## U.S. Department of Energy

#### Date of this Version

2012

#### Citation

Agricultural and Forest Meteorology 152 (2012) 135– 148;

doi:10.1016/j.agrformet.2011.09.014

#### Abstract

A fundamental equation of eddy covariance (FQEC) is derived that allows the net ecosystem exchange (NEE) *N̅s* of a specified atmospheric constituent s to be measured with the constraint of conservation of any other atmospheric constituent (e.g. N2, argon, or dry air). It is shown that if the condition │*N̅s│* ˃˃* │X̅s│ │N̅co2*│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:

*N̅ _{s} = c̅_{d}w’X’_{s}│h + ∫^{h}_{0} c̅d(z) ∂X_{s}/*

*∂t dz +*

*∫*

^{h}_{0}[X*̅*

_{s}(z)- X*̅*

_{s}(h)]*∂̅c*

*̅*

_{d}*̅/∂t dz = c̅*

_{d̅}(h) {w̅’X̅’_{s}│h + ∫^{h}_{0}∂X_{s}/*∂t dz}*.

Here *w* is vertical velocity, *c* molar density, *t* time, *h* eddy flux measurement height, *z* vertical distance and *Xs= cs/cd* molar mixing ratio relative to dry air. Subscripts *s, d* and CO2 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 CO2 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 c̅d̅w̅’X̅’s s but requires the change in storage to be calculated as the ‘*effective change in storage*’ as follows:

*∫ ^{h}_{0 }*

*∂̅c*

*̅*

_{s}*̅/*

*∂̅t*

*̅ dz – X*

*̅*

_{s}(h)*∫*

^{h}_{0}∂̅c̅_{d̅}/∂t dz = ∫^{h}_{0}c̅_{d̅}(z) - ∂Xs/*∂t dz +*

*∫*

^{h}_{0}[X*̅*

_{s}(z)- X

_{̅s}(h)]*∂̅c*

*̅*

_{d}*̅/∂t dz= c̅d (h) ∫*

^{h}_{0}∂X_{s}/*∂t dz.*

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.