No Arabic abstract
We present the condensation method that exploits the heterogeneity of the probability distribution functions (PDF) of event locations to improve the spatial information content of seismic catalogs. The method reduces the size of seismic catalogs while improving the access to the spatial information content of seismic catalogs. The PDFs of events are first ranked by decreasing location errors and then successively condensed onto better located and lower variance event PDFs. The obtained condensed catalog attributes different weights to each event, providing an optimal spatial representation with respect to the spatially varying location capability of the seismic network. Synthetic tests on fractal distributions perturbed with realistic location errors show that condensation improves spatial information content of the original catalog. Applied to Southern California seismicity, the new condensed catalog highlights major mapped fault traces and reveals possible additional structures while reducing the catalog length by ~25%. The condensation method allows us to account for location error information within a point based spatial analysis. We demonstrate this by comparing the multifractal properties of the condensed catalog locations with those of the original catalog. We evidence different spatial scaling regimes characterized by distinct multifractal spectra and separated by transition scales. We interpret the upper scale as to agree with the thickness of the brittle crust, while the lower scale (2.5km) might depend on the relocation procedure. Accounting for these new results, the Epidemic Type Aftershock Model formulation suggests that, contrary to previous studies, large earthquakes dominate the earthquake triggering process. This implies that the limited capability of detecting small magnitude events cannot be used to argue that earthquakes are unpredictable in general.
We report a similarity of fluctuations in equilibrium critical phenomena and non-equilibrium systems, which is based on the concept of natural time. The world-wide seismicity as well as that of San Andreas fault system and Japan are analyzed. An order parameter is chosen and its fluctuations relative to the standard deviation of the distribution are studied. We find that the scaled distributions fall on the same curve, which interestingly exhibits, over four orders of magnitude, features similar to those in several equilibrium critical phenomena (e.g., 2D Ising model) as well as in non-equilibrium systems (e.g., 3D turbulent flow).
The rise in the frequency of anthropogenic earthquakes due to deep fluid injections is posing serious economic, societal, and legal challenges to geo-energy and waste-disposal projects. We propose an actuarial approach to mitigate this risk, first by defining an autonomous decision-making process based on an adaptive traffic light system (ATLS) to stop risky injections, and second by quantifying a cost of public safety based on the probability of an injection-well being abandoned. The ATLS underlying statistical model is first confirmed to be representative of injection-induced seismicity, with examples taken from past reservoir stimulation experiments (mostly from Enhanced Geothermal Systems, EGS). Then the decision strategy is formalized: Being integrable, the model yields a closed-form ATLS solution that maps a risk-based safety standard or norm to an earthquake magnitude not to exceed during stimulation. Finally, the EGS levelized cost of electricity (LCOE) is reformulated in terms of null expectation, with the cost of abandoned injection-well implemented. We find that the price increase to mitigate the increased seismic risk in populated areas can counterbalance the heat credit. However this public safety cost disappears if buildings are based on earthquake-resistant designs or if a more relaxed risk safety standard or norm is chosen.
The correlation properties of the magnitudes of a time series (sometimes called volatility) are associated with nonlinear and multifractal properties and have been applied in a great variety of fields. Here, we have obtained analytically the expression of the autocorrelation of the magnitude series of a linear Gaussian noise as a function of its correlation as well as several analytical relations involving them. For both, models and natural signals, the deviation from these equations can be used as an index of non-linearity that can be applied to relatively short records and that does not require the presence of scaling in the time series under study. We apply this approach to show that the heart-beat records during rest show higher non-linearities than the records of the same subject during moderate exercise. This behavior is also achieved on average for the analyzed set of 10 semiprofessional soccer players. This result agrees with the fact that other measures of complexity are dramatically reduced during exercise and can shed light on its relationship with the withdrawal of parasympathetic tone and/or the activation of sympathetic activity during physical activity.
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.
In line of the intermediate-term monitoring of seismic activity aimed at prediction of the world largest earthquakes the seismic dynamics of the Earths lithosphere is analysed as a single whole, which is the ultimate scale of the complex hierarchical non-linear system. The present study demonstrates that the lithosphere does behave, at least in intermediate-term scale, as non-linear dynamic system that reveals classical symptoms of instability at the approach of catastrophe, i.e., mega-earthquake. These are: (i) transformation of magnitude distribution, (ii) spatial redistribution of seismic activity, (iii) rise and acceleration of activity, (iv) change of dependencies across magnitudes of different types, and other patterns of collective behaviour. The observed global scale seismic behaviour implies the state of criticality of the Earth lithosphere in the last decade.