Hybrid Fuzzy-Analytic Hierarchy Process (AHP) Model for Porphyry Copper Prospecting in Simorgh Area, Eastern Lut Block of Iran

The eastern Lut block of Iran has a high potential for porphyry copper mineralization due to the subduction tectonic regime. It is located in an inaccessible region and has harsh arid conditions for traditional mineral exploration campaigns. The objective of this study is to use Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) remote sensing data for porphyry copper exploration in Simorgh Area, eastern Lut block of Iran. Hydrothermal alteration zones such as argillic, phyllic and propylitic zones associated with porphyry copper systems in the study were identified using false color composition (FCC), band ratio (BR), principal component analysis (PCA) and minimal noise fraction (MNF). The thematic alteration layers extracted from FCC, BR, PCA and MNF were integrated using hybrid Fuzzy-AHP model to generate a porphyry copper potential map for the study area. Four high potential zones were identified in the central, western, eastern and northeastern of the study area. Fieldwork was used to validate the approach used in this study. This investigation exhibits that the use of hybrid Fuzzy-AHP model for the identification of hydrothermal alteration zones associated with porphyry copper systems that is typically applicable to ASTER data and can be used for porphyry copper potential mapping in many analogous metallogenic provinces.


Introduction
Remote sensing is a suitable tool to explore mineralization potential points [1] in the early stages of mineral exploration [2,3]. Hydrothermal fluids react with host rocks and change their mineralogical and chemical composition to form mineral deposits such as porphyry copper deposits [4,5]. These changes give rise to altered minerals in rocks while giving them unique spectral characteristics [4,6]. Geologists are used images from various satellites to detect the altered rocks associated with ore mineralization [7]. Integrating the results of image processing techniques can be used to increase the accuracy of the outputs and generate accurate potential map for reconnaissance stage of mineral exploration [3,8,9].
Analytic Hierarchy Process (AHP) is a way to help make decisions and emphasizes the importance of a decision maker as well as the value of comparing alternative options in the decision process [10]. This method is one of the popular methods in determining potential areas by combining information layers in mineral exploration [11,12]. One of its advantages is prioritizing each layer and weighting them. These weights can provide a map after applied to fuzzy layers. Given that the probability of mineralization in the whole region varies from 0 to 100, this method is called Fuzzy-AHP. Fuzzy maps show mineral anomalies based on the probability of mineralization. This means that in a region, the probability of mineralization is drawn from zero to one hundred percent. In fact, this is an advantage in mineral exploration, because the fuzzy anomalous areas are marked with a degree of probability. Due to the presence of porphyry copper ores, alterations and mineral indices, these layers can be used to make high-precision decisions [9].
Due to subduction tectonic regime documented in the eastern Lut block of Iran [13,14] (Figure 1a), the presence of porphyry copper systems can be anticipated. Numerous evidences have shown the occurrences of extensive alteration zones related to porphyry copper systems in this region [15], which unfortunately have not received much attention so far for porphyry copper exploration. Additionally, this zone is remote and has harsh arid conditions for field mineral exploration campaigns. Simorgh region is located between the two metallogenic provinces of Zahedan-Saravan and Birjand-Mirjaveh. In Simorgh region (Figure 1b), volcanic units include andesite and dacite and andesitic units are the most widespread in the region. Also, Hypabyssal intrusions interspersed with porphyry textures and penetrated volcanic rocks. This area has potential for copper mineralization as well as hydrothermal alteration zones related to porphyry copper mineralization can be identified and mapped using remote sensing techniques.
Mining 2022, 1, FOR PEER REVIEW 2 map after applied to fuzzy layers. Given that the probability of mineralization in the whole region varies from 0 to 100, this method is called Fuzzy-AHP. Fuzzy maps show mineral anomalies based on the probability of mineralization. This means that in a region, the probability of mineralization is drawn from zero to one hundred percent. In fact, this is an advantage in mineral exploration, because the fuzzy anomalous areas are marked with a degree of probability. Due to the presence of porphyry copper ores, alterations and mineral indices, these layers can be used to make high-precision decisions [9]. Due to subduction tectonic regime documented in the eastern Lut block of Iran [13,14] (Figure 1a), the presence of porphyry copper systems can be anticipated. Numerous evidences have shown the occurrences of extensive alteration zones related to porphyry copper systems in this region [15], which unfortunately have not received much attention so far for porphyry copper exploration. Additionally, this zone is remote and has harsh arid conditions for field mineral exploration campaigns. Simorgh region is located between the two metallogenic provinces of Zahedan-Saravan and Birjand-Mirjaveh. In Simorgh region (Figure 1b), volcanic units include andesite and dacite and andesitic units are the most widespread in the region. Also, Hypabyssal intrusions interspersed with porphyry textures and penetrated volcanic rocks. This area has potential for copper mineralization as well as hydrothermal alteration zones related to porphyry copper mineralization can be identified and mapped using remote sensing techniques. The purpose of this study is to use Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) remote sensing data for porphyry copper prospecting in Simorgh area, eastern Lut block of Iran. The specific objectives of this study are: (i) to identify hydrothermal alteration zones such as argillic, phyllic and propylitic zones associated with porphyry copper systems in the study area using false color composition (FCC), band ratio (BR), principal component analysis (PCA) and minimal noise fraction (MNF); (ii) to integrate the thematic layers extracted from FCC, BR, PCA and MNF using hybrid Fuzzy-AHP model; (iii) to generate porphyry copper potential map for the study area.

