Millennial and sub-millennial scale climatic variations recorded in polar ice cores over the last glacial period

Since its discovery in Greenland ice cores, the millennial scale climatic variability of the last glacial period has been increasingly documented at all latitudes with studies focusing mainly on Marine Isotopic Stage 3 (MIS 3; 28–60 thousand of years before present, hereafter ka) and characterized by short Dansgaard-Oeschger (DO) events. Recent and new results obtained on the EPICA and NorthGRIP ice cores now precisely describe the rapid variations of Antarctic and Greenland temperature during MIS 5 (73.5–123 ka), a time period corresponding to relatively high sea level. The results display a succession of abrupt events associated with long Greenland InterStadial phases (GIS) enabling us to highlight a sub-millennial scale climatic variability depicted by (i) short-lived and abrupt warming events preceding some GIS (precursor-type events) and (ii) abrupt warming events at the end of some GIS (rebound-type events). The occurrence of these sub-millennial scale events is suggested to be driven by the insolation at high northern latitudes together with the internal forcing of ice sheets. Thanks to a recent NorthGRIP-EPICA Dronning Maud Land (EDML) common timescale over MIS 5, the bipolar sequence of climatic events can be established at millennial to sub-millennial timescale. This shows that for extraordinary long stadial durations the accompanying Antarctic warming amplitude cannot be described by a simple linear relationship between the two as expected from the bipolar seesaw concept. We also show that when ice sheets are extensive, Antarctica does not necessarily warm during the whole GS as the thermal bipolar seesaw model would predict, questioning the Greenland ice core temperature records as a proxy for AMOC changes throughout the glacial period.

Integrated reconstruction of Holocene millennial-scale environmental changes in Tierra del Fuego, southernmost South America

This study presents new paleoenvironmental data obtained from sedimentary cores from Lago Fagnano, an elon- gated lake located at 54°S in southernmost South America. Data from palynomorphs (pollen, spores and algae) and associated palynofacies as well as from diatom taxa retrieved from these cores compared with other regional proxies contribute to evaluate the similarities and differences in the climate patterns based on different proxies from southernmost Patagonia. The pollen analysis reveals that a grass steppe environment existed during the early Holocene (11,300–~8000 cal a BP) followed by a major vegetation change characterized by development of forest-steppe ecotone communities between ~8000 and ~6500 cal a BP, under more humid conditions. Between ~ 6500 and ~ 4000 cal a BP, expansion and colonization by Nothofagus forests reflect an increase in effec- tive moisture levels, while openness in the forest communities characterizes the region after ~ 1100 cal a BP. The palynological organic matter combined with the algal content reflects hydrological changes occurring in the lake and its nutrient status, probably in close relation with past climate oscillations. All these past ecological changes are closely related to oscillations in precipitation and temperature as a response to the variations in the latitudinal position and/or strength of the Southern Westerlies wind belt during the Holocene.

Biological and physical controls in the Southern Ocean on past millennial-scale atmospheric CO2 changes

Millennial-scale climate changes during the last glacial period and deglaciation were accompanied by rapid changes in atmospheric CO2 that remain unexplained. While the role of the Southern Ocean as a ’control valve’ on ocean–atmosphere CO2 exchange has been emphasized, the exact nature of this role, in particular the relative contributions of physical (for example, ocean dynamics and air–sea gas exchange) versus biological processes (for example, export productivity), remains poorly constrained. Here we combine reconstructions of bottom-water [O2], export production and 14C ventilation ages in the sub-Antarctic Atlantic, and show that atmospheric CO2 pulses during the last glacial- and deglacial periods were consistently accompanied by decreases in the biological export of carbon and increases in deep-ocean ventilation via southern-sourced water masses. These findings demonstrate how the Southern Ocean’s ’organic carbon pump’ has exerted a tight control on atmospheric CO2, and thus global climate, specifically via a synergy of both physical and biological processes.

Patterns of millennial variability over the last 500 ka

Millennial variability is a robust feature of many paleoclimate records, at least throughout the last several glacial cycles. Here we use the mean signal from Antarctic climate events 1 to 4 to probe the EPICA Dome C temperature proxy reconstruction through the last 500 ka for similar millennial-scale events. We find that clusters of millennial events occurred in a regular fashion over half of the time during this with a mean recurrence interval of 21 kyr. We find that there is no consistent link between ice-rafted debris deposition and millennial variability. Instead we speculate that changes in the zonality of atmospheric circulation over the North Atlantic form a viable alternative to freshwater release from icebergs as a trigger for millennial variability. We suggest that millennial changes in the zonality of atmospheric circulation over the North Atlantic are linked to precession via sea-ice feedbacks and that this relationship is modified by the presence of the large, Northern Hemisphere ice sheets during glacial periods.

