Forward modelling in the space domain is a very important task in geodesy and other geosciences. From topographical or isostatic information in the form of digital terrain model (DTM) and density model, the effects of these parameters or their derivatives on the gravity potential can be evaluated for different applications. In most cases, height or height-layer models are in use, which are gridded with respect to spherical coordinates. This holds for global as well as regional or even local applications. The definition of the spherical gridlines leads immediately to the spherical volume element, that is, the tesseroid. Only in the specific case that the observation point is located on the symmetry axis of the spherical coordinate system does the Newton integral have a closed analytical solution. More specifically, the effect of a tesseroid can be determined by evaluating the analytical solution of a segment of a spherical zonal band. To apply this aspect in practice, the DTM must be transformed into the local spherical azimuthal system of the observation point (UNIPOL approach). In the general case, the Newton integral can be solved, for example, using a Taylor series expansion of the integral kernel and a subsequently applied term-wise integration (GIK approach). Within this contribution, the two fundamentally different tesseroid approaches, namely, the GIK and the UNIPOL approach are compared. This comparison is performed, in particular, with regard to the required computational time and the approximation error under different test scenarios. The numerical studies show that both approaches are equivalent in terms of accuracy for both the gravitational potential and gravity; however, the UNIPOL approach is more time consuming because, for each observation point, the whole DTM must be transformed. Small numerical differences exist between the compared approaches for special constellations regarding the source point and the observation point.

Comparison of two different approaches for computing the gravitational effect of a tesseroid / A.M. Marotta, K. Seitz, R. Barzaghi, T. Grombein, B. Heck. - In: STUDIA GEOPHYSICA ET GEODAETICA. - ISSN 0039-3169. - 63:3(2019 Jul), pp. 321-344.

Comparison of two different approaches for computing the gravitational effect of a tesseroid

A.M. Marotta
Primo
;
2019-07

Abstract

Forward modelling in the space domain is a very important task in geodesy and other geosciences. From topographical or isostatic information in the form of digital terrain model (DTM) and density model, the effects of these parameters or their derivatives on the gravity potential can be evaluated for different applications. In most cases, height or height-layer models are in use, which are gridded with respect to spherical coordinates. This holds for global as well as regional or even local applications. The definition of the spherical gridlines leads immediately to the spherical volume element, that is, the tesseroid. Only in the specific case that the observation point is located on the symmetry axis of the spherical coordinate system does the Newton integral have a closed analytical solution. More specifically, the effect of a tesseroid can be determined by evaluating the analytical solution of a segment of a spherical zonal band. To apply this aspect in practice, the DTM must be transformed into the local spherical azimuthal system of the observation point (UNIPOL approach). In the general case, the Newton integral can be solved, for example, using a Taylor series expansion of the integral kernel and a subsequently applied term-wise integration (GIK approach). Within this contribution, the two fundamentally different tesseroid approaches, namely, the GIK and the UNIPOL approach are compared. This comparison is performed, in particular, with regard to the required computational time and the approximation error under different test scenarios. The numerical studies show that both approaches are equivalent in terms of accuracy for both the gravitational potential and gravity; however, the UNIPOL approach is more time consuming because, for each observation point, the whole DTM must be transformed. Small numerical differences exist between the compared approaches for special constellations regarding the source point and the observation point.
tesseroid; gravity forward modelling; Newton integral
Settore GEO/10 - Geofisica della Terra Solida
Settore ICAR/06 - Topografia e Cartografia
Article (author)
File in questo prodotto:
File Dimensione Formato  
10.1007_s11200-018-0454-2.pdf

accesso riservato

Tipologia: Publisher's version/PDF
Dimensione 2.26 MB
Formato Adobe PDF
2.26 MB Adobe PDF   Visualizza/Apri   Richiedi una copia
Pubblicazioni consigliate

Caricamento pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: http://hdl.handle.net/2434/662321
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 3
  • ???jsp.display-item.citation.isi??? 2
social impact