诊断分析在数值预报产品释用中的探讨

针对地方性天气预报特点，根据现有数值预报产品及实时探空资料的接收情况，作者利用中尺度数值预报模式ＭＭ４引入中尺度信息，提出一个实际可行的夏季降水预报方案，以求提高区域范围内具有明显突发性对流降水的预报能力。作法是用ＭＭ４模式计算逐日物理量诊断场，加入同期的欧洲中心数值预报产品，共计４１个场作为预报因子场；应用典型相关分析法，并结合逐步判别和逐步回归法，计算得出９个分区和８个单站中雨以上降水预报方程。降水预报方程的因子变量多数为中尺度诊断分析场构成的典型团子，可以认为，ＭＭ４模式计算得出的中尺度信息对预报中雨以上降水是有贡献的。     

诊断分析在数值预报产品释用中的探讨

针对地方性天气预报特点，根据现有数值预报产品及实时控空资料的接收情况，作者利用中尺度数值预报模式ＭＭ４引入了中尺度信息，提出了一个实际可行的夏季降水预报方案，以求提高区域范围内具有明显突发性对流降水的预报能力。     

一次罕见的华南大暴雨过程的 诊断与数值模拟研究

首先对1994年6月中旬的华南暴雨过程进行了诊断分析,然后,利用PSU和NCAR发展起来的有限区域中尺度模拟系统MM5进行了数值模拟研究。诊断结果表明:此次过程长江以南的大部分地区具备大尺度的降水条件,而在两广地区降水更强,这可能是由地形等其他条件决定的。两广地区活动较强的对流云团是造成暴雨的主要中尺度系统,而东部地区主要是锋面造成的连续性降水,对流活动较两广地区相对较弱。数值模拟结果显示:MM5模式对这次华南的强降水过程具有较好的模拟能力,静力方案对大范围雨区的模拟比非静力方案要好一些,但对大暴雨区的模拟非静力方案似乎更具优越性。对流参数化方案的试验表明,Kuo方案能较好地模拟大范围的雨区,Grell方案对强暴雨中心的模拟似乎比其他方案好。对行星边界层方案的试验显示,边界层过程对强的对流性降水的模拟是不可忽略的。这些看法与已有的诊断结果是相符的。     

两类不同风灾个例超级单体特征对比分析

采用分钟级加密自动气象站观测资料,盐城、淮安和岳阳、荆州雷达探测数据,以及欧洲中期天气预报中心（ECMWF）高分辨率的ERA-Interim全球再分析数据,对比分析了2016年6月23日江苏阜宁龙卷灾害和2015年6月1日湖北监利下击暴流大风灾害的环境特征与超级单体的结构特征。结果表明:（1）两次强对流大风灾害发生在相似的低空环流背景下:风灾发生在低空急流出口区左侧的暖区内、850 hPa低涡中心东侧6—7个经距的位置;环境大气的对流有效位能大于2000 J/kg。但是风灾的类型不同,江苏阜宁大风灾害主要由超级单体龙卷造成,监利“东方之星”沉船事故主要是超级单体触发的下击暴流造成。短时强降水中心与风灾中心的相对位置不同:阜宁龙卷移动方向的左侧伴随着最强短时降水;湖北监利沉船事件发生期间,风灾中心与短时强降水中心基本重合。鉴于不同性质的对流大风位置与超级单体母体的中心位置对应关系上存在差异,通过比较地面观测的瞬时大风与瞬时强降水中心的相对位置将有助于区分强对流大风的性质。（2）环境风垂直切变强度对对流风暴结构、发展、维持有重要影响:阜宁龙卷发生时,其上空0—6 km风垂直切变达4×10-3 s-1,超级单体有明显的向前倾斜结构,形成有界弱回波区;而监利强对流沉船位置0—6 km风垂直切变只有2.3×10-3 s-1左右,风暴单体中的上升气流近乎于垂直。阜宁超级单体中气旋,首先出现在0—1.5 km风垂直切变和0—3 km风暴相对螺旋度带状大值区,在向抬升凝结高度更低的环境移动过程中,其底部不断下降,形成龙卷;而在监利沉船区,中低层风切变和风暴相对螺旋度相对要弱得多,对应风暴单体中的中气旋强度、持续性较弱,中气旋底部高度维持在1.6 km左右。（3）环境湿度垂直结构特征不同可能是风暴单体形成不同类型灾害大风的重要环境因子。监利下击暴流造成的风灾发生时,在地面气温迅速下降过程中,气压变化呈现快速跳升又快速下降的“尖锥”形,气压峰值比降水峰值提前4 min出现。它与对流层中高层环境大气中较为深厚的干空气卷入对流风暴中造成水物质强烈蒸发、冷却过程有关。而阜宁风灾过程中,环境大气中层仅存在非常浅薄的干层,加之低层较为深厚的饱和大气环境,对应的地面冷池效应相对较弱。      

