Dry-air boundary conditions for correction of eddy flux measurements

For measurements of eddy fluxes in the atmospheric boundary layer of gases (such as CO₂) whose average concentration is very large compared to the fluctuations, corrections for air density fluctuations are required. With the boundary condition of no flux of dry air at the surface, the evaporation correction to eddy fluxes is 2.6 times larger than has been estimated with the boundary condition of no mass flux at all at the surface. The heat flux correction is also increased by a few per cent.     

Eddy-covariance CO2 flux measurements using open- and closed-path CO2 analysers: Corrections for analyser water vapour sensitivity and damping of fluctuations in air sampling tubes

Methods of calibrating infrared CO₂ analysers for sensitivity to CO₂ and water vapour are described. Equations to correct eddy covariance CO₂ flux measurements are presented for: (i) analyser cross-sensitivity to water vapour and the effects of density fluctuations arising from atmospheric fluxes of water vapour and sensible heat, (ii) flux losses caused by signal processing and limited instrument frequency response for open- and closed-path CO₂ analysers, and (iii) flux losses resulting from damping of concentration fluctuations in a tube used to sample air for closed-path CO₂ analysers. Examples of flux corrections required for typical instruments are presented.     

Correction Of Eddy-Covariance Measurements Incorporating Both Advective Effects And Density Fluxes

Equations are presented to correct eddy-covariancemeasurements for both fluctuations in density andnon-zero mean advection, induced by convergence ordivergence of flow, and spatial source/sinkinhomogeneity, under steady-state and transientconditions. This correction collapses to theWebb–Pearman–Leuning expression ifthe mean vertical velocity is zero, and formally addsthe Webb–Pearman–Leuning expression to the correctionssuggested by Lee for conditions ofnon-zero vertical velocity and source/sink and meanscalar horizontal homogeneity. The equation requiresmeasurement of the mean vertical gradients of thescalar concentration of interest (air temperature,humidity, CO₂) as well as an accurateestimation of the mean vertical velocity, in additionto the vertical eddy covariance of the scalar. Simplemethods for the approximation of sensor tilt andcomplex terrain flow angle are presented, to allowestimation of non-zero mean vertical velocities. Theequations are applied to data from a maize crop and aforest to give examples of when the correction issignificant. In addition, a term for thethermodynamic expansion energy associated with watervapour flux is derived, which implies that the sonictemperature derived sensible heat flux will accuratelyinclude this contribution.     

Eddy correlation measurements of methane flux in a northern peatland ecosystem

A pilot study to measure methane flux using eddy correlation sensors was conducted in a peatland ecosystem in north central Minnesota. A prototype tunable diode laser spectrometer system was employed to measure the fluctuations in methane concentration.The logarithmic cospectrum of methane concentration and vertical wind velocity fluctuations under moderately unstable conditions had a peak nearf = 0.10 (wheref is the nondimensional frequency) and was quite similar to the cospectra of water vapor and sensible heat. Daytime methane flux during the first two weeks of August ranged from 120 to 270 mg m^-2 day^-1. The temporal variation in methane fluxes was consistent with changes in peat temperature and water table elevation. Our results compared well with the range of values obtained in previous studies in Minnesota peatlands.These field observations demonstrate the utility of the micrometeorological eddy correlation technique for measuring surface fluxes of methane. The current state-of-the-art in tunable diode laser spectroscopy makes this approach practical for use in key ecosystems.Published as Paper No. 9556, Journal Series, Nebraska Agricultural Research Division.     

Comparison of eddy-covariance measurements of CO2 fluxes by open- and closed-path CO2 analysers

Eddy fluxes of CO₂ estimated using a sonic anemometer and a closed-path analyser were, on average, 16% lower than those obtained with the same anemometer and an adjacent open-path CO₂ analyser. Covariances between vertical windspeed and CO₂ density from the closed-path analyser were calculated using data points for CO₂ that were delayed relative to anemometer data by the time required for a parcel of air to travel from the tube inlet to the CO₂ sensor. Air flow in the intake tube was laminar. Densities of CO₂ that had been corrected for spurious fluctuations arising from fluctuations in temperature and humidity were used in the flux calculations. Corrections for the cross-sensitivity of CO₂ analysers to water vapour were also incorporated. Spectral analysis of the corrected CO₂ signal from the closed-path analyser showed that damping of fluctuations in the sampling tube at frequencies f > 0.1 Hz caused the apparent loss in flux. The measured losses can be predicted accurately using theory that describes the damping of oscillations in a sampling tube. High-frequency response of the closed-path system can be improved substantially by ensuring turbulent flow in the tube, using a combination of high volumetric flow rate and small tube diameter. The analysis of attenuation of turbulent fluctuations in flow through tubes is applicable to the measurement of fluxes of other minor atmospheric constituents using the eddy covariance method.     

