An Agricultural Package for MODFLOW 6 Using the Application Programming Interface

An agricultural water use package has been developed for MODFLOW 6 using the MODFLOW Application Programming Interface (API). The MODFLOW API Agricultural Water Use Package (API-Ag) was based on the approach to simulate irrigation demand in the MODFLOW-NWT and GSFLOW Agricultural Water Use (AG) Package. The API-Ag Package differs from the previous approach by implementing new features and support for additional irrigation providers. New features include representation of deficit and over-irrigation, Multi-Aquifer Well and Lake Package irrigation providers, and support for structured, vertex, and unstructured grid models. Three example problems are presented that demonstrate how the API-Ag Package improves representation of highly managed systems and are further used to validate the irrigation demand and delivery formulations. Irrigation volumes simulated in the three example problems show excellent agreement with the MODFLOW-NWT AG Package.     

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).     

A Simple Method for Simulating Groundwater Interactions with Fens to Forecast Development Effects

Protection of fens–wetlands dependent on groundwater discharge–requires characterization of groundwater sources and stresses. Because instrumentation and numerical modeling of fens is labor intensive, easy-to-apply methods that model fen distribution and their vulnerability to development are desirable. Here we demonstrate that fen areas can be simulated using existing steady-state MODFLOW models when the unsaturated zone flow (UZF) package is included. In cells where the water table is near land surface, the UZF package calculates a head difference and scaled conductance at these “seepage drain” cells to generate average rates of vertical seepage to the land. This formulation, which represents an alternative to blanketing the MODFLOW domain with drains, requires very little input from the user because unsaturated flow-routing is inactive and results are primarily driven by easily obtained topographic information. Like the drain approach, it has the advantage that the distribution of seepage areas is not predetermined by the modeler, but rather emerges from simulated heads. Beyond the drain approach, it takes account of intracell land surface variation to explicitly quantify multiple surficial flows corresponding to infiltration, rejected recharge, recharge and land-surface seepage. Application of the method to a basin in southeastern Wisconsin demonstrates how it can be used as a decision-support tool to first, reproduce fen distribution and, second, forecast drawdown and reduced seepage at fens in response to shallow pumping.     

MODFLOW as a Configurable Multi-Model Hydrologic Simulator

MODFLOW 6 is the latest in a line of six “core” versions of MODFLOW released by the U.S. Geological Survey. The MODFLOW 6 architecture supports incorporation of additional hydrologic processes, in addition to groundwater flow, and allows interaction between processes. The architecture supports multiple model instances and multiple types of models within a single simulation, a flexible approach to formulating and solving the equations that represent hydrologic processes, and recent advances in interoperability, which allow MODFLOW to be accessed and controlled by external programs. The present version of MODFLOW 6 consolidates popular capabilities available in MODFLOW variants, such as the unstructured grid support in MODFLOW-USG, the Newton-Raphson formulation in MODFLOW-NWT, and the support for partitioned stress boundaries in MODFLOW-CDSS. The flexible multi-model capability allows users to configure MODFLOW 6 simulations to represent the local-grid refinement (LGR) capabilities available in MODFLOW-LGR, the multi-species transport capabilities in MT3DMS, and the coupled variable-density capabilities available in SEAWAT. This paper provides a new, holistic and integrated overview of simulation capabilities made possible by the MODFLOW 6 architecture, and describes how ongoing and future development can take advantage of the program architecture to integrate new capabilities in a way that is minimally invasive and automatically compatible with the existing MODFLOW 6 code.     

Irrigation effects in the northern lake states: Wisconsin central sands revisited

Irrigated agriculture has expanded greatly in the water-rich U.S. northern lake states during the past half century. Source water there is usually obtained from glacial aquifers strongly connected to surface waters, so irrigation has a potential to locally decrease base flows in streams and water levels in aquifers, lakes, and wetlands. During the nascent phase of the irrigation expansion, water availability was explored in works of some fame in the Wisconsin central sands by Weeks et al. (1965) on the Little Plover River and Weeks and Stangland (1971) on "headwater area" streams and lakes. Four decades later, and after irrigation has grown to a dominant landscape presence, we revisited irrigation effects on central sands hydrology. Irrigation effects have been substantial, on average decreasing base flows by a third or more in many stream headwaters and diminishing water levels by more than a meter in places. This explains why some surface waters have become flow and stage impaired, sometimes to the point of drying, with attendant losses of aquatic ecosystems. Irrigation exerts its effects by increasing evapotranspiration by an estimated 45 to 142 mm/year compared with pre-irrigated land cover. We conclude that irrigation water availability in the northern lake states and other regions with strong groundwater-surface water connections is tied to concerns for surface water health, requiring a focus on managing the upper few meters of aquifers on which surface waters depend rather than the depletability of an aquifer.     

