Estimation of the Magnetotelluric Response Function: The Path from Robust Estimation to a StableMaximum Likelihood Estimator

Estimation of the Magnetotelluric Response Function: The Path from Robust Estimation to a... The robust statistical model of a Gaussian core contaminated by outlying data in use since the 1980s, and which underlies modern estimation of the magnetotelluric (MT) response function, is re-examined from first principles. The residuals from robust estimators applied to MT data are shown to be systematically long-tailed compared to a distribution based on the Gaussian and hence inconsistent with the robust model. Instead, MT data are pervasively described by the stable distribution family for which the Gaussian is an end member, but whose remaining distributions have algebraic rather than exponential tails. The validity of the stable model is rigorously demonstrated using a permutation test. A maximum likelihood estimator (MLE), including the use of a remote reference, that exploits the stable nature of MT data is formulated, and its two-stage implementation, in which stable parameters are first fit to the residuals, and then the MT responses are solved for, with iteration between them, is described. The MLE is inherently robust, but differs from a conventional robust estimator because it is based on a statistical model derived from the data rather than being ad hoc. Finally, the covariance matrices obtained from MT data are pervasively improper as a result of weak non-stationarity, and the Cramér–Rao lower bound for the improper covariance matrix is derived, resulting in reliable second-order statistics for MT responses. The stable MLE was applied to an exemplar broadband data set from northwest Namibia. The stable MLE is shown to be consistent with the statistical model underlying linear regression and hence is unconditionally unbiased, in contrast to the robust model. The MLE is compared to conventional robust remote reference and two-stage estimators, establishing that the standard errors of the former are systematically smaller than for either of the latter, and that the standardized differences between them exhibit excursions that are both too frequent and too large to be described by Gaussian statistics. These excursions are more prevalent when the tail thickness parameter of the stable distribution is small, and are attributed to rising bias in the robust estimator that is also consistent with the Berry–Esséen theorem that defines the rate of convergence to the central limit theorem value. An explanation for weak non-stationarity of MT data is proposed, and several extensions to the present work are described. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Surveys in Geophysics Springer Journals

Estimation of the Magnetotelluric Response Function: The Path from Robust Estimation to a StableMaximum Likelihood Estimator

Loading next page...
 
/lp/springer_journal/estimation-of-the-magnetotelluric-response-function-the-path-from-gymdf6ySN4
Publisher
Springer Netherlands
Copyright
Copyright © 2017 by Springer Science+Business Media B.V.
Subject
Earth Sciences; Geophysics/Geodesy; Earth Sciences, general; Astronomy, Observations and Techniques
ISSN
0169-3298
eISSN
1573-0956
D.O.I.
10.1007/s10712-017-9422-6
Publisher site
See Article on Publisher Site

Abstract

The robust statistical model of a Gaussian core contaminated by outlying data in use since the 1980s, and which underlies modern estimation of the magnetotelluric (MT) response function, is re-examined from first principles. The residuals from robust estimators applied to MT data are shown to be systematically long-tailed compared to a distribution based on the Gaussian and hence inconsistent with the robust model. Instead, MT data are pervasively described by the stable distribution family for which the Gaussian is an end member, but whose remaining distributions have algebraic rather than exponential tails. The validity of the stable model is rigorously demonstrated using a permutation test. A maximum likelihood estimator (MLE), including the use of a remote reference, that exploits the stable nature of MT data is formulated, and its two-stage implementation, in which stable parameters are first fit to the residuals, and then the MT responses are solved for, with iteration between them, is described. The MLE is inherently robust, but differs from a conventional robust estimator because it is based on a statistical model derived from the data rather than being ad hoc. Finally, the covariance matrices obtained from MT data are pervasively improper as a result of weak non-stationarity, and the Cramér–Rao lower bound for the improper covariance matrix is derived, resulting in reliable second-order statistics for MT responses. The stable MLE was applied to an exemplar broadband data set from northwest Namibia. The stable MLE is shown to be consistent with the statistical model underlying linear regression and hence is unconditionally unbiased, in contrast to the robust model. The MLE is compared to conventional robust remote reference and two-stage estimators, establishing that the standard errors of the former are systematically smaller than for either of the latter, and that the standardized differences between them exhibit excursions that are both too frequent and too large to be described by Gaussian statistics. These excursions are more prevalent when the tail thickness parameter of the stable distribution is small, and are attributed to rising bias in the robust estimator that is also consistent with the Berry–Esséen theorem that defines the rate of convergence to the central limit theorem value. An explanation for weak non-stationarity of MT data is proposed, and several extensions to the present work are described.

Journal

Surveys in GeophysicsSpringer Journals

Published: Aug 24, 2017

References

You’re reading a free preview. Subscribe to read the entire article.


DeepDyve is your
personal research library

It’s your single place to instantly
discover and read the research
that matters to you.

Enjoy affordable access to
over 12 million articles from more than
10,000 peer-reviewed journals.

All for just $49/month

Explore the DeepDyve Library

Unlimited reading

Read as many articles as you need. Full articles with original layout, charts and figures. Read online, from anywhere.

Stay up to date

Keep up with your field with Personalized Recommendations and Follow Journals to get automatic updates.

Organize your research

It’s easy to organize your research with our built-in tools.

Your journals are on DeepDyve

Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.

All the latest content is available, no embargo periods.

See the journals in your area

DeepDyve Freelancer

DeepDyve Pro

Price
FREE
$49/month

$360/year
Save searches from
Google Scholar,
PubMed
Create lists to
organize your research
Export lists, citations
Read DeepDyve articles
Abstract access only
Unlimited access to over
18 million full-text articles
Print
20 pages/month
PDF Discount
20% off