Outline
Comptes Rendus

External Geophysics, Climate and Environment (Climate)
Rapid climatic changes of the last 90 kyr recorded on the European continent
Comptes Rendus. Géoscience, Volume 337 (2005) no. 10-11, pp. 970-982.

Abstracts

Millennial-scale events of the last glacial period (Dansgaard–Oeschger events (D/O)), characterized by great temperature changes in the ice core records, have also influenced the vegetation and the climate of the European continent. Three types of continental archives – the loess deposits, the pollen assemblage record of a continental margin core and speleothems – spread on western Europe along a SW–NE profile, give information about the humidity, temperature and vegetation changes that were due to these millennial scale climatic changes.

Les événements millénaires qui ponctuent la dernière période glaciaire – événements de Dansgaard/Oeschger (D/O) –, marqués par de fortes variations de température dans les enregistrements des glaces du Groenland, ont aussi fortement influencé la végétation et le climat sur le continent européen. Trois types d'archives continentales – les lœss, les pollens d'une carotte de la marge continentale et les spéléothèmes – répartis en Europe de l'Ouest sur un profil SW–NE, apportent des informations complémentaires sur les changements de végétation, d'humidité, de température, ainsi que sur la chronologie de cette variabilité climatique millénaire.

Metadata
Accepted:
Published online:
DOI: 10.1016/j.crte.2005.04.010
Keywords: Dansgaard/Oeschger, Mediterranean vegetation, Loess, Speleothems, Europe
Mots-clés : Dansgaard/Oeschger, Végétation méditerranéenne, Lœss, Spéléothèmes, Europe

Dominique Genty 1; Nathalie Combourieu Nebout 1; Christine Hatté 1; Dominique Blamart 1; Bassam Ghaleb 2; Lauren Isabello 2

