A Note on the Flux-Variance Similarity Relationships for Heat and Water Vapour in the Unstable Atmospheric Surface Layer

Atmospheric surface layer (ASL) experiments over the past 10 years demonstrate that the flux-variance similarity functions for water vapour are consistently larger in magnitude than their temperature counterpart. In addition, latent heat flux calculations using the flux-variance method do not compare as favorably to eddy- correlation measurements when compared to their sensible heat counterpart. These two findings, in concert with measured heat to water vapour transport efficiencies in excess of unity, are commonly used as evidence of dissimilarity between heat and water vapour transport in the unstable atmospheric surface layer. In this note, it is demonstrated that even if near equality in flux-profile similarity functions for heat and water vapour is satisfied, the flux-variance similarity functions for water vapour are larger in magnitude than temperature for a planar, homogeneous, unstably-stratified turbulent boundary-layer flow.     

Heat and water vapour fluxes and scalar roughness lengths over an Antarctic ice shelf

We present eddy-correlation measurements of heat and water vapour fluxes made during the Antarctic winter. The surface layer was stably stratified throughout the period of observation and sensible heat fluxes were always directed downwards. However, both upward and downward water vapour fluxes were observed. Their magnitude was generally small and the latent heat flux was not a significant fraction of the surface energy budget. The variation of heat and water vapour fluxes with stability is well described by Monin-Obukhov similarity theory but the scalar roughness lengths for heat and water vapour appear to be much larger than the momentum roughness length. Possible explanations of this effect are discussed.     

Testing of the bandpass eddy covariance method for a long-term measurement of water vapour flux over a forest

The bandpass eddy covariance method has been used to measure the turbulent flux of scalar quantities using a slow-responsescalar sensor. The method issimilar in principle to the traditional eddy correlation method but includes the estimation of high-frequency components of the flux on the basis of cospectral similarity in the atmospheric surface layer. In order to investigate the performance of the method, measurements of the water vapour flux over a forest with the bandpass eddy covariance method and the direct eddy correlation method were compared. The flux obtained by the bandpass eddy covariance method agreed with that by the eddy correlation method within ±20% for most cases, in spite of a rather slow sensor-response of the adopted hygrometer. This result supports its relevance to a long-term continuous operation, since a stable, low-maintenance,general-purpose sensor canbe utilized for scalar quantities. Oneweak point of the method isits difficulty in principle to measure the correct flux when the magnitude of the sensible heat flux is very small, because the method uses the sensible heat flux as a standard reference for the prediction of undetectable high-frequency components of the scalar flux. An advanced method is then presented to increase its robustness. In the new method, output signals from a slow-response sensor are corrected using empirical frequency-responsefunctions for the sensor,thereby extending the width of the bandpass frequency region where components of the flux are directly measured (not predicted). The advanced method produced correct fluxes for all cases including the cases of small sensible heat flux. The advanced bandpass eddy covariance method is thus appropriate for along-term measurement of the scalar fluxes.     

Coherent eddies and temperature structure functions for three contrasting surfaces. Part II: Renewal model for sensible heat flux

Sensible heat, latent heat, and other scalar fluxes cannot be measuredwithin short dense canopies, e.g., straw mulches, with standard approachessuch as eddy correlation, Bowen ratio-energy balance, aerodynamic, andvariance methods. However, recently developed surface renewal models, thatare based on the fact that most of the turbulent transfer within and abovecanopies is associated with large-scale coherent eddies, which are evidentas ramp patterns in scalar time series, offer a feasible solution. Wepresent a new air renewal model that calculates sensible heat flux atdifferent heights within and above a canopy from the average cubictemperature structure function, sampled at a moderate rate, and measuredaverage friction velocity. The model is calibrated and tested with datameasured above and within a Douglas-fir forest and above a straw mulch andbare soil. We show that the model describes half-hour variations ofsensible heat flux very well, both within the canopy and roughnesssublayers and in the inertial sublayer, for stable and unstable atmosphericconditions. The combined empirical coefficient that appears in the modelhas an apparently universal value of about 0.4 for all surfaces andheights, which makes application of the model particularly simple. Themodel is used to predict daytime and nighttime sensible heat flux profileswithin the straw mulch and within a small bare opening in the mulch.     

