The Importance of Water Vapor Measurements and Corrections

Printable PDF: The Importance of Water Vapor Measurements and Corrections

Application note that describes how water vapor influences the measurement of carbon dioxide by infrared gas analysis in several ways, which can lead to significant measurement errors.

Water vapor is known to influence the measurement of carbon dioxide by infrared gas analysis in several ways, which can lead to significant measurement errors. Spectral cross-sensitivity due to absorption band broadening, and inherent instrument cross sensitivity can both cause over estimations of CO2 mole fraction in samples containing water vapor when their effects are not accounted for. Dilution of samples by the addition of water vapor may not be important when measuring actual CO2 mole fractions, but can lead to significant errors in flux measurements. In this note we describe the basis of each of these three processes and discuss how each can affect the measurement of CO2.

Infrared Absorption by Gases and Absorption Band Broadening

At the sub-molecular scale the positions of atoms within molecules are not entirely fixed. As they move they stretch and bend their bonds, creating vibrations within the molecule. Energy differences between the possible vibrational states that result from this movement make it possible for the molecule to absorb infrared radiation. Changes in angular momentum as the molecule rotates about its axis can also cause infrared radiation to be absorbed. The energy differences resulting from changes in the vibrational and rotational states of the molecule cause fluctuations in its dipole moment. These oscillations interact with the alternating electrical field of electromagnetic radiation and if the frequencies of oscillation match, the radiation will be absorbed by the molecule.

Since the frequency of oscillation for electromagnetic radiation (ν) is the inverse of its wavelength (λ=1/ν) and the various energy states that lead to its absorption by a molecule are the result of the molecule's structure, the absorption of infrared radiation is both wavelength and absorber species dependent. For a given molecular species, absorption of infrared radiation will occur in bands at various wavelengths across the infrared region of the spectrum (Figure 1‑1, A). Each of these absorption bands is comprised of individual absorption lines (Figure 1‑1, B) that result from rotational transitions, and have a generally Lorentzian line shape described by

1‑1

where k(ν) is the absorption coefficient k of radiation at frequency ν, S is the line strength, (ν-ν0) is the change in frequency across the absorption line and α is the absorption line half width.

Figure 1‑1. Absorption spectrum of CO2. Carbon dioxide absorption data from HITRAN96 (Rothman et al. 1998). A: The absorption spectrum of CO2 in the infrared region of the electromagnetic spectrum. B: The absorption of CO2 in the 4.2 μm region.

For infrared gas analysis with LI-COR gas analyzers, we are interested in measuring the total absorption of infrared light across a given waveband to determine absorber species concentration. The total absorption measured by the analyzer, A is in principle the result of integrating the absorption across the entire band as shown by

1‑2

where A'(ν) is the absorption of infrared radiation at a given radiation frequency, and is given by Beer-Lambert's Law

1‑3

where w is the absorber concentration and l is the optical path length (Burch and Williams 1964). These relationships form the basic principles behind the calibration functions employed in LI-COR gas analyzers, because they allow us to relate the measured absorption to the absorber concentration. Because the relationship between w and A'(ν) is influenced by k(ν), anything that influences k(ν) will affect A and has to be accounted for in the calibration function.

Recalling from equation 1‑1, the absorption coefficient k is a function of radiation frequency, and is a product of the line strength, S, and absorption line half width, α. For the range of conditions typically encountered in IRGA measurements, we can treat S as a constant for a given frequency (Burch et al. 1962). However, since individual absorption features are the product of changes in energy states anything that induces a change in energy state will influence α. From Kinetic Theory it can be shown that α is a function of intermolecular collision frequency:

1‑4

where ρi is the number of molecules of type i in a given volume, Da,i is the sum of the optical collision diameters of the absorber molecule and a molecule of type i, and ma and mi are the masses of the absorber molecule and a molecule of type i, respectively. Equation 1‑4 indicates that there are a handful of physical factors that influence the width of spectral lines, and shows the theoretical basis for band broadening (also referred to as collisional or pressure broadening) of spectral features.

The density of component i in a gas mixture, ρi is a product of both the pressure (P) and temperature (T) in the gas mixture as given by the Ideal Gas Law:

1‑5

where R is the ideal gas constant. As a result the width of individual spectral lines changes with both T and P. However, the total absorption of an individual line does not change. As the lines become broader, tightly spaced lines comprising absorption bands overlap. This causes the total absorption of the band to increase, despite no change in the absorption of each of the individual absorption lines comprising the band.

