Detection and characterization of fluid-driven earthquake clusters

Author(s):  
Sebastian Hainzl ◽  
Tomas Fischer

<p>Natural earthquake clusters are often related to a mainshock, which triggers the sequence by its induced stress changes. These clusters are called mainshock-aftershock sequences and statistically well explained by earthquake-earthquake interactions according to the Epidemic Type Aftershock Sequence (ETAS) model. Additionally, aseismic processes such as slow slip, dike propagation or fluid flow might also play a role in the initiation and driving of the earthquake sequence. Earthquake swarms, which lacks a dominant earthquake, are often believed to indicate such transient aseismic forcing signals. However, swarm-type clusters can also occur by chance in ETAS-simulations and thus not necessarily related to aseismic drivers. Thus, more sophisticated quantification of the space-time-magnitude characteristics of earthquake sequences are required for discrimination. Migration patterns are one of those properties which can be indicative for aseismic triggering. We suggest simple measures to identify and quantify migration patterns and test those for synthetic data, data from fluid injection experiments, and natural swarm activity related to fluid flow in NW Bohemia and Long Valley caldera. We analyze their potential to discriminate from ETAS-type clusters and compare it with those of time-magnitude characteristics of the activity such as seismic moment ratios and skewness. Our results are finally used to discriminate earthquake clusters in California and elsewhere.</p>

2021 ◽  
Author(s):  
Shubham Sharma ◽  
Shyam Nandan ◽  
Sebastian Hainzl

<p>Currently, the Epidemic Type Aftershock Sequence (ETAS) model is state-of-the-art for forecasting aftershocks. However, the under-performance of ETAS in forecasting the spatial distribution of aftershocks following a large earthquake make us adopt alternative approaches for the modelling of the spatial ETAS-kernel. Here we develop a hybrid physics and statics based forecasting model. The model uses stress changes, calculated from inverted slip models of large earthquakes, as the basis of the spatial kernel in the ETAS model in order to get more reliable estimates of spatiotemporal distribution of aftershocks. We evaluate six alternative approaches of stress-based ETAS-kernels and rank their performance against the base ETAS model. In all cases, an expectation maximization (EM) algorithm is used to estimate the ETAS parameters. The model approach has been tested on synthetic data to check if the known parameters can be inverted successfully. We apply the proposed method to forecast aftershocks of mainshocks available in SRCMOD database, which includes 192 mainshocks with magnitudes in the range between 4.1 and 9.2 occurred from 1906 to 2020. The probabilistic earthquake forecasts generated by the hybrid model have been tested using established CSEP test metrics and procedures. We show that the additional stress information, provided to estimate the spatial probability distribution, leads to more reliable spatiotemporal ETAS-forecasts of aftershocks as compared to the base ETAS model.</p>


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>


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.


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.


2017 ◽  
Vol 43 (4) ◽  
pp. 1994
Author(s):  
A.C. Astiopoulos ◽  
E. Papadimitriou ◽  
V. Karakostas ◽  
D. Gospodinov ◽  
G. Drakatos

The statistical properties of the aftershock occurrence are among the main issues in investigating the earthquake generation process. Seismicity rate changes during a seismic sequence, which are detected by the application of statistical models, are proved to be precursors of strong events occurring during the seismic excitation. Application of these models provides a tool in assessing the imminent seismic hazard, oftentimes by the estimation of the expected occurrence rate and comparison of the predicted rate with the observed one. The aim of this study is to examine the temporal distribution and especially the occurrence rate variations of aftershocks for two seismic sequences that took place, the first one near Skyros island in 2001 and the second one near Lefkada island in 2003, in order to detect and determine rate changes in connection with the evolution of the seismic activity. Analysis is performed through space–time stochastic models which are developed, based upon both aftershocks clustering studies and specific assumptions. The models applied are the Modified Omori Formula (MOF), the Epidemic Type Aftershock Sequence (ETAS) and the Restricted Epidemic Type Aftershock Sequence (RETAS). The modelling of seismicity rate changes, during the evolution of the particular seismic sequences, is then attempted in association with and as evidence of static stress changes


2020 ◽  
Author(s):  
Robert Churchill ◽  
Maximilian Werner ◽  
Juliet Biggs ◽  
Ake Fagereng

