Combining rough sets and GIS techniques to assess aquifer vulnerability characteristics in the semi-arid South Texas

The coastal semi-arid region of South Texas is undergoing significant growth causing an enormous burden on its limited water resources. Understanding regional-scale vulnerability of this resource is important for sustainable water resources management and land use development. In this study, DRASTIC methodology is integrated with an information-analytic technique called rough sets to understand groundwater vulnerability characteristics in 18 different counties of South Texas. The rough set theory provides three useful metrics: the strength factor which depicts how vulnerability characteristics occur over the area; the certainty factor computes the relative probabilities for various vulnerability states within a county and the coverage factor which elucidates the fraction of a specific vulnerability state present in each county. The coupling of rough sets with GIS is particularly advantageous to cluster counties exhibiting similar vulnerability characteristics and to obtain other related insights. The application of the approach indicates that the groundwater vulnerability exhibits greater variability along the coast than in the interior sections of the area. The shallow aquifer in Aransas, DeWitt, Goliad and Gonzales counties is the most vulnerable, while the aquifer in Duval, Jim Wells, Karnes, Live Oak, Nueces and San Patricio is less vulnerable. This approach should prove useful to regional planners and environmental managers entrusted with the protection of groundwater resources.     

Determining shallow aquifer vulnerability by the DRASTIC model and hydrochemistry in granitic terrain,southern India

Shallow aquifer vulnerability has been assessed using GIS-based DRASTIC model by incorporating the major geological and hydrogeological factors that affect and control the groundwater contamination in a granitic terrain. It provides a relative indication of aquifer vulnerability to the contamination. Further, it has been cross-verified with hydrochemical signatures such as total dissolved solids (TDS), \(\hbox {Cl}^{-},\, \hbox {HCO}_{3}^{-},\, \hbox {SO}_{4}^{2-}\) and \(\hbox {Cl}^{-}/\hbox {HCO}_{3}^{-}\) molar ratios. The results show four zones of aquifer vulnerability (i.e., negligible, low, moderate and high) based on the variation of DRASTIC Vulnerability Index (DVI) between 39 and 132. About 57% area in the central part is found moderately and highly contaminated due to the 80 functional tannery disposals and is more prone to groundwater aquifer vulnerability. The high range values of TDS (2304–39,100 mg/l); \(\hbox {Na}^{+}\)(239– 6,046 mg/l) and \(\hbox {Cl}^{-}\) (532–13,652 mg/l) are well correlated with the observed high vulnerable zones. The values of \(\hbox {Cl}^{-}/\hbox {HCO}_{3}^{-}\) (molar ratios: 1.4–106.8) in the high vulnerable zone obviously indicate deterioration of the aquifer due to contamination. Further cumulative probability distributions of these parameters indicate several threshold values which are able to demarcate the diverse vulnerability zones in granitic terrain.     

An aquifer vulnerability assessment of the Benin Formation aquifer,Calabar, southeastern Nigeria,using DRASTIC and GIS approach

An aquifer vulnerability of the Benin Formation aquifer (Calabar, southern Nigeria) has been assessed using a combination of DRASTIC index and GIS technology. The assessment was necessitated by the fact that uncontrolled disposal of domestic, industrial and agricultural wastes have caused groundwater contamination. Therefore, prevention of contamination, monitoring and management of the aquifer was urgently required to increase the efficient use of the current water supplies. The DRASTIC method uses seven parameters (depth to groundwater table, net recharge, aquifer media, soil media, topography, influence of vadose zone and hydraulic conductivity), which were used to produce vulnerability maps. The drastic vulnerability index ranged between 124 and 170. The vulnerability map shows that the aquifer is highly vulnerable in southeastern parts of the area covering about 22 %. The medium vulnerability area covers about 56.8 % of Calabar extending from the southwest to northern parts. 21.2 % of the area covering the central and northern parts the area lies within the low vulnerability zone. The present industrial and activities are located in the eastern and western parts, which falls within the low-medium vulnerability areas. Documented nitrate concentration in hand-dug wells and boreholes are in agreement with vulnerability zones. Sensitivity analysis was performed to evaluate the sensitivity of each parameter between map layers such that subjectivity can be reduced to an extent and new weights computed for each DRASTIC parameter. As management options sensitive areas, especially in the southern parts of Calabar area, should be protected from future development.     

