A simple unified theory for flow in the canopy and roughness sublayer

The mean flow profile within and above a tall canopy is well known to violate the standard boundary-layer flux–gradient relationships. Here we present a theory for the flow profile that is comprised of a canopy model coupled to a modified surface-layer model. The coupling between the two components and the modifications to the surface-layer profiles are formulated through the mixing layer analogy for the flow at a canopy top. This analogy provides an additional length scale—the vorticity thickness—upon which the flow just above the canopy, within the so-called roughness sublayer, depends. A natural form for the vertical profiles within the roughness sublayer follows that overcomes problems with many earlier forms in the literature. Predictions of the mean flow profiles are shown to match observations over a range of canopy types and stabilities. The unified theory predicts that key parameters, such as the displacement height and roughness length, have a significant dependence on the boundary-layer stability. Assuming one of these parameters a priori leads to the incorrect variation with stability of the others and incorrect predictions of the mean wind speed profile. The roughness sublayer has a greater impact on the mean wind speed in stable than unstable conditions. The presence of a roughness sublayer also allows the surface to exert a greater drag on the boundary layer for an equivalent value of the near-surface wind speed than would otherwise occur. This characteristic would alter predictions of the evolution of the boundary layer and surface states if included within numerical weather prediction models.     

Comparing mixing-length models of the diabatic wind profile over homogeneous terrain

Models of the diabatic wind profile over homogeneous terrain for the entire atmospheric boundary layer are developed using mixing-length theory and are compared to wind speed observations up to 300 m at the National Test Station for Wind Turbines at Høvsøre, Denmark. The measurements are performed within a wide range of atmospheric stability conditions, which allows a comparison of the models with the average wind profile computed in seven stability classes, showing a better agreement than compared to the traditional surface-layer wind profile. The wind profile is measured by combining cup anemometer and lidar observations, showing good agreement at the overlapping heights. The height of the boundary layer, a parameter required for the wind profile models, is estimated under neutral and stable conditions using surface-layer turbulence measurements, and under unstable conditions based on the aerosol backscatter profile from ceilometer observations.     

Surface drag and turbulence over an inhomogeneous land surface

Data collected over an inhomogeneous semi-rural area are presented. The data are compared with previous surface-layer data to determine how representative the fixed-point flux measurements are of area averages. Departures from the standard surface-layer results are found to be relatively small (~10–20%), which supports the concept of a blending height above which the flow ceases to respond to variations in the underlying surface and becomes horizontally homogeneous.Effective roughness lengths are derived for different wind directions and the relationship between the effective roughness length and upwind surface is examined in the light of recent ideas on averaging surface roughness lengths. It is found that by averaging drag coefficients, realistic values of the effective roughness length can be calculated which are not very sensitive to the precise choice of the component roughness lengths.     

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.     

Wind profiles,momentum fluxes and roughness lengths at Cabauw revisited

We describe the results of an experiment focusing on wind speed and momentum fluxes in the atmospheric boundary layer up to 200 m. The measurements were conducted in 1996 at the Cabauw site in the Netherlands. Momentum fluxes are measured using the K-Gill Propeller Vane. Estimates of the roughness length are derived using various techniques from the wind speed and flux measurements, and the observed differences are explained by considering the source area of the meteorological parameters. A clear rough-to-smooth transition is found in the wind speed profiles at Cabauw. The internal boundary layer reaches the lowest k-vane (20 m) only in the south-west direction where the obstacle-free fetch is about 2 km. The internal boundary layer is also reflected in the roughness lengths derived from the wind speed profiles. The lower part of the profile (< 40 m) is not in equilibrium and no reliable roughness analysis can be given. The upper part of the profile can be linked to a large-scale roughness length. Roughness lengths derived from the horizontal wind speed variance and gustiness have large footprints and therefore represent a large-scale average roughness. The drag coefficient is more locally determined but still represents a large-scale roughness length when it is measured above the local internal boundary layer. The roughness length at inhomogeneous sites can therefore be determined best from drag coefficient measurements just above the local internal boundary layers directly, or indirectly from horizontal wind speed variance or gustiness. In addition, the momentum and heat fluxes along the tower are analysed and these show significant variation with height related to stability and possibly surface heterogeneity. It appears that the dimensionless wind speed gradients scale well with local fluxes for the variety of conditions considered, including the unstable cases.     

Estimation of sector roughness lengths and the effect on prediction of the vertical wind speed profile

