Search NASA⌕ Search

SEARCH · Search NASA

Results for “Groundwater model”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The Co-Transport of PFAS and Cr(VI) in porous media

PFAS and Cr are present at some sites as co-contaminants. Here, the objective of this research was to investigate the co-transport behavior of per- and polyfluoroalkyl substances (PFAS) and hexavalent chromium (Cr(VI)) in porous media. Miscible-displacement experiments were conducted using two soils and an aquifer sediment with different geochemical properties. Perfluorooctane sulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) were employed as model PFAS. The retardation of PFOS was decreased in the presence of Cr(VI). Conversely, the transport and retardation of PFOA was not affected by the presence of Cr(VI). The reduction of PFOS retardation caused by Cr(VI) is likely due to sorption competition for both organic-carbon and inorganic (metal-oxides and clay minerals) domains. The relative contributions of the three soil constituents to PFOS sorption and the potential for competition between PFOS and Cr(VI) is a function of the geochemical composition of the porous media (i.e., organic carbon, metal-oxides and clay minerals). The PFAS had minimal impact on the retention and transport of Cr(VI). To our knowledge, the results presented herein represent the first reported data for PFOS and Cr(VI) co-transport in porous media. The results of this study indicate that the presence of Cr(VI) has the potential to increase the migration potential of PFOS in soil and groundwater, which should be considered when characterizing electroplating facilities, leather tanning facilities, and other co-contaminated sites.

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

Predicting future well performance for environmental remediation design using deep learning

Here in this study, we developed a deep learning (DL) framework with a multi-channel three-dimensional convolutional neural network (MC3D-CNN) to predict well performance and thereby assist future environmental remediation design. Such prediction of extraction well performance at designated locations is critical for configuring pump-and-treat (P&T) well network design and operation, setting reasonable target closure dates for overall remedying, and estimating remedy costs. The framework is developed with operational and monitoring data routinely collected during P&T remedy operations, including well extraction and injection rates as well as in situ contaminant concentrations. Traditionally, the collected data were rarely used for purposes other than assessing past well performance and the accuracy of the conceptual site model. However, recent advances in data-driven computational approaches enable better use of the large datasets to inform future well performance, enhance site characterization, and improve remediation planning. In this study, we established a DL framework to integrate transient three-dimensional contaminant plumes and multiple aquifer properties (e.g., hydraulic conductivity and hydrostratigraphic maps) to identify characteristic patterns controlling and representing extraction well mass recovery, aiming at providing future mass recovery estimates for existing wells and candidate wells at any proposed locations. We evaluated our framework by using a realistic synthetic dataset generated from a well-calibrated flow and transport model used in the 200 West Area of the U.S. Department of Energy’s Hanford Site in southeastern Washington state. The multi-channel feature in our framework allows integration of various types and temporal densities of training datasets for DL model development. Overall, we found that the trained DL model achieved an accuracy of over 90% in ranking extraction well performance in validation datasets, and over 80% in predicting high-performance-ranking well locations. This data-informed approach provides a flexible tool to support adaptive site management, streamline decision-making, and potentially reduce remediation time and costs. Our DL framework can be used as a filtering tool to improve the current P&T network optimization design by reducing the number of candidate well locations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Whole metagenome sequencing and 16S rRNA gene amplicon analyses reveal the complex microbiome responsible for the success of enhanced in-situ reductive dechlorination (ERD) of a tetrachloroethene-contaminated Superfund site

