Effectiveness of data assimilation using probabilistic techniques in a model for synoptic velocity field evolution

The results of numerical experiments involving data assimilation by linear and non-linear models for synoptic ocean dynamics are presented. This paper considers the model's response to data assimilation at diverse values of the relative error in the determination of the initial fields and space-time discreteness of data assimilation. The data obtained using the optimal filtering are compared with those provided by the optimal interpolation (OI). A qualitative difference in the response of the linear and non-linear models to data assimilation is noted. The limiting values of the time-space discreteness of observations assimilation are determined, which allow correct application of OI.Translated by Vladimir A. Puchkin.     

Assessment of the effectiveness of POLYMODE data assimilation by a numerical barotropic model for synoptic motions

This paper considers the results of numerical experiments involving POLYMODE data assimilation by a barotropic model for synoptic ocean dynamics. The model's response the data assimilation for various space-time discretenesses of assimilation is studied. Results derived from the application of optimal interpolation algorithms and modified optimal filtration algorithms are compared. Qualitative similarity to the calculations carried out through the simulation modelling technique is noted. Optimal assimilation algorithms are determined, depending on the space-time discreteness used. An optimal sampling discreteness for the POLYMODE conditions is suggested.Translated by Vladimir A. Puchkin.     

Assimilation of altimetry data in the numerical model of the circulation in the Black Sea: The model numerical experiments

The numerical algorithm of the Kalman optimum filtration generalized for the case of three-dimensional fields and a multicomponent vector of the ocean state, with level measurements discrete over space and time being available, is given. The results of model numerical experiments on the assimilation of data on the Black Sea level are given. An attempt to estimate the effect of the time interval of data input on the results of field reconstruction was made.Translated by Mikhail M. Trufanov.     

基于不同背景场条件的雷达资料同化在登陆台风“桑美”中的应用研究

本文针对2006年登陆我国的超强台风“桑美”,分别采用美国国家环境预报中心的全球预报系统(Global Forecasting System, GFS)再分析资料和日本气象厅(Japan Meteorological Agency, JMA)区域客观再分析资料作为背景场,利用中尺度数值模式WRF(Weather Research and Forecasting Model)及其三维变分同化系统进行多普勒雷达资料同化和数值模拟试验,考察不同的背景场条件下雷达资料同化对台风初始场、内部结构及其随后确定性预报的影响。结果表明:GFS试验和JMA试验在同化了雷达资料之后分析出的台风700 hPa风场和500 hPa高度场相比其初始场均有所增强,JMA试验在3 h同化窗内的均方根误差和最小海平面气压的改进效果均比GFS试验显著,同时对台风动力和热力结构的改进效果也优于GFS试验;JMA试验对台风降水、路径、强度的预报均优于GFS试验,且能预报出台风前沿的降水,更加接近观测实况。     

海洋资料同化对气候季节-年际预测技巧及初始场的影响试验

集合卡尔曼滤波(Ensemble Kalman filter, EnKF)是一种国内外广泛使用的海洋资料同化方案, 用集合成员的状态集合表征模式的背景误差协方差, 结合观测误差协方差, 计算卡尔曼增益矩阵, 有效地将观测信息添加到模式初始场中。由于季节、年际预测很大程度上受到初始场的影响, 因此资料同化可以提高模式的预测性能。本文在NUIST-CFS1.0预测系统逐日SST nudging的初始化方案上, 利用EnKF在每个月末将全场(full field)海表温度(sea surface temperature, SST)、温盐廓线(in-situ temperature and salinity profiles, T-S profiles)以及卫星观测海平面高度异常(sea level anomalies, SLA)观测资料同化到模式初始场中, 对比分析了无海洋资料同化以及加入同化后初始场的区别、加入海洋资料同化后模式提前1~24个月预测性能的差异以及对于厄尔尼诺-南方涛动(El Niño-southern oscillation, ENSO)预测技巧的影响。结果表明, 加入海洋资料同化能有效地改进初始场, 并且呈现随深度增加初始场改进越显著的特征。加入同化后, 对全球SST、次表层海水温度的平均预测技巧均有一定的提高, 也表现出随深度增加预测技巧改进越明显的特征。但加入海洋资料同化后, 模式对ENSO的预测技巧有所下降, 可能是由于模式误差的存在, 使得同化后的预测初始场从接近观测的状态又逐渐恢复到与模式动力相匹配的状态, 加剧了赤道太平洋冷舌偏西、中东部偏暖的气候平均态漂移。     

