大气层结对青藏高原近地层湍流特征的影响

本文利用1986年夏季青藏高原试验（TIPMEX-86）期间获取的资料，计算分析了该年6月和7月拉萨和那曲近地层摩擦速度u*、特征温度θ*、Monin-Obukhov长度L及理查逊数Ri在不同大气层结及不同高度的分布特征。分析结果表明:不同大气层结下，湍流运动形式是不同的，湍流热力场的变化远远大于湍流动力场。就空间变化而言，θ*随高度变化较大，而u*的变化较小。对大气层结而言，在稳定条件下，u*、θ*随高度变化较大；在不稳定条件下，θ*、u*随高度变化较小。在不稳定层结相近的状况下，高原开阔草地上的各种湍流参数，接近或大于北京城郊非均匀地形的数值，而在很稳定的层结下，则高原的参数值小于北京。     

大气边界层非常定下沉急流和阵风的起沙机理

首先用动量原理导出水平沙尘通量Q的Bagnold型公式,由此可以反推出使Q达到统计平衡态所需的气流水平动量下传的作用时间.其次导出准常定下沉急流情况下的公式Q,此时公式中的摩擦速度u*应由u* *代替,其中u4* *=u4*m+u4*,u*m和u*各指由下沉平均急流和脉动的动量垂直通量所对应的特征速度(习惯上称为摩擦速度).其后指出,在研究风蚀和计算沙尘通量时,应以1 min做统计平均为宜,否则会过低估算了起沙尘的时段长度和起沙尘量.最后指出,不同频率的扰动和脉动的起沙尘效率是不同的,以适中的频率效率最高,从而在非常定急流作用下,应用u* *e代替u* *,其中u4* *e=u4*m+u4*e+(1+A4) u4*g,而u*m,u*e和u*g则各为由平均急流、湍流脉动和阵风所导致的相应摩擦速度,A为阵风增益系数.     

MM5和ETA相似理论近地层方案对农田下垫面通量模拟比较研究

近地层湍流通量计算对于中尺度数值模式有重要意义,湍流通量的参数化是当前大气边界层研究的重要课题之一.选择青藏高原东缘大理观象台边界层通量观测系统,离线测试了WRF区域模式中的两种常用的近地层参数化方案(MM5相似理论非迭代方案A和ETA相似理论迭代方案B),并将参数化方案计算结果与边界层铁塔涡动相关法的观测值进行对比分析.在大理观象台观测场不同植被随季节交替的状况下,根据边界层铁塔4层高度风速拟合,发现近地层空气动力学粗糙度随季节变化特征明显.将拟合的空气动力学粗糙度输入模式参数化方案进行通量计算.结果表明:稳定度是影响近地层参数化方案精度的重要因素,在不稳定条件下方案B低估了动量通量,方案A优于方案B,而在稳定条件下方案A低估了动量通量,方案B优于方案A,两种方案总体来看误差不大.对于大理边界层通量观测场农田植被交替的环境条件,不同季节下垫面植被类型的差异,以及植被的稀疏对近地层参数化方案湍流通量计算结果的精度有显着影响.方案B考虑了空气动力学粗糙度z₀和热量粗糙度z_0h的差异,不稳定条件下感热通量计算结果在裸土或稀少植被条件下明显优于方案A.针对方案B不稳定条件下感热通量计算结果在裸土下垫面仍出现高估的现象,使用了Zeng等1998年提出的用辐射地表温度订正裸土下垫面感热能量方法后,计算结果也有明显改善.     

海气界面动量通量算法的改进

COARE（Coupled Ocean-Atmosphere Response Experiment）算法是国际上较先进的计算海气界面通量的算法。最新的COARE 3.0算法包含TY01方案和O02方案两种考虑真实海浪状态的海面空气动力学粗糙度方案;当有效波高和谱峰周期缺省时,利用Taylor01风浪特征量参数化方案可对它们进行参数化。在此基础上引入自主设计的WHP（Wind-wave significant wave Height and dominant wave Period）风浪特征量参数化方案和三种海面空气动力学粗糙度方案,即SCOR方案、GW03方案、PYP07方案以及更简单的Andreas12、Vickers15摩擦速度算法,利用NDBC（National Data Buoy Center）浮标数据、解放军理工大学风浪流水槽实验数据,针对中高风速条件（10 m/s≤U₁₀≤25 m/s）比较上述方案对摩擦速度的预测效果。结果表明:WHP方案在COARE 3.0算法中的应用效果优于Taylor01方案,且新引入的SCOR、GW03、PYP07海面粗糙度方案对风浪特征量的观测误差敏感性更小、稳定度更高;水槽实测数据的对比结果表明,WHP方案结合SCOR海面粗糙度方案计算的摩擦速度与实测值最接近;引入其他三种实测数据的检验结果表明,原始COARE 3.0算法会低估摩擦速度,而WHP方案结合SCOR海面粗糙度方案能更准确地预测摩擦速度随10 m风速的增长趋势。      