Field Test of a Hybrid Finite‐Difference and Analytic Element Regional Model

Regional finite‐difference models often have cell sizes that are too large to sufficiently model well‐stream interactions. Here, a steady‐state hybrid model is applied whereby the upper layer or layers of a coarse MODFLOW model are replaced by the analytic element model GFLOW, which represents surface waters and wells as line and point sinks. The two models are coupled by transferring cell‐by‐cell leakage obtained from the original MODFLOW model to the bottom of the GFLOW model. A real‐world test of the hybrid model approach is applied on a subdomain of an existing model of the Lake Michigan Basin. The original (coarse) MODFLOW model consists of six layers, the top four of which are aggregated into GFLOW as a single layer, while the bottom two layers remain part of MODFLOW in the hybrid model. The hybrid model and a refined “benchmark” MODFLOW model simulate similar baseflows. The hybrid and benchmark models also simulate similar baseflow reductions due to nearby pumping when the well is located within the layers represented by GFLOW. However, the benchmark model requires refinement of the model grid in the local area of interest, while the hybrid approach uses a gridless top layer and is thus unaffected by grid discretization errors. The hybrid approach is well suited to facilitate cost‐effective retrofitting of existing coarse grid MODFLOW models commonly used for regional studies because it leverages the strengths of both finite‐difference and analytic element methods for predictions in mildly heterogeneous systems that can be simulated with steady‐state conditions.     

A simulation and prediction of agricultural irrigation on groundwater in well irrigation area of the piedmont of Mt. Taihang,North China

Declining groundwater levels caused by irrigation is the main problem for agricultural development in Northern China. Due to both economics and increased population, surface water has become almost non‐existent and groundwater is the only water resource left. Currently the groundwater is declining at a rate between 50 and 100 cm per year. Sustainable development in Northern China requires effective management of the groundwater resources. In this study, the effect of future irrigation patterns on the decline of the groundwater table is examined with the aid of MODFLOW. MODFLOW was calibrated for five observation wells in the county. The calibrated model fitted the observed data well over a 7‐year period. The simulated results showed that the groundwater decline would be decreased, and perhaps halted, by decreasing the use of irrigation. Copyright © 2005 John Wiley & Sons, Ltd.     

A MODFLOW Infiltration Device Package for Simulating Storm Water Infiltration

This article describes a MODFLOW Infiltration Device (INFD) Package that can simulate infiltration devices and their two‐way interaction with groundwater. The INFD Package relies on a water balance including inflow of storm water, leakage‐like seepage through the device faces, overflow, and change in storage. The water balance for the device can be simulated in multiple INFD time steps within a single MODFLOW time step, and infiltration from the device can be routed through the unsaturated zone to the groundwater table. A benchmark test shows that the INFD Package's analytical solution for stage computes exact results for transient behavior. To achieve similar accuracy by the numerical solution of the MODFLOW Surface‐Water Routing (SWR1) Process requires many small time steps. Furthermore, the INFD Package includes an improved representation of flow through the INFD sides that results in lower infiltration rates than simulated by SWR1. The INFD Package is also demonstrated in a transient simulation of a hypothetical catchment where two devices interact differently with groundwater. This simulation demonstrates that device and groundwater interaction depends on the thickness of the unsaturated zone because a shallow groundwater table (a likely result from storm water infiltration itself) may occupy retention volume, whereas a thick unsaturated zone may cause a phase shift and a change of amplitude in groundwater table response to a change of infiltration. We thus find that the INFD Package accommodates the simulation of infiltration devices and groundwater in an integrated manner on small as well as large spatial and temporal scales.     

Adjoint simulation of stream depletion due to aquifer pumping

If an aquifer is hydraulically connected to an adjacent stream, a pumping well operating in the aquifer will draw some water from aquifer storage and some water from the stream, causing stream depletion. Several analytical, semi-analytical, and numerical approaches have been developed to estimate stream depletion due to pumping. These approaches are effective if the well location is known. If a new well is to be installed, it may be desirable to install the well at a location where stream depletion is minimal. If several possible locations are considered for the location of a new well, stream depletion would have to be estimated for all possible well locations, which can be computationally inefficient. The adjoint approach for estimating stream depletion is a more efficient alternative because with one simulation of the adjoint model, stream depletion can be estimated for pumping at a well at any location. We derive the adjoint equations for a coupled system with a confined aquifer, an overlying unconfined aquifer, and a river that is hydraulically connected to the unconfined aquifer. We assume that the stage in the river is known, and is independent of the stream depletion, consistent with the assumptions of the MODFLOW river package. We describe how the adjoint equations can be solved using MODFLOW. In an illustrative example, we show that for this scenario, the adjoint approach is as accurate as standard forward numerical simulation methods, and requires substantially less computational effort.     

