An Effcient Implementation of the Ensemble Kalman Filter Based on Iterative Sherman Morrison Formula
Author(s) -
Elías D. Niño,
Adrian Sandu,
J. G. Anderson
Publication year - 2012
Publication title -
procedia computer science
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.334
H-Index - 76
ISSN - 1877-0509
DOI - 10.1016/j.procs.2012.04.115
Subject(s) - computer science , ensemble kalman filter , kalman filter , computation , data assimilation , algorithm , covariance , extended kalman filter , iterative method , covariance matrix , inverse , mathematical optimization , artificial intelligence , mathematics , statistics , physics , geometry , meteorology
This paper proposes an effcient implementation of the ensemble Kalman filter (EnKF) for the solution of largescale data assimilation problems. The implementation exploits the special structure of the covariance matrix and solves the analysis step by iteratively applying the Sherman Morrison formula. The iterative implementation leads to p savings in both memory and run time. The number of operations for the iterative method is O(nens2.nobs, while for p 3the standard implementation the cost is O nobs. The new implementation of the EnKF is tested using the Lorenz 96 model and shows a better performance than EnKF with the direct computation of the inverse
Accelerating Research
Robert Robinson Avenue,
Oxford Science Park, Oxford
OX4 4GP, United Kingdom
Address
John Eccles HouseRobert Robinson Avenue,
Oxford Science Park, Oxford
OX4 4GP, United Kingdom