Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “carbon dynamics”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

CO 2 adsorption mechanisms at the ZIF-8 interface in a Type 3 porous liquid

Porous liquids (PLs) are an attractive material for gas separation and carbon sequestration due to their permanent internal porosity and high adsorption capacity. PLs that contain zeolitic imidazole frameworks (ZIFs), such as ZIF-8, form PLs through exclusion of aqueous solvents from the framework pore due to its hydrophobicity. The gas adsorption sites in ZIF-8 based PLs are historically unknown; gas molecules could be captured in the ZIF-8 pore or adsorb at the ZIF-8 interface. To address this question, ab initio molecular dynamics was used to predict CO 2 binding sites in a PL composed of a ZIF-8 particle solvated in a water, ethylene glycol, and 2-methylimidazole solvent system. Further, the results show that CO 2 energetically prefers to reside inside the ZIF-8 pore aperture due to strong van der Waals interactions with the terminal imidazoles. However, the CO 2 binding site can be blocked by larger solvent molecules that have greater adsorption interactions. CO 2 molecules were unable to diffuse into the ZIF-8 pore, with CO 2 adsorption occurring due to binding with the ZIF-8 surface. Therefore, future design of ZIF-based PLs for enhanced CO 2 adsorption should be based on the strength of gas binding at the solvated particle surface.

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

Topography and canopy cover influence soil organic carbon composition and distribution across a forested hillslope in the discontinuous permafrost zone

Abstract Topography and canopy cover influence ground temperature in warming permafrost landscapes, yet soil temperature heterogeneity introduced by mesotopographic slope positions, microtopographic differences in vegetation cover, and the subsequent impact of contrasting temperature conditions on soil organic carbon (SOC) dynamics are understudied. Buffering of permafrost‐affected soils against warming air temperatures in boreal forests can reflect surface soil characteristics (e.g., thickness of organic material) as well as the degree and type of canopy cover (e.g., open cover vs. closed cover). Both landscape and soil properties interact to determine meso‐ and microscale heterogeneity of ground warming. We sampled a hillslope catena transect in a discontinuous permafrost zone near Fairbanks, Alaska, to test the small‐scale (1 to 3 m) impacts of slope position and cover type on soil organic matter composition. Mineral active layer samples were collected from backslope, low backslope, and footslope positions at depths spanning 19 to 60 cm. We examined soil mineralogical composition, soil moisture, total carbon and nitrogen content, and organic mat thickness in conjunction with an assessment of SOC composition using Fourier‐transform ion cyclotron resonance mass spectrometry (FT‐ICR‐MS). Soils in the footslope position had a higher relative contribution of lignin‐like compounds, whereas backslope soils had more aliphatic and condensed aromatic compounds as determined using FT‐ICR‐MS. The effect of open versus closed tree canopy cover varied with the slope position. On the backslope, we found higher oxidation of molecules under open cover than closed cover, indicating an effect of warmer soil temperature on decomposition. Little to no effect of the canopy was observed in soils at the footslope position, which we attributed, in part, to the strong impact of soil moisture content in SOC dynamics in the water‐gathering footslope position. The thin organic mat under open cover on the backslope position may have contributed to differences in soil temperature and thus SOC oxidation under open and closed canopies. Here, the thinner organic mat did not appear to buffer the underlying soil against warm season air temperatures and thus increased SOC decomposition as indicated by the higher oxidation of SOC molecules and a lower contribution of simple molecules under open cover than the closed canopy sites. Our findings suggest that the role of canopy cover in SOC dynamics varies as a function of landscape position and soil properties, namely, organic mat thickness and soil moisture. Condition‐specific heterogeneity of SOC composition under open and closed canopy cover highlights the protective effect of canopy cover for soils on backslope positions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2019

The Arctic-Boreal vulnerability experiment model benchmarking system

