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).     

Methodological Considerations Regarding the Measurement of Turbulent Fluxes in the Urban Roughness Sublayer: The Role of Scintillometery

We address some of the methodological challenges associated with the measurement of turbulence and use of scintillometers in the urban roughness sublayer (RSL). Two small-aperture scintillometers were located near the roof interface in a densely urbanized part of Basel, Switzerland, as part of the Basel Urban Boundary-Layer Experiment (BUBBLE) in the summer of 2002. Eddy correlation instruments were co-located near the mid-point of each scintillometer path for data verification purposes. The study presents the first values of the inner length scale of turbulence (l ₀) and the refractive index structure parameter of air for a city and demonstrates the influence of mechanical driven turbulence on dissipation. Comparison of dissipation values determined from the two approaches show large scatter that is possibly due to the spatial inhomogeneity of the turbulence statistics within the RSL. Velocity and temperature spectra display a −2/3 slope in the inertial subrange, although the spectral ratio is less than the theoretical prediction of 4/3 expected for isotropy. Conventional Monin–Obukhov equations used to calculate fluxes from the scintillometer were replaced with urban forms of the equations. The results suggest that the scintillometer may be an appropriate tool for the measurement of sensible heat flux (Q _H ) above the rooftops given a suitable determination of the effective measurement height.     

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.     

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.     

Area-Averaged Sensible Heat Flux and a New Method to Determine Zero-Plane Displacement Length over an Urban Surface using Scintillometry

Field observations of area-averagedturbulence characteristics were conducted in a densely built-up residential neighbourhood in Tokyo, Japan. In addition to eddy-correlation (EC) sensors a scintillometer was used for the first time in a city. Significant results include: (1) Scintillometer-derived sensible heat fluxes, Q_H, obtained at a height 3.5 times the building height agree well with those using the EC technique; (2) source areas for the scintillometer fluxes are larger than for the EC sensors, so that at low heights over inhomogeneous terrain scintillometry offers advantages; (3) new similarity relationships for dissipation rates are proposed for urban areas; (4) a new technique that uses simultaneous scintillation measurements at two heights to directly estimate area-averaged zero-plane displacement height, z_d, is proposed. z_d estimated in this way depends slightly on atmospheric stability (lower z_d under more unstable conditions).     

Comparison Of Upwind Land Use And Roughness Length Measured In The Urban Boundary Layer

The results of wind measurements taken in two studies at an urban sitein Birmingham (UK) are presented. Displacement heights and roughnesslengths are calculated that compare well with previous measurements inthe urban canopy. Dimensionless windshear and vertical-velocity standarddeviation are also shown to obey Monin–Obukhov similarity theory. Themean roughness lengths measured at the highest measuring position (45 m)increase with the fraction of upwind urban cover calculated using asource-area model.     

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.     

Influence of Obstacles on the Aerodynamic Roughness of the Netherlands

The aim of this study was to analyse the influence of large- and small-scale obstacles (orography, tree lines, and dikes) on the effective aerodynamic roughness of the Netherlands, a relatively flat, small-scale landscape. The roughness averaging approach was based on drag coefficients. The effective roughness was locally dominated by small-scale obstacles such as tree lines and dikes. Even at a regional scale (40,000 km2), the small-scale obstacle drag was of the same order of magnitude as the shear stress due to landuse. The neglect of those obstacles on a regional scale would result in approximately 10% overestimated averaged windspeed at 10~m above the surface. It was concluded that small-scale obstacles need to be taken into account to calculate the aerodynamic roughness of flat landscapes. Orography was of minor importance in this lowland country.     

A Roughness Sublayer Wind Profile Above A Non-Uniform Surface

In atmospheric models for different scales the underlying surface consists of patches of bare soil and plant communities with different morphological parameters. Experimental evidence indicates that there is a significant departure of the wind profile above a vegetative surface from that predicted by the logarithmic relationship, which gives values that are greater than those observed. This situation can seriously disturb the physical picture concerning the transfers of momentum, heat and water vapour from the surface into the atmosphere.The intention of this paper is to generalise the calculation of exchange of momentum between the atmosphere and a non-homogenous vegetative surface, and to derive a general equation for the wind speed profile in a roughness sublayer under neutral conditions. Furthermore, these results are extended to non-neutral cases. The suggested expression for the wind profile is compared with some earlier approaches and the observations obtained above a broad range of plant communities.     

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.     

盘锦芦苇湿地空气动力学参数动态特征及其影响因素分析

