Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “water table depth”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Development of inter-grid-cell lateral unsaturated and saturated flow model in the E3SM Land Model (v2.0)

Abstract. The lateral transport of water in the subsurface is important in modulating terrestrial water energy distribution. Although a few land surface models have recently included lateral saturated flow within and across grid cells, it is not a default configuration in the Climate Model Intercomparison Project version 6 experiments. In this work, we developed the lateral subsurface flow model within both unsaturated and saturated zones in the Energy Exascale Earth System Model (E3SM) Land Model version 2 (ELMv2.0). The new model, called ELMlat, was benchmarked against PFLOTRAN, a 3D subsurface flow and transport model, for three idealized hillslopes that included a convergent hillslope, divergent hillslope, and tilted V-shaped hillslope with variably saturated initial conditions. ELMlat showed comparable performance against PFLOTRAN in terms of capturing the dynamics of soil moisture and groundwater table for the three benchmark hillslope problems. Specifically, the mean absolute errors (MAEs) of the soil moisture in the top 10 layers between ELMlat and PFLOTRAN were within 1 %±3 %, and the MAEs of water table depth were within ±0.2 m. Next, ELMlat was applied to the Little Washita experimental watershed to assess its prediction of groundwater table, soil moisture, and soil temperature. The spatial pattern of simulated groundwater table depth agreed well with the global groundwater table benchmark dataset generated from a global model calibrated with long-term observations. The effects of lateral groundwater flow on the energy flux partitioning were more prominent in lowland areas with shallower groundwater tables, where the difference in simulated annual surface soil temperature could reach 0.3–0.4 ∘C between ELMv2.0 and ELMlat. Incorporating lateral subsurface flow in ELM improves the representation of the subsurface hydrology, which will provide a good basis for future large-scale applications.

58 GEOSCIENCES↗
DOE OSTI2025

Subsurface hydrological controls on the short-term effects of hurricanes on nitrate–nitrogen runoff loading: a case study of Hurricane Ida using the Energy Exascale Earth System Model (E3SM) Land Model (v2.1)

When the nutrient level in the soil surpasses vegetation demand, nutrient losses due to surface runoff and subsurface leaching are the major reasons for the deterioration of water quality. The lower Mississippi River basin (LMRB) is one of the sub-basins that deliver the highest nitrogen loads to the Gulf of Mexico. Potential changes in episodic events induced by hurricanes may exacerbate water quality issue in the future. However, uncertainties in modeling the hydrologic response to hurricanes may limit the modeling of nutrient losses during such events. Using a machine learning approach, we calibrated the land component of the Energy Exascale Earth System Model (E3SM), or ELM, version 2.1, based on the water table depth (WTD) of a calibrated 3D subsurface hydrology model. While the overall performance of the calibrated ELM is satisfactory, some discrepancies in WTD remain in slope areas with low precipitation due to the missing lateral flow process in ELM. Simulations including biogeochemistry performed using ELM with and without model calibration showed important influences of soil hydrology, precipitation intensity, and runoff parameterization on the magnitude of nitrogen runoff loss and the leaching pathway. Despite such sensitivities, both ELM simulations produced reduced WTD and increased runoff and accelerated nitrate–nitrogen runoff loading during Hurricane Ida in August 2021, consistent with the observations. With observations suggesting more pronounced effects of Hurricane Ida on nitrogen runoff than the simulations, we identified factors for model improvement to provide a useful tool for studying hurricane-induced nutrient losses in the LMRB region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Bedrock weathering rates, reactive nitrogen influxes and effluxes, and nitrous oxide emissions rates from the Pumphouse Hillslope, East River Watershed, Colorado

