An official website of the United States government
Abstract Mantle plumes are typically considered secondary features of mantle convection, yet their surface effects over Earth's evolution may have been significant. We use 2‐D convection models to show that mantle plumes can in fact cause the termination of a subduction zone. This extreme case of plume‐slab interaction is found when the slab is readily weakened, for example, by damage‐type rheology, and the subducting slab is young. We posit that this mechanism may be relevant, particularly for the early Earth, and a subdued version of these plume‐slab interactions may remain relevant for modern subduction zones. Such core‐mantle boundary–surface interactions may be behind some of the complexity of tomographically imaged mantle structures, for example, in South America. More generally, plume “talk back” to subduction zones may make plate tectonics more episodic.
more »
« less
Plume‐Driven Subduction Termination in 3‐D Mantle Convection Models
Abstract The effect of mantle plumes is secondary to that of subducting slabs for modern plate tectonics when considering plate driving forces. However, the impact of plumes on tectonics and planetary surface evolution may nonetheless have been significant. We use numerical mantle convection models in a 3‐D spherical chunk geometry with damage rheology to study some of the dynamics of plume‐slab interactions. Substantiating our earlier 2‐D results, we observe a range of interaction scenarios, and that the plume‐driven subduction terminations we had identified earlier persist in more realistic convective flow. We analyze the dynamics of plume affected subduction, including in terms of their geometry, frequency, and the overall effect of plumes on surface dynamics as a function of the fraction of internal to bottom heating. Some versions of such plume‐slab interplay may be relevant for geologic events, for example, for the inferred ∼183 Ma Karoo large igneous province formation and associated slab disruption. More recent examples may include the impingement of the Afar plume underneath Africa leading to disruption of the Hellenic slab, and the current complex structure imaged for the subduction of the Nazca plate under South America. Our results imply that plumes may play a significant role not just in kick‐starting plate tectonics, but also in major modifications of slab‐driven plate motions, including for the present‐day mantle.
more »
« less
- Award ID(s):
- 2121666
- PAR ID:
- 10522216
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geochemistry, Geophysics, Geosystems
- Volume:
- 25
- Issue:
- 7
- ISSN:
- 1525-2027
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
SUMMARY Tectonic plate motions predominantly result from a balance between the potential energy change of the subducting slab and viscous dissipation in the mantle, bending lithosphere and slab–upper plate interface. A wide range of observations from active subduction zones and exhumed rocks suggest that subduction interface shear zone rheology is sensitive to the composition of subducting crustal material—for example, sediments versus mafic igneous oceanic crust. Here we use 2-D numerical models of dynamically consistent subduction to systematically investigate how subduction interface viscosity influences large-scale subduction kinematics and dynamics. Our model consists of an oceanic slab subducting beneath an overriding continental plate. The slab includes an oceanic crustal/weak layer that controls the rheology of the interface. We implement a range of slab and interface strengths and explore how the kinematics respond for an initial upper mantle slab stage, and subsequent quasi-steady-state ponding near a viscosity jump at the 660-km-discontinuity. If material properties are suitably averaged, our results confirm the effect of interface strength on plate motions as based on simplified viscous dissipation analysis: a ∼2 order of magnitude increase in interface viscosity can decrease convergence speeds by ∼1 order of magnitude. However, the full dynamic solutions show a range of interesting behaviour including an interplay between interface strength and overriding plate topography and an end-member weak interface-weak slab case that results in slab break-off/tearing. Additionally, for models with a spatially limited, weak sediment strip embedded in regular interface material, as might be expected for the subduction of different types of oceanic materials through Earth’s history, the transient response of enhanced rollback and subduction velocity is different for strong and weak slabs. Our work substantiates earlier suggestions as to the importance of the plate interface, and expands the range of quantifiable links between plate reorganizations, the nature of the incoming and overriding plate and the potential geological record.more » « less
-
Abstract The Tonga‐Kermadec subduction zone exhibits the fastest observed trench retreat and convergence near its northern end. However, a paradox exists: despite the rapid trench retreat, the Tonga slab maintains a relatively steep dip angle above 400 km depth. The slab turns flat around 400 km, then steepening again until encountering a stagnant segment near 670 km. Despite its significance for understanding slab dynamics, no existing numerical model has successfully demonstrated how such a distinct slab morphology can be generated under the fast convergence. Here we run subduction models that successfully reproduce the slab geometries while incorporating the observed subduction rate. We use a hybrid velocity boundary condition, imposing velocities on the arc and subducting plate while allowing the overriding plate to respond freely. This approach is crucial for achieving a good match between the modeled and observed Tonga slab. The results explain how the detailed slab structure is highly sensitive to physical parameters including the seafloor age and the mantle viscosity. Notably, a nonlinear rheology, where dislocation creep reduces upper mantle viscosity under strong mantle flow, is essential. The weakened upper mantle allows for a faster slab sinking rate, which explains the large dip angle. Our findings highlight the utilizing rheological parameters that lead to extreme viscosity variations within numerical models to achieve an accurate representation of complex subduction systems like the Tonga‐Kermadec zone. Our study opens new avenues for further study of ocean‐ocean subduction systems, advancing our understanding of their role in shaping regional and global tectonics.more » « less
-
Abstract The tectonic configuration of the Caribbean plate is defined by inward‐dipping double subduction at its boundaries with the North American and Cocos plates. This geometry resulted from a Paleogene plate reorganization, which involved the abandonment of an older subduction system, the Great Arc of the Caribbean (GAC), and conversion into a transform margin during Lesser Antilles (LA) arc formation. Previous models suggest that a collision between the GAC and the Bahamas platform along the North American passive margin caused this event. However, geological and geophysical constraints from the Greater Antilles do not show a large‐scale compressional episode that should correspond to such a collision. We propose an alternative model for the evolution of the region where lower mantle penetration of the Farallon slab promotes the onset of subduction at the LA. We integrate tectonic constraints with seismic tomography to analyze the timing and dynamics of the reorganization, showing that the onset of LA subduction corresponds to the timing of Farallon/Cocos slab penetration. With numerical subduction models, we explore whether slab penetration constitutes a dynamically feasible set of mechanisms to initiate subduction in the overriding plate. In our models, when the first slab (Farallon/Cocos) enters the lower mantle, compressive stresses increase at the eastern margin of the upper plate, and a second subduction zone (LA) is initiated. The resulting first‐order slab geometries, timings, and kinematics compare well with plate reconstructions. More generally, similar slab dynamics may provide a mechanism not only for the Caribbean reorganization but also for other tectonic episodes throughout the Americas.more » « less
-
Abstract Pacific‐Panthalassa plate tectonics are the most challenging on Earth to reconstruct during the Mesozoic and Cenozoic eras due to extensive subduction, which has resulted in large (>9,000 km length) unconstrained gaps between the Pacific and Laurasia (now NE Asia) back to the Early Jurassic. We build four contrasted NW Pacific‐Panthalassa global plate reconstructions and assimilate their velocity fields into global geodynamic models. We compare our predicted present mantle structure, synthetic geoid and dynamic topography to Earth observations. P‐wave tomographic filtering of predicted mantle structures allows for more explicit comparisons to global tomography. Plate reconstructions that include intra‐oceanic subduction in NW Pacific‐Panthalassa fit better to the observed geoid and residual topography, challenging popular models of Andean‐style subduction along East Asia. Our geodynamic models predict significant SE‐ward lateral slab advections within the NW Pacific basin lower mantle (∼2,500 km from Mesozoic times to present) that would confound “vertical slab sinking”‐style restorations of imaged slabs and past subduction zone locations.more » « less
