Fault zone complexity naturally produces the full slip spectrum: Insights from numerical models

In addition to regular earthquakes, observations of spatiotemporally complex slip events have multiplied over the last decades. These slip events range along different time scales: from creep , slow slip events to LFEs and tremors. At present, these events are generally interpreted by imposed frictional heterogeneities along the fault plane. However, fault systems are geometrically complex in nature over different scales. We aim in this work to investigate the role of “realistic” fault geometry on the dynamics of slip events. We consider a fault system in a 2D quasi-dynamic setting. The fault system consists of a main self-similar rough fault, surrounded by a dense network of off-fault fractures. All fractures are frictionally homogeneous (rate weakening) and can potentially undergo dynamic slip. We aim to understand how the deformation in the volume is accomodated by the off-fault damage zone and the main fault. What fraction of the “supplied” moment rate is hosted by the off-fault fractures during an earthquake cycle?