SEARCH · Search NASA

Results for “subsurface flow”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2012

3-D Inversion of MT Data From the Raft River Geothermal Field - Preliminary Results

Geothermal systems are, in many cases, characterized by geoelectric structures with discernible anomalous zones, related to the presence of fluids and conductive alteration minerals. The MT method has been successfully used since the 1970’s to image subsurface electrical resistivity, from which in turn can be delineated structural controls and system characteristics such as the clay cap and fluid distributions. Toward this purpose, a 3-D magnetotelluric (MT) survey of 105 stations was carried out at the Raft River Geothermal area in southern Idaho, to further characterize the field generally, fault distributions, and paths of or barriers to fluid flow. This work was done in support of developing the EGS potential of the field, as part of a Department of Energy-funded project. An electrical resistivity model of the area was constructed from the MT data using a 3-D inversion program on a multicore workstation. The model is interpreted using available lithological logs and geochemical data. Discernible in the resistivity structure are separate low-resistivity zones, attributed to increased fluids and clay alteration, within the basement from which the geothermal field produces, and within the overlying sediments, proximal to the wells. Here, the inferred compartmentalization of the field is similar in nature to that inferred from fluid geochemistry analysis. A much deeper conductor with possible lower crustal sources is seen also which appears to correlate with a similar feature in the regional Earthscope MT project.

15 GEOTHERMAL ENERGY↗

DOE OSTI2019

Path Forward for Groundwater Compliance Strategy at the Tuba City, Arizona, Disposal Site - 19581

The Tuba City, Arizona, Disposal Site is located in the Navajo Nation and near the Hopi Reservation. A uranium mill operated at the site from 1956 to 1966. During the operational period, mill tailings were slurried to an unlined impoundment. Excess slurry water was transferred to unlined evaporation ponds. Groundwater contamination resulting from the infiltration of fluids in the mill tailings impoundment and evaporation ponds was characterized in the 1990's. The site's Groundwater Compliance Action Plan was prepared in 1999, and a groundwater extraction and treatment system was constructed in 2002. Issues with aquifer productivity and treatment plant reliability led to the suspension of operations in 2014. The limitations of the remediation system prompted a realization that the site conceptual model should be reevaluated. An integrated study to better understand the geochemical and hydrologic factors affecting contaminant transport, to support the development of a new groundwater compliance strategy, was initiated in 2017 and is ongoing. In addition to the groundwater science components, the study included the evaluation of enhanced evaporation methods. The integrated evaluation will lead to refinement of the site groundwater flow model and development of a contaminant transport model and subsequent development of remedial alternatives. Remedial alternatives will be evaluated to determine a new groundwater compliance strategy. The objectives of geochemical testing were to confirm and quantify geochemical reactions controlling uranium fate and transport. These reactions are important because they control the distribution of uranium on the solid phase; potential for release from the solid phase; mobility of uranium in the water phase; and potential for uranium redistribution onto the currently uncontaminated solid phase. A key outcome is the change in understanding of the mass of uranium in the subsurface that is potentially soluble. These results contradict the previous basis for remediation: that uranium was almost exclusively present in soluble form, within the boundaries of the contaminant plume. Testing showed that uranium sorption downgradient from the plume was measurable but relatively minimal. The movement of other mill-related constituents (nitrate, sulfate, alkalinity) results in the dissolution of naturally occurring uranium and a higher than background concentration of uranium in groundwater. These constituents are not subject to sorption, and their downgradient concentrations are reduced only by diffusion and dispersion. The geochemical testing provides conceptual and quantitative information on potential plume migration. While the testing revealed the potential availability of uranium contaminant sources, the aquifer buffering capacity, pH, and presence of oxygen in groundwater may, individually or in combination, be controlling factors for plume movement. Evaluation of evaporative treatment alternatives was undertaken to determine one or more sustainable and efficient methods capable of treating 151 liters per minute (L/min) (40 gallons per minute), the annual average flow rate. In addition to establishing the baseline capacity for evaporation using the site's existing pond, prospective methods of enhanced evaporation included preheating water before inflow to the evaporation pond, vertical circulation in the pond, dye addition, and Wind-Aided Intensified Evaporation. Consideration of a 151 L/min evaporative treatment arose from previous analysis of extraction well production rates and contaminant mass removal, which indicated that about 80% contaminant removal, relative to full-scale operation of the site's groundwater treatment plant, could be achieved with a smaller capacity system. The lower-flow evaporative option also represented a potentially more sustainable treatment system in comparison to restarting or reconfiguring the groundwater treatment plant. Aquifer pumping tests were designed to develop or improve estimates of hydraulic properties such as transmissivity, horizontal and vertical hydraulic conductivity, storativity, and storage coefficient. Testing involved pumping from one well at a constant flow rate for up to 5 days while monitoring groundwater drawdown in the pumping well and five surrounding wells and then shutting down the pumping well to observe groundwater recovery. The aquifer pumping test results provide insight regarding the hydraulics of the saturated zone. The saturated zone is in a low-conductivity formation with preferential flow in the horizontal direction. The aquifer exhibits low vertical conductance. The use of AQTESOLV software indicates that the saturated zone is a layered system due to the cross-bedded sand-grain orientation within the Navajo Sandstone. These indicators of limited vertical flow provide a potential justification for the observed depth of the contaminant plume of 15-46 meters (50-150 feet) below ground surface). Given the relatively low vertical conductance, groundwater extraction should be focused on the upper portions of the saturated zone. Aquifer response testing provided an improved understanding of aquifer characteristics. However, the primary objective, to identify optimal groundwater extraction strategies, has not been fully addressed. A list of path forward activities for hydrologic analyses was developed, with an overarching objective of providing input parameters for groundwater flow and contaminant transport modeling. Defensible predictions of contaminant transport, through modeling, are critical to the development of a new site groundwater compliance strategy. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2014

