A closer look at the relationship between slab (un)bending and double seismic zone seismicity

The origin of double seismic zones (DSZs), parallel planes of intraslab seismicity observed in many subduction zones around the globe, is still highly debated. While most researchers assume that fluid release from prograde metamorphic reactions in the slab is an important control on DSZ occurrence, the role of slab unbending is currently unclear.
Slab bending at the outer rise is instrumental in hydrating the downgoing oceanic plate through bend faulting, and is evident from earthquake focal mechanisms (prevalence of shallow normal faulting events). Observations from NE Japan show that focal mechanisms of DSZ earthquakes are downdip compressive in the upper and downdip extensive in the lower plane of the DSZ, which strongly hints at slab unbending. This coincidence of slab unbending and DSZ seismicity in NE Japan has given rise to several models in which unbending forces are a prerequisite for DSZ occurrence.

To globally test a potential correlation of slab unbending with DSZ seismicity, we derived downdip slab surface curvatures on trench-perpendicular profiles every 50 km along all major oceanic slabs using the slab2 grids of slab surface depth. We here make a steady-state assumption, i.e. we assume that the slab geometry is relatively constant with time, so that the downdip gradient of slab curvature corresponds to slab (un)bending. We compiled the loci and depth extent of all DSZ observations avalable in literature, and compare these to the slab bending or unbending estimates.

Preliminary results indicate that while there is a clear correspondence between the depth of slab unbending to DSZ seismicity in the Japan-Kurile slab, most other slabs do not show this correlation. Moreover, some DSZs deviate from the above-mentioned focal mechanism pattern and exhibit downdip extension in both planes (e.g. Northern Chile, New Zealand). It appears that the global variability of slab geometries in the depth range 50-200 km is larger than anticipated, and DSZ seismicity is not limited to slabs where unbending is prevalent at these depths. The Northern Chile case is especially interesting because focal mechanisms there not only do not fit the pattern observed in NE Japan, but also can not be explained with the current slab geometry alone. This could indicate a direct influence of ongoing metamorphic reactions on focal mechanisms (e.g. via volume reduction and densification), or it may be a hint that our steady-state assumption is invalid for the Nazca slab here (i.e. that it is in the process of changing its geometry).