Impact of Different Cloud Microphysics Parameterization Schemes on Typhoon Intensity and Structure

The impact of different cloud microphysics parameterization schemes on the intensity and structure of the Super-strong Typhoon Rammasun (1409) in 2014 is investigated using the Weather Research and Forecasting model version 3.4 with eight cloud microphysics parameterization schemes. Results indicate that the uncertainty of cloud microphysics schemes results in typhoon forecast uncertainties, which increase with forecast time. Typhoon forecast uncertainty primarily affects intensity predictions, with significant differences in predicted typhoon intensity using the various cloud microphysics schemes. Typhoon forecast uncertainty also affects the predicted typhoon structure. Greater typhoon intensity is accompanied by smaller vortex width, tighter vortex structure, stronger wind in the middle and lower troposphere, greater height of the strong wind region, smaller thickness of the eyewall and the outward extension of the eyewall, and a warmer warm core at upper levels of the eye. The differences among the various cloud microphysics schemes lead to the different amounts and distributions of water vapor and hydrometeors in clouds. Different hydrometeors have different vertical distributions. In the radial direction, the maxima for the various hydrometeors forecast by a single cloud microphysics scheme are collocated with each other and with the center of maximum precipitation. When the hydrometeor concentration is high and hydrometeors exist at lower altitudes, more precipitation often occurs. Both the vertical and horizontal winds are the strongest at the location of maximum precipitation. Results also indicate that typhoon intensities forecast by cloud microphysics schemes containing graupel processes are noticeably greater than those forecast by schemes without graupel processes. Among the eight cloud microphysics schemes investigated, typhoon intensity forecasts using the WRF Single-Moment 6-class and Thompson schemes are the most accurate.     

Uncertainty in Tropical Cyclone Intensity Predictions due to Uncertainty in Initial Conditions

Focusing on the role of initial condition uncertainty,we use WRF initial perturbation ensemble forecasts to investigate the uncertainty in intensity forecasts of Tropical Cyclone(TC)Rammasun(1409),which is the strongest TC to have made landfall in China during the past 50 years.Forecast results indicate that initial condition uncertainty leads to TC forecast uncertainty,particularly for TC intensity.This uncertainty increases with forecast time,with a more rapid and significant increase after 24 h.The predicted TC develops slowly before 24 h,and at this stage the TC in the member forecasting the strongest final TC is not the strongest among all members.However,after 24 h,the TC in this member strengthens much more than that the TC in other members.The variations in convective instability,precipitation,surface upward heat flux,and surface upward water vapor flux show similar characteristics to the variation in TC intensity,and there is a strong correlation between TC intensity and both the surface upward heat flux and the surface upward water vapor flux.The initial condition differences that result in the maximum intensity difference are smaller than the errors in the analysis system.Differences in initial humidity,and to a lesser extent initial temperature differences,at the surface and at lower heights are the key factors leading to differences in the forecasted TC intensity.These differences in initial humidity and temperature relate to both the overall values and distribution of these parameters.     

WRF中三种边界层参数化方案对新疆2.28大风过程模拟的对比分析

许多研究调查了模式预报对边界层方案的敏感性,但是这些研究基本上针对的是热力驱动的混合边界层。对于动力驱动的边界层,不同边界层方案的不同性能以及所带来的不同边界层气象要素的预报还不清楚。运用WRF3.4.1中三种边界参数化方案（YSU、MYJ、ACM2）对新疆2.28大风过程进行数值模拟分析,结果显示：三种边界参数化方案基本能模拟出发生大风的区域及大风过程中10m风速、2m温度和比湿的变化趋势;三种方案模拟的边界层内大气的温度、湿度出现差异与它们对边界层顶的夹卷过程、边界层内垂直混合的处理有关;YSU方案的模拟结果使得更多的高空动量下传,同时更多的有效位能转化为动能,MYJ方案模拟的20m/s的风场区域更大,受地形影响更明显,边界层内湍流更强。     

WRF模式边界层参数化方案对川渝盆地西南涡降水模拟的影响

