Ensembles of ETAS Models Provide Optimal Operational Earthquake Forecasting During Swarms: Insights from the 2015 San Ramon, California Swarm

2019 ◽  
Vol 109 (6) ◽  
pp. 2145-2158 ◽  
Author(s):  
Andrea L. Llenos ◽  
Andrew J. Michael

Abstract Earthquake swarms, typically modeled as time‐varying changes in background seismicity, which are driven by external processes such as fluid flow or aseismic creep, present challenges for operational earthquake forecasting. Although the time decay of aftershock sequences can be estimated with the modified Omori law, it is difficult to forecast the temporal behavior of seismicity rates during a swarm. To explore these issues, we apply the epidemic‐type aftershock sequence (ETAS) model (Ogata, 1988) to the 2015 San Ramon, California swarm, which lasted several weeks and had almost 100 2≤M≤3.6 earthquakes. We develop three‐day forecasts during the swarm based on an ETAS model fit to all prior seismicity in the region as well as an ETAS model fit only to previous swarms in the region, which is better at capturing the higher background rate during the swarm. We also explore forecasts in which the background rate is updated periodically during the swarm using data over different lookback windows and find generally these models perform better than the models in which the background rate is fixed. Finally, we construct ensemble forecasts by combining the different models weighted according to their performance. The ensemble forecasts outperform all of the individual models and allow us to avoid making arbitrary choices at the outset of a swarm as to which single model will perform the best.

2020 ◽  
Vol 91 (3) ◽  
pp. 1567-1578 ◽  
Author(s):  
Kevin R. Milner ◽  
Edward H. Field ◽  
William H. Savran ◽  
Morgan T. Page ◽  
Thomas H. Jordan

Abstract The first Uniform California Earthquake Rupture Forecast, Version 3–epidemic-type aftershock sequence (UCERF3-ETAS) aftershock simulations were running on a high-performance computing cluster within 33 min of the 4 July 2019 M 6.4 Searles Valley earthquake. UCERF3-ETAS, an extension of the third Uniform California Earthquake Rupture Forecast (UCERF3), is the first comprehensive, fault-based, epidemic-type aftershock sequence (ETAS) model. It produces ensembles of synthetic aftershock sequences both on and off explicitly modeled UCERF3 faults to answer a key question repeatedly asked during the Ridgecrest sequence: What are the chances that the earthquake that just occurred will turn out to be the foreshock of an even bigger event? As the sequence unfolded—including one such larger event, the 5 July 2019 M 7.1 Ridgecrest earthquake almost 34 hr later—we updated the model with observed aftershocks, finite-rupture estimates, sequence-specific parameters, and alternative UCERF3-ETAS variants. Although configuring and running UCERF3-ETAS at the time of the earthquake was not fully automated, considerable effort had been focused in 2018 on improving model documentation and ease of use with a public GitHub repository, command line tools, and flexible configuration files. These efforts allowed us to quickly respond and efficiently configure new simulations as the sequence evolved. Here, we discuss lessons learned during the Ridgecrest sequence, including sensitivities of fault triggering probabilities to poorly constrained finite-rupture estimates and model assumptions, as well as implications for UCERF3-ETAS operationalization.


2019 ◽  
Vol 219 (3) ◽  
pp. 2148-2164
Author(s):  
A M Lombardi

