[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Optimal Inversion of Conversion Parameters from Satellite AOD to Ground Aerosol Extinction Coefficient Using Automatic Differentiation
Next Article in Special Issue
Land Use/Land Cover Changes Impact on Groundwater Level and Quality in the Northern Part of the United Arab Emirates
Previous Article in Journal
Monitoring of Coral Reefs Using Artificial Intelligence: A Feasible and Cost-Effective Approach
Previous Article in Special Issue
Application of Probabilistic and Machine Learning Models for Groundwater Potentiality Mapping in Damghan Sedimentary Plain, Iran
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Ensemble of MCDM-Artificial Intelligence Techniques for Groundwater-Potential Mapping in Arid and Semi-Arid Regions (Iran)

by
Alireza Arabameri
1,
Saro Lee
2,3,*,
John P. Tiefenbacher
4 and
Phuong Thao Thi Ngo
5,*
1
Department of Geomorphology, Tarbiat Modares University, Tehran 36581-17994, Iran
2
Geoscience Platform Research Division, Korea Institute of Geoscience and Mineral Resources (KIGAM), 124, Gwahak-ro Yuseong-gu, Daejeon 34132, Korea
3
Korea University of Science and Technology, 217 Gajeong-ro Yuseong-gu, Daejeon 34113, Korea
4
Department of Geography, Texas State University, San Marcos, TX 78666, USA
5
Institute of Research and Development, Duy Tan University, Da Nang 550000, Vietnam
*
Authors to whom correspondence should be addressed.
Remote Sens. 2020, 12(3), 490; https://doi.org/10.3390/rs12030490
Submission received: 25 December 2019 / Revised: 30 January 2020 / Accepted: 2 February 2020 / Published: 4 February 2020
Graphical abstract
">
Figure 1
<p>Location of study area in Semnan Province and Iran (inset) and the locations of training and validation wells in the study area.</p> ">
Figure 2
<p>Flowchart of the research.</p> ">
Figure 3
<p>Location of random points in the study area.</p> ">
Figure 4
<p>Groundwater conditioning factors. (<b>a</b>) elevation, (<b>b</b>) slope, (<b>c</b>) aspect, (<b>d</b>) stream power index (SPI), (<b>e</b>) topography position index (TPI), (<b>f</b>) topography wetness index (TWI), (<b>g</b>) terrain ruggedness index (TRI), (<b>h</b>) convergence index (CI), (<b>i</b>) distance to stream, (<b>j</b>) distance to road, (<b>k</b>) distance to fault, (<b>l</b>) drainage density, (<b>m</b>) rainfall, (<b>n</b>) soil type, (<b>o</b>) land use/land cover (LULC), (<b>p</b>) lithology, (<b>q</b>) normalized difference vegetation index (NDVI).</p> ">
Figure 4 Cont.
<p>Groundwater conditioning factors. (<b>a</b>) elevation, (<b>b</b>) slope, (<b>c</b>) aspect, (<b>d</b>) stream power index (SPI), (<b>e</b>) topography position index (TPI), (<b>f</b>) topography wetness index (TWI), (<b>g</b>) terrain ruggedness index (TRI), (<b>h</b>) convergence index (CI), (<b>i</b>) distance to stream, (<b>j</b>) distance to road, (<b>k</b>) distance to fault, (<b>l</b>) drainage density, (<b>m</b>) rainfall, (<b>n</b>) soil type, (<b>o</b>) land use/land cover (LULC), (<b>p</b>) lithology, (<b>q</b>) normalized difference vegetation index (NDVI).</p> ">
Figure 5
<p>Relative importance of conditioning factors using random forest model.</p> ">
Figure 6
<p>Weighted GWCFs using the FR model for (<b>a</b>) elevation, (<b>b</b>) slope, (<b>c</b>) aspect, (<b>d</b>) stream power index (SPI), (<b>e</b>) topography position index (TPI), (<b>f</b>) topography wetness index (TWI), (<b>g</b>) terrain ruggedness index (TRI), (<b>h</b>) convergence index (CI), (<b>i</b>) distance to stream, (<b>j</b>) distance to road, (<b>k</b>) distance to fault, (<b>l</b>) drainage density, (<b>m</b>) rainfall, (<b>n</b>) soil type, (<b>o</b>) LU/LC, (<b>p</b>) lithology, (<b>q</b>) NDVI.</p> ">
Figure 6 Cont.
<p>Weighted GWCFs using the FR model for (<b>a</b>) elevation, (<b>b</b>) slope, (<b>c</b>) aspect, (<b>d</b>) stream power index (SPI), (<b>e</b>) topography position index (TPI), (<b>f</b>) topography wetness index (TWI), (<b>g</b>) terrain ruggedness index (TRI), (<b>h</b>) convergence index (CI), (<b>i</b>) distance to stream, (<b>j</b>) distance to road, (<b>k</b>) distance to fault, (<b>l</b>) drainage density, (<b>m</b>) rainfall, (<b>n</b>) soil type, (<b>o</b>) LU/LC, (<b>p</b>) lithology, (<b>q</b>) NDVI.</p> ">
Figure 7
<p>Groundwater potential map using the VIKOR-RF ensemble model.</p> ">
Figure 8
<p>Percentage of potential classes.</p> ">
Figure 9
<p>Area under curve (AUC). (<b>a</b>) Prediction rate curve (PRC), (<b>b</b>) Success rate curve (SRC).</p> ">
Figure 10
<p>Trend of frequency ratio and seed cell area index indicators.</p> ">
Versions Notes

Abstract

:
The aim of this research is to introduce a novel ensemble approach using Vise Kriterijumska Optimizacija I Kompromisno Resenje (VIKOR), frequency ratio (FR), and random forest (RF) models for groundwater-potential mapping (GWPM) in Bastam watershed, Iran. This region suffers from freshwater shortages and the identification of new groundwater sites is a critical need. Remote sensing and geographic information system (GIS) were used to reduce time and financial costs of rapid assessment of groundwater resources. Seventeen physiographical, hydrological, and geological groundwater conditioning factors (GWCFs) were derived from a spatial geo-database. Groundwater data were gathered in field surveys and well-yield data were acquired from the Iranian Department of Water Resources Management for 89 locations with high yield potential values ≥ 11 m3 h−1. These data were mapped in a GIS. From these locations, 62 (70%) were randomly selected to be used for model training, and the remaining 27 (30%) were used for validation of the model. The relative weights of the GWCFs were determined with an RF model. For GWPM, 220 randomly selected points in the study area and their final weights were determined with the VIKOR model. A groundwater potential map was created by interpolating the values at these points using Kriging in GIS. Finally, the area under receiver operating characteristic (AUROC) curve was plotted for the groundwater potential map. The success rate curve (SRC) was computed for the training dataset, and the prediction rate curve (PRC) was calculated for the validation dataset. Results of RF analysis show that land use and land cover, lithology, and elevation are the most significant determinants of groundwater occurrence. The validation results show that the ensemble model had excellent prediction performance (PRC = 0.934) and goodness-of-fit (SRC = 0.925) and reasonably high classification accuracy. The results of this study could aid management of groundwater resources and assist planners and decision makers in groundwater-investment planning to achieve sustainability.

Graphical Abstract">

Graphical Abstract

1. Introduction

Increasing demands for water due to population growth and industrialization combined with anticipated changes in aquatic resources due to global warming and regional climate change highlight the urgent need for a quantitative methodology that can model groundwater production [1,2,3,4,5,6,7]. Groundwater-potential mapping (GPM) is the first and most important step in groundwater management [8]. Since the methods to identify the parameters that influence the spatial distribution of groundwater and the ways to acquire data are evolving, new approaches to develop accurate and useful information for decision makers are needed [9]. GPM is effectively achieved by combining field study, remote sensing (RS), and geographic information system (GIS) methods [10,11]. Over the past decade, the growing use of satellite data, thematic maps, and land cover and land use (LULC) data has made it easier to study groundwater potential. As the range of data and number of variables that must be gathered and compiled spatially to simultaneously examine the critical factors grows, GIS is increasingly useful because of its capacity to process diverse spatial data from diverse sources and in different formats. The combination of GIS with new modeling methods can be a powerful tool for decision makers [12,13].
Several statistical, multi-criteria decision-making (MCDM), and machine learning approaches have been developed for GPM with thematic data layers that describe soils, topography, hydrology, lithology, and LULC [1]. The most important advantages of these methods avoid the need for highly complex or expensive data and produce results that are simple to analyze and evaluate [14]. Data needed for these models are easily acquired by integrating RS data and field surveys and can be routinely integrated in a GIS [15,16]. The extraction of spatial information from RS data has been undertaken for many tasks, including prediction of the spatial occurrence of extreme processes that yield environmental hazards (such as landslides, soil erosion, and flooding) and prediction of the distribution and dimensions of natural resources (like groundwater) [17,18,19]. Using RS and GIS, it is possible to develop a GPM with limited costs and reasonable levels of accuracy over large areas [20]. GIS enables faster spatial analysis, combines a large array of data that describe diverse spatial characteristics from diverse sources, and makes information management easier [21,22,23]. GIS is used in environmental management because it offers several benefits. High-quality data analysis can be conducted with it and it is often adaptable to the use of new and advanced analytical methods. It enables easy and fast management and analyses of large datasets. Innovations occur constantly and there are improvements in methods to create and update maps. GIS can be used to link phenomena and relate vast amounts of information in data tables. Data can be rapidly updated and integrated from continuous and scattered data acquisition activities. GIS is methodologically easy to review, and such systems can be used to model, test hypotheses, and to make predictions [24].
Many powerful statistical, machine learning, and MCDM methods have been developed and tested by researchers. These can be classified into several groups:
  • MCDM-based models, like the analytic hierarchy process (AHP) [25,26],
  • Bivariate and multivariate statistical models, like weights of evidence (WoE) [27], evidential belief function (EBF) [28,29,30], Shannon’s entropy (SE) [31], maximum entropy (ME) [32], frequency ratio (FR) [33], certainty factor (CF) [34], and logistic regression (LR) [35,36],
  • Data mining/machine learning models, like boosted regression tree (BRT) [37], classification and regression tree (CART) [38], multivariate adaptive regression splines (MARS) [37], artificial neural network (ANN) [39], random forest (RF) [40], fuzzy logic (FL) [41], Fisher’s linear discriminant function (FLDA) [42], support vector machine (SVM) [43], logistic model tree (LMT) [44], K-nearest neighbor (KNN) [45], and quadratic discriminant analysis [46].
Many approaches have been used to identify areas with high groundwater potential and the array of approaches has increased in recent years because of the growing pressure to clarify the distribution of groundwater and to identify new sources of water. High-precision maps are needed to better manage aquifers that are in use and to avoid costly exploration for new deposits, especially during droughts. Hybrid groundwater mapping techniques are a vital necessity [35]. Ensemble machine learning frameworks can be used to develop better prediction models and they easily cope with complex multi-dimensional data [40]. Kordestani et al. [30] hybridized an EBF and BRT (i.e., creating an EBF-BRT) algorithm and demonstrated a groundwater contamination-susceptibility model. Naghibi et al. [47] hybridized four data mining models and FR into a novel ensemble model for a groundwater-related study. Their ensemble model exhibited better performance than stand-alone applications of each separate model. Chen et al. [35] developed a novel hybrid approach with WoE, LR, and functional tree (FT) models that individually outperformed the models. Miraki et al. [40] used a random subspace (RS) ensemble-RF (i.e., RS-RF) classifier to model groundwater susceptibility and found that the hybrid model outperformed individual models. Razavi-Termeh et al. [48] integrated three bivariate statistical models with RF and logistic tree models to map groundwater potential and revealed that the RF statistical models improved prediction accuracy. Though the idiographic studies demonstrated hybrid models that performed strong groundwater susceptibility modeling in given study areas, there is no evidence to show that they are universally appropriate or applicable. The applicability of the novel approaches for better prediction of groundwater susceptibility models needs further evaluation.
Therefore, the main purpose of this research is to evaluate the integration of MCDM with machine-learning approaches to map groundwater potential. FR models, based on simple statistical methods, are used to examine correlations between groundwater-related factors and groundwater potential. This ensemble approach is novel and has not yet been applied in groundwater potential mapping. This study evaluated the performance of the Vise Kriterijumska Optimizacija I Kompromisno Resenje (VIKOR)-RF-FR ensemble approach for groundwater potential mapping. Statistical evaluation metrics—area under the receiver operating characteristic curve (AUC-ROC), FR, and seed cell area index (SCAI)—were used to determine the goodness-of-fit and predictive performance of the ensemble algorithm. Data from the advanced land-observing satellite (ALOS) phased array type L-band synthetic aperture radar (PALSAR) digital elevation model (DEM) with a spatial resolution of 12.5 m were used to acquire topographical and hydrological factors. ALOS-PALSAR DEM has been determined to have a higher accuracy than advanced spaceborne thermal emission and reflection radiometer (ASTER) and shuttle radar topography mission (SRTM) DEMs [49].
MCDM models are useful for solving complex problems because of their capacity to enable examination and selection of alternatives based on criteria [50,51,52,53,54,55,56,57,58,59]. MCDMs can analyze basic statistical trends in big data in various study areas [60,61]. These models are an important part of modern decision-making science which is decision-making in the face of varying decision criteria and selection from alternative outcomes. MCDM methods have been developed not only because of the diversity of processes and problems that require assessment of multiple factors, but because experts desire to improve decision making with better techniques [62,63,64,65]. The improvements are the outcomes of advances in mathematical science, computational science, and computer technology. Machine-learning models have also gained popularity among environmental managers and decision makers [66,67,68]. Machine-learning models compare favorably to other models because they do not require a defined and strict set of a priori assumptions as do many other statistical methods. These models are also used to assess the relative importance of each of the conditioning factors [66]. The FR model simplifies the probabilistic relationships between the dependent and independent variables to quantify the relationships between the reclassified factors and groundwater capacity. The relationships calculated using the FR model enabled easier interpretation of groundwater potential indicated by the rating of each factor.

2. Materials and Methods

2.1. Study Area

Bastam watershed, with an area of 1329.43 km2, is located between 36°27′02″ and 36°47′13″N and between 54°24′23″ and 55°11′08″E in Semnan Province, Iran (Figure 1). The region is topographically mountainous with an elevation range from 1357 to 3893 m.a.s.l. and an average elevation of 2078.3 m asl. The average slope in the study area is 13.5° and the maximum is 70.6°. The region’s climates are arid and semi-arid: the average annual rainfall of the study area is 262.4 mm and the average annual temperature is 12.8 °C [69]. Valley terrace deposits, stream channels, braided channels, flood plain deposits, low elevation piedmont fans, and valley terrace deposits are the most important lithological units of the study area [70]. Entisols (32.9%), entisols/inceptisols (32.4%), and mollisols (21.5%) are the most widespread soils in the study area [71].
Geologically, Devonian-Eocene rocks, including andesitic and basaltic volcanic, conglomerate, and sandstone, are distributed as mountain ranges in the northern part of the study area extending from the eastern to the western boundaries. Cretaceous granitic and volcanic rocks are intruded as mountain ranges and isolated hills. They are found in the eastern part of the study area. Moreover, Quaternary sediments, fluvial deposits along streams, and collovium, are found in the central section of study area. Groundwater in this area is stored in these Quaternary sediments, both alluvial and sedimentary formations, in cracks, fractures, faults, bedding planes, and large cavities of consolidated rocks of Permo-Carboniferous meta-sediment and granitic aquifers.
Due to Iran’s aridity, industrialization, and population growth, groundwater is a vital natural resource. More than 90% of groundwater is used for agriculture [72]. During the period 1971–2014, the number of water wells in the country increased dramatically, from 47,000 to more than 789,000. There are currently about 196,000 deep wells and 593,000 semi-deep wells in Iran and they extract an average of 0.17 and 0.02 million m3, respectively [72]. Groundwater consumption in Iran has gone from less than 18 billion m3 in 1971 to more than 61 billion m3 in 2014. In recent years, average groundwater levels in the country have dropped by about 51 cm, a decrease of 5.086 billion m3/year [48]. Approximately 20,000 people depend on the groundwater resources of the Bastam watershed for personal consumption and agriculture [73]. Mismanagement of water in this region has led to a 15 m drop of water tables over a 20-year period and to land-surface subsidence [73].

2.2. Methodology

This research involved six steps (Figure 2):
(i) Prepare the database—A groundwater-inventory map (GWIM) was created and 17 groundwater conditioning factors (GWCFs) were added as layers to a GIS (ArcGIS). The groundwater inventory locations were separated into two randomly selected groups for training and validation.
(ii) Assess the data—A multi-collinearity analysis was performed on the GWCFs using indices of tolerance (TOL) and the variance inflation factor (VIF). TOL < 0.2 and VIF > 5 indicate a multicollinearity issue among the variables [74]. The linearity of the GWCFs was assessed, and increasing linearity between the factors will reduce the accuracy of the models [65].
(iii) Calculate relative importance of GWCFs—Because a VIKOR MCDM model cannot determine the relative importance of predictive variables, a RF model was used. To do this, the same number of non-well locations were randomly selected in the study area as there were wells yielding more than 11 m3/h and they were used to train the RF model. Next, the Extract Multi Values to Point tool in ArcGIS extracted the values of each GWCF at well (coded 1) and non-well (coded 0) locations. The data were formatted in Excel and imported into R 3.1.1 software. The Rattle package in R software was used to calculate the weight of each GWCF. Given a groundwater dataset, GWD = (xi, yi), (i = 1, 2, 3, … N) with xi ∈ RM are the conditioning factors. N is the total number of samples, and RM is the total number of conditioning factors, yi ∈ (1,0) is the output which contains two classes, well and non- well locations, with the first coded “1” and the second coded “0”. The bootstrap aggregating algorithm generated n subsets, which consist of m factors (m ≤ M). A CART algorithm constructed a tree classifier for each subset. The tree classifiers were aggregated to an RF classifier. Two parameters of each the RF, n, and m, were determined when building the RF model. The parameter, n, should be large enough to ensure the robustness of the RF model (i.e., n = 500).
(iv) Prepare the decision matrix—Two hundred and twenty-two random points (Figure 3) (called alternatives) were selected in the study area and their GWCF values were extracted. A decision matrix of 17 columns (the GWCFs) and 220 rows (alternatives) was created. FR models were used to calculate the values of the alternatives’ pixels. In MCDM procedures, GWCF values should have either ascending or descending orders. The actual values did not follow this trend in this study, so FR was used to fit an ascending trend. Higher GWCF values indicated locations that were more likely to have groundwater.
(v) Run the hybrid model and map groundwater—The VIKOR was run to develop the groundwater potential model. The produced values indicated strength of association with the presence of groundwater at the 220 locations. Using the Ordinary Kriging tool in ArcGIS 10.5, the distribution of groundwater potential was mapped. Ordinary Kriging has been shown to efficiently map nonlinear relationships and has been used for a variety of real-world applications [75].
(vi) Assess performance—The receiver operating characteristics (ROC) curves and their area under the curve (AUROC), frequency ratio (FR), and seed cell area index (SCAI), were used to evaluate the accuracy of the model predictions.

2.3. Data Preparation

The most important step in groundwater-potential mapping (GWPM) is the creation of a well-inventory map (WIM) [64]. Extensive field surveys were conducted, and well-yield data were acquired from the Iranian Department of Water Resources Management to characterize the known spatial distribution of groundwater in the region. Actual pumping-test production measurements of wells expressed in m3h−1 were used to represent yield. Eighty-nine wells with a high yield (≥11 m3h−1) were selected for GWPM [1] and modeling (Figure 1). The data from these wells were randomly grouped into a training dataset of 62 (70%) wells and a validation dataset of 27 (30%) to be used to test the model [76].
Geo-environmental features were used as GWCFs to predict groundwater potential [1]. Based on a literature review [1,30,32,36,37,39,40], availability of data, environmental characteristics, scale of research, and a multicollinearity test, 17 RS-derived GWCFs that reflected topography, hydrology, and LULC were selected. Some GWCF measurements were derived from the ALOS PALSAR DEM with a spatial resolution of 12.5 m [1,52,53]. RS and thematic maps were also used to measure others. The GWCFs are elevation, slope, aspect, stream power index (SPI), topographic position index (TPI), topographic wetness index (TWI), terrain ruggedness index (TRI), convergence index (CI), distance to streams, distance to faults, distance to roads, drainage density (DD), rainfall, soil type, land use and land cover (LULC), lithology, and normalized difference vegetation index (NDVI) (Figure 4a–q).
Elevation, slope, and aspect were derived from ALOS PALSAR DEM in ArcGIS. Distances to streams, roads, and faults were computed as Euclidean distances pixel centroids to the nearest polyline. The vector-based information was then rasterized to coincide with the resolution of the DEM. Road and fault-factor maps (scales of 1:50,000 and 1:100,000, respectively) were obtained from the National Geographic Organization of Iran (www.ngo-org.ir) and the Geological Society of Iran (GSI) (http://www.gsi.ir/), respectively. Descriptions of lithological units in the study area were compiled. A factor map of streams was extracted from the DEM in ArcHydrov10.4 and ArcGISv10.4. Soil texture, lithology, and LULC were derived from three thematic maps: the soil map of the Semnan Agricultural and Natural Resources, Research and Education Center, the geological map of the Geological Society of Iran (GSI) (http://www.gsi.ir/), and the land use map created by the Iranian Soil Conservation and Watershed Management Research Institute (https://www.scwmri.ac.ir). All three of these maps were at the 1:100,000 scale.
Annual precipitation data from nine weather stations (at Mojen, Farahzad, Bastam, Abr, Karkhaneh, Semnan, Shahroud, Tarzeh, and Shah Kuh-e Bala) for the period from 1986 to 2016 were mapped and Kriging [1,24,51,54] was used to map annual rainfall patterns in ArcGIS 10.5 [69]. Land uses significantly influence the distribution of groundwater [1]. Greater vegetative cover tends to both indicate and promote the presence and extent of aquifers, as forests and denser vegetation reduce runoff rates and encourage infiltration more than do barren and sparsely vegetated lands. Equations (1)–(4) were used to calculate TWI, SPI, TRI, and TPI [77,78,79].
TWI = In   ( A S / tan β )
SPI = A S   × tan β
TRI = | X | ( max 2 min 2 )
TPI = E pixel E surrounding
where A S is the specific catchment area of the basin (m2/m), β is slope steepness in degrees, x indicate elevation of each neighbor cell to cell (0, 0) (m), and max and min show the largest and smallest elevation value among nine neighbor pixels, respectively. Epixel depicts the elevation of the cell, and Esurrounding is the mean elevation of the neighboring pixels.
TPI compares the elevation of each pixel in the DEM with the average height of the pixels around it. This factor enables classification of landscapes into morphological classes. The positive and negative values indicate that a pixel is higher or lower than the pixels that surround it [80]. The plan curvature reflects the directional variations along a curve. The effect of plan curvature on slope erosion is that it causes convergence and divergence of water along and away from a flowing stream [81]. The profile curvature represents the amount of elevation variation along the flow path. When the concavity of the surface curvature is increasing, its values are negative, and, if decreasing, they are positive [82]. This index indicates the flow velocity, erosion (in negative amounts), sedimentation (in positive values), and geomorphology of the area [83]. TPI, plan curvature, profile curvature, and convergence index were extracted from the DEM in system for automated geoscientific analyses (SAGA-GIS).
NDVI is used to evaluate vegetation types and density using quantitative analysis of albedo or reflectance of visible and NIR portions of the electromagnetic spectrum. On a scale from −1 to +1, high vegetation density usually has values between 0.05 and 1, bare soil areas are close to 0, and water bodies have values closer to −1. Vegetation depends upon moisture and promotes infiltration into soils and therefore is highly indicative of a likely presence of groundwater [1]. NDVI values were computed from a LANDSAT-8 image with 30 m spatial resolution acquired from the USGS (https://earthexplorer.usgs.gov/) (Table 1).
The GWCFs used in the GWPM were converted to raster layers at a 12.5 m resolution. The raster matrix was 2395 columns by 1460 rows.
Classes of LULC: A) Agriculture, B) dense forest, C) good range, D) agri-dry farming, E) dry farming, F) low forest, J) Woodland, H) mod-forest, I) mod-range, J) poor-range, K) rock, and L) Urban.
Classes of lithology: A) limestone, B) marl and shale, C) tuffaceous shale, Andesitic and basaltic volcanics, D) limestone, and shale, E) marl, shale and detritic limestone, F) Red conglomerate and sandstone, G) siltstone and shale, P) fluvial conglomerate, and sandstone, I) braided channel and flood plain deposits, J) dolomite, shale and sandstone.

