Shear-Related Gold Ores in the Wadi Hodein Shear Belt, South Eastern Desert of Egypt: Analysis of Remote Sensing, Field and Structural Data

: Space-borne multispectral and radar data were used to comprehensively map geological contacts, lithologies and structural elements controlling gold-bearing quartz veins in the Wadi Hodein area in Egypt. In this study, enhancement algorithms, band combinations, band math (BM), Principal Component Analysis (PCA), decorrelation stretch and mineralogical indices were applied to Landsat-8 OLI, ASTER and ALOS PALSAR following a pre-designed ﬂow chart. Together with the ﬁeld observations, the results of the image processing techniques were exported to the GIS environment and subsequently fused to generate a potentiality map. The Wadi Hodein shear belt is a ductile shear corridor developed in response to non-coaxial convergence and northward escape tectonics that accompanied the ﬁnal stages of terrane accretion and cratonization (~680–600 Ma) in the northern part of the Arabian–Nubian Shield. The evolution of this shear belt encompassed a protracted ~E–W shortening and recurrent sinistral transpression as manifested by east-dipping thrusts and high-angle reverse shear zones. Gold-mineralized shear zones cut heterogeneously deformed ophiolites and metavolcaniclastic rocks and attenuate in and around granodioritic intrusions. The gold mineralization event was evidently epigenetic in the metamorphic rocks and was likely attributed to rejuvenated tectonism and circulation of hot ﬂuids during transpressional deformation. The superposition of the NW–SE folds by NNW-trending, kilometer scale tight and reclined folds shaped the overall framework of the Wadi Hodein belt. Shallow NNW- or SSE-plunging mineral and stretching lineations on steeply dipping shear planes depict a considerable simple shear component. The results of image processing complying with ﬁeld observations and structural analysis suggest that the coincidence of shear zones, hydrothermal alteration and crosscutting dikes in the study area could be considered as a model criterion in exploration for new gold targets.


Introduction
Remote sensing data are commonly used to identify geological structures associated with extensive shear belts, particularly in vast and rugged terrains. The innovative technologies in the remote sensing promoted the application of space-borne imagery data for detailed structural geology mapping [1][2][3][4][5][6]. The application of multisensor satellite imagery can be considered as a cost-efficient exploration strategy for prospecting orogenic gold mineralization in transpression and transtension zones, which are located in harsh regions around the world [7,8]. In Egypt, structural analysis aided by the interpretations of spaceborne data has been efficiently employed to constrain the major structures in the South Figure 1. Distribution of the known gold mining sites and gold quartz vein occurrences relative to the major fault/shear structures and ophiolitic masses in the South Eastern Desert terrane (modified after [11]). Inset shows a location map of the study area.
Gold bearing quartz veins and their alteration zones occur in different geologic and structural settings in the southern Eastern Desert [10] (Figure 1). These different settings include: (1) Au-quartz veins hosted by ductile shear zones at contacts between ophiolitic and island arc terranes (e.g., Um El-Tuyor El-Foqani, Betam, Seiga, Shashoba, Um Garayat, and Haimur deposits); (2) Au-quartz veins occurring along steeply dipping anastomosing shear zones wrapped around or cutting syn-or late-orogenic granitoid intrusions (e.g., Korbiai, Madari, Romite, Egat deposits); (3) Au-quartz-carbonate veins in association with brittle-ductile shear zones separating listvenized serpentinite from the underlying successions of intercalated metavolcanic and volcaniclastic metasediments (e.g., Hutit, El-Beida, and El-Anbat deposits). In the present study, structural analysis supported by remote sensing studies was aimed at revealing what controls the distribution of the scattered goldquartz vein occurrences in the Wadi Hodein area. Multisensor satellite data were used to resolve the spatial relationship between gold mineralization and district-scale structures. The relative timing of the successive structural events was also envisaged to place the gold mineralization in the broader evolutionary framework of the South Eastern Desert.

Geological Setting
The Neoproterozoic rocks exposed in Hodein shear belt comprises variably deformed ophiolitic mélange rocks, island-arc metavolcanic and metavolcaniclastic successions ( Figure 2). These rocks are unconformably overlain by Cretaceous sandstone formations (Nubian Sandstone) and locally by Miocene marine turbidites and carbonate. Dioritic and granitic gneisses and migmatites are exposed around Wadi Beitan and Wadi Khuda, in the western part of the study area. Wadi Khuda gneisses occur as an ENE-WSW trending belt forming a major doubly plunging anticline ( Figure S1a,b), amphibolite and migmatites [35]. It is intruded by syn-tectonic tonalite and diorite as well as late-to post-tectonic syenogranites and leucogranites. The gneissic rocks of Wadi Beitan form a major NW-trending belt in the northwestern part of the Wadi Hodein shear system [27] (Figure 2). This belt com-prises fine-to coarse-grained biotite gneisses, biotite-hornblende gneisses and hornblende gneisses together with subordinate varieties of garnet-biotite-hornblende gneisses, augen gneiss, amphibolites (Figure S1c,d), migmatites and mylonites [11,27]. The contact between the Wadi Beitan gneisses belt and the volcaniclastic rocks is defined by a SW-verging thrust fault associated with discrete zones of mylonitized granites ( Figure S1d). trending belt forming a major doubly plunging anticline ( Figure S1a,b), amphiboli migmatites [35]. It is intruded by syn-tectonic tonalite and diorite as well as late-to tectonic syenogranites and leucogranites. The gneissic rocks of Wadi Beitan form a NW-trending belt in the northwestern part of the Wadi Hodein shear system [27] ( 2). This belt comprises fine-to coarse-grained biotite gneisses, biotite-hornblende gn and hornblende gneisses together with subordinate varieties of garnet-biotiteblende gneisses, augen gneiss, amphibolites (Figure S1c,d), migmatites and myl [11,27]. The contact between the Wadi Beitan gneisses belt and the volcaniclastic ro defined by a SW-verging thrust fault associated with discrete zones of mylonitized ites ( Figure S1d).  [36][37][38][39]).
The ophiolitic mafic metavolcanics to northeast of Gabal Abu Dahr are massive or pillowed ( Figure S2a) metabasalt and metabasaltic andesite. Pillowed metabasalt bodies are exposed as ellipsoidal, globular and tabular shaped bodies, ranging between 10 cm and 1 m in diameter ( Figure S2a,b). The mafic metavolcanics along Wadi Hutib and Wadi Urga El Rayan form a NNW-elongated belt and show variable degrees of shearing. The Wadi Arais ophiolites include masses of serpentinite, serpentinized peridotites ( Figure S2c), metagabbro, dunite, pyroxenite, chlorite schist, metavolcanics rocks and listwaenite along the thrust and shear zones [43,44]. The serpentinized peridotite occurs as steep mountains along NW-SE striking thrust faults ( Figure S2c). Irregular serpentinite masses are incorporated in the volcaniclastic metasediments. The ophiolitic metagabbro rocks form hills of moderate relief, which have tectonic contacts against the serpentinite and metavolcanic blocks. The Gabal Arais mass is composed chiefly of actinolite hornblende metagabbro and hornblende metagabbro and less commonly of appinitic metagabbro [45]. The arc metavolcanic rocks along Wadi Um Araka form an E-W elongate belt comprised mainly meta-andesite, metadacite, dacitic tuff and andesitic tuff with minor amounts of metabasalt [45]. The foliated varieties of the Um Araka metavolcanics are mainly schistose metatuffs. The ophiolitic mélange of Wadi Arais comprises volcaniclastic metasediments and exotic clasts of metagabbro and amphibolite.
The metagabbro-diorite-tonalite-granodiorite complex forms a large elongated belt in the Wadi Khuda area. These rocks intrude Wadi Khuda gneisses ( Figure S2d) and the eastern parts of Wadi Rahaba-Gabal Abu Dahr ophiolitic mélange and intruded by late to post-tectonic biotite and muscovite granites ( Figure S2e). Their contacts with Wadi Khuda gneisses are characterized by the strongly foliated diorite and amphibolite [46,47]. The gabbro-diorite complex forms variable degrees of deformation, from massive to strongly sheared or foliated. The metagabbro-diorite rocks of the Um Eleiga exhibit circular zonation where gabbro forms the core and the granodiorite occupies the margin [41].
The tonalite of Wadi Urga El Rayan form N-S trending, strongly weathered and weakly foliated plutons in the central part of the study area. Granodiorite forms a N-S trending pluton along the Wadi Hutib and E-W oriented small intrusions along the Wadi Beitan ( Figure S2f). The granitoid intrusions that cut the ophiolitic mélange of the Wadi Khashab-Wadi Um Karaba range in composition from tonalite to granodiorite. The granitic rocks of Gabal El Maarafay and Gabal Al Farayid include monzogranite, granodiorite, and tonalite and are locally foliated [48]. At Wadi Khuda, leucogranite forms small stocklike and dikelike bodies intruding on the core of the Khuda gneisses and the granite pluton [46]. The ophiolitic rocks of Wadi Arais are cut by monzogranite and syenogranite of the Gabal El-Nukeiba and Gabal Handusa. A small circular intrusion of younger gabbro in the southwestern part of the study area (G. Homraii) cuts metavolcanic rocks and the older granitoid terranes. The gold-bearing quartz veins in the study area cut sheared metavolcanics, volcaniclastic metasediments and gabbro-diorite complex. The mineralized quartz veins are controlled consistently by NW-and NNW-trending shear zones, and are associated with extensive pervasive hydrothermal alteration. Gold occurrences in the study area include El-Beida, Khashab, El-Anbat, Um Teneidab, Urga Ryan, Hutit and Um Eleiga.