Microbial Activity and Precipitation at Solution-Solution Mixing Zones in Porous Media -- Subsurface Biogeochemical Research

The goal for this research was to understand how best to add compounds to receptive microbial communities in porous media in order to achieve optimal calcite precipitation in a volumetrically significant space and to understand the physiological health of the cells that are responsible for the calcite precipitation. The specific objectives were to: (1) develop better tools for visually examining biofilms in porous media and calcium carbonate precipitation being mediated by microbes in porous media, and (2) demonstrate the effectiveness of using that tool within a flow cell model system.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

ParFlow Sand Tank: A tool for groundwater exploration

The ParFlow Sand Tank model is an open source application designed to allow users to interactively simulate and visualize groundwater movement through the subsurface. The app is designed for both research and education; teaching hydrogeology concepts and making it easy explore and run sophisticated groundwater simulations. Our goal is to support increased accessibility and usability of research grade hydrology tools for research and teaching. The Sand Tank application simulates groundwater and surface water fluxes as well as contaminant transport in real time using the integrated physical hydrology model ParFlow (Kollet & Maxwell, 2006; Maxwell & Miller, 2005; Osei-Kuffuor et al., 2014) and the particle tracking code EcoSlim (Maxwell et al., 2019). ParFlow is a numerical hydrology model that simulates spatially distributed groundwater and surface water flow. It is a well established research tool with more than 90 publications documenting its development use to advance our understanding of groundwater dynamics and groundwater surface water interactions from the hillslope to the continental scale e.g. (Condon et al., 2020; Condon & Maxwell, 2019; Maxwell & Condon, 2016). It is designed for efficient parallel computation and has been run on many platforms spanning from laptops to supercomputers. However, one of the challenges of ParFlow is that it requires significant training and hydrologic expertise to develop simulations. The Sand Tank application makes this model accessible to anyone for education and exploration. Our application uses ParFlow for its simulation backend and ParaView for the data loading and processing. The communication infrastructure relies on the ParaViewWeb framework. We use model templates deployed in Docker images to setup the Sand Tank framework. Users can build the application locally or interact with it through our web deployment. When interacting with a template users can interactively change model parameters like subsurface processes or pump/inject water into the subsurface and watch the system respond to their changes in real time as the simulation runs. Additionally, our template setup will allow more advanced users to build custom templates of increasing complexity for both research and educational purposes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Splitting tensile strength of shale cores: intact versus fractured and sealed with ureolysis-induced calcium carbonate precipitation (UICP)

