An official website of the United States government
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
Seismic magnitude clustering is prevalent in field and laboratory catalogs
Abstract Clustering of earthquake magnitudes is still actively debated, compared to well-established spatial and temporal clustering. Magnitude clustering is not currently implemented in earthquake forecasting but would be important if larger magnitude events are more likely to be followed by similar sized events. Here we show statistically significant magnitude clustering present in many different field and laboratory catalogs at a wide range of spatial scales (mm to 1000 km). It is universal in field catalogs across fault types and tectonic/induced settings, while laboratory results are unaffected by loading protocol or rock types and show temporal stability. The absence of clustering can be imposed by a global tensile stress, although clustering still occurs when isolating to triggered event pairs or spatial patches where shear stress dominates. Magnitude clustering is most prominent at short time and distance scales and modeling indicates >20% repeating magnitudes in some cases, implying it can help to narrow physical mechanisms for seismogenesis.
more »
« less
- Award ID(s):
- 2025073
- PAR ID:
- 10496121
- Publisher / Repository:
- Nature
- Date Published:
- Journal Name:
- Nature Communications
- Volume:
- 14
- Issue:
- 1
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
An accurate assessment of seismic hazard requires a combination of earthquake physics and statistical analysis. Because of the limitations in the investigation of the seismogenic sources and of the short temporal intervals covered by earthquake catalogs, laboratory experiments have played a crucial role in improving our understanding of earthquake phenomena. However, differences exist between acoustic emissions in the lab, events in small, regulated systems (e.g., mines) and natural seismicity. One of the most pressing issues concerns the role of mechanical parameters and how they affect seismic activity across boundary conditions and spatial-temporal scales. Here, we focus on fault friction. There is evidence inferred from geodesy, computational simulations and seismological investigations that most large faults are weak and characterized by very low static friction coefficients which are inconsistent with those of smaller faults and laboratory experiments. We support the hypothesis that static friction decreases with fault size due to the presence of fabrics, roughness, structural asperities and network geometry. We also model its scaling behavior as dependent on a few physical properties (e.g., fault fractal dimension). Conversely, dynamic coefficients are not affected by the spatial scale. Mathematical derivations are based on the hypothesis that earthquake onset results from fracture instability controlled by the extremes of fault shear strength. We validate this using a simple model for earthquake occurrence rooted in fracture mechanics, which reproduces key features of major seismicity (i.e., interevent time distribution, clustering and frequency-size relationship).more » « less
-
Abstract We introduce an algorithm for declustering earthquake catalogs based on the nearest‐neighbor analysis of seismicity. The algorithm discriminates between background and clustered events by random thinning that removes events according to a space‐varying threshold. The threshold is estimated using randomized‐reshuffled catalogs that are stationary, have independent space and time components, and preserve the space distribution of the original catalog. Analysis of catalog produced by the Epidemic Type Aftershock Sequence model demonstrates that the algorithm correctly classifies over 80% of background and clustered events, correctly reconstructs the stationary and space‐dependent background intensity, and shows high stability with respect to random realizations (over 75% of events have the same estimated type in over 90% of random realizations). The declustering algorithm is applied to the global Northern California Earthquake Data Center catalog with magnitudesm≥ 4 during 2000–2015; a Southern California catalog withm≥ 2.5, 3.5 during 1981–2017; an area around the 1992 Landers rupture zone withm≥ 0.0 during 1981–2015; and the Parkfield segment of San Andreas fault withm≥ 1.0 during 1984–2014. The null hypotheses of stationarity and space‐time independence are not rejected by several tests applied to the estimated background events of the global and Southern California catalogs with magnitude ranges Δm< 4. However, both hypotheses are rejected for catalogs with larger range of magnitudes Δm> 4. The deviations from the nulls are mainly due to local temporal fluctuations of seismicity and activity switching among subregions; they can be traced back to the original catalogs and represent genuine features of background seismicity.more » « less
-
Abstract We discuss data of three laboratory stick‐slip experiments on Westerly Granite samples performed at elevated confining pressure and constant displacement rate on rough fracture surfaces. The experiments produced complex slip patterns including fast and slow ruptures with large and small fault slips, as well as failure events on the fault surface producing acoustic emission bursts without externally‐detectable stress drop. Preparatory processes leading to large slips were tracked with an ensemble of ten seismo‐mechanical and statistical parameters characterizing local and global damage and stress evolution, localization and clustering processes, as well as event interactions. We decompose complex spatio‐temporal trends in the lab‐quake characteristics and identify persistent effects of evolving fault roughness and damage at different length scales, and local stress evolution approaching large events. The observed trends highlight labquake localization processes on different spatial and temporal scales. The preparatory process of large slip events includes smaller events marked by confined bursts of acoustic emission activity that collectively prepare the fault surface for a system‐wide failure by conditioning the large‐scale stress field. Our results are consistent overall with an evolving process of intermittent criticality leading to large failure events, and may contribute to improved forecasting of large natural earthquakes.more » « less
-
Abstract We present an analysis of magnitude clustering of microfractures inferred from acoustic emissions (AEs) during stick‐slip (SS) dynamics of faulted Westerly granite samples in frictional sliding experiments, with and without fluids, under triaxial loading with constant displacement rate. We investigate magnitude clustering in time across periods during, preceding and after macroscopic slip events on laboratory faults. Our findings reveal that magnitude clustering exists such that subsequent AEs tend to have more similar magnitudes than expected. Yet, this clustering only exists during macroscopic slip events and is strongest during major slip events in fluid‐saturated and dry samples. We demonstrate that robust magnitude clustering arises from variations in frequency‐magnitude distributions of AE events during macroscopic slip events. These temporal variations indicate a prevalence of larger AE events right after (0.3–3 s) the SS onset. Hence, magnitude clustering is a consequence of non‐stationarities.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1038/s41467-023-37782-5