Estimation of surface heat and momentum fluxes using the flux-variance method above uniform and non-uniform terrain

Eddy-correlation measurements above an uneven-aged forest, a uniform-irrigated bare soil field, and within a grass-covered forest clearing were used to investigate the usefulness of the fluxvariance method above uniform and non-uniform terrain. For this purpose, the Monin and Obukhov (1954) variance similarity functions were compared with direct measurements. Such comparisons were in close agreement for momentum and heat but not for water vapor. Deviations between measured and predicted similarity functions for water vapor were attributed to three factors: 1) the active role of temperature in surface-layer turbulence, 2) dissimilarity between sources and sinks of heat and water vapor at the ground surface, and 3) the non-uniformity in water vapor sources and sinks. It was demonstrated that the latter non-uniformity contributed to horizontal gradients that do not scale with the vertical flux. These three factors resulted in a turbulence regime that appeared more efficient in transporting heat than water vapor for the dynamic convective sublayer but not for the dynamic sublayer. The agreement between eddy-correlation measured and flux-variance predicted sensible heat flux was better than that for latent heat flux at all three sites. The flux-variance method systematically overestimated the latent heat flux when compared to eddy-correlation measurements. It was demonstrated that the non-uniformity in water vapor sources reduced the surface flux when compared to an equivalent uniform terrain subjected to identical shear stress, sensible heat flux, and atmospheric water vapor variance. Finally, the correlation between the temperature and water vapor fluctuations was related to the relative efficiency of surface-layer turbulence in removing heat and water vapor. These relations were used to assess critical assumptions in the derivation of the flux-variance formulation.     

Temperature and humidity flux-variance relations determined by one-dimensional eddy correlation

It may be possible to estimate surface fluxes of scalar quantities from measurement of their variance and mean wind speed. The flux-variance relation for temperature and humidity was investigated over prairie and desert-shrub plant communities. Fluxes were measured by one-dimensional eddy correlation, humidity by fast-response wet-bulb psychrometers and Krypton open-path hygrometers, temperature by fine-wire thermocouples, and mean windspeed by a cup anemometer. The quality of the flux-variance relation proved to be good enough for application to flux measurement. Regressions of flux estimated by the variance technique versus measured flux usually had r ² values greater than 0.97 for sensible heat flux and greater than 0.88 for water vapor flux. More uniform surfaces tended to yield the same flux-variance relations except when fluxes were small. This exception supported the hypothesis that sparse sources of flux may increase variance downwind. Nonuniform surfaces yielded flux-variance relations that were less predictable, although reasonably accurate once determined. The flux-variance relation for humidity was quite variable over dry surfaces with senescent vegetation.     

基于通量方差法对青藏高原东侧高寒草甸地表通量的估算

利用位于青藏高原东侧理塘大气综合观测站2008年观测资料,分析了高寒草甸下垫面上地表通量的时间变化特征,确定了温度、水汽和CO₂的归一化标准差在不稳定情况下随稳定度变化的通量方差关系,应用通量方差法对感热、潜热和CO₂通量进行了计算,并与涡旋相关系统的观测结果进行了比较。结果表明：地表通量月平均日变化呈较为规则的日循环特征,季节变化特征也很明显,雨季（5-9月）潜热大于感热,干季则以感热为主,CO₂通量以6-9月值最大。在不稳定条件下,温度、水汽和CO₂的归一化标准差随稳定度的变化均满足-1/3规律,其通量方差相似性常数分别为1.2,1.4和0.9。通量方差法估算出的通量值与涡旋相关观测得到的通量值有较好的一致性,但感热通量的效果优于潜热通量和CO₂通量。该方法高估了感热通量尤其是潜热通量,而低估了CO₂通量。采用直接观测的感热通量值计算潜热通量和CO₂通量可改善计算结果。     

