Articles | Volume 16, issue 2
Research article
24 Jan 2023
Research article |  | 24 Jan 2023

Laser-induced sublimation extraction for centimeter-resolution multi-species greenhouse gas analysis on ice cores

Lars Mächler, Daniel Baggenstos, Florian Krauss, Jochen Schmitt, Bernhard Bereiter, Remo Walther, Christoph Reinhard, Béla Tuzson, Lukas Emmenegger, and Hubertus Fischer

Precision, accuracy, and temporal resolution are key to making full use of atmospheric trace gas records in ice cores. These aspects will become especially crucial for ice cores that aim to extend the ice core record to the last 1.5 Myr, i.e., across the Mid-Pleistocene Transition (as currently drilled within the European project Beyond EPICA – Oldest Ice Core (BE-OIC)). The ice from this period is expected to be close to bedrock and, due to glacier flow, extremely thinned with 15 000 years of climate history contained in only 1 m of ice. Accordingly, for a century-scale resolution, the sample vertical extent must be reduced to a few centimeters containing only about 1–2 mL air STP.

We present a novel combined system for the extraction and the simultaneous measurement of CO2, CH4, and N2O concentrations, as well as δ13CO2, which achieves a vertical resolution of 1–2 cm (3.5×3.5 cm cross section) with precisions of 0.4 ppm, 3 ppb, 1 ppb, and 0.04 ‰, respectively, in sublimation tests with standard gas over gas-free ice. This is accomplished by employing a directional and continuous laser-induced sublimation followed by analysis of the sample gas by a quantum cascade laser absorption spectrometer (QCLAS). Besides the low sample volume requirements and the vertical resolution capabilities, the described method holds additional advantages over previous methods, including the immunity of the highly specific QCLAS analysis to drilling fluid contamination as well as the non-destructive nature of the spectroscopic gas analysis. The combined extraction and analysis system was extensively tested by sublimating gas-free ice with introduction of a standard gas to determine the accuracy and characterize potential artifacts. Moreover, Antarctic ice samples were measured to confirm the measurement performance, covering the range of variability expected in Pleistocene ice and highlighting the vertical resolution capabilities critical for its application within BE-OIC.

1 Introduction

Ice cores are unique natural archives preserving valuable information on atmospheric greenhouse gas concentrations (CO2, CH4, and N2O) in typically centennial resolution over hundreds of thousands of years. Reconstructing paleo-atmospheric changes in these gases and their isotopic composition has provided valuable insights into their biogeochemical cycles (Blunier and Brook2001; Loulergue et al.2008; Schilt et al.2010a; Mitchell et al.2011; Schmitt et al.2012; Marcott et al.2014; Bereiter et al.2015; Bauska et al.2018; Nehrbass-Ahles et al.2020). However, due to the time-consuming and technically challenging nature of ice core analyses and the limited sample availability, the records are still restricted in precision and resolution. This is particularly pronounced for δ13CO2, where data are only available for selected time intervals during the last glacial cycle and often lack the required precision to make full use of this parameter (Schmitt et al.2012; Schneider et al.2013; Eggleston et al.2016; Bauska et al.2018). These limitations will become even more critical with the advent of a new “Oldest Ice” ice core from East Antarctica (Beyond EPICA – Oldest Ice Core, BE-OIC) that aims to extend the ice core greenhouse gas records over the last 1.5 Myr. The deepest part of such a core will be highly thinned by glacier flow, leading to 15 000 years of atmospheric history being archived in only 1 m of ice (Fischer et al.2013). To make full use of this deep ice core archive, new analytical techniques are required that minimize sample size and increase spatial and temporal resolution and at the same time maximize precision and accuracy of the analyses to be able to reliably quantify past variabilities. A centimeter-scale resolution is a reasonable target for such a technique, as the mean age difference of 150 years between neighboring centimeter samples in such old, thinned ice is similar to the width of the gas age distribution of East Antarctic ice cores. Note that other processes counteracting signal preservation might become significant in the deep ice. For example, the increased vertical concentration gradients in highly thinned ice together with the high in situ temperature of such ice close to bedrock are acting to diminish atmospheric variability by gas diffusion in the ice (Bereiter et al.2014). To quantify and correct such diffusive smoothing will also set new standards on the required precision of the gas analyses.

In the following, we will briefly discuss the prevalent issues in ice core greenhouse gas analyses that still hinder more precise and higher-resolution records. The most common approach to determining CH4 and N2O concentrations in ice core samples is melt extraction (either with or without refreezing) followed by quantification of the concentrations with gas chromatography or laser spectroscopy (e.g., Oyabu et al.2020). Alternatively, CH4 concentrations can be measured in continuous flow analysis (CFA) campaigns with typical depth resolution of 5–10 cm (Stowasser et al.2012; Chappellaz et al.2013; Rhodes et al.2017). The discrete sample method requires a sample size of 20–40 g of ice (Schilt et al.2010a; Mitchell et al.2011; Baumgartner et al.2014; Ryu et al.2018; Schmidely et al.2021), corresponding to a typical depth resolution of several centimeters. CH4 records from melt extraction techniques of discrete ice samples can be affected by production of CH4 during extraction (in extractu) observed in dusty glacial ice from Greenland (Lee et al.2020; Muehl et al.2022). CFA records of CH4 (Stowasser et al.2012; Chappellaz et al.2013; Rhodes et al.2017) appear not to be subject to this issue but suffer from accuracy problems due to the solubility of CH4 in the meltwater stream.

For CO2, any contact with liquid water during the extraction procedure must be strictly avoided to circumvent chemical reactions (for example acid–carbonate reactions or oxidation of organic carbon compounds) because they affect the measured CO2 content. For this reason, the analysis of CO2 and CH4/N2O was mostly separated in the past despite the larger amount of sample required for both analyses. CO2 measurements are typically performed using mechanical extraction devices (crushing, milling, grating, etc.) combined with either optical gas analyzers (e.g., near-infrared (NIR) or laser spectroscopy) or gas chromatography (GC) (Lüthi et al.2010; Bereiter et al.2013; Marcott et al.2014; Nehrbass-Ahles et al.2020; Shin et al.2020). The Bern cracker, Grenoble mill, or the Bern Centrifugal Ice Microtome (CIM) achieve a vertical sample size of 2–5 cm (corresponding to as low as 15 g of ice). Typical CO2 precisions achieved using these methods are 1–3 ppm. The sample throughput ranges from a few measurements up to 20 measurements per day in the case of the Bern CIM.

These dry-extraction techniques release less than 100 % of the gases trapped in the ice (Bereiter et al.2015) and have varying extraction efficiencies for air bubbles and different clathrate sizes, which leads to systematic effects in CO2, mainly in and below the brittle ice zone (Lüthi et al.2010). The only dry-extraction technique for CO2 from ice core samples that enables 100 % extraction efficiency is sublimation (Schmitt et al.2011). While current sublimation extraction methods with a throughput of one to three samples per day are slow in comparison to the CIM, the unfractionated recovery of air also enables reliable and precise δ13CO2 analyses in partially or fully clathrated ice (Schmitt et al.2011; Schneider et al.2013). The interconnected limitation of sample size and precision is especially pronounced for δ13CO2. Previous work investigating CO2 ice core samples for their carbon isotopic composition using gas chromatography–mass spectrometry (GC–MS) or dual-inlet isotope ratio mass spectrometry (IRMS) achieved reproducibilities of 0.06 ‰–0.08 ‰ (Schmitt et al.2011; Schneider et al.2013) on samples as small as 30 g of ice. Bauska et al. (2014) obtained 0.02 ‰, however requiring samples that are 10-fold larger. All of these methods required separation of CO2 (either via GC or cryogenic trapping) from other gases for the analysis. Finally, interference by drill fluid traces in extracted air samples may affect the gas analyses. Such drill fluid contamination becomes more frequent in the brittle and the very deep ice, where mechanical fracturing enables the contaminant to penetrate the ice core.

To achieve the necessary resolution for the oldest ice while avoiding all issues outlined above, we developed a novel extraction and analysis method using continuous sublimation extraction of ice core sections coupled to a custom-made dual-quantum cascade laser absorption spectrometer (QCLAS; Bereiter et al.2020). In particular, using sublimation, the new method avoids the water phase that leads to in extractu formation of CO2 and CH4 and ensures a 100 % extraction efficiency for all gases. Furthermore, the QCLAS is insensitive to drill fluid traces and is designed to simultaneously measure CO2, CH4, and N2O concentrations as well as δ13CO2 directly from 1–2 mL of extracted air (Bereiter et al.2020). In the following, we present the technical details of the developed method comprising the laser-induced sublimation extraction (LISE) system, the vacuum system, the standard gas inlet, and the calibration procedure. While the QCLAS is described in detail by Bereiter et al. (2020), its particular use and the exact measurement procedure, including modification to the inlet and the calibration/correction procedure for the ice core air sample measurements, are described here. A detailed evaluation of the LISE using this measurement setup follows, and finally, we analyze samples from well-characterized ice cores (EPICA Dome C (EDC) and Talos Dome) to demonstrate the excellent performance of the LISE/QCLAS system.

2 Laser sublimation extraction method

In this section, we describe the laser-induced sublimation extraction (LISE) system, starting with the sublimation apparatus, followed by the vacuum system, the standard gases, and then the LISE measurement procedure itself.

