An Integrated Method for Dynamic Height Error Correction in GNSS-IR Sea Level Retrievals

Sea level is an important variable for studying water cycle and coastal hazards under global warming. Global Navigation Satellite System Interferometric Reflectometry (GNSS-IR) has emerged as a relatively new technique for monitoring sea level variations, leveraging signals from GNSS constellations. However, dynamic height errors, primarily caused by non-stationary sea surfaces, compromise the precision of GNSS-IR sea level retrievals and necessitate robust correction. In this study, we propose a new method to correct the dynamic height error by integrating the commonly used tidal analysis method and the cubic spline fitting method. The proposed method is applied to the GNSS-IR sea level retrievals from multiple systems and multiple frequency bands at two coastal GNSS stations, MAYG and HKQT. At MAYG, the results show that our method significantly reduces the Root Mean Square Error (RMSE) of the GNSS-IR sea level retrievals by 42.1% (11.4 cm) to 15.7 cm, performing better than the single tidal analysis method (16.5 cm) and the cubic spline fitting method (21.4 cm). At HKQT, our method improves the accuracy by 21.5% (3.1 cm) to 10.3 cm, which is still better than that of the tidal analysis method (11.3 cm) and the cubic spline fitting method (12.4 cm). Compared to the tidal analysis method and the cubic spline fitting method, our method maintains high retrieval retention while enhancing precision. The effectiveness of our method is further validated in the two storm surge events caused by Typhoon Hato and Typhoon Mangkhut in Hong Kong.