SEARCH · Engineering Papers

Results for “coordinated optimization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,513 records · Page 84

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Optimized Substrate Positioning Enables Switches in the C–H Cleavage Site and Reaction Outcome in the Hydroxylation–Epoxidation Sequence Catalyzed by Hyoscyamine 6β-Hydroxylase

Hyoscyamine 6β-hydroxylase (H6H) is an Fe(II)- and 2-oxoglutarate-dependent (Fe/2OG) oxygenase that catalyzes the last two steps in the biosynthesis of scopolamine, a prolifically administered anti-nausea drug. After its namesake first reaction, H6H couples the newly installed C6-bonded oxygen to C7 to form the epoxide of scopolamine. Oxoiron(IV) (ferryl) intermediates initiate both reactions by cleaving C–H bonds, but it remains unclear how the enzyme switches target site and promotes (C6)O–C7 coupling in preference to C7 hydroxylation in the second step. In one possible epoxidation mechanism, the C6 oxygen would – analogously to mechanisms proposed for the Fe/2OG halogenases and, in the preceding paper, N-acetylnorloline synthase (LolO) – coordinate as alkoxide to the C7–H-cleaving ferryl intermediate to enable alkoxyl coupling to the ensuing C7 radical. Here we provide structural and kinetic evidence that H6H instead exploits the distinct spatial dependencies of competitive C–H-cleavage (C6 vs C7) and C–O-coupling (oxygen rebound vs cyclization) steps to promote the two-step sequence without substrate coordination or repositioning for the epoxidation step. Structural comparisons of ferryl-mimicking vanadyl complexes of wild-type H6H and a variant that preferentially hydroxylates C7 of 6-hydroxyhyoscyamine suggest that only a modest (~ 10°) shift in the Fe–O–H(C7) approach angle is sufficient to change the outcome. Finally, the observation that, in wild-type H6H, 2 H 2 O solvent also increases the C7-hydroxylation:epoxidation ratio by ~ 8-fold implies that the latter outcome requires cleavage of the alcohol O-H bond, which, unlike in the LolO oxacyclization, is not accomplished in advance of C–H cleavage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Geometry-dependent distributed polarizability models for the water molecule

Geometry-dependent distributed polarizability models have been constructed by fits to ab initio calculations at the coupled cluster level of theory with up to noniterative triple excitations in an augmented triple-zeta quality basis set for the water molecule in the field of a point charge. The investigated models include (i) charge-flow polarizabilities between chemically bonded atoms, (ii) isotropic or anisotropic dipolar polarizabilities on oxygen atom or on all atoms, and (iii) combinations of models (i) and (ii). For each model, the polarizability parameters have been optimized to reproduce the induction energy of a water molecule polarized by a point charge successively occupying a grid of points surrounding the molecule. The quality of the models is ascertained by examining their ability to reproduce these induction energies as well as the molecular dipolar and quadrupolar polarizabilities. The geometry dependence of the distributed polarizability models has been explored by changing bond lengths and HOH angle to generate 125 molecular structures (reduced to 75 symmetry-unique ones). For each considered model, the distributed polarizability components have been fitted as a function of the geometry by a Taylor expansion in monomer coordinate displacements up to the sum of powers equal to 4.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Multi-Agent Deep Reinforcement Learning for Realistic Distribution System Voltage Control Using PV Inverters

Over the last few decades, the deployment of distributed solar photovoltaic (PV) systems has increased consistently. High PV penetration could cause adverse effects on the grid, such as voltage violations. This paper proposes a new distributed soft actor-critic based multi-agent deep reinforcement learning (SAC-MADRL) control solution to minimize the PV real power curtailment while keeping the grid voltage in an acceptable range. New reward functions have been designed to coordinate different agents during the learning process, yielding improved convergence. Comparison results with other control methods on a real feeder in western Colorado U.S. with 80% penetration of PVs demonstrate that the proposed method has better capability of effectively regulating voltage while minimizing the PV real power curtailment.

distribution system↗

DOE OSTI2017

Models for the dynamics of dust-like matter in the self-gravity field: The method of hydrodynamic substitutions

