农田植被层上方湍流通量输送特征分析

根据宵田植被层上方的温湿风梯度观测资料,采用基于莫宁－奥布霍夫相似理论的廓线梯度法计算了农田植被层上方的湍流通量,进而确定出拖曳系数、斯坦顿数和道尔顿数。文中检验了用于农田植被层上方湍流通量计算的普适函数的适用性,讨论了湍流通量以及拖曳系数等随稳定度参数、粗糙度高、平均风速等因子的变化规律。结果表明,农田植被层上方的湍流通量随层结构稳定度变化而变化;拖曳系数、斯坦顿数和道尔顿数也不是常数,而是层结     

北京北郊冬季大风过程湍流通量演变特征的分析研究

利用中国科学院大气物理研究所325 m气象观测塔1993年12月～1994年1月大气边界层实验资料, 计算分析了大风过境过程中47 m和120 m高度湍流通量演变特征及其影响因子, 以及与风速、 稳定度等参数的关系。结果表明: 大风过程对近地面层的物质能量输送有着重要影响, 大风之前出现短时间动量上传和热量下传; 大风过程中的湍流通量数值明显高于过境后, 水平方向湍流通量数值和能量增加幅度大于垂直方向; 当风速大于临界值5 m/s时, 湍流通量与风速、 湍流动能的相关迅速增大; 湍流谱特征表现为湍流能量的低频部分增加、 湍流谱曲线变宽; 大风能强烈影响近地面层的能量收支。     

论近地面层湍流通量观测的一些问题

本文回顾并总结了目前在近地面层内观测湍流通量的方法。指出了各种观测方法的适用条件。比较了它们的测量误差和它们的优缺点，并讨论了各种测量方法存在的问题。该文的讨论对于选择近地面层湍流通量观测方案、提高观测精度、评价观测结果并对结果进行修正都有一定参考价值。     

地表粗糙度非均匀性对模式湍流通量计算的影响

理论分析和数值试验结果表明,地表参量的不均匀性对网格区地表湍流通量的计算有重要影响。近中性大气层结条件下,次网格粗糙度长度的变差系数、网格平均粗糙度及参考高度的选取是影响网格湍流通量计算的主要因素,其中次网格粗糙度长度的变差系数对计算偏差起主要决定作用。实际计算表明,某些特定地区(如植被气候过渡带)粗糙度的地表非均匀性引起的计算相对误差可达40%以上,选取特定的参考高度能改善高网格湍流通量计算的效果。非中性大气层结条件下,由地表粗糙度不均匀性所致的平均风速、位温梯度以及近地层大气稳定度的次网格分布都对网格湍流通量(感热通量)计算产生影响,比较而言,相对误差大小对大气稳定度的次网格分布最为敏感。所以,在目前的数值模式中有必要进一步对湍流通量计算过程中由于地表不均匀性产生的计算偏差加以考虑。     

青藏高原东南部昌都地区近地层湍流输送的观测研究

利用我国1998年第2次青藏高原气象科学试验 (TIPEX) 昌都的湍流观测资料, 初步分析了三维风速、温度和湿度湍流谱、归一化方差与稳定度的关系、湍流通量和整体输送系数.结果表明, 白天标准化的湍流谱大多数满足2/3次方定律.风速方差与稳定度的关系在不稳定条件下基本服从1/3次方定律, 但水平风速方差比垂直速度方差与稳定度的关系离散.中性层结下的A、B值比平原地区大, 与起伏地形的值接近, C值与平原地区和起伏地区的差异不大.温度和湿度方差与不稳定层结的关系基本满足1/3次方定律, 而在稳定层结下的关系均无明显规律.湍流通量在干湿期转换中有突变现象, 干期的动量和热量通量均大于湿期, 热源强度以感热为主, 占80%.湿期热源强度明显增大, 其中以潜热为主, 占热源强度的60%.整体输送系数与稳定度的关系明显, 拟合得到的中性层结条件下的C_DN和C_HN分别为4.40×10-3和4.70×10-3, 其结果大于高原中、西部的同期观测值.     

近地面层中湍流热通量和长波辐射通量相互作用的数值研究

本文详细讨论了在稳定层结条件下,近地面层内有各种典型的气象要素分布及下垫面状况时,长波辐射通量与湍流热通量之间的相互关系问题.结果表明:长波辐射通量对近地面层的热力结构有相当大的影响,它常常减弱了常值湍流热通量层的厚度.在弱风和层结稳定度较大的情况下,“常值湍流热通量层”将很薄,有时甚至难以存在.     

