GIS integration of hydrogeological and geoelectrical data for groundwater potential modeling in the western part of greater Kushtia district of Bangladesh

Geographical Information System (GIS) has been used in this paper to delineate groundwater resources potential in the western part of greater Kushtia district of Bangladesh, where urgent attention for augmentation of irrigation water supply is required. Thematic maps of transmissivity, hydraulic conductivity, specific yield, net recharge, aquifer thickness, surface water bodies, aquifer resistivity, overburden aquitard thickness and its resistivity have been prepared and assigned weight according to their relative importance using Analytical Hierarchical Process for the preparation of groundwater potential model. Since the values within each thematic map vary significantly, they are classified into various ranges or types and assigned ratings. Finally, the thematic maps are integrated using GIS to prepare the groundwater potential map for the study area in terms of Ground Water Potential Index (GWPI). The evolved map indicates that 22.51% of the study area have GWPI more than 0.70 and therefore, have excellent prospective for exploitation. About 69.12% of the area with GWPI ranging from 0.50 to 0.70 is also quite promising for groundwater abstraction, while the rest 8.37% area having GWPI below 0.50 indicates moderate potential. The obtained map of groundwater potential is found in good agreement with the yields of available pumping test data.     

Estimation of anisotropic hydraulic conductivity using geophysical data in a coastal aquifer of Karnataka,India

Estimation of hydraulic parameters is essential to understand the interaction between groundwater flow and seawater intrusion. Though several studies have addressed hydraulic parameter estimation, based on pumping tests as well as geophysical methods, not many studies have addressed the problem with clayey formations being present. In this study, a methodology is proposed to estimate anisotropic hydraulic conductivity and porosity values for the coastal aquifer with unconsolidated formations. For this purpose, the one-dimensional resistivity of the aquifer and the groundwater conductivity data are used to estimate porosity at discrete points. The hydraulic conductivity values are estimated by its mutual dependence with porosity and petrophysical parameters. From these estimated values, the bilinear relationship between hydraulic conductivity and aquifer resistivity is established based on the clay content of the sampled formation. The methodology is applied on a coastal aquifer along with the coastal Karnataka, India, which has significant clayey formations embedded in unconsolidated rock. The estimation of hydraulic conductivity values from the established correlations has a correlation coefficient of 0.83 with pumping test data, indicating good reliability of the methodology. The established correlations also enable the estimation of horizontal hydraulic conductivity on two-dimensional resistivity sections, which was not addressed by earlier studies. The inventive approach of using the established bilinear correlations at one-dimensional to two-dimensional resistivity sections is verified by the comparison method. The horizontal hydraulic conductivity agrees with previous findings from inverse modelling. Additionally, this study provides critical insights into the estimation of vertical hydraulic conductivity and an equation is formulated which relates vertical hydraulic conductivity with horizontal. Based on the approach presented, the anisotropic hydraulic conductivity of any type aquifer with embedded clayey formations can be estimated. The anisotropic hydraulic conductivity has the potential to be used as an important input to the groundwater models.     

Geophysical Assessment of Groundwater Potential: A Case Study from Mian Channu Area,Pakistan

An integrated study using geophysical method in combination with pumping tests and geochemical method was carried out to delineate groundwater potential zones in Mian Channu area of Pakistan. Vertical electrical soundings (VES) using Schlumberger configuration with maximum current electrode spacing (AB/2 = 200 m) were conducted at 50 stations and 10 pumping tests at borehole sites were performed in close proximity to 10 of the VES stations. The aim of this study is to establish a correlation between the hydraulic parameters obtained from geophysical method and pumping tests so that the aquifer potential can be estimated from the geoelectrical surface measurements where no pumping tests exist. The aquifer parameters, namely, transmissivity and hydraulic conductivity were estimated from Dar Zarrouyk parameters by interpreting the layer parameters such as true resistivities and thicknesses. Geoelectrical succession of five‐layer strata (i.e., topsoil, clay, clay sand, sand, and sand gravel) with sand as a dominant lithology was found in the study area. Physicochemical parameters interpreted by World Health Organization and Food and Agriculture Organization were well correlated with the aquifer parameters obtained by geoelectrical method and pumping tests. The aquifer potential zones identified by modeled resistivity, Dar Zarrouk parameters, pumped aquifer parameters, and physicochemical parameters reveal that sand and gravel sand with high values of transmissivity and hydraulic conductivity are highly promising water bearing layers in northwest of the study area. Strong correlation between estimated and pumped aquifer parameters suggest that, in case of sparse well data, geophysical technique is useful to estimate the hydraulic potential of the aquifer with varying lithology.     