应用盘锦芦苇生态系统观测场的涡度相关通量观测系统和小气候梯度系统的观测资料,估算了芦苇下垫面的空气动力学参数,并分析了其影响因素。结果表明:芦苇湿地零平面位移d和粗糙度Z0的季节变化总体呈现先增加后减少的单峰曲线变化规律。Z0和d均在9月达到最大,分别为0.24 m和1.85 m;12月达到最小,分别为0.03 m和0.02 m。芦苇湿地的Z0和d的季节变化主要受叶面积指数(LAI)、株高(h)及风速的影响。d随h的增加而单调增加,Z0随着h的变化遵循二次曲线,d/h、Z0/h与LAI之间的复相关系数分别为0.99和0.78。研究亦表明,大气稳定度和风向也对d和Z0有一定的影响。     

Modelling The Mean Velocity Profile In The Urban Canopy Layer

A simple model originally derived for meanwind speed profiles in vegetative canopy flows ismodified for application to arrays ofthree-dimensional surface obstacles (cubes), whichcould be representative of a simple urban-typesurface. It is shown that for cube arrays that arenot too densely packed, the predicted exponentialvelocity profile provides an adequate fit to thespatially averaged velocity profile (u(z))within the obstacle canopy. Application of the model to a set of wind-tunnel dataallows for the evaluation of an empirical fittingparameter called the attenuation coefficient. This isrelated to the turbulence length scale, which can befound by manipulating the results of thegradient-diffusion model used to derive the velocityprofile. The results show a reduction of theturbulence length scale with increasing obstaclepacking density. By assuming a linear transition fromthis length scale at the top of the canopy to theclassical Prandtl length scale in the overlyinginertial sublayer, an acceptable model is obtained forthe full velocity profile within simple obstaclearrays, from the ground up to the overlyingsemi-logarithmic region.     

基于MODIS资料的中国东部时间序列空气动力学粗糙度和零平面位移高度估算

空气动力学粗糙度和零平面位移高度是很多气候模型和陆面模式中的重要参数,采用气象学方法推导这两个参数对于大范围长时间序列的计算需要大量长期的野外观测,而遥感方法可以快速经济的提供大范围数据,在本研究中采用形态学模型,以MOD IS产品数据作为数据源,计算植被冠层面积指数,估算了中国东部2001—2003年归一化到植被高度的1 km空间分辨率时间序列空气动力学粗糙度和零平面位移高度,结果表明空气动力学粗糙度和零平面位移高度均存在季节变化特征,并对所采用的模型和参数估算的结果进行了讨论。     

Spatial Variability of Both Turbulent Fluxes and Temperature Profiles in an Urban Roughness Layer

The spatial variability of both turbulent flow statistics in the roughness sublayer (RSL) and temperature profiles within and above the canopy layer (CL) were investigated experimentally in a densely built-up residential area in Tokyo, Japan. Using five towers with measuring devices, each tower isolated from the others by at least 200 m, we collected high-frequency measurements of velocity and temperature at a height z=1.8 z _H, where z _H, the mean building height in the area, is 7.3 m. Also, temperature profiles were measured from z=0.4 to 1.8 z _H. The ‘areal mean’ geometric parameters that were obtained for the areas within 200 m of each tower were fairly homogeneous among the tower sites. The main results are as follows: (1) The spatial variability of all RSL turbulent statistics, except the sensible heat flux, was comparable to that reported in a pine forest. Also, the variability decreased with increasing friction velocity. (2) The spatial variability of the RSL sensible heat flux was larger than that reported in a pine forest. Also, the variability depended on the time of the day and became larger in the morning. The difference among the sites was well related to the areal fraction of vegetation. (3) The spatial variability of the CL temperature profile depended on the time of the day and became larger in the morning. Nevertheless, the spatial standard deviation of CL temperature was always below 0.7 K. (4) It is suggested that the “warming-up” process in the morning when heat storage is dominant increases the spatial variation of RSL sensible heat flux and CL temperature according to the local properties around each tower and the variation decreases once there is further convective mixing in the midday     

基于MODIS资料的中国东部时间序列空气动力学粗糙度和零平面位移高度估算

空气动力学粗糙度和零平面位移高度是很多气候模型和陆面模式中的重要参数,采用气象学方法推导这两个参数对于大范围长时间序列的计算需要大量长期的野外观测,而遥感方法可以快速经济的提供大范围数据,在本研究中采用形态学模型,以MODIS产品数据作为数据源,计算植被冠层面积指数,估算了中国东部2001-2003年归一化到植被高度的...     

