Maximum Contaminant Concentration Research Articles

Analytical model for organic contaminant transport in a cut-off wall and aquifer dual-domain system considering barrier arrangements

This paper presents an analytical study of organic contaminants transport in a cut-off wall and aquifer dual-domain system, considering the effects of the inlet boundary conditions and cut-off structural arrangements. The comprehensive sensitivity analysis of parameters focusing on the breakthrough time, attenuation time and cumulative concentration are presented using the Monte Carlo simulation and Sobol global method. The simplified constant inlet boundary condition can lead to an excessively conservative prediction of the contaminant breakthrough compared to the ‘finite mass’ and ‘decaying source’ boundary conditions. The cut-off wall hydraulic performance can be enhanced by reducing the contaminant's head loss, shape factor, half-life and cut-off wall hydraulic conductivity while increasing the retardation factor. The contaminant's half-life can largely influence the maximum contaminant concentration, attenuation time and breakthrough time. For example, the maximum contaminant concentrations for T1/2 = 1.4 years and T1/2 = 100 years are 13 and 123 times greater than that for T1/2 = 0.1 year, respectively. The influence of the variation of shape factor on the breakthrough curve should be taken into consideration. Altering the structural arrangement of the cut-off wall to accommodate a smaller shape factor increases the contaminant breakthrough time. The proposed solution allows the analysis of a cut-off wall and aquifer system with different inlet boundary conditions and structural arrangements of the cut-off wall.

Bioaccessible metals in dust materials from non-sulfide Zn deposit and related hydrometallurgical operation

Mining and processing of ores in arid (desert) areas generates high amounts of dust, which might be enriched in potentially harmful elements. We studied dust fractions of ores, soils, and technological materials from mining and related hydrometallurgical operation at former Skorpion Zinc non-sulfide Zn deposit in southern Namibia (closed and placed under maintenance in 2020). Chemical and mineralogical investigation was combined with oral bioaccessibility testing of fine dust fractions (<48 μm and <10 μm) in simulated gastric fluid (SGF) to assess potential risk of intake of metallic contaminants (Cd, Cu, Pb, Zn) for staff operating in the area. The bulk metals concentrations were largely variable and ranked as follows: soils < tailings ≪ Skorpion ores < imported ores and dross used for feed ore blending. Maximum contaminant concentrations in the original granular materials were 927 mg Cd/kg, 9150 mg Cu/kg, 50 g Pb/kg and 706 g Zn/kg, respectively, and generally increased as a function of decreasing grain size. The highest bioaccessible concentrations of Cd and Pb yielded imported ores from Taiwan and Turkey and, together with the milled dross, these samples also exhibited the highest Zn bioaccessibilities. The exposure estimates calculated for a worker (weighing 70 kg) in this mining/ore processing operation at a dust ingestion rate of 100 mg/day indicated that most dust samples (soils, tailings, Skorpion ores) exhibited metals intake values far below tolerable daily intake limits. The overall health risk was limited in all mining and ore processing areas except for the ore blending area, where imported ores and recycled dross enriched in bioaccessible Cd, Pb and/or Zn were used for the ore blending. Safety measures required by the mine operator (wearing of masks by the operating staff) helped to prevent the staff's exposure to potentially contaminated dust even in this blending ore area.

The Contaminate Infiltration Model for Heat and Concentration Transport Within Porous Media Under Electromagnetic Fields

Abstract This paper studies the groundwater model of the influence of physical parameters, including input frequency of the electromagnetic, and input concentration of contaminants in groundwater, on the velocity pattern, temperature distribution, and concentration distribution of convective heat transfer in saturated porous media as soil. The mathematical models have solved seven equations in this simulation study, i.e., Maxwell's equation, heat transfer in fluid and solid phases, momentum, and concentration equations. The effect of frequencies and input concentrations of contaminants on the convective heat transfer and concentration distribution in porous media as soil under an electromagnetic wave is investigated. The results indicate that the electromagnetic wave frequency of 2.45 GHz has the most influence on the temperature distribution, velocity patterns, and concentration distribution of the fluid within the porous media as soil during saturated flow in groundwater. The inlet fluid concentration of the contaminant at 30 mol/dm3 has the most impact on the temperature distribution between the implementation of an electromagnetic wave of 2.45 GHz. So, this numerical model provides simple decision data based on comparing the maximum contaminant concentrations of porous media as soil samples with surface soil screening levels such as petroleum engineering and agricultural engineering. This result can be used by the engineer as a guide to determine whether further investigation is needed.

