Next Article in Journal
Comparison of Two Types of Survey Designs for Acoustic Estimates of Fish Resources in the Three Gorges Reservoir, China
Previous Article in Journal
Research on Financial Support, Technological Improvement and Marine Economic Development for China’s Coastal Regions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Hybrid Spatial–Analytical Network Process Model for Groundwater Inventory in a Semi-Arid Hard Rock Aquifer System—A Case Study

1
Centre for Advanced Research in Environment, School of Civil Engineering, SASTRA Deemed University, Thanjavur 613401, India
2
Department of Geology, Alagappa Government Arts College, Karaikudi 630003, India
3
Peter the Great St. Petersburg Polytechnic University, Saint Petersburg 195251, Russia
4
Division of Research and Innovation, Uttaranchal University, Dehradun 248007, India
*
Author to whom correspondence should be addressed.
Water 2022, 14(17), 2743; https://doi.org/10.3390/w14172743
Submission received: 14 June 2022 / Revised: 2 August 2022 / Accepted: 3 August 2022 / Published: 2 September 2022
(This article belongs to the Section Hydrogeology)

Abstract

:
Growing agricultural, industrial, and residential needs have increased the demand for groundwater resources. Targeting groundwater has become a challenging endeavour because of the complex interplay between varying climatic, geological, hydrological, and physiographic elements. This study proposes a hybrid RS, GIS, and ANP method to delineate groundwater zones. The resource was evaluated using seven surface hydrological and six subsurface aquifer parameters. The analytic network process model was used to determine the global priority vectors of each subclass. Surface and subsurface groundwater potential maps were created by assigning the resulting weights and spatially integrating them. Later, an integrated potential map was created by combining them. The validation of the obtained results using water level data demonstrates that the integrated map accurately predicted the zones. The area under study has 172.94 km2 of good groundwater potential. An area of 393.01 km2 is classified as having a moderate potential, and an area of 410 km2 is classified as having low potential. These findings will be beneficial to regional policymaking and long-term groundwater management. The results show that an integrated approach using ANP can better determine the groundwater potential zones in semi-arid zones.

1. Introduction

