Seismicity Declustering and Hazard Analysis of the Oklahoma–Kansas Region

2019 ◽  
Vol 109 (6) ◽  
pp. 2356-2366 ◽  
Author(s):  
Ganyu Teng ◽  
Jack W. Baker
Keyword(s):  
Induced Seismicity ◽  
Hazard Analysis ◽  
Aftershock Sequence ◽  
Eastern United States ◽  
Etas Model ◽  
Epidemic Type Aftershock Sequence ◽  
Window Method ◽  
Seismicity Declustering ◽  
Hazard Level ◽  
Seismicity Catalog

Abstract This study is an evaluation of the suitability of several declustering method for induced seismicity and their impacts on hazard analysis of the Oklahoma–Kansas region. We considered the methods proposed by Gardner and Knopoff (1974), Reasenberg (1985), Zaliapin and Ben‐Zion (2013), and the stochastic declustering method (Zhuang et al., 2002) based on the epidemic‐type aftershock sequence (ETAS) model (Ogata, 1988, 1998). The results show that the choice of declustering method has a significant impact on the declustered catalog and the resulting hazard analysis of the Oklahoma–Kansas region. The Gardner and Knopoff method, which is currently implemented in the U.S. Geological Survey one‐year seismic‐hazard forecast for the central and eastern United States, has unexpected features when used for this induced seismicity catalog. It removes 80% of earthquakes and fails to reflect the changes in background rates that have occurred in the past few years. This results in a slight increase in the hazard level from 2016 to 2017, despite a decrease in seismic activities in 2017. The Gardner and Knopoff method also frequently identifies aftershocks with much stronger shaking intensities than their associated mainshocks. These features are mostly due to the window method implemented in the Gardner and Knopoff method. Compared with the Gardner and Knopoff method, the other three methods are able to capture the changing hazard level in the region. However, the ETAS model potentially overestimates the foreshock effect and generates negligible probabilities of large earthquakes being mainshocks. The Reasenberg and Zaliapin and Ben‐Zion methods have similar performance on catalog declustering and hazard analysis. Compared with the ETAS method, these two methods are easier to implement and faster to generate the declustered catalog. The results from this study suggest that both Reasenberg and Zaliapin and Ben‐Zion declustering methods are suitable for declustering and hazard analysis for induced seismicity in the Oklahoma–Kansas region.

Improving Earthquake Doublet Frequency Predictions by Modified Spatial Trigger Kernels in the Epidemic-Type Aftershock Sequence (ETAS) Model

2021 ◽  
Author(s):  
Christian Grimm ◽  
Martin Käser ◽  
Sebastian Hainzl ◽  
Marco Pagani ◽  
Helmut Küchenhoff
Keyword(s):  
Hazard Analysis ◽  
Aftershock Sequence ◽  
Economic Losses ◽  
Etas Model ◽  
Independent Events ◽  
Epidemic Type Aftershock Sequence ◽  
Time Space ◽  
Spatiotemporal Clustering ◽  
Background Events ◽  
The Impact

ABSTRACT Earthquake sequences add a substantial hazard beyond the solely declustered perspective of common probabilistic seismic hazard analysis. A particularly strong driver for both social and economic losses are so-called earthquake doublets (more generally multiplets), that is, sequences of two (or more) comparatively large events in spatial and temporal proximity. Without differentiating between foreshocks and aftershocks, we hypothesize three main influencing factors of doublet occurrence: (1) the number of direct and secondary aftershocks triggered by an earthquake; (2) the occurrence of independent clusters and seismic background events in the same time–space window; and (3) the magnitude size distribution of triggered events (in contrast to independent events). We tested synthetic catalogs simulated by a standard epidemic-type aftershock sequence (ETAS) model for both Japan and southern California. Our findings show that the common ETAS approach significantly underestimates doublet frequencies compared with observations in historical catalogs. In combination with that the simulated catalogs show a smoother spatiotemporal clustering compared with the observed counterparts. Focusing on the impact on direct aftershock productivity and total cluster sizes, we propose two modifications of the ETAS spatial kernel to improve doublet rate predictions: (a) a restriction of the spatial function to a maximum distance of 2.5 estimated rupture lengths and (b) an anisotropic function with contour lines constructed by a box with two semicircular ends around the estimated rupture segment. These modifications shift the triggering potential from weaker to stronger events and consequently improve doublet rate predictions for larger events, despite still underestimating historic doublet occurrence rates. Besides, the results for the restricted spatial functions fulfill better the empirical Båth’s law for the maximum aftershock magnitude. The tested clustering properties of strong events are not sufficiently incorporated in typically used global catalog scale measures, such as log-likelihood values, which would favor the conventional, unrestricted models.

