Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsurface system”

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 613 records · Page 34

Quartz Dissolution Effects on Flow Channelization and Transport Behavior in Three‐Dimensional Fracture Networks

We perform a set of reactive transport simulations in three-dimensional fracture networks to characterize the impact of geochemical reactions on flow channelization. Flow channelization, a frequently observed phenomenon in porous and fractured subsurface rock formations, results from the spatially variable hydraulic resistance offered by a geological structure. In addition to geo-structural features such as network connectivity, geometry, and hydraulic resistance, geochemical reactions, for example, dissolution and precipitation, can dynamically inhibit or enhance flow channelization. These geochemical processes can change the fracture permeability leading to increased flow channelization, which are localized connected regions of high volumetric flow rates that are seemingly ubiquitous in the subsurface. In our simulations, fractures partially filled with quartz are gradually dissolved until quasi-steady state conditions are obtained. We compare the flow field's initial unreacted and final dissolved states in terms of flow and transport observations. We observe that the dissolved fracture networks provide less resistance to flow and exhibit increased flow channelization when compared to their unreacted counterparts. However, there is substantial variability in the magnitude of these changes which implies that the channelization strongly depends on the network structure. In turn, we identify the interplay between the particular network structure and the impact of geochemical dissolution on flow channelization. The presented results indicate that geological systems that have been weathering or reactive for longer times in older landscapes are likely to have increased flow channelization compared to their equivalent but younger counterparts, which implies a time dependence on flow channelization in fractured media.

Hyman, Jeffrey D.↗

Model Package Report: Geoframework Model of the Hanford Site 100 Area

The purpose of the 100 Area Geoframework Model (GFM) is to provide a reasonable, consistent, and defensible three dimensional representation of the hydrostratigraphic units below the River Corridor at the Hanford Site to support contaminant fate and transport models. The GFM is a three dimensional representation of the subsurface geologic structure. From this, three dimensional geologic model-exported results, in the form of points or surfaces, are used as inputs to populate and assemble the various numerical model architectures. The objective of this report is to define the process used to produce a hydrostratigraphic model for the hydrostratigraphic units beneath the Hanford Site 100 Area. The GFM may support several other CH2M HILL Plateau Remediation Company project needs and objectives, including providing geologic information to support remedial investigations and actions and as a tool to present River Corridor geology and contaminant extents. The GFM was constructed based on information through 2019 available in the Integrated Data Management System, Hanford Environmental Information System, and the Hanford Site geologic contacts (GeoContacts) datasets. Revisions to the 100 Area GFM may occur to incorporate new data and information. Each version will be maintained in configuration control using date stamps for identifying supporting databases, figures, and interpretations. Supporting data include the final three dimensional geoframework surfaces, two dimensional structure and isopach maps, and all the geologic contact inputs and interpreted geologic log data. These final products provide a supporting set of geologic information that together allow the creation of the GFM. This information is managed, updated, and maintained via CH2M HILL Plateau Remediation Company under the Environmental Model Management Archive.

58 GEOSCIENCES↗

An ensemble of 48 physically perturbed model estimates of the 1/8° terrestrial water budget over the conterminous United States, 1980–2015

Terrestrial water budget (TWB) data over large domains are of high interest for various hydrological applications. Spatiotemporally continuous and physically consistent estimations of TWB rely on land surface models (LSMs). As an augmentation of the operational North American Land Data Assimilation System Phase 2 (NLDAS-2) four-LSM ensemble, this paper describes a dataset simulated from an ensemble of 48 physics configurations of the Noah LSM with multi-physics options (Noah-MP). The 48 Noah-MP physics configurations are selected to give a representative cross-section of commonly used LSMs for parameterizing runoff, atmospheric surface layer turbulence, soil moisture limitation on photosynthesis, and stomatal conductance. The dataset spans from 1980 to 2015 over the conterminous United States (CONUS) at a monthly temporal resolution and a 1/8° spatial resolution. The dataset variables include total evapotranspiration and its constituents (canopy evaporation, soil evaporation, and transpiration), runoff (the surface and subsurface components), as well as terrestrial water storage (snow water equivalent, four-layer soil water content from the surface down to 2 m, and the groundwater storage anomaly). The dataset is available at https://doi.org/10.5281/zenodo.7109816. Evaluations carried out in this study and previous investigations show that the ensemble per forms well in reproducing the observed terrestrial water storage, snow water equivalent, soil moisture, and runoff. Noah-MP complements the NLDAS models well, and adding Noah-MP consistently improves the NLDAS es timations of the above variables in most areas of CONUS. Besides, the perturbed-physics ensemble facilitates the identification of model deficiencies. The parameterizations of shallow snow, spatially varying groundwater dynamics, and near-surface atmospheric turbulence should be improved in future model versions.