Models for the dynamics of a dust-like medium in the self-gravity field are investigated. Solutions of the corresponding problems are constructed by the method of hydrodynamic substitutions generalizing the Cole–Hopf substitutions. The method is extended to multidimensional ideal and viscous fluid flows with cylindrical and spherical symmetries for which exact solutions are constructed. Solutions for the dynamics of self-gravitating dust with arbitrary initial distributions of both fluid density and velocity are constructed using special coordinate transformations. In particular, the problem of cosmological expansion is considered in terms of Newton’s gravity theory. Models of a one-dimensional viscous dust fluid flow and some problems of gas hydrodynamics are considered. Examples of exact solutions and their brief analysis are provided.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

AI Applications to Physics Experiments at Jefferson Lab

We survey how AI/ML is being deployed across Jefferson Lab's experimental and accelerator programs. In EPSCI, Hydra applies computer vision to automate real-time data-quality monitoring across all four experimental halls, replacing manual inspection of hundreds to thousands of histograms per shift. AIEC (AI Experiment Controls) uses ML to stabilize drift chamber gains and is now part of standard CEBAF production running, while AI Optimized Polarization (AIOP) targets autonomous control of polarized targets and photon beam angular alignment. In CASA, cavity fault classification models identify faulted cavities and trip types from waveform data with ~85% and ~78% agreement to labeled data, respectively, and are deployed in production; a separate effort applies LLMs and hybrid search to make the CEBAF operations logbook AI-ready. QCD-focused work includes transformer- and GAN-based generative models for particle-level event simulation, with distributed GAN training scaling studies on Polaris. Additional efforts span ML-on-FPGA for the EIC and a new Data Science Department coordinating anomaly detection, uncertainty quantification, and HPC-scalable ML lab-wide. Collectively, these projects illustrate AI's growing role in improving efficiency across JLab's nuclear physics mission.

