Search

Abstract. Deglaciations are characterized by relatively fast and near-synchronous changes in ice sheet volume, ocean temperature, and atmospheric greenhouse gas concentrations, but glacial inceptions occur more gradually. Understanding the evolution of ice sheet, ocean, and atmospheric conditions from interglacial to glacial maximum provides important insight into the interplay of these components of our climate system. Using noble gas measurements in ancient ice samples, we reconstruct mean ocean temperature (MOT) from 74 to 59.5 ka BP, covering the Marine Isotope Stage (MIS) 5-4 boundary, MIS 4, and part of the MIS 4-3 transition. Comparing this MOT reconstruction to previously published MOT reconstructions from the last glacial cycle, we find that the majority of interglacial-glacial ocean cooling occurred across MIS 5, and MOT reached full glacial levels by MIS 4 (−2.7 ± 0.3 °C relative to the Holocene). Comparing MOT to contemporaneous records of CO2 and benthic ?18O, we find that ocean cooling and the solubility pump can explain most of the CO2 drawdown and increase in ?18O across MIS 5. The timing of ocean warming and cooling in our record indicates that millennial scale climate variability plays a crucial role in setting mean ocean temperature during this interval, as seen during other periods, such as the last deglaciation.

Abstract. Deglaciations are characterized by relatively fast andnear-synchronous changes in ice sheet volume, ocean temperature, andatmospheric greenhouse gas concentrations, but glacial inception occurs moregradually. Understanding the evolution of ice sheet, ocean, and atmosphereconditions from interglacial to glacial maximum provides insight into theinterplay of these components of the climate system. Using noble gasmeasurements in ancient ice samples, we reconstruct mean ocean temperature(MOT) from 74 to 59.7 ka, covering the Marine Isotope Stage (MIS) 5a–4boundary, MIS 4, and part of the MIS 4–3 transition. Comparing this MOTreconstruction to previously published MOT reconstructions from the last andpenultimate deglaciation, we find that the majority of the lastinterglacial–glacial ocean cooling must have occurred within MIS 5. MOTreached equally cold conditions in MIS 4 as in MIS 2 (−2.7 ± 0.3 ∘C relative to the Holocene, −0.1 ± 0.3 ∘Crelative to MIS 2). Using a carbon cycle model to quantify the CO2 solubility pump, we show that ocean cooling can explain most of theCO2 drawdown (32 ± 4 of 40 ppm) across MIS 5. Comparing MOT tocontemporaneous records of benthic δ18O, we find that ocean coolingcan also explain the majority of the δ18O increase across MIS 5 (0.7 ‰of 1.3 ‰). The timing of ocean warming and cooling inthe record and the comparison to coeval Antarctic isotope data suggest anintimate link between ocean heat content, Southern Hemisphere high-latitude climate,and ocean circulation on orbital and millennial timescales.

The energy imbalance at the top of the atmosphere determines the temporal evolution of the global climate, and vice versa changes in the climate system can alter the planetary energy fluxes. This interplay is fundamental to our understanding of Earth's heat budget and the climate system. However, even today, the direct measurement of global radiative fluxes is difficult, such that most assessments are based on changes in the total energy content of the climate system. We apply the same approach to estimate the long-term evolution of Earth's radiative imbalance in the past. New measurements of noble gas-derived mean ocean temperature from the European Project for Ice Coring in Antarctica Dome C ice core covering the last 40,000 y, combined with recent results from the West Antarctic Ice Sheet Divide ice core and the sea-level record, allow us to quantitatively reconstruct the history of the climate system energy budget. The temporal derivative of this quantity must be equal to the planetary radiative imbalance. During the deglaciation, a positive imbalance of typically +0.2 W⋅m-2 is maintained for ∼10,000 y, however, with two distinct peaks that reach up to 0.4 W⋅m-2 during times of substantially reduced Atlantic Meridional Overturning Circulation. We conclude that these peaks are related to net changes in ocean heat uptake, likely due to rapid changes in North Atlantic deep-water formation and their impact on the global radiative balance, while changes in cloud coverage, albeit uncertain, may also factor into the picture.

