松林冠层大气湍流结构观测研究

为了更好地了解大气与物质和能量交换，用二层三维超声风速温度仪测量了重庆市郊松林风速和温度脉动值。采样速率为每秒１次和１１次，数字量记录。计算了湍流动最通量和热通量的日变化，以及湍流风速的统计量和功率谱。结果表明，冠层上动量通量向下传输，而冠层内大多向上传输；冠层内湍流风速ｕ的三阶矩平均值大于零，ｗ的小于零；冠层上无因次湍流风速分量ｕ，ｗ标准差和无因次湍流动能耗散率与局地Ｍｏｎｉｎ－Ｏｂｕｋｈｏｖ长     

植冠层结构对其湍流场影响的数值试验

根据湍流应力模式的模拟结果，讨论了植被冠层中湍流能量的平衡特征，并通过数值试验就冠层结构对流场结构的影响进行了分析。     

土壤-植被-大气系统水分能量传输模拟和验证

本文建立了包括地下水的土壤－植被－大气系统水分能量传输综合模型,并对模型中冠层动量湍流交换反梯度传输现象的描述进行了改进。在此基础上,依据冠层内动量和热量的交换以及辐射传输过程,同时求解地表、冠层能量平衡方程,进而模拟饱和－非饱和土壤的水热传输。用浅地下水地区冬小麦田间试验资料对模型进行验证,结果表明,系统能量平衡各分量和土壤含水量的模拟与观测结果相当一致。模型敏感性分析发现,叶面积指数对总蒸散量的影响随叶面积指数的增加而逐渐减弱;叶片最小气孔阻力对总蒸散量的影响,在该阻力较小时更显著;地下水位对蒸散量的影响在它小于１．５ｍ时不显著,而在１．５～１．７５ｍ之间时,蒸散减小较快,主要由于土壤蒸发减小显著,冠层蒸腾稍有增加。     

植物冠层动量交换特征的实验研究

使用湍流梯度测试资料，对植物冠层动量交换特征进行了详细研究，结果表明:森林冠层内惯性副区能谱曲线仍可用幂指数描述，但斜率比-2/3更负；森林冠层内湍流尺度有变小的趋势；森林上层的耗散系数比下层大；由植被吸收引起动量及动量通量随冠层深度增加而明显减小；冠层下层的动量通量和耗散系数分别与上层的量有好的正相关；森林冠层内耗散系数和动量通量随大气稳定度有明显变化。     

植被冠层湍流代数应力模式

从植被冠层湍流时空平均雷诺应力方程出发，根据植被冠层内有关参数化的准则，建立了湍流代数应力模式，并用玉米和小麦两种作物地的实测资料对模式进行了验证。     

城市建筑动力学效应对对流边界层影响的敏感性试验

本文将大涡模拟应用于城市对流边界层(CBL)湍流结构和流场特征的研究,在大涡模式中,拖曳系数取与建筑物高度及建筑物高度标准差有关的表达式以考虑次网格建筑物对风速和湍流动能(TKE)的面积平均影响.模拟结果表明,由于城市建筑物对气流的拖曳作用,使建筑物冠层及整个CBL内风速大幅度减小,城市冠层内部风速减小尤为明显,在夹卷层内,风速有一明显的跃变.在边界层中部对流运动已经发展成为较强的热泡,城市建筑物的动力学效应使热泡的水平尺度增大,CBL内平均上升气流速度和下沉气流速度减小,同时使CBL中上升气流所占比例比平坦地面增大.城市建筑物使CBL低层热通量、动量通量、速度方差和位温方差明显增大,但对近地层高度以上的湍流量影响不大.     

山地边界层中场变量的动力学研究

采用与实测较接近的二次函数来表达Ｅｋｍａｎ层中的湍流粘性系数Ｋ，在圆形气压场条件下，求得了山地上空边界层中的风速，进而求得散度、涡度和垂直速度等场变量随高度的分布。并作图分析了这些场变量的一些动力学特征。改进了以往在求解析解时，略去运动方程中湍流粘性力项中的关于高度的一阶导数项，以及取山坡面上风速为零作下边界条件等欠合理欠精确的做法。所求得的风速、散度、涡度和垂直速度均用简单的初等函数表示出来，有助于边界层参数化和深化对边界层动力学的认识。     

