Oblique subduction and mantle flow control on upper plate deformation: 3D geodynamic modeling

Most subduction zones on Earth are oblique, i.e., the angle between the plate convergence vector and the trench notably differs from 90°. Therefore, modeling and understanding the strain partitioning in the forearc, the development of extensional basins in the back-arc region and the diachronous transition from subduction to collision require a 3D approach. Here, we assess how oblique oceanic subduction and subsequent collision and associated mantle flow around the subducted lithosphere control the thermo-mechanical evolution of active margins. We conducted a series of 3D thermo-mechanical subduction models and discuss the influence of different subduction obliquity angles, the role of mantle flow variations and their connection with sediment transport and back-arc deformation. Numerical models are complemented by scaled analogue models to visualize the mantle flow evolution. Oceanic subduction along an oblique trench results in asymmetric mantle return flow leading to the gradual decrease of the subduction obliquity angle driven by the gradual rotation of the lower plate and the along-trench variation of slab retreat. This creates laterally variable subduction velocities and slab geometries. Back-arc extension is governed by both the toroidal mantle flow along the slab edges and by the oblique subduction induced lateral mantle flow gradient. The diachronous transition from oceanic to continental subduction and collision facilitates the laterally variable trench advance and retreat and back-arc deformation. Tectonically induced lateral sediment transport in the trench and along the subduction interface decreases its strength and viscosity and can alter subduction velocities. Our model results provide critical insights into the evolution of oblique subduction and collisional systems, such as the Arabia-Eurasia convergence zone.