The Effect of Stratification on the Aerodynamic Roughness Length and Displacement Height

The roughness length, z _0u , and displacement height, d _0u , characterise the resistance exerted by the roughness elements on turbulent flows and provide a conventional boundary condition for a wide range of turbulent-flow problems. Classical laboratory experiments and theories treat z _0u and d _0u as geometric parameters independent of the characteristics of the flow. In this paper, we demonstrate essential stability dependences—stronger for the roughness length (especially in stable stratification) and weaker but still pronounced for the displacement height. We develop a scaling-analysis model for these dependences and verify it against experimental data.     

Aerodynamic Parameters of a UK City Derived from Morphological Data

Detailed three-dimensional building data and a morphometric model are used to estimate the aerodynamic roughness length z ₀ and displacement height d over a major UK city (Leeds). Firstly, using an adaptive grid, the city is divided into neighbourhood regions that are each of a relatively consistent geometry throughout. Secondly, for each neighbourhood, a number of geometric parameters are calculated. Finally, these are used as input into a morphometric model that considers the influence of height variability to predict aerodynamic roughness length and displacement height. Predictions are compared with estimations made using standard tables of aerodynamic parameters. The comparison suggests that the accuracy of plan-area-density based tables is likely to be limited, and that height-based tables of aerodynamic parameters may be more accurate for UK cities. The displacement heights in the standard tables are shown to be lower than the current predictions. The importance of geometric details in determining z ₀ and d is then explored. Height variability is observed to greatly increase the predicted values. However, building footprint shape only has a significant influence upon the predictions when height variability is not considered. Finally, we develop simple relations to quantify the influence of height variation upon predicted z ₀ and d via the standard deviation of building heights. The difference in these predictions compared to the more complex approach highlights the importance of considering the specific shape of the building-height distributions. Collectively, these results suggest that to accurately predict aerodynamic parameters of real urban areas, height variability must be considered in detail, but it may be acceptable to make simple assumptions about building layout and footprint shape.     

Momentum Transfer By A Mountain Meadow Canopy: A Simulation Analysis Based On Massman's (1997) Model

Using a mountain meadow as a case study it is the objective of the present paper todevelop a simple parameterisation for the within-canopy variation of the phytoelementdrag (C_d) and sheltering (P_m) coefficients required for Massman's model of momentum transfer by vegetation. A constant ratio between C_d and P_m is found to overestimate wind speed in the upper canopy and underestimate it in the lower canopy.Two simple parameterisations of C_d/P_m as a function of the plant area density and the cumulative plant area index are developed, using values optimised by least-squares regression between measured and predicted within-canopy wind speeds. A validation with independently measured data indicates that both parameterisations work reliably for simulating wind speed in the investigated meadow. Model predictions of the normalised zero-plane displacement height and the momentum roughness length fall only partly within the range of values given in literature, which may be explained by the accumulation of plantmatter close to the soil surface specific for the investigated canopies. The seasonal course of the normalised zero-plane displacement height and the momentum roughness length are discussed in terms of the seasonal variation of the amount and density of plant matter.     

A New Aerodynamic Parametrization for Real Urban Surfaces

This study conducted large-eddy simulations (LES) of fully developed turbulent flow within and above explicitly resolved buildings in Tokyo and Nagoya, Japan. The more than 100 LES results, each covering a 1,000 $\times $ 1,000 m $^{2}$ area with 2-m resolution, provide a database of the horizontally-averaged turbulent statistics and surface drag corresponding to various urban morphologies. The vertical profiles of horizontally-averaged wind velocity mostly follow a logarithmic law even for districts with high-rise buildings, allowing estimates of aerodynamic parameters such as displacement height and roughness length using the von Karman constant $=$ 0.4. As an alternative derivation of the aerodynamic parameters, a regression of roughness length and variable Karman constant was also attempted, using a displacement height physically determined as the central height of drag action. Although both the regression methods worked, the former gives larger (smaller) values of displacement height (roughness length) by 20–25 % than the latter. The LES database clearly illustrates the essential difference in bulk flow properties between real urban surfaces and simplified arrays. The vertical profiles of horizontally-averaged momentum flux were influenced by the maximum building height and the standard deviation of building height, as well as conventional geometric parameters such as the average building height, frontal area index, and plane area index. On the basis of these investigations, a new aerodynamic parametrization of roughness length and displacement height in terms of the five geometric parameters described above was empirically proposed. The new parametrizations work well for both real urban morphologies and simplified model geometries.     