Ureolysis-induced calcium carbonate precipitation (UICP) is a biomineral solution where the urease enzyme converts urea and calcium into calcium carbonate. The resulting biomineral can bridge gaps in fractured shale, reduce undesired fluid flow, limit fracture propagation, better store carbon dioxide, and potentially enhance well efficiency. The mechanical properties of shale cores were investigated using a modified Brazilian indirect tensile strength test. An investigation of intact shale using Eagle Ford and Wolfcamp cores was conducted at varying temperatures. Results show no significant difference between shale types (average tensile strength = 6.19 MPa). Eagle Ford displayed higher strength at elevated temperature, but temperature did not influence Wolfcamp. Comparatively, cores with a single, lengthwise heterogeneous fracture were sealed with UICP and further tested for tensile strength. UICP was delivered via a flow-through method which injected 20–30 sequential patterns of ureolytic microorganisms and UICP-promoting fluids into the fracture until permeability reduced by three orders of magnitude or with an immersion method which placed cores treated with guar gum and UICP-promoting fluids into a batch reactor, demonstrating that guar gum is a suitable inclusion and may reduce the number of flow-through injections required. Tensile results for both delivery methods were variable (0.15–8 MPa), and in some cores the biomineralized fracture split apart, possibly due to insufficient sealing and/or heterogeneity in the composite UICP-shale cores. Notably in other cores the biomineralized fracture remained intact, demonstrating more cohesion than the surrounding shale, indicating that UICP may produce a strong seal for subsurface application.

58 GEOSCIENCES↗

DOE OSTI2014

Transport of Organic Contaminants Mobilized from Coal through Sandstone Overlying a Geological Carbon Sequestration Reservoir

Column experiments were conducted using a wetted sandstone rock installed in a tri-axial core holder to study the flow and transport of organic compounds mobilized by scCO2 under simulated geologic carbon storage (GCS) conditions. The sandstone rock was collected from a formation overlying a deep saline reservoir at a GCS demonstration site. Rock core effluent pressures were set at 0, 500, or 1000 psig and the core temperature was set at 20 or 50°C to simulate the transport to different subsurface depths. The concentrations of the organic compounds in the column effluent and their distribution within the sandstone core were monitored. Results indicate that the mobility though the core sample was much higher for BTEX compounds than for naphthalene. Retention of organic compounds from the vapor phase to the core appeared to be primarily controlled by partitioning from the vapor phase to the aqueous phase. Adsorption to the surfaces of the wetted sandstone was also significant for naphthalene. Reduced temperature and elevated pressure resulted in greater partitioning of the mobilized organic contaminants into the water phase.

Zhong, Lirong↗

DOE OSTI2021

Geostatistical Mapping of Salinity Conditioned on Borehole Logs, Montebello Oil Field, California

We present a geostatistics-based stochastic salinity estimation framework for the Montebello Oil Field that capitalizes on available total dissolved solids (TDS) data from groundwater samples as well as electrical resistivity (ER) data from borehole logging. Data from TDS samples (n = 4924) was coded into an indicator framework based on falling below four selected thresholds (500, 1000, 3000, and 10,000 mg/L). Collocated TDS-ER data from the surrounding groundwater basin were then employed to produce a kernel density estimator to establish conditional probabilities for ER data (n = 8 boreholes) falling below the selected TDS thresholds within the Montebello Oil Field area. Directional variograms were estimated from these indicator coded data, and 500 TDS realizations from conditional indicator simulation were generated for the subsurface region above the Montebello Oil Field reservoir. Simulations were summarized as 3D maps of median TDS, most likely salinity class, and probability for exceeding each of the specified TDS thresholds. Results suggested TDS was below 500 mg/L in most of the study area, with a trend toward higher values (500 to 1000 mg/L) to the southwest; consistent with the average regional groundwater flow direction. Discrete localized zones of TDS greater than 1000 mg/L were observed, with one of these zones in the greater than 10,000 mg/L range; however, these areas were not prevalent. The probabilistic approach used here is adaptable and is readily modified to include additional data and types and can be employed in time-lapse salinity modeling through Bayesian updating.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Distributed Temperature Sensing (DTS) to Monitor CO 2 Migration in an Enhanced Oil Recovery Field in Northern Michigan

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. MRCSP conducted analysis of CO 2 flow patterns using distributed temperature sensing (DTS) data acquired from the Chester 16 pinnacle reef located in Otsego County, Michigan. The reef has one injection well Chester 6-16 and one monitoring well Chester 8-16 where DTS is installed. DTS provides a means to track temperature within the wellbores in real time and characterize the migration patterns of CO 2 in the subsurface when DTS data is combined with other operational data such as the injection rates, bottomhole pressures and temperature. The primary objective of this study was to determine the distribution of CO 2 inflow zones in the A1 Carbonate and Brown Niagaran formations.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2016

