Sensitivity analysis in groundwater vulnerability assessment based on GIS in the Mahdia-Ksour Essaf aquifer,Tunisia: a validation study

Abstract

The assessment of groundwater vulnerability to pollution has proved to be an effective tool for water resource management, especially in arid and semi-arid regions like Mahdia and Ksour Essaf. The main objective of this study is to assess the aquifer vulnerability by applying the DRASTIC method as well as using sensitivity analysis to evaluate the effect of each DRASTIC parameter on the final vulnerability map. An additional objective is to demonstrate the role of the GIS techniques in the vulnerability assessment. The DRASTIC method assigns a high vulnerability to the coast of the Mahdia-Ksour Essaf. The lowest values are observed in the southern part of the study area. A sensitivity analysis applied in this study suggests that net recharge, aquifer media and depth of groundwater are the key factors determining vulnerability. The model is validated with groundwater quality data and the results have shown strong relationships between modified DRASTIC Vulnerability Index and nitrate and chloride concentrations.

Citation Saidi, S., Bouri, S. & Ben Dhia, H. (2011) Sensitivity analysis in groundwater vulnerability assessment based on GIS in the Mahdia-Ksour Essaf aquifer, Tunisia: a validation study. Hydrol. Sci. J. 56(2), 288–304.     

Estimation of groundwater vulnerability to pollution based on DRASTIC in the Niipele sub-basin of the Cuvelai Etosha Basin,Namibia

Surface water is a scarce resource in Namibia with about sixty percent of Namibia's population dependent on groundwater for drinking purposes. With increasing population, the country faces water challenges and thus groundwater resources need to be managed properly. One important aspect of Integrated Water Resources Management is the protection of water resources, including protection of groundwater from contamination and over-exploitation. This study explores vulnerability mapping as a basic tool for protecting groundwater resources from pollution. It estimates groundwater vulnerability to pollution in the upper Niipele sub-basin of the Cuvelai-Etosha in Northern Namibia using the DRASTIC index. The DRASTIC index uses GIS to estimate groundwater vulnerability by overlaying different spatially referenced hydrogeological parameters that affect groundwater contamination. The study assesses the discontinuous perched aquifer (KDP) and the Ohangwena multi-layered aquifer 1 (KOH-1). For perched aquifers, point data was regionalized by a hydrotope approach whereas for KOH-1 aquifer, inverse distance weighting was used. The hydrotope approach categorized different parts of the hydrogeological system with similar properties into five hydrotopes. The result suggests that the discontinuous perched aquifers are more vulnerable than Ohangwena multi-layered aquifer 1. This implies that vulnerability increases with decreasing depth to water table because contaminants have short travel time to reach the aquifer when they are introduced on land surface. The nitrate concentration ranges between 2 and 288 mg/l in perched aquifers while in Ohangwena multi-layered aquifer 1, it ranges between 1 and 133 mg/l. It was observed that perched aquifers have high nitrate concentrations than Ohangwena 1 aquifer, which correlates well with the vulnerability results.     

Application of Groundwater Vulnerability Overlay and Index Methods to the Jijel Plain Area (Algeria)

Today, scientists are deeply concerned by the vulnerability of groundwater reservoirs to pollution. Relatively simple overlay and index methods can be used to produce groundwater vulnerability maps in geographic information system. In addition, this study deals with contamination from nonpoint sources. In this study, two such models, DRASTIC and GOD, were applied in the Jijel Plain area of northeast Algeria and compared with measured groundwater nitrate concentrations. This showed that results from DRASTIC were better than GOD, 69% correlation with nitrate compared to 56%. DRASTIC was better able to identify vulnerable zones along the river valleys. The DRASTIC model was then modified using the nitrate concentrations to optimize the rating score given within each parameter range and sensitivity analysis to change the weighting given for each parameter. These combined changes gave a final Pearson's correlation of 83% with nitrate. This showed that recharge, aquifer type, and topography were the key factors in controlling vulnerability to nitrate pollution.     

GIS based DRASTIC model for groundwater vulnerability assessment: Case study of the shallow mio-plio-quaternary aquifer (Southeastern Tunisia)