Numerical Analysis of Potential Contaminant Migration from Abandoned In Situ Coal Conversion Reactors

Abstract. In the context of a potential utilisation of coal resources located in the Mecsek mountain area in Southern Hungary, an assessment of groundwater pollution resulting from a potential water-borne contaminant pool remaining in in situ coal conversion reactors after site abandonment has been undertaken in the scope of the present study. The respective contaminants may be of organic and inorganic nature. A sensitivity analysis was carried out by means of numerical simulations of fluid flow as well as contaminant and heat transport including retardation to assess spatial contaminant migration. Hereby, the main uncertainties, e.g., changes in hydraulic gradient and hydraulic contributions of the complex regional and local fault systems in the study area, were assessed in a deterministic way to identify the relevant parameters. Overall 512 simulations of potential groundwater contamination scenarios within a time horizon exceeding the local post-operational monitoring period were performed, based on maximum contaminant concentrations, cumulative mass balances as well as migration distances of the contaminant plume. The simulation results show that regional faults represent the main contaminant migration pathway, and that contamination is unlikely assuming the given reference model parametrisation. However, contamination within a simulation time of 50 years is possible for specific geological conditions, e.g., if the hydraulic conductivity of the regional faults exceeds a maximum value of 1 × 10−5 m s−1. Further, the parameter data analysis shows that freshwater aquifer contamination is highly non-linear and has a bimodal distribution. The bivariate correlation coefficient heatmap shows slightly positive correlations for the pressure difference, the fault permeability and the simulation time, as well as a negative correlation for the retardation coefficient. The results of this sensitivity analysis have been integrated into a specific toolkit for risk assessment for that purpose.

Determination of patulin in apple juice by liquid chromatography-electrospray tandem mass spectrometry

Patulin is a mycotoxin produced by several fungi, (Penicillium, Aspergillus, Byssochlamys). The main sources of patulin intake in human diet are apples, apple juice and apple nectar, and for this reason, apple based foods are monitored for the presence of this mycotoxin. Commission Regulation EC No 1881/2006 lays down maximum residue limits (MRLs) of 50 µg/kg in apple juice and cider, 25 µg/kg in solid apple products, and 10 µg/kg in products for infants and young children. In Serbia, maximum permitted amounts of patulin in fruit juices, reconstituted concentrated fruit juices and fruit nectars, as well as in solid apple products, including apple compote and apple puree, intended for direct human consumption are prescribed in the Regulation on maximum concentrations of certain contaminants in foodstuffs. This paper presents the LC-MS/MS method for quantitative determination of patulin in apple juice. Criteria for method validation were taken from Commission Decision 2002/657/EC. Linearity was confirmed in the concentration ranges of 0-100 µg/kg, with the limit of detection (LoD) of 9.85 µg/kg. The performance of the method was successfully verified by participating in a proficiency study.

On the Maximum Concentration of Contaminants in Natural Aquifers

On the Maximum Concentration of Contaminants in Natural Aquifers

A practical approach for the determination of environmental quality standards-based discharge limits: the case of Tersakan sub-basin of Yeşilırmak River in Turkey.

The control of point source discharges to rivers has become more significant since the establishment of environmental quality standards (EQSs). Many countries, including Turkey, have set EQS values for various contaminants. One important challenge regarding these EQSs is to reconcile the effluent limits that are technically and economically achievable with the ones that are required to accomplish the EQSs. The Tersakan sub-basin of Yeşilırmak River acquires good examples of this challenge due to the industrial and agricultural discharge activities present. In this study, a new, simplistic, and less data-driven approach is developed to facilitate this compromise and implemented for all suitable discharge points within the sub-basin. The foundation of this approach is that effluent discharges may mix and become diluted within negligibly short distances from the point of discharge where exceedance of EQSs can be permissible. The approach modularly combines different analytical solutions of the advective-dispersive mass transport equation that are applicable under different mixing conditions and estimates maximum allowable discharge concentrations of contaminants. The results of the case study in the Tersakan sub-basinindicate that none of the studied discharges need load reduction to achieve EQSs. However, in various points, tridecane, nickel, bis(2-ethylhexyl) terephthalate, NH4-N, total phosphorus, and free cyanide have consumed more than 10% of their discharge quotas estimated by the mentioned approach. Therefore, for the sub-basin, these six contaminants and their corresponding two discharge points may require more attention in the future.

