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

Abstract

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.