VERIFICATION OF TROPICAL CYCLONE RAINFALL PREDICTIONS FROM CMA AND JMA GLOBAL MODELS

The number of tropical cyclone (TC) genesis over the South China Sea and the Northwest Pacific Ocean in 2009 is significantly less than the average (27.4). However, the number of landfall TC over mainland China and its associated rainfall is more than the average. This paper focuses on the performance of numerical weather prediction (NWP) of landfall TC precipitation over China in 2009. The China Meteorological Administration (CMA) and Japan Meteorological Agency (JMA) models are compared. Although the schemes of physical processes, the data assimilation system and the dynamic frame are entirely different for the two models, the results of forecast verification are similar to each other for TC rainfall and track except for TC Goni. In this paper, a day with daily rainfall amount greater than 50 mm was selected as a storm rain day when there was a TC affecting the mainland. There are 32 storm rain days related to the landing of typhoons and tropical depressions. The rainfall forecast verification methods of National Meteorological Centre (NMC) of CMA are selected to verify the models’ rainfall forecast. Observational precipitation analyses related to TCs in 2009 indicate a U-shape spatial distribution in China. It is found that the rain belt forecasted by the two models within 60 hours shows good agreement with observations, both in the location and the maximum rainfall center. Beyond 3 days, the forecasted rainfall belt shifts northward on average, and the rainfall amount of the model forecasts becomes under-predicted. The rainfall intensity of CMA model forecast is more reasonable than that of JMA model. For heavy rain, the JMA model made more missing forecasts. The TC rainfall is verified in Guangdong, Guangxi, Fujian and Hainan where rainfall amount related to TCs is relatively larger than in other regions. The results indicate that the model forecast for Guangdong and Guangxi is more skillful than that for Hainan. The rainfall forecast for Hainan remains difficult for the models because of insufficient observation data and special tropical ocean climate.     

2016年GRAPES_TYM改进及对台风预报影响

为了进一步提高国家气象中心区域模式台风数值预报系统（GRAPES_TYM）的预报能力,2016年对模式参考大气廓线以及涡旋初始化方案进行了改进:由模式初始场水平方向平均的一维参考大气代替原来的等温大气,涡旋初始化方案取消了原涡旋重定位并将涡旋强度调整半径由原来的12°减小到4°。对2014—2016年的生命史超过3 d的所有台风进行了回算,路径及近地面最大风速统计误差分析表明:参考大气的改进可以减小模式对台风预报路径预报的系统北偏和平均路径误差,尤其是140°E以东的转向台风。涡旋初始化方案中强度调整半径的减小会进一步减小模式预报路径的北偏趋势,从而进一步减小平均误差。同业务系统预报结果相比,改进后的GRAPES_TYM（包括参考大气和涡旋初始化）可以使平均路径误差分别减小10%（24 h）,12%（48 h）,16%（72 h）,14%（96 h）以及15%（120 h）。同美国NCEP全球模式路径预报相比,GRAPES_TYM在西行、西北行登陆我国的台风路径预报有一定优势。     

广州“5.7”暴雨预报的模式不确定性研究

2017年5月7日,在弱天气尺度强迫下,广州发生了暖区特大暴雨,局地发展迅速,降水强度极端,多家业务模式出现了漏报情况。为了探究此次降水过程模式预报的不确定性,采用条件非线性最优参数扰动（Conditional Nonlinear Optimal Perturbation related to Parameters,CNOP-P）方法筛选出最能体现中小尺度系统非线性误差增长特征的关键物理参数,以此构造CNOP-P-RP模式扰动方案,并基于CMA-Meso模式进行对流尺度集合预报试验,最后探究了CNOP-P关键参数影响局地对流发生、发展不同阶段的物理机理。结果显示,不同降水阶段的CNOP-P敏感参数主要与垂直扩散、云雨自动转换或其他水成物向雨滴的转换有关。与业务上常用的随机物理倾向扰动（Stochastically Perturbed Parameterization Tendencies,SPPT）方案相比,在本次降水过程中,基于CNOP-P-RP方案构造的集合预报试验具有更高的降水和地面要素的概率预报技巧,集合预报系统可靠性也占优。进一步分析发现,垂直扩散不确定性导致的山前温度梯度和...     

适用于非静力大气模式的完全平衡多矩约束有限体积方法

大气数值模式离散垂直动量方程时,一般而言气压梯度力和重力不易维持严格的静力平衡关系。为精确平衡数值离散的垂直气压梯度力和重力,基于高精度多矩约束有限体积方法引入完全平衡数值公式,即以满足静力平衡关系的热力学参考态对重力源项进行数值离散构造,发展了适用于非静力大气的完全平衡多矩约束有限体积方法。一维标准数值试验表明,完全平衡多矩约束有限体积方法能在较粗糙的计算网格点上保持静力平衡参考态的数值计算误差在计算机的单精度（10-6）和双精度（10-14）水平,在具有小量级扰动的初始条件下,完全平衡多矩约束有限体积方法能较好地模拟扰动的传播,二维非静力热泡试验进一步验证了完全平衡多矩约束有限体积方法对非静力大气运动的模拟能力。数值试验结果验证了所发展方法的完全平衡属性和适用性,这为非静力大气模式发展提供了良好参考价值。     

