Geology Reports⌕ Search

Geology topics

Eric D. Morway

Publications and source records attributed to Eric D. Morway.

At least 19 recordsLinked to original sources

Assessing nonpoint-source uranium pollution in an irrigated stream-aquifer system

Uranium (U) in rocks and soils of arid and semi-arid environments can be mobilized by irrigation and fertilization, posing environmental and health risks. Elevated U, along with selenium (Se) and nitrate (NO 3 ) co-constituents, necessitates careful monitoring and management. We developed a distributed-parameter numerical model to assess U pollution in an irrigated stream-aquifer system, applying it to a 552 km 2 region in Colorado's Lower Arkansas River Valley (LARV) over 14 years. A MODFLOW model, describing groundwater and stream flow, was coupled with an RT3D-OTIS model to portray reactive U transport. Calibration using the PESTPP-iES iterative ensemble smoother (iES) software indicated good agreement with observed U concentrations. The model revealed substantial and variable U levels across the LARV, highlighting potential hotspots and possible contributing factors, such as geological composition of the bedrock and near-surface shale and aquifer sediments derived from them, irrigation practices, and riparian landscape. U levels exceed the chronic standard (85th percentile = 30 μg/L, set by the US Environmental Protection Agency), which is the permissible regulatory threshold, in groundwater across 44 % of the region and along the river by an average factor of 2.9. Simulated average U concentrations in the non-riparian aquifer and river are 124 μg/L and 60 μg/L, respectively, compared with 112 μg/L and 62 μg/L for measured values. The average 85th percentile U concentration is 222 μg/L in the aquifer and 82 μg/L in the river. Average simulated U mass loading to the river is 0.17 kg/day per km, compared to an estimated 0.23 kg/day per km. Findings provide a baseline for comparing future simulated outcomes of alternative best management practices (BMPs) for U pollution mitigation and offer a methodology applicable to other irrigated regions.

Colorado↗

A new groundwater energy transport model for the MODFLOW hydrologic simulator

Heat transport in the subsurface is an important aspect of research related to the effects of a warming climate on ecological services (i.e., cold-water refugia); the development of geothermal resources for energy banking schemes (i.e., aquifer thermal energy storage [ATES]); and the effects of temperature on other aspects of groundwater quality, such as nutrient cycling. Historically, simulation of heat transport using the MODFLOW groundwater simulator and related codes was performed by scaling the input parameters of a solute-transport model to emulate heat transport. However, that approach required additional pre- and post-processing of input and output and could not account for the variation in effective thermal storage and transport properties during transient, unsaturated flow, for example. True heat-transport capabilities in the context of MODFLOW were first introduced in a variant called USG-Transport. More recently, a new groundwater energy-transport (GWE) model type has been added to MODFLOW 6, the core version of the MODFLOW hydrologic simulator. GWE supports the simulation of heat transport on structured or unstructured grids as well as within and between features of advanced packages that represent streams, lakes, multi-aquifer wells, and the unsaturated zone. GWE is integrated within MODFLOW 6 and is accessible through the FloPy Python package and the MODFLOW 6 application programming interface (API). An example simulation demonstrates conduction between grid cells through both the water and the solid aquifer material, including thermal bleeding from saturated overburden cells into a groundwater flow field.

Groundwater↗

Remote sensing evapotranspiration in ensemble-based framework to enhance cascade routing and re-infiltration concept in integrated hydrological model applied to support decision making

Integrated hydrological models (IHMs) help characterize the complexity of surface–groundwater interactions. The cascade routing and re-infiltration (CRR) concept, recently applied to a MODFLOW 6 IHM, improved conceptualization and simulation of overland flow processes. The CRR controls the transfer of rejected infiltration and groundwater exfiltration from upslope areas to adjacent downslope areas where that water can be evaporated, re-infiltrated back to subsurface, or discharged to streams as direct runoff. The partitioning between these three components is controlled by uncertain parameters that must be estimated. Thus, by quantifying and reducing those uncertainties, next to uncertainties of the other model parameters (e.g. hydraulic and storage parameters), the reliability of the CRR is improved and the IHM is better suited for decision support modelling, the two key objectives of this work. To this end, the remotely sensed MODIS-ET product was incorporated into the calibration process for complementing traditional hydraulic head and streamflow observations. A total of approximately 150,000 observations guided the calibration of a 13-year MODFLOW 6 IHM simulation of the Sardon catchment (Spain) with daily stress periods. The model input uncertainty was represented by grid-cell-scale parameterization, yielding approximately 500,000 unknown input parameters to be conditioned. The calibration was carried out through an iterative ensemble smoother. Incorporating the MODIS-ET data improved the CRR implementation, and reduced uncertainties associated with other model parameters. Additionally, it significantly reduced the uncertainty associated with net recharge, a critical flux for water management that cannot be directly measured and rather is commonly estimated by IHM simulations.

