Integration of a Priori Gravity Field Models in Boundary Element Formulations to Geodetic Boundary Value Problems

Abstract

Current high resolution geopotential models of the Earth are based on a combination of satellite and terrestrial data. Satellite data are well-suited to recover the long-wavelength features of the geopotential up to some degree N, whereas terrestrial gravity and height data fix the medium and short wavelengths. Usually, the recovering of the medium and short-wavelengths from terrestrial data is formulated as a boundary value problem (BVP) for the difference between the Earth’s geopotential and the long-wavelength geopotential model as derived from satellite data commonly referred to as the disturbance potential. Since a number of geopotential coefficients of the satellite model cannot be improved by terrestrial data, we should fix them when solving the BVP. Then we are faced with a constrained (overdetermined) BVP for the Laplace equation. This has implications for the representation formula and/or the choice of the trial & test space in Galerkin boundary element methods.We consider multipole representation, modified kernel functions, and modified trial spaces. The latter are the best choice for theoretical and numerical reasons. We propose a general method to construct a system of base functions that fix an a priori given set of geopotential coefficients. In addition, we address the problem of compression rates and stability, which implies the use of multiscale base functions. Various implementations are tested and compared for the altimetry-gravimetry II BVP.KeywordsBoundary Element MethodGeopotential ModelTest SpaceTerrestrial DataGeopotential CoefficientThese keywords were added by machine and not by the authors. This process is experimental and the keywords may be updated as the learning algorithm improves.