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.

SourceResearch connectionFrom yearThrough year

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

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Crustal deformation, the earthquake cycle, and models of viscoelastic flow in the asthenosphere

The crustal deformation patterns associated with the earthquake cycle can depend strongly on the rheological properties of subcrustal material. Substantial deviations from the simple patterns for a uniformly elastic earth are expected when viscoelastic flow of subcrustal material is considered. The detailed description of the deformation pattern and in particular the surface displacements, displacement rates, strains, and strain rates depend on the structure and geometry of the material near the seismogenic zone. The origin of some of these differences are resolved by analyzing several different linear viscoelastic models with a common finite element computational technique. The models involve strike-slip faulting and include a thin channel asthenosphere model, a model with a varying thickness lithosphere, and a model with a viscoelastic inclusion below the brittle slip plane. The calculations reveal that the surface deformation pattern is most sensitive to the rheology of the material that lies below the slip plane in a volume whose extent is a few times the fault depth. If this material is viscoelastic, the surface deformation pattern resembles that of an elastic layer lying over a viscoelastic half-space. When the thickness or breath of the viscoelastic material is less than a few times the fault depth, then the surface deformation pattern is altered and geodetic measurements are potentially useful for studying the details of subsurface geometry and structure. Distinguishing among the various models is best accomplished by making geodetic measurements not only near the fault but out to distances equal to several times the fault depth. This is where the model differences are greatest; these differences will be most readily detected shortly after an earthquake when viscoelastic effects are most pronounced.

Cohen, S. C.↗

DOE OSTI2023

Dissolution and reprecipitation of amorphous silica in silica Rich shales induces Non-Monotonic evolution of porosity in acidic reactive environments

Advances in sustainable subsurface energy technologies are crucial for meeting our energy and resource needs for a climate-resilient future. Novel strategies to harness subsurface shale reservoirs for recovering valuable metals and for enabling CO 2 storage are influenced by the morphological and mineralogical heterogeneities of these materials. In this context, delineating the interactions of highly acidic solutions such as wet supercritical CO 2 on shales with varying mineralogy is crucial to inform the stability of caprock seal for CO 2 storage and enhancements in permeability for fluid transport, reactivity, and storage. The feedback chemical effects associated with the interactions of acidic solutions on the morphologies and mineralogies of shales have not been extensively investigated. These insights are crucial for assessing temporal changes in the reactivity and the fate of the fluids in subsurface environments. In this study, we investigate the effect of 1M HCl solution on the chemistry and morphology of three different shale samples with varying carbonate, clay and silica contents. An increase in the amorphous content, from 37% to 41.3%, of silica-rich and carbonate/clay lean shale is noted due to reactions with an acidic solution which is attributed to the dissolution of Si-bearing phases such as clays, accompanied by SiO 2 precipitation. In shales bearing high content of clays and carbonates, significant increase in the pore volumes and surface areas are noted. Non-monotonic changes in the micron-scale porosity of silica rich–carbonate/clay lean (e.g., Mowry shale) are noted using in-situ X-ray microtomography experiments. Due to the initial mobilization of silica and dissolution of carbonate/clay phases, the total porosity slightly increases from 6.7% to 10.7% followed by a decrease to similar to 4% caused by SiO 2 reprecipitation. In conclusion, these findings suggest that even though silica is less reactive in acidic environments, the changes in the amorphous and crystalline content due to dissolution and reprecipitation alter the porosity and fluid flow paths.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2022

The ParFlow Sandtank: An interactive educational tool making invisible groundwater visible

Physical aquifer models are a highly effective teaching tool for hydrology education, however they come with inherent limitations that include the high cost to purchase, the static configuration of the model materials, the time required to visualize hydrogeological phenomena, and the effort to reset and clean them over time. To address these and other limitations, we have developed an interactive computer simulation of a physical aquifer model called the ParFlow Sandtank. In this gamified interface, users run the simulation using a familiar web-app like interface with sliders and buttons while learning real hydrologic concepts. Our user interface allows participants to dive into the world of hydrology, understanding assumptions about model parameters such as hydraulic conductivity, making decisions about inputs to groundwater aquifer systems such as pumping rates, visualizing outputs such as stream flow, transport, and saturation, and exploring various factors that impact real environmental systems such as climate change. The ParFlow Sandtank has already been used in a variety of educational settings with more than 9,000 users per year, and we feel this emerging educational tool can be used broadly in educational environments and can be scaled-up to provide greater accessibility for students and educators. Here we present the capabilities and workflow of the ParFlow Sandtank, two use cases, and additional tools and custom templates that have been developed to support and enhance the reach of the ParFlow Sandtank.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRSDate not supplied