2.1 Sublimation extraction

The core of the LISE consists of a 2 L extraction vessel, in which a rectangular block of ice is placed (Fig. 1). The vertical dimension of the ice block is typically 9 cm, which is sufficient for up to six consecutive subsamples. The cross section is variable and can be adapted to the dimensions of the ice samples, but typically a cross-section of 10–15 cm2 is required to achieve a vertical sample resolution of 1–2 cm. The air in the ice is released through sublimation of the ice block from top to bottom by vertical illumination with a 150 W near-infrared (NIR) light source. The resulting water vapor is frozen to the cooled walls of the extraction vessel, which leads to low vapor pressure. Consecutively, the released air is cryogenically trapped in dip tubes cooled to 15 K.

Figure 1Schematic view of a cross section through the extraction vessel during sublimation. The NIR laser beam induces sublimation at the top ice surface. A small fraction of the beam misses the ice sample and hits the underlying gas-free ice block, which prevents further reflection within the vessel. The resulting water vapor is efficiently removed by freezing to the walls of the vessel, while the released gases are evacuated out of the vessel via a glass tube from the bottom of the vessel.


The laser (BrightLase® Ultra-500 diode laser, QPC Lasers) emits at a wavelength of 1530 nm. At this wavelength, the penetration depth in ice is less than 1 mm (Warren and Brandt2008), which prevents subsurface melting and minimizes sublimation of the sides of the ice block (see Sect. 2.3.1). Further, the wavelength is within the transmission regime of common viewports on a borosilicate/Pyrex basis. The sublimation of a 1 cm high 3.5 × 3.5 cm2 ice cuboid requires ca. 33 kJ, which corresponds to 20 min irradiation with 25 W or 2 W cm−2. The 150 W maximum optical power of the laser therefore gives a large margin for losses of light before absorption by the ice. The laser beam is transmitted via a fiber optic patch cable (MHP550L02, Thorlabs) and initially has a Gaussian spatial power distribution. An optical homogenization unit (Bayerisches Laserzentrum GmbH) is used to widen and homogenize the beam and to reshape it into a quadratic cross section, while a specifically shaped aperture further forms it to match the ice block cross section (Figs. 1 and 2). The unit consists of a collimator, two arrays of cylindrical lenses that split the incoming beam into different beamlets, and a Fourier lens that superposes the beamlets at its focal length of 25 cm into a 4×4 cm2 cross section (e.g., Zimmermann et al.2007). The optical system is mounted on a XY linear translation stage, the position of which can be adjusted in the vertical direction. A visible (635 nm) trace laser facilitates the alignment.

Figure 2(a) Cross section of the laser beam (at 1530 nm) after passing through the homogenization unit projected on a VRC4 IR detector card showing somewhat higher power densities at the edges. (b) Mask used to further shape the beam. The beam size of the visible trace laser is about 5 % larger than that of the NIR laser.


The extraction vessel (Fig. 1) consists of a sealed glass vessel with a flat base and lid (approx. 100 mm diameter) with a Kovar glass adapter to a ConFlat flange (DN100 CF). A CF spacer flange with two DN16 CF ports connects the vessel to the rest of the vacuum line. The top flange, where the laser beam enters the vessel, is a DN100 CF viewport (Pfeiffer), which can be removed for loading ice samples.

In principle, ice sublimates at any pressure below the triple point of water (6 hPa). However, we empirically found that the pressure in the extraction vessel should be kept below 0.2 hPa to avoid the buildup of needle-like sublimation structures (further on referred to as “penitentes”; see Sect. 2.3.1). Cooling the walls of the glass vessel with a temperature-controlled (ca. −120C) aluminum ring removes most of the vapor by freezing it out onto the vessel wall. The ring is cooled by supplying liquid nitrogen (LN2) drops to its inner compartment, which is filled with copper swarfs (size ca 2–3 mm; see also Schmitt et al.2011). The swarfs increase the surface area and tortuosity, thus facilitating evaporation of the LN2 rather than pooling at the bottom of the ring. Taking into account the different thermal contraction coefficients of aluminum and glass, we left a gap of about 1 mm between ring and glass at room temperature. The ring and the glass part of the vessel are placed in a dewar, and the gap between the dewar edge and the vessel is covered by a thermal insulator (polyurethane foam). A PID controller (diaTRON 116, Jumo) regulates the LN2 supply and stabilizes the temperature of the ring within 0.2 C.

All the metal surfaces of the extraction vessel were coated with gold to reduce adsorption and desorption of gases, prevent oxidation, and in particular to suppress any CO2 formation from potential carbon impurities in the steel walls and welds. Despite the gold plating, desorption of gases still occurs when the walls get warmer during the sublimation by scattered NIR light. Accordingly, the top flanges are cooled to −15C by an air stream at −110C of pressurized air cooled via a copper heat exchanger mounted in a 2 L dewar (Schmitt et al.2011) to minimize such effects. In a similar fashion to the cooling ring, LN2 is dripped into the dewar, at a rate regulated by a PID controller (West 2300, UK) to reach and maintain the targeted temperature (measured by a thermocouple) of the exiting air stream.

The entire extraction system, consisting of the homogenizer, extraction vessel, and LN2 supply, is surrounded by an aluminum casing for laser safety. Part of the casing can be removed for loading of the sample, refilling of LN2, or readjusting of the homogenizer position.

2.2 Vacuum system and extraction line

Once the gases are released from the ice, they are pushed (similar to a diffusion pump) downwards by the vapor stream that originates from the ice surface to the cooled outer surfaces of the vessel where the water vapor freezes (Fig. 1). To prevent gas accumulation and fractionation at the bottom of the vessel downstream of this “water vapor pump”, we evacuate the vessel through an L-shaped glass tube from the bottom of the extraction vessel. The upper end of this tube is connected to one of the exit ports of the middle flange, leading to the extraction line that consists of up to V3 of 3/8 in. stainless steel tubing (bold line in Fig. 3). The gases follow the pressure gradient caused by freezing the residual water vapor flux in a second, external water trap at −100C, and finally the dry gas is collected in dip tubes cooled to 15 K in a helium exchange cryostat acting as a cryopump at the end of the line.

Figure 3Schematic diagram of the laser-induced sublimation extraction (LISE) system. Tubing size upstream of V21 and V23 is mostly 1/4 in., except in the sample pathway indicated by bold lines (3/8 in.) and upstream of V12 (1/8 in.). Individual components are described in the text.


The external water trap (Fig. 3) consists of a 10 cm long cooled aluminum body assembled around 3/8 in. tubing. Again, the temperature of the aluminum body is regulated within 0.2 C using a PID controller (diaTRON 116, Jumo) via a LN2 supply to bores filled with steel wool in the aluminum body. The cryostat (CCS-500/204, Janis) used to cool the dip tubes to 15 K is an exchange gas closed-cycle cold trap system with six dip tube ports. Each dip tube consists of a sealed off 1/4 in. SilcoNert® coated tube connected with a VCR connector to a full-metal valve (Fujikin, FUDDFM-71G-6.35). Six dip tubes are attached to a custom made aluminum manifold (hereafter called Hexatubus, photo in Appendix Fig. A1) that connects with a 1/4 in. tube to the rest of the extraction line at V3 for consecutive automated filling. The six dip tubes are consecutively filled with the gas released from the ice sample. After the extraction, the entire Hexatubus is transferred to the QCLAS system (see Sect. 3) for the gas analysis.

Besides removing the released gas from the extraction vessel, the middle flange of the extraction vessel also allows standard gases to be introduced at a rate similar to the gas release rate of the ice core samples, to mimic the extraction process and quantify potential blank contributions or fractionation processes. To this end, the second port of the extraction vessel connects via valve V1 to the outlet of the standard capillary that is used to introduce standard gases (see Sect. 2.4) to the extraction system (Fig. 3).

All valves in contact with sample or standard gas at low pressure, i.e., V20, V22, V1, V2, V3, and V4, are full-metal pneumatic valves (Fujikin FWBR-71 series). All seals in the extraction system are either CF copper or VCR stainless steel gaskets. The high vacuum is maintained by two pump stations, each equipped with a turbomolecular pump (Turbo HiPace 80, Pfeiffer) backed by a diaphragm pump (MVP 015-4, Pfeiffer). A similar diaphragm pump is used to evacuate the standard gas inlet and in a first stage the extraction vessel to a pressure below 1 hPa, which is measured using a Pirani gauge (P1, TPR 280, Pfeiffer). Both pump stations are connected to the extraction line via a LN2 trap that removes any remnant water vapor, hence improving the vacuum and assuring us that the pressure monitored at the turbomolecular pump essentially reflects dry-gas pressure. The pressures P20 and P21 at the pump stations are measured by a Pirani/cold cathode combination gauge (PKR 361, Pfeiffer). In the extraction line, we installed capacitance gauges (P2, CMR 363, and P3, CMR 364, Pfeiffer) due to their suitable materials exposed to the vacuum side, i.e ceramic and stainless steel, and good accuracy with low zero drift.

2.3 Characterization of the LISE

2.3.1 Sublimation front

