Model constraints on infiltration of surface-derived fluids deep into the brittle crust

Abstract

Summary Downward flow of surface-derived water deep into the upper crust is investigated using two dimensional coupled hydrothermal numerical models. In the models, downward flow is driven by either topographic gradients or seismic pumping, while it is facilitated by large episodic variations in fault permeability, intended to mimic fracturing and healing on a fault over repeated seismic cycles. The models show that both forcing scenarios are equally capable of driving surface-derived fluid to the base of faults at 10 km depth in several tens of thousands of years under certain conditions. Downward flow of cold fluid occurs almost exclusively during and shortly after earthquakes, while during the remaining portion of the seismic cycle fluids remain relatively stationary while they undergo thermal relaxation (i.e., heating). Rapid downward flow is favoured by a large coseismic permeability, long permeability healing time scale, and large coseismic dilatancy or high topographic relief above the fault at the surface. However, downward fluid flow is completely inhibited if fluid pressures exceeds the hydrostatic gradient, even by modest amounts, which suggests that deep fluid infiltration is unlikely to occur in every region.