Evaluation of Hydraulic Conductivity Estimates from Various Approaches with Groundwater Flow Models.

Significant efforts have been expended for improved characterization of hydraulic conductivity (K) and specific storage (Ss ) to better understand groundwater flow and contaminant transport processes. Conventional methods including grain size analyses (GSA), permeameter, slug, and pumping tests have been utilized extensively, while Direct Push-based Hydraulic Profiling Tool (HPT) surveys have been developed to obtain high-resolution K estimates. Moreover, inverse modeling approaches based on geology-based zonations, and highly parameterized Hydraulic Tomography (HT) have also been advanced to map spatial variations of K and Ss between and beyond boreholes. While different methods are available, it is unclear which one yields K estimates that are most useful for high resolution predictions of groundwater flow. Therefore, the main objective of this study is to evaluate various K estimates at a highly heterogeneous field site obtained with three categories of characterization techniques including: (1) conventional methods (GSA, permeameter, and slug tests); (2) HPT surveys; and (3) inverse modeling based on geology-based zonations and highly parameterized approaches. The performance of each approach is first qualitatively analyzed by comparing K estimates to site geology. Then, steady-state and transient groundwater flow models are employed to quantitatively assess various K estimates by simulating pumping tests not used for parameter estimation. Results reveal that inverse modeling approaches yield the best drawdown predictions under both steady and transient conditions. In contrast, conventional methods and HPT surveys yield biased predictions. Based on our research, it appears that inverse modeling and data fusion are necessary steps in predicting accurate groundwater flow behavior.

An Electron-Balance Based Approach to Predict the Decreasing Denitrification Potential of an Aquifer.

Numerical models for reactive transport can be used to estimate the breakthrough of a contaminant in a pumping well or at other receptors. However, as natural aquifers are highly heterogeneous with unknown spatial details, reactive transport predictions on the aquifer scale require a stochastic framework for uncertainty analysis. The high computational demand of spatially explicit reactive-transport models hampers such analysis, thus motivating the search for simplified estimation tools. We suggest performing an electron balance between the reactants in the infiltrating solution and in the aquifer matrix to obtain the hypothetical time of dissolved-reactant breakthrough at a receptor if the reaction with the matrix was instantaneous. This time we denote as the advective breakthrough time for instantaneous reaction (τinst ). It depends on the amount of the reaction partner present in the matrix, the mass flux of the dissolved reactant, and the stoichiometry. While the shape of the reactive-species breakthrough curve depends on various kinetic parameters, the overall timing scales with τinst . We calculate the latter by particle tracking. The effort of computing τinst is so low that stochastic calculations become feasible. We apply the concept to a two-dimensional test case of aerobic respiration and denitrification. A detailed spatially explicit reactive-transport model includes microbial dynamics. Scaling the time of local breakthrough curves observed at individual points by τinst decreased the variability of electron-donor breakthrough curves significantly. We conclude that the advective breakthrough time for instantaneous reaction is efficient in estimating the time over which an aquifer retains its degradation potential.

Delineating baseflow contribution areas for streams - A model and methods comparison.

This study addresses the delineation of areas that contribute baseflow to a stream reach, also known as stream capture zones. Such areas can be delineated using standard well capture zone delineation methods, with three important differences: (1) natural gradients are smaller compared to those produced by supply wells and are therefore subject to greater numerical errors, (2) stream discharge varies seasonally, and (3) stream discharge varies spatially. This study focuses on model-related uncertainties due to model characteristics, discretization schemes, delineation methods, and particle tracking algorithms. The methodology is applied to the Alder Creek watershed in southwestern Ontario. Four different model codes are compared: HydroGeoSphere, WATFLOW, MODFLOW, and FEFLOW. In addition, two delineation methods are compared: reverse particle tracking and reverse transport, where the latter considers local-scale parameter uncertainty by using a macrodispersion term to produce a capture probability plume. The results from this study indicate that different models can calibrate acceptably well to the same data and produce very similar distributions of hydraulic head, but can produce different capture zones. The stream capture zone is found to be highly sensitive to the particle tracking algorithm. It was also found that particle tracking by itself, if applied to complex systems such as the Alder Creek watershed, would require considerable subjective judgement in the delineation of stream capture zones. Reverse transport is an alternative and more reliable approach that provides probability intervals for the baseflow contribution areas, taking uncertainty into account. The two approaches can be used together to enhance the confidence in the final outcome.