Sardon catchment↗

Assessing potential effects of changes in water use in the middle Carson River Basin with a numerical groundwater-flow model, Eagle, Dayton, and Churchill Valleys, west-central Nevada

During the economic boom of the mid part of the first decade of the 2000s in northwestern Nevada, municipal and housing growth increased use of the water resources of this semi-arid region. In 2008, when the economy slowed, new housing development stopped, and immediate pressure on groundwater resources abated. The U.S. Geological Survey, in cooperation with the Bureau of Reclamation, began a hydrogeologic study of the middle Carson River Basin. The first half of the study reviewed and synthesized previous geologic studies and contributed new datasets that served as a foundation for a three-dimensional, transient numerical model of groundwater and surface-water flow for the middle Carson River Basin extending from Eagle Valley to Churchill Valley. The model can be used to evaluate the effects of proposed alternative management strategies on groundwater sustainability, flows in the Carson River, and routine operation of Lahontan Reservoir and can also provide a basis for basin-wide investigations seeking to quantitatively evaluate the effects of climate change or yet-to-be-determined alternative management strategies. The middle Carson model was constructed using the U.S. Geological Survey groundwater modeling software MODFLOW-NWT. MODFLOW is widely used groundwater modeling software and is well-suited for evaluating groundwater and surface-water interactions. The model uses 550-feet square grid cells that align with the previously published model for Carson Valley (adjacent upstream valley). Six grid layers with more finely resolved vertical resolution near the perimeter of the active model domain and near surface-water features, compared to other areas of the active model domain, hone the simulated groundwater and surface-water exchanges. In addition to simulating groundwater and surface-water interaction, crop and phreatophyte evapotranspiration, lake evaporation, mountain-front recharge, recharge from irrigation return flows, and groundwater pumping are also simulated. Surface-water flow entering the model domain, including the Carson River, tributary inflow from perennial streams in Eagle Valley, and trans-basin imports through the Truckee Canal (surface water diverted from the Truckee River) are specified according to U.S. Geological Survey streamgage records. Groundwater pumpage and surface-water diversions to 10 agricultural ditches and the managed release from Lahontan Reservoir, at the end of the middle Carson River Basin, are specified according to water-manager records. The model simulation period extended from 2000 through 2010 (January 1, 2000, to December 31, 2010) using 574 weekly stress periods, with a single steady-state stress period at the beginning of the simulation that establishes initial conditions by approximating average conditions during the transient simulation period. All available observations for this period were used during the model calibration process, performed using automated parameter-estimation software. Calibration targets included observations of groundwater elevations in wells, streamflow, differences in observed streamflow between successive streamgages and actual evapotranspiration from irrigated lands. Among all 5,296 simulated and observed groundwater level pairs, the mean error was 1.42 feet; the mean absolute error, 7.71 feet; and the percent bias was −0.1 percent. Three alternative management scenarios, run using the entire period of analysis (2000–10), were simulated to improve understanding of the potential effects of (1) loss of irrigated agricultural lands following conversion of water-rights to municipal groundwater rights; (2) reclaiming treated wastewater with induction wells; and (3) exercising permitted but under-utilized groundwater rights. Scenarios 2 and 3 were further explored using two and four subscenarios, respectively. Simulated scenario results ranged from having little effect on the groundwater system relative to a baseline simulation to having spatially extensive and large groundwater-level declines (10 to 20 feet) compared to the baseline simulation. None of the simulated scenarios increased delivery of river flows to Lahontan Reservoir. On the contrary, one of the subscenarios under alternative management scenario 3 led to surface-water delivery shortfalls of more than 10,000 acre-feet per year. Future model improvements may include an extension of the model simulation period backward and forward in time and directly linking it to the upstream Carson Valley groundwater model. Furthermore, converting this MODFLOW model to a GSFLOW model, which fully integrates groundwater and surface-water flows including precipitation runoff and infiltration, may provide an improved tool for comprehensive management of water-resources in the middle Carson River Basin.

Nevada↗

Long-term mercury loading and trapping dynamics in a Western North America reservoir