Progress in Simulated Water Well Performance on Mars

Recent studies of human Mars missions considered the impacts of an abundant supply of in-situ, accessible water on these mission scenarios. Discovery of exposed water ice scarps in Martian mid-latitudes has bolstered the evidence for massive amounts of almost pure water in buried deposits in regions considered candidates for these future human missions. This paper describes progress towards adapting a long-standing terrestrial technique for accessing and extracting water from these mid-latitude sources of ice. This approach relies on mechanical drills to access the ice through overlying debris. Once the ice layer has been reached, a technique known as a Rodriguez Well, or Rodwell, is used to melt the ice, store the resulting water in a subsurface cavity until needed, and then pump the water to the surface for use. Previous work by the authors utilized a computer simulation to predict the performance of a Martian Rodwell. This simulation was originally developed to predict performance in terrestrial Polar Regions. Whereas the basic approach is appropriate for a similar well on Mars, several parameters had been empirically derived and required experiments simulating the Martian environment to determine the values appropriate for a Martian Rodwell simulation. These experiments have now been completed and preliminary values have been determined for the empirical parameters. Test results are consistent with a dimensionless Sh(Ra) correlation developed by Bower and Saylor (2009) and based on lab tests of water evaporating from a pool into a large air chamber. They are also consistent with correlations derived by Ingersoll (1970) and by Hecht (2002), based on natural-convection heat transfer. Each of these earlier correlations found Sh~Ra1/3, with the 1/3rd power characteristic of large Ra where the boundary layer over the pool surface is turbulent. A somewhat surprising result was that this power-law was found to extend into a range of Ra where laminar flow would be expected. Another surprising result was found in that instead of following established correlations for natural-convection heat transfer over flat plates, the dimensionless heat-transfer rates (Nu) were much larger and approximately independent of Ra. This suggests that mass transfer from the pool enhances convective heat transfer, possibly by enhanced mixing in the boundary layer. Although this possibility has been suggested in the literature, if it is proven to be true under these conditions it could play an important role in the use of a Rodwell on Mars, where evaporative mass transfer is enhanced relative to terrestrial wells.

Human Mars Mission↗

DOE OSTI2021

Entropic Control of HD Exchange Rates over Dilute Pd-in-Au Alloy Nanoparticle Catalysts

Dilute Pd-in-Au alloy catalysts are promising materials for selective hydrogenation catalysis. Extensive surface science studies have contributed mechanistic insight on the energetic aspect of hydrogen dissociation, migration and recombination on dilute alloy systems. Yet, translating these fundamental concepts to the kinetics and free energy of hydrogen dissociation on nanoparticle catalysts operating at ambient pressures and temperatures remains challenging. Here, the effect of the Pd concentration and Pd ensemble size on the catalytic activity, apparent activation energy and rate limiting process is addressed by combining experiment and theory. Experiments in a flow reactor show that a compositional change from 4 to 8 atm% Pd of the Pd-in-Au alloy catalyst leads to strong increase in activity, even exceeding the activity per Pd atom of monometallic Pd under the same conditions, albeit with an increase in apparent activation energy. First-principles calculations show that the rate and apparent activation enthalpy for HD exchange increase when increasing the Pd ensemble size from single Pd atoms to Pd trimers in a Au surface, suggesting that the ensemble size distribution shifts from mainly single Pd atoms on the 4 atm% Pd alloy to larger Pd ensembles of at least three atoms for the 8 atm% Pd/Au catalyst. The DFT studies also indicated that the rate-controlling process is different: H 2 (D 2 ) dissociation determines the rate for single atoms whereas recombination of adsorbed H and D determines the rate on Pd trimers. Both experiment and theory suggest that the increased reaction rate with increasing Pd content and ensemble size stems from an entropic driving force. Lastly, our results support hydrogen migration between Pd sites via Au and indicate that the dilute alloy design prevents the formation of subsurface hydrogen, which is crucial in achieving high selectivity in hydrogenation catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2000

