An official website of the United States government
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
Rupture Dynamics of Cascading Earthquakes in a Multiscale Fracture Network
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
- PAR ID:
- 10496074
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 129
- Issue:
- 3
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Earthquakes vary in size over many orders of magnitude, often rupturing in complex multifault and multievent sequences. Despite the large number of observed earthquakes, the scaling of the earthquake energy budget remains enigmatic. We propose that fundamentally different fracture processes govern small and large earthquakes. We combined seismological observations with physics-based earthquake models, finding that both dynamic weakening and restrengthening effects are non-negligible in the energy budget of small earthquakes. We established a linear scaling relationship between fracture energy and fault size and a break in scaling with slip. We applied this scaling using supercomputing and unveiled large dynamic rupture earthquake cascades involving >700 multiscale fractures within a fault damage zone. We provide a simple explanation for seismicity across all scales with implications for comprehending earthquake genesis and multifault rupture cascades.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
-
Abstract The Mediterranean Hellenic Arc subduction zone (HASZ) has generated several 8 earthquakes and tsunamis. Seismic‐probabilistic tsunami hazard assessment typically utilizes uniform or stochastic earthquake models, which may not represent dynamic rupture and tsunami generation complexity. We present an ensemble of ten 3D dynamic rupture earthquake scenarios for the HASZ, utilizing a realistic slab geometry. Our simplest models use uniform along‐arc pre‐stresses or a single circular initial stress asperity. We then introduce progressively more complex models varying initial shear stress along‐arc, multiple asperities based on scale‐dependent critical slip weakening distance, and a most complex model blending all aforementioned heterogeneities. Thereby, regional initial conditions are constrained without relying on detailed geodetic locking models. Varying epicentral locations in the simplest, homogeneous model leads to different rupture speeds and moment magnitudes. We observe dynamic fault slip penetrating the shallow slip‐strengthening region and affecting seafloor uplift. Off‐fault plastic deformation can double vertical seafloor uplift. A single‐asperity model generates a 8 scenario resembling the 1303 Crete earthquake. Using along‐strike varying initial stresses results in 8.0–8.5 dynamic rupture scenarios with diverse slip rates and uplift patterns. The model with the most heterogeneous initial conditions yields a 7.5 scenario. Dynamic rupture complexity in prestress and fracture energy tends to lower earthquake magnitude but enhances tsunamigenic displacements. Our results offer insights into the dynamics of potential large Hellenic Arc megathrust earthquakes and may inform future physics‐based joint seismic and tsunami hazard assessments.more » « less
-
Previous geodetic and teleseismic observations of the 2021 Mw 7.4 Maduo earthquake imply surprising but difficult-to-constrain complexity, including rupture across multiple fault segments and supershear rupture. Here, we present an integrated analysis of multi-fault 3D dynamic rupture models, high-resolution optical correlation analysis, and joint optical-InSAR slip inversion. Our preferred model, validated by the teleseismic multi-peak moment rate release, includes unilateral eastward double-onset supershear speeds and cascading rupture dynamically triggering two adjacent fault branches. We propose that pronounced along-strike variation in fracture energy, complex fault geometries, and multi-scale variable prestress drives this event's complex rupture dynamics. We illustrate how supershear transition has signatures in modeled and observed off-fault deformation. Our study opens new avenues to combine observations and models to better understand complex earthquake dynamics, including local and potentially repeating supershear episodes across immature faults or under heterogeneous stress and strength conditions, which are potentially not unusual.more » « less