The North Railroad Avenue Plume (NRAP) Superfund site in New Mexico, USA exemplifies successful chlorinated solvent bioremediation. NRAP was the result of leakage from a dry-cleaning that operated for 37 years. The presence of tetrachloroethene biodegradation byproducts, organohalide respiring genera (OHRG), and reductive dehalogenase (rdh) genes detected in groundwater samples indicated that enhanced reductive dechlorination (ERD) was the remedy of choice. This was achieved through biostimulation by mixing emulsified vegetable oil into the contaminated aquifer. This report combines metagenomic techniques with site monitoring metadata to reveal new details of ERD. DNA extracts from groundwater samples collected prior to and at four, 23 and 39 months after remedy implementation were subjected to whole metagenome sequencing (WMS) and 16S rRNA gene amplicon (16S) analyses. The response of the indigenous NRAP microbiome to ERD protocols is consistent with results obtained from microcosms, dechlorinating consortia, and observations at other contaminated sites. WMS detects three times as many phyla and six times as many genera as 16S. Both techniques reveal abundance changes in Dehalococcoides and Dehalobacter that reflect organohalide form and availability. Methane was not detected before biostimulation but appeared afterwards, corresponding to an increase in methanogenic Archaea. Assembly of WMS reads produced scaffolds containing rdh genes from Dehalococcoides, Dehalobacter, Dehalogenimonas, Desulfocarbo, and Desulfobacula. Anaerobic and aerobic cometabolic organohalide degrading microbes that increase in abundance include methanogenic Archaea, methanotrophs, Dechloromonas, and Xanthobacter, some of which contain hydrolytic dehalogenase genes. Aerobic cometabolism may be supported by oxygen gradients existing in aquifer microenvironments or by microbes that produce O 2 via microbial dismutation. The NRAP model for successful ERD is consistent with the established pathway and identifies new taxa and processes that support this syntrophic process. This project explores the potential of metagenomic tools (MGT) as the next advancement in bioremediation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Impact of reservoir parameters and wellbore permeability uncertainties on CO 2 and brine leakage potential at the Shenhua CO 2 Storage Site, China

Here, this paper describes a system-level risk assessment for the Shenhua CO 2 storage site, China, using the National Risk Assessment Partnership Integrated Assessment Model for Carbon Storage (NRAP-IAM-CS). Here, we begin by determining the optimal number of Monte Carlo (MC) simulations to achieve CO 2 and brine leakage result convergence. Then, we calculate mass CO 2 and brine leakage to the atmosphere and a hypothetical shallow aquifer. Finally, we assess the geochemical impacts in the event of leakage as if there were a shallow freshwater aquifer at the Shenhua site. Simulation results show that leakage results tend to stabilize after 300 Monte Carlo simulations. When the three wells on site are assigned a permeability of 10 -11 m 2 (representing significantly leaking wells), moderate CO 2 and brine leakages occur, and the percentage of CO 2 leakage exceeds the threshold value we set based on the Intergovernmental Panel on Climate Change (IPCC). This is, however, unlikely to be the case for the Shenhua site. For all the other scenarios, the CO 2 leakage is trivial although there is still the possibility of CO 2 leakage into the groundwater aquifer exceeding the 1% threshold over 1000 years, assuming constant legacy wellbore permeability. For the significantly leaking scenario, there is a 10% probability to have a moderate (2.7 × 10 7 m 3 ) leakage-affected volume in the shallow aquifer with the pH below 6.5, and a large (1.4 × 10 8 m 3 ) volume with the total dissolved solids (TDS) above 500 ppm, hence pH and TDS may be considered for site monitoring plans. Based on the simulation results, there is a very low probability of significant CO 2 /brine leakage through the existing wells at the Shenhua CO 2 storage site.

58 GEOSCIENCES↗
DOE OSTI2023

Meteorological, reference evapotranspiration and estimated transpiration data, July 2020-Dec 2021, Slate River Floodplain, Crested Butte, Colorado

