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.     

Suburban-rural energy balance comparisons in summer for Vancouver,B.C.

Simultaneous energy balance observations at a rural and a suburban site in Vancouver, B.C. during the summer of 1983 are presented. The study is a follow-up to that conducted in 1980. Many of the 1980 results were unexpected and the present study seeks to assess their representativeness. The net radiant, turbulent sensible, and rural soil heat flux densities were measured directly. The suburban heat storage was parameterized and the turbulent latent heat flux densities were resolved as residuals in the energy balances. The 1983 average diurnal energy partitioning for both sites was typical of those quoted in the literature, suggesting that the 1980 results represent an extreme case. Suburban-rural differences showed the suburban area to have a 4% increase in net radiation, a 51% increase in turbulent sensible heat, and a 46% decrease in turbulent latent heat flux density. The values of the average daytime Bowen ratio were 0.46 and 1.28 for the rural and suburban areas, respectively. The sensible heat flux density exhibited relatively large values in the late afternoon and remained directed upward on many summer evenings. Large day-to-day variability in the relative magnitude of the suburban turbulent fluxes may have been due to synoptic influences. In this environment, the turbulent surface and mixed layers are closely coupled because of the low aerodynamic resistance over the rough surface.     

Evaluation of the profile and the resistance method for estimation of surface fluxes of momentum,sensible and latent heat

Summary Hourly lysimetric and micrometeorological data taken over a grass surface at the Meteorological Research Unit, Cardington U.K. have been analysed. A temperature difference and measurements of wind speed at only one height, combined with an independently estimated effective roughness length allowed sensible heat and momentum fluxes determination by the profile method on an hourly basis. The estimates are compared with direct measurements of sensible heat and friction velocity obtained by the eddy correlation method. The sensible and latent heat fluxes are also modelled by the resistance method. Equations based on the Monin—Obukhov similarity theory are used to account for stability effects through various forms of parameterization Aerodynamic and surface resistances, necessary for the Penman—Monteith equation are calculated from routinely measured meteorological data. The profile method for estimation of sensible heat flux and friction velocity is found to work excellently on the discussed daytime experimental data which correspond mainly to near neutral or slightly unstable conditions.Surface latent and sensible heat fluxes can also be described very well by the resistance method. A slightly better estimate of the sensible heat flux is achieved when stability corrections are taken into account. On the contrary Penman-Monteith equation for estimating latent heat flux is insensitive to adjustments for atmospheric stability.The comparison of the various methods leads to the establishment of empirical relationships which correlate various quantities such as soil heat flux, resistances, evapotranspiration etc. to routinely measured meteorological data.With 8 Figures     

Further Development of the Vegetated Urban Canopy Model Including a Grass-Covered Surface Parametrization and Photosynthesis Effects

The vegetated urban canopy model (VUCM), which includes parametrizations of urban physical processes for artificial surfaces and vegetated areas in an integrated system, has been further developed by including physical processes associated with grass-covered surfaces in urban pervious surfaces and the photosynthesis effects of urban vegetation. Using measurements made from three urban/suburban sites during the BUBBLE field campaign in 2002, the model’s performance in modelling surface fluxes (momentum flux, net radiation, sensible and latent heat fluxes and storage heat flux) and canopy air conditions (canopy air temperature and specific humidity) was critically evaluated for the non-precipitation and the precipitation days. The observed surface fluxes at the urban/suburban sites were significantly altered by precipitation as well as urban vegetation. Especially, the storage heat at urban surfaces and underlying substrates varied drastically depending on weather conditions while having an important role in the formation of a nocturnal urban surface layer. Unlike the nighttime canopy air temperature that was largely affected by the storage-heat release, the daytime canopy air conditions were highly influenced by the vertical turbulent exchange with the overlying atmosphere. The VUCM well reproduced these observed features in surface fluxes and canopy air conditions at all sites while performing well for both the non-precipitation and the precipitation days. The newly implemented parametrizations clearly improved the model’s performance in the simulation of sensible and latent heat fluxes at the sites, more noticeably at the suburban site where the vegetated area fraction is the largest among the sites. Sensitivity analyses for model input parameters in VUCM elucidated the relative importance of the morphological, aerodynamic, hydrological and radiative/thermal properties in modelling urban surface fluxes and canopy air conditions for daytime and nighttime periods. These results suggest that the VUCM has great potential for urban atmospheric numerical modelling for a range of cities and weather conditions in addition to having a better physical basis in the representation of urban vegetated areas and associated physical processes.     

