中尺度大气数值模式发展现状和应用前景

对国内外当前一些先进的中尺度大气数值模式的发展现状,应用前景及发展趋势作了概要综述。其内容包括：模式动力学的改进,叫就度模拟系统特征,区域谱模 发展积云参数化和显式云物理方案,行星边界层参数化,大气辐射参数、四维资料同化,区域实时数值天气预报,中尺度数值天气预报应用前景及新一代中尺度模式发展趋势。     

中尺度数值模拟中的边界层多尺度湍流参数化方案

该文在多尺度湍流理论的研究成果基础上, 将边界层湍流风谱与平均量的梯度相联系, 建立了边界层多尺度湍流参数化子模式, 之后放入MM5模式中进行了个例模拟研究, 并与MM5模式附带的M RF边界层参数化、Blackadar高分辨率边界层参数化的模拟结果进行比较和分析。结果表明, 多尺度湍流理论能够反映出实际大气边界层中热量垂直输送的规律, 将其用于中尺度数值预报模式的边界层物理过程参数化是可行的; 多尺度湍流参数化在地表层和边界层内各个层次上都着重考虑含能量最大的涡的作用以及水平热力不均匀性的影响, 因此在地形和下垫面比较复杂的区域, 对中尺度天气系统的模拟有进一步发展的前景。     

The influences of boundary layer parameterization schemes on mesoscale heavy rain system

1.Introducti0nThemesoscaleoperationalmodel-whichisoftenused,isMM4orMM5'butMM4isusedfrequentlyonlO3kmscale.Thephysicalprocessesinthismodeldevelopconstantly.FororiginalMM4,thecomputationofsurfacefluxesisnotaccurate,andKmodelfortheturbu-lencefluxesbetweenany2levelsneedstobeimprovedbynewtreatment.Inordertostudytheinfluencesofboundarylayerparameterizationschemesonmesoscaleheavyrainsystem,sur-facefluxesandKmodelinoriginalMM4areimprovedbytherecentresearchinthispaPer.Theflux-profilerelationsforv…     

不均匀地表产生的中尺度通量的数值试验

采用北京大学三维的复杂地形中尺度模式,结合陆面过程模式(SiB),模拟了草原和沙漠并存的下垫面的边界层大气运动.利用SiB模式计算了地表辐射、感热、潜热通量,并且预报地表温度.中尺度模式则模拟了沙漠地区受热抬升,形成的辐合运动,垂直速度的分布,不同高度上水平流场的变化以及中尺度动量和热量通量,把中尺度通量跟湍流通量进行了比较,以确定这种中尺度运动在GCM模式的参数化过程中的重要性.试验表明中尺度通量尤其是热量通量要比湍流通量大很多.     

物理过程参数化方案对中尺度暴雨数值模拟影响的研究

利用中尺度非静力MM 5模式和中国 2 0 0 1年 8月的 4个暴雨个例 ,研究了非绝热物理过程对中国暴雨动力和热力场预报的影响 ,深入分析了对流参数化方案在中尺度暴雨预报中的作用 ,讨论了利用模式扰动方法开展中国暴雨集合预报的可行性。结果表明 ,在短期数值预报中 ,非绝热物理过程对高度场预报影响较小 ,但边界层方案和对流参数化方案对产生暴雨的 3个基本条件即水汽通量散度、垂直速度、不稳定层结的影响很明显。不同对流参数化方案所预报的中尺度热力、动力场离差的结构特征与所预报降水的离差特征相似 ,且主要是在模式积分初期迅速增加 ,其后即趋于稳定。对中国热力场较均匀的暴雨过程 ,可以通过扰动模式的边界层和对流参数化方案 ,构造集合预报模式     

Mesoscale Vegetation-Breeze Circulations and Their Impact on Boundary Layer Structures at Night

The impact of well watered mesoscale wheat over mid-latitude arid areas on mesoscale boundary layer structures (MBLS) and climate has been investigated in the study .using a mesoscale biophysical, meteorological model (BM) developed in the current study. The BM is composed of six modules:mesoscale atmospheric module, soil module, vegetation module, snow-atmosphere interaction module, underlying surface meteorology module and subgrid scale flux parameterization module. The six modules constitute an interacting system by supplying boundary conditions to each other.The investigation indicates that a horizontal pressure gradient associated with mesoscale perturbations in temperature and humidity is created during the day, which results from more water transpired from the vegetation canopy (VC) and evaporated from underlying wet soil. Non-classical mesoscale circulations (called as vegetation-breeze) are forced by the pressure perturbations with wind speeds about 5 m / s, flowing from the VC to the adjacent     

