一种新的资料同化方法

为寻求一种快速有效的四维变分资料同化(英文缩写4DVar)作了有意义的尝试,提出了映射观测的新概念和反向四维变分资料同化的新思路,并以此为基础建立了三维变分映射资料同化(英文缩写为3DVM:3-DimensionalVariational data assimilation of Mapped observation)。该方法与传统的四维变分资料同化一样,不仅考虑了模式的动力和物理约束,使得同化后的初值与模式协调,而且通过模式方程对同化窗口中不同时刻的观测资料作了最佳拟合。与传统四维变分同化方法不同的是,由3DVM得到的初值不在同化窗口的始端,而在窗口的末端。正是所求初值时刻的改变,使得该方法的计算代价大大减少,几乎与三维变分资料同化(英文缩写3DVar)相当,这实际上是用3DVar的代价实现了4DVar的功能。同时,由于3DVM不再需要切线性和伴随近似来计算代价函数的梯度也提高了同化的精度。对具体的台风个例(Dan)用AMSU-A反演的温度场进行变分同化模拟试验,发现3DVM能比传统4DVar产生更好的初值,而且所花计算时间只需4DVar的1/7。     

利用ARGO资料改进海洋资料同化和海洋模式中的物理过程

国际ARGO计划的实施每年将可提供多达 10万个剖面 (0～ 2 0 0 0m水深 )的海水温度和盐度资料 ,这些资料的获取无疑将会大大促进海洋和大气科学的发展 ,使人们加深对海洋过程的了解 ,揭示海 气相互作用的机理 ,为长期天气预报和短期气候预测提供模式初始场 ,提高长期天气预报和短期气候预测的能力。如何利用这些资料开展研究工作以及在实际业务中应用这些资料是目前大气和海洋科学界的一个前沿课题。本研究将ARGO浮标资料引入了国家气候中心的NCC GODAS同化系统 ,结果分析表明 ,同化ARGO资料后所得到的海温场在三大洋中不仅在温度数值的大小 ,而且在分布形式方面都与观测场具有较好的一致性 ,可以很好地反映出观测到的冬季和夏季海温的分布形式以及海温的季节变化特征和异常特征。本研究还应用最新的ARGO海洋观测资料 ,通过建立新的热带西太平洋次表层海温参数化方案 ,改进了Zebiak Cane(1987)海洋模式 (ZC模式 ) ,克服了ZC模式几乎没有模拟赤道西太平洋表层和次表层海温变化能力的缺陷。在ZC模式中引用新的次表层海温参数化方案后 ,在赤道西太平洋不仅次表层海温的模拟得到了改善 ,对海面温度异常的模拟也有了较大的改进 ,不仅模拟出了赤道西太平洋表层和次表层海温异常的年际变化特征 ,也模拟出了与观     

全球海洋资料同化系统ZFL_GODAS的研制和初步评估试验

本研究发展了一个全球海洋资料同化系统ZFL_GODAS。该系统是一个短期气候数值预测业务系统的子系统,为短期气候预测海气耦合模式提供全球海洋初始场。系统能够同化的观测资料包括卫星高度计资料、卫星海表温度(SST)资料,以及Argo、XBT、TAO等各种不同来源的现场温盐廓线资料。系统使用的海洋模式为中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室开发的气候系统海洋模式LICOM1.0,同化方案为集合最优插值(EnOI)方案。系统使用一个由海洋模式自由积分得到的静态样本来估计背景场误差协方差。这样的基于集合样本的背景场误差协方差具有多变量协变、各向异性的特征,且能反映海洋物理过程固有的空间尺度特征。针对EnOI同化程序的特点,开发了一套特色鲜明、负载均衡、高效的并行化同化程序。本文通过与不同类型观测资料的比较,对同化系统的性能进行了评估。通过比较海表温度和海面高度的年际变率,海表温度异常随时间的变化,SST、海面高度异常(SLA)以及次表层温盐预报产品的均方根误差,5年平均温度偏差廓线、平均盐度廓线、平均纬向流速廓线等发现:系统工作正常、同化效果较好;经过同化以后,各变量都更加接近观测,误差更小,与观测场的相关性更好,可以为短期气候预测系统提供较好的海洋初始场,也可以为物理海洋学的研究提供有效的再分析资料。     

