Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geochemical Modeling”

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

NEWTS USGS Brackish Water Database

Data from the USGS Brackish Water Database on brackish water compositions from across the US. Data has been filtered for high salinity and completeness of data attributes. Data has been transposed and converted for ease of input into aqueous chemistry modeling software including OLI Systems and Geochemist's Workbench. Original data from: Qi, S.L., and Harris, A.C., 2017, Geochemical Database for the Brackish Groundwater Assessment of the United States: U.S. Geological Survey data release, https://doi.org/10.5066/F72F7KK1.

Brackish Groundwater↗
DOE OSTI2020

Reactive particle-tracking solutions to a benchmark problem on heavy metal cycling in lake sediments

Geochemical systems are known to exhibit highly variable spatiotemporal behavior. This may be observed both in non-smooth concentration curves in space for a single sampling time and also in variability between samples taken from the same location at different times. However, most models that are designed to simulate these systems provide only single-solution smooth curves and fail to capture the noise and variability seen in the data. We apply a recently developed reactive particle-tracking method to a system that displays highly complex geochemical behavior. When the method is made to most closely resemble a corresponding Eulerian method, in its unperturbed form, we see near-exact match between solutions of the two models. More importantly, we consider two approaches for perturbing the model and find that the spatially-perturbed condition is able to capture a greater degree of the variability present in the data. This method of perturbation is a task to which particle methods are uniquely suited and Eulerian models are not well-suited. Additionally, because of the nature of the algorithm, noisy spatial gradients can be highly resolved by a large number of mobile particles, and this incurs negligible computational cost, as compared to expensive chemistry calculations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Application of an Embedded Fracture and Borehole Modeling Approach to the Understanding of EGS Collab Experiment 1

EGS Collab is a series of meso-scale experiments and associated numerical simulation activities being funded by the United States Department of Energy, Geothermal Technologies Office (GTO) to investigate enhanced geothermal system processes under in-situ stress and slightly elevated temperature conditions. This project is designed to provide scientists and engineers with immediate access to impermeable rock at scales larger than possible in the laboratory, but generally smaller than those for commercial production. Immediate access to rock is provided via the existing drifts of the former Homestake Gold Mine, now operated as the Sanford Underground Research Facility in Lead, South Dakota. The objectives of the EGS Collab project are to develop well controlled fracture networks between injector and producer boreholes using normal and shear stimulation for permeability enhancement. The first experimental site was located off the West Access drift on the 4850 Level (4850 feet below ground surface) in phyllite of the Precambrian Poorman formation, and involved the creation of a fracture network comprising a combination of hydraulic and natural fractures. A second experimental site is now being considered near the battery alcove on the 4100 Level in amphibolite and rhyolite of the Yates Unit. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO’s flagship EGS research effort. Experiment 1 has comprised a series of successful tests, including a long-term chilled-water circulation test, but the testbed has two atypical EGS elements. Active ventilation in the adjacent drift over a 50-year period cooled the testbed rock mass and hydraulic fracturing intersected monitoring boreholes making them conduits for fluid flow. Numerical simulations executed in support of the design of the EGS Collab Experiment 1, computed a radial temperature distribution and stress gradient orthogonal to the drift axis, resulting in the forecast of an oblong hydraulic fracture geometry, extended in the direction of the drift from the stimulation borehole. This paper describes the numerical simulation of long-term chilled-water circulation test with an embedded fracture and borehole modeling approach. The principal objective of the simulation work is show agreement between experimental observations in terms of production fluid temperatures, monitoring borehole temperatures, injection and production pressures, and tracer recoveries using a single conceptual model for the fracture network, using characterization and monitoring data generally available for EGS. A secondary objective is to improve the numerical simulation result comparisons with additional data available from the broader monitoring equipment within the testbed versus deeper, hotter, and more remote EGS. One complicating factor for comparisons between the numerical simulation results and experimental observations, is the Joule-Thomson heating associated with large pressure drops across the fracture network.

numerical simulation, embedded fractures and boreh↗
DOE OSTI2022

On a Unified Core Characterization Methodology to Support the Systematic Assessment of Rare Earth Elements and Critical Minerals Bearing Unconventional Carbon Ores and Sedimentary Strata

