SEARCH · Engineering Papers

Results for “integrated modeling approach”

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 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Analysis of Multi-Output Hybrid Energy Systems Interacting with the Grid: Application of Improved Price-Taker and Price-Maker Approaches to Nuclear-Hydrogen Systems

The growing recognition of the value of hydrogen as an energy intermediate in supporting future power systems with high shares of variable renewable energy has prompted many studies to quantify the economic potential of multi-output hybrid systems, which are one type of integrated energy systems (IES). Because of the complexity of modeling multiple sectors, these studies typically use simplified modeling approaches to capture the interactions between sectors. In this study, we explore the implications of alternative modeling approaches for nuclear-hydrogen IES focusing on a power system in the Midwest United States. We combine highly resolved capacity expansion and production cost modeling tools of the power system with a detailed hydrogen system optimization tool to determine the optimal electrolyzer and storage sizing and optimal operations of the nuclear-hydrogen hybrid resource across three future study years. We compare economic and operational outcomes across a spectrum of modeling approaches, including a non-hybridized base approach; a traditional price-taker approach that does not include the impact of hydrogen production on the electricity system; a power-system-focused price-maker approach that does not account for temporal hydrogen constraints; and two improved price-taker and price-maker approaches that each address the impact of revenue-optimal levels of electricity production on the resulting power system and temporal hydrogen constraints on the overall feasible solution. Results show how a traditional price-taker approach can overestimate the economic benefits of multi-output nuclear-hydrogen IES compared to our two improved approaches that estimate both hydrogen system constraints and power system interaction. We find that hydrogen output requirements and storage size limits are key drivers to overall operations and some economic outcomes. Under our assumed constant hydrogen output requirement, storage costs, test system, and modeling approaches, our results indicate that hybridization can provide a net benefit, but results are sensitive to the treatment of hydrogen revenues and electricity prices as impacted by the power system evolution.

capacity expansion modeling↗

DOE OSTI2025

LinkML: an open data modeling framework

Background Scientific research relies on well-structured, standardized data; however, much of it is stored in formats such as free-text lab notebooks, nonstandardized spreadsheets, or data repositories. This lack of structure challenges interoperability, making data integration, validation, and reuse difficult. Findings LinkML (Linked Data Modeling Language) is an open framework that simplifies the process of authoring, validating, and sharing data. LinkML can describe a range of data structures, from flat, list-based models to complex, interrelated, and normalized models that utilize polymorphism and compound inheritance. It offers an approachable syntax that is not tied to any one technical architecture and can be integrated seamlessly with many existing frameworks. The LinkML syntax provides a standard way to describe schemas, classes, and relationships, allowing modelers to build well-defined, stable, and optionally ontology-aligned data structures. Once defined, LinkML schemas may be imported into other LinkML schemas. These key features make LinkML an accessible platform for interdisciplinary collaboration and a reliable way to define and share data semantics. Conclusions LinkML helps reduce heterogeneity, complexity, and the proliferation of single-use data models while simultaneously enabling compliance with FAIR (Findable, Accessible, Interoperable, and Reusable) data standards. LinkML has seen increasing adoption in various fields, including biology, chemistry, biomedicine, microbiome research, finance, electrical engineering, transportation, and commercial software development. In short, LinkML makes implicit models explicitly computable and allows data to be standardized at their origin. LinkML documentation and code are available at https://linkml.io/.

AI-ready data↗

DOE OSTI2025

Coarse Graining Discrete Element Method Information in Particle-in-Cell Length Scales Using a Machine Learning Approach

This report details the development of a machine learning (ML)-driven framework to coarse-grain inter-particle collision dynamics from high-fidelity Discrete Element Method (DEM) simulations to Particle-in-Cell (PIC) scales for gas-solid systems. Traditional PIC models, while computationally efficient, rely on empirical granular stress formulations that fail to capture the full complexity of collision physics, particularly the heterogeneity in particle dynamics. This study adopts a bottom-up approach, integrating insights from DEM simulations to improve the physical fidelity and interpretability of PIC-scale models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Federated Machine Learning-Based Anomaly Detection System for Synchrophasor Network Using Heterogeneous Data Sets: Preprint

