Mapping of Groundwater through the Integration of Remote Sensing and Vertical Electrical Sounding in ASALs: A Case Study of Turkana South Sub-County, Kenya

Turkana South Sub-County falls in the arid and semi-arid lands of Kenya, which are characterized temperatures of 20˚C to 41˚C with an average of 30.5˚C and precipitation in the range of 52 mm to 480 mm per year. The area has limited availability of water resources. The area has a land surface of 18,000 km 2 and lies between Longitudes 35˚10'00" and 36˚45'00" East and between Latitudes 1˚0'00" and 3˚0'00" North. The study area faces immense difficulties in trying to meet its water requirements for her rapidly increasing populace, livestock needs and other developmental programmes. Surface water in the area is becoming unobtainable, thus the choice of groundwater ex-ploitation. Indeed with the snowballing demand necessitated by energy, agricultural and livestock production needs in this area, there is an appeal into investigation for groundwater in this greatly remote and extensive area. Thus, the success in the exploration, development and management of groundwater in such a large area calls for such methods that can easily be used to zero down to exploitable targets. Given the expansiveness of the study area, remote sensing (RS) has been used to map lineaments which in turn have been used to deduce faulting affecting the geology of the area. The mapping by RS helped in identifying weathered zones, the fracture systems and fault zones sign of deeper weathered breaks were mapped and it was observed that there was a comparison of the features to the coincidence raster and aspect interpreted lineaments. The earmarked points were later investigated using vertical electrical sounding (VES) to establish their feasibility for groundwater availability and extraction. The investigated VES points in areas of high lineament density revealed decreasing resistivity with increasing depth of investigation, a confirmation of faulting. The data presented curves, type KH, type KK, and type Q. Further investigation was done in points of low lineament density, considered to be least affected by faulting and deep weathering, which turned out to be having increasing resistivity with increasing depth of investigation, an indication of layered lithology. The zones of low lineament density are represented by a type A curve, a type H curve, type HH. There is a good agreement in the results of remote sensing data and VES data models generated where the high lineaments density points concur with areas of good groundwater potential.


