Exceedance Probabilities for Large Earthquakes from DIY Local Earthquake Ensemble Nowcasting and Forecasting: Magnitude, Natural Time, and Calendar Time

In this paper, we describe a method for computing calendar time forecasts in a local area for large earthquakes of a target magnitude MT using a count of small earthquakes in the magnitude range MS to MT in the area. Using the idea that the Gutenberg–Richter (GR) relation is valid throughout the surrounding region, we define an ensemble of earthquakes in larger surrounding regions to be used in computing the forecast. What follows is simple data mining. “Local” is defined by the probability of a large earthquake occurring within a defined circle of arbitrary radius surrounding a point of interest. The main (and for that matter, the only) assumption for all these works is that the GR magnitude–frequency relation holds. The method has significant skill, as defined by the Receiver Operating Characteristic (ROC) test, which improves as the time since the last major earthquake increases. The probability is conditioned on the number of small earthquakes n(t), with M ≥ MS = 3.49, that have occurred since the last large earthquake. The probability is computed directly as the Positive Predictive Value (PPV) associated with the ROC curve. The method is compared with the UCERF3 forecasts for the UCERF3-defined geographic boxes centered on Los Angeles and San Francisco and serves as an indicative benchmark. The method is then applied to a 125 km radius circular area around Los Angeles, California, following the 17 January 1994 magnitude M6.7 Northridge earthquake, and short-term forecasts (1-year and 5-year) are computed. We further apply the method to six additional geographic regions with validation by comparison with an estimate of the time-independent conditional Poisson probability. These regions are Athens, Greece; Chengdu, China; Jakarta, Indonesia; Lima, Peru; Santiago, Chile; and Tangshan, China.