# Episodic back-arc spreading centre jumps controlled by transform fault to overriding plate strength ratio

Nicholas Schliffke, Jeroen van Hunen, Mark B. Allen, Valentina Magni,
Frédéric Gueydan

► To cite this version:

Nicholas Schliffke, Jeroen van Hunen, Mark B. Allen, Valentina Magni, Frédéric Gueydan. Episodic back-arc spreading centre jumps controlled by transform fault to overriding plate strength ratio. Nature Communications, 2022, 13, 10.1038/s41467-022-28228-5. insu-03661251

HAL Id: insu-03661251

https://insu.hal.science/insu-03661251

Submitted on 7 May 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE

https://doi.org/10.1038/s41467-022-28228-5

OPEN

<figure><img src="image_2.png" /></figure>

Check for updates

# Episodic back-arc spreading centre jumps controlled by transform fault to overriding plate strength ratio

Nicholas Schliffke ID <fnref n="1" />✉, Jeroen van Hunen ID <fnref n="1" />, Mark B. Allen ID <fnref n="1" />, Valentina Magni ID ² & Frédéric Gueydan ³

Spreading centre jumps are a common feature of oceanic back-arc basins. Jumps are con-
ventionally suggested to be triggered by plate velocity changes, pre-existing weaknesses, or
punctuated events such as the opening of slab windows. Here, we present 3D numerical
models of back-arc spreading centre jumps evolving naturally in a homogeneous subduction
system surrounded by continents without a trigger event. Spreading centres jump towards
their subduction zone if the distance from trench to spreading centre becomes too long. In
particular, jumps to a new spreading centre occur when the resistance on the boundary
transform faults enabling relative motion of back-arc and neighbouring plates is larger than
the resistance to break the overriding plate closer to trench. Time and distance of spreading
centres jumps are, thus, controlled by the ratio between the transform fault and overriding
plate strengths. Despite being less complex than natural systems, our models explain why
narrow subducting plates (e.g. Calabrian slab), have more frequent and closely-spaced
spreading jumps than wider subduction zones (e.g. Scotia). It also explains why wide back-arc
basins undergo no spreading centre jumps in their life cycle.

<footnote n="1">¹Department of Earth Sciences, Durham University, DH1 3LE Durham, UK. ²The Centre for Earth Evolution and Dynamics (CEED), University of Oslo, SemSaelandsvei 24, PO Box 1048, Blindern, NO-0316 Oslo, Norway. ³Géosciences Montpellier, Université Montpellier, place E. Bataillon 34095, Montpellier cedex 5, France. ✉email: nico.schliffke@web.de</footnote>
The cycle of an oceanic plate from the formation at a sea-floor spreading to consumption at subduction zones is fundamental to plate tectonics¹. Most subduction zone trenches retreat with respect to the lower mantle, i.e. they migrate towards the subducting plate²,³. This 'rollback' exerts tensional stresses to the overriding lithosphere and frequently forms back-arc spreading centres⁴,⁵ that may evolve into small oceanic basins. Formation and location of back-arc spreading centres is controlled by heterogeneities in the buoyancy of subducting plate or strength of the overriding plate: rupturing of the overriding plates is often enhanced by either weak arcs with melt- or fluid-weakening⁶,⁷ or by trench rotation due to subduction of buoyant material⁸⁻¹⁰.

While spreading in most back-arc basins ceases after 5–15 Myr of activity⁸, others experience long-term spreading and episodic back-arc spreading jumps⁸,¹¹,¹², when old back-arc spreading centres are abandoned and new ones form closer to the trench. Such jumps have been recorded in South East Asia (e.g. Philippine Sea plate¹³, Lau Basin¹⁴,¹⁵) and in relatively narrow subduction zones (e.g. Scotia¹⁶⁻¹⁸, Lesser Antilles¹⁹ and Tyrhenian Sea²⁰,²¹). The former regions experienced large-scale plate-reorganisation¹⁴ and subduction polarity flips¹⁴,²², possibly causing transient stress field changes that facilitate back-arc spreading jumps. The latter regions, however, show continuous trench retreat (with respect to fixed neighbouring plates and the upper mantle) without necessarily requiring major plate reorganisations, plate rotations, or significant variations in buoyancy or rheology. In those regions, the main back-arc basin characteristics are large-scale bounding transform faults that decouple them from their neighbouring plates, and STEP²,²³,²⁴ ('subduction-transform-edge-propagator') faults that are created at slab edges, thereby tearing through lithosphere at the plate boundary (Fig. 1).

