GRAPES伴随模式底层数据栈优化

GRAPES伴随模式是其四维变分同化系统的核心组成部分。由于其计算过程复杂,临时数据较多,实现中采用断点存储策略可以有效减少伴随模式的计算时间和存储空间。极限断点存储策略是在单积分步内以全存储策略实现为基础,将其中部分基态以计算代替的一种类断点存储策略。在该策略的支持下,需要一种新的数据管理结构,来保证程序的正确运行。文章提出了在已有栈基础上优化的新数据存储管理方式——嵌套多链栈,这种结构可有效满足使用极限断点存储技术实现GRAPES伴随模式的初态管理需求。试验表明:相比断点存储技术,在总内存增加不超过30%的情况下可使GRAPES的运行效率提高1倍。     

GRAPES非静力数值预报模式的三维变分资料同化系统的发展

为了减少分析变量与模式状态变量之间的插值误差,改善业务预报模式的初值质量,在GRAPES等压面三维变分资料同化系统的基础上,研究发展了针对GRAPES区域模式的非静力模式变量三维变分资料同化系统(GRAPES m3DVAR).该资料同化系统的垂直坐标及其分析变量的水平分布格式、垂直跳点方案与GRAPES预报模式保持完全一致.由于垂直坐标的变化和非静力关系,m3DVAR分析系统中设计了求解动力学约束方程的新方案.通过有效的高精度数学方案,避免了地形追随坐标下平衡方程的非线性项造成的复杂计算,有效解决了非静力平衡条件下求解平衡方程中非线性项的切线性方程和伴随方程引起的困难.重新构造各种观测算子,并考虑了质量场和风场之间的平衡约束关系、背景误差协方差结构,实现对探空、地面资料、船舶报等常规观测的同化.理想单点试验和实际资料的多变量资料同化分析结果表明,非静力模式变量三维变分资料同化系统能够正确地描写多变量之间的相互作用以及物理约束关系,分析结果合理,能够有效减少原等压面三维变分资料同化系统的分析与模式变量之间需要相互插值、变换产生的误差,在一定程度上提高了分析场质量,对预报模式的初值具有一定改善.     

GRAPES模式切线性垂直扩散方案的误差分析和改进

针对GRAPES四维变分同化系统的升级, 研究了GRAPES模式垂直扩散方案线性化问题。通过2005年8月7—27日21个个例的批量试验, 发现在GRAPES模式垂直扩散方案源代码的基础上逐句线性化得到的切线性垂直扩散方案即使能通过正确性测试试验, 在少数情况下也会存在很大误差。切线性模式计算的扰动气压场和扰动风场可能出现明显异常, 这种异常与垂直扩散方案中地表动量通量的强非线性有关。如果在切线性垂直扩散方案中忽略地表动量通量扰动, 既可以避免异常的出现, 又不影响其他正常时刻的计算精度。修改后的切线性垂直扩散方案能够在所有变量上一致地提高切线性模式的计算精度。     

复杂地形下地面观测资料同化I. 模式地形与观测站地形高度差异对地面资料同化的影响研究

中国地形复杂，模式地形与实际观测地形存在一定高度差异，因此设计合理的复杂地形下地面观测资料的同化方案有利于使我国目前仅用作探测手段的地面观测资料（常规地面观测站和地面自动站）在中尺度数值模式中得到充分利用。作者在MM5_3DVAR同化系统中利用近地层相似理论将地面观测资料进行直接三维变分同化分析，并对地面资料同化方案设计中是否需要考虑模式与实际观测站地形高度差异进行探讨研究。研究结果表明：通过近地层相似理论将地面观测资料同化到数值模式能起到一定的作用，并且地面观测资料（温度、 湿度、 风场、 地面气压）中各物理量同化到数值模式都能影响24小时降水数值结果，但各物理量起的作用大小不一样，其中影响最大的是温度，其次为湿度；地面观测资料同化方案设计有必要考虑模式地形与实际观测站地形高度差异，适当考虑这种高度差异能取得较好的结果。     

Experiments of the WRF three-/four-dimensional variational (3/4DVAR) data assimilation in the forecasting of Antarctic cyclones

