Journal cover Journal topic
Nonlinear Processes in Geophysics An interactive open-access journal of the European Geosciences Union

Journal metrics

  • IF value: 1.321 IF 1.321
  • IF 5-year<br/> value: 1.636 IF 5-year
  • SNIP value: 0.903 SNIP 0.903
  • SJR value: 0.709 SJR 0.709
  • IPP value: 1.455 IPP 1.455
  • h5-index value: 20 h5-index 20
Nonlin. Processes Geophys., 23, 391-405, 2016
© Author(s) 2016. This work is distributed
under the Creative Commons Attribution 3.0 License.
Research article
04 Nov 2016
A local particle filter for high-dimensional geophysical systems
Stephen G. Penny1,2,3 and Takemasa Miyoshi3,1 1Department of Atmospheric and Oceanic Science, University of Maryland, College Park, MD, USA
2National Centers for Environmental Prediction, College Park, MD, USA
3RIKEN Advanced Institute for Computational Science, Kobe, Japan
Abstract. A local particle filter (LPF) is introduced that outperforms traditional ensemble Kalman filters in highly nonlinear/non-Gaussian scenarios, both in accuracy and computational cost. The standard sampling importance resampling (SIR) particle filter is augmented with an observation-space localization approach, for which an independent analysis is computed locally at each grid point. The deterministic resampling approach of Kitagawa is adapted for application locally and combined with interpolation of the analysis weights to smooth the transition between neighboring points. Gaussian noise is applied with magnitude equal to the local analysis spread to prevent particle degeneracy while maintaining the estimate of the growing dynamical instabilities. The approach is validated against the local ensemble transform Kalman filter (LETKF) using the 40-variable Lorenz-96 (L96) model. The results show that (1) the accuracy of LPF surpasses LETKF as the forecast length increases (thus increasing the degree of nonlinearity), (2) the cost of LPF is significantly lower than LETKF as the ensemble size increases, and (3) LPF prevents filter divergence experienced by LETKF in cases with non-Gaussian observation error distributions.

Citation: Penny, S. G. and Miyoshi, T.: A local particle filter for high-dimensional geophysical systems, Nonlin. Processes Geophys., 23, 391-405, doi:10.5194/npg-23-391-2016, 2016.
Publications Copernicus
Short summary
Particle filters in their basic form have been shown to be unusable for large geophysical systems because the number of required particles grows exponentially with the size of the system. We have applied the ideas of localized analyses at each model grid point and use ensemble weight smoothing to blend each local analysis with its neighbors. This new local particle filter (LPF) makes large geophysical applications tractable for particle filters and is competitive with a popular EnKF alternative.
Particle filters in their basic form have been shown to be unusable for large geophysical...