Marine gravity field recovery relies heavily on satellite altimetry. Thanks to the evolution of altimetry missions and the improvements in altimeter data processing methods, the marine gravity field model has been prominently enhanced in accuracy and resolution. However, high-accuracy and high-resolution gravity field recovery from satellite altimeter data remains particularly challenging. We provide an overview of advances in satellite altimetry for marine gravity field recovery, focusing on the impact factors and available models of altimetric gravity field construction. Firstly, the evolution of altimetry missions and the contribution to gravity field recovery are reviewed, from the existing altimetry missions to the future altimetry missions. Secondly, because the methods of altimeter data processing are of great significance when obtaining high-quality sea surface height observations, these improved methods are summarized and analyzed, especially for coastal altimetry. In addition, the problems to be resolved in altimeter data processing are highlighted. Thirdly, the characteristics of gravity recovery methods are analyzed, including the inverse Stokes formula, the inverse Vening Meinesz formula, Laplace’s equation, and least squares collocation. Furthermore, the latest global marine gravity field models are introduced, including the use of altimeter data and processing methods. The performance of the available global gravity field model is also evaluated by shipboard gravity measurements. The root mean square of difference between the available global marine gravity model and shipboard gravity from the National Centers for Environmental Information is approximately 5.10 mGal in the low-middle latitude regions, which is better than the result in high-latitude regions. In coastal areas, the accuracy of models still needs to be further improved, particularly within 40 km from the coastline. Meanwhile, the SDUST2021GRA model derived from the Shandong University of Science and Technology team also exhibited an exciting performance. Finally, the future challenges for marine gravity field recovery from satellite altimetry are discussed.

The resolutions of current global altimetric gravity models and mean sea surface models are around 12 km wavelength resolving 6 km features, and for many years it has been difficult to improve the resolution further in a systematic way. For both Jason 1 and 2, a Geodetic Mission (GM) has been carried out as a part of the Extension-of-Life phase. The GM for Jason-1 lasted 406 days. The GM for Jason-2 was planned to provide ground-tracks with a systematic spacing of 4 km after 2 years and potentially 2 km after 4 years. Unfortunately, the satellite ceased operation in October 2019 after 2 years of Geodetic Mission but still provided a fantastic dataset for high resolution gravity recovery. We highlight the improvement to the gravity field which has been derived from the 2 years GM. When an Extension-of-Life phase is conducted, the satellite instruments will be old. Particularly Jason-2 suffered from several safe-holds and instrument outages during the GM. This leads to systematic gaps in the data-coverage and degrades the quality of the derived gravity field. For the first time, the Jason-2 GM was “rewound” to mitigate the effect of the outages, and we evaluate the effect of “mission rewind” on gravity. With the recent successful launch of Sentinel-6 Michael Freilich (S6-MF, formerly Jason CS), we investigate the possibility creating an altimetric dataset with 2 km track spacing as this would lead to fundamental increase in the spatial resolution of global altimetric gravity fields. We investigate the effect of bisecting the ground-tracks of existing GM to create a mesh with twice the resolution rather than starting all over with a new GM. The idea explores the unique opportunity to inject Jason-3 GM into the same orbital plane as used for Jason-2 GM but bisecting the existing Jason-2 tracks. This way, the already 2-years Jason-2 GM could be used to create a 2 km grid after only 2 years of Jason-3 GM, rather than starting all over with a new GM for Jason-3.

HaiYang-2A (HY-2A, where ‘Haiyang’ means ‘Ocean’ in Chinese) has provided reliable sea surface height observations for gravity with uniform ocean data coverage on a global scale for more than 8 years, particularly with denser across track sampling during the geodetic mission since March 2016. This paper aims at modeling and evaluating the regional marine gravity field at 1′×1′ resolution by incorporating HY-2A altimeter waveform data from 7 complete 168-day cycles in the geodetic mission phase. Initial evaluation indicates that, firstly, the measurements in the geodetic mission stay at a consistent accuracy level with observations at the start-of-life stage according to statistics of discrepancies at crossover points cycle by cycle. Secondly, range precision improvement can be achieved using a two-pass weighted least-squares retracker. Thirdly, a downsampling procedure combined with a low-pass filter is designed for HY-2A 20 Hz data to obtain 5 Hz measurements with enhanced precision. We calculate the 1′×1′ marine gravity field model over the South China Sea area by using the EGM2008 model as a reference field with the remove/restore method. The verifications with published models and shipborne gravimetric data show that HY-2A GM data is capable of improving marine gravity field modeling. Results show slightly higher accuracy than other models with similar input datasets but not including HY-2A. The accuracy is also compared with the latest DTU17 and SIO V27.1 model.

Haiyang-2 (HY-2) missions have accumulated sea surface height (SSH) observations on a global scale for more than 10 years. Four satellites, HY-2A, HY-2B, HY-2C and HY-2D, provide even but differently distributed data, which play a complementary role in marine gravity studies with other missions. Therefore, this paper evaluates the performances of HY-2 altimetric data in marine gravity modeling from the following four perspectives: SSH accuracy, geoid signal resolution ability, vertical deflections and gravity anomaly. First, the centimeter-magnitude accuracy level of HY-2 data is proved by analyzing SSH discrepancies at crossover points within a certain time limit. Second, the spectral analysis of repetitive along-track data sequences in a time domain shows a geoid resolution range from 18 to 24 km. Taking HY-2 exact repeat missions (ERM), for example, the resolution could be remarkably enhanced by stacking repetitive cycles. Third, validation with an XGM2019 model showed that vertical deflections were reliably computed for all HY-2 missions, but HY-2A performed slightly worse than the other HY-2 missions. Meanwhile, HY-2C and HY-2D with a ~66° orbital inclination obviously had an improved ability to capture east–west signals compared to HY-2A and HY-2B. Finally, we constructed global marine gravity results based on three input datasets, HY-2 dataset only, multi-satellite dataset without HY-2 and multi-satellite dataset with HY-2. Validations were performed using published models and shipborne gravimetric data. The results showed that the HY-2 dataset is capable of improving marine gravity anomaly recoveries and that the accuracy of NSOAS22 with incorporated HY-2 data is comparable to DTU21 and SS V31.1. Furthermore, HY-2 observations should not be the only input dataset to construct a 1’ × 1’ resolution marine gravity model.