2.4. Models

2.4.1. Frequency Ratio (FR)

FR is the probability of occurrence of a specific phenomenon. FR determines the relationship between GWCFs and well locations [61]. To calculate FR for a conditioning factor, Equation (5) is used [84]:
FR = A / B C / D
where A is the number of wells in a class, B is the total number of wells, C is the number of pixels in that class, and D is the number of pixels [85].

2.4.2. VIKOR (Vlse Kriterijumsk Optimizacija Kompromisno Resenje)

The VIKOR method was first introduced by Opricovic [86]. VIKOR is derived from compromising programming. This model is an MCDM method that determines the best solution for a discrete decision that has disproportionate and conflicting criteria. This method focuses on ranking and selection from alternatives and it determines a compromise solution to guide decision-making. The compromise solution is a feasible solution that is the closest solution to the ideal solution. Compatibility is also based on concessions [87]. The following six steps are involved in the VIKOR model:
(1)
Prepare the decision matrix.
(2)
Calculate the normalized matrix, as shown in Equation (6):
r ij = x ij 1 m x ij
where, xij is the ith alternative performance of the jth criterion, and m is the alternative numbers.
(3)
Calculate a weighted normalized matrix, as shown in Equation (7):
f ij = r ij × W j
where Wj is the weight of the criterion.
(4)
Identify the ideal positive (Equation (8)) and negative (Equation (9)) options:
f i + = max i f ij
f i = min i f ij
where   f i + is the positive ideal solution for the ith criterion and f i is the negative ideal solution for the ith criterion.
(5)
Calculate the utility index (Equation (10)) and incompatibility index (the distance from positive and negative ideal solution) (Equation (11)):
S i = L 1 , i = i = 1 n W j ( f i + f ij ) ( f i + f i )
R i = L , i = max [ i = 1 n W j ( f i + f ij ) / ( f i + f i ) ]
where   S i and R i are the distance of the ith alternative to the ideal positive and negative solutions, respectively.
(6)
Calculate the VIKOR index and determine the final weight of the alternatives, as shown in Equation (12):
Q i = V × ( S i S ) ( S S ) +   ( 1 V ) ×   ( R i R ) ( R R )
where V is a constant (0.5), S is min Si, S is max Si, R is min Ri, and R is max Ri. The greater the value of the VIKOR index in one alternative, the greater the importance of that alternative.

