Next Article in Journal
Workshop—Introduction
Previous Article in Journal
Physics of Shock and Physiological Effects on Biological Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Delineation of Groundwater Potential Zones Using Multi-Criteria Decision Analysis: The Case of Balkh Province, Northern Afghanistan †

by
Anayatullah Popalzai
1,
Hemayatullah Ahmadi
2,*,
Abdul Baqi Rahmani
1 and
Emrah Pekkan
2,3
1
Department of Geological Engineering and Exploration of Mines, Faculty of Geology and Mines, Kabul Polytechnic University, Kabul 1001, Afghanistan
2
Department of Remote Sensing and Geographic Information System, Institute of Graduate Programs, Eskişehir Technical University, Eskişehir 26555, Turkey
3
Institute of Earth and Space Sciences, Eskişehir Technical University, Eskişehir 26555, Turkey
*
Author to whom correspondence should be addressed.
Presented at the 4th International Electronic Conference on Geosciences, 1–15 December 2022; Available online: https://sciforum.net/event/IECG2022.
Proceedings 2023, 87(1), 41; https://doi.org/10.3390/IECG2022-14817
Published: 21 July 2023
(This article belongs to the Proceedings of The 4th International Electronic Conference on Geosciences)

Abstract

:
Without proper data and standard methods, assessing and determining groundwater resources in mountains and flat dry plains is difficult for a region. In this study, a GIS-based groundwater resource assessment is deemed a viable option, so the Analytical Hierarchy Process (AHP) is used to designate groundwater potential zones in the Balkh Province of northern Afghanistan. Several influential factors were provided for this purpose, including lithology, distance from the river, slope, drainage density, (LU/LC), lineament density, and rainfall. The groundwater potential zone map was grouped into four groundwater potential zones. They are very low (10.87%), low (35.13%), moderate (30.76%), and high (30.76%). The use of water level data was used to validate the results. The overall accuracy (71%) was found accordingly.

1. Introduction

Groundwater, often known as subterranean water, refers to all waters found beneath the earth’s surface [1]. Rainwater and/or snowmelt water seeps down through the soil and into the underlying rocks to generate groundwater, which is defined as water in the saturated zone [2,3]. It is one of the world’s most important natural resources, providing a key supply of water for towns, industry, agricultural outputs, and tourism activities [1,4]. Based on [5], there are 663 million people without access to safe water and 2.4 billion people without decent sanitation in the world. Every day, women and young girls altogether spend 266 million hours searching for a location to obtain water.
Groundwater occurrence and distribution are influenced by the climate and regional setting of the region, as well as surface and subsurface characteristics such as underlying rock fractures, land use type, geomorphic features, structural features, and their interactions with hydrological characteristics [6,7,8].
Afghanistan is situated in a semi-arid region with frequent precipitation fluctuations. The surface and groundwater are utilized for a variety of reasons, including drinking water, agriculture, and industry [9,10].
Due to low annual precipitation across the country and the disappearance of surface water, almost everyone now uses groundwater for multiple purposes, causing groundwater depletion and contamination, particularly in large cities like Kabul, Herat, and Balkh, which have grown in population in the last two decades [11,12,13,14].
At this time, remote sensing and geographic information systems (GISs) may be utilized to successfully manage groundwater and predict groundwater potential at the national, regional, and local levels [4].
According to the literature, several approaches have been widely used for the potentiality mapping of groundwater resources, including Multiple-Criteria Decision Analysis (MCDA) of Analytical Hierarchy Process (AHP), Frequency Ratio (FR), influencing factors combining Frequency Ratio, Index of Entropy, Certainty Factor (CF), etc. [4,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38]. Different influencing factors, such as geology, lithology, geomorphology, slope, land use, land cover, drainage density, lineament density, and rainfall, were evaluated in these works based on their weights and rates.
There has been relatively little research, notably in the hydrogeology sector, employing innovative technology on Afghan territory, particularly in the Balkh Province. Ref. [30] studied the paleochannel of the Balkh River (northern Afghanistan) and human occupation since the Bronze Age. In 2016, DACAAR carried out Vertical Electrical Sounding for groundwater assessment in Balkh’s new university, Nahr-I-Shahi District of Balkh Province. In recent years, ref. [32] has been evaluating the effects of climate change on field water demands using limited ground information in the Balkh Province. Recently, Shroder et al. [14] summarized the geology and water resources of Afghanistan and highlighted the main challenges in groundwater resources. Geospatial technology concerning the groundwater potentiality has thus far only been used in Kabul, the capital of Afghanistan; for example, ref. [33] assessed the potential zones by considering influencing factors, e.g., slope, geology, soil type, LU/LC, lineaments density, rainfall, and drainage density. Taking into account the present scenario in the Balkh Province, northern Afghanistan, and the absence of appropriate data and study on groundwater availability, which is the area’s primary concern, we chose to conduct this GIS-based study with the following key objectives: (1) to determine groundwater potential zones considering influencing factor using the Multiple-Criterial Decision-making Analysis (MCDA) of the Analytical Hierarchy Process (AHP) method, and (2) to uncover the effectiveness of AHP for the demarcation of groundwater reservoirs within the extent of the Balkh Province in northern Afghanistan.

2. Study Area

The study area (Balkh Province) lies in northern Afghanistan, at 36.76° N, 66.9° E, and 334 m above sea level. Balkh Province has a total area of 16,840 km2. Amudarya (Oxus River) is located in the province’s northwestern corner that separates Afghanistan from Tajikistan and Uzbekistan, with Kunduz Province in the east, Samangan and Sari Pul provinces in the south and southwest, and Jawzjan Province in the west (Figure 1).

3. Materials and Methods

