Search NASA⌕ Search

SEARCH · Search NASA

Results for “Groundwater Modeling”

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

Resilient water infrastructure partnerships in institutionally complex systems face challenging supply and financial risk tradeoffs

Abstract As regions around the world invest billions in new infrastructure to overcome increasing water scarcity, better guidance is needed to facilitate cooperative planning and investment in institutionally complex and interconnected water supply systems. This work combines detailed water resource system ensemble modeling with multiobjective intelligent search to explore infrastructure investment partnership design in the context of ongoing canal rehabilitation and groundwater banking in California. Here we demonstrate that severe tradeoffs can emerge between conflicting goals related to water supply deliveries, partnership size, and the underlying financial risks associated with cooperative infrastructure investments. We show how hydroclimatic variability and institutional complexity can create significant uncertainty in realized water supply benefits and heterogeneity in partners’ financial risks that threaten infrastructure investment partnership viability. We demonstrate how multiobjective intelligent search can design partnerships with substantially higher water supply benefits and a fraction of the financial risk compared to status quo planning processes. This work has important implications globally for efforts to use cooperative infrastructure investments to enhance the climate resilience and financial stability of water supply systems.

Science & Technology - Other Topics↗
DOE OSTI2026

Superstructure Optimization for Brine Valorization from Brackish Water Desalination

This poster presents preliminary results from a superstructure optimization framework developed to identify cost-optimal brine valorization configurations for brackish water desalination plants across diverse U.S. regional feed chemistries. The study uses brackish groundwater compositions from Arizona, California, Florida, New Mexico, and Texas. Using Pyomo Generalized Disjunctive Programming (GDP) within the WaterTAP modeling environment, the optimization framework simultaneously evaluates thousands of candidate treatment configurations, spanning nanofiltration, reverse osmosis, and chemical precipitation, to minimize the levelized cost of water (LCOW) while meeting water recovery targets and product recovery constraints. Results across eight representative feed clusters demonstrate water recovery rates of 57–87% and net LCOW values ranging from -$0.032/m³ (net revenue-positive) to $0.80/m. Notably, no single process configuration was optimal across all feed types, underscoring the necessity of feed-specific optimization. Products targeted include calcium carbonate (CaCO₃) at $0.01/kg and sodium chloride (NaCl) at $0.10/kg, both at 95% purity, with product revenues offsetting treatment costs in several scenarios. The work advances NAWI's process systems engineering capabilities for multi-configuration screening.

58 GEOSCIENCES↗
DOE OSTI2025

Pantex Plant Ogallala Aquifer and Perched Groundwater Contingency Plan

