Methods to Estimate Surface Fluxes of Momentum and Heat from Routine Weather Observations for Dispersion Applications under Stable Stratification

We examine the efficacy of two methods commonly used to estimate the vertical turbulent fluxes of momentum and sensible heat from routinely observed mean quantities in the surface layer under stable stratification. The single-level method uses mean wind speed and temperature measurements at a single height, whereas the two-level method uses mean wind speed measurements at a single height and mean temperature measurements at two heights. These methods are used in popular meteorological processors such as the U.S. Environmental Protection Agency approved AERMET and CALMET to generate inputs for dispersion simulations. We use data from a flux station of the U.K. Met Office at Cardington for comparison. It is found that the single-level method does not describe the flux variation in the weakly stable regime at all, because of its assumption that the temperature scale, i.e. the ratio of the kinematic sensible heat flux to the friction velocity, is constant, which is plausible only under strongly stable conditions. On the other hand, the two-level method provides a physically realistic variation of the fluxes with stability, but the required temperature measurements at two levels are usually not available on a routine basis. If measurements of the standard deviation of temperature are also available, in addition to the mean temperature at a single level, then they can be usefully employed in a third (single-level) method, with the consequence that the computed fluxes are very similar to those obtained from the two-level method. An improvement to the original single-level method is considered, and flux calculations under low wind conditions are also discussed.     

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     

Using Temperature Fluctuation Measurements to Estimate Meteorological Inputs for Modelling Dispersion During Convective Conditions in Urban Areas

We examine the performance of several methods to estimate meteorological inputs for modelling dispersion in urban areas during convective conditions. Sensible heat flux, surface friction velocity and turbulent velocities are estimated from measurements of mean wind speed and the standard deviation of temperature fluctuations at a single level on a tower at two suburban sites and at one urban site in Riverside, California. These estimates are compared with observations made at these sites during a field study conducted in 2007. The sensible heat flux is overestimated in the urban area, while it is underestimated at a suburban site when temperature fluctuations are used in the free convection formulation to estimate heat flux. The bias in heat flux estimates can be reduced through a correction that depends on stability. It turns out that the bias in heat flux estimates has a minor effect on the prediction of surface friction velocity and turbulent velocities. Estimates of sensible heat flux, surface friction velocity and turbulent velocities are sensitive to estimates of aerodynamic roughness length, and we suggest estimating the aerodynamic roughness length through detailed micrometeorological measurements made during a limited field study. An examination of the impact of the uncertainty in estimating surface micrometeorology on concentrations indicates that, at small distances from a surface release, ground-level concentrations computed using estimates of heat flux and surface friction compare well with the those based on observed values: the bias is small and the 95% confidence interval of the ratio of the two concentrations is 1.7. However, at distances much larger than the Obukhov length, this confidence interval is close to 2.3 because errors in both friction velocity and heat flux affect plume spread. Finally, we show that using measurements of temperature fluctuations in estimating heat flux is an improvement on that based on the surface energy balance, even when net radiation measurements are available.     

The exponential probability density function and concentration fluctuations in smoke plumes

Observations of 1-s average concentration fluctuations during two trials of a U.S. Army diffusion experiment are presented and compared with model predictions based on an exponential probability density function (pdf). The source is near the surface and concentration monitors are on lines about 30 to 100 m downwind of the source. The observed ratio of the standard deviation to the mean of the concentration fluctuations is about 1.3 on the mean plume axis and 4 to 5 on the mean plume edges. Plume intermittency (fraction of non-zero readings) is about 50%; on the mean plume axis and 10%; on the mean plume edges. A meandering plume model is combined with an exponential pdf assumption to produce predictions of the intermittency and the standard deviation of the concentration fluctuations that are within 20%; of the observations.     

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.     

Validation of the Canadian land surface scheme (class) for a subarctic open woodland

Abstract

This study reports on testing of the peatland version of the Canadian Land Surface Scheme (CLASS) for simulating the energy balance of subarctic open woodland terrain. Model results are compared against several years of measured data from a site near Churchill, Manitoba. In contrast to most forest environments, the floor of the open forest plays a large role in total ecosystem energy exchange. This behaviour presents a significant challenge for land surface models like CLASS and their simplified treatment of vegetation canopies.

Simulations of summer energy balance for seven years encompassing a wide range of meteorological conditions produced consistent results. Root mean square errors for sensible and latent heat fluxes fell between 11 and 28 W m^?2. CLASS consistently underestimated slightly the daily latent heat flux and overestimated the sensible heat flux, average mean bias errors being ‐7.6 and 9.1 W m^?2, respectively. The soil heat flux was less well represented. In general, CLASS was able to capture the diurnal and seasonal behaviour of the measured fluxes under a range of conditions with reasonable accuracy.