This data package includes a time series of meteorological data (air temperature, relative humidity, barometric pressure, precipitation, wind speed, soil heat flux, incoming shortwave radiation, incoming longwave radiation, outgoing shortwave radiation and outgoing longwave radiation) from the Slate River floodplain outside Crested Butte, Colorado, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? The package also includes reference evapotranspiration calculated via the Penman-Monteith equation and an estimate of willow (Salix spp.) evapotranspiration using the method of Irmak et al., J. Hydrology, 2013. Due to two malfunctioning sensors (incoming shortwave radiation from Oct 2020-May 2021 and incoming longwave radiation from May 2021-March 2022), observed net radiation and calculated evapotranspiration are only included from July-October 2020. However, users may be able to extend calculated evapotranspiration using estimates of these two parameters from other sources. Both sensors have been replaced and data beyond 2021 be uploaded to this repository when it is available.Within the data package, "FLMD.csv" describes file-level metadata and "dd.csv" defines column headers and universal terms across the dataset. The data package includes 3 "*data.csv" files, two for meteorological data for 2020 and 2021 and one for reference evapotranspiration and transpiration data from 2020 ("SLT_OBJ2_PenmanMonteith_2020_data.csv"). Each "*data.csv" file of meteorological data has a corresponding "*_InstallationMethods.csv" file that describes the location, sensor model, sensor serial number and other metadata corresponding to each measured parameter. We also include a Jupyter notebook ("PenmanMonteithET_Calculations.ipynb") that details the assumptions and equations used to calculate daily and hourly reference evapotranspiration and daily willow transpiration from raw meteorological data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Soil water content, matric potential, carbon dioxide and oxygen concentrations, Oct 2018-Dec 2021, Slate River Floodplain, Crested Butte, Colorado

This data package includes a time series of soil sensor data (temperature, water content, bulk electrical conductivity, porewater dissolved oxygen and porewater dissolved carbon dioxide) in a vertical profile from the Slate River floodplain outside Crested Butte, Colorado, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? The package includes: (1) soil temperature, volumetric water content and electrical conductivity at 40, 60 and 82.5 cm depth; (2) soil matric potential at 40, 60, 79 and 100 cm depth; (3) soil CO2 concentrations at 40, 60 and 82.5 cm depth; and (4) soil oxygen concentrations at 60, 82.5, 100, 135, 170 and 182 cm depth. Both the carbon dioxide and oxygen sensors are optical sensors that can measure the partial pressure of oxygen in both saturated and unsaturated conditions. Unfortunately, soil CO2 in the profile is unexpectedly high and above the sensor calibration range (0-25,000 ppm). In addition, soil CO2 sensors failed within a year of deployment, so we only report CO2 data from 2019-2020.Within the data package, "FLMD.csv" describes file-level metadata and "dd.csv" defines column headers and universal terms across the dataset. The data package includes 4 "*data.csv" files, one for each calendar year in the dataset. Each "*data.csv" file has a corresponding "*_InstallationMethods.csv" file that describes the location, sensor model, sensor serial number and other metadata corresponding for each measured parameter. Because sensors have been added over time, not every sensor has data dating back to Oct 2018. Note that there is a data gap over winter 2019-2020 due to a power outage. While this repository currently only contains data through December 2021, the dataset will be updated as additional years are collected and processed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

3D Seismic Anatomy of a Watershed Reveals Climate–Topography Coupling That Drives Water Flowpaths and Bedrock Weathering

To investigate how bedrock transforms to soil, we mapped the topography of the interface demarcating onset of weathering under an east-west trending shale watershed in the Valley and Ridge province in the U.S.A. Using wave-equation travel-time tomography from a seismic array of >4000 geophones, we obtained a 3D P-wave velocity (Vp) model that resolves structures ~20 m below land surface (mbls). The depth of mobile soil and the onset of dissolution of chlorite roughly match Vp = 600 m/s and Vp = 2700 m/s, respectively. Chlorite dissolution initiates porosity growth in the shale matrix. Depth to the 2700 m/s contour is greater under the N- as compared to S-facing hillslopes and under sub-planar as compared to concave-up land surfaces. Broadly, the geometries of the ‘soil’ and ‘chlorite’ Vp contours are consistent with the calculated potential for shear fracture opening under weak regional compression. However, this calculated fracture potential does not consistently explain observations related to N- vs. S-facing aspect nor fracture density observed by borehole televiewer. Apparently, regional compression is only a secondary influence on Vp: the primary driver of P-wave slowing in the upper layers of this catchment is topographic control of reactive water flowpaths and their integrated effects on weathering. Here, the Vp result is best explained as the long-term integrated effect of groundwater flow-induced geochemical weathering of shale in response to climate-driven patterns of micro- and macro-topography.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

