Search

Earthquake nowcasting has been proposed as a means of tracking the change in large earthquake potential in a seismically active area. The method was developed using observable seismic data, in which probabilities of future large earthquakes can be computed using Receiver Operating Characteristic methods. Furthermore, analysis of the Shannon information content of the earthquake catalogs has been used to show that there is information contained in the catalogs, and that it can vary in time. So an important question remains, where does the information originate? In this paper, we examine this question using stochastic simulations of earthquake catalogs. Our catalog simulations are computed using an Earthquake Rescaled Aftershock Seismicity (“ERAS”) stochastic model. This model is similar in many ways to other stochastic seismicity simulations, but has the advantage that the model has only 2 free parameters to be set, one for the aftershock (Omori-Utsu) time decay, and one for the aftershock spatial migration away from the epicenter. Generating a simulation catalog and fitting the two parameters to the observed catalog such as California takes only a few minutes of wall clock time. While clustering can arise from random, Poisson statistics, we show that significant information in the simulation catalogs arises from the “non-Poisson” power-law aftershock clustering, implying that the practice of de-clustering observed catalogs may remove information that would otherwise be useful in forecasting and nowcasting. We also show that the nowcasting method provides similar results with the ERAS model as it does with observed seismicity.

The last 12 months have provided further evidence of the potential for cascading ecological and socio-political crises that were warned of 12 months ago. Then a consensus statement from the Regional Action on Climate Change Symposium warned: "the Earth's climatic, ecological, and human systems are converging towards a crisis that threatens to engulf global civilization within the lifetimes of children now living." Since then, the consequences of a broad set of extreme climate events (notably droughts, floods, and fires) have been compounded by interaction with impacts from multiple pandemics (including COVID-19 and cholera) and the Russia-Ukraine war. As a result, new connections are becoming visible between climate change and human health, large vulnerable populations are experiencing food crises, climate refugees are on the move, and the risks of water, food, and climate disruption have been visibly converging and compounding. Many vulnerable populations now face serious challenges to adapt. In light of these trends, this year, RACC identifies a range of measures to be taken at global and regional levels to bolster the resilience of these populations in the face of such emerging crises. In particular, at all scales, there is a need for globally available local data, reliable analytic techniques, community capacity to plan adaptation strategies, and the resources (scientific, technical, cultural, and economic) to implement them. To date, the rate of growth of the support for climate change resilience lags behind the rapid growth of cascading and converging risks. As an urgent message to COP27, it is proposed that the time is now right to devote much greater emphasis, global funding, and support to the increasing adaptation needs of vulnerable populations.

