SEARCH · Engineering Papers

Results for “GROUND WATER”

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 757 records · Page 42

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Snowfall and snow accumulation during the MOSAiC winter and spring seasons

Data from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition allowed us to investigate the temporal dynamics of snowfall, snow accumulation and erosion in great detail for almost the whole accumulation season (November 2019 to May 2020). We computed cumulative snow water equivalent (SWE) over the sea ice based on snow depth and density retrievals from a SnowMicroPen and approximately weekly measured snow depths along fixed transect paths. We used the derived SWE from the snow cover to compare with precipitation sensors installed during MOSAiC. The data were also compared with ERA5 reanalysis snowfall rates for the drift track. We found an accumulated snow mass of 38 mm SWE between the end of October 2019 and end of April 2020. The initial SWE over first-year ice relative to second-year ice increased from 50 % to 90 % by end of the investigation period. Further, we found that the Vaisala Present Weather Detector 22, an optical precipitation sensor, and installed on a railing on the top deck of research vessel Polarstern, was least affected by blowing snow and showed good agreements with SWE retrievals along the transect. On the contrary, the OTT Pluvio 2 pluviometer and the OTT Parsivel 2 laser disdrometer were largely affected by wind and blowing snow, leading to too high measured precipitation rates. These are largely reduced when eliminating drifting snow periods in the comparison. ERA5 reveals good timing of the snowfall events and good agreement with ground measurements with an overestimation tendency. Retrieved snowfall from the ship-based K a -band ARM zenith radar shows good agreements with SWE of the snow cover and differences comparable to those of ERA5. Based on the results, we suggest the K a -band radar-derived snowfall as an upper limit and the present weather detector on RV Polarstern as a lower limit of a cumulative snowfall range. Based on these findings, we suggest a cumulative snowfall of 72 to 107 mm and a precipitation mass loss of the snow cover due to erosion and sublimation as between 47 % and 68 %, for the time period between 31 October 2019 and 26 April 2020. Extending this period beyond available snow cover measurements, we suggest a cumulative snowfall of 98–114 mm.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Synthesis of the Candidate Topological Compound Ni 3 Pb 2

Spin-orbit coupling enables the realization of topologically nontrivial ground states. As spin-orbit coupling increases with increasing atomic number, compounds featuring heavy elements, such as lead, offer a pathway toward creating new topologically nontrivial materials. Here, by employing a high-pressure flux synthesis method, we synthesized single crystals of Ni 3 Pb 2 , the first structurally characterized bulk binary phase in the Ni-Pb system. Combining experimental and theoretical techniques, we examined structure and bonding in Ni 3 Pb 2 , revealing the impact of chemical substitutions on electronic structure features of importance for controlling topological behavior. From these results, we determined that Ni 3 Pb 2 completes a series of structurally related transition-metal-heavy main group intermetallic materials that exhibit diverse electronic structures, opening a platform for synthetically tunable topologically nontrivial materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Laser Measurements and Modeling of Shock Tunnel Freestream Velocity and Multispecies Thermal Nonequilibrium

Coherent anti-Stokes Raman scattering (CARS) and nitric oxide molecular tagging velocimetry (NO-MTV) were used to characterize the freestream in Sandia’s Hypersonic Shock Tunnel (HST) using a burst-mode laser operated at 100 kHz. Experiments were performed at nominal freestream velocities of 3 and 4 km/s using both air and N 2 test gas. The CARS diagnostic provides nonequilibrium characterization of the flow by measuring vibrational and rotational temperatures of N 2 and O 2 , which are compared to NO temperatures from separate laser absorption experiments. Simultaneous, collinear freestream velocities were measured using NO MTV along with pitot pressure measurements. Furthermore, this extensive freestream dataset is compared to nonequilibrium CFD that is capable of modeling species-specific, vibrational temperatures throughout the nozzle expansion. Significant nonequilibrium between vibrational and rotational temperatures was measured at each flow condition. N 2 exhibits the most nonequilibrium, followed by O 2 and NO. The CFD model captures this trend, although it consistently overpredicts N 2 and O 2 vibrational temperatures. At 3 km/s, the modeled NO nonequilibrium is underpredicted, whereas it is overpredicted at 4 km/s. Good agreement is seen between CFD and the velocity and rotational temperature measurements. Experiments with water added to the test gas yielded no discernable difference in vibrational relaxation.

Aerodynamics↗

DOE OSTI2020

Collaborative Research: Natural Organic Matter and Microbial Controls on Mobilization/Immobilization of I and Pu in Soils and Waters Affected by Radionuclide Releases in USA and Japan

