The statistical distribution of magnetotelluric apparent resistivity and phase
Citable URI
http://hdl.handle.net/1912/2976As published
http://dx.doi.org/10.1111/j.1365246X.2007.03523.xDOI
10.1111/j.1365246X.2007.03523.xAbstract
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 noncentral chisquared with 2 degrees of freedom, with the noncentrality 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 shortertailed 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 nonparametric 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.
Description
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): 127132, doi:10.1111/j.1365246X.2007.03523.x.
Collections
Related items
Showing items related by title, author, creator and subject.

Exploring practical estimates of the ensemble size necessary for particle filters
Slivinski, Laura; Snyder, Chris (American Meteorological Society, 20151111)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 ... 
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. ... 
Correction of seafloor magnetotelluric data for topographic effects during inversion
Baba, Kiyoshi; Chave, Alan D. (American Geophysical Union, 20051214)The large contrast in electrical conductivity between seawater and the underlying seafloor accumulates boundary electric charges which can severely distort observed electric and magnetic fields. For marine magnetotelluric ...