Earthquake dynamics driven by a viscous fluid

The influence of the brittle-ductile transition zone on aftershock and foreshock occurrence

Aftershock occurrence is characterized by scaling behaviors with quite universal exponents. At the same time, deviations from universality have been proposed as a tool to discriminate aftershocks from foreshocks. Here we show that the change in rheological behavior of the crust, from velocity weakening to velocity strengthening, represents a viable mechanism to explain statistical features of both aftershocks and foreshocks. More precisely, we present a model of the seismic fault described as a velocity weakening elastic layer coupled to a velocity strengthening visco-elastic layer. We show that the statistical properties of aftershocks in instrumental catalogs are recovered at a quantitative level, quite independently of the value of model parameters. We also find that large earthquakes are often anticipated by a preparatory phase characterized by the occurrence of foreshocks. Their magnitude distribution is significantly flatter than the aftershock one, in agreement with recent results for forecasting tools based on foreshocks.

Universal avalanche statistics and triggering close to failure in a mean-field model of rheological fracture

The hypothesis of critical failure relates the presence of an ultimate stability point in the structural constitutive equation of materials to a divergence of characteristic scales in the microscopic dynamics responsible for deformation. Avalanche models involving critical failure have determined common universality classes for stick-slip processes and fracture. However, not all empirical failure processes exhibit the trademarks of criticality. The rheological properties of materials introduce dissipation, usually reproduced in conceptual models as a hardening of the coarse grained elements of the system. Here, we investigate the effects of transient hardening on (i) the activity rate and (ii) the statistical properties of avalanches. We find the explicit representation of transient hardening in the presence of generalized viscoelasticity and solve the corresponding mean-field model of fracture. In the quasistatic limit, the accelerated energy release is invariant with respect to rheology and the avalanche propagation can be reinterpreted in terms of a stochastic counting process. A single universality class can be defined from such analogy, and all statistical properties depend only on the distance to criticality. We also prove that interevent correlations emerge due to the hardening-even in the quasistatic limit-that can be interpreted as "aftershocks" and "foreshocks."

Power-law rheology controls aftershock triggering and decay

The occurrence of aftershocks is a signature of physical systems exhibiting relaxation phenomena. They are observed in various natural or experimental systems and usually obey several non-trivial empirical laws. Here we consider a cellular automaton realization of a nonlinear viscoelastic slider-block model in order to infer the physical mechanisms of triggering responsible for the occurrence of aftershocks. We show that nonlinear viscoelasticity plays a critical role in the occurrence of aftershocks. The model reproduces several empirical laws describing the statistics of aftershocks. In case of earthquakes, the proposed model suggests that the power-law rheology of the fault gauge, underlying lower crust, and upper mantle controls the decay rate of aftershocks. This is verified by analysing several prominent aftershock sequences for which the rheological properties of the underlying crust and upper mantle were established.

Interoccurrence time statistics in the two-dimensional Burridge-Knopoff earthquake model

We have numerically investigated statistical properties of the so-called interoccurrence time or the waiting time, i.e., the time interval between successive earthquakes, based on the two-dimensional (2D) spring-block (Burridge-Knopoff) model, selecting the velocity-weakening property as the constitutive friction law. The statistical properties of frequency distribution and the cumulative distribution of the interoccurrence time are discussed by tuning the dynamical parameters, namely, a stiffness and frictional property of a fault. We optimize these model parameters to reproduce the interoccurrence time statistics in nature; the frequency and cumulative distribution can be described by the power law and Zipf-Mandelbrot type power law, respectively. In an optimal case, the b value of the Gutenberg-Richter law and the ratio of wave propagation velocity are in agreement with those derived from real earthquakes. As the threshold of magnitude is increased, the interoccurrence time distribution tends to follow an exponential distribution. Hence it is suggested that a temporal sequence of earthquakes, aside from small-magnitude events, is a Poisson process, which is observed in nature. We found that the interoccurrence time statistics derived from the 2D BK (original) model can efficiently reproduce that of real earthquakes, so that the model can be recognized as a realistic one in view of interoccurrence time statistics.