Earthquake interevent time distribution for induced micro-, nano-, and picoseismicity

Self-similar properties of avalanche statistics in a simple turbulent model

In this paper, we consider a simplified model of turbulence for large Reynolds numbers driven by a constant power energy input on large scales. In the statistical stationary regime, the behaviour of the kinetic energy is characterized by two well-defined phases: a laminar phase where the kinetic energy grows linearly for a (random) time [Formula: see text] followed by abrupt avalanche-like energy drops of sizes [Formula: see text] due to strong intermittent fluctuations of energy dissipation. We study the probability distribution [Formula: see text] and [Formula: see text] which both exhibit a quite well-defined scaling behaviour. Although [Formula: see text] and [Formula: see text] are not statistically correlated, we suggest and numerically checked that their scaling properties are related based on a simple, but non-trivial, scaling argument. We propose that the same approach can be used for other systems showing avalanche-like behaviour such as amorphous solids and seismic events. This article is part of the theme issue 'Scaling the turbulence edifice (part 1)'.

Controlling crackling dynamics by triggering low-intensity avalanches

We examine the effect of small, spatially localized excitations applied periodically in different manners, on the crackling dynamics of a brittle crack driven slowly in a heterogeneous solid. When properly adjusted, these excitations are observed to radically modify avalanche statistics and considerably limit the magnitude of the largest events. Surprisingly, this does not require information on the front loading state at the time of excitation; applying it either at a random location or at the most loaded point gives the same results. Subsequently, we unravel how the excitation amplitude, spatial extent, and frequency govern the effect. We find that the excitation efficiency is ruled by a single reduced parameter, namely the injected power per unit front length; the suppression of extreme avalanches is maximum at a well-defined optimal value of this control parameter. analysis opens another way to control the largest events in crackling dynamics. Beyond fracture problems, it may be relevant for crackling systems described by models of the same universality class, such as the wetting of heterogeneous substrates or magnetic walls in amorphous magnets.

Earthquakes and entropy: Characterization of occurrence of earthquakes in southern Spain and Alboran Sea

We propose the use of entropy, H, as an indicator of the equilibrium state of a seismically active region (seismic system). The relationship between an increase in H and the occurrence of a great earthquake in a study area can be predicted by acknowledging the irreversible transition of a system. From this point of view, the seismic system evolves from an unstable initial state (due to external stresses) to another, where the stresses have dropped after the earthquake occurred. It is an irreversible transition that entails an increase in entropy. Five seismic episodes were analyzed in the south of the Iberian Peninsula, the Alboran Sea (Mediterranean Sea), and the North of Morocco: two of them of moderate-high magnitude (Al Hoceima, 2004 and 2016) and three of them of moderate-low magnitude (Adra, 1993-1994; Moron, 2007; and Torreperogil, 2012-2013). The results are remarkably in line with the theoretical forecasts; in other words: an earthquake, understood as an irreversible transition, must suppose an increase in entropy.

Interevent-time distribution and aftershock frequency in non-stationary induced seismicity

The initial footprint of an earthquake can be extended considerably by triggering of clustered aftershocks. Such earthquake-earthquake interactions have been studied extensively for data-rich, stationary natural seismicity. Induced seismicity, however, is intrinsically inhomogeneous in time and space and may have a limited catalog of events; this may hamper the distinction between human-induced background events and triggered aftershocks. Here we introduce a novel Gamma Accelerated-Failure-Time model for efficiently analyzing interevent-time distributions in such cases. It addresses the spatiotemporal variation and quantifies, per event, the probability of each event to have been triggered. Distentangling the obscuring aftershocks from the background events is a crucial step to better understand the causal relationship between operational parameters and non-stationary induced seismicity. Applied to the Groningen gas field in the North of the Netherlands, our model elucidates geological and operational drivers of seismicity and has been used to test for aftershock triggering. We find that the hazard rate in Groningen is indeed enhanced after each event and conclude that aftershock triggering cannot be ignored. In particular we find that the non-stationary interevent-time distribution is well described by our Gamma model. This model suggests that 27.0(± 8.5)% of the recorded events in the Groningen field can be attributed to triggering.

Dynamics of cascades on burstiness-controlled temporal networks

Burstiness, the tendency of interaction events to be heterogeneously distributed in time, is critical to information diffusion in physical and social systems. However, an analytical framework capturing the effect of burstiness on generic dynamics is lacking. Here we develop a master equation formalism to study cascades on temporal networks with burstiness modelled by renewal processes. Supported by numerical and data-driven simulations, we describe the interplay between heterogeneous temporal interactions and models of threshold-driven and epidemic spreading. We find that increasing interevent time variance can both accelerate and decelerate spreading for threshold models, but can only decelerate epidemic spreading. When accounting for the skewness of different interevent time distributions, spreading times collapse onto a universal curve. Our framework uncovers a deep yet subtle connection between generic diffusion mechanisms and underlying temporal network structures that impacts a broad class of networked phenomena, from spin interactions to epidemic contagion and language dynamics.