A lagrangian random-walk model for simulating water vapor,CO2 and sensible heat flux densities and scalar profiles over and within a soybean canopy

An integrated canopy micrometeorological model is described for calculating CO₂, water vapor and sensible heat exchange rates and scalar concentration profiles over and within a crop canopy. The integrated model employs a Lagrangian random walk algorithm to calculate turbulent diffusion. The integrated model extends previous Lagrangian modelling efforts by employing biochemical, physiological and micrometeorological principles to evaluate vegetative sources and sinks. Model simulations of water vapor, CO₂ and sensible heat flux densities are tested against measurements made over a soybean canopy, while calculations of scalar profiles are tested against measurements made above and within the canopy. The model simulates energy and mass fluxes and scalar profiles above the canopy successfully. On the other hand, model calculations of scalar profiles inside the canopy do not match measurements.The tested Lagrangian model is also used to evaluate simpler modelling schemes, as needed for regional and global applications. Simple, half-order closure modelling schemes (which assume a constant scalar profile in the canopy) do not yield large errors in the computation of latent heat (LE) and CO₂ (F _c ) flux densities. Small errors occur because the source-sink formulation of LE andF _c are relatively insensitive to changes in scalar concentrations and the scalar gradients are small. On the other hand, complicated modelling frames may be needed to calculate sensible heat flux densities; the source-sink formulation of sensible heat is closely coupled to the within-canopy air temperature profile.     

Monin–Obukhov Similarity Functions for the Structure Parameters of Temperature and Humidity

Monin–Obukhov similarity functions for the structure parameters of temperature and humidity are needed to derive surface heat and water vapour fluxes from scintillometer measurements and it is often assumed that the two functions are identical in the atmospheric surface layer. Nevertheless, this assumption has not yet been verified experimentally. This study investigates the dissimilarity between the turbulent transport of sensible heat and water vapour, with a specific focus on the difference between the Monin–Obukhov similarity functions for the structure parameters. Using two datasets collected over homogeneous surfaces where the surface sources of sensible heat and water vapour are well correlated, we observe that under stable and very unstable conditions, the two functions are similar. This similarity however breaks down under weakly unstable conditions; in that regime, the absolute values of the correlations between temperature and humidity are also observed to be low, most likely due to large-scale eddies that transport unsteadiness, advection or entrainment effects from the outer layer. We analyze and demonstrate how this reduction in the correlation leads to dissimilarity between the turbulent transport of these two scalars and the corresponding Monin–Obukhov similarity functions for their structure parameters. A model to derive sensible and latent heat fluxes from structure parameters without measuring the friction velocity is tested and found to work very well under moderately to strongly unstable conditions (−z/L > 0.5). Finally, we discuss the modelling of the cross-structure parameter over wet surfaces, which is crucial for correcting water vapour effects on optical scintillometer measurements and also for obtaining surface sensible and latent heat fluxes from the two-wavelength scintillometry.     

Intra-field variability of scalar flux densities across a transition between a desert and an irrigated potato field

This paper reports on measurements of sensible and latent heat and CO₂ fluxes made over an irrigated potato field, growing next to a patch of desert. The study was conducted using two eddy correlation systems. One measurement system was located within the equilibrium boundary layer 800 m downwind from the edge of the potato field. The other measurement system was mobile and was placed at various downwind positions to probe the horizontal transition of vertical scalar fluxes. Latent (LE) and sensible (H) heat fluxes, measured at 4 m above the surface, exhibited marked variations with downwind distance over the field. Only after the fetch to height ratio exceeded 75 to 1 didLE andH become invariant with downwind distance. When latent and sensible heat fluxes were measured upwind of this threshold, significant advection of humidity-deficit occurred, causing a vertical flux divergence ofH andLE.The measured fluxes of momentum, heat, and moisture were compared with predictions from a second-order closure two-dimensional atmospheric boundary layer model. There is good agreement between measurements and model predictions. A soil-plant-atmosphere model was used to examine nonlinear feedbacks between humidity-deficits, stomatal conductance and evaporation. Data interpretation with this model revealed that the advection of hot dry air did not enhance surface evaporation rates near the upwind edge of the potato field, because of negative feedbacks among stomatal conductance, humidity-deficits, andLE. This finding is consistent with results from several recent studies.     

