Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Content Model”

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

3013 Baseline Inner Container Report: Description, Imaging, and Analysis

This report presents consolidated findings from three microscopy studies of baseline Bagless Transfer Containers (BTCs), which are the inner containers used in DOE Standard 3013-compliant plutonium storage packages. The studies were conducted by Los Alamos National Laboratory (LANL), Savannah River National Laboratory (SRNL), and DNV, a subcontractor specializing in corrosion science. The baseline BTCs examined had not been exposed to corrosive environments, serving as control samples for comparison with destructively examined (DE) containers. Additionally, the Background section provides an overview of the BTC design, including manufacturing methods, material properties, and container configurations that may influence corrosion behavior. The primary objective of the LANL and SRNL studies was to establish reference characteristics of BTCs to distinguish manufacturing artifacts from corrosion-induced features observed in DE containers. Both laboratories used advanced microscopy techniques, including Wide Area Microscopy (WAMS), Laser Confocal Microscopy (LCM), Scanning Electron Microscopy (SEM), and Focused Ion Beam (FIB), to reveal that most surface anomalies were shallow and attributable to flow-forming processes, not corrosion. Subsurface impurities, identified as non-metallic inclusions, were also observed. DNV’s investigation focused on characterizing the mechanical and chemical properties of the BTC material to support crack growth modeling. Their findings indicate higher hardness and strength but reduced ductility—consistent with flow-formed 304L stainless steel and a higher-than-expected martensite content, in contrast to manufacturer claims. Ongoing and future studies, including Electron Backscatter Diffraction (EBSD), aim to further understand the microstructural factors influencing crack initiation and propagation. The combined studies provide a crucial baseline for distinguishing corrosion effects in DE containers, contributing to the safe and reliable long-term storage of plutonium-bearing materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Synthetic Streamflow Datasets to Support Emulation of Water Allocations via LSTM

This archive is the data companion to the bonney_et-al_2026_erc metarepo which generates synthetic data, trains an LSTM model, and generates performance metrics on the trained model. While the generation of the synthetic data is fully reprodicible, it is a computationally expensive process. This data archive contains the synthetic datasets needed for training and testing an LSTM model and reproduction of figures and tables. In addition, supplemenatary data products generating and visualizing results is also included, such as geospatial data for the basin. Contents There are two high level directories: `WRAP_archive/` and `repo_data/`. The `WRAP_archive` directory contains compressed intermediate dataproducts from the dataset generation workflow (marked as "I_Dataset_Generation" in the metarepo). These data products are not required by any scripts in the metarepo, but they are archived as they are expensive to generate and may have useful information for other analyses. The `repo_data` directory contains the necessary data for reproducing the workflow in the metarepo and should be decompressed and moved into the top level of the metarepo. Additional details are provided in README.md.

drought↗
DOE OSTI2022

DNS of Hydrodynamic Instabilities of Laminar H2/O2/N2 Flames at Elevated Pressure

The atypical combustion properties of hydrogen have long provided an interest for hydrogen-enriched combustion, emphasized in recent years by increasingly stringent emission regulations. Significant work has been focused on the small hydrocarbon(C1-C3)/hydrogen blends, but much less is known about the effect of hydrogen on heavier hydrocarbon combustion, of interest for the aviation and automotive industry. In this work, we perform direct numerical simulations (DNS) of spherically expanding laminar n-dodecane/H2 flames in a constant-volume vessel at elevated pressure and a range of hydrogen seeding levels. A low Mach number model is used to represent the gas and flame dynamics allowing for a temporally varying, spatially homogeneous pressure field, while also enabling numerical evolution of the system numerically at time step sizes governed by the advective CFL limit rather than acoustic processes. We use an adaptive mesh refinement (AMR) approach to tackle the large separation of scales associated with high pressure premixed flame surfaces propagating in a laboratory-scale closed domain. We incorporate detailed chemistry and transport models for the reacting flow and focus the analysis on the effects of the pressure rise on the flame propagation characteristics and morphology as the hydrogen content is increased. The presence of fast diffusive H2 and the pressure rise results in dramatic thinning of the flame and the potential onset of hydrodyanmic and theromdiffusive instabilities. We assess the potential role of these instabilities on the ability to accurately measure laminar burning speed experimentally based on the mean flame surface propagation speed.

high pressure↗
DOE OSTI2023

