Search NASA⌕ Search

SEARCH · Search NASA

Results for “Contaminant transport”

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

Integration of Large-Scale Electrical Imaging into Geological Framework Development and Refinement

Geologic framework models (GFMs) are critical to the construction of reliable simulation models of groundwater flow and contaminant transport. To support GFM development, direct information (e.g., core samples, fluid samples, hydraulic testing) tends to be sparse and separated by large distances relative to the spatial scales of aquifer heterogeneity. There are additional challenges associated with highly contaminated legacy waste sites, where drilling is particularly costly, and invasive sampling requires specialized handling and disposal of hazardous materials. At these sites in particular, non-invasive geophysical imaging can play an important role in filling spatial gaps between boreholes and reducing characterization costs by optimizing and minimizing the number of necessary boreholes. Here this paper presents a case study demonstrating the use of large-scale (> 30 km 2 ) electrical mapping to identify hydrostratigraphy and potential paleochannels at the Hanford Site, located in Washington State, USA. In two field campaigns, over 36 line-kilometers of electrical resistivity tomography (ERT) data were collected along 14 transects. ERT surveys were sited and performed to image critical aspects (e.g., paleochannels, stratigraphic contacts) of the subsurface, demonstrating a general workflow for integrating ERT with GFM development. Inconsistencies between the GFM and ERT were catalogued to provide a basis for future site characterization using complementary geophysical methods and (or) direct sampling.

58 GEOSCIENCES↗
DOE OSTI2021

Vadose and Saturated Zone Flow and Transport Model Package Report for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington

This model package report (MPR) documents the development of the integrated vadose and saturated zone flow and transport model developed for the performance assessment of Trenches 31 and 34 in the low-level burial grounds (LLBGs). This modeling capability is intended for use in addressing the analysis requirements outlined in DOE O 435.1, Chg 1, Radioactive Waste Management1. The overall objective of the modeling effort is to provide a basis for making informed disposal decisions pertinent to Trenches 31 and 34. The purpose of the MPR is to document the development of the three-dimensional numerical vadose zone and saturated flow and transport model test case and evaluate its adequacy to support the LLBG performance assessment. The purpose is not to present results for DOE 435.1 decision making. The use of the model to perform base case and sensitivity analysis for the LLBG performance assessment, including inputs and results, is documented separately in subsequent environmental calculation files. This report discusses the development and translation of the conceptual model for flow and contaminant transport into the LLBG performance assessment three-dimensional numerical flow and transport model evaluated using the Subsurface Transport Over Multiple Phases (STOMP©) simulator. The development of representative geologic framework is described along with the implementation of waste release models used to represent contaminant release from waste disposed in the trenches. The report also provides the technical basis for specific model parameters and boundary conditions, along with description of modeling assumptions. This MPR includes certain calculations that are necessary to demonstrate the soundness of the model. Results provided by the model include vadose zone and saturated zone flow fields, and estimates of the possible future concentration in groundwater of technetium-99 and iodine-129 as example test cases. As an evaluation of the test cases, the model estimates of current vadose conditions are compared to available and analogous field and laboratory data, and compared to the results from the original performance assessment model (WHC-EP-0645, Performance Assessment for the Disposal of Low-Level Waste in the 200 West Area Burial Grounds 2 ). The features, events, and processes applicable to vadose zone and saturated zone flow and transport model are identified, and representative initial estimates for various parameters are documented. Note that the parameter estimates presented in this MPR are for illustration purposes, and may or may not reflect values selected for eventual performance assessment. Several key topical discussions (i.e., basis for recharge estimates, basis for vadose zone modeling, basis for saturated zone model development, and calibration) are included, which serve as the groundwork for confidence building for the groundwater pathway modeling and results. Numerical simulation results based on an example test case are included to illustrate the use of the combined saturated-unsaturated model in performance assessment calculations. Sensitivity and uncertainty results are not included in this MPR. Those results will be included in future environmental calculation files. The inclusion of the trapezoidal trench geometry and construction details in the finite difference grid introduces some gross simplifications regarding the trench and liner systems. The model test case presented in this MPR includes the assumption that the trench liner system does not affect flow through or around the trenches after the liner system is assumed to fail. This MPR also includes results of other test cases that involve alternate assumptions about how to incorporate the hydraulic effects of the trenches into the vadose zone of the model. These alternate cases do not attempt to account for the presence of the liner system after its assumed failure either. None of these cases is considered to be the base case at this time. Analysis and alternate cases that attempt to account for the presence of the liner system in greater detail, and its effect on flow through and around it, are to be documented in subsequent environmental calculation files. Depending on the results of those analyses and alternate cases, the eventual base case may involve more detailed inclusion of the effects of the liner system hydraulics.

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

