Geoelectric Exploration of the Purı́sima-Rumicruz District , Jujuy Province , Argentina

Purı́sima-Rumicruz district consists of several polymetallic veins, rich in copper and lead. The veins have a breccia texture and are located in shear bands, with locally high concentrations of sulfides in a carbonatic gangue. The host rock of the veins consists of a black shale sequence with thin levels of fine sands, with a low degree of metamorphism, and corresponds to the Acoite formation (Lower Ordovician). The depositional environment was an open clastic platform where storm processes prevailed. An exploratory program using induced polarization and resistivity with the dipole-dipole method was carried out with the objective of establishing geoelectric anomalies in depth, which would permit the identification of exploratory targets. The group of veins was divided into three sectors: Puŕısima, El Brechón, and La Nueva. Low-resistivity values correspond to faulty zones, some of which are related to surface veins. The presence of diagenetic pyrite in the shales of Acoite formation considerably increases background chargeability (induced polarization) values. However, small chargeability anomalies associated to low-resistivity values represent exploratory targets for the Purı́sima-Rumicruz district.


Introduction
Purísima-Rumicruz is a polymetallic mining district, located on the southern section of Cerro Colorado de Cobre, some 20 km southeast of the town of Abra Pampa (Cochinoca department, Jujuy province, Argentina).
Small-scale copper exploitation took place from the 1940's to the early 1970's, with a variable production during that period; reserves totalled some 6,750 t of mineral, with 5.5% Cu content.From the mid-1980s to date, a more ambitious exploitation plan has been in order.The firm that currently holds the mining concession is Bolland S. A.
A geophysical exploratory survey based on geoelectricinduced polarization (or chargeability), and resistivity, was carried out with the purpose of identifying exploratory targets for drilling aiming to increase the district's reserves.This time-efficient methodological approach has been a cornerstone of mining exploration for decades [1][2][3][4].The research objective of this study was to detect zones with high chargeability which are spatially linked to zones of low resistivity, a combination of characteristics that has been demonstrated to be indicative of metalliferous mineralizations associated with shear zone veins [5,6].

Geologic Setting
The Purísima-Rumicruz district is located on the limit of the geologic provinces of Puna and Cordillera Oriental.Mineral deposits are typically veins and are contained in rocks of the Acoite formation [7], of Tremadocian-Arenigian age, composed of sands and shales, and deposited on a platform in which storms dominated and determined a progradational sequence, and it increases in grain and strata thickness upsection.Cretaceous units from the Salta Group (Pirgua subgroup and Lecho formation) (Figure 1) are on tectonic contact on top of this unit.The host rocks of the mineralization are ductilly deformed by folding with a north-south axis, southwardly dipping.The fragile deformation is characterized by faults mainly with NW-SE and E-W orientations.
The veins present a breccia structure, with multiepisodic characteristics.According to the mineral association, there are two distinctive typologies: with quartz gangue and with calcite and baritine gangue.The former generally presents few metailiferous mineralizations, while the latter has sulfides and arsenides, and Cu-Pb sulfur salts with trace contents of Ni, Co, As, Ag, and U.
The vein system comprises some 3,500 m of outcropping structures.Veins run for short distances and present a rosarylike arrangement (Figure 1).The mineralization has been typified as five-element deposits [7,8].

