Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hydro”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Enhanced MPM framework with multipatch isogeometric analysis for geotechnical applications

Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.

58 GEOSCIENCES↗
DOE OSTI2026

Beyond Component Optimization: Systems Level Biodesign for Lanthanide Recovery

Global demand for lanthanides (Ln) is projected to rise sharply over the next decade, while geographically concentrated supply chains and the low concentrations and matrix complexity of secondary feedstocks limit the reach of conventional hydro- and pyrometallurgical separation. Engineered biological systems offer a selective, low-energy alternative, and component-level advances in Ln-binding proteins, AI-designed selective scaffolds, and cell-surface display platforms now rival synthetic chelators in affinity and selectivity. These components, however, remain functionally isolated. Currently, there are no engineered chassis coupling recognition, intracellular trafficking, accumulation, and controlled release into an end-to-end pipeline. Here, we outline how new biodesign strategies and chassis selection must move beyond bioleaching to encompass the full recovery pathway. Achieving this requires integrating AI/ML-guided design, genome-scale build tools, high-throughput phenotyping, and biophysical transport modeling within a Design–Build–Test–Learn cycle tuned to recognition, trafficking, accumulation, and release.

Biodesign↗
DOE OSTI2025

The value of hydropower flexibility for electricity system decarbonization

Hydropower is an abundant, dispatchable, clean energy resource that will play an important role in supporting the clean energy transition. In particular, dispatchable hydropower can provide the operational flexibility that will be required in future systems with high variable renewable energy penetrations. However, the theoretical operational flexibility of hydropower can be restricted in practice by various non-power constraints. In this paper, we quantify how increasing the operational flexibility of dispatchable hydropower resources with reservoirs impacts least-cost generation portfolios and supports power system decarbonization. Specifically, we conduct a capacity expansion analysis of a two-zone system: a hydro-dominated region and a neighboring region with aggressive decarbonization targets that are represented by the United States Pacific Northwest and California respectively. We then introduce a quantifiable index for characterizing the operational flexibility of reservoir hydropower and assess how changes in this metric impact the system-optimal generation portfolio. We find that increasing hydropower flexibility leads to more investment in wind generation, less investment in natural gas generation, lower system costs, and lower system emissions. We further demonstrate a substitution effect between the grid services provided by flexible hydropower operation, increased transmission capacity on a congested line, and energy storage resources. Finally, we show that increasing the operational flexibility of hydropower increases the effective load carrying capability of both hydropower and wind resources. This research supports a more nuanced understanding of how hydropower can support electricity system decarbonization and may motivate reassessing the cost-benefit tradeoffs of non-power constraints that restrict operational flexibility.

Capacity expansion modeling↗
DOE OSTI2025

Hydropower and environmental flow management: System-level trade-offs at Glen Canyon Dam

The research focuses on the Colorado River Basin, specifically examining the Glen Canyon Dam (GCD) and its influence on surrounding aquatic ecosystems. This area is crucial due to its role in hydropower production and its impact on downstream environments, including the Grand Canyon National Park. This study explores the integration of environmental factors into hydro dispatch modeling at GCD to tackle ecological challenges posed by the invasive smallmouth bass (SMB). Utilizing the GTMax SL and SERM models, the research assesses the effects of SMB control experiments on hydropower generation, economic value, and grid stability. The study examines the financial and economic impacts of bypass flows designed to release colder water to prevent SMB spawning, which can significantly reduce hydropower output and increase costs. The research identifies that declining reservoir levels and rising water temperatures in Lake Powell have facilitated SMB spawning, posing a threat to native fish populations like the endangered humpback chub. The findings highlight the importance of adaptive management strategies to balance ecological preservation with hydropower generation amid long-term weather-related challenges. The study underscores the need for comprehensive assessments of flow options to prevent SMB establishment below GCD, considering the broader implications for sediment dynamics and ecological interactions.

Ecological impact assessment↗
DOE OSTI2024

Coupled model for liquid lithium plasma facing components

