# Breakdown of Scaling and Friction Weakening in the Critical Granular Flow

###### Abstract

The way granular materials response to an applied shear stress is of the utmost relevance to both human activities and natural environment. One of the their most intriguing and less understood behavior, is the stick-instability, whose most dramatic manifestation are earthquakes, ultimately governed by the dynamics of rocks and debris jammed within the fault gauge. Many of the features of earthquakes, i.e. intermittency, broad times and energy scale involved, are mimicked by a very simple experimental set-up, where small beads of glass under load are slowly sheared by an elastic medium. Analyzing data from long lasting experiments, we identify a critical dynamical regime, that can be related to known theoretical models used for ”crackling-noise” phenomena. In particular, we focus on the average shape of the slip velocity, observing a ”breakdown of scaling”: while small slips show a self-similar shape, large does not, in a way that suggests the presence of subtle inertial effects within the granular system. In order to characterise the crossover between the two regimes, we investigate the frictional response of the system, which we trat as a stochastic quantity. Computing different averages, we evidence a weakening effect, whose Stribeck threshold velocity can be related to the aforementioned breaking of scaling.

###### pacs:

45.70.-n 45.70.Ht 05.65.+b 07.79.Sp, 89.75.Da### .1 Motivations

The way a granular bed responses to an applied shear stress reveals many of the peculiarities of this poorly comprehended ”state” of matter Jaeger et al. (2000). When a granular bed is sheared slowly enough by an elastic medium driven at constant velocity, nor the shear stress neither the shear rate can be directly controlled from outside ann (2016). Rather, the system sets itself in a state at the edge between jamming and mobility Dalton and Corcoran (2001); Ciamarra et al. (2010); Dalton et al. (2005); Baldassarri et al. (2006); Petri et al. (2008); Geller et al. (2015); Zadeh et al. (2018), exhibiting intermittent flow also called stick-slip. This is an instance, among many others, of phenomena displaying intermittent and erratic activity, in the form of bursts, or avalanches, characterized by wild fluctuations of physical quantities, and for this reason named crackling noise Sethna et al. (2001). Examples include earthquakes Main (1996), fractures Petri et al. (1994), structural phase transitions Cannelli et al. (1993) and plastic deformation Dimiduk et al. (2006). These diverse phenomena share several common statistical features. In particular physical quantities display often long range correlations and self-similar distributions, i.e. power laws, over a wide range of values. Such properties are usually ascribed to the vicinity of some critical transition Bak et al. (1988); Sethna et al. (2001), which in granular media could be the jamming transition Dauchot et al. (2005). Consistently, critical transitions bring about the existence of universality classes: systems that are microscopically very different, can display similar and universal statistical properties in their critical dynamics. Within this spirit we have designed an experimental setup suitable to observe such an irregular granular dynamics Dalton et al. (2005), characterized by critical fluctuations and reminiscent of that displayed by the aforementioned wide class of physical systems.

In order to compare different systems exhibiting critical dynamics, several quantities can be analyzed. Recent literature witnesses a surge of interest for the average avalanche (or burst) shape (or profile). Introduced in the context of Barkhausen noise in ferromagnetic materials Kuntz and Sethna (2000), the average avalanche shape can provide a much sharper tool to test theory against experiments than the simple comparison of critical exponents characterizing probability distributions. As shown for simple stochastic processes, the geometrical and scaling properties of the average shape of a fluctuation depends on the temporal correlations of the dynamics Baldassarri et al. (2003); Colaiori et al. (2004); Colaiori (2008). Such observation has allowed, for instance, to evidence a (negative) effective mass in magnetic domain walls Zapperi et al. (2005).

Average avalanche shapes have been investigated for a variety of materials, well beyond magnetic systems Papanikolaou et al. (2011). Among the others: dislocation avalanches in plastically deformed intermetallic compounds Chrzan and Mills (1994) and in gold and niobium crystals Sparks et al. (2017); stress drop avalanches at the yielding transition in metallic glasses Antonaglia et al. (2014) and, via numerical simulations, in amorphous systems Liu et al. (2016); Lagogianni et al. (2018); bursts of load redistribution in heterogeneous materials under a constant external load Danku and Kun (2013). Many biological studies have also measured average burst shape in cortical bursting activity Roberts et al. (2014); Wikstro et al. (2015), in transport processes in living cells Wang et al. (2013), as well as in ants Gallotti and Chialvo (2018) and human Chialvo et al. (2015) activity. Many other bursty dynamics have been investigated by means of this general tool, as stellar processes Sheikh et al. (2016) or Earth’s magnetospheric dynamics Consolini et al. (2008), and earthquakes Mehta et al. (2006). The dependence of the avalanche shape from the interaction range has been studied in elastic depinning models Laurson et al. (2013).

