Open AccessNon‐linear conjugate gradient inversion for global EM induction: resolution studies
Author(s) -
Kelbert Anna,
Egbert Gary D.,
Schultz Adam
Publication year - 2008
Publication title -
geophysical journal international
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0956-540X
DOI - 10.1111/j.1365-246x.2008.03717.x
Subject(s) - magnetotellurics , conjugate gradient method , inversion (geology) , jacobian matrix and determinant , electromagnetic induction , computation , geology , geodesy , solver , mathematical analysis , geophysics , mathematics , algorithm , physics , mathematical optimization , electrical resistivity and conductivity , seismology , electromagnetic coil , tectonics , quantum mechanics
SUMMARY We develop a non‐linear conjugate gradient inversion for global long period electromagnetic induction studies. The scheme requires computation of derivatives of the regularized penalty functional. We derive analytical and numerical expressions for these derivatives, and the associated Jacobian, and show how these can be efficiently implemented by generalizing and extending an existing finite difference forward solver. Using layered spherical harmonics to parametrize the model space, we invert a range of synthetic data sets to test the inversion, and to study vertical and horizontal resolution of currently available data sets. We conclude that the currently available long‐period global geomagnetic observatory data in the period range 5–107 d can resolve large scale (300–500 km vertically, thousands of km horizontally) heterogeneities in mantle electrical conductivity reliably at depths ∼ 670–1600 km. By extending induction response to 0.2–5 d (including daily variation periods), upper‐mantle structure could also be resolved.