In this study, geographic information systems and remote sensing along with other conventional approaches were used to map groundwater potential areas by applying the Multiple-Criteria Decision Making Analysis (MCDM) of the AHP model. Seven thematic layers, including lithology, slope, distance from the river, drainage density, land cover/land use, lineament density, and rainfall, were generated and weighted considering the expert knowledge and the previous literature. The entire procedure of the methodology is depicted in Figure 2.
All of the selected thematic layers have been rasterized and then projected into UTM Projection, WGS 84 Datum Zone 42, with a spatial resolution of 50 m. The slope and drainage density map were created using NASA’s ASTER GDEM with a 30 m resolution. A hydrology tool was used to prepare the drainage density.
The geology and lineament density data were obtained from the United States Geological Survey (USGS) website. The line density tool was also used to produce lineament density from the lineament map.
The LULC map for the study area was derived using the GIS environment from the FAO-Digital LULC Map incorporating Landsat 8 OLI/TIRS data. Following that, rainfall data for the research region was collected from a CRU (Climatic Research Unit) web satellite with a high resolution of 0.5°0.5°. Reclassification was completed when all of the selected thematic layers were converted to a rasterized format with the same cell size. Then, using the MCDA-AHP approach, weights were given and prioritized based on the literature and expert knowledge.
To delineate the final groundwater potential map for the study area, a Weighted Overlay Analysis (WOA) was conducted by combining all of the raster layers in the GIS environment.

3.1. Influential Factors

3.1.1. Lithology

The lithology of aquifer rocks affects their porosity and permeability and is a fundamental indicator of the quality and amount of groundwater in a specific location [4,37].
The lithology layer was derived from the Afghanistan geological map prepared by [39] and digitized using USGS. As illustrated in Figure 3A, the following four types of units occur within the study area:
  • Sandstone;
  • Conglomerate, siltstone, gypsum;
  • Shingly and detrital sediments, gravel, sand, clay, clay sand, loess, travertine;
  • Conglomerate, shingly sediments, loam, andesite, dacite, and their tuffs, phonolite.
ArcGIS 10.5 was used to reclassify and analyze the lithologic map for study (Figure 3A). Because of the importance of geology on groundwater potential, sedimentary and metamorphic rocks cover the majority of the study area. Furthermore, sedimentary rocks cover the central (Nahri Shahi, Balkh, Mazar-i-Sharif, Dihdadi, a part of Khulm, and Char Bolak) and northern (Kaldar, Dawlat Abad, and Shortepa) parts of the study area, which have the highest score and are most important for groundwater storage due to porosity such as pores, joints, fractures, bedding planes alluvium, and sandstone. Because of the lesser penetrability of the igneous rocks and mountainous areas in the south (Sholgara, Kishindih, Charkint, and Marmul), southeast (a section of Khulm), and southwest (a part of Chimtal), it is understood that groundwater availability is lower in these places.

3.1.2. Slope

The slope’s steepness directly influences the rate of rainfall infiltration and is thus used to assess groundwater possibilities. Due to insufficient time to percolate the surface, a steeper slope results in a higher surface runoff and less groundwater recharge. There is always an inverse connection between water infiltration and slope degree [40,41]. A mild and flat slope, on the other hand, always presents tremendous potential for groundwater recharging and has good groundwater development prospects.
The slope map in degrees for the study area, constructed using the ArcGIS 10.5 environment from the ASTER DEM 30 m imagery of the area (Figure 3B), reveals five slope classes: (1) 0–4.49, (2) 4.49–11.65, (3) 11.65–20.91, (4) 20.91–33.46, and (5) 33.46–76.18. According to Figure 3B, the lowest slope (0–4.49) is found in the northwest (Dawlat Abad and Charbolak) region of the province, indicating an essentially flat topography. In the south (Chimtal, Sholgara, Kishindih, Charkint, and Marmul) and in the east (a part of Khulm), the highest slope (33.46–76.18) appears in the shape of spread strips along the mountains. As a result, groundwater recharge is mostly likely in the northern section of the state (Kaldar, Dawlat Abad, Shortepa, and a part of Khulm).

3.1.3. Distance from River

Groundwater levels are significantly affected by the distance from the river. This appears to be a good impact, particularly in flat places [42,43]. The groundwater level is greater near the river than it is further away. When it comes to groundwater risk, the distance from water bodies might be subjective. As a result (Figure 3C), all locations within 500 m of the drainage line are classified as having a very high potential for groundwater recharge, whereas areas within 2000, 4000, 7000, 11,000, and >11,000 m of the river are classified as having a very high, high, moderate, low, and very low potential for groundwater recharge, respectively.

3.1.4. Lineament Density

Lineaments are linear or curved features on the Earth’s surface that indicate weaker bedrock zones and are used as a secondary aquifer in hard rock locations. They may be seen in remote sensing images [44,45].
Lineament density might reflect a region’s groundwater potentiality and permeability [46]. Geological features such as lineaments can be identified using remote sensing methods [44,47].
By taking into account the DEM of the study area, PCI Geomatic provided the lineament dataset. The resulting map was then resampled, reclassified, and plotted in the ArcGIS environment.
There are several intersecting lineaments in the study area. A groundwater potential zone is defined as the intersection of lineaments. The density of lineaments has been created and classified into five categories: extremely poor (0–0.0004), poor (0.0004–0.0008), moderate (0.0008–0.0013), high (0.0013–0.0017), and very high (0.0017–0.002) (Figure 3D). Because high-density lineaments are better for groundwater potential than low-density lineaments, higher weights are given to high-density lineaments, and lower weights are given to low-density lineaments.
The greatest concentration of lineament density is mostly in the southern (Marmul, Chimtal, Sholgara, and Kishindih) and southeastern regions (a part of Khulm) of the study area and these places are regarded to be a likely groundwater concentration zone. However, in the northeastern section, the density is very low or nearly non-existent, resulting in limited groundwater potentiality.

3.1.5. Drainage Density

