Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “lack”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Neutron irradiation & thermomechanical experiment (NITE) - design

For the reliable long-term operation of fusion power plants, it is crucial to understand and predict the lifetime of materials in use. These materials include all structural and functional materials utilized at the first wall, blanket, magnets, and shielding. The key challenge is, that the harsh environment including high heat fluxes, high thermal stress and stress cycling, neutron irradiation, and sputtering on such materials should not be viewed separately. Currently, the synergistic loads cannot be evaluated experimentally because of the lack of adequate facilities. The purpose of that work is to design a synergetic Neutron Irradiation and Thermomechanical Experiment (NITE) for fusion materials. This design will leverage the existing Advanced-Test-Reactor (ATR), a fission reactor at the Idaho National Laboratory. We also acknowledge that with existing fission reactors the exact fusion condition cannot be created, and the limitations are critically discussed. The combination of neutron irradiation with a high heat flux is the focus. This is realized with an irradiation capsule design that includes a TRISO fueled region inside the capsule to enable a steady-state heat flux on one side of the specimen. In conclusion, the experimental design modeling showed that steady-state heat fluxes of 2.4 MW/m 2 with a thermal gradient of above 250°C can be achieved in a 5 mm thick specimen.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Tikiri—Towards a lightweight blockchain for IoT

Internet of Things (IoT) platforms have been deployed in several domains to enhance efficiency of business process and improve productivity. Most IoT platforms comprise of heterogeneous software and hardware components which can potentially introduce security and privacy challenges. Blockchain technology has been proposed as one of the solutions to realize IoT security by leveraging the (a) Immutable ledger, (b) Decentralized architecture and (c) Strong cryptography primitives. However, integrating blockchain platforms with IoT based applications presents several challenges due to lack of (a) acceptable performance on resource-constrained devices, (b) high transaction throughput, (c) keyword-based search and retrieve, (d) transaction back pressure operations, and (e) real-time response. In this paper, we propose a lightweight blockchain platform, “Tikiri”, for resource-constrained IoT devices. Tikiri uses Apache Kafka for the consensus and proposes new blockchain architecture to handle real-time transaction execution on the blockchain. Tikiri is characterized by functional programming and actor-based smart contract platform that realizes concurrent execution of transactions in the blockchain. Tikiri realizes a lightweight and scalable blockchain that can provides performance on the resource-constrained IoT devices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

CG-Kit: Code Generation Toolkit for performant and maintainable variants of source code applied to Flash-X hydrodynamics simulations

CG-Kit is a new Code Generation tool-Kit that we have developed as a part of the solution for portability and maintainability for multiphysics computing applications. The development of CG-Kit is rooted in the urgent need created by the shifting landscape of high-performance computing platforms and the algorithmic complexities of a particular large-scale multiphysics application: Flash-X. To efficiently use computing resources on a heterogeneous node, an application must have a map of computation to resources and a mechanism to move the data and computation to the resources according to the map. Most existing performance portability solutions are focussed on abstracting the expression of computations so that a unified source code can be specialized to run on different resources. However, such an approach is insufficient for a code like Flash-X, which has a multitude of code components that can be assembled in various permutations and combinations to form different instances of applications. Similar challenges apply to any code that has composability, where a single specified way of apportioning work among devices may not be optimal. Additionally, use cases arise where the optimal control flow of computation may differ for different devices while the underlying numerics remain identical. This combination leads to unique challenges including handling an existing large code base in Fortran and/or C/C++, subdivision of code into a great variety of units supporting a wide range of physics and numerical methods, different parallelization techniques for distributed and shared memory systems and accelerator devices, and heterogeneity of computing platforms requiring coexisting variants of parallel algorithms. All of these challenges demand that scientific software developers apply existing knowledge about domain applications, algorithms, and computing platforms to determine custom abstractions and granularity for code generation. There is a critical lack of tools to tackle those problems. CG-Kit is designed to fill this gap by providing a user with the ability to express their desired control flow and computation-to-resource map in the form a pseudocode-like recipe. It consists of standalone tools that can be combined into highly specific and, we argue, highly effective portability and maintainability toolchains. Here we present the design of our new tools: parametrized source trees, control flow graphs, and recipes. The tools are implemented in Python. They are agnostic to the programming language of the source code targeted for code generation. In conclusion, we demonstrate the capabilities of the toolkit with two examples, first, multithreaded variants of the basic AXPY operation, and second, variants of parallel algorithms within a hydrodynamics solver, called Spark, from Flash-X that operates on block-structured adaptive meshes.