Numerical analysis provides the design choice and operating window of liquid metal Plasma Facing Components (PFC) concepts. Coupled analysis of boundary plasma together with the surrounding boundary structures is required. Here, to achieve this goal, PPPL is developing a comprehensive multi-physics model for modeling of PFCs in fusion devices. The model includes the fluid-kinetic code SOLPS-ITER and the flow and heat transfer code CFX from ANSYS. SOLPS-ITER was augmented with a liquid metal boundary condition algorithm, allowing direct two-way coupling of the plasma analysis with the two-dimensional analytical slab flow model which includes heat convection in the liquid metal PFC. The target heat flux resulting from this coupled analysis is used as a boundary condition for detailed 3D Computational Fluid Dynamics (CFD) Magneto Hydro Dynamics (MHD) and heat transfer analysis. A new formulation of MHD equations is introduced in the numerical procedure ensuring current conservation of the discretized equations. Results of the 3D analysis are used for final validation of the coupled model. A PFC design where a porous wall is used to stabilize the liquid metal surface, while MHD drive is used to push the liquid metal flow inside the PFC, will be investigated in the regimes where vapor shielding is created for enhanced volumetric plasma heat dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A comparative analysis of numerical approaches for the description of gas flow in clay-based repository systems: From a laboratory to a large-scale gas injection test

There is nowadays a consensus among many countries that geological disposal is a favourable solution for the long-term management. Although different host formations and different barrier systems are under consideration around the world, clay-based materials form an important component for waste isolation in most national programmes. Hence, a good comprehension of the effect of gas flow on the hydro-mechanical behaviour of clay-based soils is essential, both at laboratory and field scale. Task B under the international cooperative project DECOVALEX-2023 has recently shown that, after some enhancement, models can be employed to reproduce laboratory scale tests, even with different sample geometries 37 . However, further work is required to understand whether they can be applied to simulate a large-scale experiment. Up-scaling of models for the advective transport of gas through clay-based low permeable material presents a number of problems related to the difficulty in obtaining consistent hydrogeological parameters and constitutive relationships at both laboratory and field scale. Based on a unique dataset from a large-scale gas injection test (Lasgit) performed at the Äspö Hard Rock Laboratory (Sweden), Task B within DECOVALEX-2023 has explored the refinement of these numerical strategies applied to the simulation of gas flow. Work performed within the task reveals that codes do not need to be substantially modified from the laboratory models to reproduce full-scale tests: indeed, model parameters calibrated and validated at laboratory scale have been applied to predict field scale gas flow at Lasgit, including peak gas pressure and injected cumulative gas volume. By means of (1) the introduction of interfaces between blocks to reflect the experimental configuration and the (2) adjustment of some parameters (e.g., higher permeability), the updated models are able to represent most of the key features observed in the experimental data, even at a large scale.

Tamayo-Mas, E↗
DOE OSTI2024

Leveraging automotive fuel cells can supply zero-emission peak power in the near-term

An increasingly decarbonized yet resilient power grid requires the corresponding build-out of dispatchable zero-emission resources to supply peak power. However, there is a recognized dearth of solutions which can serve multi-day peak demand events both cost-effectively and with near-term deployability. Here, we find that pairing low-cost automotive fuel cells with hydrogen storage in salt caverns can serve as a peaker plant at less than 500 US$/kW at present, a fraction of the cost of conventional fossil fuel-fired peakers. We demonstrate the peaker’s value for long duration storage by comparing it with pumped hydro and assessing its profitability within Texas’ energy-only market region. Although deployment of these peakers is constrained by the presence of salt caverns, we show that a number of sites in the United States and Europe are endowed with suitable salt formations, while utilizing hydrogen storage in pressurized containers could form a location-agnostic peak power solution.

25 ENERGY STORAGE↗
DOE OSTI2026

Structural dynamics of the renewable energy economy: A longitudinal input-output insights for a resilient transition

As countries accelerate their energy transitions, understanding how renewable energy (RE) systems structurally integrate into national economies is essential. This study presents a longitudinal economic input-output (EIO) analysis of the renewable energy sector in South Korea from 2016 to 2022. We develop a novel EIO-based framework that disaggregates the RE sector both by energy source (thermal, hydro, nuclear and renewable) and by industrial function (manufacturing, generation, and services), allowing for a detailed assessment of production dynamics, value-added creation, and import dependency. By quantifying backward and forward linkages and induced economic effects, the analysis reveals persistent structural vulnerabilities in renewable manufacturing and increasing sectoral interdependencies. Results reveal that while the renewable energy sector's production and value-added shares have increased, critical segments remain highly import-dependent, particularly in equipment manufacturing. The analysis highlights systemic gaps in domestic supply chain resilience and offers sector-specific insights for reducing vulnerability and enhancing energy security. Although applied to South Korea as a case study, the proposed framework is designed to be transferable to other national contexts where renewable energy planning requires economic structural insights. The findings offer policy-relevant guidance for enhancing domestic energy resilience and aligning industrial strategy with long-term decarbonization goals.

Economic linkage↗
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 OSTI2024

Modeling nuclear waste disposal in crystalline rocks at the Forsmark and Olkiluoto repository sites – Evaluation of potential thermal–mechanical damage to repository excavations

