## On numerical properties of the ensemble Kalman filter for data assimilation.(English)Zbl 1195.93137

Summary: 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 analysis on the numerical errors of the EnKF is conducted in a general setting. Error bounds are provided and convergence of the EnKF to the exact Kalman filter is established. The analysis reveals that the ensemble errors induced by the Monte Carlo sampling can be dominant, compared to other errors such as the numerical integration error of the underlying model equations. Methods to reduce sampling errors are discussed. In particular, we present a deterministic sampling strategy based on cubature rules (qEnKF) which offers much improved accuracy. The analysis also suggests a less obvious fact – more frequent data assimilation may lead to larger numerical errors of the EnKF. Numerical examples are provided to verify the theoretical findings and to demonstrate the improved performance of the qEnKF.

### MSC:

 93E11 Filtering in stochastic control theory 93C15 Control/observation systems governed by ordinary differential equations 93E10 Estimation and detection in stochastic control theory 93C41 Control/observation systems with incomplete information

EnKF
Full Text:

### References:

 [1] Bishop, C. H.; Etherto, B. J.; Majumdar, S. J., Adaptive sampling with the ensemble transform Kalman filter, part i: theoretical aspects, Mon. Weather Rev., 129, 420-436 (2001) [2] Burgers, G.; Leeuwen, P. V.; Evensen, G., Analysis scheme in the ensemble Kalman filter, Mon. Weather Rev., 126, 1719-1724 (1998) [3] Cohn, S. E., An introduction to estimation theory, J. Meteor. Soc. Jpn., 75, 257-288 (1997) [4] Cools, R., An encyclopaedia of cubature formulas, J. Complex., 19, 445-453 (2003) · Zbl 1061.41020 [5] Davis, P. J., Interpolation and Approximation (1975), Dover · Zbl 0111.06003 [6] Engels, H., Numerical Quadrature and Cubature (1980), Academic Press · Zbl 0435.65013 [7] Evensen, G., Sequential data assimilation with a nonlinear quasi-geostrophic model using Monte Carlo methods to forecast error statistics, J. Geophys. Res., 99, 10143-10162 (1994) [8] Evensen, G., The ensemble Kalman filter: theoretical formulation and practical implementation, Ocean Dyn., 53, 343-367 (2003) [9] Evensen, G., Sampling strategies and square root analysis schemes for the EnKF, Ocean Dyn., 54, 539-560 (2004) [10] Evensen, G., Data Assimilation, the Ensemble Kalman Filter (2007), Springer-Verlag: Springer-Verlag Berlin · Zbl 1157.86001 [11] Gelb, A., Applied Optimal Estimation (1974), MIT Press: MIT Press Cambridge [12] Haber, S., Numerical evaluation of multiple integrals, SIAM Rev., 12, 4, 481-526 (1970) · Zbl 0206.46905 [13] Jazwinski, A. H., Stochastic Processes and Filtering Theory (1970), Academic Press: Academic Press San Diego, CA · Zbl 0203.50101 [14] Kalman, R.; Bucy, R., New results in linear prediction and filter theory, Trans. AMSE J. Basic Engrg., 83D, 85-108 (1961) [15] L. Nerger, Parallel filter algorithms for data assimilation in oceanography, Reports on Polar and Marine Research 487. Ph.D. Thesis, Alfred Wegener Institute for Polar and Marine Research, University of Bremen, Bremerhaven, Germany, 2004.; L. Nerger, Parallel filter algorithms for data assimilation in oceanography, Reports on Polar and Marine Research 487. Ph.D. Thesis, Alfred Wegener Institute for Polar and Marine Research, University of Bremen, Bremerhaven, Germany, 2004. [16] Pham, D. T., Stochastic methods for sequential data assimilation in strongly nonlinear systems, Mon. Weather Rev., 129, 1194-1207 (2001) [17] Quarteroni, A.; Sacco, R.; Saleri, F., Numerical Mathematics (2000), Springer-Verlag: Springer-Verlag New York · Zbl 0943.65001 [18] Stroud, A. H., Approximate Calculation of Multiple Integrals (1971), Prentice-Hall: Prentice-Hall Englewood Cliffs, NJ · Zbl 0379.65013 [19] Tippett, M. K.; Anderson, J. L.; Bishop, C. H.; Hamill, T. M.; Whitaker, J. S., Ensemble square-root filters, Mon. Weather Rev., 131, 1485-1490 (2003) [20] Whitaker, J. S.; Hamill, T. M., Ensemble data assimilation without perturbed observations, Mon. Weather Rev., 130, 1913-1924 (2002) [21] Xiu, D., Numerical integration formulas of degree two, Appl. Numer. Math. (2007)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.