Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reproducible science”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Challenges in synthesis of heterostructures

Developing the ability to reproducibly prepare materials drives the advancement of human technology. Trends in historical advances and extrapolation of current technologies suggest that the next logical step is developing the ability to control structure and composition at near atomic scales in device manufacturing to take advantage of emergent properties discovered in 2D materials and heterostructures. This article reviews a subset of synthetic approaches used to prepare heterostructures, outlining key reaction steps, and how experimental parameters are used to control which product forms. The positive attributes of each synthesis approach are discussed along with some of their limitations. The elemental reactant (MER) approach is discussed in more depth.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Self-dual ν = 1 bosonic quantum Hall state in mixed-dimensional QED

Here, we consider a (2 + 1)-dimensional Wilson-Fisher boson coupled to a (3 + 1)-dimensional U(1) gauge field. This theory possesses a strong-weak duality in terms of the coupling constant e and is self-dual at a particular value of $\mathcal{e}$. We derive exact relations between transport coefficients for a $\mathcal{ν}$ = 1 quantum Hall state at the self-dual point. Using boson-fermion duality, we map the $\mathcal{ν}$ = 1 bosonic quantum Hall state to a Fermi sea of the dual fermion and observe that the exact relationships between transport coefficients at the bosonic self-dual point are reproduced by a simple random-phase approximation (RPA), coupled with a Drude formula, in the fermionic theory. We explain this success of the RPA by pointing out a cancellation of a parity-breaking term in the fermion theory which occurs only at the self-dual point, resulting in the fermion self-dual theory explored previously. In addition, we argue that the equivalence of the self-dual structure can be understood in terms of electromagnetic duality or modular invariance, and these features are not inherited by the nonrelativistic cousins of the present model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Hicks Pries, Sulman, et al 2018. In situ incubation of 13C-labeled litter in three soil pits at Blodgett Forest, CA from 2013 to 2016

Even though over half of the world's soil organic carbon (SOC) is stored in subsoils (>20 cm deep), and the old ages of subsoil OC indicate its cycling differs from surface SOC, there are few studies examining in situ decomposition processes in deep soils. The purpose of this dataset is to elucidate these processes. In this study we added 13C-labeled fine roots to 15, 55, and 95 cm depths of a well-characterized coniferous forest Alfisol and monitored the amount of root-derived C remaining over 6, 12, and 30 months. We recovered the root-derived C in microbial phospholipid fatty acids (PLFAs) after 6 months and in coarse (>2 mm) particulate, fine (<2 mm) particulate, and dense, mineral-associated pools after 6, 12, and 30 months. Overall, root decomposition in the first 6 months was similar among all depths but significantly diverged at 30 months with faster decomposition at 15 cm than at 95 cm. There were more fungal and Gram negative-associated PLFAs at 15 cm than at 95 cm, and 13C analysis revealed those microbial groups preferred the added root carbon to native SOC. Mineral-associations were not the cause of slower decomposition at depth because similar amounts of applied root C was recovered in the dense fraction at all depths. The largest difference among depths was in the amount of root C recovered in the coarse particulate fraction, which was greater at 95 cm (50%) than at 15 cm (15%). Slower decomposition of the particulate pool at depth likely contributed to the increase in C:N ratios and depletion of δ13C values below 60 cm depth in our soil profiles. Simulations of these soils using the CORPSE model, which incorporates microbial priming effects and mineral stabilization of SOC, reproduced patterns of particulate and mineral-associated SOC over both time and depth and suggested that a lack of priming by root exudates at depth could account for the slower decomposition rate of particulate root material. Decomposition of deep particulate SOC may increase if root exudation or dissolved OC transport to depth increases.This dataset includes (1) characterization of the soil pits prior to litter addition, (2) characterization of how the litter changed after 6, 12, and 30 months, and (3) the CORPSE model structure, parameterization, and output.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Synthesis of wafer-scale SWCNT forests with remarkably invariant structural properties in a bulk-diffusion-controlled kinetic regime