Estimating Aerodynamic Parameters of Urban-Like Surfaces with Heterogeneous Building Heights

There are many geometrical factors than can influence the aerodynamic parameters of urban surfaces and hence the vertical wind profiles found above. The knowledge of these parameters has applications in numerous fields, such as dispersion modelling, wind loading calculations, and estimating the wind energy resource at urban locations. Using quasi-empirical modelling, we estimate the dependence of the aerodynamic roughness length and zero-plane displacement for idealized urban surfaces, on the two most significant geometrical characteristics; surface area density and building height variability. A validation of the spatially-averaged, logarithmic wind profiles predicted by the model is carried out, via comparisons with available wind-tunnel and numerical data for arrays of square based blocks of uniform and heterogeneous heights. The model predicts two important properties of the aerodynamic parameters of surfaces of heterogeneous heights that have been suggested by experiments. Firstly, the zero-plane displacement of a heterogeneous array can exceed the surface mean building height significantly. Secondly, the characteristic peak in roughness length with respect to surface area density becomes much softer for heterogeneous arrays compared to uniform arrays, since a variation in building height can prevent a skimming flow regime from occurring. Overall the simple model performs well against available experimental data and may offer more accurate estimates of surface aerodynamic parameters for complex urban surfaces compared to models that do not include height variability.     

Simplified expressions for vegetation roughness length and zero-plane displacement as functions of canopy height and area index

Using a previous treatment of drag and drag partition on rough surfaces, simple analytic expressions are derived for the roughness length (z ₀) and zero-plane displacement (d) of vegetated surfaces, as functions of canopy height (h) and area index (). The resulting expressions provide a good fit to numerous field and wind tunnel data, and are suitable for applications such as surface parameterisations in atmospheric models.     

北京城市下垫面动力学参数的计算与分析

采用北京325 m铁塔2008—2012年的单层超声观测资料,基于莫宁-奥布霍夫相似理论(Monin-Obukhov similarity theory)和前人提出的最小误差分析方法,计算了铁塔周边下垫面的零平面位移高度和动力粗糙度长度。结果表明,由于铁塔位于北京市区,其周边下垫面呈现极其复杂的非均匀性,所以对应铁塔周边不同的扇区,零平面位移高度和动力粗糙度长度各有不同。平均而言,在2008—2012年间,铁塔周边下垫面的零平面位移高度为34.4 m,动力粗糙度长度为1.16 m。此外,综合前人的计算结果发现,铁塔周边的零平面位移高度和动力粗糙度长度在2001年之前呈显著增加的趋势,而在2001年以后并未增长,这一现象与铁塔周边的城市化进程相对应。     

A study of wind speed modification and internal boundary-layer heights in a coastal region

Wind profile data within the first two kilometres of a coast have been used to study the wind field modification downstream of this surface discontinuity. The land area is generally very flat, having an overall roughness length of 0.04 m. A wind model, suitable for practical applications and inexpensive to run, has been tested against the data and was found to give satisfactory results. Knowing the climatological statistics of wind and stratification, e.g., at the coast, the model may thus be used to estimate, on a climatological basis, how the wind field is modified with distance inland, at least in areas with only minor topography. This type of information is of great importance when locating wind turbines. It is in these cases also important to know the statistics of the internal boundary-layer (IBL) height, as the turbulence intensity may be quite different in and above the IBL, which in turn may influence load and fatigue calculations. Using the wind profile data, the IBL height was clearly discernible in the majority of cases. Having very unstable stratification over land, the IBL height could, however, not be determined from the wind profiles, as the wind in these cases did not decrease inland. This result was also obtained using the wind model. A simple model of the type z _IBL = a · x ^b, was instead tested, and was shown to give reasonable results.     

A simple model of drainage flow on a slope

The concept of a cold air ‘Parcel’ is introduced for describing the bulk properties of drainage flow. By means of a model based on the momentum and sensible heat transports under calm conditions, the thickness h and velocity u of the Parcel are derived in simple forms. It is shown that h and u correspond to the inversion height and maximum velocity of actual drainage flow. The governing parameters for h and u are the length and vertical drop of the slope, potential temperature difference between the ambient atmosphere and the Parcel, aerodynamic condition of the slope surface expressed by the mean bulk coefficients, and ambient stability. The mean bulk coefficients depend on the roughness lengths for the velocity and potential temperature profiles and are decreasing functions of the slope length. The Parcel Model agrees qualitatively with Manins and Sawford's (1979) model under neutral ambient stratification. But agreement is not so good under stable conditions. The thickness and velocity of drainage flow predicted by the Parcel Model agree with observations on slopes several tens of meters to several hundred kilometers long.     

