Document Type : Original Research Paper


Faculty of Mining, Petroleum and Geophysics Engineering, Shahrood University of Technology, Shahrood, Iran


In any geophysical exploration, the final goal is to achieve an accurate image of the relevant underground property. In order to achieve this, the geophysical operation using the electrical resistivity and induced polarization (IP) methods is conducted to explore the sub-surface sulfide mineralization. Considering the mineralization evidence in the Kaboudan area near the Bardeskan city, first, geophysical surveying of the polymetallic deposit is carried out using the electrical resistivity and IP methods by employing the rectangle array in order to detect the electrical anomalies in the area. Then for delineation of the identified anomalies and investigation of the mineralization in the area, the 2D resistivity and chargeability cross-sections are prepared and interpreted with the help of the geological information. This geophysical survey in the area has led to the identification of several potential areas for mineralization. Then in order to obtain a detailed picture of the sub-surface mineralization and an overview of the in-depth mineralization distribution, a 3D modeling of the acquired data is made, and the results of this modeling are shown in 3D forms. The mineralization zones are identified in the studied area from their high chargeability values as well as the low to medium electrical resistivity amounts. This can be attributed to the metal mineralization and the presence of sulfide minerals in the mineralization zones. Mineralization in many places of the studied area is determined with an approximate east-west trend as well as somewhat varying the intensities of the electrical resistivity and chargeability amounts. The geological and drilling information obtained from the area confirm the interpretations.