Millennial changes in North American wildfire and soil activity over the last glacial cycle

Climate changes in the North Atlantic region during the last glacial cycle were dominated by the slow waxing and waning of the North American ice sheet as well as by intermittent, millennial-scale Dansgaard–Oeschger climate oscillations. However, prior to the last deglaciation, the responses of North American vegetation and biomass burning to these climate variations are uncertain. Ammonium in Greenland ice cores, a product from North American soil emissions and biomass burning events, can help to fill this gap. Here we use continuous, high-resolution measurements of ammonium concentrations between 110,000 to 10,000 years ago from the Greenland NGRIP and GRIP ice cores to reconstruct North American wildfire activity and soil ammonium emissions. We find that on orbital timescales soil emissions increased under warmer climate conditions when vegetation expanded northwards into previously ice-covered areas. For millennial-scale interstadial warm periods during Marine Isotope Stage 3, the fire recurrence rate increased in parallel to the rapid warmings, whereas soil emissions rose more slowly, reflecting slow ice shrinkage and delayed ecosystem changes. We conclude that sudden warming events had little impact on soil ammonium emissions and ammonium transport to Greenland, but did result in a substantial increase in the frequency of North American wildfires.

Impact of Holocene climate changes on alpine and treeline vegetation at Sanetsch Pass, Bernese Alps, Switzerland

In order to infer reactions of treeline and alpine vegetation to climatic change, past vegetation changes are reconstructed on the basis of pollen, macrofossil and charcoal analysis. The sampled sediment cores originate from the small pond Emines, located at the Sanetsch Pass (connecting the Valais and Bern, Switzerland) at an altitude of 2288 m a.s.l. Today's treeline is at ca. 2200 m a.s.l. in the area, though due to special pass (saddle) conditions it is locally depressed to ca. 2060 m a.s.l. Our results reveal that the area around Emines was covered by treeless alpine vegetation during most of the past 12,000 years. Single individuals of Betula, Larix decidua and possibly Pinus cembra occurred during the Holocene. Major centennial to millennial-scale responses of treeline vegetation to climatic changes are evident. However, alpine vegetation composition remained rather stable between 11,500 and 6000 cal. BP, showing that Holocene climatic changes of +/− 1 °C hardly influenced the local vegetation at Emines. The rapid warming of 3–4 °C at the Late Glacial/Holocene transition (11,600 cal. BP) caused significant altitudinal displacements of alpine species that were additionally affected by the rapid upward movement of trees and shrubs. Since the beginning of the Neolithic, vegetation changes at Sanetsch Pass resulted from a combination of climate change and human impact. Anthropogenic fire increase and land-use change combined with a natural change from subcontinental to more oceanic climate during the second half of the Holocene led to the disappearance of P. cembra in the study area, but favoured the occurrence of Picea abies and Alnus viridis. The mid- to late-Holocene decline of Abies alba was primarily a consequence of human impact, since this mesic species should have benefitted from a shift to more oceanic conditions. Future alpine vegetation changes will be a function of the amplitude and rapidity of global warming as well as human land use. Our results imply that alpine vegetation at our treeline pass site was never replaced by forests since the last ice-age. This may change in the future if anticipated climate change will induce upslope migration of trees. The results of this study emphasise the necessity of climate change mitigation in order to prevent biodiversity losses as a consequence of unprecedented community and species displacement in response to climatic change.

Impact of suborbital climate changes in the North Atlantic on ice sheet dynamics at the Mid-Pleistocene Transition

[1] Early and Mid-Pleistocene climate, ocean hydrography and ice sheet dynamics have been reconstructed using a high-resolution data set (planktonic and benthicδ18O time series, faunal-based sea surface temperature (SST) reconstructions and ice-rafted debris (IRD)) record from a high-deposition-rate sedimentary succession recovered at the Gardar Drift formation in the subpolar North Atlantic (Integrated Ocean Drilling Program Leg 306, Site U1314). Our sedimentary record spans from late in Marine Isotope Stage (MIS) 31 to MIS 19 (1069–779 ka). Different trends of the benthic and planktonic oxygen isotopes, SST and IRD records before and after MIS 25 (∼940 ka) evidence the large increase in Northern Hemisphere ice-volume, linked to the cyclicity change from the 41-kyr to the 100-kyr that occurred during the Mid-Pleistocene Transition (MPT). Beside longer glacial-interglacial (G-IG) variability, millennial-scale fluctuations were a pervasive feature across our study. Negative excursions in the benthicδ18O time series observed at the times of IRD events may be related to glacio-eustatic changes due to ice sheets retreats and/or to changes in deep hydrography. Time series analysis on surface water proxies (IRD, SST and planktonicδ18O) of the interval between MIS 31 to MIS 26 shows that the timing of these millennial-scale climate changes are related to half-precessional (10 kyr) components of the insolation forcing, which are interpreted as cross-equatorial heat transport toward high latitudes during both equinox insolation maxima at the equator.