Aerodynamic roughness of vegetated surfaces

Available experimental results indicate that as the density of roughness elements over a horizontally homogeneous surface is varied, the roughness length, z ₀, varies in a manner that exhibits a maximum at intermediate density values. In an attempt to explain this behaviour, the available analytical solutions for the wind profile inside dense homogeneous canopies were reviewed. The review indicated that the variation of z ₀ with density depends on the interrelationship between the leaf density, a, and the mixing length, l. In view of this finding, a numerical model was devised based on a simple rule for constructing mixing-length profiles in the canopy. The rule states that the actual value of l is the maximum possible under the two constraints: l l _i and ¦dl/dz¦ k, where k is the von Karman constant and the intrinsic mixing length, l _i, is a function of the local internal structure of the canopy. The model which ensures a smooth transition from dense to thin canopy, was used to reproduce the observed maximum of z ₀. The model is also capable of handling vertically non-homogeneous canopies.     

A comparative study of some mathematical models of the mean wind structure and aerodynamic drag of plant canopies

A semi-analytical method for describing the mean wind profile and shear stress within plant canopies and for estimating the roughness length and the displacement height is presented. This method incorporates density and vertical structure of the canopy and includes simple parameterizations of the roughness sublayer and shelter factor. Some of the wind profiles examined are consistent with first-order closure techniques while others are consistent with second-order closure techniques. Some profiles show a shearless region near the base of the canopy; however, none displays a secondary maximum there. Comparing several different analytical expressions for the canopy wind profile against observations suggests that one particular type of profile (an Airy function which is associated with the triangular foliage surface area density distribution) is superior to the others. Because of the numerical simplicity of the methods outlined, it is suggested that they may be profitably used in large-scale models of plant-atmosphere exchanges.     

A Numerical Study on Effects of Land-Surface Heterogeneity from “Combined Approach” on Atmospheric Process Part I: Principle and Method

A method based on Giorgi (1997a, 1997b) and referred to as ‘combined approach’, which is a combination of mosaic approach and analytical-statistical-dynamical approach, is proposed. Compared with those of other approaches, the main advantage of the combined approach is that it not only can represent both interpatch and intrapatch variability, but also cost less computational time when the land surface heterogeneity is considered. Because the independent variable of probability density function (PDF) is extended to the single valued function of basic meteorological characteristic quantities, which is much more universal, the analytical expressions of the characteristic quantities, (e.g., drag coefficient, snow coverage, leaf surface aerodynamical resistance) affected by roughness length are derived, when the roughness length (and/ or the zero plane displacement) heterogeneity has been mainly taken into account with the approach. On the basis of the rule which the PDF parameters should follow, we choose a function y of the roughness length z0 as the PDF independent variable, and set different values of the two parameters width ratio αn and height ratio γ of PDF (here a linear, symmetric PDF is applied) for sensitivity experiments, from which some conclusions can be drawn, e.g., relevant characteristic terms show different sensitivities to the heterogeneous characteristic (i.e., roughness length), which suggests that we should consider the heterogeneities of the more sensitive terms in our model instead of the heterogeneities of the rest, and which also implies that when the land surface scheme is coupled into the global or regional atmospheric model, sensitivity tests against the distribution of the heterogeneous characteristic are very necessary; when the parameterαn is close to zero, little heterogeneity is represented, andαn differs with cases, which have an upper limit of about 0.6; in the reasonable range ofαn, a peak-like distribution of roughness length can be depicted by a small value ofγ, etc..     

Local advection processes in the surface layer of the marginal ice zone