[1]. Abdulsamad, F., Revil, A., Ahmed, A.S., Coperey, A., Karaoulis, M., Nicaise, S. and Peyras, L. (2019). Induced polarization tomography applied to the detection and the monitoring of leaks in embankments. Engineering Geology, 254: 89-101.
[2]. Aizebeokhai, A.P. (2009). Geoelectrical resistivity imaging in environmental studies. In Appropriate Technologies for Environmental Protection in the Developing World (pp. 297-305). Springer, Dordrecht.
[3]. Deceuster, J. and Kaufmann, O. (2012). Improving the delineation of hydrocarbon‐impacted soils and water through induced polarization tomographies: A field study at an industrial waste land. Journal of Contaminant Hydrology, 136, 25-42.
[4]. Di Maio, R., Fais, S., Ligas, P., Piegari, E., Raga, R. and Cossu, R. (2018). 3D geophysical imaging for site-specific characterization plan of an old landfill. Waste Management, 76: 629-642.
[5]. Eftekharnezhad, J., Hamzehpour, B., Aghanabati, A. and Baroyant, V. (1976). Geological map of Kashmar, 1:250000 scale. Geological Survey of Iran, Tehran, Iran.
[6]. Gallistl, J., Weigand, M., Stumvoll, M., Ottowitz, D., Glade, T. and Orozco, A.F. (2018). Delineation of sub-surface variability in clay-rich landslides through spectral induced polarization imaging and electromagnetic methods. Engineering geology. 245: 292-308.
[7]. Gazoty, A., Fiandaca, G., Pedersen, J., Auken, E., Christiansen, A.V. and Pedersen, J.K. (2012). Application of time domain-induced polarization to the mapping of litho-types in a landfill site. Hydrology and Earth System Sciences. 16 (6): 1793-1804.
[8]. Goad, R.E., Mumin, A.H., Duke, N.A., Neale, K.L., Mulligan, D.L. and Camier, W.J. (2000). The NICO and Sue-Dianne Proterozoic, iron oxide-hosted, polymetallic deposits, Northwest Territories: Application of the Olympic Dam model in exploration. Exploration and Mining Geology. 9 (2): 123-140.
[9]. Google Earth Pro (32-bit). (July 21, 2020). Kaboudan area, Razavi Khorasan Province, Iran. 40 S 588898.00 m E 3920760.45 m N, Eye alt 2.51 Km. <> (Accessed August 20, 2020).
[10]. Gurin, G., Tarasov, A., Ilyin, Y. and Titov, K. (2013). Time-domain spectral-induced polarization of disseminated electronic conductors: Laboratory data analysis through the Debye decomposition approach. Journal of Applied Geophysics. 98: 44-53.
[11]. Gurin, G.V., Tarasov, A.V., Il’in, Y.T. and Titov, K.V. (2017). Transient characteristics of induced polarization in inhomogeneous media (from results of 2D numerical simulation). Russian Geology and Geophysics. 58 (5): 624-634.
[12]. Heritiana, A.R., Riva, R., Ralay, R. and Boni, R. (2019). Evaluation of flake graphite ore using self-potential, electrical resistivity tomography (ERT) and induced polarization methods in east coast of Madagascar. Journal of Applied Geophysics.
[13]. Hördt, A., Blaschek, R., Kemna, A. and Zisser, N. (2007). Hydraulic conductivity estimation from induced polarisation data at the field scale—the Krauthausen case history. Journal of Applied Geophysics. 62 (1): 33-46.
[14]. Karimpour, M. (2010). Taknar massive sulfide, U-Pb zircon geochemistry and Sr-Nd isotopic characteristic of Bornaward granitoids, Iran. 18th symposium of crystallography and mineralogy of Iran. Tabriz, Iran.
[15]. Kearey, P., Brooks, M. and Hill, I. (2013). An introduction to geophysical exploration. John Wiley and Sons.
[16]. Kirsch, R. (2006). Groundwater geophysics (Vol. 493). Berlin: Springer.
[17]. Knödel, K., Lange, G. and Voigt, H.J. (2007). Environmental geology: handbook of field methods and case studies. Springer Science & Business Media.
[18]. Loke, M.H. (2010). Tutorial: RES2DINV ver. 3.59, Rapid 2-D Resistivity & IP inversion using the least-squares method. Malaysia, Geotomo Software.
[19]. Mary, B., Saracco, G., Peyras, L., Vennetier, M., Mériaux, P. and Camerlynck, C. (2016). Mapping tree root system in dikes using induced polarization: Focus on the influence of soil water content. Journal of Applied Geophysics. 135: 387-396.
[20]. Pelton, W.H., Ward, S.H., Hallof, P.G., Sill, W.R. and Nelson, P.H. (1978). Mineral discrimination and removal of inductive coupling with multifrequency IP. Geophysics. 43 (3): 588-609.
[21]. Qi, Y., El-Kaliouby, H., Revil, A., Ahmed, A. S., Ghorbani, A. and Li, J. (2019). 3D modeling of frequency-and time-domain electromagnetic methods with induced polarization effects. Computers & geosciences. 124: 85-92.
[22]. Revil, A., Qi, Y., Ghorbani, A., Ahmed, A.S., Ricci, T. and Labazuy, P. (2018). Electrical conductivity and induced polarization investigations at Krafla volcano, Iceland. Journal of Volcanology and Geothermal Research. 368: 73-90.
[23]. Safari, M. (2009). Prospecting report of 1:25000 Kaboudan (Bardeskan area). Geological Survey of Iran, Tehran, Iran.
[24]. Safari, M. (2010). Introducing new mineralization of massive sulfide deposits in Kaboudan area (north of Bardeskan) using feature analysis method and element carat maps. 1st Conference of Iranian Society of Economic Geology. Mashhad, Iran.
[25]. Safari, M. and Akbari-moghadam, M. (2011). Exploration of massive sulfide deposits in Kaboudan area, north of Bardeskan. 29th National Geosciences Congress, Geological Survey of Iran, Tehran, Iran.
[26]. Selley, R.C., Cocks, R. and Plimer, I. (2004). Encyclopedia of geology. Academic Press.
[27]. Shahrabi, M., Hosseini, K., Shabani, K., Massomi, R. and Haddadan, M., (2006). Geological map of Bardaskan, 1:100000 scale. Geological Survey of Iran, Tehran, Iran.
[28]. Shao, Z., Revil, A., Mao, D. and Wang, D. (2018). Finding buried metallic pipes using a non-destructive approach based on 3D time-domain induced polarization data. Journal of Applied Geophysics. 151: 234-245.
[29]. Veeken, P.C., Legeydo, P.J., Davidenko, Y.A., Kudryavceva, E.O., Ivanov, S.A. and Chuvaev, A. (2009). Benefits of the induced polarization geo-electric method to hydrocarbon exploration. Geophysics. 74 (2): B47-B59.
[30]. White, R.M.S., Collins, S. and Loke, M.H. (2003). Resistivity and IP arrays, optimized for data collection and inversion. Exploration Geophysics. 34 (4): 229-232.