Person:

Platt, John

Field observations of major earthquake fault zones show that shear deformation is often confined to principal slipping zones that may be of order 1–100 μm wide, located within a broader gouge layer of order 10–100 mm wide. This paper examines the possibility that the extreme strain localization observed may be due to the coupling of shear heating, thermal pressurization, and diffusion. In the absence of a stabilizing mechanism shear deformation in a continuum analysis will collapse to an infinitesimally thin zone. Two possible stabilizing mechanisms, studied in this paper, are rate-strengthening friction and dilatancy. For rate-strengthening friction alone, a linear stability analysis shows that uniform shear of a gouge layer is unstable for perturbations exceeding a critical wavelength. Using this critical wavelength we predict a width for the localized zone as a function of the gouge properties. Taking representative parameters for fault gouge at typical centroidal depths of crustal seismogenic zones, we predict localized zones of order 5–40 μm wide, roughly consistent with field and experimental observations. For dilatancy alone, linearized strain rate perturbations with a sufficiently large wavelength will undergo transient exponential growth before decaying back to uniform shear. The total perturbation strain accumulated during this transient strain rate localization is shown to be largely controlled by a single dimensionless parameter E, which is a measure of the dilatancy of the gouge material due to an increase in strain rate.

Field and laboratory observations indicate that at seismic slip rates most shearing is confined to a very narrow zone, just a few tens to hundreds of microns wide, and sometimes as small as a few microns. Rice et al. (2014) analyzed the stability of uniform shear in a fluid-saturated gouge material. They considered two distinct mechanisms to limit localization to a finite thickness zone, rate-strengthening friction, and dilatancy. In this paper we use numerical simulations to extend beyond the linearized perturbation context in Rice et al. (2014), and study the behavior after the loss of stability. Neglecting dilatancy we find that straining localizes to a width that is almost independent of the gouge layer width, suggesting that the localized zone width is set by the physical properties of the gouge material. Choosing parameters thought to be representative of a crustal depth of 7 km, this predicts that deformation should be confined to a zone between 4 and 44 μm wide. Next, considering dilatancy alone we again find a localized zone thickness that is independent of gouge layer thickness. For dilatancy alone we predict localized zone thicknesses between 1 and 2 μm wide for a depth of 7 km. Finally, we study the impact of localization on the shear strength and temperature evolution of the gouge material. Strain rate localization focuses frictional heating into a narrower zone, leading to a much faster temperature rise than that predicted when localization is not accounted for. Since the dynamic weakening mechanism considered here is thermally driven, this leads to accelerated dynamic weakening.