Seasonal to Decadal Variability in Focused Groundwater and Contaminant Discharge along a Channelized Stream

AbstractFocused groundwater discharge to streams is problematic at contaminated sites because high fluxes can limit natural attenuation in the hyporheic zone. However, information on location, spatial evolution, and temporal persistence of springs in unlithified sediments over multiyear time scales is limited. We examine discharge at point (~1‐m) to reach (~300‐m) scales along a stream that intercepts trichloroethene and technetium‐99 plumes from a Superfund site. During 2011 to 2012, we seasonally monitored stream and spring flow and contaminant concentrations, along with probing streambed temperatures on a grid in winter and summer, building on prior monitoring during 1999 to 2002. Baseflow measured by both gauging and dye dilution generally increased with distance downstream, and stream and spring discharge varied seasonally, from minima in October to January to maxima in February to June. Thermal anomalies identified by probing occupied approximately 3% to 6% of the reach and typically coincided with visible springs or seeps. Locations of anomalies were similar to those identified in summer 2002, although some orifices disappeared and others emerged. Vertical groundwater fluxes calculated from probing tended to be less than net fluxes calculated from stream discharge, perhaps in part because the assumption of one‐dimensional, steady‐state flow in calculating point fluxes was simplistic. Maximum contaminant concentrations and fluxes decreased between 1999 to 2001 and 2011 to 2012 as a result of partial capture by an upgradient pump‐and‐treat system. Our findings confirm that springs in unlithified sediments can remain stationary within a few meters over decadal time scales, and seasonal variability in discharge can be greater than decadal variability.

Field performance of the radon-deficit technique to detect and delineate a complex DNAPL accumulation in a multi-layer soil profile.

The performance of the radon (222Rn)-deficit technique has been evaluated at a site in which a complex DNAPL mixture (mostly hexachlorocyclohexanes and chlorobenzenes) has contaminated all four layers (from top to bottom: anthropic backfill, silt, gravel and marl) of the soil profile. Soil gas samples were collected at two depths (0.8m and 1.7m) in seven field campaigns and a total of 186 222Rn measurements were performed with a pulse ionization detector. A statistical assessment of the influence of field parameters on the results revealed that sampling depth and atmospheric pressure did not significantly affect the measurements, while the location of the sampling point and ground-level atmospheric temperature did. In order to remove the bias introduced by varying field temperatures and hence to be able to jointly interpret 222Rn measurements from different campaigns, 222Rn concentrations were rescaled by dividing each individual datum by the mean 222Rn concentration of its corresponding field campaign. Rescaled 222Rn maps showed a high spatial correlation between 222Rn minima and maximum contaminant concentrations in the top two layers of the soil profile, successfully delineating the surface trace of DNAPL accumulation in the anthropic backfill and silt layers. However, no correlation could be established between 222Rn concentrations in superficial soil gas and contaminant concentration in the deeper two layers of the soil profile. These results indicate that the 222Rn-deficit technique is unable to describe the vertical variation of contamination processes with depth but can be an effective tool for the preliminary characterization of sites in which the distance between the inlet point of the sampling probe and the contaminant accumulation falls within the effective diffusion length of 222Rn in the affected soil profile.

Method to Reduce Aerosolized Contaminant Concentration Exposure to Healthcare Workers During the COVID-19 Pandemic when Temporary Isolation Systems Are Required.

The COVID-19 pandemic has strained the healthcare system. It has led to the use of temporary isolation systems and less-then-optimum patient placement configurations because of inadequate number of isolation rooms, both of which can compromise provider safety. Three key elements require special attention to reduce the maximum and average aerosolized contaminant concentration exposure to a healthcare worker in any isolation system: flow rate; air changes per hour; and patient placement. This is important because concentration exposures of aerosolized contaminants to healthcare workers in hospitals using temporary isolation systems can reach levels 21–30 times greater than a properly engineered negative pressure isolation room. A working knowledge of these three elements can help create a safer environment for healthcare workers when isolation rooms are not available.

Groundwater contamination in mega cities with finite sources

