不同边界层参数化方案对台风“利奇马”模拟的影响研究

基于中尺度数值模式WRF,对比分析了六种大气边界层物理过程参数化方案(BouLac、MYJ、UW、YSU、ACM2、SH)对台风“利奇马”模拟结果的影响。结果表明,不同边界层方案对“利奇马”路径的模拟结果影响较小,但对其强度和结构演变的模拟结果影响显著。其中,局地闭合方案UW方案模拟的结果最强,局地闭合方案BouLac次之,而局地闭合方案MYJ和三种非局地闭合方案YSU、ACM2和SH的模拟强度都较弱。这些方案中,BouLac模拟的海平面最低气压与实况最为接近。通过对比这些边界层方案的模拟结果发现,由于台风强度的差异受到热力和动力的共同影响,边界层方案如模拟得到的地表潜热通量和边界层中湍流扩散系数较大,将导致较大的径向风和低层辐合,从而模拟得到较强的台风强度;反之,则台风强度较弱。      

边界层参数化对海南岛海风环流结构模拟的影响

利用WRF V3.7详细分析了应用8种边界层参数化方案（YSU、MYNN2.5、MYNN3、ACM2、BouLac、UW、SH、GBM）所模拟的2014年5月25日海南岛海风环流结构的差异,其中YSU、ACM2和SH为非局地闭合方案,MYNN2.5、MYNN3、BouLac、UW和GBM为局地闭合方案。结果表明:对于海风环流水平结构的模拟,15时,YSU、ACM2、BouLac、UW和SH模拟的北部海风较强,SH和GBM的内陆风速偏大。温度与海风发展强度相对应,MYNN2.5与MYNN3模拟的岛屿温度偏低,海陆温差小,海风相对较弱。对于海风环流垂直结构的模拟,09时海风开始,但强度较小,且存在残余陆风,向内陆传播距离较短,YSU、MYNN2.5和SH方案的海风相对较强。12时,海风已呈现出较为清晰的环流结构,YSU和ACM2的海风厚度及向内陆传播距离相对强于其它方案,MYNN3的环流结构则不太明显,且向内陆推进距离短,海风相对较弱。15时,海风发展强盛,MYNN2.5和MYNN3方案模拟的海风垂直强度较小,ACM2方案的海风垂直环流特征最为明显。18时,海风的强度和扰动均有所减弱,ACM2、BouLac和UW的整体海风相对强于其它方案。21时海风已基本转为陆风,BouLac与UW的陆风环流结构最为清晰。位温、水汽及海风垂直环流强度的发展变化与海风的演变过程基本一致。造成ACM2模拟海风偏强的原因是其边界层垂直混合偏强,形成了足够的湍流混合强度所致。对于边界层高度的模拟,ACM2的边界层顶最高,这与此方案所模拟的海风强度偏大相吻合,其它方案的边界层高度与海风强度并不完全一致。      

边界层方案对华北低层O3垂直分布模拟的影响

利用WRF-Chem模式,采用3种边界层参数化方案 (YSU, MYJ和ACM2),针对1个晴空、静稳日 (2013年8月26日20:00—27日20:00(北京时)) 进行模拟,着重分析不同边界层参数化方案对夜间残留层形成及日出前后O₃浓度垂直分布形式的模拟效果,并与固城站地面及垂直同步观测资料进行对比。结果表明:3种边界层参数化方案均能够模拟出温度及风速的区域分布形式以及风温垂直结构的变化特征;相比之下,MYJ方案模拟的夜间边界层高度较YSU方案和ACM2方案明显偏高,该对比结果可能是导致近地面污染物浓度模拟差异的重要原因;在夜间稳定层结至日出后稳定状态打破的边界层结构演变过程中,采用YSU方案和ACM2方案模拟的温度和风速垂直扩线形式与观测结果更为接近;同样采用非局地闭合的YSU方案和同时考虑局地和非局地闭合的ACM2方案,对于边界层高度内O₃浓度垂直分布形式的模拟效果具有明显优势。     

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

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

XJ-RUCv1.0系统不同边界层方案对地面气象要素预报的影响初探