The measurement of turbulent surface-layer fluxes by use of bichromatic scintillation

Scintillation measurements with a HeNe and a CO₂ laser were used to derive turbulent fluxes of heat and momentum in the surface layer. This was achieved by the structure constant or dissipation technique, i.e., by relating the measured structure constants and inner scales of refractive index fluctuations to structure constants of temperature fluctuations and dissipation rates of turbulent kinetic energy, respectively, and then assuming Monin-Obukhov similarity.The resulting heat fluxes agree well with measurements using the eddy correlation technique but for averaging periods of 10 min, the optical data show a much smoother and physically more plausible behaviour. The optically derived friction velocities are in good agreement with estimates derived from wind velocity and surface roughness. It was also observed that for stationary conditions, 1-min averaged optical measurements already provide good estimates for longer averaged heat and momentum fluxes.Even though some uncertainty remains about the empirical constants and Monin-Obukhov similarity expressions used, the method clearly proves to be of great value for monitoring surface-layer turbulence.     

Momentum- and Heat-Flux Parametrization at Dome C,Antarctica: A Sensitivity Study

An extensive meteorological observational dataset at Dome C, East Antarctic Plateau, enabled estimation of the sensitivity of surface momentum and sensible heat fluxes to aerodynamic roughness length and atmospheric stability in this region. Our study reveals that (1) because of the preferential orientation of snow micro-reliefs (sastrugi), the aerodynamic roughness length \(z_{0}\) varies by more than two orders of magnitude depending on the wind direction; consequently, estimating the turbulent fluxes with a realistic but constant \(z_{0}\) of 1 mm leads to a mean friction velocity bias of \(24\,\%\) in near-neutral conditions; (2) the dependence of the ratio of the roughness length for heat \(z_{0t}\) to \(z_{0}\) on the roughness Reynolds number is shown to be in reasonable agreement with previous models; (3) the wide range of atmospheric stability at Dome C makes the flux very sensitive to the choice of the stability functions; stability function models presumed to be suitable for stable conditions were evaluated and shown to generally underestimate the dimensionless vertical temperature gradient; as these models differ increasingly with increases in the stability parameter z / L, heat flux and friction velocity relative differences reached \(100\,\%\) when \(z/L > 1\); (4) the shallowness of the stable boundary layer is responsible for significant sensitivity to the height of the observed temperature and wind data used to estimate the fluxes. Consistent flux results were obtained with atmospheric measurements at heights up to 2 m. Our sensitivity study revealed the need to include a dynamical parametrization of roughness length over Antarctica in climate models and to develop new parametrizations of the surface fluxes in very stable conditions, accounting, for instance, for the divergence in both radiative and turbulent fluxes in the first few metres of the boundary layer.     

Surface fluxes estimation over agricultural areas. Comparisonof methods and the effects of land surface inhomogeneity