Limitations of an eddy-correlation technique for the determination of the carbon dioxide and sensible heat fluxes

A modified infrared CO₂ gas analyzer, a small thermocouple assembly, a heated-thermocouple anemometer for horizontal wind, and a propeller-type vertical wind sensor were used to measure the eddy fluxes of heat and CO₂ above a corn crop. Experimental results of these fluxes are discussed. The main sources of errors of the eddy fluxes using these instruments were estimated:

1. Sensors with a time constant of 0.5 s appear to be fast enough to detect most of the vertical CO₂ transfer as long as the sensors are located at least one meter above the crop surface.
2. The deviation from steady-state conditions for 10-min periods was found to have a significant effect on the eddy flux estimates.
3. Temperature fluctuations of the air sample passing through the CO₂ infrared gas analyzer were found to be non-negligible but could be easily corrected.
4. A 1° misalignment of the vertical anemometer affected these eddy fluxes by less than 10% under all circumstances studied.

     

A novel ozone sensor for direct eddy flux measurements

A small, lightweight (1.5 kg) and fast-response ozone sensor for direct eddy flux measurements has been built. The basis for detection is the chemiluminescence of an organic dye adsorbed on dry silica gel in the reaction with ozone. The chemiluminescence is monitored with a cheap and small blue-sensitive photomultiplier. At a flow rate of 100 l min^-1 the ozone sensor has a 90% response time of significantly better than 0.1 s with a detection limit lower than 50 ppt at S/N=3. There are no interferences from other atmospheric trace gases like NO_x, H₂O₂ and PAN. Water vapour and SO₂ enhance the chemiluminescence efficiency of the ozone sensor. Since their response times are 22 seconds and 30 minutes, respectively, no correlation between rapid ozone fluctuations and those of these two trace gases is noticed by the ozone sensor when operating at a frequency of 10 Hz.The ozone sensor was tested for several weeks in continuous measurements of ozone fluxes and deposition velocities over different croplands using the eddy correlation technique. Good agreement was found between ozone dry deposition velocities derived from profile measurements and by eddy correlation.     

Estimates of sensible heat flux from observations of temperature fluctuations

Comparisons between sensible heat flux measured using eddy correlation instrumentation and estimated using the temperature fluctuation method are presented for four types of surface in West Africa. Agreement between measured and estimated values is good. Regression of estimated on measured sensible heat flux gave a mean slope of 0.98 with a mean r ² of 0.94 for bare soil, mature millet, fallow savannah and tiger bush. Estimates of heat flux from temperature fluctuations measured by an instrument mounted beneath a tethered balloon are also shown to be in close agreement with eddy correlation measurements made at the surface (regression slope = 0.98, r ² = 0.84). The results provide evidence that the ratio /_×is indeed a universal function of z/L for all the surface types considered.     

On dynamically consistent eddy fluxes

The role of mesoscale oceanic eddies in driving the large-scale currents is studied in an eddy-resolving, double-gyre ocean model. The new diagnostic method is proposed, which is based on dynamical decomposition of the flow into the large-scale and eddy components. The method yields the time history of the eddy forcing, which can be used as additional, external forcing in the corresponding non-eddy-resolving model of the gyres. The main strength of this approach is in its dynamical consistency: the non-eddy-resolving solution driven by the eddy forcing history correctly approximates the original large-scale flow component. It is shown that statistical decompositions, which are based on space-time filtering diagnostics, are dynamically inconsistent. The diagnostics algorithm is formulated and tested, and the diagnosed eddies are analysed, both statistically and dynamically. It is argued that the main dynamic role of the eddies is to maintain the eastward-jet extension of the subtropical western boundary current (WBC). This is done largely by both the time–mean isopycnal-thickness flux and the relative-vorticity eddy flux fluctuations. The fluctuations drive large-scale flow through the nonlinear rectification mechanism. The relative-vorticity flux contributes mostly to the eastward jet meandering. Finally, eddy fluxes driven by both the eddies and the large-scale flow are found to be important. The latter is typically neglected in the analysis, but here it corresponds to important large-scale feedback on the eddies.     

Eddy correlation measurement of CO2 flux using a closed-path sensor: Theory and field tests against an open-path sensor