The three-/four-dimensional variational data assimilation systems (3/4DVAR) of the Weather Research and Forecasting (WRF) model were explored in the forecasting of two Antarctic synoptic cyclones, which had large influence on the Ross Sea/Ross Ice Shelf region in October 2007. A suite of variational data assimilation experiments, including regular 3DVAR, high-resolution 3DVAR, and 4DVAR experiments, were designed to evaluate their performances in weather analysis and forecasting in Antarctica. In general, both 4DVAR and high-resolution 3DVAR experiments showed better forecasting skill than regular 3DVAR experiments. High-resolution 3DVAR experiments were the most efficient in reducing the analysis errors of surface winds and temperature, and had the best performance during the first 24 h of forecasting. However, during the following forecast period, 4DVAR experiments showed either better or about comparable performance to high-resolution 3DVAR experiments. These results indicate that increasing the spatial resolution during 3DVAR is an economical approach to improving the weather analysis and forecasting over Antarctica. At the same time, the 4DVAR approach had a longer impact on forecasting than the high-resolution 3DVAR approach. Understandably, both of the variational assimilation approaches are promising techniques toward improving the regional analysis and forecasting over Antarctica.     

在四维变分同化中运用集合协方差的试验

利用浅水方程模式和模式模拟资料进行数值试验比较3种不同的背景误差协方差矩阵处理方法对四维变分(4DVAR)资料同化的影响.3种背景误差协方差矩阵分别是:(1)对单一变量将背景误差协方差矩阵简化为对角矩阵;(2)将背景误差协方差矩阵的作用简化为高斯过滤;(3)由预报集合生成背景误差协方差矩阵并利用奇异值分解技术解决矩阵的求逆.通过一系列数值试验,比较不同观测密度、不同观测误差下3种背景误差协方差处理方法对4DVAR同化效果的影响.结果表明,背景误差协方差的结构对4DVAR有重大影响.当观测资料的空间密度不够高时,采用对角矩阵得不到满意的结果.高斯过滤方案可以明显改善同化结果,但是对背景误差特征长度比较敏感.第3种方法采用的背景误差协方差矩阵是流型依赖的,而且并不以显式的方式出现在目标函数中.避免了对它求逆的复杂运算.由于做了降维处理,在观测点的密度较低和观测误差较大时可望取得较好的同化结果,同化效果较为稳定.     

复杂地形下地面观测资料同化II. 模式地形与观测站地形高度差异代表性误差

第I部分研究结果(徐枝芳等,2007)表明模式与实际观测站地形高度差异对地面观测资料同化效果有较大影响。此文在MM5_3DVAR同化系统中利用近地层相似理论将地面观测资料进行直接三维变分同化分析,考虑模式与实际观测站地形高度差异对同化效果的影响,提出在地面观测误差中增加地形代表性误差来解决这个问题。研究结果表明:地面资料同化分析时,在其观测误差中加入一项新的误差——地形代表性误差,能较好地解决地面资料同化分析中模式与观测站地形高度差异问题;地面资料参与同化分析,在观测误差中加入与模式和实际观测站地形高度差异大小相关的地形代表性误差时,地面观测值对分析值的影响随着地形高度差异代表性误差的加入而减小,同时又部分地将地面观测信息通过变分分析融进分析场,使得低层分析更接近真实场,且地面资料利用率更高,24小时降水数值预报(模拟)的效果较好。     

GRAPES混合同化方法在青藏高原区域的初步试验

针对青藏高原地区气象观测站点稀少和模式同化分析质量较低的问题,将GRAPES区域集合变分(En-3DVAR)混合同化系统应用于青藏高原地区,进行了单点理想试验和真实观测资料同化分析预报试验,分析评估青藏高原混合同化分析增量及预报误差的水平垂直结构特征及其合理性,并与中国东部平原地区进行比对。单点理想试验表明,En-3DVAR混合同化系统中背景误差协方差具备流依赖属性。真实资料混合同化试验结果表明,基于集合预报估计的分析增量具有流依赖的特征,特别是在高原地区和槽脊系统附近;青藏高原地区分析场的绝对误差总体低于3DVAR系统,其平均绝对误差在中层和高层低于平原地区,说明在青藏高原地区的改进效果略优于平原地区。需要关注的是,青藏高原地区En-3DVAR混合同化分析增量总体大于3DVAR,特别是近地面层u风分量分析增量明显偏大,这可能与青藏高原复杂地形有关。     