Describing alpine lake influence on stream network temperatures: A statistical modelling approach

Systematic variations in atmospheric heat exchange, surface residence time, and groundwater influx across montane stream networks commonly produce an increasing stream temperature trend with decreasing elevation. However, complex stream temperature profiles that differ from this common longitudinal trend also exist, suggesting that stream temperatures may be influenced by complex interactions among hydrologic and atmospheric processes. Lakes within stream networks form one potential source of temperature profile complexity due to the spatially variable contribution of lake-sourced water to stream flow. We investigated temperature profile complexity in a multi-season stream temperature dataset collected across a montane stream network containing many alpine lakes. This investigation was performed by making comparisons between multiple statistical models that used different combinations of stream and lake characteristics to represent specific hypotheses for the controls on stream temperature. The compared models included a set of models which used a topographically derived estimate of the hydrologic influence of lakes to separate and quantify the effects of stream elevation and lake source-water contributions to longitudinal stream temperature patterns. This source-water mixing model provided a parsimonious explanation for complex stream-network temperature patterns in the summer and autumn, and this approach may be further applicable to other systems where stream temperatures are influenced by multiple water sources. Simpler models that discounted lake effects were more optimal during the winter and spring, suggesting that complex patterns in stream temperature profiles may emerge and subside temporally, across seasons, in response to diversity of water temperatures from different sources.     

A Generic Closed-Loop Contaminant Treatment System Package for MODFLOW 6

We present a contaminant treatment system (CTS) package for MODFLOW 6 that facilitates the simulation of pump-and-treat systems for groundwater remediation. Using the “nonintrusive” MODFLOW 6 application programming interface (API) capability, the CTS package can balance flows between extraction and injection wells within the outer flow solution loop and applies blended concentration/mass treatment efficiency within the outer transport solution loop. The former can be important when the requested extraction rate cannot be satisfied by the current simulated groundwater system conditions, while the latter can be important for simulating incomplete/imperfect treatment schemes. Furthermore, the CTS package allows users to temporally vary all aspects of a simulated CTS system, including the configuration and location of injection and extraction wells, and the CTS efficiency. This flexibility combined with the API-based implementation provide a generic and general CTS package that can be applied across the wide range of MODFLOW 6 simulation options and that evolves in step with MODFLOW 6 code modifications and advancements without needing to update the CTS package itself.     

Embedding isochronous cells overland flow module into MODFLOW

This study integrates a simple overland flow module (isochronous cells model) with the river module of MODFLOW such that temporal and spatial interactions between stream flow and groundwater can be simulated using net rainfall data of a watershed. The isochronous cells model is an efficient travel time runoff approach based on geographic information system (GIS) that considers both spatial and temporal variations of net rainfall through hill slope of the watershed. This overland module is easily coupled with MODFLOW river routing module. Specifically, the stream flow from the isochronous cells model is directly assigned to both sides of river cells of the MODFLOW model. Such an integration of MODFLOW and isochronous cells model is especially useful in watersheds where river flow data are limited. The feasibility of this integrated model was demonstrated using a case study in the middle and downstream regions of the Yitong River watershed, China. Copyright © 2012 John Wiley & Sons, Ltd.     

Introducing environmental thresholds into water withdrawal management of mountain streams in the Kura River basin,Azerbaijan

Abstract

The study aims to set and implement environmentally relevant limits for the exploitation of mountain streams in the Kura River basin of Azerbaijan. Such streams represent the preferred spawning grounds for valuable sturgeon of the Caspian Sea, but experience continuously increasing exploitation in the form of water withdrawals for industry and irrigation. Since no detailed environmental flow assessments have been conducted on any of the Kura basin streams, an interim approach is suggested based on minimum flow, referred to as “base environmental minimum”. The latter may be estimated from the unregulated parts of observed or simulated daily flow records. Environmental flow requirements for individual months of an individual year may be calculated using correction factors related to monthly rainfall. Simple relationships are suggested for base environmental flow estimation at ungauged sites, and the implications of river pollution for monthly environmental requirements are examined. Further, definition of environmentally critical periods in a stream is proposed based on a ratio of observed to “environmental” flow as an indicator of environmental stress. It is illustrated that the conjunctive use of several closely located streams for water supply may significantly reduce the duration of, or completely eliminate, environmentally critical periods. The idea of environmentally acceptable areal water withdrawal is formulated, so that the overall approach may be applied for environmentally sustainable water withdrawal management in other small streams.     

