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

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

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

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

半干旱区不同下垫面近地层湍流通量特征分析

本文分析了2003～2005年半干旱区退化草地和农田下垫面近地层湍流通量日、季、年的变化特征, 探讨了不同年份的气候背景和下垫面土地覆被对地气交换过程的影响。结果表明: 半干旱区退化草地和农田近地层湍流通量具有明显的日变化、季变化周期; 地气间通量交换年际间的差异主要受当年的气候背景影响, 尤其是降水的影响; 同时还受到下垫面覆被的影响。潜热通量和感热通量的能量分配比率呈反位相变化, 且农田和退化草地的变化趋势相似, 在夏季潜热通量所占比例均超过感热通量; 两种下垫面的波恩比月均值变化趋势十分相似, 量级接近, 夏季低, 春、秋季高; 夏季退化草地和农田下垫面波恩比均小于或等于1。     

南京夏季城郊湍流统计特性及湍流通量对比分析

利用2010年南京夏季城市热岛三维结构试验中的湍流观测资料,对南京市委党校和南京信息工程大学观测场2个观测点的湍流统计量和湍流通量进行了对比分析。结果表明：城市平均风速大于郊区;城郊2地的湍流强度都呈现出Ix〉Iy〉Iz的趋势,水平方向上的湍强城郊相差不大,而垂直方向上,郊区的Iz要明显低于城市;城、郊的风速归一化标准差σu/u＊、σv/u＊、σw/u＊都符合MoninObukhov相似理论的＂1/3＂定律,城市的垂直风速归一化标准差大于郊区,而水平方向上的风速归一化标准差则明显小于其他下垫面;郊区风速归一化标准差与常熟农田的拟合结果较为相似,而城市风速归一化标准差与长白山森林的拟合结果更为接近;夏季城市以感热通量为主,而郊区湍流能量的输送以潜热通量为主,且城市的湍流热通量受太阳辐射的影响更大。     

基于WRF模式的海面湍流通量参数化方案的研究

中尺度数值模拟结果特别是高影响天气的精细预报对近地层动量和热量通量极为敏感,因此近地层湍流通量参数化方案一直是大气科学研究中一个十分重要的课题.以TOGA-COARE观测试验资料为基础,本文将湍流通量参数化方案模块从天气研究预报(WRF)模式中提取出来,与最新研发的湍流通量参数化方案(即LGLC方案)进行对比测试分析....     

湍流通量计算方法和误差的比较研究

近地层的风、温、湿梯度以及辐射和土壤热通量资料在过去的几十年间被广泛地用于计算陆-气感热、潜热通量.然而,不同的计算方法之间存在一些差异,随意的选用可能会造成计算结果存在较大的误差.文中利用绿洲系统能量与水分循环过程观测试验的2005年6月甘肃金塔绿洲中部观测资料,在分析了观测期间近地层微气象特征的基础上,对鲍恩比法、空气动力学法(包含两种不同计算方案)和变分法计算的绿洲农田下垫面的湍流通量进行了比较,并通过敏感性实验和理论分析探讨了不同方法间计算结果的差异和误差的来源.结果表明:不同方法计算的湍流通量特征存在明显差异.鲍恩比法从理论上满足能量平衡关系,但在昼夜交替和夜晚时出现了计算不稳定现象,因此在这些时段不适用.变分法避免了鲍恩比法计算夜间湍流通量的不稳定性,而且二者计算的白天湍流通量有较好的相关性;同时变分法由于加入了辐射和土壤热通量信息,明显改善了空气动力学法白天的能量平衡状况;变分法计算得到的潜热通量对辐射和微气象条件有较为合理的响应;对数据误差的敏感性试验表明,鲍恩比法和空气动力学法应用于绿洲下垫面上对观测误差较为敏感,而变分法相对比较稳定.进一步的分析表明,鲍恩比法在-2<β<0的范围内由于方法本身局限而适用性较差;而空气动力学方法计算结果的差异和误差来源于所应用的相似性函数.因而与鲍恩比法和空气动力学法计算的湍流通量值相比较,变分法的计算结果更加合理和稳定,在绿洲下垫面有较好的适用性.该文研究方法和结果可为这些方法的使用提供参考.     

台风“珍珠”登陆期间动量通量的多尺度分析