温湿统计平衡约束关系对GRAPES全球湿度分析的作用

为改进GRAPES全球三维变分同化系统(GRAPES-3DVar)的湿度分析,借鉴Hólm等(2002)的思想,在背景误差协方差结构中引入湿度与温度的统计平衡约束关系。通过扣除湿度变化中与温度有关的平衡部分获取非平衡拟相对湿度,并引入非线性对称变换对其做标准化处理,将处理后的变量作为新的湿度控制变量。统计结果表明,温湿统计平衡约束主要出现在中高纬度对流层中层相对湿度大于80%的区域,与大尺度抬升凝结加热有关;新的湿度控制变量能满足无偏、高斯分布特征。单点理想观测试验结果表明,新的湿度分析具备了流依赖特征,并能有效地抑制负水汽与超饱和水汽的出现。同化循环与预报试验结果表明,新方案给出的湿度分析的偏差和均方根误差均有所减小。而针对降水预报的检验结果表明,引入新方案后的0.1-10 mm降水预报,在ETS评分没有显著降低的情况下,BIAS评分更靠近1,降水空报有所减缓。然而60-84 h的25 mm以上的降水漏报现象更为明显,表明湿度同化分析方案还有改进空间。通过引入温湿统计平衡约束关系,完善了GRAPES-3DVar分析框架,为全球湿度分析的持续改进奠定了坚实基础。      

基于高光谱大气红外探测器模拟亮温的稳健变分反演研究

经典变分反演法是基于观测误差服从高斯分布的假定,对偏离均值较大的离群值较敏感。当实际观测数据包含离群值观测误差呈现非高斯分布时,如果采用经典反演法进行变分反演就会产生大的偏差,甚至导致变分反演的失败。使用经典变分反演法首先需要进行质量控制,剔除所谓的离群值,但有相当部分的离群值包含了一些亮点,如天气现象。如果对其“视而不见”,则对很多重要的信息就无法把握。基于此,研究采用稳健变分反演的思想同化这些离群值,主要思想是把M-估计法（L2、Huber、Fair和Cauchy一估计）的权重函数耦合到经典变分反演中,在每次变分反演极小化迭代过程中重新估计观测项对经典变分反演目标泛函的贡献率。采用高光谱大气红外探测器（AtmosphericInfraRedSounder,AIRS）的通道模拟亮温进行理想试验,结果表明：采用Huber一估计进行稳健变分反演对温度和湿度反演具有较好的效果;采用Cauchy-估计得到的效果反而更差,这是由Cauchy等分布函数固有的缺陷所决定的。因此,稳健变分反演观测误差非高斯分布是可行的,但依赖M-估计法权重函数的选取。     

Coupling Ensemble Kalman Filter with Four-dimensional Variational Data Assimilation

This study examines the performance of coupling the deterministic four-dimensional variational assimilation system (4DVAR) with an ensemble Kalman filter (EnKF) to produce a superior hybrid approach for data assimilation. The coupled assimilation scheme (E4DVAR) benefits from using the state-dependent uncertainty provided by EnKF while taking advantage of 4DVAR in preventing filter divergence: the 4DVAR analysis produces posterior maximum likelihood solutions through minimization of a cost function about which the ensemble perturbations are transformed, and the resulting ensemble analysis can be propagated forward both for the next assimilation cycle and as a basis for ensemble forecasting. The feasibility and effectiveness of this coupled approach are demonstrated in an idealized model with simulated observations. It is found that the E4DVAR is capable of outperforming both 4DVAR and the EnKF under both perfect- and imperfect-model scenarios. The performance of the coupled scheme is also less sensitive to either the ensemble size or the assimilation window length than those for standard EnKF or 4DVAR implementations.     

多普勒激光雷达风场反演方法研究

