Conditionality of inverse solutions to discretised integral equations in geoid modelling from local gravity data

Robert Tenzer, P. Novák

Research output: Journal article publicationJournal articleAcademic researchpeer-review

8 Citations (Scopus)

Abstract

The eigenvalue decomposition technique is used for analysis of conditionality of two alternative solutions for a determination of the geoid from local gravity data. The first solution is based on the standard two-step approach utilising the inverse of the Abel-Poisson integral equation (downward continuation) and consequently the Stokes/ Hotine integration (gravity inversion). The second solution is based on a single integral that combines the downward continuation and the gravity inversion in one integral equation. Extreme eigenvalues and corresponding condition numbers of matrix operators are investigated to compare the stability of inverse problems of the above-mentioned computational models. To preserve a dominantly diagonal structure of the matrices for inverse solutions, the horizontal positions of the parameterised solution on the geoid and of data points are identical. The numerical experiments using real data reveal that the direct gravity inversion is numerically more stable than the downward continuation procedure in the two-step approach.
Original languageEnglish
Pages (from-to)53-70
Number of pages18
JournalStudia Geophysica et Geodaetica
Volume52
Issue number1
DOIs
Publication statusPublished - 1 Jan 2008
Externally publishedYes

Keywords

  • Geoid
  • Ggravity
  • Inverse problem
  • Numerical stability

ASJC Scopus subject areas

  • Geophysics
  • Geochemistry and Petrology

Cite this