Synchrophasor technology is widely deployed in the energy management system to monitor the grid health at micro level and perform necessary corrective actions in real time; however, integrated phasor devices and data aggregators are exposed to several cybersecurity threats. This paper proposes a federated ML(FML)-based ADS to detect several data integrity attacks in the synchrophasor network. The proposed approach integrates the horizontal FML technique and consists of substation-based local models and a control center-based global model. The proposed methodology includes training local models using heterogeneous data sets that include network and grid information and updating the global model through multiple iterations by sharing model gradients. Finally, the trained global model is applied to identify cyberattacks, normal operation, and physical events. To validate the proof of concept, we used synthetic data sets generated by Mississippi State University and Oak Ridge National Laboratory for training and testing the classification models using the National Renewable Energy Laboratory's high performance computing resources. Our experimental results, computed through several performance measures, reveal that the proposed approach shows consistent performance during the binary, three-class, and multiclass classifications while ensuring privacy of synchrophasor data.

anomaly detection system↗

DOE OSTI2022

A Fast Time-Stepping Strategy for Dynamical Systems Equipped with a Surrogate Model

Simulation of complex dynamical systems arising in many applications is computationally challenging due to their size and complexity. Model order reduction, machine learning, and other types of surrogate modeling techniques offer cheaper and simpler ways to describe the dynamics of these systems but are inexact and introduce additional approximation errors. In order to overcome the computational difficulties of the full complex models, on one hand, and the limitations of surrogate models, on the other, this work proposes a new accelerated time-stepping strategy that combines information from both. This approach is based on the multirate infinitesimal general-structure additive Runge--Kutta framework. The inexpensive surrogate model is integrated with a small time step to guide the solution trajectory, and the full model is treated with a large time step to occasionally correct for the surrogate model error and ensure convergence. Here, we provide a theoretical error analysis, and several numerical experiments, to show that this approach can be significantly more efficient than using only the full or only the surrogate model for the integration.

Surrogate models↗

DOE OSTI2023

Integration of Experimental Hydroprocessing and FCC Data with Process (Aspen Plus) and Refinery Optimization (Aspen PIMS) Models

NREL's Economic, Sustainability, and Market Analysis (ESMA) team develops process models in Aspen Plus to support techno-economic analysis (TEA) and life-cycle assessment (LCA) of conversion pathways from renewable and circular resources to fuels and chemicals. These tools are being applied to assess probable operational constraints or bottlenecks associated with refinery co-processing and repurposing opportunities. The team also utilizes optimizable refinery models to quantify opportunities for biofuel and bioproduct pathways through integration with existing refining infrastructure. The goals of analysis are to model the transition of refineries to renewable and circular feedstocks and quantify the costs production and CO2 abatement through utilization of existing refineries. At the core of the ESMA team's work is the process data derived from conversion experiments on hydroprocessing and fluid catalytic cracking performed at NREL. The objective of this presentation is to highlight new analysis approaches and capabilities in pathway analysis and refinery optimization modeling with specific emphasis on integration of experimental data.

BIOMASS FUELS,ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2022

Configuration Space Integration for Adsorbate Partition Functions: The Effect of Anharmonicity on the Thermophysical Properties of CO–Pt(111) and CH 3 OH–Cu(111)

A method for computing anharmonic thermophysical properties for adsorbates on metal surfaces has been extended to include libration, or frustrated rotation. Classical phase space integration is used with Monte Carlo sampling of the configuration space to obtain the partition function of CO on Pt(111) and CH 3 OH on Cu(111). A minima-preserving neural network potential energy surrogate is used within the integration routines. Direct state counting using discrete variable representation is used to benchmark the results. We find that the phase space integration approach is in excellent agreement with the direct state counting results. Comparison with standard models such as the harmonic oscillator indicates that anharmonicity contributes significantly to the thermodynamic properties of CH 3 OH on Cu(111). We find that there is also a considerable difference between the harmonic oscillator and phase space integration for CO on Pt(111), although the discrepancy can largely be attributed to the presence of multiple binding sites within the unit cell. We demonstrate that a multisite harmonic oscillator model might be sufficient for CO-Pt(111). A more thorough description of the potential energy surface, which can be achieved with phase space integration, is necessary for weakly bound adsorbates such as CH 3 OH. In conclusion, the thermophysical properties were used to calculate free energies of adsorption on the respective metals, and subsequently the equilibrium constants and Langmuir isotherms in relevant temperature ranges. The results show that the choice of model to obtain partition functions greatly affects the resulting surface coverages in kinetic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

