Filling the Data Gaps Within GRACE Missions Using Singular Spectrum Analysis

Dozens of missing epochs in the monthly gravity product of the satellite mission Gravity Recovery and Climate Experiment (GRACE) and its follow-on (GRACE-FO) mission greatly inhibit the complete analysis and full utilization of the data. Despite previous attempts to handle this problem, a general all-purpose gap-filling solution is still lacking. Here we propose a non-parametric, data-adaptive and easy-to-implement approach-composed of the Singular Spectrum Analysis gap-filling technique, cross-validation, and spectral testing for significant components-to produce reasonable gap-filling results in the form of spherical harmonic coefficients (SHCs). We demonstrate that this approach is adept at inferring missing data from long-term and oscillatory changes extracted from available observations. A comparison in the spectral domain reveals that the gap-filling result in the same signal intensity as the product of GRACE missions for spherical harmonic degrees below 30. As the degree increases above 30, the degree variance of the gap-filling result decreases more rapidly than that of GRACE/GRACE-FO SHCs, demonstrating effective suppression of noise. As a result, our approach can reduce noise in the oceans without sacrificing resolutions on land. The quality of the gap-filling product is evaluated through comparison with a surface mass balance based estimate in Greenland, Swarm gravity solutions and a climate-driven water storage model, all of which confirm the good performance of our results. This study makes a ready-to-use gap-filling product in the form of SHCs together with proper error estimates available. Nonetheless, our method is also applicable to smoothed gridded observations.

Automated summary

This study presents a non-parametric, data-adaptive method to fill the missing epochs in the monthly gravity product of satellite missions GRACE and GRACE-FO. Comparison in the spectral domain shows that the method achieves the same signal intensity as GRACE missions for spherical harmonic degrees below 30 and effectively suppresses noise as the degree increases above 30.