Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “WATER TABLES”

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

Volumetric Soil Moisture Measurements at the Teller 27 Site, Seward Peninsula, Alaska, 2022-2023

The Teller 27 watershed on the Seward Peninsula, Alaska, has been well characterized by the NGEE Arctic project. The study site is underlain by discontinuous permafrost that is thawing as the climate warms. As a result, the site is experiencing a short-term wetting trend as a perched water table above the remaining permafrost provides plant available water during the growing season. Soil moisture patterns drive microbial activity and plant species compositions including plant density and height. This study aimed to understand how soil moisture patterns were influenced by tundra microtopography. To accomplish this, we placed a strategic network of soil moisture sensors in micro-highs, micro-lows, and control areas under different vegetation types within the Teller 27 watershed from summer of 2022 through fall of 2023. This data was used in conjunction with other soil data from the Teller 27 watershed to gain a more comprehensive understanding of soil moisture patterns in a rapidly thawing discontinuous permafrost region. This dataset includes six *.csv files: four of time series soil moisture data, one of field soil moisture data, and one of site conditions. The dataset also includes one *.kml file of the watershed and the soil moisture sensor sites as well as this user guide to provide details on data collection and processing methods. NGEE Arctic Project Summary The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy’s Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy’s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
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 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 OSTI2025

Soil oxygen dynamics: a key mediator of tile drainage impacts on coupled hydrological, biogeochemical, and crop systems

Tile drainage removes excess water and is an essential, widely adopted management practice to enhance crop productivity in the US Midwest and throughout the world. Tile drainage has been shown to significantly change hydrological and biogeochemical cycles by lowering the water table and reducing the residence time of soil water, although examining the complex interactions and feedbacks in an integrated hydrology–biogeochemistry–crop system remains elusive. Oxygen dynamics are critical to unraveling these interactions and have been ignored or oversimplified in existing models. Understanding these impacts is essential, particularly so because tile drainage has been highlighted as an adaptation under projected wetter springs and drier summers in the changing climate in the US Midwest. We used the ecosys model that uniquely incorporates first-principle soil oxygen dynamics and crop oxygen uptake mechanisms to quantify the impacts of tile drainage on hydrological and biogeochemical cycles and crop growth in corn–soybean rotation fields. The model was validated with data from a multi-treatment, multi-year experiment in Washington, IA. The relative root mean square error (rRMSE) for the corn and soybean yield in validation is 5.66 % and 12.57 %, respectively. The Pearson coefficient (r) of the monthly tile flow during the growing season is 0.78. Plant oxygen stress turns out as an emergent property of the equilibrium between the soil oxygen supply and biological demand. The impact of tile drainage on the system is achieved through a series of coupled feedback mechanisms. The model results show that tile drainage reduces the soil water content and enhances soil oxygenation. It additionally increases the subsurface discharge and elevates inorganic nitrogen leaching, with seasonal variations influenced by climate and crop phenology. The improved aerobic condition alleviates crop oxygen stress during wet springs, thereby promoting crop root growth during the early growth stage. The development of greater root density, in turn, mitigates water stress during dry summers, leading to an overall increase in the crop yield by ∼6 %. These functions indicate the potential of tile drainage in bolstering crop resilience to climate change and the use of this modeling tool for large-scale assessments of tile drainage. The model reveals the underlying causal mechanisms that drive the agroecosystem response to drainage on the coupled hydrology, biogeochemistry, and crop system dynamics.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Geophysical Characterization and Monitoring of the 200 Area Treated Effluent Disposal Facility to Support Permit Renewal; TEDF Geophysical Monitoring