AI‐Driven Defect Engineering for Advanced Thermoelectric Materials

Thermoelectric materials offer a promising pathway to directly convert waste heat to electricity. However, achieving high performance remains challenging due to intrinsic trade-offs between electrical conductivity, the Seebeck coefficient, and thermal conductivity, which are further complicated by the presence of defects. This review explores how artificial intelligence (AI) and machine learning (ML) are transforming thermoelectric materials design. Advanced ML approaches including deep neural networks, graph-based models, and transformer architectures, integrated with high-throughput simulations and growing databases, effectively capture structure-property relationships in a complex multiscale defect space and overcome the “curse of dimensionality”. This review discusses AI-enhanced defect engineering strategies such as composition optimization, entropy and dislocation engineering, and grain boundary design, along with emerging inverse design techniques for generating materials with targeted properties. Finally, it outlines future opportunities in novel physics mechanisms and sustainability, highlighting the critical role of AI in accelerating the discovery of thermoelectric materials.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Feasibility of DEIM for retrieving the initial field via dimensionality reduction

When parameter estimation is solved in a high-dimensional space, the dimensionality reduction strategy becomes the primary consideration for alleviating the tremendous computational cost. Here, the discrete empirical interpolation method (DEIM) is explored to retrieve the initial condition (IC) by combining the polynomial chaos (PC) based ensemble Kalman filter (i.e. PC-EnKF), where a non-intrusive PC expansion is considered as a surrogate model in place of the forward model in the prediction step of the ensemble Kalman filter, resulting in fewer forward model integrations but with a comparable accuracy as Monte Carlo-based approaches. The DEIM acts as a hyper-reduction tool to provide the low-dimensional input for the high-dimensional initial field, which can be reconstructed using the information on the sparse interpolation grid points that is adaptively obtained through PC-EnKF data assimilation method. Thus an innovative framework to reconstruct the IC is developed. The detailed procedure at each assimilation iteration includes: the determination of the spatial interpolation points, the estimation of the initial values on the interpolation locations using the optimal observations, and the reconstruction of IC in the full space. The current study uses the reconstruction field of initial conditions of the Navier-Stokes equations as an example to illustrate the efficacy of our method. The experimental results demonstrate the proposed algorithm achieves a satisfactory reconstruction for the initial field. The proposed method helps to extend the applicable area of DEIM in solving inverse problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Environmental Thermodynamic Effects on Future Changes of Mesoscale Convective Systems and Isolated Deep Convection Over the US: A Case Study

Mesoscale convective systems (MCSs) and isolated deep convection (IDC) significantly influence local weather conditions and the hydrological cycle across the United States. Their future changes under warming is critical to the society. However, the exact impacts of thermodynamic environmental changes on them are difficult to resolve using climate model simulations which integrate both dynamic and thermodynamic factors. This study explores a theoretical modeling approach to isolate environmental thermodynamic effects on their future changes through a case study. For a 50-day period during the 2020 summer, a pair of convection-permitting model simulations indicate contrasting changes in MCSs and IDC between inland and coastal regions under pseudo global warming. Driven by the thermodynamic environments of these simulations, a single-column parcel model indicates a decrease in the frequency of IDC occurrences, along with increases in duration and precipitation amount under warming, attributable to rising most unstable convective available potential energy (MUCAPE), convective inhibition (MUCIN), and precipitable water (PW). A multi-column parcel model reveals contrasting changes in the frequency and mean area of MCSs between inland and coastal regions, underscoring the increase in mean MUCIN over inland regions and increases in mean MUCAPE and PW over coastal regions. The increase in mean MCS area over inland regions is linked to the interplay between accelerated gust fronts and enhanced subsidence strength, which is further traced to the unchanged mean MUCIN but increased mean MUCAPE in large-scale environments. In conclusion, these case study results suggest a major role of thermodynamic environmental changes in controlling characteristics of MCSs and IDC under warming.

Mesoscale convective system↗

DOE OSTI2020

Visualization of the HIV-1 Env glycan shield across scales