Summary  The Bowen Ratio-Energy Balance (BREB) and the aerodynamic method were used to estimate turbulent fluxes of sensible and latent heat flux over an irrigated agricultural area (IAA) and over two dry agricultural areas (DAA1 and DAA2). These turbulent fluxes were analysed and particular attention paid to two specific areas. First, a quantitative analysis of sensible and latent heat fluxes obtained by the BREB method was carried out, taking into account different soil type, vegetation and surface conditions. The results showed that in IAA latent heat flux was higher than sensible heat flux, except in summer months, while in DAA1 and DAA2, sensible heat flux was higher except in the months when the vegetation was at the stage of maximum development. Second, sensible and latent heat fluxes estimates from the BREB method were compared with those obtained from the aerodynamic method. In this comparison factors such as soil type, soil vegetation cover, homogeneity or inhomogeneity of terrain and mesoscale effects such as orography and wind patterns were taken into account. The results show that in conditions of light wind, the two methods only concur if the condition of horizontal homogeniety is fulfilled. The influence of inhomogeneity seems to decrease and agreement between methods improves, if the wind is stronger and the effects of meso and synoptic scales are predominant. Received May 18, 1999/Revised March 15, 2000     

Turbulence Measurements above a Pine Forest

Eddy fluxes of momentum, sensible and latent heat, and turbulence spectra measured over the Thetford Forest during 10 days in the Spring of 1973 are described. The measured total heat flux (H + λE) for 122 20-min periods agreed closely on average with independent estimates from an energy balance method. There was evidence that the energy balance data gave small systematic overestimates of available energy during the hours before noon, compensated by slight underestimates for the remainder of the day. A comparison of measured wind speeds and friction velocities in neutral stability confirmed the validity of the aerodynamic method for estimating momentum fluxes at heights of a few roughness lengths above the canopy. In stable conditions the log-linear wind profile $U=(u_{?t}/k)({? ln}((z-d)/z_{0})+←pha(z-d-z_{0})/L)$ with α = 3.4 ± 0.4 provided a good fit to the data. Spectra in unstable conditions were generally more sharply peaked than those measured by other workers over smoother terrain: differences were less marked in the case of vertical velocity in stable conditions. Temperature spectra in these stable conditions showed high energy at relatively low wavenumbers, and wT cospectra showed a cospectral gap; both of these results were associated with an intermittent sawtooth structure in the temperature fluctuations.     

Where does all the energy go? Surface energy partitioning in suburban Christchurch under stable wintertime conditions

Summary Typical observations of surface energy fluxes in urban areas usually employ the eddy covariance approach with measurements from a tower at a height well above the roughness elements. Net radiation and the turbulent fluxes are directly measured, the anthropogenic flux may be estimated and the storage flux is calculated as an energy balance residual. This paper reports both measurement and modelling of energy fluxes during wintertime in suburban areas of Christchurch, New Zealand. Under settled anticyclonic conditions, a strong inversion can occur over the city which severely restricts turbulent mixing. Even after sunrise the turbulent fluxes are small, and if one assumes advection is negligible, this means the storage flux increases in importance to very high levels. This paper suggests that these high storage rates are physically unrealistic for this environment. Rather, it is likely that some energy, which is assumed to be dissipated as heat storage, is more likely lost through mesoscale advection or attributed to errors caused by unsuitable measurement techniques under conditions with low friction velocity. However, even these two processes cannot fully account for flux loss. A full study is recommended to resolve this issue and meanwhile, caution is advised when applying current research methodology to estimate storage fluxes in urban areas in stable wintertime conditions.     

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.     

Pilot study using SPOT satellite imagery over Apalachicola national forest to determine appropriate spatial scale for area-wide aggregation of surface fluxes

