Next Article in Journal
Bottled Mineral Waters as Unconventional Sampling in Hydro-Geological Research
Previous Article in Journal
Investigation of the Internal Flow in a Francis Turbine for Comparing the Flow Noise of Different Operation Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of GIS, Multi-Criteria Decision-Making Techniques for Mapping Groundwater Potential Zones: A Case Study of Thalawa Division, Sri Lanka

by
Dilnu Chanuwan Wijesinghe
1,2,
Prabuddh Kumar Mishra
3,*,
Neel Chaminda Withanage
1,2,
Kamal Abdelrahman
4,
Vishal Mishra
5,
Sumita Tripathi
6 and
Mohammed S. Fnais
4
1
Department of Geography, Faculty of Humanities and Social Sciences, University of Ruhuna, Matara 81000, Sri Lanka
2
School of Geographical Sciences, Southwest University, Bebei District, Chongqing 400715, China
3
Department of Geography, Shivaji College, University of Delhi, New Delhi 110027, India
4
Department of Geology and Geophysics, College of Science, King Saud University, Riyadh 11451, Saudi Arabia
5
Helmholtz Centre Potsdam, Section 1.4 ‘Remote Sensing and Geoinformatics’, GFZ German Research Centre for Geosciences, Telegrafenberg, 14473 Potsdam, Germany
6
Department of Environment Studies, Shri Lal Bahadur Shastri National Sanskrit University, New Delhi 110016, India
*
Author to whom correspondence should be addressed.
Water 2023, 15(19), 3462; https://doi.org/10.3390/w15193462
Submission received: 8 August 2023 / Revised: 5 September 2023 / Accepted: 27 September 2023 / Published: 30 September 2023

Abstract

:
Groundwater resources are depleting due to phenomena such as significant climate change and overexploitation. Therefore, it is essential to estimate water production and identify potential groundwater zones. An integrated conceptual framework comprising GIS and the analytical hierarchy process (AHP) has been applied for the present study to identify groundwater potential areas in the Thalawa division of Sri Lanka. The criteria, including rainfall, soil types, slope, stream density, lineament density, geology, geomorphology, and land use, were taken into account as the most contributing factors when identifying the groundwater zones. Weights were allocated proportionally to the eight thematic layers according to their importance. Hierarchical ranking and final normalized weighting of these determinants were performed using the pairwise comparison matrix (PCM) available in AHP. Based on the results obtained, the groundwater potential zone (GWPZ) was classified into three regions: low potentiality (33.4%), moderate potentiality (55.8%), and high potentiality (10.6%). Finally, the zoning map was compared to find consistency with field data on groundwater discharge and depth taken from 18 wells in the division. The results revealed that the GIS-multi-criteria decision-making (MCDM) approach brings about noticeably better results, which can support groundwater resource planning and sustainable use in the research area.

1. Introduction

Water is an essential resource for every organism to maintain its life. Water that stays below the earth’s surface is called groundwater. Also, the term refers to all water found beneath the surface of the ground as groundwater [1,2]. The total volume of water resources and water on earth is estimated at 1386 trillion liters, of which 97% is saltwater, 3% is freshwater, and only 0.6% is available as groundwater [3]. It is crucial for ecology, food security, and human health [4]. Due to the scarcity of surface water, the significance of subsurface water reduction is highlighted, notably in dry areas [5]. With climate change and increased use of groundwater, many groundwater sources are experiencing water depletion. Groundwater availability depends on several factors, such as porousness, permeability, storage capacity, and transmissivity. The term “Groundwater potential” can be defined as the potential for groundwater to exist in an area [6,7].
The main indicator of groundwater storage reduction is the decrease in groundwater heads in wells. Groundwater in the Sri Lankan dry zone is mainly fed by rainfall. It should also be mentioned that, due to changes in rainfall patterns and excess water withdrawal for domestic and agricultural activities, the amount of water that reaches the ground is gradually decreasing. It is proven that 72% of Sri Lanka’s rural population and 22% of its urban residents use groundwater for drinking and domestic purposes [7]. The importance of a proper assessment of groundwater potential is indicated by studies on water scarcity and drought [8,9,10,11]. The potential for groundwater in Sri Lanka is lower compared to surface water resources [12]. According to the hydrogeological conditions, about 90% of the land area of the country consists of hard rocks with low potential for groundwater, and the rest of the land consists of sedimentary rocks with high groundwater potential [13]. Sedimentary rock is limited to the north, northwest, and northeast regions of Sri Lanka [13]. There is a significant need to use groundwater in a sustainable manner in the dry zone due to reduced rainfall and an increased population [14,15]. Thalawa is one of the most important agricultural regions in the dry zone that belongs to irrigated farming. In the study area, groundwater is currently obtained from shallow and very deep wells. Sustainable development of groundwater resources as the best option to support dry zone people can contribute to improving their well-being by increasing agricultural productivity without depleting groundwater resources.
We must make decisions in our daily lives; some may be complex, while others are simple. Multi-criteria decision analysis (MCDA) ranks potential actions or alternatives in order of priority. To make the best choice feasible, MCDA is used to assess and contrast several factors, which are frequently at odds with one another [16]. Multi-criteria approaches are referred to in the academic community in a variety of ways, including multiple-criteria decision aiding (MCDA), multi-criteria decision making (MCDM), multi-objective decision making (MODM), and multi-attribute decision making (MADM) [16,17,18]. Different MCDA techniques are available, including the analytical hierarchy process (AHP), TOPSIS, PROMETHEE, ELECTRE, SWARA, WASPAS, etc. However, the AHP, TOPSIS, and VIKOR approaches were the most popular techniques that have been used in previous research works during the last four decades [16,19]. The spatial dimension of the assessment criteria, decision substitutes, and geographic data models are all factors that MCDA takes into account while assessing the criteria [20,21]. Researchers frequently use AHP to support decision making in the processes of environmental planning and natural resource management because of its capacity to recognize and balance the significance of complex aspects [22,23,24,25]. Different approaches have been used in recent years to delineate groundwater potential zones globally. Many researchers have found in their studies that AHP and MCDA methods are effective tools for identifying groundwater potential. In the case of Sri Lanka, only a limited set of studies have been carried out to identify groundwater potential areas based on the geographic information system (GIS) [3,12,14,26,27]. Groundwater research has reached a turning point with the use of remote sensing (RS) and GIS in resource discovery, which has significantly aided in groundwater resource analysis, monitoring, and protection [28,29]. Abijita et al. [30] attempted to delineate potential groundwater zones in the Ponnaniyar Basin, Tamil Nadu, using AHP and the multi-influence factor (MIF). Detection of potential groundwater zones through an appropriate modeling approach was essential in solving water problems in the drought-prone Kilinochchi district [12,14]. Kumar et al. [31] have tried to identify groundwater potential zones in the Chennai river basin using GIS and AHP in their study. Rajasekhar et al. [32] identified groundwater potential areas in the Jiledubanderu River catchment, India, using GIS, AHP, and combined fuzzy-AHP techniques. A study was undertaken by Doke et al. [33] using a systematic and scientific GIS-based AHP to prepare a groundwater potential map for the Ulhas Basin, India. The concept of using the latest techniques, such as RS and GIS, for groundwater management research is relatively new [34,35,36,37]. As time- and cost-effective methods, combined GIS and AHP demonstrate it is a useful method for defining probable groundwater zones [38,39,40,41]. Hence, based on the literature, the GIS-MCDA integrated method has been used for the present study.
In the Anuradhapura district, there is not sufficient water to meet the required volume, and the potential of groundwater needs to be explored. No related studies have been carried out in this area previously. Therefore, the study has been conducted to fill the existing research gap in the study area. The resulting groundwater potential map will provide better insights into sustainable water resource management. The current study is primarily focused on using GIS MCDA methodologies to map the potential groundwater zones using eight different criteria. The research was structured into the five sections listed below. Section 1 is devoted to explaining the research background and previous literature, and Section 2 describes the study area, materials, and methods. The results are thoroughly explained in Section 3, along with the model validation. In the discussion section, the similarities and differences of the key findings are compared with other similar research works. Section 5 gives the conclusions.

2. Materials and Methods

2.1. Study Area