利用台风"珍珠"登陆前后的近地层湍流观测资料,分析了该台风经过观测场地前后地面气象要素的变化及其动量通量特征。结果表明,台风"珍珠"经过观测场地前后的近地层气象要素发生了急剧的变化,并且在台风前部存在强的中尺度对流系统,反映在风速能谱密度结构上,频率f在3×10-4~2×10-3Hz之间的中尺度信号对能谱的贡献比平稳天气形势下的能谱贡献大很多,尤其是顺风方向风速的能谱密度的峰值与湍流信号的峰值相当;动量通量分析结果表明,台风中心经过观测场地的前后三小时,近地层通量以向下输送的中尺度通量为主,湍流通量的贡献相对于中尺度通量较小,也是向下输送的;而在其他时段,近地层通量主要以向上输送的湍流通量为主,中尺度通量量值很小,可以忽略。     

珠穆朗玛峰绒布河谷近地层大气湍流及能量输送特征分析

利用中国科学院2005年珠穆朗玛峰地区科学考察期间 (4月2日至6月7日) 收集的大气观测资料,分析了珠峰绒布河谷近地层水平风速、温度、湍流强度、湍流通量日变化及地表能量平衡特征。通过分析得出近地层三维风速方差与稳定度的关系基本满足1/3次方规律;珠峰绒布河谷近地层大气水平风速、温度、动量通量、感热通量和潜热通量均存在明显的日变化;地表获得的能量很大一部分以感热形式散失掉了,潜热所占比重很小。另外,还发现绒布河谷地区地表能量通量各分量并不满足能量平衡方程Rn=Hs+Le+G。通过对地面加热场的分析发现珠峰地表白天是强热源,晚上转变为弱冷源。     

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

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

不同平均时间对LOPEX10资料涡动相关湍流通量计算结果影响的探讨

通过处理涡动相关系统观测的近地层湍流脉动量可以获取地—气间感热和潜热通量,然而选择不同平均时间对通量计算的结果有较大影响。采用黄土高原陆面过程野外观测试验(LOPEX10)期间获得的涡动相关系统观测资料,分析了不同平均时间对湍流通量计算的影响,并采用雷诺平均和分解方法推导了平均时间引起的通量差值的数学表达式(Flux Compensation,FC)。结果表明:(1)FC公式可以说明采用不同平均时间数据之间的关系,也可以直接计算低频涡旋对湍流通量的贡献。FC公式计算的结果与直接计算的不同平均时间通量计算之差的相关系数在0.95以上,并可以确定计算湍流通量的最佳平均时间。(2)通过采用Ogive函数确定了计算LOPEX10期间通量的最佳平均时间长度为30min,印证了利用雷诺平均和分解方法计算湍流通量补偿的准确性。(3)通过进一步的数学变换,证明了平均时间对湍流通量计算的影响直接与湍流低频变化相关,FC公式可以用来确定涡动相关观测数据的最佳平均时间,并且在获得较高时间分辨率的湍流通量数据的同时,补偿因平均时间过短而遗漏的低频信息。     

六种近地层湍流动量输送系数计算方案对比分析

选取30多年来近地层湍流通量研究中具有代表性的六种参数化方案, 应用GAME/Tibet试验中那曲通量观测站的实测资料, 对比分析了各方案计算所得的湍流动量输送系数 (C_M) 之间的差异。结果表明:六种参数化方案计算得到的湍流动量输送系数之间存在较大差异。对于那曲观测站稀疏短草下垫面而言, 稳定条件下当理查孙数小于0.1时, 除Businger71方案存在显著低估以外, 其他各方案均能较好估算湍流动量输送系数; 不稳定条件下, Dyer74方案对湍流动量输送系数的估算效果最好, 其次为Wang02, Launiainen95和Louis82方案, Businger71方案误差较大。     

THE FORMATION AND THE DEVELOPMENT OF CAMS TWO-DIMENSIONAL DYNAMICAL-PHYSICAL CLIMATE MODEL

In the paper,we have developed a 2-D physical-dynamical coupled climate model.Some sensitive experiments have been done by use of this model.First of all,we have studied the effects of different results by two radiational calculation schemes on circulation variation.The calculated results have shown that the different radiation parameterization schemes give different results,therefore the variational effects of wind,temperature,and humidity field are presented on the medium-range circulation variation.Besides,we have also studied the role of the meridional eddy momentum fluxes in formation of the monsoon over East Asia.The results of study have shown that on the average,the roles of meridional eddy flux of momentum formative processes of monsoon which is added to momentum equations by using the scheme of moist process parameterization with plateau are manifest.     

RADIATION BALANCE AND TURBULENT FLUX CHARACTERISTICS OVER MIZUHO STATION IN ANTARCTICA

