A Modelling Study of Flux Imbalance and the Influence of Entrainment in the Convective Boundary Layer

It is frequently observed in field experiments that the eddy covariance heat fluxes are systematically underestimated as compared to the available energy. The flux imbalance problem is investigated using the NCAR’s large-eddy simulation (LES) model imbedded with an online scheme to calculate Reynolds-averaged fluxes. A top–down and a bottom–up tracer are implemented into the LES model to quantify the influence of entrainment and bottom–up diffusion processes on flux imbalance. The results show that the flux imbalance follows a set of universal functions that capture the exponential decreasing dependence on u _*/w _*, where u _* and w _* are friction velocity and the convective velocity scale, respectively, and an elliptic relationship to z/z _i , where z _i is the mixing-layer height. The source location in the boundary layer is an important factor controlling the imbalance magnitude and its horizontal and vertical distributions. The flux imbalance of heat and the bottom–up tracer is tightly related to turbulent coherent structures, whereas for the top–down diffusion, such relations are weak to nonexistent. Our results are broadly consistent with previous studies on the flux imbalance problem, suggesting that the published results are robust and are not artefacts of numerical schemes.     

Consequences of Incomplete Surface Energy Balance Closure for CO2 Fluxes from Open-Path CO2/H2O Infrared Gas Analysers

We present an approach for assessing the impact of systematic biases in measured energy fluxes on CO₂ flux estimates obtained from open-path eddy-covariance systems. In our analysis, we present equations to analyse the propagation of errors through the Webb, Pearman, and Leuning (WPL) algorithm [Quart. J. Roy. Meteorol. Soc. 106, 85–100, 1980] that is widely used to account for density fluctuations on CO₂ flux measurements. Our results suggest that incomplete energy balance closure does not necessarily lead to an underestimation of CO₂ fluxes despite the existence of surface energy imbalance; either an overestimation or underestimation of CO₂ fluxes is possible depending on local atmospheric conditions and measurement errors in the sensible heat, latent heat, and CO₂ fluxes. We use open-path eddy-covariance fluxes measured over a black spruce forest in interior Alaska to explore several energy imbalance scenarios and their consequences for CO₂ fluxes.     

Consequences of Uncertainties in CO<Subscript>2</Subscript> Density for Estimating Net Ecosystem CO<Subscript>2</Subscript> Exchange by Open-path Eddy Covariance

Errors in the estimation of CO₂ surface exchange by open-path eddy covariance, introduced during the removal of density terms [Webb et al. Quart J Roy Meteorol Soc 106:85–100, (1980) - WPL], can happen both because of errors in energy fluxes [Liu et al. Boundary-Layer Meteorol 120:65–85, (2006)] but also because of inaccuracies in other terms included in the density corrections, most notably due to measurements of absolute CO₂ density (ρ _c ). Equations are derived to examine the propagation of all errors through the WPL algorithm. For an open-path eddy covariance system operating in the Sierra de Gádor in south-east Spain, examples are presented of the inability of an unattended, open-path infrared gas analyzer (IRGA) to reliably report ρ _c and the need for additional instrumentation to determine calibration corrections. A sensitivity analysis shows that relatively large and systematic errors in net ecosystem exchange (NEE) can result from uncertainties in ρ _c in a semi-arid climate with large sensible heat fluxes (H _s ) and (wet) mineral deposition. When ρ_c is underestimated by 5% due to lens contamination, this implies a 13% overestimation of monthly CO₂ uptake.     

Surface energy balance measurements over a banana plantation in South China

The land surface energy exchange depends highly on the surface properties. Little is known of the energy balance over a typical banana plantation of humid tropics. In this study, we examine the characteristics of surface energy exchange over a typical banana field in South China during the period of May 2010 to April 2011 by using the eddy covariance and micrometeorological tower. The results showed that the diurnal and seasonal variations in surface latent heat flux were larger compared with those over the nearby grassland. The dominant energy partitioning varies with season. The latent heat flux was the main consumer of net radiation in summer, whereas the sensible heat flux was the main consumer in winter. The increasing cloud coverage and rain appear to control the surface energy balance with the development of the monsoon. Due to increased afternoon convective cloud systems in the monsoon active period, downward shortwave radiation was dramatically diminished around 14:00?pm. The annual mean Bowen ratio was 0.69, which fell within the range of other vegetated surfaces. The observed surface energy components were not closed, and the ratio of turbulent fluxes to the available energy was about 77 % in October–January and about 85 % in the other months after considering soil heat and air heat storage.     

