Copper, Zinc, and Lead Mineral Prospectivity Mapping in the North of Tafresh, Markazi Province, Central Iran, Using the AHP-OWA Method ()
1. Introduction
Stream sediment geochemistry is a useful technique for identifying important deposits [1] . The aim of this research is to examine the efficacy of stream sediment geochemical survey for locating copper, zinc, and lead in the area under study.
The elemental correlations were calculated by applying the Pearson correlation method, factor analysis, and classical method to elemental concentration data to study copper, zinc, and lead paragenesis and to define geochemical anomaly thresholds to find copper, zinc, and lead prospects in the region.
The use of Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) satellite imagery in mineral exploration has developed in recent years [2] . In this research, hydrothermal alteration zones were detected using ASTER data in different algorithms: Relative absorption Band Depth (RBD), Minimum Noise Fraction (MNF), and Least Squares Fit (LS-Fit).
Geospatial data, such as geochemical, geological, structural, and remotely sensed data, are particularly suitable to extract information that is useful in mineral prospecting [3] - [9] .
The application of the analytic hierarchy process (AHP) method is useful for binary index overlay modeling, as well as for multi-class index overlay modeling and fuzzy logic modeling [10] .
The ordered weighted averaging (OWA) model provides a continuous fuzzy between fuzzy intersection (AND) and union (OR) [11] [12] . The AHP is used to deduce attribute weights and the OWA operator function is used to produce a wide variety of decisions [13] . Combining information layers with AHP-OWA is an effective method to produce mineral potential mapping in the region. Field observations have confirmed the results.
2. Geology Setting
The area under study is situated between 50˚E to 50˚15'E and 34˚30'N to 35˚N in the Urmia-Dokhtar structural zone of Iran [14] . The Urumieh-Dokhtar mag- matic belt (UDMB) hosts porphyry copper deposits [15] - [21] .
The oldest outcrops in the region are the Middle Triassic sandy shale rocks (the Nayband formation). Eocene volcanic rocks (mostly tuff and andesitic rocks with limestone interlayers compositions) with a thickness of approximately 3000 meters cover most of the area. Intrusive tonalite to granodiorite has been injected into Eocene volcanic rocks and has metamorphosed in greenschist facies [22] . The Lower Red, Qom, and Upper Red formations, with sedimentary rocks, outcrop in the northern parts of the area. Eocene dykes with micro-diorite and andesitic and Neogene dykes with basaltic-andesitic compositions have greatly spread in the area under study (Figure 1).
3. Materials and Methods
520 stream sediment samples gathered and analyzed by the Geological Survey of Iran (GSI) in 2006. The samples were chemically analyzed by ICP-MS for 19 elements (Ag, As, Au, Ba, Bi, Cd, Co, Cr, Cu, Hg, Ni, Pb, Sb, Sc, Sn, Sr, V, W, Y, and Zn), with detection limit of 0.05 ppm for copper, zinc, and lead. Thompson
Figure 1. The simple geology map of study area and its location in UDMB.
control diagrams and relative error calculations were used to control the accuracy of the analysis.
The elemental correlation coefficients and paragenesis were found by applying bivariate and multivariate statistical analysis [23] [24] to the geochemical data using SPSS 18 statistical software package. Bivariate matrix of correlation coefficients was calculated by the Pearson correlation method [25] .
Hydrothermally altered minerals and lithological mapping have been successfully recognized through remote sensing instruments for mineral exploration [26] [27] [28] [29] [30] ; therefore, in this research, the ASTER image has been used and data analyses were done by ENVI 4.2 software. This work mainly focused on mapping of hydrothermal alterations in the volcano-plutonic belt; therefore, the sedimentary rocks and quaternary deposits outside the volcanic and plutonic rocks were masked. Evidence maps were produced from the data layers using ArcGIS 9.3 software and potential maps were generated by the AHP-OWA method using IDRISI 17.0 software. After data processing, the validity of the results was confirmed by field work.
4. Geochemical Data
In this study, Pearson correlation confirmed a strong correlation between Pb
Table 1. Pearson correlation coefficients for chemical analyses of samples from the area under study.
and As, Pb and Ag, and Pb and Ba (Table 1). Reduction of variables was accomplished by factor analysis, and 19 elements were classified into three factors. Table 2 shows that:
Factor 1 (43.64% of variance) is strongly positively correlated to Y, Cd, Sr, Sb, Co, Hg, and Cu and negatively correlated to Pb and Sn, with very weak correlation to Zn.
Factor 2 (21.85% of variance) is strongly positively correlated to Cr, Ni, Sc, and V and negatively correlated to Bi.
Factor 3 (19.15% of variance) is strongly positively correlated to As, Ag, W, and Ba, negatively correlated to Zn, and weakly correlated to Pb. Statistically significant correlation was not found between Pb and Zn in the area under study.
The catchment basins produced by a digital elevation model (DEM) in ARCGIS software (Figure 2).
We separated the background and anomaly populations by classical statistics method [31] as three classes for zinc, lead, and copper in the area under study (Figure 3). is normally considered as the threshold value, as a possible anomaly, and as a probable anomaly.
5. Remote Sensing
Recent years, remote sensing images are used to identify hydrothermally altered
Table 2. Varimax rotated factor matrix for 19 elements in 520 stream sediment samples of the area under study.
rocks and to produce the maps of geology, faults, and fractures [32] . In this study, an ASTER image has been used to detect alteration zones.
5.1. ASTER Data
ASTER is one of the recently designed generations of medium resolution multispectral satellite remote systems that are reported to have substantial applications in geological mapping [2] . ASTER has 14 bands: Three visible and near infrared radiation (0.52 to 0.86 μm) (VNIR, 15 m); six shortwave infrared radiation (1.6 to 2.43 μm) (SWIR, 30 m) and five thermal infrared radiations (8.125 to 11.65 μm) (TIR, 90m) [33] [34] . Many investigations have been carried out using ASTER [35] [36] . The ASTER spectral properties were used to map iron oxide, argillic, phyllic, and propylitic alteration zones by RBD, MNF, and LS-Fit methods in the present study.
5.2. Image Processing Methods on ASTER Imagery
5.2.1. Relative Absorption Band Depth (RBD)
An RBD image is produced by dividing the sum of the maximum reflectance
Figure 2. Locations of 520 stream sediment samples and their watersheds in the area under study.
Figure 3. Geochemical anomaly maps of Zn, Pb and Cu in the area under study.
bands to maximum absorption band [37] [38] . Phyllic alteration with abundance of illite and sericite has strong absorption at 2.20 µm (ASTER band 6) [39] . The argillic alteration (kaolinite and alunite), has strong absorption at 2.17 µm (ASTER band 5) [40] , and propylitic alteration dominated by chlorite, epidote and carbonates, has strong absorption in the 2.335 µm (ASTER band 8) (Figure 4). Table 3 shows the common band combinations that were used to create the maps of argillic, phyllic, and propylitic alteration zones in the present study (Figure 5).
5.2.2. Minimum Noise Fraction (MNF)
Because spectral anomalies are related to hydrothermal mineralization, this method is more interest to exploration geologists [41] [42] [43] . MNF Bands 7, 6, 4, and 3 have been used for, propylitic argillic, phyllic, and iron oxide alterations in the present study (Figure 6).
5.2.3. Least Squares Fitting (LS-Fit)
The superiority of this method compared to other methods is the reduction of noise in images [44] . To identify alteration zones, visible and near-infrared (VNIR) bands used as the input bands and band 1 as the modeled band for Iron oxide, band 4 for argillic, band 6 for phyllic and band 9 for propylitic (Figure 7).
6. Generation of Information Layers for Study Area
6.1. Fault Layer
Mineralization is more probable along a fault; therefore, faults are considered as one of the most significant elements in mineralization [45] . The presence of these geological features indicates enhanced structural permeability of the rocks in the subsurface, which facilitates upward migration of ground water that has come in contact with and has leached substances from buried deposits. These arguments suggest that the significance of multi-element stream sediment anomalies in sample catchment basins can be screened for or examined further using fault density as a factor [46] . In this research, the map of fault distributions
Figure 4. Laboratory spectra of muscovite, kaolinite, alunite, epidote, calcite, and chlorite resampled to ASTER band passes.
Table 3. Common ratio and band combinations for mapping alteration zones in the area under study.
Figure 5. Argillic, phyllic, propylitic and iron oxide images were prepared by RBD method.
Figure 6. The argillic, phyllic, propylitic and iron oxide images were prepared by MNF method.
Figure 7. Argillic, phyllic, propylitic and iron oxide images were prepared by Ls-fit method.
was derived from the 1:100,000 geological map of Tafresh (Figure 8). The fault map was used to create a fault density map (Figure 8) and reclassified into ten classes in ArcGIS. The highest number had the strongest preference (Figure 8).
6.2. Lithology Layer
This layer was derived from the 1:100,000 geological map of Taferesh. The most important component in this layer is the volcanic rocks (i.e., andesitic basalt, basaltic andesite, and andesite). The units of the geologic map were reclassified into ten classes according to their mineralization ability (Figure 9). The sedimentary rocks and younger geological units, such as the older and younger terraces, were assigned the smallest value.
6.3. Dike Layer
The abundance of dikes interpreted from geological map in the area under study was applied to create separate linear structural features. A dike evidence map was generated by reclassifying the dike buffer map (Figure 10) into nine classes.
6.4. Alteration Layer
The present research used different algorithms to excavate the most common hydrothermal alterations, namely iron oxide, phyllic, argillic, and propylitic, by processing ASTER images. The alteration evidence map was produced by combining the results from different algorithms and reclassifying to ten classes (Figure 11).
6.5. Geochemical Layer
The geochemical evidence map for three elements, zinc, lead, and copper, was produced after providing geochemical maps, combining them, and reclassifying to ten classes in ArcGIS software (Figure 12).
Figure 8. The fault layer of the area under study: fault map (left), fault density map (center) and fault density reclassified map (right).
Figure 9. The lithology layer of the area under study.
Figure 10. The dike layer of the area under study: Dike buffer map (left) and dike buffer reclassified map (right).
Figure 11. The alteration layer of the area under study.
Figure 12. (a) Pb geochemical layer. (b) Zn geochemical layer. (c) Cu geochemical layer. (d) The geochemical evidence layer of the area under study
7. Integration of the Layers and Mapping the Potential Area
The mineral potential map of the area under study multiple layers (lithology, fault, dike, alteration, and geochemistry) were produced by combining them by the AHP-OWA method.
7.1. The Analytical Hierarchy Process (AHP) Method
7.1.1. The AHP Hierarchy
The hierarchical diagram has four main levels: goals, objectives, attributes, and alternatives [47] . The hierarchical diagram was produced for information layers of the area under study (Figure 13).
7.1.2. Pairwise Comparisons
Figure 13. Hierarchical diagram of exploration layers in the area under study.
Table 4. Scale for pair-wise comparisons [48] .
Note: Adapted from Saaty (1980).
Table 5. Pairwise comparison matrix of the level of objectives and calculated weights.
7.2. The Ordered Weighted Averaging (OWA)
The OWA prepares a proper status of AND, Average and OR degree [49] [50] [51] (Figure 14).
7.3. The AHP-OWA Model
To identify the most suitable area in the present research, three outcomes have been generated by AHP-OWA procedures at AND, OR, and Average for achieving the favorable areas (Figure 15).
Figure 14. Decision making strategy space in OWA method [52] .
Figure 15. Favorable area maps which have been generated by AHP-OWA procedures in the area under study.
8. Field Checking
The accuracy of the results from AHP-OWA (Average) model was verified by selecting four accessible locations in the area under study (Figure 16).
8.1. Target 1
This area is located along a thrust fault (50˚10'38"E, 34˚53'08"N) that placed older rocks (Eocene Volcanic Rocks) over younger Oligo-Miocene sedimentary rocks. Argillic, propylitic, and iron oxide alterations affected the rocks. Copper mineralization (malachite) was observed in a barite vein that is perpendicular to the thrust fault trend. The Zn-bearing minerals were detected with Zinc-Zap, which reacts with Zn-bearing minerals to form bright scarlet stains [53] , while Fe sulfide minerals, primarily pyrite, cause a dark blue color [54] and other minerals do not react with the zap. In this site, testing the barite vein with this zap gave blue reaction paint, shows the presence of Fe-subfield minerals, and a red reaction color, and indicates the presence of Zn-bearing minerals (Photo 1).
8.2. Target 2
This area is located along the fault between Oligo-Miocene sedimentary rocks and Eocene andesitic rocks (50˚09'30"E, 34˚52'37"N). Propylitic and iron oxide alterations have been seen in this area. The barite veins with copper mineralization (malachite) have occurred in volcanic rocks. In this site, barite veins tested with Zinc-zap gave a blue reaction color, indicating the presence of Fe-subfield minerals (Photo 2).
8.3. Target 3
This site is located in 50˚04'51"E, 34˚48'27"N. The most important observations here are that numerous Eocene dioritic dikes, ranging in width from 1 to about 3 meter, have intruded into the Eocene volcanic rocks. Argillic alteration has affected a wider area around Eocene dikes. Field checking revealed that mineralization has not occurred in this site (Photo 3).
8.4. Target 4
The last site for checking was located in 50˚12'25"E, 34˚44'06"N. In this site, in-
Photo 1. The observations of the target 1 in the area under study.
Photo 2. The observations of the target 2 in the area under study
trusive igneous rock (g) (tonalite-granodiorite) has intruded into the Eocene volcanic rocks. Evidence of extensive hydrothermal alteration and mineralization has not been observed in this outcrop. Pyrite has altered to iron oxide minerals in these rocks (Photo 4).
9. Results
The area under study is located in Urmia-Dokhtar structural zone of Iran. The
Photo 3. The observations of the target 3 in the area under study.
Photo 4. The observations of the target 4 in the area under study.
Urumieh-Dokhtar magmatic belt (UDMB) is the main copper-bearing region in Iran. We have obtained the following results from this research.
9.1. Chemical Analysis Results
The chemical analysis of stream sediment samples in the area under study shows that:
・ A strong correlation exists between Pb and Ba, Pb and Ag, and Pb and As, but no statistically significant correlation is evident between Pb and Zn.
・ Factor analysis suggests that the element community can be related to parent rocks or mineralization of copper, zinc, and lead.
9.2. Remote Sensing Results
We mapped the distribution of hydrothermal alterations using ASTER data and obtained the following results in the area under study:
・ The main alteration types in the region are phyllic, argillic, propylitic, and iron oxide.
・ The development of hydrothermal alteration halos occurred around the faults and dikes.
・ The integration of the RBD, MNF, and LS-Fit algorithms create complete and accurate information regarding the distribution of hydrothermal alteration zones in the region.
9.3. Integration of the Layer Results
When several themes were combined using the AHP-OWA technique, the following results were achieved for the area under study:
・ Faults have played a significant role in mineralization in this area.
・ No mineralization has occurred in intrusive rocks.
・ Based on the final combined models generated by AHP-OWA, the “Average” model was considered for field checking.
9.4. Field Checking Results
Limited field checking was done to verify the validity of the results in this study. Based on field observation, the following results were obtained:
・ The results from remote sensing study have been verified by direct observation.
・ Copper and zinc mineralization have been discovered in the area under study.
・ Mineralization has not been found in the contact of dikes and country rocks.
・ Mineralization has not been found inside and contact of intrusion.
・ Field checking confirmed the accuracy of results from this study.
10. Conclusion
To produce the mineral exploration map for copper, zinc, and lead mineralization, we used fuzzy AHP and OWA technique on the basis of evidence layers derived from geological, geochemical, and remotely sensed data. Field anomaly checking demonstrates that the results of the fuzzy AHP-OWA model are reliable for the predicted areas of mineral potential.