The Treated Effluent Disposal Facility (TEDF), located in the 200 East Area of the Hanford Site, is a site where non-hazardous and non-radioactive liquid wastes are disposed into two state-permitted infiltration basins. In 2016, the Washington State Department of Ecology denied a permit renewal request for TEDF due to the inability to adequately assess the impact of TEDF discharge water on the underlying groundwater quality. TEDF overlies the relatively impermeable Ringold Lower Mud (RLM) unit, whose upper contact lies in the vadose zone approximately 30 m below ground surface and approximately 10 m above the water table boundary. The RLM is assumed to isolate TEDF discharge water from the natural groundwater aquifer, which is monitored using wellbores screened below the RLM. Therefore, samples collected from monitoring wells near TEDF are not considered representative of TEDF discharge water. Rather, TEDF discharge water is assumed to mound on top of the RLM to form a perched aquifer. To support permitting of the TEDF, a new monitoring well is required that can be used to sample water from the presumed perched water aquifer above the RLM. Ideally, the screened section of the well would be located at the peak of the perched water mound(s), which presumably occurs at the point of maximum vertical flux from TEDF to the RLM, or equivalently where the dominant infiltration flow paths reach the RLM. This report describes how time-lapse 3D electrical resistivity tomography (ERT) was used in conjunction with nominal TEDF discharge operations to image the dominant flow paths from each pond to the RLM. Results are summarized in Figure ES.1. Figure ES.1A shows a satellite image of the TEDF overlain by an array of surface ERT electrodes. The solid and dashed black circles denote the zones of maximum vertical flux at the RLM within the south and north ponds, respectively, and presumably the regions where perched water peaks during discharge. Figure ES.1B shows time-lapse difference imaging results approximately 15 days after switching discharge from the south pond to the north pond. Blue iso-surfaces beneath the north pond denote regions of increasing bulk electrical conductivity caused by increasing saturation due to infiltrating water. Red iso-surfaces beneath the south pond denote regions that were previously saturated during the south pond discharge and are now de-saturating, causing a decrease in bulk electrical conductivity. In both cases, the zones of maximum change mark the dominant flow paths to the RLM. Figure ES.1C shows time-lapse imaging results approximately 15 days after switching discharge from the north pond to the south pond. In this case, blue iso-surfaces mark the dominant flow paths to the RLM from the south pond. The red iso-surfaces mark the dominant flow paths that existed during discharge to the north pond. The regions of maximum change in bulk conductivity (due to saturation or desaturation) that mark the primary flow paths are equivalent in both cases and denoted by the dashed solid and black circles. If a perched water zone forms on the RLM, it is likely to mound within or near the dashed circle during discharge to the north pond, and in the solid circle during discharge to the south pond. In other words, if perched water mounding occurs, the circles mark the optimum locations suggested by the ERT imaging for monitoring boreholes to be placed, enabling samples to be collected that are representative of TEDF discharge impacts on groundwater quality.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Seismic refraction tracks porosity generation and possible CO 2 production at depth under a headwater catchment

Significance To understand groundwater in weathered bedrock aquifers, we explore (bio)geochemical reactions using geophysical measurements. The result leads to a conceptual model for the three-dimensional distribution of geochemical alteration in the subsurface—well beyond the typical approach of inferring (bio)geochemical reactions in the subsurface from pits, boreholes, and outcrops. With this geophysical technique we discovered the prevalence of gas beneath the water table at this site. This geophysical viewing of the subsurface highlights patterns of water movement, reactions of clays, and microbial and abiotic gas-producing reactions recorded in the subsurface and thus illuminates storage, flow, and reactions of groundwater in bedrock systems. Protecting the groundwater resource beneath our feet will rely on such enhanced capacities to see into the subsurface.

58 GEOSCIENCES↗
DOE OSTI2023

Attaining freshwater and estuarine-water soil saturation in an ecosystem-scale coastal flooding experiment

