The sustainable management of groundwater resources is crucial for ecological diversity, human health, and economic growth. This study employs scientific concepts and advanced techniques, including the analytic hierarchy process (AHP) and Fuzzy-AHP, to identify groundwater potential zones (GWPZs). Thematic maps representing drainage density, elevation, soil, geomorphology, slope, land use and land cover, and rainfall are used to delineate the GWPZs. Both techniques are employed to assign weights to these thematic maps based on their characteristics and water potential. The study revealed that in the investigated area, 17.76 and 18.27% of the final GWPZs (AHP and Fuzzy-AHP) can be classified as having poor potential, while 72.79 and 71.07% are categorized as having moderate potential. Moreover, 9.45 and 10.69% of the final GWPZs are identified as having high potential using the AHP and Fuzzy-AHP models, respectively. Receiver operating characteristics (ROCs) analysis is employed to validate these findings, demonstrating that the Fuzzy-AHP technique achieves an accuracy of 74% in identifying GWPZs in the region. This study utilizes the best method derived from both models to identify 26 suitable locations for artificial recharge sites. The reliable findings of this research offer valuable insights into decision-makers and water users in the Kinnerasani Watershed.

  • The novel Fuzzy-AHP method for identifying GWPZs.

  • Suitable locations for artificial recharge sites were determined by selecting the best model between AHP and Fuzzy-AHP.

  • The obtained results are valuable for making informed decisions and facilitating sustainable groundwater management planning.

Water stands as the most consumed material on the Earth due to its indispensable role in supporting life (Zade et al. 2021; Nguyen et al. 2022; Roy et al. 2022). However, a significant portion (97%) of the Earth's water is salty, rendering it unsuitable for most human purposes (Das et al. 2019). Only a mere 3% of the planet's water comprises freshwater, which is the sole source appropriate for human use (Omar et al. 2020; Goswami & Ghosal 2022). Among these freshwater resources, groundwater, which refers to water stored in underground aquifers (Post & Abarca 2010; Nyakundi et al. 2022), plays a crucial role in safeguarding and preserving water within the Earth's crust (Fitts 2013). Given its wide-ranging applications in all aspects of life, groundwater stands as an important natural resource (Jha et al. 2007). Over time, the utilization of groundwater has significantly increased due to the overuse of surface water and the unpredictability of rainfall. In order to complement surface freshwater sources, households, industries, and agriculture have predominantly relied on groundwater supplies (Andualem & Demeke 2019). Groundwater serves as a crucial water source for approximately 2.5 billion people worldwide, fulfilling their domestic water requirements (UNESCO 2015). Many nations, particularly those situated in semi-arid and arid regions, heavily rely on groundwater for various purposes such as industry, agriculture, and household use (Serele et al. 2020; Ahmed et al. 2021). However, due to the escalating global population and rapid industrialization, there has been a substantial increase in the dependence on groundwater, leading to immense strain on this invaluable resource (Das et al. 2021). Consequently, it is crucial to establish effective delineation techniques and implement sustainable groundwater management practices to meet the rising demand, not only for the present generation but also for future generations (Goswami & Ghosal 2022). Identifying potential zones for groundwater recharge becomes essential for maintaining groundwater quality and ensuring the efficiency of the groundwater management system.

It can be difficult to chart this potential recharge zone because groundwater is typically an inconspicuous source (Allafta et al. 2021). Numerous geoenvironmental factors, including geomorphology, land use and land cover (LULC), lithology, climatic conditions, the extent of fractures, soil, drainage pattern, lineament features, secondary porosity, landform, topography, geological structures, slope, etc., and also their interactions, determine the presence and movement of groundwater in a given area (Greenbaum 1985; Gupta et al. 2020; Dibs et al. 2023).

The many geoenvironmental elements that have an impact on groundwater growth deep inside the Earth's crust make it challenging to identify groundwater. The hunt for groundwater has been answered over time using conventional techniques, such as geophysical examinations (Andualem & Demeke 2019; Agyemang 2022; Shekar & Mathew 2022b). Nevertheless, the drawbacks of these conventional techniques are that they require more money, labour, and time to complete (Gumma & Pavelic 2013; Shyam et al. 2021; Mathew et al. 2022). Geospatial techniques' more recent development has opened the door for creative and effective methods of groundwater exploration (Mandal et al. 2016). These innovative techniques, in contrast to the conventional ones, are quick, economical, and time-saving, and they can even be used in inaccessible locations (Serele et al. 2020; Shekar & Mathew 2022a). In order to analyze and integrate data to provide meaningful outputs for exploring, assessing, and managing essential groundwater resources, the Geographic Information Systems (GIS) and Remote Sensing (RS) technique has proven particularly important (Serele et al. 2020; Ijumulana et al. 2021; Omar et al. 2021; Shekar & Mathew 2022c; Thacker et al. 2022). Over the past few decades, multi-criteria decision analysis (MCDA) approaches, led by remote sensing and GIS, have become more and more popular for managing water resources and the environment (Allafta et al. 2021). Delineating the GWPZs globally is aided by integrated geospatial approaches, using several multi-criteria decision-making techniques. The study conducted by Qadir et al. (2020) aimed to map groundwater potential in a specific region using remote sensing and GIS techniques. The researchers integrated various thematic layers, including lineament density, slope, Topographic Wetness Index (TWI) geology, drainage density, and LULC. Each of these layers was assigned appropriate weights to account for their relative significance in determining groundwater potential using remote sensing and GIS techniques. According to the paper by Arunbose et al. (2021), the authors focused solely on using the GIS-based analytical hierarchy process (AHP) technique method for their investigation of groundwater potential zones (GWPZs).

GIS and AHP integration is a potent tool for combining the evaluation of several criteria (Kumar et al. 2022). The subjective Saaty (1980) technique, known as the AHP, allows users to select the relative weighting of each criterion in a solution to a problem that depends on many criteria. The distribution of impact among the various thematic levels can be determined using AHP, according to Saaty (1996) and Wind and Saaty (1980). It aids in the spatially complex challenge of environmental issues for decision-makers. To find a GWPZ, several characteristics have been taken into consideration (Allafta et al. 2021). Based on their relative weights in the process of deciding where to locate the groundwater zone, each characteristic's weights were assigned. For each problem in the AHP technique, a hierarchical model with goals, criteria, sub-criteria, and options is used (Saaty 2008). The weights or priority of the pairwise comparisons are used to resolve issues with AHP (Ho 2008). When comparing the criteria in pairs, the AHP approach is susceptible to discrepancies to a certain extent.

Fuzzy logic, a branch of artificial intelligence (AI), focuses on approximate reasoning rather than strict precision. It provides a mathematical framework to handle uncertainty and vagueness, acknowledging that things can be partially true or false and allowing for degrees of truth to be expressed through mathematical expressions (Dubois & Prade 2003). Fuzzy logic finds application in various AI domains, including control systems, pattern recognition, decision-making, and data analysis. It proves particularly useful in situations with insufficient or ambiguous information, where conventional logic may not provide clear solutions (Bezdek et al. 1999; Al-Jarrah & Halawani 2001). Recent research highlights the integration of AI, fuzzy logic, and remote sensing techniques for groundwater potential mapping. These methods enable the creation of accurate groundwater potential maps that support resource management and effective groundwater management. To achieve this, it is important to examine the logical consistency of pairwise comparisons (Cheng 1997; Chang et al. 2008). To address these challenges and improve decision-making confidence, Van-Laarhoven & Pedrycz (1983) introduced the Fuzzy-AHP method. Fuzzy-AHP is an extension of the AHP designed for hierarchical fuzzy multi-criteria decision-making. This approach combines alternative strategy selection, systematic weighing, and principles from fuzzy set theory and hierarchical structure analysis. By employing fuzzy numbers to represent weight values, the Fuzzy-AHP approach generates accurate and realistic results (Aryafar et al. 2013; Lee et al. 2013).

Multidisciplinary decision-making has been recognized over the past 20 years as a good method for managing water because it gives decisions more structure, readability, transparency, and firmness (Roberts & Goodwin 2002; Flug et al. 2005). Water resource management requires proper MCDA, which should concentrate on the basis of the decision-making process, including conditional decision-making and potential decisions. Integration with GIS-based techniques, which yield more accurate results and need less computation time than conventional GWPZ field techniques, is a major benefit of these methods (Zolekar & Bhagat 2015).

