BOUNDARY LAYER GROWTH AND LAPSE RATE CHANGES DETERMINED BY LIDAR AND SURFACE HEAT FLUX IN SOFIA

In this study the results from a boundary layer experiment,conducted in autumn 1991 over a flat,build-up urbanarea in Southeast Sofia,together with some models for mixed layer growth rates are used to investigate the layered struc-ture of the vertical atmospheric stability distribution in the Sofia Valley.Lidar measurements of aerosol layer heightsand morning boundary layer development are combined with surface eddy correlation measurements of kinematic heatand moisture fluxes,profiles of temperature and humidity,wind speed and wind direction.A diagnostic method is pres-ented for determining vertical lapse rates using surface meteorological measurements and lidar returns observed duringthe transition from nighttime stable stratification to daytime convective boundary layer after the sunrise.     

BOUNDARY LAYER GROWTH AND LAPSE RATE CHANGES DETERMINED BY LIDAR AND SURFACE HEAT FLUX IN SOFIA*

In this study the results from a boundary layer experiment,conducted in autumn 1991 over a flat,build-up urban area in Southeast Sofia,together with some models for mixed layer growth rates are used to investigate the layered structure of the vertical atmospheric stability distribution in the Sofia Valley.Lidar measurements of aerosol layer heights and morning boundary layer development are combined with surface eddy correlation measurements of kinematic heat and moisture fluxes,profiles of temperature and humidity,wind speed and wind direction.A diagnostic method is presented for determining vertical lapse rates using surface meteorological measurements and lidar returns observed during the transition from nighttime stable stratification to daytime convective boundary layer after the sunrise.     

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.     

Surface and Mixed-Layer Variance Methods to Estimate Regional Sensible Heat Flux at the Surface

Turbulence measurements in the lower half of the convective boundary layer (CBL), which includes both mixed layer and surface layer, were carried out with five sonic anemometers mounted on a 213-m tower over a complex flat suburban area with patches of forest, agricultural land, houses and buildings. Also made were radiosoundings of temperature, humidity and wind speed, to determine the CBL height. The sonic anemometer data of wind speed and temperature were processed to derive the second-moment turbulent statistics and were analyzed to investigate the applicability of variance methods to estimate regional surface fluxes of sensible heat. It was found that the temperature variances in the lower mixed layer, coupled with universal functions, produced sensible heat fluxes H over the area with an rms error of the order of 40 Wm^-2 when compared with H values derived from the eddy correlation method. The variance of the vertical wind speed did not produce as good a result. In contrast, the surface-layer temperature variances yielded H values with rms error of the order of 20 Wm^-2, even though the underlying surface was non-uniform and highly non-isothermal, above which enhanced temperature variances could be suspected.     

ESTIMATES OF SEA SURFACE FLUXES DURING TOGA-COARE IOP

With the marine ship observation data set obtained by three cruise-phases of Chinese Xi-angyanghong 5 from November 5,1992 to February 19,1993 in the TOGA-COARE IOP at 2°S,156°E.the sea surface fluxes at this point are estimated by three different bulk schemes.Firstly.aquasi-linear relation is found between the neutral drag coefficients and wind speed.Then,the sta-bility-dependent drag and heat transfer coefficients are solved in the iterative method.Based onthus-derived transfer coefficients,the momentum,sensible and latent heat fluxes are calculated.In the warm pool region,the fluxes corresponding to the westeily winds are much greater thanthose of the easterly trade winds.The magnitude of sea surface fluxes depends upon the atmo-spheric stability as well,in particular in the case of weak wind condition.The estimated sea sur-face net heat budget shows that considerable amount of heat transport from ocean to atmosphere ismainly produced by the effective longwave radiation,latent and sensible heat fluxes.Among themthe value of latent heat flux is the largest and the sensible heat flux is the smallest.Finally,an ac-curacy analysis is made by direct measurements with the eddy-correlation method on the JapaneseR/V Hakuho board at the same time.It is shown that the bulk-derived fluxes are acceptable withmuch confidence.The estimated effective longwave radiation is used to compare with those by di-rect observations on Xiangyanghong 5.This research is compared with TOGA investigation in thetropical western Pacific.     

DERIVATION OF FLUXES FROM PROFILES OVER A MODERATELY HOMOGENEOUS FOREST

