Abstract Modeling of injection‐induced seismicity plays a significant role in understanding the seismic risks associated with different anthropogenic activities including sub‐surface energy harvesting and sequestration. Here we investigate multi‐cycle modeling of induced seismicity in a complex fault zone hosting a primary fault surrounded by multiple secondary faults with different orientations. The simulation results show spatially heterogeneous distribution of seismicity as well as temporal clustering. While fault criticality remains an important factor, fault network interactions modulate individual fault behavior through stress transfer mechanisms. The complex interplay between pore pressure diffusion, aseismic slip propagation, and elastic stress redistribution creates emergent spatiotemporal patterns that cannot be predicted from single fault models alone. Furthermore, the simulations show that higher confining stress leads to reduced complexity and the seismicity tends to localize on the primary fault. These results may help in assessing induced seismicity hazard and guide safer injection protocols to achieve eco‐friendly developments.
Mia et al. (Fri,) studied this question.