Effective pressure: The terms Da,i, ma and mi in equation 1‑4 vary depending on the species of absorber molecule and the interacting molecules in the gas mixture. Thus, not all molecular species are equally effective at broadening the absorption lines of another species and their total broadening effect is influenced by their abundance in the mixture. This leads to the concept of effective pressure, where the total broadening potential of a component in a gas mixture is a function of both its partial pressure (pi) and a weighting factor (ai) as given by

1‑6

where PEi is the effective pressure of component i in the mixture. The weighting factor ai, is generally expressed relative to the broadening potential of nitrogen (aN2=1), and is larger for molecules similar to the absorber species (Table 1).

Table 1‑1. Example broadening coefficients for some common atmospheric gas species. Data are taken from Bruch et. al., 1962 and are referenced to the broadening potential of N2. Band refers to the absorbing gas absorption band used to determine the broadening coefficients.
Absorbing gas Band (cm-1) Broadening gas ai
N2O 2224 O2
CH4
N2O
0.83
1.08
1.12
CO 2143 CH4
CO
1.12
1.02
CO2 2350 CO2 1.3
CH4 3020 CO2
CH4
1.25
1.3

In a gas mixture, both the species of interest and all other species in the mixture can exert a broadening effect on the measurement. For measurements of CO2 in moist air the total effective pressure (PE) of the gas mixture can be written as

1‑7

where pc is the partial pressure of CO2, ac is the self broadening coefficient for CO2, pw is the partial pressure of water vapor, aw is the broadening coefficient for water vapor on CO2, pa is the partial pressure of dry air and aa is the broadening coefficient for dry air on CO2. Over the range of partial pressures typically encountered in atmospheric measurements, the partial pressure of CO2 is so small that its self-broadening potential can be neglected. LI-COR calibrates all of its analyzers using mixtures of CO2 in air. If individual components in dry air are constant then aa can be taken as 1 and equation 1‑7 rewritten as

1‑8

Equation 1‑8 is implemented in the calibration functions used by most LI-COR gas analyzers in a slightly different form. Given that the total pressure (P) of a gas mixture is the sum of the partial pressures of its components, neglecting CO2

1‑9

By dividing by P, equation 1‑9 can be rewritten in units of mole fraction

1‑10

where Xa is the mole fraction of dry air and Xw is the mole fraction of water vapor. Substituting from equation (10), equation (8) can be rewritten in such a way that only the mole fraction of water vapor and the total pressure need to be known to correct the measured absorption for band broadening:

1‑11

where the function . Equation 1‑11 gives the form of the band broadening correction applied in the calibration functions of most LI-COR gas analyzers. The exceptions are those analyzers which do not measure water vapor, such as the LI-800 and LI-820. The LI-6400/LI- 6400XT includes an additional broadening term based on user specified oxygen content of the sample gas. This makes possible accurate measurements of photosynthesis under non-photorespiratory conditions.

Theory suggests the value of aw should be constant and could be applied to all gas analyzers. In practice however, we find that when aw is empirically derived from measurements made by a gas analyzer, its value is unique to that type of gas analyzer. This discrepancy results from differences in the optics employed in the different analyzers, and is a property of their spectral band width and transmittance. A summary of the band broadening coefficients employed in LI-COR gas analyzers is given in the appendix. It should be noted that these values are different from values of aw published in the literature, both because they are analyzer specific and because they use dry air rather than nitrogen as a reference.

The effect of band broadening on the measurement is relatively small, but is still significant. Band broadening causes an increase in the apparent CO2 mole fraction, as shown by the uncorrected data in Figure 2. This occurs because a gas mixture containing water vapor will always have a larger value for PE as compared to P. While the errors in Figure 1‑2 are only on the order of one percent, they can still be reduced significantly by applying a band broadening correction.

Figure 1‑2. The effect of band broadening corrections on CO2 measurements. Dilution corrected measurements of CO2 mole fraction were made with an LI-6262 at 31°C and 48°C, over a range of pressures (70.5 to 97.2 kPa) and vapor pressures (0.89 to 2.46 kPa) (data from McDermitt et al. 1993). Without a band broadening correction (black circles) the measured CO2 mole fraction was consistently higher than the actual mole fraction. When the band broadening correction was applied (blue circles), the differences between measured and actual values were typically one μmol/mol or less.

Instrument Cross-sensitivity