Remote Sensing Data Characteristics and Analysis
Iron oxides, clay and carbonate ± sulfate in the hydrothermal alteration zones have specific spectral signatures in the visible near infrared (VNIR), and shortwave infrared (SWIR) spectral regions, respectively [13,14]. False color composite (FCC) and band combination images are used in geological applications based on known spectral signatures of mineral phases in specific wavelength regions. Mapping the hydrothermal alteration zones is a prime focus of mineral exploration programs using the remote-sensing data [1][2][3][4][5][6][7][8]11,12]. Most of the hydrothermal alteration mineral species have distinctive features in the shortwave infrared (SWIR) region, making the Landsat-8 OLI and ASTER sensors suitable data sources for mineral mapping, particularly in arid regions. Moreover, considerable numbers of image processing techniques, such as band-rationing (BR), principal component analysis (PCA), and the spectral mineralogical indices, have been proven effective in lithological and hydrothermal alteration mapping if verified by the fieldwork [1,7,8,11,12]. The present study integrates Landsat-8 OLI/TIRS, ASTER and ALOS PALSAR data for comprehensive mapping of the lithological units and geological structures. Additionally, we aimed to detect the mineralized zones in the Wadi Hodein area. The characteristics of the used remote sensing data are presented in Tables 1 and 2.

Data and Processing Techniques
Two cloud-free ASTER scenes AST_L1T00312252006082430 and 00312252006082422 were acquired on 24 August 2006. The Landsat-8 OLIL/TIRS scene 'LC08_L1TP1730442019073001_T1 was acquired on 30 July 2019 with path/raw 173/44. Four ALOS PALSAR scenes, ALP-SRP075090450, ALPSRP075090460, ALPSRP077570450 and ALPSRP077570460, of L-Band level 1.5 images. For detailed information on the preprocessing techniques of the remote sensing data, adopted software and methodology see the supplementary data file (Section S1 in supplementary material). The remote sensing data have been processed for lithological mapping, highlighting structural elements and delineating the alteration zones. Image processing techniques, i.e., enhancement algorithms, band combinations (FCC), band math (BM), Principal Component Analysis (PCA), decorrelation stretch and mineralogical indices were applied for outlining geological mapping and the hydrothermal alteration zones in the study area. Finally, multicriteria approach was applied to produce the potentiality gold mineralization map of the study area. Figure S3 shows the flow chart of the methodology adopted for the present study.

Landsat-8 OLI and ASTER-Based Lithological Mapping
The Landsat-8 OLI data have higher radiometric resolutions (16 bits) and lower spectral resolution compared with ASTER data. Several Landsat-8 OLI and ASTER band combinations images are employed for mapping lithological units and structural elements in the study area. The Landsat-8 band combination images (RGB-753; Figure 3a) and ASTER (RGB-431; Figure S4a) discriminate well between different lithologic units and highlight the structural elements. The ophiolitic rocks appear as dark green pixels, metavolcanics and metasediments are represented by dark to medium brown pixels while the highly weathered granitoids and gneisses exhibit light brown signature. The E-W, NW-, NEtrending faults/fractures that cut these rock outcrops are conspicuous. The band math technique, which is used for reducing the effects of topography and enhancing the spectral differences between bands, divides the digital number value of one band by the digital number value of another band [49]. Band math is a useful for highlighting hydrothermal alteration minerals, as it minimizes the illumination differences caused by the topographic features. Landsat-8 band math image (6/7) was generated to map CO 3 and OH-bearing minerals. This image ratio highlights clay minerals, serpentine, and many alteration zones as bright pixels, where the metavolcanic-metasedimentary rocks are manifested by gray pixels (Figure 3b). On the other hand, the grey-scale image of ASTER band math (7 + 9/8) presented as carbonate/chlorite index. It highlights the highly sheared and carbonatized rocks (serpentinite and talc-carbonate) and metavolcanics with bright image signature ( Figure S4b). Ophiolitic gabbros and metasedimentary rocks appear as medium grey pixels, while the gneisses and granitoid rocks have black spectral signatures. Hassan and Sadek [15] used this ratio to differentiate between talc-carbonate and ophiolitic basic metavolcanic rocks in the Gerf area south of the present study area.
Segal [50] used band ratio images for the enhancement of spectral contrasts among the bands considered in the ratio operation, and they have successfully been used in the mapping of alteration zones. Therefore, Landsat-8 FCC of the Abrams ratio (RGB-6/7, 4/3, 5/4), Chica-Olma ratio (RGB-6/7, 6/5, 4/2) and Kaufmann ratio (RGB-7/5, 5/4, 6/7) (Figure 4a-c, respectively) were generated to characterize iron-oxide, hydroxyl minerals, clay minerals and containing rocks. The Abrams ratio (RGB-6/7, 4/3, 5/4) image successfully delineated mafic mineral-rich rocks, including the ophiolitic serpentinite, and highly tectonized ophiolites, which are clearly demarcated as bright red color, and the highly tectonized mafic metavolcanic rocks appear as purple and dark blue pixels ( Figure 4a). The Chica-Olma ratio (RGB-6/7, 6/5, 4/2; Figure 4b) image clearly delineates the serpentinite rocks and their sharp contacts with the surrounding rocks of metavolcanicmetasedimentary association and granitoids, where theses contacts are considered the favorable sites for gold mineralization. Serpentinite appear as bloody red pixels, metasedimentary rocks as dark green, whereas the metavolcanics and quartz-carbonate (listvenite) have a bright green signature. The altered materials and clay minerals appeared as olive green pixels. The ophiolite assemblage is highlighted successfully as bright blue and quartz-carbonate (listvenite) as violet pixels in the Kaufmann ratio (RGB-7/5, 5/4, 6/7; Segal [50] used band ratio images for the enhancement of spectral contrasts among the bands considered in the ratio operation, and they have successfully been used in the mapping of alteration zones. Therefore, Landsat-8 FCC of the Abrams ratio (RGB-6/7, 4/3, 5/4), Chica-Olma ratio (RGB-6/7, 6/5, 4/2) and Kaufmann ratio (RGB-7/5, 5/4, 6/7) ( Figure  4a-c, respectively) were generated to characterize iron-oxide, hydroxyl minerals, clay minerals and containing rocks. The Abrams ratio (RGB-6/7, 4/3, 5/4) image successfully delineated mafic mineral-rich rocks, including the ophiolitic serpentinite, and highly tectonized ophiolites, which are clearly demarcated as bright red color, and the highly tectonized mafic metavolcanic rocks appear as purple and dark blue pixels (Figure 4a). The Chica-Olma ratio (RGB-6/7, 6/5, 4/2; Figure 4b) image clearly delineates the serpentinite rocks and their sharp contacts with the surrounding rocks of metavolcanic-metasedimentary association and granitoids, where theses contacts are considered the favorable sites for gold mineralization. Serpentinite appear as bloody red pixels, metasedimentary rocks as dark green, whereas the metavolcanics and quartz-carbonate (listvenite) have a bright green signature. The altered materials and clay minerals appeared as olive green pixels. The ophiolite assemblage is highlighted successfully as bright blue and quartz-carbonate (listvenite) as violet pixels in the Kaufmann ratio (RGB-7/5, 5/4, 6/7; Figure 4c) image. The metasedimentary rocks are characterized by dark cyan color, whereas the metavolcanics appeared by reddish brown pixels. These distinctive spectral signatures helped in heightening the foliation trajectories of the highly deformed and tectonized rocks. The NW-thrusted contacts and E-W strike slip faults can be detected from the three image ratios (Figure 4a-c). Moreover, decorrelation stretching is an important image enhancement method that is used to improve visual interpretation of satellite images [51]. For bands having the maximum variance, decorrelation stretching was applied to 7-6-5 of Landsat-8 bands (Figure 4d). The dark blue pixels of serpentinite can be differentiated from the greenish brown pixels of ophiolitc metagabbros. The metasedimentary rocks and metavolcanics attained greenish red and bright green colors, respectively. Kaolinite-and illite-bearing rocks show bright red signature, whereas the gneissic and granitoid rocks appear as light violet pixels (Figure 4d).
Principal Component Analysis (PCA), as one of the spectral enhancement techniques, has been used for lithological discrimination [15]. Principal component analysis (PCA) is a multivariate statistical technique that selects uncorrelated linear combinations (eigenvector loadings) of variables in such a way that each successively tracted linear combination, or principal component (PC), has a smaller variance [52]. The PCA transformation was carried out for the VNIR and SWIR bands of Landsat-8 and ASTER to extract lithological and alteration zones information, which related to gold mineralization. The eigenvector matrices of the Landsat-8 and ASTER data derived from the PCA are given in Tables S1 and S2, respectively). The three first PCA images (PC1, PC2 and PC3), containing the highest topographical and spectral information, are found suitable for lithological discrimination. The PCA images can bear crucial information related to the alteration minerals, which could be reflected in the eigenvector loading of the absorption and reflection bands [12]. Gupta et al. [53] suggested that when the eigenvector loading is strong (as positive or negative sign) in the reflection and absorption bands of the target mineral or mineral group, the enhanced pixels related to the mineral or mineral group will manifest as bright pixels when this loading is positive in the reflection band; conversely, the pixels will manifest as dark when the loading is negative in the reflection band.  The NW-thrusted contacts and E-W strike slip faults can be detected from the three image ratios (Figure 4a-c). Moreover, decorrelation stretching is an important image enhancement method that is used to improve visual interpretation of satellite images [51]. According to the analysis of the eigenvector loadings of the seven bands of Landsat-8 PCA (Table S1), we produced the color composite images (RGB-PC2, PC1, PC4 and RGB-PC2, PC4, PC5; Figure S5a,b, respectively). Bands 2, 4, 5, 6, 7 of Landsat-8 contain valuable information for iron oxides and hydroxyl-bearing minerals' mapping [54]. The rocks consisting of high contents of Al and/or Mg-OH-bearing minerals (ophiolites) were best delineated by bright lemon-yellow pixel signatures in Figure S5a and yellowish orange pixels in Figure 5a. The two images emphasize the contact between felsic and basic metavolcanics and volcaniclastic metasediments and gneisses along Wadi Hodein and Wadi Beitan. The areas with rich altered materials show cyan color tones. Based on the eigenvector loadings of the nine bands of ASTER PCA (Table S2), the color composite images (RGB-PC1, PC2, PC3 and RGB-PC6, PC3, PC1; Figure S5a,b, respectively) were prepared. Figure  S5b shows bright magenta pixels for iron oxide/hydroxide minerals, which are mostly associated with ophiolitic serpentinite and metagabbros. The metasedimentary rocks are manifested by light magenta pixels, metavolcanics as bright green, and gneisses and granitoids as cyan pixels. Figure 5b shows the metavolcanic/metavolcaniclastic rocks appearing as reddish color pixels. The light cyan color domains within the ophiolitic mélange terranes whilst metasediments represent grass green and gneisses and granitoids are presented by yellow green pixels. The four images enhance the visualization of the E-W strike slip faults and foliation trajectories.

