Quality of stormwater runoff from an urbanised watershed

A field monitoring network was set up within the Stamford canal watershed in 1989 to study both the quantitative and qualitative aspects of storm runoff from this urbanised catchment. The data acquisition equipment comprised a continuous recording rain gauge, a water level recorder and an automatic water sampler capable of sampling storm runoff at preset intervals during rainfall events. Water samples were collected after each storm and laboratory tests were carried out on the physical and chemical properties of the storm water. Preliminary findings on the temporal variations of stormwater quality during single storms and the effects of antecedent dry weather period on the quality are presented. The average ranges of some of the significant quality parameters found in the storm runoff were also established. The quality of storm runoff from the catchment under study was found to be of an acceptable level and could potentially be developed as a water catchment area.     

Relationship between landscape characteristics and surface water quality

The effects of landscape characteristics on surface water quality were evaluated in terms of land-use condition, soil type and slope. The case area, the Chichiawan stream in the Wulin catchment in Taiwan, is Formosan landlocked salmon's natural habitat. Due to the agriculture behavior and mankind's activities, the water and environmental quality has gradually worsened. This study applied WinVAST model to predict hydrological responses and non-point source pollution (NPSP) exports in the Wulin catchment. The land-use condition and the slope of land surface in a catchment are major effect factors for watershed responses, including flows and pollutant exports. This work discussed the possible variation of watershed responses induced by the change of land-use condition, soil type and slope, etc. The results show that hydrological responses are highly relative to the value of Curve Number (CN); Pollutant exports have large relation to the average slope of the land surface in the Wulin catchment.     

Comparison of a Subjective and a Physical Approach for Identification of Priority Areas for Soil and Water Management in a Watershed – A Case Study of Nagwan Watershed in Hazaribagh District of Jharkhand, India

The present investigation was an attempt to compare the within-watershed prioritization capabilities of a physical model based SDSS with the SYI and RPI model based subjective-SDSS, conventionally devised for between-watershed prioritizations, by All India Soil Survey and Land Use Planning Division of Ministry of Agriculture, Government of India. Application of these two approaches on a test watershed situated in Damodar-Barakar catchment in India, the second most seriously eroded area in the world, demonstrated that the proposed physical model based SDSS was capable of realistically and logically mimicking the sub-watershed-scaled water and soil losses thereby suggesting its immense application potential for priority area identification within the test watershed. In contrast to the proposed physical method, the subjective approaches, which assigned totally reverse priorities to about 67–93% of the test-sub-watersheds, were observed to be incapable of realistically assessing the impact of topography and varied land use and soil types in the test watershed on their sub-watershed scaled run-off and soil loss generating potential. Besides, the physical approach could be used for assessing the annual dynamics of the total water and sediment yields, under prevailing resource management systems in the test watershed with good to moderately good correlation coefficients of 0.83 and 0.65; model efficiency coefficients of 0.54 and 0.70; mean relative errors of –4.28% and –17.97% and root mean square prediction errors of 71.8 mm and 9.63 t/ha, respectively.     

Determining storm sampling requirements for improving precision of annual load estimates of nutrients from a small forested watershed

This study sought to determine the lowest number of storm events required for adequate estimation of annual nutrient loads from a forested watershed using the regression equation between cumulative load (∑L) and cumulative stream discharge (∑Q). Hydrological surveys were conducted for 4 years, and stream water was sampled sequentially at 15-60-min intervals during 24 h in 20 events, as well as weekly in a small forested watershed. The bootstrap sampling technique was used to determine the regression (∑L-∑Q) equations of dissolved nitrogen (DN) and phosphorus (DP), particulate nitrogen (PN) and phosphorus (PP), dissolved inorganic nitrogen (DIN), and suspended solid (SS) for each dataset of ∑L and ∑Q. For dissolved nutrients (DN, DP, DIN), the coefficient of variance (CV) in 100 replicates of 4-year average annual load estimates was below 20% with datasets composed of five storm events. For particulate nutrients (PN, PP, SS), the CV exceeded 20%, even with datasets composed of more than ten storm events. The differences in the number of storm events required for precise load estimates between dissolved and particulate nutrients were attributed to the goodness of fit of the ∑L-∑Q equations. Bootstrap simulation based on flow-stratified sampling resulted in fewer storm events than the simulation based on random sampling and showed that only three storm events were required to give a CV below 20% for dissolved nutrients. These results indicate that a sampling design considering discharge levels reduces the frequency of laborious chemical analyses of water samples required throughout the year.     