Cross-sensitivity due to direct absorption: In its simplest form the optical bench of an infrared gas analyzer (IRGA) consists of a source of infrared radiation and means of detecting radiant power, with a gas sample positioned in between the two. In practice however, the optical bench of an IRGA (Figure 1‑3) designed to measure CO2 and water vapor is more complicated because both species absorb infrared radiation in the same general region of the spectrum. Instruments employing broad band sources, such as the heated filaments used in most LI-COR gas analyzers, use optical filters to isolate absorption bands specific to the absorber species of interest. If absorption features from other gases fall in the transmission band of these filters, the instrument will show a cross sensitivity to those gases. Additionally, while these optical filters have a defined bandwidth of peak transmittance, they can show some slight out of band transmission. When this out of band transmission corresponds to a region of absorption for another gas species the instrument can show cross-sensitivity to that gas species.

The effect of cross-sensitivity due to direct absorption is an over estimation of the absorber species concentration when the interfering gas is present in the sample. For instruments that show this type of cross-sensitivity it can be corrected for in the instruments' calibration function.

Figure 1‑3. IRGA design schemes. A. Mechanical differential IRGA utilizing a chopper wheel. B. An optical differential IRGA utilizing a filter wheel. C. A simple optical differential IRGA using fixed filters. 1. Infrared source. 2. Focusing lens. 3. Sample cell. 4. Reference cell 5. Mirror. 6. Water vapor inband filter. 7. CO2 in-band filter. 8. Detector. 9. Dichroic mirror. 10. Chopper wheel. 11. Drive motor. 12. Window. 13. Filter wheel. 14. CO2 out-of-band filter. 15. Water vapor out-of-band filter.

Electrical cross-sensitivity: Since absorption is not directly measureable, it is calculated from the transmittance of radiation through the sample as

1‑12

where Φi is the transmitted radiant power in the presence of the absorber species and Φo is the transmitted radiant power in the absence of the absorber species. The terms Φi and Φo can be determined from measurements at two different wavebands, one centered on an absorption band of the target species (in-band) and one centered in a region of no absorption by the target species (out-of-band), or by comparing measurements between two separate gas samples. This leads to two basic approaches to the instrument design, utilizing either a mechanical or optical differential to determine Φi and Φo (Figure 3).

In instruments using an optical differential approach, the in and out-of-band filters maybe mounted on a spinning filter wheel that moves them sequentially in line with the detector (Figure 3B). The timing of the wheel is precisely controlled such that when the appropriate filter is in line the instrument makes a measurement of Φi or Φo for either CO2 or water vapor. In this scheme the time response of the detector can cause an additional form of cross-sensitivity. The detectors used in IRGAs capture radiant energy and convert it to an electrical current. These detectors can show a small lag in the signal rise and decay from the start and stop of the radiation (i.e. the time period between when a filter comes in and out of line with the detector). While the time span of this lag is extremely small, when the detector output is sampled at high frequency the decay lag from one filter can overlap the next filter. In this situation the instrument can show a cross-sensitivity between the CO2 and H2O measurements. However, because the overlapping signal is small the magnitude of this cross-sensitivity is also small.

Some LI-COR gas analyzers are known to show a slight cross-sensitivity in their measurements between CO2 and water vapor. Those that do include a cross-sensitivity term in their calibration functions, which corrects for this effect. Because the processes leading to the cross-sensitivity are due to component-specific properties, each analyzer is individually characterized for cross-sensitivity during the calibration procedure. A summary of which LI-COR analyzers include a cross-sensitivity correction is included in the appendix. For these analyzers the cross-sensitivity correction is significantly smaller than that for band broadening.

Dilution by Foreign Gases

Unlike the broadening of spectral lines by foreign gases or instrument cross-sensitivity, dilution of samples by foreign gases does not cause errors in the actual measurement of the absorber gas. Rather the correction of concentration data for dilution effects yields a theoretical concentration of a component in a gas mixture assuming some concentration of the diluting gas. As such, adjusting for dilution effects is not important in all measurements. For measurements where the goal is simply to know the CO2 content of a given gas sample, a dilution correction may not be needed. However, in situations where comparisons are to be made between mole fractions of one component in varying gas mixtures, correcting for the dilution effect can be very important (Penning et al. 1984). A common example is when calculating carbon fluxes from measurements of CO2 mole fractions in air. If the system used to make the flux measurements operates under constant pressure and the sample is not dried before reaching the IRGA, the influx of water vapor to the system due to transpiration and/or evaporation will cause a proportionate decrease in the sample CO2 concentration (Figure 1‑4, A). Without applying a dilution correction, under these conditions carbon uptake will be overestimated and carbon losses will be underestimated. Figure 1‑4, B shows corresponding flux calculations for both diluted measurements and dilution corrected measurements. It is clear that the uncorrected measurements indicated a sizable carbon uptake, even though the sample material was inert.