The global water cycle has undergone significant changes over the last century (Milly et al., 2005 [1]). Urbanization and population growth have increased domestic, agricultural, and industrial water needs (Okello et al., 2015 [2]). The per capita actual and virtual water consumption has also significantly increased (Hoekstra and Chapagain, 2007 [3]). The present annual global water demand is around 4600 km3, and it is projected to rise between 5500 and 6000 km3 by 2050 (Wada et al., 2014 [4]; Boretti and Rosa, 2019 [5]). Increasing demand and poor water quality have made groundwater the only reliable source in many parts of the world. Consequently, extreme high groundwater extraction rates are witnessed in parts of America, Asia, Europe, and the Middle East. Especially in the Russian Federation, the water supply of more than 50% of the territory and population is based only on groundwater, but now some territorial entities including the Murmansk, Kurgan, Omsk, Novgorod, and Yaroslavl regions, certain areas of the Arkhangelsk, Rostov, and Tyumen regions, the Republic of Kalmykia and the Stavropol Territory are experiencing significant water deficits (Pykhtin, 2019 [6]). Globally 30% of the world’s largest groundwater systems are in distress.
In India, too, one of the most urgent sustainability challenges is the rapid depletion of groundwater. Several aquifers are experiencing water stress (Gleeson et al., 2012 [7]). For the past 20 years, India has been the world’s most rapidly depleting region in terrestrial water storage (TWS), with a rate of 3–4 centimetres per year (World Bank, 2012 [8]). As reported by the Ministry of Jal Shakti, the annual per capita water availability has decreased from 1816 m3 in 2001 to 1545 m3 in 2001 and is projected to drop to 1367 cubic metres in 2031 (MoWR, 2021 [9]). Unsustainable irrigation pumping and changes in land use have depleted the groundwater (Srinivasan et al., 2015 [10]; Penny et al., 2018 [11]). Environment India predicts that by 2050, six states will be water scarce, and four will be water stressed (Grasso, 2021 [12]). In 2000, nonrenewable groundwater extraction for irrigation contributed approximately 68 km3, the highest contribution in the world (Wada et al., 2012 [13]). The water table in the northwest region of India has dropped below 20 m and is continuing to fall at a rate of more than 1 m per year (MacDonald et al., 2016 [14]; Mishra et al., 2018 [15]).
Due to overexploitation, unsaturated shallow aquifer zones are expanding horizontally and vertically (Luo et al., 2019 [16]). Furthermore, future climate change will likely impact shallow unconfined aquifers (Russo and Lall, 2017 [17]). Therefore, it is essential to precisely target the deep groundwater potential zones to exploit the shallow groundwater rationally. However, due to diverse geological formations, intricate tectonic structures, climatic variations, and various hydrological conditions, groundwater behaviour is complex in deeper hard rock zones. Therefore, a large quantity of multidisciplinary data and a well-defined methodology are required to delineate the potential zone and model the groundwater.
Researchers have developed various direct and indirect methods to explore groundwater over the last few decades. The surface methods include geophysics (MacDonald et al., 2005 [18]; Porsani et al., 2005 [19]; Anomohanran et al., 2017 [20]), electromagnetic methods (Cook et al., 1989 [21]; Selvam and Sivasubramanian, 2012 [22]), and geology (Solomon and Quiel, 2006 [23]; Kudamnya et al., 2019 [24]). However, the investigations are costly, take considerable time, and require a skilled workforce.
Indirect methods, such as logistic regression (Ozdemir, 2011 [25]; Pourtaghi and Pourghasemi, 2014 [26]), decision trees (Chenini, 2010 [27]), artificial neural network models (Corsini et al., 2009 [28]), Shannon’s entropy (Naghibi et al., 2015 [29]), and machine learning techniques such as random forest and maximum entropy (Srivastava et al., 2017 [30]; Elbeltagi et al., 2021 [31]), and so on, have recently gained attention in demarcating resources. However, geographical information systems (GIS) and remote sensing are also powerful and cost-effective in demarcating the resources. Satellite data provide baseline information on groundwater influencing parameters, namely drainage, slope, land use, land cover, soil, geomorphology, lineament and lithology (Saraf and Choudhury, 1998 [32]; Thakur and Raghuwanshi, 2008 [33]; Dar et al., 2010 [34]; Oikonomidis et al., 2015 [35]; Magesh et al., 2012 [36]; Saravanan et al., 2020 [37]). GIS offers a structured framework for storing and integrating spatial and non-spatial datasets (Gumma and Pavelic, 2013 [38]; Senanayake et al., 2012 [39]), but it is important to understand that the parameters used to determine groundwater potential come from various sources. Further, the degree of influence of the parameters will also change over time and space. Such inappropriate geographic data integration from diverse sources might lead to erroneous spatial predictions. However, the analysis can minimise the unprecedented uncertainties by combining geostatistical techniques.
This study attempted to overcome this challenge by combining geographical data with a multicriteria decision-making (MCDM) tool to determine the groundwater potential zone. Multicriteria decision analysis (MCDA) offers a systematic way of evaluating multiple parameters (Kumar et al., 2014 [40]; Rahmati et al., 2015 [41]; Razandi et al., 2015 [42]; Jennifer and Jha, 2017 [43]). The easiest and most comprehensive MCDA technique is the simple additive weighting technique (SAW). Based on significance, weights are assigned qualitatively to each parameter (Kamaraju et al., 1996 [44]; Shaban et al., 2006 [45]). Other statistical methods commonly used to prospect groundwater include frequency ratio (Manap et al., 2013 [46]; Naghibi et al., 2015 [29]; Trabelsi et al., 2019 [47]), logistic regression (Pourtaghi and Pourghasemi, 2014 [26] and weights of evidence (Ozdemir, 2011 [25]). Subsequently, Saaty (1980) [48] proposed the analytical hierarchy approach (AHP). This semi-quantitative top-down decision model is used to predict weights based on the ideal relationship between the parameters through a pairwise comparison. However, all the parameters in the physical world are extremely interconnected. Thus the analytic network process (ANP), an AHP inheritor developed by Saaty (1980) [48], makes it possible to establish the complex interrelationships through a network model (Jenifer and Jha, 2017 [43]; Rajani et al., 2018 [49]; Siva et al., 2017 [50]; Sujatha, 2020 [51]; Gebru et al., 2020 [52]). Since the non–linear structure is much closer to reality, it produces better results. Although the model has extensive application in geo environments (Meade and Sarkis, 1998 [53], Sujatha and Sridhar, 2017 [54]; Sujatha, 2020 [51]), its use in groundwater studies is limited.
Tiruttani is a drought-prone zone faced with extreme water scarcity, which widely affects agricultural and drinking water uses. Besides in a few limited areas, the region’s surface water supply is inadequate to meet local requirements. Therefore, this study aims to create a groundwater prospecting model utilising the ANP approach. Furthermore, groundwater research frequently includes surface factors (Saraf et al., 2004 [55]), whereas the subsurface parameters are least investigated. In this study, an experiment was undertaken to identify prospective groundwater zones by ANP modelling using both surface and subsurface factors and perform a validation by comparing the results with field data.

2. Description of the Study Area

The areas between northern latitudes 12°56′30″ N and 13°18′30″ N and longitudes 79°30′30″ E and 79°47′ E in and around the Tiruttani block of Tamil Nadu were studied (Figure 1). The study comprises an area of 975 sq. km with typical massive rock shoots dissected by fractures. The area’s annual rainfall varies from 950 mm to 1150 mm under the influence of both northeast and southwest monsoons. The main rivers that flow within the study region are Nagari and Palar. Groundwater is found in the weathered mantle under phreatic conditions and in the fissured and fractured zones at deeper depths under semi-confined conditions. The weathered zone is between 2 and 12 m thick.

3. Material and Methods

In the present study, groundwater potential assessment in the Tituttani hard rock region was performed using seven surface and six subsurface parameters. The proposed method consists of four steps (Figure 2): (1) surface parameter-based groundwater prospecting, (2) subsurface parameter-based groundwater prospecting, (3) an integrated new method that considers both parameters, and (4) characterizing the distribution of groundwater potential and assessing the model efficiency.

3.1. Constructing Thematic Layers

The surface parameters used to identify the spatial distribution of various groundwater potential zones were extracted from multiple analogue and digital datasets such as topographical maps, satellite images and digital elevation models. The lithology map was derived from Tamil Nadu’s Geological and Mineral map published by the Geological Survey of India. After geo-rectifying and projecting, it was digitised and converted into a thematic layer using GIS software. The IRS 1D–LISS III geocoded (P142 R51, P143 R51 dated 15 December 2015) digital satellite data was processed using the supervised classification maximum likelihood method and the geomorphic units were mapped. Likewise, land use and land cover units were interpreted by adopting the NRSC level II classification system (NRSC manual, 2012 [56]). A digital elevation model (DEM) with a spatial resolution of 90 m × 90 m and an elevation accuracy of 3 m to 5 m was developed using Shuttle Radar Topographic Mission (SRTM) data. The spatial distribution of lineaments was extracted from satellite images and further updated using the DEM, and shaded relief outputs were produced from the SRTM data. Subsequently, the lineament density map was then generated by the inverse distance weighted (IDW) method. A drainage map was made by scanning, rectifying, and digitising the Survey of India topographical sheets (57O/8, 57O/11, 57O/12, 57O/16, 57P/9) of 1:50,000 scale, and a map on drainage density was then prepared as earlier. The soil map was reproduced from the soil resource map published by the soil survey and land use organisation (SSLU) of the Tamil Nadu agricultural department. The slope gradient was derived from the SRTM DEM. Ground truth verification a was performed using a GPS survey. The interpreted geomorphology, land use, land cover and major lineaments were checked in the field.
A database of subsurface parameters influencing the occurrence and movement of groundwater was created. Geophysical resistivity surveys were performed at ten locations, and the aquifer thickness map was then prepared. A total of five bore wells, five dug cum bore wells, and 15 dug wells were selected for investigation. GPS was used to geotag the locations of the wells mentioned above, and their water levels were recorded. A pumping test was conducted, and parameters such as permeability, transmissivity, specific capacity and storage coefficient were determined in the field. Later, they were spatially plotted, interpolated and classified using GIS software.
The analytical network process (ANP) model was employed in the study for delineating the groundwater potential zone (Figure 2). The ANP framework was constructed for surface and subsurface parameters separately using super decision software. Global priority vectors indicating the influence of each parameter compared with others of being a groundwater potential parameter were derived. Derived vector values were dully allocated and transformed into raster layers. By spatially adding the respective data sets, surface and sub–potential surface maps were created. In addition, an integrated groundwater potential map was also generated by combining the resulting outputs. The water level data were obtained at random intervals and classified into shallow, moderate, and deep based on their depth to validate the outcomes. Correlation coefficients were obtained by correlating them spatially with the potential surface, subsurface, and integrated groundwater maps.

3.2. Analytical Network Process Model

ANP is a multicriteria decision model that supersets the well-known Saaty (1980) [48] developed analytic hierarchy process (AHP) model. ANP provides a broad framework for dealing with choices without assuming anything. ANP reduces a multicriteria complex problem to a non-linear network (Saaty, 1980 [48]; Sujatha and Sridhar, 2017 [54]; Sujatha, 2020 [51]). More importantly, the interdependence of network elements allows better modelling of complex problems since most real-life situations are non-linear. Thus, they are closer to reality and, as a result, yield more accurate results. This study developed two ANP models to target groundwater: one based on surface parameters and the other based on subsurface parameters.
There are four stages to the ANP modeling process.
Step-1. This step involves framing the network structure accordingly. The goal is to identify a potential groundwater zone. The thematic layers form the criteria group, while the categories in each theme include the sub-criteria group. The potential zones, namely high, moderate and poor, comprise the alternatives group.
Step-2. The pairwise comparison between the elements in each group are correlated with all other elements in the network by Saaty’s 1–9 scale (Table 1). Subsequently, local priority vectors are derived.
If both the sets of comparisons then prove inconsistent, a review is necessary. As a result, verifying that they are consistent is required before using comparison matrices for evaluation and analysis. The most well-known consistency index for the pairwise comparison matrices in the ANP is the consistency ratio by Saaty. The consistency index (CI) of the derived weights of each set of pairwise comparisons is calculated by:
CI = (λmax − n)/(n − 1)
where n is the size of the pairwise matrix, an if CI is less than 0.10, it is assumed that the judgments are consistent (Saaty, 1980 [48]). Similarly, the consistency of the entire pairwise comparison matrix is checked using the consistency ratio (CR):
CR = CI/RI
where RI is a random consistency index computed for different n values is given in Table 2.
For perfect consistency, the consistency measure λmax equals n; therefore, the consistency index CI equates to zero, resulting in a CR of zero. However, this is not the case on many occasions, and the CR usually has a value greater than zero. Saaty recommends that CR be 0.1 or less, and the expert should revise the judgments if CR is significantly greater than 0.1.
Step-3. To obtain the overall priorities of the elements, a supermatrix is formed. The super matrix serves as a unifying framework, representing the overall relative dominance of one element over another in the network.
The diagram labeled “Figure 3” provides a description of the overall structure of the super matrix, where Wij is a block matrix consisting of priority weight vectors (w) of the influence of the elements in the ith cluster concerning the jth cluster, eNn signifies the nth element in the Nth cluster, and CN signifies the Nth cluster. In a situation where the ith cluster has no effect on the ith cluster, Wij is set to zero (a case of inner dependence). The super matrix created by completing this phase is referred to as the “initial supermatrix.”
There are three different kinds of super matrices: (i) the limit super matrix, (ii) the weighted super matrix, and (iii) the unweighted super matrix. No cluster comparisons were made since the current ANP model has just two clusters. Therefore, the weighted and unweighted super matrix are identical. As a result, the values of the other two super matrices remain unknown.
Step-4. By computation, global priority vectors representing a numeric value between 0 and 1 are derived, which indicate the relative importance of each element.

4. Database on Factors Influencing Groundwater Potential

4.1. Surface Parameters

The earth’s surface morphology and its characteristics determine the rate of infiltration and indirectly influence the groundwater potential of the region. The surface groundwater influencing parameters such as drainage, slope, land use and land cover, soil, geomorphology, lineament and lithology were considered (Abdalla, 2012) [57].

4.1.1. Lithology

The six major lithological units in the study area are shown in Figure 4a. The largest are the hornblende fissile-biotite gneiss and epidote-hornblende gneiss (68.25 percent). Within the gneiss, basic dykes with an EW orientation are noted. Granite patches are observed in the central, fluvial deposits in the southern, and clay in the eastern parts of the study area. Deep and moderate weathering is observed in the gneissic and charnockitic formations, respectively. The thickness of the joint/fractured zone is considerable in gneissic formations and comparatively less in granite and charnockite. Gneiss was assigned high potential in groundwater, charnokite and sandstone as moderate potential, and granite with low potential (Krishnamoorthy et al., 2000 [58]). Pairwise comparisons are established in ANP accordingly.

4.1.2. Geomorphology

Geomorphology is a field that deals with groundwater accumulation and movement in a terrain (Kale and Gupta, 2001 [59]). The pediplain geomorphic unit comprises the maximum part of the study area (738.63 sq. km), followed by the alluvial plain (125.43 sq. km). Although denudational hills, floodplain, pediment, and structural hills only form a a small portion (Figure 4b). Pediplains are characterised by a thick weathered zone and are moderately favourable for groundwater. The alluvial plain also signifies groundwater occurrence. Due to their closeness to streams and soil features, flood plains account for 3.45% of the study area and have the greatest potential for groundwater recharge [60,61,62].

4.1.3. Land Use and Land Cover

The potentiality of groundwater and the land surface characteristics are highly interrelated. Land use and land cover govern the degree of infiltration and runoff. Most of the study is dominated by cropland (598.74 sq. km), followed by water bodies (107.93 sq. km), fallow land (33.43 sq. km), and wasteland (22.97 sq. km). Smaller areas are covered by built-up land (31.70 sq. km), land with scrub (63.19 sq. km), land without scrub (33.56 sq. km), and plantation classes (78.69 sq. km) (Figure 4c). The croplands are irrigated and have a higher water holding capacity, allowing time for infiltration and are, thus, a more influential parameter in prospecting (Siva et al., 2017 [50]).

4.1.4. Lineament Density

Lineaments are the rectilinear or barely curvilinear structural components formed due to seismotectonic events. The remotely sensed data lineaments were mapped based on texture, relief, drainage, and tonal linearity. An isoline map was subsequently prepared and graded into high, medium, and low lineament density zones using the inverse distance weighted method (IDW) in ArcGIS (Figure 4d). The high-density region is approximately 460.25 sq. km. The medium-density region occupies 398.23 sq. km, and the low-density region is 117.43 sq. km. Increased secondary porosity enhances groundwater recharge in a highly faulted and fractured region (Yeh et al., 2016 [63]). Thus, the higher the lineament density, the more significant the groundwater potential.

4.1.5. Drainage Density

The drainage pattern is a reflection of the geological and topographical influences of the rocks that underlie it. The groundwater condition of an area can be inferred through its drainage pattern and density (Magesh et al., 2012 [36]). Dendritic, trellis and parallel patterns are prevalent in the study. In addition, rectangular, circular, and radial patterns are also observed. The drainage density was calculated, interpolated spatially, and categorised as high, medium and low drainage density (Figure 4e). The density areas with moderate to medium drainage would have higher surface infiltration and permeability.

4.1.6. Soil

The soil of the study area consists primarily of inceptisol (63%). The eastern portion also includes alfisol and entisol. Additionally, patches of vertisol are seen at four locations (Figure 4f). In large part, the recharge of aquifers is determined by soil infiltration and storage. Entisols are immature soils with little or no horizon and have the original parent material characteristics and, thus, have reasonable infiltration. Vertisol, alfisol, and inceptisol have the lowest infiltration rates due to texture and organic matter (Duiker et al., 2001 [64]).

4.1.7. Slope

The region’s slope map was divided into three classes based on the NRSA classification scheme. Most of the study region falls below the gentle slope category (0 to 1%), and medium slopes (1–5%) and steep slopes (>5%) are found in the north and northwest parts, respectively (Figure 4g). A higher slope suggests a higher potential for erosion, higher surface runoff and less infiltration (Rajaveni et al., 2017 [60]). In contrast, areas with gentle slopes have greater potential for groundwater recharge.

4.2. Subsurface Parameters

Aquifer properties are responsible for the percolation, transport and storage of water in hard rock terrain. A pumping test was performed in twenty-five wells to determine the hydraulic properties such as transmissivity (T), permeability (K), storage coefficient (S), and specific capacity (SC; Banton and Bangoy, 1996 [64]; Sanders, 1998 [65]; Hasan et al., 2020 [66]).

4.2.1. Water Level

Water level provides information on groundwater storativity and is, thus, a significant factor in determining potential regions. The study’s water levels range from 8.1 m below ground level (bgl) to 15.60 m. The levels were subsequently plotted, and isolines were generated using the ArcGIS IDW tool. They were then classified as high (11.31 and 15.60), medium (10.7 to 11.31), and low (10.7 to 11.31), and low (8.01 to 10.07 bgl). The map reveals that all three groups extend around the entire area of study (Figure 5a). The more favourable zone is the shallow water level accompanied by moderate and deep water (Suganthi et al., 2013 [67]).

4.2.2. Aquifer Thickness

For groundwater development, it is critical to have an adequate depth of aquifer since this provides a greater supply of resources. The formation’s resistivity was measured at ten locations up to 100 m in depth. Their actual resistivities were calculated based on the observed apparent resistivity. An isoline aquifer thickness map was developed using the ArcGIS IDW tool, based on the depth at which bedrock was encountered (Figure 5b). The high aquifer thickness zone is located in the central portion (325.58 sq. km). In contrast, the medium zone is located in the southern and northern areas (325.57 sq.km). The study’s western and eastern parts (325.67 sq. km) are predominantly in the low zone. Deeper aquifers have a higher holding capacity and are, thus, more significant in identifying zones with higher potential for groundwater recharge.

4.2.3. Permeability

Permeability indicates the movement of water within the aquifer system. Permeability depends upon the nature of the rock’s porosity and interconnected pores. The study area’s permeability varies from 2.60 to 27.91 m3/day. Based on these values, the study area was interpolated spatially and graded as low permeability zone (2.60 to 10.04 m3/day), medium zone (10.04 to 12.82 m3/day), and high zone (12.80 to 27.91 m3/day). The spatial distribution of the low, medium and high permeability zones is 325.10, 325.50, and 325.16 sq. km (Figure 5c). High permeability values are favourable for a high groundwater potential zone (da Costa et al., 2019 [68]).

4.2.4. Transmissivity

Transmissivity refers to the amount of water transferred horizontally through the aquifer’s formation or its maximum saturation depth. The assessed transmissivity was between 7.79 m2/day and 600.14 m2/day. As before, these values were divided into high transmissivity zones (144.84 to 600.14 m2/day), medium zones (70.51 to 144.84 m2/day), and low transmissivity zones (7.79 to 70.51 m2/day) (Figure 5d). Low and moderate transmissivity zones have a high potential to hold water, whereas high transmissivity implies low groundwater potential; thus, comparisons were made in the ANP model (Hassan et al., 2016 [69]).

4.2.5. Storage Coefficient

The storage coefficient in the study area varies from 0.0001 to 0.0009 and is further graded into high (0.0004 to 0.0009), medium (0.0002 to 0.0004) and low (0.0001 to 0.0002) storage coefficient zones (Figure 5e). A high storage coefficient indicates a good groundwater potential zone, a medium storage coefficient a moderate potential and a low storage coefficient a poor potential zone (Sridharan et al., 1995 [70]).

4.2.6. Specific Capacity

The yield per unit drawdown measures a specific capacity (Q/s). The tests found that the specific capacity ranges between 0.0042 and 0.522. As for the previous parameters, it was spatially divided into low (0.0042 to 0.160), medium (0.160 to 0.256), and high (0.256 to 0.0522) (Figure 5f). High specific capacity is considered a good potential groundwater zone, the moderate specific capacity area is considered medium and low has poor potential. The comparisons were based on the above relationship (Summers, 1972 [71]).

5. Results and Discussion

5.1. Mapping Potential Groundwater Zones

The complex ANP model was developed using super decisions software 2.8.0. Four clusters were primarily developed: goal, criteria, alternative, and potential. The potential groundwater zone was designated as a goal cluster. The parameters affecting the surface groundwater, i.e., drainage, slope, land use and land cover, soil, geomorphology, lineament and lithology, were grouped under the criteria cluster. The subclasses of the aforementioned parameters were grouped within the alternative cluster. A fourth cluster was created, namely potential, with good, moderate and poor nodes to extract the potential groundwater outcomes.
Then concerning the goal, each element (parameter) within the criteria cluster was compared to the other (inner dependency). For instance, concerning the goal (groundwater potential), geomorphology was correlated with drainage density, soil, slope, and within the same cluster. Similarly, their subclass elements in the alternative cluster were compared amongst themselves. For example, in addition to comparing it with its contemporaries such as pediment and alluvial plain, the flood plain of the subclass of geomorphology was also compared with subclasses of other parameters such as lineament density high, and forth and barren land in land use and land cover.
Then, the nodes (elements) of the criteria cluster were compared with all nodes in the alternative cluster (outer dependence). For example, concerning groundwater potential, the geomorphology in the alternative cluster was correlated with its subclass and the subclasses of other parameters. In the next step, each node in the alternative cluster was compared with elements in the potential cluster. For example, a pairwise analysis between the floodplain and high drainage density was established in good potential. Similarly, all subclasses in the cluster of criteria were compared concerning the potential cluster’s good, moderate, and poor elements. Likewise, a reverse comparison (feedback) was also made between the elements in various clusters.
A supermatrix was constructed by establishing 361 pairwise comparisons between nodes in the various clusters, and its resultant priority vectors are displayed in Table 3. The analytical output from ANP may be found in the Raw data column. The Normals column contains output data that has been normalized. The Ideals column contains numerical values representing the relative dominance of each parameter over the other within the network in terms of groundwater potentiality (Feibert et al., 2016 [72]). Groundwater occurrence, movement, and quality are all influenced by geological units. This area’s basement is made up of gneiss and charnokite (Figure 4a). These rocks have medium porosity and moderate groundwater potential zones (Hazell et al., 1992 [73]), earning 9 and 10 in the overall priority vector.
Sandstone has a value of 18 because of its poorly sorted and highly cemented character, making it a low groundwater potential zone (Chilton and Seiler, 2006 [74]). On the other hand, Granites and dykes have a lower porosity (Verma and Patel, 2021 [75]) and are ranked with values of 20 and 25, respectively, as low groundwater potential zones. Geomorphologically, floodplains (rank 1, Table 3) and pediplains (rank 3, Table 3) yield moderate to good groundwater (Thapa et al., 2017 [76]), earning them higher rankings. Pediment ranks 11 due to its moderate to poor groundwater potential. Because they are unfractured rock with limited infiltration, structural hill landforms (rank 24, Table 3) have very low groundwater potential (Ramaiah et al., 2012 [77]). Cropland (rank 2, Table 3) and water bodies (rank 6, Table 3) both allow for a large amount of infiltration and have a strong groundwater potential, and as a result, cropland and water bodies are awarded higher rankings in terms of land use and land cover (Varughese et al. 2012 [78], Rajaveni et al., 2017 [60]).
On the other hand, built-up land and wasteland have a limited infiltration capability (Das, 2017 [79]), earning them a grade of 16 and 21, respectively (Table 3). The slope significantly influences surface water infiltration from the ground. If the slope percentage is relatively low, surface water, mostly made up of precipitation, has a greater opportunity to percolate into the subsurface over a longer period. Runoff is more immediate on steep slopes, resulting in less time for water to remain on the ground surface and a significant loss in groundwater recharge potential (Jaafarzadeh et al., 2021 [80]). Rankings of 7, 14 and 27 were calculated for the gentle, moderate and steep slope classes, respectively (Table 3). As far as drainage density is concerned, a high drainage density stimulates high runoff while a low density promotes infiltration (Rashid et al., 2012 [81], Bagyaraj et al., 2013 [82]), whereas lineament density, the process is reversed, meaning that high density promotes infiltration and low density leads to runoff (Lentswe and Molwalefhe, 2020 [83]). The overall rankings derived from the above parameters also substantiate this (Table 3).
Subsequently, the values that were calculated as ideal were then integrated after being assigned to the appropriate variables (subclasses). The overall cumulative values range between 1.6809 and 6.5366 (Figure 6a). In order to make the results easier to understand, they were reorganized into three distinct categories: high (6.5366–4.8461), moderate (4.8461–3.2275) and low (3.2275–1.6809) potential zones (Figure 6a). High prospect zones covering 210.92 sq. km are located sporadically in the northwestern and western regions of the resulting surface groundwater potential map (Figure 6a). The moderate zones with an area of 341.2 sq. km to the west, surround the high zones mentioned above. The areas with low potential, with 423.84 sq. km, are mainly observed in the southern and eastern parts of the study area (Figure 6a).
Next, six subsurface parameters viz: storage coefficient, specific capacity, transmissivity, permeability, aquifer thickness and water level were grouped below the criteria cluster and their subclasses under the alternative cluster. The deeper the aquifer thickness, the greater the fracture density and storage potential. The overall ranking value confirmed this with a ranking of 1 for deep aquifer thickness (Table 4), 8 for moderate, and 15 for shallow aquifer thickness (Raji and Abdulkadir, 2020 [84]). The water level is a significant component in assessing the groundwater potential since it immediately offers evidence about the storativity of groundwater. A shallow water level suggests a greater potential for groundwater (Chatterjee and Dutta, 2022 [85]). Accordingly, the shallow water level category received score of 2 (Table 4), whereas the moderate and deep water level categories receive scores of 10 and 14, respectively. High permeability results in a high groundwater potential zone (Machiwal et al., 2015 [86]). The overall ranking obtained by high (3), moderate (12), and low permeability (18) categories also suggest the same (Table 4). The higher the transmissivity, the more productive the aquifer (Akhter and Hasan, 2016 [87]). The acquired rankings also confirm this with 6, 7 and 16 for high, moderate and low categories of transmissivity (Table 4). A high storage coefficient corresponds to a good groundwater potential zone (5), a medium storage coefficient area to a medium potential zone (11) and an area with a low storage coefficient to a poor potential zone (17). Similarly, a high specific capacity (Kumar et al., 2017 [88]) is recognised as a good groundwater potential zone (4), moderate specific capacity area as medium (9) and an area possessing low specific capacity as a poor potential zone (13) (Table 4).
Thus, the total rankings produced by various aquifer parameters using the ANP model corroborate the fundamental aquifer characteristics (Nsiah et al., 2018 [89]; Dhakate, 2020 [90]; Hasan et al., 2021 [91]; Masoud et al., 2022 [92]). Subsequently, the derived ideal values (Table 4) were substituted by the respective factors, and by integrating all, subsurface parameters-based groundwater potential maps were derived. The resulting cumulative value for pixels varies from 3.0285 to 0.6042. They were then reclassified as high (3.0285–2.214), moderate (2.214–1.4122), and low (1.4122–0.6041) potential zones (Figure 6b). The high potential zone seems to transact as a linear patch from the middle to the northwestern area and comprises 175.57 sq. km. The moderate potential region also has a similar linear pattern covering 431.16 sq. km. At the same time, the low potential region appears to occupy the northeastern, southeastern, and southern parts of the study area as isolated pockets spanning an area of 369.2 sq. km (Figure 6b).
Based on the surface (Figure 6a) and subsurface potential maps (Figure 6b), it is inferred that most of the area has significant groundwater prospects. In contrast, the Tirrutani block is designated by the Central Ground Water Board (CGWB, 2007 [93]) as overexploited. In addition, the above maps revealed an intricate pattern in the spatial distribution of potential regions. Hence, to contemplate the collective impact, both the unclassified surface and subsurface potential maps were combined spatially, integrated, and then an “integrated groundwater potential zone” map was prepared. The resultant pixel values range from 9.5850 to 2.2311 and were later classified into high (9.5850–7.1337), moderate (7.1337–4.6824), and low (4.6824–2.2311) potential zones (Figure 6c).
From the resultant integrated groundwater potential map, the eastern and southern parts seem to have high potential (Figure 6c). The visual correlation of the high potential zone with the surface (Figure 6a) and subsurface parameter maps (Figure 6b) prepared earlier shows that, spatially, the zone is characterised by pediplain, high lineament density, moderate drainage density, cropland, moderate water level, deep aquifer thickness, high permeability, moderate transmissivity, a moderate to high storage coefficient and a high specific capacity. The central area was revealed as a low potential zone (Figure 6c). All the surface and subsurface parameters in this low-prospect region are less conducive to groundwater (Figure 4a–g and Figure 5a–f)). The spatial assessment shows that, despite the alluvial plain, the least desirable parameters, namely land with and without scrub, plantation and fallow land, low lineament density, moderate to high drainage density, gentle to moderate slope, moderate to deep water level, low to medium permeability, medium to high transmissivity, low to medium storage coefficient, and low to medium-specific capacity are observed in this zone. The eastern and northern parts are categorised as a moderate groundwater potential zone (Figure 6c). The spatial correlation shows moderate potential parameters in this region. Thus, the results obtained corroborate the general opinion offered by many earlier workers (Garewal et al., 2015 [94]; Owolabi et al., 2020 [95]; Singha et al., 2021 [96]).