基于VDRAS的快速更新雷达四维变分分析系统

基于雷达资料快速更新四维变分同化 (RR4DVar) 技术和三维数值云模式,初步研发了一个针对对流尺度数值模拟的快速更新雷达四维变分分析系统。系统通过对京津冀6部多普勒天气雷达资料进行RR4DVar同化,并融合5 min自动气象站观测和中尺度数值模式结果,可快速分析得到12~18 min更新的低层大气三维动力、热力场的对流尺度结构特征。针对2009年7月22日发生在京津冀的一次强风暴个例,通过一系列敏感性试验,并利用局地加密资料进行检验对比,表明有效的雷达资料RR4DVar同化及自动气象站和中尺度模式资料融合方案、恰当的中尺度背景场设置和动力约束方法是获得合理结果的关键。研究也表明:恰当的系统配置能够模拟出与对流生消发展密切相关的近风暴环境特征,包括低层入流、垂直风切变、低层辐合上升和暖舌等,以及风暴自身形成的冷池、出流等与风暴演变密切相关的对流尺度结构。     

Radar-based Characteristics and Formation Environment of Supercells in the Landfalling Typhoon Mujigae in 2015

This study presents the radar-based characteristics and formation environment of supercells spawned by the tornadic landfalling Typhoon Mujigae(2015)in October 2015.More than 100 supercells were identified within a 24-hour period around the time of the typhoon’s landfall,of which three were tornadic with a rotational intensity clearly stronger than those of non-tornadic supercells.The identified supercells were concentrated within a relatively small area in the northeast quadrant beyond 140 km from the typhoon center.These supercells were found more likely to form over flat topography and were difficult to maintain in mountainous regions.During the study period,more supercells formed offshore than onshore.The mesocyclones of the identified supercells were characterized by a small diameter generally less than 5 km and a shallow depth generally less than 4 km above ground level.An environmental analysis revealed that the northeast quadrant had the most favorable conditions for the genesis of supercell in this typhoon case.The nondimensional supercell composite parameter(SCP)and entraining-SCP(E-SCP)were effective in separating supercell from non-supercell environment.Even though the atmosphere in the typhoon’s northeast quadrant was characterized by an E-SCP/SCP value supportive of supercell organization,orography was an impeditive factor for the supercell development.These findings support the use of traditional parameters obtained from midlatitude supercells to assess the supercell potential in a tropical cyclone envelope.     

The Uncertainty of Mesoscale Numerical Prediction of Heavy Rain in South China and the Ensemble Simulations?

In the context of non-hydrostatic MM5 version we have explored the impact of convective parameterization schemes on uncertainty in mesoscale numerical prediction of South China heavy rain and mesoscale heavy rainfall short-range ensemble simulation by using two kinds of physics perturbation methods through a heavy rain case occurring on June 8, 1998 in Guangdong and Fujian Provinces. The results show the physical process of impacts of convective schemes on heavy rainfall is that different latent heat of convective condensation produced by different convective schemes can make local temperature perturbation, leading to the difference of local vertical speed by the intrinsic dynamic and thermodynamic processes of atmosphere,and therefore, making difference of the timing, locations and strength of mesh scale and subgrid scale precipitation later. New precipitations become the new source of latent heat and temperature perturbation,which finally make the dynamic and thermodynamic structures different in the simulations. Two kinds of methods are used to construct different model version stochastically. The first one is using different convective parameterization and planetary boundary layer schemes, the second is adjusting different parameters of convective trigger functions in Grell scheme. The results indicate that the first ensemble simulations can provide more uncertainty information of location and strength of heavy rainfall than the second. The single determinate predictions of heavy rain are unstable; physics ensemble predictions can reflect the uncertainty of heavy rain, provide more useful guidance and have higher application value.Physics ensembles suggest that model errors should be taken into consideration in the heavy rainfall ensembles. Although the method of using different parameters in Grell scheme could not produce good results, how to construct the perturbation model or adjust the parameter in one scheme according to the physical meaning of the parameter still needs further investigation. The limitation of the current study is that it is based on a single case and more cases will be addressed in the future researches.     