应用WRF v4.0模式五种边界层参数化方案（YSU、MYJ、MYNN2、ACM2和SH）,对2016 年汛期（5～9月）在川渝盆地东部造成暴雨的所有西南涡过程进行了数值模拟,检验评估了它们对各量级降水的预报能力,并基于加密的L波段秒级探空资料对比分析了模拟与实况边界层结构的差异,结合各方案对湍流运动的算法特点探讨了其差异的原因,最后对ACM2方案进行了湍流强度调整,由此改善其对于川渝盆地边界层与西南涡降水的模拟能力。结果表明:ACM2和YSU方案TS评分表现较好,相对其它方案ACM2空报较少,这种可以根据周围环境的稳定性切换局地或非局地算法的方案更适合于盆地西南涡降水模拟,但边界层方案对西南涡降水的空报都较普遍,尤以大量级降水更明显;精细的探空资料进一步表明,所有方案模拟的白天边界层高度都偏高,湍流混合强度都偏强。通过参数调整而降低混合强度的ACM2方案,模拟的边界层温湿结构则更符合实际观测,其边界层下部温度更低、湿度更高,减少了大量级降水的空报,使盆地西南涡降水模拟有一定改善;边界层参数化方案对西南涡模拟的差别主要体现为不同的西南涡位置与降水强度,但归根到底都源于方案的局地或非局地特性、不同的混合强度这两方面原因。因此,根据不同特定区域下垫面环境与气候状况合理选择方案的特性和混合强度是准确模拟边界层结构及其降水过程的关键。     

北京地区夏季边界层结构日变化的高分辨模拟对比

使用WRF中尺度数值模式, 分别选用两种不同的边界层参数化方案 (MYJ, YSU) 和3种陆面参数化方案 (SLAB, Noah, RUC), 对2004年7月1日08:00—7月4日20:00 (北京时) 北京地区夏季边界层结构进行1 km的高分辨模拟。对比分析了近地面层风场、温度场以及边界层的日变化特征, 结果发现:WRF模式基本模拟出了北京夏季边界层的日变化特征; 在边界层方案中, MYJ方案描述的边界层结构较YSU方案合理; Noah陆面模式较好地反映了城市的热岛效应; 无降水时, 风速及边界层高度对于陆面过程不敏感, 而降水发生后, 陆面过程对于边界层结构的影响增大; 各方案模拟的城区风速明显偏大, 这是因为没有充分考虑城市建筑物的阻力作用。     

INFLUENCE OF MODEL HORIZONTAL RESOLUTION ON THE INTENSITY AND STRUCTURE OF RAMMASUN

We use the WRF (V3.4) model as the experimental model and select three horizontal resolutions of 15, 9, and 3km to research the influence of the model’s horizontal resolution on the intensity and structure of the super-strong typhoon Rammasun (1409) in 2014. The results indicate that the horizontal resolution has a very large impact on the intensity and structure of Rammasun. The Rammasun intensity increases as the horizontal resolution increases. When the horizontal resolution increases from 9km to 3km, the enhancement of intensity is more obvious, but the strongest intensity simulated by 3km horizontal resolution is still weaker than the observed strongest intensity. Along with the increase of horizontal resolution, the horizontal scale of the Rammasun vortex decreases, and the vortex gradually contracts toward its center. The vortex structure changes from loose to compact and deep. The maximum wind radius, thickness of the eye wall, and outward inclination of the eye wall with height decrease, and the wind in the inner core region, updraft motion along the eye wall, and strength of the warm core become stronger. Additionally, the pressure gradient and temperature gradient of the eye wall region increase, and the vortex intensity becomes stronger. When the horizontal resolution increases from 9km to 3km, the change in the Rammasun structure is much larger than the change when the horizontal resolution increases from 15km to 9km. When the model does not employ the method of convection parameterization, the Rammasun intensity simulated with 3km horizontal resolution is slightly weaker than the intensity simulated with 3km horizontal resolution when the Kain–Fritsch (KF) convection parameterization scheme is adopted, while the intensity simulated with 9km horizontal resolution is much weaker than the intensity simulated with 9km horizontal resolution when the KF scheme is adopted. The influence of the horizontal resolution on the intensity and structure of Rammasun is larger than the influence when the KF scheme is adopted.     