The dense array of N-linked glycans on the HIV-1 envelope glycoprotein (Env), known as the “glycan shield,” is a key determinant of immunogenicity, yet intrinsic heterogeneity confounds typical structure–function analysis. Here, we present an integrated approach of single-particle electron cryomicroscopy (cryo-EM), computational modeling, and site-specific mass spectrometry (MS) to probe glycan shield structure and behavior at multiple levels. We found that dynamics lead to an extensive network of interglycan interactions that drive the formation of higher-order structure within the glycan shield. This structure defines diffuse boundaries between buried and exposed protein surface and creates a mapping of potentially immunogenic sites on Env. Analysis of Env expressed in different cell lines revealed how cryo-EM can detect subtle changes in glycan occupancy, composition, and dynamics that impact glycan shield structure and epitope accessibility. Importantly, this identified unforeseen changes in the glycan shield of Env obtained from expression in the same cell line used for vaccine production. Finally, by capturing the enzymatic deglycosylation of Env in a time-resolved manner, we found that highly connected glycan clusters are resistant to digestion and help stabilize the prefusion trimer, suggesting the glycan shield may function beyond immune evasion.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Boron nitride nanotube precursor formation during high-temperature synthesis: kinetic and thermodynamic modelling

In this work, we performed integrated modelling of the chemical pathways of formation for boron nitride nanotube (BNNT) precursors during high-temperature synthesis in a B/N 2 mixture. Integrated modelling includes quantum chemistry, Quantum–classical molecular dynamics, thermodynamic modelling, and kinetic approaches. We demonstrate that BN compounds are formed via the interaction of molecular nitrogen with small boron clusters, rather than through interactions with less reactive liquid boron. (This process can also be described as N 2 molecule fixation.) Liquid boron evaporates to produce these boron clusters (B m with m ≤ 5), which are subsequently converted into B m N n chains. The production of such chains is crucial to the growth of BNNTs because these chains form the building blocks of bigger and longer BN chains and rings, which are in turn the building blocks of fullborenes and BNNTs. Additionally, kinetic modelling revealed that B 4 N 4 and B 5 N 4 species in particular play a major role in the N 2 molecule fixation process. Furthermore, the formation of these species via reactions with B 4 and B 5 clusters is not adequately described under the assumption of thermodynamic equilibrium, as is demonstrated in our kinetic modelling. Thus, the accumulation of both B 4 N 4 and B 5 N 4 depends on the background gas pressure and the gas cooling rate. Long BN chains and rings, which are precursors of the fullborene and BNNT growth, form via self-assembly of components B 4 N 4 and B 5 N 4 . Our modelling results—particularly the increased densities of B 4 N 4 and B 5 N 4 species at higher gas pressures—explain the experimentally observed effect of gas pressure on the yield of high-quality BNNTs. The catalytic role of hydrogen was also studied; it is shown that HBNH molecules can be the main precursor of BNNT synthesis in the presence of hydrogen.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2023

Spin-textured neutron beams with orbital angular momentum

In this work, we present a rigorous theoretical framework underpinning the technique of spin-echo modulated small-angle neutron scattering (SEMSANS), and show how the technique can be extended in order to generate spin-textured neutron beams with orbital angular momentum (OAM) via birefringent neutron spin-polarization devices known as magnetic Wollaston prisms. Neutron OAM beams are mathematically characterized by a “cork-screw” phase singularity e iℓΦ about the propagation axis where ℓ is the OAM quantum number. To understand the precise relationship between the emergent OAM state and the variety of spin textures realized by various setups, we have developed a path-integral approach that in the interferometric limit makes a judicious use of magnetic Snell's law. We show that our proposed technique produces a complex two-dimensional pattern of spin-OAM entangled states which may be useful as a probe of quantum magnetic materials. We compare our path-integral approach to the well-known single-path Larmor precession model and present a pedagogical derivation of magnetic Snell's law of refraction for both massive and massless particles based on Maupertuis's action principle.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

Adaptive Online Model Update Algorithm for Predictive Control in Networked Systems

