Publications

Export 92 results:
Sort by: Author Title Type [ Year  (Desc)]
2017
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.

Rosen, JL, Brook EJ, Severinghaus JP, Blunier T, Mitchell LE, Lee JE, Edwards JS, Gkinis V.  2014.  An ice core record of near-synchronous global climate changes at the Bolling transition. Nature Geoscience. 7:459-463.   10.1038/ngeo2147   AbstractWebsite

The abrupt warming that initiated the Bolling-Allerod interstadial was the penultimate warming in a series of climate variations known as Dansgaard-Oeschger events. Despite the clear expression of this transition in numerous palaeoclimate records, the relative timing of climate shifts in different regions of the world and their causes are subject to debate. Here we explore the phasing of global climate change at the onset of the Bolling-Allerod using air preserved in bubbles in the North Greenland Eemian ice core. Specifically, we measured methane concentrations, which act as a proxy for low-latitude climate, and the N-15/N-14 ratio of N-2, which reflects Greenland surface temperature, over the same interval of time. We use an atmospheric box model and a firn air model to account for potential uncertainties in the data, and find that changes in Greenland temperature and atmospheric methane emissions at the Bolling onset occurred essentially synchronously, with temperature leading by 4.5(-24)(+21) years. We cannot exclude the possibility that tropical climate could iag changing methane concentrations by up to several decades, if the initial methane rise came from boreal sources alone. However, because even boreal methane-producing regions lie far from Greenland, we conclude that the mechanism that drove abrupt change at this time must be capable of rapidly transmitting climate changes across the globe.

Orsi, AJ, Cornuelle BD, Severinghaus JP.  2014.  Magnitude and temporal evolution of Dansgaard-Oeschger event 8 abrupt temperature change inferred from nitrogen and argon isotopes in GISP2 ice using a new least-squares inversion. Earth and Planetary Science Letters. 395:81-90.   10.1016/j.epsl.2014.03.030   AbstractWebsite

Polar temperature is often inferred from water isotopes in ice cores. However, non-temperature effects on 3180 are important during the abrupt events of the last glacial period, such as changes in the seasonality of precipitation, the northward movement of the storm track, and the increase in accumulation. These effects complicate the interpretation of 8180 as a temperature proxy. Here, we present an independent surface temperature reconstruction, which allows us to test the relationship between delta O-18(ice) and temperature, during Dansgaard-Oeschger event 8, 38.2 thousand yrs ago using new delta N-15 and delta Ar-40 data from the GISP2 ice core in Greenland. This temperature reconstruction relies on a new inversion of inert gas isotope data using generalized least-squares, and includes a robust uncertainty estimation. We find that both temperature and delta O-18 increased in two steps of 20 and 140 yrs, with an overall amplitude of 11.80 +/- 1.8 degrees C between the stadial and interstadial centennial-mean temperature. The coefficient alpha = d delta O-18/dT changes with each time-segment, which shows that non-temperature sources of fractionation have a significant contribution to the delta O-18 signal. When measured on century-averaged values, we find that alpha = d delta O-18/dT = 0.32 +/- 0.06%(0)/degrees C, which is similar to the glacial/Holocene value of 0.328%(o)/degrees C. (C) 2014 Elsevier B.V. All rights reserved.

Buizert, C, Baggenstos D, Jiang W, Purtschert R, Petrenko VV, Lu ZT, Muller P, Kuhl T, Lee J, Severinghaus JP, Brook EJ.  2014.  Radiometric Kr-81 dating identifies 120,000-year-old ice at Taylor Glacier, Antarctica. Proceedings of the National Academy of Sciences of the United States of America. 111:6876-6881.   10.1073/pnas.1320329111   AbstractWebsite

We present successful Kr-81-Kr radiometric dating of ancient polar ice. Krypton was extracted from the air bubbles in four similar to 350-kg polar ice samples from Taylor Glacier in the McMurdo Dry Valleys, Antarctica, and dated using Atom Trap Trace Analysis (ATTA). The Kr-81 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) Kr-85 and Ar-39 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 Kr-81 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 Kr-81 analysis requires a 40-80-kg ice sample; as sample requirements continue to decrease, Kr-81 dating of ice cores is a future possibility.

2013
Arnold, T, Harth CM, Mühle J, Manning AJ, Salameh PK, Kim J, Ivy DJ, Steele PL, Petrenko VV, Severinghaus JP, Baggenstos D, Weiss RF.  2013.  Nitrogen trifluoride global emissions estimated from updated atmospheric measurements. Proceedings of the National Academy of Sciences.   10.1073/pnas.1212346110   AbstractWebsite