Figure 1‑4. An example of dilution in an open flow-through system. A: Sample chamber CO2 mole fractions (blue circles) were measured with an inert transpiring body (wet filter paper) enclosed in an O-ring sealed chamber attached to the LI-6400XT. During measurements the chamber was supplied with a gas of known CO2 mole fraction (400 μmol/ mol) and transpiration was varied by adjusting the initial moisture content of the air stream. The decrease in sample chamber CO2 mole fraction is real and becomes larger as the transpiration rate increases. B: The apparent assimilation rate (gray circles), calculated without a dilution correction, also increases with transpiration rate, while the actual assimilation rate reported by the instrument (open circles) remained near zero for all the measurements.

In practice we are often interested in knowing the mole fraction of CO2 in the absence of water vapor, or standardizing two CO2 measurements to the same water vapor content. In the first case, equation (1‑20) can be written in standard units of measurement for the two gases as

1‑13

where c is the mole fraction of CO2 (μmol/mol) in dry air, c' is the measured CO2 mole fraction (μmol/mol) of wet air, and w is the measured water vapor mole fraction (mmol/ mol). In the latter case, correcting one measurement to the water vapor concentration in another is done by multiplying one of the measured CO2 mole fractions by the ratio of the bracketed term in equation (1‑17) for the two measurements, as shown by

1‑14

where is the mole fraction of CO2 in the sample (μmol/mol) corrected to the water vapor content of the reference measurement wr (mmol/mol), is the mole fraction of CO2 measured in the sample (μmol/mol), and ws is the water vapor content in the sample (mmol/mol).

For flux measurements, equation (1‑13) or (1‑14) may be applied to raw concentration data or a dilution correction may be embedded within the flux equation. Figure 5 provides examples of some common flux measurements and how a dilution correction is applied in each. The appendix includes a summary of which LI-COR gas analyzers apply dilution corrections and how the correction is implemented.

Deriving the Dilution Correction

For any gas mixture at constant temperature and pressure, an increase in the quantity of one of its components will result in an equivalent decrease in the quantity of the remaining components. This occurs because, as shown by the Ideal Gas Law at any given temperature and pressure, only a certain number of molecules (N) can occupy a given volume (V)

1‑15

If N represents the total number of molecules in a multicomponent gas mixture, then N is equal to the sum of the quantities of its components (ni), in moles

1‑16

If some quantity of a diluting gas (nd) is added to this mixture, while pressure and temperature remain constant, then the sum of ni must be reduced by an equivalent amount as shown by

1‑17

where n’i is the diluted molar quantity of component i in the mixture. The term in brackets in equation (1‑17) represents the proportion of the total number of moles in the mixture comprised by all components excluding the diluting gas. Substituting for N from equation (1‑16) yields

1‑18

which is solvable for the number of moles of all components if the diluting gas were not present. Equation (1‑18) can be rewritten in terms of mole fraction by dividing by N to yield

1‑19

where the mole fraction of component i, Xi, is equal to ni/N. This can be written for a single component of the gas mixture as

1‑20

to give the mole fraction of component i in the absence of the diluting gas.

Table 1‑2. How dilution corrections are applied for some common flux measurements. For an open flow-through system such as the LI-6400XT, the correction is used to adjust the sample cell CO2 mole fraction to the water vapor content in the reference cell. For a closed transient system like the LI-8100 and for Eddy Correlation measurements using a closed path analyzer like the LI-6262, the correction can be applied to the measured CO2 mole fraction before the flux calculation is done. For Eddy Correlation measurements made with an open path analyzer such as the LI-7500, the dilution correction occurs as part of the Webb-Pearman-Leuning (WPL) correction which corrects for apparent fluxes due to both heat and water vapor (Webb et al. 1980). For a summary of the derivations of these equations, see Welles and McDermitt 2005.
    Flux Equation Dilution Term