A significant gap exists in our understanding and ability to predict the spatial occurrence and extent of rare earth elements (REE) and certain critical minerals (CM) in sedimentary strata. This is largely due to a lack of existing, systematic, and well-distributed REE and CM samples and analyses in United States sedimentary basins. In addition, the type of sampling and characterization performed to date has generally lacked the resolution and approach required to constrain geologic and geographic heterogeneities typical of subsurface, mineral resources. Here, we describe a robust and systematic method for collecting core scale characterization data that can be applied to studies on the contextual and spatial attributes, the geologic history, and lithostratigraphy of sedimentary basins. The methods were developed using drilled cores from coal bearing sedimentary strata in the Powder River Basin, Wyoming (PRB). The goal of this effort is to create a unified core characterization methodology to guide systematic collection of key data to achieve a foundation of spatially and geologically constrained REEs and CMs. This guidance covers a range of measurement types and methods that are each useful either individually or in combination to support characterization and delineation of REE and CM occurrences. The methods herein, whether used in part or in full, establish a framework to guide consistent acquisition of geological, geochemical, and geospatial datasets that are key to assessing and validating REE and CM occurrences from geologic sources to support future exploration, assessment, and techno-economic related models and analyses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Porewater Geochemical Assessment of Seismic Indications for Gas Hydrate Presence and Absence: Mahia Slope, East of New Zealand’s North Island

We compare sediment vertical methane flux off the Mahia Peninsula, on the Hikurangi Margin, east of New Zealand’s North Island, with a combination of geochemical, multichannel seismic and sub-bottom profiler data. Stable carbon isotope data provided an overview of methane contributions to shallow sediment carbon pools. Methane varied considerably in concentration and vertical flux across stations in close proximities. At two Mahia transects, methane profiles correlated well with integrated seismic and TOPAS data for predicting vertical methane migration rates from deep to shallow sediment. However, at our “control site”, where no seismic blanking or indications of vertical gas migration were observed, geochemical data were similar to the two Mahia transect lines. This apparent mismatch between seismic and geochemistry data suggests a potential to underestimate gas hydrate volumes based on standard seismic data interpretations. To accurately assess global gas hydrate deposits, multiple approaches for initial assessment, e.g., seismic data interpretation, heatflow profiling and controlled-source electromagnetics, should be compared to geochemical sediment and porewater profiles. A more thorough data matrix will provide better accuracy in gas hydrate volume for modeling climate change and potential available energy content.

03 NATURAL GAS↗
DOE OSTI2024

Rare earth element enrichment in coal and coal-adjacent strata of the Uinta Region, Utah and Colorado

This study aims to quantify rare earth element enrichment within coal and coal-adjacent strata in the Uinta Region of central Utah and western Colorado. Rare earth elements are a subset of critical minerals as defined by the U.S. Geological Survey. These elements are used for a wide variety of applications, including renewable energy technology in the transition toward carbon-neutral energy. While rare earth element enrichment has been associated with Appalachian coals, there has been a more limited evaluation of western U.S. coals. Here, samples from six active mines, four idle/historical mines, four mine waste piles, and seven stratigraphically complete cores within the Uinta Region were geochemically evaluated using portable X-ray fluorescence ( n = 3,113) and inductively coupled plasma-mass spectrometry ( n = 145) elemental analytical methods. Results suggest that 24%–45% of stratigraphically coal-adjacent carbonaceous shale and siltstone units show rare earth element enrichment (>200 ppm), as do 100% of sampled igneous material. A small subset (5%–8%) of coal samples display rare earth element enrichment, specifically in cases containing volcanic ash. This study proposes two multi-step depositional and diagenetic models to explain the enrichment process, requiring the emplacement and mobilization of rare earth element source material due to hydrothermal and other external influences. Historical geochemical evaluations of Uinta Region coal and coal-adjacent data are sparse, emphasizing the statistical significance of this research. These results support the utilization of active mines and coal processing waste piles for the future of domestic rare earth element extraction, offering economic and environmental solutions to pressing global demands.

Coe, Haley H.↗
DOE OSTI2024

Rare Earth Elements (REE) and Critical Minerals (CM) in Middle Pennsylvanian-Age Coals and Associated Sediments in the Central Appalachian Basin, Eastern U.S.A.