Background of the Project
As a factor of growth and development, water is a subject of great standing in many arid and semi-arid lands (Kamal & Hamid, 2011  Sub-County characteristically is defined by warm and hot climate (Sultan, 1970).
Recent study by Opiyo (2014) show that extreme drought scenarios in Turkana have gone up over the last 63 years, with 28.5 percent occurring between 1950 and 1970, up to 47.9 percent in the last two decades from 1990 to 2012. Indeed a report on Turkana county gives a general view of temperatures of 20˚C to 41˚C with an average of 30.5˚C (Government of Kenya, 2013). The region has two rainfall seasons, where in the months of April to July experiences the long rains (akiporo) and the months of October to November happens the short rains in the range of 52 mm to 480 mm per year and an average of 200 mm (Turkana County Government Report, 2015).
The main drainage in the Sub-County is controlled by Kerio and the Suguta Rivers with their tributaries (e.g. Kalabata, Kasamanang and Nasaken) which Sultan (1970) notes that these are sandy rivers that contain water at depth and D. M. Nyaberi et al. these supplies are of considerable importance in this arid region. Groundwater forms one important source of water in the Sub-County through boreholes and shallow wells although the outcome of sited and drilled boreholes is that 50 percent ends up dry. Groundwater drilling is very expensive and thus its exploration calls for the putting into practice of actual approaches of precision and that saves both money and time (Jawad & Yahya, 2013). Therefore, without proper survey the results are uneven supply of water and uneconomical programmes of water provision through boreholes of which most turnouts dry upon drilling. Jawad and Yahya (2013) notes the usability of applying remotely sensed data for establishing the most favourable locations of underground water, thus providing the best initial approach of pinpointing sites which can be pursued later by other ground methods.
Most boreholes and shallow wells in the study area are drilled proximate to rivers and stream channels, a demonstration that surface hydrology acts as a guide to most groundwater surveys rather than the pursuit of hydrogeology and D. M. Nyaberi et al. Journal of Geoscience and Environment Protection hydrogeophysics. The siting of boreholes by means of the geophysical methods only occasionally is limiting and therefore the success of the procedure calls for incorporation of several relevant survey techniques. Indeed digital elevation models (DEMs) and geological maps have been used by Mogaji, Aboyeji & Omosuyi (2011), to accomplish mapping and analysis of lineaments in the Basement Complex area of Ondo State, Nigeria. GIS's application of digital elevation model helps to improve the interpretation (Maina & Tudunwada, 2017) of lineaments in analysis of remote sensing data. In the research the application of a GIS framework by overlaying the satellite imagery with the digital elevation maps helped in the demarcation of the tile lines that eventually resulted in mapping of the lineaments.
The use of remote sensing based data (Surajit, 2014) founded on spatial layers of groundwater controlling parameters through application GIS have been effective in coming up with the groundwater potential zone maps. Surajit (2014) further remarks that remote sensed data does not stand in place of the in situ information gathering, but rather it helps to minimize the measure of field information accumulation by other techniques for it is proven to assist in targeted measurements, therefore it is still vital to affirm the correctness of remotely sensed information and unquestionably remote sensing collaboration with other techniques. Remote sensing methods have played an unlimited role in groundwater searching in hard rock areas especially on lineament identification (Anbazhagan & Manivel, 2010) and have been applied successfully to delineate areas with groundwater occurrences.
In this study, an endeavour was made towards the delineation and classification of probable underground water potential areas by using an integration of remote sensing and VES data. The correlation of faulting to changes in the course of the rivers (Mason & Gibson, 1957), common feature of faulting as evident with the presence of Trellis drainage (Fairburn & Matheson, 1970) and their association to groundwater is on large scale investigated. The research maps the groundwater potential of the Sub-County whose outcome will eventually act as a guide to siting productive boreholes and shallow wells contrary to the prevailing situation. The current approach for mapping used is able to assist in prioritizing potential regions thus offering backing for site-specific field geophysical examinations towards borehole siting.

Hydrogeologic Background of the Study Area
The Sub-County is an arid and semi-arid land and thus availability and access to water remains the biggest problem facing the people of Sub-County. The key water sources supplying the rural parts in Sub-County presently are unprotected hand dug wells, streams and boreholes. The hydrogeological situation of the region is rather difficult, with the geology of this part of the world being predominantly basement metamorphic rocks and igneous rocks, all of which naturally forms poor aquifers. The account of the geology of the area is given as Miocene basalts (Dodson, 1963) nephelinites (Champion, 1935) and amphibolites and the  McCall, 1964). The existence of volcanic rocks' source is taunted to be from a number of central rift volcanoes (Sultan, 1970). Gwynne (1969) has given an account of the occurrence of recent sedimentation in river channels which are associated with groundwater potential.
Equally, fracture and fault zones offer the best areas for groundwater movement, recharge and storage. Equally some areas are covered by sedimentary rocks both consolidated and unconsolidated including clastic sediments, limestone and other carbonate rocks.

Remote Sensing
Remote sensing denotes, the attaining of info concerning things or spaces on the Earth's surface deprived of the direct contact with the object or area, where many sensing devices record information related to an object by measuring ref- lection and radiation of electromagnetic energy transmitted from the object's surfaces (Shefali, 2013). The data used is from sentinel 1 (Zakaria et al., 2017), is more effective in restitution of lineaments as compared to Landsat-8 and ASTER. The sentinel 1's operation is in four exclusive acquisition modes: 1) and 4) Wave (WV). The one of interest is the operation in single polarization (HH or VV) and dual polarization (HH + HV or VV + VH). This is normally executed by one transmit cable (switchable to H or V) and two parallel receive cables for H and V polarization and on synthesis generated lineaments which are key in analysis of the terrain. The secondary data from the primary conflict-free mode over land, IW was used. Sentinel-1 Wave mode is similar to ERS and Evnisat wave mode imaging but with improved spatial resolution, larger vignettes and a "leap frog" acquisition pattern.
The sentinel series observes the earth at high resolution and therefore is very useful in groundwater management and other applications and information services. In this research, the assembly of a false-colour composite image (Sentinel 1a cross-polarised (VH)) image for the study area by combination of selected spectral band ratios and maps was used to develop Sentinel 1a VV Polarised image. The decision to use VH, is that they show a good contrast and enough detail permit hydrographic network mapping where minimum vegetation is available (Saint-Jean & Singhroy, 2000), whereas the VV polarization permits a better positioning of alignments in form of faults, scarps, banks and therefore vital in mapping roughness variations.

