SEARCH · Engineering Papers

Results for “Inference mechanisms”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Electrochemical Random-Access Memory: Progress, Perspectives, and Opportunities

Non-von Neumann computing using neuromorphic systems based on analogue synaptic and neuronal elements has emerged as a potential solution to tackle the growing need for more efficient data processing, but progress toward practical systems has been stymied due to a lack of materials and devices with the appropriate attributes. Recently, solid state electrochemical ion-insertion, also known as electrochemical random access memory (ECRAM) has emerged as a promising approach to realize the needed device characteristics. ECRAM is a three terminal device that operates by tuning electronic conductance in functional materials through solid-state electrochemical redox reactions. This mechanism can be considered as a gate-controlled bulk modulation of dopants and/or phases in the channel. Early work demonstrating that ECRAM can achieve nearly ideal analogue synaptic characteristics has sparked tremendous interest in this approach. More recently, the realization that electrochemical ion insertion can be used to tune the electronic properties of many types of materials including transition metal oxides, layered two-dimensional materials, organic and coordination polymers, and that the changes in conductance can span orders of magnitude has further attracted interest in ECRAM as the basis for analogue synaptic elements for inference accelerators as well as for dynamical devices that can emulate a wide range of neuronal characteristics for implementation in analogue spiking neural networks. At its core, ECRAM shares many fundamental aspects with rechargeable batteries, where ion insertion materials are used extensively for their ability to reversibly store charge and energy. Computing applications, however, present drastically different requirements: systems will require many millions of devices, scaled down to tens of nanometers, all while achieving reliable electronic-state tuning at scaled-up rates and endurances, and with minimal energy dissipation and noise. Further, in this review, we discuss the history, basic concepts, recent progress, as well as the challenges and opportunities for different types of ECRAM, broadly grouped by their primary mobile ionic charge carrier, including Li, protons, and oxygen vacancies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Optical properties and composition of viscous organic particles found in the Southern Great Plains

Abstract. Atmospheric high-viscosity organic particles (HVOPs) were observed in samples of ambient aerosols collected in April and May 2016 in the Southern Great Plains of the United States. These particles were apportioned as either airborne soil organic particles (ASOPs) or tar balls (TBs) from biomass burning based on spetro-microscopic imaging and assessments of meteorological records of smoke and precipitation data. Regardless of their apportionment, the number fractions of HVOPs were positively correlated (R2=0.85) with increased values of absorption Ångström exponent (AAE) measured in situ for ambient aerosol at the site. Extending this correlation to 100 % HVOPs yields an AAE of 2.6, similar to previous literature reports of the class of light-absorbing organic particles known as brown carbon (BrC). One out of the three samples investigated had a significant number of ASOPs, while the other two samples contained TBs. Although there are chemical similarities between ASOPs and TBs, they can be distinguished based on composition inferred from near-edge absorption X-ray fine structure (NEXAFS) spectroscopy. ASOPs were distinguished from TBs based on their average -COOH/C=C and -COOH/COH peak ratios, with ASOPs having lower ratios. NEXAFS spectra of filtered soil organic brine particles nebulized from field samples of standing water deposited after rain were consistent with ASOPs when laboratory particles were generated by bubble bursting at the air–organic brine interface. However, particles generated by nebulizing the bulk volume of soil organic brine had a particle composition different from ASOPs. These observations are consistent with the raindrop generation mechanism responsible for ASOP emissions in the area of study. In contrast, nebulized samples carry with them higher fractions of soil inorganics dissolved in the bulk volume of soil brine, which are not aerosolized by the raindrop mechanism. Our results support the bubble bursting mechanism of particle generation during rainfall resulting in the ejection of soil organics into the atmosphere. In addition, our results show that ASOPs may only be atmospherically relevant during times when suitable emission conditions are met.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Modeling of melting and erosion in high velocity microparticle impacts

