Search NASA⌕ Search

SEARCH · Search NASA

Results for “CO2 Flux”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Hydroxide Exchange Membrane Carbon Capture (HEMCC) Using Nickel Hydroxide Batteries and Flow-through Membranes

Proposed is an electrochemical nickel hydroxide based hydroxide exchange membrane carbon capture (HEMCC) device for Direct Air Capture (DAC) of CO2. DAC has been identified as one of the key net negative carbon technologies to achieve net zero carbon emissions. Net negative carbon technologies are required to offset continued emissions from dilute CO2 sources such as agriculture and construction. The majority of current DAC technologies at scale (>1 KT∙yr-1) are adsorbent based technologies with significant energy cost. The traditional DAC energy cost is primarily driven by the temperature swing required to regenerate the sorbent and has been shown to be 1.8 MWh·ton-1 at the system level. Electrochemical pH swing devices are a growing research area for carbon capture devices with the goal of lowering the energy cost required for DAC. A pH gradient is built by generating OH- at the cathode and consuming OH- at the anode. An acid-base equilibrium with CO2 allows for the capture of CO2 at the cathode and release at the anode. This extends from other electrochemical CO2 capture devices based on pKa shifts of an electrochemically active organic species allowing for the capture and release of CO2. Electrochemical CO2 capture is considered promising based on potentially low energy costs to capture CO2 in comparison with current temperature swing adsorption technologies. This work explores Ni(OH)2 electrodes to produce the pH gradient for CO2 capture and release. At the cathode NiOOH is reduced to Ni(OH)2 while at the anode Ni(OH)2 is oxidized to NiOOH. The symmetrical electrodes allow for a low voltage requirement; the thermodynamic potential difference of standard electrochemical reactions is zero. Most of the voltage observed is to produce the pH gradient with the remainder driving the polarization of the electrodes. There is a resistance component as well, but this is small in comparison due to the low current densities used in the device, nominally 2 mA·cm-1. Two similar devices are presented, a traditional MEA (membrane electrode assembly) and a flow-through MEA. The traditional MEA separates the two Ni(OH)2 electrodes with an 80μm Piperion® membrane. While the flow-through membrane separates the electrodes with a three piece membrane consisting of two 80μm Piperion® membranes with a porous membrane between them. In the traditional MEA system air is passed over the cathode for capture, while the flow-through MEA the air is passed through the porous membrane isolated from the electrodes. The traditional MEA has been used to establish a baseline performance of the device and has been shown to capture CO2 at an energy cost of 1 MWh·ton-1 at the device level. An understanding has been built around the components of that energy cost including the relationship of flux to current density, effect of a regeneration process, transient battery behavior, and gas losses coinciding with changing the polarization of the batteries. The flow-through MEA looks to address of transient battery behavior and gas losses. It allows for denser, higher capacity electrodes, which can lean on traditional Ni-MH battery technology used in alkaline batteries used today. The higher capacities, limit the transient battery effect on flux in the device. Gas losses are addressed by having a continuous inlet air stream to the device and continuous outlet product.

Buchen, James↗
DOE OSTI2022

Root system architecture and environmental flux analysis in mature crops using 3D root mesocosms

Current methods of root sampling typically only obtain small or incomplete sections of root systems and do not capture their true complexity. To facilitate the visualization and analysis of full-sized plant root systems in 3-dimensions, we developed customized mesocosm growth containers. While highly scalable, the design presented here uses an internal volume of 45 ft 3 (1.27 m 3 ), suitable for large crop and bioenergy grass root systems to grow largely unconstrained. Furthermore, they allow for the excavation and preservation of 3-dimensional root system architecture (RSA), and facilitate the collection of time-resolved subterranean environmental data. Sensor arrays monitoring matric potential, temperature and CO 2 levels are buried in a grid formation at various depths to assess environmental fluxes at regular intervals. Methods of 3D data visualization of fluxes were developed to allow for comparison with root system architectural traits. Following harvest, the recovered root system can be digitally reconstructed in 3D through photogrammetry, which is an inexpensive method requiring only an appropriate studio space and a digital camera. We developed a pipeline to extract features from the 3D point clouds, or from derived skeletons that include point cloud voxel number as a proxy for biomass, total root system length, volume, depth, convex hull volume and solidity as a function of depth. Ground-truthing these features with biomass measurements from manually dissected root systems showed a high correlation. We evaluated switchgrass, maize, and sorghum root systems to highlight the capability for species wide comparisons. We focused on two switchgrass ecotypes, upland (VS16) and lowland (WBC3), in identical environments to demonstrate widely different root system architectures that may be indicative of core differences in their rhizoeconomic foraging strategies. Finally, we imposed a strong physiological water stress and manipulated the growth medium to demonstrate whole root system plasticity in response to environmental stimuli. Hence, these new “3D Root Mesocosms” and accompanying computational analysis provides a new paradigm for study of mature crop systems and the environmental fluxes that shape them.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗
DOE OSTI2023

