This study formulates a sustainable groundwater recharge plan for Kawardha District, Chhattisgarh, India, by integrating the analytical hierarchy process (AHP), remote sensing (RS), and Geographic Information System (GIS)-based multi-criteria decision analysis (MCDA). Eight thematic layers – geology, geomorphology, land use/land cover (LULC), soil, lineament density, drainage density, rainfall, and slope – were used to evaluate groundwater recharge potential. The district's elevation varies between 273 and 1,016 m, with slopes primarily below 2%, while drainage density peaks at 3.4 km/sq.km. Pre- and post-monsoon groundwater levels range from 1.6–8.1 to 0.6–7.2 m, respectively. Forests dominate the six classified LULC categories, and sandy loam is the primary soil type. AHP weighting prioritized geology (27%) and geomorphology (24%). Groundwater potential zones (GPZs) were categorized into five classes: very low (0.50%), low (17.60%), moderate (60.46%), high (20.91%), and very high (0.53%), with a 79% accuracy validated by well data. To address irregular rainfall and complex terrain, 64 check dams, 35 percolation tanks, and 20 farm ponds are recommended. This approach offers valuable insights for groundwater management in areas with similar hydrogeological issues, aiding water resource and agricultural planning.

  • The study integrates analytical hierarchy process, remote sensing, and GIS within an multi-criteria decision analysis framework for sustainable groundwater management.

  • Eight thematic maps were generated, including geology, geomorphology, and land use, to assess groundwater potential.

  • Groundwater potential zones were classified into five categories with an accuracy of 79%.

  • Proposed 64 check dams, 35 percolation tanks, and 20 farm ponds to enhance recharge.

Groundwater is a key source of fresh water and essential for the health and well-being of all living beings. It plays a vital role in supporting economic development and meeting human needs. However, the over-exploitation of groundwater has become a significant concern. To address this issue and prevent excessive extraction, it is important for authorities to consider implementing legislative measures. Major reasons for groundwater depletion are population growth, urbanization, industrialization, deforestation, and climate change, as well as mis-management of the existing water resources (Nassery et al. 2017; Baghel et al. 2023). The need for high-quality water has greatly expanded in recent years due to a growth in agricultural, industrial and domestic activities. Groundwater is the most widely used fresh water entity due to purity and easily availability (Dinka et al. 2015; Maghrebi et al. 2021). Poor management of fresh water resources could lead to significant problems in the future, especially if groundwater is prioritized for use by households, agriculture, and industrial sectors (Kuylenstierna et al. 1997). Global food security and economic growth highly depends on groundwater resources towards qualitative aspects and sustainable resources (Tortajada & González-Gómez 2022). In India, many regions are experiencing a decline in groundwater levels at a rate of approximately 1–2 m per year (Singh & Singh 2002). India is now the largest groundwater user in the world. About 70% of groundwater withdrawn globally is used in the agriculture sector and irrigation in the arid/semi-arid counties (Maroufpoor & Shiri 2020; Mukherjee et al. 2022). The significance of groundwater resources in India is underscored by the fact that approximately 87% is used for agriculture, while the remaining 11 and 2% are allocated to domestic and industrial purposes, respectively (Mukherjee 2018). India has vast geographical features and climatic variability along with uncertainty, proper management of the water resources is quite difficult due to interstate disputes. India is blessed with a sufficient amount of rainfall as well as snowmelt fresh water from the river system (Prasad & Roy 2005). Due to inadequate management of water resources in India, much of the runoff flows from upstream to downstream areas and is excessively discharged into the ocean (Jain et al. 2007). Adverse climatic conditions like drought also occurred in various places in India. Hence, groundwater is the only source of water for survival as well as food security (Shah & Mishra 2020).

Climatic variability, particularly in terms of rainfall intensity, plays a major role in determining infiltration rates and runoff generation from a catchment area. The rainfall intensity is high as compared to the infiltration rate, maximum probability of peak runoff in shorter time span and vice versa (Hou et al. 2019). The catchment's curve number is relatively high, which could contribute to increased runoff due to insufficient infiltration rates or groundwater recharge (Sahu et al. 2010). Further, various researchers have suggested the groundwater recharge plans have benefitted from the above-mentioned problems (Rajasekhar et al. 2019; Dar et al. 2021). Historical data highlight that India has experienced over exploitation of groundwater without conservation measures. In this context various cases of groundwater depletion were observed in water shortage regions in India. Henceforth, India will be on the edge of a major upcoming water calamity (Singh et al. 2018). The status of groundwater not only depends on single parameters, it is linked with number of spatio-temporal parameters, i.e., precipitation, snow-melt water, topography, land use and land cover patterns, physio-chemical properties of soil, existing water resources, climatic variability, geographical features of the region, etc. However, more or less most important parameter of groundwater contribution is precipitation and existing water resources of the catchment (Shah 2008). In sustainable planning of groundwater resources development, excess runoff management is the major criteria for storage structure inside the catchment area (Chenini et al. 2010). In recent past, most of the conventional methods are based on historical data sets which are adopted for groundwater recharge plans by making the small check dams, ponds, artificial lakes, farm ponds, etc. (Abraham & Mohan 2019). Hydrogeological, geological, and geophysical techniques have been used in a variety of field-based conventional methods to delineate GPZs; however, these methods are often point-based, costly, laborious and time consuming in spatiotemporal information (Doke et al. 2021). Presently, GIS- and RS-based techniques are most widely adopted for groundwater exploration and recharge planning and mitigation (Jat et al. 2009; Beden et al. 2023). The frequency ratio is among the improved methodological strategies for groundwater exploration that various scientists have been working on in recent years (Amponsah et al. 2023). The decision tree model (Stumpp et al. 2016), the maximum entropy model (Rahmati et al. 2016), the fuzzy logic model (Cameron & Peloso 2001), the Dempster-Shafer model (Roy & Datta 2019), the weights of evidence model (Masetti et al. 2007; Tahmassebipoor et al. 2016), the artificial neural network model (Mohanty et al. 2010; Chitsazan et al. 2015), and the logistic regression model have all been effectively utilized towards groundwater studies on a global level.

The AHP is one of the most widely used and respected decision-making methods among various MCDA techniques. Its popularity stems from its structural simplicity, effectiveness, minimal bias, clarity, efficiency, and its ability to consider multiple factors that influence a watershed (Srdjevic & Medeiros 2008; Jhariya et al. 2016; Akinlalu et al. 2021). Saaty created the AHP technique for decision-making in the early 1980s (Saaty 1980). It has a number of benefits and uses quantitative and qualitative techniques to address issues (Chandio et al. 2013). By maintaining natural surface water flow, GPZs can be used for artificial groundwater recharge, thereby enhancing groundwater storage. In regions where aquifers have been depleted due to excessive development, artificial recharge techniques are often employed to boost sustainable productivity (Gnanachandrasamy et al. 2018). The process of exploring groundwater is essentially a geophysical and hydrogeological operation that relies on accurate interpretation of hydrological indicators and evidence (Antony Ravindran 2012). If the necessary steps are not taken promptly, India may face a severe water crisis, with many states already experiencing significant water scarcity. Although the country receives an average of 1,200 mm of precipitation annually, which is more than many comparably sized countries, the distribution of this rainfall is uneven. Consequently, different regions can experience drought conditions despite the overall high levels of precipitation (Pathak & Dodamani 2019; Biswas et al. 2023). In many areas worldwide, groundwater levels have rapidly declined as a result of the monsoon's repeated failure and overuse. Because of the general decline in the groundwater table, there is also less water per person accessible. But significant precipitation is lost because of surface drainage into the ocean (Joshi & Tyagi 1994; Briscoe 2005; Kumar et al. 2005). Areas with higher drainage density promote increased surface runoff compared to those with lower drainage density. Surface water bodies, such as rivers and ponds, can act as recharge zones, enhancing the potential for nearby groundwater replenishment. Therefore, recognizing and measuring these features is critical in developing a groundwater potential model for a specific region. In India, about 42% of its area has been covered with rocky formation. Therefore, it is considered a crucial zone, and managerial measures are very difficult in these areas. Understanding the groundwater conditions completely requires a detailed hydrogeological examination (Saraf & Choudhury 1998).

Groundwater recharge structures have become an important aspect of water management practices in current times due to the increasing scarcity of water resources. The Kawardha district in Chhattisgarh is an area that has been grappling with significant water shortages, largely due to the over-exploitation of groundwater and inadequate water management practices. In this regard, the current study aims to suggest suitable groundwater recharge structures for the Kawardha district of Chhattisgarh. The study will focus on identifying the best possible locations for the implementation of recharge structures, considering the geological, hydrological and topographical characteristics of the region.

The study will explore the implementation of recharging wells, percolation tanks, check dams, and other effective structures to enhance groundwater recharge in the region. It will also recommend management practices to ensure the proper maintenance and long-term sustainability of these recharge structures. These practices will focus on preserving the integrity of the recharge systems and mitigating damage from both natural and human-induced factors. In summary, the study aims to provide valuable insights into the groundwater recharge situation in the Kawardha district of Chhattisgarh and recommend appropriate recharge structures for the area. This will not only address water scarcity issues but also support the sustainable management of groundwater resources in the region.