This study was carried out in the Thalawa Divisional Secretariat Division (DSD), located in the Anuradhapura district of the North Central Province of Sri Lanka, covering a 218.45 km2 area and having 39 Grama Niladhari Divisions (GNDs). It lies between 8°10′63″ N and 80°36′72″ E (Figure 1). The study area belongs to the dry zone; the average annual rainfall in the study area is around 1300 mm, with a mean annual temperature of 21 to 32 °C. There is 3182 ha of agricultural land in the Thalawa DSD under the agrarian service division of Thalawa, Eppawala, and Katiyawa [42]. The primary source of income in the study region is agriculture, with paddy being the predominant crop [43]. The area is part of the Mahaweli H zone, where irrigation water is used to support agriculture. All agricultural zones in the region are covered by the irrigation canal system. The northeast monsoon, which lasts from December to February, contributes significantly to the region’s annual rainfall. Based on the annual rainfall pattern, agriculture is carried out in two seasons: ‘Yala’ from March to September and ‘Maha’ from October to February. The tank cascades spread over the area are linked to agriculture; a substantial portion of the entire area’s topography is made up of very flat terrain with heights of less than 150 feet. The area has unique climatic characteristics, and except for a few months of the year, the other months are dry. There are a variety of plants adapted to dry climates. Reddish-brown earth soils, which are commonly found in the dry zone, are also available in the Talawa area. The lithological formation in the region is an important factor in the groundwater composition, quality, and formation of aquifers.

2.2. Selecting Criteria and Data Preparation