Sensor Recommendations for Long Term Monitoring of the F-Area Seepage Basins

In mid-2018, a new paradigm for long-term monitoring was developed after of decade of applied research projects funded by the Department of Energy’s office of Environmental Management Technology Development program. The program at SRNL was focused on transitioning complex environmental waste sites from active to passive remediations strategies. A key result of these studies was that the use of enhanced attenuation approaches at radiologically contaminated sites will result in the creation of secondary source areas in the subsurface that will require monitoring for decades. Alternative monitoring approaches are being developed and tested at the Savannah River Site’s F-Area Hazardous Waste Management Facility, the new paradigm provides innovative solutions that will significantly lower costs of monitoring through the coupling of data collection, machine learning and deterministic groundwater modeling. The foundation of this approach is a well-optimized network of sensors for measuring hydrogeochemical master variables that control, and therefore act as indicators of groundwater contaminant transport. By monitoring changes in the controlling master variables over time arising from geological and environmental shifts, predictive modelling can assist with identifying new strategies for ensuring regulatory requirements are met if trends toward conditions for potential remobilization of attenuated contaminants are detected. In this report, we evaluated commercially available single parameter sensor platforms (e.g., temperature/depth) and configurable multi-parameter sensor platforms (e.g., pH, oxidation-reduction potential, temperature, depth, dissolved oxygen, and conductivity). Each was scored using an optimization function based on how well the system supports the proposed long-term monitoring paradigm, in general, and the site-specific conditions at F-Area, in particular. Several viable sensor systems were identified. Of these, a combined platform including the In-Situ Aqua TROLL 500 multi-parameter sensor platform and the In-Situ temperature/depth sensor had the highest rating and was identified as the most suitable candidate for installation and monitoring of the master variables and potentiometric surface that control groundwater contaminant plumes emanating from the F-Area Seepage Basins. The discussion of recommended potential deployment locations builds upon recommendations made by Denham et al (2019).

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

Adaptive Multirate Mass Transfer (aMMT) Model: A New Approach to Upscale Regional‐Scale Transport Under Transient Flow Conditions

Abstract The long‐term evaluation of regional‐scale groundwater quality needs efficient upscaling methods for transient flow. Upscaling techniques, such as the Multirate Mass Transfer (MRMT) method with constant upscaling parameters, have been used for transport with steady‐state flow, yet the upscaling parameters (i.e., rate coefficients) may be time dependent. This study proposed and validated an adaptive MRMT (aMMT) method by allowing the mass transfer coefficients in MRMT to change with the flow field. Advective‐dispersive contaminant transport simulated in a 3‐D heterogeneous medium was used as a reference solution. Equivalent transport under homogeneous flow conditions was evaluated by applying the MRMT and aMMT models for upscaling. The relationship between mass transfer coefficients and flow rates was fitted under steady‐state flow driven by various hydraulic gradients. A power law relationship was obtained, which was then used to update the mass transfer coefficients in each stress period under transient flow conditions in the aMMT method. Results indicated that for advection‐dominated transport, both the MRMT and aMMT methods can upscale the anomalous transport dynamics affected by subgrid heterogeneity under transient flow conditions. Whereas for diffusion‐dominated systems, the MRMT model failed to capture the tails of tracer breakthrough curves after the boundary condition changed, but the results from the aMMT model were significantly improved. However, if the overall flow direction changed, both MRMT and aMMT failed to represent the breakthrough curve tail generated by the heterogeneous system. The results point toward a promising path for upscaling transport in complex aquifers with transient flow.