Groundwater contamination due to multiple sources occurring in mega cities was modelled. One constant source contamination was considered at the source boundary, whereas other sources may join in between at various locations at different times. Initially, the aquifer was contamination-free in mega cities and was subsequently contaminated by means of different sources in due course of time. One-dimensional ADE (Advection Dispersion Equation) for modelling groundwater contamination was used and solved analytically in the semi-infinite aquifer domain for a finite number of point sources. A numerical solution was also obtained for two sources to compare analytical solutions. Results were examined for different velocity profiles to show the maximum contaminant concentration level with distance. This may be helpful to model the maximum possible number of point sources of contamination (i.e., it represents approximately what happens in the field situation). Some remedial measures may be taken to overcome these kinds of contamination problems in mega cities by treating the sources so that recharge of the aquifer is not affected.

A CFD analysis of the pollutant dispersion from cooling towers with various configurations in the lower region of atmospheric boundary layer

The mixing and merging of buoyant plumes originating from multiple small cooling towers into the atmosphere is numerically investigated. The effects of different arrangements of cooling towers as well as outlet geometries on the mixing of the plumes are examined. The side by side and tandem arrangements of two sources and also two types of multi-flue cooling towers are considered. The various ways by which the counter rotating vortex pair, as the dominant mechanism, affect the flow pattern in each aforementioned configuration are investigated. For tandem arrangement, far from sources, the outlet flow of the downstream cooling tower surrounds the plumes originating from the upstream cooling tower and in the region near the cooling towers, the pollutants are mostly originated from the upstream cooling tower. Maximum pollutant concentrations at distances 10 and 40 times the diameter downstream of the leading cooling tower increase by 67% and 29% with respect to those of a single cooling tower, respectively. For the side by side arrangement, the counter rotating vortices are stretched due to the large low pressure area created downstream of the cooling towers. Mixing of the plumes with the surrounding air is reduced as a result of contraction of vortices. Maximum contaminant concentrations at distances 10 and 40 times the diameter downstream of the cooling towers increase by 29% and 41% with respect to those of a single tower, respectively. Finally, the differences between flow fields formed around diamond and square configurations of multi-flue cooling towers are extensively discussed.

Simultaneous determination of multiresidue and multiclass emerging contaminants in waters by rotating-disk sorptive extraction–derivatization-gas chromatography/mass spectrometry

An efficient method has been developed for the multiresidue and multiclass determination of 16 emerging contaminants (parabens, hormones, anti-inflammatory drugs, triclosan and bisphenol A) in water samples using rotating-disk sorptive extraction (RDSE) and gas chromatography coupled to mass spectrometry (GC-MS). Silylation of the compounds prior to GC-MS analysis was optimized using a factorial experimental design; the optimal derivatization conditions to maximize the response of the set of analytes included 70 μL of N-methyl-N-(trimethylsilyl)trifluoroacetamide at 80 °C for 35 min. RDSE was implemented using Oasis® HLB as a sorptive phase and an extraction time of 60 min.The method was applied to Chilean environmental samples. In tap water, none of the analytes under study were detected. In the river and well waters, the concentrations of the four detected contaminants were below 0.38 μg L−1. In the effluent and influent of the wastewater treatment plant, the maximum concentrations of contaminants were 3.1 and 4.2 μg L−1, respectively.The proposed analytical strategy suggests clear improvements with respect to other methods reported in the literature, considering not only the different steps involved in the analytical process (extraction, derivatization and chromatography) but also taking into account that this method involves the determination of different families of analytes with different physicochemical and structural properties.

Cadmium, lead, and arsenic contamination in paddy soils of a mining area and their exposure effects on human HEPG2 and keratinocyte cell-lines

A mining district in south China shows significant metal(loid) contamination in paddy fields. In the soils, average Pb, Cd and As concentrations were 460.1, 11.7 and 35.1mgkg−1 respectively, which were higher than the environmental quality standard for agricultural soils in China (GB15618-1995) and UK Clea Soil Guideline Value. The average contents of Pb, Cd and As in rice were 5.24, 1.1 and 0.7mgkg−1 respectively, which were about 25, 4.5 or 2.5 times greater than the limit values of the maximum safe contaminant concentration standard in food of China (GB 2762-2012), and about 25, 10 or 1 times greater than the limit values of FAO/WHO standard. The elevated contents of Pb, Cd and As detected in soils around the factories, indicated that their spatial distribution was influenced by anthropogenic activity, while greater concentrations of Cd in rice appeared in the northwest region of the factories, indicating that the spatial distribution of heavy metals was also affected by natural factors. As human exposure around mining districts is mainly through oral intake of food and dermal contact, the effects of these metals on the viability and MT protein of HepG2 and KERTr cells were investigated. The cell viability decreased with increasing metal concentrations. Co-exposure to heavy metals (Pb+Cd) increased the metals (Pb or Cd)-mediated MT protein induction in both human HepG2 and KERTr cells. Increased levels of MT protein will lead to greater risk of carcinogenic manifestations, and it is likely that chronic exposure to metals may increase the risk to human health. Nevertheless, when co-exposure to two or more metals occur (such as As+Pb), they may have an antagonistic effect thus reducing the toxic effects of each other. CapsuleMetal contaminations in paddy soils and rice were influenced by anthropogenic activity; metal co-exposure induced MT protein in human cells.