In this paper we acquire and analyze for the first time the average shapes of slip velocity and of friction force in a sheared granular system, directly in the deep critical phase where it displays intermittent flow. Our findings also shed light on apparently contradictory recent observations Barés et al. (2017); Denisov et al. (2017), and supply new essential elements to improve the formulation of new and more effective dynamical models, with important impact on the understanding of related natural and technological issues.

### .2 Introduction

We study the stick-slip dynamics at the level of the single slip event, as illustrated in Fig. 1. The left panel reproduces the angular velocity, during a slip, of a slider that rotates while in contact with the granular bed. The middle panel shows the corresponding frictional torque experienced by the slider. The motion can be described as a function of the internal avalanche time , which starts at the beginning of the slip and ends when the system sticks. Each slip event has its own duration and its size (the grey area in Fig. 1, left panel). The average velocity shape is performed considering many slips with the same total duration as function of internal avalanche time . A similar averaging procedure is followed to obtain the average friction shape: i.e. the average friction torque exerted by the granular medium at the internal time during a slip event of total duration . The right panel of Fig. 1 shows the intricate, complex relation between friction and velocity during the intermittent, stick-slip dynamics.

In our study we observe the existence of a cross-over from small to large slips. We identify it as a breakdown of the critical scaling and show that such transition is in turn related to a change in the frictional properties of the system. Specifically, we find that the average velocity of the cross-over avalanches corresponds to a characteristic value marking a dynamical transition from weakening to thickening frictional behavior of the system. Average shape for avalanches of stress drop Denisov et al. (2017) and energy drop rate Barés et al. (2017) have been recently investigated in slow but continuous flow, where velocity never drops to zero and the stress is the relevant fluctuating quantity. While in Denisov et al. (2017) average avalanches have been found to display symmetrical and self-similar shapes, in Barés et al. (2017) these properties have been observed only in avalanches sufficiently small. Our investigations, conducted in the critical state, contribute to clarify the origin of these contradictory behavior observed in a different situation.

### .3 Experimental set up

The experimental set up (see Fig. 8) is similar to that employed in Dalton et al. (2005); Baldassarri et al. (2006); Petri et al. (2008); Leoni et al. (2010); Ciamarra et al. (2012) and described in more detail in the Appendix. The apparatus consists of an assembly of glass spheres laying in an annular channel and sheared by a horizontally rotating top plate driven by a motor. The instantaneous angular position of the plate and of the motor, respectively and are acquired by means of two optical encoders.

The plate is coupled to the motor through a soft torsion spring of elastic constant , so the instantaneous frictional torque, , exerted by the granular medium can be derived from the equation of motion for the plate:

(1) |

where is the inertia of the plate-axis system. The motor angular speed is kept constant, so that , but the interaction between the top plate and the granular medium is crucial in determining the instantaneous plate velocity, leading to different possible regimes. When both the driving speed and spring constant are low enough the critical dynamics, in which the plate performs highly irregular and intermittent motion, is approached.

### .4 Scaling analysis

We have performed long experimental runs in the critical, stick-slip, regime measuring the angular coordinate of the plate , from which we have derived the plate angular velocity . We have collected statistics for a large number of avalanches: the distribution of corresponding durations and sizes are shown in Fig. 2. Both distributions exhibit a slow decay, roughly close to a power law, terminating by a bulging cutoff for large sizes. Similar broad distributions are shared by other quantities, e.g. the plate velocity, Baldassarri et al. (2006) (not shown).

As recalled in the introduction, power law decay in distributions are generally considered the hallmark for criticality. If this is the correct scenario, one should observe self-similar scaling relations in average quantities too. In particular, we consider the average shape of velocity during an avalanche of a fixed duration, defined as:

where is the plate velocity during the observed avalanche of duration , whose total number is , and is the internal time within the slip: . Although the average velocity shape depends on both and , criticality should imply that an invariant function exists, such that it can be expressed as:

(2) |

The function determines how the average event size scales with respect to the slip duration . In fact, integrating the above equation with respect to one gets:

(3) |

(where without loss of generality we have assumed ). The function represents the average invariant pulse shape, which is expected not to depend on the slip duration and can be computed via the above equations as

(4) |

The previous scaling scenario is produced by several theoretical models for critical dynamics. One paradigmatic model for crackling noise is the so called ABBM model Alessandro et al. (1990), proposed to describe the intermittent statistics of electric noise recorded during hysteresis loops in ferromagnetic materials (Barkhausen noise). It is simple enough to allow exact analytical results Alessandro et al. (1990); Feller (1951); Colaiori (2008); Papanikolaou et al. (2011); Dobrinevski et al. (2012), and it predicts power law distributions for avalanche sizes and durations, as well as parabolic average avalanche shape, in the scaling regime. In the conclusive section we will discuss the connections between this model and what we observe in our study.

