GRAPES全球三维变分同化业务系统性能

近年来,GRAPES全球三维变分同化系统分析性能和稳定性有了长足进步。该文简要介绍了近两年GRAPES全球:三维变分同化技术的发展与改进情况,包括同化框架技术、资料同化应用技术与系统稳定性等方面。分析诊断了两年的同化循环试验结果,以探空资料作为参考,对ERA-Interim再分析场、NCEP FNL分析场和GRAPES全球三维变分分析场的统计特征进行了比较;以ERA-Interim再分析场作为参考,对NCEP FNL分析场、T639分析场和GRAPES全球三维变分分析场进行比较。结果表明:GRAPES分析场的质量明显优于T639分析场,性能上达到了业务化的要求,但相比NCEP FNL分析场还有一定差距,特别是对流层内湿度分析场的误差还比较大。     

COSMIC资料在GRAPES全球三维变分同化系统的初步研究

利用COSMIC获得的GPS无线电掩星数据所反演得到的大气温湿度廓线资料,在GRAPES全球三维变分同化系统中做连续循环同化试验,研究了加入COSMIC资料后得到的GRAPES分析场及其5天中期预报是否有所改善。研究结果表明：通过加入COSMIC资料连续循环同化得到的分析场,相对原来没有加COSMIC资料连续同化出来的GRAPES分析场有明显改善,由此GRAPES全球模式全球5天中期预报水平有明显提高,其中加入没有稀疏化后COSMIC资料对南半球GRAPES分析场的改善以及GRAPES全球模式5天中期预报水平的提高都比较显著;加入稀疏化后的COSMIC资料对北半球GRAPES分析场的改善以及GRAPES全球模式5天中期预报水平的提高都有比较显著的效果。     

GRAPES全球三维变分同化系统的检验与诊断

中国气象局数值预报中心新近升级的GRAPES全球三维变分同化系统的大气基本状态变量在物理属性与定义的网格和坐标上与预报模式保持一致,是一个完全针对GRAPES预报模式的同化系统。该系统不仅有利于减小分析误差,也是构建GRAPES四维变分同化系统的基本环节之一。该文通过与观测资料的对比、与国际其他业务中心分析场的对比,以及中期数值预报的检验,对新的GRAPES全球三维变分同化系统性能进行较全面讨论,并通过对这一系统的检验,探索资料同化系统性能的检验方法,尤其是观测资料同化效果的定量评价方法。诊断结果表明:在宏观特征上,GRAPES变分同化系统的分析场与欧洲中期数值预报中心和美国国家环境预测中心的分析场十分相似, 但细节上仍有差别。这些差别主要源自GRAPES同化系统中探空、地面报、掩星以及飞机报观测的贡献偏大,而卫星垂直探测仪观测资料的作用尚未充分发挥。从探测单要素来讲,风及湿度观测的作用发挥不够。此外,青藏高原周围地区、模式高层及赤道地区分析场偏差较大,它们与模式地形及高层的处理等有关系,这些问题有待进一步改进。     

集合变分混合同化方案在快速循环同化系统中的应用研究

基于北京市气象局快速循环同化系统RMAPS-ST以及对流尺度集合预报系统RMAPS-EN,构建了En-3DVAR集合变分混合同化系统,将该系统应用到业务快速循环同化系统中并进行试验,分别在冷启动与循环启动环境下对比了混合同化系统（Hybrid）与三维变分（3DVAR）的同化预报效果。获得的结论如下:单点试验结果表明,混合同化系统分析增量的分布与集合预报离散度分布具有较好的对应关系;在冷启动和循环启动中,三维变分的分析增量都表现出各向同性的特点,混合同化分析增量均表现出一定的流依赖特征;降水个例分析表明,在冷启动环境中,Hybrid与3DVAR效果相当,而在循环启动中,Hybrid的降水预报相对于3DVAR有较明显的改进效果;批量试验检验结果表明,冷启动中,Hybrid与3DVAR的评分大致相当,而在循环启动中,Hybrid相对于3DVAR的评分有明显改进;集合离散度和背景场误差的相关性分析表明二者在循环启动环境下具有更好的相关性。     

极轨卫星高层通道在GRAPES同化系统中的应用