Study Region The Carson River including Lahontan Reservoir in Northwestern Nevada, USA Study Focus The discovery, mining, and processing of silver and gold from the Comstock Lode in northwestern Nevada heavily contaminated the Carson River system with mercury (Hg). The river now contains some of the highest recorded water column and bed sediment Hg concentrations reported in peer-reviewed literature. Acute Hg contamination in river and reservoir bed sediment has led to elevated methylmercury (MeHg) concentrations across all trophic levels of Lahontan Reservoir, culminating in significant health risks to humans. Lahontan Reservoir, located downstream from the mills that processed Comstock ore, has served as a Hg trap since the completion of the dam in 1915. Starting in 1997 and extending through 2021, the U.S. Geological Survey has collected and analyzed hundreds of discrete water samples entering and exiting Lahontan Reservoir for unfiltered total Hg (uf.THg), thereby providing a first-of-its-kind opportunity for studying long-term Hg trapping efficiencies within a western North American reservoir. Continuous time series of uf.THg concentration and flux above and below the reservoir were estimated using the weighted regressions on time, discharge, and season with the Kalman filtering (WRTDS-K) method employed with newly developed methods to account for non-natural (i.e., managed) hydrographs below a reservoir. New hydrological insights for the region An estimated 31,650 kg (34.89 US tons) of uf.THg entered Lahontan Reservoir during the 25 year period of analysis, accounting for approximately 0.5% of the total uf.THg estimated to have been released to the Carson River system [6.8 million kg (7500 US tons)] over a multi-decade mining boom. Moreover, approximately 92% of the estimated uf.THg entering Lahontan Reservoir was trapped. On an annual basis, however, trapping efficiencies range between 34% and 98%, and are closely related to the total annual discharge. Results also indicate that flow-normalized uf.THg concentrations and loads above and below the reservoir are trending down.

Nevada↗

New capabilities in MT3D-USGS for simulating unsaturated-zone heat transport

Changes in climate and land use will alter groundwater heat transport dynamics in the future. These changes will in turn affect watershed processes (e.g., nutrient cycling) as well as watershed characteristics (e.g., distribution and persistence of cold-water habitat). Thus, groundwater flow and heat transport models at watershed scales that can characterize and quantify thermal impacts of surface temperature change on groundwater system temperatures are needed to forecast changes to groundwater-linked ecosystems in riparian zones, streams, and lakes. Including unsaturated zone processes has previously been shown to be important for properly determining the timing and magnitude of groundwater recharge (Hunt et al. 2008). Similarly, heat transport dynamics in the saturated-zone, as well as connected surface-water systems, can be appreciably influenced by unsaturated-zone processes; in this way the unsaturated zone forms an inextricable link between land surface where change occurs and the groundwater system that transmit that change. This paper presents new capabilities for the existing MT3D-USGS transport simulator by adding functionality for simulating heat transport through the unsaturated zone. New simulation capabilities are verified through comparison of simulation results with those of the variably-saturated heat transport simulator VS2DH under steady and transient conditions for both water and heat flow. The new capabilities are assessed using a number of conceptualizations and include evaluations of convective and conductive heat flow. These additional capabilities increase the utility for applied watershed-scale simulations, which in turn should facilitate more realistic characterizations of temperature change on thermally sensitive ecosystems, such as stream habitat.

Groundwater↗

Simulation of heat flow in a synthetic watershed: Lags and dampening across multiple pathways under a climate-forcing scenario

Although there is widespread agreement that future climates tend toward warming, the response of aquatic ecosystems to that warming is not well understood. This work, a continuation of companion research, explores the role of distinct watershed pathways in lagging and dampening climate-change signals. It subjects a synthetic flow and transport model to a 30-year warming signal based on climate projections, quantifying the heat breakthrough on a monthly time step along connected pathways. The system corresponds to a temperate watershed roughly 27 km on a side and consists of (a) land-surface processes of overland flow, (b) infiltration through an unsaturated zone (UZ) above an unconfined sandy aquifer overlying impermeable bedrock, and (c) groundwater flow along shallow and deep pathlines that converge as discharge to a surface-water network. Numerical simulations show that about 40% of the warming applied to watershed infiltration arrives at the water table and that the UZ stores a large fraction of the upward-trending heat signal. Additionally, once groundwater reaches the surface-water network after traveling through the saturated zone, only about 10% of the original warm-up signal is returned to streams by discharge. However, increases in the simulated streamflow temperatures are of similar magnitude to increases at the water table, due to the addition of heat by storm runoff, which bypasses UZ and groundwater storage and counteracts subsurface dampening. The synthetic modeling method and tentative findings reported here provide a potential workflow for real-world applications of climate-change modeling at the full watershed scale.

Water↗

The MODFLOW Application Programming Interface for simulationcontrol and software interoperability

