Application of Landsat-8, Sentinel-2, ASTER and WorldView-3 Spectral Imagery for Exploration of Carbonate-Hosted Pb-Zn Deposits in the Central Iranian Terrane (CIT)

The exploration of carbonate-hosted Pb-Zn mineralization is challenging due to the complex structural-geological settings and costly using geophysical and geochemical techniques. Hydrothermal alteration minerals and structural features are typically associated with this type of mineralization. Application of multi-sensor remote sensing satellite imagery as a fast and inexpensive tool for mapping alteration zones and lithological units associated with carbonate-hosted Pb-Zn deposits is worthwhile. Multiple sources of spectral data derived from different remote sensing sensors can be utilized for detailed mapping a variety of hydrothermal alteration minerals in the visible near infrared (VNIR) and the shortwave infrared (SWIR) regions. In this research, Landsat-8, Sentinel-2, Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and WorldView-3 (WV-3) satellite remote sensing sensors were used for prospecting Zn-Pb mineralization in the central part of the Kashmar–Kerman Tectonic Zone (KKTZ), the Central Iranian Terrane (CIT). The KKTZ has high potential for hosting Pb-Zn mineralization due to its specific geodynamic conditions (folded and thrust belt) and the occurrence of large carbonate platforms. For the processing of the satellite remote sensing datasets, band ratios and principal component analysis (PCA) techniques were adopted and implemented. Fuzzy logic modeling was applied to integrate the thematic layers produced by image processing techniques for generating mineral prospectivity maps of the study area. The spatial distribution of iron oxide/hydroxides, hydroxyl-bearing and carbonate minerals and dolomite were mapped using specialized band ratios and analyzing eigenvector loadings of the PC images. Subsequently, mineral prospectivity maps of the study area were generated by fusing the selected PC thematic layers using fuzzy logic modeling. The most favorable/prospective zones for hydrothermal ore mineralizations and carbonate-hosted Pb-Zn mineralization in the study region were particularly mapped and indicated. Confusion matrix, field reconnaissance and laboratory analysis were carried out to verify the occurrence of alteration zones and highly prospective locations of carbonate-hosted Pb-Zn mineralization in the study area. Results indicate that the spectral data derived from multi-sensor remote sensing satellite datasets can be broadly used for generating remote sensing-based prospectivity maps for exploration of carbonate-hosted Pb-Zn mineralization in many metallogenic provinces around the world.