Drainage density is inversely proportional to permeability; it is an important characteristic to consider when assessing the groundwater zone [48]. Also, the form and structure of the bedrock, the kind of vegetation, the rainfall absorption capacity of soils, infiltration, and slope gradient all influence an area’s drainage system [24].
ASTER DEM was used to derive the drainage density (Figure 3E). Similar to the lineament density, the drainage density is a quantitative quantity determined from the drainage pattern using comparable processes. The drainage density is defined as the ratio of the total stream lengths to the size of the grid area under consideration [9,37].
Using the equation below, the drainage density index was calculated:
D d = i = 1 i = n D i A   k m 1
where ∑Di is the total length of all streams i (km) and A is the area of the grid (km2).
As shown in Figure 3E, the resultant map for the area was divided into five classes: 0–0.034, 0.034–0.089, 0.089–0.15, 0.15–0.22, and 0.22–0.38. Low-drainage-density locations had higher ranks, whereas high-drainage-density areas received lower rankings, as seen in Figure 3E. Almost the whole northern part of the area had a high to extremely high drainage density. On the other hand, the poor drainage density covered the majority of the study area’s southern portion.

3.1.6. Rainfall

As part of the hydrological cycle, rainfall is another important component that affects groundwater potentiality [48]. This parameter greatly impacts the MCDA efficiency and groundwater potential [37]. Rainfall distribution and gradient are directly influenced by runoff penetration; this impacts the quantity of water that may percolate into the groundwater system, and it is feasible to improve groundwater potential [46]. As a result, rainfall is an important hydrologic component in our study. The rainfall dataset was retrieved for the research region from a CRU (Center for Hydrometeorology and Remote Sensing) web satellite with a high resolution of 0.5°0.5°.
The rainfall resulting map was classified into four categories, including 418.45–486.54 mm, 486.54–534 mm, 534–581.45 mm, 581.45–637.91 mm, and 637.91–763 mm. In comparison to the northern (Shortepa, Kaldar, Nahri shahi, and Shortepa), northwestern (Char Bolak and Dawlat Abad), and southern (Charkint, Marmul, Sholgara, Kishindih, and a part of Chimtal) parts of the study area, the eastern (a part of the Khulm) and southern (Charkint, Marmul, Sholgara, Kishindih) parts are shown in Figure 3F.

3.1.7. Land Use/Land Cover

The surface of the earth is described by land use and land cover. The natural layer comprises elements such as water bodies, forests, and so on, but the modified natural environment includes parameters such as croplands, fallow lands, plantations, and built-up areas, among others. These variables have an impact on groundwater either directly or indirectly.
The study area’s land use/land cover map shows six major groups of land use/land cover that were recognized and demarcated (Figure 3G). Satellite images may be used to interpret the diverse land use/land cover classifications of the research region, including agricultural, bare land, built-up land (town, village), forest land, water bodies, and sand cover.
Agricultural land almost completely surrounds the study area on all sides. In agricultural regions, runoff is typically lower, and penetrations are higher. However, the quantity of intrusion in settlement areas is usually lower. As a result, agricultural land and water bodies receive a high ranking due to their continuous recharge capacity. In contrast, built-up regions receive a poor ranking due to their negative influence on groundwater capacity.

3.2. Analytical Hierarchy Process (AHP)

AHP is one of the most appropriate decision-making techniques introduced by Saaty based on the multiple-criteria approach [49]. It assists decision-makers in resolving the spatial complexity of environmental challenges. The assignment of weights to each parameter and their normalization is a critical aspect in generating suitable results as the outcome is solely dependent on the weight assignments. However, there are numerous weight estimating strategies, but AHP is regarded as a promising methodology in groundwater modeling that may offer speedy, reliable, and cost-effective results [41].
AHP’s major features include hierarchical formulation, cost-effective and time-saving methods, and exact findings [26]. The key advantage of AHP is that it allows for expert-based factor ranking and ranking re-arrangement in the event of reduced accuracy.
The four steps of the AHP model are weight assignment, pair-wise comparison matrix, weight normalization, and consistency check. A total of seven criteria were chosen for this study to control the flow factor and storage capacity of groundwater in the study area [4,46,50,51,52,53].
The tables and final groundwater potential zones equation are calculated using the steps below (Table 1, Table 2, Table 3, Table 4 and Table 5).
Step 1: Pair-wise comparison matrix (PCM) computed using Equation (2):
X = X 11 X 12 X 1 n X 21 X 22 X 2 n X n 1 X n 1 X n n
where X is the pair-wise comparison matrix; Xnn is the indicator of the pair-wise matrix element.
Step 2: Normalizing the weights using Equation (3):
N W = G M n 1 N f G M n  
where NW is normalized weights; GMn is considered the geometric mean of the nth row of the pair-wise matrix (X).
Step 3: The consistency ratio (CR) is used to validate the pair-wise judgment matrix using Equation (4):
C R = C I R I
where CR is calculated by dividing CI (consistency index) by RCI (random consistency index) of Saaty [49] (Table 4).
Step 4: To calculate CI, the following Equation (5) was adopted:
C I = ( λ max n ) ( n 1 )
where λmax is the principal Eigen, and n represents the number of parameters selected for the study.
According to Saaty (1990) [49], a CR value of 0.10 is sufficient to continue the assessment. Furthermore, if the CR value is greater than 0.10, the analysis must be modified from the start to find the source of the matrix’s discrepancy. If the CR value is a perfect 0, then the PCM is absolutely accurate. However, the CR value in our study is less than 0.10, as calculated in Table 4, and hence the analysis can proceed.

3.3. Groundwater Potential Zone (GWPZ)

The groundwater potential zone (GWPZ) is a dimensionless variable that aids in predicting a region’s groundwater potential. The groundwater potential zone (GWPZ) was calculated using the weighted linear combination approach as follows [48]:
G W P Z = w = 1 m i = 1 n ( W i × X j )
where Wi is the normalized weight of the i thematic layer, Xj is the rank value of each class with respect to the j layer, m is the total number of thematic layers, and n is the total number of classes in a thematic layer. Equation (7) was used to determine the GWPZ for each cell:
G W P Z = i = 1 n A H P = L t W L t R + S l W S l R + D D W D D R + L C / L U W L C / L U R + L D W L D R + R f W R f R + R D W R D R
where Lt is the lithology, Sl is the slope, RD is the distance from the river, LC/LU is land cover/land use, respectively, DD is the drainage density, Rf is the rainfall, and LD is lineament density. While the subscripts W and R indicate final weights and normalized rates for each conditioning factor, respectively.
The quantile classification approach was used to divide GWPZ values into four categories: extremely low, low, moderate, and high (Figure 4). Each class in this classification technique has the same number of characteristics. Several researchers have also used the quantile approach because of its effectiveness [4,21].
Considering the geographical and environmental condition of the study area and the influence of the distance from the river, which has the second rank of weighing, we generated two groundwater potential zones based on the seven and six layers, excluding the distance-from-the-river layer as illustrated in Figure 5.

