Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsurface”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Vadose Zone Model for SALDS for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3 (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort: carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The parallel version of the Subsurface Transport Over Multiple Phases (STOMP1) simulator, officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP) is used to simulate flow and transport for the vadose models. The documentation for the STOMP code is comprehensive. The theoretical and numerical approaches applied in the STOMP code are documented in a published theory guide (PNNL-12030, STOMP Subsurface Transport Over Multiple Phases Version 2.0 Theory Guide). The code has undergone a rigorous verification procedure against analytical solutions, laboratory-scale experiments, and field-scale demonstrations. The application guide (PNNL-11216, STOMP Subsurface Transport Over Multiple Phases Application Guide) provides instructive examples in the application of the code to classical groundwater problems. The user’s guide (PNNL-15782, STOMP: Subsurface Transport Over Multiple Phases Version 4.0: User’s Guide) describes the general use, input file formatting, compilation, and execution of the code. The primary output of the vadose zone modeling is radionuclide transfer rates to the groundwater for input into the P2R model. The rates will be summed over the 100 by 100 m P2R grid cells that fall within the vadose zone model source domain. The Hanford Site Central Plateau was subdivided into 26 individual vadose zone models, with 13 in the 200 East Area and 13 in the 200 West Area. Waste sites that have a completed performance assessment (PA) or past-leak analysis were not included as sources of radionuclides. Instead the vadose zone to groundwater transfer rates of the Environmental Restoration Disposal Facility, Integrated Disposal Facility, US Ecology, and Waste Management Area C (WMA C) PAs and the past-leak analysis for WMA C were used as direct input to the P2R model. Each of the vadose zone models is documented in separate environmental calculation files (ECFs). This ECF describes the State-Approved Land Disposal Site (SALDS) model. The scope of this ECF is to document the development and results of the SALDS vadose zone model. CP-63515, Model Package Report: Central Plateau Vadose Zone Models, describes the approach, assumptions, process of determining the number of models required and domain of each model, input data, and processing common to all the models.

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

Spectral Induced Polarization-Biogeochemical Relationships for Remediation Amendment Monitoring

Geophysical tools such as electrical resistivity (ER) can indirectly monitor subsurface changes in response to remedial injections. These methods exhibit relatively low spatial resolution compared to sediment core characterization but are advantageous due to the ability to collect measurements non-intrusively over time across large volumes of the subsurface. Moreover, along with confirmatory groundwater or core sampling, geophysical methods can be used during active biogeochemical remedies to monitor short-term contaminant transformations and mobility, as well as part of an overall strategy for long-term monitoring of subsurface contaminated sites. The use of alternating current spectral induced polarization (SIP) provides significantly more information than conventional geophysical methods like direct current ER. SIP allows for monitoring changes in both solution and surface conductivity by separation of real and imaginary conductivity, respectively, as well as surface capacitance. In principle, SIP can measure indicators of remedy progression such as precipitation reactions that sequester contaminants, potentially providing a better indication of amendment delivery and reactivity as compared to conventional ER methods. However, multiple processes and material properties have overlapping (interacting) electrical responses across a range of frequencies. Hence, the purpose of this scoping study was to evaluate if SIP can be used to monitor (a) amendment delivery and (b) precipitation and reactivity of amendments under consideration at Hanford. SIP measurements were collected in fully saturated columns packed with sand and Hanford formation sediments containing (a) amendments that were highly conductive with significant capacitance (zero valent iron – ZVI, sulfur modified iron – SMI) and (b) amendments that exhibited low electrical conductivity with a small capacitance (calcite, apatite, bismuth). The sand was a quartz material with homogenous particle size that exhibited a relatively low surface conductivity. It was used as a control for comparison with the sediments from the Hanford Site, which have a greater surface conductivity because of their complex mineralogy and heterogeneous particle size distribution and may have complex interactions with amendments. The amendment mass fraction was varied to represent the different stages and subsurface locations associated with the delivery of a remedy. The primary objective was to identify the solution and solid surface changes associated with the delivery amendments and their secondary reactions within the subsurface. The figure below summarizes results for the amendments tested through this project. The SIP phase shift or imaginary conductivity change for the high conductivity amendments was more than 10 times that of the low conductivity amendments, highlighting the relative ease of detection of ZVI and SMI independent of the background signal from sand or Hanford sediments. The ZVI phase shift and imaginary conductivity changes occur primarily at high frequency (> 100 Hz) whereas SMI changes were at low frequency (0.01 to 10 Hz). SIP signals of SMI also increased over time and the maximum shifted to lower frequencies. The low conductivity amendments (calcite, abiotic and biotic apatite, bismuth subnitrate) exhibited relatively small phase and imaginary conductivity changes when added to sediments. The changes were above minimum detection limits (0.5 mrad for phase shift, 0.03 µS/cm for imaginary conductivity) for the highest concentration except for the commercial bismuth material. The lowest amendment concentration that can be detected is likely sediment specific, as minerals in sediments (clays, magnetite, Fe-oxides) have some capacitance and, therefore, exhibit a variable phase shift.

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

