Un continent qui s'est hissé tout seul au-dessus du seuil du froid
Imaginez un continent plat, baigné d'un climat doux, sans une seule calotte glaciaire à l'horizon. C'est l'Antarctique tel qu'il existait il y a 50 millions d'années. Aujourd'hui, il abrite la plus grande étendue de glace de la planète . Comment est-on passé de l'un à l'autre ? Et surtout : pourquoi l'Antarctique a-t-il basculé dans le froid permanent 29 millions d'années avant l'Arctique, alors que le climat terrestre était encore relativement clément ?
La réponse vient de paraître dans Science, sous la plume de Thomas Gernon et ses collègues des universités de Southampton, Durham, Potsdam, Utrecht et Florence, ainsi que du Centre GFZ Helmholtz pour les géosciences. Et elle est, pour le moins, inattendue : ce n'est pas d'abord la baisse du CO₂ atmosphérique qui a glacé l'Antarctique. C'est le sol lui-même qui s'est soulevé — lentement, inexorablement — jusqu'à franchir un seuil d'altitude critique. Un seuil géométrique, pas climatique.
Des vagues dans la roche, à 15 kilomètres par million d'années
Tout commence il y a entre 201 et 143 millions d'années, quand le supercontinent Gondwana se fragmente. L'Antarctique se sépare de l'Afrique. Cette rupture colossale n'est pas seulement une affaire de surface : elle déclenche en profondeur ce que les géodynamiciens appellent des ondes de manteau — non pas des ondes sismiques rapides, mais des instabilités convectives qui migrent lentement à l'intérieur du manteau terrestre, la couche visqueuse située entre la croûte et le noyau.
La mécanique est presque poétique dans sa lenteur. Sous la lithosphère (la plaque rigide qui forme le plancher continental), des « gouttes » de roche dense se détachent progressivement de la base du continent. En perdant ce lest profond, la plaque s'allège et remonte — comme un ballon à air chaud qui largue du sable, selon l'image proposée par Sascha Brune, responsable de la modélisation géodynamique au GFZ. Ces ondes se propagent à environ 15 à 20 kilomètres par million d'années, une vitesse établie dans une étude parue dans Nature qui a formalisé ce mécanisme pour les marges cratoniques issues de la dispersion du Gondwana.
Quand ces ondes atteignent l'Antarctique oriental, elles soulèvent le sol et font émerger les monts Gamburtsev — une chaîne de 1 200 kilomètres aujourd'hui enfouie sous 1 à 3 kilomètres de glace. Avant 50 millions d'années, la plupart de ces montagnes ne dépassaient pas 1 500 mètres. Vers 45 millions d'années, la grande majorité du plateau antarctique oriental franchit les 2 000 mètres d'altitude.
2 000 mètres : le seuil qui change tout
Pourquoi 2 000 mètres ? Parce que la température de l'air baisse d'environ 1 °C tous les 100 mètres d'altitude gagnés — un gradient bien connu des alpinistes. Franchir ce seuil, c'est entrer dans une zone où les neiges hivernales ne fondent plus complètement l'été. Les glaciers s'installent. Et une fois installés, ils déclenchent une mécanique implacable.
La glace et la neige fraîche réfléchissent une grande partie de la lumière solaire vers l'espace — c'est ce qu'on appelle l'albédo. Une surface enneigée renvoie environ 80 % du rayonnement solaire incident, contre moins de 10 % pour un sol sombre ou un océan. Quand la calotte s'étend, la région se refroidit davantage, ce qui permet à la glace de s'étendre encore plus. Comme le résume Philip Goodwin, physicien du climat à l'Université de Southampton : « À mesure que la calotte glaciaire s'étendait, sa surface brillante réfléchissait davantage de lumière solaire vers l'espace, refroidissant encore la région. » Ce refroidissement lié à l'effet glace-albédo est évalué à environ 1 °C supplémentaire selon le communiqué de l'Université de Southampton. L'air plus froid retient moins de vapeur d'eau — un autre gaz à effet de serre — ce qui amplifie encore le refroidissement. La boucle est bouclée.
C'est un exemple classique de rétroaction positive : chaque effet amplifie sa propre cause. Les mathématiciens du climat parlent d'un système qui « bascule » — il existe un état chaud stable et un état froid stable, et le passage de l'un à l'autre peut être brutal une fois le seuil franchi. Ce cadre formel est celui des modèles de bilan énergétique, dont l'analyse de stabilité a été posée par Gerald North, Robert Cahalan et James Coakley dans un survey fondateur publié dans Reviews of Geophysics : ces modèles montrent qu'une petite variation de paramètre peut faire basculer un système climatique entre deux équilibres distincts, l'un glaciaire, l'autre non.
Modéliser 100 millions d'années : un défi computationnel
Reconstituer l'évolution d'un continent sur 100 millions d'années n'est pas une affaire de conjecture. L'équipe a couplé trois familles de modèles : des simulations géodynamiques de propagation des ondes de manteau, des reconstructions topographiques de l'érosion et du soulèvement, et des modèles de calotte glaciaire. Les reconstructions du relief sont disponibles tous les 2 millions d'années, de 100 à 36 millions d'années avant le présent, dans l'archive ouverte déposée sur Zenodo par l'équipe.
Pour simuler l'érosion fluviale et l'évolution du paysage, les chercheurs ont utilisé FastScape, un algorithme dont la particularité est d'être de complexité O(n) — c'est-à-dire que le temps de calcul croît linéairement avec le nombre de points du terrain modélisé, et non pas de façon quadratique ou exponentielle. Cette propriété, décrite par Jean Braun et Sean Willett dans Geomorphology, rend possible le traitement de grilles pouvant atteindre 100 millions de nœuds. Sans cette efficacité algorithmique, simuler un continent entier sur des dizaines de millions d'années serait tout simplement hors de portée. La modélisation de la calotte elle-même repose sur PISM (Parallel Ice Sheet Model), dont la version PIK est décrite dans The Cryosphere : plutôt que de résoudre les équations complètes de Stokes pour l'écoulement de la glace — un problème d'une complexité prohibitive à grande échelle —, PISM superpose deux approximations adaptées aux objets très larges et très minces que sont les calottes.
Thea Hincks, chercheuse senior à Southampton, résume l'enjeu : « Nous avons constaté que nos modèles peuvent reproduire de façon réaliste l'évolution de l'escarpement côtier de deux kilomètres de haut, du plateau élevé et des montagnes intérieures, amorçant à terme la formation de l'inlandsis de l'Antarctique de l'Est. » Ce n'est pas une reconstruction narrative — c'est une validation quantitative, point par point, de formes de paysage observables aujourd'hui sous la glace.
Pourquoi l'Arctique a attendu 29 millions d'années de plus
L'asymétrie entre les deux pôles était l'un des grands mystères de la paléoclimatologie. L'explication classique — la baisse du CO₂ atmosphérique — ne suffisait pas : si c'était le seul facteur, les deux pôles auraient dû réagir de façon plus symétrique. Comme le formule Gernon lui-même : « Si les baisses de CO₂ étaient seules déterminantes, on s'attendrait à ce que les pôles réagissent de façon plus symétrique. »
La réponse est topographique. Les montagnes arctiques n'ont jamais atteint les altitudes des Gamburtsev. Sans seuil d'altitude suffisant, pas de nucléation pérenne de glace, même avec un CO₂ en baisse. L'Antarctique avait une longueur d'avance géologique — une longueur d'avance de 100 millions d'années, gravée dans la dynamique profonde de la Terre.
John Goodge, géologue à l'Université du Minnesota Duluth et extérieur à l'étude, juge la démonstration convaincante : « L'article défend de façon très convaincante l'idée que c'est cette réponse tectonique, se développant sur une période de plusieurs dizaines de millions d'années, qui conduit à une situation de seuil où ce relief de haute altitude atteint une élévation critique lui permettant de laisser se former la glace de façon permanente. » Il note cependant que des forages profonds sous la glace permettraient de mieux contraindre les propriétés de la lithosphère sous les Gamburtsev et d'affiner les modèles.
Ce que cette découverte change, au fond, c'est notre façon de lire les transitions climatiques majeures. Le CO₂ reste un acteur — mais il n'est pas le metteur en scène. Parfois, c'est la géologie profonde qui prépare la scène, silencieusement, sur des dizaines de millions d'années, avant que le climat n'entre en jeu. Une leçon qui n'est pas sans résonance pour comprendre les futurs « points de bascule » du système climatique actuel.
Concepts à emporter
La ligne de glace comme inconnue mathématique
L'un des aspects les plus élégants de la modélisation climatique mobilisée dans cette étude est la façon dont elle traite la frontière entre glace et sol nu — ce que les glaciologues appellent la ligne de glace — non pas comme une donnée fixe, mais comme une inconnue du problème.
Le cadre formel est celui des modèles de bilan énergétique (EBM pour energy balance model), introduits indépendamment par Mikhail Budyko et William Sellers à la fin des années 1960. Dans leur version la plus simple, ces modèles représentent la température moyenne d'une latitude x par une équation d'équilibre :
C · ∂T/∂t = Q·(1 − α(T)) − A − B·T + D·∂²T/∂x²
où Q est le flux solaire incident, α(T) l'albédo (qui saute brusquement quand la température passe sous un seuil de gel), A + B·T le rayonnement infrarouge émis, et D·∂²T/∂x² le transport de chaleur horizontal. La non-linéarité cruciale est dans α(T) : l'albédo de la glace (≈ 0,6) est bien supérieur à celui du sol nu (≈ 0,1), ce qui crée une discontinuité dans la réponse du système.
La position de la ligne de glace, notée η, est alors définie implicitement par T(η) = T*, où T* est la température de gel. Comme l'analyse Walsh et Rackauckas dans Discrete and Continuous Dynamical Systems, ce couplage entre l'équation de bilan et la ligne de glace mobile produit un système dynamique à équilibres multiples : il peut exister simultanément un état « Terre libre de glace », un état « calotte polaire stable » et un état « Terre entièrement glacée » (la fameuse Terre boule de neige). La stabilité de chacun dépend des paramètres — dont l'altitude du terrain, qui modifie localement T*.
C'est précisément ce que le soulèvement des Gamburtsev a modifié : en élevant le terrain de 500 mètres, il a abaissé la température effective locale d'environ 5 °C, faisant basculer la région d'un équilibre sans glace permanente vers un équilibre avec calotte stable. Un changement de paramètre géométrique qui a réécrit l'histoire climatique d'un continent entier.