In this project, the relationship between natural organic matter (NOM) and two radioactive elements that are relevant to nuclear waste disposal were studied: Plutonium (Pu) and Iodine. The human and environmental risks associated with Pu stem mainly from the very long half-lives of several of its isotopes ( 238 Pu, 88 yr; 239 Pu, 24,100 yr; 240 Pu, 6560 yr) and its radiotoxicity. Understanding Pu biogeochemical behavior in both near-field (>10-11M) and far-field scales (<10-11M) is imperative to the development of approaches for reprocessing Pu, remediation of Pu contamination and accurate assessment of risks posed by disposal practices for Pu-bearing wastes. The environmental mobility of Pu can be affected by redox potential, pH, adsorption, precipitation, complexation, colloid formation, and microbial activity, of which the first characteristic has the most profound influence. Numerous studies have shown high affinity of Pu towards NOM, as well as to mineral phases. NOM is ubiquitous in the environment, e.g., both fulvic and humic acids are able to reduce Pu(V,VI) to Pu(IV) and the redox potential of NOM is positively related to the abundance of phenolic/acidic OH groups. NOM can either facilitate or limit actinide migration, depending on specific biogeochemical conditions including pH, mineral and organic matter characteristics, etc. The other radionuclide of interest is radioiodine ( 129 I). 129I is a major by-product of nuclear fission and of serious concern to the Department of Energy (DOE) as it is among the top risk drivers at existing and potential radiowaste-contaminated sites. The risk of 129 I stems largely from its high bioconcentration factor (90% of the body’s iodine is accumulated in the thyroid), a high inventory at source terms, a very long-half life (16M years), and rapid mobility in the subsurface environment. As a consequence, 129 I has the lowest drinking water standard (1 pCi/L) among all radionuclides in the Federal Register. With a novel and sensitive gas chromatography-mass spectrometry (GC-MS) method developed in our lab, it is possible to quickly and simultaneously determine the distribution of 129 I and stable 127 I forms in environments, as low as 2 pCi/L for 129 I. This method was subsequently validated using accelerator mass spectroscopy, AMS. IO 3 - and organo-I were determined as major species in the groundwater of SRS and the Hanford Site, contrary to thermodynamic predictions that I- should be the dominant species at these sites. Mobility of 129 I was also demonstrated to depend greatly on the I species and its concentration, sediment pH, and redox state, with times to achieve equilibrium taking up to 12 weeks. Along the groundwater pathway in the F-Area of SRS, 129 I- supplied from the seepage basins was transformed to 129 IO 3 - and organo- 129 I with increasing iodine sediment sorption, causing the lower total 127 I and 129 I concentrations along the gradient transect of the waste plume. By contrast, groundwater 129 I concentrations in the wetlands (as high as 1617.3 pCi/L) were greatly elevated with respect to the source term (159.3 pCi/L). While the NOM promoted the uptake of 129 I to the wetland sediment, it also promoted the formation of soluble organic fraction. A small fraction of NOM that is bound to iodine can behave as a mobile organo-I source. Iodide was enzymatically incorporated into NOM, whereas both iodide and iodate were abiotically bound to NOM, under certain conditions. Iodate removal from the mobile aqueous phase can also occur through incorporation into carbonate (e.g., at the Hanford Site, USA). Thus immobilization and re-mobilization of iodine species were influenced by pH, Eh and the presence of NOM and metal oxides, which adds to the complexity of site remediation action. A ground-breaking result was to elucidate the products (i.e. organo-iodine moieties formed via enzymatic and non-enzymatic processes) at the molecular level by nuclear magnetic resonance (NMR) and electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICRMS). We found that iodine-NOM interactions may be influenced by NOM hydrophobic aliphatic moieties. From the perspective of ESI-FTICRMS, organo-iodine formulas were ascribed to the groups of unsaturated hydrocarbons, lignins and proteins. Iodate is likely abiotically reduced to reactive iodine species by lignin- and tannin-like compounds or carboxylic-rich alicyclic molecules (CRAM). We also investigated microbial mechanisms in iodine incorporation into NOM. We established that soil bacteria isolated from F-Area of SRS did not accumulate significant amounts of I- (0.2-2%). Intracellular uptake of I- decreases with increasing pH when pH ranged from 4 to 6. In contrast, 44 out of 84 strains isolated from the F-Area of SRS can transform I- to IO 3 - and organo-iodine. In some cases, oxidation was facilitated in the presence of H 2 O 2 . Microbes can also excrete organic acids that enhance I- oxidation by lowering the ambient pH and reacting with H 2 O 2 to form peroxy carboxylic acids. At lower pH values (≤5), H 2 O 2 hydrolysis was the driving force for iodide-oxidation; whereas, at pH ≥ 6, spontaneous decomposition of peroxy carboxylic acids, originating from H 2 O 2 and organic acids were the primary cause of iodide oxidation. Lastly, it was determined that microbial processes involved in Mn (II) are capable of directly oxidizing I- via enzymatic catalysis (i.e., multicopper oxidases), or indirectly through the formation of reactive oxygen species (ROS) and/or biogenic manganese oxides. ROS-mediated oxidation of I- was found to predominate at pH >5, whereas the enzymatic and Mn oxide pathways were more active at pH < 5. Together, this project has resulted in 9 publications in high-impact journals, and the training of 1 Ph.D and 4 undergraduate students.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