Salient features of study area

The study area is the Kawardha district in Chhattisgarh, India, located between latitudes 21°43′66″ N and 22°34′54″ N, and longitudes 80°46′2″ E and 81°35′50″ E (Figure 1). The district's topography is varied, with elevations ranging from 273 to 1,016 m above mean sea level, creating a significant elevation difference of 743 m. According to the digital elevation model (DEM), the district features steep slopes, particularly in the north-west, which contribute to high and rapid surface runoff in the downstream areas. The district spans approximately 4,447 km² and receives an average annual rainfall of 858.5 mm, with about 80% occurring during the monsoon season from June to September. Kawardha is bordered by Mungeli, Bemetara, and Rajnandgaon districts, and its western boundary connects with Madhya Pradesh. The drainage pattern directs runoff from northwest to southeast. Major land uses include forests, agriculture, and barren lands. The soil supports crops such as paddy, maize, sugarcane, and pigeonpea in the kharif season, and wheat, chickpea, soybean, and groundnuts in the rabi season. The region experiences a dry sub-humid climate with temperatures ranging from 15 °C in January to 43 °C in May.
Figure 1

Location of the study area.

Figure 1

Location of the study area.

Close modal

Data acquisition

Slope, geology, rainfall, drainage network, drainage density, soil, geomorphology, lineament, land use land cover (LULC), and relief were analysed to establish a GPZ for the study area. For the purpose, Shuttle Radar Topography Mission (SRTM) DEM was used for the delineation of a drainage network of the study area and validation of drainage pattern, survey of India toposheets have been utilized. Soil data have been acquired from the food and agriculture organization (FAO) site, lineament density map downloaded from the Bhuvan (ISRO), India website. Rainfall fine resolution (0.25° × 0.25°) gridded data have been downloaded from India Meteorological Department (IMD), Pune. Additionally, rainfall data were extracted using Python code and then converted into area-weighted average rainfall with MATLAB and ArcGIS software. The geomorphological map of the study area was obtained from the Bhukosh website, which is maintained by the Geological Survey of India (GSI). Hydrogeological data were manually collected from the Office of the Senior Geohydrologist in Raipur, Chhattisgarh. Furthermore, the LULC map was created using 30-meter resolution imagery from LANDSAT 8-9 OLI/TIRS C2 L1 for the year 2022, which was downloaded from the USGS Earth Explorer website managed by NASA. Observation well discharge data have been collected from the Central Ground Water Board (CGWB) website report. SRTM DEM data were also downloaded from the USGS earth explorer website. Detailed information on data acquisition is shown in Table 1. Processing of entire GIS data sets was done by ArcGIS 10.5 and ERDAS 2016.

Table 1

Different input parameters and their sources

S. No.DataParametersSource
1. Survey of India Top sheet no. 64 B, C, F and G Location map validation Survey of India (SOI)
Source:https://onlinemaps.surveyofindia.gov.in/FreeMapSpecification.aspx 
2. Soil data Soil texture, Soil type FAO Soil Map
Source:https://www.fao.org/soils-portal/data-hub/soil-maps-and-databases/faounesco-soil-map-of-the-world/en/ 
3. Lineament data Lineament buffer Bhuvan (ISRO)
Source:https://bhuvan-app1.nrsc.gov.in/thematic/thematic/index.php 
4. Meteorological data Rainfall India Meteorological Department (IMD), Pune (India)
IMD gridded data (0.25° × 0.25°)
Source:https://www.imdpune.gov.in/cmpg/Griddata/Rainfall_25_NetCDF.html 
5. Geomorphology data Geology, Geomorphology data Geological Survey of India (Bhukosh)
Source: https://bhukosh.gsi.gov.in/Bhukosh/Public 
6. Hydro-geological data Groundwater level data of study area Office of Senior Geohydrologist, Raipur
State Groundwater Department Chhattisgarh 
7. USGS Satellite Imageries SRTM Digital Elevation Model (DEM) of 30 m resolution USGS Earth Explorer
Source: https://earthexplorer.usgs.gov/ 
8. USGS Landsat 8 Imageries 2022 Landsat 8-9 OLI/TIRS C2 L1 of 30-m resolution USGS Earth Explorer
Source: https://earthexplorer.usgs.gov/ 
9. Well discharge data Observation well discharge data CGWB NAQUIM report of Kawardha district
Source: http://cgwb.gov.in/aquifer-mapping 
S. No.DataParametersSource
1. Survey of India Top sheet no. 64 B, C, F and G Location map validation Survey of India (SOI)
Source:https://onlinemaps.surveyofindia.gov.in/FreeMapSpecification.aspx 
2. Soil data Soil texture, Soil type FAO Soil Map
Source:https://www.fao.org/soils-portal/data-hub/soil-maps-and-databases/faounesco-soil-map-of-the-world/en/ 
3. Lineament data Lineament buffer Bhuvan (ISRO)
Source:https://bhuvan-app1.nrsc.gov.in/thematic/thematic/index.php 
4. Meteorological data Rainfall India Meteorological Department (IMD), Pune (India)
IMD gridded data (0.25° × 0.25°)
Source:https://www.imdpune.gov.in/cmpg/Griddata/Rainfall_25_NetCDF.html 
5. Geomorphology data Geology, Geomorphology data Geological Survey of India (Bhukosh)
Source: https://bhukosh.gsi.gov.in/Bhukosh/Public 
6. Hydro-geological data Groundwater level data of study area Office of Senior Geohydrologist, Raipur
State Groundwater Department Chhattisgarh 
7. USGS Satellite Imageries SRTM Digital Elevation Model (DEM) of 30 m resolution USGS Earth Explorer
Source: https://earthexplorer.usgs.gov/ 
8. USGS Landsat 8 Imageries 2022 Landsat 8-9 OLI/TIRS C2 L1 of 30-m resolution USGS Earth Explorer
Source: https://earthexplorer.usgs.gov/ 
9. Well discharge data Observation well discharge data CGWB NAQUIM report of Kawardha district
Source: http://cgwb.gov.in/aquifer-mapping 

Generation of thematic layers

To develop a sustainable groundwater recharge plan by the identification of groundwater potential zones, different thematic layers were prepared using remote sensing data, topographic maps, geological maps reports and data from different sources. The satellite data, Landsat 8-9 OLI/TIRS C2 L1 of 30 m resolution geocoded imagery pertaining to year 2022 was downloaded from USGS Earth Explorer for preparation of land use/land cover map. It was generated from the Landsat 8 data in ERDAS IMAGINE software using supervised classification techniques. The different layers such as geology, geomorphology, lineament, drainage density soil, slope, rainfall, and groundwater fluctuation maps were digitized using ERDAS IMAGINE and ArcGIS software SRTM DEM has been used to prepare the slope map of the study area. The pre-monsoon and post-monsoon groundwater level maps were prepared from the data taken from the State Groundwater Department of the bore wells and the inventory wells using interpolation techniques with the help of ArcGIS software. Rainfall maps were also prepared with the help of ArcGIS software using interpolation techniques.

Spatial data processing

The regular grid-based computation method, with the aid of GIS software, is considered a powerful way to store, process and analyze urban structure, and spatial growth patterns (Yeh & Li 2002). The basic concept of such a method is to divide the whole space into continuous square grids after settling on a predefined cell size (henceforth grid cell); however, an appropriate selection has to be justified enough, which is always a critical step for initializing the geoprocessing environmental settings. A common argument suggests that cell size should be sufficiently large to ensure accurate computational results and interpretation. To determine the optimal cell size, the coverage area of the largest physical features was considered. In Dhaka city, only a small number of buildings (8 out of 270,392) had a coverage area exceeding 10,000 m2. Therefore, a grid cell size of 100 × 100 m was chosen, which aligns with common practices in the literature (Yeh & Li 2002; Larondelle et al. 2014) and is based on Dhaka city's building feature. The spatial variability of ECD has to be localized as per one of the basic requirements in spatial analysis, so it was assumed that man-made activities can be considered almost constant within a grid cell size of (100 × 100) sq.m ground area.

Analytical hierarchy process

The analytical hierarchy process (AHP) was developed by Saaty in 1980 for pairwise comparison for understanding the priority from all the alternatives. This method indicates the relative importance between each other towards output result. This method solves the decision-making process, because of its simplicity and accuracy. This method is adopted for tangible and non-tangible criteria. The AHP process is based on three basic principles: decomposition, prioritization, and synthesis (Saaty 1980).

The decomposition is a complex problem to be decomposed I to different levels which can be manageable with few elements. It is a top to bottom process. The level of decomposing is simplifying the complex problem. A lower level of hierarchy defines detailed elements of the previous one level. The higher one is denoted as global objective.

