Efficient Parameterization of the Observation Error Covariance Matrix for Square Root or Ensemble Kalman Filters: Application to Ocean Altimetry

Abstract In the Kalman filter standard algorithm, the computational complexity of the observational update is proportional to the cube of the number y of observations (leading behavior for large y). In realistic atmospheric or oceanic applications, involving an increasing quantity of available observations, this often leads to a prohibitive cost and to the necessity of simplifying the problem by aggregating or dropping observations. If the filter error covariance matrices are in square root form, as in square root or ensemble Kalman filters, the standard algorithm can be transformed to be linear in y, providing that the observation error covariance matrix is diagonal. This is a significant drawback of this transformed algorithm and often leads to an assumption of uncorrelated observation errors for the sake of numerical efficiency. In this paper, it is shown that the linearity of the transformed algorithm in y can be preserved for other forms of the observation error covariance matrix. In particular, quit...

[1]  F. Rabier,et al.  The potential of high‐density observations for numerical weather prediction: A study with simulated observations , 2003 .

[2]  P. Houtekamer,et al.  Data Assimilation Using an Ensemble Kalman Filter Technique , 1998 .

[3]  Jeffrey L. Anderson A Local Least Squares Framework for Ensemble Filtering , 2003 .

[4]  J. Whitaker,et al.  Ensemble Square Root Filters , 2003, Statistical Methods for Climate Scientists.

[5]  Jean-Michel Brankart,et al.  Implementation of a multivariate data assimilation scheme for isopycnic coordinate ocean models: Application to a 1993–1996 hindcast of the North Atlantic Ocean circulation , 2003 .

[6]  G. Wahba,et al.  A Correspondence Between Bayesian Estimation on Stochastic Processes and Smoothing by Splines , 1970 .

[7]  P. Woodworth,et al.  Editorial (From the issue entitled "Special Issue in honour and in memory of Christian Le Provost - Modelling, observing, and forecasting sea level, ocean tides and ocean circulation: reviews and recent progress") , 2006 .

[8]  Thierry Penduff,et al.  Influence of numerical schemes on current-topography interactions in 1/4° global ocean simulations , 2007 .

[9]  P. McIntosh Oceanographic data interpolation: Objective analysis and splines , 1990 .

[10]  W. Brown,et al.  On the origins of the North Brazil Current , 1994 .

[11]  Florence Rabier,et al.  The interaction between model resolution, observation resolution and observation density in data assimilation: A one‐dimensional study , 2002 .

[12]  Ichiro Fukumori,et al.  A Partitioned Kalman Filter and Smoother , 2002 .

[13]  Istvan Szunyogh,et al.  A Local Ensemble Kalman Filter for Atmospheric Data Assimilation , 2002 .

[14]  Jean-Michel Brankart,et al.  Optimal Analysis of In Situ Data in the Western Mediterranean Using Statistics and Cross-Validation , 1996 .

[15]  David M. Fratantoni,et al.  Rings of the North Brazil Current: Their structure and behavior inferred from observations and a numerical simulation , 1995 .

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

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

[18]  Jean-Michel Brankart,et al.  Assimilation of sea-surface temperature and altimetric observations during 1992–1993 into an eddy permitting primitive equation model of the North Atlantic Ocean , 2003 .