In Force

BIPM Rapport

BIPM BIPM-2018/03 : 2018
Comparison of ozone reference standards of the CHMI and the BIPM
Joële Viallon (BIPM), Philippe Moussay (BIPM), Faraz Idrees (BIPM), Robert Wielgosz (BIPM), Jan Šilhavý (CHMI), Miroslav Vokoun (CHMI).
BIPM Rapport

In Force




Foreword

A comparison of the ozone reference standards of the Czech Hydrometeorological Institute (CHMI) and of the Bureau International des Poids et Mesures (BIPM) has been performed. Both institutes maintain Standard Reference Photometers (SRPs), developed by the National Institute of Standards and Technology (NIST), as their reference standards. The instruments were compared over a nominal ozone mole fraction range of 0nmol/mol to 860nmol/mol and the results showed good agreement.

Comparison of ozone reference standards of the CHMI and the BIPM

1.  Introduction

A comparison of the ozone reference standards of the Czech Hydrometeorological Institute (CHMI) and of the Bureau International des Poids et Mesures (BIPM) has been performed. Both institutes maintain Standard Reference Photometers (SRPs) developed by the National Institute of Standards and Technology (NIST) as their reference standards. This comparison was performed following the protocol established for the key comparison BIPM.QM-K1, adapted to the measurement range normally covered by the CHMI as part of their calibration service, which is 10nmol/mol to 870nmol/mol . The measurement protocol is described briefly in Chapter 4. A description of the standards is given in Chapter 5 of this report, together with their uncertainty budgets. The data treatment is explained in Chapter 6, and the results of the comparison are given in Chapter 7.

2.  Terms and definitions

No terms and definitions are listed in this document.

x nom

nominal ozone mole fraction in dry air supplied by the ozone generator.

x A , i

i th measurement of the nominal value x nom by the photometer A.

x ¯ A

the mean of N measurements of the nominal value x nom measured by the photometer A : x ¯ A = 1 N i = 1 N x A , i

s A

standard deviation of N measurements of the nominal value x nom measured by the photometer A : s A 2 = 1 N 1 i = 1 N x A , i x ¯ A 2

  • The result of the linear regression fit performed between two sets of data measured by photometers A and B during a comparison is written as: x A = a A,B x B + b A,B . With this notation, photometer A is compared to photometer B. a A,B is dimensionless and b A,B is expressed in units of nmol/mol .

3.  Measurement schedule

The measurements were performed on 20 September 2017 at the BIPM.

4.  Measurement protocol

This comparison was performed following the protocol established for the key comparison BIPM.QM-K1, with a modified measurement range: the instruments were compared over a nominal ozone mole fraction range of 0nmol/mol to 860nmol/mol (instead of a nominal range of 0nmol/mol to 500nmol/mol ).

