1. Introduction
Socio-economic activities, sustainability and well-being of human population rely on freshwater supply [
1]. Scarcity of surface water sources, together with a decline in quality and quantity of such resources in arid and semi-arid environments intensified reliance on groundwater for domestic and agricultural use [
2]. There is also the problem of surface water pollution due to anthropogenic activities and lack of wastewater treatment [
3]. In sub-Saharan Africa, the increasing level of freshwater demand due to the ever-growing population necessitates more attention on groundwater sources in many communities [
4,
5]. Ground-based geophysical surveys give precise information on the groundwater potential (GWP) of an area, but the method is expensive and time consuming for a large area. This calls for a cost-effective multi-criteria decision-making approach in conjunction with remote sensing technology and geographic information system (GIS) for identification of potential groundwater resources areas that could pave the way for an effective and sustainable groundwater exploration and exploitation initiatives [
6]. Groundwater occurrence and distribution in a geographical location depends on the anthropogenic, physiographic and climatic factors affecting hydrological conditions of the area. Analyzing relationship and overall influence of such factors on groundwater availability is essential for effective prediction of GWP of an area [
7].
In their quest for a suitable prediction methodology of potential groundwater areas through modeling the effects of physiographic variables on hydrologic parameters to achieve sustainable groundwater management, hydrogeologists found the method of analytic hierarchy process (AHP) as one of the multi-criteria decision analyses models useful for effective groundwater management [
8,
9]. The method provide a means for dealing with complex spatial decision problem by decomposing the spatial decision problem into a hierarchy of sub problems for easy comprehension and subjective evaluation [
10].
Nithya et al. [
11] studied the influence of geology, lineaments, geomorphology, slope, drainage density, soil and rainfall on GWP of Chittar basin India, using GIS based AHP. Results was validated using groundwater yield of 24 open wells. The GWP map classified the areas into very good, good, medium, poor and very poor GWP. Ifediegwu [
12] used integrated RS, GIS and AHP to study the influence of eight thematic maps (geology, rainfall, geomorphology, slope, drainage density, soil, land use land cover (LULC) and lineament density) on GWP in Lafia, Nigeria. Validation was done using data from 50 boreholes. Results classified 19.3, 12.9, 57.8, and 10% of the study area in to good, moderate, poor, and very poor GWP areas. Kamaraj et al. [
13] integrated remote sensing and GIS to determine the effect of 15 parameters (ie: geology, geomorphology, lineament density, topsoil resistivity and thickness, weathered zone resistivity and thickness, first fractured resistivity and thickness, second fractured resistivity and thickness, slope, LULC, drainage density and rainfall) on GWP of Mettur taluk, Salem district and Tamil Nadu, India. Result showed that top soil thickness, pediment, lithology, slope type and rainfall as the predominant factors influencing groundwater potential. Ally et al. [
14] combined seven thematic map layers of lithology, soil types, lineament, magnetic intensity, slope, drainage density and elevation to develop GWP zone map of Mpwapwa District, Tanzania. The resulting GWP zone was validated using area under the curve and overlaying method. It was shown that 19% of the area was classified as very good, 31% as good, 28% as moderate, 22% as poor and 2% as very poor GWP zones.
Researches have also been carried out to compare AHP with other approaches. Razandi et al. [
15] compared the use of AHP, frequency ratio and certainty factor to map GWP of Varamin Plain, Tehran Iran. Receiver operating characteristics was used to determine the accuracy of the GWP maps and results indicated that the frequency ratio out performed both the AHP and certainty factor methods. Kaur et al. [
16] used remote sensing and GIS to compare AHP with Catastrophe theory for demarcation of GWP zones of Panipat district, India. Results showed that the percentage of area covered by each model is almost the same thus the two models can be used for delineation of GWP areas. Shekar and Mathew [
17] combined slope, drainage density, rainfall, geomorphology, LULC, curvature, soil, topographic wetness index, distance from the river, and elevation as thematic layers and used AHP and fuzzy-AHP in GIS environment to assess GWP zones in Murredu river basin, India. The zones were classified as poor, moderate and good GWP areas
AHP has been used, in GIS environment, by many scholars in prediction of GWP areas to explore groundwater resources [
18,
19,
20,
21,
22]. The approach is proved to be a reliable means of facilitating, interpreting, analyzing large volume of diverse dataset which are then modeled through analytic procedure to produce useful information for use by authorities concerned [
23]. For this purpose, the present study employed AHP, using remote sensing and GIS, to integrate hydro-geological and climatic data to assess and map GWP areas which could serve as an exploration/exploitation guide for efficient and sustainable management of groundwater resources in Kano state, Nigeria. Little. The use of AHP in combination with remote sensing and GIS is novel in this part of the world as information is scarce in sub-Saharan Africa. The research will provide policy makers on the best way to sustainably manage the groundwater in the area.
2. Materials and Methods
2.1. Study Area
The study area is in Northern part of Nigeria, West Africa. The area lies between latitudes 10
0 23
′ 40
″ and 12
0 34
′24
″ North, longitudes7
0 41
′ 15
″ and 9
0 21
′ 21
″ East (
Figure 1). The total areal coverage of the state is estimated at 20,131 km
2 with about two-third of the area belonging to crystalline basement rocks and the remaining one-third is sedimentary formation [
24]. Mean annual rainfall is 635 mm in the north and 1000 mm in the south and occurs between May and October [
25]. Kano state is the most populous state in Nigeria with a population density of 764/km2. The 2023 population of the state is about 15,462,200 people with an annual growth rate is 3.2% [
26]. During the 2010 – 2013 period, there was a steady annual decrease of groundwater level and groundwater beneath the floodplains dropped from 9000 MCM to 5000 MCM from 1964 to 1987 in the Chad Formation area of the region [
27]. Surface water is not readily available and this puts more pressure on groundwater. There is an indiscriminate drilling of boreholes in the state and this results in the decline of water table as more boreholes and wells dry up. This underscores the need to have an assessment for the groundwater with a view to sustainably mange it.
2.2. Data Collction
The administrative boundary map, in shape file format, downloaded from DIVA-GIS (
http://www.divagis.com), was employed for the generation of study areal extent. Satellite imagery of LandSat 8 OLI and Shuttle Radar Topographical Mission Digital Elevation Model (SRTM DEM) both having 30 m spatial resolution were obtained from United States Geological Survey site (
www.glovis.usgs.com). Two hundred and fourty five (245) set of boreholes log data spread across the study area, obtained from the archives of Kano Agricultural and Rural Development Authority (KNARDA), were used to generate hydro-geological profile. The boreholes were ranked according to their yield in the order W1 to W245 with W1 having the least yield and W245 the highest. Also, rainfall data of 9 meteorological stations was obtained from Kano State Water Board.
2.3. Data Processing
The methodology adopted in this research includes the following stages as shown in
Figure 2. For the remotely sensed data, the administrative boundary data was processed in ArcGIS environment to mask down Kano State the study area as shown in
Figure 1. Landsat 8 OLI was clipped to the study area and ERDAS 9.0 were applied for the satellite image processing. Land-use-land-cover (LULC) was interpreted using visual interpretation technique and supervised classification of maximum likelihood classification was carried out for final LULC map. For this study slope steepness and stream network information were deduced from SRTM DEM having 30m spatial resolution while the stream network information was consecutively processed to produce flow direction information for the surface drainage, flow accumulation information, stream order information and finally, generating drainage density map for the study area.
Rainfall distribution map of the area under study was produced from the rainfall data of nine meteorological stations across the study area ie (Challawa, Karaye, Gajale, Riruwai, Wudil, Ranka, Joda, Wak and BUK) using inverse distance weighting (IDW) interpolation method. While available information on static water level, soil media, vadose media and aquifer media, obtained from the KNARDA, was used to produce the static water level, soil media, vadose media and aquifer media maps of the study area respectively.
2.4. Analytic Hierarchy Process (AHP)
The AHP is a multi-criteria decision-making process involving experts ranking on the relative importance of some chosen parameters. The steps involved are as follows:
Step 1: This step involves identification of thematic layers. In this work, eight thematic layers of LULC, drainage density, slope, rainfall, static water level, soil media, vadose media and aquifer media were chosen
Step 2: The experts’ opinion was used to generate pairwise comparison matrix (
Table 1) based on their relative importance using Saaty’s scale. Maps of the parameters were classified on a uniform rank of 1 – 5, where a scale of 1 denotes very low, 2 low, 3 implies moderate, 4 represents high and 5 denotes very high GWP [
4,
28].
Step 3: Computations of the normalized weights utilizing the criteria’s geometric mean (
Gm) as shown in Equation 1. The normalized weights are presented in
Table 2
where
Wn = Eigen vector of the matrix
Step 4: The principal eigen value and Consistency Index (CI) were calculated in order to capture the uncertainty in experts’ judgements using following Equation 2
where
= largest eigenvalue from the pairwise comparison matrix,
n = number of classes. The consistency of the pairwise comparison matrix was measured using Consistency Ratio (CR) given by Equation 3
where
RI = random consistency index = 1.41 when
n = 8
[29].
Step 5: The resulting thematic maps of LULC, drainage density, slope, rainfall, static water level, soil media, vadose media and aquifer media were integrated according to Equation 4 using weighted overlay method in ArcGIS 10.5 to generate an index value signifying level of GWP per each cell.
where
GWPI = the GWP index,
Wi = the determined weight of each thematic layer.
Xi = the ranks for the classes within each thematic layer.
Finally, the prediction accuracy of the model was determined and the proposed GWP map was validated with recorded borehole yield data using Pearson correlation coefficient to serve as easy guide for assessment of groundwater availability within the study area (Equation 5).
where
xi is the respective borehole yield value at a particular point,
yi is the GWP index value at that point.
4. Discussion
Globally, the problem of stress on fresh water supply is a well-documented issue. This problem is exacerbated by climate change which manifest itself in drought in some parts of the world. Highly populated areas, such as Kano, become more vulnerable arising from the over-abstraction of water sources due to high water demand from the high population of about 15,462,200 persons [
26]. Surface water pollution and grossly insufficient municipal water supply has led to the demand for groundwater in the state [
2,
27]. The United Nations Sustainable Development Goal No. 6 attempts to address this issue and Nigeria is not on track to achieving this by 2030 as only 29% of the population use safely managed drinking water sources [
35]. Other targets related to SDG No. 6 are also below 50% which is far from the expectations. The country has to fast-track the implementation of policies to achieve this aim.
The geology of Kano State is crystalline basement complex rock in two-third of the basin and the remaining is sedimentary formation. There is a sharp decline in groundwater level which calls for researches towards sustainability of groundwater. AHP was integrated in to RS and GIS to delineate and map GWP areas. Eight thematic layers affecting groundwater recharge potential including LULC, drainage density, slope, rainfall, static water level, soil media, vadose media and aquifer media were chosen in this research. Results of spatial analysis show that better groundwater prospective areas occupied total area of 2165 km
2 (10.70%) and mostly situated in Gurun, Dugol, Zago, Kumbo in the north and Dukku, Doguwa, Riruwai, Dambazau in the southern parts of Kano State (
Figure 11). These are mainly attributed to the nature of the topography, drainage density, rainfall and these areas could be suitable for groundwater development for agricultural utilization. Moderate GWP areas are widely distributed across the study area and covered 12482 km
2 (62%) of the state and could be suitable for groundwater development for domestic use. Meanwhile poor GWP areas covered 5484 km
2 (27.30%) and concentrated mostly in areas having high drainage density and low rainfall at central and southeastern parts of Kano state. These poor GWP areas could not be sufficient for groundwater development for commercial agriculture and hence the need for surface water infrastructure to augment the existing water supply system. For sustainability, efforts should be stepped-up to improve water supply, using surface water sources, in areas close to existing surface waters. This will go a long way in relieving the current pressure on groundwater.
Author Contributions
“Conceptualization, A.A.S. and S.D.; methodology, A.A.S; software, S.I.A.; validation, S.D and A.A; formal analysis, A.A.; investigation, A.A.S..; resources, A.A.; data curation, A.A.S.; writing—original draft preparation, A.A.S.; writing—review and editing, S.D. and S.I.A.; supervision, S.D.; project administration, A.A. and S.D.; funding acquisition, A.A. All authors have read and agreed to the published version of the manuscript.”