The MODFLOW API allows other programs to control MODFLOW and interactively change variables without having to modify the source code. The MODFLOW API is based on the Basic Model Interface (BMI), which is a set of conventions that define how to initialize a simulation, update the model state by advancing in time, and finalize the run. For many existing MODFLOW coupling applications, the information provided to MODFLOW must be updated multiple times in a time step. As this capability to modify variables within a time step is not defined by the BMI, an extension to BMI was developed. This eXtended Model Interface is part of the MODFLOW API and allows such a tight coupling to other models. Examples are included for a variety of use cases, including new flexibility for users to develop custom packages without modifying the MODFLOW source code and coupling MODFLOW with other models and optimization libraries.

Environmental Modelling & Software↗

Integrated hydrology and operations modeling to evaluate climate change impacts in an agricultural valley irrigated with snowmelt runoff

Applying models to developed agricultural regions remains a difficult problem because there are no existing modeling codes that represent both the complex physics of the hydrology and anthropogenic manipulations to water distribution and consumption. We apply an integrated groundwater – surface water and hydrologic river operations model to an irrigated river valley in northwestern Nevada/northern California, United States to evaluate the impacts of climate change on snow-fed agricultural systems that use surface water and groundwater conjunctively. We explicitly represent individual surface water rights within the hydrologic model and allow the integrated code to change river diversions in response to earlier snowmelt runoff and water availability. Historically under-used supplemental groundwater rights are dynamically activated within the model to offset diminished surface water deliveries. The model accounts for feedbacks between the natural hydrology and anthropogenic stresses, which is a first-of-its-kind assessment of the impacts of climate change on individual water rights, and more broadly on river basin operations. Earlier snowmelt decreases annual surface water deliveries to all water rights, not just the junior water rights, owing to a lack of surface water storage in the upper river basin capable of capturing earlier runoff. Conversely, downstream irrigators with access to reservoir storage benefit from earlier runoff flowing past upstream points of diversion prior to the start of the irrigation season. Despite regional shifts toward greater reliance on groundwater for irrigation, crop consumption (a common surrogate for crop yield) decreases due to spatiotemporal changes in water supply that preferentially impact a subset of growers in the region.

California↗

Use of the MODFLOW 6 water mover package to represent natural and managed hydrologic connections

The latest release of MODFLOW 6, the current core version of the MODFLOW groundwater modeling software, debuted a new package dubbed the “mover” (MVR). Using a generalized approach, MVR facilitates the transfer of water among any arbitrary combination of simulated features (i.e., pumping wells, stream, drains, lakes, etc.) within a MODFLOW 6 simulation. Four “rules” controlling the amount of water transferred from a providing feature to a receiving feature are currently available. In this way, MVR can represent natural connections between features, for example streams entering or exiting lakes, and perhaps more interestingly, it also can transfer water among simulated features to more accurately simulate water management. An example model representative of an agricultural setting demonstrates some of the available MVR connections. For example, an irrigation event that transfers surface water from an irrigation delivery ditch to multiple cropped areas demonstrates a “one-to-many” connection that is possible within MVR. Conversely, irrigation or precipitation runoff from multiple fields may be routed to a particular stream segment using “many-to-one” MVR connections. MVR supports many additional connection types, several of which are demonstrated by the included example problem.

Groundwater↗

Trends in nitrogen, phosphorus, and sediment concentrations and loads in streams draining to Lake Tahoe, California, Nevada, USA

Lake Tahoe, a large freshwater lake of the eastern Sierra Nevada in California and Nevada, has 63 tributaries that are sources of nutrients and sediment to the lake. The Tahoe watershed is relatively small, and the surface area of the lake occupies about 38% of the watershed area (1313 km 2 ). Only about 6% of the watershed is urbanized or residential land, and as part of a plan to maintain water clarity, wastewater is exported out of the basin. The lake's clarity has been diminishing due to algae and fine sediment, prompting development of management plans. Much of the annual discharge and nutrient load to the lake results from snowmelt in the spring and summer months. To understand the relative importance of land use, climate, forest management, and other factors affecting trends in nutrient stream concentrations and loads, a Weighted Regression on Time Discharge and Season (WRTDS) model simulated these trends over a time frame of >25 years (mid-1970s to 2017). All studied locations generally show nitrate concentration and load trending down. Ammonium concentration and load initially trended down then increased continuously after 2005. Some locations show initially decreasing orthophosphate trends, followed by small significant increases in concentration and loads starting around 2000 to 2005. Total Kjeldahl nitrogen, total phosphorus and suspended sediment mostly trended downward. Overall, the trends in various forms of nitrogen were observed at most sites irrespective of the degree of development and indicate a change in ecological conditions is affecting the nitrogen cycle throughout the watershed, most likely attributable to forest aggradation and fire suppression. Ratios of bioavailable nitrogen in the form of nitrate and ammonium to orthophosphate have also trended downward during the period of record suggesting a shift of these streams from phosphorus limited to nitrogen limited.

