Sulfide minerals are a group of compounds with the presence of sulfur. The most abundant and economically members of this group are pyrites, pyrrhotite, chalcocite, galena, sphalerite, and the group of copper sulfides minerals (Vaughan & Corkhill, 2017). It is difficult to distinguish between massive sulfides and other conductors of unknown character and to determine the direction of the orebody (Telford, Geldart, & Sheriff, 1990b). Electrical Resistivity Tomography (ERT) is the most commonly used geophysical methods for imaging subsurface features and map geological variations. It can detect subsurface sulfide mineral distribution by studying the nature of the flow of electricity in the earth, for being uniquely able to see a large range of magnitudes which can vary up to 20 orders (Morgan, 2012). This method can delineate the various sources of mineralization according to their types (Evrard et al., 2018). The resistivity method is used to map spatial variations in subsurface electrical conductivity, while the induced polarization (IP) method is used to map changes in chargeability.
A variety of applications of resistivity and Induced polarization to mineral exploration were demonstrated in many researchers and significantly prove their success by locating the best drilling points (Langore, Alikaj, & Gjovreku, 1989; Oldenburg & Li, 1994; White et al., 2001; Wilkinson et al., 2006; Qian et al., 2007; Bery et al., 2012; Loke et al., 2013; Dandi, 2014; Tavakoli et al., 2016a; Ali et al., 2020). The geological structures of occurrence of sulfide minerals show low resistivity and high chargeability (Johnson & Anderson, 1981; Langore, Alikaj, & Gjovreku, 1989; White et al., 2001; Yoshioka & Zhdanov, 2005; Ja, Oo, & Arce, 2014; Côrtes et al., 2016). The IP can detect information of small conductive rocks lost in the resistivity models (Tavakoli et al., 2016b). It is also able to detect sulfide in the presence of disseminated mineral deposit (Olowofela, Ajani, & Oladunjoye, 2008).
This review attempts to provide a clear summary of the success and problems encountered with resistivity and IP methods in sulfide mineral exploration. We focus on determining the efficiency of Resistivity and IP in exploring sulfide minerals by identifying their geophysical properties, the fundamentals of geophysics, and the host rocks by highlighting strengths and weaknesses. It also presents a brief discussion of the limitations of the electrical resistivity and IP methods.
2. Para Genesis of Sulfides
Sulfide can be formed abiotically or biogenically.
2.1. Abiogenic Sulfides Formation
Sulfide minerals can form abiologically through the crust at high temperatures above about 200˚C (Plimer & Finlow-Bates, 1978; Trudinger, 1981; McDonald et al., 2018). High-grade sulfide deposition from abiogenic sulfide exhaled in hydrothermal fluids at oceanic spreading centres and rift zones (Trudinger, 1981). According to McDonald et al. (2018), Three phases are defined: The first is an early-stage of low-temperature (<250˚C), dominated by Fe-Zn-Cu-rich mineralization, the second stage present high-temperature (>300˚C) and is dominated by Cu-Fe-rich mineralization; and the third is at low-temperature phase (<150˚C) dominated by Fe-rich mineralization (McDonald et al., 2018). However, at this low temperature (120˚C - 140˚C), pyrrhotite can be found in a stable state (Plimer & Finlow-Bates, 1978).
Figure 1 shows the example of sulfides formation at high temperatures.
The older the deposit and the deeper its site of deposition, the more likely it is to contain abundant pyrrhotite, lack barite (Plimer & Finlow-Bates, 1978). Deposits which formed at shallow depth (<500 m) where circulation was probably restricted, pyrite is far more abundant than pyrrhotite (Plimer & Finlow-Bates, 1978).
2.2. Biogenic Sulfides Formation
In biogenic sulfides formation, bacteria activities are the primary source of sulfides. Bacterial sulfate reduction is a potential source of high sulfide concentrations in typical marine sediments. However, fixation is limited by a lack of reactive metals and could lead at best to the formation of relatively low-grade ores. Biogenic sulfide might also contribute to mineralization in hydrothermal situations (Trudinger, 1981). The sulfate reduction process takes place in much lower temperatures (<150˚C) (Trudinger, 1981; McDonald et al., 2018). In hydrothermal systems, thermophilic bacteria biogenically reduce seawater sulfate ( ) to sulfide (H2S) (McDonald et al., 2018). When these react with , they can form a very fine particle black precipitates poorly crystalline mackinawite (FeS) or a mixture with greigite (Fe3S4). With time, these iron minerals transform into pyrite. Figure 2 indicates three possible pathways of transformations.
Table 1 shows examples of sulfide minerals and their host rocks from both abiotically and biogenically.
Table 1. The major types of sulfide ore deposits and host rocks.
Abbreviations used are as follows: po-pyrrhotite, pn-pentlandite, py-pyrite, cpy-chalcopyrite, viol-violarite, cass-cassiterite, sph-sphalerite, gn-galena, bn-bornite, ttd-tetrahedrite, asp-arsenopyrite, cinn-cinnabar, uran-uraninite (Fontboté et al., 2017; Vaughan & Corkhill, 2017).
3. Theory and Basic Principles of the Resistivity and IP Methods
3.1. Theory and Basic Principles of the Resistivity Method
The electrical resistivity surveys consist of measurement of the resistivity distribution of subsurface heterogeneities. Electric current is generated into the ground, creating a stationary current flow in the earth. Potentials measured, produce features of the subsurface by providing information on their electrical properties.
For simplicity, all layers are expected to be horizontal. The electrical resistance, R of a material is related to its physical dimension, cross-sectional area, A and length, L through the resistivity, ρ.
The electrical resistance of the cylindrical body R (Ω), is given by:
where V is the potential and I is the current.
The conductivity σ (Sm−1), is the reciprocal of the resistivity.
In a homogeneous and isotropic half-space, the current density J (A/m2) is:
The potential V becomes:
The measured voltage with the injected current is converted into apparent resistivity ( ) by a geometric factor which considers the air-earth interface (Telford et al., 1990a)
where: k is a geometrical factor [m].
are the distances of the real and mirror effect of the ground surface at potential points (M, N) respectively.
is the measured difference potential at points M and N; I is the applied electric current.
3.2. Theory and Basic Principles of the IP Method
IP has been widely used in geophysical surveys for many decades to provide information about the complex conductivity (chargeability) of the subsurface. IP effect can be measured by both time-domain and the frequency domain. The most common is Time-Domain IP. When the current is applied to the ground, the ground can act as a capacitor and store electric charge. The imposed voltage is switched off; it first drops to an intermediate voltage and then gradually decays (Figure 3 and Figure 4). The polarization characteristics of the subsurface
Figure 3. Integrated decay voltage used as a measure of chargeability m.
Figure 4. Time window chargeability.
are produced by switch-off the transmitted current and measure the potential decay. This is known as the apparent chargeability (M). The most common measurement of the magnitude of the IP effect in the time domain is chargeability M.
The common chargeability expression is:
where Vp is the primary voltage, Vs represents secondary voltage, and Vt is the voltage decay with a time interval between t1 and t2. Units of chargeability are millivolts per volt (mV/V) or millisecond (ms).
In the frequency domain, measurements involve measuring the resistivity at varying frequencies one higher than the other, transmitted one after the other. The IP effect is the percentage frequency effect (PFE) defined as:
where and represent apparent resistivity at high and low frequencies respectively. PFE has units of percent.
The IP effect is influenced by the effective resistivity of the host rock and, consequently, by the type of electrolyte, pore size, temperature etc. Though, the metal factor (MF) can remove this effect to some extent (Telford, Geldart, & Sheriff, 1990b).
Metal Factor is defined as the amount of sulfide minerals contained in rocks, in which the amount depends on the value FE. Metal factor is formulated as follows:
1000 is an arbitrary factor chosen to give MF values appropriate values from 10 upwards, and MF has of mhos per foot.
The concept of metal factor arose early in the development of the IP method, because it was felt that low resistivity values would damp the value of PFE more than high resistivity values.
3.3. Characteristics of Different Arrays Configurations Types
Several conventional electrode arrays are used for resistivity and IP surveying. There is no “best” method that can match with all situations encountered in the fields (Daniels, 1977) but the choice depends on the aim of prospecting; the geological aspects of the subsurface, the availability of holes; and economics (Daniels, 1977; Bing & Greenhalgh, 2000). Table 2 describes common arrays configuration depending on the measured potential by assuming that the earth’s resistivity is homogeneous and isotropic.
Surface configuration: Surface geophysical methods present a low spatial resolution, which decreases with increasing depth of investigation. Moreover, this can hinder the exposure of small-scale features such as cavities or leachate plumes. The use of electrodes in boreholes and pole-dipole arrays are recommended to extend the adequate depth of investigation (White et al., 2001).
Hole-to-surface measurements are performed by setting a pole or dipole source down a borehole and making surface dipole measurements away from the source hole (Daniels, 1977).
Borehole: Vertical profiles can detect both in-hole and off-hole features, but the enlargement depends on the configuration (Palich & Qian, 2007).
Single-hole arrays (Vertical Resistivity Profiling (VRP)) use a pole or bipole source and potential receiver in the same borehole (Daniels, 1977).
Hole-to-hole measurements (Borehole Resistivity Tomography (BRT)), this is performed by installing a pole or a dipole current source in a borehole and placing a pole or a dipole potential-receiver nearby borehole (Daniels, 1977).
The cross-borehole resistivity measurements seem to be more useful than single-borehole measurements (Yang & Ward, 1985; Wright & Ward, 1987; Qian, Milkereit, & Gräber, 2007; Ali et al., 2020). Qian, Milkereit, and Gräber (2007) reported the success of the cross-borehole method, where the imaging of the massive sulfide mineralization was more extended compared to the single borehole. However, the cross-borehole anomalies are smaller by using a pole-pole array than the dipole-dipole array. Though, with the cross-borehole mise-à-la-masse method, the anomalies are more significant than other cross borehole methods (Wright & Ward, 1987).
Table 2. The common electrode arrays.
For shallow targets hole- to-surface measurements may locate the target with a minimal field effort (Daniels & Scott, 1981; Daniels, 1983; Tsourlos et al., 2011). In the case of borehole-to-surface, Pole-dipole and dipole-dipole configurations are preferred because other arrays such as Wenner, Wenner-Schlumberger, and Gradient Schlumberger are more prone to artifact distortion (White et al., 2001; Tsourlos et al., 2011). The lateral position of the target can easily set by hole-to-surface or conventional dipole-dipole surface arrays (Daniels & Scott, 1981). Ideally, the horizontal position of the body can be established by a conventional surface survey or a hole-to-surface survey while the vertical position of the body can be determined by hole-to-hole or single-hole measurement (Daniels & Scott, 1981). For the intersection of mineralization and drill holes, mise-á-la-masse surveys appear more promising for imaging the minerals’ orientation and strike extent (Tyne, 1980; Guo, Dentith, & Zhao, 2000). However, for near-surface with borehole resistivity surveys, accurate array locations and correct geometry factors are very important to account for borehole deviation effect and for accurate imaging (Guo, Milkreit, & Qian, 2014). For all, shorter electrode intervals are needed to better define shallow anomalies and mineralized veins (Mammo, 2013).
4. 1D, 2D, and 3D Imaging Surveys and Inversion Processes
4.1. 1D, 2D, and 3D Dimensional Surveys
The Electrical Resistivity and IP methods can provide 1D, 2D or 3D images of their distribution in the subsoil.
4.1.1. 1D Dimensional Surveys
Since the early 1920s to the late 1980s, the resistivity method’s first commercial use was one-dimensional (1-D) mapping method for the profiling and sounding methods (Ghosh, 1971; Ward, 1990; Meheni et al., 1996; Loke, 2011; Sharma & Verma, 2015). The method is quick and simple in application (Ghosh, 1971; Ogunbo, 2018), it has been widely used to investigate the ground for resource management, such as mineral, petroleum, and groundwater resources (Loke, 2011). With this survey, the subsurface should consist of horizontal layers (Loke, 2011; Loke et al., 2013). In homogeneous, 1-D survey provides good results (Karim, 2015). However, 1-D survey does not provide significant lateral changes, which may bring to gives inaccurate results or misleading results (Loke, 2011; Loke et al., 2013).
4.1.2. 2D Dimensional Surveys
2-D imaging surveys are mostly used in mineral exploration. Two-dimensional multi-electrode arrays produce simultaneous and display both horizontal and vertical variations in resistivity (Loke, 2011; Loke et al., 2013). Wenner, Wenner-Schlumberger, dipole-dipole, pole-pole or pole-dipole arrays are the most common use, depending on the respective position of the potential electrodes and the current electrodes.
The sensitivity to horizontal and vertical heterogeneities, depth of investigation, data coverage and signal strength change depends on type of configurations array and the background noise level (Samouëlian et al., 2005). Commonly 2-D inversion for resistivity data was used by many of the reviewed papers (White et al., 2001; Moreira et al., 2012; Goto et al., 2013; Mammo, 2013; Côrtes et al., 2016; Tavakoli et al., 2016b). Figure 5 clearly shows the correlation of Resistivity and IP and their efficiency in 2D prospecting.
Figure 5 indicates that both resistivity and IP show a large contrast between the electrical properties of sulfide and their host rocks. The 2-D and 3-D forward interpretations are based on finite-difference, finite-element, transmission surface, integral-equation, or hybrid finite-element/integral equation formulation to provide a reasonable estimate of the subsurface structure (Zonge et al., 2014).
4.1.3. 3D Dimensional Surveys
3D modeling is now playing an important role in very complex exploration
(Loke et al., 2013; Sharma & Verma, 2015). Three-dimensional electrical resistivity can be obtained by reconstructing a two-dimensional network of parallel pseudo-sections or using a square array of four electrodes (Bentley & Gharibi, 2004). Several studies show the effectiveness of the 3D survey compared to 1 D and 2 D. The 3D surveys perform well with high accuracy in short time (Bentley & Gharibi, 2004). When the dipole-dipole or pole-dipole electrode array is more affected by spatial effects, 3D inversion can control the validity of the pseudo-3D approach (Orfanos & Apostolopoulos, 2011). 3D surveys are useful in the creation of the amplitude and geometry of the complex resistivity anomalies (Bentley & Gharibi, 2004; Milkereit et al., 2008; Orfanos & Apostolopoulos, 2011; Mammo, 2013; Tavakoli et al., 2016a). Thus, the method can give correct position and the improved delineation of the target (White et al., 2001; Orfanos & Apostolopoulos, 2011; Moreira et al., 2012; Côrtes et al., 2016; Tavakoli et al., 2016b). Figure 6 and Figure 7 hold examples that show the effectiveness of 3D in sulfide explorations. 3Dinversion survey techniques provide high data acquisition rates, increase target resolution, offer a greater depth of penetration, and hence very cost-effective (White et al., 2001; Moreira et al., 2012; Côrtes et al., 2016; Tavakoli et al., 2016b).
The above images mainly contributed to better understanding the distribution of the sulfide mineralization around the mineralization zone. For example in Tavakoli et al.’s study, 3Dresistivity/IP was capable of providing a comprehensive overview of the electrical distribution in the upper ~450 m of the crust, which is of great importance for targeting VMS ore (Tavakoli et al., 2016b).
3D inversion is suitable for detecting electrical anisotropy while 1D and 2D were less pronounced (Meheni et al., 1996; Sretenovic & Arnaut, 2019). But, many facts need to be considered before choosing the method such as subsurface heterogeneity, background noise level, characteristics of an array, etc. Thus, multiple configurations are suggested to improve the reading of different subsoil features (Hesse, Jolivet, & Tabbagh, 1986; Côrtes et al., 2016).
4.2. Inversion Process for Sulfide Investigation
To obtain an accurate picture of the subsurface, it is necessary to carry out the pseudo section’s inversion (Loke & Barkert, 1995). The inversion is performed by forward modeling operator by calculating theoretical data from input data and thereby calculating derivatives of the data with regard to the parameter (Pelton, Rijo, & Swift, 1978; Madsen et al., 2018). The massive sulfides deposits are indicated by the coincidence of high chargeability and low resistivity values after inversion (Moreira et al., 2012). It is noticeable that a suitable inversion method must simultaneously minimize the effects of data error and model parameter errors (Ja, Oo, & Arce, 2014). Therefore, to produce a better inversion, there should be re-run adjusting default parameters by manual selection of settings such as vertical and horizontal weighting (Robertson & Hart, 2013). The “least-squares method” is the most method used amongst the papers analyzed in this review. The “least-squares method” seems to be the method that can successfully remove the distortions in the apparent resistivity pseudo section because it separates the electrode array geometry’s effect on the apparent resistivity values from that which results from the subsurface resistivity (Loke & Barkert, 1995). It also separates overlapping anomalies caused by different bodies (Loke & Barkert, 1995).
5. Interpretation of Electrical Resistivity and IP Variation for Sulfide Ores and Their Host Rocks
5.1. Electrical Properties of Sulfide Host Rocks
The effectiveness of detecting sulfide deposits depends on the evolution of their hosting with adequate contrast in their physical properties between the base metals and their host rocks. Resistivity appears to be negatively correlated with porosity, Basalts indicate high resistivity and low IP along the resistivity, volcanic tuff areas are predominated by high resistivity values (Tavakoli et al., 2016b; Komori et al., 2017; Evrard et al., 2018) etc. The resistivity variation of various host rocks and sulfide minerals is summarized in Table 3.
Sulfide in a sedimentary host may be derived directly from the mantle or by reduction of sulfate either by a chemical or biological mechanism (Trudinger, 1981). Resistivity and chargeability of sulfide mineralization in sedimentary host
Table 3. Electrical resistivities range sulfide minerals and host rocks.
expects to much contrast and detectable (Meju, 2002; Airo, 2015), because of their resistivity variations: Shale and clays present low values followed by sandstones with intermediate values, while coal and limestone beds have high values (Lau, 2000). Generally, sulfide in sedimentary hosts indicates low resistivity values with a strong correlation between the low resistivity and high chargeability areas for most of sedimentary host (Côrtes et al., 2016; Evrard et al., 2018). Massive sulfide surveyed around the submarine area usually shows low resistivity values lower than the seawater value (Goto et al., 2013; Komori et al., 2017; Ishizu et al., 2019). But, the sulfide hosted by quartz and carbonates leads to act as good electrical insulators with higher resistivity values (Katsube et al., 2003; Moreira et al., 2014; Côrtes et al., 2016; Han et al., 2016; Ishizu et al., 2019).
Seafloor massive sulfide (SMS) deposits are considered a high potential source of economic minerals s such as zinc, copper, tin, gold, and silver (Emsbo, 2007; Haroon et al., 2018). The electrical resistivity methods have revealed that SMS deposits exhibit a strong IP effect and lower resistivity than the surrounding host rock (Iijima & Sayanagi, 2013; Hördt et al., 2016; Ishizu et al., 2019). SMS area deposit indicates low resistivity values similar to or lower than the seawater (Iijima & Sayanagi, 2013).
Volcanogenic Massive Sulfide (VMS) deposits appear to have been formed by concentrations of base metal mineralizations by which their host rocks are dominated by submarine volcanic (Doyle & Allen, 2003; Iijima & Sayanagi, 2013; Ishizu et al., 2019). Most VMS deposits are classified based on metal content, tectonic context, or age (Barrie & Hannington, 1999; Doyle & Allen, 2003; Iijima & Sayanagi, 2013; Ishizu et al., 2019). VMS deposit may be characterized by a distinct zonation of the ore, gangue, and hydrothermally altered minerals. But, compared to sedimentary rock, igneous and metamorphic rocks present very high resistivity values (Lau, 2000; Airo, 2015). In the Iberian Pyrite belt, in Bathurst districts in Canada and Kuroko province of Japan, sulfide observed as volcano-sedimentary stratigraphic layers for being closely associated with volcanic lavas and pyroclastic rocks. They have formed a different relative participation of sedimentary/bacterio-genic and hydrothermal processes (Boulter, 1996; Barrie & Hannington, 1999; Sáez et al., 1999; Lentz & McCutcheon, 2006; Lentz, Thorne, & Beal, 2009; Walker, 2010; Piercey, Peter, & Herrington, 2015; Almod et al., 2019). VMS deposits (Sulfide in mafic or felsic hosted or occur as intercalation of basalts) are characterized by high electrical conductivity, high chargeability bordering zones, and sediments (such as graphitic shales) are more conductive than the associated volcanic host (Tavakoli et al., 2016b). However, the geochemical processes of weathering such as humidity, porosity, and hydrothermal alteration form clay can affect the target’s resistivity. Also, major faults, fracture zones, and related structural control features in VMS prospects may be saturated with saline fluids, making them electrically conductive (Meju, 2002; Moreira et al., 2014; Tavakoli et al., 2016b). Though, when mineralized, they will be recognized by their elevated chargeability (Hawke & Brooker, 2001).
5.2. Electrical Resistivity and Chargeability of Sulfide Ores
Common sulfides except sphalerite, have high electrical conductivities (Wells, 1914; Ross, 1957; Meju, 2002; Pearce, Pattrick, & Vaughan, 2006; Aradis, Annigan, & Ewing, 2007; Airo, 2015) Resistivity and IP of Sulfide deposits depend very strongly on their quantity and the degree of connection between the various mineral grains or veins in the relevant host rock. When sulfide grains are isolated, they are not good electrical conductors because conductivity increases with the concentration of metallic particles (Revil, Florsch, & Mao, 2015). The ore bodies constituting more than 80% sulfide are considered as having “massive” texture; it is clear they are good electrical conductors (Langore et al., 1989). 20% up to 80% of the total volume is considered as “veinlet” texture. Due to electrical connections between the veinlets, such sulfide are also good electrical conductors (Langore et al., 1989). Sulfide content in the disseminated ore is 10% - 20% (Langore et al., 1989), because of the disseminated texture, sulfide is not good electrical conductors. It is also noted that the chargeability increase depending on the volumetric content of metallic particles (Scott & West, 1969; Revil, Florsch, & Mao, 2015).
The Pyrite and galena are frequently available as reasonably large single crystals. Simultaneously, chalcopyrite, pyrrhotite, and arsenopyrite occur in a compact, microcrystalline forms (Parasnis, 1956) massive sulfide vary in thickness from 2 to 10 cm, massive chalcopyrite layers vary in thickness from 2 to 5 mm (Tivey et al., 1995). Pyrrhotite laths are 20 to 30 μm long, and 3 μm wide and occur in interstices and intergrown with outer edges of chalcopyrite grain. Forty to 100 μm clusters of sphalerite are intergrown with minor amounts of Pyrite (<10 μm) (Tivey et al., 1995).
Pyrrhotite is highly conducive and can significantly influence a sulfide body’s general electrical conductivity signature (Roach & Fitzpatrick, 2003). Pyrite and chalcopyrite also give potentials higher polarized (too conductive) state than galena (Wells, 1914; Roach & Fitzpatrick, 2003), but sphalerite is generally insultator with less conductivity and lower IP signatures than other sulfide minerals (Komori et al., 2017; Evrard et al., 2018). It can be conducive when associated with other conducive/chargeable sulfide minerals such as pyrite and pyrrhotite, even in small quantities (Roach & Fitzpatrick, 2003; Evrard et al., 2018). Even though (Hawke & Brooker, 2001)’s research showed that sphalerite is the primary mineral source of the induced polarisation effect observed in their investigated samples (Hawke & Brooker, 2001). It is also noted that Cu-bearing ores are likely to be more conductive than sphalerite-rich Zn ores (Airo, 2015). Chalcopyrite is highly conductive but usually constitutes only a small part of the total sulfide volume in many deposits. Table 4 clearly describes the conductivities and chargeability of major sulfide minerals.
As mentioned in Table 4, the geophysical responses of massive sulfide usually dominated by low resistivity and high chargeability. It is not easy to make any conclusive deduction about exact resistivity and IP values for each type of sulfide.
Table 4. Resistivity and chargeability of common sulfide minerals.
Different geological settings of the areas, make identical minerals to have different geophysical signatures or different minerals to produce the identical geophysical signature (Langore, Alikaj, & Gjovreku, 1989; Moreira et al., 2014; Côrtes et al., 2016). Some of the reviewed papers did not found differences between various sulfides, i.e., pyrite and chalcopyrite, in the spectral IP parameters with laboratory analysis (Langore, Alikaj, & Gjovreku, 1989). Also, mineralization may be of a chargeable nature because of the presence of other minerals which strongly associated with them like pyrite/marcasite (Côrtes et al., 2016; Evrard et al., 2018). Thus, minerals with high metallic conductivity even for disseminated sulfide ores, can have strong IP effects (Wang & Strangwayt, 1981). Though, the presence of graphite-bearing limestone, IP surveys is less useful for mapping sulfide mineralization (Guo, Dentith, & Zhao, 2000). But, graphite rarely occurs in massive sulfide. Table 5 presents examples of resistivity and IP methods in the exploration of sulfide deposits. It describes the types of surveys used, types of minerals and types of host rocks encountered, and shows their resistivity and IP responses.
6. Factors Affecting Exploration of Sulfides
Here are some factors that affect exploration of sulfide:
• The exploration of sulfides is potential sources of pollution, of the air and of surface waters and soils (Vaughan & Corkhill, 2017).
• The 1-D survey does not provide significant lateral changes, which may bring to gives inaccurate results or misleading results.
• When using 2D, the geological structures do not change in the direction perpendicular to the survey line, creating errors in interpretation.
• To delete negative IP data may cause the loss of essential information (Dahlin & Loke, 2015).
Table 5. Examples of resistivity and IP methods in the exploration of sulfide deposits.
• Interference from grounded structures such as pipelines, power lines, and fences can distort the signals (Sternberg, 2002).
• Sometimes pyrite in fine-grained sediments such as mudstones and shales may produce false anomalies (Sternberg, 2002).
• In the interpretation of resistivity data, the author should emphasize the fact that there are other causes of low resistivity other than sulfides and explain how to constrain the interpretation to be unique to sulfides.
It is much known about resistivity and IP surveys in mineral exploration. These methods are simple, fast, cheap, and effective for sulfide mineral exploration. The reviewed papers showed that these methods have been successfully used in sulfide exploration whereby many results showed that low resistivity correlates with high chargeability and vice versa. This review is the wealth of resistivity and IP properties exhibited by sulfide minerals and their host. This review covered different distinct techniques of investigation to make the reader understand the procedures and appropriate methods that can be applied in mining exploration by showing the success and problems affecting exploration. This research could be useful as a guide for the future geophysical studies such as geophysics, mineralogy or in mineral exploration. However, many authors did not give the values attributed to each type of minerals but the range value. Some failed to include the area covered during their survey. It could be suggested therefore to present the area covered during the geophysical survey or give the number of traverse and line with their space interval from which interested readers can get the idea of the area covered for future research.
This research study is supported by Hohai University.
Consent to Participate
I solemnly declare that this manuscript is the collection of original research data and does not contain any copied contents from other sources. Proper citations of previous works were made and to the best of our knowledge.
Consent for Publication
This work, whether in the same or different form, has not been presented and will not be submitted to any journal. I am fully aware of the legal outcome of this statement from me. Therefore, the authors consent for this paper to be published by this journal after acceptance.
Availability of Data and Material
All the data and materials used are available.
All the authors contributed significantly to writing this article. The authors read and approved the final manuscript.
 Ali, M., Sun, S., Qian, W., Bohari, A. D., Claire, D., Faruwa, A. R., & Zhang, Y. (2020). Borehole Resistivity and Induced Polarization Tomography at the Canadian Shield for Mineral Exploration in North-Western Sudbury. In E3S Web of Conferences (Vol. 168, p. 00002). EDP Sciences.
 Babaei, M. et al. (2020). Geostatistical Modeling of Electrical Resistivity Tomography for Imaging Porphyry Cu Mineralization in Takht-e-Gonbad Deposit, Iran. Journal of Mining & Environment, 11, 143-159.
 Barrie, C. T., & Hannington, M. D. (1999). Chap. 1. Classification of Volcanic-Associated Massive Sulfide Deposits Based on Host-Rock Composition. In C. T. Barrie, & M. D. Hannington (Eds.), Volcanic-Associated Massive Sulfide Deposits: Processes and Examples in Modern and Ancient Settings (8th ed.). Society of Economic Geologists.
 Bentley, L. R., & Gharibi, M. (2004). Two- and Three-Dimensional Electrical Resistivity Imaging at a Heterogeneous Remediation Site Case History Two- and Three-Dimensional Electrical Resistivity Imaging at a Heterogeneous Remediation Site. Geophysics, 69, 674-680.
 Bery, A. A., Saad, R., Mohamad, E. T., Jinmin, M., Azwin, I. N., Tan, N. M. A., & Nordiana, M. (2012). Electrical Resistivity and Induced Polarization Data Correlation with Conductivity for Iron Ore Exploration. The Electronic Journal of Geotechnical Engineering, 17, 3223-3233.
 Bing, Z., & Greenhalgh, S. A. (2000). Cross-Hole Resistivity Tomography Using Different Electrode Configurations. Geophysical Prospecting, 48, 887-912.
 Boulter, C. A. (1996). Extensional Tectonics and Magmatism as Drivers of Convection Leading to Iberian Pyrite Belt Massive Sulphide Deposits? Journal of the Geological Society, 153, 181-184.
 Dahlin, T., & Loke, M. H. (2015). Negative Apparent Chargeability in Time-Domain Induced Polarisation Data. Journal of Applied Geophysics, 123, 322-332.
 Dandi, H. (2014). Application of Magnetometery, Electrical Resistivity and Induced Polarization for Exploration of Polymetal Deposits, a Case Study: Halab Dandi, Zanjan, Iran. 2nd International Conference on Advances in Engineering Sciences and Applied Mathematics, Istanbul, 4-7.
 Daniels, J. J., & Scott, J. H. (1981). Interpretation of Hole-to-Surface Resistivity Measurements at Yucca Mountain, Nevada Test Site. Report, U.S. Geological Survey.
 Emsbo, P. (2007). Review of Metalliferous Sediments of the World—Fundamental Theory of Deep-Sea Hydrothermal Sedimentation. Economic Geology, 102, 1181-1181.
 Furman, A., & Ferre, T. P. A. (2003). A Sensitivity Analysis of Electrical Resistivity Tomography Array Types Using Analytical Element Modeling. Vadose Zone Journal, 2, 416-423.
 Galley, A. G. (1993). Characteristics of Semi-Conformable Alteration Zones Associated with Volcanogenic Massive Sulphide Districts. Journal of Geochemical Exploration, 48, 175-200.
 Ghosh, D. P. (1971). The Application of Linear Filter Theory to the Direct Interpretation of Geoelectrical Resistivity Sounding Measurements. Geophysical Prospecting, 19, 192-217.
 Gibson, H. L. et al. (2007). The VMS Model: Advances and Application to Exploration Targeting. In B. Milkereit (Ed.), Proceedings of Exploration 07: Fifth Decennial International Conference on Mineral Exploration (pp. 713-730).
 Goto, T. et al. (2013). Resistivity Survey of Seafloor Massive Sulfide Areas in the Iheya North Area, off Okinawa, Japan. Proceedings of the 11th SEGJ International Symposium, Yokohama, 18-21 November 2013, 298-301.
 Han, M. et al. (2016). Induced Polarization Imaging Applied to Exploration for Low-Sulfidation Epithermal Au-Ag Deposits, Seongsan Mineralized District, South Korea. Journal of Geophysics and Engineering, 13, 817-823.
 Haroon, A. et al. (2018). Marine Dipole-Dipole Controlled Source Electromagnetic and Coincident-Loop Transient Electromagnetic Experiments to Detect Seafloor Massive Sulphides: Effects of Three-Dimensional Bathymetry. Geophysical Journal International, 215, 2156-2171.
 Hassan, E., Rai, J. K., & Anekwe, U. O. (2017). Geoelectrical Survey of Ground Water in Some Parts of Kebbi State Nigeria, a Case Study of Federal Polytechnic Bye-Pass Birnin Kebbi and Magoro Primary Health Center Fakai Local Government. Geosciences, 7, 141-149.
 Hawke, P. J., & Brooker, P. I. (2001). Can Sphalerite Be a Polarisable Mineral? An Example from the Century Zn-Pb Deposit. ASEG Extended Abstracts, 2001, 1-4.
 Hesse, A., Jolivet, A., & Tabbagh, A. (1986). New Prospects in Shallow Depth Electrical Surveying for Archaeological and Pedological Applications. Geophysics, 51, 585-594.
 Idornigie, A., Olorunfemi, M. O., & Omitogun, A. A. (2006). Electrical Resistivity Determination of Subsurface Layers, Subsoil Competence and Soil Corrositivity at an Engineering Site Location in Akungba-Akoko, South Western Nigeria. Ife. Journal of Science, 8, 159-177.
 Iijima, K., & Sayanagi, K. (2013). Resistivity Survey of Seafloor Massive Sulfide Areas in the Iheya North Area, off Okinawa, Japan. Proceedings of the 11th SEGJ International Symposium, Yokohama, 18-21 November 2013, 298-301.
 Johnson, G. R., & Anderson, L. A. (1981). A Laboratory Study of Some Physical Properties of Sulfide Ores in Igneous and Metamorphic Rocks from the Burnt Nubble Area, Somerset County, Maine. US Department of the Interior, Geological Survey.
 Karim, H. H. (2015). Implementation of Electrical Resistivity Imaging (ERI) Technique for Near Surface Investigation. The 2nd International Conference on Buildings, Construction and Environmental Engineering, October 2015, 101-108.
 Komori, S. et al. (2017). Depth Profiles of Resistivity and Spectral IP for Active Modern Submarine Hydrothermal Deposits: A Case Study from the Iheya North Knoll and the Iheya Minor Ridge in Okinawa Trough, Japan. Earth, Planets and Space, 69, Article No. 114.
 Kumar, D. et al. (2016). Integrated Geophysical and Geological Studies for Mineral Prospecting in Betul-Chhindwara Belt (BCB), Central India. Journal Geological Society of India, 87, 383-396.
 Langore, L., Alikaj, P., & Gjovreku, D. (1989). Achievements in Copper Sulphide Exploration in Albania with IP and EM Methods. Geophysical Prospecting, 37, 975-991.
 Lentz, D. R., & McCutcheon, S. R. (2006). The Brunswick No. 6 Massive Sulfide Deposit, Bathurst Mining Camp, Northern New Brunswick, Canada: A Synopsis of the Geology and Hydrothermal Alteration System. Exploration and Mining Geology, 15, 1-34.
 Madsen, L. M., Fiandaca, G., Christiansen, A. V., & Auken, E. (2018). Resolution of Well-Known Resistivity Equivalences by Inclusion of Time-Domain Induced Polarization Data. Geophysics, 83, E47-E54.
 McDonald, M. J. et al. (2018). Mineral Assemblages, Textures and in Situ Sulphur Isotope Geochemistry of Sulphide Mineralization from the Cyprus-Type Ice Volcanogenic Massive Sulphide (VMS) Deposit, Yukon, Canada. Minerals, 8, 501.
 Meheni, Y., Guérin, R., Benderitter, Y., & Tabbagh, A. (1996). Subsurface DC Resistivity Mapping: Approximate 1-D Interpretation. Journal of Applied Geophysics, 34, 255-269.
 Milkereit, B. et al. (2008). Geophysical Imaging of a “Blind” Zn-Pb-Ag Deposit. In 70th EAGE Conference and Exhibition incorporating SPE EUROPEC 2008 (pp. 9-12). European Association of Geoscientists & Engineers.
 Moreira, C. A. et al. (2012). Geoelectrical Prospection of Disseminated Sulfide Mineral Occurrences in Camaquã Sedimentary Basin, Rio Grande do Sul State, Brazil. Brazilian Journal of Geophysics, 30, 169-179.
 Moreira, C. A. et al. (2014). Geological and Geophysical Data Integration for Delimitation of Mineralized Areas in a Supergene Manganese Deposits. Geofísica Internacional, 53, 199-210.
 Olowofela, J. A., Ajani, O. O., & Oladunjoye, M. A. (2008). Application of Induced Polarization Method to Delineate Sulphide Ore Deposit in Osina Area of Benue State, Nigeria. Ife Journal of Science, 10, 137-150.
 Orfanos, C., & Apostolopoulos, G. (2011). 2D-3D Resistivity and Microgravity Measurements for the Detection of an Ancient Tunnel in the Lavrion Area, Greece. Near Surface Geophysics, 9, 449-457.
 Pearce, C. I., Pattrick, R. A. D., & Vaughan, D. J. (2006). Electrical and Magnetic Properties of Sulfides. In D. J. Vaughan (Ed.), Sulfide Mineralogy and Geochemistry. Reviews in Mineralogy and Geochemistry (Vol. 61, pp. 127-180). De Gruyter.
 Piercey, S. J., Peter, J. M., & Herrington, R. J. (2015). Zn-Rich Volcanogenic Massive Sulphide (VMS) Deposits. In S. Archibald, & S. Piercey (Eds.), Current Perspectives on Zinc Deposits (pp. 37-57). Irish Association for Economic Geology.
 Plimer, I. R., & Finlow-Bates, T. (1978). Relationship between Primary Iron Sulphide Species, Sulphur Source, Depth of Formation and Age of Submarine Exhalative Sulphide Deposits. Mineralium Deposita, 13, 399-410.
 Qian, W., Milkereit, B., & Gräber, M. (2007). Borehole Resistivity Tomography for Mineral Exploration. In 69th EAGE Conference and Exhibition incorporating SPE EUROPEC 2007 (pp. 11-14). European Association of Geoscientists & Engineers.
 Qian, W., Milkereit, B., McDowell, G., Stevens, K., & Halladay, S. (2007). Borehole Resistivity Logging and Tomography for Mineral Exploration. In B. Milkereit (Ed.), Proceedings of Exploration 07: Fifth Decennial International Conference on Mineral Exploration (pp. 1115-1118).
 Revil, A., Florsch, N., & Mao, D. (2015). Induced Polarization Response of Porous Media with Metallic Particles—Part 1: A Theory for Disseminated Semiconductors. Geophysics, 80, D525-D538.
 Sáez, R., Pascual, E., Toscano, M. & Almodóvar, G. R. (1999). The Iberian Type of Volcano-Sedimentary Massive Sulphide Deposits. Mineralium Deposita, 34, 549-570.
 Sretenovic, B., & Arnaut, F. (2019). Inadequacy of 1D, 2D and 3D Resistivity Inverse Modelling in the Presence of Electrical Anisotropy. Earth Sciences, 8, 102-116.
 Sternberg, B. K. (2002). A Review of Some Experience with the Induced-Polarization/Resistivity Method for Hydrocarbon Surveys: Successes and Limitations. Geophysics, 56, 1522-1532.
 Tavakoli, S., Bauer, T. E. et al. (2016a). Deep Massive Sulphide Exploration Using 2D and 3D Geoelectrical and Induced Polarization Data in Skellefte Mining District, Northern Sweden. Geophysical Prospecting, 1-18.
 Tavakoli, S., Thunehed, H. et al. (2016b). Targeting VMS Ore Deposits within the Upper 1.5 km of the Crust in the Central Skellefte District, Northern Sweden, Using 2D and 3D-Resistivity/IP Data. Journal of Geophysical Prospecting.
 Tivey, M. K. et al. (1995). Deducing Patterns of Fluid Flow and Mixing within the TAG Active Hydrothermal Mound Using Mineralogical and Geochemical Data. Journal of Geophysical Research, 100, 12527-12555.
 Tsourlos, P., Ogilvy, R., Papazachos, C., & Meldrum, P. (2011). Measurement and Inversion Schemes for Single Borehole-to-Surface Electrical Resistivity Tomography Surveys. Journal of Geophysics and Engineering, 8, 487-497.
 Walker, J. A. (2010). Stratigraphy and Lithogeochemistry of Early Devonian Volcano-Sedimentary Rocks Hosting the Nash Creek Zn-Pb-Ag Deposit, Northern New Brunswick. In G. L. Martin, & E. A. Smith (Eds.), Geological Investigations in New Brunswick for 2009 (pp. 52-97). New Brunswick Department of Natural Resources; Lands, Minerals and Petroleum Division, Mineral Resource Report 2010-1.
 Ward, S. H. (1988). The Resistivity and Induced Polarization Methods. In Symposium on the Application of Geophysics to Engineering and Environmental Problems (pp. 109-250). Environmental & Engineering Geophysical Society.
 Ward, S. H. (1990). Resistivity and Induced Polarization Methods. In S. H. Ward (Ed.), Geotechnical and Environmental Geophysics (pp. 147-190). Investigations in Geophysics No. 5 in SEG.
 White, R. M. S., Collins, S., Denne, R., Hee, R., & Brown, P. (2001). A New Survey Design for 3D IP Inversion Modelling at Copper Hill. Exploration Geophysics, 32, 152-155.
 Wilkinson, P. B. et al. (2006). Optimization of Array Configurations and Panel Combinations for the Detection and Imaging of Abandoned Mineshafts Using 3D Cross-Hole Electrical Resistivity Tomography. Journal of Environmental and Engineering Geophysics, 11, 213-221.
 Yoshioka, K., & Zhdanov, M. S. (2005). Three-Dimensional Nonlinear Regularized Inversion of the Induced Polarization Data Based on the Cole-Cole Model. Physics of the Earth and Planetary Interiors, 150, 29-43.
 Zonge, K., Wynn, J. et al. (2014). Chap. 9. Resistivity, Induced Polarization, and Complex Resistivity. In Near-Surface Geophysics Part 1: Concepts and Fundamentals (pp. 265-300). Society of Exploration Geophysicists.