从GTS获得的海洋温、盐资料在BCC海洋同化系统中的质量控制及同化结果初步分析

介绍了从全球电信系统(GTS)上获得的海洋温度、盐度观测资料在中国国家气候中心(BCC)新一代海洋同化系统中的应用情况。通过资料的质量控制判断温、盐观测的重复记录、观测深度、地形、极端值、气候变率、层结、空间差异,有效地过滤了错误的或不可靠的观测信息。质量控制后,将温、盐观测资料加入同化系统,有效地改进了模块化全球海洋环流模式MOM4中的全球热带、副热带海洋,尤其是太平洋地区的多年平均海表温度、盐度场分布特征;此外,同化温、盐资料对南北半球中纬度地区的海表温度分布特征也有明显的改进。对比同化前后的均方根误差(RMSE)发现,同化后大部分海区,尤其是热带海洋的海表温度/盐度的均方根误差明显降低,降幅通常在0.1—1.0℃/psu,模拟与观测的海表温、盐分布特征也更为接近。进一步分析指出,同化明显地改善了模式对Nino3、Nino4区海温时间演变特征的模拟,同化后的Nino3海温与最优插值海表温度的差异减小,但其通常在上半年改进较多(差值绝对值多在0.5℃左右),而在下半年则改进较少(差值绝对值常达1℃左右);Nino4区的海温特征则改进明显,其与最优插值海表温度的差值绝对值通常都控制在0.5℃以下。     

GRAPES集合卡尔曼滤波资料同化系统Ⅱ：区域分析及集合预报

GRAPES集合卡尔曼滤波资料同化方法能够分批同化常规观测资料,GRAPES集合卡尔曼滤波同化系统的设计及其与GRAPES三维变分同化系统的对比试验结果表明,GRAPES集合卡尔曼滤波系统能够得到合理的分析,并且具有实际运行能力。在此基础上,进行集合卡尔曼滤波区域同化分析及集合预报试验,对比区域模式面三维变分同化分析预报结果,研究表明,集合卡尔曼滤波分析比三维变分分析具有一定优势,降水预报更接近实况。考察了预报误差特征随天气形势的变化情况,表明预报误差相关场和均方差的分布随着天气形式不同而变化。     

四维变分同化技术在优化Z-C预测模式初始场中的应用试验

文中使用四维变分同化技术将海温观测资料同化到Zebiak-Cane模式中,通过优化模式的初始场提高了模式的预报技巧.通过用理想场进行检验,说明所建立的同化伴随模式是正确的 .用文中建立的四维变分同化模式以1997年1月为初始场所做的预报结果与实况相比,结果较好.这对今后ENSO预报打下了良好的基础.     

海洋模式对1995年5月南海AXBT海温资料的同化分析

利用美国普林斯顿大学海洋模式 (POM )对 1995年 5月南海AXBT海洋观测资料和MISST卫星观测的海面温度进行模式四维同化分析研究。POM模式同化处理使用最优插值分析与逼近方法 ,从模式积分到第 4年开始 ,使用COADS气候平均风应力场、净热通量场强迫海洋 ,从第 5年 5月开始加入 1995年 5月的MISST和AXBT资料同化积分 30d ,然后对同化结果进行分析。结果表明 ,资料同化可以改进模式对温度场的预报 ,也可改进模式对海流的模拟。分析结果还显示 ,1995年 5月南海浅层有暖涡和冷涡形成     

散射计风场的三维变分对海雾数值模拟的影响

为了评估同化QuikSCAT海面风场资料对海雾模拟的影响,对发生在2006年4月3—5日的海雾过程,首先通过3个敏感性试验,研究了不同边界层参数化方案对海雾模拟的影响,发现YSU边界层参数化方案更适合海雾过程的模拟。然后对2006年4月3—5日平流雾过程和2005年6月23—24日辐射雾过程利用WRF-3DVAR系统将QuikSCAT海面风场资料同化到模式中,并以未同化和同化了QuikSCAT海面风场资料的数据为初始场,应用WRF模式进行模拟预报,同时对模拟预报得到的结果与实况 (卫星云图和地面观测) 进行对比,结果表明:QuikSCAT海面风场资料的三维变分同化能够改善低层其他要素场,对海雾预报有明显的正效应,但对高层的影响相对有限。     