2.4.3. Random Forest

RF is one of the most well-known machine-learning algorithms [88]. RF is nonparametric and based on decision trees. Many decision trees are grown in the classification of RF algorithms [89]. The RF algorithm:
Considers the number of tests (N) and the number of variables (M),
Enters R(m) variables to decide on each tree node (m should be less than M),
Selects test data for the tree by using the n-times placement of N samples, and the rest of the samples are used to estimate the tree-error,
Selects M variables for each tree node, the basis for decision making in each node. The best groups are calculated on m variables in the test-sample, and
Expands each tree completely without pruning.
The GWCFs are prioritized using the mean decrease accuracy and the mean decrease Gini. Using the mean decrease accuracy, as opposed to the mean decrease Gini, determines the priority of effective factors more effectively, especially when environmental factors are being compared [68].

2.5. Validation of Results

AUC-ROC, FR, and SCAI criteria help to discern the quality of definite and probable identification, as well as the quality of prediction [65]. AUC-ROC represents a model’s predictive accuracy by describing its ability to predict whether predefined events have occurred or will occur [90]. AUC-ROC curves indicate the sensitivity of the model to the percentage of cells or unstable units correctly predicted by the model against the percentage of predicted unstable cells within the total. These values reflect the model’s ability to correctly distinguish between positive and negative observations in a validation sample. High sensitivity indicates that most predictions were correct (true positives), while specificity values of 1 indicate a high number of false positives. In the AUC-ROC, the false positive rate (1 specificity) is shown on the x axis (Equation (13)) and the true positive rate (sensitivity) on the y axis (Equation (14)):
x = 1 specificity = 1 [ TN ( TN + FP ) ]
y = sensitivity = [ TN ( TP + FN ) ]
where TN is true-negative, FP is false-positive, TP is true-positive, and FN is false-negative [85]. The quantitative–qualitative relationship classes between AUC-ROC and prediction accuracy, ranging from 0 to 1, are: excellent (0.9–1), very good (0.8–0.9), good (0.7–0.8), moderate (0.6–0.7), and weak (0.5–0.6) [80]. FR and SCAI criteria were used to evaluate the separations between classes. These criteria determine the classification accuracy of a model. The FR is the ratio of the percentage of wells in each class to the percentage area of that class [91], and SCAI is the ratio of the percentage area of each zoning class to the percentage of wells in each class [92].

3. Results

3.1. Multi-Collinearity Analysis

The multicollinearity test shows that there is no collinearity among conditioning factors (Table 2). TOL among factors ranges from 0.15 to 0.84 and VIF ranges from 1.1 to 4.8. Therefore, all GWCFs were included in the GWPM process.

3.2. Determining the Relative Importance of GWCFs using the RF Model

The RF model with an out-of-bag error rate = 8.57% was applied in R using the Caret package. This scenario implies that the accuracy of the model is equal to 91.43%, and exhibits the results of the confusion matrix. According to the results (Table 3), of the 80 non-well locations, 74 (92.5%) were predicted to be non-well and 6 (7.5%) were predicted to be locations of wells. Of the 80 well locations, 11 (13.75%) were predicted to be non-well locations and 69 (86.25%) were predicted to be wells. The determination of the relative importance of the GWCFs (using mean decrease accuracy of the RF model) (Figure 5) shows that LULC (18.5), lithology (13.9), and elevation (11.9) were the main factors that identified the locations of greater amounts of groundwater. These were followed by TPI (9.1), NDVI (7.9), TWI (7.7), distance to road (7.4), drainage density (6.3), slope aspect (5.7), soil type (5.3), distance to fault (5.03), TRI (4.9), slope degree (4.8), rainfall (4.7), CI (4.6), distance to stream (2.4), and SPI (1).

3.3. Determining the Value of Each Pixel of GWCFs using the FR Model

The FR model revealed the spatial relationships between the GWCFs and wells (Table 4). Analysis of the values of each GWCF at each well location (Figure 6a–q) shows that locations with lower elevations, slope angles, CI, SPI, distance to road, distance to stream, rainfall, TRI, and NDVI had lower FR values (i.e., these factors have strong inverse relationships with groundwater potential). Locations with higher TWI and drainage density have higher values of FR (i.e., strong positive relationships to groundwater potential). In terms of LULC and lithology, agriculture and quaternary formations, not surprisingly, were strongly associated with well locations.

