Ré-examen du processus d'initiation d'une zone de subduction dans les modèles numériques d'intéractions lithosphère-asthénosphère // Revisiting subduction initiation in numerical models of plate and mantle dynamics

Pour qu'une nouvelle zone de subduction se mette en place, une frontière de plaque de faible résistance doit réussir à se former, par la localisation de la déformation sur une (mince) interface de découplage entre deux plaques lithosphériques, plaques dont le 'coeur' peut-être très résistant résistant. Lallemand et Arcay (2021) ont répertorié des cas naturels d'initiation de subduction au cours du Cénozoïque, Arcay et al. (2023) montrant que le succès de l'initiation au niveau d'une faille transformante océanique dépendait de la profondeur de la zone faible interloque (où la convergence est susceptible de se localiser) par rapport à la profondeur où la résistance au cisaillement était maximale au sein de la plaque chevauchante. Cette profondeur critique peut être estimée dans certains cas (suivant le modèle rhéologique choisi) par la profondeur de transition entre les mécanismes de déformations fragile et ductile. Cette transition dépend de la paramétrisation rhéologique du manteau lithosphérique froid, en approximant la déformation fragile par une contrainte seuil dépendante de la profondeur et la déformation ductile par une loi d'écoulement par fluage - dislocation. Par ailleurs, des données récentes sur la dynamique des dislocations dans l'olivine à basse température (expérimentation, modélisation) prédisent une résistance maximale plus faible au voisinage de la transition fragile-ductile, ce qui pourrait favoriser la localisation de la déformation par des rétroactions non-Newtoniennes lors de l'enfoncement de la plaque ou de l'extension de la plaque (Garel+ 2020, Van Broeck+ 2026). Cependant, ceci n'a encore jamais été exploré dans un contexte d'initiation de subduction.

Ce sujet de thèse vise à évaluer comment l'implémentation du fluage- dislocation à basse température influe sur les scénarios d'initiation et de maturation d'une zone de subduction naissante. Différentes paramétrisations rhéologiques (domaines cassant, ductile, transition fragile-ductile) seront comparées dans des modèles numériques soit régionaux, soit à l'échelle plus vaste de la convection mantellique. Nous étudierons en particulier :

- si le fluage-dislocation basse température permet l'initiation d'une zone de subduction même pour une contrainte seuil plus élevée (résistance cassante)

- comment le partitionnement de la déformation, l'échauffement par cisaillement, la force tectonique compressive, la résistance à la flexion, l'affaiblissement/résistance au cisaillement de l'asthénosphère, la traction de la proto-plaque plongeante, etc. évoluent au cours des différentes étapes de du développement d'une zone de subduction

- comment, selon l'équilibre des forces, la subduction peut (ou non) devenir auto-entretenue après son initiation

- si la déformation de surface, résultant de la dynamique profonde de la plaque ou de l'évolution du découplage mécanique de l'interface, peut initier une nouvelle zone de subduction.

Les modèles régionaux permettront de comprendre les rétroactions entre les dépendances rhéologiques, le champ de déformation et la dynamique des plaques. Les modèles de convection mantellique à grande échelle permettront quant à eux d'explorer dans un cadre auto-régulé et dépendant du temps, par exemple comment une subduction peut affecter l'initiation d'une autre zone de subduction dans son voisinage. Les résultats seront comparés aux enregistrements géologiques: mode d'accommodation de la déformation, contexte géodynamique, échelle de temps de la localisation, et évolution des paramètres de subduction (vitesses de convergence, âge des plaques convergentes). Enfin, les expériences numériques seront analysées pour quantifier les rétroactions thermiques et mécaniques impliquées dans les différentes étapes du processus d'initiation de la subduction. Nous chercherons à établir des lois de comportement pertinentes (principaux paramètres de contrôle, évolution de l'équilibre des forces au cours du développement d'une zone de subduction).
------------------------------------------------------------------------------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------------------------------------------------------------------------------

During the initiation of a subduction zone, a weak plate boundary is expected to form by the localization of deformation on a decoupling interface between two strong lithosphere cores. Lallemand and Arcay (2021) have documented natural cases of subduction initiations during the Cenozoic, with Arcay+ (2023) showing that successful initiation at an oceanic transform fault depends on the weak zone depth (where convergence is likely to localize) with respect to the depth where the resistance to shear is maximum within the overriding plate. This critical depth can be approximated in some cases by the brittle-ductile transition depth. This transition depends on the rheological parameterization of the cold lithospheric mantle, approximating brittle deformation by a depth-dependent yield strength and ductile deformation by a dislocation creep flow law. Recent data on low-temperature dislocation dynamics predict a smaller peak resistance at the brittle-ductile transition, which favor deformation localization through non-Newtonian feedbacks during slab sinking or plate extension (Garel et al. 2020, Van Broeck et al. 2026). However to our knowledge the possible role of a low-temperature plasticity has never been investigated in the process of subduction initiation.

The present PhD subject aims at evaluating how the implementation of low-temperature dislocation creep impacts scenarios of subduction initiation and evolution. Different rheological parameterizations of yield-strength and dislocation creep will be compared in regional and large-scale numerical models of mantle dynamics. We will in particular investigate :

- if low-temperature creep results in subduction initiation even for larger yield-stress (brittle resistance)

- how deformation partitionning, shear heating, tectonic force, bending resistance, asthenosphere weakening/drag, slab pull, etc. evolve during the different stages of subduction initiation

- how subduction may (or not) be self-sustained after the first initial stages depending on the force balance

- if surface deformation, resulting from deep slab dynamics or variation in interface decoupling, can trigger the initiation of a new subduction zone.

Regional models will permit an understanding of feedbacks between rheological dependencies, deformation pattern and plate dynamics. Large-scale mantle convection models will on the other hand allow to explore self-consistent and time-dependent thermo-mechanical fields, e.g. how a neighbouring subduction may affect another subduction initiation. The simulation results will be compared to geological records, regarding the mode of accommodation of deformations, geological settings, timescale of localization and the evolution of subduction parameters (convergence speeds, age of the converging plates). The numerical experiments will be analyzed to decipher the thermal and mechanical feedbacks involved in the different stages of the process of subduction initiation. We expect to derive some relevant scaling laws (main controlling parameters, evolution of force balance as a subduction zone develops) to highlight the key ingredients of a successful subduction zone initiation.
------------------------------------------------------------------------------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------------------------------------------------------------------------------

Début de la thèse : 01/10/2026
WEB : https://anr.fr/en/funded-projects-and-impact/funded-projects/project/funded/project/b2d9d3668f92a3b9fbbf7866072501ef-6acf14cabe/?tx_anrprojects_funded%5Bcontroller%5D=Funded&cHash=13630af98714fb12ed2a

Concours GAIA

Université de Montpellier

584 GAIA - Biodiversité, Agriculture, Alimentation, Environnement, Terre, Eau

- Master 2 en sciences de la Terre, sciences des matériaux ou physique - Connaissances en physique (mécanique des fluides, transferts thermiques) - Connaissances en traitement de données et programmation (par exemple Python, script Linux) - Motivé, ouvert d'esprit et doté d'un esprit critique - Esprit d'équipe et bonnes aptitudes à la communication - Bonnes aptitudes à la communication écrite et orale en anglais
- Master 2 in Earth Sciences, Material Science or Physics - Knowledge in physics (fluid mechanics, heat transfer) - Knowledge in data processing and programmation (e.g. Python) - Motivated, open-minded and with critical thinking - Teamwork and communication skills - Written and oral communication skills in English - Scientific interest and curiosity