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.     

A Comment on “How Well Can We Measure the Vertical Wind Speed? Implications for Fluxes of Energy and Mass” by Kochendorfer et al.

Kochendorfer et al. (Boundary-Layer Meterol, 145:383–398, 2012) conducted an experiment to evaluate azimuth and angle-of-attack dependent errors of sonic anemometer measurements. Several questions are raised regarding the experimental design and the presented results. The finding that instruments with non-orthogonal sonic paths underestimate fluctuations of vertical wind speed and consequently also scalar fluxes by about 10 % is compared with the results of a hitherto unpublished side-by-side field comparison and other past intercomparison experiments. Scale considerations are presented that raise considerable doubts on the validity of the implicit assumption of Kochendorfer et al. (2012) that the turbulent wind vector is highly correlated across a distance of 1.2 m at a height of 2.5 m over flat grassland, which corresponds to the separation between the sonic anemometers tested in their experiment. Nevertheless, new developments in sonic anemometer design to minimize transducer-shadow effects are desirable.     

On drag coefficient parameterization with post processed direct fluxes measurements over the ocean

This study presents an evaluation of the atmospheric factors influencing the post-processing for fast-response data of horizontal momentum, vertical wind component, temperature, and water vapor to measure turbulent fluxes. They are observed at the Ieodo ocean research station over the Yellow Sea during the period of October 2004 to February 2008. The post process methods employed here are composed of quality control and tilt correction for turbulent flux measurement. The present result of quality control on the fast-response data shows that total removal ratio of the data generally depends on the factors such as a wind speed, relative humidity, significant wave height, visibility, and stability parameter (z/L). Especially, the removal ratio of water vapor data is significantly increased on light wind and strong stability conditions. The results show that the total removal ratio of water vapor data increases when wind speed is less than 3 m s^?1 and wave height is less than 1 m. The total removal ratio of water vapor data also increases with the value of the stability parameter. Three different algorithms of tilt correction methods (double rotation, triple rotation, and planar fit) are applied to correct the tilt of the sonic anemometer used in the observation. Friction velocities in near neutral state are greater than friction velocity in other states. Drag coefficients are categorized in terms of stabilities and seasons.     

Large-Scale Turbulence Structures and Their Contributions to the Momentum Flux and Turbulence in the Near-Neutral Atmospheric Boundary Layer Observed from a 213-m Tall Meteorological Tower

Large-scale turbulence structures in the near-neutral atmospheric boundary layer (ABL) are investigated on the basis of observations made from the 213-m tall meteorological tower at Tsukuba, Japan. Vertical profiles of wind speed and turbulent fluxes in the ABL were obtained with sonic anemometer-thermometers at six levels of the tower. From the archived data, 31 near-neutral cases are selected for the analysis of turbulence structures. For the typical case, event detection by the integral wavelet transform with a large time scale (180 s) from the streamwise velocity component (u) at the highest level (200 m) reveals a descending high-speed structure with a time scale of approximately 100 s (a spatial scale of 1 km at the 200-m height). By applying the wavelet transform to the u velocity component at each level, the intermittent appearance of large-scale high-speed structures extending also in the vertical is detected. These structures usually make a large contribution to the downward momentum transfer and induce the enhancement of turbulent kinetic energy. This behaviour is like that of “active” turbulent motions. From the analysis of the two-point space–time correlation of wavelet coefficients for the u velocity component, the vertical extent and the downward influence of large-scale structures are examined. Large fluctuations in the large-scale range (wavelet variance at the selected time scale) at the 200-m level tend to induce the large correlation between the higher and lower levels.     

Evaluation of the Interpretation of Ceilometer Data with RASS and Radiosonde Data

Since 2006 different remote monitoring methods for determining mixing-layer height have been operated in parallel in Augsburg (Germany). One method is based on the operation of eye-safe commercial mini-lidar systems (ceilometers). The optical backscatter intensities recorded with ceilometers provide information about the range-dependent aerosol concentration; gradient minima within this profile mark the tops of mixed layers. Special software for these ceilometers provides routine retrievals of lower atmospheric layering. A second method, based on sodar observations, detects the height of a turbulent layer characterized by high acoustic backscatter intensities due to thermal fluctuations and a high variance of the vertical velocity component. This information is extended by measurements with a radio-acoustic sounding system (RASS) that directly provides the vertical temperature profile from the detection of acoustic signal propagation and thus temperature inversions that mark atmospheric layers. Ceilometer backscatter information is evaluated by comparison with parallel measurements. Data are presented from 2 years of combined ceilometer and RASS measurements at the same site and from comparison with a nearby (60 km) radiosonde for larger-scale humidity information. This evaluation is designed to ensure mixing-layer height monitoring from ceilometer data more reliable.     

Observations of a multi-level turbulence structure in a very stable atmospheric boundary layer