### .5 Average shape of slip velocity

To investigate the properties of the average pulse shape, and to test its invariance and the scaling hypothesis, we have divided all the avalanches observed in the experiments into classes according to their duration (see Appendix). Figure 3 (main panel) shows the slip size as function of its duration for all the slips considered in the statistics, and the different colors correspond to the different classes of duration. For each class, , we have computed the average slip size and duration , and the average velocity measured as function of the internal time . According to Eq. (4), in order to obtain , this average velocity has then been normalized to the ratio .

The resulting average shapes for each class of duration are shown in Fig. 4 (light, grey points in 3, corresponding to very short slips at the limit of the system resolution, have been discarded). All classes exhibit comparable values of the rescaled maximum velocity implying that longer avalanches are also faster. However, rescaled average shapes unveil that there are two kinds of avalanches. Some of them, say short, have the shape described by a unique function , visible in Fig. 4 (left panel). That is, their size and duration are related by the well defined scaling law Eq. 2. On the contrary, the average velocity shapes of long avalanches (right panel) change with the duration and cannot be reduced to a universal form by a homogeneous rescaling of the variables. Moreover, they do not display the almost symmetric shape characterizing small avalanches.

As anticipated, Barés and coworkers Barés et al. (2017) have recently measured the average shape of stress drop rate avalanches in a bidimensional granular system driven at constant shear rate. Similarly to the present findings, they observe that larger slips develop left asymmetries. They have hypothesized a possible role of the static friction between particles and supporting glass, and of nonlinear elasticity, given by the relatively soft nature of the grain material employed in their experiments. We can however exclude these factors in our experiments, where the interface grain-wall is small with respect to the bulk and the beads are made of glass. The leftwards asymmetry observed in experiments represents a very interesting phenomenon, which in general is expected from non trivial dynamical effects, and cannot be due to the simple inertia of the moving plate (which should produce opposite asymmetry Baldassarri et al. (2003); Colaiori et al. (2004)).

In some magnetic materials, a leftwards asymmetry has been observed and related to memory effects acting as an effective negative mass of domain walls Zapperi et al. (2005). In our experiments we cannot exclude the existence of such an ”effective” inertia of the system, due to some memory introduced by the underlying granular. For instance, in some experiments Nasuno et al. (1998) researchers noted an increased inertia of the slider moving on a granular bed, due to the grains dragged by the slider itself and a similar augmented inertia has been observed also in our previous experiments Baldassarri et al. (2006). Since the quantity of grains dragged by the disk during its motion could change during the irregular motion of the system one should consider the inertia as a dynamical quantity, rather than a constant, and this could in principle be one origin of the asymmetries. A left asymmetry has been also observed in earthquakes Houston et al. (1998); Mehta et al. (2006).

### .6 Breaking of scaling

More insight into the mechanisms leading to the scaling breakdown can be gained by looking again at the plot relating and shown in Fig. 3. The first information coming from this plot is that there exists a definite statistical scaling between slip size and duration, as shown by the scattering of data. The white squared symbols in the main plot represent the average slip size and duration of each class (statistical errors are negligible on these averages). It is seen that, at least for the four lower classes, they follow an algebraic relation: (red continuous line). The value of the exponent turns out to be . This is close to, but clearly different from, the value of expected from extant models (for instance the ABBM model mentioned above Alessandro et al. (1990)).

The other information supplied by the scatter plot of Fig. 3 is that this behavior changes at large slips, where a linear dependence, looks more appropriate (yellow dashed line). Interestingly, the crossover between the two behaviors takes place around the fourth class, exactly where the scaling of the average pulse shape, shown in Fig. 4, breaks down. The inset of Fig. 3 puts into better evidence this cross-over. There, we have plotted the quantity as function of . We observe a weakly superlinear relation for small slips, followed by a plateau at large slips. Note that the ratio between and is nothing but the plate average velocity during the slip. This observation allows to identify a critical velocity, as the ratio between the average slip size, rad, and the average duration, s, of the fourth class: rad/s. We speculate that during large slips (), when the plate reaches high velocities , it could experience some sudden increase of friction. In the next section it will be seen that this increase indeed appears, as a dynamical effect.

### .7 Stochastic friction

The forces ruling the slip dynamics are the spring torque and the granular friction. While the first one just depends linearly on the instantaneous angle, the second displays a complex behavior (see Fig. 1, central and right panel) from which interesting features emerge.

The classical Mohr–Coulomb criterion predicts constant friction at low shear, and increasing values when the system enters the Bagnold’s regime Bagnold (1954, 1962), a behavior well observed experimentally at constant shear (see e.g. Savage and Sayed (1984)).