Open System
Eddy Correlation Open Path Analyzer
Closed System
Eddy Correlation Closed Path Analyzer
Fc flux of CO2; f flow; X’cr reference cell CO2 mole fraction; X’cs sample cell CO2 mole fraction; Xwr reference cell water vapor mole fraction; Xws sample cell water vapor mole fraction; s area; Xcswr dilution corrected sample cell CO2 mole fraction; dXc/dt change in dilution corrected CO2 mole fraction with time; Xc dilution corrected CO2 mole fraction; X’c measured CO2 mole fraction; Xw water vapor mole fraction; w' fluctuation in mean vertical wind speed; ρ’c fluctuation in mean CO2 density; μ ratio of molecular weights of dry air and water vapor; ρc mean CO2 density, ρa mean dry air density; ρv mean water vapor density; T' fluctuation in mean air temperature; T mean air temperature

Appendix: Summary of water vapor corrections applied in LI-COR gas analyzers

Instrument Species detected Band broadening coefficient (aw) Cross sensitivity Dilution correction Implementation
LI-800 CO2 None   No  
LI-820 CO2 None   No  
LI-840 CO2 and H2O 1.45 H2O on CO2 No  
LI-840A CO2 and H2O 1.451 H2O on CO2 No  
LI-62002 CO2 None   Yes Carbon dioxide mole fraction is not dilution corrected, but a transpiration correction is applied in the photosynthesis equation.
LI-6251 CO2 1.57   No Band broadening and dilution corrections are not applied by the instrument. They can be applied to the data post collection following the instructions in the user manual.
LI-6252 CO2 1.57   Yes (optional) Band broadening and dilution corrections can be applied if the water vapor concentration is supplied to the instrument via an external sensor. The dilution correction is the same as applied in the LI-6262.
LI-6262 CO2 and H2O 1.57   Yes Using FCT 76 a dilution correction can be applied, which adjusts the sample CO2 mole fraction to the same water vapor content as that in the reference cell.
LI-6400/XT CO2 and H2O 1.5   Yes A dilution correction is applied in the photosynthesis equation. When configured as a closed system, the instrument calculates the CO2 mole fraction relative to dry air.
LI-7000 CO2 and H2O 1.57   No  
LI-7200/RS CO2 and H2O 1.15 CO2 on H2O
H2O on CO2
Yes Calculates CO2 and water vapor mole fractions relative to dry air.
LI-7500/A/RS CO2 and H2O 1.15 CO2 on H2O
H2O on CO2
No  
LI-8100/A CO2 and H2O 1.15 CO2 on H2O
H2O on CO2
Yes Calculates the CO2 mole fraction relative to dry air for use in flux calculations.

References

1 Burch, D.E. and D. Williams. 1964. Test of theoretical band model approximations. Applied Optics 3:55-61
2 Burch, D.E., E.B. Singleton and D. Williams. 1962. Absorption line broadening in the infrared. Applied Optics 1:359-363
3 McDermitt, D.K., J.M. Welles and R.D. Eckles. 1993. Effects of temperature, pressure and water vapor on gas phase infrared absorption by CO2. LI-COR, Inc.
4 Parkinson, K.J. 1971. Carbon dioxide infra-red gas analysis. Journal of Experimental Botany 22:169-176
5 Penning de Vries, F. W. T., J. W. J. Akkersdijk and J. L. P. van Oorschot. 1984. An error in measuring respiration and photosynthesis due to transpiration. Photosynthetica 18:146-149
6 Rothman, W.C., C.P. Rinsland, A. Goldman, S.T. Massie, D.P. Edwards, J.M. Flaud. A. Perrin, C. Camy-Peyret, V. Dana, J.Y. Mandin, J. Schroeder, A. McCann, R.R. Gamache, R.B. Wattson, K. Yoshino, K.V. Chance, K.W. Jucks, L.R. Brown, V. Nemtchinov, and P. Varanasi. 1998. The HITRAN molecular spectroscopic database and HAWKS (HITRAN Atmospheric Workstation). 1996 ed. Journal of Quantitative Spectroscopy & Radiative Transfer 60:665-710
7 Webb, E.K., G.I. Pearman, and R. Leuning. 1980. Correction of flux measurements for density effects due to heat and water vapor transfer. Quarterly Journal of the Royal Meteorological Society 106:85-100
8 Welles, J.M. and D.K. McDermitt. 2005. Measuring carbon dioxide in the atmosphere. Micrometeorology in agricultural systems, Agronomy Monograph no. 47:287-320