Richard Gray Niswonger

The Santa Rosa Plain Hydrologic Model (SRPHM) was developed and published in 2014 through a collaboration between the U.S. Geological Survey (USGS) and Sonoma Water to analyze the hydrologic system in the Santa Rosa Plain watershed, help meet the increasing demand for fresh water, and prepare for future uncertainties in water resources. The original model simulated hydrological conditions and water use from water years 1975 to 2010. Recently (2023), the USGS, in cooperation with Sonoma Water and the California State Water Resources Control Board, updated the SRPHM model to extend its simulation period to the end of the 2018 calendar year, incorporate new estimates of rural and agricultural water use, and use efficient input format for climate variables. The updated model was recalibrated, and evaluation of the new model calibration is included in this report. This report presents the results of comparing the hydraulic heads, streamflow, and groundwater budget simulated by the updated model with those generated by the original model and observed data. The main difference in the simulated budget between the original and updated SRPHM is the estimates of agricultural pumping, rural domestic pumping, and return flow generated from rural water use that was not simulated in the original model. The revised agricultural pumping is simulated using the agricultural package, which constrains pumping to available groundwater. The use of the agricultural package leads to a more realistic estimation of agricultural water use, with revised agricultural pumping being one-third less than that in the original model. The revised rural pumping is about half of the pumping in the original model because of using detailed parcel data to estimate population density in rural areas instead of coarse census tracts. Overall, average total inflows for water years 2006–10 simulated by the updated model were about 2 percent less than the original model, and the average total updated outflows were nearly 5 percent less than the original model. The updated model was then used to generate stream depletion maps, simulate climate change scenarios during 2019–99, and simulate water rights allocation using the Model for Decision Support in Integrated River Basin Management (MODSIM). The results from simulating eight future climate scenarios indicated either an increase in groundwater storage or no significant change in the next 80 years, along with an increase in recharge, an increase in actual evapotranspiration in six out of eight climate projections, and an increase in surface runoff. The increases in the simulated future groundwater storage, recharge, evapotranspiration, and runoff in most climate projections are mainly driven by the projected increase in precipitation in most of the future climate scenarios. The updated model also was used to test a pilot case study demonstrating water-resource allocation among different users with different water rights using the integrated MODSIM-Groundwater and Surface-Water Flow Model (GSFLOW) platform. The updated SRPHM serves as a valuable tool for analyzing historical and future hydrologic conditions in the Santa Rosa Plain watershed and preparing for future uncertainties.

Estimation of human water withdrawals is more important now than ever due to uncertain water supplies, population growth, and climate change. Fourteen percent of the total water withdrawal in the United States is used for public supply, typically including deliveries to domestic, commercial, and occasionally including industrial, irrigation, and thermoelectric water withdrawal. Stewards of water resources in the USA require estimates of water withdrawals to manage and plan for future demands and sustainable water supplies. This study compiled the most comprehensive conterminous United States water withdrawal data set to date and developed a machine learning framework for estimating public supply withdrawals and associated uncertainty for the period 2000–2020. The modeling approach provides service area resolution estimates to allow for annual and monthly water withdrawal estimation while incorporating a complex array of driving factors that include hydroclimatic, demographic, socioeconomic, geographic, and land use factors. Model results reveal highly variable and lognormally distributed per‐capita water withdrawal, spanning from 30 to 650 gallons per capita per day (GPCD), across community, regional, and national scales, with pronounced seasonal variations. Analysis of estimated withdrawal trends indicates that the national annual average withdrawal experienced a decline at a rate of 0.58 GPCD/year during the period from 2000 to 2020. Model interpretation reveals a complex interplay between public supply withdrawal and key predictors, including population size, warm‐season precipitation, counts of large buildings and houses, and areas of urban and commercial land use. The developed models can forecast future public supply driven by various climate, demographic, and socioeconomic scenarios.

Groundwater interactions with mountain streams are often simplified in model projections, potentially leading to inaccurate estimates of streamflow response to climate change. Here, using a high-resolution, integrated hydrological model extending 400 m into the subsurface, we find groundwater an important and stable source of historical streamflow in a mountainous watershed of the Colorado River. In a warmer climate, increased forest water use is predicted to reduce groundwater recharge resulting in groundwater storage loss. Losses are expected to be most severe during dry years and cannot recover to historical levels even during simulated wet periods. Groundwater depletion substantially reduces annual streamflow with intermittent conditions predicted when precipitation is low. Expanding results across the region suggests groundwater declines will be highest in the Colorado Headwater and Gunnison basins. Our research highlights the tight coupling of vegetation and groundwater dynamics and that excluding explicit groundwater response to warming may underestimate future reductions in mountain streamflow.

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.