应用通量方差法估算戈壁绿洲下垫面湍流通量的研究

利用“绿洲系统能量与水分循环过程观测试验” 2005年绿洲、戈壁点的观测资料, 分析与讨论了温度、水汽的归一化标准差随稳定度变化的通量方差关系, 应用通量方差法对感热、 潜热通量进行了计算, 并同涡动相关系统的观测结果进行了比较。不稳定条件下, 戈壁点温度归一化标准差随稳定度变化的通量方差关系优于下垫面非均匀性更强的绿洲点, 绿洲点水汽的归一化标准差随稳定度变化的通量方差关系较温度量表现得更好。对同一站点, 归一化温度标准差的通量方差关系并不总是优于水汽的通量方差关系, 其取决于该站点的温度以及水汽的源汇分布情况; 通量方差法对两个站点的感热、 潜热通量均有较好的再现, 但戈壁点感热通量的计算效果优于非均匀性更强的绿洲点。应用通量方差法对潜热通量计算时若采用直接观测的感热通量, 则潜热通量的计算效果具有一定程度的改善。     

Estimation of Turbulent Fluxes Using the Flux-Variance Method over an Alpine Meadow Surface in the Eastern Tibetan Plateau

The flux-variance similarity relation and the vertical transfer of scalars exhibit dissimilarity over different types of surfaces,resulting in different parameterization approaches of relative transport efficiency among scalars to estimate turbulent fluxes using the flux-variance method.We investigated these issues using eddycovariance measurements over an open,homogeneous and flat grassland in the eastern Tibetan Plateau in summer under intermediate hydrological conditions during rainy season.In unstable conditions,the temperature,water vapor,and CO2 followed the flux-variance similarity relation,but did not show in precisely the same way due to different roles(active or passive) of these scalars.Similarity constants of temperature,water vapor and CO2 were found to be 1.12,1.19 and 1.17,respectively.Heat transportation was more efficient than water vapor and CO2.Based on the estimated sensible heat flux,five parameterization methods of relative transport efficiency of heat to water vapor and CO2 were examined to estimate latent heat and CO2 fluxes.The strategy of local determination of flux-variance similarity relation is recommended for the estimation of latent heat and CO2 fluxes.This approach is better for representing the averaged relative transport efficiency,and technically easier to apply,compared to other more complex ones.     

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.     

Dissimilarity of Scalar Transport in the Convective Boundary Layer in Inhomogeneous Landscapes

A land-surface model (LSM) is coupled with a large-eddy simulation (LES) model to investigate the vegetation-atmosphere exchange of heat, water vapour, and carbon dioxide (CO₂) in heterogeneous landscapes. The dissimilarity of scalar transport in the lower convective boundary layer is quantified in several ways: eddy diffusivity, spatial structure of the scalar fields, and spatial and temporal variations in the surface fluxes of these scalars. The results show that eddy diffusivities differ among the three scalars, by up to 10–12%, in the surface layer; the difference is partly attributed to the influence of top-down diffusion. The turbulence-organized structures of CO₂ bear more resemblance to those of water vapour than those of the potential temperature. The surface fluxes when coupled with the flow aloft show large spatial variations even with perfectly homogeneous surface conditions and constant solar radiation forcing across the horizontal simulation domain. In general, the surface sensible heat flux shows the greatest spatial and temporal variations, and the CO₂ flux the least. Furthermore, our results show that the one-dimensional land-surface model scheme underestimates the surface heat flux by 3–8% and overestimates the water vapour and CO₂ fluxes by 2–8% and 1–9%, respectively, as compared to the flux simulated with the coupled LES-LSM.     

Comparison of areally averaged turbulent fluxes over non-homogeneous terrain: Results from the EFEDA-field experiment