An Assessment of Potential Dose Impacts from External Contamination on Naval Reactors Facility Waste Canisters (Special Analysis)

This Special Analysis (SA) was performed to address a request by Naval Reactors Facility (NRF) Waste Programs for a permanent exception to limits of removable surface contamination on the exterior of waste canisters shipped to the Remote-Handled Low-Level Waste (RHLLW) Disposal Facility as specified in the waste acceptance criteria (WAC) (PLN-5446). The purpose of this SA is to determine if the NRF-requested levels for removable surface contamination on the exterior of all NRF waste canisters is within the bounds of the current performance assessment (PA). This was done by calculating the groundwater all-pathways dose contribution from surface contamination on the exterior of NRF canisters for the following cases: (1) the exteriors of all NRF waste canisters are contaminated to the 10 CFR 835 Appendix D allowable limits in the current WAC, and (2) the exteriors of all NRF waste canisters are contaminated to the limits requested by NRF Waste Programs. Dose impacts for each case were compared to each other and to the all pathways dose for the PA base case. Dose impacts were also compared to the all pathways dose limit specified in DOE O 435.1. A simple assessment of the potential impacts of the increase in surface contamination using the NRF-requested limits on the biotic, air, and inadvertent intruder pathways was also performed. Based on the results of this SA, the increase in canister exterior contamination limits requested by NRF are well within the bounds of the current PA and will not result in a violation of performance objectives. The results also show the increased limits do not reflect or necessitate a fundamental change to the PA conceptual model, nor a change to the way exterior contamination is not included in PA dose calculations. Therefore, it is recommended the NRF request for an exception to the current limits for external surface contamination be accepted and the revised limits for NRF-generated waste canisters be added to the WAC. The monitoring plan will also be revised to identify external canister contamination as a potential mobile source term that may be detected by monitoring earlier than potential releases from waste. The PA, composite analysis (CA), closure plan, and PA/CA maintenance plan do not require revision. Recommendations are also included for inclusion of NRF procedures used to limit water pool radionuclide variability (and thus canister surface contamination variability) to the waste generator certification process.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Quantifying local and global mass balance errors in physics-informed neural networks

Physics-informed neural networks (PINN) have recently become attractive for solving partial differential equations (PDEs) that describe physics laws. By including PDE-based loss functions, physics laws such as mass balance are enforced softly in PINN. This paper investigates how mass balance constraints are satisfied when PINN is used to solve the resulting PDEs. We investigate PINN’s ability to solve the 1D saturated groundwater flow equations (diffusion equations) for homogeneous and heterogeneous media and evaluate the local and global mass balance errors. We compare the obtained PINN’s solution and associated mass balance errors against a two-point finite volume numerical method and the corresponding analytical solution. We also evaluate the accuracy of PINN in solving the 1D saturated groundwater flow equation with and without incorporating hydraulic heads as training data. We demonstrate that PINN’s local and global mass balance errors are significant compared to the finite volume approach. Tuning the PINN’s hyperparameters, such as the number of collocation points, training data, hidden layers, nodes, epochs, and learning rate, did not improve the solution accuracy or the mass balance errors compared to the finite volume solution. Mass balance errors could considerably challenge the utility of PINN in applications where ensuring compliance with physical and mathematical properties is crucial.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Iodine Air Sampling: Installation and Preliminary Results