54 ENVIRONMENTAL SCIENCES↗

Geophysical Impacts and Spectroscopic Identification of a Hydrous Iron Sulfate on Icy Worlds

Over geologic time-scales, large volumes of exogenic sulfur ions from Io's plasma torus have been supplied to the surface of Europa and Ganymede, which, combined with recent interpretations of orbiter images, dynamical modeling, and surface-subsurface exchange, suggests further sulfur transport into the interior of the icy worlds. These observations motivate mixed-phase spectral modeling for interpreting orbiter spectroscopy data and determination of hydration states of candidate surface materials including hydrous sulfates. In this work, we present a combined experimental and theoretical study of the low temperature and high pressure vibrational spectral signature of the iron-sulfate monohydrate endmember, szomolnokite (FeSO 4 ·H 2 O). By employing synchrotron Fourier-transform infrared spectroscopy (FTIR) in the diamond anvil cell up to 23 GPa and down to 20 K, we explore the extreme range of pressure-temperature domains relevant to icy environments throughout our solar system and beyond. Combined with our density-functional theory quantum-mechanics molecular dynamics results, we demonstrate that experimentally observed infrared features in the O-H stretching region commonly associated with nH 2 O (n > 1) hydration states can be attributed to a pure monohydrate without the need for pressure-induced exsolved ice, other coexisting hydrous iron sulfates, or strong overtone and combination modes. We further discuss the possibility of lateral variations in density and shear properties on icy worlds associated with temperature variations and the high-pressure phases of kieserite group monohydrated sulfates.

Geosciences↗

Revealing complex subsurface dynamics with continuous seismic monitoring: Observations using distributed acoustic sensing and surface orbital vibrators during hydraulic fracturing

Understanding hydraulic fracturing is crucial to improving the stimulation of unconventional reservoirs and increasing fluid production. This study develops a novel seismic monitoring technology using distributed acoustic sensing (DAS) and surface orbital vibrators (SOV) to capture fracture seismic response and mechanical properties at high temporal intervals. We analyze continuous time-lapse vertical seismic profiling (VSP) data acquired every hour during the first nine days of treatment of an unconventional reservoir in the Austin Chalk/Eagle Field Laboratory. The VSP data contain clear seismic signals scattered from the activated fractures. The spatiotemporal changes of the fracture reflectivity revealed by the SOV/DAS data correlate well with the observations of fracture locations inferred from low-frequency DAS data. These results capture the fracture opening and closure processes, as well as highlight potential prestage activations of the fractures due to hydraulic connectivity with preexisting fracture systems. Therefore, analysis of the presented data set provides a unique opportunity to understand fracture initiation and subsequent evolution, not only in the context of unconventional resources but also in enhanced geothermal systems.

Correa, Julia↗

Hawai‘i Play Fairway (Final Report)