3.4. Application of VIKOR Model

Extracting the GWCF values at random locations, the constructed decision matrix, and the final weight of random locations using the VIKOR model indicate that groundwater potential is strongly indicated by weighting values. The weights of random locations range from 0 to 1. Higher values indicate greater groundwater potential.

3.5. Groundwater Potential Mapping Using the FR-RF-VIKOR Ensemble Model

Groundwater potential was mapped with the ensemble model (values varied from 0.01 to 0.92). Based on the natural break classification method, these values were divided into five classes: very low (0.01–0.14), low (0.14–0.28), moderate (0.28–0.44), high (0.44–0.65), and very high (0.65–0.92) groundwater potential (Figure 7). About 48.3% of the total area of 642.1 km2 is classified (Figure 8) as having ‘very low’ groundwater potential, 22.8% (303.3 km2) has a ‘low’ groundwater potential zone, 11.8% of the area (157.8 km2) has ‘moderate’ potential, 8.5% (113.9 km2) has high groundwater potential, and 8.4% (112.05 km2) has very high groundwater potential.

3.6. Validation of Results

Validation of the ensemble model is based on the prediction rate curve (PRC) and the success rate curve (SRC) (Figure 9). The ensemble model, with PRC = 0.93 and SRC = 0.92, has excellent groundwater potential prediction accuracy. The FR and SCAI model results (Figure 10) show ascending and descending trends, respectively. In terms of statistical assessment of the modeling values, FR increased and SCAI decreased, indicating that the classification accuracy of the ensemble model was reasonable.

4. Discussion

Groundwater resources are crucial in regions with limited surface water, particularly in arid and semi-arid climates [81]. When these regions have significant amounts of industrial and agricultural water demands, aquifers can be consumed at rates that exceed recharge rates and are at increasing risk of irreparable damage. Therefore, it is important to have accurate assessments of the spatial dimensions of groundwater in a watershed. Such assessments can improve water-management planning and use strategies [93]. Groundwater mapping can be improved and made more economical by identifying the geophysical and hydrological conditioning factors that are associated with sub-surface storage. Several methods have been developed to map groundwater potential, and each has advantages and disadvantages.
This research developed a novel state-of-the-art approach using an ensemble of statistical, machine-learning, and MCDM approaches in conjunction with RS and GIS techniques to map groundwater potential. MCDM models, statistical models, and machine-learning algorithms (like VIKOR, FR, and RF, respectively) can be used to support scientific decisions for problems that have a diversity of management criteria [94]. VIKOR can analyze both quantitative and qualitative data, uses simple mathematical formulas, produces stable results, is flexible, is highly efficient, is very accurate, it can factor uncertainty into analyses, is comprehensive in its analysis, produces visual representations of the data, enables comparison of alternatives, and can incorporate and analyze large amounts of data. FR was used for its simplicity of equations, easily interpretable results, and its high efficiency [95]. RS and GIS save time and are cost-effective ways to improve data gathering and analysis [96,97,98,99].
Machine-learning algorithms can model natural phenomena (e.g., groundwater potential and groundwater presence) with nonlinear relationships. They do not require removal of outliers, data transformations, or statistical assumptions. They can handle complex nonlinear relationships between conditioning factors and groundwater potentiality and can automatically analyze effective-factor (i.e., predictor) interactions. RF was the machine-learning model employed in this research because of its advantages. The RF method can accommodate different types of predictor variables as well as missing values. It can assess the relationships and interactions between effective factors. It can handle uncertainty, can apply several input factors without eliminating any, and can return a small set of categories to maintain high prediction accuracy [54,67]. Classification accuracy, however, is affected by the number, scale, types, and precision of inputs. Incorporating all suitable factors increases modeling accuracy. Moreover, compared to other models, RF can include more datasets [53]. It can overfit some datasets with noisy classification and regression tasks. It can be biased toward attribute data with more categories, and it can also be biased toward attributes with more levels. Therefore, the variable-importance scores from RF are not reliable for such data. There are a few drawbacks to this algorithm, however. It does require significant computational memory, analysis of large datasets consumes a lot of time, there are temporal high costs for pruning, there are a high number of overlapping end nodes, and errors accumulate at each layer as trees grow.
The analysis of the spatial relationships between GWCFs and well locations using the FR model showed that elevation, slope, CI, SPI, distances to the nearest road, distances to the nearest stream, less rainfall, TRI, and NDVI are inversely correlated to groundwater potential. TWI and drainage density are directly correlated to well locations. These relationships are consistent with previous studies [29,33]. Elevation and relief play important roles in groundwater potential because weather and climatic conditions vary greatly between elevations and cause differential soil development and vegetation growth patterns. Areas at low elevations, with gentle slopes, and with low drainage density are usually associated with highly permeable material, higher vegetative cover, and less relief. These areas have higher groundwater potential. Conversely, higher elevations, steeper slopes, and higher drainage density usually have less soil, more impermeable subsurface material, sparse vegetation, and runoff-promoting relief. TWI, CI, TRI, and SPI topographic indices predict soil moisture, groundwater flow, and slope stability. These indices have been used to reflect spatial patterns of soil-moisture [77]. In terms of LULC and lithology, agricultural areas and regions with Quaternary formations have the strongest relationships to wells. Agriculture tends not to occur in places without access to water and wells are not often drilled in places that lack demand for water. The LULC categories covary with soil conditions and subsequently reflect groundwater occurrence [100]. Lithology influences groundwater because lithological variation is associated with varying porosity and hydraulic conductivity in rock formations and subsequent soils. Soil has an impact on recharge due to permeability factors and composition, and so it can influence the potential size of an aquifer [100].
The RF model indicated that LULC, lithology, and elevation were most strongly related to groundwater flow rates. This is consistent with other research [1,8,31,34,38,100]. Rahmati et al. [8] mapped groundwater potential in the Mehran Region of Iran by combining RF and maximum entropy models using data from 163 wells. They analyzed ten GWCFs and elevation, drainage density, lithology, and land use were found to be the greatest determinants of groundwater production. Arabameri et al. [1] combined RF, WoE, binary logistic regression, and the technique for order-preference by similarity to ideal solution (TOPSIS) for groundwater potential mapping on the Shahroud Plain in Iran. Analysis of the relationships between 122 wells and 15 GWCFs revealed that LULC, soil type, and slope were key correlates to groundwater. Naghibi et al. [38] used boosted regression tree, classification and regression tree, and RF machine-learning models for GWPM in Iran. Thirteen hydrological-geological-physiographical conditioning factors were assessed at 864 springs. Elevation and drainage density were the most important correlates to spring locations.
Comparing the goodness-of-fit and performance of the ensemble model using AUROC, FR, and SCAI methods revealed that the ensemble model had excellent prediction accuracy, and this was consistent with other research [30,43]. Combined or ensemble models are often more effective than are individual, stand-alone models. Kordestani et al. [30] showed that by combining EBF and BRT for GWPM, prediction accuracy (AUROC) climbed from 75% for EBF alone to 82% for EBF-BRT. Four data-mining models and FR were combined as an ensemble model [43] and showed better performance by the ensemble model because it reduced overfitting.
Overfitting occurs due to noise, limited training data, and classification complexity. Solutions are difficult because of the width of the array of potential causes. Three aspects of this study suggest that overfitting could have occurred: the area from which the non-well locations were selected was confined to the region that contained wells even though the study region extended well beyond the region covered by wells, the ratio of 70:30 training to validation datasets was used without an assessment of the accuracy of the sampling ratio, and despite the multicollinearity assessment, GWCFs still had noise.
A lack of confidence (uncertainty) might result from doubt about the completeness of knowledge. Data could be suspected of being unclear, inaccurate, unreliable, inconclusive, or potentially false [101,102,103,104]. Machine-learning techniques like RF can account for nonlinear relationships and can handle uncertainty in data. In an RF model, the prediction interval has two main components that determine its width: one represents the uncertainty of the predicted mean (or some other parameter)—this is the confidence interval—and the other represents the variability of individual observations around that mean. BRT uses uncertainty reduction to arrive at a final number of trees [105]. There are many sources of error in datasets used for modeling: measurement errors, sampling bias, limitations in field data collection, genetic variability, etc. These errors can affect model accuracy [106]. The strength of the novel ensemble model introduced here to model natural phenomena with nonlinear relationships is confirmed by the results of this study.
There are several things that can enhance GWPM accuracy and they should be the foci of future research. First, there should be a complete analysis of the implications of limiting or extending the area from which the selection of non-well locations is made to places where wells are located or to areas beyond the distribution of wells. Second, the conventional 70:30 ratio of the datasets used for training and validation should be revisited, as proportion may be causing issues. Third, the methods used to select features may be having unseen effects on the modeling process, and fourth, there should be a comparison of the VIKOR-RF-FR ensemble to other ensemble models to predict groundwater locations.

5. Conclusions

Many methods have been used to assess regional groundwater potential. This study created an MCDM statistical machine-learning ensemble model to spatially predict groundwater in Bastam Watershed, Iran. VIKOR, FR, and RF models were combined to map and assess the characteristics of groundwater-bearing landscapes. Like many other parts of the world, the Bastam Watershed suffers from a shortage of water. Seventeen conditioning factors (elevation, slope, aspect, SPI, TPI, TWI, TRI, CI, distance to stream, distance to fault, distance to road, drainage density, rainfall, soil type, LU/LC, lithology, and NDVI) were assessed at 89 well locations to determine each conditioning factors’ statistical association with groundwater presence and productivity. RF modeling revealed that LULC, lithology, and elevation were the most important factors predicting groundwater potential and production. Validation confirmed that the ensemble model had high prediction accuracy. The VIKOR-RF-FR ensemble model is a systematic, objective, and scientific decision-making model that can be used to support groundwater resource development and management plans. A GWPM can help planners guide future groundwater exploration. The same analysis in other districts with similar topographic and geological conditions can save time and money in their quests for protecting and developing groundwater resources. A limitation of this research is that a fixed combination of GWCFs were considered for modeling. Different combinations of GWCFs should be explored in order to further improve groundwater potential modeling.

Author Contributions

