Title: Fault Network Geometry Modulates Earthquake Source Spectra Across Scales
Abstract Earthquake source spectra provide unique insights into the earthquake rupture process. Motivated by previous research suggesting that complex fault geometries enhance high‐frequency seismic radiation, we study the influence of fault network geometry on earthquake source spectra using multiple independent observations. At regional scales, we examine correlations of stress drop measurements with surface fault trace misalignment in Southern California, Japan, and Central Italy. At a global scale, we examine correlations of moment‐rate function complexity of large earthquakes with focal mechanism variability, a proxy for local fault complexity. Despite significant scatter in the observations, we find overall consistent positive correlations. The concept that elastic interactions of discrete fault structures during the earthquake rupture process generates high‐frequency ground motions offers a coherent framework for interpreting our observations. These findings suggest that variations in fault complexity explain why some earthquakes produce stronger high‐frequency ground motions than others.  more » « less
Award ID(s):
2146640
PAR ID:
10693366
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
American Geophysical Union
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
52
Issue:
12
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Understanding the generation of damaging, high‐frequency ground motions during earthquakes is essential both for fundamental science and for effective hazard preparation. Various theories exist regarding the origin of high‐frequency ground motions, including the standard paradigm linked to slip heterogeneity on the rupture plane, and alternative perspectives associated with fault complexity. To assess these competing hypotheses, we measure ground motion amplitudes in different frequency bands for 3 ≤ M ≤ 5.8 earthquakes in Southern California and compare them to empirical ground motion models. We utilize a Bayesian inference technique called the Integrated Nested Laplace Approximation (INLA) to identify earthquake source regions that produce higher or lower ground motions than expected. Our analysis reveals a strong correlation between fault complexity measurements and the high‐frequency ground motion event terms identified by INLA. These findings suggest that earthquakes on complex faults (or fault networks) lead to stronger‐than‐expected ground motions at high frequencies. 
    more » « less
  2. 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
  3. 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
  4. Abstract Earthquake rupture directivity impacts ground motions and provides insights on fault zone properties and earthquake physics. However, measuring directivity of small earthquakes is challenging due to their compact rupture sizes and complex path and site effects at high frequencies. Here, we develop a new approach that deconvolves energy envelopes of the S wave trains to remove path and site effects and robustly resolve azimuthal variations in durations of apparent source‐time functions. Our method benefits from the coherence of energy envelopes for high‐frequency seismic data, which provides more stable directivity results than conventional waveform deconvolution methods. We validate our method using both synthetic tests and real observations. We apply the algorithm to determine rupture directivities of 58 magnitude 3.5–5.5 earthquakes during the 2019 Ridgecrest earthquake sequence. The rupture directivities suggest an orthogonal interlocking fault system consistent with aftershock locations. Additionally, the rupture directivity pattern appears to correlate with spatial heterogeneity in earthquake stress drops. Our energy envelope deconvolution method enables directivity measurements at smaller magnitudes than traditional approaches and has potential for constraining small earthquake rupture dynamics. 
    more » « less
  5. Abstract High-frequency energy is often challenging to model because it is influenced by numerous factors in earthquake source processes. Hence, we investigate the factors that impact the high-frequency spectra of Mw 1.5–3 earthquakes from the 2019 Ridgecrest earthquake sequence, by clustering the earthquakes and stacking their velocity spectra to average source effects. We observe an azimuthal variation in the spectral amplitudes, where stations close to fault strike record more high-frequency energy than other stations. The increase in spectral amplitudes for stations that are near fault strikes is most apparent around 15–25 Hz, and a secondary high-frequency peak can be observed in some of the clusters. Our results are consistent across clusters with average depths of 2.0–9.7 km and average magnitudes of Mw 1.6–2.7. We evaluate the factors impacting high-frequency energy by calculating the relative site effects using regional earthquakes and modeling fault zone amplification and rupture directivity using dynamic rupture simulations. The analysis shows fault zone amplification may be the primary factor causing the observed azimuthal variation of high-frequency energy, though a combination of these factors may also contribute to the variability of high-frequency velocity spectra. 
    more » « less