Hanford Site Composite Analysis: LLBG-200-W B Vadose Zone Model

The objectives of the vadose modeling for the updated Hanford Site Composite Analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates to the CA saturated zone model (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort: carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The parallel version of the Subsurface Transport Over Multiple Phases (STOMP 1 ) simulator, officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP), is used to simulate flow and transport for the vadose models. The documentation for the STOMP code is comprehensive. The theoretical and numerical approaches applied in the STOMP code are documented in a published theory guide (PNNL-12030, STOMP Subsurface Transport Over Multiple Phases Version 2.0 Theory Guide). The code has undergone a rigorous verification procedure against analytical solutions, laboratory-scale experiments, and field-scale demonstrations. The application guide (PNNL-11216, STOMP Subsurface Transport Over Multiple Phases Application Guide) provides instructive examples in the application of the code to classical groundwater problems. The user’s guide (PNNL-15782, STOMP: Subsurface Transport Over Multiple Phases Version 4.0: User’s Guide) describes the general use, input file formatting, compilation, and execution of the code. The primary output of the vadose zone modeling is radionuclide transfer rates to the groundwater for input into the saturated zone model. The rates will be summed over the 100 by 100 m saturated zone model grid cells that fall within the vadose zone model source domain.

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

Advanced Inversion Algorithms for Scientific Data Analysis [Slides]

Accurate subsurface characterization is crucial for all subsurface energy exploration. Accurate characterization of uncertain subsurface properties is also critical for monitoring storage of CO 2 , estimating pathways of subsurface contaminant transport, and monitoring potential nuclear explosions for treaty verification. This research will advance our world-leading subsurface sensing capabilities that are crucial for LANL missions in energy security (geothermal energy, oil/gas resource exploration, geologic carbon storage) and nuclear security (facility monitoring, detonation detection).

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Colloid-Facilitated Actinide Transport in a Cementitious-Impacted SRS Groundwater

The objective of this report is to provide a literature review relevant to the question of whether colloid facilitated transport of plutonium is enhanced in a cementitious groundwater system on the Savannah River Site (SRS). Contaminant transport is commonly described as taking place in a system with a mobile aqueous phase and an immobile solid phase. There has been an increasing awareness of a third phase, a mobile solid phase, also referred to as a mobile colloidal phase. Mobile colloids consist of organic and/or inorganic submicron-particles that move with groundwater flow. When radionuclides are associated with the mobile colloids, the net effect is that radionuclides can move faster through the subsurface system than would be predicted by transport models that do not include mobile colloids. It is important to distinguish between subsurface colloids and subsurface mobile colloids. The subsurface environment includes an enormous reservoir of colloids, but only a tiny fraction, if any, are mobile. Mobile colloid formation is commonly described as involving a three-step process: genesis, stabilization, and transport. It was concluded that there is a strong likelihood that submicron colloids of plutonium exist near the source term, (i.e., the genesis step is likely completed). While such particles have not been directly detected, it is highly likely that plutonium could either attach to submicron particles in the sediment or that submicron plutonium fragments exist at the source. However, the tendency for these plutonium colloids to move in the SRS subsurface is extremely low, especially in engineered cementitious environments where the ionic strength of the solution and the elevated concentrations of divalent cations greatly curtails colloid suspension stability (i.e., the stabilization step is likely not completed). Together these data strongly indicate that plutonium-bearing colloids would likely exist in the near field of a tank closure facility, but the colloids would be unlikely to be mobile, thereby providing a vector for enhanced transport. Under cementitious leachate impacted groundwater conditions, it is reasonable to assume plutonium transport occurs primarily as a two-phase system, a mobile aqueous phase and an immobile solid phase.

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