Measurements of fluxes and profiles of wind andtemperature are performed in the roughness layer ofa moderately homogeneous forest location. Weinvestigate to what extent vertical scalar fluxescan be derived from profile measurements. Theinfluence of inhomogeneities in the upwind terrainis investigated with footprint analysis and with aninhomogeneous surface-layer model. Four methods toestimate displacement height are suggested, amongthem is a method involving the structure parameterof the vertical wind. All methods give a decrease ofdisplacement height with increasing wind speed,while roughness length is found to increase withincreasing wind speed. For near-neutral conditionsdimensionless temperature gradients are found to besubstantially lower than the surface-layer valuesfound in the literature for homogeneous terrain with lowvegetation. Dimensionless shear however iscomparable with the surface-layer value. The heightof the roughness layer is 20 times the roughnesslength. Two schemes with locally derived surfaceparameters are tested to derive friction velocityand sensible heat flux from the profilemeasurements. These site specific schemes performsatisfactorily. A third scheme based on surface parameters chosen a priorifrom the literatureperforms significantly worse especially for low windspeed and unstable cases.     

An Optimal Inverse Method Using Doppler Lidar Measurements to Estimate the Surface Sensible Heat Flux

Inverse methods are widely used in various fields of atmospheric science. However, such methods are not commonly used within the boundary-layer community, where robust observations of surface fluxes are a particular concern. We present a new technique for deriving surface sensible heat fluxes from boundary-layer turbulence observations using an inverse method. Doppler lidar observations of vertical velocity variance are combined with two well-known mixed-layer scaling forward models for a convective boundary layer (CBL). The inverse method is validated using large-eddy simulations of a CBL with increasing wind speed. The majority of the estimated heat fluxes agree within error with the proscribed heat flux, across all wind speeds tested. The method is then applied to Doppler lidar data from the Chilbolton Observatory, UK. Heat fluxes are compared with those from a mast-mounted sonic anemometer. Errors in estimated heat fluxes are on average 18 %, an improvement on previous techniques. However, a significant negative bias is observed (on average $-63\,\%$ ) that is more pronounced in the morning. Results are improved for the fully-developed CBL later in the day, which suggests that the bias is largely related to the choice of forward model, which is kept deliberately simple for this study. Overall, the inverse method provided reasonable flux estimates for the simple case of a CBL. Results shown here demonstrate that this method has promise in utilizing ground-based remote sensing to derive surface fluxes. Extension of the method is relatively straight-forward, and could include more complex forward models, or other measurements.     

Ship-board measurements and estimations of air-sea fluxes in the western tropical pacific during TOGA COARE

Direct air-sea flux measurements were made on R/V Kexue #1 at 4 ° S, 156 ° E during the Tropical Ocean Global Atmosphere (TOGA) Coupled Ocean-Atmospheric Response Experiment (COARE) Intensive Observation Period (IOP). An array of six accelerometers was used to measure the motion of the anchored ship, and a sonic anemometer and Lyman-α hygrometer were used to measure the turbulent wind vector and specific humidity. The contamination of the turbulent wind components by ship motion was largely removed by an improvement of a procedure due to Shao based on the acceleration signals. The scheme of the wind correction for ship motion is briefly outlined. Results are presented from data for the best wind direction relative to the ship to minimize flow distortion effects. Both the time series and the power spectra of the sonic-measured wind components show swell-induced ship motion contamination, which is largely removed by the accelerometer correction scheme. There was less contamination in the longitudinal wind component than in the vertical and transverse components. The spectral characteristics of the surface-layer turbulence properties are compared with those from previous land and ocean results. Momentum and latent heat fluxes were calculated by eddy correlation and compared to those estimated by the inertial dissipation method and the TOGA COARE bulk formula. The estimations of wind stress determined by eddy correlation are smaller than those from the TOGA COARE bulk formula, especially for higher wind speeds, while those from the bulk formula and inertial dissipation technique are generally in agreement. The estimations of latent heat flux from the three different methods are in reasonable agreement. The effect of the correction for ship motion on latent heat fluxes is not as large as on momentum fluxes.     

Quality analysis applied on eddy covariance measurements at complex forest sites using footprint modelling