California, Nevada↗

Managed aquifer recharge in snow-fed river basins: What, why and how?

What does climate change mean for snow-fed river basins? Climate change poses unique challenges in snow-fed river basins across the western United States because the majority of water supply originates as snow (Dettinger, Udall, & Georgakakos, 2015). In the Sierra Nevada, recent observations include changes in snow accumulation and snowmelt, and shifts in peak streamflow timing (Barnhart et al., 2016; Hatchett et al., 2017; Kim & Jain, 2010; McCabe, Wolock, & Valentin, 2018; Mote, Li, Lettenmaier, Xiao, & Engel, 2018). Such changes upstream alter surface water deliveries downstream, as well as groundwater recharge utilized as both primary and supplemental water supply (Godsey et al., 2014; Harpold, 2016; Jasechko et al., 2014). basin where snowmelt runoff produces substantial water supply to meet diverse agricultural, environmental and urban water demand (Figure 1). The East and West Forks join at the confluence of the Carson River near the north end of the Carson Valley, a rich agricultural region (40,000 acres) that grows primarily alfalfa hay. The majority of irrigators rely on surface water delivered through a network of earthen ditches constructed in the mid-19th and early 20th centuries. Flow through these earthen networks and the practice of flood irrigation contribute significantly to groundwater recharge. Because no upstream surface water reservoirs exist, snowpack that accumulates through winter and melts slowly through spring has acted as a “natural” reservoir, providing ample supply through the summer irrigati agricultural, environmental and urban water demand (Figure 1). The East and West Forks join at the confluence of the Carson River near the north end of the season. Some irrigators have permitted access to supplemental groundwater that is useful during periods of drought for augmenting shortfalls in surface water delivery. Groundwater is the primary source of municipal and industrial water supply for surrounding communities (e.g., Carson City, Minden, Gardnerville, Dayton). Across the basin, water use is highly regulated through federal, tribal, state and local water-sharing agreements based on prior appropriation doctrine (Wilds, 2014). Carson River surface water allocations follow the Alpine Decree, initiated by the United States Department of Interior in 1925 and signed into law in 1980, following 55 years of litigation, to adjudicate surface water rights to individual parties (NDWP, 1999). The Alpine Decree acknowledges return flows to lower river segments, and thus each river segment is distributed autonomously. This means that the most junior water right on an upper segment can be fulfilled before considering the most senior water right on a lower segment. Ultimately, the ruling is at the discretion of the Federal Water Master to satisfy the needs of each water right Downstream of Carson Valley, surface water flows are stored in Lahontan Reservoir, the nation’s first desert reclamation project (est. 1906), where releases are managed to meet the Newland’s Project irrigation water demand and for environmental use on the Stillwater National Wildlife Refuge. Flows from the Carson River are supplemented through diversions from the Truckee River via the Truckee Canal, resulting in a trans-basin water supply system. How is the Water for the Seasons research program informing snow-fed river basin communities? In the Truckee-Carson River System, researchers and local water managers are working together to assess climate change impacts to water supply and explore how model simulations can produce useful information to support local climate adaptation. Twelve key water managers represent agricultural, environmental, urban and regulatory water-use communities, and bring to the table diverse input and perspectives on how to adapt to climate change. Hydrologists use this input to craft scenarios and simulations that meet the information needs of local water managers. Biannual workshops provide an opportunity for information exchange, where researchers and key water managers generate new knowledge of river system function. That is, researchers share results of models that examine the physical potential, and managers validate the on-the-ground potential, further informing the research process. Coincident to this research program, the region faced a prolonged drought period (2012-2016) with historically low snowpack, followed by a historic wet year (2017) that brought winter and spring flooding as a result of atmospheric river storm events (Sterle et al., 2019). For the Carson River, an important observation made by managers was that peak streamflow that had traditionally coincided with peak irrigation demands, had shifted to earlier in the spring, with summer baseflow also decreasing (Sterle & Singletary, 2017). Managers shared with researchers concerns over potential future impacts that changing snowpack will have on surface water deliveries and reliance on groundwater, as the region’s population and economy continue to grow. During workshops that occurred over this period, local water managers and researchers discussed ways to evaluate water distribution and use that honors the existing legal framework and accounts for changing snowpack regimes (amount, rain versus snow, timing). In response to managers growing interest, researchers introduced the concept of managed aquifer recharge as one potential strategy to adapt and enhance regional water sustainability. What is managed aquifer recharge? Simply stated, managed aquifer recharge is the intentional recharge of structures to spread water over agricultural lands, allowing water to naturally infiltrate into the groundwater system (Bouwer, 1999; Niswonger et al., 2017). The latter may occur during the irrigation season by applying excess water, or during the nonirrigation season when evapotranspiration losses are low. Figure 2 illustrates managed aquifer recharge in a snow-fed river basin, where streamflow generated from snowmelt runoff is diverted to agricultural lands to recharge the aquifer. Such flood irrigation practices, including water delivery through earthen ditch networks, provide incidental but significant aquifer recharge through seepage and deep drainage beneath fields (Niswonger, Allander, & Jeton, 2014). The effects of managed aquifer recharge can vary depending on the location and intensity of practice. For example, implementing managed aquifer recharge water into the groundwater system (Dillon, 2009). This differs from the incidental recharge that may occur as part of normal irrigation practices. Managed recharge may occur by injection into the aquifer through existingwells, or by using existing conveyance adjacent to/along the river’s floodplain has the potential to enhance late-season instream flows due to increased return flows, resulting in greater downstream deliveries as well as improving ecological conditions (Niswonger et al., 2017). Implementing managed aquifer recharge away from the river’s floodplain has the potential to enhance groundwater supply which is increasingly relied upon during surface water shortage (Green et al., 2011), by storing water in available aquifer space in the deep aquifer. At the basin scale, managed aquifer recharge may lead to regional groundwater sustainability. Is the Carson River Basin a candidate for managed aquifer recharge? The physical limitations to implementing managed aquifer recharge in the Carson River Basin hinges on three key factors. The first factor relates to the physical connectivity between rivers and streams, and the irrigation delivery network of canals and ditches that divert water to agricultural lands (Niswonger et al., 2017). In the Carson River Basin the mechanisms for getting water to fields is already in place. Thus, intentionally routing high flows that occur in wet years through this system during the nonirrigation season would mimic what occurs naturally during the irrigation season. The second factor relates to the occurrence of atmospheric river storm events that deliver large amounts of precipitation to the region, much greater than average (Dettinger et al., 2015). With increased frequency and intensity projected under a warmer climate, such events have the potential to produce excess water over short periods of time that could be stored through mechanisms such as managed aquifer recharge (Niswonger et al., 2017). The third factor relates to the change in snowpack accumulation and shifts in snowmelt timing observed elsewhere in the Sierra Nevada (e.g., Godsey et al., 2014; Mote et al., 2018). Having a mechanism in place to maximize use of earlier snowmelt and shifts in streamflow timing could be advantageous and enhance regional groundwater sustainability. As part of the Water for the Seasons study, a hypothetical scenario was developed to determine the feasibility of managed aquifer recharge in the Carson River Basin, assuming no legal constraints. During “wet” or above-average water years, irrigators in the Upper Carson Valley would divert high flows and spread water over agricultural lands during the nonirrigation season. Assuming flows are abundant and “early,” diversions would begin prior to the growing season, when water would otherwise flow downstream to the Lahontan Reservoir. During “dry” years or drought periods, when surface water availability is less, irrigators in the Upper Carson Valley could augment surface water shortages with groundwater, allowing available surface water flows to flow downstream. Researchers hypothesize the amount of water has the potential to boost baseflow to support environmental instream flows, for example. What concerns have local water managers expressed? The hypothetical managed aquifer recharge scenario was presented to water managers in a workshop setting. Presentations included an overview of the hydrologic and operations modeling tools used to evaluate managed aquifer recharge by simulating the timing and distribution of water in the upper watershed. Specifically, in the Upper Carson Valley, a hydrologic model (GSFLOW) simulates streamflow driven by snowmelt, and surface and groundwater interactions, while a river basin operations model (MODSIM) allocates water according to the prior appropriation doctrine in the basin (see Figure 1) (Morway, Niswonger, & Triana, 2016; Niswonger et al., 2017). Integrating these two modeling tools advances the evaluation of climate impacts on water availability in agricultural communities and the resulting impacts of alternative management strategies (Morway et al., 2016). When asked about the viability of managed aquifer recharge, the perspectives of 11 managers varied (Figure 3). Regardless of rating, all managers questioned, “How would thisreally work?” Several managers questioned whether models could simulate the connectivity between surface and groundwater to accurately quantify changes to instream flow. Others raised concerns that managed aquifer recharge violates the Alpine Decree and Nevada Water Law. Still others requested researchers consider alternatives that could work within the confines of current (2019) water law. Managers posed specific questions that should be considered when evaluating the potential for managed aquifer recharge. For example: What triggers implementation of managed aquifer recharge?How “high” or “low” must annual flows be to initiate managed aquifer recharge? When in the water year is this determined? Where exactly in the Carson Valley is managed aquiferre charge possible? For example, what areas away from the floodplain could ensure long-term storage? Can model simulations quantify potential benefits and consequences system-wide?Would this information support decision-making, such as permitting of additional supplemental groundwater rights? How are researchers going to address managers’ research questions? Managers’ perspectives help to validate the on-the-ground potential of particular strategies and further refine alternative management scenarios. For example, understanding that managers are concerned with oversaturated fields helps researchers to define conditions in the model, such as what defines a wet versus “too” wet type of year and where to focus irrigation for managed aquifer recharge. Incorporating these nuances provides more accurate quantification of the potential benefits and consequences for users across the basin. Modeling is underway to simulate managed aquifer recharge scenarios and explore basin-wide implications. Researchers and local water managers will convene to collaboratively review results and further assess whether this or other strategies could work under the confines of existing water law. Subsequent fact sheets will present these findings.

