On temperature initialization in primitive equation forecast models

Four procedures of specifying model initial temperature were described and tested in the present study. It was found that the use of observed temperatures along with a proper vertical interpolation scheme was not only acceptable, but produced less error than the use of temperatures derived from geopotential height through the hydrostatic equation did. Use of the difference form of the hydrostatic equation would produce unacceptable errors in the initial temperatures, unrealistic horizontal and vertical distribution of temperature, and these errors would influence the calculation of the pressure gradient force, resulting in substantial, artificial disturbances within the model domain.In addition, an approach to check the initial data was described. Taking advantage of the fact that the geostrophic wind in sigma coordinates should be nondivergent, geopotential height and temperature were used to calculate the pressure gradient force terms and an initial divergence of the geostrophic wind. This approach can he used for comparing different initialization schemes for identical input data.     

GRAPES模式标准初始化方案设计与实现

文章介绍了GRAPES模式标准初始化系统的方案设计及程序结构,对3种不同方案产生的初始场进行了分析。结果表明：水平风场、温度场、位温场对插值方法的精度并不是很敏感,而高度场以及气压Exner函数则对垂直插值方法的精度要求很高;在3种模式变量产生方案中,方案一的效果较差,但它计算速度快,在地形梯度不大的地方可以选择使用;方案二与方案三的效果则比较好,用户可根据自己的需求任意选择。GRAPES模式标准初始化系统在完成其主要功能的同时,通过灵活的程序设计,对标准初始化系统的区域、分辨率、垂直层次、水平插值方法、垂直插值方法等的任意选取,使得标准初始化系统可以很容易满足模式初值的各种需求,为模式开发研究实验提供了各种便利。     

DRP-4DVar方法同化AIRS反演资料在一次江淮流域暴雨中的应用

利用经济省时的降维投影四维变分同化方法(DRP-4DVar),在2009年7月22～23日江淮流域的一次大暴雨过程中同化晴空条件下高光谱大气红外探测仪(AIRS)反演温度、湿度廓线,改进此次强降水过程的模拟。试验结果分析显示,同化AIRS反演的温度及湿度场后,基于四维变分同化系统的模式约束,能够改进湿度场、高度场、高低层散度场。从累积降水量偏差图及同化试验增量图可以看到,正降水量偏差对应于正湿度增量、负位势高度增量及低层负散度高层正散度增量,负降水量偏差则与之相反。同化试验较参照试验可更好地模拟出暴雨的天气形势、对暴雨的落区及强度有更好的反映。此外,从单次同化与连续同化的试验对比结果看出,连续同化试验结果较单次同化结果有进一步的改进,说明不断加入新的观测资料可以更好地模拟强降水过程。     

西北太平洋副热带高压强度的统计预报

本文对1982年7—9月中央气象台发布的东亚范围500百帕的48小时预报图进行了误差分析，发现40°N以北的中高纬度预报误差较小，而在35°N以南的较低纬度，尤其是副热带地区，预报误差较大。因此我们对西风带高度场和副热带高度场分别进行了车贝雪夫多项式和自然正交函数展开。分析表明，他们之间具有一定的关联。在此基础上，提出了一个用车贝雪夫系数去预报同一时刻副热带高压系数（时间系数），从而获得48小时副热带高度场的统计预报方法。通过独立样本检验并与B模式的预报结果进行对照，证明本文的预报方案具有参考和应用价值。     

Estimating the error in vertically interpolated forecast‐observation differences: U‐ and v‐wind components

Abstract

Data assimilation in numerical weather forecasting corrects current forecast values by subtracting a portion of interpolated forecast‐minus‐observation differences at the points of a three‐dimensional grid. Deviations used in updating a forecast data field are forecast errors obtained or derived from observations available at update time. When observations are missing at mandatory levels, construction of full vertical soundings by interpolation introduces extraneous errors. The present paper is concerned with determination of the error in vertical extrapolations of surface winds, and of aircraft and satellite cloud‐tracked winds. In addition it examines the effect on accuracy of using location‐specific statistics compared to averaged statistics as the basis for the interpolation weighting scheme and compares errors of one‐ and two‐variable interpolations.

Interpolation accuracy tests demonstrate the influence of the interpolation scheme on the quality of interpolated information used in forecast updating. The results show that the level of accuracy exceeds the benchmark provided by monthly mean forecast error values only with bivariate interpolation of wind components from off‐level data sources.     

Blacklist Design of AMDAR Temperature Data and Their Application in the CMA-GFS

