Seismicity in a model governed by competing frictional weakening and healing mechanisms

Observations from laboratory, field and numerical work spanning a wide range of space and time scales suggest a strain dependent progressive evolution of material properties that control the stability of earthquake faults. The associated weakening mechanisms are counterbalanced by a variety of restrengthening mechanisms. The efficiency of the healing processes depends on local material properties and on rheologic, temperature, and hydraulic conditions. We investigate the relative effects of these competing non-linear feedbacks on seismogenesis in the context of evolving frictional properties, using a mechanical earthquake model that is governed by slip weakening friction. Weakening and strengthening mechanisms are parametrized by the evolution of the frictional control variable—the slip weakening rate R—using empirical relationships obtained from laboratory experiments. In our model, weakening depends on the slip of an earthquake and tends to increase R, following the behaviour of real and simulated frictional interfaces. Healing causes R to decrease and depends on the time passed since the last slip. Results from models with these competing feedbacks are compared with simulations using non-evolving friction. Compared to fixed R conditions, evolving properties result in a significantly increased variability in the system dynamics. We find that for a given set of weakening parameters the resulting seismicity patterns are sensitive to details of the restrengthening process, such as the healing rate b and a lower cutoff time, tc , up to which no significant change in the friction parameter is observed. For relatively large and small cutoff times, the statistics are typical of fixed large and small R values, respectively. However, a wide range of intermediate values leads to significant fluctuations in the internal energy levels. The frequency-size statistics of earthquake occurrence show corresponding non-stationary characteristics on time scales over which negligible fluctuations are observed in the fixed-R case. The progressive evolution implies that -— except for extreme weakening and healing rates -— faults and fault networks possibly are not well characterized by steady states on typical catalogue time scales, thus highlighting the essential role of memory and history dependence in seismogenesis. The results suggest that an extrapolation to future seismicity occurrence based on temporally limited data may be misleading due to variability in seismicity patterns associated with competing mechanisms that affect fault stability