Groundwater resources play a vital role in the Kinnerasani River basin, located in the Telangana state of India. Encompassing an area of approximately 3,316 km2, the basin is predominantly agricultural, relying heavily on groundwater for irrigation of crops such as paddy, vegetables, cotton, and maize. However, the Kinnerasani River faces challenges, including insufficient water flow during dry seasons, which significantly impact both the river's ecosystem and the local communities that depend on it. Given the region's economic dependence on agriculture and groundwater, sustainable management of groundwater resources is crucial for ensuring the socio-economic development of the area. As per the author's knowledge, there has been a lack of investigations regarding the potential of groundwater in the current study area. Furthermore, there have been limited studies that have utilized the AHP and Fuzzy-AHP techniques. In order to bridge these gaps in knowledge and address these challenges, this article endeavours to identify the GWPZs within the Kinnerasani River basin by employing RS and GIS techniques. The objectives of this study are to develop a groundwater potential map in the Kinnerasani River basin by incorporating seven thematic layers: drainage density, elevation, soil, geomorphology, slope, LULC, and rainfall. The application of AHP and Fuzzy-AHP techniques will be employed to analyze and integrate these layers. The accuracy of the GWPZ map will be validated by comparing it with existing well data points using the receiver operating characteristic (ROC) and area under the curve (AUC) methods. Furthermore, this study aims to identify suitable locations for artificial recharge sites in the Kinnerasani River basin based on the generated groundwater potential map. By achieving these objectives, we hope to provide valuable insights into the potential distribution of groundwater resources in the study area and contribute to informed decision-making regarding groundwater management and sustainable development in the Kinnerasani River basin.

Study area

Figure 1 shows the location of the Kinnerasani Watershed. The Kinnerasani Watershed lies between latitudes 17°10′0ʺ and 18°10′0ʺ North and longitudes 80°10′0ʺ and 81°50′0ʺ East. The Kinnerasani River's mouth is situated at latitudes 17°38′15ʺ North and longitudes of 80°53′48ʺ East. The total area of the Kinnerasani River Basin is 3,316 km2. The Kinnerasani Watershed's elevation ranges from 38 to 784 m above sea level based on the Shuttle Radar Topography Mission-Digital Elevation Model (SRTM-DEM). The soil types in the current watershed are classified as clay loam and clay by the Food and Agricultural Organization (FAO) in their soil classification from 1988 (FAO 1988). Lower Triassic to Upper Carboniferous Sedimentary Rocks and Undivided Precambrian Rocks are the area's geologic classification, according to US Geological Survey geographic classification (1998). Sentinel-2 data show that crops were the most common type of land use, followed by trees, shrubs, constructed areas, water, bare ground, grass, and flooded vegetation.
Figure 1

Map of the Kinnerasani Watershed's location.

Figure 1

Map of the Kinnerasani Watershed's location.

