SEARCH · Engineering Papers

Results for “subsurface system”

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 505 records · Page 28

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Detrital U-Pb zircon and 40 Ar/ 39 Ar muscovite geochronology from Middle Pennsylvanian strata in the Anadarko Basin, Texas Panhandle, USA

The Late Mississippian-Permian Anadarko Basin formed in Texas and Oklahoma, USA as the result of inversion of Neoproterozic and Cambrian rift structures. Subsidence was driven by flexural loading of the Amarillo-Wichita Uplift, and this uplift may represent the easternmost element of the Ancestral Rocky Mountains system. The northwestern part of this basin has generally been interpreted to have been filled by sediment derived from the Ancestral Front Range Uplift, ~475 km to the northwest during the early stages of basin filling. We test this model using U-Pb detrital zircon and 40 Ar/ 39 Ar detrital muscovite results from three subsurface samples of the Morrow B sandstone in the northwestern part of the Anadarko Basin. We provide a new maximum depositional age of 310.9 ± 4.9 Ma that indicates the age of the Morrow B to be no older than late Atokan to early Desmoinesian Age (Moscovian), ~10 Myr younger than previously interpreted. In contrast to some previous interpretations, we propose that the most likely source for the sediment in the Morrow B is the Amarillo Uplift to the south. Detrital zircon and detrital muscovite data have age peaks at 900–1300 Ma, 1370 Ma and 1600–1800 Ma corresponding to derivation from Grenville, Granite-Rhyolite and Yavapai-Mazatzal basement provinces, respectively. A dominant detrital zircon peak at 1370 Ma suggests that Mesoproterozoic granites in the Amarillo Uplift were exposed by Middle Pennsylvanian time, and detritus eroded from the Amarillo Uplift dominated the lowstand sediment budget of the Texas Panhandle during this time; small volumes of sediment were likely sourced from the Ancestral Front Range. In conclusion, this study presents the first detrital geochronology data from the subsurface Anadarko Basin and the first detrital muscovite data from late Paleozoic southwestern Laurentia. The results presented here highlight the interpretive power of combined detrital zircon and muscovite datasets.

58 GEOSCIENCES↗

DOE OSTI2022

The impact of freeze-thaw history on soil carbon response to experimental freeze-thaw cycles

This dataset contains data used for the paper "The impact of freeze-thaw history on soil carbon response to experimental freeze-thaw cycles". The Related References field will be updated with a full citation when available.Freeze-thaw is a disturbance process in cold regions where permafrost soils are becoming vulnerable to temperature fluctuations above 0˚C. Freeze-thaw alters soil physical and biogeochemical properties with implications for carbon persistence and emissions in Arctic landscapes. We examined whether different freeze-thaw histories in two soil systems led to contrasting biogeochemical responses under a laboratory-controlled freeze-thaw incubation. We investigated controls on soil organic matter (SOM) composition through Fourier-transform ion cyclotron resonance mass spectroscopy (FT-ICR-MS) to identify nominal carbon oxidation states and relative abundances of aliphatic-type carbon molecules in both surface and subsurface soils. Soil cores (~ 60 cm-depth) were sampled from two sites in Alaskan permafrost landscapes with different in situ freeze-thaw characteristics: Healy (>40 freeze-thaw cycles annually) and Toolik (<150 freeze-thaw cycles annually). FT-ICR-MS was coupled with in situ temperature data and soil properties (i.e., soil texture, mineralogy) to assess (1) differences in SOM composition associated with previous freeze-thaw history and (2) sensitivity to experimental freeze-thaw in the extracted cores. Control (freeze-only) samples showed greater carbon oxidation in Healy soils compared with Toolik, even in lower mineral horizons where freeze-thaw history was comparable across both sites. Healy showed the most loss of carbon compounds following experimental freeze-thaw in the lower mineral depths, including a decrease in aliphatics. Toolik soils responded more slowly to freeze-thaw as shown by intermediary carbon oxidation distributed across multiple carbon compound classes. Variations in the response of permafrost carbon chemistry to freeze-thaw is an important factor for predicting changes in soil function as permafrost thaws in high northern latitudes.This dataset contains a compressed (.zip) archive of the data and R scripts used for this manuscript. The dataset includes files in .csv format, which can be accessed and processed using MS Excel or R. This archive can also be accessed on GitHub at https://github.com/Erin-Rooney/FTC-FTICR (DOI 10.5281/zenodo.6533064).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs: Preprint