In the prioritization, the AHP conducts comparisons of the decision elements in a systematic way to determine their priority scales. While paired comparisons are a straightforward concept (crucial for involving non-technically trained experts), their application in the AHP is innovative and unique. Only homogeneous elements are to be compared. The pairwise comparison judgments are made using the AHP's fundamental scale, which consists of absolute values from 1 to 9 and their reciprocals, as shown in Table 2.

Table 2

Saaty's scale of importance

ValueDescription
Extreme significance 
Very strong importance 
Extremely to extreme importance 
Strong plus importance 
Strong importance 
Moderate importance 
Moderate importance plus 
Weak relevance 
Equal importance 
ValueDescription
Extreme significance 
Very strong importance 
Extremely to extreme importance 
Strong plus importance 
Strong importance 
Moderate importance 
Moderate importance plus 
Weak relevance 
Equal importance 

MCDA is a highly effective tool that is widely used to solve difficult decision-making issues. The AHP involves structuring a complex decision into a hierarchical model, starting with defining the main goal and breaking it down into criteria and sub-criteria. Decision-makers then use pairwise comparisons to evaluate the relative importance of each criterion and alternative, constructing a matrix to derive weights. These weights are used to rank alternatives, with a consistency check ensuring the judgments are logical. The final step aggregates these results to identify the best alternative, and optionally, sensitivity analysis can be performed to test the robustness of the decision against variations in the inputs (Machiwal et al. 2011; Neshat et al. 2014; Roozbahani et al. 2018; Dey et al. 2021). The step by step concise outline involved in AHP is as follows:

Construct pairwise comparison matrices: For each level of the hierarchy (e.g., criteria, sub-criteria), construct a matrix where each element represents the relative importance of criterion i compared to the criterion j.

Normalize the matrix: Normalize the matrix by dividing each element by the sum of its column. For element in the matrix:
(1)
Compute the priority vector: Calculate the average of the normalized values in each row to get the priority vector. For criterion :
(2)
where is the weight of criterion i, and n is the number of criteria?

Calculate the consistency ratio (CR):

  • Consistency Index (CI) is computed and is provided in Equation (3) to evaluate the consistency of relative evaluations of all the criteria and sub-criteria.
    (3)
  • Compute the CR by dividing CI by the Random Consistency Index (RI):
    (4)
    where RI is the average consistency index for random matrices of size n (provided in AHP Table 3). The analysis can proceed with a CR of 0.10 or less, according to Saaty. The judgment must be evaluated if the consistency value is more than 0.10 in order to determine the underlying causes of the inconsistency and make the required modifications (Saaty 1980). If the CR value is 0, which suggests the pairwise comparison has perfect consistency. Because the threshold value is not greater than 0.1, the judgements matrix is comparatively consistent (Bera et al. 2019).
Table 3

Random consistency index values as per n values

N12345678910
RCI 0.58 0.90 1.12 1.24 1.32 1.41 1.45 1.49 
N12345678910
RCI 0.58 0.90 1.12 1.24 1.32 1.41 1.45 1.49 

Source: Saaty 1980.

Aggregate of the results:

Multiply the priority weights of criteria by the weights of alternatives for each criterion. For alternative with criteria weights and alternative weights :
(5)
where is the score for alternative .

These equations facilitate a systematic approach to decision-making in AHP, ensuring that decisions are based on structured comparisons and consistent judgments (Malczewski 1999; Saaty 2004; Achu et al. 2020).

Estimation of groundwater potential zone map

For the purpose, weighted overlay analysis has been performed in ArcGIS 10.5. Initially, all thematic layers in the add layer option in ArcGIS were reclassified as per required classes. Further, the weight overlay tool has been utilized using AHP-based assigned weight values. Finally, run the tool and get the result output as groundwater potential zone map. Once the analysis is complete, classify the output into zones (e.g., very high, high, moderate, low, very low potential) and create a final groundwater potential zone map. For validation of the result, observation well data (ground level and discharge data of 2021–22) has been used for the study. Detailed methodology is shown in Figure 2.
Figure 2

Flowchart for the development of sustainable groundwater recharge plan.

Figure 2

Flowchart for the development of sustainable groundwater recharge plan.

Close modal

Identification of suitable area/location for groundwater recharge

For identification of appropriate area a groundwater recharge zone map is prepared using the MCDA technique. In the groundwater potential zones maps areas showing low discharge will be given first preference for groundwater recharge planning (Koohbanani et al. 2018). In general post monsoon the groundwater level map is observed for identifying the appropriate area for groundwater recharge; the areas showing high depth of water in the post monsoon season will be selected for recharge planning.

Subsurface storage estimation

Assessment of subsurface storage is done on the basis of average depth to water level records. For determining the amount of storage space that is available, contour maps that have been created using average post-monsoon water level data and acceptable contour intervals can be employed. The overall area in which the inter-contour areas between subsequent contours are determined and the amount of subsurface storage space that is open for recharging is determined by multiplying the volume of water below a particular cut-off level by the specific yield of the aquifer material. The cut-off water level is chosen in this manner to prevent water logging in the area after recharging. For our study, 2 m is chosen as the cutoff level.

Estimation of water required for recharge

It is necessary to evaluate the real source water requirement in order to determine the subsurface storage area. It is necessary to specify the average recharge efficiency of each individual structure based on the knowledge gathered from field tests. The amount of water needed for artificial recharge is determined by multiplying the volume of subsurface storage capacity by the reciprocal of the planned structure's proposed recharge efficiency. This results in the total volume of actual source water needed at the surface. In the present study, site selection criteria adopted for check dam, percolation tank and farm pond are illustrated in Table 4.

Table 4

Site selection criteria for the structures

Type of structureMaximum water level (m)SlopeStorage lossLanduseStream order
Check Dam 4–5 <15% Low Agriculture, barren land >Third order 
Percolation Tank 6–7 <10% Moderate–high Agriculture, barren land Second to third order 
Farm Pond 2–2.5 <5% Moerate–low Agriculture, barren land First order 
Type of structureMaximum water level (m)SlopeStorage lossLanduseStream order
Check Dam 4–5 <15% Low Agriculture, barren land >Third order 
Percolation Tank 6–7 <10% Moderate–high Agriculture, barren land Second to third order 
Farm Pond 2–2.5 <5% Moerate–low Agriculture, barren land First order 

Geology map

The geologic history has a significant impact on groundwater distribution and occurrence in every terrain (Yeh & Shellnutt 2016). The geological formations in the district of Kawardha are diverse and span a wide range of ages. In the southwestern part of the district, the Nandgaon Group, which dates back to the Palaeo to Meso Proterozoic period (2,500–1,600 million years), is composed of meta-basalt and Malanjkhand Granitoids. The meta-basalt is exposed in this region, while the Granitoids range in composition from granodiorite to quartz diorite. Additionally, during the Meso-Proterozoic era (2,000–1,600 million years ago), rocks such as quartz diorite, granite porphyry, and basic dykes were inserted into the Bilaspur-Raigarh-Surguja Group. This group comprises the Chandrapur and Raipur groups. The Chandrapur group is primarily composed of arenite, ferruginous sandstone, and polymictic conglomerate, with the ferruginous sandstone also containing mudrock sections. In contrast, the Raipur zones include the Gunderdehi, Chandi, Tarenga, Hirri, and Maniyari formations, each having unique characteristics and compositions (CGWB NAQUIM Report 2022). The Deccan Trap is located in the western and northern portions of the district, which formed during the Cretaceous to Palaeogene period around 65–60 million years ago. The Deccan Trap is primarily composed of unclassified basalt, with the Linga Formation being referred to in the central-western part. Nonporphyritic basaltic lava flows make up the Linga Formation. Additionally, the elevated plateaus of the Deccan Trap often have a layer of laterite, formed during the Cainozoic era. In the present study, 12 classes are involved in the geological map of the study area. A detailed description of each class is shown in Table 5. A geology map of the study area is shown in Figure 3(a).
Table 5

Geological features of the study area

Geology classesArea (km2)Area (%)
Bilaspur-Raigarh-Surguja belt gp. 25.16 0.60045 
Chandarpurgp. (Chhattisgarh sgp.) 5.08 0.12124 
Chhotanagpur gneissic complex 67.92 1.62094 
Chilpigp. 766.89 18.3022 
Deccan trap 418.82 9.99532 
Laterite/bauxite 73.94 1.76461 
Malanjkhand granite 172.96 4.12777 
Nandgaongp. (pitepani intrusive) 230.06 5.49048 
Nandgaongp. (bijli rhyolite fm.) 410.41 9.79461 
Raipur gp. (Chhattisgarh sgp.) 1,859.78 44.3845 
Singhoragp. (Chhattisgarh sgp.) 156.5 3.73494 
Tirodi gneissic complex gp. 2.64 0.063 
Geology classesArea (km2)Area (%)
Bilaspur-Raigarh-Surguja belt gp. 25.16 0.60045 
Chandarpurgp. (Chhattisgarh sgp.) 5.08 0.12124 
Chhotanagpur gneissic complex 67.92 1.62094 
Chilpigp. 766.89 18.3022 
Deccan trap 418.82 9.99532 
Laterite/bauxite 73.94 1.76461 
Malanjkhand granite 172.96 4.12777 
Nandgaongp. (pitepani intrusive) 230.06 5.49048 
Nandgaongp. (bijli rhyolite fm.) 410.41 9.79461 
Raipur gp. (Chhattisgarh sgp.) 1,859.78 44.3845 
Singhoragp. (Chhattisgarh sgp.) 156.5 3.73494 
Tirodi gneissic complex gp. 2.64 0.063 
Figure 3