Algorithmic portability↗
DOE OSTI2026

The role of quantum computing in advancing scientific high-performance computing: A perspective from the ADAC institute

Quantum computing (QC) has gained significant attention over the past two decades due to its potential for speeding up classically demanding tasks. This transition from an academic focus to a thriving commercial sector is reflected in substantial global investments. While advancements in qubit counts and functionalities continue at a rapid pace, current quantum systems still lack the scalability for practical applications, facing challenges such as too high error rates and limited coherence times. Here, this perspective paper examines the relationship between QC and high-performance computing (HPC), highlighting their complementary roles in enhancing computational efficiency. It is widely acknowledged that even fully error-corrected QC will not be suited for all computational tasks. Rather, future compute infrastructures are anticipated to employ quantum acceleration within hybrid systems that integrate HPC and QC. While QC can enhance classical computing, traditional HPC remains essential for maximizing quantum acceleration. This integration is a priority for supercomputing centers and companies, sparking innovation to address the challenges of merging these technologies. The novelty of this work lies in its unique perspective, reflecting the collective insights of the Accelerated Data Analytics and Computing (ADAC) Institute, a global consortium of over 20 leading HPC centers. Recognizing the growing importance of QC, ADAC established a Quantum Computing Working Group in 2023 to foster collaboration and knowledge-sharing among its members. This paper synthesizes insights from the group’s collaborative efforts and incorporates findings from a member survey that captures shared experiences, ongoing projects, and strategic directions. By outlining the current landscape and challenges of QC integration into HPC ecosystems, this work offers HPC specialists practical and forward-looking guidance on the opportunities and implications of QC in computationally intensive endeavors.

Accelerated Data Analytics and↗
DOE OSTI2026

DTLMod: A simulation framework for in situ workflow optimization

In situ processing workflows have become essential for coping with the explosion in data volume and velocity in large-scale scientific computing, providing domain scientists with early insights at runtime. Multiple frameworks implement this paradigm through a data transport layer (DTL), offering different data access modes and deployment schemes, but researchers currently lack the appropriate tools to assess design and deployment options before committing to costly real experiments. We introduce DTLMod, an open-source simulated DTL that enables performance evaluation of in situ workflow configurations at scale. Built on SimGrid, it links into any SimGrid-based simulator and is available in C++ and Python. We evaluate DTLMod along four axes: scalability (tens of thousands of simulated processes across interconnected clusters in seconds, with linear memory scaling), versatility (three implementation variants trading fidelity for speed), accuracy (simulated times faithfully reflecting real behavior), and practical utility (two use cases demonstrating evidence-based workflow design decisions).

Suter, Fred [ORNL] (ORCID:0000000319021955)↗
DOE OSTI2021

The solubility of thorium in carbonate-bearing solutions at hydrothermal conditions

Thorium mineralization is frequently hosted in carbonate-bearing rocks, and thorium commonly substitutes into the structures of carbonate-bearing minerals that have precipitated from or been modified by hydrothermal fluids. Given this common association, it is reasonable to consider the hypothesis that the presence of carbonate ligands in hydrothermal solutions promotes the transport of Th through the formation of stable aqueous complexes. Our ability to evaluate this hypothesis, however, is hindered by the lack of experimental data for Th-carbonate species at conditions beyond ambient. The low-temperature data indicate that carbonate is a strong complexing agent for Th. In this contribution, we investigate the solubility of Th in carbonate-bearing fluids relevant to natural systems (0.05–0.5 m NaHCO 3 /Na 2 CO 3 ; pH T ~ 7.8–9.8) at elevated temperature (175–250 °C). We demonstrate that, in contrast to the behavior of Th at low temperature, the stability of Th-carbonate complexes is not sufficient for them to predominate at these conditions. Instead, the solubility of Th is governed by hydrolysis reactions. Under the experimental conditions investigated, the predominant hydroxyl complexes are Th(OH) 4 0 and Th(OH) 5 -. Furthermore, thermodynamic formation constants were derived for these species at the temperatures considered in our experiments (log β 4 = 43.34 and 44.31 at 175 and 200 °C, respectively, and log β 5 = 46.15 and 47.9 at 225 and 250 °C, respectively) to permit forward modeling of Th mobility in natural systems. Our study indicates that carbonate ions are unlikely to play a role in transporting Th in hydrothermal fluids. Summarizing the results of this study and our previous studies of the solubility of Th in hydrothermal fluids, we conclude that SO 4 2- is the primary ligand responsible for the hydrothermal transport of Th.

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