Summary Measuring turbulent fluxes with the eddy covariance method has become a widely accepted and powerful tool for the determination of long term data sets for the exchange of momentum, sensible and latent heat, and trace gases such as CO₂ between the atmosphere and the underlying surface. Several flux networks developed continuous measurements above complex terrain, e.g. AmeriFlux and EUROFLUX, with a strong focus on the net exchange of CO₂ between the atmosphere and the underlying surface. Under many conditions basic assumptions for the eddy covariance method in its simplified form, such as stationarity of the flow, homogeneity of the surface and fully developed turbulence of the flow field, are not fulfilled. To deal with non-ideal conditions which are common at many FLUXNET sites, quality tests have been developed to check if these basic theoretical assumptions are valid.In the framework of the CARBOEUROFLUX project, we combined quality tests described by Foken and Wichura (1996) with the analytical footprint model of Schmid (1997). The aim was to identify suitable wind sectors and meteorological conditions for flux measurements. These tools were used on data of 18 participating sites. Quality tests were applied on the fluxes of momentum, sensible and latent heat, and on the CO₂-flux, respectively. The influence of the topography on the vertical wind component was also checked. At many sites the land use around the flux towers is not homogeneous or the fetch may not be large enough. So the relative contribution of the land use type intended to be measured was also investigated. Thus the developed tool allows comparative investigations of the measured turbulent fluxes at different sites if using the same technique and algorithms for the determination of the fluxes as well as analyses of potential problems caused by influences of the surrounding land use patterns.     

论塔层风 、温廓线

本文改进了Zilitinkevich的工作,得到在塔层内动力学上合理的风廓线.并近似推求了不稳定层结下边界层高度h的表达式,使风、温廓线能用于不同层结.由近地层理论从近地层风、温求出通量后,即可推求塔层风、温分布.325m气象塔资料证明这一廓线达到一定的精度.     

Estimating heat fluxes by merging profile formulae and the energy budget with a variational technique

A variational technique (VT) is applied to estimate surface sensible and latent heat fluxes based on observations of air temperature, wind speed, and humidity, respectively, at three heights (1 m, 4 m, and 10 m), and the surface energy and radiation budgets by the surface energy and radiation system (SERBS). The method fully uses all information provided by the measurements of air temperature, wind, and humidity profiles, the surface energy budget, and the similarity profile formulae as well. Data collected at Feixi experiment station installed by the China Heavy Rain Experiment and Study (HeRES) Program are used to test the method. Results show that the proposed technique can overcome the well-known unstablility problem that occurs when the Bowen method becomes singular; in comparison with the profile method, it reduces both the sensitivities of latent heat fluxes to observational errors in humidity and those of sensible heat fluxes to observational errors in temperature, while the estimated heat fluxes approximately satisfy the surface energy budget. Therefore, the variational technique is more reliable and stable than the two conventional methods in estimating surface sensible and latent heat fluxes.     

1998年南海季风爆发期间近海面层大气湍流结构和通量输送的观测研究

1998年5月14日至6月22日,在西沙永兴岛近海铁塔上进行了一次海-气通量观测试验,观测期包括了西南季风爆发前、爆发、爆发后风速加强等几个阶段。这次试验获得该年西南季风爆发的天气特征以及由涡度相关法、廓线法计算的动量、感热、潜热通量及湍流强度等一些统计量分布。分析结果表明,观测期湍流强度σu、σv、σw与平均风速之比为0.096、0.066、0.045;在近中性条件下（z/L≈0）,各相似函数基本为常数,σu/u＊≈3,σv/u＊≈2,σw/u＊≈1.25;在稳定条件下u、v、w三方向σ/u＊近似相等;在稳定与不稳定条件下σt/t＊随稳定度参数z/L的变化趋势相似,但符号相反。摩擦速度u＊随风速的变化接近于关系式u＊=0.029U10+0.006（U10为10m高度30min平均风速）,空气动力粗糙度长度z0变化在0.01～0.35mm之间。观测期动量通量变化在0.05～0.30N/m2之间,季风爆发期明显增大;曳力系数CD平均为1.12×10-3,它随平均风速变化可以表示为关系式103CD=0.003U210+0.020U10+0.836或103CD=0.056U10+0.732。感热通量由海洋输向大气,平均值为7.8W/m2,数值在0～15W/m2之间变化,季风爆发后明显增大;潜热通量数值一般变化于50～200W/m2,白天中午数值较高,夜间较低。季风爆发后也明显增大。观测期间鲍恩比（BowenRatio）日平均为0.05左右。     

Airborne surveys of the atmospheric boundary layer above the marginal ice zone on the Newfoundland shelf

Abstract

Airborne measurements in the atmospheric boundary layer (ABL) above the marginal ice zone (MIZ) on the Newfoundland Shelf reveal strong lateral variations in mean wind, temperature and the vertical fluxes of heat and momentum under conditions of cold, off‐ice wind. Flux measurements in (and near) the surface layer indicate that the neutral 10‐m drag coefficient depends on ice concentration, ranging from 2 × 10^‐3 at 10% coverage to 5 × 10^‐3 at 90%. Furthermore, cross‐ice‐edge transects consistently show increasing wind speed, temperature and heat flux in the off‐ice direction, but the momentum flux may either increase or decrease, depending on the relative importance of surface buoyancy flux and roughness. For the conditions encountered in this experiment, it appears surface wave maturity does not have a significant influence on the drag coefficient in fetch‐limited regimes near the ice edge.     

