<figure><img src="image_1.png" /></figure>

Check for
updates

# Pollen from the Deep-Sea: A Breakthrough in the Mystery of the Ice Ages

María F. Sánchez Goñi¹,²*, Stéphanie Desprat¹,², William J. Fletcher³,
César Morales-Molino¹,², Filipa Naughton⁴,⁵, Dulce Oliveira¹,²,⁴, Dunia H. Urrego⁶ and
Coralie Zorzi⁷

¹École Pratique des Hautes Études, EPHE PSL University, Paris, France, ²Environnements et Paléoenvironnements Océaniques et Continentaux, UMR 5805, Université de Bordeaux, Pessac, France, ³Quaternary Environments and Geoarchaeology, Department of Geography, School of Environment, Education and Development, The University of Manchester, Manchester, United Kingdom, ⁴Instituto Português do Mar e da Atmosfera, Portuguese Institute of Sea and Atmosphere, Lisbon, Portugal, ⁵Center of Marine Sciences, Algarve University, Faro, Portugal, ⁶College of Life and Environmental Sciences, University of Exeter, Exeter, United Kingdom, ⁷GEOTOP, Université du Québec à Montréal, Montreal, QC, Canada

## OPEN ACCESS

**Edited by:**
Valentí Rull,
Instituto de Ciencias de la Tierra
Jaume Almera (CSIC), Spain

**Reviewed by:**
Francesc Burjachs,
Institució Catalana de Recerca i
Estudis Avançats (ICREA), Spain
Jean-Pierre Suc,
UMR 7193, Institut des Sciences
de la Terre Paris (ISTEP), France

**\*Correspondence:**
María F. Sánchez Goñi
maría.sanchez-goni@u-bordeaux.fr

**Specialty section:**
This article was submitted to
Agroecology and Land Use Systems,
a section of the journal
Frontiers in Plant Science

**Received:** 12 July 2017
**Accepted:** 09 January 2018

**Published:** 26 January 2018

**Citation:**
Sánchez Goñi MF, Desprat S, Fletcher WJ, Morales-Molino C, Naughton F, Oliveira D, Urrego DH and Zorzi C (2018) Pollen from the Deep-Sea: A Breakthrough in the Mystery of the Ice Ages. *Front. Plant Sci.* 9:38. doi: 10.3389/fpls.2018.00038

Pollen from deep-sea sedimentary sequences provides an integrated regional reconstruction of vegetation and climate (temperature, precipitation, and seasonality) on the adjacent continent. More importantly, the direct correlation of pollen, marine and ice indicators allows comparison of the atmospheric climatic changes that have affected the continent with the response of the Earth's other reservoirs, i.e., the oceans and cryosphere, without any chronological uncertainty. The study of long continuous pollen records from the European margin has revealed a changing and complex interplay between European climate, North Atlantic sea surface temperatures (SSTs), ice growth and decay, and high- and low-latitude forcing at orbital and millennial timescales. These records have shown that the amplitude of the last five terrestrial interglacials was similar above 40°N, while below 40°N their magnitude differed due to precession-modulated changes in seasonality and, particularly, winter precipitation. These records also showed that vegetation response was in dynamic equilibrium with rapid climate changes such as the Dangaard-Oeschger (D-O) cycles and Heinrich events, similar in magnitude and velocity to the ongoing global warming. However, the magnitude of the millennial-scale warming events of the last glacial period was regionally-specific. Precession seems to have imprinted regions below 40°N while obliquity, which controls average annual temperature, probably mediated the impact of D-O warming events above 40°N. A decoupling between high- and low-latitude climate was also observed within last glacial warm (Greenland interstadials) and cold phases (Greenland stadials). The synchronous response of western European vegetation/climate and eastern North Atlantic SSTs to D-O cycles was not a pervasive feature throughout the Quaternary. During periods of ice growth such as MIS 5a/4, MIS 11c/b and MIS 19c/b, repeated millennial-scale cold-air/warm-sea decoupling events occurred on the European margin superimposed to a long-term air-sea decoupling trend. Strong air-sea thermal contrasts promoted the production of water vapor that was then transported northward by the
westerlies and fed ice sheets. This interaction between long-term and shorter time-
scale climatic variability may have amplified insolation decreases and thus explain the
Ice Ages. This hypothesis should be tested by the integration of stochastic processes
in Earth models of intermediate complexity.

Keywords: vegetation, millennial-scale climate variability, Dansgaard-Oeschger cycles, Heinrich events, glaciations, interglacials, Europe, Quaternary

INTRODUCTION

“To look in the oceans for direct evidence of past continental climates seems paradoxical. However, marine sediments contain far better terrestrial paleoclimate records than most continental deposits (Heusser, 1986/1987).” This visionary statement, based on previous seminal works published in the 1960s in a special issue of Marine Geology on Marine Palynology (e.g., Groot and Groot, 1966; Traverse and Ginsburg, 1966; Zagwijn and Veenstra, 1966), anticipated the major contribution that pollen analysis from marine cores would make to Paleoclimatology in subsequent decades. It is now well-established that long and continuous deep-sea sedimentary sequences collected near the continents provide high-quality and chronologically well-constrained pollen records documenting past changes in the vegetation and climate of the adjacent landmasses. Since the pioneering work of Heusser, Hooghiemstra, Rossignol-Strick, Van Campo, and Turon in the 1970 and 1980s in the North Pacific and the Atlantic Oceans and the Mediterranean and Arabian Seas (for references see Supplementary Table S1), an increasing number of deep-sea pollen records has been published from the Iberian and South African margins, the Mediterranean Sea, and the South Pacific and Indian Oceans. They mainly span all or a part of the last 1 million years. These records, which compare without chronological ambiguity marine and terrestrial stratigraphies, are pivotal for documenting potential leads and lags between regional atmospherically-driven vegetation, oceanic conditions and ice dynamics. The relevance of this approach was shown in two seminal papers where the authors demonstrated that the Eemian interglacial, ca. 128,000 years ago (128 ka), was not exactly the terrestrial equivalent of the marine interglacial Marine Isotope Stage (MIS) 5e, and that warmth and forest persisted in south-western Europe during periods of Northern Hemisphere ice growth (Sánchez Goñi et al., 1999; Shackleton et al., 2003). Therefore, pollen-rich marine records are of prime importance for understanding the interactions between the ocean and the atmosphere leading to the orbitally-paced deglaciations and glacial inceptions. More recently, coupled analysis of pollen and marine climatic tracers at a finer resolution has revealed a unique suitability to investigate the timing and amplitude of rapid millennial-scale climate changes in different regions of the world as well (e.g., Sánchez Goñi et al., 2008). This direct sea-land comparison approach also allows testing whether vegetation response to short periods of forcing is in dynamic equilibrium with climate (Webb, 1986) or if there is a lag between climate change and vegetation response (disequilibrium hypothesis, Birks, 1981; Bennett and Willis, 1995), due to species competition or different migration rates.

In the first part of the paper, we focus on the suitability of pollen analysis in deep-sea sediments to trace vegetation and atmospheric conditions through time, highlighting its key role in understanding the mechanisms underlying climate change via the direct comparison with ocean and ice dynamics. We also present a global compilation of pollen records obtained from Quaternary deep-sea cores, and identify the regions where few marine pollen records exist despite their relevance for the study of climate variability. The second part is devoted to the responses of the European vegetation and climate to long-term and, particularly, millennial-scale climate changes. We focus on western Europe because this region is directly affected by North Atlantic millennial-scale iceberg discharges, atmospheric Greenland warming and cooling events (Shackleton et al., 2000), and the close influence of the Northern Hemisphere ice growth and decay, as well as being strongly implicated in the feedback processes bringing the Earth's system to glaciations (Ruddiman and Mcintyre, 1984). The third part of the paper aims to show that the interactions between short-term and long-term climate variability may be the potential missing piece of the Ice Age puzzle (Hodell, 2016), given that the Milankovitch astronomical theory of climate, i.e., changes in the seasonal distribution of the solar energy, alone cannot explain the ice age cycles (Milankovitch, 1941).

UNDERSTANDING PAST CLIMATE
CHANGES

Climate is the sum of meteorological phenomena that characterizes the mean state of the atmosphere (i.e., temperature, precipitation, greenhouse gasses – GHG) over 30 years that depends, in turn, on the dynamics of various components of the climate system (Claussen, 2007). Insolation is the external astronomical (orbital) forcing that determines Earth's regional climates and is defined as the amount of energy per surface unit that the Earth receives from the Sun. Insolation is controlled by the distance between the Earth and the Sun that depends on eccentricity (the shape of the Earth's orbit), obliquity (the tilt of the Earth's axis) and precession (the orientation of the Earth's axis); the latter determines the amplitude of the seasons. These orbital parameters vary over time and trigger major climate changes with cyclicities of 100,000, 40,000, and 21,000 years (e.g., Berger and Loutre, 2004). A change in insolation affects the Earth's five main climatic reservoirs-atmosphere, ocean, land surfaces, cryosphere, and vegetation- and each of them affects, in turn, the Earth's other reservoirs through feedback mechanisms that amplify or reduce the original climate change
(Ruddiman, 2001). Thus the frequency, duration and magnitude of a given climate change is the result of the interactions between orbital external forcing and internal feedback loops involving, for example, ocean currents, particularly the Atlantic Meridional Overturning Circulation (AMOC), ice-sheet and sea-ice dynamics, vegetation, volcanic eruptions, GHG concentrations, and albedo (Ruddiman, 2001).

Marine and terrestrial palaeoclimatic records and ice archives allow the reconstruction of past climate changes in the Earth's different reservoirs. However, understanding the mechanisms controlling the frequency, duration and amplitude of climate changes requires the comparison of these different records on a common timescale, a challenging task due to the fragmentary nature of certain sequences and the uncertainties inherent in different dating methods and age models used. The age models of marine and terrestrial cores are mostly based on radiometric dates ($^{14}$C, U/Th...) and isotopic stratigraphy/orbital tuning, while the ice archives are dated by mean of layer counting or physical models of ice accumulation. The three types of records have uncertainties ranging from a few decades to several centuries and even millennia (e.g., Bazin et al., 2013). Continuous pollen-rich deep-sea sedimentary sequences trace the regional vegetation history and, consequently, the climate of the landmasses nearby and represent a unique way to circumvent the aforementioned issues when comparing different records (vegetation/atmosphere, marine, and ice). Marine climatic indicators from these sequences allow the quantitative reconstruction of sea surface temperature (SST) and salinity (e.g., planktonic foraminifera, dinoflagellate cysts, calcareous nannofossil and diatom assemblages, alkenones), deep ocean conditions (e.g., Mg/Ca, benthic foraminifer assemblages, carbon isotopic ratio -δ$^{13}$C- of benthic foraminifera, Pa/Th), iceberg dynamics (Ice Rafted Debris or IRD, the coarse sediments transported by icebergs) and the ice volume stored in the ice caps [oxygen isotopic ratio -δ$^{18}$O- of benthic foraminifera; although a bias exists between this isotopic ratio and global ice volume (Skinner and Shackleton, 2006)] (Figure 1). This approach is usually named "direct land-sea correlation or comparison," although the insights extend not only to terrestrial and marine conditions but also to the cryosphere. Other indicators found in marine sediments give complementary information on terrestrial climate such as charcoal particles providing relevant information about climatically-driven fire regime prior the beginning of human impact on the environment (Daniau et al., 2007, 2013), dust concentrations (e.g., deMenocal, 2004; Jullien et al., 2007) and terrestrial organic compounds or biomarkers (δD, GDGT, δ$^{13}$C from leaf waxes) that give useful information on precipitation and, in some cases, on temperatures (e.g., Weijers et al., 2007).

This direct correlation allows the identification of possible synchronicity or time-lags between the response of these reservoirs to a given climate change and, therefore, provides insights into the feedback processes underlying this climate change. Finally, the comparison between directly correlated palaeoclimatic data and model simulations allows the assessment of the reliability of the mechanisms implicated in the different models to reproduce the observed climate change (e.g., Sánchez

Goñi et al., 2012). However, before identifying regional land-sea linkages associated with past climate changes, it is crucial to understand the present-day regional pollen signals.

## Pollen in Deep-Sea Sediments

A substantial quantity of pollen grains produced, released and dispersed from terrestrial higher plants (as well as spores of lower plants) reaches the ocean mainly by fluvial and atmospheric transport processes, with oceanic currents playing a negligible role (Koreneva, 1966; Stanley, 1966). The relative contribution of one or other transport vector is regionally-dependent (Dupont, 2011). In the Arctic, for instance, pollen is also transported by sea ice scouring and sediment transport in addition to northeasterly winds (Mudie and McCarthy, 2006). Marine sediments located further than 300 km offshore are weakly influenced by rivers and the dominant pollen transport is aeolian (e.g., Hooghiemstra et al., 1988). In contrast, the pollen preserved in marine sediments located under the influence of river plumes (nepheloid layers) is mainly of fluvial origin (e.g., Heusser and Balsam, 1977; Cambon et al., 1997; Mudie and McCarthy, 2006). Once it arrives in ocean surface waters, pollen is ingested by planktonic organisms and later integrated into their fecal pellets or agglomerated with clays (Mudie and McCarthy, 2006; Beaudouin et al., 2007). Due to these processes, pollen buoyancy decreases and is little influenced by ocean currents. It thus becomes an integral part of marine snow and crosses the water column with a relatively high speed (estimated at ~100 m/day in the Atlantic water column) before being deposited at the bottom of the oceans (Hooghiemstra et al., 1992). Furthermore, sediments under the influence of upwelling are particularly rich in pollen for two reasons: the intensification of the downward particle transport through the water column (Ratmeyer et al., 1999) and the better preservation in almost anoxic sediments (Dupont, 2011).