Robust synthesis of vertically aligned carbon nanotubes (VACNT) at large scale is required to accelerate deployment of numerous cutting-edge devices to emerging commercial applications. To address this need, this work demonstrates that the structural characteristics of single-walled CNTs (SWCNTs) produced in a growth regime dominated by bulk diffusion of the gaseous carbon precursor are remarkably invariant over a broad range of process conditions (precursor concentration increased up to 30-fold, catalyst substrate area from 1 cm 2 to 180 cm 2 , growth pressure from 20 to 790 mbar, and gas flowrates up to 8-fold). A simple growth kinetics model is derived that accounts for both reactant bulk diffusion and competing byproduct formation in the presence of excess hydrogen, which quantitatively reproduces the experimental data in the entire isothermal parameter space. The model enables critical predictions for process scale-up optimization including a potential 6-fold increase in CNT production rate, a ~90% carbon conversion efficiency in select conditions, and elimination of reaction rate decay observed at high pressures by using a hydrogen-free growth environment. This model-guided opportunity to substantially improve synthesis throughput and efficiency, in conjunction with the remarkably invariant CNT properties, is appealing for reliable VACNT device fabrication at both small and industrial scales.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Allosteric control of olefin isomerization kinetics via remote metal binding and its mechanochemical analysis

Abstract Allosteric control of reaction thermodynamics is well understood, but the mechanisms by which changes in local geometries of receptor sites lower activation reaction barriers in electronically uncoupled, remote reaction moieties remain relatively unexplored. Here we report a molecular scaffold in which the rate of thermal E-to-Z isomerization of an alkene increases by a factor of as much as 10 4 in response to fast binding of a metal ion to a remote receptor site. A mechanochemical model of the olefin coupled to a compressive harmonic spring reproduces the observed acceleration quantitatively, adding the studied isomerization to the very few reactions demonstrated to be sensitive to extrinsic compressive force. The work validates experimentally the generalization of mechanochemical kinetics to compressive loads and demonstrates that the formalism of force-coupled reactivity offers a productive framework for the quantitative analysis of the molecular basis of allosteric control of reaction kinetics. Important differences in the effects of compressive vs. tensile force on the kinetic stabilities of molecules are discussed.

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

polishCLR: A Nextflow Workflow for Polishing PacBio CLR Genome Assemblies