垂直-倾斜对流一体化参数化方案的实现及数值试验

在Kuo-Anthes垂直对流参数化方案和Nordeng倾斜对流参数化方案基础上,提出了垂直-倾斜对流一体化参数化方案,并引入MM5模式中.利用该方案对2008年1月28-29日发生在中国南方的一次暴雪过程和2005年"海棠"台风过程进行了数值模拟,模拟结果表明,此次暴雪过程在垂直方向主要表现对流稳定状态,但在对流层低层始终存在条件件对称不稳定层,并且当条件性对称不稳定区向高层发展时,伴随着强上升运动作为触发机制,引发条件性对称不稳定能量的释放,产生更多的对流降水,使模拟的总降水量与实况更加一致.条件性对称不稳定的发展加强与降雪强度、辐合辐散和上升运动变化一致,条件性对称不稳定是造成暴雪发展加强的主要机制之一.通过对"海棠"台风72 h的模拟表明,条件性对称不稳定主要发生在台风的低层,且其水平分布呈螺旋状结构.条件性对称不稳定效应对台风路径影响较小,但对台风强度影响较大,在模式中考虑垂直-倾斜对流一体化参数化方案后,与仅考虑垂直积云对流参数化方案相比,72 h模拟的平均台风中心最低气压降低了3 hPa,最大达8 hPa.在模式中考虑条件性对称不稳定的影响,可使模式台风中上层的暖心结构更加明显,上升运动和对流性降水增强,对流释放的更多凝结潜热使台风得到进一步加强.     

Urban scale modelling: The Sydney hailstorm of 14 April 1999

Summary The aim of this study is to point to the very large improvements that are taking place in a range of modelling applications in the urban areas. The particular phenomenon chosen in this study is a supercell, but it could well have been any other aspect of urban modelling. The Sydney hailstorm of 14 April 1999 was a long-lived, high precipitation supercell that produced a massive damage bill of over 2 billion Australian dollars from its hail swath. The Sydney hailstorm was poorly forecast for a number of reasons including: the severe weather season had officially ended so there were no specialist staff on duty when the hailstorm struck Sydney; the storm proved very difficult to predict and it was expected to continue heading out to sea; and the forecast guidance from all available operational numerical models was inadequate at the resolution required for a supercell simulation. Here, our interest is on the last of the problems, namely, the quality of the operational numerical model guidance, especially given the impact it had on a densely populated urban region.In this study, we compare the numerical guidance available at the time with current modeling capability which, although in research mode at present, will soon be available in real-time mode. The operational models were hydrostatic models run at horizontal resolutions of 25km at best, compared with 1km horizontal resolution for the non-hydrostatic research model. The research model also had a high-order differencing scheme and a sophisticated six phase cloud physics scheme compared with the much simpler parameterized convection in the operational models. The operational model produced very little convective precipitation and it was displaced well to the north of Sydney. The research model generated a supercell with a track and a hail size distribution that was encouragingly close to the observed.     

两次强冰雹超级单体风暴双偏振特征对比

利用S波段双偏振天气雷达资料、探空和地面常规气象观测资料及灾情调查, 对2020年6月25日河北省蠡县和2021年7月9日山东省章丘的两次特大冰雹超级单体风暴双偏振特征进行对比。结果表明:两次超级单体风暴均发生在西北气流形势下, 章丘风暴具有较强的对流有效位能、较大的湿度和较高的湿球0℃层高度。蠡县风暴强度明显大于章丘风暴, 但差分反射率柱和比差分相移柱高度明显低于章丘风暴。蠡县风暴弱回波区上方存在深厚的强度超过65 dBZ强回波悬垂, 即悬垂的冰粒子循环增长产生较大的冰雹粒子, 大的冰雹粒子进入下降通道后, 再次产生明显增长且更加不规则, 导致更强的水平极化反射率因子和更小的相关系数。湿度的垂直分布是风暴发展强度的关键环境因素之一。蠡县超级单体风暴的产生环境非常干, 章丘超级单体风暴的产生环境相对较湿。     