Biome-scale temperature sensitivity of ecosystem respiration revealed by atmospheric CO2 observations

Abstract The temperature sensitivity of ecosystem respiration regulates how the terrestrial carbon sink responds to a warming climate but has been difficult to constrain observationally beyond the plot scale. Here we use observations of atmospheric CO 2 concentrations from a network of towers together with carbon flux estimates from state-of-the-art terrestrial biosphere models to characterize the temperature sensitivity of ecosystem respiration, as represented by the Arrhenius activation energy, over various North American biomes. We infer activation energies of 0.43 eV for North America and 0.38 eV to 0.53 eV for major biomes therein, which are substantially below those reported for plot-scale studies (approximately 0.65 eV). This discrepancy suggests that sparse plot-scale observations do not capture the spatial-scale dependence and biome specificity of the temperature sensitivity. We further show that adjusting the apparent temperature sensitivity in model estimates markedly improves their ability to represent observed atmospheric CO 2 variability. This study provides observationally constrained estimates of the temperature sensitivity of ecosystem respiration directly at the biome scale and reveals that temperature sensitivities at this scale are lower than those based on earlier plot-scale studies. These findings call for additional work to assess the resilience of large-scale carbon sinks to warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Host analysis-guided selection and targeted engineering (HASTE) of Lipomyces tetrasporus for the conversion of CO2-derived feedstocks

Efficient and cost-competitive bioproduction calls for utilizing CO2-derived feedstocks, such as products from electro-reduction of CO2 and hydrolysate from lignocellulosic biomass. However, efficiently using all their carbon components, including acetate, glucose, and xylose, remains a challenge. Here, we characterize Lipomyces tetrasporus, a novel, robust yeast strain capable of effectively assimilating these carbon sources. We used an integrated systems biology approach combining ¹³C metabolic flux analysis, dynamic labeling experiments, and RNA sequencing. We conducted the first metabolic flux analysis for glucose, xylose, and acetate catabolism in this species. Dynamic labeling revealed a highly active TCA cycle during acetate metabolism, evidenced by rapid citrate and malate accumulation. The strain demonstrated strong NADH/NADPH production and acetyl-CoA synthase activity. Using insights and gene targets from this analysis, we engineered L. tetrasporus for malate production. The engineered strain produced 7.5 g/L malic acid (0.25 g/g yield) in shake flasks with glucose-acetate media and 28.8 g/L malic acid at a yield of 0.20 g/g in fed-batch mode with corn-stover hydrolysate. Together, these insights and rational strain engineering establish L. tetrasporus as a versatile, Crabtree-negative platform that is an energy-CO2-bioproduction nexus for channeling CO2 carbon into value-added bioproducts.

Xiao, Zhengyang↗
DOE OSTI2024

A multi-model assessment of the early last deglaciation (PMIP4 LDv1): a meltwater perspective

Abstract. The last deglaciation (∼20–11 ka BP) is a period of a major, long-term climate transition from a glacial to interglacial state that features multiple centennial- to decadal-scale abrupt climate variations whose root cause is still not fully understood. To better understand this time period, the Paleoclimate Modelling Intercomparison Project (PMIP) has provided a framework for an internationally coordinated endeavour in simulating the last deglaciation whilst encompassing a broad range of models. Here, we present a multi-model intercomparison of 17 transient simulations of the early part of the last deglaciation (∼20–15 ka BP) from nine different climate models spanning a range of model complexities and uncertain boundary conditions and forcings. The numerous simulations available provide the opportunity to better understand the chain of events and mechanisms of climate changes between 20 and 15 ka BP and our collective ability to simulate them. We conclude that the amount of freshwater forcing and whether it follows the ice sheet reconstruction or induces an inferred Atlantic meridional overturning circulation (AMOC) history, heavily impacts the deglacial climate evolution for each simulation rather than differences in the model physics. The course of the deglaciation is consistent between simulations except when the freshwater forcing is above 0.1 Sv – at least 70 % of the simulations agree that there is warming by 15 ka BP in most places excluding the location of meltwater input. For simulations with freshwater forcings that exceed 0.1 Sv from 18 ka BP, warming is delayed in the North Atlantic and surface air temperature correlations with AMOC strength are much higher. However, we find that the state of the AMOC coming out of the Last Glacial Maximum (LGM) also plays a key role in the AMOC sensitivity to model forcings. In addition, we show that the response of each model to the chosen meltwater scenario depends largely on the sensitivity of the model to the freshwater forcing and other aspects of the experimental design (e.g. CO2 forcing or ice sheet reconstruction). The results provide insight into the ability of our models to simulate the first part of the deglaciation and how choices between uncertain boundary conditions and forcings, with a focus on freshwater fluxes, can impact model outputs. We can use these findings as helpful insight in the design of future simulations of this time period.