How Well Can We Measure the Vertical Wind Speed? Implications for Fluxes of Energy and Mass

Sonic anemometers are capable of measuring the wind speed in all three dimensions at high frequencies (10–50 Hz), and are relied upon to estimate eddy-covariance-based fluxes of mass and energy over a wide variety of surfaces and ecosystems. In this study, wind-velocity measurement errors from a three-dimensional sonic anemometer with a non-orthogonal transducer orientation were estimated for over 100 combinations of angle-of-attack and wind direction using a novel technique to measure the true angle-of-attack and wind speed within the turbulent atmospheric surface layer. Corrections to the vertical wind speed varied from −5 to 37% for all angles-of-attack and wind directions examined. When applied to eddy-covariance data from three NOAA flux sites, the wind-velocity corrections increased the magnitude of CO₂ fluxes, sensible heat fluxes, and latent heat fluxes by ≈11%, with the actual magnitude of flux corrections dependent upon sonic anemometer, surface type, and scalar. A sonic anemometer that uses vertically aligned transducers to measure the vertical wind speed was also tested at four angles-of-attack, and corrections to the vertical wind speed measured using this anemometer were within ±1% of zero. Sensible heat fluxes over a forest canopy measured using this anemometer were 15% greater than sensible heat fluxes measured using a sonic anemometer with a non-orthogonal transducer orientation. These results indicate that sensors with a non-orthogonal transducer orientation, which includes the majority of the research-grade three-dimensional sonic anemometers currently in use, should be redesigned to minimize sine errors by measuring the vertical wind speed using one pair of vertically aligned transducers.     

Acoustic tomography and conventional meteorological measurements over heterogeneous surfaces

Summary Conventional micro-meteorological measuring methods are not particularly suitable for the investigation of the energy exchange under heterogeneous surface conditions. To consider the influence of the different surface properties, area-covered and spatially averaged meteorological measurements in combination with highly resolved simulations are necessary. In this context, the method of acoustic travel time tomography is introduced to provide information about the horizontal temperature and wind field.To check the applicability of the tomographic method for this problem and to provide data for the model initialisation a field experiment was carried out.The tomographic system and conventional meteorological equipment were utilised inside the investigation area with an extension of several hundred meters (300×700m²), which was arranged over areas with different surface properties: grassland and bare soil.The results of the field experiment show that differences between the measuring systems (in-situ and remote sensing) exist and the heterogeneity of the underlying surface is visible in the near surface temperature and wind field at a scale, which can be resolved with highly resolved numerical models. Depending on the incoming solar radiation and the local advection regional distinctions in the air temperature and wind field as well as in the vertical sensible heat fluxes were observed.The investigations demonstrate that the sensitivity of the acoustic tomography is sufficient to verify gradients in the meteorological fields even when the horizontal differences in sensible heat flux are small.     

Water vapour flux profiles in the convective boundary layer

Summary Water vapour flux profiles in the atmospheric boundary layer have been derived from measurements of water vapour density fluctuations by a ground-based Differential Absorption Lidar (DIAL) and of vertical wind fluctuations by a ground-based Doppler lidar. The data were collected during the field experiment LITFASS-2003 in May/June 2003 in the area of Lindenberg, Germany. The eddy-correlation method was applied, and error estimates of ±50 W/m² for latent heat flux were found. Since the sampling error dominates the overall measurement accuracy, time intervals between 60 and 120 min were required for a reliable flux calculation, depending on wind speed. Rather large errors may occur with low wind speed because the diurnal cycle restricts the useful interval length. In the lower height range, these measurements are compared with DIAL/radar-RASS fluxes. The agreement is good when comparing covariance and error values. The lidar flux profiles are well complemented by tower measurements at 50 and 90 m above ground and by area-averaged near surface fluxes from a network of micrometeorological stations. Water vapour flux profiles in the convective boundary layer exhibit different structures mainly depending on the magnitude of the entrainment flux. In situations with dry air above the boundary layer a positive entrainment flux is observed which can even exceed the surface flux. Flux profiles which linearly increase from the surface to the top of the boundary layer are observed as well as profiles which decrease in the lower part and increase in the upper part of the boundary layer. In situations with humid air above the boundary layer the entrainment flux is about zero in the upper part of the boundary layer and the profiles in most cases show a linear decrease.     

Turbulent characteristics and bulk transfer coefficients over the desert in the HEIFE area