The net sublimation rate can be described as a modified Hertz–Knudsen equation, and hence it is proportional to the difference between the vapor pressure above the ice surface and the equilibrium vapor pressure at ice temperature, with a proportionality coefficient that strongly decreases towards lower temperature (Kossacki et al.1999). With a heat source (here supplied by the NIR laser), the ice temperature will rise and increase the sublimation rate until a steady state is reached at the surface or it starts to melt. Hence, to sustain sublimation, the vapor needs to be transported away and frozen at the internal trap. For a top surface sublimation speed of 0.05–0.1 cm min−1, the surface of the ice needs to have a temperature of −27 to −20C, thus far below the melting point. The heat input from the laser beam necessary to sustain such sublimation rates is 2–4 W cm−2 and almost independent of ice temperature as sublimation heat only varies within 1 % over a sublimation temperature range from 0 to −100C (Feistel and Wagner2007).

The cooling process of the ice during the pump-down phase before the NIR laser is switched on is mainly caused by the heat loss due to sublimation and decreases the temperature of the ice surface until the loss reaches a steady state with the thermal radiation from the surrounding vessel. The residual flow of air liberated from the ice is about 2 orders of magnitude smaller than during the extraction, depending on both the size of the ice sample (total surface vs. irradiated surface) and of the laser power. Thus, during the pump-down process, we expect about a 5 % loss of the ice volume, which also efficiently removes any adsorbed gases from the ice sample surface.

Irradiating the top surface of the ice increases its temperature and hence the sublimation rate but also induces a temperature gradient across the ice. A temperature gradient of 10 K cm−1 in the ice block removes around 2–3 W from the sublimation front as the thermal conductivity of ice ranges from 2.2 W mK−1 at 0 C to 3.2 W mK−1 at −100C (Fukusako1990). While such a heat loss can be easily compensated for by higher irradiation power, it may lead to increased temperature on the flanks of the ice, inducing unwanted sublimation at these locations, thereby curving the ice surface and thus decreasing the effective resolution of the extraction.

To attain a uniform and horizontal sublimation front, we aimed at a homogeneous square beam shape. Despite the homogenizer beam shaping, however, the power distribution still exhibited higher density at the edges of the square (roughly 50 % higher intensity in the outer 2.5 mm, Fig. 2) and led us to use an aperture as a mask to cut off the outer energy maxima, leading to a more homogeneous illumination of the top surface. The mask with a side length of 1.5 cm was placed 9 cm below the Fourier lens. However, we empirically found that using a square mask, the ice at the corners of the ice block retreated more slowly than the main surface, resulting in 1–2 cm high ice peaks at each corner after 6 cm of sublimation. We adjusted the aperture design to increase irradiation at the corners of the ice by drilling 3 mm holes centered at each corner of the mask, which improved the flatness of the sublimation front considerably.

2.3.2 Stacked ice resolution experiments

The vertical resolution of the gas extraction that can be achieved with this setup was investigated by preparing ice consisting of stacked sections alternating between gas-free ice (FI) and bubble ice (BI). Observing the released gas flow while sublimating over the ice sections reveals the blurring of transitions and whether side sublimation of lower sections contributes to the gas flow. We prepared a 9 cm long ice sample consisting of stacked 1×3.5×3.5 cm3 sections alternating between FI and BI and another 7 cm long ice sample consisting of alternating 2 cm sections of the same lateral dimensions. The bubble ice sections were cut from the Byrd ice core (1797 m depth, tube no. 1147 C) with an air content somewhere between 0.10 and 0.11 mLSTP g−1 (Raynaud and Whillans1982). The sections were “glued” together by carefully distributing a small amount of Milli-Q water on the surface before stacking them. Even though Milli-Q water is not gas-free, its contribution to the measured air flow is minimal.

The ice was then sublimated as described in Sect. 2.5 with a laser power of 100 W down to ca. 1/4 of the bottom-most bubble ice section. The released gases increase the pressure above TP2 detected by P21 (Fig. 3), which is proportional to the gas flow as all water vapor has been removed by the traps. A rough calibration of the P21 pressure–flow relationship was achieved by admitting gas through the standard gas capillary at known flow rates into the system. The transitions from FI to BI are used as depth markers and enable conversion from time to length. Note that the conversion is assumed linear, which leads to an overestimation of distance in the last ice section due to lower laser intensity.

Figure 4Released air flow (solid lines, “measured”) during sublimation (100 W laser power) of alternating stacked gas-free ice (FI) and bubble ice (BI). Blue is alternating 1 cm ice layers, starting with FI. Green is 2 cm layers starting with BI. The dashed lines mark comparable minima during the sublimation of FI that form extraction flow “valleys”. For the sublimated distance, a conversion factor of 1 cm per 1050 s is used.


The extracted air flow vs. sublimated depth of the 1 cm stacked ice (Fig. 4) follows a slightly deformed periodic square wave with resolved peaks, demonstrating that the achieved resolution is better than 1 cm. The minima between peaks, i.e., during sublimation of FI sections, do not decrease with depth, indicating that only gases released due to side sublimation of the following BI section significantly contribute to the gas flow but not of the lower BI sections. The gas flow during the FI section is, therefore, the sum of the decreasing flow from the previous and the increasing flow from the following BI section. While in the 1 cm stack the gas flow from two BI sections does mix in the middle part of a transition through a FI section, in the 2 cm stacked ice, both contributions reach the background level before the start of the sublimation of the next BI, indicating that the two peaks are fully separated (Fig. 4). The steeper slope of flow from BI to FI than from FI to BI is a result of gas emission due to side sublimation. The emission is not present in the former as the (FI) ice sublimated on the side is gas-free. In both transitions, the non-uniform sublimation fronts contribute to the blurring of the transition. The non-uniform sublimation front results roughly in an exponential increase or decrease when sublimated over a transition with a half-value depth of 0.1–0.2 cm. The side sublimation contribution increases almost linearly within 0.5 cm to around 20 % of peak height approximately 0.2–0.3 cm before peak height, where the exponential increase from gases released at the sublimation front starts to become significant. Note that the pump time required to reach 25 % of initial pressure is around 100 s, corresponding to 0.1 cm of sublimated ice; mixing after release does, therefore, not significantly decrease resolution. In conclusion, two single pulses can be distinguished within 0.2–0.4 cm (roughly double the distance in which the signal drops to 50 % of its initial value) and are fully separated at a spatial distance of 1 cm. While we have shown that the contribution of side sublimation from lower layers is small (a few percent at most), it may still lead to some smoothing across the subsamples if the concentration gradient in a sample is large.

2.3.3 Penitentes

A faster sublimation is beneficial considering the potential contamination from gases desorbed from the surfaces. While optical power loss can be compensated for using a higher laser power, a phenomenon appearing on the sublimating surface that might impact resolution is the occurrence of penitente structures. In nature, penitentes are needle-like sublimation ice structures at much larger scale found in high-altitude dry mountain regions, e.g., in the Andes (Bergeron et al.2006). The formation of such jagged ice spikes is linked to a surface instability driven by the combination of sublimation of the ice and reflection of the sunlight by surrounding features. The redistribution of incident radiation into initially small depressions in the ice leads to locally higher sublimation rates and further increases the surface depressions. This positive feedback mechanism can result in nature in penitentes up to several meters in height with tilts toward the predominant direction of solar radiation (Bergeron et al.2006; Cathles et al.2014).

We observed similar penitente-like structures of much smaller size appearing and evolving in time that negatively impact the sublimation front and thus limit the vertical resolution of the extraction. Their appearance was observed only when the pressure measured at gauge P2 was higher than 0.2 hPa (Fig. 5). The size of such penitentes can exceed 1 cm after 3–4 cm of sublimation and “grow” further. The pressure threshold slightly depends on the gas mixture in the vessel, as in gas-free ice the transition to the buildup of penitentes is lower than for bubble ice. We speculate that a higher pressure inhibits efficient water vapor transport away from the ice surface to play a major role in the buildup of penitentes (Berisford et al.2021). Setting the internal trap temperature to −120C and restricting the sublimation speed to below 1.0 cm per 20 min keeps the pressure in the extraction vessel below 0.15 hPa. With this configuration, the formation of penitentes can be avoided.

Figure 5Ice after 4 cm of sublimation. Penitentes in panel (b) occur at pressures of more than 0.15–0.2 hPa measured at the top of the extraction vessel, while in panel (a) the pressure was kept below 0.15 hPa.


2.4 Standard gas inlet and calibration gases

The principle of identical treatment (Werner and Brand2001) of sample and reference material cannot be completely adhered to in our case, due to the lack of a suitable ice reference material for gases. Thus, similar to Schmitt et al. (2011), we approximate the identical treatment by admitting a calibration gas into the extraction vessel while sublimating gas-free ice. The flow rate is adjusted to match the gas flow released from ice core samples. Subsequent steps are kept identical to the ice core sample procedure (for details, see Sect. 2.5). The standard gases from pressurized cylinders are first introduced to the standard venting chamber via V12 up to a pressure set by the respective regulator at the cylinder (depending on the selected flow rate around 0.1 MPa over-pressure). A vent capillary (fused silica, 100 µm i.d., 0.5 m) keeps the system flushed at a rate of 1–2 mLSTP min−1 to avoid isotopic fractionation. The capillary (fused silica, 50 µm i.d., 1 m) through which the gases are admitted into the extraction line induces a laminar flow. Due to the pressure at the outlet that is 4 orders of magnitude smaller, this flow is only dependent on the pressure P10 at the high-pressure end.