Most of Hawai'i's geothermal resources are blind—their manifestations, such as hot springs and steam vents, do not appear on the ground surface because the heated water flows far below. With the exception of $K\bar{i}lauea East Rift Zone$, in most areas of Hawai'i, high lateral permeability in the first kilometer below ground surface prevents surface thermal features from developing. As a methodology for discovering these blind resources, Play Fairway Analysis (PFA) involves finding potential locations of blind hydrothermal systems and describing potential geothermal sources in rift-zone settings. Using the PFA to find Hawai'i's geothermal resources, the University of Hawai'i (UH) conducted the Hawai‘i Play Fairway Project, Hawai'i's first statewide geothermal resource assessment since 1985. Sponsored by the U.S. Department of Energy, the Hawai'i Play Fairway Project provided an updated resource assessment, a roadmap for additional exploration activities, and the identification of areas for further exploration. Benefitting from UH's core competency in earth sciences and experienced geothermal researchers, the project comprised three phases. During the first phase, the team identified, compiled, and ranked existing geologic, groundwater, and geophysical datasets relevant to subsurface heat, fluid and permeability. Using a Bayesian statistical approach, the team developed a statistical methodology to integrate these data into a resource probability map. The team evaluated the confidence in the probability value and considered development viability of areas with geothermal resources. With these analyses, the team identified 10 locations in the Hawaiian Islands for exploration activities. For the second phase, the team collected new groundwater data in 10 locations across the state and new geophysical data on $L\bar{a}na‘i, Maui$, and central Hawai'i Island and modeled topographically induced stress to better characterize subsurface permeability. Analyzing the subsurface stresses, the team evaluated the potential for fracture-induced permeability. The team inverted the MT and gravity data to produce 3D models of resistivity and density, respectively, on $L\bar{a}na‘i$, across $Haleakal\bar{a}'s$ SW rift (Maui), and surrounding Mauna Kea (Hawai‘i Island). The team developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The team incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. Through combining data from the first and second phases, the team determined locations for further exploration during the third phase. For MT and gravity surveys, the team recommended $Kaua'i's$ $L\bar{i}hu'e$ $Basin$, the east rift of $Maui's$ $Haleakal\bar{a}$ volcano, and the southwest rift of Hawai'i Island's Mauna Loa volcano. The MT and gravity surveys aimed to enable improved confidence in the resource potential in these locations. For drilling deep groundwater well(s), the team recommended Southeast Mauna Kea and $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$. During the third phase, further exploration involved drilling a groundwater well in $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$ and performing more geophysical surveys. We deepened an existing water well proximal to our target area on $L\bar{a}na'i$ due to funding constraints that precluded us from spudding a new well that would exceed 1km depth. Drilling was preceded by a number of substantial elements including: writing an Environmental Assessment and the subsequent legal process, performance of deviation logging, lowering a camera down the well, coordinating site preparation with $P\bar{u}lama$ $L\bar{a}na'i$, shipping the UH-owned rig interisland, procuring supplies, and leading 3 community meetings on $L\bar{a}na'i$. Drilling occurred 24/7 the entire month of June 2019 over which time $L\bar{a}na'i$ $Well$ 10 was deepened from 427 m to 1057 m, with continuous core collected. We measured a roughly linear temperature gradient averaging 42°C/km and a maximum bottom hole temperature of 66°C. This gradient is more than twice the background for Hawai'i and within a range of gradients measured in this depth range for some exploration wells within KERZ. We consider these results encouraging for $L\bar{a}na'i's$ resource potential and recommend following with a slim hole within $L\bar{a}na'i's$ caldera (our target zone) to ~ 2 km. Further, the positive implications such results have for the island of O‘ahu are substantial - the shield stage of O'ahu's volcanoes ended 1-2 My earlier. However, O'ahu uses more electricity than the rest of the islands combined, and the utility recently called for 500-700MW of firm, dispatchable renewable electricity on O'ahu by 2033. In Phase 3, we also collected limited new encouraging groundwater data, and updated our thoughts on the probabilities of fluid and permeability at resource depths (Pr F = 1; Pr P = mostly unconstrained). Ultimately, we advocate for using our final probability of heat, and confidence in this probability, to drive the next phase of exploration. We contend further development of geothermal in Hawai‘i will enable the state to achieve its 100% renewable policy objective and Hawai'i to transition off of fossil fuels through geothermal discovery and development. The project not only produced a large amount of data and expanded the existing knowledge of Hawai'i's geothermal resources, but also produced publications, theses, presentations, core photos, datasets, media reports, television interviews, community events, and a blog. Students and new professionals benefitted from the project's hands-on research experiences and educational opportunities and earned awards and recognition.

15 GEOTHERMAL ENERGY↗

Geophysical borehole logging data of wells ER-GLS1, ER-GUM1, ER-PLM7, and ER-PLM8 at the East River Watershed, Colorado

