Skip to main content

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

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

To refer to this page use:
Abstract: Accurate 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.
Publication Date: 22-Nov-2018
Citation: Gharti, 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.
DOI: doi:10.1093/gji/ggy495
ISSN: 0956-540X
EISSN: 1365-246X
Pages: 1364 - 1393
Type of Material: Journal Article
Journal/Proceeding Title: Geophysical Journal International
Version: Final published version. Article is made available in OAR by the publisher's permission or policy.

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