Hydrostatic and Non‐Hydrostatic Baroclinic Instability in the Dynamical Core of the DOE Global Climate Model

Abstract The dynamical core of the Department of Energy global climate model is used to understand the role of non‐hydrostatic dynamics in the simulation of dry and moist baroclinic waves. To reduce computational cost, the Diabatic Acceleration and REscaling approach is adopted. Scale analysis and numerical simulations suggest that the model solution is not distorted by the change of spherical metric terms due to the change of Earth radius, neither are the inter‐scale interactions strongly altered as indicated by the analysis of spectral flux. Compared with hydrostatic simulations, the onset of baroclinic instability is delayed under non‐hydrostatic dynamics as the associated weaker vertical motions tend to increase the critical wavelength and narrow the range of waves that can be baroclinically unstable. During the development of baroclinic waves, non‐hydrostatic dynamics tends to induce vertical motions in the upper troposphere, accelerating the eastward propagation of upper‐level ridges through their impact on local vorticity tendency. These processes reduce the westward tilt of the vertical ridge axes and suppress the conversion of mean flow available potential energy to eddy kinetic energy, leading to weaker baroclinic eddies than in hydrostatic simulations. The contrasts between hydrostatic and non‐hydrostatic settings hold in supplemental experiments that vary the background flow Rossby number, the amount of water vapor content and vertical resolution. We also find that mesoscale and smaller‐scale activities can be considerably under‐represented when the vertical resolution is limited to that typically used in global climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A leaf-level spectral library to support high-throughput plant phenotyping: predictive accuracy and model transfer

Abstract Leaf-level hyperspectral reflectance has become an effective tool for high-throughput phenotyping of plant leaf traits due to its rapid, low-cost, multi-sensing, and non-destructive nature. However, collecting samples for model calibration can still be expensive, and models show poor transferability among different datasets. This study had three specific objectives: first, to assemble a large library of leaf hyperspectral data (n=2460) from maize and sorghum; second, to evaluate two machine-learning approaches to estimate nine leaf properties (chlorophyll, thickness, water content, nitrogen, phosphorus, potassium, calcium, magnesium, and sulfur); and third, to investigate the usefulness of this spectral library for predicting external datasets (n=445) including soybean and camelina using extra-weighted spiking. Internal cross-validation showed satisfactory performance of the spectral library to estimate all nine traits (mean R2=0.688), with partial least-squares regression outperforming deep neural network models. Models calibrated solely using the spectral library showed degraded performance on external datasets (mean R2=0.159 for camelina, 0.337 for soybean). Models improved significantly when a small portion of external samples (n=20) was added to the library via extra-weighted spiking (mean R2=0.574 for camelina, 0.536 for soybean). The leaf-level spectral library greatly benefits plant physiological and biochemical phenotyping, whilst extra-weight spiking improves model transferability and extends its utility.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Disruption thermal load mitigation with shattered pellet injection on the Joint European Torus (JET)

Disruption mitigation remains a critical, unresolved challenge for ITER. To aid in addressing this challenge, a shattered pellet injection (SPI) system was installed on JET and experiments conducted at a range of thermal energy fractions and stored energies in excess of 7 MJ. The primary goals of these experiments were to investigate the efficacy of the SPI on JET and the ability of the plasma to assimilate multiple pellets. Single pellet injections produced a saturation in total radiated energy (Wrad) with increasing injected neon content, suggesting total radiation of stored thermal energy. Further increases in injected neon quantities resulted in reduced cooling times and current quench (CQ) durations, indicating higher impurity assimilation. No significant variation in CQ duration or W rad was observed when varying the deuterium content at fixed neon quantities. Additionally, higher assimilation, inferred by shorter CQ durations, was measured when a mechanical punch was used to launch the pellets and this was attributed to a lower pellet velocity leading to higher solid content in the pellet plume and larger fragments penetrating deeper into the plasma. Radiation asymmetries averaged over the cooling time were inferred from Emis3D and ranged from 1.6 to 1.9. Asymmetries averaged over the entire disruption sequence were found to increase at higher thermal energy fractions. The radiated energy fractions decreased with increasing thermal energy fractions but this trend was eliminated when toroidal asymmetries were accounted for with Emis3D. Pure deuterium pellets were able to produce cooling times of up to 75 ms with a gradual loss in thermal stored energy of up to 80%. Experiments with multiple pellet injection indicated W rad can be increased through pellet superposition and density can be increased with an additional D2 injection without a reduction in W rad . KPRAD modelling accurately reproduced the cooling times and the CQ duration at high thermal energies. Assimilation estimates from KPRAD indicated CQ rates scale strongly whilst W rad scales weakly and saturates with assimilated neon content. Comparable W rad can be achieved with lower assimilated neon quantities as longer cooling times are attained. Thus reduced neon content can be preferential in a thermal load mitigation scheme as it may reduce radiation asymmetries and prevent flash melting.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The importance of Aitken mode aerosol particles for cloud sustenance in the summertime high Arctic – a simulation study supported by observational data

