A greedy algorithm for nonlinear inverse problems with an application to nonlinear inverse gravimetry

• Published in: GEM international journal on geomathematics 2018-11, Vol.9 (2), p.167-198
• Main Authors: Kontak, Max, Michel, Volker
• Format: Article
• Language: English
• Subjects: Applications of Mathematics; Computational Science and Engineering; Earth Sciences; Mathematical Problems in Medical Imaging and Earth Sciences; Mathematics; Mathematics and Statistics; Original Paper
• ISSN: 1869-2672; 1869-2680
• DOI: 10.1007/s13137-018-0110-6

Based on the Regularized Functional Matching Pursuit (RFMP) algorithm for linear inverse problems, we present an analogous iterative greedy algorithm for nonlinear inverse problems, called RFMP_NL. In comparison to established methods for nonlinear inverse problems, the algorithm is able to combine very diverse types of basis functions, for example, localized and global functions. This is important, in particular, in geoscientific applications, where global structures have to be distinguished from local anomalies. Furthermore, in contrast to other methods, the algorithm does not require the solution of large linear systems. We apply the RFMP_NL to the nonlinear inverse problem of gravimetry, where gravitational data are inverted for the shape of the surface or inner layer boundaries of planetary bodies. This inverse problem is described by a nonlinear integral operator, for which we additionally provide the Fréchet derivative. Finally, we present two synthetic numerical examples to show that it is beneficial to apply the presented method to inverse gravimetric problems.