Subsurface data analysis, reservoir modeling, and machine learning (ML) techniques have been applied to the Brady Hot Springs (BHS) geothermal field in Nevada, USA to further characterize the subsurface and assist with optimizing reservoir management. Hundreds of reservoir simulations have been conducted in TETRAD-G and CMG STARS to explore different injection and production fluid flow rates and allocations and to develop a training data set for ML. This process included simulating the historical injection and production since 1979 and prediction of future performance through 2040. ML networks were created and trained using TensorFlow based on multilayer perceptron (MLP), long short-term memory (LSTM), and convolutional neural network (CNN) architectures. These networks took as input selected flow rates, injection temperatures, and historical field operation data and produced estimates of future production temperatures. This approach was first successfully tested on a simplified single fracture doublet system, followed by the application to the BHS reservoir. Using an initial BHS dataset with 37 simulated scenarios, the trained and validated network predicted the production temperature for 6 production wells with the mean absolute percentage error of less than 8%. In a complementary analysis effort, the principal component analysis applied to 13 BHS geological parameters revealed that vertical fracture permeability shows the strongest correlation with fault density and fault intersection density. A new BHS reservoir model was developed considering the fault intersection density as proxy for permeability. This new reservoir model helps to explore under-exploited zones in the reservoir. A data gathering plan to obtain additional subsurface data was developed; it includes temperature surveying for three idle injection wells, at which the reservoir simulations indicate high bottom-hole temperatures. The collected data assist with calibrating the reservoir model and may lead to converting these wells to producers to access under-exploited zones in the reservoir. Data gathering activities are planned for the first quarter of 2021.

40 EE - Geothermal Technologies Office (EE-4G)↗

DOE OSTI2024

Pilgrim Hot Springs: GEOPHIRES Inputs and Outputs for Direct-Use Geothermal District Heating and Cooling

This dataset includes files for a techno-economic analysis conducted using the GEOPHIRES simulator to examine the feasibility of expanding a larger district heating site in a remote location: Pilgrim Hot Springs, Alaska. Files included here are GEOPHIRES inputs and outputs for five different scenarios with varying demand, cycle, and system design characteristics to analyze. Also included is the link to the GEOPHIRES GitHub, as well as a link to the dataset that contains the energy modelling used to determine the heating demand for the district. For a list of the differences between scenarios, see the included "Input Overview.txt" file. Fields included in the input files are: subsurface technical parameters, surface technical parameters, financial parameters, capital and O&M parameters, as well as simulation parameters. The output files are case reports that summarize all equipment, reservoir characteristics, costs, and heating profiles.

15 GEOTHERMAL ENERGY↗

DOE OSTI2025

Representing lateral groundwater flow from land to river in Earth system models

Lateral groundwater flow (LGF) is an important hydrologic process in controlling water table dynamics. Due to the relatively coarse spatial resolutions of land surface models, the representation of this process is often overlooked or overly simplified. In this study, we developed a hillslope-based lateral groundwater flow model. Specifically, we first developed a hillslope definition model based on an existing watershed delineation model to represent the subgrid spatial variability in topography. Building upon this hillslope definition, we then developed a physical-based lateral groundwater flow using Darcy’s equation. This model explicitly considers the relationships between the groundwater table along the hillslope and the river water table levels. We coupled this intra-grid model to the land component (E3SM Land Model: ELM) and river component (MOdel for Scale Adaptive River Transport: MOSART) of the Energy Exascale Earth System Model (E3SM). We tested both the hillslope definition model and the lateral groundwater flow model and performed sensitivity experiments using different configurations. Simulations for a single grid cell at 0.5°×0.5° within the Amazon basin show that the definition of hillslope is the key to modeling lateral flow processes and the runoff partition between surface and subsurface can be dramatically changed using the hillslope approach. Although our method provides a pathway to improve the lateral flow process, future improvements are needed to better capture the subgrid structure to account for the spatial variability in hillslopes within the simulated grid of land surface models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Microbial Community Composition in Deep‐Subsurface Reservoir Fluids Reveals Natural Interwell Connectivity