非高斯分布观测误差资料的变分质量控制对暴雨预报的影响

观测资料的质量控制直接影响数值预报资料同化的分析质量。本文针对GRAPES区域同化和预报的三维变分资料同化系统,发展了基于观测误差为"高斯分布+均匀分布"模型的变分质量控制方案,讨论了该方案的初始启动和关键参数,并检验分析了其适用性与有效性。同时,以全球预报系统(GFS)资料作为背景场,利用探空、地面、船舶、飞机、云迹风等常规观测资料和COSMIC卫星反演资料进行同化和预报,分析了华南地区特大暴雨的个例试验和2013年8月共31 d的批量试验。试验结果表明:变分质量控制能够依据观测资料的不同质量对观测权重进行合理调整,对位势高度、气压、风、比湿的分析增量场和分析场改善显著,尤其在强降水区具有更加明显的效果;对降水落区、降水强度及中心位置的预报质量具有较好的提高,特别对暴雨、大暴雨等较大降水量级的预报能力反映出更好的改善效果,充分显示了变分质量控制在中小尺度剧烈天气过程中对同化分析和预报的重要作用。     

Argo大洋观测资料的同化及其在短期气候预测和海洋分析中的应用

国际Argo(Array for Real-time Geostrophic Oceanography)计划的实施,提供了前所未有的全球深海大洋0～2000 m水深范围内的海水温度和盐度观测资料,在大气和海洋科研业务中应用这一全新的资料,是深入认识大气和海洋变异、提高我国气候预测、海洋监测分析和预报能力的一个关键所在.通过开发非线性温—盐协调同化方案和利用同化高度计资料来调整模式的温度和盐度场,建立了可同化包括Argo等多种海洋观测资料的全球海洋资料变分同化系统,提高了对全球海洋的监测分析能力.实现了海洋资料同化系统与全球海气耦合模式的耦合,显著提高了短期气候预测水平.利用Argo资料改进了海洋动力模式中的物理过程参数化方案,有效提高了海洋模式对真实大洋的模拟能力和对厄尔尼诺/拉尼娜的预测能力.开发了利用Argo浮标漂流轨迹推算全球海洋表层和中层流的方法,提高了推算的全球表层流、中层流资料质量,有效弥补了洋流观测的匮乏.     

A New Global Four-Dimensional Variational Ocean Data Assimilation System and Its Application

A four-dimensional variational data assimilation （4DVar） system of the LASG/IAP Climate Ocean Model, version 1.0 （LICOM1.0）, named LICOM-3DVM, has been developed using the three-dimensional variational data assimilation of mapped observation （3DVM）, a 4DVar method newly proposed in the past two years. Two experiments with 12-year model integrations were designed to validate it. One is the assimilation run, called ASSM, which incorporated the analyzed weekly sea surface temperature （SST） fields from Reynolds and Smith （OISST） between 1990 and 2001 once a week by the LICOM-3DVM. The other is the control run without any assimilation, named CTL. ASSM shows that the simulated temperatures of the upper ocean （above 50 meters）, especially the SST of equatorial Pacific, coincide with the Tropic Atmosphere Ocean （TAO） mooring data, the World Ocean Atlas 2001 （WOA01） data and the Met Office Hadley Centre＇s sea ice and sea surface temperature （HadISST） data. It decreased the cold bias existing in CTL in the eastern Pacific and produced a Nifio index that agrees with observation well. The validation results suggest that the LICOM-3DVM is able to effectively adjust the model results of the ocean temperature, although it＇s hard to correct the subsurface results and it even makes them worse in some areas due to the incorporation of only surface data. Future development of the LICOM-3DVM is to include subsurface in situ observations and satellite observations to further improve model simulations.     

Numerical experiments for Typhoon Dan incorporating AMSU-A retrieved data with 3DVM

