Data Assimilation with the Ensemble Kalman Filter and the SEIK Filter applied to a Finite Element Model of the North Atlantic

Abstract Currently there are different approaches to filter algorithms based on the Kalman filter. One of the most used filter algorithms is the Ensemble Kalman Filter (EnKF). It uses a Monte Carlo approach to the filtering problem. Another approach is given by the Singular Evolutive Extended Kalman (SEEK) and Singular Evolutive Interpolated Kalman (SEIK) filters. These filters operate explicitly on a low-dimensional error space which is represented by an ensemble of model states. The EnKF and the SEIK filter have been implemented within a parallel data assimilation framework in the Finite Element Ocean Model FEOM. In order to compare the filter performances of the algorithms, several data assimilation experiments are performed. The filter algorithms have been applied with a model configuration of FEOM for the North Atlantic to assimilate the sea surface height in twin experiments. The dependence of the filter estimates on the represented error subspace is discussed. In the experiments the SEIK algorithm provides better estimates than the EnKF. Furthermore, the SEIK filter is much cheaper in terms of computing time.

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

[2]  Jens Schröter,et al.  A comparison of error subspace Kalman filters , 2005 .

[3]  P. Houtekamer,et al.  A Sequential Ensemble Kalman Filter for Atmospheric Data Assimilation , 2001 .

[4]  Jens Schröter,et al.  A finite-element ocean model: principles and evaluation , 2004 .

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

[6]  D. Pham.  Stochastic Methods for Sequential Data Assimilation in Strongly Nonlinear Systems , 2001 .

[7]  Jens Schröter,et al.  Evaluation of an eddy-permitting finite-element ocean model in the North Atlantic , 2004 .

[8]  Jeffrey L. Anderson.  An Ensemble Adjustment Kalman Filter for Data Assimilation , 2001 .

[9]  Pierre Brasseur,et al.  Assimilation of altimetric data in the mid-latitude oceans using the Singular Evolutive Extended Kalman filter with an eddy-resolving, primitive equation model , 1999 .

[10]  Dinh-Tuan Pham,et al.  A simplified reduced order Kalman filtering and application to altimetric data assimilation in Tropical Pacific , 2002 .

[11]  G. Evensen.  Sampling strategies and square root analysis schemes for the EnKF , 2004 .

[12]  Lars Nerger,et al.  Parallel filter algorithms for data assimilation in oceanography , 2003 .

[13]  K. Brusdala,et al.  A demonstration of ensemble-based assimilation methods with a layered OGCM from the perspective of operational ocean forecasting systems , 2003 .

[14]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[15]  Lars Nerger,et al.  PDAF - The Parallel Data Assimilation Framework: Experiences with Kalman Filtering , 2005 .

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