不同下垫面空气动力学参数的研究

文中利用中国科学院沙漠研究所与日本国家农业环境技术研究所合作于1990—1994年在中国内蒙古自治区奈曼市半干旱地区沙丘和植被区下垫面观测的微气象数据,根据Monin-Obukhov相似性理论,计算了重度干扰草原、中度干扰草原、轻度干扰草原、无干扰草原、沙丘、沙丘内地、草地、稻田、小麦田、大豆田和玉米田11种下垫面的空气动力学参数粗糙度长度z0,零平面位移d,摩擦速度u*,并分析了它们与水平风速u和Richardson数的关系,比较了不同人为干扰草原生态系统条件下的空气动力学特征。结果表明:地表生物量和覆盖率随着人为干扰强度的增加而减少。不同人为干扰下垫面的粗糙长度与生物量和植被高度以及地表起伏程度有着密切关系;Richardson数也是其影响因子。风速、粗糙度都与摩擦速度成正相关,但对于不同下垫面有所不同,从中可以看到草地对沙漠化有一定的防治作用。同一种下垫面不同时期的空气动力学参数也存在差异。这些结果对建立陆面过程和区域气候模式具有重要的意义。     

半干旱区榆中地表粗糙度年变化及影响机理

利用2006年6月—2010年12月兰州大学半干旱气候与环境观测站 (SACOL) 观测资料,分析了黄土高原自然植被下垫面榆中地表粗糙度时间变化特征,考虑到地形、植被物理特征以及降水和热力条件的影响,分析了东南风向和西北风向粗糙度年变化规律及其影响机理,并分别给出这两个风向归一化粗糙度与时间的拟合关系式。研究发现:对于非均一下垫面,由于地形起伏和下垫面植被差别造成的不同风向粗糙度差异显著。选取东南风向和西北风向,这两个风向的地表粗糙度无论是量级还是年变化特征都有很大差别,且由于地形和植被的差别,东南风向粗糙度年变化趋势与稳定度年变化趋势一致,粗糙度与稳定度存在一定相关关系,而西北风向粗糙度年变化趋势与降水量年变化趋势一致,粗糙度与降水量相关性较好。     

Wave-dependence of friction velocity,roughness length,and drag coefficient over coastal and open water surfaces by using three databases

The parameterization of friction velocity, roughness length, and the drag coefficient over coastal zones and open water surfaces enables us to better understand the physical processes of air-water interaction. In context of measurements from the Humidity Exchange over the Sea Main Experiment (HEXMAX), we recently proposed wave-parameter dependent approaches to sea surface friction velocity and the aerodynamic roughness by using the dimensional analysis method. To extend the application of these approaches to a range of natural surface conditions, the present study is to assess this approach by using both coastal shallow (RASEX) and open water surface measurements (Lake Ontario and Grand Banks ERS-1 SAR) where wind speeds were greater than 6.44 m s^-1. Friction velocities, the surface aerodynamic roughness, and the neutral drag coefficient estimated by these approaches under moderate wind conditions were compared with the measurements mentioned above. Results showed that the coefficients in these approaches for coastal shallow water surface differ from those for open water surfaces, and that the aerodynamic roughness length in terms of wave age or significant wave height should be treated differently for coastal shallow and open water surfaces.     

Bulk Formulation of the Surface Heat Flux

An interpretive literature survey examines different approachesfor applying the bulk aerodynamic formulato predict the surface heat flux. The surface heat flux is often predicted in terms of the surface radiation temperature, which is also used to predict the upward longwave radiation and the heat flux into the soil. In models, the thermal roughness length based on the surface radiation temperature (radiometric roughness length) is often specified to be smaller than the roughness length for momentum for a number of distinct reasons. The definition of the radiometric roughness length depends on the way that the surface temperature is measured, the choice of stability functions and displacement height and inclusion of any additional resistances.Using airborne eddy correlation data collected over eight different sites including bare soil, crops and grassland and several types of forests, the radiometric roughness length is found to vary by orders of magnitude in a manner that is difficult to formulate. Alternatively, we evaluate the approach where the thermal roughness length is equated with the better behaved roughness length for momentum and the corresponding aerodynamic surface temperature is modelled in terms of the surface radiation temperature, solar radiation, and vegetation index. The influence of wind speed and soil moisture on the difference between the aerodynamic and surface radiation temperatures is also examined.     