Consistent Controls on Trace Metal Micronutrient Speciation in Wetland Soils and Stream Sediments

Trace metal are essential for microbially-mediated biogeochemical processes occurring in anoxic wetland soils and stream bed sediments, but low availability of these elements may inhibit anaerobic element cycling and transformations. Solid-phase speciation is likely a critical control on trace metal availability but has seen limited study in anoxic systems having concentrations similar to geological background levels, where metal limitations may be most prevalent. We have investigated trace metal concentrations and solid-phase speciation in three freshwater subsurface aquatic systems: marsh wetland soils, riparian wetland soils, and the sediments of a streambed. These systems displayed low solid-phase trace metal concentrations, generally at or below geological background levels, which generally followed the trend Zn > Cu ≈ Ni > Co and showed no correlation with major element compositions. All soils and sediments were dominated by quartz but varied in clay mineralogy as well as the organic matter, total sulfur, and total iron contents. X-ray absorption near-edge structure (XANES) spectroscopy shows that sulfur speciation in both wetlands is dominated by organic sulfur. Elemental sulfur and iron sulfides together made up <25% of the sulfur in the wetland soils, but the distribution between inorganic and organic forms was reversed in the stream sediments. Ferrous and ferric iron in clay minerals were common species identified by both XANES and extended X-ray absorption fine structure (EXAFS) spectroscopies at all sites. Iron(III) oxides were substantial components in all but the marsh wetland soils. Quantitative analysis of copper, nickel, and zinc XANES spectra revealed similar metal speciation across all sites. Copper speciation was dominated by sulfides, adsorbed species, and minor amounts of copper bound to organic matter; no metallic copper was detected. Nickel speciation also varied little and was dominated by nickel in clay mineral octahedral sheets and nickel sulfide, with adsorbed species also present. Zinc speciation was slightly more varied, with the marsh wetland soils and stream bed sediments containing adsorbed species, zinc associated with clay mineral structures, and zinc bound to reduced sulfur groups on organic matter, whereas the riparian wetland soils lacked clay-associated zinc but contained zinc sulfide. Trace metals bound to reduced sulfur occurred at every site, with a greater sulfur-bound fraction for copper. The fractional abundance of sulfur-bound species showed no relationship with soil or sediment total sulfur content, which varied by two orders of magnitude. More broadly, the observations in this study suggest that trace metal speciation in freshwater wetland soils and stream sediments is consistently dominated by a small set of recurring components which are distinct for each metal. Furthermore, this may represent a general geochemical phenomenon in anoxic soils and sediments containing trace metals at background concentrations (as low as 3 µg g -1 ) that was not predicted from systems that are contaminated with or naturally-enriched in copper, nickel, or zinc.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Titanium isotope systematics of refractory inclusions: Echoes of molecular cloud heterogeneity

Calcium-aluminum-rich inclusions (CAIs) are highly refractory objects found in different chondrite groups and represent some of the oldest known solids of the Solar System. As such, CAIs provide key information regarding the conditions prevailing in the solar protoplanetary disk as well as subsequent mixing and transport processes. Many studies have investigated CAIs for their isotopic compositions and reported nucleosynthetic isotope anomalies in numerous elements, which are typically explained by the variable incorporation of isotopically highly anomalous presolar phases. However, with the exception of 54 Cr-enriched nanospinels, the exact presolar phases responsible for the isotopic heterogeneities are yet to be identified. To address this issue, we here present in-situ Ti isotopic analyses obtained on a diverse set of CAIs from various CV3 chondrites. In this work, the in-situ measurements were performed by targeting individual mineral phases of 15 CAIs with laser-ablation mass spectrometry and indicate significant inter- and intra-CAI isotopic heterogeneity in the neutron-rich isotope 50 Ti. This is particularly pronounced for primitive fine-grained CAIs, whereas coarse-grained CAIs, which have been subject to melting, exhibit smaller degrees of Ti isotopic heterogeneity. To further investigate this Ti isotopic heterogeneity, we additionally obtained Ti isotopic compositions of sequential acid leachates from two fine-grained and two coarse-grained CAIs derived from CV3 chondrites. In contrast to potential expectations from the first part of the study, we do not observe any significant intra-CAI Ti isotopic heterogeneity between the different leaching steps. The lack of intra-CAI Ti isotopic heterogeneity in the acid leachate samples of this study likely reflects that the leaching procedure is unable to efficiently separate the carriers of isotopically anomalous Ti in CAIs. By comparing the bulk CAI Ti isotope compositions with Ti isotope data for hibonite-rich objects from the literature, we find that the range of Ti isotope compositions recorded by CAIs from various chondrite groups can be accounted for by the averaging of hibonite grains. In turn, the variable Ti isotope compositions of hibonite grains can be explained by the averaging of isotopically diverse presolar grains present in the Sun’s parental molecular cloud. This effect of averaging is statistically supported by the central limit theorem, and the concept has the potential to be useful for other isotopic systems.

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