Geology of the Study Area
The Simorgh Area is approximately 40 square kilometers and located 110 km southwest of Nehbandan city and 24 km southwest of Dehsalm village in the southern part of South Khorasan province in east of Iran (Figure 1a,b). The Simorgh area is located in Dehsalm 1:100,000 geological map published by the Geological Survey of Iran (GSI), between the geographical coordinates of 59°11′00′′ to 59°15′56′′ East and 30°58′43′′ to 31°01′37.5″ north. It is placed in the central part of the Lut block and affected by igneous The purpose of this study is to use Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) remote sensing data for porphyry copper prospecting in Simorgh area, eastern Lut block of Iran. The specific objectives of this study are: (i) to identify hydrothermal alteration zones such as argillic, phyllic and propylitic zones associated with porphyry copper systems in the study area using false color composition (FCC), band ratio (BR), principal component analysis (PCA) and minimal noise fraction (MNF); (ii) to integrate the thematic layers extracted from FCC, BR, PCA and MNF using hybrid Fuzzy-AHP model; (iii) to generate porphyry copper potential map for the study area.

Geology of the Study Area
The Simorgh Area is approximately 40 square kilometers and located 110 km southwest of Nehbandan city and 24 km southwest of Dehsalm village in the southern part of South Khorasan province in east of Iran (Figure 1a,b). The Simorgh area is located in Dehsalm 1:100,000 geological map published by the Geological Survey of Iran (GSI), between the geographical coordinates of 59 • 11 00 to 59 • 15 56 East and 30 • 58 43 to 31 • 01 37.5 north. It is placed in the central part of the Lut block and affected by igneous and metamorphic processes. The study area is dry and warm with a very low annual average of precipitation. From the geomorphology point of view, the study area is characterized by several low altitude mountain hills separated by alluvial plains and dry water ways" According to the geological map, the rock units can be divided into five groups as follows [16]: (1) Pyroclastic and volcanic units including green and brown tuffs, volcanic breccia, hornblende andesite and hornblende quartz latite; (2) Semi-deep masses with porphyry texture, most of which are porphyry monzonite to porphyry diorite; (3) White dykes with a combination of porphyry granodiorite to porphyry syenite and general northeastsouthwest trend; (4) Rhyolite volcanic units with sulfide mineralization; (5) Single-source and multi-source hydrothermal breccia with an area of about 550 × 450 square meters in the central part of the outcrop (Figure 1b).

Remote Sensing Data Characteristics
ASTER was launched into orbit by the US National Aeronautics and Space Administration (NASA) in December 1999. The sensor is providing multispectral imaging of radiant and reflective electromagnetic energy from the Earth's surface and atmosphere in 14 bands. Spectral bands include: three bands in the visible and near infrared range (VNIR) with a spatial resolution of 15 m, 6 bands in the short-wave infrared (SWIR) with a spatial resolution of 30 m and finally 5 bands in the thermal infrared (TIR) range with a spatial resolution of 90 m [17] (Table 1). Due to their spectral capabilities and free accessibility, ASTER images, have found numerous applications in exploration of the mineral resources in last three decades [18,19]. Among them, studies on the ability of SWIR bands to identify hydrothermal alteration minerals such as alunite, kaolinite, calcite, dolomite, chlorite, and muscovite are noteworthy [20,21]. In addition, VNIR and TIR data can identify vegetation and iron oxides in topsoil and carbonates and silicates, respectively [22]. Given that the spatial resolution within the TIR range is 90 m, it is possible to identify rocks that are highly siliceous and have outcrops of 90 m or larger [23].

