Please use this identifier to cite or link to this item:
|Title:||Data space conjugate gradient inversion for 2-D magnetotelluric data|
Oregon State University
|Keywords:||Earth and Planetary Sciences|
|Citation:||Geophysical Journal International. Vol.170, No.3 (2007), 986-994|
|Abstract:||A data space approach to magnetotelluric (MT) inversion reduces the size of the system of equations that must be solved from M × M, as required for a model space approach, to only N × N, where M is the number of model parameter and N is the number of data. This reduction makes 3-D MT inversion on a personal computer possible for modest values of M and N. However, the need to store the N × M sensitivity matrix J remains a serious limitation. Here, we consider application of conjugate gradient (CG) methods to solve the system of data space Gauss-Newton equations. With this approach J is not explicitly formed and stored, but instead the product of J with an arbitrary vector is computed by solving one forward problem. As a test of this data space conjugate gradient (DCG) algorithm, we consider the 2-D MT inverse problem. Computational efficiency is assessed and compared to the data space Occam's (DASOCC) inversion by counting the number of forward modelling calls. Experiments with synthetic data show that although DCG requires significantly less memory, it generally requires more forward problem solutions than a scheme such as DASOCC, which is based on a full computation of J. © 2007 The Authors Journal compilation © 2007 RAS.|
|Appears in Collections:||Scopus 2006-2010|
Files in This Item:
There are no files associated with this item.
Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.