We conduct coupled thermo-hydro-mechanical modeling of a KBS-3V repository design in crystalline rocks, using data and conditions from the Forsmark in Olkiluoto repository sites in Sweden and Finland. The study focuses on repository performance related to the impact of thermal and hydraulic evolution on the potential for thermal–mechanical damage to underground repository excavations. For the designs and conditions considered at the Forsmark and Olkiluoto repository sites, the simulations show a peak temperature well under the adopted performance target of a 100°C maximum temperature, whereas there is still a high potential for thermal–mechanical damage to the KBS-3V waste deposition holes. The thermal–mechanical damage is much more likely if rock permeability is so low that it delays saturation and swelling of bentonite-clay-based backfill beyond the time for the thermal–mechanical peak, which occurs 50 to 100 years after nuclear waste deposition. We also found that sidewalls of the KBS-3V emplacement tunnels are vulnerable to tensile fracturing due to the combined effect of thermal stressing and backfill swelling. The study highlights a strong interaction between bentonite-based backfill and host rock through capillary suction along with induced rock desaturation. A careful design and selection of the bentonite-clay-based backfill materials for KBS-3V tunnels and deposition holes can facilitate a timely saturation and backfill swelling that in turn can minimize thermal–mechanical damage.

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

Catalytic Upcycling of Polyolefins

The large production volumes of commodity polyolefins (specifically, polyethylene, polypropylene, polystyrene, and poly(vinyl chloride)), in conjunction with their low unit values and multitude of short-term uses, have resulted in a significant and pressing waste management challenge. Only a small fraction of these polyolefins is currently mechanically recycled, with the rest being incinerated, accumulating in landfills, or leaking into the natural environment. Since polyolefins are energy-rich materials, there is considerable interest in recouping some of their chemical value while simultaneously motivating more responsible end-of-life management. An emerging strategy is catalytic depolymerization, in which a portion of the C–C bonds in the polyolefin backbone is broken with the assistance of a catalyst and, in some cases, additional small molecule reagents. When the products are small molecules or materials with higher value in their own right, or as chemical feedstocks, the process is called upcycling. This review summarizes recent progress for four major catalytic upcycling strategies: hydrogenolysis, (hydro)cracking, tandem processes involving metathesis, and selective oxidation. Key considerations include macromolecular reaction mechanisms relative to small molecule mechanisms, catalyst design for macromolecular transformations, and the effect of process conditions on product selectivity. Metrics for describing polyolefin upcycling are critically evaluated, and an outlook for future advances is described.

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

Integrated Carbon Capture and Storage in Hydrogen Production: A Combined Techno-Economic and Life Cycle Assessment

This paper presents a coupled techno-economic and life cycle assessment of “blue” hydrogen to be produced at a hydrogen facility through steam methane reforming (SMR) equipped with carbon capture and storage (CCS). Blue hydrogen was modeled in ChemCAD, while an integrated asset model represented the carbon capture and storage chain. An unabated carbon dioxide (CO 2 ) configuration release 11.99 kgCO 2 -eq/ kgH 2 . Capturing ≥95% of the CO 2 stream lowers the carbon footprint to 6.59 kgCO 2 -eq/kgH2 but raises the levelized cost of hydrogen (LCOH) from $\$$1.82/kgH 2 (no CO 2 capture) to $\$$3.22/kgH 2 ; the U.S. 45Q tax credit reduces it to $\$$2.59/kgH 2 . Incorporating CCS reduces the levelized net present value from $\$$0.87/kgH 2 to $\$$0.74/kgH 2 , owing to additional capture, transport, and storage costs. Supplying SMR with low-carbon electricity, especially nuclear, wind, or hydro, delivers the lowest carbon footprint relative to geothermal or grid mixes. Sensitivity analysis identifies that hydrogen sales price, internal rate of return, and CCS cost as the strongest economic levers, while electricity demand dominates residual lifecycle emissions. The results underscore a clear trade-off; substantial CO 2 reductions are achievable, but only with higher production costs, making supportive policy instruments, access to clean power, and robust hydrogen markets essential for large-scale deployment of blue hydrogen.

carbon capture↗
DOE OSTI2024

Combined Cyclohexene Dehydrogenation and Oxygenation over VO x /CeO 2 with O 2 and Vaporized H 2 O 2 as Oxidants