The Central Appalachian Basin (CAB) of Kentucky, Tennessee, West Virginia and Virginia has a long history of coal mining and oil and gas extraction that has empowered the regional and national economies, the development of infrastructure, and a highly trained energy resources work force. As our societal demands for advanced technologies have rapidly increased in recent years, coal-related materials are viewed as an important new unconventional domestic source of critical minerals (CM) that are required for telecommunications, aerospace and transportation industries, electronics, the transition to low-carbon emissions energy production, and many consumer products. Coal-related materials encompass coal, associated sediments, coal mining waste materials, produced waters, and ash residues from coal-fired power plants. An important objective of the Evolve Central Appalachia Project (Evolve CAPP), sponsored by the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL) is to assess the quantity and distribution of CM resources in the CAB region. The rare earth elements (REE) are considered highest priority, although other important CM such as niobium, gallium, and zirconium are known to occur in the Middle Pennsylvanian-age coals and sediments. Working with coal industry partners who provided access to drill cores, coal-related sediments and waste materials, over 600 samples have been collected for laboratory analysis, and over 730 materials have been scanned using portable x-ray fluorescence (pXRF) equipment. The application of pXRF provides the means for real-time semiquantitative analysis of CM content at very close spacings (typically 2-3 inch intervals) along drill core and in-situ channel samples that span the roof, coal seam, and floor rock. The comparison of pXRF geochemical data with laboratory results, geologic data, and downhole spectral gamma logs can provide high resolution input to lithologic and depositional models for future CM resource evaluations. The preliminary findings show that pXRF is capable of accurately measuring low concentrations of many of the CM with a high level of confidence (Ba, Cr, Ga, K, Nb, Rb, Sr, Th, Y), whereas for others (La, Ce, Co, Mn, Nd, Ni, Sc, Ti, V) the detection limits are very high or spectral interferences increase the uncertainty. Notably, the mean abundances of Y (34 ppm), La (86 ppm), Ga (40 ppm), and V (146 ppm) in coal underclays in the CAB region are up to 7X enriched compared with the overlying coal. These values also exceed the reported concentrations in published reference materials for upper continental crustal rocks (Rudnick and Gao, 2003), North American Shale Composite (Gromet et al., 1984), and North American coal (Finkelman, 1993). The pXRF data are in part verified by laboratory results that indicate the mean Y abundance (36 ppm) is highly correlated (R2 = 0.807) with total REE (ΣREE). The correlation is even higher (R2 = 0.957) with heavy REE (ΣHREE). Applying these correlations to the pXRF data for the coal underclays, the mean estimated values for ΣREE+Y and ΣHREE+Y are 270 ppm and 58 ppm, respectively. Although these average values are not considered high, the range of Y measured by pXRF in the coal underclays extended as high as 114 ppm, which would suggest ΣREE+Y equal to 847 ppm. The mean abundances of Zr (192 ppm) and Th (21 ppm) in coal underclays are also enriched compared with the overlying coal and these results likely reflect the presence of resistant detrital heavy minerals such as monazite, xenotime, and zircon in the underclay matrix. Several of the profiled coal seams and associated wall rocks contained thin volcanic ash layers up to 4-5 inches in thickness. The extent to which these ash fall layers provided a source for CM under the paleoenvironmental conditions that resulted in coal deposits in the CAB remains to be fully studied. Continuing investigations in the Evolve CAPP study area will include laboratory determinations of mineralogic and clay compositions, and evaluations of CM geochemical mobility in the coal and coal underclays.

Lassetter, Billy↗
DOE OSTI2026

Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0

Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.

Critical Minerals↗
DOE OSTI2025

Planetesimal Impact Vapor Plumes and Nebular Shocks Form Chondritic Mixtures