Hydrothermal Alteration Zones Detection
Shortwave infrared (SWIR) channels of ASTER can increase the accuracy of the spectral identification of minerals and rock units [55][56][57][58][59][60][61]. ASTER indices including kaolinite, clay and muscovite and carbonate alteration zones [12]. The mineralogical indices are reflectance combinations of two or more spectral bands signifying the relative abundance of the target objects. The application of mineralogical indices for the ASTER bands is efficient for lithological and hydrothermal alteration discrimination [62]. The six ASTER-SWIR bands were used to map the hydrothermal alteration minerals. Therefore, four spectral mineralogical indices (OH-bearing mineral index, OHI; kaolinite index, KLI; calcite index, CLI and alunite index, ALI) were used to characterize the alteration zone in the study area. The OHI is calculated as (band 7/band 6) × (band 4/band 6), the KLI is developed by way of (band 4/band 5) × (band 8/band 6), CLI is calculated as (band 6/band 8) ×

Hydrothermal Alteration Zones Detection
Shortwave infrared (SWIR) channels of ASTER can increase the accuracy of the spectral identification of minerals and rock units [55][56][57][58][59][60][61]. ASTER indices including kaolinite, clay and muscovite and carbonate alteration zones [12]. The mineralogical indices are reflectance combinations of two or more spectral bands signifying the relative abundance of the target objects. The application of mineralogical indices for the ASTER bands is efficient for lithological and hydrothermal alteration discrimination [62]. The six ASTER-SWIR bands were used to map the hydrothermal alteration minerals. Therefore, four spectral mineralogical indices (OH-bearing mineral index, OHI; kaolinite index, KLI; calcite index, CLI and alunite index, ALI) were used to characterize the alteration zone in the study area. The OHI is calculated as (band 7/band 6) × (band 4/band 6), the KLI is developed by way of (band 4/band 5) × (band 8/band 6), CLI is calculated as (band 6/band 8) × (band 9/band 8), while ALI is expressed as (band 7/band 5) × (band 7/band 8) [62]. The hydrothermal alteration zones are traced as bright pixels in the highly tectonized ophiolites, metasedimentary and island arc metavolcanic and metavolcaniclastic rocks (Figure 6a-d). The island arc rocks appear as bright pixels likely due to the abundant OH-minerals such as mica, amphiboles, chlorite and epidote ( Figure 6a). It is noticed that the Urga, Um Teneidab, Khashab-1, Khashab-2 and El-Anbat gold mines overlap these alteration zones. Kaolinitic and clay minerals rich rocks such as granitoids, wadi deposits and altered ophiolitic rocks show bright pixels in Figure 6b, where Um Eleiga, Um Karaba and Hutit gold mines are seen. The calcite index grey image discriminates perfectly the calcite rich rocks zone of serpentinite, talc carbonate rocks, island arc metavolcanics and highly sheared metasediments as bright and gray image pixels.
The false-color composite (FCC) ratio image (OHI, KLI, CLI) was elaborated by combination of the three gray scale mineralogical indices images in the RGB channels ( Figure 6d). This image characterizes the OH-bearing rocks (meta-ultramafic) as light blue pixels, metavolcanics and metasediments as dark blue and reddish blue pixels, respectively. The quartz-rich rocks have a reddish yellow image signature, whereas the clay minerals bearing rocks and stream sediments have lemon and bright green image signatures. The hydrothermal alteration zones in the study area are associated with the highly deformed ophiolitic and the ductile deformation zone. The results of alteration mapping derived from ASTER datasets reveal the spatial association of highly sheared rocks with gold occurrences.

PALSAR-Based Lineaments Extraction
Lineaments' extraction and analysis is considered a fast tool in geologic mapping and mineral exploration [63]. Hydrothermal ore deposits occur within or nearby fracture zones. Thus, lineament mapping is useful to infer locations of unexposed mineralization or dilatation sites [64]. High-resolution PALSAR data with full polarization and variable offnadir angle have significantly enhanced the structural mapping and lineaments extraction. The adaptive Lee and Local Sigma filters can be applied to PALSAR to eliminate speckles and to enhance structural lineaments [6,65]. In the present study, four PALSAR scenes were orthorectified, mosaicked and enhanced by Lee filter to trace the linear features and major structures. The HH, HV, HH + HV bands were transformed by PCA method to PCA RGB images (PC1, PC2, PC3; Figure 7a). The automated lineament extraction was carried out using PC1 image ( Figure 7b) by using LINE algorithm PCI Geomatica software. This algorithm is adopted with parameters, including radius of filter in pixels (RADI), threshold for edge gradient (GTHR), threshold for curve length (LTHR), threshold for line fitting error (FTHR), threshold for angular difference (ATHR), and threshold for linking distance (DTHR). Minerals 2021, 11, x FOR PEER REVIEW 12 of 39   The extracted lineaments were overlaid on the PC1 image ( Figure 7c) and exported to Rockworks software to generate the azimuth-frequency diagram (inset, Figure 7c). The lineaments layer was imported in the ArcMap (version 10.5, ESRI (Environmental Systems Research Institute, Redlands, CA, USA) environment to generate the lineament density map using the line density module in the spatial analyst toolbox (Figure 7d). The resultant lineaments layer shows predominantly NW-SE, N-S and E-W trending structural grain, which agrees mostly with the major thrusts and faults in Figure 2. Well-developed foliations, shear cleavages and strike slip faults offsets are also recognized (Figure 7a,b). The gold mining sites occur close to the high-density zones of the lineaments (Figure 7d), similar to conclusions reached by Zoheir and Emam [7], and Zoheir et al. [11,12] in the south Eastern Desert of Egypt.

Geospatial Modelling for Gold Mineralization Zones
Geographic Information System (GIS) is capable of manipulating, overlaying, integrating, and storing digital data in the geo-database as thematic layers. This approach plays an essential role in creating and integrating numerous geo-datasets from different geospatial data such as geology, structural density, etc. [66]. Remote sensing and (GIS) were integrated to identify new gold mineralization sites. For geospatial modelling of these mineralization zones, the evidence layers were integrated according to their relative contribution to the gold mineralization [67]. Where three main processes are recommended, including ranking processes for the evidence layers, weighting process are used for the classes of each layer, and data integration is used to assess gold mineralization potential zones [67,68]. Multisource datasets were obtained, analyzed and integrated to develop a geospatial model. Seven thematic maps for the study area were prepared and then converted into raster or vector form to be easily integrated with the GIS tools. A weight was assigned to each factor layer according to the different ranking methods [69][70][71][72][73] (Table 3). The weighted sum of all the factor layers was calculated ( Figure S6). The layers were then processed using ArcGIS analysis tools to place the assigned gold occurrences on a knowledge-based hierarchy, utilizing the Spatial Analyst tool of ArcGIS ( Figure S6). To perform this process three basic steps were followed: spatial database building, data analysis, and data integration. The created layers for the geospatial model were subjected to reclassification and assigned suitable weights ( Figure S6 and Table 3). Moreover, the geospatial thematic maps of the gold mineralization model were ranked from 1 to 5 (where 5 is the most favorable and 1 is the least) based on their suitability to host gold mineralization.
Gold mineralization is associated with the ophiolitic rocks and their contacts, close to contacts, and mostly along shear zones. A mineralization potential map was created and classified the study area into five relatively zones, varying from very low to very high potential zones (Figure 8). The high potential zones are distributed mainly around ophiolitic rocks contacts (serpentinite, talc carbonate and ophiolitic metagabbros) and metavolcanic and metasediments association rocks, respectively. The new promising sites occur close to the already known gold mining sites ( Figure 8). The integrated approach demarcates several potential zones along the NW-Hodein shear zone. Validation of this model was accomplished by testing ancient mining localities in the study area. Figure 8 shows that Um Eleiga, Um Karaba, Hutit and Urga occur in the very high potential zone, while Um Teneidab, Khashab 1, 2 and El-Anbat match the high potential zone. Additionally, as constrained by new fieldwork and the available literature [11,12], occurrences of gold-bearing quartz veins are mostly confined to thrust ophiolite/metavolcanic-metasediments contacts.

