Preavalanche instabilities in a granular pile

We investigate numerically the transition between static equilibrium and dynamic surface flow of a 2D cohesionless granular system driven by a continuous gravity loading. This transition is characterized by intermittent local dynamic rearrangements and can be described by an order parameter defined as the density of critical contacts, i.e., contacts where the friction is fully mobilized. Analysis of the spatial correlations of critical contacts shows the occurrence of "fluidized" clusters which exhibit a power-law divergence in size at the approach of the stability limit. The results are compatible with recent models that describe the granular system during the static/dynamic transition as a multiphase system.