Blacklist methods are used in the CMA Global Forecasting System (CMA-GFS) to improve the application of aircraft temperature data to numerical weather prediction in the Northern Hemisphere and the tropics. In this paper, the ERA5 re-analysis data are used to analyze aircraft temperature observation errors of each aircraft and a blacklist is established using pre-quality controls and threshold methods. The blacklist-filtered and blacklisted aircraft temperature data are then applied to the four-dimensional variational assimilation system, respectively, and an assimilation cycle forecast for the period from September 1 to 30, 2019 is carried out. The results show an uneven distribution in the global aircraft blacklist data. After the application of the blacklist methods, the RMSE of geopotential height and temperature analysis field decrease in the vertical direction by a maximum of ~ 1.5 gpm at 200 hPa and ~ 0.15 K at 250 hPa, respectively. Overall, the blacklist methods of aircraft temperature data improve the analysis and forecast in the CMA-GFS.     

T63模式月动力延伸预报高度场的改进实验

为克服T63模式月动力延伸预报中纬向平均环流的系统性误差较大的情形,文章利用NCEP/NCAR逐候再分析500 hPa高度场资料和非线性时空序列预测理论的局域近似法进行逐候纬向平均高度距平场预报.近30组个例的预报效果分析表明,就1～3旬总体而言,非线性时空序列预测方法对纬向平均高度距平场的预报优于持续性预报和模式动力延伸预报,体现了改善纬向平均高度场的能力.尤其是第3旬的预报,当持续性预报偏差与实况偏差明显增大、动力预报技巧相对于第1旬和第2旬降低时,相空间重构结果仍然保持一定的优势.     

Testing the Ability of Numerical Model to Simulate Climate and Its Change With 4D-EOF Analysis

The four-dimensional empirical orthogonal function (4D-EOF), which in reality is a simple combination of three-dimensional EOF (3D-EOF) and extended EOF (EEOF), is put forward in this paper to test the ability of numerical model to simulate climate and its change. The 4D-E0F analysis is able to reveal not only the horizontal characteristic pattern of analyzed variable, and its corresponding annual and inter-annual variations, but also the vertical structural characteristics. The method suggested is then used to analyze the monthly mean 100-, 500-, 70G-, and 1000-hPa geopotential height fields (4941 grids and grid spacing 60 km) and their anomaly fields in 1989-1998 simulated by the MM5V3 from the RMIP (Regional Climate Model Inter-comparison Project for East Asia)-II, as well as their counterparts (used as the observed fields) from the NCEP/NCAR re-analysis dataset in the same period. The ability of MM5V3 in simulating East Asian climate and its change is tested by comparing the 4D-EOF analysis results of the simulated and observed datasets. The comparative analyzed results show that the horizontal pattern of the first eigenvector of the observed monthly mean geopotential height fields and its vertical equivalent barotropic feature were well simulated; the simulations of the first two eigenvectors of the observed monthly mean geopotential height anomaly fields were also successful for their horizontal abnormal distributions and significant equivalent barotropic features in the vertical were well reproduced; and furthermore, the observed characteristics, such as the variation with height, the annual and inter-annual variations of the monthly mean geopotential height/anomaly fields were also well reflected in the simulation. Therefore, the 4D-EOF is able to comprehensively test numerical model's ability of simulating the climate and its change, and the simulation ability of MM5V3 for the climate and its change in East Asia in the 1990s was satisfactory.     

探空观测黑名单检查在变分同化系统中的应用

针对探空观测资料使用造成的某些区域GRAPES分析场存在虚假的高、低压系统问题,该文通过对比全球探空资料的位势高度观测与NCEP分析场,统计站点中观测质量较差的时次出现频数,确定探空位势高度观测黑名单。研究表明:500 hPa在印度地区、北大西洋和南极洲附近的探空位势高度观测与NCEP分析场的均方根误差在30 gpm以上的站点较多,且位势高度观测不可靠观测比率为20%以上的站点主要集中这些区域,以上观测站均列入黑名单。文中在GRAPES全球三维变分分析场的质量控制中加入探空位势高度观测黑名单检查,通过6 h分析预报循环试验表明:探空位势高度观测黑名单检查能有效提高分析场质量,GRAPES位势高度分析场在南极洲附近和印度地区有所改善。     

数值天气预报业务模式现状与展望

采用中国气象局2014年6月1日—30日14时加密探空资料,利用华东区域中尺度数值预报业务系统比较同化加密探空观测资料前后模式预报结果的差异。研究表明, 同化加密探空资料后,对模式初始时刻不同高度的位势高度、比湿、温度、风速等变量均有一定的影响;对位势高度、温度和风场的影响在高层100—150 hPa比较显著,而对比湿的影响主要体现在低层700—750 hPa。同化加密探空资料后模式初始场更接近实况。批量数值试验的统计检验表明,同化加密探空观测资料后对强降水及形势场预报均有不同程度改进,24 h暴雨和大暴雨量级降水的预报技巧分别提高了2.5%和8.1%。

     

GRAPES全球三维变分同化中卫星微波温度计亮温的背景误差及在质量控制中的应用

