Abstract

Residual terrain modelling (RTM) plays a key role for short-scale gravity modelling in physical geodesy, e.g. for interpolation of observed gravity and augmentation of global geopotential models (GGMs). However, approximation errors encountered in RTM computation schemes are little investigated. The goal of the present paper is to examine widely used classical RTM techniques in order to provide insights into RTM-specific approximation errors and the resulting RTM accuracy. This is achieved by introducing a new, independent RTM technique as baseline that relies on the combination of (1) a full-scale global numerical integration in the spatial domain and (2) ultra-high-degree spectral forward modelling. The global integration provides the full gravity signal of the complete (detailed) topography, and the spectral modelling that of the RTM reference topography. As a main benefit, the RTM baseline technique inherently solves the “non-harmonicity problem” encountered in classical RTM techniques for points inside the reference topography. The new technique is utilized in a closed-loop type testing regime for in-depth examination of four variants of classical RTM techniques used in the literature which are all affected by one or two types of RTM-specific approximation errors. These are errors due to the (1) harmonic correction (HC) needed for points located inside the reference topography, (2) mass simplification, (3) vertical computation point inconsistency, and (4) neglect of terrain correction (TC) of the reference topography. For the Himalaya Mountains and the European Alps, and a degree-2160 reference topography, RTM approximation errors are quantified. As key finding, approximation errors associated with the standard HC ( $$ 4\pi G\rho H_{\text{P}}^{\text{RTM}} ) $$ may reach amplitudes of ~ 10 mGal for points located deep inside the reference topography. We further show that the popular RTM approximation ( $$ 2\pi G\rho H_{\text{P}}^{\text{RTM}} - {\text{TC}} $$ ) suffers from severe errors that may reach ~ 90 mGal amplitudes in rugged terrain. As a general conclusion, the RTM baseline technique allows inspecting present and future RTM techniques down to the sub-mGal level, thus improving our understanding of technique characteristics and errors. We expect the insights to be useful for future RTM applications, e.g. in geoid modelling using remove–compute–restore techniques, and in the development of new GGMs or high-resolution augmentations thereof.

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.