Snoll, Brooke (ORCID:0000000230914394)↗
DOE OSTI2023

FUN-BioCROP model with litter decomposition parameters derived from the LIDET dataset

This repository contains the code and data necessary to run the FUN-BioCROP (Fixation and Uptake of Nitrogen-Bioenergy Carbon, Rhizosphere, Organisms, and Protection) model with litter decomposition parameters derived from a modified Monte Carlo simulation that used the Long-term Intersite Decomposition Experiment Team dataset. Related publication:Juice, S.M., Ridgeway, J.R., Hartman, M.D., Parton, W.J., Berardi, D.M., Sulman, B.N., Allen, K.E., & Brzostek, E.R. Reparameterizing litter decomposition using a simplified Monte Carlo method improves litter decay simulated by a microbial model and alters bioenergy soil carbon estimates. Description of Files: FUNBioCROP_LIDET Study.Rmd R code with FUN-BioCROP model that can be run with 10 different sets of parameters for litter decomposition (Baseline, LIDET, or eight other best parameter sets identified in the modified Monte Carlo simulation. CORPSE Functions_Bioenergy_V2.R Code with CORPSE model functions, called by FUNBioCROP_LIDET Study.Rmd Model Input Data: bulk.csv, bulk_till.csv, rhizo.csv, rhizo_till.csv, litter.csv Initial C and N (kg C or N/m2) pool values for each soil compartment, final values from spin up. All five files have the same columns: (Column - Description - Units) uFastC - Unprotected fast decomposing carbon - kg carbon/m2 uSlowC - Unprotected slow decomposing carbon - kg carbon/m2 uNecroC - Unprotected necromass carbon - kg carbon/m2 pFastC - Protected fast decomposing carbon - kg carbon/m2 pSlowC - Protected slow decomposing carbon - kg carbon/m2 pNecroC - Protected necromass carbon - kg carbon/m2 livingMicrobeC - Carbon in living microbial biomass - kg carbon/m2 uFastN - Unprotected fast decomposing nitrogen - kg nitrogen/m2 uSlowN - Unprotected slow decomposing nitrogen - kg nitrogen/m2 uNecroN - Unprotected necromass nitrogen - kg nitrogen/m2 pFastN - Protected fast decomposing nitrogen - kg nitrogen/m2 pSlowN - Protected slow decomposing nitrogen - kg nitrogen/m2 pNecroN - Protected necromass nitrogen - kg nitrogen/m2 inorganicN - Inorganic nitrogen - kg nitrogen/m2 CO2 - Carbon in carbon dioxide - kg carbon/m2 livingMicrobeN - Nitrogen in living microbial biomass - kg nitrogen/m2 Model Input Data: FluxTower_AvgSoilT.csv: Average daily soil temperature (oC) at 10 cm depth at University of Illinois Urbana-Champaign (UIUC) Energy Farm flux tower from 7/2008-3/2016. (One year of averaged data) Model Input Data: FluxTower_AvgSoilVWC.csv: Average daily soil volumetric water content (VWC) at 10 cm depth at UIUC Energy Farm flux tower from 7/2008-3/2016. (One year of averaged data) Model Input Data: input_CCS_LIDET Study.csv: This file has daily data to run FUN-BioCROP (Column - Description - Units): yr - calendar year - year doy - day of year (1 to 365) (no leap year) - day anpp - aboveground NPP (DayCent) - kg C/m2/day bnpp - belowground NPP (DayCent) - kg C/m2/day aglivc - live aboveground biomass carbon (DayCent) - kg C/m2 bglivcj - live juvenile fine root biomass carbon (DayCent) - kg C/m2 bglivcm - live mature fine root biomass carbon (DayCent) - kg C/m2 aglivn - live aboveground biomass nitrogen (DayCent) - kg N/m2 bglivnj - live juvenile fine root biomass nitrogen (DayCent) - kg N/m2 bglivnm - live mature fine root biomass nitrogen (DayCent) - kg N/m2 nyr - simulation year - year cult - indicates a cultivation event (0 or 1) crop - indicates a new crop (0 or 1) fert - indicates a fertilizer event (0 or 1) frst - indicates the first day of the growing season (0 or 1) harv - indicates a harvest event (0 or 1) last - indicates the end of the growing season (0 or 1) croptype - crop type (0=none; 1=alfalfa; 2=corn; 3=grass clover pasture; 4=soybean; 5=wheat) cropsrl - crop specific root length - mm/g root cultrhizmix - fraction of rhizosphere mixed with bulk soil during cultivation (0.0-1.0) - fraction cultlitmix - fraction of litter mixed with bulk soil during cultivation (0.0-1.0) - fraction harvremov - fraction of above ground biomass removed during harvest (0.0-1.0) - fraction fertamt - fertilization amount - g N/m2 lifehist - plant life history (0 = annual, 1 = perennial) froot_turnover_c - amount of C in fine root turnover - kg C/m2 froot_turnover_n - amount of N in fine root turnover - kg N/m2 agrd_turnover_c - amount of C in aboveground biomass turnover - kg C/m2 agrd_turnover_n - amount of N in aboveground biomass turnover - kg N/m2 leaf_litter_fastfrac - Fast decomposing fraction of leaf litter (0.0-1.0) - fraction root_litter_fastfrac - Fast decomposing fraction of root litter (0.0-1.0) - fraction root_diameter - root diameter - mm root_length - root length - mm root/m2 rhizo_frac - fraction of total soil volume that is rhizosphere (0.0 - 1.0) - fraction date - date in format YYYY-MM-DD Instructions: Save the model code ("FUN-BioCROP_LIDET Study.Rmd") and accompanything files (data streams and CORPSE function code) in the same folder. In model code "Chunk 3: Load CORPSE Data Streams" set the working directory (setwd) to the folder with the files saved in step #1. In "Chunk 5: Define LIDET parameter sets" select the litter decomposition parameter set to be used in the run, and comment out all other sets. If changing any parameter values, edit them in "Chunk 6: Load parameters." Run all chunks up to and including "Chunk 10: Prepare Data for Export." In "Chunk 11: Export Output Data" edit data frames for export and filenames, as necessary. "Chunk 12: Graph Total Soil C" makes a figure of C remaining over the model run period. Description of each model chunk (in file FUN-BioCROP_LIDET Study.Rmd): Chunk 1: Remove all functions, clear memory. Removes all functions from R environment, clears the memory. Chunk 2: Load Packages. Loads packages necessary to run the code. Chunk 3: Load CORPSE Data Streams. Sets the working directory and loads the data files necessary to run CORPSE. Chunk 4: Load CORPSE Functions. Loads the R script with CORPSE functions from the working directory, "CORPSE Functions_Bioenergy_V2.R". Chunk 5: Define LIDET parameter sets. Has ten different parameter sets for litter decomposition tested in this study: Baseline parameters, LIDET parameters, and the other 8 best performing parameter sets identified in the modified Monte Carlo. To run the model, all but one parameter set must be commented out. Chunk 6: Load Parameters. Loads all fixed parameters to run the model. Data frame with definitions of parameters is in the CORPSE function script "CORPSE Functions_Bioenergy_V2.R" Chunk 7: Prepare Data Streams. Takes data streams loaded in Chunk 3 and puts them in the format necessary to run the model. The model is coded to run at least two sites at a time, so if only one site is being run it must be run in duplicate. Individual data tables of daily values are created in this chunk from the input data file. Chunk 8: Set Initial Conditions. Creates data tables of soil C and N pools for each soil compartment (rhizo_till, rhizo, bulk_till, bulk, litter) and loads initial values into the data tables. Creates lists for each soil compartment to hold model output. Chunk 9: Load FUN Data and Set Up Matrices. Uses DayCent data to calculate FUN input data: root and leaf N demand, total N demand, plant CN, leaf N available for retranslocation, and litter production. Creates matrices for FUN model outputs. Chunk 10: Run Model. Runs the model. Chunk 11: Prepare Data for Export. Combines data from each day saved as lists into data frames for each soil compartment. Adds values from all soil compartments together to calculate total soil values, creates separate data frames for each soil C and N pool (e.g., protected slow C) for the total soil value. Adds different C and N pools together to calculate total soil C and N for all layers. Creates data frame of ratio of protected to unprotected SOC. Organizes FUN data for export. Chunk 12: Export Results. Exports CSV files of model results to the working directory. Chunk 13: Graph Total Soil C. Makes figure of C remaining over time. Related Links: Original FUN-BioCROP model: https://github.com/BrzostekEcologyLab/FUN-BioCROP LIDET dataset: https://andlter.forestry.oregonstate.edu/data/abstract.aspx?dbcode=TD023

Juice, Stephanie↗
DOE OSTI2024

Biogeochemistry simulations for the Salt Marsh Accretion Response to Temperature eXperiment (SMARTX)

Coastal ecosystems have been largely ignored in Earth system models but are important zones for carbon and nutrient processing. Interactions between water, microbes, soil, sediments, and vegetation are important for mechanistic representation of coastal processes and ecosystem function. To investigate the role of these feedbacks, we used a reactive transport model (PFLOTRAN) that has the capability to be connected to the Energy Exascale Earth System Model (E3SM). PFLOTRAN was used to incorporate redox reactions and track chemical species important for coastal ecosystems as well as define simple representations of vegetation dynamics. Our goal was to incorporate oxygen flux, salinity, pH, sulfur cycling, and methane production along with plant-mediated transport of gases and tidal flux. Using porewater profile and incubation data for model calibration and evaluation, we were able to create depth-resolved biogeochemical soil profiles for saltmarsh habitat and use this updated representation to simulate direct and indirect effects of elevated CO2 and temperature on subsurface biogeochemical cycling. We found that simply changing the partial pressure of CO2 or increasing temperature in the model did not fully reproduce observed changes in the porewater profile, but the inclusion of plant or microbial responses to CO2 and temperature manipulations was more accurate in representing porewater concentrations. This indicates the importance of characterizing tightly coupled vegetation-subsurface processes for developing predictive understanding and the need for measurement of plant-soil interactions on the same time scale to understand how hotspots or moments are generated.Included in this data package are PFLOTRAN input (PFLOTRAN input files and chemical database) files for simulating single column biogeochemistry, root, and tide interactions at the Global Change Research Wetland (Kirkpatrick Marsh; Edgewater, MD). The biogeochemical network includes soil organic matter decomposition, nitrogen, iron, and sulfur cycling, and methanogenesis. Reduced species can be oxidized and plant processes include oxygen and nutrient priming, methane release, and nutrient uptake.Inputs:TAI_database.dat - geochemical database for reactions, more information on database structure and variables can be found here https://www.pflotran.org/documentation/user_guide/cards/pages/geochemical_database.htmlswamp.in - input file for biogeochemical network in PFLOTRANswamp_eCO2.in - input file for biogeochemical network in PFLOTRAN with input gas partial pressures/concentrations adjusted for elevated CO2 treatmentsOutputs:swamp_obs_0.tec - hourly porewater concentrations from from multiple depths in the soil columnswamp_eCO2_obs_0.tec - hourly porewater concentrations from multiple depths in the soil column for elevated CO2 treatmentsPFLOTRAN code access: https://github.com/fmyuan/pflotran-elm-interface.git

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Exploring structure-function relationships governing transport by and regulation of the major cyanobacterial bicarbonate transporters SbtA and BicA

Cyanobacteria play a large and essential role in the global carbon cycle, while also representing promising microbial platforms for the direct bioproduction of fuels and chemicals from solar energy and carbon dioxide (CO 2 ). In both cases, a critical step in the overall photosynthetic process is the transport of inorganic carbon (Ci) across the cell membrane, which occurs via different mechanisms and associated transporter proteins that target either CO 2 or bicarbonate (HCO 3 - ). In the latter case, BicA (constitutively-expressed, high-flux/low-affinity Na + /HCO 3 - symporter) and SbtA (low-flux/high-affinity Na + /HCO 3 - symporter, expressed only under low CO 2 conditions) have been identified as the two main cyanobacterial HCO 3 - transporters. This project entailed combining detailed functional and high-resolution structural studies to help fill key knowledge gaps regarding the mechanisms of HCO 3 - transport by BicA and SbtA. Specifically, major goals of this project included: i) identifying key residues and structural features of BicA and SbtA that are responsible for modulating their uptake flux and affinity behaviors, and ii) obtaining high-resolution structures of full-length BicA and SbtA.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

