Consequences of Groundwater-Model Vertical Discretization in Risk-Based Decision-Making

One of the first and most important decisions facing practitioners when constructing a numerical groundwater model is vertical discretization. Several factors will influence this decision, such as the conceptual model of the system and hydrostratigraphy, data availability, resulting computational burden, and the purpose of the modeling analysis. Using a coarse vertical discretization is an attractive option for practitioners because it reduces data requirements and model construction efforts, can increase model stability, and can reduce computational demand. However, using a coarse vertical discretization as a form of model simplification is not without consequence; this may give rise to unwanted side-effects such as biases in decision-relevant simulated outputs. Given its foundational role in the modeled representation of the aquifer system, herein we investigate how vertical discretization may affect decision-relevant simulated outputs using a paired complex-simple model analysis. A Bayesian framework and decision analysis approach are adopted. Two case studies are considered, one of a synthetic, linked unsaturated-zone/surface-water/groundwater hydrologic model and one of a real-world linked surface-water/groundwater hydrologic-nitrate transport model. With these models, we analyze decisions related to abstraction-induced changes in ecologically important streamflow characteristics and differences in groundwater and surface-water nitrate concentrations and mass loads following potential land-use change. We show that for some decision-relevant simulated outputs, coarse vertical discretization induces bias in important simulated outputs, and can lead to incorrect resource management action. For others, a coarse vertical discretization has little or no consequence for resource management decision-making.     

Numerical Modeling of Nitrate Removal in Anoxic Groundwater during River Flooding of Riparian Zones

A numerical study demonstrates the effects of flooding on subsurface hydrological flowpaths and nitrate removal in anoxic groundwater in riparian zones with a top peat layer. A series of two-dimensional numerical simulations with changing conditions for flow (steady state or transient with flooding), hydrogeology, denitrification, and duration of flooding demonstrate how flowpaths, residence times, and nitrate removal are affected. In periods with no flooding groundwater flows horizontally and discharges to the river through the riverbed. During periods with flooding, shallow groundwater is forced upwards as discharge through peat layers that often have more optimal conditions for denitrification caused by the presence of highly reactive organic matter. The contrast in hydraulic conductivity between the sand aquifer and the overlying peat layer, as well as the flooding duration, have a significant role in determining the degree of nitrate removal.     

Cyclodextrin-Enhanced Vertical Flushing of a Trichloroethene Contaminated Aquifer

Pilot-scale testing of an innovative ground water remediation technology was conducted in a source zone of a trichloroethene-contaminated Superfund site in Tucson, Arizona. The technology is designed to enhance the removal of low-solubility organic contaminants from heterogeneous sedimentary aquifers by using a dual-screened vertical circulation well to inject and extract solutions containing a complexing sugar (hydroxypropyl-beta-cyclodextrin (HPCD]). Prior to initiating the pilot test, tracer tests were conducted to determine hydraulic characteristics of the vertical flow field and to evaluate trichloroethene-elution behavior during water flushing. The pilot test involved injecting approximately 4 m³ of a 20% HPCD solution into the upper screened interval of the well and extracting from the lower screened interval. The results of the pilot test indicate that the cyclodextrin solution increased the rate of trichloroethene removal from the aquifer. The concentrations of trichloroethene in the ground water extracted from the lower screened interval of the well increased by a factor of three (∼750 μg/L) in the presence of the cyclodextrin pulse, compared to concentrations obtained during previous water flushing (∼250 μg/L). Furthermore, the concentration of trichloroethene in water collected from the circulation well under static conditions was reduced to 6% of the levels measured prior to the test.     

Spectroscopic Characterization of Fulvic Acid Fractions of a Contaminated Groundwater

Two fulvic acid (FA) samples taken from a former gas production facility in the Southwest of Germany were characterized using advanced fluorescence techniques. Steady-state fluorescence (fluorescence excitation, synchronous fluorescence) as well as time-resolved fluorescence were applied. Distinct differences between the sample B22 FA taken within the contamination plume and the sample B53 FA taken downstream were found. Comparison with a model compound for metabolites and humic substances revealed that due to the downstream passage the characteristics of the dissolved organic matter became more humic-like. The assignment of single classes of compounds in the sample B22 FA is discussed in terms of their synchronous fluorescence spectra and fluorescence decay time distribution.     

