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 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS2008

Design of a Mars Airplane Propulsion System for the Aerial Regional-Scale Environmental Survey (ARES) Mission Concept

The Aerial Regional-Scale Environmental Survey (ARES) is a Mars exploration mission concept that utilizes a rocket propelled airplane to take scientific measurements of atmospheric, surface, and subsurface phenomena. The liquid rocket propulsion system design has matured through several design cycles and trade studies since the inception of the ARES concept in 2002. This paper describes the process of selecting a bipropellant system over other propulsion system options, and provides details on the rocket system design, thrusters, propellant tank and PMD design, propellant isolation, and flow control hardware. The paper also summarizes computer model results of thruster plume interactions and simulated flight performance. The airplane has a 6.25 m wingspan with a total wet mass of 185 kg and has to ability to fly over 600 km through the atmosphere of Mars with 45 kg of MMH / MON3 propellant.

Kuhl, Christopher A.↗

NASA NTRS2005

Aerobic Anoxygenic Phototrophic Bacteria in the Mid-Atlantic Bight and the North Pacific Gyre

The abundance of aerobic anoxygenic phototrophic (AM) bacteria, cyanobacteria and heterotrophs was examined in the Mid-Atlantic Bight and the central North Pacific gyre using infrared fluorescence microscopy coupled with image analysis and flow cytometry. AAP bacteria comprised 5% to 16% of total prokaryotes in the Atlantic but only 5% or less in the Pacific. In the Atlantic, AAP bacterial abundance was as much as 2-fold higher than Prochlorococcus and 10-folder higher than Synechococcus. In contrast, Prochlorococcus outnumbered AAP bacteria 5- to 50-fold in the Pacific. In both oceans, subsurface abundance maxima occurred within the photic zone, and AAP bacteria were least abundant below the 1% light depth. Concentrations of bacteriochlorophyll a (BChl a) were low (approx.1%) compared to chlorophyll a. Although the BChl a content of AAP bacteria per cell was typically 20- to 250-fold lower than the divinyl-chlorophyll a content of Prochlorococcus, in shelf break water the pigment content of AAP bacteria approached that of Prochlorococcus. The abundance of AAP bacteria rivaled some groups of strictly heterotrophic bacteria and was often higher than the abundance of known AAP genera (Erythrobacter and Roseobacter spp.). The distribution of AAP bacteria in the water column, which was similar in the Atlantic and the Pacific, was consistent with phototrophy.

Cottrell, Matthew T.↗

DOE OSTI2023

Closed Loop Geothermal Working Group: GeoCLUSTER App, Subsurface Simulation Results, and Publications

To better understand the heat production, electricity generation performance, and economic viability of closed-loop geothermal systems in hot-dry rock, the Closed-Loop Geothermal Working Group -- a consortium of several national labs and academic institutions has tabulated time-dependent numerical solutions and levelized cost results of two popular closed-loop heat exchanger designs (u-tube and co-axial). The heat exchanger designs were evaluated for two working fluids (water and supercritical CO2) while varying seven continuous independent parameters of interest (mass flow rate, vertical depth, horizontal extent, borehole diameter, formation gradient, formation conductivity, and injection temperature). The corresponding numerical solutions (approximately 1.2 million per heat exchanger design) are stored as multi-dimensional HDF5 datasets and can be queried at off-grid points using multi-dimensional linear interpolation. A Python script was developed to query this database and estimate time-dependent electricity generation using an organic Rankine cycle (for water) or direct turbine expansion cycle (for CO2) and perform a cost assessment. This document aims to give an overview of the HDF5 database file and highlights how to read, visualize, and query quantities of interest (e.g., levelized cost of electricity, levelized cost of heat) using the accompanying Python scripts. Details regarding the capital, operation, and maintenance and levelized cost calculation using the techno-economic analysis script are provided. This data submission will contain results from the Closed Loop Geothermal Working Group study that are within the public domain, including publications, simulation results, databases, and computer codes. GeoCLUSTER is a Python-based web application created using Dash, an open-source framework built on top of Flask that streamlines the building of data dashboards. GeoCLUSTER provides users with a collection of interactive methods for streamlining the exploration and visualization of an HDF5 dataset. The GeoCluster app and database are contained in the compressed file geocluster_vx.zip, where the "x" refers to the version number. For example, geocluster_v1.zip is Version 1 of the app. This zip file also contains installation instructions. **To use the GeoCLUSTER app in the cloud, click the link to "GeoCLUSTER on AWS" in the Resources section below. To use the GeoCLUSTER app locally, download the geocluster_vx.zip to your computer and uncompress this file. When uncompressed this file comprises two directories and the geocluster_installation.pdf file. The geo-data app contains the HDF5 database in condensed format, and the GeoCLUSTER directory contains the GeoCLUSTER app in the subdirectory dash_app, as app.py. The geocluster_installation.pdf file provides instructions on installing Python, the needed Python modules, and then executing the app.