In a full year simulation, CLASS reproduced the annual variations in energy balance with some discrepancies associated with snow accumulation and ablation periods. The model performance was sensitive to both snow density and specification of the surface cover. Recommendations for improving the model for subarctic woodlands and terrain types with similar features are discussed.     

A Re-examination of Density Effects in Eddy Covariance Measurements of CO2 Fluxes

Corrections of density effects resulting from air-parcel expansion/compression are important in interpreting eddy covariance fluxes of water vapor and CO2 when open-path systems are used. To account for these effects, mean vertical velocity and perturbation of the density of dry air are two critical parameters in treating those physical processes responsible for density variations. Based on various underlying assumptions, different studies have obtained different formulas for the mean vertical velocity and perturbation of the density of dry air, leading to a number of approaches to correct density effects. In this study, we re-examine physical processes related to different assumptions that are made to formulate the density effects. Specifically, we re-examine the assumptions of a zero dry air flux and a zero moist air flux in the surface layer, used for treating density variations, and their implications for correcting density effects. It is found that physical processes in relation to the assumption of a zero dry air flux account for the influence of dry air expansion/compression on density variations. Meanwhile, physical processes in relation to the assumption of a zero moist air flux account for the influence of moist air expansion/compression on density variations. In this study, we also re-examine mixing ratio issues. Our results indicate that the assumption of a zero dry air flux favors the use of the mixing ratio relative to dry air, while the assumption of a zero moist air flux favors the use of the mixing ratio relative to the total moist air. Additionally, we compare different formula for the mean vertical velocity, generated by air-parcel expansion/compression, and for density effect corrections using eddy covariance data measured over three boreal ecosystems.     

青藏高原五道梁地区湍流输送特征的研究

根据１９９４年６－７月在青藏高原五道梁地区的湍流脉动观测资料，分析了该地区近地层能量平衡、感热和谱热的日变化及湍流强度和端流谱特征。结果表明：晴天该地区近地层能量基本平衡，各能量分量的日变化与常情况相同；白天感热通量的输送占主导地位，潜热通量占次要地位，符合半干旱的一般特征。     

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.     

不同气象条件下湍流通量的研究

利用1989年12月－1990年1月在重庆市测得的三轴风速仪资料和气象资料，1990年7月－8月在天津测得的三轴风速仪资料、温度脉冲资料、气象资料、计算了不同气象条件下的动量通量、热通量、摩擦速度和特征温度。这对了解重庆和天津市郊不同气象条件下的近地面层的湍流通量特征是十分有益的。     

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.     

On the Prediction of the Kurtosis of Velocity Derivatives in a Turbulent Field

The prediction of the values of non-dimensional fourth-order moment (kurtosis) of the velocity derivative in a turbulent field is made under the assumption that the values of kurtosis depend on both the turbulence Reynolds number and the intermittency factor. The method consists of modeling a suitable probability density of the variable in a given turbulence Reynolds number and the intermittency factor. A crude model of the probability density function is derived, and the numerical calculations based on the model show excellent agreement with many of the experimental data. The analysis shows that the values of kurtosis depend strongly on the intermittency factor, and that depending on the value of the intermittency factor, it is entirely possible to have values of kurtosis as low as five in a flow with a turbulence Reynolds number of 5000.     

A simple model for describing the heat balance of a shallow lake with application to Lake Chad

A model is described allowing the computation of a shallow lake's surface temperature as a function of the meteorological conditions of the surroundings. An equilibrium temperature with no heat flux through the surface is defined. The daily cycle of the equilibrium temperature permits computation of the average and the amplitude of the lake's temperature in the form of a simple sine wave. The computed values are compared with satellite measurements of the surface temperature and good agreement is achieved. When the computed surface temperature is applied, the values of the fluxes of latent and sensible heat are in good agreement with observational data.     

An evaluation of aircraft flux measurements of CO2, water vapor and sensible heat

Ground-based flux measurements of carbon dioxide and water vapor integrate physiological processes taking place on a field scale. Aircraft flux measurements have recently been undertaken to attempt to widen the scope of applicability of such measurements. However, because of the intermittency of turbulent transfer, flux measurements must be averaged over long periods of time or long distances to give reproducible results. This requirement makes it difficult to relate aircraft flux measurements to local surface processes. Flux measurements of CO₂, latent and sensible heat obtained from repeated passes in four directions and at three elevations over a homogeneous wheat-growing area are compared with ground-based measurements. Averages based on four runs of 4 km in length gave results consistent with ground-based measurements. The largest percentage differences were in the sensible heat flux. Cospectral analyses showed no significant high frequency losses for the data from flight levels of 25 and 50 m, but an underestimation of approximately 10% resulted at 10 m. Flight direction with respect to wind direction was relatively unimportant at 10 and 25 m but some effects were observed at 50 m. It was also shown that at 25 m, over a relatively smooth and homogeneous surface, the means of either three or four runs 4 km in length were similar to the means of 12–16 km runs. This confirms that at this altitude, most of the flux contribution is contained at wavelengths less than 4 km and that the mean of 3 to 4 passes accounts for most of the intermittency of turbulent transfer.     

