A Discontinuous Galerkin Method for Simulating 3D Seismic Wave Propagation in Nonlinear Rock Models: Verification and Application to the 2015 Mw 7.8 Gorkha Earthquake

TL;DR

This paper introduces a 3D discontinuous Galerkin method incorporating nonlinear damage models to simulate seismic wave propagation, capturing co-seismic wave speed reductions and ground motion amplification during the 2015 Gorkha earthquake.

Contribution

The study develops an open-source, scalable 3D discontinuous Galerkin simulation framework with nonlinear damage rheology, advancing modeling of regional seismic wave speed changes and ground motion effects.

Findings

Wave speed reductions ranged from <0.01% to >50% near the source.

Nonlinear damage models better capture low-frequency ground motion amplification.

The method scales efficiently on high-performance computing systems.

Abstract

The nonlinear mechanical responses of rocks and soils to seismic waves play an important role in earthquake physics, influencing ground motion from source to site. Continuous geophysical monitoring, such as ambient noise interferometry, has revealed co-seismic wave speed reductions extending tens of kilometers from earthquake sources. However, the mechanisms governing these changes remain challenging to model, especially at regional scales. Using a nonlinear damage model constrained by laboratory experiments, we develop and apply an open-source 3D discontinuous Galerkin method to simulate regional co-seismic wave speed changes during the 2015 Mw7.8 Gorkha earthquake. We find pronounced spatial variations of co-seismic wave speed reduction, ranging from <0.01% to >50%, particularly close to the source and within the Kathmandu Basin. The most significant reduction occurs within the sedimentary basin and varies with basin depths, while wave speed reductions correlate with the fault slip distribution near the source. By comparing ground motions from simulations with elastic, viscoelastic, elastoplastic, and nonlinear damage rheologies, we demonstrate that the nonlinear damage model effectively captures low-frequency ground motion amplification due to strain-dependent wave speed reductions in soft sediments. We verify the accuracy of our approach through comparisons with analytical solutions and assess its scalability on high-performance computing systems. The model shows near-linear strong and weak scaling up to 2048 nodes, enabling efficient large-scale simulations. Our findings provide a physics-based framework to quantify nonlinear earthquake effects and emphasize the importance of damage-induced wave speed variations for seismic hazard assessment and ground motion predictions.