Preprocessing
In this study, the ASTER image (cloud free), level 1B, acquired on 15 May 2001, was used. After geo-referencing and geometric corrections of the image, a subset scene covering the study area was selected, cropped and pre-processed using crosstalk and Fast Line of sight Atmospheric Analysis of Spectral Hypercubes (FLAASH) methods. Crosstalk correction was performed to eliminate the effects of energy seepage from band 4 to bands 5 and 9, and FLAASH was performed to correct atmospheric effects [5].

Image Processing Techniques
All image processing tasks in this study have been done using ENVI ® Image Processing & Analysis Software version 5.1. (L3Harris Geospashtial, Boulder, Colorado, USA) After performing the necessary preprocessing steps, some widely used standard image processing methods were implemented including band composition, band ratio, principal component analysis (PCA) and minimal noise fraction (MNF). Due to interactions of different wavelengths with different rock units and also the higher sensitivity of the human eye to changes in color images compared to tone changes in black and white ones, different band combinations are used to highlight different phenomena. Assigning one of the three primary colors red, green, and blue to each wavelength creates a color combination, the Mining 2022, 2 4 most important of which is the false color composite (FCC) [24]. In the band ratio method, the ratio of same location pixels' values is calculated for two desired bands. By band ratio, some factors and parameters with adverse effects (such as topography) are removed, this method can be used to detect changes which cannot be recognized in single band images [3,25].
Principal component analysis (PCA) is a multivariate statistical method by which data are converted into new components that are least interdependent [26,27]. In this sense, spectral data is compressed into a new coordinate system, redundant data is removed, and finally new components with information from all spectral bands are created, and are used instead of the original data. This method has been used extensively by geoscientists to map lithologies as well as alterations [1,28,29]. In this study, VNIR + SWIR bands of ASTER were used. Table 2 shows statistical results derived from PCA analysis. The minimum noise fraction (MNF) method is used to determine the intrinsic dimensions of the image, to separate noise from the data, and help to reduce the amount of computational processing [30]. The MNF conversion consists of two steps: the first is noise removal, in which the principal components of the noise covariance matrix are calculated, and the second, in which the main components of noise-free data are extracted. Next, the results are divided into two groups, the group with high eigenvalues and the other group with eigenvalues close to one. The use of data with high eigenvalues causes noise separation from the data and improves spectral results. In this study, SWIR bands of ASTER were used to detect hydrothermal alterations from within igneous rocks. Table 3 shows statistical results derived from MNF analysis.

Fuzzy-AHP Model
In order to generate mineral potential map of the study area, information layers derived from FCC, BR, PCA and MNF were combined using Fuzzy-AHP model. The basis of this method is the weighting of each layer. In such a way that after fuzzification of layers, in the matrix, two-by-two weight comparisons are performed. And then using the AHP model, the maps are combined [31]. This model includes 4 steps as following: (1) Determining criteria and sub-criteria for using in the model.
(2) Determining the weight of criteria and sub-criteria using AHP model. (3) Fuzzification of information layers using fuzzy logic. (4) Integration of fuzzified layers, using the weights calculated in step 2.
All the layers used in this study were integrated by the Fuzzy-AHP method, as well as the steps of this method are presented in Figure 2. This figure shows the main criteria and methods used to obtain each of the information layers.
layers, in the matrix, two-by-two weight comparisons are performed. And then using the AHP model, the maps are combined [31]. This model includes 4 steps as following: (1) Determining criteria and sub-criteria for using in the model.
(2) Determining the weight of criteria and sub-criteria using AHP model.
(3) Fuzzification of information layers using fuzzy logic. (4) Integration of fuzzified layers, using the weights calculated in step 2.
All the layers used in this study were integrated by the Fuzzy-AHP method, as well as the steps of this method are presented in Figure 2. This figure shows the main criteria and methods used to obtain each of the information layers.

False Color Combination (FCC)
In order to detect propylitic and phyllic alteration zones in the study area, bands 4, 6 and 8 were used to produce Red (4), Green (6) and Blue (8) color combinations, respectively ( Figure 3). In the image-map, propylitic alteration zone (with high amounts of minerals containing Fe, Mg-OH) is observed in green color and phyllic alteration zone (with high amounts of minerals containing Al-OH) is observed in pink to yellow color. This is