Copper availability governs nitrous oxide accumulation in wetland soils and stream sediments

Denitrification is microbially-mediated through enzymes containing metal cofactors. Laboratory studies of pure cultures have highlighted that the availability of Cu, required for the multicopper enzyme nitrous oxide reductase, can limit N2O reduction. However, in natural aquatic systems, such as wetlands and hyporheic zones in stream beds, the role of Cu in controlling denitrification remains incompletely understood. In this study, we collected soils and sediments from three natural environments -- riparian wetlands, marsh wetlands, and a stream -- to investigate their nitrogen species transformation activity at background Cu levels and different supplemented Cu loadings. All of the systems contained solid-phase associated Cu below or around geological levels (40 - 280 nmol g-1) and exhibited low dissolved Cu (3-50 nM), which made them appropriate sites for evaluating the effect of limited Cu availability on denitrification. In laboratory incubation experiments, high concentrations of N2O accumulated in all microcosms lacking Cu amendment except for one stream sediment sample. With Cu added to provide dissolved concentrations at trace levels (10-300 nM), reduction rate of N2O to N 2 in the wetland soils and stream sediments was enhanced. A kinetic model could account for the trends in nitrogen species by combining the reactions for microbial reduction of NO 3 - to NO 2 -/N 2 O/N 2 and abiotic reduction of NO 2 - to N 2 . The model revealed that the rate of N 2 O to N 2 conversion increased significantly in the presence of Cu. For riparian wetland soils and stream sediments, the kinetic model also suggested that overall denitrification is driven by abiotic reduction of NO2- in the presence of inorganic electron donors. This study demonstrated that natural aquatic systems containing Cu at concentrations less than or equal to crustal abundances may display incomplete reduction of N 2 O to N 2 that would cause N 2 O accumulation and release to the atmosphere.

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

Morphological and crystallographic controls in the replacement of calcite and aragonite by cerussite and otavite

Mineral replacement reactions are essential solid-fluid interactions in natural and industrial processes including metasomatism, diagenesis, and metamorphism, and sequestration of metal toxins from polluted water. In this paper, we explore the morphological evolution of aragonite and calcite (polymorphs of CaCO 3 ) during their replacement by cerussite (PbCO 3 , an isomorph of aragonite) and otavite (CdCO 3 , an isomorph of calcite) in acidic aqueous solutions (with initial pH 3 to 4). Observations using scanning electron and synchrotron-based transmission X-ray microscopies reveal the formation of pseudomorphic shells (~5 to 10 µm thick) of cerussite and otavite when Pb 2+ is reacted with either calcite or aragonite and when Cd 2+ is reacted with aragonite. The formation of pore space at the substrate-precipitate interface and within the precipitate is found to be key in promoting chemical exchange between the dissolving phase and the solution. In contrast, minimal reactivity is observed when Cd 2+ is reacted with calcite. These results demonstrate that the lack of substrate-precipitate epitaxy (either because of inconsistent symmetries or due to large lattice constant mismatch) is a critical factor that enables replacement of primary CaCO 3 minerals by Pb- or Cd-bearing secondary minerals while the presence of epitaxial relationships inhibit those reactions. When calcite and aragonite coexist, the replacement of aragonite by cerussite and otavite occurs preferentially over calcite. We postulate that the relative reactivity of calcite and aragonite to Pb 2+ and Cd 2+ can be determined by interplay between morphological and crystallographic properties of CaCO 3 materials as well as solution chemistry. Our results demonstrate that coupled dissolution-precipitation processes of carbonate minerals play an important role in controlling mobility and bioavailability of lead and cadmium in environments where acidic metal-contaminated water interacts with carbonate-rich rocks and soils.