5.2. Validation of Results

The results of the surface, subsurface and integrated groundwater potential maps were validated using water level data. The water level was collected from 35 locations and based on their depths to ground level, fourteen wells were categorised as shallow (2.40 to 10 m bgl), thirteen as moderate (12 to 18 m bgl), and eight as deep (24 to 90 m bgl). This data was then spatially correlated with the above maps (Figure 6a–c). The frequency ratio method is used as an indicator for validating the potential maps (Ozdemir, 2011 [8]; Naghibi et al., 2015 [12]). The frequency ratio is an approach to bivariate statistics that may be described as FR = (W/TW)/(CP/TP), where W is some pixels of well, i.e., good, in each category, and TW is the total of good pixels in the study area. TP is the total number of pixels in the study area, and CP is the number of pixels in each potential zone category.

5.2.1. Spatial Correlation between Water Level and Surface Parameters-Based Groundwater Potential Zones

The number of pixels comprising each potential zone category was first calculated. This shows that the overall study area consists of 1,070,566 (TP) pixels. From this total, 466,483 (CP) fall in the high potential zone category, 374,608 pixels in the moderate, and 229,475 pixels in the low potential category. Then, by dividing the number of pixels in each category by total pixels, the CP/TP ratio was calculated. For example, the ratio between pixels under the high potential zone category to overall pixels is 0.3636. Subsequently, the W/TW value was calculated by dividing the number of wells in each water level category (W) by the total number of wells (TW). For instance, 13 wells fall into the shallow water level category in the high groundwater potential zone. Thus, the W/TW ratio is 0.3714. Using these derived values, the frequency ratio was calculated. Accordingly, the frequency ratio was calculated for all the water level categories. The results show that the shallow water level falling within the high potential zone category among various correlations has a clear correlation (0.8542). In contrast, all the other categories exhibit only a poor correlation (Table 5).