大气物理过程方案对PM2.5预报的敏感性试验

利用WRF-Chem大气化学模式,选择2015年12月中旬发生在我国的大范围空气污染过程,在采用同样化学方案条件下,针对模式中不同物理过程及其参数化方案开展了地面PM_(2.5)预报的敏感性试验。结果表明:该模式能较好展示此次PM_(2.5)污染的演变过程,与实况也较接近,但对青海经宁夏至内蒙的PM_(2.5)高值区出现了漏报现象,这可能是模式外边界未对污染物做更新所致。对地面PM_(2.5)的预报,各物理过程的敏感度不同,边界层(含近地面层)过程的影响要明显大于积云对流及微物理过程的影响,不同的参数化方案会造成不同的预报误差。边界层过程QNSE和与其配套的近地面层方案的组合是预报的较佳组合;而TEMF和与其配套的组合以及ACM2和Pleim-X的组合则不佳。合理的物理过程参数化方案有助于提高PM_(2.5)预报质量。模式预报对排放源也有适应过程,其Spin-up时间较气象要素长。     

大气运动的动能在湍流边界层中的耗散

在大气能量学中,大气平均运动的动能通过湍流摩擦力转化为脉动运动的能量并最终转化为内能的“变性能量”是大气动能的“汇”,正确估计其大小对研究动能的收支,以至天气系统的生衰,发展均有一定意义。此项由于与湍流通量有关,因而难于从大尺度气象参数计算。以往的研究中,或是将动能平衡方程中其余各项算出,而将此项作为余项估计出,或用地转风由粗浅的边界层模式算出。     

中尺度模式云降水物理方案介绍

介绍大气模式中湿过程的饱和凝结、对流参数化和显式云物理三种方法的适用条件和特点.中尺度大气模式的对流参数化和显式过程一般有多种方案,介绍了多种方案的主要处理方法、基本微物理过程、预报量和主要特点.对流云模式都采用显式过程,但不同模式的湿过程有一定的差别,有各自的微物理特点,对有代表性的云物理方案进行了叙述和分析.     

Representing Model Uncertainty by Multi-Stochastic Physics Approaches in the GRAPES Ensemble

To represent model uncertainties more comprehensively,a stochastically perturbed parameterization(SPP)scheme consisting of temporally and spatially varying perturbations of 18 parameters in the microphysics,convection,boundary layer,and surface layer parameterization schemes,as well as the stochastically perturbed parameterization tendencies(SPPT)scheme,and the stochastic kinetic energy backscatter(SKEB)scheme,is applied in the Global and Regional Assimilation and Prediction Enhanced System-Regional Ensemble Prediction System(GRAPES-REPS)to evaluate and compare the general performance of various combinations of multiple stochastic physics schemes.Six experiments are performed for a summer month(1-30 June 2015)over China and multiple verification metrics are used.The results show that:(1)All stochastic experiments outperform the control(CTL)experiment,and all combinations of stochastic parameterization schemes perform better than the single SPP scheme,indicating that stochastic methods can effectively improve the forecast skill,and combinations of multiple stochastic parameterization schemes can better represent model uncertainties;(2)The combination of all three stochastic physics schemes(SPP,SPPT,and SKEB)outperforms any other combination of two schemes in precipitation forecasting and surface and upper-air verification to better represent the model uncertainties and improve the forecast skill;(3)Combining SKEB with SPP and/or SPPT results in a notable increase in the spread and reduction in outliers for the upper-air wind speed.SKEB directly perturbs the wind field and therefore its addition will greatly impact the upper-air wind-speed fields,and it contributes most to the improvement in spread and outliers for wind;(4)The introduction of SPP has a positive added value,and does not lead to large changes in the evolution of the kinetic energy(KE)spectrum at any wavelength;(5)The introduction of SPPT and SKEB would cause a 5%-10%and 30%-80%change in the KE of mesoscale systems,and all three stochastic schemes(SPP,SPPT,and SKEB)mainly affect the KE of mesoscale systems.This study indicates the potential of combining multiple stochastic physics schemes and lays a foundation for the future development and design of regional and global ensembles.     

Proof of the monotonicity of grid size and its application in grid-size selection for mesoscale models