针对GRAPES(Global/Regional Assimilation and Prediction System)模式三维变分系统高层背景场温湿廓线外推方案的局限性,提出以气候垂直廓线重新构造高层温湿垂直结构,以减小外推方案的偏差。首先采用一维变分同化系统,展开模拟实验:分析目前模式中使用的外推方案误差及其对反演结果的影响,利用高层大气气候廓线构造垂直结构并分析同化偏差。最后,运用GRAPES全球分析预报系统进行同化实验并分析改进程度。结果显示:模拟研究表明采用高层背景场温湿廓线外推方案与实际观测相比最大偏差在1 h Pa附近可达数十度以上,不仅影响平流层,而且对对流层也有影响;用气候温度数据修正GRAPES高层温度数据,可以减少50%以上的偏差,证明了用气候值高层数据优化现行GRAPES模式中同化系统高层插值方案的可行性。全球GRAPES三维变分同化试验结果显示,改进方案不仅显著的改善平流层分析质量,对对流层中高层也有改进。     

GRAPeS三维变分同化系统的理想试验

文中介绍了一种新的、适合格点模式的三维变分同化方案GRAPeS 3DVAR。该方案采用相互独立的流函数、非平衡速度势函数、非平衡位势和水汽作为分析变量。通过变量变换对目标函数进行预调节 ,不仅避免了直接计算背景误差协方差逆矩阵的困难 ,而且改善了Hessian矩阵的性状 ,提高了收敛速度。采用EOF分解方法 ,将三维分析变量投影到垂直摸态上 ,分解成为二维场 ;水平方向采用数字 (递归 )滤波器代替矩阵运算 ,实现和简化了方案的求解。此外 ,还考虑了质量场和风场之间的平衡约束关系。理想试验结果表明 ,GRAPeS 3DVAR能够正确地反映多变量之间相互作用关系 ,收敛迅速 ,分析结果合理     

雷达TREC风的三维变分同化应用与试验

介绍了雷达TREC风场在中国气象科学研究院开发的三维变分同化系统GRAPeS-3DVar中的应用,研究了GRAPES-3DVar系统同化雷达TREC风场资料的一些基本特性,发现雷达TREC风场在GRAPES-3DVar中有较好的应用价值,能有效地提高分析产品的质量,并通过一个热带气旋登陆变性过程的数值试验,说明了在GRAPES-3DVar中加入雷达TREC风场也能改善预报效果。     

华北地区一次大雾过程的三维变分同化试验

文中采用WRF非静力数值预报模式及其三维变分同化系统(WRF3D-Var),对2006年1月13—14日发生在华北地区及山东半岛的一次大雾过程进行了包括GTS(Global Telecommunication System)资料、AMDAR(Aircraft Meteorological Data Relay)资料和9210资料的不同资料组合的三维变分同化试验,以及时间间隔分别为6、3和1h不同时间频率的循环同化试验,并以同化分析场为初始场进行了36h的模拟试验。对同化分析场和模拟结果进行了分析,分析结果表明,采用三维变分方法同化AMDAR等多种非常规观测资料后,分析场均有明显的改变,对雾区的模拟结果也有局部不同程度的修正。进一步分析起修正作用的原因得知同化资料后对低层的湿度和层结趋稳性有所改善。同化GTS资料对低层的增湿贡献明显,但对层结趋稳性贡献不大;而同化AMDAR资料主要使层结趋稳性明显,对增湿无贡献;9210资料对低层湿度和层结趋稳性均有贡献。不同时间间隔的循环同化试验表明,多时次的循环同化比单时次的同化分析增量要大,逐时循环同化与6和3h循环同化相比,可明显改善模拟效果。     

全球GRAPES等压面三维变分分析预报循环系统及试验