High velocity collisions can cause melting, a situation involving elasticity, plasticity, and phase transformations in tandem, which is very challenging for mechanical modeling. Here, this paper advances the use of smoothed particle hydrodynamics (SPH) simulations for this complex situation, investigating stainless steel microspheres impacting a tin substrate. After extensive validation against a range of mechanical and geometrical measurements of experimental impacts, the model provides a detailed, quantitative view of melt volume production over a broad range of impact velocities. It also enables direct decomposition of the impact energy partitioning, including analysis of the energy absorbed by plasticity in the impactor, the system kinetic energy, and the adiabatic heating of the substrate in addition to the melting itself. Since these quantities are not available experimentally, SPH modeling substantially enhances the quantitative understanding of such impacts. Finally, the model speaks to the phase characteristics of the ejecta formed during impact, confirming the experimental inference that the ejected material likely contains a significant fraction of solid material.

Constitutive modeling↗

DOE OSTI2023

Determining Stress in Metallic Conducting Layers of Microelectronics Devices Using High Resolution Electron Backscatter Diffraction and Finite Element Analysis

Delayed failure due to stress voiding is a concern with some aging microelectronics, as these voids can grow large enough to cause an open circuit. Local measurements of stress in the metallic layers are crucial to understanding and predicting this failure, but such measurements are complicated by the fact that exposing the aluminum conducting lines will relieve most of their stress. In this study, we instead mechanically thin the device substrate and measure distortions on the thinned surface using high resolution electron backscatter diffraction (HREBSD). These measurements are then related to the stresses in the metallic layers through elastic simulations. This study found that in legacy components that had no obvious voids, the stresses were comparable to the theoretical stresses at the time of manufacture (≈300 MPa). Distortion fields in the substrate were also determined around known voids, which may be directly compared to stress voiding models. In conclusion, the technique presented here for stress determination, HREBSD coupled with finite element analysis to infer subsurface stresses, is a valuable tool for assessing failure in layered microelectronics devices.

HREBSD↗

DOE OSTI2024

Lineage frequency time series reveal elevated levels of genetic drift in SARS-CoV-2 transmission in England

Genetic drift in infectious disease transmission results from randomness of transmission and host recovery or death. The strength of genetic drift for SARS-CoV-2 transmission is expected to be high due to high levels of superspreading, and this is expected to substantially impact disease epidemiology and evolution. However, we don’t yet have an understanding of how genetic drift changes over time or across locations. Furthermore, noise that results from data collection can potentially confound estimates of genetic drift. To address this challenge, we develop and validate a method to jointly infer genetic drift and measurement noise from time-series lineage frequency data. Our method is highly scalable to increasingly large genomic datasets, which overcomes a limitation in commonly used phylogenetic methods. We apply this method to over 490,000 SARS-CoV-2 genomic sequences from England collected between March 2020 and December 2021 by the COVID-19 Genomics UK (COG-UK) consortium and separately infer the strength of genetic drift for pre-B.1.177, B.1.177, Alpha, and Delta. We find that even after correcting for measurement noise, the strength of genetic drift is consistently, throughout time, higher than that expected from the observed number of COVID-19 positive individuals in England by 1 to 3 orders of magnitude, which cannot be explained by literature values of superspreading. Our estimates of genetic drift suggest low and time-varying establishment probabilities for new mutations, inform the parametrization of SARS-CoV-2 evolutionary models, and motivate future studies of the potential mechanisms for increased stochasticity in this system.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

PaleoSTeHM v1.0: a modern, scalable spatiotemporal hierarchical modeling framework for paleo-environmental data

Abstract. Geological records of past environmental change provide crucial insights into long-term climate variability, trends, non-stationarity, and nonlinear feedback mechanisms. However, reconstructing spatiotemporal fields from these records is statistically challenging due to their sparse, indirect, and noisy nature. Here, we present PaleoSTeHM, a scalable and modern framework for spatiotemporal hierarchical modeling of paleo-environmental data. This framework enables the implementation of flexible statistical models that rigorously quantify spatial and temporal variability from geological data while clearly distinguishing measurement and inferential uncertainty from process variability. We illustrate its application by reconstructing temporal and spatiotemporal paleo-sea-level changes across multiple locations. Using various modeling and analysis choices, PaleoSTeHM demonstrates the impact of different methods on inference results and computational efficiency. Our results highlight the critical role of model selection in addressing specific paleo-environmental questions, showcasing the PaleoSTeHM framework's potential to enhance the robustness and transparency of paleo-environmental reconstructions.

