Search NASA⌕ Search

SEARCH · Search NASA

Results for “Groundwater model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

Quantifying the Impact of Lagged Hydrological Responses on the Effectiveness of Groundwater Conservation

Many irrigated agricultural areas seek to prolong the lifetime of their groundwater resources by reducing pumping. However, it is unclear how lagged responses, such as reduced groundwater recharge caused by more efficient irrigation, may impact the long-term effectiveness of conservation initiatives. Here, we use a variably saturated, simplified surrogate groundwater model to: 1) analyze aquifer responses to pumping reductions, 2) quantify time lags between reductions and groundwater level responses, and 3) identify the physical controls on lagged responses. We explore a range of plausible model parameters for an area of the High Plains Aquifer (USA) where stakeholder-driven conservation has slowed groundwater depletion. We identify two types of lagged responses that reduce the long-term effectiveness of groundwater conservation, recharge-dominated and lateral-flow-dominated, with vertical hydraulic conductivity (K Z ) the major controlling variable. When high K z allows percolation to reach the aquifer, more efficient irrigation reduces groundwater recharge. By contrast, when low K z impedes vertical flow, short term changes in recharge are negligible, but pumping reductions alter the lateral flow between the groundwater conservation area and the surrounding regions (lateral-flow-dominated response). For the modeled area, we found that a pumping reduction of 30% resulted in median usable lifetime extensions of 20 or 25 years, depending on the dominant lagged response mechanism (recharge- vs. lateral-flow-dominated). These estimates are far shorter than estimates that do not account for lagged responses. To conclude, results indicate that conservation-based pumping reductions can extend aquifer lifetimes, but lagged responses can create a sizable difference between the initially perceived and actual long-term effectiveness.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of Timber Mountain Recharge and Groundwater Flow in Relation to Pahute Mesa-Oasis Valley Flow System

The purpose of this study was to evaluate the effect of Timber Mountain recharge on groundwater flow in the surrounding Pahute Mesa-Oasis Valley groundwater basin using several different methodologies. These methodologies included examining climatic and soil moisture data, estimating groundwater recharge by chloride mass balance, interpreting water chemical and isotopic data, and modeling groundwater flow along faults.

54 ENVIRONMENTAL SCIENCES↗

Sensor Recommendations for Long Term Monitoring of the F-Area Seepage Basins

In mid-2018, a new paradigm for long-term monitoring was developed after of decade of applied research projects funded by the Department of Energy’s office of Environmental Management Technology Development program. The program at SRNL was focused on transitioning complex environmental waste sites from active to passive remediations strategies. A key result of these studies was that the use of enhanced attenuation approaches at radiologically contaminated sites will result in the creation of secondary source areas in the subsurface that will require monitoring for decades. Alternative monitoring approaches are being developed and tested at the Savannah River Site’s F-Area Hazardous Waste Management Facility, the new paradigm provides innovative solutions that will significantly lower costs of monitoring through the coupling of data collection, machine learning and deterministic groundwater modeling. The foundation of this approach is a well-optimized network of sensors for measuring hydrogeochemical master variables that control, and therefore act as indicators of groundwater contaminant transport. By monitoring changes in the controlling master variables over time arising from geological and environmental shifts, predictive modelling can assist with identifying new strategies for ensuring regulatory requirements are met if trends toward conditions for potential remobilization of attenuated contaminants are detected. In this report, we evaluated commercially available single parameter sensor platforms (e.g., temperature/depth) and configurable multi-parameter sensor platforms (e.g., pH, oxidation-reduction potential, temperature, depth, dissolved oxygen, and conductivity). Each was scored using an optimization function based on how well the system supports the proposed long-term monitoring paradigm, in general, and the site-specific conditions at F-Area, in particular. Several viable sensor systems were identified. Of these, a combined platform including the In-Situ Aqua TROLL 500 multi-parameter sensor platform and the In-Situ temperature/depth sensor had the highest rating and was identified as the most suitable candidate for installation and monitoring of the master variables and potentiometric surface that control groundwater contaminant plumes emanating from the F-Area Seepage Basins. The discussion of recommended potential deployment locations builds upon recommendations made by Denham et al (2019).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Numerical investigations to identify environmental factors for field-scale reactive transport of pathogens at riverbank filtration sites

