Nowcasting Earthquakes: Imaging the Earthquake Cycle in California With Machine Learning.

We propose a new machine learning‐based method for nowcasting earthquakes to image the time‐dependent earthquake cycle. The result is a timeseries that may correspond to the process of stress accumulation and release. The timeseries are constructed by using principal component analysis of regional seismicity. The patterns are found as eigenvectors of the cross‐correlation matrix of a collection of seismicity timeseries in a coarse grained regional spatial grid (pattern recognition via unsupervised machine learning). The eigenvalues of this matrix represent the relative importance of the various eigenpatterns. Using the eigenvectors and eigenvalues, we compute the weighted correlation timeseries of the regional seismicity. This timeseries has the property that the weighted correlation generally decreases prior to major earthquakes in the region, and increases suddenly just after a major earthquake occurs. As in a previous paper (Rundle & Donnellan, 2020, https://doi.org/10.1029/2020ea001097), we find that this method produces a nowcasting timeseries that resembles the hypothesized regional stress accumulation and release process characterizing the earthquake cycle. We then address the problem of whether the timeseries contain information regarding future large earthquakes. For this, we compute a receiver operating characteristic and determine the decision thresholds for several future time periods of interest (optimization via supervised machine learning). We find that signals can be detected that can be used to characterize the information content of the timeseries. These signals may be useful in assessing present and near‐future seismic hazards. Plain Language Summary: Major earthquakes on fault systems in a tectonically active region are thought to occur in approximately repetitive cycles as a result of the buildup and release of tectonic forces (stress). Nowcasting is a technique adopted from weather, finance, and other fields that use readily observable proxy data to represent the unobservable stress accumulation process of interest. This paper presents a method that computes a timeseries representing the weighted correlation of small earthquake activity in the California region from 1950 to 2020. Prior to major magnitude M > 7 earthquakes, the timeseries trends toward lower values. Just after the earthquake occurs, the timeseries increases suddenly in association with the earthquake, before resuming its gradual trend toward lower values. Plotting the timeseries on an inverted scale, one sees a cyclic behavior that strongly resembles the hypothesized earthquake cycle. In principle, we can therefore use this timeseries for nowcasting, as a proxy for stress accumulation and release. Using methods of signal detection first developed for radar by the British in the 1940's, we find that the timeseries contain information about future large earthquakes that can be used for hazard assessment. Key Points: The current state of the earthquake cycle of tectonic stress accumulation and release is unobservable with existing methodsWe show that readily observable small earthquake correlations can be used to nowcast the current state of the earthquake cycleMachine learning techniques indicate that signals corresponding to future large earthquakes can be detected in a correlation time series [ABSTRACT FROM AUTHOR]