GIS-based DRASTIC, Pesticide DRASTIC and the Susceptibility Index (SI): comparative study for evaluation of pollution potential in the Nabeul-Hammamet shallow aquifer, Tunisia

Three vulnerability index models were applied to assess the pollution potential of Nabeul-Hammamet shallow aquifer, Tunisia: DRASTIC, Pesticide DRASTIC and the Susceptibility Index (SI). An output map layer of each one was obtained using a geographic information system (GIS). The SI layer was overlain with DRASTIC and Pesticide DRASTIC and the percentage areas of agreement and divergence in vulnerability categories were extracted. DRASTIC results suggest the aquifer has mostly low vulnerability. Pesticide DRASTIC and SI identify three vulnerability categories (low, moderate, high) in the aquifer. Published data on current chemical groundwater composition indicate that parts of the aquifer are highly contaminated, revealing that DRASTIC underestimates the risk of pollution; Pesticide DRASTIC and SI reflect this risk better. Agreement in vulnerability categories between the two last models is found for 64 % of the aquifer area. To help manage land-use allocation and prevent Nabeul-Hammamet-aquifer contamination, DRASTIC is not recommended. Pesticide DRASTIC and SI are recommended but for slightly different applications. SI helps in the monitoring of current vulnerable areas and, thus, in contamination prevention. Pesticide DRASTIC could better intervene as a criterion in a multi-criteria analysis to select the best sites for specific on-the-ground practice or future land use.     

Comparison of DRASTIC and DC resistivity modeling for assessing aquifer vulnerability in the central Nile Delta,Egypt

DRASTIC indexing and integrated electrical conductivity (IEC) modeling are approaches for assessing aquifer vulnerability to surface pollution. DRASTIC indexing is more common, but IEC modeling is faster and more cost-effective because it requires less data and fewer processing steps. This study aimed to compare DRASTIC indexing with IEC modeling to determine whether the latter is sufficient on its own. Both approaches are utilized to determine zones vulnerable to groundwater pollution in the Nile Delta. Hence, assessing the nature and degree of risk are important for realizing effective measures toward damage minimization. For DRASTIC indexing, hydrogeological factors such as depth to aquifer, recharge rate, aquifer media, soil permeability, topography, impact of the vadose zone, and hydraulic conductivity were combined in a geographical information system environment for assessing the aquifer vulnerability. For IEC modeling, DC resistivity data were collected from 36 surface sounding points to cover the entire area and used to estimate the IEC index. Additionally, the vulnerable zones identified by both approaches were tested using a local-scale resistivity survey in the form of 1D and 2D resistivity imaging to determine the permeable pathways in the vadose zone. A correlation of 0.82 was obtained between the DRASTIC indexing and IEC modeling results. For additional benefit, the obtained DRASTIC and IEC models were used together to develop a vulnerability map. This map showed a very high vulnerability zone, a high-vulnerability zone, and moderate- and low-vulnerability zones constituting 19.89, 41, 27, and 12%, respectively, of the study area. Identifying where groundwater is more vulnerable to pollution enables more effective protection and management of groundwater resources in vulnerable areas.     

A comparative study of four vulnerability mapping methods in a detritic aquifer under mediterranean climatic conditions

This paper presents the results of a comparative study relating to the application of four vulnerability mapping methods, GOD, AVI, DRASTIC and SINTACS, in a pilot detritic aquifer situated in NW Morocco, known as the Martil–Alila aquifer. The principal objective of this work is to determine the most suitable such methods for this aquifer type within a Mediterranean context, and to show the effect of the rainfall variations that are characteristic of the Mediterranean climate on the degree of vulnerability. The methods applied distinguish five classes of vulnerability, these being irregularly divided up in space, with the division varying according to the method in question. The vulnerability maps obtained by the different methods strongly suggest that the eastern half of the aquifer is more vulnerable to contamination than the western half, for all hydrological situations. The effect of climatic conditions on the degree of vulnerability is well represented by the DRASTIC, according to which the aquifer is moderately to strongly vulnerable during humid hydrological years and weakly to moderately vulnerable during dry ones. For the other methods, this climatic effect is limited to the area occupied by the two predominant classes (“High” and “Low” for GOD and “High” and “Moderate” for SINTACS) while it is null for AVI. In conclusion, DRASTIC appears the most suitable for mapping the vulnerability to contamination of Mediterranean coastal detritic aquifers such as the Martil–Alila aquifer.     

