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.     

Near Wall Flow over Urban-like Roughness

In this study, comprehensive measurements over a number of urban-type surfaces with the same area density of 25% have been performed in a wind tunnel. The experiments were conducted at a free stream velocity of 10 m s^-1 and the main instrumentation was 120 ° x-wire anemometry, but measurement accuracy was checked using laser Doppler anemometry.The results haveconfirmed the strong three-dimensionalityof the turbulent flow inthe roughness sublayer and the depths of the inertial sublayer (log-law region) and roughness sublayer for each surface have been determined. Spatial averaging has been used to remove the variability of the flow in the roughness sublayer due to individual obstacles and it is shown that the spatially averaged mean velocity in the inertial sublayer and roughness sublayer can,together, be described by a single log-law with a mean zero-plane displacement and roughness length for the surface, provided that the proper surface stress is known. The spatially averaged shear stresses in the inertial sublayer and roughness sublayer are compared with the surface stress deduced from form drag measurements on the roughness elements themselves.The dispersive stress arising from the spatial inhomogeneity in the mean flow profiles was deduced from the data and is shown to be negligible compared with the usual Reynolds stresses in the roughness sublayer. Comparisons have been made between a homogeneous (regular element array) surface and one consisting of random height elements of the same total volume. Although the upper limits of the inertial sublayer for both surfaces were almost identical at equivalent fetch, the roughness sublayer was much thicker for the random surface than for the uniform surface, the friction velocity and the roughness length were significantly larger and the `roughness efficiency' was greater. It is argued that the inertial sublayer may not exist at all in some of the more extreme rough urban areas. These results will provide fundamental information for modelling urban air quality and forecasting urban wind climates.     

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

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

The effect of uncertainty in aerodynamic resistance on evaporation estimates from the combination model

Evaporation estimates from a soybean crop calculated from the combination model are insensitive to aerodynamic resistance. The insensitivity arises from a strong link between evaporation and the vapour pressure deficit of the air and bulk stomatal resistance. The sensitivity of aerodynamic resistance to errors in surface roughness and zero-plane displacement is considered. The resistance is found to be more sensitive to errors in surface roughness than to errors in zero-plane displacement. However, large errors in these have little effect on calculated evaporation. Both surface roughness and zero-plane displacement are empirically related to crop height and leaf area index. Errors incurred by ignoring bluff-body effects and atmospheric stability are small in estimating both resistance and evaporation. Evaporation can be calculated adequately from empirical estimates of surface roughness and zero-plane displacement and single-level measurements of windspeed.     

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.     

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

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

Estimation of Roughness Parameters Within Sparse Urban-Like Obstacle Arrays

We conduct wind-tunnel experiments on three different uniform roughness arrays composed of sparsely distributed rectangular cylinders for the estimation of surface parameters. Roughness parameters such as the roughness length z ₀ and zero-plane displacement d are extracted using a best-fit approximation of the measured wind velocity. We also perform a large-eddy simulation (LES) to confirm that four sampling points are sufficient to surrogate a space average above the canopy layer of the sparse roughness arrays. We propose a new morphological model from a systematic analysis of experimental data on the arrays. The friction velocity predicted by the proposed model agrees well with the peak value of the measured Reynolds shear stress ${(-\left<\overline{u'w'}\right>)^{0.5}}${(-\left<\overline{u'w'}\right>)^{0.5}}. The proposed model is further validated in an additional wind-tunnel experiment conducted on a scaled configuration of a real urban area exposed to four wind directions. The proposed model is found to perform very well particularly in the estimation of the friction velocity, readily leading to a better estimation of turbulence, which is essential for an accurate prediction of pollutant dispersion.     

Study of Aerodynamic Parameters on Different Underling Surfaces

Aerodynamic parameters including the zero-plane displacement (d), roughness length (20), and friction velocity (u*) on the different underlying surfaces of heavy-grazing site, medium-grazing site, light-grazing site, no-grazing site, dune, inter-dune, grassland, rice paddy site, wheat site, soybean site, and maize site have been computed based on the Monin-Obukhov similarity theory by utilizing the micrometeorologically observed data of dune and vegetation in the semi-arid area at Naiman, Inner Mongolia of China, conducted jointly by the Institute of Desert Research, Chinese Academy of Sciences and the National Institute of Agro-Environmental Sciences of Japan in 1990-1994. And their relationships between wind speed and Richardson number are analyzed. The aerodynamic characteristics of different man-made disturbed grassland ecosystems are also compared. Result shows that the vegetation coverage and the above-ground biomass decrease with the increase in man-made stress of the grassland. The roughness length for different underlying surfaces is closely related to vegetation height, above-ground biomass, and ground surface undulation, and Richardson number Ri is also its influencing factor. The friction velocity varies largely on different underlying surfaces, and it is positively proportional to wind speed and roughness length. The aerodynamic parameters of various times on the same underlying surface are different, too. Above results indicate that grassland and vegetation are of significance in preventing desertification, especially in the arid and semi-arid land ecosystems. And the results of this paper are also important for constructing the land surface physical process as well as regional climate model.     

Study on Effects of Building Morphology on Urban Boundary Layer Using an Urban Canopy Model

An urban canopy model is incorporated into the Nanjing University Regional Boundary Layer Model. Temperature simulated by the urban canopy model is in better agreement with the observation, especially in the night time, than that simulated by the traditional slab model. The coupled model is used to study the effects of building morphology on urban boundary layer and meteorological environment by changing urban area, building height, and building density.It is found that when the urban area is expanded, the urban boundary layer heat flux, thermal turbulence, and the turbulent momentum flux and kinetic energy all increase or enhance, causing the surface air temperature to rise up. The stability of urban atmospheric stratification is affected to different extent at different times of the day.When the building height goes up, the aerodynamic roughness height, zero plane displacement height of urban area, and ratio of building height to street width all increase. Therefore, the increase in building height results in the decrease of the surface heat flux, urban surface temperature, mean wind speed, and turbulent kinetic energy in daytime. While at night, as more heat storage is released by higher buildings, thermal turbulence is more active and surface heat flux increases, leading to a higher urban temperature.As the building density increases, the aerodynamic roughness height of urban area decreases, and the effect of urban canopy on radiation strengthens. The increase of building density results in the decrease in urban surface heat flux, momentum flux, and air temperature, the increase in mean wind speed, and the weakening of turbulence in the daytime. While at night, the urban temperature increases due to the release of more heat storage.     

A Scale-Adaptive Approach for Spatially-Varying Urban Morphology Characterization in Boundary Layer Parametrization Using Multi-Resolution Analysis

Urban morphology characterization is crucial for the parametrization of boundary-layer development over urban areas. One complexity in such a characterization is the three-dimensional variation of the urban canopies and textures, which are customarily reduced to and represented by one-dimensional varying parametrization such as the aerodynamic roughness length $z_{0}$ and zero-plane displacement $d$ . The scope of the paper is to provide novel means for a scale-adaptive spatially-varying parametrization of the boundary layer by addressing this 3-D variation. Specifically, the 3-D variation of urban geometries often poses questions in the multi-scale modelling of air pollution dispersion and other climate or weather-related modelling applications that have not been addressed yet, such as: (a) how we represent urban attributes (parameters) appropriately for the multi-scale nature and multi-resolution basis of weather numerical models, (b) how we quantify the uniqueness of an urban database in the context of modelling urban effects in large-scale weather numerical models, and (c) how we derive the impact and influence of a particular building in pre-specified sub-domain areas of the urban database. We illustrate how multi-resolution analysis (MRA) addresses and answers the afore-mentioned questions by taking as an example the Central Business District of Oklahoma City. The selection of MRA is motivated by its capacity for multi-scale sampling; in the MRA the “urban” signal depicting a city is decomposed into an approximation, a representation at a higher scale, and a detail, the part removed at lower scales to yield the approximation. Different levels of approximations were deduced for the building height $\bar{{H}}$ and planar packing density $\lambda _\mathrm{p}$ . A spatially-varying characterization with a scale-adaptive capacity is obtained for the boundary-layer parameters (aerodynamic roughness length $z_{0}$ and zero-plane displacement $d$ ) using the MRA-deduced results for the building height and the planar packing density with a morphometric model; an attribute that is shown to be of great advantage to multi-scale and multi-resolution numerical weather prediction models.     

应用城市冠层模式研究建筑物形态对城市边界层的影响

文中将城市冠层模式耦合到南京大学城市尺度边界层模式中,通过模拟对比发现,耦合模式对城市地区气温模拟结果更接近于观测值,尤其是对城市地区夜间气温模拟的改进.运用改进耦合模式通过多个敏感性试验的模拟,从城市面积扩张、建筑物高度增加、建筑物分布密度变化等角度研究城市建筑物三维几何形态变化对城市边界层及城市气象环境的影响,试验结果表明:(1)城市面积扩张使得城市下垫面的热通量增大,热力湍流活动增强,动量通量输送增强,城市湍能增大,湍流扩散系数变大,城市气温升高,且对不同时刻城市区域大气层结稳定度均有不同程度的影响.(2)建筑物高度增加增大了城市下垫面的粗糙度和零平面位移.同时也增大了城市街渠高宽比.城市建筑物越高,白天城市地区地表热通量越小,城市上空大气温度越低,平均风速减小,湍能减小;夜间由于高大建筑物释放储热比低矮建筑物要多,其热力湍流相对活跃,地表热通量增大,使得城市区域气温较高.(3)建筑物密度增大,会减小城市下垫面的粗糙度同时增强街渠对辐射的影响.建筑物密度增大在白天会减小地表热通量和动量通量,使城市气温降低,平均风速增大,城市湍流活动能力减弱;夜间城市释放较多储热使得气温较高.     

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.     

基于近地层垂直速度能谱尺度特征确定复杂下垫面零平面位移

本文利用位于南京市郊区的南京大学仙林校区SORPES观测站多层湍流观测数据分析了湍流谱特征,以白天不稳定条件下垂直速度能谱谱峰对应的长度尺度也就是离地高度为判据,探讨了运用该方法确定复杂下垫面零平面位移的可行性。统计分析表明,该方法确定的长度尺度呈现出较为一致的概率分布形状,概率最大的长度尺度对应于离地高度,在复杂下垫面情况下这个高度就是零平面位移高度到观测高度之间的距离,将观测点的离地高度减去这个距离就能得到零平面位移。本文同时运用不稳定条件下垂直速度方差在近地层中的相似关系来确定零平面位移,并与谱方法得到的结果进行对比。结果表明,谱方法和方差法得到的零平面位移非常接近。     

Impact of surface variations on the momentum flux above the urban canopy

Based on the momentum flux–wind profile relationship of the Monin–Obukhov Similarity (MOS) theory, the observational data from the urban boundary layer field campaign in Nanjing are used to calculate the friction velocity ( $ {u_*} $ ) at the top of the urban canopy and the calculated results are evaluated. The urban surface roughness parameters (the roughness length z ₀ and zero-plane displacement height z _d) are estimated with the Ba method (Bottema’s morphological method). Two different regimes are employed for the calculations. In the homogeneous approach, z ₀ and z _d are averagely derived from the surface elements in the whole study area; while in the heterogeneous approach, z ₀ and z _d are locally derived from the surface elements in the corresponding upwind fetches (or source areas). The calculated friction velocities are compared to the measurement data. The results show that the calculated friction velocities from the heterogeneous approach are in better agreement with the observed values than those from the homogeneous approach are. This study implies that the local roughness parameters can properly represent the dynamical heterogeneity of urban surface, and its application can significantly improve the performance of parameterizations based on the MOS theory in the urban roughness sublayer.     

Aerodynamic Parameters of Urban Building Arrays with Random Geometries

It is difficult to describe the flow characteristics within and above urban canopies using only geometrical parameters such as plan area index (λ _p ) and frontal area index (λ _f ) because urban surfaces comprise buildings with random layouts, shapes, and heights. Furthermore, two types of ‘randomness’ are associated with the geometry of building arrays: the randomness of element heights (vertical) and that of the rotation angles of each block (horizontal). In this study, wind-tunnel experiments were conducted on seven types of urban building arrays with various roughness packing densities to measure the bulk drag coefficient (C _d ) and mean wind profile; aerodynamic parameters such as roughness length (z _o ) and displacement height (d) were also estimated. The results are compared with previous results from regular arrays having neither ‘vertical’ nor ‘horizontal’ randomness. In vertical random arrays, the plot of C _d and z _o versus λ _f exhibited a monotonic increase, and z _o increased by a factor of almost two for λ _f = 48–70%. C _d was strongly influenced by the standard deviation of the height of blocks (σ) when λ _p ≥ 17%, whereas C _d was independent of σ when λ _p = 7%. In the case of horizontal random arrays, the plot of the estimated C _d against λ _f showed a peak. The effect of both vertical and horizontal randomness of the layout on aerodynamic parameters can be explained by the structure of the vortices around the blocks; the aspect ratio of the block is an appropriate index for the estimation of such features.     

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

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

The Characteristics and Parameterization of Aerodynamic Roughness Length over Heterogeneous Surfaces

Aerodynamic roughness length (z0m is a key factor in surface flux estimations with remote sensing algorithms and/or land surface models. This paper calculates z0m over several land surfaces, with 3 years of experimental data from Xiaotangshan. The results show that z0m is direction-dependent, mainly due to the heterogeneity of the size and spatial distribution of the roughness elements inside the source area along different wind directions. Furthermore, a heuristic parameterization of the aerodynamic roughness length for heterogeneous surfaces is proposed. Individual z0m over each surface component (patch) is calculated firstly with the characteristic parameters of the roughness elements (vegetation height, leaf area index, etc.), then z0m over the whole experimental field is aggregated, using the footprint weighting method.     

Experimental study of the transfer velocity for urban surfaces with a water evaporation method

A major problem in urban climate modelling is determining how the heat fluxes from various canyon surfaces are affected by canyon flow. To address this problem, we developed a water evaporation method involving filter paper to study the distribution of the convective transfer velocity in urban street canyons. In this method, filter paper is pasted onto a building model and the evaporation rate from the paper is measured with an electric balance. The method was tested on 2D (two-dimensional) street canyon models and 3D model arrangements. Moreover, in this technique, it is easy to restrict the flux within an arbitrary surface in question. That is, the evaporation distribution on a surface can be studied by using several small pieces of filter paper. In the 2D case, the wall transfer velocity was strongly dependent on the canyon aspect ratio for perpendicular wind directions and it varied widely with height within both windward and leeward wall surfaces. For 3D cubic arrays, the relation to canyon aspect ratio was largely different from that of the 2D canyon. And, as a case study, the variation of wind direction was investigated for a city-like setting. The area-averaged transfer velocity was insensitive to wind direction but its local deviation was significant. Finally, we measured the transfer velocity for a clustered block array surrounded by relatively wide streets. The effect of spatial heterogeneity on the transfer velocity was significant. Moreover, for a fixed total building volume, the transfer velocity was considerably larger when the building height varied than when it was uniform. Therefore, the water evaporation method with filter paper is expected to be useful for studying the transfer velocity and ventilation rates in urban areas with various canyon shapes.