Seismic tomography of Earth’s mantle images abundant slab remnants, often located in close proximity to active subduction systems. The impact of such remnants on the dynamics of subduction remains under explored. Here, we use simulations of multi-material free subduction in a 3-D spherical shell geometry to examine the interaction between visco-plastic slabs and remnants that are positioned above, within and below the mantle transition zone. Depending on their size, negatively buoyant remnants can set up mantle flow of similar strength and length scales as that due to active subduction. As such, we find that remnants located within a few hundred km from a slab tip can locally enhance sinking by up to a factor 2. Remnant location influences trench motion: the trench advances towards a remnant positioned in the mantle wedge region, whereas remnants in the sub-slab region enhance trench retreat. These motions aid in rotating the subducting slab and remnant towards each other, reducing the distance between them, and further enhancing the positive interaction of their mantle flow fields. In this process, the trench develops along-strike variations in shape that are dependent on the remnant’s location. Slab-remnant interactions may explain the poor correlation between subducting plate velocities and subducting plate age found in recent plate tectonic reconstructions. Our results imply that slab-remnant interactions affect the evolution of subducting slabs and trench geometry. Remnant-induced downwelling may also anchor and sustain subduction systems, facilitate subduction initiation, and contribute to plate reorganisation events.

The effects of sphericity are regularly neglected in numerical and laboratory studies that examine the factors controlling subduction dynamics. Most existing studies have been executed in a Cartesian domain, with the small number of simulations undertaken in a spherical shell incorporating plates with an oversimplified rheology, limiting their applicability. Here, we simulate free-subduction of composite visco-plastic plates in 3-D Cartesian and spherical shell domains, to examine the role of sphericity in dictating the dynamics of subduction, and highlight the limitations of Cartesian models. We identify two irreconcilable differences between Cartesian and spherical models, which limit the suitability of Cartesian-based studies: (i) the presence of sidewall boundaries in Cartesian models, which modify the flow regime; and (ii) the reduction of space with depth in spherical shells, alongside the radial gravity direction, which cannot be captured in Cartesian domains. Although Cartesian models generally predict comparable subduction regimes and slab morphologies to their spherical counterparts, there are significant quantitative discrepancies. We find that simulations in Cartesian domains that exceed Earth’s dimensions overestimate trench retreat. Conversely, due to boundary effects, simulations in smaller Cartesian domains overestimate the variation of trench curvature driven by plate width. Importantly, spherical models consistently predict higher sinking velocities and a reduction in slab width with depth, particularly for wider subduction systems, enhancing along-strike slab buckling and trench curvature. Results imply that sphericity must be considered when simulating Earth’s subduction systems, and that it is essential for accurately predicting the dynamics of subduction zones of width ~2400 km or more.

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.

The crust and upper mantle structure of the Greater and Lesser Antilles Arc provides insights into key subduction zone processes in a unique region of slow convergence of old slow-spreading oceanic lithosphere. We use ambient noise tomography gathered from island broadband seismic stations and the temporary ocean bottom seismometer network installed as part of the VoiLA experiment to map crustal and upper mantle shear-wave velocity of the eastern Greater Antilles and the Lesser Antilles Arc. We find sediment thickness, based on the depth to the 2.0 km/s contour in the Grenada and Tobago basins up to 15 km in the south, with thinner sediments near the arc and to the north. We observe thicker crust, based on the depth to the 4.0 km/s velocity contour, beneath the arc platforms with the greatest crustal thickness of around 30 km, likely related to crustal addition from arc volcanism through time. There are distinct low velocity zones (4.2-4.4 km/s) in the mantle wedge (30-50 km depth), beneath the Mona Passage, Guadeloupe-Martinique, and the Grenadines. The Mona passage mantle anomaly may be related to ongoing extension there, while the Guadeloupe-Martinique and Grenadine anomalies are likely related to fluid flux, upwelling, and/or partial melt related to nearby slab features. The location of the Guadeloupe-Martinique anomaly is slightly to the south of the obliquely subducted fracture zones. This feature could be explained by either three-dimensional mantle flow, a gap in the slab, variable slab hydration, and/or melt dynamics including ponding and interactions with the upper plate.

The role of Earth’s spherical geometry in modulating the evolution of subduction zones is poorly understood. Here, we simulate multi-material free-subduction in a 3-D spherical shell domain, to investigate the effect of plate thickness, density (combined approximating age) and width on the evolution of subduction systems. To isolate the role of sphericity, we compare results with equivalent Cartesian models. The first-order predictions of our spherical cases are generally consistent with existing Cartesian studies: (i) slabs retreat more, at a shallower dip, as plate age increases, due to increased bending resistance and sinking rates; and (ii) wider slabs can develop along-strike variations in trench curvature, trending towards a ‘W’-shape, due to toroidal flow at slab edges. We find, however, that these along-strike variations are restricted to older, stronger, retreating slabs. When compared to slabs in Cartesian models, in a spherical domain: (i) slabs descend faster, due to the convergence of downwelling material with depth; (ii) these faster sinking rates reduce the time available for bending at the trench, resulting in effectively stronger slabs; (iii) the curvature of slabs increases their effective strength; and (iv) the curvature of the transition zone tends to enhance slab stagnation. These differences between spherical and Cartesian cases become more prominent as slab width increases. Taken together, our results suggest that Cartesian models are suitable for simulating narrow subduction zones, but spherical models should be utilised when investigating subduction zones wider than ~ 2000 km: at such length-scales, the consequences of Earth’s curvature cannot be ignored.

Constraints on the structure of cratonic lithosphere are essential to improve our understanding of craton formation, evolution and long-term stability. Here, we perform a joint inversion for the thermal and compositional structure of the mantle lithosphere below the South America Platform, using Rayleigh wave group velocities, elevation, and geoid height. Thick thermal lithosphere (200-300 km) is found below the southern Amazonian and São Francisco Cratons and adjoining Parecis Basin and northern Paraná Basin. The southern Rio de la Plata Craton also retains a 200-250 km thick keel. Compositionally, Amazonian, São Francisco and Rio de la Plata lithosphere has a metasomatic and possibly eclogite signature similar to that of North American Proterozoic collision belts. Parecis and northern Paraná lithosphere has likely been altered by Mesozoic plume activity throughout most of its depth, while the rest of the Paraná Basin and the Chaco and Patanal basins appear to have lost the lithospheric root below ~100 km depth that was there during intracratonic basin formation. The low elevation and high geoid of the western Paraná Basin requires a dense (eclogite) layer within the crust/shallow lithosphere, possibly associated with the NeoProterozoic western Paraná Suture Zone and/or Mesozoic plume activity, while topography and geoid of the basins further west and of the western Rio de la Plata craton seem affected by dynamic (subduction-related) topography. Thus, the variable geophysical structure of the platform lithosphere reflects a history that involves besides some stable keels, significant modification and thinning.