Evaluation of an Impulse-Response Emulator for Groundwater Contaminant Transport Modeling.

There is a significant need to develop decision support tools capable of delivering accurate representations of environmental conditions, such as ground and surface water solute concentrations, in a timely and computationally efficient manner. Such tools can be leveraged to assess a large number of potential management strategies for mitigating non-point source pollutants. Here, we assess the effectiveness of the impulse-response emulation approach to approximate process-based groundwater model estimates of solute transport from MODFLOW and MT3D over a wide range of model inputs and parameters, with the goal of assessing where in parameter space the assumptions underlying this emulation approach are valid. The impulse-response emulator was developed using the sensitivity analysis utilities in the PEST++ software suite and is capable of approximating MODFLOW/MT3D estimates of solute transport over a large portion of the parameter space tested, except in cases where the Courant number is above 0.5. Across all runs tested, the highest percent errors were at the plume fronts. These results suggest that the impulse-response approach may be suitable for emulation of solute transport models for a wide range of cases, except when high-resolution outputs are needed, or when very low concentrations at plume edges are of particular interest.

SFRmaker and Linesink-Maker: Rapid Construction of Streamflow Routing Networks from Hydrography Data.

Groundwater models have evolved to encompass more aspects of the water cycle, but the incorporation of realistic boundary conditions representing surface water remains time-consuming and error-prone. We present two Python packages that robustly automate this process using readily available hydrography data as the primary input. SFRmaker creates input for the MODFLOW SFR package, while Linesink-maker creates linesink string input for the GFLOW analytic element program. These programs can reduce weeks or even months of manual effort to a few minutes of execution time, and carry the added advantages of reduced potential for error, improved reproducibility and facilitation of step-wise modeling through reduced dependency on a particular conceptual model or discretization. Two real-world examples at the county to multi-state scales are presented.

A Semi-Structured MODFLOW-USG Model to Evaluate Local Water Sources to Wells for Decision Support.

In order to better represent the configuration of the stream network and simulate local groundwater-surface water interactions, a version of MODFLOW with refined spacing in the topmost layer was applied to a Lake Michigan Basin (LMB) regional groundwater-flow model developed by the U.S. Geological. Regional MODFLOW models commonly use coarse grids over large areas; this coarse spacing precludes model application to local management issues (e.g., surface-water depletion by wells) without recourse to labor-intensive inset models. Implementation of an unstructured formulation within the MODFLOW framework (MODFLOW-USG) allows application of regional models to address local problems. A "semi-structured" approach (uniform lateral spacing within layers, different lateral spacing among layers) was tested using the LMB regional model. The parent 20-layer model with uniform 5000-foot (1524-m) lateral spacing was converted to 4 layers with 500-foot (152-m) spacing in the top glacial (Quaternary) layer, where surface water features are located, overlying coarser resolution layers representing deeper deposits. This semi-structured version of the LMB model reproduces regional flow conditions, whereas the finer resolution in the top layer improves the accuracy of the simulated response of surface water to shallow wells. One application of the semi-structured LMB model is to provide statistical measures of the correlation between modeled inputs and the simulated amount of water that wells derive from local surface water. The relations identified in this paper serve as the basis for metamodels to predict (with uncertainty) surface-water depletion in response to shallow pumping within and potentially beyond the modeled area, see Fienen et al. (2015a).

Groundwater availability as constrained by hydrogeology and environmental flows.

Groundwater pumping from aquifers in hydraulic connection with nearby streams has the potential to cause adverse impacts by decreasing flows to levels below those necessary to maintain aquatic ecosystems. The recent passage of the Great Lakes-St. Lawrence River Basin Water Resources Compact has brought attention to this issue in the Great Lakes region. In particular, the legislation requires the Great Lakes states to enact measures for limiting water withdrawals that can cause adverse ecosystem impacts. This study explores how both hydrogeologic and environmental flow limitations may constrain groundwater availability in the Great Lakes Basin. A methodology for calculating maximum allowable pumping rates is presented. Groundwater availability across the basin may be constrained by a combination of hydrogeologic yield and environmental flow limitations varying over both local and regional scales. The results are sensitive to factors such as pumping time, regional and local hydrogeology, streambed conductance, and streamflow depletion limits. Understanding how these restrictions constrain groundwater usage and which hydrogeologic characteristics and spatial variables have the most influence on potential streamflow depletions has important water resources policy and management implications.

A new capture fraction method to map how pumpage affects surface water flow.