58 GEOSCIENCES↗

DOE OSTI2024

Diverse signatures of convergent evolution in cactus-associated yeasts

Many distantly related organisms have convergently evolved traits and lifestyles that enable them to live in similar ecological environments. However, the extent of phenotypic convergence evolving through the same or distinct genetic trajectories remains an open question. Here, we leverage a comprehensive dataset of genomic and phenotypic data from 1,049 yeast species in the subphylum Saccharomycotina (Kingdom Fungi, Phylum Ascomycota) to explore signatures of convergent evolution in cactophilic yeasts, ecological specialists associated with cacti. We inferred that the ecological association of yeasts with cacti arose independently approximately 17 times. Using a machine learning–based approach, we further found that cactophily can be predicted with 76% accuracy from both functional genomic and phenotypic data. The most informative feature for predicting cactophily was thermotolerance, which we found to be likely associated with altered evolutionary rates of genes impacting the cell envelope in several cactophilic lineages. We also identified horizontal gene transfer and duplication events of plant cell wall–degrading enzymes in distantly related cactophilic clades, suggesting that putatively adaptive traits evolved independently through disparate molecular mechanisms. Notably, we found that multiple cactophilic species and their close relatives have been reported as emerging human opportunistic pathogens, suggesting that the cactophilic lifestyle—and perhaps more generally lifestyles favoring thermotolerance—might preadapt yeasts to cause human disease. This work underscores the potential of a multifaceted approach involving high-throughput genomic and phenotypic data to shed light onto ecological adaptation and highlights how convergent evolution to wild environments could facilitate the transition to human pathogenicity.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

WHOLESCALE - Water & Hole Observations Leverage Effective Stress Calculations And Lessen Expenses (Final Technical Report 2020 - 2024)

The WHOLESCALE acronym stands for Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses. The goal of the WHOLESCALE project is to simulate the spatial distribution and temporal evolution of stress in the geothermal system at San Emidio in Nevada, United States. To reach this goal, the WHOLESCALE team has developed a methodology to incorporate and interpret data from four methods of measurement into a multi-physics model that couples thermal, hydrological, and mechanical (T H-M) processes. The WHOLESCALE team has applied this methodology at the San Emidio geothermal field, located ~100 km north of Reno, Nevada in the northwestern Basin and Range province. The WHOLESCALE team includes 30 individuals working at two universities, two national laboratories, and one industry partner. Two master-degree students and five post-doctoral researchers have gained professional experience and earned partial financial support via the WHOLESCALE project. The WHOLESCALE team has taken advantage of the perturbations created by changes in pumping operations during planned shutdowns in 2016, 2021, and 2022 to infer temporal changes in the state of stress in the geothermal system at San Emidio, Nevada, U.S. The WHOLESCALE results support the working hypothesis that increasing pore-fluid pressure reduces the effective normal stress acting across fault zones. During normal operations, pumping in deep production wells decreases fluid pressures and thus increases the effective normal stresses on faults, reducing microseismicity. During planned shutdowns, the cessation of production increases pore-fluid pressure and reduces effective normal stress. The WHOLESCALE products generated during the 4-year period between 2020 and 2024 include: three articles published in the open-access, peer-reviewed scientific literature, two master’s theses, 20 presentations or papers at scientific conferences, and 17 data sets available on public repositories. The WHOLESCALE project has been completed in two phases that included three performance periods separated by two Go/No-go Stage Gate Reviews. Tasks were classified by data type (i.e., Geologic Structure, Borehole, Geodesy, Hydrology, Seismology, and Modeling). The first phase of the project started July 31, 2020 and included ongoing project coordination (Task 1), a project kickoff (Task 2), analysis of existing data (Task 3), development of the initial stress model & deployment design (Task 4), and Go/No-go Decision Point #1 (Task 5). Phase II began with implementing the 2022 deployment (Task 6), followed by Go/No-go Decision Point #2 (Task 7) The remainder of Phase II consisted of analyzing data collected during deployment (Task 8), calibration of the stress model on all observations (Task 9), and the Final Review (August 23, 2024) & Reporting (Task 10).

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

