Intermediate-depth seismicity and intraslab stress changes due to outer-rise faulting

Intermediate-depth seismicity is common in subducting slabs and the seismicity rate shows some statistically significant yet enigmatic global positive correlation with the maximal throw of outer-rise normal faults. Here, we have simulated the formation and subduction of outer-rise faults, using 2D thermomechanical numerical models of intra-oceanic subduction with coupled brittle-ductile damage of bending plates. We observed that outer-rise faults are formed episodically during slab segmentation and their maximal throw grows with time. When been subducted to intermediate depth, these faults are locally reactivated either by i) slab unbending/bending, simultaneous to the formation of new outer-rise faults or ii) episodic interplate coupling related to the rugged morphology of the faulted downgoing plate. Faults reactivation is concurrent with a local, transient deviatoric stress increase in intraslab domains among these structures. We suggest that slab domains affected by stress increase could be the appropriate location where potential brittle deformation can occur, generating intermediate-depth intraslab earthquakes, that are predominantly localized in heterogeneous regions of dense faulting formed within slab-segments boundaries. The temporal coincidence of stress growth at intermediate depths and throw-growth of, newly-formed, outer-rise faults at the surface may possibly explain the observed global positive correlation of intermediate-depth seismicity rate with maximal fault throw.