The origin of chondrules and the chondritic sedimentary rocks that dominate the meteoritic record is a long-standing problem in planetary science. Here, we develop a physical model for the formation of chondritic mixtures as an outcome of vaporizing collisions between planetesimals that were dynamically excited by the growth and migration of planets. We present calculations of nebular shock waves generated by impact vapor plumes and focus on aspects of the plume interaction with the nebular gas and dust that have been neglected in previous studies of impact ejecta. We find that when water dominates the vapor, the plumes are relatively cool. However, the plume expansion is supersonic and can drive strong shock waves in the dusty nebular gas. Portions of these nebular shock fronts initially melt nebular dust, forming chondrules that are coupled to the moving front. As the shock front expands and cools, the chondrules solidify, while the shock front entrains additional dust. Eventually, the plume expansion stalls and then hydrodynamically collapses, turbulently mixing variably processed dust and size-sorted chondrules. For probable impact parameters and nebular conditions during giant planet growth and migration, the impact-generated mixtures have characteristics that span the range observed in chondritic meteorites, providing an environment for rapid formation of chondritic assemblages after chondrule formation. Our impact vapor and nebular shocks model, named the IVANS model, links chondrule formation to the overall context of planet formation and provides a framework for interpreting the detailed chronological and geochemical record contained in chondritic meteorites.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

DEEPEN 3D PFA Index Models for Exploration Datasets at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the development of the DEEPEN 3D play fairway analysis (PFA) methodology for magmatic plays (conventional hydrothermal, superhot EGS, and supercritical), index models needed to be developed to map values in geoscientific exploration datasets to favorability index values. This GDR submission includes those index models. Index models were created by binning values in exploration datasets into chunks based on their favorability, and then applying a number between 0 and 5 to each chunk, where 0 represents very unfavorable data values and 5 represents very favorable data values. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were created using histograms of the distributions of each exploration dataset in combination with literature and input from experts about what combinations of geophysical, geological, and geochemical signatures are considered favorable at Newberry. This is in attempt to create similar sized bins based on the current understanding of how different anomalies map to favorable areas for the different types of geothermal plays (i.e., conventional hydrothermal, superhot EGS, and supercritical). For example, an area of partial melt would likely appear as an area of low density, high conductivity, low vp, and high vp/vs. This means that these target anomalies would be given high (4 or 5) index values for the purpose of imaging the heat source. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were produced for the following datasets: - Geologic model - Alteration model - vp/vs - vp - vs - Temperature model - Seismicity (density*magnitude) - Density - Resistivity - Fault distance - Earthquake cutoff depth model

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

An integrated experimental–modeling approach to identify key processes for carbon mineralization in fractured mafic and ultramafic rocks

Abstract Controlling atmospheric warming requires immediate reduction of carbon dioxide (CO2) emissions, as well as the active removal and sequestration of CO2 from current point sources. One promising proposed strategy to reduce atmospheric CO2 levels is geologic carbon sequestration (GCS), where CO2 is injected into the subsurface and reacts with the formation to precipitate carbonate minerals. Rapid mineralization has recently been reported for field tests in mafic and ultramafic rocks. However, unlike saline aquifers and depleted oil and gas reservoirs historically considered for GCS, these formations can have extremely low porosities and permeabilities, limiting storage volumes and reactive mineral surfaces to the preexisting fracture network. As a result, coupling between geochemical interactions and the fracture network evolution is a critical component of long-term, sustainable carbon storage. In this paper, we summarize recent advances in integrating experimental and modeling approaches to determine the first-order processes for carbon mineralization in a fractured mafic/ultramafic rock system. We observe the critical role of fracture aperture, flow, and surface characteristics in controlling the quantity, identity, and morphology of secondary precipitates and present where the influence of these factors can be reflected in newly developed thermo-hydro-mechanical–chemical models. Our findings provide a roadmap for future work on carbon mineralization, as we present the most important system components and key challenges that we are overcoming to enable GCS in mafic and ultramafic rocks.

58 GEOSCIENCES↗
DOE OSTI2024

Water foraging with dynamic roots in E3SM; The role of roots in terrestrial ecosystem memory on intermediate timescales (Final Technical Report)