Assessing Suitable Geologic Carbon Storage Sites Across Utah

Utah has a wealth of potential geological reservoirs for carbon dioxide storage (CS) and a long history of geologic research resulting in an abundance of available subsurface data to evaluate CS potential. Reservoirs may include sandstone, carbonate, and basalt; these rock types are plentiful in Utah’s subsurface and the complex Phanerozoic history throughout the state requires evaluating each geologic region individually for promising reservoir-seal pairs for CO2 storage. Classifying Utah by geologic provinces (or “geo-regions”) allows for customized thinking about suitable CS reservoir and seal distribution, CO2 point sources, land use, and existing infrastructure. Preliminary results from this study highlight the geologic CS potential across 15 geo-regions. Four regions stand out as having high CS potential: the Uinta Basin, San Rafael Swell, Paradox Basin, and the southern Basin and Range Province. The Uinta Basin and San Rafael Swell geo-regions are well suited for CS and have several projects ongoing to evaluate Cretaceous Frontier and Naturita Formations, Jurassic Navajo Sandstone and Entrada Sandstone, and Permian Weber Sandstone reservoir units that lie beneath robust sealing units like the ~5000-ft-thick Mancos Shale and Carmel Formation. Reservoirs such as the Navajo and Weber Sandstones have been demonstrated to be suitable reservoirs through a long history of oil and gas exploration in Utah. New areas of interest include the southern Basin and Range in southwest Utah, where the Jurassic Navajo Sandstone is overlain by the sealing Carmel Formation at suitable depths (>3000 ft), and has good porosities based on outcrop analogue data. Just to the north (e.g., central Basin and Range), legacy wells and 2D seismic data show possible salt and subsurface basalt flows that may provide additional possible CS reservoirs and seals. This geo-region also has the advantage of being coupled with geothermal energy resources that may be used to power burgeoning direct air capture technologies. In the Paradox Basin of southeastern Utah, the Leadville Limestone is a potential storage reservoir beneath the thick (4000–5000 ft), salt-bearing Pennsylvanian Paradox Formation. Although the northern and western parts of Utah offer CS potential, these areas typically contain less infrastructure and subsurface penetrations, creating geologic uncertainty associated with subsurface seals and reservoirs due to a lack of data. Overall, this statewide assessment and ranking is the first step to aid in evaluating CS potential across Utah and provides a foundation for future research in the most favorable locations.

58 GEOSCIENCES↗
DOE OSTI2025

STILGAR End-of-Project Report

The Subsurface Tunnel Imaging LeveraGed by Analysis of Rayleigh wave ellipticity (STILGAR) project demonstrated an integrated geophysical approach for detecting, locating, and characterizing underground structural changes using dense seismic arrays and advanced inversion techniques. Field campaigns were conducted at two operational mines—the Redmond salt mine (Utah) and Graymont Pleasant Gap limestone mine (Pennsylvania)—providing real-world testbeds for monitoring anthropogenic subsurface activity. At the Redmond salt mine, seismic interferometry combined with back-projection inversion successfully identified continuous, low-amplitude signals from mining operations. The approach differentiated stationary from migrating anthropogenic sources, captured daily operational cycles, and validated the potential of passive seismic monitoring for remote detection of underground activity. At the Graymont Pleasant Gap mine, two dense seismic deployments in the spring and fall of 2023 generated over 4 TB of high-resolution data. Key outcomes included the relocation of 199 underground and 8 surface explosions with accuracies within tens of meters and the development of a 3D P-wave velocity model using the triple-difference tomography algorithm (tomoTD) that resolved major structural features such as the mine entrance, low-velocity tunnels, and roof-collapse areas. Ambient noise cross-correlation and back-projection analyses revealed persistent sources linked to ongoing mining activity, whereas horizontal-to-vertical spectral ratio (HVSR) and ellipticity studies confirmed stable site responses across seasons and identified soil thickness trends consistent with regional erosional and depositional processes. Checkerboard and sensitivity tests further validated the robustness of the tomographic results. Overall, the findings emphasize that although significant progress has been made in subsurface imaging, further work is needed to enhance the detection and localization of underground structures. Accurate imaging requires higher frequencies, yet anthropogenic sources tend to dominate the seismic record at those frequencies, and high-frequency surface waves are affected by higher modes that complicate interpretation. The improved detection and localization of human-induced signals enabled detailed temporal and spatial mapping of daily mine operations, demonstrating the feasibility of continuous anthropogenic source monitoring. Sensitivity to signals from nontraditional sources, such as fan operations, highlights the broader applicability of this approach to other industrial environments in which continuous and impulsive signals are present. The field campaigns produced a substantial volume of high-quality seismic data, supporting the development and testing of new methods for seismic source characterization and subsurface imaging. Future deployments should include sensors capable of recording lower frequencies to probe deeper structures, increase bandwidth to enhance resolution and sensitivity to both shallow and deep targets, and collect additional large-scale datasets to refine imaging and source characterization techniques. Moreover, conducting 3D modeling studies of seismic wavefields at higher frequencies will provide a better understanding of wave scattering and cavity–wavefield interactions in complex underground environments. In conclusion, the STILGAR project demonstrated that integrated seismic monitoring can effectively characterize underground operations, capturing both natural and anthropogenic signals. The approaches developed provide a foundation for improved detection, localization, and imaging of subsurface structures and are directly transferable to broader industrial monitoring applications.