Guo, Zhilin↗
DOE OSTI2020

The Influence of Aqueous Dispersivity on Transport of Radiological Material at LANL MDA G - 20450

Aqueous dispersivity is a phenomenon by which a solute moving in groundwater mixes with water that does not contain the solute and becomes diluted. Solutes take paths of varying lengths through pore networks or experience slower velocities due to interactions with pore walls. Aqueous dispersivity leads to breakthrough curves with earlier arrival than expected in purely advective systems and persistent 'tails,' or late-arriving, dilute concentrations of the solute. Dispersivity can be longitudinal (in the direction of flow) or transverse (orthogonal to flow). LANL has disposed of radioactive material at Material Disposal Area (MDA) G since the late 1950's. Most of the disposed materials at MDA G are low-level waste (LLW), but other types of radioactive wastes are disposed or stored at the site. For example, there is a collection of remote-handled transuranic waste (RH-TRU) waste in 33 shafts, which were intended to be stored retrievably. Aqueous-phase radionuclides included in models of contaminant transport at MDA G are subject to dispersive processes in porous media of the vadose zone and saturated zone. To model dispersivity at MDA G probabilistically, we created site-specific distributions of aqueous dispersivity. First, we performed a comprehensive review of peer-reviewed publications, reports from other performance assessments of MDA G, and documentation from sites that were considered similar based on expert opinion, such as Yucca Mountain. Data were included from a wide range of sources, so that distributions developed based on the data reflect the current state of knowledge in properties of the subsurface environment. Data were grouped based on value type (e.g., determined via experiment or used in a different model), rock type, zone type (saturated or unsaturated), and site, and were weighted in the distributions according to their relevance to the site and zone type. We combined Monte Carlo simulation and bootstrapping techniques to produce approximate distributions of mean aqueous dispersivity. Distributions of dispersivity for saturated and unsaturated zone transport are differentiated by data weight and the length scales over which transport occurs. We find that underestimating or disregarding transverse dispersivity may lead to overstating simulated well capture in the saturated zone at MDA G. Neglecting transverse dispersivity is a common assumption due in part to use of one-dimensional transport models that do not consider solute spreading orthogonal to flow. Also, we find that variability in longitudinal dispersivity is associated with substantial travel time alterations in the MDA G unsaturated zone and in simulated contaminant breakthrough times at a downstream observation plane in the MDA G saturated zone. Finally, at large scales (>1000 m), a linear trend between the dispersivity coefficient and length over which transport occurs may overestimate the dispersivity coefficient. Future modeling work with application to solute spreading could include interacting advective and diffusive zones in the highly fractured basalt below MDA G. (authors)

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

Contaminant Source Identification from Finite Sensor Data: Perron–Frobenius Operator and Bayesian Inference

Sensors in the built environment ensure safety and comfort by tracking contaminants in the occupied space. In the event of contaminant release, it is important to use the limited sensor data to rapidly and accurately identify the release location of the contaminant. Identification of the release location will enable subsequent remediation as well as evacuation decision-making. In previous work, we used an operator theoretic approach—based on the Perron–Frobenius (PF) operator—to estimate the contaminant concentration distribution in the domain given a finite amount of streaming sensor data. In the current work, the approach is extended to identify the most probable contaminant release location. The release location identification is framed as a Bayesian inference problem. The Bayesian inference approach requires considering multiple release location scenarios, which is done efficiently using the discrete PF operator. The discrete PF operator provides a fast, effective and accurate model for contaminant transport modeling. The utility of our PF-based Bayesian inference methodology is illustrated using single-point release scenarios in both two and three-dimensional cases. The method provides a fast, accurate, and efficient framework for real-time identification of contaminant source location.