The potential importance of Aitken mode particles (diameters ~ 25–80 nm) for stratiform mixed-phase clouds in the summertime high Arctic (>80° N) has been investigated using two large-eddy simulation models. We find that, in both models, Aitken mode particles significantly affect the simulated microphysical and radiative properties of the cloud and can help sustain the cloud when accumulation mode concentrations are low (< 10–20 cm –3 ), even when the particles have low hygroscopicity (hygroscopicity parameter – κ=0.1). However, the influence of the Aitken mode decreases if the overall liquid water content of the cloud is low, either due to a higher ice fraction or due to low radiative cooling rates. An analysis of the simulated supersaturation (ss) statistics shows that the ss frequently reaches 0.5 % and sometimes even exceeds 1 %, which confirms that Aitken mode particles can be activated. The modelling results are in qualitative agreement with observations of the Hoppel minimum obtained from four different expeditions in the high Arctic. Our findings highlight the importance of better understanding Aitken mode particle formation, chemical properties and emissions, particularly in clean environments such as the high Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

AATTENUATION—The Atmospheric Attenuation Model for CSP Tower Plants: A Look-Up Table for Operational Implementation

Attenuation of solar radiation between the receiver and the heliostat field in concentrated solar power (CSP) tower plants can reduce the overall system performance significantly. The attenuation varies strongly with time and the average attenuation at different sites might also vary strongly from each other. If no site specific attenuation data is available, the optimal plant design cannot be determined and rough estimations of the attenuation effect are required leading to high uncertainties of yield analysis calculations. The attenuation is caused mainly by water vapor content and aerosol particles in the lower atmospheric layer above ground. Although several on-site measurement systems have been developed during recent years, attenuation data sets are usually not available to be included during the plant project development. An Atmospheric Attenuation (AATTENUATION) model to derive the atmospheric transmittance between a heliostat and receiver on the basis of common direct normal irradiance (DNI), temperature, relative humidity, and barometric pressure measurements was developed and validated by the authors earlier. The model allows the accurate estimation of attenuation for sites with low attenuation and gives an estimation of the attenuation for less clear sites. However, the site-dependent coefficients of the AATTENUATION model had to be developed individually for each site of interest, which required time-consuming radiative transfer simulations, considering the exact location and altitude, as well as the pre-dominant aerosol type at the location. This strongly limited the application of the model despite its typically available input data. In this manuscript, a look-up table (LUT) is presented which enables the application of the AATTENUATION model at the site of interest without the necessity to perform the according complex radiative transfer calculations for each site individually. This enables the application of the AATTENUATION model for virtually all resource assessments for tower plants and in an operational mode in real time within plant monitoring systems around the world. The LUT also facilitates the generation of solar attenuation maps on the basis of long-term meteorological data sets which can be considered during resource assessment for CSP tower plant projects. The LUTs are provided together with this manuscript as supplementary files. The LUT for the AATTENUATION model was developed for a solar zenith angle (SZA) grid of 1°, an altitude grid of 100 m, 7 different standard aerosol types and the standard AFGL atmospheres for mid-latitudes and the tropics. The LUT was tested against the original version of the AATTENUATION model at 4 sites in Morocco and Spain, and it was found that the additional uncertainty introduced by the application of the LUT is negligible. With the information of latitude, longitude, altitude above mean sea level, DNI, relative humidity (RH), ambient temperature (Tair), and barometric pressure (bp), the attenuation can be now derived easily for each site of interest.

14 SOLAR ENERGY↗
DOE OSTI2020

Data from Unique Contributions of Chlorophyll and Nitrogen to Predict Crop Photosynthetic Capacity from Leaf Spectroscopy

