Topological inversion in geodesy-based, non-linear problems in geophysics

  • Authors:
  • Vasso Saltogianni;Stathis C. Stiros

  • Affiliations:
  • Laboratory of Geodesy and Geodetic Applications, Department of Civil Engineering, University of Patras, Patras 26504, Greece;Laboratory of Geodesy and Geodetic Applications, Department of Civil Engineering, University of Patras, Patras 26504, Greece

  • Venue:
  • Computers & Geosciences
  • Year:
  • 2013

Quantified Score

Hi-index 0.00

Visualization

Abstract

Geophysical phenomena such as volcanism and earthquake faulting are modeled on the basis of geodetic observations leading to redundant systems of highly non-linear equations. Still, such systems of equations cannot be solved using traditional analytical techniques. For this reason forward modeling or numerical inversion techniques are used, but these techniques have serious limitations. To overcome these problems, we propose a numerical/topological, grid-search based technique in the R^m space, a generalization and refinement of techniques used in some cases of low-accuracy 2-D positioning. In contrast to conventional solutions which tend to minimize a certain function and directly obtain a point solution of a system of equations, our algorithm has a different strategy. At first, an optimal R^m space which is defined by a grid and which contains the true solution is mapped as the intersection of grid spaces defined on the basis of the uncertainties of each observation. Then, the center of weight of this R^m space and its variance-covariance matrix are computed, and this point solution approximates the true solution of the system of equations. The efficiency of the proposed algorithm and its compatibility with conventional adjustment is demonstrated on the basis of two characteristic case studies.