典型特大城市下垫面空气动力学特征参数的计算与分析研究

基于2012年12月—2013年11月上海徐家汇气象铁塔风和湍流观测资料及铁塔周边500m半径范围(平均建筑高度为39.5 m)下垫面信息,分别利用温度方差方法、风速标准差方法、涡度相关法和形态学方法计算典型特大城市下垫面零平面位移和动力学粗糙度。结果表明,温度脉动方差法计算的零平面位移与建筑物高度分布存在差异,但在58~208°方向与建筑顶面积系数呈现正相关(相关系数0.73),计算公式的经验参数、稳定度阈值和风速对计算结果影响较为显著;涡度相关方法和风速标准差方法计算的动力学粗糙度随风向分布趋势相似(相关系数0.94),并与零平面位移呈负相关关系。形态学方法与动力学方法之间数值和随风向的变化趋势均存在差异,风温观测仪器源区和形态学方法分析区域具有显著影响。     

城市下垫面空气动力学参数的估算

为了定量描述北京城市下垫面的空气动力学特征,为模式提供准确的下垫面参数.利用气象塔大气湍流观测资料,结合Martano(2000)由单层超声风速、温度资料估算非均匀下垫面空气动力学参数的方法,计算了中国科学院大气物理研究所气象塔附近的下垫面空气动力粗糙度z0和零平面位移d,即z0为1.75 m和d为40.12 m.该结果与前人研究成果的比较结果表明最近的8年内,该处的零平面位移和空气动力粗糙度明显增大,这与该塔周围城市建设状况吻合.     

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

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

复杂地形条件下零平面位移和空气动力学粗糙度的计算——以珠海南亚热带常绿阔叶林地区为例

空气动力学粗糙度和零平面位移是影响陆面过程的重要参数。以往对空气动力学粗糙度和零平面位移的计算往往针对均一程度较高的下垫面。以珠海南亚热带常绿阔叶林地区为例,探索了复杂地形下垫面空气动力学粗糙度和零平面位移的计算。用拟合法和粗糙元法计算了该地区下垫面的空气动力学粗糙度和零平面位移,并分析了两种方法的误差来源。计算结果表明,拟合法在这种复杂地形下的计算准确度并不理想,由于这种方法的代表范围较小,易受到周边地形影响而使测量结果不符合对数廓线的形式,从而得不到准确的结果。与拟合法相比,粗糙元法可以给出所计算下垫面的具体范围。粗糙元法计算结果的代表范围比拟合法大得多,因此,粗糙元法比拟合法更加能够反映的是该区域下垫面空气动力学粗糙度和零平面位移的总体情况。另外,粗糙元法可以计算那些因为风速资料不足无法用风速廓线直接加以计算的区域的空气动力学粗糙度和零平面位移。根据拟合法的计算结果,冬季315~45°范围内,z_0=2.96 m,d=22.48 m;45~135°范围内,z_0=1.26 m,d=8.81 m。春夏季315~45°范围内,z_0=3.90 m,d=27.00 m;45~135°范围内,z_0=1.51 m,d=14.83 m。根据粗糙元法的计算结果,东南西北四个方向1 km×1 km范围内z_0分别为0.94 m、1.28 m、1.30 m、2.08 m;d分别为13.87 m、18.79 m、19.12 m、30.61 m。     

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

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

非单一水平均匀下垫面空气动力学参数的确定

文章介绍了一种利用单一高度风速、温度湍流资料确定空气动力学参数的方法。该方法无须进行风速廓线的测量，可以应用于非单一水平均匀下垫面和非中性层结，避免了主观性。用该方法计算了北京城市北部边缘325 m气象塔附近的零值位移d和地表粗糙 z0。结果表明:该处下垫面零值位移d和地表粗糙度z0与风向有很强的依赖关系，与气象塔周围的城市建设相对应。     

鄱阳湖地区零平面位移及动力粗糙度的计算

利用2013年7月1日至2014年6月30日鄱阳湖东岸70 m铁塔涡动相关观测资料,应用Martano方法和TVM（Temperature Variance Method）方法分别计算了该地地表零平面位移d和粗糙度z₀,通过代回Monin-Obukhov相似性理论的风廓线关系计算摩擦速度,以验证与实测摩擦速度的一致性。速度和温度标准差的归一化拟合线分别与Panosky等和Tillman给出的曲线趋势一致,表明该站观测数据总体满足近地层相似性。Martano方法计算结果随季节变化较大,春夏季的粗糙度是秋冬季的6.3倍;陆面方向零平面位移和粗糙度分别为来自湖面的2倍和10倍;Martano方法比TVM方法对季节和方向的敏感性更强。Martano方法计算得到零平面位移和粗糙度对摩擦速度造成了约9.9%的高估;而TVM方法对摩擦速度造成了约32.8%的高估;Martano方法计算的摩擦速度和观测值的一致性更好。     

深圳356米气象塔观测数据的质量控制方法和空气动力学参数研究

气象高塔数据资料弥足珍贵, 对其进行质量控制将为后续科学研究和业务工作的开展提供便利; 此外, 利用塔基观测资料计算空气动力学参数有助于校正模式空气动力学参数理论值。对2017-2018年深圳356 m气象梯度观测塔共13层的每10 s风速、风向、相对湿度、温度探测资料进行数据质量控制, 基于莫宁-奥布霍夫相似理论和数据质量控制后的气象梯度观测塔近地层(10 m、20 m、40 m、50 m和80 m) 1分钟平均的风温资料, 利用最小二乘法拟合迭代计算了近中性条件下深圳气象梯度观测塔下垫面空气动力学粗糙度(z₀)和零平面位移(d)。结果表明:深圳气象梯度观测塔的气象探测资料数据质量很高, 连续两年平均数据缺失率为1.28%, 数据错误率为0.01%。近中性边界层条件下, 深圳气象梯度观测塔下垫面空气动力学粗糙度均值为0.35 m, 零平面位移均值为5.33 m, 结果合理可信。研究表明空气动力学参数受下垫面非均匀性、植株柔软性、气流来向、风速等的共同影响。      