The Pantex Plant Ogallala Aquifer and Perched Groundwater Contingency Plan has been developed in accordance with the requirements identified in the: • Interagency Agreement for the Pantex Superfund Site, Article 8.5 Work to be Performed, • Compliance Plan Provision of Hazardous Waste Permit No. 50284, and • Record of Decision for Groundwater, Soil, and Associated Media, Pantex Plant. A Long‐Term Monitoring System Design has been designed to monitor conditions in the perched groundwater including changes in the perched aquifer as a result of implementing the response actions. Monitoring is required for verifying the effectiveness of perched groundwater response actions (i.e., conditions in the perched aquifer are being affected as intended) and for confirming that the perched aquifer and Ogallala Aquifer characterization as defined in the Resource Conservation and Recovery Act Facility Investigation Report and the Corrective Measure Studies/Feasibility Study remains accurate. If monitoring results obtained through the monitoring network identify an unexpected condition or deviation, contingent actions will be considered and implemented as necessary to ensure continued protection of the Ogallala Aquifer and human health and the environment. Potential deviations to expected technology performance may be encountered for each of the four primary response actions that compose the selected remedy for perched groundwater; Playa 1 Pump and Treat System, Southeast Area Pump and Treat System, Southeast Area In‐Situ Bioremediation System (comprised of the Southeast In‐Situ Bioremediation System Original System, Southeast Area In‐Situ Bioremediation System Extension System, Offsite In‐Situ Bioremediation System, Perchlorate/Chromium ISB, Northeast ISB and County Road 8 ISB), and Zone 11 In‐Situ Bioremediation System. Monitoring will also be conducted to determine if there are deviations to the expected characterization, e.g., contaminants not expected as a result of the RCRA Facility Investigation characterization. Deviations to expected conditions in the Ogallala Aquifer could also be encountered if the response actions in the perched groundwater are not performing as expected, i.e., preventing contaminants from migrating to the Ogallala Aquifer. Currently, Pantex has begun investigation of detections of high explosives above groundwater protection standards in wells on the Texas Tech University property and a plume that is moving to the northeast from that area. Due to those detections, this Plan recognizes the fact that future detections in the Ogallala will be focused on first‐ time detections of analytes. After a remedy is determined, this Plan will require modification to address-deviations and contingent actions. This Plan was developed to identify the contingent actions necessary to mitigate impacts resulting from deviations to site conditions or response action performance. The Plan defines the environmental problem being addressed by the response actions, clarifies the expected conditions and objectives of the response actions, and identifies the potential deviations to the response actions (due to site conditions or technology performance) that could be encountered. The deviations were evaluated to determine the likelihood of occurrence, potential impact, and time to respond to avoid impact. The Plan also identifies the monitoring outlined in the Long‐Term Monitoring System Design Report (Consolidated Nuclear Security, 2024) and Sampling Analysis Plan (PanTeXas Deterrence, 2024) that will be used to detect the deviations. Lastly, the Plan specifies the contingent actions that could be implemented in response to the deviations. Because each response focuses on a discrete portion of the perched aquifer and contaminant plume, each response action has a different set of expected conditions, and therefore differing impacts from deviations to the site and technology expectations. As a result, the contingent actions are identified for each response action and potential deviation including specific constituents, location, and conditions. If deviations are encountered that impact the ability of the response action to meet performance objectives, the contingent actions will be focused on ensuring the response action can meet the performance objective. Contingent actions may be implemented as interim actions (ISMs/removal actions) in accordance with the Record of Decision, Interagency Agreement, and Hazardous Waste Permit‐50284, if warranted by the specific circumstances. For deviations to site characterization expected conditions, the contingent action will focus on determination of the source of the deviation, determination of the appropriate response, and evaluation of additional work to be completed. However, if the deviation to characterization impacts the performance of the response action, the contingent action will again focus on ensuring performance objectives can be met. Early source term removals and cleanup actions have been implemented to protect the Ogallala Aquifer. Because of these actions and based on modeling results, the expected conditions in the Ogallala Aquifer are that constituents of concern will not be detected above the Groundwater Protection Standards (GWPSs) nor will they reach potential points of exposure above the GWPS. The primary deviation of concern for the Ogallala is if constituents are detected in the Ogallala Aquifer near or above GWPSs. If it occurs, this change in expected conditions would require further evaluation of site and contaminant characteristics to determine an appropriate course of action. The evaluation would include additional monitoring, source identification, implementation of interim protective measures (if necessary), and delineation of extent. These evaluations are necessary to determine an appropriate response action for the Ogallala. The primary goal of the Plan is to provide for the continued protection of the Ogallala Aquifer and the health of its consumers. In recognition, this Plan presents a flexible and rational approach for making future decisions associated with confirming the change in perched and Ogallala aquifer conditions and identifying a response (technical activities, changes to response actions, regulatory oversight, and public involvement).

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

Flow Directions and Ages of Subsurface Water in a Salt Marsh System Constrained by Isotope Tracing

Salt marshes are dynamic hydrologic systems where terrestrial groundwater, terrestrial surface water, and seawater mix due to bi-directional flows and pressure gradients. Due to the counteracting terrestrial and marine forcings that control these environments, we do not comprehensively understand water fluxes in these complex coastal systems. To understand the water sources, flow directions, and velocities in salt marsh porewater, we employed a combination of geochemical tracers and analytical models across a hillslope-to-salt marsh continuum in a salt marsh experiencing daily inundation of estuarine surface water (SW) from tides and mixing of fresh seasonal groundwater. We used tritium ( 3 H) as a hydrologic tracer to assess porewater ages and stable water isotope ($\delta$ 2 H and $\delta$ 18 O) analyses to separate isotopically distinct estuarine and terrestrial groundwater across different depths and landscape positions in the study transect. We employed electrical conductivity to constrain the role of source mixing and evapotranspiration in salt marsh hydrology. Salinity and stable isotopes revealed that transpiration, rather than evaporation, increased subsurface water salinity to concentrations above estuarine SW during summer. Elevated salinity at depth indicated that salt marsh subsurface water is recharged during the dry growing season. Seasonal recharge patterns drive long-term deep subsurface water dynamics across the salt marsh, with 3 H ages of 3–7 years, and daily tidal cycles drive short-term shallow porewater dynamics with 3 H ages of 0 ± 3.6 years. Furthermore, our conceptual understanding of the spatiotemporal changes in SW-subsurface water interactions at the terrestrial-marine interface quantifies the hydrological constraints we are missing to improve our understanding of biogeochemical cycles within the salt marsh.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Model Data Archive Associated with Manuscript "Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon"