An extension to the DRASTIC model to assess groundwater vulnerability to pollution: application to the Haouz aquifer of Marrakech (Morocco)

An extension to the DRASTIC model is proposed in order to assess aquifer vulnerability to pollution. In contrast to the DRASTIC model, which considers the unsaturated and saturated zones together and computes a global intrinsic vulnerability index, the suggested approach discriminates between the aquifer vertical vulnerability (a concept related to the pollutant percolation) and the groundwater susceptibility (a concept that depends on the behaviour and uses of the groundwater). This approach is applied to the Haouz aquifer (Morocco) that supplies water to the Marrakech area. This aquifer is widely overexploited and there is evidence that the groundwater quality is threatened by various sources of pollution. Evaluation of the vertical vulnerability indicates that the aquifer mainly presents a moderate-to-weak vertical vulnerability. The zones potentially most favourable to pollutant percolation are mainly located in Central Haouz, along or near the surface wadis. The aquifer susceptibility is high in places located near the N’Fis, Baaja and Issil wadis. Everywhere else, low-to-moderate susceptibility is observed. This new approach therefore enables areas of vertical vulnerability and areas of susceptibility to be delineated separately. As a result, it constitutes a valuable decision-making tool for optimising the management of aquifer water resources and land-use planning.     

Supervised committee fuzzy logic model to assess groundwater intrinsic vulnerability in multiple aquifer systems

Groundwater vulnerability modeling is an alternative approach to evaluate groundwater contamination especially in areas affected by intensive anthropogenic activities. However, the DRASTIC model as a well-known method to assess groundwater vulnerability suffers from the inherent uncertainty associated with its seven essential parameters. In this study, three different fuzzy logic (FL) models (Sugeno fuzzy logic, Mamdani fuzzy logic, and Larsen fuzzy logic) are adopted to improve the DRASTIC system to be more realistic. The vulnerability map of groundwater from multiple aquifer systems (i.e., karstic, alluvium, and complex) in Basara basin, Iraq, was created using the FL models. Validation of the FL models results using NO₃-N concentration obtained from wells and springs of the study area indicating that all of the three FL models are applicable for improving the DRASTIC model. However, each of the FL models has its own advantages for groundwater vulnerability estimation in different types of aquifer systems in the Basara basin. Therefore, this study proposes the supervised committee fuzzy logic (SCFL) as a multimodel method to combine the advantages of individual FL models. The SCFL method confirms that no water well with high NO₃-N levels would be classified as low risk and vice versa. The study suggests that this approach has provided a convenient estimation of pollution risk in the study area and therefore, a more accurate prediction of the intrinsic vulnerability to pollution in the multiple aquifer system can be achieved through SCFL method.     

Using Generic and Pesticide DRASTIC GIS-based models for vulnerability assessment of the Quaternary aquifer at Sohag, Egypt

Groundwater resources in the Sohag area, Egypt are currently threatened by contamination from municipal and industrial activities, and agricultural pesticides. To cope with the growing population, there has been development in the desert zone on both sides of the Nile Valley including agricultural investment areas, wastewater disposal sites, new urban areas, and industry. Use of agrochemicals in the old cultivated and newly reclaimed lands and wastewater disposal sites in the study area represent the most hazardous contamination sources. Prevention of contamination and management of the Quaternary aquifer is urgently needed. To address vulnerability assessment of the Quaternary aquifer, the Generic and Pesticide DRASTIC GIS-based models have been used. The Generic DRASTIC index ranged between 94 and 189, and the Pesticide DRASTIC index ranged between 94 and 226. The results showed that 83% of the Quaternary aquifer is characterized by the high and very high vulnerability classes to municipal, industrial and agricultural pesticides contamination. It was found that nearly all the development projects are located in the very high vulnerability class areas. Management alternatives for the Quaternary aquifer may be improved by application of these models, allowing sensitive groundwater sources to be protected for continuing use in the future.     

Modelling the vulnerability of groundwater to contamination in an unconfined alluvial aquifer in Pakistan

