Modification methods of the Stokes’ kernel for determining the (quasi-) geoid with the Remove-Compute-Restore technique

The geoid and quasi-geoid serve as the reference surfaces of the orthometric and normal height systems, respectively. In order to improve the accuracy of the (quasi-) geoid determined by the Stokes integral with use of the Remove-Compute-Restore (RCR) technique, various modification methods for the spherical Stokes’ kernels, including the spheroidal, cosine-, power-, and Molodensky-modified kernels, are studied in this paper. In addition to the traditional Molodensky-modified Stokes’ kernel, a more effective Molodensky-modified Stokes’ kernel is put forward. A general formula for spectral decomposition of the Stokes integral in the RCR mode is derived, followed by the spectral analysis to reveal the transfer principles of gravity data when using different Stokes’ kernels. The spheroidal and modified Stokes integrals can cause spectral leakage phenomenon, and a method to eliminate spectral leakage is presented based on spectral analysis. The research indicates the low truncation degree of the spheroidal Stokes’ kernel and the low modification degrees of the modified Stokes’ kernel affect the accuracy of the (quasi-) geoid significantly. Quantitative methods for estimating the empirical values of the parameters of the low-degree spheroidal and modified Stokes’ kernels are proposed and the effectiveness of the methods is validated through numerical tests.