Nowcasting is a term originating from economics, finance, and meteorology. It refers to the process of determining the uncertain state of the economy, markets or the weather at the current time by indirect means. In this paper, we describe a simple two-parameter data analysis that reveals hidden order in otherwise seemingly chaotic earthquake seismicity. One of these parameters relates to a mechanism of seismic quiescence arising from the physics of strain-hardening of the crust prior to major events. We observe an earthquake cycle associated with major earthquakes in California, similar to what has long been postulated. An estimate of the earthquake hazard revealed by this state variable time series can be optimized by the use of machine learning in the form of the Receiver Operating Characteristic skill score. The ROC skill is used here as a loss function in a supervised learning mode. Our analysis is conducted in the region of 5° × 5° in latitude-longitude centered on Los Angeles, a region which we used in previous papers to build similar time series using more involved methods (Rundle & Donnellan, 2020, https://doi.org/10.1029/2020EA001097; Rundle, Donnellan et al., 2021, https://doi.org/10.1029/2021EA001757; Rundle, Stein et al., 2021, https://doi.org/10.1088/1361-6633/abf893). Here we show that not only does the state variable time series have forecast skill, the associated spatial probability densities have skill as well. In addition, use of the standard ROC and Precision (PPV) metrics allow probabilities of current earthquake hazard to be defined in a simple, straightforward, and rigorous way.

Abstract: The earthquake cycle of stress accumulation and release is associated with the elastic rebound hypothesis proposed by H.F. Reid following the M7.9 San Francisco earthquake of 1906. However, observing details of the actual values of time- and space-dependent tectonic stress is not possible at the present time. In two previous papers, we have proposed methods to image the earthquake cycle in California by means of proxy variables. These variables are based on correlations in patterns of small earthquakes that occur nearly continuously in time. The purpose of the present paper is to compare these two methods by evaluating their information content using decision thresholds and Receiver Operating Characteristic methods together with Shannon information entropy. Using seismic data from 1940 to present in California, we find that both methods provide nearly equivalent information on the rise and fall of earthquake correlations associated with major earthquakes in the region. We conclude that the resulting timeseries can be viewed as proxies for the cycle of stress accumulation and release associated with major tectonic activity. Article Highlights: The current state of the earthquake cycle of tectonic stress accumulation and release is unobservableWe review two methods for visualizing the current state of the earthquake cycle from correlation in small earthquake patternsMachine learning techniques indicate that signals in a correlation time series corresponding to future large earthquakes can be detected

GeoGateway (http://geo-gateway.org) is a web-based interface for analysis and modeling of geodetic imaging data and to support response to related disasters. Geodetic imaging data product currently supported by GeoGateway include Global Navigation Satellite System (GNSS) daily position time series and derived velocities and displacements and airborne Interferometric Synthetic Aperture Radar (InSAR) from NASA's UAVSAR platform. GeoGateway allows users to layer data products in a web map interface and extract information from various tools. Extracted products can be downloaded for further analysis. GeoGateway includes overlays of California fault traces, seismicity from user selected search parameters, and user supplied map files. GeoGateway also provides earthquake nowcasts and hazard maps as well as products created for related response to natural disasters. A user guide is present in the GeoGateway interface. The GeoGateway development team is also growing the user base through workshops, webinars, and video tutorials. GeoGateway is used in the classroom and for research by experts and non-experts including by students.

A promising tsunami warning system outlined here links Total Electron Content (TEC) disturbances in the ionosphere detected by GNSS satellites to maps of inundated areas along the affected coastlines. This will be accomplished through a pipeline of earthquake, ionosphere, and tsunami simulations. The focus of this paper is to demonstrate the feasibility and validity of our tsunami simulator for use in this simulation pipeline. The simulator we have developed is called Tsunami Squares, a computationally inexpensive program that employs a unique method of propagating water according to shallow water theory. Unlike other wave simulators, Tsunami Squares can easily propagate waves from sea to land, allowing flooded areas along the coast to be mapped. To validate its accuracy, a comparison is made to the Regional Ocean Modeling System tsunami simulator, which solves the Navier Stokes equations directly using a finite-difference method. The two simulation techniques show good agreement when compared over the 2011 Japan event and the 2010 Chile event. In addition, examples of inundation plots are displayed here to demonstrate its capability for such a warning system.

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.

We present a data-driven approach to clustering or grouping Global Navigation Satellite System (GNSS) stations according to observed velocities, displacements or other selected characteristics. Clustering GNSS stations provides useful scientific information, and is a necessary initial step in other analysis, such as detecting aseismic transient signals (Granat et al., 2013, https://doi.org/10.1785/0220130039). Desired features of the data can be selected for clustering, including some subset of displacement or velocity components, uncertainty estimates, station location, and other relevant information. Based on those selections, the clustering procedure autonomously groups the GNSS stations according to a selected clustering method. We have implemented this approach as a Python application, allowing us to draw upon the full range of open source clustering methods available in Pythons scikit-learn package (Pedregosa et al., 2011, https://doi.org/10.5555/1953048.2078195). The application returns the stations labeled by group as a table and color coded KML file and is designed to work with the GNSS information available from GeoGateway (Donnellan et al., 2021, https://doi.org/10.1007/s12145-020-00561-7; Heflin et al., 2020, https://doi.org/10.1029/2019ea000644) but is easily extensible. We demonstrate the methodology on California and western Nevada. The results show partitions that follow faults or geologic boundaries, including for recent large earthquakes and post-seismic motion. The San Andreas fault system is most prominent, reflecting Pacific-North American plate boundary motion. Deformation reflected as class boundaries is distributed north and south of the central California creeping section. For most models a cluster boundary connects the southernmost San Andreas fault with the Eastern California Shear Zone (ECSZ) rather than continuing through the San Gorgonio Pass.