Abstract
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.
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 (14C, 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 -δ13C- 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 -δ18O- 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, δ13C from leaf waxes) that give useful information on precipitation and, in some cases, on temperatures (e.g., Weijers et al., 2007).
FIGURE 1
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 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.
FIGURE 2
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 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-3) 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 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).
FIGURE 3
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 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.
FIGURE 4
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 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 (δ18O) 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 (Railsback 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).
FIGURE 5
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 δ18O 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 sub-section 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).
FIGURE 6
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).
FIGURE 7
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 8
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 (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.
FIGURE 9
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, δ13C 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 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).
FIGURE 10
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 17O-excess, increase in CO2 and methane mixing ratio and heavier δD-CH4 and δ18Oatm 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 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 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).
FIGURE 11
FIGURE 12
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 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 δ13C 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.
FIGURE 13
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 δ18O 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 δ18O 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., δ18O 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.
Statements
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.
Conflict of interest
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.
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
TABLE S1List of the 129 deep-sea pollen records distributed over the world.
References
1
AllenJ. R. M.BrandtU.BrauerA.HubbertenH.-W.HuntleyB.KellerJ.et al (1999). Rapid environmental changes in southern Europe during the last glacial period.Nature400740–743. 10.1038/23432
2
Alonso-GarciaM.SierroF. J.KuceraM.FloresJ. A.CachoI.AndersenN. (2011). Ocean circulation, ice sheet growth and interhemispheric coupling of millennial climate variability during the mid-Pleistocene (ca 800–400 ka).Quat. Sci. Rev.303234–3247. 10.1016/j.quascirev.2011.08.005
3
Alvarez-SolasJ.RamsteinG. (2011). On the triggering mechanism of Heinrich events.Proc. Natl. Acad. Sci. U.S.A.108E1359–E1360. 10.1073/pnas.1116575108
4
AugustinL.BarbanteC.BarnesP. R.BarnolaJ. M.BiglerM.CastellanoE. (2004). Eight glacial cycles from an Antarctic ice core.Nature429623–628. 10.1038/nature02599
5
BaileyR. G. (1998). Ecoregions: The Ecosystem Geography of the Oceans and Continents.New York, NY: Springer.
6
BarkerS.ChenJ.GongX.JonkersL.KnorrG.ThornalleyD. (2015). Icebergs not the trigger for North Atlantic cold events.Nature520333–336. 10.1038/nature14330
7
BarkerS.KnorrG.EdwardsR. L.ParreninF.PutnamA. E.SkinnerL. C.et al (2011). 800,000 years of abrupt climate variability.Science334347–351. 10.1126/science.1203580
8
Bar-OrR.ErlickC.GildorH. (2008). The role of dust in glacial–interglacial cycles.Quat. Sci. Rev.27201–208. 10.1016/j.quascirev.2007.10.015
9
BazinL.LandaisA.Lemieux-DudonB.Toyé Mahamadou KeleH.VeresD.ParreninF.et al (2013). An optimized multi-proxy, multi-site Antarctic ice and gas orbital chronology (AICC2012): 120-800 ka.Clim. Past91715–1731. 10.5194/cp-9-1715-2013
10
BeaudouinC.SucJ.-P.EscarguelG.ArnaudM.CharmassonS. (2007). The significance of pollen signal in present-day marine terrigenous sediments: the example of the Gulf of Lions (Western Mediterranean Sea).Geobios40159–172. 10.1016/j.geobios.2006.04.003
11
BennettK. D.WillisK. J. (1995). The role of ecological factors in controlling vegetation dynamics on long temporal scales.G. Bot. Ital.129243–254. 10.1080/11263509509436132
12
BergerA. (1978). Long-term variations of daily insolation and Quaternary climatic changes.J. Atmos. Sci.352362–2367. 10.1175/1520-0469(1978)035<2362:LTVODI>2.0.CO;2
13
BergerA.LoutreM. F. (1991). Insolation values for the climate of the last 10 million years.Quat. Sci. Rev.10297–317. 10.1016/0277-3791(91)90033-Q
14
BergerA.LoutreM. F. (2004). Astronomical theory of climate change.J. Phys. IV France1211–35. 10.1051/jp4:2004121001
15
BergerA.LoutreM. F.MéliceJ. L. (2006). Equatorial insolation: from precession harmonics to eccentricity frequencies.Clim. Past2131–136. 10.5194/cp-2-131-2006
16
BirksH. J. B. (1981). “The use of pollen analysis in the reconstruction of past climates: a review,” inClimate and HistoryedsWigleyT. M. L.IngramM. J.FarmerG. (Cambridge: Cambridge University Press) 111–138.
17
Blanco CastroE.Casado GonzálezM. A.Costa TenorioM.Escribano BombínR.García AntónM.Génova FusterM.et al (1997). Los Bosques Ibéricos.Barcelona: Planeta572.
18
BondG.BroeckerW.JohnsenS.McManusJ.LabeyrieL.JouzelJ.et al (1993). Correlations between climate records from North Atlantic sediments and Greenland ice.Nature365143–147. 10.1038/365143a0
19
BondG.LottiR. (1995). Iceberg discharges into the North Atlantic on millennial time scales during the Last Glaciation.Science2671005–1009. 10.1126/science.267.5200.1005
20
BordonA.PeyronO.LézineA.-M.BrewerS.FouacheE. (2009). Pollen-inferred Late-Glacial and Holocene climate in southern Balkans (Lake Maliq).Quat. Int.20019–30. 10.1016/j.quaint.2008.05.014
21
CambonG.SucJ.-P.AloisiJ.-C.GiresseP.MonacoA.TouzaniA.et al (1997). Modern pollen deposition in the Rhone delta area (Lagoonal and marine sediments), France.Grana36105–113. 10.1080/00173139709362596
22
CarréM.CheddadiR. (2017). Seasonality in long-term climate change.Quaternaire28173–177. 10.4000/quaternaire.8018
23
ChabaudL.Sánchez GoñiM. F.DespratS.RossignolL. (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.Holocene24787–797. 10.1177/0959683614530439
24
ChapmanM. R.ShackletonN. J. (1999). Global ice-volume fluctuations, North Atlantic ice-rafted events, and deep-ocean circulation changes between 130 and 70 ka.Geology27795–798. 10.1130/0091-7613(1999)027<0795:GIVFNA>2.3.CO;2
25
ClaussenM. (2007). “Introduction to climate forcing and climate feedbacks,” inThe Climate of Past Interglacials: Developments in Quaternary ScienceVol. 7edsSirockoF.ClaussenM.Sánchez GoñiM. F.LittT. (Amsterdam: Elsevier) 3–11. 10.1016/S1571-0866(07)80026-1
26
Combourieu NeboutN.PeyronO.DormoyI.DespratS.BeaudoinC.KotthoffU.et al (2009). Rapid climatic variability in the west Mediterranean during the last 25 000 years from high resolution pollen data.Clim. Past5503–521. 10.5194/cp-5-503-2009
27
CrucifixM.LoutreM. F. (2002). Transient simulations over the last interglacial period (126-115 kyr BP): feedback and forcing analysis.Clim. Dyn.19417–433. 10.1007/s00382-002-0234-z
28
CrucifixM.MitsuiT.LenoirG. (2016). “Challenges for ice age dynamics: a dynamical systems perspective,” inNonlinear and Stochastic Climate DynamicsedsFranzkeC.O’KaneT. J. (Cambridge: Cambridge University Press) 1–32.
29
DaniauA. L.Sánchez GoñiM.-F.BeaufortL.Laggoun-DéfargeF.LoutreM.-F.DupratJ. (2007). Dansgaard-Oeschger climatic variability revealed by fire emissions in southwestern Iberia.Quat. Sci. Rev.261369–1383. 10.1016/j.quascirev.2007.02.005
30
DaniauA.-L.Sánchez GoñiM. F.MartinezP.UrregoD. H.Bout-RoumazeillesV.DespratS.et al (2013). Orbital-scale climate forcing of grassland burning in southern Africa.Proc. Natl. Acad. Sci. U.S.A.1105069–5073. 10.1073/pnas.1214292110
31
DansgaardW.JohnsenS.ClausenH. B.Dahl-JensenD.GundestrupN.HammerC. U.et al (1984). “North Atlantic climatic oscillations revealed by deep Greenland ice cores,” inClimate Processes and Climate SensitivityedsHansenJ. E.TakahashiT. (Washington, DC: American Geophysical Union) 288–298.
32
de BeaulieuJ.-L.Andrieu-PonelV.ReilleM.GrügerE.TzedakisC.SvobodovaH. (2001). An attempt at correlation between the Velay pollen sequence and the Middle Pleistocene stratigraphy from central Europe.Quat. Sci. Rev.201593–1602. 10.1016/S0277-3791(01)00027-0
33
deMenocalB. (2004). African climate change and faunal evolution during the Plio-Pleistocene.Earth Planet. Sci. Lett.20043–24. 10.1016/S0012-821X(04)00003-2
34
DespratS.Diaz FernandezM.CoulonT.EzzatL.Pessarossi-LangloisJ.GilL.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.421998–2009. 10.1111/jbi.12566
35
DespratS.OliveiraO.NaughtonF.Sánchez GoñiM. 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.Quaternaire28259–269. 10.4000/quaternaire.8102
36
DespratS.Sánchez GoñiM. F.McManusJ. F.DupratJ.CortijoE. (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. Past553–72. 10.5194/cp-5-53-2009
37
DespratS.Sánchez GoñiM. F.NaughtonF.TuronJ. L.DupratJ.MalaizéB.et al (2007). “Climate variability of the last five isotopic interglacials from direct land-sea-ice correlation,” inThe Climate of Past InterglacialsedsSirockoF.ClaussenM.Sánchez GoñiM. F.LittT. (Amsterdam: Elsevier) 375–386.
38
DespratS.Sánchez GoñiM. F.TuronJ.-L.DupratJ.MalaizéB.PeypouquetJ.-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.251010–1026. 10.1016/j.quascirev.2006.01.001
39
DespratS.Sánchez GoñiM. F.TuronJ.-L.McManusJ. F.LoutreM. F.DupratJ.et al (2005). Is vegetation responsible for glacial inception during periods of muted insolation changes?Quat. Sci. Rev.241361–1374. 10.1016/j.quascirev.2005.01.005
40
DupontL. (2011). Orbital scale vegetation change in Africa.Quat. Sci. Rev.303589–3602. 10.1016/j.quascirev.2011.09.019
41
DupontL. M.WyputtaU. (2003). Reconstructing pathways of aeolian pollen transport to the marine sediments along the coastline of SW Africa.Quat. Sci. Rev.22157–174. 10.1016/S0277-3791(02)00032-X
42
ElliotM.LabeyrieL.DokkenT.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.194151–163. 10.1016/S0012-821X(01)00561-1
43
FerrettiP.CrowhurstS. J.NaafsB. D. A.BarbanteC. (2015). The marine isotope stage 19 in the mid-latitude North Atlantic Ocean: astronomical signature and intra-interglacial variability.Quat. Sci. Rev.10895–110. 10.1016/j.quascirev.2014.10.024
44
FletcherW. J.Sánchez GoñiM. F. (2008). Orbital- and sub-orbital-scale climate impacts on vegetation of the western Mediterranean basin over the last 48,000 yr.Quat. Res.70451–464. 10.1016/j.yqres.2008.07.002
45
FletcherW. J.Sánchez GoñiM. F.AllenJ. R. M.CheddadiR.Combourieu-NeboutN.HuntleyB.et al (2010). Millennial-scale variability during the last glacial in vegetation records from Europe.Quat. Sci. Rev.292839–2864. 10.1016/j.quascirev.2009.11.015
46
FollieriM.MagriD.SadoriL. (1988). 250.000-year pollen record from valle di Castiglione (Roma).Pollen Spores30329–356.
47
GanopolskiA.RahmstorfS. (2001). Rapid changes of glacial climate simulated in a coupled climate model.Nature409153–158. 10.1038/35051500
48
GanopolskiA.WinkelmannR.SchellnhuberH. J. (2016). Critical insolation–CO2 relation for diagnosing past and future glacial inception.Nature529200–203. 10.1038/nature16494
49
GouveiaC.TrigoR. M.DaCamaraC. C.LibonatiR.PereiraJ. M. C. (2008). The North Atlantic Oscillation and European vegetation dynamics.Int. J. Climatol.281835–1847. 10.1002/joc.1682
50
GrootJ. J.GrootC. R. (1966). Marine palynology: possibilities, limitations, problems.Mar. Geol.4387–395. 10.1016/0025-3227(66)90007-7
51
GuillevicM.BazinL.LandaisA.StowasserC.Masson-DelmotteV.BlunierT.et al (2014). Evidence for a three-phase sequence during Heinrich Stadial 4 using a multiproxy approach based on Greenland ice core records.Clim. Past102115–2133. 10.5194/cp-10-2115-2014
52
GuiotJ. (1990). Methodology of the last climatic cycle reconstruction from pollen data.Palaeogeogr. Palaeoclimatol. Palaeoecol.8049–69. 10.1016/0031-0182(90)90033-4
53
GuiotJ.De VernalA. (2007). “Transfer functions: methods for quantitative paleoceanography based on microfossils,” inDevelopments in Marine GeologyVol. 1edsHillaire-MarcelC.de VernalA. (Amsterdam: Elsevier) 523–563.
54
HaysJ. D.ImbrieJ.ShackeltonN. J. (1976). Variations in the Earth’s orbit pacemaker of the ice ages.Science1941121–1132. 10.1126/science.194.4270.1121
55
HeinrichH. (1988). Origin and consequences of cyclic ice rafting in the northeast Atlantic Ocean during the past 130,000 years.Quat. Res.29142–152. 10.1016/0033-5894(88)90057-9
56
HeusserL. E.BalsamW. L. (1977). Pollen distribution in the N.E. Pacific Ocean.Quat. Res.745–62. 10.1016/0033-5894(77)90013-8
57
HeusserL. E.Van de GeerG. (1994). Direct correlation of terrestrial and marine paleoclimatic records from four glacial-interglacial cycles-DSDP site 594 Southwest Pacific.Quat. Sci. Rev.13273–282. 10.1016/0277-3791(94)90030-2
58
HeusserL. E.VmBryantJr. (1985). “Quaternary palynology of marine sediments in the northeast Pacific, northwest Atlantic, and Gulf of Mexico,” inPollen Records of Late-Quaternary North American Sedimentsed.HollowayR. G. (Dallas, TX: AASP Foundation) 385–403.
59
HeusserL. E. (1986/1987). Pollen in marine cores: evidence of past climates.Oceanus2964–70.
60
HodellD. A. (2016). The smoking gun of the ice ages.Science3541235–1236. 10.1126/science.aal4111
61
HodellD. A.LourensL.CrowhurstS.KonijnendijkT.TjallingiiR.Jiménez-EspejoF.et al (2015). A reference time scale for site U1385 (Shackleton Site) on the SW Iberian Margin.Glob. Planet. Change13349–64. 10.1016/j.gloplacha.2015.07.002
62
HoerlingM.EischeidJ.PerlwitzJ.QuanX. W.ZhangT.PegionP. (2012). On the increased frequency of Mediterranean drought.J. Clim.252146–2161. 10.1175/JCLI-D-11-00296.1
63
HooghiemstraH.LézineA.-M.LeroyS. A. G.DupontL.MarretF. (1988). Late Quaternary palynology in marine sediments: a synthesis of the understanding of pollen distribution patterns in the NW African setting.Quat. Int.14829–44. 10.1016/j.quaint.2005.11.005
64
HooghiemstraH.StallingH.AgwuC. O. C.DupontL. 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.741–53. 10.1016/0034-6667(92)90137-6
65
HooghmiestraH.AgwuC. O. C.BeugH.-J. (1986). Pollen and spore distribution in recent marine sediments: a record of NW-African seasonal wind patterns and vegetation belts.Meteor Forsch. Ergeb. C4087–135.
66
JullienE.GroussetF. E.MalaizeB.DupratJ.Sánchez GoniM. F.EynaudF.et al (2007). Tropical “dusty Heinrich-like events” vs. high-latitude “icy Heinrich events”.Quat. Res.68379–386. 10.1016/j.yqres.2007.07.007
67
KageyamaM.MignotJ.SwingedouwD.MarzinC.AlkamaR.MartiO. (2009). Glacial climate sensitivity to different states of the Atlantic Meridional Overturning Circulation: results from the IPSL model.Clim. Past5551–570. 10.5194/cp-5-551-2009
68
KnuttiR.FlückigerJ.StockerT. F.TimmermannA. (2004). Strong hemispheric coupling of glacial climate through freshwater discharge and ocean circulation.Nature430851–856. 10.1038/nature02786
69
KorenevaE. V. (1966). Marine palynological researches in the U.S.S.R.Mar. Geol.4565–574. 10.1016/0025-3227(66)90016-8
70
KorenevaE. V. (1971). “Spores and pollen in Mediterranean bottom sediments,” inThe Micropaleontology of OceansedsFunnelB. M.RiedelW. R. (Cambridge: Cambridge University Press) 361–371.
71
KühlN.GebhardtC.LittT.HenseA. (2002). Probability density functions as botanical-climatological transfer functions for climate reconstruction.Quat. Res.58381–392. 10.1006/qres.2002.2380
72
LandaisA.BarnolaJ. M.Masson-DelmotteV.JouzelJ.ChappellazJ.CaillonN.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.31101–113. 10.1029/2004GL021193
73
LevavasseurG.VracM.RocheD. M.PaillardD. (2012). Statistical modelling of a new global potential vegetation distribution.Environ. Res. Lett.7:044019. 10.1088/1748-9326/7/4/044019
74
LoulergueL.SchiltA.SpahniR.Masson-DelmotteV.BlunierT.LemieuxB.et al (2008). Orbital and millennial-scale features of atmospheric CH4 over the past 800,000 years.Nature453383–386. 10.1038/nature06950
75
LuthiD.Le FlochM.BereiterB.BlunierT.BarnolaJ.-M.SiegenthalerU.et al (2008). High-resolution carbon dioxide concentration record 650,000-800,000 years before present.Nature453379–382. 10.1038/nature06949
76
Lynch-StieglitzJ. (2017). The Atlantic Meridional Overturning Circulation and abrupt climate change.Annu. Rev. Mar. Sci.983–104. 10.1146/annurev-marine-010816-060415
77
MartratB.GrimaltJ. O.ShackletonN. J.De AbreuL.HutterliM. A.StockerT. F. (2007). Four cycles of recurring deep and surface water destabilizations on the Iberian margin.Science317502–507. 10.1126/science.1139994
78
McIntyreA.MolfinoB. (1996). Forcing of Atlantic Equatorial and subpolar millennial cycles by precession.Science2741867–1870. 10.1126/science.274.5294.1867
79
McManusJ. F.BondG. C.BroeckerW. S.JohnsenS.LabeyrieL.HigginsS. (1994). High-resolution climate records from the North Atlantic during the last interglacial.Nature371326–329. 10.1038/371326a0
80
McManusJ. F.FrancoisR.GherardiJ.-M.KeigwinL. D.Brown-LegerS. (2004). Collapse and rapid resumption of Atlantic meridional circulation linked to deglacial climate changes.Nature428834–837. 10.1038/nature02494
81
McManusJ. F.OppoD. W.CullenJ. L. (1999). A 0.5-million-year record of millennial-scale climate variability in the North Atlantic.Science283971–975. 10.1126/science.283.5404.971
82
MilankovitchM. M. (1941). Kanon der Erdbestrahlung.Beograd: Köninglich Serbische Akademie.
83
MontadeV.Combourieu-NeboutN.KisselC.HaberleS. H.SianiG.MichelE. (2013). Vegetation and climate changes during the last 22000 yr from a marine core near Taitao Peninsula, southern Chile.Palaeogeogr. Palaeoclimatol. Palaeoecol.369335–348. 10.1016/j.palaeo.2012.11.001
84
MudieJ.McCarthyF. M. G. (2006). Marine palynology: potentials for onshore–offshore correlation of Pleistocene–Holocene records.Trans. R. Soc. S. Afr.61139–157. 10.1080/00359190609519964
85
NaughtonF.Sánchez GoñiM. F.DespratS.TuronJ.-L.DupratJ.MalaizéB.et al (2007). Present-day and past (last 25000 years) marine pollen signal off western Iberia.Mar. Micropaleontol.6291–114. 10.1016/j.marmicro.2006.07.006
86
NaughtonF.Sánchez GoñiM. F.KageyamaM.BardE.DupratJ.CortijoE.et al (2009). Wet to dry climatic trend in north western Iberia within Heinrich events.Earth Planet. Sci. Lett.284329–342. 10.1016/j.epsl.2009.05.001
87
NaughtonF.Sánchez GoñiM. F.RodriguesT.SalgueiroE.CostasS.DespratS.et al (2016). Climate variability across the last deglaciation in NW Iberia and its margin.Quat. Int.4149–22. 10.1016/j.quaint.2015.08.073
88
North Greenland Ice-Core Project Members [NGRIP] (2004). High resolution climate record of the northern hemisphere reaching into the last glacial interglacial period.Nature431147–151.
89
OliveiraD. (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.
90
OliveiraD.DespratS.RodriguesT.NaughtonF.HodellD.TrigoR.et al (2016). The complexity of millennial-scale variability in southwestern Europe during MIS 11.Quat. Res.86373–387. 10.1016/j.yqres.2016.09.002
91
OliveiraD.DespratS.YinQ.NaughtonF.TrigoR.RodriguesT.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.11–20. 10.1007/s00382-017-3948-7
92
OppoD. W.McManusJ. F.CullenJ. L. (1998). Abrupt climate events 500000 to 340000 Years Ago: evidence from subpolar North Atlantic sediments.Science2791335–1338. 10.1126/science.279.5355.1335
93
OppoD. W.McManusJ. F.CullenJ. L. (2006). Evolution and demise of the Last Interglacial warmth in the subpolar North Atlantic.Quat. Sci. Rev.253268–3277. 10.1016/j.quascirev.2006.07.006
94
PeyronO.GuiotJ.CheddadiR.TarasovP.ReilleM.de BeaulieuJ.-L.et al (1998). Climatic reconstruction in Europe for 18,000 yr B.P. from Pollen Data.Quat. Res.49183–196. 10.1006/qres.1997.1961
95
PoluninO.WaltersM. (1985). A Guide to the Vegetation of Britain and Europe.New York, NY: Oxford University Press238.
96
PonsA.ReilleM. (1988). The Holocene-and Upper Pleistocene pollen record from Padul (Granada, Spain): a new study.Palaeogeogr. Palaeoclimatol. Palaeoecol.66243–263. 10.1016/0031-0182(88)90202-7
97
PrenticeI. C.GuiotJ.HuntleyB.JollyD.CheddadiR. (1996). Reconstructing biomes from palaeoecological data: a general method and its application to European pollen data at 0 and 6 ka.Clim. Dyn.12185–194. 10.1007/BF00211617
98
RailsbackL. B.GibbardL.HeadM. J.VoarintsoaN. R. G.ToucanneS. (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.11194–106. 10.1016/j.quascirev.2015.01.012
99
RasmussenS. O.BiglerM.BlockleyS. P.BlunierT.BuchardtS. L.ClausenH. 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.10614–28. 10.1016/j.quascirev.2014.09.007
100
RatmeyerV.BalzerW.BegamettiG.ChiapelloI.FischerG.WyputtaU. (1999). Seasonal impact of mineral dust on deep-ocean particle flux in the eastern subtropical Atlantic Ocean.Mar. Geol.159241–252. 10.1016/S0025-3227(98)00197-2
101
ReilleM.de BeaulieuJ.-L. (1990). Pollen analysis of a long upper Pleistocene continental sequence in a Velay maar (Massif Central, France).Palaeogeogr. Palaeoclimatol. Palaeoecol.8035–48. 10.1016/0031-0182(90)90032-3
102
ReilleM.de BeaulieuJ.-L.SvobodovaV.Andrieu-PonelV.GoeuryC. (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.15665–685. 10.1002/1099-1417(200010)15:7<665::AID-JQS560>3.0.CO;2-G
103
RisebrobakkenB.DokkenT.OtteraO. H.JansenE.GaoY.DrangeH. (2007). Inception of the Northern European ice sheet due to contrasting ocean and insolation forcing.Quat. Res.67128–135. 10.1016/j.yqres.2006.07.007
104
RocheD.PaillardD.CortijoE. (2004). Constraints on the duration and freshwater release of Heinrich event 4 through isotope modelling.Nature432379–382. 10.1038/nature03059
105
Rossignol-StrickM. (1983). African monsoons, an immediate climate response to orbital insolation.Nature30446–49. 10.1038/304046a0
106
RoucouxK. H. (2000). Millennial Scale Vegetation and Climate Variability in North-West Iberia during the Last Glacial Stage.Cambridge: University of Cambridge.
107
RoucouxK. H.TzedakisC.de AbreuL.ShackletonN. J. (2006). Climate and vegetation changes 180000 to 345000 years ago recorded in a deep-sea core off Portugal.Earth Planet. Sci. Lett.249307–325. 10.1016/j.epsl.2006.07.005
108
RuddimanW. F. (2001). Earth’s Climate: Past and Future.New York, NY: W.H. Freeman & Sons.
109
RuddimanW. F.McIntyreA. (1979). Warmth of the subpolar North Atlantic Ocean during northern hemisphere ice-sheet growth.Science204173–175. 10.1126/science.204.4389.173
110
RuddimanW. F.McintyreA. (1984). Ice-age thermal response and climatic role of the surface Atlantic Ocean, 40°N to 63°N.Geol. Soc. Am. Bull.95381–396. 10.1130/0016-7606(1984)95<381:ITRACR>2.0.CO;2
111
SadoriL.KoutsodendrisA.PanagiotopoulosK.MasiA.BertiniA.Combourieu-NeboutN.et al (2016). Pollen-based paleoenvironmental and paleoclimatic change at Lake Ohrid (south-eastern Europe) during the past 500ka.Biogeosciences131423–1437. 10.5194/bg-13-1423-2016
112
Sánchez GoñiM. F. (2014). “Rapid climatic variability in Europe over the last glacial,” inTransitions, Ruptures et Continuité en Préhistoireed.JaubertJ. (Paris: Bulletin de la Société Préhistorique de France) 185–193.
113
Sánchez GoñiM. F. (2016). Le climat au temps des peintres de Lascaux.Dossiers Archéol.37642–47.
114
Sánchez GoñiM. F.BakkerP.DespratS.CarlsonA.Van MeerbeeckC. J.PeyronO.et al (2012). European climatic optimum and enhanced Greenland ice sheet melting during the last interglacial.Geology40627–630. 10.1130/G32908.1
115
Sánchez GoñiM. F.BardE.LandaisA.RossignolL.d’ErricoF. (2013). Air-sea temperature decoupling in Western Europe during the last interglacial-glacial transition.Nat. Geosci.6837–841. 10.1038/ngeo1924
116
Sánchez GoñiM. F.CachoI.TuronJ.-L.GuiotJ.SierroF. J.PeypouquetJ.-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.1995–105. 10.1007/s00382-001-0212-x
117
Sánchez GoñiM. F.DespratS.DaniauA.-L.BassinotF. C.Polanco-MartinezJ. M.HarrisonS. 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. Data9679–695. 10.5194/essd-9-679-2017
118
Sánchez GoñiM. F.EynaudF.TuronJ.-L.ShackletonN. J. (1999). High resolution palynological record off the Iberian margin: direct land-sea correlation for the Last Interglacial complex.Earth Planet. Sci. Lett.171123–137. 10.1016/S0012-821X(99)00141-7
119
Sánchez GoñiM. F.HannonG. (1999). High altitude vegetational patterns on the Iberian Mountain chain (north-central Spain) during the Holocene.Holocene939–57. 10.1191/095968399671230625
120
Sánchez GoñiM. F.HarrisonS. (2010). Millennial-scale climate variability and vegetation changes during the Last Glacial: concepts and terminology.Quat. Sci. Rev.292823–2827. 10.1111/j.1365-2486.2012.02766.x
121
Sánchez GoñiM. F.LandaisA.CachoI.DupratJ.RossignolL. (2009). Contrasting intrainterstadial climatic evolution between high and middle latitudes: a close-up of Greenland Interstadials 8 and 12.Geochem. Geophys. Geosyst.10:Q04U04.
122
Sánchez GoñiM. F.LandaisA.FletcherW. J.NaughtonF.DespratS.DupratJ. (2008). Contrasting impacts of Dansgaard-Oeschger events over a western European latitudinal transect modulated by orbital parameters.Quat. Sci. Rev.271136–1151. 10.1016/j.quascirev.2008.03.003
123
Sánchez GoñiM. F.LoutreM. F.PeyronO.SantosL.DupratJ.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.231111–130. 10.1016/j.epsl.2004.12.010
124
Sánchez GoñiM. F.OrtuE.BanksW. E.GiraudeauJ.LeroyerC.HanquiezV. (2016a). The expansion of Central and Northern European Neolithic populations was associated with a multi-century warm winter and wetter climate.Holocene261188–1199. 10.1177/0959683616638435
125
Sánchez GoñiM. F.RodriguesT.HodellD. A.Polanco-MartínezJ. M.Alonso-GarcíaM.Hernández-AlmeidaI.et al (2016b). Tropically-driven climate shifts in southwestern Europe during MIS 19 a low eccentricity interglacial.Earth Planet. Sci. Lett.44881–93. 10.1016/j.epsl.2016.05.018
126
Sánchez GoñiM. F.TuronJ.-L.EynaudF.GendreauS. (2000). European climatic response to millenial-scale climatic changes in the atmosphere-ocean system during the Last Glacial period.Quat. Res.54394–403. 10.1006/qres.2000.2176
127
SeillèsB.Sánchez GoñiM. F.LedruM.-P.UrregoD. H.MartinezP.HanquiezV.et al (2016). Holocene land–sea climatic links on the equatorial Pacific coast (Bay of Guayaquil, Ecuador).Holocene26567–577. 10.1177/0959683615612566
128
ShackletonN. J.HallM. A.VincentE. (2000). Phase relationships between millennial-scale events 64000-24000 years ago.Paleoceanography15565–569. 10.1029/2000PA000513
129
ShackletonN. J.OpdykeN. D. (1973). Oxygen isotope and palaeomagnetic stratigraphy of Equatorial Pacific core V28-238: oxygen isotope temperatures and ice volumes on a 105 year and 106 year scale.Quat. Res.339–55. 10.1016/0033-5894(73)90052-5
130
ShackletonN. J.Sánchez GoñiM. F.PaillerD.LancelotY. (2003). Marine isotope substage 5e and the eemian interglacial.Glob. Planet. Change757151–155. 10.1016/S0921-8181(02)00181-9
131
SiddallM.RohlingE. J.Almogi-LabinA.HemlebenC.MeischnerD.SchmelzerI.et al (2003). Sea-level fluctuations during the last glacial cycle.Nature423853–858. 10.1038/nature01690
132
SiegenthalerU.StockerT. F.MonninE.LüthiD.SchwanderJ.StaufferB.et al (2005). Stable carbon cycle–climate relationship during the late pleistocene.Science3101313–1317. 10.1126/science.1120130
133
SkinnerL. C.ShackletonN. J. (2006). Deconstructing Terminations I and II: revisiting the glacioeustatic paradigm based on deep-water temperature estimates.Quat. Sci. Rev.253312–3321. 10.1016/j.quascirev.2006.07.005
134
SpahniR.ChappellazJ.StockerT. F.LoulergueL.HausammannG.KawamuraK.et al (2005). Atmospheric methane and nitrous oxide of the Late Pleistocene from Antarctic ice cores.Science3101317–1321. 10.1126/science.1120132
135
StanleyE. (1966). The application of palynology to oceanology with reference to the northwestern Atlantic.Deep Sea Res. II13921–939. 10.1016/0011-7471(76)90911-6
136
SugitaS. (2007). Theory of quantitative reconstruction of vegetation I: pollen from large sites REVEALS regional vegetation composition.Holocene17229–241. 10.1177/0959683607075837
137
TraverseA.GinsburgR. N. (1966). Palynology of the surface sediments of Great Bahama Bank, as related to water movement and sedimentation.Mar. Geol.4417–459. 10.1016/0025-3227(66)90010-7
138
TuronJ.-L. (1984a). Direct land/sea correlations in the last interglacial complex.Nature309673–676. 10.1038/309673a0
139
TuronJ.-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 I313.
140
TuronJ.-L.LézineA.-M.DenèfleM. (2003). Land-sea correlations for the last glaciation inferred from a pollen and dinocyst record from the Portuguese margin.Quat. Res.5988–96. 10.1016/S0033-5894(02)00018-2
141
TzedakisC.HooghiemstraH.PälikeH. (2006). The last 1.35 million years at Tenaghi Philippon: revised chronostratigraphy and long-term vegetation trends.Quat. Sci. Rev.253416–3430. 10.1016/j.quascirev.2006.09.002
142
UrregoD. H.Sánchez GoñiM. F.DaniauA. L.LechevrelS.HanquiezV. (2015). Increased aridity in southwestern Africa during the warmest periods of the last interglacial.Clim. Past111417–1431. 10.5194/cp-11-1417-2015
143
Van CampoM. 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 SporesXXVI497–518.
144
VoelkerA. H. L.Workshop Participant (2002). Global distribution of centennial-scale records for Marine Isotope Stage (MIS) 3: a database.Quat. Sci. Rev.211185–1212. 10.1016/S0277-3791(01)00139-1
145
WaelbroeckC.LabeyrieL.MichelE.DuplessyJ. C.McManusJ. F.LambeckK.et al (2002). Sea-level and deep water temperature changes derived from benthic foraminifer isotopic records.Quat. Sci. Rev.21295–305. 10.1016/S0277-3791(01)00101-9
146
WangX.AulerA. S.EdwardsR. L.ChengH.CristalliS.SmartL.et al (2004). Wet periods in northeastern Brazil over the past 210 kyr linked to distant climate anomalies.Nature432740–743. 10.1038/nature03067
147
WangY. J.ChengH.EdwardsR. L.KongX.ShaoX.ChenS.et al (2008). Millennial- and orbital-scale changes in the East Asian Monsoon over the past 224000 years.Nature4511090–1093. 10.1038/nature06692
148
WebbT.III (1986). Is vegetation in equilibrium with climate? How to interpret late-Quaternary pollen data.Vegetatio6775–91. 10.1007/BF00037359
149
WeijersJ. W. H.SchefußE.SchoutenS.DamstéJ. S. S. (2007). Coupled thermal and hydrological evolution of tropical Africa over the last deglaciation.Science3151701–1704. 10.1126/science.1138131
150
WeldeabS.LeaD. W.SchneiderR. R.AndersenN. (2007). 155000 years of West African monsoon and ocean thermal evolution.Science3161303–1307. 10.1126/science.1140461
151
WolffE. W.ChappellazJ.BlunierT.RasmussenS. O.SvenssonA. C. (2010). Millennial-scale variability during the last glacial: the ice core record.Quat. Sci. Rev.292828–2838. 10.1016/j.quascirev.2009.10.013
152
YinQ.BergerA. (2015). Interglacial analogues of the Holocene and its natural near future.Quat. Sci. Rev.12028–46. 10.1016/j.quascirev.2015.04.008
153
YinQ. Z.BergerA. (2012). Individual contribution of insolation and CO2 to the interglacial climates of the past 800,000 years.Clim. Dyn.38709–734. 10.1007/s00382-011-1013-5
154
ZagwijnW. H.VeenstraH. J. (1966). A pollen-analytical study of cores from the Outer Silver Pit, North Sea.Mar. Geol.4539–551. 10.1016/0025-3227(66)90014-4
155
ZorziC.Sánchez GoñiM. F.AnupamaK.PrasadS.HanquiezV.JohnsonJ.et al (2015). Indian monsoon variations during three contrasting climatic periods: the Holocene, HeinrichStadial 2 and the last interglacial-glacial transition.Quat. Sci. Rev.12550–60. 10.1016/j.quascirev.2015.06.009
Summary
Keywords
vegetation, millennial-scale climate variability, Dansgaard-Oeschger cycles, Heinrich events, glaciations, interglacials, Europe, Quaternary
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
Received
12 July 2017
Accepted
09 January 2018
Published
26 January 2018
Volume
9 - 2018
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
Updates
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.
*Correspondence: María F. Sánchez Goñi, maria.sanchez-goni@u-bordeaux.fr
This article was submitted to Agroecology and Land Use Systems, a section of the journal Frontiers in Plant Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.