Low Temperature Geothermal Play Fairway Analysis For The Appalachian Basin: Phase 1 Revised Report November 18, 2016

Geothermal energy is an attractive sustainable energy source. Yet project developers need confirmation of the resource base to warrant their time and financial resources. The Geothermal Play Fairway Analysis of the Appalachian Basin evaluated risk metrics that communicate the favorability of potential low-temperature geothermal energy resources in reservoirs more than 1000 m below the surface. This analysis is focused on the direct use of the heat, rather than on electricity production. Four risk factors of concern for direct-use geothermal plays in the Appalachian Basin portions of New York, Pennsylvania, and West Virginia are examined individually, and then in combination: 1) thermal resource quality, 2) natural reservoir quality, 3) induced seismicity, and 4) utilization opportunities. Uncertainty in the risk estimation is quantified. Based on these metrics, geothermal plays in the Appalachian Basin were identified as potentially viable for a variety of direct-use-heat applications. The methodologies developed in this project may be applied in other sedimentary basins as a foundation for low temperature (50-150 °C), direct use geothermal resource, risk, and uncertainty assessment. Three methods with which to combine the four risk factors were used. Among these, the averaging of the individual risk factors indicates the most favorable counties within the study area are the West Virginia counties of Monongalia, Harrison, Lewis (dubbed the Morgantown–Clarksburg play fairway), Putnam, and Kanawha (Charleston play fairway), the New York counties of Chemung and Steuben plus adjacent Bradford county in Pennsylvania (Corning–Ithaca play fairway), and the Pennsylvania counties of Mercer, Crawford, Erie, and Warren, and adjacent Chautauqua county in New York (together, the Meadville–Jamestown play fairway). These higher priority regions are surrounded by broader medium priority zones. Also worthy of additional exploration is a broad region near Pittsburgh Pennsylvania, for which the available geological data are insufficient to fully analyze the geological risks but yet the population is high. First, to assess the spatial variation in the depth to which one would need to drill to obtain geothermal temperatures that are useful to a future project, the project used bottom-hole temperature data from Appalachian Basin oil and gas exploration. These bottom hole temperature data are abundant but of low quality. Second, the project examined the potential for sufficient water flow rates through rocks to harvest heat from a geothermal well field, considering only natural reservoirs. This analysis provides a very incomplete picture of spatial variability of natural reservoirs because the oil and gas reservoir data lack key properties and are spatially biased toward those locations with profitable amounts of hydrocarbons in the rock pore spaces. Third, in light of the fact that earthquake activity has been induced in several states by subsurface work related to the oil and gas industry, this project examined the potential for similar activity in the Appalachian Basin. Acknowledging that data for such a task are insufficient, we utilized what was available: records of seismic activity, regional estimates of the orientations of stress in the rocks, and locations and orientations of zones of lateral change in rock properties at depths down to several kilometers below Earth’s surface. With these data, we created a first approximation of spatially variable risks for induced earthquakes. Because no data existed with which to test the reliability of these methods, the results have a high degree of uncertainty. Fourth, we examined the spatial variability of the above-the-ground factors that contribute to the economical viability of projects to tap low-temperature geothermal resources for direct-use. We worked principally with population density as a regionally known variable that would impact the cost of district heating. The resulting maps omit the costs of producing the hot water from the ground, because the below-ground costs are directly coupled to the thermal resource risk factor and natural reservoir risk factor – later analyses of those costs will be needed. The result of the district heating analysis is highly skewed: few census locations yielded a low estimated surface cost. The team also identified more than 165 prospects for high value direct-use geothermal energy opportunities such as industrial sites, university campuses, and federal facilities, among others. At the closure of this regional analysis, the most significant technical uncertainties are 1) reservoir distribution and capacities; 2) validity of thermal resource maps, and 3) the holistic estimation of Levelized Cost of Heat for favorable geological situations.