CACTI ARM Aerial Facility Measurements of Ice Nucleating Particles

The dataset comprises measures, using Colorado State University's Ice Spectrometer (IS, an immersion freezing device with a range from 0°C down to -26 to -29°C) of atmospheric ice nucleating particle (INP) concentrations taken on the Atmospheric Radiation Measurement (ARM) program Aerial Facilty (AAF) G-1 aircraft. INP measurements on the G-1 were collected from varied altitudes on different flights over the region of the Sierras de Córdoba mountain range of north-central Argentina, centred over ARM's Mobile Facility (AMF-1) near Villa Yacanto, where ground-based INP measures were being taken concurrently. Both studies took place as part of the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) Experiment. A filter sampling system was deployed on the G-1 to collect aerosol particles for post-processing of INPs collected on filters (measuring their immersion freezing ability) once returned to Colorado State University (CSU). Filter holders used were pre-cleaned aluminum in-line units loaded with pre-cleaned and pre-sterilized 47 mm dia. Nuclepore polycarbonate filters (0.2 mm pore size). Filters were drawn for varied times, resulting in varied volumes collected (51 to 1667 SL). Mass flow rate was recorded in real-time so that total sampled volume (at standard temperature and pressure) could be determined. A total of 34 sample filters were collected over the IOP, including 5 blanks. Filters were stored at -20°C freezer prior to frozen return to Colorado State University (CSU). . Processing to obtain spectra of INP number concentration active via the immersion freezing mechanism versus temperature was conducted using CSU's IS instrument (McCluskey et al., 2018). For measurment of INPs, collected aerosol particles were re-suspended in 7 mL of 0.02 µm-filtered deionized water. Aliquots of each suspension, and serial dilutions, were dispensed into trays which were fit into aluminum blocks in the IS. Samples are cooled at 0.33°C min-1 and the freezing temperatures of wells recorded automatically. Cumulative INP concentrations were determined by first calculating the INPs per mL of suspension based on Vali (1971) and then converting to concentration per standard liter of air using the proportion of the total liquid sample dispensed and the air sample volumes. Aliquots of suspensions from selected samples were also heat treated (95°C for 20 min) to denature and deactivate biological INPs, and digested in 10% H2O2 at 95°C under UV-B to remove all organic carbon INPs. McCluskey, C. S., J. Ovadnevaite, M. Rinaldi, J. Atkinson, F. Belosi, D. Ceburnis, S. Marullo, T. C. J. Hill, U. Lohmann, Z. A. Kanji, C. O’Dowd, S. M. Kreidenweis, P. J. DeMott, 2018: Marine and Terrestrial Organic Ice Nucleating Particles in Pristine Marine to Continentally-Influenced Northeast Atlantic Air Masses, Journal of Geophysical Research: Atmospheres, 123, 6196–6212, https://doi.org/10.1029/2017JD028033. Vali, G., 1971: Quantitative evaluation of experimental results on the heterogeneous freezing nucleation of supercooled liquids. J. Atmos. Sci., 28, 402–409.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

International Applications for Floating Solar PV (FPV) [Slides]

Floating solar is a comparatively new concept. Floating solar is about to form a third pillar in the market for solar energy next to ground mounted and rooftop solar. The Asia-Pacific region already dominates the fast-emerging market in floating solar PV (FPV), and it looks set to build on that position as more nations get on board and costs fall. Asia-Pacific is set to increase its market share from 74% in 2020 to 87% in 2026 in terms of global installed FPV capacity. However, since floating solar systems are still a new application, there is a lack of information on its performance, cost, technical complexity, and operations and maintenance (O&M). This seminar series seeks to partially remedy that through presenting case studies that highlight the benefits, challenges, and applications of FPV in South and Southeast Asia.