58 GEOSCIENCES↗
DOE OSTI2023

Position-specific carbon isotopes of Murchison amino acids elucidate extraterrestrial abiotic organic synthesis networks

The Murchison meteorite is a well-studied carbonaceous chondrite with relatively high concentrations of amino acids thought to be endogenous to the meteorite, in part because they are characterized by carbon isotope (δ 13 C) values higher than those typical of terrestrial amino acids. Past studies have proposed that extraterrestrial amino acids in the Murchison meteorite could have formed by Strecker synthesis (for α-amino acids), Michael addition (for β-amino acids), or reductive amination, but a lack of constraints have prevented confident discrimination among these possibilities, or assignment of specific formation pathways to each of several specific amino acids. Position-specific carbon isotope analysis differentiates amongst these mechanisms by relating molecular sites to isotopically distinct carbon sources and by constraining isotope effects associated with elementary chemical reactions. Prior measurements of the position-specific carbon isotopic composition of α-alanine from the Murchison CM chondrite demonstrated that alanine’s high δ 13 C VPDB value is attributable to the amine carbon (δ 13 C VPDB = +142 ± 20‰), consistent with Strecker synthesis drawing on 13 C-rich carbonyl groups in precursors (Chimiak et al., 2021). Here, we measured the δ 13 C composition of fragment ions generated by electron impact ionization of derivatized α-alanine, β-alanine, and aspartic acid from Murchison via gas chromatography-Fourier transform mass spectrometry. α-Alanine’s amine carbon yielded δ 13 C VPDB = +109 ± 21‰, which is consistent with the previously measured value and with formation from 13 C-rich precursors. β-Alanine’s amine carbon presents a lower δ 13 C VPDB = +33 ± 24‰, which supports formation from 13 C-rich precursors but potentially via a Michael addition mechanism rather than Strecker synthesis. Aspartic acid’s amine carbon has δ 13 C VPDB = -14 ± 5‰, suggesting synthesis from precursors distinct from those that generated the alanine isomers. Further, these measurements indicate that Murchison amino acids are a mixture of compounds made from different synthesis mechanisms, though some subsets likely drew on the same substrates; this conclusion highlights the complexity of extraterrestrial organic synthesis networks and the potential of emerging methods of isotope ratio analysis to elucidate the details of those networks.

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

Mineral surface-catalyzed oxidation of Mn(II) by bromate: Implications for the occurrence of Mn oxides on Mars

The occurrence of manganese (Mn) oxides on Mars is believed to be an indicator of an O 2 -rich paleoenvironment of Mars because Mn oxides often form through the oxidation of Mn(II) by O 2 on the surface of Earth. An alternative formation pathway was recently proposed, in which Mn(II) is oxidized by bromate (BrO 3 - ), a common oxidant in contemporary Martian regolith. However, the oxidation of Mn(II) by bromate in solution is kinetically controlled and slow unless using very high concentrations (100 mM) of reactants that may be irrelevant to the conditions of Mars. We conducted laboratory simulations to determine whether iron (Fe) oxides (hematite and goethite) and a phyllosilicate (montmorillonite), abundant minerals on the surface of Mars, could catalyze the oxidation of Mn(II) by bromate. Hematite and goethite, but not montmorillonite, dramatically accelerated the oxidation with a low concentration (1 mM) of Mn(II) and bromate under various solution conditions. Further, the reaction system was autocatalytic with Fe oxides initiating the oxidation of Mn(II) at the early stage and the subsequent catalysis mainly provided by the Mn oxide products. In contrast to producing Mn(IV)O 2 only during the homogeneous oxidation of Mn(II) by bromate in solutions, the heterogeneous mineral-surface catalyzed oxidation resulted in a mixture of Mn(III)OOH and Mn(IV)O 2 phases. Mn(III)OOH was an intermediate product and can be further oxidized by bromate to Mn(IV)O 2 . The occurrence and accumulation of the intermediate product MnOOH can be attributed to its rapid formation due to surface-enhanced nucleation and growth on Fe oxide surfaces and to its higher resistance to oxidation by bromate than Mn(III) ions or clusters. Overall, mineralsurface catalyzed oxidation of Mn(II) by bromate is favorable from both thermodynamic and kinetic perspectives, and can be a major pathway for the occurrence of Mn oxides on Mars where microorganisms are lacking to catalyze the reaction. Our study further improves our understanding of the thermodynamic and kinetic controls on Mn(II) oxidation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Equilibrium Fe isotope fractionation between olivine, pyroxene, spinel and MORB glass: Implications for mantle partial melting to generate MORBs