We have examined the potential of using a closed-path sensor to accurately measure eddy fluxes of CO₂. Five inlet tubeflow configurations were employed in the experimental setup. The fluxes of CO₂ were compared against those measured with an open-path sensor. Sampling air through an intake tube causes a loss of flux, due to the attenuation of CO₂ density fluctuations. Adjustments need to be made to correct for this loss and to account for density effects due to the simultaneous transfer of heat and water vapor. Theory quantifying these effects is discussed.The raw CO₂ flux measured with the closed-path sensor was smaller than that measured with the open-path sensor by about 15% (on average) for the turbulent tubeflow configurations with a short (3 m) intake tube, by 31% for turbulent tubeflow with a longer (6 m) intake tube and by 24% for laminar tubeflow. The difference was, in part, caused by tube attenuation of the CO₂ density fluctuations and inadequate sensor time response. The elimination of the flux adjustment for the simultaneous transfer of sensible heat (i.e., the attenuation of ambient temperature fluctuations in the intake tube) generally accounted for the rest of this difference.The raw flux measured with the closed-path sensor was corrected for frequency response and density effects. Except in the case of laminar tubeflow, the corrected closed-path flux agreed consistently with the corrected open-path flux within a few percent (<5%). These results suggest that closed-path sensors, with appropriate corrections, can be used to measure CO₂ flux accurately. Recommendations are included on selecting an optimum flow configuration to minimize the effect of sampling air through a tube.Published as Paper No. 9938, Journal Series, Nebraska Agricultural Research Division.     

北京城区二氧化碳时空分布及湍流谱特征

利用北京325 m气象塔上安装的7层CO₂涡动相关系统在2014年12月到2015年11月的观测资料,分析了北京城区不同高度上CO₂浓度、通量时空分布及湍流谱的特征。结果表明:城市CO₂浓度日变化除了冬季都呈现双峰型,冬季由于人为碳源排放的大幅增加,双峰型不明显。每层的CO₂浓度、通量都有明显的季节变化:冬季最高,春末、夏季最低。CO₂浓度整体随高度的增加而降低。北京城区是CO₂源,CO₂通量的日变化不如CO₂浓度日变化规律明显。CO₂通量在47 m以下为负,47 m以上为正。通量在140 m以下随高度的增加而增加;140m以上随高度的增加而减少。根据对CO₂时空分布的分析可知:边界层CO₂浓度、通量强烈受到碳源、下垫面植被、大气稳定度、环境温度和天气过程等因素的影响。各变量谱与Kaimal等的研究结果接近:归一化速度谱和CO₂谱在惯性子区有－2/3的斜率,在低频区与稳定度参数（Z/L）有一定的关系。这说明复杂地形的城市下垫面的湍流谱结构与平坦地形相比没有太大的实质性差异。     

Estimation of the Exchange Coefficient of Heat During Low Wind Convective Conditions

For the first time, the exchange coefficient of heat C_H has been estimated from eddy correlation of velocity and virtual temperature fluctuations using sonic anemometer measurements made at low wind speeds over the monsoon land atJodhpur (26°18' N, 73°04' E), a semi arid station. It shows strong dependence on wind speed, increasing rapidly with decreasing wind speed, and scales according to a power law C_H = 0.025U₁₀ ^-0.7 (where U₁₀ is the mean wind speed at 10-m height). A similar but more rapid increase in the drag coefficient C_Dhas already been reported in an earlier study. Low winds (<4 m s^-1) are associated with both near neutral and strong unstable situations. It is noted that C_H increases with increasing instability. The present observations best describe a low wind convective regime as revealed in the scaling behaviour of drag, sensible heat flux and the non-dimensional temperature gradient. Neutral drag and heat cofficients,corrected using Monin–Obukhov (M–O) theory, show a more uniform behaviour at low wind speeds in convective conditions, when compared with the observed coefficients discussed in a coming paper.At low wind convective conditions, M-O theory is unable to capture the observed linear dependence of drag on wind speed, unlike during forced convections. The non-dimensional shear inferred from the present data shows noticeable deviations from Businger's formulation, a forced convection similarity. Heat flux is insensitive to drag associated with weak winds superposed on true free convection. With heat flux as the primary variable, definition of new velocity scales leads to a new drag parameterization scheme at low wind speeds during convective conditionsdiscussed in a coming paper.     

Automated in-situ monitoring of atmospheric non-methane hydrocarbon concentrations and gradients

A fully automated system measuring C2–C6 hydrocarbon concentrations and vertical gradients was installed at Harvard Forest in Petersham, Massachusetts, using a gas chromatograph with dual flame ionization detectors and cryogenic sample preconcentration. Measurements were made simultaneously at two heights above the forest canopy at forty five minute intervals, continuously from July 1992 to the present. Data for concentration gradients were combined with CO₂ flux measured by eddy correlation to determine the rates of production of biogenic hydrocarbons by the forest.     