采用三维变分同化反演(3DVAR)、 四维变分同化反演(4DVAR) 对多普勒激光雷达资料反演风场的方法进行了研究, 利用车载多普勒激光雷达在2008年残奥会测试赛期间外场试验取得的数据, 反演了海面10 m高度处的风场, 并将风场反演结果与浮标资料进行了对比分析, 结果表明: 3DVAR、4DVAR风场反演方法均能实现近海面风的精细化风场反演, 并能反映出风向的变化, 反演风场与浮标数据基本一致, 在风速较大的天气情况, 3DVAR与4DVAR反演风场的一致性要好于风速较小的天气情况; 4DVAR反演方法中以浮标资料作为背景场, 使得其与浮标的符合程度要好于3DVAR方法反演风场; 反演风场的风向与浮标风向具有很好的相关关系, 反演风场的风速与浮标风速具有一定的相关关系, 反演风场的风向、风速与浮标的风向、风速之间平均均方根误差和平均绝对误差表明, 这两序列之间具有一定差别, 在风速较小的天气情况下使用时需要注意。     

Radar Data Assimilation for the Simulation of Mesoscale Convective Systems

A heavy rainfall case related to Mesoscale Convective Systems (MCSs) over the Korean Peninsula was selected to investigate the impact of radar data assimilation on a heavy rainfall forecast. The Weather Research and Forecasting (WRF) three-dimensional variational (3DVAR) data assimilation system with tuning of the length scale of the background error covariance and observation error parameters was used to assimilate radar radial velocity and reflectivity data. The radar data used in the assimilation experiments were preprocessed using quality-control procedures and interpolated/thinned into Cartesian coordinates by the SPRINT/CEDRIC packages. Sensitivity experiments were carried out in order to determine the optimal values of the assimilation window length and the update frequency used for the rapid update cycle and incremental analysis update experiments. The assimilation of radar data has a positive influence on the heavy rainfall forecast. Quantitative features of the heavy rainfall case, such as the maximum rainfall amount and Root Mean Squared Differences (RMSDs) of zonal/meridional wind components, were improved by tuning of the length scale and observation error parameters. Qualitative features of the case, such as the maximum rainfall position and time series of hourly rainfall, were enhanced by an incremental analysis update technique. The positive effects of the radar data assimilation and the tuning of the length scale and observation error parameters were clearly shown by the 3DVAR increment.     

降水资料同化在GRAPES-MESO模式中应用试验研究

利用国家气象中心中尺度业务数值预报模式GRAPES-MESO v3.0,以2010年6月1～30日为例,开展地面降水率1DVAR(one-dimensional variational assimilation)同化方案在GRAPES-3DVAR(three-dimensional variational assimilation)同化系统中的应用试验研究(ASSI试验),并以未加降水资料同化的试验为对照试验(CNTL试验),以评估全国1h加密雨量资料在模式中同化应用的效果。结果表明:1)在相对湿度背景误差和降水率观测误差范围内,1DVAR同化方案能够对湿度廓线进行有意义的调整,使分析降水向观测降水靠近;ASSI试验对初始温、压、湿、风场的修正主要为正效果;2)对2010年6月17～21日江南、华南连续性降水过程进行了分析,整体而言ASSI试验对逐日及逐时降水强度的预报普遍强于CNTL试验,与实况更加接近;3)ASSI试验对2010年6月1～30日08时起报的0～24 h模式预报的小雨、中雨、大雨、暴雨、大暴雨各个降水量级TS评分及ETS评分相比CNTL试验均有较明显提高,预报偏差也更接近于1;4)ASSI试验较CNTL试验能更好地模拟雨带的分布、雨带演变特征和降水强度的变化;5)对降水所做的典型个例和统计检验分析从不同角度说明了地面降水资料1DVAR同化方案在GRAPES-3DVAR系统中的应用改善了GRAPES-MESO v3.0的降水模拟效果。     

适应数值模式垂直坐标的GPS掩星反演资料稀疏化方案

GPS掩星反演大气温湿资料具有高垂直分辨率、高精度、受云和降水影响小等优点,针对GRAPES同化预报系统,发展设计了一种既考虑预报模式高度-地形追随垂直坐标不均匀分层特点,又结合掩星反演资料特性的新适应性垂直稀疏化方案。通过个例试验和批量试验,探索了该适应性稀疏化方案对分析预报质量的影响。试验结果表明:选取合适的稀疏参数,新稀疏化方案的GPS掩星反演资料对背景场的调整更加有效,分析场质量更高;位势高度、比湿、温度和风场等预报场的均方根误差均更小;适应性稀疏化方案对改善台风路径预报具有积极作用。批量试验则进一步证实了适应性稀疏化方案对分析场质量有明显的改善作用。     

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

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