15 GEOTHERMAL ENERGY↗

NASA NTRS2001

Seismic Study of the Dynamics of the Solar Subsurface from SoHO Observations

In collaboration with Dr. Baudin, we have developed and refined the new observational methodology for local helioseismology known as time-distance analysis. Global helioseismology study the solar oscillations as a superposition of resonant modes, whose properties (mode frequencies) reflect the global structure of the sun (sound speed stratification, rotation rate, etc). In contrast, local helioseismology look at the solar oscillations as wave packets whose propagation will be affected by perturbations of the media sampled. These local perturbations (sound speed or velocity flows) will modify the propagation time, that in turn can be used as a diagnostic tool for a given region. From a data reduction perspective, the processing of solar dopplergrams that result in time-distance maps, i.e. propagation times as a function of distance between bounces at the surface, is radically different from the methodology used for global mode analysis. We have, in a first step, further develop the programs needed to carry out such analysis. We have then applied them to NMI data set, and explore the trade-off between various averaging and filtering approaches - steps required to improve the signal-to-noise ratio of correlation maps - and the resulting stability and precision of the fitted propagation times. While excessive averaging (whether over space, propagation distance, or time) will reduce the diagnostic potential of the method, insufficient averaging lead to unstable fits, or uncertainties so large as to hide the information we seek. In a second phase, we have developed the analysis methodology required to infer local properties from perturbation in time propagation. Namely, we have developed time-distance inversion techniques, with an emphasis on inferences of velocity flows from time anomalies. Note also that during the period covered by this grant, all the investigators on this proposal (i.e., Drs. Baudin, Eff-Darwich, Korzennik, and Noyes) took part in the organization of the SOHO 6 /GONG 99 Workshop: Structure and Dynamics of the Interior of the Sun and Sunlike Stars, held on June 1-4 1999 at the Boston Park Plaza Hotel in Boston, Massachusetts, USA. it was very well attended by more than 160 participants from 26 countries from all over the world. The proceedings were published in two volumes as ESA SP-418, with Sessions I-III in Volume 1, and Sessions IV-VI in Volume 2 (1,000 pages in total). The complete contents are also included in digital form on a CD-ROM included with Volume 1. This CD-ROM also contains additional multi-media material that complements some of the contributions.

Korzennik, Sylvain G.↗

DOE OSTI2023

Developing Methods to Assess Changes in Mechanical Properties of Shale Modified by Engineered Mineral Precipitation

Fractures in subsurface shale formations serve multiple purposes, for example, in the recovery of resources in hydraulic fracturing or as potential harmful leakage passages through caprocks that may contribute undesired fluids to the atmosphere or functional groundwater aquifers. A proposed method to seal or influence fracture properties is Ureolysis-Induced Calcium Carbonate Precipitation (UICP), a bio-mineralization technology driven by the enzymatic hydrolysis of urea, resulting in the formation of calcium carbonate. Sporosarcina pasteurii is a common microbe used as the source of the urease enzyme that catalyzes the chemical reaction. The resulting calcium carbonate can bridge the gaps in fractured shale and reduce fluid flow through fractures. However, there is little information on how this process affects the mechanical properties of the resulting biomineralized shale. This study represents the first step toward determining the influence of UICP treatment on shale material and its subsequent mechanical strength properties. This methods development study aims to determine the effect that temperature has on the tensile strength of intact, unfractured shale cores (2.54 cm (1 in) diameter, 5.08 cm (2 in) long). Tensile strength was determined indirectly using a modified Brazilian test where the splitting tensile strength is attained by applying a compressive load onto the core. Shale cores from Eagle Ford and Wolfcamp formations were tested at both room temperature and 60°C to determine if increased temperature influences the tensile strength of the rock. This data will help to assess the necessity of testing biomineralized cores at temperature. Though 60°C may not mimic subsurface temperatures of the shales used in this study, it was chosen due to limitations of the UICP process while still approaching temperatures of shale formations. This project aims to evaluate what effect temperature has on the mechanical properties of intact shale cores so that engineered or natural rock fractures that are sealed by biomineralization can be better understood.

Bedey, Kayla↗