利用全球通讯系统GTS传输的常规观测资料和卫星资料,用全球与区域同化预报系统GRAPES(global/regional assimilation and prediction system)全球等压面三维变分分析和中期数值预报模式进行为期两个月6 h分析预报循环试验,并对分析结果进行了诊断分析。结果表明,位势高度在中低纬对流层顶及以上区域存在5~20 gpm负偏差,在100 hPa以下的南北极区域存在5~20 gpm正偏差;相对湿度在南极存在较大负偏差,在高层有较大正偏差。风场在赤道中高层有较明显的偏差。在500 hPa上北半球GRAPES位势高度与NCEP资料的均方根误差在10 gpm左右,南半球在15~24 gpm之间。与探空观测相比,经变分分析后,分析的偏差和均方根误差都有所减小,但位势高度背景场存在明显负偏差,而且在100 hPa以上位势高度改进不如低层明显。全球等压面三维变分分析预报循环试验表明,分析预报循环系统能稳定运行,但是分析和6 h预报场还存在系统性偏差。     

雷达资料在云南一次强降水过程中的三维变分同化试验

鉴于云南观测信息相对不足、局地强降水突出的现状,利用WRF(Weather Research and Forecasting)模式及其变分同化系统进行雷达反射率因子和反演风场的三维变分同化试验。通过对2012年9月12日00:00—13日00:00发生在云南的一次强降水过程进行数值模拟和对比分析,结果表明,同时同化雷达反演风场和基本反射率因子,对区域模式同化系统中风矢量、相对湿度、位势高度几个基本分析量都有明显影响。雷达资料的同化,有利于区域模式初始场中强降水区域的上游中低层空气湿度增加、水汽输送增强和强降水发生区域的风场辐合加强,从而改善区域模式对强降水落区、强度的预报质量。对于切变线等天气尺度系统影响下的强降水过程,雷达资料的同化持续时间选取3 h、同化间隔为1 h较适宜。另外,雷达反演风场和基本反射率因子的同化均对降水预报改善有明显贡献,且多种资料的同化效果好于单一资料同化。     

Effects of Different Initial Fields on GRAPES Numerical Prediction

In this paper,a heavy rainfall process occurring in the Huaihe River Basin during 9-10 July 2005 is studied by the new generation numerical weather prediction model system-GRAPES,from the view of different initial field effects on the prediction of the model.Several numerical experiments are conducted with the initial conditions and lateral boundary fields provided by T213 L31 and NCEP final analyses,respectively. The sensitivity of prediction products generated by GRAPES to different initial conditions,including effects of three-dimensional variational assimilation on the results,is discussed.After analyzing the differences between the two initial fields and the four simulated results,the memonic ability of the model to initial fields and their influences on precipitation forecast are investigated.Analyses show the obvious differences of sub-synoptic scale between T213 and NCEP initial fields,which result in the corresponding different simulation results,and the differences do not disappear with the integration running.It also shows that for the same initial field whether it has data assimilation or not,it only obviously influences the GRAPES model results in the initial 24 h.Then the differences reduce.In addition,both the Iocation and intensity of heavy rain forecasted by GRAPES model Further is very close to the fact,but the forecasting area of strong torrential rain has some differences from the fact.For the same initial field when it has assimilation, the 9-12-,12-24-,and 0-24-h precipitation forecasts of the model are better than those without assimilation. All these suggest that the ability of GRAPES numerical prediction depends on the different initial fields and lateral boundary conditions to some extent,and the differences of initial fields will determine the differences of GRAPES simulated results.     

地形重力波拖曳参数化对热带气旋强度和路径预报影响的研究

本文在GRAPES_TMM（Global/Regional Assimilation and Prediction System for Tropical Mesoscale Model）——中国南海台风模式版（面向南海和东南亚）中发展和引进了KA95（Kim and Arakawa,1995）地形重力波拖曳参数化方案（GWDO）,并对2012年主要的9个登陆台风进行了试验对比研究,考察了不同标准Richardson数（Ri_c）的GWDO试验对台风路径和强度预报的影响。结果表明,在引入地形重力波拖曳参数化过程后,模式对台风登陆时路径和强度的预报能力均要有提高,对台风预报时长越长,GWDO的影响也更为显著。对双台风“SAOLA”和“DAMREY”试验结果表明,GWDO对台风外围距台风中心150 km的对流层中下层风速减弱较为明显,减弱了GRAPES区域模式对台风强度预报偏强的现象,对台风强度长时间预报改善更为明显。不同标准Ri_c对重力波拖曳力的计算较为敏感,当Ri_c取1.0时,动能迅速的在低层被频散,能量无法有效地上传;Ri_c取0.25时,大部分的能量在中高层被频散。总的来说,Ri_c取0.75时对台风路径和强度预报改进更为显著,其结果可为业务预报提供指导意义。     