Previous explanations for the occurrence of spreading jumps invoked internal or external triggers, such as interaction between the oceanic slab and the 660-km discontinuity²¹, opening of lateral slab windows at depth⁹, rearranged mantle flow¹¹, or inherited weaknesses in the overriding plate¹². The role of prominent STEP and transform faults in potentially controlling the back-arc spreading jumps have never been constrained.

This work thus focuses on subduction systems exhibiting back-arc spreading jumps without significant plate reorganisation, but instead associated with STEP- and transform faults. We hypothesize that internal feedback mechanisms within a homogeneous subduction system can lead to episodic back-arc spreading centre jumps without any triggering: as back-arc basins open, the transform faults at their edge become longer and induce more resistance. Back-arc spreading centres may lock up as a consequence, and jump towards a new location closer to the trench, thereby reducing the transform fault length. To test this, we investigate the occurrence and characteristics of back-arc spreading jumps in three-dimensional numerical models with homogenous plates.

## Results

**Controls of trench retreat and back-arc spreading.** Back-arc spreading and trench rollback are intrinsically linked. In models with subducting slabs only and no overriding plate, resisting forces are limited to deformation of the slab and mantle, and in that case trench retreat is a function of slab strength²⁵ and width as well as toroidal mantle flow around slab edges²⁶. Adding a strong or thick overriding plate reduces retreat velocities²⁷,²⁸, and high retreat velocities are most likely for thin overriding plates²⁸,²⁹, or for those ruptured by back-arc spreading centres¹². In addition, the presence of neighbouring plates further reduces trench retreat velocities³⁰. The strength of interfaces to

neighbouring plates, along which the retreating slab edge must tear through lithosphere, controls trench retreat and curvature³¹, high tear resistance causes strong trench curvature and slow rollback³¹. In contrast, low tear resistance facilitates trench retreat.

The transition from a rigid overriding plate to back-arc spreading requires other processes or weaknesses in addition to extensional stresses from trench rollback alone⁹,¹¹. Melt- or fluid percolation⁶,³² or a thin volcanic arc⁷ can substantially weaken the overriding plate and support plate rupturing. Moreover, trench rotation localising tensional stresses locally after adjacent continental collision, has shown to trigger back-arc basins formation⁹,¹⁰. Mechanisms for opening back-arc basins have been explored in some detail in analogue and numerical models, but modelling the mechanisms for subsequent back-arc spreading jumps has received much less attention.

Tracing back-arc spreading centre jumps across geological timescales is difficult because abandoned back-arc basins are often deformed or destroyed by subduction. Examples of recent inactive back-arc basins that are currently subducting are the Sulu Basin (which opened between 19 and 15 Ma during subduction of the Celebes Sea⁸) and the North Loyalty Basin (opened 44–35 Ma in the New Caledonia subduction zone⁸). However, several active and recent examples do provide data for the length scales and repeat rates of ridge jumps. In the Central Mediterranean subduction zone, the southeastward rollback of the Calabrian slab caused overriding plate extension in the Tyrhenian Sea since about 14 Ma²¹,³³ (Fig. 1a). Slab rollback is enabled by a STEP fault in the south³⁴ Episodic pulses of fast trench rollback²¹ coincide with jumps from the extensional Cornaglia Basin to spreading in the Vavilov Basin (4.3–2.6 Ma)³⁵ and most recently the Marsili Basin (~2–1 Ma)³⁶,³⁷. These pulses have been suggested to be caused by slab stagnation on top of the lower mantle²¹ and/or opening of slab windows in the previously wider slab⁹. The Scotia subduction zone has retreated eastwards between Antarctica and South America since ~40 Ma¹⁶,¹⁷ with two major back-arc spreading centres (Fig. 1b). The present-day configuration was established with the opening of the East Scotia ridge at ~6–8 Ma, close to the eastward retreating trench, and the West Scotia Ridge simultaneously becoming inactive¹⁸,³⁸ (see comment in Supplementary Tables). Other than inherited weaknesses or overriding plate heterogeneities¹², no triggers for this ridge jump have been identified. Although tectonically complex with local plate rotations, in the Caribbean, an Eocene back-arc spreading jump from the Venezuelan Basin to the Grenada/Tobago Basin has been proposed to be essentially caused by continuous rollback between bounding STEP-faults¹⁹ (Fig. 1c).

