On the statistics of magnetotelluric rotational invariants
Citable URI
http://hdl.handle.net/1912/6283As published
https://doi.org/10.1093/gji/ggt366DOI
10.1093/gji/ggt366Abstract
The statistical properties of the Swift skew, the phasesensitive skew and the WAL invariants I1−I7 and Q are examined through analytic derivation of their probability density functions and/or simulation based on a Gaussian model for the magnetotelluric response tensor. The WAL invariants I1−I2 are shown to be distributed as a folded Gaussian, and are statistically well behaved in the sense that all of their moments are defined. The probability density functions for Swift skew, phasesensitive skew and the WAL invariants I3−I4, I7 and Q are derived analytically or by simulation, and are shown to have no moments of order 2 or more. Since their support is semiinfinite or infinite, they cannot be represented trigonometrically, and hence are inconsistent with a Mohr circle interpretation. By contrast, the WAL invariants I5−I6 are supported on [ − 1, 1], and are inferred to have a beta distribution based on analysis and simulation. Estimation of rotational invariants from data is described using two approaches: as the ratio of magnetotelluric responses that are themselves averages, and as averages of sectionbysection estimates of the invariant. Confidence intervals on the former utilize either Fieller's theorem, which is preferred because it is capable of yielding semiinfinite or infinite confidence intervals, or the less accurate delta method. Because sectionbysection averages of most of the rotational invariants are drawn from distributions with infinite variance, the classical central limit theorem does not pertain. Instead, their averaging is accomplished using the median in place of the mean for location and an order statistic model to bound the confidence interval of the median. An example using real data demonstrates that the ratio of averages approach has serious systematic bias issues that render the result physically inconsistent, while the average of ratios result is a smooth, physically interpretable function of period, and is the preferred approach.
Description
Author Posting. © Author, 2013. This article is posted here by permission of Oxford University Press on behalf of The Royal Astronomical Society for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 196 (2014): 111130, doi:10.1093/gji/ggt366.
Collections
Related items
Showing items related by title, author, creator and subject.

The statistical distribution of magnetotelluric apparent resistivity and phase
Chave, Alan D.; Lezaeta, Pamela F. (John Wiley & Sons, 20070729)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 ... 
Segmentation of plate coupling, fate of subduction fluids, and modes of arc magmatism in Cascadia, inferred from magnetotelluric resistivity
Wannamaker, Philip E.; Evans, Rob L.; Bedrosian, Paul A.; Unsworth, Martyn J.; Maris, Virginie; McGary, R Shane (John Wiley & Sons, 20141111)Five magnetotelluric (MT) profiles have been acquired across the Cascadia subduction system and transformed using 2D and 3D nonlinear inversion to yield electrical resistivity cross sections to depths of ∼200 km. Distinct ... 
Electrical structure beneath the northern MELT line on the East Pacific Rise at 15°45′S
Baba, Kiyoshi; Tarits, Pascal; Chave, Alan D.; Evans, Rob L.; Hirth, Greg; Mackie, Randall L. (American Geophysical Union, 20061116)The electrical structure of the upper mantle beneath the East Pacific Rise (EPR) at 15°45′S is imaged by inverting seafloor magnetotelluric data obtained during the Mantle ELectromagnetic and Tomography (MELT) experiment. ...