MM5和ETA相似理论近地层方案对农田下垫面通量模拟比较研究

近地层湍流通量计算对于中尺度数值模式有重要意义, 湍流通量的参数化是当前大气边界层研究的重要课题之一。选择青藏高原东缘大理观象台边界层通量观测系统, 离线测试了WRF区域模式中的两种常用的近地层参数化方案(MM5相似理论非迭代方案A和ETA 相似理论迭代方案B), 并将参数化方案计算结果与边界层铁塔涡动相关法的观测值进行对比分析。在大理观象台观测场不同植被随季节交替的状况下, 根据边界层铁塔4层高度风速拟合, 发现近地层空气动力学粗糙度随季节变化特征明显。将拟合的空气动力学粗糙度输入模式参数化方案进行通量计算。结果表明:稳定度是影响近地层参数化方案精度的重要因素, 在不稳定条件下方案B低估了动量通量, 方案A优于方案B, 而在稳定条件下方案A低估了动量通量, 方案B优于方案A, 两种方案总体来看误差不大。对于大理边界层通量观测场地农田植被交替的环境条件, 不同季节下垫面植被类型的差异, 以及植被的稀疏对近地层参数化方案湍流通量计算结果的精度有显著影响。方案B考虑了空气动力学粗糙度z₀和热量粗糙度z_0h的差异, 不稳定条件下感热通量计算结果在裸土或稀少植被条件下明显优于方案A。针对方案B不稳定条件下感热通量计算结果在裸土下垫面仍出现高估的现象, 在使用了(Zeng, et al, 1998)提出的对于使用辐射地表温度在裸土下垫面时的订正方法后, 计算结果也有明显改善。     

高寒草甸湍流特征量的季节变化特征

利用青藏高原东北边缘的若尔盖高寒湿地生态系统研究站玛曲高寒草甸观测场提供的2010年近地面层梯度和涡动观测资料,基于Monin-Obukhov相似理论(MOST)分析了该地区风速分量、温度、水汽和CO2的归一化标准差随稳定度参数z/L在均匀平坦地表的变化关系及其季节变化特征。结果表明,风速分量σu/u*、σv/u*和σw/u*与稳定度参数z/L呈1/3次方关系,并随|z/L|增大而增大,在近中性(|z/L|0.05)时趋于常数,且σu/u*σv/u*σw/u*;水平方向风速的归一化标准差随空气动力学粗糙度z0m的增大而减小,σw/u*无明显季节变化,可近似为常数;温度、水汽和CO2的归一化标准差在不稳定条件下与z/L呈-1/3次方关系,归一化值随|z/L|的减小而增大,在中性条件下σq/q*σc/C*σT/T*。σT/T*、σq/q*和σc/C*的季节变化分别与感热、潜热和CO2通量的季节变化相反,且主要决定于特征尺度T*、q*和C*。其中感热通量的季节变化与潜热通量、CO2通量相反,且年变化幅度最小,相应地,σT/T*的季节变化与σc/C*和σq/q*相反,且年变化幅度最小。     

北京市郊区及城区边缘的大气湍流结构特征

为了比较城市边缘地区和郊区大气边界层中湍流结构的差异,本文分析了1986年春、夏两季在北京市城市边缘地区和郊区两处大气湍流观测结果.结果表明,由于城市边缘地区的下垫表面具有较高数值的粗糙度,使其近地面层大气受到较强的动力影响,在其湍流速度分量谱曲线;无量纲化湍流速度分量标准偏差,σ_u/u*,σ_v/u*和σ_w/u*;以及无量纲化湍流耗散率φ_ε等湍流特征量在近中性层结条件下显示出与大气稳定度的关系比较微弱.     

GRAPES_Meso区域模式积云计算方案引进及预报效果检验

云量在模式中的准确模拟既是一个重点问题,也是一个难点问题。经评估发现,在GRAPES_Meso（Global/Regional Assimilation and PrEdiction System）模式中存在模拟云量偏少的问题,而ECMWF模式中的积云云量计算方案对云量模拟具有一定优势。鉴于此,本文借鉴ECMWF积云云量计算方案对GRAPES中尺度模式中的云量计算方案进行了适当改进,并通过影响试验模拟三维云量和地面温度来验证其改进效果。结果表明:较改进前,改进云量方案模拟的云量增多,单层云量更加准确,垂直方向云层结构改进效果明显;地面温度与观测值对比误差减小。综合分析模拟结果,认为改进云量计算方案在GRAPES模式中对云量和地面温度的模拟有较好的改进。     

黑河下游绿洲近地层湍流输送特征研究