古尔班通古特沙漠大气边界层参数化方案的模拟评估[1]

为准确描述我国最大的固定/半固定沙漠-古尔班通古特沙漠区域的大气边界层结构，本文利用该沙漠腹地2017年的梯度铁塔和通量观测数据，基于中尺度气象模式WRF (Weather Research and Forecast v3.7.1)，分析了5种边界层参数化方案在古尔班通古特沙漠的适用性。结果表明：1）采用WRF模拟沙漠腹地近地层内的边界层特征时，2m气温的模拟存在冷偏差，5种边界层参数化方案均能较好地模拟出四个季节2m气温的日变化特征，其中非局地方案ACM2(Asymmetric Convective Model version 2)对2m气温效果最好，局地方案BL方案的模拟偏差最大；2）5种边界层参数化方案均能够模拟出10m风速的日变化特征，其中局地方案BL(Bougeault-Lacarrere)对10m风速效果最佳；3）采用WRF模拟沙漠近地层内的地表通量特征时，感热通量存在高估现象，潜热通量存在低估现象，5种边界层参数化方案均能较好地模拟出四个季节模拟时间段内地表净辐射通量的日变化特征，其中局地方案MYJ(Mellor-Yamada-Janjie)的模拟精度最高。     

Evaluation of Two Initialization Schemes for Simulating the Rapid Intensification of Typhoon Lekima(2019)

Two different initialization schemes for tropical cyclone(TC) prediction in numerical models are evaluated based on a case study of Typhoon Lekima(2019). The first is a dynamical initialization(DI) scheme where the axisymmetric TC vortex in the initial conditions is spun up through the 6-h cycle runs before the initial forecast time. The second scheme is a bogussing scheme where the analysis TC vortex is replaced by a synthetic Rankine vortex. Results show that although both initialization schemes can help improve the simulated rapid intensification(RI) of Lekima, the simulation employing the DI scheme(DIS) reproduces better the RI onset and intensification rate than that employing the bogussing scheme(BOG).Further analyses show the cycle runs of DI help establish a realistic TC structure with stronger secondary circulation than those in the control run and BOG, leading to fast vortex spinup and contraction of the radius of maximum wind(RMW).The resultant strong inner-core primary circulation favors precession of the midlevel vortex under the moderate vertical wind shear(VWS) and thus helps vortex alignment, contributing to an earlier RI onset. Afterwards, the decreased vertical shear and the stronger convection inside the RMW support the persistent RI of Lekima in DIS. In contrast, the reduced VWS is not well captured and the inner-core convection is weaker and resides farther away from the TC center in BOG,leading to slower intensification. The results imply that the DI effectively improves the prediction of the inner-core process,which is crucial to the RI forecast.     

利用TWP-ICE试验资料对比两种边界层参数化方案

利用高分辨率WRF单气柱模式,选取了两种边界层参数化方案(YSU,MYJ),对TWP-ICE(Tropical Warm Pool International Cloud Experiment)试验期间的个例进行数值模拟,比较了两种方案对边界层结构、云和降水模拟的影响。结果表明:季风活跃期,YSU方案模拟的湍流交换系数较小,湍流混合偏弱,边界层内热通量偏小,使地表热量和水汽不易向上输送,水汽含量在近地表明显偏多,而在边界层及其以上大气层具有显著的干偏差,因此该方案模拟的云中液态水和固态水含量偏低,云量偏少,降水率偏小;MYJ方案对于季风活跃期的边界层结构具有较好的模拟能力,其模拟的云和降水更为准确。季风抑制期,MYJ方案模拟的夜间边界层结构存在较大误差,这是因为该方案模拟的夜间湍流交换系数较大,湍流混合偏强,边界层内热通量偏大,模拟的位温和水汽混合比在边界层内随高度变化较小,而观测廓线在边界层内存在较大梯度。季风抑制期两种方案模拟的云和降水均比观测值偏多,方案之间的差异较小。     

宁夏两次低涡暴雨过程边界层方案的敏感性试验