Nevada↗

Simulating selenium and nitrogen fate and transport in coupled stream-aquifer systems of irrigated regions

Elevated levels of selenium (Se) in aqueous environments can harm aquatic life and endanger livestock and human health. Although Se occurs naturally in the rocks and soils of many alluvial aquifers, mining and agricultural activities can increase its rate of mobilization and transport to surface waters. Attention is given here to regions where nonpoint source return flows from irrigated lands carry pollutant loads to aquifers and streams, contributing to concentrations that violate regulatory and performance standards. Of particular concern is the heightened level and mobilization of Se influenced by nitrate (NO 3 ), a harmful pollutant in its own right. We present a numerical model that simulates the reactive transport of Se and nitrogen (N) species in a coupled groundwater-surface water system. Building upon a conceptual model that incorporates the major processes affecting Se and NO 3 transport in an irrigated watershed, the model links the finite-difference models MODFLOW, UZF-RT3D, and OTIS, to simulate flow and reactive transport of multiple chemical species in both the aquifer and a stream network, with mass exchange between the two. The capability of the new model is showcased by calibration, testing, and application to a 500 km 2 region in Colorado’s Lower Arkansas River Valley using a rich data set gathered over a 10-yr period. Simulation of spatial and temporal distributions of Se concentration reveals conditions that exceed standards in groundwater for approximately 20% of the area. For the Arkansas River, standards are exceeded by 290%–450%. Simulation indicates that river concentrations of NO 3 alone are near the current interim standard for the total of all dissolved N species. These results indicate the need for future use of the developed model to investigate the prospects for land and water best management practices to decrease pollutant levels.

