Richard G. Niswonger

Analysis of streambed temperatures in ephemeral channels to determine streamflow frequency and duration

Spatial and temporal patterns in streamflow are rarely monitored for ephemeral streams. Flashy, erosive streamflows common in ephemeral channels create a series of operational and maintenance problems, which makes it impractical to deploy a series of gaging stations along ephemeral channels. Streambed temperature is a robust and inexpensive parameter to monitor remotely, leading to the possibility of analyzing temperature patterns to estimate streamflow frequency and duration along ephemeral channels. A simulation model was utilized to examine various atmospheric and hydrological upper boundary conditions compared with a series of hypothetical temperature-monitoring depths within the streambed. Simulation results indicate that streamflow events were distinguished from changing atmospheric conditions with greater certainty using temperatures at shallow depths (e.g., 10–20 cm) as opposed to the streambed surface. Three ephemeral streams in the American Southwest were instrumented to monitor streambed temperature for determining the accuracy of using this approach to ascertain the long-term temporal and spatial extent of streamflow along each stream channel. Streambed temperature data were collected at the surface or at shallow depth along each stream channel, using thermistors encased in waterproof, single-channel data loggers tethered to anchors in the channel. On the basis of comparisons with site information, such as direct field observations and upstream flow records, diurnal temperature variations successfully detected the presence and duration of streamflow for all sites.

MODFLOW-NWT, A Newton formulation for MODFLOW-2005

This report documents a Newton formulation of MODFLOW-2005, called MODFLOW-NWT. MODFLOWNWT is a standalone program that is intended for solving problems involving drying and rewetting nonlinearities of the unconfined groundwater-flow equation. MODFLOWNWT must be used with the Upstream-Weighting (UPW) Package for calculating intercell conductances in a different manner than is done in the Block-Centered Flow (BCF), Layer Property Flow (LPF), or Hydrogeologic-Unit Flow (HUF; Anderman and Hill, 2000) Packages. The UPW Package treats nonlinearities of cell drying and rewetting by use of a continuous function of groundwater head, rather than the discrete approach of drying and rewetting that is used by the BCF, LPF, and HUF Packages. This further enables application of the Newton formulation for unconfined groundwater-flow problems because conductance derivatives required by the Newton method are smooth over the full range of head for a model cell. The NWT linearization approach generates an asymmetric matrix, which is different from the standard MODFLOW formulation that generates a symmetric matrix. Because all linear solvers presently available for use with MODFLOW-2005 solve only symmetric matrices, MODFLOW-NWT includes two previously developed asymmetric matrix-solver options. The matrix-solver options include a generalized-minimum-residual (GMRES) Solver and an Orthomin / stabilized conjugate-gradient (CGSTAB) Solver. The GMRES Solver is documented in a previously published report, such that only a brief description and input instructions are provided in this report. However, the CGSTAB Solver (called χMD) is documented in this report. Flow-property input for the UPW Package is designed based on the LPF Package and material-property input is identical to that for the LPF Package except that the rewetting and vertical-conductance correction options of the LPF Package are not available with the UPW Package. Input files constructed for the LPF Package can be used with slight modification as input for the UPW Package. This report presents the theory and methods used by MODFLOW-NWT, including the UPW Package. Additionally, this report provides comparisons of the new methodology to analytical solutions of groundwater flow and to standard MODFLOW-2005 results by use of an unconfined aquifer MODFLOW example problem. The standard MODFLOW-2005 simulation uses the LPF Package with the wet/dry option active. A new example problem also is presented to demonstrate MODFLOWNWT’s ability to provide a solution for a difficult unconfined groundwater-flow problem.

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.

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.

Mixing effects on nitrogen and oxygen concentrations and the relationship to mean residence time in a hyporheic zone of a riffle‐pool sequence

Flow paths and residence times in the hyporheic zone are known to influence biogeochemical processes such as nitrification and denitrification. The exchange across the sediment-water interface may involve mixing of surface water and groundwater through complex hyporheic flow paths that contribute to highly variable biogeochemically active zones. Despite the recognition of these patterns in the literature, conceptualization and analysis of flow paths and nitrogen transformations beneath riffle-pool sequences often neglect to consider bed form driven exchange along the entire reach. In this study, the spatial and temporal distribution of dissolved oxygen (DO), nitrate ( NO3−) and ammonium ( NH4+) were monitored in the hyporheic zone beneath a riffle-pool sequence on a losing section of the Truckee River, NV. Spatially varying hyporheic exchange and the occurrence of multi-scale hyporheic mixing cells are shown to influence concentrations of DO and NO3− and the mean residence time (MRT) of riffle and pool areas. Distinct patterns observed in piezometers are shown to be influenced by the first large flow event following a steady 8 month period of low flow conditions. Increases in surface water discharge resulted in reversed hydraulic gradients and production of nitrate through nitrification at small vertical spatial scales (0.10–0.25 m) beneath the sediment-water interface. In areas with high downward flow rates and low MRT, denitrification may be limited. The use of a longitudinal two-dimensional flow model helped identify important mechanisms such as multi-scale hyporheic mixing cells and spatially varying MRT, an important driver for nitrogen transformation in the riverbed. Our observations of DO and NO3− concentrations and model simulations highlight the role of multi-scale hyporheic mixing cells on MRT and nitrogen transformations in the hyporheic zone of riffle-pool sequences.

Evaluating mountain meadow groundwater response to Pinyon-Juniper and temperature in a Great Basin watershed†