利用中尺度数值模式WRF选取3种边界层参数化方案（YSU、MYJ和ACM2）对2018年8月20日和2019年8月2日宁夏两次低涡暴雨过程进行敏感性试验,对比分析不同边界层参数化方案对暴雨过程模拟的影响。结果表明：YSU和MYJ方案对暴雨的预报性能均较优,ACM2方案对强降水的空报率较大;YSU方案对逐时中等强度降水的模拟效果更优。YSU方案对低涡系统的移动路径和中心强度模拟效果最好。相比ACM2方案,YSU和MYJ方案对环境场物理量和边界层内水汽混合比的模拟更接近实况,并且YSU方案能准确地反映边界层内位温分布特征。YSU与ACM2方案的湍流混合强度更强,后者维持较强湍流的时间更长。总体上,YSU方案对宁夏低涡暴雨过程的模拟效果最优。     

热带加热异常影响冬季平流层极涡强度的数值模拟

本文利用大气环流模式SAMIL/LASG,通过选择两种对流参数化方案,研究了热带加热异常对热带外平流层模拟的影响。结果表明,因不同对流参数化方案引起的热带对流加热状况的差异,可显著影响模式对北半球冬季平流层极涡强度的模拟偏差。与采用Manabe对流参数化方案相比,采用Tiedtke参数化方案可以显著改善对平流层极涡强度的模拟,使平流层极涡“过强”及极区“过冷”的模拟偏差得到明显改善。研究其中的影响过程发现,由于Manabe方案最大凝结潜热加热高度过低,在对流层中低层;而Tiedtke方案的最大凝结潜热加热位置在对流层中上层,因而Tiedtke(Manabe)方案时热带大气温度在对流层中上层较为偏暖(偏冷),在平流层低层较为偏冷(偏暖)。自秋季开始,与热带对流层高层温度的暖偏差相联系,热带外对流层高层以及热带平流层低层出现伴随的温度冷偏差;与之对应,平流层中纬度从秋季开始也出现持续的温度暖偏差。另外,随着秋冬季节平流层行星波活动的出现,Tiedtke方案时热带外地区行星波1波的强度也明显强于Manabe方案,使得秋冬季节涡动引起的向极热通量在Tiedtke方案时明显偏强,从而造成了冬季平流层极区温度偏暖、极涡强度偏弱。     

东北冷涡数值预报系统建立与业务应用

东北冷涡强降水一直是东北地区气象预报的难点和重点。从提高东北冷涡降水数值预报质量着手,选择模式中可能对东北冷涡强对流有影响的不同辐射方案、边界层方案和积云参数化方案等多种物理过程,对东北冷涡降水预报的性能进行批量试验;利用天气学检验方法,选择预报结果较好的RRTM长波辐射方案、Gayno Seaman边界层方案、Grell积云参数化方案和混合相湿过程显式方案等物理方案,组合建立了东北冷涡数值预报系统。经检验表明：系统对东北冷涡降水预报的水平较东北区域中尺度数值预报业务系统有所提高。同时利用天气学检验方法,从降水预报的强度、落区和移速等方面着手,评估东北冷涡数值预报产品对冷涡降水预报的性能和系统误差,得出降水强度易偏强、落区易偏西偏南和范围易偏小等定性结论,可直接应用于预报业务,取得了较好的应用效果。     

边界层参数化和湿过程对切变线低涡发展影响的中尺度模拟

1981年8月14至22日,在青藏高原东侧陕、甘、川毗邻区,出现了一次持续暴雨过程.对其主过程的涡度诊断表明,切变线低涡初生於青藏高原东侧低空,其后,不仅出现了高、低空涡度中心的叠加和耦合过程,而且在行星边界层(PBL)内发生了迅速加强过程.这种过程与切变低涡的持续发展相伴. 一系列48和72小时中尺度数值模拟结果揭示:只有高分辨PBL参数化的湿物理过程,才能较好地模拟出切变低涡在PBL内的突然加强及其涡度演变和垂直结构;而总体PBL参数化,虽然其涡度演变和垂直结构也能基本模拟出来,但不能模拟出切变线低涡在PBL内的迅速加强;总体PBL参数化干过程的模拟结果指出,该方案较前两种湿过程方案均差,特别是被模拟系统的强度明显减弱,位置偏离. 模拟与诊断和观测之间的比较揭示,尽可能仔细地改进PBL参数化和湿物理过程,才有可能改进对这类中尺度系统发生和发展的模拟.     