Micro-scale turbulent transport processes over the marginal ice zone have been studied by use of a two-dimensional numerical model. It has been found that internal boundary layers (IBLs) of horizontal mean velocity, temperature, and specific humidity reveal a near field and a far field. In the near field, the change in surface roughness dominates the height and growth rate of a velocity IBL. The change in surface heat flux governs the near field of a temperature and humidity IBL. In the far field, approximately x/¦L _*2 ¦ ~ 20, where L _*2 is the downstream Obukhov length, the downstream stratification more and more influences the growth rate of IBLs basically by modifying the eddy viscosity.Above more complex terrain consisting of an ensemble of ice strips and leads, a merging height h ^M develops, below which the horizontal variability of the surface modification is clearly observed; h ^M varies with the length scale L of surface modification approximately in proportion to h ^M /L ~ 1/20 – 1/10, as a rule of thumb. Above the merging height, an enveloping IBL exists, whose growth depends on the ice cover, i.e., on the integral of surface modification, but changes very little with L.Local advection of momentum, heat, and moisture clearly affects the local surface heat fluxes. Sensible and latent heat fluxes are found to show also a near and far field. However, if areally averaged surface fluxes are to be deduced from grid-averaged flow variables, then details of local advection can be neglected to a reasonably good approximation.     

Field And Wind-Tunnel Studies Of Aerodynamic Roughness Length

The aerodynamic roughness length (z₀) values of three Gobi desert surfaces were obtained by measurement of the boundary-layer wind profile in the field. To clarify the factors affecting the Gobi surface aerodynamic roughness length, a wind-tunnel experiment was conducted. The wind-tunnel simulation shows that z₀ values increase with increasingsize and coverage of roughness elements. Especially, the shape and height of roughnesselements are more important than other factors in affecting roughness length. The roughness length increases with decreasing values of the geometric parameter (the ratio of element horizontal surface area to height, ) of roughness elements. But at a higher free stream velocity, the height is more important than the shape in affecting roughness length.     

Momentum Absorption in Rough-Wall Boundary Layers with Sparse Roughness Elements in Random and Clustered Distributions

A wind-tunnel experiment has been used to investigate momentum absorption by rough surfaces with sparse random and clustered distributions of roughness elements. An unusual (though longstanding) method was used to measure the boundary-layer depth δ and friction velocity u _* and thence to infer the functional relationship z ₀/h = f(λ) between the normalised roughness length z ₀/ h and the roughness density λ (where z ₀ is the roughness length and h the mean height of the roughness elements). The method for finding u _* is based on fitting the velocity defect in the outer layer to a functional form for the dimensionless velocity-defect profile in a canonical zero-pressure-gradient boundary layer. For the conditions investigated here, involving boundary layers over sparse roughness with strong local heterogeneity, this velocity-defect-law method is found to be more robust than several alternative methods for finding u _* (uw covariance, momentum integral and slope of the logarithmic velocity profile).The experimental results show that, (1) there is general agreement in the relationship z ₀/h = f(λ) between the present experiment with random arrays and other wind-tunnel experiments with regular arrays; (2) the main effect of clustering is to increase the scatter in the z ₀/h = f(λ) relationship, through increased local horizontal heterogeneity; (3) this scatter obscures any trend in the z ₀/h = f(λ) relationship in response to clustering; and (4) the agreement between the body of wind-tunnel data (taken as a whole) and field data is good, though with scatter for which it is likely that a major contribution stems from local horizontal heterogeneity in the field.     

Aerodynamic Scaling for Estimating the Mean Height of Dense Canopies

We used an aerodynamic method to objectively determine a representative canopy height, using standard meteorological measurements. The canopy height may change if the tree height is used to represent the actual canopy, but little work to date has focused on creating a standard for determining the representative canopy height. Here we propose the ‘aerodynamic canopy height’ h _a as the most effective means of resolving the representative canopy height for all forests. We determined h _a by simple linear regression between zero-plane displacement d and roughness length z ₀, without the need for stand inventory data. The applicability of h _a was confirmed in five different forests, including a forest with a complex canopy structure. Comparison with stand inventory data showed that h _a was almost equivalent to the representative height of trees composing the crown surface if the forest had a simple structure, or to the representative height of taller trees composing the upper canopy in forests with a complex canopy structure. The linear relationship between d and z ₀ was explained by assuming that the logarithmic wind profile above the canopy and the exponential wind profile within the canopy were continuous and smooth at canopy height. This was supported by observations, which showed that h _a was essentially the same as the height defined by the inflection point of the vertical profile of wind speed. The applicability of h _a was also verified using data from several previous studies.     

Shear stress partitioning in large patches of roughness in the atmospheric inertial sublayer

