非静力AREM模式设计及其数值模拟II:数值模拟试验

在程锐等（2018）中,我们完成了非静力AREM（Advanced Regional Eta-coordinate Model）模式动力框架设计。本文将通过理想和实例试验检验其模拟能力。设计理想试验并通过与国际成熟的中尺度非静力框架比较,直接检验非静力AREM三维动力框架在细致分辨率（约1 km）下的模拟性能。可以看出,非静力AREM与ARPS（Advanced Regional Prediction System）、WRF（Weather Research and Forecasting Model）模拟出类似的积云对流结构及演变特征,从而基本确证了本文发展的非静力框架的正确性。结合原静力平衡模式的初始化和物理参数化过程,形成非静力AREM模式系统。台风实例模拟表明,粗分辨率下静力、非静力AREM模式性能接近;但在高分辨率下,非静力明显优于静力模式。我们还开展了批量降水试验检验,对非静力AREM模式性能进行了进一步的验证。     

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.     

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

尺度适应(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方案相比原方案,对降水强度及落区分布的模拟有一定的正效果,随着模式水平分辨率提高,次网格降水减少、格点降水增多、对流层中低层夹卷略有增强,对原来存在的对流层高层及低层偏冷的偏差有一定的改进,对流活跃区域的上升气流强度、云中水凝物含量更符合真实的天气系统演变。综合来看,改进后的方案更适用于高分辨率数值预报模式,该研究结果可以为尺度适应对流参数化方案的应用及数值模式强降水预报性能的优化提供有益的参考。     

Application of the hybrid convection parameterization scheme HYMACS to different meteorological situations

Convection can cause severe precipitation events and is thus of major interest for quantitative precipitation forecasts. Due to its small horizontal scale convection usually has to be parameterized in numerical weather prediction (NWP) models as a subgrid scale phenomenon. Classical mass flux convection schemes assume grid box sizes much larger than the scale of the convective circulation. Thus, the convective mass transport is closed in the local grid column and no net mass transport occurs on the grid scale.In contemporary NWP models with grid sizes of a few kilometers, where convection is already partially resolved, the classical approach leads to a conceptual problem. This can be overcome by a hybrid mass flux convection scheme (HYMACS), in which only the small scale convective updrafts and downdrafts are parameterized, whereas the treatment of the larger scale environmental subsidence is left to the grid scale equations. Different to the classical schemes, HYMACS produces a net convective mass flux exerting pressure gradient forces to the grid scale model.We discuss real cases of different meteorological situations simulated with the operational weather forecast model COSMO of the German Meteorological Service as the hosting model of HYMACS. Precipitation and cloud top pressures simulated with HYMACS are compared with results from classical convection schemes and observational data from the DWD station network and the SEVIRI (MSG) satellite instrument.     

Evaluation of cloud and precipitation parameterization using a single-column model: A TWP-ICE case study

Cloud and precipitation parameterization schemes are evaluated, and their sensitivity to the method and/or parameters used to determine cloud physical processes is examined using a singlecolumn version of the Unified Model (SCUM). In the experiment for TWP-ICE, cloud fraction is overestimated (underestimated) in the upper (lower) troposphere due to the wet (dry) bias. The precipitation rate is well simulated during the active monsoon period, but overestimated during the suppressed monsoon and clear skies periods. In the moist convection scheme, trigger condition and entrainment process affect the lower tropospheric humidity through the impact on convective occurrence frequency and intensity, respectively. Strengthening the trigger condition and using the adaptive entrainment method alleviate the low-level dry bias. In the microphysics scheme, more large-scale precipitation is produced with prognostic rain, due to rain sedimentation considering vertical velocity of rain drop, than with diagnostic rain. Less ice/snow deposition with the prognostic two-ice category results in lower ice water content and upper-level cloud fraction than with the diagnostic splitting method for the twoice category. In the cloud macrophysics scheme, the prognostic cloud fraction and cloud/ice water content scheme produces a larger cloud fraction and more cloud/ice water content than the diagnostic scheme, mainly due to detrainment from moist convection (cloud source) that surpasses the effect of convective heating and drying (cloud sink). This affects temperature by influencing the radiative, convective, and microphysical processes. The experiment with combined modifications in cloud and precipitation schemes shows that interaction between modified moist convection and cloud macrophysics schemes results in more alleviation of the cold bias not only at the lower levels but also at the upper levels.     

Diagnostic analyses of convective events — The effect of propagating gust fronts

This paper is devoted to the testing experiment of a postprocessing tool aimed at the objective analysis of propagating gust fronts in a given convective environment. The tool is being developed to be applicable in the operational mode by utilizing NWP model outputs. The experiment was carried out on two summer convective cases which occurred in the Czech Republic. The cases were numerically simulated by the limited area NWP model LM COSMO with the horizontal resolution of 2.8 km. They represent different types of convective systems, both accompanied by objectively identifiable gust fronts and causing heavy precipitation. The event from July 2000 was characterised by the development of isolated thunderstorms. The other event from July 1998 was a long-lasting and organised convective system — a squall line. The hypothesis was that the developed postprocessing tool is capable to evaluate the role which downdraft outflows played in the decay and initiation of convective cells by interaction with convective environment and thus in prolongation of convection lifetime.The procedures of the Objective Analysis of Gust Fronts (OAGF) were applied to the thermodynamic outputs of the LM COSMO. The aim was to determine the position of gust fronts within the domain and to assess their speed of movement and the potential to initiate convection according to the properties of ambient vertical shear and stability as well as humidity conditions ahead of the respective downdraft outflows. In addition, the Radar Simulation Model (RSM) was employed to monitor the simulated convective systems in arbitrary PPI and RHI scans and to verify qualitatively the LM COSMO precipitation forecasts.The case studies has confirmed the applicability of the LM COSMO–OAGF chain and RSM, which may represent the potential for improving the operational nowcasting of hazardous convection phenomena. In both simulations, the objective gust fronts moved on into the vertical shear-favourable environment for triggering new convection. In addition to the dynamical organisation, there were also favourable stability and humidity conditions in the area of forced upward motions in the simulation of the event from July 2000.     

Evaluating the impact of lightning data assimilation on mesoscale model simulations of a flash flood inducing storm

Flash floods are associated with highly localized convective storms producing heavy rainfall. Quantitative precipitation forecasting of such storms will potentially benefit from explicit representations of deep moist convection in numerical weather prediction models. However, explicit representation of moist convection is still not viable in operational mesoscale models, which rely on convective parameterizations for issuing short to medium-range forecasts. In this study we evaluate a technique that uses regional Cloud-to-Ground (CG) lightning observations to define areas of deep moist convection in thunderstorms and adjust the model-generated precipitation fields in those regions. The study focuses on a major flash flood inducing storm in central Europe (23 August 2005) that was simulated with the aid of an operational weather forecasting system (POSEIDON system based on Eta/NCEP model). The performance of the technique is assessed using as reference distributed rainfall estimates from a network of radar observations. The results indicate that CG lightning data can offer sufficient information to increase the mesoscale model skill in reproducing local convective precipitation that leads to flash floods. The model error correction is shown to be proportional to the density of lightning occurrence, making the technique potentially suitable for operational forecasting of flash flood inducing thunderstorms.     

初始场和分辨率对WRF模拟对流性降水的影响

文章利用WRF模式分析了2015年7月29日内蒙古中部地区一次典型冷涡过程带来的对流性降水,探讨了初始场和网格分辨率对对流性降水模拟结果的影响。初始场和网格分辨率均会影响模式模拟结果,对于对流性短时强降水而言,在相同的初始场条件下,网格分辨率越高,对降水强度和落区的预报越准确,而在相同的高分辨率网格条件下,初始场的精度越高,模拟不同等级降水的TS评分均较高,表明在高分辨率网格和高精度的初始场条件下,模式对不同等级的降水均有较好的模拟。     

Synoptic and mesoscale study of a severe convective outbreak with the nonhydrostatic Global Environmental Multiscale (GEM) model

Summary ?A nonhydrostatic 4-km version of the Global Environmental Multiscale (GEM) model, with detailed microphysics included, was used to forecast the initiation, development, and structure of a tornado-producing supercell storm that occurred near Pine Lake (Alberta, Canada) on 15 July 2000. Examination of observations and comparison with conceptual models indicate that this storm is a good example of supercell storms that regularly produce summertime severe weather over Alberta. It was found that the high-resolution model was able to reproduce the early initiation of convective activity along the Rocky Mountains foothills, as well as the rapid northeastward propagation towards the Pine Lake area and the subsequent intensification into a supercell storm. The general structures of the forecasted convective system correspond well with conceptual representations of such events. Overall, this high-resolution forecast of the Pine Lake supercell storm was a significant improvement over the current operational version of the GEM model (24 km), which was not able to intensify the foothills’ convection into a supercell storm. Finally, it was found that the nonhydrostatic version of the model produces better trajectory and propagation speed of the convective system, as compared with the hydrostatic one. Received March 20, 2001; revised August 24, 2001     

2001年8月23日华北强风暴动力机制的数值研究

利用非静力平衡模式(MM5 V3.5)对2001年8月23日发生在华北地区的一次强对流风暴过程进行了数值模拟,在取得合理模拟结果的基础上,着重分析了风暴发生的热力条件、太行山地形动力作用和风暴中尺度结构特征,并依据倾斜涡度发展理论,分析了强风暴的发展机制。结果表明:本次风暴是在以西北气流为主导的环流背景下产生的,低层增温,中高层降温和整层增湿是造成层结不稳定的重要原因。上游地区对流引发的干冷下沉气流沿太行山背风面下滑形成的下坡风是触发强风暴发生的直接动力机制。湿位涡分析表明,风暴发生区具有明显的等熵面倾斜,对流层中低层pm1<0区域,同时pm2>0,满足倾斜涡度发展的条件,对强风暴的发展具有一定指示性。风暴发展早期,其垂直方向次级环流可能与新雷暴的产生和雷暴的跳跃式传播有关。太行山地形引起的边界层风场的变化,包括下坡气流和边界层中尺度辐合线对风暴的触发、组织和移动发挥着重要作用。     

A Modeling Study on the Mechanism of Convective Initiation of a Squall Line over Northeastern China

High-resolution numerical simulation results of a squall line initiated along a convergence zone in northeast China on 26 June 2014 were presented in this study. The simulation was performed by a convection-permitting model with coarse and fine grids of 4 and 1.33 km, respectively, and the simulation results were validated against the observation. Results showed that the simulation adequately reproduced the life cycle of the squall line, which allowed detailed investigation of the mechanism of convective initiation in this case. The synoptic condition was favorable for convective initiation and the convection was triggered in a convergence zone, where a branch of dry and cold air and a branch of moist and warm air collided. The water vapor flux divergence was inhomogeneous and some cores of water vapor convergence existed in the convergence zone. These cores were the spots where water vapor converged intensely and the air there was forced to rise, creating favorable spots where the convection was initially triggered. A series of quasi-equally spaced vortices near the surface, which themselves were the result of horizontal shear instability, were accountable for the inhomogeneity of the surface water vapor flux divergence. These vortices rotated the moist air into their north and dry air into their south, thus creating more favorable spots for convective initiation in their north. After initiation, the updraft turned the horizontal vorticity into vertical vorticity in the mid-level. The vortices near the surface collaborated with the vorticity maxima in the mid-level and enhanced the development of convection by providing water vapor.     

积云参数化方案对热带降水年循环模态模拟的影响

本文利用中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室( LASG)发展的大气环流模式(SAMIL),采用Zhang-McFarlane (ZM)和Tiedtke (TDK)两种积云对流参数化方案,讨论了积云对流参数化方案对热带降水年循环模态模拟的影响.结果表明,两种积云对流参数化方案均能合理再现...     

苏北一次强降水超级单体风暴过程的诊断分析

利用常规观测资料、NCEP再分析资料、FY2C卫星和多普勒雷达资料,对2008年7月22日发生在苏北的一次强降水超级单体风暴过程进行诊断分析。天气分析显示,风暴发生于高湿、较低的抬升凝结高度、强对流不稳定(3 445 J/kg)和中到强的垂直风切变(0~6 km,18 m/s)环境,这种大气环境非常有利于强降水超级单体风暴的发生发展。雷达回波分析揭示,该超级单体的演化可归结为"孤立单体—经典强降水超级单体—减弱东移"三个阶段,持续时间超过2 h。强降水超级单体风暴成熟期,呈现出典型的倒"V"型缺口、中低层有界弱回波区和反射率因子大值区由低层向高层往低层入流一侧倾斜的特征,相应的雷达径向速度场显示在倒"V"型缺口附近的强降水区中存在一个成熟的中气旋。湿位涡的诊断结果表明:高层干冷空气侵入触发潜在对流不稳定能量释放,有利于对流运动的发展;中低层大气对流不稳定与条件对称不稳定共存,既有垂直对流,又有倾斜对流发生,同时边界层的偏东风入流向暴雨区提供充沛的水汽,对暴雨的发生发展起增幅作用。     

The effects of precipitation schemes and horizontal resolution on the major rainband in typhoon Flo (1990) predicted by the MRI mesoscale nonhydrostatic model

Summary ?This paper describes a numerical study of the major spiral rainband in typhoon Flo (1990) using the Meteorological Research Institute Mesoscale Nonhydrostatic Model (MRI-NHM). The effects of precipitation schemes and horizontal resolution on the representation of the simulated rainband are discussed. Dynamic and thermodynamic structures of the simulated major rainband to the north of the storm center are well represented in the model with a 5 km horizontal resolution. The structures are consistent with observational results reported for other tropical cyclones. Among the realistic features are: a cold pool and convergence on the inner side of the band; convergence above low-level inflow layers; and the outward slope of the updraft with height. The band is caused by the motion of the storm through its surroundings where horizontal wind has vertical shear. The simulation of the structure and precipitation pattern associated with the major rainband depends on the precipitation scheme rather than the horizontal resolution. The band appears more realistic when using explicit cloud microphysics as a precipitation scheme, rather than moist convective adjustment. This result is attributable to the difference in scheme triggering. In the simulation with moist convective adjustment, the elimination of vertical instability in low-level atmosphere is excessive, suppressing band formation. The overall structure of the band is also more realistic in the simulation using explicit cloud microphysics, because a cold pool exists in the lower layers and the vertical axis of upward flow tilts outward. This result suggests that prediction will partly depend on variables associated with cloud microphysics, such as the mixing ratio of cloud water. The horizontal grid distance, which varied between 5 and 20 km, quantitatively influenced the rainfall amount, although the large-scale band structure remained unchanged. The rainfall amount increased as the grid interval was reduced from 20 to 10-km but decreased as the interval was further reduced from 10 to 5 km. Received March 20, 2001; revised August 20, 2001     

Modification of a Parametrization of Shallow Convection in the Grey Zone Using a Mesoscale Model

In current operational numerical weather prediction models, the effect of shallow convection is parametrized. The grey zone of shallow convection is found between the horizontal resolutions of mesoscale numerical models (2–3 km) and large-eddy simulations (10–100 m or finer). At these horizontal scales the shallow convection is to some extent explicitly resolved by the model. The shallow-convection parametrization is still needed, but has to be regulated according to the model horizontal resolution. Here the behaviour of the non-hydrostatic mesoscale numerical weather prediction model Application of Research to Operations at Mesoscale is examined in the grey zone and a new scale-adaptive surface closure of its shallow-convection parametrization, dependent on horizontal resolution, is defined based on large-eddy simulation. The new closure is tested on a series of numerical experiments and validated on a 15-day-long real case period. Its impact on the development of deep convection is examined in detail. The idealized simulations show promising results, as the mean profiles of the subgrid and resolved turbulence change in the desired way. Based on the real case tests our modification has a low impact on model performance, but is part of a set of upgrades of the current parametrization that is aimed to treat the shallow convection grey zone.     

A method of forecasting the heavy convective precipitation zones during the cold season

A method of forecasting the heavy convective precipitation zones is based on the computation of the Q-vector divergence and moist equivalent-potential vorticity in the layer from 925 to 700 hPa. The prognostic fields of meteorological parameters from NWP models in GRIB or GRIB2 codes are used for the computation. Presented are the results of the statistical analysis. The coefficients of linear regression are computed and the working equation for forecasting heavy convective precipitation zones during the cold season is derived. Analyzed are some heavy snowfalls in the central region of the European part of Russia.     

The prediction of two cases of severe convection: implications for forecast guidance

Summary Severe convection in the Sydney basin regularly produces destructive winds, heavy rain or flash flooding, and damage from large hail. Such convective activity is a major forecasting challenge for the Sydney basin, especially during the October to April warm season. There presently is a need to provide timely, accurate and reliable numerical guidance to supplement the current probabilistic convective outlooks, issued by the operational forecasters. Initial work has been carried out that examines two cases of severe convection in the Sydney basin. The performance of a very high resolution (2 km) numerical weather prediction (NWP) model is assessed in terms of how well it performed in providing guidance on heavy rainfall and hail, as well as other key mesoscale phenomena such as low level convergence lines. The model results discriminated well between severe convection that actually did occur in the first case (1 December, 2000) and the failure of severe convection to develop in the second case (8 December, 2000). The operational forecasters predicted severe convection to occur in both cases. It is now planned to take the next step of augmenting the quasi-subjective approach of the operational forecasters with the NWP guidance to provide an enhanced capacity to forecast severe convection in the Sydney basin. Received October 20, 2001 Revised December 28, 2001     

雷暴与强对流临近天气预报技术进展

临近预报指0—6h(0—2h为重点)的高时空分辨率的天气预报,预报对象是该时段内出现明显变化的天气现象,主要包括雷暴、强对流、降水、冬季暴风雪、冻雨、沙尘暴、低能见度(雾)、天空云量等,其中,以雷暴和强对流天气的临近预报最具挑战性。综述了针对雷暴和强对流天气的以主观预报为主、结合客观算法的临近预报技术,同时讨论了高分辨率数值预报模式在临近预报中的应用。主观临近预报技术包括基于多普勒天气雷达观测数据并结合其他资料(常规高空和地面观测、气象卫星云图、快速同化循环的数值预报产品等)对雷暴生成、发展和衰减,特别是对强对流天气(包括强冰雹、龙卷、雷暴大风和对流性暴雨)的临近预报,客观算法包括几种应用最广的雷达回波或云图外推算法和强对流天气识别技术。高分辨率数值预报模式的应用包括与雷达回波外推融合延长临近预报时效,与各种观测资料融合得到快速更新的三维格点资料为雷暴和强对流近风暴环境的判断提供重要参考。     

初边界条件不确定性对AREM模拟一次华南致洪暴雨的影响

利用AREM（advanced regional E-gridη-coordinates model）模式,针对2005年6月21日发生在华南的一次特大致洪暴雨过程,研究了模式初始场和边界场的不确定性对AREM模拟大暴雨过程的影响。研究表明：模式初始场和边界场对模式模拟降水的不同时段影响存在明显差异,初始扰动误差越大模式误差也越大;误差增长先在中小尺度内伴随着湿对流不稳定,且增长极其迅速,接着向大尺度传播,由于对流有效位能的逐渐释放,大气不稳定度降低,误差在大尺度上增长缓慢;在初始场和边界场相同精度的情况下,增加边界场的中尺度信息,尤其是400 km以下尺度的信息,比增加初始场的中尺度信息更能有效抑制误差的增长。     

基于4种高分辨率模式的辽宁省雷达回波预报邻域检验结果对比

选取2017—2018年6—9月辽宁省不同降水性质, 具有2种不同特征的20次天气过程个例, 应用模糊检验邻域法中的分数技巧评分(Fraction Skill Score, FSS), 评估华东模式、华北模式、GRAPES_3km模式和睿图东北模式对辽宁省中小尺度系统的预报能力。结果表明: 区域性降水过程和局地性降水过程雷达回波强度越小, 邻域半径越大, 高分辨率模式预报技巧越高。当雷达回波大于30 dBz时, 各高分辨率模式对局地性降水的雷达回波预报FSS评分均较高。当邻域半径为3 km时, 区域性降水过程中, 华北模式预报技巧在各级别雷达回波预报中均高于其他模式, 最大FSS差值为0.031。局地性降水过程中, 华东模式预报效果较好, 最大FSS评分为0.127, 表明华东模式预报中小尺度对流系统能力更强。局地性降水过程, 睿图东北模式在08—23时预报时次中, “中间”时次的预报效果优于“两头”时次的预报, 两个时次最大FSS差值为0.121。