高大建筑物影响城市粗糙副层流场特征的数值研究

用数值模拟方法研究了高大建筑物对城市粗糙副层气流场特征的影响。数值模式采用基于雷诺平均纳维—斯托克斯方程组的应用计算流体力学FLUENT软件,次网格湍流参数化选用k-ε闭合方案。建筑物用立方体表示,并规则排列于模拟区域内。通过改变高大建筑物的数量与位置,对建筑物阵列内及其上空的气流特征进行了多个算例的数值模拟。依据模拟结果计算获得建筑物区域的面积平均风速廓线,结果表明各算例的粗糙副层风速廓线各不相同。运用动力学方法由风速廓线计算出各算例的零平面位移高度和粗糙度,并与几种计算零平面位移高度和粗糙度的形态学方法进行了比较检验。结果表明两种形态学方法(Ba、Ma)计算所得的零平面位移高度与动力学方法计算结果很接近,但对于粗糙度而言,几种形态学方法的计算结果都明显偏高。     

黄河上游玛曲地区空气动力学参数的确定及其在陆面过程模式中的应用

为了定量描述黄河上游玛曲地区草地下垫面的空气动力学特征,为模式提供准确的下垫面参数,利用"黄河源区气候与环境综合观测研究站"2006年9月的湍流观测资料,结合Martano由单层超声观测资料确定下垫面空气动力学参数的方法,计算了黄河上游玛曲地区草原下垫面空气动力学粗糙度z0和零平面位移d,即z0为0.035 m,d为0.143 m。同时,将z0和d应用于陆面过程模式CoLM中,检验其对模式模拟性能的影响,结果表明,改进陆面参数后的模式对感热通量和潜热通量的模拟均有明显改善。     

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.     

Vertical variation of roughness length at the Boulder Atmospheric Observatory

Observations near the BAO tower site have suggested a roughness length of about 0.01 m. Wind profiles from the tower yield values between 0.04 and 0.30 m. We now analyze fluctuations of vertical velocity at 10 m, and get intermediate values, differing with wind direction. Thus, inferred z ₀ values increase systematically with height.     

Regional roughness of the landes forest and surface shear stress under neutral conditions

Mean wind velocity profiles were measured by means of radio-windsondes over the Landes region in southwestern France, which consists primarily of pine forests with scattered villages and clearings with various crops. Analysis of neutral profiles indicated the existence of a logarithmic layer between approximately z – d ₀ = 67(±18)z ₀ and 128(+-32)z ₀ (z is the height above the ground, z ₀ the surface roughness and d ₀ the displacement height). The upper limit can also be given as z – d ₀ = 0.33 (±0.18)h, where h is the height of the bottom of the inversion. The profiles showed that the surface roughness of this terrain is around 1.2 m and the displacement height 6.0 m. Shear stresses derived from the profiles were in good agreement with those obtained just above the forest canopy at a nearby location with the eddy correlation method by a team from the Institute of Hydrology (Wallingford, England).     

Zero-plane displacement and roughness length for tall vegetation,derived from a simple mass conservation hypothesis

A method for the determination of the zero-plane displacement, d, and roughness length, z ₀, for tall vegetation is described. A new relationship between d and z ₀ is developed by imposing the condition of mass conservation on the logarithmic wind profile. Further, d and z ₀ can be evaluated directly if independent measurements of friction velocity are available in addition to wind profile measurements. The proposed method takes into account the existence of a transition layer immediately above the vegetation where the logarithmic wind profile law is not valid. Only one level of wind speed measurements is necessary within the inertial sub-layer.The method is applied to wind profile and eddy correlation measurements taken in and above an 18.5 m pine forest to yield d = 12.7 m and z ₀ = 1.28 m. The choice of height for the upper level of measurement and problems with measuring canopy flow are discussed.Work carried out while on leave at the Institute of Hydrology.     

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.     

An analytical method for the determination of the roughness parameters over complex regions

The values of roughness length for momentum z ₀and zero-plane displacement d ₀over a hilly rough complex region with vegetation were evaluated without any assumption concerning z ₀and d ₀.It was found that for widely scattered profile data, the method of least squares will not give a reasonable result in determining the roughness parameters. For this purpose, the method of maximum correlation was introduced instead. This method gave a fair result for captive balloon observations conducted in hilly terrain mainly covered with forest in the northwestern part of the Kanto Plain, Japan.     

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.     

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.     

Roughness length for heat over an urban canopy

The roughness length for heat z_T was evaluated over an urban canopy, using the measured sensible heat flux and radiometric temperature. To overcome thermal heterogeneity in the urban area, the measured radiometric temperature was transformed into the equivalent temperature of an upward longwave radiation flux. The equivalent temperature was found to provide an effective parameterization of the radiometric temperature. The daytime average of the resulting ln(z_T/z₀) was 10, where z₀ is the aerodynamic roughness length. This result generally agrees with previous studies; however, the anthropogenic heat is a large uncertainty, which could cause an error at least 240% in z_T.     

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.