A comprehensive strategy for stable, high productivity cultivation of microalgae with controllable biomass composition

This project was performed primarily as a partnership between three Universities – The University of Toledo (UToledo), Montana State University (MSU) and The University of North Carolina at Chapel Hill (UNC). North Carolina State University also contributed to the project through a sub-subcontract. The project developed fundamental and applied knowhow for cultivation systems that are not restricted by (i) proximate availability of flue gas or other high concentration CO 2 sources and (ii) the energy and infrastructure burden to deliver CO 2 over long distances. To overcome these constraints, this project focused on cultivation of alkaliphilic algal cultures in high pH, high alkalinity media (pH>10, alkalinity >100mEq) and demonstrated high biomass productivity, improvement of biomass composition for fuel production, and culture stability. The high media alkalinity ensured availability of sufficient HCO 3 - in solution for high rates of carbon fixation while the high pH of the media facilitated rapid scavenging of atmospheric CO 2 . Overall, the direct air capture rates (DAC) in our system were sufficient to achieve high biomass productivities using atmospheric CO 2 alone. The project also developed and validated toolkits related to metabolic flux modeling and microbiome characterization.

09 BIOMASS FUELS↗
DOE OSTI2025

Reactivity of Shale to Supercritical CO 2 : Insights from Microstructural Characterization and Mineral Phase Evolution in Caney Shales for CCUS Applications

