An official website of the United States government
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
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:
- 10533388
- Publisher / Repository:
- American Association for the Advancement of Science
- 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 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 Age-progressive seamount tracks generated by lithospheric motion over a stationary mantle plume have long been used to reconstruct absolute plate motion (APM) models. However, the basis of these models requires the plumes to move significantly slower than the overriding lithosphere. When a plume interacts with a convergent or divergent plate boundary, it is often deflected within the strong local mantle flow fields associated with such regimes. Here, we examined the age progression and geometry of the Samoa hotspot track, focusing on lava flow samples dredged from the deep flanks of seamounts in order to best reconstruct when a given seamount was overlying the mantle plume (i.e., during the shield-building stage). The Samoan seamounts display an apparent local plate velocity of 7.8 cm/yr from 0 to 9 Ma, 11.1 cm/yr from 9 to 14 Ma, and 5.6 cm/yr from 14 to 24 Ma. Current fixed and mobile hotspot Pacific APM models cannot reproduce the geometry of the Samoa seamount track if a long-term fixed hotspot location, currently beneath the active Vailulu’u Seamount, is assumed. Rather, reconstruction of the eruptive locations of the Samoan seamounts using APM models indicates that the surface expression of the plume migrated ~2° northward in the Pliocene. Large-scale mantle flow beneath the Pacific Ocean Basin cannot explain this plume migration. Instead, the best explanation is that toroidal flow fields—generated by westward migration of the Tonga Trench and associated slab rollback—have deflected the conduit northward over the past 2–3 m.y. These observations provide novel constraints on the ways in which plume-trench interactions can alter hotspot track geometries.more » « less
-
Abstract Subduction zones are fundamental features of Earth's mantle convection and plate tectonics, but mantle flow and pressure around slabs are poorly understood because of the lack of direct observational constraints on subsurface flow. To characterize the linkages between slabs and mantle flow, we integrate high‐resolution representations of Earth's lithosphere and slabs into a suite of global mantle convection models to produce physically plausible present‐day flow fields for Earth's mantle. We find that subduction zones containing wide, thick, and long slabs dominate regional mantle flow in the neighboring regions and this flow conforms to patterns predicted by simpler regional subduction models. These subduction zones, such as Kuril‐Japan‐Izu‐Bonin‐Mariana, feature prismatic poloidal flow coupled to the downgoing slab that rotates toward toroidal slab‐parallel flow near the slab edge. However, other subduction zones, such as Sumatra, deviate from this pattern because of the competing influence of other slabs or longer‐wavelength mantle flow, showing that upper mantle flow can link separate subduction zones and how flow at subduction zones is influenced by broader scale mantle flow. We find that the non‐linear dislocation creep reduces the coupling between slab motion and asthenospheric flow and increases the occurrence of non‐ideal flow, in line with inferences derived from seismological constraints on mantle anisotropy.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1029/2024GC011523