This data package supports the publication “Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon” by Li et al. (2026). The package contains processed model inputs, configuration files, restart files, simulation outputs, scripts, and visualization products used to evaluate post-fire dissolved organic carbon (DOC) dynamics in the Naches River Watershed, Washington, USA, following the 2021 Schneider Springs Fire. The modeling workflow couples ELM-BGC, the biogeochemistry-enabled Energy Exascale Earth System Model Land Model; ATS, the Advanced Terrestrial Simulator for integrated surface-subsurface hydrology; and PFLOTRAN, a reactive transport model for multicomponent aqueous geochemistry. Together, these models simulate how wildfire-induced changes in vegetation, litter, coarse woody debris, and soil organic matter influence DOC production, transport, and reaction from burned hillslopes to stream networks. The archive includes preprocessed meteorological, geospatial, hydrologic, and biogeochemical forcing data; ELM-BGC-derived DOC source terms; ATS mesh files; PFLOTRAN reactive-transport inputs; model configuration files; spin-up and transient restart files; watershed-scale diagnostic outputs; stream concentration time series; and figures or visualization files used to inspect and reproduce key results. File types include Hierarchical Data Format 5 (HDF5) files for gridded forcing and model-coupling data, model input and configuration files for ELM-BGC, ATS, and PFLOTRAN, restart and simulation-output files generated by the modeling workflow, tabular or time-series diagnostic outputs, scripts for post-processing and figure generation, and image or visualization products associated with the manuscript. Use of the package depends on the intended task. Re-running the simulations requires the relevant modeling software, including ELM-BGC, ATS, and PFLOTRAN as ATS's geochemical engine. Inspecting outputs and reproducing figures requires Python with scientific plotting libraries such as Matplotlib, and three-dimensional model outputs may be viewed with ParaView. Geographic information system files or maps may be inspected with ArcGIS Pro or comparable GIS software. The data package is intended to enable traceability, reuse, and partial reproduction of the coupled land-to-watershed hydro-biogeochemical modeling workflow used to test how wildfire disturbance affects terrestrial carbon pools and downstream DOC dynamics.

ATS↗
NASA NTRS1988

Computer modeling of the mineralogy of the Martian surface, as modified by aqueous alteration

Mineralogical constraints can be placed on the Martian surface by assuming chemical equilibria among the surface rocks, atmosphere and hypothesized percolating groundwater. A study was made of possible Martian surface mineralogy, as modified by the action of aqueous alteration, using the EQ3/6 computer codes. These codes calculate gas fugacities, aqueous speciation, ionic strength, pH, Eh and concentration and degree of mineral saturation for complex aqueous systems. Thus, these codes are also able to consider mineralogical solid solutions. These codes are able to predict the likely alteration phases which will occur as the result of weathering on the Martian surface. Knowledge of the stability conditions of these phases will then assist in the definition of the specifications for the sample canister of the proposed Martian sample return mission. The model and its results are discussed.

Zolensky, M. E.↗
DOE OSTI2026

Inputs to GCAM-USA: IM3 Phase 2 Experiments