Abstract Coastal upland forests are facing widespread mortality as sea-level rise accelerates and precipitation and storm regimes change. The loss of coastal forests has significant implications for the coastal carbon cycle; yet, predicting mortality likelihood is difficult due to our limited understanding of disturbance impacts on coastal forests. The manipulative, ecosystem-scale Terrestrial Ecosystem Manipulation to Probe the Effects of Storm Treatments (TEMPEST) experiment addresses the potential for freshwater and estuarine-water disturbance events to alter tree function, species composition, and ecosystem processes in a deciduous coastal forest in MD, USA. The experiment uses a large-unit (2000 m 2 ), un-replicated experimental design, with three 50 m × 40 m plots serving as control, freshwater, and estuarine-water treatments. Transient saturation (5 h) of the entire soil rooting zone (0–30 cm) across a 2000 m 2 coastal forest was attained by delivering 300 m 3 of water through a spatially distributed irrigation network at a rate just above the soil infiltration rate. Our water delivery approach also elevated the water table (typically ~ 2 m belowground) and achieved extensive, low-level inundation (~ 8 cm standing water). A TEMPEST simulation approximated a 15-cm rainfall event and based on historic records, was of comparable intensity to a 10-year storm for the area. This characterization was supported by showing that Hurricane Ida’s (~ 5 cm rainfall) hydrologic impacts were shorter (40% lower duration) and less expansive (80% less coverage) than those generated through experimental manipulation. Future work will apply TEMPEST treatments to evaluate coastal forest resilience to changing hydrologic disturbance regimes and identify conditions that initiate ecosystem state transitions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Vadose Zone Model for A Trenches 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. 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 parallel version of the Subsurface Transport Over Multiple Phases (STOMP) simulator, officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP), is used to simulate flow and transport for the vadose models. The scope of this ECF is to document the development and results of the A Trenches Area vadose zone model. The A Trenches Area model covers the region between the 200 East Area and southwest of the 216-B-3 Pond. The primary waste sites in this area are cribs associated with PUREX Plant (202-A).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Vadose Zone Model for B Complex 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. 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 simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The scope of this ECF is to document the development and results of the B Complex vadose zone model. The B Complex model encompasses the region around single-shell tank farms 241-B, 241-BX, and 241-BY, and adjacent cribs, trenches, and burial grounds in the northwestern part of the 200 East Area. Sources of waste to this area were primarily from B Plant (221-B) and U Plant (221-U).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Vadose Zone Model for B Plant 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. 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 simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The scope of this ECF is to document the development and results of the B Plant Area vadose zone model. The B Plant Area model encompasses the region around B Plant (221-B) in the western part of 200 East Area. The plant was used to separate plutonium from irradiated fuel using the bismuth-phosphate process from 1945 to 1956, and for strontium and cesium recovery from 1968 to 1985.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Vadose Zone Model for B-3 Pond 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. 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 simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The scope of this ECF is to document the development and results of the B-3 Pond vadose zone model. The B-3 Pond model covers the main lobe of the 216-B-3 Pond System located in a topographic low east of the 200 East Area. The pond system operated from 1945 to 1997 and received effluents from the Plutonium Uranium Extraction Plant, B Plant, 284-E Powerhouse, and other 200 East Area facilities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Vadose Zone Model for B-3C Pond 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. 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 simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The scope of this ECF is to document the development and results of the B-3C Pond vadose zone model. The B-3C Pond model covers the 216-B-3C Expansion Pond of the 216-B-3 Pond System located in a topographic low east of the 200 East Area. The pond system operated from 1945 to 1997 and received effluents from the Plutonium Uranium Extraction Plant, B Plant, 284-E Powerhouse, and other 200 East Area facilities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

BAWLD-CH 4 : a comprehensive dataset of methane fluxes from boreal and arctic ecosystems

Methane (CH 4 ) emissions from the boreal and arctic region are globally significant and highly sensitive to climate change. There is currently a wide range in estimates of high-latitude annual CH 4 fluxes, where estimates based on land cover inventories and empirical CH 4 flux data or process models (bottom-up approaches) generally are greater than atmospheric inversions (top-down approaches). A limitation of bottom-up approaches has been the lack of harmonization between inventories of site-level CH 4 flux data and the land cover classes present in high-latitude spatial datasets. Here we present a comprehensive dataset of small-scale, surface CH 4 flux data from 540 terrestrial sites (wetland and non-wetland) and 1247 aquatic sites (lakes and ponds), compiled from 189 studies. The Boreal–Arctic Wetland and Lake Methane Dataset (BAWLD-CH 4 ) was constructed in parallel with a compatible land cover dataset, sharing the same land cover classes to enable refined bottom-up assessments. BAWLD-CH 4 includes information on site-level CH 4 fluxes but also on study design (measurement method, timing, and frequency) and site characteristics (vegetation, climate, hydrology, soil, and sediment types, permafrost conditions, lake size and depth, and our determination of land cover class). The different land cover classes had distinct CH 4 fluxes, resulting from definitions that were either based on or co-varied with key environmental controls. Fluxes of CH 4 from terrestrial ecosystems were primarily influenced by water table position, soil temperature, and vegetation composition, while CH 4 fluxes from aquatic ecosystems were primarily influenced by water temperature, lake size, and lake genesis. Models could explain more of the between-site variability in CH 4 fluxes for terrestrial than aquatic ecosystems, likely due to both less precise assessments of lake CH 4 fluxes and fewer consistently reported lake site characteristics. Analysis of BAWLD-CH 4 identified both land cover classes and regions within the boreal and arctic domain, where future studies should be focused, alongside methodological approaches. Overall, BAWLD-CH 4 provides a comprehensive dataset of CH 4 emissions from high-latitude ecosystems that are useful for identifying research opportunities, for comparison against new field data, and model parameterization or validation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Continental Scale Hydrostratigraphy: Basin‐Scale Testing of Alternative Data‐Driven Approaches

