Using the Extended Kalman Filter with a Multilayer Quasi-Geostrophic Ocean Model

The formulation of the extended Kalman filter for a multilayer nonlinear quasi-geostrophic ocean circulation model is discussed. The nonlinearity in the ocean model leads to an approximative equation for error covariance propagation, where the transition matrix is dependent on the state trajectory. This nonlinearity complicates the dynamics of the error covariance propagation, and effects which are nonexistent in linear systems contribute significantly. The transition matrix can be split into two parts, where one part results in pure evolution of error covariances in the model velocity field, and the other part contains a statistical correction term caused by the nonlinearity in the model. This correction term leads to a linear unbounded instability, which is caused by the statistical linearization of the nonlinear error propagation equation. Different ways of handling this instability are discussed. Further, nonlinear small-scale instabilities also develop, since energy is accumulated at wavelengths 2Δx, owing to the numerical discretization. These small-scale oscillations are removed with a Shapiro filter, and the effect they have on the error covariance propagation is discussed. Some data assimilation experiments are performed using the full extended Kalman filter, to examine the properties of the filter. An experiment where only the first part of the transition matrix is used to propagate the error covariances is also performed. This simplified experiment actually performs better than the full extended Kalman filter because the unbounded instability associated with the statistical correction term is avoided.

[1]  J. Neumann,et al.  Numerical Integration of the Barotropic Vorticity Equation , 1950 .

[2]  E. Lorenz Deterministic nonperiodic flow , 1963 .

[3]  Ralph Shapiro,et al.  Smoothing, filtering, and boundary effects , 1970 .

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

[5]  R. Sweet A Cyclic Reduction Algorithm for Solving Block Tridiagonal Systems of Arbitrary Dimension , 1977 .

[6]  G. Bierman Factorization methods for discrete sequential estimation , 1977 .

[7]  William R. Holland,et al.  The Role of Mesoscale Eddies in the General Circulation of the Ocean—Numerical Experiments Using a Wind-Driven Quasi-Geostrophic Model , 1978 .

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

[9]  M. Ikeda Meanders and Detached Eddies of a Strong Eastward-Flowing Jet Using a Two-Layer Quasi-Geostrophic Model , 1981 .

[10]  S. Cohn,et al.  Applications of Estimation Theory to Numerical Weather Prediction , 1981 .

[11]  J. Apel,et al.  Mesoscale Eddies Detached from Spatially Growing Meanders in an Eastward-Flowing Oceanic jet Using a Two-Layer Quasi-Geostrophic Model , 1981 .

[12]  P. Smolarkiewicz The Multi-Dimensional Crowley Advection Scheme , 1982 .

[13]  Shenn-Yu Chao,et al.  Bimodality of the Kuroshio , 1984 .

[14]  J. M. Lewis,et al.  The use of adjoint equations to solve a variational adjustment problem with advective constraints , 1985 .

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

[16]  W. Budgell Nonlinear data assimilation for shallow water equations in branched channels , 1986 .

[17]  R. Stengel Stochastic Optimal Control: Theory and Application , 1986 .

[18]  W. Budgell Stochastic filtering of linear shallow water wave processes , 1987 .

[19]  W. Budgell,et al.  Ocean Data Assimilation and the Moan Filter: Spatial Regularity , 1987 .

[20]  William Carlisle Thacker,et al.  Fitting models to inadequate data by enforcing spatial and temporal smoothness , 1988 .

[21]  Andrew F. Bennett,et al.  Numerical simulation of flows with locally characteristic boundaries , 1988 .

[22]  W. Thacker,et al.  Fitting dynamics to data , 1988 .

[23]  O. Talagrand,et al.  Short-range evolution of small perturbations in a barotropic model , 1988 .

[24]  W. Budgell,et al.  The Kalman smoother for a linear quasi-geostrophic model of ocean circulation , 1989 .

[25]  Michael Ghil,et al.  Meteorological data assimilation for oceanographers. Part I: Description and theoretical framework☆ , 1989 .

[26]  William Carlisle Thacker,et al.  Data assimilation into a numerical equatorial ocean model. II. Assimilation experiments , 1989 .

[27]  Robert N. Miller,et al.  A Kalman Filter Analysis of Sea Level Height in the Tropical Pacific , 1989 .

[28]  Robert Bryan Long,et al.  Data assimilation into a numerical equatorial ocean model. I. The model and the assimilation algorithm , 1989 .

[29]  Johnny A. Johannessen,et al.  A Process Study of Mesoscale Meanders and Eddies in the Norwegian Coastal Current , 1989 .

[30]  P. Courtier,et al.  Variational assimilation of meteorological observations with the direct and adjoint shallow-water equations , 1990 .

[31]  Robert N. Miller Tropical Data Assimilation Experiments With Simulated Data' The Impact of the Tropical Ocean and Global Atmosphere Thermal Array for the Ocean , 1990 .

[32]  M. Ghil,et al.  Data assimilation in meteorology and oceanography , 1991 .

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

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

[35]  G. Evensen,et al.  Modeled and observed mesoscale circulation and wave‐current refraction during the 1988 Norwegian Continental Shelf Experiment , 1991 .

[36]  William H. Raymond,et al.  A review of recursive and implicit filters , 1991 .

[37]  Robert N. Miller,et al.  Ocean data assimilation using optimal interpolation with a quasi-geostrophic model , 1991 .