Close modal
This study aims to identify a prospective groundwater availability zone in the Indian Kinnerasani Watershed. Groundwater potential maps are produced by analysing many factors, including drainage density, elevation, soil, geomorphology, slope, LULC, and rainfall, and then using weighted AHP and Fuzzy-AHP models. The schematic design (Figure 2) provides a detailed explanation of all the techniques used in this study. Table 1 lists the data sources used in the current work. The data from the sentinel satellite image were used to generate the LULC layer (ESRI land cover – https://livingatlas.arcgis.com/landcover/). Using information from the Geological Survey of India, geomorphology was created (https://bhukosh.gsi.gov.in/Bhukosh/Public). Using CRU data, rainfall data were produced (https://crudata.uea.ac.uk/cru/data/hrg/). The drainage density and elevation layers were produced using the SRTM-DEM. Utilizing SRTM data that was collected from the NASA website, the theme layer of the slope was created (https://earthexplorer.usgs.gov/). GWPZs were created for the Kinnerasani River basin based on AHP by combining all seven theme layers using the weighted overlay analysis technique on the GIS platform. The resultant GWPZs map indicated poor, moderate, and high zones. The effectiveness of AHP approaches was validated using ROC.
Table 1

Data used for this research

S. No.Thematic layersData sources
Geomorphology Geological Survey of India 
Drainage density SRTM-DEM data 
Soil FAO soil classification (1988
Slope SRTM-DEM data 
LULC Environmental Systems Research Institute (ESRI) 2021 Land Cove 
Rainfall distribution Climatic Research Unit 
Elevation SRTM-DEM data 
Lineament Density Geological Survey of India 
S. No.Thematic layersData sources
Geomorphology Geological Survey of India 
Drainage density SRTM-DEM data 
Soil FAO soil classification (1988
Slope SRTM-DEM data 
LULC Environmental Systems Research Institute (ESRI) 2021 Land Cove 
Rainfall distribution Climatic Research Unit 
Elevation SRTM-DEM data 
Lineament Density Geological Survey of India 
Figure 2

Flowchart for identifying GWPZs and artificial recharge sites.

Figure 2

Flowchart for identifying GWPZs and artificial recharge sites.

Close modal

Analytical hierarchy process

Numerous parameters and their subclasses have varied effects on how groundwater behaves, thus they must be assessed and weighed appropriately. This study uses an MCDM approach based on the AHP. AHP, an MCDM tool, is frequently used to manage groundwater (Rahmati & Melesse 2016). Each thematic map's weight was determined using the AHP method. After considerable research and expert advice, theme layers are given priority. The AHP technique, which is based on multi-criteria decision-making, establishes the weight of each thematic layer using Saaty's (1–9) scale. By contrasting the attributes in pairs, the comparison matrix for the seven thematic levels was developed. The eigenvector is computed to normalize the weight assigned to a layer (Saaty 1980). According to the above approach, taking into consideration earlier studies and the expert's knowledge, final weights were established for each thematic layer. GWPZs were created for the Kinnerasani River basin based on AHP by combining all seven theme layers using the weighted overlay analysis technique on the GIS platform. The resultant GWPZs map indicated poor, moderate, and high zones. The effectiveness of AHP approaches was validated using ROC.

Fuzzy-AHP

In decision-making processes, discrepancies often arise when weighing quantitative and qualitative parameters. To address this, methods such as MCDM techniques play a vital role. These techniques involve assigning weights to different parameters and ranking a subset of parameters based on professional judgement within a specific thematic context (Mohanty & Walling 2008; Quadrio et al. 2015). While the AHP is commonly used in MCDM, the Fuzzy-AHP approach developed by Van-Laarhoven & Pedrycz (1983) provides a subjective assessment of weighing multiple criteria, effectively dealing with differences that may emerge during pairwise comparisons (Chang 1996; Sener et al. 2018). In the present study, the AI-based Fuzzy-AHP approach was utilized to examine the GWPZ, considering seven parameters. To construct the triangular fuzzy judgement matrix, triangular fuzzy numbers were employed as a replacement for crisp numbers in the pairwise comparison matrix (Lyu et al. 2020). This fuzzy approach enables a more nuanced representation of uncertainties and enhances the comprehensiveness of the decision-making process. The geometric mean approach of Buckley (1985) was used to obtain the fuzzy geometric mean and fuzzy weights of each criterion. Defuzzification, or converting a fuzzy number into a crisp number, was ultimately accomplished by first determining the arithmetic mean and then the normalized weight. The combined seven theme layers were processed using the weighted overlay analysis technique in the GIS platform to generate the GWPZ map for the Kinnerasani River basin, employing the Fuzzy-AHP approach. The final GWPZ map identified poor, moderate, and high zones, representing varying levels of groundwater potential. To validate the results of the Fuzzy-AHP technique, the ROC analysis was employed.

Description of the thematic maps for groundwater potential mapping

In this study, the GWPZ map was developed using GIS by integrating multiple thematic layers that contribute to the occurrence of groundwater. The following is a concise overview of the thematic maps utilized for groundwater potential mapping, along with their corresponding weights assigned through AHP and Fuzzy-AHP, including validation and the resulting final GWPZ map.

Drainage density

The drainage density (Dd), a crucial indicator of the linear scale of landform components in stream-eroded topography, was first established by Horton (1932). The Dd is defined as the ratio of the catchment's area to the total of all streams of all kinds, given in units of km/km2. It quantifies the average spacing and proximity of stream channels throughout the entire basin, providing a measurable representation of their density (Horton 1932). Drainage density observations across diverse geological and climatic conditions reveal a consistent trend: regions characterized by highly permeable subsoil, dense vegetative cover, and low relief are more prone to exhibit low drainage density. High drainage density is influenced by factors such as low vegetative cover, impermeable underlying material, and hilly topography. The drainage density is classified into different categories: very low (0–0.49 km/km²), low (0.5–0.91 km/km²), moderate (0.92–1.3 km/km²), high (1.4–1.7 km/km²), and very high (1.8–3.1 km/km²). In the zonation of groundwater potential, greater emphasis is placed on areas with lower density (higher weight) compared with those with higher density (lower weight). Figure 3 illustrates the drainage density chart for the Kinnerasani River basin.
Figure 3

Drainage density map of Kinnerasani River basin.

Figure 3

Drainage density map of Kinnerasani River basin.

Close modal

Slope

One of the most important variables affecting groundwater potential is slope. In general, groundwater potential rises as slope values decline (Delin et al. 2000). Flat areas with low slopes may store precipitation and replenish the groundwater, whereas elevated areas with steep slopes will have significant runoff and little infiltration. The slope map for the basin was created using DEM with a 30 m resolution that was obtained from SRTM data (Figure 4). The slope 8 (degrees) values were categorized into five categories: fat (0–2.9), gentle (3–7.5), medium (7.6–15), steep (16–23), and very steep (24–26). Fat and gentle slopes are given more weight. For slopes that are steep as well as very steep, a low weight is provided.
Figure 4

Slope map of Kinnerasani River basin.

Figure 4

Slope map of Kinnerasani River basin.

Close modal

Geomorphology

The percolation and recharging of water into the Earth's subsurface are significantly influenced by geomorphology (Donselaar et al. 2017). The area's geomorphic units regulate the flow of groundwater below the surface. Any surface's potential to recharge is significantly influenced by geomorphology. Waterbody-river, alluvial plain, quarry and mine dump, low dissected hills and valleys, flood plain, pediment pediplain complex, moderately dissected hills and valleys, and waterbodies-other, dam and reservoir, high dissected hills, and valleys are the principal geomorphic units found in the study region (Figure 5).
Figure 5

Geomorphology map of Kinnerasani River basin.

Figure 5

Geomorphology map of Kinnerasani River basin.

Close modal

Soil

The type, thickness, texture, and composition of the soil all affect infiltration characteristics and, in turn, groundwater recharge, making soil an essential component in assessing GWPZ (Shit et al. 2015; Achu et al. 2020; Saade et al. 2021). A specific area's soil features, kinds, and distribution are influenced by factors such as geomorphology, geology, relief, time, and others. The degree of permeability, which is established by the interaction between infiltration rates, and runoff and the properties of the soil, defines the groundwater potential (Tesfaye 2010). The specifics of the soil classifications found in the basin by FAO soil (Figure 6). The clay loam soil is most usually found in the watershed, which acts as a region with a high capacity for groundwater. and clay soil serves as a region that is poorly permeable for replenishing groundwater.
Figure 6

Soil map of Kinnerasani River basin.

Figure 6

Soil map of Kinnerasani River basin.

Close modal

Elevation

Elevation plays a significant role in identifying potential groundwater areas. To incorporate elevation information, a DEM was utilized to generate the elevation data for the research area. In this study, the contribution of the Kinnerasani Watershed to groundwater potential and recharge was analyzed based on five elevation classes (Figure 7). The SRTM-DEM was used to obtain the elevation information for the research region. The DEM was then classified into five distinct classes based on their elevation ranges: very low (38–137.5 m), low (137.6–204.8 m), moderate (204.9–286.7 m), high (286.8–403.7 m), and extremely high (403.8–784 m). This classification allows for a better understanding and analysis of the elevation variations within the study area. In the groundwater potential zoning, low elevation areas were given a higher weightage, indicating their greater potential, while high elevation areas received a lower weightage. This weighting approach allows for a comprehensive assessment of the groundwater potential across different elevation zones.
Figure 7

Elevation map of Kinnerasani River basin.

Figure 7

Elevation map of Kinnerasani River basin.

Close modal

LULC

The LULC map gives accurate details on the area's geography as well as the different land use classifications (Martin et al. 2017). The amount of vegetation in a given area affects its capacity for recharging and surface runoff. The denser the vegetation, the higher the infiltration rate because vegetation enhances infiltration and reduces runoff. The map was created using information from Sentinel-2 imagery for the year 2020. (ESRI land cover). The resolution of Sentinel-2 is 10 m. LULC categories such as flooded vegetation, scrub/shrub, crops, bare ground, range land, trees, built-up area, water, and grass were considered (Figure 8).
Figure 8

LULC map of Kinnerasani River basin.

Figure 8

LULC map of Kinnerasani River basin.

Close modal

Rainfall

Rainfall is a significant source of groundwater recharge and is also a crucial part of the water cycle (Yu & Lin 2015). Groundwater recharging is influenced by rainfall distributions in a spatial–temporal manner. The increased rainfall will probably lead to an increase in groundwater potential. Little groundwater recharge may occur in a location with little rainfall. Rainfall duration and intensity have an effect on infiltration. While quick, high-intensity downpours boost surface runoff at the expense of infiltration, low intensity, long-duration rain impacts more infiltration than runoff. High rainfall receives high weights, and vice versa. Data on monthly precipitation for the year 2016 were gathered from the Climatic Research Unit. Rainfall varies from 1,080.3 to 1,147 mm in the year 2016. Using the kriging interpolation approach, a map of the spatial distribution of rainfall was created. In Figure 9, the amount of precipitation has been divided into five categories: very low (1,126–1,149 mm), low (1,150–1,165 mm), moderate (1,166–1,180 mm), high (1,181–1,196 mm), and very high (1,197–1,219 mm).
Figure 9

Rainfall map of Kinnerasani River basin.

Figure 9

Rainfall map of Kinnerasani River basin.

Close modal

Many thematic layers that favour the presence of groundwater have been combined to build a groundwater potential map using GIS. GWPZ mapping thematic maps utilized in this work, their weights determined by AHP, Fuzzy-AHP with validation, and the final GWPZ map are briefly described in the following.

Analytical hierarchical process

After developing theme layers for groundwater modelling, the next step involved weighing the parameters considering their varying impacts on aquifer replenishment (Serele et al. 2020). The AHP is widely recognized as a reliable method and recommended by many researchers (Roy et al. 2022). Saaty (1980) introduced the AHP as a powerful decision-making approach for determining the relative weights of different variables. In this study, the AHP was employed to evaluate the relative contributions of each factor to groundwater recharge. Saaty (1980) devised a scale of relative relevance ranging from 1 to 9, which was utilized to construct the matrix. The matrix was completed by comparing each pair of criteria and assigning rankings based on their relative importance. The diagonal elements of the matrix were set to 1, indicating equal weights with other diagonal elements. Criteria that exerted a significant influence over others were assigned a rating of 9. The pairwise matrix created for this study is presented in Table 2, whereas Table 3 shows the Saaty scale of relative relevance. According to Saaty (1990), the reciprocal values were utilized to fill the lower triangular matrix and the values were placed in the upper triangle to construct a diagonal matrix. Additionally, Table 4 shows that the relative weights were uniform. By determining the fundamental eigenvalue and examining the consistency ratio (CR), the thoughts' coherence was assessed (Saaty 1990).

Table 2

AHP-based pairwise comparison matrix

 
 
Table 3

Relative importance scale

IntensityDefinition
Equal importance 
Moderate 
Strong 
Very strong 
Extreme importance 
2,4,6,8 Intermediate values 
IntensityDefinition
Equal importance 
Moderate 
Strong 
Very strong 
Extreme importance 
2,4,6,8 Intermediate values 
Table 4

Pairwise comparison matrix with normalized principal eigenvector

 
 

The pairwise matrix of the principal eigenvalue matrix and the eigenvectors from Table 4 was combined to get the primary eigenvalue (max) of Table 5. For a 7 × 7 matrix, a primary eigenvalue of 7.68 was found, and the consistency index (CI) was calculated using this value (Table 5). Also standardized were the relative weights (Table 4). This fundamental eigenvalue is used to evaluate conceptual coherence by looking at the CR (Saaty 1990). Equation (1) is used to compute CI to obtain the CR:
(1)
where n is the number of parameters and max is the primary eigenvalue. The CI for this study is thus
Table 5

Calculating the primary eigenvalue

Column sums (row 9 of Table 4)Eigenvectors (column 9 of Table 4)Parameter rank
(1)(2)(1) × (2)
Geomorphology 3.83 0.23 0.89 
Drainage density 3.43 0.28 0.97 
Rainfall 6.20 0.18 1.14 
Elevation 11.40 0.12 1.31 
Slope 13.50 0.10 1.33 
LULC 21.33 0.05 1.11 
Soil 26.00 0.04 0.92 
Sum (λmax  7.68 
Column sums (row 9 of Table 4)Eigenvectors (column 9 of Table 4)Parameter rank
(1)(2)(1) × (2)
Geomorphology 3.83 0.23 0.89 
Drainage density 3.43 0.28 0.97 
Rainfall 6.20 0.18 1.14 
Elevation 11.40 0.12 1.31 
Slope 13.50 0.10 1.33 
LULC 21.33 0.05 1.11 
Soil 26.00 0.04 0.92 
Sum (λmax  7.68 
The reliability of the output was evaluated using the CR and CI values that were calculated (Equation (2)):
(2)
Changing the order of the matrix affects the values of the random consistency index (RI) and CR (Table 6; Saaty 1980):
Table 6

For different n scores, ratio index (RI)

N345678910
RI 0.58 0.89 1.12 1.24 1.32 1.41 1.45 1.49 
N345678910
RI 0.58 0.89 1.12 1.24 1.32 1.41 1.45 1.49 
The weighted overlay analysis conducted in ArcGIS was a suitable approach to map the GWPZs by considering the weighted features. The analysis took into account various thematic layers, each assigned a weight based on its importance in determining the groundwater potential. Among the thematic layers, drainage density emerged as the most influential indicator, with a weight of 28.1%. This suggests that areas with higher drainage density are more likely to have favourable groundwater potential. Following drainage density, geomorphology, rainfall, and elevation were found to have the next highest levels of influence, with weights of 23.7, 18.3, and 11.3% respectively. On the other hand, slope, LULC, and soil were considered to have the least impact on the delineation of groundwater recharge regions, with weights of 9.8, 5.3, and 3.6%, respectively. This indicates that these factors are less significant in determining groundwater potential compared with other variables. By overlaying the influential thematic layers using the weighted overlay analysis, a final GWPZs map was generated. This map, as depicted in Figure 10, provides a spatial representation of the GWPZs in the research area. Additionally, the map includes three zones: poor, moderate, and high, representing the different levels of groundwater potential. According to the AHP model, the research region was classified as 17.76% poor potential, 72.79% moderate potential, and 9.45% high potential, as shown in Figure 11. This indicates that the majority of the study area falls under the moderate GWPZ category.
Figure 10

AHP-based GWPZ map.

Figure 10

AHP-based GWPZ map.

Close modal
Figure 11

GWPZs' areal distribution.

Figure 11

GWPZs' areal distribution.

Close modal

Fuzzy-AHP

AHP and the related MCDM approach analyze difficult decisions using numerical and pragmatic principles (Saaty 1990). While AHP is commonly used in MCDM, the Fuzzy-AHP methodology, introduced by Van-Laarhoven and Pedrycz in 1983, allows for a subjective assessment of weighing multiple criteria to address discrepancies that arise during pairwise comparisons (Chang 1996; Tan et al. 2013; Sener et al. 2018). By employing the Fuzzy-AHP technique, factors can be compared and their relative importance can be determined, facilitating the management of challenging decisions. The integration of Fuzzy-AHP into multi-criteria decisions enhances the accuracy of groundwater potential mapping, as suggested by Sener et al. (2018). Fuzzy logic, a branch of AI, deals with approximate reasoning rather than precise calculations, providing a mathematical framework to handle uncertainty and vagueness by acknowledging that things can be partially true or false and allowing degrees of truth to be expressed mathematically. In this study, seven different theme levels were considered. Experts assigned criterion weights using Saaty's scale (1–9) to compare the criteria through a questionnaire. In the context of groundwater potential, a parameter with a low weight has minimal impact, while one with a high weight exerts significant influence. The experts' comments will be taken into account to review the pairwise comparisons and arrive at a conclusion. For example, in AHP, equal importance is denoted as ‘1,’ whereas in Fuzzy-AHP, it is represented as (1,1,1) for the lower, middle, and upper triangle values (Table 7).

Table 7

Scale of Saaty is fuzzy

DefinitionWeightsFuzzy-AHP weights
Equal importance (1,1,1) 
Moderate (2,3,4) 
Strong (4,5,6) 
Very strong (6,7,8) 
Extreme importance (9,9,9) 
Intermediate values 2;4;6;8 (1,2,3); (3,4,5); (5,6,7); (7,8,9) 
DefinitionWeightsFuzzy-AHP weights
Equal importance (1,1,1) 
Moderate (2,3,4) 
Strong (4,5,6) 
Very strong (6,7,8) 
Extreme importance (9,9,9) 
Intermediate values 2;4;6;8 (1,2,3); (3,4,5); (5,6,7); (7,8,9) 

In Table 2 AHP pairwise comparison matrix, values 1, 2, 3, and 5 and 6 are converted into fuzzy numbers (using ‘1’ as an example of 1,1,1). Half, third, fifth, and sixth reciprocals are not directly translated into fuzzy numbers. To find the reciprocal fuzzy number, use Equation (3),
(3)
where l denotes lower, m denotes middle, and u denotes upper.
The reciprocal value 1/2, which is written as in fourth row and second column of Table 2, is changed into a fuzzy number using Equation (3):

Similar to this, Table 8 shows the residual reciprocal values after they have been converted into fuzzy numbers.

Table 8

Pairwise comparison matrix that has been fuzzy

MatrixGeomorphologyDrainage densityRainfallElevationSlopeLULCSoil
Geomorphology (1,1,1) (1,1,1) (1,2,3) (2,3,4) (2,3,4) (2,3,4) (2,3,4) 
Drainage density (1,1,1) (1,1,1) (1,2,3) (2,3,4) (4,5,6) (4,5,6) (4,5,6) 
Rainfall (1/3, 1/2, 1/1) (1/3, 1/2, 1/1) (1,1,1) (2,3,4) (2,3,4) (2,3,4) (4,5,6) 
Elevation (1/4, 1/3, 1/2) (1/4, 1/3, 1/2) (1/4, 1/3, 1/2) (1,1,1) (1,1,1) (4,5,6) (4,5,6) 
Slope (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/4, 1/3, 1/2) (1,1,1) (1,1,1) (3,4,5) (3,4,5) 
LULC (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/5, 1/4, 1/3) (1,1,1) (2,3,4) 
Soil (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/6, 1/5, 1/4) (1/6, 1/5, 1/4) (1/5, 1/4, 1/3) (1/4, 1/3, 1/2) (1,1,1) 
MatrixGeomorphologyDrainage densityRainfallElevationSlopeLULCSoil
Geomorphology (1,1,1) (1,1,1) (1,2,3) (2,3,4) (2,3,4) (2,3,4) (2,3,4) 
Drainage density (1,1,1) (1,1,1) (1,2,3) (2,3,4) (4,5,6) (4,5,6) (4,5,6) 
Rainfall (1/3, 1/2, 1/1) (1/3, 1/2, 1/1) (1,1,1) (2,3,4) (2,3,4) (2,3,4) (4,5,6) 
Elevation (1/4, 1/3, 1/2) (1/4, 1/3, 1/2) (1/4, 1/3, 1/2) (1,1,1) (1,1,1) (4,5,6) (4,5,6) 
Slope (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/4, 1/3, 1/2) (1,1,1) (1,1,1) (3,4,5) (3,4,5) 
LULC (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/5, 1/4, 1/3) (1,1,1) (2,3,4) 
Soil (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/6, 1/5, 1/4) (1/6, 1/5, 1/4) (1/5, 1/4, 1/3) (1/4, 1/3, 1/2) (1,1,1) 

The fuzzy geometric mean and fuzzy weights of each criterion were obtained using Buckley (1985) geometric mean technique. The next step, as shown in the following, is solving Equation (4) to obtain the fuzzy geometric mean value, ri:
(4)
The fuzzy geometric mean value in Table 9 second row and ninth column is calculated as shown in the following:
Table 9

Fuzzy geometric mean values

MatrixGeomorphologyDrainage densityRainfallElevationSlopeLULCSoilFuzzy geometric mean value (ri)
Geomorphology (1,1,1) (1,1,1) (1,2,3) (2,3,4) (2,3,4) (2,3,4) (2,3,4) (1.49, 2.07, 2.58) 
Drainage density (1,1,1) (1,1,1) (1,2,3) (2,3,4) (4,5,6) (4,5,6) (4,5,6) (2.0, 2.57, 3.07) 
Rainfall (1/3, 1/2, 1/1) (1/3, 1/2, 1/1) (1,1,1) (2,3,4) (2,3,4) (2,3,4) (4,5,6) (1.19,1.65,2.34) 
Elevation (1/4, 1/3, 1/2) (1/4, 1/3, 1/2) (1/4, 1/3, 1/2) (1,1,1) (1,1,1) (4,5,6) (4,5,6) (0.82, 0.98,1.24) 
Slope (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/4, 1/3, 1/2) (1,1,1) (1,1,1) (3,4,5) (3,4,5) (0.71,0.86,1.07) 
LULC (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/5, 1/4, 1/3) (1,1,1) (2,3,4) (0.35,0.44,0.57) 
Soil (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/6, 1/5, 1/4) (1/6, 1/5, 1/4) (1/5, 1/4, 1/3) (1/4, 1/3, 1/2) (1,1,1) (0.24,0.30,0.39) 
MatrixGeomorphologyDrainage densityRainfallElevationSlopeLULCSoilFuzzy geometric mean value (ri)
Geomorphology (1,1,1) (1,1,1) (1,2,3) (2,3,4) (2,3,4) (2,3,4) (2,3,4) (1.49, 2.07, 2.58) 
Drainage density (1,1,1) (1,1,1) (1,2,3) (2,3,4) (4,5,6) (4,5,6) (4,5,6) (2.0, 2.57, 3.07) 
Rainfall (1/3, 1/2, 1/1) (1/3, 1/2, 1/1) (1,1,1) (2,3,4) (2,3,4) (2,3,4) (4,5,6) (1.19,1.65,2.34) 
Elevation (1/4, 1/3, 1/2) (1/4, 1/3, 1/2) (1/4, 1/3, 1/2) (1,1,1) (1,1,1) (4,5,6) (4,5,6) (0.82, 0.98,1.24) 
Slope (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/4, 1/3, 1/2) (1,1,1) (1,1,1) (3,4,5) (3,4,5) (0.71,0.86,1.07) 
LULC (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/5, 1/4, 1/3) (1,1,1) (2,3,4) (0.35,0.44,0.57) 
Soil (1/4, 1/3, 1/2) (1/6, 1/5, 1/4) (1/6, 1/5, 1/4) (1/6, 1/5, 1/4) (1/5, 1/4, 1/3) (1/4, 1/3, 1/2) (1,1,1) (0.24,0.30,0.39) 

Table 9 displays the final values for the fuzzy geometric means in a manner similar to this.

The next step is to generate fuzzy weights for each criterion. Equation (5) is used to determine fuzzy weights:
(5)
All fuzzy geometric mean values must first be added. In the following Equation (6) is given for adding fuzzy geometric values:
(6)
To obtain the total of fuzzy numbers, the lower values, middle values, and upper values are put together as illustrated in the following example from Table 9:
As demonstrated in Equation (5), the estimated value (6.8, 8.87, and 11.26) above should be reciprocated. The following computation is made using reciprocal Equation (3), which is then applied,:

As indicated in Tables 10 and 11, the subsequent step entails multiplying each fuzzy geometric mean value by the reciprocal of the geometric mean summation.

Table 10

Fuzzy weights and a fuzzy geometric mean value

ParametersFuzzy geometric mean value (ri)Fuzzy weights wi
Geomorphology (1.49, 2.07, 2.58) (1.49, 2.07, 2.58) ⊗ (1/11.26, 1/8.87, 1/6.8) 
Drainage density (2.0, 2.57, 3.07) (2.0, 2.57, 3.07) ⊗ (1/11.26, 1/8.87, 1/6.8) 
Rainfall (1.19,1.65,2.34) (1.19,1.65,2.34) ⊗ (1/11.26, 1/8.87, 1/6.8) 
Elevation (0.82, 0.98,1.24) (0.82, 0.98,1.24) ⊗ (1/11.26, 1/8.87, 1/6.8) 
Slope (0.71,0.86,1.07) (0.71,0.86,1.07) ⊗ (1/11.26, 1/8.87, 1/6.8) 
LULC (0.35,0.44,0.57) (0.35,0.44,0.57) ⊗ (1/11.26, 1/8.87, 1/6.8) 
Soil (0.24,0.30,0.39) (0.24,0.30,0.39) ⊗ (1/11.26, 1/8.87, 1/6.8) 
ParametersFuzzy geometric mean value (ri)Fuzzy weights wi
Geomorphology (1.49, 2.07, 2.58) (1.49, 2.07, 2.58) ⊗ (1/11.26, 1/8.87, 1/6.8) 
Drainage density (2.0, 2.57, 3.07) (2.0, 2.57, 3.07) ⊗ (1/11.26, 1/8.87, 1/6.8) 
Rainfall (1.19,1.65,2.34) (1.19,1.65,2.34) ⊗ (1/11.26, 1/8.87, 1/6.8) 
Elevation (0.82, 0.98,1.24) (0.82, 0.98,1.24) ⊗ (1/11.26, 1/8.87, 1/6.8) 
Slope (0.71,0.86,1.07) (0.71,0.86,1.07) ⊗ (1/11.26, 1/8.87, 1/6.8) 
LULC (0.35,0.44,0.57) (0.35,0.44,0.57) ⊗ (1/11.26, 1/8.87, 1/6.8) 
Soil (0.24,0.30,0.39) (0.24,0.30,0.39) ⊗ (1/11.26, 1/8.87, 1/6.8) 
Table 11

Fuzzy weights wi

ParametersFuzzy weights wi
Geomorphology (0.13,0.23,0.38) 
Drainage density (0.18,0.29, 0.45) 
Rainfall (0.10,0.18,0.34) 
Elevation (0.07, 0.11, 0.18) 
Slope (0.06, 0.09,0.16) 
LULC 0.03, 0.05,0.08) 
Soil (0.02, 0.03, 0.06) 
ParametersFuzzy weights wi
Geomorphology (0.13,0.23,0.38) 
Drainage density (0.18,0.29, 0.45) 
Rainfall (0.10,0.18,0.34) 
Elevation (0.07, 0.11, 0.18) 
Slope (0.06, 0.09,0.16) 
LULC 0.03, 0.05,0.08) 
Soil (0.02, 0.03, 0.06) 

Defuzzification, the subsequent step, converts a fuzzy number into a precise numerical value. Using Equation (7) and as indicated in Table 12, the defuzzification process is carried out:
(7)
Table 12

Defuzzification

ParametersFuzzy weights, wiWeights
Geomorphology (0.13,0.23,0.38) 0.25 
Drainage density (0.18,0.29, 0.45) 0.31 
Rainfall (0.10,0.18,0.34) 0.21 
Elevation (0.07, 0.11, 0.18) 0.12 
Slope (0.06, 0.09,0.16) 0.10 
LULC 0.03, 0.05,0.08) 0.05 
Soil (0.02, 0.03, 0.06) 0.04 
Total  =(0.25 + 0.31 + 0.21 + 0.12 + 0.10 + 0.05 + 0.04)
=1.07 
ParametersFuzzy weights, wiWeights
Geomorphology (0.13,0.23,0.38) 0.25 
Drainage density (0.18,0.29, 0.45) 0.31 
Rainfall (0.10,0.18,0.34) 0.21 
Elevation (0.07, 0.11, 0.18) 0.12 
Slope (0.06, 0.09,0.16) 0.10 
LULC 0.03, 0.05,0.08) 0.05 
Soil (0.02, 0.03, 0.06) 0.04 
Total  =(0.25 + 0.31 + 0.21 + 0.12 + 0.10 + 0.05 + 0.04)
=1.07 

Once weights have been obtained, all weights must then be added together as illustrated in Table 12. 1.07 is an unacceptable number for the sum. To obtain the weight total of 1 as stated in Table 13, the weights are normalized.

Table 13

Normalized weights defuzzification

ParametersWeightsNormalized weights
Geomorphology 0.25 0.25/1.07 = 0.23 
Drainage density 0.31 0.31/1.07 = 0.29 
Rainfall 0.21 0.21/1.07 = 0.19 
Elevation 0.12 0.12/1.07 = 0.11 
Slope 0.10 0.10/1.07 = 0.10 
LULC 0.05 0.05/1.07 = 0.05 
Soil 0.04 0.04/1.07 = 0.03 
Total 1.07 1.00 
ParametersWeightsNormalized weights
Geomorphology 0.25 0.25/1.07 = 0.23 
Drainage density 0.31 0.31/1.07 = 0.29 
Rainfall 0.21 0.21/1.07 = 0.19 
Elevation 0.12 0.12/1.07 = 0.11 
Slope 0.10 0.10/1.07 = 0.10 
LULC 0.05 0.05/1.07 = 0.05 
Soil 0.04 0.04/1.07 = 0.03 
Total 1.07 1.00 

The normalized weights, which sum up to 1, were applied to each thematic layer included in the research. These weights indicate the relative importance of each layer in determining the GWPZ. In this present study, drainage density had the highest influence with a weight of 29%. Geomorphology was the second most important factor with a weight of 23%, followed by rainfall with a weight of 19%. On the other hand, elevation, slope, LULC, and soil were considered to have the least impact on defining groundwater recharge areas, with weights of 11, 10, 5, and 3%, respectively. This suggests that these characteristics have a relatively lower influence compared with other factors. Using the weighted overlay analysis in ArcGIS, the thematic layers representing the influential factors were overlaid to create the final GWPZs map. Figure 12 illustrates the resulting map, displaying the GWPZs in the study area. Additionally, three zones – poor, moderate, and high – were added to this map to indicate different levels of groundwater potential. Upon analysis with the Fuzzy-AHP model, Figure 13 shows that the present study area was classified as 18.27% poor potential, 71.07% moderate potential, and 10.69% high potential. This indicates that a significant portion of the study area falls within the moderate GWPZ.
Figure 12

Fuzzy-AHP-based GWPZ map.

Figure 12

Fuzzy-AHP-based GWPZ map.

Close modal
Figure 13

Percentage of the areal distribution.

Figure 13

Percentage of the areal distribution.

Close modal

Validation of groundwater potential map

Validation is crucial in scientific research. Data from 30 actual groundwater wells were used to confirm the final GWPZ obtained from the open data Telangana portal (https://data.telangana.gov.in) (Febi et al. 2020). With the use of the already-existing well data points, the study's validity was checked, and the created groundwater potential map was verified. The GWPZ map of the research region was created by combining overlay investigations of the selected parameters given above with the GIS-based AHP and Fuzzy-AHP methods. A common technique for assessing the accuracy of a model outcome is the ROC (Janizadeh et al. 2019). The AUC in ROC analysis represents the suggested model's level of accuracy. The correlation between the AUC value and the prediction accuracy can be classified into five categories based on their qualitative and quantitative relationship. These categories are: excellent (0.9–1), very good (0.8–0.9), good (0.7–0.8), average (0.6–0.7), and poor (0.5–0.6) (Yesilnacar & Topal 2005). For example, Thanh et al. (2022) conducted a validation study of an AHP-based model for GWPZs by comparing it with other models. The study found that the AHP-based model achieved an AUC value of 0.72, which indicates a relatively lower level of accuracy in predicting GWPZs compared with the alternative models. On the other side, El-Magd & Eldosouky (2021) conducted a validation study where they compared an AHP-based model with other models for GWPZs. Their results showed that the AHP model achieved an AUC value of 0.72. The present study supports the findings reported by Rahmati et al. (2015) in their research on GIS-AHP coupled groundwater potential mapping in Iran. Validation of the results showed that the AHP method exhibited fairly good prediction accuracy, with an AUC value of 73.66%. Additionally, the research conducted by Das (2019) on groundwater potential zonation in the Vaitarna basin involved the validation of an AHP-based model against other models for GWPZs. The results of the study showed that the AHP model obtained an AUC value of 0.70, indicating relatively less accuracy compared with the alternative models. The GWPZs map employing AHP had an AUC of 0.73, which means the accuracy of the model is 73% (Figure 14). On the other side, the GWPZs map employing Fuzzy-AHP had an AUC of 0.74, which means the accuracy of the model is 74% (Figure 15). The model's accuracy indicates that both AHP and Fuzzy-AHP results are good, but Fuzzy-AHP is more accurate in assessing groundwater potential.
Figure 14

Validation for GWPZs (AHP).

Figure 14

Validation for GWPZs (AHP).

Close modal
Figure 15

Validation for GWPZs (Fuzzy-AHP).

Figure 15

Validation for GWPZs (Fuzzy-AHP).

Close modal

Suitable artificial recharge sites

The practice of artificial aquifer recharge has gained significance in addressing groundwater degradation and depletion issues in various river basins worldwide (Aju et al. 2021). Artificial recharge involves the deliberate replenishment of groundwater by redirecting excess surface water into underground layers for storage and future use. This process helps replenish the water supply within aquifers. Studies have demonstrated the effectiveness of artificial recharge systems in preserving and replenishing groundwater resources. Capturing and conserving rainwater, which would otherwise seep into the ground or flow through surface streams, is an essential aspect of artificial recharge (Kumar et al. 2021).

Artificial recharge methods, such as spreading, are employed to raise the water table and enhance groundwater levels. The utilization of RS and GIS techniques plays a crucial role in identifying suitable locations for water recharge. Researchers worldwide have successfully employed GIS to delineate suitable areas for various artificial recharge technologies (Ahmadi et al. 2017). By combining scientific methods, including remote sensing, GIS, and artificial recharge techniques, sustainable management of groundwater resources can be achieved. These approaches contribute to the preservation and replenishment of aquifers, ensuring the availability of water for future generations.

In order to identify suitable sites for groundwater recharge, the Fuzzy-AHP method was employed and validated as the best approach. Figure 16 illustrates the resulting map obtained from overlaying the lineament and drainage maps, which facilitated the identification of appropriate locations for artificial recharge. A total of 26 artificial recharge sites were identified, with 2 falling within the high recharge zone, 22 within the moderate zone, and 2 within the poor zone as determined by the Fuzzy-AHP analysis.
Figure 16

Artificial recharge sites of the Kinnerasani Watershed.

Figure 16

Artificial recharge sites of the Kinnerasani Watershed.

Close modal

For the implementation of artificial recharge, it is recommended to utilize gabion structures as the primary facilities, particularly in lower order streams. These structures enhance the infiltration of water into the subsurface formations, aiding in the recharge process. It is noteworthy that a significant number of identified artificial recharge sites are located in agricultural land, which holds promising implications for fostering higher economic growth. By strategically selecting and implementing artificial recharge sites based on the findings of this study, decision-makers and water users in the research area can actively contribute to sustainable groundwater management. This approach not only improves groundwater recharge potential but also promotes economic development and ecological sustainability within the Kinnerasani Watershed.

The purpose of the current work is to define GWPZs in the Kinnerasani Watershed in India using integrated data, including remote sensing and GIS, in conjunction with the AHP and Fuzzy-AHP techniques. This study utilized seven thematic layers, namely drainage density, slope, elevation, LULC, soil, geomorphology, and rainfall, which were analyzed to delineate the GWPZs. The present study region was divided into three unique areas, each characterized by a different level of groundwater potential, as depicted in the final GWPZ map. Among the factors considered, drainage density emerged as a significant influence on groundwater potential in both the AHP and Fuzzy-AHP models. According to the AHP model, the study area was classified as 17.76% poor potential, 72.79% moderate potential, and 9.45% high potential. In comparison, the Fuzzy-AHP approach classified the area as 18.27% poor potential, 71.07% moderate potential, and 10.69% high potential. Evaluation of predictive capability through ROC calculations revealed that the Fuzzy-AHP technique exhibited better accuracy (AUC = 74%) than the AHP technique (AUC = 73%).

The GWPZ map generated by the Fuzzy-AHP technique serves as a valuable database for groundwater planning and management in the study area. Additionally, this study utilizes the best method derived from both models to identify 26 suitable locations for artificial recharge sites. Given the extensive agricultural practices in the study area, the findings hold considerable potential for improving irrigation facilities, agricultural output, and domestic water use productivity. To maintain sustainable groundwater use in the area, decision-makers can leverage the research findings to develop effective groundwater management strategies and smart development plans. By incorporating the integrated approach and methodologies employed in this study, informed decisions can be made to ensure optimal utilization of groundwater resources in the Kinnerasani Watershed.

Future studies on GWPZs can focus on exploring groundwater–surface water interactions, improving data integration, incorporating socio-economic factors, considering the impacts of climate change, and implementing long-term monitoring and modelling.

The authors would like to thank Telangana Ground Water Department Water Level Data for making public available (https://data.telangana.gov.in) as well as the US Geological Survey (USGS) for making the satellite data available (https://earthexplorer.usgs.gov/). Thank you to the Climatic Research Unit for providing rainfall data on their website at https://crudata.uea.ac.uk/cru/data/hrg/. We would especially like to thank the Geological Survey of India for providing the geomorphology and lineament data on their website at https://bhukosh.gsi.gov.in/Bhukosh/Public.

There was no funding for this project.

All relevant data are available from an online repository or repositories.

The authors declare there is no conflict.

Ahmed
N.
,
Hoque
M. A. A.
,
Pradhan
B.
&
Arabameri
A.
2021
Spatio-temporal assessment of groundwater potential zone in the drought-prone area of Bangladesh using GIS-based bivariate models
.
Nat. Resour. Res.
30
,
3315
3337
. https://doi.org/10.1007/s11053-021-09870-0.
Aju
C. D.
,
Achu
A. L.
,
Raicy
M. C.
&
Reghunath
R.
2021
Identification of suitable sites and structures for artificial groundwater recharge for sustainable water resources management in Vamanapuram River Basin, South India
.
Hydro. Res.
24
37
. https://doi.org/10.1016/j.hydres.2021.04.001.
Allafta
H.
,
Opp
C.
&
Patra
S.
2021
Identification of groundwater potential zones using remote sensing and GIS techniques: a case study of the Shatt al-Arab Basin
.
Remote Sens.
13
,
1
28
. https://doi.org/10.3390/rs13010112.
Andualem
T. G.
&
Demeke
G. G.
2019
Groundwater potential assessment using GIS and remote sensing: a case study of Guna tana landscape, upper Blue Nile Basin, Ethiopia
.
J. Hydrol.: Reg. Stud.
24
,
100610
.
Arunbose
S.
,
Srinivas
Y.
,
Rajkumar
S.
,
Nair
N. C.
&
Kaliraj
S.
2021
Remote sensing, GIS and AHP techniques based investigation of groundwater potential zones in the Karumeniyar river basin, Tamil Nadu, southern India
.
Groundwater Sustainable Dev.
14
,
100586
. https://doi.org/10.1016/J.GSD.2021.100586.
Bezdek
J. C.
,
Keller
J.
,
Krishnapuram
R.
&
Pal
N. R.
1999
Fuzzy Models and Algorithms for Pattern Recognition and Image Processing. The Handbooks of Fuzzy Sets Series
.
Kluwer Academic
,
Dordrecht
.
Buckley
J. J.
1985
Fuzzy hierarchical analysis
.
Fuzzy Set Syst.
17
,
233
247
. https://doi.org/10.1016/0165-0114(85)90090-9.
Chang
D. Y.
1996
Applications of the extent analysis method on Fuzzy AHP
.
Eur. J. Oper. Res.
95
(
3
),
649
655
.
Chang
N. B.
,
Parvathinathan
G.
&
Breeden
J. B.
2008
Combining GIS with fuzzy multicriteria decision-making for landfill siting in a fastgrowing urban region
.
J. Environ. Manage.
87
(
1
),
139
153
.
Das
B.
,
Pal
S. C.
,
Malik
S.
&
Chakrabortty
R.
2019
Modeling groundwater potential zones of Puruliya district, West Bengal, India using remote sensing and GIS techniques
.
Geol. Ecol. Landscapes
3
,
223
237
. https://doi.org/10.1080/ 24749508.2018.1555740.
Das
M.
,
Parveen
T.
,
Ghosh
D.
&
Alam
J.
2021
Assessing groundwater status and human perception in drought-prone areas: a case of Bankura-I and Bankura-II blocks, West Bengal (India)
.
Environ. Earth Sci.
80
(
18
). https://doi. org/10.1007/s12665-021-09909-8.
Delin
G. N.
,
Healy
R. W.
,
Landon
M. K.
&
Bohlke
J. K.
2000
Effects of topography and soil properties on recharge at two sites in an agricultural field
.
J. Am. Water Resour. Assoc.
36
,
1401
1416
. https://doi.org/10.1111/j.1752-1688.2000.tb05735.x.
Dibs
H.
,
Ali
A. H.
,
Al-Ansari
N.
&
Abed
S. A.
2023
Fusion Landsat-8 thermal TIRS and OLI datasets for superior monitoring and change detection using remote sensing
.
Emerging Sci. J.
7
(
2
).
doi:10.28991/ESJ-2023-07-02-09
.
Donselaar
M. E.
,
Bhatt
A. G.
&
Ghosh
A. K.
2017
On the relation between fluvio-deltaic flood basin geomorphology and the wide-spread occurrence of arsenic pollution in shallow aquifers
.
Sci. Total Environ.
574
,
901
913
. https://doi.org/10.1016/j.scitotenv.2016.09.074.
Dubois
D.
&
Prade
H.
2003
Fuzzy set and possibility theory-based methods in artificial intelligence
.
Artif. Intell.
148
(
1–2
),
1
9
. https://doi.org/10.1016/S0004-3702(03)00118-8.
El-Magd
S. A. A.
&
Eldosouky
A. H.
2021
An improved approach for predicting the groundwater potentiality in the low desert lands; El-Marashda area, Northwest Qena City, Egypt
.
J. Afr. Earth. Sci.
179
,
104200
. https://doi.org/10.1016/j.jafrearsci.2021.104200.
FAO
1988
Soil Map of the World, Revised Legend with Corrections and Updates. Report of the World Soil Resource, Food and Agricultural Organization. FAO, Rome, Italy, Report 60.
Febi
B.
,
Achu
A. L.
,
Jiminisha
K.
,
Aiysha
V. A.
&
Raghunath
R.
2020
Landslide susceptibility modeling using integrated evidential belief function based logistic regression method: a study from Sothern Western Ghat, India
.
J. RSASE
20
,
100411
.
Fitts
C. R.
2013
Groundwater Science
.
Elsevier
. https://doi.org/10.1016/C2009-0- 62950-0.
Flug
B. M.
,
Seitz
H. L. H.
,
Member
A.
&
Scott
J. F.
2005
Multicriteria decision analysis applied to Glen Canyon Dam
.
J. Water Resources Plann. Manage.
126
(
October
),
270
276
.
Goswami
T.
&
Ghosal
S.
2022
Understanding the suitability of two MCDM techniques in mapping the groundwater potential zones of semi-arid Bankura District in Eastern India
.
Groundwater Sustainable Dev.
17
. https://doi.org/10.1016/j. gsd.2022.100727 100727.
Greenbaum
D.
1985
Review of Remote Sensing Applications to Groundwater Exploration in Basement and Regolith
.
British Geological Survey
,
Nicker Hill
.
Gumma
M. K.
&
Pavelic
P.
2013
Mapping of groundwater potential zones across Ghana using remote sensing, geographic information systems, and spatial modeling
.
Environ. Monit. Assess.
185
,
3561
3579
. https://doi.org/10.1007/ s10661-012-2810-y.
Horton
R. E.
1932
Drainage-basin characteristics
.
Trans. Am. Geophys. Union
13
,
350
361
.
Ijumulana
J.
,
Ligate
F.
,
Irunde
R.
,
Bhattacharya
P.
,
Maity
J. P.
,
Ahmad
A.
&
Mtalo
F.
2021
Spatial uncertainties in fluoride levels and health risks in endemic fluorotic regions of northern Tanzania
.
Groundwater Sustainable Dev.
14
,
100618
.
Janizadeh
S.
,
Avand
M.
,
Jaafari
A.
,
Phong
T. V.
,
Bayat
M.
,
Ahmadisharaf
E.
,
Prakash
I.
,
Pham
B. T.
&
Lee
S.
2019
Prediction success of machine learning methods for flash flood susceptibility mapping in the Tafresh watershed, Iran
.
Sustainability
11
(
19
),
5426
. https://doi.org/10.3390/su11195426.
Kumar
A.
,
Kanga
S.
,
Taloor
A. K.
,
Singh
S. K.
&
Durin
B.
2021
Surface runoff estimation of Sind river basin using integrated SCS-CN and GIS techniques
.
Hydro. Res.
4
,
61
74
. ISSN 2589-7578, https://doi.org/10.1016/j.hydres.2021.08.001.
Kumar
M.
,
Singh
P.
&
Singh
P.
2022
Integrating GIS and remote sensing for delineation of groundwater potential zones in Bundelkhand Region, India, Egypt
.
J. Remote Sens. Space. Sci.
25
,
387
404
. https://doi.org/10.1016/J. EJRS.2022.03.003.
Lyu
H. M.
,
Sun
W. J.
,
Shen
S. L.
&
Zhou
A. N.
2020
Risk assessment using a new consulting process in fuzzy AHP
.
J. Constr. Eng. Manage.
146
,
04019112
.
Mandal
U.
,
Sahoo
S.
,
Munusamy
S. B.
,
Dhar
A.
,
Panda
S. N.
,
Kar
A.
&
Mishra
P. K.
2016
Delineation of groundwater potential zones of coastal groundwater basin using multi-criteria decision-making technique
.
Water Resour. Manage.
30
,
4293
4310
. https://doi.org/10.1007/s11269-016-1421-8.
Martin
S. L.
,
Hayes
D. B.
,
Kendall
A. D.
&
Hyndman
D. W.
2017
The land-use legacy effect: towards a mechanistic understanding of time-lagged water quality responses to land use/cover
.
Sci. Total Environ.
579
,
1794
1803
. https://doi.org/10.1016/j.scitotenv.2016.11.158.
Mohanty
W. K.
&
Walling
M. Y.
2008
First order seismic microzonation of Haldia, Bengal Basin (India) using a GIS platform
.
Pure Appl. Geophys.
165
(
7
),
1325
1350
.
Nguyen
T. G.
,
Phan
K. A.
&
Huynh
T. H. N.
2022
Application of integrated-weight water quality index in groundwater quality evaluation
.
Civ. Eng. J.
8
,
11
.
doi:10.28991/CEJ-2022-08-11-020
.
Nyakundi
R.
,
Nyadawa
M.
&
Mwangi
J.
2022
Effect of recharge and abstraction on groundwater levels
.
Civ. Eng. J.
8
,
05
.
doi:10.28991/CEJ-2022-08-05-05
.
Omar
P. J.
,
Dwivedi
S. B.
,
Dikshit
P. K. S.
,
2020
Sustainable development and management of groundwater in Varanasi, India
. In:
Advances in Water Resources Engineering and Management
(
AlKhaddar
R.
,
Singh
R. K.
,
Dutta
S.
&
Kumari
M.
, eds).
Springer Singapore
,
Singapore
, pp.
201
209
.
Omar
P. J.
,
Gaur
S.
,
Dwivedi
S. B.
,
Dikshit
P. K. S.
,
2021
Development and application of the integrated GIS-MODFLOW model
. In:
Fate and Transport of Subsurface Pollutants
(
Gupta
P. K.
&
Bharagava
R. N.
, eds).
Springer Singapore
,
Singapore
, pp.
305
314
. https://doi.org/10.1007/978-981-15-6564-9_16.
Qadir
J.
,
Bhat
M. S.
,
Alam
A.
&
Rashid
I.
2020
Mapping groundwater potential zones using remote sensing and GIS approach in Jammu Himalaya, Jammu and Kashmir
.
GeoJournal
85
,
487
504
. https://doi.org/10.1007/S10708-019-09981- 5/TABLES/10.
Quadrio
B.
,
Ambrosanio
M.
&
Ioannilli
M.
2015
A GIS-based tool for reliability assessment of seismic microzonation studies according to Italian instructions and criteria. In: Engineering Geology for Society and Territory, Vol. 5. Springer, pp. 1123–1126.
Rahmati
O.
,
Samani
A. N.
,
Mahdavi
M.
,
Pourghasemi
H. R.
&
Zeinivand
H.
2015
Groundwater potential mapping at Kurdistan region of Iran using analytic hierarchy process and GIS
.
Arabian J. Geosci.
8
(
9
),
7059
7071
.
Roberts
R.
&
Goodwin
P.
2002
Weight approximations in multi-attribute decision models
.
J. Multi-Criteria Decis. Anal.
11
(
6
),
291
303
. https://doi.org/10.1002/ mcda.320.
Saade
J.
,
Atieh
M.
,
Ghanimeh
S.
&
Golmohammadi
G.
2021
Modeling impact of climate change on surface water availability using SWAT model in a semi-arid basin: case of El Kalb River, Lebanon
.
Hydrology
8
(
3
),
134
. https://doi.org/10.3390/hydrology8030134.
Saaty
T. L.
1980
The Analytic Hierarchy Process
.
McGraw-Hill
,
New York
, p.
287
Saaty
T. L.
1996
Decision Making with Dependence and Feedback, The Analytic Network Process
.
RWS Publications
,
Pittsburgh, PA
.
Saaty
T. L.
2008
Decision making with the analytic hierarchy process
.
Int. J. Serv. Sci.
1
(
1
),
83
98
.
Sener
E.
,
Sehnaz
S.
&
Aysen
D.
2018
Groundwater potential mapping by combining fuzzy analytic hierarchy process and GIS in Beysehir Lake Basin, Turkey
.
Arabian J. Geosci.
11
,
187
.
doi:10.1007/s12517-018-3510-x
.
Serele
C.
,
Pérez-Hoyos
A.
&
Kayitakire
F.
2020
Mapping of groundwater potential zones in the drought-prone areas of south Madagascar using geospatial techniques
.
Geosci. Front.
11
,
1403
1413
. https://doi.org/10.1016/j. gsf.2019.11.012.
Shit
P. K.
,
Paira
R.
,
Bhunia
G.
&
Maiti
R.
2015
Modeling of potential gully erosion hazard using geo-spatial technology at Garbheta block, West Bengal in India
.
Model. Earth Syst. Environ.
1
,
2
. https://doi.org/10.1007/s40808-015-0001-x.
Shyam
G. M.
,
Taloor
A. K.
,
Singh
S. K.
&
Kanga
S.
2021
Sustainable water management using rainfall-runoff modeling: a geospatial approach
.
Groundwater Sustainable Dev.
15
,
100676
.
Tan
R. R.
,
Aviso
K. B.
,
Huelgas
A. P.
&
Promentilla
M. A. B.
2013
Fuzzy AHP approach to selection problems in process engineering involving quantitative and qualitative aspects
.
Process Saf. Environ. Prot.
doi:10.1016/j.psep.2013.11.005
.
Tesfaye
T.
2010
Ground water potential evaluation based on integrated GIS and RS techniques in Bilate river catchment, South rift valley of Ethiopia. Am. Sci. Res. J. Eng. Technol. Sci.. ISSN (Print) 2313-4410, ISSN (Online) 2313-4402. Global Society of Scientific Research and Researchers. Available from: http://asrjetsjournal.org
Thacker
H.
,
Shah
Y.
,
Borah
A. J.
,
Jadeja
Y.
,
Thakkar
M.
,
Bhimani
S.
&
Chauhan
G.
2022
Assessment of groundwater potential zones across Katrol Hill Fault, Kachchh, Western India: a remote sensing and GIS approach
.
Open J. Geol.
12
(
2
),
111
135
.
Thanh
N. N.
,
Chotpantarat
S.
,
Trung
N. H.
,
Ngu
N. H.
&
Muoi
L. N.
2022
Mapping groundwater potential zones in Kanchanaburi Province, Thailand by integrating of analytic hierarchy process, frequency ratio, and random forest
.
Ecol. Indic.
145
,
109591
. https://doi.org/10.1016/j.ecolind.2022.109591.
UNESCO
2015
The United Nations World Water Development Report 2015: Water for A Sustainable World
.
UNESCO Publishing
,
Paris
.
Van-Laarhoven
P. J. M.
&
Pedrycz
W.
1983
A fuzzy extension of Saaty's priority theory
.
Fuzzy Sets Syst.
11
,
229
241
.
Zade
N.
,
Avadhoot
B.
,
Kumar
D. P.
,
Pradip
S.
&
Robin
D.
2021
Variability of mechanical properties of cellular lightweight concrete infill and its effect on seismic safety
.
Nat. Hazard. Rev.
22
,
4021039
. https://doi.org/10.1061/(ASCE)NH.1527-6996.0000501.
Zolekar
R. B.
&
Bhagat
V. S.
2015
Multi-criteria land suitability analysis for agriculture in hilly zone: remote sensing and GIS approach
.
Comput. Electron. Agric.
118
,
300
321
. https://doi.org/10.1016/j.compag.2015.09.016.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).