高分辨有限区数值预报模式的省时显式积分

本文提出利用省时显式方案积分两个高分辨有限区数值预报模式。其结果表明，这一积分方案可以使模式的计算效率大大提高，而由时间积分方案的变化所造成的预报差异比模式空间离散化方案的不同所造成的预报差异要小。同时还讨论了将省时显式方案用于有限区模式的边界处理问题，提出了一个与蛙跃格式等价的边界处理方案。     

NUMERICAL ANALYSIS OF THE DIFFERENCES IN DIFFUSION BETWEEN p－SURFACE AND σ－SURFACE

Dimensional analysis and reduction are done to two existing schemes of 4th-order linear horizontal diffusion, and detailed control experiments between them are made using a topographyincluded mesoscale model. Horizontal diffusion is calculated on the or surface in one (known as Scheme A afterwards ), and on the p-surface in another (Scheme B). Experiments show that differences are small in smooth-terrain areas and very large in steep mountain areas, with the 24h rainfall prediction deviating by 50 mm between forecasts of the two sChemes. The reason may be that temperature and humidity are falsely diffused in Scheme A, which causes abnormal temperature and humidity, and results in the anomalies of the unstable layer and convective processes. In addition, Scheme A could also bring about circulation anomalies which assumingly have direct link to the convective anomalies in the scheme. Furthermore, perturbation may also affect surrounding areas by wave-like propagation such that precipitation anomalies may occur in the area. The analysis indicate that Scheme B is necessary and feasible for it minimizes diffusion-involved forecast abnormality in steep mountains and areas around.     

新世纪初我国数值天气预报的科技创新研究

概要介绍最近5年在国家科技攻关项目“中国气象数值预报系统科技创新研究”框架内所取得的主要成果，重点是卫星等遥感资料在变分同化中的应用，高分辨非静力数值预报模式的发展，全球资料同化与中期数值天气预报系统的发展，数值天气预报系统的模块化与并行计算，数值天气预报新技术的研究等，并扼要介绍我国新一代数值天气预报系统的业务应用试验。最后，对我国数值天气预报的进一步发展做了讨论。     

非线性平衡方案在三维变分同化系统中的应用

文中考虑气流出现明显的弯曲和辐合时,资料变分同化方案中的风压场间动力学约束问题。在强涡旋、锋面等剧烈变化的系统中,地转平衡和线性平衡方程均不能很好地描述风压场关系,非线性平衡方程考虑了流场曲率的作用,能更恰当地反映风压场的关系。把非线性平衡方案应用于三维变分同化分析系统中,可以得到和天气形势相关联的模式变量背景误差协方差,对热带气旋进行同化试验后得到的分析场更协调。结果表明,非线性平衡方程的应用,可以提高强涡旋系统中三维变分同化分析场的质量,对热带气旋强度预报有一定的改善。     

雷达反射率因子在中尺度云分辨模式初始化中的应用Ⅰ：云微物理量和垂直速度的反演

对于水平网格距小于10 km的高分辨率非静力平衡的显式云分辨模式,云微物理量的初始化以及物理量之间的相互协调是十分重要的,并且一直是云分析领域的一个难题。考虑在云雨处于定常状态的前提下,根据暖云过程,可得到云中水成物之间相互转化以及与垂直速度的约束关系,而雨水含量跟雷达反射率因子(Z)有关。因此,采用合肥新一代天气雷达2003年7月5日02时(北京时)的观测资料,针对雷达探测的特点对雷达原始数据进行了坐标转换,并进行了插值处理和简单的质量控制,然后依据Z-qr关系和定常暖云方案,反演了雨水混合比(qr)、云水混合比(qc)、水汽混合比(qv)和垂直速度(w)。结果表明,由此得到的雷达回波主要特征及量值分布与雷达站的CAPPI图像基本一致;水汽、云水、雨水和垂直速度量值大小的分布与雷达回波强度的分布也是吻合的,体现出了云中水成物和垂直速度的三维分布结构;各物理量的量值也符合梅雨锋暴雨的特点,梅雨锋积层混合云系中层状云和对流云的差别十分明显。在较强的回波区,雨水在6 km以下含量较大,最大值位于4 km左右,可超过3.0 g/kg;上升速度在5 km左右最大,最大值超过5 m/s;云水含量大值区位于5 km以上,在7 km高度上达到最大,超过3.0 g/kg;雨滴末速度虽然上下比较一致,一般几米每秒,但在5 km左右为大值区。     

