Abstract:Gravimetric measurements allow monitoring mass transport in the Earth system and one example is temporal variations in the global water cycle including the cryosphere. Coupling gravity observations to the underlying processes requires modeling of the relevant mass elements. When observing the gravity field from space (i.e., by satellites), point mass models or simple spherical surface layers (mascons) are a reasonable choice. However, for terrestrial observations taken right on the surface of the relevant masses, precise modeling of the geometry of those mass elements and the related gravitational effects becomes crucial for the quality of the modeling and the correct interpretation of the results. In this study, two different modeling techniques are investigated, the right rectangular prism and the general polyhedron. Both approximations provide the possibility to analytically calculate the corresponding gravitational signal. Although modeling in terms of polyhedrons is numerically more complex and less efficient, the representation of the actual surface is superior as the inclined polyhedral surface fits much better to the actual topographic surface than the step function implied by rectangular prisms with horizontal top plane. The case study presented here is related to the gravimetric mass balance of a mountain glacier (Vernagtferner in the Ötz valley, Austria) where on the one hand high resolution digital terrain models (DTM of resolution 1 m × 1 m) are available, while on the other hand, the accuracy requirements are very high (on the level of 10 −7 m/s 2 , i.e., 10 μGal). We evaluate first and second order derivatives of the gravitational potential induced by the melting of ice masses between epochs 2009 and 2016 at gravity sites on and off the actual glacier body. Furthermore, the relation of ground surface slope to the differences between the modeling approaches is investigated for a small part of the DTM, containing a variety of different slope regimes. Besides the original grid spacing, denser and coarser grids are also evaluated in order to examine the relation of the considered DTM resolution to the modeling differences. The modeling choice, the location of the computation points and the morphology of the mass distribution are the main factors that determine the quality of the results. The comparison of the two methods for the original 1 m resolution, provides differences up to 2.9 μGal for the vertical component of the attraction V z , 9.1 μGal/m for the gradient V zx , 4.6 μGal/m for V yz and 17.1 μGal/m for V zz . Coarser grid resolutions increase these values significantly, up to 500 μGal at specific locations, while grids denser than the 1 m resolution decrease them for about 48%.

Gravity field modeling in mountainous areas based on band-limited SRBFs

Regional Gravity Field Modeling Using Band-Limited SRBFs: A Case Study in Colorado

High-resolution regional gravity field modeling in data-challenging regions for the realization of geopotential-based height systems

SRBF_Soft: a Python-based open-source software for regional gravity field modeling using spherical radial basis functions based on the data-adaptive network design methodology

On the spectral combination of satellite gravity model, terrestrial and airborne gravity data for local gravimetric geoid computation

Experiences with the use of mass-density maps in residual gravity forward modelling

Spherical radial basis functions model: approximating an integral functional of an isotropic Gaussian random field

A High-Resolution Earth’s Gravity Field Model SGG-UGM-2 from GOCE, GRACE, Satellite Altimetry, and EGM2008

GNSS Gravity Leveling

Quasi Geoid and Geoid Modeling with the Use of Terrestrial and Airborne Gravity Data by the GGI Method—A Case Study in the Mountainous Area of Colorado

Augmented Gravity Field Modelling by Combining EIGEN_6C4 and Topographic Potential Models

Forward Gravity Modelling to Augment High-Resolution Combined Gravity Field Models

The Refined Gravity Field Models for Height System Unification in China

Analysis for Satellite Gravity Field Measurement Based on Relative Weight of Gravitational Coefficients

A hybrid approach for recovering high-resolution temporal gravity fields from satellite laser ranging

On the combined gravity gradient modeling for applied geophysics

Enhancing Regional Quasi-Geoid Refinement Precision: An Analytical Approach Employing ADS80 Tri-Linear Array Stereoscopic Imagery and GNSS Gravity-Potential Leveling

Height Variations in the Mississippi River Basin From Daily Time-Varying Satellite Gravity Data and Hydrological Model

RTM Gravity Forward-Modeling Using Topography/Bathymetry Data to Improve High-Degree Global Geopotential Models in the Coastal Zone

Integration of satellite geodetic observations for regional geoid modeling using remove-compute-restore technique

Analytical computation of local gravitational effects of mountain glacier mass change from polyhedral and prismatic modeling - test case Vernagtferner, Austrian Alps