MASNUM-WAM海浪模式集合Kalman滤波同化研究-I.风场扰动对海浪模拟影响

模式集合样本的代表性和观测信息的可靠性是制约数据同化效果的重要因素,而前者对海浪模式同化的影响尤为显著。由于海浪模式对初始场的敏感性较弱,来自大气的风输入源函数是海浪的重要能量输入,如何合理地对风输入进行扰动,构造海浪的集合模式运行,是实现和改进海浪模式集合Kalman滤波同化的关键问题。为了实现海浪模式集合运行,本文提出了风场的三种集合扰动方案,分别为：纯随机数、随机场和时间滞后的风场扰动方法。本研究利用2014年1月ECMWF全球风场,基于这三种风场扰动方法开展了集合海浪模式的集合运行实验,并统计分析了海浪特征要素（有效波高）和二维波数谱对风场扰动的响应。结果表明,随机场集合扰动方案所构造的风场集合效果最佳,所得海浪模拟结果的集合样本发散度适中,能够较为合理地反映背景误差的统计特征,可用于进一步的集合Kalman滤波海浪数据同化实验。     

TRMM海表降水率同化在热带气旋Danas(2001) 数值模拟中的应用研究

研究了TRMM/TMI海表降水率资料的四维变分同化在热带气旋(TC)数值模拟中的作用.使用中尺度气象模式MM5设计了若干数值试验模拟了TC Danas(2001)由热带低压初生到台风生成的发展过程.在满足MM5模式动力约束的前提下,将TRMM海表降水率资料直接同化进入较高分辨率(18 km)的模式初始场.结果表明,使用MM5模式的4D-VAR同化系统直接同化TRMM/TMI海表降水率资料是可行的.这种做法提高了TRMM资料的利用率,不仅在模式初始场中加入了更多实测信息,而且避免了两次同化(1DVAR+4DVAR)可能引起的误差.直接同化TRMM资料通过调整气压、温度、湿度等要素初始场,改善了模式对热带气旋结构(如暖心、涡度、散度)的描述和降水的模拟.在此基础上,同化不仅改进了对Danas强度的模拟,而且成功地模拟了热带气旋环境场的演变过程,因而改进了路径的模拟.     

Calculation of adapted black sea fields on the basis of assimilation of climatic temperature and salinity data into the model

A scheme of a periodic assimilation of climatic temperature and salinity data into a model is proposed. In accordance with the criteria chosen on the basis of numerical experiments, an optimum assimilation period of three hours is set. The proposed method allows a hydrodynamic adaptation of the fields of temperature, salinity, and flow velocity that are close to the climatic ones within the framework of the proposed criteria. The main feature of the adapted fields of vertical flows in the layer 0–700 m is their smooth character and an almost complete absence of small-scale noise. The extreme values of vertical velocity in this layer decrease by almost an order of magnitude as compared to analogous values calculated for the variant of data assimilation with a weak adjustment of fields.     

Optimization of the data assimilation algorithm in a model for synoptic dynamics of the velocity field

The results of numerical experiments involving data assimilation by a non-linear model for synoptic dynamics of the ocean are given. The results of application of simplified optimal filtration algorithms (with the correlation error function being parametrized for both linear and non-linear variants) and of the optimal interpolation are intercompared. The use of a non-linear equation to predict average values and a linear one to predict dispersion in the filtration algorithm is suggested as most appropriate for the cases, when anisotropy of the correlation function is to be included. In other cases, it is considered worth-while to apply optimal interpolation.Translated by V. Puchkin.     

Correcting the errors in the initial conditions and wind stress in storm surge simulation using an adjoint optimal technique