华南中尺度暴雨数值预报的不确定性与集合预报试验

利用非静力MM5模式,分析了不同积云对流参数化方案对华南暖区暴雨数值预报的不确定性影响,进行了中尺度暴雨模式扰动集合预报试验。不同对流参数化方案的对流凝结加热引起不同的局地温度扰动,通过大气内部的热力动力过程,导致垂直速度的差异,进而影响网格尺度和次网格尺度降水时间、地点和强度。后续降水再通过凝结潜热释放形成新的扰动源。不同积云对流参数化方案还可引起扰动源能量传播方式不同,最终使模拟大气的动力和热力结构有差异。针对物理过程的不确定性,使用两种模式扰动方法构造集合预报扰动模式,第一种方法是随机组合不同积云对流参数化方案和边界层方案,第二种方法是扰动Grell积云对流参数化方案中主要参数振幅。集合预报结果表明,第一种方法的集合预报效果优于第二种方法,仅扰动参数振幅值似乎还不足以反映华南暴雨预报的不确定性。单一的确定性预报在暴雨落区和强度方面的可信度不稳定,集合产品能给华南暴雨过程提供更有用价值的指导预报,具有较高的应用价值。     

藏北高原地表能量和边界层结构的数值模拟

利用耦合了NCAR LSM陆面过程的中尺度模式MM5V3.7和2002年8月CAMP/Tibet加强期的观测资料,对藏北高原地区地气交换过程进行了48 h模拟研究。模式较好地模拟了该地区的山谷风环流;并将模拟的地表通量在中尺度区域上与NCEP/NCAR全球大气再分析格点资料(NNRP)获得的结果进行了比较,同时也与单站的实测值进行了比较,结果显示:模拟的地表通量与NNRP得到的结果比较吻合,同时可以得到雨季时藏北、藏东地区潜热通量大于感热通量,而高原西部感热通量大于潜热通量,这与观测试验分析结果一致;与单站试验结果比较,模拟的感热通量与实测值一致,潜热通量的模拟值和实测值有一定差别。模拟的边界层位温廓线与实测值比较,模式模拟的对流混合层和夜间残留层都与实测结果吻合,但模拟的混合层高度较实测值高。由此来看,中尺度模式MM5V3.7能够较好地模拟藏北高原的地表能量和边界层结构特征,但还需要进一步完善陆面过程和物理过程参数化方案。     

不同云降水方案对一次登陆台风的降水模拟

利用美国最新发展的新一代中尺度数值预报模式(WRF),研究了不同云降水方案对一次登陆台风的降水模拟问题,结果表明,在网格距适当小的情况下,同时采用积云对流参数化方案和云微物理方案(Kessler方案)时,其降水预报优于只使用积云对流参数化方案时的预报;在台风降水模拟的初期,Kain-Fritsch方案比Betts-Miller方案产生降水更快,更接近实况降水;在台风登陆后随时间的延长,对流降水重要性逐步下降,网格尺度降水逐渐增强。     

Development of an Integrated Vertical-Slantwise Convective Parameterization Scheme and Its Associated Numerical Experiments

An integrated vertical-slantwise convective parameterization scheme,based on the vertical Kuo-Anthes and the slantwise Nordeng convective parameterization schemes,is introduced into the MM5 model.By employing the MM5 model with the proposed scheme,numerical simulations of a snowstorm event that occurred over southern China on 28 29 January 2008 and of Typhoon Haitang (2005) are conducted.The results indicate that during the snowstorm event,the atmosphere was convectively stable in the vertical direction but with conditional symmetric instability (CSI) in the lower troposphere,and when the area of CSI developed and extended to upper levels,strong rising motion occurred and triggered the release of large amount of energy,producing enhanced convective precipitation with the total precipitation much closer to the observation.The development and strengthening of CSI corresponded to changes in the intensity of snowfall,convergence,and ascending motions of air,revealing that CSI was responsible for the initiation and growth of the snowstorm.The results from a 72-h explicit simulation of Typhoon Haitang indicate that CSI occurred mainly at lower levels with a well-defined spiral structure,and it tended to have a larger impact on the intensity of typhoon than on its track.The minimum pressure at the typhoon center for the 72-h runs with the integrated vertical-slantwise convective parameterization scheme was on average 3 hPa (maximum 8 hPa) lower than that from the runs with only the vertical cumulus parameterization scheme.Introducing the influence of CSI into the model has improved the warm core structure at the middle and upper levels of the typhoon,with stronger and persistent upward motions causing increased precipitation,and the latent heat released through convection in turn made the typhoon develop further.     

