Global Gravity Field Modelling by Solving the Infinite Nonlinear Fixed Geodetic Boundary Value Problem

Abstract The aim of presented paper is to solve the nonlinear geodetic boundary value problem (BVP) by the finite element method (FEM) involving the mapped infinite elements (MIE). In comparison to our previous works, see [ M acák , M. et al .: On an iterative approach to solving the nonlinear satellite-fixed geodetic boundary-value problem . In: IAG Symp. Vol. 142 (2016), pp. 185–192.] and [M acák , M. et al .: A gravity field modelling in mountainous areas by solving the nonlinear satellite-fixed geodetic boundary value problem with the finite element method , Acta Geodaetica et Geophysica, 58 (2023), 305–320.] dealing with bounded domains, in this paper we propose and study numerical concept on unbounded domains, given as an exterior BVP for the Laplace equation outside the gravitating body, e.g. Earth, with the nonlinear boundary condition (BC) prescribed on the Earth’s surface and considering the solution regularity condition at infinity. This concept can be found in many scientific disciplines being also the most natural from physical geodesy point of view, see, e.g., [B ackus , G. E.: Application of a non-linear boundary-value problem for Laplace ’ s equation to gravity and geomagnetic intensity surveys , Q. J. Mech. Appl. Math. 2 (1968), 195–221.] and of large practical importance when we are not able to prescribe BCs on a bounded domain. The proposed concept is based on the iterative procedure, and as the numerical method we have implemented the FEM with the MIE to take into account the regularity of the disturbing potential at infinity. Since the boundary of the computational domain is the discretized real Earth’s surface considering its topography, as finite and infinite elements we have chosen the triangular prisms. We study and verify this numerical approach by a testing experiment with a homogeneous sphere, by the experiment using EGM2008, and finally, we present one detailed numerical experiment with DTU21GRA data.

Authors

Institutions

Publication Details

Journal
Tatra Mountains Mathematical Publications
Published
2025-09-12
DOI
https://doi.org/10.2478/tmmp-2025-0016
Primary Topic
Geophysics and Gravity Measurements
Type
article
Field-Weighted Citation Impact
0.00
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

Global Gravity Field Modelling by Solving the Infinite Nonlinear Fixed Geodetic Boundary Value Problem

Karol Mikula, Zuzana Minarechová, Róbert Čunderĺık, Marek Macák
Tatra Mountains Mathematical Publications
Geophysics and Gravity Measurements
article

Global Gravity Field Modelling by Solving the Infinite Nonlinear Fixed Geodetic Boundary Value Problem

Karol Mikula, Zuzana Minarechová, Róbert Čunderĺık, Marek Macák
article en

Abstract

Abstract The aim of presented paper is to solve the nonlinear geodetic boundary value problem (BVP) by the finite element method (FEM) involving the mapped infinite elements (MIE). In comparison to our previous works, see [ M acák , M. et al .: On an iterative approach to solving the nonlinear satellite-fixed geodetic boundary-value problem . In: IAG Symp. Vol. 142 (2016), pp. 185–192.] and [M acák , M. et al .: A gravity field modelling in mountainous areas by solving the nonlinear satellite-fixed geodetic boundary value problem with the finite element method , Acta Geodaetica et Geophysica, 58 (2023), 305–320.] dealing with bounded domains, in this paper we propose and study numerical concept on unbounded domains, given as an exterior BVP for the Laplace equation outside the gravitating body, e.g. Earth, with the nonlinear boundary condition (BC) prescribed on the Earth’s surface and considering the solution regularity condition at infinity. This concept can be found in many scientific disciplines being also the most natural from physical geodesy point of view, see, e.g., [B ackus , G. E.: Application of a non-linear boundary-value problem for Laplace ’ s equation to gravity and geomagnetic intensity surveys , Q. J. Mech. Appl. Math. 2 (1968), 195–221.] and of large practical importance when we are not able to prescribe BCs on a bounded domain. The proposed concept is based on the iterative procedure, and as the numerical method we have implemented the FEM with the MIE to take into account the regularity of the disturbing potential at infinity. Since the boundary of the computational domain is the discretized real Earth’s surface considering its topography, as finite and infinite elements we have chosen the triangular prisms. We study and verify this numerical approach by a testing experiment with a homogeneous sphere, by the experiment using EGM2008, and finally, we present one detailed numerical experiment with DTU21GRA data.

Tatra Mountains Mathematical PublicationsVol. 91(3)
Slovak University of Technology in Bratislava (SK)
Openalex Percentile: Top 74%
Geophysics and Gravity Measurements
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.