Colorado↗

Evaluation of bias associated with capture maps derived from nonlinear groundwater flow models

The impact of groundwater withdrawal on surface water is a concern of water users and water managers, particularly in the arid western United States. Capture maps are useful tools to spatially assess the impact of groundwater pumping on water sources (e.g., streamflow depletion) and are being used more frequently for conjunctive management of surface water and groundwater. Capture maps have been derived using linear groundwater flow models and rely on the principle of superposition to demonstrate the effects of pumping in various locations on resources of interest. However, nonlinear models are often necessary to simulate head-dependent boundary conditions and unconfined aquifers. Capture maps developed using nonlinear models with the principle of superposition may over- or underestimate capture magnitude and spatial extent. This paper presents new methods for generating capture difference maps, which assess spatial effects of model nonlinearity on capture fraction sensitivity to pumping rate, and for calculating the bias associated with capture maps. The sensitivity of capture map bias to selected parameters related to model design and conceptualization for the arid western United States is explored. This study finds that the simulation of stream continuity, pumping rates, stream incision, well proximity to capture sources, aquifer hydraulic conductivity, and groundwater evapotranspiration extinction depth substantially affect capture map bias. Capture difference maps demonstrate that regions with large capture fraction differences are indicative of greater potential capture map bias. Understanding both spatial and temporal bias in capture maps derived from nonlinear groundwater flow models improves their utility and defensibility as conjunctive-use management tools.

Groundwater↗

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 U.S., 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 multiyear cycles of extreme precipitation. An ongoing 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.

Water Resources Research↗

Long-term trends of surface-water mercury and methylmercury concentrations downstream of historic mining within the Carson River watershed