Methods
Geoelectric profiles were drafted for three zones labeled as La Nueva, Purísima, and El Brechón, according to the principal outcrops of metalliferous veins in the district.A lineal multielectrode arrangement was used to register data in a dipole-dipole fashion.The equipment used was an IPR-12 time domain IP/resistivity receiver with an IPC-9/200 W power source, both manufactured by Scintrex.The equipment produced real-time automated readings, therefore enabling the operator to process as many samples as considered appropriate, within a variable number of cycles for each station.
Eight dipole potentials were read simultaneously, and resistivity and chargeability readings were logged in the equipment's memory.Additionally, the register of Kconstant readings for each position and the mean and standard deviation enable us to perform a reliability test on the induced polarization (IP) data stored.
The power source and the receptor are automatically synchronized, and the process is controlled by the first potential dipole.The source generates a continuous current, commuted in variable duration cycles.This allows to bypass spontaneous potential compensations, given that they are done statistically.
The apparent resistivity is then calculated by the formula where K is the geometric constant of the device.Values are expressed in units of ohm•m.Apparent chargeability is calculated on the time domain after the expression: where t 1 is the initial time of measurements; t 2 is the final time of measurements; t r is the measurement interval (t 2 − t 1 ); V p is the tension during emission; V s is the tension measured during the integration period.Chargeability values are expressed in units of mV/V.
In the present study, electrodes were separated by 25 m, and measurements for each station were taken for values ranging between n = 1 and n = 8 (n being the different depth levels for data recording).The duration of each cycle was 2 sec.400-600 m long profiles were drafted with a separation of 200 m between profiles.La Nueva section was covered with three 400 m long profiles with NNW-SSE orientation.The exploratory objective was to search for the lateral continuity of a vein first observed with a small surface sounding, which evidenced a 50 cm wide structure, rich in calcosine.In order to map regional plots of resistivity, a 0-1,500 ohm•m range with a 150 ohm.m interval was chosen.In the case of chargeability, fields were plotted with a 0-180 mV/V range and a 20 mV/V interval.
With the objective of identifying outcropping veins in the area and recognizing blind structures, the Purísima section was covered with nine 600 m long profiles.In order to map resistivity, the same intervals, as in the previous case, were used.In the case of chargeability, fields were plotted with a 0-80 mV/V range and an 8 mV/V interval, due to the shallower bottom of this section.
Carrying out a prospective study of the continuity of El Brechón vein, the presence of polarizable minerals and the identification of associated blind structures, five 600-m long profiles were drafted for El Brechón section.Having found lower resistivity values than in the other sections, a range of 0 to >750 ohm•m was chosen for mapping resistivity, plotting fields with 50 ohm•m intervals.As in the previous section, chargeability fields were plotted with a 0-80 mV/V range and an 8 mV/V interval.
Geotomo's RES2DINV software was used in order to produce the 2D inversion model for resistivity and chargeability profiles [9][10][11].The 2D model consists of a number of rectangular blocks with an arrangement linked to the distribution of data points in the pseudosections.The program uses a modeling subroutine in order to compute apparent resistivity values forward, a nonlineal minimum squares optimization technique for the inversion routine [12,13].The aim of the program is to produce an apparent resistivity pseudosection that adjusts to real measurements, determining resistivity for the rectangular blocks.The optimization method basically reduces the difference between the calculated apparent resistivity values and those measured, adjusting the resistivity of the model's blocks (Figures 2,  3, and 4).A measure of this difference is determined by the quadratic mean error (RMS).However, with the lowest RMS possible, the model can sometimes show great artifact variations in resistivity values, and it may not always be the most adequate model from a geologic perspective.In general, the most cautious approach is to choose an iteration model in which the RMS does not change significantly (Table 1).
The same method used in mining prospection for determining geochemical anomalies [14] was used to obtain normal (background) chargeability values, given that chargeability values present a distribution that fits well into a lognormal model (after Kolmogorov-Smirnov's goodness of fit test, P < 0.01 in all cases).A background upper limit or threshold was determined as the mean + 2 SD; chargeability values higher than the threshold were considered anomalous.

Results
Results from the inversion of geoelectric data and true resistivities and chargeabilities indicate that rocks of the Acoite formation present resistivity values characteristic of those of sedimentary rocks, with a mean value of 319 ohm•m.Induced polarization values obtained from the inversion model are generally very high, with some values >150 mV/V and a mean value for the three sections of 23 mV/V.The analysis of chargeability data indicates that there were highly significant differences among the mean values of each section (ANOVA, F (2;4688) = 287.32;P = 0.00) (Figure 5).Table 2 lists descriptive statistics for resistivity and chargeability for each zone.
4.1.La Nueva Section. Figure 6 shows LN200 resistivity and chargeability profiles for the section.In the first case, the distribution of resistivities is irregular.In general, the low resistivity zone is placed on the southern margin of the vein; to the south there is a homogeneous zone of moderate resistivities, whereas to the north the pattern is more heterogeneous, mainly presenting low-resistivities and interspersed small shallow zones of high resistivities.The chargeability profile shows a large anomalous zone located ∼50 m south of the vein and separated by a very low chargeability band which coincides with the low resistivity band.To the north of the vein, chargeability is comparatively moderate to low.

Purísima Section.
As in the previous case, two low chargeability zones located approximately on the main valley axis separate two sections with different resistivities (i.e., low resistivity in the south and high in the north).The deep zones of the veins comprised in profile no.600 are spatially associated to one of the low-resistivity zones located southward of the valley axis (Figure 7).With respect to chargeability, the anomalous zone is located some 50 m deep in the area of veins.In general, chargeabilities southward of the valley axis are higher than those to the north.

El Brechón Section.
Figure 8 shows profile no.600 for the section.In this section, the vein coincides spatially with a tabular, very low resistivity band, which is <50 ohm•m.
Although the resistivity distribution on the rest of the profile is heterogeneous, the NE section is generally more resistive than the SW section.Chargeabilities are also heterogeneously distributed; the more developed anomalous zone is clearly unrelated to  the section's vein.Below 50 m, a small anomaly coincides spatially with the vein.veins are located in multiepisodic faults.On the one hand, fractured bands can generate low resistivity zones due to a strong increment in the porosity from rock fractures [5,[17][18][19][20].On the other hand, this result, in rocks with different resistivities placed adjacently at each side of the fault [20].

Interpretation and Discussion
In the three sections of the district comprised in this study, low resistivity zones were recognized, which might be interpreted as the result of an increased effective porosity, itself the product of faults.The tabular structure of some of these zones, and their spatial association with mineralized structures (i.e., La Nueva and El Brechón), indicates that they might as well be filled-in fault zones.In the case of Purísima section, even though the low resistivity zone coincident with the deep zone of veins on profile no.600 does not have tabular structures, its alignment with the valley's axis indicates that it might be associated with a faulty zone.Resistivity values in these zones (0-150 ohm•m) are similar to those recorded by Suzuki et al. [21] in Paleozoic shales with a high degree of fragmentation due to faults.These authors also found that resistivity declines gradually as a function of the intensity of cataclasis.The mean resistivity value in shales is 450 ohm.m, similar to that of low-degree metamorphic rocks of Acoite formation.
In La Nueva and El Brechón sections, low resistivity bands associated with the veins clearly separate zones with different geoelectric characteristics.An interpretation of this is that the faults that harbor the current mineralization brought together different lithologic types in the same formation.

Chargeability.
The chargeability values recorded in this study are not frequent in sedimentary rocks.Some chargeability values were as high as those of porphyry copper and gold [22].Even though the presence of graphite, clay minerals, and the mineralogical composition of some basic rocks can produce high chargeability [23], in the case of Acoite formation this might be due to the relatively high abundances of disseminated pyrite (probably of diagenetic origin).The majority of pyrite grains in Acoite have framboidal texture, characteristic of precipitations under abiotic conditions in an euxinic environment during diagenesis [24].Therefore, they are genetically unrelated to the polymetallic mineralization present elsewhere in the mining district.The mean grain size of framboidal pyrite at Acoite formation is in the order of 10 μm, which favors an increment in induced polarization.Grain size is an important variable in IP records; a fine granulometry of polarizing minerals translates into greater surface and, accordingly, into an increase in chargeability (with constant volume) [25].
These characteristics of the host rock explain the high background chargeability values registered, which are not necessarily related to metallic mineralization in the district.Our discriminating analysis by area indicates that there are considerable differences in chargeability.
Both in La Nueva and El Brechón sections, the larger volume and anomalous chargeability zones (>90 mV/V and >59 mV/V, resp.) are spatially unrelated to the zones where the veins outcrop.The extent of their developments might point to an association with disseminating diagenetic pyrite in the host rock, rather than being associated with a hydrothermal polymetallic mineralization.On the other hand, in La Nueva section, the presence of outcroppings stained by iron oxides on top of the anomaly (surely a byproduct of meteoric altering of pyrite) reinforces the idea that zones with anomalous chargeability are linked to a marked pyritization of the host rock.
On Purísima and El Brechón sections, below 40 and 50 m, respectively, small chargeability anomalies appear spatially associated to veins, which can be attributed to an International Journal of Geophysics epigenetic hydrothermal mineralization within a pyritized host rock.Pittard and Bourne [26] reported similar results for Centenary gold deposit (Australia), in reference to increased IP values given by the minerals coming from the host rock, unrelated to economic mineralization.In this case, chargeability is due to the pyrite associated with gold, sharing a common hydrothermal origin, and to magnetite, a hypogenic mineral in doleritic host rocks.

Conclusions
The geoelectric exploration of the Purísima-Rumicruz district revealed a marked relationship between low resistivity bands and mineralizations located in faulty zones.The fracture plane of the rock enabled fluid circulation and mineral precipitation, while the reactivation of the faults produced multiepisodic breachy structures, which characterize the district.As a consequence of fractures, the resulting increased porosity and permeability would decrease resistivity to values below 150 ohm•m in bands that stand out from the host rocks.Moreover, these faulty zones become evident by linking together lithologic varieties of different geoelectric properties.
In general, the high IP values registered seem to be the product of a stronger diagenetic pyritization of the host rock rather than an epigenetic mineralization.Anomalous zones comprise large rock volumes spatially unrelated to outcropping veins and the faulty zones that acted as mineralization traps.However, small zones with anomalous chargeability located below depths ∼50 m over veins from Purísima and El Brechón sections constitute good exploratory targets.

Figure 2 :
Figure 2: Measured and calculated apparent pseudosections and inverse models sections of resistivity and chargeability for no.200 La Nueva section.

Figure 3 :
Figure 3: Measured and calculated apparent pseudosections and inverse models sections of resistivity and chargeability for no.600 Purísima section.

Figure 4 :
Figure 4: Measured and calculated apparent pseudosections and inverse models sections of resistivity and chargeability for no.600 Brechón section.

Figure 8 :
Figure8: Geological, resistivity, and chargeability profile no.600 for El Brechón section.The vein is spatially associated with a low-resistivity tabular area.An extensive area presenting anomalous chargeability is noted on the SW section of the profile and a smaller one in the deep zone (below ∼50 m) of the vein.

Table 1 :
Absolute error for each iteration.