Atmospheric nitrous oxide contributes directly to global warming, yet models of the nitrogen cycle do not account for bedrock, the largest pool of terrestrial nitrogen, as a source of nitrous oxide. Although it is known that release rates of nitrogen from bedrock are large, there is an incomplete understanding of the connection between bedrock-hosted nitrogen and atmospheric nitrous oxide. Here, we quantify nitrogen fluxes and mass balances at a hillslope underlain by marine shale. We found that at this site bedrock weathering contributes 78% of the subsurface reactive nitrogen, while atmospheric sources (commonly regarded as the sole sources of reactive nitrogen in pristine environments) account for only the remaining 22%. About 56% of the total subsurface reactive nitrogen denitrifies, including 14% emitted as nitrous oxide. The remaining reactive nitrogen discharges in porewaters to a floodplain where additional denitrification likely occurs. We also found that the release of bedrock nitrogen occurs primarily within the zone of the seasonally fluctuating water table and suggest that the accumulation of nitrate in the vadoes zone, often attributed to fertilization and soil leaching, may also include contributions from weathered nitrogen-rich bedrock. Our hillslope study suggests that under oxygenated and moisture-rich conditions, weathering of deep, nitrogen-rich bedrock makes an important contribution to the nitrogen cycle. The data files are in Excel, which can be accessed using Microsoft Office, and consist of many data sets from the Pumphouse Hillslope PLM (Pumphouse Lower Montane) 1, 2, 3, and 4. They include soil to rock (0-10 meters) solid phase minerals and elements compositions; time- and depth-resolved pore-water chemistry and pore-gas compositions; time- and depth-resolved water table depths and water fluxes; subsurface weathering rates; nitrogen influxes and effluxes and mass balance. The attached paper is in a Word document, which can be accessed using Microsoft Office, and in pdf format.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Soil and groundwater environmental sensor data, Wax Lake Delta, Louisiana, March 2023 - March 2024

This study evaluates how environmental parameters that integrate biogeochemical processes vary with water table fluctuations in the freshwater Wax Lake Delta (WLD) in Louisiana, U.S.A. This data package contains seven *.csv files and one Excel file that compiles all the data from the individual .csv files. This dataset reports high frequency (15-min) observations of water level, soil redox potential, specific conductance, and pH made for one year along elevation transects located on the older, proximal (OT) and younger, distal (YT) ends of a deltaic island. Water depth relative to the ground surface (cm; HOBO U20L-04; error ± 0.4 cm), water pH and temperature (HOBO MX2501), and specific conductance and temperature (HOBO U24-001) sensors were installed in March 2023. Water depth was corrected for barometric pressure recorded by a separate logger secured to a platform above the highest water level. Soil redox probes (SWAP ORP-40-4-B) were also installed in March 2023. Each probe had four Pt sensors (2 mm width) placed at 10 cm, 20 cm, 30 cm, and 40 cm below the ground surface. Redox data were referenced to an external Ag0/AgCl (3M KCl) reference probe placed in saturated ground and recorded on CR1000X dataloggers (Campbell Scientific) powered by solar panels. A second reference probe was positioned near the primary reference probe for backup and data correction. The tops of the soil redox probes and soil moisture probes were flush with the soil surface so that sensors are reported at their indicated depths below ground surface. Here, we report data collected between 15 March 2023 to 15 March 2024 for all sensors, with some differences due to exact dates of sensor placement or data gaps associated with sensor malfunction. For example, water depth at OT4 was not recorded between March to November 2023. Data flags indicate whether a value is valid (1) or was excluded from data analysis in the associated manuscript (-1).

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2024

Geophysical and Environmental Monitoring Data, and Subsurface Flow Modelling Results for Chicken Bone Meadow, Mt. Snodgrass, Crested Butte, CO

This dataset includes geoelectrical monitoring data acquired between October 2021 and November 2022, soil moisture and temperature data, groundwater data obtained from borehole SNIB covering the period from June 2021 to September 2022, and hydrological modelling results. The data were acquired to investigate how variations in bedrock type and topography, and vegetation cover control subsurface flow dynamics. To provide insights into the subsurface flow dynamics and their controls, a monitoring transect was installed at the Chicken Bone Meadow, Mt. Snodgrass, Crested Butte, CO, measuring the spatio-temporal variations of soil moisture, soil and snow temperature, subsurface electrical resistivity variations, and groundwater dynamics. Field data are organized in a folder structure, with Electrical Resistivity Tomography (ERT) data being provided as one file per measurement, and data of the soil moisture and temperature sensors being provided as text files covering the entire monitoring period. The ‘Locations.csv’ file contains the location of all sensors, given in NAD83 – UTM Zone 13N. ERT monitoring data has been processed to filter data based on reciprocal errors (data with errors > 30% were removed), a linear error model was fitted to each survey, and to ensure a constant set of measurements for time-lapse inversion, filtered data were interpolated and assigned a 100% measurement error. Soil moisture and temperature data were acquired at 15 min intervals, and averaged to provide 1h data. Weather data and borehole data (groundwater depth, conductivity and temperature) were acquired at 30 min intervals, and are provided as daily measurements; all measurements are averaged, except of precipitation values, which are given as daily accumulation. The hydrological model was set up along the ERT monitoring transect, and net infiltration was used as surface boundary condition and derived from the weather data. Four different results are provided, (1) results for a parameterization using hydraulic permeability and porosity as derived from the ERT data through petrophysical relationships, and (2) three simplified model results, using 1 to 3 geological layers above the bedrock. Modelling was performed using PFLOTRAN, and for each model the PFLOTRAN input files are provided. The result files include weekly hydrological modelling results (e.g., saturation, velocities, pressures), as well as the model parameterization. The 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.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Evaluating integrated surface/subsurface permafrost thermal hydrology models in ATS (v0.88) against observations from a polygonal tundra site

