Uncertainties associated with integral-based solutions to geodetic boundary-value problems

Physical geodesy applies potential theory to study the Earth’s gravitational field in space outside and up to a few km inside the Earth’s mass. Among various tools offered by this theory, boundary-value problems are particularly popular for the transformation or continuation of gravitational field parameters across space. Traditional problems, formulated and solved as early as in the nineteenth century, have been gradually supplemented with new problems, as new observational methods and data are available. In most cases, the emphasis is on formulating a functional relationship involving two functions in 3-D space; the values of one function are searched but unobservable; the values of the other function are observable but with errors. Such mathematical models (observation equations) are referred to as deterministic. Since observed data burdened with observational errors are used for their solutions, the relevant stochastic models must be formulated to provide uncertainties of the estimated parameters against which their quality can be evaluated. This article discusses the boundary-value problems of potential theory formulated for gravitational data currently or in the foreseeable future used by physical geodesy. Their solutions in the form of integral formulas and integral equations are reviewed, practical estimators applicable to numerical solutions of the deterministic models are formulated, and their related stochastic models are introduced. Deterministic and stochastic models represent a complete solution to problems in physical geodesy providing estimates of unknown parameters and their error variances (mean squared errors). On the other hand, analyses of error covariances can reveal problems related to the observed data and/or the design of the mathematical models. Numerical experiments demonstrate the applicability of stochastic models in practice.