Table 1Reference values of custom-made standard gases prepared for the calibration of the QCLAS from Bereiter et al. (2020), no. 1–7; and reference air obtained from INSTAAR, no. B4 and B3. The concentrations and isotope ratio values cover the range of expected variations found in air samples from ice cores. The standard no. 2 is matrix air (free of CO2 and CH4) that was used to generate the various standard mixtures by adding different amounts of trace gases. The uncertainty of the values is given in parentheses as 1σ.

Download Print Version | Download XLSX

For the calibration of the LISE (see Sect. 4.1) and the QCLAS (see Sect. 3.2), and to link the CO2 isotope ratios derived by the QCLAS to the international Vienna PeeDee Belemnite (VPDB) scale, we prepared a set of calibration gases (Table 1) with specific mixtures covering the full range from the preindustrial to glacial atmospheric composition in CO2, CH4, and N2O and a range of 3.54 ‰ to 10 ‰ in δ13CO2. The production procedure and further details are discussed in Bereiter et al. (2020). Note that standard no. 1 and no. 6 have similar isotopic composition within 0.02 ‰ and CO2 concentrations, which are close to the values expected in ice core samples and are therefore used as the preferred working standards for characterization of the LISE and as QCLAS reference gases. Two additional air reference cylinders with preindustrial atmospheric composition in CO2, however, with present-day δ13CO2 signature (Table 1), were used as independent targets to determine the long-term accuracy and precision of the QCLAS. These cylinders were prepared at the Stable Isotope Lab at the Institute for Arctic and Alpine Research (INSTAAR; University of Colorado, USA).

2.5 LISE procedure

To avoid contamination of the ice core sample during processing, the extraction vessel (EV) is rinsed after each use with ethanol followed by ultrapure water and baked out at 195 C overnight in a vacuum oven. The EV is then remounted to the extraction apparatus. Prior to an extraction, the EV is evacuated again overnight combined with heating to 100 C to remove water vapor and potential contamination from the inner surface and passivate the metal surfaces once more. The six dip tubes of the Hexatubus are placed into the (warm) cryo-cooler and attached to the system. The dip tubes and the extraction line are evacuated with a turbomolecular pump overnight, reaching pressures below 10−7 hPa.

The extraction procedure starts by turning on the He cryostat, as it takes more than 2 h to reach the 15 K needed to trap the gas components in the dip tube. The EV is isolated from the extraction line and filled with argon gas (Ar; ALPHAGAZ 1) to avoid ambient air entering the system during the ice sample loading. The internal water vapor trap around the EV is cooled down to −25C before loading the EV with the ice core sample. During the cool-down time, the ice core sample is prepared.

The ice core sample is cut from an ice core section with a band saw in a cold room at −20C. To minimize the risk of contamination, a layer of at least 5 mm is removed from each surface of the bulk ice core sample. Furthermore, the surfaces of the ice core sample in contact with the band saw are further shaved with a stainless steel microtome knife. After this procedure, the weight of a 9 cm tall ice sample amounts to 70–80 g. To load the ice core sample, the viewport is unmounted, while a steady flow of Ar limits the invasion of laboratory air into the EV. First, a pre-cooled glass cube (6×5×5 cm3) is placed into the EV, followed by a layer of gas-free ice (2×5×5 cm3) that prevents the NIR laser beam from hitting the glass cube. The ice core sample is then added and aligned to the NIR laser beam with the help of the visible trace laser. For all handling steps, polyethylene gloves are worn over regular gloves to reduce potential contamination. After re-mounting the viewport and sealing the EV, the Ar flow is switched off, and a diaphragm pump is used to evacuate the extraction volume via V1 and V20. Meanwhile, the internal water vapor trap is set to −70C, and the metal flange parts are cooled to −15C with an external cooling ring. After 30 min of pumping, the (predominantly water vapor) pressure in the extraction vessel reaches 0.1 hPa, and a turbomolecular pump takes over the vessel evacuation. The internal trap temperature is then lowered to −120C, while the external water trap, in the line connecting the EV to the Hexatubus, is set to −100C, hindering any water vapor emanating from the EV from reaching the dip tubes. After another 30 min of turbomolecular pump-down, the pressure at P20 reaches 2×10-6 hPa. During the entire pump-down procedure, 1–2 mm of ice on each surface is sublimated away, which acts as an additional cleaning step (see also Sect. 2.3.1).

After finishing the pump-down, V1 is closed, and the gas stream is routed into the extraction line via V2. The sublimation is initiated by setting the NIR laser to 100 W. Note that due to radiation losses, only about 30 % of the power is available to sublimate ice (see Sect. 2.3.2). To stabilize the gas flow and equilibrate the extraction line surfaces, the liberated sample gas is wasted into the turbomolecular pump for the first 10 min of the sublimation process. A longer pumping time did not improve the results and leads to unnecessary loss of sample. Afterwards, collection of the gas sample starts by consecutive opening and closing of the six dip tubes for 25 min each (corresponding to  1 mL STP each). By integrating the readouts of the pressure gauge P3, we can estimate the total amount of gas collected in one dip tube with an accuracy of 1 %. With six collected air samples from a 9 cm ice core sample, the vertical resolution is 1.5 cm. After completing the filling procedure, the entire Hexatubus is detached from the extraction system and removed from the cryo-cooler. The trapped sample gases inside the dip tubes are left for at least 1 h to homogenize prior to the QCLAS measurement.

The standard gas over gas-free ice procedure only differs from the above in that during the sublimation process, standard gas is introduced through the capillary into the extraction vessel at a rate similar to the one released from an actual ice sample. All pressure readings and valve switches are performed automatically by a custom-made LabVIEW program. The only task that needs to be performed manually during the sublimation is to periodically lower the homogenizer unit of the NIR laser to keep the distance to the top of the sublimating ice sample comparable.

3 QCLAS gas analysis

The QCLAS system, which performs the actual measurements of greenhouse gas concentrations and the isotopic composition of CO2, is described and discussed in detail by Bereiter et al. (2020). In this section we provide a short summary of its main characteristics and the measurement procedure and describe additional modifications that were necessary to ensure the fidelity of the (isotopic) analysis of small air samples derived from ice core sublimation. A series of measurements of reference gases over 6 months were used to assess the spectrometer's long-term stability and reproducibility.

3.1 Upgraded QCLAS system

The key elements of the system are the dual-laser concept for simultaneous measurement of CO2, CH4, and N2O concentrations as well as δ13CO2, the custom-developed, low-volume multipass absorption cell (MPC), custom-made high-performance data acquisition and laser driving electronics and the vacuum manifold that allows for a non-fractionating admission of standard gases or samples into the absorption cell. A multi-port selector (6-port, Vici Valco Instruments Inc.) with individual vents allows for switching between the standard gas cylinders while maintaining an uninterrupted gas flow. A calibration scheme was implemented to derive transfer functions that relate the instrument response to externally calibrated standard gases for a range of pressures in the absorption cell. For the calibration the standard gases no. 3–no. 7 in Table 1 are used, which span the full range of expected values in CO2, CH4, and N2O concentrations as well as δ13CO2 isotopic composition in ice cores.

With this setup, we demonstrated sample-standard gas repeatabilities (1σ) for discrete gas samples of 1 mL STP amount (corresponding to 10–15 g ice) of 0.03 ppm for CO2, 2.2 ppb for CH4, 1 ppb for N2O, and 0.04 ‰ for δ13CO2 (Bereiter et al.2020). Further we showed that the laser spectroscopic approach is immune to potential contaminating compounds found in drilling fluids.

Figure 6Gas handling system used to introduce the small air samples from dip tubes into the QCLAS. The entire system is made of ultra-high-vacuum (UHV) stainless-steel tubing and uses exclusively VCR seals in the low-pressure part downstream of the flow restrictions (critical orifices and fused silica capillary) up to the multi-pass cell. The valves are pneumatically actuated, full-metal valves (Fujikin Inc., Japan).


To carry out the standard gas admission and mixing described below, the gas handling system (see Fig. 6) has been modified from Bereiter et al. (2020) and now includes two mass-flow controllers (Analyt MTC, 358 series) along with a critical orifice (20 µm hole diameter) and a capillary flow restriction, respectively. Fast switching between gases with the multi-port selector has proven inadequate since it takes more than 10 min for the trace gas concentrations to stabilize after a switch. Thus, a direct line circumventing the multi-port selector has been added, allowing for the admission of two different standard gas streams into the MPC in rapid succession. This second line enables the analysis of different standard gases within minutes of each other, as required by the measurement procedure, and also improves the determination of the calibration functions by providing an unchanging reference, whilst the calibration gases are cycled through via the multi-port selector.

3.2 QCLAS measurement procedure