Numerical simulations are essential tools for understanding the complex hydrologic response of Arctic regions to a warming climate. However, strong coupling among thermal and hydrological processes on the surface and in the subsurface and the significant role that subtle variations in surface topography have in regulating flow direction and surface storage lead to significant uncertainties. Careful model evaluation against field observations is thus important to build confidence.We evaluate the integrated surface/subsurface permafrost thermal hydrology models in the Advanced Terrestrial Simulator (ATS) against field observations from polygonal tundra at the Barrow Environmental Observatory. ATS couples a multiphase, 3D representation of subsurface thermal hydrology with representations of overland nonisothermal flows, snow processes, and surface energy balance. We simulated thermal hydrology of a 3D ice-wedge polygon with geometry that is abstracted but broadly consistent with the surface microtopography at our study site. The simulations were forced by meteorological data and observed water table elevations in ice-wedge polygon troughs. With limited calibration of parameters appearing in the soil evaporation model, the 3-year simulations agreed reasonably well with snow depth, summer water table elevations in the polygon center, and high-frequency soil temperature measurements at several depths in the trough, rim, and center of the polygon. Upscaled evaporation is in good agreement with flux tower observations. The simulations were found to be sensitive to parameters in the bare soil evaporation model, snowpack, and the lateral saturated hydraulic conductivity.Timing of fall freeze-up was found to be sensitive to initial snow density, illustrating the importance of including snow aging effects. The study provides new support for an emerging class of integrated surface/subsurface permafrost simulators.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Five Years of Dissolved Oxygen, Temperature, Salinity, Depth, Weather Data from a Transitioning Wetland at Beaver Creek, Washington, USA

Groundwater dissolved oxygen (DO) variability in coastal system remains poorly understood despite its importance for biogeochemical cycling and ecosystem modeling. Here we investigate the temporal variability in groundwater DO and its hydro-climatic drivers across hourly to seasonal timescales in a transitioning wetland at Beaver Creek, Washington, USA. The site is transitioning from a freshwater forest to a brackish tidal wetland following removal of a barrier in 2014 that prevented tides from accessing the freshwater creek. By utilizing novel optical dissolved oxygen instrumentation (Opti O2, LLC) we obtained continuous, high-frequency (5-minute), in-situ measurements of DO from the flood-plain from June 26th, 2019 through September 30th, 2024. This 63 month dataset is comprised of groundwater dissolved oxygen, temperature, water level and salinity timeseries from the floodplain. This dataset also includes rainfall, air pressure, air temperature, and solar radiation data collected with a co-located Campbell ClimaVUE50 weather sensor. All data is contained within a single csv (2019-06-26 to 2024-09-30 Beaver Creek DO, saln, BGS, temp, weather.csv) that can easily be viewed either using software such as Excel or using any text editor.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Height Above Nearest Drainage (HAND) and Hydraulic Property Table for CONUS