Mei, Xinxin [Thomas Jefferson National Accelerator↗

NASA NTRS2006

4-Vinyl-1,3-Dioxolane-2-One as an Additive for Li-Ion Cells

Electrolyte additive 4-vinyl-1,3-dioxolane-2-one has been found to be promising for rechargeable lithium-ion electrochemical cells. This and other additives, along with advanced electrolytes comprising solutions of LiPF6 in various mixtures of carbonate solvents, have been investigated in a continuing effort to improve the performances of rechargeable lithium-ion electrochemical cells, especially at low temperatures. In contrast to work by other researchers who have investigated the use of this additive to improve the high-temperature resilience of Li-ion cells, the current work involves the incorporation of 4-vinyl-1,3-dioxolane-2-one into quaternary carbonate electrolyte mixtures, previously optimized for low-temperature applications, resulting in improved low-temperature performance. The benefit afforded by 4-vinyl-1,3- dioxolane-2-one can be better understood in the light of relevant information from a number of prior NASA Tech Briefs articles about electrolytes and additives for such cells. To recapitulate: The loss of performance with decreasing temperature is attributable largely to a decrease of ionic conductivity and the increase in viscosity of the electrolyte. What is needed to extend the lower limit of operating temperature is a stable electrolyte solution with relatively small lowtemperature viscosity, a large electric permittivity, adequate coordination behavior, and appropriate ranges of solubilities of liquid and salt constituents. Whether the anode is made of graphitic or non-graphitic carbon, a film on the surface of the anode acts as a solid/electrolyte interface (SEI), the nature of which is critical to low-temperature performance. Desirably, the surface film should exert a chemically protective (passivating) effect on both the anode and the electrolyte, yet should remain conductive to lithium ions to facilitate intercalation and de-intercalation of the ions into and out of the carbon during discharging and charging, respectively. The additives investigated previously include alkyl pyrocarbonates. Those additives help to improve low-temperature performances by giving rise to the formation of SEIs having desired properties. The formation of the SEIs is believed to be facilitated by products (e.g., CO2) of the decomposition of these additives. These decomposition products are believed to react to form Li2CO3-based films on the carbon electrodes. The present additive, 4-vinyl-1,3-dioxolane-2-one, also helps to improve lowtemperature performance by contributing to the formation of SEIs having desired properties, but probably in a different manner: It is believed that, as part of the decomposition process, the compound polymerizes on the surfaces of carbon electrodes.

Smart, Marshall↗

NASA NTRSDate not supplied

Model Predictive Control Study for an Electrified Turbofan Engine

A control sensitivity study is conducted for electrified aircraft propulsion concepts using a publicly available turbofan engine simulation. A Proportional-Integral (PI) controller serves as a baseline for comparison against an advanced Model Predictive Control (MPC) approach. MPC is used in the control scheme to allow for subsystem coordination of engine and electric power loads. The electrified turbofan engine includes electric machines to augment the amount of torque being applied to the spools of the engine during transient phases of operation. Power augmentation allows for energy to be stored or exhausted in a manner that increases the fuel efficiency of the engine, reducing aircraft emissions and noise. The study compares the fan speed, high pressure compressor stall margin, and low pressure compressor stall margin against the baseline PI control system, an MPC system without electric machines, and an MPC system with electric machines. The MPC system with electric machines is optimized using the sample time, control horizon and prediction horizon. The MPC design provided a stable damped response as desired and allowed operability margins to be maintained

Electrified Aircraft Propulsion↗

DOE OSTI2015

Multi-GPU implementation of a VMAT treatment plan optimization algorithm

Purpose: Volumetric modulated arc therapy (VMAT) optimization is a computationally challenging problem due to its large data size, high degrees of freedom, and many hardware constraints. High-performance graphics processing units (GPUs) have been used to speed up the computations. However, GPU’s relatively small memory size cannot handle cases with a large dose-deposition coefficient (DDC) matrix in cases of, e.g., those with a large target size, multiple targets, multiple arcs, and/or small beamlet size. The main purpose of this paper is to report an implementation of a column-generation-based VMAT algorithm, previously developed in the authors’ group, on a multi-GPU platform to solve the memory limitation problem. While the column-generation-based VMAT algorithm has been previously developed, the GPU implementation details have not been reported. Hence, another purpose is to present detailed techniques employed for GPU implementation. The authors also would like to utilize this particular problem as an example problem to study the feasibility of using a multi-GPU platform to solve large-scale problems in medical physics. Methods: The column-generation approach generates VMAT apertures sequentially by solving a pricing problem (PP) and a master problem (MP) iteratively. In the authors’ method, the sparse DDC matrix is first stored on a CPU in coordinate list format (COO). On the GPU side, this matrix is split into four submatrices according to beam angles, which are stored on four GPUs in compressed sparse row format. Computation of beamlet price, the first step in PP, is accomplished using multi-GPUs. A fast inter-GPU data transfer scheme is accomplished using peer-to-peer access. The remaining steps of PP and MP problems are implemented on CPU or a single GPU due to their modest problem scale and computational loads. Barzilai and Borwein algorithm with a subspace step scheme is adopted here to solve the MP problem. A head and neck (H and N) cancer case is then used to validate the authors’ method. The authors also compare their multi-GPU implementation with three different single GPU implementation strategies, i.e., truncating DDC matrix (S1), repeatedly transferring DDC matrix between CPU and GPU (S2), and porting computations involving DDC matrix to CPU (S3), in terms of both plan quality and computational efficiency. Two more H and N patient cases and three prostate cases are used to demonstrate the advantages of the authors’ method. Results: The authors’ multi-GPU implementation can finish the optimization process within ∼1 min for the H and N patient case. S1 leads to an inferior plan quality although its total time was 10 s shorter than the multi-GPU implementation due to the reduced matrix size. S2 and S3 yield the same plan quality as the multi-GPU implementation but take ∼4 and ∼6 min, respectively. High computational efficiency was consistently achieved for the other five patient cases tested, with VMAT plans of clinically acceptable quality obtained within 23–46 s. Conversely, to obtain clinically comparable or acceptable plans for all six of these VMAT cases that the authors have tested in this paper, the optimization time needed in a commercial TPS system on CPU was found to be in an order of several minutes. Conclusions: The results demonstrate that the multi-GPU implementation of the authors’ column-generation-based VMAT optimization can handle the large-scale VMAT optimization problem efficiently without sacrificing plan quality. The authors’ study may serve as an example to shed some light on other large-scale medical physics problems that require multi-GPU techniques.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

Producing glucose 6-phosphate from cellulosic biomass: Structural insights into levoglucosan bioconversion

The most abundant carbohydrate product of cellulosic biomass pyrolysis is the anhydrosugar levoglucosan (1,6-anhydro-β-d-glucopyranose), which can be converted to glucose 6-phosphate by levoglucosan kinase (LGK). In addition to the canonical kinase phosphotransfer reaction, the conversion requires cleavage of the 1,6-anhydro ring to allow ATP-dependent phosphorylation of the sugar O6 atom. Using x-ray crystallography, we show that LGK binds two magnesium ions in the active site that are additionally coordinated with the nucleotide and water molecules to result in ideal octahedral coordination. To further verify the metal binding sites, we co-crystallized LGK in the presence of manganese instead of magnesium and solved the structure de novo using the anomalous signal from four manganese atoms in the dimeric structure. The first metal is required for catalysis, whereas our work suggests that the second is either required or significantly promotes the catalytic rate. Although the enzyme binds its sugar substrate in a similar orientation to the structurally related 1,6-anhydro-N-acetylmuramic acid kinase (AnmK), it forms markedly fewer bonding interactions with the substrate. In this orientation, the sugar is in an optimal position to couple phosphorylation with ring cleavage. We also observed a second alternate binding orientation for levoglucosan, and in these structures, ADP was found to bind with lower affinity. These combined observations provide an explanation for the high Km of LGK for levoglucosan. Furthermore, greater knowledge of the factors that contribute to the catalytic efficiency of LGK can be used to improve applications of this enzyme for levoglucosan-derived biofuel production.

09 BIOMASS FUELS↗

DOE OSTI2026

A Survey on the Expanding Scope and Interdisciplinary Opportunities for Processing-in-Memory Techniques

Processing-in-Memory (PIM) is emerging as a practical path to overcome the limitations of traditional von Neumann architectures. At its core, PIM systems implement computing primitives such as logic operations and multiply-accumulate acceleration through compute-in-memory, near-memory processing, or hybrid designs. The role of memory cells varies widely across technologies, acting as inputs, outputs, or analog accumulators through bit-lines and sense amplifiers. This diversity creates trade-offs in precision, bandwidth, latency, and programmability, making it difficult to build a unified understanding on the progress of the field. In this survey, we organize recent advances of PIM into three areas. First, we discuss the progress on the architectural optimizations of PIM and its integration with both DRAM and emerging non-volatile memories. Second, we examine how PIM is being used to accelerate key computing domains, including generative AI workloads and high-performance kernels, along with new approaches. Third, we highlight the growing adoption of PIM in computational sciences, where it is being applied to solve interdisciplinary problems such as genome analysis, mRNA quantification, mass spectrometry, quantum circuit simulation, wave modeling, and secure computation. Finally, we synthesize the major challenges that continue to slow PIM adoption, including manufacturing constraints, power delivery, thermal reliability, data consistency, runtime and memory-management coordination, and the difficulty of building portable software abstractions without sacrificing commercial viability. This work provides an updated, structured perspective on PIM’s potential across computing and computational sciences and the barriers that must be solved for it to reach its full impact.

Asifuzzaman, Kazi [Oak Ridge National Laboratory (↗

DOE OSTI2025

Metal–Oxide Interface Sites Created Using Atomic Layer Deposition and Tested for CO Oxidation

The performance of catalysts made out of Pt supported on TiO 2 thin films grown on SBA-15 (a silica mesoporous material) by atomic layer deposition (ALD) was characterized systematically by combining in situ infrared absorption spectroscopy (IR) with other techniques including electron microscopy and adsorption−desorption isothermal measurements. The titania films in the resulting high-surface-area catalysts were evenly distributed throughout the inner surface of the SBA-15 mesopores, and their thickness could be controlled at a submonolayer level, with 3 to 4 TiO 2 ALD cycles needed for the complete coverage of the silica sites. The titania films could be deposited either before or after adding the metal (Pt), which was dispersed in the form of small nanoparticles (NPs) approximately 4−6 nm in diameter, in order to exert some control on the density and nature of the Pt/TiO 2 interface sites. One important lesson deriving from this work is that such an order of deposition leads to significantly different catalysts in spite of the fact that most of their structural properties are similar. If the Pt is deposited on the titania films, the resulting metal NPs are slightly smaller than those grown on silica and display CO adsorption sites with lower surface Pt coordination numbers. On the other hand, when TiO 2 is deposited on the Pt/SBA-15 starting material, some titania grows on the metal and partially blocks its surface while also creating new interface sites where CO binds more weakly and displays lower C−O stretching frequencies. In terms of catalytic performance, the results from in situ IR CO site titration and kinetic measurements combined suggest a mechanism where CO first adsorbs on Pt atop sites and then migrates to Pt/TiO 2 interface sites, where oxidation takes place. Both types of sites appear to be similar in all the catalysts tested, but catalytic performance could be optimized by tuning their surface densities. Maximum catalytic activity was obtained when the TiO2 films were deposited first and with TiO 2 coverages of at least half a monolayer, that is, after at least 2 ALD cycles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Rational solvent molecule tuning for high-performance lithium metal battery electrolytes

Electrolyte engineering improved cycling of Li metal batteries and anode-free cells at low current densities; however, high-rate capability and tuning of ionic conduction in electrolytes are desirable yet less-studied. In this study, we design and synthesize a family of fluorinated-1,2-diethoxyethanes as electrolyte solvents. The position and amount of F atoms functionalized on 1,2-diethoxyethane were found to greatly affect electrolyte performance. Partially fluorinated, locally polar –CHF 2 is identified as the optimal group rather than fully fluorinated –CF 3 in common designs. Paired with 1.2 M lithium bis(fluorosulfonyl)imide, these developed single-salt-single-solvent electrolytes simultaneously enable high conductivity, low and stable overpotential, >99.5% Li||Cu half-cell efficiency (up to 99.9%, ±0.1% fluctuation) and fast activation (Li efficiency >99.3% within two cycles). Combined with high-voltage stability, these electrolytes achieve roughly 270 cycles in 50-μm-thin Li||high-loading-NMC811 full batteries and >140 cycles in fast-cycling Cu||microparticle-LiFePO 4 industrial pouch cells under realistic testing conditions. The correlation of Li + –solvent coordination, solvation environments and battery performance is investigated to understand structure–property relationships.

25 ENERGY STORAGE↗

NASA NTRS1978

A coordinated X-ray, optical, and microwave study of the flare star Proxima Centauri

Results are reported for a three-day coordinated observing program to monitor the flare star Proxima Centauri in the X-ray, optical, and radio spectrum. During this interval 30 optical flares and 12 possible radio bursts were observed. The SAS 3 X-ray satellite made no X-ray detections. An upper limit of 0.08 on the X-ray/optical luminosity ratio is derived for the brightest optical flare. The most sensitive of the radio telescopes failed to detect 6-cm emission during one major and three minor optical flares, and on this basis an upper limit on the flare radio emission (1 hundred-thousandth of the optimal luminosity) is derived.

Haisch, B. M.↗

DOE OSTI2024

Probing Basal and Prismatic Planes of Graphitic Materials for Metal Single Atom and Subnanometer Cluster Stabilization

Abstract Supported metal single atom catalysis is a dynamic research area in catalysis science combining the advantages of homogeneous and heterogeneous catalysis. Understanding the interactions between metal single atoms and the support constitutes a challenge facing the development of such catalysts, since these interactions are essential in optimizing the catalytic performance. For conventional carbon supports, two types of surfaces can contribute to single atom stabilization: the basal planes and the prismatic surface; both of which can be decorated by defects and surface oxygen groups. To date, most studies on carbon‐supported single atom catalysts focused on nitrogen‐doped carbons, which, unlike classic carbon materials, have a fairly well‐defined chemical environment. Herein we report the synthesis, characterization and modeling of rhodium single atom catalysts supported on carbon materials presenting distinct concentrations of surface oxygen groups and basal/prismatic surface area. The influence of these parameters on the speciation of the Rh species, their coordination and ultimately on their catalytic performance in hydrogenation and hydroformylation reactions is analyzed. The results obtained show that catalysis itself is an interesting tool for the fine characterization of these materials, for which the detection of small quantities of metal clusters remains a challenge, even when combining several cutting‐edge analytical methods.

Vidal, Mathieu↗

DOE OSTI2025

Synthesis of UO 2 nanoparticles via coulometric titration: influence of electrolytes and ligands on synthesized particle properties

Here, this study investigates the controlled synthesis of uranium oxide nanoparticles (UO 2 NPs) via coulometric titration under ambient conditions, focusing on the impact of electrolyte anions, salt concentrations, and strongly complexing ligands on particle formation and properties. Using dilute solutions of mineral acids (HNO 3 , HCl, and HClO 4 ), we demonstrate that the choice of electrolyte anion affects particle size, polydispersity, and surface charge. Particle characterization using dynamic light scattering and transmission electron microscopy shows that particles synthesized in perchlorate are largest and show a high degree of polydispersity. In comparison, particles synthesized in chloride are smaller and more uniform in size. Synthesis in nitrate yields a wide variety of particle sizes, with the major fraction (>65 %) having a smaller size than that obtained in the other two electrolytes. The presence of more strongly coordinating ligands such as sulfate and acetate modulate hydrolysis and condensation reactions, with sulfate suppressing nanoparticle formation across a wide concentration range ([SO 4 2- ] > 5 mM) and acetate enabling stable colloidal suspensions at [HOAc] ≤ 0.1 M. Synthesis in concentrated electrolytes (2 M NaNO 3 , NaCl, or NaClO 4 ) accelerates reaction kinetics but introduces challenges, as particles showed increased polydispersity and aggregation, and were more prone to oxidation. Electrolyte effects on actinide oxide nanoparticle formation are discussed and a short comparison to established nanoparticle syntheses is drawn. This work underscores the importance of tailoring synthesis parameters to achieve desired nanoparticle properties, providing valuable insights for optimizing UO 2 NP production for various applications.

Actinide hydrolysis↗

DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗

DOE OSTI2022

Consistent responses of vegetation gas exchange to elevated atmospheric CO 2 emerge from heuristic and optimization models

Elevated atmospheric CO 2 concentration is expected to increase leaf CO 2 assimilation rates, thus promoting plant growth and increasing leaf area. It also decreases stomatal conductance, allowing water savings, which have been hypothesized to drive large-scale greening, in particular in arid and semiarid climates. However, the increase in leaf area could reduce the benefits of elevated CO 2 concentration through soil water depletion. The net effect of elevated CO 2 on leaf- and canopy-level gas exchange remains uncertain. To address this question, we compare the outcomes of a heuristic model based on the Partitioning of Equilibrium Transpiration and Assimilation (PETA) hypothesis and three model variants based on stomatal optimization theory. Predicted relative changes in leaf- and canopy-level gas exchange rates are used as a metric of plant responses to changes in atmospheric CO 2 concentration. Both model approaches predict reductions in leaf-level transpiration rate due to decreased stomatal conductance under elevated CO 2 , but negligible (PETA) or no (optimization) changes in canopy-level transpiration due to the compensatory effect of increased leaf area. Leaf- and canopy-level CO 2 assimilation is predicted to increase, with an amplification of the CO 2 fertilization effect at the canopy level due to the enhanced leaf area. The expected increase in vapour pressure deficit (VPD) under warmer conditions is generally predicted to decrease the sensitivity of gas exchange to atmospheric CO 2 concentration in both models. The consistent predictions by different models that canopy-level transpiration varies little under elevated CO 2 due to combined stomatal conductance reduction and leaf area increase highlight the coordination of physiological and morphological characteristics in vegetation to maximize resource use (here water) under altered climatic conditions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

Towards a universal model for carbon dioxide uptake by plants

Gross primary production (GPP) - the uptake of carbon dioxide (CO 2 ) by leaves, and its conversion to sugars by photosynthesis - is the basis for life on land. Earth System Models (ESMs) incorporating the interactions of land ecosystems and climate are used to predict the future of the terrestrial sink for anthropogenic CO 2 . ESMs require accurate representation of GPP. However, current ESMs disagree on how GPP responds to environmental variations, suggesting a need for a more robust theoretical framework for modelling. Here in this work, we focus on a key quantity for GPP, the ratio of leaf internal to external CO 2 (χ). χ is tightly regulated and depends on environmental conditions, but is represented empirically and incompletely in today's models. We show that a simple evolutionary optimality hypothesis predicts specific quantitative dependencies of χ on temperature, vapour pressure deficit and elevation; and that these same dependencies emerge from an independent analysis of empirical χ values, derived from a worldwide dataset of >3,500 leaf stable carbon isotope measurements. A single global equation embodying these relationships then unifies the empirical light-use efficiency model with the standard model of C 3 photosynthesis, and successfully predicts GPP measured at eddy-covariance flux sites. This success is notable given the equation's simplicity and broad applicability across biomes and plant functional types. Finally, it provides a theoretical underpinning for the analysis of plant functional coordination across species and emergent properties of ecosystems, and a potential basis for the reformulation of the controls of GPP in next-generation ESMs.

54 ENVIRONMENTAL SCIENCES↗