Three-dimensional magnetotelluric inversion using conjugate gradients

Abstract

We have developed an inversion procedure that uses conjugate gradient relaxation methods. Although one can generalize the method to all inverse problems, we demonstrate its use to invert magnetotelluric data for 3-D earth models. This procedure allows us to bypass the actual computation of the sensitivity matrix <b>A</b> or the inversion of the <b>A</b>T<b>A</b> term. In fact, with the relaxation approach, one only needs to compute the effect of the sensitivity matrix or its transpose multiplying an arbitrary vector. We show that each of these requires one forward problem with a distributed set of sources either in the volume (for <b>A</b> multiplying a vector) or on the surface (for <b>A</b>T multiplying a vector). This significantly reduces the computational requirements needed to do a 3-D inversion. For this paper, we have simplified the boundary conditions by assuming the model is repeated in the horizontal directions, but this is not a necessary constraint of the method. The algorithm reduces data errors to the 2 per cent level for noise-free synthetic 3-D magnetotelluric data.