The photosynthetic capacity or the CO2-saturated photosynthetic rate (Vmax), chlorophyll, and nitrogen are closely linked leaf traits that determine C4 crop photosynthesis and yield. Accurate, timely, rapid, and non-destructive approaches to predict leaf photosynthetic traits from hyperspectral reflectance are urgently needed for high-throughput crop monitoring to ensure food and bioenergy security. Therefore, this study thoroughly evaluated the state-of-the-art physically based radiative transfer models (RTMs), data-driven partial least squares regression (PLSR), and generalized PLSR (gPLSR) models to estimate leaf traits from leaf-clip hyperspectral reflectance, which was collected from maize (Zea mays L.) bioenergy plots with diverse genotypes, growth stages, treatments with nitrogen fertilizers, and ozone stresses in three growing seasons. The results show that leaf RTMs considering bidirectional effects can give accurate estimates of chlorophyll content (Pearson correlation r=0.95), while gPLSR enabled retrieval of leaf nitrogen concentration (r=0.85). Using PLSR with field measurements for training, the cross-validation indicates that Vmax can be well predicted from spectra (r=0.81). The integration of chlorophyll content (strongly related to visible spectra) and nitrogen concentration (linked to shortwave infrared signals) can provide better predictions of Vmax (r=0.71) than only using either chlorophyll or nitrogen individually. This study highlights that leaf chlorophyll content and nitrogen concentration have key and unique contributions to Vmax prediction.

Biomass Analytics↗
DOE OSTI2022

Electrochemical Modeling and Experimental Verification of Lithiation Gradients in Oxide Cathodes of Lithium-Ion Cells

Lithiated nickel-cobalt-manganese oxides, such as NCM523, are used in the positive electrode (cathode) of Li-ion cells. Using operando X-ray diffraction profilometry, lithiation gradients in the cathode matrix can be observed and quantified by expansion into Legendre polynomials with time-dependent weights. These weights (referred to as gradients) increase in magnitude when electric current flows through the cell, decrease during potentiostatic hold and finally relax to zero when the current is interrupted during open circuit rest. Both physics-based electrochemical models and operando X-ray experiments suggest that the time constants for gradient growth and abatement are primarily determined by ionic diffusion in the oxide particles, which in turn depends on their lithium content. In contrast, the magnitude of gradients depends mainly on the applied current. The X-ray profilometry provides a way of directly probing the formation and disappearance of Li gradients across the cathode during fast cycling, which can help to diagnose the effects of material degradation in the cells.

25 ENERGY STORAGE↗
DOE OSTI2022

Translation and rotation equivariant normalizing flow (TRENF) for optimal cosmological analysis

ABSTRACT Our Universe is homogeneous and isotropic, and its perturbations obey translation and rotation symmetry. In this work, we develop translation and rotation equivariant normalizing flow (TRENF), a generative normalizing flow (NF) model which explicitly incorporates these symmetries, defining the data likelihood via a sequence of Fourier space-based convolutions and pixel-wise non-linear transforms. TRENF gives direct access to the high dimensional data likelihood p(x|y) as a function of the labels y, such as cosmological parameters. In contrast to traditional analyses based on summary statistics, the NF approach has no loss of information since it preserves the full dimensionality of the data. On Gaussian random fields, the TRENF likelihood agrees well with the analytical expression and saturates the Fisher information content in the labels y. On non-linear cosmological overdensity fields from N-body simulations, TRENF leads to significant improvements in constraining power over the standard power spectrum summary statistic. TRENF is also a generative model of the data, and we show that TRENF samples agree well with the N-body simulations it trained on, and that the inverse mapping of the data agrees well with a Gaussian white noise both visually and on various summary statistics: when this is perfectly achieved the resulting p(x|y) likelihood analysis becomes optimal. Finally, we develop a generalization of this model that can handle effects that break the symmetry of the data, such as the survey mask, which enables likelihood analysis on data without periodic boundaries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Physics-Informed Machine-Learning Prediction of Curie Temperatures and Its Promise for Guiding the Discovery of Functional Magnetic Materials