Overview This dataset contains XML input files for the IM3 Phase 2 version of GCAM-USA. The files are organized into two categories: Scenario-specific inputs represent hydroclimate and socioeconomic effects on water availability, heating and cooling degree-hours, and agricultural productivity. They support eight IM3 canonical scenarios: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 Model-improvement inputs extend GCAM-USA v5.3 with updated representations of coal and nuclear power plant retirements, electricity trade among U.S. interconnections, offshore carbon storage costs, and groundwater depletion constraints. Data structure Scenario-specific inputs rcp45_runoff/ and rcp85_runoff/XML files describing water availability by HUC2 basin under the RCP 4.5 and RCP 8.5 scenarios. rcp45_hdcd/ and rcp85_hdcd/XML files containing monthly-day and monthly-night heating and cooling degree-hours at the U.S. state level for different RCP-SSP combinations. rcp45_agyields/ and rcp85_agyields/XML files describing changes in agricultural productivity at the intersection of GCAM regions and HUC2 water basins for different RCP-SSP combinations. rcp45_emissions_pathway/The emissions-constraint XML file used to represent the RCP 4.5 pathway. Model-improvement inputs core_retire/Updates coal-fired power plant retirement schedules based on New England ISO. GCAMUSA_IM3_elec_trade_interconnect.xmlRestricts electricity trade to occur within the ERCOT, WECC, and IE interconnections. nuclear_USA.xmlUpdates the retirement schedules of the Diablo Canyon and Palisades nuclear power plants. high_cost_offshore_carbon.xmlUpdates the assumed cost of offshore carbon storage. water_supply_constrained_gleeson_5pct.xmlReplaces WaterGAP historical groundwater-depletion estimates with data from the Gleeson dataset and limits groundwater extraction to 5% of the available groundwater in each Superwell grid cell. How to use the data This dataset is designed for use with the IM3 version of GCAM-USA. Download or clone GCAM-USA from the IM3 GCAM GitHub repository at https://github.com/IMMM-SFA/gcam-core and check out the gcam-usa-im3 branch. Place the downloaded folder im3scenarios in the gcam-core/input directory while preserving the provided folder structure.

Energy↗
DOE OSTI2024

A Deep Learning Based Framework to Identify Undocumented Orphaned Oil and Gas Wells from Historical Maps: A Case Study for California and Oklahoma

Undocumented Orphaned Wells (UOWs) are wells without an operator that have limited or no documentation with regulatory authorities. An estimated 310,000 to 800,000 UOWs exist in the United States (US), whose locations are largely unknown. These wells can potentially leak methane and other volatile organic compounds to the atmosphere, and contaminate groundwater. In this study, we developed a novel framework utilizing a state-of-the-art computer vision neural network model to identify the precise locations of potential UOWs. The U-Net model is trained to detect oil and gas well symbols in georeferenced historical topographic maps, and potential UOWs are identified as symbols that are further than 100 m from any documented well. A custom tool was developed to rapidly validate the potential UOW locations. We applied this framework to four counties in California and Oklahoma, leading to the discovery of 1301 potential UOWs across >40,000 km 2 . We confirmed the presence of 29 UOWs from satellite images and 15 UOWs from magnetic surveys in the field with a spatial accuracy on the order of 10 m. This framework can be scaled to identify potential UOWs across the US since the historical maps are available for the entire nation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Pumped-Storage Hydropower using Abandoned Underground Mines (PSH-AUM) as an Innovative Energy Storage Technology for Fossil-Integrated Systems

Pumped-storage hydropower (PSH) provides around 95% of all utility-scale energy storage in the U.S. and globally. Additional deployment of PSH is hampered by excessively long permitting and commissioning requirements and is constrained to locations for which natural topography provides suitable elevation relief between the upper and lower reservoirs (the ΔH challenge). The purpose of this research was to evaluate Pumped-Storage Hydropower using Abandoned Underground Mines (PSH-AUM) as a means to solve the ΔH challenge and initiate the commercialization pathway for a promising new energy storage technology. Four primary tasks were conducted: Screening and ranking of candidate sites for project development; multiphase reservoir modeling to evaluate mine performance; techno-economic analysis and preliminary designs for PSH-AUM systems integrated with fossil-fuel power units; and stakeholder engagement to identify pathways to commercialization of this new technology. Key results were achieved in each of the four primary tasks. Candidate site screening determined that nearly 10,000 underground mines were spatially locatable, of which more than 100 sites appear suitable for integration with existing fossil power plants that are expected to remain in longer-term operation. Mine reservoir models were developed using PNNL’s STOMP simulator and parameterized using candidate site data to evaluate interactions with the surrounding groundwater system and confirm the potential for some mines to accommodate inflows and outflows on the order of 100 m3/s over 8-hour durations (1.6 GWh system) without excessive aqueous pressures. Techno-economic analysis resulted in a project cost optimization scheme to identify key sensitivities and the development of preliminary designs that can minimize overall costs of deployment. Finally, stakeholder engagement with industry, state government, and local economic development leaders confirmed the viability of PSH-AUM as a promising new technology. Our Phase I project results suggest that PSH-AUM technology has a domestic market potential on the order of $100+ Billion with ample space for technology development to commercialization within the next decade.

13 HYDRO ENERGY↗
NASA NTRS2009

