An official website of the United States government
SUMMARY Earthquakes come in clusters formed of mostly aftershock sequences, swarms and occasional foreshock sequences. This clustering is thought to result either from stress transfer among faults, a process referred to as cascading, or from transient loading by aseismic slip (pre-slip, afterslip or slow slip events). The ETAS statistical model is often used to quantify the fraction of clustering due to stress transfer and to assess the eventual need for aseismic slip to explain foreshocks or swarms. Another popular model of clustering relies on the earthquake nucleation model derived from experimental rate-and-state friction. According to this model, earthquakes cluster because they are time-advanced by the stress change imparted by the mainshock. This model ignores stress interactions among aftershocks and cannot explain foreshocks or swarms in the absence of transient loading. Here, we analyse foreshock, swarm and aftershock sequences resulting from cascades in a Discrete Fault Network model governed by rate-and-state friction. We show that the model produces realistic swarms, foreshocks and aftershocks. The Omori law, characterizing the temporal decay of aftershocks, emerges in all simulations independently of the assumed initial condition. In our simulations, the Omori law results from the earthquake nucleation process due to rate and state friction and from the heterogeneous stress changes due to the coseismic stress transfers. By contrast, the inverse Omori law, which characterizes the accelerating rate of foreshocks, emerges only in the simulations with a dense enough fault system. A high-density complex fault zone favours fault interactions and the emergence of an accelerating sequence of foreshocks. Seismicity catalogues generated with our discrete fault network model can generally be fitted with the ETAS model but with some material differences. In the discrete fault network simulations, fault interactions are weaker in aftershock sequences because they occur in a broader zone of lower fault density and because of the depletion of critically stressed faults. The productivity of the cascading process is, therefore, significantly higher in foreshocks than in aftershocks if fault zone complexity is high. This effect is not captured by the ETAS model of fault interactions. It follows that a foreshock acceleration stronger than expected from ETAS statistics does not necessarily require aseismic slip preceding the mainshock (pre-slip). It can be a manifestation of a cascading process enhanced by the topological properties of the fault network. Similarly, earthquake swarms might not always imply transient loading by aseismic slip, as they can emerge from stress interactions.
more »
« less
A deep Gaussian process model for seismicity background rates
SUMMARY The spatio-temporal properties of seismicity give us incisive insight into the stress state evolution and fault structures of the crust. Empirical models based on self-exciting point processes continue to provide an important tool for analysing seismicity, given the epistemic uncertainty associated with physical models. In particular, the epidemic-type aftershock sequence (ETAS) model acts as a reference model for studying seismicity catalogues. The traditional ETAS model uses simple parametric definitions for the background rate of triggering-independent seismicity. This reduces the effectiveness of the basic ETAS model in modelling the temporally complex seismicity patterns seen in seismic swarms that are dominated by aseismic tectonic processes such as fluid injection rather than aftershock triggering. In order to robustly capture time-varying seismicity rates, we introduce a deep Gaussian process (GP) formulation for the background rate as an extension to ETAS. GPs are a robust non-parametric model for function spaces with covariance structure. By conditioning the length-scale structure of a GP with another GP, we have a deep-GP: a probabilistic, hierarchical model that automatically tunes its structure to match data constraints. We show how the deep-GP-ETAS model can be efficiently sampled by making use of a Metropolis-within-Gibbs scheme, taking advantage of the branching process formulation of ETAS and a stochastic partial differential equation (SPDE) approximation for Matérn GPs. We illustrate our method using synthetic examples, and show that the deep-GP-ETAS model successfully captures multiscale temporal behaviour in the background forcing rate of seismicity. We then apply the results to two real-data catalogues: the Ridgecrest, CA 2019 July 5 Mw 7.1 event catalogue, showing that deep-GP-ETAS can successfully characterize a classical aftershock sequence; and the 2016–2019 Cahuilla, CA earthquake swarm, which shows two distinct phases of aseismic forcing concordant with a fluid injection-driven initial sequence, arrest of the fluid along a physical barrier and release following the largest Mw 4.4 event of the sequence.
more »
« less
- Award ID(s):
- 2034167
- PAR ID:
- 10400984
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 234
- Issue:
- 1
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- p. 427-438
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The development of new earthquake forecasting models is often motivated by one of the following complementary goals: to gain new insights into the governing physics and to produce improved forecasts quantified by objective metrics. Often, one comes at the cost of the other. Here, we propose a question-driven ensemble (QDE) modeling approach to address both goals. We first describe flexible epidemic-type aftershock sequence (ETAS) models in which we relax the assumptions of parametrically defined aftershock productivity and background earthquake rates during model calibration. Instead, both productivity and background rates are calibrated with data such that their variability is optimally represented by the model. Then we consider 64 QDE models in pseudoprospective forecasting experiments for southern California and Italy. QDE models are constructed by combining model parameters of different ingredient models, in which the rules for how to combine parameters are defined by questions about the future seismicity. The QDE models can be interpreted as models that address different questions with different ingredient models. We find that certain models best address the same issues in both regions, and that QDE models can substantially outperform the standard ETAS and all ingredient models. The best performing QDE model is obtained through the combination of models allowing flexible background seismicity and flexible aftershock productivity, respectively, in which the former parameterizes the spatial distribution of background earthquakes and the partitioning of seismicity into background events and aftershocks, and the latter is used to parameterize the spatiotemporal occurrence of aftershocks.more » « less
-
Foreshocks are the most obvious signature of the earthquake nucleation stage and could, in principle, forewarn of an impending earthquake. However, foreshocks are only sometimes observed, and we have a limited understanding of the physics that controls their occurrence. In this work, we use high-resolution earthquake catalogs and estimates of source properties to understand the spatiotemporal evolution of a sequence of 11 foreshocks that occurred ~ 6.5 hours before the 2020 Mw 4.8 Mentone earthquake in west Texas. Elevated pore-pressure and poroelastic stressing from subsurface fluid injection from oil-gas operations is often invoked to explain seismicity in west Texas and the surrounding region. However, here we show that static stresses induced from the initial ML 4.0 foreshock significantly perturbed the local shear stress along the fault and could have triggered the Mentone mainshock. The majority (9/11) of the earthquakes leading up to the Mentone mainshock nucleated in areas where the static shear stresses were increased from the initial ML 4.0 foreshock. The spatiotemporal properties of the 11 earthquakes that preceded the mainshock cannot easily be explained in the context of a preslip or cascade nucleation model. We show that at least 6/11 events are better classified as aftershocks of the initial ML 4.0. Together, our results suggest that a combination of physical mechanisms contributed to the occurrence of the 11 earthquakes that preceded the mainshock, including static-stressing from earthquake-earthquake interactions, aseismic creep, and stress perturbations induced from fluid injection. Our work highlights the role of earthquake-earthquake triggering in induced earthquake sequences, and suggests that such triggering could help sustain seismic activity following initial stressing perturbations from fluid injection.more » « less
-
Abstract In areas of induced seismicity, earthquakes can be triggered by stress changes due to fluid injection and static deformation from fault slip. Here we present a method to distinguish between injection‐driven and earthquake‐driven triggering of induced seismicity by combining a calibrated, fully coupled, poroelastic stress model of wastewater injection with interpretation of a machine learning algorithm trained on both earthquake catalog and modeled stress features. We investigate seismicity from Paradox Valley, Colorado as an ideal test case: a single, high‐pressure injector that has induced thousands of earthquakes since 1991. Using feature importance analysis, we find that injection‐driven earthquakes are approximately 225% of the total catalog but act as background events that can trigger subsequent aftershocks. Injection‐driven events also have distinct spatiotemporal clustering properties with a larger b‐value, closer proximity to the well, and earlier occurrence in the injection history. Generalization of our technique can help characterize triggering processes in other regions where induced seismicity occurs.more » « less
-
Abstract Microseismicity associated with fluid pressurization in the subsurface occurs during fluid injection but can also be triggered after injection shut‐in. Understanding the extent and duration of the post‐injection microseismicity is critical to limit the risk of fluid‐induced seismicity and insure the safe utilization of the subsurface. Using theoretical and numerical techniques, we investigated how aseismic slip on a fault plane evolves and stops after a fluid pressurization event. We found that the locking mechanisms controlling the arrest of aseismic slip highly depend on the initial fault stress criticality and the pressurization duration. The absolute arrest time of fault aseismic slip after injection shut‐in is proportional to the pressurization duration and increases significantly with the initial fault stress criticality. Given that microseismicity can be triggered by aseismic slip, these results provide insights into the mechanics controlling the arrest of microseismicity after fluid pressurization as a milestone toward induced seismicity mitigation strategies.more » « less