New ice cores retrieved from the Taylor Glacier (Antarctica) blue ice area contain ice and air spanning the Marine Isotope Stage (MIS) 5-4 transition, a period of global cooling and ice sheet expansion. We determine chronologies for the ice and air bubbles in the new ice cores by visually matching variations in gas- and ice-phase tracers to preexisting ice core records. The chronologies reveal an ice age-gas age difference (Δage) approaching 10 ka during MIS 4, implying very low snow accumulation in the Taylor Glacier accumulation zone. A revised chronology for the analogous section of the Taylor Dome ice core (84 to 55 ka), located to the south of the Taylor Glacier accumulation zone, shows that Δage did not exceed 3 ka. The difference in Δage between the two records during MIS 4 is similar in magnitude but opposite in direction to what is observed at the Last Glacial Maximum. This relationship implies that a spatial gradient in snow accumulation existed across the Taylor Dome region during MIS 4 that was oriented in the opposite direction of the accumulation gradient during the Last Glacial Maximum.

Abstract. A new ice core retrieved from the Taylor Glacier blue ice area contains ice and air spanning the Marine Isotope Stage (MIS) 5/4 transition (74 to 65 ka), a period of global cooling and glacial inception. Dating the ice and air bubbles in the new ice core reveals an ice age-gas age difference (Δage) approaching 10 ka during MIS 4, implying very low accumulation at the Taylor Glacier accumulation zone on the northern flank of Taylor Dome. A revised chronology for the Taylor Dome ice core (80 to 55 ka), situated to the south of the Taylor Glacier accumulation zone, shows that Δage did not exceed 2.5 ka at that location. The difference in Δage between the new Taylor Glacier ice core and the Taylor Dome ice core implies a spatial gradient in snow accumulation across Taylor Dome that intensified during the last glacial inception and through MIS 4.

Little is known about the ocean temperature's long-term response to climate perturbations owing to limited observations and a lack of robust reconstructions. Although most of the anthropogenic heat added to the climate system has been taken up by the ocean up until now, its role in a century and beyond is uncertain. Here, using noble gases trapped in ice cores, we show that the mean global ocean temperature increased by 2.57 ± 0.24 degrees Celsius over the last glacial transition (20,000 to 10,000 years ago). Our reconstruction provides unprecedented precision and temporal resolution for the integrated global ocean, in contrast to the depth-, region-, organism- and season-specific estimates provided by other methods. We find that the mean global ocean temperature is closely correlated with Antarctic temperature and has no lead or lag with atmospheric CO2, thereby confirming the important role of Southern Hemisphere climate in global climate trends. We also reveal an enigmatic 700-year warming during the early Younger Dryas period (about 12,000 years ago) that surpasses estimates of modern ocean heat uptake.

Glacial-state greenhouse gas concentrations and Southern Hemisphere climate conditions persisted until ∼17.7 ka, when a nearly synchronous acceleration in deglaciation was recorded in paleoclimate proxies in large parts of the Southern Hemisphere, with many changes ascribed to a sudden poleward shift in the Southern Hemisphere westerlies and subsequent climate impacts. We used high-resolution chemical measurements in the West Antarctic Ice Sheet Divide, Byrd, and other ice cores to document a unique, ∼192-y series of halogen-rich volcanic eruptions exactly at the start of accelerated deglaciation, with tephra identifying the nearby Mount Takahe volcano as the source. Extensive fallout from these massive eruptions has been found >2,800 km from Mount Takahe. Sulfur isotope anomalies and marked decreases in ice core bromine consistent with increased surface UV radiation indicate that the eruptions led to stratospheric ozone depletion. Rather than a highly improbable coincidence, circulation and climate changes extending from the Antarctic Peninsula to the subtropics-similar to those associated with modern stratospheric ozone depletion over Antarctica-plausibly link the Mount Takahe eruptions to the onset of accelerated Southern Hemisphere deglaciation ∼17.7 ka.