An Experimental Study on the Application of Ultrasonic Technology for Demulsifying Crude Oil and Water Emulsions

An emulsion is the mixture of two immiscible fluids, where one fluid appears as droplets within another. In the oil and gas industry, produced crude oil generally comes with an appreciable amount of water within it in an emulsified form. Before produced crude oil can be prepared for purchase, the water associated with it must be removed. A process known as demulsification is required in order to separate an emulsion into its two phases. In the industry, a number of demulsification techniques are already present; these include thermal, mechanical, chemical, and electrical techniques. Crude oil and gas produced from wells originally come with water, salts, and volatile gases such as oxygen, carbon dioxide, and sometimes hydrogen sulfide, etc. Hence, the petroleum mixture needs to be refined-water, salt, and non-hydrocarbon gases to be separated from the mixture, in order to meet certain oil and gas specifications (which state the maximum concentrations of such contaminants) and make it ready for purchase and transportation. Sonication provides a cheap, simple, and harmless (as it involves mainly the propagation of sound waves) way of separating crude oils from water droplets via demulsification. In addition, if needed, it can be used for emulsification processes as well. Hence, a study of sonification as a way for crude refinement or chemical mixing has important implications for the oil and gas. This investigation proposes the use of ultrasonication as a new and cost-effective technique to aid in the demulsification of crude oil emulsion. The effectiveness of this technique was gauged through its comparison to the already present methods in the industry. Based on the investigation it was found that centrifuge served as the best demulsification method for it reduced the turbidity by 86%. In addition, the reduced turbidity achieved with proposed ultrasonication method ranges from 20%-60%.

Using turbidity for designing water networks

Some methods to design water networks with minimum fresh water consumption are based on the selection of a key contaminant. In most of these “single contaminant methods”, a maximum allowable concentration of contaminants must be established in water demands and water sources. Turbidity is not a contaminant concentration but is a property that represents the “sum” of other contaminants, with the advantage that it can be cheaper and easily measured than biological oxygen demand, chemical oxygen demand, suspended solids, dissolved solids, among others. The objective of this paper is to demonstrate that turbidity can be used directly in the design of water networks just like any other contaminant concentration. A mathematical demonstration is presented and in order to validate the mathematical results, the design of a water network for a guava fudge production process is performed. The material recovery pinch diagram and nearest neighbors algorithm were used for the design of the water network. Nevertheless, this water network could be designed using other single contaminant methodologies. The maximum error between the expected and the real turbidity values in the water network was 3.3%. These results corroborate the usefulness of turbidity in the design of water networks.

Determination of Hexavalent Chromium (Cr(VI)) Concentrations via Ion Chromatography and UV-Vis Spectrophotometry in Samples Collected from Nacogdoches Wastewater Treatment Plant, East Texas (USA)

The concentration of hexavalent chromium (Cr(VI)), a toxic environmental pollutant and carcinogen, was determined in samples collected from Nacogdoches Wastewater Treatment Plant (NWWTP) using ion chromatography and UV-visible spectrophotometry (IC, UV-Vis). On reaction with 1,5-diphenylcarbazide (DPC) Cr+6 forms a 1,5-diphenylcarbazide-Cr(VI) complex, which is then analyzed at 530 nm and 540 nm, respectively. Via ion chromatography Cr(VI) concentrations were in the range of 0.00190±0.0020 and 0.0010±0.0006 ppm at the influent and effluent, respectively. With the use of standard addition wastewater samples were spiked with a 0.5 ppm Cr(VI) standard of various amounts and subsequently analyzed with UV-Vis spectrophotometry. The spiked concentrations gave Cr(VI) concentrations in the range of 0.0090±0.0060 ppm and 0.0040±0.0061 ppm at the influent and influent wastewater, respectively. The determined Cr(VI) concentrations through the ion chromatography and UV-Vis spectrophotometry are below the maximum USEPA contaminant concentration of 0.1 ppm. From the analysis, the NWWTP efficiently removes Cr(VI) before discharge into the environment through La Nana Creek. The removal efficiency for Cr(VI) was determined to be ≥92.8% along the wastewater treatment stages from the influent (aeration stage) to the effluent stages prior to discharge into the La Nana Creek.

