Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

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

At least 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Modeling the Effects of Artificial Drainage on Agriculture-dominated Watersheds using a Fully Distributed Integrated Hydrology Model: Datasets, scripts, model files

This model-data archive supports the research paper that demonstrates the integration of agricultural drainage features—specifically, narrow engineered ditches and tile drains—into a fully distributed, basin-scale integrated surface-subsurface hydrology model (ISSHM), Amanzi-ATS. The model employs innovative computational meshes aligned with agricultural ditches and incorporates the physically based Hooghoudt's drainage equation to simulate tile drainage, offering a novel strategy that enhances the accuracy of hydrological simulations.The archived dataset includes input parameters, model configurations, and select simulation outputs for the Amanzi-ATS model that successfully captured the streamflow patterns in the Portage River Watershed as validated by USGS gauge readings. Jupyter notebook for the preparation of model inputs and post-processing of outputs are also included. The model's predictive performance achieved a normalized Kling-Gupta Efficiency (KGE) of 0.81, surpassing SWAT without the necessity for site-specific calibration.The Amanzi-ATS model presented in this modeL-data archive allows for numerical experiments to explore the shifts in the flow structure under different drainage scenarios. As a tool for advancing the understanding of distributed hydrological responses and nutrient cycling, this archived model provides valuable insights for researchers, modelers, and decision-makers involved in watershed management and environmental modeling.The Watershed Workflow package is implemented in Python3. The Jupyter notebooks can be executed through multiple open-source tools, for example, Anaconda Jupyter Lab, VS Studio Code, etc. Other data files include CSV and HDF5 files, which can be read through Python scripts. The input files for the ATS model, open-source integrated hydrology, and transport model, are in XML format and can be edited in any commonly used text editors.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Phyllosilicate Formation on Early Mars Via Open-System Acid Alteration of Basaltic Glass

Smectites are widespread on Mars, but neutral/alkaline conditions favorable for smectite formation have surprisingly not led to abundant carbonates. Smectite formation on Mars could also occur in acidic environments unfavorable for carbonate formation. Acidic smectite formation has been demonstrated in batch (closed hydrologic) systems, however, the mechanisms and octahedral composition of smectite forming in acidic flow-through (open hydrologic) systems are not fully understood. Stapafell basaltic glass was hydrothermally altered (190◦C) at two flow rates corresponding to low and high water to rock ratio (W/R) and initial pH (pH_0) values of 2, 3, 4 and 6, and a batch low W/R experiment was conducted at pH_0 2. Kaolinite, montmorillonite and chlorite formed at pH_0 2 at low W/R; no phyllosilicates formed at pH_0 2 at high W/R; and lizardite formed at pH_0 ≥3 at both W/R ratios. Lizardite, kaolinite, and montmorillonite in these experiments formed by precipitation from solution and chlorite likely formed through alteration of montmorillonite and/or basalt. Saponite formed at pH_0 2 in batch conditions by alteration of basaltic glass. Comparison of experimental data with martian phyllosilicate assemblages indicated that smectite formation on Mars likely occurred under water-limited conditions. Al-rich smectite could form in open-system low W/R subsurface environments under a narrow range of pH (pH <3) while saponite could form in closed low W/R systems under acidic to alkaline conditions. Combined open and closed hydrological regimes could be responsible for development of clay mineral stratigraphies observed on Mars. The acidic conditions required for formation of Al-rich smectite were unfavorable for carbonate precipitation, but carbonate precipitation could occur together with Fe/Mg-smectite in closed systems at pH >4. The lack of carbonates occurring together with Fe/Mg-smectite could be cause by low pCO2 in subsurface closed environments and/or by overall lack of carbonate precipitation despite significant oversaturation.

S J Ralston↗
DOE OSTI2021

2D-to-3D image translation of complex nanoporous volumes using generative networks