4. Results and Discussion

The major goal of this research is to determine groundwater potential mapping for the Balkh area in northern Afghanistan using the AHP model.
As a result, the study area’s low to very low groundwater potential zone is mostly found in the south and southeastern parts. High mountains with varied terrain, strong to severe slopes, and coarse to gravelly loam soil dominate these locations.
Furthermore, because the majority of precipitation falls as surface runoff, the area’s high relief and irregular topography result in limited infiltration capability. As a result, the groundwater level in this region of the study area is always low. As a result, this region of very low and low groundwater potential encompasses roughly 7644.2 km2 (46%). A moderate and high groundwater potential zone, taken together, covers an area of around 8973.3 km2 (54%).
The existence of high to moderate groundwater potential zones in the specified area might be influenced by the presence of high rainfall, gravel, and sandy gravel lithology, all of which indicate a good aquifer, linked lineament, lower slope, and low drainage density.
In fact, the presence of these structures in the study area will increase the likelihood of groundwater penetration and storage. The presence of considerable rains will amplify this effect. Moderate to low rainfall, lithology type (clay, shale, marl, and limestone), low lineament with no connections, steep slope, and high drainage density describes low to very low groundwater potential zones. As a result, there is less probability of groundwater penetration, and more rainfall is likely to flow as surface runoff.
Furthermore, agricultural fields in most locations have a high groundwater potential, indicating that the groundwater level in such areas is shallow. Geophysical studies and well logs might help determine the ideal locations for drilling groundwater discharge wells in this zone.
The area is divided into four groundwater potential zones based on seven thematic layers (lithology, slope, drainage density, lineament density, land use/land cover, distance from the river, and rainfall). They are very low (10.87%), low (35.13%), moderate (30.76%), and high(23.24%). The study area’s extremely good (high) potential zones are often found in the north, whereas low potential zones are found in the south, as summarized in Table 6.
The low and very low zones near the river are illustrated in Figure 4, and this result is due to the dryness of the river channel throughout the year. It is worth noting that after analyzing the layers, the distance from the river presents good results in those areas where the river channels have water throughout the year.
As a result, we ignored only the distance between layers when generating the groundwater potential area, instead of basing it on six thematic layers (lithology, slope, drainage density, lineament density, land use/land cover, and rainfall), and the result (Table 7) shows that the area is divided into four groundwater potential zones. They are extremely low (11.6%), low (19.9%), moderate (32.4%), and high (36%). The research area’s extremely excellent potential zones are often found in the north, whereas poor potential zones are found in the south. Meanwhile, the output of the created map, which is based on six thematic layers, is more accurate and dependable.
Validation is the most significant modeling procedure from the standpoint of scientific importance [4,54]. As a result, assessing the resulting GWPZ is critical. The groundwater well locations in the validation datasets were used to verify the GWPZ defined in this work.
The result was verified using data from 5776 monitoring wells collected over nine years (2005–2022) by the MRRD (Ministry of Rural Rehabilitation and Development) and AGS (Afghanistan Geological Department), with support from the USGS, DACAAR (Danish Committee for Aid to Afghan Refugees), and the Ministry of Energy and Water. The collected data were normalized based on the final result, and the validation procedure was applied to 5608 well data. Water levels in the monitoring wells ranged from 3.15 to 85.03 m on average. The mentioned wells were classified into three classes based on the static water level in the wells, 3.15–16.38 m, 16.38–31.46 m, and 31.46–85.03 m, which were attributed to high, moderate, and low groundwater levels, respectively, and inserted on the final GWPZ map. We then extracted the values for the selected well data and created a confusion matrix as described in Table 8 and Table 9. When it comes to tackling classification difficulties, the confusion matrix is a prominent tool. It can be used for both binary and multiclass classification issues.
According to six thematic layers, the total accuracy values of both layers were 71 percent and 54 percent, respectively. The accuracy evaluation of these approaches based on six thematic levels reveals a strong link between prospective groundwater zones and well record values. Consequently, the current study’s findings show the technique’s validity for groundwater potential zones in Afghanistan’s Balkh Province, and it may be extended to other locations with comparable geological settings.

5. Conclusions

For a long time, groundwater potential mapping has been carried out in this country by utilizing several traditional and remotely based methodologies. However, remote sensing technology and GIS studies have limitations, even though these studies make it simple and accessible for specialists to undertake prospective mapping at minimal prices and in a short amount of time.
In this study, the groundwater potential of different regions within the study area was predicted using the Analytical Hierarchy Process (AHP) model in the ArcGIS environment by considering seven thematic layers: lithology, slope, distance from the river, drainage density, land use/land cover, lineament density, and rainfall.
The weights and rankings assigned to each class and sub-class determine the consistency of the AHP model for groundwater potential delineation. As a result, understanding the elements that influence the targeted item is essential for improved results. According to the final groundwater potential map, 11.6% of the research region is classified as “very low potential”, 19.9% as “low potential”, 32.4 percent as “moderate”, and 35.9% as “high groundwater potential”.
High groundwater potential zones have few restrictions for intensive groundwater management. However, due to climate change and population growth, we must address the limits of urbanization and agricultural development in this area. For sustainable surface and groundwater management, regions categorized as moderate groundwater potential zones require good management planning addressing agriculture, plantation, and appropriate economic activities.
The areas with very low and low groundwater potential are largely covered by built-up areas, steep to extremely steep slopes, barren terrain, and thin soils with mountains, necessitating careful attention and planning for the management of surface and groundwater water resources. Artificial recharge packages, rainwater collection, the extension of agricultural areas, and the creation of forest areas are all suitable solutions for low and intermediate areas.
As a result, the groundwater potential map developed by this study’s robust AHP approach might be a useful resource for groundwater planning and management not only in this region, but also in agricultural areas. Furthermore, because the approach used in this study was efficient and precise, it could be used to perform groundwater exploration in diverse places with similar geologic, meteorological, and environmental characteristics.
The resulting maps show that the groundwater potential model created using AHP in GIS is good. The findings of this study may be used as a guide for responsible departments to perform correct groundwater management, distribution, and consumption planning in the future. Furthermore, to achieve reliable groundwater potential mapping at a wide scale over local districts and villages, the findings should be followed up with more extensive fieldwork and other related investigations.