**Modelling results.** We present spreading centre jumps in 3D-models with a simple visco-plastic rheology without any compositional strength variation (see Methodology). A narrow oceanic basin (along-strike width varied between 400–1500 km) is situated between two continental plates and subducting below an overriding continental plate (Fig. 2a). To establish gravitationally driven subduction, a slab segment is initially placed in the mantle and between short prescribed transform faults.

Initially, subduction and rollback can occur with little resistance between the pre-existing transform faults in all models due to decoupling of subducting and overriding plate from neighbouring plates (see Supplementary Figures for initial evolution). To study natural plate interaction, we limit the extent of the prescribed transform faults and allow for their spontaneous propagation. Once rollback reaches the end of the prescribed transform faults, further rollback of the subducting slab is only feasible by tearing of lithosphere between subducting and
<figure><img src="image_3.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_6.png" /></figure>

**Fig. 1 Mechanism sketch.** **a–c** Extinct and active back-arc spreading centre in the narrow Tyrrhenian Sea (adapted from Gueydan et al.³³), Scotia (adapted from Eagles & Jokat¹⁸, and Lesser Antilles at ~40 Ma (adapted from Allen et al.¹⁹. All subduction zones retreat within the bounds of STEP- or similar transform faults. **d** Sketch of a spreading centre jump in a retreating subduction zone, bound laterally by transform faults. The subducting slab tears along STEP (subduction-transform-edge-propagator) faults, whereas the back-arc moves relative to the neighbouring along transform faults (grey shaded regions). Within a short period is old spreading centre and the associate volcanic arc extinct, and a new spreading centre forms closer to the retreating trench.

neighbouring plates at STEP faults (Fig. 1d, Supplementary Fig. 1). At this point, the first back-arc basin is formed: high tear resistance at the STEP faults reduces lateral rollback at the plate interface³¹,³⁹, and the resulting trench curvature localizes stresses and ruptures the upper plate. The initial opening of a back-arc basin is facilitated and simplified by this ending of the prescribed transform faults, but subsequent back-arc spreading will involve STEP fault propagation that will form an essential ingredient for back-arc spreading jumps.

Models with oceanic plate widths between 600 and 1000 km show ongoing rollback with back-arc spreading (Fig. 2b). High strain rates at the STEP faults reduce the effective viscosity due to the nonlinear rheology (see Methodology) and allow for tearing

of the slab edge through the lithosphere. Resulting continued trench retreat increases the distance between trench and back-arc spreading centre. Stress and strain rate concentrate at the plate edges and lead to natural plate-like behaviour and decoupling of overriding and neighbouring plates along transform faults. (Fig. 2b). The resulting low frictional resistance of a young and small back-arc plate allows for rollback rates to increase from the initial 3 cm/yr up to 7 cm/yr (Fig. 3a) in this first 5–10 Myrs of back-arc spreading (Fig. 3a).

All deformation in the model, both in the mantle and along plate boundaries at the surface, is entirely driven internally by subduction and localises in the bending downgoing plate, at the trench, at the STEP faults, and at the transform faults between the
<figure><img src="image_7.png" /><figcaption>**Fig. 2 Model evolution.** Model setup (a) and evolution (b-e) of the second invariant of the strain rate ($\dot{\epsilon}_{ll} = \sqrt{0.5\epsilon_{ij}^2}$) and temperature for an 800 km wide slab which is used as reference model. Arrows represent absolute motion with reference to the fixed model edge. Model evolution plots show an inset area of panel (a) of surface strain rate in the centre of the model domain and the vertical-temperature profile up to 200 km depth along the red line. **a** A narrow slab subducts between two buoyant continents and below an overriding continental plate. The slab is initially placed up to 150 km deep in the mantle and between two prescribed transform faults ending after 660 km in x-direction (dashed red lines) to ensure subduction initiation. **b** After opening of an initial back-arc basin, transform faults form, connecting the trench and spreading centre. **c** At a critical distance, increased coupling between back-arc and neighbouring plates forces the transform faults to become inactive and localization of stresses closer to the trench. **d** A rapid ridge jump occurs. **e** The slab width is decreased and new transform faults are formed.</figcaption></figure>

back-arc spreading centre and neighbouring plates. Of these, only deformation along the transform faults increases with time as the back-arc basin grows, and, as a result, rollback rates reduce (Fig. 3a shows a decrease from 7 to 2–3 cm/yr). Given that the total amount of available energy to overcome resistance against deformation is limited, the transform faults eventually reach a threshold length after which they lock up (Fig. 2c). Ongoing rollback then localises stresses closer to the trench where a new spreading centre fully develops within 2–4 Myr (Fig. 2d, e), while the old spreading centre is abandoned rapidly. This new spreading centre forms where the toroidal flow around the slab edges causes the highest stresses at the base of the overriding plate (see Supplementary Fig. 3)$^{40}$. After the spreading centre jump, the slab width shortens and total deformation at surface is reduced (Fig. 2d), which causes a pulse of increased rollback rate by about 1 cm/yr (Fig. 3a).

The timing of the back-arc spreading jump is governed by the relative strengths of the transform faults and the overriding plate. Wider overriding plates require more energy to break than narrow ones, so that old back-arc spreading centres remain active for longer. Indeed, within the set of models with back-arc spreading jumps, the maximum transform fault length and lifespan of a

back-arc spreading centre both linearly increase with overriding plate width (Fig. 3). A first-order theoretical derivation of supporting this relationship is presented in the Supplementary Discussion. Our conclusion that limited deformation along STEP- and transform faults is crucial for back-arc spreading jumps is supported by model cases that do not develop jumps, and models for a 400 km and a 1500 km wide slab, as well as a model with a weak plate interface along the entire neighbouring plate, are presented in Supplementary Figs. 4–6. Each model has a different reason for the absence of jumps: a 400 km wide slab is unable to overcome tear resistance at STEP faults, i.e. the strain rate induced by the slab pull is too low to weaken the lithosphere and prevents rollback and subduction completely (Fig. 3 and Supplementary Fig. 4). In contrast, the 1500 km wide slab retreats continuously until it reaches the model edge (after >4000 km retreat) without a spreading jump (Supplementary Fig. 5). Although a jump might occur in a longer model domain, such lengthy uninterrupted rollback is unlikely to happen within a typical lifetime of an oceanic plate.

No spreading jumps develop in models with prescribed weak zones (i.e. transform faults) along neighbouring plates (Supplementary Fig. 6). This demonstrates the role of STEP fault tearing
<figure><img src="image_8.png" /><figcaption>**Fig. 3** Temporal evolution. **a** Trench velocities and **b** transform fault lengths for models with various plate widths. Model times are relative to the first back-arc basin opening. Opening of back-arc basins increase rollback rates (i.e. negative trench velocities) temporally. While the overriding plate and transforms grow with ongoing rollback, rollback rates decrease. Spreading jumps are displayed in the suddenly decreased transform length. Subduction ceases in the 400 km model, as well as in the 600 km model after the back-arc spreading centre jump, due to the inability to tear through lithosphere at plates interfaces (see Supplementary Information).</figcaption></figure>

in the development of spreading jumps. With prescribed weaknesses, which are weak and require no energy to form, both the slab and overriding plate can retreat along neighbouring plate effortlessly, making the total transform fault distance between trench and back-arc basin irrelevant. Subducting and overriding plate are basically independent from neighbouring plates and thus spreading centres jumps are not feasible.

Discussion

Our models show that the ratio of transform fault versus over-
riding plate strengths governs the wavelength and frequency of
back-arc spreading jumps. Narrow plates can only maintain
decoupling of back-arc and neighbouring plates along transform
faults of finite length. Furthermore, we show why the widest
subduction zones are unlikely to undergo any spreading jumps
with this mechanism and possibly require external factors or
heterogeneities in the plate (e.g., plate re-organisation, lower
mantle slab penetration, nearby collision, along-trench variations
of trench retreat velocities). Moreover, much wider plates
(>2000 km) are likely to have a different trench curvature and
velocity along strike due to mantle flow associated with slab
rollback<sup>26,41</sup>. To test the applicability of our models, we compare
the results to natural back-arcs that have experienced spreading
jumps and to those with a single spreading centre (Fig. 1 and 4).
Two characteristics are used: (1) The distance between the trench
and the centre of the back-arc extension (Fig. 4a, see our Sup-
plementary Tables), which, for our models and natural cases, is
equivalent to the transform fault length, and (2) the back-arc
rifting/spreading duration (Fig. 4b), which is relatively well con-
strained in nature by dating of appropriate sedimentary rocks and
magnetic seafloor lineations.

<figure><img src="image_9.png" /><figcaption>**Fig. 4** Comparison to data. Comparison of models and observations of back-arc spreading with and without back-arc spreading jumps. Timing in the natural data comprises the entire rifting and spreading duration. **a** The largest distance between trench and back-arc spreading centre. **b** Spreading duration of back-arcs. For models without spreading centre jumps (non-filled red circles), we measure the spreading period of the first back-arc. Models with back-arc spreading jumps form a linear trend, while the ones without spreading jumps plot below this trend. Observational data of back-arc spreading jumps plot close to the regression line, while all other active or recently inactive basins plot below it.</figcaption></figure>

Figure 4 shows that plate widths are linearly correlated to both
characteristics in all models with back-arc spreading jumps, while
models without jumps (400 km and 1500 km) plot clearly below
the regression line. Data points from observations show a pri-
marily short back-arc spreading duration of 3–15 Myrs which
also explains short distances between trench and spreading centre
(<800 km). Exceptions are subduction systems retreating between
continental margins and with spreading jumps, i.e. Tyrrhenian
Sea (Calabria), Lesser Antilles and Scotia, which all plot close to
the regression line. The regression line thus represents an
important threshold: given that we expect spreading jumps to
open a new spreading centre and abandon an old spreading
centre close to this line, we infer that no active back-arc spreading
centre should plot above the line, consistent with natural data.

Our models suggest that the longest subduction zones with
back-arc basins, such as the New Hebrides (~1600 km) are not
likely to jump, as they would require long-term, uninterrupted
rollback of >40 Myrs, assuming our extrapolation is valid at such
subduction zone lengths. According to our models, the narrowest
plate widths for spreading jumps in the central Mediterranean
(<400 km) should not be subducting at all, as lithospheric tearing
at STEPs is not feasible. Observations of such subducting narrow
plates show, however, that the minimum plate width is lower in
nature than in our models. Additional processes reducing tear
resistance in nature, and hence also enabling back-arc spreading
jumps in narrow subduction systems, could be weaker plate
boundaries³¹,³⁹, 'true' brittle deformation at the STEP-faults (see Methodology for rheology), reactivation of old-faults or strain-weakening. Very narrow subduction systems have lower energy dissipation, in which case the exact details of the plate boundary friction become more important, and models could potentially benefit from some of these more complex rheological features².

In conclusion, back-arc spreading centres remain active for as long as slip along the lengthening transform faults remains energetically favourable over rifting and rupturing the overriding plate in a new location closer to the subduction zone. Once a threshold length is reached by ongoing rollback, the old spreading centre is locked and a new spreading centre is formed. Given that the available potential energy is proportional to the plate width⁹, frictional forces along transform faults might be small compared to the driving buoyancy forces for wide subduction zones, but they are important and non-negligible for narrow subduction zones.

**Data availability**

See Supplementary Information for Data used in this study. No new Data was created.

**Code availability**

CITCOM is available upon request.

Received: 13 July 2020; Accepted: 10 January 2022;

Published online: 31 January 2022

**References**

1. Crameri, F., Conrad, C. P., Montési, L. & Lithgow-Bertelloni, C. R. The dynamic life of an oceanic plate. *Tectonophysics* **760**, 107–135 (2019).

2. Munch, J., Gerya, T. & Ueda, K. Oceanic crust recycling controlled by weakening at slab edges. *Nat. Commun.* **11**, 1–6 (2020).

3. Garfunkel, Z., Anderson, C. A. & Schubert, G. Mantle circulation and the lateral migration of subducted slabs. *J. Geophys. Res. Solid Earth* **91**, 7205–7223 (1986).

4. Dewey, J. F. Episodicity, sequence, and style at convergent plate boundaries. *Cont. crust its Miner. Depos.* **20**, 553–573 (1980).

5. Dvorkin, J., Nur, A., Mavko, G. & Ben-Avraham, Z. Narrow subducting slabs and the origin of backarc basins. *Tectonophysics* **227**, 63–79 (1993).

6. Baitsch-Ghirardello, B., Gerya, T. V. & Burg, J.-P. Geodynamic regimes of intra-oceanic subduction: Implications for arc extension vs. shortening processes. *Gondwana Res.* **25**, 546–560 (2014).

7. Boutelier, D. & Cruden, A. Slab rollback rate and trench curvature controlled by arc deformation. *Geology* **41**, 911–914 (2013).

8. Wallace, L. M., Ellis, S. & Mann, P. Collisional model for rapid fore-arc block rotations, arc curvature, and episodic back-arc rifting in subduction settings. *Geochem. Geophys. Geosyst.* **10**, Q05001 (2009).

9. Magni, V., Faccenna, C., van Hunen, J. & Funiciello, F. How collision triggers backarc extension: insight into Mediterranean style of extension from 3-D numerical models. *Geology* **42**, 511–514 (2014).

10. Menant, A., Sternai, P., Jolivet, L., Guillou-Frottier, L. & Gerya, T. 3D numerical modeling of mantle flow, crustal dynamics and magma genesis associated with slab roll-back and tearing: The eastern Mediterranean case. *Earth Planet. Sci. Lett.* **442**, 93–107 (2016).

11. Sdrolias, M. & Müller, R. D. Controls on back-arc basin formation. *Geochem. Geophys. Geosyst.* **7**, Q04016 (2006).

12. Clark, S. R., Stegman, D. & Müller, R. D. Episodicity in back-arc tectonic regimes. *Phys. Earth Planet. Inter.* **171**, 265–279 (2008).

13. Deschamps, A. et al. Propagators and ridge jumps in a back-arc basin, the West Philippine Basin. *Terra Nov.* **20**, 327–332 (2008).

14. Schellart, W. P., Lister, G. S. & Toy, V. G. A Late Cretaceous and Cenozoic reconstruction of the Southwest Pacific region: tectonics controlled by subduction and slab rollback processes. *Earth-Sci. Rev.* **76**, 191–233 (2006).

15. Parson, L. M., Pearce, J. A., Murton, B. J. & Hodkinson, R. A. Role of ridge jumps and ridge propagation in the tectonic evolution of the Lau back-arc basin, southwest Pacific. *Geology* **18**, 470–473 (1990).

16. Maldonado, A. et al. A model of oceanic development by ridge jumping: opening of the Scotia Sea. *Glob. Planet. Change* **123**, 152–173 (2014).

17. Maldonado, A., Dalziel, I. W. D. & Leat, P. T. The global relevance of the Scotia Arc: an introduction. *Glob. Planet. Change* **125**, A1–A8 (2015).

18. Eagles, G. & Jokat, W. Tectonic reconstructions for paleobathymetry in Drake Passage. *Tectonophysics* **611**, 28–50 (2014).

19. Allen, R. W. et al. The role of arc migration in the development of the Lesser Antilles: a new tectonic model for the Cenozoic evolution of the eastern Caribbean. *Geology* **47**, 891–895 (2019).

20. Spadini, G. & Podladchikov, Y. Spacing of consecutive normal faulting in the lithosphere: a dynamic model for rift axis jumping (Tyrrenian Sea). *Earth Planet. Sci. Lett.* **144**, 21–34 (1996).

21. Faccenna, C., Funiciello, F., Giardini, D. & Lucente, P. Episodic back-arc extension during restricted mantle convection in the Central Mediterranean. *Earth Planet. Sci. Lett.* **187**, 105–116 (2001).

22. Teng, L. S., Lee, C. T., Tsai, Y. B. & Hsiao, L.-Y. Slab breakoff as a mechanism for flipping of subduction polarity in Taiwan. *Geology* **28**, 155–158 (2000).

23. Govers, R. & Wortel, M. J. R. Lithosphere tearing at STEP faults: Response to edges of subduction zones. *Earth Planet. Sci. Lett.* **236**, 505–523 (2005).

24. Rosenbaum, G., Gasparon, M., Lucente, F. P., Peccerillo, A. & Miller, M. S. Kinematics of slab tear faults during subduction segmentation and implications for Italian magmatism. *Tectonics* **27**, TC2008 https://doi.org/10.1029/2007TC002143 (2008).

25. Stegman, D. R., Farrington, R., Capitanio, F. A. & Schellart, W. P. A regime diagram for subduction styles from 3-D numerical models of free subduction. *Tectonophysics* **483**, 29–45 (2010).

26. Schellart, W. P., Freeman, J., Stegman, D. R., Moresi, L. & May, D. Evolution and diversity of subduction zones controlled by slab width. *Nature* **446**, 308 (2007).

27. Sharples, W., Jadamec, M. A., Moresi, L.-N. & Capitanio, F. A. Overriding plate controls on subduction evolution. *J. Geophys. Res. Solid Earth* **119**, 6684–6704 (2014).

28. Holt, A. F., Becker, T. W. & Buffett, B. A. Trench migration and overriding plate stress in dynamic subduction models. *Geophys. J. Int.* **201**, 172–192 (2015).

29. Hertgen, S. et al. Influence of the thickness of the overriding plate on convergence zone dynamics. *Geochem., Geophys. Geosyst.* **21**, e2019GC008678 (2020).

30. Yamato, P., Husson, L., Braun, J., Loiselet, C. & Thieulot, C. Influence of surrounding plates on 3D subduction dynamics. *Geophys. Res. Lett.* **36**, L07303 https://doi.org/10.1029/2008GL036942 (2009).

31. Hale, A. J. et al. Dynamics of slab tear faults: Insights from numerical modelling. *Tectonophysics* **483**, 58–70 (2010).

32. Nakakuki, T. & Mura, E. Dynamics of slab rollback and induced back-arc basin formation. *Earth Planet. Sci. Lett.* **361**, 287–297 (2013).

33. Gueydan, F., Brun, J.-P., Phillippon, M. & Noury, M. Sequential extension as a record of Corsica Rotation during Apennines slab roll-back. *Tectonophysics* **710**, 149–161 (2017).

34. Polonia, A. et al. The Ionian and Alfeo-Etna fault zones: New segments of an evolving plate boundary in the central Mediterranean Sea? *Tectonophysics* **675**, 69–90 (2016).

35. Kastens, K. A., Mascle, J. & others. The geological evolution of the Tyrrenian Sea: an introduction to the scientific results of ODP Leg 107. in *Proc. Ocean Drilling Program, Scientific Results* **107**, 3–26 (1990).

36. Guillaume, B., Funiciello, F., Faccenna, C., Martinod, J. & Olivetti, V. Spreading pulses of the Tyrrenian Sea during the narrowing of the Calabrian slab. *Geology* **38**, 819–822 (2010).

37. Nicolosi, I., Speranza, F. & Chiappini, M. Ultrafast oceanic spreading of the Marsili Basin, southern Tyrrenian Sea: evidence from magnetic anomaly analysis. *Geology* **34**, 717–720 (2006).

38. Barker, P. F. & Hill, I. A. Back-arc extension in the Scotia Sea. *Philos. Trans. R. Soc. Lond. Ser. A, Math. Phys. Sci.* **300**, 249–262 (1981).

39. Nijholt, N. & Govers, R. The role of passive margins on the evolution of Subduction-Transform Edge Propagators (STEPs). *J. Geophys. Res. Solid Earth* **120**, 7203–7230 (2015).

40. Schellart, W. P. & Moresi, L. A new driving mechanism for backarc extension and backarc shortening through slab sinking induced toroidal and poloidal mantle flow: Results from dynamic subduction models with an overriding plate. *J. Geophys. Res. Solid Earth* **118**, 3221–3248 (2013).

41. Strak, V. & Schellart, W. P. Control of slab width on subduction-induced upper mantle flow and associated upwellings: insights from analog models. *J. Geophys. Res. Solid Earth* **121**, 4641–4654 (2016).

**Acknowledgements**

This work has been supported by EU FP7 Marie Curie Initial Training Network 'Subi-top', grant agreement No. 674899. JvH acknowledges funding from ERC (grants NE/K010824/1 and NE/M000281/1); V.M. acknowledges support from the Research Council of Norway through its Centres of Excellence funding scheme, Project Number 223272. MBA acknowledges ERC grant NE/H021620/1. We thank Robert Allen for valuable input and Carmen Gaina for fruitful discussions. This work made use of the computational facilities of Hamilton HPC at Durham University.
## Author contributions

N.S. set up the project and the numerical models, and, wrote the manuscript. J.v.H. provided advice on the numerical modelling and contributed to writing the manuscript. M.A. provided the sketches and geological advice. He further provided support in writing the manuscript. V.M. provided the base model on which these models were set up and contributed in writing the manuscript. F.G. provided data and advice on geological implications and contributed in writing the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-022-28228-5.

**Correspondence** and requests for materials should be addressed to Nicholas Schliffke.

**Peer review information** *Nature Communications* thanks Taras Gerya, Douwe van Hinsbergen, and Vincent Strak for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022