False Color Combination (FCC)
In order to detect propylitic and phyllic alteration zones in the study area, bands 4, 6 and 8 were used to produce Red (4), Green (6) and Blue (8) color combinations, respectively ( Figure 3). In the image-map, propylitic alteration zone (with high amounts of minerals containing Fe, Mg-OH) is observed in green color and phyllic alteration zone (with high amounts of minerals containing Al-OH) is observed in pink to yellow color. This is due to the greater reflection of Al-OH minerals (alunite, kaolinite and muscovite) in the range of band 4 compared to bands 6 and 8 [29]. Many parts of recent alluvium are appeared as white color due to high amount of clay minerals (Figure 3). On the other hand, mafic to intermediate lithologies such as pyroxene diorite porphyry, syenite and diorite porphyry are manifested in green to black color due to high amount of mafic minerals (chlorite and epidote) in their content (Figure 3).
due to the greater reflection of Al-OH minerals (alunite, kaolinite and muscovite) in the range of band 4 compared to bands 6 and 8 [29]. Many parts of recent alluvium are appeared as white color due to high amount of clay minerals (Figure 3). On the other hand, mafic to intermediate lithologies such as pyroxene diorite porphyry, syenite and diorite porphyry are manifested in green to black color due to high amount of mafic minerals (chlorite and epidote) in their content (Figure 3). Band ratios of 5/6 and 8/9 have great capability to detect minerals with Al-Si-(OH) composition such as kaolinite (argillic zone), and Mg-Si-(OH) such as chlorite and epidote (propylitic zone), respectively [30]. Band ratio of 7/6 has high potential to identify muscovite-rich zone [31]. Muscovite mineralization can be considered as indicator of phyllic alteration zone, which is high potential zone of economical porphyry copper mineralization [32]. A RGB color combination of 7/6, 5/6 and 8/9 (R: 7/6, G: 5/6, B: 8/9) was developed to map phyllic, argillic and propylitic zones in the study area. The phyllic zone (with high economic potential) is represented in yellow to reddish yellow, the argillic zone is depicted as green and the propylitic zone is typically appeared in blue around other alteration zones (Figure 4). Band ratio technique identified the alteration zones specifically compared to FCC method. Propylitic alteration zone appears in green color, while phyllic alteration is depicted in pink to yellow.

Band Ratio (BR).
Band ratios of 5/6 and 8/9 have great capability to detect minerals with Al-Si-(OH) composition such as kaolinite (argillic zone), and Mg-Si-(OH) such as chlorite and epidote (propylitic zone), respectively [30]. Band ratio of 7/6 has high potential to identify muscovite-rich zone [31]. Muscovite mineralization can be considered as indicator of phyllic alteration zone, which is high potential zone of economical porphyry copper mineralization [32]. A RGB color combination of 7/6, 5/6 and 8/9 (R: 7/6, G: 5/6, B: 8/9) was developed to map phyllic, argillic and propylitic zones in the study area. The phyllic zone (with high economic potential) is represented in yellow to reddish yellow, the argillic zone is depicted as green and the propylitic zone is typically appeared in blue around other alteration zones (Figure 4). Band ratio technique identified the alteration zones specifically compared to FCC method.

Principal Component Analysis (PCA)
Principal component analysis was performed using a covariance matrix on nine bands of ASTER image covering the study area. The resulting eigenvector matrix is shown in Table 2. the PC that contains the target spectral information shows the highest eigen-

Principal Component Analysis (PCA)
Principal component analysis was performed using a covariance matrix on nine bands of ASTER image covering the study area. The resulting eigenvector matrix is shown in Table 2. the PC that contains the target spectral information shows the highest eigenvector loadings from the ASTER bands, coinciding with the target's most diagnostic features, but with opposite signs (+ or −). Component 1 (PC1) contains positive coefficients of all nine VNIR and SWIR bands and, due to its specific value, covers 91.6% of the total variance of the image. The overall brightness of the image creates a strong correlation among the bands. In the present analysis, this overall brightness is fully found in the first component (PC1). Vegetation has high reflectance and absorption in bands 3 and 2 of ASTER, respectively. The largest loading difference among the mentioned bands is observed in the PC4, and therefore the vegetation is found in PC4 (Figure 5a). The low percentage of variability of this component (0.94) indicates the dispersion of vegetation in the study area.