scholarly journals Probabilistic Forecasting of Hydraulic Fracturing-Induced Seismicity Using an Injection-Rate Driven ETAS Model

2021 ◽  
Author(s):  
Simone Mancini ◽  
Maximilian Jonas Werner ◽  
Margarita Segou ◽  
Brian Baptie
Keyword(s):  
Hydraulic Fracturing ◽  
Shale Gas ◽  
Induced Seismicity ◽  
Risk Mitigation ◽  
Injection Rate ◽  
Aftershock Sequence ◽  
Probabilistic Forecasting ◽  
Etas Model ◽  
Epidemic Type Aftershock Sequence ◽  
Out Of Sample

Abstract The development of robust forecasts of human-induced seismicity is highly desirable to mitigate the effects of disturbing or damaging earthquakes. We assess the performance of a well-established statistical model, the epidemic-type aftershock sequence (ETAS) model, with a catalog of ∼93,000 microearthquakes observed at the Preston New Road (PNR, United Kingdom) unconventional shale gas site during, and after hydraulic fracturing of the PNR-1z and PNR-2 wells. Because ETAS was developed for slower loading rate tectonic seismicity, to account for seismicity caused by pressurized fluid, we also generate three modified ETAS with background rates proportional to injection rates. We find that (1) the standard ETAS captures low seismicity between and after injections but is outperformed by the modified model during high-seismicity periods, and (2) the injection-rate driven ETAS substantially improves when the forecast is calibrated on sleeve-specific pumping data. We finally forecast out-of-sample the PNR-2 seismicity using the average response to injection observed at PNR-1z, achieving better predictive skills than the in-sample standard ETAS. The insights from this study contribute toward producing informative seismicity forecasts for real-time decision making and risk mitigation techniques during unconventional shale gas development.

Operational Earthquake Forecasting during the 2019 Ridgecrest, California, Earthquake Sequence with the UCERF3-ETAS Model

2020 ◽  
Vol 91 (3) ◽  
pp. 1567-1578 ◽  
Author(s):  
Kevin R. Milner ◽  
Edward H. Field ◽  
William H. Savran ◽  
Morgan T. Page ◽  
Thomas H. Jordan
Keyword(s):  
High Performance ◽  
Aftershock Sequence ◽  
Lessons Learned ◽  
Ease Of Use ◽  
Earthquake Forecasting ◽  
Earthquake Rupture ◽  
Etas Model ◽  
Epidemic Type Aftershock Sequence ◽  
Aftershock Sequences ◽  
High Performance Computing Cluster

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.

Testing of the foreshock hypothesis within an epidemic like description of seismicity

2020 ◽  
Author(s):  
G Petrillo ◽  
E Lippiello
Keyword(s):  
Seismic Activity ◽  
Good Description ◽  
Aftershock Sequence ◽  
Short Term ◽  
Etas Model ◽  
Temporal Clustering ◽  
Epidemic Type Aftershock Sequence ◽  
Preparatory Phase ◽  
Spatio Temporal ◽  
The One

Summary The Epidemic Type Aftershock Sequence (ETAS) model provides a good description of the post-seismic spatio-temporal clustering of seismicity and is also able to capture some features of the increase of seismic activity caused by foreshocks. Recent results, however, have shown that the number of foreshocks observed in instrumental catalogs is significantly much larger than the one predicted by the ETAS model. Here we show that it is possible to keep an epidemic description of post-seismic activity and, at the same time, to incorporate pre-seismic temporal clustering, related to foreshocks. Taking also into-account the short-term incompleteness of instrumental catalogs, we present a model which achieves very good description of the southern California seismicity both on the aftershock and on the foreshock side. Our results indicate that the existence of a preparatory phase anticipating mainshocks represents the most plausible explanation for the occurrence of foreshocks.

Point process modelling of reservoir-induced seismicity

2001 ◽  
Vol 38 (A) ◽  
pp. 232-242 ◽  
Author(s):  
Masajiro Imoto
Keyword(s):  
Point Process ◽  
Induced Seismicity ◽  
Earthquake Hazard ◽  
Information Criterion ◽  
Aftershock Sequence ◽  
Water Levels ◽  
Tectonic Activity ◽  
Reservoir Induced Seismicity ◽  
Epidemic Type Aftershock Sequence ◽  
Background Seismicity

