Currently, seafloor topography inversion based on satellite altimetry gravity data provides the principal means to predict the global seafloor topography. Researchers often use sea surface geoid height or gravity anomaly to predict sea depth in the space domain. In this paper, a comprehensive discussion on seafloor topography inversion formulas in the space domain is presented using sea surface geoid height, gravity anomaly and introduces an approach that uses vertical gravity gradient. This would be the first study to estimate seafloor topography by vertical gravity gradient in the space domain. Further, a nonlinear iterative least-square inversion process is discussed. Using the search area for the Malaysia Airlines Flight MH370 as study site, we used the DTU17 gravity anomaly model and SIO V29.1 vertical gravity gradient to generate the seafloor topography. The results of the proposed bathymetric models were analyzed and compared with the DTU18 and SIO V20.1 bathymetric models. The experimental results show that the gravity anomaly and vertical gravity gradient in the study area are strongly correlated with the seafloor topography in the 20–200 km wavelength range. The optimal initial iteration values for seafloor topography variance and correlation length are 0.6365 ${\mathrm{km}}^2$ and 10.5′, respectively. Shipborne measurements from SONAR data were used as external checkpoints to evaluate the bathymetric models. The results show that the RMS for BAT_VGG_ILS (inversion model constructed by vertical gravity gradient) is smaller than for BAT_GA_ILS (inversion model constructed by gravity anomaly) and BAT_GA_VGG_ILS (inversion model constructed by gravity anomaly and vertical gravity gradient). The relative accuracy of the DTU18 bathymetry model was 9.27%, while the relative accuracy of the proposed seafloor models was higher than 4%. Within the 200 m difference range, the proportion of checkpoints for BAT_VGG_ILS was close to 95%, about 80% for BAT_GA_ILS and BAT_GA_VGG_ILS, and less than 50% for the DTU18. The results show that the nonlinear iterative least square method in the space domain is feasible.

We introduce an iterative inversion method to address the problems in high-order seafloor topography inversion using gravity data (gravity anomaly and vertical gravity gradient anomaly), such as the difficulty in computing the equation and the uniqueness of the calculation results. A part of the South China Sea is selected as the experimental area. Considering the coherence and admittance function of gravity topography and vertical gravity gradient topography, the inversion band of the gravity anomaly and vertical gravity gradient anomaly in the study area is 30 km–120 km. Seafloor topography models of different orders are constructed using an iterative method, and the performance of each seafloor topography model is analyzed against ETOPO1 and other seafloor topography models. The experimental results show that as the inversion order increases, the clarity and richness of seafloor topographic expression continuously improve. However, the accuracy of seafloor topography inversion does not improve significantly when the inversion order exceeds a certain value, which is related to the contribution of high-order seafloor topography to gravity information. The results show that the accuracy of BGT4 (inversion model constructed by the gravity anomaly) is slightly poorer than that of BVGGT4 (inversion model constructed by the vertical gravity gradient anomaly) in areas with complex topography, such as multi-seamounts and trenches, and the results are generally better in areas with flat seafloor topography.

One of the key constraints for the accelerometer of GRACE-type gravity satellites to accurately measure the non-gravitational accelerations acting on the satellite is that the center of mass of the satellite and the proof mass of the accelerometer should maintain a coincidence. In addition, the accuracy requirement is that the center of mass offset (CM-offset) in the three directions is less than 100 microns. Since the center of mass (CoM) of the satellite will change with the consumption of cold-gas fuel in the tanks, it is necessary to regularly carry out the CoM calibration maneuver. Firstly, the observation equations consisting of the accelerometer linear acceleration, angular acceleration, and the CM-offset vector are established in order to estimate the amount of CM-offset. Then, according to the estimated CM-offset, the satellite mass trim mechanisms are used to change the satellite’s CoM, so that the satellite’s CoM always approaches the proof mass of the accelerometer, with an accuracy of 100 μm per axis. The CM-offset of the satellite of GRACE-FO is estimated by using the accelerometer, star camera, magnetic torquer, magnetometer, and the precision orbit data during the GRACE-C CM-offset calibration period on 1 February 2020. Four kinds of CM-offset results are obtained by four different angular accelerations as follows: the angular acceleration based on the attitude dynamics (“MTQ angular acceleration”), the accelerometer angular acceleration calibrated by MTQ, the accelerometer angular acceleration, and the angular acceleration calculated by the star camera. By comparing the four kinds of CM-offset results that are estimated by the four different methods, all four of the results are shown to have the same level of accuracy. Based on the accelerometer (calibrated) angular acceleration, the difference with the JPL result is 0.5 μm, while the difference between the conventional method and the JPL result is 6.0 μm. All four of the methods can achieve the requirement of 50 μm accuracy and using four CM-offset estimation methods simultaneously can improve the integrity of the calibration results. Subsequently, the CM-offset results of GRACE-C since its launch are estimated here. The calibration algorithm that is proposed in this paper can be used as a reference in the calibration of gravity satellites carrying an accelerometer payload.