Minimum Noise Fraction (MNF)
MNF conversion was performed on SWIR bands of the ASTER image. The statistical results of this conversion are shown in Table 3. The eigenvalues of each converted MNF band represent the amount of information that the target band contains from the original image. So that the bands that contain the greatest noise have eigenvalues close to 1. Accordingly, bands with low eigenvalues will have very limited spatial integrity, which indicates the predominance of noise in those bands. By increasing the MNF principal number, the image quality decreases.
According to the Table 3, bands 4-6 are not suitable for performing RGB color combinations due to their eigenvalues close to 1. Therefore, only MNF bands with an eigenvalue percentage greater than 1 such as bands 1-3 were considered to produce the RGB color combination in this study. Figure 6 shows the RGB color combination of MNF bands 1-3 for the study area. According to this image, the areas that have undergone hydrothermal alteration are yellowish-brown in color, and most of the lithological boundaries are typically visible, which is in good agreement with the geological map of the area ( Figure  1). There is also a promising area in the central part of the image. Iron oxides have the highest absorption in bands 1 and 2 of ASTER, while they show the high reflection in band 4 [2]. In PC4, the eigenvectors related to bands 2 and 4 have moderate to high loadings with the opposite sign. The eigenvectors of bands 2, and 4 in the PC4 are −0.415 and 0.610, respectively. Iron oxides are shown as black pixels in Figure 5a.
Minerals containing Al (OH) (kaolinite, illite, muscovite and alunite) have very high absorption in band 6 of ASTER. Also, minerals containing Fe, Mg (OH) such as chlorite, epidote as well as carbonates such as calcite and dolomite have specific absorption in bands 8 of ASTER [1,3,33]. Considering the amount and sign of loadings obtained for bands 4 and 6 in the PC5, it can be concluded that this PC is suitable for showing argillic and phyllic alterations (Figure 5b). Also, considering the characteristics of bands 8 and 9 for the PC8, it can be said that this PC represents propylitic alteration (Figure 5c). Finally, RGB color combination of PC4, PC5 and PC8 (R: PC4, G: PC5, B: PC8) was created for a better display of the iron oxides and alterations in the study area. In Figure 5d, iron oxides, phyllic and argillic alterations are shown with light red, green to yellow, and blue, respectively. These alterations are well observed around the target areas.

Minimum Noise Fraction (MNF)
MNF conversion was performed on SWIR bands of the ASTER image. The statistical results of this conversion are shown in Table 3. The eigenvalues of each converted MNF band represent the amount of information that the target band contains from the original image. So that the bands that contain the greatest noise have eigenvalues close to 1. Accordingly, bands with low eigenvalues will have very limited spatial integrity, which indicates the predominance of noise in those bands. By increasing the MNF principal number, the image quality decreases.
According to the Table 3, bands 4-6 are not suitable for performing RGB color combinations due to their eigenvalues close to 1. Therefore, only MNF bands with an eigenvalue percentage greater than 1 such as bands 1-3 were considered to produce the RGB color combination in this study. Figure 6 shows the RGB color combination of MNF bands 1-3 for the study area. According to this image, the areas that have undergone hydrothermal alteration are yellowish-brown in color, and most of the lithological boundaries are typically visible, which is in good agreement with the geological map of the area (Figure 1). There is also a promising area in the central part of the image.

Integration of Geological and Alteration Layers
Using Fuzzy-AHP model, all of alteration layers resulted from FCC, BR, PCA and MNF were integrated. In order to increase the accuracy of the integration of information layers, a geological layer containing intrusive masses was also added. Then, according to the experts of mineral exploration, each of the criteria was weighed. The matrix of twoby-two criteria weighting is presented in Table 4. The final rank of the importance of each layer as well as the percentage of their impact in the combination is presented in Table 5. In order to compare binary and preferential criteria, which are: 1-Geology 2-PCA 3-MNF 4-FCC 5-BR, based on the opinion of mineral exploration experts as well as data processing, a comparison was made. The aim of this study was to investigate the alterations associated with porphyry copper mineralization [31][32][33][34].
Based on the available information layers, the lithology of the region was recognized as the main factor in creating and controlling hydrothermal alterations in the region. Therefore, in binary comparison, it is preferable to all information layers and criteria [1]. In the case of methods PCA, MNF, FCC and BR, which have all been done to identify hydrothermal alterations in the region by analyzing ASTER satellite images, a binary comparison was made based on the validity and accuracy of the methods [35][36][37][38][39][40]. Finally, based on the AHP method, the final rank of each of these criteria was calculated in Table  5.

