Geological mapping and mineral prospectivity using remote sensing and GIS in parts of Hamissana, Northeast Sudan

In recent years, various geological activities and different mineral prospecting and exploration programs have been intensified along the Red Sea hills in order to elucidate the geological maps and to evaluate the mineral potentials. This study is therefore aimed at testing the viability of using remote sensing and geographic information system (GIS) techniques for geological mapping and prospecting for gold mineralization in the area. The study area is located in northeast Sudan and covers an area of about 1379 km2. Different digital image processing techniques were applied to Landsat 8 Operational Land Imager image in order to increase the discrimination between various lithological units and to delineate wall rock alteration which represents target zones for gold mineralization. Image sharpening was performed to enhance the spatial resolution of the images for more detailed information. Contrast stretching was applied after the various digital processing procedures to produce more interpretable images. The principal component analysis transformations yielded saturated images and resulted in more interpretable images than the original data. Several ratio images were prepared, combined together and displayed in RGB color composite ratio images. This process revealed the existence of alteration zones in the study area. These zones extend from the northeast to the southwest in the acid meta-volcanic and silica barite rocks. The enhanced satellite images were implemented in the GIS environment to facilitate the final production of the geological map at scale 1:400,000. X-ray fluorescence analyses prove that selected samples taken from the wall rock alteration zones are gold-bearing.