Primitive mid-ocean ridge basalts (MORBs) exhibit Fe isotopic compositions heavier than the upper mantle by +0.074 ± 0.028 ‰ for δ 56 Fe. The processes responsible for this isotopic difference remain unclear. Modeling of Fe isotope fractionation during mantle partial melting requires reliable equilibrium Fe isotope fractionation factors between minerals and melts, for which consistent data are still lacking. Here, in this study, we used Nuclear Resonant Inelastic X-ray Scattering (NRIXS) technique to measure Fe force constants for a MORB glass (ALV 519-4-1) and natural mantle minerals (olivine, orthopyroxene, clinopyroxene, and spinel) to determine the equilibrium Fe isotope fractionation factors between them. The force constants determined in this study, in increasing order, are 167 ± 26 N/m for spinel, 175 ± 17 N/m for olivine, 176 ± 20 N/m for MORB glass, 205 ± 26 N/m for clinopyroxene, and 219 ± 36 N/m for orthopyroxene. We evaluated the previously proposed mechanisms for the heavy Fe isotopic composition of MORBs, including (i) mantle partial melting, (ii) mantle lithological heterogeneity, with pyroxenite in the source, (iii) mantle metasomatism by low-degree melts, and (iv) fractional crystallization of olivine from melts. For (i), we used the pMELTS program to simulate adiabatic decompression melting of mantle peridotites, and calculated Fe isotope fractionation based on Fe 3+ –Fe 2+ equilibrium-controlled fractionation, where Fe 3+ forms stronger bonds and is more incompatible than Fe 2+ . At 10 wt% peridotite melting, corresponding to MORB generation, only +0.03 ‰ Fe isotope fractionation between the melt and the original bulk composition (Δ 56 Fe = δ 56 Fe melt - δ 56 Fe 0 ) was produced, insufficient to account for the observed MORB-upper mantle difference. For (ii), melting of pyroxenites yields smaller Fe isotope fractionation than melting of peridotites, making it unlikely the cause for the MORB-upper mantle isotopic difference. For (iii), both the Fe 3+ /ΣFe ratio and the δ 56 Fe of melts increase with the degree of partial melting, indicating that low-degree melts are not isotopically heavy enough to significantly alter the isotopic composition of lithospheric mantle through metasomatism. For (iv), equilibrium isotope fractionation between olivine and melt is near zero. These results suggest that equilibrium Fe isotope fractionation alone cannot explain the MORB isotopic signature, highlighting the potential role of kinetic isotope fractionation. Using a diffusion model, we calculated kinetic Fe and Mg isotope fractionations associated with (iv) olivine crystallization from a melt, and found that the predicted Fe and Mg isotope fractionations were inconsistent with observations in MORBs. Qualitatively, two processes could have induced kinetic Fe isotope fractionation during MORB generation: (a) Fe-Mg interdiffusion between melt and solid during melt migration and (b) reactive melt-rock interactions during melt focusing. However, a quantitative understanding of their role in modifying the melt isotopic composition remains limited and requires further investigation.

Fe isotopes↗
DOE OSTI2026

The solubility and speciation of REE phosphate endmembers (CePO 4 and YPO 4 ) in Cl-rich aqueous fluids from 350 to 450 °C and implications for natural systems

