Modeling dynamic ruptures on extended faults for microearthquakes induced by fluid injection

Understanding the dynamics of microearthquakes is a timely challenge to solve current paradoxes in earthquake mechanics, such as the stress drop and fracture energy scaling with seismic moment. Dynamic modeling of microearthquakes induced by fluid injection is also relevant for studying rupture propagation following a stimulated nucleation. We study the main features of unstable dynamic ruptures caused by fluid injection on a target preexisting fault (50m x 50m) generating a M_w=1 event. The selected fault is located in the Bedretto Underground Laboratory (Swiss Alps) at ≈1000m depth. We perform fully dynamic rupture simulations and model seismic wave propagation in 3D by adopting a linear slip-weakening law. We use the distributed multi-GPU implementation of SeisSol on the supercomputer Marconi100.

 Stress field and fault geometry are well constrained by in-situ observations, allowing us to minimize the a priori imposed parameters. We investigate the scaling relations of stress drop, slip-weakening distance (D_c) and fracture energy (G_c) focusing on their role in governing dynamics of rupture propagation and arrest for a target M_w=1 induced earthquake. We explore different homogenous conditions of frictional parameters, and we show that the spontaneous arrest of the rupture is possible in the modeled stress regime, by assuming a high ratio between stress excess and dynamic stress drop (the fault strength parameter S), characterizing the fault before the fluid pressure change. The rupture arrest of modeled induced earthquakes depends on the heterogeneity of dynamic parameters due to the spatially variable effective normal stress. Moreover, for a fault with high S values (not ready to slip), small variations of D_c (≈0.5÷1.2 mm) can drive the rupture from self-arrested to run-away. Studying dynamic interactions (stress transfer) among slipping points on the rupturing fault provides insight on the breakdown process zone and shear stress evolution at the crack tip leading to failure. The inferred spatial dimension of the cohesive zone in our models is nearly ~0.3-0.4m, with a maximum slip of ~0.6 cm. Finally, we compare stress drop and fracture energy estimated from synthetic waveforms with assumed dynamic parameters. Our results suggest that meso-scale processes near the crack-tip affect rupture dynamics of micro-earthquakes.