Response of Nebraska Sand Hills Natural Vegetation to Drought, Fire, Grazing, and Plant Functional Type Shifts as Simulated by the Century Model

The Nebraska Sand Hills exist in a semi-arid climatic environment and the land surface is grassland growing on sandy soils. These soils have been periodically active throughout the Holocene, but are currently stabilized by the vegetation. However, a shift in climate could cause grassland death and eventual sand dune remobilization. Our studies used the CENTURY nutrient cycling and ecosystem model to investigate the impacts of drought, plant functional type, fire, grazing, and erosion on Nebraska Sand Hills vegetation and dune stability. Fire and grazing alone had little impact on the vegetation, but when combined with mild drought, biomass decreased. Overall biomass increased if one plant functional type was allowed to dominate the ecosystem. Addition of as little as 1 mm of erosion per year under current climate conditions decreases vegetation as much as a drought 20 percent drier than the worst of the 1930s drought years in Nebraska.     

Aerodynamic drag coefficient and roughness length for three seasons over a tropical western Indian station

Land Surface Processes Experiment (LASPEX) was conducted over semi-arid region of western India in 1997. As a part of this program, wind and temperature observations were taken using slow as well as fast response sensors over a semi-arid station Anand (22°35′N, 72°55′E) situated in Gujarat state of India. Turbulent parameters such as drag coefficient and sensible heat flux were estimated using eddy correlation method and aerodynamic roughness length was estimated using wind profiles. The analysis has been carried out for the data representing summer, monsoon and winter seasons. It was found that the wind speed does not exceed 5 ms^− 1 during the observational period considered in this study. Relationship of aerodynamic drag coefficient and roughness length with wind speed and stability has been investigated. Aerodynamic roughness length was greater in the stable conditions when the wind speed was low and it reduced drastically during convective conditions. The resulting values of aerodynamic roughness length and drag coefficient for the monsoon period agree well with values reported in literature over Indian subcontinent for homogeneous grass covered surfaces.     

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.     

Diurnal Variations of Thermal Roughness Height over a Grassland

The thermal roughness height associated with the surface radiation temperature has been previously found to vary between different surface types. This study finds that the thermal roughness height varies diurnally even over a homogeneous senescent grassland. The corresponding roughness length for momentum is relatively constant.Both the aerodynamic temperature and the surface radiation temperature are found to be closely related to the air temperature in the middle of the grass canopy. However, the aerodynamic temperature is strongly influenced by the horizontally integrated heat transfer, while the surface radiation temperature represents the integrated thermal emission through the grass depth within the field of view of the radiometer. The aerodynamic temperature is less sensitive to variations and measurement errors in sensible heat flux, wind speed, and air temperature than the thermal roughness height. We find that formulating the aerodynamic temperature in terms of the surface radiation temperature is better posed for use in the bulk formula than using the surface radiation temperature directly and adjusting the thermal roughness length.     

中国北方四类典型下垫面能量分配特征及其环境影响因子研究

中国北方疆域辽阔,多种植被类型交错分布,陆面过程较为复杂。研究该区域地表能量平衡及分配过程,有助于进一步了解其生态系统的生产力及气候变化。本文使用"中国北方协同观测实验"观测资料,选取中国北方四种植被类型——高寒草甸、荒漠草地、玉米农田及半干旱草地,在深入认识地表能量分配特征的基础上,讨论了环境因子对能量分配的影响。结果表明,上述四类典型下垫面能量闭合度在71%~91%,其中半干旱草地最大,玉米农田最小。从各能量分量与净辐射的比率来看,荒漠草地、半干旱草地净辐射主要通过感热形式加热大气,感热通量占净辐射的百分比分别为29%和36%;高寒草甸、玉米农田潜热通量占主导地位,其值分别占净辐射的42%和37%。影响各下垫面水热综合参数Bowen比的气象因子中,半干旱草地受地-气温差(R2=0.64)和饱和水汽压差(R2=0.56)的影响;玉米农田则主要受地-气温差(R2=0.52)的影响;荒漠草地地-气温差(R2=0.37)和平均风速(R2=0.23)是其主要影响因子。Bowen比与NDVI之间存在显著的指数关系(R2=0.74),且随NDVI的增大而减小。     

