Search NASA⌕ Search

SEARCH · Search NASA

Results for “interface kinetics”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Rigor and Reproducibility in Electrocatalysis: Best Practices for Operando Studies

Operando measurements have rapidly expanded the scope of electrocatalysis by enabling direct observation of catalytic interfaces under working conditions and by linking structural, compositional, and spectroscopic observables to activity and selectivity. However, the growth of operando methods has outpaced the adoption of broadly shared experimental standards, creating persistent challenges in reproducibility, interpretation, and comparison across laboratories and platforms. This perspective synthesizes discussions from the 2025 National Science Foundation Workshop on Rigor and Reproducibility in Electrocatalysis and outlines a practical framework for the rigorous use of operando measurements in electrocatalysis. We highlight three recurring needs: careful implementation of complex methods to avoid overinterpretation; recognition that (subtle) differences in reactor architecture, hydrodynamics, and electrical boundary conditions can alter apparent kinetics and selectivity; and transparent reporting standards that enable meaningful cross-comparison without constraining measurement-specific cell innovation. Focusing on widely used techniques (including X-ray and vibrational spectroscopies, mass spectrometry, and electron microscopy), we discuss technique-specific pitfalls, cross-validation strategies, and recurring platform-agnostic considerations such as mass transport, current distribution, temporal-resolution mismatches, and catalyst evolution. This Perspective aims to strengthen the mechanistic inference and improve the reproducibility, comparability, and predictive value of operando electrocatalysis research.

X-ray absorption spectroscopy↗
DOE OSTI2026

Elucidating key reducing species beyond ions in hydrogen plasma smelting reduction of iron ore

Hydrogen plasma smelting reduction (HPSR) of iron ore has attracted significant attention over the past decade due to its high-temperature operation, rapid plasma mediated reduction kinetics, and simpler density-based separation of molten iron product, compared to H2-based solid-state reduction. All of these attributes enable processing of low-grade ores for downstream use in electric-arc furnaces, as virgin iron with low gangue content is required for high quality steel and improved furnace operation. While positive ions exist within the plasma arc, this work demonstrates that near the anodic ore surface, hydrogen radicals and vibrationally excited hydrogen species dominate and their densities correlate well with observed reduction rates. Species concentrations in the transferred plasma arc and at the plasma-ore interface are evaluated using coupled thermal plasma and near-wall non-equilibrium plasma models. The thermal plasma model is validated against experimental voltage data and spectroscopic measurements of plasma temperature and density for varying current inputs. Modeling of the near surface thermochemical non-equilibrium and micrometer scale anode sheath layer reveals, in addition to the expected H + , significant concentrations of ArH + and H$^+_3$ ions, typically not observed in thermal plasmas under thermodynamic equilibrium. Our results show that the inverted sheath structure at the anodic ore surface strongly suppresses reactive positive ion fluxes, while non-equilibrium electron-impact processes generate abundant hydrogen radicals and vibrationally excited species. These findings highlight the critical role of non-equilibrium effects in hydrogen arc-driven iron ore reduction and advance understanding beyond prevailing hypotheses centered on hydrogen ion-driven mechanisms.

08 HYDROGEN↗
DOE OSTI2025

Charge-Transfer States Enable Spin-Selective Formation of Quartet State Qudits in Luminescent Radical-Chromophore Dyads

Luminescent tris(2,4,6-trichlorophenyl)methyl (TTM) radicals are promising doublet emitters that have been used to generate spin-optical interfaces in molecular electron spin qudits. In particular, photoexcitation of covalent TTM-chromophore dyads can generate quartet spin states that can be detected via photoluminescence. However, the mechanism of quartet spin initialization is complicated by competing charge-transfer (CT) and energy-transfer (EnT) pathways between the radical and the chromophore. Here, we demonstrate the role of CT intermediates by engineering a covalent dyad of TTM and naphthalene-(1,4:5,8)-bis(dicarboximide) (NDI), a well-known electron acceptor chromophore. Transient absorption spectroscopy indicates that photoexcitation of TTM results in ultrafast electron transfer from 2* TTM to NDI, which kinetically outcompetes EnT. Electron paramagnetic resonance spectroscopy reveals that charge recombination proceeds via spin–orbit charge-transfer intersystem crossing to generate 3* NDI, followed by spin mixing with 2TTM to form the quartet state. Furthermore, this process uniquely spin-polarizes the quartet state, providing fundamental design principles to purify the initial wave function for quantum information processing.