(a) Geology map; (b) geomorphology map; (c) LULC map; (d) stream order map; (e) drainage density map; (f) rainfall map; (g) slope map; (h) lineament density map; (i) pre-monsoon GWL map; (j) post-monsoon GWL map; (k) soil map; and (l) DEM map.

Figure 3

(a) Geology map; (b) geomorphology map; (c) LULC map; (d) stream order map; (e) drainage density map; (f) rainfall map; (g) slope map; (h) lineament density map; (i) pre-monsoon GWL map; (j) post-monsoon GWL map; (k) soil map; and (l) DEM map.

Close modal

Geomorphology map

Geomorphology, which depicts a region's topography and landforms, is one of the crucial factors usually utilized to define groundwater potential zones. Additionally, it provides information on processes including temperature variations, geochemical reactions, and the distribution of different landform features (Kumar et al. 2016a, 2016b; Thapa et al. 2017). In terms of landforms, the area can be divided into two separate regions: (1) the hilly western part and (2) the flat eastern region. The Maikal and Mangata ranges are located in the western region of the district. The district slopes generally in the direction of the east. Landforms such as hills, valleys, plains, plateaus brought about by erosion, and floodplains created by filled river beds may be found in the district's center and southern regions. There are plateaus, hills, and high-level plateaus of volcanic origin in the district's northern portion. The western region of the district is a part of the Banjara sub-basin, while the majority of it is a part of the Seonath sub-basin. The surface drainage system is formed by the tributaries of the Seonath and Banjara rivers. Except for the southwest portion of the Banjara sub-basin, where the gradient is towards the north, the area's overall gradient is toward the southeast. The Maikal mountain ranges of the Satpura border the district's northern and western regions. The southeast has a minimum elevation of 300 m above sea level, while the tallest mountain, Kesmarda, is situated in the northwest at an elevation of 931 m above sea level. In the present study, 12 classes are involved in the geomorphological map of the study area. A detailed description of each class is shown in Table 6. A geomorphology map of the study area is shown in Figure 3(b).

Table 6

Geomorphological features of the study area

Geomorphologic unitsArea (sq.km)Area (%)
Active flood plain 4.65 0.001 
Dam and reservoir 7.49 0.002 
Highly dissected structural hills and valleys 552.92 0.137 
Low dissected denudational hills and valleys 69.42 0.017 
Low dissected denudational upper plateau 27.01 0.007 
Low dissected structural hills and valleys 13.19 0.003 
Moderately dissected denudational hills and valleys 3.75 0.001 
Moderately dissected denudational upper plateau 47.35 0.012 
Moderately dissected structural hills and valleys 78.34 0.019 
Moderately dissected structural upper plateau 691.06 0.171 
Pediment pediplain complex 2,503.15 0.620 
Pond 12.84 0.003 
River 20.26 0.005 
Waterbodies 6.47 0.002 
Geomorphologic unitsArea (sq.km)Area (%)
Active flood plain 4.65 0.001 
Dam and reservoir 7.49 0.002 
Highly dissected structural hills and valleys 552.92 0.137 
Low dissected denudational hills and valleys 69.42 0.017 
Low dissected denudational upper plateau 27.01 0.007 
Low dissected structural hills and valleys 13.19 0.003 
Moderately dissected denudational hills and valleys 3.75 0.001 
Moderately dissected denudational upper plateau 47.35 0.012 
Moderately dissected structural hills and valleys 78.34 0.019 
Moderately dissected structural upper plateau 691.06 0.171 
Pediment pediplain complex 2,503.15 0.620 
Pond 12.84 0.003 
River 20.26 0.005 
Waterbodies 6.47 0.002 

LULC map

LULC provides groundwater requirements as well as critical data on infiltration, soil moisture, surface water, groundwater, etc. In this study, the LULC map is prepared using supervised classification from Landsat 8 imageries of the year 2022. Six classes of LULC, viz. water bodies, forest, barren land, settlement, agricultural and fallow land were identified and demarcated as shown in Table 7 and Figure 3(c).

Table 7

LULC map of the study area

No.ClassArea (sq.km)Area (%)
Water body 76.47 1.8279 
Forest 1,138.17 27.2062 
Barren land 112.69 2.69368 
Settlement 206.6 4.93846 
Agriculture 1,975.13 47.2125 
Fallow land 674.43 16.1212 
No.ClassArea (sq.km)Area (%)
Water body 76.47 1.8279 
Forest 1,138.17 27.2062 
Barren land 112.69 2.69368 
Settlement 206.6 4.93846 
Agriculture 1,975.13 47.2125 
Fallow land 674.43 16.1212 

Drainage density map

Drainage density is primarily derived from the demarcation of streams and stream order shown in Figure 3(d). Drainage density has a significant impact on both groundwater availability and contamination (Ganapuram et al. 2009). The drainage system depends on the lithology and is a significant measure of infiltration rate. Permeability and drainage density have the opposite relationships. Therefore, it is essential in determining the groundwater potential zone. The drainage density of a drainage basin is determined by dividing the total length of its rivers by the drainage basin's area. The area's potential for groundwater is less favourable due to low infiltration caused by high drainage density. The potential of the groundwater is increased more by excessive infiltration caused by low drainage density. The drainage density was prepared and reclassified in five classes, viz. 0–0.706 km/sq.km, 0.707–1.2 km/sq.km, 1.21–1.71 km/sq.km, 1.72–2.28 km/sq.km and 2.29–3.4 km/sq.km shown in Figure 3(e).

Rainfall map

Rainfall is the main source of water in the hydrologic cycle and the main determinant of a region's groundwater (Nouaceur & Murărescu 2016). For the current study IMD gridded rainfall data were used for the year 2022 for rainfall map preparation. The spatial distribution of rainfall was mapped using the IDW interpolation method (Earls & Dixon 2007). The rainfall map is categorizing into five classes (1,200–1,400 mm), (1,400–1,500 mm), (1,500–1,600 mm), (1,600–1,700 mm) and (1,700–1,800 mm) shown in Figure 3(f).

Slope map

A crucial topographical component is the slope, which describes how steep the ground surface is. Slope gives important information on the types of geologic and geodynamic processes occurring at the level of the region (Riley et al. 1999). Overall, the topography in the district varies between 273 and 1,016 m (743 m elevation difference) above mean sea level. Slope of the surface primarily influences surface run-off and infiltration rate (Singh et al. 2013). Larger slopes create less recharge because the water from precipitation does not have enough time to dwell in the saturated zone and rushes fast down a steep hill when it rains (ref). The slope map was prepared using SRTM 30 m DEM and reclassified into five classes (0–2%), (2–4%), (4–6%), (6–12%) and (12–25%) shown in Figure 3(g).

Lineament density map

Lineaments are structurally regulated linear or curved characteristics. It can be recognized by their comparatively linear alignments in the satellite image (Nampak et al. 2014; Nag & Kundu 2016). Lines represent the faulting and fracture zones leading to increased secondary porosity and permeability. By adding a Web Mapping Service (WMS) layer to ArcMap 10.5, extracting the research region's lineaments, and using the line density tool to create a lineament density map, lineaments of the study area were taken from Bhuvan NRSC data. The lineament density is reclassified in five classes 0–0.065 km/Sq.km, 0.065–0.195 km/Sq.km, 0.195–0.344 km/Sq.km, 0.344–0.532 km/Sq.km and 0.532–0.924 km/Sq.km. The ranks for lineament density are determined by how closely spaced the lineaments are to one another. The strength of the groundwater potential is shown to diminish as the distance from the lineaments rises, shown in Figure 3(h).

Soil map

The type of soil has a big impact on how much water may seep into underground formations, which also has an impact on groundwater recharge (Das 2017). The hydraulic characteristics of the soil and its texture are the main factors considered for determining the rate of infiltration. In our study the soil map was prepared by using FAO soil data and a total of four types of soil are found under the study area: clay, sandy loam, loam and cay loam, as shown in Figure 3(k). The elevation difference of the study area is about 743 m, which is estimated by using the DEM shown in Figure 3(l).

Groundwater fluctuations

Groundwater fluctuation is dependent on pre-monsoon and post-monsoon data sets of the study area. The pre- and post-monsoon groundwater levels vary from 1.6–8.3 m to 0.6–7.2 m, respectively, shown in Figure 3(i) and 3(j). Whereas the groundwater fluctuation of the study area falls under 0.3–5.4 m, shown in Figure 4(b). The results of the research indicate that the groundwater level in the study area is largely steady throughout the post-monsoon season while showing a minor upward tendency in the Bodla block during the pre-monsoon season. These results are in line with previous studies, which showed that groundwater levels in the area were steady throughout the post-monsoon season and slightly rising during the pre-monsoon season. The study indicates the need for monitoring groundwater levels in order to ensure sustainable use of this essential resource and offers useful information for the management of water resources in the region, shown in Figure 4(b).
Figure 4