58 GEOSCIENCES↗
DOE OSTI2025

Produced Water DNA Database (PW-DNA): Utilizing KBase to generate an environmental specific curated molecular database

The deep subsurface is estimated to host the majority of Earth’s microbial biomass yet remains one of the most challenging environments to access and study. One common approach to investigate these microbial communities is through the analysis of produced water from subsurface reservoirs, where researchers can assess water and gas chemistry along with molecular (DNA/RNA) sequence data. Advances in high-throughput sequencing have greatly expanded our understanding of these environments and their biotechnological potential. However, further progress requires large-scale, integrative meta-analyses across diverse datasets. To address this need, we developed the Produced Water-DNA (PW-DNA) Database, a curated, publicly available resource that consolidates microbial DNA/RNA sequences, geochemical data, and relevant metadata from in situ hydrocarbon environments such as coal beds, oil reservoirs, and natural gas systems. The PW-DNA database delivers three core benefits to the research community: (1) it improves data sharing by linking environmental microbial datasets with corresponding geochemical parameters, enabling more robust filtering and analysis; (2) it connects with complementary research databases to promote broader dissemination and interoperability; and (3) it supports technological innovation by serving as a resource for identifying microbial trends and exploring genetic potential. While individual studies have highlighted basin-specific microbial communities and functional redundancy in biogeochemical cycling, a comprehensive, system-wide perspective is needed to better understand connectivity and novelty across subsurface ecosystems. By designing the PW-DNA in the KBase platform, we provide a reproducible, visual framework for integrating large-scale genomic and geochemical data, enabling researchers to perform more informed analyses and experimental design. Ultimately, this resource enhances the ability to identify, characterize, and interpret microbial functions across diverse subsurface environments, thereby accelerating discovery in subsurface microbiology and biotechnology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Numerical Modeling of Anisotropic Elastic-Wave Sensitivity Propagation for Optimal Design of Time-Lapse Seismic Surveys

Reliable subsurface time-lapse seismic monitoring is crucial for many geophysical applications, such as enhanced geothermal system characterization, geologic carbon utilization and storage, and conventional and unconventional oil/gas reservoir characterization, etc. In this work, we develop an elastic-wave sensitivity propagation method for optimal design of cost-effective time-lapse seismic surveys considering the fact that most of subsurface geologic layers and fractured reservoirs are anisotropic instead of isotropic. For anisotropic media, we define monitoring criteria using qP- and qS-wave sensitivity energies after decomposing qP- and qS-wave components from the total elastic-wave sensitivity wavefield using a hybrid time- and frequency-domain approach. Geophones should therefore be placed at locations with significant qP- and qS-wave sensitivity energies for cost-effective time-lapse seismic monitoring in an anisotropic geology setting. Our numerical modeling results for a modified anisotropic Hess model demonstrate that, compared with the isotropic case, subsurface anisotropy changes the spatial distributions of elastic-wave sensitivity energies. Consequently, it is necessary to consider subsurface anisotropies when designing the spatial distribution of geophones for cost-effective time-lapse seismic monitoring. This finding suggests that it is essential to use our new anisotropic elastic-wave sensitivity modeling method for optimal design of time-lapse seismic surveys to reliably monitor the changes in subsurface reservoirs, fracture zones or target monitoring regions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Development of inter-grid-cell lateral unsaturated and saturated flow model in the E3SM Land Model (v2.0)