Vertical Electrical Sounding Method
The The theoretical assumption is that the medium penetrated by the current in the sub-surface of the earth is being homogeneous and isotropic and thus the resistivity p, and the potentials V M and V N measure happening at M and N, respectively. The total current I flows away from or toward each electrode across the surface of a half sphere with area ( ) 2 1 4 2 r π and with the distances between the electrodes given by r 11 , r 12 , etc., and V = 0 infinitely far from the current source, the potentials at M and N are given by (Rhett, 2001): Combination of the two equations gives a complete resultant in field measurement where the outcome equation is given as: subsequently, the resistivity of the homogeneous subsurface is resolved by: and where K is referred to as the geometric factor as applied to the Schlumberger configuration and is given from the relationship: otherwise normally, the geology is heterogeneous thus its resistivity is presented as apparent resistivity: and in cases where r 11 , r 12 , r 21 , r 22 values are given in metres, and ∆v and I are measured in millivolts and milliamperes respectively, then the pa would be in ohm-meters.
The VES data was collected in points of interest in the study area pursuant to

Results
The remote sensing data was processed into giving drainage lines and lineaments extraction. Lineaments were extracted automatically by considering two fundamental calculations: foremost, was the usage of a filter in the detection of   The analyzed data presented in Figure 6 and Figure 7 paints a picture of two separate scenarios in the areas investigated. The profiles in Figure 6 are VES results of points of low-density lineaments while the profiles in Figure 7 are VES results of points of high-density lineament.

Discussion
The The VES curve for point AA is a type KK which means that there are two combined K-type graphs with apparent resistivities ρ 1 < ρ 2 > ρ 3 < ρ 4 . This presents a four-layer model in which the information indicates geoeletrical layers of decreasing values that indicate more potential for groundwater potential as one goes deeper. The VES curve for point AC is a type KH graph with apparent resistivities ρ 1 < ρ 2 > ρ 3 > ρ 4 which is also indicative of a four-layer model. It portends two aquifers where there is a shallower aquifer at about 10 metres deep which is associated with the sediments deposited in the area and a deeper one at about 40 to 80 metres associated with a concealed fault system in the area. The VES at point AD represents a type KK curve which means that there are two combined K-type graphs with apparent resistivities ρ 1 < ρ 2 > ρ 3 < ρ 4 . This  The layered resistivity model for point AB represents a type A curve which is a three-layer model with the apparent resistivities ρ 1 < ρ 2 < ρ 3 . This shows an increasing resistivity with depth and that signifies a lithological setting which obeys the Law of Superposition where in its interpretation given the way of overlaying of geology the outcome is reduced porosity and permeability and therefore reduced electrical conductivity in the rocks with depth. This curve affords an interstitial aquifer about 20 metres deep which might be a boundary between two geological formations overlaying one another. These zones at times de- This is also a three-layer model with apparent resistivities ρ 1 < ρ 2 > ρ 3 . This presents a shallow aquifer up to 15 metres associated with sediments but otherwise the deeper regions are not viable. The presented resistivity also shows a continuous lithology obeying law of superposition. The curve type for layered resistivity model for point AF is type HH, with apparent resistivities ρ 1 > ρ 2 < ρ 3 , ρ 1 < ρ 2 < ρ 3 , representing a three-layer model. This Remote sensing data is epitome for finding an overall overview of the morphological characteristics of a large research area at the commencement of the exploratory process, thus defining a smaller area as priority targets for more localize follow-up studies. This has effectively been applied to map lineaments in Turkana South Sub-County resulting in targeted points which have been confirmed through VES to be of potential for groundwater exploitation. Electrical resistivity through VES mode of deployment was handy in determining the resistivity of the geological media for the purpose of determining groundwater prospective in fractured unconfined aquifers in the predetermined points through remote sensing.

Conclusion
The remote sensing method has been successfully used to map regions of high D. M. Nyaberi et al.
and low-density lineaments in the study area. The regions of high-density lineaments were interpreted to refer to areas affected by faulting and deep weathering. In consideration of the interpretation of remote sensing data, the researcher concludes that the area is affected by different fault elements trending NE-SW, N-S and minor in the E-W reflecting the orientation of different lithologies and as an indication of past tectonic activities. Further use of vertical electrical sounding resistivity data resulted in the conclusion that the occurrence of groundwater in the study area is considered to exist either in zones of weaknesses necessitated by faulting and subsequent deep weathering, or in the interstitial openings in the beddings.
Thus it is evident that with integration of remote sensing and VES data, there is a better understanding of groundwater availability of the area. Although with the studies of using an integration of various method of survey is commended, still holding is the fact that groundwater potential of an area involves establishing its quantity and quality. Thus this is achieved by undertaking drilling and well test pumping which helps determine the yield and the subsequent quality analysis of the water to establish its portability.