Sandpile-based model for capturing magnitude distributions and spatiotemporal clustering and separation in regional earthquakes

Abstract. We propose a cellular automata model for earthquake occurrences patterned after the sandpile model of self-organized criticality (SOC). By incorporating a single parameter describing the probability to target the most susceptible site, the model successfully reproduces the statistical signatures of seismicity. The energy distributions closely follow power-law probability density functions (PDFs) with a scaling exponent of around −1. 6, consistent with the expectations of the Gutenberg–Richter (GR) law, for a wide range of the targeted triggering probability values. Additionally, for targeted triggering probabilities within the range 0.004–0.007, we observe spatiotemporal distributions that show bimodal behavior, which is not observed previously for the original sandpile. For this critical range of values for the probability, model statistics show remarkable comparison with long-period empirical data from earthquakes from different seismogenic regions. The proposed model has key advantages, the foremost of which is the fact that it simultaneously captures the energy, space, and time statistics of earthquakes by just introducing a single parameter, while introducing minimal parameters in the simple rules of the sandpile. We believe that the critical targeting probability parameterizes the memory that is inherently present in earthquake-generating regions.

[1]  Keisuke Ito,et al.  Earthquakes as self-organized critical phenomena , 1990 .

[2]  Christensen,et al.  Self-organized criticality in a continuous, nonconservative cellular automaton modeling earthquakes. , 1992, Physical review letters.

[3]  Tang,et al.  Self-Organized Criticality: An Explanation of 1/f Noise , 2011 .

[4]  Stefan Boettcher,et al.  Interoccurrence times in the Bak-Tang-Wiesenfeld sandpile model: a comparison with the observed statistics of solar flares. , 2005, Physical review letters.

[5]  R. Batac Statistical Properties of the Immediate Aftershocks of the 15 October 2013 Magnitude 7.1 Earthquake in Bohol, Philippines , 2016, Acta Geophysica.

[6]  P. Bak,et al.  Earthquakes as a self‐organized critical phenomenon , 1989 .

[7]  L. Knopoff,et al.  Model and theoretical seismicity , 1967 .

[8]  Ilya Zaliapin,et al.  Clustering analysis of seismicity and aftershock identification. , 2007, Physical review letters.

[9]  Bruce D. Malamud,et al.  Cellular-automata models applied to natural hazards , 2000, Comput. Sci. Eng..

[10]  S. Lübeck Large-scale simulations of the Zhang sandpile model , 1997 .

[12]  H. Kantz,et al.  Observing spatio-temporal clustering and separation using interevent distributions of regional earthquakes , 2014 .

[13]  Christopher Monterola,et al.  Loss of criticality in the avalanche statistics of sandpiles with dissipative sites , 2015, Commun. Nonlinear Sci. Numer. Simul..

[14]  John McCloskey,et al.  Heterogeneity in a self-organized critical earthquake model , 1996 .

[15]  D. E. Juanico,et al.  Avalanche Statistics of Driven Granular Slides in a Miniature Mound , 2008 .

[16]  Zhang,et al.  Scaling theory of self-organized criticality. , 1989, Physical review letters.

[17]  F. Landes,et al.  Scaling laws in earthquake occurrence: Disorder, viscosity, and finite size effects in Olami-Feder-Christensen models. , 2016, Physical Review E.

[18]  M. Baiesi,et al.  Scale-free networks of earthquakes and aftershocks. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[19]  I. Zaliapin,et al.  Artefacts of earthquake location errors and short-term incompleteness on seismicity clusters in southern California , 2015 .

[20]  E A Jagla,et al.  Forest-fire analogy to explain the b value of the Gutenberg-Richter law for earthquakes. , 2013, Physical review letters.

[21]  D Sornette,et al.  "Universal" distribution of interearthquake times explained. , 2006, Physical review letters.

[22]  Drossel,et al.  Self-organized critical forest-fire model. , 1992, Physical review letters.

[23]  A scale‐invariant cellular‐automata model for distribited seismicity , 1991 .

[24]  T. V. McEvilly Seismicity of the earth and associated phenomena , 1967 .

[25]  Shlomo Havlin,et al.  Memory in the occurrence of earthquakes. , 2005, Physical review letters.

[26]  Mark Naylor,et al.  Origin and nonuniversality of the earthquake interevent time distribution. , 2009, Physical review letters.

[27]  Stefan Hergarten,et al.  Foreshocks and aftershocks in the Olami-Feder-Christensen model. , 2002, Physical review letters.