Impact of storm water on groundwater quality below retention/detention basins

Groundwater from 33 monitoring of peripheral wells of Karachi, Pakistan were evaluated in terms of pre- and post-monsoon seasons to find out the impact of storm water infiltration, as storm water infiltration by retention basin receives urban runoff water from the nearby areas. This may increase the risk of groundwater contamination for heavy metals, where the soil is sandy and water table is shallow. Concentration of dissolved oxygen is significantly low in groundwater beneath detention basin during pre-monsoon season, which effected the concentration of zinc and iron. The models of trace metals shown in basin groundwater reflect the land use served by the basins, while it differed from background concentration as storm water releases high concentration of certain trace metals such as copper and cadmium. Recharge by storm water infiltration decreases the concentration and detection frequency of iron, lead, and zinc in background groundwater; however, the study does not point a considerable risk for groundwater contamination due to storm water infiltration.     

Concentrations of dissolved herbicides and pharmaceuticals in a small river in Luxembourg

Urban and agricultural areas affect the hydraulic patterns as well as the water quality of receiving drainage systems, especially of catchments smaller than 50 km(2). Urban runoff is prone to contamination due to pollutants like pesticides or pharmaceuticals. Agricultural areas are possible sources of nutrient and herbicide contamination for receiving water bodies. The pollution is derived from leaching by subsurface flow, as well as wash-off and erosion caused by surface runoff. In the Luxembourgish Mess River catchment, the pharmaceutical and pesticide concentrations are comparable with those detected by other authors in different river systems worldwide. Some investigated pesticide concentrations infringe current regulations. The maximum allowable concentration for diuron of 1.8 μg l(?-?1) is exceeded fourfold by measured 7.41 μg l(?-?1) in a flood event. The load of dissolved pesticides reaching the stream gauge is primarily determined by the amount applied to the surfaces within the catchment area. Storm water runoff from urban areas causes short-lived but high-pollutant concentrations and moderate loads, whereas moderate concentrations and high loads are representative for agricultural inputs to the drainage system. Dissolved herbicides, sulfonamides, tetracyclines, analgesics and hormones can be used as indicators to investigate runoff generation processes, including inputs from anthropogenic sources. The measurements prove that the influence of kinematic wave effects on the relationship between hydrograph and chemographs should not be neglected in smaller basins. The time lag shows that it is not possible to connect analysed substances of defined samples to the corresponding section of the hydrograph.     

An experimental method to verify soil conservation by check dams on the Loess Plateau, China

A successful experiment with a physical model requires necessary conditions of similarity. This study presents an experimental method with a semi-scale physical model. The model is used to monitor and verify soil conservation by check dams in a small watershed on the Loess Plateau of China. During experiments, the model–prototype ratio of geomorphic variables was kept constant under each rainfall event. Consequently, experimental data are available for verification of soil erosion processes in the field and for predicting soil loss in a model watershed with check dams. Thus, it can predict the amount of soil loss in a catchment. This study also mentions four criteria: similarities of watershed geometry, grain size and bare land, Froude number (Fr) for rainfall event, and soil erosion in downscaled models. The efficacy of the proposed method was confirmed using these criteria in two different downscaled model experiments. The B-Model, a large scale model, simulates watershed prototype. The two small scale models, D_a and D_b, have different erosion rates, but are the same size. These two models simulate hydraulic processes in the B-Model. Experiment results show that while soil loss in the small scale models was converted by multiplying the soil loss scale number, it was very close to that of the B-Model. Obviously, with a semi-scale physical model, experiments are available to verify and predict soil loss in a small watershed area with check dam system on the Loess Plateau, China.     

Agricultural Chemicals in the Alluvial Aquifer of a Typical County of the Arkansas Delta

Statistical methods and a Geographic Information System (GIS) were used to investigate potential indicators of ground water vulnerability to agricultural chemical contamination in a representative area of the Mississippi River alluvial aquifer. A total of 47 wells were sampled for analysis of nitrate, phosphorus, potassium, and 13 pesticides commonly-used in the area. Ten soil and hydrogeologic variables and five ground water vulnerability indices were examined to explain the variations of chemical concentrations. The results showed that no individual soil or hydrogeologic variables or their linear combinations could explain more than 25% of the variation of the chemical concentrations. A quadratic response surface model with the values of confining unit thickness, slope, soil permeability, depth to ground water, and recharge rate accounted for 62% of the variation of nitrate, 43% of P, and 83% of K, suggesting that the interactions among soil and hydrogeologic variables were significant. Observed trends of decreasing nitrate and P concentrations with increasing well depth and/or depth to ground water seemed to correlate with carbonate equilibrium in the aquifer and more reduced environment with depth. In view of uncertainties involved, it was recognized that the limitations associated with input data resolution used in GIS and the formulation of leaching indices limited their use for predicting ground water vulnerability. Misuse of pesticides could be another factor that would complicate the relationships between pesticide concentrations and the vulnerability indices.     