The F-Area wetlands are a known sink for radioiodine generated at the Savannah River Site (SRS). Wastewater containing iodine-129 ( 129 I) was discharged into the F-Area Seepage Basins for disposal, but the acidity of the wastewater mobilized some of the contaminants from the basin soil through the vadose zone and into the groundwater; ultimately creating a groundwater plume that extends towards Fourmile Branch and its associated wetland areas. Dissolved radioiodine specifically has migrated to the nearby F-Area wetlands, where it persists as a contaminant of concern. Accumulation and resurfacing of 129 I in the wetlands and nearby areas creates a potential secondary source zone where changes in biogeochemical conditions could result in the release of this sequestered 129 I to the air. Monitoring data has indicated that 129 I has a seasonal behavior at several surface water stations, with high concentrations during summertime and low concentrations during wintertime. Still, it is not known which factors control the seasonal behavior of 129 I. Further, it has been reported that 129 I contributed approximately 23% to the 2019 SRS offsite air pathway dose to the representative person for that year (SRNS, 2020), but air measurements of 129 I have never been taken in the F-Area wetlands. Because air emissions of 129 I from the wetlands have never been investigated, it is not known if this could be a potential release pathway that needs to be included in the site conceptual model. To address these unknowns, the Savannah River National Laboratory (SRNL) has an applied research task to investigate whether 129 I in the wetlands is emitted to the air and, if so, whether air concentrations exhibit seasonal variation. This report presents information about the system recently installed in the F-Area wetlands for measuring 129 I, as well as analysis results for samples collected June through August 2023. Quantification of this constituent will provide a better understanding of the contribution of airborne 129 I from the F-Area wetlands to the offsite dose and better inform ongoing and future remedial actions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Response of Subsurface Nitrogen-Cycling Microbial Communities to Environmental Fluctuations (Final Technical Report)

Riparian floodplains are dynamic ecosystems linking terrestrial and riverine systems. These floodplains experience hydrological shifts such as changes in water table height, flooding, and drought and can be ‘hotspots’ of biogeochemical cycling due to shifting sediment moisture (and saturation) and subsurface exchanges of water, nutrients, and other compounds across different sediment layers. Subsurface microbial communities are the primary drivers of biogeochemical processes in floodplains, and thus their structure and function can directly influence both surface and groundwater quality. The microbial nitrogen (N) cycle is particularly important in floodplains as it affects nutrient availability and removal. Two functional guilds of chemoautotrophic (i.e. CO2-fixing) microorganisms are responsible for the first oxidative step of the N cycle, nitrification: ammonia-oxidizing archaea (AOA) and bacteria (AOB) catalyze the oxidation of ammonia to nitrite, while nitrite-oxidizing bacteria (NOB) oxidize nitrite to nitrate. Despite the critical role nitrification plays in N-cycling in both terrestrial and aquatic ecosystems, our understanding of the diversity, ecophysiology, and activity of nitrifying organisms in subsurface floodplain soils/sediments is extremely limited. To help address this critical knowledge gap, the overarching goal of this project was to determine how shifts in key environmental parameters and gradients impact microbial N-cycling communities/processes, with particular emphasis on nitrification, within hydrologically-variable floodplain sediments in the Wind River Basin near Riverton, Wyoming. The three specific objectives of this project were to: (1) to associate in situ environmental drivers of N cycling with distinct functional guilds; (2) determine the guild response to variation in key ecosystem drivers; and (3) develop a dynamic ecosystem model of the microbial N cycle with the Riverton subsurface using community genomic and biogeochemical data collected in the first two objectives. Over the course of this project, we employed both 16S rRNA gene amplicon sequencing and genome-resolved metagenomics to examine the phylogenetic diversity and metabolic potential of subsurface nitrifier communities within 68 samples collected across multiple sites, depths, and time points within the Riverton floodplain, allowing for both spatial and temporal investigations at different scales. This project benefitted tremendously from recent advances in high-throughput sequencing technologies coupled with dramatic improvements in the computational tools and algorithms available for analyzing such large, complex genomic datasets. By pairing these cutting-edge genomic approaches with depth-resolved sampling and detailed geochemical analyses of the Riverton floodplain, we have gained novel insights into the structure and function of subsurface nitrifier communities in relation to both hydrology and biogeochemistry. This project resulted in the most detailed and comprehensive characterization of N-cycling floodplain microbial communities to date and will hopefully inspire and pave the way for future studies using similar approaches in other floodplains. Indeed, such information is critical for understanding subsurface biogeochemical cycling and how elemental stores are altered from perturbations initiated by the water cycle within floodplains. Finally, because of the terrestrial-aquatic nature of the Riverton floodplain, results from this project are also of relevance to disciplines such as soil science, estuarine science, limnology & oceanography, biogeochemistry, geobiology, environmental engineering, as well as genomics and data science.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Energy‐Water Asynchrony Principally Determines Water Available for Runoff From Snowmelt in Continental Montane Forests

