Assimilation of altimetric data in the mid-latitude oceans using the Singular Evolutive Extended Kalman filter with an eddy-resolving, primitive equation model

A new data assimilation scheme has been elaborated for ocean circulation models based on the concept of an evolutive, reduced-order Kalman filter. The dimension of the assimilation problem is reduced by expressing the initial error covariance matrix as a truncated series of orthogonal perturbations. This error sub-space evolves during the assimilation so as to capture the growing modes of the estimation error. The algorithm has been formulated in quite a general fashion to make it tractable with a large variety of ocean models and measurement types. In the present paper, we have examined three possible strategies to compute the evolution of the error subspace in the so-called Singular Evolutive Extended Kalman (SEEK) filter: the steady filter considers a time-independent error sub-space, the apprentice filter progressively enriches the error sub-space with the information learned from the innovation vector after each analysis step, and the dynamical filter updates the error sub-space according to the model dynamics. The SEEK filter has been implemented to assimilate synthetic observations of the surface topography in a non-linear, primitive equation model that uses density as vertical coordinate. A simplified box configuration has been adopted to simulate a Gulf Stream-like current and its associated eddies and gyres with a resolution of 20 km in the horizontal, and 4 levels in the vertical. The concept of twin experiments is used to demonstrate that the conventional SEEK filter must be complemented by a learning mechanism in order to model the misrepresented tail of the error covariance matrix. An approach based on the vertical physics of the isopycnal model, is shown particularly robust to control the velocity field in deep layers with surface observations only. The cost of the method makes it a suitable candidate for large-size assimilation problems and operational applications.

[1]  Keith Haines,et al.  Altimetric assimilation with water property conservation , 1996 .

[2]  Ichiro Fukumori,et al.  Fitting Dynamic Models to the Geosat Sea Level Observations in the Tropical Pacific Ocean. Part II: A Linear, Wind-driven Model , 1993 .

[3]  Arthur Gelb,et al.  Applied Optimal Estimation , 1974 .

[4]  Jacques Verron,et al.  A singular evolutive extended Kalman filter for data assimilation in oceanography , 1998 .

[5]  Jacques Verron,et al.  An extended Kalman filter to assimilate satellite altimeter data into a nonlinear numerical model of the tropical Pacific Ocean: Method and validation , 1999 .

[6]  C. Wunsch,et al.  How well does a 1/4° global circulation model simulate large-scale oceanic observations? , 1996 .

[7]  Geir Evensen,et al.  Comparison of sequential updating, Kalman filter and variational methods for assimilating Rossby waves in the simulated Geosat altimeter data into a quasi-geostrophic model , 1995 .

[8]  Ichiro Fukumori,et al.  Assimilation of TOPEX sea level measurements with a reduced‐gravity, shallow water model of the tropical Pacific Ocean , 1995 .

[9]  John D. McCalpin,et al.  Phenomenology of the low-frequency variability in a reduced-gravity, quasigeostrophic double-gyre model , 1996 .

[10]  Isabelle Blanchet,et al.  A Comparison of Adaptive Kalman Filters for a Tropical Pacific Ocean Model , 1997 .

[11]  G. Evensen,et al.  Assimilation of Geosat altimeter data for the Agulhas current using the ensemble Kalman filter with , 1996 .

[12]  Alexey Kaplan,et al.  Mapping tropical Pacific sea level : Data assimilation via a reduced state space Kalman filter , 1996 .

[13]  Robert N. Miller.  Toward the Application of the Kalman Filter to Regional Open Ocean Modeling , 1986 .

[14]  Robert M. Chervin,et al.  Ocean general circulation from a global eddy‐resolving model , 1992 .

[15]  R. Bleck.  Simulation of coastal upwelling frontogenesis with an isopycnic coordinate model , 1978 .

[16]  R. E. Young,et al.  A comparison of two methods for the assimilation of altimeter data into a shallow-water model , 1993 .

[17]  G. Evensen.  Using the Extended Kalman Filter with a Multilayer Quasi-Geostrophic Ocean Model , 1992 .

[18]  G. Evensen.  Sequential data assimilation with a nonlinear quasi‐geostrophic model using Monte Carlo methods to forecast error statistics , 1994 .

[19]  François Bouttier,et al.  The dynamics of error covariances in a barotropic model , 1993 .

[20]  Eric P. Chassignet,et al.  The Influence of Boundary Conditions on Midlatitude Jet Separation in Ocean Numerical Models , 1991 .

[21]  P. Courtier,et al.  Assimilation of Simulated Wind Lidar Data with a Kalman Filter , 1993 .

[22]  Rainer Bleck,et al.  A wind-driven isopycnic coordinate model of the north and equatorial Atlantic Ocean 2. The Atlantic Basin Experiments , 1990 .

[23]  Philippe Courtier,et al.  Unified Notation for Data Assimilation : Operational, Sequential and Variational , 1997 .

[24]  Andreas Oschlies,et al.  Assimilation of Geosat altimeter data into an eddy-resolving primitive equation model of the North Atlantic Ocean , 1996 .

[25]  P. Malanotte‐Rizzoli,et al.  An approximate Kaiman filter for ocean data assimilation: An example with an idealized Gulf Stream model , 1995 .

[26]  Rainer Bleck,et al.  Wind‐driven spin‐up in eddy‐resolving ocean models formulated in isopycnic and isobaric coordinates , 1986 .

[27]  Jean-Marc Molines,et al.  Assimilation of TOPEX/POSEIDON altimeter data into a circulation model of the North Atlantic , 1994 .

[28]  Allan R. Robinson,et al.  Assimilation of altimeter Eddy fields in a limited-area quasi-geostrophic model , 1987 .

[29]  D. Dee.  Simplification of the Kalman filter for meteorological data assimilation , 1991 .

[30]  Predictability experiments in the North Atlantic Ocean: Outcome of a quasi-geostrophic model with assimilation of TOPEX/POSEIDON altimeter data , 1996 .

[31]  M. Gavart,et al.  Isopycnal EOFs in the Azores Current Region: A Statistical Tool forDynamical Analysis and Data Assimilation , 1997 .

[32]  Eric Blayo,et al.  Complementarity of ERS 1 and TOPEX/POSEIDON altimeter data in estimating the ocean circulation: Assimilation into a model of the North Atlantic , 1997 .

[33]  A. Jazwinski.  Stochastic Processes and Filtering Theory , 1970 .

[34]  George L. Mellor,et al.  A Gulf Stream model and an altimetry assimilation scheme , 1991 .

[35]  Rainer Bleck,et al.  A Wind-Driven Isopycnic Coordinate Model of the North and Equatorial Atlantic Ocean 1 , 1990 .

[36]  G. Evensen,et al.  Analysis Scheme in the Ensemble Kalman Filter , 1998 .

[37]  Stephen E. Cohn,et al.  The Behavior of Forecast Error Covariances for a Kalman Filter in Two Dimensions , 1991 .

[38]  Anthony Rosati,et al.  The sea surface pressure formulation of rigid lid models. Implications for altimetric data assimilation studies , 1995 .