Risk-Based Wellhead Protection Decision Support: A Repeatable Workflow Approach.

Environmental water management often benefits from a risk-based approach where information on the area of interest is characterized, assembled, and incorporated into a decision model considering uncertainty. This includes prior information from literature, field measurements, professional interpretation, and data assimilation resulting in a decision tool with a posterior uncertainty assessment accounting for prior understanding and what is learned through model development and data assimilation. Model construction and data assimilation are time consuming and prone to errors, which motivates a repeatable workflow where revisions resulting from new interpretations or discovery of errors can be addressed and the analyses repeated efficiently and rigorously. In this work, motivated by the real world application of delineating risk-based (probabilistic) sources of water to supply wells in a humid temperate climate, a scripted workflow was generated for groundwater model construction, data assimilation, particle-tracking and post-processing. The workflow leverages existing datasets describing hydrogeology, hydrography, water use, recharge, and lateral boundaries. These specific data are available in the United States but the tools can be applied to similar datasets worldwide. The workflow builds the model, performs ensemble-based history matching, and uses a posterior Monte Carlo approach to provide probabilistic capture zones describing source water to wells in a risk-based framework. The water managers can then select areas of varying levels of protection based on their tolerance for risk of potential wrongness of the underlying models. All the tools in this workflow are open-source and free, which facilitates testing of this repeatable and transparent approach to other environmental problems.

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.

Risk-Based Decision-Support Groundwater Modeling for the Lower San Antonio River Basin, Texas, USA.

A numerical surface-water/groundwater model was developed for the lower San Antonio River Basin to evaluate the responses of low base flows and groundwater levels within the basin under conditions of reduced recharge and increased groundwater withdrawals. Batch data assimilation through history matching used a simulation of historical conditions (2006-2013); this process included history-matching to groundwater levels and base-flow estimates at several gages, and was completed in a high-dimensional (highly parameterized) framework. The model was developed in an uncertainty framework such that parameters, observations, and scenarios of interest are envisioned stochastically as distributions of potential values. Results indicate that groundwater contributions to surface water during periods of low flow may be reduced from 6% to 25% with a corresponding 25% reduction in recharge and a 25% increase in groundwater pumping over an 8-year planning period. Furthermore, results indicate groundwater-level reductions in some hydrostratigraphic units are more likely than in other hydrostratigraphic units over an 8-year period under drought conditions with the higher groundwater withdrawal scenario.

Active thermal tracer tests for improved hydrostratigraphic characterization.

Subsurface heterogeneity in hydraulic properties and processes is a fundamental challenge in hydrogeology. We have developed an improved method of borehole dilution testing for hydrostratigraphic characterization, in which distributed temperature sensing (DTS) is used to monitor advective heat movement. DTS offers many advantages over conventional technologies including response times in the order of seconds rather than minutes, the ability to profile temperature synoptically in a well without disturbing the fluid column, sensitivity to a wider range of flow rates than conventional spinner and heat pulse flow meters, and the ease of interpretation. Open-well thermal dilution tests in two multiaquifer wells near Madison, Wisconsin, provided detailed information on the borehole flow regimes, including flow rates and the locations of inflows from both fractures and porous media. The results led to an enhanced understanding of flow in a hydrostratigraphic unit previously conceptualized as homogenous and isotropic.