Finite volume method: a good match to airborne gravimetry?

Abstract

Numerical methods, like the finite element method (FEM) or finite volume method (FVM), are widely used to provide solutions in many boundary value problems. In previous studies, these numerical methods have also been applied in geodesy but demanded extensive computations because the upper boundary condition was usually set up at the satellite orbit level, hundreds of kilometers above the Earth. The relatively large distances between the lower boundary of the Earth's surface and the upper boundary exacerbate the computation loads because of the required discretization in between. Considering that many areas, such as the US, have uniformly distributed airborne gravity data just a few kilometers above the topography, we adapt the upper boundary from the satellite orbit level to the mean flight level of the airborne gravimetry. The significant decrease in the domain of solution dramatically reduces the large computation demand for FEM or FVM. This paper demonstrates the advantages of using FVM in the decreased domain in simulated and actual field cases in study areas of interest. In the simulated case, the FVM numerical results show that precision improvement of about an order of magnitude can be obtained when moving the upper boundary from 250 to 10 km, the upper altitude of the GRAV-D flights. A 2–3 cm level of accurate quasi-geoid model can be obtained for the actual datasets depending on different schemes used to model the topographic mass. In flat areas, the FVM solution can reach to about 1 cm precision, which is comparable with the counterparts from classical methods. The paper also demonstrates how to find the upper boundary if no airborne data are available. Finally, the numerical method provides a 3D discrete representation of the entire local gravity field instead of a surface solution, a (quasi) geoid model.