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 a unified vertical datum system is to determine the vertical offset between the local vertical datum and the global vertical datum. For this purpose, the geodetic boundary value problem (GBVP) approach based on the remove-compute-restore (RCR) technique is used to determine the vertical datum parameters in this paper. In the RCR technique, a global geopotential model (GGM) is required to remove and restore the long wavelengths of the gravity field. The satellite missions of the GRACE (Gravity Recovery and Climate Experiment) and GOCE (Gravity field and steady-state Ocean Circulation Exploration) offer high accuracy medium–long gravity filed information, but GRACE/GOCE-based GGMs are restricted to medium–long wavelengths because the maximum degree of their spherical harmonic representation is limited, which is known as an omission error. To compensate for the omission error of GRACE/GOCE-based GGM, a weighting method is used to determine the combined GGM by combining the high-resolution EGM2008 model (Earth Gravitational Model 2008) and GRACE/GOCE-based GGM to effectively bridge the spectral gap between satellite and terrestrial data. An additional consideration for the high-frequency gravity signals is induced by the topography, and the residual terrain model (RTM) is used to recover the omission errors effect of the combined GGM. In addition, to facilitate practical implementation of the GBVP approach, the effects of the indirect bias term, the spectral accuracy of the GGM, and the systematic levelling errors and distortions in estimations of the vertical datum parameters are investigated in this study. Finally, as a result of the GBVP solution based on the combined DIR_R6/EGM2008 model, RTM, and residual gravity, the geopotential values of the North American Vertical Datum of 1988 (NAVD88), the Australian Height Datum (AHD), and the Hong Kong Principal Datum (HKPD) are estimated to be equal to 62636861.31 ± 0.96, 62653852.60 ± 0.95 and 62636860.55 ± 0.29 m²s⁻², respectively. The vertical offsets of NAVD88, AHD, and HKPD with respect to the global geoid are estimated as −0.809 ± 0.090, 0.082 ± 0.093, and −0.731 ± 0.030 m, respectively.

This study presents a methodology for constructing a quasi-geoid model with millimeter-level accuracy over the Shangyu area in China, following the guidelines of the International Association of Geodesy Joint Working Group 2.2.2, known as “The 1 cm geoid experiment”. Our approach combines two steps to ensure exceptional accuracy. First, we employ Molodensky’s theory to model the gravity field, accounting for non-level surfaces and considering complex terrain effects. Through an exhaustive analysis of these influential factors, we implement a comprehensive suite of applicable formulae within Molodensky’s series solution, enabling a thorough assessment of their impacts on height anomalies within the gravimetric quasi-geoid model. Second, we utilize a hybrid method that involves a multi-surface function using the least-squares method and a robust estimation technique. This approach enables the interpolation of quasi-geoid heights by incorporating ellipsoidal and leveling normal heights, as well as gravimetric quasi-geoid data. Through a numerical example, we demonstrate the efficiency of our solution concept, achieving an accuracy of 0.79 cm compared to independent global navigation satellite system (GNSS)/leveling measurements. By developing this methodology, our study contributes to the advancement of geodesy research and provides a valuable methodology for creating highly precise quasi-geoid models in geodetic applications.

The accurate estimation of the zero-height geopotential level in a local vertical datum (LVD) is critical for linking traditional height reference systems to a global height system. In this paper, we investigate the theoretical and practical challenges involved in determining the offset between the North American vertical datum (NAVD) and the global vertical datum (GVD). Drawing on the classical theory of the vertical system in physical geodesy, we define the vertical datum offset and derive rigorous formulas for its calculation. We examine various factors that affect the determination of the offset, including the global gravitational models (GGMs), geodetic reference system, tide system, tilt error, and omission error. Using terrestrial gravity data and gravity anomalies from multiple GGMs in conjunction with Global Navigation Satellite System (GNSS) and orthometric heights, we estimate the vertical offset between the NAVD and GVD. Our results indicate that the geopotential difference approach and the geodetic boundary value problem (GBVP) approach yield consistent results. When the normal gravity geopotential of the geodetic reference system is selected as the gravity geopotential of the global height datum, the NAVD is approximately 0.04 m higher than the GVD relative to the GRS80 ellipsoid, and 0.97 cm higher than the GVD relative to the WGS84 ellipsoid. When the Gauss–Listing geopotential value is chosen as the gravity geopotential of the global height datum, the NAVD is roughly 1.45 m higher than the GVD relative to the GRS80 ellipsoid, and approximately 0.52 m higher than the GVD relative to the WGS84 ellipsoid.

