ABSTRACT We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress-drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude-spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress-drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two-week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. In this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite-fault modeling, ground-motion-based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between-method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress-drop values.
more »
« less
Improving Stress Drop Estimation through Point-Wise Spectral Ratio Stacking
Stress drop, a crucial source parameter in earthquake studies, significantly influences ground motion prediction and seismic hazard assessment. Despite several existing methods to estimate stress drops, the resulting stress drop estimates often exhibit a wide variation of up to 3-4 orders of magnitude. In this study, we address the robustness of stress drop estimation by introducing a point-wise spectral ratio stacking approach based on empirical Green’s functions (eGfs). Conventional trace-wise stacking can lead to data exclusion due to high signal-to-noise ratio requirements across a wide range of frequency. By adopting point-wise stacking, we maximize the utilization of useful recording information, leading to more accurate stress drop estimates. We applied the point-wise spectral ratio stacking method to a comprehensive dataset comprising global earthquakes from 1990 to 2020 with magnitude larger than Mw5.5 and depth shallower than 50 km. We first verified the moment magnitudes of earthquakes estimated from the resulting seismic moment ratios. We found that the moment magnitude of master events best consistent with catalog magnitudes when the magnitude difference between master and their eGfs differs by about 0.5. Our analysis indicates that stress drop of shallow earthquakes exhibits no depth dependence, while showing a slight increase with magnitude. The results obtained through our optimized stacking process shed new light on stress drop estimate of shallow earthquakes and have the potential to enhance the understanding of earthquake mechanics.
more »
« less
- Award ID(s):
- 2019379
- PAR ID:
- 10535372
- Publisher / Repository:
- AGU
- Date Published:
- Format(s):
- Medium: X
- Location:
- San Francisco, CA
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT Earthquake stress drop—a key parameter for describing the energetics of earthquake rupture—can be estimated in several different, but theoretically equivalent, ways. However, independent estimates for the same earthquakes sometimes differ significantly. We find that earthquake source complexity plays a significant role in why theoretically (for simple rupture models) equivalent methods produce different estimates. We apply time- and frequency-domain methods to estimate stress drops for real earthquakes in the SCARDEC (Seismic source ChAracteristics Retrieved from DEConvolving teleseismic body waves, Vallée and Douet, 2016) source time function (STF) database and analyze how rupture complexity drives stress-drop estimate discrepancies. Specifically, we identify two complexity metrics—Brune relative energy (BRE) and spectral decay—that parameterize an earthquake’s complexity relative to the standard Brune model and strongly correlate with the estimate discrepancies. We find that the observed systematic magnitude–stress-drop trends may reflect underlying changes in STF complexity, not necessarily trends in actual stress drop. Both the decay and BRE parameters vary systematically with magnitude, but whether this magnitude–complexity relationship is real remains unresolved.more » « less
-
ABSTRACT Although the Brune source model describes earthquake moment release as a single pulse, it is widely used in studies of complex earthquakes with multiple episodes of high moment release (i.e., multiple subevents). In this study, we investigate how corner frequency estimates of earthquakes with multiple subevents are biased if they are based on the Brune source model. By assuming complex sources as a sum of multiple Brune sources, we analyze 1640 source time functions of Mw 5.5–8.0 earthquakes in the seismic source characteristic retrieved from deconvolving teleseismic body waves catalog to estimate the corner frequencies, onset times, and seismic moments of subevents. We identify more subevents for strike-slip earthquakes than dip-slip earthquakes, and the number of resolvable subevents increases with magnitude. We find that earthquake corner frequency correlates best with the corner frequency of the subevent with the highest moment release (i.e., the largest subsevent). This suggests that, when the Brune model is used, the estimated corner frequency and, therefore, the stress drop of a complex earthquake is determined primarily by the largest subevent rather than the total rupture area. Our results imply that, in addition to the simplified assumption of a radial rupture area with a constant rupture velocity, the stress variation of asperities, rather than the average stress change of the whole fault, contributes to the large variance of stress-drop estimates.more » « less
-
Abstract Elevated seismic noise for moderate‐size earthquakes recorded at teleseismic distances has limited our ability to see their complexity. We develop a machine‐learning‐based algorithm to separate noise and earthquake signals that overlap in frequency. The multi‐task encoder‐decoder model is built around a kernel pre‐trained on local (e.g., short distances) earthquake data (Yin et al., 2022,https://doi.org/10.1093/gji/ggac290) and is modified by continued learning with high‐quality teleseismic data. We denoise teleseismic P waves of deep Mw5.0+ earthquakes and use the clean P waves to estimate source characteristics with reduced uncertainties of these understudied earthquakes. We find a scaling of moment and duration to beM0 ≃ τ4, and a resulting strong scaling of stress drop and radiated energy with magnitude ( and ). The median radiation efficiency is 5%, a low value compared to crustal earthquakes. Overall, we show that deep earthquakes have weak rupture directivity and few subevents, suggesting a simple model of a circular crack with radial rupture propagation is appropriate. When accounting for their respective scaling with earthquake size, we find no systematic depth variations of duration, stress drop, or radiated energy within the 100–700 km depth range. Our study supports the findings of Poli and Prieto (2016,https://doi.org/10.1002/2016jb013521) with a doubled amount of earthquakes investigated and with earthquakes of lower magnitudes.more » « less
-
The 2019 Ridgecrest, CA earthquake sequence has provided a unique opportunity and a rich dataset to understand earthquake source properties and near-fault structure. Using the high-quality seismic data provided by the SCEC Stress Drop Validation group, we first estimate the corner frequency of M2.0-4.5 earthquakes by applying the spectral ratio method based on empirical Green’s function (Liu et al., 2020). We relate corner frequency estimates to stress drops assuming the Brune source model and circular cracks. Our preliminary results show increasing median stress drops with magnitude for both P and S waves, from 1 MPa for M2.0 events to 10 MPa for M4.0 events, though the limited frequency bandwidth may cause underestimation for small events. The estimated moment magnitude is proportional to the catalog magnitude by a factor of 0.72, which is close to 0.74 estimated by Trugman (2020) for the Ridgecrest earthquake sequence. In the second part of the study, we examine the impact of fault zone structure on the azimuthal variation of the source spectra. Using kinematic simulations and observations of the 2003 Big Bear earthquake sequence, Huang et al. (2016) showed that fault damage zones can act as an effective wave guide and cause high-frequency wave amplification along directions close to fault strike. We use clusters of M1.5-3 earthquakes in the Ridgecrest region to further examine the azimuthal variation of the stacked source spectra and investigate if the near-source structure can affect our corner frequency estimates. We aim to develop robust methods that utilize high-quality seismic data to illuminate earthquake source processes and fault zone properties.more » « less
An official website of the United States government

