Off-fault plasticity in three-dimensional dynamic rupture simulations using a modal Discontinuous Galerkin method on unstructured meshes: Implementation, verification and application

Stephanie Wollherr, Alice Agnes Gabriel, Carsten Uphoff

Research output: Contribution to journalArticlepeer-review

46 Scopus citations


The dynamics and potential size of earthquakes depend crucially on rupture transfers between adjacent fault segments. To accurately describe earthquake source dynamics, numerical models can account for realistic fault geometries and rheologies such as nonlinear inelastic processes off the slip interface. We present implementation, verification and application of off-fault Drucker-Prager plasticity in the open source software SeisSol ( SeisSol is based on an arbitrary high-order derivative modal Discontinuous Galerkin method using unstructured, tetrahedral meshes specifically suited for complex geometries. Two implementation approaches are detailed, modelling plastic failure either employing subelemental quadrature points or switching to nodal basis coefficients. At fine fault discretizations, the nodal basis approach is up to six times more efficient in terms of computational costs while yielding comparable accuracy. Both methods are verified in community benchmark problems and by 3-D numerical h- and p-refinement studies with heterogeneous initial stresses. We observe no spectral convergence for on-fault quantities with respect to a given reference solution, but rather discuss a limitation to low-order convergence for heterogeneous 3-D dynamic rupture problems. For simulations including plasticity, a high fault resolution may be less crucial than commonly assumed, due to the regularization of peak slip rate and an increase of the minimum cohesive zone width. In large-scale dynamic rupture simulations based on the 1992 Landers earthquake, we observe high rupture complexity including reverse slip, direct branching and dynamic triggering. The spatiotemporal distribution of rupture transfers are altered distinctively by plastic energy absorption, correlated with locations of geometrical fault complexity. Computational cost increases by 7 per cent when accounting for off-fault plasticity in the demonstrating application. Our results imply that the combination of fully 3-D dynamic modelling, complex fault geometries and off-fault plastic yielding is important to realistically capture dynamic rupture transfers in natural fault systems.
Original languageEnglish (US)
Pages (from-to)1556-1584
Number of pages29
JournalGeophysical Journal International
Issue number3
StatePublished - May 28 2018
Externally publishedYes

Bibliographical note

KAUST Repository Item: Exported on 2022-06-09
Acknowledged KAUST grant number(s): ORS-2016-CRG5-3027-04, OSR-CRG2017-3389
Acknowledgements: The work presented in this paper was supported by the German Research Foundation (DFG; project no. KA 2281/4-1, AOBJ 584936/TG-92); by the Bavarian Competence Network for Technical and Scientific High Performance Computing (KONWIHR), project GeoPF (Geophysics for PetaFlop Computing); by the Volkswagen Foundation (project ASCETE - Advanced Simulation of Coupled Earthquake-Tsunami Events, grant no. 88479); by the European Union's Horizon 2020 research and innovation program under grant agreement no. 671698 and by King Abdullah University of Sciene and Technology (KAUST) in Thuwal, Saudi Arabia, under grant ORS-2016-CRG5-3027-04 and OSR-CRG2017-3389. Computing resourceswere provided by the Leibniz Supercomputing Centre (LRZ, projects no. h019z, pr63qo and pr45fi on SuperMUC). The authors thank Dave A. May and Kenneth C. Duru for fruitful discussions. Furthermore, we thank LMU Geophysics' system administrator team for their support (Oeser, 2009). The manuscript benefited from valuable reviews and comments by R. Ando and an anonymous reviewer.
This publication acknowledges KAUST support, but has no KAUST affiliated authors.


Dive into the research topics of 'Off-fault plasticity in three-dimensional dynamic rupture simulations using a modal Discontinuous Galerkin method on unstructured meshes: Implementation, verification and application'. Together they form a unique fingerprint.

Cite this