The fundamental equation of eddy covariance and its application in flux measurements

A fundamental equation of eddy covariance (FQEC) is derived that allows the net ecosystem exchange (NEE) Ns of a specified atmospheric constituent s to be measured with the constraint of conservation of any other atmospheric constituent (e.g. N₂, argon, or dry air). It is shown that if the condition Nsχs NCO2 is true, the conservation of mass can be applied with the assumption of no net ecosystem source or sink of dry air and the FQEC is reduced to the following equation and its approximation for horizontally homogeneous mass fluxes: Here w is vertical velocity, c molar density, t time, h eddy flux measurement height, z vertical distance and χ_s=c_s/c_d molar mixing ratio relative to dry air. Subscripts s, d and CO₂ are for the specified constituent, dry air and carbon dioxide, respectively. Primes and overbars refer to turbulent fluctuations and time averages, respectively. This equation and its approximation are derived for non-steady state conditions that build on the steady-state theory of Webb, Pearman and Leuning (WPL; Webb et al., 1980. Quart. J. R. Meteorol. Soc. 106, 85-100), theory that is widely used to calculate the eddy fluxes of CO₂ and other trace gases. The original WPL constraint of no vertical flux of dry air across the EC measurement plane, which is valid only for steady-state conditions, is replaced with the requirement of no net ecosystem source or sink of dry air for non-steady state conditions. This replacement does not affect the 'eddy flux' term cd w'χ's but requires the change in storage to be calculated as the 'effective change in storage' as follows: Without doing so, significant diurnal and seasonal biases may occur. We demonstrate that the effective change in storage can be estimated accurately with a properly designed profile of mixing ratio measurements made at multiple heights. However further simplification by using a single measurement at the EC instrumentation height is shown to produce substantial biases. It is emphasized that an adequately designed profile system for measuring the effective change in storage in proper units is as important as the eddy flux term for determining NEE. When the EC instrumentation measures densities rather than mixing ratios, it is necessary to use:Ns w'c's h+χs w'c'v +c̄w'T' T̄h+cd (h)∫0hχs t dz. Here T is temperature and c_v and c are the molar densities of water vapor and moist air, respectively. For some atmospheric gas species such as N₂ and O₂, the condition Ns χs NCO2 is not satisfied and additional information is needed in order to apply the EC technique with the constraint of conservation of dry air.