(a) Groundwater potential zone map; (b) groundwater level fluctuation map; and (c) artificial groundwater recharge structures points map for the study area.

Figure 4

(a) Groundwater potential zone map; (b) groundwater level fluctuation map; and (c) artificial groundwater recharge structures points map for the study area.

Close modal

Groundwater potential zone

Planning and development of the region's water resources must consider and manage the region's hard rock geology (CGWB NAQUIM Report). Similar to other hard rock areas, Kawardha's groundwater potential is greatly influenced by geology. In our study, we assigned the highest preference to geology for delineating groundwater potential, recognizing its importance. Conversely, in our carefully constructed pairwise comparison matrix, rainfall obtained the lowest rank and weightage depending on the suitability of the research location. By considering both Chhattisgarh and Kawardha's hard rock geology, we can enhance water resource management strategies in the region. The comparison matrix was created based on how significant each component was and the way it influenced the groundwater potential shown in Table 8. After preparation of the pairwise comparison matrix normalization of the weights was done. Following this, the important calculation which is required for the acceptability of the matrix were calculated (Table 9) as per the standard formula used by many researchers (Saranya & Saravanan 2020; Dar et al. 2021; Pande et al. 2021).

Random consistency index = 1.41 (value taken from Table 2)

Table 8

Pairwise comparison matrix

ParametersGeologyGeomorphologyLULCSoilLineamentDrainage densitySlopeRainfall
Geology 
Geomorphology 
LULC 0.33 0.5 
Soil 0.33 0.33 0.33 
Lineament 0.33 0.33 0.33 0.33 
Drainage Density 0.33 0.33 0.2 0.33 0.5 
Slope 0.11 0.11 0.2 0.25 0.5 0.5 
Rainfall 0.14 0.25 0.14 0.25 0.25 0.33 0.33 
ParametersGeologyGeomorphologyLULCSoilLineamentDrainage densitySlopeRainfall
Geology 
Geomorphology 
LULC 0.33 0.5 
Soil 0.33 0.33 0.33 
Lineament 0.33 0.33 0.33 0.33 
Drainage Density 0.33 0.33 0.2 0.33 0.5 
Slope 0.11 0.11 0.2 0.25 0.5 0.5 
Rainfall 0.14 0.25 0.14 0.25 0.25 0.33 0.33 
Table 9

Normalization of the comparison matrix

ParameterGeologyGeomorphologyLULCSoilLineamentDrainage DensitySlopeRainfallNPE
Geology 0.28 0.26 0.42 0.27 0.21 0.17 0.28 0.21 0.26 
Geomorphology 0.28 0.26 0.28 0.27 0.21 0.17 0.28 0.12 0.23 
LULC 0.09 0.13 0.14 0.27 0.21 0.28 0.15 0.21 0.19 
Soil 0.09 0.09 0.05 0.09 0.21 0.17 0.12 0.12 0.12 
Lineament 0.09 0.09 0.05 0.03 0.07 0.11 0.06 0.12 0.08 
Drainage Density 0.09 0.09 0.03 0.03 0.04 0.06 0.06 0.09 0.06 
Slope 0.03 0.03 0.03 0.02 0.04 0.03 0.03 0.09 0.04 
Rainfall 0.04 0.06 0.02 0.02 0.02 0.02 0.01 0.03 0.03 
ParameterGeologyGeomorphologyLULCSoilLineamentDrainage DensitySlopeRainfallNPE
Geology 0.28 0.26 0.42 0.27 0.21 0.17 0.28 0.21 0.26 
Geomorphology 0.28 0.26 0.28 0.27 0.21 0.17 0.28 0.12 0.23 
LULC 0.09 0.13 0.14 0.27 0.21 0.28 0.15 0.21 0.19 
Soil 0.09 0.09 0.05 0.09 0.21 0.17 0.12 0.12 0.12 
Lineament 0.09 0.09 0.05 0.03 0.07 0.11 0.06 0.12 0.08 
Drainage Density 0.09 0.09 0.03 0.03 0.04 0.06 0.06 0.09 0.06 
Slope 0.03 0.03 0.03 0.02 0.04 0.03 0.03 0.09 0.04 
Rainfall 0.04 0.06 0.02 0.02 0.02 0.02 0.01 0.03 0.03 

According to Saaty (1984), it is adequate to continue the analysis with a CR of 0.10 or below. The judgment must be evaluated if the consistency value is more than 0.10 in order to determine the underlying causes of the inconsistency and make the required modifications. If the CR value is 0, this suggests the pairwise comparison has perfect consistency. Because the threshold value is no higher than 0.1, the judgements matrix is essentially consistent. Furthermore, influencing percentage is decided by the normalized principal eigenvector, which is shown in Table 10.

Table 10

Percentage distributions of parameters using AHP-based analysis in GIS overlay analysis

S. No.ParameterInfluence (%)Feature classesWeights (1–9)
1. Geology 27 Bilaspur-Raigarh-Surguja belt gp. 
Chandarpurgp. (Chhattisgarh sgp.) 
Chhotanagpur gneissic complex 
Chilpigp. 
Deccan trap 
Laterite/bauxite 
Malanjkhand granite 
Nandgaongp. (pitepani intrusive) 
Nandgaongp.(bijli rhyolite fm.) 
Raipur Gp.(Chhattisgarh sgp.) 
SinghoraGp. (Chhattisgarh sgp.) 
Tirodi gneissic complex gp. 
2. Geomorphology 24 Active Flood Plain 
Dam and Reservoir 
Highly dissected structural hills and valleys 
Low dissected denudational hills and valleys 
Low dissected denudational upper plateau 
Low dissected structural hills and valleys 
Moderately dissected denudational hills and valleys 
Moderately dissected denudational upper plateau 
Moderately dissected structural hills and valleys 
Moderately dissected structural upper plateau 
Pediment pediplain complex 
Pond 
River 
3. LULC 18 Agriculture 
Forest-cover 
Fallow-land 
Settlement 
Barren-land 
Water body 
4. Soil 12 Clay loam 
Loam 
Sandy-loam 
Clay 
5. Lineament (km/sq.km) 0–0.065 
0.06–0.195 
0.195–0.344 
0.344–0.5327 
0.5327–0.924 
6. Drainage density (km/sq.km) 0–0.68 
0.68–1.4 
1.5–2 
2.1–2.7 
2.8–3.4 
7. Rainfall (mm) 1,200–1,400 
1,400–1,500 
1,500–1,600 
1,600–1,700 
1,700–1,800 
8. Slope (%) 0–2 
2–4 
4–8 
8–12 
12–25 
S. No.ParameterInfluence (%)Feature classesWeights (1–9)
1. Geology 27 Bilaspur-Raigarh-Surguja belt gp. 
Chandarpurgp. (Chhattisgarh sgp.) 
Chhotanagpur gneissic complex 
Chilpigp. 
Deccan trap 
Laterite/bauxite 
Malanjkhand granite 
Nandgaongp. (pitepani intrusive) 
Nandgaongp.(bijli rhyolite fm.) 
Raipur Gp.(Chhattisgarh sgp.) 
SinghoraGp. (Chhattisgarh sgp.) 
Tirodi gneissic complex gp. 
2. Geomorphology 24 Active Flood Plain 
Dam and Reservoir 
Highly dissected structural hills and valleys 
Low dissected denudational hills and valleys 
Low dissected denudational upper plateau 
Low dissected structural hills and valleys 
Moderately dissected denudational hills and valleys 
Moderately dissected denudational upper plateau 
Moderately dissected structural hills and valleys 
Moderately dissected structural upper plateau 
Pediment pediplain complex 
Pond 
River 
3. LULC 18 Agriculture 
Forest-cover 
Fallow-land 
Settlement 
Barren-land 
Water body 
4. Soil 12 Clay loam 
Loam 
Sandy-loam 
Clay 
5. Lineament (km/sq.km) 0–0.065 
0.06–0.195 
0.195–0.344 
0.344–0.5327 
0.5327–0.924 
6. Drainage density (km/sq.km) 0–0.68 
0.68–1.4 
1.5–2 
2.1–2.7 
2.8–3.4 
7. Rainfall (mm) 1,200–1,400 
1,400–1,500 
1,500–1,600 
1,600–1,700 
1,700–1,800 
8. Slope (%) 0–2 
2–4 
4–8 
8–12 
12–25 