Mapping groundwater recharge potential zones in arid region using GIS and Landsat approaches,southeast Tunisia

The rapid growth of population and agricultural and industrial activities has caused an increase in demand on the Jeffara aquifers of Gabes (southeast Tunisia). In fact, the over-pumping of this aquifer system has resulted in water-level declines ranging from 0.25 to 1 m/year during the past three decades. The aim of this study is to identify favourable artificial recharge sites of this aquifer system based on the combined use of remotely sensed data, a geographic information system (GIS), the Shuttle Radar Topography Mission (SRTM) product and a multiple criteria decision making (MCDM) technique. The delineation of artificial recharge zones shows high to moderate potential for groundwater recharge (40%) in the Gabes region, with high precision of good potential proposed sites. Recharge processes need to account for natural conditions and ecosystems.     

Spatial delineation of riparian groundwater within alluvium deposit of mountainous region using Laplace equation

The interaction between surface and groundwater plays a key role in a riparian ecosystem while the size of riparian groundwater has not been typically incorporated into hydrological modelling systems. An extensive geophysical survey composed of 25 individual DC electrical resistivity profiles was conducted at the Blair–Wallis site in Wyoming. The observed resistivity images show a near‐surface aquifer interpreted as the saturated alluvium deposit along the channel, rather than the geological bedrock. Based on the electrical resistivity images, it can be inferred that only the near‐surface portion of the groundwater actively interacts with the stream flow in the mountainous and hilly watershed. This study attempted the spatial extrapolation of the measured riparian aquifer depths by means of fitting functions based on the surface topography. The analysis indicated that the boundary of the riparian aquifer well corresponds to the topographical inflexion point of the hill slope profile. It was also demonstrated that the extent of alluvium deposit, where the area of riparian aquifer is indicated, can be delineated using the slope and curvature maps in the geographic information system. Then, the parabolic and biharmonic functions were tested for the groundwater depth estimation using the developed alluvium deposit map. The proposed methodology was effective if geological diffusion processes by wind and water dominated the topography. The spatial map of the active aquifer will be useful in hydrological drought analysis because it is considered to be a main source of baseflow during dry seasons.     

2D and 3D resistivity inversion of Schlumberger vertical electrical soundings in Wadi El Natrun,Egypt: A case study

The Wadi El Natrun area is characterized by a very complicated geological and hydrogeological system. 45 vertical electrical soundings (Schlumberger array) were measured in the study area to elucidate the peculiarity of this unique regime, specifically the nature of waterless area. 2D and 3D resistivity inversion based on the finite element technique and regularization method were applied on the data set. 2D and 3D model resolution was investigated through the use of the Depth and Volume of Investigation Indexes. A very good matching was found between the zones of high resistivity, the waterless area, and the non-productive wells. The low resistivity zones (corresponding to Lower Pliocene clay) were also identified. The middle resistivity fresh water aquifer zones were recognized. Available results can assist in the aquifer management by selecting the most productive zone of groundwater.     

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.     

Modelling groundwater levels in an urban coastal aquifer using artificial neural networks

The prediction of groundwater levels in a basin is of immense importance for the management of groundwater resources, especially in coastal regions where the water table fluctuations are to be limited to avoid seawater intrusion. In this paper, an Artificial Neural Network (ANN) methodology is presented to predict groundwater levels in individual wells with one month lead. Groundwater levels were also predicted in neighboring wells using model parameters from the best network of a well. This methodology is applied to an urban coastal aquifer in Andhra Pradesh state, India. The results suggest that the feed forward neural network with Levenberg Marquardt (LM) algorithm is a good choice for predicting groundwater levels in individual wells. Bayesian Regularization (BR) model parameters of Balaji Nagar well are also used successfully to predict groundwater levels in the study area. It was observed that the ANN‐based algorithms were a better choice for the prediction of groundwater levels with limited hydrological parameters. Copyright © 2007 John Wiley & Sons, Ltd.     

