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 and near-synchronous changes in ice sheet volume, ocean temperature, and atmospheric greenhouse gas concentrations, but glacial inception occurs more gradually. Understanding the evolution of ice sheet, ocean, and atmosphere conditions from interglacial to glacial maximum provides insight into the interplay of these components of the climate system. Using noble gas measurements in ancient ice samples, we reconstruct mean ocean temperature (MOT) from 74 to 59.7 ka, covering the Marine Isotope Stage (MIS) 5a–4 boundary, MIS 4, and part of the MIS 4–3 transition. Comparing this MOT reconstruction to previously published MOT reconstructions from the last and penultimate deglaciation, we find that the majority of the last interglacial–glacial ocean cooling must have occurred within MIS 5. MOT reached equally cold conditions in MIS 4 as in MIS 2 (−2.7 ± 0.3 ∘C relative to the Holocene, −0.1 ± 0.3 ∘C relative to MIS 2). Using a carbon cycle model to quantify the CO2 solubility pump, we show that ocean cooling can explain most of the CO2 drawdown (32 ± 4 of 40 ppm) across MIS 5. Comparing MOT to contemporaneous records of benthic δ18O, we find that ocean cooling can also explain the majority of the δ18O increase across MIS 5 (0.7 ‰ of 1.3 ‰). The timing of ocean warming and cooling in the record and the comparison to coeval Antarctic isotope data suggest an intimate link between ocean heat content, Southern Hemisphere high-latitude climate, and ocean circulation on orbital and millennial timescales.

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.

We present successful (81)Kr-Kr radiometric dating of ancient polar ice. Krypton was extracted from the air bubbles in four ∼350-kg polar ice samples from Taylor Glacier in the McMurdo Dry Valleys, Antarctica, and dated using Atom Trap Trace Analysis (ATTA). The (81)Kr radiometric ages agree with independent age estimates obtained from stratigraphic dating techniques with a mean absolute age offset of 6 ± 2.5 ka. Our experimental methods and sampling strategy are validated by (i) (85)Kr and (39)Ar analyses that show the samples to be free of modern air contamination and (ii) air content measurements that show the ice did not experience gas loss. We estimate the error in the (81)Kr ages due to past geomagnetic variability to be below 3 ka. We show that ice from the previous interglacial period (Marine Isotope Stage 5e, 130-115 ka before present) can be found in abundance near the surface of Taylor Glacier. Our study paves the way for reliable radiometric dating of ancient ice in blue ice areas and margin sites where large samples are available, greatly enhancing their scientific value as archives of old ice and meteorites. At present, ATTA (81)Kr analysis requires a 40-80-kg ice sample; as sample requirements continue to decrease, (81)Kr dating of ice cores is a future possibility.

Here we use high-precision carbon isotope data (δ¹³C-CO₂) to show atmospheric CO₂ during Marine Isotope Stage 4 (MIS 4, ~70.5-59 ka) was controlled by a succession of millennial-scale processes. Enriched δ¹³C-CO₂ during peak glaciation suggests increased ocean carbon storage. Variations in δ¹³C-CO₂ in early MIS 4 suggest multiple processes were active during CO₂ drawdown, potentially including decreased land carbon and decreased Southern Ocean air-sea gas exchange superposed on increased ocean carbon storage. CO₂ remained low during MIS 4 while δ¹³C-CO₂ fluctuations suggest changes in Southern Ocean and North Atlantic air-sea gas exchange. A 7 ppm increase in CO₂ at the onset of Dansgaard-Oeschger event 19 (72.1 ka) and 27 ppm increase in CO₂ during late MIS 4 (Heinrich Stadial 6, ~63.5-60 ka) involved additions of isotopically light carbon to the atmosphere. The terrestrial biosphere and Southern Ocean air-sea gas exchange are possible sources, with the latter event also involving decreased ocean carbon storage.

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.

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.

Methane (CH₄) 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 (¹⁴CH₄) 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.

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.