Conceptualization, A.A.; Data curation, A.A.; Formal analysis, A.A.; Investigation, A.A. S.L., and P.T.T.N.; Methodology, A.A., S.L., J.P.T., and P.T.T.N.; Resources, A.A.; Software, A.A.; Supervision, A.A., S.L., and P.T.T.N.; Validation, A.A.; Writing—original draft, A.A.; Writing—review and editing, A.A., J.P.T., S.L., and P.T.T.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This research was conducted by the Basic Research Project of the Korea Institute of Geoscience and Mineral Resources (KIGAM) funded by the Ministry of Science and ICT.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arabameri, A.; Rezaei, K.; Cerda, A.; Lombardo, L.; Rodrigo-Comino, J. GIS-based groundwater potential mapping in Shahroud plain, Iran. A comparison among statistical (bivariate and multivariate), data mining and MCDM approaches. Sci. Total Environ. 2019, 658, 160–177. [Google Scholar] [CrossRef]
  2. Foster, S. Groundwater assessing vulnerability and promoting protection of a threatened resource. In Proceedings of the 8th Stockholm Water Symposium, Stockholm, Sweden, 10–13 August 1998. [Google Scholar]
  3. UN. Water for people, water for life. In The UN World Water Development Report (WWDR), UNESCO; Publishing and Berghahn Books: New York, NY, USA, 2003; p. 34. [Google Scholar]
  4. ElNaqa, A.; AlShayeb, A. Groundwater protection and management strategy in jordan. Water Resour. Manag. 2008, 23, 2379–2394. [Google Scholar] [CrossRef]
  5. Schematization and Management Organ of Iran (SMOI). 2004. Available online: http://www.ncc.org.ir/ (accessed on 12 August 2018).
  6. Magesh, N.S.; Chandrasekar, N.; John, 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] [Green Version]
  7. Lee, S.; Hyun, Y.; Lee, M.J. Groundwater Potential Mapping Using Data Mining Models of Big Data Analysis in Goyang-si, South Korea. Sustainability 2019, 11, 1678. [Google Scholar] [CrossRef] [Green Version]
  8. Rahmati, O.; Pourghasemi, H.R.; Melesse, A.M.N. Application of GIS-based data driven random forest and maximum entropy models for groundwater potential mapping: A case study at Mehran region, Iran. Catena 2016, 137, 360–372. [Google Scholar] [CrossRef]
  9. 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]
  10. Mondal, S. Remote sensing and GIS based ground water potential mapping of kangshabati irrigation command area, west bengal. Geogr. Nat. Disasters 2012, 1, 1–8. [Google Scholar] [CrossRef]
  11. Oikonomidis, D.; Dimogianni, S.; Kazakis, N.; Voudouris, K. A GIS/Remote Sensing-based methodology for groundwater potentiality assessment in Tirnavos area, Greece. J. Hydrol. 2015, 525, 197–208. [Google Scholar] [CrossRef]
  12. Thomas, A.; Sharma, P.K.; Sharma, M.K.; Sood, A. Hydrogeomorphological mapping in assessing groundwater by using remote sensing datada case study in Lehra Gage Block, Sangrur district, Punjab. J. Indian Soc. Remote Sens. 1999, 27, 31–42. [Google Scholar] [CrossRef]
  13. Muralidhar, M.; Raju, K.R.K.; Raju, K.S.V.P.; Prasad, J.R. Remote sensing applications for the evaluation of water resources in rainfed area, Warangal district, Andhra Pradesh. Indian Miner. 2003, 34, 33–40. [Google Scholar]
  14. Pinto, D.; Shrestha, S.; Babel, M.S.; Ninsawat, S. Delineation of groundwater potential zones in the Comoro watershed, Timor Leste using GIS, remote sensing and analytic hierarchy process (AHP) technique. Appl. Water Sci. 2017, 7, 503–519. [Google Scholar] [CrossRef] [Green Version]
  15. Magaia, L.A.; Goto, T.N.; Masoud, A.A.; Koike, K. Identifying groundwater potential in crystalline basement rocks using remote sensing and electromagnetic sounding techniques in Central Western Mozambique. Nat. Resour. Res. 2018, 27, 275–298. [Google Scholar] [CrossRef]
  16. Luís, A.M. Development of Regional Exploration Techniques for Groundwater Resources in Semiarid Areas Through Integration of Remote Sensing and Geophysical Survey. Ph.D. Thesis, Kyoto University, Kyoto, Japan, 2018. [Google Scholar]
  17. Pradhan, B. Groundwater potential zonation for basaltic watersheds using satellite remote sensing data and GIS techniques. Open Geosci. 2009, 1, 120–129. [Google Scholar] [CrossRef]
  18. Bera, K.; Bandyopadhyay, J. Ground water potential mapping in Dulung watershed using remote sensing & GIS techniques, West Bengal, India. Int. J. Sci. Res. Publ. 2012, 2, 1–7. [Google Scholar]
  19. Al-Ruzouq, R.; Shanableh, A.; Merabtene, T.; Siddique, M.; Khalil, M.A.; Idris, A.; Almulla, E. Potential groundwater zone mapping based on geo-hydrological considerations and multi-criteria spatial analysis: North UAE. CATENA 2019, 173, 511–524. [Google Scholar] [CrossRef]
  20. Shi, Z.H.; Cai, C.F.; Ding, S.W.; Wang, T.W.; Chow, T.L. Soil conservation planning at the small watershed level using RUSLE with GIS: A case study in the Three Gorge area of China. Catena 2004, 55, 33–48. [Google Scholar] [CrossRef]
  21. Sharma, S.; Kumar Mahajan, A. Comparative evaluation of GIS-based landslide susceptibility mapping using statistical and heuristic approach for Dharamshala region of Kangra Valley, India. Geoenviron. Disasters 2018, 5, 4. [Google Scholar] [CrossRef] [Green Version]
  22. Saro, L. Current and future status of GIS-based landslide susceptibility mapping: A literature review. Korea J. Remote Sens. 2019, 35, 179–193. [Google Scholar]
  23. Saro, L.; Oh, H.J. Landslide Susceptibility Prediction using Evidential Belief Function, Weight of Evidence and Artificial Neural Network Models. Korea J. Remote Sens. 2019, 35, 299–316. [Google Scholar]
  24. Arabameri, A.; Pourghasemi, H.R.; Yamani, M. Applying different scenarios for landslide spatial modeling using computational intelligence methods. Environ. Earth Sci. 2017, 76, 832. [Google Scholar] [CrossRef]
  25. Singh, S.K.; Srivastava, K.; Gupta, M.; Thakur, K.; Mukherjee, S. Appraisal of land use/land cover of mangrove forest ecosystem using support vector machine. Environ. Earth Sci. 2014, 71, 2245–2255. [Google Scholar] [CrossRef]
  26. Mahato, S.; Pal, S. Groundwater Potential Mapping in a Rural River Basin by Union (OR) and Intersection (AND) of Four Multi-criteria Decision-Making Models. Nat. Resour. Res. 2019, 28, 523–545. [Google Scholar] [CrossRef]
  27. Lee, S.; Kim, Y.S.; Oh, H.J. Application of a weights of evidence method and GIS to regional groundwater productivity potential mapping. J. Environ. Manag. 2012, 96, 91–105. [Google Scholar] [CrossRef] [PubMed]
  28. Park, I.; Kim, Y.; Lee, S. Groundwater productivity potential mapping using evidential belief function. Groundwater 2014, 52, 201–207. [Google Scholar] [CrossRef]
  29. Mogaji, K.A.; Omosuyi, G.O.; Adelusi, A.O.; Lim, H.S. Application of GIS-based evidential belief function model to regional groundwater recharge potential zones mapping in hardrock geologic terrain. Environ. Process. 2016, 3, 93–123. [Google Scholar] [CrossRef]
  30. Kordestani, M.D.; Naghibi, S.A.; Hashemi, H.; Ahmadi, K.; Kalantar, B.; Pradhan, B. Groundwater potential mapping using a novel data-mining ensemble model. Hydrogeol. J. 2019, 27, 211–224. [Google Scholar] [CrossRef] [Green Version]
  31. Naghibi, S.A.; Pourghasemi, H.R.; Pourtaghie, Z.S.; Rezaei, A. Groundwater qanat potential mapping using frequency ratio and Shannon’s entropy models in the Moghan Watershed, Iran. Earth Sci. Inf. 2015, 8, 171–186. [Google Scholar] [CrossRef]
  32. Golkarian, A.; Naghibi, S.A.; Kalantar, B.; Pradhan, B. Groundwater potential mapping using C5.0, random forest, and multivariate adaptive regression spline models in GIS. Environ. Monit. Assess. 2018, 190, 149. [Google Scholar] [CrossRef]
  33. Trabelsi, F.; Lee, S.; Khlifi, S.; Arfaoui, A. Frequency Ratio Model for Mapping Groundwater Potential Zones Using GIS and Remote Sensing; Medjerda Watershed Tunisia. In Advances in Sustainable and Environmental Hydrology, Hydrogeology, Hydrochemistry and Water Resources; Chaminé, H., Barbieri, M., Eds.; Springer: Berlin/Heidelberg, Germany, 2019. [Google Scholar]
  34. Razandi, Y.; Pourghasemi, H.R.; SamaniNeisani, N.; Rahmati, O. Application of analytical hierarchy process, frequency ratio, and certainty factor models for groundwater potential mapping using GIS. Earth Sci. Inf. 2015, 8, 867–883. [Google Scholar] [CrossRef]
  35. Chen, W.; Li, H.; Hou, E.; Wang, S.; Wang, G.; Panahi, M. GIS-based groundwater potential analysis using novel ensemble weights-of-evidence with logistic regression and functional tree models. Sci. Total Environ. 2018, 1, 853–867. [Google Scholar] [CrossRef] [Green Version]
  36. Lee, S.; Lee, C.W.; Kim, J.C. Groundwater Productivity Potential Mapping Using Logistic Regression and Boosted Tree Models: The Case of Okcheon City in Korea. In Advances in Remote Sensing and Geo Informatics Applications; El-Askary, H., Lee, S., Eds.; Springer: Berlin/Heidelberg, Germany, 2019. [Google Scholar]
  37. Golkarian, A.; Rahmati, O. Use of a maximum entropy model to identify the key factors that influence groundwater availability on the Gonabad Plain, Iran. Environ. Earth Sci. 2018, 77, 369. [Google Scholar] [CrossRef]
  38. Naghibi, S.A.; Pourghasemi, H.R.; Dixon, B. GIS-based groundwater potential mapping using boosted regression tree, classification and regression tree, and random forest machine learning models in Iran. Environ. Monit. Assess. 2016, 188, 44. [Google Scholar] [CrossRef] [PubMed]
  39. Lee, S.; Song, K.Y.; Kim, Y.; Park, I. Regional groundwater productivity potential mapping using a geographic information system (GIS) based artificial neural network model. Hydrogeol. J. 2012, 20, 1511–1527. [Google Scholar] [CrossRef]
  40. Miraki, S.; Hedayati Zanganeh, S.; Chapi, K.; Singh, V.P.; Shirzadi, A.; Shahabi, H.; Thai Pham, B. Mapping Groundwater Potential Using a Novel Hybrid Intelligence Approach. Water Resour. Manag. 2019, 33, 281–302. [Google Scholar] [CrossRef]
  41. Shahid, S.; Nath, S.K.; Kamal, A.S. GIS integration of remote sensing and topographic data using fuzzy logic for ground water assessment in Midnapur district, India. Geocarto Int. 2014, 17, 69–74. [Google Scholar] [CrossRef]
  42. Chen, W.; Pradhan, B.; Li, S.; Shahabi, H.; Mojaddadi Rizeei, H.; Hou, E.; Wang, S. Novel Hybrid Integration Approach of Bagging-Based Fisher’s Linear Discriminant Function for Groundwater Potential Analysis. Nat. Resour. Res. 2019, 28, 1239–1258. [Google Scholar] [CrossRef] [Green Version]
  43. Naghibi, S.A.; Ahmadi, K.; Daneshi, A. Application of support vector machine, random forest, and genetic algorithm optimized randomforest models in groundwater potential mapping. Water Resour. Manag. 2017, 31, 2761–2775. [Google Scholar] [CrossRef]
  44. Rahmati, O.; Naghibi, S.A.; Shahabi, H.; Tien Bui, D.; Pradhan, B.; Aareh, A. Groundwater spring potential modelling: Comprising the capability and robustness of three different modeling approaches. Hydrology 2018, 565, 248–261. [Google Scholar] [CrossRef]
  45. Naghibi, S.A.; Pourghasemi, H.R.; Abbaspour, K. A comparison between ten advanced and soft computing models for groundwater qanat potential assessment in Iran using R and GIS. Theor. Appl. Clim. 2018, 131, 967–984. [Google Scholar] [CrossRef]
  46. Naghibi, S.A.; Dashtpagerdi, M.M. Evaluation of four supervised learning methods for groundwater spring potential mapping in Khalkhal region (Iran) using GIS-based features. Hydrogeol. J. 2017, 25, 169–189. [Google Scholar] [CrossRef]
  47. Naghibi, S.A.; Moghaddam, D.D.; Kalantar, B.; Pradhan, B.; Kisi, O. A comparative assessment of GIS-based data mining models and a novel ensemble model in groundwater well potential mapping. J. Hydrol. 2017, 548, 471–483. [Google Scholar] [CrossRef]
  48. Razavi-Termeh, S.V.; Sadeghi-Niaraki, A.; Choi, S.-M. 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] [Green Version]
  49. Alganci, U.; Besol, B.; Sertel, E. Accuracy Assessment of Different Digital Surface Models. ISPRS Int. J. Geo Inf. 2018, 7, 114. [Google Scholar] [CrossRef] [Green Version]
  50. Arabameri, A.; Pourghasemi, H.R.; Cerda, A. Erodibility prioritization of sub-watersheds using morphometric parameters analysis and its mapping: A comparison among TOPSIS, VIKOR, SAW, and CF multi-criteria decision making models. Sci. Total Environ. 2017, 613, 1385–1400. [Google Scholar]
  51. Arabameri, A.; Pradhan, B.; Rezaei, K.; Yamani, M.; Pourghasemi, H.R.; Lombardo, L. Spatial modelling of gully erosion using Evidential Belief Function, Logistic Regression and a new ensemble EBF–LR algorithm. Land Degrad. Dev. 2018, 29, 4035–4049. [Google Scholar] [CrossRef]
  52. Arabameri, A.; Rezaei, K.; Pourghasemi, H.R.; Lee, S.; Yamani, M. GIS-based gully erosion susceptibility mapping: A comparison among three data-driven models and AHP knowledge-based technique. Environ. Earth Sci. 2018, 77, 628. [Google Scholar] [CrossRef]
  53. Arabameri, A.; Pradhan, B.; Pourghasemi, H.R.; Rezaei, K. Identification of erosion-prone areas using different multi-criteria decision-making techniques and GIS. Geomat. Nat. Hazards Risk 2018, 9, 1129–1155. [Google Scholar] [CrossRef] [Green Version]
  54. Arabameri, A.; Rezaei, K.; Cerda, A.; Conoscenti, C.; Kalantari, Z. A comparison of statistical methods and multi-criteria decision making to map flood hazard susceptibility in Northern Iran. Sci. Total Environ. 2019, 660, 443–458. [Google Scholar] [CrossRef]
  55. Arabameri, A.; Pradhan, B.; Rezaei, K.; Conoscenti, C. Gully erosion susceptibility mapping using GIS-based multi-criteria decision analysis techniques. CATENA 2019, 180, 282–297. [Google Scholar] [CrossRef]
  56. Arabameri, A. Application of the Analytic Hierarchy Process (AHP) for locating fire stations: Case Study Maku City. Merit Res. J. ArtSoc. Sci. Humanit. 2014, 2, 1–10. [Google Scholar]
  57. Arabameri, A.; Ramesht, M.H. Site Selection of Landfill with emphasis on Hydrogeomorphological–environmental parameters Shahrood-Bastam watershed. Sci. J. Manag. Syst. 2017, 16, 55–80. [Google Scholar]
  58. Yamani, M.; Arabameri, A. Comparison and evaluation of three methods of multi attribute decision making methods in choosing the best plant species for environmental management (Case study: Chah Jam Erg). Nat. Environ. Chang. 2015, 1, 49–62. [Google Scholar]
  59. Arabameri, A. Zoning Mashhad Watershed for Artificial Recharge of Underground Aquifers Using Topsis Model and GIS Technique. Glob. J. Hum. Soc. Sci. B Geogr. Geo Sci. Environ. Disaster Manag. 2014, 14, 45–53. [Google Scholar]
  60. Mousavi, S.M.; Golkarian, A.; Naghibi, S.A.; Kalantar, B.; Pradhan, B. GIS-based groundwater spring potential mapping using data mining boosted regression tree and probabilistic frequency ratio models in Iran. Aims Geosci. 2017, 3, 91–115. [Google Scholar]
  61. Siahkamari, S.; Haghizadeh, A.; Zeinivand, H.; Tahmasebipour, N.; Rahmati, O. Spatial prediction of flood-susceptible areas using frequency ratio and maximum entropy models. Geocarto Int. 2018, 33, 927–941. [Google Scholar] [CrossRef]
  62. Arabameri, A.; Pradhan, B.; Rezaei, K.; Sohrabi, M.; Kalantari, Z. GIS-based landslide susceptibility mapping using numerical risk factor bivariate model and its ensemble with linear multivariate regression and boosted regression tree algorithms. J. Mt. Sci. 2019, 16, 595–618. [Google Scholar] [CrossRef]
  63. Arabameri, A.; Pradhan, B.; Rezaei, K.; Saro, L.; Sohrabi, M. An Ensemble Model for Landslide Susceptibility Mapping in a Forested Area. Geocarto Int. 2019, 1–26. [Google Scholar] [CrossRef]
  64. Arabameri, A.; Pradhan, B.; Rezaei, K. Spatial prediction of gully erosion using ALOS PALSAR data and ensemble bivariate and data mining models. Geosci. J. 2019, 24, 669–686. [Google Scholar] [CrossRef]
  65. Arabameri, A.; Pradhan, B.; Rezaei, K.; Lee, C.-W. Assessment of Landslide Susceptibility Using Statistical- and Artificial Intelligence-Based FR–RF Integrated Model and Multiresolution DEMs. Remote Sens. 2019, 11, 999. [Google Scholar] [CrossRef] [Green Version]
  66. Arabameri, A.; Pradhan, B.; Pourghasemi, H.R.; Rezaei, K.; Kerle, N. Spatial Modelling of Gully Erosion Using GIS and R Programing: A Comparison among Three Data Mining Algorithms. Appl. Sci. 2018, 8, 1369. [Google Scholar] [CrossRef] [Green Version]
  67. Arabameri, A.; Pourghasemi, H.R. Spatial Modeling of Gully Erosion Using Linear and Quadratic Discriminant Analyses in GIS and R. In Spatial Modeling in GIS and R for Earth and Environmental Sciences, 1st ed.; Pourghasemi, H.R., Gokceoglu, C., Eds.; Elsevier: Amsterdam, The Netherlands, 2019; p. 796. [Google Scholar]
  68. Arabameri, A.; Pradhan, B.; Rezaei, K. Gully erosion zonation mapping using integrated geographically weighted regression with certainty factor and random forest models in GIS. J. Environ. Manag. 2019, 232, 928–942. [Google Scholar] [CrossRef]
  69. IRIMO. Summary Reports of Iran’s Extreme Climatic Events. In Ministry of Roads and Urban Development; Iran Meteorological Organization: Tehran, Iran, 2012; Available online: http://www.cri.ac.ir (accessed on 12 August 2018).
  70. GSI. Geology Survey of Iran. 1997. Available online: http://www.gsi.ir/Main/Lang_en/index.html (accessed on 12 August 2018).
  71. IUSSWorking Group WRB14. World Reference Base for Soil Resources 2014, World Soil Resources Report; FAO: Rome, Italy, 2014. [Google Scholar]
  72. Noor, H. Analysis of Groundwater Resource Utilization and Their Current Condition in Iran. Iran. J. Rainwater Catchment Syst. 2018, 5, 29–38. [Google Scholar]
  73. Ramesht, M.H.; Arabameri, A. Shahrood-Bastam Basin Zoning for the Purpose of Artificial Underground Aquifer Recharge by Using Linear Assignment Method and GIS Technique. Geogr. Space 2013, 12, 134–149. [Google Scholar]
  74. Chen, W.; Xie, X.; Wang, J.; Pradhan, B.; Hong, H.; Bui, D.T. A comparative study of logistic model tree, random forest, and classification and regression tree models for spatial prediction of landslide susceptibility. Catena 2017, 151, 147–160. [Google Scholar] [CrossRef] [Green Version]
  75. Liang, C.-P.; Chen, J.-S.; Chien, Y.-C.; Chen, C.-F. Spatial analysis of the risk to human health from exposure to arsenic contaminated groundwater: A kriging approach. Sci. Total Environ. 2018, 627, 1048–1057. [Google Scholar] [CrossRef]
  76. Haghizadeh, A.; Davoudi Moghadam, D.; Pourghasemi, H.R. GIS-based bivariate statistical techniques for groundwater potential analysis (an example of Iran). J. Earth Syst. Sci. 2017, 126, 109. [Google Scholar] [CrossRef] [Green Version]
  77. Moore, I.D.; Grayson, R.B.; Ladson, A.R. Digital terrain modelling: A review of hydrological, geomorphological, and biological applications. Hydrol. Process. 1991, 5, 3–30. [Google Scholar] [CrossRef]
  78. De Reu, J.; Bourgeois, J.; Bats, M.; Zwertvaegher, A.; Gelorini, V.; De Smedt, P.; Chu, W.; Antrop, M.; De Maeyer, P.; Finke, P.; et al. Application of the topographic position index to heterogeneous landscapes. Geomorphology 2013, 186, 39–49. [Google Scholar] [CrossRef]
  79. Althuwaynee, O.F.; Pradhan, B.; Park, H.J.; Lee, J.H. A novel ensemble bivariate statistical evidential belief function with knowledge-based analytical hierarchy process and multivariate statistical logistic regression for landslide susceptibility mapping. Catena 2014, 114, 21–36. [Google Scholar] [CrossRef]
  80. Jenness, J. Surface Areas and Ratios from Elevation Grid; Jenness Enterprises: Flagstaff, AZ, USA, 2012. [Google Scholar]
  81. Oh, H.J.; Pradhan, B. Application of a neuro-fuzzy model to landslide-susceptibility mapping for shallow landslides in a tropical hilly area. Comput. Geosci. 2011, 37, 1264–1276. [Google Scholar] [CrossRef]
  82. Hengl, T.; Gruber, S.; Shrestha, D.P. Digital terrain analysis in ILWIS. International Institute for Geo-Information Science and Earth Observation Enschede. Int. Inst. Geoinf. Sci. Earth Obs. Enschede Neth. 2003, 62, 1–56. [Google Scholar]
  83. Yesilnacar, E.K. The Application of Computational Intelligence to Landslide Susceptibility Mapping in Turkey. Ph.D. Thesis, Department of Geomatics the University of Melbourne, Melbourne, Australia, 2005; p. 423. [Google Scholar]
  84. Bonham-Carter, G.F. Geographic information systems for geoscientists: Modeling with GIS. In Computer Methods in the Geosciences; Pergamon: Bergama, Turkey, 1994. [Google Scholar]
  85. Shafapour Tehrany, M.; Kumar, L. The application of a Dempster–Shafer-based evidential belief function in flood susceptibility mapping and comparison with frequency ratio and logistic regression methods. Environ. Earth Sci. 2018, 77, 490. [Google Scholar] [CrossRef]
  86. Opricovic, S. Multicriteria optimization of civil engineering systems. Fac. Civ. Eng. Belgrade 1998, 2, 5–21. [Google Scholar]
  87. Opricovic, S.; Tzeng, G.H. Multicriteria planning of post-earthquake sustainable reconstruction. Comput. Aided Civ. Infrastruct. Eng. 2002, 17, 211–220. [Google Scholar] [CrossRef]
  88. Breiman, L. Random forests. J. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  89. Hastie, T. The Elements of Statistical Learning: Data Mining, Inference, and Prediction, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2001; p. 533. [Google Scholar]
  90. Swets, J.A. Measuring the accuracy of diagnostic systems. Science 1988, 240, 1285–1293. [Google Scholar] [CrossRef] [Green Version]
  91. Pradhan, B. Remote sensing and GIS-based landslide hazard analysis and cross-validation using multivariate logistic regression model on three test areas in Malaysia. Adv. Space Res. 2010, 45, 1244–1256. [Google Scholar] [CrossRef]
  92. Yilmaz, C.; Topal, T.; Suzen, M.L. GIS-based landslide susceptibility mapping using bivariate statistical analysis in Devrek (Zonguldak-Turkey). Environ. Earth Sci. 2012, 65, 2161–2178. [Google Scholar] [CrossRef]
  93. 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. 2012, 7, 711–724. [Google Scholar] [CrossRef]
  94. Kiker, G.A.; Bridges, T.S.; Varghese, A.; Seager, T.P.; Linkovjj, I. Application of multi-criteria decision analysis in environmental decision making. Integr. Environ. Assess. Manag. 2005, 1, 95–108. [Google Scholar] [CrossRef] [PubMed]
  95. Aditian, A.; Kubota, T.; Shinohara, Y. Comparison of GIS-based landslide susceptibility models using frequency ratio, logistic regression, and artificial neural network in a tertiary region of Ambon, Indonesia. Geomorphology 2018, 318, 101–111. [Google Scholar] [CrossRef]
  96. 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]
  97. Arkoprovo, B.; Adarsa, J.; Prakash, S.S. Delineation of groundwater potential zones using satellite remote sensing and geographic information system techniques: A case study from Ganjam district, Orissa, India. Res. J. Recent. Sci. 2012, 1, 59–66. [Google Scholar]
  98. Hammouri, N.; ElNaqa, A.; Barakat, M. An integrated approach to groundwater exploration using remote sensing and geographic information system. J. Water Resour. Prot. 2012, 4, 717–724. [Google Scholar] [CrossRef] [Green Version]
  99. Davoodi, M.D.; Rezaei, M.; Pourghasemi, H.R.; Pourtaghi, Z.S.; Pradhan, B. Groundwater spring potential mapping using bivariate statistical model and GIS in the Taleghan watershed Iran. Arab. J. Geosci. 2015, 8, 913–929. [Google Scholar]
  100. 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]
  101. Paris, G.; Robilliard, D.; Fonlupt, C. Exploring Overfitting in Genetic Programming. In Proceedings of the International Conference Evolution Artificielle, Marseilles, France, 17 October 2004; pp. 267–277. [Google Scholar]
  102. Arabameri, A.; Chen, W.; Lombardo, L.; Blaschke, T.; Tien Bui, D. Hybrid Computational Intelligence Models for Improvement Gully Erosion Assessment. Remote Sens. 2020, 12, 140. [Google Scholar] [CrossRef] [Green Version]
  103. Arabameri, A.; Blaschke, T.; Pradhan, B.; Pourghasemi, H.R.; Tiefenbacher, J.P.; Bui, D.T. Evaluation of Recent Advanced Soft Computing Techniques for Gully Erosion Susceptibility Mapping: A Comparative Study. Sensors 2020, 20, 335. [Google Scholar] [CrossRef] [Green Version]
  104. Arabameri, A.; Chen, W.; Blaschke, T.; Tiefenbacher, J.P.; Pradhan, B.; Tien Bui, D. Gully Head-Cut Distribution Modeling Using Machine Learning Methods—A Case Study of N.W. Iran. Water 2020, 12, 16. [Google Scholar] [CrossRef] [Green Version]
  105. Pourtaghi, Z.S.; Pourghasemi, H.R. GIS-based groundwater spring potential assessment and mapping in the Birjand Township, southern Khorasan Province, Iran. Hydrogeol. J. 2015, 22, 643–662. [Google Scholar] [CrossRef]
  106. Sahoo, S.; Munusamy, S.B.; Dhar, A.; Kar, A.; Ram, P. Appraising the accuracy of multi-class frequency ratio and weights of evidence method for delineation of regional groundwater potential zones in canal command system. Water Resour. Manag. 2017, 31, 4399–4413. [Google Scholar] [CrossRef]
