Sediment buoyancy controls the effective slab pull force and deviatoric stress along trenches: Insights from 3D free-subduction model

Abstract

In the complicated force equilibrium of subduction dynamics, the slab pull resulting from the negative buoyancy of the lithosphere is considered the primary driving force of plate tectonics. In contrast, the isostatic restoring force due to lighter materials covering the oceanic crust (e.g., sediment) has been recognized as a relatively minor contributor. However, as the restoring force counteracts the gravitational pull, the resultant interaction has the potential to control the force balance governing the subduction system. Here, we conducted a series of three-dimensional viscoelastic-free subduction simulations with varying sediment distributions, thicknesses, and densities. We found that a higher restoring force suppressed the trench retreat and decreased the trench velocity. The trench curvatures increased with a high differential restoring force along the trench strike over time. Our results consistently showed a negative correlation between the effective slab pull and the isostatic restoring forces. Moreover, the magnitude and distribution of the isostatic restoring force along the trench strike determined the deviatoric stress level and location of the stress concentration within the subducting slab. We also found that locally abundant sediment margins may lead to a high deviatoric stress in the trench, which is associated with the occurrence of large earthquakes. We argue that sediment buoyancy can significantly affect the subduction kinematics and dynamics to a greater degree than previously suggested. Our findings provide fresh insights into the development of a global predicted velocity model, that reflects the effect of the isostatic restoring force, and provides a better understanding of subduction earthquakes.