SUMMARY The operational earthquake forecasting (OEF) is a procedure aimed at informing communities on how seismic hazard changes with time. This can help them live with seismicity and mitigate risk of destructive earthquakes. A successful short-term prediction scheme is not yet produced, but the search for it should not be abandoned. This requires more research on seismogenetic processes and, specifically, inclusion of any information about earthquakes in models, to improve forecast of future events, at any spatio-temporal-magnitude scale. The short- and long-term forecast perspectives of earthquake occurrence followed, up to now, separate paths, involving different data and peculiar models. But actually they are not so different and have common features, being parts of the same physical process. Research on earthquake predictability can help to search for a common path in different forecast perspectives. This study aims to improve the modelling of long-term features of seismicity inside the epidemic type aftershock sequence (ETAS) model, largely used for short-term forecast and OEF procedures. Specifically, a more comprehensive estimation of background seismicity rate inside the ETAS model is attempted, by merging different types of data (seismological instrumental, historical, geological), such that information on faults and on long-term seismicity integrates instrumental data, on which the ETAS models are generally set up. The main finding is that long-term historical seismicity and geological fault data improve the pseudo-prospective forecasts of independent seismicity. The study is divided in three parts. The first consists in models formulation and parameter estimation on recent seismicity of Italy. Specifically, two versions of ETAS model are compared: a ‘standard’, previously published, formulation, only based on instrumental seismicity, and a new version, integrating different types of data for background seismicity estimation. Secondly, a pseudo-prospective test is performed on independent seismicity, both to test the reliability of formulated models and to compare them, in order to identify the best version. Finally, a prospective forecast is made, to point out differences and similarities in predicting future seismicity between two models. This study must be considered in the context of its limitations; anyway, it proves, beyond argument, the usefulness of a more sophisticated estimation of background rate, inside short-term modelling of earthquakes.


2019 ◽  
Vol 71 (1) ◽  
Author(s):  
Takao Kumazawa ◽  
Yosihiko Ogata ◽  
Hiroshi Tsuruoka

AbstractWe applied the epidemic type aftershock sequence (ETAS) model, the two-stage ETAS model and the non-stationary ETAS model to investigate the detailed features of the series of earthquake occurrences before and after the M6.7 Hokkaido Eastern Iburi earthquake on 6 September 2018, based on earthquake data from October 1997. First, after the 2003 M8.0 Tokachi-Oki earthquake, seismic activity in the Eastern Iburi region reduced relative to the ETAS model. During this period, the depth ranges of the seismicity were migrating towards shallow depths, where a swarm cluster, including a M5.1 earthquake, finally occurred in the deepest part of the range. This swarm activity was well described by the non-stationary ETAS model until the M6.7 main shock. The aftershocks of the M6.7 earthquake obeyed the ETAS model until the M5.8 largest aftershock, except for a period of several days when small, swarm-like activity was found at the southern end of the aftershock region. However, when we focus on the medium and larger aftershocks, we observed quiescence relative to the ETAS model from 8.6 days after the main shock until the M5.8 largest aftershock. For micro-earthquakes, we further studied the separated aftershock sequences in the naturally divided aftershock volumes. We found that the temporal changes in the background rate and triggering coefficient (aftershock productivity) in respective sub-volumes were in contrast with each other. In particular, relative quiescence was seen in the northern deep zones that includes the M5.8 largest aftershock. Furthermore, changes in the b-values of the whole aftershock activity showed an increasing trend with respect to the logarithm of elapsed time during the entire aftershock period, which is ultimately explained by the spatially different characteristics of the aftershocks.


Author(s):  
Yue Liu ◽  
Jiancang Zhuang ◽  
Changsheng Jiang

Abstract The aftershock zone of the 1976 Ms 7.8 Tangshan, China, earthquake remains seismically active, experiencing moderate events such as the 5 December 2019 Ms 4.5 Fengnan event. It is still debated whether aftershock sequences following large earthquakes in low-seismicity continental regions can persist for several centuries. To understand the current stage of the Tangshan aftershock sequence, we analyze the sequence record and separate background seismicity from the triggering effect using a finite-source epidemic-type aftershock sequence model. Our results show that the background rate notably decreases after the mainshock. The estimated probability that the most recent 5 December 2019 Ms 4.5 Fengnan District, Tangshan, earthquake is a background event is 50.6%. This indicates that the contemporary seismicity in the Tangshan aftershock zone can be characterized as a transition from aftershock activity to background seismicity. Although the aftershock sequence is still active in the Tangshan region, it is overridden by background seismicity.