Abstract Integrated hydrological modeling is an effective method for understanding interactions between parts of the hydrologic cycle, quantifying water resources, and furthering knowledge of hydrologic processes. However, these models are dependent on robust and accurate datasets that physically represent spatial characteristics as model inputs. This study evaluates multiple data‐driven approaches for estimating hydraulic conductivity and subsurface properties at the continental‐scale, constructed from existing subsurface dataset components. Each subsurface configuration represents upper (unconfined) hydrogeology, lower (confined) hydrogeology, and the presence of a vertical flow barrier. Configurations are tested in two large‐scale U.S. watersheds using an integrated model. Model results are compared to observed streamflow and steady state water table depth (WTD). We provide model results for a range of configurations and show that both WTD and surface water partitioning are important indicators of performance. We also show that geology data source, total subsurface depth, anisotropy, and inclusion of a vertical flow barrier are the most important considerations for subsurface configurations. While a range of configurations proved viable, we provide a recommended Selected National Configuration 1 km resolution subsurface dataset for use in distributed large‐and continental‐scale hydrologic modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Biogeochemical controls on iron speciation and cycling across upland to shoreline gradients in freshwater and estuarine coastal soils (Lake Erie and Chesapeake Bay, United States)

Coastal environments are dynamic interfaces that mediate carbon and nutrient exchanges between terrestrial landscapes and open waters, and understanding the biogeochemical factors controlling these exchanges, particularly iron (Fe) redox transformations, is crucial for predicting coastal ecosystem functions. Here, we investigated the mechanisms controlling Fe speciation changes across upland-to-shoreline gradients in freshwater and estuarine soils using Fe K-edge X-ray absorption spectroscopy, solid and porewater composition analysis, and 16S rRNA sequencing analysis. We show that Fe transformations depend primarily on inundation patterns. In unsaturated uplands, Fe occurs as Fe(III) oxyhydroxides, mainly goethite (9–35 %), Fe(II,III)-phyllosilicates (39–89 %), and Fe(III)-organic species (0–61 %). Soils influenced by estuarine waters exhibit porewater sulfide concentrations reaching up to 221 μM, Fe- and S-cycling bacteria, and up to 81 % pyrite (FeS 2 ), indicating that sulfur-driven redox dynamics control Fe transformations. In lacustrine wetlands, Fe(III) reduction is indicated by porewater Fe(II) concentrations increasing to 1.0–2.1 mM, and ~10–15 % of Fe as Fe(II,III)-(hydr)oxides (green rust), vivianite (Fe 3 (PO 4 ) 2 ·8H 2 O), and/or adsorbed Fe(II) species. EXAFS data also indicate reduction of structural Fe(III) to Fe(II) in phyllosilicates. The presence of Fe- and S-cycling bacteria, as well as sulfide (0–10 μM), suggests that Fe-cycling is microbially driven and potentially coupled with cryptic S-cycling. Fe(II) oxidation was indicated above/near the water table by the presence of Fe(III) oxyhydroxides (ferrihydrite, lepidocrocite). Furthermore, negligible Fe(III) or sulfate reduction was observed at some water-saturated sites located at the upland-wetland transition, likely due to oxic (sub-)surface water inputs. Overall, our results highlight the importance of considering both Fe-speciation and hydro-biogeochemical dynamics when predicting Fe-cycling at coastal interfaces.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Response of Subsurface Nitrogen-Cycling Microbial Communities to Environmental Fluctuations (Final Technical Report)