Charge transfer↗
DOE OSTI2025

Synthesis challenges, thermodynamic stability, and growth kinetics of La–Si–P ternary compounds

Although many new compounds have been recently predicted with the help of machine learning, the successful experimental synthesis of these compounds remains challenging. Computational insights about the thermodynamic stability and phase formation kinetics among the ground state and competing metastable phases are highly desirable to rationalize and attempt to overcome synthesis challenges experimentally. In this work, we explore synthetic challenges within ternary La–Si–P compounds through feedback between experimental and computational studies. We discuss the experimental challenges in forming three computationally predicted ternary phases (La 2 SiP, La 5 SiP 3 , and La 2 SiP 3 ). To understand the synthetic challenges, we performed molecular dynamics (MD) simulations using an accurate and efficient artificial neural network machine learning (ANN-ML) interatomic potential. We study the phase stability and formation kinetics of these ternary phases in relation to the reported and synthesized La 2 SiP 4 phase. While the growth of the La 2 SiP 4 phase can be reproduced by our MD simulation, our results indicate that the rapid formation of a Si-substituted LaP crystalline phase is a major barrier to the synthesis of the predicted La 2 SiP, La 5 SiP 3 , and La 2 SiP 3 ternary compounds, agreeing well with experimental observations. Our simulations also suggest that there is a narrow temperature window in which the La 2 SiP 3 phase can be grown from the solid–liquid interface.

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

In Situ Characterization of Surface Recombination in p-Si/SiO x Based Photoelectrochemical Cells

Time-resolved infrared (TRIR) and electrochemical impedance spectroscopy (EIS) were utilized to quantify surface states present at silicon oxide (SiO x ) protected crystalline p-Si electrolyte interfaces. The primary goal was to identify p-Si/SiO x photoelectrodes with both low surface recombination rates and efficient multi-electron transfer to an acceptor present in the external electrolyte. Three SiO x layers were investigated: native oxide (nOx), chemical oxide (cOx), and rapid thermal annealed (RTA) thermal oxide (tOx). Comparative study with [Ru(bpy) 3 ](PF 6 ) 2 as the electron acceptor indicated that tOx was most optimal with a small effective recombination rate, multi-electron transfer capability, and photovoltage of 500 ± 50 mV. A secondary goal was to analyze the surface recombination rates with the Shockley–Read–Hall (SRH) kinetic model. Two surface states were identified from this analysis, one closer to the CB edge (V t,1 ) and the other near the midgap (V t,2 ). EIS and SRH analyses revealed that a forming gas (5% H 2 /N 2 ) anneal (FGA) decreased surface recombination for tOx and nOx through a lower density of surface states. In the case of tOx, the infrared data indicated that V t,2 was completely removed. Here, the energetic positions of the band edges were correlated with the surface state density; low densities corresponded to more favorable potentials for inversion layer formation, which is expected to be most optimal for photocatalysis. Collectively this study indicates that the free carrier dynamics provided by TRIR represent a powerful in situ probe of the band edge and the surface state energetics in silicon based photoelectrochemical cells.

Electrochemical Impedance Spectroscopy↗
DOE OSTI2025

Single-Crystallization of O3-Type Layered Oxide Cathode for Na-Ion Battery

The development of high-energy-density Na-ion batteries places significant demands on single-crystal layered oxide cathodes, especially for further high-voltage, solid-state battery scenarios. In the O3-type structure, due to the original sluggish Na ion diffusion kinetics (approximately 1 order of magnitude lower than that of Li-ion), and further hindrance against diffusion kinetics caused by single-crystal architecture, these inherent defects lead to the decline in the electrochemical performance. Herein, we demonstrated that the single crystallization of O3-type NaNi 1/3 Fe 1/3 Mn 1/3 O 2 cathode (d 50 = 5.04 μm) aggravates surface-to-bulk phase inhomogeneity distribution, which is attributed to the uneven Na ions extraction. Moreover, the Na-depletion of the surface/shell region not only aggravates Na ion diffusion resistance but also leads to a higher valence state of transition-metal elements (e.g., Ni/Fe) near the surface of the single-crystal particle, which further compromises the cathode-electrolyte interface stability. Furthermore, not limited to revealing the challenges, tuning the particle size and moderating quasi-single-crystal strategies have been proven to effectively mitigate the negative uneven distributions of Na ions, phases, and valence/oxidative states, resulting in efficient modification for single crystallization of Na-layered oxide cathodes.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Equivalent equilibrium ion distributions in ion sources and beams