植物冠层内物质交换特征的实验研究

使用湍流和SO2通量梯度测试资料，对植被冠层内物质交换特征进行了详细研究，揭示了一些有意义的结果: 植物的生物过程对物质的吸收作用非常明显； 在冠层内物质通量是随冠层深度加大而明显减小；对高的植被，上层沉积速度(Vgu)大于下层 (Vgd)； 冠层内的沉积速度(Vg) 表现出明显的日变化，冠层内物质的生物吸收作用与太阳总辐射量有直接联系；冠层的Vg与速度尺度V＊和平均风成正变关系；森林的Vg比麦地小；新的Vg理论公式能更好地预测重庆森林和麦地的结果。     

城市冠层上下大气湍流特征分析

利用兰州市榆中县城市冠层架设的3台涡动相关仪(EC)观测大气湍流资料,分析了城市冠层上下不同下垫面湍流通量和不同风向范围内湍流动能变化特征,之后对城市冠层上下3台EC观测湍流动能最小、最大分布方向上分别对应的最小最大湍流动能的风速谱进行了研究,进一步检验了局地相似理论在城市冠层上下的适用性。结果表明:(1)城市冠层之上水泥、砖石等构成下垫面和城市冠层之内草坪下垫面观测感热通量、摩擦速度较为接近,观测潜热通量、CO2通量在白天差异明显。(2)城市冠层之上的湍流动能总体上大于冠层之内,冠层之上气流来向的上风向较为开阔时湍流动能较大,而冠层之内气流来向的上风向为街道口时湍流动能较大。(3)城市冠层之上的湍涡尺度大于冠层之内,城市冠层小尺度湍涡风速谱在惯性副区基本符合-2/3次方关系,且准各向同性,大尺度湍涡风速谱在惯性副区不符合-2/3次方关系,且各向异性。(4)不稳定层结下,城市冠层上下无量纲速度方差与稳定度基本满足1/3次方局地相似关系,稳定层结下不满足;城市冠层上下无量纲温度、湿度、CO2浓度方差在所有层结下均不满足-1/3次方局地相似关系。(5)近中性层结下,城市冠层上下u、v、w方向无量纲速度方差分别为3.52,3.03,1.49和2.62,2.22,1.50。     

城市冠层中湍流运动的统计特征

对１９９７ 年夏天和冬天北京湍流运动的各种统计特征量进行了初步的统计分析。结果表明,城市冠层中湍流运动的各种统计特征量与平坦下垫面条件下边界层湍流运动的相比, 有不同的地方也有相似的地方; 无论白天还是夜晚, 垂直方向的湍流强度和湍流脉动风速标准差均小于水平方向的, 水平方向的相应湍流特征量则总是接近相等; 城市冠层中湍流脉动强度和标准差几乎均大于平坦下垫面边界层的; 平均风速ｕ≥１ ｍ ／ｓ 时的湍流统计特征量与ｕ＜ １ ｍ ／ｓ 时的有所不同; 城市冠层的阻力系数较大, 可达００６２５,Ｐａｎｏｆｓｋｙ 等提出的公式σｗ ／ｕ＊ ＝ １３ （１－ ３ｚ／Ｌ）１／３在城市冠层中并不适用。     

Edge Flow and Canopy Structure: A Large-Eddy Simulation Study