In this article, we introduce an adaptive on-line model update algorithm designed for predictive control applications in networked systems, particularly focusing on power distribution systems. Unlike traditional methods that depend on historical data for offline model identification, our approach utilizes real-time data for continuous model updates. This method integrates seamlessly with existing online control and optimization algorithms and provides timely updates in response to real-time changes. This methodology offers significant advantages, including a reduction in the communication network bandwidth requirements by minimizing the data exchanged at each iteration and enabling the model to adapt after disturbances. Furthermore, our algorithm is tailored for non-linear convex models, enhancing its applicability to practical scenarios. The efficacy of the proposed method is validated through a numerical study, demonstrating improved control performance using a synthetic IEEE test case.

data-driven model predictive control↗

DOE OSTI2024

Labeling sequential data from noisy annotations

Crowdsourcing algorithms often work under the assumption that the data samples are independent. Recent work has shown that data dependence, such as temporal correlations in sequential data, can be leveraged to improve the label quality. Existing methods that exploit this special structure rely on third-order statistics of the annotator outputs to ensure the identifiability of key latent parameters, which are costly to acquire. This work proposes an approach for integrating crowdsourced annotations under the Dawid-Skene/Hidden Markov Model (DS-HMM) for sequential data based on second-order statistics, which naturally enjoys a lower sample complexity. An effective algorithm is proposed to tackle the challenging optimization problem associated with the proposed estimator. Numerical experiments showcase the effectiveness of the data labeling paradigm.

Marrinan, Timothy P.↗

DOE OSTI2020

Preliminary description of a new creep-fatigue design method that reduces over conservatism and simplifies the high temperature design process

The report provides the initial description of a new creep-fatigue design method for structural components in high temperature nuclear service. The new method is based on an integrated elastic-perfectly plastic (EPP) analysis and Simplified Model Test (SMT) approach that reduces over conservatism, improves the treatment of elastic follow up, and simplifies the design procedure, when compared with the current creep-fatigue design methods in ASME Boiler and Pressure Vessel Code. Developing the design charts for the EPP-SMT design method requires extrapolating SMT test data as a function of hold time and follow up factor. The report develops the preliminary design charts for Alloy 617 at temperatures between 800°C and 950°C by combining two extrapolation approaches developed in a previous work. The report also presents a comparative analysis between the EPP-SMT design method and the current ASME creep-fatigue design methods by evaluating design life of two sample geometries under different loading conditions. Results from the comparative analysis verify the EPP-SMT design charts but suggest the requirement of additional test data in the low strain range regime for improving the extrapolation procedure that will further reduce the over conservatism in the creep-fatigue damage evaluation. The report also concludes that the EPP-SMT design procedure can account for effect of primary load on creep-fatigue life by using a fixed, bounding value of follow up in constructing the design charts. The conclusions to this report describe the future work required to complete this new design method so it can be codified through a nuclear Code Case.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Interactions between molecular-scale biogeochemical processes and hyporheic exchange for understanding coupled Fe-S-C cycling in iron-rich riparian wetlands (Final Technical Report)

Riparian wetlands are dynamic interfaces that exert strong control over water quality, contaminant mobility, and greenhouse gas emissions. These systems are characterized by hyporheic exchange between oxic surface water and anoxic groundwater, which generates steep redox gradients and promotes spatially and temporally variable microbial activity. Despite growing recognition of tightly coupled iron (Fe), sulfur (S), and carbon (C) cycling in these environments, the mechanisms governing these interactions—particularly under low-sulfate freshwater conditions—remain poorly constrained. This project developed a mechanistic understanding of how hydrologic variability and microbial processes interact to control Fe–S–C cycling in iron-rich riparian wetlands. Using a multi-scale and multi-method approach integrating field observations, geochemical and spectroscopic analyses, metagenomics, and reactive transport modeling, we demonstrate that “cryptic” sulfur cycling—rapid sulfur transformations involving intermediate-valence species—plays a dominant and previously underrecognized role in freshwater wetlands. These processes persist despite low sulfate concentrations and significantly influence iron reduction, carbon mineralization, and methane dynamics. The results show that cryptic sulfur cycling enhances dissolved Fe 2+ production, regulates methane concentrations, and is strongly controlled by climate-driven hydrologic fluxes. By linking hydroclimate, subsurface flow, and biogeochemical reactions, this work provides a predictive framework for understanding how wetland systems respond to environmental change, with direct implications for water quality and carbon cycling.

54 ENVIRONMENTAL SCIENCES↗