Figure 1. Location of study area in Semnan Province and Iran (inset) and the locations of training and validation wells in the study area.
Figure 1. Location of study area in Semnan Province and Iran (inset) and the locations of training and validation wells in the study area.
Remotesensing 12 00490 g001
Figure 2. Flowchart of the research.
Figure 2. Flowchart of the research.
Remotesensing 12 00490 g002
Figure 3. Location of random points in the study area.
Figure 3. Location of random points in the study area.
Remotesensing 12 00490 g003
Figure 4. Groundwater conditioning factors. (a) elevation, (b) slope, (c) aspect, (d) stream power index (SPI), (e) topography position index (TPI), (f) topography wetness index (TWI), (g) terrain ruggedness index (TRI), (h) convergence index (CI), (i) distance to stream, (j) distance to road, (k) distance to fault, (l) drainage density, (m) rainfall, (n) soil type, (o) land use/land cover (LULC), (p) lithology, (q) normalized difference vegetation index (NDVI).
Figure 4. Groundwater conditioning factors. (a) elevation, (b) slope, (c) aspect, (d) stream power index (SPI), (e) topography position index (TPI), (f) topography wetness index (TWI), (g) terrain ruggedness index (TRI), (h) convergence index (CI), (i) distance to stream, (j) distance to road, (k) distance to fault, (l) drainage density, (m) rainfall, (n) soil type, (o) land use/land cover (LULC), (p) lithology, (q) normalized difference vegetation index (NDVI).
Remotesensing 12 00490 g004aRemotesensing 12 00490 g004b
Figure 5. Relative importance of conditioning factors using random forest model.
Figure 5. Relative importance of conditioning factors using random forest model.
Remotesensing 12 00490 g005
Figure 6. Weighted GWCFs using the FR model for (a) elevation, (b) slope, (c) aspect, (d) stream power index (SPI), (e) topography position index (TPI), (f) topography wetness index (TWI), (g) terrain ruggedness index (TRI), (h) convergence index (CI), (i) distance to stream, (j) distance to road, (k) distance to fault, (l) drainage density, (m) rainfall, (n) soil type, (o) LU/LC, (p) lithology, (q) NDVI.
Figure 6. Weighted GWCFs using the FR model for (a) elevation, (b) slope, (c) aspect, (d) stream power index (SPI), (e) topography position index (TPI), (f) topography wetness index (TWI), (g) terrain ruggedness index (TRI), (h) convergence index (CI), (i) distance to stream, (j) distance to road, (k) distance to fault, (l) drainage density, (m) rainfall, (n) soil type, (o) LU/LC, (p) lithology, (q) NDVI.
Remotesensing 12 00490 g006aRemotesensing 12 00490 g006b
Figure 7. Groundwater potential map using the VIKOR-RF ensemble model.
Figure 7. Groundwater potential map using the VIKOR-RF ensemble model.
Remotesensing 12 00490 g007
Figure 8. Percentage of potential classes.
Figure 8. Percentage of potential classes.
Remotesensing 12 00490 g008
Figure 9. Area under curve (AUC). (a) Prediction rate curve (PRC), (b) Success rate curve (SRC).
Figure 9. Area under curve (AUC). (a) Prediction rate curve (PRC), (b) Success rate curve (SRC).
Remotesensing 12 00490 g009
Figure 10. Trend of frequency ratio and seed cell area index indicators.
Figure 10. Trend of frequency ratio and seed cell area index indicators.
Remotesensing 12 00490 g010
Table 1. Minimum and maximum values and classes of groundwater conditioning factors (GWDFs).
Table 1. Minimum and maximum values and classes of groundwater conditioning factors (GWDFs).
FactorsMinMaxClassesMethods
Elevation (m)135738931) <1659, 2) 1659–1990, 3) 1990–2342, 4) 2342–2724, 5) 2724–3163, 6) >3163Natural break (Jenks)
Slope (°)070.661) <5, 2) 5–10, 3) 10–15, 4) 15–20, 5) 20–30, 6) >30Natural break (Jenks)
Aspect−1337.51) Flat (−1), 2) North (0–22.5), 3) Northeast (22.5–67.5), 4) East (67.5–112.5), 5) Southeast (112.5–157.5), 6) South (157.5–202.5), 7) Southwest (202.5–247.5), 8) West (247.5–292.5), 9) Northwest (292.5–337.5)Directional units
SPI6.2723.051) <9.16, 2) 9.16–11.07, 3) 11.07–12.85, 4) 12.85–15.35, 5) >15.35Natural break (Jenks)
TPI−63.270.61) <−10.74, 2) −10.74–−3.3, 3) −3.3–2.9, 4) 2.9 0 11.8, 5) >11.8Natural break (Jenks)
TWI1.2021.041) <4.86, 2) 4.86–7.27, 3) 7.27–10.85, 4) >10.85Natural break (Jenks)
TRI063.11) <3.22, 2) 3.22–7.43, 3) 7.43–11.64, 4) 11.64–17.59, 5) >17.59Natural break (Jenks)
CI (100/m)−1001001) <−39.6, 2) −39.6–−11.3, 3) −11.3–9.01, 4) 9.01–38.03, 5) >38.03Natural break (Jenks)
Distance to stream (m)02317.51) <100, 2) 100–200, 3) 200–300, 4) 300–400, 5) >400Manual
Distance to road (m)011,5971) <500, 2) 500–1000, 3) 1000–1500, 4) 1500–2000, 5) >2000Manual
Distance to fault(m)016,1461) <500, 2) 500–1000, 3) 1000–1500, 4) 1500–2000, 5) >2000Manual
Drainage density (km/km2)0.153.181) <0.85, 2) 0.85–1.49, 3) 1.49–2.19, 4) >2.19Natural break (Jenks)
Rainfall (mm)1577221) <248.5, 2) 248.5–350.3, 3) 350.3–461.02, 4) 461.02–582.7, 5) >582.7Natural break (Jenks)
Soil type--1) Entisols, 2) Alfisols, 3) Aridisols, 4) Entisols/Inceptisols, 5) Inceptisols, 6) MollisolsSoil types/Orders
LULC--1) Agriculture, 2) Dense forest, 3) Good range, 4) Agri-dry farming, 5) Dry farming, 6) Low forest, 7) Woodland, 8) Mod-forest, 9) Mod-range, 10) Poor-range, 11) Rock, 12) UrbanSupervised
Classification
Lithology--1) A, 2) B, 3) C, 4) D, 5) E, 6) F, 7) G, 8) H, 9) I, 10) J.Lithological
Units
NDVI−0.240.541) <−0.01, 2) −0.01–0.07, 3) 0.07–0.12,
4) 0.12–0.21, 5) 0.21–0.32, 6) >0.21
Natural break (Jenks)
Table 2. Multi-collinearity test among groundwater conditioning factors.
Table 2. Multi-collinearity test among groundwater conditioning factors.
FactorsCollinearity StatisticsFactorsCollinearity Statistics
ToleranceVIF *ToleranceVIF
NDVI0.61.4Rainfall0.24.3
Distance to fault0.71.4Distance to road0.42.1
TRI0. 24.7Slope0.14.8
Slope aspect0.81.1Soil type0.32.5
CI0.71.3SPI0.51.7
Drainage density0.32.7Distance to stream0.61.4
Elevation0.33.2TPI0.32.7
Lithology0.23.4TWI0.32.7
LU/LC0.23.9
* VIF is variance inflation factor.
Table 3. Confusion matrix from the RF model.
Table 3. Confusion matrix from the RF model.
No Well (0)Well (1)
no well (0)746
well (1)1169
Table 4. Spatial relation between groundwater conditioning factors and wells using the frequency ratio model.
Table 4. Spatial relation between groundwater conditioning factors and wells using the frequency ratio model.
FactorClassNo. of Pixels%No. of Wells%FR
Elevation (m)<1659479,81432.485385.482.63
1659–1990309,67120.96914.520.69
1990–2342239,85416.2400.000.00
2342–2724209,10914.1600.000.00
2724–3163146,9129.9500.000.00
>316391,7926.2100.000.00
Slope (°)< 5588,36339.8362100.002.51
5–10149,25910.1000.000.00
10–15124,1598.4100.000.00
15–20143,8749.7400.000.00
20–30280,85919.0100.000.00
>30190,63712.9100.000.00
AspectF62620.4223.237.61
N132,3048.9600.000.00
NE152,50910.321422.582.19
E245,19816.601117.741.07
SE341,07123.091727.421.19
S273,88018.541016.130.87
SW144,9419.8158.060.82
W89,7186.0723.230.53
NW91,2696.1811.610.26
SPI<9.16322,59921.843962.902.88
9.16–11.07437,27129.601524.190.82
11.07–12.85455,11330.8134.840.16
12.85–15.35203,52113.7846.450.47
>15.3558,6483.9711.610.41
TPI<−10.7469,1204.6800.000.00
−10.74–−3.3233,43515.8000.000.00
−3.3–2.9914,95661.9462100.001.61
2.9 0 11.8201,14413.6200.000.00
>11.858,4973.9600.000.00
TWI<4.86582,46439.4300.000.00
4.86–7.27585,13039.613861.291.55
7.27–10.85237,43616.071829.031.81
>10.8572,1224.8869.681.98
TRI<3.22692,26646.8662100.001.61
3.22–7.43286,70819.4100.000.00
7.43–11.64285,45319.3200.000.00
11.64–17.59174,86411.8400.000.00
>17.5937,8602.5600.000.00
CI (100/m)<−39.691,4626.221930.654.92
−39.6–−11.3289,80719.721219.350.98
−11.3–9.01663,25545.121219.350.43
9.01–38.03341,55323.24914.520.62
>38.0383,7515.701016.132.83
Distance to stream (m)<100408,79727.672946.771.69
100–200298,96620.241625.811.28
200–300245,50416.62914.520.87
300–400157,38310.6569.680.91
>400366,50224.8123.230.13
Distance to road (m)<500224,54515.202133.872.23
500–1000193,76813.121524.191.84
1000–1500166,95711.301320.971.86
1500–2000147,4319.98914.521.45
>2000744,45150.4046.450.13
Distance to fault (m)<500102,3536.9300.000.00
500–1000103,0216.9700.000.00
1000–1500101,8306.8900.000.00
1500–2000101,9386.9000.000.00
>2000103,8467.0334.840.69
Drainage density (km/km2)<0.85398,19526.9600.000.00
0.85–1.49497,68433.6946.450.19
1.49–2.19368,97324.983454.842.20
>2.19212,30014.372438.712.69
Rainfall<248.5284,46619.265283.874.36
248.5–350.3535,31936.241016.130.45
350.3–461.02321,49121.7600.000.00
461.02–582.7208,01814.0800.000.00
>582.7127,8578.6600.000.00
Soil typeEntisols486,28232.9223.230.10
Alfisols70860.4800.000.00
Aridisols161,52610.935080.657.37
Entisols/Inceptisols479,69832.471016.130.50
Inceptisols24,7771.6800.000.00
Mollisols317,78321.5100.000.00
LU/LCAgriculture281,15919.0362100.005.25
Dense forest1550.0100.000.00
Good range86600.5900.000.00
Agri-dry farming2430.0200.000.00
Dry farming16,2351.1000.000.00
Low forest337,47122.8500.000.00
Woodland233,61915.8200.000.00
Mod-forest83,5275.6500.000.00
Mod-range105,5927.1500.000.00
Poor-range382,11425.8700.000.00
Rock23,1541.5700.000.00
Urban52200.3500.000.00
LithologyA76,4755.1800.000.00
B131,6738.9100.000.00
C114,7487.7758.061.04
D94,1496.3700.000.00
E33,7222.2800.000.00
F117,9077.9800.000.00
G31,5642.1400.000.00
H134,0599.0800.000.00
I594,53140.255791.942.28
J148,32410.0400.000.00
NDVI<−0.0118,5081.261829.0323.12
−0.01–0.0744,1452.991524.198.08
0.07–0.1289,9476.101320.973.44
0.12–0.21362,49124.591016.130.66
0.21–0.32538,55336.5369.680.26
>0.32420,47928.5200.000.00

