Spatially and temporally varying adaptive covariance inflation for ensemble filters

Abstract

Ensemble filters are used in many data assimilation applications in geophysics. Basic implementations of ensemble filters are trivial but are susceptible to errors from many sources. Model error, sampling error and fundamental inconsistencies between the filter assumptions and reality combine to produce assimilations that are suboptimal or suffer from filter divergence. Several auxiliary algorithms have been developed to help filters tolerate these errors. For instance, covariance inflation combats the tendency of ensembles to have insufficient variance by increasing the variance during the assimilation. The amount of inflation is usually determined by trial and error. It is possible, however, to design Bayesian algorithms that determine the inflation adaptively. A spatially and temporally varying adaptive inflation algorithm is described. A normally distributed inflation random variable is associated with each element of the model state vector. Adaptive inflation is demonstrated in two low-order model experiments. In the first, the dominant error source is small ensemble sampling error. In the second, the model error is dominant. The adaptive inflation assimilations have better mean and variance estimates than other inflation methods.