Groundwater is the main source of water in arid regions. Thus, groundwater pollution becomes a major issue due to the increasing contamination, which poses serious and harmful risk to the environment. Groundwater vulnerability maps can be used as a tool to help decision makers to protect groundwater resources from contamination. The vulnerability of the Mio-Plio-Quaternary shallow aquifer (Southeast Tunisia) has been assessed using a DRASTIC model based on Geographic Information System (GIS). The different parameters of the model were collected from several sources and converted into thematic maps using ArcGis©. Each DRASTIC parameter was assigned a weight and rating based on a range of information within the parameter. Groundwater vulnerability map shows a large area (48%) with high risk of pollution. It indicates that the Southern part of the aquifer and the wadi beds are the most susceptible to contamination. The measured nitrate concentration is coherent with the DRASTIC model results.     

A supervised committee machine artificial intelligent for improving DRASTIC method to assess groundwater contamination risk: a case study from Tabriz plain aquifer,Iran

Vulnerability maps are designed to show areas of greatest potential for groundwater contamination on the basis of hydrogeological conditions and human impacts. The objective of this research is (1) to assess the groundwater vulnerability using DRASTIC method and (2) to improve the DRASTIC method for evaluation of groundwater contamination risk using AI methods, such as ANN, SFL, MFL, NF and SCMAI approaches. This optimization method is illustrated using a case study. For this purpose, DRASTIC model is developed using seven parameters. For validating the contamination risk assessment, a total of 243 groundwater samples were collected from different aquifer types of the study area to analyze \( {\text{NO}}_{ 3}^{ - } \) concentration. To develop AI and CMAI models, 243 data points are divided in two sets; training and validation based on cross validation approach. The calculated vulnerability indices from the DRASTIC method are corrected by the \( {\text{NO}}_{3}^{ - } \) data used in the training step. The input data of the AI models include seven parameters of DRASTIC method. However, the output is the corrected vulnerability index using \( {\text{NO}}_{3}^{ - } \) concentration data from the study area, which is called groundwater contamination risk. In other words, there is some target value (known output) which is estimated by some formula from DRASTIC vulnerability and \( {\text{NO}}_{3}^{ - } \) concentration values. After model training, the AI models are verified by the second \( {\text{NO}}_{3}^{ - } \) concentration dataset. The results revealed that NF and SFL produced acceptable performance while ANN and MFL had poor prediction. A supervised committee machine artificial intelligent (SCMAI), which combines the results of individual AI models using a supervised artificial neural network, was developed for better prediction of vulnerability. The performance of SCMAI was also compared to those of the simple averaging and weighted averaging committee machine intelligent (CMI) methods. As a result, the SCMAI model produced reliable estimates of groundwater contamination risk.     

Correlation between nitrate contamination and ground water pollution potential

AQUIPRO, a PC-based method, was used to assess aquifer vulnerability using digital water well logs. The AQUIPRO model is a parameter/factor weighting system for rating the pollution potential of an aquifer. This method uses the well depth, as well as the clay and partial clay thickness in a well, to generate pollution potential scores. In this model, aquifer protection increases as the AQUIPRO vulnerability scores increase and ground water pollution potential decreases. Computerized water well records of 2435 domestic wells with partial chemistry data were used to determine the ground water pollution potential of Kalamazoo County, Michigan. Theoretically, low AQUIPRO pollution potential scores should have more frequent occurrences of ground water contamination events than areas with high AQUIPRO scores with similar land-use, well construction, and well densities. The relative AQUIPRO scores were compared with the frequency of occurrences of nitrate-N in ground water wells. The average nitrate-N concentrations within each relative AQUIPRO vulnerability scores category were also compared. The results indicate that domestic wells containing 5 mg/L or more nitrate-N showed a positive correlation between the frequency of occurrences of nitrate-N and relative decrease of AQUIPRO (r2 = 0.99) vulnerability scores. In other words, as the ground water pollution potential increases, the occurrence frequency of nitrate-N also increases. Furthermore, the results show that as the relative AQUIPRO (r2 = 0.96) vulnerability scores decrease, the mean nitrate-N concentrations also increase.     

An approach to catchment‐scale groundwater nitrate risk assessment from diffuse agricultural sources: a case study in the Upper Bann,Northern Ireland