Groundwater is a resource that can be restored, but in the previous four to five decades, diverse anthropogenic activities and unbalanced development have drastically reduced the recharge of this important life-sustaining resource. A thorough understanding of groundwater potential is essential for effective planning and sustainable development of an area. The planning and building of the remedial structures that will improve the processes of groundwater recharge require such information. The eight parameters listed above were taken into consideration while determining groundwater potential zones. Based upon the pairwise comparison matrix, weights were assigned in accordance with these characteristics when using the AHP technique. The values, 0.11 and 0.074, respectively, were found to be the CR and CI for high GPZs. The Kawardha district's groundwater potential zones map, which was created using ArcGIS, is a useful resource for learning how groundwater potential is distributed over various places. According to the analysis, 60.46% of the district's total area, or the bulk of the district, has moderate groundwater potential. This suggests a moderate capacity for groundwater availability in these areas. Following the moderate zone, the high potential zone covers 20.91% of the district, encompassing parts of Bodla, Pandariya, and Kawardha blocks, where groundwater availability is generally favorable. On the other side, 17.60% of the district is made up of low potential areas, including the Bodla region and a few blocks in Pandariya that have lower groundwater potential. Smaller portions of the district are represented by the extremely low and high potential zones, which make up 0.50 and 0.53% of the district's total area, respectively. While the very low potential areas indicate regions with highly constrained groundwater supply, the very high potential areas emphasize the zones with the district's greatest groundwater potential shown in Figure 4(a). The distribution of groundwater potential zones according to area and percentage is shown in Table 11.

Table 11

Groundwater potential zones and area distribution

ClassArea (km2)Area (%)
Very low 20.97 0.50 
Low 727.62 17.60 
Moderate 2,498.72 60.46 
Good 864.3 20.91 
Very good 20.79 0.53 
ClassArea (km2)Area (%)
Very low 20.97 0.50 
Low 727.62 17.60 
Moderate 2,498.72 60.46 
Good 864.3 20.91 
Very good 20.79 0.53 

Validation of groundwater potential zones map

The examination of eight contributing elements was added into the groundwater potential zones map produced by ArcGIS. Discharge data from the research region were used to verify the accuracy of the map that was generated. In the present study, 19 observation wells data (2021–22) were utilized for validation of groundwater potential zones. These observation well locations were overlaid onto the groundwater potential map, allowing for a comparative analysis between the predicted potential zones and the observed discharge values. By placing the positions of the observation wells across the groundwater potential map, a comprehensive assessment of the map's accuracy and reliability was conducted (Maghrebi et al. 2023). The efficacy of the model and the validity of the identified potential zones can both be determined by comparing the projected potential zones with the actual groundwater discharge data. This validation procedure (Table 12) aids in ensuring that the generated groundwater potential zones map corresponds to the actual groundwater dynamics that have been observed (Yeganeh et al. 2024). Following is the accuracy assessment table.

Table 12

Geographical entities of the station along with discharge

S No.StationsLatitudeLongitudeDischarge (lps)Condition
Birendranagar 21.82 81.18 Agree 
Charbhata 21.88 81.23 10.85 Agree 
Dasrangpur 21.85 81.38 7.4 Agree 
Dasranpur OW 21.85 81.38 9.31 Agree 
Indauri 21.95 81.37 2.6 Disagree 
Kawardha 22.03 81.22 6.4 Disagree 
Kawardha OW 22.03 81.22 23 Disagree 
SahaspurLohara 21.83 81.13 Agree 
Ranvirpur 21.83 81.18 2.6 Agree 
10 Pondi 22.12 81.27 5.76 Agree 
11 Kishangarh 22.22 81.50 3.5 Agree 
12 Pandaria 22.22 81.40 8.7 Agree 
13 Kunda 22.12 81.45 3.5 Agree 
14 Damapur 22.05 81.40 5.8 Agree 
15 Bagarra 22.15 81.48 7.14 Agree 
16 Pandatarai 22.19 81.33 14.5 Agree 
17 Bodla 22.15 81.22 14.5 Agree 
18 Raveli 22.10 81.33 0.2 Disagree 
19 Bhaluchwa 22.09 81.22 7.14 Agree 
S No.StationsLatitudeLongitudeDischarge (lps)Condition
Birendranagar 21.82 81.18 Agree 
Charbhata 21.88 81.23 10.85 Agree 
Dasrangpur 21.85 81.38 7.4 Agree 
Dasranpur OW 21.85 81.38 9.31 Agree 
Indauri 21.95 81.37 2.6 Disagree 
Kawardha 22.03 81.22 6.4 Disagree 
Kawardha OW 22.03 81.22 23 Disagree 
SahaspurLohara 21.83 81.13 Agree 
Ranvirpur 21.83 81.18 2.6 Agree 
10 Pondi 22.12 81.27 5.76 Agree 
11 Kishangarh 22.22 81.50 3.5 Agree 
12 Pandaria 22.22 81.40 8.7 Agree 
13 Kunda 22.12 81.45 3.5 Agree 
14 Damapur 22.05 81.40 5.8 Agree 
15 Bagarra 22.15 81.48 7.14 Agree 
16 Pandatarai 22.19 81.33 14.5 Agree 
17 Bodla 22.15 81.22 14.5 Agree 
18 Raveli 22.10 81.33 0.2 Disagree 
19 Bhaluchwa 22.09 81.22 7.14 Agree 

On the basis of Table 12, the accuracy assessment was performed by the following calculations (6):

Total agree points = 15

Total disagree points = 4
(6)

Groundwater recharge plan for the study area

Subsurface storage area is determined using the thickness of the accessible unsaturated zone (below 2 mbgl) in the post-monsoon season and the precise yield of the phreatic aquifer. The restriction against saturating the vadose zone below 2 m is maintained in order to avoid water logging and soil salinity. The total volume of unsaturated strata has been multiplied by the area's specific yield, or 0.03%, to get the real amount of water required to replace the aquifer up to 2 m (Kumar et al. 2016a, 2016b). A 150.21 M Cu m vadose zone that is suitable for artificial recharge exists in the research region. The actual volume of water required to saturate the vadose zone was calculated, and the net quantity of source water available was calculated by accounting for the artificial recharge structure's 75% efficiency. The value was multiplied by the reciprocal of 75% efficiency, which is 1.33. According to study, 199.5 Mm3 of water is needed for artificial recharge. 30% of the volume from each block is used to estimate the amount of water needed for recharge, and the corresponding runoff measurements for each block are 396.06, 381.01, 419.49, and 438.39 mm. Percolation tanks and check dams are acceptable for the research region for groundwater recharge, according to the study. These structures were proposed after several thematic layers were overlaid. The intersection point of lineament and drainage line was selected for groundwater recharging through check dam, similarly the sink point with lineament was selected for groundwater recharging through percolation tank. The results of the study made it abundantly evident that the district of Kawardha's groundwater recharge plans would be enough for the growth of sustainable agriculture (Bahraseman et al. 2024). Percolation tank (PT), check dam (CD), and farm pond (FP) are mainly suitable for artificial recharge facilities in the areas suggested (Figure 5) to be built and the number of the particular structure block wise (Table 13) are in the following Figure 4(c).
Table 13

Artificial groundwater recharge structures block wise distribution

S No.BlockCheck damPercolation tankFarm pond
Kawardha 11 
Pandariya 29 10 
Bodla 12 21 
Sahaspur Lohara 12 
S No.BlockCheck damPercolation tankFarm pond
Kawardha 11 
Pandariya 29 10 
Bodla 12 21 
Sahaspur Lohara 12 
Figure 5

Adopted soil and water conservation cum recharge structures along with description.

Figure 5

Adopted soil and water conservation cum recharge structures along with description.

Close modal

The present study provides a comprehensive approach to evaluating and implementing sustainable groundwater recharge strategies. It is an in-depth discussion considering specific conditions: High Slope: The mountainous region of Kawardha, with its varied topography, poses significant challenges for groundwater recharge. High slopes typically result in rapid runoff and less water infiltration, which reduces groundwater recharge potential. Rainfall characteristics: Heavy rainfall (1,200–1,800 mm): The region experiences substantial annual rainfall, which, while beneficial for recharge, also poses risks of erosion and runoff if not managed properly. Rainfall intensity: High-intensity rainfall increases surface runoff and erosion, reducing the effectiveness of traditional recharge methods and increasing the need for robust and adaptive solutions. Erratic distribution: Irregular rainfall distribution can create periods of drought followed by intense rainfall events, complicating recharge efforts and leading to variable groundwater levels (Das et al. 2019a, 2019b; Mirboluki et al. 2024). Lack of conservation measures: The absence of existing conservation measures exacerbates issues related to runoff and erosion, leading to reduced groundwater recharge and potentially impacting local water availability and agricultural productivity.

Solutions for sustainable groundwater recharge

Based on the characteristics and challenges outlined, the study suggests implementing various low-cost conservation measures. Following are the measures that can address the specific issues.

Farm ponds: Design criteria: Suitable for areas with slopes less than 5% and a maximum water level of 2–2.5 m. Advantages: Farm ponds can capture and store runoff during heavy rainfall, helping to reduce soil erosion and increase infiltration rates. They are especially useful for small-scale, localized water storage and can be integrated into agricultural practices to improve water availability during dry periods. Percolation ponds: Design criteria: Designed for slopes less than 10% and with a maximum water level of 6–7 m. Advantages: These ponds facilitate groundwater recharge by allowing water to percolate into the soil over a larger area. They are effective in regions with moderate slopes and can help enhance groundwater recharge by storing water for longer periods. Check dams: Design criteria: Effective on slopes less than 15% with a maximum water level of 4–5 m. Advantages: Check dams can slow down runoff, reduce soil erosion, and promote water infiltration. They are particularly useful in mountainous regions to manage high runoff and facilitate groundwater recharge in areas with steeper slopes.