基于沙漠绿洲戈壁区域同化预报系统（Desert Oasis Gobi Regional Assimilation Forecast System,简称DOGRAFSv1.0）不同边界层方案（YSU,ACM2）对2013年1月和7月的预报结果和实况资料,利用客观检验方法,初步探究不同边界层方案对地面气象要素预报的影响,结果表明：ACM2方案对地面气象要素预报效果总体优于YSU方案,尤其对1月份ACM2方案预报结果优势显著;对于2 m温度,ACM2方案对低温预报偏高的问题有一定改善,但对高温预报偏低的问题有所加重;对于10 m风速,ACM2方案的预报误差明显低于YSU方案;对于降水,ACM2方案在一定程度上改善了1月漏报和空报的问题,但对7月改进不明显,这主要由于ACM2方案在稳定或中性层结下采用的局地闭合算法削弱了边界层顶卷挟和边界层内混合作用,进而抑制了边界层物理过程与对流过程的耦合作用。     

Sensitivity of the regional climate model RegCM4.2 to planetary boundary layer parameterisation

This study investigates the performance of two planetary boundary layer (PBL) parameterisations in the regional climate model RegCM4.2 with specific focus on the recently implemented prognostic turbulent kinetic energy parameterisation scheme: the University of Washington (UW) scheme. When compared with the default Holtslag scheme, the UW scheme, in the 10-year experiments over the European domain, shows a substantial cooling. It reduces winter warm bias over the north-eastern Europe by 2 °C and reduces summer warm bias over central Europe by 3 °C. A part of the detected cooling is ascribed to a general reduction in lower tropospheric eddy heat diffusivity with the UW scheme. While differences in temperature tendency due to PBL schemes are mostly localized to the lower troposphere, the schemes show a much higher diversity in how vertical turbulent mixing of the water vapour mixing ratio is governed. Differences in the water vapour mixing ratio tendency due to the PBL scheme are present almost throughout the troposphere. However, they alone cannot explain the overall water vapour mixing ratio profiles, suggesting strong interaction between the PBL and other model parameterisations. An additional 18-member ensemble with the UW scheme is made, where two formulations of the master turbulent length scale in unstable conditions are tested and unconstrained parameters associated with (a) the evaporative enhancement of the cloud-top entrainment and (b) the formulation of the master turbulent length scale in stable conditions are systematically perturbed. These experiments suggest that the master turbulent length scale in the UW scheme could be further refined in the current implementation in the RegCM model. It was also found that the UW scheme is less sensitive to the variations of the other two selected unconstrained parameters, supporting the choice of these parameters in the default formulation of the UW scheme.     

边界层参数化方案对“莫兰蒂”台风（1614）登陆阶段影响的数值模拟研究

利用中尺度数值模式WRF v3.8中的YSU、MYJ、QNSE、ACM2、UW、GBM、Boulac七种不同边界层参数化方案,采用高分辨率（1.33 km）数值试验的方法研究了不同边界层方案对模拟台风“莫兰蒂”（1614）登陆减弱阶段的移动路径、强度、结构、降水量、近地层有关物理量场分布等方面的影响,结果表明:（1）“莫兰蒂”台风登陆减弱阶段,不同边界层方案对台风路径、强度、降水量模拟影响显著,24 h内模拟台风路径、最低气压、最大风速及24 h累积降水量极值的最大差异分别达80 km、11 hPa、27 m s?1及241 mm;（2）Boulac方案模拟台风路径与实况最为接近,GBM、YSU和MYJ方案分别次之,ACM2和UW方案再次之,而QNSE方案最差;UW和QNSE方案模拟的最低气压以及MYJ和QNSE方案模拟的最大风速与观测最为接近;不同边界层方案均模拟出台风登陆阶段最低气压逐渐升高以及其升高速率在台风登陆后大于登陆前的特征,这与实况一致,但台风登陆前各方案模拟最低气压升高速度均大于实况,而台风登陆后却又不及实况;（3）Boulac方案模拟的24 h降水分布、强降水落区、结构、强度和各量级降水TS评分均最优,MYJ方案次之;而QNSE、UW和ACM2方案雨带向西北方向推进过快,各量级降水TS评分均较差;（4）综合台风路径、强度和降水模拟,Boulac和MYJ方案相对最优,其中Boulac方案在台风路径和降水模拟上更优,而MYJ方案在台风强度模拟上更优;YSU和GBM方案次之,而QNSE、UW和ACM2方案相对较差;（5）不同边界层方案计算的近地层潜热通量、感热通量显著不同,进而影响台风路径、强度、降水量模拟存在显著差异。比较而言,QNSE方案潜热通量相对异常偏高,MYJ和Boulac方案量值适中,其余方案相对偏低;QNSE方案感热通量相对略偏高,MYJ方案适中,其他方案则相对显著偏低;（6）不同边界层方案模拟降水区边界层热、动力结构显著不同,其中Boulac方案具有较明显优势,尤其是对日间边界层结构的模拟。     

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

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

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

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