NASA's Arctic-Boreal Vulnerability Experiment (ABoVE) integrates field and airborne data into modeling and synthesis activities for understanding Arctic and Boreal ecosystem dynamics. The ABoVE Benchmarking System (ABS) is an operational software package to evaluate terrestrial biosphere models against key indicators of Arctic and Boreal ecosystem dynamics, i.e.: carbon biogeochemistry, vegetation, permafrost, hydrology, and disturbance. The ABS utilizes satellite remote sensing data, airborne data, and field data from ABoVE as well as collaborating research networks in the region, e.g.: the Permafrost Carbon Network, the International Soil Carbon Network, the Northern Circumpolar Soil Carbon Database, AmeriFlux sites, the Moderate Resolution Imaging Spectroradiometer, the Orbiting Carbon Observatory 2, and the Soil Moisture Active Passive mission. The ABS is designed to be interactive for researchers interested in having their models accurately represent observations of key Arctic indicators: a user submits model results to the system, the system evaluates the model results against a set of Arctic-Boreal benchmarks outlined in the ABoVE Concise Experiment Plan, and the user then receives a quantitative scoring of model strengths and deficiencies through a web interface. This interactivity allows model developers to iteratively improve their model for the Arctic-Boreal Region by evaluating results from successive model versions. We show here, for illustration, the improvement of the Lund–Potsdam–Jena-Wald Schnee und Landschaft (LPJwsl) version model through the ABoVE ABS as a new permafrost module is coupled to the existing model framework. The ABS will continue to incorporate new benchmarks that address indicators of Arctic-Boreal ecosystem dynamics as they become available.

ABoVE↗
DOE OSTI2023

Topography and canopy cover influence soil organic carbon composition and distribution across a forested hillslope in the discontinuous permafrost zone

This dataset contains data used for the paper "Topography and canopy cover influence soil organic carbon composition and distribution across a forested hillslope in the discontinuous permafrost zone". The Related References field will be updated with a full citation when available. Topography and canopy cover influence ground temperature in warming permafrost landscapes, yet soil temperature heterogeneity introduced by meso-topographic slope positions, microtopographic differences in vegetation cover, and the subsequent impact of contrasting temperature conditions to soil organic carbon (SOC) dynamics are understudied. Buffering of permafrost-affected soils against warming air temperatures in boreal forests can reflect surface soil characteristics (e.g., thickness of organic material) as well as the degree and type of canopy cover (e.g., open cover vs closed cover). Both landscape and soil properties interact to determine meso- and micro-scale heterogeneity of ground warming. We sampled a hillslope catena transect in a discontinuous permafrost zone near Fairbanks, Alaska to test the small-scale (1 to 3 meter) impacts of slope position and cover type on soil organic matter composition. Mineral active layer samples were collected from backslope, low backslope, and footslope positions at depths spanning 19 to 60 cm. We examined soil mineralogical composition, soil moisture, total carbon and nitrogen content, and organic mat thickness in conjunction with an assessment of SOC composition using Fourier-transform ion Cyclotron Resonance Mass Spectrometry (FT-ICR-MS). Soils in the footslope position had a higher relative contribution of lignin-like compounds while backslope soils had more aliphatic and condensed aromatic compounds as determined by FT-ICR-MS. The effect of open versus closed tree canopy cover varied with slope position. On the backslope, we found higher oxidation of molecules under open cover compared with closed cover, indicating an effect of warmer soil temperature on decomposition. Little to no effect of canopy was observed for soils at the footslope position, which we attributed, in part, to the strong impact of soil moisture content in SOC dynamics in the water-gathering footslope position. The thin organic mat under open cover on the backslope position may have contributed to differences in soil temperature and thus SOC oxidation under open and closed canopy. Here, the thinner organic mat did not appear to buffer the underlying soil against warm season air temperatures and thus increased SOC decomposition as indicated by higher oxidation of SOC molecules and a lower contribution of simple molecules under open cover compared with the closed canopy sites. Our findings suggest that the role of canopy cover in SOC dynamics varies as a function of landscape position and soil properties, namely organic mat thickness and soil moisture. Condition-specific heterogeneity of SOC composition under open and closed canopy cover highlights the protective effect of canopy cover for soils on backslope positions. This dataset contains a compressed (.zip) archive of the data and R scripts used for this manuscript. The dataset includes files in .csv format, which can be accessed and processed using MS Excel or R. This archive can also be accessed on GitHub at https://github.com/Erin-Rooney/Y1_fairbanks (DOI: 10.5281/zenodo.8071247).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2018

Impacts of Land Carbon Flux Variability on Variations in Atmospheric CO2