Sharp heterogeneities in forest structure, such as edges, are often responsible for wind damage. In order to better understand the behaviour of turbulent flow through canopy edges, large-eddy simulations (LES) have been performed at very fine scale (2 m) within and above heterogeneous vegetation canopies. A modified version of the Advanced Regional Prediction System (ARPS), previously validated in homogeneous conditions against field and wind-tunnel measurements, has been used for this purpose. Here it is validated in a simple forest-clearing-forest configuration. The model is shown to be able to reproduce accurately the main features observed in turbulent edge flow, especially the “enhanced gust zone” (EGZ) present around the canopy top at a few canopy heights downwind from the edge, and the turbulent region that develops further downstream. The EGZ is characterized by a peak in streamwise velocity skewness, which reflects the presence of intense intermittent wind gusts. A sensitivity study of the edge flow to the forest morphology shows that with increasing canopy density the flow adjusts faster and turbulent features such as the EGZ become more marked. When the canopy is characterized by a sparse trunk space the length of the adjustment region increases significantly due to the formation of a sub-canopy wind jet from the leading edge. It is shown that the position and magnitude of the EGZ are related to the mean upward motion formed around canopy top behind the leading edge, caused by the deceleration in the sub-canopy. Indeed, this mean upward motion advects low turbulence levels from the bottom of the canopy; this emphasises the passage of sudden strong wind gusts from the clearing, thereby increasing the skewness in streamwise velocity as compared with locations further downstream where ambient turbulence is stronger.     

Characteristics of canopy turbulence over a deciduous forest on complex terrain

Canopy turbulence plays an important role in mass and energy exchanges at the canopy-atmosphere interface. Despite extensive studies on canopy turbulence over a flat terrain, less attention has been given to canopy turbulence in a complex terrain. The purpose of this study is to scrutinize characteristics of canopy turbulence in roughness sublayer over a hilly forest terrain. We investigated basic turbulence statistics, conditionally sampled statistics, and turbulence spectrum in terms of different atmospheric stabilities, wind direction and vertical structures of momentum fluxes. Similarly to canopy turbulence over a homogeneous terrain, turbulence statistics showed coherent structure. Both quadrant and spectrum analysis corroborated the role of intermittent and energetic eddies with length scale of the order of canopy height, regardless of wind direction except for shift of peak in vertical wind spectrum to relatively high frequency in the down-valley wind. However, the magnitude of the momentum correlation coefficient in a neutral condition was smaller than typical value over a flat terrain. Further scrutiny manifested that, in the up-valley flow, temperature skewness was larger and the contribution of ejection to both momentum and heat fluxes was larger compared to the downvalley flow, indicating that thermal instability and weaker wind shear in up-valley flow asymmetrically affect turbulent transport within the canopy.     

Application of transilient turbulent theory to study interactions between the atmospheric boundary layer and forest canopies

The new Forest-Land-Atmosphere ModEl called FLAME is presented. The first-order, nonlocal turbulence closure called transilient turbulence theory (Stull, 1993) is applied to study the interactions between a forested land-surface and the atmospheric boundary layer (ABL). The transilient scheme is used for unequal vertical grid spacing and includes the effects of drag, wake turbulence, and interference to vertical mixing by plant elements. Radiation transfer within the vegetation and the equations for the energy balance at the leaf surface have been taken from Norman (1979). Among others, the model predicts profiles of air temperature, humidity and wind velocity within the ABL, sensible and latent heat fluxes from the soil and the vegetation, the stomata and aerodynamic resistances, as well as profiles of temperature and water content in the soil. Preliminary studies carried out for a cloud free day and idealized initial conditions are presented. The canopy height is 30 m within a vertical domain of 3 km. The model is able to capture some of the effects usually observed within and above forested areas, including the relative wind speed maximum in the trunk space and the counter gradient-fluxes in the lower part of the plant stand. Of special interest is the determination of the location and magnitude of the turbulent mixing between model layers, which permits one to identify the effects of large eddies transporting momentum and scalar quantities into the canopy. A comparison between model simulations and field measurements will be presented in a future paper.     

Turbulence structure in a deciduous forest

Three-dimensional wind velocity components were measured at two levels above and at six levels within a fully-leafed deciduous forest. Greatest shear occurs in the upper 20% of the canopy, where over 70% of the foliage is concentrated. The turbulence structure inside the canopy is characterized as non-Gaussian, intermittant and highly turbulent. This feature is supported by large turbulence intensities, skewness and kurtosis values and by the large infrequent sweeps and ejections that dominate tangential momentum transfer. Considerable day/night differences were observed in the vertical profiles of the mean streamwise wind velocity and turbulence intensities since the stability of the nocturnal boundary layer dampens turbulence above and within the canopy.     