模式变量背景误差在观测空间的投影,也即观测变量的背景误差包含了变分同化系统的重要信息,其在诊断和分析变分同化系统中资料的影响等方面具有重要作用,特别是在背景场检查质量控制中。在GRAPES全球三维变分同化(3DVar)系统中仅给定了控制变量的背景误差,并未直接给定观测变量的背景误差。为了能够对GRAPES全球3DVar进行全面的诊断和分析,改进卫星微波温度计资料的质量控制,推导出GRAPES全球3DVar同化系统控制变量随机扰动方法估计观测变量的背景误差的公式,为分析和改进GRAPES全球3DVar提供了一个有力工具,并进而估计了AMSU-A亮温的背景误差,分析了AMSU-A不同通道亮温的背景误差特征,将其应用于GRAPES全球3DVar的AMSU-A亮温的背景场检查质量控制中。结果表明,控制变量随机扰动方法估计的GRAPES全球3DVar同化系统AMSU-A亮温的背景误差正确合理。同化循环预报试验结果表明,亮温的背景误差在背景场检查中的应用显著提高了GRAPES全球3DVar同化的亮温资料的数量,显著提高了GRAPES南半球对流层中高层位势高度场的预报技巧。在GRAPES全球3DVar同化系统中推导和实现的控制变量扰动方法为诊断和分析GRAPES全球3DVar观测资料同化效果提供了有力工具。      

基于T213集合预报的延伸期产品释用方法及初步试验

基于T213集合预报系统2008年1月中国及附近区域500 hPa高度和850 hPa温度的1~15 d预报资料,构建延伸期产品释用方法,通过对逐日11~15 d预报资料做集合平均和后向衰减权重系数滑动平均,进而得到延伸期(11~15 d)候平均和候距平预报,并对预报效果进行检验,结果表明:对11~15 d预报场做集合平均和后向衰减权重系数滑动平均均能降低预报误差,改善整体预报效果。由此得到的500 hPa高度场和850 hPa温度场11~1 5 d候平均预报误差与逐日控制预报第5d的水平相当,候距平相关系数均接近0.6,整体而言效果较好,具备一定的应用价值。500 hPa高度场和850 hPa温度场11~15 d候距平预报在中国大陆地区位相准确率均较高,东南沿海和东北部分地区稍差,且850 hPa温度场的位相准确率整体高于500 hPa高度场。候距平预报对延伸期(11~15 d)的大范围持续性异常距平具有较强的捕捉能力,对异常距平出现的范围和分布、强距平中心的位置的预报均较好,但强度整体偏弱。     

GRAPES全球集合预报初始条件及模式物理过程不确定性方法研究

为描述GRAPES全球模式初始条件的不确定性,基于适合集合预报应用的GRAPES全球奇异向量技术,依据大气初始误差符合正态分布的特征,采用高斯取样奇异向量来构造全球集合预报初始扰动,在此基础上建立了GRAPES全球集合预报系统(GRAPES-GEPS)。利用GRAPES全球同化分析场,对采用初始扰动的GRAPES-GEPS连续试验预报结果进行检验和分析。结果表明:GRAPES-GEPS中高度场、风场及温度场预报的集合离散度能有效快速增加,集合平均均方根误差与集合离散度的关系合理;相对控制预报的均方根误差,集合平均的预报优势在预报中期非常显著。为进一步体现GRAPES-GEPS中模式物理过程的不确定性,发展了模式物理过程倾向随机扰动技术(SPPT)。试验结果表明:SPPT方案的应用有效提高了GRAPES-GEPS在南、北半球和热带地区等压面要素预报的集合离散度,同时一定程度减小了集合平均误差,进而改进了集合平均误差与集合离散度的关系,其中SPPT方案在热带地区的改进最为显著。本文发展的基于奇异向量的初始扰动方法和模式扰动SPPT方案在中国气象局2018年12月业务化运行的GRAPES-GEPS中得到了应用。     

利用MODIS红外资料反演大气温湿度廓线的研究

概述了利用特征向量统计回归反演算法,从EOS/MODIS的红外通道资料反演大气温湿度垂直分布的过程,并与美国国家环境预报中心NCEP(National Centers for Environmental Prediction)的等压面再分析场资料按照纬度和气压高度进行了真实性检验。结果表明：由MODIS资料反演得到的大气温湿度参数能够揭示大气温湿度的垂直分布。在各个等压面上均方根误差平均值在中纬度地区为3.39K,低纬度地区为1.40K,近地面层、对流层顶附近及下垫面地形复杂的区域误差较大,总体上低纬度地区要好于中纬度地区。反演的水汽误差也为低纬度地区小于中纬度地区,且随高度升高,中、高纬度误差都逐渐减小并逐渐接近。     