While induced bank filtration is a proven method for facilitating sustainable drinking water production, it is at risk from surface water contaminations (e.g., pathogens). Induced bank filtration and pathogen transport in groundwater have been studied extensively. However, long-term studies that consider real-world conditions are missing. These conditions include seasonal changes to environmental conditions and waterworks operations. Therefore, to analyze the effect of seasonal changes on the transport of human pathogenic viruses and their indicators in induced bank filtration, concentrations of adenoviruses and pathogen indicators were monitored over 16 months at an active bank filtration plant at the Rhine River, in Düsseldorf (Germany). Based on this data, a 2D groundwater model was created in PFLOTRAN that simulated flow, heat transport, conservative transport of chloride and the resulting electrical conductivity, reactive transport of oxygen and nitrate, and colloid-based transport of coliforms, somatic coliphages, and adenoviruses. The results show that reduced travel time was the key factor determining periods with a low removal of coliforms and somatic coliphages in the aquifer. Travel time was controlled by river level variations during rainy seasons, and the waterworks extraction rates during dry seasons. Further, for adenovirus transport, travel distance in the subsurface appeared to be the key factor, while travel time had no significant impact. Coliform removal increased when the colmation layer permeability decreased, while coliphage and adenovirus removal was unaffected by the colmation layer permeability. Seasonal changes in temperature and oxygen content did not significantly impact the removal of coliphages and adenoviruses in groundwater. Denitrifying conditions correlated with a lowered coliform removal, but the modelling could not establish a connection between denitrifying conditions and coliform removal. Our study showed that removal of pathogens and pathogen indicators at induced bank filtration plants varies greatly in time and space (e.g., for coliforms from 1 to 4 log-levels at 20m travel distance), and that adenovirus transport differs considerably from transport of coliforms and somatic coliphages.

59 BASIC BIOLOGICAL SCIENCES↗

MODFLOW6 models used to evaluate potential stresses and hydrologic conditions driving water-level fluctuations in well ER-5-3-2, Frenchman Flat, Southern Nevada

The hydrograph for well ER-5-3-2 in Frenchman Flat, southern Nevada, has previously unexplained water-level fluctuations. Four, three-dimensional, groundwater models (MODFLOW 6) were developed to evaluate potential stresses and hydrologic conditions affecting the well ER-5-3-2 hydrograph. Four model scenarios were developed that simulated: (1) wellbore leakage without recharge, (2) wellbore leakage with recharge, (3) shallow (low transmissivity) and deep (high transmissivity) carbonate rocks, and (4) lateral heterogeneity of carbonate rocks. Input and output files for the four model scenarios are in the model and output directories, respectively. Hydraulic conductivity, specific storage, and wellbore-leakage rates (when simulated) were estimated with parameter estimation (PEST) by minimizing a weighted composite, sum-of-squares objective function. The objective function was informed by measurement and Tikhonov regularization observations. Measurement observations included drawdowns from the constant-rate aquifer test and water-level altitudes measured in well ER-5-3-2 from 2001-2021. Tikhonov regularization informed hydraulic conductivity and specific storage parameters that were insensitive to measurement observations, where homogeneity was the preferred relation. Batch files, executables, and MODFLOW 6, PEST, and post-processing utilities are in the ancillary directory. Supplementary data also are included in the ancillary directory, including site information, high-frequency water-level and aquifer-test data, transmissivity estimates, water-chemistry data, and water-temperature analyses. This USGS data release contains data, analyses, and model files for the simulations and analysis results described in U.S. Geological Survey Scientific Investigations Report (https://doi.org/10.3133/sir20225132).

54 ENVIRONMENTAL SCIENCES↗

Application of the Recharge Estimation Tool (RET) to Prepare Spatially and Temporally Variable Recharge Boundary Conditions for Hanford Site Composite Analysis Vadose Zone Models

This environmental calculation file (ECF) describes the development of a tool for translating recharge estimates into readable input for STOMP© (Subsurface Transport Over Multiple Phases) vadose zone models primarily supporting the vadose zone (VZ) facets of the updated Hanford Site Composite Analysis (CA) and the Hanford Site Cumulative Impact Evaluation (CIE). The recharge estimates are spatiotemporally variable and are produced by the Recharge Evolution Tool (RET) described in Hanford Site-wide Natural Recharge Boundary Condition for Groundwater Models (ECF-HANFORD-15-0019). Outputs from the RET are given in Esri’s™ feature class format with yearly estimates and associated metadata encapsulated in file geodatabase objects. For STOMP models, the translated output is a text file in the format of an input boundary condition card, consistent with STOMP software requirements. The text file contains assimilated spatiotemporal recharge estimates produced by the RET. The tool discussed in this document will be referred to as the “RET2STOMP” tool.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Bayesian Approach to Estimation of Water Table Elevations Using Historical Rasters as Prior Information 2019 - 20430

In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)