14 SOLAR ENERGY↗

DOE OSTI2023

LA100 Equity Strategies. Chapter 2: Procedural Justice

The LA100 Equity Strategies project synthesizes community guidance with robust research, modeling, and analysis to identify strategy options that can increase equitable outcomes in Los Angeles' clean energy transition. Grounded in the analysis of past and ongoing energy inequities and engagement with underserved communities, the project presents community-guided and community-tailored strategies that aim to operationalize recognition and procedural justice. This chapter focuses on procedural justice, examining priorities identified during the community engagement portion of the LA100 Equity Strategies project. This process, and our approach to partnering with community-based organizations (CBOs) and the communities they serve, was developed from the baseline analysis in Chapter 1, which centers on recognition justice, examining past and current inequities in LA. Recognition, procedural, and distributional justice are the three tenets of energy justice around which the LA100 Equity Strategies project is organized. Procedural justice prioritizes fair, equitable, and inclusive participation in the decision-making process. This tenet's practical application entails who is invited and able to participate, whose voices are considered as decisions are made, the co-development of procedures to inform this deliberative process, and who has access to formal measures of regulation and accountability. Engaging with Angelenos to examine the causes of inequities and identify impact areas and priorities that center community experiences, values, and goals represents a first step in developing energy equity strategies for the Los Angeles Department of Water and Power (LADWP) to achieve distributional justice - the just and equitable distribution of energy benefits and burdens in LA's energy transition. Equity, as community members insisted, is about making - and following through with - a commitment to prioritize historically underserved and overburdened communities. Community engagement is a principal method for operationalizing this commitment, guiding our analytic approach and potential equity strategies. Using procedural justice as an analytical tool, this chapter presents the LA100 Equity Strategies approach to community engagement from 2021 to 2023 and the results of this process in relation to community-identified barriers and burdens impacting procedural justice outcomes. We analyze the procedural elements of reaching the equity goals Angelenos set to inform future LADWP decision-making and program development.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2026

A compact, low-power epithermal neutron counter for lunar water detection

The detection and characterization of lunar water are critical for enabling sustainable human and robotic exploration of the Moon. Orbital neutron spectrometers, such as instruments on Lunar Prospector and the Lunar Reconnaissance Orbiter, have revealed hydrogen-rich regions near the poles but are limited by coarse spatial resolution and low counting efficiency. We present a compact, lightweight, and low-power epithermal neutron detector based on boron-coated silicon imagers, designed to probe subsurface hydrogen at decimeter scales from mobile platforms such as lunar rovers. This instrument leverages the high neutron capture cross-section of 10 B to convert epithermal neutrons into detectable α and 7 Li ions in a fully-depleted silicon imager, providing a unique event topology to identify neutrons while suppressing backgrounds. Monte Carlo simulations demonstrate that a 3 μm boron layer achieves optimal neutron detection efficiency, further enhanced with polyethylene moderation to improve sensitivity to the 0.4 eV–500 keV epithermal energy range. For a 10 cm 2 active area, the detector achieves sensitivity to H 2 O weight fractions as low as 0.01 wt% in a 15 minute measurement. This scalable, portable, low-mass design is well-suited for integration into upcoming Artemis and commercial lunar rovers, providing a transformative capability for in-situ resource prospecting and ground-truth validation of orbital measurements.