Introduction
Remote sensing and GIS techniques have been used for geologic interpretations with resounding success. Remote sensing techniques are generally employed because they provide a cost-effective approach resulting from their ability to access difficult terrains and landforms (such as some mountains and forest terrains) and data collection can be done rapidly at frequent intervals on a large scale. This thus opened a new era in the mapping of lithologies and alteration zones. Previously, geological maps were developed from field observations derived from the conventional ground survey studies. The plotting of these maps and extrapolation of details were quite rigorous, leading to some unavoidable errors and inaccuracies. The advent of remote sensing and GIS operations has resulted in continuous procedural changes in mapping activities. Now remote sensing techniques play an important role in mapping programs (Farina et al. 2005). Remote sensing is the science of acquiring, processing and interpreting images and related data, acquired from aircraft and satellites that record the interaction between matter and electromagnetic energy (Sabins 1997).
Using data acquired from remote sensing, geologists have focussed on global experiences in several fields of earth science such as environmental, mineral ore and hydrocarbon exploration (Aminzadeh and Samani 2006;Elsayed Zeinelabdein and Albiely 2008;Elsayed and El Nadi 2014;Omer and Elsayed Zeinelabdein 2018;Petrovic et al. 2012). The mapping of hydrothermal alteration is based on the examination of the spectral signatures of the mineral veins. It is now known that certain minerals associated with hydrothermal alterations have unique spectral features that allow them to be remotely identified and detected (Modabberi et al. 2017). Lithological mapping and the recognition of hydrothermally altered minerals through the application of remotely sensed data have given encouraging results in the exploration of gold, magnetite and other ore deposits (Amer et al. 2012;Azizi et al. 2010;Rajendran et al. 2012Van Ruitenbeek et al. 2012.
Remote sensing provides a great merit of synoptic view that is unachievable with most other geophysical prospecting methods. There is also the advantage associated with giving a regional and integrated perspective of the relationships among various land features.
Nowadays, multispectral and high-resolution data are available with digital image processing techniques with advanced capabilities. These have further enlarged the potential of remote sensing in delineating geological structures and lithologies by generating enhanced and interpretable images with a better accuracy level (Drury 1987). The existing multispectral satellite systems are designed and programmed to study natural resources with special focus on vegetative cover, lithology and mineral/ore exploration (Crippen and Blom 2001;Yousif and Shedid 1999;Abrams 1984).
Several applications of remote sensing in geological studies involve delineating structures/structural features, mineral resource exploration rock and soil type discrimination. (Kruse and Dietz 1991). The launch of Landsat Thematic Mapper (TM) in 1982 gave earth scientists access to better spatial (30 m) and spectral resolution (Abrams 1984), compared to the multispectral scanner (MSS) used, for detailed geological studies (Podwysocki et al. 1983). The Landsat TM and ETM+data have been used in various geological surveys for lithological discrimination, lineaments and mineral mapping by using hyper-spectral laboratory (Abrams 1984). Multispectral remote sensing is essential in lithological mapping and zones of alteration, especially in regions of exposed bedrock and bedrock surfaces.
Studies conducted in Southwest USA and elsewhere (Prost 1980;Rowan and Kahle 1982;Podwysocki et al. 1983;Abrams 1984) confirmed that the use of the ratio of TM bands 5 and 7 could make areas of hydrothermal alterations distinguishable.
Other studies utilized band ratioing method to delineate alteration zones around eastern Red Sea Hills (Elsayed Zeinelabdein and Albiely 2008), while lithological mapping was accomplished through image processing techniques in North Sudan (Ali et al. 2012). Remote sensing, GIS and geochemical techniques were used in prospecting for gold mineralization in North Kordofan, Sudan, with results revealing favourable hydrothermal mineral alteration signals, and the presence of auriferous host rocks (Elsayed Zeinelabdein et al. 2020).
These studies reveal that the spectral characteristics of TM bands make them well suited for lithological mapping.

Location, accessibility and geography
The area studied is in the west of Hamissana, Wadi Edom, which is located in the northwestern flanks of the Red Sea Hills bounded between latitudes 20° 22`N and 20° 50`N and longitudes 34° 00`E and 34° 45`E. The total area studied is approximately 1379 square kilometres (Fig. 1). It forms a part of the Gabgaba terrain that is known to be created among the other four terrains of Northeast Sudan (Tokar, Heya, Gebeit and Gerf) during the Pan-African event (950-500 Ma) due to subduction-related processes and collisional tectonics sutured together by dismembered-ophiolite belts and major shear zones.

Topography and drainage system
The study area (Wadi Eldoum) lies in the northern part of the Red Sea hills topographic terrain, which generally ranges in height between 500 and 1500 m above sea level with some isolated peaks that reach up to 2000 m.
The drainage system of the area is irregular, having dendritic pattern (Fig. 2), and seems to be structurally controlled by shear zone, as well as faults, major joint and zones of weakness such as foliations. It is represented by many Wadies and Khors (seasonally intermittent streams). They have many tributaries which come from hills and high relief areas.

Climate and vegetation
The prevailing climate is the semi-desert type, with very hot and dry summer (March-October), during which day time temperatures can reach up to 46 °C. Winter season is from November to February when day temperature averages 25 °C.
The rainy season in the Red Sea Hills's inland is usually between July and September. The annual rainfall of the whole region is between 100 and 190 mm, Evaporation rate is very high and runoff is rapid. Geological field investigations are best undertaken during the winter season.  Vegetation is poor and scarce. Bushes and small trees such as Acacia Nilotica, Acacia Etabacia and Xerophytes are confined to the seasonal water courses and some of the sandy plains. Short-lived grass flourishes following the occasional rains on the hill sides and in the plains.

Materials and method
Different data, materials and methods of data processing have been utilized in the study, and they are illustrated as follows.
Operational land imager (OLI) and thermal infrared sensor (TIRS), Landsat-8 (path 173/Raw 46) acquired on April 30, 2015, were downloaded from the United States Geological Survey (USGS) in a digital format with 16-bit pixel values characteristics. It contains 11 bands. Nine bands with OLI (from band 2 to band 9) were carried in order to define the main lithological units, major structural elements, drainage pattern and to delineate the zones of gold mineralization. Digital Terrain Model (DTM, 90 m resolution) of the area was obtained to enhance the interpretation process and to delineate the drainage system and lineaments. Digital image processing was carried out using the ERDAS IMAGINE ® version 14.00 build 350, while ArcGIS 10.3 was used for data integration and geological mapping.
Visual interpretation of digitally enhanced images was done; different data types were integrated with the ERDAS IMAGINE ® environment and then with the geographic information system (ArcGIS 10.3), in order to generate the final geological map and delineate zone of mineralization.
Thirty-six chip samples were collected from quartz veins and mylonitic rocks. Other rock types were collected for geochemical investigation and analyzed using the XRF method. Twenty-five rock hand specimens were collected for thin and polish sectioning of representative samples for petrographic work (using an optical microscope). This is in addition to 52 rock samples collected which represent different selected localities for petrographic interpretation.
Collections of ground control points, rock sampling and structural elements were recorded by GPS receiver during these geo-traverses.

Digital image processing
This is the manipulation of digital data by computer programs in order to improve the appearance of an image (Gibson 2000). The data used in digital image processing of the present study are the Operational Land Imager (OLI) and the Thermal Infrared Sensor (TIRS). The spatial resolution of OLI multispectral bands is 30 m, the panchromatic band is 15 m, and the thermal band is 100 m. The thermal infrared band has not been used because of its coarse spatial resolution. The software ERDAS IMAGINE is a revolutionary image processing system used in processing the abovementioned satellite data to enhance the quality of their raw digital data and produce images suitable for geological interpretation.
Image interpretation entails distinguishing features based mainly on tonal and textural variations. The following digital image processing techniques outlined below were used.

Image preparation
The image preparation for the present study included the following steps:

Image subseting
Using ERDAS IMAGINE ® software, spatial subseting was performed on the raw digital images including band 2, 3, 4, 5, 6 and 7 in order to spatially resize the image to the following coordinates: longitude 34° 00ʹE-34°45ʹE and latitude 20° 22′N-20°50`N. This subset covers the study area and helped in tracing major structural elements over a wider regional area.

Resolution enhancement
The resolution of an image refers to the possible detail reflected by the imagery. This could be in the form of spatial resolution, spectral resolution and temporal resolution. In this study, image sharpening was used to automatically merge the low-resolution multispectral band with the highresolution greyscale band (with resampling to the high-resolution pixel size).

Contrast stretching
Contrast stretching modifies the grayscale to produce a more interpretable image (Sabins 2000). The number of actually recorded intensity levels in the scene is rather low, and the full dynamic range of the digital image (256 levels) is not fully utilized. It is worth mentioning that contrast stretching was applied in this study after various digital processing procedures to enhance the resulting image, and was not applied as a separate enhancement process. This process basically increases the spatial frequency of the imagery by enhancing such high-frequency features like lineaments.
linearly uncorrelated variables, called principal components (El Zalaky et al. 2018). The transformation of the raw data using PCA can result in new principal component images that are more interpretable than the original data (Jensen 1996). The PCA is used to compress the information content of a number of bands of imagery or to reduce the dimensionality from a number of bands in two or three PCs (Jensen 1996). In continuation, PCA is extensively applied in remote sensing to determine alterations in rocks and vegetation in a given area. It is also used to reduce the redundancy of information in highly correlated image set.
The flowchart diagram of the adopted methodology for preparation of the geological map of the study area and map of gold mineralization zones is shown in Fig. 3.

Color composite images
Different color composites were constructed during the course of the study. High spectral resolution is important when producing color composite images. The rule of color composites is to set the most informative band for a particular purpose in the red, the next in green and the least informative band in blue filters (Drury 1993). The different RGB combinations discriminate the rock types which is useful in geological application. Any three bands of the six OLI bands can produce color composite image.
A false color composite of bands 7, 5, 3 in RGB, respectively, was assembled in this work as shown in Fig. 4.
On this image, the sheared acid meta-volcanic appears in dark brown color and has linear shapes with the general trend of northeast-southwest (NE-SW), while the basic meta-volcanic appears in light brown color. The silica barite is light green in color and the intrusive granite appears brighter and has spot shape which makes it easily discernible from other rock units and has sharp contact with the adjacent rock units. Another color composite was prepared using the infrared bands of the image set. It was composed by employing bands 7, 6, 5 in RGB, respectively, as shown in Fig. 5. This band combination is termed the infrared image.
Using the six reflective OLI multispectral bands, principal component analysis (PCA) was performed over the study area.
New axes of principal component analysis were computed statistical analyses. This procedure resulted in the production of six principal component images illustrated in Fig. 6 The first PC image contains the largest percentage of data variance, and by comparison, the variance in last PCs image is seen to be negligible due to the fact that the image appears totally noisy and contains no information.
Principal components produced more colorful color composite images than spectral color composite images because of the uncorrelated nature of the obtained PC images.
PC1 and PC2 display more lithological contrast, the topographic expression is better, and it is known that they discriminate well between the visible and near-infrared (VNIR) and shortwave infrared (SWIR) bands. Although PC3 and PC4 have low variance values, they still display fair lithological contrast, whereas PC5 and PC6, with very low variance values, are less informative and show more noise.
The principal component color composite (Fig. 7) was created by assaying the PC1, PC2 and PC3 to the RGB filters, respectively. Composite images of PC1, PC2 and PC3 in RGB and PC2, PC3 and PC4 in RGB have provided much lithological information and discrimination between rock units.

Rationing technique
Image ratioing is one of the change detection techniques used in remote sensing. Change detection is the process of identifying differences in the state of an object or phenomenon by observing it at different times (Singh 1989).
Spectral band rationing is a proven technique which allows identification of geological materials based on the reorganization of diagnostic absorption bands. It minimizes the effect of topographic slope, aspect, and Albedo differences between rocks and enhances the subtle differences in reflectivity between bands which are diagnostic of various surface materials (Estes 1983). Ratio images are prepared by dividing the DN (digital number) value in one spectral band by the corresponding DN value in another band for each pixel (Sabins 2000;Afify 2011). Band rationing has been used successfully with Landsat data for identifying the hydrothermally altered rocks, particularly volcanic rocks, Ratio images are known for enhancement of spectral contrasts among the bands considered in the rationing and have successfully been used in mapping alteration zones (Segal 1983;Kenea 1997). During this study, several ratio images were prepared for the purpose of geological mapping and alteration zones delineation. Generally, the OLI 4/2 image ratio shows high values for iron stained hydrothermally altered rock and ferric oxides because iron minerals have low reflectance in band 2 and high reflectance in band 4, while the OLI6/7 ratio is used to enhance the hydroxyl-bearing minerals and distinguishes altered rocks containing clays and alunite from unaltered rocks. Several OLI bands ratio have been computed: 6/7, 6/5, 6/2, 4/2, 5/4 and 4/5.
Based on the spectral features of ferric and hydroxylbearing minerals, two ratio false color composite images were produced by combining three ratio images at a time in RGB. The different combinations of ratio images are shown in Table 1, and this provides the computation suitable for use with Landsat 7 ETM+images and Landsat 8 OLI images. Sultan et al. (1987) used Thematic Mapper (TM) 5/1 ratio (equivalent to OLI 6/2 ratio images) because opaque minerals have reflectance features in band 5 and absorption features in band 1. Sultan's color composite ratio image was obtained by using OLI image ratios which corresponds to TM image ratios as 6/7, 6/2 and 4/5*6/5 in R, G and B, respectively. The resultant image (Fig. 8) shows that the acid meta-volcanic rocks appear bluish in color, the metasediment reddish in color, the sheared granite in deep blue color, the post-tectonic granite appears brighter, and the clay minerals appear in pink.
Mineral's composite ratio image was obtained using OLI band ratios 6/7, 7/5 and 4/2 in RGB, respectively. The obtained image (Fig. 9) shows meta-volcanic in reddish, while the meta-sediment is in light green color and easily discernible from other rock units. The alteration zone is mapped in reddish yellow.

Geology
The study area is covered by different lithological units with predominant acid meta-volcanic rocks and syn-tectonic granite in addition to various post tectonic rocks. There is also the presence of few scattered low laying outcrops of  sedimentary rocks and superficial deposits. Figure 10 is a map showing the regional geology of the study area.
From the field description and geochronological result, the study area is commonly characterized by Mylonites rocks that are called generally meta-volcano sedimentary/   Figure 11 shows photomicrographs of quartzite and marble, representing the meta-sediments. In the same vein,

Structure
The most pronounced structural features in the study area are the presence of faults. They have the general trends of northeast-southwest (NE-SW), northwest-southeast (NW-SE) and east-west (E-W) (Fig. 15). Most of these faults are strike-slip where the dominant displacement is horizontal and parallel to the strike of the fault. The evidences of these faults are the rocks brecciation, shearing, ferruginization and slickensides.

Alteration zones
The breakdown of feldspars and ferromagnesian minerals to a variety of clays and other hydroxyl bearing minerals is the commonest type of alteration. Such minerals can be detected by remote sensing techniques since the shortwave infrared (SWIR) range of their spectra exhibit absorption (Drury 1993). It has been well established that ferric iron exhibits pronounced absorption features at around 0.82 µm and 1.3 µm (infrared region), whereas ferrous iron has absorption at 1.0 µm, 4.8-2.0 µm and 0.55-0.45 µm (Ultra-violet and violet range). Hydroxyl-bearing minerals such as clays have a major absorption feature around 1.9 µm, 2.35 µm and 2.5 µm (band 7) (Gupta 2003). Detection of the above mentioned minerals has been used as prospecting guide for ore deposits.
To fulfil the objectives of the present study, alteration zone mapping was conducted using Landsat OLI image applying different digital image procedures, such as rationing technique and feature oriented principal component analysis. The results of the different techniques utilized to this end are presented in this section.

Mapping alteration zones using rationing techniques
Extensive hydrothermal alteration zones and weathering of the sulphide mineralization within the acid volcanic represent a significant mineral province in the Red Sea Hills of Sudan (Robertson Research International, RRI 1985). False color composite image called the hydrothermal composite because of its ability to enhance hydrothermal wall rock alteration was produced using combination of bands 5/7, 3/1, 4/3 in R, G, B, respectively (Fig. 16).
On this image, the altered rocks appear in reddish-pink colors. These alteration zones are observed along linear structures trending NE-SW in the acid meta-volcanic rocks.

Mapping alteration zones using feature oriented principal component analysis
This method was created by Crosta and McMoore (1989) with the purpose of using remote sensing satellite images for geological exploration in an environment lacking bedrock exposure and whereby selected bands believed to contain Fig. 16 Ratio image color composite produced using the following ratios 6/7, 4/2, 5/4 in R, G, B, respectively. Red circles are samples collected for geochemical analysis For enhancement of both iron-oxides and hydroxyl-bearing zones, two sets of four TM bands were selected: bands 1, 3, 4 and 7, and 1, 4, 5 and 7. This selection is based on the fact that the former have contrasting signatures in TM bands 1 and 3, and the latter in TM bands 5 and 7 (Loughlin 1991). The above described two PC transformations were conducted for the study area. Table 3 summarizes the basic statistics of the bands involved in the transformation and gives the eigenvector loadings of both band sets. As it is clear from Table 3, PC4 has mapped iron-oxides as negative loading (− 0.77) in TM-3 and as (0.13) in TM-1. For hydroxylated minerals (Table 3), PC4 has relatively more contrasting signature loading (− 0.69) in TM-5 and (0.70) in TM-7.
PC3 was selected from the first transformation to map the iron oxide. PC4 was selected from the second transformation to map the hydroxyl-bearing minerals. A linear combination of the two PCs was produced to map the hydroxyl-bearing minerals and iron oxides. A composite image was prepared by combining the "hydroxyl," hydroxyl + iron-oxide and iron oxide images in RGB, respectively. A low pass spatial filtering of 3 × 3 kernel was run in order to remove the incorporated significant noise and improving the hue by removing isolated pixels smaller than the applied window size. The obtained image is displayed in Fig. 17. On this image, all alterations have been mapped in brownish golden color.

X-ray fluorescence (XRF) analysis
Groundtruthing and field check-up were made to authenticate the results obtained by the alteration zone mapping through satellite image interpretation. Geochemical samples (Table 4) were collected from quartz vein and alteration host rocks.
The geochemical analysis results show that the highest anomaly value for gold is 43.913 ppm and the lowest is 0.009 ppm. Results from the work of Elsayed (2010) from five outcrop surfaces around the Red Sea Hills showed gold content ranging from 1.1 ppm to 12.2 ppm in sample AZ2. It is clear from the results of the geochemical analysis as well as geological mapping and alteration zone mapping that the area is highly prospective for gold. Therefore, this area is extremely encouraging for further ground survey and exploration work.

Conclusion
Remote sensing technique has proved a powerful method for large scale geological mapping in relatively short time and with a reasonable degree of exactness. This mapping is of great importance since the exploration of any given area must start with geological mapping. However, care must be taken in the usage of satellite images as base maps Fig. 17 Principal component analysis color composite obtained by assigning hydroxyl," hydroxyl + ironoxide and iron-oxide images in R, G, and B, respectively (Red circles are samples collected for geochemical analysis) for geological mapping as different rock types have similar reflectance properties.
The use of remote sensing in such areas as the Red Sea Hills must be considered as a major and efficient tool as these areas are characterized by hilly topography and harsh environment. This makes the conventional field survey difficult if not impossible. Remote sensing results obtained in this study have emphasized the viability of remote sensing technique in geological mapping and mineral exploration. Alteration zones which are normally rich in mineral deposits were successfully mapped using transformation techniques such as band ratioing and Feature Oriented Principal Components Analysis (FOPCA). The outcome of these transformations shows that the area has potential for mineral exploration and exploitation. This process revealed the existence of alteration zones within the study area.
Geochemical analysis proves the altered rocks to be gold (Au) bearing. This is in general, an encouraging result for further exploration in the outlined area, since it is known from literature that the mineralization zones in the Red Sea Hills have the same linear structures, in addition to their very existence in the acid meta-volcanic rocks as the alteration zones in the study area.
For future studies, more detailed fieldwork should be conducted especially in the outlined alteration zones to verify the remote sensing results and to collect good amounts of field measurements to enhance the geologic understanding of the area and to re-correct the final geological maps. Also, a detailed structural analysis of the lineaments of the study area should be carried out. Future research could incorporate high (spectral and spatial) resolution satellite data to map certain windows of geologic and/or target localities of economic significance in the area. This will support the production of more detailed geological maps of suitable scales.
Funding The authors received no special funding for this research.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.