Numerical simulations of heavy rainfall over central Korea on 21 September 2010 using the WRF model

On 21 September 2010, heavy rainfall with a local maximum of 259 mm d-1occurred near Seoul, South Korea. We examined the ability of the Weather Research and Forecasting(WRF) model in reproducing this disastrous rainfall event and identified the role of two physical processes: planetary boundary layer(PBL) and microphysics(MPS) processes. The WRF model was forced by 6-hourly National Centers for Environmental Prediction(NCEP) Final analysis(FNL) data for 36 hours form 1200 UTC 20 to 0000 UTC 22 September 2010. Twenty-five experiments were performed, consisting of five different PBL schemes—Yonsei University(YSU), Mellor-Yamada-Janjic(MYJ), Quasi Normal Scale Elimination(QNSE),Bougeault and Lacarrere(Bou Lac), and University of Washington(UW)—and five different MPS schemes—WRF SingleMoment 6-class(WSM6), Goddard, Thompson, Milbrandt 2-moments, and Morrison 2-moments. As expected, there was a specific combination of MPS and PBL schemes that showed good skill in forecasting the precipitation. However, there was no specific PBL or MPS scheme that outperformed the others in all aspects. The experiments with the UW PBL or Thompson MPS scheme showed a relatively small amount of precipitation. Analyses form the sensitivity experiments confirmed that the spatial distribution of the simulated precipitation was dominated by the PBL processes, whereas the MPS processes determined the amount of rainfall. It was also found that the temporal evolution of the precipitation was influenced more by the PBL processes than by the MPS processes.     

基于多种边界层方案的天津PM2.5集合预报试验

基于大气化学模式WRF/Chem采用4种边界层方案（YSU,BL,MYJ和MYN3）模拟2015年全年天津地区细颗粒物质量浓度演变趋势,分析不同边界层方案对天津地区PM_2.5质量浓度模拟和预报的影响,并构建多种边界层方案的集合预报产品,以期提高天津地区PM_2.5质量浓度预报效果。结果表明:大气化学模式的4种边界层方案在空气质量模拟中均有较好的适用性,PM_2.5模拟值与观测值相关系数达到0.76左右,晴空和大风天气BL方案表现较优,阴天和小风天气YSU和MYJ方案表现优于其他两种方案,综合所有过程没有任何单一方案具有显著优势。基于上述原因,该文开展多种边界层方案的天津空气质量集合预报试验,经过对2015年全年模拟值分析,多边界层方案和多气溶胶机制扰动的集合预报可以减小PM_2.5质量浓度预报的相对误差和均方根误差,降低重污染天气预报的漏报率,提升空气质量等级预报能力,在计算资源得到保证的基础上,是一种可以提升数值模式PM_2.5质量浓度预报能力的有效手段。     

利用激光雷达结合数值模式估算兰州远郊榆中地区夏季边界层高度

用激光雷达资料,采用小波变化法反演兰州远郊榆中地区兰州大学半干旱气候与环境观测站(SACOL)的边界层高度,并利用WRF中尺度数值模式,选取两种不同边界层参数化方案(YSU、MYJ)模拟了该地区边界层及其高度.分析表明激光雷达反演边界层高度与WRF模拟边界层高度结果基本一致;WRF选用YSU方案能较好反映热对流边界层,而MYJ方案对于动力作用边界层模拟较好.日出后08:00(北京时间,下同)SACOL不稳定边界层开始发展,17:00达到最大高度.热对流边界层可以达到2 km;动力作用边界层可达到1.5 km,之后热对流边界层下降速度明显高于动力作用边界层.     

Improvement of the K-profile Model for the Planetary Boundary Layer based on Large Eddy Simulation Data