Whole-genome sequencing distinguishes the two most common giant kelp ecomorphs

Abstract Giant kelp, Macrocystis pyrifera, exists as distinct morphological variants—or “ecomorphs”—in different populations, yet the mechanism for this variation is uncertain, and environmental drivers for either adaptive or plastic phenotypes have not been identified. The ecomorphs Macrocystis “pyrifera” and M. “integrifolia” are distributed throughout temperate waters of North and South America with almost no geographic overlap and exhibit an incongruous, non-mirrored, distribution across the equator. This study evaluates the degree of genetic divergence between M. “pyrifera” and M. “integrifolia” across 18 populations in Chile and California using whole-genome sequencing and single-nucleotide polymorphism markers. Our results based on a principal component analysis, admixture clustering by genetic similarity, and phylogenetic inference demonstrate that M. “pyrifera” and M. “integrifolia” are genetically distinguishable. Analyses reveal separation by Northern and Southern Hemispheres and between morphs within hemispheres, suggesting that the convergent “integrifolia” morphology arose separately in each hemisphere. This is the first study to use whole-genome sequencing to understand genetic divergence in giant kelp ecomorphs, identifying 83 potential genes under selection and providing novel insights about Macrocystis evolution that were not evident with previous genetic techniques. Future studies are needed to uncover the environmental forces driving local adaptation and presumed convergent evolution of these morphs.

Environmental Sciences & Ecology↗

DOE OSTI2023

Pressure Transient Analyses and Poroelastic Modeling of Hydraulic Fracture Dilation for Multiple Injections at the Devine Fracture Pilot Site

Our team has conducted electromagnetic (EM) surveys for the past six years to monitor hydraulic-fracture behavior at the Devine Fracture Pilot Site (DFPS). The sub-horizontal orientation of a shallow hydraulic fracture at the DFPS provides uniform access to the fracture area for interrogation and data collection. Ahmadian et al. (2023) suggested a possible correlation between spatiotemporal changes in the flow rate, bottomhole pressure (BHP), and the observed surface recorded electric field at the DFPS. In this paper, we present the development of poroelastic forward models and pressure transient analyses (PTAs) to support the development of a multiphysics inverse model for these EM surveys. First, we conducted PTAs of the shut-in periods after six injections out of 10 to determine the fracture closure pressure (FCP) or the overburden pressure used in a poroelastic fracture reopening model. Second, we developed a finite-element poroelastic model throughout five injection cycles to include the effect of the cumulative injected volumes due to the previous injections on current fracture dilation in the presence of highly permeable unpropped and propped zones adjacent to the cohesive layer that models fracture reopening. Fracture reopening in this poroelastic model is based on a calibrated traction-separation response using the bottomhole pressure collected in two injection campaigns in 2020 and 2022. We used the outcomes of a previous simulation study of the primary hydraulic-fracturing stimulation to define the dimension of an unpropped fracture zone ahead of the propped fracture area. The PTAs led to FCPs consistent with those obtained using the injection data collected at the DFPS in 2020. Further, these analyses showed that at later injections, the fracture closure occurred at a later time with respect to the shut-in time, inferring the effect of cumulative injected volumes in previous injections. The simulation results show that considering the propped and unpropped fracture zones improves our poroelastic model in predicting the injection-well BHP. The numerical simulation results demonstrate a significant excess pore pressure near the fracture because of the preceding formation loadings by the previous injections. The obtained fracture dilation area and fluid pressure distribution provide a basis to improve the development of a multiphysics inverse model. Furthermore, in an iteratively coupled scheme, this pressure distribution can be introduced into EM models to render a holistic view of the causative mechanisms for the surface signal anomalies.

02 PETROLEUM↗

DOE OSTI2023

Microphysical Modeling of Fault Slip and Stability Transition in Hydrothermal Conditions