The Moho interface is an important parameter for describing the structure and morphology of the Earth's crust, and it is of significant importance in the study of the formation and evolution of the crust and mantle, as well as deep-seated dynamic processes (Stern et al., 2018). Existing Moho models derived from seismic data often suffer from inaccuracies due to irregular distribution and regional imbalances of seismic data. However, with the development of gravity satellite technology, high-precision satellite gravity data has injected new vitality into the study of lithospheric tectonic features and crustal evolution. In this study, constrained by seismic data (Li et al., 2013; Zhang et al., 2021), we utilized an improved regularized Bott method (Uieda et al., 2017) to invert high-precision satellite gravity data and obtained a high-precision unified Moho depth model for the East Asian region, encompassing both land and sea areas. The research results show that the Moho depth model exhibits a continuous increase in depth from east to west, and its overall distribution in the horizontal direction is non-uniform, displaying distinct regional block features. This paper provides a high-resolution and high-precision Moho model for studying the evolution of the East Asian continental tectonics and plate interactions, and further discusses the macrostructural framework and geological implications of East Asia.ReferencesLi Y Gao M, Wu Q. Crustal Thickness Map of the Chinese Mainland from Teleseismic Receiver Functions [J]. Tectonophysics, 2013, 611.Stern, Robert, J, et al. Continental crust of China: A brief guide for the perplexed [J]. Earth Science Reviews the International Geological Journal Bridging the Gap Between Research Articles & Textbooks, 2018.Uieda L, Barbosa V. Fast nonlinear gravity inversion in spherical coordinates with application to the South American Moho [J]. Geophysical Journal International, 2016.Zhang J , Yang G , Tan H , et al. Mapping the Moho depth and ocean-continent transition in the South China Sea using gravity inversion [J]. Journal of Asian Earth Sciences, 2021, 218(3–4):104864.

China is located in the southeast of the Eurasian Plate and is subject to the effects of subducting, squeezing and collision by the Pacific Plate to the east, Philippine Plate to the southeast and Indian Ocean Plate to the southwest. It has exceptional geotectonic structure. Based on the satellite gravity data with high precision, high resolution and ample geophysical information, combined with geological data, by using satellite gravity potential field and its full tensor gradient, this paper studies the distribution characteristics of gravity anomalies and the identification of tectonic boundaries in China and surrounding regions. Results suggest that the Bouguer gravity anomaly in eastern China reduces gradually from east to west, mostly in the direction of NNE; in the western, it reduces gradually in a wave mode from north to south, mainly in the directions of NW and NWW. In general, the stress field reduces gradually from west to east, and the tectonic of stress field in western China is complex. The maximum principal compressive stress in Xinjiang exists in SN direction and that in Qinghai-Tibet Plateau mostly changes gradually from NNE to SSE; the change in eastern China is relatively simple, and the maximum principal compressive stress direction gradually changes from NE to WE and then to SE. In addition to the above study results, by comprehensively referencing the previous studies by other people and by using the boundary identification methods based on the satellite gravity full-tensor gradient data and its combinations, we update the extension route of Red River fault zone and deduce the tectonic unit boundary between the North China and South China active tectonic block regions. This paper identifies in China and surrounding regions 6 primary active tectonic blocks, 22 secondary active tectonic blocks, 3 tertiary active tectonic blocks and the 20 active tectonic block boundary zones constituted of deformation belts and active tectonic belts with various geometric structures and width variations. The results of this study can improve the understanding of gravity anomalies and boundary structures in China and surrounding regions, and provide certain geophysical supports for geological structure analysis and crustal dynamic process.