An estimate of roughness length is required by some atmospheric models and is also used in the logarithmic profile to determine the increase of wind speed with height under neutral conditions. The choice of technique for determining roughness lengths is generally constrained by the available input data. Here, we compare sets of roughness lengths derived by different methods for the same site and evaluate their impact on the prediction of the vertical wind speed profile.Wind speed and direction data have been collected at four heights over a three-year period at the North Norfolk Wind Monitoring Site. Wind speed profiles were used to generate sector roughness lengths based on the logarithmic profile formula. This is the only direct way of determining roughness lengths. The simplest and cheapest method is to use maps with published tables giving roughness length estimates for different terrain types. Alternatively Wieringa (1976, 1986) and Beljaars (1987) give formulae for determining roughness lengths from wind speed gusts or standard deviations.The four sets of estimated roughness lengths vary considerably. They were used to estimate 34 m wind speeds from 12.7 m observations. The profile-derived roughnesses are used simply as a check on the prediction of the wind speed profiles. The terrain-derived roughness lengths give reasonable results. Gust-derived and standard deviation roughnesses both predict wind speeds which are lower than the observed ones. The error is greater in the case of standard deviation roughnesses. If stability corrections are applied in the prediction of the vertical wind speed profile, the results are considerably improved.     

Intercomparison of Planetary Boundary-Layer Parametrizations in the WRF Model for a Single Day from CASES-99

This study compares five planetary boundary-layer (PBL) parametrizations in the Weather Research and Forecasting (WRF) numerical model for a single day from the Cooperative Atmosphere-Surface Exchange Study (CASES-99) field program. The five schemes include two first-order closure schemes—the Yonsei University (YSU) PBL and Asymmetric Convective Model version 2 (ACM2), and three turbulent kinetic energy (TKE) closure schemes—the Mellor–Yamada–Janjić (MYJ), quasi-normal scale elimination (QNSE), and Bougeault–Lacarrére (BouLac) PBL. The comparison results reveal that discrepancies among thermodynamic surface variables from different schemes are large at daytime, while the variables converge at nighttime with large deviations from those observed. On the other hand, wind components are more divergent at nighttime with significant biases. Regarding PBL structures, a non-local scheme with the entrainment flux proportional to the surface flux is favourable in unstable conditions. In stable conditions, the local TKE closure schemes show better performance. The sensitivity of simulated variables to surface-layer parametrizations is also investigated to assess relative contributions of the surface-layer parametrizations to typical features of each PBL scheme. In the surface layer, temperature and moisture are more strongly influenced by surface-layer formulations than by PBL mixing algorithms in both convective and stable regimes, while wind speed depends on vertical diffusion formulations in the convective regime. Regarding PBL structures, surface-layer formulations only contribute to near-surface variability and then PBL mean properties, whereas shapes of the profiles are determined by PBL mixing algorithms.     

Neutral surface-layer flow over isolated roughness strips

Numerical simulations are used to study neutral surface-layer flow which passes over isolated or widely separated strips of modified roughness embedded in an otherwise homogeneous surface.Simple power laws are given for the maximum height and horizontal extent of turbulent momentum and horizontal mean velocity wakes with an assessment of the range of validity of this proposal. Furthermore, it is indicated, how vertical profiles of energy budget and horizontal mean velocity are distorted by a roughness strip.     

风廓线雷达在高空风场分析中的应用

利用大连风廓线雷达高时空分辨率风场观测资料,统计2011年雷达站上空各层水平及垂直风速的分布特征.通过分析发现:最大水平风速通常出现在12 km上下,受高空急流的影响,各季节高空最大水平风速出现高度不同,4 km以下高空水平风速随高度的变化各月份存在一定差异,4 km以上至最大风速层,水平风速随高度的升高而增大,最大风速层以上至雷达测量的上限水平风速随高度增加先减小后增大;高空垂直风速在夏季较为明显,秋季次之,冬春季节最小;6月是全年月均垂直风速最大的月份,在500～1300 m高度层存在一个上升气流中心,平均风速大于0.6 m/s,2月各高度平均垂直风速全年最小.     

Analysis And Simulation Of Surface-Layer Winds Using Multiplicative Cascade Models With Self-Similar Probability Densities

