During subduction, the oceanic lithosphere of one tectonic plate dives beneath another at a convergent margin and is recycled into Earth's mantle (e.g., Kearey et al., 2009;Stern, 2002). As subducted slabs descend, their negative buoyancy provides a key driving force for plate tectonics, and they continue to influence surface processes in a number of ways (e.g.,

Abstract. Subducting slabs are an important driver of plate motions, yet the relative importance of different forces in governing subduction motions and styles remains incompletely understood. Basal drag has been proposed to be a minor contributor to subduction forcing because of the lack of correlation between plate size and velocity in observed and reconstructed plate motions. Furthermore, in single subduction system models, low basal drag leads to subduction behaviour most consistent with the observation that trench migration velocities are generally low compared to convergence velocities. By contrast, analytical calculations and global mantle flow models indicate basal drag can be substantial. In this study, we revisit this problem by examining the drag at the base of the lithosphere, for a single subduction system, in 2D models with a free trench and composite non-linear rheology. We compare the behaviour of short and long plates for a range of asthenospheric and lithospheric rheologies. We reproduce results from previous modelling studies, including low ratios of trench over plate motions. However, we also find that any combination of asthenosphere and lithosphere viscosity that produces Earth-like subduction behaviour leads to a correlation of velocities with plate size, due to the role of basal drag. By examining Cenozoic plate motion reconstructions, we find that slab age and plate size are positively correlated: higher slab pull for older plates tends to be offset by higher basal drag below these larger plates. This, in part, explains the lack of plate velocity–size correlation in observations, despite the important role of basal drag in the subduction force balance.

The subduction of positively buoyant features has been implicated in the development of flat and shallow dipping slabs, the formation of cusps in trench geometry, and the cessation of associated arc magmatism. However, how such buoyant anomalies influence subduction dynamics to produce these different tectonic expressions remains debated. In this paper, using a series of multi-material 3-D simulations of free subduction, we investigate how linear buoyant ridges modify subduction dynamics, in particular downgoing plate velocities, trench motions and slab morphology. We examine the sensitivity of results to downgoing plate age (affecting buoyancy and strength), ridge buoyancy and ridge location along the trench, finding that buoyant ridges can locally change slab sinking and trench retreat rates, in turn modifying the evolution of slab morphology at depth and trench shape at the surface. In all cases examined, trench retreat is reduced, or switches to trench advance, where the ridge subducts. These effects depend strongly on downgoing plate age: on young, weak plates, the change in trench shape is more localised than on old, strong plates. Slab shallowing at the ridge only occurs for young plates, while the stronger and more negatively buoyant older plates pull down the ridge at a steeper angle than the rest of the slab. On old plates, ridges located near regions of trench stagnation or advance, which typically develop in wide slabs, have a stronger effect on trench and slab shape. The combined effects of buoyant feature location, subducting plate age and overriding plate properties can result in a range of responses: from mainly trench deformation, through local slab shallowing, to the formation of a flat slab, a variation in expressions also observed on Earth.

Many of the factors expected to control the dynamics and evolution of Earth's subduction zones are under-explored in an Earth-like spherical geometry. Here, we simulate multi-material free-subduction of a complex rheology slab in a 3-D spherical shell domain, to investigate the effect of plate age (simulated by covarying plate thickness and density) and width on the evolution of subduction systems. We find that the first-order predictions of our spherical cases are generally consistent with existing Cartesian studies: (i) as subducting plate age increases, slabs retreat more and subduct at a shallower dip angle, due to increased bending resistance and sinking rates; and (ii) wider slabs can develop along-strike variations in trench curvature due to toroidal flow at slab edges, trending towards a 'W'-shaped trench with increasing slab width. We find, however, that these along-strike variations are restricted to older, stronger, retreating slabs:. Younger slabs that drive minimal trench motion remain relatively straight along the length of the subduction zone. We summarise our results into a regime diagram, which highlights how slab age modulates the effect of slab width, and present examples of the evolutionary history of subduction zones that are consistent with our model predictions.

Subduction is the process where oceanic lithosphere descends into the underlying mantle at a convergent plate boundary, providing the major buoyancy force that drives plate tectonics and mantle convection (e.g., Forsyth