54 ENVIRONMENTAL SCIENCES↗

Carbon Tetrachloride Degradation Results for 200-ZP-1 Operable Unit

Carbon tetrachloride (CT) contamination in the 200-ZP-1 Operable Unit (OU) at the Hanford Site originated from large-volume discharges to the subsurface during plutonium production operations between 1955 and 1973. Contamination migrated through more than 70 meters of unsaturated sediment to reach the underlying unconfined aquifer, where it persists as a large and complex groundwater plume. The 200-ZP-1 OU Record of Decision (ROD) requires that groundwater CT concentrations be reduced to 3.4 µg/L within 125 years. Current groundwater modeling projections estimate that the existing pump-and-treat, even when combined with monitored natural attenuation (specifically hydrolysis), will not achieve this target within the designated timeframe. A fundamental contributor to this shortfall is the extremely slow rate of CT hydrolysis under Hanford aquifer conditions, which has been estimated to have a half-life of 630 years. If faster-acting biotic and abiotic degradation processes are operating within the aquifer, their contribution to CT mass reduction could have a meaningful impact. However, site-specific measurements of these processes and their rates have not previously been performed. This report documents the results of a two-phase laboratory investigation designed to characterize and quantify the capacity of site-specific 200-ZP-1 OU sediments and groundwater to support natural attenuation of CT through biotic and abiotic pathways. In this context, degradation capacity is defined as the intrinsic potential of the subsurface matrix to transform CT under optimized, controlled conditions. System capacity is evaluated in two ways: (1) as rate-limited capacity, which establishes the maximum kinetic velocity of CT transformation and is measured using half-lives and first order rate constants; and (2) as mass limited capacity, which defines the total contaminant mass the batch experimental system can degrade before reactants are exhausted, representing the maximum amount of contaminant the microbial community and reactive mineral phases can transform under the experimental conditions.

abiotic degradation↗

3D hydrogeophysical characterization of managed aquifer recharge basins

Aquifers are increasingly stressed. Managed aquifer recharge provides a potential solution to mitigate this stress and provide sustainable groundwater resources. Subsurface properties are known to have a strong control on the infiltration rates that can be achieved. However, these properties are often highly heterogeneous and difficult to assess with conventional probing techniques. Here, we show the application of 3D geophysical imaging to assess the recharge potential and its variation across several basins used for managed aquifer recharge. We link in-situ measurements of saturated hydraulic conductivity with the electrical resistivity of the subsurface to establish petrophysical relationships and use those relationships to estimate the distribution of hydraulic conductivity throughout the five recharge basins. Our results show a considerable variability in the hydraulic properties, i.e., soil texture and saturated hydraulic conductivity, that have a direct impact on potential infiltration rates. We use the 3D hydraulic property distributions to model groundwater recharge and provide estimates for infiltration rates and volumes, and use this approach to assess the impact of management activities on groundwater recharge performance. Having such data not only enables us to predict infiltration rates, but also provides means for optimizing such water infrastructure.

54 ENVIRONMENTAL SCIENCES↗

Characterization of C-Reactor (105-C) Moderator Tanks 204 and 205