42 ENGINEERING↗
DOE OSTI2020

Smart Ventilation for Advanced California Homes

This project investigated smart ventilation approaches to minimize energy use for providing indoor air quality (IAQ) in high performance new California homes. Evaluation criteria included annual ventilation-related energy, peak energy and time-of-use savings, and the indoor air quality relative to a minimally code-compliant ventilation system. The simulations used CONTAM’s air flow and contaminant transport model, combined with the EnergyPlus building loads model. House types representing the default California Energy Code compliance homes were investigated for four California climate zones, covering a wide range of climate types. Both single and multi-zone smart ventilation controls were investigated. Contaminant sources included contaminants emitted continuously and varying with time, temperature and relative humidity, episodic emissions from occupant activities and outdoor particles. Single-zone ventilation controls that varied ventilation depending on outdoor temperatures were able to consistently save half of ventilation-related energy without compromising long-term IAQ. Ventilation strategies that tracked occupancy were less successful, because this work included generic contaminants with constant background emission rates. Energy performance for occupancy controls improved with a one-hour pre-occupancy flush out strategy. The addition of zoning ventilation controls did not offer significant IAQ to energy improvements compared to non-zonal versions of the same ventilation system type. The best controls had HVAC energy savings of 10-20%, with individual cases reaching up to 40% savings. However, these savings cannot be achieved without worsening personal exposures for at least one contaminant. A metric is needed to assess the competing changes in exposure to different contaminants in order to determine the net-health impacts of a control strategy. Controls that directly sensed contaminants and controlled them to acceptable levels showed that the California OEHHA limit for formaldehyde completely dominates system performance, with homes not able to meet the limit even with continuous operation of a fan sized to twice the current code minimum.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Assessing Effects of Climate Change on Legacy Waste at the Enewetak Atoll

The Republic of the Marshall Islands (RMI) is in the central Pacific Ocean ~4,500 km west of Hawaii. The Enewetak Atoll, located in the northwest part of the RMI, was the site for 43 nuclear weapon tests between 1948 and 1958. Fallout and deposition from the tests contaminated the island surfaces, lagoon waters and sediment, and nearby ocean waters at the atoll. In the 1970s, a cleanup effort collected radioactive waste and placed it in the Cactus Crater on Runit Island (also called the Runit Dome). In December 2021, Congress directed the U.S. Department of Energy to study the impacts of climate change on the Runit Dome nuclear waste disposal site. Pacific Northwest National Laboratory (PNNL) assembled a multidisciplinary team of climate scientists, ocean modelers, environmental scientists, and health physicists to assess the likely effects of remaining radionuclides at the Enewetak Atoll. PNNL’s approach focused on effects of tropical cyclones that were postulated to mobilize and transport contaminated lagoon sediments and result in human and biota exposure. PNNL’s study estimated (1) the radionuclide source term, (2) the effects of climate change on severe storms, (3) mobilization and transport of radionuclides, and (4) radiation dose to humans and biota. Radionuclides in the lagoon and/or ocean waters of the Enewetak Atoll were characterized by the U.S. Atomic Energy Commission (AEC) in 1972, Woods Hole Oceanographic Institution in 2015, and Lawrence Livermore National Laboratory in 2018. The RMI Nationwide Radiological Study was conducted in the early 1990s for radionuclides remaining in island soils. The 1972 AEC survey remains the most comprehensive source of radionuclide data on lagoon sediments. Climate change modeling at a regional scale in the central Pacific Ocean is limited. PNNL climate scientists simulated severe historical storms postulated to occur both in a recent climate (2015) and in the future (2090) using the Advanced Research Weather Research and Forecasting (WRF-ARW) model, employing a pseudo-global-warming technique. A postulated complete, future failure of the Runit Dome was also considered. PNNL developed a high-resolution regional ocean hydrodynamics model covering the entire RMI extended economic zone using the Finite Volume Coastal Ocean Model (FVCOM). The FVCOM model was run using global reanalysis data for current climate and WRF-ARW simulation for the future climate. PNNL also developed a radionuclide fate and transport model using the FVCOM Integrated Compartment Model (FVCOM-ICM) to simulate the current and future mobilization and transport of radionuclides sorbed to lagoon sediments and the exchange of radionuclides between the water and sediment. FVCOM-ICM-predicted radionuclide concentrations were then used to estimate radiation dose to humans and biota at all islands of the Enewetak Atoll. Under current climate conditions, annual radiation exposures for the southern islands including Enewetak (Fred) and Medren (Elmer) were below the current U.S. standards. Radiation doses were somewhat elevated starting at Runit Island northward and westward to Enjebi Island (Janet). The islands in the northwest quadrant, particularly Bokoluo (Alice) and Bokombako (Belle), remain relatively contaminated. The islands in the southwestern quadrant have low contamination. The highest contribution to radiation doses comes from consumption of locally grown foods. Two radionuclides, 90Sr and 137Cs, contributed the greatest fraction for most terrestrial foods. In current climate conditions, the storms temporarily increased radionuclide concentrations in the lagoon waters, increasing the radiation dose slightly. In future conditions, doses are expected to be smaller, primarily because of the radioactive decay of the shorter-lived radioisotopes of 90Sr and 137Cs. This could make all islands in the far northwest of the atoll – except Bokombako (Belle) and perhaps Bokoluo (Alice) – suitable for residency. For the f

