Thermal pressurization and onset of melting in fault zones

Abstract

We examine how frictional heating drives the evolution of temperature, strength, and fracture energy during earthquake slip. For small slip distances, heat and pore fluid are unable to escape the shearing fault core, and the behavior is well approximated by simple analytical models that neglect any transport. Following large slip distances, the finite width of the shear zone is small compared to the thicknesses of the thermal and hydrological boundary layers, and the fault behavior approaches that predicted for the idealized case of slip on a plane. To evaluate the range in which the predictions of these two sets of approximations are valid, we develop a model that describes how frictional dissipation within a finite shear zone drives heat and mass transport through the surrounding static gouge. With realistic parameter values and slips greater than a few centimeters, the subsequent evolution of strength and fracture energy are approximated well by the planar slip model. However, the temperature evolution is much more sensitive to the finite shear zone thickness, and the ultimate temperature rise tends to be intermediate between that predicted for the two simplified cases. We explore the range of conditions necessary for melting to begin and focus in particular on the potential role of fault zone damage in facilitating fluid transport and promoting larger temperature increases. We discuss how the apparent scarcity of exhumed pseudotachylytes places constraints on some of the more uncertain fault zone parameters.