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.

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.

Simultaneous stable isotope analysis of methane and nitrous oxide on ice core samples

Methane and nitrous oxide are important greenhouse gases which show a strong increase in atmospheric mixing ratios since pre-industrial time as well as large variations during past climate changes. The understanding of their biogeochemical cycles can be improved using stable isotope analysis. However, high-precision isotope measurements on air trapped in ice cores are challenging because of the high susceptibility to contamination and fractionation. Here, we present a dry extraction system for combined CH4 and N2O stable isotope analysis from ice core air, using an ice grating device. The system allows simultaneous analysis of δD(CH4) or δ13C(CH4), together with δ15N(N2O), δ18O(N2O) and δ15N(NO+ fragment) on a single ice core sample, using two isotope mass spectrometry systems. The optimum quantity of ice for analysis is about 600 g with typical "Holocene" mixing ratios for CH4 and N2O. In this case, the reproducibility (1σ ) is 2.1‰ for δD(CH4), 0.18‰ for δ13C(CH4), 0.51‰ for δ15N(N2O), 0.69‰ for δ18O(N2O) and 1.12‰ for δ15N(NO+ fragment). For smaller amounts of ice the standard deviation increases, particularly for N2O isotopologues. For both gases, small-scale intercalibrations using air and/or ice samples have been carried out in collaboration with other institutes that are currently involved in isotope measurements of ice core air. Significant differences are shown between the calibration scales, but those offsets are consistent and can therefore be corrected for.

Methodical study of nitrous oxide eddy covariance measurements using quantum cascade laser spectrometery over a Swiss forest

Nitrous oxide fluxes were measured at the Lägeren CarboEurope IP flux site over the multi-species mixed forest dominated by European beech and Norway spruce. Measurements were carried out during a four-week period in October–November 2005 during leaf senescence. Fluxes were measured with a standard ultrasonic anemometer in combination with a quantum cascade laser absorption spectrometer that measured N2O, CO2, and H2O mixing ratios simultaneously at 5 Hz time resolution. To distinguish insignificant fluxes from significant ones it is proposed to use a new approach based on the significance of the correlation coefficient between vertical wind speed and mixing ratio fluctuations. This procedure eliminated roughly 56% of our half-hourly fluxes. Based on the remaining, quality checked N2O fluxes we quantified the mean efflux at 0.8±0.4 μmol m−2 h−1 (mean ± standard error). Most of the contribution to the N2O flux occurred during a 6.5-h period starting 4.5 h before each precipitation event. No relation with precipitation amount could be found. Visibility data representing fog density and duration at the site indicate that wetting of the canopy may have as strong an effect on N2O effluxes as does below-ground microbial activity. It is speculated that above-ground N2O production from the senescing leaves at high moisture (fog, drizzle, onset of precipitation event) may be responsible for part of the measured flux.

Can the carbon isotopic composition of methane be reconstructed from multi-site firn air measurements?

Methane is a strong greenhouse gas and large uncertainties exist concerning the future evolution of its atmospheric abundance. Analyzing methane atmospheric mixing and stable isotope ratios in air trapped in polar ice sheets helps in reconstructing the evolution of its sources and sinks in the past. This is important to improve predictions of atmospheric CH4 mixing ratios in the future under the influence of a changing climate. The aim of this study is to assess whether past atmospheric δ13C(CH4) variations can be reliably reconstructed from firn air measurements. Isotope reconstructions obtained with a state of the art firn model from different individual sites show unexpectedly large discrepancies and are mutually inconsistent. We show that small changes in the diffusivity profiles at individual sites lead to strong differences in the firn fractionation, which can explain a large part of these discrepancies. Using slightly modified diffusivities for some sites, and neglecting samples for which the firn fractionation signals are strongest, a combined multi-site inversion can be performed, which returns an isotope reconstruction that is consistent with firn data. However, the isotope trends are lower than what has been concluded from Southern Hemisphere (SH) archived air samples and high-accumulation ice core data. We conclude that with the current datasets and understanding of firn air transport, a high precision reconstruction of δ13C of CH4 from firn air samples is not possible, because reconstructed atmospheric trends over the last 50 yr of 0.3–1.5 ‰ are of the same magnitude as inherent uncertainties in the method, which are the firn fractionation correction (up to ~2 ‰ at individual sites), the Kr isobaric interference (up to ~0.8 ‰, system dependent), inter-laboratory calibration offsets (~0.2 ‰) and uncertainties in past CH4 levels (~0.5 ‰).