Understanding mineral–fluid interactions in shale under supercritical CO 2 (scCO 2 ) conditions is relevant for assessing long-term geochemical containment. This study characterizes mineralogical transformations and elemental redistribution in five Caney Shale samples serving as proxies for reservoir (R1, R2, R3) and caprock (D1, D2) facies, subjected to 30-day static exposure to pure scCO 2 at 60 °C and 17.23 MPa (2500 psi), with no brine or impurities introduced. SEM-EDS analyses were conducted before and after exposure, with mineral phases classified into silicates, carbonates, sulfides, and organic matter. Initial compositions were dominated by quartz (38–47 wt.%), illite (16–23 wt.%), carbonates (12–18 wt.%), and organic matter (8–11 wt.%). Post-exposure, carbonate loss ranged from 15 to 40% in reservoir samples and up to 20% in caprock samples. Illite and K-feldspar showed depletion of Fe 2+ , Mg 2+ , and K + at grain edges and cleavages, while pyrite underwent oxidation with Fe redistribution. Organic matter exhibited scCO 2 -induced surface alteration and apparent sorption effects, most pronounced in R2 and R3. Elemental mapping revealed Ca 2+ , Mg 2+ , Fe 2+ , and Si 4+ mobilization near reactive interfaces, though no secondary mineral precipitates formed. Reservoir samples developed localized porosity, whereas caprock samples retained more structural clay integrity. The results advance understanding of mineral reactivity and elemental fluxes in shale-based CO 2 sequestration.

