Anomalous Power Law Distribution of Total Lifetimes of Aftershock Sequences: Application to Earthquake Aftershock Sequences

Alexander I. Saichev, & Didier Sornette

Published October 2004, SCEC Contribution #858

We consider a general stochastic branching process, which is relevant to earthquakes, and study the distributions of global lifetimes of the branching processes. In the earthquake context, this amounts to the distribution of the total durations of aftershock sequences including aftershocks of arbitrary generation number. Our results extend previous results on the distribution of the total number of offspring (direct and indirect aftershocks in seismicity) and of the total number of generations before extinction. We consider a branching model of triggered seismicity, the epidemic-type aftershock sequence model, which 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 (“productivity” or “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 μ is characterized by a power law ∼1∕μ1+γ and the bare Omori law for the memory of previous triggering mothers decays slowly as ∼1∕t1+θ , with 0<θ<1 relevant for earthquakes. Using the tool of generating probability functions and a quasistatic approximation which is shown to be exact asymptotically for large durations, we show that the density distribution of total aftershock lifetimes scales as ∼1∕t1+θ∕γ when the average branching ratio is critical (n=1) . The coefficient 1<γ=b∕α<2 quantifies the interplay between the exponent b≈1 of the Gutenberg-Richter magnitude distribution ∼10−bm and the increase ∼10αm of the number of aftershocks with mainshock magnitude m (productivity), with 0.5<α<1 . The renormalization of the bare Omori decay law ∼1∕t1+θ into ∼1∕t1+θ∕γ stems from the nonlinear amplification due to the heavy-tailed distribution of fertilities and the critical nature of the branching cascade process. In the subcritical case n<1 , the crossover from ∼1∕t1+θ∕γ at early times to ∼1∕t1+θ at longer times is described. More generally, our results apply to any stochastic branching process with a power-law distribution of offspring per parent and a long memory.

Key Words
foreshocks, models

Citation
Saichev, A. I., & Sornette, D. (2004). Anomalous Power Law Distribution of Total Lifetimes of Aftershock Sequences: Application to Earthquake Aftershock Sequences. Physical Review E, 70(4), 046123. doi: 10.1103/PhysRevE.70.046123.