Two sets of assimilation experiments on a landfalling typhoon—Typhoon Dan(1999)over the western North Pacific were designed to compare the performances of two kinds of variational data assimilation schemes that are the 3-Dimensional Variational data assimilation of Mapped observation(3DVM)and the 4-dimensional variational data assimilation(4DVar).Results show that:(1)both the 3DVM and 4DVar successfully improved the simulations of typhoon intensity and track incorporating the satellite AMSU-A retrieved temperature and wind data into the initial conditions,and the 3DVM more significantly due to the flow-dependent of background error covariance matrix and observation error covariance matrix like 3-dimensional variational data assimilation(3DVar)circle;(2)inclusions of extra model integration iterations at each observation time in the 3DVM make it more consistent with prediction model;(3)the 3DVM is much more time-saving due to the exclusion of the adjoint technique in it.     

A New Approach to Data Assimilation

A significant attempt to design a timesaving and efficient four-dimensional variational data assimilation (4DVar) has been made in this paper, and a new approach to data assimilation, which is noted as 'three-dimensional variational data assimilation of mapped observation (3DVM)' is proposed, based on the new concept of mapped observation and the new idea of backward 4DVar. Like the available 4DVar, 3DVM produces an optimal initial condition (IC) that is consistent with the prediction model due to the inclusion of model constraints and best fits the observations in the assimilation window through the model solution trajectory. Different from the 4DVar, the IC derived from 3DVM is located at the end of the assimilation window rather than at the beginning conventionally. This change greatly reduces the computing cost for the new approach, which is almost the same as that of the three-dimensional variational data assimilation (3DVar). Especially, such a change is able to improve assimilation accuracy because it does not need the tangential linear and adjoint approximations to calculate the gradient of cost function. Therefore, in numerical test, the new approach produces better IC than 4DVar does for 72-h simulation of TY9914 (Dan), by assimilating the three-dimensional fields of temperature and wind retrieved from the Advanced Microwave Sounding Unit-A (AMSU-A) observations. Meanwhile, it takes only 1/7 of the computing cost that the 4DVar requires for the same initialization with the same retrieved data.     

Initialization and simulation of a landfalling typhoon using a variational bogus mapped data assimilation (BMDA)

Summary Recently, a new data assimilation method called “3-dimensional variational data assimilation of mapped observation (3DVM)” has been developed by the authors. We have shown that the new method is very efficient and inexpensive compared with its counterpart 4-dimensional variational data assimilation (4DVar). The new method has been implemented into the Penn State/NCAR mesoscale model MM5V1 (MM5_3DVM). In this study, we apply the new method to the bogus data assimilation (BDA) available in the original MM5 with the 4DVar. By the new approach, a specified sea-level pressure (SLP) field (bogus data) is incorporated into MM5 through the 3DVM (for convenient, we call it variational bogus mapped data assimilation – BMDA) instead of the original 4DVar data assimilation. To demonstrate the effectiveness of the new 3DVM method, initialization and simulation of a landfalling typhoon – typhoon Dan (1999) over the western North Pacific with the new method are compared with that with its counterpart 4DVar in MM5. Results show that the initial structure and the simulated intensity and track are improved more significantly using 3DVM than 4DVar. Sensitivity experiments also show that the simulated typhoon track and intensity are more sensitive to the size of the assimilation window in the 4DVar than that in the 3DVM. Meanwhile, 3DVM takes much less computing cost than its counterpart 4DVar for a given time window.     

A global ocean temperature and altimeter data assimilation system for studies of climate variability

 An ocean data assimilation (ODA) system which can assimilate both temperature and altimeter observations has been applied to the global ocean and tested between January 1993–October 1996. A statistical method has been used to convert sea surface height (SSH) anomalies observations from TOPEX/POSEIDON into synthetic temperature profiles. The innovative aspect of this method is the introduction of time dependency in the correlations used to transform the altimeter observations into temperature corrections. The assimilation system is based on a univariate variational optimal interpolation scheme applied to assimilate both in situ and synthetic temperature profiles. In addition, a longer global analysis for the upper-ocean temperature starting from January 1979 and ending November 1997, has been produced to examine the skill of sea temperature assimilation with a rather simple and practical method. The temperature analysis shows encouraging improvement over a corresponding ocean simulation when compared to independent (not assimilated) temperature data both at seasonal and interannual time scales. However, the univariate data assimilation of hydrographic data does not result in an improvement of the velocity field. In fact the assimilation of sparse in situ data can introduce unrealistic spatial variability in the temperature field which affects the velocity field in a negative way. This deficiency is partially overcome when we also assimilate altimeter observations since the coverage is complete and uniform for this data. In particular, our study shows that temperature corrections due to the altimeter signal have a positive impact on the current system in the tropical Pacific. Received: 28 May 2000 / Accepted: 6 November 2000     

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

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

