Abstract
When applied to strongly nonlinear chaotic dynamics the extended Kalman filter (EKF) is prone to divergence due to the difficulty of correctly forecasting the forecast error probability density function. In operational forecasting applications ensemble Kalman filters circumvent this problem with empirical procedures such as covariance inflation. This paper presents an extension of the EKF that includes nonlinear terms in the evolution of the forecast error estimate. This is achieved starting from a particular square-root implementation of the EKF with assimilation confined in the unstable subspace (EKF-AUS), that is, the span of the Lyapunov vectors with non-negative exponents. When the error evolution is nonlinear, the space where it is confined is no more restricted to the unstable and neutral subspace causing filter divergence. The algorithm presented here, denominated EKF-AUS-NL, includes the nonlinear terms in the error dynamics: These result from the nonlinear interaction among the leading Lyapunov vectors and account for all directions where the error growth may take place. Numerical results show that with the nonlinear terms included, filter divergence can be avoided. We test the algorithm on the Lorenz96 model, showing very promising results.
- Received 4 December 2014
DOI:https://doi.org/10.1103/PhysRevE.91.042905
©2015 American Physical Society