A modeling approach for assessing the effect of multiple alpine lakes in sequence on nutrient transport

The effects of a single lake on downstream water chemistry may be compounded by the presence of additional lakes within the watershed, augmenting or negating the effects of the first lake. Multiple, linked lakes are a common feature of many watersheds and these resemble reactors in series often studied in engineering. The effects of multiple lakes in series on nutrient transport are largely unexplored. We populated and calibrated a simple lake model to investigate the role of a sub-alpine lake (Bull Trout Lake (BTL), Rocky Mountains, USA) on the transport of the macronutrients during the summer of 2008. Further, we developed a sequential model in which four identical lakes (copies of the BTL model) were connected in series. All lakes in the sequence retarded the flux of nutrients, thus slowing their transport downstream. The first lake in the sequence dramatically altered stream water chemistry and served as a sink for C and P and a source of N, while additional lakes downstream became sources of C, N and P. Although additional downstream lakes resulted in important changes to water chemistry and nutrient transport, the nature of the changes were similar from Lakes 2 to 4 and the magnitude of the changes diminished with distance downstream. Our lake model served as an effective tool for assessing the nutrient budget of the lake and the hypothetical effect of multiple lakes in sequence in a landscape limnology framework.     

A Hybrid Finite-Difference and Analytic Element Groundwater Model

Regional finite-difference models tend to have large cell sizes, often on the order of 1–2 km on a side. Although the regional flow patterns in deeper formations may be adequately represented by such a model, the intricate surface water and groundwater interactions in the shallower layers are not. Several stream reaches and nearby wells may occur in a single cell, precluding any meaningful modeling of the surface water and groundwater interactions between the individual features. We propose to replace the upper MODFLOW layer or layers, in which the surface water and groundwater interactions occur, by an analytic element model (GFLOW) that does not employ a model grid; instead, it represents wells and surface waters directly by the use of point-sinks and line-sinks. For many practical cases it suffices to provide GFLOW with the vertical leakage rates calculated in the original coarse MODFLOW model in order to obtain a good representation of surface water and groundwater interactions. However, when the combined transmissivities in the deeper (MODFLOW) layers dominate, the accuracy of the GFLOW solution diminishes. For those cases, an iterative coupling procedure, whereby the leakages between the GFLOW and MODFLOW model are updated, appreciably improves the overall solution, albeit at considerable computational cost. The coupled GFLOW–MODFLOW model is applicable to relatively large areas, in many cases to the entire model domain, thus forming an attractive alternative to local grid refinement or inset models.     

Ground water flow in a desert basin: challenges of simulating transport of dissolved chromium

A large chromium plume that evolved from chromium releases in a valley near the Mojave River was studied to understand the processes controlling fate and migration of chromium in ground water and used as a tracer to study the dynamics of a basin and range ground water system. The valley that was studied is naturally arid with high evapotranspiration such that essentially no precipitation infiltrates to the water table. The dominant natural hydrogeologic processes are recharge to the ground water system from the Mojave River during the infrequent episodes when there is flow in the river, and ground water flow toward a playa lake where the ground water evaporates. Agricultural pumping in the valley from the mid-1930s to the 1970s significantly altered ground water flow conditions by decreasing water levels in the valley by more than 20 m. This pumping declined significantly as a result of dewatering of the aquifer, and water levels have since recovered modestly. The ground water system was modeled using MODFLOW, and chromium transport was simulated using MT3D. Several innovative modifications were made to these modeling programs to simulate important processes in this ground water system. Modifications to MODFLOW include developing a new well package that estimates pumping rates from irrigation wells at each time step based on available drawdown. MT3D was modified to account for mass trapped above the water table when the water table declines beneath nonirrigated areas and to redistribute mass to the system when water levels rise.     

Incorporating Groundwater-Surface Water Interaction into River Management Models