The rare earth elements (REE) are important metals used increasingly in advanced technologies. Within the crust, the elements Ce and Y are commonly more abundant compared to other lanthanides and comprise important end-member constituents of REE-bearing minerals. Specifically, Ce is part of the light (L) REE which have larger ionic radii than the heavy (H) REE, which are grouped together with Y. These differences in ionic radius can lead to important physico-chemical trends within the lanthanide group. Despite a recent increase in experimental and thermodynamic data for the REE at high temperature and pressure, there is still a significant lack of these data at supercritical conditions. In this study we conducted batch-type experiments to measure the solubility of REE phosphates (CePO 4 and YPO 4 ) at varying starting pH (1.5–10), and salinity (0.01–1.4 mol/kg NaCl) at 350 and P sat , and from 400 to 450 °C at 700 bar. Results show that the solubility of Ce (33–0.14 ppb) is generally higher than Y (13–0.13 ppb) and that Ce complexes more strongly with both chloride and hydroxyl ligands compared to Y. The solubilities of both REE phosphates are highly pH-dependent and, to a lesser extent, depend on salinity at the studied conditions. The solubility data from this study were implemented into the GEMSFITS program to optimize the thermodynamic properties of Ce and Y hydroxyl and chloride species. The updated standard partial molal Gibbs energies of formation (Δ f G 0 T,P ) are used within the experimental temperature and pressure range to accurately predict the CePO 4 and YPO 4 solubility and Ce and Y speciation behavior. Based on the updated thermodynamic properties we also provide formation constants (log β n Cl,OH ) for Ce and Y hydroxyl and chloride species. Updated thermodynamic properties are applied to model REE-apatite dissolution and REE mobility based on the Pea Ridge iron oxide apatite deposit in Missouri, USA. The apatite dissolution model replicates natural observations including the replacement of monazite and xenotime after apatite and is an example of the utility of the new thermodynamic constants applied to supercritical crustal fluids. Furthermore, the findings of this study advance the predictive capabilities of geochemical models, our understanding of the behavior of individual REE, and permit modeling the overarching fractionation trends between LREE and HREE in supercritical crustal fluids.

58 GEOSCIENCES↗
DOE OSTI2023

Effects of head-starting on multi-year space use and survival of an at-risk tortoise

A major challenge in the recovery of long-lived at-risk taxa like turtles is low juvenile recruitment. Head-starting—the raising of juveniles to larger sizes to improve survival—is one tool that can be used in circumstances where juvenile recruitment is limited. Due to declining populations and difficulty detecting juveniles, however, lack of knowledge of the ecology of juveniles can hinder efforts to develop and evaluate head-starting programs for many turtle species. We sought to inform recovery efforts of Mojave desert tortoises by quantifying multi-year space use and survival of head-started juveniles after release. We radio-tracked tortoises head-started under three different husbandry treatments that varied in rearing duration (from two to over six years) and whether head-starting included an indoor rearing component the first year. We compared postrelease space use and survival as a function of treatment, release size, and time since release. We found that space use, including home range size and site fidelity, varied by husbandry treatment, with smaller and younger tortoises having smaller home ranges and higher site fidelity. Additionally, home range size decreased and site fidelity increased with time since release across treatments. Tortoises with an indoor-rearing component experiencing increased risk of mortality as movement increased compared to tortoises reared solely outdoors. Nevertheless, survival did not differ among treatments or with tortoise age or size. Regardless of husbandry treatment, head-started tortoises exhibited similar space-use and survival overall. Our study provides insight into juvenile tortoise behavior and head-starting as a tool for tortoise conservation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

How to estimate soil organic carbon stocks of agricultural fields? perspectives using ex-ante evaluation

Estimating soil organic carbon (SOC) stocks of agricultural fields has a range of important applications from development of sustainable management practices to monitoring carbon stocks. There are many estimation strategies with the potential for more reliable estimates of SOC stock and more efficient use of soil sampling and analysis resources, especially by leveraging readily available auxiliary information such as remote sensing. However, concrete guidance for strategy selection is lacking. This study narrows this gap with a comparison of strategies for estimating deep SOC stock (0–60 cm) in a prototypical field. Using high density SOC stock measurements and simulation, we built on past studies by 1) ex-ante evaluating a large number of strategy options, 2) using a Bayesian approach to quantify the uncertainty of the comparison, and 3) considering multiple Bayesian models to assess sensitivity to this modeling choice. We found that, using readily available auxiliary information, both balanced and stratified sampling offer substantial improvements over simple random sampling. The auxiliary information most important for this improvement is a Sentinel-2 SOC index = blue / (green × red), followed by the topographic wetness index. We found that these results are robust to the choice of mapping method, but that there is uncertainty in the magnitude of improvement. Here, we recommend future studies implement this Bayesian approach for simulated ex-ante evaluation of SOC stock estimation strategies across more fields to investigate the generalizability of these findings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Uncertainty quantification of the convolutional neural networks on permeability estimation from micro-CT scanned sandstone and carbonate rock images