Abstract Field and laboratory observations indicate that the frictional behaviors of faults depend on hydrothermal conditions. We extend the microphysical Chen‐Niemeijer‐Spiers (CNS) model to hydrothermal conditions by using the observed temperature variation of indentation hardness to infer the temperature dependence of a microphysical parameter . This parameter is assumed constant in previous versions of the CNS model. A simple spring‐slider system is used to simulate the fault system and investigate the steady‐state frictional behaviors of wet granite gouges. Our numerical results quantitatively reproduce experimental data showing the frictional‐plastic transition. The results also describe the transition from velocity‐strengthening at low temperatures (<160°C), to velocity‐weakening at intermediate temperatures (160°C–370°C), then back to velocity‐strengthening at high temperatures (>370°C). In our extended CNS model, these results suggest that the dominant shear deformation mechanism does transition from frictional granular flow to fully plastic creep with increasing temperature.

58 GEOSCIENCES↗

DOE OSTI2020

Probabilistic modeling and prediction of out-of-plane unidirectional composite lamina properties

Computational simulation provides an efficient means to predict the behavior of customized hybrid material configurations using validated, physics-based models. One limitation to this approach is the quality and quantity of available data to characterize the many constituent input properties. Therefore, a systematic approach to identify the most influential parameters on the hybrid behavior and quantify the corresponding uncertainty in predictive capabilities is required. Here, an approach using Bayesian multimodel inference and imprecise global sensitivity analysis is presented to investigate the effects of sparse constituent data on the prediction of composite material properties. The methodology allows the identification, using quantified uncertainties, of the most influential constituent material parameters for specified homogenized properties. This sensitivity analysis further enables a dimension reduction when assessing the influence of uncertainties on material properties and can be used to inform testing programs of the constituent properties that require additional testing/data collection in order to minimize uncertainty in macro-scale composite properties. The methodology is specifically demonstrated on the prediction and sensitivity analysis of out-of-plane mechanical properties of a unidirectional lamina.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Soil and soil fraction radiocarbon, C and N from the Enriched Background Isotope Study, Oak Ridge National Laboratory, Tennessee, USA

A large accidental 14CO2 release at the Oak Ridge Reservation in Tennessee, USA provided a strong label pulse into adjacent, well-studied oak forests, resulting in highly elevated Δ14C values in leaf litter (~1000‰) and roots (~260–450‰). A four-year manipulative study was conducted to provide insight into mechanisms underlying soil Organic matter dynamics, specifically to determine the relative contribution of litter versus roots to the bulk mineral soil C pool, as well as to free light, occluded light and heavy fractions. The heavy fraction was further split into fractions with densities of 1.7–2.4 g cm-3 and >2.4 g cm-3 to test the homogeneity of the mineral-associated fraction of C. Substantial concentrations of label were detected in all soil fractions within a year of the 14CO2 release, indicating rapid incorporation of newly fixed photosynthates in all fractions of soil organic C, regardless of differences in stability inferred by previous work. This rapid incorporation of label occurred only in treatments where roots were labeled, indicating that roots are the major source of inputs to mineral soil C stocks at these sites. Separation of the heavy fraction into subfractions of intermediate (1.7–2.4 g cm -3) and high (>2.4 g cm-3) density indicated that both subfractions incorporated label at similar rates, despite significant differences in degree of microbial processing.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Edge supercurrent reveals competition between condensates in a Weyl superconductor

In topological materials, the edge states are readily distinguished from the bulk states. The situation where a topological semimetal becomes superconducting so that Cooper pairs occupy both the bulk and the edge states is not well understood. In particular, we do not know if we can force their pairing symmetries to be different. Here we show that, when supercurrent is injected into the superconducting Weyl semimetal MoTe 2 from Nb contacts, the invasive s-wave pairing potential from Nb is incompatible with the intrinsic Cooper pair condensate in MoTe 2 . This incompatibility leads to strong stochasticity in the switching current and an unusual anti-hysteretic behaviour in the current–voltage loops. There is also an asymmetry in the edge oscillations where, as the magnetic field crosses zero, the phase noise switches from one with a noisy spectrum to one that is noise free. Furthermore, using the noise spectrum as a guide, we track the anomalous features to field-induced switching of the device gap function between s-wave symmetry and the unconventional symmetry intrinsic to MoTe 2 . We infer that the behaviour of the gap function along the edges is different from that in the bulk.