This research highlights development and application of an integrated hydrologic model (GSFLOW) to a semiarid, snow-dominated watershed in the Great Basin to evaluate Pinyon-Juniper (PJ) and temperature controls on mountain meadow shallow groundwater. The work used Google Earth Engine Landsat satellite and gridded climate archives for model evaluation. Model simulations across three decades indicated that the watershed operates on a threshold response to precipitation (P) > 400 mm y-1 to produce a positive yield (P-ET; 9%) resulting in stream discharge and a rebound in meadow groundwater levels during these wetter years. Observed and simulated meadow groundwater response to large P correlates with above average predicted soil moisture and with a normalized difference vegetation index (NDVI) threshold value > 0.3. A return to assumed pre-expansion PJ conditions or an increase in temperature to mid-21st century shifts yielded by only ±1% during the multi-decade simulation period; but changes of approximately ±4% occurred during wet years. Changes in annual yield were largely dampened by the spatial and temporal redistribution of evapotranspiration (ET) across the watershed. Yet, the influence of this redistribution and vegetation structural controls on snowmelt altered recharge to control water table depth in the meadow. Even a small-scale removal of PJ (0.5 km2) proximal to the meadow will promote a stable, shallow groundwater system resilient to droughts, while modest increases in temperature will produce a meadow susceptible to declining water levels and a community structure likely to move toward dry and degraded conditions. This article is protected by copyright. All rights reserved.

Groundwater exchanges near a channelized versus unmodified stream mouth discharging to a subalpine lake

The terminus of a stream flowing into a larger river, pond, lake, or reservoir is referred to as the stream-mouth reach or simply the stream mouth. The terminus is often characterized by rapidly changing thermal and hydraulic conditions that result in abrupt shifts in surface water/groundwater (sw/gw) exchange patterns, creating the potential for unique biogeochemical processes and ecosystems. Worldwide shoreline development is changing stream-lake interfaces through channelization of stream mouths, i.e., channel straightening and bank stabilization to prevent natural meandering at the shoreline. In the central Sierra Nevada (USA), Lake Tahoes shoreline has an abundance of both “unmodified” (i.e., not engineered though potentially impacted by broader watershed engineering) and channelized stream mouths. Two representative stream mouths along the lakes north shore, one channelized and one unmodified, were selected to compare and contrast water and heat exchanges. Hydraulic and thermal properties were monitored during separate campaigns in September 2012 and 2013 and sw/gw exchanges were estimated within the stream mouth-shoreline continuum. Heat-flow and water-flow patterns indicated clear differences in the channelized versus the unmodified stream mouth. For the channelized stream mouth, relatively modulated, cool-temperature, low-velocity longitudinal streambed flows discharged offshore beneath warmer buoyant lakeshore water. In contrast, a seasonal barrier bar formed across the unmodified stream mouth, creating higher-velocity subsurface flow paths and higher diurnal temperature variations relative to shoreline water. As a consequence, channelization altered sw/gw exchanges potentially altering biogeochemical processing and ecological systems in and near the stream mouth.

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.

Feedback of land subsidence on the movement and conjunctive use of water resources

The dependency of surface- or groundwater flows and aquifer hydraulic properties on dewatering-induced layer deformation is not available in the USGSs groundwater model MODFLOW. A new integrated hydrologic model, MODFLOW-OWHM, formulates this dependency by coupling mesh deformation with aquifer transmissivity and storage and by linking land subsidence/uplift with deformation-dependent flows that also depend on aquifer head and other flow terms. In a test example, flows most affected were stream seepage and evapotranspiration from groundwater (ETgw). Deformation feedback also had an indirect effect on conjunctive surface- and groundwater use components: Changed stream seepage and streamflows influenced surface-water deliveries and returnflows. Changed ETgw affected irrigation demand, which jointly with altered surface-water supplies resulted in changed supplemental groundwater requirements and pumping and changed return runoff. This modeling feature will improve the impact assessment of dewatering-induced land subsidence/uplift (following irrigation pumping or coal-seam gas extraction) on surface receptors, inter-basin transfers, and surface-infrastructure integrity. We develop a method to simulate deformation-dependent flows for MODFLOW.We demonstrate the significance of linking subsidence to conjunctive water use.The linkages affect flows across the landscape, surface water, and groundwater.Linked flows are relevant to resource issues that include conjunctive water use.

Managed aquifer recharge through off‐season irrigation in agricultural regions:

Options for increasing reservoir storage in developed regions are limited and prohibitively expensive. Projected increases in demand call for new long-term water storage to help sustain agriculture, municipalities, industry, and ecological services. Managed aquifer recharge (MAR) is becoming an integral component of water resources around the world. However, MAR faces challenges, including infrastructure costs, difficulty in enhancing recharge, water quality issues, and lack of available water supplies. Here we examine, through simulation modeling of a hypothetical agricultural subbasin in the western US, the potential of agricultural managed aquifer recharge (Ag-MAR) via canal seepage and off-season field irrigation. Weather phenomenon in many regions around the world exhibit decadal and other multi-year cycles of extreme precipitation. An on-going challenge is to develop approaches to store greater amounts of water during these events. Simulations presented herein incorporate Ag-MAR programs and demonstrate that there is potential to enhance regional recharge by 7-13%, increase crop consumptive use by 9-12%, and increase natural vegetation consumption by 20-30%, where larger relative increases occur for lower aquifer hydraulic conductivity and higher specific yield values. Annual Increases in groundwater levels were 7 m, and sustained levels following several years of drought were greater than 2 m. Results demonstrate that Ag-MAR has great potential to enhance long-term sustainability of water resources in agricultural basins.