Including Vertical Fault Structures in Layered Groundwater Flow Models

Accurate representation of groundwater flow and solute transport requires a sound representation of the underlying geometry of aquifers. Faults can have a significant influence on the structure and connectivity of aquifers, which may allow permeable units to connect, and aquifers to seal when juxtaposed against lower permeability units. Robust representation of groundwater flow around faults remains challenging despite the significance of faults for flow and transport. We present a methodology for the inclusion of faults utilizing the unstructured grid features of MODFLOW-USG and MODFLOW 6. The method focuses on the representation of fault geometries using non-neighbor connections between juxtaposed layers. We present an illustration of the method for a synthetic fluvial aquifer. The combined impact of the heterogeneous aquifer and fault offset is clearly visible where channel features at different depths in the aquifer were connected at the fault. These results highlight the importance of representing fault features in groundwater flow models.     

Importance of Spatial Resolution in Global Groundwater Modeling

Global-scale gradient-based groundwater models are a new endeavor for hydrologists who wish to improve global hydrological models (GHMs). In particular, the integration of such groundwater models into GHMs improves the simulation of water flows between surface water and groundwater and of capillary rise and thus evapotranspiration. Currently, these models are not able to simulate water table depth adequately over the entire globe. Unsatisfactory model performance compared to well observations suggests that a higher spatial resolution is required to better represent the high spatial variability of land surface and groundwater elevations. In this study, we use New Zealand as a testbed and analyze the impacts of spatial resolution on the results of global groundwater models. Steady-state hydraulic heads simulated by two versions of the global groundwater model G³M, at spatial resolutions of 5 arc-minutes (9 km) and 30 arc-seconds (900 m), are compared with observations from the Canterbury region. The output of three other groundwater models with different spatial resolutions is analyzed as well. Considering the spatial distribution of residuals, general patterns of unsatisfactory model performance remain at the higher resolutions, suggesting that an increase in model resolution alone does not fix problems such as the systematic overestimation of hydraulic head. We conclude that (1) a new understanding of how low-resolution global groundwater models can be evaluated is required, and (2) merely increasing the spatial resolution of global-scale groundwater models will not improve the simulation of the global freshwater system.     

太湖风生流数值模式及其对总磷分布的影响

A new up-winding finite element numerical model, which is two-dimensional and suitable for modeling lake current and the distribution of total phosphorus(TP) in shallow lakes, is derived. Moreover, it is used in the study of wind-driven current in Taihu Lake, and the impact of lake current field on the distribution of TP is also discussed.     

Examining the Extraction Efficiency of Petroleum-Derived Dissolved Organic Matter in Contaminated Groundwater Plumes

The extraction efficiency of petroleum-derived dissolved organic matter (DOM) was examined for groundwater samples from an aquifer contaminated with crude oil. Five different types of extraction techniques were investigated to determine which method is best suited for the analysis of potentially toxic petroleum-derived DOM. The five types were a liquid-liquid extraction (LLE) with dichloromethane (DCM) and total petroleum hydrocarbons-diesel range (TPH_d) with DCM (EPA method 3510C), and three solid-phase extraction (SPE) stationary phases that are routinely used for extraction of polar analytes from water. For the LLE and TPH_d, that is selective for nonpolar compounds, the extraction efficiency of petroleum-derived DOM decreased downgradient as the petroleum-derived DOM becomes increasingly polar due to biodegradation. In contrast, the average extraction efficiency by the SPE methods was greater than 65% across the gradient. The results showed that SPE is more efficient for extracting petroleum-derived DOM at hydrocarbon-contaminated sites. The use of a method with greater extraction efficiency for partially oxidized hydrocarbons may prove useful in determining relationships between their composition and structure and potential for risks to human health or the environment.     

The Impact of Vertical Resolution in the Explicit Numerical Forecasting of Radiation Fog: A Case Study