The selection of criteria for ranking is an important step in the suitability assessment process in any area for potential groundwater zones. Table 1 illustrates the criteria considered in earlier research when analyzing suitability for potential groundwater zones. Necessary data were obtained from various institutes, and then the steps required for ground-water potential zoning were followed by GIS MCDA [44].
The preparation of thematic layers (criteria) included RS data extraction, digitization of existing maps, and the collection of institutional data. Using preliminary investigation as a basis, the eight thematic layers were developed as follows: rainfall, geology, geomorphology, land use, soil type, stream density, lineament density, and slope [12,14,16,23,36]. Thematic layers were created once all the data were prepared, and these layers were then converted into raster datasets [35]. Finally, utilizing GIS-MCDA-integrated approaches, potential groundwater zones have been identified. The summary of the data sources is described in Table 2.
The methodological flowchart for the groundwater potential zones is illustrated in Figure 2. To map groundwater potentiality, the control factors for groundwater movement, storage, and occurrence may be investigated [62,63]. A rainfall map was produced using the inverse distance weighted (IDW) interpolation technique using point data sources employing Arc GIS 10.8 software. GNU Octave 7.3 software was used to calculate the criteria and factors for AHP weight using an algorithm developed by Mathew, 2020 (https://mathewmanoj.wordpress.com/mul (20 July 2023). Because the point data were scattered and sparse, IDW methods were selected instead of distance thresholding. The dependable IDW model is utilized in this work to interpolate geographical information based on the idea of weighting distance [12]. The geological and geomorphological layers of the study area were prepared using an existing map of the geological survey and mines bureau (GSMB) under the scale of 1:100,000. The soil map was created by digitizing the resource map with the help of the Irrigation Department. Shuttle radar topography mission (SRTM)—digital elevation model (DEM) was used to create the slope, stream, and lineament density layers of the area. Land use information was collected from the Survey Department of Sri Lanka. Finally, these thematic layers underwent raster data conversion and weighted overlay analysis in the Arc GIS environment.

2.3. Assignment of Weights and Criteria Normalization

In integrated analysis, the weight assignment of each feature class is most important since the weight assignment depends on the output. The AHP proposed by Saaty [63,64] was applied to allocate weight to each of the criteria used in the study. The AHP is a multi-criteria decision-making technique widely used in the field of groundwater studies as well as in environmental and other geospatial contexts [57]. Multi-criteria decision analysis can be identified as a commonly accepted and important technique for solving complicated problems [12]. Saaty’s scale was used in allocating standard weights [65]. Establishing weights for each criterion was the next stage. The weight computation determined the relative weights of the various criteria [66]. The AHP has the benefit of reducing pairwise comparisons of complicated judgments and aiding in determining the weight of the criterion [67,68].
For the chosen theme levels, a pairwise comparison matrix (PCM) was initially created to assess the scale weight of the relevant layers in relation to their contribution to groundwater potentiality. After that, a pairwise comparison matrix (M) was prepared, as in Equation (1). If n is the number of criteria, the size of M is n × n [61].
M = 1 p q 1 / p 1 r 1 / q 1 / r 1
Here, each component of M uses a number ranging from 1 to 9 (Table 3) to represent the relative weight of the two criteria [69]. Typically, the ratings range from 1 (equal importance) to 9 (extreme importance).
Based on earlier research and expert opinions obtained by referring a semi-structured questionnaire to six experts from different fields, the Saaty scale was used to assign weights to the selected criteria and determine their relationship with characteristics and influence on groundwater potential. Three GIS experts (University academics), one hydrologist, one geologist, and a land use planning director participated in the semi-structured questionnaire survey within their busy schedules and the limited pool of GIS and hydrology experts in the Sri Lankan context. The proportion of influence of the thematic layers and the categorization of the constraints were computed using the PCM, the relative weight matrix, and the normalized primary eigenvalue (Table 4).
The order of layer incentive on groundwater potential is expressed by the eigenvector [70]. The criteria of high groundwater potential are given more weight, whereas the criteria of low groundwater potential are given less weight. The column elements should be divided by the sum of the elements of the same column to normalize M and determine the weight of each criterion. The necessary relative test weights are provided by averaging the rows of the new matrix. Hence, all data are prepared as thematic layers and weighted overlay analysis using spatial analyst tools.

2.4. Normalized Weights and Identification Groundwater Potentiality

The relationship between the layers and their relative importance for the generation of the 8 × 8 pairwise matrix and the groundwater potential preparation determine how the eight thematic layers were integrated. They are displayed as, slope (SP), rainfall (RF), geomorphology (GM), soil (SL), geology (GL), stream density (SD), land use (LU), and lineament density (LD). Consequently, the result is a continuous mapping of suitability to produce a composite suitability map. The overlay tool creates raster layers using a common measurement range and weights each one based on its significance, which gives the final layer values ranges of 1–5 [12,34,71]. The last step was to prepare the composite map for groundwater potential. Weighted criteria are integrated to produce the potential map. This combination was performed by the weighted linear combination (WLC) method. The ability to achieve the relationship between the eight thematic maps using AHP with different classes is remarkable. Based on the PCM, the relative weight matrix and normalized weights were assigned to estimate the importance of the thematic layers on groundwater potential (Table 5).
The groundwater potential map was created by superimposing all of the criterion maps and utilizing the WLC technique with Arc GIS software 10.8 [71]. To calculate the groundwater potential index (GWPI), groundwater potential areas were produced using eight layers inserted into the GIS. The following formula describes the WLC process [72,73,74].
GWPI   =   w = 1 m 1 n   ( W j × X i )
Here, GWPI is the groundwater potential index, Wj is the normalized weight of the j-th thematic layer, Xi refers to the weight of the I class of the criteria, m represents the number of criteria, and n denotes the total number of classes. Table 6 illustrates the weights and ranks for each of the eight impact factors.
The relative weight (Wi) of the slope is 0.0732, and the ranking was fifth among all criteria. As the Thalawa area is on a flat surface belonging to the dry zone, very high rough, slope (deep) features cannot be identified. Rainfall became the most important criterion, gaining a high relative weight of 0.3109. Geomorphology was the second most important criterion in groundwater potential zoning, derived at 0.2231 relative weight. According to its significance, the soil criterion gained 0.1061 relative weight and became the fourth important factor. Geology was the third most important factor in groundwater potential zoning and derived 0.1650 relative weight from AHP. The land use criterion is the least important among all others according to the relative weight assigned by AHP derived from 0.0352. Stream density gained 0.0495 relative weight (sixth most important criterion) according to the AHP results. Lineament density gained a 0.0369 relative weight and was reported as the seventh most important factor in groundwater potential zoning. To perform GIS overlay analysis, each criterion was assigned a ranking order that ranges from 1–5 (1—very low, 2—low, 3—moderate, 4—high, 5—very high). These ranks were allocated based on the ranking orders collected from experts’ opinions using the questionnaire survey.

3. Results

3.1. Groundwater Potentiality for Major Criteria

The outputs of the AHP estimation of eight criteria and sub-criteria used in the research and standardized rating values (ri) are shown in Table 5 and Table 6. As well as the thematic layers that are produced through the reclassification of main criteria using ri values are shown in Figure 3. Details of all these criteria and their spatial distribution are described below.

3.1.1. Slope

The flat terrain (Figure 3a) of the study area may reveal the possibility of high groundwater potential. The slope of the area was categorized into five groups very high (0.018–1.5), high (1.6–3.9), moderate (4–8.3), low (8.4–15), and very low (16–22). Towards the center of the area, there is a slightly steeper slope. Due to the flat terrain and high infiltration rate, the research area with a slope of 0.018–1.5 is classified as having “very high” groundwater storage. The slope gradient between 16 and 22 has considered the groundwater storage as ‘very low’. While the very high-potential areas comprised 69% the high-potential areas were covered by 19.4%. Moderate, low, and very low were covered by 0.20%, 0.28%, and 7.9%, respectively.

3.1.2. Rainfall

Rainfall has a dominant effect on the hydrological cycle of the area and is directly related to the groundwater capacity. As a result of the dry climate of the study area, there are constant variations in the groundwater potential. The area receives rainfall with the activation of the northeast monsoon from September to December. The annual rainfall of the area has been classified into five classes (Figure 3b). The maximum and minimum rainfall in the area are 116 mm and 88 mm, respectively. Compared with the northern portion of the area, a higher trend of rainfall intensity can be detected in the southeast quarter. Rainfall distribution and slope gradient greatly influence surface water runoff, which also contributes to the determination of groundwater potential. As per the reclassification results, 67.6% of the area is of low potential, while moderate, high, and very high potential areas are covered by 18%, 9.5, and 4.7%, respectively.

3.1.3. Geomorphology

In the study area, two prominent landforms, namely, lower intermediate plantation surfaces and lower plant surfaces, have been observed. Conspicuous landforms were reclassified into two classes (Figure 3c). Porous and permeable zones are well explained in geomorphology and can be considered an essential phenomenon in groundwater recharge. High weight was allocated to the lower intermediate plantation surfaces, and low weight was assigned to the low plantation surfaces geomorphological unit that has moderate groundwater potential zones. 71.3% of the total area is in the moderate potential zone, while 28.7% is in the high potential zone.

3.1.4. Soil Types

Analysis of soil types revealed that the study area is mainly covered by four major soil types. Namely, alluvial, low-humic gley, red-yellow podzolic, and reddish-brown earth. This reddish-brown soil, which is typical of the dry zone, is spread over a large area. The brown texture and drainage of this soil are also widespread. The majority of the study area consists of reddish-brown earth soils (165 km2). Alluvial soils are covered in the south and northwest areas, and low-humic gley soils can be identified towards the center of the area. Red-yellow podzolic soils are spread over an area of 11.5 km2. In determining the influence of soil types on the occurrence of groundwater potential in the area, it can be identified that alluvial soils and low-humic gley soil contribute to being considered “very high” and “high”, respectively. Red-yellow podzolic soil was assigned a moderate weight because it was generally more conducive to stabilizing groundwater potential than reddish-brown soil. In depicting the groundwater potential, reddish-brown soil was assigned low weight due to its fine surface nature. When 74.1% of the area is in a low potential zone, high, very high, and moderate potential zones are covered at 13.7%, 6.7%, and 5.2%, respectively.

3.1.5. Geology

In the study area, the geological features are formed in relation to three types of rocks, and the determination of groundwater potential varies according to the geological types. Biotite Gneiss/Hornblende, Calciphyre/Minor Marble, Carbonatite, Charnockitic Gneiss, Granitic Gneiss with Pinkish Microcline, and Quartzite/Quartz Schist are the seven geological forms found in the area of study. Biotite Gneiss/Hornblende represents 47.5% and covers mainly the Northern and southwest parts, while other formations such as Granitic Gneiss with Pinkish Microcline, Quartz Schist, Carbonatite, and Calciphyre/Minor Marble are mostly identified in the Eastern portion of the area. Charnockitic Gneiss also represents 2.3% of the study area. The geology of the area indicates that the possible high groundwater holding formation is only Calciphyre/Minor Marble. Also, Calciphyre/Minor Marble is formed under the metamorphic rocks, which greatly influence the consistent groundwater capacity of the area. Calciphyre/minor marble is assigned high weights, while biotite gneiss/hornblende, carbonatite, charnokitic gneiss, and quartzite are given low weights, respectively. Granitic gneiss with pinkish Microcline, the second most widespread geological type in the area, has been found to have moderate suitability for groundwater potential. According to the thematic rating map (Figure 3e), only 1.4% of the area is in the high potential zone, while 42.9% and 51.6% are in the moderate and very low potential zones, respectively.

3.1.6. Land Use

Prominent land use/land covers include, paddy, homesteads, water bodies, forests, roads, and scrubs. Around 40% of the total area is rice paddy land. There is a tendency for water to infiltrate more in the vicinity of agricultural lands and forest areas. This is because a vegetated area has the ability to retain water and also allows for proper drainage. In the study, the water bodies and the forest cover have been identified as land covers that have a high groundwater potential. As the drainage system of the area is mainly connected with tanks and canals, water percolation occurs regularly. Therefore, it has been suggested in the study to give a higher weight to those land uses. Both homestead and paddy land uses have moderate suitability for groundwater potential (Figure 3f). The road network and scrubs are categorized as very low and low, respectively, for groundwater potential. After the rating assignment, 58.6% of the area was identified as having moderate potential, while 8.9% represented very high potential. Low and very low potential areas represented 26.4% and 4.1%, respectively.

3.1.7. Stream Density

The stream density has been categorized into five categories, and it varies in the range of <0.7 to >3.68 km2. A higher weight is given to regions with low stream density, and a lower weight is given to areas with very high stream density. In the study area, it is not possible to identify many regions with significantly high stream densities, and the stream density is more concentrated in the Southern part of the research area. Areas with low stream density have high levels of infiltration from rainfall. Hence, permeability works inversely and plays an important role in stream density for flow distribution and infiltration rates. Accordingly, 36.8% of the area (Figure 3g) is in the very high potential zone, while 22.9% is in the high potential zone. 5.2% and 12.1% of the area are in the very low and low potential zones, while the moderate potential zone is covered by 22.8%.

3.1.8. Lineament Density

The lineament map of the study area was derived from SRTM DEM using USGS. The lineament density of the study area varies from 0.7 km2 to 3.6 km2. In the study area, there is no specific place where the density of lines is high and the lines are distributed in a spreading manner. The linear features were classified into five groups, namely: 0–0.7 km2, 0.7–1.4 km2, 1.4–2.2 km2, 2.2–2.9 km2, and 2.9–3.6 km2. Based on the rating assignment for the thematic layer, 56.9% of the area was covered by a very low potential zone, while very high and high potential zones covered 17% together (Figure 3h).

3.2. Distribution of the Groundwater Potential Zones in Thlawa

In the study, every main criterion was individually rated and multiplied by AHP weights. Knowledge-based rating and weighting of various classes for each thematic layer have been assigned through the weighted overlay analysis process based on experts’ judgment. The weighted linear combination (WLC) output map displays three distinct groups, including low, moderate, and high potential for groundwater (Figure 4). High (23.34 km2), moderate (122.04 km2), and low (73.07 km2) groundwater potential zones cover about 10.6%, 55.8%, and 33.4% of the area, respectively. The factors of rainfall and geology have directly contributed to the presence of high groundwater potential conditions in the east, southeast, and south regions.

3.3. Validation of Groundwater Potential Zones in Thalawa

Any suitability assessment must be cross-checked with actual ground-truth information to maintain the reliability of the results. The delineated groundwater potential zone map was validated using secondary data and well-discharge statistics gathered during the field survey. To validate the resulting groundwater potential zone map, water discharge and depth data were integrated. Below-ground level (bgl) data were collected from the water resource board and Mahaweli Authority, Sri Lanka. Accordingly, the data was collected from 18 groundwater wells in three different potential zones (Table 7). Water discharge from ground wells ranges from 0.9 Ls−1 to 92 Ls−1. Also, the depth of the water level is between 0–77 m.
According to the validation results, four groundwater wells are located in the high potential area, while nine and five wells are located in moderate and low groundwater potential zones, respectively. The water discharge of the wells in the high potential zones is between 12 Ls−1 and 92 Ls−1 with a mean of 26 Ls−1 and the depth of groundwater ranges from 0.7 m to 4.2 m with a mean of 2.1 m. The mean water discharge of the wells in moderate potential zones is 44, with 14 Ls−1 minimum and 73 Ls−1 maximum. The mean groundwater depth of the wells in moderate zones is 21.2 m (min: 0.7 m max: 32 m). The water discharge of the wells in low potential zones ranges from 0.7 Ls−1 to 0.3 Ls−1 with a 0.92 Ls−1 mean. The mean water depth is 68 m (min: 37 m, and max: 77 m). The groundwater potential zone produced from the AHP approach demonstrated satisfactory levels of results when predicting the groundwater potential zone in Thalawa DS, Anuradhapura district, according to the verification results. The outcomes further demonstrated that groundwater potential zones might be identified using the methods presented here. Since the research area is essentially an agricultural-dominant DS, this will be more beneficial as a time-cost-efficient approach to choosing and finding groundswells for agricultural uses. This will reduce overexploitation and help conserve the area’s scarce groundwater resources. But for this purpose, spatial data at a finer spatial scale may increase the accuracy of the results. The study has limitations since we did not analyze the relationship between groundwater potential and well yield data.

4. Discussion

4.1. The Potentiality of Groundwater, Water Depth, and Discharge

The study area was divided into three distinct groundwater potential zones, namely high, moderate, and low. These zones encompassed 10.6%, 55.8%, and 33.4% of the total area, respectively. The Southern portion of the division is predominantly occupied by high-potential areas, while moderate-potential areas are scattered along the Eastern portion. Low-potential areas, on the other hand, are found in the Western and Northern parts of the division, as depicted in Figure 4. The acceptability of the derived groundwater potential map, generated using GIS-MCDA, was determined based on the groundwater depth and discharge data obtained from a set of representative ground wells. These data were utilized for validation. The existing body of research indicates a considerable correlation between groundwater potential zones and the discharge and water yield seen in the wells within the study locations [75,76,77].

4.2. The Impact of Slope, Rainfall, and Elevation on Groundwater Potential

The slope is a significant topographic element that influences the process of surface runoff. The impact of slope on groundwater recharge in aquifers is significant [78]. The relationship between slope gradient and surface water percolation is a significant factor in the delineation of groundwater potential zones [79,80]. Previous research has also indicated that flat terrain exhibits significant potential for groundwater accumulation as a result of the extended duration required for water to percolate [75,81,82,83].
The precipitation patterns exert a significant influence on the hydrological processes within the region and exhibit a direct correlation with the groundwater storage capacity [67,78]. The study area exhibits a maximum rainfall of 116 mm and a minimum rainfall of 88 mm. Several previous studies have demonstrated a strong positive association between rainfall and groundwater [75,84,85]. A study conducted on the Beshilo River watershed in Ethiopia revealed that precipitation had a positive impact on groundwater storage capacity. However, it was observed that the soil structure exhibited variations throughout the catchment area. In regions characterized by the presence of clay deposits, reduced rates of infiltration were found to correspond with decreased groundwater potential. In the context of sandy loam soil, it is observed that regions exhibit rapid absorption and subsequent contribution to the groundwater reservoir. In contrast to higher elevations, lower altitudes exhibited a higher likelihood of groundwater occurrence [86].

4.3. Effect of Geological Factors on Groundwater Potentiality

The utilization of geomorphology as a factor in various studies pertaining to the assessment of groundwater potentiality is of great significance. This is due to its ability to depict the landform and topography of a specific geographical region. Drought, classified as a natural calamity, manifests in regions characterized by arid climatic conditions, leading to the proliferation of vegetation patterns that have evolved to withstand such environmental constraints. The study area is characterized by low plantation surfaces, inselbergs, and a thin soil cover. The eroded relics observed in this context are identified as the oldest plain composed of Vijayan gneiss and quartzite [12]. The moderate and high potentiality of the bulk of the area can be attributed to the high permeability of water on low-level and intermediate plantation surfaces. Prior research has also demonstrated that areas characterized by rocks and sediments exhibiting a high degree of permeability experience rapid infiltration of water into the underlying soil [75,87].
Approximately 40% of the overall land surface consists of rice paddy fields, characterized by a higher degree of porosity that facilitates enhanced water percolation. Consequently, a significant portion of the land falls within the high and moderate potential classifications. Previous research has indicated a positive correlation between arable and agricultural land with moderate and high groundwater potential [75,88,89].
The stream density in a particular area governs the permeability and the percolation rate of precipitation. In the study area, the identification of locations with notably high stream densities that have had a discernible positive impact on water infiltration and movement, particularly in the Southern portion of the research area, is not feasible. Prior research has demonstrated that areas with low stream densities exhibit elevated groundwater potentials [75,90,91]. The hydrological processes of runoff and groundwater penetration in the area are significantly influenced by the presence of linear and curved structural elements [72,79,81]. Water movements have greater intensity in regions characterized by a higher density of lineaments. Previous investigations have also documented similar findings, indicating a favorable correlation between groundwater output and lineament features [75,92,93].
The spatial arrangement of groundwater within a given area is contingent upon the geological attributes of that region, which in turn impact the processes of infiltration and percolation [32,94,95]. Unconsolidated sediments, such as alluvium, exhibit a deficiency in the partitioning process, resulting in a significant proportion (55%) of the region being classified as zones with low and very low groundwater potential. Granitic gneiss, which ranks as the second most prevalent geological formation within the region under study, has been determined to possess a moderate level of appropriateness in terms of its potential for groundwater resources. Other studies utilizing GIS and multi-criteria MCDA have yielded comparable results in the mapping of groundwater potential zones [75,95]. The study area exhibits a limited extent of high water infiltration soil types, namely red-yellow podzolic soils and alluvial soils, which are found in very high and high potential zones, accounting for only 19% of the total area. In contrast, a significant portion of the study area, covering 165 km2, is characterized by reddish-brown earth with low infiltration capabilities, contributing to 74% of the area classified as having low potential for water infiltration.

5. Conclusions

The current study tries to delineate groundwater potential zones in the Thalawa division using GIS and AHP-MCDA techniques. Weights were assigned to eight potential criteria using AHP. A final groundwater potential map was generated by zoning the area into three classes: high, moderate, and low, using overlay analysis by the WLC method. The final map revealed that 10.6% of the study area has high groundwater potential, 55.08% has moderate potential, and 33.4% has low groundwater potential. Finally, the derived potential zone map was compared with the water discharge and depth data taken from representative groundwater wells in the study area. The mean discharge and mean depth of the groundwater wells in high-potential zones are 26 Ls−1 and 2.1 m. In the moderate zone, the mean discharge was 44 Ls−1 and the mean depth was 21.2 m, while ground wells in the low potential zone reported 0.92 Ls−1 mean discharge and 68 m mean depth. Thus, it revealed that the GIS-AHP integrated zoning map is acceptable and accurate. High groundwater potential areas are located in the Southern portion of the division and moderate potential zones are distributed over the Eastern and Southeastern portions. Low-potential areas are mainly distributed along the Western portion. Due to a lack of data, the relationship between groundwater potential and water yield was not used to validate the result. Incorporating those bivariate analyses into future research will be more helpful in deriving accurate and reliable results. This information will essentially support groundwater management planning decisions. The zoning map will provide policy guidelines for local planning authorities, especially for their agricultural water management. The study will pave the way for more advanced research that incorporates more criteria in national and international contexts.

Author Contributions

Conceptualization, D.C.W. and N.C.W.; Data curation, D.C.W.; Formal analysis, D.C.W.; Funding acquisition, K.A. and M.S.F.; Methodology, P.K.M. and S.T.; Resources, D.C.W., V.M., and M.S.F.; Software, D.C.W.; Visualization, P.K.M. and K.A.; Writing—original draft, D.C.W. and N.C.W.; Writing—review and editing, P.K.M. All authors have read and agreed to the published version of the manuscript.

Funding

Deep thanks and gratitude to the Researchers Supporting Project (RSP2023R249), King Saud University, Riyadh, Saudi Arabia, for funding the research article.

Data Availability Statement

Not Applicable.

Acknowledgments

Our sincere thanks go to the Survey Department, Geological Survey & Mines Bureau, and Meteorology Department of Colombo, Sri Lanka, which helped us succeed in the research. The authors also thankfully acknowledge USGS. Also, we are thankful to the anonymous reviewers, who helped us immensely to improve the quality of the manuscript through their constructive comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ibrahim-Bathis, K.; Ahmed, S.A. Geospatial technology for delineating groundwater potential zones in Doddahalla watershed of Chitradurga district, India. Egypt. J. Remote Sens. Space Sci. 2016, 19, 223–234. [Google Scholar] [CrossRef]
  2. Nampak, H.; Pradhan, B.; Manap, M.A. Application of GIS-based data driven evidential belief function model to predict groundwater potential zonation. J. Hydrol. 2014, 513, 283–300. [Google Scholar] [CrossRef]
  3. Senanayake, I.P.; Dissanayake, D.M.D.O.K.; Mayadunna, B.B.; Weerasekera, W.L. An approach to delineate groundwater recharge potential sites in Ambalantota, Sri Lanka using GIS techniques. Geosci. Front. 2016, 7, 115–124. [Google Scholar] [CrossRef]
  4. Gleeson, T.; Befus, K.M.; Jasechko, S.; Luijendijk, E.; Cardenas, C.M.B. The global volume and distribution of modern groundwater. Nat. Geosci. 2016, 9, 161–167. [Google Scholar] [CrossRef]
  5. Kumari, M.; Rai, S.R.; Mishra, P.K. Impact of Land-use/cover Changes on Groundwater Level Dynamics in Semi-Arid Region of India. Ann. Nat. Assoc. Geogr. India 2023, 43, 152–173. [Google Scholar] [CrossRef]
  6. Jha, M.K.; Chowdary, V.M.; Chowdhury, A. Groundwater assessment in Salboni Block, West Bengal (India) using remote sensing, geographical information system and multi-criteria decision analysis techniques. Hydrogeol. J. 2010, 18, 1713–1728. [Google Scholar] [CrossRef]
  7. Villholth, K.G.; Rajasooriyar, L.D. Groundwater Resources and Management Challenges in Sri Lanka—An Overview. Water Resour. Manag. 2010, 24, 1489–1513. [Google Scholar] [CrossRef]
  8. Yu, H.L.; Chu, H.J. Recharge signal identification based on groundwater level observations. Environ. Monit. Assess. 2011, 184, 5971–5982. [Google Scholar] [CrossRef] [PubMed]
  9. Chu, H.J. Drought detection of regional nonparametric standardized groundwater index. Water Resour. Manag. 2018, 32, 3119–3134. [Google Scholar] [CrossRef]
  10. Mishra, A.; Rai, A.; Mishra, P.K.; Rai, S.R. Evaluation of hydro-chemistry in a phreatic aquifer in the Vindhyan Region, India, using entropy weighted approach and geochemical modelling. Acta Geochim. 2023, 42, 648–672. [Google Scholar] [CrossRef]
  11. Pasalari, H.; Nodehi, R.N.; Mahvi, A.H.; Yaghmaeian, K.; Charrahi, Z. Landfill site selection using a hybrid system of AHP-fuzzy in GIS environment: A case study in Shiraz city, Iran. MethodsX 2019, 6, 1454–1466. [Google Scholar] [CrossRef] [PubMed]
  12. Pathmanandakumar, V.; Thasarathan, N.; Ranagalage, M. An Approach to Delineate Potential Groundwater Zones in Kilinochchi District, Sri Lanka, Using GIS Techniques. ISPRS Int. J. Geo-Inf. 2021, 10, 730. [Google Scholar] [CrossRef]
  13. Premathalika, K.M. Groundwater, Environment and Management. 2022. Available online: https://divaina.lk (accessed on 31 May 2023).
  14. Pankaj, K.; Herath, S.; Avtar, R.; Takeuchi, K. Mapping of groundwater potential zones in Killinochi area, Sri Lanka, using GIS and remote sensing techniques. Sustain. Water Resour. Manag. 2016, 2, 419–430. [Google Scholar]
  15. Alahacoon, N.; Edirisinghe, M. Spatial variability of rainfall trends in Sri Lanka from 1989 to 2019 as an indication of climate change. ISPRS Int. J. Geo-Inf. 2021, 10, 84. [Google Scholar] [CrossRef]
  16. Basílio, M.P.; Pereira, V.; Costa, H.G.; Santos, M.; Ghosh, A. A systematic review of the applications of multi-criteria decision aid methods (1977–2022). Electronics 2022, 11, 1720. [Google Scholar] [CrossRef]
  17. Sałabun, W.; Watróbski, J.; Shekhovtsov, A. Are MCDA Methods Benchmarkable? A Comparative Study of TOPSIS, VIKOR, COPRAS, and PROMETHEE II Methods. Symmetry 2020, 12, 1549. [Google Scholar] [CrossRef]
  18. Withange, W.K.N.C.; Gunathilaka, M.D.K.L.; Mishra, P.K.; Wijesinghe, W.M.D.C.; Thripathi, S. Indexing habitat suitability and human-elephant conflicts using GIS-MCDA in a human-dominated landscape. Geogr. Sustain. 2023, 4, 343–355. [Google Scholar] [CrossRef]
  19. Behzadian, M.; Otaghsara, S.K.; Yazdani, M.; Ignatius, J. A state-of the-art survey of TOPSIS applications. Expert Syst. Appl. 2012, 39, 13051–13069. [Google Scholar] [CrossRef]
  20. Greene, R.; Devillers, R.; Luther, J.E.; Eddy, B.G. GIS-based multiple-criteria decision analysis. Geogr. Compass 2011, 5, 412–432. [Google Scholar] [CrossRef]
  21. Wijesinghe, W.M.D.C.; Mishra, P.K.; Tripathi, S.; Abdelrahman, K.; Tiwari, A.; Fnais, M.S. Integrated flood hazard vulnerability modeling of Neluwa (Sri Lanka) using analytical hierarchy process and geospatial techniques. Water 2023, 15, 1212. [Google Scholar] [CrossRef]
  22. Malczewski, J. GIS-based multi-criteria decision analysis: A survey of the literature. Int. J. Geogr. Inf. 2006, 20, 703–726. [Google Scholar] [CrossRef]
  23. Mondal, B.K.; Kumari, S.; Ghosh, A.; Mishra, P.K. Transformation and risk assessment of the East Kolkata Wetlands (India) using fuzzy MCDM method and geospatial technology. Geogr. Sustain. 2022, 3, 191–203. [Google Scholar] [CrossRef]
  24. Jayasinghe, A.; Withanage, W. A geographical information system-based multi-criteria decision analysis of potato cultivation land suitability in Welimada divisional secretariat, Sri Lanka. Potato J. 2020, 47, 126–134. [Google Scholar]
  25. Wijesinghe, W.M.D.C.; Withanage, N.C. Application of the GIS-MCDA to identify land suitability for legume crops cultivation: An empirical investigationbased on Thalawa DSD, Sri Lanka. Int. J. Spat. Temporal Multimed. Inf. Syst. 2022, in press. [CrossRef]
  26. Kariyawasam, T.; Basnayake, V.; Wanniarachchi, S.; Sarukkalige, R.; Rathnayake, U. Application of GIS Techniques in Identifying Artificial Groundwater Recharging Zones in Arid Regions: A Case Study in Tissamaharama, Sri Lanka. Hydrology 2022, 9, 224. [Google Scholar] [CrossRef]
  27. Ramya, R.; Nanthakumaran, A.; Senanayake, I.P. Identification of artificial groundwater recharge zones in Vavuniya district using remote sensing and GIS. AGRIEAST J. Agric. Sci. 2019, 13, 44–55. [Google Scholar] [CrossRef]
  28. Shao, Z.; Huq, E.; Cai, B.; Altan, O.; Li, Y. Integrated remote sensing and GIS approach using Fuzzy-AHP to delineate and identify groundwater potential zones in semi-arid Shanxi. Environ. Model. Softw. 2020, 134, 104868. [Google Scholar] [CrossRef]
  29. Silwal, C.B.; Pathak, D. Review on practices and state of the art methods on delineation of ground water potential using GIS and remote sensing. Bull. Dep. Geol. 2018, 20, 7–20. [Google Scholar] [CrossRef]
  30. Abijitha, D.; Saravanan, S.; Singha, L.; Jennifer, J.J.; Saranya, T.; Parthasarathy, K.S.S. GIS-based multi-criteria analysis for identification of potential groundwater recharge zones—A case study from Ponnaniyaru watershed, Tamil Nadu, India. HydroResearch 2020, 3, 1–14. [Google Scholar] [CrossRef]
  31. Kumar, M.; Singh, P.; Singh, P.; Region, C. Fuzzy AHP based GIS and remote sensing techniques for the groundwater potential zonation for Bundelkhand Craton Region, India Fuzzy AHP based GIS and remote sensing techniques for. Geocarto Int. 2021, 37, 6671–6694. [Google Scholar] [CrossRef]
  32. Rajasekhar, M.; Sudarsana Raju, G.; Sreenivasulu, Y.; Siddi Raju, R. Delineation of groundwater potential zones in semi-arid region of Jilledubanderu river basin, Anantapur District, Andhra Pradesh, India using fuzzy logic, AHP and integrated fuzzy-AHP approaches. HydroResearch 2019, 2, 97–108. [Google Scholar] [CrossRef]
  33. Doke, A.B.; Zolekar, R.B.; Patel, H.; Das, S. Geospatial mapping of groundwater potential zones using multi-criteria decision-making AHP approach in a hardrock basaltic terrain in India. Eco. Indic. 2021, 127, 107685. [Google Scholar] [CrossRef]
  34. Chowdhury, A.; Jha, M.K.; Chowdary, V.M. Delineation of groundwater recharge zones and identification of artificial recharge sites in West Medinipur district, West Bengal, using RS, GIS and MCDM techniques. Environ. Earth Sci. 2010, 59, 1209–1222. [Google Scholar] [CrossRef]
  35. Mandal, B.; Dolui, G.; Satpathy, S. Land suitability assessment for potential surface irrigation of river catchment for irrigation development in Kansai watershed, Purulia, West Bengal, India. Sustain. Water Resour. Manag. 2018, 4, 699–714. [Google Scholar] [CrossRef]
  36. Berhanu, K.G.; Hatiye, S.D. Identification of groundwater potential zones using proxy data: Case study of Megech watershed, Ethiopia. J. Hydrol. Reg. Stud. 2020, 28, 100676. [Google Scholar] [CrossRef]
  37. Tamiru, H.; Wagari, M.; Tadese, B. An integrated Artificial Intelligence and GIS spatial analyst tools for Delineation of Groundwater Potential Zones in complex terrain: Fincha Catchment, Abay Basi, Ethiopia. Air Soil. Water Res. 2022, 15, 11786221211045972. [Google Scholar] [CrossRef]
  38. Andualem, T.G.; Demeke, G.G. Groundwater potential assessment using GIS and remote sensing: A case study of Guna tana landscape, upper blue Nile Basin, Ethiopia. J. Hydrol. Reg. Stud. 2019, 24, 100610. [Google Scholar] [CrossRef]
  39. Aslan, V.; Çelik, R. Integrated GIS-based multi-criteria analysis for groundwater potential mapping in the euphrates’s sub-basin, harran basin, turkey. Sustainability 2021, 13, 7375. [Google Scholar] [CrossRef]
  40. Achu, A.L.; Thomas, J.; Reghunath, R. Groundwater for Sustainable Development Multi-criteria decision analysis for delineation of groundwater potential zones in a tropical river basin using remote sensing, GIS and analytical hierarchy process (AHP). Groundw. Sustain. Dev. 2020, 10, 100365. [Google Scholar] [CrossRef]
  41. Sajil Kumar, P.J.; Elango, L.; Schneider, M. GIS and AHP Based Groundwater Potential Zones Delineation in Chennai River Basin (CRB), India. Sustainability 2022, 14, 1830. [Google Scholar] [CrossRef]
  42. Divisional Secretariat Division. Resource Profile—2020, Divisional secretariat Division: Thalawa, Anuradhapura, Sri Lanka. 2020. Available online: http://www.thalwa.ds.lk (accessed on 1 June 2023).
  43. Wijesinghe, W.M.D.C.; Withanage, W.K.N.C. Detection of the changes in land use and land cover using remote sensing and GIS in Thalawa DS Division. Prathimana J. 2021, 14, 72–86. [Google Scholar]
  44. Sarwar, A.; Ahmad, S.R.; Rehmani, M.I.A.; Asif Javid, M.; Gulzar, S.; Shehzad, M.A.; Shabbir Dar, J.; Baazeem, A.; Iqbal, M.A.; Rahman, M.H.U.; et al. Mapping groundwater potential for irrigation, by geographical information system and remote sensing Techniques: A case study of district lower Dir, Pakistan. Atmosphere 2021, 12, 669. [Google Scholar] [CrossRef]
  45. Pal, S.; Kundu, S.; Mahato, S. Groundwater potential zones for sustainable management plans in a river basin of India and Bangladesh. J. Clean. Prod. 2020, 257, 120311. [Google Scholar] [CrossRef]
  46. Verma, N.; Patel, R.K. Delineation of groundwater potential zones in lower Rihand River Basin, India using geospatial techniques and AHP. Egypt. J. Remote Sens. Space Sci. 2020, 24, 559–570. [Google Scholar] [CrossRef]
  47. Senthilkumar, M.; Gnanasundar, D.; Arumugam, R. Identifying groundwater recharge zones using remote sensing & GIS techniques in Amaravathi aquifer system, Tamil Nadu, South India. Sustain. Environ. Res. 2019, 29, 15. [Google Scholar] [CrossRef]
  48. Arulbalaji, P.; Padmalal, D.; Sreelash, K. GIS and AHP techniques based delineation of groundwater potential zones: A case study from southern western Ghats, India. Sci. Rep. 2019, 9, 2082. [Google Scholar] [CrossRef] [PubMed]
  49. Arefin, R. Groundwater for sustainable development groundwater potential zone identification at Plio-Pleistocene elevated tract, Bangladesh: AHP-GIS and remote sensing approach. Groundw. Sustain. Dev. 2020, 10, 100340. [Google Scholar] [CrossRef]
  50. Yıldırım, Ü. Identification of groundwater potential zones using GIS and multi-criteria decision-making techniques: A case study upper Coruh River Basin (NE Turkey). ISPRS Int. J. Geo-Inf. 2021, 10, 396. [Google Scholar] [CrossRef]
  51. Jhariya, D.C.; Mondal, K.C.; Kumar, T.; Indhulekha, K.; Khan, R.; Singh, V.K. Assessment of groundwater potential zone using GIS-based multi-influencing factor (MIF), multi-criteria decision analysis (MCDA) and electrical resistivity survey techniques in Raipur city, Chhattisgarh, India. Aqua Water Infrastruct. Ecosyst. Soc. 2021, 70, 375–400. [Google Scholar] [CrossRef]
  52. Benjmel, K.; Amraoui, F.; Boutaleb, S.; Ouchchen, M.; Tahiri, A.; Touab, A. Mapping of groundwater potential zones in crystalline terrain using remote sensing, GIS techniques, and multicriteria data analysis (Case of the Ighrem Region, Western Anti-Atlas, Morocco). Water 2020, 12, 471. [Google Scholar] [CrossRef]
  53. Singh, S.K.; Zeddies, M.; Shankar, U.; Griffiths, G.A. Potential groundwater recharge zones within New Zealand. Geosci. Front. 2019, 10, 1065–1072. [Google Scholar] [CrossRef]
  54. Tiwari, A.; Ahuja, A.; Vishwakarma, B.D.; Jain, K. Groundwater potential zone (GWPZ) for urban development site suitability analysis in Bhopal, India. J. Indian Soc. Remote Sens. 2019, 47, 1793–1815. [Google Scholar] [CrossRef]
  55. Pradhan, S.; Kumar, S.; Kumar, Y.; Sharma, H.C. Assessment of groundwater utilization status and prediction of water table depth using different heuristic models in an Indian interbasin. Soft Comput. 2019, 23, 10261–10285. [Google Scholar] [CrossRef]
  56. Metrereological Department of Sri Lanka. Annual Rainfall Data; Metrereological Department of Sri Lanka: Colombo, Sri Lanka, 2022.
  57. Geological Survey and Mines Bureau of Sri Lanka. Geological Map of Sri Lanka; Geological Survey and Mines Bureau of Sri Lanka: Sri Jayawardenepura Kotte, Sri Lanka, 1992.
  58. Geological Survey and Mines Bureau of Sri Lanka. Geomorphological Map of Sri Lanka; Geological Survey and Mines Bureau of Sri Lanka: Sri Jayawardenepura Kotte, Sri Lanka, 1982.
  59. Department of Irrigation. Soil Map of Sri Lanka; Department of Cartography, ITC: Colombo, Sri Lanka, 1988.
  60. Survey Department of Sri Lanka. Land Use Data; Survey Department of Sri Lanka: Colombo, Sri Lanka, 2022.
  61. USGS United States Geological Survey (USGS). Available online: https://earthexplorer.usgs.gov (accessed on 22 May 2023).
  62. Gunarathna, M.H.J.P.; Nirmanee, K.G.S.; Kumari, M.K.N. Geostatistical analysis of spatial and seasonal variation of groundwater level: A comprehensive study in Malwathu Oya cascade-I, I, Anuradhapura, Sri Lanka. Int. Res. J. Environ. Sci. 2016, 5, 29–36. [Google Scholar]
  63. Kolli, M.; Opp, C.; Groll, M. Mapping of potential groundwater recharge zones in the Kolleru Lake catchment, India, by using remote sensing and GIS techniques. Nat. Resour. 2020, 11, 127–145. [Google Scholar] [CrossRef]
  64. Tolche, A. Groundwater potential mapping using geospatial techniques: A case study of Dhungeta-Ramis sub-basin, Ethiopia. Geol. Ecol. Landsc. 2020, 5, 65–80. [Google Scholar] [CrossRef]
  65. Saaty, T.L. The Analytic Hierarchy Process: Planning, Priority Setting, Resource Allocation; McGraw-Hill: New York, NY, USA, 1980. [Google Scholar]
  66. Das, N.; Mukhopadhyay, S. Application of multi-criteriadecision making technique for the assessment of groundwater potential zones: A study on Birbhum district, west Bengal, India. Environ. Dev. Sustain. 2020, 22, 931–955. [Google Scholar] [CrossRef]
  67. Rahmati, O.; Samani, A.N.; Mahdavi, M.; Pourghasemi, H.R.; Zeinivand, H. Groundwater potential mapping at Kurdistan Region of Iran using analytic hierarchy process and GIS. Arab. J. Geosci. 2015, 8, 7059–7071. [Google Scholar] [CrossRef]
  68. Saaty, T.L. A scaling method for priorities in hierarchical structures. J. Math. Psychol. 1977, 15, 234–281. [Google Scholar] [CrossRef]
  69. Saaty, T.L. Fundamentals of the Analytic Network Process. In Proceedings of the International Symposium of the Analytic Hierarchy Process (ISAHP), Kobe, Japan, 12–14 August 1999. [Google Scholar]
  70. Saaty, T. Decision-making with the AHP: Why is the principal eigenvector necessary. Eur. J. Oper. Res. 2003, 145, 85–91. [Google Scholar] [CrossRef]
  71. Arya, S.; Subramani, T.; Karunanidhi, D. Delineation of groundwater potential zones and recommendation of artificial recharge structures for augmentation of groundwater resources in Vattamalaikarai Basin, South India. Environ. Earth Sci. 2020, 79, 102. [Google Scholar] [CrossRef]
  72. Malczewski, J. GIS and Multicriteria Decision Analysis; Wiley: Hoboken, NJ, USA, 1999. [Google Scholar]
  73. Suthakar, K. Assessment of land suitability potential for selected field crops using GIS based multicriteria analysis: Evaluating the case for Jaffna Peninsula, Sri Lanka’. Sri Lanka J. South Asian Stud. 2015, 1, 12–27. [Google Scholar]
  74. Sener, E.; Davraz, A.; Ozcelik, M. An integration of GIS and remote sensing in groundwater investigations: A case study in Burdur, Turkey. Hydrogeol. J. 2005, 13, 826–834. [Google Scholar] [CrossRef]
  75. Abrar, H.; Kura, A.L.; Dube, E.E.; Beyene, D.L. AHP based analysis of groundwater potential in the western escarpment of the Ethiopian rift valley. Geol. Ecol. Landsc. 2021, 7, 175–188. [Google Scholar] [CrossRef]
  76. Murmu, P.; Kumar, M.; Lal, D.; Sonker, I.; Singh, S.K. Delineation of groundwater potential zones using geospatial techniques and analytical hierarchy process in Dumka district, Jharkhand, India. Groundw. Sustain. Dev. 2019, 9, 100239. [Google Scholar] [CrossRef]
  77. Magesh, N.; Chandrasekar, N.; Soundranayagam, J. Delineation of groundwater potential zones in Theni district, Tamil Nadu, using remote sensing, GIS and MIF techniques. Geosci. Front. 2012, 3, 189–196. [Google Scholar] [CrossRef]
  78. Gupta, M.; Srivastava, P. Integrating GIS and remote sensing for identification of groundwater potential zones in the hilly terrain of Pavagarh, Gujarat, India. Water Int. 2010, 35, 233–245. [Google Scholar] [CrossRef]
  79. Satapathy, I.; Syed, T. Characterization of groundwater potential and artificial recharge sites in Bokaro district, Jharkhand (India), using remote sensing and GIS-based techniques. Environ. Earth Sci. 2015, 74, 4215–4232. [Google Scholar] [CrossRef]
  80. Qadir, J.; Bhat, M.S.; Alam, A.; Rashid, I. Mapping groundwater potential zones using remote sensing and GIS approach in Jammu Himalaya, Jammu and Kashmir. GeoJournal 2020, 85, 487–504. [Google Scholar] [CrossRef]
  81. Allafta, H.; Opp, C.; Patra, S. Identification of groundwater potential zones using remote sensing and GIS techniques: A case study of the Shatt Al-Arab Basin. Remote Sens. 2021, 13, 112. [Google Scholar] [CrossRef]
  82. Razavi-Termeh, S.; Sadeghi-Niaraki, A.; Choi, S. Groundwater potential mapping using an integrated ensemble of three Bivariate statistical models with random forest and logistic model tree models. Water 2019, 11, 1596. [Google Scholar] [CrossRef]
  83. Agarwal, R.; Garg, P.K. Remote sensing and GIS based groundwater potential & recharge zones mapping using multi-criteria decision making technique. Water Resour. Manag. 2016, 30, 243–260. [Google Scholar] [CrossRef]
  84. Lilienfeld, A.; Asmild, M. Estimation of excess water use in irrigated agriculture: A Data Envelopment Analysis approach. Agric. Water Manag. 2007, 94, 73–82. [Google Scholar] [CrossRef]
  85. Akale, A.T.; Dagnew, D.C.; Moges, M.A.; Tilahun, S.A.; Steenhuis, T.S. The effect of landscape interventions on groundwater flow and surface runoff in a watershed in the Upper Reaches of the Blue Nile. Water 2019, 11, 2188. [Google Scholar] [CrossRef]
  86. Kumar, V.; Mondal, N.; Ahmed, S. Identification of groundwater potential zones using RS, GIS and AHP techniques: A case study in a part of Deccan Volcanic Province (DVP), Maharashtra, India. J. Indian Soc. Remote Sens. 2020, 48, 497–511. [Google Scholar] [CrossRef]
  87. Yeh, H.F.; Cheng, Y.S.; Lin, H.I.; Lee, C.H. Mapping groundwater recharge potential zone using a GIS approach in Hualian River, Taiwan. Sustain. Environ. Res. 2016, 26, 33–43. [Google Scholar] [CrossRef]
  88. Fenta, A.A.; Kifle, A.; Gebreyohannes, T.; Hailu, G. Spatial analysis of groundwater potential using remote sensing and GIS-based multi-criteria evaluation in Raya Valley, northern Ethiopia. Hydrogeol. J. 2015, 23, 195–206. [Google Scholar] [CrossRef]
  89. Ansari, Z.R.; Rao, L.A.K.; Yusuf, A. GIS based morphometric analysis of Yamuna drainage network in parts of Fatehabad area of Agra district, Uttar Pradesh. J. Geol. Soc. India 2012, 79, 505–514. [Google Scholar] [CrossRef]
  90. Etikala, B.; Golla, V.; Li, P.; Renati, S. Deciphering groundwater potential zones using MIF technique and GIS: A study fromTirupati area, Chittoor District, Andhra Pradesh, India. HydroResearch 2019, 1, 1–7. [Google Scholar] [CrossRef]
  91. Magowe, M.; Carr, J.R. Relationship between lineaments and groundwater occurrence in western Botswana. GroundWater 1999, 37, 282–286. [Google Scholar] [CrossRef]
  92. Rashid, M.; Lone, M.A.; Ahmed, S. Integrating geospatial and ground geophysical information as guidelines for groundwater potential zones in hard rock terrains of south India. Environ. Monit. Assess. 2012, 184, 4829–4839. [Google Scholar] [CrossRef]
  93. Mondal, B.K.; Sahoo, S.; Das, R.; Mishra, P.K.; Abdelrahman, K.; Acharya, A.; Lee, M.-A.; Tiwari, A.; Fnais, M.S. Assessing groundwater dynamics and potentiality in the lower Ganga plain, India. Water 2022, 14, 2180. [Google Scholar] [CrossRef]
  94. Mishra, A.K.; Upadhyay, A.; Mishra, P.K.; Srivastava, A.; Rai, S.C. Evaluating geo-hydrological environs through morphometric aspects using geospatial techniques: A case study of Kashang Khad watershed in the Middle Himalayas, India. Quat. Sci. Adv. 2023, 11, 100096. [Google Scholar] [CrossRef]
  95. Ahmad, I.; Dar, M.A.; Andualem, T.G.; Teka, A.H. GIS-based multi-criteria evaluation of groundwater potential of the Beshilo River basin, Ethiopia. J. Afr. Earth Sci. 2020, 164, 103747. [Google Scholar] [CrossRef]
Figure 1. Geographical location of the study area: (a) Location of Sri Lanka in the Indian Ocean; (b) Location of Thalawa DSD; (c) Thalawa DSD in Landsat 8 false color (5,4,3) composite.
Figure 1. Geographical location of the study area: (a) Location of Sri Lanka in the Indian Ocean; (b) Location of Thalawa DSD; (c) Thalawa DSD in Landsat 8 false color (5,4,3) composite.
Water 15 03462 g001
Figure 2. Methodological flowchart for delineating groundwater potential zones.
Figure 2. Methodological flowchart for delineating groundwater potential zones.
Water 15 03462 g002
Figure 3. Distribution of rating values for eight criteria in groundwater potential zones in Thalawa: (a) Slope; (b) Rainfall; (c) Soil; (d) Geomorphology; (e) Geology; (f) Land use; (g) Stream density; (h) Lineament density.
Figure 3. Distribution of rating values for eight criteria in groundwater potential zones in Thalawa: (a) Slope; (b) Rainfall; (c) Soil; (d) Geomorphology; (e) Geology; (f) Land use; (g) Stream density; (h) Lineament density.
Water 15 03462 g003
Figure 4. Groundwater potential zones of the study area.
Figure 4. Groundwater potential zones of the study area.
Water 15 03462 g004
Table 1. Criteria used in previous studies to determine the potentiality of groundwater.
Table 1. Criteria used in previous studies to determine the potentiality of groundwater.
ReferencesRFGMGLSLSPLUDSLDASGLTWI
Ibrahim-Bathis and Ahmed [1]×× ××××
Pathmanandakumar et al. [12]××××××××
Aslan & Celik [39]××××××××
Kumar et al. [31]××××××××××
Sarwar et al. [44]××××××××
Pal et al. [45]× × ×× ××
Verma & Patel [46]××××××××
Senthilkumar et al. [47] × ××× ×
Arulbalaji et al. [48]×××××××× ×
Arefin [49]×××× ××
Yıldırım [50]×××××××× ×
Jhariya et al. [51]××××××××
Benjmel et al. [52] × ××× ×
Singh et al. [53] ×××××
Tiwari et al. [54]××××××××
Pradhan et al. [55]× ××× ×
Notes: RF: Rainfall, GM: Geomorphology, GL: Geology, SL: Soil, SP: Slope, LU: Land Use, DS: Drainage Density, LD: Lineament Density, AS: Aspect, GL: Groundwater Level, TWI: Topographic Wetness Index (TWI).
Table 2. Data sources used for mapping groundwater potential mapping.
Table 2. Data sources used for mapping groundwater potential mapping.
VariablesDataResolutionSource Locations
RainfallRainfall Data Department of Meteorology [56]
GeologyGeological map 1:100,000Geological Survey and Mines Bureau [57,58]
GeomorphologyGeomorphological map 1:100,000Shuttle Radar Topography Mission [58]
Soil typeSoil map 1:100,000Irrigation Department [59]
Land useLand use Data Survey Department of Sri Lanka [60]
SlopeShuttle Radar Topography Mission (SRTM)30 mUnited States Geological Survey [61]
Stream DensityShuttle Radar Topography Mission (SRTM)30 mUnited States Geological Survey [61]
Lineament DensityShuttle Radar Topography Mission (SRTM)30 mUnited States Geological Survey [61]
Table 3. AHP Scale.
Table 3. AHP Scale.
Scale123456789
ImportanceEquallyWeakModeratelyModerate PlusStrongStrong plusVery StrongVery very StrongExtreme
Table 4. Normalized pairwise comparison matrix analysis for the AHP Process.
Table 4. Normalized pairwise comparison matrix analysis for the AHP Process.
CriteriaSlopeRainfallGeomorphologySoilGeologyStream DensityLand UseLineament Density
Slope1.000.3330.20.3330.23.003.003.00
Rainfall3.001.005.003.005.005.003.005.00
Geomorphology5.000.21.005.003.005.005.005.00
Soil3.000.3330.21.000.25.005.005.00
Geology5.000.20.3335.001.005.003.005.00
Stream Density0.3330.20.20.20.21.003.001.00
Land use0.3330.3330.20.20.3330.3331.005.00
Lineament Density0.3330.20.20.20.21.000.21.00
Table 5. Normalized pairwise comparison matrix and weights were obtained for each criterion.
Table 5. Normalized pairwise comparison matrix and weights were obtained for each criterion.
CriteriaSPRFGMSLGLSDLULDNormalized Weight%
Slope (SP)0.060.120.030.020.020.120.130.100.07327.33%
Rainfall (RF)0.170.360.680.200.490.200.130.170.310931.09%
Geomorphology (GM)0.280.070.140.330.300.200.220.170.223122.31%
Soil (SL)0.170.120.030.070.020.200.220.170.106110.61%
Geology (GL)0.280.070.050.330.100.200.130.170.165016.51%
Stream Density (SD)0.020.070.030.010.020.040.130.030.04954.9%
Land use (LU)0.020.120.030.010.030.010.040.170.03523.53%
Lineament Density (LD)0.020.070.030.010.020.040.010.030.03693.7%
Table 6. Weight and ranking for different criteria.
Table 6. Weight and ranking for different criteria.
CriteriaWeightFeatureRank (ri)Potentiality Level
Slope
(Degrees)
70.018–1.55Very High
1.6–3.94High
4–8.33Moderate
8.4–152Low
16–221Very Low
Rainfall
(mm per month)
3188.9–94.42Low
94.5–99.92Low
100–1053Moderate
106–1114High
112–1165Very High
Geomorphology22Lower Levels of Intermediate Plantation Surfaces4High
Lower Plant Surfaces, Inselbergs, and thin Soil (Dry zone)3Moderate
Soil Types10Alluvial Soils5Very High
Low Humic Gley Soils4High
Red-Yellow Podzoilc Soils3Moderate
Redish-Brown Earth Soils2Low
Geology16Biotite Gneiss/Hornblende1Very Low
Calciphyre/Minor Marble4High
Carbonatite1Very Low
Charnockitic Gneiss2Low
Granitic Gneiss with Pinkish Microcline3Moderate
Quartzite/Quartz Schist2Low
Stream Density
(km2)
40–0.7355Very High
0.736–1.474High
1.48–2.213Moderate
2.22–2.942Low
2.95–3.681Very Low
Lineament Density (km2)30–0.3861Very Low
0.387–0.7722Low
0.773–1.163Moderate
1.17–1.544High
1.55–1.935Very High
Land use3Paddy3Moderate
Homestead3Moderate
Water Bodies4High
Forest3High
Road Network1Very Low
Scrubs2Low
Table 7. The data of verified sample wells in the Thalawa division.
Table 7. The data of verified sample wells in the Thalawa division.
ID of WellX and Y Coordinates (WGS 84)Depth of GW (m bgl)Water Discharge (Ls−1)GW Potential Zone
XY
GW1324,502164,5414.274.4High
GW2323,444164,2242.892High
GW3322,809162,2131.643.6High
GW4321,327160,8310.712High
GW5323,126157,34528.564.2Moderate
GW6324,396159,88524.673Moderate
GW7325,349161,8952958.3Moderate
GW8328,841163,8002134Moderate
GW9330,217164,4351954.7Moderate
GW10331,381164,0120.728.6Moderate
GW11331,910162,5301614Moderate
GW12338,366155,96914.329.5Moderate
GW13340,012147,1854.836.2Moderate
GW14329,370155,22872.848.6Low
GW15331,805154,80564.577Low
GW16331,180153,9587764.5Low
GW17335,297152,9003737Low
GW18335,297150,5714654.6Low
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wijesinghe, D.C.; Mishra, P.K.; Withanage, N.C.; Abdelrahman, K.; Mishra, V.; Tripathi, S.; Fnais, M.S. Application of GIS, Multi-Criteria Decision-Making Techniques for Mapping Groundwater Potential Zones: A Case Study of Thalawa Division, Sri Lanka. Water 2023, 15, 3462. https://doi.org/10.3390/w15193462

AMA Style

Wijesinghe DC, Mishra PK, Withanage NC, Abdelrahman K, Mishra V, Tripathi S, Fnais MS. Application of GIS, Multi-Criteria Decision-Making Techniques for Mapping Groundwater Potential Zones: A Case Study of Thalawa Division, Sri Lanka. Water. 2023; 15(19):3462. https://doi.org/10.3390/w15193462

Chicago/Turabian Style

Wijesinghe, Dilnu Chanuwan, Prabuddh Kumar Mishra, Neel Chaminda Withanage, Kamal Abdelrahman, Vishal Mishra, Sumita Tripathi, and Mohammed S. Fnais. 2023. "Application of GIS, Multi-Criteria Decision-Making Techniques for Mapping Groundwater Potential Zones: A Case Study of Thalawa Division, Sri Lanka" Water 15, no. 19: 3462. https://doi.org/10.3390/w15193462

APA Style

Wijesinghe, D. C., Mishra, P. K., Withanage, N. C., Abdelrahman, K., Mishra, V., Tripathi, S., & Fnais, M. S. (2023). Application of GIS, Multi-Criteria Decision-Making Techniques for Mapping Groundwater Potential Zones: A Case Study of Thalawa Division, Sri Lanka. Water, 15(19), 3462. https://doi.org/10.3390/w15193462

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

Article Metrics

Back to TopTop