GRAPES的新初始化方案

四维变分同化由于引入预报模式作为一项约束,理论上它的分析场已经具有较好的平衡性,但实施时还会有诸多因重力波导致的高频振荡过程,因此,四维变分同化(4DVar)分析仍需要初始化。文中描述了GRAPES全球四维变分同化系统(GRAPES-4DVar)的新初始化方案的科学设计、公式演绎以及试验结果。GRAPES-4DVar的新初始化方案采用数字滤波方案作为代价函数的一项约束控制重力波引发的不平衡结构,约束强加在分析增量上与极小化迭代过程同步进行。新的初始化方案是变分同化系统的一部分,数字滤波的积分时间与4DVar的同化时间窗一致,不会对4DVar产生额外的计算资源消耗;并能适应长时间窗的同化,不会因为时间窗的延长而削弱慢波过程。新初始化方案中,模式轨迹的光滑程度可在变分同化中通过重力波控制项的权重系数方便控制。GRAPES全球四维变分同化的理想和循环同化批量试验都表明,在四维变分同化中,重力波的控制依然非常重要,具有初始化的GRAPES试验,无论分析还是预报技巧都较无初始化的有明显优势。与以前分析和滤波独立实施的旧初始化方案相比,新方案的分析和预报效果略优,同时有效地节省循环同化系统的运行时间,这对四维变分同化来说非常重要。     

资料同化对2017年登陆广东沿海台风的短期降水与路径预报影响

利用华南精细数值天气预报模式,设计了无同化资料(CTRL)、同化雷达反演水汽(EXP1)以及同化雷达反演水汽、地面和探空资料(EXP2)三个试验,对2017年登陆广东沿海的四个台风降水预报与路径预报进行模拟,以评估资料同化对登陆台风短期降水预报、路径预报的影响。分析结果如下:雷达反演水汽同化后对未来24小时降水预报技巧均有正的改善,对台风路径预报影响不大;在此基础上同化地面、探空资料后对台风路径预报有改进,对降水预报改进不明显(与EXP1比)。通过诊断分析台风“玛娃”,发现模式初值场水汽的增量配合对流上升区有利于短时间内成云致雨,从而提高短时降水预报;地面及探空资料同化有利于登陆台风的短时路径预报。      

A STUDY ON THE APPLICATION OF FY-2E CLOUD DRIFT WIND HEIGHT REASSIGNMENT IN NUMERICAL FORECAST OF TYPHOON CHANTHU (1003) TRACK

In this paper, we first analyzed cloud drift wind(CDW) data distribution in the vertical direction, and then reassigned the height of every CDW in the research domain in terms of background information, and finally, conducted contrast numerical experiments of assimilating the CDW data before and after reassignment to examine the impacts on the forecast of the track of Typhoon Chanthu(1003) from 00:00(Coordinated Universal Time) 21 July to 00:00 UTC23 July, 2010. The analysis results of the CDW data indicate that the number of CDWs is mainly distributed in the midand upper-troposphere above 500 h Pa, with the maximum number at about 300 h Pa. The height reassigning method mentioned in this work may update the height effectively, and the CDW data are distributed reasonably and no obvious contradiction occurs in the horizontal direction after height reassignment. After assimilating the height-reassigned CDW data, especially the water vapor CDW data, the initial wind field around Typhoon Chanthu(1003) became more reasonable, and then the steering current leading the typhoon to move to the correct location became stronger. As a result, the numerical track predictions are improved.