Geophysical techniques to aquifer locating and monitoring for industrial zones in North Hanoi, Vietnam

Geophysical methods were applied for hydrogeological targets in many countries including Vietnam. This paper presents results of using complex geophysical techniques as well as 2D electrical resistivity imaging (ERI), vertical electrical sounding (VES), very low frequency (VLF), and seismic refraction for geological structure investigation for locating the aquifers and assessing the hydrogeological conditions for groundwater potential in industrial zones of North Hanoi, Vietnam. The locations of two aquifers are determined by their depth and thickness on the basis of resistivity and seismic velocity values which were proved by stratifications of three boreholes to 40–60 m of depth on the study area. There are connections from surface water to shallow aquifer by hydraulic windows, as follows from VLF data. The deeper aquifer can be considered as a potential groundwater supply, but the water level is descending in time, as shown by hydrological monitoring. However, with careful use and by reducing sources of pollution, groundwater can continue to be an important natural resource for future.     

Geoelectrical exploration for groundwater in Al Maha Forest,Ain Jouhra,Morocco

An area of about 30 km² located in Ain Jouhra, south of Rabat, Morocco, was the subject of a geoelectric resistivity investigation. The main goal of the investigation was the assessment of the groundwater potential of the uppermost aquifer. The aquifer conditions such as depth, thickness and boundaries were also investigated. The obtained apparent resistivity curves were first analysed qualitatively and classified using simple curve shapes. Thereafter, the data were converted to resistivity and thickness pairs semi‐quantitatively by means of master curves and then quantitatively by computer modelling using ATO and Winsev software (Zohdy, 1989; Zohdy and Bisdrof, 1989). Lithological control from the available single well with a stratigraphic log aided in the correlation of the resistivity values to different rock units. Three different AB‐spacing iso‐resistivity maps, an isopach map of the main groundwater‐bearing horizon, the depth to the aquifer substratum map and five geoelectric cross‐sections were constructed. The interpretation of these soundings indicates the presence of an unconfined to semi‐confined sandy aquifer with relatively important extent and varying thickness. The maximal thickness of the aquifer is recorded in the central part of the investigated area and is thinning southwards to pinch out farther to the south. Geophysical as well as field data indicate a hydraulic connection between the upper and deeper aquifers. Indeed, the two aquifers are separated from each other by a marly substratum that is indicated throughout the area by the lowest values of the interpreted true resistivity. The value of this resistivity varies laterally, most likely due to the lateral variation in the shale‐to‐sand ratio. The altitude of the substratum decreases towards the north, and increases southwards. Regarding the availability of the groundwater in the study area, zones with high potential are theoretically expected to occur in the central part where the transversal resistance is greatest. However, sufficient water supply and high flow rates from wells intended to produce restrictively from the most upper aquifer are not likely to exist. This conclusion, which seems to be very pessimistic, is evidenced from two real field and experimental observations. The first is the rapid fall of the level of Gharnoug lake, despite the ongoing feeding by three wells. Hence, the amount of water level drop cannot be accounted for by the evaporation alone. That means that the deeper aquifer is continuously draining the upper aquifer at a high flow rate. Very low rates are recorded in all the wells that penetrated only the upper aquifer, the exception being the well that reached deeper into the lower aquifer. The flow rate in this lower aquifer measured 18 litre s^?1. Copyright © 2007 John Wiley & Sons, Ltd.     

Geophysical Investigation of Fresh‐Saline Water Interface: A Case Study from South Punjab,Pakistan