Accurate representation of groundwater-surface water interactions is critical to modeling low river flows in the semi-arid southwestern United States. Although a number of groundwater-surface water models exist, they are seldom integrated with river operation/management models. A link between the object-oriented river and reservoir operations model, RiverWare, and the groundwater model, MODFLOW, was developed to incorporate groundwater-surface water interaction processes, such as river seepage/gains, riparian evapotranspiration, and irrigation return flows, into a rule-based water allocations model. An explicit approach is used in which the two models run in tandem, exchanging data once in each computational time step. Because the MODFLOW grid is typically at a finer resolution than RiverWare objects, the linked model employs spatial interpolation and summation for compatible communication of exchanged variables. The performance of the linked model is illustrated through two applications in the Middle Rio Grande Basin in New Mexico where overappropriation impacts endangered species habitats. In one application, the linked model results are compared with historical data; the other illustrates use of the linked model for determining management strategies needed to attain an in-stream flow target. The flows predicted by the linked model at gauge locations are reasonably accurate except during a few very low flow periods when discrepancies may be attributable to stream gaging uncertainties or inaccurate documentation of diversions. The linked model accounted for complex diversions, releases, groundwater pumpage, irrigation return flows, and seepage between the groundwater system and canals/drains to achieve a schedule of releases that satisfied the in-stream target flow.     

MODFLOW/MT3DMS-based simulation of variable-density ground water flow and transport

This paper presents an approach for coupling MODFLOW and MT3DMS for the simulation of variable-density ground water flow. MODFLOW routines were modified to solve a variable-density form of the ground water flow equation in which the density terms are calculated using an equation of state and the simulated MT3DMS solute concentrations. Changes to the MODFLOW and MT3DMS input files were kept to a minimum, and thus existing data files and data files created with most pre- and postprocessors can be used directly with the SEAWAT code. The approach was tested by simulating the Henry problem and two of the saltpool laboratory experiments (low- and high-density cases). For the Henry problem, the simulated results compared well with the steady-state semianalytic solution and also the transient isochlor movement as simulated by a finite-element model. For the saltpool problem, the simulated breakthrough curves compared better with the laboratory measurements for the low-density case than for the high-density case but showed good agreement with the measured salinity isosurfaces for both cases. Results from the test cases presented here indicate that the MODFLOW/MT3DMS approach provides accurate solutions for problems involving variable-density ground water flow and solute transport.     

Patterns of hydrologic connectivity in the McMurdo Dry Valleys,Antarctica: a synthesis of 20 years of hydrologic data

Streams in the McMurdo Dry Valleys (MDVs) of Antarctica moderate an important hydrologic and biogeochemical connection between upland alpine glaciers, valley‐bottom soils, and lowland closed‐basin lakes. Moreover, MDV streams are simple but dynamic systems ideal for studying interacting hydrologic and ecological dynamics. This work synthesizes 20 years of hydrologic data, collected as part of the MDVs Long‐Term Ecological Research project, to assess spatial and temporal dynamics of hydrologic connectivity between glaciers, streams, and lakes. Long‐term records of stream discharge (Q), specific electrical conductance (EC), and water temperature (T) from 18 streams were analysed in order to quantify the magnitude, duration, and frequency of hydrologic connections over daily, annual, and inter‐annual timescales. At a daily timescale, we observe predictable diurnal variations in Q, EC, and T. At an annual timescale, we observe longer streams to be more intermittent, warmer, and have higher median EC values, compared to shorter streams. Longer streams also behave chemostatically with respect to EC, whereas shorter streams are more strongly characterized by dilution. Inter‐annually, we observe significant variability in annual runoff volumes, likely because of climatic variability over the 20 record years considered. Hydrologic connections at all timescales are vital to stream ecosystem structure and function. This synthesis of hydrologic connectivity in the MDVs provides a useful end‐member template for assessing hydrologic connectivity in more structurally complex temperate watersheds. Copyright © 2016 John Wiley & Sons, Ltd.     

Robust representation of dry cells in single-layer MODFLOW models

Dewatered or "dry" grid cells in the USGS ground water modeling software MODFLOW may cause nonphysical artifacts, trigger convergence failures, or interfere with parameter estimation. These difficulties can be avoided in two dimensions by modifying the spatial differencing scheme and the iterative procedure used to resolve nonlinearities. Specifically, the spatial differencing scheme is modified to use the water level on the upstream side of a pair of adjacent cells to calculate the saturated thickness and hence intercell conductance for the pair. This makes it possible to explicitly constrain the water level in a cell to be at or above the cell bottom elevation without introducing nonphysical artifacts. Thus constrained, all initially active cells will remain active throughout the simulation. It was necessary to replace MODFLOW's Picard iteration method with the Newton-Raphson method to achieve convergence in demanding applications involving many dry cells. Tests using a MODFLOW variant based on the new method produced results nearly identical to conventional MODFLOW in situations where conventional MODFLOW converges. The new method is extremely robust and converged in scenarios where conventional MODFLOW failed to converge, such as when almost all cells dewatered. An example application to the Edwards Aquifer in south-central Texas further demonstrates the utility of the new method.