Statistical analysis techniques based on multiplicative cascades are investigated for use with surface-layer wind data sets collected in the atmospheric boundary layer over flat farm land. The data were found to exhibit multiscaling statistics, allowing the surface-layer winds to be simulated with the use of multiplicative random cascades. The study found evidence that, for the surface-layer at least, these cascade models (andhence the methods of multifractal analysis) should be applied in separate ways to the microscale inertial range, and the mesoscale. This is at odds with the view found in the existing literature, which proposes a `universal multifractal' model to replace the widely held view that there exists separate microscale, mesoscale and synoptic scales for which the processes governing each are different. At least two separate ranges of scaling are suggested for surface-layer wind data, corresponding to the microscale inertial range and the mesoscale. For the case of the mesoscale range, a self-similar distribution of weighting factors was found for the wind speed data themselves, rather than for an intermediate (dissipation) field, as is required for themicroscale data.     

边界层方案对华北低层O3垂直分布模拟的影响

利用WRF-Chem模式,采用3种边界层参数化方案 (YSU, MYJ和ACM2),针对1个晴空、静稳日 (2013年8月26日20:00—27日20:00(北京时)) 进行模拟,着重分析不同边界层参数化方案对夜间残留层形成及日出前后O₃浓度垂直分布形式的模拟效果,并与固城站地面及垂直同步观测资料进行对比。结果表明:3种边界层参数化方案均能够模拟出温度及风速的区域分布形式以及风温垂直结构的变化特征;相比之下,MYJ方案模拟的夜间边界层高度较YSU方案和ACM2方案明显偏高,该对比结果可能是导致近地面污染物浓度模拟差异的重要原因;在夜间稳定层结至日出后稳定状态打破的边界层结构演变过程中,采用YSU方案和ACM2方案模拟的温度和风速垂直扩线形式与观测结果更为接近;同样采用非局地闭合的YSU方案和同时考虑局地和非局地闭合的ACM2方案,对于边界层高度内O₃浓度垂直分布形式的模拟效果具有明显优势。     

Use of a High-Resolution Sodar to Study Surface-layer Turbulence at Night

Measurements in the atmospheric surface layer are generally made with point sensors located in the first few tens of metres. In most cases, however, these measurements are not representative of the whole surface layer. Standard Doppler sodars allow a continuous display of the turbulent thermal structure and wind profiles in the boundary layer up to 1000 m, with a few points, if any, in the surface layer. To overcome these limitations a new sodar configuration is proposed that allows for a higher resolution in the surface layer. Because of its capabilities (echo recording starting at 2 m, echo intensity vertical resolution of approximately 2 m, temporal resolution of 1 s) this sodar is called the surface-layer mini-sodar (SLM-sodar). Features and capabilities of the SLM-sodar are described and compared with the sodar. The comparison of the thermal vertical structure given by the SLM-sodar and the sodar provides evidence that, in most cases, the surface layer presents a level of complexity comparable to that of the entire boundary layer. Considering its high vertical resolution, the SLM-sodar is a promising system for the study of the nocturnal surface layer. The nocturnal SLM-sodar measurements have shown that, depending on wind speed, the structure of the surface layer may change substantially within a short time period. At night, when the wind speed is greater than 3 m s⁻¹, mechanical mixing destroys the wavy structure present in the nocturnal layer. Sonic anemometer measurements have shown that, in such cases, also the sensible heat flux varies with height, reaching a peak in correspondence with the wind speed peak. Under these conditions the assumption of horizontal homogeneity of the surface layer and the choice of the averaging time need to be carefully treated.     

Sensitivity of Turbine-Height Wind Speeds to Parameters in Planetary Boundary-Layer and Surface-Layer Schemes in the Weather Research and Forecasting Model

We evaluate the sensitivity of simulated turbine-height wind speeds to 26 parameters within the Mellor–Yamada–Nakanishi–Niino (MYNN) planetary boundary-layer scheme and MM5 surface-layer scheme of the Weather Research and Forecasting model over an area of complex terrain. An efficient sampling algorithm and generalized linear model are used to explore the multiple-dimensional parameter space and quantify the parametric sensitivity of simulated turbine-height wind speeds. The results indicate that most of the variability in the ensemble simulations is due to parameters related to the dissipation of turbulent kinetic energy (TKE), Prandtl number, turbulent length scales, surface roughness, and the von Kármán constant. The parameter associated with the TKE dissipation rate is found to be most important, and a larger dissipation rate produces larger hub-height wind speeds. A larger Prandtl number results in smaller nighttime wind speeds. Increasing surface roughness reduces the frequencies of both extremely weak and strong airflows, implying a reduction in the variability of wind speed. All of the above parameters significantly affect the vertical profiles of wind speed and the magnitude of wind shear. The relative contributions of individual parameters are found to be dependent on both the terrain slope and atmospheric stability.     

An Improvement of the Louis Scheme for the Surface Layer in an Atmospheric Modelling System

The Louis scheme for calculating the vertical eddy fluxes within the atmospheric surface layer is improved by broadening the original assumptions. In our approach, the momentum and heat transfer roughness lengths (z₀ and z_T respectively) can be different, and z₀ need not be negligibly small compared with the lowest height (z) in modelling. For these conditions, we choose more consistent wind and potential temperature profile forms, then derive new algorithms for calculating fluxes. Improvement is demonstrated for a wide range of z/L (L is the Obukhov length), z/z₀ and z₀ z_T, by comparing these fluxes with those derived from a theoretical surface-layer model. The improved algorithms can be used in atmospheric modelling systems for more varied surfaces and a wide range of atmospheric stability.     

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.     

Observed Wind Profiles and Turbulence Fluxes over an ice Surface with Changing Surface Roughness

Wind profile and eddy-correlation data obtained at two sites on a melting glacier surface in Iceland during the summer of 1996 are presented. Throughout the experiment the surface roughness increased rapidly from smooth to very rough, with the largest roughness element height obtained being about 1.7 m. In a layer close to the rough surface we find that the wind speed profiles were disturbed showing horizontal inhomogeneities as in a roughness sublayer. Its height was approximately two times the height of the main roughness elements (h) at both sites throughout the experiment. From the wind profiles and eddy-correlation data we calculated corrections for the displaced zero plane as a function of time and compared these with results obtained from a drag partitioning model. In general, the agreement was reasonable considering the ranges of uncertainty but the results indicate that the increasing horizontal anisotropy of the surface probably limits the use of the model. The values obtained for the roughness lengths are in good agreement with those calculated from a simple linear model, i.e., z₀/h = 0.5 with the frontal area index. Above the roughness sublayer the wind profiles, normalised standard deviations of wind speed, and the balance of the turbulence kinetic energy budget behaved as over an ideal homogeneous surface thereby confirming similarity of the flow.     

Weibull-k Revisited: “Tall” Profiles and Height Variation of Wind Statistics

The Weibull distribution is commonly used to describe climatological wind-speed distributions in the atmospheric boundary layer. While vertical profiles of mean wind speed in the atmospheric boundary layer have received significant attention, the variation of the shape of the wind distribution with height is less understood. Previously we derived a probabilistic model based on similarity theory for calculating the effects of stability and planetary boundary-layer depth upon long-term mean wind profiles. However, some applications (e.g. wind energy estimation) require the Weibull shape parameter (k), as well as mean wind speed. Towards the aim of improving predictions of the Weibull- \(k\) profile, we develop expressions for the profile of long-term variance of wind speed, including a method extending our probabilistic wind-profile theory; together these two profiles lead to a profile of Weibull-shape parameter. Further, an alternate model for the vertical profile of Weibull shape parameter is made, improving upon a basis set forth by Wieringa (Boundary-Layer Meteorol, 1989, Vol. 47, 85–110), and connecting with a newly-corrected corollary of the perturbed geostrophic-drag theory of Troen and Petersen (European Wind Atlas, 1989, Risø National Laboratory, Roskilde). Comparing the models for Weibull-k profiles, a new interpretation and explanation is given for the vertical variation of the shape of wind-speed distributions. Results of the modelling are shown for a number of sites, with a discussion of the models’ efficacy and applicability. The latter includes a comparative evaluation of Wieringa-type empirical models and perturbed-geostrophic forms with regard to surface-layer behaviour, as well as for heights where climatological wind-speed variability is not dominated by surface effects.     

The wind field in the nonlinear multilayer planetary boundary layer

By use of the small parameter expansion method, the nonlinear planetary boundary layer (PBL) is studied in this paper. The PBL is divided into the surface layer and the Ekman layer, which is divided into several sublayers. In the surface-layer, the eddy coefficient K is taken as a linear function of height; in the Ekman layer, different constant K values are taken within different sublayers: these values are determined from O'Brien's formula (O'Brien, 1970) approximately. Under the upper and lower boundary conditions and the continuity conditions of the wind velocities and turbulent stresses at each boundary between sublayers, analytical expressions for wind velocity in all sublayers and the vertical velocity at the top of the PBL are obtained. A specific example of steady axisymmetrical circular high and low pressure areas is analysed, and some new conclusions are obtained. The results are in better agreement with reality than previous results. This example also shows that the vertical velocity at the top of the PBL caused by friction approaches zero near the center of a high or low pressure system for this model, but attains its maximum absolute values near the center of the high or low pressure area for Wu's (1984) model. This is due to the fact that in our model, the geostrophic wind speed near the center of this specific vortex approaches zero, which causes the wind shear and the friction effect to be very weak. Therefore the wind distribution in the PBL is very sensitive to the type of eddy coefficient.