How cementation and fluid flow influence slip behavior at the subduction interface

Much of the complexity of subduction-zone earthquake size and temporal patterns owes to linkages among fluid flow, stress, and fault healing. To investigate these linkages, we introduce a novel numerical model that tracks cementation and fluid flow within the framework of an earthquake simulator. In the model, there are interseismic increases in cohesion across the plate boundary and decreases in porosity and permeability caused by cementation along the interface. Seismogenic slip is sensitive to the effective stress and therefore fluid pressure; in turn, slip events increase porosity by fracturing. The model therefore accounts for positive and negative feedbacks that modify slip behavior through the seismic cycle. The model produces temporal clustering of earthquakes in the seismic record of the Aleutian margin, which has well-documented along-strike variations in locking characteristics. Model results illustrate how physical, geochemical, and hydraulic linkages can affect natural slip behavior.