Author Contributions

Conceptualization, A.P. and H.A.; methodology, A.P. and H.A.; software, A.P. and H.A.; validation, H.A., E.P. and A.B.R.; formal analysis, A.P.; investigation, A.P.; resources, H.A.; data curation, A.B.R.; writing—original draft preparation, A.P.; writing—review and editing, H.A.; visualization, H.A.; supervision, E.P.; project administration, H.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bear, J.J.; Cheng, H.-D.A. Groundwater and Aquifers. In Modeling Groundwater Flow and Contaminant Transport; Springer: Dordrecht, The Netherlands, 2010; pp. 65–80. [Google Scholar]
  2. Dunne, T. Hydrologyy mechanics, and geomorphic implications of erosion by subsurface flow. Spec. Pap. Geol. Soc. Am. 1990, 252, 1–28. [Google Scholar] [CrossRef]
  3. Adinolfi Falcone, R.; Falgiani, A.; Parisse, B.; Petitta, M.; Spizzico, M.; Tallini, M. Chemical and isotopic (δ18O‰, δ2H‰, δ13C‰, 222Rn) multi-tracing for groundwater conceptual model of carbonate aquifer (Gran Sasso INFN underground laboratory—Central Italy). J. Hydrol. 2008, 357, 368–388. [Google Scholar] [CrossRef]
  4. Ahmadi, H.; Kaya, O.; Babadagi, E.; Savas, T.; Pekkan, E. GIS—Based groundwater potentiality mapping using AHP and FR models in central Antalya, Turkey. Environ. Sci. Proc. 2021, 5, 8741. [Google Scholar]
  5. Graham, J.P.; Hirai, M.; Kim, S.S. An analysis of water collection labor among women and children in 24 sub-Saharan African countries. PLoS ONE 2016, 11, e0155981. [Google Scholar] [CrossRef] [PubMed]
  6. Jasrotia, A.S.; Kumar, R.; Saraf, A.K. Delineation of groundwater recharge sites using integrated remote sensing and GIS in Jammu district, India Delineation of groundwater recharge sites using integrated remote sensing and GIS in Jammu district, India. Int. J. Remote Sens. 2007, 28, 5019–5036. [Google Scholar] [CrossRef]
  7. Kumar, T.; Gautam, A.K.; Kumar, T. Appraising the accuracy of GIS-based Multi-criteria decision making technique for delineation of Groundwater potential zones. Water Resour. Manag. 2014, 28, 4449–4466. [Google Scholar] [CrossRef]
  8. Ibrahim-Bathis, K.; Ahmed, S.A. Geospatial technology for delineating groundwater potential zones in Doddahalla watershed of Chitradurga district, India. Egypt. J. Remote Sens. Sp. Sci. 2016, 19, 223–234. [Google Scholar] [CrossRef]
  9. Rahmati, O.; Nazari Samani, A.; Mahdavi, M.; Pourghasemi, H.R.; Zeinivand, H. Groundwater potential mapping at Kurdistan region of Iran using analytic hierarchy process and GIS. Arab. J. Geosci. 2015, 8, 7059–7071. [Google Scholar] [CrossRef]
  10. Flem, B.; Reimann, C.; Birke, M.; Banks, D.; Filzmoser, P.; Frengstad, B. Inorganic chemical quality of European tap-water: 2. Geographical distribution. Appl. Geochem. 2015, 59, 211–224. [Google Scholar] [CrossRef]
  11. Manap, M.A.; Nampak, H.; Pradhan, B.; Lee, S.; Sulaiman, W.N.A.; Ramli, M.F. Application of probabilistic-based frequency ratio model in groundwater potential mapping using remote sensing data and GIS. Arab. J. Geosci. 2014, 7, 711–724. [Google Scholar] [CrossRef]
  12. Goes, B.J.M.; Howarth, S.E.; Wardlaw, R.B.; Hancock, I.R.; Parajuli, U.N. Integrated water resources management in an insecure river basin: A case study of Helmand River Basin, Afghanistan. Int. J. Water Resour. Dev. 2016, 32, 3–25. [Google Scholar] [CrossRef]
  13. Azami, A.; Sagin, J.; Sadat, S.H.; Hejran, H. Sustainable Irrigation: Karez System in Afghanistan. Cent. Asian J. Water Res. 2020, 6, 65–86. [Google Scholar] [CrossRef]
  14. Shroder, J.F.; Eqrar, N.; Waizy, H.; Ahmadi, H.; Weihs, B.J. Review of the Geology of Afghanistan and its water resources. Int. Geol. Rev. 2021, 64, 1009–1031. [Google Scholar] [CrossRef]
  15. Ahmad, A.Y.; Al-Ghouti, M.A. Approaches to achieve sustainable use and management of groundwater resources in Qatar: A review. Groundw. Sustain. Dev. 2020, 11, 100367. [Google Scholar] [CrossRef]
  16. Esteban, E.; Dinar, A. Cooperative Management of Groundwater Resources in the Presence of Environmental Externalities. Environ. Resour. Econ. 2013, 54, 443–469. [Google Scholar] [CrossRef]
  17. Singh, A. Groundwater resources management through the applications of simulation modeling: A review. Sci. Total Environ. 2014, 499, 414–423. [Google Scholar] [CrossRef]
  18. Gaur, S.; Ch, S.; Graillot, D.; Chahar, B.R.; Kumar, D.N. Application of Artificial Neural Networks and Particle Swarm Optimization for the Management of Groundwater Resources. Water Resour. Manag. 2013, 27, 927–941. [Google Scholar] [CrossRef]
  19. Goldman, M.; Neubauer, F.M. Groundwater exploration using integrated geophysical techniques. Surv. Geophys. 1994, 15, 331–361. [Google Scholar] [CrossRef]
  20. Hasan, M.; Shang, Y.; Jin, W. Delineation of weathered/fracture zones for aquifer potential using an integrated geophysical approach: A case study from South China. J. Appl. Geophys. 2018, 157, 47–60. [Google Scholar] [CrossRef]
  21. Razandi, Y.; Pourghasemi, H.R.; Neisani, N.S.; Rahmati, O. Application of analytical hierarchy process, frequency ratio, and certainty factor models for groundwater potential mapping using GIS. Earth Sci. Inform. 2015, 8, 867–883. [Google Scholar] [CrossRef]
  22. Hossein, A.; Ardakani, H.; Ekhtesasi, M.R. Groundwater potentiality through Analytic Hierarchy Process (AHP) using remote sensing and Geographic Information System (GIS). JGeope 2016, 6, 75–88. [Google Scholar]
  23. Ozdemir, A. GIS-based groundwater spring potential mapping in the Sultan Mountains (Konya, Turkey) using frequency ratio, weights of evidence and logistic regression methods and their comparison. J. Hydrol. 2011, 411, 290–308. [Google Scholar] [CrossRef]
  24. Manap, M.A.; Sulaiman, W.N.A.; Ramli, M.F.; Pradhan, B.; Surip, N. A knowledge-driven GIS modeling technique for groundwater potential mapping at the Upper Langat Basin, Malaysia. Arab. J. Geosci. 2013, 6, 1621–1637. [Google Scholar] [CrossRef]
  25. Arnous, M.O.; El-Rayes, A.E.; Geriesh, M.H.; Ghodeif, K.O.; Al-Oshari, F.A. Groundwater potentiality mapping of tertiary volcanic aquifer in IBB basin, Yemen by using remote sensing and GIS tools. J. Coast. Conserv. 2020, 24, 1–20. [Google Scholar] [CrossRef]
  26. Das, S. Comparison among influencing factor, frequency ratio, and analytical hierarchy process techniques for groundwater potential zonation in Vaitarna basin, Maharashtra, India. Groundw. Sustain. Dev. 2019, 8, 617–629. [Google Scholar] [CrossRef]
  27. Zghibi, A.; Mirchi, A.; Msaddek, M.H.; Merzougui, A.; Zouhri, L.; Taupin, J.D.; Chekirbane, A.; Chenini, I.; Tarhouni, J. Using analytical hierarchy process and multi-influencing factors to map groundwater recharge zones in a semi-arid mediterranean coastal aquifer. Water 2020, 12, 2525. [Google Scholar] [CrossRef]
  28. Naghibi, S.A.; Pourghasemi, H.R.; Pourtaghi, Z.S.; Rezaei, A. Groundwater qanat potential mapping using frequency ratio and Shannon’s entropy models in the Moghan watershed, Iran. Earth Sci. Inform. 2015, 8, 171–186. [Google Scholar] [CrossRef]
  29. Kumari, A.; Kumar, D.; Warwade, P. Application of multi-criteria decision making (MCDM) and electrical resistivity tomography (ERT) techniques for identification of groundwater recharge zone(s) in granitic hard rock aquifer. J. Earth Syst. Sci. 2021, 130, 81. [Google Scholar] [CrossRef]
  30. Fouache, E.; Besenval, R.; Cosandey, C.; Coussot, C.; Ghilardi, M.; Huot, S.; Lamothe, M. Palaeochannels of the Balkh river (northern Afghanistan) and human occupation since the Bronze Age period. J. Archaeol. Sci. 2012, 39, 3415–3427. [Google Scholar] [CrossRef]
  31. Mahaqi, A.; Moheghi, M.M.; Mehiqi, M.; Moheghy, M.A. Hydrogeochemical characteristics and groundwater quality assessment for drinking and irrigation purposes in the Mazar-i-Sharif city, North Afghanistan. Appl. Water Sci. 2018, 8, 1–10. [Google Scholar] [CrossRef]
  32. Jami, A.; Tasumi, M.; Mosier, T.M.; Somura, H.; Konishi, T. Evaluation of the effects of climate change on field-water demands using limited ground information: A case study in Balkh province, Afghanistan. Irrig. Sci. 2019, 37, 583–595. [Google Scholar] [CrossRef]
  33. Nasir, M.J.; Khan, S.; Ayaz, T.; Khan, A.Z.; Ahmad, W.; Lei, M. An integrated geospatial multi-influencing factor approach to delineate and identify groundwater potential zones in Kabul Province, Afghanistan. Environ. Earth Sci. 2021, 80, 453. [Google Scholar] [CrossRef]
  34. Singh, S.K.; Noori, A.R. Delineation of groundwater recharge potential zones for its sustainable development utilizing GIS approach in Kabul basin, Afghanistan. Arab. J. Geosci. 2022, 15, 213. [Google Scholar] [CrossRef]
  35. Szeliga, W.; Bilham, R.; Kakar, D.M.; Lodi, S.H. Interseismic strain accumulation along the western boundary of the Indian subcontinent. J. Geophys. Res. Solid Earth 2012, 117, 8822. [Google Scholar] [CrossRef]
  36. Muradi, S.; Phien-Wej, N.; Giao, P.H. Depletion of Water Resources, Issues and Challenges of Water Supply Management in Mazar-i-Sharif City, Afghanistan. Res. J. Environ. Earth Sci. 2013, 5, 242–251. [Google Scholar] [CrossRef]
  37. Adiat, K.A.N.; Nawawi, M.N.M.; Abdullah, K. Assessing the accuracy of GIS-based elementary multi criteria decision analysis as a spatial prediction tool—A case of predicting potential zones of sustainable groundwater resources. J. Hydrol. 2012, 440–441, 75–89. [Google Scholar] [CrossRef]
  38. Mohammadi-Behzad, H.R.; Charchi, A.; Kalantari, N.; Nejad, A.M.; Vardanjani, H.K. Delineation of groundwater potential zones using remote sensing (RS), geographical information system (GIS) and analytic hierarchy process (AHP) techniques: A case study in the Leylia–Keynow watershed, southwest of Iran. Carbonates Evaporites 2019, 34, 1307–1319. [Google Scholar] [CrossRef]
  39. Abdullah, S.H.; Chmyriov, V.M.; Dronov, V.I. Geology and Mineral Resources of Afghanistan; British Geological Survey: Nottingham, UK, 2008; Volume 2. [Google Scholar]
  40. 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]
  41. 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]
  42. Kuriqi, A.; Ali, R.; Pham, Q.B.; Montenegro Gambini, J.; Gupta, V.; Malik, A.; Linh, N.T.T.; Joshi, Y.; Anh, D.T.; Nam, V.T.; et al. Seasonality shift and streamflow flow variability trends in central India. Acta Geophys. 2020, 68, 1461–1475. [Google Scholar] [CrossRef]
  43. 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. Ecol. Indic. 2021, 127, 107685. [Google Scholar] [CrossRef]
  44. Ahmadi, H.; Pekkan, E. Fault-Based Geological Lineaments Extraction Using Remote Sensing and GIS—A Review. Geosciences 2021, 11, 183. [Google Scholar] [CrossRef]
  45. Sharma, M.; Kumar, S.; Tiwari, M.; Shukla, Y.; Pandey, S.K.; Srivastava, P.; Banerjee, S. Palaeobiological constraints and the precambrian biosphere: Indian evidence. Proc. Indian Natl. Sci. Acad. 2012, 78, 407–422. [Google Scholar]
  46. Magesh, N.S.; Chandrasekar, N.; Soundranayagam, J.P. Delineation of groundwater potential zones in Theni district, Tamil Nadu, using remote sensing, GIS and MIF techniques. Geosci. Front. 2012, 3, 189–196. [Google Scholar] [CrossRef]
  47. Chaabouni, R.; Bouaziz, S.; Peresson, H.; Wolfgang, J. Lineament analysis of South Jenein Area (Southern Tunisia) using remote sensing data and geographic information system. Egypt. J. Remote Sens. Space Sci. 2012, 15, 197–206. [Google Scholar] [CrossRef]
  48. Shekhar, S.; Pandey, A.C. Delineation of groundwater potential zone in hard rock terrain of India using remote sensing, geographical information system (GIS) and analytic hierarchy process (AHP) techniques. Geocarto Int. 2015, 30, 402–421. [Google Scholar] [CrossRef]
  49. Saaty, T. Decision Making for Leaders: The Analytic Hierarchy Process for Decisions in a Complex World; RWS Publications: Pittsburgh, PA, USA, 1990. [Google Scholar]
  50. Agarwal, E.; Agarwal, R.; Garg, R.D.; Garg, P.K. Delineation of groundwater potential zone: An AHP/ANP approach. J. Earth Syst. Sci. 2013, 122, 887–898. [Google Scholar] [CrossRef]
  51. Saha, D.; Dhar, Y.R.; Vittala, S.S. Delineation of groundwater development potential zones in parts of marginal Ganga Alluvial Plain in South Bihar, Eastern India. Environ. Monit. Assess. 2010, 165, 179–191. [Google Scholar] [CrossRef]
  52. Prasad, R.K.; Mondal, N.C.; Banerjee, P.; Nandakumar, M.V.; Singh, V.S. Deciphering potential groundwater zone in hard rock through the application of GIS. Environ. Geol. 2008, 55, 467–475. [Google Scholar] [CrossRef]
  53. 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]
  54. Fabbri, A.; Chung, C.-J.F.; Fabbri, A.G. Validation of Spatial Prediction Models for Landslide Hazard Mapping. Nat. Hazards 2003, 30, 451–472. [Google Scholar] [CrossRef]