Introduction
Remote sensing has provided tools for geological exploration for almost four decades. Nowadays, many satellite remote sensing datasets are accessible freely and can be extensively used for mineral exploration projects [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15]. Pb-Zn sulfide mineralization is typically associated with hydrothermal alteration zones, the contact boundaries of lithological units and structural features such as faults and fractures [5,[16][17][18][19]. According to World Bank Commodities Price Forecast (WBCPF), the price and consumption of Pb and Zn are increasing annually [20]. Pb and Zn are a necessity for the steady development of many countries around the world [21]. Accordingly, the exploration of Pb-Zn deposits using remote sensing satellite imagery as an available and inexpensive tool is of practical and economic interest.
Carbonate-hosted Pb-Zn deposits are some of the most significant sources of Pb and Zn [22]. The major hydrothermal alteration zones associated with carbonate-hosted Pb-Zn deposits are: (i) dissolution and hydrothermal brecciation, (ii) dolomite and calcite alteration, (iii) silicification and (iv) clay, mica, and feldspar diagenesis [22][23][24]. Besides, gossans as oxidation products of sulphide mineralized rocks are documented with carbonate-hosted Pb-Zn mineralization [25]. Particularly, studies on carbonate-hosted Pb-Zn deposits in Iran, India, China and Greenland have shown the possibility to identify hydrothermal alteration and iron oxides associated with Pb-Zn deposits using Landsat-8 and ASTER satellite imagery [5,[16][17][18][19]. The application of multi-sensor remote sensing satellite imagery and fusing the most informative alteration thematic layers using geostatistical models can provide a low-cost exploration approach for generating remote sensing-based prospectivity maps [15]. Multiple sources of spectral data derived from different remote sensing sensors can be utilized for detailed mapping a variety of hydrothermal alteration minerals in the VNIR and the SWIR regions.
Landsat-8 imagery contains nine bands (0.433 to 2.290 µm; 30 m spatial resolution) in the VNIR and SWIR regions ( Table 1). The VNIR spectral bands are particularly sensitive for mapping iron oxides/hydroxides, while SWIR spectral bands are responsive for detecting hydroxyl-bearing minerals and carbonates. These spectral bands have been broadly used for mapping hydrothermal alteration zones associated with hydrothermal ore mineralizations [3,5,12,15]. Sentinel-2 has thirteen spectral bands in the VNIR and the SWIR regions (0.433 to 2.280 µm; spatial resolutions from 10 to 60 m) ( Table 1) which are useful to identify iron oxides/hydroxides and hydroxyl-bearing minerals [6]. Six spectral bands in the SWIR range (1.600 to 2.430 µm; 30 m spatial resolution) allow the ASTER sensor to map clay and carbonate minerals (Table 1). Detailed detection and discrimination of hydroxyl-bearing minerals Figure 1. Geology map of the study area (modified from [36]).

Multi-Sensor Data Characteristics
In this investigation, multi-sensor satellite imagery, including Landsat-8, Sentinel-2, ASTER, and WV-3, was utilized for prospecting the carbonate-hosted Pb-Zn mineralization in the central part of the KKTZ region. Table 1 summarizes the technical performance and attributes of the Landsat-8, Sentinel-2, ASTER and WV-3 sensors. Broad spectral bands of Landsat-8 and Sentinel-2 were selected and used to discriminate lithological units and alteration zones for regional mapping of the study area. Narrow spectral bands of ASTER and WV-3 were used for detailed mapping of clay and carbonate minerals and Fe 2+ and Fe 3+ mineral groups.  19 100 TIRS2 (11) 11.50-12.51 Figure 1. Geology map of the study area (modified from [36]).

PrePprocessing of the Remote Sensing Datasets
The Fast Line-of-sight Atmospheric Analysis of Hypercubes (FLAASH) algorithm [54] were applied to Landsat-8 (Operational Land Imager (OLI) bands) by implementing the Mid-Latitude Summer (MLS) and Rural aerosol models [55]. Sentinel-2 data layer stacked of VNIR+SWIR bands (bands 2, 3, 4, 8, 11 and 12) with 10 m spatial dimension was generated to obtain a six-band dataset. The QUick Atmospheric Correction (QUAC) was performed on this dataset by using mud filtering to eliminate highly structured materials such as shallow water, mud and vegetation [55]. ASTER data layer stacked of VNIR+SWIR bands with 30-meter spatial dimensions was generated by using Pan-sharpening method [56]. Internal Average Relative Reflectance (IARR) calibration [57] was applied to Crosstalk corrected [58] ASTER data for atmospheric correction. The conversion to the Top-of-Atmosphere (TOA) spectral radiance and absolute radiometric correction are required for the WV-3 relative radiometrically corrected images [59]. Therefore, the corrections were executed and FLAASH algorithm was applied to the WV-3 data.

Image Processing Techniques
Band Ratios and Principal Components Analysis (PCA) image processing procedures were executed to extract key information related to alteration minerals and lithological units from the pre-processed remote sensing datasets. Successively, the most rational thematic layers of the alteration zones were fused using fuzzy logic modeling to generate mineral prospectivity maps of the study area. Finally, field reconnaissance, laboratory analysis and confusion matrix were carried out for verifying the remote sensing results. An overview of the methodological flowchart used in this study is displayed in Figure 2

Band Ratios
Band ratios method is broadly used for mapping hydrothermal alteration minerals and lithological units [2][3][4][60][61][62][63]. By ratioing bands that correspond to certain absorptions and reflectance, the pixels with particular mineral or mineral groups are highlighted [64][65][66][67]. Furthermore, this technique is proficient in reducing the topographic effects generated by slope orientations and solar illumination angles [68]. Several mathematical expressions were used for detecting alteration minerals or mineral groups, which are generally called Relative Absorption Band Depth (RBD) [69]. It includes three-point ratio formulation for revealing the mineral spectral intensities attributed to Fe 2+ , Fe 3+ , Fe-OH, Al-OH, Fe, Mg-OH, Si-OH, SO, CO 3 and SiO 2 [60,[70][71][72]. For a particular absorption or emissivity distinction, the numerator is the sum of the bands indicating the shoulders and the denominator is the band placed nearby the absorption or emissivity feature minimum [69]. To characterize iron oxide/hydroxides (hematite, goethite and jarosite), the VNIR spectral bands contain the most important information due to electronic transitions of Fe 3+ /Fe 2+ in the VNIR region from 0.45 to 1.2 µm [73,74] (Figure 3). In this study for detecting iron oxide at a regional scale, 4/2 band ratio of Landsat-8, Sentinel-2 and ASTER was selected to highlight iron oxide/hydroxides. Hydroxyl-bearing (Al-OH) alteration and carbonates (muscovite, kaolinite, gypsum, calcite and dolomite) show spectral absorption features in the 2.1-2.5 µm region due to overtones and combinations of the fundamental vibrations [75], whereas their spectral reflectance typically occur in 1.55-1.75 µm in the SWIR regions ( Figure 3). These characteristics are matched with band 7 (2.11-2.29 µm) and band 6 (1.57-1.65 µm) of Landsat-8, as well as band 12 (2.100-2.280 µm) and band 11 (1.565-1.655 µm) of Sentinel-2, respectively ( Figure 3). Therefore, the 6/7 band ratio of Landsat-8 and 11/12 band ratio of Sentinel-2 were used to map hydroxyl bearing alteration minerals and carbonates in this study at the regional scale.

Principal Components Analysis (PCA)
PCA is a mathematical technique that transforms a quantity of correlated variables into a number of uncorrelated linear variables called principal components (PCs) [83][84][85]. The PCA is usually implemented on a square symmetric matrix. It can be based on the covariance matrix (scaled sums of squares and cross products) or correlation matrix (sums of squares and cross products from standardized data) [83,86,87]. The PCA is broadly used to map hydrothermal alteration minerals and lithological units using spectral bands of remote sensing sensors [8,61,62,84,[88][89][90][91][92]. The uncorrelated linear combinations (eigenvector loadings) contain indicative information allied to spectral characteristics of alteration minerals that can be expected from the specific spectral bands in the VNIR and SWIR regions [84,89]. Accordingly, a PC contains strong eigenvector loadings for indicative bands such as reflective and absorptive bands of an alteration mineral or mineral group with opposite signs enhances that mineral or mineral group as bright or dark pixels in the PC image. Positive loading in a reflective band enhances the alteration mineral as bright pixels, while negative loading is in a reflective band depicts the alteration mineral as dark pixels [84,89]. In this analysis, the PCA method was implemented based on covariance matrix to the selected bands of Landsat-8 (OLI bands), Sentinel-2 (bands 2, 3, 4, 8, 11 and 12) and ASTER (VNIR+SWIR bands) and WV-3 (VNIR bands) for identifying hydrothermal alteration mineral assemblages in the study area. Tables 2-5 show the eigenvector matrix for the selected bands of the remote sensing datasets.

Fuzzy Logic Modeling
Fuzzy logic modeling is based on the fuzzy set theory, which was proposed by Zadeh [93]. It is a form of many-valued logic in which the true values of variables may be any real number between 0 and 1 both inclusive [94]. A fuzzy set of A is a set of ordered pairs: where µ_A (x) is termed the membership function or membership grade of x in A. µ_A (x) maps x to membership space (M), when M contains only the two points 0 and 1. The range of µ_A (x) is [0, 1], where zero expresses non-membership and one expresses full membership [93]. Fuzzy logic modeling has been successfully applied for mineral prospectivity mapping in metallogenic provinces [29,[95][96][97][98].
The application of fuzzy logic modeling for mineral prospectivity mapping normally includes three main feed-forward stages: (i) fuzzification of evidential data; (ii) logical integration of fuzzy evidential maps with the aid of an inference network and appropriate fuzzy set operations; and (iii) defuzzification of fuzzy mineral prospectivity output to aid its interpretation [96]. A set of fuzzy membership values is expressed in a continuous series from 0 to 1. This 0-1 scale, however, does not constitute a probability density function. Function-member values are established for each evidence map that will be integrated.
In the fuzzy logic method, a total of sheet maps (fuzzy membership) based on the significance distance of features are weighted (for each pixel or spatial position, a particular weight between 0 and 1 is appointed) [96,99]. Five operators that are useful for combining mineral exploration datasets include the fuzzy AND, fuzzy OR, fuzzy algebraic product, fuzzy algebraic sum and fuzzy gamma [95,[100][101][102]. In this analysis, the multiclass evidential image-maps were reclassified in 10 classes of equal interval and then were fuzzified using the linear membership function. The fuzzy gamma operator was used for mapping the prospective areas using Landsat-8, Sentinel-2, ASTER and WV-3 alteration thematic input layers. After testing several values for the γ parameter, it was adjusted 0.70 for ensuring a flexible compromise between the fuzzy algebraic sum and the fuzzy algebraic product [99]. Table 6 shows the fuzzification parameters for the input layers used in this analysis.

Fieldwork Data and Laboratory Analysis
Geological survey and laboratory analysis were carried out to confirm the image processing results and mineral prospectivity mapping for the central part of the KKTZ region. A global positioning system (GPS) survey was conducted in the study area for verifying the spatial distribution of alteration zones and lithological units using a handheld Monterra GPS (average accuracy of 3 m; Garmin, New Taipei City, Taiwan). Additionally, numerous photos were taken from alteration zones and lithological units during the field surveys. Thirty hand specimens were collected from alteration zones, ore mineralization and lithological units for laboratory analysis. Polish sections of ore mineralization and thin sections of alteration zones and lithological units were prepared. For a detailed mineralogical study of alteration zones, X-ray diffraction (XRD) analysis was implemented using an X'pert Pro XRD diffractometer (Philips, Amsterdam, The Netherlands) located at the Iran Mineral Processing Research Center (IMPRC, Tehran, Iran). Moreover, analytical spectral devices (ASD) spectroscopy was performed to the samples collected from the main lithological units exposed in the study area using a FieldSpec3 ® spectroradiometer (Malvern Panalytical Ltd., Malvern, UK, operating from 0.35 µm to 2.5 µm) located at the University of Kerman Institute of Science and High Technology (Kerman, Iran).

Lithological and Alteration Mapping Using Landsat-8, Sentinel-2 and ASTER
For generating a regional view of lithological units in the study region, the Red-Green-Blue (RGB) false color composite of bands 2, 5 and 7 for Landsat-8 and bands 2, 8 and 12 for Sentinel-2 were considered, respectively. The resultant images show most of the lithological units having spectral features related to Fe 3+ and Fe 3+ /Fe 2+ iron oxides and clay and carbonate minerals. Regarding the geological map of the study area (see Figure 1), the identification and lithological discrimination of the units in Landsat-8 and Sentinel-2 resultant images were almost similar. Therefore, the RGB false color composite image of the Sentinel-2 spectral bands was selected and presented herein ( Figure 4).
The lithological units such as the rhyolite and acidic tuff interbedded with crystalline limestone (PE r ) unit, alternation of limestone, dolomite and shale (the Shemshak, Espahk and Surkh Formations) (TR), the alternation of light green sandstone, slate and phyllite (C ph.s ) and the Zarigan granite (Zrg) show a mixture of Fe 3+ iron oxides and clay and carbonates (brown, blue and purple shades). The marble, garnet schist, biotite schist and quartzite (PE 1.s ) unit, light brown thin bedded dolomite and limestone, red shale and siltstone and white to red quartzite (Kuhbanan Formation) (E 1 ), marl interbedded with gypsum and conglomerate (Ngc) unit and the massive to thick bedded limestone (K 1 ) unit mostly contain Fe 3+ /Fe 2+ iron oxides mixed with clay and carbonates (light green and cyan hues). The sphilite, basalt, tuff and grey to light brown limestone (PE v ) unit is characterized as dark brown color due to high content of Fe 3+ iron oxides. The northeastern part of the scene (yellow polygon) shows a variety of colors related to absorption features of target alteration minerals including Fe 3+ and Fe 3+ /Fe 2+ iron oxides and clay and carbonate minerals. The massive to medium bedded limestone and dolomite (C 1d ) unit, gypsum interbedded with grey limestone (DCg), the C ph.s unit and the PE r unit are the main lithological units in this zone. A high level of Fe 3+ iron oxides (orange, pink to rose blush shades), mixture of Fe 3+ and Fe 3+ /Fe 2+ iron oxides (yellow tone) and clay and carbonate minerals (light blue hue) in the zone is notable (Figure 4). Note that the occurrence of carbonate-hosted Pb-Zn mineralization is reported in this zone (yellow polygon) and detailed alteration mapping is presented by authors during this study.
Remote Sens. 2020, 12, x FOR PEER REVIEW 12 of 35 false color composite image of the Sentinel-2 spectral bands was selected and presented herein ( Figure 4). The lithological units such as the rhyolite and acidic tuff interbedded with crystalline limestone (PE r ) unit, alternation of limestone, dolomite and shale (the Shemshak, Espahk and Surkh Formations) (TR), the alternation of light green sandstone, slate and phyllite (C ph.s ) and the Zarigan granite (Zrg) show a mixture of Fe 3+ iron oxides and clay and carbonates (brown, blue and purple shades). The marble, garnet schist, biotite schist and quartzite (PE 1.s ) unit, light brown thin bedded dolomite and limestone, red shale and siltstone and white to red quartzite (Kuhbanan Formation) (E 1 ), marl interbedded with gypsum and conglomerate (Ngc) unit and the massive to thick bedded limestone (K 1 ) unit mostly contain Fe 3+ /Fe 2+ iron oxides mixed with clay and carbonates (light green and cyan hues). The sphilite, basalt, tuff and grey to light brown limestone (PE v ) unit is characterized as dark brown color due to high content of Fe 3+ iron oxides. The northeastern part of the scene (yellow polygon) shows a variety of colors related to absorption features of target alteration minerals including Fe 3+ and Fe 3+ /Fe 2+ iron oxides and clay and carbonate minerals. The massive to medium bedded limestone and dolomite (C 1d ) unit, gypsum interbedded with grey limestone (DCg), the C ph.s unit and the PE r unit are the main lithological units in this zone. A high level of Fe 3+ iron oxides (orange, pink to rose blush shades), mixture of Fe 3+ and Fe 3+ /Fe 2+ iron oxides (yellow tone) and clay and carbonate minerals (light blue hue) in the zone is notable (Figure 4). Note that the occurrence of carbonate-hosted Pb-Zn mineralization is reported in this zone (yellow polygon) and detailed alteration mapping is presented by authors during this study.  The RGB false color composite of bands 6, 2 and 8 of ASTER discriminate most of the lithological units containing Al/Fe-OH, Fe +2 /Fe 3+ and Mg-Fe-OH/CO3 absorption properties ( Figure 5). The PE r unit, the Zarigan granite (Zrg), thin bedded limestone and dolomite (PE Ed ) unit and some parts of the Kuhbanan Formation (E 1 ) appear in purple due to Al/Fe-OH minerals that is slightly mixed with Mg-Fe-OH/CO 3 minerals. The DCg unit is depicted as a green tone because of Fe +2 /Fe 3+ minerals. The C ph.s unit is characterized in cyan hue attributed to mixture of Fe +2 /Fe 3+ and Mg-Fe-OH/CO 3 minerals. The Ngc unit and the K 1 unit are represented as bright yellow shade probably because of combination between Mg-Fe-OH/CO 3 , Al/Fe-OH and Fe +2 /Fe 3+ minerals. The PE v unit is considered as dark shade due to high content of Fe 3+ iron oxides. The TR unit having Al/Fe-OH and Mg-Fe-OH/CO 3 minerals is presented in a dark purple tone. The C 1d unit is manifested in brown-golden color due to Mg-Fe-OH/CO 3 minerals mixed with Fe +2 /Fe 3+ minerals. The E unit, PE r and some parts of the Kuhbanan Formation (E 1 ) appear as dark green hue because of Fe +2 /Fe 3+ minerals ( Figure 5).
Remote Sens. 2020, 12, x FOR PEER REVIEW 13 of 35 The RGB false color composite of bands 6, 2 and 8 of ASTER discriminate most of the lithological units containing Al/Fe-OH, Fe +2 /Fe 3+ and Mg-Fe-OH/CO3 absorption properties ( Figure  5). The PE r unit, the Zarigan granite (Zrg), thin bedded limestone and dolomite (PE Ed ) unit and some parts of the Kuhbanan Formation (E 1 ) appear in purple due to Al/Fe-OH minerals that is slightly mixed with Mg-Fe-OH/CO3 minerals. The DCg unit is depicted as a green tone because of Fe +2 /Fe 3+ minerals. The C ph.s unit is characterized in cyan hue attributed to mixture of Fe +2 /Fe 3+ and Mg-Fe-OH/CO3 minerals. The Ngc unit and the K 1 unit are represented as bright yellow shade probably because of combination between Mg-Fe-OH/CO3, Al/Fe-OH and Fe +2 /Fe 3+ minerals. The PE v unit is considered as dark shade due to high content of Fe 3+ iron oxides. The TR unit having Al/Fe-OH and Mg-Fe-OH/CO3 minerals is presented in a dark purple tone. The C 1d unit is manifested in brown-golden color due to Mg-Fe-OH/CO3 minerals mixed with Fe +2 /Fe 3+ minerals. The E unit, PE r and some parts of the Kuhbanan Formation (E 1 ) appear as dark green hue because of Fe +2 /Fe 3+ minerals ( Figure 5). show the results of band ratios for mapping iron oxide/hydroxides and clay and carbonate minerals derived from Landsat-8 and Sentinel-2 spectral bands. The 4/2 band ratio of Landsat-8 and Sentinel-2 shows the spatial distribution of iron oxide/hydroxide minerals (red pixels), which are mostly mapped in the PE 1.s unit, the Ngc unit, the K 1 unit, the PE r unit, the DCg unit and Quaternary deposits (Q t ) ( Figures 6A and B). The 6/7 band ratio of Landsat-8 and 11/12 band ratio of Sentinel-2 map the surface distribution of hydroxyl-bearing alteration minerals and  Figure 6A-D shows the results of band ratios for mapping iron oxide/hydroxides and clay and carbonate minerals derived from Landsat-8 and Sentinel-2 spectral bands. The 4/2 band ratio of Landsat-8 and Sentinel-2 shows the spatial distribution of iron oxide/hydroxide minerals (red pixels), which are mostly mapped in the PE 1.s unit, the Ngc unit, the K 1 unit, the PE r unit, the DCg unit and Quaternary deposits (Q t ) ( Figure 6A,B). The 6/7 band ratio of Landsat-8 and 11/12 band ratio of Sentinel-2 map the surface distribution of hydroxyl-bearing alteration minerals and carbonates (green pixels) ( Figure 6C,D). Accordingly, the PE 1.s unit, tuff and thick bedded of red silty shale (E) unit, the K 1 unit, the PE r unit, the pink to white granite (Zarigan granite) (Zrg), the C 1d unit and the Q t deposits are mapped due to high content of clay and carbonate minerals in their composition. Clay and carbonate minerals show higher surface abundance in the Landsat-8 ratio-image ( Figure 6C) compared to the Sentinel-2 ratio-image ( Figure 6D). The K 1 unit is particularly mapped in the Landsat-8 ratio-image ( Figure 6C), while it shows the smaller spatial distribution in the Sentinel-2 ratio-image ( Figure 6D). It is because the band 12 (2.100-2.280 µm) of Sentinel-2 mostly covers the absorption features of hydroxyl-bearing minerals (2.10-2.20 µm), while the absorption features of carbonates typically concentrate around 2.350 to 2.450 µm [69,82]. Hence, carbonates could not be mapped using 11/12 band ratio of Sentinel-2, properly.
Remote Sens. 2020, 12, x FOR PEER REVIEW 14 of 35 carbonates (green pixels) ( Figures 6C and D). Accordingly, the PE 1.s unit, tuff and thick bedded of red silty shale (E) unit, the K 1 unit, the PE r unit, the pink to white granite (Zarigan granite) (Zrg), the C 1d unit and the Q t deposits are mapped due to high content of clay and carbonate minerals in their composition. Clay and carbonate minerals show higher surface abundance in the Landsat-8 ratio-image ( Figure 6C) compared to the Sentinel-2 ratio-image ( Figure 6D). The K 1 unit is particularly mapped in the Landsat-8 ratio-image ( Figure 6C), while it shows the smaller spatial distribution in the Sentinel-2 ratio-image ( Figure 6D). It is because the band 12 (2.100-2.280 µm) of Sentinel-2 mostly covers the absorption features of hydroxyl-bearing minerals (2.10-2.20 µm), while the absorption features of carbonates typically concentrate around 2.350 to 2.450 µm [69,82]. Hence, carbonates could not be mapped using 11/12 band ratio of Sentinel-2, properly. The 4/2 band ratio of ASTER shows the iron oxide/hydroxide minerals in red pixels, which are typically mapped in the PE 1.s unit, the K 1 unit, the PE r unit, the DCg unit, the PE Ed unit, the Kuhbanan The 4/2 band ratio of ASTER shows the iron oxide/hydroxide minerals in red pixels, which are typically mapped in the PE 1.s unit, the K 1 unit, the PE r unit, the DCg unit, the PE Ed unit, the Kuhbanan Formation (E 1 ) and the C 1d unit ( Figure 7A). The Al/Fe-OH minerals are specifically mapped as mustard pixels in the RBD1 image of ASTER ( Figure 7B). Comparison to the geology map of the study area, the high abundance of Al/Fe-OH minerals (mustard pixels) are associated with the lithological units of the C ph.s , the PE r , the Kuhbanan Formation (E 1 ), the PE v unit, the PE Ed unit, the TR and Quaternary deposits (Qt) ( Figure 7B). The RBD2 image of ASTER shows the surface abundance of Mg-Fe-OH/CO3 minerals (green pixels; Figure 7C), which are typically associated with the K 1 unit and the C 1d unit. Although the Kuhbanan Formation (E 1 ), the PE 1.s unit, the Ngc unit and the PE r unit are also partially mapped in Figure 7C. Dolomitic units are specifically mapped in the RBD3 image of ASTER as yellow pixels ( Figure 7D). The C 1d lithological unit is strongly represented in yellow pixels. Additionally, some parts of the Kuhbanan Formation (E 1 ), the TR unit, PE Ed unit, the PE 1.s unit and the PE r unit are mapped in Figure 7D due to the high content of dolomite in their lithological composition. Formation (E 1 ) and the C 1d unit ( Figure 7A). The Al/Fe-OH minerals are specifically mapped as mustard pixels in the RBD1 image of ASTER ( Figure 7B). Comparison to the geology map of the study area, the high abundance of Al/Fe-OH minerals (mustard pixels) are associated with the lithological units of the C ph.s , the PE r , the Kuhbanan Formation (E 1 ), the PE v unit, the PE Ed unit, the TR and Quaternary deposits (Qt) ( Figure 7B). The RBD2 image of ASTER shows the surface abundance of Mg-Fe-OH/CO3 minerals (green pixels; Figure 7C), which are typically associated with the K 1 unit and the C 1d unit. Although the Kuhbanan Formation (E 1 ), the PE 1.s unit, the Ngc unit and the PE r unit are also partially mapped in Figure 7C. Dolomitic units are specifically mapped in the RBD3 image of ASTER as yellow pixels ( Figure 7D). The C 1d lithological unit is strongly represented in yellow pixels. Additionally, some parts of the Kuhbanan Formation (E 1 ), the TR unit, PE Ed unit, the PE 1.s unit and the PE r unit are mapped in Figure 7D due to the high content of dolomite in their lithological composition.  Analyzing the eigenvector loadings derived from the PCA method for Landsat-8 and Sentinel-2 indicates that the PC4 and PC5 contain key information related to the hydroxyl-bearing and carbonate and iron oxide/hydroxide minerals. Table 2 shows the eigenvector matrix for Landsat-8 selected bands. The PC4 has strong loadings of bands 6 (0.667216) and 7 (−0.48894) with opposite signs. Thus, the PC4 image identify hydroxyl-bearing minerals and carbonates as bright pixels due to a positive sign in the reflection band (band 6). Figure 8A) shows the PC4 image of Landsat-8. The spatial distribution of hydroxyl-bearing minerals and carbonates (green pixels) is mostly identified with the units of the PE 1.s , the PE r , the K 1 , the C ph.s , the DCg, the C 1d and Q t deposits. The Zarigan granite (Zrg) and Kuhbanan Formation (E 1 ) also show some hydroxyl-bearing minerals in the PC4 image of Landsat-8 ( Figure 8A). The PC5 of Landsat-8 has moderate positive contribution in band 2 (0.318038), negative strong loadings in band 4 (−0.579504) and strong positive contribution (0.516661) in band 5 (Table 2). Accordingly, iron oxide/hydroxide (Fe 3+ /Fe 2+ ) minerals manifest as dark pixels in the PC5 image because of negative contribution of band 4 (reflection band). The dark pixels were inverted to bright pixels by negation (multiplication to −1). Figure 8B displays the spatial distribution of iron oxide/hydroxides (red pixels) in the study area. The units of the PE 1.s , the C ph.s , the Ngc, the PE r , the Zarigan granite (Zrg), the TR and the DCg show a high abundance of oxide/hydroxides in the PC5 image of Landsat-8 ( Figure 8B). Comparison of the Landsat-8 PCA results to the band ratio indicates that the spatial distribution of iron oxide/hydroxides minerals and hydroxyl-bearing minerals and carbonates in the PC images is generally less widespread. Moreover, the alteration zones mapped by the PCA method show better spatial relationship with the documented ore mineral occurrences in the study area. The PC4 has strong positive loadings in band 2 (0.812776), while it shows strong negative loadings in band 4 (−0.429509) (see Table 4). Iron oxide/hydroxides minerals characterize by high absorption features about 0.40 to 1.10 µm and high reflection around 1.60 µm [75]. Considering the spectral location of bands 2 (0.63-0.69 µm) and 4 (1.60-1.70 µm) of ASTER and the eigenvector loadings in PC4, it is discernable that iron oxide/hydroxides minerals depict as dark pixels in the PC4 image, which consequently negated to bright pixels and overlain by red color ( Figure 9B). The Ngc unit, the Kuhbanan Formation (E 1 ) and the PE Ed unit are strongly mapped in the image. However, some zones located in the K 1 unit, the PE 1.S unit, the PE r unit and Qt deposits are also distinguishable ( Figure 9B). The PC5 shows a strong contribution in band 4 (-0.668838) and band 6 (0.425904) and band 9 (0.507113) with opposite signs (see Table 4). The OH/S-O/CO3-bearing minerals exhibit diagnostic absorption features at 2.20 to 2.50 µm, which are coincident with bands 6 to 9 of ASTER [74]. Therefore, OH/S-O/CO3-bearing minerals can be mapped as dark pixels in the PC5 image. Figure 9C shows the PC5 image that dark pixels are converted to bright and overlain by mustard color. Spatial distribution of OH/S-O/CO3-bearing minerals is clearly observable in the C ph.s unit, the DCg unit, the PE Ed unit the TR unit the PE r unit and Qt deposits. Some small parts of the Kuhbanan Formation (E 1 ) and the PE 1.S unit are also appeared as mustard pixels ( Figure 9C). The results of PC5 image is almost matched to the RBD1 image-map (see Figure 7). Considering eigenvector loadings for selected bands of Sentinel-2 (see Table 3), the PCA4 includes a great contribution of band 11 (−0.641949) with a negative sign and band 12 (0.567634) with a positive sign. Hence, hydroxyl-bearing minerals represent in dark pixels of the PC4 image ( Figure 8C). The dark pixels were converted to bright pixels by negation, subsequently. Results show that the spatial distribution of hydroxyl-bearing minerals is typically associated with the units of the PE 1.s , the C ph.s , the Kuhbanan Formation (E 1 ), the PE r , the K 1 , the Zarigan granite (Zrg), the Ngc, the C 1d and Q t deposits ( Figure 8C). Iron oxide/hydroxide (Fe 3+ /Fe 2+ ) minerals appear as bright pixels in the PCA5 image of Sentinel-2 due to a strong positive loading in band 4 (0.698634) and negative loadings of band 2 (−0.454829) and band 8 (−0.523627), respectively (see Table 3). Iron oxide/hydroxide minerals are represented in bright pixels as a result of a positive loading in the reflection band (band 4) ( Figure 8D). The high surface abundance of iron oxide/hydroxide (red pixels) is highlighted with the PE 1.s unit, the C ph.s unit, the Ngc unit, the PE r unit, the DCg unit, the Zarigan granite (Zrg), the TR unit and Q t deposits ( Figure 8D). The spatial distribution of iron oxide/hydroxides minerals in the PCA5 image of Sentinel-2 is less prevalent compared to the PCA5 of Landsat-8 (see Figure 8B,D). Nevertheless, the PC4 image of Landsat-8 and Sentinel-2 are almost identical in many parts (see Figure 8A,C).
The PCA technique was also implemented on ASTER VNIR+SWIR bands for mapping the target alteration minerals. The eigenvector matrix for ASTER data is shown in Table 4. The PC3 contains strong positive loading in band 4 (0.564435) and strong negative loading in band 8 (−0.896111). The Mg-Fe-OH/CO3 has high reflectance about 1.6 µm that is coincident with band 4 of ASTER. The Fe-Mg-OH and CO3 minerals exhibit diagnostic absorption features near 2.350 µm, which is matched with bands 8 of ASTER [60,75]. So, the PC3 image maps Mg-Fe-OH/CO3 minerals as bright pixels because of positive loading in band 4 (reflection band). Figure 9A shows the PC3 image-map of the study area that overlain by green color. The green pixels show a high concentration of Mg-Fe-OH/CO 3 minerals. Referring to the local geology map of the study zone, the Mg-Fe-OH/CO 3 minerals are typically identified in the K 1 unit, the Ngc unit and the C 1d . The Kuhbanan Formation (E 1 ), the PE r unit, the Zarigan granite (Zrg) and the PE 1.S unit are also weakly mapped in the PC3 image ( Figure 9A). These lithological units contain a high abundance of carbonates and Mg-Fe-OH minerals, which can be detected by band 8 (2.295-2.365 µm) of ASTER. The detected pixels in the PC3 image-map are comparable to the RBD2 image-map (see Figure 7C).
The PC4 has strong positive loadings in band 2 (0.812776), while it shows strong negative loadings in band 4 (−0.429509) (see Table 4). Iron oxide/hydroxides minerals characterize by high absorption features about 0.40 to 1.10 µm and high reflection around 1.60 µm [75]. Considering the spectral location of bands 2 (0.63-0.69 µm) and 4 (1.60-1.70 µm) of ASTER and the eigenvector loadings in PC4, it is discernable that iron oxide/hydroxides minerals depict as dark pixels in the PC4 image, which consequently negated to bright pixels and overlain by red color ( Figure 9B). The Ngc unit, the Kuhbanan Formation (E 1 ) and the PE Ed unit are strongly mapped in the image. However, some zones located in the K 1 unit, the PE 1.S unit, the PE r unit and Qt deposits are also distinguishable ( Figure 9B). The PC5 shows a strong contribution in band 4 (-0.668838) and band 6 (0.425904) and band 9 (0.507113) with opposite signs (see Table 4). The OH/S-O/CO3-bearing minerals exhibit diagnostic absorption features at 2.20 to 2.50 µm, which are coincident with bands 6 to 9 of ASTER [74]. Therefore, OH/S-O/CO3-bearing minerals can be mapped as dark pixels in the PC5 image. Figure 9C shows the PC5 image that dark pixels are converted to bright and overlain by mustard color. Spatial distribution of OH/S-O/CO3-bearing minerals is clearly observable in the C ph.s unit, the DCg unit, the PE Ed unit the TR unit the PE r unit and Qt deposits. Some small parts of the Kuhbanan Formation (E 1 ) and the PE 1.S unit are also appeared as mustard pixels ( Figure 9C). The results of PC5 image is almost matched to the RBD1 image-map (see Figure 7).  Figure 10 shows the WV-3 scene covering the selected subset of the study area contains Zn-Pb mineralization. The band ratio image-map of 5/3 for mapping Fe 3+ oxides (A), the band ratio image-map of 5/2 for identifying iron-stained alteration (B), the band ratio image-map of 7/4 for detecting dolomite/Fe 2+ (C), and RGB false color composite image-map of 7/4, 5/3 and 5/2 for discriminating lithological units (D) are shown in Figure 10. The surface distribution of Fe 3+ oxides (bright pixels) is generally associated with the DCg unit, the PE r unit and PE Ed unit ( Figure 10A). Checking the PC6 and PC7 images indicate that they are mostly noisy and uninformative. Considering the laboratory reflectance spectra of dolomite resampled to response functions of VINR+SWIR bands of ASTER (see Figure 4) reveals that the PC8 image feasibly contains vital information related to the spatial distribution of dolomite. Bands 4 and 5 can be considered reflection bands, whereas band 7 can be deliberated absorption bands for detecting dolomite. The PCA8 has strong positive loadings in band 4 (0.404433) and band 5 (0.558407) and great contribution of band 7 (−0.623420) with a negative sign (see Table 4). The PC8 image-map of the study zone is shown in Figure 9D. Regarding the geology map of the study zone, the C 1d unit clearly appears as bright pixels in the PC8 image. Moreover, the Kuhbanan Formation (E 1 ), the TR, the PE r unit, the PE 1.S unit and the PE Ed unit are partially characterized in yellow pixels. The resultant image-map of the PC8 is similar to the RBD3. Figure 10 shows the WV-3 scene covering the selected subset of the study area contains Zn-Pb mineralization. The band ratio image-map of 5/3 for mapping Fe 3+ oxides (A), the band ratio image-map of 5/2 for identifying iron-stained alteration (B), the band ratio image-map of 7/4 for detecting dolomite/Fe 2+ (C), and RGB false color composite image-map of 7/4, 5/3 and 5/2 for discriminating lithological units (D) are shown in Figure 10. The surface distribution of Fe 3+ oxides (bright pixels) is generally associated with the DCg unit, the PE r unit and PE Ed unit ( Figure 10A). The band ratio image-map of 5/2 generates analogous results for mapping all iron oxides (iron-stained alteration) associated with the above-mentioned lithological units (Figure 10 B). This similarity of iron oxide mapping results might be related to extensive iron-stained occurrences in the lithological units. The spatial distribution of dolomite/Fe 2+ (bright pixels) is typically mapped in the eastern part of the selected subset (Figure 10 C). The PE r unit, PE Ed unit and the C 1d unit are highlighted in the image-map of the 7/4 band ratio ( Figure 10C). Figure 10D shows RGB false color The band ratio image-map of 5/2 generates analogous results for mapping all iron oxides (iron-stained alteration) associated with the above-mentioned lithological units ( Figure 10B). This similarity of iron oxide mapping results might be related to extensive iron-stained occurrences in the lithological units. The spatial distribution of dolomite/Fe 2+ (bright pixels) is typically mapped in the eastern part of the selected subset ( Figure 10C). The PE r unit, PE Ed unit and the C 1d unit are highlighted in the image-map of the 7/4 band ratio ( Figure 10C). Figure 10D shows RGB false color composite of 7/4, 5/3 and 5/2 band ratio images for the selected subset. The discrimination of lithological units is characteristically discernable based on different composition of iron oxide/hydroxide minerals.

Detailed Detection of Iron Oxide/Hydroxides and Dolomit Using WV-3
The lithological units with a high abundance of iron oxide/hydroxide minerals such as the DCg unit appear in cyan. The exposures of the PE r , PE Ed and C 1d units contain the high abundance of dolomite/Fe 2+ represent as whitish-yellow, golden yellow and light brown. On the other hand, the Kuhbanan Formation (E 1 ), Q t deposits and the C ph.s unit having a low abundance of iron oxides show recognizable colors (shades of gray) and boundaries with other lithologies (see Figure 10D).
Analyzing the PCA statistical results for the WV-3 selected subset shows the PC3, PC4 and PC6 contain essential information for mapping iron-stained alteration, dolomite/Fe 2+ and Fe 3+ oxides, respectively. The PC3 shows high negative loading in band 2 (−0.598506) and strong positive loading in band 5 (0.661001) (see Table 5). Thus, iron-stained alteration can be mapped as bright pixels because of band 5 that is assumed as a reflection band. Figure 11A shows the DCg unit and some small parts of the PE r unit as bright pixels. The PC4 has moderate positive loading in band 2 (0.345377), strong positive loading in band 4 (0.631659) and high negative loading in band 7 (−0.509311) (see Table 5). Dolomite/Fe 2+ can be mapped as dark pixels due to the negative sign in the reflection band (band 7). Figure 11B shows the negated image-map of the PC4. The PE r unit, the PE Ed unit and the C 1d unit are represented as bright pixels. The image-map of the PC4 is identical to the 7/4 band ratio resultant (see Figure 10C). The PC6 contains a strong contribution in band 2 (−0.427799) and band 3 (−0.515988) with a negative sign and strong loading in band 5 (0.646312) with a positive sign (see Table 5). Accordingly, Fe 3+ oxides manifest in bright pixels because of positive sign in band 5, which is assumed as a reflection band ( Figure 11C). The DCg unit, the PE r unit and the PE Ed unit are mapped in the PC6 image-map. Figure 11D shows RGB false color composite image-map of the PC4, PC6 and PC3, respectively. Lithological units are differentiated stronger than RGB false color composite of band ratios (see Figure 10D). The lithological boundaries of the PE r unit, the PE Ed unit and the C 1d unit (represented as golden to orange-yellow) with other lithological units such as the DCg, the C ph.s and the E 1 units (depicted in shades of gray and purple) are distinguishable ( Figure 11D).

Generating Mineral Prospectivity Maps for the Study Area
The fuzzy-logic model was utilized to produce mineral prospectivity maps of the favorable areas for ore mineralizations in the study region using most rational alteration thematic layers derived from image processing techniques (see Table 6). In this analysis, the PCA output was considered more informative compared to band ratios output. It is due to the fact that the PCA is statistically based algorithm and uses uncorrelated linear combinations (eigenvector loadings) to map spectral characteristics of alteration minerals. Therefore, the PC4 and PC5 thematic layers of Landsat-8 and Sentinel-2 and PC4, PC5 and PC8 thematic layers of ASTER were selected to be integrated by application of the fuzzy gamma operator (γ = 0.7) for generating mineral prospectivity map at a regional scale. Figure 12 shows the regional mineral prospectivity map of the study area. Evaluating the fuzzy membership indicates that the high favorability index is associated with some of the lithological units in the study area. The PE1s unit, the PE r unit, the Kuhbanan Formation (E 1 ), the TR unit, the C ph.s unit, the Ngc unit, the Zarigan granite (Zrg) and the DCg unit show high value (0.7 to 1.0) of the favorability index. Most of the mines and ore indications are located in the high value zones of the favorability index and associated with fault systems in the study area ( Figure 12). Results demonstrate that the PE 1s unit, the PE r unit and the Zarigan granite (Zrg) show the highest value (0.9 to 1.0) of the favorability index. The alteration zones associated with these lithological units are the most favorable/prospective areas for ore mineralizations at the regional scale. Some of the high prospective zone demarcated using dashed black ellipsoids and circles in Figure 12, which can be considered for future mineral exploration in the study region.
Remote Sens. 2020, 12, x FOR PEER REVIEW 21 of 35 composite of 7/4, 5/3 and 5/2 band ratio images for the selected subset. The discrimination of lithological units is characteristically discernable based on different composition of iron oxide/hydroxide minerals. The lithological units with a high abundance of iron oxide/hydroxide minerals such as the DCg unit appear in cyan. The exposures of the PE r , PE Ed and C 1d units contain the high abundance of dolomite/Fe 2+ represent as whitish-yellow, golden yellow and light brown. On the other hand, the Kuhbanan Formation (E 1 ), Q t deposits and the C ph.s unit having a low abundance of iron oxides show recognizable colors (shades of gray) and boundaries with other lithologies (see Figure 10D).
Analyzing the PCA statistical results for the WV-3 selected subset shows the PC3, PC4 and PC6 contain essential information for mapping iron-stained alteration, dolomite/Fe 2+ and Fe 3+ oxides, respectively. The PC3 shows high negative loading in band 2 (−0.598506) and strong positive loading in band 5 (0.661001) (see Table 5). Thus, iron-stained alteration can be mapped as bright pixels  Figure 13 shows the local mineral prospectivity map of the study area derived from the PC3, PC4 and PC6 thematic layers (most rational alteration thematic layers) for WV-3 data. The fuzzy fuzzy gamma operator (γ = 0.7) was used to fuse the selected alteration thematic layers (see Table 3). The highest value of (0.8 to 1.0) the favorability index is obtained for the PE r unit, the PE ED unit and the C 1d unit. In addition, the DCg unit shows a high value (0.6 to 0.9) of the favorability index in some parts of the study area. The Pb-Zn mineralization zones have moderate to high favorability index value (0.6 to 1.0) and are also adjoining to fault systems ( Figure 13). Accordingly, the most favorable/prospective zones for Pb-Zn mineralization in the study area are alterations (especially dolomitic zone) associated with the PE r unit, the PE ED unit and the C 1d unit, mainly in fault contact zones with impermeable lithological units. Black polygons, ellipsoids and circles show some of the high prospective zones for future mineral prospecting in the study region ( Figure 13). These high prospective zones were selected to check during filed reconnaissance in this study. Locations of in situ observation are shown in Figure 13.
Remote Sens. 2020, 12, x FOR PEER REVIEW 23 of 35 Figure 12. Mineral prospectivity map of the study area at regional scale derived from Landsat-8, Sentinal-2 and ASTER selected alteration thematic layers. Magenta polygon delineates WV-3 imagery. Figure 13 shows the local mineral prospectivity map of the study area derived from the PC3, PC4 and PC6 thematic layers (most rational alteration thematic layers) for WV-3 data. The fuzzy fuzzy gamma operator (γ=0.7) was used to fuse the selected alteration thematic layers (see Table 3). The highest value of (0.8 to 1.0) the favorability index is obtained for the PE r unit, the PE ED unit and the C 1d unit. In addition, the DCg unit shows a high value (0.6 to 0.9) of the favorability index in some parts of the study area. The Pb-Zn mineralization zones have moderate to high favorability index value (0.6 to 1.0) and are also adjoining to fault systems ( Figure 13). Accordingly, the most favorable/prospective zones for Pb-Zn mineralization in the study area are alterations (especially dolomitic zone) associated with the PE r unit, the PE ED unit and the C 1d unit, mainly in fault contact zones with impermeable lithological units. Black polygons, ellipsoids and circles show some of the high prospective zones for future mineral prospecting in the study region ( Figure 13). These high prospective zones were selected to check during filed reconnaissance in this study. Locations of in situ observation are shown in Figure 13.

Verifying the Spatial Distribution of Alteration Zones and High Prospective Areas
Field reconnaissance was conducted to verify the occurrence of alteration zones and high prospective locations of Pb-Zn mineralization in the study area. The global positioning system (GPS) survey was undertaken for estimating the overall accuracy of the image processing techniques. Several field photos were taken to record the location, contact and characteristics of alteration zones and lithological units. Rock sample collection was carried out in several parts of the study regions and high prospective zones for laboratory analysis such as polish and thin sections and XRD analysis and the ASD spectroscopy. Surface manifestations of mineral assemblages such as iron oxide/hydroxides, dolomite, shale, calcite and gypsum are widespread in the study area. Some of the highly prospective zones, mineralogically interesting alteration zones and lithological units in the

Verifying the Spatial Distribution of Alteration Zones and High Prospective Areas
Field reconnaissance was conducted to verify the occurrence of alteration zones and high prospective locations of Pb-Zn mineralization in the study area. The global positioning system (GPS) survey was undertaken for estimating the overall accuracy of the image processing techniques. Several field photos were taken to record the location, contact and characteristics of alteration zones and lithological units. Rock sample collection was carried out in several parts of the study regions and high prospective zones for laboratory analysis such as polish and thin sections and XRD analysis and the ASD spectroscopy. Surface manifestations of mineral assemblages such as iron oxide/hydroxides, dolomite, shale, calcite and gypsum are widespread in the study area. Some of the highly prospective zones, mineralogically interesting alteration zones and lithological units in the study area were particularly visited. The argillic alteration, sericitic zones, iron oxides and dolomitization were found in the PE r unit, the PE ED unit, the E 1 unit and the C 1d unit, which showed some surface expression of hematite, malachite, pyrite, galena and sphalerite. Some old open-pit quarries were also found in the alteration zones of the lithological units. In many parts of the study area, surface expression of Pb-Zn mineralization was typically observed in the fault contact of dolomite with other lithological units ( Figure 14A-F).
Polish section study shows the presence of chalcopyrite, pyrite, malachite, smithsonite, sphalerite, galena, hematite and limonite. Thin section study typically displays the association of iron oxide/hydroxides with dolomite that mostly concentrated in the fractures. The XRD analysis of the samples collected from alteration zones inside old open-pit quarries and surrounding areas shows the presence of quartz, dolomite, calcite, muscovite, chlorite, gypsum, albite, illite, jarosite and malachite. The ASD spectroscopy analysis for main lithological units such as shale, gypsum, dolomite and calcite were measured in this study ( Figure 15). The laboratory reflectance spectra from the shale sample display three distinguishable absorption features about 1.40 µm attributed to OH/H 2 O stretches, 1.90 µm related to H 2 O stretches and 2.20 µm due to combination of the OH-stretching fundamental with Al-OH bending mode, respectively ( Figure 15). These absorption features occur in shale because of the high content of clay minerals (Al-rich phyllosilicates) [103]. The reflectance spectra derived from gypsum exhibits three absorption features, which are identical to the shale sample ( Figure 15). But, the absorption feature at 1.90 µm (due to H 2 O stretches) is stronger than the shale sample and the absorption feature near 2.20 µm is related to S-O bending mode [104]. The reflectance spectra of dolomite contain two main absorption features related to Fe 2+ at 0.9 to 1.2 µm and CO 3 in 2.35 µm ( Figure 15). Calcite absorption properties typically concentrated about 2.35 µm (Figure 15), which is attributed to the vibrational processes of the CO 3 bending mode [82]. Furthermore, the confusion matrix and Kappa Coefficient were calculated for alteration mapping results versus field GPS surveys ( Table 7). The overall accuracy for Landsat-8 and Sentinel-2 datasets is 86.66% and 83.33%, respectively. The Kappa Coefficient of 0.83 for Landsat-8 and 0.81 for Sentinel-2 is also assessed. For the ASTER dataset, the overall accuracy and Kappa Coefficient are 70% and 0.68, respectively. The calculation of the confusion matrix for WV-3 shows the overall accuracy of 83.33% and Kappa Coefficient of 0.82 (Table 7). zones were identified (see Figure 12), which are mostly associated with alteration zones in the PE 1s unit, the PE r unit, the Zarigan granite (Zrg), the Ngc, the C ph.s and PE Ed lithological units. Most of the prospective zones are located adjacent to the NW-SE and N-S fault systems, which likely acted as fluid pathways for hydrothermal ore mineralizations. Documented mineral occurrences show also close spatial relationship to the fault systems in the study area. The ASTER dataset was used for detailed mapping iron oxide/hydroxides, Al/Fe-OH minerals, Mg-Fe-OH/CO3 minerals and dolomite associated with carbonate-hosted Pb-Zn deposits. The spatial distribution of these alteration minerals was comprehensively mapped. Results derived from the PCA for the study area show that calcite, gypsum, hematite and goethite are the main spectrally-spatially minerals. However, dolomite, jarosite, kaolinite and muscovite are minor spectrally-spatially minerals. The analysis of confusion matrix shows that calcite has the highest user's accuracy (83%) and producer's accuracy (83%). However, the lowest user's accuracy (57%) was obtained for argillic alteration + Shale (kaolinite and muscovite) (see Table 7). It means that the spectral mixing and confusion of calcite with other alteration minerals is low, while for argillic alteration + Shale is high. Iron oxide/hydroxides shows the user's accuracy of 80% and producer's accuracy of 67%. Gypsum and dolomitization have similar user's accuracy of 67% and producer's accuracy of 67%. Thus, spectral mixing with other alteration minerals especially dolomite could be Moreover, dolomite has some spectral mixing with calcite. ASTER dataset shows the overall accuracy of 70% and Kappa Coefficient of 0.68. Accordingly, ASTER can map and discriminate different alteration minerals appropriately, but some spectral mixing and confusion are also associated with alteration mapping. ASTER PC image-maps show iron oxide/hydroxides, OH/S-O/CO3-bearing minerals and dolomite associated with the C 1d unit, the PE r unit, the PE ED unit, the Kuhbanan Formation (E 1 ) and the C ph.s are high favorable/prospective zones for carbonate-hosted Pb-Zn mineralization (see Figure 12). The N-S faults are dominant structural features associated with the high potential zones; however, the NW-SE and NE-SW faults are also associated with few of the prospective zones.
More detailed surface distribution of Fe 3+ oxides, iron-stained alteration and dolomite in the highly prospective zones were identified using WV-3 dataset. Dolomite is mostly identified associated with the PE r unit, PE Ed unit and the C 1d unit. Calculation of user's accuracy (83%) and producer's accuracy (83%) for iron oxide/hydroxides and dolomite shows a low rate of spectral mixing and confusion. The overall accuracy of 83.33% and Kappa Coefficient of 0.82 show a very good rate of accuracy and agreement for WV-3 dataset (see Table 7). WV-3 mineral prospectivity map (see Figure 13) shows the zones of most favorable/prospective zones for the Pb-Zn mineralization in the study area. The dolomitic zones in the PE r unit, the PE ED unit and the C 1d unit are typically indicated as the highest potential zone. In the vicinity of N-S, NW-SE and NE-SW fault systems, several of the most prospective zones are identified in the study area (see Figure 13).
Fieldwork was conducted comprehensively in the highly prospective zones for observing surface expression of Pb-Zn mineralization and related alteration zones and lithological units. Iron oxide/hydroxides and dolomite as well as weak argillic/sericitic zones was found with surface expression of hematite, malachite, pyrite, galena and sphalerite in the PE r unit, the PE ED unit, the E 1 unit and the C 1d unit. Surface manifestation of Pb-Zn mineralization was typically recorded in the fault contact of dolomite with impermeable lithological units. The presence of quartz, dolomite, calcite, muscovite, chlorite, gypsum, albite, illite, jarosite and malachite is verified by XRD analysis. These evidences emphasized that the fault systems provided fluid conduits for Pb-Zn mineralization that hydrothermally altered the host lithologies and afterwards oxidized during supergene

Discussion
Remote sensing satellite imagery has been successfully utilized to detect the major hydrothermal alteration minerals associated with the carbonate-hosted Pb-Zn deposits in metallogenic provinces around the world [5,[16][17][18][19]. In this investigation, multiple sources of spectral data derived from Landsat-8, Sentinel-2, ASTER and WV-3 sensors were utilized for detailed mapping a variety of hydrothermal alteration minerals in the central part of the KKTZ region, the CIT, Iran. Band ratios and PCA image processing techniques were used to produce thematic maps of hydrothermal alteration minerals for indicating the high prospective zones. Lithological units and alteration mineral zones were mapped based on spectral absorption characteristics of Fe 3+ /Fe 2+ and Al-OH/CO3 minerals by implementing band ratios and PCA techniques to spectral bands of Landsat-8 and Sentinel-2. The spatial distribution of iron oxide/hydroxide minerals was mapped in the lithological units of the PE 1s , the Ngc, the K 1 , the PE r , the DCg and Q t deposits using the 4/2 band ratio and PC5 of Landsat-8 and Sentinel-2. Analyzing confusion matrix for mapping iron oxide/hydroxide minerals using these sensors (see Table 7) indicates that the highest the user's accuracy (71%) and producer's accuracy (83%) were obtained for Sentinel-2 dataset. It shows that the higher spatial resolution of Sentinel-2 (10 m) clearly impacts on mapping minerals. The 6/7 band ratio of Landsat-8 and 11/12 band ratio of Sentinel-2 and the PC4 images were used to identify the surface distribution of hydroxyl-bearing alteration minerals and carbonates in the study region. The alteration zones associated with lithological units of the PE 1s , the E, the K 1 , the PE r , the Zrg, the C 1d and the Q t were highlighted. The Landsat-8 show user's accuracy (91%) and producer's accuracy (91%) for hydroxyl-bearing alteration minerals and carbonates, which are higher than user's accuracy (90%) and producer's accuracy (83%) of Sentinel-2 (see Table 7). It probably can be attributed to spectral coverage of band 12 (2.100-2.280 µm) of Sentinel-2, which is not able to map carbonates, properly. The overall accuracy of Landsat-8 for mapping target alteration minerals is 86.66% and higher than Sentinel-2 (83.33%) because of better spectral band placement in Landsat-8 for detection of hydroxyl-bearing alteration minerals and carbonates. Fuzzy logic modeling was used to fuse the most informative thematic alteration layers (the PC4 and PC5 images). Subsequently, mineral prospectivity map for the study area was generated. Several prospective zones were identified (see Figure 12), which are mostly associated with alteration zones in the PE 1s unit, the PE r unit, the Zarigan granite (Zrg), the Ngc, the C ph.s and PE Ed lithological units. Most of the prospective zones are located adjacent to the NW-SE and N-S fault systems, which likely acted as fluid pathways for hydrothermal ore mineralizations. Documented mineral occurrences show also close spatial relationship to the fault systems in the study area.
The ASTER dataset was used for detailed mapping iron oxide/hydroxides, Al/Fe-OH minerals, Mg-Fe-OH/CO3 minerals and dolomite associated with carbonate-hosted Pb-Zn deposits. The spatial distribution of these alteration minerals was comprehensively mapped. Results derived from the PCA for the study area show that calcite, gypsum, hematite and goethite are the main spectrally-spatially minerals. However, dolomite, jarosite, kaolinite and muscovite are minor spectrally-spatially minerals.
The analysis of confusion matrix shows that calcite has the highest user's accuracy (83%) and producer's accuracy (83%). However, the lowest user's accuracy (57%) was obtained for argillic alteration + Shale (kaolinite and muscovite) (see Table 7). It means that the spectral mixing and confusion of calcite with other alteration minerals is low, while for argillic alteration + Shale is high. Iron oxide/hydroxides shows the user's accuracy of 80% and producer's accuracy of 67%. Gypsum and dolomitization have similar user's accuracy of 67% and producer's accuracy of 67%. Thus, spectral mixing with other alteration minerals especially dolomite could be feasible for iron oxide/hydroxides. The spectral mixing and confusion between gypsum and argillic alteration + shale can be expected.
Moreover, dolomite has some spectral mixing with calcite. ASTER dataset shows the overall accuracy of 70% and Kappa Coefficient of 0.68. Accordingly, ASTER can map and discriminate different alteration minerals appropriately, but some spectral mixing and confusion are also associated with alteration mapping. ASTER PC image-maps show iron oxide/hydroxides, OH/S-O/CO3-bearing minerals and dolomite associated with the C 1d unit, the PE r unit, the PE ED unit, the Kuhbanan Formation (E 1 ) and the C ph.s are high favorable/prospective zones for carbonate-hosted Pb-Zn mineralization (see Figure 12). The N-S faults are dominant structural features associated with the high potential zones; however, the NW-SE and NE-SW faults are also associated with few of the prospective zones.
More detailed surface distribution of Fe 3+ oxides, iron-stained alteration and dolomite in the highly prospective zones were identified using WV-3 dataset. Dolomite is mostly identified associated with the PE r unit, PE Ed unit and the C 1d unit. Calculation of user's accuracy (83%) and producer's accuracy (83%) for iron oxide/hydroxides and dolomite shows a low rate of spectral mixing and confusion. The overall accuracy of 83.33% and Kappa Coefficient of 0.82 show a very good rate of accuracy and agreement for WV-3 dataset (see Table 7). WV-3 mineral prospectivity map (see Figure 13) shows the zones of most favorable/prospective zones for the Pb-Zn mineralization in the study area. The dolomitic zones in the PE r unit, the PE ED unit and the C 1d unit are typically indicated as the highest potential zone. In the vicinity of N-S, NW-SE and NE-SW fault systems, several of the most prospective zones are identified in the study area (see Figure 13).
Fieldwork was conducted comprehensively in the highly prospective zones for observing surface expression of Pb-Zn mineralization and related alteration zones and lithological units. Iron oxide/hydroxides and dolomite as well as weak argillic/sericitic zones was found with surface expression of hematite, malachite, pyrite, galena and sphalerite in the PE r unit, the PE ED unit, the E 1 unit and the C 1d unit. Surface manifestation of Pb-Zn mineralization was typically recorded in the fault contact of dolomite with impermeable lithological units. The presence of quartz, dolomite, calcite, muscovite, chlorite, gypsum, albite, illite, jarosite and malachite is verified by XRD analysis. These evidences emphasized that the fault systems provided fluid conduits for Pb-Zn mineralization that hydrothermally altered the host lithologies and afterwards oxidized during supergene processes. The ASD spectroscopy analysis of main lithological units showed distinct absorption features related to Fe 2+ , OH/H 2 O, H 2 O, Al-OH, S-O and CO 3 for collected rock samples such as dolomite, shale, gypsum and calcite. Generally, confusion matrix and Kappa Coefficient calculated for alteration mapping results versus field GPS survey show reasonable overall accuracy (70% to 86.66%) and good rate of agreement (0.68 to 0.83). Landsat-8 and Sentinel-2 generally mapped iron oxide/hydroxides and hydroxyl-bearing and carbonate minerals and indicated favorable/prospective zones, while ASTER and WV-3 comprehensively detected and discriminated hematite, goethite, jarosite, gypsum, calcite, dolomite, kaolinite and muscovite and the highest prospective zones.

Conclusions
This investigation demonstrates the application of multi-source spectral data in the range of VNIR and SWIR wavelengths provided by Landsat-8, Sentinel-2, ASTER and WV-3 for indicating the highly prospective zones of carbonate-hosted Pb-Zn deposits in the KKTZ of the CIT region of Iran. Results show that significant information related to the iron oxide/hydroxide and hydroxyl-bearing and carbonate minerals can be easily obtained by implementing some specialized band ratio (e.g., 4/2, 6/7 and 11/12) and PCA technique to Landsat-8 and Sentinel-2 datasets. Fusing of the most informative alteration thematic layers by the fuzzy-logic model is a reliable approach for generating remote sensing-based mineral prospectivity map. Landsat-8/Sentinel-2/ASTER mineral prospectivity map for the KKTZ indicated the potential zones are mostly located in the altered zones of the PE 1s unit, the PE r unit and the Zarigan granite (Zrg), which normally are near to the NW-SE and N-S fault systems. Detailed surface distribution of Al/Fe-OH minerals, Mg-Fe-OH/CO3 minerals and dolomite was detected using RBD and PCA methods to ASTER dataset. The alteration zones (especially gossan and dolomite) associated with C 1d , the PE r , the PE ED , the Kuhbanan Formation (E 1 ) and the C ph.s lithological units and close to N-S, NW-SE and NE-SW faults are considered highly favorable/prospective zones for carbonate-hosted Pb-Zn mineralization in the study area. Detailed identification of Fe 3+ oxides, iron-stained alteration and dolomite in the highly prospective zones was obtained using WV-3 VNIR spectral bands processing (e.g., band ratios of 5/3, 5/2 and 7/4; PCA techniques). The most favorable/prospective zones for Pb-Zn mineralization in the study area are dolomite and gossan alteration zones located in the PE r , the PE ED and the C 1d lithological unit, which are exclusively placed in the fault contact zones of dolomitic occurrences with impermeable lithological units. The N-S, NW-SE and NE-SW trends fault systems provided fluid conduits for Pb-Zn mineralization and deposition and subsequent alteration zones in the study area. Therefore, it is recommended that detailed structural analysis of lineaments using Synthetic Aperture Radar (SAR) remote sensing data such as the Phased Array type L-band Synthetic Aperture Radar (PALSAR) is also required for future mineral prospection in the CIT. Additionally, SWIR bands of WV-3 can provide high spectral and spatial data for detailed alteration mapping. On the other hand, thermal infrared (TIR) data from Landsat-8 and ASTER are also valuable to map minerals like quartz, andradite, gypsum, calcite, dolomite, diopside and albite, comprehensively. In conclusion, the application of multi-sensor remote sensing satellite imagery and fusing the most informative alteration thematic layers using the fuzzy-logic model can provide a low-cost exploration approach for prospecting carbonate-hosted Pb-Zn mineralization in the KKTZ of CIT and other metllogenic provinces around the world. The results demonstrated in this investigation represent a significant contribution of space-borne multispectral systems to generate mineral prospectivity maps at various scales (regional, district and local scale). This approach could be very interesting for stakeholders and mining/exploration companies to use different types of space-borne multispectral data for distinct phases of mineral exploration.