Selective oxidation of hydrocarbons over metal oxide catalysts frequently invokes Mars-van Krevelen (MvK) mechanisms involving lattice oxygens. At typical reaction temperatures, oxidative dehydrogenation reactions predominate and compete with total oxidation; production of more useful oxygenates is challenging. However, because regeneration of lattice oxygens has been previously hypothesized to involve transient formation of surface hydroperoxides or peroxides, MvK cycles should be capable of forming products derived from both lattice oxygen and surface (hydro)peroxide species, under certain conditions. In this report, we study the low temperature reaction of cyclohexene with O 2 over CeO 2 -supported VO x catalysts. Below 150°C and for higher VO x loading, cyclohexene converts to benzene and trans-1,2-cyclohexanediol with no appreciable CO 2 formation. Control experiments show similar behavior in other supported vanadia catalysts. Conversely, using vaporized H 2 O 2 as oxidant exclusively gives 1,2-cyclohexanedione, suggesting that MvK reoxidation cycles do not directly involve H 2 O 2 . Furthermore, these findings contribute to the current understanding of selective oxidation mechanisms, and demonstrate routes to the formation of industrially useful oxygenates with O 2 or vaporized H 2 O 2 over metal oxide catalysts.

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

Integrated Effects of Site Hydrology and Vegetation on Exchange Fluxes and Nutrient Cycling at a Coastal Terrestrial‐Aquatic Interface

Abstract The complex interactions among soil, vegetation, and site hydrologic conditions driven by precipitation and tidal cycles control the biogeochemical transformations and bi‐directional exchange of carbon and nutrients across the terrestrial–aquatic interfaces (TAIs) in coastal regions. This study uses a highly mechanistic model, Advanced Terrestrial Simulator (ATS)‐PFLOTRAN, to explore how these interactions affect exchanges of materials and carbon and nitrogen cycling. We used a transect in the Chesapeake Bay region that spans zones of open water, coastal wetland, transition, and upland forest. We designed several simulation scenarios to parse the effects of the individual controlling factors and the sensitivity of carbon cycling to reaction rate parameters derived from laboratory experiments. Our simulations reveal an active zone for carbon cycling under the transition zones between the wetland and the upland. Evapotranspiration is found to enhance the exchange fluxes between the surface and subsurface domains, resulting in a higher dissolved oxygen concentration in the TAIs. The transport of organic carbon derived from plant leaves and roots provide an additional source of organic carbon needed for the aerobic respiration and denitrification processes in the TAIs. The variability in reaction rate parameters associated with microbial activities is also found to play a dominant role in controlling the heterogeneity and dynamics of the simulated redox conditions. This modeling‐focused exploratory study enabled us to better understand the complex interactions among soil, water and microbes that govern the hydro‐biogeochemical processes at the TAIs, which is an important step toward representing coastal ecosystems in larger‐scale Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Modeling the Effects of Artificial Drainage on Agriculture-Dominated Watersheds Using a Fully Distributed Integrated Hydrology Model

In agriculture-dominated watersheds where natural drainage is poor, agricultural ditches (narrow engineered channels) and tile drains (perforated pipes) are widely employed to enhance surface and subsurface drainage, respectively. Despite their relatively small scale, these features exert substantial control over the hydro-biogeochemical function of watersheds and their effects need to be represented in the models. We introduce a novel strategy to incorporate the effects of artificial agricultural drainage into a fully distributed basin-scale integrated surface-subsurface hydrology models. In our approach, narrow agriculture ditches for surface drainage are resolved efficiently using ditch-aligned computational meshes that are hydrologically conditioned to ensure connectivity in the stream/ditch network. For tile drainage in the subsurface, we use the physically based Hooghoudt's drainage equation as a subgrid model and route the water drained through tiles to the nearest ditch. Without site-specific calibration, this model reproduced observed streamflow in the Portage River Watershed (>1,000 km 2 ) as recorded by a USGS gauge with good accuracy (normalized KGE = 0.81) and outperformed a calibrated SWAT model (normalized KGE = 0.68). Numerical experiments confirm that artificial drainage reduces surface inundations and effectively controls the water table. At the watershed scale, artificial drainage increases baseflow but has little effect on watershed discharges above the 90th percentile. The strong physical underpinnings and reduced need for calibration allow us to study the impacts of artificial drainage on distributed hydrological response in terms of fluxes and states and provide a platform for investigating watershed-scale nutrient transport.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantifying Streambed Grain Size, Uncertainty, and Hydrobiogeochemical Parameters Using Machine Learning Model YOLO

