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.

Performance Analysis of the χMD Matrix Solver Package for MODFLOW-USG.

The χMD matrix solver package is incorporated into USGS groundwater modeling software, such as MODFLOW-NWT, MODFLOW-USG, and MT3D. The solver is used to solve matrices assembled through numerical discretization of the groundwater flow equation, and solute transport equations. χMD has demonstrated its higher robustness, faster execution speed, and more efficient memory usage compared to the existing solvers for many types of groundwater flow problems. χMD uses preconditioned iterative Krylov-subspace methods and consists of preconditioning and acceleration modules. Because the solver package uses a variety of preconditioning features including level-based incomplete lower-upper (ILU) factorization method with a drop tolerance scheme, users must choose optimal preconditioning parameters to improve execution speed and robustness. In order to examine how the preconditioning parameters, ILU factorization level, and drop tolerance values affect the overall performance of the matrix solver, we evaluated five different groundwater model applications using MODFLOW-USG that include different numerical complexities. For those five cases, the number of discretization nodes varied from 10,000 cells to 730,300 cells. From the analysis, we found that the preconditioning parameters greatly affect execution times and memory usage of the preconditioning and acceleration procedures. In addition, a combination of the ILU level between five to seven and the drop tolerance value between 10-2 and 10-3 usually resulted in shorter overall execution time. Our study suggests that the users can elicit higher performance and robustness of the χMD matrix solver using this combination of the parameters and enhance computational efficiency of solving groundwater and solute transport problems.

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.

Simulating the effect of climate extremes on groundwater flow through a lakebed.

Groundwater exchanges with lakes resulting from cyclical wet and dry climate extremes maintain lake levels in the environment in ways that are not well understood, in part because they remain difficult to simulate. To better understand the atypical groundwater interactions with lakes caused by climatic extremes, an original conceptual approach is introduced using MODFLOW-2005 and a kinematic-wave approximation to variably saturated flow that allows lake size and position in the basin to change while accurately representing the daily lake volume and three-dimensional variably saturated groundwater flow responses in the basin. Daily groundwater interactions are simulated for a calibrated lake basin in Florida over a decade that included historic wet and dry departures from the average rainfall. The divergent climate extremes subjected nearly 70% of the maximum lakebed area and 75% of the maximum shoreline perimeter to both groundwater inflow and lake leakage. About half of the lakebed area subject to flow reversals also went dry. A flow-through pattern present for 73% of the decade caused net leakage from the lake 80% of the time. Runoff from the saturated lake margin offset the groundwater deficit only about half of that time. A centripetal flow pattern present for 6% of the decade was important for maintaining the lake stage and generated 30% of all net groundwater inflow. Pumping effects superimposed on dry climate extremes induced the least frequent but most cautionary flow pattern with leakage from over 90% of the actual lakebed area.

Modeling variably saturated multispecies reactive groundwater solute transport with MODFLOW-UZF and RT3D.

A numerical model was developed that is capable of simulating multispecies reactive solute transport in variably saturated porous media. This model consists of a modified version of the reactive transport model RT3D (Reactive Transport in 3 Dimensions) that is linked to the Unsaturated-Zone Flow (UZF1) package and MODFLOW. Referred to as UZF-RT3D, the model is tested against published analytical benchmarks as well as other published contaminant transport models, including HYDRUS-1D, VS2DT, and SUTRA, and the coupled flow and transport modeling system of CATHY and TRAN3D. Comparisons in one-dimensional, two-dimensional, and three-dimensional variably saturated systems are explored. While several test cases are included to verify the correct implementation of variably saturated transport in UZF-RT3D, other cases are included to demonstrate the usefulness of the code in terms of model run-time and handling the reaction kinetics of multiple interacting species in variably saturated subsurface systems. As UZF1 relies on a kinematic-wave approximation for unsaturated flow that neglects the diffusive terms in Richards equation, UZF-RT3D can be used for large-scale aquifer systems for which the UZF1 formulation is reasonable, that is, capillary-pressure gradients can be neglected and soil parameters can be treated as homogeneous. Decreased model run-time and the ability to include site-specific chemical species and chemical reactions make UZF-RT3D an attractive model for efficient simulation of multispecies reactive transport in variably saturated large-scale subsurface systems.

Modeling variably saturated subsurface solute transport with MODFLOW-UZF and MT3DMS.

The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated-zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume-averaged approach in which Lagrangian-based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF-MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably-Saturated Two-Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two- and three-dimensional simulations also were investigated to ensure unsaturated-saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large-scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF-MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three-dimensional variably saturated flow and transport simulations revealed UZF-MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide-spread use of both MODFLOW and MT3DMS, the added capability of unsaturated-zone transport in this familiar modeling framework stands to benefit a broad user-ship.

Approaches to the simulation of unconfined flow and perched groundwater flow in MODFLOW.

Various approaches have been proposed to manage the nonlinearities associated with the unconfined flow equation and to simulate perched groundwater conditions using the MODFLOW family of codes. The approaches comprise a variety of numerical techniques to prevent dry cells from becoming inactive and to achieve a stable solution focused on formulations of the unconfined, partially-saturated, groundwater flow equation. Keeping dry cells active avoids a discontinuous head solution which in turn improves the effectiveness of parameter estimation software that relies on continuous derivatives. Most approaches implement an upstream weighting of intercell conductance and Newton-Raphson linearization to obtain robust convergence. In this study, several published approaches were implemented in a stepwise manner into MODFLOW for comparative analysis. First, a comparative analysis of the methods is presented using synthetic examples that create convergence issues or difficulty in handling perched conditions with the more common dry-cell simulation capabilities of MODFLOW. Next, a field-scale three-dimensional simulation is presented to examine the stability and performance of the discussed approaches in larger, practical, simulation settings.

River-aquifer interactions, geologic heterogeneity, and low-flow management.

Low river flows are commonly controlled by river-aquifer exchange, the magnitude of which is governed by hydraulic properties of both aquifer and aquitard materials beneath the river. Low flows are often important ecologically. Numerical simulations were used to assess how textural heterogeneity of an alluvial system influences river seepage and low flows. The Cosumnes River in California was used as a test case. Declining fall flows in the Cosumnes River have threatened Chinook salmon runs. A ground water-surface water model for the lower river basin was developed, which incorporates detailed geostatistical simulations of aquifer heterogeneity. Six different realizations of heterogeneity and a homogenous model were run for a 3-year period. Net annual seepage from the river was found to be similar among the models. However, spatial distribution of seepage along the channel, water table configuration and the level of local connection, and disconnection between the river and aquifer showed strong variations among the different heterogeneous models. Most importantly, the heterogeneous models suggest that river seepage losses can be reduced by local reconnections, even when the regional water table remains well below the riverbed. The percentage of river channel responsible for 50% of total river seepage ranged from 10% to 26% in the heterogeneous models as opposed to 23% in the homogeneous model. Differences in seepage between the models resulted in up to 13 d difference in the number of days the river was open for salmon migration during the critical fall months in one given year.