To achieve the highest accuracy for the dip tube analysis, we use a pointwise offset correction approach. This is required to account for spectrometer-related artifacts causing dependencies of the retrieved δ13CO2 values on gas pressure and CO2 amount fraction. We decided to minimize these effects rather than just correct for them, mainly because we observed a variability of the dependencies over time. The dip tube measurement is preceded and followed by standard gases in the following manner (Fig. 7): first an aliquot of standard gas no. 6 is admitted into the absorption cell at 5.5 hPa, which serves as a conditioning step and control measurement. Then the MPC is fully evacuated, and the sample gas is expanded from the dip tube into the MPC. After another full evacuation, the MPC is filled with a standard gas mixture (no. 2 and no. 7) that is concentration-matched within 1 ppm and pressure-matched within 1 Pa to the preceding sample. This standard is used later to pointwise-correct the concentration-specific offset in δ13CO2. The concentration matching is achieved by mixing a constant gas flow (approximately 5 mL min−1) of standard no. 7 with a variable gas flow of the CO2-free standard no. 2 (Fig. 6). The parameters of the mass-flow controller for standard no. 2 are calibrated beforehand and cover a range from 150 to 300 ppm. Finally, a pressure-matched aliquot of standard no. 6 is admitted into the cell. This standard is used later to pointwise-correct any offsets in the concentration values. The entire cycle of standard and sample sequences for each dip tube takes less than 20 min.

Figure 7Measurement procedure of a single dip tube sample as recorded in the pressure time-series inside the absorption cell. The sample is flanked by several standard gases in a bracketing approach. The shaded blocks indicate measurement periods. The “shoulder” of approx. 1.7 hPa before every standard gas (STD) period indicates a 30 s purge to equilibrate all surfaces before the cell is filled to a specified pressure.


Each measurement block lasts for 80 s, with the signal being averaged over the entire period. These raw data are calibrated using previously determined calibration curves (according to Bereiter et al.2020) and are denoted as [X]cal. The next step is to apply a drift correction to the sample measurement using the value of the immediately following standard measurements. The drift-corrected values of the sample concentrations [X]SA,cal,pdc (with SA for sample and pdc for pointwise drift correction), using the calibrated concentrations of the standard no. 6 pressure matched measurements (ST6), are calculated in the following way:

(1) [ X ] SA , cal , pdc = [ X ] SA , cal + [ X ] ST 6 , ref - [ X ] ST 6 , cal ,

where [X] stands for the CO2, CH4, and N2O concentrations, and ST6,ref refers to the externally established reference concentrations of standard no. 6 from Table 1.

For δ13CO2, the correction takes a similar form:

(2) δ 13 C SA , cal , pdc = δ 13 C SA , cal + δ 13 C STcm , ref - δ 13 C STcm , cal ,

where STcm stands for the concentration-matched standard gas mixture (the third block in Fig. 7).

3.3 Long-term stability and precision

The long-term stability of the QCLAS system is verified by repeatedly measuring two standard gases (B3 and B4 in Table 1) over the period from September 2021 to March 2022 in approximately 2-week intervals for a total of 49 individual measurements of each gas. These standards are well suited to mimic ice-core-derived ancient air samples, as they are not used in the calibration, making them truly independent targets, and their CO2 concentrations span the range expected in actual samples (see Sect. 2.4). The results are displayed in Fig. 8. For CO2, the drift correction improves the repeatability by about a factor of 2, i.e., from 0.2 ppm to better than 0.1 ppm. In absolute terms, the stated externally calibrated values of B3 and B4 are 0.1 to 0.2 ppm higher than our measurements, which is acceptable for the intended ice core analyses.

Figure 8Repeated biweekly measurements (n=49) of B3 and B4 over 6 months. The two standards differ significantly in all measured concentrations but are essentially identical in δ13CO2 (Table 1). Each measurement set, separated by the dashed gray lines, is composed of three aliquots of each standard. Data are shown as calibrated only (cal) as well as calibrated and drift-corrected (cal and pdc) according to Eqs. (1) and (2). For δ13CO2, drift-corrected values using the concentration-matched standard (cal pdc matched) are shown additionally. On the lower left of each subplot, the standard deviation (1σ) achieved with the different corrections is shown. For δ13CO2, the standard deviation is shown for the combined measurements of both standards. The circled symbol denotes an outlier identified in the secondary (“reference”) channel (more details in the main text). Outliers are not removed from the standard deviation, and hence our precision estimates based on these measurements are conservative. The horizontal magenta lines depict the externally determined values as listed in Table 1.


For CH4, the drift correction has no clear effect on the repeatability. The long-term repeatability is about 3.5 ppb but could be better if not for a few outliers. These outliers can be identified, however, using the secondary beam path in the QCLAS that is directed through a reference cell and monitored with an additional IR detector. (Bereiter et al.2020). As B3 and B4 were not independently referenced for CH4 or N2O, the accuracy of our measurements cannot be assessed with them.

For N2O, the drift correction improves the precision for B3 and B4 to 1 to 1.5  ppb but is again impacted by outliers. Especially for B4, there is a conspicuous bias in the first half of the measurement series whose cause we have not been able to identify.

The largest impact of the pointwise drift correction is on the δ13CO2 values. The repeatability improves from 0.49 ‰ to 0.03 ‰ for the concentration-matched standard. For the two standards B3 and B4, we found within our precision an identical value in δ13CO2 of 8.713 ‰, which is in good agreement with the externally determined isotopic composition from Table 1.

Overall, the repeatability of the measurements, while slightly worse for concentrations than those of Bereiter et al. (2020), is sufficient for the intended ice core analyses. All reported measurements in Sect. 4 include the drift correction.

4 Ice measurements

4.1 Standard over ice experiments

In this section we examine sublimation experiments where calibration gases (see Table 1) are continuously admitted into the extraction vessel while sublimating a piece of gas-free ice (Fig. 9). This mimics the extraction with real ice samples closely and is supposed to reveal anomalies in the studied gas concentrations that arise from the process of sublimation such as (i) production by reflected or stray laser radiation hitting the metal flange, (ii) co-deposition of trace gases with the freezing water vapor on the side of the extraction vessel, (iii) fractionation in the water vapor flux from the sublimating ice sample, (iv) ad- or desorption at the glass and metal walls in the presence of water vapor, and (v) potentially other unknown effects.

Figure 9The standard-over-ice experiments characterize the small offsets that are introduced through the extraction process and gas transfer from the extraction vessel to the MPC. The standards no. 3 (orange) and no. 7 (red) were measured once each, while no. 1 (blue) was repeated five times. Error bars show the long-term reproducibility of discrete standards measured by the QCLAS only, as determined in Fig. 8. For standard no. 1, the error bars are omitted to improve legibility. For CO2, the error bars are smaller than the marker size.


We found that for all target parameters the biases to the expected values are small, with average anomalies of +0.9 ppm for CO2, 0.11 ‰ for δ13CO2, +2.9 ppb for CH4, and +0.3 ppb for N2O. Furthermore, the anomalies are not significantly different for the three standard gases, which span the range of concentrations expected from real ice samples and for δ13CO2 even exceed that. The magnitude of these biases suggests that the aforementioned processes alter the composition of the sample gas only marginally. The standard deviation calculated over all data shown in Fig. 9 amounts to 0.42 ppm, 3.1 ppb, 1.1 ppb, and 0.043 ‰, respectively, for CO2, CH4, N2O, and δ13CO2. These represent a conservative estimate of the precision for a single ice-core-derived measurement.

Considering the CO2 measurements for which a very high precision is achieved, a consistent and reproducible pattern is observed: an elevated concentration in the first dip tube, followed by a lower concentration value in the second dip tube. For δ13CO2, a trend towards more depleted values with increasing sample number was found. It is not clear yet what causes this behavior, but apparently it is inherent in our system. A contamination by insufficient pumping has been ruled out. In Fig. 11 we show results for the Talos Dome ice samples with and without a finger-specific bias correction, which clearly show an improvement of the homogeneity of the results for CO2 and δ13CO2.

When accounting for individual dip tube average biases, the standard deviations of the standard-over-ice measurements improve to 0.29 ppm and 0.031 ‰ for CO2 and δ13CO2. As we are using averages of seven (CO2, CH4, and N2O) and five (δ13CO2) individual experiments, the biases are reasonably confined, but their quantification can be further improved in the future.

4.2 Ice core samples

As the final aim of our development is to measure ice core samples, we used individual EPICA Dome C (EDC) samples of glacial and interglacial ages covering the full range of concentrations expected in Pleistocene ice and one continuous 57 cm long ice section from the Talos Dome ice core to assess the resolution capabilities of the LISE extraction method.

4.2.1 EDC

A total of seven EDC ice core samples were measured with the LISE/QCLAS system, three from a glacial and interglacial period and one additional fully clathrated sample from 1518 m depth (see Table 2 and Fig. 10). For CO2, the measured values agree very well with previous data obtained with different extraction techniques. The variability of 1–2 ppm within a 10 cm long ice sample is noticeably higher than in our standard-over-ice experiments, reflecting true variability in the concentration in the ice, as to be expected from the stochastic nature of bubble trapping at the firn–ice transition. The sample at 768.35 m depth is an outlier in this regard with a much higher intra-sample variability, likely due to the fact that this sample lies in the bubble-to-clathrate transition zone (Neff2014), where layered early clathratization occurs, and differently fast permeation rates between bubbly layers and clathrate layers for different gas species lead to the centimeter-scale variability in the gas composition (Lüthi et al.2010; Oyabu et al.2021).

Table 2Analyzed EDC samples. Each sample was  9 cm long, yielding four to five subsamples. Standard deviations (SDs) are 1σ and given in parentheses.

Download Print Version | Download XLSX

