Skip to main content

Spectral-infinite-element simulations of coseismic and post-earthquake deformation

Author(s): Gharti, Hom Nath; Langer, Leah; Tromp, Jeroen

Download
To refer to this page use: http://arks.princeton.edu/ark:/88435/pr1zc7rt5g
Full metadata record
DC FieldValueLanguage
dc.contributor.authorGharti, Hom Nath-
dc.contributor.authorLanger, Leah-
dc.contributor.authorTromp, Jeroen-
dc.date.accessioned2022-01-25T14:59:29Z-
dc.date.available2022-01-25T14:59:29Z-
dc.date.issued2018-11-22en_US
dc.identifier.citationGharti, Hom Nath, Leah Langer, and Jeroen Tromp. "Spectral-infinite-element simulations of coseismic and post-earthquake deformation." Geophysical Journal International 216, no. 2 (2019): 1364-1393. doi:10.1093/gji/ggy495.en_US
dc.identifier.issn0956-540X-
dc.identifier.urihttp://arks.princeton.edu/ark:/88435/pr1zc7rt5g-
dc.description.abstractAccurate and efficient simulations of coseismic and post-earthquake deformation are important for proper inferences of earthquake source parameters and subsurface structure. These simulations are often performed using a truncated half-space model with approximate boundary conditions. The use of such boundary conditions introduces inaccuracies unless a sufficiently large model is used, which greatly increases the computational cost. To solve this problem, we develop a new approach by combining the spectral-element method with the mapped infinite-element method. In this approach, we still use a truncated model domain, but add a single outer layer of infinite elements. While the spectral elements capture the domain, the infinite elements capture the far-field boundary conditions. The additional computational cost due to the extra layer of infinite elements is insignificant. Numerical integration is performed via Gauss–Legendre–Lobatto and Gauss–Radau quadratures in the spectral and infinite elements, respectively. We implement an equivalent moment-density tensor approach and a split-node approach for the earthquake source, and discuss the advantages of each method. For post-earthquake deformation, we implement a general Maxwell rheology using a second-order accurate and unconditionally stable recurrence algorithm. We benchmark our results with the Okada analytical solutions for coseismic deformation, and with the Savage & Prescott analytical solution and the PyLith finite-element code for post-earthquake deformation.en_US
dc.format.extent1364 - 1393en_US
dc.language.isoen_USen_US
dc.relation.ispartofGeophysical Journal Internationalen_US
dc.rightsFinal published version. Article is made available in OAR by the publisher's permission or policy.en_US
dc.titleSpectral-infinite-element simulations of coseismic and post-earthquake deformationen_US
dc.typeJournal Articleen_US
dc.identifier.doidoi:10.1093/gji/ggy495-
dc.identifier.eissn1365-246X-
pu.type.symplectichttp://www.symplectic.co.uk/publications/atom-terms/1.0/journal-articleen_US

Files in This Item:
File Description SizeFormat 
Spectral-infinite-element_simulations_coseismic_post-earthquake_deformation.pdf9.57 MBAdobe PDFView/Download


Items in OAR@Princeton are protected by copyright, with all rights reserved, unless otherwise indicated.