In June 1991 the EFEDA-field experiment (ECHIVAL Field Experiment in a Desertification-Threatened Area) was carried out in the Spanish province Castilla-La Mancha, to improve the understanding of the interactions between the soil, the vegetation and the atmosphere.Here results of energy balance studies at the Barrax site are given, one of the three intensively studied experimental sites within Castilla-La Mancha. This area is characterized by a large fraction of irrigated fields (40%) while the remaining 60% was fallow land at the end of June 1991. The energy balances over these two characteristic land-use classes totally differ. While for the irrigated fields the evapotranspiration is dominant, for the non-irrigated fields the sensible and the soil heat fluxes dominate and the latent heat flux is nearly negligible.In order to achieve areally averaged turbulent fluxes, surface, SODAR and aircraft data have been used. Comparing the surface fluxes from all three facilities, it can be found that:The sensible heat flux estimation from SODAR (_w-method) gives reliable results when a calibration of _w is done with another independent system (e.g. aircraft).Aircraft measurements in conjunction with energy budget methods yield surface fluxes of sensible heat, which are about 20% lower than the areally averaged values calculated by the surface measurements. The areally averaged latent heat fluxes from aircraft and surface measurements agree better than the sensible heat fluxes.     

Regional-Scale Heat and Water Vapour Fluxes in an Agricultural Landscape: An Evaluation of CBL Budget Methods at OASIS