The continental flood inundation mapping (CFIM) framework is a high-performance computing (HPC)-based computational framework for the Height Above Nearest Drainage (HAND)-based inundation mapping methodology. Using the 10m Digital Elevation Model (DEM) data produced by U.S. Geological Survey (USGS) 3DEP (the 3-D Elevation Program) and the NHDPlus hydrography dataset produced by USGS and the U.S. Environmental Protection Agency (EPA), a hydrological terrain raster called HAND is computed for HUC6 units in the conterminous U.S. (CONUS). The value of each raster cell in HAND is an approximation of the relative elevation between the cell and its nearest water stream. Derived from HAND, a hydraulic property table is established to calculate river geometry properties for each of the 2.7 million river reaches covered by NHDPlus (5.5 million kilometers in total length). This table is a lookup table for water depth given an input stream flow value. Such lookup is available between water depth 0m and 25m at 1-foot interval. The flood inundation map is then computed by using HAND and this lookup table based on the near real-time water forecast from the National Water Model (NWM) at the National Oceanic and Atmospheric Administration (NOAA). HAND and the Hydraulic Property Table version 0.2. is created to correspond to data updates in the USGS 1/3 arcsec DEM, the USGS National Hydrography Dataset, including its Water Boundary Dataset, and the NHDPlus medium resolution dataset. This dataset comprises 331 HUC6 units for CONUS (excluding the five great lakes units), each is a downloadable zip file. Version 0.1 was computed at the National Center for Supercomputing Applications (NCSA) at the University of Illinois at Urban-Champaign in 2016 and is currently hosted at the Texas Advanced Computing Center (TACC). Please cite this data DOI and the following publications and data DOIs when you use HAND and the hydraulic property table: Liu, Yan Y., David R. Maidment, David G. Tarboton, Xing Zheng, and Shaowen Wang. 'A CyberGIS integration and computation framework for high-resolution continental-scale flood inundation mapping.' JAWRA Journal of the American Water Resources Association 54, no. 4 (2018): 770-784. DOI: 10.1111/1752-1688.12660 Zheng, Xing, David G. Tarboton, David R. Maidment, Yan Y. Liu, and Paola Passalacqua. 'River channel geometry and rating curve estimation using height above the nearest drainage.' JAWRA Journal of the American Water Resources Association 54, no. 4 (2018): 785-806. DOI: 10.1111/1752-1688.12661 Liu, Y. (2018). Height Above Nearest Drainage (HAND) for CONUS - v0.1, HydroShare, https://doi.org/10.4211/hs.69f7d237675c4c73938481904358c789 LICENSE FOR USE -- MAPS AND DATA DISCLAIMER This resource is shared under the Creative Commons Attribution CC BY, http://creativecommons.org/licenses/by/4.0/ MAPS AND DATA DISCLAIMER The Oak Ridge National Laboratory (ORNL) shall not be held liable for improper or incorrect use of the data described or information contained on this map or associated series of maps. The data and related map graphics are not legal, land survey or engineering documents and are not intended to be used as such. ORNL gives no warranty, express or implied, as to the accuracy, reliability, utility or completeness of this information. The user of these maps and data assumes all responsibility and risk for the use of the maps and data. ORNL disclaims all warranties, representations or endorsements either express or implied, with regard to the information contained in this map product, including, but not limited to, all implied warranties of merchantability, fitness for a particular purpose or non-infringement. This preliminary map product is for research and review purposes only. It is not intended to be used for emergency management operational or life safety decisions at the local or regional governmental level or by the general public. Users requiring information regarding hazardous conditions or meteorological conditions for specific geographic areas should consult directly with their city or county emergency management office.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Height Above Nearest Drainage (HAND) and Hydraulic Property Table for CONUS - Version 0.21 (20200601)

