Articles | Volume 14, issue 6
Atmos. Meas. Tech., 14, 4279–4304, 2021
Atmos. Meas. Tech., 14, 4279–4304, 2021

Research article 09 Jun 2021

Research article | 09 Jun 2021

Simultaneous measurement of δ13C, δ18O and δ17O of atmospheric CO2 – performance assessment of a dual-laser absorption spectrometer

Simultaneous measurement of δ13C, δ18O and δ17O of atmospheric CO2 – performance assessment of a dual-laser absorption spectrometer
Pharahilda M. Steur1, Hubertus A. Scheeren1, Dave D. Nelson2, J. Barry McManus2, and Harro A. J. Meijer1 Pharahilda M. Steur et al.
  • 1Centre for Isotope Research, University of Groningen, Nijenborgh 6, 9747 AG Groningen, the Netherlands
  • 2Aerodyne Research Inc., 45 Manning Road, Billerica, MA 01821-3976, USA

Correspondence: Pharahilda M. Steur (


Using laser absorption spectrometry for the measurement of stable isotopes of atmospheric CO2 instead of the traditional isotope ratio mass spectrometry method decreases sample preparation time significantly, and uncertainties in the measurement accuracy due to CO2 extraction and isobaric interferences are avoided. In this study we present the measurement performance of a new dual-laser instrument developed for the simultaneous measurement of the δ13C, δ18O and δ17O of atmospheric CO2 in discrete air samples, referred to as the Stable Isotopes of CO2 Absorption Spectrometer (SICAS). We compare two different calibration methods: the ratio method, based on the measured isotope ratio and a CO2 mole fraction dependency correction, and the isotopologue method, based on measured isotopologue abundances. Calibration with the ratio method and isotopologue method is based on three different assigned whole-air references calibrated on the VPDB (Vienna Pee Dee Belemnite) and the WMO 2007 (World Meteorological Organization) scale for their stable isotope compositions and their CO2 mole fractions, respectively. An additional quality control tank is included in both methods to follow long-term instrument performance. Measurements of the quality control tank show that the measurement precision and accuracy of both calibration methods is of similar quality for δ13C and δ18O measurements. During one specific measurement period the precision and accuracy of the quality control tank reach WMO compatibility requirements, being 0.01 ‰ for δ13C and 0.05 ‰ for δ18O. Uncertainty contributions of the scale uncertainties of the reference gases add another 0.03 ‰ and 0.05 ‰ to the combined uncertainty of the sample measurements. Hence, reaching WMO compatibility for sample measurements on the SICAS requires reduction of the scale uncertainty of the reference gases used for calibration. An intercomparison of flask samples over a wide range of CO2 mole fractions has been conducted with the Max Planck Institute for Biogeochemistry, resulting in a mean residual of 0.01 ‰ and −0.01 ‰ and a standard deviation of 0.05 ‰ and 0.07 ‰ for the δ13C measurements calibrated using the ratio method and the isotopologue method, respectively. The δ18O could not be compared due to depletion of the δ18O signal in our sample flasks because of storage times being too long. Finally, we evaluate the potential of our Δ17O measurements as a tracer for gross primary production by vegetation through photosynthesis. Here, a measurement precision of <0.01 ‰ would be a prerequisite for capturing seasonal variations in the Δ17O signal. Lowest standard errors for the δ17O and Δ17O of the ratio method and the isotopologue method are 0.02 ‰ and 0.02 ‰ and 0.01 ‰ and 0.02 ‰, respectively. The accuracy results show consequently results that are too enriched for both the δ17O and Δ17O measurements for both methods. This is probably due to the fact that two of our reference gases were not measured directly but were determined indirectly. The ratio method shows residuals ranging from 0.06 ‰ to 0.08 ‰ and from 0.06 ‰ to 0.1 ‰ for the δ17O and Δ17O results, respectively. The isotopologue method shows residuals ranging from 0.04 ‰ to 0.1 ‰ and from 0.05 ‰ to 0.13 ‰ for the δ17O and Δ17O results, respectively. Direct determination of the δ17O of all reference gases would improve the accuracy of the δ17O and thereby of the Δ17O measurements.

1 Introduction

As atmospheric CO2 (atm-CO2) is the most important contributor to anthropogenic global warming, keeping track of its sources and sinks is essential for understanding and predicting the consequences of climate change for natural systems and societies and for assessing and quantifying the possible mitigating measures. The stable isotope (si) composition of atm-CO2 is often used as an additional tool to distinguish between anthropogenic emissions and the influence of the biosphere on varying CO2 mole fractions (Pataki et al.2003; Zhou et al.2005). For this reason, the si composition of atm-CO2 is monitored at a considerable number of atmospheric measurement stations around the globe. Due to the large size of the carbon reservoir of the atmosphere and the long lifetime of CO2 in the atmosphere, the effects of sources and sinks on the atmospheric composition are heavily diluted. Changes in the isotope composition of atm-CO2 are therefore relatively small compared to the actual changes in carbon fluxes (IAEA2002). Hence, current climate change and meteorological research, as well as the monitoring of CO2 emissions, require accurate and precise greenhouse gas measurements that can meet the WMO/GAW (Global Atmosphere Watch Program of the World Meteorological Organization) inter-laboratory compatibility goals of 0.01 ‰ for δ13C and 0.05 ‰ for δ18O of atm-CO2 for the Northern Hemisphere (Crotwell et al.2020).

Traditionally, high-precision stable isotope measurements are done using isotope ratio mass spectrometry (IRMS) (Roeloffzen et al.1991; Trolier et al.1996; Allison and Francey1995), which requires extraction of CO2 from the air sample before a measurement is possible. This is a time-consuming process wherein very strict, 100 % extraction procedures need to be applied to avoid isotope fractionation and to prevent isotope exchange of CO2 molecules with other gases or water. Extraction of CO2 from air is a major contributor to both random and systematic-scale differences between laboratories and thus complicates the comparison of measurements (Wendeberg et al.2013). Further, due to the isobaric interferences of both different CO2 isotopologues and N2O molecules, which are also trapped with the (cryogenic) extraction of CO2 from air, corrections need to be applied for the determination of the δ13C and δ18O values. Due to the mass interference of the 12C17O16O isotopologue with 13C16O16O (and to a lesser extent 13C17O16O and 12C17O17O with 12C18O16O), the δ13C results need a correction (usually referred to as “ion correction”) that builds upon an assumed fixed relation between δ17O and δ18O. This assumed relation has varied in the past (Santrock et al.1985; Allison et al.1995; Assonov and Brenninkmeijer2003; Brand et al.2010), again giving rise to systematic differences (and confusion) between laboratories. Determination of the δ17O of CO2 samples itself using IRMS is extremely complex, due to the mass overlap of the 13C and 17O containing isotopologues, and can only be done using very advanced techniques, restricted to just a few laboratories at the moment (see Adnew et al.2019, and references therein). As the δ17O in addition to the δ18O values in atmospheric CO2 have the potential to be a tracer for gross primary production and anthropogenic emissions (Laskar et al.2016; Luz et al.1999; Koren et al.2019), a less labor-intensive method that would enable all three stable isotopologues of atm-CO2 to be analyzed at a sufficient precision would be an asset.

Optical (infrared) spectroscopy now offers this possibility following strong developments in recent years in Fourier-transform infrared (FTIR) spectroscopy and, especially for the laser light sources, to perform isotopologue measurements showing precision close to, or even surpassing, IRMS measurements (Tuzson et al.2008; Vogel et al.2013; McManus et al.2015). The technique was developed in the 1990s to a level where useful isotope signals could be measured, first on pure compounds such as water vapor (Kerstel et al.1999) and soon also directly on CO2 in dry whole-air samples (Becker et al.1992; Murnick and Peer1994; Erdélyi et al.2002; Gagliardi et al.2003). Extraction of CO2 from the air can therefore be avoided, and smaller sample sizes suffice. Finally, optical spectroscopy is truly isotopologue-specific and is thus free of isobaric interferences, hence giving the possibility to directly measure the δ17O in addition to the δ13C and δ18O. Recent studies already showed the effectiveness of optical spectroscopy for the measurement of δ17O in pure CO2 for various applications (Sakai et al.2017; Stoltmann et al.2017; Prokhorov et al.2019).

In this paper we present the performance, in terms of precision and accuracy, of an Aerodyne dual-laser optical spectrometer (CW-IC-TILDAS-D) in use since September 2017, for the simultaneous measurement of δ13C, δ18O and δ17O of atm-CO2, which we refer to as the Stable Isotopes of CO2 Absorption Spectrometer (SICAS). The instrument performance over time is discussed, followed by an analysis of the CO2 mole fraction dependency of the instrument. We report CO2 mole fractions in micromoles per mole (µmol mol−1), also referred to as parts per million (ppm). The actual ways of performing a calibrated measurement using either individual isotopologue measurements or isotope ratios is discussed, and whole-air measurement results of both calibration methods are evaluated for their compatibility with IRMS stable isotope measurements. Conclusively, the usefulness of the triple oxygen isotope measurements for capturing signals of atmospheric CO2 sources and sinks is evaluated.

2 Instrument description

2.1 Instrumental setup

The optical bench as depicted in Fig. 1 includes the two lasers, several mirrors to combine and deflect the laser beams, the optical cell and two detectors. The two interband cascade lasers (Nanoplus GmbH, Germany) operate in the mid-infrared region. The isotopologues that are measured are 12C16O2, 13C16O2, 12C16O18O and 12C16O17O, which from now on will be indicated as 626, 636, 628 and 627 respectively, following the HITRAN database notation (Gordon et al.2017). Application of a small current ramp causes small frequency variations, so the lasers are swept (with a sweep frequency of 1.7 kHz) over a spectral range in which ro-vibrational transitions of the isotopologues occur with similar optical depths (Tuzson et al.2008). Laser 1 operates in the spectral range of 2350 cm−1 (4.25 µm) for measurement of 627 (and 626), and laser 2 operates around 2310 cm−1 (4.33 µm) for the measurement of 626, 636 and 628. The lasers are thermoelectrically cooled and stabilized to temperatures of −1.1 and 9.9 C, respectively. The beams are introduced in a multi-pass aluminum cell with a volume of 0.16 L, in which an air sample is present at low pressure (∼50 mbar). The total path length of the laser light in the optical cell is 36 m.

Figure 1Scheme of the optical board of the SICAS (figure adapted from McManus et al.2015). Two pathways are shown, both consisting of signals from both lasers: the sample measurement beam in red and the reference beam in blue. The reference pathway is in our case only used for fitting purposes. RBS stands for reference beam splitter. One of the detectors is used to read the signal of the sample beam, the other for the reference beam. The red trace laser is co-aligned with the sample path to visualize the sample pathway to ease alignment.

After passing the cell, the lasers are led to a thermoelectrically cooled infrared detector, measuring the signal from the lasers in the spectral range (Fig. 2). The lasers, optical cell and detectors are all in a housing that is continuously flushed with N2 gas to avoid any other absorption by CO2 than from gas in the optical cell. The temperature within the housing is controlled using a recirculating liquid chiller set at a temperature of 20 C to keep the temperature in the cell stable. Within a measurement sequence (∼12 h) the temperature does not normally fluctuate more than 0.05 C. The absorption spectra are derived by the software TDLWintel (McManus et al.2005) that fits the measured signal based on known molecular absorption profiles from the HITRAN database (Rothman et al.2013). On the basis of the integration of the peaks at the specific wavelengths, measured pressure and temperature in the optical cell and the constant path length, the isotopologue mole fractions are calculated by the TDLWintel software with an output frequency of 1 Hz. For convenience, the default output for the isotopologue mole fractions is scaled for the “natural abundances” of the 626, 636, 628 and 627 as defined in Rothman et al. (2013), but for obtaining the raw mole fractions, this scaling is avoided.

Figure 2Typical absorption spectrum, transmission spectrum and residual for laser 1 (a) for measurement of 627 and 626 and laser 2 (b) for measurement of 626, 628 and 636. The residuals show systematic deviations at the line positions. These deviations are primarily due to the use of the Voigt line shape function in the spectral fitting model, rather than a more complex line shape function such as Hartmann–Tran. Careful analysis has shown that the use of the more convenient Voigt line shape function does not add noise, drift or calibration error as implemented in the isotope analyzer.


The gas inlet system, depicted in Fig. 3, is designed to measure discrete air samples in static mode, such that one can quickly switch between measurements of different samples. The system consists of Swagelok stainless steel tubing and connections and pneumatic valves. There are two inlet ports (11 and 14) which are connected to the sample cross at the heart of system (from now on indicated as inlet volume), where a sample is collected at the target pressure of 200±0.25 mbar before it is connected to the optical cell. One of the inlet ports (11) is connected to a 1/8” VICI multi-valve (Valco Instruments) with 15 potential positions for flask samples or cylinders. The cylinders depicted in Fig. 3 will be defined in Sects. 2.2 and 3.2. When the VICI valve switches from position, the volume between port 10 and 9 is flushed seven times with the sample gas to prevent memory effects due to the dead volume of the VICI valve. A sample gas is led into the inlet volume at reduced flow, as a critical orifice is placed right before the inlet valve, while another gas is being measured inside the optical cell. Since the closing and opening of the valves are controlled by the TDLWintel software, it also controls the duration of the flow into the inlet volume. The target pressure is reached using input from a pressure sensor placed inside the inlet volume. After evacuation of the optical cell (opening valve 22 and 23), the gas from the inlet volume can immediately be brought into the optical cell (opening valve 19 and 23), thereby reducing the sample pressure to ∼50 mbar.

Figure 3Gas inlet system of the SICAS with one VICI multi-valve inlet port, connected to three high-pressure natural-air tanks and 12 free ports for samples. The includes an extra inlet port for the working gas, also a high-pressure natural-air tank.


The gas handling procedures are different for measurements of air from cylinders or flasks. For the cylinders, single-stage pressure regulators are in use (Rotarex, model SMT SI220), set at an outlet pressure of ∼600–1000 mbar (absolute). If measurements are started after more than 2 d of inactivity, the internal volume of the regulators is flushed 10 times to prevent fractionation effects. To open and close the flasks, we use a custom-built click-on electromotor valve system (Neubert et al.2004), making it possible to open the flasks automatically before the measurement. Before opening the flask, the volume between valve 9 and the closed flask is evacuated, so there is no need to flush extensively, and less sample gas is lost. The actions described above are all steered by a command program developed by Aerodyne Research Inc. called the Switcher program. A bespoke script writing program developed in FileMaker Pro enables us to quickly write scripts for measurement sequences and to directly link those measurements to an internal database.

Figure 4The Allan variance as a function of the integration time in seconds for a single gas measurement plotted for both the measured isotopologue abundance (a, b) and the isotope ratios (c, d) in September 2017 (a, c) and July 2019 (b, d). The best achieved precision and corresponding integration times are shown as a table in the plots.


2.2 Instrument performance

The SICAS measurement performance was evaluated by determining the Allan variance of the four measured isotopologue abundances and the three isotope ratios as a function of measurement time on a single whole-air sample in the sealed optical cell. The isotope ratios, defined here as the ratio of the rare isotopologue (636, 628 and 627) and the most abundant 626 isotopologue, are r636, r628 and r6271. This experiment was first done in September 2017 and repeated in July 2019 to see the development in time of the measurement precision (Fig. 4). In all cases, drifts outweigh the averaging process after time periods ranging from 16 to 75 s, and this is short compared to the duration of the normal measurement sequences. This is a firm indication that continuous drift correction using gas from a high-pressurized cylinder, of comparable CO2 concentration and isotope composition as atmospheric samples, is necessary for optimal results. The cylinder used for drift correction which we define as the working gas contains natural air, of which the isotope composition and the CO2 concentration is known.

The precision became significantly worse for all species but isotopologue 627 in the time period between September 2017 and July 2019. In this same period a gradual but significant decrease (of about 50 %) in the measured laser intensity was observed. For most species this led to an increase of the optimal integration time, which is logical given the fact that the minimal precision was higher, such that the increase due to drift influences the acquired precision at a higher integration time and also at a higher variance level. Figure 4 shows the rapid variance increase due to drift for all isotopologues after less than 1 min for the September 2017 measurements, and the same happens for the July 2019 measurements, only less visible due to the higher minimal variance levels.

The decreased laser intensity, potentially leading to a deteriorated signal-to-noise ratio, was caused by contamination of the mirrors in the optical cell, most likely due to precipitation of ultra-fine salt-based aerosols from the sample air occurring during evacuation of the cell. The majority of flask samples measured on the SICAS are from the atmospheric measurement station Lutjewad, which is located on the northern coast of the Netherlands in a rural area dominated by cropland and grassland mainly used for dairy cows. The aerosol composition at Lutjewad is therefore expected to be dominated by sea salt and ammonium nitrate from agricultural emissions. Hence, we were able to clean the mirrors and retrieve ∼80 % of the original laser signal by flushing the mirrors with ultrapure water and ethanol (in that order). This procedure, performed on 31 October 2019, deviates from the recommended mirror cleaning instructions in which it is advised to use ethanol only to clean the mirrors. The additional use of ultrapure water was in our case necessary since the precipitated aerosols were not dissolved in ethanol and were therefore not removed when we used ethanol only. Despite the increase of the laser signal due to the cleaning procedure, precision did not improve as a consequence of it. This indicates that other, still unidentified, issues played a role in the decrease of measurement precision.

To reduce short-term instrumental drift, all sample measurements needed to be alternated with measurements of the working gas, as then the drift-corrected signal can be expressed as

(1) I S ( t ) = M S ( t ) M WG ( t ) ,

where M stands for measurement which can be either the measured isotope ratio or isotopologue abundance, S stands for sample, WG stands for working gas, t stands for time of the sample gas measurement and I stands for the index referring to the drift-corrected sample measurement, of either the isotope ratio or the isotopologue abundance. WG(t) is the measured working gas at time t derived from the time-dependent linear regression of the measurements of the working gas bracketing the sample gas measurement. The effectiveness of this drift correction method was tested for the measured isotope ratios only. Although one of the tested calibration methods uses isotopologue abundances for the initial calibration, the isotope composition is expressed as a delta value and will therefore eventually be calculated using isotope ratios (Sect. 3.3.3). The precision of the isotope ratios will therefore always determine the measurement precision.

A tank was measured >10 times alternately with the working gas. The relative standard deviations were calculated for n=5 and n=10, both with and without drift correction (Table 1). It is expected that, if the drift correction is effective, the standard deviations of the uncorrected values are higher than the standard deviations of the corrected values. The drift correction is effective as the standard deviations of the corrected values are always lower than of the uncorrected values. Although the drift correction procedure is not perfect, as we see a small increase of the standard deviation between n=5 and n=10 between 0.005 ‰ and 0.012 ‰, we can still conclude that the drift correction will result in a better repeatability of the isotope ratios.

Table 1Relative standard deviations for n=5 and n=10 of uncorrected (uncor) and corrected (cor) isotope ratio sample measurements. Sample measurements were always bracketed by measurements of the working gas. Standard deviations of the uncorrected measurements only use the sample measurements; standard deviations of the corrected measurements use drift-corrected (Eq. 1) sample measurements using the working gas measurements.

Download Print Version | Download XLSX

Cross-contamination, being the dilution of a small volume of the working gas in the sample aliquot that is being measured, and vice versa, as described for a dual-inlet IRMS in Meijer et al. (2000), will occur in the SICAS due to the continuous switching between sample and machine working gas. If cross-contamination is not corrected for, DI-IRMS measurement inaccuracy can occur when samples of a highly deviating isotope composition are measured. On the SICAS, only atmospheric samples are measured that are of very similar isotope values. The CO2 mole fraction of the samples can deviate quite strongly from the machine working gas, so effects of cross-contamination will have an influence on the CO2 mole fraction in the optical cell. From experimental data, we quantified the fraction of the preceding sample that affects a sample measurement to be max 0.01 %. A sensitivity analysis was performed using this fraction and showed that this is such a small amount that scale effects due to cross-contamination are well below the precision found in this study (for a detailed description of the analysis, see Appendix E). If samples of CO2 concentrations outside the range of atmospheric samples are measured, it will be essential to also take into account the surface adsorption effects of the aluminum cell which is known to adsorb CO2 (Leuenberger et al.2015). CO2 adsorption in the cell of the SICAS was clearly visible as a drop of measured CO2 concentration when an atmospheric sample was let into the cell right after the cell was flushed with a CO2-free flush gas (hence stripped from CO2 molecules sticking to the cell surface).

3 Calibration experiments

3.1 The CO2 mole fraction dependency

The stable isotope composition of atmospheric CO2 is expressed as a delta value on the VPDB (Vienna Pee Dee Belemnite) (13C)/VPDB-CO2 (17O and 18O) scales, which are realized by producing CO2 gas (using phosphoric acid under well-defined circumstances) from the IAEA-603 marble primary reference material (successor to the now obsolete NBS-19) (IAEA2016). A complication when compared to classical DI-IRMS isotope measurements (or to optical measurements of pure CO2 for that matter) is that in the practice of laser absorption spectroscopy the mole fraction of CO2 in a gas affects the measured stable isotope ratios (and thus delta values) of CO2. Quantification, let alone elimination of this CO2 mole fraction dependence (CMFD), is difficult (McManus et al.2015), but two sources of CMFD were identified by Wen et al. (2013) and related to different calibration strategies. In the first place, CMFD results from nonideal fitting of the absorption spectra, which will to some extent always occur. Capturing the true absorption spectrum is very complicated, due to, among others, line broadening effects of the various components of the air, far wing overlap of distant but strong absorptions, temperature and pressure variability and the choice of line shape function (see Fig. 2 and caption). Secondly, a more “trivial” CMFD is introduced when calibration is done on measured isotopologue ratios, and the intercept of the relation between the isotopologues and the CO2 mole fraction is nonzero (Griffith et al.2012). This effect can be explained by expressing the calculation of the isotopologue ratio by

(2) r * = X * X 626 ,

in which X* is the measured isotopologue mole fraction, and * indicates which of the rare isotopologues is used. When the relation of the measured isotopologue mole fraction and the CO2 mole fraction is linear, this can be described by

(3) X * = X CO 2 α + β .

When Eq. (3) is brought into 2 for either or both of the rare and the abundant isotopologue mole fraction, and β is nonzero for one of those, this leads to an approximate inverse dependence of the measured ratios on the concentration (Griffith et al.2012).

3.1.1 Experiment description

Three experiments have been conducted over the last 2 years to determine the CMFD and to assess its stability over time. These experiments were conducted in December 2017 (experiment 1), in December 2018 (experiment 2) and in May 2019 (experiment 3). Experiment 1 was conducted in cooperation with the Institute for Marine and Atmospheric research Utrecht (IMAU) and served as the initial determination of the CMFD on the SICAS. Experiments 2 and 3 were meant to assess the stability of the CMFD over time. A methodology to determine the CMFD of the r636 for a comparable dual-laser instrument has been described by McManus et al. (2015). In their study, a pure CO2 working gas was diluted back to different CO2 mole fractions using a setup including computer-controlled valves connected to a flow of air without CO2 (“zero air”). CO2 and zero-air mixtures were led directly into the continuous-flow dual-laser instrument. In this way, it was possible to measure the CMFD over a wide range of CO2 mole fractions, from ∼0 to 1000 ppm. The CMFD correction function for the isotope ratios was derived by applying a fourth-order polynomial fit to these measurements.

For determination of the CMFD on the SICAS, this approach was used with some adjustments. The SICAS is designed for the measurement of atmospheric samples of which the relevant range of CO2 mole fractions is ∼370–500 ppm, and CMFD experiments were therefore for the most part conducted in this range. The SICAS measures discrete air samples; hence air mixtures were manually prepared in sample flasks by back-diluting a well-known pure CO2 in-house reference gas to different CO2 mole fractions in the ambient range. Air samples for experiment 1 were prepared at the IMAU, Utrecht University.

Air samples for experiments 2 and 3 were prepared manually in our own laboratory; for the detailed procedure, see Appendix A. The dilutor gas consists of natural air scrubbed of CO2 and H2O using Ascarite® (sodium-hydroxide-coated silica, Sigma-Aldrich) and Sicapent® (phosphoric anhydride, phosphorus(V) oxide), which results in dry, CO2-free natural air. For experiment 2, additional samples were prepared using synthetic air mixtures with and without 1 % Argon as dilutor gas for evaluation of the effect of air composition on the CMFD (see also Sect. 3.1.6). With our manual preparation system, we were able to prepare 10, 12 (with the dilutor being whole air) and 7 flasks for experiments 1, 2 and 3 respectively, that were within our relevant range of atmospheric CO2 mole fractions. McManus et al. (2015) applied a polynomial curve fit on the isotope ratio as a function of the CO2 mole fraction. In this study we focus on a narrower range of CO2 mole fractions, and therefore we expect that a linear or quadratic relationship is sufficient to describe the measured ratios as a function of the CO2 mole fraction. We therefore considered the lower number of samples that were used for the three experiments in comparison to the continuous-flow experiment by McManus et al. (2015) to be sufficient. Griffith (2018) showed that a combination of a linear and inverse relationship to the CO2 mole fraction is theoretically expected, and this relationship fitted the data used in his study in practice. As we expect to have a relation of the measured delta values and the CO2 mole fraction which is close to linear, we use a quadratic relation which approximates this expected theoretical relation closely.

In the next two paragraphs we will discuss the results of the above-described experiments for evaluation of the two sources of CMFD according to Wen et al. (2013) for the SICAS.

3.1.2 Spectroscopic nonlinearities of measured isotopologues

The first source described by Wen et al. (2013), nonlinearity of the relation between the measured isotopologue mole fraction and the CO2 mole fraction, is determined by analysis of the linear fits of the measured rare isotopologue mole fractions (X636, X628 and X627) as a function of the measured X626. We used the IS(t) from Eq. (1) for both the rare isotopologue and the abundant isotopologue mole fractions. The 626 mole fraction is calculated by multiplying I626(t)dc by the known 626 mole fraction of the working gas. The residuals of the linear fits are manipulated such that residuals of the lowest mole fractions are zero (Fig. 5). A linear relation would result in residuals scattering around zero, without a pattern, while systematic nonlinearities would result in a significant pattern, recurring for the different experiments. From the results in Fig. 5 we can conclude that nonlinearities occur; however, these are only clearly visible in experiment 1 for the X636 and the X627 isotopologue and to a lesser degree in experiment 2 for the X636 isotopologue. The maximum residuals of both the X636 and the X627 are highest in experiment 1, which is also the experiment covering the highest range of CO2 mole fractions. From these experiments we can therefore conclude that nonlinearities of the measured rare isotopologue mole fractions and the X626 isotopologue occur but are only significant if the range of CO2 mole fraction is higher than 100 ppm. For the X628 we do not see significant nonlinearities, even if the CO2 mole fraction is much higher than 100 ppm. The maximum residuals of the X628 are not influenced by the CO2 mole fraction, and we therefore conclude that nonlinearities are below the level of detection in these experiments.

Figure 5Residuals (expressed in ‰ relative to the measured amount fraction) of the linear fit of the rare isotopologue abundances as a function of the X626 and the quadratic fit on the residuals. From top to bottom: experiment 1, experiment 2 and experiment 3. The colors red, dark blue and light blue are used for the isotopologues 636, 628 and 627 respectively. Error bars are the combined standard deviations of the 626 and rare isotopologue measurements. Per isotopologue the R2 of the quadratic fit on the residuals is indicated in the tables on the right, as well as the maximum residual (in ‰) on the linear fit of the rare isotopologue as a function of the X626.


Figure 6Measured δ636 of three experiments; black points are experiment 1, red points are experiment 2 and green points are experiment 3.


3.1.3 Introduced dependency on measured delta values

The second source for CMFD, described by Wen et al. (2013), is the introduced dependency on measured isotope ratios if intercepts of the different isotopologues of the analyzer's signal are nonzero, or as in our case for some experiments, if different isotopologues of the analyzer's signal are nonlinear in a different way. In this paragraph we look into the different possibilities to correct for the CMFD of the measured deltas based on observations of the experiments that were described in the section above.

Isotope ratios are susceptible to instrumental drift, but delta values are drift-corrected as the uncalibrated delta value δS is calculated by

(4) δ S * = r S ( t ) * r WG ( t ) * - 1 ,

where S(t) and WG(t) stand for sample and working gas at the time of the sample measurement, respectively and * stands for the rare isotopologue of which the delta is calculated. The rWG(t)* is calculated using the same method as MWG(t) is calculated in Eq. (1). The CMFDs of the deltas are determined by conducting a linear fit on the measured delta values as a function of the measured CO2 mole fraction.

Table 2Slopes derived from the linear fits of the three measured deltas and CO2 mole fractions and the standard errors of the slopes. Delta values calculated with Eq. (4).

Download Print Version | Download XLSX

Table 3Mean residuals for correction of the CMFD of the three measured deltas using three different scenarios; lin and q are calculated relations, using the linear and quadratic fit, respectively, of the rare isotopologue as a function of the abundant isotopologue. Fit and fit q are the linear and quadratic fit, respectively, of the measured delta values as a function of the CO2 mole fraction. The minimum and maximum CO2 mole fractions that were used per experiment are shown in the first column.

Download Print Version | Download XLSX

The results for δ636 are shown in Fig. 6, and slopes of all deltas and the standard errors of the slopes are shown in Table 2. Note that in some cases the standard error of the slope is close to the slope itself, and it is therefore questionable whether a significant CMFD is measured at all. As the CO2 used for the different experiments was not of similar isotope composition, the δ636 measurements in Fig. 6 were normalized such that at the CO2 mole fraction of 400 ppm all ratios are 1. Only the calculated slope is therefore of importance when considering the CMFD of the different experiments. From Table 2 it is clear that the δ636 shows the strongest CMFD. The results show that the CMFD varies for the three different experiments for all measured deltas. Changing instrumental conditions can be an explanation for this change in the CMFD. A drop in measured laser intensity, for instance, was observed over the period between experiment 1 and experiment 3. We should, however, also consider the different range of CO2 mole fractions of the different experiments.

Although most of the variance occurring in the observed CMFD of the deltas (especially of the δ636) can be explained by the linear relationship we found with the measured CO2 mole fraction, we can, from the observed nonlinearities of the measured isotopologues, expect that these relations are better explained by a polynomial relation. We compare therefore both linear and quadratic fits of the measured deltas with calculated relations derived from the fits of the rare isotopologues as a function of the measured 626 isotopologue mole fraction. The theoretically expected combination of a linear and inverse relationship as described in Griffith (2018) showed very similar results as the quadratic fit results, so we consider the quadratic fit to be a good approximation of the theoretically expected relationship. Two relations are calculated: assuming a linear dependency of the rare isotopologue on the abundant isotopologue and assuming a quadratic dependency of the rare isotopologue on the abundant isotopologue. To compare all four scenarios (assuming a linear or quadratic CMFD of the measured deltas and calculation of the CMFD of the deltas assuming a linear and a quadratic dependency of the rare isotopologues on the abundant isotopologue), the mean of the absolute residuals of the observations was calculated for all three experiments and is shown in Table 3. The quadratic fit of the deltas (fit q) shows the lowest mean residuals (except the δ13C in experiment 3), followed by the calculated relation of the deltas when using a quadratic relation of the individual isotopologues and the CO2 mole fraction (q). From these results it can therefore be concluded that determination of the quadratic CMFD of the deltas will give the most accurate results in most cases. It is, however, the question whether this is feasible in practice, as we also know that the CMFD can change through time due to changing instrumental conditions. Determination of a (accurate) quadratic relation requires at least three measurement points (but preferably more) of atm-CO2 of the same isotope composition. In our lab, CO2 in air samples of the same isotope composition but deviating CO2 mole fractions are prepared manually, again introducing uncertainties, and doing these experiments regularly is therefore labor- and time-intensive. Note as well that the range of the CO2 mole fractions in the three experiments is quite high, considering the range of CO2 mole fractions in atmospheric samples. The differences between the four scenarios are significantly smaller in experiment 3 (covering 96 ppm) than in experiment 1 (covering 621 ppm). In the daily procedure of the SICAS, there are at least two CO2-in-air reference gases (in short “reference gases”), high pressurized cylinders containing gas of known isotope composition and CO2 mole fraction, measured bracketing most of the CO2 mole fractions (covering 82 ppm) that occur in atmospheric samples. As all sample and reference measurements are divided by measurements of the working gas when the delta values are calculated, the measured delta value of the working gas should always be zero. The two reference cylinders, together with the zero point for the working gas, provide us with three points to determine a quadratic CMFD of the measured deltas. In this way it is possible to apply a quadratic CMFD correction on the measured deltas. It should be noted that tests showed that the improvements of a quadratic fit (in this form) compared to a linear fit were very small within the narrow range of CO2 mole fractions occurring in the atmosphere, in line with the results of Table 3. However, when samples of very deviating CO2 mole fractions are measured, a quadratic fit will certainly improve the accuracy of the measurement.

3.2 Standard materials and reference scales

Four high-pressure gas tanks (40 L Luxfer aluminum, alloy 6061, max. pressure of 200 bar) containing reference gases are used in the daily measurement procedure of the SICAS: a working gas used for drift correction and possibly for a first calibration step, a quality control tank that is being treated as a sample and two tanks containing a high mole fraction reference gas and a low mole fraction reference gas, from now defined as the high reference and the low reference, which can thus be used for CMFD corrections. The high and low reference cover a great part of the CO2 mole fraction range occurring in atmospheric samples.

It is known that for laser spectroscopy the composition of the sample air affects the absorption line profiles by pressure broadening effects (“matrix effects”), with non-negligible consequences (Nakamichi et al.2006; Nara et al.2012; Harris et al.2020). Hence, it is likely that air composition affects CO2 isotope measurements for the SICAS as well. The possible effects of air composition on the CMFD have been tested by measurement of samples of the same CO2, mixed to different CO2 mole fractions, prepared according to the method described in Sect. 3.1.1 and Appendix A, using three different dilutor gases. The gases that have been used in addition to the CO2-free natural air (whole air) were synthetic air (20 % O2 and 80 % N2, purity is >=99.99 %) and the same synthetic air with addition of 1% of Argon, both prepared by Linde Gas. Linear fits on the measured r636 as a function of the CO2 mole fraction show a small but significant difference of the resulting slopes of 0.0014 ‰ ppm−1 (Table 4) between the synthetic air and whole-air samples. For the r628 and r627, the slope was much smaller, and the standard error of the slope was too large to determine a significant difference between the use of the synthetic dilutors and whole air. Nevertheless, to avoid inaccuracy due to a different CMFD of r636 of samples and references, we solely use gas consisting of natural, dried air as then the effects of the (very small) variability in air composition are negligible.

Table 4CMFD for samples of the same CO2 diluted back with different dilutors. Per dilutor the slopes, resulting from the linear fits of measured r636 and 626 isotopologue mole fraction (ppm), and the standard errors of the slopes are indicated.

Download Print Version | Download XLSX

The gas tanks were produced in-house from dry, compressed natural air collected on the roof of our institute using a RIX compressor (model SA-3). The high and low reference were produced as follows: the high-reference cylinder was filled up to ∼150 bar in winter at 15 January 2018, so the resulting CO2 mole fraction is relatively high (423.77±0.07 ppm). The low-reference cylinder was subsequently produced by transferring air from the high-reference cylinder to an empty cylinder, using the pressure difference, while completely removing CO2 from the air as it flew through a tube filled with Ascarite®. After the low-reference cylinder was filled up to ∼13 bar with CO2-free air, the Ascarite® filled tube was removed, and the filling was continued until the pressure of both cylinders was ∼70 bar. In this way the CO2 mole fraction of the low-reference cylinder was reduced in comparison with the high-reference cylinder, without influencing the CO2 isotope ratios. The resulting CO2 mole fraction of the low reference was 342.81±0.01 ppm. A scheme of the whole setup and detailed description of the procedure can be found in Appendix B.

The CO2 mole fraction of the tanks was measured on a PICARRO G2401 gas mole fraction analyzer and calibrated using in-house working standards, linked to the WMO 2007 scale for CO2 with a suite of four primary standards provided by the Earth System Research Laboratory (ESRL) of the National Oceanic and Atmosphere Administration (NOAA). The uncertainty of the WMO 2007 scale was estimated to be 0.07 µmol mol−1. The typical measurement precision of the PICARRO G2401 measurements is 0.01 µmol mol−1, resulting in a combined uncertainty of 0.07 µmol mol−1 for the assigned CO2 mole fraction values of the calibration tanks, while difference between the two cylinders is known with a much lower uncertainty. The PICARRO analysis is based on the 626 isotopologue mole fraction, not on whole CO2. This is a potential source of error if the isotope composition of different reference gases varies significantly. As the isotope compositions of the used reference gases are close (see Table 5), the variation is not significant for this error (Griffith2018).

Table 5Calibrated whole-air working standards used in daily operation of the SICAS measurements. CO2 measurements were conducted in our lab on a PICARRO G2401 gas mole fraction analyzer, and the δ13C and δ18O values were measured at the MPI-BGC with a MAT-252 dual-inlet IRMS. The δ17O values of the working gas and the quality control tank were measured at the IMAU, while the δ17O values of the low and high references were indirectly determined using our own measurements on the SICAS. Errors are all combined errors, including measurement precision, measurement accuracy and scale uncertainty.

Download Print Version | Download XLSX

Aliquots of all four tanks have been analyzed at the MPI-BGC in Jena by IRMS to link the δ13C and δ18O directly to the JRAS-06 scale (Jena Reference Air Set for isotope measurements of CO2 in air (VPDB/VPDB-CO2 scale)) (Wendeberg et al.2013). The JRAS-06 scale uses calcites mixed into CO2-free whole air to link isotope measurements of atm-CO2 to the VPDB scale. An overview of our reference gases measured at the MPI-BGC and their final propagated error is presented in Table 5, and it can be seen that the low and high reference are very close in isotope composition but seem to differ slightly in their δ13C composition (by 0.05 ‰).

Aliquots of the working gas and quality control gas were analyzed for their δ18O and δ17O values at the IMAU in Utrecht. These values were related to the VSMOW (Vienna Standard Mean Ocean Water) scale using two pure in-house reference gases. The δ17O values are converted to the VPDB-CO2 scale using the known relations between the reference materials VSMOW and VPDB. As the low reference and high reference were not measured at the IMAU, the δ17O values were calculated from experimental results in which a linear CMFD correction was conducted using the measured δ17S (as in Eq. 4) of the low and high reference, assuming that the δ17O values of both gases are similar. Subsequently another linear fit is conducted on the CMFD corrected δ17O values using the known values of the working gas and quality control gas, deriving the calibrated δ17O values of the low reference and high reference. Note that for measurement of our reference gases by the MPI-BGC and IMAU, aliquots were prepared using the “sausage” method, meaning that several (in this case five) flasks are connected and flushed with the sample gas, resulting in a similar air sample in all flasks. However, deviations of the sampled air and the air in reference cylinders due to small leakages or other gas handling problems might be introduced.

3.3 Calibration methods

Two different calibration strategies are discussed in this section. The calibration strategies are based on the two main approaches for calibration of isotope measurements, as also described by Griffith et al. (2012) and, more recently, by Griffith (2018), being (1) to determine the isotopologue ratios, and calibrate them, taking the introduced CMFD into account, from now on defined as the ratio method (RM), and (2) to first calibrate the absolute isotopologue mole fractions individually and then calculate the isotopologue ratios, from now on defined as the isotopologue method (IM). We give a brief introduction of the two calibration methods, as described in literature, and we describe the measurement procedure that is used for both calibration methods. This section ends with a detailed description of both methods as applied for the SICAS measurements.

The RM, being very similar to calibration strategies applied by isotope measurements using DI-IRMS (Meijer2009), is usually based on reference gases covering delta values of a range which is similar to the range of the measured samples. Determination of the CMFD can be done by measuring different tanks of varying CO2 mole fractions or by dynamical dilution of pure CO2 with CO2-free air (Braden-Behrens et al.2017; Sturm et al.2012; Griffith et al.2012; McManus et al.2015; Tuzson et al.2008), again covering the CO2 mole fraction range of the measured samples.

The IM has the advantage that there is no need to take the introduced CMFD into account (Griffith2018). As all isotopologues are calibrated independently, it is only necessary to use reference gases covering the range of isotopologue abundances as occurring in the samples. This can be realized by using reference gases containing CO2 of similar isotope composition but varying CO2 mole fractions, as described in Griffith (2018) and successfully implemented in Griffith et al. (2012), Flores et al. (2017), Wehr et al. (2013) and Tans et al. (2017). The range of delta values that is measured in samples of atmospheric background air is limited (range in unpolluted troposphere is −9.5 ‰ to −7.5 ‰ and −2 ‰ to +2 ‰ for δ13C and δ18O, respectively; Crotwell et al.2020); hence this also applies to the range of delta values that should be covered by the reference gases when applying the RM. We decided therefore to use the same reference gases to test both calibration methods, varying mainly in the CO2 mole fraction (342.81–424.52 µmol mol−1).

3.3.1 Measurement procedure

The measurement procedure that is used for both calibration methods is based on the alternating measurements of samples/reference gases and the working gas, so the drift-corrected measurement value can be calculated as in Eq. (1). Per sample/reference gas measurement, there are nine iterations of successive sample and working gas measurements, from now on called a measurement series, before switching to the next sample/reference gas measurement series. One measurement series lasts ∼30 min. Sample series are conducted once, while the reference gases series (low and high reference) are repeated four times throughout a measurement sequence. The quality control gas, a gas of known isotope composition which is not included in the calibration procedure, is also measured four times throughout the measurement sequence. One measurement sequence in which 12 samples are measured lasts therefore ∼12 h. For the nine measurement values of each measurement series, outliers are determined using the outlier identification method for very small samples by Rousseeuw and Verboven (2002), and the mean values of the measurement series are calculated.

3.3.2 Ratio method

In the RM, measured isotopologue mole fractions are used for the estimation of isotope ratios (Eq. 1), which are calibrated to the international VPDB-CO2 scale by measurement of several in-house CO2-in-air references within the same measurement sequence. The working gas is used both for drift correction and a first calibration step, and the uncalibrated delta value δS is calculated by

(5) δ S * = r S * r WG * - 1 ,

where S and WG stand for sample and working gas respectively. The calibrated δ13C, δ18O and δ17O based on the working gas that is used are then derived by

(6) δ SCal = 1 + δ WG δ S + δ WG ,

in which δWG is the known delta value of the working gas on the VPDB-CO2 scale.

Up to this point, the procedures are more or less identical to those for IRMS measurements (but without the ion correction and N2O correction unnecessary here). CMFD correction is specific for laser absorption spectroscopy and is crucial (as can be concluded from Sect. 3.1.3) to derive accurate measurement results when calibration is done using the isotope ratios. We developed a calibration method based on the idea that including the measurement of two reference gases covering the CO2 mole fraction range of the measured samples (in our case the low- and high-reference gas) enables the correction of the measured isotope ratios. These two reference gases are measured several times throughout the measurement sequence and a quadratic fit of the mean of the residuals (measured δSCal – assigned δVPDB), including the residual of zero for the (hypothetical) working gas measurement, as a function of the CO2 mole fraction is done, so the following calibration formula can then be determined:

(7) δ VPDB = δ SCal - CO 2 2 a + CO 2 b + c ,

in which a and b are the second- and first-order coefficients respectively, c is the intercept of the quadratic fit of the residuals and the CO2 mole fractions of the two reference gases, [CO2] is the measured CO2 mole fraction and δVPDB is the calibrated δ value on the VPDB scale.

3.3.3 Isotopologue method

The IM as described by Flores et al. (2017) following methods earlier described by Griffith et al. (2012) will be briefly explained here for clarity, before explaining the application of the IM on the SICAS. Basically, the method treats the CO2 isotopologues as if they were independent species, calibrates their mixing ratios individually and only then combines the results to build isotope ratios and delta values. The mole fractions (X) of the four most abundant isotopologues of a measured CO2 sample are determined using a suite (in our case the working gas and the high- and low-reference gas) of references gases with known CO2 mole fractions and isotope compositions. The CO2 mole fractions are ideally chosen such that normally occurring CO2 mole fractions in atmospheric air are bracketed by the two reference gases. The low- and high-reference gases cover the range between 324.81 and 424.52 ppm, meaning that this method is only valid for samples within that range of CO2 concentrations. The actual (or assigned) mole fractions (Xa) of the four most abundant isotopologues of the reference gases can be calculated using calculations 1–11 in Flores et al. (2017), which are listed in Appendix B. Although the nonlinearity of isotopologues as a function of the absolute CO2 mole fraction has not been investigated in this study, it is very likely that nonlinearities occur, according to the results discussed in Sect. 3.1.2. The broad range of CO2 mole fractions that are covered by the reference gases, together with a hypothetical measurement of the working gas (of which the normalized isotopologue abundance will always be 1), enable us to do a quadratic fit of the measured isotopologue abundance as a function of the assigned isotopologue mole fractions, by

(8) X a = X m 2 c + X m d + e ,

in which c and d are the second- and first-order coefficients, respectively, and e is the intercept of the quadratic fit of Xm as a function of Xa of the reference gases. The resulting Xa values are used to calculate the isotope composition using calculation 1–11 in Appendix B. The introduced CMFD due to calibration on measured isotope ratios will not occur with this method, and a CMFD correction is therefore not necessary to yield accurate results.

A complete overview of all calculation steps of both the RM and IM can be found in Appendix C.

Figure 7Quality control gas δ13C (a, b) and δ18O (c, d) measurements for both the RM (a, c) and IM (b, d). The assigned value of the quality control gas is indicated by the dotted black line, and the WMO compatibility goals are indicated by the dotted gray lines. The error bars show the standard error of the measurements. Color of the points indicates whether the measurements were performed in a high-quality (green), medium-quality (black) or low-quality (red) measurement period (see Sect. 4.1 for definitions).


4 Results and discussion

4.1 Monitoring measurement quality and comparison of calibration methods of δ13C and δ18O

To capture the very small signals in time series of the isotope composition of atm-CO2, it is crucial to keep track of the instrument's performance over the course of longer measurement periods. Variations in precision and accuracy of the isotope measurements on the SICAS are monitored by measurement of a quality control gas in every measurement sequence. Since the quality control gas measurement is not used for any correction or calibration procedures, it can be considered a known sample measurement that gives an indication of the overall instrument performance. Based on the WMO compatibility goals required for isotope measurements of atm-CO2 we categorized (high quality – H, medium quality – M – and low quality – L) three measurement periods for both the RM and IM. A period is rated as H if both the mean accuracy and the mean precision (expressed as the standard error) of the quality control gas measurements over that period are within the WMO compatibility goals (0.01 ‰ for δ13C and 0.05 ‰ for δ18O; Crotwell et al.2020); if the accuracy or precision is within the requirements but the other one is not, it is rated M, and if both accuracy and precision do not fulfil the requirements, it is rated L. Measurements of the quality control gas done over the period of 20 November 2019 until 4 February 2020 are shown in Fig. 7, and we assigned three distinct measurement periods based on the quality of the measurements. The mean residuals and standard errors of all quality control gas measurements during the three periods are shown in Table 6.

Table 6Mean residuals and standard errors of the quality control measurements in the three different measurement periods.

Download Print Version | Download XLSX

From the results we learn that the difference in performance between the two methods is minimal. The precision of the quality control gas measurements shows the same results, while the accuracy shows small differences between the methods for the different periods. High-quality performances are reached in period 1 for the δ13C measurements, but in periods 2 and 3 both the precision and the accuracy are worse than 0.01 ‰; hence the measurement quality is low. The δ18O measurements show high-quality performance over the whole period.

4.2 Uncertainty budget

A combined uncertainty consisting of measurement uncertainties and scale uncertainties is calculated for the sample measurements. Measurement uncertainties include the standard error of the sample measurement, the repeatability of all (usually four) measurements of the quality control gas throughout the measurement sequence and the residual of the mean of the quality control gas measurements from the assigned value. The measurement uncertainties will therefore vary with each measurement/measurement sequence. We observe a high repeatability in all sequences included in the analysis of Fig. 7 (eight in total), with standard errors ranging between 0.008 ‰ and 0.03 ‰ and a mean of 0.02 ‰ for δ13C and standard errors ranging between 0.007 ‰ and 0.01 ‰ and a mean of 0.008 ‰ for δ18O, for both methods. The residuals in these sequences show a higher contribution to the combined uncertainty and a small difference between the two calibration methods. The absolute residuals of the RM range between 0.002 ‰ and 0.04 ‰ with a mean of 0.024 ‰ for δ13C and between 0.016 ‰ and 0.05 ‰ with a mean of 0.04 ‰ for δ18O. For the IM the residuals range between 0.006 ‰ and 0.03 ‰ with a mean of 0.02 ‰ for δ13C and between 0.012 ‰ and 0.04 ‰ with a mean of 0.03 ‰ for δ18O. Hence, the RM shows slightly higher contributions to the combined uncertainty as a result of the accuracy of the quality control gas measurements.

The scale uncertainties, which are fixed for all measurement sequences in which the working gas, low reference and high reference are used for the sample calibration, were simulated using the Monte Carlo method. Input values were generated by choosing random numbers of normal distribution with the assigned value and uncertainty as in Table 5, being the mean and the standard deviation around the mean, respectively. As the RM and IM follow different calibration schemes, the Monte Carlo simulations are discussed separately; for the RM the scale uncertainties of the assigned delta values result in an uncertainty in the calculated residuals which are quadratically fitted against the measured CO2 mole fraction. The average uncertainties in the calibrated delta values of the five simulations are 0.03 ‰ and 0.05 ‰ for δ13C and δ18O, respectively.

Besides the uncertainties introduced by the scale uncertainties of the delta values, the calibrated measurements of the IM are also affected by the scale uncertainties of the CO2 mole fractions. Both the uncertainties in the delta values and in the CO2 mole fractions affect the calculated assigned isotopologue abundances, which are quadratically fitted against the measured isotopologue abundances. The uncertainties in the assigned delta values result in average uncertainties of 0.03 ‰ and 0.06 ‰ for δ13C and δ18O, respectively. The uncertainties in the assigned CO2 mole fractions result in uncertainties of 0.005 ‰ and 0.018 ‰ for δ13C and δ18O, respectively, and are small compared to the uncertainties of the assigned delta values.

Reducing the combined uncertainty of the δ13C and δ18O measurements of the SICAS will be most effective by determining the isotope composition of the reference gases with a lower uncertainty on the VPDB-CO2 scale.

4.3 Intercomparison flask measurements

To test the accuracy of SICAS flask measurements over a wide range of CO2 mixing ratios, as well as testing the lab compatibility of the SICAS measurements, we measured flask samples that are part of an ongoing lab intercomparison of atmospheric trace gas measurements including the δ13C and δ18O of CO2 (Levin et al.2004). Since 2002, the sausage flask intercomparison program (from now on defined as ICP) has provided aliquots of three high-pressure cylinders containing natural air covering a CO2 mixing ratio range of 340–450 µmol mol−1 every 2 to 3 months (occasionally longer periods). Participating laboratories send six flasks to the ICOS-CAL lab in Jena where these are filled with air from the three cylinders (two flasks per cylinder) with the sausage method. The ICP provides therefore the opportunity to compare flask measurements on the SICAS with IRMS flask measurements of the MPI-BGC and other groups. We measured sausage series 90–94, which were filled between April 2018 and January 2020, and calibrated the isotope measurements both with the RM and the IM. SICAS measurements took place in the period from December 2019 to April 2020, with the consequence that the storage time of the flasks varies between 3 and 20 months. To place these results in the context of intercomparison results of well-established isotope and measurement laboratories, the ICP results of the Earth System Research Laboratory of the National Oceanic and Atmospheric Administration (NOAA) (Trolier et al.1996) were also compared to the MPI-BGC results for the same sausage series.

The lab intercomparison is presented in the usual way: the mean and standard deviation of the differences between our SICAS δ13C and δ18O results (both RM and IM calibrated) and the MPI-BGC ones are shown in Table 7, along with the NOAA-MPI-BGC differences. The mean values of the differences for the SICAS RM and NOAA results are both below 0.01 ‰, while the standard deviations of the differences are 0.05 ‰ and 0.07 ‰, respectively. The SICAS results calibrated with the IM show an offset with MPI-BGC of −0.013 ‰ and a standard deviation of the differences of 0.07 ‰. We can therefore conclude that the differences in performance between the RM and the IM are minimal, and both methods show comparable results for the measured differences between MPI-BGC as for the differences between the NOAA and the MPI-BGC.

Table 7Lab intercomparison of ICP sausage 90–94 results, only including data points within the CO2 mole fraction range of the used calibration tanks (342.81–424.52 µmol mol−1). Differences between the SICAS, of both calibration methods, as well as the NOAA IRMS results and the MPI-BGC IRMS results are shown. The mean difference as well as the standard deviation of the differences of the δ13C and δ18O are shown.

Download Print Version | Download XLSX

Table 8Results of isotope measurements of quality control gas from the tank and quality control gas air in flasks (calibrated with the RM) at different periods after the flask filling procedure. The last column shows the number of cylinder measurements or the number of flasks that were used to calculate the average and the standard deviation.

Download Print Version | Download XLSX

When we compare the δ18O measurements, we find that the SICAS results are consequently significantly more depleted with an average difference of −0.4 ‰ compared to the MPI-BGC results and that the differences vary strongly with a standard deviation of 0.16 ‰. δ18O results of the ICP program show in general a larger scatter among the labs than δ13C results (Levin et al.2004), as is also visible in Table 7 for the NOAA-MPI-BGC differences. The differences between the SICAS- and the MPI-BGC results, however, are far larger than those (or than in fact all differences in the ICP program). The reason for this signal being too depleted is presumably equilibration of CO2 with water molecules on the glass surface inside the CIO-type sample flasks during storage. Earlier (unpublished) results from our CO2 extraction system indicated that the water content of our dried atmospheric air samples increased as a function of time inside the flasks. Our atmospheric samples are stored at atmospheric pressure or lower (down to ∼800 mbar) when part of the sample has been consumed by different measurement devices. The CIO flasks are sealed with two Louwers–Hapert valves and Viton O-rings of which it is known that permeation of water vapor (as well as other gases) occurs over time (Sturm et al.2004). Both the pressure gradient and the water vapor gradient between the lab atmosphere and the dry sample air inside the flask lead to permeation of water molecules through the valve seals. To check this hypothesis, an experiment was conducted in which CIO flasks were filled with quality control gas and were measured the same day of the filling procedure and 1 week and 3 months later (see Table 8). The results show no significant change in the δ13C, while for the δ18O there is a strong depletion of the flask measurements after 3 months, deviating more than −0.2 ‰ in comparison to the cylinder measurements. After 1 week there is no change in the δ18O, indicating that depletion of the δ18O in the CIO flasks occurs over longer time periods. As the flasks from the ICP were measured at the SICAS after relatively long storage times, sometimes almost 2 years, this is likely the explanation of the values being too depleted in comparison to the MPI-BGC results. A depletion twice as small as for δ18O is observed in the δ17O values, as one would expect for isotopic exchange with water. Further investigations about the changing oxygen isotope signal in CIO-sample flasks are being conducted with the aim to be able to make reliable assessments of the quality of δ18O and δ17O flask measurements on the SICAS.

To check the performance of the SICAS for both the IM and RM over the wide CO2 range that is covered by the ICP sausage samples, the differences between the MPI-BGC and the SICAS results are plotted in Fig. 8 against the measured CO2 mole fraction. Shown is that for both methods the highest differences are seen at the higher end of the CO2 mole fraction above 425 ppm and therefore far out of the range that is covered by the high and low references (∼343–425 ppm). Extrapolation of the calibration methods outside the CO2 mole fraction range of the reference gases yields worse compatibility with MPI-BGC, possibly due to the nonlinear character of both the isotopologue CO2 dependency and the ratio CO2 dependency. It should therefore be concluded that, to achieve highly accurate results of isotope measurements over the whole range of CO2 mole fractions found in atmospheric samples, the range covered by the reference gases would ideally be changed to ∼380–450 ppm. The results of the IM are slightly better in the CO2 range above 425 ppm; specifically, the point closest to 440 ppm shows a significantly smaller residual (∼0.1 ‰ less) than the RM. The better result of extrapolation of the determined calibration curves for the IM method could be due to the lesser degree of nonlinearity of the measured isotopologue abundances as a function of the assigned isotopologue abundances, in comparison to the nonlinearity of the measured isotope ratios as a function of the CO2 mole fraction. More points in this higher range are needed, however, to draw any further conclusions on this matter.

Figure 8Results of the intercomparison of δ13C measurements on the SICAS and on the IRMS facility at the MPI-BGC for both the RM (a) and IM (b). The MPI-BGC results were subtracted from the SICAS results, and the error bars show the combined uncertainty of the SICAS measurements. The dotted gray lines show the 0.03 ‰ range of residuals. Red data points are outside of the CO2 mole fraction range of the reference gases.


4.4 Potential of SICAS Δ17O measurements for atmospheric research

With the direct measurement of δ17O in addition to δ18O (triple oxygen isotope composition) of atm-CO2, the δ17O excess (Δ17O) can be calculated. Δ17O measurements can be a tracer for biosphere activity (Hoag et al.2005), atmospheric circulation patterns (Mrozek et al.2016) and different combustion processes (Horváth et al.2012). The Δ17O is usually defined as

(9) Δ 17 O = ln 1 + δ 17 O - λ ln 1 + δ 18 O .

Variations in the Δ17O signal in the troposphere mainly depend on biosphere activity and the influx of stratospheric CO2 (Koren et al.2019; Hofmann et al.2017; Hoag et al.2005). High measurement precision and accuracy of both the δ18O and the δ17O are needed to capture spatial gradients and seasonal cycles in the Δ17O, of which seasonal variations of 0.13 ‰ (Hofmann et al.2017) and 0.211 ‰ (Liang et al.2017) have been reported. So far it has been an extremely complex and time intensive process to measure δ17O of CO2 using DI-IRMS (Hofmann and Pack2010; Barkan and Luz2012; Mahata et al.2013; Adnew et al.2019). Dual-laser absorption spectroscopy as presented in this paper does not require any sample preparation and would therefore be a great step forward in the use of Δ17O as a tracer for atm-CO2. Here we present the measurement precision and stability of the δ17O, as well as the Δ17O measurements of our quality control tank in Fig. 9 and Table 9, and we evaluate the potential for contribution in the field of triple oxygen isotope composition studies.

Figure 9Quality control gas δ17O (a, b) and δ17O (c, d) measurement averages for the three measurement periods for both the RM (a, c) and the IM (b, d). The averages are indicated by the dotted black line, and the dotted gray line shows the 0.05 ‰ range around the average. The error bars show the standard error of the measurements.


Table 9Average of the residuals from the assigned value and mean of the standard error of the quality control gas δ17O and δ17O measurements per period for both the RM and the IM.

Download Print Version | Download XLSX

All results show too enriched values according to the assigned values, which is probably due to the fact that the assigned δ17O values of the low and high references have been determined indirectly, as discussed in Sect. 3.2. A direct determination of the δ17O values of our low and high references would supposedly improve the accuracy of both methods. The Δ17O accuracy is dependent on both the δ17O and δ18O results, where Δ17O values will deviate more if those results deviate in opposite directions and vice versa. Furthermore, it is striking that the mean standard errors of measurement periods 2 and 3 are twice as low for the IM than for the RM. The r627, used for the RM, is calculated by dividing X627, derived from laser 1, by X626 derived from laser 2. It can be that the two lasers do not drift in the same direction, and the advantage of canceling out these drifts by dividing the two measured values will not apply. The outlier analysis of the IM might in that case be more effective as it is performed on both the measured 16O and 17O abundances, while for the RM it is only performed on the r627. A comparison of the correlation coefficients between the 627 peak results and the 626 peak results from both lasers shows no significant difference (and a value of ∼0.65), meaning that using the 626 peak of laser 1 for the δ17O calibration will not improve the precision of the RM results.

Due to the lower seasonal variations of the Δ17O values, even higher measurement precision is a prerequisite, and in Hofmann et al. (2017) it is stated that a measurement precision of 0.01 ‰ or better is required to capture these variations and to use the Δ17O value as a potential tracer for gross primary productivity (GPP). This precision has now not yet been achieved, but the results of the IM calibrated values show that small improvements in the measurement precision of the SICAS can bring the Δ17O measurements close to the 0.01 ‰ precision. This could for instance be accomplished by deciding to conduct more iterations per measurement series, if sample size allows for this. In Sect. 2.2 the contamination of the mirrors was discussed as the potential cause for the decreased signal-to-noise ratio over the period September 2017–July 2019. Placing new mirrors in the optical cell might therefore improve the quality of the measurements further. As the quality of the Δ17O measurements depends directly on the quality of the δ18O and the δ17O measurements, it will be important to monitor the measurement quality of both isotope values over time using the measurements of the quality control gas. If SICAS measurements are to be used for comparison with Δ17O measurements from other labs or measurement devices, it is necessary to add the error introduced by the scale uncertainties of the reference gases as well. For both the δ17O and Δ17O these uncertainties are 0.08 ‰, as calculated with a Monte Carlo simulation, described in Sect. 4.2. As long as only measurements from this device are used, seasonal and diurnal cycles are measured with much lower uncertainties. The high residuals found for the quality control gas measurements of the δ17O and Δ17O show that these uncertainties are probably an underestimation, as the assigned values of the low and high reference, which were not directly measured at the IMAU, are not known with high accuracy. For reducing the combined uncertainty it is therefore crucial to have all reference gases directly determined for their δ17O values, as well as reducing the scale uncertainties of both the δ17O and δ18O values of the reference tanks.

5 Conclusions and outlook

In this study we discuss the measurement performances of our Aerodyne dual-laser absorption spectrometer in static mode for stable isotope measurements of atm-CO2 in dry whole-air samples. We implemented two different calibration methods based on the same measurement procedure, the ratio method (RM) and the isotopologue method (IM). Short-term instrumental drift can effectively be corrected by continuously alternating sample measurements with measurements of a machine working gas. Nine aliquots are measured per sample/reference gas, and two reference gases covering a wide range of CO2 mole fractions, as well as a quality control tank serving as a known unknown, are measured four times throughout a measurement sequence. The RM is based on calibration of measured isotopologue ratios (or delta values, calculated directly from the measured sample and working gas ratios), including correction for a nonlinear CO2 mole fraction dependency. This correction is determined by doing a quadratic fit of the residuals of the calibrated delta values of the reference gases as a function of the measured CO2 mole fractions. The IM is based on calibration of measured isotopologue abundances, using a quadratic fit of the measured values of the reference gases as a function of the assigned isotopologue values. In optimal measurement conditions, precision and accuracy of <0.01 ‰ and <0.05 ‰ for δ13C and δ18O are reached for measurements of the quality control tank for both calibration methods. The combined uncertainty of the measurements also includes the repeatability of the four quality control gas measurements throughout the measurement sequence, with mean values of 0.014 ‰ and 0.012 ‰. The last components in the combined uncertainty calculation are caused by scale uncertainties of the reference gases used for the sample calibration, which are 0.03 ‰ and 0.05 ‰ for δ13C and δ18O of the RM, respectively, and 0.03 ‰ and 0.06 ‰ for δ13C and δ18O of the IM, respectively.

A comparison of SICAS results, for both calibration methods, with results from the MPI-BGC from the sausage ICP shows that sample results within the CO2 mole fraction range of both methods are of similar quality when compared to the MPI-BGC results. Better results were achieved for the IM for samples outside the CO2 mole fraction range, but more measurements are needed to determine whether the IM is indeed less vulnerable to extrapolation of the calibration. As extrapolation should at any time be avoided, using reference gases that cover the range of atmospheric samples is advisable for more reliable measurement results. We found that δ18O measurements were consequently too depleted due to storage times of the CIO flasks being too long before measurement. Future investigations will give more insight into the stability of the oxygen isotopes within the CIO flasks, and we will evaluate the possibility of a correction based on storage time.

δ17O and Δ17O results of the quality control gas show values that are consequently too enriched, which is probably caused due to the indirect determination of the δ17O values of two of the reference gases. The measurement precision is significantly better for the IM, with standard errors not higher than 0.02 ‰, while the measurement precision of the RM shows standard errors ranging between 0.02 ‰ and 0.05 ‰. Results of the IM come close to the required 0.01 ‰ precision to capture seasonal variations of the atmospheric Δ17O signal. For a combined uncertainty of the δ17O and Δ17O values, an additional uncertainty of 0.08 ‰ must be added due to effects of the scale uncertainties of the reference gases, indicating that improved determination of the oxygen stable isotope values of the reference gases will be essential to reach high-precision Δ17O measurements that are compatible with measurements from other labs. We will show the actual achievements of Δ17O measurements with this instrument for a record of atmospheric CO2 samples of our atmospheric monitoring station Lutjewad in a forthcoming paper.

Appendix A: Preparation procedure for CO2-in-air samples

The pure CO2 aliquots were prepared by connecting a 20 mL flask containing a pure CO2 local reference gas to a calibrated adjustable volume. The required amount of CO2 in the adjustable volume could be determined by measuring the pressure at a resolution of 1 mbar using a pressure sensor (Keller LEO 2). Both the sample flask and adjustable volume were connected to a vacuum (3.3×10-5 mbar) glass line. The CO2 in the adjustable volume was transferred cryogenically (using liquid nitrogen) into a small glass tube shape attachment on the side of the evacuated sample flask which was custom-made for this purpose, and subsequently the zero-air dilutor gas was added. The dilutor gas consists of natural air scrubbed of CO2 and H2O using Ascarite® (sodium hydroxide coated silica, Sigma-Aldrich) and Sicapent® (phosphoric anhydride, phosphorus(V) oxide), which results in dry, CO2-free natural air. For experiment 2, additional samples were prepared using synthetic air mixtures with and without 1 % Argon as dilutor gas for evaluation of the effect of air composition on the CMFD (see also Sect. 3.1.6). After closing the flask, the mixture was put to rest for at least one night before measurement to ensure the CO2 and the dilutor were completely mixed.

Appendix B: Equations for calculation of isotopologue mole fractions

Individual isotopologues of standards of known CO2 mole fractions and isotope composition are calculated for the IM calibration method by the equations below, according to Flores et al. (2017), starting with equations for the atomic abundances X in each of the calibration gas mixtures (Eqs. B1–B5):




(B8) R 18 = R VPDB - CO 2 18 1 + δ 18 O

and δ13C and δ18O are the delta values. RVPDB-CO213 (0.011180), RVPDB-CO217 (0.0003931) and RVPDB-CO218 (0.00208835) values were taken from Brand et al. (2010) for VPDB-CO2. Then each carbon dioxide isotopologue mole fraction in the reference gas was calculated according to its composition using Eqs. (B9) and (B10):


For a more elaborated explanation of these equations, see Flores et al. (2017).

Appendix C: Step-by-step calculation steps

C1 Ratio method

  1. Calculate ratios from the measured isotopologue abundances:

    (C1) r * = X * X 626 ,

    with r being the ratio, X the measured isotopologue abundance as the default output, scaled in this way for the natural abundance, * stands for one of the three rare isotopologues (636, 628 or 627) and 626 standard for the abundant isotopologue. The CO2 mole fraction is calculated by

    (C2) CO 2 = A 626 + A 636 + A 628 + A 627 ,

    with A being the actual measured abundance, calculated back in this way using the natural abundance values for the isotopologues as defined in Rothman et al. (2013).

  2. Use only the relevant interval (in our case 30–60 s) from measured ratio and [CO2] per measurement.

  3. Do a drift correction and calculate the uncalibrated delta value by

    (C3) δ S * ( t ) dc = r S ( t ) r WG ( t ) - 1 ,

    with S for sample, t for time of the measurement, dc for drift-corrected and WG for working gas and with rWG(t)* derived from applying a time-dependent linear fit of the rWG(t-1)* and rWG(t+1)*, following

    (C4) r WG ( t ) * = α + β t .

    The [CO2] is also drift-corrected by

    (C5) CO 2 S ( t ) dc = CO 2 S ( t ) CO 2 WG ( t ) .
  4. Group all δ* values and [CO2] values per measurement series, and do an outlier analysis per series. We adapted the method as described in Rousseeuw and Verboven (2002).

    • a.

      Define the variable “sborder” (sborder = 2), defining the strictness of filtering.

    • b.

      Calculate for all values in the series the (absolute) deviation from the median of the series, resulting in a new series containing the distance from the median (DM).

    • c.

      Calculate the MAD (median absolute deviation), by

      (C6) MAD = 1.483 median ( DM ) .
    • d.

      Calculate per value of the series the deviation with the following equation:

      (C7) deviation = abs x 1 , 2 N - median sborder MAD ,

      with x1,2 … N the measurement values from the measurement series.

    • e.

      If the deviation of a value is higher than 1, the value is identified as an outlier.

  5. Calculate the mean and standard error per measurement series, excluding the identified outliers.

  6. Do first a one-point-calibration on all mean values using the known values of the working gas, by

    (C8) δ * A = ( 1 + WG * A ) δ * + WG * A ,

    with A standing for atom (C or O) and WG* being the assigned isotope value of the working gas. The [CO2] is calibrated by

    (C9) CO 2 c = CO 2 S CO 2 WG ,

    with [CO2]c the calibrated [CO2] value, [CO2]S the mean of the sample [CO2] measurement and [CO2]WG the assigned CO2 mole fraction value of the working gas.

  7. Calculate the means of the δ*A values and the [CO2]c values of the high- and the low-reference gas measurements that were done throughout the measurement sequence (we normally do four measurements of both reference cylinder) and calculate the residual of the means from their assigned δ*A values.

  8. Use the two calculated residuals, together with a residual of 0 for a hypothetical working gas measurement, to do a quadratic fit (ax2+bx+c) of the residuals as a function of the [CO2]c, and calculate the final δ*A on the VPDB scale by

    (C10) δ * A VPDB = δ * A - CO 2 c 2 a + CO 2 c b + c .
  9. Calculate the combined uncertainty by

    (C11) cu δ * A = sud 2 + QC SE 2 + QC res 2 + SE m 2 ,

    in which sud is the scale uncertainty of delta values (derived from a Monte Carlo simulation), QCSE is the standard error of all (usually four) quality control gas measurements throughout the measurement sequence, QCres is the residual of the mean of all quality control gas measurements and SEm is the standard error of the measurement.

C2 Isotopologue method

  1. Use only the relevant interval (30–60 s) from measured isotopologue abundances per measurement.

  2. Do a drift correction by

    (C12) a S ( t ) dc * = a S ( t ) a WG ( t ) .

    The aWG(t) is derived with the same method as for the RM.

  3. Group all a* per measurement series, and do an outlier analysis per sample. The same method as for the RM is used.

  4. Calculate the mean per measurement series, excluding the identified outliers.

  5. Calculate the quadratic calibration curves (αx2+βx+γ) for all four isotopologues, by fitting the mean of all low- and high-reference measurements (usually four per measurement sequence) and an additional value of 1 for the hypothetical working gas measurement as a function of the assigned isotopologue mole fraction.

  6. Calculate the calibrated isotopologue mole fractions of all four isotopologues for all measurements, so not for the mean of the grouped measurements but for all drift-corrected a* values from step 2:

    (C13) a S ( cal ) * = α a * 2 + β a * + γ .
  7. Calculate the calibrated delta values using the calibrated isotopologue abundances for all sample measurements.

  8. Group all δ*A per measurement series, and do an outlier analysis per sample, again using the same method as described in the RM.

  9. Calculate the mean and standard error of all δ*A per measurement series, excluding the identified outliers.

  10. Calculate the combined uncertainty of the measurement:

    (C14) cu δ * A = sud 2 + suc 2 + QC SE 2 + QC res 2 + SE m 2 ,

    in which suc is uncertainty introduced by the scale uncertainty of the CO2 mole fractions.

Appendix D: Setup for preparation of low reference

The setup is as follows: a high reference, filled up to ∼150 bar with dry natural air, is connected to a similar, empty cylinder. Half of the air in the high-reference tank will be transferred (passive transfer using the pressure difference) into the empty cylinder to produce the low reference. The CO2 mole fraction in the low reference is reduced by leading part of the air over an Ascarite®-filled cartridge that removes all CO2 from the air, so no isotope fractionation will occur. Successively it is led over a magnesium-perchlorate-filled cartridge to remove water from the air that is potentially stored in the hydrophilic Ascarite®. A needle valve installed before the cartridges creates a low flow to ensure the complete removal of the CO2 from the air. The pressure sensor installed after the repercussion valve enables us to estimate when the low-reference cylinder is filled with the amount of CO2-free air needed to obtain the preferred CO2 mole fraction. When the preferred amount of CO2-free air is transferred into the low-reference cylinder, the cartridges are decoupled from the system to transfer the rest of the air from the high-reference cylinder to the low-reference cylinder.

Figure D1Setup for the preparation of the low reference. The cylinder on the left is the high reference, filled until ∼150 bar at start time of preparation. The right cylinder is empty at start time of preparation.


Appendix E: Analysis cross-contamination

To determine whether cross-contamination has the potential to affect isotope measurements on the SICAS, a simulation was conducted in which we use the measurement procedure described in this paper. Input in the simulation is an experimentally derived value, which expresses how much a measured sample is affected by the sample that was measured in the optical cell before. The experiment was conducted as follows: the high reference was measured eight times in a row, each time letting in a new aliquot, followed by the low reference which was also measured eight times in a row, and this procedure was repeated three times. The usual flushing procedure was applied every time there was a switch between the cylinders. It can be expected that the first measurement of a series of eight of the low reference is affected the most by the preceding measurement of the high-reference gas. The last measurements of a series of eight will be affected less and will be closer to the “true” value. We quantified this effect by applying the following equation to all series of measured isotopologues:

(E1) CC = M 1 * - mean M 6 , 7 , 8 * mean M 6 , 7 , 8 * 100 ,

in which CC is the cross-contamination in percent, and M stands for measurement, with * being the isotopologue and the number indicated at the underscore the number of the measurement. The CC values we observed were low, ranging from undetectable up to 0.01 % at most. We used this highest value for our simulation.

A simulation for a measurement sequence was set up in Excel, following the measurement procedure as described in this paper, only using three sample measurements per measurement series instead of nine. Included in the simulation are measurements of the low- and high-reference gas and two hypothetical samples with CO2 concentrations of 480 and 340 ppm, a δ13C values of −7 ‰ and −11 ‰ and δ18O values of −1 ‰ and −4 ‰, respectively. All measurements are alternated with measurements of the working gas, according to the measurement procedure described in this paper. We use the actual values for CO2 concentration and isotope composition of all reference gases in the simulation. The measurements were simulated by

(E2) M t * = M t - 1 * 0.01 × 10 - 2 + 1 - 0.01 × 10 - 2 Tr * ,

with Mt* being the simulated measurement at time t and * indicating which isotopologue measurement is simulated, Mt−1 being the preceding simulated measurement and Tr being the true isotopologue abundance of the sample or reference gas that is being measured at time t. The first value that is put in the simulation contains the true values for all measured isotopologue abundances. For all sample measurements, a normalized measurement is calculated by dividing Mt* by Mt-1* (being the working gas measurement).

These simulated, normalized measurements of the low- and high-reference gases are used to do a linear fit as a function of the true value and in doing so calculate the calibration curves. These curves are used to calculate the calibrated sample measurements, and the measured δ13C and δ18O measurements can be calculated. We find deviations from the measured simulation values of maximum 0.0002 ‰ for both δ13C and δ18O, so 2 orders of magnitude lower than the measurement precision.

Data availability

All data that have been used for this study which were measured at the SICAS can be found in the Supplement.


The supplement related to this article is available online at:

Author contributions

PMS, HAS and HAJM conceived the experiments, which were conducted by PMS and HAS, and PMS carried out the data analysis. DDN and JBM optimized the fit and contributed with technical assistance for development of the gas handling system as well as solving problems with the instrumentation. The manuscript was written by PMS, and HAS and HAJM contributed with discussions and comments throughout the writing process.

Competing interests

Authors Dave D. Nelson and J. Barry McManus work for Aerodyne Research Inc., which is the company that developed the instrument described in this study.


We would like to thank Heiko M. Moossen, as well as his colleagues at the MPI-BGC, for measuring our reference cylinders and providing us with the data that we required for the intercomparison. We greatly acknowledge the help of Getachew Adnew from IMAU who prepared the samples for the first CMFD experiment and helped measuring them. He also measured the δ17O composition of our working gas and quality control gas. We also want to thank Marcel de Vries, Bert A. M. Kers and Ramon Ritchie for helping us to develop the sampling system as well as the required software development and Dipayan Paul for general assistance. We thank the anonymous reviewer, David Griffith and Edgar Flores for their valuable comments, as well as the associate editor Tim Arnold for his contributions.

Financial support

This research has been supported by the European Metrology Programme for Innovation and Research (grant no. 16ENV06SIRS).

Review statement

This paper was edited by Tim Arnold and reviewed by Edgar Flores, David Griffith, and one anonymous referee.


Adnew, G. A., Hofmann, M. E., Paul, D., Laskar, A., Surma, J., Albrecht, N., Pack, A., Schwieters, J., Koren, G., Peters, W., and Röckmann, T.: Determination of the triple oxygen and carbon isotopic composition of CO2 from atomic ion fragments formed in the ion source of the 253 Ultra high‐resolution isotope ratio mass spectrometer, Rapid Commun. Mass Spectrom., 33, 1363–1380,, 2019. a, b

Allison, C. and Francey, R.: High precision stable isotope measurements of atmospheric trace gases, Reference and intercomparison materials for stable isotopes of light elements, IAEA-TECDOC, Vienna, 131–154, 1995. a

Allison, C., Francey, R., and Meijer, H.: Recommendations for the reporting of stable isotope measurements of carbon and oxygen in CO2 gas, IAEA-TECDOC, Vienna, 155–162, 1995. a

Assonov, S. S. and Brenninkmeijer, C. A. M.: On the 17O correction for CO2 mass spectrometric isotopic analysis, Rapid Commun. Mass Spectrom., 17, 1007–1016,, 2003. a

Barkan, E. and Luz, B.: High-precision measurements of 17O/16O and 18O/16O ratios in CO2, Rapid Commun. Mass Spectrom., 26, 2733–2738,, 2012. a

Becker, J. F., Sauke, T. B., and Loewenstein, M.: Stable isotope analysis using tunable diode laser spectroscopy, Appl. Optics, 31, 1921–1927, 1992. a

Braden-Behrens, J., Yan, Y., and Knohl, A.: A new instrument for stable isotope measurements of 13C and 18O in CO2 – Instrument performance and ecological application of the Delta Ray IRIS analyzer, Atmos. Meas. Tech., 10, 4537–4560,, 2017. a

Brand, W. A., Assonov, S. S., and Coplen, T. B.: Correction for the 17O interference in δ13C measurements when analyzing CO2 with stable isotope mass spectrometry (IUPAC Technical Report), Pure Appl. Chem., 82, 1719–1733,, 2010. a, b

Crotwell, A., Lee, H., and Steinbacher, M.: 20th WMO/IAEA Meeting on Carbon Dioxide, Other Greenhouse Gases and Related Tracers Measurement Techniques (GGMT-2019), GAW Report 255, 2–5 September 2019, Jeju Island, South Korea, 2020. a, b, c

Erdélyi, M., Richter, D., and Tittel, F.: 13CO2/12CO2 isotopic ratio measurements using a difference frequency-based sensor operating at 4.35 µm, Appl. Phys. B, 75, 289–295,, 2002. a

Flores, E., Viallon, J., Moussay, P., Griffith, D. W. T., and Wielgosz, R. I.: Calibration Strategies for FT-IR and Other Isotope Ratio Infrared Spectrometer Instruments for Accurate δ13C and δ18O Measurements of CO2 in Air, Anal. Chem., 89, 3648–3655,, 2017. a, b, c, d, e

Gagliardi, G., Castrillo, A., Iannone, R. Q., Kerstel, E. R., and Gianfrani, L.: High-precision determination of the 13CO2/12CO2 isotope ratio using a portable 2.008-µm diode-laser spectrometer, Appl. Phys. B, 77, 119–124,, 2003. a

Gordon, I. E., Rothman, L. S., Hill, C., Kochanov, R. V., Tan, Y., Bernath, P. F., Birk, M., Boudon, V., Campargue, A., Chance, K. V., Drouin, B. J., Flaud, J. M., Gamache, R. R., Hodges, J. T., Jacquemart, D., Perevalov, V. I., Perrin, A., Shine, K. P., Smith, M. A., Tennyson, J., Toon, G. C., Tran, H., Tyuterev, V. G., Barbe, A., Császár, A. G., Devi, V. M., Furtenbacher, T., Harrison, J. J., Hartmann, J. M., Jolly, A., Johnson, T. J., Karman, T., Kleiner, I., Kyuberis, A. A., Loos, J., Lyulin, O. M., Massie, S. T., Mikhailenko, S. N., Moazzen-Ahmadi, N., Müller, H. S., Naumenko, O. V., Nikitin, A. V., Polyansky, O. L., Rey, M., Rotger, M., Sharpe, S. W., Sung, K., Starikova, E., Tashkun, S. A., Auwera, J. V., Wagner, G., Wilzewski, J., Wcisło, P., Yu, S., and Zak, E. J.: The HITRAN2016 molecular spectroscopic database, J. Quant. Spectrosc. Ra., 203, 3–69,, 2017. a

Griffith, D. W.: Calibration of isotopologue-specific optical trace gas analysers: A practical guide, Atmos. Meas. Tech., 11, 6189–6201,, 2018. a, b, c, d, e, f

Griffith, D. W., Deutscher, N. M., Caldow, C., Kettlewell, G., Riggenbach, M., and Hammer, S.: A Fourier transform infrared trace gas and isotope analyser for atmospheric applications, Atmos. Meas. Tech., 5, 2481–2498,, 2012. a, b, c, d, e, f

Harris, S. J., Liisberg, J., Xia, L., Wei, J., Zeyer, K., Yu, L., Barthel, M., Wolf, B., Kelly, B. F., Cendón, D. I., Blunier, T., Six, J., and Mohn, J.: N2O isotopocule measurements using laser spectroscopy: Analyzer characterization and intercomparison, Atmos. Meas. Tech., 13, 2797–2831,, 2020. a

Hoag, K. J., Still, C. J., Fung, I. Y., and Boering, K. A.: Triple oxygen isotope composition of tropospheric carbon dioxide as a tracer of terrestrial gross carbon fluxes, Geophys. Res. Lett., 32, 1–5,, 2005. a, b

Hofmann, M. and Pack, A.: Development of a technique for high-precision analysis of triple oxygen isotope ratios in carbon dioxide, Anal. Chem., 82, 4357–4361, 2010. a

Hofmann, M. E. G., Horváth, B., Schneider, L., Peters, W., Schützenmeister, K., and Pack, A.: Atmospheric measurements of Δ17O in CO2 in Göttingen, Germany reveal a seasonal cycle driven by biospheric uptake, Geochim. Cosmochim. Ac., 199, 143–163,, 2017. a, b, c

Horváth, B., Hofmann, M. E., and Pack, A.: On the triple oxygen isotope composition of carbon dioxide from some combustion processes, Geochim. Cosmochim. Ac., 95, 160–168,, 2012. a

IAEA: Stable isotope measurement techniques for atmospheric greenhouse gases, IAEA-TECDOC, Vienna, 2002. a

IAEA: Reference Sheet: Certified Reference Material: IAEA-603 (calcite) – Stable Isotope Reference Material, available at: /IAEA-603.htm (last access: 28 May 2021), 2016. a

Kerstel, E. R., Van Trigt, R., Dam, N., Reuss, J., and Meijer, H. A. J.: Simultaneous determination of the 2H/1H, 17O/16O, and 18O/16O isotope abundance ratios in water by means of laser spectrometry, Anal. Chem., 71, 5297–5303,, 1999. a

Koren, G., Schneider, L., van der Velde, I. R., van Schaik, E., Gromov, S. S., Adnew, G. A., Mrozek Martino, D. J., Hofmann, M. E. G., Liang, M.-C., Mahata, S., Bergamaschi, P., van der Laan-Luijkx, I. T., Krol, M. C., Röckmann, T., and Peters, W.: Global 3-D Simulations of the Triple Oxygen Isotope Signature Δ17O in Atmospheric CO2, J. Geophys. Res.-Atmos., 124, 8808–8836,, 2019. a, b

Laskar, A. H., Mahata, S., and Liang, M. C.: Identification of Anthropogenic CO2 Using Triple Oxygen and Clumped Isotopes, Environ. Sci. Technol., 50, 11806–11814,, 2016. a

Leuenberger, M. C., Schibig, M. F., and Nyfeler, P.: Gas adsorption and desorption effects on cylinders and their importance for long-term gas records, Atmos. Meas. Tech., 8, 5289–5299,, 2015. a

Levin, I., Facklam, C., Schmidt, M., Ramonet, M., Ciais, P., Xueref, I., Langenfelds, R., Allison, C., Francey, R., Jordan, A., Rothe, M., Brand, W. A., Neubert, R. E., Meijer, H. A. J., Machida, T., and Mukai, H.: Results of inter-comparison programme for analysis of “sausage” flask air samples, Special report 2, Terrestrial and Atmospheric Carbon Observing System Infrastructure, Gif-Sur-Yvette, France, 2004. a, b

Liang, M. C., Mahata, S., Laskar, A. H., and Bhattacharya, S. K.: Spatiotemporal variability of oxygen isotope anomaly in near surface air CO2 over urban, semi-urban and ocean areas in and around Taiwan, Aerosol Air Qual. Res., 17, 706–720,, 2017. a

Luz, B., Barkan, E., and Bender, M. L.: Triple-isotope composition of atmospheric oxygen as a tracer of biosphere productivity, Lett. Nat., 400, 547–550, 1999. a

Mahata, S., Bhattacharya, S. K., Wang, C. H., and Liang, M. C.: Oxygen isotope exchange between O2 and CO2 over hot platinum: An innovative technique for measuring Δ17O in CO2, Anal. Chem., 85, 6894–6901,, 2013. a

McManus, J. B., Nelson, D. D., Shorter, J. H., Jimenez, R., Herndon, S., Saleska, S., and Zahniser, M.: A high precision pulsed quantum cascade laser for measurements of stable isotopes of carbon dioxide, J. Modern Opt., 52, 2309–2321,, 2005. a

McManus, J. B., Nelson, D. D., and Zahniser, M. S.: Design and performance of a dual-laser instrument for multiple isotopologues of carbon dioxide and water, Opt. Exp., 23, 6569,, 2015. a, b, c, d, e, f, g

Meijer, H. A. J.: Stable isotope quality assurance using the `Calibrated IRMS' strategy, Isotop. Environ. Health Stud., 45, 150–163,, 2009. a

Meijer, H. A. J., Neubert, R. E. M., and Visser, G. H.: Cross contamination in dual inlet isotope ratio mass spectrometers, Int. J. Mass Spectrom., 198, 45–61,, 2000. a

Mrozek, D. J., Van Der Veen, C., Hofmann, M. E., Chen, H., Kivi, R., Heikkinen, P., and Röckmann, T.: Stratospheric Air Sub-sampler (SAS) and its application to analysis of Δ17O(CO2) from small air samples collected with an AirCore, Atmos. Meas. Tech., 9, 5607–5620,, 2016. a

Murnick, D. E. and Peer, B. J.: Laser-based analysis of carbon isotope ratios, Science, 263, 945–947, 1994. a

Nakamichi, S., Kawaguchi, Y., Fukuda, H., Enami, S., Hashimoto, S., Kawasaki, M., Umekawa, T., Morino, I., Suto, H., and Inoue, G.: Buffer-gas pressure broadening for the (3001)III (0 0 0) band of CO2 measured with continuous-wave cavity ring-down spectroscopy, Chem. Phys., 8, 364–368,, 2006. a

Nara, H., Tanimoto, H., Tohjima, Y., Mukai, H., Nojiri, Y., Katsumata, K., and Rella, C. W.: Effect of air composition (N2, O2, Ar, and H2O) on CO2 and CH4 measurement by wavelength-scanned cavity ring-down spectroscopy: calibration and measurement strategy, Atmos. Meas. Tech., 5, 2689–2701,, 2012. a

Neubert, R. E., Spijkervet, L. L., Schut, J. K., Been, H. A., and Meijer, H. A. J.: A computer-controlled continuous air drying and flask sampling system, J. Atmos. Ocean. Tech., 21, 651–659,<0651:ACCADA>2.0.CO;2, 2004. a

Pataki, D. E., Bowling, D. R., and Ehleringer, J. R.: Seasonal cycle of carbon dioxide and its isotopic composition in an urban atmosphere: Anthropogenic and biogenic effects, J. Geophys. Res., 108, 4735,, 2003. a

Prokhorov, I., Kluge, T., and Janssen, C.: Laser Absorption Spectroscopy of Rare and Doubly Substituted Carbon Dioxide Isotopologues, Anal. Chem., 91, 15491–15499,, 2019. a

Roeloffzen, J. C., Mook, W. G., and Keeling, C. D.: Trend and variations in stable carbon isotopes of atmospheric carbon dioxide, Stable isotopes in plant nutrition, soil fertility and environmental studies, IAEA, Vienna, 601–618, 1991. a

Rothman, L. S., Gordon, I. E., Babikov, Y., Barbe, A., Chris Benner, D., Bernath, P. F., Birk, M., Bizzocchi, L., Boudon, V., Brown, L. R., Campargue, A., Chance, K., Cohen, E. A., Coudert, L. H., Devi, V. M., Drouin, B. J., Fayt, A., Flaud, J. M., Gamache, R. R., Harrison, J. J., Hartmann, J. M., Hill, C., Hodges, J. T., Jacquemart, D., Jolly, A., Lamouroux, J., Le Roy, R. J., Li, G., Long, D. A., Lyulin, O. M., Mackie, C. J., Massie, S. T., Mikhailenko, S., Müller, H. S., Naumenko, O. V., Nikitin, A. V., Orphal, J., Perevalov, V., Perrin, A., Polovtseva, E. R., Richard, C., Smith, M. A., Starikova, E., Sung, K., Tashkun, S., Tennyson, J., Toon, G. C., Tyuterev, V. G., and Wagner, G.: The HITRAN2012 molecular spectroscopic database, J. Quant. Spectrosc. Ra., 130, 4–50,, 2013. a, b, c

Rousseeuw, P. J. and Verboven, S.: Robust estimation in very small samples, Comput. Stat. Data Anal., 40, 741–758,, 2002. a, b

Sakai, S., Matsuda, S., Hikida, T., Shimono, A., McManus, J. B., Zahniser, M., Nelson, D., Dettman, D. L., Yang, D., and Ohkouchi, N.: High-Precision Simultaneous 18O/16O, 13C/12C, and 17O/16O Analyses for Microgram Quantities of CaCO3 by Tunable Infrared Laser Absorption Spectroscopy, Anal. Chem., 89, 11846–11852,, 2017. a

Santrock, J., Studley, S. A., and Hayes, J. M.: Isotopic Analyses Based on the Mass Spectra of Carbon Dioxide, Anal. Chem., 57, 1444–1448,, 1985. a

Stoltmann, T., Casado, M., Daëron, M., Landais, A., and Kassi, S.: Direct, Precise Measurements of Isotopologue Abundance Ratios in CO2 Using Molecular Absorption Spectroscopy: Application to Δ17O, Anal. Chem., 89, 10129–10132,, 2017.  a

Sturm, P., Leuenberger, M., Sirignano, C., Neubert, R. E. M., Meijer, H. A. J., Langenfelds, R., Brand, W. A., and Tohjima, Y.: Permeation of atmospheric gases through polymer O-rings used in flasks for air sampling, J. Geophys. Res.-Atmos., 109, 1–9,, 2004. a

Sturm, P., Eugster, W., and Knohl, A.: Eddy covariance measurements of CO2 isotopologues with a quantum cascade laser absorption spectrometer, Agr. Forest Meteorol., 152, 73–82,, 2012. a

Tans, P. P., Crotwell, A. M., and Thoning, K. W.: Abundances of isotopologues and calibration of CO2 greenhouse gas measurements, Atmos. Meas. Tech., 10, 2669–2685,, 2017. a

Trolier, M., White, J. W., Tans, P. P., Masarie, K. A., and Gemery, P. A.: Monitoring the isotopic composition of atmospheric CO2: Measurements from the NOAA global air sampling network, J. Geophys. Res.- Atmos., 101, 25897–25916,, 1996. a, b

Tuzson, B., Mohn, J., Zeeman, M. J., Werner, R. A., Eugster, W., Zahniser, M. S., Nelson, D. D., McManus, J. B., and Emmenegger, L.: High precision and continuous field measurements of δ13C and δ18O in carbon dioxide with a cryogen-free QCLAS, Appl. Phys. B, 92, 451–458,, 2008. a, b, c

Vogel, F., Huang, L., Ernst, D., Giroux, L., and Worthy, D.: Evaluation of a cavity ring-down spectrometer for in situ observations of 13CO2, Atmos. Meas. Tech., 6, 301–308,, 2013. a

Wehr, R., Munger, J. W., Nelson, D. D., McManus, J. B., Zahniser, M. S., Wofsy, S. C., and Saleska, S. R.: Long-term eddy covariance measurements of the isotopic composition of the ecosystem-atmosphere exchange of CO2 in a temperate forest, Agr. Forest Meteorol., 181, 69–84,, 2013. a

Wen, X. F., Meng, Y., Zhang, X. Y., Sun, X. M., and Lee, X.: Evaluating calibration strategies for isotope ratio infrared spectroscopy for atmospheric 13CO2/12CO2 measurement, Atmos. Meas. Tech., 6, 1491–1501,, 2013. a, b, c, d

Wendeberg, M., Richter, J. M., Rothe, M., and Brand, W. A.: Jena Reference Air Set (JRAS): A multi-point scale anchor for isotope measurements of CO2 in air, Atmos. Meas. Tech., 6, 817–822,, 2013. a, b

Zhou, L., Conway, T. J., White, J. W. C., Mukai, H., Zhang, X., Wen, Y., Li, J., and Macclune, K.: Long-term record of atmospheric CO2 and stable isotopic ratios at Waliguan Observatory: Background features and possible drivers, 1991–2002, Global Biogeochem. Cy., 19, 1–9,, 2005. a


Note that the r628 and r627 differ strictly speaking from the isotope ratios (r18 and r17) by a factor of 2.

Short summary
For understanding the sources and sinks of atmospheric CO2, measurement of stable isotopes has proven to be highly valuable. We present a new method using laser absorption spectroscopy to simultaneously conduct measurements of three CO2 isotopes, directly on dry-air samples. This new method reduces sample preparation time significantly, compared to the conventional method in which measurements are conducted on pure CO2, and avoids measurement biases introduced by CO2 extraction.