Detector modelling and simulations I (interaction ↗

DOE OSTI2020

ACCESS: Confirmation of No Potassium in the Atmosphere of WASP-31b

We present a new optical (400–950 nm) transmission spectrum of the hot Jupiter WASP-31b (M = 0.48 M {sub J}; R = 1.54 R {sub J}; P = 3.41 days), obtained by combining four transit observations. These transits were observed with IMACS on the Magellan Baade Telescope at Las Campanas Observatory as part of the ACCESS project. We investigate the presence of clouds/hazes in the upper atmosphere of this planet, as well as the contribution of stellar activity on the observed features. In addition, we search for absorption features of the alkali elements Na i and K i, with particular focus on K i, for which there have been two previously published disagreeing results. Observations with Hubble Space Telescope (HST)/STIS detected K i, whereas ground-based low- and high-resolution observations did not. We use equilibrium and nonequilibrium chemistry retrievals to explore the planetary and stellar parameter space of the system with our optical data combined with existing near-IR observations. Our best-fit model is that with a scattering slope consistent with a Rayleigh slope (α=5.3{sub −3.1}{sup +2.9}), high-altitude clouds at a log cloud top pressure of −3.6{sub −2.1}{sup +2.7} bars, and possible muted H{sub 2}O features. We find that our observations support other ground-based claims of no K i. Clouds are likely why signals like H{sub 2}O are extremely muted and Na or K cannot be detected. We then juxtapose our Magellan/IMACS transmission spectrum with existing VLT/FORS2, HST/WFC3, HST/STIS, and Spitzer observations to further constrain the optical-to-infrared atmospheric features of the planet. We find that a steeper scattering slope (α = 8.3 ± 1.5) is anchored by STIS wavelengths blueward of 400 nm and only the original STIS observations show significant potassium signal.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Probabilistic Multi-Hazard Performance Assessment of Concrete Structures in Nuclear Installations

Concrete structures in nuclear installations are subject to time-dependent degradation mechanisms that can deteriorate their physical and mechanical properties, potentially exacerbating the risk of structural failure under external forces such as a seismic event. Previous research has extensively investigated the seismic response of nuclear concrete structures and the associated risk, as well as their effect on structural components safety margins. However, substantial work is still necessary to incorporate concrete aging effects into such evaluations. In fact, most models in the literature assume pristine concrete conditions and do not account for the impact of aging on the structural components’ fragility curves. This work identifies relevant time-dependent degradation mechanisms and provides simplified models to predict the the evolution of key material properties based on data from the literature. Namely, this work focuses on the aging effects of corrosion, alkali–silica reaction (ASR), and irradiation on reinforced concrete within US Department of Energy (DOE) nuclear facilities and nuclear power plants (NPP) structures. Furthermore, degradation models based on literature data are presented that define the relationship between probabilistic material properties and the concrete’s age. In this work, sampled material properties served as input for a simplified finite element model (FEM) of a critical nuclear structural system, with the output of the FEM being the seismic response for a given ground motion. The results of the FEM were then used within a probabilistic performance assessment with a statistically significant number of samples. The research presented herein addresses the detrimental effects of hazards caused by natural phenomena on deteriorated concrete elements of nuclear installations. This work directly benefits the safety analysis performed on US DOE/ National Nuclear Security Administration (NNSA) nuclear facilities located in areas prone to seismic activity. The results presented herein could aid in the improvement of DOE-STD-1020, the DOE Standard that addresses seismic risk analysis and capacity evaluation in DOE facilities. DOE-STD-1020 refers to the requirements in American Society of Civil Engineers (ASCE) 4-98, now superseded by ASCE 4-16, that shall be met in performing dynamic response analyses and generating in-structure response spectra, provided that such requirements are consistent with the requirements of ASCE/Structural Engineering Institute (SEI) 43-05. Moreover, the results presented herein could also aid in the updating of section C3.1.1. of ASCE 4-16 to account for the effects of aging on the stiffness of reinforced elements and American Concrete Institute (ACI) 349.3R-18, “Report on Evaluation and Repair of Existing Nuclear Safety-Related Concrete Structures.” Ultimately, this work can assist the risk assessment of potential lifetime extension of the existing US commercial nuclear fleet (light water reactors) and the safety analysis of the emerging advanced nuclear reactors. The proposed proof-of-concept methodology employs open-source DOE computational tools and is transferable to commercial software commonly used by engineering firms.

42 ENGINEERING↗

DOE OSTI2023

Intact mass spectrometry screening to optimize hydroxyl radical dose for protein footprinting

Hydroxyl radical protein footprinting (HRPF) using synchrotron radiation is a well-validated method to assess protein structure in the native solution state. In this method, X-ray radiolysis of water generates hydroxyl radicals that can react with solvent accessible side chains of proteins, with mass spectrometry used to detect the resulting labeled products. An ideal footprinting dose provides sufficient labeling to measure the structure but not so much as to influence the results. The optimization of hydroxyl radical dose is typically performed using an indirect Alexa488 fluorescence assay sensitive to hydroxyl radical concentration, but full evaluation of the experiment's outcome relies upon bottom-up liquid chromatography mass spectrometry (LC-MS) measurements to directly determine sites and extent of oxidative labeling at the peptide and protein level. A direct evaluation of the extent of labeling to provide direct and absolute measurements of dose and “safe” dose ranges in terms of, for example, average numbers of labels per protein, would provide immediate feedback on experimental outcomes prior to embarking on detailed LC-MS analyses. To this end, we describe an approach to integrate intact MS screening of labeled samples immediately following exposure, along with metrics to quantify the extent of observed labeling from the intact mass spectra. Intact MS results on the model protein lysozyme were evaluated in the context of Alexa488 assay results and a bottom-up LC-MS analysis of the same samples. In conclusion, this approach provides a basis for placing delivered hydroxyl radical dose metrics on firmer technical grounds for synchrotron X-ray footprinting of proteins, with explicit parameters to increase the likelihood of a productive experimental outcome. Further, the method directs approaches to provide absolute and direct dosimetry for all types of labeling for protein footprinting.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Quantitative imaging and automated fuel pin identification for passive gamma emission tomography

Compliance of member States to the Treaty on the Non-Proliferation of Nuclear Weapons is monitored through nuclear safeguards. The Passive Gamma Emission Tomography (PGET) system is a novel instrument developed within the framework of the International Atomic Energy Agency (IAEA) project JNT 1510, which included the European Commission, Finland, Hungary and Sweden. The PGET is used for the verification of spent nuclear fuel stored in water pools. Advanced image reconstruction techniques are crucial for obtaining high-quality cross-sectional images of the spent-fuel bundle to allow inspectors of the IAEA to monitor nuclear material and promptly identify its diversion. In this work, we have developed a software suite to accurately reconstruct the spent-fuel cross sectional image, automatically identify present fuel rods, and estimate their activity. Unique image reconstruction challenges are posed by the measurement of spent fuel, due to its high activity and the self-attenuation. While the former is mitigated by detector physical collimation, we implemented a linear forward model to model the detector responses to the fuel rods inside the PGET, to account for the latter. The image reconstruction is performed by solving a regularized linear inverse problem using the fast-iterative shrinkage-thresholding algorithm. We have also implemented the traditional filtered back projection (FBP) method based on the inverse Radon transform for comparison and applied both methods to reconstruct images of simulated mockup fuel assemblies. Higher image resolution and fewer reconstruction artifacts were obtained with the inverse-problem approach, with the mean-square-error reduced by 50%, and the structural-similarity improved by 200%. We then used a convolutional neural network (CNN) to automatically identify the bundle type and extract the pin locations from the images; the estimated activity levels finally being compared with the ground truth. The proposed computational methods accurately estimated the activity levels of the present pins, with an associated uncertainty of approximately 5%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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 OSTI2022

What explains the year-to-year variation in growing season timing of boreal black spruce forests?

Amplified climate warming in high latitudes is expected to affect growing season timing of the vast boreal biome. It is unclear whether the presence of permafrost (perennially frozen ground) might have an influence on changes in growing season timing. Here, this study examined how different environmental variables explained, either directly or indirectly, the variation in growing season timing of boreal forest stands with and without permafrost. We expected that environmental variables explaining the variation in growing season timing differed or had different explanatory power depending on permafrost presence or absence. The growing season was delineated from daily gross primary productivity (GPP) time series derived from 40 site-year data of net ecosystem carbon dioxide exchange measured with eddy covariance techniques over five black spruce (Picea mariana [Mill.])-dominated boreal forest stands in North America. In permafrost-free forest stands, a combination of start in canopy ‘green-up’ in spring and the timing of air and soil temperature increasing above freezing explained the start-of-season (SOS GPP ). Results from commonality analysis and structural equation modeling suggest that canopy ‘green-up’ and air temperature directly affected SOS GPP in permafrost-free forest stands. In addition, soil temperature acted as mediator for an indirect effect of air temperature on SOS GPP . In contrast, none of the environmental variables, or their combination, explained the variation in SOS GPP in forest stands with permafrost. The explanatory power of environmental variables was more consistent regarding the end-of-season (EOS GPP ). In both, forest stands with and without permafrost, EOS GPP was directly explained by mean soil water content in the fall and the first day of continuous snowpack formation. A better understanding how environmental variables control SOS GPP and EOS GPP in forest stands with and without permafrost will help to refine parameterizations of the boreal biome in Earth system models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Evaluating cloud liquid detection against Cloudnet using cloud radar Doppler spectra in a pre-trained artificial neural network

Detection of liquid-containing cloud layers in thick mixed-phase clouds or multi-layer cloud situations from ground-based remote-sensing instruments still poses observational challenges, yet improvements are crucial since the existence of multi-layer liquid layers in mixed-phase cloud situations influences cloud radiative effects, cloud lifetime, and precipitation formation processes. Hydrometeor target classifications such as from Cloudnet that require a lidar signal for the classification of liquid are limited to the maximum height of lidar signal penetration and thus often lead to underestimations of liquid-containing cloud layers. Here we evaluate the Cloudnet liquid detection against the approach of Luke et al. (2010) which extracts morphological features in cloud-penetrating cloud radar Doppler spectra measurements in an artificial neural network (ANN) approach to classify liquid beyond full lidar signal attenuation based on the simulation of the two lidar parameters particle backscatter coefficient and particle depolarization ratio. We show that the ANN of Luke et al. (2010) which was trained under Arctic conditions can successfully be applied to observations at the mid-latitudes obtained during the 7-week-long ACCEPT field experiment in Cabauw, the Netherlands, in 2014. In a sensitivity study covering the whole duration of the ACCEPT campaign, different liquid-detection thresholds for ANN-predicted lidar variables are applied and evaluated against the Cloudnet target classification. Independent validation of the liquid mask from the standard Cloudnet target classification against the ANN-based technique is realized by comparisons to observations of microwave radiometer liquid-water path, ceilometer liquid-layer base altitude, and radiosonde relative humidity. In addition, a case-study comparison against the cloud feature mask detected by the space-borne lidar aboard the CALIPSO satellite is presented. Three conclusions were drawn from the investigation. First, it was found that the threshold selection criteria of liquid-related lidar backscatter and depolarization alone control the liquid detection considerably. Second, all threshold values used in the ANN framework were found to outperform the Cloudnet target classification for deep or multi-layer cloud situations where the lidar signal is fully attenuated within low liquid layers and the cloud radar is able to detect the microphysical fingerprint of liquid in higher cloud layers. Third, if lidar data are available, Cloudnet is at least as good as the ANN. The times when Cloudnet outperforms the ANN in liquid detections are often associated with situations where cloud dynamics smear the imprint of cloud microphysics on the radar Doppler spectra.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Computationally Accelerated Discovery and Experimental Demonstration of High-Performance Materials for Advanced Solar Thermochemical Hydrogen Production

This project achieved its overarching goal of accelerating the discovery and validation of solar thermochemical hydrogen (STCH) materials through a tightly integrated approach that combined high-throughput computational screening, advanced machine learning (ML), and experimental testing. Guided by the objectives outlined in the Statement of Project Objectives (SOPO), our work fulfilled all major milestones across four technical tasks and delivered scientific breakthroughs and practical tools that significantly exceeded the original scope of the project. We began by addressing the challenge of predicting material phase stability through machine learning. A novel Python module was developed to generate thousands of meaningful features from composition, structure, and electronic properties, enabling rapid and reproducible ML model development. Using these tools, we trained a model to predict temperature-dependent Gibbs energies (G(T)) for inorganic crystalline materials with near-chemical accuracy—roughly 40 meV/atom—marking the first such descriptor of its kind. We also introduced a new machine-learned tolerance factor, τ, that accurately predicted perovskite formability with over 90% success, outperforming traditional heuristic models, such as the Goldschmidt tolerance factor. These capabilities allowed for rapid and accurate predictions of phase stability across a vast oxide composition space, setting the stage for high-throughput thermodynamic screening. Building on this foundation, we conducted an extensive computational screening of candidate STCH oxide materials. Over 1.1 million perovskite compositions were evaluated using the τ descriptor, leading to the identification of more than 27,000 predicted stable structures. Using density functional theory (DFT), we refined over 68,000 multinary perovskite structures and computed oxygen vacancy formation energies for over 1,300 ternary and double perovskites. These calculations enabled us to isolate compounds with redox behavior consistent with STCH requirements and resulted in a public dataset now hosted on the Materials Project. Recognizing that thermodynamic screening alone is insufficient, we addressed kinetic limitations by developing a suite of tools to estimate transition state (TS) energies for key redox reactions. We implemented a novel bounding approach that provides lower and upper estimates of TS energies with dramatically reduced computational cost, requiring less than 10% of the CPU time of a full nudged elastic band (NEB) calculation while maintaining high accuracy. This enabled rapid evaluation of over 200 reaction pathways across 90 materials. To further accelerate screening, we developed a SISSO-based ML model to predict diffusion barriers with a 96.7% success rate in classifying fast vs. slow materials, supporting a robust, data-driven framework for assessing redox kinetics. Experimental validation was critical to confirming the predictive power of our models. We synthesized and tested a wide array of candidate materials, including Mn-doped hercynite and several Gd- and La-based perovskites. Notably, Sr 0.4 Gd 0.6 Mn 0.6 Al 0.4 O 3 (SGMA) and Gd 0.5 La 0.5 Co 0.5 Fe 0.5 O 3 (GLCF) emerged as leading STCH materials, exhibiting robust redox cycling and high hydrogen yields exceeding 150 µmol H 2 /g per cycle. These materials also retained over 50% of their hydrogen productivity under high-conversion conditions (H 2 O:H 2 = 1333:1), demonstrating strong thermodynamic favorability and promising performance under industrially relevant scenarios. Additional candidates, such as La 2 MnNiO 6 (L2MN), were found to produce even higher yields than ceria under standard STCH conditions. Our collaborators at Sandia National Laboratories confirmed these findings using high-temperature X-ray diffraction and thermogravimetric analysis, observing stable phase evolution and reversible redox activity. In several respects, the project went beyond the goals initially outlined in the SOPO. We published 17 peer-reviewed articles, including a large dataset of over 66,000 theoretical perovskites and a new structure prediction method (SPuDS-DFT) that accurately identifies ground-state structures at a fraction of the cost of traditional DFT. We demonstrated that our machine-learned G(T) model offers accuracy rivaling quasiharmonic calculations while being orders of magnitude faster. In partnership with the Materials Project, we made our datasets openly available, providing a powerful new resource for the broader materials science community. The combined computational and experimental advances of this project represent a significant advance in STCH materials discovery. By creating a robust, generalizable, and open workflow for thermodynamic and kinetic screening, and validating key findings through synthesis and reactor testing, we have provided a practical and scalable pathway for the rapid identification of new redox-active materials. The tools, data, and materials developed under this project are already supporting ongoing research and have laid the groundwork for the next generation of solar fuel technologies.

08 HYDROGEN↗

DOE OSTI2023

Evaluation of aerosol–cloud interactions in E3SM using a Lagrangian framework

A Lagrangian framework is used to evaluate aerosol–cloud interactions in the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) version 1 (E3SMv1) for measurements taken at Graciosa Island in the Azores where a U.S. Department of Energy Atmosphere Radiation Measurement (ARM) site is located. This framework uses direct measurements of cloud condensation nuclei (CCN) concentration (instead of relying on satellite retrievals of aerosol optical depth) and incorporates a suite of ground-based ARM measurements, satellite retrievals, and meteorological reanalysis products that when applied to over a 1500 trajectories provides key insights into the evolution of low-level clouds and aerosol radiative forcing that is not feasible from a traditional Eulerian analysis framework. Significantly lower concentrations (40 %) of surface CCN concentration are measured when precipitation rates in 48 h back trajectories average above 1.2 mm d -1 in the Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (IMERG) product. The depletion of CCN concentration when precipitation rates are elevated is nearly twice as large in the ARM observations compared to E3SMv1 simulations. The model CCN concentration bias remains significant despite modifying the autoconversion and accretion rates in warm clouds. As the clouds in trajectories associated with larger surface-based CCN concentration advect away from Graciosa Island, they maintain higher values of droplet number concentrations (N d ) over multiple days in observations and E3SM simulations compared to trajectories that start with lower CCN concentrations. The response remains robust even after controlling for meteorological factors such as lower troposphere stability, the degree of cloud coupling with the surface, and island wake effects. E3SMv1 simulates a multi-day aerosol effect on clouds and a Twomey radiative effect that is within 30 % of the ARM and satellite observations. However, the mean cloud droplet concentration is more than 2–3 times larger than in the observations. While Twomey radiative effects are similar amongst autoconversion and accretion sensitivity experiments, the liquid water path and cloud fraction adjustments are positive when using a regression model as opposed to negative when using the present-day minus pre-industrial aerosol emissions approach. This result suggests that tuning the autoconversion and accretion alone is unlikely to produce the desired aerosol susceptibilities in E3SMv1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Community Geothermal: Soil Conductivity, Borehole Design, Energy Models, and Load Data for a Residential System Development - Hinesburg, VT

This dataset contains materials from the Coalition for Community-Supported Affordable Geothermal Energy Systems (C2SAGES) project, which evaluated the techno-economic feasibility of a community geothermal system for a residential development in Hinesburg, VT. The dataset includes detailed soil conductivity test reports, energy models, borehole design reports, hourly energy loads for heating, cooling, and hot water, and design layouts. EnergyPlus was used to model building energy loads, and Modelica software was applied for geothermal loop sizing based on these loads and soil conductivity results. Python scripts for network design further refined the models. Key files include PDF reports on borehole design (with projections for 1-year, 15-year, and 30-year systems), soil conductivity test results, EnergyPlus modeling outputs, and 2D/3D design drawings in PDF, DWG, and DXF formats. Python notebooks for network design and OnePipe model files are also provided, with Modelica required for viewing certain files. Outputs and modeling data are in various formats including CSV, JPG, HTML, and IDF, with units and data clearly labeled to support understanding of system design and performance for the proposed geothermal solution.

15 GEOTHERMAL ENERGY↗