对流层中层中尺度涡旋在台风榴莲(2001)生成中的作用—数值模拟及验证

西太平洋热带气旋（TC）的生成和季风槽及中尺度对流系统（MCS）的活动有密切关系,但以往这方面的实例数值模拟很少。为了进一步探讨由MCS对流强迫产生的对流层中层中尺度涡旋（MCV）在TC生成中的作用,作者利用非静力平衡的中尺度模式PSU/NCAR MM5对台风榴莲（2001）的生成过程进行了高分辨率（6 km）数值模拟和比较验证。结果表明:模式成功地模拟出榴莲的生成地点,其与MCS的相对位置关系与以往的观测研究结果一致;模拟的TC移动路径、强度变化与最优观测报告比较接近,准确反映了TC未来登陆地点,以及中心气压缓慢下降和迅速下降两个阶段;对云系演变的模拟,成功模拟出了TC初生时的涡旋云系和季风槽中MCS云系的分离现象,以及在TC登陆前达到成熟阶段时出现的台风眼和螺旋云带。此外,模式还成功模拟出中层MCV,它的水平尺度约200 km,位于800~400 hPa之间,具有暖心结构等,均与已有观测结果相近。模式初始场中包含有充分的MCS信息,是模拟取得成功的关键因素之一。     

尺度适应的对流参数化方案对一次华南飑线模拟的影响

尺度适应(scale-aware)的物理过程是现代数值预报发展的一种趋势,本文针对GRAPES_Meso(Global/Regional Assimilation and Prediction Enhanced System)模式没有考虑物理过程尺度适应的现状,首先在KFeta(Kain-Fritsch Eta)积云对流参数化方案中引进了尺度适应过程,对该方案的对流时间尺度、格点垂直速度以及夹卷率进行了基于尺度适应参数化的改进。为研究尺度适应KFeta方案与原KFeta方案对不同分辨率模式模拟结果的影响,选取了一次华南飑线过程进行数值模拟和影响分析。结果表明:在3 km、5 km、10 km、20 km水平分辨率的GRAPES_Meso模式中,尺度适应KFeta方案相比原方案,对降水强度及落区分布的模拟有一定的正效果,随着模式水平分辨率提高,次网格降水减少、格点降水增多、对流层中低层夹卷略有增强,对原来存在的对流层高层及低层偏冷的偏差有一定的改进,对流活跃区域的上升气流强度、云中水凝物含量更符合真实的天气系统演变。综合来看,改进后的方案更适用于高分辨率数值预报模式,该研究结果可以为尺度适应对流参数化方案的应用及数值模式强降水预报性能的优化提供有益的参考。     

The Uncertainty of Mesoscale Numerical Prediction of Heavy Rain in South China and the Ensemble Simulations

1. IntroductionHeavy rain is a kind of severe natural calamitythat influences South China. After decades of years oftests and theoretical exploration by Chinese scientists,significant progresses have been achieved in its predic-tion and basic theoretical studies (Huang, 1986; Xue,1999; Zhou et al., 2003). Currently, the mesoscale nu-merical model has already been employed as one of themajor tools in the prediction and research on heavyrain in South China, promoting considerably the ac-curac…     

Development of an integrated vertical-slantwise convective parameterization scheme and its associated numerical experiments

An integrated vertical-slantwise convective parameterization scheme, based on the vertical Kuo-Anthes and the slantwise Nordeng convective parameterization schemes, is introduced into the MM5 model. By employing the MM5 model with the proposed scheme, numerical simulations of a snowstorm event that occurred over southern China on 28-29 January 2008 and of Typhoon Haitang (2005) are conducted. The results indicate that during the snowstorm event, the atmosphere was convectively stable in the vertical directi...     