The purpose of this data set is to provide a means of characterizing the bedrock variability both within a single borehole location and across the East River Watershed. The East River is part of the Watershed Function Scientific Focus Area (WFSFA) located in the Upper Colorado River Basin, United States, and all boreholes are located within Mancos shale bedrock. The data set contains geophysical borehole logging data acquired in October 2018 in boreholes ER-GUM1 (East River- Gothic Upper Montane Borehole #1) and ER-GLS1 (East River- Gothic Lower Subalpine Borehole #1), and acquired in October 2019 in boreholes ER-PLM7 (East River- Pumphouse Lower Montane Borehole #7) and ER-PLM8 (East River- Pumphouse Lower Montane Borehole #8). It contains depth, natural gamma ray emission, spectral gamma ray emission(Th, K, U), magnetic susceptibility, fluid temperature and electrical conductivity, electromagnetic (EM) induction derived rock electrical resistivity, heat pulse and spinner flow measurements to estimate groundwater flow, and Nuclear Magnetic Resonance (NMR) derived water content and hydraulic conductivity. Data for each borehole are provided as comma delimited text files. For each borehole accurate location information is provided (based on RTK GPS (Real Time Kinematic Global Positioning System) surveying). Borehole logging data were acquired using Mount Sopris borehole tools, and Vista Clara Javelin for NMR data.

54 ENVIRONMENTAL SCIENCES↗

Microseismic Event Denoising: Removal of Borehole Waves

The vertical borehole array at Farnsworth Unit, TX is used to monitor microseismic activity in the subsurface around the Carbon Capture and Sequestration (CCS) reservoir. The array consists of 16 3-component seismometers spaced vertically in a single borehole. Tube or borehole waves traveling up or down the borehole can corrupt signals of interest, such as microseismic events. A denoising convolutional neural network (DCNN) was trained to remove borehole waves from seismic waveforms of microseismic events for the purpose of reducing unwanted signal detections and better characterizing events of interest. This R&D leverages the work of Sandia colleague Rigo Tibi, who used a DCNN developed by Greg Beroza at Stanford University to improve the signal-to-noise ratio (SNR) of teleseismic events detected by the International Monitoring System.

58 GEOSCIENCES↗

Information content of spatially distributed ground-based measurements for hydrologic-parameter calibration in mixed rain-snow mountain headwaters

Parameters in hydrologic models used in mixed rain-snow regions are often uncertain to calibrate and overfitted on streamflow. To contribute addressing these challenges, we used an algorithm that assesses modeling performances through time (Dynamic Identifiability Analysis) to quantify the information content of spatially distributed ground-based measurements for identifying optimal parameter values in the Precipitation Runoff Modeling System (PRMS) model. Including spatially distributed ground-based measurements in Identifiability Analysis allowed us to unambiguously estimate more parameter values than only using streamflow (seven parameters instead of two out of a pool of thirty-three). Peaks in information gain were obtained when using dew-point temperature to identify precipitation phase-partitioning parameters. Multi-attribute identifiability analysis also yielded optimal parameter values that were temporally less variable than those estimated using streamflow alone. Overall, identifying parameter values using ground-based measurements improved the simulation of key drivers of the surface-water budget, such as air temperature and precipitation-phase partitioning. However, parameters simulating surface-to-subsurface mass fluxes like snow accumulation and melt or evapotranspiration were poorly identified by any attribute and so emerged as key sources of predictive uncertainty for this distributed-parameter hydrologic model. In conclusion, this work demonstrates the value of expanded ground-based measurements for identifying parameters in distributed-parameter hydrologic models and so diagnosing their conceptual uncertainty across the water budget.

54 ENVIRONMENTAL SCIENCES↗

In situ precipitation of hydrous ferric oxide (HFO) for remediation of subsurface iodine contamination

A practical approach for in-situ hydrous ferric oxide (HFO) precipitation was developed for iodine immobilization under field-scale conditions at the Hanford Site. A series of 1D meter-long bench top column experiments packed with Hanford sediments was conducted with a single acidic ferric solution (0.1 M, pH = 1.5) injection. Because carbonate and clay minerals are widely present in sediments, self-pH buffering of the injected acidic ferric solution occurred due to mineral dissolution, leading to HFO precipitation under a neutral condition. Up to ~12 mg/g HFO successfully precipitated and evenly distributed in the column sediments, and the remobilization of the neo-formed HFO precipitates was limited (< 4 wt.% after 100 pore volumes of flushing). The transport of IO3- in the HFO-amended sediments was strongly retarded through both adsorption and co-precipitation processes. However, reversible adsorption of iodine on HFO was observed, which might limit its application to slow moving groundwater systems.

Wang, Guohui↗

Simulations of Water Flow in Relation to a Steady Criticality in an Unsaturated Alluvial Repository

Oak Ridge National Laboratory (ORNL) is examining the conditions under which a criticality event could be ruled out for large packages placed in an unsaturated alluvial repository. This work specifically examines the water level in a stylized waste package that is assumed (1) to fail on the top, and (2) to empty as water is driven away—not by boiling—but by fluid flow phenomena driven by thermal gradients. This analysis uses the PFLOTRAN code package. Specifically, this year’s work includes examination of the heat rate at which water is driven away for various percolation rates using the same assumptions listed in the 2022 annual report. The STANDARDS (formerly UNF-ST&DARDS) package has predicted dual-purpose canister (DPC) neutron multiplication factors (k-effectives) ranging above 1.1 (unitless) with failed components. The DPC designs were not created with the intent to use for disposal purposes. Emplacing such packages may be worth consideration if it can be shown that criticality is either extremely unlikely or of very low consequence to the repository’s safety case. This work aims to help establish hydrologic boundary conditions for a highly hypothetical, highly stylized subsurface criticality event. The approach to critical is assumed to be slow, producing steady heating rates under a few kilowatts. For a DPC containing capillary media, the simulated power required for dry-out was proportional to the percolation rate but was unaffected by the size of a large breech on top of the DPC. Results for single cases with no capillary media in the DPC suggest that higher infiltration rates may be required in those cases to cause flooding. Bounding mechanisms and assumptions must be identified and characterized further by performing additional perturbations to various input parameters to demonstrate robustness of the results. This might be accomplished by examining responses to power pulses to determine the system transfer function required to enable a reduced order feedback model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

In Situ Non-Destructive Temporal Measurements of the Rhizosphere Microbiome ‘Hot-Spots’ Using Metaproteomics

Rhizosphere arguably embodies the most diverse microbial ecosystem on the planet, yet it is largely a functional ‘black box’ of belowground plant-microbiome interactions. The rhizosphere is the primary site of entry for subsurface injection of fixed carbon (C) into soil with impacts on local to global scale C biogeochemistry and ultimately Earth’s climate. While spatial organization of rhizosphere is central to its function, small scale and steep microbial and geochemical gradients within this dynamic region make it easily disrupted by sampling. The significant challenge presented by sampling blocks elucidation of discreet functions, drivers, and interactions within rhizosphere ecosystems. Here, we describe a non-destructive sampling method linked to metaproteomic analysis in order to measure temporal shifts in the microbial composition and function of rhizosphere. A robust, non-destructive method of sampling microbial hotspots within rhizosphere provides an unperturbed window into the elusive functional interactome of this system over time and space.

59 BASIC BIOLOGICAL SCIENCES↗

Ecosystems and Networks Integrated with Genes and Molecular Assemblies (ENIGMA). Final report

The goal of the Stahl lab is to use a combination of isotopic, genomic, and activity-based analyses in the field, of selected field isolates, and defined synthetic communities to identify sources and sinks of different nitrogen species in the field, focusing on the FRC at the ORNL. The most general research objective is to characterize variants in the pathway for denitrification among field isolates and synthetic microbial assemblies as a foundation for predicting the biotic and abiotic controls of nitrogen and carbon transformations in the field. Specifically, we have focused on characterizing environmental controls of nitrous oxide emissions, an important greenhouse gas and ozone depleting molecule. Studies have mapped the flux of nitrous oxide and controlling genetic elements within the FRC subsurface, indicating the importance of different enzyme variants in nitrous oxide production and consumption in controlling surface emissions. We have developed methods to quantify and catalog the isotopic signatures imparted on nitrous oxide by different nitrous oxide reductase variants in isolated organisms and model synthetic communities. We have developed methods to quantify and catalog the nitrous oxide affinity of field relevant isolates and their biomass yield when respiring nitrous oxide on different carbon substrates. We continue to expand the collection of field-relevant microbial isolates within ENIGMA by focused enrichment and isolation from field material. These objectives have been facilitated by the development of new technologies, including: 1) the development of a multiplexed automated optical density measurement system for quantifying growth kinetics of microorganisms dependent on gaseous substrates in monoculture or gas exchange in synthetic communities and 2) the development of a microcalorimetric analysis format to quantify the impact of environmental stress on microbial maintenance energy. Since field relevant processes are mediated by assemblies of organisms interacting with each other and their local environment, we have also collaborated with other ENIGMA investigators on the assembly and characterization of synthetic communities. These studies have evaluated the role of organisms lacking a complete pathway for denitrification in the control of field processes, measuring process rates of assemblies of organisms that individually lack the full denitrification pathway but mediate the complete denitrification when growing together. These efforts have so far resulted in the publication of 18 peer reviewed papers ranging from ecological and physiological drivers of pathway optimization to field-scale process characterization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multi‐Layered Systems for Permanent Geologic Storage of CO 2 at the Gigatonne Scale