Processing and quality control of flux data during LITFASS-2003

Different aspects of the quality assurance and quality control (QA/QC) of micrometeorological measurements were combined to create a comprehensive algorithm which was then applied to experimental data from LITFASS-2003 (Lindenberg Inhomogeneous Terrain—Fluxes between Atmosphere and Surface: a long term Study). Eddy-covariance measurements of the latent heat flux were the main focus of the QA/QC efforts. The results of a turbulence sensor intercomparison experiment showed deviations between the different eddy-covariance systems on the order of 15%, or less than 30 W m⁻², for the latent heat flux and 5%, or less than 10 W m⁻², for the sensible heat flux. In order to avoid uncertainties due to the post-processing of turbulence data, a comprehensive software package was used for the analysis of experimental data from LITFASS-2003, including all necessary procedures for corrections and quality control. An overview of the quality test results shows that for most of the days more than 80% of the available latent heat flux data are of high quality so long as there are no instrumental problems. The representativeness of a flux value for the target land-use type was analysed using a stochastic footprint model. Different methods to calculate soil heat fluxes at the surface are discussed and a sensitivity analysis is conducted to select the most robust method for LITFASS-2003. The lack of energy balance closure, which was found for LITFASS-2003, can probably be attributed to the presence of low-frequency flux contributions that cannot be resolved with an averaging time of 30 min. Though the QA/QC system has been developed for the requirements of LITFASS-2003, it can also be applied to other experiments dealing with similar objectives.     

Towards Closing the Surface Energy Budget of a Mid-latitude Grassland

Observations for May and August, 2005, from a long-term grassland meteorological station situated in central Netherlands were used to evaluate the closure of the surface energy budget. We compute all possible enthalpy changes, such as the grass cover heat storage, dew water heat storage, air mass heat storage and the photosynthesis energy flux, over an averaging time interval. In addition, the soil heat flux was estimated using a harmonic analysis technique to obtain a more accurate assessment of the surface soil heat flux. By doing so, a closure of 96% was obtained. The harmonic analysis technique appears to improve closure by 9%, the photosynthesis for 3% and the rest of the storage terms for a 3% improvement of the energy budget closure. For calm nights (friction velocity u _* < 0.1 m s⁻¹) when the eddy covariance technique is unreliable for measurement of the vertical turbulent fluxes, the inclusion of a scheme that calculates dew fluxes improves the energy budget closure significantly.     

On Monin–Obukhov Scaling in and Above the Atmospheric Surface Layer: The Complexities of Elevated Scintillometer Measurements

In scintillometry Monin–Obukhov similarity theory (MOST) is used to calculate the surface sensible heat flux from the structure parameter of temperature (C_T²){(C_{T^2})} . In order to prevent saturation a scintillometer can be installed at an elevated level. However, in that case the observation level might be located outside the atmospheric surface layer (ASL) and thus the validity of MOST questioned. Therefore, we examine two concepts to determine the turbulent surface sensible heat flux from the structure parameter at elevated levels with data obtained at 60-m height on the Cabauw tower (the Netherlands). In the first concept (MOSTs) C_T²{C_{T^2}} is still scaled with the surface flux, whereas in the second (MOSTl) C_T²{C_{T^2}} is scaled with the local sensible heat flux. The C_T²{C_{T^2}} obtained from both concepts is compared with direct observations of C_T²{C_{T^2}} using a sonic anemometer/thermometer. In the afternoon (when the measurement height is located within the ASL) both concepts give results that are comparable to the directly observed values of C_T²{C_{T^2}} . In the morning (data outside the ASL), our data do not unequivocally support either of the two concepts. First, the peak in C_T²{C_{T^2}} that occurs when the measurement height is located in the entrainment zone disqualifies the use of MOST. Second, during the morning transition, local scaling shows the correct pattern (zero flux and a minimum in C_T²{C_{T^2}}) but underestimates C_T²{C_{T^2}} by a factor of ten. Third, from the best linear fit a we found that the slope of MOSTl gave better results, whereas the offset is closer to zero for MOSTs. Further, the correlation between the direct observations and MOST-scaled results is low and similar for the two concepts. In the end, we conclude that MOST is not applicable for the morning hours when the observation level is above the ASL.     

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.     