However, it is doubtful whether this behavior could be relevant to the stick-slip dynamics observed in the critical regime. More generally, friction in granular systems is usually measured under controlled shear strain or stress, but its properties can be dramatically different when observed in the self-organized state, as exemplified in Fig. 1, right panel. Some statistical features of friction in this state have been investigated in Dalton et al. (2005); Baldassarri et al. (2006); Petri et al. (2008), but despite this quantity plays a crucial role for the system dynamics, it has never been systematically measured to date during stick slip.

In the critical regime friction is a random quantity. It depends on the details of the network of contacts between particles in the granular bed. Fluctuations in the frictional response of the granular medium result from the stress propagation on the evolving network of grain contacts, and are at the very origin of the motion stochasticity. This fact has a number of consequences and some subtleties. A random friction force, as a stochastic quantity, can be described by statistical estimators like averages, moments, correlators, etc. Nevertheless, several averages can be defined, which depend on the driving protocol and can be very different from each other. More specifically, one can consider the time average of the friction over the full dynamics, but this is not always really meaningful, especially in the critical regime. As shown in Dalton et al. (2005) the statistical distribution of friction in this regime is characterized by fat tails, as opposite with the continuous sliding where it is normal. Another possible average, Baldassarri et al. (2006); Leoni et al. (2010) is the average friction conditioned to the (instantaneous) plate velocity:

(5) |

In Fig. 5 we plot such conditioned average friction during the stick slip critical regime. As noted in Baldassarri et al. (2006), an interesting Stribeck-shaped (that is, a shear weakening followed by a thickening) friction curve appears, featuring weakening for small velocities and recovering the Bagnold behavior at high velocities. However, this velocity weakening arises as a dynamical effect. In fact, a different driving protocol can give different results: For instance, at constant shear Savage and Sayed (1984) the average friction is constant at low and intermediate speeds.

The analysis of Fig. 5 allows to identify a velocity corresponding to the position of the minimum of the average friction . Our experiments clearly indicate that the position of this minimum does not depends on the drive velocity (Fig. 13 in Appendix) and it is always attained near the velocity rad/s. This value is very close to the value marking the crossover in the scaling of vs (see Fig. 3), which in turn is related to the breakdown of scaling of the average avalanche shape shown in Fig. 4. This corroborates the previous interpretation of the crossover phenomena and the breaking of the critical scaling of the dynamics as due to the weakening followed by the increase of friction experienced by the plate during larger, faster avalanches Fig. 5.

In order to better investigate whether and how friction dynamical behavior can influence the average velocity shape we have also analyzed the average shape of friction along the slip. In an analogous fashion to what done for computing the velocity shapes, one can define as the average frictional torque for slips of the same duration . In practice, we have computed the average value of the friction torque over slips of similar duration , according to the same classes of duration adopted for velocities (see Fig. 3 and Appendix). The results, presented in Fig. 6, show that the breaking of scaling of the velocity shapes corresponds to a change in the frictional properties of . For small avalanches, i.e. those corresponding to the cases in which average velocity shape obeys scaling (curves plotted in the left graph of Fig. 6), the average friction maintains an almost constant value along the whole slip, whose value is independent fron the slip duration. On the contrary, the curves corresponding to longer slips (shown in the right t plot of Fig. 6) display different shapes that, as in the case of velocity (Fig. 4), strongly depend on , and cannot be collapsed. Note also that higher frictions are experienced during longer slips.

Let us stress here the difference between the two average procedures considered in this work. The average shown in Fig. 6, are performed over slips of similar duration, at the same internal avalanche time . Instead, the (conditional) average , defined in Eq. 5 and shown in Fig. 5, mixes events of any duration, and it depends on the instantaneous plate velocity . The two quantities give different aspects of the same (stochastic) physical phenomenon. Nevertheless, the combination of the two analysis suggests that the quite complex friction weakening behavior of is mainly due to large slips, which show a non constant average friction in time (see Fig. 6, right panel), in contrast with small avalanches, where the average keeps mainly constant.

By combining the analysis of friction and velocity shapes, one can consider the curves resulting from plotting the average friction as a function of the average velocity, in slips of similar duration, as shown in Fig. 7. They show that while, as anticipated, friction has a low velocity dependence in small slips (left panel), in large ones it splits into a two-valued function (right panel), displaying an hysteresis, with well different dependencies on the (average) plate velocity in the accelerating and decelerating phases of the slips. This evolution is very similar to that observed in periodic stick-slip Nasuno et al. (1997, 1998), where all slips have identical extension, duration, and velocity profile.

### .8 Considerations and conclusions

Our experiments show a good scaling of the average velocity shape for small avalanches, with an almost symmetric average shape. For larger avalanches however, scaling (Eq. 2) is broken: for large slips the shape takes a clear leftwards shape in agreement with what observed in seismic data Houston et al. (1998); Barés et al. (2017) (and recently in Barés et al. (2014)).

