Asthenospheric flow-driven lithospheric deformation in analogue models – a novel methodological approach and implications for natural systems

In subduction systems, asthenospheric flow, generated by subducting slabs, is considered as one of the key forces contributing to the deformation of the overlying lithosphere. Previous analogue modelling studies predominantly focused on understanding the kinematics and dynamics of subduction roll-back-driven asthenospheric flow, without looking at the influence of that flow on upper-plate deformation due to the modelling setups or methodological limitations. We developed a novel analogue modelling approach where gravity-driven asthenospheric flow represents the main driver for upper plate deformation. Volume-constant flow within the deformation box is achieved by an inlet-outlet system. In the models, we gradually increase the setup complexity from single-layer asthenosphere-only models to 4-layer asthenosphere-lithosphere models to test flow velocity distribution and its sensitivity to the outlet size, model thickness and rheological stratification of the model, as well as the transfer of deformation from the asthenosphere to the overlying lithosphere. Furthermore, we study the effects of the inherited lithospheric structures, such as weak zones representing old sutures, on deformation transfer. The results are compared with the Pannonian-Carpathians system of south-eastern Europe, where the large Pannonian back-arc basin formed during the Miocene retreat of the Carpathians slab.

For the methodological approach, the results show that asthenospheric flow can be fully controlled by the inlet-outlet system by adjusting the outlet size, which provides an efficient mechanism for the deformation of the overlying mechanically stratified lithosphere. The models also demonstrate that the back-arc extension is initiated farther away from the asthenospheric flow origin (i.e., the outlet in the models or slab-roll back in nature). The subsequent deformation propagates in two directions, towards the flow origin, and farther away from it, both directions controlled by the shape of an indenter located laterally to the subduction zone. Most of the back-arc extension and the lithospheric thinning are accommodated in the area farther to the “slab” due to the strain shadow effect of the indenter. The indenter also contributes significantly to the strain partitioning in its closer proximity where a complex pattern of bi-directional extension, transtensional, strike-slip and transpressional deformation forms. The weak zones accommodate the onset of back-arc extension or act as transfer zones between areas with different extension rates, depending on their orientation relative to the asthenospheric flow. These models show several similarities with the Pannonian-Carpathians system, where most of the Pannonian lithospheric thinning is located at a significant distance from the subducting Carpathians slab, bypassing the Transylvanian-Apuseni area. This extension started by reactivation of the Neotethys suture zone, while the Mid-Hungarian Fault zone transferred the deformation between areas of higher extension to the south and lower extension to the north. Furthermore, several triangular-shaped sub-basins within and at the margin of the Pannonian Basin are radially located around the Moesian NW corner, similar to our modelling results. The complex pattern of the bi-directional extension and strike-slip observed in the models were recorded by the Carpathians-Balkanides orocline in the vicinity of the Moesian indenter.