Trench migration, net rotation and slab–mantle coupling

Laboratory models have been conducted to improve our understanding of the role that the resistance of the slab to bending and its coupling to the ambient mantle play in subduction dynamics over geological time scales. Our models are set up with a viscous plate of silicone (lithosphere) subducting under negative buoyancy in a viscous layer of glucose syrup (mantle). For our study, the lithosphere/upper mantle viscosity contrast has been systematically varied, from ~ 10 to ~ 105 in order to explore the parameter space between weak and strong slab dynamics. We found that subduction is characterized by a retreating mode for viscosity ratios > 104, by the coexistence of a retreating mode and an advancing mode for viscosity ratios between ~ 104 and ~ 102, and quasi-stationary, Rayleigh–Taylor like behaviour for ratios < 102. By combining our experimental results and kinematic data from current subduction zones in four reference frames which differ in the amount of net rotation, we infer that a lithosphere/upper mantle viscosity contrast of 150–500 is necessary to obtain realistic trench/subducting plate velocity ratios as well as the variability of subduction styles observed in nature.