An integrated approach for addressing uncertainty in the delineation of groundwater management areas.

Uncertainty is a pervasive but often poorly understood factor in the delineation of wellhead protection areas (WHPAs), which can discourage water managers and practitioners from relying on model results. To make uncertainty more understandable and thereby remove a barrier to the acceptance of models in the WHPA context, we present a simple approach for dealing with uncertainty. The approach considers two spatial scales for representing uncertainty: local and global. At the local scale, uncertainties are assumed to be due to heterogeneities, and a capture zone is expressed in terms of a capture probability plume. At the global scale, uncertainties are expressed through scenario analysis, using a limited number of physically realistic scenarios. The two scales are integrated by using the precautionary principle to merge the individual capture probability plumes corresponding to the different scenarios. The approach applies to both wellhead protection and the mitigation of contaminated aquifers, or in general, to groundwater management areas. An example relates to the WHPA for a supply well located in a complex glacial aquifer system in southwestern Ontario, where we focus on uncertainty due to the spatial distributions of recharge. While different recharge scenarios calibrate equally well to the same data, they result in different capture probability plumes. Using the precautionary approach, the different plumes are merged into two types of maps delineating groundwater management areas for either wellhead protection or aquifer mitigation. The study shows that calibrations may be non-unique, and that finding a "best" model on the basis of the calibration fit may not be possible.

A simple method to assess unsaturated zone time lag in the travel time from ground surface to receptor.

In contaminant travel from ground surface to groundwater receptors, the time taken in travelling through the unsaturated zone is known as the unsaturated zone time lag. Depending on the situation, this time lag may or may not be significant within the context of the overall problem. A method is presented for assessing the importance of the unsaturated zone in the travel time from source to receptor in terms of estimates of both the absolute and the relative advective times. A choice of different techniques for both unsaturated and saturated travel time estimation is provided. This method may be useful for practitioners to decide whether to incorporate unsaturated processes in conceptual and numerical models and can also be used to roughly estimate the total travel time between points near ground surface and a groundwater receptor. This method was applied to a field site located in a glacial aquifer system in Ontario, Canada. Advective travel times were estimated using techniques with different levels of sophistication. The application of the proposed method indicates that the time lag in the unsaturated zone is significant at this field site and should be taken into account. For this case, sophisticated and simplified techniques lead to similar assessments when the same knowledge of the hydraulic conductivity field is assumed. When there is significant uncertainty regarding the hydraulic conductivity, simplified calculations did not lead to a conclusive decision.

Nitrate controls methyl mercury production in a streambed bioreactor.

Organic carbon bioreactors provide low-cost, passive treatment of a variety of environmental contaminants but can have undesirable side effects in some cases. This study examines the production of methyl mercury (MeHg) in a streambed bioreactor consisting of 40 m³ of wood chips and designed to treat nitrate (NO3) in an agricultural drainage ditch in southern Ontario (Avon site). The reactor provides 30 to 100% removal of NO3-N concentrations of 0.6 to 4.4 mg L(-1), but sulfate (SO4(2-)) reducing conditions develop when NO3 removal is complete. Sulfate reducing conditions are known to stimulation the production of MeHg in natural wetlands. Over one seasonal cycle, effluent MeHg ranged from 0.01 to 0.76 ng L(-1) and total Hg ranged from 1.3 to 3.4 ng L(-1). During all sampling events when reducing conditions were only sufficient to promote NO3(-) reduction (or denitrification) ( = 5, late fall 2009, winter 2010), MeHg concentrations decreased in the reactor and it was a net sink for MeHg (mean flux of -5.1 µg m(-2) yr(-1)). During all sampling events when SO4(2-) reducing conditions were present ( = 6, early fall 2009, spring 2010), MeHg concentrations increased in the reactor and it was a strong source of MeHg to the stream (mean flux of 15.2 µg m(-2) yr(-1)). Total Hg was consistently removed in the reactor (10 of 11 sampling events) and was correlated to the total suspended sediment load ( r² = 0.69), which was removed in the reactor by physical filtration. This study shows that organic carbon bioreactors can be a strong source of MeHg production when SO4(2-) reducing conditions develop; however, maintaining NO3-N concentrations > 0.5 mg L suppresses the production of MeHg.