Prasad, Rajiv↗
DOE OSTI2023

Single-Well Push–Pull Tracer Test Analyses to Determine Aquifer Reactive Transport Parameters at a Former Uranium Mill Site (Grand Junction, Colorado)

At a former uranium mill site where tailings have been removed, prior work has determined several potential ongoing secondary uranium sources. These include locations with uranium sorbed to organic carbon, uranium in the unsaturated zone, and uranium associated with the presence of gypsum. To better understand uranium mobility controls at the site, four single-well push–pull tests (with a drift phase) were completed with the goal of deriving aquifer flow and contaminant transport parameters for inclusion in a future sitewide reactive transport model. This goes beyond the traditional use of a constant sorption distribution coefficient (K d ) and allows for the evaluation of alternative remedial injection fluids, which can produce variable K d values. Dispersion was first removed from the resulting data to determine possible reactions before conducting reactive transport simulations. These initial analyses indicated the potential need to include cation exchange, uranium sorption, and gypsum dissolution. A reactive transport model using multiple layers to account for partially penetrating wells was completed using the PHT-USG reactive transport modeling code and calibrated using PEST. The model results quantify the hydraulic conductivity and dispersion parameters using the injected tracer concentrations. Uranium sorption, cation exchange, and gypsum dissolution parameters were quantified by comparing the simulated versus observed geochemistry. All simulations required some cation exchange and calcite equilibrium, and one simulation required gypsum dissolution to improve the model fit for calcium and sulfate. Uranium sorption parameters were not strongly influenced by the other parameter values but were highly influenced by uranium concentrations during the drift phase, with possible kinetic rate limitations. Thus, a future recommendation for such push–pull tests is to collect more geochemical data during the drift phase. The final uranium sorption parameters were within the range of values determined from prior column testing. The flow and transport parameters derived from these single-well push–pull tests will provide initial parameters for any future sitewide reactive transport model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Predicting multiphase flow and tracer transport for an underground chemical explosive test

