Error Covariance Modeling in the GMAO Ocean Ensemble Kalman Filter

Abstract

In practical applications of the ensemble Kalman filter (EnKF) for ocean data assimilation, the computational burden and memory limitations usually require a trade-off between ensemble size and model resolution. This is certainly true for the NASA Global Modeling and Assimilation Office (GMAO) ocean EnKF used for ocean climate analyses. The importance of resolution for the adequate representation of the dominant current systems means that small ensembles, with their concomitant sampling biases, have to be used. Hence, strategies have been sought to address sampling problems and to improve the performance of the EnKF for a given ensemble size. Approaches assessed herein consist of spatiotemporal filtering of background-error covariances, improving the system-noise representation, imposing a steady-state error covariance model, and speeding up the analysis by performing the most expensive operation of the analysis on a coarser computational grid. A judicious combination of these approaches leads to significant performance improvements, especially with very small ensembles.