Our analyses show that the breakdown of velocity scaling goes along with changes in the friction behavior, pointing out a strict relation between the two phenomena. On the opposite, spring-block models with only Coulomb friction generate symmetric slips Aharonov and Spark (2004); Bizzarri and Petri (2016). Effective friction laws accounting for elapsed time and/or space have been incorporated in solid-on-solid interface models, through the dependence on so called state variables Ruina (1983); Dieterich (1994); Baumberger and Caroli (2006). These rate-and-state laws are often adopted for studying and modeling co-seismic fault shearing, together with their other simpler forms Scholz (1998); Bizzarri and Cocco (2003); Kawamura et al. (2012); Bizzarri and Petri (2016). They are essentially phenomenological and can describe both velocity weakening and hardening, depending on the adopted parameters (which are not derivable from microscopic principles). These laws have shown to work to some extent also for interstitial granular matter Marone (1998), but with some inconsistencies Mair and Marone (1999). Moreover they have been drawn from experiments where velocity is forced to change in sudden steps and they don’t seem to have been never investigated in the critical stick-slip. Attempts to do this with smoothly varying velocity have been done in Leoni et al. (2010). In a very recent work DeGiuli and Wyart (2017) both friction weakening and hysteresis have been numerically investigated during granular shear cycles, showing that these are due to contact instabilities induced by the acoustic waves generated during granular fluidization. It is thus clear that granular flow cannot be effectually modeled without the inclusion of more refined and realistic friction laws.

An effective modeling approach to the critical granular dynamics cannot as well exclude a stochastic description of friction, which generates the slip unpredictability and their range of variability, with the following change in the slip shapes. To our knowledge, the only few attempts in this direction Baldassarri et al. (2006); Leoni et al. (2010); Dahmen et al. (2011) are inspired to the aforementioned ABBM model Alessandro et al. (1990), which represents the mean field approximation for the motion of a driven elastic interface in a random environment Zapperi et al. (1998); Le Doussal and Wiese (2012). From the dynamical point of view, it describes a spring–slider model subjected to a friction where both viscous and a random pinning components are present, in the overdamped (i.e. negligible inertia) approximation. At small, but finite driving rate, the ABBM model predicts an intermittent, critical dynamics for the block motion. Similarly to our observations, avalanche statistics show a scaling regime for short slips, whose average velocity has parabolic shapes. However, an exponent relates to , which is different from what observed in our experiment. Moreover, for longer slips, ABBM predicts flatter symmetical shapes, witnessing a cut-off in the velocity correlation. No inertial effects are present, due to the overdamped approximation.

A variant of the ABBM model for critical granular dynamics has been introduced in Baldassarri et al. (2006), where, based on empirical observations, a simple Stribeck-like rate dependence, showing a minimum, of the average granular friction was adopted. Moreover, more physically, a cut off in the spatial correlation of the random force was considered and, at odds with the original model, inertia was taken in account. Later on, attempts to introduce in the model a state dependent weakening friction have been done Leoni et al. (2010)), and further investigations are in progress.

We think that the insights provided by the present study can explain the contradictory recent observations in Barés et al. (2017); Denisov et al. (2017) and can be useful to advance such efforts to improve models. In particular, they show that inertia can play an important role in both weakening-hysteresis DeGiuli and Wyart (2017) and in the determining scaling exponent (an inertial ABBM model has been studied in Le Doussal et al. (2012)). Even at the microscopic level, grain inertia can influence the avalanche statistics. For instance, in sandpile models, largely studied in the context of SOC (Self Organized Criticality), the tendency of real sand grains to keep moving once they start facilitate the emergence of huge avalanches. Recent theoretical developments propose, in the presence of such facilitation effects, a scenario called Self-Organised Bistability Di Santo et al. (2016), where again a breaking of scaling is associated to the appearance of large avalanches (”kings”).

We conclude that weakening is a genuine property exhibited by granular dynamics at variable shear rate, and that randomness and memory are a general features of friction that cannot be overlooked in the formulation of effectual models. Such models can have impact on the understanding of many phenomena occurring in the large realm of granular systems, and in particular of self organized natural phenomena like landslides, and earthquake, where it is not yet clear the way different mechanisms can contribute to the shear weakening observed in coseismic fault shearing Di Toro et al. (2011). Further investigation on theoretical models incorporating more realistic, specifically memory dependent friction laws, and new experiments will allow to better understand the mechanism for which criticality breaks down.

## Acnokwlegments

This work has been supported by the grant FIRB RBFR081IUK_003 from the Italian Ministry for Education and Research

## I Appendix

### i.1 The experimental set up