DRASTIC has drawbacks in groundwater risk assessment that are important in guiding activities to prevention agricultural diffuse groundwater pollution. This paper presents an improved and GIS‐based D‐DRASTIC approach for groundwater nitrate risk assessment from diffuse agricultural sources based on DRASTIC. D‐DRASTIC considers the risk concept, nitrate loading, pollutant transport with runoff, depth to water, net recharge, aquifer media, soil media, topography, impact of the vadose zone media, and the hydraulic conductivity of the aquifer. D‐DRASTIC was developed within an ArcGIS environment and applied to the Upper Bann Catchment, Northern Ireland as a case study. D‐DRASTIC shows that ‘very high’ and ‘high’ zones of groundwater nitrate risk occupy 5% and 11% of the case study area, respectively. When considering groundwater pollution sources and pathways, the results using D‐DRASTIC are helpful in guiding the activities of groundwater pollution prevention at the catchment scale in the context of better implementation of the EU Water Framework Directive. Copyright © 2008 John Wiley & Sons, Ltd.     

Catchment-scale vulnerability assessment of groundwater pollution from diffuse sources using the DRASTIC method: a case study

Abstract

The catchment-scale groundwater vulnerability assessment that delineates zones representing different levels of groundwater susceptibility to contaminants from diffuse agricultural sources has become an important element in groundwater pollution prevention for the implementation of the EU Water Framework Directive (WFD). This paper evaluates the DRASTIC method using an ArcGIS platform for assessing groundwater vulnerability in the Upper Bann catchment, Northern Ireland. Groundwater vulnerability maps of both general pollutants and pesticides in the study area were generated by using data on the factors depth to water, net recharge, aquifer media, soil media, topography, impact of vadose zone, and hydraulic conductivity, as defined in DRASTIC. The mountain areas in the study area have “high” (in 4.5% of the study area) or “moderate” (in 25.5%) vulnerability for general pollutants due to high rainfall, net recharge and soil permeability. However, by considering the diffuse agricultural sources, the mountain areas are actually at low groundwater pollution risk. The results of overlaying the maps of land use and the groundwater vulnerability are closer to the reality. This study shows that the DRASTIC method is helpful for guiding the prevention practices of groundwater pollution at the catchment scale in the UK.

Citation Yang, Y. S. & Wang, L. (2010 Yang, Y. S. and Wang, L. 2010. A review of modelling tools for implementation of the EU Water Framework Directive in handling diffuse water pollution. Water Resour. Manage., 24: 1819–1843.  [Google Scholar]) Catchment scale vulnerability assessment of groundwater pollution from diffuse sources using the DRASTIC method: a case study. Hydrol. Sci. J. 55(7), 1206–1216.     

Vulnerability evaluation of a coastal aquifer via GALDIT model and comparison with DRASTIC index using quality parameters

In recent years, environmental assessments of groundwater resources have resulted in the development of models that help identify the vulnerable zones. An aquifer is investigated using both GALDIT and DRASTIC indices. The GALDIT model is developed to determine the vulnerability of coastal aquifers in terms of saltwater intrusion whereas the DRASTIC model is generally applicable to all aquifers. Having compared the results of both the GALDIT and DRASTIC models with quality parameters, the salinity model proved to be more appropriate in identifying the vulnerability of coastal aquifers. The results show a Pearson correlation coefficient between TDS and the GALDIT vulnerability map of 0.58 while the corresponding value for the DRASTIC index is 0.48.

EDITOR D. Koutsoyiannis

ASSOCIATE EDITOR A. Fiori     

综合指数影响评价模型在地下水防污性能评价中的应用

针对DRASTIC模型在地下水防污性能评价中存在的不足,提出模型的改进方法,在此基础上构建一套基于DRASTIC模型的综合指数影响评价（DRAICQ）模型,并应用于淮河流域中部某市地下水防污性能评价的实例研究中,绘制了研究区地下水防污性能评价分区图。结果表明：DRAICQ模型的计算结果更符合研究区的实际情况,评价效果较好;研究区地下水防污性能评价分区图对该区制定地下水防污措施具有重要的参考指导作用。     

SOM-DRASTIC: using self-organizing map for evaluating groundwater potential to pollution