geochemical CO2 sequestration↗
DOE OSTI2021

Spatial distributions of X CO 2 seasonal cycle amplitude and phase over northern high-latitude regions

Satellite-based observations of atmospheric carbon dioxide (CO 2 ) provide measurements in remote regions, such as the biologically sensitive but undersampled northern high latitudes, and are progressing toward true global data coverage. Recent improvements in satellite retrievals of total column-averaged dry air mole fractions of CO 2 (X CO 2 ) from the NASA Orbiting Carbon Observatory 2 (OCO-2) have allowed for unprecedented data coverage of northern high-latitude regions, while maintaining acceptable accuracy and consistency relative to ground-based observations, and finally providing sufficient data in spring and autumn for analysis of satellite-observed X CO 2 seasonal cycles across a majority of terrestrial northern high-latitude regions. Here, we present an analysis of X CO 2 seasonal cycles calculated from OCO-2 data for temperate, boreal, and tundra regions, subdivided into 5° latitude by 20° longitude zones. We quantify the seasonal cycle amplitudes (SCAs) and the annual half drawdown day (HDD). OCO-2 SCAs are in good agreement with ground-based observations at five high-latitude sites, and OCO-2 SCAs show very close agreement with SCAs calculated for model estimates of X CO 2 from the Copernicus Atmosphere Monitoring Services (CAMS) global inversion-optimized greenhouse gas flux model v19r1 and the CarbonTracker2019 model (CT2019B). Model estimates of X CO 2 from the GEOS-Chem CO 2 simulation version 12.7.2 with underlying biospheric fluxes from CarbonTracker2019 (GC-CT2019) yield SCAs of larger magnitude and spread over a larger range than those from CAMS, CT2019B, or OCO-2; however, GC-CT2019 SCAs still exhibit a very similar spatial distribution across northern high-latitude regions to that from CAMS, CT2019B, and OCO-2. Zones in the Asian boreal forest were found to have exceptionally large SCA and early HDD, and both OCO-2 data and model estimates yield a distinct longitudinal gradient of increasing SCA from west to east across the Eurasian continent. In northern high-latitude regions, spanning latitudes from 47 to 72° N, longitudinal gradients in both SCA and HDD are at least as pronounced as latitudinal gradients, suggesting a role for global atmospheric transport patterns in defining spatial distributions of X CO2 seasonality across these regions. GEOS-Chem surface contact tracers show that the largest X CO 2 SCAs occur in areas with the greatest contact with land surfaces, integrated over 15–30d. The correlation of X CO 2 SCA with these land surface contact tracers is stronger than the correlation of X CO 2 SCA with the SCA of CO 2 fluxes or the total annual CO 2 flux within each 5° latitude by 20° longitude zone. This indicates that accumulation of terrestrial CO 2 flux during atmospheric transport is a major driver of regional variations in X CO 2 SCA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Real-Time Partitioning of Diurnal Stem CO 2 Efflux into Local Stem Respiration and Xylem Transport Processes