Radar Reflectivity of the Martian Polar Regions

Radar experiments provide a unique method of probing the surfaces and subsurfaces of planetary bodies. Information on surface and subsurface structures and properties can be extracted from radar data. There is a well developed history of radar investigations of the planet Mars, beginning with the first reports of variations of scattering properties as a function of martian longitude in the mid-1960's. Because of the rapid rotation of Mars, the standard technique of delay-doppler mapping cannot be used as effectively as it is on Mercury and Venus, making it a more difficult radar target. Techniques have been developed to overcome this difficulty, including an inversion technique to combine many doppler-only (CW) experiments from different viewing geometries into a map of surface reflectivity, and new random long-code techniques. These techniques still suffer from ambiguities, however. A third technique which does not suffer from these ambiguities is the combination of the powerful Goldstone transmitter with the VLA as the receiving instrument to create a combined radar imaging instrument. We have used this combined radar instrument to image the surface of Mars in 3.5-cm radar reflectivity during the 1988, 1992/93, and 1999 oppositions. During the 1988 experiments, the residual south polar ice cap (RSPIC) was the brightest radar reflector on the planet - intrinsically brighter than even the Tharsis lava flows. This was quite an unexpected result. In contrast, during the 1992/93 experiments, the residual north polar ice cap (RNPIC) was not nearly so bright, and in fact showed no enhancement at all, This was puzzling, given the 1988 results for the RSPIC. We attributed the lack of a radar reflection enhancement to a combination of three effects: the geometry was different; the season was different; and the intrinsic scattering from the residual ice caps was different. The 1999 experiments provided a chance to test the relative importance of these three effects, since the RNPIC would be visible in a geometry very similar to the RSPIC in 1988, and the season would also be very similar. Preliminary reduction of the data taken in 1999 shows that the RNPIC did in fact show a radar reflectivity enhancement, but that it was still not as bright as the RSPIC in 1988. Details regarding the radar reflectivity of the residual ice caps, as well as the surrounding polar layered terrains (in light of our radar reflectivity data) will be discussed. Additional information is contained in the original extended abstract.

Butler, B. J.↗

DOE OSTI2023

Low-Temperature Geothermal Geospatial Datasets: An Example from Alaska: Preprint

This project is part of a broader effort focused on geothermal heating and cooling with the purpose of demonstrating the multi-faceted value of integrating GHC/GHX into national decarbonization plans and community energy plans. Currently, there is a need to better define and evaluate low-temperature geothermal resource potential to provide the basis for supporting private sector investment. For shallow GHC/GHX fields, there is no formal compilation of subsurface characteristics shared among industry practitioners that can improve system design and operations. Alaska is specifically noted in this work, because heretofore, it has not received a similar focus in geothermal potential evaluations as the contiguous United States. The methodology consists of leveraging relevant data generate a baseline low-temperature resources (< 150 degrees C) geospatial dataset to compare and analyze information accessible to anyone trying to understand the potential of GHC/GHX and small-scale low-temperature geothermal power in Alaska (e.g., energy modelers, communities, planners, and policymakers). Importantly, this project is identifying data related to (1) the evaluation of GHC/GHX in the shallow subsurface, and (2) the evaluation of low-temperature geothermal resource availability. Additionally, data is being compiled to assess repurposing of oil and gas wells to contribute co-produced fluids toward the geothermal direct use and heating and cooling resource potential. The goal of this project is facilitating future deployment of GHC/GHX analysis and community-led programs and update the low-temperature geothermal resources assessment of Alaska. A better understanding of shallow potential for geothermal heat exchange (GHX) will improve design and operations of highly efficient GHC systems. The deployment and impact that can be achieved for low temperature geothermal resources will contribute to decarbonization goals and facilitate widespread electrification by shaving and shifting grid loads.

bottom-hole temperature↗

NASA NTRS2020

Innovative Thermal Management & Control to Surmount Challenges of Exploring Ocean Worlds on Europa and Enceladus