A point process procedure can be used to study reservoir-induced seismicity (RIS), in which the intensity function representing earthquake hazard is a combination of three terms: a constant background term, an ETAS (epidemic-type aftershock sequence) term for aftershocks, and a time function derived from observation of water levels of a reservoir. This paper presents the results of such a study of the seismicity in the vicinity of the Tarbela reservoir in Pakistan. Making allowance for changes in detection capability and the background seismicity related to tectonic activity, earthquakes of magnitude ≥ 2.0, occurring between May 1978 and January 1982 and whose epicentres were within 100 km of the reservoir, were used in this analysis. Several different intensities were compared via their Akaike information criterion (AIC) values relative to those of a Poisson process. The results demonstrate that the seismicity within 20 km of the reservoir correlates with water levels of the reservoir, namely, active periods occur about 250 days after the appearance of low water levels. This suggests that unloading the reservoir activates the seismicity beneath it. Seasonal variations of the seismicity in an area up to 100 km from the reservoir were also found, but these could not be adequately interpreted by an appropriate RIS mechanism.

Improving earthquake doublet rate predictions in ETAS by using modified spatial trigger distributions

2021 ◽  
Author(s):  
Christian Grimm ◽  
Martin Käser ◽  
Sebastian Hainzl ◽  
Marco Pagani ◽  
Helmut Küchenhoff
Keyword(s):  
Hazard Analysis ◽  
Aftershock Sequence ◽  
Model Fit ◽  
Occurrence Rate ◽  
Economic Losses ◽  
Anisotropic Models ◽  
Independent Events ◽  
Epidemic Type Aftershock Sequence ◽  
Time Space ◽  
The Impact

<p>Earthquake sequences add significant hazard beyond the solely declustered perspective of common probabilistic seismic hazard analysis (PSHA). A particularly strong driver for both social and economic losses are so-called earthquake doublets (more generally multiplets), i.e. sequences of two (or more) comparatively large events in spatial and temporal proximity. Not differentiating between foreshocks and aftershocks, we hypothesize three main drivers of doublet occurrence: (1) the number of direct aftershocks triggered by an earthquake; (2) the underlying, independent background seismicity in the same time-space window; and (3) the magnitude size distribution of triggered events (in contrast to independent events). We tested synthetic catalogs simulated by a common, isotropic epidemic type aftershock sequence (ETAS) model for both Japan and Southern California. Our findings show that the standard ETAS approach dramatically underestimates doublet frequencies compared to observations in historical catalogs. Among others, the results partially smooth out pronounced peaks of temporal and spatial event clustering. Focusing on the impact on direct aftershock productivity, we propose two modifications of the ETAS spatial kernel in order to improve doublet rate predictions: (a) a restriction of the spatial function to a maximum distance of 2.5 estimated rupture lengths; (b) an anisotropic function with contour lines constructed by a box with two semicircular ends around the estimated rupture line. The restriction of the spatial extent shifts triggering potential from weaker to stronger events and in consequence improves doublet rate predictions for larger events. However, this improvement goes at the cost of a weaker overall model fit according to AIC. The anisotropic models improve the overall model fit, but have minor impact on doublet occurrence rate predictions.</p>

How Do Statistical Parameters of Induced Seismicity Correlate with Fluid Injection? Case of Oklahoma

2021 ◽  
Author(s):  
Hideo Aochi ◽  
Julie Maury ◽  
Thomas Le Guenan
Keyword(s):  
Induced Seismicity ◽  
Aftershock Sequence ◽  
Statistical Parameters ◽  
Exponential Equation ◽  
Seismicity Rate ◽  
Epidemic Type Aftershock Sequence ◽  
Background Seismicity ◽  
Mechanical Interpretation ◽  
Independent Variable ◽  
Grid Points

Abstract The seismicity evolution in Oklahoma between 2010 and 2018 is analyzed systematically using an epidemic-type aftershock sequence model. To retrieve the nonstationary seismicity component, we systematically use a moving window of 200 events, each within a radius of 20 km at grid points spaced every 0.2°. Fifty-three areas in total are selected for our analysis. The evolution of the background seismicity rate μ is successfully retrieved toward its peak at the end of 2014 and during 2015, whereas the triggering parameter K is stable, slightly decreasing when the seismicity is activated. Consequently, the ratio of μ to the observed seismicity rate is not stationary. The acceleration of μ can be fit with an exponential equation relating μ to the normalized injected volume. After the peak, the attenuation phase can be fit with an exponential equation with time since peak as the independent variable. As a result, the evolution of induced seismicity can be followed statistically after it begins. The turning points, such as activation of the seismicity and timing of the peak, are difficult to identify solely from this statistical analysis and require a subsequent mechanical interpretation.