2020 ◽  
Author(s):  
Christian Grimm ◽  
Martin Käser ◽  
Helmut Küchenhoff

<p>While Probabilistic Seismic Hazard Assessment is commonly based on earthquake catalogues in a declustered form, ongoing seismicity in aftershock sequences is known to be able to add significant hazard, which can also increase the damage potential to already affected structures in risk assessment. Especially so-called earthquake doublets (multiplets), i.e. a cluster mainshock being followed or preceded by one (or more) events with a similarly strong magnitude occurring within pre-defined temporal and spatial limits, can cause loss multiplication effects to the insurance industry, which therefore has a pronounced interest in investigating the frequency of earthquake doublets to happen worldwide. A widely used method to analyse and simulate the triggering process of earthquake sequences is the Epidemic Type Aftershock Sequence (ETAS) model. We estimate the ETAS model parameters for some regional areas and produce synthetic catalogues, which are then analysed particularly with respect to the occurrence of earthquake doublets and compared to the observed history. Also, different seismic subduction-type regions in the world are pointed out to have shown differing relative frequencies of earthquake doublets. Regression models are used to study whether certain mainshock and local, geophysical properties such as magnitude, dip and rake angle, depth, distance to subduction plate interface and velocity of converging subduction plates nearby show explanatory power for the probability of a cluster containing an earthquake doublet.</p>


2021 ◽  
Author(s):  
Christian Grimm ◽  
Sebastian Hainzl ◽  
Martin Käser ◽  
Helmut Küchenhoff

Abstract Strong earthquakes cause aftershock sequences that are clustered in time according to a power decay law, and in space along their extended rupture, shaping a typically elongate pattern of aftershock locations. A widely used approach to model seismic clustering is the Epidemic Type Aftershock Sequence (ETAS) model, that shows three major biases: First, the conventional ETAS approach assumes isotropic spatial triggering, which stands in conflict with observations and geophysical arguments for strong earthquakes. Second, the spatial kernel has unlimited extent, allowing smaller events to exert disproportionate trigger potential over an unrealistically large area. Third, the ETAS model assumes complete event records and neglects inevitable short-term aftershock incompleteness as a consequence of overlapping coda waves. These three effects can substantially bias the parameter estimation and particularly lead to underestimated cluster sizes. In this article, we combine the approach of Grimm (2021), which introduced a generalized anisotropic and locally restricted spatial kernel, with the ETAS-Incomplete (ETASI) time model of Hainzl (2021), to define an ETASI space-time model with flexible spatial kernel that solves the abovementioned shortcomings. We apply different model versions to a triad of forecasting experiments of the 2019 Ridgecrest sequence, and evaluate the prediction quality with respect to cluster size, largest aftershock magnitude and spatial distribution. The new model provides the potential of more realistic simulations of on-going aftershock activity, e.g.~allowing better predictions of the probability and location of a strong, damaging aftershock, which might be beneficial for short term risk assessment and desaster response.


2019 ◽  
Vol 71 (1) ◽  
Author(s):  
Yosihiko Ogata ◽  
Koich Katsura ◽  
Hiroshi Tsuruoka ◽  
Naoshi Hirata

Abstract We propose an extended 3D space (longitude, latitude, and depth) epidemic-type aftershock sequence (ETAS) model for seismicity forecasts beneath the greater Tokyo area (the Kanto region), which also takes into account the effects induced by the M9 Tohoku-Oki earthquake of 2011. The model is characterized by a number of 3D location-dependent parameters, such as the background seismicity rates, and the productivity rate induced by the Tohoku earthquake. These allow production of high-resolution predictive mappings in zones where hypocenters are densely populated. The optimally inverted 3D spatial images of the characterizing parameters effectively discriminate seismicity features in the crust and near the plate boundaries. The success of the model is demonstrated using short-, intermediate- and long-term probability forecasts of intermediate and large earthquake occurrences beneath the Kanto region.


