An official website of the United States government
Abstract There is a strong need to model potential rupture behaviors for the next Cascadia megathrust earthquake. However, there exists significant uncertainty regarding the extent of downdip rupture and rupture speed. To address this problem, we study how the transition region (i.e., the gap), which separates the locked from slow‐slip regions, influences coseismic rupture propagation using 2‐D dynamic rupture simulations governed by a slip‐weakening friction law. We show that rupture propagation through the gap is strongly controlled by the amount of accumulated tectonic initial shear stress and gap friction level. A large amplitude negative dynamic stress drop is needed to arrest downdip rupture. We also observe downdip supershear rupture when the gradient in effective normal stress from the locked to slow‐slip regions is dramatic. Our results justify kinematic rupture models that extend below the gap and suggests the possibility of high‐frequency energy radiation during the next Cascadia megathrust earthquake.
more »
« less
Working with Dynamic Earthquake Rupture Models: A Practical Guide
Abstract Dynamic rupture models are physics-based simulations that couple fracture mechanics to wave propagation and are used to explain specific earthquake observations or to generate a suite of predictions to understand the influence of frictional, geometrical, stress, and material parameters. These simulations can model single earthquakes or multiple earthquake cycles. The objective of this article is to provide a self-contained and practical guide for students starting in the field of earthquake dynamics. Senior researchers who are interested in learning the first-order constraints and general approaches to dynamic rupture problems will also benefit. We believe this guide is timely given the recent growth of computational resources and the range of sophisticated modeling software that are now available. We start with a succinct discussion of the essential physics of earthquake rupture propagation and walk the reader through the main concepts in dynamic rupture model design. We briefly touch on fully dynamic earthquake cycle models but leave the details of this topic for other publications. We also highlight examples throughout that demonstrate the use of dynamic rupture models to investigate various aspects of the faulting process.
more »
« less
- Award ID(s):
- 1943742
- PAR ID:
- 10380873
- Date Published:
- Journal Name:
- Seismological Research Letters
- Volume:
- 93
- Issue:
- 4
- ISSN:
- 0895-0695
- Page Range / eLocation ID:
- 2096 to 2110
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The 2023 M7.8 Kahramanmaraş/Pazarcik earthquake was larger and more destructive than what had been expected. Here we analyzed nearfield seismic records and developed a dynamic rupture model that reconciles different currently conflicting inversion results and reveals spatially non-uniform propagation speeds in this earthquake, with predominantly supershear speeds observed along the Narli fault and at the southwest (SW) end of the East Anatolian Fault (EAF). The model highlights the critical role of geometric complexity and heterogeneous frictional conditions in facilitating continued propagation and influencing rupture speed. We also constrained the conditions that allowed for the rupture to jump from the Narli fault to EAF and to generate the delayed backpropagating rupture towards the SW. Our findings have important implications for understanding earthquake hazards and guiding future response efforts and demonstrate the value of physics based dynamic modeling fused with near-field data in enhancing our understanding of earthquake mechanisms and improving risk assessment.more » « less
-
Physics-based dynamic rupture simulations are valuable for assessing the seismic hazard in the Cascadia subduction zone (CSZ), but require assumptions about fault stress and material properties. Geodetic slip deficit models (SDMs) may provide information about the initial stresses governing megathrust earthquake dynamics. We present a unified workflow linking SDMs to 3D dynamic rupture simulations, and 22 rupture scenarios to unravel the dynamic trade-offs of assumptions for SDMs, rigidity, and pore fluid pressure. We find that margin-wide rupture, an earthquake that ruptures the entire length of the plate boundary, requires a large slip deficit in the central CSZ. Comparisons between Gaussian and smoother, shallow-coupled SDMs show significant differences in stress distributions and rupture dynamics. Variations in depth-dependent rigidity cause competing effects, particularly in the near-trench region. Higher overall rigidity can increase fault slip but also result in lower initial shear stresses, inhibiting slip. The state of pore fluid pressure is crucial in balancing SDM-informed initial shear stresses with realistic dynamic rupture processes, especially assuming small recurrence time scaling factors. This study highlights the importance of self-consistent assumptions for rigidity and initial stresses between geodetic, structural, and dynamic rupture models, providing a foundation for future simulations of ground motions and tsunami generation.more » « less
-
Abstract Dynamic rupture simulations generate synthetic waveforms that account for nonlinear source and path complexity. Here, we analyze millions of spatially dense waveforms from 3D dynamic rupture simulations in a novel way to illuminate the spectral fingerprints of earthquake physics. We define a Brune-type equivalent near-field corner frequency (fc) to analyze the spatial variability of ground-motion spectra and unravel their link to source complexity. We first investigate a simple 3D strike-slip setup, including an asperity and a barrier, and illustrate basic relations between source properties and fc variations. Next, we analyze >13,000,000 synthetic near-field strong-motion waveforms generated in three high-resolution dynamic rupture simulations of real earthquakes, the 2019 Mw 7.1 Ridgecrest mainshock, the Mw 6.4 Searles Valley foreshock, and the 1992 Mw 7.3 Landers earthquake. All scenarios consider 3D fault geometries, topography, off-fault plasticity, viscoelastic attenuation, and 3D velocity structure and resolve frequencies up to 1–2 Hz. Our analysis reveals pronounced and localized patterns of elevated fc, specifically in the vertical components. We validate such fc variability with observed near-fault spectra. Using isochrone analysis, we identify the complex dynamic mechanisms that explain rays of elevated fc and cause unexpectedly impulsive, localized, vertical ground motions. Although the high vertical frequencies are also associated with path effects, rupture directivity, and coalescence of multiple rupture fronts, we show that they are dominantly caused by rake-rotated surface-breaking rupture fronts that decelerate due to fault heterogeneities or geometric complexity. Our findings highlight the potential of spatially dense ground-motion observations to further our understanding of earthquake physics directly from near-field data. Observed near-field fc variability may inform on directivity, surface rupture, and slip segmentation. Physics-based models can identify “what to look for,” for example, in the potentially vast amount of near-field large array or distributed acoustic sensing data.more » « less
-
null (Ed.)SUMMARY Earthquake ruptures are generally considered to be cracks that propagate as fracture or frictional slip on pre-existing faults. Crack models have been used to describe the spatial distribution of fault offset and the associated static stress changes along a fault, and have implications for friction evolution and the underlying physics of rupture processes. However, field measurements that could help refine idealized crack models are rare. Here, we describe large-scale laboratory earthquake experiments, where all rupture processes were contained within a 3-m long saw-cut granite fault, and we propose an analytical crack model that fits our measurements. Similar to natural earthquakes, laboratory measurements show coseismic slip that gradually tapers near the rupture tips. Measured stress changes show roughly constant stress drop in the centre of the ruptured region, a maximum stress increase near the rupture tips and a smooth transition in between, in a region we describe as the earthquake arrest zone. The proposed model generalizes the widely used elliptical crack model by adding gradually tapered slip at the ends of the rupture. Different from the cohesive zone described by fracture mechanics, we propose that the transition in stress changes and the corresponding linear taper observed in the earthquake arrest zone are the result of rupture termination conditions primarily controlled by the initial stress distribution. It is the heterogeneous initial stress distribution that controls the arrest of laboratory earthquakes, and the features of static stress changes. We also performed dynamic rupture simulations that confirm how arrest conditions can affect slip taper and static stress changes. If applicable to larger natural earthquakes, this distinction between an earthquake arrest zone (that depends on stress conditions) and a cohesive zone (that depends primarily on strength evolution) has important implications for how seismic observations of earthquake fracture energy should be interpreted.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1785/0220220022