Terrestrial ecosystems can show sustained responses to stress events that may last for years after the initial perturbation. This phenomenon is referred to as legacy or memory and it emerges from a set of ecosystem processes that are poorly captured by Earth System and Land Surface Models. Consequently, these models struggle to predict the impact that extreme climate events have on surface energy, carbon and hydrological exchange once a stressor such as drought has relaxed or in response to repeated exposure to stress. In this project, we set out to test how the addition of dynamic root profiles in land surface models impact the capacity for Earth System Models, such as the Department of Energy’s E3SM, to capture realistic legacy effects. Observations have shown that vegetation shifts their root profiles during stress events to forage for water and these altered root profiles may take years to relax back to the initial state. We hypothesized that the alteration of the belowground root structures may be a source of terrestrial ecosystem legacy that is missing from models. To test this idea we undertook three core activities. (1) We developed a global-scale analysis of the impact of dynamic roots on ecosystem legacy building from a recently developed dynamic root module. (2) We developed new root dynamics in the DOE’s Energy Land Model (ELM) that allowed not only root profiles to shift their profile but also simultaneously alter carbon allocation to fine root pools. We then tested the impact of these new dynamics with intensive sensitivity analysis at four long-term AmeriFlux sites. (3) We undertook detailed isotopic analysis of tree rings from these AmeriFlux sites to assess ecophysiological and ecohydrological legacy to stress events that can be used to benchmark the sensitivity experiments with ELM. From these core activities, we report the following key findings. Firstly, on a global scale, the addition of dynamic roots led to chronically water-stressed ecosystems recovering faster to climate stress while wetter ecosystem showed enhanced legacy. This is because across ecosystems, stress events were almost universally associated with water shortages that led to the development of deeper root profiles. These deeper root profiles proved beneficial for recovery from drought stress. While the root dynamics did not universally improve the modeled representation of legacy it showed complex transient dynamics that emerge from the addition of dynamic roots. Secondly, the sensitivity analysis illustrated long term shifts in rooting depths away from the prescribed default profile suggesting that initializing of root profiles could benefit from spin-up simulations that converge on locally optimized root profiles. In addition, by enabling dynamic allocation some of the sites predicted unrealistically low allocation to roots (and high allocation to leaves). While these changes did not dramatically alter modeled gross primary productivity, they illustrated that without more sophisticated root processes in models, there is little penalty to dramatically disinvest in roots. Thirdly, the isotopic analysis of tree rings showed highly distinct legacy responses across species and sites. For example, T. canadensis showed reduced transpiration the year after stress events illustrated by sustained elevated $\delta^{18}$O. In contrast, A saccharum displayed elevated $\delta^{13}$C associated with reduced stomatal conductance in response to the previous years’ stress event. These geochemical signatures of legacy were present despite tree growth returning to normal the year after the stress. The results show the importance of species-level dynamics in legacy that are absent in modeling that assumes common traits within plant functional types. In summary, the work here established new avenues to explore root dynamics in models while illustrating how additional root processes are needed before dynamic carbon allocation can be implemented. Lastly, this project provided mentorship to a postdoctoral fellow, training for an early career scientist and multiple undergraduate students recruited from a minority serving institution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Grain detachment and transport clogging during mineral dissolution in carbonate rocks with permeable grain boundaries

In this paper, we explore the effect of permeable grain boundaries on carbonate rock dissolution in order to improve our understanding of grain detachment and migration during reactive flow. To do so we investigated the effects of grain size and enhanced reactivity along grain boundaries on global dissolution kinetics. Variations of permeability and porosity were calculated, and exponential relationships were observed in fractured rocks. Our model employed a reactive transport framework based on the Darcy-Brinkman equation to simulate calcite dissolution in carbonate rocks composed of microporous grains and, included fluid transport along grain boundaries. The model includes fluid flow, solute transport by advection–diffusion, heterogeneous reaction between fluid and minerals and grain detachment with subsequent grain transport in macropores. The migration of solid particles due to dissolution was based on cluster analysis and local movement, and the results show that grain detachment can lead to significant decreases in permeability due to the clogging of transport pathways. Microporous media with smaller grain sizes (fine grains in this study) showed a higher average reaction rate than those with coarser grains. In addition, as the overall rates of geochemical processes are commonly affected by the presence of texture heterogeneities such as fractures, a single fracture (macropore) introduced into a microporous matrix composed of permeable grains and grain boundaries was modelled as a large channel connecting the inlet and outlet of the simulation domain. It was found that macropore clogging by grain detachment temporary decreased permeability, lowering the long-term global reaction rates. Lastly, we observed that increases in flow rate can reduce detachment and fracture clogging by reducing local dissolution along grain boundaries.