Area Completion Project (ACP)/Savannah River Nuclear Solutions (SRNS) has requested Savannah River National Laboratory (SRNL) to perform chemical and radiological analyses on moderator heavy water samples from C-Reactor (105-C) in Moderator Tanks 204 and 205. The heavy water samples have been characterized using SRNL analytical methods. Information generated from the characterization of these heavy water samples provide a basis for determining the residual contamination remaining in the tanks. These analyses are needed in order to update the C-Reactor Contaminant migration groundwater model and report, to help determine a disposition pathway for the moderator water, and to ship samples off-site to the Southwest Research Institute (SWRI). This report presents characterization results for the March 2024 C-Reactor Moderator Tank 204 and 205 heavy water samples. SRNL results are critical to allow ACP to complete this Performance Based Incentive project on time. Based on the characterization data provided in this report, results of these samples meet customer expectations and satisfy the Technical Assistance Request (TAR).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Groundwater Volume Loss in Mexico City Constrained by InSAR and GRACE Observations and Mechanical Models

Abstract Groundwater withdrawal can cause localized and rapid poroelastic subsidence, spatially broad elastic uplift of low amplitude, and changes in the gravity field. Constraining groundwater loss in Mexico City, we analyze data from the Gravity Recovery and Climate Experiment and its follow‐on mission (GRACE/FO) and Synthetic Aperture Radar (SAR) Sentinel‐1A/B images between 2014 and 2021. GRACE/FO observations yield a groundwater loss of 0.85–3.87 km 3 /yr for a region of ∼300 × 600 km surrounding Mexico City. Using the high‐resolution interferometric SAR data set, we measure >35 cm/yr subsidence within the city and up to 2 cm/yr of uplift in nearby areas. Attributing the long‐term subsidence to poroelastic aquifer compaction and the long‐term uplift to elastic unloading, we apply respective models informed by local geology, yielding groundwater loss of 0.86–12.57 km 3 /yr. Our results suggest Mexico City aquifers have been depleting at faster rates since 2015, exacerbating the socioeconomic and health impacts of long‐term groundwater overdrafts.

58 GEOSCIENCES↗

Integrating Analytical Solutions and U-Net Model for Predicting Groundwater Contaminant Plumes in Pump-and-Treat Systems

Pump-and-treat (P&T) is a common technique for groundwater remediation involving the extraction and treatment of contaminated water above ground. Optimizing the design and operation of the P&T well network is essential for maximizing the system’s effectiveness and efficiency. However, this optimization often necessitates many model evaluations, leading to computationally demanding tasks. This study introduces a novel approach that integrates analytical solutions for groundwater dynamics with the U-Net (Ronneberger et al., 2015) deep learning framework to predict groundwater contaminant plume migration under dynamic pumping conditions. By incorporating the Thiem equation (Thiem, 1906) into the input preprocessing, the U-Net model transforms sparse well data into a continuous spatial field that captures the hydraulic impacts of pumping activities. This integration enables the model to leverage both deep learning capabilities and classical physics-based groundwater theories, enhancing prediction accuracy and computational efficiency. These advancements can facilitate rapid, large-scale evaluations of P&T optimization simulations, allowing for timely and effective decision-making in well placement and system management. We demonstrate the model's robust performance across both simplified transient 2D models and a more complex 3D heterogeneous site model at the 200 West P&T facility at the Hanford Site. The U-Net-based model offers substantial computational advantages, reducing simulation times significantly compared to full physics-based models and providing a powerful tool for rapid site evaluation and P&T system optimization, such as evaluating alternative P&T well network designs. Our findings highlight the potential of advanced machine learning models to significantly enhance the efficiency and sustainability of groundwater remediation efforts, offering a novel application of U-Net architecture in environmental science.

Pump-and-treat↗

Verification and Validation of a Modified Numerical Algorithm for Simulation of Transient Unconfined Groundwater Flow

This report extends verification and validation testing of the NUFT package of codes (Nitao 2000a,b) to the US1P module of code designed for simulation of single-phase mixed or “coupled” saturated and unsaturated (or “variably saturated”) groundwater flow. Importantly, the verification and validation testing of US1P in this report includes performance evaluation of a modified numerical algorithm that was not included in prior quality assurance (QA) of the NUFT package by Carle et al. (2014) for Underground Test Area (UGTA) project activities, which are directed at assessment of radionuclide contamination in groundwater sourced from underground nuclear test locations at the Nevada National Security Site (NNSS). An immediate purpose of this report is to provide QA for ongoing large-scale, three-dimensional (3-D) groundwater flow modeling of transient water levels associated with longterm water supply pumping and underground nuclear testing at the NNSS (Jackson and Fenelon, 2018; Jackson et al., 2021). For UGTA activities, the QA refers to the standards of ASTM (1996) for verification and validation testing of groundwater modeling codes.