15 GEOTHERMAL ENERGY↗

DOE OSTI2026

Data for Soil Oxygen Dynamics: A Key Mediator of Tile Drainage Impacts on Coupled Hydrological, Biogeochemical, and Crop Systems

Tile drainage removes excess water and is an essential, widely adopted management practice to enhance crop productivity in the US Midwest and throughout the world. Tile drainage has been shown to significantly change hydrological and biogeochemical cycles by lowering the water table and reducing the residence time of soil water, although examining the complex interactions and feedbacks in an integrated hydrology–biogeochemistry–crop system remains elusive. Oxygen dynamics are critical to unraveling these interactions and have been ignored or oversimplified in existing models. Understanding these impacts is essential, particularly so because tile drainage has been highlighted as an adaptation under projected wetter springs and drier summers in the changing climate in the US Midwest. We used the ecosys model that uniquely incorporates first-principle soil oxygen dynamics and crop oxygen uptake mechanisms to quantify the impacts of tile drainage on hydrological and biogeochemical cycles and crop growth in corn–soybean rotation fields. The model was validated with data from a multi-treatment, multi-year experiment in Washington, IA. The relative root mean square error (rRMSE) for the corn and soybean yield in validation is 5.66 % and 12.57 %, respectively. The Pearson coefficient (r) of the monthly tile flow during the growing season is 0.78. Plant oxygen stress turns out as an emergent property of the equilibrium between the soil oxygen supply and biological demand. The impact of tile drainage on the system is achieved through a series of coupled feedback mechanisms. The model results show that tile drainage reduces the soil water content and enhances soil oxygenation. It additionally increases the subsurface discharge and elevates inorganic nitrogen leaching, with seasonal variations influenced by climate and crop phenology. The improved aerobic condition alleviates crop oxygen stress during wet springs, thereby promoting crop root growth during the early growth stage. The development of greater root density, in turn, mitigates water stress during dry summers, leading to an overall increase in the crop yield by ∼6 %. These functions indicate the potential of tile drainage in bolstering crop resilience to climate change and the use of this modeling tool for large-scale assessments of tile drainage. The model reveals the underlying causal mechanisms that drive the agroecosystem response to drainage on the coupled hydrology, biogeochemistry, and crop system dynamics.

Modeling↗

DOE OSTI2014

Scale-Dependent Fracture-Matrix Interactions and Their Impact on Radionuclide Transport: Development of efficient particle-tracking methods