5.2.2. Spatial Correlation between Water Level and Subsurface Parameters-Based Groundwater Potential Zones

Similarly, the spatial correlation between groundwater level categories and subsurface parameters based on potential groundwater zones shows that eleven out of fourteen shallow wells fall in the high potential zone. Out of 13 moderate water level wells, five fall in the high potential zone, four in the moderate potential zone and four in the low potential zone. Two wells in the deep water level well category fall in the high potential zone, one in moderate and five in low potential zones.
Using frequency ratios, the coincidence between water level categories and subsurface parameter-based groundwater potential was calculated (Table 6). This shows a good correlation between shallow water level and high potential zone (0.820193) and between deep water level and low potential zone (0.78834). In contrast, the other parameters have very poor correlations.

5.2.3. Spatial Correlation between Water Level and Integrated Groundwater Potential Zones

The spatial correlation between water level categories and integrated subsurface parameters shows that 13 shallow wells are in the high potential zone. One is in the moderate potential zone, which is similar to the results obtained from the correlation between water level and surface parameters. However, for the moderate category, out of 13 wells, 12 wells are in the moderate category, and one well is in the deep category. In the deep water level well category, all eight wells fall exactly in the low potential zone category.
The frequency ratio values between water level categories and integrated groundwater potential zone categories (Table 7) showed a very good correlation in all categories. The frequency ratio value between shallow water level and the high potential category is 1.021087. Similarly, the frequency ratio between moderate water level and the moderate potential zone is 0.756871, thus indicating a correlation. Likewise, deep water level wells and low potential zones also exhibited a good correlation (1.247328).
The water level and the surface potential map show a good correlation for the high potential category, whereas the other classes only show a weak correlation. However, the subsurface parameter strongly correlates between shallow water level and high potential zone and between deep water level and low potential zone. The other parameters, however, showed a very weak correlation. Interestingly, the frequency ratio values obtained for the integrated groundwater potential zone show a precise correlation in all classes. Thus, from the above, it can be surmised that the integrated map comprising both surface and subsurface parameters yields better results, followed by subsurface and surface parameter-based potential maps.