Detecting radionuclide gas seepage from clandestine underground nuclear tests is central to nonproliferation explosion monitoring research. Yet, early-time (<6 day) gas transport driven by the explosive pressure wave remains poorly constrained due to scarcity of field data. We simulate multi-phase gas transport in the vadose zone using pre-shot data from a recent chemical explosion in P-Tunnel at the Nevada National Security Site, USA. Despite using a simplified 2D-radial model, predictions of tracer arrival matched observations within one order-of-magnitude. Our results show how transient blast forcing rapidly mobilizes gases from the cavity into surrounding rock – critical for optimizing sensor placement and test planning. This unique integration of field data and modeling represents a significant improvement in our ability to predict gas migration from underground explosions. More broadly, it offers insights into the coupled dynamics of pressure waves and contaminant transport in the vadose zone, with implications for monitoring and hazard assessment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evaluation of the Radioactive Material Released in the Harborview Research and Training Building and Some Implications for Emergency Response

On 2 May 2019, during the 137 Cs source recovery operation, a source capsule in a research irradiator containing approximately 77.1 TBq was breached. Based on a geometric reconstruction analysis of the damage to the capsule, approximately 46.3 GBq (0.04%) was impacted by the chop saw (grinder) inside a mobile hot cell on the loading dock at the University of Washington Harborview Research and Training (HRT) Building. A very small fraction of the material impacted, less than 1%, was released from the mobile hot cell and then to the rest of the HRT Building. The objectives of this project were to assess the accidental release of 137 CsCl and its implications related to emergency response methods and the ramifications of 137 CsCl transport. The phenomenology of this event was also compared with past alkali halide dispersal events. The vast number of measurements and samples collected by the remediation contractors, the Department of Energy’s Nuclear Emergency Support Team, and the small number of retrospective samples collected by the authors informed the analysis. The techniques included (1) autoradiography and electron microscopy of samples collected from the HRT Building and the irradiator, (2) 3D visualization of deposition on surfaces and within the ventilation system, and (3) a study of the damage to the source capsule to evaluate the Cs particle size and particle composition due to the grinding accident. Subsequently, the cesium contaminant transport through the numerous pathways in the building was reconstructed to assess the deposition on surfaces as a function of particle size. Furthermore, the implications for emergency response are relevant to data quality and management. A Data Quality Objective guides data collection methods so that they have appropriate accuracy and precision for the intended application. Recommendations were made with respect to the sample collection protocols and archiving of samples.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2024

Evaluation of the Radioactive Material Release in the Harborview Research and Training Building and Implications for Emergency Response

On May 2, 2019, during the 137 Cs source recovery operation, a source capsule in a research irradiator containing approximately 77.1 TBq was breached. Based on a geometric reconstruction analysis of the damage to the capsule, approximately 46.3 GBq (0.04%) was impacted by the chop saw (grinder) inside a Mobile Hot Cell (MHC) on the loading dock at the University of Washington Harborview Research and Training (HRT) Building. A very small fraction of the material impacted, less than 1%) was released from the Mobile Hot Cell and then to the rest of the HRT Building. The objectives of this project were to assess the accidental release of 137 CsCl and its implications related to emergency response methods and the ramifications of 137 CsCl transport. The phenomenology of this event was also compared with past alkali halide dispersal events. The vast number of measurements and samples collected by the remediation contractors, the Department of Energy's Nuclear Emergency Support Team, and the small number of retrospective samples collected by the authors informed the analysis. The techniques included (1) autoradiography and electron microscopy of samples collected from the HRT Building and the irradiator, (2) 3D visualization of deposition on surfaces and within the ventilation system, and (3) a study of the damage to the source capsule to evaluate the Cs particle size and particle composition due to the grinding accident. Subsequently, the cesium contaminant transport through the numerous pathways in the building was reconstructed to assess the deposition on surfaces as a function of particle size. The implications for emergency response are relevant to data quality and management. A Data Quality Objective (DQO) guides data collection methods so that they have appropriate accuracy and precision for the intended application. Recommendations were made with respect to the sample collection protocols and sample archival.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2024

