Geochemical and Hydrothermal Alteration Patterns of the Abrisham-Rud Porphyry Copper District, Semnan Province, Iran

: In this study, the zonality method has been used to separate geochemical anomalies and to calculate erosional levels in the regional scale for porphyry-Cu deposit, Abrisham-Rud (Semnan province, East of Iran). In geochemical maps of multiplicative haloes, the co-existence of both the supra-ore elements and sub-ore elements local maxima implied blind mineralization in the northwest of the study area. Moreover, considering the calculated zonality indices and two previously presented geochemical models, E and NW of the study have been introduced as ZDM and BM, respectively. For comparison, the geological layer has been created by combining rock units, faults, and alterations utilizing the K-nearest neighbor (KNN) algorithm. The rock units and faults have been identiﬁed from the geological map; moreover, alterations have been detected by using remote sensing and ASTER images. In the geological layer map related to E of the study area, many parts have been detected as high potential areas; in addition, both geochemical and geological layer maps only conﬁrmed each other at the south of this area and suggested this part as high potential mineralization. Therefore, high potential areas in the geological layer map could be related to the mineralization or not. Due to the incapability of the geological layer in identifying erosional levels, mineralogy investigation could be used to recognize this level; however, because of the high cost, mineralogy is not recommended for application on a regional scale. The ﬁndings demonstrated that the zonality method has successfully distinguished geochemical anomalies including BM and ZDM without dependent on alteration and was able to predict erosional levels. Therefore, this method is more powerful than the geological layer.