Quantifying Parameter Sensitivity, Interaction and Transferability in Hydrologically Enhanced Versions of Noah-LSM over Transition Zones

We use sensitivity analysis to identify the parameters that are most responsible for shaping land surface model (LSM) simulations and to understand the complex interactions in three versions of the Noah LSM: the standard version (STD), a version enhanced with a simple groundwater module (GW), and version augmented by a dynamic phenology module (DV). We use warm season, high-frequency, near-surface states and turbulent fluxes collected over nine sites in the US Southern Great Plains. We quantify changes in the pattern of sensitive parameters, the amount and nature of the interaction between parameters, and the covariance structure of the distribution of behavioral parameter sets. Using Sobol s total and first-order sensitivity indexes, we show that very few parameters directly control the variance of the model output. Significant parameter interaction occurs so that not only the optimal parameter values differ between models, but the relationships between parameters change. GW decreases parameter interaction and appears to improve model realism, especially at wetter sites. DV increases parameter interaction and decreases identifiability, implying it is overparameterized and/or underconstrained. A case study at a wet site shows GW has two functional modes: one that mimics STD and a second in which GW improves model function by decoupling direct evaporation and baseflow. Unsupervised classification of the posterior distributions of behavioral parameter sets cannot group similar sites based solely on soil or vegetation type, helping to explain why transferability between sites and models is not straightforward. This evidence suggests a priori assignment of parameters should also consider climatic differences.

Rosero, Enrique↗
DOE OSTI2025

Evaluating How Climate Adaptation Measures Affect the Interconnected Water‐Energy Resource Systems of the Western United States

Abstract The Western US faces increasing water stress from the impacts of climate change, making it difficult to meet water demands for the region's cities, agriculture, and hydropower generators. Existing literature suggests that climate adaptation measures such as water conservation, cropland retirement, wastewater recycling, and managed aquifer recharge can alleviate some of these challenges. Few analyses, however, compare the relative efficacy and system‐wide effects of these adaptations under different climate projections across the entire Western United States. Here we use a Western US‐wide water systems model to evaluate, by sector and sub‐region, how the widespread implementation of these adaptive measures impacts water demands, water deliveries, and electricity use related to the water system for three different climate projections. We find that wastewater recycling has greater potential to lower unmet indoor water demands than urban indoor water conservation measures. However, when implemented at scale, indoor water conservation reduces electricity use by an average of 683 Terawatt hours while wastewater recycling increases energy use by an average of 721 Terawatt hours, cumulatively from 2020 to 2070. Cropland retirement and aquifer recharge adaptations increase the ability to meet agricultural water demand, increase groundwater storage, and reduce summertime electricity use. While most of these findings are consistent across different climate projections, the benefits of aquifer recharge are sensitive to spatial variation of precipitation. Given the limitations and tradeoffs of each individually, the results suggest that a portfolio of adaptation measures will be needed for a climate‐resilient water and energy future in the Western US. Plain Language Summary The Western US faces increased water stress from the impacts of climate change, making it difficult to meet demands for cities, agriculture, and hydropower facilities. Adaptation measures like water conservation, retiring agricultural lands, recycling wastewater, and storing water underground can address these challenges. However, there is little modeling to understand the impact of these adaptations if they were implemented across the entire Western US while also considering climate change. We use a Western US‐wide water systems model to evaluate how implementing adaptation measures impacts groundwater levels, the ability to meet water demands, and electricity use related to water in three different possible climate futures. We find that recycling wastewater to drinking water standards does a better job of meeting urban water demands than water conservation measures. However, indoor water conservation reduces energy use for water while recycling wastewater increases energy use for water. Retiring agricultural lands and storing water underground both increase the ability to meet agricultural water demands, increase underground water storage, and reduce summertime energy use. No adaptation measures provides benefits across every metric we track. Therefore, multiple adaptation measures will likely be needed to achieve a climate‐resilient future for energy and water in the Western US. Key Points Water‐related climate change adaptations have different benefits and tradeoffs for interconnected water‐energy systems in the Western U.S. Wastewater recycling most increases urban water coverage, but has a significant tradeoff in the form of greater electricity use Cropland retirement and aquifer recharge benefit water and electricity systems most in the summer when those systems are most stressed

Singhal, A↗
DOE OSTI2020

The Role of Subsurface Flow on Evapotranspiration: A Global Sensitivity Analysis