ABSTRACT Changes in the volume, rate, and timing of the snowmelt water pulse have profound implications for seasonal soil moisture, evapotranspiration (ET), groundwater recharge, and downstream water availability, especially in the context of climate change. Here, we present an empirical analysis of water available for runoff using five eddy covariance towers located in continental montane forests across a regional gradient of snow depth, precipitation seasonality, and aridity. We specifically investigated how energy‐water asynchrony (i.e., snowmelt timing relative to atmospheric demand), surface water input intensity (rain and snowmelt), and observed winter ET (winter AET) impact multiple water balance metrics that determine water available for runoff (WAfR). Overall, we found that WAfR had the strongest relationship with energy‐water asynchrony (adjusted r 2 = 0.52) and that winter AET was correlated to total water year evapotranspiration but not to other water balance metrics. Stepwise regression analysis demonstrated that none of the tested mechanisms were strongly related to the Budyko‐type runoff anomaly (highest adjusted r 2 = 0.21). We, therefore, conclude that WAfR from continental montane forests is most sensitive to the degree of energy‐water asynchrony that occurs. The results of this empirical study identify the physical mechanisms driving variability of WAfR in continental montane forests and are thus broadly relevant to the hydrologic management and modelling communities.

Webb, Ryan William↗
DOE OSTI2026

BSEC ecohydrological and water quality fluxes from RHESSys Simulations in USGS gauged watersheds

Baltimore Environmental Social Collaborative (BSEC) Water and Water Quality Simulations from RHESSys Model The repository contains RHESSys (Tague & Band, 2004; source code) simulated ecohydrological and nutrient (nitrogen only) fluxes at daily, basin-average (RHESSys_basin_output) and monthly, grid (RHESSys_patch_output) levels. We currently simulated the following 8 watersheds in Baltimore: Dead Run Baisman Run Scotts Level Branch Moores Run Powder Mill Run Maidens Choice Run Stony Run The watershed boundaries of all studied watersheds are stored in Watershed_Boundary folder. Variables and their units are listed in the metadata. Spatial projection, NAD83 / UTM zone 18N (EPSG:26918) is used for patch-level, netCDF-format files. For more information, please contact Ruoyu Zhang (rz3jr@virginia.edu).

Baltimore MD↗
DOE OSTI2023

Application of unsupervised deep learning to image segmentation and in-situ contact angle measurements in a CO 2 -water-rock system