In this paper,by using measurements in micrometeorology and radiation balance in the surface layer(Yamanouchiet al.,1981;Wada et al.,1981;Ishikawa et al.,1982;Ohata et al.,1983),we have analyzed the diurnal and annual varia-tion characteristics of radiation balance in spring,summer,autumn and winter,and calculated the momentum flux andsensible heat flux with the aerodynamic method and profile gradient alternate method in different seasons.We have alsoobtained the diurnal variation characters of the latent heat flux from the equation for energy balance.The results fromthe calculation are compared.Finally,the relationship between the turbulent heat and momentum exchange coefficientand the diurnal variation of the Richardson number is discussed.     

Energy- and Flux-Budget Turbulence Closure Model for Stably Stratified Flows. Part II: The Role of Internal Gravity Waves

We advance our prior energy- and flux-budget (EFB) turbulence closure model for stably stratified atmospheric flow and extend it to account for an additional vertical flux of momentum and additional productions of turbulent kinetic energy (TKE), turbulent potential energy (TPE) and turbulent flux of potential temperature due to large-scale internal gravity waves (IGW). For the stationary, homogeneous regime, the first version of the EFB model disregarding large-scale IGW yielded universal dependencies of the flux Richardson number, turbulent Prandtl number, energy ratios, and normalised vertical fluxes of momentum and heat on the gradient Richardson number, Ri. Due to the large-scale IGW, these dependencies lose their universality. The maximal value of the flux Richardson number (universal constant ≈0.2–0.25 in the no-IGW regime) becomes strongly variable. In the vertically homogeneous stratification, it increases with increasing wave energy and can even exceed 1. For heterogeneous stratification, when internal gravity waves propagate towards stronger stratification, the maximal flux Richardson number decreases with increasing wave energy, reaches zero and then becomes negative. In other words, the vertical flux of potential temperature becomes counter-gradient. Internal gravity waves also reduce the anisotropy of turbulence: in contrast to the mean wind shear, which generates only horizontal TKE, internal gravity waves generate both horizontal and vertical TKE. Internal gravity waves also increase the share of TPE in the turbulent total energy (TTE = TKE + TPE). A well-known effect of internal gravity waves is their direct contribution to the vertical transport of momentum. Depending on the direction (downward or upward), internal gravity waves either strengthen or weaken the total vertical flux of momentum. Predictions from the proposed model are consistent with available data from atmospheric and laboratory experiments, direct numerical simulations and large-eddy simulations.     

WRF模式中微物理和积云参数化方案的对比试验

为了研究微物理参数化方案对珠江三角洲(简称珠三角)降水模拟的影响,利用WRF中尺度数值预报模式,在3 km模式分辨率下,在微物理方案为WSM6方案条件下,选用KF、BMJ、GD以及G3等四种积云参数化方案对2010年5月14日广东珠三角地区的一次暴雨过程进行了模拟试验。结果显示,KF方案对于降水带和降水量的模拟与实况较为一致。在积云参数化方案为KF条件下,分别选用Kessler、Lin et al、WSM 3、WSM5、Ferrier(New eta)和WSM6等6种微物理方案再次对这次暴雨过程进行模拟试验,模拟结果的对比分析表明:选用Lin et al微物理方案时,模式较好地模拟出了强降水雨带的位置和降水强度;而其他5种参数方案的模拟效果均不好,降水量明显偏小,雨带位置偏差较大;同时对低空急流、K指数和上升速度等物理量分析可知,Lin et al方案能较好地模拟出降水实况。     

湍流通量参数化方案的非迭代方法研究

基于Högström (1996) 和Beljaars et al.(1991) 的研究工作, 沿用Louis et al.(1982) 和Launiainen (1995) 的思路, 本文采用多元回归分析方法, 研发了一种采用非迭代方法的湍流通量参数化方案。该方案直接用整体理查森数、 空气动力学粗糙度长度和热力学粗糙度长度对稳定度参数进行参数化, 从而避免了通过循环迭代计算Monin－Obukhov长度。该方案不仅有效地节省了CPU计算时间, 而且其计算结果与迭代方案 (BHH方案) 的计算结果非常接近。     

ECMWF预报模式在不同地形方案下的系统误差检验

Wallace等的研究表明,由于引进了包络地形方案,欧洲中期天气预报中心(ECMWF)的中期天气预报有所改善。本文是上述工作的继续。所用的资料为原作者的数值试验结果。检验的对象包括纬向平均风场、温度场和经圈环流,及不同波数域中的涡旋和经圈环流对角动量的平均经向输送。结果表明,在冬季,包络地形方案对低纬地区的温度场和风场的预报略有改进,对总涡旋动量通量的模拟也较好。然而对天气尺度系统的预报和对热带对流及积云摩擦的参数化却有待更进一步改善。     