Integration with AHP, RS, and GIS

AHP can be used to prioritize various recharge strategies based on factors such as cost, effectiveness, and ease of implementation. This multi-criteria decision analysis method helps in evaluating and selecting the most appropriate conservation measures for specific sites. RS technologies can provide detailed topographic, hydrological, and land use data. This data is crucial for identifying suitable locations for farm ponds, percolation ponds, and check dams. It also helps in monitoring changes in land cover and assessing the impact of conservation measures. GIS integrates and analyzes spatial data to identify optimal locations for implementing recharge measures. It can help visualize slope gradients, rainfall patterns, and land use, facilitating better planning and management of groundwater recharge projects.

Implementation strategy

Site selection

Use GIS to map out suitable locations for each type of conservation measure based on slope, rainfall intensity, and soil characteristics. Incorporate AHP to evaluate and prioritize sites based on various factors like cost, feasibility, and potential impact.

Community involvement

Engage local communities in the planning and implementation phases to ensure the measures align with local needs and practices. Community involvement can also enhance maintenance and sustainability. Monitoring and Evaluation: Implement a robust monitoring system using RS and GIS to track the effectiveness of conservation measures over time. Regularly assess water levels, infiltration rates, and soil erosion to adjust strategies as needed. Capacity building: Train local stakeholders in the construction and maintenance of farm ponds, percolation ponds, and check dams. Providing education on water conservation practices and the benefits of these measures can foster long-term sustainability.

In summary, addressing groundwater recharge in Kawardha requires a multifaceted approach that considers the unique topographical and climatic challenges. By integrating AHP, RS, and GIS with practical, low-cost conservation measures, it's possible to develop a sustainable groundwater management strategy that enhances water availability and mitigates the impacts of erratic rainfall and high runoff. As for the limitation of the study, initial cost is slightly higher then the agronomical measures.

The study employs a robust methodology integrating the AHP, RS, and GIS to assess and identify optimal groundwater recharge strategies. The use of these tools enables a detailed analysis of various parameters, including geology, geomorphology, LULC, soil, lineament, drainage density, slope, and rainfall. The key findings involved the groundwater potential zones, water deficit zones and recharge structures (farm pond, percolation pond and check dam) and their locations. The study's recharge plan estimates a vadose zone of 150 MCM suitable for artificial recharge, requiring 199.5 Mm³ of water. The plan considers various factors such as post-monsoon water levels, subsurface storage, surface water needs, and runoff estimation. The thematic maps and proposed structures are validated against observed data, demonstrating an overall map accuracy of approximately 79%. From the analysis and findings, targeted recharge measures, local community engagement, ongoing monitoring (groundwater levels) and holistic approach (integrate groundwater recharge plan) are recommended for sustainable groundwater recharge planning and management. This study provides a systematic approach for sustainable groundwater management in Kawardha District by identifying appropriate recharge zones and structures, thereby addressing the critical issues of water scarcity and poor management. The proposed plan, supported by detailed analyses and practical recommendations, lays the groundwork for improving groundwater recharge and ensuring more effective water resource management in the region.

In conclusion, the proposed groundwater recharge plan involves identifying suitable areas, estimating subsurface storage, determining surface water needs, calculating runoff, and recommending appropriate recharge structures. Effective implementation, ongoing monitoring, and community involvement are essential for sustainable groundwater management in Kawardha. The study provides valuable insights into regional water resource planning and offers a systematic approach to evaluating groundwater potential.

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

