Robust multiple-station magnetotelluric data processing

被引:445
作者
Egbert, GD
机构
[1] Coll. of Oceanic and Atmosph. Sci., Oregon State University, Corvallis
关键词
coherent noise; magnetotellurics; multivariate statistics; robust estimation; transfer functions;
D O I
10.1111/j.1365-246X.1997.tb05663.x
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
Although modern magnetotelluric (MT) data are highly multivariate (multiple components, recorded at multiple stations), commonly used processing methods are based on univariate statistical procedures. Here we develop a practical robust processing scheme which is based on multivariate statistical methods. With this approach we use data from all channels to improve signal-to-noise ratios, and to diagnose possible biases due to coherent noise. To illustrate our approach we use data from two-and three-station wide-band MT arrays from an area south of San Francisco, California, where contamination of the MT signal by spatially coherent cultural electromagnetic noise is severe at some periods. To deal with such coherent noise we adopt a two-stage procedure. In the first stage we focus on reducing the effects of incoherent noise, and testing for the presence of coherent noise. To this end we have developed a robust multivariate errors-in-variables (RMEV) estimator, which estimates background noise levels, cleans up outliers in all channels, and determines the 'coherence dimension' of the array data. In the absence of coherent noise, the coherence dimension of the data will be two (corresponding to two polarizations of the plane-wave MT source fields). In this case the RMEV estimator provides direct estimates of MT impedances and inter-station transfer functions. We show, with synthetic and real data examples, that in some cases these estimates can be significantly better than those obtained with more standard robust remote reference estimators. When MT data is severely contaminated by coherent noise (as for our example arrays for periods of 4-50 s) the coherence dimension of the data will exceed two. The RMEV estimate thus provides a clear warning of coherent noise contamination. Although there appears to be no completely general automatic way to deal with this circumstance, useful results can be obtained from severely contaminated data in some cases. We show in particular how the RMEV estimator can be adapted to separate the MT signal from coherent noise for two special cases: when at least one site is unaffected by coherent noise, and when coherent noise sources are intermittent. We give examples of significant improvements in MT impedance estimates obtained with the RMEV estimate for each of these cases.
引用
收藏
页码:475 / 496
页数:22
相关论文
共 35 条
[1]   ESTIMATING LINEAR STATISTICAL RELATIONSHIPS [J].
ANDERSON, TW .
ANNALS OF STATISTICS, 1984, 12 (01) :1-45
[2]  
ANDERSON TW, 1956, 3RD P BERK S MATH ST, V5, P111
[3]  
[Anonymous], 1977, MULTIVARIATE STAT IN
[4]   INTRODUCTION TO THE SPECIAL SECTION ON THE EMSLAB DEFUCA,JUAN EXPERIMENT [J].
BOOKER, JR ;
CHAVE, AD .
JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH AND PLANETS, 1989, 94 (B10) :14093-14098
[5]   SOME COMMENTS ON MAGNETOTELLURIC RESPONSE FUNCTION ESTIMATION [J].
CHAVE, AD ;
THOMSON, DJ .
JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH AND PLANETS, 1989, 94 (B10) :14215-14225
[6]   ON THE ROBUST ESTIMATION OF POWER SPECTRA, COHERENCES, AND TRANSFER-FUNCTIONS [J].
CHAVE, AD ;
THOMSON, DJ ;
ANDER, ME .
JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH AND PLANETS, 1987, 92 (B1) :633-648
[7]   ON THE SYNTHESIS OF A LARGE GEOMAGNETIC ARRAY FROM SMALL OVERLAPPING ARRAYS [J].
EGBERT, GD .
GEOPHYSICAL JOURNAL INTERNATIONAL, 1991, 106 (01) :37-51
[8]   MULTIVARIATE-ANALYSIS OF GEOMAGNETIC ARRAY DATA .1. THE RESPONSE SPACE [J].
EGBERT, GD ;
BOOKER, JR .
JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH AND PLANETS, 1989, 94 (B10) :14227-14247
[9]   ROBUST ESTIMATION OF GEOMAGNETIC TRANSFER-FUNCTIONS [J].
EGBERT, GD ;
BOOKER, JR .
GEOPHYSICAL JOURNAL OF THE ROYAL ASTRONOMICAL SOCIETY, 1986, 87 (01) :173-194
[10]  
EGBERT GD, 1987, THESIS U WASHINGTON