Making sound inferences from geomagnetic sounding

Abstract We examine the nonlinear inverse problem of electromagnetic induction to recover electrical conductivity. As this is an ill-posed problem based on inaccurate data, there is a critical need to find the reliable features of the models of electrical conductivity. We present a method for obtaining bounds on Earth’s average conductivity that all conductivity profiles must obey. Our method is based completely on optimization theory for an all-at-once approach to inverting frequency-domain electromagnetic data. The forward modeling equations are constraints in an optimization problem solving for the electric fields and the conductivity simultaneously. There is no regularization required to solve the problem. The computational framework easily allows additional inequality constraints to be imposed, allowing us to further narrow the bounds. We draw conclusions from a global geomagnetic depth sounding data set and compare with laboratory results, inferring temperature and water content through published Boltzmann–Arrhenius conductivity models. If the upper mantle is assumed to be volatile free we find it has an average temperature of 1409–1539  ° C. For the top 1000 km of the lower mantle, we find an average temperature of 1849–2008  ° C. These are in agreement with generally accepted mantle temperatures. Our conclusions about water content of the transition zone disagree with previous research. With our bounds on conductivity, we calculate a transition zone consisting entirely of Wadsleyite has 0.27  wt.% water and as we add in a fraction of Ringwoodite, the upper bound on water content decreases proportionally. This water content is less than the 0.4 wt.% water required for melt or pooling at the 410 km seismic discontinuity.

[1]  D. Bercovici,et al.  Whole-mantle convection and the transition-zone water filter , 2002, Nature.

[2]  Philip E. Gill,et al.  Practical optimization , 1981 .

[3]  R. Parker Ideal bodies for Mars magnetics , 2003 .

[4]  D. Dobson,et al.  The electrical conductivity of the lower mantle phase magnesiowustite at high temperatures and pressures , 2000 .

[5]  G. Backus,et al.  The Resolving Power of Gross Earth Data , 1968 .

[6]  R. Parker The existence of a region inaccessible to magneto‐telluric sounding , 1982 .

[7]  P. Stark,et al.  Strict Bounds on Seismic Velocity in the Spherical Earth , 1986 .

[8]  Robert L. Parker,et al.  A theory of ideal bodies for seamount magnetism , 1991 .

[9]  R. Parker,et al.  Assigning uncertainties in the inversion of NMR relaxation data. , 2005, Journal of magnetic resonance.

[10]  R. Parker The inverse problem of electromagnetic induction: Existence and construction of solutions based on incomplete data , 1980 .

[11]  Michael A. Saunders,et al.  SNOPT: An SQP Algorithm for Large-Scale Constrained Optimization , 2002, SIAM J. Optim..

[12]  R. Parker Geophysical Inverse Theory , 1994 .

[13]  W. D. Du Frane,et al.  Anisotropy of electrical conductivity in dry olivine , 2005 .

[14]  S. Karato,et al.  Water content in the transition zone from electrical conductivity of wadsleyite and ringwoodite , 2005, Nature.

[15]  D. Oldenburg,et al.  Inversion of Magnetotelluric Data for a One-Dimensional Conductivity , 1992 .

[16]  G. Egbert,et al.  Robust estimation of geomagnetic transfer functions , 1986 .

[17]  F. Guyot,et al.  Microstructures and iron partitioning in (Mg,Fe)SiO3 perovskite‐(Mg,Fe)O magnesiowüstite assemblages: An analytical transmission electron microscopy study , 1997 .

[18]  D. Oldenburg Funnel functions in linear and nonlinear appraisal , 1983 .

[19]  M. Hirschmann,et al.  Storage capacity of H2O in nominally anhydrous minerals in the upper mantle , 2005 .

[20]  Donald L. Turcotte,et al.  Mantle Convection in the Earth and Planets , 2001 .

[21]  G. Backus,et al.  Uniqueness in the inversion of inaccurate gross Earth data , 1970, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[22]  Michael A. Saunders,et al.  USER’S GUIDE FOR SNOPT 5.3: A FORTRAN PACKAGE FOR LARGE-SCALE NONLINEAR PROGRAMMING , 2002 .

[23]  M. S. Paterson,et al.  The determination of hydroxyl by infrared absorption in quartz, silicate glasses and similar materials , 1982 .

[24]  H. Keppler,et al.  Solubility of water in the α, β and γ phases of (Mg,Fe)2SiO4 , 1996 .

[25]  E. Haber,et al.  Inversion of 3D electromagnetic data in frequency and time domain using an inexact all-at-once approach , 2004 .

[26]  E. Haber,et al.  On optimization techniques for solving nonlinear inverse problems , 2000 .

[27]  S. Constable,et al.  Global mapping of the electrically conductive lower mantle , 1996 .

[28]  P. Weidelt The inverse problem of geomagnetic induction , 1973 .

[29]  S. Constable FAST TRACK PAPER: SEO3: A new model of olivine electrical conductivity , 2006 .

[30]  Michael A. Saunders,et al.  SNOPT: An SQP Algorithm for Large-Scale Constrained Optimization , 2005, SIAM Rev..

[31]  Xu,et al.  The effect of alumina on the electrical conductivity of silicate perovskite , 1998, Science.

[32]  William Rodi,et al.  Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion , 2001 .

[33]  Kathryn A. Whaler,et al.  Numerical methods for establishing solutions to the inverse problem of electromagnetic induction , 1981 .

[34]  R. Parker The Inverse Problem of Electrical Conductivity in the Mantle , 1971 .

[35]  S. Karato,et al.  The role of hydrogen in the electrical conductivity of the upper mantle , 1990, Nature.

[36]  R. Parker,et al.  Occam's inversion; a practical algorithm for generating smooth models from electromagnetic sounding data , 1987 .

[37]  N. Olsen Long-period (30 days-1 year) electromagnetic sounding and the electrical conductivity of the lower mantle beneath Europe , 1999 .

[38]  S. Constable Constraints on Mantle Electrical Conductivity from Field and Laboratory Measurements , 1993 .