Image-based characterization offers a powerful approach to studying geological porous media at the nanoscale and images are critical to understanding reactive transport mechanisms in reservoirs relevant to energy and sustainability technologies such as carbon sequestration, subsurface hydrogen storage, and natural gas recovery. Nanoimaging presents a trade off, however, between higher-contrast sample-destructive and lower-contrast sample-preserving imaging modalities. Furthermore, high-contrast imaging modalities often acquire only 2D images, while 3D volumes are needed to characterize fully a source rock sample. In this work, we present deep learning image translation models to predict high-contrast focused ion beam-scanning electron microscopy (FIB-SEM) image volumes from transmission X-ray microscopy (TXM) images when only 2D paired training data is available. We introduce a regularization method for improving 3D volume generation from 2D-to-2D deep learning image models and apply this approach to translate 3D TXM volumes to FIB-SEM fidelity. We then segment a predicted FIB-SEM volume into a flow simulation domain and calculate the sample apparent permeability using a lattice Boltzmann method (LBM) technique. Results show that our image translation approach produces simulation domains suitable for flow visualization and allows for accurate characterization of petrophysical properties from non-destructive imaging data.

58 GEOSCIENCES↗
NASA NTRS2013

The Origin of Carbon-bearing Volatiles in Surprise Valley Hot Springs in the Great Basin: Carbon Isotope aud Water Chemistry Characterizations

There are numerous hydrothermal fields within the Great Basin of North America, some of which have been exploited for geothermal resources. With methane and other carbon-bearing compounds being observed, in some cases with high concentrations, however, their origins and formation conditions remain unknown. Thus, studying hydrothermal springs in this area provides us an opportunity to expand our knowledge of subsurface (bio)chemical processes that generate organic compounds in hydrothermal systems, and aid in future development and exploration of potential energy resources as well. While isotope measurement has long been used for recognition of their origins, there are several secondary processes that may generate variations in isotopic compositions: oxidation, re-equilibration of methane and other alkanes with CO2, mixing with compounds of other sources, etc. Therefore, in addition to isotopic analysis, other evidence, including water chemistry and rock compositions, are necessary to identify volatile compounds of different sources. Surprise Valley Hot Springs (SVHS, 41º32'N, 120º5'W), located in a typical basin and range province valley in northeastern California, is a terrestrial hydrothermal spring system of the Great Basin. Previous geophysical studies indicated the presence of clay-rich volcanic and sedimentary rocks of Tertiary age beneath the lava flows in late Tertiary and Quaternary. Water and gas samples were collected for a variety of chemical and isotope composition analyses, including in-situ pH, alkalinity, conductivity, oxidation reduction potential (ORP), major and trace elements, and C and H isotope measurements. Fluids issuing from SVHS can be classified as Na-(Cl)-SO4 type, with the major cation and anion being Na+ and SO4 2-, respectively. Thermodynamic calculation using ORP and major element data indicated that sulfate is the most dominant sulfur species, which is consistent with anion analysis results. Aquifer temperatures at depth estimated by both dissolved SiO2 and Na-K-Ca geothermometers are in the range of 125.0 to 135.4 oC, and higher than the values measured at orifices (77.3 to 90.0 oC). CO2 and homologs of straight chain alkanes (C1-C5) were identified in gas samples. Carbon isotope values of alkanes increase with carbon numbers. The 13C fractionation between CO2 and dissolved inorganic carbon suggests they are out of carbon isotope equilibrium. The hypothesis regarding the formation of carbon-bearing compounds in SVHS may involve two processes: 1) Under high heat flow conditions which are caused by regional faulting and crustal extension, original high molecular weight organic compounds (kerogens) in clay-rich rocks decomposed to generate methane and other alkane homologs. 2) The SVHS area is associated with outflow structures, and distant from the heat source. Anaerobic oxidation of methane (AOM) with sulfate at shallow depth (< 90 oC) is suggested as being responsible for the generation of CO2 in SVHS.

Fu, Qi↗
DOE OSTI2022

Pore-Scale Microenvironments Control Anthropogenic Carbon Mineralization Outcomes in Basalt

Thin sections and hand samples from 50 sidewall cores from the Wallula Basalt Pilot Demonstration, a basaltic carbon sequestration demonstration, provided the opportunity for the in-depth analysis of carbon mineralization induced by the injection of supercritical CO 2 . In this study, we used optical petrography and scanning electron microscopy to characterize the physical and chemical characteristics of the basalt components influenced by carbon mineralization reactions from all available hand samples and thin sections within the three CO 2 injection zones and caprock flow interiors. We found extensive carbonate mineralization, mostly in the form of nodules that were shown to be chemically zoned: Ca-dominant in the core regions and Ca-bearing Fe-dominant in the outer regions. Carbonate mineralization also took the form of fracture-filling carbonate cement, and acicular aragonite was also observed. Overall, we clarified the structural and paragenetic relationships between newly formed minerals, identifying a new fibro-palagonite-like, poorly crystalline silicate phase that grew on the carbonate nodules and pore-lining zeolites. Here we observed Fe-dominant carbonate precipitates surrounding acicular aragonite and rhombohedral Ca-carbonate cores, whereas previous studies of these zoned nodules did not observe these structures. A comprehensive accounting of the carbon mineralization products is vital to understand and predict the behavior of supercritical CO 2 in the subsurface given both the diversity of the host rock between and within injection zones, especially considering that the morphology and chemistry of the diverse precipitates are influenced by the pore-scale microenvironments of the basalt.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2022