Figure 1. The geographical location of Balkh Province in northern Afghanistan: (A) location of Balkh Province in country-scale watersheds, (B) location of Balkh and surrounding provinces, and (C) the geomorphological settings of Balkh Province.
Figure 1. The geographical location of Balkh Province in northern Afghanistan: (A) location of Balkh Province in country-scale watersheds, (B) location of Balkh and surrounding provinces, and (C) the geomorphological settings of Balkh Province.
Proceedings 87 00041 g001
Figure 2. Flow chart of the methodology.
Figure 2. Flow chart of the methodology.
Proceedings 87 00041 g002
Figure 3. Thematic maps of study area: (A) lithology, (B) slope, (C) distance form river, (D) lineament density, (E) drainage density, (F) rainfall, and (G) land use/land cover.
Figure 3. Thematic maps of study area: (A) lithology, (B) slope, (C) distance form river, (D) lineament density, (E) drainage density, (F) rainfall, and (G) land use/land cover.
Proceedings 87 00041 g003
Figure 4. The groundwater potential zones map generated from the AHP model based on seven thematic layers for Balkh Province.
Figure 4. The groundwater potential zones map generated from the AHP model based on seven thematic layers for Balkh Province.
Proceedings 87 00041 g004
Figure 5. The groundwater potential zones map generated from the AHP model based on six thematic layers for Balkh Province.
Figure 5. The groundwater potential zones map generated from the AHP model based on six thematic layers for Balkh Province.
Proceedings 87 00041 g005
Table 1. Pair-wise comparison matrix for AHP.
Table 1. Pair-wise comparison matrix for AHP.
CriteriaLithologyDistance from RiverSlopeDrainage DensitylU/LCLineament DensityRainfallWeights
Lithology1.0002.0004.0005.0005.0007.0006.0000.355
Distance from river0.5001.0003.0006.0004.0007.0007.0000.281
Slope0.2500.3331.0002.0003.0004.0005.0000.137
Drainage density0.2000.1670.5001.0002.0003.0004.0000.091
LU/LC0.2000.2500.3330.5001.0002.0003.0000.067
Lineament density0.1430.1430.2500.3330.5001.0001.0000.036
Rainfall0.1670.1430.2000.2500.3331.0001.0000.034
Sum2.4604.0369.28315.08315.83325.00027.0001.000
Table 2. Normalized pair-wise matrix.
Table 2. Normalized pair-wise matrix.
CriteriaLithologyDistance from RiverSlopeDrainage DensitylU/LCLineament DensityRainfallWeights%
Lithology0.40660.49560.43090.33150.31580.280.22220.354635.4650
Distance from river0.20330.24780.32320.39780.25260.280.25930.280628.0560
Slope0.10160.08260.10770.13260.18950.160.18520.137013.7030
Drainage density0.08130.04130.05390.06630.12630.120.14810.09109.1030
LU/LC0.08130.06190.03590.03310.06320.080.11110.06676.6660
Lineament density0.05810.03540.02690.02210.03160.040.03700.03593.5880
Rainfall0.06780.03540.02150.01660.02110.040.03700.03423.4200
Sum1111111 100
Table 3. CI and CR worksheets.
Table 3. CI and CR worksheets.
CriteriaLithologyDistance from RiverSlopeDrainage DensitylU/LCLineament DensityRainfallWeighted Sum ValueSum/Weights
Lithology0.35460.56110.54810.45520.33330.25110.20522.70867.6375
Distance from river0.17730.28060.41110.54620.26660.25110.23942.17237.7427
Slope0.08870.09350.13700.18210.20000.14350.17101.01577.4124
Drainage density0.07090.04680.06850.09100.13330.10760.13680.65507.1947
LU/LC0.07090.07010.04570.04550.06670.07180.10260.47337.1000
Lineament density0.05070.04010.03430.03030.03330.03590.03420.25877.2124
Rainfall0.05910.04010.02740.02280.02220.03590.03420.24167.0664
λmax7.3380
Table 4. The consistency status.
Table 4. The consistency status.
Number of CriteriaRICIInconsistencyConsistent?
71.320.05630.0427Yes
Table 5. Weights and scores of criteria and their sub-criteria.
Table 5. Weights and scores of criteria and their sub-criteria.
No.FactorsSub-ClassesRatingNormalized RatesWeights
1LithologySandstone, siltstone, shale (Lagar and Argandab Tectonic Zones); limestone, sandstone, siltstone,30.3000.35465
Red sandstone, conglomerate, siltstone, gypsum, clay20.200
Shingly and detrital sediments, gravel, sand, clay, clay sand, loess, travertine40.400
Conglomerate, shingly sediments, loam; andesite, dacite, and their tuffs, phonolite and soevite tuff10.100
Sum10
2Distance from river<2258.3450.3330.2806
2258.34–4839.30040.267
4839.300–7742.88030.200
7742.880–11,453.0120.133
11,453.01–20,567.0210.067
Sum15
3Slope<4.48150.3330.137
4.481–11.65140.267
11.651–20.91330.200
20.931–33.46120.133
33.461–76.18510.067
Sum15
4Drainage density<0.034350.3330.091
0.0343–0.089540.267
0.0895–0.150630.200
0.1506–0.226720.133
0.2267–0.380310.067
Sum15
5LU/LCAgricultural land50.1520.0667
Bare Land40.121
Sand Cover60.182
Forest Land70.212
Bulit Up30.091
Waterbody80.242
Sum33
6Lineament density<0.193010.0670.0359
0.1930–0.524020.133
0.5240–0.820530.200
0.8205–1.123940.267
1.1239–1.758350.333
Sum15
7Rainfall418.4545–484.661310.0670.0342
484.6613–531.981820.133
531.9518–580.593530.200
580.5935–637.342240.267
637.3422–76350.333
Sum15
Table 6. The percentage of each potential zones according to seven thematic layers.
Table 6. The percentage of each potential zones according to seven thematic layers.
Province ClassesArea (km2)Area (%)
Balkhvery low1807.110.87
low5837.135.13
moderate5111.630.76
high3861.723.24
Table 7. The percentage of each potential zones according to six layers.
Table 7. The percentage of each potential zones according to six layers.
ProvinceClassesArea (km2)Area (%)
Balkhvery low192111.56
low332219.99
moderate539232.45
high598236.00
Table 8. The confusion matrix based on six thematic layers.
Table 8. The confusion matrix based on six thematic layers.
Well Data
Classified DataClassesHighModerateLowGrand TotalUser Accuracy
High2917414133345584%
Moderate5541026341192153%
Low46163142236%
Total351716034885608
Producer Accuracy83%64%3% 71%
Table 9. The confusion matrix based on seven thematic layers.
Table 9. The confusion matrix based on seven thematic layers.
ProvinceClassesArea (km2)Area (%)
Balkhvery low1807.110.87
low5837.135.13
moderate5111.630.76
high3861.723.24
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