Turbulent carbon exchange in very stable conditions

Turbulent fluxes obtained using the conventional eddy covariance approach result in erratic results with large time fluctuations in extremely stable conditions. This can limit efforts to estimate components of the nocturnal energy budget and respiratory CO₂ fluxes. Well-organized fluxes that show a clear dependence on turbulent intensity were obtained when multiresolution decomposition was used to estimate turbulent exchanges. CO₂, heat and water vapour fluxes were observed at a site in the eastern Amazon basin that had been cleared for agricultural purposes. Temporal scales of the carbon transfer were determined and shown to be similar to those of latent heat, but as much as three times larger than those of sensible heat. CO₂ eddy diffusivities at the temporal scales on which most of the vertical CO₂ exchange occurs are shown to be 50 times larger than the eddy diffusivity for heat. A process associated with the vertical scale of the scalar accumulation layer is suggested to explain these different scales and turbulent diffusivities of carbon and sensible heat transfer. For an appreciable range of turbulence intensities, the observed vertical turbulent carbon exchange is insufficient to account for the locally respired CO₂ estimated independently. Evidence that shallow drainage currents may account for this is given.     

高频瞬变涡动反馈强迫对东亚/太平洋事件演变过程的作用

本文通过位势倾向方程研究了高频瞬变涡动反馈强迫(TEFF)对东亚/太平洋型在低频时间尺度上(8天以上)的演变过程(EAP事件)的影响。一个重要的发现是,无论正位相还是负位相的EAP事件,在对流层中低层上,TEFF对其中、高纬度两个异常中心的演变均起着重要的作用。尤其是高纬度中心,在对流层中层,TEFF异常在盛期前可解释其增强部分的50%左右;盛期之后则起着减缓其减弱的作用。TEFF异常可分解为由高频涡动热量通量的辐合辐散引起的TEFF_heat异常和高频涡动涡度通量引起的TEFF_vor异常。两者在对流层上层相互抵消、下层相互叠加,导致TEFF异常在对流层上层(低层)的作用要弱(强)于对流层中层。就TEFF_heat和TEFF_vor的相对贡献而言,正位相的EAP事件中,TEFF_vor的贡献较大;但在负位相EAP事件中,TEFF_heat的贡献则增大至与TEFF_vor相当的程度。这种TEFF_heat作用在负事件中的增强可能是与贝加尔湖以东地区对流层低层上斜压性的加强有关。     

On the relationship between the eddy covariance,the turbulent flux,and surface exchange for a trace gas such as CO<Subscript>2</Subscript>

In the context of CO₂ surface exchange estimation, an analysis combining the basic principles of diffusion and scalar conservation shows that the mixing ratio is the appropriate variable both for defining the (eddy covariance) turbulent flux and also for expressing the relationship between the turbulent flux and surface exchange in boundary-layer budget equations. Other scalar intensity variables sometimes chosen, both the CO₂ density and molar fraction, are susceptible to the influence of surface exchange of heat and water vapour. The application of a hypsometric analysis to the boundary-layer “control volume” below the tower measurement height reveals flaws in previously applied approaches: (a) incompressibility cannot be assumed to simplify mass conservation (the budget in terms of CO₂ density); (b) compressibility alone makes the analysis of mass conservation vulnerable to uncertainties associated with resultant non-zero vertical velocities too small to measure or model over real terrain; and (c) the WPL (Webb et al. (1980) Quart J Roy Meteorol Soc 106:85–100) “zero dry air flux” assumption is invalidated except at the surface boundary. Nevertheless, the definition and removal of the WPL terms do not hinge upon this last assumption, and so the turbulent CO₂ flux can be accurately determined by eddy covariance using gas analysers of either open- or closed-path design. An appendix discusses the necessary assumptions and appropriate interpretations for deriving the WPL terms.     

Unsteadiness as a cause of non-equality of eddy diffusivities for heat and vapour at the base of an advective inversion