Comparison and Analysis of Two Different Calculation Schemes of SLEVE- Hybrid Coordinate in GRAPES_MESO Model

The calculation scheme of the smoothed-level and hybrid (SLEVE-hybrid for short) coordinates in numerical forecasting model is not limited to one. It is divided into the semi-analytical scheme and the finite differential scheme in terms of the various differential methods of the coordinate deformation variables. Comparing the dynamic equation and the long-time batch simulation results of the two schemes, the present study draws the following conclusions. The first- order finite difference accuracy of the coordinate deformation variables in the finite differential scheme is theoretically lower than that in the semi-analytical scheme. The larger the vertical gradient of the layer thickness is, the larger the relative errors of the finite differential scheme are. The long-time batch simulation test in the GRAPES model dynamic core demonstrates that the bias of the temperature and the geopotential height in the semi-analytical scheme is smaller under the default layering, while the simulation difference of the two schemes is greatly reduced when the layering is more uniform.     

两种地形追随坐标对一次强降水过程模拟影响对比分析

利用GRAPES_Meso 区域中尺度模式,选择传统追随坐标(简称Gal.C.S 坐标)和新的平缓追随坐标(简称SLEVE1坐标)两种高度地形追随坐标,对一次强降水天气过程进行数值模拟,对比分析两种不同地形追随坐标的影响。个例模拟结果表明,SLEVE1 对降水场、风场、温度场、高度场的模拟预报均有改进;与Gal.C.S 坐标相比,SLEVE1 坐标通过平滑坐标面上的地形作用,减小虚假垂直速度进而减少虚假降水;SLEVE1 坐标模拟的各要素场剖面图等值线上的小尺度噪音有不同程度衰减,预报要素均方根误差减小,相关系数增大;SLEVE1 坐标临界衰减高度越小,模拟预报结果越好,该高度受到最大地形高度的约束。     

2014年夏季长江流域持续性低温的大尺度环流异常

采用中国地面站气温逐日观测资料、NOAA全球逐日海表温度资料,及NCEP/NCAR的全球日均再分析资料,研究了2014年持续性低温的三维结构及大尺度环流异常。结果表明,2014年的低温异常,除了在陆地上区域性地出现在长江流域,还以大尺度带状的形式、从陆地延伸到海洋上。这种带状异常不只出现在近地面,在大气各层（925～500 hPa）都能看到。分析指出,大气中的这个低温带主要由高纬大气环流异常造成。在位势高度场上,最重要的异常出现在高纬60°N,有呈带状的位势高度正距平,它引导（距平意义上的）偏北气流从正北和东北偏东方向侵入,在其南侧形成一带状的偏低温区。大气各层均呈现出这种在高纬有位势高度正距平、相应地在稍南的低纬（40°N）有位势高度负距平、两者之间为低温区的分布特征。从低层往上,这种配置型式整体表现出由南往北的倾斜,其垂直剖面表现为距平意义上的、大尺度、类似于锋面的倾斜结构。文中用简单的概念模型对此进行理解,认为这种结构是由大气动力异常和热力异常相互影响、共同作用而形成的。     

华东冬季异常冷暖与大气环流和海温的关系

利用1951-2007年华东地区14个代表站冬季(12-2月)温度资料和北半球500 hPa高度及北太平洋海温资料,通过合成分析、相关分析等方法,研究了华东地区冬季气温的气候变化及其与北半球500 hPa高度场、北太平洋海温场的关系.结果表明:华东地区冬季气温具有明显的年代际气候变化特征;前期夏季北半球500 hPa高度距平场和前期春季北太平洋海温距平场分布可作为华东冬季异常冷暖年的前兆信号;夏季北太平洋中部地区500 hPa高度场变化及前期10月西太平洋副高强弱变化,对华东地区冬季气温变化具有很好的指示性;春季南赤道海流区和西风漂流区海温异常变化,对华东地区冬季气温变化也具有很好的指示意义.     

Assimilation of Ocean Surface Wind Data by the HY-2B Satellite in GRAPES: Impacts on Analyses and Forecasts

The ocean surface wind(OSW) data retrieved from microwave scatterometers have high spatial accuracy and represent the only wind data assimilated by global numerical models on the ocean surface, thus playing an important role in improving the forecast skills of global medium-range weather prediction models. To improve the forecast skills of the Global/Regional Assimilation and Prediction System Global Forecast System(GRAPES＿GFS), the HY-2B OSW data is assimilated into the GRAPES＿GFS four-dimensio...     

Preliminary Comparison of the CMA, ECMWF, NCEP, and JMA Ensemble Prediction Systems

Based on The Observing System Research and Predictability Experiment (THORPEX) Interactive Grand Global Ensemble (TIGGE) dataset,using various verification methods,the performances of four typical ense...