Ion sources and ion beams often contain ions of more than one charge-to-mass ratio, due to the presence of multiple ionization states and/or multiple atomic and molecular ion species in the source. We show that, under certain conditions, analysis of these systems can be simplified using the notion of “equivalent equilibrium ion distributions” (EEIDs). These distributions have identical ion trajectories, space-charge distributions, and electric-field distributions, but can be composed of ion species with arbitrary charge-to-mass ratios. We use this result to obtain multi-species generalizations for the perveance of an ion diode, which relates the diode voltage to the extracted beam current, and for the perveance of a drifting ion beam, which determines the axial length-scale for beam expansion due to the self-electric field. Applied to a quasi-neutral plasma ion source containing cold ions and hot electrons, we find that the location of the beam-plasma interface (meniscus), where ions are extracted from the plasma, depends on the electron density and ion energy per unit charge at the interface, independent of the ion composition. For certain types of numerical simulations, EEIDs allow one to calculate the beam optics of an intense multi-species ion source by doing a simulation for a single ion species with an artificial charge-to-mass ratio chosen to minimize computation time. We illustrate these results using the 3D particle-in-cell, finite-element code Aleph. The simulations show how ions in a plasma source self-consistently accelerate to the supersonic speed needed for steady, radially expanding flow. Multi-species equilibria are subject to streaming instabilities because, while the different species have the same kinetic energy per unit charge, their velocities scale with the square root of the charge-to-mass ratio. The simulations show that a strong ion streaming instability can occur in a multi-species plasma ion source.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Chemical and Structural Insights into Solid Electrolyte Interphase Evolution for Sodium Metal Electrodes

The solid electrolyte interphase (SEI) critically governs the reversibility of sodium metal batteries, through dynamically mediating ion transport and interfacial reactions. However, its kinetic evolution under operating conditions, and how it influences interfacial stability, remains poorly understood. Here, in this study, we reveal that the SEI undergoes coupled chemical and mechanical changes during sodium plating and stripping, leading to spatial and temporal heterogeneity that drives interfacial degradation. Synchrotron operando grazing-incidence wide-angle X-ray scattering and soft X-ray absorption spectroscopy capture the sequential formation and dissolution of inorganic SEI phases (NaF, NaH, NaOH, Na 2 PO 3 F), accompanied by depth-dependent alterations in organic SEI components. Mesoscale modeling connects this evolving SEI heterogeneity to localized current density fluctuations and stress accumulation at the Na interface, identifying pathways to electrically isolated sodium formation. These findings show that SEI instability fundamentally limits reversibility in sodium metal batteries, and that controlling SEI chemistry–mechanics coupling is essential to achieving its durability.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Influence of zwitterionic amphiphilic copolymers on heterogeneous gypsum formation: A promising approach for scaling resistance

