skip to main content

Title: Invariant Galton–Watson branching process for earthquake occurrence

We propose a theoretical modelling framework for earthquake occurrence and clustering based on a family of invariant Galton–Watson (IGW) stochastic branching processes. The IGW process is a rigorously defined approximation to imprecisely observed or incorrectly estimated earthquake clusters modelled by Galton–Watson branching processes, including the Epidemic Type Aftershock Sequence (ETAS) model. The theory of IGW processes yields explicit distributions for multiple cluster attributes, including magnitude-dependent and magnitude-independent offspring number, cluster size and cluster combinatorial depth. Analysis of the observed seismicity in southern California demonstrates that the IGW model provides a close fit to the observed earthquake clusters. The estimated IGW parameters and derived statistics are robust with respect to the catalogue lower cut-off magnitude. The proposed model facilitates analyses of multiple quantities of seismicity based on self-similar tree attributes, and may be used to assess the proximity of seismicity to criticality.

more » « less
Award ID(s):
1723033 2122191 1722561 2122168
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Geophysical Journal International
Medium: X Size: p. 567-583
["p. 567-583"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Precipitation clusters are contiguous raining regions characterized by a precipitation threshold, size, and the total rainfall contained within—termed the cluster power. Tropical observations suggest that the probability distributions of both cluster size and power contain a power-law range (with slope ~ −1.5) bounded by a large-event “cutoff.” Events with values beyond the cutoff signify large, powerful clusters and represent extreme events. A two-dimensional stochastic model is introduced to reproduce the observed cluster distributions, including the slope and the cutoff. The model is equipped with coupled moisture and weak temperature gradient (WTG) energy equations, empirically motivated precipitation parameterization, temporally persistent noise, and lateral mixing processes, all of which collectively shape the model cluster distributions. Moisture–radiative feedbacks aid clustering, but excessively strong feedbacks push the model into a self-aggregating regime. The power-law slope is stable in a realistic parameter range. The cutoff is sensitive to multiple model parameters including the stochastic forcing amplitude, the threshold moisture value that triggers precipitation, and the lateral mixing efficiency. Among the candidates for simple analogs of precipitation clustering, percolation models are ruled out as unsatisfactory, but the stochastic branching process proves useful in formulating a neighbor probability metric. This metric measures the average number of nearest neighbors that a precipitating entity can spawn per time interval and captures the cutoff parameter sensitivity for both cluster size and power. The results here suggest that the clustering tendency and the horizontal scale limiting large tropical precipitating systems arise from aggregate effects of multiple moist processes, which are encapsulated in the neighbor probability metric.

    more » « less
  2. Abstract

    Foreshocks can provide valuable information about possible nucleation process of a mainshock. However, their physical mechanisms are still under debate. In this study, we present a comprehensive analysis of the earthquake sequence preceding the 2010 Mw7.2 El Mayor‐Cucapah mainshock, including waveform detection of missing smaller events, relative relocation, and source parameter analysis. Based on a template matching method, we find a tenfold increase in the number of earthquakes than reported in the Southern California Seismic Network catalog. The entire sequence exhibits nearly continuous episodes of foreshocks that can be loosely separated into two active clusters. Relocated foreshocks show several seismicity streaks at depth, with a consistently active cluster at depths between 14 and 16 km where the mainshock was nucleated. Stress drop measurements from a spectral ratio approach based on empirical Green's functions show a range between 3.8 and 41.7 MPa with a median of 13.0 MPa and no clear temporal variations. The relocation results, together with the source patches estimated from earthquake corner frequencies, revealed a migration front toward the mainshock hypocenter within last 8 hr and a chain of active burst immediately 6 min prior to the mainshock. Our results support combined effects of aseismic slip and cascading failure on the evolution of foreshocks.

    more » « less
  3. null (Ed.)
    SUMMARY We examine localization processes of low magnitude seismicity in relation to the occurrence of large earthquakes using three complementary analyses: (i) estimated production of rock damage by background events, (ii) evolving occupied fractional area of background seismicity and (iii) progressive coalescence of individual earthquakes into clusters. The different techniques provide information on different time scales and on the spatial extent of weakened damaged regions. Techniques (i) and (ii) use declustered catalogues to avoid the occasional strong fluctuations associated with aftershock sequences, while technique (iii) examines developing clusters in entire catalogue data. We analyse primarily earthquakes around large faults that are locked in the interseismic periods, and examine also as a contrasting example seismicity from the creeping Parkfield section of the San Andreas fault. Results of analysis (i) show that the M > 7 Landers 1992, Hector Mine 1999, El Mayor-Cucapah 2010 and Ridgecrest 2019 main shocks in Southern and Baja California were preceded in the previous decades by generation of rock damage around the eventual rupture zones. Analysis (ii) reveals localization (reduced fractional area) 2–3 yr before these main shocks and before the M > 7 Düzce 1999 earthquake in Turkey. Results with technique (iii) indicate that individual events tend to coalesce rapidly to clusters in the final 1–2 yr before the main shocks. Corresponding analyses of data from the Parkfield region show opposite delocalization patterns and decreasing clustering before the 2004 M6 earthquake. Continuing studies with these techniques, combined with analysis of geodetic data and insights from laboratory experiments and model simulations, might improve the ability to track preparation processes leading to large earthquakes. 
    more » « less

    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
  5. Abstract

    We introduce an algorithm for declustering earthquake catalogs based on the nearest‐neighbor analysis of seismicity. The algorithm discriminates between background and clustered events by random thinning that removes events according to a space‐varying threshold. The threshold is estimated using randomized‐reshuffled catalogs that are stationary, have independent space and time components, and preserve the space distribution of the original catalog. Analysis of catalog produced by the Epidemic Type Aftershock Sequence model demonstrates that the algorithm correctly classifies over 80% of background and clustered events, correctly reconstructs the stationary and space‐dependent background intensity, and shows high stability with respect to random realizations (over 75% of events have the same estimated type in over 90% of random realizations). The declustering algorithm is applied to the global Northern California Earthquake Data Center catalog with magnitudesm≥ 4 during 2000–2015; a Southern California catalog withm≥ 2.5, 3.5 during 1981–2017; an area around the 1992 Landers rupture zone withm≥ 0.0 during 1981–2015; and the Parkfield segment of San Andreas fault withm≥ 1.0 during 1984–2014. The null hypotheses of stationarity and space‐time independence are not rejected by several tests applied to the estimated background events of the global and Southern California catalogs with magnitude ranges Δm< 4. However, both hypotheses are rejected for catalogs with larger range of magnitudes Δm> 4. The deviations from the nulls are mainly due to local temporal fluctuations of seismicity and activity switching among subregions; they can be traced back to the original catalogs and represent genuine features of background seismicity.

    more » « less