Occurrence and behaviour of organophosphorus insecticides in the River Wuchuan,southeast China

Organophosphorus (OP) insecticides are used worldwide, yet their persistence in the environment is not well understood. This paper summarises recent research on 17 OP insecticides in the River Wuchuan, an important catchment in China, by determining their concentrations in water, soil, sediment and plants by gas chromatography with flame photometric detection (GC-FPD). The concentrations of the total OP insecticides ranged from 92.77 to 229 ng l-1 in river water, 1.61 to 9.93 ng g-1 dry weight in soil, 1.24 to 7.56 ng g-1 dry weight in sediment and 75.28 to 326 ng g-1 dry weight in plants. There was a relatively high abundance of methamidophos, dichlorvos, malathion, omethoate and dimethoate in all water, soil, sediment and plant samples. On average, these five compounds collectively accounted for 64%, 71%, 71% and 54% of the total OP insecticide concentration in water, soil, sediment and plants, respectively, which was similar to the composition of OP insecticides in application formulations used in our study area and in China as a whole. The results therefore reflect the application pattern to some extent, and suggest that the sources of the OP insecticides are mainly from current usage. When compared with other areas, the levels of OP insecticides in the River Wuchuan were moderate. The animal and human risks from exposure to OP insecticides in water and plants were evaluated against relevant threshold values. The results therefore provide important information on the current contamination status of a key agricultural watershed in China, and point to the need for urgent action to control the use of some excessively applied and potentially persistent agrochemicals, such as methamidophos, omethoate and dimethoate.     

Transport and fate of organochlorine pesticides in the River Wuchuan,southeast China

Persistent organic pollutants (POPs) such as chlorinated pesticides are of global concern due to their widespread occurrence, persistence, bioaccumulation and toxicity to animals and human. This paper summarises recent research on 18 chlorinated pesticides in an important catchment in China, by determining their concentrations and behaviour in water, sediment, soil and plants. The concentrations of the total pesticides were in the ranges 187-893 ng l(-1) in river water, 8.53-210 ng g(-1) dry weight in soil, 2.66-13.45 ng g(-1) dry weight in river sediment, and 651-2823 ng g(-1) dry weight in plants. The predominance of beta-HCH as the major isomer of HCHs in all water, soil, sediment and plant samples was clearly observed, due to beta-HCH's resistance to biodegradation. On average beta-HCH accounted for 44%, 53%, 50%, and 46% of the total HCH concentration in water, soil, sediment and plant, respectively. Of the DDTs, DDE accounted for 48%, 43%, 53%, 55% of the total DDT, which suggested that DDT had been transformed to its metabolites, DDE and DDD, of which DDE was the more stable. The chlorinated pesticide levels in the River Wuchuan were generally below the guideline values in China, but some sites displayed levels in excess of EC Environmental Quality Standards for HCHs and DDTs. The results therefore provide important information on the current contamination status of a key agricultural watershed in China, and point to the need for urgent actions to evaluate the long-term fate and toxicity of such persistent compounds and an appropriate remediation strategy.     

Spatial distribution and transport of heavy metals in soil, ponded-surface water and grass in a pb-contaminated watershed as related to land-use practices

The aim of this study was to investigate the spatial distribution of heavy metal in soil and evaluate the dissolution of metal from soil to ponded-surface water, leaching through soil profiles and metal uptake in grass as related to different land-use practices. The data provided a scientific basis for best-management practices for land use in Khli Ti watershed. The watershed has a Pb-contamination problem from the previous operation of a Pb-ore concentrator and abandoned Zn–Pb mine. Sampling sites were selected from a land-use map, with land-use types falling into the following four categories: forest, agricultural land, residential area and road. Soil, ponded-surface water, grass samples and soil profiles were collected. The study related soil characteristics from different land-use practices and locations with observed metal concentrations in ponded-surface water and soil. High enrichment factors of Pb and As in soil were found. Partitioning coefficient, K_d values were in the order: Cr > Pb > Ni > Cu > Cd > Zn. Soil disturbance from land-use activities including tillage and traffic increased leaching of trace metal from soils. Pb in soil was significantly taken up by grass even though the Transfer Factor, TF values were rather low. Agricultural activities in the watershed must be limited. Moreover, land encroachments in the upper and middle part of the watershed which have high potential of Pb must be strictly controlled in order to reduce the Pb contamination from non-point sources.     