Figure 10Comparison of selected EDC ice core samples (Table 2) with published records. Reconstructed atmospheric CO2 and CH4 data compilations from Köhler et al. (2017). Individual δ13CO2 data points and a smoothing spline are from Eggleston et al. (2016). N2O data from Spahni et al. (2005). All LISE/QCLS data are corrected for the extraction biases established by the standard-over-ice experiments. δ13CO2 is additionally corrected for gravitational fractionation using interpolated δ15N2 data from Haeberli et al. (2021).​​​​​​​


Similarly, also for CH4, the absolute values agree well with the published record. The intra-sample standard deviation is with 4–16 ppb also exceeding the analytical uncertainty determined above. This variability can likely be attributed to differential layer trapping, which produces noise on the centimeter scale (Fourteau et al.2017). Interestingly, the CH4 concentrations do not seem affected by fractionation in the transition from bubbles to clathrates represented by the 768.35 m sample. This agrees with the high-resolution CH4 measurements by Fourteau et al. (2020), where sections from the transition zone likewise do not show inflated variability.

N2O is known to be affected by (variable) in situ production during periods of high dust content (Schilt et al.2010a), which confounds the comparison to other N2O data during glacial maxima. Note that the Talos Dome ice core is least affected by this process due to its lower dust content (Schilt et al.2010b; Fischer et al.2019). Generally, our EDC samples agree well with the published record, but there is a tendency to a positive bias of  7 ppb, even for interglacial samples not affected by in situ production. A similarly high bias compared to previous literature values has been observed recently by Schmidely et al. (2021) using a refined, more efficient melt extraction technique and may point to non-quantitative extraction of highly soluble N2O from the meltwater in older studies. The 100 % extraction efficiency of our sublimation method avoids any such solubility issues. The centimeter-scale variability also observed in interglacial samples is elevated compared to the established analytical performance as for CO2 and CH4 and likely for the same reasons.

The δ13CO2 has the worst signal-to-noise ratio of our measured parameters, complicating the comparison to existing records. In addition, δ13CO2 needs to be corrected for gravitational fractionation in the firn column (Craig et al.1988; Sowers et al.1989; Schwander et al.1993) to be able to compare reconstructed atmospheric values from different ice cores. For our samples, the correction is done using interpolated δ15N2 data from Haeberli et al. (2021). After this correction, the atmospheric δ13CO2 values agree well with previous records. The variability is with 0.02 ‰ to 0.07 ‰ only slightly elevated with respect to the analytical precision, except for the sample from the bubble-clathrate transition zone (768.35 m depth, 0.11 ‰ SD).

4.2.2 Talos Dome

A second set of validation measurements was performed using ice core from Talos Dome (TALDICE; Urbini et al.2006). A 57 cm long section from 1548 m depth, less than 100 m above bedrock, was analyzed continuously, leading to 27 individual data points over the entire section. Although initially undated, the deep TALDICE ice core below 1438 m has recently been re-examined using δ18Oatm, δD, and 81Kr combined with an inverse ice chronology model to produce the TALDICE-deep1 chronology (Crotti et al.2021). This close to bedrock, the stratigraphic layers are highly thinned. The thinning leads to the climatic signal getting compressed, and atmospheric variability is expected to occur on the decimeter scale. The deep TALDICE ice is thus similar in terms of thinning to the BE-OI core currently being drilled at Little Dome C, where a highly compressed stratigraphy similarly requires high-resolution analysis techniques to fully recover the paleo-atmospheric signal (Fischer et al.2013). Therefore, these TALDICE ice samples provide a benchmark test for the new LISE extraction, though the stratigraphy expected for the deep BE-OIC will be even more compressed.

Figure 11A 57 cm section of the Talos Dome ice core analyzed in  2 cm resolution with the LISE/QCLS system described. Faded colors show the raw data (with calibration and drift correction), while solid colors denote the same measurements but with a dip-tube-specific correction for average offsets deduced from the standard-over-ice experiments. Vertical gray lines separate individual,  10 cm long ice samples, which were sublimated yielding typically five subsamples. The shorter sample with only two subsamples is due to a visible fracture at 1548.25 cm depth.


The results (Fig. 11) confirm that the studied section contains carbon dioxide jump (CDJ) 10a, where CO2 concentrations increase from 210 to 219 ppm in less than 250 years, as previously documented in the EDC ice core (Nehrbass-Ahles et al.2020). In our study, however, this short jump is likely resolved by 15 individual data points (see discussion on age scale below), showing a relatively constant rise rate over this interval except for one slightly elevated sample after the first quarter of the jump. For both CH4 and N2O, there is variability on short timescales that was previously not resolved with very rapid swings in CH4 concentrations of 50 ppb over 5 cm, similar in size to the CH4 changes observed for some of the smaller Dansgaard–Oeschger events during Marine Isotope Stage 3 (Baumgartner et al.2014; Rhodes et al.2015). The events in CH4 are also observed in N2O, although more gradual in nature, as expected from the atmospheric lifetime of N2O that is more than 10 times longer and with an earlier onset of the N2O increase, as previously observed for Dansgaard–Oeschger events (Schilt et al.2013). Signs of in situ production such as high uncorrelated scatter of the measurements are not observed during this 57 cm sequence.

Regarding possible corrections for the dip-tube-specific biases documented in the standard-over-ice experiments (Sect. 4.1), there is strong evidence in the TALDICE measurements that indeed both CO2 and δ13CO2 are affected by the same patterns that were established in Sect. 4.1. For CO2, the first measurement of every subset (separated by vertical lines in Fig. 11) is elevated compared to the following measurements. After subtracting the established standard-over-ice biases, this feature disappears. Similarly for δ13CO2, there is a decreasing trend in each subset in the uncorrected data, which is also the case in the standard-over-ice experiments. After applying the standard-over-ice correction, these artifactual trends are removed, leaving a gradual enrichment over the entire 57 cm section. Accounting for a linear trend over the entire 57 cm sequence, the remaining variability amounts to 0.05 ‰ (1σ), just slightly more than the expected statistical noise due to the measurement precision.

Figure 12Comparing our high-resolution Talos Dome data to EDC records. EDC CO2 and CH4 from Nehrbass-Ahles et al. (2020) including carbon dioxide jump 10a and N2O (Spahni et al.2005) are plotted on the AICC2012 age scale (Veres et al.2013), which is also the basis for TALDICE-deep1 (Crotti et al.2021). The Talos Dome data are shown both on the TALDICE-deep1 age scale and on an arbitrary scaling that optimizes the fit to the EDC record.


Comparison with previous studies is complicated due to the large discrepancy in vertical resolution and thus time (Fig. 12). Only CO2 has a reference record from EDC (Nehrbass-Ahles et al.2020) that has enough resolution to allow for a direct comparison. For CH4 and N2O, the absolute concentrations broadly agree, but due to the additional variability resolved in our new data as well as the uncertainties in the age scale (see below), it is difficult to make any further interpretation. It does appear that, just like for EDC ice, our N2O concentrations are again slightly elevated compared to the published data. δ13CO2 has not previously been measured in this time interval at all.

We point out that on the TALDICE-deep1 age scale (Crotti et al.2021), our CO2 concentration data do not fit the published record very well. However, there are several caveats to consider: (i) TALDICE-deep1 ends at 1548 m depth, which means the age of our data is unconstrained towards older ages; (ii) the stated uncertainty of TALDICE-deep1 in this age interval is almost 2000 years, and (iii) the tie points at this age are rather poor (their Fig. 3b). All of these points lead us to speculate that the gas age–depth curve of TALDICE-deep1 is too steep and slightly too old at this depth compared to EDC on AICC2012 (Veres et al.2013). To match our CO2 data to the EDC record, we need to flatten the TALDICE age–depth curve (by manually adjusting tie points), which also aligns better with the 81Kr-dated horizons (their Fig. 6b).

5 Conclusions

We described a coupled extraction and analysis system capable of simultaneous, high-precision measurements of CO2, CH4, and N2O concentrations, as well as δ13CO2, on ice core samples. Compared to existing techniques, our LISE/QCLAS system has the following advantages: (i) no preprocessing of the sample (separation of individual trace gases) is necessary as the QCLAS measures absorption spectra directly in whole air; (ii) the measurement of the carbon isotopic composition of CO2 and the greenhouse gas concentrations takes place on the same sample of extracted air without isobaric interferences common to mass spectrometry, eliminating the need for splitting the sample; (iii) sublimation achieves 100 % extraction efficiency, which is especially crucial for deeper, semi-clathrated and clathrated ice; (iv) controlled sublimation from top to bottom of a sample reaches a vertical resolution (1–2 cm) needed to resolve atmospheric variability in highly thinned ice; (v) equal treatment of ice samples and standard gas over gas-free ice allows us to characterize the entire system and correct for alterations to the studied gases due to the extraction procedure; (vi) there is the possibility to reuse the sample gas due to the non-destructive nature of the QCLAS measurement. This is especially helpful for the gravitational correction of δ13CO2, as δ15N2 could be measured on the same sample. This recycling of sample gas will be tested in the near future. (vii) The analysis is unaffected by drilling fluid contamination of the sample gas, and last but not least (viii) the precision of CO2 measurements but also of δ13CO2 for such small samples sizes (15 g of ice) is significantly improved over previous approaches. While some of these qualities were achieved by previous measurement techniques, the combination of all of them is unique and makes the LISE/QCLAS approach especially suited for the analysis of deep, highly thinned ice as expected from the BE-OI drilling at Little Dome C.