Pollen studies on modern deep ocean surface sediments from the Iberian margin (Turon, 1984a; Naughton et al., 2007), the Bay of Biscay (Turon, 1984b), the Mediterranean (Koreneva, 1971; Beaudouin et al., 2007), the African margin (Hooghmiestra et al., 1986; Dupont and Wyputta, 2003), and the western and eastern margins of North America and the Gulf of Mexico (Heusser and Balsam, 1977; Heusser, 1985; Heusser and Van de Geer, 1994) showed that pollen assemblages from the sediment surface of the ocean floor reflect an integrated image of the regional vegetation of the adjacent continent and, consequently, the climatic parameters under which this vegetation has developed. Palynological richness, i.e., the number of pollen types per sample standardized to a constant pollen sum, is similar in fossil samples from deep-sea cores, modern terrestrial and deep-sea surface samples, and fossil assemblages from lakes and estuaries, oscillating between 15 and 25 morphotypes with respect to a pollen sum of at least 100 terrestrial pollen grains (e.g., Sánchez Goñi and Hannon, 1999; Naughton et al., 2007). An ongoing study on the pollen representation of modern vegetation in moss-polsters along the Tajo basin shows that palynological richness ranges from 13 to 41 with respect to a pollen sum of 500 terrestrial pollen grains. In particular, Naughton et al. (2007) demonstrated by comparison of modern pollen spectra from the deep ocean, estuaries and continent that the pollen signal of the
<figure><img src="image_2.png" /><figcaption>**FIGURE 1** | Illustration of a sediment core from a continental margin, which captures multiple terrestrial and oceanic tracers, including: transport and deposition of pollen grains, coarse sediments derived from the discharge and melting of icebergs [ice-rafted detritus (IRD)], planktonic (near surface) and benthic (bottom-dwelling) foraminifer, and dinoflagellates. Adapted from Sánchez Goñi (2016).</figcaption></figure>

Iberian margin is similar to that found in estuarine sediments from western Iberia that reliably represents the broad regional vegetation of the related hydrographic basins (Figure 2). On the one hand, pollen assemblages from the northern part of the margin reflect the composition of temperate forest that develops in northwestern Iberia dominated by deciduous oaks. On the other hand, pollen samples from the southern part of the margin capture the composition of Mediterranean forest, i.e., warm-temperate forest dominated by sclerophyllous trees and shrubs that characterizes southern Iberia (Figure 2). In this region, pollen is mainly transported by the rivers because the dominant winds come from the north-west and the hydrographic Iberian basins are large and thus favor the transport of sediment load, including pollen (Dupont and Wyputta, 2003; Naughton et al., 2007). From a botanist's point of view, this type of study has therefore a limitation for reconstruction of ancient local plant communities (Groot and Groot, 1966). However, the similarity of western European terrestrial pollen sequences and eastern North Atlantic deep-sea pollen records (e.g., Naughton et al., 2007 for the last 18,000 years in north-west Iberia; Fletcher and Sánchez Goñi, 2008 for the last 50,000 years in south-east Iberia; Sánchez Goñi et al., 2012 for the Last Interglacial in western

France) demonstrate once more that marine pollen records provide a reliable image of the vegetation history of the adjacent landmasses (Groot and Groot, 1966). Moreover, the regional vegetation is directly linked to climate conditions as the present-day distribution of the major biomes is governed by climatic parameters (Bailey, 1998). For western Europe Gouveia et al. (2008) have recently shown the rapidity with which vegetation in this region responds to North Atlantic atmospheric processes, i.e., the westerlies.

The comparison between terrestrial and marine modern pollen samples to better interpret marine pollen records has been further improved and applied to south-western Africa (Urrego et al., 2015). In particular, modern pollen spectra were used to assess the distribution of Poaceae pollen abundance and other pollen taxa with potential indicator value for large biomes in southern Africa, and therefore climatic zones in this region. For instance, pollen percentages of Poaceae up to 70% reflect the Nama-karoo and its transition with the fined-leaved savanna. These results can substantially change the interpretation of the marine pollen fossil record collected off northwestern South Africa. Contrary to the studies that interpret the increase of Poaceae pollen percentages as the result of
<figure><img src="image_3.png" /><figcaption>**FIGURE 2** | Map of the Iberian Peninsula with the main types of potential vegetation (Braun-Blanquet and Walter in Blanco Castro et al., 1997), present-day pollen signature in surface sedimentary samples from the Ría de Vigo and Barreiro (Tajo estuary), and samples from the top samples of cores MD99-2331 and MD95-2042 covering the last centuries and located in the Atlantic and Mediterranean region, respectively.</figcaption></figure>

humidity-driven savanna expansions in southwestern Africa, the new pollen calibration study offers an alternative interpretation indicating that fine-leaved savanna developed in this region due to the aridity increase during the warm and humid Northern Hemisphere periods of the last interglacial (Urrego et al., 2015). This approach should be applied to other regions to improve the interpretation of marine pollen assemblages in terms of vegetation cover and composition. However, the quantification of the vegetation abundance from marine pollen assemblages to infer, among others, changes in albedo and their influence on climate change, remains a challenge. The REVEALS model was conceived to obtain quantitative reconstructions of regional vegetation cover around large lakes from pollen data, which resembles the spatial scale of marine pollen samples (Sugita, 2007). Nevertheless, this model only considers atmospheric pollen dispersal and deposition, thus disregarding the importance of pollen input from inlet streams and surface run-off, which are often the main vectors delivering pollen to the marine environment (Sugita, 2007).

The concentration (usually between 1,000 and 50,000 grains cm⁻³) and taxonomical diversity of marine pollen assemblages are in general high enough to allow the quantitative estimation of annual and, more importantly, seasonal temperatures and precipitation (e.g., Sánchez Goñi et al., 2002, 2005, 2016a;

Desprat et al., 2007; Combourieu Nebout et al., 2009) (Figure 3). Knowledge about past climate seasonality is crucial to understand natural climate variability (Carré and Cheddadi, 2017) and the sensitivity of vegetation to changes in this parameter makes pollen data particularly well-suited to track shifts in seasonality through time. Quantitative pollen-based climatic estimates are based on different techniques. The most frequently used is the Modern Analog Technique (MAT) based on a comparison of past assemblages to modern pollen assemblages through the calculation of the shortest weighted distance (Guiot, 1990). This method requires high-quality, taxonomically consistent modern pollen and climate datasets. The modern database includes different continental pollen spectra (moss polsters, soil surface samples, and lacustrine sediments) from Eurasia and northern Africa (Peyron et al., 1998; Bordon et al., 2009). This continental pollen database is, however, slightly modified to estimate climatic parameters from marine pollen samples. Since *Pinus* pollen is overrepresented in marine sediments (Heusser and Balsam, 1977; Hooghiemstra et al., 1988), this pollen type is excluded from the calculation of pollen percentages for marine pollen spectra as well as from the continental pollen database. As potential errors would be similar for adjacent samples of the same record, we consider that this approach is still valid for the quantification of climate
<figure><img src="image_4.png" /></figure>

**FIGURE 3** | Maps of western Europe with pollen-based quantitative climate estimates, winter temperature and annual precipitation, and sea surface temperature (SST) on its margin during Heinrich stadial (HS) 4 and Greenland interstadial (GI) 8. (A) Air and sea surface winter temperatures during HS 4, (B) air and sea surface winter temperature during GI 8, (C) annual precipitation during HS 4, and (D) annual precipitation during GI 8. These maps interpolate the temperatures and precipitation estimated from four cores strategically located in the Atlantic and Mediterranean regions (Sánchez Goñi et al., 2002 and unpublished data). Black circles indicate the location of cores MD95-2043, MD95-2042, MD99-2331, and MD04-2845 (Naughton et al., 2007; Fletcher and Sánchez Goñi, 2008; Sánchez Goñi et al., 2008).

change anomalies. The application of the MAT to estimate quantitatively climatic parameters for marine pollen records should be further improved by integrating marine modern pollen assemblages in the database, and by application of other reconstruction approaches that rely less strongly on the availability of suitable modern analogs, such as indicator approaches (Kühl et al., 2002) and biomisation (Prentice et al., 1996).

## A Global Compilation of Deep-Sea Pollen Records

During the last 60 years, several oceanographic cruises organized in the framework of the international programs IMAGES (International Marine Global change Study) and ODP/IODP (International Ocean Drilling Program) have retrieved a large number of pollen-rich marine cores along different oceanic margins. Most of these cores come from oceanic rises and abyssal plains that are the most favorable areas for marine palynological work as terrigenous sediments are predominant

(Groot and Groot, 1966). We have performed a literature survey to compile pollen records from deep-sea cores spanning the last 1 million years, finding 129 sites all over the world (Supplementary Table S1 and Figure 4). This compilation shows that there are 74 high-resolution pollen sequences (better than 1000 years between adjacent samples) and highlights the paucity of pollen records from the Indian (e.g., Rossignol-Strick, 1983) and the South Pacific (e.g., Montade et al., 2013; Seillès et al., 2016) Oceans. Likewise, almost half of the records (n = 57) do not cover periods older than the Last Glacial Maximum, and few sequences (n = 19) record several orbital climatic cycles. Therefore, the regional expression of long-term and rapid, millennial-scale, global climate changes is far from being well-documented, particularly in the Indian Ocean where only two out of the six available records cover the last climatic cycle (last 150,000 years) with a coarse time resolution (worse than 1000 years between consecutive samples). To fill this gap, IODP expedition 353 "Indian Monsoon Rainfall" collected a sedimentary sequence
<figure><img src="image_5.png" /></figure>

**FIGURE 4** | Global map with the main types of potential vegetation (Levavasseur et al., 2012) and the marine pollen records covering different intervals of the last 2.58 Myrs. We have restricted the selection to pollen records from sections cored in the present-day marine realm. Black circles indicate the pollen records with a temporal resolution better than 1000 years (<1000 years) between adjacent samples. White circles indicate the pollen records with a temporal resolution worse than 1000 years (>1000 years) between adjacent samples. Gray circles indicate records where no information about the resolution is given. Left bottom corner: histogram indicating the number of marine pollen records by region and the number of high resolution records among the 129 records compiled in this work (see **Supplementary Table S1** for details and references). Warm-temperate forest biome includes the Mediterranean forest; Temperate forest biome includes the Atlantic forest; Grassland and dry shrubland include the Mediterranean plants, i.e., sclerophyllous trees and shrubs such as *Quercus ilex* type (holm oak, kermes oak), *Q. suber* type (cork oak), *Olea* (olive tree), *Phillyrea*, *Pistacia*, and *Cistus* (rockroses); Desert and Tundra-dry tundra are the present-day closest vegetation biomes to the past semi-desert and steppe environments inferred from the pollen analysis of last glacial sedimentary sequences, respectively.

in 2015 on the eastern Indian margin, site U1446, which will enable to trace for the first time past vegetation and Indian summer monsoon variability over the last 1 million years.

# EUROPEAN VEGETATION AND CLIMATE RESPONSE TO LONG-TERM AND RAPID CLIMATE VARIABILITY

## General Setting

Since 1995, the Marion Dufresne and Joides Resolution oceanographic ships have retrieved several high-resolution pollen-rich deep-sea sedimentary sequences from the eastern North Atlantic and the Mediterranean Sea within the framework of the IMAGES and IODP programs. These sequences are strategically located in the present-day Mediterranean and Atlantic regions, characterized by a markedly seasonal climate, with warm and dry summers and wet and mild winters, and a year-round humid climate with warm summers and cool winters, respectively. In winter, both regions are under the influence of the westerlies, but in summer the Azores high affects southern Iberia inducing a pronounced summer drought while the westerlies still play a prominent role from northern Iberia northward. Additionally, the region located below 40°N is affected by

precession, and the subtropical ocean water currents (Ruddiman and McIntyre, 1984).

At present, the seasonal distribution of precipitation, with a relatively abundant amount of rainfall from autumn to spring and a pronounced summer drought, determines the dominance of Mediterranean forest (Polunin and Walters, 1985; Gouveia et al., 2008). This forest comprises evergreen sclerophyllous woodland (evergreen *Quercus*, *Olea*, *Pistacia*, *Phillyrea*, *Cistus*) at low altitudes, and deciduous forest (mostly deciduous *Quercus*) at higher elevation. On the contrary, the occurrence of Atlantic forest, dominated by deciduous *Quercus* (oak), is today controlled by winter length and harshness, given that annual rainfall is high and summer drought absent or negligible (Polunin and Walters, 1985). Therefore, deciduous *Quercus* is included in both the 'Atlantic forest' and the 'Mediterranean forest' assemblages.

Pollen analysis of these European margin sequences has allowed the reconstruction of the responses of western European vegetation and climate to long- and short-term global climate changes of the Ice Ages.

## European Vegetation and Climate Response to Past Global Long-term Climate Changes

Geological archives show that the Earth's climate has experienced large changes during the last 2.58 million years (Shackleton
<figure><img src="image_6.png" /></figure>

**FIGURE 5** | Pollen and isotopic (δ¹⁸O₄) records from the southwestern Iberian margin (Sánchez Goñi et al., 1999, 2016b; Roucoux et al., 2006; Chabaud et al., 2014; Hodell et al., 2015; Oliveira et al., 2016; and unpublished data). The Mediterranean forest is mainly composed of deciduous oaks and Mediterranean plants, i.e., sclerophyllous trees and shrubs such as *Quercus ilex* type (holm oak, kermes oak), *Q. suber* type (cork oak), *Olea* (olive tree), *Phillyrea*, *Pistacia*, and *Cistus* (rockroses). Pollen data are compared to CO₂ and CH₄ records from ice cores (Siegenthaler et al., 2005; Spahni et al., 2005; Louergue et al., 2008; Luthi et al., 2008) and astronomical forcing (Berger, 1978). Adapted from Desprat et al. (2017).

and Opdyke, 1973; Hays et al., 1976), with alternating warm (interglacial) and cold (glacial) periods. During the glacial periods ice occupied large regions of the Northern Hemisphere resulting in low sea level stands while during interglacials the ice sheets only covered Greenland and sea level was similar or higher than at present. These changes in ice volume are identified in marine sedimentary sequences by measuring the oxygen isotope composition (δ¹⁸O) of benthic foraminifer carbonate shells, providing the basis of the marine isotopic stratigraphy. Depleted values define the interglacial periods *sensu lato* (s.l.) or odd numbered MIS, while heavy values mark the glacial periods corresponding to even-numbered MIS. Within interglacial periods, an alternation of ice retreats, called “e,” “c” and “a,” and ice advances, “d” and “b,” is observed, also resulting from orbital forcing. The interglacial period *sensu stricto* (s.s.) corresponds to the substage with minimal ice volume, often located just after deglaciation (Railback et al., 2015).