Boundary-layer measurements conducted at the Marsta site in Sweden from a winter-time situation (23–25 Feb.) with stable stratification have been analysed. The data comprise wind and temperature profile measurements up to 30 m, turbulence measurements at 2, 6 and 30 m and Doppler acoustic sounder data up to about 150 m. The upwind fetch at the site is flat and free from obstacles to a distance of ca 5 km for the particular sector chosen for the experiment.During the night, a two-layer vertical structure developed. Analysis of power spectra, co-spectra and variances in a shallow and very stable turbulent boundary layer near the ground show that the turbulence is fully developed and follow the universal behaviour.Above, at a height of 30 m, another turbulent layer is produced by increased wind shear near a low-level jet. This turbulent upper layer can be regarded as a layer of free shear flow. At this height, there also exist wave-turbulence interactions at low frequencies which sometimes cause a countergradient heat flux.     

Estimating the components of the sensible heat budget of a tall forest canopy in complex terrain

Ultrasonic wind measurements, sonic temperature and air temperature data at two heights in the advection experiment MORE II were used to establish a complete budget of sensible heat including vertical advection, horizontal advection and horizontal turbulent flux divergence. MORE II took place at the long-term Carbo-Europe IP site in Tharandt, Germany. During the growing period of 2003 three additional towers were established to measure all relevant parameters for an estimation of advective fluxes, primarily of CO₂. Additionally, in relation to other advection experiments, a calculation of the horizontal turbulent flux divergence is proposed and the relation of this flux to atmospheric stability and friction velocity is discussed. In order to obtain a complete budget, different scaling heights for horizontal advection and horizontal turbulent flux divergence are tested. It is shown that neglecting advective fluxes may lead to incorrect results. If advective fluxes are taken into account, the sensible heat budget based upon vertical turbulent flux and storage change only, is reduced by approximately 30%. Additional consideration of horizontal turbulent flux divergence would in turn add 5–10% to this sum (i.e., the sum of vertical turbulent flux plus storage change plus horizontal and vertical advection). In comparison with available energy horizontal advection is important at night whilst horizontal turbulent flux divergence is rather insignificant. Obviously, advective fluxes typically improve poor nighttime energy budget closure and might change ecosystem respiration fluxes considerably.     

Triple-Hot-Film Anemometer Performance in Cases-99 and A Comparison With Sonic Anemometer Measurements

Two levels of triple-hot-film and sonic anemometers were deployed on a 5.5-m towerduring the Cooperative Atmospheric Surface Exchange Study (CASES-99) in October1999. Each triple-hot-film probe was collocated 50 mm from the sonic sensing path ona common boom. Various problems with using triple-hot-films in the atmosphere toresolve wind components are addressed including the derivation of a yaw angle correction using the collocated sensors. It was found that output voltage drift due to changes in environmental temperature could be monitored and corrected using an automated system. Non-unique solutions to heat transfer equations can be resolved using a collocated sonic anemometer. Multi-resolution decomposition of the hot-film data was used to estimate appropriate day and night averaging periods for turbulent flux measurements in and near the roughness sub-layer. Finally, triple-hot-film measurements of mean wind magnitude (M), turbulent kinetic energy (TKE), sensible heat flux (H), and local friction velocity (u_*) are compared to those of the collocated CSAT3 sonic anemometers. Overall, the mean wind magnitudes measured by the triple-hot-film and the collocated sonic sensorswere close, consistent and independent of stability or proximity to the ground. The turbulent statistics, TKE, u_*, and H, measured by the two sensor systems were reasonably close together at z = 5 m. However, the ratio of sonic measurement/hot-film measurement decreased toward the ground surface, especially during stable conditions.     

Measurements and Modelling of the Wind Speed Profile in the Marine Atmospheric Boundary Layer

We present measurements from 2006 of the marine wind speed profile at a site located 18 km from the west coast of Denmark in the North Sea. Measurements from mast-mounted cup anemometers up to a height of 45 m are extended to 161 m using LiDAR observations. Atmospheric turbulent flux measurements performed in 2004 with a sonic anemometer are compared to a bulk Richardson number formulation of the atmospheric stability. This is used to classify the LiDAR/cup wind speed profiles into atmospheric stability classes. The observations are compared to a simplified model for the wind speed profile that accounts for the effect of the boundary-layer height. For unstable and neutral atmospheric conditions the boundary-layer height could be neglected, whereas for stable conditions it is comparable to the measuring heights and therefore essential to include. It is interesting to note that, although it is derived from a different physical approach, the simplified wind speed profile conforms to the traditional expressions of the surface layer when the effect of the boundary-layer height is neglected.     

Modification of air flow over an escarpment — Results from the Hjardemål experiment