Groundwater is considered as the most important water resource, especially in arid and semi-arid regions, so it is crucial to impede this source of water to be contaminated. One of the most common methods to assess groundwater vulnerability is DRASTIC method. However, the subjectivity existing in defining DRASTIC weights and ratings as well as inadaptability of the parameters involved in this method with special geology, hydrogeology, land use and climatic conditions have urged researchers to modify this method. In this paper, a new method combining a special type of the neural networks called Self-Organizing Map (SOM) and the traditional DRASTIC model resulting in the hybrid SOM-DRASTIC model is applied to modify and improve DRASTIC Model. The traditional DRASTIC method holds a summation among all negative effects of different factors contributing to vulnerability, while the proposed hybrid method is able of classifying the groundwater vulnerability and deriving the real relation existing between the DRASTIC parameters as the inputs and the vulnerability class as the output of the method. The vulnerability assessment process was performed on the Zayandeh-Rud river basin aquifers in Iran. The SOM-DRASTIC identified the northern parts of the study area as the most vulnerable areas with a drastically fractured structure, while the traditional DRASTIC ranked the western parts as the most vulnerable regions with a high rate of net recharge. The results demonstrate that the proposed method can be used by managers and decision-makers as an alternative robust tool for vulnerability-based classification and land use planning.     

A GIS-Based GOD Model and Hazard Index Analysis: The Quaternary Coastal Collo Aquifer (NE-Algeria)

The quaternary coastal Collo aquifer in northeast Algeria (NE Algeria) marks an important local water resource supporting domestic, industrial and agricultural activities. The aquifer shows signs of contamination due to the existence of various pollution sources, especially nitrogen compounds. Focusing the local identification of key vulnerable zones and related main hazard types for wise future water management, the present study highlights results from a coupled analysis of the well-established Geographical Information System (GIS)-based GOD (groundwater occurrence, overall aquifer class, depth to groundwater) hazard index analysis and the COST Action 620 plan. Most prevalent hazard types in the study area were identified as the urban/residential areas without public sewage systems, landfill and agricultural/pasturing areas. Regarding the vulnerability analysis particularly the northern aquifer region is endangered, dominated by high (22.4%) and moderate (27.4%) vulnerability classes. Central, western and southern aquifer regions are characterized by low (23.3%) and very low (26.9%) vulnerability classes. Overall, these GOD-derived results are in good agreement with earlier results obtained by the more complex DRASTIC approach. Final risk assessment and validation related to 2014/2015 nitrate sampling campaigns indicate that “high risk” and “very high risk” classes only apply to a small part of the study area in the northern sector (8%), whereas the main part (>60%) broadly affecting the central, western and southern sector only bears a low to very low risk of water pollution. Apart from a future-oriented groundwater abstraction strategy it is recommended to update the evaluation regularly to effectively consider dynamic changes of local anthropogenic activities and hazards.     

A field demonstration of the entropy-weighted fuzzy DRASTIC method for groundwater vulnerability assessment

Abstract

Groundwater vulnerability assessment based on the DRASTIC index has been widely used since the 1980s to map potential risks of groundwater contamination. However, its applicability and usefulness are affected by two uncertain and subjective factors. One is the discretization of continuous input variables and the other is the assignment of different weights to the index variables. In this study, an entropy-weighted fuzzy-optimization approach was developed to augment and improve the classic DRASTIC method by reducing the uncertainties associated with variable discretization and weight assignment. The modified DRASTIC method was applied to a study site in Shandong, north China. The entropy-weighted fuzzy-optimization approach is shown to provide a more rigorous delineation of the relative vulnerability distribution. Meanwhile, the new approach does not require the use of more parameters. The results suggest that this approach significantly improves and enhances the ability of the classic DRASTIC method in a more systematic and rigorous way.

Editor D. Koutsoyiannis

Citation Yu, C., Zhang, B.X., Yao, Y.Y., Meng, F.H., and Zheng, C.M., 2012. A field demonstration of the entropy-weighted fuzzy DRASTIC method for groundwater vulnerability assessment. Hydrological Sciences Journal, 57 (7), 1420–1432.     

Assessment of groundwater pollution risk as a characteristic of the stability of its quality

Definitions of hazard and risk of groundwater pollution are given. A deterministic method for the assessment of groundwater pollution risk using estimates of groundwater protection against and vulnerability to pollution and stability indicators of groundwater quality is considered. Also presented are the principal methodological approaches to the assessment of groundwater protection against pollution and the formation of the structure of indicators and indices characterizing the stability of groundwater quality. The structure of hazard, risks, and damages associated with groundwater pollution is shown. Expert appraisal method is used for the assessment of groundwater pollution risks.     

