Mapping Groundwater Potential Zones Using a Knowledge-Driven Approach and GIS Analysis

Despite the Sahara being one of the most arid regions on Earth, it has experienced rainfall conditions in the past and could hold plentiful groundwater resources. Thus, groundwater is one of the most precious water resources in this region, which suffers from water shortage due to the limited rainfall caused by climatic conditions. This article will assess the knowledge-driven techniques employed to develop a model to integrate the multicriteria derived from geologic, geomorphic, structural, seismic, hydrologic, and remotely sensed data. This model was tested on the defunct Kom Ombo area of Egypt’s Nile river basin in the eastern Sahara, which covers ~28,200 km2, to reveal the promising areas of groundwater resources. To optimize the output map, we updated the model by adding the automated depression resulting from a fill-difference approach and seismic activity layers combined with other evidential maps, including slope, topography, geology, drainage density, lineament density, soil characteristics, rainfall, and morphometric characteristics, after assigning a weight for each using a Geographic Information System (GIS)-based knowledge-driven approach. The paleochannels and soil characteristics were visualized using Advanced Land Observing Satellite (ALOS)/Phased Array type L-band Synthetic Aperture Radar (PALSAR) data. Several hydromorphic characteristics, sinks/depressions, and sub-basin characteristics were extracted using Shuttle Radar Topography Mission (SRTM) data. The results revealed that the assessed groundwater potential zones (GPZs) can be arranged into five distinctive groups, depending on their probability for groundwater, namely very low (6.56%), low (22.62%), moderate (30.75%), high (29.71%), and very high (10.34%). The downstream areas and Wadi Garara have very high recharge and storage potential. Interferometry Synthetic Aperture Radar (InSAR) coherence change detection (CCD) derived from Sentinel-1 data revealed a consistency between areas with high InSAR CCD (low change) that received a plausible amount of surface water and those with very low InSAR CCD values close to 0 (high change). Landsat data validated the areas that received runoff and are of high potentiality. The twenty-nine groundwater well locations overlaid on the GPZs, to assess the predicted model, indicated that about 86.17% of the wells were matched with very good to moderate potential zones.