Superconducting properties and materials↗

DOE OSTI2026

SODAs: sparse optimization for the discovery of differential and algebraic equations

Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.

DAE↗

DOE OSTI2021

Short-range thermal magnon diffusion in magnetic garnet

Here, using the spin Seebeck effect (SSE), we study the propagation distance of thermally induced spin currents inside a magnetic insulator thin film in the short-range regime. We disambiguate spin currents driven by temperature and chemical potential gradients by comparing the SSE signal before and after adding a heat-sinking capping layer on the same device. We report that the measured spin decay behavior near the heat source is well accounted for by a diffusion model where the magnon diffusion length is in submicron range, in other words, two orders of magnitude smaller than previous estimates inferred from the long-range behavior. Our results highlight the caveat in applying a diffusive theory to describe thermally generated magnon transport, where a single decay length may not capture the behavior on all length scales.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Visible-light Phase Curves from the Second Year of the TESS Primary Mission

We carried out a systematic study of full-orbit phase curves for known transiting systems in the northern ecliptic sky that were observed during Year 2 of the TESS primary mission. We applied the same methodology for target selection, data processing, and light-curve fitting as we did in our Year 1 study. Out of the 15 transiting systems selected for analysis, seven—HAT-P-7, KELT-1, KELT-9, KELT-16, KELT-20, Kepler-13A, and WASP-12—show statistically significant secondary eclipses and day–night atmospheric brightness modulations. Small eastward dayside hot-spot offsets were measured for KELT-9b and WASP-12b. KELT-1, Kepler-13A, and WASP-12 show additional phase-curve variability attributed to the tidal distortion of the host star; the amplitudes of these signals are consistent with theoretical predictions. We combined occultation measurements from TESS and Spitzer to compute dayside brightness temperatures, TESS-band geometric albedos, Bond albedos, and phase integrals for several systems. The new albedo values solidify the previously reported trend between dayside temperature and geometric albedo for planets with 1500 K < T {sub day} < 3000 K. For Kepler-13Ab, we carried out an atmospheric retrieval of the full secondary eclipse spectrum, which revealed a noninverted temperature–pressure profile, significant H{sub 2}O and K absorption in the near-infrared, evidence for strong optical atmospheric opacity due to sodium, and a confirmation of the high geometric albedo inferred from our simpler analysis. We explore the implications of the phase integrals (ratios of Bond to geometric albedos) for understanding exoplanet clouds. We also report updated transit ephemerides for all of the systems studied in this work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Bayesian Physics Informed Spatio-Temporal Network for Streamflow Data Imputation

Reliable reconstruction of incomplete streamflow records is critical for improving hydrological forecasting, flood preparedness, and water resource management. However, large observational gaps and uncertainties in governing physical parameters limit the accuracy of traditional statistical and machinelearning imputation frameworks. To address these challenges, we develop a Bayesian Physics-Informed Spatio-Temporal Network (BPI-STNet) that jointly captures spatial and temporal dependencies while enforcing hydrologic consistency through embedded physical constraints. The framework integrates a GraphSAGE-LSTM architecture to model spatial connectivity across gauges and temporal flow dynamics, coupled with a Bayesian update mechanism to estimate uncertain parameters in a simplified water-balance framework. Unlike conventional physics-informed networks that rely on sampling-based posterior estimation, BPI-STNet derives an analytic solution to the inverse problem, allowing closed-form Bayesian updates of uncertain parameters Λ={α,β,k} using Gaussian priors and likelihoods. Applied to daily observations from the Susquehanna River Basin (1980-2022), BPI-STNet achieves substantial improvements over a purely data-driven RGNN baseline, which reduced RMSE by 23 % and MAE by 9 %, and achieving an average NSE values up to 0.96. The results demonstrate that coupling Bayesian inference with physics-informed learning yields physically consistent, uncertainty-aware reconstructions that preserve the temporal persistence and statistical distribution of observed flows. The proposed framework establishes a generalizable paradigm for data-sparse hydrologic systems where both data fidelity and physical interpretability are essential.

Krishnan Kutty Ambika, Anukesh [ORNL] (ORCID:00000↗