Quasistatic kinetic avalanches and self-organized criticality in deviatorically loaded granular media

Sheared granular matter and the empirical relations of seismicity

The frictional instability associated with earthquake initiation and earthquake dynamics is believed to be mainly controlled by the dynamics of fragmented rocks within the fault gauge. Principal features of the emerging seismicity (e.g., intermittent dynamics and broad time and/or energy scales) have been replicated by simple experimental setups, which involve a slowly driven slider on top of granular matter, for example. Yet these setups are often physically limited and might not allow one to determine the underlying nature of specific features and, hence, the universality and generality of the experimental observations. Here, we address this challenge by a numerical study of a spring-slider experiment based on two-dimensional discrete element method simulations, which allows us to control the properties of the granular matter and of the surface of the slider, for example. Upon quasistatic loading, stick-slip-type behavior emerges which is contrasted by a stable sliding regime at finite driving rates, in agreement with experimental observations. Across large parameter ranges for damping, interparticle friction, particle polydispersity, etc., the earthquake-like dynamics associated with the former regime results in several robust scale-free statistical features also observed in experiments. At first sight, these closely resemble the main empirical relations of tectonic seismicity at geological scales. This includes the Gutenberg-Richter distribution of event sizes, the Omori-Utsu-type decay of aftershock rates, as well as the aftershock productivity relation and broad recurrence time distributions. Yet, we show that the correlations associated with tectonic aftershocks are absent such that the origin of the Omori-Utsu relation, the aftershock productivity relation, and Båth's relation in the simulations is fundamentally different from the case of tectonic seismicity. This, we believe, is mainly due to a lack of macroscale relaxation processes that are closely tied to the generation of real aftershocks. We argue that the same is true for previous laboratory experiments.

Signatures of the spatial extent of plastic events in the yielding transition in amorphous solids

Amorphous solids are yield stress materials that flow when a sufficient load is applied. Their flow consists of periods of elastic loading interrupted by rapid stress drops, or avalanches, coming from microscopic rearrangements known as shear transformations (STs). Here we show that the spatial extent of avalanches in a steadily sheared amorphous solid has a profound effect on the distribution of local residual stresses that in turn determines the stress drop statistics. As reported earlier, the most unstable sites are located in a flat "plateau" region that decreases with system size. While the entrance into the plateau is set by the lower cutoff of the mechanical noise produced by individual STs, the departure from the usually assumed power-law (pseudogap) form of the residual stress distribution comes from far field effects related to spatially extended rearrangements. Interestingly, we observe that the average residual stress of the weakest sites is located in an intermediate power-law regime between the pseudogap and the plateau regimes, whose exponent decreases with system size. Our findings imply a new scaling relation linking the exponents characterizing the avalanche size and residual stress distributions.