An official website of the United States government
A precise understanding of earthquake magnitudes is vital for accurate calculations of magnitude exceedance probabilities and seismic hazard assessment. However, characterization of earthquake magnitude, particularly for small events, is complicated by differences in network capabilities and procedures. Furthermore, the use of differing magnitude scales for events of various sizes introduces additional challenges and produces disparate magnitude estimates for the same events. To address the need for a consistent magnitude estimation procedure that can accurately estimate magnitude across a wide magnitude range and in diverse tectonic environments, we investigate the use of the relative magnitude method. This approach utilizes amplitude ratios of highly correlated waveforms among numerous interlinked event pairs to compute magnitude for a group of events. While the relative magnitude method is advantageous because it can be applied uniformly in various regions and does not require distance or attenuation corrections, there are several parameters that currently require human decision which may introduce bias. These include acceptable thresholds for signal-to-noise ratios and cross-correlation, filtering procedures, sampling windows, and station selection. Our research focuses on computing new relative magnitudes for events in southern California, including the 2019 Ridgecrest sequence. We investigate the uncertainty that human decision may impose on the resulting magnitudes and compare our results to other magnitude estimation methods. Finally, we present our recommendations for routine procedures that minimize uncertainty and variability in the relative magnitude method, aiming to enhance the utility of this method for future users.
more »
« less
This content will become publicly available on December 12, 2025
Quantifying Magnitude Uncertainty of the 2019 Ridgecrest Earthquake Sequence Through a Sensitivity Study of the Relative Magnitude Method
ABSTRACT Precise knowledge of earthquake magnitudes is vital for accurate characterization of seismic hazards. However, the estimation of earthquake magnitude, particularly for small events, is complicated by differences in network procedures and completeness. This produces disparate magnitude estimates for the same event and emphasizes the need for a consistent and transportable magnitude estimation procedure. Here, we investigate the use of the relative magnitude method, which measures earthquake magnitude from a least-squares inversion of interlinked waveform amplitude ratios. Our results show that that the relative magnitude method can establish both local and moment magnitudes for many events in the 2019 Ridgecrest sequence. The method also provides constraints on moment magnitude estimates for M <3 events, which are not routinely available using current methods. Although the relative magnitude method is advantageous because it can be applied uniformly in various regions and does not require empirical distance or attenuation corrections, there are several parameters that require subjective decision making and may introduce bias in the resulting magnitude estimates. These include acceptable thresholds for signal-to-noise ratios and cross correlation, filtering procedures, sampling windows, and station selection. Here, we not only calculate magnitude but also investigate how the subjective decision making affects the resulting magnitudes. Based on our analysis, we present recommendations to enhance the utility of this method for future users.
more »
« less
- Award ID(s):
- 2315814
- PAR ID:
- 10631906
- Publisher / Repository:
- GeoScienceWorld
- Date Published:
- Journal Name:
- Bulletin of the Seismological Society of America
- Volume:
- 115
- Issue:
- 3
- ISSN:
- 0037-1106
- Page Range / eLocation ID:
- 1294 to 1307
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Accurate estimates of earthquake magnitude are necessary to improve our understanding of seismic hazard. Unbiased magnitudes for small earthquakes are especially important because magnitude exceedance probabilities for large earthquakes are derived from the behavior of small earthquakes. Also, accurate characterization of small events is becoming increasingly important for ground motion models. However, catalog magnitudes may vary for the same event depending on network procedures and capabilities. In addition, different magnitude scales are often used for events of varying sizes. For example, moment magnitude (Mw) is the widely preferred estimate for earthquake size but it is often not available for small earthquakes (M < 3.5). As a result, statistical measures such as magnitude frequency distribution (MFD) and b-value can be biased depending on magnitude type and uncertainties that arise during the measurement process. In this research we demonstrate the capability of the relative magnitude method to provide a uniform and accurate estimate of earthquake magnitude in a variety of regions, while only requiring the use of waveform data. The study regions include the Permian Basin in Texas, central Oklahoma, and southern California. We present results in which only relative magnitudes are used to estimate MFD and b-value as well as relative magnitudes that are benchmarked to an absolute scale using a coda-envelope derived Mw calibration for small events. We also discuss potential sources of uncertainty in the relative magnitude method such as acceptable signal-to-noise ratios, cross-correlation thresholds, and choice of scaling constant, as well as our attempts to mitigate those uncertainties.more » « less
-
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
-
Abstract The magnitude‐frequency distribution (MFD) describes the relative proportion of earthquake magnitudes and provides vital information for seismic hazard assessment. Theb‐value, derived from the MFD, is commonly used to estimate the probability that a future earthquake will exceed a specified magnitude threshold. Improved MFD andb‐value estimates are of great importance in the central and eastern United States where high volumes of fluid injection have contributed to a significant rise in seismicity over the last decade. In this study, we recalculate the magnitudes of 8,775 events for the 2011 Prague, Oklahoma sequence using a relative magnitude approach that depends only on waveform data to calculate magnitudes. We also compare the distribution of successive magnitude differences to the MFD and show that a combination of the magnitude difference distribution (MDFD) and relative magnitudes yields a reliable estimate ofb‐value. Using the MDFD and relative magnitudes, we examine the temporal and spatial variations in theb‐value and show thatb‐value ranges between ∼0.6 and 0.85 during the aftershock sequence for at least 5 months after theM5.7 mainshock, though areas surrounding the northeast part of the sequence experience higherb‐values (0.7–0.85) than the southwestern part of the Meeker‐Prague fault whereb‐value is the lowest (0.6–0.7). We also identify a cluster of off‐fault events with the highestb‐values in the catalog (0.85). These new estimates of MFD andb‐value will contribute to understanding of the relations between induced and tectonic earthquake sequences and promote discussion regarding the use ofb‐value in induced seismic hazard estimation.more » « less
-
The clustering of earthquake magnitudes is poorly understood compared to spatial and temporal clustering. Better understanding of correlations between earthquake magnitudes could provide insight into the mechanisms of earthquake rupture and fault interactions, and improve earthquake forecasting models. In this study we present a novel method of examining how seismic magnitude clustering occurs beyond the next event in the catalog and evolves with time and space between earthquake events. We first evaluate the clustering signature over time and space using double-difference located catalogs from Southern and Northern California. The strength of magnitude clustering appears to decay linearly with distance between events and logarithmically with time. The signature persists for longer distances (more than 50km) and times (several days) than previously thought, indicating that magnitude clustering is not driven solely by repeated rupture of an identical fault patch or Omori aftershock processes. The decay patterns occur in all magnitude ranges of the catalog and are demonstrated across multiple methodologies of study. These patterns are also shown to be present in laboratory rock fracture catalogs but absent in ETAS synthetic catalogs. Incorporating magnitude clustering decay patterns into earthquake forecasting models such as ETAS could improve their accuracy.more » « less