The continental flood inundation mapping (CFIM) framework is a high-performance computing (HPC)-based computational framework for the Height Above Nearest Drainage (HAND)-based inundation mapping methodology. Using the 10m Digital Elevation Model (DEM) data produced by U.S. Geological Survey (USGS) 3DEP (the 3-D Elevation Program) and the NHDPlus hydrography dataset produced by USGS and the U.S. Environmental Protection Agency (EPA), a hydrological terrain raster called HAND is computed for HUC6 units in the conterminous U.S. (CONUS). The value of each raster cell in HAND is an approximation of the relative elevation between the cell and its nearest water stream. Derived from HAND, a hydraulic property table is established to calculate river geometry properties for each of the 2.7 million river reaches covered by NHDPlus (5.5 million kilometers in total length). This table is a lookup table for water depth given an input stream flow value. Such lookup is available between water depth 0m and 25m at 1-foot interval. The flood inundation map is then computed by using HAND and this lookup table based on the near real-time water forecast from the National Water Model (NWM) at the National Oceanic and Atmospheric Administration (NOAA). HAND and the Hydraulic Property Table version v0.21 is an upgrade to v0.2. Version 0.21 applies NHD HR etching to DEMs for all HUC6 units before running the HAND workflow. This new version also fixed the issue of void areas along HUC boundary. This dataset comprises 331 HUC6 units for CONUS (excluding the five great lakes units), each is a downloadable zip file. Version 0.1 was computed at the National Center for Supercomputing Applications (NCSA) at the University of Illinois at Urban-Champaign in 2016 and is currently hosted at the Texas Advanced Computing Center (TACC). Please cite this data DOI and the following publications and data DOIs when you use HAND and the hydraulic property table: Liu, Yan Y. Height Above Nearest Drainage (HAND) and Hydraulic Property Table for CONUS - Version 0.2. (20200301). Oak Ridge Leadership Computing Facility. (2020). DOI: 10.13139/ORNLNCCS/1608331 Liu, Yan Y., David R. Maidment, David G. Tarboton, Xing Zheng, and Shaowen Wang. 'A CyberGIS integration and computation framework for high-resolution continental-scale flood inundation mapping.' JAWRA Journal of the American Water Resources Association 54, no. 4 (2018): 770-784. DOI: 10.1111/1752-1688.12660 Zheng, Xing, David G. Tarboton, David R. Maidment, Yan Y. Liu, and Paola Passalacqua. 'River channel geometry and rating curve estimation using height above the nearest drainage.' JAWRA Journal of the American Water Resources Association 54, no. 4 (2018): 785-806. DOI: 10.1111/1752-1688.12661 Liu, Y. (2018). Height Above Nearest Drainage (HAND) for CONUS, HydroShare, https://doi.org/10.4211/hs.69f7d237675c4c73938481904358c789 LICENSE FOR USE -- MAPS AND DATA DISCLAIMER This resource is shared under the Creative Commons Attribution CC BY, http://creativecommons.org/licenses/by/4.0/ MAPS AND DATA DISCLAIMER The Oak Ridge National Laboratory (ORNL) shall not be held liable for improper or incorrect use of the data described or information contained on this map or associated series of maps. The data and related map graphics are not legal, land survey or engineering documents and are not intended to be used as such. ORNL gives no warranty, express or implied, as to the accuracy, reliability, utility or completeness of this information. The user of these maps and data assumes all responsibility and risk for the use of the maps and data. ORNL disclaims all warranties, representations or endorsements either express or implied, with regard to the information contained in this map product, including, but not limited to, all implied warranties of merchantability, fitness for a particular purpose or non-infringement. This preliminary map product is for research and review purposes only. It is not intended to be used for emergency management operational or life safety decisions at the local or regional governmental level or by the general public. Users requiring information regarding hazardous conditions or meteorological conditions for specific geographic areas should consult directly with their city or county emergency management office.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Weathering of chlorite from grain to watershed: The role and distribution of oxidation reactions in the subsurface

The reaction mechanism of weathering of chlorite, an important rock-forming phyllosilicate, is not well understood in natural settings. Here, in this work we investigated the weathering of Fe-rich chlorite from deep protolith to saprock to soil across a small shale-underlain watershed in the Appalachian Mountains, USA (Shale Hills). We found that oxidation of Fe(II) in chlorite always occurs prior to dissolution of the interlayers of the mineral. The oxidation of pyrite and chlorite commence near the water table across narrow depth intervals under the upper-catchment ridges, but well below the water table across wide depth intervals under the valley. We hypothesize that these patterns can be explained by hydrological and geochemical differences between the ridge and the valley: oxygenated water descends sub-vertically (1D flow) under the ridge, while under the valley, oxygen-depleted water moves upward to the stream and laterally out of the watershed in the subsurface (3D flow). Geochemical and mineralogical characterization indicates that the transformation of Fe-rich chlorite at Shale Hills is initiated by the oxidation of Fe(II). Next, the interlayer hydroxide sheet dissolves to form hydroxy-interlayered vermiculite and then vermiculite. During the transformation, Mg and Fe are released into solution and Fe is reprecipitated as goethite in pore space. Delivery of oxygen to the deep subsurface by infiltration of meteoric water is thought to control the initial transformation of chlorite at Shale Hills. It is possible that weathering of many Fe(II)-rich minerals is initiated by oxidation as mediated by rates of subsurface oxygen delivery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Groundwater elevation data for monitoring wells within the East and Taylor River basins, Colorado (USA)

