A Systematic Comparison of Particle Filter and EnKF in Assimilating Time‐Averaged Observations. Issue 24 (20th December 2017)

• Record Type: Journal Article
• Authors: Liu, Huaran; Liu, Zhengyu; Lu, Feiyu
• Is Part Of: Journal of geophysical research. Volume 122:Issue 24(2017)
• Page Start: 13, 155
• Page End: 13, 173
• Publication Date: 2017-12-20
• Subjects: paleoclimate data assimilation -- particle filter -- ensemble Kalman filter; Atmospheric physics -- Periodicals; Geophysics -- Periodicals; 551.5
• DOI: 10.1002/2017JD026798
• Languages: English
• ISSNs: 2169-897X

Abstract

Abstract: The particle filter (PF) and the ensemble Kalman filter (EnKF) are two promising and popularly adopted types of ensemble‐based data assimilation methods for paleoclimate reconstruction. However, no systematic comparison between them has been attempted. We compare these two uncertainty based methods in pseudoproxy experiments where synthetic seasonal mean sea surface temperature observations are assimilated. Their skills are evaluated with regards to local, hemispherically averaged and globally averaged analysis error, and their ability to capture large‐scale modes of variability. It is found that the EAKF (Ensemble Adjustment Kalman filter, a variant of EnKF) performs better than the PF with only one third of the ensemble size, despite PF's theoretical superiority in allowing for non‐Gaussian statistics and nonlinear dynamics. The success of the EAKF is attributed to the facts that (1) Gaussian assumption is somewhat appropriate for this application; (2) The EAKF is less sensitive to sampling errors than the PF due to the different methodological natures. Sixteen members are enough to estimate accurate covariance for the EAKF, but 48 (even 96) members still underrepresent the state space of high‐dimensional system for the PF. Our study highlights the importance of a large localization radius in the application of the EnKF to paleoclimate reconstruction due to the sparse proxy network and suggests that additional techniques, such as localization or clustered particle filter, are needed to improve the PF for paleoclimate reconstruction, in addition to the simple importance resampling currently adopted by most research. Key Points: Pseudoproxy experiments are conducted using particle filter with simple importance resampling and ensemble Kalman filter to evaluate the performance for paleoclimate reconstruction We evaluated the strengths and weaknesses of these two filters for paleoclimate data assimilation In sparse observation network, large localization radius is crucial to EnKF