Europa Clipper Mission: System Integration Review Report

Jupiter’s icy moon Europa is a prime target in our exploration of potentially habitable worlds beyond Earth. The combination of a subsurface liquid water layer in contact with a rocky seafloor may yield an ocean rich in the elements and energy needed for the emergence of life, and for potentially sustaining life through time. Europa may hold the clues to one of NASA’s long-standing quests – to determine whether or not we are alone in the universe. The Europa Clipper mission will characterize Europa’s habitability as the first step in the search for potential life at Europa by conducting approximately four dozen flybys. The project is preparing for a System Integration Review in November of 2021. This paper will summarize changes from the mission’s Preliminary Design Review (PDR) baseline, including science objectives, launch vehicle, interplanetary trajectory, science tour, flight system, payload and mission phases. Since PDR several of the Level 1 science requirements were descoped and the science requirement flow was realigned, in order to control cost growth. The realignment deemphasized science redundancy, which had resulted in an over-constrained science tour, in favor of a more targeted linkage to the primary responding instrument. Additionally, changes to the PDR flight system design baseline include changes to the sounding radar antenna design, the magnetometer architecture, and the heat redistribution system pumps. All changes have been incorporated into a reference science tour used to demonstrate that the requirements and design of the flight system and payload are mature enough to support flight builds.

Cook, Kendra↗
NASA NTRS2024

Dielectric Permittivity Measurements of Geologic Samples at Venus Surface Temperatures

Synthetic-aperture radar (SAR) and radar sounding are two powerful remote sensing techniques used to examine the surface and subsurface of planets. Their value is critical in the case of Venus due to this planet’s global and visually opaque cloud layer. The recently selected VERITAS and EnVision orbital missions, slated to study the planet in the early 2030’s, employ radar instruments to image the surface (VERITAS and EnVision) and probe the subsurface (EnVision). For both techniques, the complex dielectric permittivity of the surface and subsurface materials control radar backscatter intensity and penetration depth. Since the 1970’s, many studies have characterized the dielectric permittivity of geologic materials collected both on Earth and on the Moon to constrain the radar reflectivity and transparency at a variety of radar frequencies [e.g., 1]. Two consistent findings from those studies are that complex permittivity depends on both the density (or porosity) and the composition of the material. Denser materials and more metal-rich compositions tend to have higher permittivity values and to incur greater propagation path loss. Water content also elevates permittivity of otherwise dry rocks/soils/regolith [e.g., 2]. However, permittivity is independent of frequency for the 0.5-18 GHz range, within which SAR imaging commonly lies. The one aspect poorly explored thus far is the effect of elevated temperatures on dielectric permittivity. Previous work suggests that increased lattice vibrations and flow of free electrons associated with an increase in temperature elevate both the real (reflectivity) and imaginary (path loss) parts of the complex permittivity. We are currently conducting a 3-year campaign to measure the complex dielectric permittivity of Venusrelevant geologic materials at Venus surface temperatures (380°C - 460°C) and with frequencies relevant to VERITAS, EnVision, and Magellan imaging synthetic aperture radars.

M Barmatz↗
DOE OSTI2024

Sub-Diurnal Methane Variations on Mars Driven by Barometric Pumping and Planetary Boundary Layer Evolution