This dataset is comprised of temporal variations in groundwater elevation data for the 24 monitoring wells located throughout the East River watershed. Seasonal to annual variations in groundwater elevations are a critical property of mountainous watersheds needed to understand both hydrological and below ground biogeochemical processes. Such data serve as a critical constraint for numerical models describing coupled groundwater-surface water behavior within the watershed. Additionally, the offset between the maximum and minimum groundwater elevations defines the extent of the bedrock weathering zone, with annual excursions in the groundwater hydrographic (i.e., the rising and falling hydrographic limbs) imposing primary controls on bedrock saturation state and redox conditions that govern biogeochemical reactions impacting nitrogen, carbon, and metals cycling. Manufacturer-specific software is used to download pressure data from each transducer, with broadly available spreadsheet software (e.g. Microsoft Excel) used to convert temporal variations in water pressure to elevations in units of meters above mean sea level. As additional monitoring wells are installed within the East River watershed and new groundwater monitoring wells are installed in the Taylor River watershed, temporal groundwater elevation data will be included as a part of this master dataset. Details regarding the metadata associated with each monitoring well location, including well depths, screened intervals, well location coordinates, and bedrock type, are included, as is a standard operating procedure for generating groundwater elevation data from water pressure values recorded by the pressure transducers. This dataset includes: (1) a zip file (East_River_Watershed_Compiled_Groundwater_Elevation_Data_Plots.zip), containing (a) PNG of groundwater hydrographs, (b) a CSV file with groundwater elevation data, and (c) CSV file containing metadata organized by location; (2) an Excel file (East_River_Watershed_Compiled_Groundwater_Elevation_Data_Plots.xlsx) with the groundwater elevation data, groundwater hydrographs, and metadata organized by location; (3) a Word file (Groundwater_elevation_data_protocols.docx) and a PDF file version (Groundwater_elevation_data_protocols.pdf) containing field protocols and methods; (4) a location metadata (locations.csv) file; (5) a file level metadata (flmd.csv); and (6) data dictionary (dd.csv) file.This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Utah FORGE: Phase 3 Ground Water Level Monitoring Data

Groundwater data around the Utah FORGE site has been collected to determine the piezometric levels and compositional variability. Field measurements in 2020 are designed to obtain and survey new groundwater wells that have been drilled on the distal periphery of the Utah FORGE project area. The spreadsheet included in this submission contains ground water level data from the Utah FORGE site as monitored from wells WOW-2 and WOW-3. The well locations are also provided as UTM, Zone 12, NAD83 coordinates.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Utah FORGE Groundwater Levels: Updated 2021

Contains Utah FORGE ground water data for wells WOW2 and WOW3. The data was updated on October 5th, 2021 and contains legacy data. Groundwater data includes the level, offset, date, and time for each measurement. Temperature, drift, water elevation and other parameters are recorded. Figures in the data include a water elevation over time plot. Legacy data ranges from year 1976 to 2019 and updated data ranges from year 2019 to 2021.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Model Scripts for "Old-Aged Groundwater Contributes to Mountain Hillslope Hydrologic Dynamics"