scholarly journals Weighted complex networks applied to seismicity: the Itoiz dam (Northern Spain)

2011 ◽  
Vol 18 (4) ◽  
pp. 477-487 ◽  
Author(s):  
A. Jiménez ◽  
F. Luzón
Keyword(s):  
Complex Networks ◽  
Induced Seismicity ◽  
Aftershock Sequence ◽  
Northern Spain ◽  
Topological Features ◽  
Epidemic Type Aftershock Sequence ◽  
Main Shocks ◽  
Random Models ◽  
Seismic Clusters ◽  
Weighted Complex Networks

Abstract. On 18 September 2004, an earthquake of magnitude mbLg = 4.6 was recorded near the Itoiz dam (Northern Spain). It occurred after the first impoundment of the reservoir and has been catalogued by some authors as induced seismicity. We analyzed the seismicity in the region as weighted complex networks and tried to differentiate this event from others that occurred nearby. We calculated the main topological features of the networks formed by the seismic clusters and compared them. We compared the results with a series of simulations, and showed that the clusters were better modelled with the Epidemic-Type Aftershock Sequence (ETAS) model than with random models. We found that the properties of the different clusters are grouped according to the magnitude of the main shocks and the number of events in each cluster, and that no distinct feature could be obtained for the 18 September 2004 series. We found that the nodes with the highest strength are the most important in the networks' traffic, and are associated with the events with the highest magnitude within the clusters.

scholarly journals The Relevance of Foreshocks in Earthquake Triggering: A Statistical Study

2018 ◽  
Author(s):  
Eugenio Lippiello ◽  
Cataldo Godano ◽  
Lucilla De Arcangelis
Keyword(s):  
Statistical Study ◽  
Aftershock Sequence ◽  
Statistical Features ◽  
Earthquake Triggering ◽  
Etas Model ◽  
Epidemic Type Aftershock Sequence ◽  
The One ◽  
Triggering Process ◽  
Large Earthquakes ◽  
Good Agreement

An increase of seismic activity is often observed before large earthquakes. Events responsible for this increase are usually named foreshock and their occurrence probably represents the most reliable precursory pattern. Many foreshocks statistical features can be interpreted in terms of the standard mainshock-to-aftershock triggering process and are recovered in the Epidemic Type Aftershock Sequence ETAS model. Here we present a statistical study of instrumental seismic catalogs from four different geographic regions. We focus on some common features of foreshocks in the four catalogs which cannot be reproduced by the ETAS model. In particular we find in instrumental catalogs a significantly larger number of foreshocks than the one predicted by the ETAS model. We show that this foreshock excess cannot be attributed to catalog incompleteness. We therefore propose a generalized formulation of the ETAS model, the ETAFS model, which explicitly includes foreshock occurrence. Statistical features of aftershocks and foreshocks in the ETAFS model are in very good agreement with instrumental results.

Exploring probabilistic seismic risk assessment accounting for seismicity clustering and damage accumulation: Part I. Hazard analysis

2020 ◽  
pp. 875529302095733
Author(s):  
Athanasios N Papadopoulos ◽  
Paolo Bazzurro ◽  
Warner Marzocchi
Keyword(s):  
Seismic Risk ◽  
Hazard Analysis ◽  
Central Italy ◽  
Building Design ◽  
Companion Paper ◽  
Aftershock Sequence ◽  
Risk Models ◽  
Epidemic Type Aftershock Sequence ◽  
One Year ◽  
Compare And Contrast

Probabilistic seismic hazard analysis (PSHA), as a tool to assess the probability that ground motion of a given intensity or larger is experienced at a given site and time span, has historically comprised the basis of both building design codes in earthquake-prone regions and seismic risk models. The PSHA traditionally refers solely to mainshock events and typically employs a homogeneous Poisson process to model their occurrence. Nevertheless, recent disasters, such as the 2010–2011 Christchurch sequence or the 2016 Central Italy earthquakes, to name a few, have highlighted the potential pitfalls of neglecting the occurrence of foreshocks, aftershocks, and other triggered events, and pinpointed the need to revisit the current practice. Herein, we employ the epidemic-type aftershock sequence (ETAS) model to describe seismicity in Central Italy, investigate the model’s capability to reproduce salient features of observed seismicity, and compare ETAS-derived one-year hazard estimates with ones obtained with a standard mainshock-only Poisson-based hazard model. A companion paper uses the hazard models derived herein to compare and contrast loss estimates for the residential exposure of Umbria in Central Italy.

Sign in / Sign up

Export Citation Format

Share Document