1- Study Area
Nayband village is situated at approximately 29°54' longitude and 15°32' latitude on the western margin of the Lut Desert in eastern Iran, one of the southernmost points in Tabas County, South Khorasan. The study area is recognized as a desert region with a dry climate. This region is characterized by extreme temperature fluctuations, with high temperatures often exceeding 45°C in summer and dropping to around − 5°C in winter. Rainfall is generally scarce, occurring mostly in winter and early spring. In this area, there are no rivers with constant flow, but seasonal rains can lead to the formation of floods in the rivers. The study area is located at an elevation of approximately 1095 meters above sea level. Regarding the statistics of floods in this region, there is no available information. However, as part of the South Khorasan province, it should be noted that the rivers in this province mostly have seasonal flows, and the occurrence of occasional floods is minimal. The annual surface water flow volume is estimated to be around 400 million cubic meters. These water resources are utilized through traditional methods (South Khorasan Regional Water Authority, 2004). Also the primary source of water resources in this province comes from atmospheric precipitation, totaling over 12 billion cubic meters annually.
About 90% of the total precipitation, roughly 10.6 billion cubic meters, evaporates from the earth's surface, and only 1.4 billion cubic meters or 10% infiltrates in aquifers, contributing to groundwater recharge. Additionally, approximately 2 to 3% of the precipitation, equivalent to 400 million cubic meters, flows as surface water in rivers (South Khorasan Regional Water Authority, 2004). Furthermore, Table 1 and Table 2 present the details respectively regarding the Qantas and springs utilized in this study.
Table 1
Information on the utilized qanats for validation (South Khorasan Regional Water Company, 2022)
Average discharge (L/S)
|
year
|
name
|
city
|
Qanat
|
7.76
|
1398
|
Nayband
|
Tabas
|
Zardgah Nayband
|
0.09
|
1398
|
Nayband
|
Tabas
|
Ab Reza
|
3.51
|
1398
|
Nayband
|
Tabas
|
Dehnonayband
|
Table 2
Information on the utilized springs for validation (South Khorasan Regional Water Company, 2022)
Average discharge (L/S)
|
year
|
name
|
city
|
spring
|
2.275
|
1398
|
Nayband
|
Tabas
|
Dig Rostam2
|
2.15
|
1398
|
Nayband
|
Tabas
|
Dig Rostam1
|
0.045
|
1398
|
Nayband
|
Tabas
|
Ab Reza
|
2- Data
To conduct the site selection process, it is necessary to first collect diverse spatial information from various sources and then organize them in the form of different information layers within the coverage of GIS-related software. Subsequently, using GIS tools and techniques, analyze this spatial information. These steps ensure that the data is interpreted accurately and logically. Then, based on the results of the analyses, select and determine suitable locations. This method, aided by spatial data and generated maps, contributes to improving decision-making and area management.
In order to prepare the necessary information layers for this research, satellite images of the area were initially collected through the Landsat system. These satellite images then underwent pre-processing stages, including geometric and radiometric corrections. After performing these corrections and pre-processing steps using ENVI software, satellite data were processed and analyzed to extract patterns and information required for the research. This process, from obtaining satellite images to processing data to obtain spatial information and patterns, ultimately contributed to the research.
In this study, remote sensing data were collected from various sources. ETM images with a resolution of 30 meters and SRTM images were the primary sources used. Additionally, geological maps of the region were used to extract information related to lithology and fault lines. SRTM images were also used to create land-use maps. SRTM data were considered as the base data in the research. Slope and topography analyses were obtained from the Digital Elevation Model (DEM) and SRTM Digital Surface Model (DSM) maps.
In this research, eight factors influencing the determination of groundwater potential have been utilized including:
- Slope of the land.
- Geomorphology: Geological features of the land surface.
- Soil: Type and characteristics of the soil in the region.
- Lineament density (Fault lines): The presence of fault lines or land fractures.
- Land cover: Type and vegetation cover in the area.
- Precipitation: Amount and pattern of precipitation in the region.
- Aquifer density and proximity to surface water: Proximity to rivers or surface water sources.
These factors can have a significant impact on the potential of groundwater resources, and based on the analysis of these factors, maps of groundwater potential have been prepared.
The density of drainage patterns and proximity to surface water:
In studies related to aquifers and groundwater resources, the characteristics of drainage patterns, especially permeability, are of paramount importance. Drainage patterns play a crucial role in nourishing groundwater resources. The density and features of drainage patterns have a significant impact on aquifer permeability. Increasing the density of drainage patterns leads to improved aquifer permeability, signifying enhanced ability for aquifers to draw from sedimentary water. Therefore, analyzing drainage patterns can contribute to a better understanding and management of groundwater resources.
Additionally, a detailed examination of drainage patterns in different regions improves our understanding of the hydrogeological characteristics of groundwater resources, providing more precise information about permeability, yield, and other aquifer properties. This information can be effective in planning and optimizing the management of groundwater resources for various purposes.
In reality, the drainage pattern of a region reveals important features of the land surface and subsurface structures. The competence of groundwater potential is associated with the drainage pattern, as it indicates the distribution of surface runoff and land permeability. The density of the drainage pattern can significantly affect the yield of aquifers. Drainage patterns contribute to the proper distribution and management of surface water resources, guiding rainwater towards aquifers and groundwater. This interaction between drainage patterns and groundwater potential highlights the importance of detailed analysis of drainage patterns and management measures in different regions (Manikandan, Kiruthika, & Sureshbabu, 2014). In this analysis, in addition to drainage pattern density, the distance from drainage patterns has been examined as an essential parameter. Using the "Distance" command in ArcGIS software, the study area has been classified into five categories based on the distance from drainage patterns, including very close, close, moderate, weak, and very weak. This classification reflects the distance of different areas from drainage patterns and their proximity or distance to drainage pattern routes, serving as a basis for analyzing groundwater potential.
Land Slope:
The slope of the land plays a crucial role in controlling surface runoff, groundwater recharge, and the movement of surface water. Indeed, this layer is recognized as one of the fundamental factors influencing groundwater recharge and water movement in different regions. The topography of an area, especially the slope of the land surface, can significantly impact the permeability of water in soils and rocks.
Areas with lower slopes are typically classified as highly suitable for groundwater recharge because a gentle slope can enhance the permeability of soils and rocks, allowing surface water to infiltrate into groundwater. Conversely, areas with moderate slopes are often considered suitable for groundwater recharge due to their relatively stable topography and typically good permeability.
Regions with steeper slopes are generally classified as less favorable for groundwater recharge. These areas have relatively high surface runoff and lower infiltration of surface water into groundwater, resulting in less groundwater storage. Therefore, land slope can be crucial in terms of its impact on groundwater recharge and surface water flow in the region (Elubid et al., 2020). In this study, Digital Elevation Model (DEM) data in the Arc Map software was utilized to generate slope maps for the region.
Geomorphology:
Geomorphology and the identification of surface features in a region are of special importance in studying groundwater and karsts. Geomorphology refers to the shapes and structures of the Earth's surface that can have a significant impact on the flow and storage of groundwater. Satellites and satellite images can provide useful information about these features (Asgari Moghaddam et al., 2007). By analyzing satellite images and identifying surface features in the area, regions with high groundwater potential can be determined. Features that may be examined in satellite images include elevations, shape and nature of the land surface, the presence of karst structures, changes in vegetation density, and rock and sediment density.
Using this information, researchers can identify areas with favorable conditions for groundwater storage and recharge and utilize these resources for exploration and optimal management. This data can contribute to improving water resource plans and decision-making processes (NRSC, 2010). In this study, the region has been divided into three categories using the geomorphology map:
- Plains: These areas are typically flat with a low to nearly zero slope. The low slope can make these areas suitable for groundwater recharge.
- Grasslands with moderate slope: This category includes areas with a moderate slope, usually grasslands and meadows. These areas may also be suitable for groundwater recharge.
- Lands with steep slope: This category includes areas with a steep and mountainous slope. These types of lands are generally not suitable for groundwater recharge because a steep slope can lead to surface runoff and a decrease in the density of water channels.
With this classification, more information about the geomorphological features of the study areas is available, and this information can be useful in potential assessment and groundwater resource management.
Soil:
Soil is a fundamental parameter in groundwater potential and can significantly impact permeability and groundwater recharge. In this regard, two important soil characteristics are discussed as follows:
Permeability: Permeability, or soil transmissivity, indicates the soil's ability to allow water to pass and infiltrate underground. Soils with high permeability can easily feed water to aquifers. Therefore, soils with high permeability are generally considered good for groundwater recharge.
Lineament density (Faults):
In summary, the presence of faults and fault zones can play a crucial role in the distribution and recharge of groundwater. These areas may have a high density of faults, contributing to the accumulation of groundwater. This information can be utilized in planning and managing groundwater resources in the region (Haridas, Aravindan, & Girish, 1998).
These maps and analyses aid in understanding the factors influencing surface water infiltration and groundwater distribution. The density of faults and the rock fracture system are among the critical factors that can play a fundamental role in the recharge and distribution of groundwater. The information derived from these analyses assists water resource managers and decision-makers in planning and managing water resources, aiming to enhance efficiency and address water scarcity issues. The fault density map for the Naein Dasht basin was prepared and classified into five categories. The northwestern part of the region exhibits a higher fault density compared to other areas.
Land Cover:
The presence of vegetation cover can have significant positive effects on water resource management. This vegetative cover not only prevents soil erosion and surface runoff but also facilitates water infiltration and absorption, contributing to the improvement of groundwater quality and quantity. Therefore, preserving and increasing vegetation cover is of great importance in water resource conservation. The Naein Dasht region map was classified into six categories: barren areas, saline lands, agricultural lands, poor pastures, foothills, and rocky areas.
Precipitation:
Groundwater sources are often recharged through subsurface rainfall, where water infiltrates into underground layers and nourishes them. Factors such as topography (land slope), lithology (rock characteristics), vegetation cover, and soil properties play crucial roles in controlling this process. These factors can regulate how water is absorbed and permeates into the soil and subsurface layers, thereby exerting different effects on groundwater resources. Information obtained from studying these factors can be highly significant for better water resource management and determining the potential of areas for groundwater storage and supply.
Therefore, the precipitation layer of the region is one of the influential layers in groundwater potential assessment. This map was prepared using Kriging tools in ArcMap, based on the total precipitation over the past 20 years. The precipitation map of the Naein Dasht region was divided into 5 classes based on the amount of precipitation, and higher scores were assigned to areas with higher precipitation.
Distance map from ground water level:
The Distance map from ground water level map was obtained from the groundwater level data of the South Khorasan province. Then, using the Euclidean Distance tool in ArcMap 10.8, a map of the distance from the groundwater level was created and classified into 5 categories. The scoring criterion was based on the idea that the farther the distance from the groundwater level, the lower the score assigned.
In continuation, maps have been classified based on these values, as shown in the Table 3 and Fig. 2.
Table 3
Assigned Values to Each Layer
Layer
|
Distance(m)
|
Value
|
Distance from the waterway
|
1000-0
|
5
|
2000 − 1000
|
4
|
3000 − 2000
|
3
|
4000 − 3000
|
2
|
4000<
|
1
|
Layer
|
range
|
Value
|
Slope
|
5 − 0
|
5
|
10 − 5
|
4
|
15 − 10
|
3
|
20 − 15
|
2
|
20<
|
1
|
Geomorphology
|
class
|
Value
|
plain
|
3
|
Pastures with moderate slope
|
2
|
Lands with steep slopes
|
1
|
Soil
|
soil pattern
|
Value
|
very light
|
4
|
style
|
3
|
Medium (loamy silty).
|
2
|
heavy (silty clay)
|
1
|
|
class (km/km2)
|
Value
|
Lineament density(fault)
|
0.2-0
|
1
|
0.4 − 0.2
|
2
|
0.4–0.6
|
3
|
0.8 − 0.6
|
4
|
1.8–0.8
|
5
|
|
class
|
Value
|
land use
|
Salt marshes and barren areas
|
5
|
agriculture
|
4
|
Poor pastures
|
3
|
mid-range
|
2
|
Rock
|
1
|
|
Precipitation
|
Value
|
rainfall
|
50
|
1
|
75
|
2
|
100
|
3
|
150
|
4
|
200
|
5
|
Weighting of Information Layers:
As mentioned, in this study, the Fuzzy Analytical Hierarchy Process (Fuzzy AHP) method has been utilized. Fuzzy AHP was introduced by two Dutch researchers, Van Laarhoven and Pedrycz (1983). They presented this method as an improved approach to address the problems existing in the Analytical Hierarchy Process (AHP), including scale inequalities, uncertainty, difficulties in precise pairwise comparisons, and incomplete understanding of human thought elements. This method, using fuzzy numbers in a triangular form in the pairwise comparison matrices and employing the least squares logarithmic method, has been developed (Habibi et al., 2014; Anabestani and Vasal, 2016).
Fuzzy-AHP is a mathematical theory designed to model uncertainty in decision-making processes related to human decision-makers. This method allows decision-makers to freely choose their preferred values. The uncertainty of decision-makers is expressed using fuzzy numbers, so Fuzzy-AHP uses a range of values to represent this uncertainty. This method uses fuzzy numbers for comparing options and employs the geometric mean method for weighting and prioritization. It can easily be extended to fuzzy cases and helps determine a unique solution for the pairwise comparison matrix.
In this study, first, a questionnaire named 'Pairwise Comparison of Criteria (information layers discussed in the previous section)' is designed and prepared. For each of these criteria, the opinions of 20 experts and specialists from the province are collected, and a ranking is assigned for each criterion. Then, using the Fuzzy-AHP technique (Fuzzy Analytical Hierarchy Process), weighting is performed for these information layers, followed by normalization. In the final stage, a linear combination of weighting in a GIS environment is used to integrate thematic maps to produce a groundwater potential map (Machiwal, Jha, & Mal, 2011). In the Table 4, the phase spectrum equivalent to the scale of nine hours in the AHP technique is shown.
Table 4
Fuzzy spectrum equivalent to nine-hour scale in AHP technique (Habibi et al., 2013)
Fuzzy Inverse
|
Fuzzy numbers
|
Comparative verbal phrase i compared to j
|
(1,1,1)
|
(1,1,1)
|
Preferred Equally
|
(0.333,0.5,1)
|
(1,2,3)
|
|
(0.25,0.333,0.5)
|
(2,3,4)
|
Preferred moderately
|
(0.2,0.25,0.333)
|
(3,4,5)
|
|
(0.166,0.2,0.25)
|
(4,5,6)
|
Preferred Strongly
|
(0.142,0.16,0.2)
|
(5,6,7)
|
|
(0.125,0.142,0.166)
|
(6,7,8)
|
Very strongly Preferred
|
(0.111,0.125,0.142)
|
(7,8,9)
|
|
(0.111,0.111,0.111)
|
(9,9,9)
|
Extremely Preferred
|
In the fuzzy hierarchical method, to achieve the final map of groundwater potential, each criterion is weighted, and then they are combined based on these weights. To compare the criteria and weight them, the opinions of experts and specialists in this field are considered as criteria.
The binary comparison matrix of the criteria is displayed in EXCEL software, and the pairwise criteria comparison chart, which is visible in the Fig. 3 and Fig. 4. This process is carried out to determine the relative importance of the criteria in forming the groundwater potential map.
By conducting the criteria comparison, it was revealed that the precipitation criterion with a normalized weight of 0.313 has the most significant impact, while the geomorphology criterion with a weight of 0.051 has the least impact. In other words, according to the fuzzy AHP analysis, the precipitation criterion is more important than the other criteria in determining the potential of groundwater resources in this study.
Calculation of Inconsistency Rate:
The inconsistency rate is an index that indicates the degree of compatibility between pairwise comparisons. Pairwise comparisons are made by experts or specialists, and if these comparisons are compatible with each other, the results can be trusted. In fact, the inconsistency rate reflects the extent of disagreement or lack of consensus among the members of the study group regarding pairwise comparisons. If the inconsistency rate is very low, it indicates greater agreement among the members on pairwise comparisons, and there is more confidence in the results. However, if the inconsistency rate is high, it may indicate a lack of agreement and significant differences of opinion within the group, which may lead to uncertainty in the comparisons.
Weight assignments will be compatible only when the consistency ratio is equal to or less than 10%; otherwise, weights need to be re-evaluated and estimated to reduce inconsistency (Saaty, 1990).
The assigned weights will be compatible only when the consistency ratio (CR) is equal to or less than 10%. Otherwise, the weights should be re-evaluated and estimated to reduce inconsistency (Saaty, 1990).
The inconsistency rate of each matrix can be determined through the following steps:
1. Form the pairwise comparison matrix A.
2. Determine the weight vector (W).
3. Calculate the maximum eigenvalue of A (\(\:{{\lambda\:}}_{max}\)).
Estimate the values of \(\:{{\lambda\:}}_{max}\)by dividing the values obtained for \(\:{{\lambda\:}}_{max}\).W by the corresponding W.
In this study, the value of \(\:{{\lambda\:}}_{max}\)was found to be 8.95%.
Table 5
Consistency Index (Qudsiehpour, 2008)
n
|
1
|
2
|
3
|
4
|
5
|
6
|
7
|
8
|
9
|
10
|
R
|
0
|
0
|
0.58
|
0.9
|
1.12
|
1.24
|
1.32
|
1.41
|
1.45
|
1.45
|
Combining Information Layers and Creating Groundwater Potential Map:
To generate a groundwater potential map, you need to first obtain maps corresponding to each sub-criterion based on their normalized weights.
Table 6
Weights of Criteria and Sub-criteria, and Final Weights
Fuzzy-AHP
|
sub- criteria
|
Criterion
|
Final weight
|
sub-criteria weight
|
Criterion weight
|
0.026
|
0.098
|
0.313
|
50
|
Rainfall
|
0.040
|
0.152
|
75
|
0.052
|
0.195
|
100
|
0.064
|
0.240
|
150
|
0.084
|
0.315
|
200
|
0.063
|
0.319
|
0.233
|
1000
|
Distance from the waterway
|
0.050
|
0.253
|
2000
|
0.034
|
0.172
|
3000
|
0.027
|
0.137
|
4000
|
0.023
|
0.120
|
> 4000
|
0.032
|
0.360
|
0.106
|
Very light texture
|
Soil
|
0.024
|
0.270
|
Light texture
|
0.019
|
0.210
|
Medium texture
|
0.014
|
0.160
|
Heavy texture
|
0.014
|
0.318
|
0.054
|
0–5
|
Slope
|
0.010
|
0.228
|
5–10
|
0.009
|
0.198
|
10–15
|
0.006
|
0.137
|
20 − 15
|
0.005
|
0.120
|
> 20
|
0.007
|
0.099
|
0.084
|
0.2-0
|
Fault
|
0.010
|
0.140
|
0.4 − 0.2
|
0.012
|
0.172
|
0.4–0.6
|
0.020
|
0.278
|
0.8 − 0.6
|
0.022
|
0.311
|
1.8–0.8
|
0.048
|
0.359
|
0.158
|
Barren and salt marsh areas
|
Land use
|
0.027
|
0.207
|
agriculture
|
0.026
|
0.199
|
Poor pastures
|
0.016
|
0.125
|
mid-range
|
0.014
|
0.109
|
Rock
|
0.019
|
0.443
|
0.051
|
plain
|
Geomorphology
|
0.013
|
0.317
|
Medium slope pastures
|
0.010
|
0.240
|
Steep land
|
Table 6 shows the weights of criteria, sub-criteria, and the final weight for each. In this table, the final weight for each sub-criterion is obtained by multiplying the weight of the main criterion by the weight of the sub-criterion. To generate the groundwater potential map, the maps corresponding to each sub-criterion are combined, and then, using the fuzzy hierarchical method, the map is classified into 5 quantitative categories, ranging from very good to very poor. As a result of this classification, approximately 16% of the area is classified as very poor, 22% as poor, 24% as moderate, 22% as good, and 14% as very good.
Validation
Through the analysis of the locations of wellsqanats, and springs in the area, the results indicate that the region is very good in terms of groundwater potential and possesses significant groundwater resources. These findings demonstrate that the use of remote sensing techniques, Geographic Information System (GIS), and fuzzy-AHP modeling are effective in estimating the potential of groundwater resources and decision-making in water management. According to the information in Fig. 9, there are 10 well rings, 3 qanats, and 3 springs in the region, with a significant portion of the resources derived from groundwater placed in the very good zone in terms of groundwater potential. This reflects the high accuracy of the research results in employing GIS, and fuzzy-AHP techniques in groundwater potential assessment.