Methane (CH4) is a powerful greenhouse gas and plays a key part in global atmospheric chemistry. Natural geological emissions (fossil methane vented naturally from marine and terrestrial seeps and mud volcanoes) are thought to contribute around 52 teragrams of methane per year to the global methane source, about 10 per cent of the total, but both bottom-up methods (measuring emissions) and top-down approaches (measuring atmospheric mole fractions and isotopes) for constraining these geological emissions have been associated with large uncertainties. Here we use ice core measurements to quantify the absolute amount of radiocarbon-containing methane (14CH4) in the past atmosphere and show that geological methane emissions were no higher than 15.4 teragrams per year (95 per cent confidence), averaged over the abrupt warming event that occurred between the Younger Dryas and Preboreal intervals, approximately 11,600 years ago. Assuming that past geological methane emissions were no lower than today, our results indicate that current estimates of today's natural geological methane emissions (about 52 teragrams per year) are too high and, by extension, that current estimates of anthropogenic fossil methane emissions are too low. Our results also improve on and confirm earlier findings that the rapid increase of about 50 per cent in mole fraction of atmospheric methane at the Younger Dryas-Preboreal event was driven by contemporaneous methane from sources such as wetlands; our findings constrain the contribution from old carbon reservoirs (marine methane hydrates, permafrost and methane trapped under ice) to 19 per cent or less (95 per cent confidence). To the extent that the characteristics of the most recent deglaciation and the Younger Dryas-Preboreal warming are comparable to those of the current anthropogenic warming, our measurements suggest that large future atmospheric releases of methane from old carbon sources are unlikely to occur.

Abstract. Old ice for paleo-environmental studies, traditionally accessed through deep core drilling on domes and ridges on the large ice sheets, can also be retrieved at the surface from ice sheet margins and blue ice areas. The practically unlimited amount of ice available at these sites satisfies a need in the community for studies of trace components requiring large sample volumes. For margin sites to be useful as ancient ice archives, the ice stratigraphy needs to be understood and age models need to be established. We present measurements of trapped gases in ice from Taylor Glacier, Antarctica, to date the ice and assess the completeness of the stratigraphic section. Using δ18O of O2 and methane concentrations, we unambiguously identify ice from the last glacial cycle, covering every climate interval from the early Holocene to the penultimate interglacial. A high-resolution transect reveals the last deglaciation and the Last Glacial Maximum (LGM) in detail. We observe large-scale deformation in the form of folding, but individual stratigraphic layers do not appear to have undergone irregular thinning. Rather, it appears that the entire LGM-deglaciation sequence has been transported from the interior of the ice sheet to the surface of Taylor Glacier relatively undisturbed. We present an age model that builds the foundation for gas studies on Taylor Glacier. A comparison with the Taylor Dome ice core confirms that the section we studied on Taylor Glacier is better suited for paleo-climate reconstructions of the LGM due to higher accumulation rates.

Abstract. Marine sediments, speleothems, paleo-lake elevations, and ice core methane and δ18O of O2 (δ18Oatm) records provide ample evidence for repeated abrupt meridional shifts in tropical rainfall belts throughout the last glacial cycle. To improve understanding of the impact of abrupt events on the global terrestrial biosphere, we present composite records of δ18Oatm and inferred changes in fractionation by the global terrestrial biosphere (ΔεLAND) from discrete gas measurements in the WAIS Divide (WD) and Siple Dome (SD) Antarctic ice cores. On the common WD timescale, it is evident that maxima in ΔεLAND are synchronous with or shortly follow small-amplitude WD CH4 peaks that occur within Heinrich stadials 1, 2, 4, and 5 – periods of low atmospheric CH4 concentrations. These local CH4 maxima have been suggested as markers of abrupt climate responses to Heinrich events. Based on our analysis of the modern seasonal cycle of gross primary productivity (GPP)-weighted δ18O of terrestrial precipitation (the source water for atmospheric O2 production), we propose a simple mechanism by which ΔεLAND tracks the centroid latitude of terrestrial oxygen production. As intense rainfall and oxygen production migrate northward, ΔεLAND should decrease due to the underlying meridional gradient in rainfall δ18O. A southward shift should increase ΔεLAND. Monsoon intensity also influences δ18O of precipitation, and although we cannot determine the relative contributions of the two mechanisms, both act in the same direction. Therefore, we suggest that abrupt increases in ΔεLAND unambiguously imply a southward shift of tropical rainfall. The exact magnitude of this shift, however, remains under-constrained by ΔεLAND.