This paper examines the effect of non-stationarity of the wind on similarity of the eddy diffusivities for heat and vapour within a stable layer at the bottom of an internal boundary layer formed downwind of a dry-to-wet transition. First, we present some experimental data taken above a rice crop downwind of very extensive dry range lands at Warrawidgee, NSW, Australia. These data establish that periods of higher wind speed were associated with periods of higher saturation deficit in the canopy of the rice crop, and lower Bowen ratio. It is shown that Bowen ratios calculated for 30-second sub-intervals varied three-fold within a single 20-minute averaging period. Thus periods of higher wind speed corresponded to periods of higher moisture flux and smaller sensible heat flux.An idealized situation is then analysed theoretically. It is assumed that the time scale of the slow variations of the wind is long compared with the surface-layer time scale and that fetch is sufficient that the air near the ground is in continuous equilibrium with the surface. Using a two-scale Reynolds decomposition of the fluctuating wind and scalar variables into active and inactive components, it is shown that unsteadiness can lead to an eddy diffusivity for saturation deficit, calculated as the ratio of average flux to average gradient, that is larger than that for total energy calculated in a similar way. Using this ratio to calculate the ratio of diffusivities for temperature and humidity, K_T/K_q, it is found that the latter can be much larger than one if the Bowen ratio is small and negative. Despite this, assuming K_T = K_q and using the Bowen ratio method to calculate surface energy fluxes will usually incur only minor errors.     

A Simple Method of Estimating Scalar Fluxes Over Forests

A simple aerodynamic-variance method is proposed to fill gaps in continuous CO₂ flux measurements in rainy conditions, when open-path analysers do not function. The method requires turbulent conditions (friction velocity greater than 0.1 ms^–1), and uses measurements of mean wind speed, and standard deviations of temperature and CO₂ concentration fluctuations to complement, and at times replace, eddy-covariance measurements of friction velocity, sensible heat flux and CO₂ flux. Friction velocity is estimated from the mean wind speed with a flux-gradient relationship modified for the roughness sublayer. Since normalised standard deviations do not follow Monin-Obukhov similarity theory in the roughness sublayer, a simple classification scheme according to the scalar turbulence scale was used. This scheme is shown to produce sensible heat and CO₂ flux estimates that are well correlated with the measured values.     

A Simple Model for the Afternoon and Early Evening Decay of Convective Turbulence Over Different Land Surfaces

A simple model to study the decay of turbulent kinetic energy (TKE) in the convective surface layer is presented. In this model, the TKE is dependent upon two terms, the turbulent dissipation rate and the surface buoyancy fluctuations. The time evolution of the surface sensible heat flux is modelled based on fitting functions of actual measurements from the LITFASS-2003 field campaign. These fitting functions carry an amplitude and a time scale. With this approach, the sensible heat flux can be estimated without having to solve the entire surface energy balance. The period of interest covers two characteristic transition sub-periods involved in the decay of convective boundary-layer turbulence. The first sub-period is the afternoon transition, when the sensible heat flux starts to decrease in response to the reduction in solar radiation. It is typically associated with a decay rate of TKE of approximately t ⁻² (t is time following the start of the decay) after several convective eddy turnover times. The early evening transition is the second sub-period, typically just before sunset when the surface sensible heat flux becomes negative. This sub-period is characterized by an abrupt decay in TKE associated with the rapid collapse of turbulence. Overall, the results presented show a significant improvement of the modelled TKE decay when compared to the often applied assumption of a sensible heat flux decreasing instantaneously or with a very short forcing time scale. In addition, for atmospheric modelling studies, it is suggested that the afternoon and early evening decay of sensible heat flux be modelled as a complementary error function.     

Comparative Measurements Of Co2flux Over A Forest Using Closed-Path And Open-Path Co2analysers

An eddy covariance system using a closed-path CO₂analyser was constructed for long-term CO₂flux measurements above a forest, and its total frequency response was valuated experimentally. The amplitude and phase responses of the system wereexamined through a preliminary test, in which a prescribed pattern of CO₂fluctuation was input to the system. The result showed that the amplitude of the output from the system was attenuated as frequency increased, with a half-power frequency of 0.3 Hz. The phase was delayed by the air sampling through a long tube, but the delay in phase decreased asfrequency increased. We then presented a new technique for the correction of flux loss due to the inadequate system response for the eddy covariance measurements of CO₂ flux. Using the present system and the correction technique, diurnal variations in CO₂ flux were measured over a temperate deciduous forest on three days in 1997. The results were compared with the CO₂fluxes measured with a fast response open-path gas analyser. The CO₂fluxes from the closed-path system agreed with those from the open-path system after the Webb, Pearman and Leuning correction was made for the latter. In the present test over a forest, the contribution of the frequency-response correction to the CO₂fluxes was small and its averaged percentage was only 3%in the daytime. However, the percentage would likely increase, if the system were applied to a shorter vegetation site where high frequency components are more important. The comparison confirmed that we can obtain correct measurements of CO₂flux using the present closed-path system and the correction technique.