Rock permeability is one of the most crucial properties affecting subsurface fluid flow behaviors. To accurately and robustly estimate the permeability, Digital Rock Physics, including micro-CT scanning technology and direct flow simulations on scanned images, has prevailed in recent years. Besides, machine learning techniques such as convolutional neural networks (CNNs) have been widely adopted and achieved success in permeability estimations directly from rock images. However, existing ML methods used for permeability estimation from rock images lack uncertainty quantification that causes unreliable predictions and overconfident estimations on out-of-distribution (OOD) samples. Here, in this work, we propose a PI3NN-CNN framework to address this problem. PI3NN-CNN consists of a CNN model for absolute permeability estimation and a PI3NN method to quantify the estimation uncertainty. It is able to quantify the uncertainty for in-distribution (InD) data with a desired confidence level, and identify OOD samples to avoid overconfident predictions. We demonstrate the method using micro-CT scanned images from two sandstone and two carbonate rocks. We found that PI3NN-CNN generates accurate predictions for InD samples, while producing high-quality prediction uncertainties regardless of the prediction accuracy. Meanwhile, PI3NN-CNN identifies OOD samples using its special network initialization scheme. The unique feature of PI3NN-CNN makes it applicable to more complex real-world image-based data for robust learning and predictions without overconfident estimations when the ground-truth information is unavailable.

58 GEOSCIENCES↗
DOE OSTI2023

Stability study of aqueous foams under high-temperature and high-pressure conditions relevant to Enhanced Geothermal Systems (EGS)

In Enhanced Geothermal System (EGS) hydraulic fracturing is carried out by injecting cold water into deep Hot Dry Rocks (HDR) under carefully controlled conditions to create new or reopen existing fractures. Water-based fracturing fluids demonstrate some challenges including immense quantity of water usage, water sensitivity of the formations, water blocking, and lack of proppant carrying capacity and transportation. Thus, an alternative is to use foam-based fracturing fluid which offers potential advantage over conventional water-based fracturing fluid such as minimum water usage, reduced wellbore damage, high proppant carrying capacity, and less environmental damage. However, foams are complex mixture of gaseous phase and liquid phase which are thermodynamically unstable at downhole conditions, and their stability can decrease over time due to liquid drainage, bubble coarsening, and coalescence. This paper shows laboratory experiments executed to study foam stability at high temperature (200°C) and high pressure (6.9MPa) conditions which simulates the geothermal environment. Foam stability was characterized by half-life of foam, which is defined as the time taken by the foam to decreases by 50% of its original height due to drainage. In this paper, two types of gaseous phases, nitrogen (N 2 ) and carbon dioxide (CO 2 ) were investigated. Also, based on successful practice of foam-based fracturing fluid in oil and gas industries, four surfactants, including Alpha olefin sulfonate (AOS), Sodium dodecyl sulfonate (SDS), Tergitol™ (NP – 40), and Cetyltrimethylammonium chloride (CTAC) at optimum concentration of 1wt.% were tested for best stability performance. In addition, different stabilizing agents including guar gum, bentonite clay, crosslinker, silicon dioxide nanoparticles (SiO 2 ), graphene oxide (GO) were also studied. Experimental results showed that N 2 foams were more stable than CO 2 foams. It was observed that foam half-life decreased with the increase in temperature. Among all the surfactants, AOS foams showed the most promising thermal stability at high temperatures. Moreover, with the addition of stabilizing agents, foam's half-life was enhanced. Stabilizing agents such as crosslinker and GO dispersion showed the most stable foams with half-life recorded at 20min and 17min, respectively, at 200°C and 6.9MPa. Finally, pressure also showed a positive effect on foam stability; with increased pressure, foam half-life was increased. Based on the experimental data, analytical models for the effect of temperature and pressure were developed, considering foam degradation is a first-order kinetic reaction that linearly depends on the foam drainage mechanism. The effect of temperature on foam half-life was studied as an exponential decay model. In this model, foam half-life is a function of drainage rate constant (D A ) and activation energy (E a ) of the foam system. The effect of pressure on foam half-life was found to obey a power-law model where an increase in pressure showed an increase in foam half-life. Furthermore, a linear relation was studied for the effect of pressure on foam activation energy and drainage rate. Then the combined effects of temperature and pressure were studied, which yielded an analytical model to predict the foam stabilities in terms of half-life for different foam compositions. In conclusion, this research indicates that with an appropriate selection of surfactants and stabilizing agents, it is possible to obtain stable foams, which could replace conventional water fracturing fluid under EGS conditions.

15 GEOTHERMAL ENERGY↗