A Discontinuous Galerkin Method for Simulating 3D Seismic Wave Propagation in Nonlinear Rock Models: Verification and Application to the M<sub><i>w</i></sub> 7.8 2015 Gorkha, Nepal Earthquake

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 M w 7.8 Gorkha earthquake. We find pronounced spatial variations of co‐seismic wave speed reduction, ranging from &lt;0.01% to &gt;50%, particularly close to the source and within the Kathmandu Basin, while disagreement with observations remains. The most significant reduction occurs within the sedimentary basin and varies with basin depths, whereas 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 2,048 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.

Authors

Institutions

Publication Details

Journal
Journal of Geophysical Research Solid Earth
Published
2025-07-01
DOI
https://doi.org/10.1029/2025jb031378
Citations
2
Primary Topic
Seismic Waves and Analysis
Type
article
Field-Weighted Citation Impact
3.73

Funders

Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

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

Alice‐Agnes Gabriel, Zihua Niu, Thomas Ulrich, Vladimir Lyakhovsky et al.
2 citations
Journal of Geophysical Research Solid Earth
Seismic Waves and Analysis
3.73
article

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

Alice‐Agnes Gabriel, Zihua Niu, Thomas Ulrich, Vladimir Lyakhovsky, Heiner Igel, Sebastian Wolf
article en
2 citations

Abstract

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 M w 7.8 Gorkha earthquake. We find pronounced spatial variations of co‐seismic wave speed reduction, ranging from &lt;0.01% to &gt;50%, particularly close to the source and within the Kathmandu Basin, while disagreement with observations remains. The most significant reduction occurs within the sedimentary basin and varies with basin depths, whereas 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 2,048 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.

Journal of Geophysical Research Solid EarthVol. 130(7)
Scripps Institution of Oceanography (US), University of California San Diego (US), Geological Survey of Israel (IL), Technical University of Munich (DE), Ludwig-Maximilians-Universität München (DE)
National Science Foundation, National Aeronautics and Space Administration, HORIZON EUROPE European Research Council
Openalex Percentile: Top 8%
Seismic Waves and Analysis
3.73
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.