Plan
Comptes Rendus

Géophysique externe, climat et environnement
Insolation, CO2 et précipitations en période interglaciaire
Comptes Rendus. Géoscience, Volume 337 (2005) no. 1-2, pp. 69-78.

Résumés

Des modèles climatiques de complexités différentes ont prouvé qu'une partie importante de la variation climatique à long terme (de la dizaine à la centaine de milliers d'années) pouvait s'expliquer par les variations des paramètres orbitaux aux mêmes échelles de temps. Un de ces modèles climatiques de complexité intermédiaire (LLN 2-D NH) [Gallée et al., J. Geophys. Res. 96 (1991) 13161] a permis de simuler les variations climatiques à long terme, dont la succession des cycles glaciaires–interglaciaires durant le Quaternaire. Ayant démontré ses capacités à simuler les climats du passé, le modèle a été utilisé pour simuler le climat futur. Dans la plupart des hypothèses, le modèle a suggéré la prolongation de notre interglaciaire sur près de 50 000 ans. Un modèle atmosphérique régional (MAR) [Gallée, Mon. Weather Rev. 123 (7) (1995) 2051–2069 ; Gallée, Schayes, Bound-Lay. Meteorol. 59 (1994) 141–161] a ensuite permis de préciser les implications locales des changements climatiques simulés par LLN 2-D NH, notamment pour une augmentation de la concentration en CO2 dans l'atmosphère et pour un changement du forçage astronomique.

Climate models of different complexity confirmed that the long-term variations of the astronomical parameters are driving the climatic variations during the Quaternary (and other geological times) at the timescale of tens to hundreds of thousands of year. An Earth system Model of Intermediate Complexity (LLN 2-D NH) [Gallée et al., J. Geophys. Res. 96 (1991) 13161] was able to simulate long-term climatic variations such as the recurrence of glacial interglacial cycles during the Quaternary. This model is used here to simulate the future climate. In most of the hypotheses, the model simulates an extended present interglacial covering the next 50 000 years. A regional atmospheric model (MAR) [Gallée, Mon. Weather Rev. 123 (7) (1995) 2051–2069; Gallée, Schayes, Bound-Lay. Meteorol. 59 (1994) 141–161] allows the simulation of the climatic changes at a more local scale, in particular in the case of an increase in atmospheric CO2 concentration and of a change in the astronomical forcing.

Métadonnées
Reçu le :
Accepté le :
Publié le :
DOI : 10.1016/j.crte.2004.07.009
Mot clés : Climat futur, Modèle climatique, Modèle climatique de complexité intermédiaire (EMIC), Modèle atmosphérique régional
Keywords: Future climate, Climate model, Earth system Model of Intermediate Complexity (EMIC), Regional atmospheric model

Marie-France Loutre 1 ; André Berger 1

