Une approche numérique pour comprendre les vitesses d'englacement du Quaternaire

Au cours du Quaternaire (depuis 2.6 Ma), les calottes glaciaires connaissent différents épisodes d'avancée-retrait correspondant aux cycles glaciaires-interglaciaires. L'étude de ces événements épisodiques permet de mieux comprendre les mécanismes à l'origine de l'évolution de la Terre et d'améliorer les prévisions dans le contexte du réchauffement climatique actuel.La simulation des interactions entre les calottes polaires et le climat sur des échelles de temps aussi longues nécessite des approches de modélisation numérique qui représentent suffisamment le système réel tout en maintenant des coûts de calcul faibles. Dans la première partie de cette thèse, j'utilise un modèle système terre système terrestre de complexité intermédiaire (iLOVECLIM) couplé au modèle 3D de calotte polaires GRISLI pour simuler l'avancée abrupte de la calotte glaciaire au début du dernier cycle glaciaire (120-115 kaBP). Les résultats indiquent que les débuts de glaciation ne peuvent pas être expliqués uniquement par la théorie astronomique (en résponse aux forçages orbitaux). Les rôles de la biosphère et de l'océan par le biais de différents mécanismes de rétroaction doivent être inclus pour expliquer la localisation et l'étendue de l'avancée de la calotte glaciaire. De plus, une simulation appropriée du processus d'accumulation de la calotte glaciaire est essentielle pour obtenir des résultats corrects.Dans la deuxième partie de la thèse, j'étudie les comportements d'un modèle de neige multicouche BESSI afin de fournir une simulation de bilan de masse de surface (SMB) davantage basée sur la physique pour iLOVECLIM-GRISLI. Le modèle de neige présente de bons résultats par rapport à un modèle climatique régional MAR de pointe pour le climat actuel dans différentes conditions de calotte glaciaire. Pour le dernier interglaciaire (130-116 kaBP), BESSI forcé par iLOVECLIM montre une plus grande sensibilité aux forçages climatiques que la paramétrisation SMB existante d'iLOVECLIM-GRISLI. En outre, l'évolution du SMB simulée par BESSI-iLOVECLIM se situe également dans une fourchette acceptable par rapport à d'autres études. Cependant, comme ce modèle de neige est davantage fondé sur la physique que la paramétrisation existante, l'influence des biais d'iLOVECLIM est plus importante pour BESSI, ce qui nuit à ses performances. Moyennant des travaux à venir sur la correction de biais et la méthode de couplage, mon étude ouvre la voie à l'utilisation de BESSI dans le cadre du couplage entre le modèle de climat iLOVECLIM et le modèle de calottes glaciaires GRISLI.