Terrain characteristics can be accurately represented in spectrum space. Terrain spectra can quantitatively reflect the effect of topographic dynamic forcing on the atmosphere. In wavelength space, topographic spectral energy decreases with decreasing wavelength, in spite of several departures. This relationship is approximated by an exponential function. A power law relationship between the terrain height spectra and wavelength is fitted by the least-squares method, and the fitting slope is associated with grid-size selection for mesoscale models. The monotonicity of grid size is investigated, and it is strictly proved that grid size increases with increasing fitting exponent, indicating that the universal grid size is determined by the minimum fitting exponent. An example of landslide-prone areas in western Sichuan is given, and the universal grid spacing of 4.1 km is shown to be a requirement to resolve 90% of terrain height variance for mesoscale models, without resorting to the parameterization of subgrid-scale terrain variance. Comparison among results of different simulations shows that the simulations estimate the observed precipitation well when using a resolution of 4.1 km or finer. Although the main flow patterns are similar, finer grids produce more complex patterns that show divergence zones, convergence zones and vortices.Horizontal grid size significantly affects the vertical structure of the convective boundary layer. Stronger vertical wind components are simulated for finer grid resolutions. In particular, noticeable sinking airflows over mountains are captured for those model configurations.     

Mesoscale convection from a large-scale perspective

This essay concerns precipitating convective cloud systems and convectively-driven mesoscale circulations (“mesoscale convection”) and their role in the large-scale structure of the atmosphere. Mesoscale convection is an important and ubiquitous process on scales of motion spanning a few kilometers to many hundreds of kilometers. It plays a role in the input of energy to the climate system through the radiative effect of upper-tropospheric cloud and water vapor, and enhanced surface fluxes. This is in addition to its important effect on energy, heat and momentum transport within the atmosphere. However, mesoscale convection is neither parameterized nor adequately resolved in atmospheric general circulation models. Its representation in mean-flow terms raises issues that are quite distinct from classical approaches to sub-grid scale convection parameterization.Cloud-resolving modeling and theoretical concepts pertinent to the transport properties and mean-flow effects of organized convection are summarized, as are the main convective parameterization techniques used in global models. Two principal themes that are relevant to the representation of organized mesoscale systems are discussed. First, mesoscale transports and their sub-grid scale approximation with emphasis on dynamical approaches. Second, long time-scale modeling of mesoscale cloud systems that involves the collective effect of convection, boundary and surface layers, radiation, microphysics acting under the influence of large-scale forcing.Finally, major research programs that address the role of precipitating convection and mesoscale processes in global models are summarized.     

Recent Progress in Numerical Atmospheric Modeling in China

This review summarizes the scientific and technical progress in atmospheric modeling in China since 2011,including the dynamical core,model physics,data assimilation,ensemble forecasting,and model evaluation strategies.In terms of the dynamical core,important efforts have been made in the improvement of the existing model formulations and in exploring new modeling approaches that can better adapt to massively parallel computers and global multiscale modeling.With regard to model physics,various achievements in physical representations have been made,especially a trend toward scale-aware parameterization for accommodating the increase of model resolution.In the field of data assimilation,a 4D-Var system has been developed and is operationally used by the National Meteorological Center of China,and its performance is promising.Furthermore,ensemble forecasting has played a more important role in operational forecast systems and progressed in many fundamental techniques.Model evaluation strategies,including key performance metrics and standardized experimental protocols,have been proposed and widely applied to better understand the strengths and weaknesses of the systems,offering key routes for model improvement.The paper concludes with a concise summary of the status quo and a brief outlook in terms of future development.     

适用于GRAPES模式C-P边界层方案的设计和实现

基于K廓线闭合方案,通过考虑不稳定边界层和稳定边界层中热量交换系数在半层上求取及下边界条件的设置,将温湿倾向在整层上直接计算,设计了Charney-Phillips跳点（简称C-P跳点）的边界层方案,使之与GRAPES全球模式的C-P跳点相协调,解决了Lorenz跳点物理过程与C-P跳点动力框架耦合时插值造成的不协调问题,同时避免了耦合时反复插值造成的误差,提高了边界层物理过程参数化方案及其反馈的准确性和合理性。试验表明:C-P跳点边界层方案因为避免了温度和湿度在垂直方向上的插值,消除了温湿变量在垂直方向上的锯齿状抖动,使温湿廓线分布更合理,减小了模式预报误差,形势场的预报效果也得到一定改善。C-P边界层方案的应用提升了GRAPES全球模式的总体预报性能。     