An adjoint data assimilation methodology is applied to the Princeton Ocean Model and is evaluated by obtaining “optimal” initial conditions, sea surface forcing conditions, or both for coastal storm surge modelling. By prescribing different error sources and setting the corresponding control variables, we performed several sets of identical twin experiments by assimilating model-generated water levels. The experiment results show that, when the forecasting errors are caused by the initial (or surface boundary) conditions, adjusting initial (or surface boundary) conditions accordingly can significantly improve the storm surge simulation. However, when the forecasting errors are caused by surface boundary (or initial) conditions, data assimilation targeting improving the initial (or surface boundary) conditions is ineffective. When the forecasting errors are caused by both the initial and surface boundary conditions, adjusting both the initial and surface boundary conditions leads to the best results. In practice, we do not know whether the errors are caused by initial conditions or surface boundary conditions, therefore it is better to adjust both initial and surface boundary conditions in adjoint data assimilation.     

直接同化卫星辐射数据的黄海海雾数值试验研究

基于中尺度大气模式WRF及其3DVAR.模块,采用循环3DVAR数据同化方案,针对6次明显的黄海海雾过程,实施了一系列直接同化ATOVS卫星辐射数据数值试验.在试验中设计了不同化任何观测数据、仅同化GTS常规数据、仅直接同化辐射数据,同时同化二者,以及同化不同疏密程度辐射数据的对比研究方案.利用地面水平能见度与卫星云图对模拟的海雾雾区进行了评估,并比较了各种同化方案所形成初始场的差异.对试验结果的统计分析表明:同化试验较好地再现了影响海雾的天气系统,模拟雾区与实际观测较为吻合,并且初始温度场和湿度场对比不同化任何观测数据的试验有明显的改善;仅同化辐射数据的结果略优于仅同化常规数据的结果,疏化或者只同化海上辐射数据几乎不影响模拟的雾区,但却可以大幅节约计算资源;同时同化常规数据与辐射数据的结果为单独同化它们所得结果的综合体现,总体效果最好.     

An Operational Circulation Modeling System for the Gulf of Maine/Georges Bank Region—Part I: Basic Elements

The basic elements of a prototype operational data assimilation modeling system that can provide near-real-time information on the ocean water property and circulation environment in the Gulf of Maine (GOM)/Georges Bank (GB) region are described in this paper. This application of the Harvard Ocean Prediction System (HOPS, Harvard University, Cambridge, MA) model includes development of protocols for the following: 1) the production of model initial fields from an objective blending of climatological and feature model (FM) hydrographic data with fishing-boat-measured bottom temperature data, 2) the ldquowarm startrdquo of the model to produce reasonably realistic initial model fields, 3) converting real-time Fleet Numerical Meteorological and Oceanographic Center (FNMOC, Monterey, CA) model nowcast and forecast winds and/or National Data Buoy Center (NDBC, Stennis Space Center, MS) operational wind measurements to model wind stress forcing fields, and 4) the assimilation of satellite-derived sea surface temperature (SST). These protocols are shown herein to evolve the initial model fields, which were dominated by climatological data, toward more dynamically balanced, realistic fields. Thus, the model nowcasts, with the assimilation of one SST field, are well positioned to produce reasonably realistic ocean fields within a few model days (MDs).     

Response of a coupled ocean-ice-atmosphere model to data assimilation in the tropical zone of the Pacific Ocean

This paper deals with numerical experiments based on the coupled ECHAM-HOPE model. The results of experiments are analyzed. The initial fields for the calculations over time periods from one month to one year are constructed based on the results of the data assimilation of temperature profiles from TOGA-TAO moorings. The perturbations of the initial fields and the propagation of these perturbations with specific computational time intervals are analyzed on the basis of the results of these experiments. It is shown that the strongest impacts of the perturbations are localized in specific regions of the World Ocean corresponding to the energetically active zones of the Earth. The mechanism of the transition of these perturbations is also studied. Different statistical properties of the ensemble of experiments are presented.     

云迹风资料同化对东亚海域风场数值模拟的影响

