Estimation and calibration of the water isotope differential diffusion length in ice core records

Palaeoclimatic information can be retrieved from the diffusion of the stable water isotope signal during firnification of snow. The diffusion length, a measure for the amount of diffusion a layer has experienced, depends on the firn temperature and the accumulation rate. We show that the estimation of the diffusion length using power spectral densities (PSDs) of the record of a single isotope species can be biased by uncertainties in spectral properties of the isotope signal prior to diffusion. By using a second water isotope and calculating the difference in diffusion lengths between the two isotopes, this problem is circumvented. We study the PSD method applied to two isotopes in detail and additionally present a new forward diffusion method for retrieving the differential diffusion length based on the Pearson correlation between the two isotope signals. The two methods are discussed and extensively tested on synthetic data which are generated in a Monte Carlo manner. We show that calibration of the PSD method with this synthetic data is necessary to be able to objectively determine the differential diffusion length. The correlation-based method proves to be a good alternative for the PSD method as it yields precision equal to or somewhat higher than the PSD method. The use of synthetic data also allows us to estimate the accuracy and precision of the two methods and to choose the best sampling strategy to obtain past temperatures with the required precision. In addition to application to synthetic data the two methods are tested on stable-isotope records from the EPICA (European Project for Ice Coring in Antarctica) ice core drilled in Dronning Maud Land, Antarctica, showing that reliable firn temperatures can be reconstructed with a typical uncertainty of 1.5 and 2 °C for the Holocene period and 2 and 2.5 °C for the last glacial period for the correlation and PSD method, respectively.

The role of seasonality of mineral dust concentration and size on glacial/interglacial dust changes in the EPICA Dronning Maud Land ice core

