An official website of the United States government
Abstract Using a novel high‐performance computing implementation of a nonlinear continuum damage‐breakage model, we explore interactions between 3D co‐seismic off‐fault damage, seismic radiation, and rupture dynamics. Our simulations demonstrate that off‐fault damage enhances high‐frequency wave radiation above 1 Hz, reduces rupture speed and alters the total kinetic energy. We identify distinct damage regimes separated by solid‐granular transition, with smooth distributions under low damage conditions transitioning to localized, mesh‐independent shear bands upon reaching brittle failure. The shear band orientations depend systematically on the background stress and agree with analytical predictions. The brittle damage inhibits transitions to supershear rupture propagation and the rupture front strain field results in locally reduced damage accumulation during supershear transition. The dynamically generated damage yields uniform and isotropic ratios of fault‐normal to fault‐parallel high‐frequency ground motions. Co‐seismic damage zones exhibit depth‐dependent width variations, becoming broader near the Earth's surface consistent with field observations, even under uniform stress conditions. We discover a new delayed dynamic triggering mechanism in multi‐fault systems, driven by reductions in elastic moduli and the ensuing stress heterogeneities in 3D tensile fault step‐overs. This mechanism affects the static and dynamic stress fields and includes the formation of high shear‐traction fronts around localized damage zones. The brittle damage facilitates rupture cascading across faults, linking delay times directly to damage rheology and fault zone evolution. Our results help explain near‐fault high‐frequency isotropic radiation and delayed rupture triggering, improving our understanding of earthquake processes, seismic wavefields and fault system interactions.
more »
« less
A Discontinuous Galerkin Method for Simulating 3D Seismic Wave Propagation in Nonlinear Rock Models: Verification and Application to the M w 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 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, 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.
more »
« less
- PAR ID:
- 10633045
- Publisher / Repository:
- AGU
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 130
- Issue:
- 7
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Abstract We measure pseudospectral and peak ground motions from 44 intermediate‐depth Mw≥4.9 earthquakes in the Cook Inlet region of southern Alaska, including those from the 2018 Mw 7.1 earthquake near Anchorage, to identify regional amplification features (0.1–5 s period). Ground‐motion residuals are computed with respect to an empirical ground‐motion model for intraslab subduction earthquakes, and we compute bias, between‐, and within‐event terms through a linear mixed‐effects regression. Between‐event residuals are analyzed to assess the relative source characteristics of the Cook Inlet earthquakes and suggest a difference in the scaling of the source with depth, relative to global observations. The within‐event residuals are analyzed to investigate regional amplification, and various spatial patterns manifest, including correlations of amplification with depth of the Cook Inlet basin and varying amplifications east and west of the center of the basin. Three earthquake clusters are analyzed separately and indicate spatial amplification patterns that depend on source location and exhibit variations in the depth scaling of long‐period basin amplification. The observations inform future seismic hazard modeling efforts in the Cook Inlet region. More broadly, they suggest a greater complexity of basin and regional amplification than is currently used in seismic hazard analyses.more » « less
-
ABSTRACT We simulate shaking in Tacoma, Washington, and surrounding areas from Mw 6.5 and 7.0 earthquakes on the Tacoma fault. Ground motions are directly modeled up to 2.5 Hz using kinematic, finite-fault sources; a 3D seismic velocity model considering regional geology; and a model mesh with 30 m sampling at the ground surface. In addition, we explore how adjustments to the seismic velocity model affect predicted shaking over a range of periods. These adjustments include the addition of a region-specific geotechnical gradient, surface topography, and a fault damage zone. We find that the simulated shaking tends to be near estimates from empirical ground-motion models (GMMs). However, long-period (T = 5.0 s) shaking within the Tacoma basin is typically underpredicted by the GMMs. The fit between simulated and GMM-derived short-period (T = 0.5 s) shaking is significantly improved with the addition of the geotechnical gradient. From comparing different Mw 6.5 earthquake scenarios, we also find that the response of the Tacoma basin is sensitive to the azimuth of incoming seismic waves. In adding surface topography to the simulation, we find that average ground motion is similar to that produced from the nontopography model. However, shaking is often amplified at topographic highs and deamplified at topographic lows, and the wavefield undergoes extensive scattering. Adding a fault damage zone has the effect of amplifying short-period shaking adjacent to the fault, while reducing far-field shaking. Intermediate-period shaking is amplified within the Tacoma basin, likely due to enhanced surface-wave generation attributable to the fault damage zone waveguide. When applied in the same model, the topography and fault damage zone adjustments often enhance or reduce the effects of one another, adding further complexity to the wavefield. These results emphasize the importance of improving near-surface velocity model resolution as waveform simulations progress toward higher frequencies.more » « less
-
ABSTRACT Fault zones exhibit geometrical complexity and are often surrounded by multiscale fracture networks within their damage zones, potentially influencing rupture dynamics and near-field ground motions. In this study, we investigate the ground-motion characteristics of cascading ruptures across damage zone fracture networks of moderate-size earthquakes (Mw 5.5–6.0) using high-resolution 3D dynamic rupture simulations. Our models feature a listric normal fault surrounded by more than 800 fractures, emulating a major fault and its associated damage zone. We analyze three cases: a cascading rupture propagating within the fracture network (Mw 5.5), a non-cascading main-fault rupture with off-fault fracture slip (Mw 6.0), and a main-fault rupture without a fracture network (Mw 6.0). Cascading ruptures within the fracture network produce distinct ground-motion signatures with enriched high-frequency content, arising from simultaneous slip of multiple fractures and parts of the main fault, resembling source coda-wave-like signatures. This case shows elevated near-field characteristic frequency (fc) and stress drop, approximately an order of magnitude higher than the estimation directly on the fault of the dynamic rupture simulation. The inferred fc of the modeled vertical ground-motion components reflects the complexity of the radiation pattern and rupture directivity of fracture-network cascading earthquakes. We show that this is consistent with observations of strong azimuthal dependence of corner frequency in the 2009–2016 central Apennines, Italy, earthquake, sequence. Simulated ground motions from fracture-network cascading ruptures also show pronounced azimuthal variations in peak ground acceleration (PGA), peak ground velocity, and pseudospectral acceleration, with average PGA nearly double that of the non-cascading cases. Cascading ruptures radiate high-frequency seismic energy, yield nontypical ground-motion characteristics including coda-wave-like signatures, and may result in a significantly higher seismologically inferred stress drop and PGA. Such outcomes emphasize the critical role of fault-zone complexity in affecting rupture dynamics and seismic radiation and have important implications for physics-based seismic hazard assessment.more » « less
-
Abstract The San Fernando Valley (SFV), a densely populated region in Southern California, has high earthquake hazard due to a complex network of active faults and the amplifying effects of the sedimentary basin. Since the devastating 1994 Mw 6.7 Northridge earthquake, numerous studies have examined its structure using various geological and geophysical datasets. However, current seismic velocity models still lack the resolution to accurately image the near-surface velocity structure and concealed or blind faults, which are critical for high-frequency wavefield simulations and earthquake hazard modeling. To address these challenges, we develop a 3D high-resolution shear-wave velocity model for the SFV using ambient noise data from a dense array of 140 seismic nodes and 10 Southern California Seismic Network stations. We also invert gravity data to map the basin geometry and integrate horizontal-to-vertical spectral ratios and aeromagnetic data to constrain interfaces and map major geological structures. With a lateral resolution of 250 m near the basin center, our model reveals previously unresolved geological features, including the detailed geometry of the basin and previously unmapped structure of faults at depth. The basin deepens from the Santa Monica Mountains in the south to approximately 4 km near its center and 7 km in the Sylmar sub-basin at the basin’s northern margin. Strong velocity contrasts are observed across major faults, at the basin edges, and in the basin’s upper 500 m, for which we measure velocities as low as 200 m/s. Our high-resolution model will enhance ground-motion simulations and earthquake hazard assessments for the SFV and has implications for other urban areas with high seismic risk.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1029/2025JB031378
