Dynamic–gravimetric preparation of metrologically traceable primary calibration standards for halogenated greenhouse gases
- 1Federal Institute of Metrology METAS, 3003 Berne, Switzerland
- 2Laboratory for Air Pollution and Environmental Technology, Empa, Swiss Federal Laboratories for Materials Science and Technology, Dübendorf, Switzerland
- anow at: Federal Office for the Environment FOEN, Air Pollution and Chemicals Division, 3003 Berne, Switzerland
Correspondence: Myriam Guillevic (firstname.lastname@example.org)
For many years, the comparability of measurements obtained with various instruments within a global-scale air quality monitoring network has been ensured by anchoring all results to a unique suite of reference gas mixtures, also called a “primary calibration scale”. Such suites of reference gas mixtures are usually prepared and then stored over decades in pressurised cylinders by a designated laboratory. For the halogenated gases which have been measured over the last 40 years, this anchoring method is highly relevant as measurement reproducibility is currently much better (< 1 %, k = 2 or 95 % confidence interval) than the expanded uncertainty of a reference gas mixture (usually > 2 %). Meanwhile, newly emitted halogenated gases are already measured in the atmosphere at pmol mol−1 levels, while still lacking an established reference standard. For compounds prone to adsorption on material surfaces, it is difficult to evaluate mixture stability and thus variations in the molar fractions over time in cylinders at pmol mol−1 levels.
To support atmospheric monitoring of halogenated gases, we create new primary calibration scales for SF6 (sulfur hexafluoride), HFC-125 (pentafluoroethane), HFO-1234yf (or HFC-1234yf, 2,3,3,3-tetrafluoroprop-1-ene), HCFC-132b (1,2-dichloro-1,1-difluoroethane) and CFC-13 (chlorotrifluoromethane). The preparation method, newly applied to halocarbons, is dynamic and gravimetric: it is based on the permeation principle followed by dynamic dilution and cryo-filling of the mixture in cylinders. The obtained METAS-2017 primary calibration scales are made of 11 cylinders containing these five substances at near-ambient and slightly varying molar fractions. Each prepared molar fraction is traceable to the realisation of SI units (International System of Units) and is assigned an uncertainty estimate following international guidelines (JCGM, 2008), ranging from 0.6 % for SF6 to 1.3 % (k = 2) for all other substances. The smallest uncertainty obtained for SF6 is mostly explained by the high substance purity level in the permeator and the low SF6 contamination of the matrix gas. The measured internal consistency of the suite ranges from 0.23 % for SF6 to 1.1 % for HFO-1234yf (k=1). The expanded uncertainty after verification (i.e. measurement of the cylinders vs. each others) ranges from 1 to 2 % (k = 2).
This work combines the advantages of SI-traceable reference gas mixture preparation with a calibration scale system for its use as anchor by a monitoring network. Such a combined system supports maximising compatibility within the network while linking all reference values to the SI and assigning carefully estimated uncertainties.
For SF6, comparison of the METAS-2017 calibration scale with the scale prepared by SIO (Scripps Institution of Oceanography, SIO-05) shows excellent concordance, the ratio METAS-2017 / SIO-05 being 1.002. For HFC-125, the METAS-2017 calibration scale is measured as 7 % lower than SIO-14; for HFO-1234yf, it is 9 % lower than Empa-2013. No other scale for HCFC-132b was available for comparison. Finally, for CFC-13 the METAS-2017 primary calibration scale is 5 % higher than the interim calibration scale (Interim-98) that was in use within the Advanced Global Atmospheric Gases Experiment (AGAGE) network before adopting the scale established in the present work.
Since the 1970s, atmospheric measurements of CFCs (chlorofluorocarbons) and HCFCs (hydrochlorofluorocarbons), used as refrigerants and blowing agents, have evidenced their role in stratospheric ozone layer depletion (Molina and Rowland, 1974; WMO, 1981). The reduction of CFC use (and later HCFCs) has been under strict regulations of the Montreal Protocol on Substances that Deplete the Ozone Layer since entering into force in 1989 (WMO, 2014). While the molar fractions of major CFCs are now declining in the atmosphere, some longer-lived minor CFCs are still increasing (Laube et al., 2014; Vollmer et al., 2018; WMO, 2014). HFCs (hydrofluorocarbons) were introduced as replacement for CFCs and HCFCs. Their emissions, though not harmful to the ozone layer, are still increasing and contributing to global warming due to their high radiative forcing (Harris and Wuebbles, 2014; Velders et al., 2009). For this reason the recent Kigali Amendment (October 2016) added these HFCs to the Montreal Protocol. As a consequence and in compliance with European Directive 2006/40/EC, replacement compounds were introduced, foremost the HFOs (Vollmer et al., 2015).
Non-Article 5 (developed) countries of the Montreal Protocol are bound to report their CFC, HCFC and HFC production and consumption based on so-called “bottom-up” inventories. Within the Kyoto Protocol, aiming at limiting climate change, Annex-1 (developed) countries report inventories of greenhouse gas emissions to the atmosphere, for which reduction targets are set. Reviews of the success of such reduction targets and projections of future developments are so far still based on those inventories, while ozone layer recovery and climate change depend on real atmospheric molar fractions.
Atmospheric measurements of halogenated compounds are currently provided by several networks such as AGAGE (Advanced Global Atmospheric Gases Experiment), NOAA (National Oceanic and Atmospheric Administration) and GAW (Global Atmospheric Watch). Such measurements, used to precisely estimate atmospheric molar fraction of these halogenated substance together with associated trends, are crucial to understand and predict the evolution of stratospheric ozone and estimate their radiative forcing thereby refining future climatic projections. Furthermore, based on these measurements and using atmospheric transport modelling, emissions can be quantified (Brunner et al., 2017; Prinn et al., 2000; Rigby et al., 2010). The comparison of top-down reconstructions with bottom-up inventories shows agreement for some gases but also discrepancies that can be considerable for others (Hu et al., 2016; Lunt et al., 2015; Sherry et al., 2017; Simmonds et al., 2016; Weiss and Prinn, 2011). The top-down approach thus is a complementing and independent way to review production–consumption–emission inventories and compliance with reduction targets, while assessments of climate forcing and stratospheric ozone rely on observations of atmospheric composition.
To quantify atmospheric molar fractions at specific sites, detect gradients between monitoring stations, evaluate data consistency and lack of biases and thereby attribute emissions to specific geographical areas, results have to be linked to accurate calibration scales. Currently for most monitored halogenated gases, measurement precision is as low as 0.4 % (Miller et al., 2008). This is much better than the expanded uncertainty of reference gas mixtures, estimated to be on the order of 2 to 4 % for SIO standards (Prinn et al., 2000), 0.3 to 3 % for NIST (Rhoderick et al., 2015) and 0.6 to < 2 % for NOAA (Hall et al., 2007; Lim et al., 2017; Montzka et al., 2015), as examples. Given this technical state of the art, monitoring networks have developed a so-called “primary calibration scale” system, in which all stations of a specific network are anchored to the same suite of primary reference gas mixtures with a calibration chain as short as possible. To ensure long-term continuity of atmospheric composition data, such calibration scales have to be maintained over years and only replaced if substantial necessity arises, with a scheme to properly back-calibrate all past measurements if a new calibration scale is defined. Accurate and stable calibration scales directly impact the quality of emission quantification, trend estimation and the calculation of values relevant for future climate projections, such as radiative forcing, which is derived from atmospheric composition.
The calibration scale approach enables a high degree of consistency, but still requires detecting and documenting systematic offsets as an indicator for potentially existing systematic biases. This QA–QC procedure includes regular intercomparisons between instruments installed at the same sites and reference gas mixtures or air sample flask exchanges (Hall et al., 2014; Rhoderick et al., 2015). In addition, geographically close monitoring stations anchored to different calibration scales regularly assess the potential development of a bias over time, which can dependent on the molar fraction (O'Doherty et al., 2004; Rigby et al., 2010; Simmonds et al., 2017; Vollmer et al., 2016).
The production of a robust and accurate primary calibration scale at low molar fraction level, which is usually created through a set of reference gas mixtures, is a major task (Dlugokencky et al., 2005; Hall et al., 2007; Prinn et al., 2000; Zhao and Tans, 2006; Zhao et al., 1997). In particular for compounds prone to adsorption on material surfaces, it is difficult to evaluate mixture stability and thus variations of molar fractions over time in cylinders. As a consequence, some halogenated gases are measured in the atmosphere against reference standards that lack conventional calibration (relative calibration scales) or with only limited sets of primary reference standards (Vollmer et al., 2015).
To support atmospheric monitoring of halogenated gases, we present here a method to produce reference gas mixtures at near-atmospheric molar fractions for SF6 (sulfur hexafluoride), HFC-125 (pentafluoroethane), HFO-1234yf (HFC-1234yf, 2,3,3,3-tetrafluoroprop-1-ene, newly emitted compound), HCFC-132b (1,2-dichloro-1,1-difluoroethane) and CFC-13 (chlorotrifluoromethane). While SF6 and HFC-125 have already widely used calibration scales (see Sect. 4.1), HCFC-132b and CFC-13 have been measured for many years in AGAGE yet not reported on a conventional calibration scale (Vollmer et al., 2018). We apply a preparation technique combining dynamic gravimetry (ISO 6145-10, 2002) and dynamic dilution (ISO 6145-7, 2009), followed by cryo-filling in cylinders. This is an alternative to the prevailing preparation method for such compounds, which is static gravimetry (Hall et al., 2007; Lim et al., 2017; Prinn et al., 2000; Rhoderick et al., 2015). The produced suite of reference gas mixtures is SI-traceable; i.e. all measured or relevant quantities of each preparation step are linked to the realisation of SI units (the International System of Units) through an unbroken chain of calibrations. Finally, the prepared suite of mixtures is assigned a carefully quantified expanded uncertainty, taking into account all known relevant potential sources of uncertainties, following international recommendations (JCGM, 2008). The aim is to provide independent calibration scales and compare them to other available scales. This would provide existing calibration scales, which have so far been used on a relative basis, with a link to the SI.
In this paper, we present in detail the method developed to prepare SI-traceable reference gas mixtures for the mentioned halogenated gases at near-atmospheric background levels, i.e. as low as 1 pmol mol−1 (Sect. 2). The calculations to assign prepared values are described together with the method. The associated uncertainty budgets, established following JCGM (2008), are presented and discussed in Sect. 2.5. The internal consistencies of the calibration scales are determined in Sect. 3.2. These SI-traceable primary standards are then compared to other standards currently in use in the AGAGE network, as reported in Sect. 4. In addition, we determine conversion factors which compare the METAS-2017 suite to other primary calibration scales.
2.1 Method overview: dynamic–gravimetric generation process
The suite is designed to consist of one master cylinder (hereafter MP-001, for METAS primary cylinder no. 001) containing all components at near-ambient molar fractions but none below 1 pmol mol−1, in order to not exceed preparation uncertainties of 2 %. Ten additional cylinders are filled with molar fractions bracketing those filled in cylinder MP-001 over the range from 20 % less to 30 % more. The resulting prepared molar fraction range covered by this suite varies between the five compounds, with a range of 0.9–1.5 pmol mol−1 for HFO-1234yf with the lowest molar fractions to 26–42 pmol mol−1 for HFC-125 with the highest molar fractions (see details for each substance in Table 3). This allows the later determination of the internal consistency of the suite.
The generation process consists of three successive steps, starting from pure halocarbon substances diluted to molar fractions as low as pmol mol−1 in synthetic air. In a first step a matrix gas is spiked with one pure halocarbon substance using a permeation device (Sect. 2.2). In a second step, the high molar fraction mixture is dynamically diluted to pmol mol−1 level mixture using thermal mass flow controllers (MFCs, Sect. 2.3). In the final step the mixture is successively transferred into the 11 cylinders by cryo-filling (Sect. 2.4). In order to generate multi-component mixtures, the permeation device is changed and all steps are repeated with the mixture containing a new substance being added to the same suite of cylinders.
Reference gas generation by applying the permeation method combined with dynamic dilution is an established, standardised technique, particularly for reactive gases (Brewer et al., 2011; Flores et al., 2012; Haerri et al., 2017; O'Keeffe and Ortman, 1966; Scaringelli et al., 1970). It is routinely used at METAS following a procedure in compliance with international standards (ISO 6145-10, 2002; ISO 6145-7, 2009). The permeation method is based on constant transfer of the substance of interest from a permeation device (or permeator), resulting in a mass loss which can be continuously monitored. Permeation devices are placed in a permeation chamber, i.e. a controlled atmosphere in terms of temperature and pressure, continuously flushed by a carrier gas stream. Permeators used here consist of a stainless-steel reservoir containing the pure substance as a liquid, sealed with a cap containing a polymer membrane permeable for the specific substance. All permeators used in this study (Fine Metrology, Italy) were filled with substances of purity 99 % or higher (Synquest Laboratories, Florida, USA). Substance purity levels were determined by the manufacturer using flame ionisation detector analysis. Permeators are filled under laminar flow to avoid potential contamination.
The permeation rate depends exponentially on temperature; secondary influences are carrier gas composition and pressure (Brito and Zahn, 2011; Haerri et al., 2017; Jost, 2004; Lucero, 1971; Moosbach and Hartkamp, 1993). For mass loss determination, each permeator is placed individually in our magnetic suspension balance (MSB) “Violetta” (hereafter MSB-Violetta, model FLUIDIFF MP, installed in 2014, Rubotherm, Germany, Fig. 1). A MSB system allows for continuous and unperturbed mass measurements as the temperature-controlled chamber, where the permeator is suspended, is physically decoupled from the balance itself. The stainless-steel permeation chamber of MSB-Violetta allows for precisely controlling temperature (±0.02 ∘C, measured with a Pt100 sensor), pressure (±0.1 hPa, Bronkhorst El-Press, the Netherlands) and carrier gas flow (±0.1 % of the flow, Red-y series, Vögtlin, Switzerland). In this study, permeation typically occurs at 45 ∘C and 3500 hPa. Inside the chamber, the permeator is coupled to a permanent magnet also placed in the chamber. On the outside of the chamber, an electromagnet connected to the balance (Sartorius ME66S) exerts a force over the permanent magnet thus coupling the permeation unit to the balance. To minimise balance noise, total weight and position of the permeation device are adjusted, as well as the vertical position of the electromagnet. The absolute mass of the permeation unit is measured in 3-minute intervals. To correct for balance drift, after every three measurement points a calibration mass (CM1) is automatically placed on the balance plate. Additionally, a second calibration mass (CM2) of same volume but different mass is measured every six calibration points to correct for any potential buoyancy change affecting the measurement of CM1. The masses of CM1 and CM2 are traceable to the Swiss realisation of the kilogram (Fuchs et al., 2012; Marti, 2017; Marti et al., 2015). All weight measurements and associated corrections realised with MSB-Violetta are fully automated. After each opening and closing of the permeation chamber, its tightness is checked by closing the input gas from the synthetic air cylinder as well as all exits and checking the absence of pressure decrease over time.
After inserting a permeation device, a stabilisation period is required mainly depending on chamber temperature, pressure and permeator membrane properties, before the mass loss becomes linear over time. This linear mass loss vs. time is then determined for at least 8000 min to minimise the standard deviation of the measured mass loss due to balance noise. The time window during which the mass data are used is determined so that the residuals of the fit to the mass loss over time are centred around zero and randomly distributed (see example for CFC-13 in Fig. 2 and Sect. S1 in the Supplement). For each substance i, the molar fraction xperm,i of the mixture exiting the permeation chamber can be calculated as
, the mass difference between the beginning and end of linear mass loss (g);
, corresponding time difference (min);
purityi, purity fraction of the substance filled in permeator (mol mol−1);
and Vm,carrier, molar volume of carrier gas, here synthetic air, (L mol−1). All volumes in this work are given at standard temperature and pressure (STP), i.e. 0 ∘C and 1013.25 hPa. Values are from the NIST Chemistry WebBook, assuming real gas;
qV1,i, the volumetric flow of carrier gas regulated by MFC1 (Fig. 1, L min−1);
and xres,i, the residual molar fraction of substance i in carrier gas (pmol mol−1).
Permeation rates were determined at 3500 hPa and temperatures between 36 and 45 ∘C. We observed a particularly long stabilisation time (i.e. 6000 min) for the permeator containing CFC-13 (Fig. 2). The low mass loss for HFO-1234yf (88.5 ng min−1) required a continuous measurement during 22 days to reach the required standard deviation (< 0.4 % for Stabbalance,i).
2.3 Dynamic dilution
The mixture exiting the permeation chamber, at µmol mol−1 level, is diluted to pmol mol−1 levels over two successive, dynamic dilution steps (Fig. 1). The flows are piloted by thermal MFCs (ISO 6145-7, 2009). First, the mixture exiting the permeation chamber is diluted with a first dilution flow (Fig. 1 MFC2, up to 5 L min−1). The total flow passing through the permeation chamber (MFC1) and diluting this flow (MFC2) is measured by a mass flow meter (MFM in Fig. 1, Vögtlin, Switzerland), so that only this MFM needs calibration to calculate the resulting molar fraction. After this first dilution step, xperm,i can be rewritten as
Second, a small flow of this resulting mixture is sampled by MFC3 (10 mL min−1) and is diluted by another larger flow (MFC4, up to 5 L min−1). The dilution factor fdilution,i to obtain the second dilution step is calculated as
with qVn,i, the volumetric flow of MFCn (L min−1), set for substance i.
After this second dilution step, the prepared molar fraction xfilled,i that will be filled in all cylinders can be calculated as
At this stage, the generated mixture has a molar fraction approximately 10 times higher than atmospheric levels. Note that most metal surfaces in contact with the carrier gas and the produced gas mixture are passivated by applying SilcoNert2000 coating. This includes all metal tubing, all metal surfaces of the MFCs and MFM in contact with the gas and most of the permeation chamber.
The generated mixture after the second dilution step is then transferred into 11 cylinders, named MP-001 to MP-011. The technical characteristics of these cylinders are summarised in Table 1. In brief, we used seven cylinders made of electropolished stainless steel (Essex Industries, USA) and four cylinders made of Silconert2000-coated stainless steel (Swagelok). To detect potential systematic biases due to adsorption on cylinder surfaces, this set of cylinders presents four different surface ∕ volume ratios (see details in Table 1).
A set flow of the reference gas mixture (qV5,j, 3 L min−1 for the 34 L Essex cylinders, 0.5 L min−1 for all other smaller volume cylinders) is sampled by MFC5 (CMOSsens series, Sensirion, Switzerland, Fig. 1) and then directed to a tee. The two exiting paths of this tee are piloted by pneumatic valves (Swagelok), one being normally closed (to the cylinder, V1) and the other one being normally open (to the pump, V2).
All cylinders are cleaned beforehand, being evacuated to approximately 6 Pa and filled with nitrogen at 2000 hPa (purity grade 99.999 % or better), three times. Each cylinder was then evacuated one last time to 6 Pa or lower just before being connected to the filling system.
After being connected, the tubing between the cylinder valve and the pneumatic valve V1 (Fig. 1) is evacuated to 1 hPa and filled with synthetic air, three times. After leak checking, it is filled with synthetic air at a pressure of 800 hPa ± 20 hPa, because it is the residual pressure that is lost after each filling, when the cylinder valve is manually closed (see Supplement Sect. S2).
After this preparation step, the bath is (re)filled with liquid nitrogen. Once the cylinder content is re-liquified, its valve is manually opened, and, piloted by a LabVIEW program, the pneumatic valve positions are switched for a precisely set duration in order to fill a precisely controlled volume of the reference gas mixture in each cylinder. The filling durations last from 17 min to 6 h depending on substance and target molar fraction. To avoid freezing, the cylinder valve is intermittently heated with a heat gun during filling. Once the pneumatic valves are back at their default setting, the cylinder valve is kept open for 1 min, while the tubing between the cylinder valve and V1 is heated to force all potential remaining substance of interest to be cryo-trapped in the cylinder still in liquid nitrogen. The cylinder valve is then manually closed, and the cylinder is disconnected and left standing vertically outside the building to warm up. A new cylinder is placed in the liquid nitrogen bath, connected to the filling system, and the filling procedure starts again. The state of the filling system is checked at least every 30 min during filling. Once the filling of each cylinder with one given mixture is completed, the permeation chamber of the balance is opened, a new permeator installed, the chamber closed, and a new mass loss measurement starts.
After the fillings of all five substances of interest are completed, additional precisely determined volumes of synthetic air are cryo-filled into the cylinders in order to dilute the mixtures to atmospheric molar fractions. For this special filling, synthetic air is also directed through the permeation chamber where a glass reaction tube filled with deionised water is placed, evaporating water, in order to slightly humidify each cylinder. The resulting water vapour molar fraction ranges from 20 to 70 µmol mol−1. Note that this added water is not included in the calculations of molar fractions for the halogenated compounds, which are therefore expressed in dry synthetic air. Adding water vapour to each cylinder was motivated by the fact that losses due to adsorption are known to occur for some halogenated compounds (Prinn et al., 2000). This has been evidenced by Yokohata et al. (1985) for CCl4 and CH3CCl3 in very dry mixtures (i.e. likely less than 750 nmol mol−1 of water vapour), who also experimentally showed that adding water vapour to the cylinder annulled this adsorption, water vapour being an excellent competitor for adsorption sites on a metal surface (Pogàny et al., 2016; Vaittinen et al., 2013).
Reference gas mixtures for each of the five compounds were filled in this order: SF6, HFC-125, HFO-1234yf, HCFC-132b and CFC-13. All cylinders are homogenised for a minimum of 6 h each before measurement using an automated rolling system, with cylinders lying horizontally, and alternating directions. A minimum of 1 week elapsed between the final cylinder fillings and the measurements.
Filling different mixtures successively in each cylinder j with different filling duration results in an additional dilution factor for each substance i. The durations are chosen in order to attain the same dilution factor for all substances in one given cylinder. The resulting molar fraction for substance i in cylinder j is
with the total filling duration in each cylinder j being
where Δti,j is the filling duration of mixture containing substance i in cylinder j and qV5,j is the flow of MFC5 used to regulate the flow into the cylinder during cryo-filling. Equation (5) can be simplified by removing qV5,j and re-arranging into
However, the stability component of the flow Stab has to be taken into account in the uncertainty budget (see Eq. 9).
Note that for a given substance, the molar fractions in cylinders vary only due to varying filling durations. This is designed to maximise the correlation between cylinders and therefore improve the resulting internal consistency (defined in Sect. 3.2.3) of the prepared suite of reference gas mixtures. Correlation coefficients between cylinders for one given substance therefore range from 0.96 to 0.99. This set constitutes the METAS-2017 primary calibration scale for SF6, HFC-125, HFO-1234yf, HCFC-132b and CFC-13. All computed prepared values are reported in Table 3.
2.5 Uncertainty of preparation
2.5.1 Uncertainty of prepared values in pmol mol−1
We estimate the uncertainty of the assigned molar fraction, for each substance in each cylinder, following JCGM (2008) by measuring (type A uncertainty) or assigning (type B uncertainty) an uncertainty estimate to each input quantity used in the equations presented in Sect. 2. Expanded uncertainties, noted U, are then calculated with k = 2 for a coverage probability of approximately 95 %. All uncertainty computations for the prepared molar fractions are made using the GUMWorkBench program. We describe hereafter the most important contributions to the uncertainty.
Balance measurement. The uncertainty of the weighing is estimated by fitting a linear function through the measured masses over time and using the standard deviation between this linear fit and the residuals for each point as uncertainty estimate (k = 1, e.g. Fig. 2). This measurement noise level, noted Stabbalance,i in Eq. (8), is on the order of 0.2 % of the mass loss (for SF6, whose mass loss was higher making the associated balance noise relatively smaller) to 0.6 %. In addition, we take into account potential biases due to mass calibration of the two used calibration masses (CM1, CM2) and time uncertainty, even if these two contributions are extremely small (see also Sect. S1 in the Supplement).
Permeation chamber temperature stability StabT. Once carrier gas flow and pressure are kept constant, the permeation rate varies only with temperature. The stability of the permeation chamber temperature is 0.02 ∘C over 20 min (k = 2). Based on our experience measuring temperature sensitivity of permeation rate, this corresponds to approximately 0.1 % change in permeation rate. Stabbalance,i and StabT are given a value of 1 and included in Eq. (2) in order to take into account their uncertainty contributions:
Substance purity purityi. We use the certificate provided by the substance manufacturer, i.e. purityi=0.999 for SF6 and purityi=0.99 for all other substances. As uncertainty we choose a triangular distribution with 1 as upper boundary, and 1–2 ⋅ (1−purityi) as lower boundary. This is a conservative approach.
Impurity in the carrier gas xres,i. For all fillings, a total of three 50 L cylinders of the same type of synthetic air (Pangas synthetic air 5.6) have been used. Absolute impurity levels in one cylinder were measured at Empa on the Medusa GC-MS system (Miller et al., 2008). The similarity of the impurity level in the other cylinders was checked at METAS on a similar preconcentration GC-MS system by trapping a total of 6 L of synthetic air. For each measured impurity molar fraction xres,i (ranging from 6 to 30 fmol mol−1), we use a triangular distribution centred in xres,i, with zero as lower boundary and as upper boundary (see Table S4).
Calibrated values of volumetric flows. , qV3,i and qV4,i (MFM, MFC3 and MFC4 in Fig. 1) were calibrated using a SI-traceable primary reference standard applying a system of pistons with known volume (Niederhauser and Barbe, 2002). All flows are given at standard temperature and pressure (0 ∘C, 1013.25 hPa). The actual flow has been calibrated using the same gas type, at the same input parameter values (pressure set points before and after each MFM/MFC, MFC set points). A minimum of four replicate measurements for each flow set point have been made and the average was taken as best estimate. The stability of each flow over the measurement duration (1–5 min) is better than 1 ‰. The obtained relative expanded uncertainty (k = 2) for each flow is U=0.3 % for MFC3 (due to smaller flows) and U=0.2 % for MFM and MFC4.
Filling durations Δti. Individual filling durations range from 1060 to 24 110 s (18 min to 6.7 h, see Supplement Table S3). The associated uncertainty is fixed at U = 1.8 s and takes into account the response time of each pneumatic valve as well as the computer time clock uncertainty. In percentage this uncertainty therefore decreases with increasing filling duration.
Stability of filling flow Stab. The flow stability of MFC5 was challenged due to the relatively small pressure gradient of approximately 800 hPa. We therefore take into account a flow stability component in the uncertainty of U=0.1 %. This stability component plays a role to explain the internal consistency between cylinders (Sect. 3.2.3).
The resulting, combined uncertainty is then expanded using a coverage factor k = 2 (representing a 95 % confidence interval for a Gaussian distribution). The obtained values are documented in Table 3. As an example for cylinder MP-001, we present in Fig. 3 the uncertainty contribution for the most important contributors as pie charts, for each substance. Expanded uncertainties range from 0.6 % (for SF6) to 1.3 %. The smallest uncertainty obtained for SF6 is mostly explained by the high substance purity level inside the permeator (0.999 pure, 10 times better than for the other substances), as well as low SF6 contamination of the carrier gas.
2.5.2 Uncertainty of prepared ratios
To calculate the expected internal consistency of the prepared suite of mixtures, we also calculate ratios of assigned values, using cylinder MP-001 as master cylinder. The assigned value of a ratio for substance i between cylinder j and cylinder MP-001 (marked by subscript 1 hereafter) can be calculated with a very good approximation by
The terms “Stab” represent the stability of each filling temperature and each filling flow of MFC5, respectively. Each “Stab” term is assigned a value of 1 and an expanded uncertainty of Stab.1 % for the temperature stability and Stab.1 % for the flow stability (as discussed in Sect. 2.5.1). The standard uncertainty (k = 1) of , calculated using Eq. (9), is hereafter noted .
Ratio values range from 0.8 to 1.3 and the corresponding expanded uncertainty has an actually constant value of 0.3 % over this limited ratio range (Table 3). Due to the elimination of many common factors when working in such a ratio space, the correlation coefficient between ratios is approximately 0.4.
3.1 Measurement method
The relative molar fractions of the five compounds in the 11 samples MP-001–MP-011 were determined using Medusa gas chromatography mass spectrometry (GC-MS) methods (Arnold et al., 2012; Miller et al., 2008). Medusa GC-MS systems have been in use in AGAGE for hourly field and laboratory measurements of more than 50 halogenated compounds. In brief, it consists of a multi-port inlet, a sample drying system (Nafion driers), a custom-designed preconcentration system, a capillary GC column (CP-PoraBOND Q, 0.32 mm ID × 25 m, 5 µm film thickness, Varian Chrompack) and a quadrupole mass spectrometer (Agilent Technologies, 5975).
Measurements were conducted on the Empa laboratory Medusa GC-MS and consisted of 2 L samples measured alternatingly with a reference gas measurement. MP-001 was used as the reference sample. Hence all results of the MP-002–MP-011 samples are expressed as ratios by dividing the chromatographic peak area of the sample by the mean of the area of the bracketing MP-001 measurements. Minor analytical modifications compared to the routine field measurements were adopted. In particular, to enhance precision of the measurements, compounds, which chromatographically elute near those of interest, were omitted from acquisitions, and the electron multiplier was elevated to enhance the signal size. Integration of the chromatographic peaks was done using commercial software (GCWerks). Detection limits, defined here as 3 times the noise level, are 0.015 pmol mol−1 for SF6, 0.02 pmol mol−1 for HFC-125, 0.01 pmol mol−1 for HFO-1234yf, 0.015 pmol mol−1 for HCFC-132b and 0.07 pmol mol−1 for CFC-13. All results are expressed as dry-air molar fraction. The system was shown to be free of blanks and linear in the range of molar fractions of our samples (Vollmer et al., 2015). The repeatability of the measurements was calculated as the standard deviation (1σ) of the area ratios. The standard deviation of the mean was obtained by dividing by the square root of n, n being the number of measurements.
3.2 Measurement results
Quantifying analyser response, correcting results and identifying potential outliers is done in an iterative manner. First, all measurement results with their standard deviation are used to calculate a first analyser response function. Using this function, all measurements results are corrected. The corrected values are then compared to the assigned value. Based on the results of this test, potential outliers are excluded, and an adjusted analyser response is calculated, etc. Details of the calculations are presented hereafter.
3.2.1 Analyser response calibration
Following the approach already developed at SIO, and similar to methods already used for isotopic studies (Craig, 1957; Dansgaard, 1953), we compare measured and assigned values in the ratio space, for three reasons. First, the mass spectrometer used for analysis is naturally drifting over time and to correct for this effect, measured areas are expressed as area ratios relative to the bracketing mixture used as standard, here MP-001. The most precise measurement result given by the MS is therefore a ratio of areas. Second, due to the preparation design made to maximise correlation between cylinders for one given substance, here again the prepared ratio has an uncertainty much smaller than each molar fraction separately, because when calculating these ratios many constant factors cancel out (Eq. 9). The uncertainty components that still have to be taken into account are mostly those related to the stability over time of the preparation system, i.e. the stability of the permeation temperature, the stability of each MFC flow (negligible, except for MFC5) and the filling duration uncertainty. Expanded uncertainties in ratio space are therefore 0.3 %, compared to expanded uncertainties of 0.6 to 1.3 % in the molar fraction space (Table 3). Third, the correlation between values of assigned ratios is much smaller (on the order of 0.4) than the correlation between assigned molar fraction values, and therefore using ratios is again better indicated to estimate the fit between assigned values and instrument response (ISO 6143, 2001).
To calibrate the measured value, we thus determine the analyser response function in ratio space, i.e. measured ratio vs. assigned ratio. We can therefore compare a maximum of 10 ratio values, using 11 cylinders. An example is given in Fig. 4.
The analyser response is calculated using a linear fit due to the relatively small number of measured values as well as the good linear response of the MS over this limited range:
where is the measured area ratio for substance i in cylinder j.
The fit coefficients ai and bi are computed using a bivariate weighted linear fit, following the York algorithm (York et al., 2004) as described in Cantrell (2008) (Supplement Sect. S3), coded in Octave (results in Table 4). As an additional test we ran the same fitting algorithm forcing bi=0. Interestingly, for each substance the obtained slope can then not be distinguished from a slope of 1 within uncertainties (Table 4). This suggests that the analyser response is linear within the tested range and within stated uncertainties. The deviation in the corrected values for cylinder MP-001 is no more than 0.2 %, but varies up to 0.5 % for MP-005 and 1 % for MP-006, which are at the upper and lower end of the scale. We also calculated ai and bi using a simpler weighted linear fit considering only weights associated with the measured values, giving very similar results (difference < 0.03 %).
The measured and corrected ratios and values in pmol mol−1 are calculated using
where is the prepared value in pmol mol−1 for substance i in cylinder MP-001.
Hereafter we refer to as measured ratios and to as measured molar fractions.
3.2.2 Verification test and exclusion of outliers
To compare the measured results to the prepared results, we use the measured ratios and their associated uncertainties calculated according to Eq. (11) (Supplement Table S5). We then calculate the verification criteria (ISO 6143, 2001):
Using this procedure, after the first iteration cylinder MP-008 is excluded as outlier for SF6, HFC-125 and HFO-1234yf, the measured value being systematically too low. Because this represents already the majority of all substances, we decided to exclude this cylinder for HCFC-132b and CFC-13 as well. After the second iteration cylinder MP-002 is excluded for HFC-125, MP-010 for HFO-1234yf and MP-003 for HCFC-132b. After the third iteration, cylinder MP-010 is identified as an outlier for HCFC-132b.
For cylinder MP-008 (a 4.5 L Essex cylinder filled at 24 bar, Table 1), we note that the difference between measured value and assigned value is large for SF6 (the first filled substance) and more or less decreases until showing no particular offset for CFC-13 (being filled last). An explanation for such a time varying offset could be a potential leak over time, affecting more the substance that was filled first.
During filling of cylinder MP-002 for HFC-125, we noted that the mass flow controller sampling the flow into the cylinder (MFC5, Fig. 1) was suffering from large instability during the 15 s towards the end of the filling. This particularly high instability was very likely due to both the small pressure gradient over MFC5 and the small flow from MFC4 (4 L min−1) compared to the flow sampled by MFC5 (3 L min−1), making the pressure control before MFC5 (done by PrC2) particularly challenging. To improve the pressure regulation by PrR2 and therefore limit MFC5 instability, an additional buffer volume (approximately 1 L, stainless steel) was added on the flow path just before PrC2. This instability did not occur during subsequent fillings.
We tentatively explain cylinder MP-003 being an outlier for HCFC-132b by the fact that this cylinder was filled first after a synthetic air cylinder exchange. Additional tests showed that HCFC-132b is affected by regulator contamination and takes time to clear out. The regulator was purged three times before using the synthetic air but perhaps this was not sufficient. The system being then continuously running, the purge was already sufficient for the second filling not to be affected.
We unfortunately do not have a specific explanation for cylinder MP-010 being an outlier for HCFC-132b and HFO-1234yf, potentially pointing towards operator error, and highlighting the need for a verification step. We observe, however, that most outliers are cases of substance loss (Fig. 5) and affect cylinders having the smallest total amount of gas filled and the highest surface / volume ratio (i.e. Essex cylinder 4.5 L, 24 bar). We would therefore in the future favour filling in cylinders of larger volume and pressure and further automatising the cryo-filling process to limit human intervention as much as possible and to increase the safety of the procedure. Furthermore, we plan to investigate whether the observed substance losses occurred by adsorption on cylinders walls or beforehand in the preparation system. To do so, a comparison of the molar fraction in the mixture exiting the magnetic suspension balance–dynamic dilution system with the same mixture filled in cylinders will be performed. The recent installation of a measurement system for halogenated gases at METAS in the same laboratory makes this possible. If the adsorption indeed occurs in the cylinders, we will test whether adding water vapour earlier in the sequence of fillings helps to limit this adsorption.
After excluding the outliers, we can observe no systematic bias between mixtures filled in electropolished cylinders from those filled in SilcoNert-coated cylinders. This result gives us confidence that there is also no identifiable loss of the halogenated substances on cylinder surfaces.
3.2.3 Internal consistency of the suite
To give an estimation of the preparation reproducibility, we calculate the so-called “internal consistency” of the suite of mixture (Prinn et al., 2000). This parameter quantifies the difference between assigned values and measured values for each substance.
First, we calculate the difference d for each substance, in each cylinder, as
and the associated uncertainty as
We then calculate the weighted mean difference for each substance, over the set of cylinders (without outliers):
We use the corresponding weighted standard deviation of dW,i as estimator of the internal consistency:
The calculated internal consistencies for each substance are reported in Table 4. This measured internal consistency is due to measurement reproducibility of the Medusa GC-MS system as well as how stable the preparation system is (i.e. including all potential sources of random noise in the preparation system, any systematic bias being on the contrary cancelled out when working in ratio space).
3.3 Final assigned uncertainties
3.3.1 According to ISO-6142-1:2015
This formula includes the uncertainty of the preparation as well as the uncertainty of the verification, with equal weights. The resulting uncertainties range from 0.7 % for SF6 to 1.5 % for HFO-1234yf (excluding outliers). In particular, here the uncertainties for HFC-125 and HCFC-132b become smaller (1 %) than the prepared uncertainty (1.3 %), because in this formula uncertainties from preparation and measurement are arbitrarily given equal weights.
3.3.2 According to preparation and measurement equations
We calculate the uncertainty of using the uncertainty of each component defining it, according to Eq. (12). This method therefore includes both the preparation uncertainty (through ) and the measurement uncertainty (through ). The expanded uncertainty (k = 2) ranges from 1 % for SF6 to 2 % for HCFC-132b and CFC-13 (Table 5). The disadvantage of using this method is that is included twice, once in through ai and bi and once in , which includes the same stability factors as in (Eqs. 2 and 7). To correct for this effect we would need to remove the stability factors in . However, the uncertainty budgets of the prepared mixtures (Fig. 3) suggest that this contribution is overall minor (from less than 1 to 5 % of the total), and removing one occurrence of would be a marginal change.
We favour this method to calculate the uncertainty, where the sensitivity of the final value to the preparation and verification uncertainties are taken into account, through Eq. (12), and where the resulting uncertainty is the largest – it is also the most conservative approach.
Cylinder MP-001 was used to compare the METAS-2017 calibration scale to other scales, using Empa's Medusa GC-MS as comparator. We report hereafter a brief description of each calibration scale followed by the results of the comparisons, for SF6, HFC-125, HFO-1234yf and CFC-13. For HCFC-132b, there was no other scale available for comparison.
4.1 Description of calibration scales
SIO-05. SIO developed in the 1990s a preparation method using a bootstrap technique to make a suite of standards for halogenated compounds at pmol mol−1 levels (Prinn et al., 2000; Weiss et al., 1981). This technique relies on the preparation of a first primary standard for CO2 at ppm levels (work of C. D. Keeling at SIO). CO2 is then used as bootstrap gas to prepare a second mixture containing CO2 and N2O with a known prepared N2O ∕ CO2 ratio, close to atmospheric conditions. The CO2 ∕ N2O ratio is prepared gravimetrically, by filling pure CO2 and N2O in individual glass vials, flame-sealing the vials, weighing them and mixing their content into an aliquot. In this second standard the N2O molar fraction is assigned using this known, gravimetrically prepared N2O ∕ CO2 ratio and a CO2 molar fraction calibration vs. the CO2 calibration scale. In a third step N2O is used as bootstrap gas to produce a new mixture containing N2O and halogenated gases, again with ratios close to ambient conditions. Knowing these halogen ∕ N2O ratios and by measuring the N2O molar fraction against the N2O primary standard, each halogen molar fraction is determined. This preparation system therefore combines gravimetric preparation of ratios and measurement of one of them vs. another suite of standards.
For comparison with the METAS-2017 calibration scale, we use so-called “tertiary tanks” filled with real air at SIO and calibrated vs. an SIO “secondary standard”, itself calibrated vs. the suite of gravimetrically prepared, primary reference gas mixture defining the calibration scale for each compound. It is therefore necessary to take into account the uncertainty due to scale propagation, i.e. the measurement repeatability of the secondary vs. the primary tank and of the tertiary vs. the secondary tank. The same principle applies to all other substances on a SIO calibration scale measured from a tertiary tank vs. the METAS-2017 scale.
When preparing the primary calibration scales for halogenated gases, SIO did not assign an uncertainty following JCGM (2008), but still the internal consistency of the scale was precisely determined as being u=0.4 % for SF6 (Ray F. Weiss, personal communication, October 2017). Regular comparison with NOAA, comparing results at co-located monitoring stations (Rigby et al., 2010) or through cylinder exchanges (Hall et al., 2014), shows agreement within 0.2 % or better for SF6. NOAA recently determined the uncertainty of its SF6 reference gas mixtures following JCGM (2008) as 0.062 pmol mol−1 (k = 2) for molar fractions in the range 7–10 pmol mol−1, equivalent to 0.6 to 0.9 % (Lim et al., 2017). The preparation method followed by NOAA has similarities with to the one developed at SIO, being based on static gravimetry as well. We therefore conservatively use an uncertainty for the SIO-05 assigned molar fraction of 1 % (k = 2).
METAS-2016. A set of two primary standards was prepared at METAS in 2016 (METAS-2016 calibration scale) to participate in an intercomparison for SF6 in air at atmospheric molar fractions organised by the World Calibration Centre for SF6 (Lee et al., 2017). These two standards were prepared using a similar method as for the METAS-2017 scale using permeation, dynamic dilution and cryo-filling of a nmol mol−1 molar fraction mixture containing SF6 only in synthetic air (Lee et al., 2017). After homogenisation this mother mixture was dynamically diluted into two daughter mixtures at 8 and 10 pmol mol−1, themselves transferred in cylinders by cryo-filling. The expanded uncertainty of the prepared standards is U = 1.3 %.
NOAA-2014. We use for comparison the known SF6 conversion factor between the NOAA-2014 and SIO-05 calibration scales of NOAA-2014 / SIO-05 = 1.002 ± 0.002 (Krummel et al., 2017), based on measurements at co-located stations and tank exchanges. This conversion factor is depicted by a green dashed line in Fig. 6.
METAS-2015. A first HFC-125 reference gas mixture was produced by METAS in 2015, with expanded uncertainty U = 2 %. The preparation method consisted of a permeation step using MSB-Violetta, generating a mixture at approximately 85 nmol mol−1, transferred to a SilcoNert-coated cylinder by cryo-filling. This mother mixture was then diluted into two daughter mixtures both at 17 pmol mol−1, using dynamic dilution steps. Diluted daughter mixtures were then directly injected into the Medusa GC-MS (Supplement Sect. S4 and Figs. S1, S2 and S3).
NOAA-2008. We use the published conversion factor NOAA-2008 / SIO-14 of 0.946 ± 0.008 (Simmonds et al., 2017). This factor was determined by comparing AGAGE continuous measurements by Medusa GC-MS and NOAA flask samples taken from three co-located sites.
UB-98. Before using the SIO-14 scale for HFC-125 within AGAGE, a primary calibration scale prepared by University of Bristol was in use (O'Doherty et al., 2004, 2009). The known conversion factor UB-98 / SIO-14 is 0.9237 (Christina M. Harth and Ray F. Weiss, personal communication, 2018).
Empa-2013. In 2013 Empa prepared a first calibration scale for a set of newly emitted compounds, including HFO-1234yf at ≈ 2 pmol mol−1, using volumetric dilution (Vollmer et al., 2015). The molar fraction uncertainty is likely no more than U ≤ 30 %.
METAS-2015. Two reference gas mixtures for HFO-1234yf were produced in 2015 at 2 pmol mol−1 following the same method as for the METAS-2015 standards for HFC-125 (Supplement Sect. S4). The resulting expanded uncertainty was U = 2.5 %.
Interim-98. For CFC-13, a preliminary calibration scale was developed at the University of Bristol based on dilution of a high molar fraction reference gas mixture purchased from a commercial manufacturer (O'Doherty et al., 2004). This preliminary scale has been used within AGAGE until 2017. It would be very difficult to assign an uncertainty to this mixture, and it should be noted that the aim of this Interim-98 standard was to serve as intermediate anchor in order to be able to report CFC-13 internally within AGAGE. We tentatively assign an uncertainty of 5 %. The Interim-98 scale was transferred to the SIO suite of secondary standards by measurement comparisons, and therefore the tertiary standards used at Empa have a CFC-13 assigned value on the Interim-98 scale. We provide comparison to this calibration scale to contribute to the documentation of the scale transfer from Interim-98 to METAS-2017 for CFC-13 rather than to realise a new intercomparison. Additional details on the scale transfer are given in Vollmer et al. (2018).
4.2 Results of comparisons
We express the results as ratio values, i.e. the value expressed on the METAS calibration scale divided by the value expressed on the historical calibration scale (SIO-05, SIO-14, Empa-2013 and Interim-98). A ratio of 1 would correspond to a perfect agreement between the two compared scales.
The average ratio METAS-2017 / SIO-05 using all available results is 1.002, i.e. the deviation from 1 is clearly within the ratio uncertainty (Fig. 6). This demonstrates that the two calibration scales are concordant with each other. SF6 is also the substance for which the NOAA / SIO ratio is closest to 1 (NOAA-2014 / SIO-05 ratio of 1.002). Combining these two ratios, one obtains a METAS-2017 / NOAA-2014 ratio for SF6 of 1.000. Such excellent agreements, in particular between standards produced by dynamic and static methods, can in addition to the reliability of both preparation methods be explained by the stability and non-reactivity (for instance low adsorptivity) of this substance.
The METAS-2017 calibration scale for HFC-125 is 7 % lower than the SIO-14 scale (METAS-2017 / SIO-14 = 0.929). For HFC-125, the value assigned on the SIO-14 scale is not corrected for potential impurities in the pure HFC-125 substance used for the preparation (Christina M. Harth, personal communication, 2017), while a 1 % correction is used for the METAS-2017 scale. Assuming HFC-125 sources were similar, and both scales would apply the same procedure, the disagreement would be reduced to 6 %. We plan for future reference gas mixture preparation to check the presence of substance impurities in permeators in a systematic way, to get a better estimate of the purity fraction and to quantify any potential cross-contamination, if any. For the METAS-2017 scale, we checked in particular the absence of HFC-132b as impurity in the HFC-125 permeator (see Sect. S5 in the Supplement).
Comparison of the METAS-2015 and SIO-14 calibration scales also showed a METAS value lower than the SIO value, by 4 % (METAS-2015 / SIO-14 = 0.963). The ratio NOAA-2008 / SIO-14 for HFC-125 is 0.946, one of the largest discrepancies observed between SIO and NOAA for halogenated gases (: Paul B. Krummel and Bradley D. Hall, personal communication, January 2018). Thus, comparing the SIO-14 calibration scale for HFC-125 to these three other scales (NOAA-08, METAS-2015, METAS-2017) points to a probable overestimation of the SIO-14 value (due e.g. to substance losses by adsorption) or that the values on the other calibration scales are underestimated (due to e.g. unaccounted for contamination). Cases of gas standards produced by dynamic methods yielding results lower than those produced by static method have already been observed several times with reactive substances prone to adsorption on surfaces, such as ammonia on stainless steel (van der Veen et al., 2010). This is due to the fact that when applying dynamic methods, potential losses by adsorption on surfaces can be cancelled out when the generation process reaches equilibrium, after a sufficiently long stabilisation time. Interestingly, the METAS-2017 calibration scale is even lower than the METAS-2015 scale, by 3 %. Significant improvements in the generation process were made for the METAS-2017 scale to considerably minimise the total exposition to metal surfaces, compared to the METAS-2015 scale. This improvement is potentially the cause of the observed 3 % shift towards lower values.
For HFO-1234yf, the two METAS calibration scales are lower than the Empa-2013 scale with, on average, METAS-2017 / Empa-2013 = 0.910 and METAS-2015 / Empa-2013 = 0.971. METAS-2017 is thus lower than METAS-2015 by 6 %. As with HFC-125, this latter ratio is significantly lower than 1. Within all halogenated substances studied here, this is the largest discrepancy observed between static (Empa-2013) and dynamic (METAS-2017) preparation methods, as well as the largest offset between different METAS scales. In addition, within the METAS-2017 suite of cylinders the largest offsets for outliers are also observed for HFO-1234yf (Fig. 5, HFO-1234yf in MP-010 and MP-008 is ≈ 5.5 % too low). All these observations suggest that preparing primary reference gas mixtures for HFO-1234yf with U ≤ 2 % may require dynamic generation methods, with additional minimisation of contact with surfaces.
Three cylinders on the Interim-98 calibration scale for CFC-13 have been compared to cylinder MP-001 on the METAS-2017 scale. For documentation purposes, we report the average ratio METAS-2017 / Interim-98 measured as 1.055 (Fig. 6). The comparison was extended to additional cylinders in Vollmer et al. (2018) to ensure a reliable scale transfer from Interim-98 to METAS-2017 within the AGAGE network.
We have developed a suite of primary, SI-traceable reference gas mixtures in 11 pressurised cylinders for SF6, HFC-125, HFO-1234yf, HCFC-132b and CFC-13 in synthetic air, at atmospheric molar fractions. This suite constitutes the METAS-2017 primary calibration scales for these five halogenated compounds. This work therefore combines the advantages of SI-traceable reference gas mixture preparation with a primary calibration scale system for its use as anchor by a monitoring network. Such a combined system allows us to maximise the compatibility (as defined by GAW) within the network while linking all reference values to the International System of Units and assigning carefully estimated uncertainties following international guidelines (JCGM 100:2008).
Expanded uncertainties of the METAS-2017 calibration scale after verification ranges from 1 to 2 % at a 95 % confidence interval. Such molar fractions at the pmol mol−1 level with associated expanded uncertainties of no more than 2 % clearly mark a step beyond the state of the art for dynamic methods. We have demonstrated the applicability of dynamic gravimetric generation methods coupled to cryo-filling in cylinders to prepare primary reference gas mixtures for halogenated compounds as low as 1 pmol mol−1. For stable compounds for which static gravimetric methods are also applicable (e.g. SF6), these latter methods perform better in terms of expanded uncertainties (Lim et al., 2017), but we emphasise that using a completely independent preparation method may always help to detect potential systematic biases affecting one method or the other. From a metrological point of view, this preparation exercise is therefore highly valuable, ensuring comparability and redundancy of prepared reference gas mixtures.
Comparison of the METAS-2017 calibration scale for SF6 with the scale prepared by SIO (SIO-05) leads to a conversion factor METAS-2017 / SIO-05 of 1.002, illustrating the concordance of the two scales within uncertainties. An indirect comparison with the NOAA calibration scale also yields agreeing results (METAS-2017 / NOAA-2014 = 1.000). The excellent concordance obtained for SF6 gives confidence in the reliability of the presented dynamic–gravimetric method to prepare standards for other, more reactive compounds, e.g. HFC-125.
For HFC-125, known to be more reactive than SF6, the METAS-2017 calibration scale is measured as 7 % lower than SIO-14. In addition the METAS-2017 scale for HFO-1234yf is measured 9 % lower than Empa-2013. Such an offset towards lower values for standards prepared using dynamic generation methods by contrast to methods using static gravimetry or static volumetry has been previously observed for other reactive compounds such as ammonia. This underlines the risk of substance losses by, for example, adsorption on surfaces for HFC-125 and HFO-1234yf (and potentially other reactive substances). Dynamic generation methods and/or minimisation of contact on surfaces should therefore be favoured when preparing primary reference gas mixtures for such reactive substances.
All data used to prepare the METAS-2017 suite and results of cylinder measurements within the suite are available in this article and its Supplement.
The supplement related to this article is available online at: https://doi.org/10.5194/amt-11-3351-2018-supplement.
The authors declare that they have no conflict of interest.
This article is part of the special issue “The 10th International Carbon Dioxide Conference (ICDC10) and the 19th WMO/IAEA Meeting on Carbon Dioxide, other Greenhouse Gases and Related Measurement Techniques (GGMT-2017) (AMT/ACP/BG/CP/ESD inter-journal SI)”. It is a result of the 19th WMO/IAEA Meeting on Carbon Dioxide, Other Greenhouse Gases, and Related Measurement Techniques (GGMT-2017), Empa Dübendorf, Switzerland, 27–31 August 2017.
Myriam Guillevic would like to thank members of the AGAGE community, and in
particular Ray Weiss, for their interest in this work, which was a strong
motivation to bring this work to the presented scientific level. Many thanks
to the colleagues from METAS Gas Analysis laboratory, including Hans-Peter Haerri, now retired, and from other
laboratories and technical departments, particularly the METAS workshop, who
supported this work with their skills and their equipment. We thank
collaborators of the HIGHGAS project for their valuable inputs: Fine
Metrology Srls, SilcoTek GmbH, SAES Getters, Arbor Fluidtec AG and VICI
Metronics AG. This project has been funded by METAS project AtmoChemECV and
the EMRP project HIGHGAS. The EMRP is jointly funded by the EMRP
participating countries within EURAMET and the European Union. Free softwares
Inkscape and Octave were used for data analysis, diagrams and
plots. Finally, we thank two anonymous reviewers for their time and excellent, constructive comments.
Edited by: Markus Leuenberger
Reviewed by: two anonymous referees
Arnold, T., Mühle, J., Salameh, P. K., Harth, C. M., Ivy, D. J., and Weiss, R. F.: Automated measurement of nitrogen trifluoride in ambient air, Anal. Chem., 84, 4798–4804, https://doi.org/10.1021/ac300373e, 2012. a
Brunner, D., Arnold, T., Henne, S., Manning, A., Thompson, R. L., Maione, M., O'Doherty, S., and Reimann, S.: Comparison of four inverse modelling systems applied to the estimation of HFC-125, HFC-134a, and SF6 emissions over Europe, Atmos. Chem. Phys., 17, 10651–10674, https://doi.org/10.5194/acp-17-10651-2017, 2017. a
Cantrell, C. A.: Technical Note: Review of methods for linear least-squares fitting of data and application to atmospheric chemistry problems, Atmos. Chem. Phys., 8, 5477–5487, https://doi.org/10.5194/acp-8-5477-2008, 2008. a
Craig, H.: Isotopic standards for carbon and oxygen and correction factors for mass-spectrometric analysis of carbon dioxide, Geochim. Cosmochim. Ac., 12, 133–149, 1957. a
Dansgaard, W.: Comparative measurements of standards for carbon isotopes, Geochim. Cosmochim. Ac., 3, 253–256, 1953. a
Dlugokencky, E. J., Myers, R. C., Lang, P. M., Masarie, K. A., Crotwell, A. M., Thoning, K. W., Hall, B. D., Elkins, J. W., and Steele, L. P.: Conversion of NOAA atmospheric dry air CH4 mole fractions to a gravimetrically prepared standard scale, J. Geophys. Res.-Atmos., 110, D18306, https://doi.org/10.1029/2005JD006035, 2005. a
Flores, E., Viallon, J., Moussay, P., Idrees, F., and Wielgosz, R. I.: Highly accurate nitrogen dioxide (NO2) in nitrogen standards based on permeation, Anal. Chem., 84, 10283–10290, https://doi.org/10.1021/ac3024153, 2012. a
Fuchs, P., Marti, K., and Russi, S.: New instrument for the study of “the kg, mise en pratique”: first results on the correlation between the change in mass and surface chemical state, Metrologia, 49, 607–614, https://doi.org/10.1088/0026-1394/49/6/607, 2012. a
Haerri, H.-P., Macé, T., Waldén, J., Pascale, C., Niederhauser, B., Wirtz, K., Stovcik, V., Sutour, C., Couette, J., and Waldén, T.: Dilution and permeation standards for the generation of NO, NO2 and SO2 calibration gas mixtures, Meas. Sci. Technol., 28, 035801, https://doi.org/10.1088/1361-6501/aa543d, 2017. a, b
Hall, B. D., Dutton, G. S., and Elkins, J. W.: The NOAA nitrous oxide standard scale for atmospheric observations, J. Geophys. Res.-Atmos., 112, D09305, https://doi.org/10.1029/2006JD007954, 2007. a, b, c
Hall, B. D., Engel, A., Mühle, J., Elkins, J. W., Artuso, F., Atlas, E., Aydin, M., Blake, D., Brunke, E.-G., Chiavarini, S., Fraser, P. J., Happell, J., Krummel, P. B., Levin, I., Loewenstein, M., Maione, M., Montzka, S. A., O'Doherty, S., Reimann, S., Rhoderick, G., Saltzman, E. S., Scheel, H. E., Steele, L. P., Vollmer, M. K., Weiss, R. F., Worthy, D., and Yokouchi, Y.: Results from the International Halocarbons in Air Comparison Experiment (IHALACE), Atmos. Meas. Tech., 7, 469–490, https://doi.org/10.5194/amt-7-469-2014, 2014. a, b
Harris, N. R. and Wuebbles, D. J.: Scientific assessment of ozone depletion: 2014, chap. 5: Scenarios and information for policymakers, World Meteorological Organization (WMO), Global Ozone Research and Monitoring Project, Geneva, Switzerland, 5.1–5.58, available at: https://www.wmo.int/pages/prog/arep/gaw/ozone_2014/documents/Full_report_2014_Ozone_Assessment.pdf (last access: 16 January 2018), 2014. a
Hu, L., Montzka, S. A., Miller, B. R., Andrews, A. E., Miller, J. B., Lehman, S. J., Sweeney, C., Miller, S. M., Thoning, K., Siso, C., Atlas, E. L., Blake, D. R., de Gouw, J., Gilman, J. B., Dutton, G., Elkins, J. W., Hall, B., Chen, H., Fischer, M. L., Mountain, M. E., Nehrkorn, T., Biraud, S. C., Moore, F. L., and Tans, P.: Continued emissions of carbon tetrachloride from the United States nearly two decades after its phaseout for dispersive uses, P. Natl. Acad. Sci. USA, 113, 2880–2885, https://doi.org/10.1073/pnas.1522284113, 2016. a
ISO 6142-1: Gas analysis – Preparation of calibration gas mixtures – Part 1: Gravimetric method for Class I mixtures, available at: https://www.iso.org/standard/59631.html (last access: 16 January 2018), 2015. a
ISO 6143: Gas analysis – Comparison methods for determining and checking the composition of calibration gas mixtures, available at: https://www.iso.org/standard/24665.html (last access: 16 January 2018), 2001. a, b
ISO 6145-10: Gas analysis – Preparation of calibration gas mixtures using dynamic volumetric methods – Part 10: Permeation method, available at: https://www.iso.org/standard/25916.html (last access: 16 January 2018), 2002. a, b
ISO 6145-7: Gas analysis – Preparation of calibration gas mixtures using dynamic volumetric methods – Part 7: Thermal mass-flow controllers, available at: https://www.iso.org/standard/45471.html (last access: 16 January 2018), 2009. a, b, c
JCGM: JCGM 100:2008: Evaluation of measurement data – Guide to the expression of uncertainty in measurement, available at: https://www.bipm.org/utils/common/documents/jcgm/JCGM_100_2008_E.pdf (last access: 8 June 2018), 2008. a
Krummel, P., Montzka, S., Harth, C., Miller, B., Mühle, J., Dlugokencky, E., Salameh, P., Hall, B., O'Doherty, S., Steele, L., Dutton, G., Young, D., Nance, J., Langenfelds, R., Elkins, J., Loh, Z., Lunder, C., Fraser, P., Derek, N., Mitrevski, B., Weiss, R., and Prinn, R.: An update of comparisons of non-CO2 trace gas measurements between AGAGE and NOAA at common sites, in: 19th WMO/IAEA Meeting on Carbon Dioxide, Other Greenhouse Gases, and Related Measurement Techniques (GGMT-2017), Dübendorf, Switzerland, available at: http://www.wmo.int/pages/prog/arep/gaw/documents/GGMT2017_T03_Krummel.pdf (last access: 16 January 2018), 2017. a, b
Laube, J. C., Newland, M. J., Hogan, C., Brenninkmeijer, C. A. M., Fraser, P. J., Martinerie, P., Oram, D. E., Reeves, C. E., Röckmann, T., Schwander, J., Witrant, E., and Sturges, W. T.: Newly detected ozone-depleting substances in the atmosphere, Nat. Geosci., 7, 266–269, https://doi.org/10.1038/NGEO2109, 2014. a
Lee, H., Kim, S., woo Cheo, H., Lee, S.-B. R. H., Kim, S., Yoo, H., Choe, H., Han, S.-O., Ryoo, S.-B., Hall, B., Forster, G., Camille, Y. K., Meinhardt, F., Levin, I., Ries, L., Couret, C., Guillevic, M., Wyss, S. A., Rzesanke, D., Jordan, A., Moss, R., Lee, J., and Yu, J.: The report of the result on SF6 Inter-Comparison Experiment, 2016–2017, Tech. rep., World Meteorological Organization, available at: http://www.wmo.int/pages/prog/arep/gaw/documents/FinalreportofSICE_2016.pdf (last access: 16 January 2018), 2017. a, b
Lim, J. S., Lee, J., Moon, D., Kim, J. S., Lee, J., and Hall, B. D.: Gravimetric standard gas mixtures for global monitoring of atmospheric SF6, Anal. Chem., 89, 12068–12075, https://doi.org/10.1021/acs.analchem.7b02545, 2017. a, b, c, d
Lunt, M. F., Rigby, M., Ganesan, A. L., Manning, A. J., Prinn, R. G., O'Doherty, S., Mühle, J., Harth, C. M., Salameh, P. K., Arnold, T., Weiss, R. F., Saito, T., Yokouchi, Y., Krummel, P. B., Steele, L. P., Fraser, P. J., Li, S., Park, S., Reimann, S., Vollmer, M. K., Lunder, C., Hermansen, O., Schmidbauer, N., Maione, M., Arduini, J., Young, D., and Simmonds, P. G.: Reconciling reported and unreported HFC emissions with atmospheric observations, P. Natl. Acad. Sci. USA, 112, 5927–5931, https://doi.org/10.1073/pnas.1420247112, 2015. a
Marti, K.: Dissemination of the Kilogramm at METAS: Extended Lagrange multiplier method, in: IMEKO 2017, IMEKO, available at: https://www.metas.ch/dam/data/metas/ForschungEntwicklung/publikationen/IMEKO_Poster_Marti.pdf (last access: 16 January 2018), 2017. a
Meija, J., Coplen, T. B., Berglund, M., Brand, W. A., Bievre, P. D., Gröning, M., Holden, N. E., Irrgeher, J., Loss, R. D., Walczyk, T., and Prohaska, T.: Atomic weights of the elements 2013 (IUPAC Technical Report), Pure Appl. Chem., 88, 265–291, https://doi.org/10.1515/pac-2015-0305, 2016. a
Miller, B. R., Weiss, R. F., Salameh, P. K., Tanhua, T., Greally, B. R., Mühle, J., and Simmonds, P. G.: Medusa: A sample preconcentration and GC/MS detector system for in situ measurements of atmospheric trace halocarbons, hydrocarbons, and sulfur compounds, Anal. Chem., 80, 1536–1545, https://doi.org/10.1021/ac702084k, 2008. a, b, c
Molina, M. J. and Rowland, F. S.: Stratospheric sink for chlorofluoromethanes: Chlorine atom catalyzed destruction of ozone, Nature, 249, 810–812, 1974. a
Montzka, S. A., McFarland, M., Andersen, S. O., Miller, B. R., Fahey, D. W., Hall, B. D., Hu, L., Siso, C., and Elkins, J. W.: Recent trends in global emissions of hydrochlorofluorocarbons and hydrofluorocarbons: Reflecting on the 2007 adjustments to the Montreal Protocol, J. Phys. Chem. A, 119, 4439–4449, https://doi.org/10.1021/jp5097376, 2015. a
Moosbach, E. and Hartkamp, H.: Influence of carrier gas pressure on permeation rates of the “Hartmann & Braun CGP”-NO2-test gas generator, Fresenius J. Anal. Chem., 347, 475–477, https://doi.org/10.1007/BF00324236, 1993. a
Niederhauser, B. and Barbe, J.: Bilateral comparison of primary low-gas-flow standards between the BNM-LNE and METAS, Metrologia, 39, 573–578, 2002. a
O'Doherty, S., Cunnold, D. M., Manning, A., Miller, B. R., Wang, R. H. J., Krummel, P. B., Fraser, P. J., Simmonds, P. G., McCulloch, A., Weiss, R. F., Salameh, P., Porter, L. W., Prinn, R. G., Huang, J., Sturrock, G., Ryall, D., Derwent, R. G., and Montzka, S. A.: Rapid growth of hydrofluorocarbon 134a and hydrochlorofluorocarbons 141b, 142b, and 22 from Advanced Global Atmospheric Gases Experiment (AGAGE) observations at Cape Grim, Tasmania, and Mace Head, Ireland, J. Geophys. Res.-Atmos., 109, D06310, https://doi.org/10.1029/2003JD004277, 2004. a, b, c
O'Doherty, S., Cunnold, D. M., Miller, B. R., Mühle, J., McCulloch, A., Simmonds, P. G., Manning, A. J., Reimann, S., Vollmer, M. K., Greally, B. R., Prinn, R. G., Fraser, P. J., Steele, L. P., Krummel, P. B., Dunse, B. L., Porter, L. W., Lunder, C. R., Schmidbauer, N., Hermansen, O., Salameh, P. K., Harth, C. M., Wang, R. H. J., and Weiss, R. F.: Global and regional emissions of HFC-125 (CHF2CF3) from in situ and air archive atmospheric observations at AGAGE and SOGE observatories, J. Geophys. Res.-Atmos., 114, D23304, https://doi.org/10.1029/2009JD012184, 2009. a
Pogàny, A., Balslev-Harder, D., Braban, C. F., Cassidy, N., Ebert, V., Ferracci, V., Hieta, T., Leuenberger, D., Martin, N. A., Pascale, C., Peltola, J., Persijn, S., Tiebe, C., Twigg, M. M., Vaittinen, O., van Wijk, J., Wirtz, K., and Niederhauser, B.: A metrological approach to improve accuracy and reliability of ammonia measurements in ambient air, Meas. Sci. Technol., 27, 115012, https://doi.org/10.1088/0957-0233/27/11/115012, 2016. a
Prinn, R. G., Weiss, R. F., Fraser, P. J., Simmonds, P. G., Cunnold, D. M., Alyea, F. N., O'Doherty, S., Salameh, P., Miller, B. R., Huang, J., Wang, R. H. J., Hartley, D. E., Harth, C., Steele, L. P., Sturrock, G., Midgley, P. M., and McCulloch, A.: A history of chemically and radiatively important gases in air deduced from ALE/GAGE/AGAGE, J. Geophys. Res.-Atmos., 105, 17751–17792, https://doi.org/10.1029/2000JD900141, 2000. a, b, c, d, e, f, g, h
Rhoderick, G. C., Hall, B. D., Harth, C. M., Kim, J. S., Lee, J., Montzka, S. A., Mühle, J., Reimann, S., Vollmer, M. K., and Weiss, R. F.: Comparison of halocarbon measurements in an atmospheric dry whole air sample, Elementa, 3, 000075, https://doi.org/10.12952/journal.elementa.000075, 2015. a, b, c
Rigby, M., Mühle, J., Miller, B. R., Prinn, R. G., Krummel, P. B., Steele, L. P., Fraser, P. J., Salameh, P. K., Harth, C. M., Weiss, R. F., Greally, B. R., O'Doherty, S., Simmonds, P. G., Vollmer, M. K., Reimann, S., Kim, J., Kim, K.-R., Wang, H. J., Olivier, J. G. J., Dlugokencky, E. J., Dutton, G. S., Hall, B. D., and Elkins, J. W.: History of atmospheric SF6 from 1973 to 2008, Atmos. Chem. Phys., 10, 10305–10320, https://doi.org/10.5194/acp-10-10305-2010, 2010. a, b, c
Scaringelli, F. P., O'Keeffe, A. E., Rosenberg, E., and Bell, J. P.: Preparation of known concentrations of gases and vapors with permeation devices calibrated gravimetrically, Anal. Chem., 42, 871–876, https://doi.org/10.1021/ac60290a012, 1970. a
Sherry, D., McCulloch, A., Liang, Q., Reimann, S., and Newman, P. A.: Current sources of carbon tetrachloride (CCl4) in our atmosphere, Environ. Res. Lett., 13, 024004, https://doi.org/10.1088/1748-9326/aa9c87, 2017. a
Simmonds, P. G., Rigby, M., Manning, A. J., Lunt, M. F., O'Doherty, S., McCulloch, A., Fraser, P. J., Henne, S., Vollmer, M. K., Mühle, J., Weiss, R. F., Salameh, P. K., Young, D., Reimann, S., Wenger, A., Arnold, T., Harth, C. M., Krummel, P. B., Steele, L. P., Dunse, B. L., Miller, B. R., Lunder, C. R., Hermansen, O., Schmidbauer, N., Saito, T., Yokouchi, Y., Park, S., Li, S., Yao, B., Zhou, L. X., Arduini, J., Maione, M., Wang, R. H. J., Ivy, D., and Prinn, R. G.: Global and regional emissions estimates of 1,1-difluoroethane (HFC-152a, CH3CHF2) from in situ and air archive observations, Atmos. Chem. Phys., 16, 365–382, https://doi.org/10.5194/acp-16-365-2016, 2016. a
Simmonds, P. G., Rigby, M., McCulloch, A., O'Doherty, S., Young, D., Mühle, J., Krummel, P. B., Steele, P., Fraser, P. J., Manning, A. J., Weiss, R. F., Salameh, P. K., Harth, C. M., Wang, R. H. J., and Prinn, R. G.: Changing trends and emissions of hydrochlorofluorocarbons (HCFCs) and their hydrofluorocarbon (HFCs) replacements, Atmos. Chem. Phys., 17, 4641–4655, https://doi.org/10.5194/acp-17-4641-2017, 2017. a, b, c
Vaittinen, O., Metsälä, M., Persijn, S., Vainio, M., and Halonen, L.: Adsorption of ammonia on treated stainless steel and polymer surfaces, Appl. Phys. B, 115, 1–12, https://doi.org/10.1007/s00340-013-5590-3, 2013. a
van der Veen, A. M. H., Nieuwenkamp, G., Wessel, R. M., Maruyama, M., Heo, G. S., doo Kim, Y., Moon, D. M., Niederhauser, B., Quintilii, M., Milton, M. J. T., Cox, M. G., Harris, P. M., Guenther, F. R., Rhoderick, G. C., Konopelko, L. A., Kustikov, Y. A., Pankratov, V. V., Selukov, D. N., Petrov, V. A., and Gromova, E. V.: International comparison CCQM-K46: Ammonia in nitrogen, Tech. rep., avialable at: http://stacks.iop.org/0026-1394/47/i=1A/a=08023 (last access: 16 January 2018), 2010. a
Velders, G. J. M., Fahey, D. W., Daniel, J. S., McFarland, M., and Andersen, S. O.: The large contribution of projected HFC emissions to future climate forcing, P. Natl. Acad. Sci. USA, 106, 10949–10954, https://doi.org/10.1073/pnas.0902817106, 2009. a
Vollmer, M. K., Reimann, S., Hill, M., and Brunner, D.: First observations of the fourth generation synthetic halocarbons HFC-1234yf, HFC-1234ze(E), and HCFC-1233zd(E) in the atmosphere, Environ. Sci. Technol., 49, 2703–2708, https://doi.org/10.1021/es505123x, 2015. a, b, c, d
Vollmer, M. K., Mühle, J., Trudinger, C. M., Rigby, M., Montzka, S. A., Harth, C. M., Miller, B. R., Henne, S., Krummel, P. B., Hall, B. D., Young, D., Kim, J., Arduini, J., Wenger, A., Yao, B., Reimann, S., O'Doherty, S., Maione, M., Etheridge, D. M., Li, S., Verdonik, D. P., Park, S., Dutton, G., Steele, L. P., Lunder, C. R., Rhee, T. S., Hermansen, O., Schmidbauer, N., Wang, R. H. J., Hill, M., Salameh, P. K., Langenfelds, R. L., Zhou, L., Blunier, T., Schwander, J., Elkins, J. W., Butler, J. H., Simmonds, P. G., Weiss, R. F., Prinn, R. G., and Fraser, P. J.: Atmospheric histories and global emissions of halons H-1211 (CBrClF2), H-1301 (CBrF3), and H-2402 (CBrF2CBrF2), J. Geophys. Res.-Atmos., 121, 3663–3686, https://doi.org/10.1002/2015JD024488, 2016. a
Vollmer, M. K., Young, D., Trudinger, C. M., Mühle, J., Henne, S., Rigby, M., Park, S., Li, S., Guillevic, M., Mitrevski, B., Harth, C. M., Miller, B. R., Reimann, S., Yao, B., Steele, L. P., Wyss, S. A., Lunder, C. R., Arduini, J., McCulloch, A., Wu, S., Rhee, T. S., Wang, R. H. J., Salameh, P. K., Hermansen, O., Hill, M., Langenfelds, R. L., Ivy, D., O'Doherty, S., Krummel, P. B., Maione, M., Etheridge, D. M., Zhou, L., Fraser, P. J., Prinn, R. G., Weiss, R. F., and Simmonds, P. G.: Atmospheric histories and emissions of chlorofluorocarbons CFC-13 (CClF3), SCFC-114 (C2Cl2F4), and CFC-115 (C2ClF5), Atmos. Chem. Phys., 18, 979–1002, https://doi.org/10.5194/acp-18-979-2018, 2018. a, b, c, d
Weiss, R. F. and Prinn, R. G.: Quantifying greenhouse-gas emissions from atmospheric measurements: a critical reality check for climate legislation, Philos. T. R. Soc. A, 369, 1925–1942, https://doi.org/10.1098/rsta.2011.0006, 2011. a
WMO: The stratosphere 1981: Theory and measurements, Tech. rep., World Meteorological Organization, Geneva, Switzerland, 1981. a
WMO: Scientific assessment of ozone depletion: 2014, Report 55, World Meteorological Organization (WMO), Global Ozone Research and Monitoring Project, Geneva, Switzerland, available at: https://www.wmo.int/pages/prog/arep/gaw/ozone_2014/documents/Full_report_2014_Ozone_Assessment.pdf (last access: 16 January 2018), 2014. a, b
Yokohata, A., Makide, Y., and Tominaga, T.: A new calibration method for the measurement of CCl4 concentration at 10−10 v/v level and the behavior of CCl4 in the atmosphere, B. Chem. Soc. Jpn., 58, 1308–1314, https://doi.org/10.1246/bcsj.58.1308, 1985. a
York, D., Evensen, N. M., Martínez, M. L., and Delgado, J. D. B.: Unified equations for the slope, intercept, and standard errors of the best straight line, Am. J. Phys., 72, 367–375, https://doi.org/10.1119/1.1632486, 2004. a
Zhao, C. L., Tans, P. P., and Thoning, K. W.: A high precision manometric system for absolute calibrations of CO2 in dry air, J. Geophys. Res.-Atmos., 102, 5885–5894, https://doi.org/10.1029/96JD03764, 1997. a