Abstract The identification of natural fractures and the wells they connect is crucial for the development of geological reservoirs because it may have an important impact on reservoir model construction and hydraulic fracture propagation. In this study we investigated the use of a novel data source, the microbial community composition in the reservoir formation fluids, for identification of interwell connectivity caused by natural fractures. We verified this concept at a newly developed mesoscale enhanced geothermal system (EGS) field testbed located 4,850 ft (1,478 m) beneath ground surface at the Sanford Underground Research Facility in Lead, SD. Fluids produced at or near the EGS testbed were sampled and subjected to high‐throughput 16S rRNA gene amplicon sequencing to analyze the microbial community profile therein. Despite the typically substantial heterogeneity across the community profiles of samples spatially distributed (10 m to 1.9 km apart) throughout the site, samples from two wells at the EGS testbed showed highly similar microbial community composition, suggesting the two wells intersected the same natural fracture. This evidence of natural connectivity between the two wells at the EGS testbed was later corroborated by core log analysis and sewer camera surveys into the boreholes. Besides the field case described in this study, microbial community composition as a reservoir diagnostic tool would be applicable in a much broader context such as unconventional hydrocarbon exploitation, groundwater reservoir characterization, and environmental remediation, adding valuable “hard” data capable of pinpointing the origins of fluids unambiguously.

Zhang, Yuran↗

DOE OSTI2022

Discovering hidden geothermal signatures using non-negative matrix factorization with customized k-means clustering

Discovery of hidden geothermal resources is challenging. It requires the mining of large datasets with diverse data attributes representing subsurface hydrogeological and geothermal conditions. The commonly used play fairway analysis approach typically incorporates subject-matter expertise to analyze regional data to estimate geothermal characteristics and favorability. We demonstrate an alternative approach based on machine learning (ML) to process a geothermal dataset from southwest New Mexico (SWNM). The study region includes low- and medium-temperature hydrothermal systems. Several of these systems are not well characterized because of insufficient existing data and limited past explorative work. This study discovers hidden patterns and relations in the SWNM geothermal dataset to improve our understanding of the regional hydrothermal conditions and energy-production favorability. This understanding is obtained by applying an unsupervised ML algorithm based on non-negative matrix factorization coupled with customized k-means clustering (NMFk). NMFk can automatically identify (1) hidden signatures characterizing analyzed datasets, (2) the optimal number of these signatures, (3) the dominant data attributes associated with each signature, and (4) the spatial distribution of the extracted signatures. Here, in this study, NMFk is applied to analyze 18 geological, geophysical, hydrogeological, and geothermal attributes at 44 locations in SWNM. Using NMFk, we find data patterns and identify the spatial associations of hydrothermal signatures within two physiographic provinces (Colorado Plateau and Basin and Range) and two sub-regions of these provinces (the Mogollon-Datil volcanic field and the Rio Grande rift) in SWNM. The ML algorithm extracted five hydrothermal signatures in the SWNM datasets that differentiate between low (<90°C) and medium (90-150°C)-temperature hydrothermal systems. The algorithm also suggests that the Rio Grande rift and northern Mogollon-Datil volcanic field are the most favorable regions for future geothermal resource discovery. NMFk also identified critical attributes to identify medium-temperature hydrothermal systems in the study area. The resulting NMFk model can be applied to predict geothermal conditions and their uncertainties at new SWNM locations based on limited data from unexplored regions. The code to execute the performed analyses as well as the corresponding data can be found at https://github.com/SmartTensors/GeoThermalCloud.jl.

15 GEOTHERMAL ENERGY↗

DOE OSTI2025

Surface Energy Balance Responses to Radiative Forcing in the Central Arctic From MOSAiC and Models

The Arctic surface energy budget (SEB) couples the atmosphere with the sea ice, making it useful for both studying surface processes as well as evaluating models. Improved understanding of atmosphere-ice interactions is required to improve models, requiring year-round observations to address seasonally dependent biases. This work uses novel observations from the MOSAiC expedition to quantify the responses of surface fluxes to radiative forcing over sea ice throughout a complete annual cycle. We identify two primary regimes of flux response: an ice growth regime in winter and an ice melt regime in summer. In the growth regime, changes in radiative forcing impact upwelling longwave, sensible heat, and subsurface heat fluxes, whereas in the melt regime changes in radiative forcing primarily alter the amount of melt and subsurface transmission because the surface temperature is fixed. These observed responses of surface fluxes to radiative forcing are used to evaluate seven weather forecast models during the ice growth regime. In most models, the responses of surface fluxes to radiative forcing do not match observations. Many models also have biased downwelling longwave. One model (the Coupled Arctic Forecast System; CAFS) adequately captures both the mean radiative forcing and the flux responses in winter. CAFS is further evaluated against observations spanning the full MOSAiC year, demonstrating sufficient agreement to provide a more generalized understanding of these SEB process relationships across the Arctic.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Geochemistry and Strontium Isotopes for Coal Creek Watershed, Colorado, 2021-2022

