Approximate Data Assimilation Schemes for Stable and Unstable Dynamics

Two suboptimal data assimilation schemes for stable and unstable dynamics are introduced. The first scheme, the partial singular value decomposition filter, is based on the most dominant singular modes of the tangent linear propagator. The second scheme, the partial eigendecomposition filter, is based on the most dominant eigenmodes of the propagated analysis error covariance matrix. Both schemes rely on iterative procedures like the Lanczos algorithm to compute the relevant modes. The performance of these schemes is evaluated for a shallow-water model linearized about an unstable Bickley jet. The results are contrasted against those of a reduced resolution filter, in which the gains used to update the state vector are calculated from a lower-dimensional dynamics than the dynamics that evolve the state itself. The results are also contrasted against the exact results given by the Kalman filter. These schemes are validated for the case of stable dynamics as well. The two new approximate assimilation schemes are shown to perform well with relatively few modes computed. Adaptive tuning of a modeled trailing error covariance for all three of these low-rank approximate schemes enhances performance and compensates for the approximation employed.

[1]  M. Yanai,et al.  Finite Difference Approximation for the Barotropic Instability Problem , 1968 .

[2]  J. Pedlosky.  Geophysical Fluid Dynamics , 1979 .

[3]  G. J. Haltiner.  Numerical Prediction and Dynamic Meteorology , 1980 .

[4]  Brian F. Farrell,et al.  Optimal Excitation of Baroclinic Waves , 1989 .

[5]  A. Girard.  A fast ‘Monte-Carlo cross-validation’ procedure for large least squares problems with noisy data , 1989 .

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

[7]  Anne E. Trefethen,et al.  Hydrodynamic Stability Without Eigenvalues , 1993, Science.

[8]  Eugenia Kalnay,et al.  Ensemble Forecasting at NMC: The Generation of Perturbations , 1993 .

[9]  P. Brasseur,et al.  Data Assimilation: Tools for Modelling the Ocean in a Global Change Perspective , 1994 .

[10]  Michael Ghil,et al.  Tracking Atmospheric Instabilities with the Kalman Filter. Part 1: Methodology and One-Layer Resultst , 1994 .

[11]  Andrew M. Moore,et al.  Using Adjoint Models for Stability and Predictability Analysis , 1994 .

[12]  Roberto Buizza,et al.  Sensitivity of optimal unstable structures , 1994 .

[13]  S. CohnData.  Assessing the Eeects of Data Selection with Dao's Physical-space Statistical Analysis System , 1994 .

[14]  R. Todling,et al.  Suboptimal Schemes for Atmospheric Data Assimilation Based on the Kalman Filter , 1994 .

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

[16]  D. Dee.  On-line Estimation of Error Covariance Parameters for Atmospheric Data Assimilation , 1995 .

[17]  P. Courtier.  Introduction to Numerical Weather Prediction Data Assimilation , 1996 .

[18]  F. Molteni,et al.  The ECMWF Ensemble Prediction System: Methodology and validation , 1996 .