Integration of Geological and Alteration Layers
Using Fuzzy-AHP model, all of alteration layers resulted from FCC, BR, PCA and MNF were integrated. In order to increase the accuracy of the integration of information layers, a geological layer containing intrusive masses was also added. Then, according to the experts of mineral exploration, each of the criteria was weighed. The matrix of two-by-two criteria weighting is presented in Table 4. The final rank of the importance of each layer as well as the percentage of their impact in the combination is presented in Table 5.  In order to compare binary and preferential criteria, which are: 1-Geology 2-PCA 3-MNF 4-FCC 5-BR, based on the opinion of mineral exploration experts as well as data processing, a comparison was made. The aim of this study was to investigate the alterations associated with porphyry copper mineralization [31][32][33][34].
Based on the available information layers, the lithology of the region was recognized as the main factor in creating and controlling hydrothermal alterations in the region. Therefore, in binary comparison, it is preferable to all information layers and criteria [1]. In the case of methods PCA, MNF, FCC and BR, which have all been done to identify hydrothermal alterations in the region by analyzing ASTER satellite images, a binary comparison was made based on the validity and accuracy of the methods [35][36][37][38][39][40]. Finally, based on the AHP method, the final rank of each of these criteria was calculated in Table 5.
As can be seen in Table 5, the layers of information, in order of importance high to low, are geology, PCA, MNF, BR and FCC respectively. Each of the information layers was fuzzified and then combined by the AHP method. Based on the results of weighting and determining the degree of importance of each information layer, a combined map of hydrothermal alterations and geological units affecting porphyry copper mineralization was presented in Figure 7. Four high potential zones were identified in the central, western, eastern and northeastern of the study area.  As can be seen in Table 5, the layers of information, in order of importance high to low, are geology, PCA, MNF, BR and FCC respectively. Each of the information layers was fuzzified and then combined by the AHP method. Based on the results of weighting and determining the degree of importance of each information layer, a combined map of hydrothermal alterations and geological units affecting porphyry copper mineralization was presented in Figure 7. Four high potential zones were identified in the central, western, eastern and northeastern of the study area.

Field Verification
Based on the potential map obtained from the Fuzzy-AHP model, three points of high potential zones were selected for ground control. Argillic alteration was observed at the first control point (Figure 8a). Iron oxide alteration was observed at the second control point (Figure 8b). At the third control point, there was a significant extent of Phyllic alteration (Figure 8c). The locations of all mentioned controlling points are shown in Figure 7. They have clear expansion at the surface with 500 to 1000 m length outcrop. According to the mineralogical studies conducted on samples taken from the target areas, illite/kaolinite, muscovite, goethite, and epidote/ carbonate were the prevalent minerals representing the argillic, phyllic, Fe-Oxide and propylitic alterations, respectively, concluding the compatibility of the field works with results obtained via remote sensing.

Field Verification
Based on the potential map obtained from the Fuzzy-AHP model, three points of high potential zones were selected for ground control. Argillic alteration was observed at the first control point (Figure 8a). Iron oxide alteration was observed at the second control point (Figure 8b). At the third control point, there was a significant extent of Phyllic alteration (Figure 8c). The locations of all mentioned controlling points are shown in Figure 7. They have clear expansion at the surface with 500 to 1000 m length outcrop. According to the mineralogical studies conducted on samples taken from the target areas, illite/kaolinite, muscovite, goethite, and epidote/ carbonate were the prevalent minerals representing the argillic, phyllic, Fe-Oxide and propylitic alterations, respectively, concluding the compatibility of the field works with results obtained via remote sensing.

Conclusions
In this study, a Hybrid Fuzzy-AHP model was developed for porphyry copper pr specting in Simorgh Area, Eastern Lut block of Iran. FCC, band ratio, PCA and MNF tec niques were used to identify alteration areas such as phyllic, argillic and propylitic zon associated with copper mineralization. Alteration thematic layers derived from the tec niques were weighted and fused to produce a potential map for the study area. Cons quently, a potential map for the study area was generated showing four main high pote tial zones in the central, western, eastern and northeastern sectors of the study area. T approach used in this investigation can be applied to other regions of Eastern Lut blo of Iran for porphyry copper prospecting.

Conclusions
In this study, a Hybrid Fuzzy-AHP model was developed for porphyry copper prospecting in Simorgh Area, Eastern Lut block of Iran. FCC, band ratio, PCA and MNF techniques were used to identify alteration areas such as phyllic, argillic and propylitic zones associated with copper mineralization. Alteration thematic layers derived from the techniques were weighted and fused to produce a potential map for the study area. Consequently, a potential map for the study area was generated showing four main high potential zones in the central, western, eastern and northeastern sectors of the study area. The approach used in this investigation can be applied to other regions of Eastern Lut block of Iran for porphyry copper prospecting.