A Discrete Elements Study of the Frictional Behavior of Fault Gouges

A series of discrete elements simulations is presented for the study of fault gouges' frictional response. The gouge is considered to have previously undergone ultra‐cataclastic flow and long‐time consolidation loading. We explore the effect of different particle characteristics such as size, polydispersity, and also shearing velocities on gouge's response under the conditions met in the seismogenic zone. Monte‐Carlo analyses suggest that the local stick‐slip events disappear when averaging over a large number of numerical samples. Moreover, the apparent material frictional response remains almost unaffected by the spatial randomness of particles' position and by the particle's size distribution. On the contrary, the mean particle size controls the formation and thickness of the observed shear bands, which appear after the peak friction is met. Furthermore, the apparent friction evolution fits well to an exponential decay law with slip, which involves a particle size dependent critical slip distance. For the studied conditions and depth, the shearing velocity is found to play a secondary role on the apparent frictional response of the gouge, which highlights the importance of analyses involving multiphysics for studying the rheology of fault gouges. Besides improving the understanding of the underlying physics of the problem, the above findings are also useful for deriving pertinent constitutive models in the case of modeling with continuum theories.