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.

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.

The coordinates of the KBR (K-band ranging system) antenna phase center of GRACE-type gravity satellites in the satellite Science Reference Frame should be precisely known, and the determination accuracy should reach 0.3 mrad in the Y (pitch) and Z (yaw) directions. Due to the precision limitation of ground measurement and the change of space environment during orbit, the KBR antenna phase center changes. In order to obtain more accurate KBR antenna phase center coordinates, it is necessary to maneuver the satellite to achieve the on-orbit calibration of the KBR antenna phase center. Based on the in-orbit calibration data of KBR of GRACE-FO satellites, a new method is proposed to estimate the antenna phase center of KBR using the inter-satellite range acceleration as the observation value. The antenna phase center of KBR is solved by the robust estimation method, and the obtained calibration results are better than 72 μm in the Y and Z directions and better than 1.3 mm in the X direction, which is 50% better than the least squares estimation algorithm. The accuracy of KBR calibration results obtained by using the data of positive maneuvers or mirror (negative) maneuvers, respectively, does not meet 0.3 mrad. It is shown that mirror maneuvers are required for KBR calibration of a GRACE-type gravity satellite to obtain antenna phase center estimation results that meet the accuracy requirements. The calibration algorithm proposed in this paper can provide reference for KBR antenna phase center calibration of Chinese GRACE-type gravity satellites.

Unification of the global vertical datum has been a key problem to be solved for geodesy over a long period, and the main challenge for unifying the global vertical datum system is to determine the geopotential value W₀ of the geoid and to calculate the vertical offset between the local vertical datum and the global vertical datum W₀. The geopotential value W₀ can be calculated using the grid mean sea surface (GMSS) data and the global geopotential model (GGM). In this study, this GMSS data was measured with adjustment methods and 24 years of merged multi-satellite altimetry data. The data of HaiYang-2A (HY-2A) and Jason-3 were first used to calculate W₀. The geopotential value W₀ was determined to be 62,636,856.82 m²s⁻² by combining the EIGEN-6C4 (European Improved Gravity Model of the Earth by New Techniques) and the GMSS data. Then, the geopotential difference approach and geodetic boundary value problem (GBVP) approach were used to determine the vertical datum parameters in this study. To compensate for the omission error of the GGM, this study utilized the remove–compute–restore (RCR) technique and the residual terrain model (RTM)-recovered high-frequency gravity signals. Finally, as a result of the GBVP solution, the geopotential value of the Australian Height Datum (AHD) was 62,636,851.935 m²s⁻², and the vertical offset of the AHD relative to the global vertical datum W₀ was 0.4885 m. As a result of the geopotential difference approach, the geopotential value of the Chinese Height datum was 62636861.412 m²s⁻², and the vertical offset of the Chinese Height datum was −0.4592 m.