An understanding of the mechanisms that control CO2 change during glacial-interglacial cycles remains elusive. Here we help to constrain changing sources with a high-precision, high-resolution deglacial record of the stable isotopic composition of carbon in CO2(δ(13)C-CO2) in air extracted from ice samples from Taylor Glacier, Antarctica. During the initial rise in atmospheric CO2 from 17.6 to 15.5 ka, these data demarcate a decrease in δ(13)C-CO2, likely due to a weakened oceanic biological pump. From 15.5 to 11.5 ka, the continued atmospheric CO2 rise of 40 ppm is associated with small changes in δ(13)C-CO2, consistent with a nearly equal contribution from a further weakening of the biological pump and rising ocean temperature. These two trends, related to marine sources, are punctuated at 16.3 and 12.9 ka with abrupt, century-scale perturbations in δ(13)C-CO2 that suggest rapid oxidation of organic land carbon or enhanced air-sea gas exchange in the Southern Ocean. Additional century-scale increases in atmospheric CO2 coincident with increases in atmospheric CH4 and Northern Hemisphere temperature at the onset of the Bølling (14.6-14.3 ka) and Holocene (11.6-11.4 ka) intervals are associated with small changes in δ(13)C-CO2, suggesting a combination of sources that included rising surface ocean temperature.

The last glacial period exhibited abrupt Dansgaard-Oeschger climatic oscillations, evidence of which is preserved in a variety of Northern Hemisphere palaeoclimate archives. Ice cores show that Antarctica cooled during the warm phases of the Greenland Dansgaard-Oeschger cycle and vice versa, suggesting an interhemispheric redistribution of heat through a mechanism called the bipolar seesaw. Variations in the Atlantic meridional overturning circulation (AMOC) strength are thought to have been important, but much uncertainty remains regarding the dynamics and trigger of these abrupt events. Key information is contained in the relative phasing of hemispheric climate variations, yet the large, poorly constrained difference between gas age and ice age and the relatively low resolution of methane records from Antarctic ice cores have so far precluded methane-based synchronization at the required sub-centennial precision. Here we use a recently drilled high-accumulation Antarctic ice core to show that, on average, abrupt Greenland warming leads the corresponding Antarctic cooling onset by 218 ± 92 years (2σ) for Dansgaard-Oeschger events, including the Bølling event; Greenland cooling leads the corresponding onset of Antarctic warming by 208 ± 96 years. Our results demonstrate a north-to-south directionality of the abrupt climatic signal, which is propagated to the Southern Hemisphere high latitudes by oceanic rather than atmospheric processes. The similar interpolar phasing of warming and cooling transitions suggests that the transfer time of the climatic signal is independent of the AMOC background state. Our findings confirm a central role for ocean circulation in the bipolar seesaw and provide clear criteria for assessing hypotheses and model simulations of Dansgaard-Oeschger dynamics.