Assimilation of Global Navigation Satellite Radio Occultation Observations in GRAPES: Operational Implementation

This paper presents the design of an observation operator for assimilation of global navigation satellite system (GNSS) radio occultation (RO) refractivity and the related operational implementation strategy in the global GRAPES variational data assimilation system. A preliminary assessment of the RO data assimilation effect is performed. The results show that the RO data are one of the most important observation types in GRAPES, as they have a significant positive impact on the analysis and forecast at all ranges, especially in the Southern Hemisphere and the global stratosphere where in-situ measurements are lacking. The GRAPES model error cannot be controlled in the Southern Hemisphere without RO data being assimilated. In addition, it is found that the RO data play a key role in the stable running of the GRAPES global assimilation and forecast system. Even in a relatively simple global data assimilation experiment, in which only the conventional and RO data are assimilated, the system is able to run for more than nine months without drift compared with NCEP analyses. The analysis skills in both the Northern and Southern Hemispheres are still relatively comparable even after nine-month integration, especially in the stratosphere where the number of conventional observations decreases and RO observations with a uniform global coverage dominate gradually.     

Assimilation of Global Navigation Satellite Radio Occultation Observations in GRAPES:Operational Implementation

This paper presents the design of an observation operator for assimilation of global navigation satellite system(GNSS) radio occultation(RO) refractivity and the related operational implementation strategy in the global GRAPES variational data assimilation system.A preliminary assessment of the RO data assimilation effect is performed.The results show that the RO data are one of the most important observation types in GRAPES,as they have a significant positive impact on the analysis and forecast at all ranges,especially in the Southern Hemisphere and the global stratosphere where in-situ measurements are lacking.The GRAPES model error cannot be controlled in the Southern Hemisphere without RO data being assimilated.In addition,it is found that the RO data play a key role in the stable running of the GRAPES global assimilation and forecast system.Even in a relatively simple global data assimilation experiment,in which only the conventional and RO data are assimilated,the system is able to run for more than nine months without drift compared with NCEP analyses.The analysis skills in both the Northern and Southern Hemispheres are still relatively comparable even after nine-month integration,especially in the stratosphere where the number of conventional observations decreases and RO observations with a uniform global coverage dominate gradually.     

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观测资料同化效果提供了有力工具。      

Simulation of heavy rainfall events over Indian monsoon region using WRF-3DVAR data assimilation system

We present the results of the impact of the 3D variational data assimilation (3DVAR) system within the Weather Research and Forecasting (WRF) model to simulate three heavy rainfall events (25–28 June 2005, 29–31 July 2004, and 7–9 August 2002) over the Indian monsoon region. For each event, two numerical experiments were performed. In the first experiment, namely the control simulation (CNTL), the low-resolution global analyses are used as the initial and boundary conditions of the model. In the second experiment (3DV-ANA), the model integration was carried out by inserting additional observations in the model’s initial conditions using the 3DVAR scheme. The 3DVAR used surface weather stations, buoy, ship, radiosonde/rawinsonde, and satellite (oceanic surface wind, cloud motion wind, and cloud top temperature) observations obtained from the India Meteorological Department (IMD). After the successful inclusion of additional observational data using the 3DVAR data assimilation technique, the resulting reanalysis was able to successfully reproduce the structure of convective organization as well as prominent synoptic features associated with the mid-tropospheric cyclones (MTC). The location and intensity of the MTC were better simulated in the 3DV-ANA as compared to the CNTL. The results demonstrate that the improved initial conditions of the mesoscale model using 3DVAR enhanced the location and amount of rainfall over the Indian monsoon region. Model verification and statistical skill were assessed with the help of available upper-air sounding data. The objective verification further highlighted the efficiency of the data assimilation system. The improvements in the 3DVAR run are uniformly better as compared to the CNTL run for all the three cases. The mesoscale 3DVAR data assimilation system is not operational in the weather forecasting centers in India and a significant finding in this study is that the assimilation of Indian conventional and non-conventional observation datasets into numerical weather forecast models can help improve the simulation accuracy of meso-convective activities over the Indian monsoon region. Results from the control experiments also highlight that weather and regional climate model simulations with coarse analysis have high uncertainty in simulating heavy rain events over the Indian monsoon region and assimilation approaches, such as the 3DVAR can help reduce this uncertainty.