The partitioning of water inputs between deep and shallow groundwater flow paths is a fundamental processes, yet is challenging to observe. Numerical models provide a valuable tool to further develop insights on these groundwater mixing processes. This repository contains scripts to run the ParFlow-CLM and EcoSLIM integrated hydrologic models along the Pumphouse Hillslope in the East River Watershed and associated python scripts to process model outputs. Files includes input decks for the associated models, python scripts, and text and css files to support the model runs and processing. The 2-D hillslope model simulates from 2000 to 2021 using transient forcing conditions from the publicly available NLDAS-2 dataset. We develop an ensemble of models with variable hydrogeologic and soil parameters. Model outputs from the ensemble of runs are compared to to water level observations at the PLM1 and PLM6 wells published at (https://data.ess-dive.lbl.gov/view/doi:10.15485/1866836) and groundwater mean ages published at (https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1960042. The model results are used to evaluate hillslope water mass-balance transience and mixing dynamics between groundwater flow paths with young (<10 year) and old (>10 year) ages.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Data From: "Warming and snow loss increase reliance on old groundwater in a Colorado River headwater"

This repository contains the data and code associated with the paper titled "Warming and snow loss increase reliance on old groundwater in a Colorado River headwater," published in Nature Geoscience, 2026. This study seeks to answer how various ages of groundwater interact with mountainous streamflow in mountainous headwaters such as the East River. It includes various model-data processing scripts, primarily for ParFlow-CLM analysis of simulated water years 2015-2021, and two numerical warming experiments (+2.5 and +4.0 degrees C), including run scripts, forcing scripts, and post-processing, as well as comparison to observation datasets, detailed below. This data requires the use of R (.r, .rmd), Python (.py), Jupyter Notebook or Jupyter Lab (.ipynb), ParFLOW-CLM, EcoSLIM. Further information on the use of all file formats mentioned below (e.g. .tff. .nc) are provided within the associated scripts and directory where the files are located. Contents & Usage ASO/: ​​Contains the bash and python scripts used to convert airborne snow observatory (ASO) data (ASO, 2023) in various data formats (georeferenced tiff file, NetCDF, UTM, and to latitude/longitude) then regrided to the ParFlow equivalent grid. Output data are in regrid_regll_data.zip and subsequently visualized and analyzed in plot_and_compare.py for Supplementary Figures A14 and A15. The wksht_ASO_comparison.xlsx spreadsheet is used to calculate the data for Supplementary Figure A16. EcoSLIM/: Contains the scripts and input files to run the EcoSLIM particle tracking simulations (/run_scripts) and the post-processing python script (/plot_scripts/eco_agedist_plots.ipynb). Jasechko et al./: Contains the jupyter notebook (Extract_Elevation.ipynb) to determine the outlet elevations of the 260 watersheds used in Jasechko et al. (2016), and the corresponding table, Table_S1_Watersheds_alt.csv. Used to create Supplementary Information Figure A2. PLM_Wells/: Contains the QA/QC-ed groundwater level time series of the PLM-1 and PLM-6 Monitoring Wells from Faybishenko et al. (2023), reformatted to water years used for Supplementary Figures A19 and and A20. ParFlow/: Contains the input files and run scripts to run ParFlow-CLM (/run_scripts), the python and tool command language (Tcl) scripts to create and distribute the ParFlow forcing simulation files (/forcing), and various scripts and intermediary files to analyze the model outputs (/post_process). SQUIRE/: Contains the processing scripts and intermediary files for the Surface QUantitatIve pRecipitation Estimation (SQUIRE) data (Grover, 2023) used to generate Supplementary Figure A18. USGS_Streamflow/: Contains the raw and gap-filled United States Geological Survey streamflow data (U.S. Geological Survey, 2026) used at the Almont station (site number 09112500). Gap-filling is performed in the R script with data from the Taylor station (site number 09110000). (/USGS_09112500_EAST_RIVER_AT_ALMONT_GAP_FILLED/code_almont_streamflow_gap_fill.Rmd). discharge/: Contains the gap-filled discharge data at the Watershed Function SFA East River pumphouse site (Newcomer et al., 2022) used to generate Supplementary Figure A13 and to compute hourly Nash-Sutcliffe model efficiency coefficients (NSE) in Table A4. snotel_and_flux_tower/: Contains the snow telemetry data (U.S. Department of Agriculture, 2024) from the Butte (site ID 380) and Schofield (site ID 737) stations, reformatted by water year, accessed with the snotelr R package. Used to create Supplementary Figure A17. Also contains the flux tower observational data (FluxTower_Pumphouse_ESS-DIVE.ET_only.h.txt) from Ryken et al. (2022) and sap flux transpiration data (MaxB_Transpiration_5Sites.daily_sums.h.txt) from Ryken (2021), used to create Supplementary Figures A22 and A23, respectively. Raw EcoSLIM model outputs are in excess of 24TB, and are stored on National Energy Research Scientific Computing Center (NERSC) and publicly available via the external link provided in the paper.

atmospheric warming↗
DOE OSTI2025

Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1

Abstract. Groundwater plays a key role in meeting water demands, supplying over 40 % of irrigation water globally, with this role likely to grow as water demands and surface water variability increase. A better understanding of the future role of groundwater in meeting sectoral demands requires an integrated hydro-economic evaluation of its cost and availability. Yet substantial gaps remain in our knowledge and modeling capabilities related to groundwater availability, recharge, feasible locations for extraction, extractable volumes, and associated extraction costs, which are essential for large-scale analyses of integrated human–water system scenarios, particularly at the global scale. To address these needs, we developed Superwell, a physics-based groundwater extraction and cost accounting model that operates at sub-annual temporal and at the coarsest 0.5° (≈50 km × 50 km) gridded spatial resolution with global coverage. The model produces location-specific groundwater supply–cost curves that provide the levelized cost to access different quantities of available groundwater. The inputs to Superwell include recent high-resolution hydrogeologic datasets of permeability, porosity, aquifer thickness, depth to water table, recharge, and hydrogeological complexity zones. It also accounts for well capital and maintenance costs, as well as the energy costs required to lift water to the surface. The model employs a Theis-based scheme coupled with an amortization-based cost accounting formulation to simulate groundwater extraction and quantify the cost of groundwater pumping. The result is a spatiotemporally flexible, physically realistic, economics-based model that produces groundwater supply–cost curves. We show examples of these supply–cost curves and the insights that can be derived from them across a set of scenarios designed to explore model outcomes. The supply–cost curves produced by the model show that most (90 %) nonrenewable groundwater in storage globally is extractable at costs lower than USD 0.57 m−3, while half of the volume remains extractable at under USD 0.108 m−3. The global unit cost is estimated to range from a minimum of USD 0.004 m−3 to a maximum of USD 3.971 m−3. We also demonstrate and discuss examples of how these cost curves could be used by linking Superwell's outputs with other models to explore coupled human–environmental system challenges, such as water resources planning and management, or broader analyses of multisectoral feedbacks.

Global Change Analysis Model (GCAM)↗
DOE OSTI2021

Convergent Hydraulic Redistribution and Groundwater Access Supported Facilitative Dependency Between Trees and Grasses in a Semi‐Arid Environment

Abstract Hydraulic redistribution is the transport of water from wet to dry soil layers, upward or downward, through plant roots. Often in savanna and woodland ecosystems, deep‐rooted trees, and shallow‐rooted grasses coexist. The degree to which these different species compete for or share soil‐water derived from precipitation or groundwater, as well as how these interactions are altered by hydraulic redistribution, is unknown. We use a multilayer canopy model and field observations to examine how the presence of deep, but tree‐root accessible, groundwater impacts seasonal patterns of hydraulic redistribution, and interaction between coexisting vegetation species in a semiarid riparian woodland (US‐CMW). Based on the simulation, trees absorb moisture at the water table (∼10 m depth) and release it in the shallow soil depth (0–3 m) during the dry pre‐monsoon season. We observed the occurrence of a new convergent hydraulic redistribution pattern during the monsoon season, where moisture is transported from both the near‐surface (0–0.5 m) and the water table to intermediate soil layers (1–5 m) through tree roots. We found that hydraulic redistribution demonstrates a growth facilitation effect at this site, supporting 49% of growing season tree transpiration and 14% of the grass transpiration. Compared to a similarly structured upland savanna without accessible groundwater, the riparian site shows an increased amount of hydraulically redistributed water and more facilitative water use between coexisting grasses and trees. These results shed light on the linkage between accessible groundwater and the role of hydraulic redistribution on the interaction between deep‐rooted and shallow‐rooted vegetation.

Lee, E.↗
DOE OSTI2021

Influence of Agricultural Managed Aquifer Recharge (AgMAR) and Stratigraphic Heterogeneities on Nitrate Reduction in the Deep Subsurface

Agricultural managed aquifer recharge (AgMAR) is a strategy whereby surface water is used to intentionally flood croplands to recharge underlying aquifers. However, nitrate (NO 3 - ) contamination in agriculturally intensive regions poses a threat to groundwater resources under AgMAR. We use a reactive transport model to understand the effects of AgMAR management strategies (i.e., by varying the frequency, duration between flooding events, and amount of water) on NO 3 - leaching to groundwater under different stratigraphic configurations and antecedent moisture conditions. We examine the potential of denitrification and nitrogen retention in deep vadose zone sediments (~15 m) using variable AgMAR application rates on two-dimensional representations of differently textured soils, soils with discontinuous bands/channels, and with preferential flow paths characteristic of agricultural fields. Simulations indicate finer textured sediments, alone or embedded within/adjacent to high flow regions, are important reducing zones providing conditions needed for denitrification. Simulation results suggest that applying water all-at-once rather than in increments transports higher concentrations of NO 3 - deeper into the profile, which may exacerbate groundwater quality. This transport into deeper depths can be aggravated by wetter antecedent soil moisture conditions. However, applying water all-at-once also increases denitrification within the vadose zone by promoting anoxic conditions. We conclude that AgMAR can be designed to enhance denitrification in the subsurface and reduce NO 3 - leaching to groundwater, while specifically accounting for lithologic heterogeneity, antecedent soil moisture conditions, and depth to the water table. Our findings are potentially relevant to other systems that experience flooding inundation such as floodplains and dedicated recharge basins.

54 ENVIRONMENTAL SCIENCES↗