Copper , Zinc , and Lead Mineral Prospectivity Mapping in the North of Tafresh , Markazi Province , Central Iran , Using the AHP-OWA Method

The reason of this research is to identify the favorable areas for copper, zinc, and lead mineralization in the western part of the 1:100,000 Tafresh geological Sheet in the Urmia-Dokhtar structural zone of Iran. Effective data layers for mineralization, such as geology, geochemistry, structures, and satellite images, were analyzed and then integrated using the AHP-OWA method to identify favorable areas. Geochemical stream samples were analyzed by univariate, multivariate, and classical statistical methods and revealed the first, second, and third class anomalies for copper, zinc, and lead in the study region. Detection of hydrothermal alteration zones by Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) satellite imagery in various algorithms, such as Relative absorption Band Depth (RBD), Minimum Noise Fraction (MNF), and Least Square Fit (LS-Fit), shows that argillic, phyllic, propylitic, and iron oxide alterations develop around the faults in the area under study. The favorable areas for copper, zinc, and lead mineralization have been identified by a combination of evidence maps of lithology, faults, dikes, geochemistry, and alteration data layers. Field observations in the area under study have confirmed the results.


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.
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).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].

Materials and Methods
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.

Geochemical Data
In this study, Pearson correlation confirmed a strong correlation between Pb 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).X S + is normally considered as the threshold value, 2 X S + as a possible anomaly, and 3 X S + as a probable anomaly.

Remote Sensing
Recent years, remote sensing images are used to identify hydrothermally altered 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.

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.

Relative Absorption Band Depth (RBD)
An RBD image is produced by dividing the sum of the maximum reflectance  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).

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).

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     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).

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.

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.

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).

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).

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.

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).

Pairwise Comparisons
In the AHP method, the pairwise comparison is the fundamental measurement mode [47].Table 4 shows the values from 1 to 9 that were used to rate the relative preferences for two elements of the hierarchy in the area under study.The pairwise comparison matrix and calculated weights for the objectives of the present study area were calculated using IDRISI software (Table 5).

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).

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).

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).

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).

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).

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.

Results
The area under study is located in Urmia-Dokhtar structural zone of Iran.The Urumieh-Dokhtar magmatic belt (UDMB) is the main copper-bearing region in Iran.We have obtained the following results from this research.

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.

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.

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.

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.

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 de- 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 .
Figure 1.The simple geology map of study area and its location in UDMB.

Figure 2 .
Figure 2. Locations of 520 stream sediment samples and their watersheds in the area under study.

Figure 3 .
Figure 3. Geochemical anomaly maps of Zn, Pb and Cu in the area under study.

Figure 8 .
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 .
Figure 9.The lithology layer of the area under study.

Figure 10 .
Figure 10.The dike layer of the area under study: Dike buffer map (left) and dike buffer reclassified map (right).

Figure 11 .
Figure 11.The alteration layer of the area under study.

Figure 12 .
Figure 12.(a) Pb geochemical layer.(b) Zn geochemical layer.(c) Cu geochemical layer.(d) The geochemical evidence layer of the area under study

Figure 13 .
Figure 13.Hierarchical diagram of exploration layers in the area under study.

Figure 15 .
Figure 15.Favorable area maps which have been generated by AHP-OWA procedures in the area under study.

Figure 16 .
Figure 16.Targets in the area under study.

Photo 2 .
Photo 2. The observations of the target 2 in the area under study

Photo 3 .
Photo 3. The observations of the target 3 in the area under study.

Table 1 .
Pearson correlation coefficients for chemical analyses of samples from the area under study.

Table 2 .
Varimax rotated factor matrix for 19 elements in 520 stream sediment samples of the area under study.

Table 3 .
Common ratio and band combinations for mapping alteration zones in the area under study.

Table 5 .
Pairwise comparison matrix of the level of objectives and calculated weights.