大规模风电场建成后对风能资源影响的研究

考虑了千万千瓦级风电基地建成后风电机群对近地面层风速的影响,采用Frandsen研究了大规模风电场内部风速损失时所使用的方法,在内边界层已经充分发展成新边界层的区域内,对轮毂高度65m处风速Uh进行了计算。结果表明,风电场建成后研究区内,风速Uh与未建场时的65m风速U0相比变小,存在风速损失,该风速损失随着U0的增大而减小,与风电机的推力系数CT性质有关;大规模风电场建成后,Uh在3～20m.s-1范围内的平均风功率密度与未建场时U0在此范围内的平均风功率密度相比损失约为58.45%,这与建场地区建场前65m处风速值大小以及各风速值出现的概率有关。     

Influence of a Two-scale Surface Roughness on a Neutral Turbulent Boundary Layer

Flow in the urban boundary layer is strongly influenced by the surface roughness, which is composed principally of isolated buildings or groups of buildings. Previous research has shown that the flow regime depends on the characteristic height of these obstacles (H), and the spacing between them (W). In reality, the urban boundary layer contains roughness elements with a wide range of length scales; in many practical situations these can be classified into large-scale roughness—buildings, or groups of buildings—and small-scale roughness, such as street furniture and elements on the façades and roofs. It is important to understand how the small-scale roughness might modify mass and momentum transfer in the urban boundary layer, but relatively little information is available concerning the potential interaction between large- and small-scale roughness elements in the different flow regimes. This problem has been studied using wind-tunnel experiments, by measuring vertical velocity profiles over a two-dimensional obstacle array, adding small-scale roughness elements to the top of larger parallel square bars. The experiments were performed for different cavity aspect ratios: the results show that the small-scale roughness increases the turbulence intensities and the momentum transfer when the large-scale obstacles are closely packed (H/W > 1) but it has very little effect for more widely-spaced obstacles (H/W < 1).     

城市下垫面空气动力学参数确定方法研究综述

地表粗糙度和零平面位移是2个重要的空气动力学参数,对于城市下垫面,粗糙元的形态非常复杂,空间分布上存在较大的非均匀性,如何确定城市下垫面地表粗糙度和零平面位移目前尚无最佳方案。概述了这2个参数适用于城市非均匀下垫面的参数化方法,包括风温观测方法和粗糙元形态学方法。探讨了各自在应用过程中存在的问题及进一步的研究工作,并且提出,通过数学模型的改进和卫星遥感技术的引入,粗糙元形态学方法会有更好的发展前景。     

Equations for the Drag Force and Aerodynamic Roughness Length of Urban Areas with Random Building Heights

We use a conceptual model to investigate how randomly varying building heights within a city affect the atmospheric drag forces and the aerodynamic roughness length of the city. The model is based on the assumptions regarding wake spreading and mutual sheltering effects proposed by Raupach (Boundary-Layer Meteorol 60:375?C395, 1992). It is applied both to canopies having uniform building heights and to those having the same building density and mean height, but with variability about the mean. For each simulated urban area, a correction is determined, due to height variability, to the shear stress predicted for the uniform building height case. It is found that u _*/u _*R , where u _* is the friction velocity and u _*R is the friction velocity from the uniform building height case, is expressed well as an algebraic function of ?? and ?? _h /h _m , where ?? is the frontal area index, ?? _h is the standard deviation of the building height, and h _m is the mean building height. The simulations also resulted in a simple algebraic relation for z ₀/z _0R as a function of ?? and ?? _h /h _m , where z ₀ is the aerodynamic roughness length and z _0R is z ₀ found from the original Raupach formulation for a uniform canopy. Model results are in keeping with those of several previous studies.     

城市下垫面空气动力学参数确定方法综述

地表粗糙度和零平面位移是2个重要的空气动力学参数,对于城市下垫面,粗糙元的形态非常复杂,空间分布上存在较大的非均匀性,如何确定城市下垫面地表粗糙度和零平面位移目前尚无最佳方案。概述了这2个参数适用于城市非均匀下垫面的参数化方法,包括风温观测方法和粗糙元形态学方法。探讨了各自在应用过程中存在的问题及进一步的研究工作,并且提出,通过数学模型的改进和卫星遥感技术的引入,粗糙元形态学方法会有更好的发展前景：