Magnetotelluric data, stable distributions and impropriety: an existential combination
Chave, Alan D.
MetadataShow full item record
KeywordTime series analysis; Numerical approximations and analysis; Fractals and multifractals; Probability distributions; Magnetotellurics
The robust statistical model of a Gaussian core contaminated by outlying data that underlies robust estimation of the magnetotelluric (MT) response function has been re-examined. The residuals from robust estimators are systematically long tailed compared to a distribution based on the Gaussian, and hence are inconsistent with the robust model. Instead, MT data are pervasively described by the alpha stable distribution family whose variance and sometimes mean are undefined. A maximum likelihood estimator (MLE) that exploits the stable nature of MT data is formulated, and its two-stage implementation in which stable parameters are first fit to the data and then the MT responses are solved for is described. The MLE is shown to be inherently robust, but differs from the conventional robust estimator because it is based on a model derived from the data, while robust estimators are ad hoc, being based on the robust model that is inconsistent with actual data. Propriety versus impropriety of the complex MT response was investigated, and a likelihood ratio test for propriety and its null distribution was established. The Cramér-Rao lower bounds for the covariance matrix of proper and improper MT responses were specified. The MLE was applied to exemplar long period and broad-band data sets from South Africa. Both are shown to be significantly stably distributed using the Kolmogorov–Smirnov goodness of fit and Ansari-Bradley non-parametric dispersion tests. Impropriety of the MT responses at both sites is pervasive, hence the improper Cramér-Rao bound was used to estimate the MLE covariance. The MLE is shown to be nearly unbiased and well described by a Gaussian distribution based on bootstrap simulation. The MLE was compared to a conventional robust estimator, establishing that the standard errors of the former are systematically smaller than for the latter and that the standardized differences between them exhibit excursions that are both too frequent and too large to be described by a Gaussian model. This is ascribed to pervasive bias of the robust estimator that is to some degree obscured by their systematically large confidence bounds. Finally, a series of topics for further investigation is proposed.
Author Posting. © Author, 2014. This article is posted here by permission of The Royal Astronomical Society for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 198 (2014): 622-636, doi: 10.1093/gji/ggu121.
Showing items related by title, author, creator and subject.
Chave, Alan D.; Lezaeta, Pamela F. (John Wiley & Sons, 2007-07-29)The marginal distributions for the magnetotelluric (MT) magnitude squared response function (and hence apparent resistivity) and phase are derived from the bivariate complex normal distribution that describes the distribution ...
Uncertainty analysis of vegetation distribution in the northern high latitudes during the 21st century with a dynamic vegetation model Jiang, Yueyang; Zhuang, Qianlai; Schaphoff, Sibyll; Sitch, Stephen; Sokolov, Andrei P.; Kicklighter, David W.; Melillo, Jerry M. (John Wiley & Sons, 2012-02-13)This study aims to assess how high-latitude vegetation may respond under various climate scenarios during the 21st century with a focus on analyzing model parameters induced uncertainty and how this uncertainty compares ...
Further analysis of target strength measurements of Antarctic krill at 38 and 120 kHz : comparison with deformed cylinder model and inference of orientation distribution Chu, Dezhang; Foote, Kenneth G.; Stanton, Timothy K. (Acoustical Society of America, 1993-05)Data collected during the krill target strength experiment [J. Acoust. Soc. Am. 87, 16–24 (1990)] are examined in the light of a recent zooplankton scattering model where the elongated animals are modeled as deformed finite ...