Glacial–interglacial dynamics of Antarctic firn columns: comparison between simulations and ice core air-δ15N measurements

Correct estimation of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice core studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: outputs of a firn densification model, and measurements of δ15N of N2 in air trapped in ice core, assuming that δ15N is only affected by gravitational fractionation in the firn column. Since the firn densification process is largely governed by surface temperature and accumulation rate, we have investigated four ice cores drilled in coastal (Berkner Island, BI, and James Ross Island, JRI) and semi-coastal (TALDICE and EPICA Dronning Maud Land, EDML) Antarctic regions. Combined with available ice core air-δ15N measurements from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rates and temperature conditions. Our δ15N profiles reveal a heterogeneous response of the firn structure to glacial–interglacial climatic changes. While firn densification simulations correctly predict TALDICE δ15N variations, they systematically fail to capture the large millennial-scale δ15N variations measured at BI and the δ15N glacial levels measured at JRI and EDML – a mismatch previously reported for central East Antarctic ice cores. New constraints of the EDML gas–ice depth offset during the Laschamp event (~41 ka) and the last deglaciation do not favour the hypothesis of a large convective zone within the firn as the explanation of the glacial firn model–δ15N data mismatch for this site. While we could not conduct an in-depth study of the influence of impurities in snow for firnification from the existing datasets, our detailed comparison between the δ15N profiles and firn model simulations under different temperature and accumulation rate scenarios suggests that the role of accumulation rate may have been underestimated in the current description of firnification models.

Estimated strength of the Atlantic overturning circulation during the last deglaciation

The Atlantic meridional overturning circulation affects the latitudinal distribution of heat, and is a key component of the climate system. Proxy reconstructions, based on sedimentary Pa-231/Th-230 ratios and the difference between surface-and deep-water radiocarbon ages, indicate that during the last glacial period, the overturning circulation was reduced during millennial-scale periods of cooling(1-5). However, much debate exists over the robustness of these proxies(6-8). Here e combine proxy reconstructions of sea surface and air temperatures and a global climate model to quantitatively estimate changes in the strength of the Atlantic meridional overturning circulation during the last glacial period. We find that, relative to the Last Glacial Maximum, the overturning circulation was reduced by approximately 14 Sv during the cold Heinrich event 1. During the Younger Dryas cold event, the overturning circulation was reduced by approximately 12 Sv, relative to the preceding warm interval. These changes are consistent with qualitative estimates of the overturning circulation from sedimentary Pa-231/Th-230 ratios. In addition, we find that the strength of the overturning circulation during the Last Glacial Maximum and the Holocene epoch are indistinguishable within the uncertainty of the reconstruction.

High-resolution glacial and deglacial record of atmospheric methane by continuous-flow and laser spectrometer analysis along the NEEM ice core

The Greenland NEEM (North Greenland Eemian Ice Drilling) operation in 2010 provided the first opportunity to combine trace-gas measurements by laser spectroscopic instruments and continuous-flow analysis along a freshly drilled ice core in a field-based setting. We present the resulting atmospheric methane (CH4) record covering the time period from 107.7 to 9.5 ka b2k (thousand years before 2000 AD). Companion discrete CH4 measurements are required to transfer the laser spectroscopic data from a relative to an absolute scale. However, even on a relative scale, the high-resolution CH4 data set significantly improves our knowledge of past atmospheric methane concentration changes. New significant sub-millennial-scale features appear during interstadials and stadials, generally associated with similar changes in water isotopic ratios of the ice, a proxy for local temperature. In addition to the midpoint of Dansgaard–Oeschger (D/O) CH4 transitions usually used for cross-dating, sharp definition of the start and end of these events brings precise depth markers (with ±20 cm uncertainty) for further cross-dating with other palaeo- or ice core records, e.g. speleothems. The method also provides an estimate of CH4 rates of change. The onsets of D/O events in the methane signal show a more rapid rate of change than their endings. The rate of CH4 increase associated with the onsets of D/O events progressively declines from 1.7 to 0.6 ppbv yr−1 in the course of marine isotope stage 3. The largest observed rate of increase takes place at the onset of D/O event #21 and reaches 2.5 ppbv yr−1.