skip to main content

Title: Cascading foreshocks, aftershocks and earthquake swarms in a discrete fault network

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
Author(s) / Creator(s):
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Geophysical Journal International
Page Range / eLocation ID:
p. 831-852
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Foreshocks are the only currently widely identified precursory seismic behavior, yet their utility and even identifiability are problematic, in part because of extreme variation in behavior. Here, we establish some global trends that help identify the expected frequency of foreshocks as well the type of earthquake most prone to foreshocks. We establish these tendencies using the global earthquake catalog of the U.S. Geological Survey National Earthquake Information Center with a completeness level of magnitude 5 and mainshocks with Mw≥7.0. Foreshocks are identified using three clustering algorithms to address the challenge of distinguishing foreshocks from background activity. The methods give a range of 15%–43% of large mainshocks having at least one foreshock but a narrower range of 13%–26% having at least one foreshock with magnitude within two units of the mainshock magnitude. These observed global foreshock rates are similar to regional values for a completeness level of magnitude 3 using the same detection conditions. The foreshock sequences have distinctive characteristics with the global composite population b-values being lower for foreshocks than for aftershocks, an attribute that is also manifested in synthetic catalogs computed by epidemic-type aftershock sequences, which intrinsically involves only cascading processes. Focal mechanism similarity of foreshocks relative to mainshocks is more pronounced than for aftershocks. Despite these distinguishing characteristics of foreshock sequences, the conditions that promote high foreshock productivity are similar to those that promote high aftershock productivity. For instance, a modestly higher percentage of interplate mainshocks have foreshocks than intraplate mainshocks, and reverse faulting events slightly more commonly have foreshocks than normal or strike-slip-faulting mainshocks. The western circum-Pacific is prone to having slightly more foreshock activity than the eastern circum-Pacific.

    more » « less
  2. Abstract

    We conduct a detailed study of the foreshock sequence preceding the 2010Mw6.7 Yushu, Qinghai earthquake in the Tibetan plateau by examining continuous waveforms recorded at a seismic station near the mainshock rupture zone. By using a deep learning phase picker—EQTransformer and a matched‐filter technique, we identify 120 foreshocks with magnitude ranging from −0.7 to 1.6, starting with aMw4.6 foreshock approximately 2 hr before theMw6.7 Yushu mainshock. Our analyses show that the foreshock sequence follows a typical Omori's law decay with ap‐value of 0.73 and the Gutenberg‐Richer frequency‐magnitudeb‐value of 0.66. We do not find any evidence of accelerating events leading up to the Yushu mainshock. Hence, they could be considered as aftershocks of theMw4.6 earthquake. We further invert for the focal mechanisms and rupture directions for both the largest foreshock and the mainshock. TheMw4.6 foreshock likely occurred on a NE‐SW trending fault conjugating to the NW‐SE trending fault of the mainshock. Coulomb stress analysis suggests theMw4.6 foreshock induces negative stress on the mainshock source area. These observations do not support either the pre‐slip or the cascade triggering model for foreshock generation. The occurrence of the foreshock, mainshock and large aftershocks appear to be modulated by the Earth's tidal forces, likely reflecting the role of high pore‐fluid pressures. Our observations, together with other recent studies, suggest that extensional step‐overs and conjugate faults along major strike‐slip faults play an important role in generating short‐term foreshock sequences.

    more » « less
  3. Abstract

    The 2021MW6.0 Yangbi, Yunnan strike‐slip earthquake occurred on an unmapped crustal fault near the Weixi‐Qiaoho‐Weishan Fault along the southeast margin of the Tibetan Plateau. Using near‐source broadband seismic data from ChinArray, we investigate the spatial and temporal rupture evolution of the mainshock using apparent moment‐rate functions (AMRFs) determined by the empirical Green's function (EGF) method. Assuming a 1D line source on the fault plane, the rupture propagated unilaterally southeastward (∼144°) over a rupture length of ∼8.0 km with an estimated rupture speed of 2.1 km/s to 2.4 km/s. A 2D coseismic slip distribution for an assumed maximum rupture propagation speed of 2.2 km/s indicates that the rupture propagated to the southeast ∼8.0 km along strike and ∼5.0 km downdip with a peak slip of ∼2.1 m before stopping near the largest foreshock, where three bifurcating subfaults intersect. Using the AMRFs, the radiated energy of the mainshock is estimated as ∼. The relatively low moment scaled radiated energyof 1.5 × 10−5and intense foreshock and aftershock activity might indicate reactivation of an immature fault. The earthquake sequence is mainly distributed along a northwest‐southeast trend, and aftershocks and foreshocks are distributed near the periphery of the mainshock large‐slip area, suggesting that the stress in the mainshock slip zone is significantly reduced to below the level for more than a few overlapping aftershock to occur.

    more » « less
  4. The recent 2019 Ridgecrest earthquake sequence in Southern California jostled the seismological community by revealing a complex and cascading foreshock series that culminated in a M7.1 mainshock. But the central Garlock fault, despite being located immediately south of this sequence, did not coseismically fail. Instead, the Garlock fault underwent post-seismic creep and exhibited a sizeable earthquake swarm. The dynamic details of the rupture process during the mainshock is largely unknown, as is the amount of stress needed to bring the Garlock fault to failure. We present an integrated view of how stresses changed on the Garlock fault during and after the mainshock using a combination of tools including kinematic slip inversion, Coulomb stress change, and dynamic rupture modeling. We show that positive Coulomb stress changes cannot easily explain observed aftershock patterns on the Garlock fault, but are consistent with where creep was documented on the central Garlock fault section. Our dynamic model is able to reproduce the main slip asperities and kinematically estimated rupture speeds (≤ 2 km/s) during the mainshock, and suggests the temporal changes in normal and shear stress on the Garlock fault were greatest near the end of rupture. The largest static and dynamic stress changes on the Garlock fault we observe from our models coincide with the creeping region, suggesting that positive stress perturbations could have caused this during or after the mainshock rupture. This analysis of near-field stress change evolution gives insight into how the Ridgecrest sequence influenced the local stress field of the northernmost Eastern California Shear Zone. 
    more » « less
  5. Oceanic transform faults display a unique combination of seismic and aseismic slip behavior, including a large globally averaged seismic deficit, and the local occurrence of repeating magnitude (M)6earthquakes with abundant foreshocks and seismic swarms, as on the Gofar transform of the East Pacific Rise and the Blanco Ridge in the northeast Pacific Ocean. However, the underlying mechanisms that govern the partitioning between seismic and aseismic slip and their interaction remain unclear. Here we present a numerical modeling study of earthquake sequences and aseismic transient slip on oceanic transform faults. In the model, strong dilatancy strengthening, supported by seismic imaging that indicates enhanced fluid-filled porosity and possible hydrothermal circulation down to the brittle–ductile transition, effectively stabilizes along-strike seismic rupture propagation and results in rupture barriers where aseismic transients arise episodically. The modeled slow slip migrates along the barrier zones at speeds ∼10 to 600 m/h, spatiotemporally correlated with the observed migration of seismic swarms on the Gofar transform. Our model thus suggests the possible prevalence of episodic aseismic transients in M6rupture barrier zones that host active swarms on oceanic transform faults and provides candidates for future seafloor geodesy experiments to verify the relation between aseismic fault slip, earthquake swarms, and fault zone hydromechanical properties.

    more » « less