Correlating the force network evolution and dynamics in slider experiments

The experiments involving a slider moving on top of granular media consisting of photoelastic particles in two dimensions have uncovered elaborate dynamics that may vary from continuous motion to crackling, periodic motion, and stick-slip type of behavior. We establish that there is a clear correlation between the slider dynamics and the response of the force network that spontaneously develop in the granular system. This correlation is established by application of the persistence homology that allows for formulation of objective measures for quantification of time-dependent force networks. We find that correlation between the slider dynamics and the force network properties is particularly strong in the dynamical regime characterized by well-defined stick-slip type of dynamics.

On interevent time distributions of avalanche dynamics

Physical systems characterized by stick-slip dynamics often display avalanches. Regardless of the diversity of their microscopic structure, these systems are governed by a power-law distribution of avalanche size and duration. Here we focus on the interevent times between avalanches and show that, unlike their distributions of size and duration, the interevent time distributions are able to distinguish different mechanical states of the system. We use experiments on granular systems and numerical simulations of emulsions to show that systems having the same probability distribution for avalanche size and duration can have different interevent time distributions. Remarkably, these interevent time distributions look similar to those for earthquakes and, if different from an exponential, are indirect evidence of non trivial space-time correlations among avalanches. Our results therefore indicate that interevent time statistics are essential to characterise the dynamics of avalanches.

Seismiclike organization of avalanches in a driven long-range elastic string as a paradigm of brittle cracks

Crack growth in heterogeneous materials sometimes exhibits crackling dynamics, made of successive impulselike events with specific scale-invariant time and size organization reminiscent of earthquakes. Here, we examine this dynamics in a model which identifies the crack front with a long-range elastic line driven in a random potential. We demonstrate that, under some circumstances, fracture grows intermittently, via scale-free impulse organized into aftershock sequences obeying the fundamental laws of statistical seismology. We examine the effects of the driving rate and system overall stiffness (unloading factor) onto the scaling exponents and cutoffs associated with the time and size organization. We unravel the specific conditions required to observe a seismiclike organization in the crack propagation problem. Beyond failure problems, implications of these results to other crackling systems are finally discussed.

Seismicity in sheared granular matter

We report on experiments investigating the dynamics of a slider that is pulled by a spring across a granular medium consisting of a vertical layer of photoelastic disks. The motion proceeds through a sequence of discrete events, analogous to seismic shocks, in which elastic energy stored in the spring is rapidly released. We measure the statistics of several properties of the individual events: the energy loss in the spring, the duration of the movement, and the temporal profile of the slider motion. We also study certain conditional probabilities and the statistics of mainshock-aftershock sequences. At low driving rates, we observe crackling with Omori-Utsu, Båth, and waiting time laws similar to those observed in seismic dynamics. At higher driving rates, where the sequence of events shows strong periodicity, we observe scaling laws and asymmetrical event shapes that are clearly distinguishable from those in the crackling regime.

Crackling to periodic dynamics in granular media

We study the local and global dynamics of sheared granular materials in a stick-slip experiment, using a slider and a spring. The system crackles, with intermittent slip avalanches, or exhibits irregular or periodic dynamics, depending on the shear rate and loading stiffness. The global force on the slider during shearing captures the transitions from the crackling to the periodic regime. We deduce a dynamic phase diagram as a function of the shear rate and the loading stiffness and associated scaling laws. Using photoelastic particles, we also capture the grain-scale stress evolution, and investigate the microscopic behavior in the different regimes.

Crack growth in heterogeneous brittle solids: intermittency, crackling and induced seismicity

Crack growth is the basic mechanism leading to the failure of brittle materials. Engineering addresses this problem within the framework of continuum mechanics, which links deterministically the crack motion to the applied loading. Such an idealization, however, fails in several situations and in particular cannot capture the highly erratic (earthquake-like) dynamics sometimes observed in slowly fracturing heterogeneous solids. Here, we examine this problem by means of innovative experiments of crack growth in artificial rocks of controlled microstructure. The dynamical events are analysed at both global and local scales, from the time fluctuation of the spatially averaged crack speed and the induced acoustic emission, respectively. Their statistics are characterized and compared with the predictions of a recent approach mapping fracture onset to the depinning of an elastic interface. Finally, the overall time-size organization of the events is characterized to shed light on the mechanisms underlying the scaling laws observed in seismology.This article is part of the theme issue 'Statistical physics of fracture and earthquakes'.

Aftershock sequences and seismic-like organization of acoustic events produced by a single propagating crack

Brittle fractures of inhomogeneous materials like rocks, concrete, or ceramics are of two types: Nominally brittle and driven by the propagation of a single dominant crack or quasi-brittle and resulting from the accumulation of many microcracks. The latter goes along with acoustic noise, whose analysis has revealed that events form aftershock sequences obeying characteristic laws reminiscent of those in seismology. Yet, their origin lacks explanation. Here we show that such a statistical organization is not only specific to the multi-cracking situations of quasi-brittle failure and seismology, but also rules the acoustic events produced by a propagating crack. This simpler situation has permitted us to relate these laws to the overall scale-free distribution of inter-event time and energy and to uncover their selection by the crack speed. These results provide a comprehensive picture of how acoustic events are organized upon material failure in the most fundamental of fracture states: single propagating cracks.