Study on ETKF-Based Initial Perturbation Scheme for GRAPES Global Ensemble Prediction

Initial perturbation scheme is one of the important problems for ensemble prediction. In this paper, ensemble initial perturbation scheme for Global/Regional Assimilation and PrEdiction System (GRAPES) global ensemble prediction is developed in terms of the ensemble transform Kalman filter (ETKF) method.A new GRAPES global ensemble prediction system (GEPS) is also constructed. The spherical simplex 14-member ensemble prediction experiments, using the simulated observation network and error characteristics of simulated observations and innovation-based in ation, are carried out for about two months. The structure characters and perturbation amplitudes of the ETKF initial perturbations and the perturbation growth characters are analyzed, and their qualities and abilities for the ensemble initial perturbations are given. The preliminary experimental results indicate that the ETKF-based GRAPES ensemble initial perturbations could identify main normal structures of analysis error variance and reflect the perturbation amplitudes.The initial perturbations and the spread are reasonable. The initial perturbation variance, which is approximately equal to the forecast error variance, is found to respond to changes in the observational spatial variations with simulated observational network density. The perturbations generated through the simplex method are also shown to exhibit a very high degree of consistency between initial analysis and short-range forecast perturbations. The appropriate growth and spread of ensemble perturbations can be maintained up to 96-h lead time. The statistical results for 52-day ensemble forecasts show that the forecast scores ofensemble average for the Northern Hemisphere are higher than that of the control forecast. Provided that using more ensemble members, a real-time observational network and a more appropriate inflation factor,better effects of the ETKF-based initial scheme should be shown.     

对流天气系统自动站雨量资料同化对降雨预报的影响

利用GRAPES(Global and Regional Assimilation and Prediction Enhanced System,全球/区域同化预报系统)三维变分同化系统,针对对流天气系统特点,用改进的郭晓岚对流参数化方案作为观测算子,同化广东省自动站记录的对流天气系统的雨量资料,并且与同化探空资料进行了比较.在雨带有明显改进的区域,分别同化这两种资料都可以调整大气低层水汽辐合增加(或辐散),对流层中下层增暖增湿(或变冷变干),从而增加(或减少)降水,表明降水的同化方案对初始场的调整在一定     

GFS物理过程包在GRAPES区域模式中的实施及改进:单柱试验

将NCEP的中期全球预报系统（GFS）整套物理过程包引入GRAPES区域模式,并利用美国ARM计划（Atmospheric Radiation Measurement Program）南部大平原外场1997年夏加强观测资料,对GRAPES模式中原有物理过程和GFS物理过程进行了单柱试验的测试和比较,评估了两组物理过程的模拟性能。单柱试验结果表明GFS试验模拟的位温和水汽混合比演变更为准确;到地短波和长波辐射的误差更小;其模拟的地面向上潜热通量与观测更为接近,而原物理过程试验的感热偏大;GFS试验对2 m气温最高和最低值的模拟也比原有物理过程更接近观测。进一步分析降水的模拟发现,两组物理过程对强降水过程均有很好的预报能力,但对中等强度降水存在漏报或者发生时间偏晚的倾向,此外还存在对小雨的空报。在3个降水事件的模拟中,GFS试验模拟的降水误差小于原物理过程试验。最后还对两组试验中明显的降水过程漏报和降水延迟进行了分析,对GFS物理过程包中对流参数化的对流触发方式的修改,可改进此次降水的模拟,进而通过更准确的反馈过程,改善对大尺度场的模拟。     

Possible sources of forecast errors generated by the global/regional assimilation and prediction system for landfalling tropical cyclones. Part I: Initial uncertainties