Rock surface wettability is a critical property that regulates multiphase flows in porous media, which can be quantified using the surface contact angle (CA). X-ray micro-computed tomography (μCT) provides an effective approach to in-situ measurements of surface CAs. However, the CA measurement accuracy depends significantly on the quality of CT image segmentation, which is the clustering of CT pixels into separate phases. Inspired by this, we developed a deep learning (DL)-based CA measurement workflow. Motivated by the recent tremendous progress in unsupervised learning techniques and aiming to avoid expensive manual data annotations, an unsupervised DL pipeline for CT image segmentation was proposed and implemented, which includes unsupervised model training and post-processing. The unsupervised model training was driven by a novel loss function constrained with feature similarity and spatial continuity and implemented by iterative forward and backward paths; the former clustered the pixel-wise feature vectors extracted by convolution neural networks, whereas the latter updated the parameters using gradient descent. An over-segmentation strategy was adopted for model training. The post-processing steps based on agglomerative hierarchical clustering (AHC) were implemented to further merge the over-segmented model output to the desired cluster number, which is intended to improve the efficiency of image segmentation. The developed unsupervised DL pipeline was compared with other commonly-used image segmentation methods using pixel-wise and physics-based evaluation metrics on a synthetic raw-image dataset, which had a known ground truth. The unsupervised DL pipeline showed the best performance. Next, the segmented images were input to an automatic CA measurement tool, and the results were validated by comparisons with manual measurements. The CA values from the manual and automatic measurements showed similar distributions and statistical properties. The automatic measurement demonstrated a wider spectrum because of the much larger number of measurement data points. The primary novelty of the unsupervised DL pipeline developed in this study lies in the novel loss function and the over-segmentation strategy associated with AHC post-processing. Finally, the workflow has been proven an efficient tool for pore-scale wettability characterization, which has a wide range of applications in fundamental studies of multiphase flows in natural porous media, which have critical implications to geological carbon sequestration, hydrocarbon energy recovery, and contaminant transport in groundwater.

42 ENGINEERING↗
DOE OSTI2020

Analysis of Hydrology and Interim Measure Performance in the Chromium Plume at Los Alamos National Laboratory - 20498

Hexavalent chromium, Cr(VI), is present in the regional aquifer beneath Los Alamos National Laboratory (LANL) at concentrations greater than the New Mexico groundwater standard (50 ppb). The source of the Cr is blowdown from a power plant where potassium dichromate was used as a corrosion inhibitor from 1956 to 1972. Blowdown comingled with effluent from another source and migrated approximately 3 miles down a canyon, primarily as surface flow, until reaching an alluvial infiltration zone. The infiltrating water moved vertically through a thick (approximately 900 ft), complex stratigraphic sequence of unsaturated tuffs, basalt, and unconsolidated alluvial deposits. Horizontal transport also took place within spatially limited perched groundwater horizons in the vadose zone (VZ). Breakthrough at the regional aquifer water table is thought to have occurred in multiple hydraulic windows or 'drip points' located approximately 3 miles from the initial release site at the power plant. The regional aquifer is highly heterogeneous and anisotropic, and displays unconfined behavior in an upper zone and confined behavior at depth. The Cr(VI) contamination has so far remained confined to the upper approximately 25 m of the aquifer. The ambient hydraulic gradient in the regional aquifer is toward the east/southeast, but the water table is extremely flat in the Cr plume area compared to surrounding areas. There is some local variability in flow direction, potentially related to effects of seasonal pumping from nearby Los Alamos county water-supply wells, which draw from the regional aquifer several hundred feet below the water table. An interim measure (IM) is currently operating at the site to achieve hydraulic control of plume migration using a series of extraction wells, ion-exchange treatment, and injection along the downgradient portion of the plume. The migration of Cr contamination in the aquifer has been previously evaluated using a calibrated probabilistic numerical model, which has also been used to design the spatial configuration of extraction and injection wells for the IM, and to identify locations for monitoring wells. The model is built using the Finite Element Heat and Mass (FEHM) transfer code, and the Model Analysis and Decision Support (MADS) software is used for analyses. To date, the portion of the IM system that has operated in a continuous manner has been very effective at hydraulic control. Modifications to the flow field associated with the influence of IM pumping and injection are evaluated here, with modeling constraints to be provided by recent observations of tracers injected into both monitoring wells and IM injection wells, as well as opportunistic observations of geochemical signatures of treated ion-exchange effluent injected into IM injection wells. The tracer and geochemical data provide insights into modified flow directions within the Cr plume area that cannot be otherwise obtained. Results of the modeling work are being used to inform adaptive management of IM operations and will be a valuable tool for development and refinement of long-term remediation strategies. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