Numerical experiments are performed with a comprehensive one-dimensional boundary layer/fog model to assess the impact of vertical resolution on explicit model forecasts of an observed fog layer. Two simulations were performed, one using a very high resolution and another with a vertical grid typical of current high-resolution mesoscale models. Both simulations were initialized with the same profiles, derived from observations from a fog field experiment. Significant differences in the onset and evolution of fog were found. The results obtained with the high-resolution simulation are in overall better agreement with available observations. The cooling rate before the appearance of fog is better represented, while the evolution of the liquid water content within the fog layer is more realistic. Fog formation is delayed in the low resolution simulation, and the water content in the fog layer shows large-amplitude oscillations. These results show that the numerical representation of key thermo-dynamical processes occurring in fog layers is significantly altered by the use of a grid with reduced vertical resolution.     

Distributed Parallel Computing in Stochastic Modeling of Groundwater Systems

Stochastic modeling is a rapidly evolving, popular approach to the study of the uncertainty and heterogeneity of groundwater systems. However, the use of Monte Carlo‐type simulations to solve practical groundwater problems often encounters computational bottlenecks that hinder the acquisition of meaningful results. To improve the computational efficiency, a system that combines stochastic model generation with MODFLOW‐related programs and distributed parallel processing is investigated. The distributed computing framework, called the Java Parallel Processing Framework, is integrated into the system to allow the batch processing of stochastic models in distributed and parallel systems. As an example, the system is applied to the stochastic delineation of well capture zones in the Pinggu Basin in Beijing. Through the use of 50 processing threads on a cluster with 10 multicore nodes, the execution times of 500 realizations are reduced to 3% compared with those of a serial execution. Through this application, the system demonstrates its potential in solving difficult computational problems in practical stochastic modeling.     

Monitoring and Mathematical Modeling of Contaminated Ground-Water Plumes in Fluvial Environments

Ground-water monitoring to delineate a contaminant plume in fluvial hydrostratigraphic units often is uncertain. Fluvial deposits consist typically of interbedded layers of sands, silts and clays, with buried stream channel deposits of sands or gravels. The channel deposits are often interpreted erroneously to be discontinuous between test holes and in cross section due to their sinuosity. Erroneous conclusions pertaining to the areal continuity of these geometrically complex deposits are inevitable unless the investigator thoroughly understands the depositional environment(s). The hydraulic conductivity of buried stream channel deposits may be several orders of magnitude higher than the matrix materials in which they are enclosed. The higher hydraulic conductivity of buried stream channel deposits has potentially significant ramifications with respect to ground-water monitoring to delineate the geometry of a contaminant plume migrating through these deposits. Ground-water monitoring at uranium mill waste disposal sites located in fluvial environments began on a significant scale in about 1977. A uranium mill tailing disposal site located in such an environment in central Wyoming is among the first sites monitored. Thirty-seven monitor wells were constructed at the site to delineate a seepage plume originating from one of the tailing ponds. This case history illustrates the need for a detailed under—standing of the hydrostratigraphy at a waste disposal site in order to interpret the meaning of ground-water quality data effectively. Water quality data from monitor wells located on a hit or miss basis often are misleading. The hydrostratigraphic horizon from which a water quality sample is collected must be well defined before the sample analyses can be interpreted quantitatively.     

Field Treatment of MTBE‐Contaminated Groundwater Using Ozone/UV Oxidation

Methyl‐tertiary butyl ether (MTBE) is often found in groundwater as a result of gasoline spills and leaking underground storage tanks. An extrapolation of occurrence data in 2008 estimated at least one detection of MTBE in approximately 165 small and large public water systems serving 896,000 people nationally (United States Environmental Protection Agency [U.S. EPA] 2008). The objective of this collaborative field study was to evaluate a small groundwater treatment system to determine the effectiveness of ultraviolet (UV)/ozone treatment in removing MTBE from contaminated drinking water wells. A pilot‐scale advanced oxidation process (AOP) system was tested to evaluate the oxidation efficiency of MTBE and intermediates under field conditions. This system used ozone as an oxidizer in the presence of UV light at hydraulic retention times varying from 1 to 3 min. MTBE removal efficiencies approaching 97% were possible with this system, even with low retention times. The intermediate t‐butyl alcohol (TBA) was removed to a lesser extent (71%) under the same test conditions. The main intermediate formed in the oxidation process of the contaminated groundwater in these studies was acetone. The concentrations of the other anticipated intermediates t‐butyl formate (TBF), isopropyl alcohol (IPA), methyl acetate (MAc), and possible co‐occurring aromatics (BTEX) in the effluent were negligible.     

