Crust-mantle delamination enables continental subduction and flake tectonics: insights from numerical modelling

Continental collision occurs when two continents are dragged towards each other by the pull of the attached subducting oceanic lithosphere. Previous geodynamic modeling studies of collisional systems focused on first-order processes (such as coupled/decoupled regimes, continental delamination, slab break-off dynamics) and regional or even local scale dynamics (e.g., exhumation of HP/UHP rocks, surface topography). However, continuous subduction of continental lithospheric mantle after the onset of collision and long-term dynamics of continental subduction remains poorly constrained. Long-term continental subduction bears major geodynamic implications for the evolution of past and present collision zones.

To this aim, we use the geodynamic code LaMEM to perform high-resolution (2048 × 512) 2D buoyancy-driven numerical models, coupled with phase diagrams to account for density changes, of continued continental subduction with conditions that favor flake tectonics. We investigate the role of lower crust rheology to assess which rheological scenarios allow continental flaking and, thus, continued subduction of continental lithospheric mantle.

Our preliminary results exhibit long-term continental subduction, due to decoupling of the lower crust from the subducting continental mantle and/or density changes. This separation allows the deformation to be transmitted onto the overriding plate, with the emplacement of the subducting plate crust onto the overriding plate spanning more than 350 km and lasting over 100 Myr.

This work was funded by the Portuguese Fundação para a Ciência e a Tecnologia (FCT) I.P./MCTES through national funds (PIDDAC) – UIDB/50019/2020 (https://doi.org/10.54499/UIDB/50019/2020), UIDP/50019/2020 (https://doi.org/10.54499/UIDP/50019/2020) and LA/P/0068/2020 (https://doi.org/10.54499/LA/P/0068/2020), and through scholarship UI/BD/154679/2023.