The experimental apparatus utilized for this research consists of a circular PPMI channel of outer and inner radii = 19.2 cm and = 12.5 cm respectively. The channel is 12 cm height and is almost filled with a bidisperse mixture 50%-50% of glass beads, with radii =1.5 mm 10% and =2 mm 10%.

A top plate, fitting the channel, can be rotated and has a few layer of grains glued to its lower face in order to better drag the underlying granular medium. The plate has mass = 1200 g and moment of inertia = 0.026 kg m, and it is free to move vertically, implying that in our experiments the medium can change volume under a nominal pressure of 176 Pa. The plate is connected to an end of torsion spring, of elastic constant 0.36 Nm/rad, while the other end of the spring is rotated by a motor at constant angular velocity . The angular positions of motor and plate are supplied by two optical encoders positioned on either side of the torsion spring, each one having a spatial resolution of rad and being sampled at 50 Hz. These measures provides the plate instantaneous position and velocity, and , as well as the friction torque, which is proportional to the angular difference between motor and plate(see Eq. (1)).

### i.2 Experimental analysis

In principle each single slip event, or avalanche, begins when starts to differ from zero and ends when goes back to zero. However, in practice it is necessary to choose a threshold value to cross, in order to get rid of the instrumental noise. This choice is to some extent arbitrary, however all the results have been observed to be independent from the chosen threshold, as long as it is small and different from 0. For our analysis we have set rad/s, and considered the seven time series reported in table 1.

series | duration | # of points | driving | # of slips |
---|---|---|---|---|

(minutes) | (rad/s) | used in analysis | ||

3900 | 5849962 | 0.0015 | 6014 | |

673 | 2020079 | 0.0022 | 1625 | |

1200 | 3600060 | 0.0044 | 5826 | |

4080 | 12240020 | 0.0055 | 2451 | |

360 | 1079977 | 0.011 | 3725 | |

240 | 720007 | 0.021 | 3973 | |

210 | 630014 | 0.033 | 4300 |

The avalanches of each series have been grouped into classes on the base of their duration, according to the first column of table 2. For each class the average avalanche duration and size have been evaluated, and instantaneous average velocity has been computed at a set of discrete times, (see main text).

Avalanches at the extremes of distributions have been dropped out. For example duration and size for avalanches from the series (EA) are plotted in Fig. 2 of the main text, with different colors for each interval. Avalanches in gray, shorter than 0.31 s, are too small to perform meaningful analysis (less than 15 points at Hz of sampling rate). The total number of avalanches employed in this statistics has then been 6014, distributed according to the second column of table 2.

The main text presents results from the series . The results from the other datasets, with the different drives reported in Table 1, display similar behaviors and are shown in Figs. 9-14 of this Appendix, to be compared with the correspondig Figs. 2-7 in the main text. Analogous results were obtained adopting different sampling frequencies and threshold values.

duration | # of avalanches |
---|---|

0.309 | 929 |

0.489 | 866 |

0.772 | 987 |

1.219 | 1694 |

1.925 | 1380 |

3.04 | 158 |

## References