The apparent respiratory quotient (ARQ) of tree stems, defined as the ratio of net stem CO 2 efflux (E S_CO2 ) to net stem O 2 influx (E S_O2 ), offers insights into the balance between local respiratory CO 2 production and CO 2 transported via the xylem. Traditional static chamber methods for measuring ARQ can introduce artifacts and obscure natural diurnal variations. Here, we employed an open flow-through stem chamber with ambient air coupled with cavity ring-down spectrometry, which uses the molecular properties of CO 2 and O 2 molecules to continuously measure E S_CO2 , E S_O2 , and ARQ, at the base of a California cherry tree (Prunus ilicifolia) during the 2024 growing season. Measurements across three stem chambers over 3–11-day periods revealed strong correlations between E S_CO2 and E S_O2 and mean ARQ values ranging from 1.3 to 2.9, far exceeding previous reports. Two distinct diurnal ARQ patterns were observed: daytime suppression with nighttime recovery, and a morning peak followed by gradual decline. Partitioning E S_CO2 into local respiration and xylem-transported CO 2 indicated that the latter can dominate when ARQ exceeds 2.0. Furthermore, transported CO 2 exhibited a higher temperature sensitivity than local respiration, with both processes showing declining temperature sensitivity above 20 °C. These findings underscore the need to differentiate stem CO 2 flux components to improve our understanding of whole-tree carbon cycling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Rational engineering to enhance C8-fatty acid biosynthesis in Picosynechococcus sp. PCC 7002

Medium-chain fatty acids (e.g., C8-C14) are important biofuel precursors that can be produced from CO2 by cyanobacteria and other photoautotrophs. However, cyanobacteria naturally direct a relatively small fraction of fixed carbon to lipid synthesis, primarily producing long-chain membrane-bound fatty acids. We investigated whether mitigating kinetic bottlenecks within the fatty acid biosynthesis (FAB) pathway could enhance flux to free fatty acid (FFA) production in Picosynechococcus sp. PCC 7002. Previous in vitro studies proposed that the FAB initiating enzyme FabH is the primary rate-controlling enzyme in PCC 7002. We hypothesized that enhancing fatty acid initiation could increase in vivo FFA production while shifting the kinetic bottleneck further downstream in the pathway. We enabled C8-FFA accumulation by knocking out the native acyl-acyl carrier protein synthetase gene (aas) and expressing the highly active Cuphea palustris-derived mutant thioesterase CpFatB1.2-M4-287 (CupTE), which selectively catalyzes C8 chain termination. We then expressed the diatom-derived Chaetoceros sp. GSL56 FabH ortholog (chKASIII), to enhance initiation, and a medium-chain selective E. coli ketosynthase ecFabF[I108F] (ecFabF∗), to enhance elongation. When expressed individually or in combination with ecFabF∗, chKASIII slowed growth and decreased net carbon fixation rates relative to the parental Δaas-CupTE strain. However, co-expression of chKASIII, ecFabF∗, and CupTE redirected a larger fraction of fixed carbon toward FFA production, increasing relative carbon flux to C8-FFA and decreasing the projected minimum selling price by four-fold. This work demonstrates how systems metabolic engineering can be applied to enhance C8-FFA production in cyanobacteria while highlighting the unpredictable physiological consequences of host metabolic burden.