All groundwater pumped is balanced by removal of water somewhere, initially from storage in the aquifer and later from capture in the form of increase in recharge and decrease in discharge. Capture that results in a loss of water in streams, rivers, and wetlands now is a concern in many parts of the United States. Hydrologists commonly use analytical and numerical approaches to study temporal variations in sources of water to wells for select points of interest. Much can be learned about coupled surface/groundwater systems, however, by looking at the spatial distribution of theoretical capture for select times of interest. Development of maps of capture requires (1) a reasonably well-constructed transient or steady state model of an aquifer with head-dependent flow boundaries representing surface water features or evapotranspiration and (2) an automated procedure to run the model repeatedly and extract results, each time with a well in a different location. This paper presents new methods for simulating and mapping capture using three-dimensional groundwater flow models and presents examples from Arizona, Oregon, and Michigan.

Linking MODFLOW with an agent-based land-use model to support decision making.

The U.S. Geological Survey numerical groundwater flow model, MODFLOW, was integrated with an agent-based land-use model to yield a simulator for environmental planning studies. Ultimately, this integrated simulator will be used as a means to organize information, illustrate potential system responses, and facilitate communication within a participatory modeling framework. Initial results show the potential system response to different zoning policy scenarios in terms of the spatial patterns of development, which is referred to as urban form, and consequent impacts on groundwater levels. These results illustrate how the integrated simulator is capable of representing the complexity of the system. From a groundwater modeling perspective, the most important aspect of the integration is that the simulator generates stresses on the groundwater system within the simulation in contrast to the traditional approach that requires the user to specify the stresses through time.

Directed site exploration for permeable reactive barrier design.

Permeable reactive barriers (PRBs) are being employed for in situ site remediation of groundwater that is typically flowing under natural gradients. Site characterization is of critical importance to the success of a PRB. A design-specific site exploration approach called quantitatively directed exploration (QDE) is presented. The QDE approach employs three spatially related matrices: (1) covariance of input parameters, (2) sensitivity of model outputs, and (3) covariance of model outputs to identify the most important location to explore based on a specific design. Sampling at the location that most reduces overall site uncertainty produces a higher probability of success of a particular design. The QDE approach is demonstrated on the Kansas City Plant, Kansas City, MO, a case study where a PRB was installed and failed. It is shown that additional quantitatively directed site exploration during the design phase could have prevented the remedial failure that was caused by missing a geologic body having high hydraulic conductivity at the south end of the barrier. The most contributing input parameter approach using head uncertainty clearly indicated where the next sampling should be made toward the high hydraulic conductivity zone. This case study demonstrates the need to include the specific design as well as site characterization uncertainty when choosing the sampling locations.

MODFLOW 2000 head uncertainty, a first-order second moment method.

A computationally efficient method to estimate the variance and covariance in piezometric head results computed through MODFLOW 2000 using a first-order second moment (FOSM) approach is presented. This methodology employs a first-order Taylor series expansion to combine model sensitivity with uncertainty in geologic data. MODFLOW 2000 is used to calculate both the ground water head and the sensitivity of head to changes in input data. From a limited number of samples, geologic data are extrapolated and their associated uncertainties are computed through a conditional probability calculation. Combining the spatially related sensitivity and input uncertainty produces the variance-covariance matrix, the diagonal of which is used to yield the standard deviation in MODFLOW 2000 head. The variance in piezometric head can be used for calibrating the model, estimating confidence intervals, directing exploration, and evaluating the reliability of a design. A case study illustrates the approach, where aquifer transmissivity is the spatially related uncertain geologic input data. The FOSM methodology is shown to be applicable for calculating output uncertainty for (1) spatially related input and output data, and (2) multiple input parameters (transmissivity and recharge).

Axi-symmetric simulation of soil vapor extraction influenced by soil fracturing.

Fracturing, either pneumatic or hydraulic, is a method to improve the performance of soil vapor extraction (SVE) in relatively low permeability soils (< 10(-5) cm/s). A two-dimensional model is presented to simulate trichloroethylene (TCE) soil vapor extraction modified by fracturing. Flow and transport is modeled using mobile macropore and micropore networks, which also have been identified in the literature as dual porosity, dual permeability, or heterogeneous flow models. In this model, fluids can flow in both the macropore and micropore networks. This represents a more general model compared to immobile micropore, mobile macropore models presented thus far in the literature for vapor flow and transport in two dimensions. The model considers pressure- and concentration-driven exchange between the macropore and micropore networks, concentration-driven exchange between the gas and sorbed phases within each network, and equilibrium exchange between the gas and water and a sorbed phase within each network. The parameters employed in an example simulation are based on field measurements made at a fractured site. Considered in the simulations were the influence of the volume percentage of fractures, the length of fractures, the relative location of the water table, and the influence of pulsed pumping. For these simulations, internetwork concentration-driven exchange most significantly affected mass removal. The volume percentage of fractures more significantly influence flow and mass removal than the length of fractures. The depth of the water table below the contamination plume only significantly influenced flow and mass removal when the water table was within 60 cm of the bottom of the contaminated soil in the vadose zone for the parameters considered in this study. Pulsed pumping was not found to increase the amount of mass removed in this study.