Probes to penetrate the thick ice shells of our solar system’s Ocean Worlds have been studied for nearly 20 years, since scientific evidence strongly suggested a subsurface ocean on the Jupiter moon called Europa. There is keen scientific interest in exploring four significant themes on such proposed missions: 1) Geodynamics, 2) Geochemistry, 3) Habitability, and 4) Life Detection. The ice shells of Ocean Worlds are predicted to be up to 40 km thick; they exhibit extreme thermal environments, with ice temperatures from 100 K to 270 K, and extreme pressure environments from vacuum to 53 MPa. Jet Propulsion Laboratory has conducted a broad-look investigation of proposed mission concepts to Europa to identify the significant technology and operational challenges of Europa ice-penetration. The thermal-mechanical system (TMS) of an ice penetration probe (IPP) designed to access the ocean of an icy moon using radioisotope thermoelectric generators for heat and power faces technological hurdles exacerbated by severe thermal and volume constraints. This study identified thermal management and control (TMC) challenges that are strongly linked to: ice penetration start-up, mobility and navigation in the ice, communications while in the ice sheet, and detecting and avoiding in-ice hazards. The major objectives of the TMC system are: 1) Absorb internal thermal energy from the IPP radioisotope power source, 2) Maintain liquid water conditions around the IPP at all times, 3) Manage and control thermal flows from probe nose to tail, and 4) Provide pressure containment for all internal probe components. This work discusses the baseline TMC system architecture and design developed to accomplish these objectives, and survive and transit the extreme ice thicknesses in pursuit of Icy/Ocean Worlds science goals. The proposed TMC system consisting of an internal pumped two-phase fluid loop “thermal bus” for thermal energy capture, variable conductance heat pipe system for passively adaptive thermal energy transport around the probe, and water jetting system for ice cutting is described and discussed. Critical testing performed to date is described.

Lee, Kuan-Lin↗

NASA NTRS1994

Sinuosity of Martian rampart ejecta deposits

The sinuosities of 2213 Martian rampart ejecta craters are quantified through measurement of the ejecta flow front perimeter and ejecta area. This quantity, called lobateness, was computed for each complete lobe of the 1582 single lobe (SL), 251 double lobe (DL), and 380 multiple lobe (ML) craters included in this study. A lobateness value of 1 indicates a circular ejecta blanket, whereas more sinuous ejecta perimeters have lobateness values greater than 1. Although resolution does have an effect on the absolute values of lobateness, the general relationships between lobateness and morphology exist regardless of resolution. Evaluation of the lobateness values reveals that the outer lobes of DL and ML craters have higher median lobateness values (i.e., are more sinuous) than the inner lobes. The outermost lobe of ML craters displays higher lobateness values than the outer lobe of DL craters or the single lobe of SL craters. Previous reports of lobateness-diameter, lobateness-latitude, and lobateness-terrain relationships for rampart craters are not supported by this study. Many of the differences between the results of this study and the previous lobateness analyses can be attributed to the inclusion of resolution effects and the distinction between different ejecta morphologies in this study. The results of this study taken together with a previous analysis of the distribution and diameter dependence of different ejecta morphologies are most consistent with the theory that Martian lobate ejecta morphologies form from impact into subsurface volatiles.

Barlow, Nadine G.↗

DOE OSTI2024

Changes in environmental and engineered conditions alter the plasma membrane lipidome of fractured shale bacteria