To better understand the role of atmospheric CO2 in the global carbon cycle, it is important to understand the factors that control its variability in time and space. Here we investigate these factors using a coupled land-atmosphere modeling system (the NASA GEOS-5 model) fitted with land carbon physics and atmospheric CO2 transport. To separate the influences of land carbon flux variability and atmospheric transport variability on the variability of CO2 concentrations, we conducted two carbon-fitted GEOS-5 AGCM simulations run in replay mode (a technique that guides the model's weather to match that of the MERRA-2 reanalysis): (i) a control simulation of 15-year carbon cycle dynamics and climate, and (ii) a simulation in which the climatological seasonal cycles of net biosphere production (NBP), as determined from diagnostics produced in the control experiment, are applied at the surface to the atmosphere instead. The impact of land flux variability on atmospheric CO2 variability is then isolated by subtracting the variability inherent in the second simulation (as induced, e.g., by strong El Nino conditions) from that in the control. These results are also presented in the context of supplemental simulations that have examined the impact of imposed drought in a free running AGCM simulation on carbon fluxes and atmospheric CO2 variability.

Lee, Eunjee↗
DOE OSTI2025

Multiphase computational fluid dynamics modeling of reacting flows in absorption columns for carbon capture

First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.

Absorption columns↗
DOE OSTI2025

Hydrology controls thermokarst and alters carbon cycling and methane emissions in peatlands near the southern limit of permafrost: Model Inputs and Outputs

This dataset contains model inputs and simulation outputs associated with the study “Hydrology controls thermokarst and alters carbon cycling and methane emissions in peatlands near the southern limit of permafrost.” The ecosys model was extended to represent thermokarst processes in laterally coupled two-dimensional transects and applied to the Scotty Creek basin (Northwest Territories, Canada). The ecosys model and documentation is available for download at https://github.com/jinyun1tang/ECOSYS. The dataset contains the inputs used to run the model for this study (meteorological forcing data (1950–2018) and site-specific soil and vegetation parameters) and model outputs used to generate Figure 3 and Figure 4 in the paper (simulated ground thermal and hydrologic states, vegetation dynamics, and carbon fluxes).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Uncertain Spatial Pattern of Future Land Use and Land Cover Change and Its Impacts on Terrestrial Carbon Cycle Over the Arctic–Boreal Region of North America

Land use and land cover change (LULCC) represents a key process of human-Earth system interaction and has profound impacts on terrestrial ecosystem carbon cycling. As a key input for ecosystem models, future gridded LULCC data is typically spatially downscaled from regional LULCC projections by integrated assessment models, such as the Global Change Analysis Model (GCAM). The uncertainty associated with the different spatial downscaling methods and its impacts on the subsequent model projections have been historically ignored and rarely examined. This study investigated this problem using two representative spatial downscaling methods and focused on their impacts on the carbon cycle over the Arctic-Boreal Vulnerability Experiment (ABoVE) domain, where extensive LULCC is expected. Specifically, we used the Future Land Use Simulation model (FLUS) and the Demeter model to generate 0.25° gridded LULCC data (i.e., LULCC FLUS and LULCC Demeter , respectively) with the same input of regional LULCC projections from GCAM, under both the low (i.e., SSP126) and high (i.e., SSP585) greenhouse gas emission scenarios. The two sets of downscaled LULCC were used to drive the Community Land Model version 5 and prognostically simulate the terrestrial carbon cycle dynamics over the 21st century. The results suggest large spatial-temporal differences between LULCC FLUS and LULCC Demeter , and the spatial distributions of the needleleaf evergreen boreal tree, broadleaf deciduous boreal tree, broadleaf deciduous boreal shrub, and C3 arctic grass are particularly different under both SSP126 and SSP585. Additionally, the spatiotemporal differences are larger under SSP126 than SSP585, due to more intensive LULCC under SSP126 than SSP585 from GCAM projection. The differences in LULCC further lead to large discrepancies in the spatial patterns of projected gross primary productivity, ecosystem respiration, and net ecosystem exchange, which represent more than 79% of the contributions of future LULCC in 2100. Additionally, the difference in carbon flux under SSP126 is generally larger than those under SSP585. This study highlights the importance of considering the uncertainties induced by the spatial downscaling process in future LULCC projections and carbon cycle simulations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Development of Carbon Flux Model Using ABI Data Over the Conterminous US

The satellite-driven carbon flux estimation has been playing important role to estimate continental-scale carbon budget. One of the biggest recent advances in the satellite-driven carbon flux modeling is utilization of high-frequent geostationary satellites to estimate diurnal cycle in carbon fluxes. The satellite based carbon flux estimation used the polar orbiting satellite sensors (e.g., Moderate Resolution Imaging Spectroradiometer (MODIS)), which allow us to observe target regions only once during the day. The new generation of geostationary satellite sensors provide frequent observations, often less than every 10 minutes. Here, we use GOES Advanced Baseline Imager (ABI) data to estimate hourly NEE over the conterminous US. We used the Terrestrial Observation Prediction System (TOPS) model for estimating hourly NEE. TOPS is a diagnostic ecosystem process model that simulates the fluxes of carbon and water through vegetation in response to climate variability. For the climate input, we developed hourly climate data using the same algorithm with NASA Earth Exchange Gridded Daily Meteorology (NEX-GDM) datasets based on machine learning techniques. The hourly climate data includes precipitation, maximum temperature, minimum temperature, dew point temperature, and solar radiation were derived from the Geostationary observations. The spatial patterns of ecosystem parameters used in TOPS are optimized using satellite Solar Induced Fluorescence (SIF) data. The high frequency GPP estimations from geostationary satellite sensors make it comparable to the instantaneous SIF data than daily GPP. We also used Ameriflux data for optimization of model parameters and the validation of the output. The derived data addresses the diurnal dynamics of carbon cycling at large scales and should help in reducing the uncertainties in carbon budget studies.

geostationary satellite↗
NASA NTRSDate not supplied

Hourly Carbon Fluxes Estimation Using the GOES Advanced Baseline Imager (ABI) Data Over the Conterminous USA

Tremendous efforts by Fluxnet scientists over the past few decades have made thousands of site-years of carbon flux observations available for advancing our understanding of carbon cycling in terrestrial ecosystems. One of key Fluxnet measurements is net ecosystem exchange (NEE) as it is directly related to carbon budget of terrestrial ecosystems. However, carbon flux estimation studies using satellite remote sensing have focused mainly on daily Gross Primary Production (GPP). The satellite based carbon flux estimation used the polar orbiting satellite sensors (e.g., Moderate Resolution Imaging Spectroradiometer (MODIS)), which allow us to observe target regions only once during the day. Because daily NEE is close to zero value, the carbon flux models using the polar orbiting satellite data have not been well used for NEE estimation. The new generation of geostationary satellite sensors (e.g., GOES Advanced Baesline Imager (ABI) and Himawari Advanced Himawari Imager (AHI)) provide frequent observations, often less than every 10 minutes. Here, we use GOES ABI data to estimate hourly NEE over the conterminous USA. We used the Terrestrial Observation Prediction System (TOPS) model for estimating hourly NEE. TOPS is a diagnostic ecosystem process model that simulates the fluxes of carbon and water through vegetation in response to climate variability. For the climate input, we developed hourly climate data using the same algorithm with NASA Earth Exchange Gridded Daily Meteorology (NEX-GDM) datasets based on machine learning techniques. The hourly climate data includes precipitation, maximum temperature, minimum temperature, dew point temperature, and, in particular, solar radiation that is directly derived from the Geostationary observations. The spatial patterns of ecosystem parameters used in TOPS are optimized using satellite Solar Induced Fluorescence (SIF) data. The high frequency GPP estimations from geostationary satellite sensors make it comparable to the instantaneous SIF data than daily GPP. We also used Fluxnet data for optimization of model parameters and the validation of the output. The derived data addresses the diurnal dynamics of carbon cycling at large scales and should help in reducing the uncertainties in carbon budget studies.

Geostationary satellite↗
NASA NTRS2005

Higher Order Chemistry Models in the CFD Simulation of Laser-Ablated Carbon Plumes

Production of single-walled carbon nanotubes (SWNT) has taken place for a number of years and by a variety of methods such-as laser ablation, chemical vapor deposition, and arc-jet ablation. Yet, little is actually understood about the exact chemical kinetics and processes that occur in SWNT formation. In recent time, NASA Johnson Space Center has devoted a considerable effort to the experimental evaluation of the laser ablation production process for SWNT originally developed at Rice University. To fully understand the nature of the laser ablation process it is necessary to understand the development of the carbon plume dynamics within the laser ablation oven. The present work is a continuation of previous studies into the efforts to model plume dynamics using computational fluid dynamics (CFD). The ultimate goal of the work is to improve understanding of the laser ablation process, and through that improved understanding, refine the laser ablation production of SWNT. Fig. 1 shows a basic schematic of the laser-ablation oven at NASA-JSC. Construction of the facility is simple in concept. Two concentric quartz tubes of 1.5 mm thickness form the inner and outer tubes with inside diameters of 2.2 and 5.08 cm respectively. At one end of the inner tube are located two 60 Hz pulsed lasers operating at 1064 nm and 532 nm wavelength with beam diameters of 5 mm aligned coaxially with the longitudinal axis of the inner quartz tube. For standard nanotube production runs, a 10 ns 532 nm pulse is followed 50 ns later by a 10 ns 1064 nm pulse. Each pulse is of 300 mJ energy. A target of carbon graphite with approximately 1% nickel and cobalt catalysts is located at the other end of the inner quartz tube. In the ordinary processing of SWNT, a base flow of 100 sccm of argon is maintained from the laser location and exits past the carbon target at a pressure of 66.7 kPa. These conditions yield a baseline mass flow through the chamber of 2.723x10(exp -6)kg/s of argon. The whole oven facility is heated to a temperature of 1473 K prior to nanotube production runs. Upon laser irradiation, part of the carbon target ablates immediately and forms a carbon vapor plume that penetrates into the argon base flow towards the laser initially at supersonic velocities. In the previous studies either a single carbon species, C3, was used to model the plume development, or a simplified 11 species carbon model reduced down from the Krestinin and Moravsky model for full fullerene chemistry was used. While both of these options yielded significant results, it was felt that the actual chemistry occ&g in the carbon plume might have a greater affect on the plume than assumed. Indeed, in the earlier attempts at modeling the carbon plume, several thermophysical characteristics could never be matched to experimental observations of plume development - mainly the propagation distance of the plume itself. In the present study, two additional chemistry models will be used to duplicate the previous studies simulations of the carbon plume. The first chemistry model used in this study is again a reduced form of the Krestinin and Moravsky rates. However the highest order carbon species allowed has been increased from C6 to C30 - therefore allowing the simulation of up to a half of the standard C60 fullerene. The second chemistry model investigated is a reduced form of a full carbon nanotube model developed at NASA-JSC. The C30 studies have already been accomplished at the present time, and the reduced SWNT model studies are currently underway. To pursue the current study, one sacrifice had to be made in that the simulation grid spacing had to be increased from 0.5 mm spacing to 1 mm spacing for the sake of computational efficiency since computational effort is proportional to the square of the number of grid points multiplied by the number of species considered. propagation that is far more in line with the experimental results observed by Puretzky et al as shown in Fig. 2. e C6 studies had yielded a far greater propagation in previous studies. In addition, chemical species development with the C30 model indicates that many higher order carbon species are produced outside of the plume proper (indicated by plotting contours of the background argon concentrations in Fig. 3) - this result was not observed in previous studies. In fact, some species primarily occurred outside of the plume itself - as shown for C27 in Fig. 4 when it is compared to Fig. 3. It could be asserted that this has occurred because all of the C27 in the plume had already been consumed in the formation of C30, but this does not seem to be indicated over time. Several other factors that arose in the previous studies have also been made more clear by the use of the higher order chemical models - one being that the use of c6 as an indicator species was mistaken. C6 is the only carbon species in the previous studies that was not injected into the flowfield as a boundary condition; it was therefore hoped that this species would provide insight into the formation of higher order carbon species for comparison to full SWNT production. But, when the plot of total mass in the plume is examined on a species by species basis in Fig. 5, it is seen that Cg was a fairly insignificant contributor to the total carbon mass in the plume and would not provide information on higher order carbon formation. the thermophysical characteristics of the carbon plume as well as simulate the carbon plume using the reduced SWNT model to provide an even better simulation of full chemistry effects upon plume propagation.

Scott, C. D.↗
DOE OSTI2024

Central transcriptional regulator controls photosynthetic growth and carbon storage in response to high light

Carbon capture and biochemical storage are some of the primary drivers of photosynthetic yield and productivity. To elucidate the mechanisms governing carbon allocation, we designed a photosynthetic light response test system for genetic and metabolic carbon assimilation tracking, using microalgae as simplified plant models. The systems biology mapping of high light-responsive photophysiology and carbon utilization dynamics between two variants of the same Picochlorum celeri species, TG1 and TG2 elucidated metabolic bottlenecks and transport rates of intermediates using instationary 13 C-fluxomics. Simultaneous global gene expression dynamics showed 73% of the annotated genes responding within one hour, elucidating a singular, diel-responsive transcription factor, closely related to the CCA1/LHY clock genes in plants, with significantly altered expression in TG2. Transgenic P. celeri TG1 cells expressing the TG2 CCA1/LHY gene, showed 15% increase in growth rates and 25% increase in storage carbohydrate content, supporting a coordinating regulatory function for a single transcription factor.

09 BIOMASS FUELS↗
DOE OSTI2025

Variational AutoEncoders Reveal Intensifying GPP Extremes in Continental United States based on CESM2 Simulations

Climate extremes significantly impact terrestrial carbon cycle dynamics, necessitating robust methods for detecting and analyzing anomalous behavior in plant productivity. This study presents a novel application of variational autoencoders (VAE) for identifying extreme events in gross primary productivity (GPP) from Community Earth System Model version 2 simulations across four AR6 regions in the Continental United States. We compare VAE-based anomaly detection with traditional singular spectral analysis (SSA) methods across three time periods: 1850-80, 1950-80, and 2050-80 under SSP5-8.5 scenario. The VAE architecture employs three dense layers and a latent space with input sequence length of 12 months, training on normalized GPP time series to reconstruct the GPP and identify anomalies based on reconstruction errors. Extreme events are defined using 5th percentile thresholds applied to both VAE and SSA anomalies. Results demonstrate strong regional agreement between VAE and SSA methods in spatial patterns of extreme event frequencies, despite VAE consistently producing higher threshold values (179-756 GgC for VAE vs. 100-784 GgC for SSA across regions and periods). Both methods reveal increasing magnitudes and frequencies of negative carbon cycle extremes toward 2050-80, particularly in Western and Central North America. The VAE approach shows comparable performance to established SSA techniques while offering computational advantages and enhanced capability for capturing non-linear temporal dependencies in carbon cycle variability. This research demonstrates the potential of deep learning approaches for extremes detection and provides a foundation for improved understanding of future carbon cycle risks under future conditions.

Sharma, Bharat [ORNL] (ORCID:0000000266982487)↗
DOE OSTI2024

Spin-Controllable Dynamics in Defect-Engineered Carbon Nanotubes as Single Photon Emitters: Data-Driven Modeling and Computations

Quantum technologies, such as quantum computing and sensing, require efficient single-photon emission (SPE) sources that operate at room temperature in telecom wavelengths. While several materials can serve as SPE sources, no single platform meets all the criteria for efficiency, ambient operation, and scalability. Single-walled carbon nanotubes (SWCNTs) with covalently attached molecules offer a promising solution. Their SPE can be easily tuned via modifications of the SWCNT's diameter, chirality, and bonded molecules, enabling emission across near-IR to telecom wavelengths at ambient conditions. However, to fully realize the potential of SWCNTs and unlock their quantum capabilities, a deeper understanding of how structural defects from molecular adducts affect their emission and competing photoexcited processes is essential. To address this gap in our knowledge, this project combined quantum chemistry calculations with data-driven methods of cheminformatics (QSAR) and machine learning (ML). The developed computational approaches have provided several design strategies for covalent functionalization of SWCNTs to improve their optical response. The collaboration with Los Alamos National Lab (LANL) enabled direct comparison of computational and experimental data, facilitating method validation. This partnership was enhanced through access to LANL's Center for Integrated Nanotechnologies (CINT) utilizing User Facility Program and summer internships, which provided three NDSU graduate students with hands-on experience at LANL. The outcomes of this project included (1) Advancing the current stage of computational methods in accurate modeling of non-adiabatic spin-dependent photoexcited dynamics and its applicability to nanosystems consisting of thousands of atoms, realized as open-access codes linked to existing DFT-based software; (2) Establishing the relationship between the structure of adducts and SWCNTs and intrinsic excitonic and spin properties of defect states for guiding novel synthetic strategies and experimental probes of chemically functionalized SWCNTs as near-IR emitting materials; (3) Generating virtual libraries of hypothetical functionalized SWCNTs for virtual screening of their chemical structures and optical properties, leveraging new functionalities of SWCNTs; (4) Offering a unique experience for NDSU graduate students that prepared them for future scientific careers related to materials modeling and big data processing. These results were summarized in 12 published journal papers and 3 recently submitted papers. One of a key finding is that the position of defect sites on the SWCNT surface primarily drives the emission redshift (up to 100 meV), while the polarity of the defect-inducing molecules has a much smaller effect (~10 meV). However, the electron-donating or withdrawing properties of a molecule influence selecting reactivity of defect sites. These insights important for optimizing synthetic protocols for desired emissions in SWCNTs. We also revealed that the interaction between two defects at various positions on the SWCNT enhances the redshift and optical activity of states, favoring strong near-IR emission. This suggests that manipulations in defect concentrations is a promising strategy for controlling efficient emission. Mostly important, the defect position was found controllable by the spin states of photoexcited intermediates: Excited aromatic molecules form ortho defects with SWCNTs at their singlet states in the presence of oxygen, while oxygen-free conditions favor para defects via the triplet-state mechanism. Additionally, a heat-activated [2+2] cycloaddition reaction facilitates divalent defect formation with fewer bonding positions that narrows emission bands. These groundbreaking findings have been experimentally validated and significantly advance our understanding of defect chemistry in SWCNTs. Using a novel encoding technique and 3D-MoRSE descriptors, we developed highly accurate ML/QSAR models to predict both the 3D structure and optical properties of SWCNTs with chemical defects. This model enabled the creation of a virtual library of 125,556 structures, providing new insights into the relationship between SWCNT-defect structure and emission.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS2020

Spatial Heterogeneity in CO2, CH4, and Energy Fluxes: Insights from Airborne Eddy Covariance Measurements Over the Mid-Atlantic Region

The exchange of carbon between the Earth’s atmosphere and biosphere influences the atmospheric abundances of carbon dioxide (CO2) and methane (CH4). Airborne eddy covariance can quantify surface-atmosphere exchange from landscape-to-regional scales, offering a unique perspective on carbon cycle dynamics. We use extensive airborne measurements to quantify fluxes of sensible heat, latent heat, CO2, and CH4across multiple ecosystems in the Mid-Atlantic region during September 2016 and May 2017. In conjunction with footprint analysis and land cover information, we use the airborne dataset to explore the effects of landscape heterogeneity on measured fluxes. Our results demonstrate large variability in CO2 uptake over mixed agricultural and forested sites, with fluxes ranging from -3.4 ± 0.7 to -11.5 ± 1.6μmol m-2s-1 for croplands and -9.1 ± 1.5 to -22.7 ± 3.2μmol m-2s-1for forests. We also report substantial CH4emissions of 32.3 ± 17.0to 76.1 ± 29.4nmol m-2s-1 from a brackish herbaceous wetland and 58.4± 12.0 to 181.2 ± 36.8nmol m-2s-1 from a freshwater forested wetland. Comparison of ecosystem-specific aircraft observations with measurements from eddy covariance flux towers along the flight path demonstrate that towers capture ~30–75% of the regional variability in ecosystem fluxes. Diel patterns measured at the tower sites suggest that peak, midday flux measurements from aircraft accurately predicts net daily CO2exchange. We discuss next steps in applying airborne observations to evaluate bottom-up flux models and improve understanding of the biophysical processes that drive carbon exchange from landscape-to-regional scales.

airborne measurements to quantify fluxes of sensib↗
DOE OSTI2022

Spatial turnover of soil viral populations and genotypes overlain by cohesive responses to moisture in grasslands

Viruses shape microbial communities, food web dynamics, and carbon and nutrient cycling in diverse ecosystems. However, little is known about the patterns and drivers of viral community composition, particularly in soil, precluding a predictive understanding of viral impacts on terrestrial habitats. To investigate soil viral community assembly processes, here we analyzed 43 soil viromes from a rainfall manipulation experiment in a Mediterranean grassland in California. We identified 5,315 viral populations (viral operational taxonomic units [vOTUs] with a representative sequence ≥10 kbp) and found that viral community composition exhibited a highly significant distance–decay relationship within the 200-m 2 field site. This pattern was recapitulated by the intrapopulation microheterogeneity trends of prevalent vOTUs (detected in ≥90% of the viromes), which tended to exhibit negative correlations between spatial distance and the genomic similarity of their predominant allelic variants. Although significant spatial structuring was also observed in the bacterial and archaeal communities, the signal was dampened relative to the viromes, suggesting differences in local assembly drivers for viruses and prokaryotes and/or differences in the temporal scales captured by viromes and total DNA. Despite the overwhelming spatial signal, evidence for environmental filtering was revealed in a protein-sharing network analysis, wherein a group of related vOTUs predicted to infect actinobacteria was shown to be significantly enriched in low-moisture samples distributed throughout the field. Overall, our results indicate a highly diverse, dynamic, active, and spatially structured soil virosphere capable of rapid responses to changing environmental conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Enriched Molecular-Level View of Saline Wetland Soil Carbon by Sensitivity-Enhanced Solid-State NMR

Soil organic matter (SOM) plays a major role in mitigating greenhouse gas emission and regulating earth’s climate, carbon cycle, and biodiversity. Wetland soils account for one-third of all SOM; however, globally, coastal wetland soils are eroding faster due to increasing sea-level rise. Our understanding of carbon sequestration dynamics in wetlands lags behind that of upland soils. Here, we employ solid-state nuclear magnetic resonance (ssNMR) to investigate the molecular-level structure of biopolymers in wetland soils spanning 11 centuries. High-resolution multidimensional spectra, enabled by dynamic nuclear polarization (DNP), demonstrate enduring preservation of molecular structures within herbaceous plant cores, notably condensing aromatic motifs and carbohydrates, even over a millennium, with the preserved cores constituting a decreasing minority among molecules from decomposition and repolymerization with depth and age. Such preserved cores occur alongside molecules from the decomposition of loosely packed parent biopolymers. These findings emphasize the relative vulnerability of coastal wetland SOM when exposed to oxygenated water due to geological and anthropogenic changes.

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

Predicted annual biomass input to maintain soil organic carbon under contrasting management

Corn (Zea mays L.) residue is an important source of lignocellulosic feedstocks for meeting the US renewable energy goals. However, there is uncertainty in how much biomass can be sustainably harvested. The objectives of this study were to (i) simulate soil organic carbon (SOC) dynamics in the 0–20, 20–30, and 0–30 cm soil depths for the 2005–2015 study period using CQESTR, a process-based C model; (ii) establish trends in SOC with residue removal and tillage management over the predictive period of 2016–2035; and (iii) recommend practices to maintain SOC stocks in a corn–soybean (Glycine max L. [Merr.]) rotation with corn residue removal in western Minnesota. The residue removal treatments in the conventional tillage since 1995 (CT 95 ) and no tillage since 1995 (NT 95 ) were initiated in 2005, and in 2006 for the no tillage since 2005 (NT 05 ). Corn residue removal treatments were none, cob removal, ~50% removal, and ~75% removal (R0, RCob, R50, and R75, respectively). Only R0 in NT 05 increased SOC (0.32 Mg C ha –1 yr –1 ) stocks in the top 30 cm soil depth by 2035 with annual aboveground biomass input of 6.30 Mg ha –1 yr –1 . It was concluded that only the NT 05 field has the potential to maintain SOC at the current rate of biomass addition, and a minimal annual aboveground C (MSC) input of 3.75 Mg C ha –1 yr –1 is required to maintain SOC until 2035. These simulation results indicate that caution should be used in removing residue especially under conventional tillage.

54 ENVIRONMENTAL SCIENCES↗