1 Institut d'astronomie et de géophysique Georges Lemaître, université catholique de Louvain, 2, chemin du Cyclotron, B-1348 Louvain-la-Neuve, Belgique
@article{CRGEOS_2005__337_1-2_69_0,
     author = {Marie-France Loutre and Andr\'e Berger},
     title = {Insolation, {CO\protect\textsubscript{2}} et pr\'ecipitations en p\'eriode interglaciaire},
     journal = {Comptes Rendus. G\'eoscience},
     pages = {69--78},
     publisher = {Elsevier},
     volume = {337},
     number = {1-2},
     year = {2005},
     doi = {10.1016/j.crte.2004.07.009},
     language = {fr},
}
TY  - JOUR
AU  - Marie-France Loutre
AU  - André Berger
TI  - Insolation, CO2 et précipitations en période interglaciaire
JO  - Comptes Rendus. Géoscience
PY  - 2005
SP  - 69
EP  - 78
VL  - 337
IS  - 1-2
PB  - Elsevier
DO  - 10.1016/j.crte.2004.07.009
LA  - fr
ID  - CRGEOS_2005__337_1-2_69_0
ER  - 
%0 Journal Article
%A Marie-France Loutre
%A André Berger
%T Insolation, CO2 et précipitations en période interglaciaire
%J Comptes Rendus. Géoscience
%D 2005
%P 69-78
%V 337
%N 1-2
%I Elsevier
%R 10.1016/j.crte.2004.07.009
%G fr
%F CRGEOS_2005__337_1-2_69_0
Marie-France Loutre; André Berger. Insolation, CO2 et précipitations en période interglaciaire. Comptes Rendus. Géoscience, Volume 337 (2005) no. 1-2, pp. 69-78. doi : 10.1016/j.crte.2004.07.009. https://comptes-rendus.academie-sciences.fr/geoscience/articles/10.1016/j.crte.2004.07.009/

Version originale du texte intégral

1 Introduction

Le climat a varié dans le passé, de manière continue et à différentes échelles de temps. Les glaciations se sont intensifiées au cours du Quaternaire moyen et supérieur. Plusieurs phases glaciaires ont pu être identifiées, notamment à partir de l'analyse des sédiments marins, des calottes glaciaires et des enregistrements continentaux. Une analogie entre ces enregistrements à l'échelle mondiale a confirmé qu'il s'agissait bien de changements climatiques globaux et synchrones. Le dernier de ces cycles glaciaires–interglaciaires a commencé avec l'Éémien, il y a quelque 125 000 ans. Cet interglaciaire a été suivi d'une glaciation progressive, qui a culminé avec le dernier maximum glaciaire, vers 21 ka (ka : millier d'années) BP (BP : Before Present ; avant le présent). Une déglaciation rapide a ensuite fait place à l'interglaciaire actuel, l'Holocène. Ce dernier a commencé il y a 10 000 ans et l'optimum climatique fut atteint il y a environ 6000 ans. Le climat du dernier maximum glaciaire était fondamentalement différent du climat actuel. Des inlandsis de 2 à 3 km d'épaisseur couvraient une partie importante des continents dans les hautes latitudes nord. L'étendue de la banquise était largement plus grande qu'actuellement. Le niveau des mers était inférieur à l'actuel d'au moins 115 m. La végétation était, elle aussi, très différente. La température globale de l'air en surface était de 5 °C plus basse que maintenant en moyenne annuelle. Le cycle hydrologique fut affecté par tous ces changements, avec un impact important sur l'environnement à l'échelle régionale.

L'étude des climats du passé est indispensable pour comprendre les mécanismes de la variabilité naturelle du climat. Cette variabilité reflète l'ensemble de toutes les interactions complexes entre le forçage externe du climat, la dynamique des différentes composantes du système climatique (atmosphère, océan, biosphère, hydrosphère, cryosphère) et les rétroactions au sein même de ce système. Les modèles mis au point pour reproduire les variations climatiques observées sont de plus en plus performants, ce qui autorise leur utilisation pour tenter de prédire le climat futur.

La théorie astronomique des paléoclimats, rendue populaire par Milankovitch [36], permet d'expliquer les variations climatiques aux échelles de 10 000 à 100 000 ans, en particulier la récurrence des cycles glaciaires–interglaciaires du Quaternaire. C'est la variation à long terme des distributions saisonnières et en latitudes de l'énergie solaire, qui joue un rôle clé dans cette explication. Elle est elle-même fonction de trois paramètres astronomiques : l'excentricité (dont les périodes principales se situent autour de 400 et 100 ka), la précession climatique (avec quatre périodes principales autour de 23 et 19 ka) et l'obliquité (dont la période principale est proche de 41 ka) [6].

Un modèle climatique de complexité intermédiaire (LLN 2-D NH) a été développé à Louvain-la-Neuve [22], pour montrer que ces variations de l'insolation induisent bien la récurrence des cycles glaciaires-interglaciaires. Il s'agit d'un modèle latitude–altitude pour l'hémisphère nord. Il représente les principales composantes du système climatique, c'est-à-dire l'atmosphère, la couche de mélange de l'océan, les banquises, les continents, les champs de neige et les inlandsis, ainsi que leurs interactions. Comme le modèle ne possède pas encore de cycle du carbone, la concentration en CO2 dans l'atmosphère y est considérée comme un forçage extérieur, en plus de l'insolation [7,8,12,22]. Validé sur le climat présent, ce modèle s'est montré capable de reproduire les principales caractéristiques de la circulation atmosphérique générale et du cycle saisonnier de la couche de mélange de l'océan, de la banquise et du champ de neige [22]. De plus, de nombreuses expériences de sensibilité ont mis en évidence le rôle des différents forçages, ainsi que l'importance des processus de rétroaction au sein du système climatique.

2 Sensibilité du modèle climatique à la concentration en CO2

Au cours des 400 000 dernières années, la concentration en CO2 atmosphérique a varié entre 180 et 300 ppmv, avec une valeur moyenne de 225 ppmv [39]. Si on impose une faible concentration en CO2, 210 ppmv, correspondant aux valeurs typiques des périodes glaciaires, le modèle est capable de simuler une succession de périodes glaciaires et interglaciaires plus ou moins comparables aux reconstructions géologiques. Dans ces conditions, la concentration en CO2 est, en effet, suffisamment faible pour permettre la mise en place de rétroactions qui amplifient la réponse climatique directe au forçage astronomique. Il s'agit principalement des rétroactions liées à la vapeur d'eau, à l'albédo et aux inlandsis. À l'opposé, une concentration élevée en CO2 ne permet pas la formation de grandes calottes de glace. La réponse climatique au CO2 est cependant loin d'être linéaire. Ainsi, pour une concentration en CO2 intermédiaire (250 ppmv), la réponse du modèle peut être plus proche de celle due à une faible concentration ou de celle due à une concentration élevée et ce, en fonction de l'insolation et de l'état, glacé ou non, dans lequel se trouve la Terre [13].

Alors que les variations de l'insolation sont capables, à elles seules, d'entraîner des variations climatiques de type glaciaire–interglaciaire, les variations de la concentration en CO2 seules ne le sont pas. Pour une distribution d'insolation fixée, le système climatique se maintient, soit dans un état glaciaire, soit dans un état interglaciaire et ce, malgré les variations de la concentration en CO2 [31].

3 Le climat des trois derniers millions d'années

Les premiers inlandsis importants de l'hémisphère nord se sont formés il y a environ 3 Ma (Ma : million d'années). Différentes hypothèses, comme l'effet combiné de l'augmentation de l'amplitude de l'obliquité et de celle de la précession [35], ont été proposées pour expliquer l'intensification des glaciations de l'hémisphère nord. Li et al. [28], utilisant le modèle LLN 2-D NH, ont suggéré une influence complémentaire de la diminution de la concentration en CO2. Le modèle ainsi forcé a simulé trois intervalles de glaciation entre 3 et 2 Ma BP. Entre 2 et 1 Ma BP, le spectre du volume de glace est dominé par une période de 41 ka et ce, en conformité avec les enregistrements de δ18O. Enfin, après 1 Ma BP, simulation et enregistrements géologiques sont tous deux caractérisés par des spectres où, cette fois, la période de 100 ka domine.

Au cours du Pliocène supérieur, la concentration en CO2 était élevée et les calottes de l'hémisphère nord étaient petites. Dans ces conditions, toute augmentation de l'insolation, même modeste, peut entraîner une déglaciation, un âge glaciaire « accidentel » ne pouvant se développer que si les étés sont particulièrement froids. Cette condition n'est satisfaite que pour une combinaison adéquate de la précession et de l'obliquité et uniquement si l'excentricité est grande. Au contraire, durant le Pléistocène supérieur, la concentration en CO2 était plus faible, ce qui a permis le développement de calottes importantes. Dans ce monde « glacé », les interglaciaires sont à présent accidentels et ne peuvent se développer que lorsque l'insolation est très élevée. Cette condition est atteinte pour des valeurs maximales de l'excentricité, pour une obliquité importante et lorsque le solstice d'été est au périhélie.

4 Les 400 000 dernières années

Le modèle climatique, dans sa version actuelle, ne contient pas de cycle du carbone. Nous avons donc dû considérer le CO2 comme un forçage externe et non comme la rétroaction qu'il est en réalité. Des résultats récents [16,19,27] montrent qu'une telle hypothèse de travail ne biaise pas nécessairement les résultats du modèle. Au contraire, l'accélération des changements climatiques, due aux rétroactions du cycle de carbone à laquelle ils arrivent, permet de supposer que notre hypothèse de travail sous-estime en fait l'impact astronomique. De manière générale, la réponse du modèle au forçage astronomique [3] et au CO2 [39] reproduit assez bien les enregistrements géologiques du climat des 400 000 dernières années [5,13,38] (Fig. 1) et, en particulier, le volume de glace estimé par SPECMAP [25,34]. En revanche, le modèle a tendance à surestimer la fonte des calottes durant les interglaciaires, à moins que le volume de glace continentale dans l'hémisphère sud, non simulé dans le modèle, n'ait changé de manière importante à ces époques. Cette fonte des calottes dans le modèle LLN n'empêche toutefois pas la glaciation suivante de se développer.

Fig. 1

Variation au cours des 500 000 dernières années (de haut en bas) : – de l'insolation à 65°N en juin [3] ; – de la concentration en CO2 atmosphérique (ppmv) d'après [39] pour les valeurs postérieures à 414 ka BP et d'après [29] pour celles antérieures à 414 ka BP ; – du volume de glace continentale simulé par le modèle LLN 2-D NH (106 km3) ; – de la courbe composite de 18O dans la chronologie de SPECMAP [25,34] ; – de la courbe composite de δ18O des carottes MD900963 et ODP site 677 [2] ; – du niveau des mers (en ‰) [41]. Les variations de δ18O de l'eau de mer sont dues à des changements de température, de volume de glace et du niveau des mers. La courbe présentée ici montre uniquement la partie du signal liée au changement du niveau des mers. La chronologie de SPECMAP [25,34] est utilisée pour définir les bornes des stades isotopiques.

Variations over the last 500 ka of (from top to bottom): – the mid-month insolation (W m−2) at 65°N in June [3]; – the atmospheric CO2 concentration (ppmv) (from [39] after 414 ka BP and from [29] prior to 414 ka BP); – the simulated northern hemisphere's continental ice volume (106 km3); – the stacked, smoothed oxygen isotope record as a function of age in the SPECMAP timescale [25,34]; – the stacked δ18O record of core MD900963 and ODP site 677 [2]; – the sea level (in ‰) [41]. Changes in ocean water δ18O are attributable to temperature, ice volume, and sea level. This curve shows only the part of the signal related to sea level change. The SPECMAP time scale [25,34] is used for the stage boundaries.

De plus, la comparaison du climat des différents interglaciaires permet de mettre en évidence le rôle qu'y joue l'insolation. En effet, durant les deux stades isotopiques 5 et 11, la concentration en CO2 reste élevée assez longtemps. En revanche, le volume de glace reste faible durant le stade 11 (voir section suivante), alors qu'il varie de manière significative pendant le stade 5. Ce comportement différent entre les stades 5 et 11 ne peut dès lors s'expliquer que par une différence fondamentale entre les distributions d'insolation qui les caractérisent. Au stade 11, l'insolation reste presque constante dans le temps par suite d'une valeur de l'excentricité proche de zéro. Au stade 5, les valeurs relativement élevées de l'excentricité induisent des variations importantes de l'insolation, dont l'amplitude atteint 20 % de la valeur moyenne à 65°N en juin. Finalement, à l'Holocène, le modèle montre une légère dégradation du climat depuis 6000 ans, qualitativement semblable au refroidissement observé au Groenland [26].

Ces expériences confirment que les variations orbitales de la Terre et, par conséquent, de l'insolation jouent un rôle de métronome sur les âges glaciaires [24]. De plus, elles montrent que le CO2 renforce le cycle de 100 ka et améliore les valeurs simulées de la température de surface [23,33]. Elles soulignent aussi toute l'importance de plusieurs processus gouvernant la réponse du système climatique aux divers forçages. Il s'agit des rétroactions albédo–vapeur d'eau–température [10], végétation–albédo [4], de l'effet de continentalité et d'altitude sur les calottes, de la réponse isostatique de leur socle rocheux [17] et de la déstabilisation de leur front sud lors de la fonte [9,11].

5 Un analogue du climat futur : le stade isotopique 11

L'Éémien a souvent été considéré comme un analogue du climat futur. Toutefois, une analyse détaillée montre que les forçages astronomiques à l'œuvre durant l'Holocène et l'Éémien sont fortement différents. En revanche, le climat autour de 400 ka BP montre une meilleure analogie avec le présent ; les paramètres astronomiques y ont un comportement fort similaire à celui du passé récent et du futur proche (prochaines dizaines de milliers d'années). L'excentricité atteindra son prochain minimum dans 27 ka et l'orbite de la Terre sera alors presque circulaire, une situation similaire à celle qui s'est produite il y a 366 ka. Tant durant l'Holocène et le proche futur que durant le stade isotopique 11 (autour de 400 ka BP), l'amplitude des variations de la précession climatique (modulée par l'excentricité) et de l'insolation (dominée par la précession, pour la plupart des latitudes et périodes de l'année) est faible [6]. En tenant compte de ces conditions, le modèle LLN 2-D NH a permis de confirmer que le stade isotopique 11 constituait un bon analogue du climat présent et futur [32]. Cette conclusion vient d'être confirmée en matière de température pour la carotte de Vostok. La déglaciation conduisant au stade 11 y est extrêmement similaire à celle qui a conduit à notre interglaciaire. Ceci laisse supposer qu'il y a des chances que le long interglaciaire du stade 11 (28 000 ans) se répète dans l'avenir [1].

Une analyse de sensibilité a montré que la longueur de l'interglaciaire simulé par le modèle au stade 11 était fortement contrainte par la phase entre l'insolation et la concentration en CO2 atmosphérique [30]. Plus précisément, c'est uniquement si l'insolation et la concentration en CO2 décroissent simultanément que le modèle simule un début de glaciation précoce (autrement dit, un court interglaciaire). Dans les autres cas, la période interglaciaire se prolonge sur plusieurs dizaines de milliers d'années. Or, pour l'interglaciaire actuel, nous savons que l'insolation diminue depuis 10 ka, alors que la concentration en CO2 dans l'atmosphère a continué d'augmenter pendant toute cette période (même hors activités humaines, sauf si l'hypothèse de Ruddiman [40] d'un anthropocène qui aurait déjà commencé il y a 8000 ans, s'avérait être réelle). Ceci explique pourquoi, dans nos simulations, notre interglaciaire persiste pendant plusieurs dizaines de milliers d'années et ce, pour presque tous les scénarios utilisés (Fig. 2). Ces scénarios supposent, (1) soit une concentration en CO2 constante, faible (210 ppmv), (2) soit une concentration variable, identique à celle des 131 000 dernières années, (3) soit un scénario similaire sur lequel se superpose une perturbation liée aux activités humaines augmentant la concentration naturelle jusqu'à 750 ppmv au cours du prochain millénaire. Le volume de glace continentale simulé sous le scénario (1) augmente dès le début de la simulation. Il reste inférieur à 5×106 km3 pendant les 15 000 prochaines années. Ensuite, il augmente plus rapidement et un premier maximum est atteint à 61 ka AP (AP : After Present ; dans le futur) (31×106 km3). Le prochain maximum glaciaire est bien atteint 40 ka plus tard, avec un volume de glace de 33×106 km3. Dans le cas des scénarios (2) et (3), le volume de glace continentale de l'hémisphère nord reste inférieur à 5×106 km3 jusqu'à 50 ka AP environ. La calotte groenlandaise disparaît même complètement pendant quelques milliers d'années dans le cas du scénario (3). Ensuite, les calottes se reforment rapidement pour atteindre un premier maximum vers 65 ka AP, et enfin le prochain maximum glaciaire est atteint à 101 ka AP avec un volume de 34×106 km3. Une faible différence entre le volume de glace simulé dans le cas du scénario naturel (2) et du scénario anthropique (3) subsiste jusque 75 ka AP (Fig. 2).

Fig. 2

Variations au cours des 100 000 prochaines années et pour l'hémisphère nord de la température moyenne annuelle de l'air en surface (à gauche en rouge), du volume de glace continentale (à gauche en bleu) et des précipitations (à droite ; en rouge pour juillet, en bleu pour janvier et en noir en moyenne annuelle) simulés par le modèle LLN 2-D NH, selon trois scénarios de concentration en CO2 atmosphérique : soit (1) une concentration constante (210 ppmv) (tirets), soit (2) une concentration variable, identique à celle des 131 000 dernières années selon la reconstruction issue de la carotte de glace de Vostok (trait plein), soit (3) le scénario (2), avec une perturbation additionnelle liée aux activités humaines, qui augmente la concentration naturelle jusqu'à 750 ppmv au cours du prochain millénaire (pointillés).

Variations over the next 100 ka of the northern hemisphere's annual mean air temperature (left, red curves), continental ice volume (left, blue curves) and precipitation (right; July in red, January in blue and annual mean in black), according to three different CO2 scenarios ((1) a constant 210 ppmv scenario is used (dashed line), (2) the last 131 ka of the Vostok reconstruction are used for the future (full line) and (3) an anthropogenic increase during the next millennium is added on the top of the natural scenario (dotted line)).

Au cours des 100 000 prochaines années, la température moyenne annuelle de l'hémisphère nord montre une tendance générale à la décroissance. Cette tendance est plus faible dans le scénario (1) (0,012 °C ka−1) que dans le scénario (2) (0,027 °C⋅ka−1). Pendant les 50 000 prochaines années, la température simulée sous le scénario (1) est 1 °C plus basse que celle simulée par le scénario (2). La réponse directe à l'augmentation de la concentration en CO2 au cours du prochain millénaire est de 2,8 °C. Bien que cette augmentation de température soit de courte durée, elle reste de 0,6 °C pendant plus de 15 000 ans, avant de reprendre sa valeur naturelle (Fig. 2). Il faut souligner que ce scénario (3) est relativement conservatoire. En effet, il suppose un retour à des concentrations naturelles en CO2 après 1000 ans, autrement dit, on suppose que le système climatique est capable d'assimiler très rapidement des émissions d'origine anthropique. Ce n'est très probablement pas le cas. Une simulation complémentaire, pour laquelle la concentration en CO2 atmosphérique diminue lentement (jusqu'à 270 ppmv dans 100 ka), montre que l'absence de glace continentale dans l'hémisphère nord persiste durant plus de 100 ka. De plus, l'augmentation de température, en moyenne annuelle, est encore supérieure à 2 °C dans 100 ka.

Le modèle climatique LLN 2-D NH contient une représentation simple du cycle hydrologique : en moyenne hémisphérique, la précipitation totale égale l'évaporation. La distribution en latitude est proportionnelle à la distribution actuelle et une correction pour les précipitations sur les calottes est appliquée [22]. Dès lors, les variations des précipitations suivent celles de la température, plus faibles dans le cas du scénario (1) et plus importantes lors d'une augmentation de la concentration en CO2 (scénario (3) ; Fig. 2).

6 Changements climatiques à l'échelle régionale

Le modèle atmosphérique régional (MAR) [20,21], couvrant l'Europe occidentale, a ensuite été utilisé pour étudier les variations climatiques liées à une augmentation de la concentration en CO2 et à un changement des paramètres orbitaux. Il s'agit d'un modèle aux équations primitives, dans lequel la coordonnée verticale est la pression normalisée σ. Les nuages et les précipitations sont représentés par un schéma de type Kessler. Les variables pronostiques et les équations de conservation sont incluses pour les différents hydrométéores (nuages, pluie, neige, grêle) [20]. Le modèle atmosphérique est couplé à un modèle de sol et comprend des schémas radiatifs solaire et infrarouge [18,37] semblables à ceux des modèles de circulation générale (GCM). Le schéma de fermeture de la turbulence utilisé est celui d'ordre 1,5 de Therry et Lacarrère [42]. Une paramétrisation de la convection profonde a également été introduite [15]. Les informations météorologiques aux frontières du domaine sont fournies par un modèle de circulation atmosphérique générale, forcé par des conditions aux frontières observées ou simulées par le modèle de complexité intermédiaire (LLN 2-D NH). Les simulations numériques couvrent une année. Nous présenterons ici uniquement les résultats pour l'hiver (décembre–janvier–février : DJF) et pour l'été (juin–juillet–août : JJA).

La première simulation porte sur les conséquences climatiques d'un triplement de la concentration en CO2 par rapport à la période actuelle (plus précisément une augmentation de 345 ppmv à 1100 ppmv) [14] (Fig. 3). Comme on peut s'y attendre, la température augmente fortement, davantage sur les continents que sur les océans. En été, la région méditerranéenne subit les plus fortes augmentations de température (+7°C en Espagne, +4,5°C en France et en Angleterre). En revanche, en hiver, les augmentations de température sont les plus importantes dans le Nord de l'Europe (supérieure à 7 °C dans le Nord de la France). Le réchauffement des hautes latitudes en hiver est lié au changement d'étendue de la glace de mer sur l'océan et de la neige sur le continent (formation tardive, fonte hâtive, diminution de la couverture estivale). En moyenne globale annuelle, les précipitations augmentent légèrement. Toutefois, la plupart des régions d'Europe occidentale subissent une assez forte réduction des précipitations estivales. Par exemple, la réduction est supérieure à 1 mm j−1 au Nord-Est de la France, la moyenne actuelle y étant de 1,9 à 2,7 mm j−1 (région Meuse–Haute-Marne). En hiver, pour les concentrations actuelles en CO2, les précipitations en Europe sont essentiellement de type stratiforme. En revanche, l'augmentation de la concentration en CO2 et le réchauffement qui en découle modifient cette caractéristique, une partie importante des précipitations hivernales prenant un caractère convectif. La comparaison entre ces deux états devient dès lors plus délicate, de même que l'établissement de conclusions générales à partir d'une seule année de simulation.

Fig. 3

Résultat du modèle MAR. Différence entre la simulation forcée par une concentration en CO2 atmosphérique de 1100 ppmv et la simulation forcée par la valeur actuelle (345 ppmv). Dans les deux cas, on a utilisé la configuration orbitale actuelle. Différence de température de surface (°C ; à gauche) et de précipitation totale (mm j−1 ; à droite) en DJF.

MAR results. Difference between the experiment using an 1100 ppmv atmospheric CO2 concentration and the experiment with a 345 ppmv atmospheric CO2 concentration. The present-day orbital parameters are used in both experiments. Difference in surface temperature (K; left) and in total precipitation (mm d−1; right) in DJF.

La seconde simulation porte sur l'effet d'un changement du forçage astronomique [14]. L'insolation considérée est celle de 67 ka AP. Elle est caractérisée par une augmentation dans l'hémisphère nord, de mars à août, pouvant atteindre 40 W m−2 en juin à 65°N, et par une diminution de l'insolation dans l'hémisphère nord le reste de l'année, inférieure à 4 W m−2 en septembre à 65°N. Le changement de température qui en découle reflète le changement de forçage. Les températures d'hiver diminuent, de manière plus marquée aux hautes latitudes à cause des rétroactions liées à l'augmentation de la couverture neigeuse, alors que les températures d'été augmentent, de 6 °C en Espagne, de 3,5 °C en Angleterre et de 2,5 °C dans le Nord de la France. En hiver, les précipitations augmentent significativement sur les îles Britanniques (+1,15 mmj−1), sur l'Ouest de la France, sur la péninsule Ibérique et en région méditerranéenne. Cette augmentation des précipitations se propage depuis le Groenland et est issue des conditions aux frontières tirées du modèle global (GCM). En été, il y a une nette fracture entre le nord, où les précipitations augmentent (+0,09 mmj−1 en Angleterre, +0,88 mmj−1 dans le Nord de l'Allemagne), et le sud, où elles diminuent (0,6 mmj−1 en Espagne).

Enfin, nous avons étudié le cas d'une augmentation de la concentration en CO2 atmosphérique (de 345 ppmv à 550 ppmv) sous les conditions de forçage astronomique de 67 ka AP [14]. Tout comme dans le cas d'une augmentation de la concentration en CO2 sous le forçage astronomique actuel, la température augmente, tant en été qu'en hiver. En hiver, elle augmente essentiellement sur le continent, du Sud de la France au Nord de l'Allemagne. Bien que l'augmentation du CO2 atmosphérique soit plus faible que dans l'expérience « période actuelle », l'amplitude des augmentations de température est semblable dans les deux expériences (de l'ordre de +10°C en France, en République tchèque et en Allemagne). Toutefois, il pourrait s'agir d'un artefact lié au choix des conditions aux frontières. En revanche, l'augmentation de température en été est réduite par rapport à l'expérience « période actuelle » et est principalement localisée dans le Sud de l'Europe.

Une caractéristique commune aux deux expériences d'augmentation de la concentration en CO2 concerne l'effet sur les changements de climat en montagne. En hiver, le réchauffement est réduit en altitude par rapport à la plaine, probablement parce que l'altitude permet à la neige de se maintenir et ainsi à la rétroaction température–albédo de continuer à jouer son rôle. En revanche, en été, l'augmentation de température est plus importante en montagne qu'en plaine. En montagne, le réchauffement initial entraîne la fonte des neiges et permet la mise en place de la rétroaction positive albédo–température, qui amplifie le réchauffement initial.

En hiver, le modèle suggère une forte diminution des précipitations dans le Bassin méditerranéen (1,82 mmj−1 en Espagne), alors que les précipitations augmentent dans le nord (+1,22 mmj−1 en France, +2,01 mmj−1 en Angleterre, +1,38 mmj−1 dans le Nord de l'Allemagne). En été, comme pour l'époque actuelle, une augmentation de la concentration en CO2 atmosphérique entraîne une diminution des précipitations sur pratiquement toute l'Europe (0,16 mmj−1 en Angleterre, 0,11 mmj−1 en Espagne), sauf sur quelques régions très localisées (+0,41 mmj−1 en République tchèque, +0,35 mmj−1 dans le Nord de l'Allemagne).

7 Conclusion

Les variations de l'énergie solaire reçue par la Terre, induites par les variations de l'orbite de la Terre autour du Soleil et de l'inclinaison de son axe de rotation, jouent un rôle fondamental dans les variations climatiques du Quaternaire, aux échelles de temps allant de la dizaine à plusieurs centaines de milliers d'années. Un modèle climatique, tel que LLN 2-D NH [22], permet de mieux comprendre le fonctionnement du système climatique et les mécanismes de rétroaction qui entrent en jeu lors des changements climatiques. Outre le rôle prépondérant de l'insolation, le CO2 joue un rôle d'amplificateur du signal initié par ces variations de l'insolation. Toutefois, le rôle du CO2 augmente significativement lorsque les variations de l'insolation sont faibles, comme c'est le cas au stade 11 et durant l'Holocène. Le modèle LLN 2-D NH, appliqué au futur, permet également d'établir des projections pour les grandes lignes du climat futur, suggérant que le climat interglaciaire puisse se maintenir pendant 50 000 ans. Il ne faut toutefois pas oublier qu'un nombre important d'hypothèses sous-tend ce modèle et qu'il ne représente que l'hémisphère nord.

L'utilisation d'une hiérarchie de modèles (modèle global de complexité intermédiaire, modèle de circulation générale, modèle régional) permet d'estimer la réponse du système climatique, à l'échelle régionale, pour des états climatiques futurs. Les résultats obtenus ici devront encore être confirmés, notamment par l'application de la même procédure d'imbrication par d'autres groupes, par des analyses de sensibilité plus poussées, par des modèles plus complets de la circulation océanique et du cycle du carbone et par l'étude de l'impact de la variabilité naturelle sur les résultats.

Acknowledgments

Merci à ceux qui, par les nombreuses discussions qu'ils ont eues avec nous, ont permis de faire avancer ce travail, en particulier Dominique Raynaud, Michel Crucifix, Philippe Marbaix, Xavier Fettweis. Ce travail a été partiellement réalisé dans le cadre du projet EU Bioclim (contrat FIKW-CT-2000–00024). Nous remercions également le coordinateur (Andra), le secrétariat technique (Quantisci) et les participants au projet Bioclim.


Bibliographie

[1] L. Augustin; C. Barbante; P.R.F. Barnes; J.-M. Barnola; M. Bigler; E. Castellano; O. Cattani; J. Chappellaz; D. Dahl-Jensen; B. Delmonte; G. Dreyfus; G. Durand; S. Falourd; H. Fischer; J. Flückiger; M.E. Hansson; P. Huybrechts; G. Jugie; S.J. Johnsen; J. Jouzel; P. Kaufmann; J. Kipfstuhl; F. Lambert; V.Y. Lipenkov; G.C. Littot; A. Longinelli; R. Lorrain; V. Maggi; V. Masson-Delmotte; H. Miller; R. Mulvaney; J. Oerlemans; H. Oerter; G. Orombelli; F. Parrenin; D.A. Peel; J.-R. Petit; D. Raynaud; C. Ritz; U. Ruth; J. Schwander; U. Siegenthaler; R. Souchez; B. Stauffer; J.P. Steffensen; B. Stenni; T.F. Stocker; I.E. Tabacco; R. Udisti; R.S.W. Van de Wal; M. van den Broeke; J. Weiss; F. Wilhelms; J.-G. Winther; E.W. Wolff; M. Zucchelli Eight glacial cycles from an Antarctic ice core, Nature, Volume 429 (2004), pp. 623-628 | DOI

[2] F.C. Bassinot; L.D. Labeyrie; E. Vincent; X. Quidelleur; N.J. Shackleton; Y. Lancelot The astronomical theory of climate and the age of the Brunhes–Matuyama magnetic reversal, Earth Planet. Sci. Lett., Volume 126 (1994), pp. 91-108

[3] A. Berger Long-term variations of daily insolation and Quaternary climatic changes, J. Atmos. Sci., Volume 35 (1978) no. 12, pp. 2362-2367

[4] A. Berger The role of CO2, sea level, and vegetation during the Milankovitch-forced glacial–interglacial cycles (L.O. Bengtssen; C.U. Hammer, eds.), Geopshere–Biosphere interactions and climate, Cambridge University Press, Cambridge, 2001

[5] A. Berger; M.-F. Loutre Modelling the climate response to astronomical and CO2 forcings, C. R. Acad. Sci., Ser. IIa, Volume 323 (1996) no. 1, pp. 1-16

[6] A. Berger; M.-F. Loutre Théorie astronomique des paléoclimats, C. R. Geoscience, Volume 336 (2004), pp. 701-709

[7] A. Berger; T. Fichefet; H. Gallée; C. Tricot; I. Marsiat; J.-P. van Ypersele Astronomical forcing of the last glacial–interglacial cycle (P.J. Crutzen; J.-C. Gérard; R. Zander, eds.), Our Changing Atmosphere, université de Liège, institut d'astrophysique, Cointe-Ougrée, 1989, pp. 353-382

[8] A. Berger; T. Fichefet; H. Gallée; I. Marsiat; C. Tricot; J.-P. van Ypersele Physical interactions within a coupled climate model over the last glacial–interglacial cycle, Trans. R. Soc. Edinb.: Earth Sci., Volume 82 (1990), pp. 357-369

[9] A. Berger; T. Fichefet; H. Gallée; C. Tricot; J.-P. van Ypersele Entering the glaciation with a 2-D coupled climate model, Q. Sci. Rev., Volume 11 (1992) no. 4, pp. 481-493

[10] A. Berger; C. Tricot; H. Gallée; M.-F. Loutre Water vapour, CO2 and insolation over the last glacial–interglacial cycles, Phil. Trans. R. Soc. Lond. B, Volume 341 (1993), pp. 253-261

[11] A. Berger; H. Gallée; C. Tricot Glaciation and deglaciation mechanisms in a coupled 2-D climate-ice sheet model, J. Glaciol., Volume 39 (1993) no. 131, pp. 45-49

[12] A. Berger; C. Tricot; H. Gallée; T. Fichefet; M.-F. Loutre The last two glacial–interglacial cycles simulated by the LLN model, NATO ASI Ser. Ser. I Global Environ. Change, Volume 22 (1994), pp. 411-452

[13] A. Berger; M.-F. Loutre; H. Gallée Sensitivity of the LLN climate model to the astronomical and CO2 forcings over the last 200 ka, Clim. Dyn., Volume 14 (1998) no. 9, pp. 615-629

[14] , 2003 http://www.andra.fr/bioclim (BIOCLIM Report D6a Regional climatic characteristics for the European sites at specific times: the dynamical downscaling, Workpackage 2, available on)

[15] O. Brasseur, Mesoscale modeling of extreme meteorological events over Belgium, dissertation doctorale, fac. sci. appl., université catholique de Louvain-la-Neuve, Belgique, 1999

[16] P.M. Cox; R.A. Betts; C.D. Jones; S.A. Spall; I.J. Totterdell Acceleration of global warming due to carbon-cycle feedbacks in a coupled climate model, Nature, Volume 408 (2000) no. 6809, pp. 184-187

[17] M. Crucifix; M.-F. Loutre; K. Lambeck; A. Berger Effect of isostatic rebound on modelled ice volume variations during the last 200 ka, Earth Planet. Sci. Lett., Volume 184 (2001), pp. 623-633

[18] Y. Fouquart; B. Bonnel Computation of the solar heating of the Earth's atmosphere: A new parameterization, Beitr. Phys. Atmos., Volume 53 (1980), pp. 35-62

[19] P. Friedlingstein; L. Bopp; P. Ciais; J.-L. Dufresne; L. Fairhead; H. Le Treut; P. Monfray; J. Orr Positive feedback between future climate change and the carbon cycle, Geophys. Res. Lett., Volume 28 (2001) no. 8, pp. 1543-1546

[20] H. Gallée Simulation of mesocyclonic activity in the Ross Sea, Antarctica, Mon. Weather Rev., Volume 123 (1995) no. 7, pp. 2051-2069

[21] H. Gallée; G. Schayes Dynamical aspects of katabatic winds evolution in the Antarctic coastal zone, Bound-Lay. Meteorol., Volume 59 (1994), pp. 141-161

[22] H. Gallée; J.-P. van Ypersele; T. Fichefet; C. Tricot; A. Berger Simulation of the last glacial cycle by a coupled sectorially averaged climate – ice-sheet model. I. The climate model, J. Geophys. Res., Volume 96 (1991), pp. 13139-13161

[23] H. Gallée; J.-P. van Ypersele; T. Fichefet; I. Marsiat; C. Tricot; A. Berger Simulation of the last glacial cycle by a coupled, sectorially averaged climate – ice-sheet model. II. Response to insolation and CO2 variation, J. Geophys. Res., Volume 97 (1992) no. D14, pp. 15713-15740

[24] J.D. Hays; J. Imbrie; N.J. Shackleton Variations in the Earth's orbit: pacemaker of the ice ages, Science, Volume 194 (1976), pp. 1121-1132

[25] J. Imbrie; J.D. Hays; D.G. Martinson; A. Mcintyre; A.C. Mix; J.J. Morley; N.G. Pisias; W.L. Prell; N.J. Shackleton The orbital theory of Pleistocene climate: support from a revised chronology of the marine 18O record (A. Berger et al., eds.), Milankovitch and Climate, D. Reidel, Hingham Mass, 1984, pp. 269-305

[26] S.J. Johnsen; D. Dahl-Jensen; W. Dansgaard; N. Gundestrup Greenland palaeotemperatures derived from GRIP bore hole temperature and ice core isotope profiles, Tellus B, Volume 47 (1995) no. 5, pp. 624-629

[27] C.D. Jones; P.M. Cox; R.L.H. Essery; D.L. Roberts; M.J. Woodage Strong carbon cycle feedbacks in a climate model with interactive CO2 and sulphate aerosols, Geophys. Res. Lett., Volume 30 (2003) no. 9, p. 1479 (32.1–32.4) | DOI

[28] X.S. Li; A. Berger; M.-F. Loutre; M.A. Maslin; G.H. Hang; R. Tiedemann Simulating Late Pliocene Northern Hemisphere climate with the LLN 2-D model, Geophys. Res. Lett., Volume 25 (1998) no. 6, pp. 915-918

[29] X.S. Li; A. Berger; M.-F. Loutre CO2 and northern hemisphere ice volume variations over the Middle and Late Quaternary, Clim. Dyn., Volume 14 (1998) no. 7–8, pp. 537-544

[30] M.-F. Loutre Clues from MIS 11 to predict the future climate – a modeling point of view, Earth Planet. Sci. Lett., Volume 212 (2003), pp. 213-224

[31] M.-F. Loutre; A. Berger No glacial–interglacial cycle in the ice volume simulated under a constant astronomical forcing and a variable CO2, Geophys. Res. Lett., Volume 27 (2000) no. 6, pp. 783-786

[32] M.-F. Loutre; A. Berger Stage 11 as an analogue for the present interglacial, Global Planet. Change, Volume 36 (2003), pp. 209-217

[33] M.-F. Loutre; A. Berger; A. Dutrieux; H. Gallée The response of the LLN climate model to the astronomical forcing over the last glacial–interglacial cycle, Terra Nostra, Schriften der Alfred-Wegener-Stiftung, Volume 1 (1994), pp. 11-15

[34] D.G. Martinson; N.G. Pisias; J.D. Hays; J. Imbrie; T.C. Moore; N.J. Shackleton Age dating and the orbital theory of the ice ages: development of a high-resolution 0 to 300 000-year chronostratigraphy, Quat. Res., Volume 27 (1987) no. 1, pp. 1-29

[35] M.A. Maslin; G.H. Haug; M. Sarnthein; R. Tiedemann; H. Erlenkeuser; R. Stax Northwest Pacific Site 882: The initiation of northern hemisphere glaciation, Proc. Ocean Drilling Program, Scientific Results, Volume 145 (1995), pp. 315-329

[36] M. Milankovitch Kanon der Erdbestrahlung und seine Anwendung auf des Eiszeitenproblem, Royal Serbian Sciences, Spec. publ. 132, vol. 33, Section of Mathematical and Natural Sciences, Belgrade, 1941, p. 633 (‘Canon of Insolation and the Ice Age Problem’, English Translation by the Israël Program for the U.S. Department of Commerce and the National Science Foundation, Washington DC, 1969, and by Zavod za Udzbenike I nastavna Sredstva, in cooperation with Muzej nauke I technike Srpske akademije nauka I umetnosti, Beograd, 1998)

[37] J. Morcrette, Sur la paramétrisation du rayonnement dans les modèles de la circulation générale atmosphérique, thèse d'État, USTL, Lille, 1984

[38] L. Pépin; D. Raynaud; J.-M. Barnola; M.-F. Loutre Hemispheric roles of climate forcings during glacial–interglacial transitions as deduced from the Vostok record and LLN-2D model experiments, J. Geophys. Res., Volume 106 (2001) no. D23, pp. 31885-31892

[39] J.R. Petit; J. Jouzel; D. Raynaud; N.I. Barkov; J.-M. Barnola; I. Basile; M. Bender; J. Chappellaz; M. Davis; G. Delaygue; M. Delmotte; V.M. Kotlyakov; M. Legrand; V.Y. Lipenkov; C. Lorius; L. Pépin; C. Ritz; E. Saltman; M. Stievenard Climate and atmospheric history of the past 420 000 years from Vostok ice core, Antarctica, Nature, Volume 399 (1999) no. 6735, pp. 429-436

[40] W.F. Ruddiman The anthropogenic greenhouse era began thousands of years ago, Clim. Change, Volume 61 (2003), pp. 261-293

[41] N.J. Shackleton The 100 000 year ice-age cycle identified and found to lag temperature, carbon dioxide and orbital eccentricity, Science, Volume 289 (2000), pp. 1897-1902

[42] G. Therry; P. Lacarrère Improving the eddy kinetic energy model for the planetary boundary-layer description, Bound-Lay. Meteorol., Volume 25 (1982), pp. 63-88


Commentaires - Politique