为了揭示云迹风资料同化对于东亚海域内各层模拟风场的影响,采用由日本GMS卫星红外和水汽通道资料反演得到的1999年6月底至7月初的云迹风资料,首先对这些云迹风资料进行质量评估,评估结果表明云迹风u,v分量的均方根误差和偏差都在所容许的误差范围之内.然后针对10个个例进行5种不同的四维同化模拟试验,多种同化试验和无同化控制试验的检验对比表明:在经过初始12 h常规高空地面资料和云迹风资料同化过程以后,模拟结果质量有了比较明显的提升;各种同化试验与控制试验间的最大差异发生在12h,也即同化结束的时刻,此后随着模拟时间的增长,两者间的差异逐渐变小;但24~36h同化试验对控制试验的改进效果又开始提升,36h后效果再次下降.由于海面上缺少常规测站,因此无法实现海面部分的测站同化,而这必将影响到海面部分的模拟质量,而云迹风资料恰可以祢补这个不足.试验表明云迹风资料对u,v风场的改进则主要集中在对流层高层(300~200hPa),对中低层的影响不很明显,这与云迹风资料主要集中在高层是相一致的.另外,敏感性试验表明经过质量控制的云迹风同化模拟的性能相对不采用质量控制的同化模拟有了一定程度的提高.     

在Zebiak-Cane模式中观测MJO强迫对ENSO可预报性的影响, 第1部分:对最大预报误差的影响

With the observational wind data and the Zebiak-Cane model, the impact of Madden-Julian Oscillation (MJO) as external forcing on El Ni(n)o–Southern Oscillation (ENSO) predictability is studied. The obs...     

三维斜压台风模式 Ⅱ.预报试验

一种斜压多重移动套网格台网模式在国家海洋环境预报中心已被应用于海洋环境预报。本文第一部分已描述了模式方程组和数值方法。本文继续概述模式网格、变分调整初始化和预报试验结果。最外粗网格域固定,内部细网格域随台风中心轨迹移动。模式中,粗细网格变量采用双向耦合。平衡方程和方程,理想台风场和变分调整方案被用于台风模式初始化。一种简单而有效的资料同化方法,即用第6h台风报和弱约束变分原理调整初始场,被应用于提高预报结果。最后本文给出预报试验结果。预报误差统计显示本模式对台风路径预报具有相当能力,而且可以提供海面风和气压场较好的预报。本模式已经与海浪模式联结,得到满意的波高预报结果。     

Assimilation of TOPEX Altimeter Data in a Reduced-Gravity Model of the Northwestern Indian Ocean

This study deals with the assimilation of TOPEX altimeter-derived sea level variability in a reduced-gravity model of the northwestern Indian Ocean. The assimilation has been done using variational assimilation method with adjoint technique. A cost function representing the misfit between the model and the data is minimized with the model equations acting as constraints. The initial conditions of the model are used as control parameters and the best-fit initial conditions are determined as a result of minimization carried out using a variable-storage quasi-Newton method. Experiment has been done with 10 days' cycle-averaged data of TOPEX. Assimilation for 10 days and 20 days has been performed and it has been seen that 20 days' assimilation of satellite data provides better results.     

Numerical analysis of hydrophysical fields on the northwestern shelf of the Black Sea

A numerical experiment with assimilation of hydrological observational data from a survey in October 2007 on the northwestern shelf of the Black Sea was carried out using the hydrodynamic model with nonlinear equations of motion, equations of heat and salt advection, and data assimilation. The results of this calculation are compared with thermohydrodynamic fields obtained without taking into account temperature and salinity measurements. It is shown that allowance for the observation data leads to qualitative and quantitative differences in the structure of the hydrophysical fields. Mesoscale eddies and intense jet streams that agree with satellite observations were found in the field of currents and were investigated. These eddies are not resolved in low-resolution field experiments.     

黄海海雾短时临近预报中云水路径的EnKF同化研究

在海雾的短时临近预报中,初始场的水汽凝结状态扮演着重要角色。为了改进初始场的云水含量,本文提出直接同化雾体云水信息的思路。针对2011年5月一次大范围的黄海海雾,借助EnKF (Ensemble Kalman Filter)方法,尝试进行了极轨卫星反演云水路径数据的同化试验。结果表明:(1)通过利用EnKF将云水混合比增加到背景场和分析场的控制变量中,构建云水观测数据与背景场之间的关系,实现云水路径数据的直接同化是可行的;(2)同化云水路径可显著改善海面气温与湿度状态,大幅提高海雾预报效果;(3)EnKF能够基于集合体动态统计流依赖的背景误差协方差是其取得良好同化效果的主要原因。值得指出的是,受集合样本误差的影响,需要特别关注云水含量与风之间的相关关系。