Abstract. The lateral transport of water in the subsurface is important in modulating terrestrial water energy distribution. Although a few land surface models have recently included lateral saturated flow within and across grid cells, it is not a default configuration in the Climate Model Intercomparison Project version 6 experiments. In this work, we developed the lateral subsurface flow model within both unsaturated and saturated zones in the Energy Exascale Earth System Model (E3SM) Land Model version 2 (ELMv2.0). The new model, called ELMlat, was benchmarked against PFLOTRAN, a 3D subsurface flow and transport model, for three idealized hillslopes that included a convergent hillslope, divergent hillslope, and tilted V-shaped hillslope with variably saturated initial conditions. ELMlat showed comparable performance against PFLOTRAN in terms of capturing the dynamics of soil moisture and groundwater table for the three benchmark hillslope problems. Specifically, the mean absolute errors (MAEs) of the soil moisture in the top 10 layers between ELMlat and PFLOTRAN were within 1 %±3 %, and the MAEs of water table depth were within ±0.2 m. Next, ELMlat was applied to the Little Washita experimental watershed to assess its prediction of groundwater table, soil moisture, and soil temperature. The spatial pattern of simulated groundwater table depth agreed well with the global groundwater table benchmark dataset generated from a global model calibrated with long-term observations. The effects of lateral groundwater flow on the energy flux partitioning were more prominent in lowland areas with shallower groundwater tables, where the difference in simulated annual surface soil temperature could reach 0.3–0.4 ∘C between ELMv2.0 and ELMlat. Incorporating lateral subsurface flow in ELM improves the representation of the subsurface hydrology, which will provide a good basis for future large-scale applications.

58 GEOSCIENCES↗
DOE OSTI2022

Multi-frequency electrical impedance tomography

Apparatus includes a plurality of geological subsurface electrical line sensors spaced apart from each other proximate a predetermined geological subsurface region of interest, with at least one of the electrical line sensors situated as a line source to produce a multi-frequency electrical impedance tomography source signal, and with at least one of the electrical line sensors situated as a line detector to receive the multi-frequency electrical impedance tomography response signal associated with the source signal that propagates through the predetermined geological subsurface region of interest, and a controller including a processor and a memory configured with instructions that, when executed by the processor, cause the processor to determine an electrical mapping over the predetermined geological subsurface region of interest based on the multi-frequency electrical impedance tomography source signal, response signal, and the spatial positions of the geological subsurface electrical line sensors.

Karra, Satish↗
DOE OSTI2024

Evolution of a Geological Model for Co-Producing Electricity at the Blackburn Oil Field, Nevada: Preprint

Oil fields characterized by significant proportions of hot water in the reservoir offer a unique opportunity for geothermal energy co-production and utilization. Development of an accurate geological model is essential in understanding the subsurface. However, this is a challenge at the preliminary stages of resource assessment when data on the subsurface parameters and subterranean fault network is scarce. Our focus in this paper is to explore a methodology for developing an accurate conceptual geological model of the complexly faulted reservoir at the Blackburn Field, Nevada, to create a numerical simulation that will eventually be used to assess the potential of co-producing electricity. There are four steps to the evolution of this model: (1) collection of public data and research on the geological setting of the system; (2) digitization of seismic time-interpretation and cross-sections to create the most representative model of the subsurface; and (3) creation of a natural-state numerical model to simulate the thermal gradient through the subsurface. Utilizing publicly available data of the geological system and Leapfrog conceptual modeling software, we refined a preliminary model to be more representative of the subsurface geology of an active geothermal exploration stage project. The conceptual geological model developed has been imported into VOLSUNG reservoir modeling software to create a natural-state model, further characterizing the Blackburn Field reservoir's thermal gradient. In the future, new seismic data will be utilized to further evolve the conceptual geological model for numerical simulations.