大理近地层山谷盆地湖陆风及湍流特征分析

研究湖陆风特征不仅能够为提高天气气候的预测能力奠定基础,而且对风能资源的开发利用等具有重要的实用意义。利用大理国家气候观象台近地面通量观测系统的2007年3月－2008年5月资料,采用涡动相关法等分析了大理近地层中湖陆风、峡谷风特征及形成原因和影响因素。结果表明：大理地区白天以东风和东南风为主,夜间以西风和西南风为主。进一步对湍流和湍流通量特征分析发现,大理地区白天不稳定层结多于夜间;湍流强度白天强于夜间,并且随着风速的增大而减小;湍流通量具有明显的日变化特征,热量交换形式以潜热为主。     

那曲高寒草地上四种地表通量计算方法的对比

利用中国科学院那曲高寒气候环境观测研究站2013年9月至2014年8月自动气象站(AWS)和涡动相关系统(EC)的观测资料,基于空气动力学法、地表能量平衡组合法、总体输送法以及涡动相关法,计算了高寒草地下垫面的湍流通量,并对不同方法计算结果间的一致性和差异性进行了分析。结果表明,不同方法计算的湍流通量特征具有明显的差异。地表能量平衡组合法满足能量平衡关系,但在早晨和傍晚层结转换期间,计算的湍流通量出现异常不稳定值;空气动力学法计算的湍流通量在整个观测期与涡动相关法计算的湍流通量相关性最高,但在大气稳定度参数接近0,计算结果不稳定;总体输送法计算的通量数据在地气温差为负值时发散明显,但该方法原理简单,适合在只有常规观测项目的业务气象站或在气象观测项目不全的野外台站使用。空气动力学法和地表能量平衡组合法与涡动相关法的湍流通量的平均偏差相对较小,而总体输送法平均偏差相对较大。研究方法和结果除了为这些方法的使用提供参考外,也为建立长时间通量序列提供了一个较为合理的依据,有助于深入了解高原地气相互作用。     

组合法确定近地面层湍流通量和通用函数

本文发展了确定近地面层湍流通量的组合法,得到了它们的一般表达式。并用这种方法确定近地面层MO相似性函数。利用ITCE的观测资料作为实例计算,结果表明这种方法是可行的。     

稳定大气近地面层的风和温度廓线

本文不做风和温度廓线相似性假设,使用最小二乘法分析位于粗糙、不均匀地形上北京气象专用塔在稳定条件下的平均风速和温度测量资料。结果表明,在近地面层,从弱稳定到强逆温的较宽的稳定层结范围内,对数加线性规律都是适用的;廓线参数β和β不是常数,它们的数值随大气稳定度而变化。文中给出的经验公式表明:廓线参数随R_i的增大而减少,随L的增大而增大。在临界理查逊数的两侧,廓线参数随R_i的变化遵循不同的规律。文中还表明,廓线相似性假设仅对大气层结稳定度较弱时才是真实的;当稳定度增大时,廓线变成不相似,P值随R_i的增大而增大。     

纳木错(湖)地区湍流数据质量控制和湍流通量变化特征

利用中国科学院纳木错多圈层综合观测研究站2009年全年的大气湍流观测资料,应用Foot-print模型分析了青藏高原非均匀下垫面湍流观测数据的数据质量、质量评价及不同下垫面对湍流通量的贡献。结果表明:纳木错(湖)地区因不同土地利用类型的差别,导致地表通量分布不均匀,草地对地表通量的贡献最大;对不同大气层结状态,观测站周围200m范围内的地表通量贡献各不相同,上风向通量贡献源区较大,湍流发展较充分。在不稳定状态和中性状态下,纳木错地区地表通量数据质量较高,即白天观测的通量数据质量较高;在稳定状态下数据质量较低,即夜间的通量数据质量较差;纳木错地区的湍流通量受湖陆风和大气稳定性影响较大。     

Convective scaling of the average dissipation rate of temperature variance in the atmospheric surface layer

The flux of sensible heat from the land surface is related to the average rate of dissipation of temperature fluctuations in the atmospheric surface layer through the temperature variance budget equation. In many cases it is desirable to estimate the heat flux from measurement or inference of the dissipation rate. Here we study how the dissipation rate scales with atmospheric stability, using three inertial range methods to calculate the dissipation rate: power spectra, second order structure functions, and third order structure functions. Experimental data are analyzed from a pair of field experiments, during which turbulent fluctuations of velocity and temperature were measured over a broad range of neutral and unstable atmospheric flows. It is shown that the temperature dissipation rate scales with a single convective power law continuously from near-neutral to strongly unstable stratification. The dissipation scaling is found to nearly match production in the near-neutral region, but to be consistently lower than production in the more convective regimes. The convective scaling is shown to offer a simplified means of computing sensible heat flux from the dissipation rate of temperature variance.Also at Johns Hopkins University, Baltimore, MarylandAlso at Los Alamos National Laboratory, Los Alamos, New Mexico.     

