Local Enhancement of Marine Gravity Field over the Spratly Islands by Combining Satellite SAR Altimeter-Derived Gravity Data

: The marine gravity ﬁeld recovery close to land/island is challenging owing to the scarcity of measured gravimetric observations and sorely contaminated satellite radar altimeter-derived data. The satellite missions that carried the synthetic aperture radar (SAR) altimeters supplied data with improved quality compared to that retrieved from the conventional radar altimeters. In this study, we combine the satellite altimeter-derived gravity data for marine gravity ﬁeld augmentation over island areas; in particular, the feasibility for regional augmentation by incorporating the SAR altimeter-derived gravity data is investigated. The gravity ﬁeld modeling results over the Spratly Islands demonstrate that the marine gravity ﬁeld is augmented by the incorporation of newly published satellite altimeter-derived gravity data. By merging the gravity models computed with the Sentinel-3A/B SAR altimetry data, the quasi-geoid and mean dynamic topography are dramatically improved, by a magnitude larger than 4 cm around areas close to islands, in comparison with the results directly derived from a combined global geopotential model alone. Further comparison of regional solutions computed from heterogeneous gravity models shows that the ones modeled from the SAR-based gravity models have better performances, the errors of which are reduced by a magnitude of 2~4 cm over the regions close to islands, in comparison with the solutions modeled with the gravity models developed without SAR altimetry data. These results highlight the superiority of using the SAR-based gravity data in marine gravity ﬁeld recovery, especially over the regions close to land/island.


Introduction
High-resolution gravity field determination at seas is a basic task in geodesy. Thanks to the Gravity Recovery and Climate Experiment (GRACE) [1,2] and Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) missions [3,4], the global gravity field has been prominently strengthened at long wavelength up to hundreds of kilometers [5][6][7][8].
On the other hand, by incorporating ground-based data at short-wavelength, the derived Global Geopotential Models (GGMs) (known as high-degree or combined GGMs) can map the gravity signals at a mean spatial resolution of 5 arc-minutes (~10 km) in global scale [9][10][11].
Despite the tremendous progresses made in global geopotential model computation over decades, the lack of globally distributed gravity data is still a major barrier to improve the combined GGMs. For regions inland, such as most areas in Asia and Africa, the measured land/airborne gravimetric observations are limited; accessible and fill-in measurements were involved in the computation of combined GGMs [9,12,13], whereas satellite altimetric gravity data was usually used over oceans as shipborne/airborne gravimetric measurements were usually inaccessible or confidentially kept due to political reasons [14]. However, the satellite altimeter-derived data is notoriously known to be of low quality close to land/island, owing to the contamination of radar altimeter waveforms and degraded quality of geophysical models used for data corrections [15][16][17][18]. Inevitably, the errors in satellite altimeter-derived gravimetric data were spread to the combined GGMs. As a result, the errors in the combined GGMs reached 6~10 centimeters over most oceans [9], and even a magnitude of decimeter level or larger over polar areas and coastal regions [19,20]. The errors in a GGM may bring about dramatic disturbances in the investigation of the mean ocean state in detail [21,22].
The improved satellite altimetry techniques result in the augmentation of global marine gravity field, by a factor of 2~4, in comparison with the gravity field models developed with old altimeter-derived records [23][24][25]. In particular, the CryoSat-2 and Sentinel-3A/B missions that carried synthetic aperture radar (SAR) altimeters derived precise records of sea surface height up to several kilometers from the coast [17,26]. The CryoSat-2 applies the low-resolution mode (LRM), SAR mode, and SAR interferometry (SARIn) mode over different areas [27][28][29]. Sentinel-3A/B inherited the SAR altimeter from CryoSat-2, and it retrieved observations with denser spatial coverage in comparison with the conventional altimetry (or operated in the LRM mode). Moreover, the data derived from the SAR waveforms have higher signal-to-noise ratio and lower speckle noise [30][31][32], and the accuracy of SAR altimetry data reaches several centimeters to decimeter level close to coast and lakes [33][34][35]. By incorporating recent altimetry observations, the derived gravimetric products have improved precision versus the ones developed without these observations [36].
The launch of the satellite altimetry missions carrying the SAR altimeter provides solid basis for local gravity field recovery close to land/island; however, little attention has been paid to gravity field enhancement over coastal or island areas by combing the SAR altimetry data, especially for the use of recently released data from Sentinel-3A/B. To our best knowledge, no existing literature has investigated and quantified the additional signals introduced from the SAR-based data retrieved from Sentinel-3A/B on regional gravity field modeling. This study aims at strengthening the marine gravity field over regions close to land/island on a regional scale based on gravity data derived from satellite altimeters. In particular, we study the feasibility for local augmentation based on SARbased gravity data. Moreover, we compare the performances of different altimetric gravity models in marine gravity field recovery. In the following, the study area and data sources are included in Section 2. Section 3 displays the gravity field modeling results, where the added signals retrieved from the newly published satellite altimeter-derived gravity data are quantified and validated; in particular, the possibility for local enhancement based on SAR-based gravity data is investigated. Section 4 contains a summary of the study and the main conclusions.