Share and Cite

MDPI and ACS Style

Arabameri, A.; Lee, S.; Tiefenbacher, J.P.; Ngo, P.T.T. Novel Ensemble of MCDM-Artificial Intelligence Techniques for Groundwater-Potential Mapping in Arid and Semi-Arid Regions (Iran). Remote Sens. 2020, 12, 490. https://doi.org/10.3390/rs12030490

AMA Style

Arabameri A, Lee S, Tiefenbacher JP, Ngo PTT. Novel Ensemble of MCDM-Artificial Intelligence Techniques for Groundwater-Potential Mapping in Arid and Semi-Arid Regions (Iran). Remote Sensing. 2020; 12(3):490. https://doi.org/10.3390/rs12030490

Chicago/Turabian Style

Arabameri, Alireza, Saro Lee, John P. Tiefenbacher, and Phuong Thao Thi Ngo. 2020. "Novel Ensemble of MCDM-Artificial Intelligence Techniques for Groundwater-Potential Mapping in Arid and Semi-Arid Regions (Iran)" Remote Sensing 12, no. 3: 490. https://doi.org/10.3390/rs12030490

APA Style

Arabameri, A., Lee, S., Tiefenbacher, J. P., & Ngo, P. T. T. (2020). Novel Ensemble of MCDM-Artificial Intelligence Techniques for Groundwater-Potential Mapping in Arid and Semi-Arid Regions (Iran). Remote Sensing, 12(3), 490. https://doi.org/10.3390/rs12030490

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