Dynamics of double-polarity subduction: application to the Western Mediterranean

The evolution of the Western Mediterranean is a highly debated question by geologists and geophysicists. Even though most scientists agree in considering slab roll-back to be the driving mechanism of the tectonic evolution of this area, there is still no consensus about the initial setup and its time evolution. A recent model proposed by Vergés and Fernàndez (2012) suggests a lateral change in subduction polarity of the Ligurian-Thetys oceanic domain to explain the formation and evolution of the Betic-Rif orogenic system and the associated Alboran back-arc basin. Such geodynamic scenario is also proposed for different converging regions. The aim of this study is to analyze the dynamic evolution of a double-polarity subduction process and its consequences in order to test the physical feasibility of this interaction and provide geometries and evolutions comparable to those proposed for the Western Mediterranean. The 3D numerical model of double-polarity subduction is carried out via the Underworld framework. Tectonic plate behavior is described by equations of fluid dynamics in the presence of several different phases. Underworld solves a non-linear Stokes flow problem using Finite Elements combined with particle-in-cell approach, thus the discretization combines a standard Eulerian Finite Element mesh with Lagrangian particles to track the location of the phases. The final model consists of two oceanic plates with viscoplastic rheology subducting into the upper mantle and the problem is driven by Rayleigh-Taylor instability. The main factors to be studied are the interaction between the two plates, the poloidal and toroidal mantle fluxes, the velocity variations of slabs, the stress distribution and the variations in the trench morphology.