DOE OSTI2021

Striving to translate shale physics across ten orders of magnitude: What have we learned?

Shales will play an important role in the successful transition of energy from fossil-based resources to renewables in the coming decades. Aside from being a significant source of low-carbon intensity fuels, like natural gas, they also serve as geologic seals of subsurface formations that may be used to isolate nuclear waste, sequester CO 2 , or store intermittent energy (e.g., solar hydrogen). Despite their importance, shales pose significant engineering and environmental challenges due to their nanoporous structure and extreme heterogeneity that spans at least ~10 orders of magnitude in spatial scale. Two challenges inhibit a system-level understanding: (1) the physics of fluid flow and phase behavior in shales are poorly understood due to the dominant molecular interactions between minerals and fluids under confinement, and (2) the apparent lack of scale separation that prevents a reliable (closed) description of the physics at any single scale of observation. In this review, we focus on the latter issue and discuss scale translation, which in its broadest sense is transforming data or simulations from one spatiotemporal scale to another. While effective scale translation is not exclusive to shales, but all geologic porous media, the need for it is especially acute in shales given their high degree of heterogeneity. Classical theories like homogenization, while indispensable, fail when scales are not separated. Other methods, like numerical upscaling, scale-translate in only one direction: small to large, but not the reverse, called downscaling. However, the confluence of advances in three areas are bringing challenging problems such as shales within reach: increased computational power and scalable algorithms; high-resolution imaging and multi-modal data acquisition; and machine learning to process massive amounts of data. While these advances equip geoscientists with a wide array of experimental and computational tools, no individual tool can probe the entire gamut of heterogeneity in shales. Their effective use, therefore, requires an ability to bridge between various data types obtained at different scales. The aim of this review is to present a coherent account of computational and experimental methods that may be used to achieve just that, i.e., to perform scale translation. We provide a broader definition of scale translation, one that transcends classical homogenization and upscaling methods, but is consistent with them and accommodates notions like downscaling and data translation. After a brief introduction to homogenization, we review hybrid methods, numerical upscaling and its recent extensions, multiscale computing, high-resolution imaging, and machine learning. We place particular emphasis on multiscale computing and propose an algorithmic framework to bridge between the pore (micro) and Darcy (macro) scales. Throughout the paper, we draw comparisons between the various methods and highlight their (often hidden) similarities, differences, benefits, and pitfalls. We finally conclude with two case studies on shales that exemplify some of the methods presented.

58 GEOSCIENCES↗

DOE OSTI2018

Dual in-aquifer and near surface processes drive arsenic mobilization in Cambodian groundwaters

Millions of people globally, and particularly in South and Southeast Asia, face chronic exposure to arsenic from reducing groundwaters in which. Arsenic release to is widely attributed largely to reductive dissolution of arsenic-bearing iron minerals, driven by metal reducing bacteria using bioavailable organic matter as an electron donor. However, the nature of the organic matter implicated in arsenic mobilization, and the location within the subsurface where these processes occur, remains debated. In a high resolution study of a largely pristine, shallow aquifer in Kandal Province, Cambodia, we have used a complementary suite of geochemical tracers (including 14 C, 3 H, 3 He, 4 He, Ne, δ 18 O, δD, CFCs and SF 6 ) to study the evolution in arsenic-prone shallow reducing groundwaters along dominant flow paths. The observation of widespread apparent 3 H- 3 He ages of <55years fundamentally challenges some previous models which concluded that groundwater residence times were on the order of hundreds of years. Surface-derived organic matter is transported to depths of >30m, and the relationships between age-related tracers and arsenic suggest that this surface-derived organic matter is likely to contribute to in-aquifer arsenic mobilization. A strong relationship between 3 H- 3 He age and depth suggests the dominance of a vertical hydrological control with an overall vertical flow velocity of ~0.4±0.1m·yr-1 across the field area. A calculated overall groundwater arsenic accumulation rate of ~0.08±0.03μM·yr -1 is broadly comparable to previous estimates from other researchers for similar reducing aquifers in Bangladesh. Although apparent arsenic groundwater accumulation rates varied significantly with site (e.g. between sand versus clay dominated sequences), rates are generally highest near the surface, perhaps reflecting the proximity to the redox cline and/or depth-dependent characteristics of the OM pool, and confounded by localized processes such as continued in-aquifer mobilization, sorption/desorption, and methanogenesis.

54 ENVIRONMENTAL SCIENCES↗

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 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 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 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 OSTI2025

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

60 APPLIED LIFE SCIENCES↗

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↗