Structural Setting and Analysis
The Wadi Hodein shear belt ( Figure 9) is a high angle NW-oriented transcurrent shear zone in the South Eastern Desert. It exhibits a sinistral sense of shear manifested by various kinematic indicators such as biotite fishes, asymmetric porphyroclasts, S-C struc-

Structural Setting and Analysis
The Wadi Hodein shear belt ( Figure 9) is a high angle NW-oriented transcurrent shear zone in the South Eastern Desert. It exhibits a sinistral sense of shear manifested by various kinematic indicators such as biotite fishes, asymmetric porphyroclasts, S-C structures and deflected markers [34,40,74]. Evidence for a dextral sense of shear superimposed on the main sinistral sense of shear was also reported within the W. Kharit-W. Hodein shear zone [74]. Abdel-Karim et al. [75] described high-grade regional metamorphism contemporaneous with the regional foliation and intrafolial folds. A later phase was mainly a dynamic metamorphism associated with the emplacement of ophiolites onto the gneissic rocks. The Wadi Rahaba-Gabal Abu Dahr and Gabal Arais ophiolites include imbricate thrust-bounded sheets and slices of serpentinized ultramafic rocks, amphibolite and metabasalt, embedded in a highly tectonized matrix of talc carbonate, schists and metasiltstones. According to several authors, the Wadi Hodein area has evolved throughout four phases of deformation, D1-D4 [11,[37][38][39][40]61,76]. The D1 deformation phase was an early NNE-SSW crustal shortening related to arc-arc collision and is manifested by S1 axial planar foliations, WNW-trending tight intrafolial and overturned folds (F1) and SSEverging thrusts [37,61]. D2 is expressed in NNW-SSE folds and NNW-SSE crenulation cleavage (S2). F2 folds have en-echelon geometry and verge towards the WSW [37][38][39]61]. The NNW-SSE crenulations and kink folds (F2) might have been developed through oblique non-coaxial deformation of cleaved rocks. D3 was a sinistral transpression along the NNW-SSE ductile shear zones due to NW-ward nappe stacking and thrusting [37]. D4 was a brittle deformation event that led to formation of the ENE-WSW dextral strike-slip faults that deformed the preexisting rocks and dislocated the earlier structures [37,61].

Wadi Khashab Shear Zone
The NNW-trending Wadi Khashab shear zone ( Figure 9) and related splays cut across the southwest part of the area, and it dips moderately or steeply to the east. This shear zone is cut by WNW-ESE and ENE-WSW strike-slip faults. The Wadi Khashab shear zone hosts several occurrences of gold-bearing quartz veins that are commonly associated with sericitized and silicified rocks. Of these occurrences, the El-Beida and Wadi Khashab occurrences are heavily worked out by artisanal miners at present. The Wadi Urga occurrence is confined to highly sheared metavolcanic rocks with signs of hydrothermal alteration along their contacts with a large granodiorite-tonalite intrusion. The~E-W Urga fault is a major, fault that dissects the Wadi Khashab shear zone and apparently displaces its tracẽ 1-2 km right-laterally. Similar parallel structures are abundant north and south of this major fault. A number of prospective targets, including Wadi Beitan, Wadi Urga Rayan, Um Teneidab and north Wadi Hutib, have been identified along similar,~E-W trending faults, or where the main NNW structural trend is cut by WNW-ESE strike slip faults. To the north, it seems that this structure may host other unexplored occurrences, based on the widespread carbonate alteration in the footwalls of the ophiolitic blocks.

Wadi Rahaba Shear Zone
The Wadi Rahaba shear zone ( Figure 9) represents an intense NNW-trending, steeply dipping, brittle-ductile shear zone cutting a generally carbonated mafic-ultramafic rocks and volcaniclastic successions in the central and northern parts of the study area. The shear zone also extends to the south, where it hosts the Anbat mine and other scattered surface mining sites. Towards the east of the Wadi Rahaba shear zone, a major thrust fault boundary marks the base of a hanging wall block of metagabbro and meta-granodiorite; and varies significantly from NW-trending and moderately steeply east-dipping in the central-east to shallow south-dipping and easterly trending in its northeastern part. A number of parallel, NW-trending, shallow to moderately easterly dipping thrust faults (or splays) occur immediately west of the Wadi Rahaba shear zone. These thrust faults mark the western boundary of the gneissic domains and are spatially associated with a number of prospective targets. Numerous significant, E-W to ENE trending late cross-faults appear to offset the major thrust-faults within this area. D4 was a brittle deformation event that led to formation of the ENE-WSW dextral strikeslip faults that deformed the preexisting rocks and dislocated the earlier structures [37,61].

Wadi Khashab Shear Zone
The NNW-trending Wadi Khashab shear zone ( Figure 9) and related splays cut across the southwest part of the area, and it dips moderately or steeply to the east. This shear zone is cut by WNW-ESE and ENE-WSW strike-slip faults. The Wadi Khashab shear zone hosts several occurrences of gold-bearing quartz veins that are commonly associated with sericitized and silicified rocks. Of these occurrences, the El-Beida and Wadi Khashab occurrences are heavily worked out by artisanal miners at present. The Wadi Urga occurrence is confined to highly sheared metavolcanic rocks with signs of hydrothermal alteration along their contacts with a large granodiorite-tonalite intrusion. The

Wadi Beitan Fault Zone
The Wadi Beitan fault zone ( Figure 9) is a major, steeply dipping, WNW-trending silicified breccia zone with a strike length of some +25 km and is traced across the Um Eleiga mine area and eastwards. This structure cuts rocks of variable types in the area, including the Um Eleiga gabbro, Abu Dahr ophiolites and Rahaba gabbro-diorite complex. This zone converges with a NNW-SSE shear zone to the east and delineates the eastern structure contact of the island arc and ophiolitic rocks. Hydrothermal breccia and silica alteration of the host rocks are common along this structure, and abundant zones of sulfide mineralization with local malachite occurrences are observed. In the vicinity of the Um Eleiga gold mine, albitization of the gabbroic host rocks are pervasive. Iron oxides and copper hydro-carbonate are observed in some localities within the Um Eleiga mine area. Noteworthy, the Wadi Beitan fault zone could have been genetically linked with the extensive dike swarms spread over the large are between G. Handusa and G. Nukeiba northwest of the Abu Dahr ophiolitic massif. These dike swarms overprint the earlier NNW-SSE structural trend. Numerous basic dikes along the WNW-ESE direction cut the Um Eleiga gabbro-diorite intrusion and extend eastward to the Abu Dahr ophiolites.

Deformation Events
The Wadi Hodein shear belt (Figure 9) is a NNW-striking belt (~105 km-long), deformed greenschist metamorphosed ophiolites and island arc rocks. It was developed via four phases of deformation; (i) D1, NNE-SSW crustal shortening, (ii) D2: NE-SW oblique convergence and transpression, (iii) D3, E-W compressional regime and (iv) D4: extension tectonics. Lithological layering represents the primary bedding surface (S0) which is still locally recognizable in the ophiolitic metagabbros and bedded metamudstone as cm-scale alternating lighter and darker layers.

D1: NNE-SSW Crustal Shortening
N-S shortening and terrane accretion occurred as a result of a collision between the Gabgaba and Gerf terranes over a~N-dipping subduction zone [77]. D1 structures are mainly preserved in the Wadi Khuda gneisses and in the Wadi Khashab-Gabal Sirsir ophiolitic mélange. Wadi Khuda gneisses represent a small structural window later cut by tonalite-granodiorite and biotite granite intrusions. In the Wadi Khuda gneisses, quartzfeldspar-micas gneissosity (S1) strikes ENE-WSW, and dips toward the NNW and SSE (Figure 10a). The axial plane gneissosity is only evident in the hinges of F1 folds which represented by metrewide, high-amplitude isoclinal, near recumbent to shallowly dipping folds. Boudins (B1) are developed on the quartzo-feldspathic layers in banded amphibolites. They have a symmetric shape and their long axes are parallel to the gneissosity (S1), suggesting a coaxial strain. The largest F1 fold is an asymmetric, doubly plunging anticline extending for up to 8 km in the Wadi Khuda gneisses (Figure 9). The axis of F1 major anticline is displaced by a NNW-striking sinistral strike-slip fault and the southern limb is superimposed by a NNW-tending anticline (F2). The early thrusts are SSE verging and are folded during the D2 event. These thrusts are preserved in the Arais ophiolitic mélange, marking the contacts between Beitan gneisses belt and the ophiolite blocks, as well as the contact between serpentinites and volcaniclastic metasediments ( Figure 9). These thrusts converge downwards onto a basal S-dipping decollement.