Nitrogen trifluoride (NF3) 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 NF3 in 2011 were 1.18 ± 0.21 Gg⋅y−1, or ∼20 Tg CO2-eq⋅y−1 (carbon dioxide equivalent emissions based on a 100-y global warming potential of 16,600 for NF3). 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 CO2 equivalents, current NF3 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 NF3 over hexafluoroethane (C2F6) in electronics manufacture is significant—emissions of between 53 and 220 Tg CO2-eq⋅y−1 were avoided during 2011. Despite these savings, total NF3 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 NF3 emissions reduction strategies to keep pace with its increasing use and to slow its rising contribution to anthropogenic climate forcing.

Petrenko, VV, Martinerie P, Novelli P, Etheridge DM, Levin I, Wang Z, Blunier T, Chappellaz J, Kaiser J, Lang P, Steele LP, Hammer S, Mak J, Langenfelds RL, Schwander J, Severinghaus JP, Witrant E, Petron G, Battle MO, Forster G, Sturges WT, Lamarque JF, Steffen K, White JWC.  2013.  A 60 yr record of atmospheric carbon monoxide reconstructed from Greenland firn air. Atmospheric Chemistry and Physics. 13:7567-7585.   10.5194/acp-13-7567-2013   AbstractWebsite

We present the first reconstruction of the Northern Hemisphere (NH) high latitude atmospheric carbon monoxide (CO) mole fraction from Greenland firn air. Firn air samples were collected at three deep ice core sites in Greenland (NGRIP in 2001, Summit in 2006 and NEEM in 2008). CO records from the three sites agree well with each other as well as with recent atmospheric measurements, indicating that CO is well preserved in the firn at these sites. CO atmospheric history was reconstructed back to the year 1950 from the measurements using a combination of two forward models of gas transport in firn and an inverse model. The reconstructed history suggests that Arctic CO in 1950 was 140-150 nmol mol(-1), which is higher than today's values. CO mole fractions rose by 10-15 nmol mol(-1) from 1950 to the 1970s and peaked in the 1970s or early 1980s, followed by a approximate to 30 nmol mol(-1) decline to today's levels. We compare the CO history with the atmospheric histories of methane, light hydrocarbons, molecular hydrogen, CO stable isotopes and hydroxyl radicals (OH), as well as with published CO emission inventories and results of a historical run from a chemistry-transport model. We find that the reconstructed Greenland CO history cannot be reconciled with available emission inventories unless unrealistically large changes in OH are assumed. We argue that the available CO emission inventories strongly underestimate historical NH emissions, and fail to capture the emission decline starting in the late 1970s, which was most likely due to reduced emissions from road transportation in North America and Europe.

Fudge, TJ, Steig EJ, Markle BR, Schoenemann SW, Ding QH, Taylor KC, McConnell JR, Brook EJ, Sowers T, White JWC, Alley RB, Cheng H, Clow GD, Cole-Dai J, Conway H, Cuffey KM, Edwards JS, Edwards RL, Edwards R, Fegyveresi JM, Ferris D, Fitzpatrick JJ, Johnson J, Hargreaves G, Lee JE, Maselli OJ, Mason W, McGwire KC, Mitchell LE, Mortensen N, Neff P, Orsi AJ, Popp TJ, Schauer AJ, Severinghaus JP, Sigl M, Spencer MK, Vaughn BH, Voigt DE, Waddington ED, Wang XF, Wong GJ, Members WDP.  2013.  Onset of deglacial warming in West Antarctica driven by local orbital forcing. Nature. 500:440-+.   10.1038/nature12376   AbstractWebsite

The cause of warming in the Southern Hemisphere during the most recent deglaciation remains a matter of debate(1,2). Hypotheses for a Northern Hemisphere trigger, through oceanic redistributions of heat, are based in part on the abrupt onset of warming seen in East Antarctic ice cores and dated to 18,000 years ago, which is several thousand years after high-latitude Northern Hemisphere summer insolation intensity began increasing from its minimum, approximately 24,000 years ago(3,4). An alternative explanation is that local solar insolation changes cause the Southern Hemisphere to warm independently(2,5). Here we present results from a new, annually resolved ice-core record from West Antarctica that reconciles these two views. The records show that 18,000 years ago snow accumulation in West Antarctica began increasing, coincident with increasing carbon dioxide concentrations, warming in East Antarctica and cooling in the Northern Hemisphere(6) associated with an abrupt decrease in Atlantic meridional overturning circulation(7). However, significant warming in West Antarctica began at least 2,000 years earlier. Circum-Antarctic sea-ice decline, driven by increasing local insolation, is the likely cause of this warming. The marine-influenced West Antarctic records suggest a more active role for the Southern Ocean in the onset of deglaciation than is inferred from ice cores in the East Antarctic interior, which are largely isolated from sea-ice changes.