On the prediction of the kurtosis of velocity derivatives in a turbulent field

The prediction of the values of non-dimensional fourth-order moment (kurtosis) of the velocity derivative in a turbulent field is made under the assumption that the values of kurtosis depend on both the turbulence Reynolds number and the intermittency factor. The method consists of modeling a suitable probability density of the variable in a given turbulence Reynolds number and the intermittency factor. A crude model of the probability density function is derived, and the numerical calculations based on the model show excellent agreement with many of the experimental data. The analysis shows that the values of kurtosis depend strongly on the intermittency factor, and that depending on the value of the intermittency factor, it is entirely possible to have values of kurtosis as low as five in a flow with a turbulence Reynolds number of 5000.     

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     

The influence of water vapor fluctuations on turbulent fluxes

The influence of water vapor fluctuations on vertical turbulent fluxes is examined. It is shown that effects on density and, consequently, buoyancy are insignificant. However, because specific heat of air is a function of specific humidity, these fluctuations are found to influence sensible heat flux significantly. The critical parameter is the Bowen ratio, and a formulation relating the heat flux assuming dry air to the true value is given. The implications of this analysis to flux-gradient relationships in the surface layer are commented on.     

A CASE STUDY OF NUMERICAL SIMULATION OF SEA FOG ON THE SOUTHERN CHINA COAST

This study uses numerical simulations to examine a case of sea fog that was observed from 20 to 22 March2011 on the southern China coast. The observation dataset includes observatory data, cloud-top temperature from MODIS, GPS sonde, and data from the Integrated Observation Platform for Marine Meteorology(IOPMM). The simulations are based on the Weather Research and Forecasting(WRF) model with four distinct parameter settings.Both the observations and simulations focus on the characteristics of the fog extent, boundary layer structure, and meteorological elements near the air-sea interface. Our main results are as follows:(1) The extent of mesoscale sea fog can be well simulated when the sea surface temperature has at least 0.5 ×0.5 horizontal resolution.(2) To accurately model the vertical structure of the sea fog, particularly the surface-based inversion, vertical levels must be added in the boundary layer.(3) When these model conditions are met, the simulations faithfully reproduce the measured downward shortwave radiation, downward longwave radiation, and surface sensible heat flux during the sea fog period.     

Turbulent Flux Measurements Above and Below the Overstory of a Boreal Aspen Forest

Turbulent flux measurements both above and beneath the canopy of a boreal aspen forest are described. Velocity skewness showed that, beneath the aspen canopy, turbulence was dominated by intermittent, downward penetrating gusts. Eulerian horizontal length scales calculated from integration of the autocorrelation function or spectral peaks were 9.0 and 1.4 times the mean aspen height of 21.5 m respectively. Above-canopy power spectral slopes for all velocity components followed the -2/3 power law, whereas beneath-canopy slopes were closer to -1 and showed a spectral short cut in the horizontal and vertical components. Cospectral patterns were similar both above and beneath the canopy. The Monin–Obukhov similarity function for the vertical wind velocity variance was a well-defined function of atmospheric stability, both above and beneath the canopy. Nocturnal flux underestimation and departures of this similarity function from that expected from Monin–Obukhov theory were a function of friction velocity. Energy balance closure greater than 80% was achieved at friction velocities greater than 0.30 and 0.10 m s-1, above and below the aspen canopy, respectively. Recalculating the latent heat flux using various averaging periods revealed a minimum of 15 min were required to capture 90% of the 30-min flux. Linear detrending reduced the flux at shorter averaging periods compared to block averaging. Lack of energy balance closure and erratic flux behaviour led to the recalculation of the latent and sensible heat fluxes using the ratio of net radiation to the sum of the energy balance terms.     

Flux decomposition into coherent structures

This study examines the intermittency of the momentum flux near the surface and the relation of such intermittency to coherent structures. Toward this goal, variances and covariances are decomposed into coherent structures and less coherent activity. The sampled structures are identified using the Haar transform and then decomposed into eigenvectors of the lagged covariance matrix.The methodology is applied to the momentum flux for a relatively stationary 50-h period of strong winds measured from a 45 m tower in the Lammefjord Experiment. Events of sinking motion with strong horizontal momentum account for the majority of the flux. Such sweeping motions arrive as gust microfronts. The large momentum flux is associated with strong coherent fluctuations of the longitudinal wind component and high correlation with relatively modest fluctuations of vertical motion. In the heated case (HAPEX), a phase lag between the vertical and horizontal velocity fluctuations leads to less efficient momentum transport by the main coherent structures.The event nature of the flux is used to formulate an expression for the flux error due to sampling problems. Estimation of the momentum flux requires a significantly longer record than for the heat flux. Modulation of the flux by mesoscale variations also affects the sampling strategy.