Matrix Diffusion and Adsorption within a rock matrix are important mechanisms for retarding transport of radionuclides in fractured rock. Due to computational limitations and difficulties in characterizing complex subsurface systems, diffusive exchange between a fracture network and surrounding rock matrix is often modeled using simplified conceptual representations. There is significant uncertainty in “effective” parameters used in these models, such as the “effective matrix diffusivity”. Often, these parameters are estimated by fitting sparse breakthrough data, and estimated values fall outside meaningful ranges, because simplified interpretive models do not consider complex three-dimensional flow. There is limited understanding of the relationship between the effective parameters and rock mass characteristics including network structure and matrix properties. There is also evidence for an apparent scale-dependence in “effective matrix diffusion” coefficients. These observations raise questions on whether fracture-matrix interaction parameters estimated from small-scale tracer tests can be used for predicting radionuclide fate and transport at the scale of DOE field sites. High-resolution three-dimensional Discrete-Fracture-Network-Matrix (DFNM) models based on well-defined local scale transport equations can help to address some of these questions. Due to tremendous advances in computational technology over the last 10 years, DFNM modeling in relatively large domains is now feasible. The overarching objective of our research is to use DFNM modeling to improve fundamental understanding of how effective parameters in conceptual models are related to fracture network structure and matrix properties. An advanced three-dimensional DFNM model is being developed, which combines upscaled particle-tracking algorithms for fracture-matrix interaction and a parallel fracture-network flow simulator. The particle-tracking algorithms allow complexity in flow fields at different scales, and track transport across fracture-matrix interfaces based on rigorous local approximations to the transport equations. This modeling approach can incorporate aperture variability, multi-scale preferential flow and matrix heterogeneity. We developed efficient particle-tracking methods for handling matrix diffusion and adsorption on fracture walls and demonstrated their efficiency for use within the context of large-scale complex fracture network models with variability in apertures across a network of fractures and within individual fractures.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Soil oxygen dynamics: a key mediator of tile drainage impacts on coupled hydrological, biogeochemical, and crop systems

60 APPLIED LIFE SCIENCES↗

DOE OSTI2013

Self-Assembling Sup-porosity: The Effect On Fluid Flow And Seismic Wave Propagation

Fractures and joints in the field often contain debris within the void spaces. Debris originates from many different mechanisms: organic and/or inorganic chemical reactions/mineralization, sediment transport, formation of a fracture, mechanical weathering or combinations of these processes. In many cases, the presence of debris forms a sub-porosity within the fracture void space. This sub-porosity often is composed of material that differs from the fracture walls in mineralogy and morphology. The sub-porosity may partially fill voids that are on the order of hundreds of microns and thereby reduce the local porosity to lengths scales on the order of sub-microns to tens of microns. It is quite clear that a sub-porosity affects fracture porosity, permeability and storativity. What is not known is how the existence/formation of a sub-porosity affects seismic wave propagation and consequently our ability to probe changes in the subsurface caused by the formation or alteration of a sub-porosity. If seismic techniques are to be developed to monitor the injection and containment of phases in sequestration reservoirs or the propping of hydraulically induced fracture to enhance oil & gas production, it is important to understand how a sub-porosity within a fracture affects macroscopic seismic and hydraulic measurements. A sub-porosity will directly affect the interrelationship between the seismic and hydraulic properties of a fracture. This reports contains the results of the three main topics of research that were performed (1) to determine the effect of a sub-porosity composed of spherical grains on seismic wave propagation across fractures, (2) to determine the effect of biofilm growth in pores and between grains on seismic wave propagation in sediment, and (3) to determine the effect of the scale of observation (field-of-view) on monitoring alteration the pore space within a fracture caused by reactive flow. A brief summary of the results for each topic is contained in the report and the full details of the research and approach are contained in the publications found in the Attachment section of this report. A list of presentation and publications of all work associated with this grant is also provided.

58 GEOSCIENCES↗

DOE OSTI2019

Radionuclide transport in brackish water through chalk fractures

Mobility of radionuclides originating from geological repositories in the subsurface has been shown to be facilitated by clay colloids. In brackish water, however, colloids may flocculate and act to immobilize radionuclides associated with them. Furthermore, little research has been conducted on radionuclide interactions with carbonate rocks. In this work, the impact of bentonite colloid presence on the transport of a cocktail of U(VI), Cs, Ce and Re through fractured chalk was investigated. Flow-through experiments were conducted with and without bentonite colloids, present as a mixture of bentonite and Ni-altered montmorillonite colloids. Ce was used as an analogue for reactive actinides in the (III) and (VI) redox states, and Re was considered an analogue for Tc. Filtered brackish groundwater (ionic strength = 170 mM) pumped from a fractured chalk aquitard in the northern Negev Desert of Israel, was used as a solution matrix. Rhenium transport was identical to that of the conservative tracer, uranine. The sorption coefficient (K d ) of U(VI), Cs and Re, calculated from batch experiments with crushed chalk, proved to be a good predictor of mass recovery in transport experiments conducted without bentonite colloids. A meaningful K d value for Ce could not be calculated due to its precipitation as a Ce-carbonate colloids. Transport of both U(VI) and Cs was indifferent to the presence of bentonite colloids. However, the addition of bentonite in the injection solution effectively immobilized Ce, decreasing its recovery from 17-41% to 0.8–1.4%. This indicates that radionuclides which interact with clay colloids that undergo flocculation and deposition may effectively be immobilized in brackish aquifers. Lastly, the results of this study have implications for the prediction of potential mobility of radionuclides in safety assessments for future geological repositories to be located in fractured carbonate rocks in general and in brackish groundwater in particular.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Impact of salinity origin on microbial communities in saline springs within the Illinois Basin, USA

