Publications

Export 102 results:
Sort by: Author Title Type [ Year  (Desc)]
2018
Wagner, TJW, Dell RW, Eisenman I, Keeling RF, Padman L, Severinghaus JP.  2018.  Wave inhibition by sea ice enables trans-Atlantic ice rafting of debris during Heinrich events. Earth and Planetary Science Letters. 495:157-163.   10.1016/j.epsl.2018.05.006   AbstractWebsite

The last glacial period was punctuated by episodes of massive iceberg calving from the Laurentide Ice Sheet, called Heinrich events, which are identified by layers of ice-rafted debris (IRD) in ocean sediment cores from the North Atlantic. The thickness of these IRD layers declines more gradually with distance from the iceberg sources than would be expected based on present-day iceberg drift and decay. Here we model icebergs as passive Lagrangian particles driven by ocean currents, winds, and sea surface temperatures. The icebergs are released in a comprehensive climate model simulation of the last glacial maximum (LGM), as well as a simulation of the modern climate. The two simulated climates result in qualitatively similar distributions of iceberg meltwater and hence debris, with the colder temperatures of the LGM having only a relatively small effect on meltwater spread. In both scenarios, meltwater flux falls off rapidly with zonal distance from the source, in contrast with the more uniform spread of IRD in sediment cores. To address this discrepancy, we propose a physical mechanism that could have prolonged the lifetime of icebergs during Heinrich events. The mechanism involves a surface layer of cold and fresh meltwater formed from, and retained around, large densely packed armadas of icebergs. This leads to wintertime sea ice formation even in relatively low latitudes. The sea ice in turn shields the icebergs from wave erosion, which is the main source of iceberg ablation. We find that sea ice could plausibly have formed around the icebergs during four months each winter. Allowing for four months of sea ice in the model results in a simulated IRD distribution which approximately agrees with the distribution of IRD in sediment cores. (C) 2018 Elsevier B.V. All rights reserved.

Birner, B, Buizert C, Wagner TJW, Severinghaus JP.  2018.  The influence of layering and barometric pumping on firn air transport in a 2-D model. Cryosphere. 12:2021-2037.   10.5194/tc-12-2021-2018   AbstractWebsite

Ancient air trapped in ice core bubbles has been paramount to developing our understanding of past climate and atmospheric composition. Before air bubbles become isolated in ice, the atmospheric signal is altered in the firn column by transport processes such as advection and diffusion. However, the influence of low-permeability layers and barometric pumping (driven by surface pressure variability) on firn air transport is not well understood and is not readily captured in conventional one-dimensional (1-D) firn air models. Here we present a two-dimensional (2-D) trace gas advection-diffusion-dispersion model that accounts for discontinuous horizontal layers of reduced permeability. We find that layering or barometric pumping individually yields too small a reduction in gravitational settling to match observations. In contrast, when both effects are active, the model's gravitational fractionation is suppressed as observed. Layering focuses airflows in certain regions in the 2-D model, which acts to amplify the dispersive mixing resulting from barometric pumping. Hence, the representation of both factors is needed to obtain a realistic emergence of the lock-in zone. In contrast to expectations, we find that the addition of barometric pumping in the layered 2-D model does not substantially change the differential kinetic fractionation of fast-and slow-diffusing trace gases. Like 1-D models, the 2-D model substantially underestimates the amount of differential kinetic fractionation seen in actual observations, suggesting that further subgrid-scale processes may be missing in the current generation of firn air transport models. However, we find robust scaling relationships between kinetic isotope fractionation of different noble gas isotope and elemental ratios. These relationships may be used to correct for kinetic fractionation in future high-precision ice core studies and can amount to a bias of up to 0.45 degrees C in noble-gas-based mean ocean temperature reconstructions at WAIS Divide, Antarctica.

Bereiter, B, Kawamura K, Severinghaus JP.  2018.  New methods for measuring atmospheric heavy noble gas isotope and elemental ratios in ice core samples. Rapid Communications in Mass Spectrometry. 32:801-814.   10.1002/rcm.8099   AbstractWebsite