SENSITIVITY TESTS OF INTERACTION BETWEEN LAND SURFACE PHYSICAL PROCESS AND ATMOSPHERIC BOUNDARY LAYER*

In this paper,an interactive model between land surface physical process and atmosphere boundary layer is established,and is used to simulate the features of soil environmental physics,surface heat fluxes,evaporation from soil and evapotranspiration from vegetation and structures of atmosphere boundary layer over grassland underlying.The sensitivity experiments are engaged in primary physics parameters.The results show that this model can obtain reasonable simulation for diurnal variations of heat balance,soil volumetric water content,resistance of vegetation evaporation,flux of surface moisture,and profiles of turbulent exchange coefficient,turbulent momentum,potential temperature,and specific humidity.The model developed can be used to study the interaction between land surface processes and atmospheric boundary layer in city regions,and can also be used in the simulation of regional climate incorporating a mesoscale model.     

MM5模式中湿方案的不确定性对“03·7”南京暴雨数值模拟不确定性的影响

利用有限区域非静力MM5模式, 分析了显式降水方案对于2003年7月4—5日南京暴雨数值模拟的不确定性影响。采用混合方案模拟此次暴雨时, 这种不确定性决定于显式和隐式方案的相互协调性及敏感性; 隐式方案基本决定了雨带的整体的空间分布, 而显式方案对于降水型及降水量起到一定的调节作用, 调节的程度与选择的参数化方案有关; 采用隐式方案Grell和KF2模拟此次暴雨时, 应考虑不同的显式方案对于降水模拟的不确定性的影响。     

The Major Research Advances of Mesoscale Weather Dynamics in China Since 2003

This paper reviews the main theoretical progress of mesoscale weather dynamics since 2003, including： （1） The dynamic mechanisms of balanced and unbalanced flow are applied to study the genesis and development problems of mesoscale circulation. The symmetric instability and transverse-wave instability are analyzed in line and vortex atmosphere convection, and further research has been done on nonlinear convective symmetric instability. The interaction between forced convection and unstable convection and the wave characteristics of mesoscale motion are also discussed. （2） Intermediate atmosphere dynamic boundary layer models are developed. The complicated nonlinear interaction is analyzed theoretically between the atmospheric boundary layer and the free atmosphere. The structure of the topography boundary layer, atmospheric frontogenesis, the structure and circulation of the low-level front and other boundary layer dynamic problems are discussed. （3） The formation and development of meso-β-scale rainstorms under the background of the East-Asia atmosphere circulation are diagnosed with the variation of MPV （moist potential vorticity） anomalies. And some physical vectors are modified and applied in the moist atmosphere.     

半干旱区植被覆盖度对边界层气候热力影响的数值模拟

在陆－气相互作用的中小尺度系统研究中，水平非均匀下垫面的强迫作用是主要的物理过程。本文利用能量闭合二维陆面过程与大气边界层耦合模式，研究了我国西北半干旱地区（３８°Ｎ，１０５°Ｅ）夏季下垫面物理特征的变化对区域边界层气候的影响。结果表明：土壤湿度、植被覆盖度对局地环流和区域边界层气候的形成起着决定性的作用。模拟结果揭示了在半干旱地区大面积植树造林、提高植被覆盖度，可涵养土壤水分，改善局地生态环境，是人工持续改造干旱、半干旱荒漠地区局地气候的重要途径。     

正压大气中尺度半平衡和准平衡动力学模式

应用描写正压大气运动的基本方程组，分析了中尺度大气运动的物理特征，指出非平衡强迫运动是引起中尺度重要天气演变的根本原因。中尺度动力学方程组是中尺度动力学理论研究的基础，因此，结合中尺度大气运动的基本特征，依据严格的尺度分析理论和摄动理论，简化基于流体力学和热力学的大气动力学方程组使之能够恰当地描述出中尺度运动的基本特征，对于中尺度动力学的发展是极为必要的。基于非线性平衡方程所得到的半平衡和准平衡动力学模式分别与半地转和准地转模式极为相似，它们可以较精确地描述中尺度大气运动的基本特征，因而，可作为中尺度动力学研究的理论基础。将准平衡动力学模式应用于中尺度涡旋系统的研究，结论表明中尺度平衡涡旋系统主要是受梯度风控制，其流场和气压场的发展演变则由一个演化方程来描写，获得了较为理想的结果。     

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

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