Here, this study aims to investigate the influence of zwitterionic amphiphilic copolymers (ZACs) in the nucleation and growth of heterogeneous CaSO4 at the zwitterion-water interface, which is crucial for the prevention of mineral scaling and consequent downtime or suboptimal performance in industries like membrane desalination, heat exchangers, and pipeline transportation. In situ grazing incidence small angle X-ray Scattering (GISAXS), and quartz crystal microbalance with dissipation (QCM-D) techniques were used to analyze the evolution of CaSO 4 particles on two new ZAC coatings: poly-(trifluoroethyl methacrylate-random-sulfobetaine methacrylate) (PTFEMA-r-SBMA, or PT:SBMA) and poly(trifluoroethyl methacrylate-random-2-methacryloyloxyethyl phosphorylcholine) (PTFEMA-r-MPC, or PT:MPC). The results showed that PT:MPC coatings promoted nucleation but inhibited crystal growth, resulting in slower overall reaction kinetics on PT:MPC coatings compared to PT:SBMA coatings. Interfacial interactions involving the substrates, sulfate minerals, and ions were examined, revealing that calcium ion adsorption, primarily governed by electrostatic attraction, played a crucial role in the nucleation and growth processes on both ZAC coatings. The crystal characterization revealed a phase transition from bassanite to gypsum on both ZAC coatings, suggesting that these zwitterionic materials can influence the mineral phase of heterogeneously formed CaSO 4 crystals. These findings enhance our understanding of the fundamental mechanisms underlying heterogeneous CaSO 4 scaling in the presence of zwitterionic materials.

42 ENGINEERING↗
DOE OSTI2026

Effects of Polymer Morphology on Solvent and Catalyst Accessibility during Polyethylene and Polystyrene Autoxidation

Efficient catalytic deconstruction of plastics requires facile solvent and catalyst access to polymer substrates to minimize mass transfer effects. Autoxidation using Co(II) acetate, Mn(II) acetate, and a radical carrier in acetic acid is a promising strategy to deconstruct mixed plastic waste, yet the role of polymer morphology in governing solvent and catalyst accessibility remains poorly understood. Here, in situ simultaneous small- and wide-angle X-ray scattering (SAXS/WAXS), complemented by X-ray fluorescence (XRF) imaging and high-pressure differential scanning calorimetry (DSC), were used to elucidate interactions between acetic acid, a Co/Mn catalyst solution, and semicrystalline polyethylene (PE) and amorphous polystyrene (PS) from room temperature to 160 °C. In PE, acetic acid and catalyst access were confined to amorphous regions and cryomilled particle interfaces at room temperature, while crystalline lamellae remained intact after soaking for up to 34 h. Increasing temperature enabled solvent uptake into PE, followed by solvent-assisted softening above 100 °C, and a modest melting-point depression that removed lamellar transport barriers upon melting. Conversely for PS, acetic acid penetrated the glassy polymer without inducing chain mobility until the glass transition was reached, above which the observed structural changes were consistent with enhanced segmental mobility which enabled bulk penetration. These results suggest that polymer morphology and thermally activated physical transitions arising from diffusion and polymer–solvent interactions can influence whether autoxidation of plastics is transport-limited or kinetically controlled, providing a framework for aligning reaction conditions with reaction outcomes.

09 BIOMASS FUELS↗
DOE OSTI2025

Decoding anomalous grain growth at room temperature during pressure-induced phase transformations

Significant grain growth is observed during the high-pressure phase transformations (PTs) at room temperature in various materials. The main focus here is grain growth from a few hundred nanometers to 10 μm within an hour during α → ω PT in Zr. No existing theory explains this phenomenon since without PT, Zr nanocrystals do not grow at room temperature even for up to 10 years. Here, in this study, a multistep mechanism for the grain growth during α → ω PT in Zr is suggested. Phase interfaces (PI) and grain boundaries (GBs) coincide and move together as a combined PI-GBs under the action of the combined thermodynamic driving force. Such a combined motion changes the diffusional grain growth mechanism to the transformational one and the martensitic mechanism of PT to a reconstructive one via an intermediate disordered phase. The primary condition is that the GB energy of the ω phase is smaller than that of the α phase, which promotes the nucleation of ω-Zr and is consistent with the absence of the reverse PT and reduction in the PT pressure with reducing grain size. Several intermediate steps for such motion are suggested and justified kinetically. Nonhydrostatic stresses due to volume reduction in the growing ω grain promote continuous growth of the existing ω grain instead of a new nucleation at other GBs. In situ synchrotron Laue diffraction experiments confirm the main predictions of the theory. The suggested mechanism provides a new insight into synergistic interaction between PTs and microstructure evolution.

anomalous grain growth during phase transformation↗
DOE OSTI2024

Explosively driven Richtmyer–Meshkov instability jet suppression and enhancement via coupling machine learning and additive manufacturing

