No abstract

Subduction of oceanic lithosphere occurs through two modes: subducting plate motion and trench migration. Using a global subduction zone data set and three-dimensional numerical subduction models, we show that slab width (W) controls these modes and the partitioning of subduction between them. Subducting plate velocity scales with W(2/3), whereas trench velocity scales with 1/W. These findings explain the Cenozoic slowdown of the Farallon plate and the decrease in subduction partitioning by its decreasing slab width. The change from Sevier-Laramide orogenesis to Basin and Range extension in North America is also explained by slab width; shortening occurred during wide-slab subduction and overriding-plate-driven trench retreat, whereas extension occurred during intermediate to narrow-slab subduction and slab-driven trench retreat.

Earthquake moment tensors in eastern Pacific (ePac) slabs typically show downdip tensional (DT) axes, whereas in the western Pacific (wPac), they typically show downdip compressional (DC) axes or have mixed orientations indicative of unbending. Prevailing conceptual models emphasize uniform stress/deformation modes, that is, bulk slab stretching or shortening, as the dominant control on intermediate depth seismic expression. In contrast, we propose that a diversity of seismic expression, including DT‐ and DC‐dominated regions, is consistent with expectations of flexural strain accumulation, based on systemic differences in slab geometry. Our analysis reveals two largely unrecognized features of ePac intraslab seismicity. First, earthquake clusters consistent with slab unbending are present in ePac slabs, albeit at much shallower depths than typical of wPac slabs. Second, intermediate depth ePac DT seismicity is strongly localized to the upper half of zones undergoing curvature increase, such as flat slab segments. Our study highlights how the seismic expression of slab flexure is impacted by the relative contribution of brittle and ductile deformation. The strongly asymmetric temperature structure that is preserved in sinking slabs means that seismicity disproportionately records the deformation regime in the colder part of the slab, above the neutral plane of bending. The expression of in‐plane stress may be discernible in terms of a systematic modifying effect on the seismic expression of flexure.

[1] Subduction of oceanic lithosphere occurs through both trenchward subducting plate motion and trench retreat. We investigate how subducting plate velocity, trench velocity and the partitioning of these two velocity components vary for individual subduction zone segments as a function of proximity to the closest lateral slab edge (D SE ). We present a global compilation for 207 trench segments from 17 active subduction zones on Earth and three-dimensional numerical models of progressive free subduction of a single oceanic plate that subducts into a stratified mantle. The results show that the subducting plate velocity is always high (≥5.1 cm/yr (models) and ≥4.2 cm/yr (nature)) and trench velocity is always low (≤2.5 cm/yr (models) and ≤1.7 cm/yr (nature)) in the center of wide subduction zones (D SE > 2200 km). Only in regions close to lateral slab edges (D SE < 1000 km), be it for narrow or wide subduction zones, can the trench velocity exceed 4 cm/yr (models) and 6 cm/yr (nature) and can the subducting plate velocity go below 4 cm/yr (models) and 2 cm/yr (nature). In general, plate velocities, trench velocities and subduction partitioning are much more variable near slab edges than in the center of wide subduction zones owing to other parameters that affect subduction kinematics. We conclude that subduction kinematics can vary considerably along individual subduction zones and that the upper bound values for trench velocity and lower bound values for subducting plate velocity and subduction partitioning at individual subduction zone segments depend critically on D SE .

No abstract