Analysis of surface energy balance closure over heterogeneous surfaces

Surface energy balance closure has been examined using eddy covariance measurements and other observations at one industrial and three agricultural sites near the Nakdong River during daytime. Energy balance closure was evaluated by calculating the long-term averaged energy balance ratio (EBR), the ratio of turbulent energy fluxes to available energy, and the statistical regression of turbulent energy fluxes against available energy using half-hourly data. The EBR of all sites ranges from 0.46 to 0.83 while the coefficient of determination (R ²) ranges from 0.37 to 0.77. The energy balance closure was relatively poor compared to homogeneous sites, indicating the influence of surface heterogeneity. Unmeasured heat storage terms also seem to play a role in the surface energy budget at the industrial and irrigated sites. The energy balance closure was better in conditions of high wind speed, low downward short wave radiation, and high friction velocity, which suggests the role of heat storage term and surface heterogeneity in surface energy balance at these sites. Spectrum analysis shows a sharp roll-off at the low frequency in co-spectrum, which indicates that low-frequency motions do not significantly contribute to turbulent fluxes. Both the spectra and cospectra in unstable conditions show a broad peak indicating the influence of multiple sizes of large eddies over heterogeneous sites. Most of ogive curves for the kinematic latent and sensible heat fluxes reach an asymptote within 30 minutes regardless of the EBR value, indicating that low frequency motion is not a main factor for energy imbalance. However, stationary eddies due to landscape heterogeneity still remains as a possible cause for energy imbalance.     

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.     

Buoyancy and The Sensible Heat Flux Budget Within Dense Canopies

In contrast to atmospheric surface-layer (ASL) turbulence, a linear relationship between turbulent heat fluxes (F_T) and vertical gradients of mean air temperature within canopies is frustrated by numerous factors, including local variation in heat sources and sinks and large-scale eddy motion whose signature is often linked with the ejection-sweep cycle. Furthermore, how atmospheric stability modifies such a relationship remains poorly understood, especially in stable canopy flows. To date, no explicit model exists for relating F_T to the mean air temperature gradient, buoyancy, and the statistical properties of the ejection-sweep cycle within the canopy volume. Using third-order cumulant expansion methods (CEM) and the heat flux budget equation, a “diagnostic” analytical relationship that links ejections and sweeps and the sensible heat flux for a wide range of atmospheric stability classes is derived. Closure model assumptions that relate scalar dissipation rates with sensible heat flux, and the validity of CEM in linking ejections and sweeps with the triple scalar-velocity correlations, were tested for a mixed hardwood forest in Lavarone, Italy. We showed that when the heat sources (S_T) and F_T have the same sign (i.e. the canopy is heating and sensible heat flux is positive), sweeps dominate the sensible heat flux. Conversely, if S_T and F_T are opposite in sign, standard gradient-diffusion closure model predict that ejections must dominate the sensible heat flux.     

Effects of the soil heat flux estimates on surface energy balance closure over a semi-arid grassland

Soil heat flux is important for surface energy balance (SEB), and inaccurate estimation of soil heat flux often leads to surface energy imbalance. In this paper, by using observations of surface radiation fluxes and soil temperature gradients at a semi-arid grassland in Xilingguole, Inner Mongolia, China from June to September 2008, the characters of the SEB for the semi-arid grassland were analyzed. Firstly, monthly averaged diurnal variations of SEB components were revealed. A 30-min forward phase displacement of soil heat flux (G) observed by a fluxplate at the depth of 5-cm below the soil surface was conducted and its effect on the SEB was studied. Secondly, the surface soil heat flux (G _s) was computed by using harmonic analysis and the effect of the soil heat storage between the surface and the fluxplate on the SEB was examined. The results show that with the 30-min forward phase displacement of observed G, the slope of the ordinary linear regression (OLR) of turbulent fluxes (H+LE) against available energy (R _n-G) increased from 0.835 to 0.842, i.e., the closure ratio of SEB increased by 0.7%, yet energy imclosure of 15.8% still existed in the SEB. When G _s, instead of G was used in the SEB equation, the slope of corresponding OLR of (H+LE) against (R _n-G _s) reached 0.979, thereby the imclosure ratio of SEB was reduced to only 2.1%.     