The ability to control the behavior of fluid instabilities at material interfaces, such as the shock-driven Richtmyer–Meshkov instability, is a grand technological challenge with a broad number of applications ranging from inertial confinement fusion experiments to explosively driven shaped charges. In this work, we use a linear-geometry shaped charge as a means of studying methods for controlling material jetting that results from the Richtmyer–Meshkov instability. A shaped charge produces a high-velocity jet by focusing the energy from the detonation of high explosives. The interaction of the resulting detonation wave with a hollowed cavity lined with a thin metal layer produces the unstable jetting effect. By modifying the characteristics of the detonation wave prior to striking the lined cavity, the kinetic energy of the jet can be enhanced or reduced. Modifying the geometry of the liner material can also be used to alter jetting properties. We apply optimization methods to investigate several design parameterizations for both enhancing or suppressing the shaped-charge jet. This is accomplished using 2D and 3D hydrodynamic simulations to investigate the design space that we consider. We also apply new additive manufacturing methods for producing the shaped-charge assemblies, which allow for the experimental testing of complicated design geometries obtained through computational optimization. We present a direct comparison of our optimized designs with experimental results carried out at the High Explosives Application Facility at Lawrence Livermore National Laboratory.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simulating the Phonon Collection Efficiency in Kinetic Inductance Phonon-Mediated Detectors

Kinetic inductance phonon-mediated (KIPM) detectors use microwave kinetic inductance detectors (MKIDs) to read out phonon signals in the substrate. They are a promising class of detectors to be used in light dark matter (DM) searches due to their potential eV-scale sensitivity and native frequency-domain multiplexability. In order to improve upon the design of these detectors, simulations are needed to understand the effects of detector design on measurable physical parameters including the phonon collection efficiency, $\eta_{ph}$, which is defined as the ratio of phonon energy detected by the sensitive target element to the incident energy deposited in the substrate. This work simulates the phonon collection efficiency for two existing detector architectures to study phonon absorption at substrate-sensor interfaces and phonon loss to non-sensitive detector elements.

Dang, Stella Q.↗
DOE OSTI2024

Simulating the Phonon Collection Efficiency in KIPMDs

Kinetic inductance phonon-mediated (KIPM) detectors are superconducting microcalorimeters that use microwave kinetic inductance detectors (MKIDs) to read out phonon signals in the device substrate. In order to improve the design of these detectors, we need to understand the effect of various detector design elements on the physical processes that take place within the detector through simulation efforts. One figure of merit for KIPM detectors is the phonon collection efficiency, $\eta_{ph]$, defined as the ratio of phonon energy detected by the sensitive element in the detector and the energy deposited by incident phonons in the substrate. Comparing the phonon collection efficiency obtained in simulation and in experiment for the same detector geometry provides insight to physical processes in the detector, such as phonon absorption at substrate-sensor interfaces and phonon loss to non-sensitive detector elements. This work models the phonon collection efficiency in two KIPM detectors with different geometries, one located at the Fermi National Accelerator Laboraory, the other located at the SLAC National Accelerator Laboratory.

Dang, Stella Q.↗
DOE OSTI2026

Unveiling the High‐Voltage Reactivity and Gas Evolution With Aluminum‐Based Chloride and Oxychloride Catholytes in Solid‐State Sodium Batteries

All-solid-state sodium batteries (ASSBs) employing halide solid electrolytes (SEs) offer a cost-effective and energy-dense alternative to conventional liquid electrolyte systems. However, their high voltage (>4 V vs. Na/Na + ) performance remains limited by interfacial instability between the cathode active material (CAM) and the SE. We present here the electrochemical and interfacial behaviors of crystalline NaAlCl 4 and amorphous sodium–aluminum–oxychloride (NACO) SEs when combined with NaNi 0.5 Mn 0.5 O 2 cathode. While oxygen incorporation in NACO enhances ionic conductivity by nearly three orders of magnitude relative to NaAlCl 4 , it does not improve high-voltage cycling stability. Cells employing NACO exhibit accelerated capacity fade, increased cell impedance growth, and intrinsic oxygen evolution above 4.5 V vs. Na 3 Sn, as revealed by operando electrochemical mass spectrometry. In contrast, the NaAlCl 4 -based cells show no detectable gas release, underscoring their superior high-voltage stability and safety. Time-of-flight secondary-ion mass spectrometry confirms the formation of Al─O and Ni/Mn─Cl species, respectively, in the SE and CAM, indicating redox-driven anion exchange that contributes to kinetic hindrance of high-voltage phase transitions. The findings establish that while oxygen incorporation enhances ionic transport, it can compromise interfacial stability, suggesting pure chloride SEs may offer a more robust and intrinsically safer pathway for developing high-energy ASSBs.