The importance of the study of fresh‐saline water incursion cannot be over‐emphasized. Borehole techniques have been widely used, but they are quite expensive, intrusive, and time consuming. The electrical resistivity method has proved very successful in groundwater assessment. This advanced technique uses the calculation of Dar‐Zarrouk (D‐Z) parameters, namely longitudinal unit conductance, transverse unit resistance, and longitudinal resistivity has been employed by using 50 vertical electrical sounding points to assess the groundwater and delineate the fresh‐saline water interface over 1045 km² area of Khanewal in Southern Punjab of Pakistan. The x‐y plots and maps of D‐Z parameters were produced to establish a decipherable vision for the occurrence and distribution of different water‐bearing formations of fresh‐saline water aquifers through a complicated situation of intermixing of different resistivity ranges for fresh‐saline water bodies. This technique is useful to reduce the ambiguity produced by the process of equivalence and suppression which cause intermixing in differentiating fresh, brackish, and saline aquifers during interpretation. The fresh‐saline water interface is correlated very well with the previous studies of water quality analysis carried out in Khanewal area. The results suggest that the D‐Z parameters are useful for demarcating different aquifer zones. The behavior and pattern of D‐Z parameters with respect to occurrence and distribution of different water‐bearing formations were effectively identified and delineated in the study area.     

Groundwater Potential Mapping Using GIS-Based Hybrid Artificial Intelligence Methods

Groundwater is one of the major valuable water resources for the use of communities, agriculture, and industries. In the present study, we have developed three novel hybrid artificial intelligence (AI) models which is a combination of modified RealAdaBoost (MRAB), bagging (BA), and rotation forest (RF) ensembles with functional tree (FT) base classifier for the groundwater potential mapping (GPM) in the basaltic terrain at DakLak province, Highland Centre, Vietnam. Based on the literature survey, these proposed hybrid AI models are new and have not been used in the GPM of an area. Geospatial techniques were used and geo-hydrological data of 130 groundwater wells and 12 topographical and geo-environmental factors were used in the model studies. One-R Attribute Evaluation feature selection method was used for the selection of relevant input parameters for the development of AI models. The performance of these models was evaluated using various statistical measures including area under the receiver operation curve (AUC). Results indicated that though all the hybrid models developed in this study enhanced the goodness-of-fit and prediction accuracy, but MRAB-FT (AUC = 0.742) model outperformed RF-FT (AUC = 0.736), BA-FT (AUC = 0.714), and single FT (AUC = 0.674) models. Therefore, the MRAB-FT model can be considered as a promising AI hybrid technique for the accurate GPM. Accurate mapping of the groundwater potential zones will help in adequately recharging the aquifer for optimum use of groundwater resources by maintaining the balance between consumption and exploitation.     

Groundwater potential zones using a combination of geospatial technology and geophysical approach: case study in Dehradun,India

ABSTRACT

A combination of geospatial, geophysical and statistical models using satellite data, the weighted index overlay (WIO) method and two-dimensional electrical resistivity tomography (2D-ERT) is applied to generate the highest potential groundwater area and to further explore the groundwater in Dehradun, India. The results show that of 19.7 km² total basin area, 0.26% falls under the “poor” category as a prospect zone for groundwater, 4.3% is “moderate”, 10.10% “moderately good”, 4.9% “good” and 0.17% “very good”. In addition, the demonstration of the geophysical survey is presented, in which Purkal Youth Society Division (PYSD) site is categorized as a shallow aquifer zone and the Guru Nanak Fifth Centenary School (GNFCS) site is a deeper aquifer zone. Our study emphasizes remote sensing and geographic information system integrated with a geophysical survey to support prospecting the most probable area and confirm the existence of groundwater.     

Characterization of seasonal local recharge using electrical resistivity tomography and magnetic resonance sounding