The comparison protocol is summarized in this section. The complete version can be downloaded from the BIPM website (http://www.bipm.org/utils/en/pdf/BIPM.QM-K1_protocol.pdf).

This comparison was performed following protocol A, corresponding to a direct comparison between the CHMI national standard SRP17 and the common reference standard BIPM-SRP27 maintained at the BIPM. A comparison between two (or more) ozone photometers consists of producing ozone-air mixtures at different amount-of-substance fractions over the required range, and measuring these with the photometers.

4.1.  Ozone generation

The same source of purified air is used for all the ozone photometers being compared. This air is used to provide reference air as well as the ozone–air mixture to each ozone photometer. Ambient air is used as the source for reference air. The air is compressed with an oil-free compressor. It is dried and scrubbed using a commercial purification system so that the mole fraction of ozone and nitrogen oxides remaining in the air is below detectable limits. The relative humidity of the reference air is monitored and the mole fraction of water in air is typically less than 3µmol/mol . The mole fraction of volatile organic hydrocarbons in the reference air was measured (November 2002), with no mole fraction of any detected component exceeding 1nmol/mol .

A common dual external manifold in Pyrex is used to supply the necessary flows of reference air and ozone–air mixtures to the ozone photometers. The two columns of this manifold are vented to atmospheric pressure.

4.2.  Comparison procedure

All the instruments were switched on and allowed to stabilize for at least 8 hours before the comparison commenced. The pressure and temperature measurement systems of the instruments were checked at this time. If any adjustments were required, these were noted. For this comparison, no adjustments were necessary.

One comparison run normally includes ten different amount-of-substance fractions distributed to cover the range, together with the measurement of reference air at the beginning and end of each run. The nominal amount-of-substance fractions were measured in a sequence imposed by the BIPM.QM-K1 protocol, plus three additional values to extend the range (0, 220, 80, 420, 120, 320, 30, 370, 170, 500, 270, 0, 570, 840, 700 and 0) nmol/mol . Each of these points is an average of ten single measurements.

For each nominal value of the ozone mole fraction x nom supplied by the ozone generator, the standard deviation s SRP27 on the set of ten consecutive measurements x SRP27 , i recorded by BIPM-SRP27 was calculated. The measurement results were considered as valid if s SRP27 was less than 1nmol/mol , which ensures that the photometers were measuring a stable ozone concentration. If not, another series of ten consecutive measurements was performed.

4.3.  Comparison repeatability

The comparison procedure was repeated continuously for BIPM.QM-K1 to evaluate its repeatability. The participant and the BIPM jointly decided when both instruments were stable enough to start recording a set of measurement results to be considered as the official comparison results. The measurements reported here were recorded in an additional run that was performed immediately afterwards.

4.4.  SRP27 stability check

A second ozone reference standard, BIPM-SRP28, was included in the comparison to verify its agreement with BIPM-SRP27 and thus follow its stability over the period of the on-going key comparison.

5.  Measurement standards

All instruments included in this comparison were Standard Reference Photometers (SRP) built by the NIST. More details on the instrument’s operating principle and capabilities can be found in [1]. The following section describes the SRP measurement principle and uncertainty budget.

5.1.  Measurement equation of a NIST SRP

The measurement of ozone mole fraction by an SRP is based on the absorption of radiation at 253.7nm by ozonized air in the gas cells of the instrument. One unique characteristic of the instrument design is the use of two gas cells to overcome the instability of the light source. The measurement equation is derived from the Beer–Lambert and ideal gas laws. The number concentration ( C ) of ozone is calculated from:

C = 1 2 σ L opt T T std P std P ln D   (1)

where

σ

is the absorption cross-section of ozone at 253.7nm under standard conditions of temperature and pressure, 1.147 6×1017cm2 /molecule [2];

L opt

is the optical path length of one of the cells;

T

is the measured temperature of the cells;

T std

is the standard temperature (273.15K );

P

is the measured pressure of the cells;

P std

is the standard pressure (101.325kPa );

D

is the product of transmittances of two cells, with the transmittance ( T r ) of one cell defined as

T r = I ozone I air   (2)

where

I ozone

is the UV radiation intensity measured from the cell when containing ozonized air; and

I air

is the UV radiation intensity measured from the cell when containing pure air (also called reference or zero air).

Using the ideal gas law, Equation (1) can be reformulated in order to express the measurement results as a mole fraction ( x ) of ozone in air:

x = 1 2 σ L opt T P R N A ln D   (3)

where

N A

is the Avogadro constant, 6.022 142×1023mol1 , and

R

is the gas constant, 8.314 472Jmol1K1 .

The formulation implemented in the SRP software is:

x = 1 2 α x L opt T T std P std P ln D   (4)

where

α x

is the linear absorption coefficient under standard conditions, expressed in cm1 , linked to the absorption cross–section with the relation:

α x = σ N A R P std T std   (5)

5.2.  Absorption cross–section for ozone

The linear absorption coefficient at standard conditions α x used within the SRP software algorithm is 308.32cm1 . This corresponds to a value for the absorption cross–section σ of 1.147 6×1017cm2 /molecule, rather than the more often quoted 1.147×1017cm2 /molecule. In the comparison of two SRP instruments, the absorption cross section can be considered to have a conventional value and its uncertainty can be set to zero. However, in the comparison of different methods or when considering the complete uncertainty budget of the method the uncertainty of the absorption cross–section should be taken into account. A consensus value of 2.12 % at a 95 % level of confidence for the uncertainty of the absorption cross–section has been proposed by the BIPM and the NIST in a recent publication [3].

5.3.  Condition of the BIPM SRPs

Compared to the original design described in [1], SRP27 and SRP28 have been modified to deal with two biases revealed by a study conducted by the BIPM and the NIST [3]. In 2009, an “SRP upgrade kit” was installed in the instruments, as described in [4].

5.4.  Uncertainty budget of the common reference BIPM-SRP27

The uncertainty budget for the ozone mole fraction in dry air ( x ) measured by the instruments BIPM-SRP27 and BIPM-SRP28 in the nominal range 0nmol/mol to 900nmol/mol is given in Table 1.

Table 1.  Uncertainty budget for the SRPs maintained by the BIPM

Component ( y )Uncertainty u y Sensitivity coefficient c i = x y contribution to u x ciuynmol/mol
SourceDistributionStandard UncertaintyCombined standard uncertainty u y
Optical Path L opt Measurement scaleRectangular0.000 6cm 0.52cm x L opt 2.89 × 10 3 x
RepeatabilityNormal0.01cm
Correction factorRectangular0.52cm
Pressure P Pressure gaugeRectangular0.029kPa 0.034kPa x P 3.37 × 10 4 x
Difference between cellsRectangular0.017kPa
Temperature T Temperature probeRectangular0.03K 0.07K x T 2.29 × 10 4 x
Temperature gradientRectangular0.058K
Ratio of intensities D Scaler resolutionRectangular 8 × 10 6 1.4 × 10 5 x D ln D 0.28
RepeatabilityTriangular 1.1 × 10 5
Absorption Cross section σ Hearn value

1.22×1019cm2 /molecule

1.22×1019cm2 /molecule

x α 1.06 × 10 2 x

As explained in the comparison protocol, following this budget the standard uncertainty associated with the ozone mole fraction measurement with the BIPM SRPs can be expressed as a numerical equation (numerical values expressed as nmol/mol ):

u x = 0.28 2 + 2.92 10 3 x 2   (6)

5.5.  Covariance terms for the common reference BIPM-SRP27

Correlations in-between the results of two measurements, performed at two different ozone amount-of-substance fractions with the BIPM–SRP27, were taken into account in the OzonE software. More details on the covariance expression can be found in the protocol. The following expression was applied:

u x i , x j = x i x j u b 2   (7)

where:

u b 2 = u 2 T T 2 + u 2 P P 2 + u 2 L opt L opt 2   (8)

The value of u b is given by the expression of the measurement uncertainty: u b = 2.92 × 10 3 .

5.6.  Condition of the CHMI SRP17

Compared to the original design, the CHMI SRP17 has been modified to deal with the two biases revealed in [3]. In August 2007, an “SRP upgrade kit” was installed by NIST, as already described in the previous comparison report [5].

5.7.  Uncertainty budget of the CHMI SRP17

The uncertainty budget for the ozone mole fraction in dry air x measured by the CHMI standard SRP17 in the nominal range 0nmol/mol to 900nmol/mol is given in Table 2.

Following this budget, as explained in the protocol of the comparison, the standard uncertainty associated with the ozone mole fraction measurement with the CHMI SRP17 can be expressed as a numerical equation (numerical values expressed as nmol/mol ):

u x = 0.28 2 + 2.92 10 3 x 2   (9)

No covariance term for the CHMI standard SRP17 was included in the calculations.

Table 2.  SRP17 uncertainty budget

Component ( y )Uncertainty u y Sensitivity coefficient c i = x y contribution to u x ciuynmol/mol
SourceDistributionStandard UncertaintyCombined standard uncertainty u y
Optical Path L opt Measurement scaleRectangular0.005cm 0.52cm x L opt 2.89 × 10 3 x
VariabilityRectangular0.004cm
DivergenceRectangular0.52cm
Pressure P Pressure gaugeRectangular0.029kPa 0.034kPa x P 3.37 × 10 4 x
Difference between cellsRectangular0.017kPa
Temperature T Temperature probeRectangular0.03K 0.07K x T 2.29 × 10 4 x
Temperature gradientRectangular0.058K
Ratio of intensities D Scaler resolutionRectangular 8 × 10 6 1.4 × 10 5 x D ln D 0.28
RepeatabilityTriangular 1.1 × 10 5
Absorption Cross section α Hearn value

1.22×1019cm2 /molecule

1.22×1019cm2 /molecule

x α 1.06 × 10 2 x

6.  Analysis of the measurement results by generalized least-squares regression

The relationship between the national and reference standards was first evaluated with a generalized least-squares regression fit, using the OzonE software. This software, which is documented in a publication [6], is an extension of the previously used software B_Least recommended by ISO standard 6143:2001 [7]. It includes the possibility to take into account correlations between measurements performed with the same instrument at different ozone mole fractions. It also facilitates the use of a transfer standard, by the handling of unavoidable correlations, which arise since this instrument needs to be calibrated by the reference standard.

In a direct comparison, a linear relationship between the ozone amount-of-substance fractions measured by SRPn and SRP27 is obtained:

x SRP n = a 0 + a 1 x SRP27   (10)

The associated uncertainties on the slope u a 1 and the intercept u a 0 are given by OzonE, as well as the covariance between them and the usual statistical parameters to validate the fitting function.

7.  Results of the comparison

SRP17, maintained by the CHMI, was compared with the SRPs maintained by the BIPM following the general procedure outlined above. A cycle of twelve comparison runs between SRP27, SRP28 and SRP17 were performed. Ozone was generated using the Environics 6100 generator with airflow of 10L/min . The nominal ozone mole fraction range covered during this exercise was 0nmolmol1 to 500nmolmol1 for the eleven first runs, then 0nmolmol1 to 840nmolmol1 for the last run. The result of the last comparison run is presented in detail in the Section 7.1. The repeatability of the results over the twelve runs is shown in Section 7.2, and the stability of the BIPM standards in Section 7.3.

7.1.  Measurement results

The measurement results of the last of the twelve recorded runs are shown in Table 3. For each nominal ozone mole fraction, the standard deviation is reported on the ten successive measurements that were recorded. The values reported here show that both instruments were in a stable regime.

An interesting way to look at these results is to display the difference between the ozone amount-of-substance fractions measured by SRP17 and SRP27 ( x 17 x 27 ) versus the ozone mole fraction measured by SRP27. Figure 1, shows that measurement results from both instruments differ little from each other over the entire concentration range of the comparison. Considering a coverage factor of k = 2 (95 % confidence interval), the difference between both standards is much lower than their combined measurement uncertainties.

Table 3.  measurement results of the comparison between the CHMI standard SRP17 and the BIPM reference standard SRP27

Reference standard BIPM-SRP27 (RS)National standard CHMI SRP17 (NS)
x RS nmol/mol s RS nmol/mol u x RS nmol/mol x NS nmol/mol s NS nmol/mol u x NS nmol/mol
0.180.350.280.040.320.28
201.040.250.65201.710.140.65
72.000.230.3572.260.220.35
386.310.241.16387.340.211.17
108.730.170.42109.170.150.42
292.780.200.90293.600.270.90
25.940.260.2926.040.220.29
339.560.471.03340.710.491.03
154.230.250.53154.800.270.53
465.610.221.39466.930.381.39
247.910.320.78248.570.170.78
0.120.170.28-0.050.270.28
576.50.51.71578.120.381.71
842.10.42.47844.580.332.48
697.90.22.06700.080.222.06
−0.100.100.28−0.040.280.28

Figure 1 — difference between the ozone amount-of-substance fractions measured by SRP17 and SRP27 versus the ozone mole fraction measured by SRP27.

The relationship between SRP17 and SRP27 is given by the result of the generalized least-square regression performed following the method described in Chapter 6:

x SRP17 = 1.003 4 x SRP27 0.06   (11)

The standard uncertainties on the parameters of the regression are u a 1 = 0.003 2 for the slope and ua0=0.21nmol/mol for the intercept. The covariance between the two parameters is cov a 0 , a 1 = 1.66 × 10 4 .

The least-squares regression results confirm that a linear fit is appropriate, with a sum of the the squared deviations (SSD) of 0.30 and a goodness of fit (GoF) equal to 0.17.

To assess the agreement of the standards using Equation (10), the difference between the calculated slope value and unity, and the intercept value and zero, together with their measurement uncertainties, need to be considered. In this comparison, the value of the intercept is consistent with an intercept of zero, considering the uncertainty in the value of this parameter; i.e a 0 < 2 u a 0 , and the value of the slope is consistent with a slope of 1; i.e. 1 a 1 < 2 u a 1 .

7.2.  Repeatability

The results of the twelve comparison runs performed successively between SRP17 and SRP27 are reported in Table 4. With a standard deviation of 0.04nmol/mol on the intercept and a relative standard deviation of 0.03 % on the slope, these data show a good repeatability.

Table 4.  results of the ten comparison runs repeated successively

RunSlope a 1 u a 1 Intercept a 0 / (nmolmol1 ) u a 0 / (nmolmol1 ) cov a 0 , a 1 GoF
11.00270.00330.090.22−2.05 10−40.395
21.00300.0033−0.040.22−2.05 10−40.273
31.00300.00330.040.22−2.04 10−40.391
41.00250.00330.050.22−2.04 10−40.363
51.00250.00330.050.22−2.04 10−40.400
61.00220.00330.090.22−2.04 10−40.756
71.00320.0033−0.010.22−2.04 10−40.277
81.00300.00330.090.22−2.04 10−40.333
91.00220.00330.080.22−2.04 10−40.381
101.00320.00330.000.22−2.05 10−40.473
111.00260.00330.070.22−2.04 10−40.629
121.00320.0032−0.020.18−1.19 10−40.260

7.3.  History of comparisons between BIPM SRP27, SRP28 and CHMI SRP17

Results of the previous comparison performed with CHMI, on the same measurement range and since the first one in December 2002, are shown in Figure 2 together with the results of this comparison. The slopes a 1 of the linear relation x SRP n = a 0 + a 1 x SRP27 are represented together with their associated uncertainties calculated at the time of each comparison. Results of previous comparisons have been corrected to take into account the changes in the reference BIPM-SRP27 described in [4], which explains the larger uncertainties associated with the corresponding slopes.

Figure 2 shows that all standards included in these comparisons remained in close agreement.

Figure 2 — Results of previous comparisons between SRP27, SRP28 and CHMI-SRP17 realised at the BIPM. Uncertainties are calculated at k = 2 , with the uncertainty budget in use at the time of each comparison.

8.  Conclusion

For the sixth time in fifteen years, a direct comparison between the CHMI national standard SRP17 and the BIPM reference standard SRP27 has been conducted at the BIPM. The measurement range of the SRP17 at the CHMI is 10nmol/mol to 870nmol/mol . This entire range was examined during the comparison, and the relative difference between the CHMI SRP17 and BIPM SRP27 measurement results was found to be much smaller than their combined measurement uncertainties.


References

[1]  Paur R.J., Bass A.M., Norris J.E. and Buckley T.J., 2003, Standard reference photometer for the assay of ozone in calibration atmospheres Env. Sci. Technol. NISTIR 6369, 25 pp

[2]  ISO 13964:1996 ISO. ISO 13964 : 1996 Ambient air — Determination of ozone — Ultraviolet photometric method (International Organization for Standardization). 1996.

[3]  Viallon J., Moussay P., Norris J.E., Guenther F.R. and Wielgosz R.I., 2006, A study of systematic biases and measurement uncertainties in ozone mole fraction measurements with the NIST Standard Reference Photometer Metrologia 43 441-450

[4]  Viallon J., Moussay P., Idrees F. and Wielgosz R.I., 2010, Upgrade of the BIPM Standard Reference Photometers for Ozone and the effect on the on-going key comparison BIPM.QM-K1, Rapport BIPM-2010/07, 16pp

[5]  Viallon J., Moussay P., Wielgosz R.I., Novák J. and Vokoun M., 2008, Final report of the on-going key comparison BIPM.QM-K1: Ozone at ambient level, comparison with CHMI, 2007 Metrologia 45, Tech. Supl. 08005

[6]  Bremser W., Viallon J. and Wielgosz R.I., 2007, Influence of correlation on the assessment of measurement result compatibility over a dynamic range Metrologia 44 495-504

[7]  ISO 6143:2001 International Organization for Standardization. Gas analysis — Comparison methods for determining and checking the composition of calibration gas mixtures. 2nd edition. June 2001. Geneva. https://www.iso.org/standard/24665.html.

Document Control

Authors:Philippe Moussay (BIPM), Faraz Idrees (BIPM), Robert Wielgosz (BIPM), Jan Šilhavý (CHMI), and Miroslav Vokoun (CHMI)