This paper evaluates convective boundary layer (CBL) budget methods as a tool for estimating regionally averaged sensible and latent heat fluxes for the study region used in OASIS (Observations at Several Interacting Scales). This is an agricultural region of mixed cropping and grazing extending about 100 km west of the town of Wagga Wagga, NSW, Australia.The analysis proceeds in three stages: first, a simpleone-dimensional model of the well-mixed layer (the CBL slab model), forced with measurements of the surface heat and evaporation fluxes, is evaluated by comparing measured and modelled CBL temperature, humidity and depths. A comparison of several entrainment schemes shows that a simple model, where the entrainment kinetic energy is parameterised as a fraction (₃) of the surface sensible heat flux, works well if is set to 0.5. Second, the slab model is coupled to a Penman–Monteith model of surface evaporation to predict regional scale evaporation and thence heat fluxes. Finally, the integral CBL budget approach, which is an inverse method using theone-dimensional slab model, is used to infer regional heat and evaporation fluxes from measured time series of CBL temperature and humidity.We find that the simple CBL slab model works reasonably well for predicting CBL depth and very well for CBL temperature, especially if approximate estimates of subsidence velocity and warming due to advection are included. Regional sensible heat fluxes estimated from the integral CBL method match those measured, although the method is very sensitive to measurement errors. Measurement-model differences were larger for short integration times, because the well-mixed assumptions are violated at particular times of the day. The corollary is that `whole-day' (0530–1530 h) estimates are in reasonable agreement with measured values. Integral methods could not be used to infer the regional evaporation flux directly because CBL humidity profiles were complex and often not well mixed until mid-afternoon. We recommend that regional evaporation fluxes be predicted either from a coupled Penman–Monteith – CBL slab model, or inferred as a residual term from estimates of the regionally averaged available energy and sensible heat flux. Furthermore, we show that inferring fluxes via integral methods will always be difficult when the scalar concentrations have either a large surface source and free atmosphere sink (in the case of water vapour and methane), or a large surface sink and upper level source (in the case of CO₂).     

The Spatial Variability of Energy and Carbon Dioxide Fluxes at the Floor of a Deciduous Forest

Fluxes of carbon dioxide, water and sensible heat were measured using three different eddy covariance systems above the forest floor of a closed deciduous forest (leaf area index approx 6). The primary objective was to examine the representativeness of a single eddy covariance system in estimating soil respiration for time scales ranging from one-half hour to more than one week. Experiments were conducted in which the eddy covariance sensors were in one of three configurations: i) collocated, ii) separated horizontally or iii) separated vertically. A measure of the variation between the three systems (CV',related to the coefficient of variation) for half-hour carbon dioxide fluxes was 0.14 (collocated systems), 0.34 (vertically separated systems at 1, 2 and 4 m above the surface), and 0.57 (systems horizontally separated by 30 m). A similar variation was found for other scalar fluxes (sensible and latent heat). Variability between systems decreased as the number of half-hour sampling periods used to obtain mean fluxes was increased. After forty-eight hours (means from ninety-six half-hour samples), CV' values for carbon dioxide fluxes were 0.07, 0.09 and 0.16 in the collocated, vertically separated and horizontally separated experiments, respectively. The time dependence of variability has implications on the appropriateness of using short-term measurements in modelling validation studies. There are also implications concerning the appropriate number of half-hour samples necessary to obtain reliable causal relationships between flux data and environmental parameters. Based on the longer-term measurements, we also discuss the representativeness of a single eddy covariance system in long-term monitoring of soil respiration and evaporation beneath forest canopies using the eddy covariance method.     

Surface renewal analysis for sensible and latent heat flux density

High frequency temperature measurements were recorded at five heights and surface renewal (SR) analysis was used to estimate sensible heat flux density (H) over 0.1 m tall grass. Traces of the temperature data showed ramp-like structures, and the mean amplitude and duration of these ramps were used to calculate H using structure functions. Data were compared with H values measured with a sonic anemometer. Latent heat flux density (E) was calculated using an energy balance and the results were compared with E computed from the sonic anemometer data. SR analysis provided good estimates of H for data recorded at all heights but the canopy top and at the highest measurement level, which was above the fully adjusted boundary layer.     

Surface turbulent flux measurements over the Loess Plateau for a semi-arid climate change study

In order to provide high quality data for climate change studies, the data quality of turbulent flux measurements at the station of SACOL (Semi-Arid Climate & Environment Observatory of Lanzhou University), which is located on a semi-arid grassland over the Loess Plateau in China, has been analyzed in detail. The effects of different procedures of the flux corrections on CO₂, momentum, and latent and sensible heat fluxes were assessed. The result showed that coordinate rotation has a great influence on the momentum flux but little on scalar fluxes. For coordinate rotation using the planar fit method, different regression planes should be determined for different wind direction sectors due to the heterogeneous nature of the ground surface. Sonic temperature correction decreased the sensible heat flux by about 9%, while WPL correction (correction for density fluctuations) increased the latent heat flux by about 10%. WPL correction is also particularly important for CO₂ fluxes. Other procedures of flux corrections, such as the time delay correction and frequency response correction, do not significantly influence the turbulent fluxes. Furthermore, quality tests on stationarity and turbulence development conditions were discussed. Parameterizations of integral turbulent characteristics (ITC) were tested and a specific parameterization scheme was provided for SACOL. The ITC test on turbulence development conditions was suggested to be applied only for the vertical velocity. The combined results of the quality tests showed that about 62%–65% of the total data were of high quality for the latent heat flux and CO₂ flux, and as much as about 76% for the sensible heat flux. For the momentum flux, however, only about 35% of the data were of high quality.     

Design of land-air parameterization scheme (LAPS) for modelling boundary layer surface processes

Summary A land-air parametrization scheme (LAPS) describes mass, energy and momentum transfer between the land surface and the atmosphere. The scheme is designed as a software package which can be run as part of an atmospheric model or a stand-alone scheme. A single layer approach is chosen for the physical and biophysical scheme background. The scheme has six prognostic variables: two temperatures (one for the canopy vegetation and one for soil surface), one interception storage, and three soil moisture storage variables. The scheme's upper boundary conditions are: air temperature, water vapour pressure, wind speed, radiation and precipitation at some reference level within the atmospheric boundary layer. The sensible and latent heat are calculated using resistance representation. The evaporation from the bare soil is parametrized using the scheme. The soil part is designed as a three-layer model which is used to describe the vertical transfer of water in the soil.The performances of the LAPS scheme were tested using the results of meteorological measurements over a maize field at the experimental site De Sinderhoeve (The Netherlands). The predicted partitioning of the absorbed radiation into sensible and latent heat fluxes is in good agreement with observations. Also, the predicted leaf temperature agrees quite well with the observed values.With 9 Figures