The geochemistry and strontium isotope data for Coal Creek Watershed, Colorado, consists of cation, anion, and 87Sr/87Sr isotope values from samples collected at 8 stream location along Coal Creek, samples from two groundwater springs within the watershed, and a shallow subsurface piezometer. All stream and spring samples were collected between June and October, 2021, and the shallow, near stream piezometer sample was collected in July of 2022. These data were collected to evaluate how groundwater contributions to Coal Creek originating from shallow vs deep flow paths respond seasonal drying. Understanding of groundwater-surface water interactions in montane systems in critical for the future of water availability in the Western US as groundwater contributions are expected to become more important for sustaining summer stream flows. This data package contains: (1) a csv of all cation samples; (2) a csv of all anion samples; (3) a csv of all 87Sr/87Sr isotope samples; and (4) a csv of locations for each sampling site. The dataset additionally includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Potential Challenges of HyBlend Storage in a Methane Reservoir Located in Southwestern United States

Hydrogen has been identified as a flexible energy carrier with zero or negative emissions across multiple energy systems. It is possible to utilize hydrogen by storing Hyblend, or hydrogen gas blended with methane, in existing natural gas infrastructure. However, the compatibility of adapting the current CH4 storage strategies to include H2 injection has not been fully demonstrated. It is essential that we understand the impact of H2 gas on the naturally occurring microbial community of subsurface storage reservoirs before deploying large-scale H2-CH4 storage. We designed a series of experiments that allowed us to identify potential geochemical and microbial challenges of HyBlend Storage in existing methane reservoirs. First, we collected and characterized field fluid samples from a methane reservoir located in southwestern United States. Next, we used these field fluid samples to complete a series of short-term reactor experiments at reservoir conditions (80 °C and ~1,000 psi) for a natural gas (100% CH4) and HyBlend(80% CH4/20% H2) storage environment to measure the transformation of gas content. We conducted both biotic and abiotic (sterilized) measurements to accurately understand and decouple abiotic and microbially driven processes. Overall, we found that our field sample was characterized by a diverse microbial community with the metabolic capacity for sulfur reduction, iron reduction, and acetogenesis. Across our reactors, there was minimal change in geochemistry.

hydrogen storage↗

DOE OSTI2024

Influence of Antarctic and Greenland Continental Shelf Circulation on High‐Latitude Oceans in E3SM