This paper investigates the possible sources of errors associated with tropical cyclone (TC) tracks forecasted using the Global/Regional Assimilation and Prediction System (GRAPES). The GRAPES forecasts were made for 16 landfalling TCs in the western North Pacific basin during the 2008 and 2009 seasons, with a forecast length of 72 hours, and using the default initial conditions (“initials”, hereafter), which are from the NCEP-FNL dataset, as well as ECMWF initials. The forecasts are compared with ECMWF forecasts. The results show that in most TCs, the GRAPES forecasts are improved when using the ECMWF initials compared with the default initials. Compared with the ECMWF initials, the default initials produce lower intensity TCs and a lower intensity subtropical high, but a higher intensity South Asia high and monsoon trough, as well as a higher temperature but lower specific humidity at the TC center. Replacement of the geopotential height and wind fields with the ECMWF initials in and around the TC center at the initial time was found to be the most efficient way to improve the forecasts. In addition, TCs that showed the greatest improvement in forecast accuracy usually had the largest initial uncertainties in TC intensity and were usually in the intensifying phase. The results demonstrate the importance of the initial intensity for TC track forecasts made using GRAPES, and indicate the model is better in describing the intensifying phase than the decaying phase of TCs. Finally, the limit of the improvement indicates that the model error associated with GRAPES forecasts may be the main cause of poor forecasts of landfalling TCs. Thus, further examinations of the model errors are required.     

CCM3模式中LSM积雪方案的改进研究(Ⅰ):修改方案介绍及其单点试验

为了改进美国NCARCCM3全球模式中LSM陆面模型中的积雪方案的模拟效果,在Sun等[1]SAST积雪模型的基础上,作了部分修改后,加进CCM3模式LSM模型中.该方案根据格点区域平均积雪深度的不同,把地面雪盖划分为1到3层不等,能在积雪表层和中间层更好地描述温度的日变化和季节变化;较详细地考虑了雪的热传导、太阳辐射的穿透吸收、雪的融化、液态水的储存、渗透和再冻结等积雪内部的主要物理过程;根据Nimbus-7卫星实测雪深资料修改了积雪覆盖度和雪面反照率的计算方案.利用前苏联6个台站1978-1983年的实测积雪资料和大气强迫数据,进行了单点模拟试验,结果表明,新的积雪参数化方案能够较好地再现积雪深度和雪水当量的逐日和季节变化特征,部分提高了积雪参数化方案对积雪的模拟能力.     

The GRAPES Variational Bias Correction Scheme and Associated Preliminary Experiments

The variational assimilation theory is generally based on unbiased observations. In practice, however, almost all observations suffer from biases arising from observational instruments, radiative transfer operator, precondition of data, and so on. Therefore, a bias correction scheme is indispensable. The current scheme for radiance bias correction in the GRAPES 3DVar system is an offline scheme. It is actually a static correction for the radiance bias before the process of cost function minimization. In consideration of its effects on forecast results, this kind of scheme has some shortcomings. Thus, this study provides a variational bias correction (VarBC) scheme for the GRAPES 3DVar system following Dee’s idea. In the VarBC scheme, the observation operator is modified and a new control variable is defined by taking the predictor coefficients as the control parameters. According to the feature of the GRAPES-3DVAR, an incremental formulation is applied and the original bias correction scheme is maintained in the actual process of observations. The VarBC is designed to co-exist with the original scheme, because it is a dynamic revision to the observational operator on the basis of the old method, i.e., it adjusts the model state vector along with the control parameters to an unbiased state in the process of minimization and the assimilation system remains consistent with available information automatically. Preliminary experimental results show that the mean departures of background-minus-observation and analysis-minus-observation are reduced as expected. In a case study of the heavy rainfall that happened in South China on 11–13 June 2008, the 500-hPa geopotential height is better simulated using the analyzed field from the VarBC as the initial condition.     

气象观测设备测试与试验系统设计与实现

为了实现气象观测设备测试评估与试验业务的规范化、标准化和自动化处理,设计实现针对常规地面气象观测设备测试评估与试验的业务系统。系统基于主流Web开发技术和Oracle数据库,以灵活配置、分层设计、模块化应用的方式,实现了测试评估与试验业务中的流程管理、数据配置、设备运行及状态管理、专家分析计算、数据查询显示、测评报告生成导出和系统管理等功能模块,实现了测试评估与试验业务的统一规范管理、提高了工作效率,并在雪深观测仪等设备的测试评估与试验中得到应用,计算结果和评估结论与人工处理一致。后期将以此为基础,进一步设计开发探空设备、雷达设备等适用的测试评估与试验平台系统。