Summary High resolution radiances from SPOT satellite imagery converted to Normalized Difference Vegetation Indices (NDVI) over a 16×16 km² mixed ground cover study-area in the Apalachicola National Forest in northwest Florida, along with in situ measurements from a Bowen ratio surface flux monitoring system and physical modeling techniques, are used to determine the length manifold beyond which degraded resolution satellite imagery fails to capture flux variability over the scene. The investigation is relevant to an understanding of how bias error is generated in methods designed to produce scale-invariant surface flux estimates from satellite measurements. Error estimates are based on assigning characteristic NDVI values to the four predominant types of ground cover found within the study-area. An open site near the center of the study-area, which satisfies the conditions for surface flux monitoring, is used for obtaining input data for a biosphere-atmosphere exchange model designed to calculate representative fluxes for the different ground covers. Continuous 6-minute meteorological and surface flux measurements were made at the monitoring site for a period of 22 days. These measurements are used in conjunction with surface layer theory to provide surface layer profile estimates of wind speed, temperature, and relative humidity at the tops of the forested sites. The measured and derived meteorological parameters, together with representative biophysical parameters, are used as input to the biosphere-atmosphere exchange model. By representing sensible and latent heat flux distributions due to the variable ground cover with characteristic NDVI values at 20-m resolution, baseline area-wide sensible and latent heat flux quantities are calculated. Error-growth curves as a function of spatial resolution for the fluxes are found by degrading the resolution of the SPOT radiances used to calculate NDVI, and rationing the associated area-wide fluxes to the baseline values. The point at which an error-growth curve becomes invariant represents the edge of a length manifold beyond which the satellite input no longer contains information on surface flux variability, even though NDVI variability continues at all scales up to that of the complete SPOT scene. The error-growth curves are non-linear, with all the error build-up taking place between 20 m and 1.6 km. Decreasing the spatial resolution of the NDVI information down to or below 1.6 km, introduces bias errors in the area-wide surface flux estimates of 10% for sensible heat and 8% for latent heat. The underlying assumptions and modeling produce uncertainty in estimating the manifold limits, however, the principal objective is to show that in using satellite data for scale-invariant surface flux retrieval, there is an optimal spatial resolution factor that can be objectively quantified.With 10 Figures     

青藏高原不同地区夏季近地层能量输送与微气象特征比较分析

利用2009年7月在青藏高原理塘、林芝、海北、拉萨获得的气象观测资料,对比分析了这些地区近地层气象要素、辐射收支及湍流通量日变化特征。结果表明:无论是高原东部、中部还是北部,无论是高原台地还是高山峡谷区,7月份近地层各气象要素、湍流通量、辐射收支都有明显的日变化。各地区的地表辐射、感热、潜热等最高值都出现在中午,最低值出现在早晨。地表反照率日变化均呈早晚高中午低的“U”型分布。地面热源强度在白天均为热源,正午为强热源,在夜间表现为弱的冷、热源交替出现。7月份近地层地气热量交换中,感热输送作用小,潜热输送占主导地位。动量通量和摩擦速度均在风速较大的下午较大,风速小的早晨小。      

Radiative surface temperature and energy balance of a wheat canopy

The sensible heat loss from a stand of winter wheat was calculated from radiometric measurements of crop surface temperature, measurements of air temperature, and an atmospheric resistance to momentum transfer; corresponding latent heat flux was obtained through the energy balance equation. These estimates of sensible and latent heat were compared with fluxes from the Bowen Ratio method. When radiative temperature was derived using a measured canopy emissivity of 0.98, calculations of sensible heat flux were systematically 50–100 W m^-2 less than Bowen Ratio values. The two techniques agreed more closely when an apparent emissivity of 0.96 was used with an apparent reflectivity of 0.03. The mean difference between the estimates of latent heat flux was then -16 ± 32 W m^-2.The surface temperature method showed less systematic error in comparison with the Bowen Ratio values than did estimates using the aerodynamic method.On leave from: University of Nottingham, School of Agriculture, Loughborough LE12 5RD.     

A scheme for computing surface fluxes from mean flow observations

A computational scheme is developed for estimating turbulent surface stress, sensible heat flux and humidity flux from mean velocity, temperature and humidity at a single height in the atmospheric surface layer; conditions at this reference level are presumed known from observations or from a numerical atmospheric circulation model. The method is based on coupling a Monin-Obukhov similarity profile to a force-restore formulation for the evolution of surface soil temperature to yield the local values of shear stress, heat flux and surface temperature. A self-contained formulation is presented including parameterizations for solar and infrared radiant flux at the surface.In addition to reference-level mean flow properties, the parameters needed to implement the scheme are thermal heat capacity of the soil, surface aerodynamic roughness, latitude, solar declination, surface albedo, surface emissivity and atmospheric transmissivity.Sample calculations are presented for (a), constant atmospheric forcing at the reference level, and (b) variable atmospheric forcing corresponding to Kahle's (1977) measurements of windspeed, air temperature and radiometer soil surface temperature under dry vegetatively sparse conditions in the Mohave Desert in California. The latter case simulated the observed diurnal variations resonably well for the parameters used.Consultant, Atmospheric Sciences Division, Department of Energy and Environment, Brookhaven National Laboratory, Upton, N.Y., pc11973, U.S.A.     

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.     

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.     