绿洲-沙漠复合地表条件下的局地和有效粗糙度

地表空气动力学粗糙度是研究水平非均匀条件下陆面过程的基础,文中对黑河试验区各不同下垫面上地表空气动力学粗糙度进行了估算和分析,分别得到了地表空气动力学粗糙度ｚ０ｍ的有效值和局地值,发现各不同典型下垫面的局地ｚ０ｍ值明显不同,尤其在有作物存在时,绿洲的局地ｚ０ｍ值比沙漠、戈壁的值明显要大;而对同一测站,有效空气动力学粗糙度比局地空气动力学粗糙度值明显要大,运用应力分解理论对这种差异进行了解释,认为有效空气动力学粗糙度与复合表面的总应力有关,其代表了较大水平尺度内的表面和较高障碍物的综合作用,当涉及对陆表通量参数化问题时,必须首先确定水平尺度及相应的地表粗糙度参数     

西北地区典型下垫面空气动力学阻抗估算方法改进及其验证

利用2010年古浪近地层野外观测试验取得的干旱区典型下垫面（荒滩、沙漠和农田）的观测资料,对常用的3种空气动力学阻抗算法——Verma—R算法、Liu算法、Thom算法进行了适用性检验,并通过分离Thom算法中的动量和热量粗糙度对其进行改进（下称Thom-2算法）。对比不同下垫面空气动力学阻抗估算结果发现：在沙漠和荒滩下垫面,动力因素较热力因素影响大,空气动力学阻抗与风速的指数关系较农田下垫面明显;与其他3种算法相比,考虑了热量粗糙度的Thom-2算法在各种典型下垫面上的估算结果与观测值相比偏差最小,且在较为均匀的荒滩和沙漠下垫面上估算的空气动力学阻抗精度较农田高;通过比较4种算法估算干旱区不同典型下垫面空气动力学阻抗结果,说明干旱区空气动力学阻抗的估算应考虑热力粗糙度的影响。     

Improvement and Application of the Similarity Saltation Model: Wind-Tunnel Experimental Investigation and Numerical Simulation of the Vertical Sand Mass Flux Distribution in the Saltation Layer

Based upon comparisons between published experimental data and simulated results on the vertical sand flux distribution in the saltation layer, Shao’s similarity saltation model has been greatly improved by correcting the average vertical particle lift-off velocity and using a more suitable universal roughness length. By the improved model, the vertical sand flux profile over the bare, dry and loose uniform sandy surface, which is quite representative of real desert surfaces, can be reproduced very well. Meanwhile, the surface transport rate and the characteristic and average saltation heights have been simulated and analyzed in detail, disclosing their relationships with friction velocity, particle size and roughness length, and the possible underlying mechanisms. Besides, the average particle lift-off velocity and the average mean vertical aerodynamic action upon the ascending particle, which determine the saltation process, are explicitly expressed by parameters involved in the similarity model, and their relationships with friction velocity, particle size and roughness length are also described concisely. The corrected average particle lift-off velocity makes it possible to investigate the characteristic particle trajectory, whose initial velocity equals the average lift-off velocity, so as to estimate the average particle against surface impacting velocity and the average aerodynamic action upon the saltation process.     

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.     

Using Temperature Fluctuation Measurements to Estimate Meteorological Inputs for Modelling Dispersion During Convective Conditions in Urban Areas