Continuous CO2/CH4/CO measurements (2012–2014) at Beromünster tall tower station in Switzerland

The understanding of the continental carbon budget is essential to predict future climate change. In order to quantify CO₂ and CH₄ fluxes at the regional scale, a measurement system was installed at the former radio tower in Beromünster as part of the Swiss greenhouse gas monitoring network (CarboCount CH). We have been measuring the mixing ratios of CO₂, CH₄ and CO on this tower with sample inlets at 12.5, 44.6, 71.5, 131.6 and 212.5 m above ground level using a cavity ring down spectroscopy (CRDS) analyzer. The first 2-year (December 2012–December 2014) continuous atmospheric record was analyzed for seasonal and diurnal variations and interspecies correlations. In addition, storage fluxes were calculated from the hourly profiles along the tower. The atmospheric growth rates from 2013 to 2014 determined from this 2-year data set were 1.78 ppm yr⁻¹, 9.66 ppb yr⁻¹ and and -1.27 ppb yr⁻¹ for CO₂, CH₄ and CO, respectively. After detrending, clear seasonal cycles were detected for CO₂ and CO, whereas CH₄ showed a stable baseline suggesting a net balance between sources and sinks over the course of the year. CO and CO₂ were strongly correlated (r² > 0.75) in winter (DJF), but almost uncorrelated in summer. In winter, anthropogenic emissions dominate the biospheric CO₂ fluxes and the variations in mixing ratios are large due to reduced vertical mixing. The diurnal variations of all species showed distinct cycles in spring and summer, with the lowest sampling level showing the most pronounced diurnal amplitudes. The storage flux estimates exhibited reasonable diurnal shapes for CO₂, but underestimated the strength of the surface sinks during daytime. This seems plausible, keeping in mind that we were only able to calculate the storage fluxes along the profile of the tower but not the flux into or out of this profile, since no Eddy covariance flux measurements were taken at the top of the tower.

Population mixing and the risk of childhood leukaemia in Switzerland: a census-based cohort study.

Childhood leukaemia (CL) may have an infectious cause and population mixing may therefore increase the risk of CL. We aimed to determine whether CL was associated with population mixing in Switzerland. We followed children aged <16 years in the Swiss National Cohort 1990-2008 and linked CL cases from the Swiss Childhood Cancer Registry to the cohort. We calculated adjusted hazard ratios (HRs) for all CL, CL at age <5 years and acute lymphoblastic leukaemia (ALL) for three measures of population mixing (population growth, in-migration and diversity of origin), stratified by degree of urbanisation. Measures of population mixing were calculated for all municipalities for the 5-year period preceding the 1990 and 2000 censuses. Analyses were based on 2,128,012 children of whom 536 developed CL. HRs comparing highest with lowest quintile of population growth were 1.11 [95 % confidence interval (CI) 0.65-1.89] in rural and 0.59 (95 % CI 0.43-0.81) in urban municipalities (interaction: p = 0.271). Results were similar for ALL and for CL at age <5 years. For level of in-migration there was evidence of a negative association with ALL. HRs comparing highest with lowest quintile were 0.60 (95 % CI 0.41-0.87) in urban and 0.61 (95 % CI 0.30-1.21) in rural settings. There was little evidence of an association with diversity of origin. This nationwide cohort study of the association between CL and population growth, in-migration and diversity of origin provides little support for the population mixing hypothesis.