<p>Aftershock sequences following large tectonic earthquakes exhibit considerable spatio-temporal complexity and suggest causative mechanisms beyond co-seismic, elasto-static Coulomb stress changes in the crust. Candidate mechanisms include dynamic triggering and postseismic processes such as viscoelastic relaxation, poroelastic rebound and aseismic afterslip, which has garnered particular interest recently. Aseismic afterslip – whereby localized frictional sliding within velocity-strengthening rheologies acts to redistribute lithospheric stresses in the postseismic phase – has been suggested by numerous studies to exert dominant control on aftershock sequence evolution, including productivity, spatial distribution and temporal decay.</p><p>As evidence is based overwhelmingly on individual case study analysis, we wish to systematically compare key metrics of aseismic afterslip and corresponding aftershock sequences to investigate this relationship. We specifically look for any empirical relationship between the seismic-equivalent moment of aseismic afterslip episodes and the corresponding aftershock sequence productivity. We first compile published afterslip models into a database containing moment estimates over varying time periods, as well as spatial distributions, temporal decays and modelling methodology as a supplementary resource. We then identify the corresponding aftershock sequence from the globally comparable USGS PDE catalog. As expected, coseismic moment exerts an obvious control on both afterslip moment and aftershock productivity – an effect we control for by normalising by mainshock moment and expected productivity (the Utsu-Seki law) respectively. Preliminary results suggest broad variability of both afterslip moment and aftershock productivity with no obvious control of afterslip on aftershocks beyond the scaling with mainshock size, including when separated by mainshock mechanism or region. As this study is insensitive to spatial and temporal distributions, we cannot rule out the potential influence afterslip exerts in these but find no evidence that afterslip drives overall productivity of aftershock sequences.</p>


2017 ◽  
Vol 50 (3) ◽  
pp. 1283
Author(s):  
K.A. Adamaki ◽  
R.G. Roberts

We investigate temporal changes in seismic activity observed in the West Corinth Gulf and North-West Peloponnese during 2008 to 2010. Two major earthquake sequences took place in the area at that time (in 2008 and 2010). Our aim is to analyse Greek seismicity to attempt to confirm the existence or non-existence of seismic precursors prior to the strongest earthquakes. Perhaps because the area is geologically and tectonically complex, we found that it was not possible to fit the data well using a consistent Epidemic Type Aftershock Sequence (ETAS) model. Nor could we unambiguously identify foreshocks to individual mainshocks. Therefore we sought patterns in aggregated foreshock catalogues. We set a magnitude threshold (M3.5) above which all the earthquakes detected in the study area are considered as “mainshocks”, and we combined all data preceding these into a single foreshock catalogue. This reveals an increase in seismicity rate not robustly observable for individual cases. The observed effect is significantly greater than that consistent with stochastic models, including ETAS, thus indicating genuine foreshock activity with potential useful precursory power, if sufficient data is available, i.e. if the magnitude of completeness is sufficiently low.


2020 ◽  
Author(s):  
José Ángel López-Comino ◽  
Martin Galis ◽  
P. Martin Mai ◽  
Xiaowei Chen ◽  
Daniel Stich

<p>Exploring the connections between injection wells and seismic migration patterns is key to understanding processes controlling growth of fluid-injection induced seismicity. Numerous seismic clusters in Oklahoma have been associated with wastewater disposal operations, providing a unique opportunity to investigate migration directions of each cluster with respect to the injection-well locations. We introduce new directivity migration parameters to identify and quantify lateral migration toward or away from the injection wells. We take into account cumulative volume and injection rate from multiple injection wells. Our results suggest a weak relationship between migration direction and the cluster-well distances. Migration away from injection wells is found for distances shorter than 5-13 km, while an opposite migration towards the wells is observed for larger distances, suggesting an increasing influence of poroelastic stress changes. This finding is more stable when considering cumulative injected volume instead of injection rate. We do not observe any relationship between migration direction and injected volume or equivalent magnitudes.</p>


2001 ◽  
Vol 38 (A) ◽  
pp. 213-221 ◽  
Author(s):  
Li Ma ◽  
Jiancang Zhuang

This paper uses the epidemic-type aftershock sequence (ETAS) point process model to study certain seismicity features of the Jiashi swarm of certain earthquakes, investigating in particular whether there is relative quiescence prior to the quite big events within the Jiashi sequence. The seven earthquake sequences studied occurred in the region of Jiashi, south of Tianshan Mountain, Xinjiang, China. The particular ETAS model that is developed is consistent with the reality of seismic activity. The various features of Jiashi swarm activity can be described as focusing in different stages. There is obvious precursory quiescence prior to most big events with Ms ≥ 6.0 within the Jiashi swarm. Thus, checking for relative quiescence can be use for earthquake prediction.


2001 ◽  
Vol 38 (A) ◽  
pp. 213-221 ◽  
Author(s):  
Li Ma ◽  
Jiancang Zhuang

This paper uses the epidemic-type aftershock sequence (ETAS) point process model to study certain seismicity features of the Jiashi swarm of certain earthquakes, investigating in particular whether there is relative quiescence prior to the quite big events within the Jiashi sequence. The seven earthquake sequences studied occurred in the region of Jiashi, south of Tianshan Mountain, Xinjiang, China. The particular ETAS model that is developed is consistent with the reality of seismic activity. The various features of Jiashi swarm activity can be described as focusing in different stages. There is obvious precursory quiescence prior to most big events with Ms ≥ 6.0 within the Jiashi swarm. Thus, checking for relative quiescence can be use for earthquake prediction.


Sign in / Sign up

Export Citation Format

Share Document