The representer method, the ensemble Kalman filter and the ensemble Kalman smoother: A comparison study using a nonlinear reduced gravity ocean model

This paper compares contending advanced data assimilation algorithms using the same dynamical model and measurements. Assimilation experiments use the ensemble Kalman filter (EnKF), the ensemble Kalman smoother (EnKS) and the representer method involving a nonlinear model and synthetic measurements of a mesoscale eddy. Twin model experiments provide the “truth” and assimilated state. The difference between truth and assimilation state is a mispositioning of an eddy in the initial state affected by a temporal shift. The systems are constructed to represent the dynamics, error covariances and data density as similarly as possible, though because of the differing assumptions in the system derivations subtle differences do occur. The results reflect some of these differences in the tangent linear assumption made in the representer adjoint and the temporal covariance of the EnKF, which does not correct initial condition errors. These differences are assessed through the accuracy of each method as a function of measurement density. Results indicate that these methods are comparably accurate for sufficiently dense measurement networks; and each is able to correct the position of a purposefully misplaced mesoscale eddy. As measurement density is decreased, the EnKS and the representer method retain accuracy longer than the EnKF. While the representer method is more accurate than the sequential methods within the time period covered by the observations (particularly during the first part of the assimilation time), the representer method is less accurate during later times and during the forecast time period for sparse networks as the tangent linear assumption becomes less accurate. Furthermore, the representer method proves to be significantly more costly (2–4 times) than the EnKS and EnKF even with only a few outer iterations of the iterated indirect representer method.     

A singular evolutive extended Kalman filter to assimilate ocean color data in a coupled physical–biochemical model of the North Atlantic ocean

Within the European DIADEM project, a data assimilation system for coupled ocean circulation and marine ecosystem models has been implemented for the North Atlantic and the Nordic Seas. One objective of this project is to demonstrate the relevance of sophisticated methods to assimilate satellite data such as altimetry, surface temperature and ocean color, into realistic ocean models. In this paper, the singular evolutive extended Kalman (SEEK) filter, which is an advanced assimilation scheme where three-dimensional, multivariate error statistics are taken into account, is used to assimilate ocean color data into the biological component of the coupled system. The marine ecosystem model, derived from the FDM model [J. Mar. Res. 48 (1990) 591], includes 11 nitrogen and carbon compartments and describes the synthesis of organic matter in the euphotic zone, its consumption by animals of upper trophic levels, and the recycling of detritic material in the deep ocean. The circulation model coupled to the ecosystem is the Miami isopycnic coordinate ocean model (MICOM), which covers the Atlantic and the Arctic Oceans with an enhanced resolution in the North Atlantic basin. The model is forced with realistic ECMWF ocean/atmosphere fluxes, which permits to resolve the seasonal variability of the circulation and mixed layer properties. In the twin assimimation experiments reported here, the predictions of the coupled model are corrected every 10 days using pseudo-measurements of surface phytoplankton as a substitute to chlorophyll concentrations measured from space. The diagnostics of these experiments indicate that the assimilation is feasible with a reduced-order Kalman filter of small rank (of order 10) as long as a sufficiently good identification of the error structure is available. In addition, the control of non-observed quantities such as zooplankton and nitrate concentrations is made possible, owing to the multivariate nature of the analysis scheme. However, a too severe truncation of the error sub-space downgrades the propagation of surface information below the mixed layer. The reduction of the actual state vector to the surface layers is therefore investigated to improve the estimation process in the perspective of sea-viewing wide field-of-view sensor (SeaWiFS) data assimilation experiments.     

Weak and strong constraint data assimilation in the inverse Regional Ocean Modeling System (ROMS): Development and application for a baroclinic coastal upwelling system

We describe the development and preliminary application of the inverse Regional Ocean Modeling System (ROMS), a four dimensional variational (4DVAR) data assimilation system for high-resolution basin-wide and coastal oceanic flows. Inverse ROMS makes use of the recently developed perturbation tangent linear (TL), representer tangent linear (RP) and adjoint (AD) models to implement an indirect representer-based generalized inverse modeling system. This modeling framework is modular. The TL, RP and AD models are used as stand-alone sub-models within the Inverse Ocean Modeling (IOM) system described in [Chua, B.S., Bennett, A.F., 2001. An inverse ocean modeling system. Ocean Modell. 3, 137–165.]. The system allows the assimilation of a wide range of observation types and uses an iterative algorithm to solve nonlinear assimilation problems. The assimilation is performed either under the perfect model assumption (strong constraint) or by also allowing for errors in the model dynamics (weak constraints). For the weak constraint case the TL and RP models are modified to include additional forcing terms on the right hand side of the model equations. These terms are needed to account for errors in the model dynamics.Inverse ROMS is tested in a realistic 3D baroclinic upwelling system with complex bottom topography, characterized by strong mesoscale eddy variability. We assimilate synthetic data for upper ocean (0–450 m) temperatures and currents over a period of 10 days using both a high resolution and a spatially and temporally aliased sampling array. During the assimilation period the flow field undergoes substantial changes from the initial state. This allows the inverse solution to extract the dynamically active information from the synthetic observations and improve the trajectory of the model state beyond the assimilation window. Both the strong and weak constraint assimilation experiments show forecast skill greater than persistence and climatology during the 10–20 days after the last observation is assimilated.Further investigation in the functional form of the model error covariance and in the use of the representer tangent linear model may lead to improvement in the forecast skill.     

珠江河口对外海潮波响应的固有频率研究——采用系统控制论方法

采用受控的自回归滑动平均模型(CARMA)分别研究了珠江河口洪、枯季对外海潮波响应的固有频率。枯季采用单输入单输出模型,洪季采用双输入单输出模型,通过验证、检验表明,模型(CARMA)可用于寻求复杂网河对外海潮波响应的固有频率,方法可靠,易操作。通过频率增益响应分析,发现枯季珠江河口的固有频率以浅水分潮频段为主,对应周期介于6~10.0 h之间,枯季固有频率未出现在全日潮波对应的频段。洪季珠江河口的固有频率也以浅水分潮频段为主,珠江河口对浅水分潮的系统增益响应幅值最大。固有频率接近浅水分潮频段的站点,枯季较洪季固有频率减小,对应周期延长。