54 ENVIRONMENTAL SCIENCES↗

Groundwater and Surface Water Flow (GSFLOW) model files to explore bedrock circulation depth and porosity in Copper Creek, Colorado

This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.

54 ENVIRONMENTAL SCIENCES↗

Effects of groundwater pumping on ground surface temperature: A regional modeling study in the North China Plain

Over-exploitation of groundwater (GW) in the North China Plain (NCP) since the 1960s has many environmental consequences. However, understanding of the dominant mechanisms remains limited, particularly at regional scale. In this study, the coupled ParFlow.CLM model representing subsurface and land-surface processes and their interactions was applied in the NCP at high spatio-temporal resolutions. The model was validated using the water and energy fluxes reported in previous studies and from the JRA-55 reanalysis. Numerical experiments were designed to examine the relative impacts of GW pumping and irrigation on the ground surface temperature (GST). Results show significant effects of GW pumping on GST in the NCP. Generally, the subsurface acts as a buffer to temporal variations in heat fluxes at the land-surface, but long-term pumping can gradually weaken this buffer, resulting in increases in the spatio-temporal variability of GST, as exemplified by hotter summers and colder winters. Considering that changes of water table depth (WTD) can significantly affect land surface heat fluxes when WTD ranges roughly between 0.01–10 m, the 0.5 m/year increase of WTD simulated by the model due to pumping can continue to increase the average WTD and hence, disturb GST, for about 20 years from the pre-pumping WTD in the NCP, before the WTD exceeds 10 m. The variations of GST are expected to be faster initially and gradually slow down due to the nonlinear behaviors of GST with WTD. The findings from this study in the NCP may also have implications for other regions with GW depletion.

Yang, Chen↗

Contaminant Transport Parameters for the Groundwater Flow and Contaminant Transport Model of Corrective Action Units 101 and 102: Central and Western Pahute Mesa, Nye County, Nevada, Revision 1

This document presents a summary and framework of available transport data and other information directly relevant to the development of the groundwater transport model for Central and Western Pahute Mesa (PM) corrective action units (CAU). This model is referred to as the “PM CAU transport model” in this document. Central and Western PM are two areas of the Nevada National Security Site (NNSS) that historically were used for underground nuclear testing. These nuclear tests resulted in groundwater contamination in the vicinity of the underground test areas. As a result, the U.S. Department of Energy (DOE), Environmental Management (EM) Nevada Program is currently conducting a corrective action investigation (CAI) of the PM underground test areas. This work is a part of the Underground Test Area (UGTA) Activity in accordance with the Federal Facility Agreement and Consent Order (FFACO) (1996, as amended).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Vadose Zone Model for U-10 West Area for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3 (CP-57037, Model Package Report: Plateau to River Groundwater Transport Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort: carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Groundwater and Surface Water Flow (GSFLOW) model files for the East River, Colorado

The data package contains model input files and executables for the East River, Colorado (750 km2) located in the headwaters of the Upper Colorado River Basin. The code applied is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution to represent topographic complexity and a daily timestep accounts for energy and water partitioning between the snowpack, vegetation, soil zone and bedrock. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model accounts for nine stratigraphic units and extends 400 m below land surface. Using this modeling framework we explore historical water budgets (water year 1987-2022) and the influence of seasonal warming and associated mechanisms driving groundwater declines and streamflow loss. A Readme_062424.txt file provides instructions on how to download all files. Input and output files are provided for the historical simulation representative of water years 1987 to 2022 (i.e. baseline) and the all-year +4C warming scenario. Instructions are provided to run the seasonal warming scenarios with warming applied only to the autumn, winter, spring or summer months. Modeled output used in figures for Carroll et al., 2024 are also provided with metadata describing where these data were obtained. This dataset additionally includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type.Updated on 06-24-2024: The dataset was updated to include model outputs (updated base.zip and 4C.zip files), figure source files (Figure_Source_Files.zip), an updated Readme file (Readme_062424.txt), climate input files (climate_input_files.zip), and the adoption of ESS-DIVE File Level Metadata and CSV reporting formats (inclusion of flmd.csv and dd.csv files).

54 ENVIRONMENTAL SCIENCES↗