We examine the performance of several methods to estimate meteorological inputs for modelling dispersion in urban areas during convective conditions. Sensible heat flux, surface friction velocity and turbulent velocities are estimated from measurements of mean wind speed and the standard deviation of temperature fluctuations at a single level on a tower at two suburban sites and at one urban site in Riverside, California. These estimates are compared with observations made at these sites during a field study conducted in 2007. The sensible heat flux is overestimated in the urban area, while it is underestimated at a suburban site when temperature fluctuations are used in the free convection formulation to estimate heat flux. The bias in heat flux estimates can be reduced through a correction that depends on stability. It turns out that the bias in heat flux estimates has a minor effect on the prediction of surface friction velocity and turbulent velocities. Estimates of sensible heat flux, surface friction velocity and turbulent velocities are sensitive to estimates of aerodynamic roughness length, and we suggest estimating the aerodynamic roughness length through detailed micrometeorological measurements made during a limited field study. An examination of the impact of the uncertainty in estimating surface micrometeorology on concentrations indicates that, at small distances from a surface release, ground-level concentrations computed using estimates of heat flux and surface friction compare well with the those based on observed values: the bias is small and the 95% confidence interval of the ratio of the two concentrations is 1.7. However, at distances much larger than the Obukhov length, this confidence interval is close to 2.3 because errors in both friction velocity and heat flux affect plume spread. Finally, we show that using measurements of temperature fluctuations in estimating heat flux is an improvement on that based on the surface energy balance, even when net radiation measurements are available.     

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.     

Transport rate of drifting snow and the mean wind speed profile

Transport rates, measured by weighing snow blown into a filter fabric trap, were greater over hard snow or ice than for the same wind speed over soft, fresh snow surfaces. Analysis of wind speed profiles from nine blizzards showed that friction between moving particles and the surface was less, and particle speeds were greater over hard surfaces. Transport rates at a given wind speed increased rapidly as aerodynamic roughness decreased in the rough-smooth transition region. Bagnold's theory for bed load transport provided a useful framework for the analysis.     

不同下垫面大孔径闪烁仪观测数据处理与分析

大孔径闪烁仪是近年兴起的测量大尺度（500 m~10 km）地表通量的仪器。北京师范大学等单位分别于2002年、2004年在北京昌平小汤山开展了大孔径闪烁仪短期观测实验, 2006年6月又在北京密云建立了长期观测站。利用这些数据, 对大孔径闪烁仪观测数据进行处理与分析, 结果表明:闪烁仪光径高度和风速是影响观测显热通量的关键因子。当地表粗糙元的高度变化相对于光径高度不可忽略时, 零平面位移需要精确确定。波文比在湿润地表需要准确确定, 而气温、气压和动力学粗糙度则为不敏感因子。计算中所需的大气稳定度可用理查孙数判断, 也可借助日出日落时间或净辐射观测值确定。稳定条件下的普适函数目前无统一表达式, 可采用仪器说明书推荐的函数。通过几个站点闪烁仪观测显热通量与涡动相关仪测量值的比较表明:大孔径闪烁仪在均匀和非均匀地表都能得到合理的显热通量观测值。     

渤海西岸空气动力学粗糙度特征分析

利用2009年12月至2010年11月渤海西岸大港风能塔观测得到的风速、风向和温度梯度资料,计算了风能塔周围空气动力学粗糙度Z0。对比分析了风廓线法和风速标准差法的计算结果,讨论了风速标准差法的适用性以及下垫面空间非均一性对粗糙度的影响。结果表明：风能塔三个主风向上的平均粗糙长度为：0.1319m(N),0.0386m(SE)和0.0182m(SW)。通过严格条件限制,利用风速标准差法可以得到与风廓线法相同的计算效果。同一方位上,利用10-30 m资料计算得到的粗糙度长度的季节差异相对较小,而50-70 m的计算结果差异相对较大。在相同高度上,粗糙度长度计算结果亦存在明显季节差异,反映了随季节变化的植被对地表粗糙度的影响。     

奈曼沙漠化对草地微气象特性影响的研究

为探索沙漠化发生机制，在内蒙古奈曼的流动沙丘，干草地和湿草地进行了微气象观测，并用热量平衡法分析观测数据。观测结果表明，沙漠化发生后，反射率增加，净 辐射减少；而且热量平衡中土壤热通量所占份额增加，潜热通量所占份额减少。干、湿草地的风速廓线在白天和夜间都近于对数分布。而沙丘的风速廓线只在晚上近于对数分布。白天沙丘近地表风速垂直切变较大，从而对沙丘表面产生较强切应力。曳力系数白天低于夜里，沙丘低于干、湿草地。     

深圳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, 结果合理可信。研究表明空气动力学参数受下垫面非均匀性、植株柔软性、气流来向、风速等的共同影响。