A Simplified Method for the Assessment of Groundwater Vulnerability to Contamination ()
1. Introduction
The maps of Groundwater vulnerability to contamination are important tools used to assist the landuse planners in addressing the problems that groundwater might have as a result of changing the landuse in a certain area. It is also used to predict the pollutants movement in the soil. This will allow planners to modify the potential harm to groundwater before serious impacts occur [1] .
There are several methods used to assess the groundwater vulnerability inclu- ding process-based mathematical models, statistical methods and overlay and in- dexing.
The process-based mathematical models are analytical or numerical solutions of mathematical equations to predict the contaminant transport in both space and time, which distinguish them from the other methods. These methods are constrained by the lack of data and the computational difficulties [2] .
The statistical models are usually used in evaluating, determining, and quanti- fying the association between measures of vulnerability and various types of information that could be related to vulnerability. The use of statistical methods is limited by the requirement for high quality data, cost and time constrains [2] .
The overlay and index methods are used for combining maps of parameters considered to be important in contaminant, where each attribute is assigned a numerical score based on its perceived importance [2] [3] . The overlay and index methods rely on simple mathematical representations of expert opinion rather than process representation or empirical data. Their major negative aspect is the fact that assigning numerical values to the descriptive entities and relative weights for the different attributes is subjective.
Examples of the developed overlay and index methods include the following (Table 1):
1. DRASTIC ( [4] - [11] ).
2. GOD ( [12] [13] [14] ).
3. EPIK ( [15] [16] [17] ).
4. SINTACS ( [18] [19] [20] ).
5. PI ( [21] ).
6. COP ( [11] , [17] and [22] [23] ).
It appears from Table 1 that these methods require subsurface data including
Table 1. Methods and used parameters for groundwater vulnerability assessment.
the unsaturated zone (GOD, DRASTIC, SINTACS, EPIK, PI and COP), hydraulic conductivity (DRASTIC and SINTACS) and aquifer thickness (SINTACS). Also, these methods require data about net recharge (DRASTIC, SINTACS, EPIK, PI and COP). These data requirement might restrict researchers and lead them either to estimate or exclude one or more of these parameters. [9] excluded the hydrau- lic conductivity and estimated the net recharge based on an equation derived by [8] .
This research is an attempt to derive a Simplified Groundwater Vulnerability Index (SGVI) based on experts’ opinions using available and inexpensive data to assist researchers in assessing groundwater vulnerability as a part of their Environmental Impact Assessment (EIA).
2. Methodology
2.1. Index Development
The adopted methodology for developing a simplified groundwater vulnerability index (SGVI) was based on the Analytical Hierarchy Process (AHP). AHP is ba- sed upon the construction of a series of Pair-Wise Comparison Matrices (PCMs), which compare each criterion to one another. The scale of PWC is divided from 1 to 9 for PCM elements to estimates the relationship between the criteria on the scale (9, 7, 5, 3, 1) where the value of 1 suggests that the criteria are equally important and a value of 9 leads one to infer that the criterion under consideration is extremely important in relation to the other criterion with which the compa- rison is made [24] (Table 2).
The processes of pairwise comparison matrix that are computed to achieve and check the pairwise comparison matrix is acceptable or not are illustrated in Figure 1 ( [25] [26] [27] ).
Based on [25] [26] [27] , the following Equations were used to calculate λmax, CI, and CR
(1)
λ represents the sum of consistency vectors divided by number of the consistency vector.
Where A is the known judgment matrix and n is the order of the matrix B.
Table 2. Scale of relative importance [24] .
Figure 1. The processes of pairwise comparison matrix.
(2)
CI represents the consistency index based on the observations of λ.
Where λmax is the average value of the consistency vector, and n is the number of criteria.
The consistency Ratio (RI) is calculated using Equation (3):
(3)
Table 3 lists the consistency index of a randomly generated pairwise compa- rison matrix [27] [28] .
In this research, nine experts from various Jordanian universities and organisations in the fields of groundwater, surface water, soil, landuse and GIS were interviewed to suggest a Simplified Groundwater Vulnerability Index (SGVI) that requires inexpensive data. The experts agreed that the following parameters should be included in the groundwater vulnerability index:
Depth of Groundwater (DG): Depth determines the depth which a contaminant must travel before reaching the aquifer.
Soil Texture (ST): Soil has a significant impact on the amount of recharge water which infiltrates into the groundwater.
Lineament Density (LD): Higher values of lineament density might pose more threat to groundwater by allowing contamination to reach the water table.
Rainfall (RF): Rainfall is the major carrier of contamination to groundwater through infiltration.
Topographic Slope (TS): Slope controls the likelihood that a pollutant will run off or remain on the surface long enough to infiltrate.
Drainage Density (DD): Drainage density is a useful measure for runoff potential. High drainage density means higher potential for runoff and hence less infiltration.
Land Use/Land Cover (LL): Landuse/Land cover affect the potential for groundwater recharge and hence impact the groundwater vulnerability to contamination.
The experts were also asked to put a weight to each parameter based on the supplied simple questionnaire listed in Table 4. The pair-wise comparison
matrix of the experts opinions was then calculated (Table 5). Table 6 lists the computed values of weights, λmax, CI, RI and CR for experts opinions. Experts were
Table 3. Average random consistency indices (RI) for different number of criteria.
Table 4. A sample of the questionnaire used to determine the relative importance of index parameters.
also asked to suggest the ratings for each parameter based a scale 1 (Lowest) to 4 (Highest). Table 7 summarises the ratings and weights for the Simplified Groundwater Vulnerability Index parameters. The final classification of the SGVI is listed in Table 8. Figure 2 illustrates the methodology adopted to calculate the final weights for the SGVI parameters based on the experts’ opinions.
2.2. Investigated Area
In order to test the SGVI, an area in the Northern part of Jordan was selected as shown in Figure 3.
The investigated has an area of 3200 km2 which comprises 3.59% of the total area of Jordan. The climate within the investigated area is characterised by hot dry summers and wet cold winters. Rainfall varies between 50 mm in the South East to 150 mm in the North West (Figure 4).
Elevation in the area varies between 530 m above sea level in the South and South East to 1224 m above sea level in the North near the Syrian border (Figure 5).
Surface water flows from the North towards the South and South West. Most of the Drainages is coming from Syria and flow towards the Azraq Oasis to the South of the investigated area (Figure 6).
The soil texture within the investigated area is classified into 4 classes: Loam, Sandy Loam, Silt Loam and Silty Clay Loam as shown in Figure 7. The shallow groundwater aquifer lies at a depth that varies between 160 m in the South to
more than 300 m in the North. This variation is due to the high variation in ground topography that characterise the area.
Table 5. The pair-wise comparison matrix of the experts opinions.
Table 6. The computed values of weights (priority vector), CI, RI and CR for experts opinions.
Table 7. Ratings and weights for the SGVI.
Table 8. Simplified Groundwater Vulnerability Index Classification.
Figure 2. The adopted methodology for weight calculation.
Figure 3. The investigated area within Jordan.
2.3. Data Collection
The data needed for this research were collected at no cost from various resources within Jordan and from internet resources. Table 9 lists the data used in this research and their sources.
3. Data Analysis and Results
The following flowchart (Figure 8) shows the methodology adopted within GIS environment to calculate the SGVI. This Figure illustrates the GIS data and tools used to map the final SGVI.
The soil texture map shown in Figure 7 was subjected to GIS manipulation by adding the appropriate ratings according to Table 7. The soil map was then con- verted into raster format and multiplied by the parameter weight (Figure 9). The ST values varied between 0.18 to 0.72.
The density function in ArcGIS was used to calculate the lineament density for the investigated area, which then classified according to Table 7. The outcome of this operation was then multiplied by the parameter weight (Figure 10). LD varied between 0.162 to 0.648.
The rainfall isohyets map shown in Figure 4 was subjected to GIS manipulation by adding the appropriate ratings according to Table 7. The rainfall isohyets map was then converted into raster format and multiplied by the parameter weight (Figure 11). RF values varied between 0.136 to 0.408.
The slope map for the investigated area was extracted from the ASTER DEM (30 m) for the investigated area (Figure 5). The slope was classified according to Table 7 and multiplied by the parameter weight (Figure 12). TS values varied between 0.099 to 0.396.
ASTER DEM was used to generate the drainage map for the investigated area (Figure 6) using the Flow Direction and Flow Accumulation tools in ArcGIS. The density function in ArcGIS was then used to calculate the drainage density for the
Table 9. Used data and their sources.
Figure 8. Analysis methodology to calculate the final SGVI.
investigated area, which then classified according to Table 7. The outcome of this operation was then multiplied by the parameter weight (Figure 13). DS values va- ried between 0.09 in to 0.36.
Landsat 8 imagery (30 m) of August 2016 was used to generate the landu- se/land cover map for the investigated area using unsupervised classification technique. The resulted map was subjected to GIS manipulation by adding the appropriate ratings according to Table 7. The landuse/land cover map was then converted into raster format and multiplied by its weight (Figure 14). It appears that LL values varied between 0.09 in most of the investigated area to 0.36 in sporadic patches in the North and North West of the investigated area. Most of the investigated area where given 0.09 due the fact that the dominant land cover in the area is basaltic rocks.
The final SGVI map (Figure 15) is the result of adding all parameters in GIS and then classifying the output based on Table 8. The summary of the SGVI for the investigated area (Table 10) showed that the area with moderate-low vulnerability has an area of 1929.7 km2 (60.3% of the total area), while the low vulnerability area has an area of 1161.42 km2 (36.3% of total area). The remaining parts of the investigated area are 102.3 km2 (3.2% of total area) as a moderate-high vulnerability and 6.58 km2 (0.2% of total area) has no risk.
4. Sensitivity Analysis of the SGVI
4.1. Statistical Analysis of the SGVI Parameters
The statistical analysis for the SGVI parameters (Table 11) indicated that ST and LL parameters contribute the highest risk of contamination with high mean values of 0.384 and 0.358 respectively. The parameters DG, LD, RF, TS and DD
contribute moderate risk of contamination with mean values of 0.244, 0.298, 0.269, 0.277 and 0.277 respectively.
The ST parameter is highly variable as its coefficient of variation (CV) value is
Figure 15. Simplified Groundwater Vulnerability Index (SGVI).
Table 10. A summary of the SGVI outcome for the investigated area.
Table 11. The statistical summary of the SGVI parameters.
50%. The DG parameter has no variation within the investigated area as it has only a single value for the entire area (0.244). LD, TS, RF and DD are moderately variable with CV values 37.92% and 32%, 21.19% and 19.61% respectively. The LL co- ver parameter is the least variable parameter (CV = 6.98%). The sensitivity is calculated based on the rating and weights designated to the feature classes of each parameter.
4.2. Map Removal Analysis
The statistics on the removal of one statistically significant parameter on the SGVI values (Table 12) showed that the most important parameters to contamination were DG, DD, RF and TS followed by LD, LL, and ST. The highest values were associated with DG (1.84), DD (1.829), RF (1.815) and TS (1.807). The parameter (ST) showed the lowest sensitivity value (1.701).
4.3. Map Removal Sensitivity Analysis
Map removal sensitivity analysis is a test developed by [29] . This test was used in this research to identify the sensitivity of each parameter in the SGVI map. [29] developed an equation (Equation (4)) to calculate the sensitivity index (S) for the parameters used in the vulnerability index.
(4)
With: S is the sensitivity index of the parameter.
V is the intrinsic vulnerability index of the method; N is the total number of parameters used to calculate V; Vxi represents the intrinsic vulnerability index obtained after removal of the parameter X.
Based on the determined partial indexes (Table 12) and Equation (4), the sensitivity index was calculated for each parameter of the SGVI (Table 13). Table 13 indicates that all parameters except the DG parameter have a strong influence on the SGVI map. DG with respective sensitivity index of 0.009 has the lowest impact on the SGVI. This is due to the fact that the entire investigated area was given a fixed number of 0.244 to represent this parameter (Depth to groundwater > 150 m).
5. Conclusions and Recommendations
Based on experts’ opinions, a new simplified groundwater vulnerability index was developed to investigate groundwater vulnerability to contamination. Seven parameters including depth to groundwater, soil texture, lineament density, rainfall, topographic slope, drainage density, and landuse/land cover were suggested by the experts. The common factor between these parameters is inexpensive data sources to build each one of them. This is in contradiction to the parame- ters found in other groundwater vulnerability indexes such as GOD, DRASTIC, SINTACS, EPIK, PI and COP.
Based on the outcomes of this research and the statistical tests conducted on the SGVI parameters, it is established that:
Soil texture and landuse/land cover contribute the highest risk of contamina-
Table 12. The partial index calculated by removing one parameter of the SGVI.
Table 13. Sensitivity index according to the map removal sensitivity analysis test for SGVI vulnerability map.
tion to groundwater while depth to groundwater, lineament density, rainfall, topographic slope and drainage density contribute moderate risk of contami- nation.
Depth to groundwater, drainage density, rainfall and topographic slope have the highest importance in the index.
All parameters have strong impacts on the index except the depth to ground- water.
Based on that, it is concluded that the new simplified groundwater vulnerability index (SGVI) is applicable to investigate groundwater vulnerability to conta- mination in areas where subsurface data is limited or not available.
Based on this conclusion, it is recommended to test the SGVI in combination with other indexes such DRASTIC, SINTACS or EPIK especially in areas where data for these indexes are available to compare its outcomes with their outcomes. It is also recommended to look for other parameters that might contribute to groundwater vulnerability to contamination under the condition which these parameters are affordable and inexpensive.
Acknowledgements
The authors of this research would like to acknowledge and appreciate the contribution of the Jordanian experts who assisted the research team in the development of the SGVI index.