25 ENERGY STORAGE↗
DOE OSTI2024

Oxidation Temperature-Dependent Electrochemical Doping of WO 3 Deposited via Atomic Layer Deposition

Silicon-based photoelectrochemical devices show promise for the performance of light-driven CO 2 reduction but suffer from instability under photoelectrochemical conditions relevant to CO 2 reduction. Coating silicon electrodes with thin layers of metal oxides has shown promise to passivate unstable silicon surfaces, and many different metal oxides can be deposited on silicon using various techniques. In this study, we investigate the fundamental photoelectrochemical performance of WO 3 -coated silicon photoelectrodes, which were generated by oxidation of W-metal films deposited via atomic layer deposition on both degenerately doped (nSi + ) and low-doped (pSi) silicon. Two different oxidation temperatures were investigated (400 and 600 °C), and it was found that the monoclinic phase of WO 3 predominates at both temperatures but that more grain boundaries are present in the 600 °C film. From X-ray photoelectron spectroscopy, the stoichiometry of both films was found to be 1:3 W:O, and low electron energy loss experiments indicate band gaps of 3.0 and 3.1 eV for 400 and 600 °C films, respectively. Cyclic voltammetry experiments showed that the electron transfer kinetics increased after continued redox cycling, particularly for the material produced at 400 °C. X-ray photoelectron spectra suggest that the observed increase in electrode conductivity is due to the formation of oxygen vacancies in the film. Electrochemical impedance spectroscopy indicated that charge transport through the films was impacted by the grain boundaries that formed during oxidation of the film. Photoelectrochemical studies on pSi/WO 3 electrodes were highly variable, only producing a photocurrent and photovoltage with some samples. Our best sample, formed at 400 °C, produced a photovoltage of 180 mV, which is lower than what has previously been reported for WO 3 -coated silicon (500 mV). We hypothesize that the variability in photoelectrochemical experiments arose from a roughened WSiO x interface that is generated during film preparation. WO 3 shows promise as a metal oxide coating for silicon, but our results suggest that formation of a high-quality interface between Si and WO 3 is vital for best performance.

Charge transfer↗
DOE OSTI2025

Polyethylene Glycol Surface Modification and Polythiophene Side-Chain Chemistry: A Combined Strategy toward High-Capacity Lithium-Ion Battery Anodes

In the development of high-capacity lithium-ion batteries (LIBs), the combined optimization of active material interfaces and polymer binder chemistry plays a critical role in improving electrode performance and longevity. This work explores a dual design strategy incorporating polyethylene glycol (PEG) surface modification and carboxylated polythiophene side-chain tailoring to enhance the electrochemical behavior of magnetite (Fe 3 O 4 )-based anodes. PEG is employed to improve interfacial stability, while carboxylated polythiophene binders with varying alkyl side-chain lengths─poly[3-(potassium-4-butanoate)thiophene-2,5-diyl] (P3KBT), poly[3-(potassium-5-pentanoate)thiophene-2,5-diyl] (P3KPT), and poly[3-(potassium-6-hexanoate)thiophene-2,5-diyl] (P3KHT)─are used to modulate molecular interactions and ion transport. Among these three analogs, the PEG–Fe 3 O 4 –P3KHT electrode exhibits superior ion-transfer kinetics, the highest capacity retention, and the lowest charge-transfer resistance after extended cycling. Compared to their non-PEG analogs, PEG-coated electrodes demonstrate enhanced structural integrity and electrochemical behavior, emphasizing the synergistic effects of surface modification and side-chain chemistry. These findings highlight the importance of interfacial interactions and molecular design in achieving robust and high-performance composite anodes for next-generation LIBs.

Fe3O4↗
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↗