2021 ◽  
Author(s):  
Jordi Baro

<p>Earthquake catalogs exhibit strong spatio-temporal correlations. As such, earthquakes are often classified into clusters of correlated activity. Clusters themselves are traditionally classified in two different kinds: (i) bursts, with a clear hierarchical structure between a single strong mainshock, preceded by a few foreshocks and followed by a power-law decaying aftershock sequence, and (ii) swarms, exhibiting a non-trivial activity rate that cannot be reduced to such a simple hierarchy between events. </p><p>The Epidemic Aftershock Sequence (ETAS) model is a linear Hawkes point process able to reproduce earthquake clusters from empirical statistical laws [Ogata, 1998]. Although not always explicit, the ETAS model is often interpreted as the outcome of a background activity driven by external forces and a Galton-Watson branching process with one-to-one causal links between events [Saichev et al., 2005]. Declustering techniques based on field observations [Baiesi & Paczuski, 2004] can be used to infer the most likely causal links between events in a cluster. Following this method, Zaliapin and Ben‐Zion (2013) determined the statistical properties of earthquake clusters characterizing bursts and swarms, finding a relationship between the predominant cluster-class and the heat flow in seismic regions.</p><p>Here, I show how the statistical properties of clusters are related to the fundamental statistics of the underlying seismogenic process, modeled in two point-process paradigms [Baró, 2020].</p><p>The classification of clusters into bursts and swarms appears naturally in the standard ETAS model with homogeneous rates and are determined by the average branching ratio (nb) and the ratio between exponents α and b characterizing the production of aftershocks and the distribution of magnitudes, respectively. The scale-free ETAS model, equivalent to the BASS model [Turcotte, et al., 2007], and usual in cold active tectonic regions, is imposed by α=b and reproduces bursts. In contrast, by imposing α<0.5b, we recover the properties of swarms, characteristic of regions with high heat flow. </p><p>Alternatively, the same declustering methodology applied to a non-homogeneous Poisson process with a non-factorizable intensity, i.e. in absence of causal links, recovers swarms with α=0, i.e. a Poisson Galton-Watson process, with similar statistical properties to the ETAS model in the regime α<0.5b.</p><p>Therefore, while bursts are likely to represent actual causal links between events, swarms can either denote causal links with low α/b ratio or variations of the background rate caused by exogenous processes introducing local and transient stress changes. Furthermore, the redundancy in the statistical laws can be used to test the hypotheses posed by the ETAS model as a memory‐less branching process. </p><p>References:</p><ul><li> <p>Baiesi, M., & Paczuski, M. (2004). <em>Physical Review E</em>, 69, 66,106. doi:10.1103/PhysRevE.69.066106.</p> </li> <li> <p>Baró, J. (2020).  <em>Journal of Geophysical Research: Solid Earth,</em> 125, e2019JB018530. doi:10.1029/2019JB018530.</p> </li> <li> <p>Ogata, Y. (1998) <em>Annals of the Institute of Statistical Mathematics,</em> 50(2), 379–402. doi:10.1023/A:1003403601725.</p> </li> <li> <p>Saichev, A., Helmstetter, A. & Sornette, D. (2005) <em>Pure appl. geophys.</em> 162, 1113–1134. doi:10.1007/s00024-004-2663-6.</p> </li> <li> <p>Turcotte, D. L., Holliday, J. R., and Rundle, J. B. (2007), <em>Geophys. Res. Lett.</em>, 34, L12303, doi:10.1029/2007GL029696.</p> </li> <li> <p>Zaliapin, I., and Ben‐Zion, Y. (2013), <em>J. Geophys. Res. Solid Earth</em>, 118, 2865– 2877, doi:10.1002/jgrb.50178.</p> </li> </ul>


2015 ◽  
Vol 57 (6) ◽  
Author(s):  
Maura Murru ◽  
Jiancang Zhuang ◽  
Rodolfo Console ◽  
Giuseppe Falcone