Abstract Streambed grain sizes control river hydro‐biogeochemical (HBGC) processes and functions. However, measuring their quantities, distributions, and uncertainties is challenging due to the diversity and heterogeneity of natural streams. This work presents a photo‐driven, artificial intelligence (AI)‐enabled, and theory‐based workflow for extracting the quantities, distributions, and uncertainties of streambed grain sizes from photos. Specifically, we first trained You Only Look Once, an object detection AI, using 11,977 grain labels from 36 photos collected from nine different stream environments. We demonstrated its accuracy with a coefficient of determination of 0.98, a Nash–Sutcliffe efficiency of 0.98, and a mean absolute relative error of 6.65% in predicting the median grain size of 20 ground‐truth photos representing nine typical stream environments. The AI is then used to extract the grain size distributions and determine their characteristic grain sizes, including the 10th, 50th, 60th, and 84th percentiles, for 1,999 photos taken at 66 sites within a watershed in the Northwest US. The results indicate that the 10th, median, 60th, and 84th percentiles of the grain sizes follow log‐normal distributions, with most likely values of 2.49, 6.62, 7.68, and 10.78 cm, respectively. The average uncertainties associated with these values are 9.70%, 7.33%, 9.27%, and 11.11%, respectively. These data allow for the computation of the quantities, distributions, and uncertainties of streambed HBGC parameters, including Manning's coefficient, Darcy‐Weisbach friction factor, top layer interstitial velocity magnitude, and nitrate uptake velocity. Additionally, major sources of uncertainty in grain sizes and their impact on HBGC parameters are examined.

58 GEOSCIENCES↗
DOE OSTI2025

Co‐Evolution of Specific Stiffness and Permeability of Rock Fractures Offset in Shear

Abstract Fractures and faults represent planes of weakness and compliance in rock masses that serve as focal points for both microearthquakes and fluid transport, with seismicity and permeability evolution closely linked. Contact stiffness is highly stress‐sensitive and directly influences permeability. We explore the co‐evolution of specific stiffness and permeability of rough fractures under normal stress and shear offset using numerical simulations. Individual rough fractures are represented by variable amplitude (Root mean square) and wavelength ( λ ) using a granular mechanics model. Contacting rough surfaces are mated, offset in shear, and then compacted in displacement mode. The compacting fractures generate stress‐dependent changes in contact porosity, which govern both permeability and stiffness evolution. We establish a universal dimensionless relationship linking specific stiffness and permeability that inherently incorporates the effects of surface roughness, shear offset, and microcracking. The observed cracking effect—where local stress redistribution and pressure‐driven microcrack propagation dynamically alter the aperture field—introduces a nonlinear permeability response at high stress. Increased roughness amplitude and larger shear offsets reduce stiffness while dampening permeability sensitivity to stress, demonstrating a strong interplay between surface texture and hydro‐mechanical behavior. While the model captures this behavior effectively, deviations emerge at very low porosities due to extreme aperture sensitivity in this limit.

He, Xinxin [Energy and Mineral Engineering Pennsyl↗
DOE OSTI2024

BioRT‐HBV 1.0: A Biogeochemical Reactive Transport Model at the Watershed Scale

Abstract Reactive Transport Models (RTMs) are essential tools for understanding and predicting intertwined ecohydrological and biogeochemical processes on land and in rivers. While traditional RTMs have focused primarily on subsurface processes, recent watershed‐scale RTMs have integrated ecohydrological and biogeochemical interactions between surface and subsurface. These emergent, watershed‐scale RTMs are often spatially explicit and require extensive data, computational power, and computational expertise. There is however a pressing need to create parsimonious models that require minimal data and are accessible to scientists with limited computational background. To that end, we have developed BioRT‐HBV 1.0, a watershed‐scale, hydro‐biogeochemical RTM that builds upon the widely used, bucket‐type HBV model known for its simplicity and minimal data requirements. BioRT‐HBV uses the conceptual structure and hydrology output of HBV to simulate processes including advective solute transport and biogeochemical reactions that depend on reaction thermodynamics and kinetics. These reactions include, for example, chemical weathering, soil respiration, and nutrient transformation. The model uses time series of weather (air temperature, precipitation, and potential evapotranspiration) and initial biogeochemical conditions of subsurface water, soils, and rocks as input, and output times series of reaction rates and solute concentrations in subsurface waters and rivers. This paper presents the model structure and governing equations and demonstrates its utility with examples simulating carbon and nitrogen processes in a headwater catchment. As shown in the examples, BioRT‐HBV can be used to illuminate the dynamics of biogeochemical reactions in the invisible, arduous‐to‐measure subsurface, and their influence on the observed stream or river chemistry and solute export. With its parsimonious structure and easy‐to‐use graphical user interface, BioRT‐HBV can be a useful research tool for users without in‐depth computational training. It can additionally serve as an educational tool that promotes pollination of ideas across disciplines and foster a diverse, equal, and inclusive user community.

Sadayappan, Kayalvizhi↗