In recent years, the Tunable Laser Spectrometer within the Sample Analysis at Mars (TLS-SAM) instrument on board the Mars Science Laboratory (MSL) Curiosity rover has detected methane variations in the atmosphere at Gale crater. Methane concentrations appear to fluctuate seasonally as well as sub-diurnally, which is difficult to reconcile with an as-yet-unknown transport mechanism delivering the gas from underground to the atmosphere. To potentially explain the fluctuations, we consider barometrically induced transport of methane from an underground source to the surface, modulated by temperature-dependent adsorption. The subsurface fractured-rock seepage model is coupled to a simplified 1-D atmospheric mixing model to provide insights on the pattern of atmospheric methane concentrations in response to transient surface methane emissions, as well as to predict sub-diurnal variation in methane abundance for the northern summer period, which is a candidate time frame for a MSL Curiosity sampling campaign. Our analysis suggests that there is a lower limit to the subsurface fracture density that can produce the observed methane patterns, below which the atmospheric methane variations would be out of phase with the observations. The best-performing model scenarios indicate a significant, short-lived methane pulse just prior to sunrise, the detection of which by TLS-SAM would be a potential indicator of the contribution of barometric pumping to Mars' atmospheric methane variations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Detecting Organic Compounds Released from Iron Oxidizing Bacteria using Sample Analysis at Mars (SAM) Like Instrument Protocols

Mars is a planet of great interest for Astrobiology since its past environmental conditions are thought to have been favourable for the emergence life. At present, the Red Planet is extremely cold and dry and the surface is exposed to intense UV and ionizing radiation, conditions generally considered to be incompatible with life as we know it on Earth. It was proposed that the shallow subsurface of Mars, where temperatures can be above freezing and liquid water can exist on rock surfaces, could harbor chemolithoautotrophic bacteria such as the iron oxidizing microorganism Pseudomonas sp. HerB. The Mars Science Laboratory (MSL) mission will provide the next opportunity to carry out in situ measurements for organic compounds of possible biological origin on Mars. One instrument onboard MSL, called the Sample Analysis at Mars (SAM) instrument suite, will carry out a broad and sensitive search for organic compounds in surface samples using either high temperature pyrolysis or chemical extraction followed by gas chromatography mass spectrometry. We present gas chromatograph mass spectrometer (GC/MS) data on crushed olivine rock powders that have been inoculated with Pseudomonas sp. HerB at different concentrations ranging from approx 10(exp 2) to 10(exp 7) cells per gram. The inoculated olivine samples were heated under helium carrier gas flow at 500 C and the pyrolysis products concentrated using a SAM-like hydrocarbon trap set at -20 C followed by trap heating and analysis by GC/Ms. In addition, the samples were also extracted using a low temperature "one-pot" chemical extraction technique using N-methyl, N-(tert-butyldimethylsilyl) trifluoroacetamide (MTBSTFA) as the silylating agent prior to GC/MS analysis. We identified several aldehydes, thiols, and alkene nitriles after pyrolysis GC/MS analysis of the bacteria that were not found in the olivine control samples that had not been inoculated with bacteria. The distribution of pyrolysis products extracted from the bacteria was clearly distinct from similar GC/MS analyses of the carbonaceous meteorite Murchison that was dominated by sulfur containing aromatic compounds. A similar comparison, if organic compounds are detected by SAM on Mars, could be useful to help discriminate between meteoritic or biological origins.

Glavin, D. P.↗
NASA NTRS2013

The Origin of Carbon-bearing Volatiles in Surprise Valley Hot Springs in the Great Basin: Carbon Isotope and Water Chemistry Characterizations