ABSTRACT Microorganisms that persist in fractured shale reservoirs cause several problems including secreting foul gases and forming biofilms. Current biocontrol measures often fail due to limited knowledge of their in situ activities. The plasma membrane protects the cell, mediates many of its critical functions, and responds to intracellular cues and ecological perturbations through physicochemical modifications. As such, it provides valuable insight into the physiological adaptation of microorganisms in disturbed environmental systems. Here, we (i) demonstrate how changes in salinity and hydraulic retention time (HRT) influence the plasma membrane intact polar lipid (IPL) chemistry of model bacterium, Halanaerobium congolense WG10, and mixed microbial consortia enriched from shale-produced fluids and (ii) elucidate adjustments in membrane IPL chemistry during biofilm growth relative to planktonic cells. We incubated H. congolense WG10 in chemostats under three salinities (7%, 13%, and 20% NaCl), operated under three HRTs (19.2, 24, and 48 h), and in drip flow biofilm reactors under the same salinity gradients. Also, mixed microbial consortia in produced fluids were enriched in triplicate chemostat vessels under three HRTs (19.2, 24, and 72 h) and biofilm reactors. Lipids were analyzed by ultra high performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). Our results show that phosphatidylglycerols, cardiolipins, and phosphatidylethanolamines were predominantly enriched in planktonic H. congolense WG10 cells grown at hypersalinity (20%) compared to optimum (13%). In addition, several zwitterionic phosphatidylcholines and phosphatidylethanolamines were higher in abundance during biofilm growth. These observations suggest that microbial adaptation and biofilm formation in fractured shale are enabled by strategic plasma membrane IPL chemistry adjustments. IMPORTANCE Microorganisms inadvertently introduced into the shale reservoir during fracturing face multiple stressors including brine-level salinities and starvation. However, some anaerobic halotolerant bacteria adapt and persist for long periods of time. They produce hydrogen sulfide, which sours the reservoir and corrodes engineering infrastructure. In addition, they form biofilms on rock matrices, which decrease shale permeability and clog fracture networks. These reduce well productivity and increase extraction costs. Under stress, microbes remodel their plasma membrane to optimize its roles in protection and mediating cellular processes such as signaling, transport, and energy metabolism. Hence, by observing changes in the membrane lipidome of model shale bacteria, Halanaerobium congolense WG10, and mixed consortia enriched from produced fluids under varying subsurface conditions and growth modes, we provide insight that advances our knowledge of the fractured shale biosystem. We also offer data-driven recommendations for improving biocontrol efficacy and the efficiency of energy recovery from unconventional formations.

03 NATURAL GAS↗

NASA NTRS2015

Noctis Landing: A Proposed Landing Site/Exploration Zone for Human Missions to the Surface of Mars

The proposed Noctis Landing Site/Exploration Zone (LS/EZ) is shown in Figure 1. Our preliminary study suggests that the proposed site meets all key Science and Resources (incl. Civil Engineering) requirements. The site is of significant interest, as the EZ not only offers a large number and wide range of regions of interest (ROIs) for short-term exploration, it is also located strategically at the crossroads between Tharsis and Valles Marineris, which are key for long-term exploration. The proposed site contains Regions of Interest (ROIs) that meet the following Science requirements: -‐ Access to (1) deposits with a high preservation potential for evidence of past habitability and fossil biosignatures and (2) sites that are promising for present habitability. The site presents a wide variety of ROIs qith likely aqueous features and deposits, including sinous channels and valleys, slope gullies, lobate debris aprons, impact craters with lobate ejecta flows, and "bathtub ring" deposits. Neutron spectrometry also suggests hydrogen is present within the topmost 0.3 m or so of 4 to 10 wt% WEH (Water Equivalent Hydrogen). -‐ Noachian and/or Hesperian rocks in a stratigraphic context that have a high likelihood of containing trapped atmospheric gases. Collapsed canyon rim material with preserved stratigraphy is abundantly present and accessible. -‐ Exposures of at least two crustal units that have regional or global extents, that are suitable for radiometric dating, and that have relative ages that sample a significant range of martian geological time. Canyons floors in Ius Chasma, Tithonium Chasma, and plateau tops on Tharsis and in Sinai Planum offer access to distinct crustal units of regional extent. -‐ Access to outcrops with linked morphological and/or geochemical signatures indicative of aqueous or groundwater/ mineral interactions. Iron and sulfur-bearing deposits on canyon floors in Noctis Labyrinthus, and in Ius Chasma (IC) and Tithonium Chasma (TC) offer many such outcrop options. -‐ Identifiable stratigraphic contacts and cross-cutting relationships from which relative ages can be determined. In place and collapsed canyon walls in NL, TC, and IC offer such opportunities. -‐ Other types of ROIs include access points to surrounding plateau top areas for longer term regional exploration. A key attribute of the proposed Noctic Landing site is its strategic location to allow the shortest possible surface excusions to Tharsis and Valles Marineris (VM). VM is the feature and region on Mars that exposes the longest record of Mars' geology and evolution through time. Tharsis is the region of Mars that has experienced the longest and most extensive volcanic history, and might still be volcanically active. Some of the youngest lava flows on Mars have been identified on the western flanks of the Tharsis Bulge, i.e., within driving range of future longrange (500 - 1000 km) pressurized rover traverses. The proposed site also contains ROIs that offer the following Resources (incl. Civil Engineering) characteristics: -‐ Access to raw material that exhibits the potential to (1) be used as feedstock for water-generating in situ resource utilization (ISRU) processes and (2) yield significant quantities (greater than 100 MT) of water. The raw material is likely in the form of hydrated minerals, and possibly ice/regolith mix. The top of the raw material deposit is at the surface. -‐ Access to a region where infrastructure can be emplaced or constructed. This region is less than 5 km from the LS and contains flat, stable terrain. The region exhibits evidence for an abundant source of loose regolith. Several deep pits in the area combined with the availability of sand suggests that some natural terrain features can be adapted for construction purposes. -‐ Access to raw material that exhibits the potential to be used as metal feedstock for ISRU and construction purposes. Iron and sulfur-rich mineral surface deposits have been identified in CRISM data in many locations in this area. Noctis Landing is the lowest-altitude location on Mars that straddles both the Tharsis region (above average geothermal gradients) and Valles Marineris (minimal crustal thickness from surface (valley floor) to a subsurface liquid water table. Noctis Landing has the potential for being an ideal site for eventual deep drilling on Mars to access deep subsurface liquid water and potentially encountering extant life. Available data remains insufficient to fully qualify the Noctis Landing site. Additional remote sensing data (visible, Near and Mid-IR, and radar) and surface reconnaissance via a high-mobility robotic rover are recommended. In particular, it will be important to assess the trafficability of the site, and its potential for yielding water and metals as a resource. Access to plateau tops from the Noctis Landing site on the canyon floor should be demonstrated. Future exploration of the site would also be enhanced significantly by the availability of robotic (tele-operatable) surveying and sample-collecting drones. Testing of the use of such collaborative science and exploration technologies should be conducted at terrestrial sites such as the Haughton-Mars Project site on Devon Island, High Arctic, among others. Note: Noctis Landing is not an official Mars nomenclature name for this location. Because the area of the proposed LS/EZ had no name, and because it is close to Noctis Labyrinthus to the West while being distinct from it, the provisional name Noctis Landing is proposed. Noctis means night in Latin.