Abstract Water resources are impacted by water‐energy balance fluxes at the land surface, most notably evapotranspiration ( ET ), the largest component of surface energy balance. While integrated hydrologic models show promise in quantifying the nonlinear dynamics at this interface, the model results are plagued by parametric uncertainties. Given the high computational demand of running multiple parameter spaces of these models, little is known about how these uncertainties propagate into land surface processes. In this work, we perform a global sensitivity analysis by computing Sobol and AMAE indices using a surrogate model constructed with a polynomial chaos expansion to assess the impacts of the subsurface physical properties on ET . We do so by modeling a semisynthetic test case. Our results show that the effects of vertical hydraulic conductivity, porosity, and the water retention curve parameter Van Genuchten α mainly control ET . However, we note that while evaporation ( E ) shows behavior similar to the ET with high sensitivity to the parameters controlling the flow in the unsaturated zone, transpiration ( T ) is very sensitive to the saturated zone parameters and groundwater flow, especially during periods without rain. Our results also show that heterogeneities in land cover impact the flow in the saturated zone. Furthermore, our work demonstrates that reduced‐order models can be developed to make integrated models more accessible for rigorous sensitivity analyses and calibration purposes.

Maina, Fadji Zaouna↗
DOE OSTI2022

Delineating Variabilities of Groundwater Level Prediction Across the Agriculturally Intensive Transboundary Aquifers of South Asia

We report the groundwater depletion in South Asia’s Himalayan, transboundary Indus-Ganges-Brahmaputra-Meghna (IGBM) rivers basin is among the highest globally. Given the high irrigation demand and population, groundwater sustainability requires an improved understanding of groundwater systems for the accurate prediction of groundwater levels (GWLs). However, the prediction of groundwater system behaviors is a significant challenge since it is dominated by spatiotemporal and subsurface depth-dependent drivers. Earlier studies that address the challenges are mainly based on the short spatial and temporal extent and/or do not separate the renewable (i.e., shallow) vs nonrenewable (i.e., deeper) groundwater signals. Here, we first identified the variable importance of spatial and depth-dependent drivers on GWL in the IGBM basin. Our results indicate a greater influence of anthropogenic factors (i.e., widespread pumping and increased population) in most parts of the IGBM basin, except in the precipitation-dominated basin of the Brahmaputra. Our next purpose was to delineate a multifactorial approach for GWL prediction using the two most used machine learning models (i.e., support vector machine and feed-forward neural network) in the literature. In general, the machine learning model outputs show a good match in comparison to the GWL from the observation wells (n = 2303 distributed across India and Bangladesh) with some limitations in areas with increased groundwater irrigation. We separately compared the results from shallow (<35 m) and deep (>35 m) observation wells, emphasizing the significance of deep groundwater pumping. Our approach highlights the importance of spatiotemporal to multidepth factors in GWL prediction and can be adopted in other parts of the globe to predict GWLs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Model‐Based Interpretation of Solute Exports and Carbon Partitioning During Shale Weathering in a Mountainous Hillslope

The weathering of sedimentary rocks in high-elevation catchments influences freshwater quality and the global carbon cycle. While individual biogeochemical mechanisms involved in this process are relatively well understood, quantifying their contributions to solute export and carbon fluxes under natural, transient conditions remains challenging. Here, we implement a numerical multidimensional and multiphase model to simulate coupled hydrological and biogeochemical processes in a shale-underlain, snow-dominated hillslope in the Rocky Mountains, Colorado. The model captures the dynamic interplay between soil respiration, mineral weathering, and climate-driven hydrological forcing, reproducing observed soil CO 2 dynamics, groundwater chemistry, and subsurface flow. Our results reveal that seasonal snowmelt enhances carbonate weathering by promoting the infiltration of CO 2 -rich water to depth, while pyrite oxidation is primarily sensitive to low water saturation that facilitates O 2 diffusion through the regolith. Topography modulates the spatial distribution of shale weathering, as steeper slopes enhance lateral drainage, favoring the delivery of reactants to greater depths. While shale weathering at our site acts as a transient carbon sink, with silicates and carbonates buffering acidity and promoting atmospheric CO 2 consumption (1% of soil-derived CO 2 ), the exported dissolved inorganic carbon is predominantly geogenic (∼73%). Consequently, when accounting for long-term marine carbonate precipitation. The current weathering regime represents a net source of carbon to the atmosphere. The oxidation of pyrite and petrogenic organic carbon together release approximately 0.9 mol·m −2 ·yr −1 of CO 2 . Our findings highlight the role of topography, hydroclimate, and the coupling between acid-base reactions in shaping the carbon balance and the solute exports in mountainous critical zones.