Reuse of Drinking Water Treatment Waste for Remediation of Heavy Metal Contaminated Groundwater

Lime softening produces an estimated 10,000 metric tons of dry drinking water treatment wastes (DWTW) per year, costing an estimated one billion dollars annually for disposal worldwide. Lime softening wastes have been investigated for reuse as internal curing agents or supplementary cementitious materials in concrete as well as a high-capacity sorbent for heavy metal removal. Lead, cadmium, and zinc are common heavy metals in groundwater contaminated by mine tailings. Cement-based filter media (CBFM) are a novel material-class for heavy metal remediation in groundwater. This study investigated the incorporation of DWTW as a recycled, low-cost additive to CBFM for the removal of lead, cadmium, and zinc. Jar testing at three different metal concentrations and breakthrough column testing using synthetic groundwater were performed to measure removal capacity and reaction kinetics. Jar testing results show as DWTW content increases at low concentrations, removal approaches 100% but at high metal concentrations removal decreases due to saturation or exhaustion of the removal mechanisms. Removal occurs through the formation of metal carbonate precipitates, surface sorption, and ion exchange with calcium according to the preferential series Pb⁺² > Zn⁺² > Cd²⁺. Removal kinetics were also measured through column testing and exceeded estimated calculations derived from batch jar testing isotherms due to the large formation of oolitic metal carbonates. Lead, cadmium, and zinc was concentrated in the column precipitates from 0.29, 0.23, and 20.0 μg/g in the influent solution to approximately 200, 130, 14,000 μg/g in the reacted DWTW-CBFM. The control and DWTW-CBFM columns had statically similar removal for zinc and lead. In the DWTW-CBFM, cadmium had decreased removal of approximately 25% due to proportionately decreased hydroxide content from cement replacement with 25% DWTW. This study shows the potential for DWTW as an enhancement to CBFM, thereby valorizing an otherwise waste material. Furthermore, the concentrative abilities of CBFM through precipitate and oolitic mineral formation could provide a minable waste product and close the waste-product cycle for DWTW.     

Modeling the Factors Impacting Pesticide Concentrations in Groundwater Wells

This study examines the effect of pumping, hydrogeology, and pesticide characteristics on pesticide concentrations in production wells using a reactive transport model in two conceptual hydrogeologic systems; a layered aquifer with and without a stream present. The pumping rate can significantly affect the pesticide breakthrough time and maximum concentration at the well. The effect of the pumping rate on the pesticide concentration depends on the hydrogeology of the aquifer; in a layered aquifer, a high pumping rate resulted in a considerably different breakthrough than a low pumping rate, while in an aquifer with a stream the effect of the pumping rate was insignificant. Pesticide application history and properties have also a great impact on the effect of the pumping rate on the concentration at the well. The findings of the study show that variable pumping rates can generate temporal variability in the concentration at the well, which helps understanding the results of groundwater monitoring programs. The results are used to provide guidance on the design of pumping and regulatory changes for the long‐term supply of safe groundwater. The fate of selected pesticides is examined, for example, if the application of bentazone in a region with a layered aquifer stops today, the concentration at the well can continue to increase for 20 years if a low pumping rate is applied. This study concludes that because of the rapid response of the pesticide concentration at the drinking water well due to changes in pumping, wellhead management is important for managing pesticide concentrations.     

Accelerating Streamline Tracking in Groundwater Flow Modeling on GPUs

Streamline simulation in groundwater flow modeling is a time-consuming process when a large number of streamlines are analyzed. We develop a parallelization method on graphics processing units (GPUs) for the semi-analytical particle tracking algorithm developed by Pollock (1988). Compute Unified Device Architecture was used to implement the parallel method. Forward and backward tracking of a streamline is handled by an individual thread. A GPU includes a grid of blocks where a block handles 32 threads. We use multi-GPUs to accelerate streamline tracking in a flow model with millions of particles. The method was examined to simulate streamlines for identifying three-dimensional (3D) flow systems in a Tóthian basin. The speedup exceeds 1000 when 8 NVIDIA GPUs are used to simulate 5 million or more streamlines.     

Ecohydrologic Process Modeling of Mountain Block Groundwater Recharge

