## Abstract It is well recognized that discontinuous analysis increments of sequential data assimilation systems, such as ensemble Kalman filters, might lead to spurious high‐frequency adjustment processes in the model dynamics. Various methods have been devised to spread out the analysis increment
A robust formulation of the ensemble Kalman filter
✍ Scribed by S. J. Thomas; J. P. Hacker; J. L. Anderson
- Publisher
- John Wiley and Sons
- Year
- 2009
- Tongue
- English
- Weight
- 338 KB
- Volume
- 135
- Category
- Article
- ISSN
- 0035-9009
- DOI
- 10.1002/qj.372
No coin nor oath required. For personal study only.
✦ Synopsis
Abstract
The ensemble Kalman filter (EnKF) can be interpreted in the more general context of linear regression theory. The recursive filter equations are equivalent to the normal equations for a weighted least‐squares estimate that minimizes a quadratic functional. Solving the normal equations is numerically unreliable and subject to large errors when the problem is ill‐conditioned. A numerically reliable and efficient algorithm is presented, based on the minimization of an alternative functional. The method relies on orthogonal rotations, is highly parallel and does not ‘square’ matrices in order to compute the analysis update. Computation of eigenvalue and singular‐value decompositions is not required. The algorithm is formulated to process observations serially or in batches and therefore easily handles spatially correlated observation errors. Numerical results are presented for existing algorithms with a hierarchy of models characterized by chaotic dynamics. Under a range of conditions, which may include model error and sampling error, the new algorithm achieves the same or lower mean square errors as the serial Potter and ensemble adjustment Kalman filter (EAKF) algorithms. Published in 2009 by John Wiley and Sons, Ltd.
📜 SIMILAR VOLUMES
## Abstract Ensemble Kalman filter techniques are widely used to assimilate observations into dynamical models. The phase‐space dimension is typically much larger than the number of ensemble members, which leads to inaccurate results in the computed covariance matrices. These inaccuracies can lead,
## Abstract The Met Office has been routinely running a short‐range ensemble prediction system since the summer of 2005. This system consists of two component ensembles, a global ensemble that provides lateral boundary conditions to a regional ensemble. The global ensemble calculates the initial co
## Abstract Ensemble Kalman Filter (EnKF) may have a longer spin‐up time to reach its asymptotic level of accuracy than the corresponding spin‐up time in variational methods (3D‐Var or 4D‐Var). During the spin‐up EnKF has to fulfill two independent requirements, namely that the ensemble mean be clo
Ensemble Kalman filter (EnKF) has been widely used as a sequential data assimilation method, primarily due to its ease of implementation resulting from replacing the covariance evolution in the traditional Kalman filter (KF) by an approximate Monte Carlo ensemble sampling. In this paper rigorous ana