The science objectives of this project are to simulate and understand the impacts of both deep-basin warm-water intrusions and land-ice melt on the continental shelf circulations and sea-ice distributions around the margins of Greenland and Antarctica. As well, the role of subsurface ocean heat from the Atlantic on declining sea-ice cover in the Arctic is explored. Mesoscale processes and fine bathymetry are implicated in cross-shelf property transports around both Greenland and Antarctica. Therefore, we configured and ran an atmospheric reanalysis-forced global ocean/sea-ice simulation on a grid that reduces from 8 km at the Equator to 2 km at the poles (UH8to2) with 60 vertical levels. It was produced using the Energy Exascale Earth System Model ‘‘HiLAT’’ code (E3SMv0-HiLAT) that uses the Parallel Ocean Program (POP) and CICE5 as its ocean and sea-ice components, respectively. Two main UH8to2 simulations were carried out: one for 1975-2009 and the other for July 2016-2020 after it was initialized from a 1/25° data-assimilative ocean/sea-ice prediction system ocean/sea-ice state. The UH8to2 is not coupled to an active land-ice model. Rather, land-ice melt is represented by observationally informed freshwater fluxes (FWFs). Short (multi-year) UH8to2 simulations were conducted to understand sensitivities when Greenland ice sheet (GrIS) melt is released only at the ocean surface or when it is distributed over the upper water column in accordance with fjord melt plume behavior; these cases were compared with a no GrIS melt case. West Greenland continental shelf currents were fastest in the vertical distribution case and an increase in baroclinic conversion at the shelf break associated with increased eddy kinetic energy was found relative to the surface release case. Further, salinity is lower and meltwater volume greater in the eastern Labrador Sea in the vertical distribution case. For the Arctic, the veracity of the UH8to2 was evaluated for 2017-2020 using available observations. Simulated seasonal sea-ice thickness and concentration are realistic, but the ice is unrealistically thin in the central and eastern Arctic in the fall. Comparisons of vertical sections of ocean temperature, salinity, and buoyancy collected from Ice-Tethered Profilers (ITPs) in the eastern Arctic in the fall and winter of 2019/2020 and co-located/concurrent UH8to2 fields show the stratification over the top 100 m of the water column is too low in the model, the simulated mixed layer too deep, and the simulated subsurface Atlantic Water (AW) too warm; these biases may contribute to the sea-ice biases. A model intercomparison study using the UH8to2 and a forced 1/25° regional Arctic ocean/sea-ice (uses the HYbrid Coordinate Ocean Model and CICE5) simulation further investigates the relationship between AW and sea-ice in the eastern Arctic. The models show a mesoscale-rich pulse of Atlantic Water extending into the eastern basin that reaches maximum intensity in late winter of 2018, after which it decreases in strength. Concurrent and co-located sea-ice melt or the inhibition of sea-ice growth is seen and is attributed to halocline mesoscale eddies doming into the mixed layer with convection bringing this heat into the vicinity of the sea-ice.

58 GEOSCIENCES↗

DOE OSTI2025

Depth-dependent links between microbial taxa and nitrous oxide emissions in a long-term cotton cropping system employing soil health practices

Long-term management practices can shape soil microbial communities in ways that influence nitrogen (N) dynamics and nitrous oxide (N 2 O) emissions. We leverage a 41-year continuous cotton cropping experiment with contrasting tillage, cover cropping, and N fertilization regimes to investigate how these long-term strategies influence soil microbial communities and their associations with N 2 O fluxes during the cotton growing season. Using 16S rRNA gene metabarcoding, we assessed microbial composition in surface and subsurface soils and evaluated its relationship with temporal N 2 O emissions. Among the management practices, N fertilization – a known driver of N 2 O emissions – had the strongest effect on microbial community composition and was linked to a greater number of taxa correlated to N 2 O emissions, particularly in surface soils. Soil pH emerged as a key variable influencing microbial structure across depth and was negatively associated with both N 2 O emissions and microbial composition in the surface layers of fertilized soils. In total, 57 archaeal/bacterial taxa were correlated with N 2 O fluxes, but only seven were shared across depths, suggesting distinct microbial contributors in surface and subsurface soils. Several of these taxa have been previously reported to be associated with N and C cycling processes such as nitrate respiration or carbon turnover, indicating functional context to their correlation with N 2 O fluxes. Temporal shifts in the abundance of key taxa aligned with seasonal peaks in N 2 O emissions, notably in early and late August, and were most pronounced under conventional tillage, hairy vetch cover cropping, and N fertilization. While 16S-based associations cannot confirm functional gene presence or activity, these findings demonstrate that long-term fertilization and associated soil acidification are dominant drivers of microbial shifts linked to N 2 O emissions and highlight the importance of accounting for depth-specific and seasonal microbial dynamics when evaluating management impacts on greenhouse gas emissions.

16S rRNA gene sequencing↗

DOE OSTI2024

Linear-frequency conversion with time-varying metasurfaces

Frequency conversion is a hallmark of nonlinearity. The spectral manifestations, emergent within a system, can typically be attributed to a marked nonlinearity within the material properties, complex geometric configurations, and/or the unique functional form of interactions taking place in the constitutive subsystems. These phenomena, irrespective of their origins, have been harnessed and exploited in applications ranging from the generation of entangled photons, a cornerstone in quantum technologies, to nanomechanical frequency mixing, advancing subsurface scanning probe microscopy. Here, we propose a frequency conversion mechanism based on time-varying metasurfaces, an emerging frontier in metamaterial research. We show how temporal properties of metasurfaces can effectively emulate a nonlinear medium, thereby facilitating frequency conversion. The proposed material configuration has the potential not only to advance integrated photonics and quantum optics, but also to create opportunities in quantum sensing, quantum materials, and crucially quantum communications. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