Simulation of soil loss processes based on rainfall runoff and the time factor of governance in the Jialing River Watershed,China

Jialing River is the largest tributary in the catchment area of Three Gorges Reservoir, and it is also one of the important areas of sediment yield in the upper reaches of the Yangtze River. In recent years, significant changes of water and sediment characteristics have taken place. The "Long Control" Project implemented since 1989 had greatly changed the surface appearance of the Jialing River Watershed (JRW), and it had made the environments of the watershed sediment yield and sediment transport change significantly. In this research, the Revised Universal Soil Loss Equation was selected and used to predict the annual average amount of soil erosion for the special water and sediment environments in the JRW after the implementation of the "Long Control" Project, and then the rainfall–runoff modulus and the time factor of governance were both considered as dynamic factors, the dynamic sediment transport model was built for soil erosion monitoring and forecasting based on the average sediment yield model. According to the dynamic model, the spatial and temporal distribution of soil erosion amount and sediment transport amount of the JRW from 1990 to 2007 was simulated using geographic information system (GIS) technology and space-grid algorithm. Simulation results showed that the average relative error of sediment transport was less than 10% except for the extreme hydrological year. The relationship between water and sediment from 1990 to 2007 showed that sediment interception effects of the soil and water conservation projects were obvious: the annual average sediment discharge reduced from 145.3 to 35 million tons, the decrement of sediment amount was about 111 million tons, and decreasing amplitude was 76%; the sediment concentration was also decreased from 2.01 to 0.578 kg/m³. These data are of great significance for the prediction and estimation of the future changing trends of sediment storage in the Three Gorges Reservoir and the particulate non-point source pollution load carried by sediment transport from watershed surface.     

Linking on-farm dairy management practices to storm-flow fecal coliform loading for California coastal watersheds

How and where to improve water quality within an agricultural watershed requires data at a spatial scale that corresponds with individual management decision units on an agricultural operation. This is particularly true in the context of water quality regulations, such as Total Maximum Daily Loads (TMDLs), that identify agriculture as one source of non-point source pollution through larger tributary watershed scale and above and below water quality investigations. We have conducted a systems approach study of 10 coastal dairies and ranches to document fecal coliform concentration and loading to surface waters at the management decision unit scale. Water quality samples were collected on a storm event basis from loading units that included: manure management systems; gutters; storm drains; pastures; and corrals and lots. In addition, in-stream samples were collected above and below the dairy facilities and from a control watershed, managed for light grazing and without a dairy facility or human residence and corresponding septic system. Samples were analyzed for fecal coliform concentration by membrane filtration. Instantaneous discharge was measured for each collected sample. Storm runoff was also calculated using the curve number method (SCS, 1985). Results for a representative dairy as well as the entire 10 dairy data set are presented. Fecal coliform concentrations demonstrate high variability both within and between loading units. Fecal coliform concentrations for pastures range from 206 to 2,288,888 cfu/100 ml and for lots from 1,933 to 166,105,000 cfu/100 ml. Mean concentrations for pastures and lots are 121,298 (SE=62,222) and 3,155,584 (SE=1,902,713) cfu/100 ml, respectively. Fecal coliform load from units of concentrated animals and manure are significantly more than units such as pastures while storm flow amounts were significantly less. Compared with results from earlier tributary scale studies in the watershed, this systems approach has generatedwater quality data that is beneficial for management decisions because of its scale and representation of current management activities. These results are facilitating on-farm changes through the cooperative efforts of dairy managers, regulatory agency staff, and sources of technical and financial assistance.     

A water quality monitoring network design methodology for the selection of critical sampling points: Part II