内蒙古奈曼流动沙丘下垫面湍流输送特征初步研究

分析了2000年7月21日～8月10日在内蒙古奈曼流动沙丘下垫面取得的近地面层湍流及辐射观测资料,讨论了无量纲湍流方差与稳定度参数z/L的关系,发现无量纲速度分量方差及无量纲温度、湿度脉动方差在不稳定层结下,均满足莫宁-奥布霍夫相似理论;同时还讨论近地面层能量的收支.发现流动沙丘下垫面感热通量与地表热通量最大值分别为170 W m-2及100 W m-2;潜热通量通常小于50 W m-2;流动沙丘净辐射最大值为400 W m-2左右.(Hs+Hl),与(Rn-G)的比值在晴天白天平均值为0.78左右,流动沙丘下垫面近地面层能量不平衡现象存在,其原因有待进一步研究.     

稻麦轮作农田近地层湍流通量计算方案对比研究

稻麦轮作农田是我国典型农田类型,其模拟效果对我国农田气候模拟具有重要参考价值。气候中尺度模拟结果对近地层通量极为敏感,选择合适的通量计算方案对模拟效果至关重要。因此,对比分析稻麦轮作农田下不同的通量计算方案具有重要意义。本文选取了8种具有代表性的近地层湍流通量计算方案,采用寿县国家气候观象台实测资料对比分析了各方案在稻麦轮作农田的计算特征和差异。结果表明,在不同稳定度和不同风速情况下,各方案的误差特征各异。本文基于归一化标准差综合评价了各方案的准确度,总体而言,所有方案的动量通量总体平均归一化平均差为0.536,其中SS14(Sharan和Srivastava,2014年)方案最大为0.575,SS20(Sharan和Srivastava等,2020年)方案最小为0.517;所有方案的感热通量总体平均归一化标准差为0.638,其中GLGS20(Gryanik等,2020年)方案最大为0.871,SS14方案最小为0.476。此外,本研究还给出了稻麦轮作农田不同稳定层结和不同风速情况下各通量计算方案的误差特征。本文的研究结论,可为准确计算近地层湍流通量提供支撑。     

Aircraft Observations of Sea-Surface Turbulent Fluxes Near the California Coast

Aircraft turbulence data from the Autonomous Ocean Sampling Network project were analyzed and compared to the Coupled Ocean–Atmosphere Response Experiment (COARE) bulk parametrization of turbulent fluxes in an ocean area near the coast of California characterized by complex atmospheric flow. Turbulent fluxes measured at about 35 m above the sea surface using the eddy-correlation method were lower than bulk estimates under unstable and stable atmospheric stratification for all but light winds. Neutral turbulent transfer coefficients were used in this comparison because they remove the effects of mean atmospheric conditions and atmospheric stability. Spectral analysis suggested that kilometre-scale longitudinal rolls affect significantly turbulence measurements even near the sea surface, depending on sampling direction. Cross-wind sampling tended to capture all the available turbulent energy. Vertical soundings showed low boundary-layer depths and high flux divergence near the sea surface in the case of sensible heat flux but minimal flux divergence for the momentum flux. Cross-wind sampling and flux divergence were found to explain most of the observed discrepancies between the measured and bulk flux estimates. At low wind speeds the drag coefficient determined with eddy correlation and an inertial dissipation method after corrections were applied still showed high values compared to bulk estimates. This discrepancy correlated with the dominance of sea swell, which was a usually observed condition under low wind speeds. Under stable atmospheric conditions measured sensible heat fluxes, which usually have low values over the ocean, were possibly affected by measurement errors and deviated significantly from bulk estimates.     

Observing and Simulating Wind-Turbine Wakes During the Evening Transition

