An official website of the United States government
Abstract Fault‐damage zones comprise multiscale fracture networks that may slip dynamically and interact with the main fault during earthquake rupture. Using 3D dynamic rupture simulations and scale‐dependent fracture energy, we examine dynamic interactions of more than 800 intersecting multiscale fractures surrounding a listric fault, emulating a major listric fault and its damage zone. We investigate 10 distinct orientations of maximum horizontal stress, probing the conditions necessary for sustained slip within the fracture network or activating the main fault. Additionally, we assess the feasibility of nucleating dynamic rupture earthquake cascades from a distant fracture and investigate the sensitivity of fracture network cascading rupture to the effective normal stress level. We model either pure cascades or main fault rupture with limited off‐fault slip. We find that cascading ruptures within the fracture network are dynamically feasible under certain conditions, including: (a) the fracture energy scales with fracture and fault size, (b) favorable relative pre‐stress of fractures within the ambient stress field, and (c) close proximity of fractures. We find that cascading rupture within the fracture network discourages rupture on the main fault. Our simulations suggest that fractures with favorable relative pre‐stress, embedded within a fault damage zone, may lead to cascading earthquake rupture that shadows main fault slip. We find that such off‐fault events may reach moment magnitudes up toMw ≈ 5.5, comparable to magnitudes that can be otherwise hosted by the main fault. Our findings offer insights into physical processes governing cascading earthquake dynamic rupture within multiscale fracture networks.
more »
« less
Ground-Motion Characteristics of Cascading Earthquakes in a Multiscale Fracture Network
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
- PAR ID:
- 10633051
- Publisher / Repository:
- AGU
- Date Published:
- Journal Name:
- Bulletin of the Seismological Society of America
- ISSN:
- 0037-1106
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The 2023 Turkey earthquake sequence involved unexpected ruptures across numerous fault segments. We present 3D dynamic rupture simulations to illuminate the complex dynamics of the earthquake doublet. Our models are constrained by observations available within days of the sequence and deliver timely, mechanically consistent explanations of the unforeseen rupture paths, diverse rupture speeds, multiple slip episodes, heterogeneous fault offsets, locally strong shaking, and fault system interactions. Our simulations link both earthquakes, matching geodetic and seismic observations and reconciling regional seismotectonics, rupture dynamics, and ground motions of a fault system represented by 10 curved dipping segments and embedded in a heterogeneous stress field. The Mw 7.8 earthquake features delayed backward branching from a steeply branching splay fault, not requiring supershear speeds. The asymmetrical dynamics of the distinct, bilateral Mw 7.7 earthquake are explained by heterogeneous fault strength, prestress orientation, fracture energy, and static stress changes from the previous earthquake. Our models explain the northward deviation of its eastern rupture and the minimal slip observed on the Sürgü fault. 3D dynamic rupture scenarios can elucidate unexpected observations shortly after major earthquakes, providing timely insights for data-driven analysis and hazard assessment toward a comprehensive, physically consistent understanding of the mechanics of multifault systems.more » « less
-
Abstract Seismic faults are surrounded by damaged rocks with reduced rigidity and enhanced attenuation. These damaged fault zone structures can amplify seismic waves and affect earthquake dynamics, yet they are typically omitted in physics‐based regional ground motion simulations. We report on the significant effects of a shallow, flower‐shaped fault zone in foreshock‐mainshock 3D dynamic rupture models of the 2019 Ridgecrest earthquake sequence. We find that the fault zone structure both amplifies and reduces ground motions not only locally but at distances exceeding 100 km. This impact on ground motions is frequency‐ and magnitude‐dependent, particularly affecting higher frequency ground motions from the foreshock because its corner frequency is closer to the fault zone's fundamental eigenfrequency. Within the fault zone, the shallow transition to a velocity‐strengthening frictional regime leads to a depth‐dependent peak slip rate increase of up to 70% and confines fault zone‐induced supershear transitions mostly to the fault zone's velocity‐weakening roots. However, the interplay of fault zone waves, free surface reflections, and rupture directivity can generate localized supershear rupture, even in narrow velocity‐strengthening regions, which are typically thought to inhibit supershear rupture. This study demonstrates that shallow fault zone structures may significantly affect intermediate‐ and far‐field ground motions and cause localized supershear rupture penetrating into velocity‐strengthening regions, with important implications for seismic hazard assessment.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
-
Fault damage zones can influence various aspects of the earthquake cycle, such as the recurrence intervals and magnitudes of large earthquakes. The properties and structure of fault damage zones are often characterized using dense arrays of seismic stations located directly above the faults. However, such arrays may not always be available. Hence, our research aims to develop a novel method to image fault damage zones using broadband stations at relatively larger distances. Previous kinematic simulations and a case study of the 2003 Big Bear earthquake sequence demonstrated that fault damage zones can act as effective waveguides, amplifying high-frequency waves along directions close to fault strike via multiple reflections within the fault damage zone. The amplified high-frequency energy can be observed by stacking P-wave spectra of earthquake clusters with highly-similar waveforms (Huang et al., 2016), and the frequency band which is amplified may be used to estimate the width and velocity contrast of the fault damage zone. We attempt to identify the high-frequency peak associated with fault zone waves in stacked spectra by conducting a large-scale study of small earthquakes (M1.5–3). We use high quality broadband data from seismic stations at hypocentral distances of 20-80 km in the 2019 Ridgecrest earthquake regions. First, we group the Ridgecrest earthquakes in clusters by their locations and their waveform similarity, and then stack their velocity spectra to average the source effects of individual earthquakes. Our results show that the stations close to the fault strike record more high-frequency energies around the characteristic frequency of fault zone reflections. We find that the increase in the amount of high-frequencies is consistent across clusters with average magnitudes ranging from 1.6-2.4, which suggests that the azimuthal variation in spectra is caused by fault zone amplification rather than rupture directivity. We will apply our method to other fault zones in California, in order to search for fault damage zone structures and estimate their material properties.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1785/0120240281