Interpretation of ice core trace gas records depends on an accurate understanding of the processes that smooth the atmospheric signal in the firn. Much work has been done to understand the processes affecting air transport in the open pores of the firn, but a paucity of data from air trapped in bubbles in the firn-ice transition region has limited the ability to constrain the effect of bubble closure processes. Here we present high-resolution measurements of firn density, methane concentrations, nitrogen isotopes, and total air content that show layering in the firn-ice transition region at the West Antarctic Ice Sheet (WAIS) Divide ice core site. Using the notion that bubble trapping is a stochastic process, we derive a new parameterization for closed porosity that incorporates the effects of layering in a steady state firn modeling approach. We include the process of bubble trapping into an open-porosity firn air transport model and obtain a good fit to the firn core data. We find that layering broadens the depth range over which bubbles are trapped, widens the modeled gas age distribution of air in closed bubbles, reduces the mean gas age of air in closed bubbles, and introduces stratigraphic irregularities in the gas age scale that have a peak-to-peak variability of ~10 years at WAIS Divide. For a more complete understanding of gas occlusion and its impact on ice core records, we suggest that this experiment be repeated at sites climatically different from WAIS Divide, for example, on the East Antarctic plateau.

Nitrous oxide (N2O) is an important greenhouse gas and ozone-depleting substance that has anthropogenic as well as natural marine and terrestrial sources. The tropospheric N2O concentrations have varied substantially in the past in concert with changing climate on glacial-interglacial and millennial timescales. It is not well understood, however, how N2O emissions from marine and terrestrial sources change in response to varying environmental conditions. The distinct isotopic compositions of marine and terrestrial N2O sources can help disentangle the relative changes in marine and terrestrial N2O emissions during past climate variations. Here we present N2O concentration and isotopic data for the last deglaciation, from 16,000 to 10,000 years before present, retrieved from air bubbles trapped in polar ice at Taylor Glacier, Antarctica. With the help of our data and a box model of the N2O cycle, we find a 30 per cent increase in total N2O emissions from the late glacial to the interglacial, with terrestrial and marine emissions contributing equally to the overall increase and generally evolving in parallel over the last deglaciation, even though there is no a priori connection between the drivers of the two sources. However, we find that terrestrial emissions dominated on centennial timescales, consistent with a state-of-the-art dynamic global vegetation and land surface process model that suggests that during the last deglaciation emission changes were strongly influenced by temperature and precipitation patterns over land surfaces. The results improve our understanding of the drivers of natural N2O emissions and are consistent with the idea that natural N2O emissions will probably increase in response to anthropogenic warming.

Author(s): Vasilii V. Petrenko (corresponding author) [1]; Andrew M. Smith [2]; Hinrich Schaefer [3]; Katja Riedel [3]; Edward Brook [4]; Daniel Baggenstos [5, 6]; Christina Harth [5]; Quan Hua [2]; [...]

Precision, accuracy, and temporal resolution are key to make full use of atmospheric trace gas records in ice cores. These aspects will become especially crucial for ice cores that aim to extend the ice core record to the last 1.5 Myr, i.e., across the Mid Pleistocene Transition (as currently drilled within the European project Beyond EPICA - Oldest Ice Core (BE-OIC)). The ice from this period is expected to be close to bedrock and, due to glacier flow, extremely thinned with 15,000 years of climate history contained in only one meter of ice. Accordingly, for a century-scale resolution, the sample vertical extent must be reduced to a few cm containing only about 1-2mL air STP. We present a novel combined system for the extraction and the simultaneous measurement of CO2, CH4, and N2O concentrations, as well as δ13CO2, which achieves a vertical resolution of 1-2 cm with precisions of 0.4 ppm, 3 ppb, 1 ppb and 0.04 ‰, respectively. This is accomplished by employing a directional and continuous laser induced sublimation followed by analysis of the sample gas by quantum cascade laser absorption spectroscopy (QCLAS). Besides the low sample volume requirements and the vertical resolution capabilities, the described method holds additional advantages over previous methods, including the immunity of the highly specific QCLAS analysis to drilling fluid contamination as well as the non-destructive nature of the spectroscopic gas analysis. The combined extraction and analysis system was extensively tested by sublimating gas-free ice with introduction of a standard gas to determine the accuracy and characterize potential artefacts. Moreover, Antarctic ice samples were measured to confirm the measurement performance, covering the range of variability expected in Pleistocene ice and to highlight the vertical resolution capabilities critical for its application within BE-OIC. [ABSTRACT FROM AUTHOR]