Petrenko, VV, Severinghaus JP, Smith AM, Riedel K, Baggenstos D, Harth C, Orsi A, Hua Q, Franz P, Takeshita Y, Brailsford GW, Weiss RF, Buizert C, Dickson A, Schaefer H.  2013.  High-precision C-14 measurements demonstrate production of in situ cosmogenic (CH4)-C-14 and rapid loss of in situ cosmogenic (CO)-C-14 in shallow Greenland firn. Earth and Planetary Science Letters. 365:190-197.   10.1016/j.epsl.2013.01.032   AbstractWebsite

Measurements of radiocarbon (C-14) in carbon dioxide (CO2), methane (CH4) and carbon monoxide (CO) from glacial ice are potentially useful for absolute dating of ice cores, studies of the past atmospheric CH4 budget and for reconstructing the past cosmic ray flux and solar activity. Interpretation of C-14 signals in ice is complicated by the fact that the two major C-14 components-trapped atmospheric and in situ cosmogenic-are present in a combined form, as well as by a very limited understanding of the in situ component. This study measured (CH4)-C-14 and (CO)-C-14 content in glacial firn with unprecedented precision to advance understanding of the in situ C-14 component. (CH4)-C-14 and (CO)-C-14 were melt-extracted on site at Summit, Greenland from three very large (similar to 1000 kg each) replicate samples of firn that spanned a depth range of 3.6-5.6 m. Non-cosmogenic C-14 contributions were carefully characterized through simulated extractions and a suite of supporting measurements. In situ cosmogenic (CO)-C-14 was quantified to better than +/- 0.6 molecules g(-1) ice, improving on the precision of the best prior ice (CO)-C-14 measurements by an order of magnitude. The (CO)-C-14 measurements indicate that most (>99%) of the in situ cosmogenic C-14 is rapidly lost from shallow Summit firn to the atmosphere. Despite this rapid C-14 loss, our measurements successfully quantified (CH4)-C-14 in the retained fraction of cosmogenic C-14 (to +/- 0.01 molecules g(-1) ice or better), and demonstrate for the first time that a significant amount of (CH4)-C-14 is produced by cosmic rays in natural ice. This conclusion increases the confidence in the results of an earlier study that used measurements of (CH4)-C-14 in glacial ice to show that wetlands were the likely main driver of the large and rapid atmospheric CH4 increase approximately 1 1.6 kyr ago. (C) 2013 Elsevier B.V. All rights reserved.

Rasmussen, SO, Abbott PM, Blunier T, Bourne AJ, Brook E, Buchardt SL, Buizert C, Chappellaz J, Clausen HB, Cook E, Dahl-Jensen D, Davies SM, Guillevic M, Kipfstuhl S, Laepple T, Seierstad IK, Severinghaus JP, Steffensen JP, Stowasser C, Svensson A, Vallelonga P, Vinther BM, Wilhelms F, Winstrup M.  2013.  A first chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core. Climate of the Past. 9:2713-2730.   10.5194/cp-9-2713-2013   AbstractWebsite

A stratigraphy-based chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core has been derived by transferring the annual layer counted Greenland Ice Core Chronology 2005 (GICC05) and its model extension (GICC05modelext) from the NGRIP core to the NEEM core using 787 match points of mainly volcanic origin identified in the electrical conductivity measurement (ECM) and dielectrical profiling (DEP) records. Tephra horizons found in both the NEEM and NGRIP ice cores are used to test the matching based on ECM and DEP and provide five additional horizons used for the timescale transfer. A thinning function reflecting the accumulated strain along the core has been determined using a Dansgaard-Johnsen flow model and an isotope-dependent accumulation rate parameterization. Flow parameters are determined from Monte Carlo analysis constrained by the observed depth-age horizons. In order to construct a chronology for the gas phase, the ice age-gas age difference (Delta age) has been reconstructed using a coupled firn densification-heat diffusion model. Temperature and accumulation inputs to the Delta age model, initially derived from the water isotope proxies, have been adjusted to optimize the fit to timing constraints from delta N-15 of nitrogen and high-resolution methane data during the abrupt onset of Greenland interstadials. The ice and gas chronologies and the corresponding thinning function represent the first chronology for the NEEM core, named GICC05modelext-NEEM-1. Based on both the flow and firn modelling results, the accumulation history for the NEEM site has been reconstructed. Together, the timescale and accumulation reconstruction provide the necessary basis for further analysis of the records from NEEM.