The utility of microspheres as surrogates for the transport of E. coli RS2g in partially saturated agricultural soil.

Polystyrene latex microspheres are widely used as surrogates for biocolloid transport in porous media; however, relatively few studies directly compare microsphere transport with that of the microorganism it is intended to represent, particularly at the field scale. Here, we compared the transport behaviour of a bacterium (Escherichia coli RS2g; 1.2 microm in diameter) and three different sized microspheres (1.1, 3.9, and 4.8 microm in diameter) within undisturbed agricultural field soil following infiltration under partially saturated conditions. The soil contained significant macroporosity. A tension infiltrometer was used to control the application of a transport solution containing Brilliant Blue FCF dye to two plots. A >2 log reduction in the concentration of all colloids was observed from the soil surface to 5 cm depth in both plots. The concentration of colloids in the soil was generally proportional to the intensity of soil dye staining; however, both the E. coli RS2g bacterium and the 1.1 microm microspheres appeared to be transported deeper than the other colloids and the visible dye along root holes at the bottom of the profile in both plots. The similarities in size and zeta potential of the 1.1 microm microspheres and the E. coli RS2g likely contributed to that outcome. Colloid concentrations in dyed soil by depth were similar between the two plots, despite differences in soil properties and infiltration patterns. The properties of the colloids and macropore density were the most important factors affecting colloid transport. These results suggest that microspheres with size and surface properties similar to the microbe of interest are useful surrogates to trace potential pathways of transport in the subsurface.

Influence of macroporosity on preferential solute and colloid transport in unsaturated field soils.

Transport of solutes and colloids in soils, particularly those subject to preferential flow along macropores, is important for assessing the vulnerability of shallow groundwater to contamination. The objective of this study was to investigate flow and transport phenomena for dissolved and colloid tracers during large infiltration events in partially saturated, macroporous soils. Controlled tracer infiltration tests were completed at two field sites in southern Ontario. A tension infiltrometer (TI) was used to infiltrate water with dissolved Brilliant Blue FCF dye simultaneously with 3.7 microm and 0.53 microm diameter fluorescent microspheres. Infiltration was conducted under maximum infiltration pressure heads ranging from -5.2 to -0.4 cm. All infiltration test sites were excavated to examine and photograph dye-stained flow patterns, map soil features, and collect samples for microsphere enumeration. Results indicated that preferential transport of dye and microspheres via macropores occurred when maximum pressure heads were greater than -3.0 cm, and the corresponding infiltration rates exceeded 2.0 cm h(-1). Dye and microspheres were detected at depths greater than 70 cm under the highest infiltration rates from both sites. Microsphere concentrations in the top 5-10 cm of soil decreased by more than two orders of magnitude relative to input concentrations, yet remained relatively constant with depth thereafter. There was some evidence for increased retention of the 3.7 microm microspheres relative to the 0.53 microm microspheres, particularly at lower infiltration pressures where straining and attachment mechanisms are most prevalent. Microspheres were observed within dye stained soil matrix surrounding individual macropores, illustrating the significance of capillary pressures in controlling the vertical migration of both tracers in the vicinity of the macropores. Overall, microsphere distributions closely followed the dye patterns, with microsphere concentrations at all depths directly related to the intensity (or concentration) of dye staining. It is concluded that the flow system influenced transport to a much greater degree than differences between dissolved and colloidal species, and hence a dye tracer could serve as a reasonable surrogate for colloid distributions in the vadose zone following individual infiltration events.