Influence of non-flat terrain and wind direction shear on canopy turbulence

We have analyzed eddy covariance data collected within open canopy to investigate the influence of non-flat terrain and wind direction shear on the canopy turbulence. The study site is located on non-flat terrain with slopes in both south-north and east-west directions. The surface elevation change is smaller than the height of roughness element such as building and tree at this site. A variety of turbulent statistics were examined as a function of wind direction in near-neutral conditions. Heterogeneous surface characteristics results in significant differences in measured turbulent statistics. Upwind trees on the flat and up-sloping terrains yield typical features of canopy turbulence while upwind elevated surface with trees yields significant wind direction shear, reduced u and w skewness, and negligible correlation between u and w. The directional dependence of turbulence statistics is due that strong wind blows more horizontally rather than following terrain, and hence combination of slope related momentum flux and canopy eddy motion decreases the magnitude of Sk _w and r _uw for the downslope flow while it enhances them for the upslope flow. Significant v skewness to the west indicates intermittent downdraft of northerly wind, possibly due to lateral shear of wind in the presence of significant wind direction shear. The effects of wind direction shear on turbulent statistics were also examined. The results showed that correlation coefficient between lateral velocities and vertical velocity show significant dependence on wind direction shear through change of lateral wind shear. Quadrant analysis shows increased outward interaction and reduced role of sweep motion for longitudinal momentum flux for the downslope flow. Multi-resolution analysis indicates that uw correlation shows peak at larger averaging time for the upslope flow than for the downslope flow, indicating that large eddy plays an active role in momentum transfer for the upslope flow. On the other hand, downslope flow shows larger velocity variances than other flows despite similar wind speed. These results suggest that non-flatness of terrain significantly influences on canopy-atmosphere exchange.     

台风莫拉克 (0908) 影响期间近地层风特性

通过对台风莫拉克 (0908) 影响范围内的33座测风塔观测资料的分析可知:台风莫拉克越靠近陆地,风场的非对称性越明显,其行进方向的左侧测风塔风向呈逆时针旋转,右侧测风塔风向顺时针旋转。在远离台风莫拉克的地方风向稳定,湍流强度变化较平稳;在台风莫拉克登陆点附近,风向、风速和湍流强度均会出现突变。台风莫拉克影响期间,湍流强度与风速的关系未出现IEC标准曲线那样随风速增大稳定减小,其I₁₅达B级和A级及以上的平均湍流强度会在风速7~17 m·s-1形成一个峰值;无论南风或北风,风速越大,各层湍流强度差异趋于减小,同等风速、高度的湍流强度偏南大风均大于偏北大风。位于台风莫拉克登陆点北侧测风塔湍流强度随风速的增加先减小后增大,最终各高度全部超过IEC标准A级曲线,而位于南侧测风塔湍流强度随风速的变化比北侧小,并随风速增大趋于标准A级曲线;另外北侧测风塔湍流强度大于南侧,且各高度偏北大风湍流强度之间的差异比南侧相应风向明显,表明北侧垂直方向的扰动更强。台风莫拉克阵风系数为1.2~1.7,其随高度变化与地形有关,一般情况下随高度升高而减小,在复杂地形条件下不符合随高度升高减小的规律。     

Large-Eddy Simulation of Flow Over a Vegetation-Like Canopy Modelled as Arrays of Bluff-Body Elements

Turbulent flow over a vegetation canopy under neutral atmospheric conditions is investigated using large-eddy simulation. Each model tree, which consists of a sphere-shaped tree crown and a cylindrical trunk, is fully resolved. The resulting turbulence statistics and the drag force on the vegetation agree well with measurements from the corresponding wind-tunnel experiment described by Böhm et al. (Boundary-Layer Meteorol, 146:393–419, 2013). Statistically, this kind of model canopy exhibits both vegetation and bluff-body-flow characteristics. The time-averaged flow skims over the top of the underlying canopy, forming a low-momentum recirculation zone on the lee-side of the bluff elements, which causes significant dispersive stress within the canopy layer. Two other numerical representations of vegetation canopies, referred to as the drag-element and drag-crown approaches, have also been developed to assess the performance of simulations. Turbulence statistics suggest that the canopy shear layer interferes with wakes behind stems and crowns. The drag-crown approach yields better agreement between numerical results and experimental measurements than does the traditional drag-element approach, thus providing a promising numerical model for simulating canopy turbulence.     

