Improving particle filter performance by smoothing observations

TL;DR

This paper proposes smoothing observations by increasing small-scale observation variance to prevent particle filter collapse, thereby enhancing uncertainty quantification and accuracy in spatially-extended systems.

Contribution

It introduces a novel observation error model that smooths the posterior mean without requiring ensemble smoothing, improving particle filter stability and performance.

Findings

Reduces particle filter collapse significantly.

Improves continuous ranked probability scores by up to 25%.

Maintains accuracy while enhancing uncertainty quantification.

Abstract

This article shows that increasing the observation variance at small scales can reduce the ensemble size required to avoid collapse in particle filtering of spatially-extended dynamics and improve the resulting uncertainty quantification at large scales. Particle filter weights depend on how well ensemble members agree with observations, and collapse occurs when a few ensemble members receive most of the weight. Collapse causes catastrophic variance underestimation. Increasing small-scale variance in the observation error model reduces the incidence of collapse by de-emphasizing small-scale differences between the ensemble members and the observations. Doing so smooths the posterior mean, though it does not smooth the individual ensemble members. Two options for implementing the proposed observation error model are described. Taking discretized elliptic differential operators as an observation error covariance matrix provides the desired property of a spectrum that grows in the approach to small scales. This choice also introduces structure exploitable by scalable computation techniques, including multigrid solvers and multiresolution approximations to the corresponding integral operator. Alternatively the observations can be smoothed and then assimilated under the assumption of independent errors, which is equivalent to assuming large errors at small scales. The method is demonstrated on a linear stochastic partial differential equation, where it significantly reduces the occurrence of particle filter collapse while maintaining accuracy. It also improves continuous ranked probability scores by as much as 25%, indicating that the weighted ensemble more accurately represents the true distribution. The method is compatible with other techniques for improving the performance of particle filters.