Basin and Range province↗
DOE OSTI2021

Modeling framework for evaluating the impacts of hydrodynamic pressure on hydrologic exchange fluxes and residence time for a large-scale river section over a long-term period

Quantifying hydrologic exchange fluxes (HEF) at the river and subsurface interface and their residence times (RT) in subsurface are important for managing the water quality and ecosystem health in dynamic river corridor systems. In this study, a modeling framework is developed for coupling the three-dimensional (3D) multi-phase surface, subsurface flow transport, and numerical tracer model for RT in a large-scale river section over a long period. The framework is utilized to evaluate the impacts of hydrodynamic pressure on HEFs and RT for a 30 km section of the Columbia River in Washington State over a three-year period. Based on comparisons between model simulations with and without considering hydrodynamic pressure, we found that hydrodynamic pressure increases the net HEFs by 7% with river gaining water from the subsurface domain, and leads to slight reduction of RT.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Kilometer-Scale Hydrologic Exchange Flows in a Gravel Bed River Corridor and Their Implications to Solute Migration

In this work, a well-characterized field site along a major river corridor was used to investigate the dynamic pathways and impacts of subsurface hydrogeologic structures on kilometer-scale hydrologic exchange flows between river water and groundwater. A aqueous uranium (U aq ) plume exists at the site that discharges to the Columbia River. We performed temporally intensive monitoring of specific conductance (SpC) and U aq concentrations within the plume at increasing distances from the river shoreline and at locations both within and outside a presumed subsurface pathway of lateral hydrologic exchange. Specific conductance and U aq were utilized as in-situ tracers of hydrologic exchange and associated groundwater-surface water mixing. Monitoring data was collected for an approximate two year period during which seasonal river elevations varied by more than 2 m causing distinct events of river water intrusion and retreat from the near-shore aquifer, and large temporal variations in well-water analyte concentrations. Statistical clustering analysis was applied to identify well groupings exhibiting common time series behaviors. PFLOTRAN was implemented within a modeling domain of 3.3 x 1.9 km with dynamic flow boundary conditions and accurate subsurface topography for the basal aquitard underlying the aquifer. In-silico calculations of river water tracer behavior were performed to assess predominant flow directions and velocities influenced by river stage variations and variable aquitard topography, and their collective influences on SpC and U aq concentrations. We demonstrate complicated multidirectional flow behaviors at the kilometer scale resulting from hydrologic exchange and subsurface structure that influence plume dynamics and determine the overall system behavior.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The transformation of U(VI) and V(V) in carnotite group minerals during dissimilatory respiration by a metal reducing bacterium

Recent results from laboratory and field studies support that dissimilatory metal reducing (DMR) bacteria influence the fate and transport of uranium in anaerobic subsurface environments. To date, most research efforts have focused on the reduction of soluble U(VI) by DMR bacteria to form insoluble uraninite (UO 2 ). Subsurface environments harbor, however, large reservoirs of U(VI) in solid or mineral form. Uranium that is structure-bound in minerals is expected to be more refractory to microbial reduction than soluble U, based on analogy with Fe respiration. The reducibility of U(VI) could impact the fate of U(IV) by controlling mineral precipitation reactions, which has implications for the long-term immobilization of U in subsurface environments. Here, we studied anoxic cultures of Shewanella putrefaciens CN32 incubated with natural carnotite-group minerals by X-ray diffraction, electron microscopy, scanning transmission X-ray microscopy (STXM). Near-edge X-ray absorption fine structure (NEXAFS) spectroscopy measurements at U-N 4,5 , V-L 2,3 , and O-K edges on cultures incubated up to 10 months show that V(V) was reduced to V(IV), whereas U was not reduced. In contrast, V(V) and U(VI) in solution were both completely reduced to lower oxidation states by CN32, as precipitates within the exopolymer surrounding the bacteria. Assays for the toxicity of U and V to CN32 showed that biofilm formation was stimulated at 0.001 M U(VI), and growth was inhibited at concentrations of U(VI) greater than 0.001 M. Vanadium did not inhibit growth or stimulate biofilm formation at any concentration tested. Investigations of the bacteria-mineral and bacteria-metal interface at the nanometer and molecular scales provide new insights into the co-respiration of V and U that help explain their biogeochemical cycling and have implications for subsurface bioremediation of these elements.

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