The CO2 mineralization potential of New Mexico basalts and associated critical element behavior and mobility.

This presentation reports on laboratory experiments assessing the carbon dioxide (CO₂) mineralization potential of New Mexico basalts, including trachybasalt and basalt lithologies. Using simulated groundwater and controlled conditions, the study evaluates the kinetics of CO₂ sequestration, mineral saturation indices, and secondary mineral formation. It also investigates the mobility of critical and trace elements such as Ni, Zn, Cu, and Li during basalt dissolution. Results indicate that New Mexico basalts can effectively sequester CO₂ through carbonate formation, with grain size and mineral composition influencing reaction rates and elemental behavior. These findings support the viability of geologic carbon storage in mafic volcanic terrains and inform strategies for critical element recovery.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Pore-resolved investigation of turbulent open channel flow over a randomly packed permeable sediment bed

Pore-resolved direct numerical simulations are performed to investigate the interactions between streamflow turbulence and groundwater flow through a randomly packed porous sediment bed for three permeability Reynolds numbers, Re K = 2.56 , 5.17 and 8.94, representative of natural stream or river systems. Time–space averaging is used to quantify the Reynolds stress, form-induced stress, mean flow and shear penetration depths, and mixing length at the sediment–water interface (SWI). Here, the mean flow and shear penetration depths increase with Re K and are found to be nonlinear functions of non-dimensional permeability. The peaks and significant values of the Reynolds stresses, form-induced stresses, and pressure variations are shown to occur in the top layer of the bed, which is also confirmed by conducting simulations of just the top layer as roughness elements over an impermeable wall. The probability distribution functions (p.d.f.s) of normalized local bed stress are found to collapse for all Reynolds numbers, and their root-mean-square fluctuations are assumed to follow logarithmic correlations. The fluctuations in local bed stress and resultant drag and lift forces on sediment grains are mainly a result of the top layer; their p.d.f.s are symmetric with heavy tails, and can be well represented by a non-Gaussian model fit. The bed stress statistics and the pressure data at the SWI potentially can be used in providing better boundary conditions in modelling of incipient motion and reach-scale transport in the hyporheic zone.

turbulence simulation↗
DOE OSTI2021

Numerical and Analytical Modeling of Flow Partitioning in Partially Saturated Fracture Networks

Infiltration processes in fractured-porous media remain a crucial, yet not very well understood component of recharge and vulnerability assessment. Under partially-saturated conditions flows in fractures, percolating fracture networks, and fault zones contribute to the fastest spectrum of infiltration velocities via preferential pathways. Specifically, the partitioning dynamics at fracture intersections determine the magnitude of flow fragmentation into vertical and horizontal components and hence the bulk flow velocity and dispersion of fracture networks. In this work, we derive an approximate analytical solution for the partitioning process and validate it using smoothed particle hydrodynamics simulations. The developed transfer function is conceptually based on simulation results and laboratory experiments carried out in previous works. It allows to efficiently simulate flow through fracture networks with simple cubic structure and arbitrary number of fractures and aperture sizes via linear response theory and convolution of a given input signal. We derive a non-dimensional bulk flow velocity ($\widetilde{v}$) and dispersion coefficient ($\widetilde{D}$) to characterize fracture networks in terms of dimensionless horizontal and vertical time scales $\tau_m$ and $\tau_0$. The dispersion coefficient is shown to strongly depend on the horizontal time scale and converges towards a constant value of $0.08$ within reasonable fluid and geometrical parameter ranges, while the non-dimensional velocity exhibits a characteristic $\widetilde{v} \sim \tau_m^{-1/2}$ scaling. Given that hydraulic information is often only available at limited places within (fractured-porous) aquifer system, such as boreholes or springs, our study intends to provide a rudimentary analytical concept to potentially reconstruct internal fracture network geometries from external boundary information, e.g., the dispersive properties of discharge (groundwater level fluctuations).

42 ENGINEERING↗