A model for sensible heat flux probability density function for near-neutral and slightly-stable atmospheric flows

The probability density function for sensible heat flux was measured above a uniform dry lakebed (Owens lake) in Owens Valley, California. It was found that for moderately stable to near neutral atmospheric stability conditions, the probability density function exhibits well defined exponential tails. These exponential tails are consistent with many laboratory boundarylayer measurements and numerical simulations. A model for the sensible heat flux probability density function was developed and tested. A key assumption in the model derivation was the near Gaussian statistics of the vertical velocity and temperature fluctuations. This assumption was verified from time series measurements of temperature and vertical velocity. The parameters for the sensible heat flux probability density function model were also derived from mean meteorological and surface conditions using surface-layer similarity theory. It was found that the best agreement between modeled and measured sensible heat flux probability density function was at the tails. Finally, a relation between the intermittency parameter, the probability density function, and the mean meteorological conditions was derived. This relation rigorously links the intermittency parameter to mean meteorological conditions.     

成都精细下垫面信息对城市气象影响的模拟试验

为了提高成都市精细化天气预报水平,使用成都地区精细下垫面土地利用资料,在WRF中耦合了单层城市冠层模式,对2008年7月6 日晴空背景下的成都城市气象特征进行了模拟,并和使用旧土地利用资料、slab模式的模拟结果进行了对比分析.模拟结果表明城区因为不透水下垫面的增加,使得地表蒸发和地表水汽通量显著减小,潜热通量减小,感...     

Aircraft Observations of Sea-Surface Turbulent Fluxes Near the California Coast

Aircraft turbulence data from the Autonomous Ocean Sampling Network project were analyzed and compared to the Coupled Ocean–Atmosphere Response Experiment (COARE) bulk parametrization of turbulent fluxes in an ocean area near the coast of California characterized by complex atmospheric flow. Turbulent fluxes measured at about 35 m above the sea surface using the eddy-correlation method were lower than bulk estimates under unstable and stable atmospheric stratification for all but light winds. Neutral turbulent transfer coefficients were used in this comparison because they remove the effects of mean atmospheric conditions and atmospheric stability. Spectral analysis suggested that kilometre-scale longitudinal rolls affect significantly turbulence measurements even near the sea surface, depending on sampling direction. Cross-wind sampling tended to capture all the available turbulent energy. Vertical soundings showed low boundary-layer depths and high flux divergence near the sea surface in the case of sensible heat flux but minimal flux divergence for the momentum flux. Cross-wind sampling and flux divergence were found to explain most of the observed discrepancies between the measured and bulk flux estimates. At low wind speeds the drag coefficient determined with eddy correlation and an inertial dissipation method after corrections were applied still showed high values compared to bulk estimates. This discrepancy correlated with the dominance of sea swell, which was a usually observed condition under low wind speeds. Under stable atmospheric conditions measured sensible heat fluxes, which usually have low values over the ocean, were possibly affected by measurement errors and deviated significantly from bulk estimates.     

北京城市下垫面边界层湍流统计特征

分析了1998年10月在北京325 m气象塔47 m、120m及280m 3层取得的湍流资料,讨论了湍流宏观统计特征量等随稳定度变化的规律.结果表明,47 m及120 m观测到的感热通量的日变化基本一致,城市边界层中近地层高度可达100m左右,这主要是由于下垫面粗糙度很大所致.同时发现当冷空气来临时,边界层中的感热通量会有明显增大.