Energy Fluxes of an Open Water Area in a mid-Latitude Prairie Wetland

Concurrent measurements of the surface energy balance components (net radiation, heat storage, and sensible and latent heat fluxes) were made in three communities (open water, Phragmites australis, Scirpus acutus) in a wetland in north-central Nebraska, U.S.A., during May-October, 1994. The Bowen ratio – energy balance method was used to calculate latent and sensible heat fluxes. This paper presents results from the open water area. The heat stored in water (G) was found to play a major role in the energy exchange over the water surface. During daytime, G consumed 45–60% of R _n , the net radiation (seasonally averaged daytime G was about 127 W m^–2). At night, G was a significant source of energy (seasonally averaged nighttime G was about -135 Wm^–). The diurnal pattern of latent heat flux ( E) did not follow that of R _n . On some days, E was near zero during midday periods with large R _n . The diurnal variability in E seemed to be significantly affected by temperature inversions formed over the cool water surface. The daily evaporation rate (E) ranged from 2 to 8 mm during the measurement period, and was generally between 70 and 135% of the equilibrium rate.     

EVALUATION OF CORRECTIONS ON TURBULENT FLUXES OBTAINED BY EDDY COVARIANCE METHOD IN HIGH WINDS

This paper investigates processing of fast-response data and corrections of turbulent fluxes obtained by using eddy covariance method based on data collected at an offshore observation tower during three Cold-intrusion(CI)events in the South China Sea in 2010. This study presents the data processing procedure in detail and compares frictional velocities(u*), sensible heat fluxes(H) and latent heat fluxes(LE) yielded by using different averaging periods and different coordinate rotation methods; evaluates the sonic temperature correction for sensible heat flux and the Webb correction for latent heat flux as a function of 10 m wind speed(u10) during the CIs. The results show(1) that the different averaging periods of 30 min and 10 min cause biases of u*(H, LE) within 5%(15%, 62%). The values of u*(H,LE) averaged from 30 mins are mostly larger than those averaged from 10 mins. We suggest that the averaging period of 10 min is not sufficiently long to capture all scale eddies and recommend 30 min averaging period in calculating turbulent fluxes using eddy covariance method during CIs;(2) that the values of u*(H, LE) obtained from double rotation(DR2) and those obtained from planar fit rotation(PF) have good agreements and correlation coefficients between them are larger than 0.99. Because PF method requires unchanged environment and it is easier to apply DR2 method, we suggest DR2 coordinate rotation method in processing fast-response data; and(3) that the median values of frictional velocity(sensible heat flux and latent heat flux) binned according to 2 m s~(-1) intervals of u_(10) increase(decrease,increase) by less than 9%(4%, 10%) by Coriolis corrections(sonic temperature corrections, Webb corrections), which decreases(decreases, increases) with increasing u10 when u10 are 5-17 m s~(-1).     

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.     

Vertical divergence of fogwater fluxes above a spruce forest

Two almost identical eddy covariance measurement setups were used to measure the fogwater fluxes to a forest ecosystem in the “Fichtelgebirge” mountains (Waldstein research site, 786 m a.s.l.) in Germany. During the first experiment, an intercomparison was carried out with both setups running simultaneously at the same measuring height on a meteorological tower, 12.5 m above the forest canopy. The results confirmed a close agreement of the turbulent fluxes between the two setups, and allowed to intercalibrate liquid water content (LWC) and gravitational fluxes. During the second experiment, the setups were mounted at a height of 12.5 and 3 m above the canopy, respectively. For the 22 fog events, a persistent negative flux divergence was observed with a greater downward flux at the upper level. To extrapolate the turbulent liquid water fluxes measured at height z to the canopy of height h_c, a conversion factor 1/[1+0.116(z−h_c)] was determined. For the fluxes of nonvolatile ions, no such correction is necessary since the net evaporation of the fog droplets appears to be the primary cause of the vertical flux divergence. Although the net evaporation reduces the liquid water flux reaching the canopy, it is not expected to change the absolute amount of ions dissolved in fogwater.     

The Simulation of the Opposing Fluxes of Latent Heat and CO<Subscript>2</Subscript> over Various Land-Use Types: Coupling a Gas Exchange Model to a Mesoscale Atmospheric Model