6. Conclusions

The current study was undertaken to identify and classify groundwater potential zones in the Tiruttani region using the RS/GIS-based ANP model. Three groundwater potential maps were created using (1) surface parameters, (2) subsurface parameters, and (3) a combined map by integrating them. The validation with ground truth data reveals that the surface parameter-based output accurately predicted the shallow potential zones, whereas the subsurface-based model foresaw the deep zones. On the other hand, the integrated method successfully predicted the shallow, moderate, and deep potential zones. The findings suggest the need of an integrated approach, especially under the current global groundwater depletion scenario.
Considering the integrated map, the good groundwater potential region tends to cover only a small area of 172.94 sq.km (17.72%), indicating that the larger portion of the study are is experiencing severe groundwater scarcity. At the same time, moderate potential zones comprise 393.01 sq. km (40.27%) and the low potential zones are 410 sq. km (42.01%). As per the CGWB report, the water table is drastically depleting in and around the study area. Interestingly, the area receives 1104 mm of mean annual rainfall and has a very good aquifer thickness. By augmenting site-specific recharge mechanisms, the groundwater resource in the Tiruttani area could be enhanced.
Our findings demonstrate that the ANP-based GIS model is a cost-effective and dependable method for delineating groundwater potential zones in hard rock terrain. Though ANP is ideal for establishing complex interrelationships between real-world factors, it is recommended that relevant criteria and influencing factors be brainstormed and selected based on the hydrological conditions of the region. The results can be further improved by using appropriate high resolution satellite images and DEM. Incorporating sizable pump test data derived from various depths and across several thematic sub-classes can yield even better results. In terms of future methodology, developing techniques based on machine learning is an exciting prospect.

Author Contributions

Conceptualization, S.R. and R.R.; methodology, S.R. and R.R.; software, S.R. and R.R.; validation, S.R. and R.R.; formal analysis, S.R., R.R., G.M. and N.I.V.; investigation, S.R. and R.R.; resources, S.R. and R.R.; data curation, S.R. and R.R.; writing—original draft preparation, S.R., R.R., G.M. and N.I.V.; writing—review and editing S.R., R.R., G.M. and N.I.V.; and visualization, S.R. and R.R.; supervision, S.R. and R.R.; project administration, S.R. and R.R.; funding acquisition, G.M. and N.I.V. All authors have read and agreed to the published version of the manuscript.

Funding

