Power Law Distributions of Offspring and Generation Numbers in Branching Models of Earthquake Triggering

A. Saichev,A. Helmstetter,D. Sornette

DOI: https://doi.org/10.1007/s00024-004-2663-6

2004-03-20

Abstract:We consider a general stochastic branching process, which is relevant to earthquakes as well as to many other systems, and we study the distributions of the total number of offsprings (direct and indirect aftershocks in seismicity) and of the total number of generations before extinction. We apply our results to a branching model of triggered seismicity, the ETAS (epidemic-type aftershock sequence) model. The ETAS model assumes that each earthquake can trigger other earthquakes (``aftershocks''). An aftershock sequence results in this model from the cascade of aftershocks of each past earthquake. Due to the large fluctuations of the number of aftershocks triggered directly by any earthquake (``fertility''), there is a large variability of the total number of aftershocks from one sequence to another, for the same mainshock magnitude. We study the regime where the distribution of fertilities mu is characterized by a power law ~1/\mu^(1+gamma). For earthquakes, we expect such a power-law distribution of fertilities with gamma = b/alpha based on the Gutenberg-Richter magnitude distribution ~10^(-bm) and on the increase ~10^(alpha m) of the number of aftershocks with the mainshock magnitude m. We derive the asymptotic distributions p_r(r) and p_g(g) of the total number r of offsprings and of the total number g of generations until extinction following a mainshock. In the regime \gamma<2 relevant for earhquakes, for which the distribution of fertilities has an infinite variance, we find p_r(r)~1/r^(1+1/gamma) and p_g(g)~1/g^(1+1/(gamma -1)). These predictions are checked by numerical simulations.

Statistical Mechanics

What problem does this paper attempt to address?

### Problems the paper attempts to solve

This paper aims to study the problems of the number of descendants and the distribution of generations in the branching process in the earthquake triggering model. Specifically, the author focuses on the distribution of all direct and indirect aftershocks (i.e., the total number of descendants of all generations) triggered by the main shock in the branching process and the total number of generations before extinction.

#### Research background and motivation
1. **Earthquake triggering model**:
- Earthquakes are usually followed by a series of aftershocks, and these aftershocks can further trigger more earthquakes, forming a cascading process.
- The ETAS (Epidemic - Type Aftershock Sequence) model assumes that each earthquake can trigger other earthquakes (called "aftershocks"), and the aftershock sequence is caused by the aftershock cascades of each past earthquake.

2. **Statistical laws**:
- The magnitude distribution of earthquakes follows the Gutenberg - Richter law: \[ p_m(m)\sim10^{-bm} \]
- There is a relationship between the magnitude \( m \) of the main shock and the number \( N_m \) of its directly triggered aftershocks: \[ N_m = K10^{\alpha m} \]

3. **Core of the problem**:
- Since the number of aftershocks triggered by each earthquake (i.e., "fertility") has great volatility, the total number of aftershocks after different main shocks varies greatly.
- The author pays special attention to the case where the distribution of fertility \( \mu \) is a power - law distribution, that is, \( p_\mu(\mu)\sim\frac{1}{\mu^{1 + \gamma}} \), where \( \gamma=\frac{b}{\alpha} \).

#### Research objectives
- **Derive the distribution**: The author's goal is to provide a general and exact derivation to determine the distribution of the total number \( r \) of aftershocks of all generations and the total number \( g \) of generations after the main shock.
- **Anomalous scaling under the power - law distribution**: When the distribution of fertility has an infinite variance (i.e., \( \gamma < 2 \)), the author finds that the distributions of the total number of descendants and the total number of generations present new power - law forms:
\[
p_r(r)\sim\frac{1}{r^{1+\frac{1}{\gamma}}}, \quad p_g(g)\sim\frac{1}{g^{1+\frac{1}{\gamma - 1}}}
\]
This is different from the distribution in the case of finite variance in the standard branching process.

#### Application scope
- **Seismology**: These results can be directly applied to the ETAS model, especially in the case of \( \alpha = 0.8 \) and \( b = 1 \), where the distribution of fertility has an infinite variance.
- **Wider applications**: These results are also applicable to any stochastic branching process with a power - law distributed fertility, such as chain reactions, material fractures, economic and social cascades and other phenomena.

Through the above analysis, the problem that this paper attempts to solve is to understand and quantify the statistical characteristics of the number of all aftershocks triggered by the main shock and their generation distribution in the earthquake triggering model, especially in the case where the fertility distribution has power - law characteristics.