Modifications of the widely used K-profile model of the planetary boundary layer (PBL), reported by Troen and Mahrt (TM) in 1986, are proposed and their effects examined by comparison with large eddy simulation (LES) data. The modifications involve three parts. First, the heat flux from the entrainment at the inversion layer is incorporated into the heat and momentum profiles, and it is used to predict the growth of the PBL directly. Second, profiles of the velocity scale and the Prandtl number in the PBL are proposed, in contrast to the constant values used in the TM model. Finally, non-local mixing of momentum was included. The results from the new PBL model and the original TM model are compared with LES data. The TM model was found to give too high PBL heights in the PBL with strong shear, and too low heights for the convection-dominated PBL, which causes unrealistic heat flux profiles. The new PBL model improves the predictability of the PBL height and produces profiles that are more realistic. Moreover, the new PBL model produces more realistic profiles of potential temperature and velocity. We also investigated how each of these three modifications affects the results, and found that explicit representation of the entrainment rate is the most critical.     

春季黄海海雾WRF参数化方案敏感性研究

利用2005—2011年10次春季黄海海雾个例开展WRF模式参数化方案敏感性研究。结果表明:边界层方案对WRF模式雾区模拟结果起决定作用,而微物理方案影响较小,它主要影响海雾浓度和高度。边界层与微物理方案的最佳组合为YSU与Lin方案,最差为Mellor-Yamada与WSM5方案;Mellor-Yamada和QNSE方案模拟的近海面湍流过强,导致边界层过高,不利于海雾的发展与维持;而MYNN与YSU方案刻画的湍流强度与边界层高度合适,有利于海雾发展与维持。MYNN方案虽与YSU方案相当,但在大多数海雾个例中,后者明显优于前者,而在有些个例中却刚好相反。因此对于某一具体海雾个例而言,所用边界层方案仍需在它们之中选择最优者。这些信息可为黄海海雾WRF模式边界层与微物理方案的选择与改进提供参考。     

WRF模式多种边界层参数化方案对四川盆地不同量级降水影响的数值试验

利用中尺度模式WRF三种边界层参数化方案(MYJ、YSU和ACM2),对2012年四川盆地夏季连续40天逐日降水量进行数值试验,并检验评估了不同边界层参数化方案下模式对分级降水量和边界层结构的模拟能力,分析了各参数化方案对降水量模拟差异的可能原因。结果表明:三种边界层参数化方案对较小量级(小雨和中雨)降水量的模拟,24 h时效优于48 h,ACM2方案效果较好;对较大量级(大雨和暴雨)降水的模拟,48 h时效优于24 h,YSU方案模拟效果较好。对比分析温江站加密探空观测与模式模拟的大气边界层结构表明,ACM2方案对小量级降水时边界层结构的模拟较为准确,而YSU方案更适合于温江站大量级降水时边界层结构的模拟。不同边界层参数化方案对各量级降水量模拟差异的可能原因是边界层湍流混合强度的不同,MYJ方案湍流混合作用较弱,导致底层大量水汽积聚,不稳定性强,容易产生虚假降水,因此对各量级降水模拟能力均有限;YSU方案具有强烈的垂直混合强度,有利于局地水汽的向上输送,更易达到大量级降水发生发展的条件,适用于盆地较大量级降水的模拟;ACM2方案在保证足够湍流混合强度的同时,在较稳定条件下会关闭非局地输送,不致于产生过强降水,适合盆地较小量级降水的数值模拟     

WRF模式对青藏高原那曲地区大气边界层模拟适用性研究

采用WRF（Weather Research and Forecasting）模式4种边界层参数化方案对青藏高原那曲地区边界层特征进行了数值模拟,并利用"第三次青藏高原大气科学试验"在青藏高原那曲地区5个站点的观测资料对模拟结果进行验证,分析不同参数化方案在那曲地区的适用性。研究表明,YSU、MYJ、ACM2和BouLac方案对2 m气温和地表温度的模拟偏低。BouLac方案模拟的地表温度偏差较小。通过对能量平衡各分量的对比分析发现,温度模拟偏低可能是向下长波辐射模拟偏低以及感热通量和潜热通量交换过强导致的。对于边界层风、位温和相对湿度垂直结构的模拟,局地方案的模拟效果均优于非局地方案。BouLac方案对那曲地区近地层温度、边界层内位温和相对湿度的垂直分布模拟效果较好。      

边界层参数化方案在“灰色区域”尺度下的适用性评估