Abstract The effectiveness of Carbon Capture and Storage (CCS) as an imperative decarbonization technology relies on the sealing capacity of a fine‐grained caprock to permanently store CO 2 deep underground. Uncertainties in assessing the caprock sealing capacity increase with the spatial and temporal scales and may delay CCS deployment at the gigatonne scale. We have developed a computationally efficient transport model to capture the dynamics of basin‐wide upward CO 2 migration in a multi‐layered setting over geological time scales. We find that massive capillary breakthrough and viscous flow of CO 2 , even through pervasively fractured caprocks, are unlikely to occur and compromise the storage security. Potential leakage from the injection reservoir is hampered by repetitive layering of overlying caprocks. This finding agrees with geologic intuition and should be understandable by the public, contributing to the development of climate policies around this technology with increased confidence that CO 2 will be indefinitely contained in the subsurface.

58 GEOSCIENCES↗

Investigation of mechanochemical interactions of hydrogen with earth materials in a subsurface gas storage

The goal of this project is to investigate the molecular interactions of H 2 with earth materials (EMs) that may potentially affect economics and safety of H 2 geological storage (HGS). We investigated (1) the H 2 intercalation into interlayers of phyllosilicates, (2) the competitive adsorption of H 2 /CH 4 onto porous materials, and (3) solubility of H 2 in interfacial and confined hydrocarbons. Our results indicate that (i) H 2 intercalation into hydrated interlayers is thermodynamically unfavorable and H 2 solubility in hydrated clay interlayers is in the same order of magnitude as that in bulk water, (ii) CH 4 outcompetes H 2 in adsorption onto kerogen, due to stronger CH 4 -kerogen interactions than H 2 -kerogen interactions, (iii) H 2 tends to dissolve more in oil than in water, and the introduction of CO 2 as a cushion gas reduces H 2 partitioning near the kaolinite surfaces. The outcomes provide foundational knowledge for preparing the USA for future storage site selection and storage system design, supporting DOE missions in clean and secured energy.