分析2003年9月5-27日在黑河下游额济纳绿洲取得的近地面层湍流观测资料,讨论无量纲湍流方差与稳定度参数Z/L的关系.无量纲风速分量σu/u*、σv/u*和σw/u*在不稳定层结与稳定度Z/L满足1/3次方时,符合Monin-Obkuhov相似理论.在近中性时,σu/u*、σw/u*值比平坦下垫面的略小些,σv/u*值比平坦下垫面的略大.无量纲温度脉动方差和湿度脉动方差在对流状态下与Z/L满足-1/3次方,在稳定和不稳定层结下,随|Z/L|减小而增大.地面热量以感热输送为主,感热通量峰值约为250 W/m2,潜热通量的峰值为170 W/m2;夜间潜热通量较小,感热通量则出现负值.动量通量日间平均在0.2-0.3 N/m2,峰值为0.31 N/m2.黑河下游绿洲柽柳林动力学粗糙度比HEIFE戈壁及其他绿洲下垫面的大,中性时z0m=0.47551 m.中性状态下动量整体输送系数CD=22.2×10-3,不稳定层结时平均值CD=10.3×10-3,稳定层结时平均值CD=7.2 × 10-3.热量整体输送系数CH,中性状态下CH=3.2 × 10-3,不稳定层结时平均值为CH=3.3×10-3,稳定层结时平均值为CH=2.7×10-3.     

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

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

一次长江流域梅雨降水中三种云量计算方案的对比研究

采用NCEP分析场,选取2010年梅雨期长江流域的一次降水过程,分别基于Slingo方案、NCAR方案和钱氏方案,利用相对湿度计算云量,并以LAPS(Local Analysis and Prediction System)系统输出的云量分析场作为观测值,分别在高层(400 h Pa)与低层(850 h Pa),从宏观比较与统计分析的角度,与计算结果进行云量大小与区域分布的对比分析。结果表明,三个云量计算方案对云量中心位置的把握均较为准确,但对云量值的计算存在大小不等的误差。NCAR方案计算结果和LAPS输出场最为吻合,能够体现出云量大值区,但区域一般偏大;Slingo方案相较NCAR方案来说略差,但也能较好地描述云带分布;此外,钱氏方案计算出的云量值始终偏小,但其能够较好地描述云带轮廓与云量的分布特征。综合对比结果,NCAR云量计算方案比其余两者更优,且在低层(850 h Pa)表现尤为明显。     

Turbulence Measurements above a Pine Forest

Eddy fluxes of momentum, sensible and latent heat, and turbulence spectra measured over the Thetford Forest during 10 days in the Spring of 1973 are described. The measured total heat flux (H + λE) for 122 20-min periods agreed closely on average with independent estimates from an energy balance method. There was evidence that the energy balance data gave small systematic overestimates of available energy during the hours before noon, compensated by slight underestimates for the remainder of the day. A comparison of measured wind speeds and friction velocities in neutral stability confirmed the validity of the aerodynamic method for estimating momentum fluxes at heights of a few roughness lengths above the canopy. In stable conditions the log-linear wind profile $U=(u_{?t}/k)({? ln}((z-d)/z_{0})+←pha(z-d-z_{0})/L)$ with α = 3.4 ± 0.4 provided a good fit to the data. Spectra in unstable conditions were generally more sharply peaked than those measured by other workers over smoother terrain: differences were less marked in the case of vertical velocity in stable conditions. Temperature spectra in these stable conditions showed high energy at relatively low wavenumbers, and wT cospectra showed a cospectral gap; both of these results were associated with an intermittent sawtooth structure in the temperature fluctuations.     

On the profiles of wind velocity in the roughness sublayer above a coniferous forest

Wind velocity and temperature measurements from a 200 m tower, locatedin a forest near Karlsruhe were used to investigatethe modified profile function of the wind velocity in theroughness sublayer.To avoid determination of the friction velocity we introduced analternative analysis with the expression instead of From the observed Fm* profiles we evaluated the profile function m*. The wind profiles observed under neutral conditions were well representedby a modified non-dimensional profile function with physically based boundary values at the top and at the bottom of theroughness sublayer.The results of our analysis can be used to take into consideration themomentum exchange between the atmosphere and a forest in mesoscaleatmospheric models in a refined way.     

Aerodynamic Roughness Length of Fresh Snow

