High-resolution numerical modeling of stress distribution in visco-elasto-plastic subducting slabs

We employ high-resolution 2D thermomechanical modelling with complex elasto-visco-plastic rheology to analyse stresses in subducting and overriding plates. The model contains a dynamic subduction channel, defined as a few km thick zone of specific (weak) visco-plastic rheology between the slab and the overriding plate. Our high-resolution model corresponds to the zoom-in of the large-scale thermomechanical model of the Central Andes, with the 4-fold resolution increase. Initial model configuration (material, temperature and stress distributions) was interpolated from the parent model. Boundary conditions including slab pull and push velocities as well as velocity of the overriding plate, were also re-mapped from the large scale model. We demonstrate that, for shallow-dipping overridden slabs, the dominant stress is unbending stress at 50–100 km depth. This stress reaches 1–1.5 GPa and strongly prevails over the slab bending stress near the trench, which is limited by the frictional plastic yield. Unbending creates double (5–10 km distance between maximums) compression zone in the upper part of the slab and a zone of extension some 20–30 km below. Pressure in these zones differs from the lithostatic pressure to as much as plus/minus 0.8 GPa. Details of the distribution and magnitude of bending and unbending stresses in the slab are mainly controlled by the overriding velocity and direction of the slab-pull force. Overpressure in the subduction channel in all models does not exceed 100 MPa. Large over- and under-pressures in the slab, do not, however, lead to significant (more than 10 km) offsets of main dehydration reactions in crust and mantle of the slab.