In order to resolve the spatial component of the design of a water quality monitoring network, a methodology has been developed to identify the critical sampling locations within a watershed. This methodology, called Critical Sampling Points (CSP), focuses on the contaminant total phosphorus (TP), and is applicable to small, predominantly agricultural-forested watersheds. The CSP methodology was translated into a model, called Water Quality Monitoring Station Analysis (WQMSA). It incorporates a geographic information system (GIS) for spatial analysis and data manipulation purposes, a hydrologic/water quality simulation model for estimating TP loads, and an artificial intelligence technology for improved input data representation. The model input data include a number of hydrologic, topographic, soils, vegetative, and land use factors. The model also includes an economic and logistics component. The validity of the CSP methodology was tested on a small experimental Pennsylvanian watershed, for which TP data from a number of single storm events were available for various sampling points within the watershed. A comparison of the ratios of observed to predicted TP loads between sampling points revealed that the model's results were promising.     

Groundwater vulnerability assessment using DRASTIC and Pesticide DRASTIC models in intense agriculture area of the Gangetic plains,India

Delineating areas susceptible to contamination from anthropogenic sources form an important component of sustainable management of groundwater resources. The present research aims at estimating vulnerability of groundwater by application of DRASTIC and Pesticide DRASTIC models in the southern part of the Gangetic plains in the state of Bihar. The DRASTIC and Pesticide DRASTIC models have considered seven parameters viz. depth to water level, net recharge, aquifer material, soil material, topography, impact of vadose zone and hydraulic conductivity. A third model, Pesticide DRASTIC LU has been adopted by adding land use as an additional parameter, to assess its impact on vulnerability zonation. The DRASTIC model indicated two vulnerable categories, moderate and high, while the Pesticide DRASTIC model revealed moderate, high and very high vulnerable categories. Out of the parameters used, depth to water level affected the vulnerability most. The parameter caused least impact was topography in DRASTIC, while in case of Pesticide DRASTIC and Pesticide DRASTIC LU models, the parameter was hydraulic conductivity. A linear regression between groundwater NO₃ concentrations and the vulnerability zonation revealed better correlation for Pesticide DRASTIC model, emphasising the effectiveness of the model in assessing groundwater vulnerability in the study region. Considering all three models, the most vulnerable areas were found to be concentrated mainly in two zones, (i) in the south-western part along Ekangarsarai-Islampur patch and (ii) around Biharsharif-Nagarnausa area in the central part. Both zones were characterised by intensive vegetable cultivation with urban areas in between.     

Atmospheric deposition and storm induced runoff of heavy metals from different impermeable urban surfaces

Contaminants deposited on impermeable surfaces migrate to stormwater following rainfall events, but accurately quantifying their spatial and temporal yields useful for mitigation purposes is challenging. To overcome limitations in current sampling methods, a system was developed for rapid quantification of contaminant build-up and wash-off dynamics from different impervious surfaces. Thin boards constructed of concrete and two types of asphalt were deployed at different locations of a large carpark to capture spatially distributed contaminants from dry atmospheric deposition over specified periods of time. Following experimental exposure time, the boards were then placed under a rainfall simulator in the laboratory to generate contaminant runoff under controlled conditions. Single parameter effects including surface roughness and material composition, number of antecedent dry days, rain intensity, and water quality on contaminant build-up and wash-off yields could be investigated. The method was applied to quantify spatial differences in deposition rates of contaminants (TSS, zinc, copper and lead) at two locations varying in their distance to vehicle traffic. Results showed that boards exposed at an unused part of the carpark >50 m from vehicular traffic captured similar amounts of contaminants compared with boards that were exposed directly adjacent to the access route, indicating substantial atmospheric contaminant transport. Furthermore, differences in contaminant accumulation as a function of surface composition were observed. Runoff from asphalt boards yielded higher zinc loads compared with concrete surfaces, whereas runoff from concrete surfaces resulted in higher TSS concentrations attributed to its smoother surfaces. The application of this method enables relationships between individual contaminant behaviour and specific catchment characteristics to be investigated and provides a technique to derive site-specific build-up and wash-off functions required for modelling contaminant loads from impermeable surfaces.     

Land degradation assessment by geo-spatially modeling different soil erodibility equations in a semi-arid catchment