Lee, Pascal↗

DOE OSTI2025

Hydrologic Model Data for the East Fork Poplar Creek Watershed Simulated with the Advanced Terrestrial Simulator (ATS): Streamflow and Network Expansion–Contraction Dynamics

This dataset supports hydrologic modeling and stream network expansion–contraction analysis for the East Fork Poplar Creek (EFPC) Watershed in Tennessee. It includes a Jupyter notebook for model setup, model configuration files, simulation outputs, and derived products used to evaluate model performance and investigate stream dynamics under varying hydrologic conditions. The dataset was generated using the Watershed Workflow Python package and the Advanced Terrestrial Simulator (ATS), enabling integrated surface–subsurface hydrologic simulations using a stream-aligned mesh. Outputs include high-resolution time series of streamflow, active network length, water table depth, and related hydrologic variables. Also included are spatially explicit stream persistency indices and classifications of reaches as perennial or non-perennial. These data facilitate reproducibility and support further research on stream intermittency and variability in network extent.The model data archive is organized in following directories:1) model_setup_inputsContains the Watershed Workflow Jupyter notebooks (accessed through any open source code editor), selected input datasets, and resulting ATS input files, including XML files (access through any open source code editor), computational mesh (.exo files can be viewed using Paraview), and meteorological forcing files (.h5 files can be accessed through h5py python package and HDFView open source software). 2) model_outputsIncludes ATS simulation outputs relevant to this study. Time series of spatially integrated or averaged variables (e.g., streamflow, water table depth) are provided as CSV files. Select spatial fields (e.g., ponded depth and water table depth) are saved as pickled Python objects to reduce file size, and can be accessed through pickle package in Python. Key geometry objects from Watershed Workflow—such as the surface mesh and river tree—are also included to support analysis of streamflow persistency and expansion–contraction dynamics. These files can also be accessed through Watershed Workflow Python package.3) model_evaluationProvides observed streamflow time series and field survey-based flow regime classifications used to evaluate model performance. Jupyter notebooks for processing ATS outputs and comparing model predictions with observations to build confidence in the model prior to scientific analysis are also included.4) Q_L_relationshipsContains workflows for generating time series of discharge, active network length, and related hydrologic variables used in the stream network expansion–contraction analysis. Includes routines for delineating baseflow-dominated periods. For each catchment, notebooks and processed data (as pickled DataFrames accessed through Pandas Python package) are provided. 5) figure_scriptsProvides the Jupyter notebooks used to generate the figures presented in the paper.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2005