<div class="page" title="Page 1"><div class="layoutArea"><div class="column"><p>In this paper, we compare the forecasting performance of several statistical models, which are used to describe the occurrence process of earthquakes in forecasting the short-term earthquake probabilities during the L’Aquila earthquake sequence in central Italy in 2009. These models include the Proximity to Past Earthquakes (PPE) model and two versions of the Epidemic Type Aftershock Sequence (ETAS) model. We used the information gains corresponding to the Poisson and binomial scores to evaluate the performance of these models. It is shown that both ETAS models work better than the PPE model. However, in comparing the two types of ETAS models, the one with the same fixed exponent coefficient (<span>alpha)</span> = 2.3 for both the productivity function and the scaling factor in the spatial response function (ETAS I), performs better in forecasting the active aftershock sequence than the model with different exponent coefficients (ETAS II), when the Poisson score is adopted. ETAS II performs better when a lower magnitude threshold of 2.0 and the binomial score are used. The reason is found to be that the catalog does not have an event of similar magnitude to the L’Aquila mainshock (M<sub>w</sub> 6.3) in the training period (April 16, 2005 to March 15, 2009), and the (<span>alpha)</span>-value is underestimated, thus the forecast seismicity is underestimated when the productivity function is extrapolated to high magnitudes. We also investigate the effect of the inclusion of small events in forecasting larger events. These results suggest that the training catalog used for estimating the model parameters should include earthquakes of magnitudes similar to the mainshock when forecasting seismicity during an aftershock sequence.</p></div></div></div>


Solid Earth ◽  
2019 ◽  
Vol 10 (5) ◽  
pp. 1519-1540
Author(s):  
Marisol Monterrubio-Velasco ◽  
F. Ramón Zúñiga ◽  
José Carlos Carrasco-Jiménez ◽  
Víctor Márquez-Ramírez ◽  
Josep de la Puente

Abstract. Earthquake aftershocks display spatiotemporal correlations arising from their self-organized critical behavior. Dynamic deterministic modeling of aftershock series is challenging to carry out due to both the physical complexity and uncertainties related to the different parameters which govern the system. Nevertheless, numerical simulations with the help of stochastic models such as the fiber bundle model (FBM) allow the use of an analog of the physical model that produces a statistical behavior with many similarities to real series. FBMs are simple discrete element models that can be characterized by using few parameters. In this work, the aim is to present a new model based on FBM that includes geometrical characteristics of fault systems. In our model, the faults are not described with typical geometric measures such as dip, strike, and slip, but they are incorporated as weak regions in the model domain that could increase the likelihood to generate earthquakes. In order to analyze the sensitivity of the model to input parameters, a parametric study is carried out. Our analysis focuses on aftershock statistics in space, time, and magnitude domains. Moreover, we analyzed the synthetic aftershock sequences properties assuming initial load configurations and suitable conditions to propagate the rupture. As an example case, we have modeled a set of real active faults related to the Northridge, California, earthquake sequence. We compare the simulation results to statistical characteristics from the Northridge sequence determining which range of parameters in our FBM version reproduces the main features observed in real aftershock series. From the results obtained, we observe that two parameters related to the initial load configuration are determinant in obtaining realistic seismicity characteristics: (1) parameter P, which represents the initial probability order, and (2) parameter π, which is the percentage of load distributed to the neighboring cells. The results show that in order to reproduce statistical characteristics of the real sequence, larger πfrac values (0.85<πfrac<0.95) and very low values of P (0.0<P≤0.08) are needed. This implies the important corollary that a very small departure from an initial random load configuration (computed by P), and also a large difference between the load transfer from on-fault segments than by off-faults (computed by πfrac), is required to initiate a rupture sequence which conforms to observed statistical properties such as the Gutenberg–Richter law, Omori law, and fractal dimension.


Sign in / Sign up

Export Citation Format

Share Document