How initial conditions-, structural-, and parameter-based model uncertainty interact and influence predictions in permafrost ecosystems: Modeling Archive

This dataset contains model output and input data, as well as source code examples for the Terrestrial Ecosystem Model with the Dynamic Vegetation Model and Dynamic Organic Soil (DVM-DOS-TEM) for the field sites Imnavait creek and the Bonanza creek Long Term Ecological Research Network (LTER). The data covers simulations from the last glacial maximum (LGM) until 2100 for a selection of paleo scenarios, setting the mean temperature of the LGM up to 10°C lower than pre-industrial conditions. The model structure was modulated to represent various model versions, and this dataset contains the relevant changes in the source code. The raw output data, the processed statistical data, the setup and processing scripts as well as parameter value distribution files from a parameter sensitivity analysis are included as well. Model outputs include active layer depth, organic soil carbon, soil layer depths, gross primary productivity (GPP) with and without nitrogen limitation, net primary productivity (NPP), soil liquid water content, heterotrophic, maintenance, and growth respiration, soil temperature, and vegetation carbon (*.nc files). The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research.Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Characterizing Reactive Transport Behavior in a Three-Dimensional Discrete Fracture Network

While several studies have linked network and in-fracture scale properties to conservative transport behavior in subsurface fractured media, studies on reactive transport cases remain relatively underdeveloped. In this study, we explore the behavior of an irreversible kinetic reaction during the interaction of two solute plumes, one consisting of species A and the other species B. When the plumes converge, these species react kinetically to form a new species C via A+B→kC. This reactive system is studied using a three-dimensional discrete fracture network (DFN) model coupled with reactive Lagrangian particle tracking. We find that the interplay of network topology and chemical properties of the reactive solutes controls reactive transport processes. The network topology drives species A and B together, and the chemical properties dictate whether and how quickly a reaction occurs. Results demonstrate that reactions are most likely to occur in high-velocity fractures that make up the network backbone. The interplay between species’ chemical properties and transport is characterized by a non-dimensional Damköhler (Da) number. We show that the spatial distribution of reactions is sensitive to Da, which subsequently influences late-time tailing behavior in outlet breakthrough time distributions. The results of this study provide initial insights into how an irreversible reaction occurs during transport in a fracture network, using a methodology that can be applied to study reactive transport in a wide range of fractured media environments and contexts.

58 GEOSCIENCES↗

DOE OSTI2022

Dissolved pyrogenic carbon leaching in soil: Effects of soil depth and pyrolysis temperature

We report Pyrogenic organic carbon (PyC) is found in soils as a heterogeneous mixture of thermally altered plant residues that range in their susceptibility to losses via mineralization, leaching, and erosion. Leaching of PyC as DOC (DPyC) within the soil profile is likely influenced by the chemical composition of solid PyC and the dominant soil processes and properties at a particular depth. Here we report the results of a 2-year laboratory decomposition and leaching study designed to investigate the interactive effects of pyrolysis temperature levels (no pyrolysis, 300 °C, and 450 °C) and soil depth (surface soil, 0–10 cm depth; subsurface soil, 50–70 cm depth) on losses of 13 C-labeled jack pine (Pinus banksiana) wood, wood pyrolyzed at 300 °C (PyC300) and wood pyrolyzed at 450 °C (PyC450) as dissolved organic carbon (DOC). Losses of wood and PyC in the 13 C- DOC pool were measured in leachates drained from soils once a month over the course of 1 year, and at one single leaching event after 2 years. We found that pyrolysis temperature levels interacted with time and soil depth to affect losses of wood and PyC as DOC throughout the 1-year incubation and leaching study, with greater DOC losses from wood in surface soil (0.73 ± 0.076% of added C) than in subsurface soil (0.40 ± 0.063% of added C) averaged across sampling time. Monthly DOC losses from PyC were not affected by soil depth. Cumulative data indicated a small contribution of wood (2.86 ± 0.07%), PyC300 (0.40 ± 0.04%), and PyC450 (0.16 ± 0.01%) to total DOC leached from soils. DOC losses from wood and PyC300 (as a proportion of added C) were greater in surface soil than in subsurface soil, whereas DOC losses from PyC450 were unaffected by differences between surface and subsurface soil. Losses of DOC were greater from PyC300 than from PyC450 in surface soil, with no significant differences in DOC losses between PyC300 and PyC450 in subsurface soil. One single leaching event at the end of the 2-year decomposition study resulted in higher DOC losses from wood in subsurface soil than in surface soil likely due to desorption, and no differences in DOC losses between PyC300 and PyC450 regardless of soil depth. Our results suggest strong interactions between the initial physicochemical composition of organic C inputs and soil properties (soil depth as a proxy) that control the mobility and transport of PyC and should be better represented in Earth system models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Co-Design of Charge Transport Superhighways to Connect Catalytic Sites in Soft Photoelectrochemical Systems