The risk of river pollution due to washout from contaminated floodplain water bodies during periods of high magnitude floods

Summary The risk of river pollution due to washout (removal of pollutants) from contaminated floodplain water bodies (floodplain lakes and quarries whose origin is related to the large-scale mining of nonmetallic building materials in the floodplain zone) during high magnitude flood periods is analyzed using a combination of one-, two- and three-dimensional hydrodynamic modeling and in situ measurements. The modeling performed for the floodplain water bodies contaminated by N compounds shows that during large magnitude floods washout occurs. The washout process consists of two stages: an initial rapid stage lasting about two hours during which the upper (3–4 m thick) layer is washed out, followed by a second stage when the concentration of NH4-N in the floodplain water body remains nearly constant. The maximum contaminant concentration in the river in the vicinity of a water intake for drinking water located 21 km downstream is attained about 9 h from the beginning of the flood; concentration of NH4-N can reach values several times larger than acceptable concentration guidelines. The initial primary peak in contaminant concentration at the water intake is followed by a slight decrease in contaminant concentration; a second peak related to the contaminant transport through the inundated floodplain subsequently occurs, after which the concentration slowly decreases, reaching acceptable values after 30–40 h. Contaminated floodplain water bodies located near drinking water supply systems are not significant sources of contamination during small and moderate floods, but during high magnitude floods, they can become sources of water pollution. Operational measures that can decrease potential health risks are discussed.

Analysis of regulatory guidance values for six commonly regulated chlorophenols

Regulatory guidance values (RGVs) are developed by regulators to protect human health by establishing maximum allowable contaminant concentrations. Contamination of residential surface soil is of particular concern because of its potential impact on the general population, especially children. As part of an ongoing examination of soil RGVs, this study analyses guidance values applied to chlorophenols, a class of chemicals used as pesticides, herbicides, and disinfectants. Chlorophenol RGVs have been promulgated by at least 45 US states, 38 United Nations member states, four Native American tribes, five federal agencies, and several city, county, and territorial jurisdictions. An analysis of the six most commonly regulated chlorophenols (2-chlorophenol, 2,4-dichlorophenol, 2,4,5-trichlorophenol, 2,4,6-trichlorophenol, 2,3,4,6-tetrachlorophenol, and pentachlorophenol) demonstrates that there is great variability in the guidance values, which vary by six to eight orders of magnitude and have distributions similar to those of lognormal random variables. Clearly, there is little universal agreement on the appropriate maximum contamination levels for chlorophenols.

Evaluation de la vulnérabilité d’un puits pour la protection d’un captage d’eau souterraine basée sur un modèle de transport de soluté: étude de cas au niveau de la ville de Jilin, Nord-Est de la Chine

Well vulnerability assessment is essential for groundwater source protection. A quantitative approach to assess well vulnerability in a well capture zone is presented, based on forward solute transport modeling. This method was applied to three groundwater source areas (Jiuzhan, Hadawan and Songyuanhada) in Jilin City, northeast China. The ratio of the maximum contaminant concentration at the well to the released concentration at the contamination source (c max/c 0) was determined as the well vulnerability indicator. The results indicated that well vulnerability was higher close to the pumping well. The well vulnerability in each groundwater source area was low. Compared with the other two source areas, the cone of depression at Jiuzhan resulted in higher spatial variability of c max/c 0 and lower minimum c max/c 0 by three orders of magnitude. Furthermore, a sensitivity analysis indicated that the denitrification rate in the aquifer was the most sensitive with respect to well vulnerability. A process to derive a NO3−N concentration at the pumping well is presented, based on determining the maximum nitrate loading limit to satisfy China’s drinking-water quality standards. Finally, the advantages, disadvantages and prospects for improving the precision of this well vulnerability assessment approach are discussed.