The impact of AIRS atmospheric temperature and moisture profiles on hurricane forecasts: Ike (2008) and Irene (2011)

Atmospheric Infra Red Sounder(AIRS) measurements are a valuable supplement to current observational data, especially over the oceans where conventional data are sparse. In this study, two types of AIRS-retrieved temperature and moisture profiles, the AIRS Science Team product(Sci Sup) and the single field-of-view(SFOV) research product, were evaluated with European Centre for Medium-Range Weather Forecasts(ECMWF) analysis data over the Atlantic Ocean during Hurricane Ike(2008) and Hurricane Irene(2011). The evaluation results showed that both types of AIRS profiles agreed well with the ECMWF analysis, especially between 200 h Pa and 700 h Pa. The average standard deviation of both temperature profiles was approximately 1 K under 200 h Pa, where the mean AIRS temperature profile from the AIRS Sci Sup retrievals was slightly colder than that from the AIRS SFOV retrievals. The mean Sci Sup moisture profile was slightly drier than that from the SFOV in the mid troposphere. A series of data assimilation and forecast experiments was then conducted with the Advanced Research version of the Weather Research and Forecasting(WRF) model and its three-dimensional variational(3DVAR) data assimilation system for hurricanes Ike and Irene. The results showed an improvement in the hurricane track due to the assimilation of AIRS clear-sky temperature profiles in the hurricane environment. In terms of total precipitable water and rainfall forecasts, the hurricane moisture environment was found to be affected by the AIRS sounding assimilation.Meanwhile, improving hurricane intensity forecasts through assimilating AIRS profiles remains a challenge for further study.     

台风数值模拟中bogus涡旋四维变分同化的研究

研究用MM5四维变分同化系统同化bogus涡旋变量对台风数值模拟的影响,对2002年13号台风Phanfone进行48h预报,结果表明：用四维变分方法同化bogus涡旋变量效果优于在背景场中直接植入模型涡旋的传统bogus技术。同时同化bogus涡旋的风、气压、温度和比湿效果不如只同化气压和其中任一变量。同化bogus涡旋的气压和比湿,能够调整模式初始场的其他变量,得到结构更清晰、与环境场更协调的初始涡旋,从而在一定程度上提高了台风Phanfone路径和强度预报的准确性。     

Effect of 3DVAR assimilation of MODIS temperature and humidity profiles on the dynamic and thermodynamic features of three monsoon depressions over the Bay of Bengal

The three-dimensional variational data assimilation (3DVAR) technique in the advanced weather research and forecast model is used to study the impact of assimilating Moderate Resolution Spectroradiometer (MODIS) retrieved temperature and humidity profiles on the dynamic and thermodynamic features for three monsoon depressions over the Bay of Bengal, India. For better understanding of the role of various physical processes in the evolution of monsoon depression, a detailed diagnostic study is performed on all the three depression cases. Numerical experiments were conducted in a system of two-way nested domains with a horizontal resolution of 36 and 12 km, respectively. The assimilation of MODIS data did improve the mean sea level pressure patterns and spatial distribution of rainfall patterns in all the three monsoon depression cases studied. Higher values of equitable threat score and lower bias values are seen consistently for the entire rainfall threshold range and for all the three depression cases with 3DVAR assimilation of MODIS temperature and humidity profiles. The current operational regional models in India do not ingest the MODIS temperature and humidity profiles and hence the present study is particularly relevant to the operational forecasting community in India in their ongoing efforts to improve weather forecasting over India.