2D interpretation of the magnetotelluric data on Mighan plain of ARAK


1 Associate Professor, Earth Physics Department, Institute of Geophysics, University of Tehran, Iran

2 M.Sc. in Geophysics, Earth Physics Department, Institute of Geophysics, University of Tehran

3 Associate Professor, Physics Department of Geophysics, Faculty of Science, University of Arak , Iran


Magnetotelluric (MT) method is a passive electromagnetic technique that uses the natural, time varying electric and magnetic field components measured at right angles at the surface of the earth to make inferences about the earth’s electrical structure which, in turn, can be related to the geology tectonics and subsurface conditions. Reflection and refraction of electromagnetic (EM) signals at both horizontal and vertical interfaces separate media of different electrical parameters. Electromagnetic methods have been developed and employed to recognize the geological features and particularly fault zones in many regions. To achieve higher lateral resolution and also greater depth penetration, the MT method is one of the most effective electromagnetic techniques to create image of the subsurface structures electrically.
    In 2011 wide frequency range of magnetotelluric measurements were carried out at Mighan plain in the southern part of the Markazi province in Iran to understand the crustal electrical conductivity of the region by putting emphasis on locating the geological structures and recognizing the bedrock and probable fault. The electric and magnetic field components were acquired along a profile at 6 stations with a 1000-meter distance between stations using GMS05 (Metronix, Germany) systems. Three magnetometers and two pairs of non-polarizable electrodes were connected to this five-channel data logger. The experimental setup included four electrodes distributed at a distance of 100 m in north–south (Ex) and east–west (Ey) direction.
    Measurements of the horizontal components of the natural electromagnetic field were used to construct the full complex impedance tensor, Z, as a function of frequency. Using the effective impedance, determinant apparent resistivities and phases were computed and used for the inversion. MT data were processed using a code from Smirnov (2003) aiming at a robust single site estimate of electromagnetic transfer functions. As the area of study is populated and close to noise sources, the recorded data has not good quality which justifies the low coherency between the electric and magnetic channels. We performed 1D inversion of the determinant data using a code from Pedersen (2004) for all sites. The 2D modeling was applied to the data to explain the data if their responses fitted the measured data within their errors. Generally, the better the fit between measured and predicted data, the more reliable model. The 2D inversion of the TE and DET-mode data using a code from Siripunvaraporn and Egbert (2000) were performed. The data were calculated as apparent resistivities and phases. Apparent resistivity and phase data exhibited fairly different characteristics in the TE and DET -modes.  we used the model obtained from the TE -mode data as an interpretation model. The resistivity model obtained from the TE -mode is consistent with the geological model of the Mighan region down to five kilometers.
The 2D models significantly illustrate two conductive blocks and a fault structure and resolved layers with sharp resistivity contrasts. As significant results, in collaboration with geological information about the presence of the Tabarteh fault, the conductivity features can be attributed to the fault. Besides, a probable hidden fault is also recognizable. The bedrock was also detected with high apparent resistivity by the two dimensional model.


Main Subjects

میرزایی، م. وقدیمی، ف، 2006، گزارش زمین‌شناسی: بررسی گسل‌های تبرته و تلخاب، دانشگاه اراک.
Cantwell, T., and Madden, T. R., 1960, Preliminary report on crustal magnetotelluric measurements: Geophysics, 65(12), 4202-4205.
Jones, A.G. & I. Dumas (1993). Electromagnetic images of a volcanic zone. Phys. Earth
Planet. Inter., 81, 289–314.
Oklahoma, 2(B), 641-711.Wannamaker, P.E., 1986. Electrical conductivity of water-undersuturated crustal melting:J Geophys. Res. 91(B6), 6321-6327.
Oskooi, B., 2004. "A broad view on the interpretation of electromagnetic data (VLF, RMT, MT, CSTMT)". PhD Thesis, Uppsala University, Sweden.
Pedersen, L. B., 2004.” Determination of the regularization level of truncated singular-valubhe decomposition inversion”, The case of 1D inversion of MT data: Geophys. Prospect, 52, 261-270.
Pedersen, L. B., and Engels, M., 2005, Routine 2D inversion of magnetotelluric data using the determinant of the impedance tensor: Geophysics, 70, G33-G41.
Siripunvaraporn, W., Egbert, G., 2000. “An efficient data-subspace inversion method for 2-D magnetotelluric data”, Geophysics, 65, 791-803.
Smirnov, M. Yu., 2003. ” Magnetotelluric data processing with a robust statistical procedure having a high breakdown point”. Geophysics. J. Int, 152, 1-7.
Swift, C. M., 1967. "A magnetotelluric investigation of electrical conductivity anomaly in the southwestern United States" .PhD Thesis Massachusetts Institute of Technology, Cambridge, MA.
Vozoff, K., 1991. "The magnetotelluric method, in Electromagnetic methods in applied geophysics.". M. N. Nabighian, Ed., Society of Exploration Geophysicists, Tulsa, Oklahoma, 2(B), 641-711.