Using Wind Driven Tumbleweed Rovers to Explore Martian Gully Features

Gully features have been observed on the slopes of numerous Martian crater walls, valleys, pits, and graben. Several mechanisms for gully formation have been proposed, including: liquid water aquifers (shallow and deep), melting ground ice, snow melt, CO2 aquifers, and dry debris flow. Remote sensing observations indicate that the most likely erosional agent is liquid water. Debate concerns the source of this water. Observations favor a liquid water aquifer as the primary candidate. The current strategy in the search for life on Mars is to "follow the water." A new vehicle known as a Tumbleweed rover may be able to conduct in-situ investigations in the gullies, which are currently inaccessible by conventional rovers. Deriving mobility through use of the surface winds on Mars, Tumbleweed rovers would be lightweight and relatively inexpensive thus allowing multiple rovers to be deployed in a single mission to survey areas for future exploration. NASA Langley Research Center (LaRC) is developing deployable structure Tumbleweed concepts. An extremely lightweight measurement acquisition system and sensors are proposed for the Tumbleweed rover that greatly increases the number of measurements performed while having negligible mass increase. The key to this method is the use of magnetic field response sensors designed as passive inductor-capacitor circuits that produce magnetic field responses whose attributes correspond to values of physical properties for which the sensors measure. The sensors do not need a physical connection to a power source or to data acquisition equipment resulting in additional weight reduction. Many of the sensors and interrogating antennae can be directly placed on the Tumbleweed using film deposition methods such as photolithography thus providing further weight reduction. Concepts are presented herein for methods to measure subsurface water, subsurface metals, planetary winds and environmental gases.

Antol, Jeffrey↗

NASA NTRS2004

Using Wind Driven Tumbleweed Rovers to Explore Martian Gully Features

Antol, Jeffrey↗

DOE OSTI2025

Phase II Field Demonstration at Lansing Smith Generating Plant, Southport, Florida (Final Report)

The Final Technical Report: Field Test Design and Pressure Management Strategies for Phase II Field Demonstration of Optimal Pressure Control, Plume Management, and Produced Water Strategies presents the culmination of multi-year efforts under the U.S. Department of Energy’s Brine Extraction and Storage Test (BEST) program, led by EPRI in partnership with Gulf Power at the Plant Smith site near Panama City, Florida. The project was designed to evaluate and demonstrate the technical feasibility of managing subsurface pressures and fluid movement associated with large-scale CO₂ injection, using low-salinity water as a proxy fluid. Through a combination of field injection testing, reservoir modeling, and optimization studies, the research team developed and refined active and passive brine extraction strategies aimed at controlling injection-induced pressure buildup, mitigating risks of fault activation, and managing plume migration. The field demonstration incorporated a new injection well (TIW-2), a new monitoring/extraction well (TEMW-A), and an existing well (TIW-1) repurposed for passive pressure relief. Complementary geophysical monitoring was designed to track plume development and evaluate the effectiveness of pressure management strategies. The report details the integrated workflow encompassing wellfield development, regulatory permitting, model calibration, and survey design. It includes a comprehensive electromagnetic (EM) modeling and inversion study used to develop a cost-effective, time-lapse geophysical monitoring plan capable of imaging the evolving low-salinity plume within the highly saline Lower Tuscaloosa formation. Reservoir simulation results guided the design of the 17-month injection program and the timing of active extraction to maintain formation pressures below a hypothetical fault reactivation threshold. Supporting analyses evaluated potential injectivity risks related to clay fines migration and geochemical incompatibility, identifying practical mitigation measures such as salinity and pH control. The final design recommends a series of crosswell EM surveys—conducted before, during, and after injection—as the most effective and economical approach for plume imaging, supplemented by continuous downhole pressure and flow monitoring. Collectively, the study provides a field-ready framework for cost-effective pressure management and monitoring in support of future CO₂ storage projects.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2025