- Jaeger et al. (2000) H. M. Jaeger, T. Shinbrot, and P. B. Umbanhowar, Proceedings of the National Academy of Sciences 97, 12959 (2000).
- ann (2016) Journal of Physics: Condensed Matter 28, 134001 (2016).
- Dalton and Corcoran (2001) F. Dalton and D. Corcoran, Physical Review E 63, 061312 (2001).
- Ciamarra et al. (2010) M. P. Ciamarra, E. Lippiello, C. Godano, and L. de Arcangelis, Physical review letters 104, 238001 (2010).
- Dalton et al. (2005) F. Dalton, F. Farrelly, A. Petri, L. Pietronero, L. Pitolli, and G. Pontuale, Phys. Rev. Lett. 95, 138001 (2005).
- Baldassarri et al. (2006) A. Baldassarri, F. Dalton, A. Petri, S. Zapperi, G. Pontuale, and L. Pietronero, Physical Review Letters 96, 118002 (2006).
- Petri et al. (2008) A. Petri, A. Baldassarri, F. Dalton, G. Pontuale, L. Pietronero, and S. Zapperi, The European Physical Journal B - Condensed Matter and Complex Systems 64, 531 (2008).
- Geller et al. (2015) D. A. Geller, R. E. Ecke, K. A. Dahmen, and S. Backhaus, Phys. Rev. E 92, 060201 (2015).
- Zadeh et al. (2018) A. A. Zadeh, J. Barés, J. E. S. Socolar, and R. P. Behringer, preprint (2018), arXiv:1810.12243v1 [cond-mat.soft] .
- Sethna et al. (2001) J. P. Sethna, K. A. Dahmen, and C. R. Myers, Nature 410, 242 (2001).
- Main (1996) I. Main, Reviews of Geophysics , 433 (1996).
- Petri et al. (1994) A. Petri, G. Paparo, A. Vespignani, A. Alippi, and M. Costantini, Phys. Rev. Lett. 73, 3423 (1994).
- Cannelli et al. (1993) G. Cannelli, R. Cantelli, and F. Cordero, Phys. Rev. Lett. 70, 3923 (1993).
- Dimiduk et al. (2006) D. M. Dimiduk, C. Woodward, R. LeSar, and M. D. Uchic, Science 312, 1188 (2006).
- Bak et al. (1988) P. Bak, C. Tang, and K. Wiesenfeld, Phys. Rev. A 38, 364 (1988).
- Dauchot et al. (2005) O. Dauchot, G. Marty, and G. Biroli, Phys. Rev. Lett. 95, 265701 (2005).
- Kuntz and Sethna (2000) M. C. Kuntz and J. P. Sethna, Physical Review B - Condensed Matter and Materials Physics 62, 11699 (2000), arXiv:9911207 [cond-mat] .
- Baldassarri et al. (2003) A. Baldassarri, F. Colaiori, and C. Castellano, Physical Review Letters 90, 60601 (2003).
- Colaiori et al. (2004) F. Colaiori, A. Baldassarri, and C. Castellano, Physical Review E (Statistical, Nonlinear, and Soft Matter Physics) 69, 41105 (2004).
- Colaiori (2008) F. Colaiori, Advances in Physics 57, 287 (2008).
- Zapperi et al. (2005) S. Zapperi, C. Castellano, F. Colaiori, and G. Durin, Nature Physics 1, 46 (2005).
- Papanikolaou et al. (2011) S. Papanikolaou, F. Bohn, R. L. Sommer, G. Durin, S. Zapperi, and J. P. Sethna, Nature Physics 7, 316 (2011).
- Chrzan and Mills (1994) D. C. Chrzan and M. J. Mills, Physical Review B 50, 30 (1994).
- Sparks et al. (2017) G. Sparks, J. Sickle, K. Dahmen, and R. Maas, arXiv:1705.06636 (2017).
- Antonaglia et al. (2014) J. Antonaglia, W. J. Wright, X. Gu, R. R. Byer, T. C. Hufnagel, M. LeBlanc, J. T. Uhl, and K. A. Dahmen, Phys. Rev. Lett. 112, 155501 (2014).
- Liu et al. (2016) C. Liu, E. E. Ferrero, F. Puosi, J.-L. Barrat, and K. Martens, Phys. Rev. Lett. 116, 065501 (2016).
- Lagogianni et al. (2018) A. E. Lagogianni, C. Liu, K. Martens, and K. Samwer, The European Physical Journal B 91, 104 (2018).
- Danku and Kun (2013) Z. Danku and F. Kun, Phys. Rev. Lett. 111, 084302 (2013).
- Roberts et al. (2014) J. A. Roberts, K. K. Iyer, S. Finnigan, S. Vanhatalo, and M. Breakspear, The Journal of Neuroscience 34, 6557 LP (2014).
- Wikstro et al. (2015) S. Wikstro, K. K. Iyer, J. A. Roberts, L. Hellstro, I. H. Pupp, D. Ley, S. Vanhatalo, and M. Breakspear, Brain , 2206 (2015).
- Wang et al. (2013) B. Wang, J. Kuo, and S. Granick, Phys. Rev. Lett. 111, 208102 (2013).
- Gallotti and Chialvo (2018) R. Gallotti and D. R. Chialvo, Journal of The Royal Society Interface 15 (2018).
- Chialvo et al. (2015) D. Chialvo, A. G. Torrado, E. Gudowska-Nowak, J. Ochab, P. Montoya, M. Nowak, and E. Tagliazucchi, Papers in Physics 7 (2015).
- Sheikh et al. (2016) M. A. Sheikh, R. L. Weaver, and K. A. Dahmen, Phys. Rev. Lett. 117, 261101 (2016).
- Consolini et al. (2008) G. Consolini, P. De Michelis, and R. Tozzi, Journal of Geophysical Research: Space Physics 113, 1 (2008).
- Mehta et al. (2006) A. P. Mehta, K. A. Dahmen, and Y. Ben-Zion, Phys. Rev. E 73, 056104 (2006).
- Laurson et al. (2013) L. Laurson, X. Illa, S. Santucci, K. Tore Tallakstad, K. J. Måloy, and M. J. Alava, Nat. Commun. 4 (2013).
- Barés et al. (2017) J. Barés, D. Wang, D. Wang, T. Bertrand, C. S. O’Hern, and R. P. Behringer, Phys. Rev. E 96, 052902 (2017), 1709.01012 .
- Denisov et al. (2017) D. V. Denisov, K. A. Lörincz, W. J. Wright, T. C. Hufnagel, A. Nawano, X. Gu, J. T. Uhl, K. A. Dahmen, and P. Schall, Scientific Reports 7, 43376 (2017).
- Leoni et al. (2010) F. Leoni, A. Baldassarri, F. Dalton, A. Petri, G. Pontuale, and S. Zapperi, Journal of Non-Crystalline Solids 357, 749 (2010).
- Ciamarra et al. (2012) M. P. Ciamarra, F. Dalton, L. de Arcangelis, C. Godano, E. Lippiello, and A. Petri, Tribology Letters 48, 89 (2012).
- Alessandro et al. (1990) B. Alessandro, C. Beatrice, G. Bertotti, and A. Montorsi, Journal of Applied Physics 68, 2901 (1990).
- Feller (1951) W. Feller, Annals of mathematics 54, 173 (1951).
- Dobrinevski et al. (2012) A. Dobrinevski, P. Le Doussal, and K. J. Wiese, Physical Review E - Statistical, Nonlinear, and Soft Matter Physics 85, 1 (2012), 1112.6307 .
- Nasuno et al. (1998) S. Nasuno, A. Kudrolli, A. Bak, and J. P. Gollub, Phys. Rev. E 58, 2161 (1998).
- Houston et al. (1998) H. Houston, H. M. Benz, and J. E. Vidale, Journal of Geophysical Research 103, 29895 (1998).
- Bagnold (1954) R. A. Bagnold, Proc. Roy. Soc. A 225, 49 (1954).
- Bagnold (1962) R. A. Bagnold, Proc. Roy. Soc. A 295, 219 (1962).
- Savage and Sayed (1984) S. B. Savage and M. Sayed, J. Fluid Mech. 142, 391 (1984).
- Nasuno et al. (1997) S. Nasuno, A. Kudrolli, and J. P. Gollub, Phys. Rev. Lett. 79, 949 (1997).
- Barés et al. (2014) J. Barés, M. L. Hattali, D. Dalmas, and D. Bonamy, Phys. Rev. Lett. 113, 264301 (2014).
- Aharonov and Spark (2004) E. Aharonov and D. Spark, Journal of Geophysical Research: Solid Earth 109 (2004).
- Bizzarri and Petri (2016) A. Bizzarri and A. Petri, Physics of the Earth and Planetary Interiors 261, 36 (2016).
- Ruina (1983) A. Ruina, Journal of Geophysical Research: Solid Earth 88, 10359 (1983).
- Dieterich (1994) J. Dieterich, Journal of Geophysical Research: Solid Earth 99, 2601 (1994).
- Baumberger and Caroli (2006) T. Baumberger and C. Caroli, Advances in Physics 55, 279 (2006).
- Scholz (1998) C. H. Scholz, Nature 391, 37 (1998).
- Bizzarri and Cocco (2003) A. Bizzarri and M. Cocco, Journal of Geophysical Research: Solid Earth 108, 2373 (2003).
- Kawamura et al. (2012) H. Kawamura, T. Hatano, N. Kato, S. Biswas, and B. K. Chakrabarti, Rev. Mod. Phys. 84, 839 (2012).
- Marone (1998) C. Marone, Ann. Revs. Earth & Plan. Sci. 26, 643 (1998).
- Mair and Marone (1999) K. Mair and C. Marone, Journal Geophysical Research 104, 899 (1999).
- DeGiuli and Wyart (2017) E. DeGiuli and M. Wyart, Proceedings of the National Academy of Sciences 114, 9284 (2017), http://www.pnas.org/content/114/35/9284.full.pdf .
- Dahmen et al. (2011) K. A. Dahmen, Y. Ben-Zion, and J. Uhl, Nature Phyisics , 1957 (2011).
- Zapperi et al. (1998) S. Zapperi, P. Cizeau, G. Durin, and H. E. Stanley, Physical Review B 58, 6353 (1998), 9803253 [cond-mat] .
- Le Doussal and Wiese (2012) P. Le Doussal and K. J. Wiese, Europhys. Lett. 97, 46004 (2012), 1104.2629 .
- Le Doussal et al. (2012) P. Le Doussal, A. Petković, and K. J. Wiese, Physical Review E 85, 061116 (2012), 1203.5620 .
- Di Santo et al. (2016) S. Di Santo, R. Burioni, A. Vezzani, and M. A. Muñoz, Physical Review Letters 116, 1 (2016), arXiv:1605.05161 .
- Di Toro et al. (2011) G. Di Toro, R. Han, T. Hirose, N. De Paola, S. Nielsen, K. Mizoguchi, F. Ferri, M. Cocco, and T. Shimamoto, Nature 471, 7339 (2011).