Method for obtaining response functions from noisy magnetotelluric data using frequency-domain independent component analysis

Abstract

The magnetotelluric (MT) method has been used for visualizing subsurface resistivity structures and more recently for monitoring resistivity changes. However, electromagnetic data often include cultural noise, which can cause errors in the estimation of MT response functions and subsurface resistivity structure analysis. Frequency-domain independent component analysis (FDICA) offers advantages for MT data processing particularly because this method can extract hidden components in the observed data. These components can be decomposed into natural MT signals and cultural noise so that the noise effect in the recovered MT data is reduced. FDICA is applied to MT data acquired at the Kakioka Magnetic Observatory in Japan. The apparent resistivity and phase curves are obtained with small estimated errors between periods of 7 and 12,000 s, although the length of the time-series data is limited. The curves are smoother than those obtained using a conventional method. Various types of synthetic noise are added to the time series at Kakioka to test the noise-reduction performance of FDICA for MT data with high noise contamination. The results demonstrate that FDICA can be used to estimate MT response functions with high accuracy even under conditions in which more than half of the time-series data are contaminated by noise.