Land degradation by soil erosion is one of the most serious problems and environmental issues in many ecosystems of arid and semi-arid regions. Especially, the disturbed areas have greater soil detachability and transportability capacity. Evaluation of land degradation in terms of soil erodibility, by using geostatistical modeling, is vital to protect and reclaim susceptible areas. Soil erodibility, described as the ability of soils to resist erosion, can be measured either directly under natural or simulated rainfall conditions, or indirectly estimated by empirical regression models. This study compares three empirical equations used to determine the soil erodibility factor of revised universal soil loss equation prediction technology based on their geospatial performances in the semi-arid catchment of the Saraykoy II Irrigation Dam located in Cankiri, Turkey. A total of 311 geo-referenced soil samples were collected with irregular intervals from the top soil layer (0-10?cm). Geostatistical analysis was performed with the point values of each equation to determine its spatial pattern. Results showed that equations that used soil organic matter in combination with the soil particle size better agreed with the variations in land use and topography of the catchment than the one using only the particle size distribution. It is recommended that the equations which dynamically integrate soil intrinsic properties with land use, topography, and its influences on the local microclimates, could be successfully used to geospatially determine sites highly susceptible to water erosion, and therefore, to select the agricultural and bio-engineering control measures needed.     

Use of cluster and principal component analyses to profile areas in California where ground water has been contaminated by pesticides

An empirical approach to profiling areas of ground water contamination by pesticides was devised that did not rely upon determining the level of vulnerability between land areas and that did not assume any particular pathway for ground water contamination. Climatic and soil data were obtained for 1-square mile sections of land in California where pesticide residues had been found in well water samples and the detection was attributed to legal agricultural applications. These sections were designated as known contaminated (KC) sections. Climate and soil data were also obtained for sections which lacked either well sampling data or a positive pesticide detection. These sections were designated as candidate sections. Statistical procedures were used to cluster groups of KC sections first with respect to climate characteristics and then with respect to soil characteristics. Principal components analysis (PCA) was used to construct a statistical profile of soil variables for each cluster of KC sections. A method based on the PCA was developed to compare the similarity of soil profiles derived for each KC section cluster to individual candidate sections. Since the profiling scheme was based only on data from KC sections, candidate sections that did not match any KC cluster profile could only be considered dissimilar to contaminated sections, receiving a status of not-classified. This profiling method is flexible and it can be revised to incorporate updated well sampling information.     

红壤丘陵区典型流域地表水水环境特征研究——以浏阳河流域为例

选取浏阳河流域为例,根据近十年河流断面监测数据,采用改进的综合污染指数法来评价水质污染程度,研究了该红壤丘陵区典型河流水质的时空变化特征,并结合土地利用和土壤特征等分析地表水质变化原因。结果表明,改进的综合污染指数法有较好的适用性;从时间特征上看,由于面源污染加剧,使得浏阳河近十年的水质污染呈增长趋势;从空间特征上看,浏阳河从上游到下游,河流污染呈增长趋势,上游水质较好,中游表现为重金属铅和汞的污染较大,而下游则是氨氮污染加剧。     

An analytical model for the assessment of pesticide exposure levels in soils and groundwater

The movement and degradation of pesticide residues in soils and groundwater are complex processes affected by soil physical, (bio)chemical, and hydrogeological properties, climatic conditions, and agricultural practices. This work presents a physically-based analytical model suitable for long-term predictions of pesticide concentrations in groundwater. The primary interest is to investigate the impact of soil environment, related physical and (bio)chemical processes, especially, volatilization, crop uptake, and agricultural practices on long-term vulnerability of groundwater to contamination by pesticides. The soil is separated into root and intermediate vadose zones, each with uniform properties. Transport in each soil zone is modeled on the basis of complete mixing, by spatial averaging the related point multiphase-transport partial differential equation (i.e., linear-reservoir models). Transport in the aquifer, however, is modeled by a two-dimensional advection-dispersion transport equation, considering adsorption and first-order decay rate. Vaporization in the soil is accounted for by assuming liquid-vapor phase partitioning using Henry's law, and vapor flux (volatilization) from the soil surface is modeled by diffusion through an air boundary layer. Sorption of liquid-phase solutes by crops is described by a linear relationship which is valid for first-order (passive) crop uptake. The model is applied to five pesticides (atrazine, bromacil, chlordane, heptachlor, and lindane), and the potential for pesticide contamination of groundwater is investigated for sandy and clayey soils. Simulation results show that groundwater contamination can be substantially reduced for clayey soil environments, where bio(chemical) degradation and volatilization are most efficient as natural loss pathways for the pesticides. Also, uptake by cross can be a significant mechanism for attenuating exposure levels in ground-water especially in a sandy soil environment, and for relatively persisting pesticides. Further, simulations indicate that changing agricultural practices can have a profound effect on vulnerability of groundwater to mobile and relatively persisting pesticides.