Search

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Permafrost and methane hydrates are large, climate-sensitive old carbon reservoirs that have the potential to emit large quantities of methane, a potent greenhouse gas, as the Earth continues to warm. We present ice core isotopic measurements of methane (Δ14C, δ13C, and δD) from the last deglaciation, which is a partial analog for modern warming. Our results show that methane emissions from old carbon reservoirs in response to deglacial warming were small (

The Last Interglacial (129–116 thousand years ago (ka)) represents one of the warmest climate intervals of the past 800,000 years and the most recent time when sea level was metres higher than today. However, the timing and magnitude of the peak warmth varies between reconstructions, and the relative importance of individual sources that contribute to the elevated sea level (mass gain versus seawater expansion) during the Last Interglacial remains uncertain. Here we present the first mean ocean temperature record for this interval from noble gas measurements in ice cores and constrain the thermal expansion contribution to sea level. Mean ocean temperature reached its maximum value of 1.1 ± 0.3 °C warmer-than-modern values at the end of the penultimate deglaciation at 129 ka, which resulted in 0.7 ± 0.3 m of thermosteric sea-level rise relative to present level. However, this maximum in ocean heat content was a transient feature; mean ocean temperature decreased in the first several thousand years of the interglacial and achieved a stable, comparable-to-modern value by ~127 ka. The synchroneity of the peak in mean ocean temperature with proxy records of abrupt transitions in the oceanic and atmospheric circulation suggests that the mean ocean temperature maximum is related to the accumulation of heat in the ocean interior during the preceding period of reduced overturning circulation.

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.

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.

The West Antarctic Ice Sheet Divide (WAIS Divide, WD) ice core is a newly drilled, high-accumulation deep ice core that provides Antarctic climate records of the past ~68 ka at unprecedented temporal resolution. The upper 2850 m (back to 31.2 ka BP) have been dated using annual-layer counting. Here we present a chronology for the deep part of the core (67.8-31.2 ka BP), which is based on stratigraphic matching to annual-layer-counted Greenland ice cores using globally well-mixed atmospheric methane. We calculate the WD gas age-ice age difference (Δage) using a combination of firn densification modeling, ice-flow modeling, and a data set of δ15N-N2, a proxy for past firn column thickness. The largest Δage at WD occurs during the Last Glacial Maximum, and is 525 ± 120 years. Internally consistent solutions can be found only when assuming little to no influence of impurity content on densification rates, contrary to a recently proposed hypothesis. We synchronize the WD chronology to a linearly scaled version of the layer-counted Greenland Ice Core Chronology (GICC05), which brings the age of Dansgaard-Oeschger (DO) events into agreement with the U/Th absolutely dated Hulu Cave speleothem record. The small Δage at WD provides valuable opportunities to investigate the timing of atmospheric greenhouse gas variations relative to Antarctic climate, as well as the interhemispheric phasing of the "bipolar seesaw".

Abstract. The West Antarctic Ice Sheet (WAIS)-Divide ice core (WAIS-D) is a newly drilled, high-accumulation deep ice core that provides Antarctic climate records of the past ∼68 ka at unprecedented temporal resolution. The upper 2850 m (back to 31.2 ka BP) have been dated using annual-layer counting. Here we present a chronology for the deep part of the core (67.8–31.2 ka BP), which is based on stratigraphic matching to annual-layer-counted Greenland ice cores using globally well-mixed atmospheric methane. We calculate the WAIS-D gas age-ice age difference (Δage) using a combination of firn densification modeling, ice flow modeling, and a dataset of δ15N-N2, a proxy for past firn column thickness. The largest Δage at WAIS-D occurs during the last glacial maximum, and is 525 ± 100 years. Internally consistent solutions can only be found when assuming little-to-no influence of impurity content on densification rates, contrary to a recently proposed hypothesis. We synchronize the WAIS-D chronology to a linearly scaled version of the layer-counted Greenland Ice Core Chronology (GICC05), which brings the age of Dansgaard-Oeschger (DO) events into agreement with the U/Th absolutely dated Hulu speleothem record. The small Δage at WAIS-D provides valuable opportunities to investigate the timing of atmospheric greenhouse gas variations relative to Antarctic climate, as well as the interhemispheric phasing of the bipolar "seesaw".

The atmospheric abundance of 81Kr is a global integrator of cosmic rays. It is insensitive to climate shifts, geographical variations, and short‐term solar cycle activity, making it an ideal standard to test models of cosmic‐ray flux on the time scale of 105 years. Here we present the first calculation of absolute 81Kr production rates in the atmosphere, and a measurement of the atmospheric 81Kr/Kr abundance via the Atom Trap Trace Analysis method. The measurement result significantly deviates from previously reported values. The agreement between measurement and model prediction supports the current understanding of the production mechanisms. Additionally, the calculated 81Kr atmospheric inventory over the past 1.5 Myr provides a more accurate input function for radiokrypton dating. Plain Language Summary: Krypton‐81 is a long‐lived radioactive isotope produced in the Earth's atmosphere by cosmic rays. It stays in the atmosphere as a noble gas for hundreds of thousands of years until its eventual nuclear decay. As a result, its abundance uniquely reflects the long‐term accumulation record of cosmic rays across the entire globe. We performed the first precise measurement of the atmospheric abundance of krypton‐81. The result agrees with the prediction of a realistic isotope production model, thus confirming the current understanding of the cosmic‐ray flux, isotope production mechanisms, and the past terrestrial and space magnetic field environment. Key Points: Conducted new precision measurement of 81Kr isotopic abundance in the atmosphere; new result differs significantly from previous valuesFirst ever cosmic‐ray flux models of 81Kr production in the atmosphere are presented; models successfully predict new measured valueMeasurement and models provide more accurate input function for radiokrypton dating over the past 1.5 Myr [ABSTRACT FROM AUTHOR]

Abstract: Changes in atmospheric CO2 on millennial‐to‐centennial timescales are key components of past climate variability during the last glacial and deglacial periods (70–10 ka), yet the sources and mechanisms responsible for the CO2 fluctuations remain largely obscure. Here we report the 13C/12C ratio of atmospheric CO2 during a key interval of the last glacial period at submillennial resolution, with coeval histories of atmospheric CO2, CH4, and N2O concentrations. The carbon isotope data suggest that the millennial‐scale CO2 variability in Marine Isotope Stage 3 is driven largely by changes in the organic carbon cycle, most likely by sequestration of respired carbon in the deep ocean. Centennial‐scale CO2 variations, distinguished by carbon isotope signatures, are associated with both abrupt hydrological change in the tropics (e.g., Heinrich events) and rapid increases in Northern Hemisphere temperature (Dansgaard‐Oeschger events). These events can be linked to modes of variability during the last deglaciation, thus suggesting that drivers of millennial and centennial CO2 variability during both periods are intimately linked to abrupt climate variability. [ABSTRACT FROM AUTHOR]

The 36Ar(n,γ)37Ar (t1/2=35 d) and 38Ar(n,γ)39Ar (269 yr) reactions were studied for the first time with a quasi-Maxwellian (kT~47 keV) neutron flux for Maxwellian average cross section (MACS) measurements at stellar energies. Gas samples were irradiated at the high-intensity Soreq applied research accelerator facility-liquid-lithium target neutron source and the 37Ar/36Ar and 39Ar/38Ar ratios in the activated samples were determined by accelerator mass spectrometry at the ATLAS facility (Argonne National Laboratory). The 37Ar activity was also measured by low-level counting at the University of Bern. Experimental MACS of 36Ar and 38Ar, corrected to the standard 30 keV thermal energy, are 1.9(3) and 1.3(2) mb, respectively, differing from the theoretical and evaluated values published to date by up to an order of magnitude. The neutron-capture cross sections of 36,38Ar are relevant to the stellar nucleosynthesis of light neutron-rich nuclides; the two experimental values are shown to affect the calculated mass fraction of nuclides in the region A=36-48 during the weak s process. The new production cross sections have implications also for the use of 37Ar and 39Ar as environmental tracers in the atmosphere and hydrosphere. [ABSTRACT FROM AUTHOR]

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., Competing Interests: The authors declare no conflict of interest.

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 CO 2 , 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., Competing Interests: The authors declare no conflict of interest.

Methane (CH 4 ) 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 ( 14 CH 4 ) 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.

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.

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.

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.

Nitrogen trifluoride (NF(3)) has potential to make a growing contribution to the Earth's radiative budget; however, our understanding of its atmospheric burden and emission rates has been limited. Based on a revision of our previous calibration and using an expanded set of atmospheric measurements together with an atmospheric model and inverse method, we estimate that the global emissions of NF(3) in 2011 were 1.18 ± 0.21 Gg⋅y(-1), or ∼20 Tg CO(2)-eq⋅y(-1) (carbon dioxide equivalent emissions based on a 100-y global warming potential of 16,600 for NF(3)). The 2011 global mean tropospheric dry air mole fraction was 0.86 ± 0.04 parts per trillion, resulting from an average emissions growth rate of 0.09 Gg⋅y(-2) over the prior decade. In terms of CO(2) equivalents, current NF(3) emissions represent between 17% and 36% of the emissions of other long-lived fluorinated compounds from electronics manufacture. We also estimate that the emissions benefit of using NF(3) over hexafluoroethane (C(2)F(6)) in electronics manufacture is significant-emissions of between 53 and 220 Tg CO(2)-eq⋅y(-1) were avoided during 2011. Despite these savings, total NF(3) emissions, currently ∼10% of production, are still significantly larger than expected assuming global implementation of ideal industrial practices. As such, there is a continuing need for improvements in NF(3) emissions reduction strategies to keep pace with its increasing use and to slow its rising contribution to anthropogenic climate forcing.