Wildland–urban interface wildfire increases metal contributions to stormwater runoff in Paradise, California

The 2018 Camp Fire was a large late-year (November) wildfire that produced an urban firestorm in the Town of Paradise, California, USA, and destroyed more than 18 000 structures. Runoff from burned wildland areas is known to contain ash, which can transport contaminants including metals into nearby watersheds. However, due to historically infrequent occurrences, the effect of wildland-urban interface (WUI) fires, such as the Camp Fire, on surface water quality has not been well-characterized. Therefore, this study investigated the effects of widespread urban burning on surface water quality in major watersheds of the Camp Fire area. Between November 2018 and May 2019, 140 surface water samples were collected, including baseflow and stormflow, from burned and unburned watersheds with varying extent of urban development. Samples were analyzed for total and filter-passing metals, dissolved organic carbon, major anions, and total suspended solids. Ash and debris from the Camp Fire contributed metals to downstream watersheds via runoff throughout the storm season. Increases in concentration up to 200-fold were found for metals Cr, Cu, Ni, Pb, and Zn in burned watersheds compared to pre-fire values. Total concentrations of Al, Cd, Cu, Pb, and Zn exceeded EPA aquatic habitat acute criteria by up to 16-fold for up to five months after the fire. To assess possible transport mechanisms and bioavailability, a subset of 18 samples was analyzed using four filters with nominal pore sizes ranging from 0.22 to 1.2 μm to determine the particulate size distribution of metals. Trace and major metals (Al, Ba, Co, Cr, Cu, Fe, Hg, Mn, Ni, Pb, and Zn) were found mostly associated with larger grain sizes (>0.45 μm), and some metals (Al, Cr, Fe, and Pb) also included a substantial colloidal phase (0.22 to 0.45 μm). This study suggests that fires in the wildland–urban interface increase metal concentrations, mainly through particulate driven transport. The metals with the largest increases are likely from anthropogenic disaster materials, though biomass ash also is a major contributor to water quality. In conclusion, the increase in metals following WUI burning may have adverse ecological impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Review: Pre-Darcy flows in low-permeability porous media

The widely used Darcy’s law specifies a linear relation between the Darcy velocity of fluid flow and the pressure gradient that drives the flow. However, studies have shown that Darcy velocity can exhibit a nonlinear dependence on the pressure gradient in low-permeability porous media such as clay and shale when the pressure gradient is adequately low. This phenomenon is referred to as low-velocity non-Darcian flow or pre-Darcy flow. This paper provides a comprehensive review of the theories, experimental data, and modeling methods for pre-Darcy flow in low-permeability porous media. The review begins by outlining the fundamental mechanisms underlying pre-Darcy flow that regulate the unique characteristics such as nonlinear dependence of the Darcy velocity on the pressure gradient and its relevance to fluid–rock interactions. The review then proceeds to present a thorough compilation of experimental investigations performed in various low-permeability geomaterials including tight sandstones, shales, and clays. Next, empirical and theoretical models and simulation methods that have been developed to fit and interpret experimental data are reviewed. Finally, the review underscores the challenges encountered in conducting and interpreting pre-Darcy flow experiments and suggests future research directions. By analyzing previous experimental investigations, this review aims to offer a valuable resource for researchers and practitioners seeking to enhance their understanding of fluid dynamics in low-permeability geomaterials. This provides insights into the application of pre-Darcy flow in numerous natural and engineered processes such as shale oil and gas recovery, contaminant transport in low-permeability aquifers, and geological disposal of nuclear waste.

58 GEOSCIENCES↗
DOE OSTI2020

Benzene homologues contaminants in a former herbicide factory site: distribution, attenuation, risk, and remediation implication