随着数值预报模式分辨率的提高,当模式网格距与含能湍涡的长度尺度相当时,模式动力过程可解析一部分湍流运动,而剩余的湍流运动仍需参数化,此时便产生了湍流参数化的“灰色区域”问题。对传统的PBL（Planetary Boundary Layer）方案在“灰色区域”下的适用性评估,是改进PBL方案以使其能够适应分辨率变化的前提和基础。本研究基于干对流边界层的大涡模拟试验,比较了WRF（Weather Research and Forecast Model）模式中四种常用的边界层参数化方案[YSU（Yonsei University）、MYJ（Mellor-Yamada-Janjic）、MYNN2.5（Mellor-Yamada-Nakanishi-Niino Level 2.5）、MYNN3）]在“灰色区域”尺度下的表现。研究表明,混合层内总热通量对所使用的参数化方案和水平分辨率均不敏感。不同参数化方案中次网格与网格通量的比例表现出对水平网格距不同的依赖性。局地PBL方案（MYJ、MYNN2.5）在混合层内的平均位温随网格距减小而增大,次网格通量随网格距减小而减小,较参考湍流场对次网格通量有所低估。YSU方案的非局地项几乎不随水平格距改变而变化,对次网格通量的表征并未表现出较强的分辨率依赖性,且过强的非局地次网格输送使混合层内温度层结呈弱稳定,抑制了可分辨湍流输送,不易于激发次级环流。MYNN3方案的非局地次网格通量（负梯度输送项）随网格距减小而减小,使其对次网格通量的表征具有较好的分辨率依赖性。PBL方案在“灰色区域”尺度下的适用性与具体分辨率有关。以分辨率500 m为例,四种PBL方案中不存在一种最佳方案,能对边界层的热力结构和湍流统计特征均有准确的描述。     

Intercomparison of Planetary Boundary-Layer Parametrizations in the WRF Model for a Single Day from CASES-99

This study compares five planetary boundary-layer (PBL) parametrizations in the Weather Research and Forecasting (WRF) numerical model for a single day from the Cooperative Atmosphere-Surface Exchange Study (CASES-99) field program. The five schemes include two first-order closure schemes—the Yonsei University (YSU) PBL and Asymmetric Convective Model version 2 (ACM2), and three turbulent kinetic energy (TKE) closure schemes—the Mellor–Yamada–Janjić (MYJ), quasi-normal scale elimination (QNSE), and Bougeault–Lacarrére (BouLac) PBL. The comparison results reveal that discrepancies among thermodynamic surface variables from different schemes are large at daytime, while the variables converge at nighttime with large deviations from those observed. On the other hand, wind components are more divergent at nighttime with significant biases. Regarding PBL structures, a non-local scheme with the entrainment flux proportional to the surface flux is favourable in unstable conditions. In stable conditions, the local TKE closure schemes show better performance. The sensitivity of simulated variables to surface-layer parametrizations is also investigated to assess relative contributions of the surface-layer parametrizations to typical features of each PBL scheme. In the surface layer, temperature and moisture are more strongly influenced by surface-layer formulations than by PBL mixing algorithms in both convective and stable regimes, while wind speed depends on vertical diffusion formulations in the convective regime. Regarding PBL structures, surface-layer formulations only contribute to near-surface variability and then PBL mean properties, whereas shapes of the profiles are determined by PBL mixing algorithms.     

Parameterization of sheared entrainment in a well-developed CBL. Part I: Evaluation of the scheme through large-eddy simulations

The entrainment flux ratio A _e and the inversion layer (IL) thickness are two key parameters in a mixed layer model. A _e is defined as the ratio of the entrainment heat flux at the mixed layer top to the surface heat flux. The IL is the layer between the mixed layer and the free atmosphere. In this study, a parameterization of A _e is derived from the TKE budget in the firstorder model for a well-developed CBL under the condition of linearly sheared geostrophic velocity with a zero value at the surface. It is also appropriate for a CBL under the condition of geostrophic velocity remaining constant with height. LESs are conducted under the above two conditions to determine the coefficients in the parameterization scheme. Results suggest that about 43% of the shear-produced TKE in the IL is available for entrainment, while the shear-produced TKE in the mixed layer and surface layer have little effect on entrainment. Based on this scheme, a new scale of convective turbulence velocity is proposed and applied to parameterize the IL thickness. The LES outputs for the CBLs under the condition of linearly sheared geostrophic velocity with a non-zero surface value are used to verify the performance of the parameterization scheme. It is found that the parameterized A _e and IL thickness agree well with the LES outputs.