Abstract

The ensemble Kalman filter (EnKF) is susceptible to losing track of observations, or ‘diverging’, when applied to large chaotic systems such as atmospheric and ocean models. Past studies have demonstrated the adverse impact of sampling error during the filter’s update step. We examine how system dynamics affect EnKF performance, and whether the absence of certain dynamic features in the ensemble may lead to divergence. The EnKF is applied to a simple chaotic model, and ensembles are checked against singular vectors of the tangent linear model, corresponding to short-term growth and Lyapunov vectors, corresponding to long-term growth. Results show that the ensemble strongly aligns itself with the subspace spanned by unstable Lyapunov vectors. Furthermore, the filter avoids divergence only if the full linearized long-term unstable subspace is spanned. However, short-term dynamics also become important as nonlinearity in the system increases. Non-linear movement prevents errors in the long-term stable subspace from decaying indefinitely. If these errors then undergo linear intermittent growth, a small ensemble may fail to properly represent all important modes, causing filter divergence. A combination of long and short-term growth dynamics are thus critical to EnKF performance. These findings can help in developing practical robust filters based on model dynamics.

Highlights

  • The ensemble Kalman filter (EnKF) is a flexible data assimilation tool for various geophysical applications

  • This work uses a simplified model of a chaotic attractor, the Lorenz-98 model (Lorenz and Emanuel, 1998), to clarify the dynamic features that affect EnKF performance

  • We find that the non-divergent filter ensemble only aligns weakly with singular vectors, while it aligns strongly with leading characteristic or backward Lyapunov vectors

Read more

Summary

Introduction

The ensemble Kalman filter (EnKF) is a flexible data assimilation tool for various geophysical applications. ‘inbreeding’, or the interaction of sampling error with other terms in the EnKF update These fixes often prove crucial for generating robust (non-divergent) estimates with the EnKF in various geophysical applications. This work examines an aspect influencing EnKF performance that has been typically overlooked in most filter divergence studies: the dynamic properties of the forecast model. We are interested in how EnKF errors, including those that arise from using small ensemble sizes, relate to dynamic features of the model, and how that correspondence impacts filter performance. We adopt a more diagnostic objective: we assess the impact of dynamic properties on EnKF performance for chaotic systems It is not readily apparent how the ensemble relates to different types of growth modes in the system, and how filter divergence depends on the ensemble filter’s ability to track these modes.

The ensemble Kalman filter
Chaos dynamics review
Singular vectors
Lyapunov vectors
Dynamical analysis and data assimilation
Experiment set-up
Results for 40-D Lorenz-98 model
Results for 6-D Lorenz-98 model
Implications for larger-scale models
Conclusion
Full Text
Paper version not known

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call