There are numerous hydrothermal fields within the Great Basin of North America, some of which have been exploited for geothermal resources. With methane and other carbon-bearing compounds being observed, in some cases with high concentrations, however, their origins and formation conditions remain unknown. Thus, studying hydrothermal springs in this area provides us an opportunity to expand our knowledge of subsurface (bio)chemical processes that generate organic compounds in hydrothermal systems, and aid in future development and exploration of potential energy resources as well. While isotope measurement has long been used for recognition of their origins, there are several secondary processes that may generate variations in isotopic compositions: oxidation, re-equilibration of methane and other alkanes with CO2, mixing with compounds of other sources, etc. Therefore, in addition to isotopic analysis, other evidence, including water chemistry and rock compositions, are necessary to identify volatile compounds of different sources. Surprise Valley Hot Springs (SVHS, 41 deg 32'N, 120 deg 5'W), located in a typical basin and range province valley in northeastern California, is a terrestrial hydrothermal spring system of the Great Basin. Previous geophysical studies indicated the presence of clay-rich volcanic and sedimentary rocks of Tertiary age beneath the lava flows in late Tertiary and Quaternary. Water and gas samples were collected for a variety of chemical and isotope composition analyses, including in-situ pH, alkalinity, conductivity, oxidation reduction potential (ORP), major and trace elements, and C and H isotope measurements. Fluids issuing from SVHS can be classified as Na-(Cl)-SO4 type, with the major cation and anion being Na+ and SO4(2-), respectively. Thermodynamic calculation using ORP and major element data indicated that sulfate is the most dominant sulfur species, which is consistent with anion analysis results. Aquifer temperatures at depth estimated by both dissolved SiO2 and Na-K-Ca geothermometers are in the range of 125.0 to 135.4 C, and higher than the values measured at orifices (77.3 to 90.0 C). CO2 and homologs of straight chain alkanes (C1-C5) were identified in gas samples. Carbon isotope values of alkanes increase with carbon numbers. The C-13 fractionation between CO2 and dissolved inorganic carbon suggests they are out of carbon isotope equilibrium. The hypothesis regarding the formation of carbon-bearing compounds in SVHS may involve two processes: 1) Under high heat flow conditions which are caused by regional faulting and crustal extension, original high molecular weight organic compounds (kerogens) in clay-rich rocks decomposed to generate methane and other alkane homologs. 2) The SVHS area is associated with outflow structures, and distant from the heat source. Anaerobic oxidation of methane (AOM) with sulfate at shallow depth (< 90 C) is suggested as being responsible for the generation of CO2 in SVHS.

Fu, Qi↗
DOE OSTI2021

Calendar Year 2020 Post-Closure Monitoring Report for Corrective Action Unit 99: Rainier Mesa/Shoshone Mountain, Underground Test Area, Nevada National Security Site, Nevada (Rev. 1)

This report presents the results of monitoring conducted for water quality, water levels, and institutional controls in calendar year (CY) 2020 by the U.S. Department of Energy (DOE), Environmental Management (EM) Nevada Program’s Underground Test Area (UGTA) Activity at Corrective Action Unit (CAU) 99, Rainier Mesa/Shoshone Mountain (RM/SM), at the Nevada National Security Site (NNSS), Nevada. This report also presents analytical laboratory results (including results of quality assurance/quality control samples, such as field duplicates [FDs]), and verification of use restrictions (URs), institutional controls, and water usage. Groundwater samples were collected; water levels were measured; and well site surveillance was conducted in support of the Underground Test Area Closure Report (CR) for Corrective Action Unit 99: Rainier Mesa/Shoshone Mountain, Nevada National Security Site, Nevada (DOE/EMNV, 2020). Rainier Mesa (RM) was the site of 61 underground nuclear tests, and Shoshone Mountain (SM) was the site of 6 underground nuclear tests. As a result of these activities, which took place from 1957 to 1992, radionuclides (RNs) were released in the subsurface in the vicinity of the detonations. These 67 underground nuclear tests are associated with 66 specific corrective action site (CAS) numbers, as listed in the Federal Facility Agreement and Consent Order (FFACO) (1996, as amended). Two tests, HURON LANDING and DIAMOND ACE (conducted simultaneously), are included within one CAS. The CR (DOE/EMNV, 2020) establishes the regulatory boundaries and regulatory boundary objectives, monitoring program, UR boundaries and URs, and other institutional controls agreed to by EM Nevada Program and the Nevada Division of Environmental Protection (NDEP) for closure of CAU 99. Regulatory boundaries for the RM/SM CAU were established to protect receptors of groundwater from RN contamination within the three downgradient groundwater basins that receive recharge from RM. The SM regulatory boundary was established to verify that RN contamination does not reach the lower carbonate aquifer below SM. The UR boundaries were established based on contaminant simulation forecasts from the flow and transport simulations with consideration of site-specific topographical controls. The URs associated with the UR boundaries were identified to protect onsite workers and the public from inadvertent exposure to contaminated groundwater as forecasted to occur assuming current conditions. The UR is intended to restrict activities that might expose workers to contaminated groundwater within the UR area. In addition to the URs, other institutional controls are established to monitor and limit access to groundwater. These include federal ownership and management in perpetuity, controlled access of the NNSS and surrounding areas, and reporting water use on the NNSS and surrounding hydrographic basins.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Editorial: Linking Hydrological and Biogeochemical Processes in Riparian Corridors