58 GEOSCIENCES↗
DOE OSTI2022

Quantifying pH buffering capacity in acidic, organic-rich Arctic soils: Measurable proxies and implications for soil carbon degradation

Dynamic pH change promoted by biogeochemical reactions in Arctic tundra soils can be a major control on the production and release of CO 2 and CH 4 , which contribute to rising global temperatures. Large quantities of soil organic matter (SOM) in these soils are susceptible to microbial decomposition, leading to pH changes during permafrost thaw. Soil pH buffering capacity (β) modulates the extent of pH change but has not been thoroughly studied and represented in predictive ecosystem scale biogeochemical models in Arctic tundra soils. In this study, we generated titration curves for 21 acidic tundra soils from three Arctic sites across northern Alaska, United States of America. Geochemical and hydrological soil properties were evaluated, and correlations with β were developed. Strong correlations between β and both gravimetric water content (Θ g ) (R 2 = 0.847, p < 0.001) and soil water retention (SWR) (R 2 = 0.849, p = 0.001) indicate that the ability of soil to retain water could be associated with its buffering properties. Correlations between β and soil organic carbon (SOC) and cation exchange capacity (CEC) were also explored, and relationships to SWR are discussed. These correlations were then used with existing soil databases reporting SOC, CEC, and SWR to estimate β across Alaska soils. We further demonstrated the quantitative relationships between β and the simulated rates of biogeochemical reactions and show that lower β leads to higher soil pH and more CH 4 production. Overall, our study provides simple proxies for β in Arctic soils and highlights the importance and implications of representing soil buffering in predictive models, thereby enabling quantitative coupling between pH dynamics associated with biogeochemical reactions. Integrating β into predictive models of Arctic biogeochemical cycling may reduce model uncertainty and further our understanding of permafrost SOM degradation accelerated by warming.

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

Exploring an ‘ideal hill': how lithology and transport mechanisms affect the possibility of a steady state during weathering and erosion

We introduce a model of chemical reaction within hills to explore how evolving porosity (and by inference, permeability) affects flow pathways and weathering. The model consists of hydrologic and reactive-transport equations that describe alteration of ferrous minerals and feldspar. These reactions were chosen because previous work emphasized that oxygen- and acid-driven weathering affects porosity differently in felsic and mafic rocks. A parameter controlling the order of the fronts is presented. In the absence of erosion, the two reaction fronts move at different velocities and the relative depths depend on geochemical conditions and starting composition. In turn, the fronts and associated changes in porosity drastically affect both the vertical and lateral velocities of water flow. For these instances, estimates of weathering advance rates based on simple models that posit unidirectional constant-velocity advection do not apply. In the model hills, weathering advance rates diminish with time as the Darcy velocities decrease with depth. The system can thus attain a dynamical steady state at any erosion rate where the regolith thickness is constant in time and velocities of both fronts become equal to one another and to the erosion rate. The slower the advection velocities in a system, the faster it attains a steady state. For example, a massive rock with relatively fast-dissolving minerals such as diabase -- where solute transport across the reaction front mainly occurs by diffusion -- can reach a steady state more quickly than granitoid rocks in which advection contributes to solute transport. The attainment of a steady state is controlled by coupling between weathering and hydrologic processes that force water to flow horizontally above reaction fronts where permeability changes rapidly with depth and acts as a partial barrier to fluid flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Surface Complexation/Ion Exchange Hybrid Model for Radionuclide Sorption to Clay Minerals (M4SF-23LL010301062)

