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.
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
- NSF-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
-
more » « less
SUMMARY -
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
-
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.
-
Earthquake swarms attributed to subsurface fluid injection are usually assumed to occur on faults destabilized by increased pore-fluid pressures. However, fluid injection could also activate aseismic slip, which might outpace pore-fluid migration and transmit earthquake-triggering stress changes beyond the fluid-pressurized region. We tested this theoretical prediction against data derived from fluid-injection experiments that activated and measured slow, aseismic slip on preexisting, shallow faults. We found that the pore pressure and slip history imply a fault whose strength is the product of a slip-weakening friction coefficient and the local effective normal stress. Using a coupled shear-rupture model, we derived constraints on the hydromechanical parameters of the actively deforming fault. The inferred aseismic rupture front propagates faster and to larger distances than the diffusion of pressurized pore fluid.more » « less
-
Abstract The recent Indios, Puerto Rico earthquake sequence has drawn attention, as the increased seismicity rate in this area was unprecedented. The sequence began on 28 December 2019, caused a 6.4 magnitude earthquake on 7 January 2020, and remained active over a year later. This sequence fits the nominal definition of an earthquake swarm in that it had an abrupt onset, a sustained high rate of seismicity without a clear triggering mainshock or evidence for Omori decay, and a lack of adherence to Bath’s law. However, the sequence also had several prominent mainshock–aftershock (MS–AS) sequences embedded within it. We applied three-station waveform cross correlation to the early part of this sequence using the Puerto Rico Seismic Network (PRSN) catalog as templates, which confirmed the mixture of swarm and MS–AS patterns. In an effort to place this intriguing sequence in the context of the previous seismicity in Puerto Rico, we investigated the existence of swarms and MS–AS sequences recorded by the PRSN since 1987 by identifying sequences with increased seismicity rate when compared to the background rate. About 59 sequences were manually verified and characterized into swarms or MS–AS. We found that 58% of the sequences follow traditional swarm patterns and 14% adhere to traditional MS–AS behavior, whereas 29% of the sequences have a mixture of both swarm and MS–AS behaviors. These findings suggest that it is not unusual for the Indios sequence to have a mixture of both the characteristics. In addition, the detection of many swarms distributed over a broad area of the subduction interface indicates stress heterogeneity and low-coupling consistent with prior studies indicating that the potential for a magnitude ∼8 megathrust earthquake along the Puerto Rico trench is unlikely.more » « less