The riparian corridor is a key component of the critical zone and an essential component of watershed systems. According to Merriam-Webster Dictionary, the word riparian is derived from the Latin word riparius, meaning “existing alongside a river.” Riparian corridors typically extend from a few meters to hundreds of meters adjacent to a river and are marked by rich biodiversity, vegetation, and intense biogeochemical activity. They act as integrators of watershed processes and constitute the primary pathways for the subsurface geochemical exports from the watershed. Although riparian corridors comprise only 2–10% of a watershed's area, as much as 90–98% of biogeochemical processing in watersheds occurs in this region, thereby affecting the subsurface geochemical exports and downstream river water quality. Indeed, the riparian corridor is a good example of the Pareto principle. This outsized contribution occurs at the interface between aquatic (river) and terrestrial (land) environments, where interactions between hydrologic and biogeochemical processes are intensified. Variations in the river corridor over time can thus also have outsize impacts. Therefore, it is important to understand the hydrological and biogeochemical linkages in riparian corridors to determine water availability and quality for sustainable management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Water Observations of Flow/No-Flow for the East-Taylor Watershed, Colorado (June-July 2025 and 2026)

This dataset provides multi-year, ground-truth visual observations of surface water flow/no-flow conditions within the East-Taylor Watershed, Colorado, collected during June and July of 2025 and 2026. In June and July 2025, on-the-ground visual observations of flow/no-flow were collected as part of the Watershed Function Scientific Focus Area (SFA) and Rocky Mountain Biological Laboratory (RMBL) Colorado Headwaters Ecological Spectroscopy Study (CHESS) campaign (further details are provided within the CHESS Project Description). We obtained 377 water observations of flow/no-flow within the East-Taylor Watershed, Colorado. These ground-truth observations were collected to validate classification maps from remote sensing data and model results within the East-Taylor Watershed. In 2025, flow/no-flow measurements were collected using a field-based app for the CHESS Campaign (Zerion iForm). Within the field app, a water observation form was created to collect coordinates and metadata about the observation. Information collected for the water observation points included information about visually-assessed streamflow presence/absence (standard question obtained from Colorado State University’s StreamTracker project), flow estimate, stream or ponded area width, canopy cover, manganese films, iron seeps, and beaver activity. For 2025 water observations, this dataset contains: (1) a data file with the water observations and coordinates (2025_Water_Observations.csv); (2) a Keyhole Markup Language Zipped (KMZ) with the water observation locations and metadata (2025_Water_Observations_Locations.kmz); (3) photos (.jpg and .jpeg) of the water observation points, organized by location, contained within 2025_Water_Observations_FieldPhotographs.zip file; and (4) water observation protocols and figures (2025_Water_Observation_Protocols.pdf). In June and July 2026, on-the-ground visual observations of flow/no-flow were collected as part of the Watershed Function SFA project. We obtained 365 water observations of flow/no-flow within the East-Taylor Watershed, Colorado. The 2026 observations focused on collecting repeat measurements at the 2025 flow/no-flow observation locations conducted as part of the CHESS campaign. These ground-truth observations were collected to understand differences in flow/no-flow in 2026, given the unprecedented 2026 drought in Colorado. In 2026, flow/no-flow measurements were collected using ArcGIS (Geographic Information System) Survey123. Within the field app, a water observation form was created to collect coordinates and metadata about the observation. Information collected for the water observation points included repeat information from the 2025 water observation effort, including visually-assessed streamflow presence/absence (standard question obtained from Colorado State University’s StreamTracker project), flow estimate, stream or ponded area width, canopy cover, manganese films, iron seeps, beaver activity, and a new metadata component of estimated stream depth (for select locations). For 2026 water observations, this dataset contains: (1) a data file with the water observations and coordinates (2026_Water_Observations.csv); (2) a Keyhole Markup Language Zipped (KMZ) with the water observation locations and metadata (2026_Water_Observations_Locations.kmz); (3) photos (.jpg) of the water observation points, organized by location, contained within 2026_Water_Observations_FieldPhotographs.zip file; and (4) water observation protocols and figures (2026_Water_Observation_Protocols.pdf). For 2025 and 2026 water observations, this dataset contains: (1) a location metadata file (locations.csv); (6) a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and (7) a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. 2026-09-02: This dataset was updated to include 2026 water observation measurements. The 2025 observation files were also updated to ensure a consistent file naming convention across water observation years.

2018 NEON and 2025 CHESS Campaigns↗
DOE OSTI2020

Nanoparticle Tracking to Probe Transport in Porous Media