海洋飞沫对台风“Morakot”结构影响的数值模拟研究

将海洋飞沫参数化引入到高分辨率、非静力中尺度模式中,并对0908号台风"Morakot"进行了数值模拟,研究了海洋飞沫对台风"Morakot"结构和强度的影响。结果表明:(1)不论是否考虑海洋飞沫作用,模式均能较好地模拟出台风"Morakot"的移动路径,说明海洋飞沫对其移动路径影响不大;(2)引入海洋飞沫参数化后,台风眼墙区域的切向风速、径向风速、垂直速度、涡度、云水混合比、雨水混合比等物理量均增强,表明飞沫对台风结构变化的影响明显;(3)海洋飞沫对台风"Morakot"演变的直接影响是在对流层低层,低层风速明显增大,大风速区的影响尤为显著;(4)飞沫的蒸发使台风范围内的潜热和感热通量明显增强,尤其是潜热通量,其大值区对应着台风中心附近的最大风速区。由于水汽和热量输送的增强,使台风眼壁附近的云水量与雨水量增多,因此降水强度明显增加。     

Numerical Simulation of the Rapid Intensification of Hurricane Katrina(2005):Sensitivity to Boundary Layer Parameterization Schemes

Accurate forecasting of the intensity changes of hurricanes is an important yet challenging problem in numerical weather prediction. The rapid intensification of Hurricane Katrina(2005) before its landfall in the southern US is studied with the Advanced Research version of the WRF(Weather Research and Forecasting) model. The sensitivity of numerical simulations to two popular planetary boundary layer(PBL) schemes, the Mellor–Yamada–Janjic(MYJ) and the Yonsei University(YSU) schemes, is investigated. It is found that, compared with the YSU simulation, the simulation with the MYJ scheme produces better track and intensity evolution, better vortex structure, and more accurate landfall time and location. Large discrepancies(e.g.,over 10 hPa in simulated minimum sea level pressure) are found between the two simulations during the rapid intensification period. Further diagnosis indicates that stronger surface fluxes and vertical mixing in the PBL from the simulation with the MYJ scheme lead to enhanced air–sea interaction, which helps generate more realistic simulations of the rapid intensification process. Overall, the results from this study suggest that improved representation of surface fluxes and vertical mixing in the PBL is essential for accurate prediction of hurricane intensity changes.     

边界层参数化方案对南海秋季台风“莎莉嘉”(2016)模拟的影响

利用美国中尺度数值模式WRF,选取两个局地(QNSE、MYJ)闭合和两个非局地(YSU、ACM2)闭合的边界层参数化方案对台风“莎莉嘉”(2016)进行了4组模拟试验,结果表明,不同边界层方案对台风路径影响较小,但对台风强度和结构有明显的影响,就本个例研究而言,非局地闭合边界层方案明显优于局地闭合边界层方案。台风强度的差异是热力和动力共同作用的影响。局地闭合方案模拟的地表焓通量、水汽通量和动量通量更大,台风偏强;局地闭合方案模拟的边界层高度更高、边界层顶的夹卷过程更强、垂直混合更强、台风暖心结构更强,从而台风也更强。台风强度的差异和台风结构的变化密切相关。      

A note on surface layer parameterizations in the weather research and forecast model

