Influence of trench width on subduction hinge retreat rates in 3‐D models of slab rollback
Top Cited Papers
- 29 March 2006
- journal article
- Published by American Geophysical Union (AGU) in Geochemistry, Geophysics, Geosystems
- Vol. 7 (3)
- https://doi.org/10.1029/2005gc001056
Abstract
Subduction of tectonic plates limited in lateral extent and with a free‐trailing tail, i.e., “free subduction,” is modeled in a three‐dimensional (3‐D) geometry. The models use a nonlinear viscoplastic rheology for the subducting plate and exhibit a wide range of behaviors depending on such plate characteristics as strength, width, and thickness. We investigate the time evolution of this progressive rollback subduction, measure the accompanying return flow in the upper mantle, and quantify the plate kinematics. Due to the 3‐D geometry, flow is allowed to accompany slab rollback around the lateral edges of the slab (the toroidal component), as opposed to 2‐D geometry, where material is forced to flow underneath the slab tip (the poloidal component). A simple force balance is provided which relates the speed of backward trench migration to the resistive forces of generating flow and weakening the plate. Our results indicate most of the gravitational energy of the system (i.e., the negative buoyancy of the subducting slab) is converted into a toroidal flow (∼69%), a much smaller amount goes into weakening the plate (∼18%), and the remaining amount goes into driving flow parallel to displacement of the slab (∼13%). For the trench widths (W) we investigate (≤1500 km), a maximum trench retreat rate occurs for trenches 600 km wide, which is attributed to the interaction between a plate of finite width and the induced flow (which has a lengthscale in the horizontal direction). These numerical results quantitatively agree with comparable 3‐D laboratory experiments using analogue models with a purely viscous plate material (Schellart, 2004a, 2004b), including correlations between increasing retreat rate with increasing plate thickness, trench width for maximum retreat rate (500 km), and estimated amount of slab buoyancy used to drive rollback‐induced flow (∼70%). Several implications for plate tectonics on Earth result from these models such as rollback subduction providing a physical mechanism for ephemeral slab graveyards situated above the more viscous lower mantle (and endothermic phase transition) prior to a flushing event into the lower mantle (mantle avalanche).Keywords
This publication has 62 references indexed in Scilit:
- Plate motions, slab dynamics and back-arc deformationPhysics of the Earth and Planetary Interiors, 2005
- Evolving force balance during incipient subductionGeochemistry, Geophysics, Geosystems, 2004
- Kinematics of subduction and subduction‐induced flow in the upper mantleJournal of Geophysical Research, 2004
- Variability in flow and temperatures within mantle subduction zonesGeochemistry, Geophysics, Geosystems, 2004
- On the role of subducting oceanic plateaus in the development of shallow flat subductionTectonophysics, 2002
- Predicting plate velocities with mantle circulation modelsGeochemistry, Geophysics, Geosystems, 2001
- Geodynamic models of deep subductionPhysics of the Earth and Planetary Interiors, 2001
- The influence of olivine metastability on deep subduction of oceanic lithospherePhysics of the Earth and Planetary Interiors, 2000
- Subduction influence on magma supply at the East Scotia RidgeEarth and Planetary Science Letters, 1997
- The accuracy of finite element solutions of Stokes's flow with strongly varying viscosityPhysics of the Earth and Planetary Interiors, 1996