Mean flow, turbulence, and surface pressure measurements over an escarpment are presented. The speed-up in the mean wind field shows the known dependence on stratification. Cross-sections of the standard deviation of horizontal and vertical wind components and of the friction velocity are derived from the data and compare favorably with the numerical model of Zeman and Jensen (1987). The modification of turbulent power spectra at intermediate frequencies can be explained by rapid distortion theory. At very low frequencies, there is a quasi-stationary response to the disturbance. Except for speed-up and standard deviations of the wind components, which are also shown for downslope wind, all results in this paper refer to upslope winds.An analysis of the vertical momentum flux reveals that upstream of the escarpment, most of the flux is transported in sweeps of fast, sinking motion to the ground. Downstream of the escarpment, ejections of slow, rising motion dominate the turbulent transport.     

Boundary-layer effects on mountain waves: a new look at some historical studies

Early studies of mountain waves reported various results that have rarely been investigated since. These include: large-amplitude mountain waves above an unstable boundary layer much higher than the mountains; a repeated downwind drift and upwind jump of mountain waves; and larger vertical wind magnitude near sunrise and/or sunset. These are investigated using over 3,000 radiosondes and meso-strato-troposphere (MST) radar. Superadiabatic temperature gradients are found beneath mountain waves, explainable by convection which appears to raise the mountain-wave launching height. Movement of mountain-wave patterns is studied by a new method using height–time vertical wind data. A swaying motion of mountain waves, with period of a few minutes, appears to be equally upwind and downwind, rather than asymmetric at the heights measurable. Also, vertical wind shows no change in mean, variance or extreme values near sunrise and sunset, despite the expected diurnal changes of boundary-layer structure. An explanation for differences between MST radar and other measurements and models of mountain waves is suggested in terms of more than one variety of mountain wave. Type 1 has stable air near the ground; type 2 is above a convective/turbulent boundary layer of significant height as compared to the mountains.     

Observations of Coherent Turbulence Structures in the Near-Neutral Atmospheric Boundary Layer

Turbulence structures of high Reynolds number flow in the near-neutral atmospheric boundary layer (ABL) are investigated based on observations at Shionomisaki and Shigaraki, Japan. A Doppler sodar measured the vertical profiles of winds in the ABL. Using the integral wavelet transform for the time series of surface wind data, the pattern of a descending high-speed structure with large vertical extent (from the surface to more than 200-m level) is depicted from the Doppler sodar data. Essentially this structure is a specific type of coherent structure that has been previously shown in experiments on turbulent boundary-layer flows. Large-scale high-speed structures in the ABL are extracted using a long time scale (240 s) for the wavelet transform. The non-dimensional interval of time between structures is evaluated as 3.0–6.2 in most cases. These structures make a large contribution to downward momentum transfer in the surface layer. Quadrant analyses of the turbulent motion measured by the sonic anemometer (20-m height) suggest that the sweep motion (high-speed downward motion) plays a substantial role in the downward momentum transfer. In general, the contribution of sweep motions to the momentum flux is nearly equal to that of ejection motions (low-speed upward motions). This contribution of sweep motions is related to the large-scale high-speed structures.     

Observations Of Nocturnal Drainage Flow In A Shallow Gully

The formation of cold air drainage flows in a shallow gully is studied during CASES-99 (Cooperative Atmosphere-Surface Exchange Study). Fast and slow response wind and temperature measurements were obtained on an instrumented 10-m tower located in the gully and from a network of thermistors and two-dimensional sonic anemometers, situated across the gully. Gully flow formed on clear nights even with significant synoptic flow. Large variations in surface temperature developed within an hour after sunset and in situ cooling was the dominant factor in wind sheltered locations. The depth of the drainage flow and the height of the down-gully wind speed maximum were found to be largest when the external wind speed above the gully flow is less than 2 m s^-1. The shallow drainage current is restricted to a depth of a few metres, and is deepest when the stratification is stronger and the external flow is weaker. During the night the drainage flow breaks down, sometimes on several occasions, due to intermittent turbulence and downward fluxes of heat and momentum. The near surface temperature may increase by 6 ° C in less than 30 min due to the vertical convergence of downward heat flux. The mixing events are related to acceleration of the flow above the gully flow and decreased Richardson number. These warming events also lead to warming of the near surface soil and reduction of the upward soil heat flux. To examine the relative importance of different physical mechanisms that could contribute to the rapid warming, and to characterize the turbulence generated during the intermittent turbulent periods, the sensible heat budget is analyzed and the behaviour of different turbulent parameters is discussed.     

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.     

Surface energy balance,parameterizations of boundary-layer heights and the application of resistance laws near an Antarctic Ice Shelf front

