Net accumulation rates derived from ice core stable isotope records of Pío XI glacier, Southern Patagonia Icefield

Pío XI, the largest glacier of the Southern Patagonia Icefield, reached its neoglacial maximum extent in 1994 and is one of the few glaciers in that area which is not retreating. In view of the recent warming it is important to understand glacier responses to climate changes. Due to its remoteness and the harsh conditions in Patagonia, no systematic mass balance studies have been performed. In this study we derived net accumulation rates for the period 2000–2006 from a 50 m (33.2 4 m weq) ice core collected in the accumulation area of Pío XI (2600 m a.s.l., 49°16'40"S, 73°21'14"W). Borehole temperatures indicate near temperate ice, but the average melt percent is only 16 ± 14%. Records of stable isotopes are well preserved and were used for identification of annual layers. Net accumulation rates range from 3.4–7.1 water equivalent (m weq) with an average of 5.8 m weq, comparable to precipitation amounts at the Chilean coast, but not as high as expected for the Icefield. Ice core stable isotope data correlate well with upper air temperatures and may be used as temperature proxy.

Glacial–interglacial dynamics of Antarctic firn columns: comparison between simulations and ice core air-δ15N measurements

Correct estimation of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice core studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: outputs of a firn densification model, and measurements of δ15N of N2 in air trapped in ice core, assuming that δ15N is only affected by gravitational fractionation in the firn column. Since the firn densification process is largely governed by surface temperature and accumulation rate, we have investigated four ice cores drilled in coastal (Berkner Island, BI, and James Ross Island, JRI) and semi-coastal (TALDICE and EPICA Dronning Maud Land, EDML) Antarctic regions. Combined with available ice core air-δ15N measurements from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rates and temperature conditions. Our δ15N profiles reveal a heterogeneous response of the firn structure to glacial–interglacial climatic changes. While firn densification simulations correctly predict TALDICE δ15N variations, they systematically fail to capture the large millennial-scale δ15N variations measured at BI and the δ15N glacial levels measured at JRI and EDML – a mismatch previously reported for central East Antarctic ice cores. New constraints of the EDML gas–ice depth offset during the Laschamp event (~41 ka) and the last deglaciation do not favour the hypothesis of a large convective zone within the firn as the explanation of the glacial firn model–δ15N data mismatch for this site. While we could not conduct an in-depth study of the influence of impurities in snow for firnification from the existing datasets, our detailed comparison between the δ15N profiles and firn model simulations under different temperature and accumulation rate scenarios suggests that the role of accumulation rate may have been underestimated in the current description of firnification models.

Diffusive equilibration of N2, O2 and CO2 mixing ratios in a 1.5-million-years-old ice core

In the framework of the International Partnerships in Ice Core Sciences, one of the most important targets is to retrieve an Antarctic ice core that extends over the last 1.5 million years (i.e. an ice core that enters the climate era when glacial–interglacial cycles followed the obliquity cycles of the earth). In such an ice core the annual layers of the oldest ice would be thinned by a factor of about 100 and the climatic information of a 10 000 yr interval would be contained in less than 1 m of ice. The gas record in such an Antarctic ice core can potentially reveal the role of greenhouse gas forcing on these 40 000 yr cycles. However, besides the extreme thinning of the annual layers, also the long residence time of the trapped air in the ice and the relatively high ice temperatures near the bedrock favour diffusive exchanges. To investigate the changes in the O2 / N2 ratio, as well as the trapped CO2 concentrations, we modelled the diffusive exchange of the trapped gases O2, N2 and CO2 along the vertical axis. However, the boundary conditions of a potential drilling site are not yet well constrained and the uncertainties in the permeation coefficients of the air constituents in the ice are large. In our simulations, we have set the drill site ice thickness at 2700 m and the bedrock ice temperature at 5–10 K below the ice pressure melting point. Using these conditions and including all further uncertainties associated with the drill site and the permeation coefficients, the results suggest that in the oldest ice the precessional variations in the O2 / N2 ratio will be damped by 50–100%, whereas CO2 concentration changes associated with glacial–interglacial variations will likely be conserved (simulated damping 5%). If the precessional O2 / N2 signal will have disappeared completely in this future ice core, orbital tuning of the ice-core age scale will be limited.

Where to find 1.5 million yr old ice for the IPICS "Oldest-Ice" ice core

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 2500 m would be required to find 1.5 Myr old ice (i.e., more than 700 m 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.

High-resolution glacial and deglacial record of atmospheric methane by continuous-flow and laser spectrometer analysis along the NEEM ice core

The Greenland NEEM (North Greenland Eemian Ice Drilling) operation in 2010 provided the first opportunity to combine trace-gas measurements by laser spectroscopic instruments and continuous-flow analysis along a freshly drilled ice core in a field-based setting. We present the resulting atmospheric methane (CH4) record covering the time period from 107.7 to 9.5 ka b2k (thousand years before 2000 AD). Companion discrete CH4 measurements are required to transfer the laser spectroscopic data from a relative to an absolute scale. However, even on a relative scale, the high-resolution CH4 data set significantly improves our knowledge of past atmospheric methane concentration changes. New significant sub-millennial-scale features appear during interstadials and stadials, generally associated with similar changes in water isotopic ratios of the ice, a proxy for local temperature. In addition to the midpoint of Dansgaard–Oeschger (D/O) CH4 transitions usually used for cross-dating, sharp definition of the start and end of these events brings precise depth markers (with ±20 cm uncertainty) for further cross-dating with other palaeo- or ice core records, e.g. speleothems. The method also provides an estimate of CH4 rates of change. The onsets of D/O events in the methane signal show a more rapid rate of change than their endings. The rate of CH4 increase associated with the onsets of D/O events progressively declines from 1.7 to 0.6 ppbv yr−1 in the course of marine isotope stage 3. The largest observed rate of increase takes place at the onset of D/O event #21 and reaches 2.5 ppbv yr−1.

High-resolution mineral dust and sea ice proxy records from the Talos Dome ice core

In this study we report on new non-sea salt calcium (nssCa2+, mineral dust proxy) and sea salt sodium (ssNa+, sea ice proxy) records along the East Antarctic Talos Dome deep ice core in centennial resolution reaching back 150 thousand years (ka) before present. During glacial conditions nssCa2+ fluxes in Talos Dome are strongly related to temperature as has been observed before in other deep Antarctic ice core records, and has been associated with synchronous changes in the main source region (southern South America) during climate variations in the last glacial. However, during warmer climate conditions Talos Dome mineral dust input is clearly elevated compared to other records mainly due to the contribution of additional local dust sources in the Ross Sea area. Based on a simple transport model, we compare nssCa2+ fluxes of different East Antarctic ice cores. From this multi-site comparison we conclude that changes in transport efficiency or atmospheric lifetime of dust particles do have a minor effect compared to source strength changes on the large-scale concentration changes observed in Antarctic ice cores during climate variations of the past 150 ka. Our transport model applied on ice core data is further validated by climate model data. The availability of multiple East Antarctic nssCa2+ records also allows for a revision of a former estimate on the atmospheric CO2 sensitivity to reduced dust induced iron fertilisation in the Southern Ocean during the transition from the Last Glacial Maximum to the Holocene (T1). While a former estimate based on the EPICA Dome C (EDC) record only suggested 20 ppm, we find that reduced dust induced iron fertilisation in the Southern Ocean may be responsible for up to 40 ppm of the total atmospheric CO2 increase during T1. During the last interglacial, ssNa+ levels of EDC and EPICA Dronning Maud Land (EDML) are only half of the Holocene levels, in line with higher temperatures during that period, indicating much reduced sea ice extent in the Atlantic as well as the Indian Ocean sector of the Southern Ocean. In contrast, Holocene ssNa+ flux in Talos Dome is about the same as during the last interglacial, indicating that there was similar ice cover present in the Ross Sea area during MIS 5.5 as during the Holocene.

Dating the Siple Dome (Antarctica) Ice Core By Manual and Computer Interpretation of Annual Layering

The Holocene portion of the Siple Dome (Antarctica) ice core was dated by interpreting the electrical, visual and chemical properties of the core. The data were interpreted manually and with a computer algorithm. The algorithm interpretation was adjusted to be consistent with atmospheric methane stratigraphic ties to the GISP2 (Greenland Ice Sheet Project 2) ice core, (BE)-B-10 stratigraphic ties to the dendrochronology C-14 record and the dated volcanic stratigraphy. The algorithm interpretation is more consistent and better quantified than the tedious and subjective manual interpretation.

CO(2) Diffusion in Polar Ice: Observations from Naturally Formed CO(2) Spikes in the Siple Dome (Antarctica) Ice Core

One common assumption in interpreting ice-core CO(2) records is that diffusion in the ice does not affect the concentration profile. However, this assumption remains untested because the extremely small CO(2) diffusion coefficient in ice has not been accurately determined in the laboratory. In this study we take advantage of high levels of CO(2) associated with refrozen layers in an ice core from Siple Dome, Antarctica, to study CO(2) diffusion rates. We use noble gases (Xe/Ar and Kr/Ar), electrical conductivity and Ca(2+) ion concentrations to show that substantial CO(2) diffusion may occur in ice on timescales of thousands of years. We estimate the permeation coefficient for CO(2) in ice is similar to 4 x 10(-21) mol m(-1) s(-1) Pa(-1) at -23 degrees C in the top 287 m (corresponding to 2.74 kyr). Smoothing of the CO(2) record by diffusion at this depth/age is one or two orders of magnitude smaller than the smoothing in the firn. However, simulations for depths of similar to 930-950m (similar to 60-70 kyr) indicate that smoothing of the CO(2) record by diffusion in deep ice is comparable to smoothing in the firn. Other types of diffusion (e.g. via liquid in ice grain boundaries or veins) may also be important but their influence has not been quantified.

Seasonal Deuterium Excess in a Tien Shan Ice Core: Influence of Moisture Transport and Recycling in Central Asia

Stable water isotope (delta(18)O, deltaD) data from a high elevation (5100 masl) ice core recovered from the Tien Shan Mountains, Kyrgyzstan, display a seasonal cycle in deuterium excess (d = deltaD - 8* delta(18)O) related to changes in the regional hydrologic cycle during 1994 - 2000. While there is a strong correlation (r(2) = 0.98) between delta(18)O and dD in the ice core samples, the regression slope (6.9) and mean d value (23.0) are significantly different than the global meteoric water line values. The resulting time-series ice core d profile contains distinct winter maxima and summer minima, with a yearly d amplitude of similar to 15 - 20parts per thousand. Local-scale processes that may affect d values preserved in the ice core are not consistent with the observed seasonal variability. Data from Central Asian monitoring sites in the Global Network of Isotopes in Precipitation (GNIP) have similar seasonal d changes. We suggest that regional-scale hydrological conditions, including seasonal changes in moisture source, transport, and recycling in the Caspian/ Aral Sea region, are responsible for the observed spatial and temporal d variability.

A 12,000 Year Record of Explosive Volcanism in the Siple Dome Ice Core, West Antarctica

Air mass trajectories in the Southern Hemisphere provide a mechanism for transport to and deposition of volcanic products on the Antarctic ice sheet from local volcanoes and from tropical and subtropical volcanic centers. This study extends the detailed record of Antarctic, South American, and equatorial volcanism over the last 12,000 years using continuous glaciochemical series developed from the Siple Dome A (SDMA) ice core, West Antarctica. The largest volcanic sulfate spike ( 280 mu g/L) occurs at 5881 B. C. E. Other large signals with unknown sources are observed around 325 B. C. E. ( 270 mu g/L) and 2818 B. C. E. ( 191 mu g/L). Ages of several large equatorial or Southern Hemisphere volcanic eruptions are synchronous with many sulfate peaks detected in the SDMA volcanic ice chemistry record. The microprobe "fingerprinting'' of glass shards in the SDMA core points to the following Antarctic volcanic centers as sources of tephra found in the SDMA core: Balenny Island, Pleiades, Mount Berlin, Mount Takahe, and Mount Melbourne as well as Mount Hudson and possibly Mount Burney volcanoes of South America. Identified volcanic sources provide an insight into the poorly resolved transport history of volcanic products from source volcanoes to the West Antarctic ice sheet.

Recent Increase in Black Carbon Concentrations from a Mt. Everest Ice Core Spanning 1860-2000 AD

A Mt. Everest ice core spanning 1860-2000 AD and analyzed at high resolution for black carbon (BC) using a Single Particle Soot Photometer (SP2) demonstrates strong seasonality, with peak concentrations during the winter-spring, and low concentrations during the summer monsoon season. BC concentrations from 1975-2000 relative to 1860-1975 have increased approximately threefold, indicating that BC from anthropogenic sources is being transported to high elevation regions of the Himalaya. The timing of the increase in BC is consistent with BC emission inventory data from South Asia and the Middle East, however since 1990 the ice core BC record does not indicate continually increasing BC concentrations. The Everest BC and dust records provide information about absorbing impurities that can contribute to glacier melt by reducing the albedo of snow and ice. There is no increasing trend in dust concentrations since 1860, and estimated surface radiative forcing due to BC in snow exceeds that of dust in snow. This suggests that a reduction in BC emissions may be an effective means to reduce the effect of absorbing impurities on snow albedo and melt, which affects Himalayan glaciers and the availability of water resources in major Asian rivers. Citation: Kaspari, S. D., M. Schwikowski, M. Gysel, M. G. Flanner, S. Kang, S. Hou, and P. A. Mayewski (2011), Recent increase in black carbon concentrations from a Mt. Everest ice core spanning 1860-2000 AD, Geophys. Res. Lett., 38, L04703, doi: 10.1029/2010GL046096.

Recent Increases in Atmospheric Concentrations of Bi, U, Cs, S and Ca from a 350-Year Mount Everest Ice Core Record

High-resolution major and trace elements (Sr, Cs, Ba, La, Ce, Pr, Nd, Sm, Eu, Tb, Dy, Ho, Er, Tm, Yb, Lu, Bi, U, Tl, Al, S, Ca, Ti, V, Cr, Mn, Fe, and Co) quantified in a Mount Everest ice core ( 6518 m above sea level) spanning the period 1650-2002 AD provides the first Asian record of trace element concentrations from the pre-industrial era, and the first continuous high-resolution Asian record from which natural baseline concentrations and subsequent changes due to anthropogenic activities can be examined. Modern concentrations of most elements remain within the pre-industrial range; however, Bi, U, and Cs concentrations and their enrichment factors (EF) have increased since the similar to 1950s, and S and Ca concentrations and their EFs have increased since the late 1980s. A comparison of the Bi, U, Cs, S, and Ca data with other ice core records and production data indicates that the increase in atmospheric concentrations of trace elements is widespread, but that enrichment varies regionally. Likely sources for the recent enrichment of these elements include mining, metal smelting, oil and coal combustion, and end uses for Bi, and mining and refinement for U and Cs. The source of the synchronous enrichment of Ca and S is less certain, but may be related to land use and environmental change.

Topographic Control of Regional Accumulation Rate Variability at South Pole and Implications for Ice-Core Interpretation

Snow-accumulation rates are known to be sensitive to local changes in ice-sheet surface slope because of the effect of katabatic winds. These topographic effects can be preserved in ice cores that are collected at non-ice-divide locations. The trajectory of an ice-core site at South Pole is reconstructed using measurements of ice-sheet motion to show that snow was probably deposited at places of different surface slope during the past 1000 years. Recent accumulation rates, derived from shallow firn cores, vary along this trajectory according to surface topography, so that on a relatively steep flank mean annual accumulation is similar to 18% smaller than on a nearby topographic depression. These modern accumulation rates are used to reinterpret the cause of accumulation rate variability with time in the long ice-core record as an ice-dynamics effect and not a climate-change signal. The results highlight the importance of conducting ancillary ice-dynamics measurements as part of ice-coring programs so that topographic effects can be deconvolved from potential climate signals.

Ice Core Record of Rising Lead Pollution in the North Pacific Atmosphere

A high-resolution, 8000 year-long ice core record from the Mt. Logan summit plateau (5300 m asl) reveals the initiation of trans-Pacific lead (Pb) pollution by ca. 1730, and a > 10-fold increase in Pb concentration (1981-1998 mean = 68.9 ng/l) above natural background (5.6 ng/l) attributed to rising anthropogenic Pb emissions from Asia. The largest rise in North Pacific Pb pollution from 1970-1998 (end of record) is contemporaneous with a decrease in Eurasian and North American Pb pollution as documented in ice core records from Greenland, Devon Island, and the European Alps. The distinct Pb pollution history in the North Pacific is interpreted to result from the later industrialization and less stringent abatement measures in Asia compared to North America and Eurasia. The Mt. Logan record shows evidence for both a rising Pb emissions signal from Asia and a trans-Pacific transport efficiency signal related to the strength of the Aleutian Low.