From the granular and fractured subsurface environment to highly engineered polymer membranes used in pharmaceutical purification, porous materials are ubiquitous in nature and industrial applications. In particular, porous media are used extensively in processes including water treatment, pharmaceutical sterilization, food/beverage processing and heterogeneous catalysis, where hindered mass transport is either essential to the process or a necessary, but undesirable limitation. Unfortunately, there are currently no universal models capable of predicting mass transport based on a description of the porous material because real porous materials are complex, and because many coupled dynamic mechanisms (e.g. adsorption, steric effects, hydrodynamic effects, electrostatic interactions, etc.) give rise to the observed macroscopic transport phenomena. While classical techniques, like nuclear magnetic resonance and dynamic light scattering, provide useful information about mass transport in porous media at the ensemble level, they provide limited insight into the microscopic mechanisms that give rise to complex phenomena such as anomalous diffusion, hindered pore-space accessibility, and unexpected retention under flow, among many others. To address this issue, we have developed refractive index matching imaging systems, combined with single-particle tracking methods, allowing the direct visualization of single particle motion within a variety of porous materials. Here, we summarize our recent efforts to advance the understanding of nanoparticle transport in porous media using single-particle tracking methods in both fundamental and applied scenarios. First, we describe the basic principles for two-dimensional and three-dimensional single-particle tracking in porous materials. Then, we describe concrete examples of nanoparticle transport in porous materials from two perspectives: (1) understanding fundamental elementary particle transport processes in porous media, including pore accessibility and cavity escape, which limit transport in porous media; (2) facilitating applications in industrial processes, e.g. by understanding the mechanisms of particle fouling and remobilization in filtration membranes. Finally, we provide an outlook of opportunities associated with investigating other types of mass transport in confined environments using single-particle tracking methods, including electrophoretic and self-propelled motion.

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

Prediction and uncertainty quantification of shale well performance using multifidelity Monte Carlo

Uncertainty quantification is an integral component of reservoir management, especially considering the inherent uncertainty in subsurface systems. While a standard practice to estimate the uncertainty, Monte Carlo (MC) simulation is computationally intense when the sampling population comprises high-fidelity simulations. Alternatively, the Multi-fidelity Monte Carlo (MFMC) simulation overcomes this computational intensity by integrating low- and high-fidelity simulations. Our goal is to minimize the number of expensive high-fidelity simulations while maintaining accuracy and using numerous fast and cheap low-fidelity simulations to efficiently sample to input parameter space of interest. We selected gas production from unconventional wells to demonstrate the potential speedups and accuracy of the MFMC approach. The model fidelity usually determines the trade-off between accuracy and efficiency. While the high-fidelity model is more accurate, the low-fidelity model is more efficient. Our high-fidelity simulation includes reservoir simulations of a hydraulically fractured well. On the other hand, our low-fidelity model comprises the parallel-plate flow model. We used differential programming to efficiently solve the 1D flow model, where automatic differentiation is used to efficiently compute the gradients. We matched the production profile of high-fidelity simulations with our low-fidelity simulations. Then, we used a support vector regression to map the high- and low-fidelity input parameters. The mapping function is essential to tune the low-dimensional parameter space of the low-fidelity model to the high-dimensional parameter space of the high-fidelity model. We found that we can use a combination of 9 high fidelity and 10,000 low fidelity simulations to efficiently and accurately simulate pressure management. This method is at least two orders of magnitude faster than only using high-fidelity simulations. Finally, from a broader perspective, MFMC could efficiently estimate the uncertainty of various systems and models, integrating low- and high-fidelity models.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2023

Data, figures, animations, and scripts associated with the manuscript "Integrated Effects of Site Hydrology and Vegetation on Exchange Fluxes and Carbon Cycling at the Coastal Terrestrial Aquatic Interface"

