Abstract Gulia and Wiemer (2019; hereafter, GW2019) proposed a near-real-time monitoring system to discriminate between foreshocks and aftershocks. Our analysis (Dascher-Cousineau et al., 2020; hereinater, DC2020) tested the sensitivity of the proposed Foreshock Traffic-Light System output to parameter choices left to expert judgment for the 2019 Ridgecrest Mw 7.1 and 2020 Puerto Rico Mw 6.4 earthquake sequences. In the accompanying comment, Gulia and Wiemer (2021) suggest that at least six different methodological deviations lead to different pseudoprospective warning levels, particularly for the Ridgecrest aftershock sequence which they had separately evaluated. Here, we show that for four of the six claimed deviations, we conformed to the criteria outlined in GW2019. Two true deviations from the defined procedure are clarified and justified here. We conclude as we did originally, by emphasizing the influence of expert judgment on the outcome in the analysis.
more »
« less
Two Foreshock Sequences Post Gulia and Wiemer (2019)
Abstract Recognizing earthquakes as foreshocks in real time would provide a valuable forecasting capability. In a recent study, Gulia and Wiemer (2019) proposed a traffic-light system that relies on abrupt changes in b-values relative to background values. The approach utilizes high-resolution earthquake catalogs to monitor localized regions around the largest events and distinguish foreshock sequences (reduced b-values) from aftershock sequences (increased b-values). The recent well-recorded earthquake foreshock sequences in Ridgecrest, California, and Maria Antonia, Puerto Rico, provide an opportunity to test the procedure. For Ridgecrest, our b-value time series indicates an elevated risk of a larger impending earthquake during the Mw 6.4 foreshock sequence and provides an ambiguous identification of the onset of the Mw 7.1 aftershock sequence. However, the exact result depends strongly on expert judgment. Monte Carlo sampling across a range of reasonable decisions most often results in ambiguous warning levels. In the case of the Puerto Rico sequence, we record significant drops in b-value prior to and following the largest event (Mw 6.4) in the sequence. The b-value has still not returned to background levels (12 February 2020). The Ridgecrest sequence roughly conforms to expectations; the Puerto Rico sequence will only do so if a larger event occurs in the future with an ensuing b-value increase. Any real-time implementation of this approach will require dense instrumentation, consistent (versioned) low completeness catalogs, well-calibrated maps of regionalized background b-values, systematic real-time catalog production, and robust decision making about the event source volumes to analyze.
more »
« less
- NSF-PAR ID:
- 10197709
- Date Published:
- Journal Name:
- Seismological Research Letters
- Volume:
- 91
- Issue:
- 5
- ISSN:
- 0895-0695
- Page Range / eLocation ID:
- 2843 to 2850
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)ABSTRACT The 2019 Ridgecrest earthquake sequence culminated in the largest seismic event in California since the 1999 Mw 7.1 Hector Mine earthquake. Here, we combine geodetic and seismic data to study the rupture process of both the 4 July Mw 6.4 foreshock and the 6 July Mw 7.1 mainshock. The results show that the Mw 6.4 foreshock rupture started on a northwest-striking right-lateral fault, and then continued on a southwest-striking fault with mainly left-lateral slip. Although most moment release during the Mw 6.4 foreshock was along the southwest-striking fault, slip on the northwest-striking fault seems to have played a more important role in triggering the Mw 7.1 mainshock that happened ∼34 hr later. Rupture of the Mw 7.1 mainshock was characterized by dominantly right-lateral slip on a series of overall northwest-striking fault strands, including the one that had already been activated during the nucleation of the Mw 6.4 foreshock. The maximum slip of the 2019 Ridgecrest earthquake was ∼5 m, located at a depth range of 3–8 km near the Mw 7.1 epicenter, corresponding to a shallow slip deficit of ∼20%–30%. Both the foreshock and mainshock had a relatively low-rupture velocity of ∼2 km/s, which is possibly related to the geometric complexity and immaturity of the eastern California shear zone faults. The 2019 Ridgecrest earthquake produced significant stress perturbations on nearby fault networks, especially along the Garlock fault segment immediately southwest of the 2019 Ridgecrest rupture, in which the coulomb stress increase was up to ∼0.5 MPa. Despite the good coverage of both geodetic and seismic observations, published coseismic slip models of the 2019 Ridgecrest earthquake sequence show large variations, which highlight the uncertainty of routinely performed earthquake rupture inversions and their interpretation for underlying rupture processes.more » « less
-
more » « less
Abstract Foreshocks are the only currently widely identified precursory seismic behavior, yet their utility and even identifiability are problematic, in part because of extreme variation in behavior. Here, we establish some global trends that help identify the expected frequency of foreshocks as well the type of earthquake most prone to foreshocks. We establish these tendencies using the global earthquake catalog of the U.S. Geological Survey National Earthquake Information Center with a completeness level of magnitude 5 and mainshocks with Mw≥7.0. Foreshocks are identified using three clustering algorithms to address the challenge of distinguishing foreshocks from background activity. The methods give a range of 15%–43% of large mainshocks having at least one foreshock but a narrower range of 13%–26% having at least one foreshock with magnitude within two units of the mainshock magnitude. These observed global foreshock rates are similar to regional values for a completeness level of magnitude 3 using the same detection conditions. The foreshock sequences have distinctive characteristics with the global composite population b-values being lower for foreshocks than for aftershocks, an attribute that is also manifested in synthetic catalogs computed by epidemic-type aftershock sequences, which intrinsically involves only cascading processes. Focal mechanism similarity of foreshocks relative to mainshocks is more pronounced than for aftershocks. Despite these distinguishing characteristics of foreshock sequences, the conditions that promote high foreshock productivity are similar to those that promote high aftershock productivity. For instance, a modestly higher percentage of interplate mainshocks have foreshocks than intraplate mainshocks, and reverse faulting events slightly more commonly have foreshocks than normal or strike-slip-faulting mainshocks. The western circum-Pacific is prone to having slightly more foreshock activity than the eastern circum-Pacific.
-
more » « less
SUMMARY The spatio-temporal properties of seismicity give us incisive insight into the stress state evolution and fault structures of the crust. Empirical models based on self-exciting point processes continue to provide an important tool for analysing seismicity, given the epistemic uncertainty associated with physical models. In particular, the epidemic-type aftershock sequence (ETAS) model acts as a reference model for studying seismicity catalogues. The traditional ETAS model uses simple parametric definitions for the background rate of triggering-independent seismicity. This reduces the effectiveness of the basic ETAS model in modelling the temporally complex seismicity patterns seen in seismic swarms that are dominated by aseismic tectonic processes such as fluid injection rather than aftershock triggering. In order to robustly capture time-varying seismicity rates, we introduce a deep Gaussian process (GP) formulation for the background rate as an extension to ETAS. GPs are a robust non-parametric model for function spaces with covariance structure. By conditioning the length-scale structure of a GP with another GP, we have a deep-GP: a probabilistic, hierarchical model that automatically tunes its structure to match data constraints. We show how the deep-GP-ETAS model can be efficiently sampled by making use of a Metropolis-within-Gibbs scheme, taking advantage of the branching process formulation of ETAS and a stochastic partial differential equation (SPDE) approximation for Matérn GPs. We illustrate our method using synthetic examples, and show that the deep-GP-ETAS model successfully captures multiscale temporal behaviour in the background forcing rate of seismicity. We then apply the results to two real-data catalogues: the Ridgecrest, CA 2019 July 5 Mw 7.1 event catalogue, showing that deep-GP-ETAS can successfully characterize a classical aftershock sequence; and the 2016–2019 Cahuilla, CA earthquake swarm, which shows two distinct phases of aseismic forcing concordant with a fluid injection-driven initial sequence, arrest of the fluid along a physical barrier and release following the largest Mw 4.4 event of the sequence.
-
Abstract The recent Indios, Puerto Rico earthquake sequence has drawn attention, as the increased seismicity rate in this area was unprecedented. The sequence began on 28 December 2019, caused a 6.4 magnitude earthquake on 7 January 2020, and remained active over a year later. This sequence fits the nominal definition of an earthquake swarm in that it had an abrupt onset, a sustained high rate of seismicity without a clear triggering mainshock or evidence for Omori decay, and a lack of adherence to Bath’s law. However, the sequence also had several prominent mainshock–aftershock (MS–AS) sequences embedded within it. We applied three-station waveform cross correlation to the early part of this sequence using the Puerto Rico Seismic Network (PRSN) catalog as templates, which confirmed the mixture of swarm and MS–AS patterns. In an effort to place this intriguing sequence in the context of the previous seismicity in Puerto Rico, we investigated the existence of swarms and MS–AS sequences recorded by the PRSN since 1987 by identifying sequences with increased seismicity rate when compared to the background rate. About 59 sequences were manually verified and characterized into swarms or MS–AS. We found that 58% of the sequences follow traditional swarm patterns and 14% adhere to traditional MS–AS behavior, whereas 29% of the sequences have a mixture of both swarm and MS–AS behaviors. These findings suggest that it is not unusual for the Indios sequence to have a mixture of both the characteristics. In addition, the detection of many swarms distributed over a broad area of the subduction interface indicates stress heterogeneity and low-coupling consistent with prior studies indicating that the potential for a magnitude ∼8 megathrust earthquake along the Puerto Rico trench is unlikely.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1785/0220200082