A mesoscale meteorological model (FOOT3DK) is coupled with a gas exchange model to simulate surface fluxes of CO₂ and H₂O under field conditions. The gas exchange model consists of a C3 single leaf photosynthesis sub-model and an extended big leaf (sun/shade) sub-model that divides the canopy into sunlit and shaded fractions. Simulated CO₂ fluxes of the stand-alone version of the gas exchange model correspond well to eddy-covariance measurements at a test site in a rural area in the west of Germany. The coupled FOOT3DK/gas exchange model is validated for the diurnal cycle at singular grid points, and delivers realistic fluxes with respect to their order of magnitude and to the general daily course. Compared to the Jarvis-based big leaf scheme, simulations of latent heat fluxes with a photosynthesis-based scheme for stomatal conductance are more realistic. As expected, flux averages are strongly influenced by the underlying land cover. While the simulated net ecosystem exchange is highly correlated with leaf area index, this correlation is much weaker for the latent heat flux. Photosynthetic CO₂ uptake is associated with transpirational water loss via the stomata, and the resulting opposing surface fluxes of CO₂ and H₂O are reproduced with the model approach. Over vegetated surfaces it is shown that the coupling of a photosynthesis-based gas exchange model with the land-surface scheme of a mesoscale model results in more realistic simulated latent heat fluxes.     

Forest-Air Fluxes Of Carbon, Water And Energy Over Non-Flat Terrain

A field study of surface-air exchange of carbon, water, and energy was conducted at a mid-latitude, mixed forest on non-flat terrain to investigate how to best interpret biological signals from the eddy flux data that may be subject to advective influences. It is shown that during periods of Southwest winds (sector with mild topography), the eddy fluxes are well-behaved in terms of energy balance closure, the existence of a constant flux layer, consistency with chamber observations and the expected abiotic controls on the fluxes. Advective influences are evident during periods with wind from a steep (15%) slope to the Northeast of the tower. These influences appear more severe on CO₂ flux, particularly in stable air, than on the energy fluxes. Large positive flux of CO₂ (> 23 mol m^-2 s^-1) occurs frequently at night. The annual sum of the carbon flux is positive, but the issue about whether the forest is a source of atmospheric carbon remains inconclusive.Attempts are made to assess vertical advectionusing the data collected on a single tower. Over the Southwestsector, vertical advection makes a statistically significant but small contribution to the 30-min energy imbalance and CO₂ flux variations. Contributions by horizontal advection may be larger but cannot be verified directly by the current experimental method.     

简单生物圈模式 (SiB2) 中湍流能量通量对近地层两类阻抗系数的敏感性研究

陆—气之间的能量交换是通过近地层湍流热量和水汽通量来实现的.以往的研究表明近地层阻抗对湍流能量通量有着不同程度的影响,但是关于陆面模式中阻抗系数的取值范围却始终没有统一的标准.为了深入了解简单生物圈模式 (SiB₂) 中近地层阻抗系数取值变化对湍流能量通量的影响,我们以那曲站为例,分别采用传统的逐个因子分析法和考虑参数间相互作用的部分因子分析法定量地研究了夏季该观测站近地层湍流能量通量分别对冠层阻抗系数C₁和地表阻抗系数C₂的敏感性响应.结果表明,感热通量对地表阻抗系数C₂更为敏感,而潜热通量则对冠层阻抗系数C₁较为敏感;感热通量随C₁增加而增大,随C₂增加而减小,而潜热通量则随C₁或C₂的增加而减小;不管是感热通量还是潜热通量,它们对阻抗系数的敏感度随阻抗系数的增大而减小,而对阻抗系数的相对敏感度则随阻抗系数的增大而增大.最后,结合那曲站夏季下垫面稀疏短草的分布特点分析了造成感热通量和潜热通量敏感变化各异的原因.     

高温和台风影响下的苏州城市粗糙子层湍流通量特征

利用苏州地区2011 年12月20 日-2012年8 月13 日的湍流观测资料对不同季节、高温、台风强天气过程下的湍流特征进行分析.结果表明:城市地区不同季节动量通量、感热通量、潜热通量日变化明显,各通量的夏季平均值、最大值均高于冬春季,夏季感热通量日最大值为160.2 W·m-2,感热在城市地表能量平衡中的作用大于潜...