Abstract Saline springs within the Illinois Basin result from the discharge of deep‐seated evaporated seawater (brine) and likely contain diverse and complex microbial communities that are poorly understood. In this study, seven saline/mineral springs with different geochemical characteristics and salinity origins were investigated using geochemical and molecular microbiological analyses to reveal the composition of microbial communities inhabiting springs and their key controlling factors. The 16S rRNA sequencing results demonstrated that each spring harbours a unique microbial community influenced by its geochemical properties and subsurface conditions. The microbial communities in springs that originated from Cambrian/Ordovician strata, which are deep confined units that have limited recharge from overlying formations, share a greater similarity in community composition and have a higher species richness and more overlapped taxa than those that originated from shallower Pennsylvanian strata, which are subject to extensive regional surface and groundwater recharge. The microbial distribution along the spring flow paths at the surface indicates that 59.8%–94.2% of total sequences in sedimentary samples originated from spring water, highlighting the role of springs in influencing microbiota in the immediate terrestrial environment. The results indicate that the springs introduce microbiota with a high biodiversity into surface terrestrial or aquatic ecosystems, potentially affecting microbial reservoirs in downstream ecosystems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Evaluating the Importance of Barometric Pumping for Subsurface Gas Transport Near an Underground Nuclear Test Site

Core Ideas Gas transport modeling can estimate ground surface breakthrough after an underground nuclear explosion. Comparison tracer‐injection experimental and model results reveals parameter sensitivity. Barometric pumping may explain gas breakthrough following an underground nuclear explosion. Amplitude and period of barometric pressure signal are key in controlling breakthrough. An underground nuclear explosion (UNE) generates and distributes radioactive gases that can be transported to the ground surface though preexisting and explosion‐induced fractures over timescales of hours to months. If detected, the presence of short‐lived radionuclides in gas is evidence of a recent UNE. Numerical modeling can provide estimates of surface arrival times that can help inform gas monitoring strategies at suspected foreign test sites. Efforts are underway at historic US UNE sites to better understand subsurface gas‐transport processes following a UNE by geologic characterization of the near‐field damage structures, field‐scale tracer experiments, and subsurface air pressure monitoring. The development of numerical models using historical and experimental datasets from former UNE sites can improve predictions by testing conceptual models, highlighting key processes, and constraining parameter ranges. The models developed in this study represent the U20az site at the Nevada National Security Site where the Barnwell device was expended in December 1989. A two‐phase (water and air), dual‐permeability flow and transport model of the U20az site was built to investigate gas transport processes under recent experimental conditions and following the Barnwell nuclear event. Results indicate that the model can explain both the lack of arrival of radioactive gas tracers in a 2013 field experiment as well as the observed arrival of radioactive gases following the 1989 Barnwell event using barometric pressure records from the respective periods, even when additional advective processes associated with the Barnwell detonation are ignored. The results demonstrate that the character of the barometric records may be a key factor in explaining the differences in transport behavior.

58 GEOSCIENCES↗

DOE OSTI2023

Radionuclide Transport Simulations Supporting Proposed Borehole Waste Disposal in Israel

