The concept of vulnerability of drinking water sources is reviewed, and a quantitative approach for assessing well vulnerability for complex three-dimensional ground water systems is developed. The approach focuses on the relative expected impact of potential contaminant sources at unknown locations within a well capture zone, providing relative measures of intrinsic well vulnerability, including the expected times of arrival of a contaminant, the dispersion-related reduction in concentration, the time taken to breach a certain quality objective, and the corresponding exposure times. Thus, the result of the analysis includes the usual advective travel time information used in conventional wellhead protection analysis, plus a set of selected quantitative measures expressing the expected impact. The technique is based on adjoint theory and combines forward- and backward-in-time transport modeling using a standard numerical flow and transport code. The methodology is demonstrated using the case study of a complex glacial multiaquifer system in Ontario. The new approach will be useful in helping water managers develop more physically based and quantitative wellhead protection strategies.

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers. L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

The process that transfers mass from a subsurface source zone of residual dense nonaqueous phase liquid (DNAPL) to the flowing groundwater is a controlling factor in determining the time required to dissolve the source by noninvasive means. While mass transfer can be kinetic or equilibrium under laboratory conditions, aqueous concentrations in the field are generally found to be below equilibrium levels. To gain insight into the mass transfer process under field conditions, we simulated the dissolution of the emplaced DNAPL source at the Canadian Forces Base Borden, Ontario, which contains a mixture of three DNAPLs. The simulations clearly show that mass transfer at this site is equilibrium‐controlled during the 1000‐day observation period and that apparent tailing of one of the organic components is due to its declining solubility, rather than mass transfer kinetics. Flow lines passing through the source are focused in a narrow streamtube downstream of the source, and equilibrium concentrations are therefore observed only at the center of the effluent plume. Since the concentration peaks can be easily missed in the sampling, streamline focusing can explain the low concentrations observed in the field.

Abstract. The spatial distribution of groundwater fluxes through a streambed can be highly variable, most often resulting from a heterogeneous distribution of aquifer and streambed permeabilities along the flow pathways. Using a groundwater flow and heat transport model, we defined four scenarios of aquifer and streambed permeability distributions to simulate and assess the impact of subsurface heterogeneity on the distribution of groundwater fluxes through the streambed: (a) a homogeneous low-K streambed within a heterogeneous aquifer; (b) a heterogeneous streambed within a homogeneous aquifer; (c) a well connected heterogeneous low-K streambed within a heterogeneous aquifer; and (d) a poorly connected heterogeneous low-K streambed within a heterogeneous aquifer. The simulation results were compared with a base case scenario, in which the streambed had the same properties as the aquifer, and with observed data. The results indicated that the aquifer has a stronger influence on the distribution of groundwater fluxes through the streambed than the streambed itself. However, a homogeneous low-K streambed, a case often implemented in regional-scale groundwater flow models, resulted in a strong homogenization of fluxes, which may have important implications for the estimation of peak mass flows. The flux distributions simulated with heterogeneous low-K streambeds were similar to the flux distributions of the base case sceCorrespondence to: E. Kalbus (edda.kalbus@web.de) nario, despite the lower permeability. The representation of heterogeneous distributions of aquifer and streambed properties in the model has been proven to be beneficial for the accuracy of flow simulations.

MINTRAN is a new model for simulating transport of multiple thermodynamically reacting chemical substances in groundwater systems. It consists of two main modules, a finite element transport module (PLUME2D), and an equilibrium geochemistry module (MINTEQA2). Making use of the local equilibrium assumption, the inherent chemical nonlinearity is confined to the chemical domain. This linearizes the coupling between the physical and chemical processes and leads to a simple and efficient two‐step sequential solution algorithm. The advantages of the coupled model include access to the comprehensive geochemical database of MINTEQA2 and the ability to simulate hydrogeological systems with realistic aquifer properties and boundary conditions under complex geochemical conditions. The model is primarily targeted toward groundwater contamination due to acidic mine tailings efiïuents but is potentially also applicable to the full range of geochemical scenarios covered by MINTEQA2. The model is tested with respect to ion exchange chemistry and with respect to precipitation/dissolution chemistry involving multiple sharp fronts. The companion paper presents two‐dimensional simulations of heavy metal transport in an acidic mine tailings environment, focusing on environmental implications.