The research is partially funded by the Ministry of Science and Higher Education of the Russian Federation under the strategic academic leadership program ‘Priority 2030’ (Agreement 075-15-2021-1333 dated 30 September 2021).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Milly, P.C.; Dunne, K.A.; Vecchia, A.V. Global pattern of trends in streamflow and water availability in a changing climate. Nature 2005, 438, 347–350. [Google Scholar] [CrossRef]
  2. Okello, C.; Tomasello, B.; Greggio, N.; Wambiji, N.; Antonellini, M. Impact of Population Growth and Climate Change on the Freshwater Resources of Lamu Island, Kenya. Water 2015, 7, 1264–1290. [Google Scholar] [CrossRef]
  3. Hoekstra, A.Y.; Chapagain, A.K. Water footprints of nations: Water use by people as a function of their consumption pattern. In Integrated Assessment of Water Resources and Global Change; Springer: Dordrecht, The Netherlands, 2006; pp. 35–48. [Google Scholar] [CrossRef]
  4. Wada, Y.; Bierkens, M.F.; De Roo, A.; Dirmeyer, P.A.; Famiglietti, J.S.; Hanasaki, N.; Wheater, H. Human–water interface in hydrological modelling: Current status and future directions. Hydrol. Earth Syst. Sci. 2017, 21, 4169–4193. [Google Scholar] [CrossRef]
  5. Boretti, A.; Rosa, L. Reassessing the projections of the World Water Development Report. NPJ Clean Water 2019, 2, 15. [Google Scholar] [CrossRef]
  6. Pykhtin, A.I. Problems of rational use and protection of groundwater within the Russian Federation. J. Appl. Eng. Sci. 2019, 17, 3. [Google Scholar]
  7. Gleeson, T.; Wada, Y.; Bierkens, M.F.; Van Beek, L.P. Water balance of global aquifers revealed by groundwater footprint. Nature 2012, 488, 197–200. [Google Scholar] [CrossRef] [PubMed]
  8. World Bank. India Groundwater: A Valuable but Diminishing Resource. 2012. Available online: https://www.worldbank.org/en/news/feature/2012/03/06/india-groundwater-critical-diminishing (accessed on 6 November 2020).
  9. MoWR. Annual Report 2020-21. 2021. Available online: https://dpiit.gov.in/sites/default/files/annualReport-English2020-21.pdf (accessed on 10 April 2021).
  10. Srinivasan, V.; Thompson, S.; Madhyastha, K.; Penny, G.; Jeremiah, K.; Lele, S. Why is the Arkavathy River drying? A multiple-hypothesis approach in a data-scarce region. Hydrol. Earth Syst. Sci. 2015, 19, 1905–1917. [Google Scholar] [CrossRef]
  11. Penny, G.; Srinivasan, V.; Dronova, I.; Lele, S.; Thompson, S. Spatial characterization of long-term hydrological change in the Arkavathy watershed adjacent to Bangalore, India. Hydrol. Earth Syst. Sci. 2018, 22, 595–610. [Google Scholar] [CrossRef]
  12. Grasso, V. 2020 State of Climate Services Report: Risk Information and Early Warning Systems. In Proceedings of the EGU General Assembly Conference Abstracts 2021, EGU21-74, Vienna, Austria, 19–30 April 2021. [Google Scholar]
  13. Wada, Y.; van Beek, L.P.; Bierkens, M.F. Nonsustainable groundwater sustaining irrigation: A global assessment. Water Resour. Res. 2012, 48, W00L06. [Google Scholar] [CrossRef]
  14. MacDonald, A.M.; Bonsor, H.C.; Ahmed, K.M.; Burgess, W.G.; Basharat, M.; Calow, R.C.; Dixit, A.; Foster, S.S.D.; Gopal, K.; Lapworth, D.J.; et al. Groundwater quality and depletion in the Indo-Gangetic Basin mapped from in situ observations. Nat. Geosci. 2016, 9, 762–766. [Google Scholar] [CrossRef]
  15. Mishra, V.; Asoka, A.; Vatta, K.; Lall, U. Groundwater depletion and associated CO2 emissions in India. Earth’s Future 2018, 6, 1672–1681. [Google Scholar] [CrossRef]
  16. Luo, P.; Kang, S.; Zhou, M.; Lyu, J.; Aisyah, S.; Binaya, M.; Regmi, R.K.; Nover, D. Water quality trend assessment in Jakarta: A rapidly growing Asian megacity. PLoS ONE 2019, 14, e0219009. [Google Scholar] [CrossRef]
  17. Russo, T.A.; Lall, U. Depletion and response of deep groundwater to climate-induced pumping variability. Nat. Geosci. 2017, 10, 105–108. [Google Scholar] [CrossRef]
  18. MacDonald, A.; Davies, J.; Calow, R.; Chilton, J. Developing groundwater: A guide for rural water supply. In NERC Open Research Archive; ITDG Publishing: Rugby, UK, 2005. [Google Scholar]
  19. Porsani, J.L.; Elis, V.R.; Hiodo, F.Y. Geophysical investigations for the characterization of fractured rock aquifers in Itu, SE Brazil. J. Appl. Geophys. 2005, 57, 119–128. [Google Scholar] [CrossRef]
  20. Anomohanran, O.; Ofomola, M.O.; Okocha, F.O. Investigation of groundwater in parts of Ndokwa District in Nigeria using geophysical logging and electrical resistivity methods: Implications for groundwater exploration. J. Afr. Earth Sci. 2017, 129, 108–116. [Google Scholar] [CrossRef]
  21. Cook, P.G.; Walker, G.R.; Jolly, I.D. Spatial variability of groundwater recharge in a semiarid region. J. Hydrol. 1989, 111, 195–212. [Google Scholar] [CrossRef]
  22. Selvam, S.; Sivasubramanian, P. Groundwater potential zone identification using geoelectrical survey: A case study from Medak district, Andhra Pradesh, India. Int. J. Geomat. Geosci. 2012, 3, 55. [Google Scholar]
  23. Solomon, S.; Quiel, F. Groundwater study using remote sensing and geographic information systems (GIS) in the central highlands of Eritrea. Hydrogeol. J. 2006, 14, 1029–1041. [Google Scholar] [CrossRef]
  24. Kudamnya, E.A.; Edet, A.E.; Ekwere, A.S. Advances in Sustainable and Environmental Hydrology, Hydrogeology, Hydrochemistry and Water Resources. In Proceedings of the 1st Springer Conference of the Arabian Journal of Geosciences (CAJG-1), Tunisia, North Africa, 25–28 November 2019; pp. 249–252. [Google Scholar]
  25. Ozdemir, A. GIS-based groundwater spring potential mapping in the Sultan Mountains (Konya, Turkey) using frequency ratio, weights of evidence and logistic regression methods and their comparison. J. Hydrol. 2011, 411, 290–308. [Google Scholar] [CrossRef]
  26. Pourtaghi, Z.S.; Pourghasemi, H.R. GIS-based groundwater spring potential assessment and mapping in the Birjand Township, southern Khorasan Province, Iran. Hydrogeol. J. 2014, 22, 643–662. [Google Scholar] [CrossRef]
  27. Chenini, I.; Ben Mammou, A. Groundwater recharge study in arid region: An approach using GIS techniques and numerical modeling. Comput. Geosci. 2010, 36, 801–817. [Google Scholar] [CrossRef]
  28. Corsini, A.; Cervi, F.; Ronchetti, F. Weight of evidence and artificial neural networks for potential groundwater spring mapping: An application to the Mt. Modino area (Northern Apennines, Italy). Geomorphology 2009, 111, 79–87. [Google Scholar] [CrossRef]
  29. Naghibi, S.A.; Pourghasemi, H.R.; Pourtaghi, Z.S.; Rezaei, A. Groundwater qanat potential mapping using frequency ratio and Shannon’s entropy models in the Moghan watershed, Iran. Earth Sci. Inform. 2015, 8, 171–186. [Google Scholar] [CrossRef]
  30. Srivastava, A.; Sahoo, B.; Raghuwanshi, N.S.; Singh, R. Evaluation of Variable-Infiltration Capacity Model and MODIS-Terra Satellite-Derived Grid-Scale Evapotranspiration Estimates in a River Basin with Tropical Monsoon-Type Climatology. J. Irrig. Drain. Eng. 2017, 143, 04017028. [Google Scholar] [CrossRef]
  31. Elbeltagi, A.; Kumari, N.; Dharpure, J.K.; Mokhtar, A.; Alsafadi, K.; Kumar, M.; Mehdinejadiani, B.; Ramezani Etedali, H.; Brouziyne, Y.; Towfiqul Islam, A.R.M.; et al. Prediction of combined terrestrial evapotranspiration index (Ctei) over large river basin based on machine learning approaches. Water 2021, 13, 547. [Google Scholar] [CrossRef]
  32. Saraf, A.K.; Choudhury, P.R. Integrated remote sensing and gis for groundwater exploration and identification of artificial recharge sites. Int. J. Remote Sens. 1998, 19, 1825–1841. [Google Scholar] [CrossRef]
  33. Thakur, G.S.; Raghuwanshi, R.S. Perspect and assessment of groundwater resources using remote sensing techniques in and around Choral river basin, Indore and Khargone districts, M.P. J. Indian Soc. Remote Sens. 2008, 36, 217–225. [Google Scholar] [CrossRef]
  34. Dar, I.A.; Sankar, K.; Dar, M.A. Deciphering groundwater potential zones in hard rock terrain using geospatial technology. Environ. Monit. Assess. 2011, 173, 597–610. [Google Scholar] [CrossRef]
  35. Oikonomidis, D.; Dimogianni, S.; Kazakis, N.; Voudouris, K. A GIS/Remote Sensing-based methodology for groundwater potentiality assessment in Tirnavos area, Greece. J. Hydrol. 2015, 525, 197–208. [Google Scholar] [CrossRef]
  36. Magesh, N.S.; Chandrasekar, N.; Soundranayagam, J.P. Delineation of groundwater potential zones in Theni district, Tamil Nadu, using remote sensing, GIS and MIF techniques. Geosci. Front. 2012, 3, 189–196. [Google Scholar] [CrossRef]
  37. Saravanan, S.; Pitchaikani, S.; Venkatesan, G. Assessment and evaluation of groundwater vulnerability index maps of Upper Palar River Basin, Tamilnadu, India. J. Earth Syst. Sci. 2020, 129, 1–13. [Google Scholar] [CrossRef]
  38. Gumma, M.K.; Pavelic, P. Mapping of groundwater potential zones across Ghana using remote sensing, geographic information systems, and spatial modeling. Environ. Monit. Assess. 2013, 185, 3561–3579. [Google Scholar] [CrossRef] [PubMed]
  39. Senanayake, I.P.; Dissanayake, D.M.D.O.K.; Puswewala, U.G.A. Analysis of the abundance of abandoned tanks in Hambantota District, Sri Lanka using GIS techniques. Egypt. J. Remote Sens. Space Sci. 2012, 15, 143–150. [Google Scholar] [CrossRef]
  40. Kumar, T.; Gautam, A.K.; Kumar, T. Appraising the accuracy of GIS-based Multi-criteria decision making technique for delineation of Groundwater potential zones. Water Resour. Manag. 2014, 28, 4449–4466. [Google Scholar] [CrossRef]
  41. Rahmati, O.; Nazari Samani, A.; Mahdavi, M.; Pourghasemi, H.R.; Zeinivand, H. Groundwater potential mapping at Kurdistan region of Iran using analytic hierarchy process and GIS. Arab. J. Geosci. 2015, 8, 7059–7071. [Google Scholar] [CrossRef]
  42. Razandi, Y.; Pourghasemi, H.R.; Neisani, N.S.; Rahmati, O. Application of analytical hierarchy process, frequency ratio, and certainty factor models for groundwater potential mapping using GIS. Earth Sci. Inform. 2015, 8, 867–883. [Google Scholar] [CrossRef]
  43. Jenifer, M.A.; Jha, M.K. Comparison of Analytic Hierarchy Process, Catastrophe and Entropy techniques for evaluating groundwater prospect of hard-rock aquifer systems. J. Hydrol. 2017, 548, 605–624. [Google Scholar] [CrossRef]
  44. Kamaraju, M.V.V.; Bhattacharya, A.; Reddy, G.S.; Rao, G.C.; Murthy, G.S.; Rao, T.C.M. Ground–Water Potential Evaluation of West Godavari District, Andhra Pradesh State, India—A GIS Approach. Ground Water 1996, 34, 318–325. [Google Scholar] [CrossRef]
  45. Shaban, A.; Khawlie, M.; Abdallah, C. Use of remote sensing and GIS to determine recharge potential zones: The case of Occidental Lebanon. Hydrogeol. J. 2011, 14, 433–443. [Google Scholar] [CrossRef]
  46. Manap, M.A.; Sulaiman, W.N.A.; Ramli, M.F.; Pradhan, B.; Surip, N. A knowledge-driven GIS modeling technique for groundwater potential mapping at the Upper Langat Basin, Malaysia. Arab. J. Geosci. 2013, 6, 1621–1637. [Google Scholar] [CrossRef]
  47. Trabelsi, F.; Lee, S.; Khlifi, S.; Arfaoui, A. Frequency Ratio Model for Mapping Groundwater Potential Zones Using GIS and Remote Sensing; Medjerda Watershed Tunisia. In Advances in Sustainable and Environmental Hydrology, Hydrogeology, Hydrochemistry and Water Resources; Springer: Berlin/Heidelberg, Germany, 2019; pp. 341–345. [Google Scholar] [CrossRef]
  48. Saaty, T.L. The Analytic Hierarchy Process. In European Journal of Operational Research; McGraw-Hill: New York, NY, USA, 1982; Volume 9, pp. 97–98. [Google Scholar]
  49. Rajani, R.; Palanivel, K.; Selvakumar, R. Delineation of Groundwater Potential Zones Using Subsurface Parameters: An ANP Approach in and Around Tiruttani Area, Tamilnadu, India. Int. J. Pure Appl. Math. 2018, 119, 575–587. [Google Scholar]
  50. Siva, G.; Nasir, N.; Selvakumar, R. Delineation of Groundwater Potential Zone in Sengipatti for Thanjavur District using Analytical Hierarchy Process. In Proceedings of the IOP Conference Series: Earth and Environmental Science, Sanya, China, 20–22 November 2017; Volume 80. [Google Scholar]
  51. Sujatha, E.R. A spatial model for the assessment of debris flow susceptibility along the Kodaikkanal-Palani traffic corridor. Front. Earth Sci. 2020, 14, 326–343. [Google Scholar] [CrossRef]
  52. Gebru, H.; Gebreyohannes, T.; Hagos, E. Identification of Groundwater Potential Zones Using Analytical Hierarchy Process (AHP) and GIS-Remote Sensing Integration, the Case of Golina River Basin, Northern Ethiopia. Int. J. Adv. Remote Sens. GIS 2020, 9, 3289–3311. [Google Scholar] [CrossRef]
  53. Meade, L.; Sarkis, J. Strategic analysis of logistics and supply chain management systems using the analytical network process. Transp. Res. Part E Logist. Transp. Rev. 1998, 34E, 201–215. [Google Scholar] [CrossRef]
  54. Sujatha, E.R.; Sridhar, V. Mapping debris flow susceptibility using analytical network process in Kodaikkanal Hills, Tamil Nadu (India). J. Earth Syst. Sci. 2017, 126, 116. [Google Scholar] [CrossRef]
  55. Saraf, A.K.; Choudhury, P.R.; Roy, B.; Sarma, B.; Vijay, S.; Choudhury, S. GIS based surface hydrological modelling in identification of groundwater recharge zones. Int. J. Remote Sens. 2004, 25, 5759–5770. [Google Scholar] [CrossRef]
  56. NRSC, I. Manual of National Land Use/Land Cover Mapping (Second Cycle) Using Multi Temporal Satellite Data; Department of Space, Hyderabad: Hyderabad, India, 2012.
  57. Abdalla, F. Mapping of groundwater prospective zones using remote sensing and GIS techniques: A case study from the Central Eastern Desert, Egypt. J. Afr. Earth Sci. 2012, 70, 8–17. [Google Scholar] [CrossRef]
  58. Krishnamurthy, J.; Mani, A.; Jayaraman, V.; Manivel, M. Groundwater resources development in hard rock terrain—An approach using remote sensing and GIS techniques. ITC J. 2000, 2, 204–215. [Google Scholar] [CrossRef]
  59. Kale, V.S.; Gupta, A. Introduction to Geomorphology; Orient Longman: Hyderabad, India, 2001. [Google Scholar]
  60. Rajaveni, S.P.; Brindha, K.; Elango, L. Geological and geomorphological controls on groundwater occurrence in a hard rock region. Appl. Water Sci. 2017, 7, 1377–1389. [Google Scholar] [CrossRef]
  61. Yeh, H.F.; Cheng, Y.S.; Lin, H.I.; Lee, C.H. Mapping groundwater recharge potential zone using a GIS approach in Hualian River, Taiwan. Sustain. Environ. Res. 2016, 26, 33–43. [Google Scholar] [CrossRef]
  62. Duiker, S.W.; Flanagan, D.C.; Lal, R. Erodibility and infiltration characteristics of five major soils of southwest Spain. Catena 2001, 45, 103–121. [Google Scholar] [CrossRef]
  63. Ballukraya, P.N.; Sharma, K.K. Estimation of Storativity from Recovery Data. Groundwater 1991, 29, 495–498. [Google Scholar] [CrossRef]
  64. Banton, O.; Bangoy, L.M. A New Method to Determine Storage Coefficient from Pumping Test Recovery Data. Ground Water 1996, 34, 772–777. [Google Scholar] [CrossRef]
  65. Sanders, L.L. Field Hydrogeology; Prentice Hall: Englewood Cliffs, NJ, USA, 1998. [Google Scholar]
  66. Hassan, H.; Waru, S.M.; Bukar, G.A.; Abdullahi, K.M. Groundwater Potentials Estimation of a Basement Terrain Using Pumping Test Data for Parts of Sanga Local Government Area, Kaduna State, Northwestern Nigeria. Open J. Mod. Hydrol. 2016, 06, 222–229. [Google Scholar] [CrossRef]
  67. Suganthi, S.; Elango, L.; Subramanian, S.K. Groundwater potential zonation by remote sensing and GIS techniques and its relation to the groundwater level in the coastal part of the Arani and Koratalai river basin, Southern India. Earth Sci. Res. J. 2013, 17, 87–95. [Google Scholar]
  68. Da Costa, A.M.; de Salis, H.H.C.; Viana, J.H.M.; Pacheco, F.A.L. Groundwater recharge potential for sustainable water use in urban areas of the Jequitiba River Basin, Brazil. Sustainability 2019, 11, 2955. [Google Scholar] [CrossRef]
  69. Hasan, M.; Shang, Y.; Jin, W.; Akhter, G. Estimation of hydraulic parameters in a hard rock aquifer using integrated surface geoelectrical method and pumping test data in southeast Guangdong, China. Geosci. J. 2021, 25, 223–242. [Google Scholar] [CrossRef]
  70. Sridharan, K.; Mohan Kumar, M.S.; Sekhar, M. Groundwater flow and storage in hard rocks aquifers. In Proceedings of the Workshop on Artificial Recharge of Groundwater in Granite Terrain, Bangalore, India; 1995. [Google Scholar]
  71. Summers, W.K. Specific Capacities of Wells in Crystalline Rocks. Groundwater 1972, 10, 37–47. [Google Scholar] [CrossRef]
  72. Feibert, D.; Sørup, C.M.; Jacobsen, P. Using the Analytic Network Process (ANP) to assess the distribution of pharmaceuticals in hospitals—A comparative case study of a Danish and American hospital. In Proceedings of the 5th World Conference on Production and Operations Management, Havana, Cuba, 6–10 September 2016; pp. 1–10. [Google Scholar]
  73. Hazell, J.R.T.; Cratchley, C.R.; Jones, C.R.C. The hydrogeology of crystalline aquifers in northern Nigeria and geophysical techniques used in their exploration. Geol. Soc. Lond. Spec. Publ. 1992, 66, 155–182. [Google Scholar] [CrossRef]
  74. Chilton, J.; Chorus, I. Protecting groundwater for health. In Managing the Quality of Drinking Water Sources; IWA Publishers: London, UK, 2006; pp. 20–24. [Google Scholar]
  75. Verma, N.; Patel, R.K. Delineation of groundwater potential zones in lower Rihand River Basin, India using geospatial techniques and AHP. Egypt. J. Remote Sens. Space Sci. 2021, 24, 559–570. [Google Scholar] [CrossRef]
  76. Thapa, R.; Gupta, S.; Guin, S.; Kaur, H. Assessment of groundwater potential zones using multi-influencing factor (MIF) and GIS: A case study from Birbhum district, West Bengal. Appl. Water Sci. 2017, 7, 4117–4131. [Google Scholar] [CrossRef]
  77. Ramaiah, S.N.; Gopalakrishna, G.S.; Srinivasa Vittala, S.; Najeeb, K.M. Geomorphological mapping for identification of ground water potential zones in hard rock areas using geo-spatial information—A case study in Malur Taluk, Kolar District, Karnataka, India. Nat. Environ. Pollut. Technol. 2012, 11, 369–376. [Google Scholar]
  78. Varughese, A.; Suhail, A.; Chitra, M.G.; Jiji, P.S.; Deepthy, C. Identification of shallow groundwater potential zones using GIS—A case study. Int. J. Adv. Eng. Appl. 2012, 1, 65–70. [Google Scholar]
  79. Das, S. Delineation of groundwater potential zone in hard rock terrain in Gangajalghati block, Bankura district, India using remote sensing and GIS techniques. Modeling Earth Syst. Environ. 2017, 3, 1589–1599. [Google Scholar] [CrossRef] [Green Version]
  80. Jaafarzadeh, M.S.; Tahmasebipour, N.; Haghizadeh, A.; Pourghasemi, H.R.; Rouhani, H. Groundwater recharge potential zonation using an ensemble of machine learning and bivariate statistical models. Sci. Rep. 2021, 11, 5587. [Google Scholar] [CrossRef]
  81. Rashid, M.; Lone, M.A.; Ahmed, S. Integrating geospatial and ground geophysical information as guidelines for groundwater potential zones in hard rock terrains of south India. Environ. Monit. Assess. 2012, 184, 4829–4839. [Google Scholar] [CrossRef]
  82. Bagyaraj, M.; Ramkumar, T.; Venkatramanan, S.; Gurugnanam, B. Application of remote sensing and GIS analysis for identifying groundwater potential zone in parts of Kodaikanal Taluk, South India. Front. Earth Sci. 2013, 7, 65–75. [Google Scholar] [CrossRef]
  83. Lentswe, G.B.; Molwalefhe, L. Delineation of potential groundwater recharge zones using analytic hierarchy process-guided GIS in the semi-arid Motloutse watershed, eastern Botswana. J. Hydrol. Reg. Stud. 2020, 28, 100674. [Google Scholar] [CrossRef]
  84. Raji, W.O.; Abdulkadir, K.A. Evaluation of groundwater potential of bedrock aquifers in Geological Sheet 223 Ilorin, Nigeria, using geo-electric sounding. Appl. Water Sci. 2020, 10, 220. [Google Scholar] [CrossRef]
  85. Chatterjee, S.; Dutta, S. Assessment of groundwater potential zone for sustainable water resource management in south-western part of Birbhum District, West Bengal. Appl. Water Sci. 2022, 12, 40. [Google Scholar] [CrossRef]
  86. Machiwal, D.; Rangi, N.; Sharma, A. Integrated knowledge-and data-driven approaches for groundwater potential zoning using GIS and multi-criteria decision making techniques on hard-rock terrain of Ahar catchment, Rajasthan, India. Environ. Earth Sci. 2015, 73, 1871–1892. [Google Scholar] [CrossRef]
  87. Akhter, G.; Hasan, M. Determination of aquifer parameters using geoelectrical sounding and pumping test data in Khanewal District, Pakistan. Open Geosci. 2016, 8, 630–638. [Google Scholar] [CrossRef]
  88. Kumar, T.; Balasubramanian, A.; Kumar, R.S.; Dushiyanthan, C.; Thiruneelakandan, B.; Suresh, R.; Karthikeyan, K.; Davidraju, D. Assessment of groundwater potential based on aquifer properties of hard rock terrain in the Chittar–Uppodai watershed, Tamil Nadu, India. Appl. Water Sci. 2016, 6, 179–186. [Google Scholar] [CrossRef]
  89. Nsiah, E.; Appiah-Adjei, E.K.; Adjei, K.A. Hydrogeological delineation of groundwater potential zones in the Nabogo basin, Ghana. J. Afr. Earth Sci. 2018, 143, 1–9. [Google Scholar] [CrossRef]
  90. Dhakate, R. Distribution of aquifer characteristics in different geomorphologic units in a granitic terrain. Appl. Water Sci. 2020, 10, 231. [Google Scholar] [CrossRef]
  91. Hasan, M.; Shang, Y.; Jin, W.; Akhter, G. Joint geophysical prospecting for groundwater exploration in weathered terrains of South Guangdong, China. Environ. Monit. Assess. 2021, 193, 1–24. [Google Scholar] [CrossRef]
  92. Masoud, A.M.; Pham, Q.B.; Alezabawy, A.K.; El-Magd, S.A.A. Efficiency of geospatial technology and multi-criteria decision analysis for groundwater potential mapping in a Semi-Arid region. Water 2022, 14, 882. [Google Scholar] [CrossRef]
  93. Central Ground Water Board. Manual on Artificial Recharge of Groundwater; Ministry of Water Resource Board, Government of India: Delhi, India, 2007.
  94. Garewal, S.K.; Vasudeo, A.D.; Landge, V.S.; Ghare, A.D. A GIS-based Modified DRASTIC (ANP) method for assessment of groundwater vulnerability: A case study of Nagpur city, India. Water Qual. Res. J. 2015, 52, 121–135. [Google Scholar] [CrossRef]
  95. Owolabi, S.T.; Madi, K.; Kalumba, A.M.; Orimoloye, I.R. A groundwater potential zone mapping approach for semi-arid environments using remote sensing (RS), geographic information system (GIS), and analytical hierarchical process (AHP) techniques: A case study of Buffalo catchment, Eastern Cape, South Africa. Arab. J. Geosci. 2020, 13, 1–17. [Google Scholar] [CrossRef]
  96. Singha, S.S.; Pasupuleti, S.; Singha, S.; Singh, R.; Venkatesh, A.S. Analytic network process based approach for delineation of groundwater potential zones in Korba district, Central India using remote sensing and GIS. Geocarto Int. 2021, 36, 1489–1511. [Google Scholar] [CrossRef]