A study of the surface energy balance with turbulent fluxes obtained by the Monin-Obukhov similarity theory and a comparison with results for resistance laws are presented for the strong baroclinic conditions in the vicinity of the Filchner/Ronne Ice Shelf front. The data are taken from a field experiment in the Antarctic summer season 1983/84. For the first time in the coastal Antarctic region, this data set comprises synchronous energy balance measurements over the polynya and the ice shelf together with soundings of the boundary layer, yielding vertical profiles of the wind velocity and temperature over the ice shelf, at the ice shelf front and over the polynya.Over the ice shelf, the radiation balance is the largest component of the energy fluxes and is mainly compensated by the subsurface energy flux and the turbulent heat flux in the daily mean. Over the polynya, turbulent fluxes of sensible and latent heat lead to large energy losses of the water surface in the night-time and in situations of very low air temperatures.Different parameterizations for boundary-layer height are compared using tethered sonde and energy balance measurements. With the height of the inversion base over the polynya and the height of the critical bulk Richardson number over the ice shelf, external parameters for the application of resistance laws were determined. The comparison of turbulent surface fluxes obtained by the energy balance measurements and by the resistance laws shows good agreement for the convective conditions over the polynya. For the stably stratified boundary layer over the ice shelf with small amounts of the turbulent heat flux, the deviation is large for the case of a cold air outflow with a superposed inertial oscillation.     

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.     

多普勒声雷达测风与直接测量仪器的观测比较

文中对多普勒声雷达测风和325 m铁塔上风速、风向仪及超声风速仪的测量值进行了观测比较。结果表明,当信噪比大于1时,两种方法测量的平均风速、风向的相关较好,风速、风向廓线较为接近;声雷达测量的垂直风速平均值和均方根差比超声风速仪的测量值略偏大。这与国外进行的观测比较结果一致。 需要指出,声雷达和直接测量仪器的观测比较中存在一定差异,特别是对比每一次的观测值就可以看出,在某些情况下还存在明显差异。进一步分析声雷达测风产生偏差的原因,对误差做出合理的订正是十分必要的。     

理想城市建筑群风场的大涡模拟

为了解城市建筑群对周围风场结构的影响,本文利用大涡模拟方法对不同建筑群分布类型下的风场状况进行了模拟,并利用通风指数对风场进行定量评估。模拟结果表明:PALM(Parallelized LES model)能够反映出建筑物影响下风场的结构特征。建筑群风场受到建筑高度、高度方差、分布、初始入流的影响。风场衰减主要体现在水平分量上,尤其是平均高度较高、高度方差大、交错不规整分布使得水平方向风场衰减、通风效果变差,而同时又会使得垂直风速分量w增大,增大垂直方向的空气交换速率。初始入流不会影响风场风廓线形态,在低层风场中平均风速大小以及总体空气交换速率基本与入流成正比关系。     

The Høvsøre Tall Wind-Profile Experiment: A Description of Wind Profile Observations in the Atmospheric Boundary Layer

We present an analysis of data from a nearly 1-year measurement campaign performed at Høvsøre, Denmark, a coastal farmland area where the terrain is flat. Within the easterly sector upstream of the site, the terrain is nearly homogenous. This topography and conditions provide a good basis for the analysis of vertical wind-speed profiles under a wide range of atmospheric stability, turbulence, and forcing conditions. One of the objectives of the campaign was to serve as a benchmark for flow over flat terrain models. The observations consist of combined wind lidar and sonic anemometer measurements at a meteorological mast. The sonic measurements cover the first 100 m and the wind lidar measures above 100 m every 50 m in the vertical. Results of the analysis of observations of the horizontal wind-speed components in the range 10–1200 m and surface turbulence fluxes are illustrated in detail, combined with forcing conditions derived from mesoscale model simulations. Ten different cases are presented. The observed wind profiles approach well the simulated gradient and geostrophic winds close to the simulated boundary-layer height during both barotropic and baroclinic conditions, respectively, except for a low-level jet case, as expected. The simulated winds are also presented for completeness and show good agreement with the measurements, generally underpredicting the turning of the wind in both barotropic and baroclinic cases.     

北京冬季大风过程大气边界层特征分析

利用北京中国科学院大气物理研究所325 m气象观测塔的气象梯度资料和湍流资料,分析了2014年11月29日至12月5日北京两次大风过程中气象要素和湍流输送特征的变化。第一次大风过程的强度和持续时间均高于第二次大风过程。强烈的风速垂直切变主要集中在距地面100 m高度范围内,最强风速垂直切变达到0.31 s~(-1)。大风过程中,阵风系数呈现随高度减小的趋势,越接近地面,阵风系数愈大。阵风强度的变化与阵风系数相似,100 m以下高度时,阵风强度随高度增大而减小。大风过程自上而下改变边界层结构,平均动能、湍流动能和摩擦速度最先从上层(280 m)发生变化且迅速增加。近地层由于风速垂直梯度的显著差异,近地层垂直方向的湍流强度最大。大风时各功率谱在低频区(0.01 s~(-1))达到峰值,大风过后各高度的能量都有所下降。