Groundwater flowpath characteristics drive variability in per- and polyfluoroalkyl substances (PFAS) loading across a stream-wetland system

Groundwater dependent ecosystems in areas with industrial and military land use are at risk of direct exposure to a wide range of contaminants, including PFAS chemicals. Glaciated terrain often has mixed high and low permeability sediments coupled with groundwater flow-through lake features. These hydrogeologic attributes create highly complex ‘source to seep’ dynamics that make spatiotemporal contaminant transport patterns difficult to predict. We investigated one such system in detail using a suite of heat-tracing and chemical methods. Numerous (n=57) preferential groundwater discharge zones (vertical flux rates ranging 0.2 to 3.2 m/d) were identified across the upper Quashnet River stream-wetland system in Mashpee, MA, USA, adjacent to an Air Force Base with several known PFAS source areas. Surface-water and groundwater samples were collected and analyzed (n=145) for precursors and terminal PFAS compounds between March and September 2022. Samples were collected at identified seeps along the Quashnet River (n=59), from wells upgradient from the stream-wetland system (n= 44), from contributing flow-through kettle lakes (n=8), and at multiple locations along the Quashnet River (n=34). Samples from seeps and wells had measured PFAS concentrations ranging from non-detect to approximate 3,500 ng/L (mean= 1,650 ng/L), and a range of deuterium excess values (3.2 to 15.9 per mil) indicative of varying degrees of groundwater-lake interaction prior to emergence at the discharge zones. Groundwater-lake interaction along flowpaths that sourced the sampled seeps was farther supported by significant correlations (p < 0.01) between deuterium excess and %PFAS precursors, and between %PFAS precursors and multiple terminal PFAS compounds (e.g., PFPeS, PFBS, PFHxS). However, some sampled seeps contributing groundwater to the stream-wetland system had much higher total PFAS concentrations (>1000 ng/L) than the upgradient kettle lakes, despite showing lake (evaporative) isotopic signatures, indicating the potential for groundwater flowpath convergence at wetland discharge zones and the influence of lakebed PFAS precursor reactions. PFAS compounds and water isotopic composition at sampled multilevel groundwater wells, rivers, lakes, and seeps suggest that a complex mixture of source groundwater and flowpath characteristics are responsible for diverse observed PFAS mixtures at preferential discharge zones across the stream-wetland system. Further, total PFAS loading patterns to the Quashnet River via groundwater discharge remained remarkably similar from winter to summer to fall conditions, despite a regional dry period in late summer 2022 with the upper river channel completely drying. This work addresses gaps in the existing PFAS literature by demonstrating the importance of subsurface fate and transport on PFAS compound concentrations in controlling contaminant mass loading in preferential groundwater discharge zones and presents a transferrable field toolkit for efficient characterization of spatially preferential PFAS transport dynamics.

Contaminant transport↗

DOE OSTI2022

Comparison of Microbial Profiling and Tracer Testing for the Characterization of Injector-Producer Interwell Connectivities

Insufficient understanding of the microbial communities and associated microbial processes in geological reservoirs hinders the utilization of this rich data source for improved resource management. In this study, along with four interwell tracer tests at a 1478-m deep fractured crystalline-rock aquifer, we analyzed the microbial communities in the injected and produced water by high-throughput sequencing. The microbial community similarities across boreholes during an interwell flow scenario frequently encountered in reservoir development was explored. Despite the significant tracer recoveries (~30%) in all tracer tests and the cumulatively >100,000 L of exogenous water (carrying exogenous microbes) injected into the 10-m-scale reservoir, the overall structure of produced-fluid microbiome did not increasingly resemble that of the injectate. However, producers with better connectivity with the injector (based on tracer test results) did have more amplicon sequence variants (ASVs) that overlapped with those in the injectate. We identified possible drivers behind our observations and verified the practicality of repeated microbial sampling in the context of reservoir characterization and long-term monitoring. We highlight that injector-producer microbial profiling could provide insights on the relative connectivities across different producers with a given injector, and that the associated logistical needs may be comparable or even less than that of classic tracer tests.

59 BASIC BIOLOGICAL SCIENCES↗