Appendix A: Additional photos

Figure A1Custom aluminum vacuum manifold, termed Hexatubus, with six 1/4 in. dip tubes sitting in the cryostat and connected to the extraction line at the top.


Data availability

The data used in this paper are available from the corresponding author upon request.

Author contributions

LM, DB, and FK designed and developed the LISE and the sample and standard introduction system for the QCLAS, under the guidance of and with input from HF and JS. RW and CR provided mechanical and electronic support. LM, DB, and FK tested the system and performed the measurements. BB developed the QCLAS with input by BT, LE, JS, and HF. HF managed and supervised the project. DB and LM prepared the manuscript with contributions from all authors.

Competing interests

The contact author has declared that none of the authors has any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


This study is part of the ERC Advanced Grant “deepSLice” (667507) to Hubertus Fischer. We gratefully acknowledge the reference measurements of the standard gases by Peter Nyfeler and Markus Leuenberger at the University of Bern.

Financial support

This research has been supported by the H2020 European Research Council (grant no. 667507) and the Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (grant no. 200020B_200328).

Review statement

This paper was edited by Thomas Röckmann and reviewed by two anonymous referees.


Baumgartner, M., Kindler, P., Eicher, O., Floch, G., Schilt, A., Schwander, J., Spahni, R., Capron, E., Chappellaz, J., Leuenberger, M., Fischer, H., and Stocker, T. F.: NGRIP CH4 concentration from 120 to 10 kyr before present and its relation to a δ15N temperature reconstruction from the same ice core, Clim. Past, 10, 903–920,, 2014. a, b

Bauska, T. K., Brook, E. J., Mix, A. C., and Ross, A.: High-precision dual-inlet IRMS measurements of the stable isotopes of CO2 and the N2O / CO2 ratio from polar ice core samples, Atmos. Meas. Tech., 7, 3825–3837,, 2014. a

Bauska, T. K., Brook, E. J., Marcott, S. A., Baggenstos, D., Shackleton, S., Severinghaus, J. P., and Petrenko, V. V.: Controls on Millennial-Scale Atmospheric CO2 Variability During the Last Glacial Period, Geophys. Res. Lett., 45, 7731–7740,, 2018. a, b

Bereiter, B., Stocker, T. F., and Fischer, H.: A centrifugal ice microtome for measurements of atmospheric CO2 on air trapped in polar ice cores, Atmos. Meas. Tech., 6, 251–262,, 2013. a

Bereiter, B., Fischer, H., Schwander, J., and Stocker, T. F.: Diffusive equilibration of N2, O2 and CO2 mixing ratios in a 1.5-million-years-old ice core, The Cryosphere, 8, 245–256,, 2014. a

Bereiter, B., Eggleston, S., Schmitt, J., Nehrbass-Ahles, C., Stocker, T. F., Fischer, H., Kipfstuhl, S., and Chappellaz, J.: Revision of the EPICA Dome C CO2 record from 800 to 600 kyr before present, Geophys. Res. Lett., 42, 542–549,, 2015. a, b

Bereiter, B., Tuzson, B., Scheidegger, P., Kupferschmid, A., Looser, H., Mächler, L., Baggenstos, D., Schmitt, J., Fischer, H., and Emmenegger, L.: High-precision laser spectrometer for multiple greenhouse gas analysis in 1 mL air from ice core samples, Atmos. Meas. Tech., 13, 6391–6406,, 2020. a, b, c, d, e, f, g, h, i, j, k

Bergeron, V., Berger, C., and Betterton, M. D.: Controlled Irradiative Formation of Penitentes, Phys. Rev. Lett., 96, 098502,, 2006. a, b

Berisford, D. F., Foster, J. T., Kosberg, J. A., Furst, B. I., Poston, M. J., Daimaru, T., Lang, M., Backman, L., Macias, A., and Hand, K. P.: Erosion of Penitentes Under Experimental Conditions Relevant to Ice-Covered Airless Worlds, J. Geophys. Res.-Planet., 126, e2021JE006955,, 2021. a

Blunier, T. and Brook, E. J.: Timing of millennial-scale climate change in Antarctica and Greenland during the last glacial period, Science, 291, 109–112,, 2001. a

Cathles, L. M., Abbot, D. S., and MacAyeal, D. R.: Intra-surface radiative transfer limits the geographic extent of snow penitents on horizontal snowfields, J. Glaciol., 60, 147–154,, 2014. a

Chappellaz, J., Stowasser, C., Blunier, T., Baslev-Clausen, D., Brook, E. J., Dallmayr, R., Faïn, X., Lee, J. E., Mitchell, L. E., Pascual, O., Romanini, D., Rosen, J., and Schüpbach, S.: High-resolution glacial and deglacial record of atmospheric methane by continuous-flow and laser spectrometer analysis along the NEEM ice core, Clim. Past, 9, 2579–2593,, 2013. a, b

Craig, H., Horibe, Y., and Sowers, T.: Gravitational Separation Of Gases And Isotopes In Polar Ice Caps, Science, 242, 1675–1678,, 1988. a

Crotti, I., Landais, A., Stenni, B., Bazin, L., Parrenin, F., Frezzotti, M., Ritterbusch, F., Lu, Z.-T., Jiang, W., Yang, G.-M., Fourré, E., Orsi, A., Jacob, R., Minster, B., Prié, F., Dreossi, G., and Barbante, C.: An extension of the TALDICE ice core age scale reaching back to MIS 10.1, Quaternary Sci. Rev., 266, 107078,, 2021. a, b, c

Eggleston, S., Schmitt, J., Bereiter, B., Schneider, R., and Fischer, H.: Evolution of the stable carbon isotope composition of atmospheric CO2 over the last glacial cycle, Paleoceanography, 31, 434–452,, 2016. a, b

Feistel, R. and Wagner, W.: Sublimation pressure and sublimation enthalpy of H2O ice Ih between 0 and 273.16 K, Geochim. Cosmochim. Ac., 71, 36–45,, 2007. a

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., and Wilhelms, F.: Where to find 1.5 million yr old ice for the IPICS “Oldest-Ice” ice core, Clim. Past, 9, 2489–2505,, 2013. a, b

Fischer, H., Schmitt, J., Bock, M., Seth, B., Joos, F., Spahni, R., Lienert, S., Battaglia, G., Stocker, B. D., Schilt, A., and Brook, E. J.: N2O changes from the Last Glacial Maximum to the preindustrial – Part 1: Quantitative reconstruction of terrestrial and marine emissions using N2O stable isotopes in ice cores, Biogeosciences, 16, 3997–4021,, 2019. a

Fourteau, K., Faïn, X., Martinerie, P., Landais, A., Ekaykin, A. A., Lipenkov, V. Ya., and Chappellaz, J.: Analytical constraints on layered gas trapping and smoothing of atmospheric variability in ice under low-accumulation conditions, Clim. Past, 13, 1815–1830,, 2017. a

Fourteau, K., Martinerie, P., Faïn, X., Ekaykin, A. A., Chappellaz, J., and Lipenkov, V.: Estimation of gas record alteration in very low-accumulation ice cores, Clim. Past, 16, 503–522,, 2020. a

Fukusako, S.: Thermophysical properties of ice, snow, and sea ice, Int. J. Thermophys., 11, 353–372,, 1990. a

Haeberli, M., Baggenstos, D., Schmitt, J., Grimmer, M., Michel, A., Kellerhals, T., and Fischer, H.: Snapshots of mean ocean temperature over the last 700 000 years using noble gases in the EPICA Dome C ice core, Clim. Past, 17, 843–867,, 2021. a, b

Köhler, P., Nehrbass-Ahles, C., Schmitt, J., Stocker, T. F., and Fischer, H.: A 156 kyr smoothed history of the atmospheric greenhouse gases CO2, CH4, and N2O and their radiative forcing, Earth Syst. Sci. Data, 9, 363–387,, 2017. a

Kossacki, K. J., Markiewicz, W. J., Skorov, Y., and Kömle, N.: Sublimation coefficient of water ice under simulated cometary-like conditions, Planet. Space Sci., 47, 1521–1530,, 1999. a

Lee, J. E., Edwards, J. S., Schmitt, J., Fischer, H., Bock, M., and Brook, E. J.: Excess methane in Greenland ice cores associated with high dust concentrations, Geochim. Cosmochim. Ac., 270, 409–430,, 2020. a

Loulergue, L., Schilt, A., Spahni, R., Masson-Delmotte, V., Blunier, T., Lemieux, B., Barnola, J.-M., Raynaud, D., Stocker, T. F., and Chappellaz, J.: Orbital and millennial-scale features of atmospheric CH4 over the past 800 000 years, Nature, 453, 383–386,, 2008. a

Lüthi, D., Bereiter, B., Stauffer, B., Winkler, R., Schwander, J., Kindler, P., Leuenberger, M., Kipfstuhl, S., Capron, E., Landais, A., Fischer, H., and Stocker, T. F.: CO2 and O2/ N2 variations in and just below the bubble–clathrate transformation zone of Antarctic ice cores, Earth Planet. Sc. Lett., 297, 226–233,, 2010. a, b, c