A groundwater recharge process of heterogeneous hard rock aquifer in the Moole Hole experimental watershed, south India, is being studied to understand the groundwater flow behaviour. Significant seasonal variations in groundwater level are observed in boreholes located at the outlet area indicating that the recharge process is probably taking place below intermittent streams. In order to localize groundwater recharge zones and to optimize implementation of boreholes, a geophysical survey was carried out during and after the 2004 monsoon across the outlet zone. Magnetic resonance soundings (MRS) have been performed to characterize the aquifer and measure groundwater level depletion. The results of MRS are consistent with the observation in boreholes, but it suffers from degraded lateral resolution. A better resolution of the regolith/bedrock interface is achieved using electrical resistivity tomography (ERT). ERT results are confirmed by resistivity logging in the boreholes. ERT surveys have been carried out twice—before and during the monsoon—across the stream area. The major feature of recharge is revealed below the stream with a decrease by 80% of the calculated resistivity. The time‐lapse ERT also shows unexpected variations at a depth of 20 m below the slopes that could have been interpreted as a consequence of a deep seasonal water flow. However, in this area time‐lapse ERT does not match with borehole data. Numerical modelling shows that in the presence of a shallow water infiltration, an inversion artefact may take place thus limiting the reliability of time‐lapse ERT. A combination of ERT with MRS provides valuable information on structure and aquifer properties respectively, giving a clue for a conceptual model of the recharge process: infiltration takes place in the conductive fractured‐fissured part of the bedrock underlying the stream and clayey material present on both sides slows down its lateral dissipation. Copyright © 2007 John Wiley & Sons, Ltd.     

Groundwater dynamics of Fongafale Islet, Funafuti Atoll, Tuvalu

Geoelectric and hydrologic surveys during spring tides revealed the spatiotemporal distribution of groundwater quality produced by tidal forcing in Fongafale Islet, Funafuti Atoll, Tuvalu. The observed low resistivity showed that saline water largely immersed the surficial Holocene aquifer, indicating that there is no thick freshwater lens in Fongafale Islet, unlike in other atoll islands of comparable size. Half of the islet was constructed by reclaiming the original swamp with porous, highly permeable coral blocks; this reclaimed area should not be considered as part of the islet width for calculation of the expected thickness of the freshwater lens. The degree of aquifer salinization depends on the topographic characteristics and the hydrologic controls on the inland propagation of the tidal forcing. Large changes in bulk resistivity and the electrical conductivity of groundwater from wells indicate that periodic salinization in phase with the semidiurnal tides was occurring widely, especially in areas at lower elevation than the high-tide level and in reclaimed areas with high permeability. Thin sheets of nearly fresh and brackish water were observed in the surficial aquifer in areas above the high-tide level and in taro swamps, respectively. The thinness of the brackish and freshwater sheets suggests that the taro swamps and the fresh groundwater resources of the islet are highly vulnerable to salinization from anticipated sea-level rise. An understanding of the inherent geologic and topographic features of an atoll is necessary to evaluate the groundwater resources of the atoll and assess the vulnerability of its water resources to climate change.     

Assimilation of historical head data to estimate spatial distributions of stream bed and aquifer hydraulic conductivity fields

Management of water resources in alluvial aquifers relies mainly on understanding interactions between hydraulically connected streams and aquifers. Numerical models that simulate this interaction often are used as decision support tools for water resource management. However, the accuracy of numerical predictions relies heavily on unknown system parameters (e.g., streambed conductivity and aquifer hydraulic conductivity), which are spatially heterogeneous and difficult to measure directly. This paper employs an ensemble smoother to invert groundwater level measurements to jointly estimate spatially varying streambed and alluvial aquifer hydraulic conductivity along a 35.6‐km segment of the South Platte River in Northeastern Colorado. The accuracy of the inversion procedure is evaluated using a synthetic experiment and historical groundwater level measurements, with the latter constituting the novelty of this study in the inversion and validation of high‐resolution fields of streambed and aquifer conductivities. Results show that the estimated streambed conductivity field and aquifer conductivity field produce an acceptable agreement between observed and simulated groundwater levels and stream flow rates. The estimated parameter fields are also used to simulate the spatially varying flow exchange between the alluvial aquifer and the stream, which exhibits high spatial variability along the river reach with a maximum average monthly aquifer gain of about 2.3 m³/day and a maximum average monthly aquifer loss of 2.8 m³/day, per unit area of streambed (m²). These results demonstrate that data assimilation inversion provides a reliable and computationally affordable tool to estimate the spatial variability of streambed and aquifer conductivities at high resolution in real‐world systems.     

Geophysical Investigation of Seepage Beneath an Earthen Dam