Drag partition measurements were made in the atmospheric inertial sublayer for six roughness configurations made up of solid elements in staggered arrays of different roughness densities. The roughness was in the form of a patch within a large open area and in the shape of an equilateral triangle with 60 m long sides. Measurements were obtained of the total shear stress (τ) acting on the surfaces, the surface shear stress on the ground between the elements (τ_S) and the drag force on the elements for each roughness array. The measurements indicated that τ_S quickly reduced near the leading edge of the roughness compared with τ, and a τ_S minimum occurs at a normalized distance (x/h, where h is element height) of (downwind of the roughness leading edge is negative), then recovers to a relatively stable value. The location of the minimum appears to scale with element height and not roughness density. The force on the elements decreases exponentially with normalized downwind distance and this rate of change scales with the roughness density, with the rate of change increasing as roughness density increases. Average τ_S : τ values for the six roughness surfaces scale predictably as a function of roughness density and in accordance with a shear stress partitioning model. The shear stress partitioning model performed very well in predicting the amount of surface shear stress, given knowledge of the stated input parameters for these patches of roughness. As the shear stress partitioning relationship within the roughness appears to come into equilibrium faster for smaller roughness element sizes it would also appear the shear stress partitioning model can be applied with confidence for smaller patches of smaller roughness elements than those used in this experiment.     

Heat flux parameterization for sparse and dense grasslands with the analytical land-atmosphere radiometer model (alarm)

The dependence of radiometric surface temperature (_s) on view angle and the unclear definition of the aerodynamic temperature, which is the temperature that gives the correct sensible heat flux estimate at a specified roughness length, bring about a challenge in estimating sensible heat flux from _s. An analytical-land-atmosphere-radiometer model (ALARM) has been developed to convert _s taken at any zenith view angle to a clearly defined equivalent isothermal surface temperature, _i, at a defined scalar roughness length. ALARM is an analytical model based on K-theory that links the foliage temperature profile to the radiometric surface temperature and the temperature felt by the turbulent lower atmosphere. ALARM has previously been applied with slightly different values of its parameters to several grassland sites of varying canopy density. Our objective in this study was to apply ALARM to these and to one additional dataset with a single parameterization. When compared to the reference (measured) values of sensible heat flux H, ALARM estimates of H had root mean square errors of about 35 W m^-2. These results were comparable to those from two other simple canopy models also tested with these datasets.     

Estimates of Roughness Parameters for Arrays of Obstacles

Some methods are evaluated and extended to estimate roughness length and zero plane displacement height for atmospheric flow over arrays of obstacles, typically buildings. It appears that the method proposed by Bottema, with an extension to account for low density obstacle arrays, performs best. Procedures are proposed to represent irregular obstacle arrangements by a representative regular array to which Bottema's method can be applied. It is shown that this can be done without loss of accuracy, in general, roughness length can be predicted within a factor of two in more than 74% of the cases (95% reliability estimate). The methods proposed by Lettau and Raupach have been included in the evaluation. Lettau's model, which only requires input on the frontal area density, predicts roughness length unbiassed for frontal area densities up to 0.3, but predictions will be within a factor of two in more than 59% of the cases only (95% reliability estimate).     

Aerodynamic Parameters of Regular Arrays of Rectangular Blocks with Various Geometries

The aerodynamic effects of various configurations of an urban array were investigated in a wind-tunnel experiment. Three aerodynamic parameters characterising arrays—the drag coefficient (C _d ), roughness length (z _o) and displacement height (d)—are used for analysis. C _d is based on the direct measurement of the total surface shear using a floating element, and the other two parameters are estimated by logarithmic fitting of the measured wind profile and predetermined total drag force. The configurations of 63 arrays used for measurement were designed to estimate the effects of layout, wind direction and the height variability of the blocks on these parameters for various roughness packing densities. The results are summarised as follows: (1) The estimated C _d and z _o of the staggered arrays peak against the plan area index (λ _p ) and frontal area index (λ _f ), in contrast with values for the square arrays, which are less sensitive to λ _p and λ _f . In addition, the square arrays with a wind direction of 45° have a considerably larger C _d , and the wind direction increases z _o/H by up to a factor of 2. (2) The effect of the non-uniformity of roughness height on z _o is more remarkable when λ _f exceeds 20%, and the discrepancy in z _o is particularly remarkable and exceeds 200%. (3) The effect of the layout of tall blocks on C _d is stronger than that of short blocks. These results indicate that the effects of both wind direction and the non-uniformity of the heights of buildings on urban aerodynamic parameters vary greatly with λ _p and λ _f ; hence, these effects should be taken into account by considering the roughness packing density.