ALBERT

Description: On Earth, oceanic plates subduct beneath a variety of overriding plate (OP) styles, from relatively thin and negatively buoyant oceanic OPs to thick and neutrally/positively buoyant continental OPs. The inclusion of an OP in numerical models of self-consistent subduction has been shown to reduce the rate that subducting slabs roll back relative to the equivalent single plate models. We use dynamic, 2-D subduction models to investigate how the mechanical properties, namely viscosity, thickness, and density, of the OP modify the slab rollback rate and the state of stress of the OP. In addition, we examine the role of the subducting plate (SP) viscosity. Because OP deformation accommodates the difference between the slab rollback rate and the far-field OP velocity, we find that the temporal variations in the rollback rate results in temporal variations in OP stress. The slabs in our models roll back rapidly until they reach the lower mantle viscosity increase, at which point the rollback velocity decreases. Concurrent with this reduction in rollback rate is a switch from an OP dominated by extensional stresses to a compressional OP. As in single plate models, the viscosity of the SP exerts a strong control on subducting slab kinematics; weaker slabs exhibit elevated sinking velocities and rollback rates. The SP viscosity also exerts a strong control on the OP stress regime. Weak slabs, either due to reduced bulk viscosity or stress-dependent weakening, have compressional OPs, while strong slabs have dominantly extensional OPs. While varying the viscosity of the OP alone does not substantially affect the OP stress state, we find that the OP thickness and buoyancy plays a substantial role in dictating the rate of slab rollback and OP stress state. Models with thick and/or negatively buoyant OPs have reduced rollback rates, and increased slab dip angles, relative to slabs with thin and/or positively buoyant OPs. Such elevated trench rollback for models with positively buoyant OPs induces extensional stresses in the OP, while OPs that are strongly negatively buoyant are under compression. While rollback is driven by the negative buoyancy of the subducting slab in such models of free subduction, we conclude that the physical properties of the OP potentially play a significant role in modulating both rollback rates and OP deformation style on Earth.

Description: The subduction of lithospheric plates is partitioned between subducting plate motion and lateral slab migration (i.e. trench retreat and advance). We use 3-D, dynamic models of subduction to address the role of a power-law mantle viscosity on subduction dynamics and, in particular, rates of trench retreat. For all numerical models tested, we find that a power-law rheology results in reduced rates of trench retreat, and elevated slab dip angles, relative to the equivalent isoviscous mantle model. We analyse the asthenospheric pressure distribution and the style of mantle flow, which exhibits only limited variability as a function of mantle rheology, in order to compute estimates of the mantle forces associated with subduction. The inclusion of a power-law rheology reduces the mantle shear force (which resists subducting plate motion) to a greater degree than it reduces the dynamic pressure gradient across the slab (which resists trench retreat). Therefore, the inclusion of a power-law mantle rheology favours a shift towards a subduction mode with a reduced trench retreat component, typically a relative reduction of order 25 per cent in our 3-D models. We suggest that this mechanism may be of importance for reducing the high trench retreat rates observed in many previous models to levels more in line with the average subduction partitioning observed on Earth at present (i.e. trench velocity ≤ plate velocity), for most absolute plate motion reference frames.