Kawamura, K, Severinghaus JP, Albert MR, Courville ZR, Fahnestock MA, Scambos T, Shields E, Shuman CA.  2013.  Kinetic fractionation of gases by deep air convection in polar firn. Atmospheric Chemistry and Physics. 13:11141-11155.   10.5194/acp-13-11141-2013   AbstractWebsite

A previously unrecognized type of gas fractionation occurs in firn air columns subjected to intense convection. It is a form of kinetic fractionation that depends on the fact that different gases have different molecular diffusivities. Convective mixing continually disturbs diffusive equilibrium, and gases diffuse back toward diffusive equilibrium under the influence of gravity and thermal gradients. In near-surface firn where convection and diffusion compete as gas transport mechanisms, slow-diffusing gases such as krypton (Kr) and xenon (Xe) are more heavily impacted by convection than fast diffusing gases such as nitrogen (N-2) and argon (Ar), and the signals are preserved in deep firn and ice. We show a simple theory that predicts this kinetic effect, and the theory is confirmed by observations using a newly-developed Kr and Xe stable isotope system in air samples from the Megadunes field site on the East Antarctic plateau. Numerical simulations confirm the effect's magnitude at this site. A main purpose of this work is to support the development of a proxy indicator of past convection in firn, for use in ice-core gas records. To this aim, we also show with the simulations that the magnitude of the kinetic effect is fairly insensitive to the exact profile of convective strength, if the overall thickness of the convective zone is kept constant. These results suggest that it may be feasible to test for the existence of an extremely deep (similar to 30-40 m) convective zone, which has been hypothesized for glacial maxima, by future ice-core measurements.

Fischer, H, Severinghaus J, Brook E, Wolff E, Albert M, Alemany O, Arthern R, Bentley C, Blankenship D, Chappellaz J, Creyts T, Dahl-Jensen D, Dinn M, Frezzotti M, Fujita S, Gallee H, Hindmarsh R, Hudspeth D, Jugie G, Kawamura K, Lipenkov V, Miller H, Mulvaney R, Parrenin F, Pattyn F, Ritz C, Schwander J, Steinhage D, van Ommen T, Wilhelms F.  2013.  Where to find 1.5 million yr old ice for the IPICS "Oldest-Ice" ice core. Climate of the Past. 9:2489-2505.   10.5194/cp-9-2489-2013   AbstractWebsite

The recovery of a 1.5 million yr long ice core from Antarctica represents a keystone of our understanding of Quaternary climate, the progression of glaciation over this time period and the role of greenhouse gas cycles in this progression. Here we tackle the question of where such ice may still be found in the Antarctic ice sheet. We can show that such old ice is most likely to exist in the plateau area of the East Antarctic ice sheet (EAIS) without stratigraphic disturbance and should be able to be recovered after careful pre-site selection studies. Based on a simple ice and heat flow model and glaciological observations, we conclude that positions in the vicinity of major domes and saddle position on the East Antarctic Plateau will most likely have such old ice in store and represent the best study areas for dedicated reconnaissance studies in the near future. In contrast to previous ice core drill site selections, however, we strongly suggest significantly reduced ice thickness to avoid bottom melting. For example for the geothermal heat flux and accumulation conditions at Dome C, an ice thickness lower than but close to about 2500m would be required to find 1.5 Myr old ice (i.e., more than 700m less than at the current EPICA Dome C drill site). Within this constraint, the resolution of an Oldest-Ice record and the distance of such old ice to the bedrock should be maximized to avoid ice flow disturbances, for example, by finding locations with minimum geothermal heat flux. As the geothermal heat flux is largely unknown for the EAIS, this parameter has to be carefully determined beforehand. In addition, detailed bedrock topography and ice flow history has to be reconstructed for candidates of an Oldest-Ice ice coring site. Finally, we argue strongly for rapid access drilling before any full, deep ice coring activity commences to bring datable samples to the surface and to allow an age check of the oldest ice.