<p>Earthquake catalogs exhibit strong spatio-temporal correlations. As such, earthquakes are often classified into clusters of correlated activity. Clusters themselves are traditionally classified in two different kinds: (i) bursts, with a clear hierarchical structure between a single strong mainshock, preceded by a few foreshocks and followed by a power-law decaying aftershock sequence, and (ii) swarms, exhibiting a non-trivial activity rate that cannot be reduced to such a simple hierarchy between events.&#160;</p><p>The Epidemic Aftershock Sequence (ETAS) model is a linear Hawkes point process able to reproduce earthquake clusters from empirical statistical laws [Ogata, 1998]. Although not always explicit, the ETAS model is often interpreted as the outcome of a background activity driven by external forces and a Galton-Watson branching process with one-to-one causal links between events [Saichev et al., 2005]. Declustering techniques based on field observations [Baiesi & Paczuski, 2004] can be used to infer the most likely causal links between events in a cluster. Following this method, Zaliapin and Ben&#8208;Zion (2013) determined the statistical properties of earthquake clusters characterizing bursts and swarms, finding a relationship between the predominant cluster-class and the heat flow in seismic regions.</p><p>Here, I show how the statistical properties of clusters are related to the fundamental statistics of the underlying seismogenic process, modeled in two point-process paradigms [Bar&#243;, 2020].</p><p>The classification of clusters into bursts and swarms appears naturally in the standard ETAS model with homogeneous rates and are determined by the average branching ratio (nb) and the ratio between exponents &#945; and b characterizing the production of aftershocks and the distribution of magnitudes, respectively. The scale-free ETAS model, equivalent to the BASS model [Turcotte, et al., 2007], and usual in cold active tectonic regions, is imposed by &#945;=b and reproduces bursts. In contrast, by imposing &#945;<0.5b, we recover the properties of swarms, characteristic of regions with high heat flow.&#160;</p><p>Alternatively, the same declustering methodology applied to a non-homogeneous Poisson process with a non-factorizable intensity, i.e. in absence of causal links, recovers swarms with &#945;=0, i.e. a Poisson Galton-Watson process, with similar statistical properties to the ETAS model in the regime &#945;<0.5b.</p><p>Therefore, while bursts are likely to represent actual causal links between events, swarms can either denote causal links with low &#945;/b ratio or variations of the background rate caused by exogenous processes introducing local and transient stress changes. Furthermore, the redundancy in the statistical laws can be used to test the hypotheses posed by the ETAS model as a memory&#8208;less branching process.&#160;</p><p>References:</p><ul><li>
<p>Baiesi, M., & Paczuski, M. (2004). <em>Physical Review E</em>, 69, 66,106. doi:10.1103/PhysRevE.69.066106.</p>
</li>
<li>
<p>Bar&#243;, J. (2020).&#160; <em>Journal of Geophysical Research: Solid Earth,</em> 125, e2019JB018530. doi:10.1029/2019JB018530.</p>
</li>
<li>
<p>Ogata, Y. (1998) <em>Annals of the Institute of Statistical Mathematics,</em> 50(2), 379&#8211;402. doi:10.1023/A:1003403601725.</p>
</li>
<li>
<p>Saichev, A., Helmstetter, A. & Sornette, D. (2005) <em>Pure appl. geophys.</em> 162, 1113&#8211;1134. doi:10.1007/s00024-004-2663-6.</p>
</li>
<li>
<p>Turcotte, D. L., Holliday, J. R., and Rundle, J. B. (2007), <em>Geophys. Res. Lett.</em>, 34, L12303, doi:10.1029/2007GL029696.</p>
</li>
<li>
<p>Zaliapin, I., and Ben&#8208;Zion, Y. (2013), <em>J. Geophys. Res. Solid Earth</em>, 118, 2865&#8211; 2877, doi:10.1002/jgrb.50178.</p>
</li>
</ul>
Spherical K-Means method to determine earthquake clusters