A scientific collaboration between the U.S. and Israel is underway to assess the suitability of a potential site for subsurface radioactive waste disposal in the Negev Desert, Israel. The Negev Desert has several favorable attributes for geologic disposal, including an arid climate, a deep vadose zone, interlayered low-permeability lithologies, and carbonate rocks with high uranium-sorption potential. These features may provide a robust natural barrier to radionuclide migration. Geologic and laboratory characterization data from the Negev Desert are incorporated into multiphase flow and transport models, solved using PFLOTRAN, to aid in site characterization and risk analysis that will support decision-making for waste disposal in an intermediate-depth borehole design. The lithology with the greatest uranium sorption potential at the site is phosphorite. We use modeling to evaluate the ability of this layer to impact uranium transport around a proposed disposal borehole. The current objective of the simulations is focused on characterizing hypothetical leakage from waste canisters and subsequent uranium migration under three infiltration scenarios. Here, we describe a hydrogeologic model based on data from a local exploratory borehole and present results for uranium flow and transport simulations under varying infiltration scenarios. We find that under the current climate conditions, it is likely that uranium will remain in the near-field of the borehole for thousands of years. However, under a hypothesized extreme climate scenario representing an increase in infiltration by a factor of 300x above present-day values, uranium may break through the phosphorite layer and exit the base of the model domain (~200 m above the water table) within 1000 years. Simulation results have direct implications for the planning of nuclear waste disposal in the Negev Desert, and specifically in intermediate-depth boreholes.

58 GEOSCIENCES↗

DOE OSTI2016

Migration and Distribution of Natural Organic Matter Injected into Subsurface Systems - 16523

The Savannah River Site (SRS) was one of the major U.S. Department of Energy facilities that produced large amounts of radioactive and hazardous waste during the Cold War era, a legacy of material production for defense purposes. The unlined F-Area seepage basins received approximately 1.8 billion gallons of acidic waste solutions containing low-level radioactivity from numerous isotopes generated by uranium slug and irradiated fuel processing. Sufficient quantities of uranium isotopes, I-129, Tc-99, and tritium migrated into the groundwater, creating an acidic plume with a pH between 3 and 5.5. Humic substances (HS) are complex heterogeneous mixtures of poly-dispersed materials formed by biochemical and chemical reactions during the decay and transformation of plant and microbial remains. They are known for their excellent binding capacity for metals, making them a strong candidate for remediation efforts to reduce the mobility of uranium (VI) in the subsurface. Studies showed that HS are an important ion exchange and metal-complexing ligand, carrying a large number of functional groups with high complexing capacity that can greatly affect the mobility behavior of actinides in natural systems. Hence, it is a potential in situ amendment for treating groundwater contaminated with these constituents. Column experiments were conducted to address the complexities of HS in flow-through porous media representative of the acidic aquifer at the SRS F-Area Seepage Basins. The humate source used, Huma-K, is high in humic and fulvic compounds. The column studies were done using the soil obtained from SRS's FAW-1 core-hole at a depth interval of 60'-70'. A bromide tracer test was performed to obtain the pore volume (PV), variance and Peclet number. The columns were saturated with an artificial groundwater solution prepared to mimic the groundwater at SRS that was pH adjusted to 3.5 and 5 until the effluent pH reached equilibrium. Approximately one (1) PV of 5000 ppm of Huma-K was injected into the columns at 2 ml/min to study the sorption of Huma-K on the sediment. Desorption of Huma-K was studied by injecting approximately four (4) PVs of the pH adjusted artificial groundwater solution. A faster increase in pH was observed in column 1 (saturated with pH 3.5 AGW) compared to column 2 (saturated with pH 5 AGW). After injection of 2 PV of AGW, the pH of the columns reached 6.5 and 7, respectively, for columns 1 and 2. After injecting an additional PV of AGW, the pH remained relatively stable. With an increase in pH from 3.5 to 5.0, the overall retention of HA increased from 461 mg/kg for the column with pH 3.5 to 642 mg/kg for the column with pH 5. This suggests that retention is not due to adsorption alone, but adsorption coupled with precipitation and re-dissolution of Huma-K as it flows through the columns is a possible explanation for greater retention at pH 5. This negates the use of a simple Langmuir isotherm to explain humate behavior in the acidic environment of column 1. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