Including Covariates in the ETAS Model Triggered Seismicity

2020 ◽  
Author(s):  
marcello chiodi ◽  
Giada Adelfio
2021 ◽  
Vol 11 (19) ◽  
pp. 9143
Author(s):  
Marcello Chiodi ◽  
Orietta Nicolis ◽  
Giada Adelfio ◽  
Nicoletta D’Angelo ◽  
Alex Gonzàlez

Chilean seismic activity is one of the strongest in the world. As already shown in previous papers, seismic activity can be usefully described by a space–time branching process, such as the ETAS (Epidemic Type Aftershock Sequences) model, which is a semiparametric model with a large time-scale component for the background seismicity and a small time-scale component for the triggered seismicity. The use of covariates can improve the description of triggered seismicity in the ETAS model, so in this paper, we study the Chilean seismicity separately for the North and South area, using some GPS-related data observed together with ordinary catalog data. Our results show evidence that the use of some covariates can improve the fitting of the ETAS model.


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>


2021 ◽  
Author(s):  
Grzegorz Lizurek ◽  
Konstantinos Leptokaropoulos ◽  
Jan Wiszniowski ◽  
Izabela Nowaczyńska ◽  
Nguyen Van Giang ◽  
...  

<p>Reservoir-triggered seismicity (RTS) is the longest known anthropogenic seismicity type. It has the potential to generate seismic events of M6 and bigger. Previous studies of this phenomenon have proved that major events are triggered on preexisting major discontinuities, forced to slip by stress changes induced by water level fluctuations and/or pore-pressure changes in the rock mass in the vicinity of reservoirs. Song Tranh 2 is an artificial water reservoir located in Central Vietnam. Its main goal is back up the water for hydropower plant. High seismic activity has been observed in this area since the reservoir was first filled in 2011. The relation between water level and seismic activity in the Song Tranh area is complex, and the lack of clear correlation between water level and seismic activity has led to the conclusion that ongoing STR2 seismic activity is an example of the delayed response type of RTS. However, the first phase of the activity observed after impoundment has been deemed a rapid response type. In this work, we proved that the seismicity recorded between 2013 and 2016 manifested seasonal trends related to water level changes during wet and dry seasons. The response of activity and its delay with respect to water level changes suggest that the main triggering factor is pore pressure change due to the significant water level changes observed. A stress orientation difference between low and high water periods is also revealed. The findings indicate that water load and related pore pressure changes influence seismic activity and stress orientation in this area.</p><p>This work was partially supported by research project no. 2017/27/B/ST10/01267, funded by the National Science Centre, Poland, under agreement no. UMO-2017/27/B/ST10/01267.</p>


2021 ◽  
Author(s):  
Kwang-Il Kim ◽  
Hwajung Yoo ◽  
Seheok Park ◽  
Juhyi Yim ◽  
Linmao Xie ◽  
...  

<p>Hydraulic stimulation for the creation of an Enhanced Geothermal System (EGS) reservoir could potentially reactivate a nearby fault and result in man-made earthquakes. In November 15, 2017, an M<sub>w</sub> 5.5 earthquake, the second largest after the initiation of the South Korean national instrumental monitoring system, occurred near an EGS project in Pohang, South Korea. The earthquake occurred on a previously unmapped fault, that is here denoted the M<sub>w</sub> 5.5 Fault. A number of previous studies to model the hydraulic stimulation in the Pohang EGS project have been carried out to identify the mechanism of seismic events. Those previous studies focused on coupled hydro-mechanical processes without the consideration of pre-existing fractures and thermal effects. This study presents an investigation of the mechanisms of induced and triggered seismicity in the Pohang EGS project through three-dimensional coupled thermo-hydro-mechanical numerical simulations. Fractures intersecting the open-hole sections of two deep boreholes, PX-1 and PX-2, clearly indicated by field observations are modeled along with the M<sub>w</sub> 5.5 Fault. Models of stress-dependent permeability models are calibrated based on the numerical reproduction of the pressure-time evolution during the field hydraulic stimulations. The Coulomb failure stress change at the M<sub>w</sub> 5.5 Fault is calculated to quantify the impact of five hydraulic stimulations. In the case of PX-2 stimulations, the pore pressure buildup results in a volumetric expansion of the reservoir and thereby the perturbation of stresses is transferred to the M<sub>w</sub> 5.5 Fault. The volumetric contraction of the reservoir by the temperature reduction could slightly perturb the stress distribution at the M<sub>w</sub> 5.5 Fault. In the case of PX-1 stimulations, shear slip of the PX-1 fracture is explicitly modeled. The modeling shows that transfer of the shear stress drop by the shear slip stabilizes the M<sub>w</sub> 5.5 Fault, which is consistent with the field observation that the seismicity was not induced at the M<sub>w</sub> 5.5 Fault by the PX-1 stimulations. The cooling-induced thermal stress additionally reduces the effective normal stress of PX-1 fracture. Thus, some additional shear slip of the PX-1 fracture is induced by the thermal effect. However, the modeling shows that for both PX-1 and PX-2 stimulations, thermally-induced stress perturbations are very small compared to pressure-induced stress perturbations.</p>


2012 ◽  
Vol 39 (10) ◽  
pp. n/a-n/a ◽  
Author(s):  
Hector Gonzalez-Huizar ◽  
Aaron A. Velasco ◽  
Zhigang Peng ◽  
Raul R. Castro
Keyword(s):  

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.


Sign in / Sign up

Export Citation Format

Share Document