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 magnitudes
This content will become publicly available on December 1, 2024
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
- NSF-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
-
more » « less
Abstract m ≥ 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 Megathrust earthquakes and their associated tsunamis cause some of the worst natural disasters. In addition to earthquakes, a wide range of slip behaviors are present at subduction zones, including slow earthquakes that span multiple orders of spatial and temporal scales. Understanding these events may shed light on the stress or strength conditions of the megathrust fault. Out of all types of slow earthquakes, very low frequency earthquakes (VLFEs) are most enigmatic because they are difficult to detect reliably, and the physical nature of VLFEs are poorly understood. Here we show three VLFEs in Cascadia that were dynamically triggered by a 2009 Mw 6.9 Canal de Ballenas earthquake in the Gulf of California. The VLFEs likely locate in between the seismogenic zone and the Cascadia episodic tremor and slip (ETS) zone, including one event with a moment magnitude of 5.7. This is the largest VLFE reported to date, causing clear geodetic signals. Our results show that the Cascadia megathrust fault might slip rapidly at some spots in this gap zone, and such a permissible slip behavior has direct seismic hazard implications for coastal communities and perhaps further inland. Further, the observed seismic sources may represent a new class of slip events, whose characteristics do not fit current understandings of slow or regular earthquakes.
-
more » « less
Abstract The magnitude‐frequency distribution (MFD) describes the relative proportion of earthquake magnitudes and provides vital information for seismic hazard assessment. The
b ‐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 theM 5.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
Abstract Machine learning can predict the timing and magnitude of laboratory earthquakes using statistics of acoustic emissions. The evolution of acoustic energy is critical for lab earthquake prediction; however, the connections between acoustic energy and fault zone processes leading to failure are poorly understood. Here, we document in detail the temporal evolution of acoustic energy during the laboratory seismic cycle. We report on friction experiments for a range of shearing velocities, normal stresses, and granular particle sizes. Acoustic emission data are recorded continuously throughout shear using broadband piezo‐ceramic sensors. The coseismic acoustic energy release scales directly with stress drop and is consistent with concepts of frictional contact mechanics and time‐dependent fault healing. Experiments conducted with larger grains (10.5 μm) show that the temporal evolution of acoustic energy scales directly with fault slip rate. In particular, the acoustic energy is low when the fault is locked and increases to a maximum during coseismic failure. Data from traditional slide‐hold‐slide friction tests confirm that acoustic energy release is closely linked to fault slip rate. Furthermore, variations in the true contact area of fault zone particles play a key role in the generation of acoustic energy. Our data show that acoustic radiation is related primarily to breaking/sliding of frictional contact junctions, which suggests that machine learning‐based laboratory earthquake prediction derives from frictional weakening processes that begin very early in the seismic cycle and well before macroscopic failure.
-
more » « less
Abstract The Húsavík‐Flatey Fault Zone (HFFZ) is the largest strike‐slip fault in Iceland and poses a high seismic risk to coastal communities. To investigate physics‐based constraints on earthquake hazards, we construct three fault system models of varying geometric complexity and model 79 3‐D multi‐fault dynamic rupture scenarios in the HFFZ. By assuming a simple regional prestress and varying hypocenter locations, we analyze the rupture dynamics, fault interactions, and the associated ground motions up to 2.5 Hz. All models account for regional seismotectonics, topo‐bathymetry, 3‐D subsurface velocity, viscoelastic attenuation, and off‐fault plasticity, and we explore the effect of fault roughness. The rupture scenarios obey earthquake scaling relations and predict magnitudes comparable to those of historical events. We show how fault system geometry and segmentation, hypocenter location, and prestress can affect the potential for rupture cascading, leading to varying slip distributions across different portions of the fault system. Our earthquake scenarios yield spatially heterogeneous near‐field ground motions modulated by geometric complexities, topography, and rupture directivity, particularly in the near‐field. The average ground motion attenuation characteristics of dynamic rupture scenarios of comparable magnitudes and mean stress drop are independent of variations in source complexity, magnitude‐consistent and in good agreement with the latest regional empirical ground motion models. However, physics‐based ground motion variability changes considerably with fault‐distance and increases for unilateral compared to bilateral ruptures. Systematic variations in physics‐based near‐fault ground motions provide important insights into the mechanics and potential earthquake hazard of large strike‐slip fault systems, such as the HFFZ.