Marcott, S. A., Bauska, T. K., Buizert, C., Steig, E. J., Rosen, J. L., Cuffey, K. M., Fudge, T. J., Severinghaus, J. P., Ahn, J., Kalk, M. L., McConnell, J. R., Sowers, T., Taylor, K. C., White, J. W. C., and Brook, E. J.: Centennial-scale changes in the global carbon cycle during the last deglaciation, Nature, 514, 616–619,, 2014. a, b

Mitchell, L. E., Brook, E. J., Sowers, T., McConnell, J. R., and Taylor, K.: Multidecadal variability of atmospheric methane, 1000–1800 C.E., J. Geophys. Res.-Biogeo., 116, G02007,, 2011. a, b

Mühl, M., Schmitt, J., Seth, B., Lee, J. E., Edwards, J. S., Brook, E. J., Blunier, T., and Fischer, H.: Excess methane, ethane and propane production in Greenland ice core samples and a first isotopic characterization of excess methane, EGUsphere [preprint],, 2022. a

Neff, P. D.: A review of the brittle ice zone in polar ice cores, Ann. Glaciol., 55, 72–82,, 2014. a

Nehrbass-Ahles, C., Shin, J., Schmitt, J., Bereiter, B., Joos, F., Schilt, A., Schmidely, L., Silva, L., Teste, G., Grilli, R., Chappellaz, J., Hodell, D., Fischer, H., and Stocker, T. F.: Abrupt CO2 release to the atmosphere under glacial and early interglacial climate conditions, Science, 369, 1000–1005,, 2020. a, b, c, d, e

Oyabu, I., Kawamura, K., Kitamura, K., Dallmayr, R., Kitamura, A., Sawada, C., Severinghaus, J. P., Beaudette, R., Orsi, A., Sugawara, S., Ishidoya, S., Dahl-Jensen, D., Goto-Azuma, K., Aoki, S., and Nakazawa, T.: New technique for high-precision, simultaneous measurements of CH4, N2O and CO2 concentrations; isotopic and elemental ratios of N2, O2 and Ar; and total air content in ice cores by wet extraction, Atmos. Meas. Tech., 13, 6703–6731,, 2020. a

Oyabu, I., Kawamura, K., Uchida, T., Fujita, S., Kitamura, K., Hirabayashi, M., Aoki, S., Morimoto, S., Nakazawa, T., Severinghaus, J. P., and Morgan, J. D.: Fractionation of O2/ N2 and Ar / N2 in the Antarctic ice sheet during bubble formation and bubble–clathrate hydrate transition from precise gas measurements of the Dome Fuji ice core, The Cryosphere, 15, 5529–5555,, 2021. a

Raynaud, D. and Whillans, I. M.: Air content of the Byrd core and past changes in the West Antarctic Ice Sheet, Ann. Glaciol., 3, 269–273,, 1982. a

Rhodes, R. H., Brook, E. J., Chiang, J. C. H., Blunier, T., Maselli, O. J., McConnell, J. R., Romanini, D., and Severinghaus, J. P.: Enhanced tropical methane production in response to iceberg discharge in the North Atlantic, Science, 348, 1016–1019,, 2015. a

Rhodes, R. H., Brook, E. J., McConnell, J. R., Blunier, T., Sime, L. C., Faïn, X., and Mulvaney, R.: Atmospheric methane variability: Centennial-scale signals in the Last Glacial Period, Global Biogeochem. Cy., 31, 575–590,, 2017. a, b

Ryu, Y., Ahn, J., and Yang, J.-W.: High-Precision Measurement of N2O Concentration in Ice Cores, Environ. Sci. Technol., 52, 731–738,, 2018. a

Schilt, A., Baumgartner, M., Blunier, T., Schwander, J., Spahni, R., Fischer, H., and Stocker, T. F.: Glacial–interglacial and millennial-scale variations in the atmospheric nitrous oxide concentration during the last 800,000 years, Quaternary Sci. Rev., 29, 182–192,, 2010a. a, b, c

Schilt, A., Baumgartner, M., Schwander, J., Buiron, D., Capron, E., Chappellaz, J., Loulergue, L., Schüpbach, S., Spahni, R., Fischer, H., and Stocker, T. F.: Atmospheric nitrous oxide during the last 140 000 years, Earth Planet. Sc. Lett., 300, 33–43,, 2010b. a

Schilt, A., Baumgartner, M., Eicher, O., Chappellaz, J., Schwander, J., Fischer, H., and Stocker, T.: The response of atmospheric nitrous oxide to climate variations during the last glacial period, Geophys. Res. Lett., 40, 1888–1893,, 2013. a

Schmidely, L., Nehrbass-Ahles, C., Schmitt, J., Han, J., Silva, L., Shin, J., Joos, F., Chappellaz, J., Fischer, H., and Stocker, T. F.: CH4 and N2O fluctuations during the penultimate deglaciation, Clim. Past, 17, 1627–1643,, 2021. a, b

Schmitt, J., Schneider, R., and Fischer, H.: A sublimation technique for high-precision measurements of δ13CO2 and mixing ratios of CO2 and N2O from air trapped in ice cores, Atmos. Meas. Tech., 4, 1445–1461,, 2011. a, b, c, d, e, f

Schmitt, J., Schneider, R., Elsig, J., Leuenberger, D., Lourantou, A., Chappellaz, J., Köhler, P., Joos, F., Stocker, T. F., Leuenberger, M., and Fischer, H.: Carbon Isotope Constraints on the Deglacial CO2 Rise from Ice Cores, Science, 336, 711–714,, 2012. a, b

Schneider, R., Schmitt, J., Köhler, P., Joos, F., and Fischer, H.: A reconstruction of atmospheric carbon dioxide and its stable carbon isotopic composition from the penultimate glacial maximum to the last glacial inception, Clim. Past, 9, 2507–2523,, 2013. a, b, c

Schwander, J., Barnola, J.-M., Andrié, C., Leuenberger, M., Ludin, A., Raynaud, D., and Stauffer, B.: The age of the air in the firn and the ice at Summit, Greenland, J. Geophys. Res.-Atmos., 98, 2831–2838,, 1993. a

Shin, J., Nehrbass-Ahles, C., Grilli, R., Chowdhry Beeman, J., Parrenin, F., Teste, G., Landais, A., Schmidely, L., Silva, L., Schmitt, J., Bereiter, B., Stocker, T. F., Fischer, H., and Chappellaz, J.: Millennial-scale atmospheric CO2 variations during the Marine Isotope Stage 6 period (190–135 ka), Clim. Past, 16, 2203–2219,, 2020. a

Sowers, T., Bender, M., and Raynaud, D.: Elemental and isotopic composition of occluded O2 and N2 in polar ice, J. Geophys. Res.-Atmos., 94, 5137–5150,, 1989. a

Spahni, R., Chappellaz, J., Stocker, T. F., Loulergue, L., Hausammann, G., Kawamura, K., Flückiger, J., Schwander, J., Raynaud, D., Masson-Delmotte, V., and Jouzel, J.: Atmospheric methane and nitrous oxide of the Late Pleistocene from Antarctic ice cores, Science, 310, 1317–1321,, 2005. a, b

Stowasser, C., Buizert, C., Gkinis, V., Chappellaz, J., Schüpbach, S., Bigler, M., Faïn, X., Sperlich, P., Baumgartner, M., Schilt, A., and Blunier, T.: Continuous measurements of methane mixing ratios from ice cores, Atmos. Meas. Tech., 5, 999–1013,, 2012. a, b

Urbini, S., Cafarella, L., Zirizzotti, A., Bianchi, C., Tabacco, I., and Frezzotti, M.: Location of a new ice core site at Talos Dome (East Antarctica), Ann. Geophys.-Italy, 49, 4–5,, 2006. a

Veres, D., Bazin, L., Landais, A., Toyé Mahamadou Kele, H., Lemieux-Dudon, B., Parrenin, F., Martinerie, P., Blayo, E., Blunier, T., Capron, E., Chappellaz, J., Rasmussen, S. O., Severi, M., Svensson, A., Vinther, B., and Wolff, E. W.: The Antarctic ice core chronology (AICC2012): an optimized multi-parameter and multi-site dating approach for the last 120 thousand years, Clim. Past, 9, 1733–1748,, 2013. a, b

Warren, S. G. and Brandt, R. E.: Optical constants of ice from the ultraviolet to the microwave: A revised compilation, J. Geophys. Res.-Atmos., 113, 1–10,, 2008. a

Werner, R. A. and Brand, W. A.: Referencing strategies and techniques in stable isotope ratio analysis, Rapid Commun. Mass Spectrom., 15, 501–519,, 2001. a

Zimmermann, M., Lindlein, N., Voelkel, R., and Weible, K. J.: Microlens laser beam homogenizer: from theory to application, in: Laser Beam Shaping VIII, edited by: Dickey, F. M. and Shealy, D. L., International Society for Optics and Photonics, SPIE, vol. 6663, 9–21,, 2007. a

Short summary
We present a new method to extract the gases from ice cores and measure their greenhouse gas composition. The ice is sublimated continuously with a near-infrared laser, releasing the gases, which are then analyzed on a laser absorption spectrometer. The main advantage over previous efforts is a low effective resolution of 1–2 cm. This capability is crucial for the analysis of highly thinned ice, as expected from ongoing drilling efforts to extend ice core history further back in time.