Asperity characteristics of the Olami-Feder-Christensen model of earthquakes

Model for tectonic tremors: Enduring events, moment rate spectrum, and moment-duration scaling

Numerous attempts have been made to obtain earthquake statistics from a theoretical-physics perspective, but these studies mostly involve regular earthquakes. In recent years, a new category of earthquakes, referred to as slow earthquakes, has been discovered. Slow earthquakes emit only weak or no seismic signals and have different statistics than regular earthquakes. Here we propose a physical model for the tremor, which is a type of slow earthquake, introducing two competing timescales in a cellular automaton model. The proposed model reproduces some observation results for tremors, such as enduring events, moment-duration scaling, size distribution, and the power spectrum of the moment rate function.

Mapping heterogeneities through avalanche statistics

Avalanche statistics of various threshold-activated dynamical systems are known to depend on the magnitude of the drive, or stress, on the system. Such dependences exist for earthquake size distributions, in sheared granular avalanches, laboratory-scale fracture and also in the outage statistics of power grids. In this work, we model threshold-activated avalanche dynamics and investigate the time required to detect local variations in the ability of model elements to bear stress. We show that the detection time follows a scaling law where the scaling exponents depend on whether the feature that is sought is either weaker, or stronger, than its surroundings. We then look at earthquake data from Sumatra and California, demonstrate the trade-off between the spatial resolution of a map of earthquake exponents (i.e. the b-values of the Gutenberg-Richter Law) and the accuracy of those exponents, and suggest a means to maximize both.This article is part of the theme issue 'Statistical physics of fracture and earthquakes'.

Statistical properties of the one-dimensional Burridge-Knopoff model of earthquakes obeying the rate- and state-dependent friction law

Statistical properties of the one-dimensional spring-block (Burridge-Knopoff) model of earthquakes obeying the rate- and state-dependent friction law are studied by extensive computer simulations. The quantities computed include the magnitude distribution, the rupture-length distribution, the main shock recurrence-time distribution, the seismic-time correlations before and after the main shock, the mean slip amount, and the mean stress drop at the main shock, etc. Events of the model can be classified into two distinct categories. One tends to be unilateral with its epicenter located at the rim of the rupture zone of the preceding event, while the other tends to be bilateral with enhanced "characteristic" features resembling the so-called "asperity." For both types of events, the distribution of the rupture length L_{r} exhibits an exponential behavior at larger sizes, ≈exp[-L_{r}/L_{0}] with a characteristic "seismic correlation length" L_{0}. The mean slip as well as the mean stress drop tends to be rupture-length independent for larger events. The continuum limit of the model is examined, where the model is found to exhibit pronounced characteristic features. In the continuum limit, the characteristic rupture length L_{0} is estimated to be ∼100 [km]. This means that, even in a hypothetical homogenous infinite fault, events cannot be indefinitely large in the exponential sense, the upper limit being of order ∼10^{3} kilometers. Implications to real seismicity are discussed.

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.

Record-breaking avalanches in driven threshold systems

Record-breaking avalanches generated by the dynamics of several driven nonlinear threshold models are studied. Such systems are characterized by intermittent behavior, where a slow buildup of energy is punctuated by an abrupt release of energy through avalanche events, which usually follow scale-invariant statistics. From the simulations of these systems it is possible to extract sequences of record-breaking avalanches, where each subsequent record-breaking event is larger in magnitude than all previous events. In the present work, several cellular automata are analyzed, among them the sandpile model, the Manna model, the Olami-Feder-Christensen (OFC) model, and the forest-fire model to investigate the record-breaking statistics of model avalanches that exhibit temporal and spatial correlations. Several statistical measures of record-breaking events are derived analytically and confirmed through numerical simulations. The statistics of record-breaking avalanches for the four models are compared to those of record-breaking events extracted from the sequences of independent and identically distributed (i.i.d.) random variables. It is found that the statistics of record-breaking avalanches for the above cellular automata exhibit behavior different from that observed for i.i.d. random variables, which in turn can be used to characterize complex spatiotemporal dynamics. The most pronounced deviations are observed in the case of the OFC model with a strong dependence on the conservation parameter of the model. This indicates that avalanches in the OFC model are not independent and exhibit spatiotemporal correlations.