RationaleThe global ocean constitutes the largest heat buffer in the global climate system, but little is known about its past changes. The isotopic and elemental ratios of heavy noble gases (krypton and xenon), together with argon and nitrogen in trapped air from ice cores, can be used to reconstruct past mean ocean temperatures (MOTs). Here we introduce two successively developed methods to measure these parameters with a sufficient precision to provide new constraints on past changes in MOT. MethodsThe air from an 800-g ice sample - containing roughly 80mL STP air - is extracted and processed to be analyzed on two independent dual-inlet isotope ratio mass spectrometers. The primary isotope ratios (N-15, Ar-40 and Kr-86 values) are obtained with precisions in the range of 1 per meg (0.001) per mass unit. The three elemental ratio values Kr/N-2, Xe/N-2 and Xe/Kr are obtained using sequential (non-simultaneous) peak-jumping, reaching precisions in the range of 0.1-0.3. ResultsThe latest version of the method achieves a 30% to 50% better precision on the elemental ratios and a twofold better sample throughput than the previous one. The method development uncovered an unexpected source of artefactual gas fractionation in a closed system that is caused by adiabatic cooling and warming of gases (termed adiabatic fractionation) - a potential source of measurement artifacts in other methods. ConclusionsThe precisions of the three elemental ratios Kr/N-2, Xe/N-2 and Xe/Kr - which all contain the same MOT information - suggest smaller uncertainties for reconstructed MOTs (+/- 0.3-0.1 degrees C) than previous studies have attained. Due to different sensitivities of the noble gases to changes in MOT, Xe/N-2 provides the best constraints on the MOT under the given precisions followed by Xe/Kr, and Kr/N-2; however, using all of them helps to detect methodological artifacts and issues with ice quality.

Bertler, NAN, Conway H, Dahl-Jensen D, Emanuelsson DB, Winstrup M, Vallelonga PT, Lee JE, Brook EJ, Severinghaus JP, Fudge TJ, Keller ED, Baisden WT, Hindmarsh RCA, Neff PD, Blunier T, Edwards R, Mayewski PA, Kipfstuhl S, Buizert C, Canessa S, Dadic R, Kjaer HA, Kurbatov A, Zhang DQ, Waddington ED, Baccolo G, Beers T, Brightley HJ, Carter L, Clemens-Sewall D, Ciobanu VG, Delmonte B, Eling L, Ellis A, Ganesh S, Golledge NR, Haines S, Handley M, Hawley RL, Hogan CM, Johnson KM, Korotkikh E, Lowry DP, Mandeno D, McKay RM, Menking JA, Naish TR, Noerling C, Ollive A, Orsi A, Proemse BC, Pyne AR, Pyne RL, Renwick J, Scherer RP, Semper S, Simonsen M, Sneed SB, Steig EJ, Tuohy A, Venugopal AU, Valero-Delgado F, Venkatesh J, Wang FT, Wang SM, Winski DA, Winton VHL, Whiteford A, Xiao CD, Yang J, Zhang X.  2018.  The Ross Sea Dipole - temperature, snow accumulation and sea ice variability in the Ross Sea region, Antarctica, over the past 2700 years. Climate of the Past. 14:193-214.   10.5194/cp-14-193-2018   AbstractWebsite

High-resolution, well-dated climate archives provide an opportunity to investigate the dynamic interactions of climate patterns relevant for future projections. Here, we present data from a new, annually dated ice core record from the eastern Ross Sea, named the Roosevelt Island Climate Evolution (RICE) ice core. Comparison of this record with climate reanalysis data for the 1979-2012 interval shows that RICE reliably captures temperature and snow precipitation variability in the region. Trends over the past 2700 years in RICE are shown to be distinct from those in West Antarctica and the western Ross Sea captured by other ice cores. For most of this interval, the eastern Ross Sea was warming (or showing isotopic enrichment for other reasons), with increased snow accumulation and perhaps decreased sea ice concentration. However, West Antarctica cooled and the western Ross Sea showed no significant isotope temperature trend. This pattern here is referred to as the Ross Sea Dipole. Notably, during the Little Ice Age, West Antarctica and the western Ross Sea experienced colder than average temperatures, while the eastern Ross Sea underwent a period of warming or increased isotopic enrichment. From the 17th century onwards, this dipole relationship changed. All three regions show current warming, with snow accumulation declining in West Antarctica and the eastern Ross Sea but increasing in the western Ross Sea. We interpret this pattern as reflecting an increase in sea ice in the eastern Ross Sea with perhaps the establishment of a modern Roosevelt Island polynya as a local moisture source for RICE.

Bereiter, B, Shackleton S, Baggenstos D, Kawamura K, Severinghaus J.  2018.  Mean global ocean temperatures during the last glacial transition. Nature. 553:39-+.   10.1038/nature25152   AbstractWebsite

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.