Benzene homologues often used as organic raw materials or as detergents in chemical industry are prone to accidental release into the environment which can cause serious long-term soil pollutions. In a large former herbicide factory site, we investigated 43 locations for benzene homologues contaminations in soil, soil gas, and groundwater; and studied the hydrogeological conditions. An inverse distance weighted interpolation method was employed to determine the pollutants three-dimensional spatial distribution in the soils. Results showed that benzene homologues residues were mainly originated from the herbicide production workshop; and that the pollution had horizontally expanded at the deeper soil layer. Contaminants had already migrated 15m downward from ground surface. Contaminant phase distribution study showed that NAPL was the primary phase (>99%) for the pollutants accumulated in the unsaturated zone; while it had not migrated to groundwater. The primary mechanism for contaminant transport and attenuation included dissolution of “occluded” NAPL into pore water and pollutant volatilization into soil pore space. Risk assessment revealed that the pollutants brought unacceptable high carcinogenic and non-carcinogenic risks to public health. In order to convert this former chemical processing factory site into a residential area, a remediation to the polluted production workshop sites is urgently required.

Yang, Shuo↗
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

Assessment of peak power demand reduction available via modulation of building ventilation systems.

Peak power demand strains electrical grids and increases cost of electricity generation, transmission and distribution infrastructure. Many studies have examined ways of reducing peak power demand, including modification of room air temperature setpoints or reduction of lighting levels. However, very few or no studies have examined the peak power reduction resource offered by temporary curtailment of building ventilation systems. For this reason, we conducted a simulation campaign in which we examined the resource offered by temporary ventilation curtailment in commercial buildings of different use types across the United States and in residences in the state of California, with the essential constraint that any changes resulted in air quality acceptable to occupants through additional ventilation to compensate for the curtailment. To do this, we employed previously validated building models implemented in the airflow and contaminant transport tool CONTAM and building thermal and systems modeling tool EnergyPlus, in some cases co-simulated. Overall, results show savings are highly dependent on building type and climate but range from 0–2 W/ft 2 and up to 40% of total peak building power demand. Depending on building type, this power shed can be conducted for 1.5–8 h before acute exposure or odor concerns are expected, assuming a safety factor of 2. This reduction is of the same magnitude as that offered by thermal control strategies such as setpoint increase, or from lighting reduction strategies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Optimizing Cryo-Focused Pyrolysis GC/MS for Tracing Soil Organic Matter Across Diverse Ecosystems

The cycling of organic matter in terrestrial soils and sediments is central to a range of biogeochemical processes that regulate nutrient cycling, crop productivity, trace gas emissions, and contaminant transport. Pyrolysis-gas chromatography/mass spectrometry (py-GC/MS) is a powerful tool for characterizing bulk soil organic matter (SOM) at the molecular level. In this study, we used a cryo-focused py-GC/MS system to analyze soil samples from seven diverse ecosystems: vernal pool, prairie pothole, temperate forest, tropical forest, tundra, wildfire-affected boreal forest, and grassland. We addressed a key bottleneck in molecular-level SOM characterization by developing an automated data analysis pipeline to optimize py-GC/MS and complementary evolved gas analysis/mass spectrometry (EGA/MS) methods, incorporating advanced tools for peak deconvolution, developing a custom compound class library, and implementing fragmentation spectrum-based molecular networking for the first time. This improved workflow was applied to soil samples from all seven ecosystems, including multiple depths and density fractions. Our findings demonstrate that ecosystem type plays a dominant role in shaping compositional differences in SOM. We also identified trends in the source of SOM compounds (e.g., microbial vs plantderived) across soil depth and density fractions, which are critical for understanding persistence and turnover of SOM. Our molecular networking analysis indicated that although many compounds are widespread across ecosystems, others are restricted to specific environments, such as wetlands. This underscores the utility of molecular-level data in elucidating the complexity of SOM composition and the environmental drivers that shape it. Such molecular-level insights can deepen our knowledge of biogeochemical SOM cycles.

54 ENVIRONMENTAL SCIENCES↗