D2: NE-SW Oblique Convergence and Transpression
The D2 deformation phase is manifested by regional foliation (S2), mylonitic foliation (S2m), major and minor folds (F2), stretching and intersection lineations (L2) and major thrusts (T2). The structures include small-scale, rootless intrafolial folds (Figure 10b), shearband boudins (Figure 10b), recumbent (Figure 10c,d), tight isoclinal and sheath folds along the first-generation ductile shear zone. These structures are largely preserved as decimetre-scale folds in a transposing and intrafolial foliation (S2) in the Beitan gneisses belt, schistose metavolcanics and volcaniclastic metasediments. The S2 in foliated metavolcanics and volcaniclastic metasediments is axial plane schistosity for the F2 folds and is defined by a strong preferred alignment of actinolite, chlorite and biotite. In the bedding-parallel quartz veins, asymmetrical "shearband-type" boudins (B2; Figure 10e) and asymmetrical "domino-boudinaged" quartz veins are common. The shearband boudins (B2) were only observed in shear zone close to the gold occurrences. Ptygmatic folding is widespread in the Beitan gneisses (Figure 10f). S2 planes have moderate to high dips predominantly toward the SW and NE. Around the synorogenic gneissose granites, ophioltic metagabbros and serpentinite mountainous blocks, the S2 foliations show important changes along strike or define different geometric patterns. In the central segment of Hodein shear belt, S1 traces strike NW-SE to NNW-SSE to NNE-SSW, whereas they curve back into a NW-SE orientation to the southwest of Gabal Dahr serpentinite's peridotites.
anticline extending for up to 8 km in the Wadi Khuda gneisses (Figure 9). The axis of F1 major anticline is displaced by a NNW-striking sinistral strike-slip fault and the southern limb is superimposed by a NNW-tending anticline (F2). The early thrusts are SSE verging and are folded during the D2 event. These thrusts are preserved in the Arais ophiolitic mélange, marking the contacts between Beitan gneisses belt and the ophiolite blocks, as well as the contact between serpentinites and volcaniclastic metasediments (Figure 9). These thrusts converge downwards onto a basal S-dipping decollement.

