The statistical distribution of magnetotelluric apparent resistivity and phase
Chave, Alan D.
Lezaeta, Pamela F.
MetadataShow full item record
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 of response function estimates when the Gauss–Markov theorem is satisfied and the regression random errors are normally distributed. The distribution of the magnitude squared response function is shown to be non-central chi-squared with 2 degrees of freedom, with the non-centrality parameter given by the squared magnitude of the true MT response. The standard estimate for the magnitude squared response function is biased, with the bias proportional to the variance and hence important when the uncertainty is large. The distribution reduces to the exponential when the expected value of the MT response function is zero. The distribution for the phase is also obtained in closed form. It reduces to the uniform distribution when the squared magnitude of the true MT response function is zero or its variance is very large. The phase distribution is symmetric and becomes increasingly concentrated as the variance decreases, although it is shorter-tailed than the Gaussian. The standard estimate for phase is unbiased. Confidence limits are derived from the distributions for magnitude squared response function and phase. Using a data set taken from the 2003 Kaapvaal transect, it is shown that the bias in the apparent resistivity is small and that confidence intervals obtained using the non-parametric delta method are very close to the true values obtained from the distributions. Thus, it appears that the computationally simple delta approximation provides accurate estimates for the confidence intervals, provided that the MT response function is obtained using an estimator that bounds the influence of extreme data.
Author Posting. © The Authors, 2007. This article is posted here by permission of John Wiley & Sons for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 171 (2007): 127-132, doi:10.1111/j.1365-246X.2007.03523.x.
Showing items related by title, author, creator and subject.
Slivinski, Laura; Snyder, Chris (American Meteorological Society, 2015-11-11)Particle filtering methods for data assimilation may suffer from the “curse of dimensionality,” where the required ensemble size grows rapidly as the dimension increases. It would, therefore, be useful to know a priori ...
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, 2014-11-11)Five magnetotelluric (MT) profiles have been acquired across the Cascadia subduction system and transformed using 2-D and 3-D nonlinear inversion to yield electrical resistivity cross sections to depths of ∼200 km. Distinct ...
Baba, Kiyoshi; Tarits, Pascal; Chave, Alan D.; Evans, Rob L.; Hirth, Greg; Mackie, Randall L. (American Geophysical Union, 2006-11-16)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. ...