Here we present a newly developed ice core gas-phase proxy that directly samples a component of the large-scale atmospheric circulation: synoptic-scale pressure variability. Surface pressure variability weakly disrupts gravitational isotopic settling in the firn layer, which is recorded in krypton-86 excess (86Krxs). We validate 86Krxs using late Holocene ice samples from eleven Antarctic and one Greenland ice core that collectively represent a wide range of surface pressure variability in the modern climate. We find a strong correlation (r = -0.94, p < 0.01) between site-average 86Krxs and site synoptic variability from reanalysis data. The main uncertainties in the method are the corrections for gas loss and thermal fractionation, and the relatively large scatter in the data. We show 86Krxs is linked to the position of the eddy-driven subpolar jet (SPJ), with a southern position enhancing pressure variability. We present a 86Krxs record covering the last 24 ka from the WAIS Divide ice core. West Antarctic synoptic activity is slightly below modern levels during the last glacial maximum (LGM); increases during the Heinrich Stadial 1 and Younger Dryas North Atlantic cold periods; weakens abruptly at the Holocene onset; remains low during the early and mid-Holocene, and gradually increases to its modern value. The WAIS Divide 86Krxs record resembles records of monsoon intensity thought to reflect changes in the meridional position of the intertropical convergence zone (ITCZ) on orbital and millennial timescales, such that West Antarctic storminess is weaker when the ITCZ is displaced northward, and stronger when it is displaced southward. We interpret variations in synoptic activity as reflecting movement of the South Pacific SPJ in parallel to the ITCZ migrations, which is the expected zonal-mean response of the eddy-driven jet in models and proxy data. Past changes to Pacific climate and the El Niño Southern Oscillation (ENSO) may amplify the signal of the SPJ migration. Our interpretation is broadly consistent with opal flux records from the Pacific Antarctic zone thought to reflect wind-driven upwelling. We emphasize that 86Krxs is a new proxy, and more work is called for to confirm, replicate and better understand these results; until such time, our conclusions regarding past atmospheric dynamics remain tentative. Current scientific understanding of firn air transport and trapping is insufficient to explain all the observed variations in 86Krxs. [ABSTRACT FROM AUTHOR]

Together with the latent heat stored in glacial ice sheets the ocean heat uptake carries the lion's share of glacial/interglacial changes in the planetary heat content but little direct information on the global mean ocean temperature (MOT) is available to constrain the ocean temperature response to glacial/interglacial climate perturbations. Using ratios of noble gases and molecular nitrogen trapped in the Antarctic EPICA Dome C ice core we are able to reconstruct MOT for peak glacial and interglacial conditions during the last 700,000 years and explore the differences between these extrema. To this end, we have to correct the noble gas ratios for gas transport effects in the firn column and gas loss fractionation processes of the samples after ice core retrieval using the full elemental matrix of N2, Ar, Kr and Xe in the ice and their individual isotopic ratios. The reconstructed MOT in peak glacials is consistently about 3.3 ± 0.4 °C cooler compared to the Holocene. Lukewarm interglacials before the Mid Brunhes event 450 kyr ago are characterized by 1.6 ± 0.4 °C lower temperatures in the ocean than the Holocene, thus, glacial/interglacial amplitudes were only about 50 % of those after the Mid Brunhes event, in line with the reduced radiative forcing by lower greenhouse gas concentrations and their Earth system feedbacks. Moreover, we find significantly increased MOTs at the onset of Marine Isotope Stage 5.5 and 9.3, which are coeval with CO2 and CH4 overshoots at that time. We link these CO2 and CH4 overshoots to a resumption of the Atlantic Meridional Overturning Circulation, which is also the starting point of the release of heat previously accumulated in the ocean during times of reduced overturning. [ABSTRACT FROM AUTHOR]