A hydrogeophysical survey is performed at small earthen dam that overlies a confined aquifer. The structure of the dam has not shown evidence of anomalous seepage internally or through the foundation prior to the survey. However, the surface topography is mounded in a localized zone 150 m downstream, and groundwater discharges from this zone periodically when the reservoir storage is maximum. We use self‐potential and electrical resistivity tomography surveys with seismic refraction tomography to (1) determine what underlying hydrogeologic factors, if any, have contributed to the successful long‐term operation of the dam without apparent indicators of anomalous seepage through its core and foundation; and (2) investigate the hydraulic connection between the reservoir and the seepage zone to determine whether there exists a potential for this success to be undermined. Geophysical data are informed by hydraulic and geotechnical borehole data. Seismic refraction tomography is performed to determine the geometry of the phreatic surface. The hydro‐stratigraphy is mapped with the resistivity data and groundwater flow patterns are determined with self‐potential data. A self‐potential model is constructed to represent a perpendicular profile extending out from the maximum cross‐section of the dam, and self‐potential data are inverted to recover the groundwater velocity field. The groundwater flow pattern through the aquifer is controlled by the bedrock topography and a preferential flow pathway exists beneath the dam. It corresponds to a sandy‐gravel layer connecting the reservoir to the downstream seepage zone.     

Estimation of spatial variability of aquifer parameters from geophysical methods: a case study of Sindhudurg district,Maharashtra, India

Accurate estimation of aquifer parameters, especially from crystalline hard rock area, assumes a special significance for management of groundwater resources. The aquifer parameters are usually estimated through pumping tests carried out on water wells. While it may be costly and time consuming for carrying out pumping tests at a number of sites, the application of geophysical methods in combination with hydro-geochemical information proves to be potential and cost effective to estimate aquifer parameters. Here a method to estimate aquifer parameters such as hydraulic conductivity, formation factor, porosity and transmissivity is presented by utilizing electrical conductivity values analysed via hydro-geochemical analysis of existing wells and the respective vertical electrical sounding (VES) points of Sindhudurg district, western Maharashtra, India. Further, prior to interpolating the distribution of aquifer parameters of the study area, variogram modelling was carried out using data driven techniques of kriging, automatic relevance determination based Bayesian neural networks (ARD-BNN) and adaptive neuro-fuzzy neural networks (ANFIS). In total, four variogram model fitting techniques such as spherical, exponential, ARD-BNN and ANFIS were compared. According to the obtained results, the spherical variogram model in interpolating transmissivity, ARD-BNN variogram model in interpolating porosity, exponential variogram model in interpolating aquifer thickness and ANFIS variogram model in interpolating hydraulic conductivity outperformed rest of the variogram models. Accordingly, the accurate aquifer parameters maps of the study area were produced by using the best variogram model. The present results suggest that there are relatively high value of hydraulic conductivity, porosity and transmissivity at Parule, Mogarne, Kudal, and Zarap, which would be useful to characterize the aquifer system over western Maharashtra.     

Spatial prediction of nitrate pollution in groundwaters using neural networks and GIS: an application to South Rhodope aquifer (Thrace,Greece)

Neural network techniques combined with Geographical Information Systems (GIS), are used in the spatial prediction of nitrate pollution in groundwaters. Initially, the most important parameters controlling groundwater pollution by nitrates are determined. These include hydraulic conductivity of the aquifer, depth to the aquifer, land uses, soil permeability, and fine to coarse grain ratio in the unsaturated zone. All these parameters were quantified in a GIS environment, and were standardized in a common scale. Subsequently, a neural network classification was applied, using a multi‐layer perceptron classifier with the back propagation (BP) algorithm, in order to categorize the examined area into categories of groundwater nitrate pollution potential. The methodology was applied to South Rhodope aquifer (Thrace, Greece). The calculation was based on information from 214 training sites, which correspond to monitored nitrate concentrations in groundwaters in the area. The predictive accuracy of the model developed reached 86% in the training samples, 74% in the overall sample and 71% in the test samples. This indicates that this methodology is promising to describe the spatial pattern of nitrate pollution. Copyright © 2008 John Wiley & Sons, Ltd.     

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.