Two surface layer parameterization schemes along with five planetary boundary layer (PBL) schemes in the Weather Research and Forecasting model (WRF) are analyzed in order to evaluate the performance of the WRF model in simulating the surface variables and turbulent fluxes over an Indian sub-continent region. These surface layer schemes are based on the fifth-generation Pennsylvania State University—National Center for Atmospheric Research Mesoscale Model (MM5) parameterization; (a) Old MM5 scheme having Businger-Dyer similarity functions and (b) revised MM5 scheme utilizing the functions that are valid for full ranges of atmospheric stabilities. The study suggests that each PBL scheme can reproduce the diurnal variation of 2 m temperature, momentum flux and sensible heat flux irrespective of the surface layer scheme used for the simulations. However, a comparison of model-simulated values of surface variables and turbulent fluxes with observed values suggests that each PBL scheme is found to systematically over-estimate the nocturnal 2 m temperature and 10 m wind speed with both the revised and old schemes during stable conditions.     

初始涡旋结构对热带气旋强度变化影响的数值研究

本文利用中尺度WRF模式,通过构造3个位于不同高度上强度相同的初始涡旋暖心中心(分别称为Low试验、Mid试验和High试验),认识暖心垂直结构对热带气旋发展的影响。理想数值试验发现,在积分6 h后在Low试验和Mid试验中最大风速半径开始收缩,眼墙内对流发展,高层暖心发展明显比High试验强,强度增加明显快于High试验,达到快速增强的标准。进一步诊断发现,暖心偏低的试验中初始CAPE值较大,低层风垂直切变较强,有利于眼墙内对流发展,非绝热加热在高层暖心形成过程中起重要作用,最大风速半径收缩比High试验快,热带气旋强度快速增加。本研究清楚表明,数值预报模式中构造初始涡旋的暖心高度对模拟热带气旋的强度发展有重要影响。     

2006年超级台风“桑美”强度与结构变化的数值模拟研究

使用一个高分辨率、非静力数值模式WRF模式对2006年超级台风Saomei强度和结构进行了数值模拟研究.首先,评估了Makin的粗糙度长度公式对台风Saomei强度和结构变化的影响,结果表明,采用新参数后,使得模拟的台风强度变化与实况最佳路径资料的强度变化更一致,对超级台风Saomei强度预报有改进;但对台风路径的影响不大.通过QuikSCAT、雷达和TRMM非常规资料的验证,进一步表明模拟的台风Saomei的结构与实况很接近,可以再现台风内核区域的部分"双眼墙"和"Annular"结构.其次,通过对台风Saomei边界层过程模拟的改进,表明在平均风速大于40 m/s时边界层各物理量明显改善,使得模式最大强度比传统的简单外推插值方案有显著改进,特别是在台风最强阶段,当台风Saomei眼墙区域的海表面拖曳系数C_d的相对变小,使得其眼墙区域的平均切向风速、径向风速、垂直风速、温度距平、涡旋动能和绝对角动量等物理量均有增强.表明台风Saomei眼墙氏域(20-40 km)各物理量的贡献对其强度和结构变化的影响十分重要.最后,在此基础上进一步分析模式海温和大尺度环境垂直风切变对台风Saomei强度和结构变化的可能影响,讨论了台风Saomei在其增强和消弱阶段中,大尺度环境垂直风切变对其强度变化的负反馈作用.     

夏季区域极涡异常对塔里木河流域降水的影响

采用1961-2015年夏季大西洋-欧洲极涡面积指数和塔里木河流域43个站降水资料,研究该区域极涡面积异常对该流域降水的影响。结果表明:两者年际变化呈显著的反位相关系;在极涡面积异常偏小（大）年,西风急流在西亚和中亚减弱（增强）,在东亚则相反;500 hPa欧洲中部和贝加尔湖地区的高压脊偏强（偏弱）,东亚低压槽偏弱（偏强）,中亚经向环流增强（减弱）;700 hPa塔里木河流域天气扰动活跃（不活跃）,东风和西南风（西北风）的水汽输送增强,西部和北部等主要降水区水汽辐合（辐散）,该流域降水偏多（偏少）。在极涡面积异常偏小年的夏季,塔里木河流域水平风场和垂直运动从高纬度到低纬度的经向变化分布与大西洋-欧洲区相似,大西洋-欧洲极涡区与塔里木河流域之间存在西北-东南的环流异常分布,大西洋-欧洲极涡面积异常可能通过该环流异常分布影响塔里木河流域及周边风场、水汽输送和垂直运动,进而影响到该流域降水。