Assessment of aquifer vulnerability to industrial waste water using resistivity measurements. A case study, along El-Gharbyia main drain, Nile Delta, Egypt

1D resistivity sounding and 2D resistivity imaging surveys were integrated with geological and hydrochemical data to assess the aquifer vulnerability and saltwater intrusion in the north of Nile Delta, Egypt. In the present study, the El-Gharbyia main drain was considered as a case study to map the sand bodies within the upper silt and clay aquitard. Twenty Schlumberger soundings and six 2D dipole-dipole profiles were executed along one profile close to the western side of the main drain. In addition, 14 groundwater samples and 4 surface water samples from the main drain were chemically analyzed to obtain the major and trace elements concentrations.The results from the resistivity and hydrochemical data were used to assess the protection of the groundwater aquifer and the potential risk of groundwater pollution. The inverted resistivities and thicknesses of the layers above the aquifer layer were used to estimate the integrated electrical conductivity (IEC) that can be used for quantification of aquifer vulnerability. According to the aquifer vulnerability assessment of an underlying sand aquifer, the southern part of the area is characterized by high vulnerability zone with slightly fresh to brackish groundwater and resistivity values of 11-23 Ω.m below the clay cap. The resistivity sections exhibit some sand bodies within the clay cap that lead to increase the recharging of surface waste water (650 mg/l salinity) and flushing the upper part of underlying saltwater aquifer. The region in the north has saltwater with resistivity less than 6 Ω.m and local vulnerable zones within the clay cap. The inverted 2D dipole-dipole profiles in the vulnerable zones, in combination with drilling information have allowed the identification of subsoil structure around the main drain that is highly affected by waste water.     

Modification of the DRASTIC Framework for Mapping Groundwater Vulnerability Zones

The DRASTIC technique is commonly used to assess groundwater vulnerability. The main disadvantage of the DRASTIC method is the difficulty associated with identifying appropriate ratings and weight assignments for each parameter. To mitigate this issue, ratings and weights can be approximated using different methods appropriate to the conditions of the study area. In this study, different linear (i.e., Wilcoxon test and statistical approaches) and nonlinear (Genetic algorithm [GA]) modifications for calibration of the DRASTIC framework using nitrate (NO₃) concentrations were compared through the preparation of groundwater vulnerability maps of the Meshqin-Shahr plain, Iran. Twenty-two groundwater samples were collected from wells in the study area, and their respective NO₃ concentrations were used to modify the ratings and weights of the DRASTIC parameters. The areas found to have the highest vulnerability were in the eastern, central, and western regions of the plain. Results showed that the modified DRASTIC frameworks performed well, compared to the unmodified DRASTIC. When measured NO₃ concentrations were correlated with the vulnerability indices produced by each method, the unmodified DRASTIC method performed most poorly, and the Wilcoxon–GA–DRASTIC method proved optimal. Compared to the unmodified DRASTIC method with an R² of 0.22, the Wilcoxon–GA–DRASTIC obtained a maximum R² value of 0.78. Modification of DRASTIC parameter ratings was found to be more efficient than the modification of the weights in establishing an accurately calibrated DRASTIC framework. However, modification of parameter ratings and weights together increased the R² value to the highest degree.     

Predicting Pesticide Attenuation in a Fractured Aquifer Using Lumped‐Parameter Models

Neighboring springs draining fractured‐rock aquifers can display large differences in water quality and flow regime, depending on local variations of the connectivity and the aperture size distribution of the fracture network. Consequently, because homogeneous equivalent parameters cannot be assumed a priori for the entire regional aquifer, the vulnerability to pollution of such springs has to be studied on a case by case basis. In this paper, a simple lumped‐parameter model usually applied to estimate the mean transit time of water (or tracer) is presented. The original exponential piston‐flow model was modified to take land‐use distribution into account and applied to predict the evolution of atrazine concentration in a series of springs draining a fractured sandstone aquifer in Luxembourg, where despite a nationwide ban in 2005, atrazine concentrations still had not begun to decrease in 2009. This persistence could be explained by exponentially distributed residence times in the aquifer, demonstrating that in some real world cases, models based on the groundwater residence time distribution can be a powerful tool for trend reversal assessments as recommended for instance by current European Union guidelines.     

