A novel Slepian approach for determining mass-term sea level from GRACE over the South China Sea

Treatment of land-to-ocean leakage is crucial in modeling the mass-term sea level (MSL) using the Gravity Recovery and Climate Experiment (GRACE) satellite gravity measurements. In this study, we utilized the Spherical Slepian functions (SSFs) to determine MSL variations in the South China Sea. A sensitivity simulation in terms of trend, annual amplitude, and phase revealed that the SSF solution with a 1° coastal buffer zone provides a better land-to-ocean leakage correction than traditionally used Spherical Harmonics (SH). It was also found that an additional smoothing procedure for SSF with low concentrating energy could significantly reduce the high-frequency noise in GRACE (e.g., the north–south strips). The spatiotemporal characteristics of the true MSL were further compared among GRACE SH and Mascon solutions, model-predicted ocean bottom pressure, and steric-corrected altimetric data (i.e., satellite-altimetric sea level minus steric effect). Results revealed that despite the SSF-inverted regional MSL solution being generally similar to other results, this technique notably recovered the realistic magnitude of detailed features. The apparent MSL signal was, for instance, precisely recovered in the central part of the Gulf of Thailand and the Sunda Shelf by the SSF, while the SH significantly underestimated it due to smoothing. In addition to the seasonality, the interannual signal decomposed from the SSF-inverted residual MSL (by removing the trend and periodic signals) was the strongest and well cross-correlated with the Southern Oscillation Index. With more detailed spatial patterns revealed by MSL from GRACE SSF, our findings demonstrate that SSF is more suitable for regional-scale studies.