Popalzai, A.; Ahmadi, H.; Rahmani, A.B.; Pekkan, E. Delineation of Groundwater Potential Zones Using Multi-Criteria Decision Analysis: The Case of Balkh Province, Northern Afghanistan. Proceedings 2023, 87, 41. https://doi.org/10.3390/IECG2022-14817

AMA Style

Popalzai A, Ahmadi H, Rahmani AB, Pekkan E. Delineation of Groundwater Potential Zones Using Multi-Criteria Decision Analysis: The Case of Balkh Province, Northern Afghanistan. Proceedings. 2023; 87(1):41. https://doi.org/10.3390/IECG2022-14817

Chicago/Turabian Style

Popalzai, Anayatullah, Hemayatullah Ahmadi, Abdul Baqi Rahmani, and Emrah Pekkan. 2023. "Delineation of Groundwater Potential Zones Using Multi-Criteria Decision Analysis: The Case of Balkh Province, Northern Afghanistan" Proceedings 87, no. 1: 41. https://doi.org/10.3390/IECG2022-14817

APA Style

Popalzai, A., Ahmadi, H., Rahmani, A. B., & Pekkan, E. (2023). Delineation of Groundwater Potential Zones Using Multi-Criteria Decision Analysis: The Case of Balkh Province, Northern Afghanistan. Proceedings, 87(1), 41. https://doi.org/10.3390/IECG2022-14817

Article Metrics

Back to TopTop