This package contains the data, figures, animations, and scripts used in "Integrated Effects of Site Hydrology and Vegetation on Exchange Fluxes and Carbon Cycling at the Coastal Terrestrial Aquatic Interface" (Li et al., 2023). This study examines the interactions between soil, vegetation and hydrologic conditions in coastal areas, with a focus on the Chesapeake Bay region. The researchers used a mechanistic model called ATS-PFLOTRAN to explore how these interactions affect exchange process in different regions and the carbon/nitrogen cycle across the terrestrial-aquatic interface (TAI). Simulation scenarios isolate the effects of control factors and reaction constants derived from laboratory experiments. The results show a carbon cycle "hot zone" in coastal wetlands and in the transition zone between wetlands and uplands. Transpiration enhances fluxes between the surface and subsurface domains and increases dissolved oxygen in TAI. The decomposition of leaf-derived organic carbon provides an additional source of carbon for aerobic respiration and denitrification in the TAI. Microbial activity plays a key role in controlling redox conditions and their variability. This modeling study improves the understanding of complex TAI interactions and facilitates the representation of coastal ecosystems in larger-scale Earth system models.Several files can be found from this data package.1. Readme,md: This file describe the Title, Target Journal, Target submission date, Co-author, Science Questions, Hypotheses, Key words, Key message, Model and Data, Repo Structure. The user can read this file first and then go to details.2. mesh.zip: This file contains the mesh file for simulation cases. 3. Figure.zip: This file contains figures used in the manuscript . 4. animation.zip: This file contains animations used in the manuscript . 5. Data.zip: This file contains all required input data, such as concentration and flow boundaries. Also, it contains the DEM and processes results on area fraction. 6. Simulation_setup_and_results.zip: This file contains setup of cases used in the manuscript . The result files are too large, if you need that data, please contact to the Author. 7. notebooks.zip: This file contains the Jupyter notebooks for performing sensitivity analysis and other pre-process and post-process analyses.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1997

Geosat Data Assimilation with Application to the Eastern North Atlantic

An attempt is made to determine the three-dimensional ocean circulation from satellite altimeter measurements by assimilating Geosat sea surface height data into an eddy-resolving QuasiGeostrophic (QG) model of the eastern North Atlantic Ocean. Results are tested against independent information from hydrographic field observations and moored current meter data collected during the Geosat ERM. The comparison supports the concept of inferring aspects of the three-dimensional flow field from sea surface height observations by combining altimetric measurements with the dynamics of ocean circulation models. A Holland-type QG model with open boundaries is set up on a 2000 km X 2000 km domain of the eastern North Atlantic between 25 deg. and 45 deg. N, 32 deg. and 8 deg. W. By using a simple nudging technique, about two years of Geosat altimeter data are assimilated into the model every five days as space-time objective analyses on the model grid. The error information resulting from the analysis is used during the assimilation procedure to account for data uncertainties. Results show an intense eddy field, which in the surface layer interacts with a meandering Azores Front. Compared to Geosat, the model leads to smoothed fields that follow the observations. Model simulations are significantly correlated with hydrographic data from March 1988 and June 1989, both close to the surface and in the subsurface. Good agreement is also found between the model velocity fields and moored current meter data in the top two model layers. The agreement is visually weak in the bottom layer, although a coherence analysis reveals an agreement between the model simulation and current meter data over the full water column at periods exceeding 80 days.

Stammer, Detlef↗
DOE OSTI2024

3D Quantification of Pore Networks and Anthropogenic Carbon Mineralization in Stacked Basalt Reservoirs

Basalt formations are promising candidates for the geologic storage of anthropogenic CO 2 due to their storage capacity, porosity, permeability, and reactive geochemical trapping ability. The Wallula Basalt Carbon Storage Pilot Project demonstrated that supercritical CO 2 injected into >800 m deep Columbia River Basalt Group stacked reservoir flow tops mineralizes to ankerite-siderite-aragonite on month-year timescales, with 60% of the 977 metric tonnes of CO 2 converted within two years. The potential impacts of mineral precipitation and consequent changes on rock porosity, pore structure, pore size, and pore size distributions have likely been underestimated hitherto. Herein, we address these knowledge gaps using X-ray Micro Computed Tomography (XMT) to evaluate the pore network architecture of sidewall cores recovered two years after CO 2 injection. In the present study, we performed a detailed quantitative analysis of the CO 2 -reacted basalt cores by XMT imaging. Reconstructed 3D images were analyzed to study the distribution and volumetric details of porosity and carbonate nodules in the cores along with the various other phases, providing insight into paragenesis and carbonate growth mechanisms, including mineralogic/chemical zonation. Finally, these findings are being used to parametrize multiphase reactive transport models to predict the fate and transport of subsurface CO 2 , enabling scale-up to commercial-scale geologic carbon storage in basalts and other reactive mafic-ultramafic formations.

54 ENVIRONMENTAL SCIENCES↗