08 HYDROGEN↗

Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY24 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., 200 DV-1 Operable Unit), particles need to be delivered a sufficient radial distance from an injection well and retained at concentrations high enough for effective treatment. An accurate description of the particle radius of influence (ROI) is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capabilities and a lack of experimental data to validate and parameterize particle transport models. To help build toward field-scale deployment, this fiscal year (FY) we have (1) developed a pre screening tool to estimate particle transport, (2) implemented particle transport models within PFLOTRAN, and (3) conducted preliminary estimations of particle ROI. While field-scale numerical simulations will ultimately be necessary before remedy design and field implementation, we have developed a pre-screening tool that offers valuable estimations of expected particle injectability and ROI in a 1-D system. The advantage of the tool is that it does not require extensive laboratory experiments and instead makes predictions based solely on routine laboratory measurements. This tool can assist in down-selection and decision-making by identifying which particle amendment systems are worth pursuing in future laboratory experiments, such as 1-D column tests and beyond. With any system, scaling up from the lab to the field presents challenges. Currently, there is no field data available for model calibration or validation. However, the theoretical particle models being developed herein are the best tools available to guide progress toward field deployment. To help bridge this gap and verify model predictions, larger-scale lab experiments are being proposed. To advance simulation capabilities, six particle transport models are being integrated into the reactive transport simulator PFLOTRAN. These include colloid filtration theory (CFT) and five additional particle transport models (M1-M5). Each model, from M1 to M5, progressively incorporates additional particle transport and retention processes. Ultimately, the simplest model capable of accurately describing 1-D column data will be selected and parameterized. During FY24, the CFT and M1 model have been fully implemented within PFLOTRAN. Using an existing 1 D column experiment, the two currently implemented particle transport models (CFT and M1), and associated parameters, were fit to this experiment. While simpler model formulations are helpful for estimations, these formulations could not fully describe particle transport and retention behavior in the previous 1-D column experiment. Thus, additional complexities will need to be considered, which will be accounted for in the M2-M5 model formulations. Additionally, because a viscous, shear thinning fluid was required to keep particles in suspension, considerations for flow will also need to also be accounted for. Therefore, a new immiscible two-phase flow mode is currently being implemented in PFLOTRAN. With some modifications, this new flow module could also support simulation of non-Newtonian liquid amendments, foams, and emulsions. We also estimated the expected ROI of solid amendments using 1-D simulations. The average predicted ROI was approximately 15 ft for micron-sized zero valent iron (mZVI) suspended in xanthan gum (XG). Using the pre screening tool and ROI estimates, additional amendment-delivery laboratory characterization and experiments are proposed. The results from additional experiments can be used to validate and parametrize particulate transport model formulations, which will ultimately provide predictive capabilities for field amendment-delivery systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Advancing ocean monitoring and knowledge for societal benefit: the urgency to expand Argo to OneArgo by 2030