Two long and continuous pollen-rich sedimentary sequences are available from the southwestern and northwestern Iberian margin covering the last 1.5 million years and 425,000 years, respectively. These Iberian margin sequences are each composed of three cores: the southern sequence is a composite of cores MD95-2042, MD01-2443, and IODP U1385 located at around 37°N (**Figure 5**; Sánchez Goñi et al., 1999, 2016b; Roucoux et al., 2006; Chabaud et al., 2014; Oliveira et al., 2016), and the northern one is a composite of cores MD99-2331, MD01-2447, and MD03-2697 collected at 42°N (Desprat et al., 2005, 2006, 2007, 2009; Sánchez Goñi et al., 2005, 2008; Naughton et al., 2007, 2009, 2016). The pollen analysis of the 425,000-year northern sequence is completed (Desprat et al., 2017) while that of the southern site is in progress (**Figure 5**). Further north,

the sequence MD04-2845 collected on the European margin, off western France, spans the last 140,000 years (Sánchez Goñi et al., 2008, 2012).

Both Iberian margin pollen sequences show an alternation of glacial periods dominated by Mediterranean steppe in southern Iberia (that we call semi-desert although it rarely included true semi-desert taxa such as *Lygeum*) and steppe/heathlands to the north, and interglacial periods when Mediterranean and Atlantic forests expanded over southern and northern Iberia, respectively. Tree pollen percentages are a good estimate for the magnitude of the interglacial in both Iberian regions. The magnitude of past Mediterranean forest expansions is interpreted as an indicator of the amount of winter precipitation in southern Iberia, given that relatively warm temperatures and summer drought occurred. The magnitude of the Atlantic forest expansion during different interglacials will mainly reflect the amplitude of the warming above 40°N, as forest development in the Atlantic domain requires a climate with a growth period of 4-6 months, and a mild winter period of 3-4 months (Polunin and Walters, 1985). Two or even three major phases of forest expansion occurred during all these interglacials, related with ice volume minima during sub-stages “a,” “c,” and “e” within each MIS (**Figure 5**). However, forest development in general reached its maximum during the earliest phase of each stage, concomitant with the highest sea level, indicating that warmest and wettest conditions occurred then, and subsequently identifying this phase with the terrestrial interglacial *sensu stricto*.

Pollen records from the Iberian margin show that the magnitude of forest development substantially differed from one interglacial to another in the south (MIS 19 to MIS 1, Mediterranean forest pollen = 40 to 80%, **Figure 5**), while in
the north the magnitude of the forest expansion was similar (Atlantic forest pollen = 80%; Desprat et al., 2017). These data suggest that temperatures were similarly warm during the last five interglacials in northwestern Iberia, whereas the amount of winter precipitation in the south was quite variable. The Holocene (MIS 1, ~10 ka), the last interglacial (MIS 5e, ~128 ka) and MIS 17c (~690 ka) would have been the wettest followed by MIS 9e (~335 ka), MIS 11c (~415 ka) and MIS 19c (~784 ka), and MIS 7e (~240 ka) being the driest. Excluding the Lake Orhid pollen record (Sadori et al., 2016), these results coincide with what is observed in the European sequences of Tenaghi Philippon and Praclaux located above 40°N (Reille et al., 2000; de Beaulieu et al., 2001; Tzedakis et al., 2006), which show similar expansions of the temperate forest during the interglacials of the last 400,000 years. So far no long terrestrial pollen sequence exists below 40°N for comparison with the southwestern Iberian margin pollen record.

The comparison of the Iberian observations with those at higher latitudes illustrates the regional variability of the magnitude and impacts of climate change. Due to the combined forcing of insolation and GHG, climatic models indicate that MIS 5e, MIS 9e, and MIS 11c were the warmest at the highest latitudes, associated with the strongest melting of ice sheets (Yin and Berger, 2012, 2015). This scenario is in marked contrast with data from southwestern Europe showing that these interglacials were not particularly warm. Our pollen data suggest that substantial differences in precipitation occurred across the different interglacials. Nevertheless, it is difficult to propose a relationship between warming at high latitudes due to GHG increase and dryness in the Mediterranean region, contrary to what climate models project for the end of this century (Hoerling et al., 2012). Identifying forcing factors responsible for the interglacial diversity in terms of precipitation in the Mediterranean region is a subject of debate. A good candidate is the orbital parameter of precession that influences regions below 40°N (Ruddiman and McIntyre, 1984). Precession controls the amplitude of the seasonality and, therefore, the amount of winter precipitation in subtropical regions as recently suggested by models (Yin and Berger, 2012) and observations (e.g., Sánchez Goñi et al., 2008; Oliveira et al., 2016, 2017).

## Rapid Climate Changes: A Focus on the Last Glacial Period

More recently, geological archives have also shown that millennial-scale changes punctuated this long-term climatic variability at the global scale as these changes are registered in North Atlantic SSTs (Voelker and Workshop Participant, 2002; Martrat et al., 2007) as well as in air temperatures over Greenland (North Greenland Ice-Core Project Members [NGRIP], 2004; Barker et al., 2011) and Antarctica (Augustin et al., 2004). These changes occur independently of the boundary climate state, i.e., in glacial and interglacial periods, and are generally associated with changes in ice volume (Siddall et al., 2003), GHG concentration (Knutti et al., 2004; Spahni et al., 2005; Loulergue et al., 2008), and AMOC

(Oppo et al., 2006; Lynch-Stieglitz, 2017). The magnitude and frequency of these short-term climatic oscillations were larger when ice caps reached a critical mass (McManus et al., 1999) even if they were not the trigger (Barker et al., 2015).

So far, the best studied millennial-scale climate variability is that observed during the last glacial period (MIS 4, MIS 3 and MIS 2, i.e., ~73-14.7 ka) characterized by a series of warming and cooling events called Dansgaard-Oeschger (D-O) cycles that were first identified in the atmosphere of Greenland and usually lasted 500-2,000 years (Dansgaard et al., 1984). These cycles recorded in the δ¹⁸O of the ice cores were characterized by a large (7 to 16°C) and rapid (within a few decades) warming event followed by a progressive decrease in temperature and a final abrupt cooling (Wolff et al., 2010). The warming and progressive cooling phase is termed Greenland Interstadial (GI), and the final cooling leading to the cold phase termed Greenland Stadial (GS). The GI phases lasted between 100 and 2,600 years (Wolff et al., 2010). Also during the last glacial period, repeated iceberg discharges cooled the surface of the North Atlantic (Heinrich, 1988; Bond et al., 1993). On the one hand, massive iceberg discharges from the Laurentide ice sheet, the so-called Heinrich events (HE), occurred with a cyclicity of 7,000-10,000 years while, on the other hand, weaker discharges coincided with iceberg fragmentation from the British-Icelandic-Scandinavian (BIS) ice cap (Elliot et al., 2001). We define a HE as the time period synchronous with the deposition of the coarse sediment Heinrich layer in a given region following the iceberg discharge while the Heinrich Stadial (HS) is the cold interval associated with the HE (Sánchez Goñi and Harrison, 2010). Some GSs encompass the HSs while the others are associated with the BIS minor iceberg discharges. These cold intervals are related to decreases in the AMOC (McManus et al., 2004; Lynch-Stieglitz, 2017) but the cause of the ice sheet collapse remains a subject of debate (e.g., Alvarez-Solas and Ramstein, 2011; Barker et al., 2015). The duration of the iceberg discharge has been simulated to be between 50 and 200 years (Roche et al., 2004) and data show that the HSs lasted longer, up to 3000 years (Sánchez Goñi and Harrison, 2010). HE and HS time intervals coincide in the Ruddiman belt, the preferential zone of IRD deposition that extends roughly between 45 and 50°N, while this is not the case outside (Naughton et al., 2009). Since the discovery of these millennial-scale climate changes, D-O cycles and HEs, the paleoclimatic community has dedicated major efforts to investigate the regional expression of this variability, the oceanic and atmospheric mechanisms involved in its transmission, and the interaction with long-term climate forcing such as the orbital parameters and ice volume. The following subsection summarizes the great contribution that palynological research conducted on marine sediments has done so far to increase current understanding of these processes in western Europe.

### D-O Cycles

Until the end of the 1990s, investigations into European pollen records reported no changes in the dominant steppe
vegetation through the last glacial period (e.g., Follieri et al., 1988; Pons and Reille, 1988; Reille and de Beaulieu, 1990). Allen et al. (1999) published a pollen record from southern Italy showing for the first time millennial-scale changes between forest expansion and contraction that corresponded to the rapid climate variability detected in Greenland and the North Atlantic Ocean. However, the chronological uncertainties between ice, terrestrial and marine sequences precluded the demonstration of whether forest expansions corresponded to North Atlantic and Greenland warming or, inversely, to cooling events. Marine pollen records have unequivocally shown that cooling events in the North Atlantic and high latitudes corresponded in western Europe with herbaceous community expansion and warming events with enhanced forest development (Roucoux, 2000; Sánchez Goñi et al., 2000, 2002, 2008). The aforementioned marine pollen records from the European margin and an additional record, MD95-2043, from the Alboran Sea (western Mediterranean, Figure 4) (Sánchez Goñi et al., 2002; Fletcher and Sánchez Goñi, 2008) also show that vegetation quickly, within 100 years, responded to the D-O cycles and HEs, and that there was a dynamical equilibrium between vegetation response and climate change for short periods of forcing. Cold episodes in North Atlantic surface temperatures related to GS, not only those associated with HEs, were synchronous with forest reductions in Iberia and western France, while forest expansions correlated with increases in SST (Figure 6) (Sánchez Goñi et al., 2008).

Notably, below 40°N, Mediterranean forest reached its maximum development at the onset of the Eemian (MIS 5), the D-O 24, D-O 21, D-O 17-16, D-O 8-7, and D-O 1 warming events, and the beginning of the Holocene (MIS 1), indicating the occurrence of enhanced hot/dry summers and wet/cool winters (Figure 7). The Atlantic sites, above 40°N, showed a contrasting pattern: during GI 12 and GI 14 the Atlantic forest experienced a strong development while the impact of D-O 17-16 and D-O 8-7 warming was rather limited. The comparison of the Mediterranean and Atlantic palaeoclimatic records reveal that there was a spatial variability in the amplitude of the forest expansions for any given D-O warming of the last glacial period (Sánchez Goñi et al., 2008) (Figure 7). The comparison of these changes in vegetation with the evolution of orbital parameters shows that the maxima in the Mediterranean forest were always synchronous with low precession values. This observation suggests that precession minima strengthened Mediterranean climate through promoting marked seasonality. The floristic composition of the Mediterranean forest corroborates that Mediterranean climate reached its best expression during precession minima (Sánchez Goñi et al., 2008). For example, pollen data from the Alboran Sea core show that Mediterranean sclerophylls such as evergreen Quercus, Olea, Cistus, Phillyrea, Coriaria myrtifolia, and Pistacia were particularly abundant during GI 8, under enhanced seasonality associated to a minimum in precession (Fletcher and Sánchez Goñi, 2008). On the contrary, the same sequence records that the abundance of the less drought-tolerant Ericaceae (heather) was higher during GI 12,

when precession reached a maximum and seasonality was reduced (Fletcher and Sánchez Goñi, 2008). In contrast, the maximum expansion of the Atlantic forest coincides with maxima in obliquity. For northern latitudes above 40°N obliquity seems to play a major role in modulating the increase of annual temperatures (Figure 7). The amplitude of temperature changes in Greenland was strong at D-O 17-16, D-O 12 and D-O 8, and weak at D-O 14 (Sánchez Goñi et al., 2008), thus also differing from the climatic patterns at lower latitudes. This means that the strongest warming events in Greenland were not necessarily particularly warm in other regions of the world. The Iberian Peninsula and wider Mediterranean region are particularly interesting from a palaeoclimatic point of view because they appear to straddle the transition zone from obliquity- to precession-influence on the climatic signal of millennial-scale warming events. Preliminary comparison between European margin and terrestrial pollen records across Europe roughly confirms the contrasting latitudinal response over Europe to the D-O warming events despite the independent and sometimes uncertain chronologies of individual records (Fletcher et al., 2010). A more in-depth comparison of these pollen records is in progress based on the harmonized chronology recently developed in the framework of the INQUA International Focus Group (IFG) ACER project (Sánchez Goñi et al., 2017).

A zoom into D-O 7 and D-O 8 cycles, 41,000–34,000 ka, shows that semi-desert plants (mainly *Artemisia*, *Chenopodiaceae* and *Ephedra*) characterized the vegetation of southern Iberia during GSs, with open Mediterranean woodland establishing during GIs (Sánchez Goñi et al., 2008; Figure 8). Simultaneously, Ericaceae and Poaceae dominated the vegetation of north-western Iberia during GSs (Figure 8), being partially replaced by Atlantic woodlands with *Pinus* and deciduous *Quercus* during warmer periods (GIs). The composition of the pinewoods of south-western/central Iberia during the last glacial also varied in response to abrupt climate changes according to the quantitative study of *Pinus* pollen grains by Desprat et al. (2015). Thus, although *P. nigra* was the dominant pine species throughout the last glacial period, *P. sylvestris* was more abundant during the colder GSs and HSs and Mediterranean pines (*P. pinaster*, *P. pinea*, *P. halepensis*) during the warmer GIs and early-to mid-Holocene. Further north, *Artemisia*, *Cyperaceae* and *Calluna* dominated in western France during the cold periods, whereas *Betula*, deciduous *Quercus* and conifers (*Pinus*, *Abies*, *Picea*) expanded slightly in the landscape during warming episodes (Sánchez Goñi et al., 2008; Figure 8). These sequences illustrate the spatial floristic diversity of western Europe in response to the D-O cycles. Interestingly, the latitudinal boundary between the Atlantic and the Mediterranean vegetation during the last glacial period seems to have been similar to that at present-day, i.e., below 40°N. Finally, it is worth highlighting that temperate forest expanded when summer SST in the North Atlantic crossed the threshold of 12°C (Sánchez Goñi et al., 2008) (Figure 8), which is the same threshold value that explains the distribution of the temperate forest in both sides of the North Atlantic Ocean at present (Van Campo, 1984).
<figure><img src="image_7.png" /></figure>

