^{1}

^{2}

^{2}

^{1}

The increase in the concentration of carbon dioxide (CO2) in the atmosphere has led the scientific community to investigate the adverse effects on humanity and nature, including the greenhouse effect, which contributes to global warming and can lead to climate change, besides the risks associated with human health. Due to the importance of metrological issues in the current scenario, the Laboratory of Gas Analysis (Lanag) of the National Institute of Metrology, Quality and Technology (Inmetro) is developing the methodology of preparation of such primary standard gas mixtures through gravimetry. For the preparation of these standards, an analysis verification step is necessary, carried out by means of an analytical technique of comparison, at low levels of concentration, by cavity ringdown (CRDS) and by gas chromatography using a flame ionization detector coupled to a methaniser catalyst (GC-FIDmeth). This work presents the method validation of both methodologies developed to analyze a range of concentration of atmospheric standards of CO2 at a matrix of synthetic clean dry air (SCDA) by CRDS and GC-FIDmeth, of which the analysis results can compromise the measurements on atmospheric air quality. The objective is to compare results of method validation of both analytical methods for low CO2 concentration through the use of primary reference mixtures developed and also by certified reference material. The validation measurement results were analyzed according to the requirements of ISO 5725: 1994 parts 1 and 2, and the calculation of measurement uncertainty followed the methodology described in ISO 6143:2001, with results showing satisfactory consistent between both selected techniques. According to the work presented here, the obtained validation results for CRDS are better than the GC results, such as the relative uncertainty of samples evaluated for CRDS was 0.4% and 2.7% for GC_FID meth.

The growing requirements concerning the accuracy of measurements and their traceability and worldwide comparability constitute a considerable challenge for the Brazilian National Metrology Institute (NMI), Inmetro, with its responsibility for ensuring the scientific background for the consistency and accuracy of all measurements in Brazil and South America society.

Withal, implementing efficient provisions in order to protect biodiversity and to act against climate change requires an effective and efficient, demand- oriented quality infrastructure (QI). This way, QI with its infrastructure can fulfil obligations laid down in international conventions such as the United Nations Framework Convention on Climate Change (UNFCCC) and use the instruments provided within the scope of the World Trade Organization (WTO). Regional research laboratories over wide will benefit from reliable, accredited measurement and testing capacities which have to be provided to ensure sustainable management.