High-performance permanent magnets with a high Curie temperature, containing less critical materials, are integral to zero-carbon energy solutions. We built a machine-learning model trained over available experimentally measured Curie temperature values to predict the T C of multicomponent magnetic materials. We chose two compositions from a pseudo-binary (Zr 1–x Ce x )Fe 2 system, namely, (Zr 0.16 Ce 0.84 )Fe 2 and (Zr 0.94 Ce 0.06 )Fe 2 , to experimentally validate the ability of our model to predict the Curie temperature of novel compounds. We also provided a detailed discussion on the correlation of the Curie temperature with the de Gennes scaling factor in rare-earth intermetallic compounds and its breakdown below a certain rare-earth content. The electronic structure calculations (density of states and Fermi surface) were performed using the density functional theory on selected compounds (Zr 0.16 Ce 0.84 )Fe 2 and (Zr 0.94 Ce 0.06 )Fe 2 to understand the electronic origin of a strong magnetic exchange. We found that the change in the electronic density of states and electron/hole fillings at the Fermi level directly correlate with the Curie temperature. Notably, our model was able to capture these key electronic structure trends, which show that physics-informed machine learning can play a crucial role in designing new high-performance magnets with improved properties for environmentally sustainable applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Implications of the correlation between bulge-to-total baryonic mass ratio and the number of satellites for SAGA galaxies

We searched for correlations between the number of satellites and fundamental galactic properties for the Milky Way-like host galaxies in order to better understand their diverse satellite populations. We specifically aim to understand why galaxies that are very similar in stellar mass content, star formation rate, and local environment have very different numbers of satellites. Deep and extensive spectroscopic observations are needed to characterize the complete satellite luminosity function beyond the Local Group. One such endeavor is an ongoing Satellites of Galactic Analogs (SAGA) spectroscopic survey that has completed spectroscopic observations of 36 Milky Way-like galaxies within their virial radii down to the luminosity of Leo I dwarf galaxy. We correlated the number of satellites of SAGA galaxies with several fundamental properties of their hosts – including total specific angular momentum, which is considered to be well preserved throughout galaxy lifetime – in an attempt to identify the main driver of their diverse satellite populations. We aim to reveal some intrinsic galactic property decisive in making more or less satellites irrespective of baryonic mass or the environment in which galaxies reside. We modeled Spitzer Heritage Archive images of SAGA host galaxies at 3.6 and 4.5 microns with GALFIT code to obtain their stellar masses. We also searched the Extragalactic Database for information on their gas content and rotation velocities. Empirical correlations, like the baryonic Tully–Fisher relation and the stellar mass–size relation were used to exclude outliers. All the available galactic properties from the literature along with measured stellar masses were correlated with the number of satellites and no significant correlation was found. However, when we considered the “expected” number of satellites based on the correlation between the baryonic bulge-to-total ratio and the number of satellites confirmed for several nearby galaxies then strong correlations emerge between this number and (1) the mass of the bulge, and (2) the total specific angular momentum. The first correlation is positive, implying that galaxies with more massive bulges have more satellites, as already confirmed. Furthermore, the second correlation with the angular momentum is negative, meaning that, the smaller the angular momentum, the greater the number of expected satellites. This would imply that either satellites cannot form if galaxy angular momentum is too high, or that satellites form inside-out, so that angular momentum is being transferred to the outer parts of the galaxies. However, deeper spectroscopic observations are needed to confirm these findings, because they rely on the expected rather than detected number of satellites. There was a luminosity limit to the SAGA survey equivalent to the luminosity of Leo I dwarf satellite of the Milky Way galaxy (the SAGA limit). In particular, correlations found in this work are very susceptible to the total number of satellites of the NGC 4158 galaxy. This galaxy is predicted to have many more satellites than detected up to the SAGA limit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Surface soil temperature and water content from warming experiment located at the Rocky Mountain Biological Laboratory, Colorado, 2015 to 2019

This data package consists of soil temperature and soil water content sensor data from the warming experiment near the Rocky Mountain Biological Laboratory. Data range is from summer 2015 to summer 2019, when the warming manipulation was terminated. The location of the warming experiment is in the upper montane life zone and consists of meadow vegetation with a mixture of sagebrush (Artemisia tridentata), forb species (Erigeron speciosus, Helianthella quinquenervis, Lathyrus lanszwertii), and graminoid species (Festuca thurberi, Carex spps). Early studies showed that the heating treatment led to earlier snowmelt as well as warmer and drier surface soils, and changes in the functional composition of aboveground biomass favoring sagebrush over forb species (Harte and Shaw, 1995). These changes were mainly observed in the upper, naturally drier zone for the experimental plots. More recent observations of aboveground biomass, soil carbon, and microbial composition require concurrent information about the abiotic environment for interpretation and extrapolation to process-based models. This data package provides the plot- and zone-level abiotic context for the final years of the experiment and consists of csv files with hourly soil temperature (WM_ST.csv) and soil water content (WM_SWC.csv).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Strangeness enhancement at its extremes: multiple (multi-)strange hadron production in pp collisions at \(\sqrt{s}=5.02\) TeV