The Carson River is a vital water resource for local municipalities and migratory birds travelling the Pacific Flyway. Historic mining practices that used mercury (Hg) to extract gold from Comstock Lode ore has left much of the river system heavily contaminated with Hg, a practice that continues in many parts of the world today. Between 1998 and 2013, the United States Geological Survey (USGS) collected and analyzed Carson River water for Hg and methylmercury (MeHg) concentrations resulting in a sixteen year record of unfiltered total mercury (uf.THg), filtered (dissolved) Hg (f.THg), total methylmercury (uf.MeHg), filtered MeHg (f.MeHg), and particulate-bound THg (p.THg) and MeHg (p.MeHg) concentrations. This represents one of the longest continuous records of Hg speciation data for any riverine system, thereby providing a unique opportunity to evaluate long-term trends in concentrations and annual loads. During the period of analysis, uf.THg concentration and load trended downward at rates of −0.85% and −1.8% per year, respectively. Conversely, the f.THg concentration increased at a rate of 1.7% per year between 1998 and 2005, and 4.9% per year between 2005 and 2013. Trends in flow-normalized partition coefficients for both Hg and MeHg suggest a statistically significant shift from the particulate to the filtered phase. The upwardly accelerating f.THg concentration and observed shift from the solid phase to the aqueous phase among the pools of Hg and MeHg within the river water column signals an increased risk of deteriorating ecological conditions in the lower basin with respect to Hg contamination. More broadly, the 16-year trend analysis, completed 140 years after the commencement of major Hg releases to the Carson River, provides a poignant example of the ongoing legacy left behind by gold and silver mining techniques that relied on Hg amalgamation, and a cautionary tale for regions still pursuing the practice in other countries.

Nevada↗

MT3D-USGS version 1: A U.S. Geological Survey release of MT3DMS updated with new and expanded transport capabilities for use with MODFLOW

MT3D-USGS, a U.S. Geological Survey updated release of the groundwater solute transport code MT3DMS, includes new transport modeling capabilities to accommodate flow terms calculated by MODFLOW packages that were previously unsupported by MT3DMS and to provide greater flexibility in the simulation of solute transport and reactive solute transport. Unsaturated-zone transport and transport within streams and lakes, including solute exchange with connected groundwater, are among the new capabilities included in the MT3D-USGS code. MT3D-USGS also includes the capability to route a solute through dry cells that may occur in the Newton-Raphson formulation of MODFLOW (that is, MODFLOW-NWT). New chemical reaction Package options include the ability to simulate inter-species reactions and parent-daughter chain reactions. A new pump-and-treat recirculation package enables the simulation of dynamic recirculation with or without treatment for combinations of wells that are represented in the flow model, mimicking the above-ground treatment of extracted water. A reformulation of the treatment of transient mass storage improves conservation of mass and yields solutions for better agreement with analytical benchmarks. Several additional features of MT3D-USGS are (1) the separate specification of the partitioning coefficient ( K d ) within mobile and immobile domains; (2) the capability to assign prescribed concentrations to the top-most active layer; (3) the change in mass storage owing to the change in water volume now appears as its own budget item in the global mass balance summary; (4) the ability to ignore cross-dispersion terms; (5) the definition of Hydrocarbon Spill-Source Package (HSS) mass loading zones using regular and irregular polygons, in addition to the currently supported circular zones; and (6) the ability to specify an absolute minimum thickness rather than the default percent minimum thickness in dry-cell circumstances. Benchmark problems that implement the new features and packages test the accuracy of new code through comparison to analytical benchmarks, as well as to solutions from other published codes. The input file structure for MT3D-USGS adheres to MT3DMS conventions for backward compatibility: the new capabilities and packages described herein are readily invoked by adding three-letter package name acronyms to the name file or by setting input flags as needed. Memory is managed in MT3D-USGS using FORTRAN modules in order to simplify code development and expansion.

Techniques and Methods↗

Toward improved simulation of river operations through integration with a hydrologic model

Advanced modeling tools are needed for informed water resources planning and management. Two classes of modeling tools are often used to this end–(1) distributed-parameter hydrologic models for quantifying supply and (2) river-operation models for sorting out demands under rule-based systems such as the prior-appropriation doctrine. Within each of these two broad classes of models, there are many software tools that excel at simulating the processes specific to each discipline, but have historically over-simplified, or at worse completely neglected, aspects of the other. As a result, water managers reliant on river-operation models for administering water resources need improved tools for representing spatially and temporally varying groundwater resources in conjunctive-use systems. A new tool is described that improves the representation of groundwater/surface-water (GW-SW) interaction within a river-operations modeling context and, in so doing, advances evaluation of system-wide hydrologic consequences of new or altered management regimes.

Environmental Modelling and Software↗