Efficient photon-to-electron-to-molecule conversion requires multi-length scale control over charge transport pathways, where electronic charges are delivered to catalytic sites under high mass transport flux. A fundamental question is how can we co-design charge transport pathways to promote efficient charge transfer to/from catalytic sites in complex three-dimensional architectures? Soft conducting polymer systems offer exceptional promise to provide three-dimensional charge transport networks, where electrolyte (ion and solvent) can interdiffuse to promote long-lived charge carriers and the molecular nature allows for strategic synthetic design of catalytic sites. Herein we combine theoretical and experimental approaches to investigate the earliest stages of photoelectrochemical deposition of near-surface catalytic sites (Pt) on soft bulk heterojunction polymeric semiconductors composed of a prototype donor (PTB7-Th) and a prototype acceptor (N2200) as a model system towards better understanding molecular catalyst-polymer site interactions. We focus initially on photoelectrochemical deposition of low Pt loadings, nanoparticle sizes (formed by progressive nucleation) below 20 nm, for both density functional theory (DFT) modeling studies and for spectroscopic characterization using surface-sensitive X-ray and UV-photoemission (XPS/UPS). DFT modeling of “n-type” N2200 slabs reveal for the first time that sulfur atoms in the thiophene units serve as the lowest-energy adsorption sites for single Pt atoms, while larger Pt clusters engage more complexly with both thiophene and naphthalene diimide (NDI) core sites. Changes in chemical composition observed by X-ray photoelectron spectroscopy (XPS) support the DFT predictions, and the angle-resolved measurements reveal that Pt nucleation initiates at subsurface sites which appear to be localized active domains that promote charge transport/transfer and enable vertical growth toward the surface. These results suggest that light-activated Pt nanoparticle deposition decorates energetically distinct sites, where photoactivity is dictated by the local energetics of those sites, and the fact that they represent the termini of charge transport “super-highways” – a small percentage of the total volume of the donor/acceptor polymeric active layer which carries most of the photocurrent generated during both Pt deposition and photoelectrochemical HER. We posit that these initial studies provide a foundational strategy for design of catalytic sites in the near surface regions of complex polymeric materials and advancing soft semiconductor-based photoelectrochemical systems. Achieving a nanometer-scale understanding of catalyst deposition and the impact of local composition and energetics on that placement, should ultimately provide the design guidelines (co-design) for a broad array of catalysts at sites that optimize that efficiency and maximize platform durability.

14 SOLAR ENERGY↗

DOE OSTI2020

Coupled Time-Lapse Full-Waveform Inversion for Subsurface Flow Problems Using Intrusive Automatic Differentiation

We describe a novel framework for estimating subsurface properties, such as rock permeability and porosity, from time-lapse observed seismic data by coupling full-waveform inversion (FWI), subsurface flow processes, and rock physics models. For the inverse modeling, we handle the back propagation of gradients by an intrusive automatic differentiation strategy that offers three levels of user control: (1) At the wave physics level, we adopted the discrete adjoint method in order to use our existing high-performance FWI code; (2) at the rock physics level, we used built-in automatic differentiation operators from the TensorFlow backend; (3) at the flow physics level, we implemented customized partial differential equation (PDE) operators for the multiphase flow equations. The three-level coupled inversion strategy strikes a good balance between computational efficiency and programming efforts, and when the gradients are chained together, it constitutes a coupled inverse system. Our numerical experiments demonstrate that the three-level coupled inverse problem is superior in terms of accuracy to a traditional decoupled inversion strategy. Additionally, our method is able to simultaneously invert for parameters in empirical relationships such as the rock physics models. Our proposed inverted model can be used for reservoir performance prediction and reservoir management/optimization purposes.

54 ENVIRONMENTAL SCIENCES↗