Wind-turbine-wake evolution during the evening transition introduces variability to wind-farm power production at a time of day typically characterized by high electricity demand. During the evening transition, the atmosphere evolves from an unstable to a stable regime, and vertical stratification of the wind profile develops as the residual planetary boundary layer decouples from the surface layer. The evolution of wind-turbine wakes during the evening transition is examined from two perspectives: wake observations from single turbines, and simulations of multiple turbine wakes using the mesoscale Weather Research and Forecasting (WRF) model. Throughout the evening transition, the wake’s wind-speed deficit and turbulence enhancement are confined within the rotor layer when the atmospheric stability changes from unstable to stable. The height variations of maximum upwind-downwind differences of wind speed and turbulence intensity gradually decrease during the evening transition. After verifying the WRF-model-simulated upwind wind speed, wind direction and turbulent kinetic energy profiles with observations, the wind-farm-scale wake evolution during the evening transition is investigated using the WRF-model wind-farm parametrization scheme. As the evening progresses, due to the presence of the wind farm, the modelled hub-height wind-speed deficit monotonically increases, the relative turbulence enhancement at hub height grows by 50%, and the downwind surface sensible heat flux increases, reducing surface cooling. Overall, the intensifying wakes from upwind turbines respond to the evolving atmospheric boundary layer during the evening transition, and undermine the power production of downwind turbines in the evening.     

A multi-limit formulation for the equilibrium depth of a stably stratified boundary layer

Currently no expression for the equilibrium depth of the turbulent stably-stratified boundary layer is available that accounts for the combined effects of rotation, surface buoyancy flux and static stability in the free flow. Various expressions proposed to date are reviewed in the light of what is meant by the stable boundary layer. Two major definitions are thoroughly discussed. The first emphasises turbulence and specifies the boundary layer as a continuously and vigorously turbulent layer adjacent to the surface. The second specifies the boundary layer in terms of the mean velocity profile, e.g. by the proximity of the actual velocity to the geostrophic velocity. It is shown that the expressions based on the second definition are relevant to the Ekman layer and portray the depth of the turbulence in the intermediate regimes, when the effects of static stability and rotation essentially interfere. Limiting asymptotic regimes dominated by either stratification or rotation are examined using the energy considerations. As a result, a simple equation for the depth of the equilibrium stable boundary layer is developed. It is valid throughout the range of stability conditions and remains in force in the limits of a perfectly neutral layer subjected to rotation and a rotation-free boundary layer dominated by surface buoyancy flux or stable density stratification at its outer edge. Dimensionless coefficients are estimated using data from observations and large-eddy simulations. Well-known and widely used formulae proposed earlier by Zilitinkevich and by Pollard, Rhines and Thompson are shown to be characteristic of the above interference regimes, when the effects of rotation and static stability (due to either surface buoyancy flux, or stratification at the outer edge of the boundary layer) are roughly equally important.     

Determination of turbulent momentum and heat fluxes by spectral methods

Methods are studied which permit one to evaluate turbulent fluxes from the results of spectral measurements in turbulent laboratory flows and an unstable atmospheric surface layer. The well known dissipation method of flux measurements, which uses spectral data related to the inertial range, is reanalyzed. New theoretical ideas and the latest experimental data are used to specify this method in cases of moderately and very strongly unstable thermal stratifications.Moreover, it is also explained how to estimate momentum and heat fluxes from data in the low frequency parts of the velocity and temperature spectra in the low frequency ranges beyond the lower limit of the inertial range. This permits one to estimate fluxes using rather simple and cheap instruments (e.g., Pilot-tubes and thermocouples in laboratory flows or cup anemometers and crude resistance thermometers in meteorological studies). The equations for flux determination are based in such cases on the recent models by Kader (1987, 1988) and Kader and Yaglom (1990, 1991) of spectral shapes at mesoscale wave numbers; these models agree quite satisfactorily with many (though not all) data of direct spectral measurements. It is shown that estimated momentum and heat fluxes in the laboratory and in an unstably stratified atmospheric surface layer obtained by the method suggested in this paper agree satisfactorily with direct flux measurements.     

天津市秋季典型环境污染过程个例分析

选择天津市秋季典型PM10污染过程2010年10月3—12日环境空气质量监测资料和常规气象资料、探空资料及NCEP资料,研究大气环境天气背景场、大气层结稳定度的特征及其对污染过程的影响。结果表明：高低空环流背景场与污染过程密切相关。在污染上升阶段,层结稳定度迅速增加,500 hPa高空处于槽前,地面在华北地形槽中,高低空风速辐合;在污染峰值阶段,层结稳定,逆温层加强,环流场稳定少变,地面风力微弱;在污染下降阶段,层结稳定程度骤降,地面冷锋和高空槽过境,降水出现,高低空偏北风增大。同时,PM10污染过程与多项层结稳定度参数显著相关,与对流凝结高度单相关系数为0.84,因此,高低空环流背景场的配置和层结稳定度变化是PM10污染出现的主要原因。     

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

分析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.