D2: NE-SW Oblique Convergence and Transpression
The D2 deformation phase is manifested by regional foliation (S2), mylonitic foliation (S2m), major and minor folds (F2), stretching and intersection lineations (L2) and major thrusts (T2). The structures include small-scale, rootless intrafolial folds (Figure 10b  Other F2 axes may have resulted from later reworking by F3. Axial planar to these folds is a spaced crenulation cleavage (S2), which is commonly well developed. F2 folds range from centimeter-to kilometer-scale and are common in Beitan gneisses belt and the ophiloitic mélanges (Figure 11a The ophiolitic mélanges in the Wadi Hodein-Wadi Beitan shear belt are thrusted over each other and over the Beitan gneisses belt (infracrustal rocks) along three major imbricate thrusts (T2). Three imbricate thrusts are named in this study: the Sirsir-Rahaba thrust, the Arais imbricate thrust and the Hodein-Beitan imbricate thrusts (Figure 9). The geometry of Wadi Khashab-Gabal Sirsir ophiolitic mélange is controlled by the oppositely dipping Sirsir-Rahaba and the southern extension of Wadi Hodein-Wadi Beitan sinistral reverse shear zone, respectively (Figure 12a-c). Gently to moderately ENE-and WSW-dipping foliation (S2) forms the dominant fabric that intensifies into belts of mylonites associated with D2 ductile thrusts. The Sirsir-Rahaba thrust dips ENE and separates ophiolitic blocks of Gabal Sirsir and the volcaniclastic metasediments of Gabal Anbat (Figures 9 and  12a-c). The Wadi Hodein-Wadi Beitan imbricate thrust forms an imbricated fan ( Figure  12a-c). Axial planar to tight-to-isoclinal, reclined F2 sheath folds have axes concentrated around the L2 mineral lineation. To the north of Gebel of Abu Dahr the thrusts mark the tectonic contacts between the ophiolitic slabs (serpentinites and talc carbonates) and the overlying metavolcanics and volcaniclastic metasediments. Additionally, Arais imbricate thrusts control the contact between Arais opoilitic mélange and the Beitan gneisses belt ( Figure 9) and also mark the tectonic contacts between the ophiolitic slices. These NE- The ophiolitic mélanges in the Wadi Hodein-Wadi Beitan shear belt are thrusted over each other and over the Beitan gneisses belt (infracrustal rocks) along three major imbricate thrusts (T2). Three imbricate thrusts are named in this study: the Sirsir-Rahaba thrust, the Arais imbricate thrust and the Hodein-Beitan imbricate thrusts (Figure 9). The geometry of Wadi Khashab-Gabal Sirsir ophiolitic mélange is controlled by the oppositely dipping Sirsir-Rahaba and the southern extension of Wadi Hodein-Wadi Beitan sinistral reverse shear zone, respectively (Figure 12a-c). Gently to moderately ENE-and WSW-dipping foliation (S2) forms the dominant fabric that intensifies into belts of mylonites associated with D2 ductile thrusts. The Sirsir-Rahaba thrust dips ENE and separates ophiolitic blocks of Gabal Sirsir and the volcaniclastic metasediments of Gabal Anbat (Figures 9 and 12a-c). The Wadi Hodein-Wadi Beitan imbricate thrust forms an imbricated fan (Figure 12a-c). Axial planar to tight-to-isoclinal, reclined F2 sheath folds have axes concentrated around the L2 mineral lineation. To the north of Gebel of Abu Dahr the thrusts mark the tectonic contacts between the ophiolitic slabs (serpentinites and talc carbonates) and the overlying metavolcanics and volcaniclastic metasediments. Additionally, Arais imbricate thrusts control the contact between Arais opoilitic mélange and the Beitan gneisses belt (Figure 9) and also mark the tectonic contacts between the ophiolitic slices. These NE-dipping thrusts constitute a backthrust to the WSW-dipping Hodein-Beitan imbricate thrusts.  Most of the small-scale, tight to isoclinal, asymmetric F2 folds in the Wadi Hodein-Wadi Beitan imbricate thrusts show top-to-the-ENE thrusting. Many shear sense indicators, including mineral lineations, asymmetric shear folds, asymmetric porphyroblasts, shearband boudins, S-C fabrics and shear bands indicate that the ophiolite blocks along Sirsir-Rahaba thrust and Arais thrusts are moved top to WSW and those along Hodein-Beitan thrusts are moved up to ENE (Figure 13a-c). Each thrust zone comprises several thrust segments showing anastomosing morphologies and form WSW-dipping tectonic duplexes. S2 gneissosity in the Wadi Beitan gneisses occur as asymmetrical folds around these boudins. Migmatitization is common in the shadows of the rotated boudins. The boudins do not display any foliations but are surrounded by the S2 gneissosity, suggesting that they formed during the D2 deformation. Shearband boudins show a top-to the-WSW movement. Most of the small-scale, tight to isoclinal, asymmetric F2 folds in the Wadi Hodein-Wadi Beitan imbricate thrusts show top-to-the-ENE thrusting. Many shear sense indicators, including mineral lineations, asymmetric shear folds, asymmetric porphyroblasts, shearband boudins, S-C fabrics and shear bands indicate that the ophiolite blocks along Sirsir-Rahaba thrust and Arais thrusts are moved top to WSW and those along Hodein-Beitan thrusts are moved up to ENE (Figure 13a-c). Each thrust zone comprises several thrust segments showing anastomosing morphologies and form WSW-dipping tectonic duplexes. S2 gneissosity in the Wadi Beitan gneisses occur as asymmetrical folds around these boudins. Migmatitization is common in the shadows of the rotated boudins. The boudins do not display any foliations but are surrounded by the S2 gneissosity, suggesting that they formed during the D2 deformation. Shearband boudins show a top-to the-WSW movement.

D3: E-W Shortening and Sinistral Shearing
D3 structures are related to transposition and reorientation of the D2 structures. S3 schistosity overprinted the pre-existing S2 schistosity. S3 foliation strikes N20-30° W and dips 60-70° to the NE or SW (Figure 12d). Mylonitic foliation in C3 shear zones is superimposed on the S2 foliation (Figure 13d). The spacing of the strike slip shear zones varies from a few centimeters to some meters. Kilometer scale NNW-oriented strike-slip shear zones are commonly developed in the Wadi Khashab-Gabal Sirsir ophiolitic mélange, especially along Wadi Khashab, along the main planes of Hodein-Beitan shear zone and Sirsir-Rahaba shear zone (Figure 9). The ophiolitic and metavolcanic rocks on both sides of the shear zone show S-shaped fault-drag folds (F3), consistent with the sinistral movement. Quartz veins showing signs of ductile shearing are sulfide-bearing and are associated with green malachite alteration in places.
The strike-slip shear zones in Hodein shear belt separate the Wadi Beitan gneisses from the overlying ophiolitic mélange (Figure 9). The mylonitized rocks are characterized by stretched quartz ribbons with an anastomosing pattern of intervening variably deformed lenses. The angle between S3 and C3 is often clearly visible and can reach 5°, indicating strong shear values. Asymmetric boudins are developed in the Wadi Beitan The strike-slip shear zones in Hodein shear belt separate the Wadi Beitan gneisses from the overlying ophiolitic mélange (Figure 9). The mylonitized rocks are characterized by stretched quartz ribbons with an anastomosing pattern of intervening variably deformed lenses. The angle between S3 and C3 is often clearly visible and can reach 5 • , indicating strong shear values. Asymmetric boudins are developed in the Wadi Beitan banded gneisses and amphibolite. Mylonitic rocks exhibit subvertical to steeply dipping foliations (N20 • W and N40 • W) and dip steeply to the SW-dipping high-angle S3 foliations in the high strain zones along the main planes of the Beitan imbricate thrusts and Sirsir-Rahaba thrust. The S3 foliation is defined by kink bands or spaced cleavage that are parallel to the F3 axial planes. S3 kink bands are sharp deflections of S2. The high-strain zones with subvertical foliation planes are subparallel and are oriented NNW-SSW except in the central part of the Wadi Hodein-Wadi Beitan shear belt where they are deflected into NNE-SSW. In these high-strain zones, isoclinal and intrafolial folds with shallow plunging fold hinges and a subvertical axial planar foliation are abundant. Shallow NNW-or SSEplunging mineral and stretching lineations are present on the steep foliation planes in the NNW-SSE striking high-strain zones.
Two distinct types of mineral lineations L3 on the S3 foliations have been recognized by the measurement of the minerals that grew during D3. The first type, here defined as L3a, is a SW and SE-plunging mineral stretching lineation along the ENE-verging thrust (T3). L3a is marked by platy quartz, amphibole fish and feldspar rods. The second type, here referred to as L3b, is a b-type lineation orientated subparallel to the F3 fold axes. The L3b is a crenulation lineation with a subhorizontal plunge to the NNW or SSE (Figure 13e). The D3 deformation produced regional-scale NNW-trending folds (F3) defined by open folds with nearly vertical axial planes. Superposition of F3 on F2 folds resulted in extensive NNW-trending tight and reclined folds and NE-trending open folds. The most common macroscopic shear sense indicators include asymmetric tight and reclined folds (Figure 13f), sigmoidal porphyroclasts (Figure 13g,h), S-C fabrics, folded and en-echelon quartz lenses (Figure 13i), shearband quartz boduins and pop up structures. The lineation indicates that the tectonic transport took place while far-field compression progressed and switched to oblique thrusting.

D4: Extension and Terrane Exhumation
D4 was a prolonged phase of brittle deformation associated with intrusion and emplacement of late-orogenic granitoid intrusions and terrane exhumation. During D4, the Wadi Hodein-Wadi Beitan shear belt was dissected by a number of dextral strike-slip faults, deformed the post-orogenic granites and the D1-D3 structures. These ENE-WSW right-lateral strike-slip faults occur in the central part of the Wadi Hodein-Wadi Beitan shear belt and along Wadi Khuda (Figures 9 and 12d-f). Another major sinistral strike-slip fault striking NNE-SSW controls the course of Wadi Urga El-Rayan. NW-SE oriented gold-bearing quartz veins were originated during D2 and were subsequently deformed by D3-D4 events.

Gold Occurrences in the Wadi Hodein Shear Belt
A number of gold occurrences in the Wadi Hodein shear belt are mainly associated with the high strain zones. The latter are confined to fault-bounded imbricated thrust sheets of allochthonous ophiolitic blocks. The Roman and Ptolemaic gold occurrences are confined to deformed granitoid and metagabbro-diorite complexes in the northern and central parts of the Hodein shear belt. In the following sections, we describe the internal structures of the mineralized quartz veins and the structural controls of several occurrences, namely: Wadi Khashab, the El-Anbat mine, the Um Teneidab mine, the Urga El-Ryan occurrence, the Hutit mine and the Um Eleiga deposit.

The Wadi Khashab Occurrence
In the Wadi Khashab area (Figure 14a), the ophiolite slices are overthrust on the arc's metavolcanic rocks and volcaniclastic metasediments. The serpentinites of Gabal El Beida form a NW-elongated large belt and slices of different sizes tectonically incorporated in the volcaniclastic metasediments. The serpentinites form dissected elongated massive masses with sheared and foliated peripheries. These serpentinite bodies are altered to talc-carbonate rocks along the shear zones. Blocks and masses of ophiolitic metagabbros usually occur within the volcaniclastic metasediments. They locally exhibit weakly developed foliation and compositional layering. Ultramafic rocks are commonly associated with pillow and amygdaloid metabasalts. The metavolcanic rocks are represented by strongly foliated metaandesite, metabasaltic andesite and intermediate to acidic metatuffs [39].
Gold-sulfide mineralization is confined to discrete shear zones of highly silicified, ferruginated metavolcanics and volcaniclastic metasediments commonly associated with sulfide-bearing granophyric dikes and quartz± carbonate veins [39]. The Wadi Khashab zone is a major sinistral strike-slip shear zone, striking NNW-SSE and dipping steeply (65-75 • ) towards the ENE. Rocks along this shear zone are carbonatized and consist mainly of highly sheared metavolcanics, serpentinites and volcaniclastic metasedimentary rocks. The wall rock alteration is present as strongly foliated flakes in quartz veins and is dominated by the silicification and chloritization associated with carbonate [39,78]. The silicified shear zone (SiO 2 > 80%) of the pillow metabasalts is relatively enriched in Au (25-45 ppm) whereas the carbonatized shear zone (SiO 2 < 50%) within the obducted area of the serpentinite rocks has a very low Au content (0.24-0.45 ppm) [78]. The samples analyzed by Zoheir [39] yielded 2.1-6.6 ppm Au, and the sulfidized granophyric dikes are also gold-bearing (0.2-2.1 ppm Au).
In the El-Anbat mine area (23 • 17 01 N, 35 • 09 21 E), the alteration zones occur as sheets or lenses along the NW-and NNW-striking thrusts, between the carbonatized serpentinite and foliated intermediate metavolcanics (Figure 14b). S2 is a regional foliation in El-Anbat area striking NNW-SSW and dipping moderately or steeply to the ENE and less commonly to the WSW. The F1 folds are still preserved in ophiolitic metagabbros and mafic metavolocanics. The F2 folds are major NNW-trending asymmetric left-stepping folds in the volcaniclastic metasediments. Different styles of minor folds are recognizedopen concentric folds, closed tight folds, symmetric and asymmetric folds, isoclinal folds and chevron folds. During D3, Conjugate sinistral and dextral strike-slip faults and shear zones are apparently developed as a component of a major transpression system.
The mineralized shear zones and the bleached sheared host rocks are confined to the steeply dipping thrusts between carbonatized serpentinites and the underlying volcaniclastic metasediments. These zones are characterized by arrays of recrystallized quartz pods, fracture-filling ankerite, sulfidized granitic offshoots, boudinage of the granitic bodies, silicified listvenite masses and intermingling carbonate sericite and malachite veinlets [38], suggesting significant gold grades in the hydrothermally altered, sheared host rocks. The analyses indicate that gold is related to mineralized quartz veins (with up to 7.5 g/t Au) and to pervasively altered wall rocks (with up to 5 g/t Au) [38]. mine down to a level of 15 m below the surface. These shear zones (up to 40 m in length and 5-30 cm in thickness) are strongly mylonitized and marked by recrystallized quartz lenses and asymmetric boudinaged quartz veins [12]. The narrow zones of sheared wall rocks are characterized variable degrees of hydrothermal alteration. The gold grade is 1 to 7g/t Au in the quartz veins from the Urga El-Ryan occurrence [84].

The Um Teneidab Deposit
The Um Teneidab mine (23 • 17 01 N, 35 • 09 21 E) is located~13 km to the southwest of the Hutit mine and to the west of Gabal Um Teneidab peak. The Um Teneidab area is underlain by gneisses and migmatites, serpentinite, volcaniclastic metasediments, and deformed mafic metavolcanic rocks. These rocks are cut by metagabbro-diorite complex (Figure 14c), as well as unmappable masses and offshoots of post-orogenic leucogranites [12,81]. Serpentinite constitutes the main mass of Gabal Um Teneidab and forms blocks and slices of variable sizes incorporated in the volcaniclastic metasediments. The mélange blocks are randomly distributed, varying from cobbles to mountain size and usually oriented in the NWSE or NNW-SSE directions. The gabbroic rocks in the Um Teneidab mine area were mapped as island arc, synorogenic metagabbro-diorite complex [12,82] and as post-orogenic younger gabbros [81]. In the present work, we agree with Hassan and El-Manakhly [82] and Zoheir et al. [12] and consider the gabbroic rocks in the Um Teneidab mine area as a variably foliated, metagabbro-diorite complex.
The Um Teneidab metagabbro-diorite complex is dissected by a number of ENE-WSW and N-S dextral strike-slip faults and WNW-ESE major fractures with no obvious lateral displacement. The zones in which the leucogranites send offshoots into the gabbroic rocks (Figure 14c) are characterized by intense shearing and alteration as well as the presence of abundant quartz veins and felsic dikes. Shearing in the area is predominantly brittle, with little ductile deformation being experienced, with alteration zones and quartz veins that are generally controlled by NW-SE shear/fault sets. The subvertical NW-and NNWstriking milky quartz veins extend up to 250 m with 8 cm to 45 cm thickness and lensoid morphology. These veins are in part or completely recrystallized. The main lode is a zone of stockwork of veinlets (70 cm wide) bordered by hydrothermally altered wall rocks forming together a~2-m-wide mineralization zone [12]. The gold in these veinlets and the hydrothermally altered wall rocks occurs as fillings in the microfractures and is scattered as sliver in altered pyrite and galena. The gold content ranges from 1 to 30 g/t in quartz veins, whereas the altered wall rocks locally contain~8 g/t Au [12].

The Urga El-Ryan Occurrence
The Urga Ryan gold deposit (35 • 5 0 E and 23 • 21 23 N) is located at~3 km to the south of the intersection between Wadi Urga Ryan and Wadi Hutib. The area around Wadi Urga El-Ryan (Figure 14d) is underlain by gneisses, volcaniclastic metasediments and mafic to intermediate metavolcanic rocks. The volcaniclastic metasediments include wide varieties of metagreywacke, metamudstone, slate, pelitic schists and metaconglomerate in decreasing order of abundance. They are strongly NW-SE foliated and less commonly NNE-SSW cleaved [48].
The metavolcanics are composed of repeatedly alternating and interfingered successions of basaltic flows together with layered andesite and subordinate dacitic flows and their corresponding metatuffs [48] as well as chlorite schists. The metavolcanic rocks are variably deformed by km-scale shear zones that led to intense shearing in the NNW-SSE direction (S2 and S3), overprinting the WNW-ESE schistosity (S1). The metagabbro is emplaced into both metavolcanics from the northern and western parts and volcaniclastic metasediments from the southern part [48]. The metagabbroic rocks are foliated and sheared, particularly along the shear zones. The syn-orogenic granites are slightly foliated and form an elongated intrusion (10 m × 2.5 m width) trending roughly N-S, which forms discontinuous hills. It is mainly tonalite in composition and is bounded from the west by the Wadi Beitan gneisses, and by metagabbro from the north. Volcaniclastic metasediments and metavolcanics border these rocks from the southeastern side. Xenoliths and roof pendants from metagabbro are enclosed within the syn-orogenic tonalite.
In the Urga El-Ryan occurrence (Figure 14d), the old mining houses spread along the main Wadi Urga Ryan and their tributaries, possibly reflecting significant mine activities. The mineralized regions are limited only to NNW-SSE striking and westward dipping zones of intense shearing in metavolcanic rocks. The quartz vein strikes N-S (355 • ) concordantly with the shear zone system and dips at 80 • E [83]. It was exploited in a deep mine down to a level of 15 m below the surface. These shear zones (up to 40 m in length and 5-30 cm in thickness) are strongly mylonitized and marked by recrystallized quartz lenses and asymmetric boudinaged quartz veins [12]. The narrow zones of sheared wall rocks are characterized variable degrees of hydrothermal alteration. The gold grade is 1 to 7g/t Au in the quartz veins from the Urga El-Ryan occurrence [84].

Hutit Deposit
The Hutit gold deposit occurs in the Wadi Huzama, which is a small tributary of the larger Wadi Rahaba (Figure 15a). The Hutit mine area is underlain by serpentinite, mafic metavolcanics, volcaniclastic metasediments, ophiolitic metagabbros, syn-orogenic metagabbros and post-orogenic granites (Figure 15a,b) [12,48,82,84,85]. Structurally, the Hutit deposit is controlled by the Rahaba-Sirsir sinistral shear zone and their kilometerscale imbricate thrust system (Figure 15b). The SW-directed thrusts bound the ophioltic blocks and dip moderately or steeply to the NE. Imbricate slices and sheets of serpentinite are tectonically overlying metavolcanics and volcaniclastic metasediments. The oblique thrust planes accommodate sinistral displacement and asymmetrical anastomosing and kinked fabrics. Asymmetric quartz lenses and drag folds indicate a sinistral sense of shear, whereas subvertical slickensides along the quartz vein walls document the reverse slip of the hanging wall block. Joints and fractures are abundant in the rock units mainly with NE-SW and NW-SE directions. Mafic to felsic dikes (basalt, basaltic andesite, andesite, rhyodacite and dacite) as well as granodiorite dikelike bodies strike NW-SE and less commonly NE-SW.

Hutit Deposit
The Hutit gold deposit occurs in the Wadi Huzama, which is a small tributary of the larger Wadi Rahaba (Figure 15a). The Hutit mine area is underlain by serpentinite, mafic metavolcanics, volcaniclastic metasediments, ophiolitic metagabbros, syn-orogenic metagabbros and post-orogenic granites (Figure 15a,b) [12,48,82,84,85]. Structurally, the Hutit deposit is controlled by the Rahaba-Sirsir sinistral shear zone and their kilometer-scale imbricate thrust system (Figure 15b). The SW-directed thrusts bound the ophioltic blocks and dip moderately or steeply to the NE. Imbricate slices and sheets of serpentinite are tectonically overlying metavolcanics and volcaniclastic metasediments. The oblique thrust planes accommodate sinistral displacement and asymmetrical anastomosing and kinked fabrics. Asymmetric quartz lenses and drag folds indicate a sinistral sense of shear, whereas subvertical slickensides along the quartz vein walls document the reverse slip of the hanging wall block. Joints and fractures are abundant in the rock units mainly with NE-SW and NW-SE directions. Mafic to felsic dikes (basalt, basaltic andesite, andesite, rhyodacite and dacite) as well as granodiorite dikelike bodies strike NW-SE and less commonly NE-SW.    (Figure 15b) were mapped at the 1:1000 scale by Hassan and El-Manakhly [82]. The ruins of the separation and grinding plants are still observed in the northern old mine, whereas the loading station, leaching basins and crusher stages were observed in the southern mine. In both mines, a main entrance through a horizontal~E-W adit leads to the veins at a distance of 20 m or 35 m [12]. The mineralized quartz veins occur along the thrust contact between the serpentinite masses and the mafic metavolcanics and volcaniclastic metasediments. Two types of anastomosing and undulating gold-bearing quartz veins were reported in the mine area, including bluish-gray (common in the northern mine) and milky quartz veins (dominant in the southern mine). Concerning the ore grade, Gabra [84] reported 1-36 g/t in samples from quartz veins the southern mine and 1-40 g/t in quartz veins from the northern mine. Takla et al. [85] reported an average of 20 g/t in samples from the two different types of quartz veins and 8 g/t Au on average in altered wall rocks.

The Um Eleiga Deposit
The Um Eleiga mine (24 • 36 44" N, 35 • 03 19" E) occurs~50 km west of the Red Sea coast. Traces of placer gold workings during Roman-Byzantine and early Islamic times include hundreds of shallow pits, dumps, ancient mining camps, stone anvils, hammers, and grinding mills [12,41,83]. The geology of Um Eleiga area is dominated by a 7 km wide intrusive complex with a gabbroid core cutting the serpentinite-chromitite of Gabal Abu Dahr and the northern part of Wadi Rahaba-Gabal Abu Dahr ophiolitc mélange. The latter composed masses and blocks of serpentinized ultramafics, pillow metabasalts and mafic metavolcanics incorporated within highly tectonised matrix of pelitic and carbonaceous metasediments. The ophiolite fragments occur as NE-dipping slabs and sheets thrusted southwestward forming imbricate thrust system. Contacts of the serpentinized ultamfics with Um Eleiga metagabbros-diorite complex are sharp and irregular. In the upper part of the Abu Dahr massif, serpentinite contains large masses of serpentinized dunite, sills or layers of wehrlite, pyroxenite and gabbro [42,86].
The mafic metavolcanic rocks comprise blocks and thrusted slices of metabasalt and basaltic andesite [42,86]. These rocks structurally overlay the Um Eleiga intrusive metagabbros-diorite complex and Abu Dahr serpentinized ultramafics (Figure 15c). Pillow metabasalt occurs in the core of a major fault-propagation fold (Figure 15c) within the highly tectonized matrix with schists and sheared serpentinite. Pillow metabasalt forms ellipsoidal, globular and tabular shaped bodies, ranging between 10 cm to 1 m in diameter. The Um Eleiga intrusive metagabbro-diorite complex is approximately 7 km long and 4 km across, and encompasses a compositional continuum from gabbro to granodiorite through diorite and subordinate tonalite (Figure 15d). This complex exhibits a circular zonation, in which gabbro occupies the core and subordinate granodiorite forms the margin [41]. The Um Eleiga metagabbro and metabasalt have a tholeiitic to calc-alkaline affinity and are interpreted to have been formed in a forearc setting [42].
In the Um Eleiga mine, the mineralized quartz veins trend mainly NE-SW or ENE-WSW and cut the gabbroic rocks in the central part of the complex and extend beyond the gabbro-diorite boundary [12]. The fault/joint intersections are the main structural control of intensely hydrothermal alteration zones and high gold contents in the central part of the Um Eleiga complex [41]. Late barren quartz veins are confined to intersections of fault and tension gashes in the brecciated gabbro core and trend mainly N-S, NW-SE and E-W (Figure 15d). Sulfide-bearing quartz veins (5-40 cm thick) cut the metagabbro-diorite complex and are partially stopped out [12]. Analyses of quartz dumps in several pits gave gold values up to 28 g/t [82]. According to Takla et al. [85], most of the ancient workings in the mine area were confined to the old terraces (lithified wadi alluvium) within and adjacent to the intensely kaolinitised gabbro. Atomic absorption analyses of heavy concentrates of these old terraces indicated that gold contents were variable from barren to 10 g/t [85].

Remote Sensing Targeting of New Gold Occurrences
The Nubian Shield is a typical example of well-exposed crystalline basement rocks, and mafic-carbonate-hydrous mineralogical indices were used to extract the representative pixels in the satellite images. Zoheir et al. [12] used Landsat-8 OLI and ASTER band ratio images and successfully benefited from the effective absorption features of the mafic rock-forming minerals and their metasomatic products. Therefore, Landsat-8 OLI and ASTER data with comprehensive fieldwork were used for mapping geological structures, lithological units and alteration zones associated with orogenic gold mineralization in the study area. The image processing techniques of these remotely sensed data maximize the relatively small differences in the spectral responses of rock forming minerals within a specific wavelength range [87]. The applied approach includes false band combination (FCC), band rationing, principal component and mineralogical indices to the Landsat-8 OLI and ASTER data. It is evident that the band math of Landsat-8 (6/7) distinguished clay minerals, serpentine, and many alteration zones, whereas ASTER (7 + 9/8) demarcated carbonatized rocks, metavolcanics and metasedimentary rocks.
The mineralogical indices of Ninomiya [62] (OHI, KLI, CLI and ALI) were used to characterize the alteration zone in the study area using ASTER data. The resultant images highlighted effectively the zone of the hydrothermal alterations. In addition, they stressed the fact that the alteration minerals zones are found in both ductile and brittle deformation zones where the integration between the structural segments (density and direction) associated with the presence of the altered minerals may help to localize and predict new occurrences of gold mineralization. A false-color composite (FCC) ratio image (OHI, KLI, CLI) was generated. This RGB FCC image shows the OH-bearing rocks (metaultramafic-mafic ophiolites and island arc meta-basic rocks) and the clay-minerals (see Figure 6). The evaluation of the remote sensing alteration mapping results with field data using an error matrix approach and Kappa Coefficient shows a very good match, which indicates the overall accuracy of 88.33% and the Kappa Coefficient of 0.77 for Landsat-8 data and the overall accuracy of 74.10% and the Kappa Coefficient of 0.65 for ASTER data, respectively (Table 4A,B). It is noticed that the old working gold occurrences are associated with the high deformation and fracturing zones bounding the sheared ophiolitic belt and faulted contacts such as El-Beida, Urga Ryan, Hutit and Um Eleiga gold mines. Multisource datasets were selected to control factor layers indicating the distribution of gold mineralization potentiality zones such lithology, alteration zones, proximity to existing gold mines, proximity to favorable contacts, proximity to major structures (especially thrusts), faults density maps and lineaments density maps. The seven evidence layers were derived, analyzed, and integrated in a GIS platform to develop the gold mineralization model (see Figure 8).

Transpressional Tectonics in the Evolution of the Wadi Hodein Shear Belt
The Wadi Hodein shear belt has evolved throughout a multistage deformation history (D1-D4), in which D2 and D3 were the most dominant. NNE-SSW crustal shortening (D1) led to development of early ENE-WSW striking S1 gneissosity and F1 doubly plunging anticline in the Wadi Khuda gneisses. D1 structures in the Wadi Khashab, Wadi El Beida and Gabal Anbat and Gabal Arais areas exhibit WNW-striking S1 schistosity, F1 minor tight, recumbent and overturned folds and T1 thrusts. NE-SW oblique convergence and transpression (D2) led to the development of NW-SE striking S2 regional foliation, S2m mylonitic foliation, F2 major kilometer scale anticlines and synclines (see  and minor open, asymmetrical to isoclinal folds, NW-plunging stretching and mineral lineations (L2) and WSW-verging major thrusts (T2). The major D2 fabrics are represented by a series of large-scale folds (F2) associated with northwest-vergent foliation (S2), indicating that the F2 folds formed in response to a NW-SE-oriented stress field. The Wadi Hodein-Beitan shear belt is cut by a suite of variably deformed granitic intrusions. Many of these intrusions occur in the vicinity of the major thrusts and are considered to have emplaced during the regional D2 deformation event. E-W shortening (D3) The Wadi Hodein-Beitan shear belt is cut by a suite of variably deformed granitic intrusions. Many of these intrusions occur in the vicinity of the major thrusts and are considered to have emplaced during the regional D2 deformation event. E-W shortening (D3) produced subvertical NNW-Striking foliation (S3), NNW-SSE thrust dominated strike slip shear zones (C3), NNW-and ENE-trending folds (F3). L3a is a SW and SE-plunging stretching lineation along the NNE-SSW T3 minor thrusts. L3b lineation is a crenulation lineation with a subhorizontal plunge. During D3, the major NNW-oriented shear zones rooted into a low-angle, deep crustal detachment, delineating the point at which oblique thrusting is partitioned into separate transpressional strike-slip and reverse-sense shear zones, i.e., Sirsir-Rahaba, Hodien-Beitan and Um Teneidab-Hutit shear zones. The sinistral porphyroclasts and NNW-plunging reclined and vertical folds as well as NE-trending open folds in sheared metavolcanics and volcaniclastic metasediments indicate a transpressive stress regime including both pure and simple compressional components. The latest extensional structures (D4) in Hodein shear belt include major brittle strike slip faults and microfaults and dikes associated with the intrusion of late to post-tectonic granites. Structures assigned to this deformation phase include ENE-WSW dextral and NNE-SSW sinistral strike-slip faults.
The newly produced structures map (see Figure 9) shows that the overall structural setting of the Wadi Hodein shear belt bears a resemblance to an asymmetric flower structures in which the ophiolitic mélange and the ophiolitie blocks are cut and bounded on both sides by reverse-slip faults, and strike-slip shear zones. The eastern (Sirsir-Rahaba and Um Teneidab-Hutit shear zones) and western boundary faults (Hodein-Beitan and Arais shear zones) correspond to a reverse-sinistral thrusts. The transpressional imbricate thrusts forming the Sirsir-Rahaba, Hodein-Beitan and Arais shear zones dip toward ENE, whereas those pertaining to the Um Teneidab-Hutit shear zones dip to WSW (see Figures 9 and 16d). All the shear zones show a moderately sinuous trace in map-view, up to 2 km-thick band of fine-laminated mylonites, developed in the footwall gneisses and volcaniclastic metasediments. The high-angle, NNW-oriented thrust-dominated, strike slip shear zones marked a near-vertical ductile deformation zone rather than an imbricate structure [88][89][90]. The thrusts show a typical sequence of propagation in the footwall of the previous thrust in a piggy-back model. Many kinematic indicators, including asymmetric porphyroblasts such as "σ-type" augens, and pressure shadows, indicate a sinistral sense of shearing. The predominance of reverse slip movement over the along strike shearing suggests a thrust-dominated transpression along the Wadi Hodein shear zone. The NW-SE to NNW-SSE subparallel Sirsir-Rahaba, Hodein-Beitan, Um Teneidab-Hutit and Arais strike-slip shear zones exhibit a consistently sinistral sense of shear.

Structures Controlling Gold Occurrences
The distribution of the gold occurrences is confined to the margins of the Wadi Hodein shear belt (see Figure 9) especially along the thrust-dominated transpressional major zones and their splays formed during D3. Conjugate sinistral (N) NW-(S) SE and dextral (E)NE-(W)SW strike-slip faults and shear zones apparently developed During D3 as a component of a major transpression system. These conjugate shear zones control the occurrence of gold in Hodein shear belt. Twisting and anastomosing S2 and T2 thrusts led to NE-trending F3 open folds and oblique sinistral shear zones on the pre-existing NW-striking thrusts. The S-C fabrics and asymmetric quartz lenses replicate left-lateral shearing, which resulted in widespread NE-SW-striking dilation bents within and around the shear zones. At the Hutit and El-Anbat mine areas, the mineralized quartz veins and the alteration zones occur along the NW-and NNW-striking thrusts between the carbonatized serpentinite and foliated metavolcanics and the metavolcanics and volcaniclastic metasediments. In the Um Eleiga mine, the mineralized quartz veins trend mainly NE-SW or ENE-WSW mainly along the fault/joint intersections. In the Urga El-Ryan occurrence, the mineralized regions are limited only to NNW-SSE striking and westward dipping zones of intense shearing in metavolcanic rocks. In the Um Teneidab deposit alteration zones and quartz veins are generally controlled by NW-SE shear/fault sets. In the Wadi Khashab occurrence, the quartz veins occur as closely spaced swarms running parallel and subparallel to NNWtrending sinistral strike-slip shear zone of Wadi Khashab.

Conclusions
Detailed geological mapping of the Wadi Hodein area was successfully accomplished by employing different processing techniques, i.e., band combinations, band math (BM), Principal Component Analysis (PCA), decorrelation stretch and mineralogical indices, to Landsat-8 OLI, ASTER and ALOS PALSAR data. Field data were integrated with results of the remote sensing studies, which enabled an updated understanding of the structural evolution of the Wadi Hodein-Wadi Beitan shear belt. The structural framework of the Wadi Hodein-Wadi Beitan shear belt was shaped by the superposition of the NW-SE folds (F2) by the NNW-trending, km-scale tight and reclined folds (F3). The overall configuration of the investigated shear belt is identical to that of an asymmetrical flower structure in which the ophiolitic mélange and the ophiolitic blocks are cut and bounded on both sides by reverse slip and strike-slip shear zones. The field studies showed that most of the gold mineralization is mainly associated with quartz veins within extensive sinistral transpressional shear zones assigned to the last ductile deformation event (D3) in the evolution of the entire shear belt. This event was characterized by strong stretching and mylonitic lineation, abundant undulatory shear zones and shear foliation disturbance. The gold-bearing quartz veins are hosted by narrow, steeply SW-and NE-dipping NNW-and NW-trending mylonitic zones that are characterized by a steeply NNW-plunging mineral lineation. This may imply the important role of oblique convergence in the formation of the auriferous quartz veins.
The processing flow chart introduced in this study has identified areas with the potential for unexplored gold occurrences and helped to illustrate the controls of the distribution of gold mineralization sites. The processing results were classified into five categories, from very high to very low potentiality zones. Gold mineralization is here interpreted to show preferential distribution along the Wadi Hodein-Wadi Beitan shear belt. The most important conclusion of this study is that gold mineralization in the Wadi Hodein-Wadi Beitan shear belt was mostly confined to steeply dipping strike-slip shear zones in the marginal parts of the shear belt. The occurrence of gold-bearing quartz veins in central shear zone is not evident so far. A conceptual model for testing new and unexplored targets of gold-quartz veins in the region translates the coincidence of shear zones, alteration and acid dikes to highly potential zones for new and possibly important gold occurrences in the study area.  Figure S3: Flowchart of the remote sensing and GIS methodologies adopted in the present study; Figure S4: Lithological discrimination using (a) ASTER (RGB-431) and (b) grey scale ASTER band ratio (7 + 9/8); Figure S5: False-color composite of principal component analysis (PCA) of (a) Landsat-8 RGB-PC2, PC1, PC4 and (b) ASTER RGB-PC1, PC2, PC3; Figure S6: Model Builder of the geospatial thematic maps used to locate the high potential zones for gold mineralization in the study area; Table S1: Eigenvector loadings of principal component analysis for Landsat-8 OLI data; Table S2: Eigenvector loadings of principal component analysis for ASTER data; Section S1: Data and Processing Techniques Funding: The present work has received no funding or grant.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.