Modification of the SWAT model to simulate regional groundwater flow using a multicell aquifer

The soil and water assessment tool (SWAT) has been widely used and thoroughly tested in many places in the world. The application of the SWAT model has pointed out that 2 of the major weaknesses of SWAT are related to the nonspatial reference of the hydrologic response unit concept and to the simplified groundwater concept, which contribute to its low performance in baseflow simulation and its inability to simulate regional groundwater flow. This study modified the groundwater module of SWAT to overcome the above limitations. The modified groundwater module has 2 aquifers. The local aquifer, which is the shallow aquifer in the original SWAT, represents a local groundwater flow system. The regional aquifer, which replaces the deep aquifer of the original SWAT, represents intermediate and regional groundwater flow systems. Groundwater recharge is partitioned into local and regional aquifer recharges. The regional aquifer is represented by a multicell aquifer (MCA) model. The regional aquifer is discretized into cells using the Thiessen polygon method, where centres of the cells are locations of groundwater observation wells. Groundwater flow between cells is modelled using Darcy's law. Return flow from cell to stream is conceptualized using a non‐linear storage–discharge relationship. The SWAT model with the modified aquifer module, the so‐called SWAT‐MCA, was tested in 2 basins (Wipperau and Neetze) with porous aquifers in a lowland area in Lower Saxony, Germany. Results from the Wipperau basin show that the SWAT‐MCA model is able (a) to simulate baseflow in a lowland area (where baseflow is a dominant source of streamflow) better than the original model and (b) to simulate regional groundwater flow, shown by the simulated groundwater levels in cells, quite well.     

Using principal component analysis in the investigation of groundwater hydrochemistry of Upper Jezireh Basin,Syria

Principal component analysis (PCA) was applied to hydrochemical and isotopic data of 34 groundwater samples. This allowed the reduction of 20 variables to four significant PCs that explain 81.9% of the total variance; F1 (47.1%) explains the groundwater mineralization, whereas F2 (17%) shows isotopic enrichment and nitrate pollution. Based on an iso-factor scores map of F1, three water zones were delineated: Zone A (F1 < ?1), with fresh groundwater from the unconfined aquifer; Zone B (1 > F1 > ?1), with moderate mineralization from the confined–unconfined aquifer boundary; and Zone C (F1 > 1), with the most mineralized hot water from the confined aquifer. The iso-factor scores map of F2 delineates positive values representing samples from the unconfined aquifer, with freshwater and nitrate contamination associated with stable isotope enrichment, whereas negative values represent samples from the confined aquifer. The results clearly demonstrate the usefulness of PCA in groundwater hydrochemistry investigations.     

Characterization of Flow Dynamics and Vulnerability in a Coastal Aquifer System

Traditional aquifer vulnerability techniques primarily rely on spatial property data for a region and are limited by their ability to directly or indirectly assess flow and transport processes occurring from the surface to depth within an aquifer system. The main objective of this study was to investigate groundwater vulnerability in terms of aquifer interconnectivity and flow dynamics. A combination of stable isotopes, groundwater age‐dating (radiocarbon), and geomorphic/geogenic spatial analyses was applied to a regional, highly developed coastal aquifer to explain the presence of nitrate at depth. The average δ¹³C value (?17.3 ± 2‰ VPDB, n = 27) is characteristic of groundwater originating from locally infiltrated precipitation through extensively cultivated soils. The average δ¹⁸O and δD values (?4.0 ± 0.1‰ VSMOW, n = 27; δD: ?19.3 ± 1‰ VSMOW, n = 27, respectively) are similar to precipitation water derived from maritime sources feeding the region's surface water and groundwater. Stable and radioactive isotopes reveal significant mixing between shallow and deep aquifers due to high velocities, hydraulic connection, and input of local recharge water to depths. Groundwater overdevelopment enhances deeper and faster modern water downward flux, amplifying aquifer vulnerability. Therefore, aquifer vulnerability is a variable, dependent on the type and degree of stress conditions experienced by a groundwater system as well as the geospatial properties at the near surface.