This progress report (Level 4 Milestone Number M4SF-23LL010301062) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Argillite International Collaborations Activity Number SF-23LL01030106. The activity is focused on our long-term commitment to engaging our partners in international nuclear waste repository research. The focus of this milestone is the establishment of international collaborations for surface complexation modeling and the associated impacts of unlocking larger, community-based datasets. More specifically, we are developing a database framework for Spent Fuel and Waste and Science Technology (SFWST) that is aligned with the Helmholtz Zentrum Dresden Rossendorf (HZDR) sorption database development group in support of the database needs of the SFWST program. In our FY22 effort, we described a detailed analysis of U(VI) sorption to quartz through both traditional surface complexation modeling and through a hybrid ML framework. In FY23, effort was placed on publication of these results and expansion of the LLNL surface complexation and ion exchange database (L-SCIE) in order to assess mineral-based radionuclide retardation under a wider variety of geochemical conditions (e.g., ionic strength, varying electrolyte compositions). Efforts were initiated to expand L-SCIE to include radionuclide surface complexation and ion exchange to clays that are relevant to subsurface geochemical processes occurring at nuclear waste repositories. In particular, a large source of sorption data for clays resides at the Paul Scherrer Institute (PSI) (work primarily by Bradbury and Baeyens) and we initiated discussions on how to retrieve those data and apply FAIR principles to those datasets. In addition to L-SCIE development, two hybrid models that incorporate AI/ML were investigated and compared to discern the most promising approaches for accurate and precise estimations of radionuclide retardation. Key considerations for future model development include (1) the ability to reduce computational burden on determining retardation coefficients for PA and (2) the ability to quantify and predict radionuclide-mineral partitioning at a more efficient, rapid pace due to automated workflows. Upon the careful consideration of the most effective modeling approaches, we are identifying ways to implement these approaches into PA. Ultimately, the data science-based workflows will provide a major incentive for other institutions to adopt a FAIR-formatted, interoperable database. LLNL will play a key role in disseminating sorption data and acting as good data stewards by updating the database in a consistent format and assessing the quality of the newly assimilated data in an organized fashion. To this end, all data and workflows are open access and made available on the LLNL Seaborg research website (https://seaborg.llnl.gov/resources/geochemical-databases-modeling-codes).

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Multistage nucleation pathway in LiF molten salt mirrors the crystal–melt interface structure

Despite over a century of studies, fundamental questions remain about the processes governing crystal nucleation from melts or solutions. Research over the past three decades has presented mounting evidence for kinetic pathways of crystal nucleation that are more complex than envisioned by the simplest forms of classical theory. Such observations have been presented for colloidal and elemental systems with covalent and metallic bonding. Despite the technological and geochemical importance of molten salts, similar studies for these ionically bonded systems are currently lacking. Here we develop a machine learning interatomic potential for a model ionic system: LiF. The potential features quantum-level accuracy for both liquid and multiple solid polymorphs over wide temperature and pressure ranges and accurately reproduces experimentally measured properties. Thanks to the efficiency of the potential, which enables microsecond-scale molecular dynamics simulations, induction times for nucleation of LiF solids from their melts are computed over a range of undercoolings. With the aid of a set of robust local order parameters established here, the simulations reveal that homogeneous crystal nucleation in undercooled melts preferentially initiates from liquid regions showing slow dynamics and high bond orientational order simultaneously, and the second-shell order of both precritical nuclei and the surface of postcritical nuclei is dominated by hexagonal close packing and body-centered cubic local structure, even though the nucleus core is dominated by face-centered cubic structure corresponding to the stable rocksalt crystal structure. Finally, we establish a connection between the crystallization pathway and the equilibrium crystal-melt interface structure.

Applied Physical Sciences↗
DOE OSTI2022

A Brief Overview of Reactive Transport Codes Used in CO 2 Applications

Reactive transport modelling consists of computational and numerical models that describe the coupled physical, chemical, mechanical, and biological processes interacting with each other over a broad range of spatial and temporal scales. This report reviews some of the reactive transport codes available in the literature for carbon dioxide (CO 2 ) applications and can assist the scientific community with the applicability of each code to solve a specific problem. Approaches implemented in these codes can be categorized in three groups: (i) continuum-scale, (ii) pore-scale, and (iii) hybrid-scale approaches. Theoretical foundations, numerical implementations, and application examples using the models are described. In the continuum approach, flow and transport are formulated in terms of a representative control volume of the medium and coupled with geochemical reactions. The medium is characterized by bulk parameters such as porosity, permeability, or reactive surface area. In the pore-scale approach, each point of space in the pore network is occupied by either a fluid or solid phase. The pore-scale approach requires an exact knowledge of the spatial and temporal phase distribution. In the hybrid-scale approach, media are described by multiple characteristic length-scales, with some regions using a pore-scale approach while others used a continuum approach.

01 COAL, LIGNITE, AND PEAT↗