2017
Seltzer, AM, Buizert C, Baggenstos D, Brook EJ, Ahn J, Yang JW, Severinghaus JP.  2017.  Does delta O-18 of O-2 record meridional shifts in tropical rainfall? Climate of the Past. 13:1323-1338.   10.5194/cp-13-1323-2017   AbstractWebsite

Marine sediments, speleothems, paleo-lake elevations, and ice core methane and delta O-18 of O-2 (delta O-18(atm)) 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 delta O-18(atm) and inferred changes in fractionation by the global terrestrial biosphere (Delta epsilon(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 Delta epsilon(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 delta(OatmO)-O-18 of terrestrial precipitation (the source water for atmospheric O-2 production), we propose a simple mechanism by which Delta epsilon(LAND) tracks the centroid latitude of terrestrial oxygen production. As intense rainfall and oxygen production migrate northward, Delta epsilon(LAND) should decrease due to the underlying meridional gradient in rainfall delta O-18. A southward shift should increase Delta epsilon(LAND). Monsoon intensity also influences delta O-18 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 Delta epsilon(LAND) unambiguously imply a southward shift of tropical rainfall. The exact magnitude of this shift, however, remains under-constrained by Delta epsilon(LAND).

Hamme, RC, Emerson SR, Severinghaus JP, Long MC, Yashayaev I.  2017.  Using noble gas measurements to derive air-sea process information and predict physical gas saturations. Geophysical Research Letters. 44:9901-9909.   10.1002/2017gl075123   AbstractWebsite

Dissolved gas distributions are important because they influence oceanic habitats and Earth's climate, yet competing controls by biology and physics make gas distributions challenging to predict. Bubble-mediated gas exchange, temperature change, and varying atmospheric pressure all push gases away from equilibrium. Here we use new noble gas measurements from the Labrador Sea to demonstrate a technique to quantify physical processes. Our analysis shows that water-mass formation can be represented by a quasi steady state in which bubble fluxes and cooling push gases away from equilibrium balanced by diffusive gas exchange forcing gases toward equilibrium. We quantify the rates of these physical processes from our measurements, allowing direct comparison to gas exchange parameterizations, and predict the physically driven saturation of other gases. This technique produces predictions that reasonably match N-2/Ar observations and demonstrates that physical processes should force SF6 to be approximate to 6% more supersaturated than CFC-11 and CFC-12, impacting ventilation age calculations. Plain Language Summary Gases dissolved in the ocean are important because they influence oceanic habitats and Earth's climate. Physics and biology combine to control the amounts of gases like carbon dioxide, oxygen, and nitrogen in the ocean. Our research seeks to disentangle and quantify the competing effects of physics and biology on dissolved gases. We present very precise measurements of dissolved noble gas concentrations (neon, argon, and krypton) in the Labrador Sea, one of the few places on Earth where the surface and deep ocean communicate with each other. Because noble gases have no biological function, responding only to physical processes in the ocean, we use these measurements to discover the amounts of physical processes that affect gases during the winter at this site, like rapid cooling of the water or bubbles injected by breaking waves. From these amounts of physical processes, we calculate the concentrations of nitrogen and chlorofluorocarbons if only physical processes affected these gases. Our work will allow oceanographers to better estimate the rate that bioavailable nutrients are being removed from the ocean (a process that biologically creates nitrogen gas) and to better determine how the ocean moves from observations of changing chlorofluorocarbons in the ocean.

McConnell, JR, Burke A, Dunbar NW, Kohler P, Thomas JL, Arienzo MM, Chellman NJ, Maselli OJ, Sigl M, Adkins JF, Baggenstos D, Burkhart JF, Brook EJ, Buizert C, Cole-Dai J, Fudge TJ, Knorr G, Graf HF, Grieman MM, Iverson N, McGwire KC, Mulvaney R, Paris G, Rhodes RH, Saltzman ES, Severinghaus JP, Steffensen JP, Taylor KC, Winckler G.  2017.  Synchronous volcanic eruptions and abrupt climate change similar to 17.7 ka plausibly linked by stratospheric ozone depletion. Proceedings of the National Academy of Sciences of the United States of America. 114:10035-10040.   10.1073/pnas.1705595114   AbstractWebsite

Glacial-state greenhouse gas concentrations and Southern Hemisphere climate conditions persisted until similar to 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, similar to 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 similar to 17.7 ka.

Baggenstos, D, Bauska TK, Severinghaus JP, Lee JE, Schaefer H, Buizert C, Brook EJ, Shackleton S, Petrenko VV.  2017.  Atmospheric gas records from Taylor Glacier, Antarctica, reveal ancient ice with ages spanning the entire last glacial cycle. Climate of the Past. 13:943-958.   10.5194/cp-13-943-2017   AbstractWebsite

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 delta O-18 of O-2 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.

Orsi, AJ, Kawamura K, Masson-Delmotte V, Fettweis X, Box JE, Dahl-Jensen D, Clow GD, Landais A, Severinghaus JP.  2017.  The recent warming trend in North Greenland. Geophysical Research Letters. 44:6235-6243.   10.1002/2016gl072212   AbstractWebsite

The Arctic is among the fastest warming regions on Earth, but it is also one with limited spatial coverage of multidecadal instrumental surface air temperature measurements. Consequently, atmospheric reanalyses are relatively unconstrained in this region, resulting in a large spread of estimated 30 year recent warming trends, which limits their use to investigate the mechanisms responsible for this trend. Here we present a surface temperature reconstruction over 1982-2011 at NEEM (North Greenland Eemian Ice Drilling Project, 51 degrees W, 77 degrees N), in North Greenland, based on the inversion of borehole temperature and inert gas isotope data. We find that NEEM has warmed by 2.7 +/- 0.33 degrees C over the past 30 years, from the long-term 1900-1970 average of -28.55 +/- 0.29 degrees C. The warming trend is principally caused by an increase in downward longwave heat flux. Atmospheric reanalyses underestimate this trend by 17%, underlining the need for more in situ observations to validate reanalyses.

Seltzer, AM, Severinghaus JP, Andraski BJ, Stonestrom DA.  2017.  Steady state fractionation of heavy noble gas isotopes in a deep unsaturated zone. Water Resources Research. 53:2716-2732.   10.1002/2016WR019655   AbstractWebsite

To explore steady state fractionation processes in the unsaturated zone (UZ), we measured argon, krypton, and xenon isotope ratios throughout a ∼110 m deep UZ at the United States Geological Survey (USGS) Amargosa Desert Research Site (ADRS) in Nevada, USA. Prior work has suggested that gravitational settling should create a nearly linear increase in heavy-to-light isotope ratios toward the bottom of stagnant air columns in porous media. Our high-precision measurements revealed a binary mixture between (1) expected steady state isotopic compositions and (2) unfractionated atmospheric air. We hypothesize that the presence of an unsealed pipe connecting the surface to the water table allowed for direct inflow of surface air in response to extensive UZ gas sampling prior to our first (2015) measurements. Observed isotopic resettling in deep UZ samples collected a year later, after sealing the pipe, supports this interpretation. Data and modeling each suggest that the strong influence of gravitational settling and weaker influences of thermal diffusion and fluxes of CO2 and water vapor accurately describe steady state isotopic fractionation of argon, krypton, and xenon within the UZ. The data confirm that heavy noble gas isotopes are sensitive indicators of UZ depth. Based on this finding, we outline a potential inverse approach to quantify past water table depths from noble gas isotope measurements in paleogroundwater, after accounting for fractionation during dissolution of UZ air and bubbles.

2016
Cuffey, KM, Clow GD, Steig EJ, Buizert C, Fudge TJ, Koutnik M, Waddington ED, Alley RB, Severinghaus JP.  2016.  Deglacial temperature history of West Antarctica. Proceedings of the National Academy of Sciences of the United States of America. 113:14249-14254.   10.1073/pnas.1609132113   AbstractWebsite

The most recent glacial to interglacial transition constitutes a remarkable natural experiment for learning how Earth's climate responds to various forcings, including a rise in atmospheric CO2. This transition has left a direct thermal remnant in the polar ice sheets, where the exceptional purity and continual accumulation of ice permit analyses not possible in other settings. For Antarctica, the deglacial warming has previously been constrained only by the water isotopic composition in ice cores, without an absolute thermometric assessment of the isotopes' sensitivity to temperature. To overcome this limitation, we measured temperatures in a deep borehole and analyzed them together with ice-core data to reconstruct the surface temperature history of West Antarctica. The deglacial warming was 11.3 +/- 1.8 degrees C, approximately two to three times the global average, in agreement with theoretical expectations for Antarctic amplification of planetary temperature changes. Consistent with evidence from glacier retreat in Southern Hemisphere mountain ranges, the Antarctic warming was mostly completed by 15 kyBP, several millennia earlier than in the Northern Hemisphere. These results constrain the role of variable oceanic heat transport between hemispheres during deglaciation and quantitatively bound the direct influence of global climate forcings on Antarctic temperature. Although climate models perform well on average in this context, some recent syntheses of deglacial climate history have underestimated Antarctic warming and the models with lowest sensitivity can be discounted.

Goodge, JW, Severinghaus JP.  2016.  Rapid Access Ice Drill: a new tool for exploration of the deep Antarctic ice sheets and subglacial geology. Journal of Glaciology. 62:1049-1064.   10.1017/jog.2016.97   AbstractWebsite

A new Rapid Access Ice Drill (RAID) will penetrate the Antarctic ice sheets in order to create borehole observatories and take cores in deep ice, the glacial bed and bedrock below. RAID is a mobile drilling system to make multiple long, narrow boreholes in a single field season in Antarctica. RAID is based on a mineral exploration-type rotary rock-coring system using threaded drill pipe to cut through ice using reverse circulation of a non-freezing fluid for pressure-compensation, maintenance of temperature and removal of ice cuttings. Near the bottom of the ice sheet, a wireline latching assembly will enable rapid coring of ice, the glacial bed and bedrock below. Once complete, boreholes will be kept open with fluid, capped and available for future down-hole measurement of temperature gradient, heat flow, ice chronology and ice deformation. RAID is designed to penetrate up to 3300 m of ice and take cores in <200 hours, allowing completion of a borehole and coring in similar to 10 d at each site. Together, the rapid drilling capability and mobility of the system, along with ice-penetrating imaging methods, will provide a unique 3-D picture of interior and subglacial features of the Antarctic ice sheets.

Buizert, C, Severinghaus JP.  2016.  Dispersion in deep polar firn driven by synoptic-scale surface pressure variability. Cryosphere. 10:2099-2111.   10.5194/tc-10-2099-20160   AbstractWebsite

Commonly, three mechanisms of firn air transport are distinguished: molecular diffusion, advection, and near-surface convective mixing. Here we identify and describe a fourth mechanism, namely dispersion driven by synoptic-scale surface pressure variability (or barometric pumping). We use published gas chromatography experiments on firn samples to derive the along-flow dispersivity of firn, and combine this dispersivity with a dynamical air pressure propagation model forced by surface air pressure time series to estimate the magnitude of dispersive mixing in the firn. We show that dispersion dominates mixing within the firn lock-in zone. Trace gas concentrations measured in firn air samples from various polar sites confirm that dispersive mixing occurs. Including dispersive mixing in a firn air transport model suggests that our theoretical estimates have the correct order of magnitude, yet may overestimate the true dispersion. We further show that strong barometric pumping, such as at the Law Dome site, may reduce the gravitational enrichment of delta N-15-N-2 and other tracers below gravitational equilibrium, questioning the traditional definition of the lock-in depth as the depth where delta N-15 enrichment ceases. Last, we propose that Kr-86 excess may act as a proxy for past synoptic activity (or paleo-storminess) at the site.

Bauska, TK, Baggenstos D, Brook EJ, Mix AC, Marcott SA, Petrenko VV, Schaefer H, Severinghaus JP, Lee JE.  2016.  Carbon isotopes characterize rapid changes in atmospheric carbon dioxide during the last deglaciation. Proceedings of the National Academy of Sciences of the United States of America. 113:3465-3470.   10.1073/pnas.1513868113   AbstractWebsite

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 (delta C-13-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 delta C-13-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 delta C-13-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 delta C-13-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 Bolling (14.6-14.3 ka) and Holocene (11.6-11.4 ka) intervals are associated with small changes in delta C-13-CO2, suggesting a combination of sources that included rising surface ocean temperature.

Petrenko, VV, Severinghaus JP, Schaefer H, Smith AM, Kuhl T, Baggenstos D, Hua Q, Brook EJ, Rose P, Kulin R, Bauska T, Harth C, Buizert C, Orsi A, Emanuele G, Lee JE, Brailsford G, Keeling R, Weiss RF.  2016.  Measurements of 14C in ancient ice from Taylor Glacier, Antarctica constrain in situ cosmogenic 14CH4 and 14CO production rates. Geochimica et Cosmochimica Acta. 177:62-77.   10.1016/j.gca.2016.01.004   Abstract

Carbon-14 (14C) is incorporated into glacial ice by trapping of atmospheric gases as well as direct near-surface in situ cosmogenic production. 14C of trapped methane (14CH4) is a powerful tracer for past CH4 emissions from “old” carbon sources such as permafrost and marine CH4 clathrates. 14C in trapped carbon dioxide (14CO2) can be used for absolute dating of ice cores. In situ produced cosmogenic 14C in carbon monoxide (14CO) can potentially be used to reconstruct the past cosmic ray flux and past solar activity. Unfortunately, the trapped atmospheric and in situ cosmogenic components of 14C in glacial ice are difficult to disentangle and a thorough understanding of the in situ cosmogenic component is needed in order to extract useful information from ice core 14C. We analyzed very large (≈1000 kg) ice samples in the 2.26–19.53 m depth range from the ablation zone of Taylor Glacier, Antarctica, to study in situ cosmogenic production of 14CH4 and 14CO. All sampled ice is >50 ka in age, allowing for the assumption that most of the measured 14C originates from recent in situ cosmogenic production as ancient ice is brought to the surface via ablation. Our results place the first constraints on cosmogenic 14CH4 production rates and improve on prior estimates of 14CO production rates in ice. We find a constant 14CH4/14CO production ratio (0.0076 ± 0.0003) for samples deeper than 3 m, which allows the use of 14CO for correcting the 14CH4 signals for the in situ cosmogenic component. Our results also provide the first unambiguous confirmation of 14C production by fast muons in a natural setting (ice or rock) and suggest that the 14C production rates in ice commonly used in the literature may be too high.

2015
Orsi, AJ, Kawamura K, Fegyveresi JM, Headly MA, Alley RB, Severinghaus JP.  2015.  Differentiating bubble-free layers from melt layers in ice cores using noble gases. Journal of Glaciology. 61:585-594.   10.3189/2015JoG14J237   AbstractWebsite

Melt layers are clear indicators of extreme summer warmth on polar ice caps. The visual identification of refrozen meltwater as clear bubble-free layers cannot be used to study some past warm periods, because, in deeper ice, bubbles are lost to clathrate formation. We present here a reliable method to detect melt events, based on the analysis of Kr/Ar and Xe/Ar ratios in ice cores, and apply it to the detection of melt in clathrate ice from the Eemian at NEEM, Greenland. Additionally, melt layers in ice cores can compromise the integrity of the gas record by dissolving soluble gases, or by altering gas transport in the firn, which affects the gas chronology. We find that the easily visible 1 mm thick bubble-free layers in the WAIS Divide ice core do not contain sufficient melt to alter the gas composition in the core, and do not cause artifacts or discontinuities in the gas chronology. The presence of these layers during winter, and the absence of anomalies in soluble gases, suggests that these layers can be formed by processes other than refreezing of meltwater. Consequently, the absence of bubbles in thin crusts is not in itself proof of a melt event.

Rhodes, RH, Brook EJ, Chiang JCH, Blunier T, Maselli OJ, McConnell JR, Romanini D, Severinghaus JP.  2015.  Enhanced tropical methane production in response to iceberg discharge in the North Atlantic. Science. 348:1016-1019.   10.1126/science.1262005   AbstractWebsite

The causal mechanisms responsible for the abrupt climate changes of the Last Glacial Period remain unclear. One major difficulty is dating ice-rafted debris deposits associated with Heinrich events: Extensive iceberg influxes into the North Atlantic Ocean linked to global impacts on climate and biogeochemistry. In a new ice core record of atmospheric methane with ultrahigh temporal resolution, we find abrupt methane increases within Heinrich stadials 1, 2, 4, and 5 that, uniquely, have no counterparts in Greenland temperature proxies. Using a heuristic model of tropical rainfall distribution, we propose that Hudson Strait Heinrich events caused rainfall intensification over Southern Hemisphere land areas, thereby producing excess methane in tropical wetlands. Our findings suggest that the climatic impacts of Heinrich events persisted for 740 to 1520 years.

Buizert, C, Adrian B, Ahn J, Albert M, Alley RB, Baggenstos D, Bauska TK, Bay RC, Bencivengo BB, Bentley CR, Brook EJ, Chellman NJ, Clow GD, Cole-Dai J, Conway H, Cravens E, Cuffey KM, Dunbar NW, Edwards JS, Fegyveresi JM, Ferris DG, Fitzpatrick JJ, Fudge TJ, Gibson CJ, Gkinis V, Goetz JJ, Gregory S, Hargreaves GM, Iverson N, Johnson JA, Jones TR, Kalk ML, Kippenhan MJ, Koffman BG, Kreutz K, Kuhl TW, Lebar DA, Lee JE, Marcott SA, Markle BR, Maselli OJ, McConnell JR, McGwire KC, Mitchell LE, Mortensen NB, Neff PD, Nishiizumi K, Nunn RM, Orsi AJ, Pasteris DR, Pedro JB, Pettit EC, Price PB, Priscu JC, Rhodes RH, Rosen JL, Schauer AJ, Schoenemann SW, Sendelbach PJ, Severinghaus JP, Shturmakov AJ, Sigl M, Slawny KR, Souney JM, Sowers TA, Spencer MK, Steig EJ, Taylor KC, Twickler MS, Vaughn BH, Voigt DE, Waddington ED, Welten KC, Wendricks AW, White JWC, Winstrup M, Wong GJ, Woodruff TE, Members WDP.  2015.  Precise interpolar phasing of abrupt climate change during the last ice age. Nature. 520:661-U169.   10.1038/nature14401   AbstractWebsite

The last glacial period exhibited abrupt Dansgaard-Oeschger climatic oscillations, evidence of which is preserved in a variety of Northern Hemisphere palaeodimate 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(4-6). 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 sigma a) for DansgaardOeschger events, including the Bolling 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.

Mitchell, LE, Buizert C, Brook EJ, Breton DJ, Fegyveresi J, Baggenstos D, Orsi A, Severinghaus J, Alley RB, Albert M, Rhodes RH, McConnell JR, Sigl M, Maselli O, Gregory S, Ahn J.  2015.  Observing and modeling the influence of layering on bubble trapping in polar firn. Journal of Geophysical Research-Atmospheres. 120:2558-2574.   10.1002/2014jd022766   AbstractWebsite

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.

Buizert, C, Cuffey KM, Severinghaus JP, Baggenstos D, Fudge TJ, Steig EJ, Markle BR, Winstrup M, Rhodes RH, Brook EJ, Sowers TA, Clow GD, Cheng H, Edwards RL, Sigl M, McConnell JR, Taylor KC.  2015.  The WAIS Divide deep ice core WD2014 chronology - Part 1: Methane synchronization (68-31 kaBP) and the gas age-ice age difference. Climate of the Past. 11:153-173.   10.5194/cp-11-153-2015   AbstractWebsite

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 similar to 68 ka at unprecedented temporal resolution. The upper 2850m (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 (Delta age) using a combination of firn densification modeling, ice-flow modeling, and a data set of delta N-15-N-2, a proxy for past firn column thickness. The largest Delta 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 Delta 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".

2014
Seierstad, IK, Abbott PM, Bigler M, Blunier T, Bourne AJ, Brook E, Buchardt SL, Buizert C, Clausen HB, Cook E, Dahl-Jensen D, Davies SM, Guillevic M, Johnsen SJ, Pedersen DS, Popp TJ, Rasmussen SO, Severinghaus JP, Svensson A, Vinther BM.  2014.  Consistently dated records from the Greenland GRIP, GISP2 and NGRIP ice cores for the past 104 ka reveal regional millennial-scale delta O-18 gradients with possible Heinrich event imprint. Quaternary Science Reviews. 106:29-46.   10.1016/j.quascirev.2014.10.032   AbstractWebsite

We present a synchronization of the NGRIP, GRIP and GISP2 ice cores onto a master chronology extending back to 104 ka before present, providing a consistent chronological framework for these three Greenland records. The synchronization aligns distinct peaks in volcanic proxy records and other impurity records (chemo-stratigraphic matching) and assumes that these layers of elevated impurity content represent the same, instantaneous event in the past at all three sites. More than 900 marker horizons between the three cores have been identified and our matching is independently confirmed by 24 new and previously identified volcanic ash (tephra) tie-points. Using the reference horizons, we transfer the widely used Greenland ice-core chronology, GICC05modelext, to the two Summit cores, GRIP and GISP2. Furthermore, we provide gas chronologies for the Summit cores that are consistent with the GICC05modelext timescale by utilizing both existing and new gas data (CH4 concentration and delta N-15 of N-2). We infer that the accumulation contrast between the stadial and interstadial phases of the glacial period was -10% greater at Summit compared to at NGRIP. The delta O-18 temperature-proxy records from NGRIP, GRIP, and GISP2 are generally very similar and display synchronous behaviour at climate transitions. The 1180 differences between Summit and NGRIP, however, changed slowly over the Last Glacial Interglacial cycle and also underwent abrupt millennial-to-centennial-scale variations. We suggest that this observed latitudinal delta O-18 gradient in Greenland during the glacial period is the result of 1) relatively higher degree of precipitation with a Pacific signature at NGRIP, 2) increased summer bias in precipitation at Summit, and 3) enhanced Rayleigh distillation due to an increased source-to-site distance and a potentially larger source-to-site temperature gradient. We propose that these processes are governed by changes in the North American Ice Sheet (NAIS) volume and North Atlantic sea-ice extent and/or sea-surface temperatures (SST) on orbital timescales, and that changing sea-ice extent and SSTs are the driving mechanisms on shorter timescales. Finally, we observe that maxima in the Summit NGRIP delta O-18 difference are roughly coincident with prominent Heinrich events. This suggests that the climatic reorganization that takes place during stadials with Heinrich events, possibly driven by a southward expansion of sea ice and low SSTs in the North Atlantic, are recorded in the ice-core records. (C) 2014 Elsevier Ltd. All rights reserved.

Schilt, A, Brook EJ, Bauska TK, Baggenstos D, Fischer H, Joos F, Petrenko VV, Schaefer H, Schmitt J, Severinghaus JP, Spahni R, Stocker TF.  2014.  Isotopic constraints on marine and terrestrial N2O emissions during the last deglaciation. Nature. 516:234-+.   10.1038/nature13971   AbstractWebsite

Nitrous oxide (N2O) is an important greenhouse gas and ozone-depleting substance that has anthropogenic as well as natural marine and terrestrial sources(1). The tropospheric N2O concentrations have varied substantially in the past in concert with changing climate on glacial-interglacial and millennial timescales(2-8). 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(4,9,10). 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(11).

Marcott, SA, Bauska TK, Buizert C, Steig EJ, Rosen JL, Cuffey KM, Fudge TJ, Severinghaus JP, Ahn J, Kalk ML, McConnell JR, Sowers T, Taylor KC, White JWC, Brook EJ.  2014.  Centennial-scale changes in the global carbon cycle during the last deglaciation. Nature. 514:616-+.   10.1038/nature13799   AbstractWebsite

Global climate and the concentration of atmospheric carbon dioxide (CO2) are correlated over recent glacial cycles(1,2). The combination of processes responsible for a rise in atmospheric CO2 at the last glacial termination(1,3) (23,000 to 9,000 years ago), however, remains uncertain(1-3). Establishing the timing and rate of CO2 changes in the past provides critical insight into the mechanisms that influence the carbon cycle and helps put present and future anthropogenic emissions in context. Here we present CO2 and methane (CH4) records of the last deglaciation from a new high-accumulation West Antarctic ice core with unprecedented temporal resolution and precise chronology. We show that although low-frequency CO2 variations parallel changes in Antarctic temperature, abrupt CO2 changes occur that have a clear relationship with abrupt climate changes in the Northern Hemisphere. A significant proportion of the direct radiative forcing associated with the rise in atmospheric CO2 occurred in three sudden steps, each of 10 to 15 parts per million. Every step took place in less than two centuries and was followed by no notable change in atmospheric CO2 for about 1,000 to 1,500 years. Slow, millennial-scale ventilation of Southern Ocean CO2-rich, deep-ocean water masses is thought to have been fundamental to the rise in atmospheric CO2 associated with the glacial termination(4), given the strong covariance of CO2 levels and Antarctic temperatures(5). Our data establish a contribution from an abrupt, centennial-scale mode of CO2 variability that is not directly related to Antarctic temperature. We suggest that processes operating on centennial timescales, probably involving the Atlantic meridional overturning circulation, seem to be influencing global carbon-cycle dynamics and are at present not widely considered in Earth system models.

Buizert, C, Gkinis V, Severinghaus JP, He F, Lecavalier BS, Kindler P, Leuenberger M, Carlson AE, Vinther B, Masson-Delmotte V, White JWC, Liu ZY, Otto-Bliesner B, Brook EJ.  2014.  Greenland temperature response to climate forcing during the last deglaciation. Science. 345:1177-1180.   10.1126/science.1254961   AbstractWebsite

Greenland ice core water isotopic composition (delta O-18) provides detailed evidence for abrupt climate changes but is by itself insufficient for quantitative reconstruction of past temperatures and their spatial patterns. We investigate Greenland temperature evolution during the last deglaciation using independent reconstructions from three ice cores and simulations with a coupled ocean-atmosphere climate model. Contrary to the traditional delta O-18 interpretation, the Younger Dryas period was 4.5 degrees +/- 2 degrees C warmer than the Oldest Dryas, due to increased carbon dioxide forcing and summer insolation. The magnitude of abrupt temperature changes is larger in central Greenland (9 degrees to 14 degrees C) than in the northwest (5 degrees to 9 degrees C), fingerprinting a North Atlantic origin. Simulated changes in temperature seasonality closely track changes in the Atlantic overturning strength and support the hypothesis that abrupt climate change is mostly a winter phenomenon.