Hydrological models are usually calibrated to in‐situ streamflow observations with reasonably long and uninterrupted records. This is challenging for poorly gage or ungaged basins where such information is not available. Even for gaged basins, the single‐objective calibration to gaged streamflow cannot guarantee reliable forecasts because, as has been documented elsewhere, the inverse problem is mathematically ill‐posed. Therefore, the inclusion of other observations, and the reproduction of other hydrological variables beyond streamflow, become critical components of accurate hydrological forecasting. In this study, six single‐ and multi‐objective model calibration schemes based on different combinations of gaged streamflow, global‐scale gridded soil moisture, actual evapotranspiration (ET), and runoff products are used for the calibration of a process‐based hydrological model for 20 catchments located within the Lake Michigan watershed, of the Laurentian Great Lakes. Results show that the addition of gridded soil moisture to gaged streamflow in model calibration improves the ET simulation performance for most of the catchments, leading to the overall best‐performing models. The monthly streamflow simulation performance for the experiments using gridded runoff products to inform the model is outperformed by those using the gaged streamflow, but the discrepancy is mitigated with increasing catchment scale. A new visualization method that effectively synthesizes model performance for the simulations of streamflow, soil moisture, and ET was also proposed. Based on the method, it is revealed that the streamflow simulation performance is relatively weak for baseflow‐dominated catchments; overall, the 20 catchment models simulate streamflow and ET better than soil moisture.

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.

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.

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 ⁻² and 10 ⁻³ 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.

Deeper flows through bedrock in mountain watersheds could be important, but lack of data to characterize bedrock properties limits understanding. To address data scarcity, we combine a previously published integrated hydrologic model of a snow‐dominated, headwater basin of the Colorado River with a new method for dating baseflow age using dissolved gas tracers SF₆, CFC‐113, N₂, and Ar. The original flow model predicts the majority of groundwater flow through shallow alluvium (6 concentrations. To match gas data, bedrock permeability is increased to allow a larger fraction of deeper and older groundwater flow (median 112 m). The updated hydrologic model indicates interannual variability in baseflow age (3–12 years) is controlled by the volume of seasonal interflow and tightly coupled to snow accumulation and monsoon rain. Deeper groundwater flow remains stable (11.7 ± 0.7 years) as a function mean historical recharge to bedrock hydraulic conductivity (R/K). A sensitivity analysis suggests that increasing bedrock K effectively moves this alpine basin away from its original conceptualization of hyperlocalized groundwater flow (high R/K) with groundwater age insensitive to changes in water inputs. Instead, this basin is situated close to the precipitation threshold defining recharge controlled groundwater flow conditions (low R/K) in which groundwater age increases with small reductions in precipitation. Work stresses the need to explore alternative methods characterizing bedrock properties in mountain basins to better quantify deeper groundwater flow and predict their hydrologic response to change.

Deeper flows through bedrock in mountain watersheds could be important but lack of data to characterize bedrock properties and link flow paths to snow-dynamics limits understanding. To address data scarcity, we combine a previously published integrated hydrologic model of a snow-dominated, headwater basin with a new method for dating baseflow age using dissolved gas tracers SF, N, Ar. The original flow model produces shallow groundwater flow (median depth 6 m), very young stream water and is unable to reproduce observed SF concentrations. To match the observed gas data, bedrock permeability is increased to allow a larger fraction of deeper groundwater flow (median depth 110 m). Results indicate that interannual variability in baseflow age (3-12 y) is dictated by the volume of seasonal interflow. Deeper groundwater flow remains stable (11.7±0.7 y) as a function of the ratio of recharge to bedrock hydraulic conductivity (R/K), where recharge is dictated by long-term climate and land use. With sensitivity experiments, we show that information gleaned from gas tracer data to increase bedrock hydraulic conductivity effectively moves this alpine basin away from shallow, topographically controlled groundwater flow with baseflow age relatively insensitive to water inputs (high R/K), and closer toward recharge-controlled conditions, in which a small shift toward a drier future with less snow accumulation will alter the groundwater flow system and increase baseflow age (low R/K). Work stresses the need to explore alternative methods characterizing bedrock properties in mountain basins to better quantify deeper groundwater flow and predict their hydrologic response to change.