Long-read sequencing has revolutionized genome assembly, yielding highly contiguous, chromosome-level contigs. However, assemblies from some third generation long read technologies, such as Pacific Biosciences (PacBio) continuous long reads (CLR), have a high error rate. Such errors can be corrected with short reads through a process called polishing. Although best practices for polishing non-model de novo genome assemblies were recently described by the Vertebrate Genome Project (VGP) Assembly community, there is a need for a publicly available, reproducible workflow that can be easily implemented and run on a conventional high performance computing environment. Here, we describe polishCLR (https://github.com/isugifNF/polishCLR), a reproducible Nextflow workflow that implements best practices for polishing assemblies made from CLR data. PolishCLR can be initiated from several input options that extend best practices to suboptimal cases. It also provides re-entry points throughout several key processes, including identifying duplicate haplotypes in purge_dups, allowing a break for scaffolding if data are available, and throughout multiple rounds of polishing and evaluation with Arrow and FreeBayes. PolishCLR is containerized and publicly available for the greater assembly community as a tool to complete assemblies from existing, error-prone long-read data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Magnetic Damping in Polycrystalline Thin-Film Fe-V Alloys

Here, we report on the magnetic damping properties of polycrystalline Fe-V alloy thin films that are deposited at room temperature. By varying the concentration of V in the alloy, the saturation magnetization can be adjusted from that of Fe to near zero.We show that exceptionally low values of the damping parameter can be maintained over the majority of this range, with a minimum damping at approximately 15%–20% V concentration. Such a minimum is qualitatively reproduced with ab initio calculations of the damping parameter, although at a concentration closer to 10% V. The measured intrinsic damping has a minimum value of (1.53±0.08)×10 –3 , which is approximately a factor of 3 higher than our calculated value of 0.48×10 –3 . From first-principles theory, we outline the factors that are mainly responsible for the trend of the damping parameter in these alloys. In particular, the band structure and resulting damping mechanism is shown to change at V concentrations greater than approximately 35% V content.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evolving Dark Sector and the Dark Dimension Scenario

String theory naturally leads to the expectation that dark energy is not stable, and may be evolving as captured by the Swampland de Sitter conjectures. Moreover, motivated by the distance conjecture, a unification of dark sector has been proposed, where the smallness of dark energy leads to one extra dimension of micron size with dark matter being the Kaluza--Klein graviton excitations in this extra dimension. We consider the natural possibility that the radius of the dark dimension varies as the dark energy decreases, leading to the variation of the dark matter mass. This correlates the decrease of the dark energy with the variation of the dark matter mass as they depend on the variations of a scalar field $ϕ$ controlling the radius of the extra dimension. A simple realization of this idea for small range of $ϕ$ is adequately captured by choosing a potential which is locally of the form $V=V_0\ {\rm exp}(-cϕ)$ and dark matter mass $m_{\rm DM}=m_0\ {\rm exp}(-c' ϕ)$ where the sign of $ϕ$ is chosen such that $c'\geq 0$ while we have two choices for the sign of $c$ depending on whether the dark dimension expands or shrinks when the dark energy dominates. We find excellent agreement with recent experimental data from DESI DR2 combined with SN measurements (from DES, Union3 or Pantheon+) and reproduces the same significance as CPL parametrization with the added benefit of providing a natural explanation for the apparent phantom behavior ($w<-1$) reported by DESI and DES based on a physical model. DESI and SN datasets independently favor non-zero values of $c'$ and $c$, respectively, both lying within the expected $\mathcal{O}(1)$ range suggested by the Swampland criteria. Moreover, our best fit value $c'\simeq 0.05 \pm 0.01$ is remarkably consistent with the experimental upper bound of $c'\lesssim 0.2$ demanded by the lack of detection of fifth force in the dark sector.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗
DOE OSTI2023

The neutron number probability distribution in coupled lumped assemblies

Here, the validity of the gamma distribution in describing the neutron number probability distribution function for both isolated and coupled multiplying assemblies when constrained to reproduce the true mean and variance is investigated in lumped geometry by numerical comparison with kinetic Monte Carlo simulations. The mean and variance are obtained from numerical solution of moment equations constructed from the relevant forward Master equation with assembly coupling coefficients obtained from a view factor method. Numerical results for a two-group, two coupled assemblies model, with static and dynamic reactivity insertion, show that except for subcritical assemblies, the gamma distribution well-approximates the number distribution. Differences in the fast and thermal neutron population shapes are explained in terms of effective source strengths due to downscatter and coupling.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

NanoPSD: A software for automatic detection of Nano-Particle Shape Distribution in electron microscopy images

Accurate quantification of the size and morphology of nanoparticles from electron microscopy (EM) images is essential to understand growth mechanisms, surface reactivity, and functional behavior in nanoscale materials. Manual analysis remains slow, subjective, and difficult to reproduce in large datasets. We introduce NanoPSD (Nano-Particle Shape Distribution), an open-source and fully automated framework for quantitative particle detection and morphology analysis from EM images. NanoPSD integrates adaptive contrast enhancement, polarity-agnostic scale-bar detection, Optical Character Recognition (OCR)-based calibration, and classical segmentation via Otsu thresholding with morphological refinement. Particle contours are used to extract geometric descriptors, including equivalent circular diameter, aspect ratio, circularity, and solidity, enabling automated classification into spherical, rod-like, and aggregate morphologies. The framework supports both single-image and batch processing, generating publication-quality visualizations, LaTeX-ready tables, and structured comma-separated values (CSV) datasets. As a demonstration, we applied NanoPSD to plasma-synthesized nanoparticle samples diagnosed via transmission electron microscopy (TEM). The code produced statistically robust size and morphology distributions spanning a few to tens of nanometers with minimal user supervision. The pipeline demonstrates high reproducibility and scalability, processing large image collections with consistent calibration and output formatting. Its modular design enables seamless integration of future deep-learning-based segmentation models, providing a pathway toward intelligent, data-driven electron microscopy analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Many-Body Potential for Simulating Self-Assembly of Polymer-Grafted Nanoparticles in a Polymer Matrix

Many-body interactions between polymer-grafted nanoparticles (NPs) play a key role in promoting their assembly into low-dimensional structures within polymer melts, even when the particles are spherical and isotropically grafted. However, capturing such interactions in simulations of NP assembly is very challenging because explicit modeling of the polymer grafts and melt chains is highly computationally expensive, even using coarse-grained models. Here, we develop a many-body potential for describing the effective interactions between spherical polymer-grafted NPs in a polymer matrix through a machine-learning approach. The approach involves using permutationally invariant polynomials to fit two- and three-body interactions derived from the potential of mean force calculations. The potential developed here reduces the computational cost by several orders of magnitude, thereby, allowing us to explore assembly behavior over large length and time scales. We show that the potential not only reproduces previously known assembled phases such as 1D strings and 2D hexagonal sheets, which generally cannot be achieved using isotropic two-body potentials, but can also help discover interesting phases such as networks, clusters, and gels. We demonstrate how each of these assembly morphologies intrinsically arises from a competition between two- and three-body interactions. Our approach for deriving many-body effective potentials can be readily extended to other colloidal systems, enabling researchers to make accurate predictions of their behavior and dissect the role of individual interaction energy terms of the overall potential in the observed behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Polysulfide entrapment and retardation in gel electrolyte Li–S batteries: experiments and modeling

Gelling the electrolyte in lithium–sulfur (Li–S) batteries provides advantages in terms of safety and mitigating polysulfide shuttling leading to reduced side reactions and dendrite growth on the anode side. In the current study, organic and ceramic crosslinkers have been used to improve the entrapment of polysulfides and prevent them from shuttling to the anode side. However, the performance of gel electrolytes is still far from being on par with liquid electrolytes in terms of raw capacity. The discharge curves exhibit an extended first plateau for gel electrolyte batteries compared to liquid electrolytes, showing that polysulfides are better trapped in the gel system. Yet, in the conversion of soluble polysulfides into insoluble polysulfides, the second plateau was delayed. Qualitative and quantitative analysis of the polysulfide diffusion in the liquid electrolyte and gel electrolyte systems were also conducted to predict the relative diffusivity ratio between the liquid electrolyte and gel electrolyte for simulation. To probe what happens at the end of the discharge, where the mass transport in the gel electrolyte appears to be slow, a modified discharge procedure was used where the cell is first paused after discharge and then discharged a second time before recharging normally. Coin cell tests for two gel electrolytes in addition to liquid electrolyte show up to 7% capacity recovery after the pause, which indicates that the diffusion of species during the pause is responsible for the failure recovery. Additionally, the recovery is higher for gel electrolytes compared to liquid, for higher C-rates, and for longer pause times. To better understand this behavior, a Li–S numerical model with mass transport limited reactions was used to examine the polysulfide diffusion coefficients expected in different electrolyte systems. The model is able to reproduce the trends seen in experiments and yields a higher recovery for smaller Li 2 S 4 diffusion coefficients, suggesting that insufficient Li 2 S 4 mass transport is responsible for failure at the end of discharge. By identifying the effect that gel electrolytes have on the state of the battery during discharge using both experimental and simulation methods, this work is a step towards leveraging the advantages of gel electrolytes and mitigating their weaknesses.

25 ENERGY STORAGE↗
DOE OSTI2021

MetFish: a Metabolomics Pipeline for Studying Microbial Communities in Chemically Extreme Environments

Metabolites have essential roles in microbial communities, including as mediators of nutrient and energy exchange, cell-to-cell communication, and antibiosis. However, detecting and quantifying metabolites and other chemicals in samples having extremes in salt or mineral content using liquid chromatography-mass spectrometry (LC-MS)-based methods remains a significant challenge. Here, we report a facile method based on in situ chemical derivatization followed by extraction for analysis of metabolites and other chemicals in hypersaline samples, enabling for the first time direct LC-MS-based exometabolomics analysis in sample matrices containing up to 2 M total dissolved salts. The method, MetFish, is applicable to molecules containing amine, carboxylic acid, carbonyl, or hydroxyl functional groups, and it can be integrated into either targeted or untargeted analysis pipelines. In targeted analyses, MetFish provided limits of quantification as low as 1 nM, broad linear dynamic ranges (up to 5 to 6 orders of magnitude) with excellent linearity, and low median interday reproducibility (e.g., 2.6%). MetFish was successfully applied in targeted and untargeted exometabolomics analyses of microbial consortia, quantifying amino acid dynamics in the exometabolome during community succession; in situ in a native prairie soil, whose exometabolome was isolated using a hypersaline extraction; and in input and produced fluids from a hydraulically fractured well, identifying dramatic changes in the exometabolome over time in the well.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Towards Automated and High-Throughput Quantitative Sizing and Isotopic Analysis of Nanoparticles via Single Particle-ICP-TOF-MS

The work described herein assesses the ability to characterize gold nanoparticles (Au NPs) of 50 and 100 nm, as well as 60 nm silver shelled gold core nanospheres (Au/Ag NPs), for their mass, respective size, and isotopic composition in an automated and unattended fashion. Here, an innovative autosampler was employed to mix and transport the blanks, standards, and samples into a high-efficiency single particle (SP) introduction system for subsequent analysis by inductively coupled plasma–time of flight–mass spectrometry (ICP-TOF-MS). Optimized NP transport efficiency into the ICP-TOF-MS was determined to be >80%. This combination, SP-ICP-TOF-MS, allowed for high-throughput sample analysis. Specifically, 50 total samples (including blanks/standards) were analyzed over 8 h, to provide an accurate characterization of the NPs. This methodology was implemented over the course of 5 days to assess its long-term reproducibility. Impressively, the in-run and day-to-day variation of sample transport is assessed to be 3.54 and 9.52% relative standard deviation (%RSD), respectively. The determination of Au NP size and concentration was of <5% relative difference from the certified values over these time periods. Isotopic characterization of the 107 Ag/ 109 Ag particles (n = 132,630) over the course of the measurements was determined to be 1.0788 ± 0.0030 with high accuracy (0.23% relative difference) when compared to the multi-collector–ICP-MS determination.

36 MATERIALS SCIENCE↗
NASA NTRS2014

Extreme Environment Simulation - Current and New Capabilities to Simulate Venus and Other Planetary Bodies

Science, technology, and planetary mission communities have a growing interest in components and systems that are capable of working in extreme (high) temperature and pressure conditions. Terrestrial applications range from scientific research, aerospace, defense, automotive systems, energy storage and power distribution, deep mining and others. As the target environments get increasingly extreme, capabilities to develop and test the sensors and systems designed to operate in such environments will be required. An application of particular importance to the planetary science community is the ability for a robotic lander to survive on the Venus surface where pressures are nearly 100 times that of Earth and temperatures approach 500C. The scientific importance and relevance of Venus missions are stated in the current Planetary Decadal Survey. Further, several missions to Venus were proposed in the most recent Discovery call. Despite this interest, the ability to accurately simulate Venus conditions at a scale that can test and validate instruments and spacecraft systems and accurately simulate the Venus atmosphere has been lacking. This paper discusses and compares the capabilities that are known to exist within and outside the United States to simulate the extreme environmental conditions found in terrestrial or planetary surfaces including the Venus atmosphere and surface. The paper then focuses on discussing the recent additional capability found in the NASA Glenn Extreme Environment Rig (GEER). The GEER, located at the NASA Glenn Research Center in Cleveland, Ohio, is designed to simulate not only the temperature and pressure extremes described, but can also accurately reproduce the atmospheric compositions of bodies in the solar system including those with acidic and hazardous elements. GEER capabilities and characteristics are described along with operational considerations relevant to potential users. The paper presents initial operating results and concludes with a sampling of investigations or tests that have been requested or expected.

Extreme Environment↗
NASA NTRS2015

2015 Space Radiation Standing Review Panel

The 2015 Space Radiation Standing Review Panel (from here on referred to as the SRP) met for a site visit in Houston, TX on December 8 - 9, 2015. The SRP met with representatives from the Space Radiation Element and members of the Human Research Program (HRP) to review the updated research plan for the Risk of Radiation Carcinogenesis Cancer Risk. The SRP also reviewed the newly revised Evidence Reports for the Risk of Acute Radiation Syndromes Due to Solar Particle Events (SPEs) (Acute Risk), the Risk of Acute (In-flight) and Late Central Nervous System Effects from Radiation Exposure (CNS Risk), and the Risk of Cardiovascular Disease and Other Degenerative Tissue Effects from Radiation (Degen Risk), as well as a status update on these Risks. The SRP would like to commend Dr. Simonsen, Dr. Huff, Dr. Nelson, and Dr. Patel for their detailed presentations. The Space Radiation Element did a great job presenting a very large volume of material. The SRP considers it to be a strong program that is well-organized, well-coordinated and generates valuable data. The SRP commended the tissue sharing protocols, working groups, systems biology analysis, and standardization of models. In several of the discussed areas the SRP suggested improvements of the research plans in the future. These include the following: It is important that the team has expanded efforts examining immunology and inflammation as important components of the space radiation biological response. This is an overarching and important focus that is likely to apply to all aspects of the program including acute, CVD, CNS, cancer and others. Given that the area of immunology/inflammation is highly complex (and especially so as it relates to radiation), it warrants the expansion of investigators expertise in immunology and inflammation to work with the individual research projects and also the NASA Specialized Center of Research (NSCORs). Historical data on radiation injury to be entered into the Watson “big data” study must be used with caution. The general scientific issues of reproducibility, details of experimental methods and data analysis from preclinical and basic research laboratories have been raised broadly over the last few years (not specific to this work) and indicate that caution must be applied in the ways these data are used. This pertains to preclinical data and also to phase 3 clinical trials in radiation oncology and medical oncology. Of course, appropriate use and analysis of these “big-data” sets also offer the potential of pinpointing limitations and extracting remaining useful information. Emphasis should be placed on the latter possibility. A key target is risk reduction from radiation exposure. Progress of the entire space program, now moving towards the Mars mission, requires timely answers to key components of human risk, which are known to be complex. Periodic review of progress should be conducted with additional resources directed into achieving critical milestones. Turning the long red bars to yellow and green (or for some risks such as CNS possibly to grey) must be high priority. That such progress will require new science and not engineering means that it should be viewed in a knowledge-based light. The technology-based aspects of engineering issues are certainly as important, however, science and knowledge-based problems are solved in a different way than engineering. Timelines for engineering are more predictable, while for science, progress can be methodical with occasional major incremental findings that can rapidly change the rate of progress. As opportunities for rapid incremental changes arise, periodic enhancement of investment is strongly recommended to enable such new knowledge to be quickly and efficiently exploited. Collaborations and linkages with National Institute of Allergy and Infectious Diseases (NIAID), the Biomedical Advanced Research and Development Authority (BARDA) and the Department of Defense (DoD) are in place and more are encouraged, where possible, with the radiation injury and medical countermeasure studies. This could include utilizing some of their animal model testing contracts to facilitate obtaining results using common platforms. Such approach will facilitate the comparison of results among laboratories, and will facilitate and accelerate the development of medical countermeasures. It is particularly noteworthy that the NASA Space Radiation Element is reaching out to the Multidisciplinary European Low Dose Initiative (MELODI) platform coordinating low dose radiation risk research, and to other international agencies that are studying low dose radiation effects in an effort to fill the void generated by the cancelation of the Department of Energy (DOE) low dose radiation program. While NASA is working actively with NIAID and BARDA to integrate their relevant findings of radiation mitigator investigations to NASA programs, the committee notes its disappointment that the United States currently lacks a dedicated low dose radiation program with clear mechanistic orientation and aimed at the quantification and mitigation of human radiation risk on Earth. This void gives to the NASA Space Radiation Program Element special societal value, but also makes its overall design more challenging.

Steinberg, Susan↗
DOE OSTI2021

Quantifying carbon budget, crop yields and their responses to environmental variability using the ecosys model for U.S. Midwestern agroecosystems

As one of the major agricultural production areas in the world, the United States (U.S.) Midwest plays a vital role in the global food supply and agricultural ecosystem services. Although significant efforts have been made in modeling the carbon cycle dynamics over this area, large uncertainty still exists in the previous simulations in terms of reproducing individual components of the carbon cycle and their responses to environmental variability. Here we evaluated the performance of an advanced agroecosystem model, ecosys, in simulating carbon budgets over the U.S. Midwest, considering both the magnitude of carbon flux/yield and its response to environmental (climate and soil) variability. We conducted model simulations and evaluations at 7 cropland eddy-covariance sites as well as over 293 counties of Illinois, Indiana, and Iowa in the U.S. Midwest. The site-level simulations showed that ecosys captured both the magnitude and seasonal patterns of carbon fluxes (i.e., net ecosystem carbon exchange (NEE), ecosystem gross primary production (GPP), and ecosystem respiration (Reco)), leaf area index (LAI), and dynamic plant carbon allocation processes, with R2 equal to 0.92, 0.87, 0.87, and 0.78 for GPP, NEE, Reco, and LAI, respectively across all the sites compared with the observations. For regional scale simulations, ecosys reproduced the spatial distribution and interannual variability of corn and soybean yields with the constraints of observed yields and a new remotely sensed GPP product, with R2 of multi-year averaged simulated and observed yield equal 0.83 and 0.80 for corn and soybean, respectively. The simulated responses of carbon cycle dynamics to environmental variability were consistent with that from the empirical observations at both site and regional scales. Furthermore, our results demonstrated the applicability of ecosys in simulating the carbon cycle and soil carbon dynamics of the U.S. Midwestern agroecosystems under different climate and soil conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Insights into storm direction effect on flood response

In this study, we investigate the directional influence of storm movement on catchment flood peak response using the synthetic circular basin. Due to the complexity in defining storm movements that require meteorological modeling, we adopt a novel approach of combining the basin rotation method (BRM) with a circular basin construct. A systematic basin rotation approach provides a proxy for the storm direction relative to the river network topology to study its hydrologic response. Using Stage-IV rainfall with 4-km by 1-hour resolution for 16 years period from 2003 to 2018 in Iowa, U.S., we systematically analyze consequent flood peak response due to storm directions (basin rotations) using a distributed hydrologic model called the Hillslope-Link Model (HLM). HLM has demonstrated the ability to reproduce observed streamflow for real river basins studied in Iowa. The hydrologic simulation results show that the BRM can quantify the relationship between the rainstorm direction and catchment hydrologic response. Also, the maximum flood peak response is strongly dependent on the interaction of the peak of the rainstorm with the peak of the width function. The results indicate significant differences in runoff volume and runoff peak and their interannual variability due to rainstorm direction. The study has important implications for flood frequency predictions and developing flood resilience and watershed management strategies, particularly due to changing storm tracks under a warming climate.

42 ENGINEERING↗