Amponsah
P. O.
,
Forson
E. D.
,
Sungzie
P. S.
&
Loh
Y. S. A.
(
2023
)
Groundwater prospectivitymodeling over the Akatsi Districts in the Volta Region of Ghana using the frequency ratio technique
,
Modeling Earth Systems and Environment
,
9
(
1
),
937
955
.
Antony Ravindran
A.
(
2012
)
Azimuthal square array resistivity method and goundwater exploration in Sanganoor, Coimbatore District, Tamilnadu, India
,
Res J Recent Sci ISSN
,
2277
,
2502
.
Bahraseman
S. E.
,
Firoozzare
A.
,
Zhang
C.
,
Yousefian
N.
,
Skominas
R.
,
Barati
R.
&
Azadi
H.
(
2024
)
Reviving the forgotten legacy: Strategies for reviving qanats as sustainable solutions for agricultural water supply in arid and semi-arid regions
,
Water Research
,
265
,
122138
.
Biswas
T.
,
Pal
S. C.
,
Ruidas
D.
,
Islam
A. R. M. T.
,
Saha
A.
,
Costache
R.
&
Shit
M.
(
2023
)
Modelling of groundwater potential zone in hard rock-dominated drought-prone region of Eastern India using integrated geospatial approach
,
Environmental Earth Sciences
,
82
(
3
),
81
.
Briscoe
J.
(
2005
)
India's Water Economy, Bracing for A Turbulent Future
.
Cameron
E.
&
Peloso
G.
(
2001
)
An application of fuzzy logic to the assessment of aquifers’ pollution potential
,
Environmental Geology
,
40
,
1305
1315
.
CGWB Master Plan for Artificial Recharge to Groundwater in India, Ministry of Water Resources, Govt. of India (2002)
.
Chandio
I. A.
,
Matori
A. N. B.
,
WanYusof
K. B.
,
Talpur
M. A. H.
,
Balogun
A. L.
&
Lawal
D. U.
(
2013
)
GIS-based analytic hierarchy process as a multicriteria decision analysis instrument: A review
,
Arabian Journal of Geosciences
,
6
,
3059
3066
.
Chitsazan
M.
,
Rahmani
G.
&
Neyamadpour
A.
(
2015
)
Forecasting groundwater level by artificial neural networks as an alternative approach to groundwater modeling
,
Journal of the Geological Society of India
,
85
,
98
106
.
Dar
T.
,
Rai
N.
&
Bhat
A.
(
2021
)
Delineation of potential groundwater recharge zones using analytical hierarchy process (AHP)
,
Geology, Ecology, and Landscapes
,
5
(
4
),
292
307
.
Das
B.
,
Pal
S. C.
,
Malik
S.
&
Chakrabortty
R.
(
2019a
)
Modeling groundwater potential zones of Puruliya district, West Bengal, India using remote sensing and GIS techniques
,
Geology, Ecology, and Landscapes
,
3
(
3
),
223
237
.
Das
U. K.
,
Roy
P.
&
Ghose
D. K.
(
2019b
)
Modeling water table depth using adaptive Neuro-Fuzzy Inference System
,
ISH Journal of Hydraulic Engineering
,
25
(
3
),
291
297
.
Dinka
M. O.
,
Loiskandl
W.
&
Ndambuki
J. M.
(
2015
)
Hydrochemical characterization of various surface water and groundwater resources available in Matahara areas, FantalleWoreda of Oromiya region
,
Journal of Hydrology: Regional Studies
,
3
,
444
456
.
Earls
J.
&
Dixon
B.
(
2007
). '
Spatial interpolation of rainfall data using ArcGIS: A comparative study
',
Proceedings of the 27th Annual ESRI International User Conference
, Vol.
31
, pp.
1
9
.
Ganapuram
S.
,
Kumar
G. T. V.
,
Krishna
I. V. M.
,
Kahya
E.
&
Demirel
M. C.
(
2009
)
Mapping of groundwater potential zones in the Musibasin using remote sensing data and GIS
,
Advances Engineering Software
,
40
,
506
518
.
Gnanachandrasamy
G.
,
Zhou
Y.
,
Bagyaraj
M.
,
Venkatramanan
S.
,
Ramkumar
T.
&
Wang
S.
(
2018
)
Remote sensing and GIS based groundwater potential zone mapping in Ariyalur District, Tamil Nadu
,
Journal of the Geological Society of India
,
92
,
484
490
.
Hou
J.
,
Zhang
Y.
,
Tong
Y.
,
Guo
K.
,
Qi
W.
&
Hinkelmann
R.
(
2019
)
Experimental study for effects of terrain features and rainfall intensity on infiltration rate of modelled permeable pavement
,
Journal of Environmental Management
,
243
,
177
186
.
Jain
S. K.
,
Agarwal
P. K.
,
Singh
V. P.
,
Jain
S. K.
,
Agarwal
P. K.
&
Singh
V. P.
(
2007
)
Problems related to water resources management in India
,
Hydrology and Water Resources of India
,
1
,
871
936
.
Jhariya
D. C.
,
Kumar
T.
,
Gobinath
M.
,
Diwan
P.
&
Kishore
N.
(
2016
)
Assessment of groundwater potential zone using remote sensing, GIS and multi criteria decision analysis techniques
,
Journal of the Geological Society of India
,
88
,
481
492
.
Joshi
P. K.
&
Tyagi
N. K.
(
1994
)
Salt affected and water logged soils in India: A review
,
Strategic Change in Indian Irrigation
,
237
252
.
Koohbanani
H.
,
Barati
R.
,
Yazdani
M.
,
Sakhdari
S.
&
Jomemanzari
R.
(
2018
)
Groundwater recharge by selection of suitable sites for underground dams using a GIS-based fuzzy approach in semi-arid regions
,
Progress in River Engineering & Hydraulic Structures
,
1
,
11
32
.
Kumar
R.
,
Singh
R. D.
&
Sharma
K. D.
(
2005
)
Water resources of India
,
Current Science
,
89
,
794
811
.
Kumar
P.
,
Herath
S.
,
Avtar
R.
&
Takeuchi
K.
(
2016a
)
Mapping of groundwater potential zones in Killinochi area, Sri Lanka, using GIS and remote sensing techniques
,
Sustainable Water Resources Management
,
2
,
419
430
.
Kumar
T.
,
Tripathi
M. P.
,
Katre
P.
&
Tiwari
P.
(
2016b
)
Groundwater recharge plan for Durg district of Chhattisgarh using satellite data and GIS technique
,
Agricultural Research Journal
,
53
(
2
),
234
242
.
Kuylenstierna
J. L.
,
Björklund
G.
&
Najlis
P.
(
1997
)
Future sustainable water use: Challenges and constraints
,
Journal of Soil and Water Conservation
,
52
(
3
),
151
156
.
Larondelle
N.
,
Haase
D.
&
Kabisch
N.
(
2014
)
Mapping the diversity of regulating ecosystem services in European cities
,
Global Environmental Change
,
26
,
119
129
.
Maghrebi
M.
,
Noori
R.
,
Partani
S.
,
Araghi
A.
,
Barati
R.
,
Farnoush
H.
&
Torabi Haghighi
A.
(
2021
)
Iran's groundwater hydrochemistry
,
Earth and Space Science
,
8
(
8
),
e2021EA001793
.
Maghrebi, M., Noori, R., Sadegh, M., Sarvarzadeh, F., Akbarzadeh, A. E., Karandish, F., Barati, H. & Taherpour, H.
(
2023
)
Anthropogenic decline of ancient, sustainable water systems: Qanats
,
Groundwater
,
61
(
1
),
139
146
.
Malczewski
J.
(
1999
)
GIS and Multicriteria Decision Analysis
.
New York, NY:
John Wiley & Sons
.
Mirboluki
A.
,
Mehraein
M.
,
Kisi
O.
,
Kuriqi
A.
&
Barati
R.
(
2024
)
Groundwater level estimation using improved deep learning and soft computing methods
,
Earth Science Informatics
,
17
(
3
),
2587
2608
.
Mohanty
S.
,
Jha
M. K.
,
Kumar
A.
&
Sudheer
K. P.
(
2010
)
Artificial neural network modeling for groundwater level forecasting in a river island of eastern India
,
Water Resources Management
,
24
,
1845
1865
.
Mukherjee
A.
(
2018
)
Overview of the Groundwater of South Asia
.
Singapore
:
Springer
, pp.
3
20
.
Nassery
H. R.
,
Adinehvand
R.
,
Salavitabar
A.
&
Barati
R.
(
2017
)
Water management using system dynamics modeling in semi-arid regions
,
Civil Engineering Journal
,
3
(
9
),
766
778
.
Neshat
A.
,
Pradhan
B.
&
Dadras
M.
(
2014
)
Groundwater vulnerability assessment using an improved DRASTIC method in GIS
,
Resources, Conservation and Recycling
,
86
,
74
86
.
Nouaceur
Z.
&
Murărescu
O.
(
2016
)
Rainfall variability and trend analysis of annual rainfall in North Africa
,
International Journal of Atmospheric Sciences
,
2016
,
1
12
.
Pande
C. B.
,
Moharir
K. N.
,
Panneerselvam
B.
,
Singh
S. K.
,
Elbeltagi
A.
,
Pham
Q. B.
, Varade, A. M. &
Rajesh
J.
(
2021
)
Delineation of groundwater potential zones for sustainable development and planning using analytical hierarchy process (AHP), and MIF techniques
,
Applied Water Science
,
11
(
12
),
186
.
Prasad
V. H.
&
Roy
P. S.
(
2005
)
Estimation of snowmelt runoff in Beas Basin, India
,
Geocarto International
,
20
(
2
),
41
47
.
Rajasekhar
M.
,
Sudarsana Raju
G.
,
Imran Basha
U.
,
Siddi Raju
R.
,
Pradeep Kumar
B.
&
Ramachandra
M.
(
2019
)
Identification of suitable sites for artificial groundwater recharge structures in semi-arid region of Anantapur District: AHP approach
,
Hydrospatial Analysis
,
3
(
1
),
1
11
.
Riley
S. J.
,
DeGloria
S. D.
&
Elliot
R.
(
1999
)
Index that quantifies topographic heterogeneity
,
Intermountain Journal of Sciences
,
5
(
1–4
),
23
27
.
Roozbahani
A.
,
Ebrahimi
E.
&
Banihabib
M. E.
(
2018
)
A framework for ground water management based on Bayesian network and MCDM techniques
,
Water Resources Management
,
32
(
15
),
4985
5005
.
Saaty
T. L.
(
1980
)
The analytic hierarchy process (AHP)
,
The Journal of the Operational Research Society
,
41
(
11
),
1073
1076
.
Saaty
T. L.
(
1984
)
The analytic hierarchy process: Decision making in complex environments. In Quantitative assessment in arms control: mathematical modeling and simulation in the analysis of arms control problems. Boston, MA: Springer US, pp. 285–308
.
Saaty
T. L.
(
2004
)
Decision making – the analytic hierarchy and network processes (AHP/ANP)
,
Journal of Systems Science and Systems Engineering
,
13
,
1
35
.
Sahu
R. K.
,
Mishra
S. K.
&
Eldho
T. I.
(
2010
)
An improved AMC-coupled runoff curve number model
,
Hydrological Processes
,
24
(
20
),
2834
2839
.
Saraf
A. K.
&
Choudhury
P. R.
(
1998
)
Integrated remote sensing and GIS for groundwater exploration and identification of artificial recharge sites
,
International Journal of Remote Sensing
,
19
(
10
),
1825
1841
.
Shah
T.
(
2008
)
India's master plan for groundwater recharge: An assessment and some suggestions for revision
,
Economic and Political Weekly
,
43
(
51
),
41
49
.
Shah
D.
&
Mishra
V.
(
2020
)
Integrated Drought Index (IDI) for drought monitoring and assessment in India
,
Water Resources Research
,
56
(
2
),
e2019WR026284
.
Singh
D. K.
&
Singh
A. K.
(
2002
)
Groundwater situation in India: Problems and perspective
,
International Journal of Water Resources Development
,
18
(
4
),
563
580
.
Singh
P.
,
Thakur
J. K.
&
Kumar
S.
(
2013
)
Delineating groundwater potential zones in a hard-rock terrain using geospatial tool
,
Hydrological Sciences Journal
,
58
(
1
),
213
223
.
Srdjevic
B.
&
Medeiros
Y. D. P.
(
2008
)
Fuzzy AHP assessment of water management plans
,
Water Resources Management
,
22
,
877
894
.
Stumpp
C.
,
Żurek
A. J.
,
Wachniew
P.
,
Gargini
A.
,
Gemitzi
A.
,
Filippini
M.
&
Witczak
S.
(
2016
)
A decision tree tool supporting the assessment of groundwater vulnerability
,
Environmental Earth Sciences
,
75
,
1
7
.
Tahmassebipoor
N.
,
Rahmati
O.
,
Noormohamadi
F.
&
Lee
S.
(
2016
)
Spatial analysis of groundwater potential using weights-of-evidence and evidential belief function models and remote sensing
,
Arabian Journal of Geosciences
,
9
,
1
18
.
Tortajada
C.
&
González-Gómez
F.
(
2022
)
Agricultural trade: Impacts on food security, groundwater and energy use
,
Current Opinion in Environmental Science and Health
,
27
,
1
7
.
Yeganeh
A.
,
Ahmadi
F.
,
Wong
Y. J.
,
Shadman
A.
,
Barati
R.
&
Saeedi
R.
(
2024
)
Shallow vs. deep learning models for groundwater level prediction: A multi-piezometer data integration approach
,
Water, Air, & Soil Pollution
,
235
(
7
),
1
24
.
Yeh
A. G. O.
&
Li
X.
(
2002
)
A cellular automata model to simulate development density for urban planning
,
Environment and Planning B: Planning and Design
,
29
(
3
),
431
450
.
Yeh
M. W.
&
Shellnutt
J. G.
(
2016
)
The initial break-up of Pangæa elicited by Late Palæozoic deglaciation
,
Scientific Reports
,
6
(
1
),
31442
.
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/).