Revisiting the latent heat nudging scheme for the rainfall assimilation of a simulated convective storm

Summary Next-generation, operational, high-resolution numerical weather prediction models require economical assimilation schemes for radar data. In the present study we evaluate and characterise the latent heat nudging (LHN) rainfall assimilation scheme within a meso-γ scale NWP model in the framework of identical twin simulations of an idealised supercell storm. Consideration is given to the model’s dynamical response to the forcing as well as to the sensitivity of the LHN scheme to uncertainty in the observations and the environment. The results indicate that the LHN scheme is well able to capture the dynamical structure and the right rainfall amount of the storm in a perfect environment. This holds true even in degraded environments but a number of important issues arise. In particular, changes in the low-level humidity field are found to affect mainly the precipitation amplitude during the assimilation with a fast adaptation of the storm to the system dynamics determined by the environment during the free forecast. A constant bias in the environmental wind field, on the other hand, has the potential to render a successful assimilation with the LHN scheme difficult, as the velocity of the forcing is not consistent with the system propagation speed determined by the wind. If the rainfall forcing moves too fast, the system propagation is supported and the assimilated storm and forecasts initialised therefrom develop properly. A too slow forcing, on the other hand, can decelerate the system and eventually disturb the system dynamics by decoupling the low-level moisture inflow from the main updrafts during the assimilation. This distortion is sustained in the free forecast. It has further been found that a sufficient temporal resolution of the rainfall input is crucial for the successful assimilation of a fast moving, coherent convective storm and that the LHN scheme, when applied to a convective storm, appears to necessitate a careful tuning.     

Meridional propagation of large-scale monsoon convective zones

Summary Observational studies indicate that the convective activity of the monsoon systems undergo intraseasonal variations with multi-week time scales. The zone of maximum monsoon convection exhibits substantial transient behavior with successive propagating from the North Indian Ocean to the heated continent. Over South Asia the zone achieves its maximum intensity. These propagations may extend over 3000 km in latitude and perhaps twice the distance in longitude and remain as coherent entities for periods greater than 2–3 weeks. Attempts to explain this phenomena using simple ocean-atmosphere models of the monsoon system had concluded that the interactive ground hydrology so modifies the total heating of the atmosphere that a steady state solution is not possible, thus promoting lateral propagation. That is, the ground hydrology forces the total heating of the atmosphere and the vertical velocity to be slightly out of phase, causing a migration of the convection towards the region of maximum heating. Whereas the lateral scale of the variations produced by the Webster (1983) model were essentially correct, they occurred at twice the frequency of the observed events and were formed near the coastal margin, rather than over the ocean.Webster's (1983) model used to pose the theories was deficient in a number of aspects. Particularly, both the ground moisture content and the thermal inertia of the model were severely underestimated. At the same time, the sea surface temperatures produced by the model between the equator and the model's land-sea boundary were far too cool. Both the atmosphere and the ocean model were modified to include a better hydrological cycle and ocean structure. The convective events produced by the modified model possessed the observed frequency and were generated well south of the coastline.The improved simulation of monsoon variability allowed the hydrological cycle feedback to be generalized. It was found that monsoon variability was constrained to lie within the bounds of a positive gradient of aconvective intensity potential (I). The function depends primarily on the surface temperature, the availability of moisture and the stability of the lower atmosphere which varies very slowly on the time scale of months. The oscillations of the monsoon perturb the mean convective intensity potential causing local enhancements of the gradient. These perturbations are caused by the hydrological feedbacks, discussed above, or by the modification of the air-sea fluxes caused by variations of the low level wind during convective events. The final result is the slow northward propagation of convection within an even slower convective regime. The ECMWF analyses show very similar behavior of the convective intensity potential. Although it is considered premature to use the model to conduct simulations of the African monsoon system, the ECMWF analysis indicates similar behavior in the convective intensity potential suggesting, at least, that the same processes control the low frequency structure of the African monsoon. The implications of the hypotheses on numerical weather prediction of monsoon phenomenon are discussed.