1 Laboratoire des sciences du climat et de l'environnement, UMR CEA/CNRS 1572, bât. 709, L'Orme des Merisiers, CEA Saclay, 91191 Gif-sur-Yvette cedex, France
2 GEOTOP, UQUAM, C.P. 8888, succursale Centre-Ville, H3C 3P8, Montréal, Canada
@article{CRGEOS_2005__337_10-11_970_0,
     author = {Dominique Genty and Nathalie Combourieu Nebout and Christine Hatt\'e and Dominique Blamart and Bassam Ghaleb and Lauren Isabello},
     title = {Rapid climatic changes of the last 90 kyr recorded on the {European} continent},
     journal = {Comptes Rendus. G\'eoscience},
     pages = {970--982},
     publisher = {Elsevier},
     volume = {337},
     number = {10-11},
     year = {2005},
     doi = {10.1016/j.crte.2005.04.010},
     language = {en},
}
TY  - JOUR
AU  - Dominique Genty
AU  - Nathalie Combourieu Nebout
AU  - Christine Hatté
AU  - Dominique Blamart
AU  - Bassam Ghaleb
AU  - Lauren Isabello
TI  - Rapid climatic changes of the last 90 kyr recorded on the European continent
JO  - Comptes Rendus. Géoscience
PY  - 2005
SP  - 970
EP  - 982
VL  - 337
IS  - 10-11
PB  - Elsevier
DO  - 10.1016/j.crte.2005.04.010
LA  - en
ID  - CRGEOS_2005__337_10-11_970_0
ER  - 
%0 Journal Article
%A Dominique Genty
%A Nathalie Combourieu Nebout
%A Christine Hatté
%A Dominique Blamart
%A Bassam Ghaleb
%A Lauren Isabello
%T Rapid climatic changes of the last 90 kyr recorded on the European continent
%J Comptes Rendus. Géoscience
%D 2005
%P 970-982
%V 337
%N 10-11
%I Elsevier
%R 10.1016/j.crte.2005.04.010
%G en
%F CRGEOS_2005__337_10-11_970_0
Dominique Genty; Nathalie Combourieu Nebout; Christine Hatté; Dominique Blamart; Bassam Ghaleb; Lauren Isabello. Rapid climatic changes of the last 90 kyr recorded on the European continent. Comptes Rendus. Géoscience, Volume 337 (2005) no. 10-11, pp. 970-982. doi : 10.1016/j.crte.2005.04.010. https://comptes-rendus.academie-sciences.fr/geoscience/articles/10.1016/j.crte.2005.04.010/

Original version of the full text

Version française abrégée

1 Introduction

Avec le récent carottage de NGRIP, ce sont maintenant vingt-cinq événements climatiques de type Dansgaard/Oeschger (D/O) qui sont dénombrés clairement dans les carottes de glace du Groenland [36]. Plusieurs d'entre eux ont atteint des amplitudes thermiques de plus de 10 °C au Groenland [25,30,31] et ont été clairement enregistrés dans les carottes marines de l'Atlantique et de la Méditerranée, démontrant un couplage important entre l'atmosphère et l'océan, au moins pour l'hémisphère nord [5–7]. La conséquence de ces événements sur le continent est, en revanche, beaucoup moins bien connue, du fait de la rareté d'archives bien datées avec un signal climatique bien préservé. Ce que l'on sait, grâce à l'analyse palynologique de plusieurs lacs du Massif central et d'Italie, est que l'influence de ces variations climatiques millénaires a été ressentie par la végétation jusqu'en Méditerranée centrale [1,16]. Les D/O ont aussi eu un impact important sur la circulation atmosphérique et la mousson d'Asie du Sud-Est en particulier, comme l'a montré le remarquable enregistrement des spéléothèmes de la grotte de Hulu (Chine) [51]. Les enregistrements montrant clairement l'impact des D/O sur l'environnement (climat, végétation) en Europe au cours de la dernière période glaciaire sont rares. Nous présentons ici les résultats de l'analyse de trois types d'archives paléoenvironnementales, qui apportent des informations complémentaires sur l'impact des D/O en Europe de l'Ouest : (1) l'enregistrement palynologique de la carotte ODP 976 en mer l'Alboran, qui nous renseigne directement sur la végétation du Sud de l'Espagne et de l'Afrique du Nord ; (2) un enregistrement isotopique (δ13C) de quatre stalagmites de la grotte de Villars (Sud-Ouest de la France), qui montre la variation de l'activité végétale et biopédologique liée aux D/O et qui donne une chronologie absolue précise pour plusieurs de ces événements, avec toutefois quatre hiatus liés à des seuils climatiques ; (3) la séquence de lœss de Nußloch (Ouest de l'Allemagne), où une modélisation inverse utilisant BIOME4 et le δ13C de la matière organique (MO) a permis une reconstruction des paléoprécipitations, mais qui a malheureusement une résolution temporelle inégale.

2 Sites et chronologies

Les trois sites s'alignent selon une direction SW–NE depuis le Sud de l'Espagne pour la carotte marine ODP 976 (36°12′N, 4°18′W, 1,108 m sous le niveau de la mer [8]) jusqu'à l'Allemagne de l'Ouest, où se situe la séquence de lœss de Nußloch (49°19′N, 8°43′E, alt.=221 m ; [20]) ; la grotte de Villars (Dordogne : 45°30′N, 0°50′E, altitude 175 m ; [11]) se situant entre ces deux sites.

La chronologie du carottage ODP 976 est basée sur 13 analyses 14C (AMS) [8]. Les âges 14C, corrigés de l'effet réservoir et calibrés, concernent la partie récente de l'enregistrement ; la partie plus ancienne est calée, d'une part, avec les stades isotopiques du site ODP 976 et, d'autre part, avec la chronologie de la carotte MD 95–2042, située au large du Portugal [8,14,43–45]. Enfin, les âges ont été ajustés en comparant les courbes palynologiques avec les stades isotopiques de NGRIP [36].

Soixante-neuf datations par U/Th par thermo-ionisation par spectrométrie de masse (TIMS) ont été effectuées sur les quatre stalagmites de la grotte de Villars ([12] et données non publiées ; points datés en bas de la Fig. 2d). Pour la plupart des échantillons, l'erreur sur l'âge à 2σ est inférieure à 2 %, donnant ainsi une des meilleures chronologies absolues sur le continent européen pour cette période.

Fig. 2

Palaeoclimatic records of the last 90 kyr in western Europe. (a) Blue: NGRIP ice core δ18O with the ss09sea chronology [24,36], black: GISP2 δ18O [15]. (b) Nußloch loess sequence: red, δ13C of the OM; black: mean annual rainfall with 2 σ errors (grey curves); at the bottom of the curve: 14C and OSL dated points with 2 σ and 1 σ error bars, respectively. (c) Pollen records of the ODP 976 marine core; green: temperate taxa (see text); red: semi-desertic taxa. (d) Villars stalagmite δ13C record with dated points at the bottom with 2 σ error bars. Hulu cave speleothem δ18O record [51] (black dashed line).

Enregistrements paléoclimatiques des derniers 90 ka en Europe de l'Ouest. (a) En bleu : δ18O de la glace de NGRIP [36], avec la chronologie ss09sea [24] ; en noir : δ18O de GISP2 [15]. (b) Séquence loessique de Nußloch. En rouge : δ13C de la MO ; en noir : précipitations moyennes annuelles (en mm), avec les erreurs à 2 σ (gris clair) ; en bas, points datés par 14C et OSL, avec les erreurs à 2 σ pour 14C et 1 σ pour OSL. (c) Enregistrement palynologique de la carotte ODP 976 de la mer d'Alboran. En vert : assemblage des taxons tempérés (arbres etc., voir texte), en rouge: assemblage des taxons semi-désertiques. (d) δ13C des stalagmites de la grotte de Villars avec, en bas, points datés en UTh (TIMS) avec leur incertitude à 2 σ. δ18O des spéléothèmes de la grotte de Hulu, Chine [51] (courbe noire).

La chronologie de la séquence loessique de Nußloch est basée sur 22 datations 14C sur matière organique lœssique (15), bois (4), os (1), carbonate de coquille de gastéropodes (1), conchioline (1) de 3 à 42 kyr BP et 21 datations IRSL/OSL de 19 à 122 kyr [21,32]. Les ages 14C postérieurs à 20 kyr BP sont calibrés selon Calib 4.3 [47], et les plus anciens sont corrigés selon la proposition de Kitagawa et Van der Plicht [29]. La chronologie utilisée ici est obtenue par extrapolation linéaire entre les valeurs médianes des intervalles de datation.

3 Résultats

Les proxy présentés ici sont les suivants :

  • – pour le site ODP 976 : assemblages polliniques : (1) de taxons tempérés composés d'arbres tels que Quercus, Fagus, Carpinus, Corylus, Alnus, Betula, Tilia, Ulmus... , associés avec Ericaceae, qui reflètent des conditions relativement chaudes et humides, typiques d'un interstade ; (2) de taxons semi-désertiques composés de Artemisia, Amaranthaceae–Chenopodiaceae et Ephedra, indicateurs de conditions froides et sèches, typiques des stadiaires ;
  • – pour le site de Villars : le δ13C de la calcite des stalagmites, qui est lié à la densité de végétation et à l'activité biopédologique, elles-mêmes dépendantes du climat ; ainsi, un δ13C faible (env. 10) est caractéristique de conditions tempérées relativement chaudes et humides, alors qu'un δ13C élevé (env. 6) est lié à des épisodes froids où l'activité végétale est fortement ralentie [12] ;
  • – pour la séquence de Nußloch : le δ13C de la MO, qui reflète la signature isotopique originale de la végétation. L'amplitude de variation des δ13C, de −23 à 25, montre que la végétation est restée de type photosynthétique C3 [20]. Aussi, les variations de δ13C observées sont-elles interprétées en termes de changement environnemental faisant intervenir la composition isotopique et la concentration du CO2 de l'atmosphère, la disponibilité en eau, la température, le type de sol et l'insolation. Une modélisation inverse avec Biome4 [19,22,26] permet de déconvoluer le signal, d'éliminer les composantes globales et, ainsi, de simuler les conditions climatiques ayant généré de tels δ13C et donc de reconstituer les paléoprécipitations [19].

Grâce à la nature différente de ces trois enregistrements paléoclimatiques, on obtient un ensemble cohérent et complémentaire des conditions environnementales qui prévalaient en Europe occidentale au cours des derniers 90 ka (Fig. 2). Alors que les assemblages polliniques du site ODP 976 montrent de façon remarquable les événements froids (Heinrich H1 à H6) grâce aux taxons semi-désertiques, les événements chauds et humides sont clairement enregistrés dans les taxons tempérés de la même carotte, mais aussi dans le δ13C et les phases de croissance des stalagmites de Villars (D/O #7 à D/O #20 ; Fig. 2). Pendant les épisodes de froid extrême, qui correspondent globalement aux stades MIS 4 et MIS 2, la croissance des stalagmites de Villars s'est arrêtée (hiatus D2 et D4 ; car l'activité végétale s'est fortement ralentie ainsi que l'écoulement), alors que la sédimentation s'est poursuivie à Nußloch (MIS 2) et au site ODP 976 permettant un enregistrement clair des variations climatiques pendant ces périodes.

3.1 Assemblages polliniques de ODP 976

Au cours des derniers 90 000 ans, les courbes polliniques obtenues sur le site ODP 976 montrent une rythmicité de grande amplitude, qui traduit des modifications cycliques de la végétation entre des forêts tempérées décidues à Quercus (chêne) et des semi-déserts à Artemisia (armoise) (Fig. 2c). Ces changements de végétation évoquent des variations climatiques entre des conditions chaudes/humides et froides/sèches. Ils sont corrélés avec les oscillations de Dansgaard/Oeschger répertoriées dans les carottes de glace du Groenland. Ainsi, les interstades correspondent au développement de la forêt tempérée, tandis que les stadiaires sont exprimés par l'extension du semi-désert. Six stadiaires sont particulièrement drastiques, avec un semi-désert très développé et une forêt réduite à son minimum. Ils traduisent une intense aridification de l'Ouest méditerranéen et sont associés, dans les environnements marins, à des périodes de refroidissement des eaux de surface qui sont corrélées avec les événements de Heinrich 1 à 6 [8,27].

3.2 Isotopes des spéléothèmes de Villars

Le profil en δ13C de Villars, fait à partir de quatre stalagmites, possède quatre hiatus, qui correspondent à des arrêts de croissance liés à des conditions froides (D2, D4 ; Fig. 2) ou à d'autres causes toujours incertaines (froid, inondations, D1, D3). Les périodes d'arrêt (67,4–61,2 ka et 31,6–15,9 ka) sont cernées par une dérive progressive du δ13C, démontrant le changement graduel des conditions climatiques à leurs abords. La phase froide correspondant au hiatus D2 (ou Villars cold phase) est un élément majeur de cette séquence, qui montre qu'un climat froid a empêché l'infiltration et la croissance des stalagmites pendant plus de six mille ans. On peut l'associer au MIS 4 et au H6, lui aussi bien enregistré, d'une part, dans la carotte ODP 976 par un pic important (en amplitude et en durée) de taxons semi-désertiques (Fig. 2c) et, d'autre part, dans le δ18O de NGRIP, avec des valeurs minimales proches de 44 (Fig. 2a). Dans la séquence de Villars, les D/O montrent une variation typique de 34 dans le δ13C (Fig. 2d), qui peuvent être mises en correspondance avec les δ18O élevés des carottes de glace (Fig. 2a). Le D/O #8 apparaît comme le dernier événement chaud majeur avant le maximum glaciaire ; il s'inscrit dans une tendance générale au refroidissement qui semble avoir débuté dès l'optimum du D/O #12 il y a environ 45 ka et s'est terminé avec les arrêts synchrones des deux stalagmite Vil-stm9 et Vil-stm27 il y a 31,7±0,1 ka (Fig. 2d). Les phases froides qui entourent les D/O sont en accord avec les pics des taxons semi-désertiques de ODP 976 et correspondent aux événements de Heinrich : H4 et H5 sont bien visibles, avec des valeurs de δ13C élevées (attention à l'échelle inversée sur la Fig. 2d), H6 correspondant au hiatus D2, H2 et H3 au hiatus D4. Mais l'événement chaud le plus marqué est certainement le D/O #12, qui a atteint son maximum il y a 45,3±0,4 ka et dont les valeurs de δ13C sont proches de celles des dépôts de calcite moderne (10,5 à 11,2), ce qui suggère des conditions climatiques voisines de l'Actuel. Enfin, la partie la plus récente de la séquence de Villars montre clairement la déglaciation : le réchauffement, marqué par la diminution du δ13C et par la reprise du concrétionnement, est daté de 15,9±0,35 ka ; le δ13C diminue alors de plus de 5, jusqu'à un optimum daté à 9,9±0,17 ka.

3.3 Séquence de Nußloch

La modélisation en mode inverse utilisant le δ13C de la séquence lœssique de Nußloch simule des précipitations variant de 280−60+120 à 460−160+80 mman−1 tout au long de la dernière glaciation, bien en dessous de la valeur actuelle de 800 mman−1. Cette aridité est en accord avec les études pédologiques et malacologiques [2,35].

Avant 28 kyr BP, les paléoprécipitations présentent de grandes oscillations de 150 à 200 mman−1. Ces fluctuations sont similaires et contemporaines des événements de D/O. Aussi, l'événement proche du D/O 14, associé à une augmentation de température de 10 à 12 °C au centre de la calotte du Groenland [9,30,33], semble correspondre à une augmentation des précipitations de 200 mman−1 (+70 %). Cette similarité entre l'enregistrement climatique de Nußloch et celui des glaces du Groenland souligne le lien étroit entre le régime de précipitation en Europe du Nord-Ouest, la taille des calottes de glace polaire et continentales, le niveau marin et par conséquent, la position du front polaire. Toutefois, de par les incertitudes sur la chronologie ancienne de la séquence de loess, la mise en coïncidence des évènements de D/O entre les lœss et les enregistrements glaciaires doit rester prudente.

Alors que le niveau marin est au plus bas (90 m sous le niveau actuel) entre 28 000 et 17 000 ans [49], le régime des précipitations est relativement constant, autour de 280 mm an−1. Dans ces conditions particulières où la Manche est émergée, l'influence des flux d'air océaniques sur Nußloch est considérablement réduite. Nußloch est donc moins sensible aux changements de position du front polaire pendant cette période de temps.

La déglaciation est caractérisée par des valeurs de précipitation plus importantes, 960−390+80 mman−1 autour de 9000 ans BP, décroissant jusqu'à 860340+160 mman−1 vers 8000 ans BP puis rejoignant le niveau actuel de 800 mm an−1.

1 Introduction

Thanks to the recent NGRIP core completing the GRIP and GISP2 cores, 25 Dansgaard/Oeschger (D/O) events have now been enumerated in Greenland records [36]. Several of these millennial-scale climatic events have reached air-temperature changes greater than 10 °C in the Greenland sites [25,30,31]; they are also detectable in the water temperature and sediments of North-Atlantic and Mediterranean marine cores, demonstrating their great link between the atmosphere and the ocean, at least in the Northern Hemisphere [5–7]. However, their consequences on the continent are still not well known, because of the scarcity of well-preserved and well-dated records. What we know is that the ocean–atmosphere system was closely coupled as far as in the central Mediterranean region, where the vegetation, reconstructed thanks to pollen lake records, reacted very quickly and sensitively to these abrupt climate changes [1,16]. The D/O events had also a large impact on the atmospheric circulation, and on the monsoon intensity specifically, as far as eastern China, as the remarkable Hulu Cave δ18O records have shown [50]. But up to now, there has been very little information about their impact on the European continent. Here we present results from three kinds of palaeoenvironmental records that give an initial pattern of the D/O impact on western Europe: (1) a pollen assemblage record from the ODP 976 marine core, in the Alboran Sea, that is directly controlled by the nearby continental vegetation of Southern Spain and North Africa; (2) a δ13C record of stalagmites from the Villars Cave, southwestern France, that recorded past vegetation and soil activity and gives an accurate chronology of several of these events with four hiatuses; (3) a loess-deposit sequence whose δ13C inverse modelling brings a useful rainfall reconstruction around the Nußloch area, in western Germany, with an irregular temporal resolution.

2 Studied sites

The three studied records are located along a SW–NE profile from southern Spain (core ODP 976) to western Germany (Nußloch loess sequence), the Villars Cave in South-France being situated between theses two sites (Fig. 1).

Fig. 1

Localisation of the studied sites: Alboran Sea marine core ODP 976 (36°12′N, 4°18′W, 1108 m water depth), Villars Cave (southwestern France; 45°30′N, 0°50′E, elevation 175 m), Nußloch loess section (49°19′N, 8°43′E, 221 m asl).

Position des sites étudiés, du sud au nord : carotte marine ODP 976 de la mer d'Alboran (36°12′N, 4°18′W, 1,108 m profondeur), grotte de Villars (Sud-Ouest de la France; 45°30′N, 0°50′E, altitude 175 m), Nußloch loess section (49°19′N, 8°43′E, altitude 221 m).

Ocean Drilling Program (ODP) Site 976 has been drilled in the Alboran Sea, the westernmost section of the Mediterranean Sea (36°12N, 4°18W, 1108-m water depth), close to the Strait of Gibraltar. This record provides new insight into the evolution of terrestrial palaeoenvironments during the past 90 000 years, which enables a direct correlation with climates in the larger northern-Europe area.

The Villars Cave (southwestern France; 45°30′N, 0°50′E, elevation: 175 m), developed in Bajocian limestone at depths from 10 to 40 m from the surface. Its environmental setting is detailed in [3,11]. The four stalagmites studied here (Vil-stm9, Vil-stm10, Vil-stm11, Vil-stm27) developed in the deepest part of the cave, about 30 to 40 m below the surface. They were all dated by uranium–thorium method by thermal ionisation mass spectrometry (TIMS UTh), and stable isotope profiles were performed along their axis with sampling spaced at about 3 mm, giving a high-resolution continental record of the last 82 kyr [12].

The Nußloch loess section is located on the right bank of the Upper Rhine graben in southern Germany (49°19′N, 8°43′E, 221 m asl, 800 mm yr−1 of mean annual precipitation). The geomorphological setting is characterized by a wide alluvial plain and an abrupt increase in elevation to the east. During the Weichselian Pleniglacial, loess accumulated on the slope and the adjacent hill country, producing a series of large loess dunes. The last interglacial/glacial climatic cycle is represented by the uppermost 20 m of loess, with a temporal resolution higher than 1 m for 1000 years for the typical loess of the pleniglacial Weichselian. An erosion hiatus eliminates the 33.846 kyr BP phase out of the deposit. This sequence constitutes the key-section of western Europe [19,20,40].

3 Methods and chronologies

3.1 Dating methods

The chronology is a key point in the interpretation and the comparison of several palaeoclimatic records. The three records that are presented here have great differences in their timescale accuracy and in their palaeoclimate signals. Taken alone, none of these archive types possesses all the best ‘qualities’ (i.e., best chronology, best climatic signal); they are however complementary: while speleothems, thanks to UTh ages, have the most accurate chronology, they do not have the easiest interpretable palaeoclimate signal, which is the opposite for the pollen record. At the same time, hiatuses due to too cold and dry periods that are marked in the speleothems by a thin brown layer, are extensively recorded in the loess deposits.

3.1.1 ODP 976 chronology

The chronology of the ODP site 976 is first based on 13 accelerator mass spectrometry (AMS) radiocarbon ages [8]. This chronology has been completed by other control points obtained according to the correlation between the Oxygen Isotope curve of ODP site 976 [8,14] and the reference curve of the core MD 95-2042 off Portugal [43–45], especially for the older part of the series. Because of the similarity in the pollen curves and the Greenland records [15,36], after being established by linear interpolation between the dates and oxygen isotope control points, the age model has been refined by correlation between pollen and North Grip records down to 120 000 yr [36].

3.1.2 Villars cave speleothems chronology

Sixty-nine UTh ages by thermal ionization mass spectrometry (TIMS) have been performed on the four studied stalagmites from the Villars Cave ([12] and unpublished data; Fig. 2). Details of the method are explained in Genty et al. [12,13]. When necessary, a 232Th correction has been performed in order to take into account the detrital origin of the 230Th. Final error is given at 2σ, which is, for most samples, lower than 2% of the age (see error bars at the bottom of Fig. 2) generating among the best-dated continental sequence.

3.1.3 Nußloch loess deposit chronology

The chronology of the Nußloch loess sequence is based on 22 14C datings performed on bulk organic matter (15), wood (4), bones (1), carbonate gastropod shells (1), conchioline (1) ranging from 3 to 42 kyrs BP and on 21 IRSL/OSL datings ranging from 19 to 122 kyr [21,32]. Ages younger than 20 kyr BP were calibrated using the Calib 4.18 program [47] and older ages were ‘corrected’ according to the Kitagawa and van der Plicht study [29]. The chronology used here was established by linear interpolation between median values of the dating ranges.

3.2 Measurement methods and proxy significance

3.2.1 ODP 976 palynological method

Pollen process follow a standard protocol already described in Combourieu Nebout et al. [8]. Palaeoenvironmental interpretation of the down-core pollen assemblage fluctuations is based on the assumption that the primary pollen contribution to Alboran-Sea sediments comes from west-Mediterranean borderlands. Modern environments range from a thermo-Mediterranean belt with Olea, Pistacia and some semi-desert representatives (Artemisia, Chenopodiaceae, Ephedra) to a meso-Mediterranean belt, represented by a sclerophyllous oak forest then a humid-temperate oak forest (Euro-Siberian trees as Quercus, Betula... , and Ericaceae), and a supra-Mediterranean belt, with a cold-temperate coniferous forest (Pinus, Abies, Cedrus) at the higher altitudes [37,39] (Fig. 1). Today, the area experiences long dry summers and mild rainy winters. The pollen spectrum, based on 120 pollen taxa, ranges from semi-desert to mountain deciduous and coniferous forest.

Interpretation of the pollen spectrum follows the modern climatic-plant relationships in Eurasia and northern Africa [38,52]. Here, we present the variations of pollen percentages of two main associations (Fig. 2):

  • – the temperate association composed of Euro-Siberian trees like Quercus, Fagus, Carpinus, Corylus, Alnus, Betula, Tilia, Ulmus... , associated with Ericaceae, which reflects warmer and moist climate characteristic of interstadials;
  • – the steppe to semi-desert association, composed of Artemisia, Amaranthaceae-Chenopodiaceae and Ephedra, whose high representations in pollen spectrum indicate dry and cold climatic conditions of stadials.

3.2.2 Stable isotope measurements on Villars cave speleothems

Samples have been taken every 1 to 3 mm, along the centre of the growth axis and analysed with the standard procedures on a VG OPTIMA leading to analytical error ±0.08 [12]. The isotopic equilibrium, checked in two ways [23], leads to the conclusion that great kinetic fractionation is not likely, and, at places where it may have occurred (R2>0.5 and slight enrichment toward the edges), the overall δ13C signal goes in the same direction as the climate changes.

Speleothem calcite δ18O is controlled by several environmental factors that can act in opposite ways, leading to complex interpretations where rainfall amount or temperature intervene [4,28,46]. On the other hand, the speleothem δ13C changes can be interpreted as evidence of vegetation density and soil activity changes for the following reasons: (1) the main sources of the carbon in speleothem calcite is the soil CO2 as demonstrated in former studies [13]; (2) δ13C decreases when the soil activity and the vegetation respiration is efficient and increases when the climate degradation has slowed down this soil and vegetation activity. Changes in the C4/C3 plant ratio suggested in other studies [10] cannot be invoked here as, until now, no evidence of C4 plants has been found in this area during the last Glacial.

3.2.3 Nußloch Loess deposit isotope measurements and BIOME 4 modelling methods

For the isotopic study, the Nußloch sequence was sampled by taking 1 to 2 g of sediments at intervals of 510 cm and analysing bulk carbon isotopic composition. The δ13C values vary from −25.72 to 23.52±0.15. The absence of pedogenesis and the dry glacial environment inducing the degradation of organic matter without distortion of the isotopic signal make typical loess suitable for an organic geochemical study. The carbon isotopic composition (δ13C) of loess organic matter reflects nicely the original isotopic signature of vegetation and informs about palaeoenvironmental conditions. In western Europe, δ13C values ranging from −23 to 26 indicate that, during the last climatic cycle, the vegetation maintained a unique presence of plants following the C3 photosynthetic pathway [20]. Thus δ13C variations were interpreted as changes in environmental conditions, related, firstly, to concentration and isotopic composition of atmospheric CO2 and water availability and, secondly, to temperature, soil type and texture, and insolation.

Vegetation models introduce such great complexity by considering first- and second-order parameters. Recently, one of them, Biome3 [22] has been modified into Biome4, including a carbon isotopes routine [19,26]. This model simulates the biome that would growth in the given climatic conditions, some production parameters (net and growth primary production, leaf area index, etc.) and the mean δ13C of the plant community. This version of Biome4 is used in an inverse mode to reconstruct the most probable climate, i.e., starting from the outputs, to infer some of the inputs of the model. We follow the procedure of Guiot et al. [17] based on a Monte-Carlo sampling of monthly precipitation, and monthly temperature in a pre-defined range (Metropolis–Hastings algorithm, [18]). The output biome and the plant NPP (Net primary Production) are constrained by temperature (annual distribution and absolute value) as well as by precipitation regime. The carbon isotopic composition, greatly dependent on precipitation, narrows the precipitation-variation range.

4 Results and discussion

Thanks to the different nature of the three studied records, we obtain a coherent and complementary pattern of the western-Europe continental environment variations of the last 90 kyr (Fig. 2). While pollen assemblages from the ODP 976 core display obviously well the cold episodes (Heinrich events H1 to H6) thanks to semi-desert taxa, the warm and more humid periods are clearly recorded in the temperate taxa of the same core and also in the Villars stalagmite growth periods and δ13C (D/O #7 to D/O #20). During extremely cold episodes, which broadly coincide with MIS 4 and MIS 2, speleothem deposition stopped because of vegetation and soil activity slowing down and/or a stop in the seepage, which is not the case for the loess deposit at Nußloch nor for the ODP 976 core, where climate variations are still visible during these stages. Finally, thanks to inverse modelling, the Nußloch loess deposit data show by taking off global components (CO2 and δ13C), the high variability of the rainfalls during this period.

4.1 ODP 976 palynological record

High-amplitude shifts during the last 90 000 years are displayed in the pollen record, which indicates quasi-cyclic oscillations of environmental conditions in the westernmost Mediterranean (Fig. 2). High down-core variability of pollen assemblages reflects rhythmic alterations from temperate forest (mainly Quercus) very similar to today's vegetation in the mountain ranges of the western Mediterranean borderlands to semi-desert vegetation as observed today in North Africa and southwestern Europe [50]. These vegetation changes are indicative of repetitive and abrupt climatic changes in the Alboran Sea's borderlands between warm/moist and cold/dry conditions. The pattern of temperate taxa oscillations mirrors the D/O climatic oscillations seen in the Greenland ice core records [1,34] in that 20 increases in temperate vegetation taxa abundance are correlated with the D/O interstadials 1 through 21 (Fig. 2). The colder stadials are characterized by abundant increases in semi-desert taxa. Six of them are particularly extreme, with the lowest occurrences of temperate taxa and the highest representation of semi-desert represented by Artemisia, Ephedra and Chenopodiaceae. Such associations, characteristic of cool steppe to cool desert [48], suggest enhanced aridity on the adjacent westernmost Mediterranean borderlands. Combined marine-terrestrial records at Site 976 [8] and in other cores in the western Mediterranean and off Portugal [42] and results obtained in continental sites [52] have already shown that during the last glacial, aridity increased noticeably in the westernmost Mediterranean borderlands far beyond its present state, in association with a clear cooling of the sea surface waters associated to Heinrich events. Then, the six dry terrestrial episodes recorded in ODP Site 976 seem tightly linked to Heinrich cold events H1 to H6 [8,27].

4.2 Villars Cave δ13C record

The Villars δ13C record, made by the stack of four stalagmites, possesses four hiatuses, which correspond to extreme climatic events that crossed thresholds and prevented stalagmite growth (Fig. 2). Hiatuses D2 and D4 that occurred at 67.4–61.2 kyr (‘Villars cold phase’) and 31.6–15.9 kyr, respectively, are the consequence of a very cold climate that progressively occurred at Villars, as the regular δ13C towards these hiatuses and also the steady increase that is observed just after show. They coincide with extremely low temperatures in the NGRIP record (low δ18O values; Fig. 2), and to cold marine stages MIS 4 and MIS 2. Because D1 (78.8–75.5 kyr) and D3 (55.7–51.8 kyr) hiatuses occurred during stable or low δ13C and δ18O levels, their cause is possibly due to flooding in the Villars Cave, due to very humid periods or, like former hiatuses, to cold periods that seem to correlate with NGRIP δ18O and OPD 976 pollen record; in this case, the transition toward the cold episode was too rapid to have been recorded in the δ13C record.

In the Villars Cave's δ13C record, D/O events #8, #12, #16, #17, #18 and #19 are particularly visible with typically 24 variations (Fig. 2; [12]). A comparison with δ18O Greenland records allows their identification, low δ13C values corresponding to high ice δ18O and thus to warmer conditions and, at least in Greenland, to higher precipitation (Fig. 2). Most D/O ages from Greenland records are within the 2σ error limits of our chronology. However, the best fit occurs with the older GISP2 chronology [15] and significant differences are observed for D/O #12, #19 and #20 with the new ss09sea Greenland chronology [24]. The fact that we made double measurements of all the UTh ages on a second stalagmite for the 52–31-kyr period, and that we found similar growth curve and similar δ13C signal, strengthen our UTh-based chronology, at least for this period (orange and red curves, Fig. 2d). D/O #8 displays a regular δ13C decrease from 39.4 to 36.2 ka, whose shape is very similar to another speleothem record from China (Hulu Cave record; [51]). It is the last important warm event before the pleniglacial period recorded in Villars stalagmites, but is already part of a climatic decline trend that started about 45 kyr ago up to the stop of the stalagmites Vil-stm9 and Vil-stm27 about 31.7±0.1 kyr ago (Fig. 2d). Cold phases surrounding the D/O events are in agreement with the semi-desert pollen taxa peaks found in the ODP 976 record, and likely correspond to the Heinrich events: H4 and H5 are well visible with high δ13C values (isotope reversed scale on the figure), H6 corresponds to the D2 hiatus and H2 and H3 to the D4 hiatus. The most pronounced warm event of the studied period is certainly the D/O #12, which reached its optimum (lower δ13C: −10.5 to 11.2), close to present-day values, at 45.3±0.4 kyr. The very high growth rate observed during this stage demonstrates also the high humidity of this period. Finally, the H5-D/O#12 transition, marked by a δ13C decrease by more than 4, appears to have occurred in less than 1 kyr.

The most recent part of the Villars record displays clearly the last deglaciation: the beginning of the warming is dated to 15.9±0.35 kyr and the Vil-stm11 δ13C decreased by more than 5 until a climatic optimum dated at 9.9±0.17 kyr. The Younger–Dryas cooling is well marked by a δ13C increase and a growth rate slowing down.

4.3 Nußloch Loess deposit isotope record and rainfall modelling

The inverse modelling simulates palaeoprecipitation ranging from 280−60+120 to 460−160+80 mmyr−1 along the glaciation well below the modern value of about 800 mm yr−1. This aridity is in agreement with pedological and malacological bio-indicators [2,35]. The palaeoprecipitation reconstruction obtained by inverse modelling can be divided into three parts. Prior to 28 kyr BP, the precipitation varies from 280−60+120 mmyr−1 to 460−160+80 mmyr−1, and is relatively constant in the ca 28–15-kyr BP interval. It finally increases during the deglaciation.

The period before 28 kyr BP presents large oscillations of about 150–200 mm yr−1. These are very similar and contemporaneous to some D/O events, as recorded by Greenland ice-cores (NGRIP). The D/O 14, 15, 16, 18 as well as 3 or 4 are particularly well recognizable. The D/O 14, which is associated with a temperature shift of 10 to 12 °C in Greenland [9,30,33], seems to be characterized by a precipitation increase of 200 mm yr−1, i.e. + ca 70%.

The similarity expressed in both ice cores and Nußloch underlines a narrow link between precipitation regime, polar and continental ice-sheets size, sea level and consequently the polar front position. At the first order, ice δ18O changes in ice are linked to those of the atmospheric temperature and depend on the ice-sheet size: the smaller the ice sheet, the more positive the δ18O and the higher the temperature. This pulls the polar front to the North Atlantic high latitudes, favouring a mild oceanic regime over western Europe and over the Nußloch site. Alternatively, during stadials, the polar front moved southwards [41], limiting the entrance of wet air over Europe.

A nearly constant level of precipitation at around 280 mm yr−1 is observed between 28 and 17 kyr BP during full glacial conditions. At this time, the English Channel was mostly emerged (sea level ca 90 m below present conditions, [49]). Such a geographical configuration is consistent with a great reduction in the entrance of the wet oceanic air fluxes. Nußloch is thus less sensitive to the changes of the polar front position during this time interval.

The deglaciation and Holocene onset are characterized by higher precipitation, reaching 960−390+80 mmyr−1 at ca 9000 years BP and decreasing to 860−340+160 mmyr−1 around 8000 yr BP, framing the present 800 mm yr−1 (Fig. 2b). The steady increase of the precipitation associated with the deglaciation is interrupted between 14 and 12.5 kyr BP, which it is tempting to relate to the Younger Dryas.

5 Conclusion

This study is an example of a multi-proxy approach to better describe and understand the abrupt millennial scale events that occurred in western Europe during the last 90 kyr. It shows the limits and complementarities of the proxy, each one bringing independent information that can be confirmed or specified by another one. For example, palaeoprecipitations have been reconstructed over western Germany thanks to loess deposits and to modelling; pollen assemblages of a marine core and especially semi-desert and arboreal assemblages have allowed the identification of cold and warm stages, leading to hypotheses on the past vegetation and on the past atmospheric circulation; finally, thanks to speleothems, an accurate absolute chronology on the past vegetation and temperature changes can be established, which is the indispensable step toward any lead and lag comparison with very distant records. These three records show that, during a transition toward a D/O event, the precipitation increases, semi-desert plants are replaced by trees, and vegetation and soil activity increase. From UTh age interpolations, all these changes appear to have been rapid, i.e. from a few decades to a few hundred years. As an example, 46.6±0.5 kyr ago, one of the most important transitions toward a warm climatic event occurred (D/O event #12). It is characterized by a rainfall increase of about 200 mm yr−1 (from 250 to 450 mm yr−1) in the Rhine Valley, by a drastic drop in semi-desert taxa in the Western Mediterranean basin, demonstrating the disappearance of the steppe and by a clear decrease in the δ13C of the speleothem calcite and by a high calcite growth rate, giving evidence of a return to an efficient plant and soil activity in southwestern France. The whole transition, which almost reached the amplitude of a full glacial–interglacial transition, occurred in less than 1.2 kyr between 46.6 kyr and 45.4 kyr; moreover, its last part, whose amplitude covers almost half of the total variation, occurred extremely rapidly in less than the 2-σ age error (i.e. <450 yr) between 45.5 and 45.2 kyr (Fig. 2c and d). Most of the other D/O transitions show a similar rapidity (D/O 9, 10, 19 and 20) and the future challenge will be, first, to get a better accuracy for these climate transitions, thanks, for example, to annual laminae in speleothems in Villars cave and others, and also to get a more complete climate quantification (P,T), in a larger spatial cover, thanks to loess and pollen assemblages. Further analyses of other similar records will permit to obtain a global idea of the timing and amplitude of the D/O events during the last glacial period, to check any time lag between distant records (i.e. south/north, coastal/continental) and then to test hypothetical mechanisms with climatic models.

Acknowledgements

We thank CNRS (France) programs ECLIPSE, PNEDC, the French–Italian Galilee project and EC programs for their financial support. D.G. thanks Hubert Versaveau for access to the Villars Cave, and T. Baritaud for his help in sampling and caving. Vil-stm10 U/Th dates were performed by M. Bar-Matthews (Geological Survey of Israel). Vil-stm9 U/Th analyses were performed by Mabs Gilmour; all other U/Th were performed by B. Ghaleb and L. Isabello at the GEOTOP. N.C.N. thanks the Ocean Drilling Program for giving us the samples for the study. The 14C AMS dates have been measured at the Leibniz Laboratory at Kiel University under the supervision of M. Nadeau and P. Grootes. We acknowledge M.H. Castera and J.-P. Cazet for their efficient technical assistance. C.H. would like to thank Denis-Didier Rousseau and Pierre Antoine, leaders of the EC Program (BIMACEL), and Joël Guiot, who opened the door of the inverse modelling world.


References

[1] J.R.M. Allen; U. Brandt; A. Brauer; H.-W. Hubberten; B. Huntley; J. Keller; M. Kraml; A. Mackensen; J. Mingram; J.F.W. Negendank; N.R. Nowaczyk; H. Oberhänsli; W.A. Watts; S. Wulf; B. Zolitschka Rapid environmental changes in southern Europe during the last glacial period, Nature, Volume 400 (1999), pp. 740-743

[2] P. Antoine; D.D. Rousseau; L. Zöller; A. Lang; A.-V. Munaut; C. Hatté; M.R. Fontugne High-resolution of the last interglacial–glacial cycle in the loess paleosol sequences of Nußloch (Upper Rhine Area, Germany), Quat. Int., Volume 76–77 (2001), pp. 211-229

[3] A. Baker; D. Genty; I. Fairchild Hydrological characterization of stalagmite dripwaters at Grotte de Villars, Dordogne, by the analysis of inorganic species and luminescent organic matter, Hydrol. Earth Syst. Sci., Volume 4 (2000), pp. 439-449

[4] M. Bar-Matthews; A. Ayalon; A. Kaufman Late Quaternary paleoclimate in the eastern Mediterranean region from stable isotope analysis of speleothems at Soreq Cave, Israel, Quat. Res., Volume 47 (1997), pp. 155-168

[5] G. Bond; W. Broecker; S. Johnsen; J. McManus; L. Labeyrie; J. Jouzel; G. Bonani Correlations between climate records from North Atlantic sediments and Greenland ice, Nature, Volume 365 (1993)

[6] I. Cacho; J.O. Grimalt; C. Pelejero; M. Canals; F. Sierro; J.A. Flores; N. Shackelton Dansgaard–Oeschger and Heinrich event imprints in Alboran Sea paleotemperatures, Paleoceanography, Volume 14 (1999), pp. 698-705

[7] O. Cayre; Y. Lancelot; E. Vincent Paleoceanographic reconstructions from planktonic foraminifera of the Iberian Margin: Temperature, salinity, and Heinrich events, Paleoceanography, Volume 14 (1999), pp. 384-396

[8] N. Combourieu Nebout; J.-L. Turon; R. Zahn; L. Capotondi; L. Londeix; K. Pahnke Enhanced aridity and atmospheric high-pressure stability over the western Mediterranean during the North Atlantic cold event of the past 50 kyr, Geology, Volume 30 (2002) no. 10, pp. 863-866

[9] D. Dahl-Jensen; K. Mosegaard; G.D. Gunderstrup; G.D. Clow; S.J. Johnsen; A.W. Hansen; N. Balling Past temperatures directly from the Greenland ice sheet, Science, Volume 282 (1998), pp. 268-271

[10] J.A. Dorale; R.L. Edwards; E. Ito; L.A. Gonzales Climate and vegetation history of the Midcontinent from 75 to 25 ka: A speleothem record from Crevice Cave, Missouri, USA, Science, Volume 282 (1998), pp. 1871-1874

[11] D. Genty; A. Baker; B. Vokal Intra- and inter-annual growth rate of modern stalagmites, Chem. Geol., Volume 176 (2001), pp. 191-212

[12] D. Genty; D. Blamart; R. Ouahdi; M. Gilmour; A. Baker; J. Jouzel; S. Van-Exter Precise dating of Dansgaard–Oeschger climate oscillations in western Europe from stalagmite data, Nature, Volume 421 (2003), pp. 833-837

[13] D. Genty; M. Massault Carbon transfer dynamics from bomb-14C and δ13C time series of a laminated stalagmite from SW France – Modeling and comparison with other stalagmite records, Geochim. Cosmochim. Acta, Volume 63 (1999) no. 10, pp. 1537-1548

[14] R. von Grafenstein; R. Zahn; R. Tiedemann; A. Murat Planktonic 18O records at Sites 976 and 977, Alboran Sea: stratigraphy, forcing, and paleoceanographic implications, Proc. ODP, Sci. Results, College Station, TX (Ocean Drilling Program), Volume 161 (1999), pp. 469-479

[15] P.M. Grootes; M. Stuiver; J.W.C. White; S. Johnsen; J. Jouzel Comparison of oxygen isotope records from the GISP2 and GRIP Greenland ice cores, Nature, Volume 366 (1993), pp. 552-554

[16] J. Guiot; J.-L. de Beaulieu; R. Cheddadi; F. David; P. Ponel; M. Reille The climate in Western Europe during the last Glacial/Interglacial cycle derived from pollen and insect remains, Palaeogeogr. Palaeoclimatol. Palaeoecol., Volume 103 (1993), pp. 73-93

[17] J. Guiot; F. Torre; D. Jolly; O. Peyron; J.-J. Boreux; R. Cheddadi Inverse vegetation modeling by Monte Carlo sampling to reconstruct palaeoclimates under changed precipitation seasonality and CO2 conditions: application to glacial climate in Mediterranean region, Ecol. Model., Volume 127 (2000), pp. 119-140

[18] W.K. Hastings Monte-Carlo sampling methods using Markov chains and their application, Biometrika, Volume 57 (1970), pp. 97-109

[19] C. Hatté, J. Guiot, Paleoprecipitation reconstruction by inverse modeling using the isotopic signal of loess organic matter: application to the Nußloch loess sequence (Rhine Valley, Germany), Clim. Dynam. (in press)

[20] C. Hatté; M.R. Fontugne; D.-D. Rousseau; P. Antoine; N. Tisnérat-Laborde δ13C variations of loess organic matter as a record of the vegetation response to climatic changes during the Weichselian, Geology, Volume 26 (1998), pp. 583-586

[21] C. Hatté; L.C.R. Pessenda; A. Lang; M. Paterne Development of an accurate and reliable 14C chronology for loess sequences, Application to the loess sequence of Nußloch (Rhine valley, Germany), Radiocarbon, Volume 43 (2001) no. 2B, pp. 611-618

[22] A. Haxeltine; I.C. Prentice BIOME 3: An equilibrium terrestrial biosphere model based on ecophysiological constraints, resource availability and competition among plant functional types, Global Biogeochem. Cycles, Volume 10 (1996) no. 4, pp. 693-709

[23] C.H. Hendy The isotopic geochemistry of speleothems – I. The calculation of the effects of different modes of formation on the isotopic composition of speleothems and their applicability as palaeoclimatic indicators, Geochim. Cosmochim. Acta, Volume 35 (1971), pp. 801-824

[24] S.J. Johnsen; D. Dahl-Jensen; N. Gundestrup; J.P. Steffensen; H.B. Clausen; H. Miller; V. Masson-Delmotte; A.E. Sveinbjörnsdottir; J. White Oxygen isotope and palaeotemperature records from six Greenland ice-core stations: Camp Century, Dye-3, GRIP, GIPS2, Renland and NorthGRIP, J. Quat. Sci., Volume 16 (2001) no. 4, pp. 229-307

[25] J. Jouzel Towards a calibration of the isotopic paleothermometer, Science, Volume 286 (1999) no. 5441, pp. 910-911

[26] J.O. Kaplan; I.C. Prentice; N. Buchmann The stable carbon isotope composition of the terrestrial biosphere: modeling at scales from the leaf to the globe, Global Biogeochem. Cycles, Volume 16 (2002) no. 4, p. 1060

[27] M. Kageyama; N. Combourieu Nebout; P. Sepulchre; O. Peyron; G. Krinner; G. Ramstein; J.-P. Cazet The Last Glacial Maximum and Heinrich Event 1 in terms of climate and vegetation around the Alboran Sea: a preliminary model-data comparison, C. R. Geoscience, Volume 337 (2005) | DOI

[28] S.-T. Kim; J.R. O'Neil Equilibrium and nonequilibrium oxygen isotope effect in synthetic carbonates, Geochim. Cosmochim. Acta, Volume 61 (1997), pp. 3461-3475

[29] H. Kitagawa; J. van der Plicht Atmospheric Radiocarbon Calibration to 45 000 yr BP: Late Glacial fluctuations and cosmogenic isotope production, Science, Volume 279 (1998), pp. 1187-1190

[30] A. Landais; N. Caillon; J. Severinghaus; J.M. Barnola; C. Goujon; J. Jouzel; V. Masson-Delmotte Isotopic measurements of air trapped in ice to quantify temperature changes, C. R. Geoscience, Volume 336 (2004) no. 11, pp. 963-970

[31] A. Landais; J. Chappellaz; M. Delmotte; J. Jouzel; T. Blunier; C. Bourg; N. Caillon; S. Cherrier; B. Malaize; V. Masson-Delmotte; D. Raynaud; J. Schwander; J.-P. Steffensen A tentative reconstruction of the last interglacial and glacial inception in Greenland based on new gas measurements in the Greenland Ice Core Project (GRIP) ice core, J. Geophys. Res. – Atmos., Volume 108 (2003) no. D18

[32] A. Lang; C. Hatté; D.-D. Rousseau; P. Antoine; M.R. Fontugne; L. Zöller; U. Hambach High-resolution chronologies for loess: comparing AMS 14C and optical dating results, Quat. Sci. Rev., Volume 22 (2003), pp. 953-959

[33] C. Lang; M. Leuenberger; J. Schwander; S. Johnsen 16 °C Rapid temperature variation in central Greenland 70 000 years ago, Science, Volume 286 (1999), pp. 934-937

[34] D.A. Meese; A.J. Gow; R.B. Alley; G.A. Zielensky; P.M. Grootes The Greenland Ice project 2 depth age scale, Methods and results, J. Geophys. Res., Volume 12 (1997), pp. 411-423

[35] O. Moine; D.-D. Rousseau; P. Antoine; C. Hatté Mise en évidence d'événements climatiques rapides par les faunes de mollusques terrestres des loess weichséliens de Nußloch (Allemagne), Quaternaire, Volume 13 (2002) no. 3–4, pp. 209-217

[36] North GRIP members, High-resolution climate record of the Northern Hemisphere reaching into the Last Interglacial period, Nature, Volume 43 (2004), pp. 147-151

[37] P. Ozenda Sur les étages de végétation dans les montagnes du bassin méditerranéen, Doc. Cartogr. Écol., Volume 16 (1975), pp. 1-32

[38] O. Peyron; J. Guiot; R. Cheddadi; P. Tarasov; M. Reille; J.-L. de Beaulieu; S. Bottema; V. Andrieu Climatic reconstruction in Europe for 18 000 yr BP from pollen data, Quat. Res., Volume 49 (1998), pp. 183-196

[39] S. Rivas Martinez Étages bioclimatiques, secteurs chorologiques et séries de végétation de l'Espagne méditerranéenne, Ecol. Méditerr., Volume 8 (1982), pp. 275-288

[40] D.-D. Rousseau; P. Antoine; C. Hatté; A. Lang; L. Zöller; M.R. Fontugne; D. Ben Othman; J.-M. Luck; O. Moine; M. Labonne; I. Bentaleb; D. Jolly Abrupt millennial climatic changes from Western European eolian records during the last glaciation, Quat. Sci. Rev., Volume 21 (2002), pp. 1577-1582

[41] W. Ruddiman; A. McIntyre The North-Atlantic Ocean during the last deglaciation, Palaeogeogr. Palaeoclimatol. Palaeoecol., Volume 35 (1981) no. 2–4, pp. 145-214

[42] M.F. Sanchez Goni; I. Cacho; J.-L. Turon; J. Guiot; F.J. Sierro; J.-P. Peypouquet; J.-O. Grimalt; N.J. Shackleton Synchroneity between marine and terrestrial responses to millennial scale climatic variability during the last glacial period in the Mediterranean region, Clim. Dynam., Volume 19 (2002), pp. 95-105

[43] N.J. Shackleton; R.G. Fairbanks; T. Chiu; F. Parrenin Absolute calibration of the Greenland time scale: implications for Antarctic time scales and for δ14C, Quat. Sci. Rev., Volume 23 (2004), pp. 1513-1522

[44] N.J. Shackleton; M. Hall; E. Vincent Phase relationships between millennial-scale events 64 000–24 000 years ago, Paleoceanography, Volume 15 (2000) no. 6, pp. 565-569

[45] N.J. Shackleton; M.F. Sanchez Goñi; D. Pailler; Y. Lancelot Marine isotope substage 5e and the Eemian interglacial, Global Planet Change, Volume 36 (2003), pp. 151-155

[46] C. Spötl; A. Mangini Stalagmite from the Austrian Alps reveals Dansgaard–Oeschger events during isotope stage 3: implications for the absolute chronology of Greenland ice cores, Earth Planet Sci. Lett. 203 (2002), pp. 507-518

[47] M. Stuiver; P.J. Reimer; E. Bard; W. Beck; G.S. Burr; K.A. Hughen; B. Kromer; F.G. Mc Cormac; J. van der Plicht; M. Spurk INTCAL98 Radiocarbon age calibration, 24 000 cal BP, Radiocarbon, Volume 40 (1998), pp. 1041-1083

[48] P.E. Tarasov; R. Cheddadi; J. Guiot; S. Bottema; O. Peyron; J. Belmonte; V. Ruiz-Sanchez; F. Saadi; S. Brewer A method to determine warm and cool steppe biomes from pollen data; application to the Mediterranean and Kazakhstan regions, J. Quat. Sci., Volume 13 (1998), pp. 335-344

[49] C. Waelbroeck; L. Labeyrie; E. Michel; J.-C. Duplessy; J.F. McManus; K. Lambeck; E. Balbon; M. Labracherie Sea-level and deep water temperature changes derived from benthic foraminifera isotopic records, Quat. Sci. Rev., Volume 21 (2002), pp. 295-305

[50] W.A. Walter Die Vegetation Ost-Europas, Nord- und Zentralasians, Gustav Fisher Verlag, Stuttgart Germany, 1974 (452 p)

[51] Y.J. Wang; H. Cheng; R.L. Edwards; Z.S. An; J.Y. Wu; C.C. Shen; J.A. Dorale A high-resolution absolute-dated Late Pleistocene Monsoon Record from Hulu Cave, China, Science, Volume 294 (2001), pp. 2345-2348

[52] W.A. Watts; J.R.M. Allen; B. Huntley Vegetation history and palaeoclimate of the last glacial period at Lago Grande di Monticchio, Southern Italy, Quat. Sci. Rev., Volume 15 (1996), pp. 133-153


Comments - Policy