Breaking Up: Comminution Mechanisms in Sheared Simulated Fault Gouge

The microstructural state and evolution of fault gouge has important implications for the mechanical behaviour, and hence the seismic slip potential of faults. We use 3D discrete element (DEM) simulations to investigate the fragmentation processes operating in fault gouge during shear. Our granular fault gouge models consist of aggregate grains, each composed of several thousand spherical particles stuck together with breakable elastic bonds. The aggregate grains are confined between two blocks of solid material and sheared under a given normal stress. During shear, the grains can fragment in a somewhat realistic way leading to an evolution of grain size, grain shape and overall texture. The ‘breaking up’ of the fault gouge is driven by two distinct comminution mechanisms: grain abrasion and grain splitting. The relative importance of the two mechanisms depends on applied normal stress, boundary wall roughness and accumulated shear strain. If normal stress is sufficiently high, grain splitting contributes significantly to comminution, particularly in the initial stages of the simulations. In contrast, grain abrasion is the dominant mechanism operating in simulations carried out at lower normal stress and is also the main fragmentation mechanism during the later stages of all simulations. Rough boundaries promote relatively more grain splitting whereas smooth boundaries favor grain abrasion. Grain splitting (plus accompanying abrasion) appears to be an efficient mechanism for reducing the mean grain size of the gouge debris and leads rapidly to a power law size distribution with an exponent that increases with strain. Grain abrasion (acting alone) is an effective way to generate excess fine grains and leads to a bimodal distribution of grain sizes. We suggest that these two distinct mechanisms would operate at different stages of a fault’s history. The resulting distributions in grain size and grain shape may significantly affect frictional strength and stability. Our results therefore have implications for the earthquake potential of seismically active faults with accumulations of gouge. They may also be relevant to the susceptibility of rockslides since non-cohesive basal shear zones will evolve in a similar way and potentially control the dynamics of the slide.