Search for: All records

Seismology is witnessing explosive growth in the diversity and scale of earthquake catalogs. A key motivation for this community effort is that more data should translate into better earthquake forecasts. Such improvements are yet to be seen. Here, we introduce the Recurrent Earthquake foreCAST (RECAST), a deep‐learning model based on recent developments in neural temporal point processes. The model enables access to a greater volume and diversity of earthquake observations, overcoming the theoretical and computational limitations of traditional approaches. We benchmark against a temporal Epidemic Type Aftershock Sequence model. Tests on synthetic data suggest that with a modest‐sized data set, RECAST accurately models earthquake‐like point processes directly from cataloged data. Tests on earthquake catalogs in Southern California indicate improved fit and forecast accuracy compared to our benchmark when the training set is sufficiently long (>10⁴events). The basic components in RECAST add flexibility and scalability for earthquake forecasting without sacrificing performance.

 

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. 

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. 

The number of aftershocks increases with mainshock size following a well‐defined scaling law. However, excursions from the average behavior are common. This variability is particularly concerning for large earthquakes where the number of aftershocks varies by factors of 100 for mainshocks of comparable magnitude. Do observable factors lead to differences in aftershock behavior? We examine aftershock productivity relative to the global average for all mainshocks () from 1990 to 2019. A global map of earthquake productivity highlights the influence of tectonic regimes. Earthquake depth, lithosphere age, and plate boundary type correspond well with earthquake productivity. We investigate the role of mainshock attributes by compiling source dimensions, radiated seismic energy, stress drop, and a measure of slip heterogeneity based on finite‐fault source inversions for the largest earthquakes from 1990 to 2017. On an individual basis, stress drop, normalized rupture width, and aspect ratio most strongly correlate with aftershock productivity. A multivariate analysis shows that a particular set of parameters (dip, lithospheric age, and normalized rupture area) combines well to improve predictions of aftershock productivity on a cross‐validated data set. Our overall analysis is consistent with a model in which the volumetric abundance of nearby stressed faults controls the aftershock productivity rather than variations in source stress. Thus, we suggest a complementary approach to aftershock forecasts based on geological and rupture properties rather than local calibration alone.

 

On 4 and 6 July 2019, two large strike‐slip earthquakes withW‐phase moment magnitudesM_WW6.5 (foreshock) andM_WW7.1 (mainshock) struck the Eastern California Shear Zone, northeast of Ridgecrest. The faulting geometry and kinematic coseismic slip distribution of both events are determined by jointly inverting seismological and geodetic observations guided by aftershock and surface rupture locations. The foreshock ruptured two orthogonal faults with a prominent L‐shaped geometry with maximum slip of ~1.1 m on the NE‐SW segment. The mainshock faulting extended NW‐SE along several primary fault segments that straddle the foreshock slip. The surface rupture and slip model indicate mostly near‐horizontal strike‐slip motion with maximum slip of ~3.7 m, but there is a localized vertical dip‐slip motion. Both the foreshock and mainshock ruptures terminate in regions of complex surface offsets. High aftershock productivity and low rupture velocity may be the result of rupture of a relatively immature fault system.