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
The effects of plate interface rheology on subduction kinematics and dynamics
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
- PAR ID:
- 10340949
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 230
- Issue:
- 2
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- 796 to 812
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Simulating present‐day solid Earth deformation and volatile cycling requires integrating diverse geophysical data sets and advanced numerical techniques to model complex multiphysics processes at high resolutions. Subduction zone modeling is particularly challenging due to the large geographic extent, localized deformation zones, and the strong feedbacks between reactive fluid transport and solid deformation. Here, we develop new workflows for simulating 3‐dimensional thermal‐mechanical subduction and patterns of volatile dehydration at convergent margins, adaptable to include reactive fluid transport. We apply these workflows to the Hikurangi margin, where recent geophysical investigations have offered unprecedented insight into the structure and processes coupling fluid transport and solid deformation across broad spatiotemporal scales. Geophysical data sets constraining the downgoing and overriding plate structure are collated with the Geodynamic World Builder, which provides the initial conditions for forward simulations using the open‐source geodynamic modeling software code ASPECT. We systematically examine how plate interface rheology and hydration of the downgoing plate and upper mantle influence Pacific–Australian convergence and seismic anisotropy. Models prescribing a plate boundary viscosity of Pa s best reproduce observed plate velocities, and changing the configuration of the Pacific–Australia plate boundary directly influences the modeled plate motions. Models considering hydrated olivine fabrics best reproduce observations of seismic anisotropy. Predicted patterns of slab dehydration and mantle melting correlate well with observations of seismic attenuation and arc volcanism. These results suggest that hydration‐related rheological heterogeneity and related fluid weakening may strongly influence slab dynamics. Future investigations integrating coupled fluid transport and global mantle flow will provide insight into the feedbacks between subduction dynamics, fluid pathways, and arc volcanism.more » « less
-
New Zealand's Hikurangi margin is known for recurring shallow slow slip, numerous forearc seeps, and a productive volcanic arc. Fluids derived from the subducting slab are implicated in these processes. However, prior studies lacked evidence of basic crustal structure of the slab, or of its water content that would allow an assessment of fluid budgets. We review several recent studies that place bounds on the fluid reservoirs within the subducting Hikurangi Plateau that could be released between the forearc and backarc regions. Subducting sediments are thickest (> 1 km) in the southern Hikurangi margin, where there is a unit of turbidites beneath the regional proto decollement. These subducting sediments begin draining near the deformation front, resulting in a 20-30 % loss of volumetric fluid content. In contrast, the central and northern Hikurangi margins lack a continuous unit of subducting sediment. Here, lenses of poorly drained sediment underthrust the forearc in the wakes of seamount collisions. The Hikurangi Plateau's crustal structure resembles normal oceanic crust with a doubled upper crust of basalt and diabase. Above this upper crust is a ~1.5 km thick unit of hydrated volcaniclastic conglomerates. Seamounts can locally increase the upper crust's thickness by an extra ~1-3 km, raising the amount of porous, altered volcanic material. Finally, P-wave velocity models of the slab's upper mantle show velocity changes that could indicate moderate differences in serpentinization. Active bend-faults that could circulate fluids to the upper mantle are sparse prior to subduction. However, upon subduction the upper mantle seismic velocities of the Hikurangi Plateau are significantly less in the north compared to the south, possibly due to enhanced slab faulting beneath the forearc. Separate thermo-petrologic models for the shallow forearc and deeper subduction system suggests that fluid release from volcaniclastic units and the thickened Hikurangi Plateau upper crust is expected to occur over a range of depths extending from ~12 km to ~130 km, providing fluids for onshore seep systems and hydrous melting of the mantle wedge, whereas dehydration of serpentinite is greatest beyond the arc front. Subducting sediments and volcaniclastic units are the most readily available source of fluids for shallow slow slip.more » « less
-
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
-
Variable viscosity in Earth’s mantle exerts a fundamental control on mantle convection and plate tectonics, yet rigorously constraining the underlying parameters has remained a challenge. Inverse methods have not been sufficiently robust to handle the severe viscosity gradients and nonlinearities (arising from dislocation creep and plastic failure) while simultaneously resolving the megathrust and bending slabs globally. Using global plate motions as constraints, we overcome these challenges by combining a scalable nonlinear Stokes solver that resolves the key tectonic features with an adjoint-based Bayesian approach. Assuming plate cooling, variations in the thickness of continental lithosphere, slabs, and broad scale lower mantle structure as well as a constant grain size through the bulk of the upper mantle, a good fit to global plate motions is found with a nonlinear upper mantle stress exponent of 2.43 0.25 (mean SD). A relatively low yield stress of 151 19 MPa is required for slabs to bend during subduction and transmit a slab pull that generates asymmetrical subduction. The recovered long-term strength of megathrusts (plate interfaces) varies between different subduction zones, with South America having a larger strength and Vanuatu and Central America having lower values with important implications for the stresses driving megathrust earthquakes.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/gji/ggac075