An Analytical Model for Mean Wind Profiles in Sparse Canopies

Existing analytical models for mean wind profiles within canopies are applicable only in dense canopy scenarios, where all momentum is absorbed by canopy elements and, hence, the effect of the ground on turbulent mixing is not important. Here, we propose a new analytical model that can simulate mean wind profiles within sparse canopies under neutral conditions. The model adopts a linearized canopy-drag parametrization and a first-order turbulence closure scheme taking into account the effects of both the ground and canopy elements on turbulent mixing. The resulting wind profile within a sparser canopy appears to be more like a logarithmic form, with the no-slip condition at the ground being satisfied. The analytical solution converges exactly to the standard surface-layer logarithmic wind profile in the case of zero canopy density (i.e., no-canopy scenario) and tends to be an exponential wind profile for a dense canopy; this feature is unique compared with existing analytical models for canopy wind profiles. Results from the new model are in good agreement with those from laboratory experiments and numerical simulations.     

Particle Image Velocimetry Measurements of Turbulent Flow Within Outdoor and Indoor Urban Scale Models and Flushing Motions in Urban Canopy Layers

We investigate the spatial characteristics of urban-like canopy flow by applying particle image velocimetry (PIV) to atmospheric turbulence. The study site was a Comprehensive Outdoor Scale MOdel (COSMO) experiment for urban climate in Japan. The PIV system captured the two-dimensional flow field within the canopy layer continuously for an hour with a sampling frequency of 30 Hz, thereby providing reliable outdoor turbulence statistics. PIV measurements in a wind-tunnel facility using similar roughness geometry, but with a lower sampling frequency of 4 Hz, were also done for comparison. The turbulent momentum flux from COSMO, and the wind tunnel showed similar values and distributions when scaled using friction velocity. Some different characteristics between outdoor and indoor flow fields were mainly caused by the larger fluctuations in wind direction for the atmospheric turbulence. The focus of the analysis is on a variety of instantaneous turbulent flow structures. One remarkable flow structure is termed ‘flushing’, that is, a large-scale upward motion prevailing across the whole vertical cross-section of a building gap. This is observed intermittently, whereby tracer particles are flushed vertically out from the canopy layer. Flushing phenomena are also observed in the wind tunnel where there is neither thermal stratification nor outer-layer turbulence. It is suggested that flushing phenomena are correlated with the passing of large-scale low-momentum regions above the canopy.     

A spectral and lag-correlation analysis of turbulence in a deciduous forest canopy

The processes influencing turbulence in a deciduous forest and the relevant length and time scales are investigated with spectral and cross-correlation analysis. Wind velocity power spectra were computed from three-dimensional wind velocity measurements made at six levels inside the plant canopy and at one level above the canopy. Velocity spectra measured within the plant canopy differ from those measured in the surface boundary layer. Noted features associated with the within-canopy turbulence spectra are: (a) power spectra measured in the canopy crown peak at higher wavenumbers than do those measured in the subcanopy trunkspace and above the canopy; (b) peak spectral values collapse to a relatively universal value when scaled according to a non-dimensional frequency comprised of the product of the natural frequency and the Eulerian time scale for vertical velocity; (c) at wavenumbers exceeding the spectral peak, the slopes of the power spectra are more negative than those observed in the surface boundary layer; (d) Eulerian length scales decrease with depth into the canopy crown, then increase with further depth into the canopy; (e) turbulent events below crown closure are more correlated with turbulent events above the canopy than are those occurring in the canopy crown; and (f) Taylor's frozen eddy hypothesis is not valid in a plant canopy. Interactions between plant elements and the mean wind and turbulence alter the processes that produce, transport and remove turbulent kinetic energy and account for the noted observations.