The ocean plays an essential role in regulating Earth’s climate, influencing weather conditions, providing sustenance for large populations, moderating anthropogenic climate change, encompassing massive biodiversity, and sustaining the global economy. Human activities are changing the oceans, stressing ocean health, threatening the critical services the ocean provides to society, with significant consequences for human well-being and safety, and economic prosperity. Effective and sustainable monitoring of the physical, biogeochemical state and ecosystem structure of the ocean, to enable climate adaptation, carbon management and sustainable marine resource management is urgently needed. The Argo program, a cornerstone of the Global Ocean Observing System (GOOS), has revolutionized ocean observation by providing real-time, freely accessible global temperature and salinity data of the upper 2,000m of the ocean (Core Argo) using cost-effective simple robotics. For the past 25 years, Argo data have underpinned many ocean, climate and weather forecasting services, playing a fundamental role in safeguarding goods and lives. Argo data have enabled clearer assessments of ocean warming, sea level change and underlying driving processes, as well as scientific breakthroughs while supporting public awareness and education. Building on Argo’s success, OneArgo aims to greatly expand Argo’s capabilities by 2030, expanding to full-ocean depth, collecting biogeochemical parameters, and observing the rapidly changing polar regions. Providing a synergistic subsurface and global extension to several key space-based Earth Observation missions and GOOS components, OneArgo will enable biogeochemical and ecosystem forecasting and new long-term climate predictions for which the deep ocean is a key component. Driving forward a revolution in our understanding of marine ecosystems and the poorly-measured polar and deep oceans, OneArgo will be instrumental to assess sea level change, ocean carbon fluxes, acidification and deoxygenation. Emerging OneArgo applications include new views of ocean mixing, ocean bathymetry and sediment transport, and ecosystem resilience assessment. Implementing OneArgo requires about $100 million annually, a significant increase compared to present Argo funding. OneArgo is a strategic and cost-effective investment which will provide decision-makers, in both government and industry, with the critical knowledge needed to navigate the present and future environmental challenges, and safeguard both the ocean and human wellbeing for generations to come.

ARGO↗

Machine Learning Analysis of Hydrologic Exchange Flows and Transit Time Distributions in a Large Regulated River

Hydrologic exchange between river channels and adjacent subsurface environments is a key process that influences water quality and ecosystem function in river corridors. High-resolution numerical models were often used to resolve the spatial and temporal variations of exchange flows, which are computationally expensive. In this study, we adopt Random Forest (RF) and Extreme Gradient Boosting (XGB) approaches for deriving reduced order models of hydrologic exchange flows and associated transit time distributions, with integrated field observations (e.g., bathymetry) and hydrodynamic simulation data (e.g., river velocity, depth). The setup allows an improved understanding of the influences of various physical, spatial, and temporal factors on the hydrologic exchange flows and transit times. The predictors also contain those derived using hybrid clustering, leveraging our previous work on river corridor system hydromorphic classification. The machine learning-based predictive models are developed and validated along the Columbia River Corridor, and the results show that the top parameters are the thickness of the top geological formation layer, the flow regime, river velocity, and river depth; the RF and XGB models can achieve 70% to 80% accuracy and therefore are effective alternatives to the computational demanding numerical models of exchange flows and transit time distributions. Each machine learning model with its favorable configuration and setup have been evaluated. The transferability of the models to other river reaches and larger scales, which mostly depends on data availability, is also discussed.

97 MATHEMATICS AND COMPUTING↗