Introduction
Water is a valuable resource necessary for life and sustainable development [1]. Affording assured water resources for agricultural and domestic activities is a seemingly tough task in areas of scarce water sources and uneven rainfall distribution [2,3]. Groundwater is an important source of freshwater, utilized for sustainable development on Earth, particularly in areas with arid/hyperarid zones. A report from UNESCO [4] stated that about 2.5 billion people around the world depend merely on groundwater sources in many activities, such as domestic use and industrial and agricultural purposes [5]. The need for freshwater resources is likely to grow and, by 2030, the level of water deficiency will reach about 40% globally [6]. In Egypt, irrigation with groundwater is an important source for farmers in desert fringes, a vulnerable area due to the risk of water scarcity. Although building massive water infrastructures yields significant benefits and is a priority in several regions around the world [7], the construction of the Ethiopia Dam has made the situation worse, and it may cause a significant water crisis for the downstream countries of the Nile River, as the populations are largely dependent on its water. Moreover, the demand for water resources has increased in the past few decades due to corresponding increases in population and economic and agricultural activities.
Groundwater is a hidden natural resource beneath the Earth's surface [8] that requires integrated approaches to predict and find [9][10][11][12][13]. Conventional field methods are often expensive and require a long time to implement. Since the advent of satellite remotely sensed images and Geographic Information Systems (GIS) procedures, developments in groundwater reconnaissance and monitoring studies [2,11,13,14] have occurred in arid/hyperarid environments. Remote sensing data have provided significantly rapid, cost-effective information on the geologic, structural, and hydrologic features that have guided us to groundwater conditions and their occurrences [9][10][11][12][13]. Developments in data analysis techniques have also promoted imaging in progressively finer image resolutions. Radar images have provided useful information and helped in delineating structural features and stream networks, along with estimating the rainfall rate (e.g., [2,13]). GIS techniques have allowed for the combination of multi-evidential layers that have contributed to the revealing of favorable areas for groundwater recharging.
Groundwater accumulation, infiltration, and movement depend on several factors, such as lithology, recharge conditions, structures, geomorphology, climatic situations, streams, and geologic sorts. Several studies (e.g., [1,5,9,12,13,[15][16][17][18][19][20][21][22][23]) utilized remotely sensed images/data and GIS for groundwater infiltration, recharging, and exploration, as well as the characterization of artificial recharge zones in many areas all over the world. Crucially, remotely sensed data is crucial to the characterization of hydrogeological properties and conditions in many environments. However, even in well-mapped areas, satellite images reveal lithologic units as well as geomorphologic and structural features that were not recognized via conventional approaches. These data enhance our information on hydrogeological conditions [24,25]. Knowledge-driven approaches that allowed for the weighting and classifying of evidential maps, based on prior knowledge and information on the infiltration capacity of the utilized layers, have been successfully performed in delineating groundwater potential zones [5,13,20]. The aforementioned studies used common evidential layers (n = 4 to 8), such as geology, lineaments, slope, drainage density, rainfall, and topography. In the present model, we selected 10 thematic layers and added two layers to the model, including seismic activities and depressions, rather than analogous previous studies.
The desert fringes of the Kom Ombo area are promising areas for land cultivation and population growth. They are distinguished by excessive water claim and pressure on existing water stores. The objectives of the present investigation are as follows: (a) to develop a model to test the ability of knowledge-driven techniques through GIS, to delineate prospective areas of groundwater resources after ranking the multicriteria; (b) to test the extracted depressions and seismic activity for groundwater potential of the present study.

Study Area
The study area is located east of the Kom Ombo area; it presently drains westward to Kom Ombo town and receives surface water from the highlands of the Red Sea in random precipitation events. The wadis are geomorphic characteristics in the eastern desert of Egypt, covering the~28,200 km 2 that extends between latitudes 22 • 4 47.88" to 25 • 16 23.48" N and longitudes 32 • 57 28.24" to 35 • 9 13" E ( Figures 1 and 2). The upper reach collects water from the basement catchment and drains to the downstream of the Nile basin, at Kom Ombo city. The district is comprised of igneous and metamorphic varieties to the east, and the Upper Cretaceous/Lower Tertiary succession to the west ( Figure 2). The Precambrian rocks, consisting of older metamorphosed rock units, include gneisses, metasediments, metavolcanics, and serpentinites, along with metagabbro-diorites complexes [26] (Figure 2). They are intruded by older granites, younger granites, and younger volcanic rocks. The sedimentary succession comprises Nubian Sandstone that covers a wide area and overlays the basement rocks [27]. Upper Cretaceous and Paleocene rocks form the western part of the present catchment. The area is affected by the NW-SE, NE-SW, NNW-SSE, and E-W faults. The N-S trend represents the alignment of the basement extension along the contact. Based on hydrologic conditions, the area is located in an arid/hyperarid zone, with little precipitation during the rainy season, much like northern African countries (Figure 1a).

Data Used and Methods
Remote-sensing data of various spatial, spectral, radiometric, and temporal resolutions have the capability to afford accurate, cost-effective, automated, near-real-time information, even in areas on the Earth that are difficult to access [1,2,7,9]. Technological advances in image processing and analysis have allowed for the extraction and combination of information in a fast manner that can help improve decision-making. It is a very significant approach to reveal the geologic, structural, and hydrologic conditions for better assessment, planning, and observation of water resources under any conditions [2].
Remote-sensing data, comprised of multispectral and radar, plays a significant role today in mapping hydrologic conditions in different climatic conditions. Multispectral remote-sensing data, comprised of information from within the visible-near-infrared and shortwave-infrared regions, utilize the land-use/cover for mapping. The data includes microwaves and visible/IR waves that are processed and integrated through the Environment for Visualizing Images (ENVI; Harris Geospatial Solutions, Broomfield, CO, USA) v5.3 and ArcGIS 10.5 software packages (ESRI Company, Redlands, CA, USA). In addition to remotely sensed data, the geological map [26] was digitized to prepare a geologic map of the present study. This integration assists in understanding the geological and structural features involved in recharging the aquifers. The data adopted for the present study are given in Figure 2.
Microwave data are significant in regional studies, such as the study of sand-buried drainages and in structural features (e.g. [28,29]). Microwave data of ALOS/PALSAR, Sentinel-1, and SRTM data are beneficial here to understand the geometry, geomorphic, and morphotectonic features. In addition to precipitation derived from TRMM records for estimating the rainfall, precipitation will be utilized in the present study between 1 January 1998 and 30 November 2013. The data were accessed from the NASA website (https://giovanni.gsfc.nasa.gov/giovanni/, accessed on 10 February 2021).
The Japanese active microwave satellite sensor ALOS/PALSAR is an advanced landobserving satellite remote-sensing instrument, frequently utilized to monitor and observe the land at any time, in totally meteoric circumstances. Such a sensor employs the L-band with a spatial resolution of 25 m and a swath width of 70 km; a mosaic of Jaxa "2017" was employed in this study. Although the data are available with the HH and HV polarizations, the present investigation employed the HH polarization.
The SRTM data are part of an international joint study conducted by National Aeronautics and Space Administration (NASA), National Geospatial-Intelligence Agency (NGA), and the Italian Space Agency (ASI) and the German Aerospace Center (Deutsches Zentrum für Luft-und Raumfahrt; DLR). Such data were used to present the topographic features of the Kom Ombo area. The data can be accessed from the United States Geological Survey (USGS) with a~30-m horizontal resolution. The digital elevation model (DEM) derived from the SRTM data was processed and visualized using GIS approaches. Stream networks were extracted by using surface flow routing constructed by the 8D [30], using the spatial analyst approach. The extraction of drainage and hill-shade relief is suitable for identifying the structural features controlling the terrain and groundwater aquifers. The "Fill" tool in ArcGIS v. 10.5 allowed for the defining of depressions/sinks after subtracting from the original DEM to reveal the "fill-difference". So, this can be represented in the matrix form as the original DEM layer, indicated by the "O" matrix, and the filled-DEM layer indicated by the "F" matrix. Let O and F be two matrices of the same order (m*n). Let O = O ij and F = F ij . Therefore, (O) − (F) is a matrix of the same order as O and F; its elements are obtained by subtracting the elements of F from the corresponding elements of O. Then, the result matrix S ij can be estimated as follows: Two scenes of Sentinel-1 were acquired on 19 September 2020 and 18 September 2018. VV polarizations were utilized to characterize an image, revealing the changes in land features. The phase and intensity data of the single look complex (SLC) SAR format were used to characterize the coherence change detection (CCD) information and define the activities of the active streams and surface land features [31]. The InSAR CCD has revealed land surface changes and activities in the spatiotemporal span [31,32]. In this method, every pixel of a SAR includes intensity and phase information. Mathematically, the interferometric coherence (G) of the two SAR scenes can be computed as in [33] (Equation (2)).
where Sc1 and Sc2 are the complex signals of the "master" and "slave" scenes. The means of the kernels with L × M pixels can be used to estimate the coherence of an output pixel [31]. The result of the coherence (G) map extends from 0, for complete change, to 1, which defines no change [34]. The Landsat Operational Land Imager (OLI) optical sensor with 9 spectral bands carried by Landsat-8 was launched on 11 February 2013. Two scenes of Level 1 (path/row 174/042 and path/row 174/043) OLI data were acquired on 9 September 2020. The obtained scenes were mosaicked and linked to the coordinate system of UTM Zone 36N. The OLI bands 2, 3, 4, 5, 6, and 7 were utilized here to conduct image transformation and enhancement techniques. The data were radiometrically calibrated using ENVI v.5.3. A scene of a Landsat-7 Enhanced Thematic Mapper (ETM+) image, acquired on 21 January 2010, was also used to visualize the flood event in 2010 over the study area. To monitor the changes in vegetation, the Normalized Difference Vegetation Index (NDVI) was carried out using Landsat series, Landsat-OLI (8 September 2020), Landsat ETM+ (10 September 2000), and Landsat-5 (15 November 1986) of the same area covered by path/row 174/043.
The GIS technique was used for analyzing the different spatial data by applying various logical conditions. Lineaments were derived from high-resolution imagery by digitization using the GIS technique and analyzed using Rockworks-16 software to extract the rose diagram. The geological map was digitized using a map from [26], prepared from a thematic layer of geology. Spatial modeling of prospecting groundwater potential zones (GPZs) using multicriteria was applied using knowledge-driven techniques.
The knowledge-driven predictive overlay model [35,36] was applied using a geographic information system (GIS). This model is based on known data, and the relative importance of each observation is decided on by the user. The GIS approach uses the raster overlay, in which each pixel of each layer has the same geographic location. This makes the overlay process more appropriate for integrating characteristics for several data sets into an output layer. The numeric values given for each layer promote the mathematical combination process and give a new grade to each pixel in the output cell. The weighted values are then given to the spatial predictor maps and their categories.
Based on the knowledge-driven technique of combining multicriteria [37], the scores of each evidential map are assigned to different maps and their intermap features. In this model, the evidential maps are weighted from 1 to 4 and the given weight was normalized to control the occurrence and infiltrate the groundwater. The intermap feature classes were also assigned weights and normalized in every evidential map. The GPZ, which represents the weighted average of the combined data-based maps [37], can be delineated using the following equation: where Wi refers to the evidential layer weight and Fi corresponds to the weight of the intermap class.

Results
GIS-based approaches, remotely sensed data, and supplementary geologic-geomorphic models were merged to delineate potential areas of plausible water resources of the selected area in Egypt. A knowledge-based method to combine the multicriteria of several factors (e.g., geology/geomorphology, topography, slope, lineaments, depression, seismic activity, drainage density, rainfall, and alluvial deposits) was derived from remotely sensed data, as the existence, infiltration, accumulation, and movements of groundwater are controlled by the aforementioned parameters that serve as inputs to the thematic layers.
The information resulting from the remotely sensed data and geologic-geomorphic features can be integrated with the GIS technique to find the GPZ using Equation (2). Groundwater cannot be perceived directly from the Earth's surface, so a variety of approaches are required to provide information regarding its potential occurrence and characteristics. Such sources of water are necessary because of their clear color and the absence of pathogenic organisms and turbidity. Furthermore, their stable chemical composition and storage capacity are significant.

Geology/Geomorphology
Lithologic characteristics are significant, as the infiltration of the precipitated water to the groundwater aquifer is governed by the porosity and permeability of the surface layer [1,5]. Based on processing and interpretation of the remotely sensed data, the area of east Kom Ombo consists of different rock units. The crystalline rocks occupy the eastern part of the present study. The Cretaceous/Tertiary succession occupies the western part of the present study. The basement rocks, sedimentary succession, Nubian sandstone, and Quaternary deposits cover about 51.23, 0.87, 34.20, and 13.68% of the study area, respectively. The wadi bed built-up of sand and gravel deposits increases the porosity, which is a good factor for groundwater recharging [38], versus the basement rocks comprised of impermeable rocks of granitoids, metavolcanics, and serpentinites ( Figure 3). The abundance of sand deposits related to the sandstone exposures in the study area support the long-term goal of infiltration capacity in the basin, controlled by structural features. Based on the geologic map, the evidential geologic layer is categorized into four classes (Table 1)

Topography
The topography layer is an important layer that governs the direction of water flow over the land; it also controls the occurrence of groundwater and recharge potential [12,13]. The wadi is surrounded by elevated hills and the main catchment is the Red Sea Highlands. In this image, the low topographic areas collect, hold, and store much more water than the high topographic areas during the rainfall season. The low altitudes yield water recharge and movement. The DEM of the study area was acquired from the SRTM data (Figure 4a). Analysis of the SRTM DEMs data revealed that the land-surface height of the area varies from 91 m to 1710 m. The topographic values revealed high runoff in an east-west direction. The area of high elevation occupies the eastern, northern, and southern parts of the studied basins. Interpretation of the DEM data revealed a network of fracture systems consistent with the structural elements that were probably incised in the plateau. Based on the infiltration capabilities, the elevations were sorted into five groups based on their tendency of water storing (Figure 4b). The areas of very low, low, moderate, high, and very high topography cover about 1.50, 11.70, 31.85, 37, and 17.90% of the area, respectively ( Table 1). The low-topography areas promote water accumulation and infiltration capabilities.

Slope
The slope is likely to affect the water flow, occurrence, movement, and infiltration capacity of groundwater [19,39]. The slope has an inverse relationship with the infiltration of the surface water to the aquifers underneath [13,40]. The slope can often be influenced by the rock types, rate of erosion, climate conditions, and tectonic activities. The slope map is derived by the SRTM DEM. It varies between 0 • and 50.84 • (Figure 4c) and is grouped into five classes. The low slope degree (flat areas) considers a promising area of surface water infiltration and recharge potential; thus, it is given a high rank, as the slope controls the runoff and vertical infiltration of precipitation [19]. However, steep-slope areas are assigned a low rank and high runoff capacity [9]. The slope direction reflects the dominant westward runoff to the Nile Valley. Based on the slope degree map, the majority of the area (~65%) extends under the gentle slope (<5 • ) surface level that represents the suitable zones for recharging capability (Figure 4d), which allow for a holding time for surface water infiltration ( Table 1). The very high potential degree (5), high (4), moderate (3), low (2), and very low (1) slope cover about 62.15, 24.48, 8.84, 3.33, and 1.18% of the studied basin area, respectively (Figure 4d).

Depressions/Sinks
Depressions are low surface elevations that must be filled with water until the water elevation reaches a level at which the flow would pour out of the depression. Using the Fill tool in ArcMap v. 10.5, the depressions/sinks were filled to their spill elevation ( Figure 5). Subsequently, the "fill-difference" could be calculated by subtracting the initial DEM from the filled-DEM. This allowed for the generation of a raster that signified the depth of the depressions in the original surface [41]. These depressions capture the flow of water, which is accumulated and stored until it infiltrates into the aquifers underneath. In arid regions, such a process can happen during rainy storms. The results of the extracted depression layers are classified into four categories: −82 to −8 (high), −7.99 to -5 (moderate), −4.99 to −2 (low), and -1.99 to 0 (very low), based on the capability of long-term capturing and storing of accumulated water (Figure 6a). The implementation of such data would benefit the GPZ model. The plausible promise area of groundwater infiltration covers about 0.90% of the entire area (Table 1). A large number of depressions were revealed along the streambed channels. These depressions are most likely consistent with specific areas (e.g., upstream of W. Beizah, a tributary of W. Shait (Figure 6b)) that most likely carved the sandstone, and are low-topography-controlled by faults along W. Shait (Figure 6c,d), behind the natural barriers or stream-neck at W. Garara (Figure 6e,f). The extracted depressions are validated using the analogous areas in Landsat-8, ALOS/PALSAR, and the DEM (Figure 6g-i), to understand the uncertainties. The results of validation revealed more than 82% of accuracy.

Stream-Networks
The nature of the stream-networks and catchment area control the harvesting of the precipitation on the highly elevated areas (Figure 7). Catchments that capture excessive amounts of rainfall have the capability to infiltrate and recharge the groundwater aquifers. Several studies have revealed that stream-network density is inversely connected to the rate of recharge processes (e.g., [5,42]), but others [2,23,40,43] have indicated that the high drainage density reveals a high infiltration capability because of the highly dissected land surface. The results of using the D8 method and stream networks were predicted in the DEM. In this article, we related the areas of high density (Figure 7a,b) to having greater suitability for recharging and infiltration capacity. As the majority of the area is a gentle slope, the areas with high drainage density are influenced by high groundwater and low runoff [23]. Therefore, the studied basins can be graded into five categories: 5.14-29.52 (very low), 29

Morphometric Characteristics (Runoff)
Based on morphometric analysis, the perimeter (P) of the studied basin is 1537 km, with a bounded area (A) of about 28,200 km 2 . The wadi was subdivided into 11 sub-basins (Tables 2 and 3; Figure 8). The expressed formula used in the morphometric analysis of the studied basins is referenced in Table. 2 [24]. The A of the delineated 11 sub-catchments varies from 678 to 8632 km 2 . Several parameters were calculated to characterize the prominent expanses of runoff (Table 2), using the linear equations referenced in [24]. Using such equations, the summation of the runoff rate of the sub-basins ranges from 25.02 (low runoff) to 49.71 (high runoff), categorized into five ranks that extend from the lowest to the highest (5). Based on that, the runoff was found to be inversely related to the infiltration capability and the runoff potential map (Figure 8b) was rearranged to infiltration capacity (Figure 8c). The potential areas of infiltration reach 30% of the study area (Table 1).

Radar Intensity
The ALOS/PALSAR data covering the study area (Figure 9a) reveal the variation in the backscatter characteristics of the soil and bedrock units. Such radar intensity data show the nature of the deposits in the study basin, which are observed with difficulty in the Landsat data [43]. The fine deposits depicted in radar data reveal fine sediments and alluvial sand sediments. The coarse-grained texture reveals bedrocks that demonstrate massive characteristics. This is because of the low, diffuse backscattered beams. The wadi bed has backscattered fine deposits. A fused ALOS/PALSAR and Landsat-8 discriminated the geologic and structural features (Figure 9b) and vegetated areas in the downstream, along with patches within the sedimentary rocks (Figure 9c). Discrimination of the fine deposits from the bedrocks and coarse-grained deposits was conducted using the logical operation (classification) of ArcGIS tools. Based on this operation, the fine deposits occupy the majority of the wadi bed and stream areas on the bed rocks of the study area. The dry sand of the wadi deposits, derived mainly from sandstone, provides optimum conditions for the used L-band to penetrate deeply, because of its transparency in radar images. These sand deposits revealed an alluvial origin caused by the erosion of the surrounding rocks. Several fractures were probably hidden throughout these deposits and connected to the groundwater aquifer. ALOS/PALSAR data revealed several areas of dark color representing fine sediments or areas rich in sand deposits; this is because the low backscattered signal, compared to the rough surface, spreads in many directions, and we know a diffuse reflection results in a bright tone on the surface. The radar map was categorized into 5 grades of categories ranging from 1 (low backscattered) to 5 (high backscattered). The results from the radar map revealed that the high values indicate locations for low groundwater recharging and the low values displayed the high GPZs (Figure 9b), which reach 26.85% of the entire area ( Table 1).
The active (ephemeral) streams during rainfall storms have a spectrally bright view compared to the surrounding pixels of Landsat-8 images (see Figure 9c). A fused ALOS/ PALSAR and Landsat-8 (Figure 9d) system probes the sand cover and shows a near-surface view of the paleochannel (cyan arrow) that was drained under sand and is barely visible in the Landsat-8 images (Figure 9c). These images revealed that the NW flow (cyan arrow) was captured by the present SW. However, in the downstream area of W. Kharit, the SW flow is captured by the NW flow (Figure 9e,f), such as Wadi Abu Suberah [44]. The filling wadi deposits appear dark in the ALOS/PALSAR images, versus the intermediate to bright outcrops of surface bedrocks reflecting textural information based on surface roughness (e.g. [9,45,46]. The radar and optical data also provide evidence of the stream capture process of Wadi Abu Suberah (west) by Wadi Kharit to the northwest, resulting from active tectonic movements (Figure 9g,h).

Lineaments
Lineaments are the linear or curved geologic structural features (faults/fractures/joints) that play a significant factor in characterizing the promising areas for groundwater potential. Recently, several studies have focused on using structural features in recognizing groundwater possibilities and explained the role of fracture systems prospecting the availability of groundwater resources (e.g., [9,12,15,47]). Such investigations revealed that groundwater potential is positively related to areas of higher lineament intensity, where using structural features controls the aquifer system, as the damaged zones of faults can be represented by channels or blocks to the flow of the groundwater [48,49]. The permeability along the fault planes is delimited by the faulted rock and time of the fracture/fault systems [49]. The faults have a tendency to be occupied with broken or crushed material (fault gouge) of high permeability that facilitates water movement (e.g., [2,50]).
The fracture/fault zones connected to the aquifers are also noteworthy. Optical/ remotely sensed images provided evidence of such a connection (Figure 10a). Lineaments were extracted using remote sensing data, such as the shaded relief map derived from SRTM data, fused ALOS/PALSAR and Landsat-8, stream-networks, and Landsat-8 ( Figure 10a). Integration of lineaments [26] and shaded relief of DEM data allowed for the preparation of the lineaments map (Figure 10b).
The E-W trend captured the NW-SE trend caused by the tectonic activities of the Red Sea. A fused ALOS/PALSAR and OLI image allowed for the visualization of the geomorphic features and structural continuity (Figure 10c) of W. Garara. This image highlighted the geologic structures and fault system control of the drainage system along the basement rocks. The structural types and movements guided the formation of the porous and permeable zone for possible recharging processes and controlling groundwater flow. The area under exploration was modified by tectonics in the Miocene, which allowed for the breaking of the plateau, contemporaneous with the Red Sea events. Moreover, high-resolution data revealed that several fracture/fault zones control the granitoids of the Arabian-Nubian Shield (Figure 10b).
The lineaments are dominated NW-SE in the extension of the Arabian-Nubian Shield crossing the main wadi channel that rends NE-SW (Figure 10a). The fused images (Figure 10d,f) improved the form of the lineaments' visualization, rather than using individual images of Landsat (Figure 10c,e). The lineaments' density map was grouped into five numerical categories (5 (highest), 4 (high), 3 (moderate), 2 (low), and 1 (lowest); Figure 10b), covering about 23.97, 35.40, 26.91, 12.74, and 0.963% of the study area, respectively (Table 1). These five categories depend on their relative recharging properties. The area of intense density is a promising site for groundwater penetration and recharging comparable to the low-density zone. The rose diagram revealed the major trends, including the NW and SE directions that affect the area; however, NE-SW faults are less abundant. The NW-SE and sub-E-W dominate the basement rocks, and the NW-SE characterize the sedimentary rocks of the area. Noteworthily, fused ALOS/PALSAR and OLI data successfully visualize the structurally controlled paleochannels that are difficult to recognize from OLI data alone (Figure 10g,h).

Rainfall Data (TRMM)
Rainfall is a crucial factor of the hydrological characteristics of a region, as it is the source of runoff and surface water during storms [5]. The average of the rainfall intensity data (mm/day) derived from the TRMM satellite covering the study area (years 1997-2015) provides important precipitation information on the conditions and anomalies of precipitation in different areas (Figure 11a). This data estimates the runoff of surface water and excessive rainfall in the catchments. The excessive rainfall intensities allow for a high recharge capability to the groundwater aquifer, as infiltration depends on the intensity and duration of precipitation. The low value of precipitation reveals the low recharge of the precipitated water to groundwater aquifers. Subsequently, the rainfall data was sorted into four numerical categories (Table 1). These categories (Figure 11b) range from 1 (low) to 4 (high). The results of the precipitation anomalies map explained that the low class of rainfall intensities indicate a lower possibility of infiltration characteristics and recharging groundwater, and the high class displays a high probability of the infiltration and recharging process (Figure 11b). The potential areas of the rainfall precipitation cover 4.77% of the study area, mostly consisting of the downstream and upstream areas. During storms, the precipitated water collected and contributed to runoff on the streams (Figure 11c,d), so the W. Shait received a high amount of precipitation during storms on 17-18 January 2010, which was validated by a bright cyan spectral signature in the band composite 7, 4, and 2 of Landsat-7 (Figure 11e,f). Therefore, these streams that receive precipitation would hold groundwater, as the majority of surface water can infiltrate the sediments below.

Earthquake (Seismicity)
Earthquakes are natural processes that can yield rupturing and deformations in the Earth's crust. In this work, we assume that the deformation of the brittle rocks yields the opening of fractures/fault systems [51] and changes in the bulk permeability and porosity [52,53]. Therefore, areas frequently subjected to earthquakes are probably dominated by ruptures reaching the subsurface aquifers ( Figure 12). The response of the groundwa-ter recharge/discharge processes driven by the occurrence of the fracture/fault systems caused by earthquakes causes an enhancement of permeability by widening fractures [54]. Therefore, areas with low earthquake intensity reveal the low recharge/discharge processes versus areas of high seismicity. The magnitudes of these earthquakes range from 0.1 to 3.9, according to the record of the National Research Institute of Astronomy and Geophysics (NRIAG). The seismic activity density map of the study area ranges between zero and 1251 ( Figure 12b). Thus, the seismic instability map was categorized into 3 rank classifications, ranging from 1 (low density) to 3 (high density). The areas of high intensity cover about 1.36% of the study area and promote groundwater infiltration characteristics (Table 1).

Groundwater Prospect Map
To understand where groundwater recharge is, several evidential maps were compiled, as water infiltration to the aquifers is controlled by geologic and hydrologic factors and landuse/cover factors. The ten (10) evidential maps are geology/geomorphology, topography, slope, depressions/sinks, stream-networks, runoff, radar backscattered, lineaments, rainfall data, and earthquakes. Figure 13 describes the GPZ map that is quantitatively categorized into several zones of analogous sorts. The delineated recharge sectors are termed very low, low, moderate, high, and very high, involving areas of about 6.56, 22.62, 30.75, 29.71, and 10.34%, respectively ( Figure 13a; Table 4). The spatial variations of the very high to high zones cover the downstream areas and fit areas of high wadi deposits on the main streams. This area consists of porous, thick sediments and a low topography, is highly fractured, and contains depressions with a flat to gentle slope, promoting the increase of water infiltration and recharging beneath aquifers that could reveal the groundwater prospectively. Very low 110-487 6.56% The twenty-nine well locations that were overlaid on the GPZs to assess the predicted model indicated that about 13, 12, and 4 wells are consistent with the very high to high, moderate, and low potential zones, respectively (Table 4). About 86.17% of the wells are matched with the very good to moderate potential zones. The well data located in the eastern side belong to the fractured Precambrian aquifer.
The InSAR coherence map reveals the change in the land-cover along wadis as a result of the runoff of active streams during rainy storms. Landsat-8 reveals agricultural activities. The InSAR CCD image (Figure 13c,g,e) provides evidence of water accumulation, as indicated by the Landsat image in panel (d), following the rainfall storm, which is consistent with the active stream marked in yellow. Most of these active streams are most likely holding much water in some ponds (Figure 13c,d) and are consistent with areas of low topography (Figure 13h,i). The basement areas that consist of igneous and metamorphic rocks, particularly in the eastern section of the area, are illustrated by absent vegetation with a low porosity and high topography, with a steep slope that aids the demarcation of poor recharging zones. Remarkably, the evidence of surface water is also detected in the Landsat image (Figure 13j) acquired on 21 January 2010, following a rainy storm (17-18 January 2010), which also matched the rainfall data in Figure 11c,d.

Discussion
Based on the analysis of the optical/radar remotely sensed data, the study area is comprised of crystalline rocks to the east and Phanerozoic sorts to west, which represent the main catchment, and are predominantly affected by NW-SE and E-W lineaments of the fracture/fault systems. Such fault systems could have controlled the groundwater movements and infiltration capacity [50]. Wadis east of Kom Ombo probably experienced intense past hydrological processes, revealed in the denudation of the Cretaceous/Tertiary succession [9,25]. Radar data support that these were affected by E-W lineaments. Abu Subeirah is captured by the NW-SE trend of W. Kharit (Figure 9e-h), due to structural activities in the Kom Ombo area.
Crucially, the present defunct wadi east of Kom Ombo has been an active stream during the past pluvial conditions, like the rest of Egypt's wadis [24,57,58]. Such conditions were likely predominant during the Late Tertiary period or later [9,59] and experienced structural activities related to Red Sea tectonics. This is judged by the existence of paleochannels and the fact that the position of the Equator at that time was more north (most likely at the present-day latitudes of Sudan and Chad) than its current latitude. Such a geometry would have allowed heavy precipitation over Egypt [60] and aided in the removal of the wide part of the Upper Cretaceous/Lower Tertiary succession to the south, as noticed in the Kom Ombo area. Importantly, the exposed Nubian Sandstone covers a wide area of the basin, and such deposits promote surface water infiltration to the groundwater aquifers and are more suitable for storing groundwater [61].
Today, the discontinuous rainfall showers on the Red Sea highlands, as the main catchment, yield surface water flow and water accumulation in ponds, as depicted by the Landsat and InSAR CCD images that allow for the recharging of groundwater aquifers. The Red Sea highlands are the most elevated section of the study area and are comprised of crystalline rocks. Such runoff recharges the aquifers underneath, as evidenced by the existence of surface water during storms.
Groundwater data analysis [56] revealed that the water wells to the east (Figure 13a) that annually recharge, mainly from rainfall, are consistent with fractured and weathered Precambrian aquifers, and their depths range from 3.8 to 31.5 m. The water wells are consistent with zones of high-very high groundwater potential to the west, and are most likely consistent with Quaternary aquifers comprised of sand and gravel, with a thickness from 40 to 120 m and depths ranging from 2 to 8.4 m. This aquifer recharges mainly from the Nile and the surface runoff resulting from the precipitation, as the hydraulic activity of the middle part of the present study is~80 m/day [62].
Based on our analysis of Landsat data, the majority of the downstream area is cultivated using water resources from surface-and groundwater. The existence of the highly vegetated areas led to the existence of the surface-and groundwater and appropriate soil properties that allow for infiltration capability [63]. The areas of agricultural activities are considered the promising areas for groundwater recharge capacity [64], but bare areas of impermeable rocks are recognized as having low recharge capacity [6].
The InSAR CCD image revealed that the highly changed areas are consistent with the CCD image, with low values (~0) that matched the vegetated areas and active streams (Figure 14a,b). The signature of the water is in cyan color in Landsat images (Figure 14c-e). Applying the vegetation index NDVI = (NIR − R)/(NIR + R) to the Landsat data, a series of years (1986, 2000, and 2020) revealed major changes in vegetation, covering areas of 86, 159, and 342 km 2 , respectively (Figure 14c-h). From 1986 to 2000 to 2020, land reclamation increased eastward (Figure 14). The output of the NDVI values extends from −1 to +1; the areas of maximum vegetation land-cover are close to +1. The expansion of agriculture is very extensive, and the presence of the surface water of the Nile and groundwater allowed for the cultivation of the majority of the downstream area.

Conclusions
The integration of ALOS/PALSAR, Landsat-8, SRTM DEM, and TRMM data successfully characterized the geologic, geomorphic, structural, and hydrologic settings of the area east of Kom Ombo and potentially demarcated the groundwater. The present study successfully compiled geologic/geomorphic and hydrologic data for revealing the potential recharge zones. The evidential maps, namely, geology, lineaments, stream-networks, topography, slope, depressions, radar backscattered, rainfall data, seismicity, and morphometric data, were integrated through a GIS-based knowledge-driven approach to reveal GPZs. The resultant map revealed five potential groundwater zones that are very low (6.56%), low (22.62%), moderate (30.75%), high (29.71%), and very high (10.34%). Such zones were verified by matching with agricultural activities depicted from satellite images, well data, and stream activities derived from the InSAR CCD images. Well data revealed that the obtained results will be beneficial for maintaining water supplies, groundwater abstraction, and understanding relationships between the surface water infiltrations in arid regions. The overall outcomes present the fact that combining different data sets derived from satellite sensors as well as geologic, structural, and hydrologic factors is a substantial tool for detecting geologic features and exploring groundwater in arid environments.

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