Insights from a coupled thermo-hydro-mechanical analysis of a layered high-temperature thermal energy storage reservoir

Coupled thermal-hydraulic-mechanical (THM) modeling is applied to investigate the performance of a seasonal high-temperature aquifer thermal energy storage operation based on data and conditions from current site investigations at the Geostorage Forsthaus pilot project in Bern (Switzerland). The model includes subhorizontal sand lenses of various lengths and dips that are embedded in a low permeability clay matrix. Thermal energy storage is simulated by seasonal injection and withdrawal of hot (up to 90 °C) water from a main well, with reservoir pressure regulated by two auxiliary wells at a distance of about 70 m from the main well. The results show how targeted injection into deeper permeable storage formations, along with active deep well pressure control, can effectively minimize geomechanical impact and the potential risk of damaging subsurface storage and sealing formations, or even surface facilities. With such pressure control, the subsurface mechanical responses are dominated by thermal strain and stress, which can be monitored with subsurface fiber optics. The study demonstrates how coupled THM modeling can be applied for the design of a safe and efficient thermal energy storage operation, and how subsurface fiber optic monitoring can be applied for performance confirmation, allowing for more confident operational forecasting.

Rutqvist, Jonny↗
DOE OSTI2021

Stratigraphic Identification with Airborne Electromagnetic Methods at the Hanford Site, Washington

Stratigraphic units can influence the fate and transport of subsurface contaminants within groundwater. Units having coarse-grained sediments act as preferential flow pathways, and therefore can accelerate the transport of contaminants to reach human and ecological receptors. At legacy waste sites, detailed knowledge of subsurface stratigraphy can be used for effective monitoring and remediation planning to help minimize risk to human health and the environment. Airborne electromagnetic (AEM) methods can non-invasively provide information on kilometer-scale or larger subsurface stratigraphic features and fill informational gaps in directly sampled data from sparsely located boreholes. In this paper, we present inversion results of a 412 line-km frequency-domain AEM survey to delineate subsurface stratigraphic features at the Hanford Site, located in southeastern Washington State. The inversion was performed using a massively parallel 3D electromagnetic modeling and inversion code, where the modeling is based on solving frequency-domain Maxwell’s equations using an unstructured-mesh finite-element method and the inversion employs a Gauss-Newton optimization scheme. The results are compared to an underlying geologic framework model (GFM), built by interpolating contact depths of stratigraphic units interpreted from site borehole datasets. In areas with good borehole coverage, the inversion results show a good match with the GFM to a depth of about 60 m. Outside of these areas, the inversion results exhibit inconsistencies from the assumptions made to create the GFM, demonstrating that the AEM survey results can be used to improve the understanding of the geological conceptual model.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Self-Assembly of Silica Nanoparticles at Water–Hydrocarbon Interfaces: Insights from In Operando Small-Angle X-ray Scattering Measurements and Molecular Dynamics Simulations

To achieve tunable controls on the interactions of siliceous materials in subsurface environments bearing water and hydrocarbons, it is essential to determine the influence of hydrocarbon–water interfaces on the self-assembly of the silica particles. The hydrophilic silica particle have a tendency to aggregate on the water front of the interface. Self-assembly of silica nanoparticles proceeds via the migration of these nanoparticles to the water–hydrocarbon interface followed by aggregation at the interface. Fractal-like morphologies of assembled silica nanoparticles at water–hydrocarbon interfaces are observed. Rapid assembly of the hydrophilic silica nanoparticles at water–hydrocarbon interfaces corresponds to an overall reduction in the surface tension of water–toluene and water–heptane systems. These studies demonstrate the silica aggregation that is the precursor of silica polymerization, or nucleation and growth is influenced by the presence of hydrocarbons in subsurface geologic environments. These insights were derived from in operando ultrasmall and small-angle X-ray scattering (USAXS/SAXS) measurements, cryo-scanning electron microscopy (Cryo-SEM) imaging, and classical molecular dynamics simulations. As a result, these studies are intended to inform current and future efforts aimed at tuning silica reactivity in subsurface geologic environments to enhance permeability, the design and use of silica-based proppants to enhance fractures, and the development of effective strategies to control silica-based scaling behavior in subsurface reservoirs.

42 ENGINEERING↗