09 BIOMASS FUELS↗
DOE OSTI2022

AmeriFlux FLUXNET-1F US-ARM ARM Southern Great Plains site- Lamont

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-ARM ARM Southern Great Plains site- Lamont. This is the FLUXNET version of the carbon flux data for the site US-ARM ARM Southern Great Plains site- Lamont produced by applying the standard ONEFlux (1F) software. Site Description - Central facility tower crop field (winter wheat, corn, soy, alfalfa). The site also has continuous measurements of precise mixing ratios of CO2, CH4, CO, N2O, and isotopic ratio of 13CO2. In addition to 4 m system, there are sonic anemometers at 25 m and 60 m.

Biraud, Sebastien↗
DOE OSTI2022

Underground sources of drinking water chemistry changes in response to potential CO 2 leakage

The purpose of this study was to quantify changes to underground sources of drinking water (USDW) quality in response to potential CO 2 leakage from geologic CO 2 sequestration (GCS) reservoirs. Here, the developed a framework of combined laboratory experiments and reactive transport simulations and used this framework to evaluate the Ogallala aquifer overlying the Farnsworth Unit (FWU), an active GCS site, as a case study. Using chemical reaction parameters obtained from laboratory experiments and numerical simulations, site-specific mechanisms of CO 2 -water-sediment interactions at the USDW aquifer were interpreted. Long-term risks of potential CO 2 leakage were then evaluated with field-scale numerical models using the regional hydrogeological characteristics and reaction parameters obtained from our experiments and simulations. Results suggest that carbonate mineral impurity and cation exchange are key mechanisms for interactions between CO 2 and the aquifer sediment. Additionally, for a large leakage rate of 0.1 % injection from one leaky well, the leakage plume might impact an area of 300 m in diameter and significantly affect the local water quality by changing pH and cation concentrations (e.g., Zn, Ba and Sr). After leakage ceases, the zone of impacted fluids would not migrate significantly in subsequent decades due to a low regional groundwater flowrate (for this case study). The relatively small area of impact might not be detected in a monitoring well given the broader spacing in a typical field scenario. Effective early leakage detection may require additional tools, e.g., borehole CO 2 movement, four-dimensional seismicity, CO 2 soil flux, samples from deeper aquifers, etc., to ensure effective leakage detection and long-term safety of GCS projects.

58 GEOSCIENCES↗
DOE OSTI2024

Evaluation of Additively Manufactured Monolithic SiC and SiC-Matrix Ceramic Matrix Composites for Concentrating Solar Receiver Applications - Fabrication and Testing of Receiver Design Feature Specimens in a Simulating Lab Test via Laser Heating

Increasing operating temperatures of solar receivers is paramount to the efficiency of concentrated solar thermal (CST) and solar power (CSP) systems. Owing to its high temperature stability combined with excellent thermal and optical properties, SiC has been the material of choice for application in high-temperature solar receivers. We report the results of our study of the effective heat transfer characteristics of several candidate SiC structure motifs, or feature geometries, which are fabricated via additive manufacturing. The SiC structure motifs studied include different permutations of three-dimensional periodic lattices and defined shapes. A solar-thermal simulating laboratory test setup is constructed using a 4kW CO2 laser system with beam shaping optics to apply concurrent radiative heating power on one face of 2”-diameter cylindrical feature specimens, representing the structure motifs of interest for receiver element design, while flowing through the sample as heat transfer fluid. Using the test setup, simulative test conditions representative of a concentrated solar flux of up to ~2000 suns could be achieved in the lab tests under varying air flow through the test structure. A simple 1D numerical analysis scheme is developed to extract an effective or compound heat transfer coefficient representative of the test structure under steady-state heat flow conditions. The test results and their use to guide the selection and optimization of SiC material and structure motifs for the receiver element design fabrication are discussed.

14 SOLAR ENERGY↗