The area of Thal Doab is located in the Indus Basin and is underlain by a thick alluvial aquifer called the Thal Doab aquifer (TDA). The TDA is undergone intense hydrological stress owing to rapid population growth and excessive groundwater use for livestock and irrigated agricultural land uses. The potential impact of these land uses on groundwater quality was assessed using a DRASTIC model in a Geographic Information System environment. Seven DRASTIC thematic maps were developed at fixed scale and then combined into a groundwater vulnerability map. The resultant vulnerability index values were grouped into four zones as low, moderate, high and very high. The study has established that 76% of the land area that is underlain by the TDA has a high to very high vulnerability to groundwater contamination mainly because of a thin soil profile, a shallow water table and the presence of soils and sediments with high hydraulic conductivity values. In addition, only 2 and 22% of the total area lie in low and moderate vulnerability zones, respectively. The outcomes of this study can be used to improve the sustainability of the groundwater resource through proper land-use management.     

Vulnerability evaluation of a coastal plain sand aquifer with a case example from Calabar,southeastern Nigeria

A point count index method using a well drillers log and field measurements has been developed following the DRASTIC and SINTACS procedures to map and evaluate the vulnerability of a coastal plain aquifer to surface and near surface contamination. The input parameters with the acronym CALOD include clay layer thickness (C), aquifer media character (A), lateritic layer thickness (L), overlying layer character (O) and the depth to groundwater level (D). The CALOD vulnerability potential index (CALOD index) is computed as the sum of the products of weights and ratings assigned to each of the input parameters. The CALOD index, varying between 15 and 75, is divided into four classes: high (>60), high-medium, (40–60), low-medium (20–40) and low (<20). The CALOD index is then used to produce a vulnerability potential map for the area. From the map, areas of high, high-medium and medium-low are consistent with the upper gravelly aquiferous zone while areas of medium-low and low are restricted to the deeper lower sandy aquiferous layer. The most important parameters affecting groundwater vulnerability to pollution in coastal areas include saturated thickness of the aquifer, depth to groundwater level, lateritic layer thickness and the aquifer media character. The concentration of some chemical pollution indicators (electrical conductivity, K, NO₃, Cl and metal load) are relatively higher for the highly vulnerable shallower upper gravelly unit in comparison to the less vulnerable deeper sandy unit. This method is very suitable for coastal plain sand aquifers especially, where data is scare.     

Groundwater vulnerability assessment in shallow aquifer of Kathmandu Valley using GIS-based DRASTIC model

In this paper, groundwater aquifer vulnerability map has been developed by incorporating the major geological and hydro-geological factors that affect and control the groundwater contamination using GIS based DRASTIC model. This work demonstrates the potential of GIS to derive a map by overlying various spatially referenced digital data layers that portrays cumulative aquifer sensitivity ratings across the Kathmandu Valley, Nepal, providing a relative indication of groundwater vulnerability to contamination. In fact, the groundwater is the major natural resources in Kathmandu for drinking purpose. The decline in groundwater levels due to the over exploitation and thus extracted water from shallow aquifer has been contaminated by the infiltration of pollutants from polluted river and land surface is continuous and serious. As the demand for water for human and industrial use has escalated and at the same time, the engineering and environmental costs are much higher for new water supplies than maintaining the existing sources already in use. Management of groundwater source and protecting its quality is therefore essential to increase efficient use of existing water supplies. Aquifer vulnerability maps developed in this study are valuable tools for environmental planning and predictive groundwater management. Further, a sensitivity analysis has been performed to evaluate the influence of single parameters on aquifer vulnerability assessment such that some subjectivity can be reduced to some extent and then new weights have been computed for each DRASTIC parameters.     

Assessment of aquifer vulnerability based on GIS and DRASTIC methods: a case study of the Senirkent-Uluborlu Basin (Isparta, Turkey)

Aquifer vulnerability has been assessed in the Senirkent-Uluborlu Basin within the Egirdir Lake catchment (Turkey) using the DRASTIC method, based on a geographic information system (GIS). There is widespread agriculture in the basin, and fertilizer (nitrate) and pesticide applications have caused groundwater contamination as a result of leaching. According to hydrogeological data from the study area, surface water and groundwater flow are towards Egirdir Lake. Hence, aquifer vulnerability in the basin should be determined by water quality in Egirdir Lake. DRASTIC layers were prepared using data such as rainfall, groundwater level, aquifer type, and hydraulic conductivity. These data were obtained from hydrogeological investigations and literature. A regional-scale aquifer-vulnerability map of the basin was prepared using overlay analysis with the aid of GIS. A DRASTIC vulnerability map, verified by nitrate in groundwater data, shows that the defined areas are compatible with land-use data. It is concluded that 20.8% of the basin area is highly vulnerable and urgent pollution-preventions measures should be taken for every kind of relevant activity within the whole basin.     

Comparative approach of three popular intrinsic vulnerability methods: case of the Beni Amir groundwater (Morocco)

The study of the intrinsic vulnerability of groundwater resources to pollution is an effective tool to control their quality degradation and contribute to their protection. It is used to delimit the vulnerable zones which do not withstand a large flow of pollutants introduced from the soil surface. Three methods of assessing the intrinsic vulnerability of groundwater: DRASTIC (Depth to water table, Recharge, Aquifer, Soil type, Topography, Impact of zone vadose, Hydraulic conductivity), DRSTI, and GOD (Groundwater occurrence, Overall aquifer class and Depth of water table) coupled with a geographic information system (GIS) are applied to the groundwater of Beni Amir, and they are compared in order to adopt the method which better characterizes the vulnerability of the aquifer to pollution. The validation of this application was made by measurements of the nitrate levels in the aquifer. Because the pollution of groundwater, in this plain, is a direct consequence of agricultural activities characterized by an intensive fertilizer application. The results clearly show that the rate of the coincidence, between the measured nitrate concentrations and the different classes of vulnerability of three methods, is 81.81, 54.54, 72.72, and 27.27%, respectively, for methods DRASTIC (classification of Engel et al. (Int Assoc Hydrol Sci Publi 235:521–526, 1996)), DRASTIC (classifications of Aller et al. (1987)), DRSTI, and GOD. Of this rate of coincidences, the DRASTIC method, with the classification of Engel et al. (Int Assoc Hydrol Sci Publi 235:521–526, 1996), allows a finer assessment and turns out the most representative of the study area.     

Evaluation of groundwater vulnerability in El-Bahariya Oasis,Western Desert,Egypt, using modelling and GIS techniques: A case study

The Nubian Sandstone Aquifer (NSSA) is the main groundwater resource of the El-Bahariya Oasis, which is located in the middle of the Western Desert of Egypt. This aquifer is composed mainly of continental clastic sediments of sandstone with shale and clay intercalations of saturated thickness ranging between 100 and 1500 m. Vulnerability assessment to delineate areas that are more susceptible to contamination from anthropogenic sources has become an important element for sustainable resources management and land use planning. Accordingly, this research aims to estimate the vulnerability of NSSA by applying the DRASTIC model as well as utilising sensitivity analyses to evaluate the relative importance of the model parameters for aquifer vulnerability in the study area. The main objective is to demonstrate the combined use of the DRASTIC and the GIS techniques as an effective method for groundwater pollution risk assessment, and mapping the areas that are prone to deterioration of groundwater quality and quantity. Based on DRASTIC index (DI) values, a groundwater vulnerability map was produced using the GIS. The aquifer analysis in the study area highlighted the following key points: the northeastern and western parts of the NSSA were dominated by ‘High’ vulnerability classes while the northwestern and southeastern parts were characterised by ‘Medium’ vulnerability classes. The elevated central part of the study area displayed ‘Low’ aquifer vulnerability. The vulnerability map shows a relatively greater risk imposed on the northeastern part of the NSSA due to the larger pollution potential of intensive vegetable cultivation. Depth-to-water, topography and hydraulic conductivity parameters were found to be more effective in assessing aquifer vulnerability.     

Vulnerability mapping of shallow groundwater aquifer using SINTACS model in the Jordan Valley area, Jordan

Jordan Valley is one of the important areas in Jordan that involves dense agricultural activities, which depend on groundwater resources. The groundwater is exploited from an unconfined shallow aquifer which is mainly composed of alluvial deposits. In the vicinity of the Kafrein and South Shunah, the shallow aquifer shows signs of contamination from a wide variety of non-point sources. In this study, a vulnerability map was created as a tool to determine areas where groundwater is most vulnerable to contamination. One of the most widely used groundwater vulnerability mapping methods is SINTACS, which is a point count system model for the assessment of groundwater pollution hazards. SINTACS model is an adaptation for Mediterranean conditions of the well-known DRASTIC model. The model takes into account several environmental factors: these include topography, hydrology, geology, hydrogeology, and pedology. Spatial knowledge of all these factors and their mutual relationships is needed in order to properly model aquifer vulnerability using this model. Geographic information system was used to express each of SINTACS parameters as a spatial thematic layer with a specific weight and score. The final SINTACS thematic layer (intrinsic vulnerability index) was produced by taking the summation of each score parameter multiplied by its specific weight. The resultant SINTACS vulnerability map of the study area indicates that the highest potential sites for contamination are along the area between Er Ramah and Kafrein area. To the north of the study area there is a small, circular area which shows fairly high potential. Elsewhere, very low to low SINTACS index values are observed, indicating areas of low vulnerability potential.     

Groundwater vulnerability and risk mapping of the Angad transboundary aquifer using DRASTIC index method in GIS environment

The protection and preservation of groundwater resources are compulsory, particularly in the arid and semiarid countries where the waters are scarce. The effects of increasing urbanization, economic development, and agricultural activities, along with the erratic and scarce rainfall, contribute to the quantitative and qualitative deterioration of these resources. This paper attempts to produce groundwater vulnerability and risk maps for the Angad transboundary aquifer using DRASTIC model. The data which correspond to the seven parameters of the model were collected and converted to thematic maps in Geographic Information System environment. The modified DRASTIC map, which is the summation of the DRASTIC index and the network fractures maps, shows two degreed of vulnerability: medium and high. This map is then integrated with a land use map to assess the potential risk of groundwater to pollution in the Angad transboundary aquifer. There are three risk zones that are identified: moderate, high, and very high.     

基于地下水流数值模型的改进DRASTIC方法

地下水脆弱性评价作为地下水资源保护和地下水开发利用规划的一个重要工具,被广泛的应用于实际工作中。尝试利用地下水数值模型为改进的DRASTIC方法提供数据支持,并以北京市平原区为例探讨地下水脆弱性评价方法。评价结果与传统方法在高值区和低值区具有很好的对应性,而基于模型的方法在地下水水位计算、含水层介质和水力传导系数确定上较传统方法更具优势,如地下水位的计算上较传统方法更为客观地体现含水介质对地下水运动的影响,且能够方便地获得模拟期内任意时间的流场数据;经由模型调试后的含水层参数数据,较传统方法更为准确。评价结果分区之间的变化较传统方法更为平滑,更符合水文地质条件渐变的特性。     

Groundwater pollution risk mapping for the Eocene aquifer of the Oum Er-Rabia basin, Morocco

Sustainable development requires the management and preservation of water resources indispensable for all human activities. When groundwater constitutes the main water resource, vulnerability maps therefore are an important tool for identifying zones of high pollution risk and taking preventive measures in potential pollution sites. The vulnerability assessment for the Eocene aquifer in the Moroccan basin of Oum Er-Rabia is based on the DRASTIC method that uses seven parameters summarizing climatic, geological, and hydrogeological conditions controlling the seepage of pollutant substances to groundwater. Vulnerability maps were produced by using GIS techniques and applying the “generic” and “agricultural” models according to the DRASTIC charter. Resulting maps revealed that the aquifer is highly vulnerable in the western part of the basin and areas being under high contamination risk are more extensive when the “agricultural” model was applied.     

A GIS-based DRASTIC model for assessing groundwater vulnerability in the Ordos Plateau, China

Groundwater plays a key role in arid regions as the majority of water is supplied by it. Groundwater pollution is a major issue, because it is susceptible to contamination from land use and other anthropogenic impacts. A study was carried out to build a vulnerability map for the Ordos Plateau using the DRASTIC model in a GIS environment. The map was designed to show the areas of the highest potential for groundwater pollution based on hydrogeological conditions. Seven environmental parameters, such as depth to water table, net recharge, aquifer media, soil media, topography, impact of the vadose zone media, and hydraulic conductivity of the aquifer, were incorporated into the DRASTIC model and GIS was used to create a groundwater vulnerability map by overlaying the available data. The results of this study show that 24.8 % of the study area has high pollution potential, 24.2 % has moderate pollution potential, 19.7 % has low pollution potential, and the remaining 31.3 % of the area has no risk of groundwater pollution. The regional distribution of nitrate is well correlated with the DRASTIC vulnerability index. In contrast to this, although the DRASTIC model indicated that the western part had no risk, nitrate concentrations were higher in some of these areas. In particular, higher nitrate concentrations were recorded along river valleys and around lakes, such as the Mulin River valley. This is mainly caused by the intensive agricultural development and favorable conditions for recharge along river valleys.