Search

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.

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.

Interferometric synthetic aperture radar (InSAR) interferograms contain valuable information about the fault systems hidden beneath the surface of the Earth. In a new approach, we aim to fit InSAR ground deformation data using a distribution of multiple seismic point sources whose parameters are found by a genetic algorithm. The resulting source distribution could provide another useful tool in solving the difficult problem of accurately mapping earthquake faults. We apply the algorithm to an ALOS-2 InSAR interferogram and perform a multifractal analysis on the resulting distribution, finding that it exhibits multifractal properties. We report first results and discuss advantages and disadvantages of this approach.

We use UAVSAR interferograms to characterize fault slip, triggered by the Mw 7.2 El Mayor-Cucapah earthquake on the 1 San Andreas Fault in the Coachella Valley providing comprehensive maps of short-term geodetic surface deformation that complement in situ measurements. Creepmeters and geological mapping of fault offsets on Durmid Hill recorded 4 and 8 mm of average triggered slip respectively on the fault, in contrast to radar views that reveal significant off-fault dextral deformation averaging 20 mm. Unlike slip in previous triggered slip events on the southernmost San Andreas fault, dextral shear in 2010 is not confined to transpressional hills in the Coachella valley. Edge detection and gradient estimation applied to the 50-m-sampled interferogram data identify the location (to 20 m) and local strike (to

Seismic nowcasting uses counts of small earthquakes as proxy data to estimate the current dynamical state of an earthquake fault system. The result is an earthquake potential score that characterizes the current state of progress of a defined geographic region through its nominal earthquake "cycle." The count of small earthquakes since the last large earthquake is the natural time that has elapsed since the last large earthquake (Varotsos et al., 2006, https://doi.org/10.1103/PhysRevE.74.021123). In addition to natural time, earthquake sequences can also be analyzed using Shannon information entropy ("information"), an idea that was pioneered by Shannon (1948, https://doi.org/10.1002/j.1538-7305.1948.tb01338.x). As a first step to add seismic information entropy into the nowcasting method, we incorporate magnitude information into the natural time counts by using event self-information. We find in this first application of seismic information entropy that the earthquake potential score values are similar to the values using only natural time. However, other characteristics of earthquake sequences, including the interevent time intervals, or the departure of higher magnitude events from the magnitude-frequency scaling line, may contain additional information.

Abstract 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.

Tectonic motion across the Los Angeles region is distributed across an intricate network of strike-slip and thrust faults that will be released in destructive earthquakes similar to or larger than the 1933 M6.4 Long Beach and 1994 M6.7 Northridge events. Here we show that Los Angeles regional thrust, strike-slip, and oblique faults are connected and move concurrently with measurable surface deformation, even in moderate magnitude earthquakes, as part of a fault system that accommodates north-south shortening and westerly tectonic escape of northern Los Angeles. The 28 March 2014 M5.1 La Habra earthquake occurred on a northeast striking, northwest dipping left-lateral oblique thrust fault northeast of Los Angeles. We present crustal deformation observation spanning the earthquake showing that concurrent deformation occurred on several structures in the shallow crust. The seismic moment of the earthquake is 82% of the total geodetic moment released. Slip within the unconsolidated upper sedimentary layer may reflect shallow release of accumulated strain on still-locked deeper structures. A future M6.1-6.3 earthquake would account for the accumulated strain. Such an event could occur on any one or several of these faults, which may not have been identified by geologic surface mapping.

Abstract 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.

Abstract 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 Python's 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.

Abstract Interferometric synthetic aperture radar (InSAR) interferograms contain valuable information about the fault systems hidden beneath the surface of the Earth. In a new approach, we aim to fit InSAR ground deformation data using a distribution of multiple seismic point sources whose parameters are found by a genetic algorithm. The resulting source distribution could provide another useful tool in solving the difficult problem of accurately mapping earthquake faults. We apply the algorithm to an ALOS‐2 InSAR interferogram and perform a multifractal analysis on the resulting distribution, finding that it exhibits multifractal properties. We report first results and discuss advantages and disadvantages of this approach.

Abstract We use UAVSAR interferograms to characterize fault slip, triggered by the Mw 7.2 El Mayor‐Cucapah earthquake on the 1 San Andreas Fault in the Coachella Valley providing comprehensive maps of short‐term geodetic surface deformation that complement in situ measurements. Creepmeters and geological mapping of fault offsets on Durmid Hill recorded 4 and 8 mm of average triggered slip respectively on the fault, in contrast to radar views that reveal significant off‐fault dextral deformation averaging 20 mm. Unlike slip in previous triggered slip events on the southernmost San Andreas fault, dextral shear in 2010 is not confined to transpressional hills in the Coachella valley. Edge detection and gradient estimation applied to the 50‐m‐sampled interferogram data identify the location (to 20 m) and local strike (to

Abstract Because cross‐polarized radar returns are highly associated with volume scatter, radar polarimetry returns tend to show strong evidence of wildfire scars and recovery in forest and chaparral. We focus on the polarimetry images from UAVSAR (PolSAR) line SanAnd_08525, which covers a roughly 20 km wide swath over the Transverse Range including parts of the Santa Monica, San Gabriel and San Bernardino Mountains. We select images from four acquisition dates from October 2009 to September 2020, very roughly 4 years apart. These are compared to fire perimeters from the national Geospatial Multi‐Agency Coordination and NIFC databases for years 2003–2020, which shows the areas affected by the major fires (west to east) Springs2013, Woolsey2018, Topanga2005, LaTuna2017, Station2009, BlueCut2016, Pilot2016, Slide2007, Butler2007, and many smaller fires. PolSAR images are shown to be helpful in identifying types and boundaries of fire, 50‐meter scale details of vegetation loss, and variability of vegetation recovery in post‐fire years.

Key Point The editors thank the 2020 peer reviewers

Abstract Seismic bursts in Southern California are sequences of small earthquakes strongly clustered in space and time and include seismic swarms and aftershock sequences. A readily observable property of these events, the radius of gyration (RG), allows us to connect the bursts to the temporal occurrence of the largest M ≥ 7 earthquakes in California since 1984. In the Southern California earthquake catalog, we identify hundreds of these potentially coherent space‐time structures in a region defined by a circle of radius 600 km around Los Angeles. We compute RG for each cluster then filter them to identify those bursts with large numbers of events closely clustered in space, which we call “compact” bursts. Our basic assumption is that these compact bursts reflect the dynamics associated with large earthquakes. Once we have filtered the burst catalog, we apply an exponential moving average to construct a time series for the Southern California region. We observe that the RG of these bursts systematically decreases prior to large earthquakes, in a process that we might term “radial localization.” The RG then rapidly increases during an aftershock sequence, and a new cycle of “radial localization” then begins. These time series display cycles of recharge and discharge reminiscent of seismic stress accumulation and release in the elastic rebound process. The complex burst dynamics we observe are evidently a property of the region as a whole, rather than being associated with individual faults. This new method allows us to improve earthquake nowcasting, which is a technique to evaluate the current state of hazard in a seismically active region.

Abstract This paper describes the methods used to estimate positions, velocities, breaks, and seasonal terms from daily Global Navigation Satellite System (GNSS) measurements. Break detection and outlier removal have been automated so that decades of daily measurements from thousands of stations can be processed in a few hours. New measurements are added, and parameters are updated every week. Model parameters allow separation of interseismic, annual, coseismic, and postseismic signals. Tools available through GeoGateway (http://geo-gateway.org) allow rapid visualization and analysis of these terms for results that can be subsetted in time or space. Results show highly variable and nonlinear motion for GPS stations in southern California. The variable motion is related to seasonal motions, distributed tectonic motion, earthquakes, and postseismic motions that can continue for years. In some areas results suggest that additional processes are responsible for the observed motions. In general, following earthquakes, stations return to their long‐term motions after 2–3 years, though some exceptions occur. The use of the tools shows nonlinear motion in the Salton Trough of southern California related to the 2010 M7.2 El Mayor‐Cucapah earthquake, 2012 Brawley earthquake swarm, and a creep event on the Superstition Hills fault in 2017.

Key Points The editors thank the 2019 peer reviewers