This study presents the results from a series of wind-tunnel experiments designed to investigate the aerodynamic roughness length z ₀ of fresh snow under no-drift conditions. A two-component hot-film anemometer was employed to obtain vertical profiles of velocity statistics in a zero pressure gradient turbulent boundary layer for flow over naturally deposited snow surfaces. The roughness of these snow surfaces was measured by means of digital photography to capture characteristic length scales that can be related to z ₀. Our results show that, under aerodynamically rough conditions, the mean value of the roughness length for fresh snow is \({\langle{z}_{0}\rangle= 0.24}\) mm with a standard deviation σ(z ₀) = 0.05 mm. In this study, we show that variations in z ₀ are associated with variations in the roughness geometry. The roughness measurements suggest that the estimated values of z ₀ are consistent with the presence of irregular roughness structures that develop during snowfalls that mimic ballistic deposition processes.     

A Simple Method of Estimating Scalar Fluxes Over Forests

A simple aerodynamic-variance method is proposed to fill gaps in continuous CO₂ flux measurements in rainy conditions, when open-path analysers do not function. The method requires turbulent conditions (friction velocity greater than 0.1 ms^–1), and uses measurements of mean wind speed, and standard deviations of temperature and CO₂ concentration fluctuations to complement, and at times replace, eddy-covariance measurements of friction velocity, sensible heat flux and CO₂ flux. Friction velocity is estimated from the mean wind speed with a flux-gradient relationship modified for the roughness sublayer. Since normalised standard deviations do not follow Monin-Obukhov similarity theory in the roughness sublayer, a simple classification scheme according to the scalar turbulence scale was used. This scheme is shown to produce sensible heat and CO₂ flux estimates that are well correlated with the measured values.     

Applicability of Temperature Discrete Equation to NMRF Boundary Layer Scheme in GRAPES Model

By deriving the discrete equation of the parameterized equation for the New Medium-Range Forecast (NMRF) boundary layer scheme in the GRAPES model, the adjusted discrete equation for temperature is obviously different from the original equation under the background of hydrostatic equilibrium and adiabatic hypothesis. In the present research, three discrete equations for temperature in the NMRF boundary layer scheme are applied, namely the original (hereafter NMRF), the adjustment (hereafter NMRF-gocp), and the one in the YSU boundary-layer scheme (hereafter NMRF-TZ). The results show that the deviations of height, temperature, U and V wind in the boundary layer in the NMRF-gocp and NMRF-TZ experiments are smaller than those in the NMRF experiment and the deviations in the NMRF-gocp experiment are the smallest. The deviations of humidity are complex for the different forecasting lead time in the three experiments. Moreover, there are obvious diurnal variations of deviations from these variables, where the diurnal variations of deviations from height and temperature are similar and those from U and V wind are also similar. However, the diurnal variation of humidity is relatively complicated. The root means square errors of 2m temperature (T2m) and 10m speed (V10m) from the three experiments show that the error of NMRF-gocp is the smallest and that of NMRF is the biggest. There is also a diurnal variation of T2m and V10m, where T2m has double peaks and V10m has only one peak. Comparison of the discrete equations between NMRF and NMRF-gocp experiments shows that the deviation of temperature is likely to be caused by the calculation of vertical eddy diffusive coefficients of heating, which also leads to the deviations of other elements.     

Aerodynamic Properties of Rough Surfaces with High Aspect-Ratio Roughness Elements: Effect of Aspect Ratio and Arrangements

We examine the effect of varying roughness-element aspect ratio on the mean velocity distributions of turbulent flow over arrays of rectangular-prism-shaped elements. Large-eddy simulations (LES) in conjunction with a sharp-interface immersed boundary method are used to simulate spatially-growing turbulent boundary layers over these rough surfaces. Arrays of aligned and staggered rectangular roughness elements with aspect ratio >1 are considered. First the temporally- and spatially-averaged velocity profiles are used to illustrate the aspect-ratio effects. For aligned prisms, the roughness length (\(z_\mathrm{o}\)) and the friction velocity (\(u_*\)) increase initially with an increase in the roughness-element aspect ratio, until the values reach a plateau at a particular aspect ratio. The exact value of this aspect ratio depends on the coverage density. Further increase in the aspect ratio changes neither \(z_\mathrm{o}\), \(u_*\) nor the bulk flow above the roughness elements. For the staggered cases, \(z_\mathrm{o}\) and \(u_*\) continue to increase for the surface coverage density and the aspect ratios investigated. To model the flow response to variations in roughness aspect ratio, we turn to a previously developed phenomenological volumetric sheltering model (Yang et al., in J Fluid Mech 789:127–165, 2016), which was intended for low to moderate aspect-ratio roughness elements. Here, we extend this model to account for high aspect-ratio roughness elements. We find that for aligned cases, the model predicts strong mutual sheltering among the roughness elements, while the effect is much weaker for staggered cases. The model-predicted \(z_\mathrm{o}\) and \(u_*\) agree well with the LES results. Results show that the model, which takes explicit account of the mutual sheltering effects, provides a rapid and reliable prediction method of roughness effects in turbulent boundary-layer flows over arrays of rectangular-prism roughness elements.