Regional mountain block recharge (MBR) is a key component of alluvial basin aquifer systems typical of the western United States. Yet neither water scientists nor resource managers have a commonly available and reasonably invoked quantitative method to constrain MBR rates. Recent advances in landscape-scale ecohydrologic process modeling offer the possibility that meteorological data and land surface physical and vegetative conditions can be used to generate estimates of MBR. A water balance was generated for a temperate 24,600-ha mountain watershed, elevation 1565 to 3207 m, using the ecosystem process model Biome-BGC (BioGeochemical Cycles) ( Running and Hunt 1993 ). Input data included remotely sensed landscape information and climate data generated with the Mountain Climate Simulator (MT-CLIM) ( Running et al. 1987 ). Estimated mean annual MBR flux into the crystalline bedrock terrain is 99,000 m³/d, or approximately 19% of annual precipitation for the 2003 water year. Controls on MBR predictions include evapotranspiration (radiation limited in wet years and moisture limited in dry years), soil properties, vegetative ecotones (significant at lower elevations), and snowmelt (dominant recharge process). The ecohydrologic model is also used to investigate how climatic and vegetative controls influence recharge dynamics within three elevation zones. The ecohydrologic model proves useful for investigating controls on recharge to mountain blocks as a function of climate and vegetation. Future efforts will need to investigate the uncertainty in the modeled water balance by incorporating an advanced understanding of mountain recharge processes, an ability to simulate those processes at varying scales, and independent approaches to calibrating MBR estimates.     

Artificial Sweeteners Identify Spatial Patterns of Historic Landfill Contaminated Groundwater Discharge in an Urban Stream

Leachate-contaminated groundwater from historical municipal landfills, typically lacking engineered liners and leachate collection systems, poses a threat to nearby urban streams, particularly to benthic ecosystems. Effective monitoring and assessment of such sites requires understanding of the spatial patterns (i.e., two-dimensional footprint) of contaminated groundwater discharge and associated controlling factors. However, discharges from groundwater contaminated by modern wastewater can complicate site assessments. The objectives of this study were to (1) demonstrate the use of artificial sweeteners (AS): saccharin (SAC), cyclamate (CYC), acesulfame (ACE), and sucralose (SUC), to distinguish groundwater discharge areas influenced by historic landfill leachate (elevated SAC and sometimes CYC; low ACE and SUC concentrations) from those influenced by wastewater (high ACE and SUC concentrations), and (2) investigate contaminant discharge patterns for two gaining urban stream reaches adjacent historic landfills at base flows. Contaminant discharge patterns revealed by the AS were strongly controlled by hyporheic flow (low AS concentrations), particularly for the straight reach, and stream sinuosity, particularly for the meandering reach. These patterns were different and the contaminant footprint coverage (<25% of streambed area) much less than most past studies (typically >50% coverage), likely due to the homogeneous streambed-aquifer conditions and shallow, narrow landfill plume in this setting.     

Reinterpreting the Importance of Oxygen‐Based Biodegradation in Chloroethene‐Contaminated Groundwater

Chlororespiration is common in shallow aquifer systems under conditions nominally identified as anoxic. Consequently, chlororespiration is a key component of remediation at many chloroethene‐contaminated sites. In some instances, limited accumulation of reductive dechlorination daughter products is interpreted as evidence that natural attenuation is not adequate for site remediation. This conclusion is justified when evidence for parent compound (tetrachloroethene, PCE, or trichloroethene, TCE) degradation is lacking. For many chloroethene‐contaminated shallow aquifer systems, however, nonconservative losses of the parent compounds are clear but the mass balance between parent compound attenuation and accumulation of reductive dechlorination daughter products is incomplete. Incomplete mass balance indicates a failure to account for important contaminant attenuation mechanisms and is consistent with contaminant degradation to nondiagnostic mineralization products like CO₂. While anoxic mineralization of chloroethene compounds has been proposed previously, recent results suggest that oxygen‐based mineralization of chloroethenes also can be significant at dissolved oxygen concentrations below the currently accepted field standard for nominally anoxic conditions. Thus, reassessment of the role and potential importance of low concentrations of oxygen in chloroethene biodegradation are needed, because mischaracterization of operant biodegradation processes can lead to expensive and ineffective remedial actions. A modified interpretive framework is provided for assessing the potential for chloroethene biodegradation under different redox conditions and the probable role of oxygen in chloroethene biodegradation.