The probability to observe a specific number of strange and multi-strange hadrons (nS), denoted as P(nS), is measured by ALICE at midrapidity (|y| < 0.5) in $$\sqrt{s}=5.02$$ TeV proton-proton (pp) collisions, dividing events into several multiplicity-density classes. Exploiting, for the first time, a technique based on counting the number of strange-particle candidates event-by-event, this measurement allows one to extend the study of strangeness production beyond the mean of the distribution. This constitutes a new test bench for production mechanisms, probing events with a large imbalance between strange and non-strange content. The analysis of a large-statistics data sample makes it possible to extract P(nS) up to a maximum nS of 7 for $${\text{K}}_{\text{S}}^{0}$$, 5 for Λ and $$\overline{\Lambda }$$, 4 for Ξ− and $${\overline{\Xi } }^{+}$$, and 2 for Ω− and $${\overline{\Omega } }^{+}$$. From this, the probability of producing strange hadron multiplets per event is calculated, thereby enabling the extension of the study of strangeness enhancement to extreme situations where several strange quarks hadronize in a single event at midrapidity. Moreover, comparing hadron combinations with different u and d quark compositions and equal overall s quark content, the contribution to the enhancement pattern coming from non-strangeness related mechanisms is isolated. The results are compared with state-of-the-art phenomenological models implemented in commonly used Monte Carlo event generators, including PYTHIA 8 Monash 2013, PYTHIA 8 with QCD-based Color Reconnection and Rope Hadronization (QCD-CR + Ropes), and EPOS LHC, which incorporates both partonic interactions and hydrodynamic evolution. These comparisons show that the new approach dramatically enhances the sensitivity to the different underlying physics mechanisms modeled by each generator.

Abualrob, I J↗
DOE OSTI2022

Receptance coupling substructure analysis and chatter frequency-informed machine learning for milling stability

This paper describes a milling stability identification approach that simultaneously considers: physics-based models for the tool tip frequency response functions and stability predictions; the binary result from a milling test (automatically labeled as stable or unstable based on frequency content); chatter frequency when an unstable result is obtained; and user risk tolerance. The algorithm applies probabilistic Bayesian machine learning with adaptive, parallelized Markov Chain Monte Carlo sampling to update the probability of stability with each milling test. Furthermore, the result is a robust solution for rapid convergence to optimized milling parameters for maximum metal removal rate using all available information.

42 ENGINEERING↗
DOE OSTI2020

Carbon Limitation Leads to Thermodynamic Regulation of Aerobic Metabolism

Organic matter (OM) metabolism in freshwater ecosystems is a critical source of uncertainty in global biogeochemical cycles, yet the processes regulating aerobic respiration in aquatic environments remain poorly understood. Aerobic respiration is typically predicted through kinetic controls such as organic carbon and oxygen concentrations, while assuming thermodynamic regulation has no influence. The role of OM thermodynamics has been mainly explored under anaerobic conditions; however, recent observations have suggested that thermodynamic properties of OM may influence aerobic respiration rates. Here, we explicitly test OM thermodynamics as a key regulator of aerobic respiration. We pair controlled microcosm experiments with ultrahigh-resolution OM characterization to show that aerobic respiration increases with OM thermodynamic favorability under carbon limitation. We also demonstrate a shift in the regulation of aerobic respiration from OM thermodynamics to nitrogen content when carbon is in excess, highlighting a central role for OM thermodynamics particularly in carbon-limited ecosystems. This illuminates a structural gap in aquatic biogeochemical models by demonstrating that both kinetic and thermodynamic constraints can modulate aerobic respiration. Our work therefore proposes a new conceptual model in which thermodynamic and nutrient limitations dually control aerobic respiration.

54 ENVIRONMENTAL SCIENCES↗