The multiple microcracking events underlying damage in inhomogeneous brittle materials form characteristic aftershocks sequences obeying similar laws to those in seismology. Here, Barés et al. evidence and explain the same organization in the acoustic noise produced by a single propagating crack.

Triggering Processes in Rock Fracture

We study triggering processes in triaxial compression experiments under a constant displacement rate on sandstone and granite samples using spatially located acoustic emission events and their focal mechanisms. We present strong evidence that event-event triggering plays an important role in the presence of large-scale or macrocopic imperfections, while such triggering is basically absent if no significant imperfections are present. In the former case, we recover all established empirical relations of aftershock seismicity including the Gutenberg-Richter relation, a modified version of the Omori-Utsu relation and the productivity relation-despite the fact that the activity is dominated by compaction-type events and triggering cascades have a swarmlike topology. For the Gutenberg-Richter relations, we find that the b value is smaller for triggered events compared to background events. Moreover, we show that triggered acoustic emission events have a focal mechanism much more similar to their associated trigger than expected by chance.

Avalanche structural rearrangement through cracking-healing in weakly stressed cold dusty plasma liquids

We experimentally investigate the spatiotemporal dynamical behaviors of the avalanche structural rearrangement through micro-cracking-healing in weakly stressed cold dusty plasma liquids, and the kinetic origins for their different spatial and temporal classifications. The crystalline ordered domains can be cracked or temporarily sustain and transfer the weak stress to remote regions for cracking-healing. It is found that cracking sites form a fractal network with cluster size following power law distribution in the xyt space. The histograms of the persistent times for sustaining regional ordered and disordered structure, the temporal cracking burst width, and quiescent time between two bursts all follow power law decays with fast descending tails. Cracking can be classified into a single temporal burst with simple line like spatial patterns and the successive cracking fluctuation with densely packed cracking clusters. For an ordered region, whether the Burgers vectors of the incoming dislocations from the boundary allow direct dislocation reduction is the key for the above two classifications through cracking a large ordered domain into medium scale corotating ordered domains or small patches. The low regional structural order at the end of a cracking burst can be regarded as an alarm for predicting the short quiescent period before the next cracking burst.

Predicting failure: acoustic emission of berlinite under compression

Acoustic emission has been measured and statistical characteristics analyzed during the stress-induced collapse of porous berlinite, AlPO4, containing up to 50 vol% porosity. Stress collapse occurs in a series of individual events (avalanches), and each avalanche leads to a jerk in sample compression with corresponding acoustic emission (AE) signals. The distribution of AE avalanche energies can be approximately described by a power law p(E)dE = E(-ε)dE (ε ~ 1.8) over a large stress interval. We observed several collapse mechanisms whereby less porous minerals show the superposition of independent jerks, which were not related to the major collapse at the failure stress. In highly porous berlinite (40% and 50%) an increase of energy emission occurred near the failure point. In contrast, the less porous samples did not show such an increase in energy emission. Instead, in the near vicinity of the main failure point they showed a reduction in the energy exponent to ~ 1.4, which is consistent with the value reported for compressed porous systems displaying critical behavior. This suggests that a critical avalanche regime with a lack of precursor events occurs. In this case, all preceding large events were 'false alarms' and unrelated to the main failure event. Our results identify a method to use pico-seismicity detection of foreshocks to warn of mine collapse before the main failure (the collapse) occurs, which can be applied to highly porous materials only.

Unveiling the complex organization of recurrent patterns in spiking dynamical systems

Complex systems displaying recurrent spike patterns are ubiquitous in nature. Understanding the organization of these patterns is a challenging task. Here we study experimentally the spiking output of a semiconductor laser with feedback. By using symbolic analysis we unveil a nontrivial organization of patterns, revealing serial spike correlations. The probabilities of the patterns display a well-defined, hierarchical and clustered structure that can be understood in terms of a delayed model. Most importantly, we identify a minimal model, a modified circle map, which displays the same symbolic organization. The validity of this minimal model is confirmed by analyzing the output of the forced laser. Since the circle map describes many dynamical systems, including neurons and cardiac cells, our results suggest that similar correlations and hierarchies of patterns can be found in other systems. Our findings also pave the way for optical neurons that could provide a controllable set up to mimic neuronal activity.

Contextual analysis framework for bursty dynamics

To understand the origin of bursty dynamics in natural and social processes we provide a general analysis framework in which the temporal process is decomposed into subprocesses and then the bursts in subprocesses, called contextual bursts, are combined to collective bursts in the original process. For the combination of subprocesses, it is required to consider the distribution of different contexts over the original process. Based on minimal assumptions for interevent time statistics, we present a theoretical analysis for the relationship between contextual and collective interevent time distributions. Our analysis framework helps to exploit contextual information available in decomposable bursty dynamics.