**FIGURE 6** | Vegetation response to the millennial-scale variability of the last glacial period in southwestern Iberia. Terrestrial and marine palaeoclimatic record from the twin cores MD95-2042 and SU81-18 (SW Iberian margin) for the last 138 ka compared with Greenland temperature variations. (A) Greenland temperature curve (100 years resolution) derived from air and water isotopic measurements. Numbers 1–25 refer to the D–O warming events, (B) pollen percentage curve of the Mediterranean forest (deciduous and evergreen *Quercus*, *Olea*, *Phillyrea*, *Pistacia*, and *Cistus*), (C) percentage curve of the polar planktonic foraminifer *Neogloboquadrina pachyderma* (s), (D) concentration curve of the ice-rafted debris (IRD, in grains per gram of dry sediment) reflecting the main episodes of iceberg melting in the Iberian margin. Thin line indicates exaggeration x10, (E) curve of the alkenone-derived (UK'37) sea-surface temperatures (SST), (F) planktic foraminifer δ¹⁸O curve on Globigerina bulloides, (G) benthic foraminifer δ¹⁸O curve indicating the marine isotopic stratigraphy. Blue intervals indicate the Heinrich stadials associated with ice-rafting events of the last climatic cycle: HS1–HS6 and C19–C24 (Bond et al., 1993; Chapman and Shackleton, 1999). N.B.: for the last 25,000 years, the IRD, *N. pachyderma* (s) and pollen data are from the twin core SU81-18 (Turon et al., 2003) that are similar to that obtained from core MD95-2942 (Chabaud et al., 2014). Adapted from Sánchez Goñi et al. (2008).

## A Close-up on HS and GI Phases

Marine pollen records unequivocally show that HSs were associated with the expansion of herbaceous vegetation reflecting cold and dry conditions in western Europe. However, a closer observation of HS 4 (40.2–38.3 ka), HS 2 (26.5–24.3 ka), and HS 1 (18–15.6 ka) in core MD99-2331 retrieved off northwestern Iberia suggests that two major phases occurred

within these stadials (Figure 9) (Naughton et al., 2009). The first phase was characterized by almost no IRD in the north-western Iberian margin but SSTs were low and large reductions of *Pinus* and Atlantic forest pollen indicate very cold atmospheric conditions as well. The first episode was the wettest one according to the maximum percentages of *Calluna*, typical at present of central European mires
<figure><img src="image_8.png" /><figcaption>**FIGURE 7** | Evolution of Greenland temperature, western European temperate forest development (37–45°N) and SSTs (37°N) over the last climatic cycle compared with iceberg discharges (37°N), precession index, obliquity and ice volume variations. (A) Greenland temperature reconstruction (100 years resolution). Numbers 1 to 24 indicates D-O warming events, (B) Atlantic forest pollen percentages in western France from deep-sea core MD04-2845 (136–35 ka), (C) Atlantic forest pollen percentages in northwestern Iberia from deep-sea cores MD99-2331 (131–18 ka) and MD03-2697 (14–0 ka), (D) Mediterranean forest pollen percentages in southwestern Iberia from deep-sea cores MD95-2042 (27–134 ka) and SU8118 (0–27 ka), (E) precession index and obliquity, (F) benthic δ18O curve from deep-sea core MD95-2042. Gray rectangles indicate the North Atlantic ice-rafting events of the last climatic cycle: HS1 to HS6 and C19 to C24 (Bond and Lotti, 1995; Chapman and Shackleton, 1999). Orange rectangles indicate the maximum expansion of the Mediterranean forest. Adapted from Sánchez Goñi et al. (2008).</figcaption></figure>

(Polunin and Walters, 1985), and high pollen concentration (not shown). The second phase of these stadials was drier and warmer as shown by the increase of semi-desert plants, and *Pinus* and temperate forest pollen percentages, respectively.

A synthesis of all available North Atlantic sites recording the last four HS shows complex climatic and IRD deposition patterns that are summarized in **Figure 10** (Naughton et al., 2009). In the first phase, the wettest and coldest conditions in the Iberian Peninsula were coetaneous with low IRD deposition and low SST in the Iberian margin. A contrasting pattern is detected in the Ruddiman belt (45–50°N) with high IRD deposition and cold SST while in Florida a dry climate was associated with cool SST off this peninsula. During the second phase, dryness in Iberia

and high IRD deposition in its margin are observed at the same time as wetness in Florida and moderate IRD deposition in the Ruddiman belt. A simple oceanographic mechanism related to changes in the strength of the AMOC alone cannot explain this complex scenario but rather an atmospheric mechanism must be involved. During the first phase, strong SST cooling of the North Atlantic as far as the southern Azores and the Gulf of Cádiz was associated with relatively high precipitation in south-western Europe. At the same time, δ¹³C measurements on benthic foraminifer and Pa/Th ratios indicate a gradual decrease of deep sea ventilation and a slowdown of the AMOC, respectively (Naughton et al., 2009). The southward displacement of the thermal front produced a strong mid-latitude meridional temperature gradient that led to the southward displacement
<figure><img src="image_9.png" /></figure>

MD04-2845
off western France

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

HS 4

<figure><img src="image_12.png" /></figure>

MD99-2331
off north-western Iberia

<figure><img src="image_13.png" /></figure>

MD95-2042
off south-western Iberia

<figure><img src="image_14.png" /></figure>

HS 4

<figure><img src="image_15.png" /></figure>

**FIGURE 8** | A "zoom" on the interval 41–34 ka, from the three western European margin cores. The location of these cores is shown in **Figure 3**. GI (Greenland Interstadial), HS (Heinrich Stadial), GS (Greenland Stadial). In contrast with the Atlantic region, *Pinus* pollen percentages in core MD95-2042 (not shown) are similar through all the considered interval, during GI and GS, and oscillate between 60 and 90%. The pictures named HS 4 and GI 8 suggest analog present-day landscapes: semi-desert (*Ephedra distachya*-type, *E. fragilis*-type, Chenopodiaceae that now is included in the Amaranthaceae pollen morphotype, *Artemisia*) and open forest in the Mediterranean region and steppe (*Poaceae*, *Cyperaceae*, *Asteraceae*, *Ericaceae* and semi-desert plants) and *Pinus-Betula* forest in central Europe, respectively. SST: sea surface temperatures estimate from planktonic foraminifer assemblages using the Modern Analog Technique (MAT) (Guiot and De Vernal, 2007). Blue rectangles indicate the cold intervals, HS 4 and GS 8. Adapted from Sánchez Goñi (2014).
<figure><img src="image_16.png" /></figure>

**FIGURE 9** | Direct land-sea correlation of MD99-2331 record. From bottom to top: (A) ice-rafted debris (IRD) concentrations (brown line), percentages of polar planktonic foraminifer [N. pachyderma (s)] (gray line), (B) winter SST estimates based on planktonic foraminifer associations, (C) pollen percentages of temperate trees (Acer, *Alnus*, *Betula*, *Corylus*, *Juniperus-Cupressus*-type, deciduous and evergreen *Quercus*, *Fraxinus excelsior*-type, *Salix*, *Tilia*, and *Ulmus*), and *Pinus*, (D) pollen percentages of (E) semi-desert plants and (F) *Calluna vulgaris*. Gray bars represent the HS 4, HS 3 (~GS 5), HS 2 (~GS 3) and HS 1 and GS 8, GS 7, GS 6, GS 4. HS 4, HS 3 and HS 1 are divided in two phases, a and b. HS 4 (a: ~39.5–38.6 ka; b: ~38.6–38 ka); HS 2 (a: ~26–25.5 ka; b: ~25.5–24 ka); HS 1 (a: ~18.8–17.5 ka; b: ~17.5–15.8 ka). Adapted from Naughton et al. (2009).
<figure><img src="image_17.png" /><figcaption>**FIGURE 10** | Mapping oceanic and atmospheric patterns of Heinrich events in the North Atlantic region within: (A) the first phase and (B) the second phase of HS 4, HS 2 and HS 1. The gray zone in the middle of the North Atlantic represents the Ruddiman belt. Adapted from Naughton et al. (2009).</figcaption></figure>

of the westerlies and increased precipitation in southwestern Europe. The scenario highlighted by our synthesis suggests that a mechanism similar to that of the present-day NAO in negative mode could have played a role (Naughton et al., 2009).

During the second phase of HS4, HS2 and HS1, climatic proxies indicate a northward migration of the thermal front associated with the northward shift of the westerlies and dryness in western Europe. This is accompanied in simulations by the southward displacement of the Inter Tropical Convergence Zone (ITCZ) (Kageyama et al., 2009) as detected from the increased of trade winds offshore Senegal (Jullien et al., 2007). As a whole, this configuration is similar to that currently produced by the positive mode of the NAO (Naughton et al., 2009). A marine pollen record from the Bay of Bengal indicates that HS 2 was associated with extremely weak summer and winter Indian monsoons, thus adding support for the southward migration of the ITCZ (Zorzi et al., 2015) during this period as also shown by South American and Asian speleothem records (e.g., Wang et al., 2004, 2008) and other marine sequences off western Africa (e.g., Weldeab et al., 2007).

The multi-phase structure of HS 4, HS 2 and HS 1 is supported by changes in ¹⁷O-excess, increase in CO₂ and methane mixing ratio and heavier δD-CH₄ and δ¹⁸O_atm observed in Greenland ice cores (Guillevic et al., 2014). In particular, in the second part of GS 9 (39.9–38.1 ka) the HS 4 imprint in Greenland would be characterized by a lower-latitude signal (without no changes in high-latitude temperature) associated with the southward migration of the ITCZ and a cold Europe (Guillevic et al., 2014). These data evidence a decoupling during GS 9 between stable cold Greenland temperature and low-latitude climate variability.

A decoupling between high and lower latitudes is also observed during GI 12 (46.8–44.3 ka) and GI 8 (38.2–36.6 ka)

(Sánchez Goñi et al., 2009). After the northern hemispheric rapid warming at the GS-GI transition, the trend during the first part of the GI is a Greenland cooling and an Iberian warming. This increase of the North Atlantic climatic gradient led to moisture transportation to Greenland from mid latitudes (lightest d-excess) and to a drying episode in Iberia. The subsequent temperature decrease in Greenland and Iberia associated with the precipitation increase in the latter region occurred when the major source of Greenland precipitation shifted to lower latitudes (d-excess increase). These examples of decoupling during GS and GI provide new targets for benchmarking climate model simulations and testing mechanisms associated with millennial variability (Guillevic et al., 2014).

## INTERACTIONS BETWEEN MILLENNIAL AND ORBITAL CLIMATE VARIABILITY: THE MISSING PIECE OF THE ICE AGE PUZZLE?

Changes in insolation roughly explain the timing of deglaciation and glacial inception (Milankovitch, 1941; Shackleton and Opdyke, 1973; Hays et al., 1976). However, internal feedback processes are needed to explain the ice age cycles, and previous research has mostly focused on those loops related to atmospheric and oceanic circulation (e.g., Risebrobakken et al., 2007; Bar-Or et al., 2008), GHG concentrations (e.g., Ganopolski et al., 2016), and vegetation-albedo feedback (Crucifix and Loutre, 2002; Desprat et al., 2005; Sánchez Goñi et al., 2005). So far little attention has been paid to the interactions between shorter (millennial-scale) and longer (10,000 to 100,000 years) term climate variability that may amplify the original orbital forcing (Hodell, 2016). Another important missing piece of the Ice Age puzzle involves the still debated origin of the
<figure><img src="image_18.png" /></figure>

**FIGURE 11** | Insolation changes versus marine and terrestrial palaeoclimatic records from core MD04-2845 for the interval 85–50 ka (Sánchez Goñi et al., 2013): (A) insolation curve at 65°N in June (Berger, 1978), and δ¹⁸O benthic foraminifer isotopic record (Cibicides wuellestorfi > 150 µm), (B) pollen records of Atlantic and Boreal forests in western France, (C) foraminifer-based summer SST (red line) in the Bay of Biscay, and percentages of temperate and subtropical foraminifer assemblages (blue line), and (D) IRD concentration and *N. pachyderma* (s) percentage (black solid line) records.

short-term climatic variability at high versus low latitudes, which has implications for the start of glaciations (Berger et al., 2006).

## Internal Mechanisms Linking Changes in Earth's Orbit to Ice Ages

In this section, we firstly focus on the specific mechanisms linking changes in Earth's orbit to Ice Ages taking as a test bed the onset of the last glaciation (i.e., MIS 5a/4 transition: ~80–70 ka), when one of the largest ice accumulations of the last 250,000 years occurred.

Ruddiman and McIntyre (1979) showed unequivocally that the ice growth in the northern hemisphere was contemporaneous with persistent warmth and high salinity in the subpolar North Atlantic (44–54° N) when boreal summer insolation was decreasing. This apparent paradox was interpreted as the development of a strong thermal gradient between a warm subpolar North Atlantic and a cold nearby land leading to an increase of moisture. This moisture was transported toward the north by storm tracks and fell as snow following the decrease in insolation. However, this theoretical model involving a warm ocean-cold land thermal contrast was not confirmed at that time due to the lack of data. Also, the interglacial/glacial
<figure><img src="image_19.png" /></figure>

**FIGURE 12** | (Top) A “zoom” on the interval 84.2–68 ka. Greenland temperature record compared with land-sea palaeoclimatic records from the two western European margin cores during the MIS 5a/4 transition (Sánchez Goñi et al., 2013). (A) Greenland temperature record (Landais et al., 2004), (B) δ¹⁸O curve of benthic foraminifer from core MD04-2845 (dark blue), and reconstructed sea level changes (light blue) (Waelbroeck et al., 2002), (C) core MD04-2845 (Bay of Biscay): foraminifer-based summer SST curve (red) with the minima and maxima values found in the set of the five selected analogs (gray surface), Atlantic forest pollen percentage curves (green), (D) MD99-2331 (NW Iberian margin) the same as core MD04-2845, and UK’37-based SST (black). Gray bands indicate warm intervals in western France. C20, C19, and C18’: cooling events in the western North Atlantic Ocean. Gl 21, Gl 20, Gl 19, and Gl 18: warm phases in Greenland. (Bottom) Schematic representation of the cold air-warm ocean contrast in the western European margin at the time of moderate iceberg discharges in the western North Atlantic Ocean.

transitions were then thought to be gradual but we know now that they were punctuated by millennial-scale events. Thus, one may wonder whether sub-orbital climatic variability, i.e., D-O cycles and HEs, affected this orbitally-controlled ice growth.

During the MIS 5a/4 transition, the direct land-sea correlation of the core MD04-2845, located at 45°N in the northern subtropical Atlantic, shows that the increase in ice volume

at orbital scale was contemporaneous with a progressive cooling in western Europe, as indicated by the replacement of temperate forest by boreal forest (mainly *Abies* and *Picea*). The simultaneous stabilization of warm SST in the Bay of Biscay demonstrates for the first time a long-term increase of the thermal warm sea-cold land gradient toward the MIS 4 glacial (Sánchez Goñi et al., 2013) (Figure 11). At sub-orbital scale, a decoupling with cooling on land synchronous with warming in SST is furthermore observed (Figure 12). Three successive episodes of strong land-sea thermal gradient and another three with weak gradients occurred during the MIS 5a/4 transition, superimposed on the long-term increasing trend in the thermal gradient (Sánchez Goñi et al., 2013).

The strong land-sea thermal gradients coincided with the North Atlantic cold events C20, C19, and C18' marked by moderate iceberg discharges in the western side of the North Atlantic (McManus et al., 1994) and related to Greenland cold episodes GS 21, GS 20, and GS 19 (Rasmussen et al., 2014) while the weak gradients coincided with the warm episodes in the North Atlantic and Greenland, GI 20, GI 19, and GI 18. Whereas the events C20, C19, and C18' were cold in the subpolar gyre and cold conditions installed in western Europe, the northern subtropical gyre remained warm. Increased snowfall in northern Europe and subsequent ice growth resulted from the high rates of moisture production resulting from the abovementioned strong land-sea thermal gradient and its transport by northward tracking storms. Southward displacement of tundra by 10° in latitude during cold phases C20, C19 and C18', as suggested by boreal forest colonization of western Europe, probably amplified ice growth owing to the increase in surface albedo (Crucifix and Loutre, 2002). Weak gradients slowed down the process but still allowed ice accumulation. The marine and pollen palaeoclimatic records from core MD99-2331, off north-western Iberia, show the same orbital and sub-orbital increases in the land-sea thermal gradient, excluding during C19 (Figure 12). In contrast, further south, in the south-western Iberian margin, the thermal gradient remained weak throughout the MIS 5a/4 transition. In summary, the direct comparison between marine and terrestrial records from the European margin demonstrates for the first time a long-term increase in the thermal gradient between the cold air and warmer sea, and three short intervals of even more pronounced thermal gradients during the last entering in glaciation. This synergy between orbital and millennial-scale variability provided a substantial source of moisture that was transported, through northward-tracking storms, to feed ice sheets in colder Greenland, northern Europe, and the Arctic.

## The Origin of the Millennial-Scale Variability and the Implication for the Ice Ages

The origin of the short-lived climate changes is still a subject of debate. On the one hand, some authors propose that high-latitude processes such as iceberg discharges and changes
<figure><img src="image_20.png" /></figure>

**FIGURE 13** | (Left) Vegetation changes versus the δ¹⁸O_b record from Site U1385 for MIS 19. (A) Percentages of *Isoetes* spores (cyan) and of semi-desert plants (*Ephedra distachya*-type, *E. fragilis*-type, Chenopodiaceae that now is included in the Amaranthaceae pollen morphotype, *Artemisia*) in orange, pollen percentages of the Mediterranean taxa (evergreen *Quercus*, *Olea*, *Pistacia*, *Phillyrea*, *Cistus*) in red and Mediterranean forest in light green (mainly deciduous *Quercus*, and Mediterranean taxa). Dark green curve: three-point weigθed-average smoothing, (B) Uk'37-based SST record, (C) δ¹⁸O record from benthic foraminifer (δ¹⁸Ob) *Cibicidoides wuellerstorfii*, Black curve: three-point weighted-average smoothing, (D) changes in insolation at 65°N in July (in black), in the precession index (e*·sin ω) (in red) and obliquity (light blue) (Berger and Loutre, 1991). Orange and gray vertical panels indicate the Tajo interglacial and the MIS 19b, respectively. Long dashed lines indicate the MIS 19 boundaries. Black arrows indicate Mediterranean forest contractions. Light gray bands indicate the intervals of relatively low Mediterranean forest cover that coincide with heavier benthic δ¹⁸O values. (Right) Schematic representation of the vegetation in the Iberian Peninsula, atmospheric and surface ocean circulations and temperature gradient in the North Atlantic during millennial events of MIS 19c and MIS 11c. Adapted from Desprat et al. (2017).

in the AMOC would have forced the rapid climate changes observed in the mid and low latitudes of the North Atlantic Ocean throughout the Pleistocene (Ganopolski and Rahmstorf, 2001). On the other hand, other authors put forward a low-latitude forcing mechanism via the harmonics of precession. Actually, the tropical latitudes receive over the course of the year a daily irradiance characterized by a double maximum that originates in the equatorial insolation cyclicities at 5,500 and 11,000 years (Berger et al., 2006). A direct consequence of this process would be a larger latitudinal thermal gradient and the enhanced transport of warmth and moisture by either atmospheric (westerlies) or oceanic circulation (subtropical gyre) from equatorial to high latitudes in the North Atlantic (Berger et al., 2006).

This low latitude forcing has also been suggested during periods of low eccentricity, when precession changes are

muted, such as the last 45,000 years and MIS 19, centered at 800,000 years ago, based only on paleoceanographic evidence (e.g., McIntyre and Molfino, 1996; Ferretti et al., 2015). Our recent land-sea direct correlation studies at the IODP Site U1385 (south-western Iberian margin) have provided firm evidence on the transport of energy from equatorial regions to the high latitudes by means of the westerlies during the MIS 19 and MIS 11, the astronomically closest analogs to the present interglacial. During MIS 19 and MIS 11, Mediterranean forest pollen percentages indicate two and three long-term major Mediterranean forest expansions, respectively, following maxima in insolation (Figure 13 for MIS 19 as an example). These major forest episodes point to the occurrence of a well-established Mediterranean climate, with higher winter precipitations and warmer and drier summers. Forest expansion in south-western Iberia therefore indicates that the prevailing
zonal configuration of the westerlies during these interglacials brought precipitation to this region (Oliveira et al., 2016; Sánchez Goñi et al., 2016b). During MIS 19c and MIS 11c, a long-term decrease in the Mediterranean forest cover indicating a progressive cooling and drying trend contrasts with the SST record showing quite stable and warm conditions (around 18°C) in the subtropical gyre. During these intervals, freshwater input is not detected on this margin and the progressive increase in the local benthic foraminifer δ¹³C indicates good deep water ventilation (Oliveira et al., 2016; Sánchez Goñi et al., 2016b). These long-term trends were punctuated by major contractions of the Mediterranean forest coeval with increases in Mediterranean semi-desert plants, which may indicate intervals when the westerlies were slightly deflected toward the north with a subsequently reduced influence in south-western Iberia. Repeated meridional shifts in the westerlies during MIS 19c and MIS 11c would have implied successive large masses of warm and moist air reaching the high latitudes of the North Atlantic via the western boundary current of the North Atlantic subtropical gyre. At sub-orbital time scales a decoupling also occurred during MIS 19c and MIS 11c with several drying and cooling events on land concomitant with stable warm SST.

The direct comparison of marine and pollen proxies shows therefore a clear air-sea decoupling during MIS 19c and MIS 11c at orbital and millennial time-scales. In contrast, the other cool and dry events during MIS 19b-19a and MIS 11b-11a were contemporaneous with SST cooling in the subtropical gyre and IRD deposition at in the subpolar gyre (IODP Site U1314 and ODP Site 980), which evidences the input of large amounts of icebergs in the subpolar North Atlantic (Oppo et al., 1998; Alonso-Garcia et al., 2011). Interestingly, the atmospheric cool and dry events during MIS 19c and MIS 11c were not associated with freshwater pulses in either the subtropical or subpolar gyres indicating that the cause of these events was not primarily related to high latitude ice-sheet dynamics. Fourier spectral analysis applied to the pollen and marine records, excluding the benthic foraminifer δ¹⁸O record, show dominant 5,000- and 10,000-year cyclicities for MIS 19 and MIS 11, likely related with the fourth and second harmonics of precession, respectively (Oliveira et al., 2016; Sánchez Goñi et al., 2016b; Oliveira, 2017). For MIS 19c, the comparison of the Mediterranean forest pollen record and its 5-kyr bandpass filter output curve with the variations in the largest amplitude of the seasonal cycle reconstructed at the equator shows good correspondence between the two records although the magnitude of changes was larger during MIS 19b-19a than during MIS 19c (Sánchez Goñi et al., 2016b). This good correspondence gives support to the low latitude origin of the observed repeated shift of the westerlies in southwestern Iberia at 37°N likely related to the harmonics of precession.

The stronger forest contractions during MIS 19b-19a and MIS 11b-11a would be explained by the regional SST cooling linked to the freshwater fluxes arriving at the subpolar gyre and reducing the AMOC. These data show once more that cold and dry episodes occurred independently of the amount

of ice volume (Oppo et al., 1998; Desprat et al., 2009). In contrast, their intensity and duration are modulated by positive feedback mechanisms on the AMOC associated with ice dynamics and, particularly, with freshwater pulses from iceberg discharges and melting (Barker et al., 2015; Oliveira et al., 2016).

During the MIS 19c/19b transition the three cold land-warm sea decoupling events coincide with three increases in the benthic δ¹⁸O values concomitant with decreases in sea level. This observation suggests that for a period of low amplitude changes in insolation a tropically-driven strong warm sea-cold land contrast at millennial-timescale in the western European margin likely contributed to the increase of moisture. Similarly to that described previously for the MIS 5a/4 transition, these repeated millennial-scale episodes of moisture increase fed the ice caps through northward storm tracks and triggered the onset of the successive glacial period, at around 774,000 years ago.

## CONCLUSION

Pollen analysis from deep-sea sedimentary sequences constitutes a powerful tool for reconstructing regional vegetation and climate changes that, in turn, influenced the global climate. One of the main advantages of pollen records from marine sediment sequences is that they allow direct correlations of climate change over land, in the ocean and in the ice domain, with minimum chronological uncertainty. This is because marine proxies (e.g., δ¹⁸O from foraminifer, dinoflagellate cysts, IRD) and terrestrial tracers (e.g., pollen) can be analyzed from the same sediment sample so ocean-ice-atmosphere paleoreadings can be obtained for the same time slice. A second crucial advantage is that besides pollen records, very few other palaeoclimate markers exist that are sensitive to the seasonality of temperature and precipitation. The pollen richness of marine sediments allows the quantitative estimation of temperatures, precipitation and, more importantly, their seasonality, one of the key parameters for understanding climate change. However, further research is necessary for improving the quantification of key climatic parameters, particularly seasonality in subtropical and tropical regions, and of vegetation cover and composition to better estimate albedo changes and their impact on global climate.

This overview has shown that European vegetation and climate responded to long-term and shorter-term climate changes. There was a dynamic equilibrium between vegetation and climate for short periods of forcing such as the D-O cycles and HEs that were similar in magnitude and velocity to the present-day global warming. However, the magnitude of the millennial-scale changes of the last glacial period was regionally-specific. Regions below 40°N appear to have been imprinted by precession, through changes in the amplitude of seasonal contrasts and particularly winter precipitation, while the amplitude of warming in northern regions seems to have been modulated by obliquity that determines the annual temperatures. Within GS and GI, particularly during HS 4, GI 12 and GI 8, a
decoupling between high and lower latitudes is observed, a decoupling that has also been identified for the amplitude of the European interglacials. In the mid-latitude regions above 40°N the forest expansion maxima, i.e., the warmest peaks, of the last five terrestrial interglacials is similar while below 40°N their magnitude differed. In this latter region, the main factor for forest expansion is winter precipitation, and this climatic parameter depends on precession.

The overview presented here also suggests that millennial scale climatic variability likely played an important role in glacial inception whatever the boundary climatic conditions were. During one of the largest and fastest ice growth phases at the MIS 5a/4 transition, the millennial variability associated with moderate iceberg discharges in the western North Atlantic Ocean pushed the warm and saline Gulf Stream toward the European margin and created optimal conditions to develop repeated and pronounced air-sea thermal contrasts in this region. The synergy between low and high frequency climate changes amplified moisture production and snow fall in the high latitudes of the Northern Hemisphere. Within interglacials MIS 11 and MIS 19, the cyclicity of the harmonics of precession, repeatedly warming the low latitudes and increasing the latitudinal thermal gradient, created similar pronounced increases of the air-sea thermal gradient in the European margin that favored moisture production and the onset of ice growth. The theoretical model proposed here, and therefore the physical mechanisms involving rapid changes in the westerlies, in the direction of the subtropical gyre and in the air-sea thermal contrast superimposed to the long-term change, needs to be evaluated by model experiments. In particular, some authors (e.g., Crucifix et al., 2016) propose the integration of stochastic models in the Earth Model of Intermediate Complexity with the aim to analyze the effect of the rapid, millennial scale, variability on the long-term glacial-interglacial climate changes.

## REFERENCES

Allen, J. R. M., Brandt, U., Brauer, A., Hubberten, H.-W., Huntley, B., Keller, J., et al. (1999). Rapid environmental changes in southern Europe during the last glacial period. *Nature* 400, 740–743. doi: 10.1038/23432

Alonso-Garcia, M., Sierro, F. J., Kucera, M., Flores, J. A., Cacho, I., and Andersen, N. (2011). Ocean circulation, ice sheet growth and interhemispheric coupling of millennial climate variability during the mid-Pleistocene (ca 800–400 ka). *Quat. Sci. Rev.* 30, 3234–3247. doi: 10.1016/j.quascirev.2011.08.005

Alvarez-Solas, J., and Ramstein, G. (2011). On the triggering mechanism of Heinrich events. *Proc. Natl. Acad. Sci. U.S.A.* 108, E1359–E1360. doi: 10.1073/pnas.1116575108

Augustin, L., Barbante, C., Barnes, P. R., Barnola, J. M., Bigler, M., and Castellano, E. (2004). Eight glacial cycles from an Antarctic ice core. *Nature* 429, 623–628. doi: 10.1038/nature02599

Bailey, R. G. (1998). *Ecoregions: The Ecosystem Geography of the Oceans and Continents*. New York, NY: Springer.

Barker, S., Chen, J., Gong, X., Jonkers, L., Knorr, G., and Thornalley, D. (2015). Icebergs not the trigger for North Atlantic cold events. *Nature* 520, 333–336. doi: 10.1038/nature14330

Barker, S., Knorr, G., Edwards, R. L., Parrenin, F., Putnam, A. E., Skinner, L. C., et al. (2011). 800,000 years of abrupt climate variability. *Science* 334, 347–351. doi: 10.1126/science.1203580

## AUTHOR CONTRIBUTIONS

MFSG designed the synthesis and wrote the manuscript. All the co-authors contributed to the writing of the manuscript.

## FUNDING

Financial support was provided by WarmClim, a LEFE-INSU, IMAGO project. CM-M was supported by an IdEx Bordeaux post-doctoral fellowship, FN by Climhol (PTDC/AAC-CLI/100157/2008) and Ultimatum (IF/01489/2015), DO by the Portuguese Foundation for Science and Technology (FCT) through the doctoral grant (SFRH/BD/9079/2012).

## ACKNOWLEDGMENTS

A substantial part of this review is based on marine sedimentary cores that were collected during the IMAGES I, GEOSCIENCES, and ALIENOR oceanographic cruises, and IODP Expedition 339. The authors are grateful to the coring and logistic teams inboard of the R/V Marion Dufresne II and Joides Resolution. They are grateful to Vincent Hanquiez for helping with the drawing of the figures.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.00038/full#supplementary-material

<caption>TABLE S1 | List of the 129 deep-sea pollen records distributed over the world.</caption>

Bar-Or, R., Erlick, C., and Gildor, H. (2008). The role of dust in glacial-interglacial cycles. *Quat. Sci. Rev.* 27, 201–208. doi: 10.1016/j.quascirev.2007.10.015

Bazin, L., Landais, A., Lemieux-Dudon, B., Toyé Mahamadou Kele, H., Veres, D., Parrenin, F., et al. (2013). An optimized multi-proxy, multi-site Antarctic ice and gas orbital chronology (AICC2012): 120-800 ka. *Clim. Past* 9, 1715–1731. doi: 10.5194/cp-9-1715-2013

Beaudouin, C., Suc, J.-P., Escarguel, G., Arnaud, M., and Charmasson, S. (2007). The significance of pollen signal in present-day marine terrigenous sediments: the example of the Gulf of Lions (Western Mediterranean Sea). *Geobios* 40, 159–172. doi: 10.1016/j.geobios.2006.04.003

Bennett, K. D., and Willis, K. J. (1995). The role of ecological factors in controlling vegetation dynamics on long temporal scales. *G. Bot. Ital.* 129, 243–254. doi: 10.1080/11263509509436132

Berger, A. (1978). Long-term variations of daily insolation and Quaternary climatic changes. *J. Atmos. Sci.* 35, 2362–2367. doi: 10.1175/1520-0469(1978)035<2362: LTVODI>2.0.CO;2

Berger, A., and Loutre, M. F. (1991). Insolation values for the climate of the last 10 million years. *Quat. Sci. Rev.* 10, 297–317. doi: 10.1016/0277-3791(91)90033-Q

Berger, A., and Loutre, M. F. (2004). Astronomical theory of climate change. *J. Phys. IV France* 121, 1–35. doi: 10.1051/jp4:2004121001

Berger, A., Loutre, M. F., and Mélice, J. L. (2006). Equatorial insolation: from precession harmonics to eccentricity frequencies. *Clim. Past* 2, 131–136. doi: 10.5194/cp-2-131-2006
Birks, H. J. B. (1981). "The use of pollen analysis in the reconstruction of past climates: a review," in *Climate and History*, eds T. M. L. Wigley, M. J. Ingram, and G. Farmer (Cambridge: Cambridge University Press), 111–138.

Blanco Castro, E., Casado González, M. A., Costa Tenorio, M., Escribano Bombin, R., García Antón, M., Génova Fuster, M., et al. (1997). *Los Bosques Ibéricos*. Barcelona: Planeta, 572.

Bond, G., Broecker, W., Johnsen, S., McManus, J., Labeyrie, L., Jouzel, J., et al. (1993). Correlations between climate records from North Atlantic sediments and Greenland ice. *Nature* 365, 143–147. doi: 10.1038/365143a0

Bond, G., and Lotti, R. (1995). Iceberg discharges into the North Atlantic on millennial time scales during the Last Glaciation. *Science* 267, 1005–1009. doi: 10.1126/science.267.5200.1005

Bordon, A., Peyron, O., Lézine, A.-M., Brewer, S., and Fouache, E. (2009). Pollen-inferred Late-Glacial and Holocene climate in southern Balkans (Lake Maliq). *Quat. Int.* 200, 19–30. doi: 10.1016/j.quaint.2008.05.014

Cambon, G., Suc, J.-P., Aloisi, J.-C., Giresse, P., Monaco, A., Touzani, A., et al. (1997). Modern pollen deposition in the Rhone delta area (Lagoonal and marine sediments), France. *Grana* 36, 105–113. doi: 10.1080/0017319709362596

Carré, M., and Cheddadi, R. (2017). Seasonality in long-term climate change. *Quaternaire* 28, 173–177. doi: 10.4000/quaternaire.8018

Chabaud, L., Sánchez Goñi, M. F., Desprat, S., and Rossignol, L. (2014). Land-sea climatic variability in the eastern North Atlantic subtropical region over the last 14,200 years: atmospheric and oceanic processes at different timescales. *Holocene* 24, 787–797. doi: 10.1177/0959683614530439

Chapman, M. R., and Shackleton, N. J. (1999). Global ice-volume fluctuations, North Atlantic ice-rafted events, and deep-ocean circulation changes between 130 and 70 ka. *Geology* 27, 795–798. doi: 10.1130/0091-7613(1999)027<0795:GIVFNA>2.3.CO;2

Claussen, M. (2007). "Introduction to climate forcing and climate feedbacks," in *The Climate of Past Interglacials: Developments in Quaternary Science*, Vol. 7, eds F. Sirocko, M. Claussen, M. F. Sánchez Goñi, and T. Litt (Amsterdam: Elsevier), 3–11. doi: 10.1016/S1571-0866(07)80026-1

Combourieu Nebout, N., Peyron, O., Dormoy, I., Desprat, S., Beaudoin, C., Kothoff, U., et al. (2009). Rapid climatic variability in the west Mediterranean during the last 25 000 years from high resolution pollen data. *Clim. Past* 5, 503–521. doi: 10.5194/cp-5-503-2009

Crucifix, M., and Loutre, M. F. (2002). Transient simulations over the last interglacial period (126-115 kyr BP): feedback and forcing analysis. *Clim. Dyn.* 19, 417–433. doi: 10.1007/s00382-002-0234-z

Crucifix, M., Mitsui, T., and Lenoir, G. (2016). "Challenges for ice age dynamics: a dynamical systems perspective," in *Nonlinear and Stochastic Climate Dynamics*, eds C. Franzke and T. J. O'Kane (Cambridge: Cambridge University Press), 1–32.

Daniau, A. L., Sánchez Goñi, M.-F., Beaufort, L., Laggoun-Défarge, F., Loutre, M.-F., and Duprat, J. (2007). Dansgaard-Oeschger climatic variability revealed by fire emissions in southwestern Iberia. *Quat. Sci. Rev.* 26, 1369–1383. doi: 10.1016/j.quascirev.2007.02.005

Daniau, A.-L., Sánchez Goñi, M. F., Martinez, P., Urrego, D. H., Bout-Roumazeilles, V., Desprat, S., et al. (2013). Orbital-scale climate forcing of grassland burning in southern Africa. *Proc. Natl. Acad. Sci. U.S.A.* 110, 5069–5073. doi: 10.1073/pnas.1214292110

Dansgaard, W., Johnsen, S., Clausen, H. B., Dahl-Jensen, D., Gundestrup, N., Hammer, C. U., et al. (1984). "North Atlantic climatic oscillations revealed by deep Greenland ice cores," in *Climate Processes and Climate Sensitivity*, eds J. E. Hansen and T. Takahashi (Washington, DC: American Geophysical Union), 288–298.

de Beaulieu, J.-L., Andrieu-Ponel, V., Reille, M., Grüger, E., Tzedakis, C., and Svobodova, H. (2001). An attempt at correlation between the Velay pollen sequence and the Middle Pleistocene stratigraphy from central Europe. *Quat. Sci. Rev.* 20, 1593–1602. doi: 10.1016/S0277-3791(01)00027-0

deMenocal, B. (2004). African climate change and faunal evolution during the Plio-Pleistocene. *Earth Planet. Sci. Lett.* 2004, 3–24. doi: 10.1016/S0012-821X(04)00003-2

Desprat, S., Diaz Fernandez, M., Coulon, T., Ezzat, L., Pessarossi-Langlois, J., Gil, L., et al. (2015). Pinus nigra (European black pine) as the dominant species of the last glacial pinewoods in south-western to central Iberia: a morphological study of modern and fossil pollen. *J. Biogeogr.* 42, 1998–2009. doi: 10.1111/jbi. 12566

Desprat, S., Oliveira, O., Naughton, F., and Sánchez Goñi, M. F. (2017). L'étude du pollen des séquences sédimentaires marines pour la compréhension du climat: l'exemple des périodes chaudes passées. *Quaternaire* 28, 259–269. doi: 10.4000/quaternaire.8102

Desprat, S., Sánchez Goñi, M. F., McManus, J. F., Duprat, J., and Cortijo, E. (2009). Millennial-scale climatic variability between 340 000 and 270 000 years ago in SW Europe: evidence from a NW Iberian margin pollen sequence. *Clim. Past* 5, 53–72. doi: 10.5194/cp-5-53-2009

Desprat, S., Sánchez Goñi, M. F., Naughton, F., Turon, J. L., Duprat, J., Malaizé, B., et al. (2007). "Climate variability of the last five isotopic interglacials from direct land-sea-ice correlation," in *The Climate of Past Interglacials*, eds F. Sirocko, M. Claussen, M. F. Sánchez Goñi, and T. Litt (Amsterdam: Elsevier), 375–386.

Desprat, S., Sánchez Goñi, M. F., Turon, J.-L., Duprat, J., Malaizé, B., and Peypouquet, J.-P. (2006). Climatic variability of Marine Isotope Stage 7: direct land-sea-ice correlation from a multiproxy analysis of a north-western Iberian margin deep-sea core. *Quat. Sci. Rev.* 25, 1010–1026. doi: 10.1016/j.quascirev. 2006.01.001

Desprat, S., Sánchez Goñi, M. F., Turon, J.-L., McManus, J. F., Loutre, M. F., Duprat, J., et al. (2005). Is vegetation responsible for glacial inception during periods of muted insolation changes? *Quat. Sci. Rev.* 24, 1361–1374. doi: 10.1016/j.quascirev.2005.01.005

Dupont, L. (2011). Orbital scale vegetation change in Africa. *Quat. Sci. Rev.* 30, 3589–3602. doi: 10.1016/j.quascirev.2011.09.019

Dupont, L. M., and Wyputta, U. (2003). Reconstructing pathways of aeolian pollen transport to the marine sediments along the coastline of SW Africa. *Quat. Sci. Rev.* 22, 157–174. doi: 10.1016/S0277-3791(02)00032-X

Elliot, M., Labeyrie, L., Dokken, T., and Manthé, S. (2001). Coherent patterns of ice-rafted debris deposits in the Nordic regions during the last glacial (10-60 ka). *Earth Planet. Sci. Lett.* 194, 151–163. doi: 10.1016/S0012-821X(01)00561-1

Ferretti, P., Crowhurst, S. J., Naafs, B. D. A., and Barbante, C. (2015). The marine isotope stage 19 in the mid-latitude North Atlantic Ocean: astronomical signature and intra-interglacial variability. *Quat. Sci. Rev.* 108, 95–110. doi: 10.1016/j.quascirev.2014.10.024

Fletcher, W. J., and Sánchez Goñi, M. F. (2008). Orbital- and sub-orbital-scale climate impacts on vegetation of the western Mediterranean basin over the last 48,000 yr. *Quat. Res.* 70, 451–464. doi: 10.1016/j.yqres.2008.07.002

Fletcher, W. J., Sánchez Goñi, M. F., Allen, J. R. M., Cheddadi, R., Combourieu-Nebout, N., Huntley, B., et al. (2010). Millennial-scale variability during the last glacial in vegetation records from Europe. *Quat. Sci. Rev.* 29, 2839–2864. doi: 10.1016/j.quascirev.2009.11.015

Follieri, M., Magri, D., and Sadori, L. (1988). 250.000-year pollen record from valle di Castiglione (Roma). *Pollen Spores* 30, 329–356.

Ganopolski, A., and Rahmstorf, S. (2001). Rapid changes of glacial climate simulated in a coupled climate model. *Nature* 409, 153–158. doi: 10.1038/ 35051500

Ganopolski, A., Winkelmann, R., and Schellnhuber, H. J. (2016). Critical insolation-CO₂ relation for diagnosing past and future glacial inception. *Nature* 529, 200–203. doi: 10.1038/nature16494

Gouveia, C., Trigo, R. M., DaCamara, C. C., Libonati, R., and Pereira, J. M. C. (2008). The North Atlantic Oscillation and European vegetation dynamics. *Int. J. Climatol.* 28, 1835–1847. doi: 10.1002/joc.1682

Groot, J. J., and Groot, C. R. (1966). Marine palynology: possibilities, limitations, problems. *Mar. Geol.* 4, 387–395. doi: 10.1016/0025-3227(66) 90007-7

Guillevic, M., Bazin, L., Landais, A., Stowasser, C., Masson-Delmotte, V., Blunier, T., et al. (2014). Evidence for a three-phase sequence during Heinrich Stadial 4 using a multiproxy approach based on Greenland ice core records. *Clim. Past* 10, 2115–2133. doi: 10.5194/cp-10-2115-2014

Guiot, J. (1990). Methodology of the last climatic cycle reconstruction from pollen data. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* 80, 49–69. doi: 10.1016/0031-0182(90)90033-4

Guiot, J., and De Vernal, A. (2007). "Transfer functions: methods for quantitative paleoceanography based on microfossils," in *Developments in Marine Geology*, Vol. 1, eds C. Hillaire-Marcel and A. de Vernal (Amsterdam: Elsevier), 523–563.

Hays, J. D., Imbrie, J., and Shackelton, N. J. (1976). Variations in the Earth's orbit pacemaker of the ice ages. *Science* 194, 1121–1132. doi: 10.1126/science.194. 4270.1121
Heinrich, H. (1988). Origin and consequences of cyclic ice rafting in the northeast Atlantic Ocean during the past 130,000 years. *Quat. Res.* 29, 142–152. doi: 10.1016/0033-5894(88)90057-9

Heusser, L. E., and Balsam, W. L. (1977). Pollen distribution in the N.E. Pacific Ocean. *Quat. Res.* 7, 45–62. doi: 10.1016/0033-5894(77)90013-8

Heusser, L. E., and Van de Geer, G. (1994). Direct correlation of terrestrial and marine paleoclimatic records from four glacial-interglacial cycles-DSDP site 594 Southwest Pacific. *Quat. Sci. Rev.* 13, 273–282. doi: 10.1016/0277-3791(94) 90030-2

Heusser, L. E., and Vm, Bryant Jr. (1985). "Quaternary palynology of marine sediments in the northeast Pacific, northwest Atlantic, and Gulf of Mexico," in *Pollen Records of Late-Quaternary North American Sediments*, ed. R. G. Holloway (Dallas, TX: AASP Foundation), 385–403.

Heusser, L. E. (1986/1987). Pollen in marine cores: evidence of past climates. *Oceanus* 29, 64–70.

Hodell, D. A. (2016). The smoking gun of the ice ages. *Science* 354, 1235–1236. doi: 10.1126/science.aal4111

Hodell, D. A., Lourens, L., Crowhurst, S., Konijnendijk, T., Tjallingii, R., Jiménez-Espejo, F., et al. (2015). A reference time scale for site U1385 (Shackleton Site) on the SW Iberian Margin. *Glob. Planet. Change* 133, 49–64. doi: 10.1016/j.gloplacha.2015.07.002

Hoerling, M., Eischeid, J., Perlwitz, J., Quan, X. W., Zhang, T., and Pegion, P. (2012). On the increased frequency of Mediterranean drought. *J. Clim.* 25, 2146–2161. doi: 10.1175/JCLI-D-11-00296.1

Hooghiemstra, H., Lézine, A.-M., Leroy, S. A. G., Dupont, L., and Marret, F. (1988). Late Quaternary palynology in marine sediments: a synthesis of the understanding of pollen distribution patterns in the NW African setting. *Quat. Int.* 148, 29–44. doi: 10.1016/j.quaint.2005.11.005

Hooghiemstra, H., Stalling, H., Agwu, C. O. C., and Dupont, L. M. (1992). Vegetational and climatic changes at the northern fringe of the Sahara 250,000–5000 years BP: evidence from 4 marine pollen records located between Portugal and the Canary Islands. *Rev. Palaeobot. Palynol.* 74, 1–53. doi: 10.1016/0034-6667(92)90137-6

Hooghmiestra, H., Agwu, C. O. C., and Beug, H.-J. (1986). Pollen and spore distribution in recent marine sediments: a record of NW-African seasonal wind patterns and vegetation belts. *Meteor Forsch. Ergeb.* C 40, 87–135.

Jullien, E., Grousset, F. E., Malaize, B., Duprat, J., Sánchez Goni, M. F., Eynaud, F., et al. (2007). Tropical “dusty Heinrich-like events” vs. high-latitude “icy Heinrich events”. *Quat. Res.* 68, 379–386. doi: 10.1016/j.yqres.2007.07.007

Kageyama, M., Mignot, J., Swingedouw, D., Marzin, C., Alkama, R., and Marti, O. (2009). Glacial climate sensitivity to different states of the Atlantic Meridional Overturning Circulation: results from the IPSL model. *Clim. Past* 5, 551–570. doi: 10.5194/cp-5-551-2009

Knutti, R., Flückiger, J., Stocker, T. F., and Timmermann, A. (2004). Strong hemispheric coupling of glacial climate through freshwater discharge and ocean circulation. *Nature* 430, 851–856. doi: 10.1038/nature02786

Koreneva, E. V. (1966). Marine palynological researches in the U.S.S.R. *Mar. Geol.* 4, 565–574. doi: 10.1016/0025-3227(66)90016-8

Koreneva, E. V. (1971). "Spores and pollen in Mediterranean bottom sediments," in *The Micropaleontology of Oceans*, eds B. M. Funnel and W. R. Riedel (Cambridge: Cambridge University Press), 361–371.

Kühl, N., Gebhardt, C., Litt, T., and Hense, A. (2002). Probability density functions as botanical-climatological transfer functions for climate reconstruction. *Quat. Res.* 58, 381–392. doi: 10.1006/qres.2002.2380

Landais, A., Barnola, J. M., Masson-Delmotte, V., Jouzel, J., Chappellaz, J., Caillon, N., et al. (2004). A continuous record of temperature evolution over a whole sequence of Dansgaard-Oeschger during Marine Isotopic Stage 4 (76 to 62 kyr BP). *Geophys. Res. Lett.* 31, 101–113. doi: 10.1029/2004GL021193

Levavasseur, G., Vrac, M., Roche, D. M., and Paillard, D. (2012). Statistical modelling of a new global potential vegetation distribution. *Environ. Res. Lett.* 7:044019. doi: 10.1088/1748-9326/7/4/044019

Loulergue, L., Schilt, A., Spahni, R., Masson-Delmotte, V., Blunier, T., Lemieux, B., et al. (2008). Orbital and millennial-scale features of atmospheric CH4 over the past 800,000 years. *Nature* 453, 383–386. doi: 10.1038/nature06950

Luthi, D., Le Floch, M., Bereiter, B., Blunier, T., Barnola, J.-M., Siegenthaler, U., et al. (2008). High-resolution carbon dioxide concentration record 650,000–800,000 years before present. *Nature* 453, 379–382. doi: 10.1038/nature06949

Lynch-Stieglitz, J. (2017). The Atlantic Meridional Overturning Circulation and abrupt climate change. *Annu. Rev. Mar. Sci.* 9, 83–104. doi: 10.1146/annurev-marine-010816-060415

Martrat, B., Grimalt, J. O., Shackleton, N. J., De Abreu, L., Hutterli, M. A., and Stocker, T. F. (2007). Four cycles of recurring deep and surface water destabilizations on the Iberian margin. *Science* 317, 502–507. doi: 10.1126/science.1139994

McIntyre, A., and Molfino, B. (1996). Forcing of Atlantic Equatorial and subpolar millennial cycles by precession. *Science* 274, 1867–1870. doi: 10.1126/science.274.5294.1867

McManus, J. F., Bond, G. C., Broecker, W. S., Johnsen, S., Labeyrie, L., and Higgins, S. (1994). High-resolution climate records from the North Atlantic during the last interglacial. *Nature* 371, 326–329. doi: 10.1038/371326a0

McManus, J. F., Francois, R., Gherardi, J.-M., Keigwin, L. D., and Brown-Leger, S. (2004). Collapse and rapid resumption of Atlantic meridional circulation linked to deglacial climate changes. *Nature* 428, 834–837. doi: 10.1038/nature02494

McManus, J. F., Oppo, D. W., and Cullen, J. L. (1999). A 0.5-million-year record of millennial-scale climate variability in the North Atlantic. *Science* 283, 971–975. doi: 10.1126/science.283.5404.971

Milankovitch, M. M. (1941). *Kanon der Erdbestrahlung*. Beograd: Köninglich Serbsche Akademie.

Montade, V., Combourieu-Nebout, N., Kissel, C., Haberle, S. H., Siani, G., and Michel, E. (2013). Vegetation and climate changes during the last 22000 yr from a marine core near Taitao Peninsula, southern Chile. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* 369, 335–348. doi: 10.1016/j.palaeo.2012.11.001

Mudie, J., and McCarthy, F. M. G. (2006). Marine palynology: potentials for onshore-offshore correlation of Pleistocene-Holocene records. *Trans. R. Soc. S. Afr.* 61, 139–157. doi: 10.1080/00359190609519964

Naughton, F., Sánchez Goñi, M. F., Desprat, S., Turon, J.-L., Duprat, J., Malaizé, B., et al. (2007). Present-day and past (last 25000 years) marine pollen signal off western Iberia. *Mar. Micropaleontol.* 62, 91–114. doi: 10.1016/j.marmicro.2006.07.006

Naughton, F., Sánchez Goñi, M. F., Kageyama, M., Bard, E., Duprat, J., Cortijo, E., et al. (2009). Wet to dry climatic trend in north western Iberia within Heinrich events. *Earth Planet. Sci. Lett.* 284, 329–342. doi: 10.1016/j.epsl.2009.05.001

Naughton, F., Sánchez Goñi, M. F., Rodrigues, T., Salgueiro, E., Costas, S., Desprat, S., et al. (2016). Climate variability across the last deglaciation in NW Iberia and its margin. *Quat. Int.* 414, 9–22. doi: 10.1016/j.quaint.2015.08.073

North Greenland Ice-Core Project Members [NGRIP] (2004). High resolution climate record of the northern hemisphere reaching into the last glacial interglacial period. *Nature* 431, 147–151.

Oliveira, D. (2017). *Comprendre les Périodes Chaudes Pendant et Après la Transition du Pléistocène Moyen (MIS 31 et MIS 11) Dans la Péninsule Ibérique*. Bordeaux: Université de Bordeaux.

Oliveira, D., Desprat, S., Rodrigues, T., Naughton, F., Hodell, D., Trigo, R., et al. (2016). The complexity of millennial-scale variability in southwestern Europe during MIS 11. *Quat. Res.* 86, 373–387. doi: 10.1016/j.yqres.2016.09.002

Oliveira, D., Desprat, S., Yin, Q., Naughton, F., Trigo, R., Rodrigues, T., et al., (2017). Unraveling the forcings controlling the vegetation and climate of the best orbital analogues for the present interglacial in SW Europe. *Clim. Dyn.* 1, 1–20. doi: 10.1007/s00382-017-3948-7

Oppo, D. W., McManus, J. F., and Cullen, J. L. (1998). Abrupt climate events 500,000 to 340,000 Years Ago: evidence from subpolar North Atlantic sediments. *Science* 279, 1335–1338. doi: 10.1126/science.279.5355.1335

Oppo, D. W., McManus, J. F., and Cullen, J. L. (2006). Evolution and demise of the Last Interglacial warmth in the subpolar North Atlantic. *Quat. Sci. Rev.* 25, 3268–3277. doi: 10.1016/j.quascirev.2006.07.006

Peyron, O., Guiot, J., Cheddadi, R., Tarasov, P., Reille, M., de Beaulieu, J.-L., et al. (1998). Climatic reconstruction in Europe for 18,000 yr B.P. from Pollen Data. *Quat. Res.* 49, 183–196. doi: 10.1006/qres.1997.1961

Polunin, O., and Walters, M. (1985). *A Guide to the Vegetation of Britain and Europe*. New York, NY: Oxford University Press, 238.

Pons, A., and Reille, M. (1988). The Holocene-and Upper Pleistocene pollen record from Padul (Granada, Spain): a new study. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* 66, 243–263. doi: 10.1016/0031-0182(88)90202-7

Prentice, I. C., Guiot, J., Huntley, B., Jolly, D., and Cheddadi, R. (1996). Reconstructing biomes from palaeoecological data: a general method and its
application to European pollen data at 0 and 6 ka. *Clim. Dyn.* 12, 185–194. doi: 10.1007/BF00211617

Railsback, L. B., Gibbard, L., Head, M. J., Voarintsoa, N. R. G., and Toucanne, S. (2015). An optimized scheme of lettered marine isotope substages for the last 1.0 million years, and the climatostratigraphic nature of isotope stages and substages. *Quat. Sci. Rev.* 111, 94–106. doi: 10.1016/j.quascirev.2015.01.012

Rasmussen, S. O., Bigler, M., Blockley, S. P., Blunier, T., Buchardt, S. L., Clausen, H. B., et al. (2014). A stratigraphic framework for abrupt climatic changes during the Last Glacial period based on three synchronized Greenland ice-core records: refining and extending the INTIMATE event stratigraphy. *Quat. Sci. Rev.* 106, 14–28. doi: 10.1016/j.quascirev.2014.09.007

Ratmeyer, V., Balzer, W., Begametti, G., Chiapello, I., Fischer, G., and Wyputta, U. (1999). Seasonal impact of mineral dust on deep-ocean particle flux in the eastern subtropical Atlantic Ocean. *Mar. Geol.* 159, 241–252. doi: 10.1016/S0025-3227(98)00197-2

Reille, M., and de Beaulieu, J.-L. (1990). Pollen analysis of a long upper Pleistocene continental sequence in a Velay maar (Massif Central, France). *Palaeogeogr. Palaeoclimatol. Palaeoecol.* 80, 35–48. doi: 10.1016/0031-0182(90)90032-3

Reille, M., de Beaulieu, J.-L., Svobodova, V., Andrieu-Ponel, V., and Goeury, C. (2000). Pollen analytical biostratigraphy of the last five climatic cycles from a long continental sequence from the Velay region (Massif Central, France). *J. Quat. Sci.* 15, 665–685. doi: 10.1002/1099-1417(200010)15:7<665>:AID-JQS560>3.0.CO;2-G

Risebrobakken, B., Dokken, T., Ottera, O. H., Jansen, E., Gao, Y., and Drange, H. (2007). Inception of the Northern European ice sheet due to contrasting ocean (and insolation) forcing. *Quat. Res.* 67, 128–135. doi: 10.1016/j.qyres.2006.07.007

Roche, D., Paillard, D., and Cortijo, E. (2004). Constraints on the duration and freshwater release of Heinrich event 4 through isotope modelling. *Nature* 432, 379–382. doi: 10.1038/nature03059

Rossignol-Strick, M. (1983). African monsoons, an immediate climate response to orbital insolation. *Nature* 304, 46–49. doi: 10.1038/304046a0

Roucoux, K. H. (2000). *Millennial Scale Vegetation and Climate Variability in North-West Iberia during the Last Glacial Stage*. Cambridge: University of Cambridge.

Roucoux, K. H., Tzedakis, C., de Abreu, L., and Shackleton, N. J. (2006). Climate and vegetation changes 180,000 to 345,000 years ago recorded in a deep-sea core off Portugal. *Earth Planet. Sci. Lett.* 249, 307–325. doi: 10.1016/j.epsl.2006.07.005

Ruddiman, W. F. (2001). *Earth's Climate: Past and Future*. New York, NY: W.H. Freeman & Sons.

Ruddiman, W. F., and McIntyre, A. (1979). Warmth of the subpolar North Atlantic Ocean during northern hemisphere ice-sheet growth. *Science* 204, 173–175. doi: 10.1126/science.204.4389.173

Ruddiman, W. F., and McIntyre, A. (1984). Ice-age thermal response and climatic role of the surface Atlantic Ocean, 40°N to 63°N. *Geol. Soc. Am. Bull.* 95, 381–396. doi: 10.1130/0016-7606(1984)95<381:ITRACR>2.0.CO;2

Sadori, L., Koutsodendris, A., Panagiotopoulos, K., Masi, A., Bertini, A., Combourieu-Nebout, N., et al. (2016). Pollen-based paleoenvironmental and paleoclimatic change at Lake Ohrid (south-eastern Europe) during the past 500ka. *Biogeosciences* 13, 1423–1437. doi: 10.5194/bg-13-1423-2016

Sánchez Goñi, M. F. (2014). "Rapid climatic variability in Europe over the last glacial," in *Transitions, Ruptures et Continuité en Préhistoire*, ed. J. Jaubert (Paris: Bulletin de la Société Préhistorique de France), 185–193.

Sánchez Goñi, M. F. (2016). Le climat au temps des peintres de Lascaux. *Dossiers Archéol.* 376, 42–47.

Sánchez Goñi, M. F., Bakker, P., Desprat, S., Carlson, A., Van Meerbeeck, C. J., Peyron, O., et al. (2012). European climatic optimum and enhanced Greenland ice sheet melting during the last interglacial. *Geology* 40, 627–630. doi: 10.1130/G32908.1

Sánchez Goñi, M. F., Bard, E., Landais, A., Rossignol, L., and d'Errico, F. (2013). Air-sea temperature decoupling in Western Europe during the last interglacial-glacial transition. *Nat. Geosci.* 6, 837–841. doi: 10.1038/ngeo1924

Sánchez Goñi, M. F., Cacho, I., Turon, J.-L., Guiot, J., Sierro, F. J., Peypoquet, J.-P., et al. (2002). Synchroneity between marine and terrestrial responses to millennial escale climatic variability during the last glacial period in the Mediterranean region. *Clim. Dyn.* 19, 95–105. doi: 10.1007/s00382-001-0212-x

Sánchez Goñi, M. F., Desprat, S., Daniau, A.-L., Bassinot, F. C., Polanco-Martinez, J. M., and Harrison, S. P. (2017). The ACER pollen and charcoal database: a

global resource to document vegetation and fire response to abrupt climate changes during the last glacial period. *Earth Syst. Sci. Data* 9, 679–695. doi: 10.5194/essd-9-679-2017

Sánchez Goñi, M. F., Eynaud, F., Turon, J.-L., and Shackleton, N. J. (1999). High resolution palynological record off the Iberian margin: direct land-sea correlation for the Last Interglacial complex. *Earth Planet. Sci. Lett.* 171, 123–137. doi: 10.1016/S0012-821X(99)00141-7

Sánchez Goñi, M. F., and Hannon, G. (1999). High altitude vegetational patterns on the Iberian Mountain chain (north-central Spain) during the Holocene. *Holocene* 9, 39–57. doi: 10.1191/095968399671320625

Sánchez Goñi, M. F., and Harrison, S. (2010). Millennial-scale climate variability and vegetation changes during the Last Glacial: concepts and terminology. *Quat. Sci. Rev.* 29, 2823–2827. doi: 10.1111/j.1365-2486.2012.02766.x

Sánchez Goñi, M. F., Landais, A., Cacho, I., Duprat, J., and Rossignol, L. (2009). Contrasting intranterstadial climatic evolution between high and middle latitudes: a close-up of Greenland Interstadials 8 and 12. *Geochem. Geophys. Geosyst.* 10:Q04U04.

Sánchez Goñi, M. F., Landais, A., Fletcher, W. J., Naughton, F., Desprat, S., and Duprat, J. (2008). Contrasting impacts of Dansgaard-Oeschger events over a western European latitudinal transect modulated by orbital parameters. *Quat. Sci. Rev.* 27, 1136–1151. doi: 10.1016/j.quascirev.2008.03.003

Sánchez Goñi, M. F., Loutre, M. F., Peyron, O., Santos, L., Duprat, J., Malaizé, B., et al. (2005). Increasing vegetation and climate gradient in Western Europe over the Last Glacial Inception (122-110 ka): data-model comparison. *Earth Planet. Sci. Lett.* 231, 111–130. doi: 10.1016/j.epsl.2004.12.010

Sánchez Goñi, M. F., Ortu, E., Banks, W. E., Giraudeau, J., Leroyer, C., and Hanquiez, V. (2016a). The expansion of Central and Northern European Neolithic populations was associated with a multi-century warm winter and wetter climate. *Holocene* 26, 1188–1199. doi: 10.1177/0959683616638435

Sánchez Goñi, M. F., Rodrigues, T., Hodell, D. A., Polanco-Martinez, J. M., Alonso-García, M., Hernández-Almeida, I., et al. (2016b). Tropically-driven climate shifts in southwestern Europe during MIS 19, a low eccentricity interglacial. *Earth Planet. Sci. Lett.* 448, 81–93. doi: 10.1016/j.epsl.2016.05.018

Sánchez Goñi, M. F., Turon, J.-L., Eynaud, F., and Gendreau, S. (2000). European climatic response to millennial-scale climatic changes in the atmosphere-ocean system during the Last Glacial period. *Quat. Res.* 54, 394–403. doi: 10.1006/qres.2000.2176

Seillés, B., Sánchez Goñi, M. F., Ledru, M.-P., Urrego, D. H., Martinez, P., Hanquiez, V., et al. (2016). Holocene land-sea climatic links on the equatorial Pacific coast (Bay of Guayaquil, Ecuador). *Holocene* 26, 567–577. doi: 10.1177/0959683615612566

Shackleton, N. J., Hall, M. A., and Vincent, E. (2000). Phase relationships between millennial-scale events 64,000-24,000 years ago. *Paleoceanography* 15, 565–569. doi: 10.1029/2000PA000513

Shackleton, N. J., and Opdyke, N. D. (1973). Oxygen isotope and palaeomagnetic stratigraphy of Equatorial Pacific core V28-238: oxygen isotope temperatures and ice volumes on a 10⁵ year and 10⁶ year scale. *Quat. Res.* 3, 39–55. doi: 10.1016/0033-5894(73)90052-5

Shackleton, N. J., Sánchez Goñi, M. F., Pailler, D., and Lancelot, Y. (2003). Marine isotope substage 5e and the eemian interglacial. *Glob. Planet. Change* 757, 151–155. doi: 10.1016/S0921-8181(02)00181-9

Siddall, M., Rohling, E. J., Almogi-Labin, A., Hemleben, C., Meischner, D., Schmelzer, I., et al. (2003). Sea-level fluctuations during the last glacial cycle. *Nature* 423, 853–858. doi: 10.1038/nature01690

Siegenthaler, U., Stocker, T. F., Monnin, E., Lüthi, D., Schwander, J., Stauffer, B., et al. (2005). Stable carbon cycle-climate relationship during the late pleistocene. *Science* 310, 1313–1317. doi: 10.1126/science.1120130

Skinner, L. C., and Shackleton, N. J. (2006). Deconstructing Terminations I and II: revisiting the glacioeustatic paradigm based on deep-water temperature estimates. *Quat. Sci. Rev.* 25, 3312–3321. doi: 10.1016/j.quascirev.2006.07.005

Spahni, R., Chappellaz, J., Stocker, T. F., Loulergue, L., Hausmann, G., Kawamura, K., et al. (2005). Atmospheric methane and nitrous oxide of the Late Pleistocene from Antarctic ice cores. *Science* 310, 1317–1321. doi: 10.1126/science.1120132

Stanley, E. (1966). The application of palynology to oceanology with reference to the northwestern Atlantic. *Deep Sea Res. II* 13, 921–939. doi: 10.1016/0011-7471(76)90911-6
Sugita, S. (2007). Theory of quantitative reconstruction of vegetation I: pollen from large sites REVEALS regional vegetation composition. *Holocene* 17, 229–241. doi: 10.1177/0959683607075837

Traverse, A., and Ginsburg, R. N. (1966). Palynology of the surface sediments of Great Bahama Bank, as related to water movement and sedimentation. *Mar. Geol.* 4, 417–459. doi: 10.1016/0025-3227(66)90010-7

Turon, J.-L. (1984a). Direct land/sea correlations in the last interglacial complex. *Nature* 309, 673–676. doi: 10.1038/309673a0

Turon, J.-L. (1984b). *Le Palynoplancton dans l'Environnement Actuel de l'Atlantique Nord-oriental: Evolution Climatique et Hydrologique Depuis le Dernier Maximum Glaciaire*. Bordeaux: Université de Bordeaux I, 313.

Turon, J.-L., Lézine, A.-M., and Denéfle, M. (2003). Land-sea correlations for the last glaciation inferred from a pollen and dinocyst record from the Portuguese margin. *Quat. Res.* 59, 88–96. doi: 10.1016/S0033-5894(02)00018-2

Tzedakis, C., Hooghiemstra, H., and Pálike, H. (2006). The last 1.35 million years at Tenaghi Philippon: revised chronostratigraphy and long-term vegetation trends. *Quat. Sci. Rev.* 25, 3416–3430. doi: 10.1016/j.quascirev.2006.09.002

Urrego, D. H., Sánchez Goñi, M. F., Daniau, A. L., Lechevrel, S., and Hanquiez, V. (2015). Increased aridity in southwestern Africa during the warmest periods of the last interglacial. *Clim. Past* 11, 1417–1431. doi: 10.5194/cp-11-1417-2015

Van Campo, M. M. (1984). Relations entre la végétation de l'Europe et les températures de surface océaniques après le dernier maximum glaciaire. *Pollen Spores* XXVI, 497–518.

Voelker, A. H. L., and Workshop Participant (2002). Global distribution of centennial-scale records for Marine Isotope Stage (MIS) 3: a database. *Quat. Sci. Rev.* 21, 1185–1212. doi: 10.1016/S0277-3791(01)00139-1

Waelbroeck, C., Labeyrie, L., Michel, E., Duplessy, J. C., McManus, J. F., Lambeck, K., et al. (2002). Sea-level and deep water temperature changes derived from benthic foraminifer isotopic records. *Quat. Sci. Rev.* 21, 295–305. doi: 10.1016/S0277-3791(01)00101-9

Wang, X., Auler, A. S., Edwards, R. L., Cheng, H., Cristalli, S., Smart, L., et al. (2004). Wet periods in northeastern Brazil over the past 210 kyr linked to distant climate anomalies. *Nature* 432, 740–743. doi: 10.1038/nature03067

Wang, Y. J., Cheng, H., Edwards, R. L., Kong, X., Shao, X., Chen, S., et al. (2008). Millennial- and orbital-scale changes in the East Asian Monsoon over the past 224,000 years. *Nature* 451, 1090–1093. doi: 10.1038/nature06692

Webb, T. III (1986). Is vegetation in equilibrium with climate? How to interpret late-Quaternary pollen data. *Vegetatio* 67, 75–91. doi: 10.1007/BF00037359

Weijers, J. W. H., Schefuß, E., Schouten, S., and Damsté, J. S. S. (2007). Coupled thermal and hydrological evolution of tropical Africa over the last deglaciation. *Science* 315, 1701–1704. doi: 10.1126/science.1138131

Weldeab, S., Lea, D. W., Schneider, R. R., and Andersen, N. (2007). 155,000 years of West African monsoon and ocean thermal evolution. *Science* 316, 1303–1307. doi: 10.1126/science.1140461

Wolff, E. W., Chappellaz, J., Blunier, T., Rasmussen, S. O., and Svensson, A. C. (2010). Millennial-scale variability during the last glacial: the ice core record. *Quat. Sci. Rev.* 29, 2828–2838. doi: 10.1016/j.quascirev.2009.10.013

Yin, Q., and Berger, A. (2015). Interglacial analogues of the Holocene and its natural near future. *Quat. Sci. Rev.* 120, 28–46. doi: 10.1016/j.quascirev.2015.04.008

Yin, Q. Z., and Berger, A. (2012). Individual contribution of insolation and CO₂ to the interglacial climates of the past 800,000 years. *Clim. Dyn.* 38, 709–734. doi: 10.1007/s00382-011-1013-5

Zagwijn, W. H., and Veenstra, H. J. (1966). A pollen-analytical study of cores from the Outer Silver Pit, North Sea. *Mar. Geol.* 4, 539–551. doi: 10.1016/0025-3227(66)90014-4

Zorzi, C., Sánchez Goñi, M. F., Anupama, K., Prasad, S., Hanquiez, V., Johnson, J., et al. (2015). Indian monsoon variations during three contrasting climatic periods: the Holocene, HeinrichStadial 2 and the last interglacial-glacial transition. *Quat. Sci. Rev.* 125, 50–60. doi: 10.1016/j.quascirev.2015.06.009

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2018 Sánchez Goñi, Desprat, Fletcher, Morales-Molino, Naughton, Oliveira, Urrego and Zorzi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.