To address the limitations in global seafloor topography model construction, a scheme is proposed that takes into account the efficiency of seafloor topography prediction, the applicability of inversion methods, the heterogeneity of seafloor environments, and the inversion advantages of sea surface gravity field element. Using the South China Sea as a study area, we analyzed and developed the methodology in modeling the seafloor topography, and then evaluated the feasibility and effectiveness of the modeling strategy. Based on the proposed modeling approach, the STO_IEU2020 global bathymetry model was constructed using various input data, including the SIO V29.1 gravity anomaly (GA) and vertical gravity gradient anomaly (VGG), as well as bathymetric data from multiple sources (single beam, multi-beam, seismic, Electronic Navigation Chart, and radar sensor). Five evaluation areas located in the Atlantic and Indian Oceans were used to assess the performance of the generated model. The results showed that 79%, 89%, 72%, 92% and 93% of the checkpoints were within the ±100 m range for the five evaluation areas, and with average relative accuracy better than 6%. The generated STO_IEU2020 model correlates well with the SIO V20.1 model, indicating that the proposed construction strategy for global seafloor topography is feasible.

Compared with airborne gravimetry, a technique frequently used to infer the seafloor topography at places inaccessible to ship soundings due to the presence of ice shelf or ice mélange, airborne gravity gradiometry inherently could achieve higher spatial resolution, thus it is promising for improved inference of seafloor topography. However, its estimation capability has not been demonstrated by real projects. Theoretical analysis through admittance shows that compared with gravity disturbance, gravity gradient is more sensitive to the short-wavelength seafloor topography but diminishes faster with the increase of the distance between the seafloor and airplane, indicating its superiority is recovering short-wavelength topographic features over shallow waters. We present the first numerical experiment that estimates seafloor topography from a 0.4-km resolution, real airborne gravity gradients. It is shown that airborne gravity gradiometry can recover smaller topographic features than typical airborne gravimetry, but the estimation accuracy is only $\pm 17$ m due to the presence of subsurface density variations. The long-wavelength effect of the subsurface density variations can be removed with the aid of constraining bathymetry inside the study area, whereas the short wavelengths cannot. This study expands the applications of airborne gravity gradiometry, and helps glaciologists understand its performance in seafloor topography estimation.

A high-precision marine gravity gradient reference map is key to enabling underwater gravity gradient matching navigation. At present, the construction of the reference maps is based on quadrilateral geographic grids. However, quadrilateral grids lead to detriangulation at high latitudes, which limits the global applicability of such maps to underwater gravity navigation. To circumvent the limitations of quadrilateral grids, a hexagonal grid is introduced for constructing the reference map. This paper analyzes the characteristics of the icosahedral Snyder equal area aperture 4 hexagon (ISEA4H) and H3 grid systems and selects an appropriate grid system. In addition, we calculate and analyze the grid model errors and matching positioning errors of hexagonal and quadrilateral grids at the same resolution. The experimental results show that the grid model and matching positioning errors of a hexagonal grid system are more than 14% and 15% less than those of a quadrilateral grid system, respectively, indicating the feasibility and effectiveness of applying hexagonal grids to gravity gradient matching navigation. Given the low construction efficiency of a marine hexagonal grid gravity gradient reference map, we propose an efficient CPU+GPU hybrid parallel scheme. A global total tensor marine hexagonal grid gravity gradient reference map model is then constructed.