Observations of surface-layer turbulence and turbulent fluxes were made over a desert in northwestern China as a part of HEIFE (HEIhe river Field Experiment). These show that the normalized variations of the vertical wind component and of the air temperature obey Monin-Obukhov similarity well, especially in free convective conditions. However, the variations of specific humidity do not obey Monin-Obukhov similarity.Mean bulk transfer coefficients of sensible heat and momentum flux are obtained as functions of stability over a wide stability range from the observed data of turbulent fluxes and mast profiles. However, the bulk transfer coefficient for water vapor could not be obtained because of the large scatter of the data. In free convective conditions, the sensible heat flux was found to be approximately proportional to the 1.4 power of temperature difference between the surface and 20m. The bulk transfer coefficient of sensible heat is also obtained as a function of the bulk Richardson number for practical convenience.     

Cross‐shore variations of near‐surface wind velocity and atmospheric turbulence at the land‐sea boundary during CASP

Abstract

Airborne measurements of mean wind velocity and turbulence in the atmospheric boundary layer under wintertime conditions of cold offshore advection suggest that at a height of 50 m the mean wind speed increases with offshore distance by roughly 20% over a horizontal scale of order 10 km. Similarly, the vertical gust velocity and turbulent kinetic energy decay on scales of order 3.5 km by factors of 1.5 and 3.2, respectively. The scale of cross‐shore variations in the vertical fluxes of heat and downwind momentum is also 10 km, and the momentum flux is found to be roughly constant to 300 m, whereas the heat flux decreases with height. The stability parameter, z/L (where z = 50 m and L is the local Monin‐Obukhov length), is generally small over land but may reach order one over the warm ocean. The magnitude and horizontal length scales associated with the offshore variations in wind speed and turbulence are reasonably consistent with model results for a simple roughness change, but a more sophisticated model is required to interpret the combined effects of surface roughness and heat flux contrasts between land and sea.

Comparisons between aircraft and profile‐adjusted surface measurements of wind speed indicate that Doppler biases of 1–2 m s^?1 in the aircraft data caused by surface motions must be accounted for. In addition, the wind direction measurements of the Minimet anemometer buoy deployed in CASP are found to be in error by 25 ± 5°, possibly due to a misalignment of the anemometer vane. The vertical fluxes of heat and momentum show reasonably good agreement with surface estimates based on the Minimet data.     

A METHOD FOR DETERMINING TURBULENT TRANSFER IN THE ATMOSPHERIC SURFACE LAYER

Derivation of bulk transport coefficients helps solving land surface processes. A similarity-based method for determining the turbulent transfer (including the flux exchange, the vertical distribution of wind and potential temperature) in the atmospheric surface layer is presented. Comparisons with iterative schemes (Businger, 1971) are given to demonstrate the advantages of the calculation methods.     

Horizontal and Vertical Turbulent Fluxes Forced by a Gravity Wave Event in the Nocturnal Atmospheric Surface Layer Over the Amazon Forest

A nocturnal gravity wave was detected over a south-western Amazon forest during the Large-Scale Biosphere–Atmosphere experiment in Amazonia (LBA) in the course of the dry-to-wet season campaign on October 2002. The atmospheric surface layer was stably stratified and had low turbulence activity, based on friction velocity values. However, the passage of the wave, an event with a period of about 180–300 s, caused negative turbulent fluxes of carbon dioxide (CO₂) and positive sensible heat fluxes, as measured by the eddy-covariance system at 60 m (≈30 m above the tree tops). The evolution of vertical profiles of air temperature, specific humidity and wind speed during the wave movement revealed that cold and drier air occupied the sub-canopy space while high wind speeds were measured above the vegetation. The analysis of wind speed and scalars high frequency data was performed using the wavelet technique, which enables the decomposition of signals in several frequencies allowed by the data sampling conditions. The results showed that the time series of vertical velocity and air temperature were −90° out of phase during the passage of the wave, implying no direct vertical transport of heat. Similarly, the time series of vertical velocity and CO₂ concentration were 90° out of phase. The wave was not directly associated with vertical fluxes of this variable but the mixing induced by its passage resulted in significant exchanges in smaller scales as measured by the eddy-covariance system. The phase differences between horizontal velocity and both air temperature and CO₂ concentration were, respectively, zero and 180°, implying phase and anti-phase relationships. As a result, the wave contributed to positive horizontal fluxes of heat and negative horizontal fluxes of carbon dioxide. Such results have to be considered in nocturnal boundary-layer surface-atmosphere exchange schemes for modelling purposes.