A Comparison Study of the Climatological Air-Sea Heat Fluxes Estimated by Different Computational Schemes of Bulk Formula

In this paper the different results of estimated air-sea heat fluxes by different computational schemes of bulk formula are researched. The varying results may bring another source of the uncertainties of flux climatology.In addition to the classical scheme and the sampling scheme, a revised sampling scheme named semi-sampling scheme is proposed in this paper. In this revised version, the monthly means of and are cal?culated from each simultaneous measurement as that in the sampling method. But the transfer coefficients in this version are defined based on monthly means of and instead of those in each measurement in older to use the files of monthly mean variables and in COADS (Woodruff et al. 1987).The results of a comparison study show that the covariance between and and that between and are less than 2 W / m2 and 5 W / m2 respectively in most areas, which are about one order of magni?tude less than the standard deviation of the fluxes. Therefore the difference between two schemes is not statistically significant on average. However, in the regions off coast of Asian and North American continents, the contribution of covariance between and can be as large as 15 W / m2, which is about one third to a half of the standard deviation. Since there are files of monthly mean and for each year in the COADS, we pro?pose to use these files directly in calculation of climatology. It is not necessary to use every individual report.The flux climatology calculated directly from the long-term mean values of basic meteorological elements shows an significant difference from that calculated from the mean of each year’s flux, particularly for the latent heat flux which is on the order of 10 to 15 W / m2, about 3-4 times of the difference between semi-sampling and classical schemes and about one third of the standard deviation. Therefore it indicates that the time interval of months is per?haps a better one in the calculation of flux climatology. It is not recommended to use the long-term means of the me?teorological variables.     

Implementation of a Stable PBL Turbulence Parameterization for the Mesoscale Model MM5: Nocturnal Flow in Complex Terrain

The difficulties associated with the parameterization of turbulence in the stable nocturnal planetary boundary layer (PBL) have been a great challenge for the nighttime predictions from mesoscale meteorological models such as MM5. As such, there is a general consensus on the need for better stable boundary-layer parameterizations. To this end, two new turbulence parameterizations based on the measurements of the Vertical Transport and Mixing (VTMX) field campaign were implemented and evaluated in MM5. A unique aspect of this parameterization is the use of a stability-dependent turbulent Prandtl number that allows momentum to be transported by the internal waves, while heat diffusion is impeded by the stratification. This improvement alleviates the problem of over-prediction of heat diffusion under stable conditions, which is a characteristic of conventional atmospheric boundary-layer schemes, such as the Medium Range Forecast (MRF) and Blackadar schemes employed in MM5. The predictions made with the new PBL scheme for the complex terrain airshed of Salt Lake City were compared with those made with a default scheme of MM5, and with observations made during the VTMX campaign. The new schemes showed an improvement in predictions, particularly for the nocturnal near-surface temperature. Surface wind predictions also improved slightly, but not to the extent of temperature predictions. The default MRF scheme showed a significantly higher surface temperature than observed, which could be attributed to the enhanced vertical heat exchange brought about by its turbulence parameterization. The modified parameterizations reduced the surface sensible heat flux, thus enhancing the strength of the near-surface inversion and lowering the temperature towards the observed values.     

Comparison of atmospheric turbulence characteristics and turbulent fluxes from two urban sites in Essen, Germany

From September 2006 to September 2007, the intersite variability of turbulence characteristics and turbulent heat fluxes was analysed at two urban stations in Essen, Germany. One site was situated within an urban residential setting while the other was located at the border of an urban park and suburban/urban residential housing. Therefore, the surroundings at both sites contributing to surface–atmosphere exchange differed in terms of surface cover and surface morphology. During the 1-year measurement period, 19% of data were characterised by stable atmospheric stratification. Since observations of urban turbulence characteristics under stable stratification are scarce, so far, this work adds additional input to this discussion. Turbulence characteristics, i.e. normalised standard deviations of wind components, were in agreement to empirical fits from other urban observations under both instable and stable atmospheric stratification. However, differences in magnitude of turbulence characteristics between sites were observable. Comparison of turbulent heat fluxes indicated typical urban features in the site located in the urban setting with increased surface heating and higher surface heat fluxes by about 30%. Also the temporal evolution of heat fluxes on the diurnal course was affected. Differences in momentum flux were of minor magnitude with about 6% variation on average between sites. Findings indicate that multiple urban flux measurements within one city may be characterised by general similarities in terms of turbulent characteristics but are still significantly influenced by differences in the surface cover of the flux footprint.