Figure 1. Map showing study area and geo-position of rain gauge stations, geophysical survey stations and pump test locations.
Figure 1. Map showing study area and geo-position of rain gauge stations, geophysical survey stations and pump test locations.
Water 14 02743 g001
Figure 2. Flow chart showing methodology.
Figure 2. Flow chart showing methodology.
Water 14 02743 g002
Figure 3. General structure of the super matrix.
Figure 3. General structure of the super matrix.
Water 14 02743 g003
Figure 4. Thematic layers on surface groundwater influencing parameters.
Figure 4. Thematic layers on surface groundwater influencing parameters.
Water 14 02743 g004
Figure 5. (af) Thematic layers on subsurface groundwater influencing parameters.
Figure 5. (af) Thematic layers on subsurface groundwater influencing parameters.
Water 14 02743 g005
Figure 6. Groundwater potential zone maps based on (a) surface parameters, (b) subsurface parameters, and (c) amalgamated parameters along with geo-positions of inventory wells.
Figure 6. Groundwater potential zone maps based on (a) surface parameters, (b) subsurface parameters, and (c) amalgamated parameters along with geo-positions of inventory wells.
Water 14 02743 g006
Table 1. Saaty’s scale of relative importance.
Table 1. Saaty’s scale of relative importance.
ImportanceScale
Extreme importance9
Very, very strong8
Very strong importance7
Strong6
Strong importance5
Moderate plus4
Moderate importance3
Weak2
Equal importance1
Table 2. Random consistency index.
Table 2. Random consistency index.
n123456789101112131415
RI000.580.901.121.241.321.411.451.491.511.531.561.571.59
Table 3. Ideal values and overall ranking of surface parameters derived through supermatrix computation.
Table 3. Ideal values and overall ranking of surface parameters derived through supermatrix computation.
Sl.NoAlternativesTotalNormalIdealRanking
1Floodplain0.05730.131111
2Cropland0.04460.1020.77822
3Pediplain0.04230.09670.73773
4DD *—low0.03210.07340.55994
5LD *—high0.03120.07130.54415
6Water bodies0.02380.05450.41546
7Gentle slope0.02310.0530.40387
8Entisols0.02270.05190.39588
9Gneiss0.02180.050.38129
10Charnokite0.01480.03390.258610
11Pediment0.01270.0290.221111
12Fallow land0.01260.02890.220512
13Vertisols0.01210.02770.211413
14Medium slope0.0110.02510.191714
15LD medium0.00920.02110.16115
16Built-up land0.00870.01980.151216
17DD medium0.00840.01910.145917
18Sandstone0.00760.01730.132218
19Alfisols0.0060.01370.104819
20Dyke0.00530.01220.093420
21Wasteland0.00530.01210.092121
22Inceptisols0.00520.0120.091622
23LD—low0.00510.01160.088223
24Structural hill0.00440.010.076624
25Granite0.00330.00760.057825
26DD—high0.00330.00750.05726
27Steep slope0.00320.00720.055327
* LD—lineament density; DD—drainage density.
Table 4. Ideal values and the overall ranking of subsurface parameters derived through supermatrix computation.
Table 4. Ideal values and the overall ranking of subsurface parameters derived through supermatrix computation.
Sl.No.AlternativesTotalNormalIdealRanking
1Aquifer thickness—deep0.08870.208511
2Water level—shallow0.06140.14440.69292
3Permeability—high0.04860.11430.54843
4Specific capacity—high0.03630.08530.40944
5Storage coefficient—high0.0320.07530.36125
6Transmissivity—high0.02870.06750.32376
7Transmissivity—moderate0.01930.04530.21757
8Aquifer thickness—moderate0.01630.03830.18388
9Specific capacity—moderate0.01560.03670.17619
10Water level—moderate0.01410.03310.158610
11Storage coefficient—moderate0.01290.03040.14611
12Permeability—moderate0.01160.02730.130812
13Specific capacity—low0.00970.02280.109213
14Water level—deep0.00840.01980.095114
15Aquifer thickness—shallow0.00710.01670.080315
16Transmissivity—low0.00510.01210.05816
17Storage coefficient—low0.00480.01130.054217
18Permeability—low0.00460.01090.052118
Table 5. Frequency ratio between water level and surface parameter-based groundwater potential zones.
Table 5. Frequency ratio between water level and surface parameter-based groundwater potential zones.
Potential Zone CategoryNo. of PixelsCP/TPNo. of WellsW/TWFR
Shallow water level
High466,4830.435734929130.3714290.852419
Moderate374,6080.34991583910.0285710.081652
Low229,4750.214349232000
Moderate water level
High466,4830.43573492950.1428570.327853
Moderate374,6080.34991583970.20.571566
Low229,4750.21434923210.0285710.133294
Deep water level
High466,4830.43573492930.0857140.196712
Moderate374,6080.34991583930.0857140.244957
Low229,4750.21434923220.0571430.266588
Total1,070,566 35
Table 6. Frequency ratio between water level and subsurface parameter-based groundwater potential zones.
Table 6. Frequency ratio between water level and subsurface parameter-based groundwater potential zones.
Potential Zone CategoryNo. of PixelsCP/TPNo. of WellsW/TWFR
Shallow water level
High410,2250.383185110.3142860.820193
Moderate466,3410.43560210.0285710.065591
Low194,0000.18121320.0571430.315336
Moderate water level
High410,2250.38318550.1428570.372815
Moderate466,3410.43560240.1142860.262363
Low194,0000.18121340.1142860.630672
Deep water level
High410,2250.38318520.0571430.149126
Moderate466,3410.43560210.0285710.065591
Low194,0000.18121350.1428570.78834
Total1,070,566 35
Table 7. Frequency ratio between water level and amalgamated groundwater potential zones.
Table 7. Frequency ratio between water level and amalgamated groundwater potential zones.
Potential Zone CategoryNo. of PixelsCP/TPNo. of WellsW/TWFR
Shallow water level
High389,4270.363758130.3714291.021087
Moderate484,9590.45299310.0285710.063073
Low196,1800.183249000
Moderate water level
High389,4270.363758000
Moderate484,9590.452993120.3428570.756871
Low196,1800.18324910.0285710.155916
Deep water level
High389,4270.363758000
Moderate484,9590.452993000
Low196,1800.18324980.2285711.247328
Total1,070,566 35
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Radhakrishnan, S.; Ramachandran, R.; Murali, G.; Vatin, N.I. A Hybrid Spatial–Analytical Network Process Model for Groundwater Inventory in a Semi-Arid Hard Rock Aquifer System—A Case Study. Water 2022, 14, 2743. https://doi.org/10.3390/w14172743

AMA Style

Radhakrishnan S, Ramachandran R, Murali G, Vatin NI. A Hybrid Spatial–Analytical Network Process Model for Groundwater Inventory in a Semi-Arid Hard Rock Aquifer System—A Case Study. Water. 2022; 14(17):2743. https://doi.org/10.3390/w14172743

Chicago/Turabian Style

Radhakrishnan, Selvakumar, Rajani Ramachandran, Gunasekaran Murali, and Nikolai Ivanovich Vatin. 2022. "A Hybrid Spatial–Analytical Network Process Model for Groundwater Inventory in a Semi-Arid Hard Rock Aquifer System—A Case Study" Water 14, no. 17: 2743. https://doi.org/10.3390/w14172743

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop