Groundwater vulnerability assessment using an improved DRASTIC method in GIS

Groundwater management can be effectively conducted by using groundwater contamination map assessment. In this study, a modified DRASTIC approach using geographic information system (GIS) was applied to evaluate groundwater vulnerability in Kerman plain (Iran). The Wilcoxon rank-sum nonparametric statistical test was applied to modify the rates of DRASTIC. In addition, the analytic hierarchy process (AHP) method was employed to evaluate the validity of the criteria and sub criteria of all the parameters of the DRASTIC model, which proposed as an alternative treatment of the imprecision demands. The GIS offers spatial analysis in which the multi index evaluation can be effectively conducted through the AHP. The non-point source pollution was effectively determined by the modified DRASTIC method compared with the traditional method. The regression coefficient revealed the relationship between the vulnerability index and the nitrate concentration. The best result was obtained by using AHP–AHP, followed by DRASTIC–AHP, modified DRASTIC–AHP, and AHP–DRASTIC models. In this study, the DRASTIC method failed to provide satisfactory result. Additionally, by using both the original DRASTIC and the modified DRASTIC methods in the study area, AHP–AHP performed highly in the Kerman plain, suggesting that the southern and south east parts of the area considerably calls for conservation against contamination.     

Groundwater vulnerability assessment for the Banyas Catchment of the Syrian coastal area using GIS and the RISKE method

Vulnerability assessment to delineate areas that are more susceptible to contamination from anthropogenic sources has become an important element for sensible resource management and landuse planning. This contribution aims at estimating aquifer vulnerability by applying the RISKE model in Banyas Catchment Area (BCA), Tartous Prefecture, west Syria. An additional objective is to demonstrate the combined use of the RISKE model and a geographical information system (GIS) as an effective method for groundwater pollution risk assessment. The RISKE model uses five environmental parameters (Rock of aquifer media, Infiltration, Soil media, Karst, and Epikarst) to characterize the hydro-geological setting and evaluate aquifer vulnerability. The elevated eastern and low western part of the study area was dominated by high vulnerability classes, while the middle part was characterized by moderate vulnerability classes. Based on the vulnerability analysis, it was found that 2% and 39% of BCA is under low and high vulnerability to groundwater contamination, respectively, while more than 52% and 5% of the area of BCA can be designated as an area of moderate and very high vulnerability to groundwater contamination, respectively. The GIS technique has provided an efficient environment for analyses and high capabilities of handling a large amount of spatial data.     

Gaining forests but losing ground: A GIS evaluation in a Himalayan watershed

GIS overlay techniques were used to provide a quantitative historic documentation of deforestation and land-use dynamics in the Middle Mountains of Nepal between 1947 and 1990. Deforestation was most critical in the 1960s, but active afforestation programs in the 1980s have reversed the process. In spite of these trends, the degradation problem is more complex. The GIS evaluation showed that 86% of the recently afforested land is now under pine plantations located primarily at lower elevations and moderately steep slopes. In contrast, rainfed agricultural expansion is most pronounced on acidic soils and steeper, upper elevation sites, suggesting marginalization of agriculture. Agricultural expansion coupled with major losses of grazing land to pine forests are the key processes pointing towards major animal feed deficits. An alternative animal feed source is suggested through GIS using a topographically based microclimatic classification to generate a tree-planting map where the optimum ecological conditions for selective native fodder tree species are identified.     

An Aquifer Vulnerability Assessment of the Paluxy Aquifer, Central Texas, USA, Using GIS and a Modified DRASTIC Approach

D epth to water, net Recharge, Aquifer media, Soil media, Topography, Impact of the vadose zone, and Conductivity of the aquifer). Using such an approach allows one to investigate the potential for groundwater contamination on a regional, rather than site-specific, scale. Based upon data from variables such as soil permeability, depth to water, aquifer hydraulic conductivity, and topography, subjective numerical weightings have been assigned according to the variable's relative importance in regional groundwater quality. The weights for each variable comprise a GIS map layer. These map layers are combined to formulate the final groundwater pollution potential map. Using this method of investigation, the pollution potential map for the study area classifies 47% of the area as having low pollution potential, 26% as having moderate pollution potential, 22% as having high pollution potential, and 5% as having very high pollution potential.     

Denitrification in the shallow ground water of a tile-drained, agricultural watershed

Nonpoint-source pollution of surface water by N is considered a major cause of hypoxia. Because Corn Belt watersheds have been identified as major sources of N in the Mississippi River basin, the fate and transport of N from midwestern agricultural watersheds have received considerable interest. The fate and transport of N in the shallow ground water of these watersheds still needs additional research. Our purpose was to estimate denitrification in the shallow ground water of a tile-drained, Corn Belt watershed with fine-grained soils. Over a 3-yr period, N was monitored in the surface and ground water of an agricultural watershed in central Illinois. A significant amount of N was transported past the tile drains and into shallow ground water. The ground water nitrate was isotopically heavier than tile drain nitrate, which can be explained by denitrification in the subsurface. Denitrifying bacteria were found at depths to 10 m throughout the watershed. Laboratory and push-pull tests showed that a significant fraction of nitrate could be denitrified rapidly. We estimated that the N denitrified in shallow ground water was equivalent to 0.3 to 6.4% of the applied N or 9 to 27% of N exported via surface water. These estimates varied by water year and peaked in a year of normal precipitation after 2 yr of below average precipitation. Three years of monitoring data indicate that shallow ground water in watersheds with fine-grained soils may be a significant N sink compared with N exported via surface water.     

Stream transport of herbicides and metabolites in a tile-drained, agricultural watershed

The occurrence of metabolites of many commonly used herbicides in streams has not been studied extensively in tile-drained watersheds. We collected water samples throughout the Upper Embarras River watershed [92% corn, Zea mays L., and soybean, Glycine max (L.) Merr.] in east-central Illinois from March 1999 through September 2000 to study the occurrence of atrazine (2-chloro-4-ethylamino-6-isopropylamino-s-triazine), metolachlor 12-chloro-N-(2-ethyl-6-methylphenyl)-N-(methoxy-1-methylethyl) acetamide], alachlor [2-chloro-N-(2,6-diethylphenyl)-N-(methoxymethyl) acetamide], acetochlor [2-chloro-N-(ethoxymethyl)-N-(2-ethyl-6-methylphenyl) acetamide], and their metabolites. River water samples were collected from three subwatersheds of varying tile density (2.8-5.3 km tile km(-2)) and from the outlet (United States Geological Survey [USGS] gage site). Near-record-low totals for stream flow occurred during the study, and nearly all flow was from tiles. Concentrations of atrazine at the USGS gage site peaked at 15 and 17 microg L(-1) in 1999 and 2000, respectively, and metolachlor at 2.7 and 3.2 microg L(-1); this was during the first significant flow event following herbicide applications. Metabolites of the chloroacetanilide herbicides were detected more often than the parent compounds (evaluated during May to July each year, when tiles were flowing), with metolachlor ethanesulfonic acid [2-[(2-ethyl-6-methylphenyl)(2-methoxy-1-methylethyl)amino]-2-oxoethanesulfonic acid] detected most often (> 90% from all sites), and metolachlor oxanilic acid [2-[(2-ethyl-6-methylphenyl)(2-methoxy-1-methylethyl)amino]-2-oxoacetic acid] second (40-100% of samples at the four sites). When summed, the median concentration of the three chloroacetanilide parent compounds (acetochlor, alachlor, and metolachlor) at the USGS gage site was 3.4 microg L(-1), whereas it was 4.3 microg L(-1) for the six metabolites. These data confirm the importance of studying chloroacetanilide metabolites, along with parent compounds, in tile-drained watersheds.     

Evaluation of the ADAPT model for simulating nitrogen dynamics in a tile-drained agricultural watershed in central Illinois

Assessing the accuracy of agronomic and water quality simulation models in different soils, land-use systems, and environments provides a basis for using and improving these models. We evaluated the performance of the ADAPT model for simulating riverine nitrate-nitrogen (NO3-N) export from a 1500-km2 watershed in central Illinois, where approximately 85% of the land is used for maize-soybean production and tile drainage is common. Soil chemical properties, crop nitrogen (N) uptake coefficient, dry matter ratio, and a denitrification reduction coefficient were used as calibration parameters to optimize the fit between measured and simulated NO3-N load from the watershed for the 1989 to 1993 period. The applicability of the calibrated parameter values was tested by using these values for simulating the 1994 to 1997 period on the same watershed. Willmott's index of agreement ranged from 0.91 to 0.97 for daily, weekly, monthly, and annual comparisons of riverine nitrate N loads. Simulation accuracy generally decreased as the time interval decreased. Willmott's index for simulated crop yields ranged from 0.91 to 0.99; however, observed crop yields were used as input to the model. The partial N budget results suggested that 52 to 72 kg N ha(-1) yr(-1) accumulated in the soil, but simulated biological N fixation associated with soybeans was considerably greater than literature values for the region. Improvement of the N fixation algorithms and incorporation of mechanisms that describe soybean yield in response to environmental conditions appear to be needed to improve the performance of the model.     

Using the soil and water assessment tool to estimate achievable water quality targets through implementation of beneficial management practices in an agricultural watershed

Runoff from crop production in agricultural watersheds can cause widespread soil loss and degradation of surface water quality. Beneficial management practices (BMPs) for soil conservation are often implemented as remedial measures because BMPs can reduce soil erosion and improve water quality. However, the efficacy of BMPs may be unknown because it can be affected by many factors, such as farming practices, land-use, soil type, topography, and climatic conditions. As such, it is difficult to estimate the impacts of BMPs on water quality through field experiments alone. In this research, the Soil and Water Assessment Tool was used to estimate achievable performance targets of water quality indicators (sediment and soluble P loadings) after implementation of combinations of selected BMPs in the Black Brook Watershed in northwestern New Brunswick, Canada. Four commonly used BMPs (flow diversion terraces [FDTs], fertilizer reductions, tillage methods, and crop rotations), were considered individually and in different combinations. At the watershed level, the best achievable sediment loading was 1.9 t ha(-1) yr(-1) (89% reduction compared with default scenario), with a BMP combination of crop rotation, FDT, and no-till. The best achievable soluble P loading was 0.5 kg ha(-1) yr(-1) (62% reduction), with a BMP combination of crop rotation and FDT and fertilizer reduction. Targets estimated through nonpoint source water quality modeling can be used to evaluate BMP implementation initiatives and provide milestones for the rehabilitation of streams and rivers in agricultural regions.     

Nitrogen and phosphorus attenuation within the stream network of a coastal, agricultural watershed

Streams alter the concentration of nutrients they transport and thereby influence nutrient loading to estuaries downstream; however, the relationship between in-stream uptake, discharge variability, and subsequent nutrient export is poorly understood. In this study, in-stream N and P uptake were examined in the stream network draining a row-crop agricultural operation in coastal North Carolina. The effect of in-stream nutrient uptake on estuarine loading was examined using continuous measurements of watershed nutrient export. From August to December 2003, 52 and 83% of the NH4+ and PO4(3-) loads were exported during storms while concurrent storm flow volume was 34% of the total. Whole-ecosystem mass transfer velocities (Vf) of NH4+ and PO4(3-), measured using short-term additions of inorganic nutrients, ranged from 0.1 to 25 mm min(-1). Using a mass balance approach, this in-stream uptake was found to attenuate 65 to 98% of the NH4+ flux and 78 to 98% of the PO4(3-) flux in small, first-order drainage ditches. For the larger channel downstream, an empirical model based on Vf and discharge was developed to estimate the percentage of the nutrient load retained in-stream. The model predicted that all of the upstream NH4+ and PO4(3-) load was retained during base flow, while 65 and 37% of the NH4+ and PO4(3-) load was retained during storms. Remineralization from the streambed (vs. terrestrial sources) was the apparent source of NH4+ and PO4(3-) to the estuary during base flow. In-stream uptake reduced the dissolved inorganic N to dissolved inorganic P ratio of water exported to the N-limited estuary, thus limiting the potential for estuarine phytoplankton growth.     

Groundwater risk assessment at a heavily industrialised catchment and the associated impacts on a peri-urban wetland

Industrial and agricultural activities often impose significant pressures to the groundwater quality and consequently degrade wetland ecosystems that depend mostly on subsurface water flow. Groundwater vulnerability and risk mapping is a widely used approach to assess the natural protection of aquifers and the associated pollution potential from human activities. In the particular study, the relatively new Pan-European methodology (COP method) has been applied in a highly industrialized peri-urban wetland catchment, located close to Athens city, to map the intrinsic vulnerability of the aquifer and evaluate the risk potential originating from local land uses. Groundwater analysis results for various parameters, including Phenols, PCBs and nutrients, have been used to validate the vulnerability and risk estimations while a biological assessment occurred to associate the mapping results with the wetland's ecological status. The results indicated that even though the natural protection of the aquifer is relatively high due to the dominant hydrogeologic and geomorphologic conditions, the groundwater pollution risk is considerable, mainly because of the existing hazardous land uses. The water quality of the groundwater accredited these findings and the ecological status of this peri-urban wetland also indicated significant impacts from industrial effluents.     

Use of agricultural land evaluation and site assessment in Linn County,Oregon, USA

Oregon state law requires each county in the state to identify agricultural land and enact policies and regulations to protect agricultural land use. State guidelines encourage the preservation of large parcels of agricultural land and discourage partitioning of agricultural land and construction of nonfarm dwellings in agricultural areas. A land evaluation and site assessment (LESA) system was developed in Linn County to aid in the identification of agricultural land and provide assistance to decision makers concerning the relative merits of requests to partition existing parcels of ricultural land and introduce nonagricultural uses.Land evaluation was determined by calculating soil potential ratings for each agricultural soil in the county based on the soil potentials for winter wheat, annual ryegrass, permanent pasture, and irrigated sweet corn. Soil potential ratings were expressed on a scale of 0 to 150 points. The land evaluation score for a parcel consists of the weighted average soil potential rating for all of the soils in the parcel, weighted by the percentage of each soil present in the parcel.Site assessment was based on the size of a parcel and on the amount of existing conflict between agricultural and nonagricultural uses, particularly rural residential uses, both adjacent to and in the vicinity of a parcel. Parcel size refers to both size in relation to a typical field and size in relation to a typical farm unit. Conflict takes into account the number of nonfarm dwellings within 1/4 mile (0.4 km) of a parcel, the amount of the perimeter that adjoins conflicting land uses, and the residential density adjacent to the parcel. Empirical scales were derived for assigning points to each of the site assessment factors. Both parcel size and conflict were worth 75 points in the model. For parcel size, 45 points were allocated to field size and 30 points to farm-unit size. For conflict, 30 points were allocated to nonfarm dwellings within 1/4 mile and 45 points to perimeter conflicts.The LESA model was validated by testing on 23 parcels in Linn County for which requests to partition and/or convert to nonagricultural uses had been received by the County Planning Department. This testing was an essential part of the development process, as it pointed out inconsistencies and errors in the model and allowed continuous adjustment of factors and point scales. The results of application of the final model to three of the case studies are presented to illustrate the concepts.Three possible uses of the information generated by the LESA system include determining the relative agricultural value of a parcel, determining grades of agricultural land suitability, and determining the impacts of changing land use on other parcels in the vicinity. Relative agricultural value is a direct outcome of application of the evaluation criteria in the LESA model. Good, marginal, and nonagricultural grades of agricultural suitability were determined by examining the data from all 23 test cases and establishing threshold point values for soil quality, conflict, parcel size, and total LESA score. Impact analyses were not done in this study, but could be achieved by calculating LESA scores for all parcels possibly affected by a land-use change both before and after a proposed change. All three applications fall short of making a specific land-use decision, but they do provide information that should be of value to the local jurisdiction charged with making such decisions.     

Use of agricultural land evaluation and site assessment in Whitman County,Washington, USA

Whitman County, Washington, USA, one of the largest wheat producing counties in the country, developed a land evaluation and site assessment system for evaluating proposed conversions of agricultural land to heavy commercial uses. The system uses soil potential indices to determine a land evaluation score and a nine-factor site assessment system to weigh suitability for development. Details on each of these are provided along with results for four sample sites.Scientific Paper no. 7165, College of Argiculture and Home Economics Research Paper, Washington State University, Pullman, Washington, project no. 0010.     

Groundwater contamination by nitrates in the city of Konya, (Turkey): a GIS perspective

Groundwater is an essential drinking water source in the city of Konya, Turkey. Approximately 75% of the city's water consumption has been supplied from 198 groundwater wells for the last six years. Nitrate (NO(3)(-)) is one of the important water quality parameters and was measured in the water samples taken from 139 wells in 1998 and from 156 wells in 2001 within the study area of 427.5 km(2). To evaluate the nitrate data, a vector-based GIS software package ArcView GIS 3.2 was used. A hardcopy map of the city was digitized in the UTM projection system. The locations of the wells were obtained by a hand-held Global Positioning System (GPS) receiver. According to the maps produced, nitrate concentrations generally tend to increase in the city center, the average concentrations being 2.2 and 16.1mg/L for the years of 1998 and 2001, respectively. A statistical correlation procedure was also applied to well depths and nitrate concentrations. As a result, correlation coefficients of 0.259 and 0.261 were obtained for data collected in 1998 and 2001. It is concluded that the distribution of nitrate concentrations is not correlated with well depths within the study area.     

Identifying crop vulnerability to groundwater abstraction: modelling and expert knowledge in a GIS

Water use is expected to increase and climate change scenarios indicate the need for more frequent water abstraction. Abstracting groundwater may have a detrimental effect on soil moisture availability for crop growth and yields. This work presents an elegant and robust method for identifying zones of crop vulnerability to abstraction. Archive groundwater level datasets were used to generate a composite groundwater surface that was subtracted from a digital terrain model. The result was the depth from surface to groundwater and identified areas underlain by shallow groundwater. Knowledge from an expert agronomist was used to define classes of risk in terms of their depth below ground level. Combining information on the permeability of geological drift types further refined the assessment of the risk of crop growth vulnerability. The nature of the mapped output is one that is easy to communicate to the intended farming audience because of the general familiarity of mapped information. Such Geographic Information System (GIS)-based products can play a significant role in the characterisation of catchments under the EU Water Framework Directive especially in the process of public liaison that is fundamental to the setting of priorities for management change. The creation of a baseline allows the impact of future increased water abstraction rates to be modelled and the vulnerability maps are in a format that can be readily understood by the various stakeholders. This methodology can readily be extended to encompass additional data layers and for a range of groundwater vulnerability issues including water resources, ecological impacts, nitrate and phosphorus.     

Artificial Groundwater Recharge Zones Mapping Using Remote Sensing and GIS: A Case Study in Indian Punjab

Artificial groundwater recharge plays a vital role in sustainable management of groundwater resources. The present study was carried out to identify the artificial groundwater recharge zones in Bist Doab basin of Indian Punjab using remote sensing and geographical information system (GIS) for augmenting groundwater resources. The study area has been facing severe water scarcity due to intensive agriculture for the past few years. The thematic layers considered in the present study are: geomorphology (2004), geology (2004), land use/land cover (2008), drainage density, slope, soil texture (2000), aquifer transmissivity, and specific yield. Different themes and related features were assigned proper weights based on their relative contribution to groundwater recharge. Normalized weights were computed using the Saaty’s analytic hierarchy process. Thematic layers were integrated in ArcGIS for delineation of artificial groundwater recharge zones. The recharge map thus obtained was divided into four zones (poor, moderate, good, and very good) based on their influence to groundwater recharge. Results indicate that 15, 18, 37, and 30 % of the study area falls under “poor,” “moderate,” “good,” and “very good” groundwater recharge zones, respectively. The highest recharge potential area is located towards western and parts of middle region because of high infiltration rates caused due to the distribution of flood plains, alluvial plain, and agricultural land. The least effective recharge potential is in the eastern and middle parts of the study area due to low infiltration rate. The results of the study can be used to formulate an efficient groundwater management plan for sustainable utilization of limited groundwater resources.     

A GIS-based DRASTIC vulnerability and net recharge reassessment in an aquifer of a semi-arid region (Metline-Ras Jebel-Raf Raf aquifer, Northern Tunisia)

This paper aims to elaborate new generic DRASTIC aquifer vulnerability maps of the coastal aquifer of Metline-Ras Jebel-Raf Raf (Northeast of Tunisia) using the GIS technique, making the data analyses easier to handle and providing better capabilities of dealing with large spatial data. A similar study was carried out in 1999 in the same aquifer using a method based on the soil water balance equation to determine the net recharge parameter. Unfortunately, the lack of data in the study area made the results unsatisfactory. By applying the Williams and Kissel equation and the Rao relationship, we intend to demonstrate that we could correctly evaluate the net recharge parameter. Moreover, new data related to the aquifer hydraulic conductivity, the soil cover and the vadose zone lithology have become available, and allowed us to elaborate suitable DRASTIC maps.     

Water quality trends and changing agricultural practices in a midwest U.S. watershed, 1994-2006

Sediment and nutrient concentrations in surface water in agricultural regions are strongly influenced by agricultural activities. In the Corn Belt, recent changes in farm management practices are likely to affect water quality, yet there are few data on these linkages at the landscape scale. We report on trends in concentrations of N as ammonium (NH(4)) and nitrate (NO(3)), soluble reactive phosphorus (SRP), and suspended sediment (SS) in three Corn Belt streams with drainage areas of 12 to 129 km(2) for 1994 through 2006. During this period, there has been an increase in conservation tillage, a decline in fertilizer use, and consolidation of animal feeding operations in our study watersheds and throughout the Corn Belt. We use an autoregressive moving average model to include the effects of discharge and season on concentrations, LOWESS plots, and analyses of changes in the relation between discharge and concentration. We found significant declines in mean monthly concentrations of NH(4) at all three streams over the 13-yr period, declines in SRP and SS in two of the three streams, and a decline in NO(3) in one stream. When trend coefficients are converted to percent per year and weighted by drainage, area changes in concentration are -8.5% for NH(4), -5.9% for SRP, -6.8% for SS, and -0.8% for NO(3). Trends in total N and P are strongly tied to trends in NO(3), SRP, and SS and indicate that total P is declining, whereas total N is not.     

Soil phosphorus, management practices, and their relationship to phosphorus delivery in the Iowa Clear Lake agricultural watershed

Clear Lake is on Iowa's list of impaired water bodies because of high P concentration. This study assessed soil-test phosphorus (STP), management practices, and P loads from its agricultural watershed. Management practice histories and STP for eight basins were surveyed in 1999. Soil samples (15-cm depth) were analyzed for STP with agronomic [Bray P1 (BP), Olsen (OP), Mehlich 3 (M3P) and environmental [iron oxide-impregnated paper (FeP) and water extraction (WP)] tests. Total phosphorus (TP) concentrations in water discharge from five basins were measured during two years, and TP loads were measured for two basins. The agronomic P tests showed that 46 to 83% (depending on the test) of the area tested above optimum for crops. Correlations among tests were high for OP, M3P, and FeP (r > 0.96) and lower for BP and WP (r = 0.88-0.93). Moldboard- and chisel-plow tillage predominated (82% of the area). Applied P (mainly fertilizer) averaged 15 kg P ha(-1) yr(-1), and 40% of the high-testing area (M3P test) was being fertilized. The mean annual water TP concentration across five basins was 275 to 474 microg L(-1). The two-year mean TP loads for the two gauged basins were 1504 and 1510 g P ha(-1) yr(-1). Water TP concentration increased linearly with increasing STP. Relationships were stronger for M3P and FeP (R2 = 0.96-0.97 for annual means and 0.77-0.79 for storm-flow events) than for BP or WP (R2 = 0.88-0.91 and 0.59-0.69, respectively). Improving P and soil conservation practices in high-testing areas could reduce P loads to the lake.