Carbon dioxide (CO_{2}) the largest contributor to the greenhouse effect, and its concentration in the atmosphere has been monitored using high accuracy techniques, such as gas chromatography, FTIR and cavity ring-down spectroscopy, instead of the established technique used at this type of monitoring―nondispersive infrared gas analyzers (ND-IR). This is need as measurements of CO_{2} in the atmosphere require that uncertainty in the concentration of reference gas mixtures be lower than 0.1 μmol∙mol^{−1} [

Aiming to integrate Inmetro’s mission to promote citizens’ quality of life and the challenges of global metrology with regard to air quality analysis, the Gas Analysis Laboratory (Lanag) has been working for the development and certification of mixtures of standard gases destined to atmospheric and air quality monitoring [_{2} at atmospheric air standard analysis by this specific spectroscopy technique. In this way, this study is proposed for the quantification of the concentration of carbon dioxide in synthetic clean dry air at a range of 370 to 835 μmol∙mol^{−1}, analyzed by both techniques: GC-FIDmeth and CRDS, of which the methods were validated and compared, in order to provide reliability of this type of measurements regarding analysis of greenhouse gases.

Low level concentration of atmospheric carbon dioxide must be accurately quantified using methods with the highest sensitivity and the lowest possible detection limit. The flame ionization detector (FID) detects CH_{4} but CO and CO_{2} have to be reduced to CH_{4} to be detected. The reduction is achieved by using a methanizer composed of a catalyst, usually metallic nickel, at high temperature in presence of hydrogen, placed between the analytical column outlet and the detector [

Online catalytic reduction of carbon monoxide to methane for detection by FID was previously described [

2 CO + 4 H 2 → CH 4 + O 2 (1)

CO 2 + 2 H 2 ↔ 2 CH 4 + O 2 (2)

The catalyst consists of a 2% coating of Ni in the form of nickel nitrate on Chromsorb G. A 1/2” long bed is packed around the bend of an 8” × 1/8” SS U-tube. The tube is clamped in a block so that the ends protrude down into the column oven for easy connection between column or TCD outlet and FID base. Heat is provided by a pair of cartridge heaters and controlled by a temperature controller. Hydrogen for the reduction is provided by adding it via a tee at the inlet to the catalyst.

Methanizer efficiency is important to monitor when analyzing trace gases with GC-FID. It can be calculated by comparing the CO_{2} and CH_{4} peak ratio against concentration ratio and can be compromised by high concentrations of hydrocarbons blocking the active sites on the catalyst. The degradation of efficiency is said to be reversible using air as a carrier gas (as shown in Equation (3)) to burn off the carbon film from the catalyst by bypassing the column first [

3 O 2 + 2 C + 2 H 2 → 2 CO 2 + 2 H 2 O (3)

The GC system used was a Varian GC3800 equipment equipped with a column Carbobond (Varian) and a flame ionization detector (FID) coupled to a methanizer catalyst. The method parameters for analysis of carbon dioxide scheme of the GC at atmospheric level range are presented on

Nearly every small gas-phase molecule (e.g., CO_{2}) has a unique near-infrared absorption spectrum. At sub-atmospheric pressure, this consists of a series of narrow, well-resolved, sharp lines, each at a characteristic wavelength. Because

Injectors | |
---|---|

Split ratio | 1 |

Heater Temperature | 150˚C |

Column specifications | |

Column Flow | 20 mL/min |

Column Pressure | 5 psi |

Column Temperature | 50˚C |

Column type | CP-Carbobond Fused silica |

Lengh | 50 m |

Inside diameter | 530 µm |

DF | 10 µm |

Carrier gas | Helium |

FID | |

Temperature | 250˚C |

N2 makeup flow | 25 mL/min |

H2 flow | 30 mL/min |

Air flow | 300 mL/min |

Methanizer temperature | 380˚C |

Running time | 10 min |

CO_{2} Retention time | 6.7 min |

these lines are well-spaced and their wavelength is well-known, the concentration of any species can be determined by measuring the strength of this absorption, i.e. the height of a specific absorption peak. But, in conventional infrared spectrometers, trace gases provide far too little absorption to measure, typically limiting sensitivity to the parts per million at best. Cavity Ring-Down Spectroscopy (CRDS) avoids this sensitivity limitation by using an effective pathlength of many kilometers. It enables gases to be monitored in seconds or less at the parts per billion level [

Thus, CRDS works by attuning light rays to the unique molecular fingerprint of the sample species. By measuring the time it takes the light to fade or “ring- down”, you receive an accurate molecular count in milliseconds. The time of light decay, in essence, provides an exact, non-invasive, and rapid means to detect contaminants in the air.

Regarding the principle of the CRDS to generate the result expected, a computer-controlled system tunes the laser off the absorption peak for the sample species to determine the τ empty value, equivalent to a zero baseline correction. It tunes back to the absorption peak to determine the τ value, dependent on the sample species concentration. Based on Beer’s Law, this value constitutes an absolute measurement and is unaffected by losses outside the ring-down cavity. The graph presented at

Science should be used to improve quality of life. In this regard, scientists must guarantee the quality of their results. These results must provide relevant information, not just for the scientific community, but also for citizens. Demonstration of the ability of an analytical method to provide reliable results is of great importance to ensure quality, safety and efficacy in the analysis. Consequently, before an analytical method is implemented for routine use, it must be previously validated to demonstrate that it is suitable for its intended purpose. The

final goal of the validation of an analytical method is to ensure that every future measurement in routine analysis will give a true assessment of the actual value of the content of an analyte in the sample [

All accredited Brazilian test laboratories shall comply with the requirements of ABNT NBR ISO/IEC 17025 [

With the methods developed for the two techniques evaluated, calibration curves were drawn for each selected range in a given technique, using the primary standards produced in-house to compose the curve and as samples, as well as, certified reference materials to verify the adjust of the proposed calibration curve [_{i}, y_{i}), which the known concentrations used are x_{i}, are plotted against the responses of the instrument, y_{i}, obtained at independent and repetitive conditions. A linear relationship between the concentrations and the measurement results is adjusted, obtained by the mathematical model of correlation: least squares numerical method [

The validation of the CRDS analysis was performed using independent results under conditions of repeatability (5 replicates) and reproducibility (3 different days). It was selected 03 (three) primary standard gaseous mixtures produced by Lanag, which range of concentration adopted was from 370 to 420 μmol∙mol^{−1}, due to the specifications from the equipment used. Those standards were analyzed to be fitted in a linear model of calibration curve, and two standards were selected as samples: one developed by Lanag at a mole fraction of 380 μmol∙mol^{−1}, and the other, a CRM from the renowned Dutch Metrology Institute, VSL, at a mole fraction of 300 μmol∙mol^{−1}, both at the lower limit of the range, in order to determine the limits of quantification of the technique evaluated.

Regarding the validation analysis of the GC using a methanizer catalyst before the FID detector was performed using independent results under conditions of repeatability (4 replicates) and reproducibility (2 different days). As it’s known that chromatographic areas results are non linear, specially if a wide range is fitted, it was selected 05 (five) primary standard gaseous mixtures from Lanag, which range of concentration varied from 450 to 835 μmol∙mol^{−1}, as the specifications from the equipment allows to go further at the higher limit of the range. Those standards were analyzed to be fitted in a quadratic model of analytical calibration curve, and two standards were selected as samples: one produced by Lanag at a mole fraction of 555 μmol∙mol^{−1}, and the other, a CRM, also from VSL, at a mole fraction of 600 μmol∙mol^{−1}, both at the middle of the range in order to prove the good adjustment of the function model specified.

The linearity was assessed by repeated injections of five primary standards produced at different concentrations englobing the selected range. The linearity is evaluated by the following statistical approach: the coefficient of correlation and the goodness of fit, derived by the part of validation of the response model at ISO 6143. To effectively test the compatibility of a prospective analysis function, calculate the measure of goodness-of-fit (GOF), defined as the maximum value of the weighted differences, | x ^ i − x i | / u ( x i ) and | y ^ i − y i | / u ( y i ) , between the coordinates of measured and adjusted calibration points ( i = 1 , 2 , 3 , ⋯ , n ). A function is admissible if GOF < 2, as well, as r2 > 0.99.

CRDS | ||||
---|---|---|---|---|

Primary standard mole fraction (μmol∙mol^{−1}) | Gravimetric uncertainty (μmol∙mol^{−1}) | Day 1―CV (%) | Day 2―CV (%) | Day 3―CV (%) |

370.55 | 0.08 | 0.21 | 0.06 | 0.05 |

390.61 | 0.08 | 0.23 | 0.60 | 0.02 |

421.32 | 0.08 | 0.03 | 0.16 | 0.08 |

GC-FIDmeth | ||||

Primary standard mole fraction (μmol∙mol^{−1}) | Gravimetric uncertainty (μmol∙mol^{−1}) | Day 1―CV (%) | Day 2―CV (%) | Day 3―CV (%) |

453.33 | 0.09 | 1.07 | 1.45 | - |

552.33 | 0.11 | 0.86 | 0.77 | - |

789.11 | 0.25 | 0.69 | 0.99 | - |

806.33 | 0.26 | 0.79 | 0.43 | - |

835.22 | 0.24 | 0.32 | 0.62 | - |

Results from the primary standards used as samples at the curve of calibration fitted with the responses from the primary standards obtained from each technique can be seen at _{c}) obtained by the regression from the gravimetric mole fraction (x_{g}); the relative expanded uncertainty (U) calculated by the calibration curve fitted; and the number of response repetitions (N) of each day analyzed at each technique evaluated.

The regression analysis of all data obtained from both techniques was made using the trend line regression tool from Microsoft Excel. And the GOF was obtained by applying the regression based on ISO6143 [

In order to verify if the regression from the calibration curve adopted is significant, it was performed an analysis of variance (ANOVA) that derived the results for the Ftest that is determined by the relation of the quadratic average of the regression and the quadratic average of the residuals. This F value is compared with the value of Fcritic tabulated at the selected confidence level (95%). If F > Fcritic, it is accepted, at the selected confidence level, that a ¹ 0, which means that the slope of the regression line is not zero, that is, the regression is significant. If F £ Fcritic, there is no indication of a linear relationship between the variables x (concentration values) and y (measurement responses) [

Another evaluation made from the linear regression of data analysis, presented at

Precision evaluates the dispersion of results between independent and repeated assays. Repeatability is the degree of agreement between the results of successive measurements under the same measurement conditions. Reproducibility is the degree of agreement between the results of successive measurements

CRDS | ||||||
---|---|---|---|---|---|---|

Day | x_{g} (μmol∙mol^{−1}) | x_{c} (μmol∙mol^{−1}) | Δ (%) | U (%) | CV (%) | N |

1 | 380.09 | 378.82 | 0.32 | 0.42 | 0.14 | 05 |

2 | 380.09 | 380.12 | −0.02 | 0.48 | 0.25 | 05 |

3 | 380.09 | 380.12 | −0.02 | 0.22 | 0.21 | 05 |

GC-FIDmeth | ||||||

Day | x_{g} (μmol∙mol^{−1}) | x_{c} (μmol∙mol^{−1}) | Δ (%) | U (%) | CV (%) | N |

1 | 555.20 | 554.40 | 0.14 | 2.31 | 1.06 | 04 |

2 | 555.20 | 552.53 | 0.48 | 3.17 | 1.08 | 04 |

Parameters of Methods Validation results | CRDS | GC-FIDmeth | |||
---|---|---|---|---|---|

Day 1 | Day 2 | Day 3 | Day 1 | Day 2 | |

GOF | 1.58 | 0.23 | 2.16 | 1.33 | 1.91 |

r2 | 0.9991 | 0.9997 | 0.9991 | 0.9972 | 0.9960 |

Standard error | 5.92 | 3.41 | 7.67 | 201.97 | 220.02 |

Observations | 43 | 43 | 39 | 20 | 20 |

F | 4783.9 | 3230.6 | 32115.0 | 3913.9 | 2968.5 |

Fcritic | 4.5 × 10^{−18} | 5.7 × 10^{−17} | 1.9 × 10^{−23} | 1.6 × 10^{−22} | 1.9 × 10^{−21} |

RSD (%) | 0.14 | 0.25 | 0.21 | 1.06 | 1.08 |

s_{r} (%) | 0.20 | 1.07 | |||

s_{R} (%) | 0.53 | 2.07 | |||

ER (%) | 0.09 | 0.3 | |||

EN (%) | 0.51 | 0.10 | |||

Z-score | 0.24 | 0.11 | |||

LOD (μmol∙mol^{−1}) | 16 | 40 | |||

LOQ (μmol∙mol^{−1}) | 50 | 121 |

under varying measurement conditions. Reproducibility was determined internally, only by Inmetro’s laboratory, by injecting patterns and samples by different analysts on different days (intermediate precision). The repeatability and the intermediate precision were evaluated by assessing the following statistical approaches: relative standard deviation (RSD %) or coefficient of variation (CV), relative deviation standard for repeatability (s_{r}), relative deviation standard for reproducibility (s_{R}). The acceptance criteria for the parameter quoted before are RSD < 1%, s_{r} < 2%, and s_{R} < 5%.

Another parameter evaluated was the repeatability limit (r), which is given by 2.8 times the relative standard deviation of the repeatability (s_{r}), for a significant level of 95%. The reproducibility limit (R) is given by 2.8 times the square root of the relative standard deviation of the reproducibility (s_{R}).

Accuracy is the agreement between the result of a test and the reference value accepted as true, and for its evaluation the value of the relative difference between the gravimetric and the analytical concentration derived from the calibration curve (Δ) from the results through the Xgenline software, as it can be seen on

・ 90% < Recovery < 110%

・ E R ( % ) = x c − x g x g ⋅ 100 (4)

・ E N ( % ) = ( x c − x g ) U c 2 + U g 2 < 1 (5)

・ Z -score = ( x c − x g ) s < 2 is satisfactory, (6)

where, x_{c} is the analytical and x_{g} the gravimetric concentration or mole fraction, U_{c} and U_{g} the expanded uncertatinty of analysis and gravimetry respectively, and s can be the standard deviation from the analysis responses or even the combined uncertainty from the samples evaluated.

Limit of detection (LOD) is the lowest concentration of analyte that can only be detected. Limit of quantification (LOQ) is the lowest concentration of analyte that can be quantified as an exact value with precision and accuracy. There are two main methods used to determine the LOD and LOQ: the signal-noise relation method and the method based on the standard deviation of response and slope of the calibration curve [

・ LOD = 3.3 × Stddeviation slope ; and (7)

・ LOQ = 10 × Stddeviation slope (8)

As observed from both techniques most of the days analyzed, r2 > 0.99 and the GOF < 2, which means all calibration curves were adjustable and satisfactory. The third day of CRDS analysis had a bad adjustment of the curve fitted. The best function adjustment for GC was a quadratic model, as presented better GOF results when it was adjusted by ISO 6143 regression than the r2 results obtained by the linear regression fit.

According to

Regarding the repeatability results from the GC analysis days are higher than the established criteria, i.e., relative standard deviation higher that 1%, which means that both process and/or method should be optimized and more replicates should be taken. Nevertheless, results from repeatability and reproducibility among days are lower than the established criteria, although results of GC are much higher than CRDS data.

According to all days of data for both methods analyzed recovery results were between 99.5% and 100.5%, which means the expected calibrated results and the real ones from gravimetric concentration agree on accuracy, as well as, it was satisfactory all other parameters tested.

The mole fraction obtained for LOD and LOQ evaluation results are considered acceptable, as are much lower than the range of concentration selected for each technique.

In summary, most of parameters evaluated presented an accepted acceptance criteria and the method validation for both techniques were considered satisfactory.

Finally, the results obtained on different days were also compared through ANOVA single factor test, which performs the analysis of simple data variance of two or more samples. The acceptance criteria used was the following: if F < Fcritic and P-value > 0.05, the average through different days of the sample are equivalent. This analysis tests the hypothesis that each sample comes from the same baseline probability distribution. Anova single factor results are presented on

According to ANOVA single factor criteria, results of average for GC-FID- meth are slightly not equivalent, as F < Fcritic and P-value a little bit higher than 0.05. On the other hand, CRDS Anova results are considered satisfactory for both criteria evaluated.

All data obtained from GC results can also be compared and validated to Lanag’s previous results on the key-comparisons participated among other national metrology institutes for the parameters evaluated [

Method validation using primary standard gas mixture was applied to carbon dioxide at atmospheric air levels. Preliminary results of key-comparisons of this gas component at such nominal concentration presented different estimations for uncertainty calculation, but usually GC-FIDmeth was used. There are no previous studies published so far for the analysis of this standard by CRDS, which is being used by known national metrology institutes nowadays. So, this paper presented the comparison of the results from a validation of method for GC-FIDmeth and CRDS applied to a range of 370 to 835 μmol∙mol^{−1} of carbon dioxide in synthetic clean dry air primary standard mixture. When compared linearity and calibration curve adjustment, it can be seen that CRDS and GC are equivalent, with GOF for both lower than the criteria value of 2. The relative

CRDS | |||||||
---|---|---|---|---|---|---|---|

Group | Counting | Sum | Average | Variance | F | Fcritic | P-value |

Day 1 | 5 | 1977630.6 | 395526.1 | 297412.1 | 54.2 | 3.9 | 9.8 × 10^{−7} |

Day 2 | 5 | 1995414.1 | 399082.8 | 997612.9 | |||

Day 3 | 5 | 2003868.3 | 400773.7 | 689556.7 | |||

GC-FIDmeth | |||||||

Group | Counting | Sum | Average | Variance | F | Fcritic | P-value |

Day 1 | 4 | 41401.6 | 10350.4 | 11975.7 | 5.4 | 6.0 | 0.06 |

Day 2 | 4 | 40681.2 | 10170.3 | 12139.4 |

standard deviation results for CRDS is four times lower than GC. The repeatability and reproducibility results for CRDS are better than the GC as well. Another parameter evaluated was the limit of quantification (LOQ), which for CRDS was obtained 50 μmol∙mol^{−1} and for GC 120 μmol∙mol^{−1}. According to ANOVA single factor criteria, results of average for GC-FIDmeth are slightly not equivalent, as F < Fcritic and P-value a little bit higher than 0.05. On the other hand, CRDS Anova results are considered satisfactory for both criteria evaluated. Finally, when comparing relative expanded uncertainty from a sample applied at the calibration curve developed for each method, GC-FIDmeth presented an estimation of 2.7%, while CRDS presented around 0.4%. It can be seen that the GC-FIDmeth method developed shall be optimized in order to get better uncertainty estimations, especially if compared this GC method evaluation to preliminary GC results from key-comparisons, which average relative uncertainty obtained was 10 times lower than the one presented at this study of 2.7%. In this way, the results obtained by comparing both validations of the methods developed for each technique conclude that CRDS, based on the robustness of the spectroscopy technique, has better results than GC-FIDmeth, especially when uncertainty is the criteria of decision. This conclusion of this study of comparison is confirmed considering the high application of CRDS when atmospheric monitoring is applied, allowing the achievement of low detection limits allied to high level of uncertainty at a good speed of response.

Fioravante, A.L., Augusto, C.R., Cunha, V.S. and d’Avila, L.A. (2018) Validation and Comparison of Calibration Techniques for Measurements of Carbon Dioxide in Atmospheric Air Standards. Journal of Environmental Protection, 9, 34-49. https://doi.org/10.4236/jep.2018.91004