We present a record of particulate dust concentration and size distribution in subannual resolution measured on the European Project for Ice Coring in Antarctica (EPICA) Dronning Maud Land (EDML) ice core drilled in the Atlantic sector of the East Antarctic plateau. The record reaches from present day back to the penultimate glacial until 145,000 years B.P. with subannual resolution from 60,000 years B.P. to the present. Mean dust concentrations are a factor of 46 higher during the glacial (~850–4600 ng/mL) compared to the Holocene (~16–112 ng/mL) with slightly smaller dust particles during the glacial comparedto the Holocene and with an absolute minimum in the dust size at 16,000 years B.P. The changes in dust concentration are mainly attributed to changes in source conditions in southern South America. An increase in the modal value of the dust size suggests that at 16,000 years B.P. a major change in atmospheric circulation apparently allowed more direct transport of dust particles to the EDML drill site. We find a clear in-phase relation of the seasonal variation in dust mass concentration and dust size during the glacial (r(conc,size) = 0.8) but no clear phase relationship during the Holocene (0 ice sheet using the size as an indicator for transport intensity, we find that the effect of the changes in the seasonality of the source emission strength and the transport intensity on the dust decrease over Transition 1 can significantly contribute to the large decrease of dust concentration from the glacial to the Holocene.

Chemical compositions of solid particles present in the Greenland NEEM ice core over the last 110,000 years

This study reports the chemical composition of particles present along Greenland’s North Greenland Eemian Ice Drilling (NEEM) ice core, back to 110,000 years before present. Insoluble and soluble particles larger than 0.45 μm were extracted from the ice core by ice sublimation, and their chemical composition was analyzed using scanning electron microscope and energy dispersive X-ray spectroscopy and micro-Raman spectroscopy. We show that the dominant insoluble components are silicates, whereas NaCl, Na₂SO₄, CaSO ₄, and CaCO₃ represent major soluble salts. For the first time, particles of CaMg(CO₃)₂ and Ca(NO₃)₂ 4H₂O are identified in a Greenland ice core. The chemical speciation of salts varies with past climatic conditions. Whereas the fraction of Na salts (NaCl + Na₂SO₄) exceeds that of Ca salts (CaSO₄+ CaCO₃) during the Holocene (0.6–11.7 kyr B.P.), the two fractions are similar during the Bølling-Allerød period (12.9–14.6 kyr B.P.). During cold climate such as over the Younger Dryas (12.0–12.6 kyr B.P.) and the Last Glacial Maximum (15.0–26.9 kyr B.P.), the fraction of Ca salts exceeds that of Na salts, showing that the most abundant ion generally controls the salt budget in each period. High-resolution analyses reveal changing particle compositions: those in Holocene ice show seasonal changes, and those in LGM ice show a difference between cloudy bands and clear layers, which again can be largely explained by the availability of ionic components in the atmospheric aerosol body of air masses reaching Greenland.

Oxygen isotopes of lake marl at Gerzensee and Leysin (Switzerland), covering the Younger Dryas and two minor oscillations, and their correlation to the GRIP ice core

The ratio of oxygen isotopes is a temperature proxy both in precipitation and in the calcite of lacustrine sediments. The very similar oxygen-isotope records from Greenland ice cores and European lake sediments during the Last Glacial Termination suggest that the drastic climatic changes occurred quasi-simultaneously on an extra-regional, probably hemispheric scale. In order to study temporal relations of the different parameters recorded in lake sediments, for example biotic response times to rapid climatic changes, a precise chronology is required. In unlaminated lake sediments there is not yet available a method to provide a high-resolution chronology, especially for periods with radiocarbon plateaux. Alternatively, an indirect time scale can be constructed by linking the lake stratigraphy with other well-dated climate records. New oxygen-isotope records from Gerzensee and Leysin, with an estimated sampling resolution of between 15 and 40 years, match the Greenlandic isotope record in many details. Under the assumption that the main variations in temperature and thus in oxygen isotopes occurred about simultaneously in Greenland and Switzerland, we have assigned a time scale to the lake sediments of Gerzensee and Leysin by wiggle-matching their stable-isotope records with those of Greenland ice cores, which are among the best dated climatic archives. We estimate a precision of 20 to 100 years during the Last Glacial Termination.

Geochemical and Microbiological Studies of Nitrous Oxide Variations within the New NEEM Greenland Ice Core During the Last Glacial Period

Deep polar ice cores provide atmospheric records of nitrous oxide (N₂O) and other trace gases reflecting climate history along with a parallel archive of microbial cells transported with mineral dust, marine and volcanic aerosols from around the globe. Our interdisciplinary study of 32 samples from different depths of the recently drilled NEEM Greenland ice core addressed the question whether the identified microorganisms were capable of post-depositional biological production of N₂O in situ. We used high-resolution geochemical and microbiological approaches to examine the N₂O concentrations, the quantitative distributions of dust, Ca⁺², NH₄⁺ and NO₃⁻ ¡ons related to N cycle pathways, the microbial abundance and diversity at specific NEEM core depths from 1758 m to 1867.8 m. Results showed varying concentrations of N₂O (220 –271.5 ppb). Microbial abundance fluctuated between 3.3 x 10⁴ and 3.3 x 10⁶ cells mL⁻¹ in direct correlation with dust and Ca²⁺ concentrations with higher cell numbers deposited during colder periods. The average values of NH₄⁺ and NO₃⁻ indicated that substrates were available for the microorganisms capable of utilizing them. PCR amplification of selected functional genes involved in bacterial and archaeal nitrification and denitrification was not successful. Sanger and Illumina MiSeq sequence analyses of SSU rRNA genes showed variable representation of Alpha-, Beta- and Gammaproteobacteria, Firmicutes, Actinobacteria, chloroplasts and fungi. The metabolic potential of the dominant genera of Proteobacteria and Firmicutes as possible N₂O producers suggested that denitrification activity may have led to in-situ production and accumulation of N₂O.

Millennial and sub-millennial scale climatic variations recorded in polar ice cores over the last glacial period

Since its discovery in Greenland ice cores, the millennial scale climatic variability of the last glacial period has been increasingly documented at all latitudes with studies focusing mainly on Marine Isotopic Stage 3 (MIS 3; 28–60 thousand of years before present, hereafter ka) and characterized by short Dansgaard-Oeschger (DO) events. Recent and new results obtained on the EPICA and NorthGRIP ice cores now precisely describe the rapid variations of Antarctic and Greenland temperature during MIS 5 (73.5–123 ka), a time period corresponding to relatively high sea level. The results display a succession of abrupt events associated with long Greenland InterStadial phases (GIS) enabling us to highlight a sub-millennial scale climatic variability depicted by (i) short-lived and abrupt warming events preceding some GIS (precursor-type events) and (ii) abrupt warming events at the end of some GIS (rebound-type events). The occurrence of these sub-millennial scale events is suggested to be driven by the insolation at high northern latitudes together with the internal forcing of ice sheets. Thanks to a recent NorthGRIP-EPICA Dronning Maud Land (EDML) common timescale over MIS 5, the bipolar sequence of climatic events can be established at millennial to sub-millennial timescale. This shows that for extraordinary long stadial durations the accompanying Antarctic warming amplitude cannot be described by a simple linear relationship between the two as expected from the bipolar seesaw concept. We also show that when ice sheets are extensive, Antarctica does not necessarily warm during the whole GS as the thermal bipolar seesaw model would predict, questioning the Greenland ice core temperature records as a proxy for AMOC changes throughout the glacial period.

Interhemispheric coupling, the West Antarctic Ice Sheet and warm Antarctic interglacials

Ice core evidence indicates that even though atmospheric CO2 concentrations did not exceed ~300 ppm at any point during the last 800 000 years, East Antarctica was at least ~3–4 °C warmer than preindustrial (CO2~280 ppm) in each of the last four interglacials. During the previous three interglacials, this anomalous warming was short lived (~3000 years) and apparently occurred before the completion of Northern Hemisphere deglaciation. Hereafter, we refer to these periods as "Warmer than Present Transients" (WPTs). We present a series of experiments to investigate the impact of deglacial meltwater on the Atlantic Meridional Overturning Circulation (AMOC) and Antarctic temperature. It is well known that a slowed AMOC would increase southern sea surface temperature (SST) through the bipolar seesaw and observational data suggests that the AMOC remained weak throughout the terminations preceding WPTs, strengthening rapidly at a time which coincides closely with peak Antarctic temperature. We present two 800 kyr transient simulations using the Intermediate Complexity model GENIE-1 which demonstrate that meltwater forcing generates transient southern warming that is consistent with the timing of WPTs, but is not sufficient (in this single parameterisation) to reproduce the magnitude of observed warmth. In order to investigate model and boundary condition uncertainty, we present three ensembles of transient GENIE-1 simulations across Termination II (135 000 to 124 000 BP) and three snapshot HadCM3 simulations at 130 000 BP. Only with consideration of the possible feedback of West Antarctic Ice Sheet (WAIS) retreat does it become possible to simulate the magnitude of observed warming.

Contamination of Alpine snow and ice at Colle Gnifetti, Swiss/Italian Alps, from nuclear weapons tests

Plutonium is present in the environment as a consequence of atmospheric nuclear tests, nuclear weapons production and industrial releases over the past 50 years. To study temporal trends, a high resolution Pu record was obtained by analyzing 52 discrete samples of an alpine firn/ice core from Colle Gnifetti (Monte Rosa, 4450 m a.s.l.), dating from 1945 to 1990. The 239Pu signal was recorded directly, without decontamination or preconcentration steps, using an Inductively Coupled Plasma - Sector Field Mass Spectrometer (ICP-SFMS) equipped with an high efficiency sample introduction system, thus requiring much less sample preparation than previously reported methods. The 239Pu profile reflects the three main periods of atmospheric nuclear weapons testing: the earliest peak lasted from 1954/55 to 1958 and was caused by the first testing period reaching a maximum in 1958. Despite a temporary halt of testing in 1959/60, the Pu concentration decreased only by half with respect to the 1958 peak due to long atmospheric residence times. In 1961/62 Pu concentrations rapidly increased reaching a maximum in 1963, which was about 40% more intense than the 1958 peak. After the signing of the "Limited Test Ban Treaty" between USA and USSR in 1964, Pu deposition decreased very sharply reaching a minimum in 1967. The third period (1967-1975) is characterized by irregular Pu concentrations with smaller peaks (about 20-30% of the 1964 peak) which might be related to the deposition of Saharan dust contaminated by the French nuclear tests of the 1960s. The data presented are in very good agreement with Pu profiles previously obtained from the Col du Dome ice core (by multi-collector ICP-MS) and Belukha ice core (by Accelerator Mass Spectrometry, AMS). Although a semi-quantitative method was employed here, the results are quantitatively comparable to previously published results.

A sublimation technique for high-precision measurements of d13CO2 and mixing ratios of CO2 and N2O from air trapped in ice cores

In order to provide high precision stable carbon isotope ratios (δ13CO2 or δ13C of CO2) from small bubbly, partially and fully clathrated ice core samples we developed a new method based on sublimation coupled to gas chromatography-isotope ratio mass spectrometry (GC-IRMS). In a first step the trapped air is quantitatively released from ~30 g of ice and CO2 together with N2O are separated from the bulk air components and stored in a miniature glass tube. In an off-line step, the extracted sample is introduced into a helium carrier flow using a minimised tube cracker device. Prior to measurement, N2O and organic sample contaminants are gas chromatographically separated from CO2. Pulses of a CO2/N2O mixture are admitted to the tube cracker and follow the path of the sample through the system. This allows an identical treatment and comparison of sample and standard peaks. The ability of the method to reproduce δ13C from bubble and clathrate ice is verified on different ice cores. We achieve reproducibilities for bubble ice between 0.05 ‰ and 0.07 ‰ and for clathrate ice between 0.05 ‰ and 0.09 ‰ (dependent on the ice core used). A comparison of our data with measurements on bubble ice from the same ice core but using a mechanical extraction device shows no significant systematic offset. In addition to δ13C, the CO2 and N2O mixing ratios can be volumetrically derived with a precision of 2 ppmv and 8 ppbv, respectively.

Continuous measurements of methane mixing ratios from ice cores

This work presents a new, field-deployable technique for continuous, high-resolution measurements of methane mixing ratios from ice cores. The technique is based on a continuous flow analysis system, where ice core samples cut along the long axis of an ice core are melted continuously. The past atmospheric air contained in the ice is separated from the melt water stream via a system for continuous gas extraction. The extracted gas is dehumidified and then analyzed by a Wavelength Scanned-Cavity Ring Down Spectrometer for methane mixing ratios. We assess the performance of the new measurement technique in terms of precision (±0.8 ppbv, 1σ), accuracy (±8 ppbv), temporal (ca. 100 s), and spatial resolution (ca. 5 cm). Using a firn air transport model, we compare the resolution of the measurement technique to the resolution of the atmospheric methane signal as preserved in ice cores in Greenland. We conclude that our measurement technique can resolve all climatically relevant variations as preserved in the ice down to an ice depth of at least 1980 m (66 000 yr before present) in the North Greenland Eemian Ice Drilling ice core. Furthermore, we describe the modifications, which are necessary to make a commercially available spectrometer suitable for continuous methane mixing ratio measurements from ice cores.

First continuous phosphate record from Greenland ice cores

A continuous and highly sensitive absorption method for detection of soluble phosphate in ice cores has been developed using a molybdate reagent and a 2m liquid waveg- uide (LWCC). The method is optimized to meet the low concentrations of phosphate in Greenland ice, it has a detection limit of around 0.1ppb and a depth resolution of approximately 2cm. The new method has been applied to obtain phosphate concen- trations from segments of two Northern Greenland ice cores: from a shallow firn core covering the most recent 120yr and from the recently obtained deep NEEM ice core in which sections from the late glacial period have been analysed. Phosphate con- centrations in 20th century ice are around 0.32ppb with no indication of anthropogenic influence in the most recent ice. In the glacial part of the NEEM ice core concentra- tions in the cold stadial periods are significantly higher, in the range of 6–24ppb, while interstadial ice concentrations are around 2ppb. In the shallow firn core, a strong cor- relation between concentrations of phosphate and insoluble dust suggests a similar deposition pattern for phosphate and dust. In the glacial ice, phosphate and dust also correlate quite strongly, however it is most likely that this correlation originates from the phosphate binding to dust during transport, with only a fraction coming directly from dust. Additionally a constant ratio between phosphate and potassium concentrations shows evidence of a possible biogenic land source.

Microgram level radiocarbon (14C) determination on carbonaceous particles in ice

In climate research the interest on carbonaceous particles has increased over the last years because of their influence on the radiation balance of the earth. Nevertheless, there is a paucity of available data regarding their concentrations and sources in the past. Such data would be important for a better understanding of their effects and for estimating their influence on future climate. Here, a technique is described to extract carbonaceous particles from ice core samples with subsequent separation of the two main constituents into organic carbon (OC) and elemental carbon (EC) for analysis of their concentrations in the past. This is combined with further analysis of OC and EC 14C/12C ratios by accelerator mass spectrometry (AMS), what can be used for source apportionment studies of past emissions. We further present how 14C analysis of the OC fraction could be used in the future to date any ice core extracted from a high-elevation glacier. Described sample preparation steps to final analysis include the combustion of micrograms of water–insoluble carbonaceous particles, primary collected by filtration of melted ice samples, the graphitisation of the obtained CO2 to solid AMS target material and final AMS measurements. Possible fractionation processes were investigated for quality assurance. Procedural blanks were reproducible and resulted in carbon masses of 1.3 ± 0.6 μg OC and 0.3 ± 0.1 μg EC per filter. The determined fraction of modern carbon (fM) for the OC blank was 0.61 ± 0.13. The analysis of processed IAEA-C6 and IAEA-C7 reference material resulted in fM = 1.521 ± 0.011 and δ13C = −10.85 ± 0.19‰, and fM = 0.505 ± 0.011 and δ13C = −14.21 ± 0.19‰, respectively, in agreement with consensus values. Initial carbon contents were thereby recovered with an average yield of 93%.