carbon cycling↗
DOE OSTI2025

Accelerating multilevel Markov Chain Monte Carlo using machine learning models

Here, this work presents an efficient approach for accelerating multilevel Markov Chain Monte Carlo (MCMC) sampling for large-scale problems using low-fidelity machine learning models. While conventional techniques for large-scale Bayesian inference often substitute computationally expensive high-fidelity models with machine learning models, thereby introducing approximation errors, our approach offers a computationally efficient alternative by augmenting high-fidelity models with low-fidelity ones within a hierarchical framework. The multilevel approach utilizes the low-fidelity machine learning model (MLM) for inexpensive evaluation of proposed samples thereby improving the acceptance of samples by the high-fidelity model. The hierarchy in our multilevel algorithm is derived from geometric multigrid hierarchy. We utilize an MLM to accelerate the coarse level sampling. Training machine learning model for the coarsest level significantly reduces the computational cost associated with generating training data and training the model. We present an MCMC algorithm to accelerate the coarsest level sampling using MLM and account for the approximation error introduced. We provide theoretical proofs of detailed balance and demonstrate that our multilevel approach constitutes a consistent MCMC algorithm. Additionally, we derive the expression for cost reduction due to machine learning model to facilitate cost analysis of the hierarchical sampling algorithm. Our technique is demonstrated on a standard benchmark inference problem in groundwater flow, where we estimate the probability density of a quantity of interest using a four-level MCMC algorithm. Our proposed algorithm accelerates multilevel sampling by a factor of two while achieving similar accuracy compared to sampling using the standard multilevel algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Metabolic Constraints of Organic Matter Mineralization and Metal Cycling

Floodplains are poorly understood and dynamic components of the global carbon cycle that not are well represented in Earth system models. Further, they have a dominant influence on the cycling of important metals within critical transport conduits between surface waters and groundwater. The physical characteristics of floodplains make the hydrology and associated coupled biology and geochemistry particularly responsive to ongoing and impending changes in climate, river management, and land development

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The influence of burn severity on dissolved organic carbon concentrations across a stream network differs based on seasonal wetness conditions

Large, high-severity wildfires in many regions across the globe have increased concerns about their impacts on carbon cycling in watersheds. Altered sources of carbon and changes in catchment hydrology after wildfire can lead to shifts in dissolved organic carbon (DOC) concentrations in streams, which can have negative impacts on aquatic ecosystem health and downstream drinking-water treatment. Despite its importance, post-fire DOC responses remain relatively unconstrained in the literature, and we lack critical knowledge of how burn severity, landscape elements, and climate interact to affect DOC concentrations. To improve our understanding of the impact of burn severity on DOC concentrations, we measured DOC at 129 sites across a stream network extending upstream, within, and downstream of a large, high-severity wildfire in Oregon, USA. We collected samples across the study sub-basin during four distinct seasonal wetness conditions. We used our high-spatial-resolution data to develop spatial stream network (SSN) models to predict DOC across the stream network and to improve our understanding of the controls on DOC concentrations. Spatially, we found no obvious wildfire signal – instead, we observed a pattern of increasing DOC concentrations from the high-elevation headwaters to the sub-basin outlet, while the mainstem maintained consistently low DOC concentrations. This suggests that effects from large wildfires may be “averaged” out at higher stream orders and larger spatial scales. When we grouped DOC concentrations by burn severity group, we observed a significant decrease in the variability of DOC concentrations in the moderate and high burn severity sub-catchments. However, our SSN models were able to predict decreases in DOC concentrations with increases in burn severity across the stream network. Decreases in DOC concentrations were also highly variable across seasonal wetness conditions, with the greatest (–1.40 to –1.64 mg L –1 ) decrease occurring in the high-severity group during the wetting season. Additionally, our models indicated that in all seasons, baseflow index was more influential in predicting DOC concentrations than burn severity was, indicating that groundwater discharge can obscure the impacts of wildfire in a stream network. Overall, our results suggested that landscape characteristics can regulate the DOC response to wildfire. Moreover, our results also indicated that the seasonal timing of sampling can influence the observed response of DOC concentrations to wildfire.

54 ENVIRONMENTAL SCIENCES↗