Implementation of CMA-GFS 3D-Var System on the Yin-Yang Grid

The global three-dimensional variational (3D-Var) data assimilation is implemented on a new quasi-uniform overset (Yin-Yang) grid on sphere. As a quasi-uniform spherical grid, it covers the sphere by overlapping two perpendicularly oriented grid components which is nothing but low latitude region of the usual latitude-longitude grid. Based on this characteristic of the Yin-Yang grid, it enables us to implement the regional 3D-Var system efficiently and accurately on the Yin or Yang component grid, respectively. The global analysis could update directly from the regional analysis since they have the same configurations like the precondition of eigenvalue decomposition for vertical direction, recursive filtering for horizontal direction, minimization method and observation operator and so on. However, the balance equation and vector wind are needed to be paid more attention on the Yin grid analysis due to its coordinate transformation. How to spread the observation information near the boundary of Yin and Yang grid is a key to the 3D- Var analysis. Extending double the horizontal correlation length distance in the overset boundary of Yin and Yang grid has successfully solved the problem. The results show that the analysis on the Yin-Yang grid is reasonable and similar to the result on the latitude-longitude (LAT-LON) grid. This paper provides a promising strategy for the development of a 3D-Var global system for overset grids.     

预估-校正的半隐式半拉格朗日时间积分方案及其在CMA-GFS模式中的应用

CMA-GFS采用的是传统的二时间层半隐式半拉格朗日时间积分方案(SISL).拉格朗日平流速度和非线性项需要采用时间外插进行计算,在急流轴附近等梯度大值区会造成计算不稳定,甚至积分中断现象.文中通过构造预估-校正半隐式半拉格朗日时间积分方案(SISL/P-C),以减少时间外插的影响;半隐式系数由原来的0.72减小到0....     

下一代大气模式中的数值方法综述

数值预报是逐日天气预报、气候预测和气象防灾减灾的核心科技支撑。为进一步提高预报预测的准确度和延长预见期,甚高分辨率、多圈层耦合、多尺度嵌套、多尺度集合、数值地球系统模拟技术等是下一代数值预报的重要发展方向。异构众核高性能计算机和E级计算的高速发展为这一发展提供了契机,但也对现有业务数值预报中采用的数值方法提出了挑战。此文仅对国内外下一代大气模式涉及到的数值方法进行综述,着重于数值算法、准均匀球面网格和时间积分方案等3个方面,期望为相关研究者提供参考。     

陆面过程对2007年淮河流域强降水数值预报的影响分析

基于GRAPES中尺度数值天气预报模式,研究了不同复杂程度陆面过程和不同土壤初始条件对2007年中国夏季淮河流域降水数值预报的影响。重点分析了下垫面非均匀性对造成我国夏季强降水的中小尺度对流系统的启动和发展的影响,探讨了中小尺度对流发生发展中陆气相互作用的可能影响。对2007年7月8日发生在淮河流域的一次对流过程采用两组敏感性试验进行研究,结果表明对流启动对陆面过程有很强的敏感性。Noah陆面模式能合理地模拟出对流的启动,而Slab陆面模式模拟的对流会延迟1～2小时,这与Noah陆面模式能合理描述地表感、潜热通量有关。GLDAS初始土壤资料能更加合理地描述实际土壤温湿的分布状况,有利于更加准确地模拟出对流降水的启动和分布位置。另外,对2007年7月淮河流域的持续性强降水天气过程进行了对比模拟,研究表明:不同复杂程度陆面过程对连续强降水数值预报有明显的影响。使用Noah陆面模式能提高对强降水预报的能力,而且随着模式分辨率的提高,降水的TS评分也在提高,Noah的TS评分总体上高于Slab。     

一种快速预处理算法及其在浅水波方程中的应用

为了提高数值求解大气方程的计算速度,我们研究了稀疏逼近逆预处理方法及其在数值求解浅水波方程中的应用,这是一种求解大型线性方程组的快速算法,该算法的核心内容是稀疏逼近逆非零元模式的选取。首先导出了一种稀疏逼近逆的非零元模式及其确定方法,然后以浅水方程的差分格式为例,借助于GMRES迭代算法,对这种预处理快速算法应用前后的计算速度进行了比较,发现该快速算法能大幅度提高运算速度。另外,该预处理快速算法简单、易于并行,是一种值得在大气方程中推广应用的方法。