Introduction
The utilization of geochemical methods for ore deposit exploration dates back to 1930. Fersman (1939) carried out the first survey of such an exploration [1]. Since then, further studies on the theory and application of geochemical exploration methods have been carried out, and these techniques have been increasingly modified and improved. Mining geochemistry is a branch of applied geochemistry, which is based on the utilization of geochemical methods that helps increase the ore reserves of known mines by assessing the ore potential of deep horizons. In other words, local and mine scale exploration models for anomaly recognition are created and developed by using mining geochemistry. Recent experiences in the application of mining geochemistry illustrate its efficiency in the discovery of blind and zone-dispersed mineralization (BM and ZDM) within areas of active and abandoned mines. Due to increasing ore reserves and mining income, this trend in geochemical exploration is very important [2]. The recognition of various alteration zones is a qualitative method that cannot help geochemists in separating BM from ZDM at a local scale. In mining geochemistry, the alteration has no basic role in separating anomalies at a local scale [3,4]. Optimal drilling points were determined by using mining geochemistry as a quantitative method without being time consuming and inducing high costs. In the past decades, several models and methods based on geochemistry have been developed for predicting geochemical anomalies as well as the locations of hidden orebodies [2,[5][6][7]. Most of these models and methods are focused on the identification of geochemical anomalies reflecting the presence of hidden orebodies and the prediction of horizons of erosional surfaces [8]. In active mines, vertical geochemical zonality is the most important feature of primary halos because of the relation to the direction of the ore-bearing fluid [9][10][11]. Beus and Grigorian (1977) used vertical geochemical zonality to predict hidden mineralization at the mine scale [12]. Grigorian (1985Grigorian ( , 1992) presented a zonality model to identify BM from ZDM [5,6]. Since then, the zonality method has been used in many studies [2][3][4][13][14][15][16][17][18][19][20][21]. Solovov (1987) used metallometric methods for the identification of geochemical anomaly (IGA) and the quantitative evaluation of ore reserves [22]. Baranov (1987) introduced a model in which horizons of the erosional surface were computed for geochemical associations [23]. Solovov (1990) introduced different relations to predict hidden orebodies by using metallometric exploration [24]. Liu and Peng (2004) presented a model to predict hidden orebodies by the synthesis of geological, geophysical, and geochemical information based on a dynamic approach [25].
In most mineral exploration methods (e.g., porphyry-Cu), a mineral potential map is obtained by using one layer or a combination of different layers [26][27][28][29][30][31], which includes field geological surveys, geochemical surveys, field geophysical surveys, and remote sensing [32][33][34]. This map consists of shallow to the deep layers, which poses a problem when these layers are not associated with mineralization. Each of these layers has a value for mapping the areas with the potential of mineralization, and a number of these layers are generally surficial and cannot be useful for identifying BM. Hydrothermally altered rocks result from chemical attacks of pre-existing rocks by hydrothermal fluids. The spatial distribution of hydrothermally altered rocks is a key to locating the main outflow zones of hydrothermal systems, which may result in the recognition of mineral deposits [35]. Minerals associated with alteration can be detected by remote sensing. These tasks are achieved by using the analysis of the spectral signatures recorded in the visible-near infrared (VNIR), short wave infrared (SWIR), and thermal infrared (TIR) regions of the electromagnetic spectrum with this spectral signature constituting the key mineral identification criterion [36]. Furthermore, the mineral deposits are spatially and genetically associated with the various types of geological structures including faults or fractures [37]. Faults and fractures, which transport magmatic, meteoric, and metal carrying hydrothermal fluids, subsequently deposit metals [38]. Zarasvandi [38][39][40][41][42]. The strong advance in remote sensing allows exploiting a variety of sources and methods in the characterization of lineaments [43]. Remote sensing is a valuable technical resource for mineral exploration when it is properly employed [44][45][46][47][48][49][50][51][52][53].
In this study, the results of the zonality method were compared with the geological layer including rock units, faults, and alterations. For this goal, a part of the 1:100,000 scale map of Abrisham-Rud (Semnan province, East of Iran) was examined. This area is a part of the Troud Range in the Khorasan porphyry tract, which few studies on porphyry mineralization have been conducted [54]. Orojnia (2003) studied the lithology and provenance of Eocene volcanic rocks in this area, which suggested identifying the economic potential of Cu [55]. Mars (2014) suggested the potential of hydrothermal alteration in the Torud Range that could be associated with an unidentified porphyry system [56]. Thus, study data are concentrated in the province of Abrisham-Rud, where it has high propects for porphyry copper mineralization.

Geological Setting of the Study Area
The study area is a part of the 1:100,000 scale map of Abrisham-Rud (Semnan province, East of Iran), which is located in the north of the Central Iran zone [57], and it is a part of Troud Range in Khorasan porphyry tract (Figure 1a). Khorasan porphyry tracts are delimited by permissive units of island arc setting of Late Cretaceous to Early Miocene age [58]. This tract includes four main ranges: Taknar-Kashmar, Kuh Mish, Sabzevar, and Torud [54]. Igneous units of this tract are shown in Figure 1a, along with the location of known porphyry-related mineral occurrences and other geologic features mentioned in this section. In the Torud Range (west of Khorasan porphyry tract), middle Eocene volcanosedimentary rocks are overlain by Eocene-Oligocene calc-alkaline and alkaline volcanic rocks, which are interlayered with shallow marine, lacustrine, and subaerial sedimentary successions. These successions are intruded by basic tholeiitic dikes and calc-alkaline quartz monzodioritic to granodioritic stocks [59,60]. The aluminum content of hornblendes indicates shallow emplacement depths [61]. In this range, Chah Shirin and Chah Mussa have been introduced as porphyry/porphyry-related deposits [58]. The study area has a dry climate, mountainous topography, and poor vegetation cover. The oldest geological units are metamorphic rocks for which its bedrock was metamorphosed in the Late Triassic. The latest units are associated with a quaternary that consists of gravel fan, terrace, clay and salt deposits, and channel deposits (Figure 1b). Most volcanic activity occurred in semi-arid settings in the form of lava. The middleupper Eocene rock units are the most extended in the study area and mostly include volcanic breccia-agglomerates and tuffs, intermediate lavas, basic and acidic rocks, and pyroclastic-sedimentary rocks. These units have outcrops and are sometimes dispersed or indistinguishable, and there are many faults in them. The volcanic breccia-agglomerate units that sometimes had intermediate lavas observed in red-brown and sometimes had dusty colors (argillic alteration), as well as the presence of green minerals (chlorite, epidote) in the volcanic fragments of this breccia, indicate lava explosive eruption in a semi-arid and shallow setting with the effect of water on alterations [55].

Geological Data
The 1:100,000 scale map of Abrisham-Rud was purchased from the Geological Survey of Iran (GSI). As shown in Figure 1b, a part of this geological map was selected.

Geochemical Data
The geochemical database of the 1:100,000 scale map of Abrisham-Rud was also purchased from GSI. The 364 samples have been surveyed in a systematic network with intervals of 1400 × 1400 m 2 (Figure 2), and these soil samples have been analyzed using the ICP-MS method. The concentrations of Cu, Mo, Pb, and Zn were considered for this study.

ASTER Data
Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) sensors are one of the multi-spectral sensors that have been installed on the Terra satellite. ASTER measures reflected radiation in three bands between 0.52 and 0.86 µm (VNIR) and in six bands between 1.60 and 2.43 µm (SWIR) and five bands of emitted radiation in the 8.125 µm to 11.65 µm (TIR) region with 15 m, 30 m, and 90 m resolution, respectively [66][67][68]. In this study, an ASTER-Level 1T image was utilized. ASTER-Level 1T data contain calibrated at-sensor radiance, which corresponds with ASTER-Level 1B that has been geometrically corrected and rotated to the north up UTM projection.

Data Preparation
Atmospheric correction was used to minimize influences of atmospheric factors in multispectral data. The Internal Average Relative Reflection (IARR) method was applied to Sentinel-2 and ASTER data. The IARR technique for mineral mapping requires no prior knowledge of geological features [69]. The average radiance for each band of the image was calculated; therefore, an average spectrum was created, and this average spectrum was divided into actual radiance for each band of each pixel to create an image of apparent reflectance. It has been suggested as the best method for arid areas with no vegetation cover [70].

Zonality Method
A zonation of a geochemical halo has a spatial nature and vectorial context that can be defined by the three parameters of dimension (space), direction, and element concentration [4]. Recognition of zonality of geochemical halos associated with BM can be achieved using four cases of complementary analyses [2]: (1) analysis of element associations representing supra-ore and sub-ore halos of mineral deposits; (2) analysis of a single component, implying false anomaly; (3) analysis of mean values of indicator elements outside significant geochemical anomalies to eliminate background noise in data analysis; and (4) mapping of multiplicative geochemical anomalies.
One of the most important indices in porphyry-Cu deposits is the ratio of Pb and Zn to Cu and Mo, which is often defined as a zonality index [2,12,71]. This index represents different exhumation levels of mineral deposits [6]. Input variables can be subdivided into the following: supra-ore, upper-ore, ore, lower-ore, and sub-ore [24]. Ziaii et al. (2012) and Ziaii et al. (2009) showed that these groups provide the necessary information to separate BM from ZDM in porphyry-Cu mineralization [2,3].

Anomaly Separation
For mapping the multiplicative index of supra-ore (Pb × Zn) and sub-ore (Cu × Mo) elements and calculating the zonality index, the threshold value is calculated by using Equation (1) for each element: where C A is the anomaly value, and C x = C o is the geometric mean of the elements contained within the background area, which is calculated by using Equation (2): where C i is the element concentration of samples, N is the number of samples, and ε is generally called a standard factor, which is calculated from Equation (3). Due to the above consideration as to the selection of the value of t, the lowest anomalous content, when trying to detect weak anomalies determined from one sampling point, in geochemical prospecting it is taken equal to C Al ≥ C o ε 3 . This relationship corresponds to a "three standard deviations" criterion extensively used in many engineering disciplines to determine quantities falling outside the probable values of a random anomaly distributed quantity. Thus, weak anomalies formed by a sequence of adjacent sampling points with increased pathfinder elements below C Al can be detected, and it is conventional to lower the threshold value according to the criterion C Am ≥ C o ε 3/ √ m , where m = 2, 3, 4 . . . 9 is the number points that may be joined, which can show a common anomaly in the geochemical map [22].

Erosional Surface
In order to predict the erosional level, Solovov (1987) suggested using areal productivity, and Beus and Grigorian (1977) suggested using the coefficient of mineralization to eliminate the syngenetic parameters of the halos, which increase anomaly detection [12,22]. It should be noted that both the areal productivity and mineralization coefficients are used to calculate the vertical geochemical zonality index.
In the systematic sampling network, linear productivity is calculated according to Equation (4): where M is the linear productivity, ∆x is the distance between the samples in each profile, C x is the values greater than the anomaly concentration, and n is the number of anomal samples. If the values of M i are preliminarily estimated in each of m profiles across the anomaly, P is determined according to Equation (5) [22]: where 2l is the distance between profiles. Therefore, the zonality index introduced by Solovov (1987) is calculated for zones I and II by using Equation (6) [22].
The zonality index introduced by Beus and Grigorian (1977) is calculated for each zone by using Equation (7) [12]: where CA is the arithmetic mean of element contents, and η(α) is the coefficient of mineralization calculated for each element by using Equation (8): where η(α) ore is the number of anomalous samples, and η A (α) is the total number of samples in each zone [12].

Lineaments Extraction
O'Leary et al. (1976) defined the term "lineaments" as a simple or composite linear feature of a surface for which its parts are aligned in a rectilinear or slightly curvilinear relationship and differs from the pattern of adjacent features and probably reflects some subsurface phenomena [72]. Faults, fractures, and large crush zones are formed by extension or compression processes and are considerable and fundamental factors on ore mineral deposition. Areas with concentrations or intersections of these structures could be suitable for the penetration of magma, ore-forming solutions, and, afterward, mineralization [73]. In other words, a detailed geological study imperatively means acquiring knowledge of present structural information, principally the lineaments [74]. The significance of lineaments is also manifested by their localization often close to several mineralogical deposits [75].
Lineament extraction methods can be conducted by using manual photointerpretation by an expert, semi-automatic detection using computer vision techniques, and automatic methods. Automatic methods have resulted in a more efficient lineament extraction process [43,[80][81][82][83]. The main steps of lineaments extraction are mentioned below: Principal Component Analysis PCA (Pearson, 1901) is a statistical method that has the advantage of compressing information contained in initial bands into new bands called principal components (PCs) [84][85][86]. This method has been commonly used in lithological mapping and lineament extraction [79]. The PCA method can reduce redundancy in different bands, which can obtain aimed dimension reduction [87], isolation of noise, and enhancement of the targeted information in the image [88]. Each PC can reflect the maximum information of the original variable, and the information contained therein is not repeated.

Filter Operations
Filter operations were used to emphasize or de-emphasize spatial frequency in the image. This frequency can be attributed to the presence of the lineaments in the area. In other words, the filtering operator can sharpen the boundary that exists between adjacent units. A directional filter is a first-order derivative edge enhancement filter that selectively enhances image features possessing specific direction components (gradients) [89]. Directional filters are used strictly for the structural analysis. These filters improve the perception of lineaments, causing an optical effect of shade worn on the image as if it was illuminated by light grazing [78].

PCI Geomatica Software
The automatic extraction of the lineaments was carried out by algorithm LINE EX-TRACTION of the PCI Geomatica software [80], which is a widely used module for automatic lineament extraction. The LINE module of PCI Geomatica software extracts linear features from an image and records polylines in vector segments by two main steps, namely edge detection and line detection, and using six parameters [43]: Edge detection

Iron Mineralization and Alteration Detection
Lowell and Guilbert (1970) described the San Manuel-Kalamazoo deposit and compared the results with 27 other porphyry-Cu deposits. According to this model, four alteration zones were introduced, which are often used for porphyry-Cu exploration. As shown in Figure 3, the zones in this model from the center to the outside are potassic, phyllic, argillic, and propylitic zones [90,91]. After presenting the Lowell-Guilbert model, some porphyry-Cu deposits were recognized to be associated with intrusive rocks posessing low silica. Hollister (1975) introduced this model and called it the diorite model, although the host pluton may be syenite, monzonite, and diorite [92]. Diorite's model differs from the Lowell-Guilbert model. In the Diorite model, sulfur concentrations were relatively low in mineralizing fluids. As a result, not all the iron oxides in the host rocks were converted to pyrite, and there are many iron remains in chlorites and biotites. Excess iron tends to occur as magnetite, which may be present in all alteration zones. Therefore, phyllic and argillic alteration zones are usually absent so that the potassic zone is surrounded by the propylitic zone [93,94].

Band Ratio
Band ratios are a very useful method for highlighting certain features or materials that cannot be seen in raw bands [95]. This method was applied to the Sentinel-2 MSI image to detect iron mineralization.

Color Composite
Colors provide more visual and conceptual information of the image. The combination of three black and white images creates a new image that can provide a better interpretation of surface features [96]. This method was used in the ASTER image for better visual interpretation of the alteration areas.
Logical Operator Algorithm Mars and Rowan (2006) developed two logical operator algorithms based on ASTERdefined band ratios for regional mappings of argillic and phyllic-altered rocks in the Zagros magmatic arc, Iran [97]. Mars (2013) used thermal images in-band ratios to map hydrothermal alterations [98]. The logical operator algorithm presented by Mars (2013) can be the best suited for hydrothermal alteration associated with porphyry-Cu mineralization on a regional scale. The logical operator algorithm performs a series of band ratios for each pixel. Each logical operator determines a true (one) or false (zero) value for each ratio by comparing the band ratio to a predetermined range of threshold values. All of the ratios in the algorithm have to be true for a value of one to be assigned to the byte image; otherwise, a zero value is produced. Thus, a byte image consisting of zeros and ones is produced with each algorithm [97]. Due to the geological settings in the study area, the logical operator algorithm was applied to the Aster image to map alterations.

Generation of The Geological Layer
The geological layer is created by combining rock units, faults, and alterations using the K-nearest neighbor (KNN) algorithm. The 1:100,000 scale map of Abrisham-Rud has complex and multistage geological settings, and some units may not have outcrop, whereas this scale does not pose a problem for this study. The layer of rock units was created by using the high value for intrusive rock and intermediate to base units, as well as the low value for other units. In the study area, faults are dense and intersect in different directions. Fault layers were created by buffers at 100 to 500 at 100 m intervals. The alteration layers for each argillic, phyllic, and propylitic (epidote-chlorite) alteration were created by buffers at 100 to 300 at 100 m intervals around the alteration.

K-Nearest Neighbor Algorithm
KNN is a non-parametrically supervised algorithm designed to solve regression and classification problems [99]. KNN is the fundamental and the simplest classification technique when there is little or no prior knowledge about the distribution of data [100]. This algorithm is quite successful when a large training data set [101] and many geological studies are provided [102][103][104][105][106].
KNN classifies objects based on the closest training examples in the feature space [107,108]. In order to classify or predict a new case, KNN relies on finding similar cases in training data. These cases are classified by voting for neighbor classes [109]. The optimal choice of the number of neighbors "K" depends on the metrics used for classification and regression purposes [109]. Thus, KNN algorithm predicts the target class through three steps [101,110]: (1) preparing the dataset consists of training, test, and feature; (2) measuring the distances between each test data and all training data depending on the weight values of each individual; and (3) finding "K" the neighbors nearest to the test data from training data based on distance and weight measurements.
The most common and simple distance metrics are Euclidean, Manhattan, and Minkowski. The Minkowski distance is generally a more complete form of distance metrics and is calculated based on Equation (9): where x and y are points to calculate the distance, k is the number of neighbors, and λ is the order of the Minkowski distance, which contains values greater than zero. Thus, where λ = 2, the Minkowski distance is equivalent to the Manhattan distance, and where λ = 1, it is equivalent to the Euclidean distance. The Manhattan distance is usually preferred over the more common Euclidean distance when there is high dimensionality in the data set [111].
In this study, the Manhattan distance was used as the nearest neighbor classifier, and weights were calculated on Equation (10) based on the distance from the target to predict in a neighborhood: where W is the weights for each nearest neighbor, 0 < W i < 1 and ∑ k i=1 Ww i = 1. Table 1 shows the calculation of background and threshold values, as well as the Clark values for supra-ore (Pb and Zn) and sub-ore (Cu and Mo) elements in the study area.  Figure 4 shows the anomaly map of multiplicative geochemical halos of the supra-ore and sub-ore (Cu × Mo) elements. As shown in Figure 4, zones I and II have been considered to predict erosional levels, and they create the geological layer. As shown in Figure 4, zone I mostly implied sub-ore (Cu and Mo) elements, and zone II included both supra-ore (Pb and Zn) and sub-ore (Cu and Mo) elements. The co-existence of both supra-ore and sub-ore elements' local maxima implies blind mineralization [20].

Zonality Method
The geochemical and geometrical similarity of genetically similar orebodies, the uniformity of ore and halo, and the tentative nature of geological and economic boundaries are all crucial for considering mining geochemical models. In addition, the vertical geochemical zonality index and their spatial associations with particular geological and geochemical factors are important aspects of mineral distributions for exploration and insight into ore geometry. The vertical geochemical zonality index could be used to estimate the erosional level of porphyry-Cu deposits [2,3]. In order to identify the erosional level in zones I and II, the presented models by Ziaii (1996) and Ziaii et al. (2009) were used [2,67]. Ziaii (1996) introduced the vertical zonality model for porphyry-Cu deposits using areal productivity and the zonality index (Equation (6)) based on porphyry-Cu deposits in Kazakhstan, Bulgaria, Armenia, and Iran (Figure 5a). The vertical variations in three zonality indices associated with porphyry-Cu deposits in areas of the same landscapegeochemical conditions in different countries are shown in Figure 5a. Values of each zonality index decrease downward uniformly despite considerable differences in local geological settings of individual porphyry-Cu deposits, suggesting the existence of uniform vertical zonality in primary halos of porphyry-Cu deposits [2,6,71]. Therefore, vertical variations in the indices allow the distinction of mineralization levels and their primary halos (supra ore, upper-ore, ore, lower-ore, and sub-ore) [6,22,24]. The values of areal productivity, mineralization coefficient, and zonality indices were introduced by Solovov (KS) and Grigorian (KG) in two zones ( Table 2). The values of (Ks) in zones I and II are equal to 0.18 and 26.57, respectively. Considering the presented model by Ziaii (1996) [71], zones I and II are ZDM and BM, respectively (Figure 5a). Table 2. Areal productivity, mineralization coefficient, and zonality index introduced by Solovov and zonality index introduced by Grigorian for Cu, Mo, Pb, and Zn elements in zones I and II. Moreover, it can be deduced from Figure 5a that similar values of the zonality index imply similar depths of mineralization and primary halos within an ore field. Thus, primary halos of mineral deposits at different depths are characterized by specific values of the zonality index. The practical exploration significance of the zonality index is for the recognition of erosional surfaces representing vertical levels of geochemical anomalies. Concerning the present erosional level, high values of the zonality index imply the presence of sub-cropping to BM, whereas low values of the index imply outcropping or already eroded deposits [14].
The values of areal productivity, mineralization coefficient, and zonality indices were introduced by Solovov (K S ) and Grigorian (K G ) in two zones ( Table 2). The values of (Ks) in zones I and II are equal to 0.18 and 26.57, respectively. Considering the presented model by Ziaii (1996) [71], zones I and II are ZDM and BM, respectively (Figure 5a). Table 2. Areal productivity, mineralization coefficient, and zonality index introduced by Solovov and zonality index introduced by Grigorian for Cu, Mo, Pb, and Zn elements in zones I and II.

Zone
Elements P, m 2 % η(α) K S K G Ziaii et al. (2009) introduced the geochemical model for porphyry-Cu deposits in Aktogy (Kazakhstan), Asarel (Bulgaria), Tekhut (Armenia), and Sungun (Iran) using mineralization coefficient and zonality index (Figure 5b) [2]. This plot shows the depth of mineralization versus the zonality index (KG). Despite considerable differences in geological settings, the linear relationship suggests the existence of a quantitatively uniform vertical geochemical zonality in the structure of primary halos of the deposits.
In this study, the values of (KG) in zones I and II are equal to 0.71 and 28.23, respectively ( Table 2). According to the presented model by Ziaii et al. (2009), erosional levels in zone I and II are nearly similar to Astamal and Songun 2 areas, respectively (Figure 5b). Based on the previous studies, Sungun 2 and Astamal areas have been recognized as BM and ZDM, respectively [2]. Therefore, the results of the introduced models by Ziaii (1996) and Ziaii et al. (2009) confirm each other in identifying the erosional level in each zone [2,71].

Remote Sensing Lineaments Extraction
The procedure of lineament extraction was accomplished in this manner: the PCA image of six Sentinel-2 MSI bands (bands 2, 3, 4, 8, 11, and 12), as shown in Figure 6a. PC1 explains the largest amount of eigenvalue among six bands. PC1 with the loading of the same signs represents overall brightness in all bands [14], and it shows that the albedo is related largely to the topographic features [76] (Figure 6b). Then, a directional filter was applied using 3 × 3 kernels in four directions with azimuths of 0 • , 45 • , 90 • , and 135 • . By using these azimuths, this filter visually enhances edges striking N-S, NE-SW, E-W, and NW-SE, respectively.
The sum of the directional filter kernel arrays is zero. The result is that areas with uniform pixel values are zeroed in the output image, while those with variable pixel values are presented as bright edges. In PCI Geomatica software, differences in the values of six parameters of edge and line detection indicate the differences of opinion among researchers. The values proposed by Adiri et al. (2017) were used [43]. Finally, the lineaments obtained from azimuths of 0 • , 45 • , 90 • , and 135 • were merged, and the repetitive segments and non-geological lineaments (river, road, etc.) were deleted.
Lineaments and lineaments density maps in two zones were demonstrated in Figure 7. Lineament density was used to find the correlation between the concentration of lineaments and the distribution of existing faults in the study area (Figure 7). In two zones, this comparison proved that the fault is well related to lineament density in most areas. These areas were generally recognized in the middle-upper Eocene units and intrusive rocks. Lineament orientation allows identifying the most frequent directions of lineaments, and they can be compared with directions related to the existing faults [43]. As shown in Figure 8, in zones I and II, the directions of the lineaments correspond to existing faults.
Lineaments may be formed, for example, by structural alignment, geomorphologic consequences, structural weaknesses, faults, valleys, rivers, the boundaries between the different lithological units, vegetation cover, and artificial objects (road, bridge, etc.) [43]. In this study, due to the importance of faults in porphyry-Cu mineralization [6], the faults obtained from the geological map were used to create a geological layer.

Iron Mineralization and Alteration Detection
The results of the band ratio applied to Sentinel-2 MSI were shown in Figure 9. According to the rock units in the studied zones and Porphyry-Cu alterations zones, iron mineralization is dispersed and dense in two zones (Figure 9), which is well identified in areas containing volcanic and intrusive rocks. In the south of zone I, which includes intrusive rock consisting of monzodiorite and monzogabbro, iron mineralization is well recognized.
In ASTER, false-color composites of SWIR bands were used for better visual interpretation of the alteration areas. Empirical combinations have shown that an image with a false-color composite (band 4 in red, band 6 in green, and band 8 in blue) is the most suitable color composite for identifying alteration areas in porphyry-Cu deposits. As shown in Figure 10, areas with the propylitic alteration are shown in green to dark green based on the alteration intensity, and areas with the argillic and phyllic alteration are shown in white and pink to red. This is due to the high reflectivity of alunite, kaolinite, and muscovite minerals in band 4 compared to bands 6 and 8.
Geological settings in the study area have made it difficult to identify some alterations, especially in zone II. Rocks containing hydrous quartz, chalcedony, opal, and amorphous silica (hydrothermal silica-rich rocks); calcite-dolomite and epidote-chlorite (propylitic); alunite-pyrophyllite-kaolinite (argillic); and sericite-muscovite (phyllic) were mapped using ASTER and logical operator algorithms ( Figure 11). It is observed that Mars (2013) used images that have to differ in correction levels from the images used in this paper. Mars (2013) used "and (b4 gt 260)" to remove the black pixel [94], but in the images that have been used in this paper, the pixel has no value higher than 260. Thus, "and" in the algorithm causes the result to be zero (Table 3). Although not all alterations are associated with ore bodies and not all ore bodies are accompanied by alteration, the presence of altered rocks is a valuable indicator of possible deposits [112]. Kaolinite is mostly related to weathering feldspars, and epidote can be related to regional metamorphism. Kaolinite and epidote anomalies can have genesis related to deposit when they have a close relationship with muscovite anomalies [65].
Due to intrusive rocks, there is a possibility of the diorite model in the south of zone I. In this area, propylitic alteration (epidote-chlorite) was identified less than other parts of zone I, and argillic and phyllic alterations are well-identified around these intrusive rocks. It should be noted that all intrusive rocks have not shown alterations, which could be due to erosion or geological settings. Furthermore, alterations have been identified in other areas of zone I, which could imply mineralization areas or presence of minerals associated with alterations. In zone II, propylitic alteration (epidote-chlorite) was not identified, but argillic and phyllic alterations were detected in part of the studied area. Therefore, because of the absence of the alterations, it is not reasonable to create a geological layer in this zone.

Geological Layer
Rock units, faults, and alterations layers were combined by using the KNN algorithm and the geological layer, as shown in Figure 12. In this procedure, faults and alterations have an important role in mapping high potential areas. The density and intersection of faults and the extent of alterations represent these areas. Considering the geological layer, the detected areas as high potential could be related to mineralization or not. These areas can be compared to anomalous areas obtained by using the zonality method. As a result, more parts of the geological layer map were detected as high potential in comparison to the zonality method. However, the results of both methods confirm each other in the south of zone I. The geological layer and mineralogical investigation cannot identify erosional levels. In addition, mineralogy is not economical for application on a regional scale. In the study area, it is not recommended.    Table 3. The logical operator algorithms used with ASTER to map hydrothermally altered rocks in zones I and II (modified after Reference [98]). (b: band; float: floating point; le: less than or equal to; lt: less than; ge: greater than or equal to; gt: greater than.)

Zone Hydrothermal Alteration Algorithm
Zone I Hydrothermal silica-rich (hydrous silica, chalcedony, opal) Propylitic Unlike the geological layer, the application of the zonality method in mineral prospecting allows further interpretation about whether delineated desirable areas are attractive for the exploration of ZDM or BM deposits. This comparison demonstrates that the zonality method for detecting anomalous areas is more powerful than the geological layer.

Conclusions
The traditional zonality method has been used in the exploration of porphyry-Cu deposits for many years and is an effective method for the distinction between sub-ore and supra-ore halos, prediction of the erosional level of mineralization, and exploration of blind mineral deposits.
Utilizing the zonality method, the geochemical maps of multiplicative haloes were mapped. In the east of the study area, multiplicative haloes of sub-ore elements (Cu and Mo) were observed, and these represent zone dispersed mineralization. In the northwest, both multiplicative supra-ore (Pb and Zn) and sub-ore (Cu and Mo) element haloes existed, and these imply blind mineralization. Thus, zones I and II, which are located in the east and northwest of the study area, were selected for calculating erosional levels and for creating the geological map. Zonality indices introduced by Solovov (1987) (Ks) and Beus and Grigorian (1977) (KG) were calculated in zones I and II [12,22]. The (Ks) values were equal to 0.18 and 26.57 in zones I and II, respectively; moreover, (KG) values in zones I and II were equal to 0.71 and 28.23, respectively. Then, the presented models by Ziaii (1996) and Ziaii et al. (2009) were considered for identifying erosional levels in these zones [2,67]. Due to these models, zones I and II were recognized as ZDM and BM, respectively. Therefore, the zonality method was successfully applied in the identification of anomalous areas, separate BM from ZDM, and predicted erosional levels.
The results of the zonality method were compared to the geological layer, which was created by rock units, faults, and alterations by using the KNN algorithm. Each of these layers plays an important role in prospecting and exploring mineral deposits. Thus, high potential areas can be identified by combining these layers. For zones I and II, rock units and faults were identified from the geological map, and the alterations were detected using ASTER images and logical operator algorithms. It was observed that the alterations layers had a significant contribution in constructing the geological layer. The alteration zones of porphyry-Cu deposits include propylitic (chlorite and epidote), argillic (alunite and kaolinite), and phyllic (sericite and muscovite), and they are important for identifying possible areas associated with porphyry-Cu systems. These alterations were detected in zone I, especially around the intrusive rocks in the S of this zone. In zone II, only argillic and phyllic alterations were identified in part of the studied area. Due to the lack of alteration in zone II, the geological layer was created only in zone I. Comparing the results of the two methods showed that more parts of the geological layer map were highlighted as having high potential. These high potential areas could be related to mineralization or not; in other words, the geological layer cannot separate BM from ZDM. However, the results of both methods correspond to each other in the south of zone I. In other, the geological layer is unable to recognize erosional levels. Therefore, mineralogy investigation is required, which is not recommended to apply on a regional scale because of its high cost. It could be concluded that the geological layer, which is based on alteration, cannot help geochemists in separating BM from ZDM and in predicting erosional levels. Data Availability Statement: All information is collected by the researchers of this article and is not accessible anywhere.