Riparian floodplains are dynamic ecosystems linking terrestrial and riverine systems. These floodplains experience hydrological shifts such as changes in water table height, flooding, and drought and can be ‘hotspots’ of biogeochemical cycling due to shifting sediment moisture (and saturation) and subsurface exchanges of water, nutrients, and other compounds across different sediment layers. Subsurface microbial communities are the primary drivers of biogeochemical processes in floodplains, and thus their structure and function can directly influence both surface and groundwater quality. The microbial nitrogen (N) cycle is particularly important in floodplains as it affects nutrient availability and removal. Two functional guilds of chemoautotrophic (i.e. CO2-fixing) microorganisms are responsible for the first oxidative step of the N cycle, nitrification: ammonia-oxidizing archaea (AOA) and bacteria (AOB) catalyze the oxidation of ammonia to nitrite, while nitrite-oxidizing bacteria (NOB) oxidize nitrite to nitrate. Despite the critical role nitrification plays in N-cycling in both terrestrial and aquatic ecosystems, our understanding of the diversity, ecophysiology, and activity of nitrifying organisms in subsurface floodplain soils/sediments is extremely limited. To help address this critical knowledge gap, the overarching goal of this project was to determine how shifts in key environmental parameters and gradients impact microbial N-cycling communities/processes, with particular emphasis on nitrification, within hydrologically-variable floodplain sediments in the Wind River Basin near Riverton, Wyoming. The three specific objectives of this project were to: (1) to associate in situ environmental drivers of N cycling with distinct functional guilds; (2) determine the guild response to variation in key ecosystem drivers; and (3) develop a dynamic ecosystem model of the microbial N cycle with the Riverton subsurface using community genomic and biogeochemical data collected in the first two objectives. Over the course of this project, we employed both 16S rRNA gene amplicon sequencing and genome-resolved metagenomics to examine the phylogenetic diversity and metabolic potential of subsurface nitrifier communities within 68 samples collected across multiple sites, depths, and time points within the Riverton floodplain, allowing for both spatial and temporal investigations at different scales. This project benefitted tremendously from recent advances in high-throughput sequencing technologies coupled with dramatic improvements in the computational tools and algorithms available for analyzing such large, complex genomic datasets. By pairing these cutting-edge genomic approaches with depth-resolved sampling and detailed geochemical analyses of the Riverton floodplain, we have gained novel insights into the structure and function of subsurface nitrifier communities in relation to both hydrology and biogeochemistry. This project resulted in the most detailed and comprehensive characterization of N-cycling floodplain microbial communities to date and will hopefully inspire and pave the way for future studies using similar approaches in other floodplains. Indeed, such information is critical for understanding subsurface biogeochemical cycling and how elemental stores are altered from perturbations initiated by the water cycle within floodplains. Finally, because of the terrestrial-aquatic nature of the Riverton floodplain, results from this project are also of relevance to disciplines such as soil science, estuarine science, limnology & oceanography, biogeochemistry, geobiology, environmental engineering, as well as genomics and data science.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Spatio–Temporal Machine Learning for Regional to Continental Scale Terrestrial Hydrology

Integrated hydrologic models can simulate coupled surface and subsurface processes but are computationally expensive to run at high resolutions over large domains. Here we develop a novel deep learning model to emulate subsurface flows simulated by the integrated ParFlow–CLM model across the contiguous US. We compare convolutional neural networks like ResNet and UNet run autoregressively against our novel architecture called the Forced SpatioTemporal RNN (FSTR). The FSTR model incorporates separate encoding of initial conditions, static parameters, and meteorological forcings, which are fused in a recurrent loop to produce spatiotemporal predictions of groundwater. We evaluate the model architectures on their ability to reproduce 4D pressure heads, water table depths, and surface soil moisture over the contiguous US at 1 km resolution and daily time steps over the course of a full water year. The FSTR model shows superior performance to the baseline models, producing stable simulations that capture both seasonal and event–scale dynamics across a wide array of hydroclimatic regimes. The emulators provide over 1,000× speedup compared to the original physical model, which will enable new capabilities like uncertainty quantification and data assimilation for integrated hydrologic modeling that were not previously possible. Our results demonstrate the promise of using specialized deep learning architectures like FSTR for emulating complex process–based models without sacrificing fidelity.

54 ENVIRONMENTAL SCIENCES↗
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)↗