Study Area and Data
We chose the study area as the Spratly Islands located in the southern part of the South China Sea (SCS) (see the area inside the red rectangle in Figure 1a), and the bathymetry is retrieved from the General Bathymetric Chart of the Oceans (GEBCO) [37]. The Spratly Islands are the largest archipelago in the SCS and are also the southernmost archipelago, where small islands, sandbanks, shoals, and coral reefs are the predominant structures [38] (see Figure 1b). The Spratly Islands archipelago is rich in natural resources, and is also a disputed archipelago with complicate governances. As a result, it is hard to implement shipborne/airborne gravimetric surveys over this area. The lack of measured gravimetric observations brings about the difficulty for marine gravity field recovery over the Spratly Islands, and the current combined GGMs were computed by using satellite altimetric Remote Sens. 2022, 14, 474 3 of 19 gravity data. Moreover, small islands, atolls, sandbanks, shoals, and reefs are in abundance in the Spratly Islands, and the return waveforms from radar altimeters have been severely contaminated; as such, there exist multiple challenges in local marine gravity refinement. However, this offers a chance to research the feasibility of using SAR-based gravimetric data in regional augmentation.
Remote Sens. 2022, 14, x FOR PEER REVIEW 3 of hard to implement shipborne/airborne gravimetric surveys over this area. The lack measured gravimetric observations brings about the difficulty for marine gravity fie recovery over the Spratly Islands, and the current combined GGMs were computed using satellite altimetric gravity data. Moreover, small islands, atolls, sandbanks, shoa and reefs are in abundance in the Spratly Islands, and the return waveforms from rad altimeters have been severely contaminated; as such, there exist multiple challenges local marine gravity refinement. However, this offers a chance to research the feasibil of using SAR-based gravimetric data in regional augmentation. As there were no measured gravimetric observations available over the Spratly lands, we use the satellite altimeter-derived data for marine gravity field modeling. F this purpose, several recently released altimetric gravity models are used. First, t widely used models developed at the Technical University of Denmark (DTU) space a Scripps Institution of Oceanography (SIO) are introduced. For the series of models d veloped at DTU space, the latest model, namely DTU21GRA, and several predecesso i.e., DTU13GRA, DTU15GRA, and DTU17GRA [36,39], are used. Typically, for the mo els that were developed in the DTU space, the updated gravity models had improv accuracy versus the previous versions, which were mainly due to the use of newly leased altimeter data and the updated data preprocessing methods. For instance, t comparison with an airborne survey over North Greenland showed that the standa deviation (SD) of the misfits between DTU17GRA and the airborne data was 3.78 mG This value increased to 8.81, 5.91, and 5.45 mGal when DTU10GRA, DTU13GRA, a DTU15GRA were evaluated, respectively. The major improvement of DTU17GRA ov DTU15GRA is that the computation of DTU17GRA contained more CryoSat-2 data a SARAL/AltiKa data from 2016 to 2017 in the geodetic phase, while the improvement DTU21GRA over DTU17GRA is that the former was computed by combining 5 years Sentinel-3A and 3 years of Sentinel-3B and reprocessed Cryosat-2 data (processed w the SAMOSA+ physical retracker).
For the models developed in SIO, the latest model, called SIO V31.1, and its previo version, namely, SIO V30.1 (hereinafter referred to as SIO31 and SIO30, respective [23,24,40], are introduced. In addition, two recently published regional models over t SCS, i.e., SCSGA V1.0 [41] and HY-2A V1.0 [42], are used. In short, we call these two gional models SCSGA and HY2A, respectively. All these satellite altimetric grav models have the spatial resolutions of one arc-minute, and the datasets used in the models' development and other associated information are shown in Table 1. As there were no measured gravimetric observations available over the Spratly Islands, we use the satellite altimeter-derived data for marine gravity field modeling. For this purpose, several recently released altimetric gravity models are used. First, the widely used models developed at the Technical University of Denmark (DTU) space and Scripps Institution of Oceanography (SIO) are introduced. For the series of models developed at DTU space, the latest model, namely DTU21GRA, and several predecessors, i.e., DTU13GRA, DTU15GRA, and DTU17GRA [36,39], are used. Typically, for the models that were developed in the DTU space, the updated gravity models had improved accuracy versus the previous versions, which were mainly due to the use of newly released altimeter data and the updated data preprocessing methods. For instance, the comparison with an airborne survey over North Greenland showed that the standard deviation (SD) of the misfits between DTU17GRA and the airborne data was 3.78 mGal. This value increased to 8.81, 5.91, and 5.45 mGal when DTU10GRA, DTU13GRA, and DTU15GRA were evaluated, respectively. The major improvement of DTU17GRA over DTU15GRA is that the computation of DTU17GRA contained more CryoSat-2 data and SARAL/AltiKa data from 2016 to 2017 in the geodetic phase, while the improvement of DTU21GRA over DTU17GRA is that the former was computed by combining 5 years of Sentinel-3A and 3 years of Sentinel-3B and reprocessed Cryosat-2 data (processed with the SAMOSA+ physical retracker).
For the models developed in SIO, the latest model, called SIO V31.1, and its previous version, namely, SIO V30.1 (hereinafter referred to as SIO31 and SIO30, respectively) [23,24,40], are introduced. In addition, two recently published regional models over the SCS, i.e., SCSGA V1.0 [41] and HY-2A V1.0 [42], are used. In short, we call these two regional models SCSGA and HY2A, respectively. All these satellite altimetric gravity models have the spatial resolutions of one arc-minute, and the datasets used in these models' development and other associated information are shown in Table 1. The modeling of DTU13GRA additionally involved data from CryoSat-2 and Jason-1, compared to its predecessor, i.e., DTU10GRA, which was derived by combining data from Topex/Poseidon (T/P), Geosat, ICESat, Jason-1, GFO, Envisat, and retracked data from ERS-1. This model was developed by incorporating data from HY-2A and data from T/P, Geosat, Jason-1, Envisat, ERS-1/2, CryoSat-2, and SARAL/AltiKa.

Yes
For the models developed at the DTU space, it is noticeable that although DTU13GRA/ DTU15GRA/DTU17GRA was computed with CryoSat-2 data, CryoSat-2 operated in the LRM over the SCS, and no SAR altimetry data was used in the computation of these models. For the same reason, SCSGA and HY2A are the models that were computed without combining SAR altimetry data, whereas the development of SIO30, SIO31, and DTU21GRA incorporated Sentinel-3A/B SAR altimeter measurements over the SCS, and the altimeter in the Sentinel-3A/B mission operated in the SAR mode all around the world.
The comparisons with an airborne survey over North Greenland and the marine gravity data from the NGA showed that DTU17GRA had improved accuracy, compared to the previous versions, e.g., DTU10GRA, DTU13GRA, and DTU15GRA [36]. Moreover, the validation with independent shipborne gravity data over the SCS displayed that the SCSGA model had comparable precision to SIO V27.1 but had better performance than DTU13GRA/DTU17GRA [41], whereas HY2A had degraded precision compared with SIO V27.1 and DTU17GRA [42].

Local Refinement with SAR-Based Altimetric Gravity Data
The remove-compute-restore approach is applied for modeling [43,44], and we use the residual terrain model (RTM) to reduce the high-frequency topographical signals [45]. XGM2019e_2159, with a full degree and order (d/o) of 2190/2159, is chosen as the reference model [46]. A satellite-only GGM, namely, GOCO06S, was used to represent the longwavelength component of XGM2019e_2159, and the gravimetric data sources from National Geospatial-Intelligence Agency (NGA) and DTU13GRA were used to compute the shortwavelength component of this model [46]. The validation against local airborne gravimetric observations demonstrated that XGM2019e_2159 had improved precision, by a magnitude of~1 mGal, compared to the GGMs that have similar expansion degrees [47].
We model the residual gravity field by using Poisson wavelets, and the long-and short-wavelength components of gravity field are recovered from the reference model and RTM, respectively [48,49]. To compute the residual gravity data, the gravity anomalies synthesized from XGM2019e_2159 up to degree 2190 and the associated RTM corrections are subtracted from the satellite altimeter-derived gravity data. Then, we parameterized the residual gravity field based on Poisson wavelets, and the function that links the gravity anomaly to the disturbing potential is seen in Klee et al. (2008) [50]. The area extends 7 • N to 12 • N in the latitudinal direction, and 111.0 • E to 116 • E in the longitudinal direction is chosen as the computational region. We locate the Poisson wavelets 5 km beneath the terrain, and we set the mean resolution of Poisson wavelets as 10 km. Moreover, we compute the quasi-geoid model with a spatial resolution of~2 km, in correspondence with the mean resolution of the satellite altimeter-derived gravity models.
To highlight the use of SAR altimetry data in local augmentation, we investigate the performances of various solutions modeled with different altimetric gravity models that are computed with and without SAR altimetry data. To do this, four different versions of altimetric gravity models that were developed at DTU space, i.e., DTU13GRA, DTU15GRA, DTU17GRA, and DTU21GRA, are used. It is of note that only the DTU21GRA model was developed with the SAR altimetry data retrieved from Sentinel-3A/B over the Spratly Islands, while the other three models were developed without SAR altimetry data. The residual quasi-geoids computed with various satellite altimeter-derived gravity models are seen in Figure 2 and reach a magnitude of several centimeters (see the statistics in Table 2). The magnitude of the contents retrieved from DTU13GRA is slightly larger than 3 cm, and the standard deviation value is~0.6 cm. Although the reference model in local gravity field modeling, i.e., XGM2019e_2159, was computed based on the DTU13GRA data over oceans, there may still be the additional contents in the DTU13GRA data that were unresolved in XGM2019e_2159, whereas the minimum/maximum value changes to −3.0/4.0 cm (−3.4/4.8 cm) when DTU15GRA (DTU17GRA) is applied for modeling. The additional signals derived from these satellite altimeter-derived gravity models concentrate over the areas close to islands, e.g., close to the Zhongye Reefs The added signals retrieved from DTU17GRA have stronger structures than those derived from the DTU13GRA data, which may be due to that the computation of DTU17GRA involved more high-quality data, including more data from Jason-1, CryoSat-2, and SARAL/AltiKa, as well as the updated data preprocessing strategies. As a result, DTU17GRA had improved accuracy compared to DTU13GRA [36]. timetry data retrieved from Sentinel-3A/B may further contribute to the marine gravity field augmentation, especially over island areas, compared to the altimetric gravity models computed without SAR altimeter data, as the SAR waveforms have the relatively high signal-to-noise ratio and can alleviate the well-known coast problem. Moreover, the SAR altimeters supply observations with denser spatial coverage than conventional radar altimeters, and the incorporation of DTU21GRA is beneficial to recover the short-wavelength component of local gravity field.

Max
Min Further, the geodetic mean dynamic topography (MDT) solutions based on the quasi-geoids derived from different altimetric gravity data are computed and compared. The MDT represents the departure of the MSS from the geoid [51][52][53][54], and we use the quasi-geoid rather than the geoid, due to fact that the quasi-geoid coincides with the geoid at seas [55]. We choose DTU21MSS as the MSS model. For the derivation of DTU21MSS, more than 5 years of Sentinel-3A data and 2 years of Sentinel-3B data were combined [56]. Moreover, an updated waveform retracker, i.e., the SAMOSA+ physical retracker, was used to preprocess the Cryosat-2 data. This further improves the quality of  Table 2. Statistics of the residual quasi-geoid heights derived from various satellite altimeter-derived gravity models developed at DTU space (units: cm).

Max
Min Mean SD In contrast, the added signals range from -5.2 to 5.3 cm (the SD value is~1.0 cm) when DTU21GRA is used. The contents retrieved from DTU21GRA have stronger structures than those derived from the three models above, particularly over the northern part and the regions close to islands. The DTU21GRA data computed with the SAR altimetry data retrieved from Sentinel-3A/B may further contribute to the marine gravity field augmentation, especially over island areas, compared to the altimetric gravity models computed without SAR altimeter data, as the SAR waveforms have the relatively high signal-to-noise ratio and can alleviate the well-known coast problem. Moreover, the SAR altimeters supply observations with denser spatial coverage than conventional radar altimeters, and the incorporation of DTU21GRA is beneficial to recover the short-wavelength component of local gravity field.
Further, the geodetic mean dynamic topography (MDT) solutions based on the quasigeoids derived from different altimetric gravity data are computed and compared. The MDT represents the departure of the MSS from the geoid [51][52][53][54], and we use the quasigeoid rather than the geoid, due to fact that the quasi-geoid coincides with the geoid at seas [55]. We choose DTU21MSS as the MSS model. For the derivation of DTU21MSS, more than 5 years of Sentinel-3A data and 2 years of Sentinel-3B data were combined [56]. Moreover, an updated waveform retracker, i.e., the SAMOSA+ physical retracker, was used to preprocess the Cryosat-2 data. This further improves the quality of the MSS model over coastal regions compared to its previous version, e.g., DTU18MSS [57].
The raw MDTs are seen in Figure 3, and we find that the dramatic oscillations exist in the model derived from XGM2019e_2159 over island areas. The errors of the quasi-geoid computed from a combined GGM spread to the MDT, which possibly brings about errors larger than several centimeters at seas [9], and even to decimeter level over coastal regions and polar areas [19,20]. In comparison, smaller disturbances are seen in the solutions computed by using the quasi-geoids that are augmented by incorporating the satellite altimeter-derived gravity data (see Figure 3b By using the newly published satellite altimeter-derived gravity data, the quasi-geoid and mean dynamic topography are augmented, compared to the GGM-derived solutions. The mutual comparisons demonstrate that the MDT computed with the quasi-geoid enhanced by DTU21GRA shows relatively smooth structures, especially in the regions close to islands. This is probably owing to that the computation of DTU21GRA involved the Sentinel-3A/3B SAR altimetry data, which brings about regional augmentation. the MSS model over coastal regions compared to its previous version, e.g., DTU18MSS [57]. The raw MDTs are seen in Figure 3, and we find that the dramatic oscillations exist in the model derived from XGM2019e_2159 over island areas. The errors of the quasi-geoid computed from a combined GGM spread to the MDT, which possibly brings about errors larger than several centimeters at seas [9], and even to decimeter level over coastal regions and polar areas [19,20]. In comparison, smaller disturbances are seen in the solutions computed by using the quasi-geoids that are augmented by incorporating the satellite altimeter-derived gravity data (see Figure 3b 8°N, 114.5°E). By using the newly published satellite altimeter-derived gravity data, the quasi-geoid and mean dynamic topography are augmented, compared to the GGM-derived solutions. The mutual comparisons demonstrate that the MDT computed with the quasi-geoid enhanced by DTU21GRA shows relatively smooth structures, especially in the regions close to islands. This is probably owing to that the computation of DTU21GRA involved the Sentinel-3A/3B SAR altimetry data, which brings about regional augmentation. The raw MDTs are further investigated along latitude 10.4°N and longitude 114.6°E, and the profiles are seen in Figure 3e. In Figure 4a, we see that the MDT derived from the quasi-geoid augmented by merging DTU21GRA has smoother properties than the one directly derived from XGM2019e_2159 and the solutions computed from other altimetric the quasi-geoid augmented by merging DTU21GRA has smoother properties than the one directly derived from XGM2019e_2159 and the solutions computed from other altimetric gravity models. This is particularly the case when this profile crosses the Zhenghe Reefs (114.6 • E) and Wufang Reef (115.7 • E), where the spike-like patterns appear in the MDTs. These spikes are usually identified as the errors due to the degraded quality of satellite altimeter-derived data close to land/island [58]. Over these regions, the XGM2019e_2159derived solution has the most prominent disturbances; and the inconsistencies between this model and the one strengthened by combining DTU21GRA is greater than 3 cm close to the Zhenghe Reefs. The mutual comparison shows that the MDTs computed from DTU13GRA, DTU15GRA, and DTU17GRA almost have consistent structures, although the computation of DTU17GRA incorporated more CryoSat-2 and SARAL/AltiKa data that was not included in the derivation of DTU13GRA/DTU15GRA. In contrast, the application of the DTU21GRA data reduces these spike-like errors, by a magnitude exceeding 2 cm close to the Zhenghe Reefs and Wufang Reef, compared to the associated MDTs modeled with the altimetric gravity data computed without SAR altimetry data. This also corresponds with the results shown in Figure 3, where the MDT modeled with the DTU21GRA data has relatively smooth structures. These results demonstrate the superiority of using the SAR-based gravity data in the determination of marine quasi-geoid and mean dynamic topography over the areas close to land/island. These spikes are usually identified as the errors due to the degraded quality of satellite altimeter-derived data close to land/island [58]. Over these regions, the XGM2019e_2159-derived solution has the most prominent disturbances; and the inconsistencies between this model and the one strengthened by combining DTU21GRA is greater than 3 cm close to the Zhenghe Reefs. The mutual comparison shows that the MDTs computed from DTU13GRA, DTU15GRA, and DTU17GRA almost have consistent structures, although the computation of DTU17GRA incorporated more CryoSat-2 and SARAL/AltiKa data that was not included in the derivation of DTU13GRA/DTU15GRA. In contrast, the application of the DTU21GRA data reduces these spike-like errors, by a magnitude exceeding 2 cm close to the Zhenghe Reefs and Wufang Reef, compared to the associated MDTs modeled with the altimetric gravity data computed without SAR altimetry data. This also corresponds with the results shown in Figure 3, where the MDT modeled with the DTU21GRA data has relatively smooth structures. These results demonstrate the superiority of using the SAR-based gravity data in the determination of marine quasi-geoid and mean dynamic topography over the areas close to land/island. The properties of the MDTs along longitude 114.6°N are seen in Figure 4b. We also find that the mean dynamic topography computed with the augmentation of The properties of the MDTs along longitude 114.6 • N are seen in Figure 4b. We also find that the mean dynamic topography computed with the augmentation of DTU21GRA has relatively small variations, especially when this profile crosses the Zhenghe Reefs (10.4 • E) and Lesi Shoal (11.35 • E). Moreover, the discrepancy between the XGM2019e_2159-derived MDT and the one strengthened by using the DTU21GRA data exceeds 3 cm (4 cm) over the Zhenghe Reefs (Lesi Shoal). Similarly, compared to the MDT modeled with the altimetric gravity model computed without SAR altimeter data, i.e., DTU13GRA, DTU15GRA, and DTU17GRA, the use of DTU21GRA data can reduce these spike-like errors by a magnitude exceeding 2 cm (3 cm) over the Zhenghe Reefs (Lesi Shoal).

Performances of Heterogeneous Altimetric Gravity Models
Moreover, we investigate the performances of heterogeneous gravity models released from different institutes in marine gravity field modeling. To achieve this, five representative models, i.e., HY2A, SCSGA, SIO30, SIO31, and DTU21GRA, are selected for the case study. The magnitudes of the added signals in terms of quasi-geoid heights derived from different gravity models are not consistent (see Figure 5). The maximum/minimum of the additional signals retrieved from HY2A is -5.5/7.0 cm (see Table 3). In comparison, the signals derived from SCSGA have weaker patterns, and the maximum/minimum value is -4.2/4.4 cm, whereas the added signals computed from the altimetric gravity data modeled with the Sentinel-3A/B SAR altimeter data, i.e., SIO30/SIO31/DTU21GRA, reach a magnitude greater than 5 cm. The differences among different MDTs result from the different data preprocessing techniques and weighting methods, as well as the different datasets used in model development. DTU21GRA has relatively small variations, especially when this profile crosses the Zhenghe Reefs (10.4°E) and Lesi Shoal (11.35°E). Moreover, the discrepancy between the XGM2019e_2159-derived MDT and the one strengthened by using the DTU21GRA data exceeds 3 cm (4 cm) over the Zhenghe Reefs (Lesi Shoal). Similarly, compared to the MDT modeled with the altimetric gravity model computed without SAR altimeter data, i.e., DTU13GRA, DTU15GRA, and DTU17GRA, the use of DTU21GRA data can reduce these spike-like errors by a magnitude exceeding 2 cm (3 cm) over the Zhenghe Reefs (Lesi Shoal).

Performances of Heterogeneous Altimetric Gravity Models
Moreover, we investigate the performances of heterogeneous gravity models released from different institutes in marine gravity field modeling. To achieve this, five representative models, i.e., HY2A, SCSGA, SIO30, SIO31, and DTU21GRA, are selected for the case study. The magnitudes of the added signals in terms of quasi-geoid heights derived from different gravity models are not consistent (see Figure 5). The maximum/minimum of the additional signals retrieved from HY2A is -5.5/7.0 cm (see Table  3). In comparison, the signals derived from SCSGA have weaker patterns, and the maximum/minimum value is -4.2/4.4 cm, whereas the added signals computed from the altimetric gravity data modeled with the Sentinel-3A/B SAR altimeter data, i.e., SIO30/SIO31/DTU21GRA, reach a magnitude greater than 5 cm. The differences among different MDTs result from the different data preprocessing techniques and weighting methods, as well as the different datasets used in model development.   The heterogeneous data preprocessing strategies, as well as data combination methods used in model development, partially cause these discrepancies among different MDTs. However, the most important reason may be due to that the computation of SIO30, SIO31, and DTU21GRA involved the Sentinel-3A/B SAR altimeter data over the study area. These results agree with the validation of heterogeneous gravity models against independent gravimetric observations over the SCS, which showed that HY2A had degraded quality versus SIO V27.1/DTU17GRA [42].    4°N, 115.7°E). The heterogeneous data preprocessing strategies, as well as data combination methods used in model development, partially cause these discrepancies among different MDTs. However, the most important reason may be due to that the computation of SIO30, SIO31, and DTU21GRA involved the Sentinel-3A/B SAR altimeter data over the study area. These results agree with the validation of heterogeneous gravity models against independent gravimetric observations over the SCS, which showed that HY2A had degraded quality versus SIO V27.1/DTU17GRA [42].  Although the MDT derived from the quasi-geoid enhanced by HY2A/SCSGA has smaller turbulences than that derived from XGM2019e_2159, the spike-like errors are still dramatic close to the Zhenghe Reefs (114.6 • E) and Wufang Reef (115.7 • E). The spikes are reduced in the solutions strengthened by using SIO31/DTU21GRA, with a reduction greater than 2 cm over the Zhenghe Reefs and 3 cm over the Wufang Reef, compared to the MDT modeled with the HY2A data. It agrees well with the modeling results in Figure 6, where the MDT modeled with the SIO31/DTU21GRA data has relatively smooth structures close to islands. The MDTs computed with the SIO30, SIO31, and DTU21GRA data almost show consistent structures along this profile, and the mutual discrepancies among these three solutions are within 2 cm. The properties of different MDTs along longitude 114.6 • N show similar results to the profile along latitude 10.4 • E (see Figure 7b). The MDTs from the altimetric gravity data computed with the SAR altimetry data, i.e., SIO30, SIO31, and DTU21GRA, have relatively small variations compared with the one derived from XGM2019e_2159 alone and the solution derived from the HY2A/SCSGA data. The application of the SAR altimetry data reduces the spike-like errors by a magnitude larger than 2 cm (4 cm) over the Zhenghe Reefs (Lesi Shoal), compared with the MDT modeled with the HY2A data. These results show that the performances of heterogeneous gravity models are not consistent, and the SAR-based models may be preferable.

Max
Remote Sens. 2022, 14, x FOR PEER REVIEW 11 of 19 pressed in the Plate Carrée projection. Note: only SIO30, SIO31, and DTU21GRA were computed with the SAR altimetry data derived from Sentinel-3A/B over the study area. Figure 7a shows the different MDTs along latitude 10.4°N. Although the MDT derived from the quasi-geoid enhanced by HY2A/SCSGA has smaller turbulences than that derived from XGM2019e_2159, the spike-like errors are still dramatic close to the Zhenghe Reefs (114.6°E) and Wufang Reef (115.7°E). The spikes are reduced in the solutions strengthened by using SIO31/DTU21GRA, with a reduction greater than 2 cm over the Zhenghe Reefs and 3 cm over the Wufang Reef, compared to the MDT modeled with the HY2A data. It agrees well with the modeling results in Figure 6, where the MDT modeled with the SIO31/DTU21GRA data has relatively smooth structures close to islands. The MDTs computed with the SIO30, SIO31, and DTU21GRA data almost show consistent structures along this profile, and the mutual discrepancies among these three solutions are within 2 cm. The properties of different MDTs along longitude 114.6° N show similar results to the profile along latitude 10.4°E (see Figure 7b). The MDTs from the altimetric gravity data computed with the SAR altimetry data, i.e., SIO30, SIO31, and DTU21GRA, have relatively small variations compared with the one derived from XGM2019e_2159 alone and the solution derived from the HY2A/SCSGA data. The application of the SAR altimetry data reduces the spike-like errors by a magnitude larger than 2 cm (4 cm) over the Zhenghe Reefs (Lesi Shoal), compared with the MDT modeled with the HY2A data. These results show that the performances of heterogeneous gravity models are not consistent, and the SAR-based models may be preferable.

The MDTs Modeled with Different Reference Models
In addition, we research the performances of various satellite altimeter-derived gravity models in marine gravity field modeling based on heterogeneous reference models, where three combined models that have similar expansion degrees to XGM2019e_2159, i.e., EGM2008, SGG-UGM-1, and GECO, are introduced. EGM2008 (d/o 2190/2159) was computed by merging GRACE observations with ground-based gravity data [9]. Moreover, two models computed by combining GOCE gravity gradients, namely GECO (d/o 2190/2159) [11] and SGG-UGM-1 (d/o 2159/2159) [59], are used. The application of GOCE gradients enhanced the global gravity field at low-frequency bands, approximately from degree 30 to 220 [60].
Three representative altimetric gravity models, i.e., HY2A, SIO31, and DTU21GRA, are used, where HY2A represents the altimetric gravity model computed by only using the conventional radar altimeter data, and SIO31 and DTU21GRA represent the SAR-based altimetric gravity models. Figure 8 shows the MDTs modeled from different reference models; for all the MDTs directly modeled from the GGMs alone, we find dramatic variations over areas close to islands. The associated solution augmented by using HY2A demonstrates smaller turbulences in comparison with the MDTs directly derived from these GGMs. This is mainly due to the incorporation of heterogeneous satellite altimeter-derived gravity data in model development, where EGM2008/SGG-UGM-1/GECO was computed by involving DNSC07GRA data at short-wavelength over the SCS. However, limited highquality altimeter data was used in developing DNSC07GRA, where no SARAL/AltiKa or CryoSat-2 data was included. In comparison, five years of CryoSat-2 data, two years of SARAL/AltiKa, and data from the HY-2A geodetic mission were included in the computation of the HY2A model, and, consequently, the HY2A model had improved accuracy versus DTU10GRA/DNSC07GRA [42]. The MDTs computed from the SIO31/DTU21GRA data further reduce these disturbances over island areas, by a magnitude of several centimeters, compared with the solutions modeled from the HY2A data. Moreover, the inconsistencies among the MDTs modeled from different reference models are prominent over the northern and eastern study area, owing to the heterogeneous data sources and modeling approaches adopted in those reference models' computation. Figures 9 and 10 show the MDTs computed from different reference models along latitude 10.4 • N and longitude 114.6 • E, respectively. For all the MDTs, we see that the GGMderived MDT has the largest variations, whereas the MDTs modeled additionally with the HY2A data reduces these disturbances: see the Zhenghe Reefs and Wufang Reef along latitude 10.4 • N and the Zhenghe Reefs and Lesi Shoal along longitude 114.6 • E. Moreover, the MDTs modeled from the SIO31/DTU21GRA data computed with the SAR altimetry data have less oscillations than those computed from the HY2A data. This is especially a case close to the Lesi Shoal along longitude 114.6 • E, where the use of DTU21GRA data reduces the spikes significantly, by a magnitude larger than 10 cm and 4 cm, compared to the results computed from these GGMs alone and ones modeled by additionally combining the HY2A data, respectively. This also demonstrates that the use of altimetric gravity data computed with the SAR altimetry data can significantly enhance the marine gravity field, no matter which global geopotential model is applied for modeling. Remote Sens. 2022, 14, x FOR PEER REVIEW 13 of 19

Conclusions
In this study, the satellite altimeter-derived gravity models are used in marine gravity field augmentation over island areas; in particular, we study the feasibility of regional enhancement by using SAR-based gravity models. The gravity field modeling results in the Spratly Islands demonstrate that the incorporation of the newly published

Conclusions
In this study, the satellite altimeter-derived gravity models are used in marine gravity field augmentation over island areas; in particular, we study the feasibility of regional enhancement by using SAR-based gravity models. The gravity field modeling results in the Spratly Islands demonstrate that the incorporation of the newly published satellite altimeter-derived gravity observations strengthens the marine gravity field. By involving the gravity data computed with the Sentinel-3A/B SAR altimetry data, the accuracies of regional solutions (including quasi-geoid and MDT models) are dramatically improved, by a magnitude larger than 4 cm around areas close to islands, in comparison with the results directly derived from a newly released high-degree GGM, namely, XGM2019e_2159. This suggests that the newly published altimetric gravity models contain additional signals that were unresolved in the currently available GGMs. Moreover, we model local gravity field based on heterogeneous combined GGMs. For all the solutions derived from the GGMs alone, significant variations are observed, which may bring about errors that are greater than 10 cm. These results show that the application of a GGM alone is inadequate in MDT modeling with centimeter-level accuracy over the areas close to land/island.
Further comparison of the solutions modeled from heterogeneous gravity models show that the performances of different models are heterogeneous, where the quasi-geoid/MDT solutions augmented with the SAR-based gravity models have better results. The use of the SAR-based gravity models in local augmentation improves the quasi-geoid/MDT solutions, by a magnitude of 2~4 cm, versus the solutions computed from the gravity models developed without SAR altimetry data. These results highlight the superiority of using the SAR-based altimetric gravity data in marine gravity field recovery, especially over the regions close to land/island. Moreover, the use of SAR altimetry data can alleviate the well-known coast problem, which also allays the issue of data scarcity in narrow/shallow waters. It is noticeable that the use of SAR-based altimetry data cannot fully overcome the coast problem; moreover, it also cannot replace the role of ground-based gravity data in gravity field modeling. By combing the ground-based gravity data, the local gravity field can be further improved.
Future work involves the further improvement of local gravity field model and the application of the model computed in this study in research fields such as geodesy and oceanography. By merging SAR altimeter records from the missions such as Sentinel-6 and using the improved waveform retracking methods, the local gravity field may be improved. Moreover, the quasi-geoid and mean dynamic topography models computed with the SAR altimetric gravity data are beneficial in understanding sea level change, ocean currents, and water exchanges with the surrounding regions close to the Spratly Islands. More importantly, these models are useful for establishing a unified vertical height datum over the Spratly Islands, where extensive coastlines and many islands exist, which make the traditionally used methods, such as the hydrostatic leveling, oceanic dynamic leveling, and trigonometric leveling methods, difficult to use for height datum unification.  Data Availability Statement: The satellite altimetric gravity models and mean sea surface models that developed at DTU space are available on https://ftp.space.dtu.dk/pub, accessed on 28 November 2021, and the altimetric gravity models computed by the Scripps Institution of Oceanography are freely accessible at https://topex.ucsd.edu/pub/global_grav_1min, accessed on 28 November 2021.
All the global geopotential models can be publicly accessed from http://icgem.gfz-potsdam.de/home, accessed on 28 November 2021.