No Arabic abstract
The aftershock productivity law, first described by Utsu in 1970, is an exponential function of the form K=K0.exp({alpha}M) where K is the number of aftershocks, M the mainshock magnitude, and {alpha} the productivity parameter. The Utsu law remains empirical in nature although it has also been retrieved in static stress simulations. Here, we explain this law based on Solid Seismicity, a geometrical theory of seismicity where seismicity patterns are described by mathematical expressions obtained from geometric operations on a permanent static stress field. We recover the exponential form but with a break in scaling predicted between small and large magnitudes M, with {alpha}=1.5ln(10) and ln(10), respectively, in agreement with results from previous static stress simulations. We suggest that the lack of break in scaling observed in seismicity catalogues (with {alpha}=ln(10)) could be an artefact from existing aftershock selection methods, which assume a continuous behavior over the full magnitude range. While the possibility for such an artefact is verified in simulations, the existence of the theoretical kink remains to be proven.
Declustering aims to divide earthquake catalogs into independent events (mainshocks), and dependent (clustered) events, and is an integral component of many seismicity studies, including seismic hazard assessment. We assess the effect of declustering on the frequency-magnitude distribution of mainshocks. In particular, we examine the dependence of the b-value of declustered catalogs on the choice of declustering approach and algorithm-specific parameters. Using the catalog of earthquakes in California since 1980, we show that the b-value decreases by up to 30% due to declustering with respect to the undeclustered catalog. The extent of the reduction is highly dependent on the declustering method and parameters applied. We then reproduce a similar effect by declustering synthetic earthquake catalogs with known b-value, which have been generated using an Epidemic-Type Aftershock Sequence (ETAS) model. Our analysis suggests that the observed decrease in b-value must, at least partially, arise from the application of the declustering algorithm on the catalog, rather than from differences in the nature of mainshocks versus fore- or aftershocks. We conclude that declustering should be considered as a potential source of bias in seismicity and hazard studies.
A likely source of earthquake clustering is static stress transfer between individual events. Previous attempts to quantify the role of static stress for earthquake triggering generally considered only the stress changes caused by large events, and often discarded data uncertainties. We conducted a robust two-fold empirical test of the static stress change hypothesis by accounting for all events of magnitude M>=2.5 and their location and focal mechanism uncertainties provided by catalogs for Southern California between 1981 and 2010, first after resolving the focal plane ambiguity and second after randomly choosing one of the two nodal planes. For both cases, we find compelling evidence supporting the static triggering with stronger evidence after resolving the focal plane ambiguity above significantly small (about 10 Pa) but consistently observed stress thresholds. The evidence for the static triggering hypothesis is robust with respect to the choice of the friction coefficient, Skemptons coefficient and magnitude threshold. Weak correlations between the Coulomb Index (fraction of earthquakes that received positive Coulomb stress change) and the coefficient of friction indicate that the role of normal stress in triggering is rather limited. Last but not the least, we determined that the characteristic time for the loss of the stress change memory of a single event is nearly independent of the amplitude of the Coulomb stress change and varies between ~95 and ~180 days implying that forecasts based on static stress changes will have poor predictive skills beyond times that are larger than a few hundred days on average.
Seeing the Earth crust as crisscrossed by faults filled with fluid at close to lithostatic pressures, we develop a model in which its elastic modulii are different in net tension versus compression. In constrast with standard nonlinear effects, this ``threshold nonlinearity is non-perturbative and occurs for infinitesimal perturbations around the lithostatic pressure taken as the reference. For a given earthquake source, such nonlinear elasticity is shown to (i) rotate, widen or narrow the different lobes of stress transfer, (ii) to modify the $1/r^2$ 2D-decay of elastic stress Green functions into the generalized power law $1/r^{gamma}$ where $gamma$ depends on the azimuth and on the amplitude of the modulii asymmetry. Using reasonable estimates, this implies an enhancement of the range of interaction between earthquakes by a factor up to 5-10 at distances of several tens of rupture length. This may explain certain long-range earthquake triggering and hydrological anomalies in wells and suggest to revisit the standard stress transfer calculations which use linear elasticity. We also show that the standard double-couple of forces representing an earthquake source leads to an opening of the corresponding fault plane, which suggests a mechanism for the non-zero isotropic component of the seismic moment tensor observed for some events.
The ETAS model is widely employed to model the spatio-temporal distribution of earthquakes, generally using spatially invariant parameters. We propose an efficient method for the estimation of spatially varying parameters, using the Expectation-Maximization (EM) algorithm and spatial Voronoi tessellation ensembles. We use the Bayesian Information Criterion (BIC) to rank inverted models given their likelihood and complexity and select the best models to finally compute an ensemble model at any location. Using a synthetic catalog, we also check that the proposed method correctly inverts the known parameters. We apply the proposed method to earthquakes included in the ANSS catalog that occurred within the time period 1981-2015 in a spatial polygon around California. The results indicate a significant spatial variation of the ETAS parameters. We find that the efficiency of earthquakes to trigger future ones (quantified by the branching ratio) positively correlates with surface heat flow. In contrast, the rate of earthquakes triggered by far-field tectonic loading or background seismicity rate shows no such correlation, suggesting the relevance of triggering possibly through fluid-induced activation. Furthermore, the branching ratio and background seismicity rate are found to be uncorrelated with hypocentral depths, indicating that the seismic coupling remains invariant of hypocentral depths in the study region. Additionally, triggering seems to be mostly dominated by small earthquakes. Consequently, the static stress change studies should not only focus on the Coulomb stress changes caused by specific moderate to large earthquakes but also account for the secondary static stress changes caused by smaller earthquakes.
The giant Sumatra-Andaman earthquake of December 26 2004 caused permanent deformations effects in a region of previously never observed extension. The GPS data from the world wide network of permanent IGS sites show significant coseismic displacements in an area exceeding 107 km^2. The effects of the permanent residual deformation field could be detected as far as Australia, the Phillipines and Japanese archipelagos, and, on the West, as far as the indian continent. The synthetic simulations of the coseismic displacement field obtained by means of a spherical model support the hypothesis that a considerable amount of energy has been released aseismically.