SEARCH · Search NASA

Results for “subsurface flow”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

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 OSTI2019

Subtask 1.1 – Advanced Characterization of Unconventional Oil and Gas Reservoirs to Enhance CO2 Storage Resource Estimates

The Energy & Environmental Research Center (EERC) developed methods and approaches for advanced characterization of unconventional oil and gas reservoirs to enhance CO2 storage resource estimates using samples collected from organic-rich shales of the Bakken Formation. The project was funded under U.S. Department of Energy (DOE) Cooperative Agreement No. DE-FE0024233, with nonfederal cost share provided by Hitachi High Technologies America. The project participants included the EERC, Hitachi High Technologies America, the team at the National Energy Technology Laboratory’s (NETL’s) computed tomography (CT) scanning laboratory in Morgantown, West Virginia, and several NETL staff who have worked on CO2 storage resource estimation. Investigations on the use of subsurface geologic formations for CO2 storage have been ongoing for many years. In the last decade, improvements in directional drilling technology and hydraulic fracturing of organic-rich shale and other tight rock formations have opened new opportunities for CO2 storage, with the potential cobenefit of enhancing gas and oil recovery. Gas-bearing, organic-rich shales have been of particular interest as a carbon storage target because their ability to store large volumes of CO2 is implied by their ability to store methane; however, estimating the CO2 storage resource of organic-rich shales and other unconventional reservoirs is challenging because of the lack of detailed geologic and petrophysical data that are needed to develop improved volumetric equations. According to NETL’s method for estimating regional-scale CO2 storage resources, the key factors that control CO2 storage in geologic formations include void space within induced fractures, natural fractures and matrix pores, as well as the occurrence and distribution of organic matter (OM) and clay minerals. A key goal of this effort was to use field emission scanning electron microscopy (FESEM) analysis of fine-grain unconventional reservoir rocks coupled with advanced image analysis and machine learning techniques to quantify the occurrence and distribution of porosity, OM, and clay minerals. By better understanding the potential migration pathways of CO2 within tight rocks and quantifying the occurrence of sorptive components (OM and clays) along potential flow pathways, methods to estimate the CO2 storage resource of such rocks can be refined. The Bakken petroleum system (Bakken) was used as the source of rock samples for this effort because it includes organic-rich shales with contrasting degrees of thermal maturity as well as tight, nonshale clastic rocks with varying clay mineral content. To better understand the key components of interest in tight rock samples, the EERC worked closely with Hitachi to optimize the FESEM settings necessary to resolve the sample mineralogy, porosity and OM content at micro- to nanoscale resolution. The EERC also worked closely with Hitachi to apply and improve its advanced mineral identification and characterization system (AMICS) software so that it is better suited for mineral identification at the high resolution necessary to resolve features of interest in fine-grained rocks such as organic-rich shales. The EERC also applied advanced image analysis software packages, including Ilastik, to classify and quantify the key components of interest within the samples using FESEM imagery. These advanced techniques allowed for quantification and differentiation of the void space within fractures, pores and OM, as well as a technique to determine and quantify the key mineral and OM components surrounding pore spaces. Being able to better quantify the void space and the potential accessibility of CO2 to sorptive minerals or OM within tight rocks allows for improved estimation of the CO2 storage resource potential of organic-rich shales and other unconventional formations. The information derived from the image analysis techniques was used to refine NETL’s method for estimating regional-scale CO2 storage resources in collaboration with staff from NETL. An additional goal of this work was to better understand rates of CO2 permeation, flow pathways, and sorption in organic-rich shales and nonshale reservoir rocks and to evaluate potential changes in CO2 flow over time. The EERC conducted long-duration (several weeks) CO2 flow-through testing at reservoir conditions using Bakken shale and nonshale samples to assess CO2 permeation rates and retention/sorption within core plugs to represent bulk matrix behavior. Conventional isothermal adsorption testing over a range of pressure conditions was also performed on pulverized Bakken samples to better assess potential maximum sorptive capacity of the samples. The results of this testing illustrated the higher sorptive capacity of the Bakken shales when compared to the sorptive capacity of the nonshale reservoir rock, which contains very little OM when compared to the Bakken shales. The results also illustrate the differences in permeation behavior as a function of fracture flow versus matrix flow. In an attempt to visualize CO2 permeation and flow behavior in Bakken shale samples over time, the EERC collaborated with NETL’s CT Lab in Morgantown, West Virginia. A variety of tests were performed at the NETL lab on unfractured and fractured samples of the Bakken shale. The unfractured Bakken shale samples did not result in CO2 migration within the samples over the testing duration; thus imagery of CO2 migration pathways within the matrix of Bakken samples was not obtained. The results of the flow-through testing within fractured samples illustrate changes in CO2 permeation rates over time, with a trend of decreasing CO2 flow over time. The results of this effort suggest that FESEM analysis coupled with advanced image analysis and machine learning techniques are effective in characterizing and quantifying the mineralogy, porosity, fractures, and OM content of 2-D imagery at the high resolution needed to identify features of interest in tight rocks. In addition, the image analysis techniques allow for identification and quantification of different porosity types as well as the key mineral components and/or OM that CO2 would contact in the reservoir as it migrates through pores and fracture networks. This is key information needed to better estimate the CO2 storage resource potential of organic-rich shales and other tight reservoirs. Additional work is needed to determine the sorptive capacity of the individual sample components (OM and clay minerals) at reservoir conditions so that this approach and the refined CO2 storage resource equation can be applied to other shale and unconventional reservoirs. In addition, subsequent work should be performed to better understand the number of images needed to use stereological approaches to estimate the volume of sample components using 2-D imagery, to explore machine learning techniques to differentiate between induced vs. natural fractures, and to estimate the connectivity of pores and fracture networks using 2-D imagery. Being able to derive this information at the microscale would facilitate scaling up of the sample characterization data used for reservoir modeling and simulation to estimate the CO2 storage resource potential of tight formations.

03 NATURAL GAS↗

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↗

DOE OSTI2019

Infiltration Data Package for the E-Area Low-Level Waste Facility Performance Assessment

The infiltration data package contains the input parameters, the cap design and material properties assumptions, and the modeling results for the Hydrologic Evaluation of Landfill Performance (HELP) infiltration model simulations performed in support of the E-Area Low-Level Waste Facility Performance Assessment (PA). The infiltration estimates establish the upper boundary condition for the PORFLOW vadose-zone model and GoldSim model simulations for the following E-Area disposal unit types: Slit and Engineered Trenches, Low-Activity Waste Vault (LAWV), Intermediate-Level Vault (ILV), Componentin- Grout (CIG) special waste form trench segments, and the Naval Reactor Component Disposal Areas (NRCDAs). The infiltration data package builds upon relevant, foundational PA technical reports and memoranda from the past 15 years and is supported by three important components: the HELP model input parameter datasheets, HELP model input and output filenames and directory structure, and infiltration rates as a function of time for each scenario for each disposal unit type. For Slit and Engineered Trenches, a single best-estimate Bahia-grass case of two-percent slope and 585- foot slope length represents an upper bound on intact closure-cap infiltration rates over a 10,000-year period. Building upon the single, bounding, and intact case, two different approaches were considered for incorporating the effect of localized cap subsidence on infiltration rates in the PORFLOW flow and transport model simulations. In the first approach, a Monte Carlo probabilistic model was developed to calculate infiltration rates for 0.54%, 2%, 3.6%, and 4.9% subsidence scenariosa reflecting historical and future non-crushable content for the Slit and Engineered Trenches. The resulting blended (spatially and slope-length averaged) infiltration rates become an inflow upper boundary condition for the PORFLOW vadose zone model simulations. In the second and preferred approach, a Monte Carlo rectangle packing algorithm (Danielson, 2019) and a simplified equation for the total mass flux of water draining into a subsided compartment (Equation 1) will be employed to implement a weighted blending of radionuclide fluxes to the water table using the results of deterministic PORFLOW vadose zone simulations representative of specific subsidence cases. Finally, HELP-model results from sensitivity studies of the intact infiltration case were fit to a log-logistic function to generate infiltration profiles over a 10,000-year period for most-optimistic, more-optimistic, best-estimate, more-pessimistic, and most-pessimistic cases, for both the intact and four subsidence scenarios. For the LAWV and IL vaults, infiltration rate profiles for a 10,000-year period were developed for both an on-vault (above the concrete/metal vault roof) and an off-vault (10-foot soil zone adjacent to vault walls) scenario. The purpose of the off-vault simulations was to confirm that subsurface runoff from the concrete vault roof will adequately drain through the lowermost backfill layers adjacent to the vault walls. The LAWV and ILV infiltration profiles assume collapse of the concrete roof at relative Years 2,905 and 7,100, respectively. The identical Bahia-grass case with two-percent slope and 585-foot slope length assumed for the Slit and Engineered Trenches will represent an upper bound on intact infiltration rates for the NRCDAs and CIG special waste form as well. While the NRCDA casks are assumed to remain structurally intact (i.e., non-assumed to collapse 300 years after the end of operations. Because both rows of the existing CIG trench segments in ST23 (formerly CIG01) extend across the planned cap crest, the subsided CIG trench segments are treated conservatively in the HELP model as a 100%-subsidence case with zero subsurface run-on. The CIG special waste form trench segment infiltration profile, therefore, consists of an intact-cap period from relative Year 100 to Year 300, followed by a 100%-subsidence period from relative Year 300 to Year 10,100.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

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 OSTI2018

Development of 3D Geological Model of Tuscarora Sandstone for Feasibility of Deep Direct-Use Geothermal at West Virginia University Main Campus

The subsurface uncertainty at West Virginia University Main Campus is dominated by the uncertainty in the projections of geofluid flowrate in the target formation, the Tuscarora Sandstone. In this paper, three cores from the heterogeneous reservoir, available through West Virginia Geologic and Economic Survey, are analyzed by performing core analysis using CT scanning and permeability measurements via minipermeameter. Additional geological data are collected through cores, published literature, seismic data, and nearby, existing wells to estimate thickness, fracture network configuration and geothermal gradient to minimize the uncertainty of well deliverability. Using these estimated reservoir properties; a 3D conceptual model for the proposed geothermal site is developed. This dataset includes a GRC conference publication along with the data used to produce results explained in the paper including minipermeability measurement data for Preston -119 core and thin section analysis photos and data for Clay-513 core.

15 GEOTHERMAL ENERGY↗

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↗

DOE OSTI2021

Continental-Scale Geographic Trends in Barometric-Pumping Efficiency Potential: A North American Case Study

Barometric pumping is a gas transport mechanism that has important implications for many applications involving subsurface gas seepage processes. Here, this study provides the first continental-scale analysis of barometric-pumping efficiency potential based on meteorology. We quantified the barometric-pumping efficiency potential at 1,257 locations across the continental US and Canada. The results provide continental-scale geographic dependencies of barometric-pumping efficiency potential, indicating a significant correlation with latitude and a nonlinear dependence on longitude. The analysis also indicates that variability in barometric-pumping efficiency potential decreases with distance from the coast and as elevation increases. Locations far from the coastline are more likely to have upper mid-range potentials, while higher elevation locations are more likely to have low potentials. The highest barometric-pumping efficiency potentials are mostly found around the Gulf of St. Lawrence around 50°N. Locations along the Atlantic coast exhibit large-scale variations in potentials with a clear increasing trend with latitude.

58 GEOSCIENCES↗

DOE OSTI2020

Distribution of Infiltration in the 216-U-10 and 216-B-3 Pond Systems 1944-1997

A modeling system to simulate fate and transport of contaminant plumes in groundwater beneath the Central Plateau to support the Cumulative Impact Evaluation (CIE) and Composite Analysis (CA) is being developed. The CIE will evaluate the cumulative effect of multiple Central Plateau radiological and chemical sources on groundwater quality (DOE/RL-2018-69, Cumulative Impact Evaluation Technical Approach Document ). This will enable cleanup decisions to be evaluated in the context of all Central Plateau waste sites and existing groundwater conditions. The CA evaluates the cumulative effect of all sources of radioactive materials on radiological dose to future members of the public (DOE/RL-2018-60, Annual Status Report (FY 2018): Composite Analysis of Low Level Waste Disposal in the Central Plateau of the Hanford Site ). Fate and transport modeling will be used to forecast future groundwater concentrations for both the CIE and CA. The CIE/CA subsurface modeling system is divided into source (inventory), vadose zone, and saturated zone facets. Source inventory is a necessary input to the vadose zone facet, and releases from the vadose zone to groundwater are inputs to the saturated zone facet. The vadose zone facet simulates historical and future releases from the vadose zone to groundwater for both liquid waste and solid waste disposal sites. Effluent release volumes and radionuclide inventory estimates for liquid waste sites are contained in ECFHANFORD- 17-0079, Hanford Soil Inventory Model (SIM-v2) Calculated Radionuclide Inventory of Direct Liquid Discharges to Soil in the Hanford Site’s 200 Areas , herein referred to as SIM-v2. The liquid waste sites in SIM-v2 include the 216-B-3 Pond near 200 East Area and the 216-U-10 Pond in 200 West Area. The 216-B-3 Pond system consisted of a main pond and 3 expansion lobes that were operational during overlapping time periods. It also included influent ditches, but these are not listed as liquid waste sources in SIM-v2 and are not included in the CIE/CA modeling. The total effluent volume released to the pond system was measured, but the partitioning of this volume to the various expansion lobes was not determined and is not accounted for in SIM-v2. Regarding the 216-U-10 Pond, it received water from several influent ditches. The pond and the influent ditches are listed as liquid waste sites in SIM-v2 and are included as sources in the CIE/CA modeling. Because water in the ditches flowed into the pond, only a portion of the effluent volumes and contaminant inventories assigned to the ditches in SIM-v2 infiltrated from the ditches. Thus, partitioning of infiltration and contaminant inventories between the ditches and the pond is not fully accounted for in SIM-v2. Further, some of the water in the 216-U-10 Pond entered overflow ditches and the amount of this overflow was not recorded and is not accounted for in SIM-v2. Thus, the purpose of this Environmental Calculation File (ECF) is to partition the infiltration of wastewater and contaminant inventory between the main pond and each expansion lobe of the 216-B-3 Pond system, and to partition the infiltration of wastewater and contaminant inventory between the ditches and pond of the 216-U-10 Pond system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1999

Lateral and Vertical Heterogeneity of Thorium in the Procellarum KREEP Terrane: As Reflected in the Ejecta Deposits of Post-Imbrium Craters

The Procellarum KREEP Terrane displays the highest concentrations of Th on the Moon. However, locations of elevated Th in this region appear to be random. As observed in the 5 deg per pixel equal-area Th data, and made more evident in the preliminary 2 deg data, Th is enhanced around the craters Aristillus, Aristarchus, Kepler, Mairan, the Apennine Bench formation, and the Fra Mauro region, while noticeably and unexpectedly lower in other locations (e.g., Archimedes, Copernicus, Eratosthenes, and Plato). We have examined the composition of the materials present in these regions with the goal of understanding the patchy nature to the distribution of Th and ultimately to decipher the geologic processes that have concentrated the Th. At present time, the published resolution of the Lunar Prospector Th gamma-ray data is low (5 deg per pixel), but this will soon be superceded by significantly higher-resolution data (2 deg per pixel). Even at this improved resolution, however, it is difficult to resolve the units that are the major source of Th. In an attempt to circumvent this problem, we employ the higher-resolution Clementine multispectral data for those regions mentioned above. We use the UV-VIS-derived compositional information and the spectral properties of craters, and their ejecta as drill holes through the mare-basalt surface to investigate the thickness and composition of underlying material. With this information we attempt to piece together the stratigraphy and geologic history of the Imbrium-Procellanim region. We processed the five-band multispectral data from the Clementine Mission (415, 750,900,950, and 1000nm) using ISIS software and calibration parameters developed by the USGS, Flagstaff, Arizona. Final image mosaics are in equal-area sinusoidal projection, and have a resolution of 250 m/pixel. Using the method of we produced maps of FeO and Ti02 composition. Here we examine the Th, FeO, and Ti02 composition and spectral properties of the craters discussed above and their ejecta, with the goal of describing the materials they excavate. One interpretation for the origin of the high-Th material is that subsurface KREEPy materials have been excavated by impact craters. The material excavated may be either volcanic KREEP (e,g., Apennine Bench Formation), KREEPy impact-melt breccia formed by the Imbrium impact (e.g., Fra Mauro Formation), or other KREEP-rich crustal material. Determining which type of material is responsible for the elevated Th and its extent is important to understanding the premare and possibly the prebasin stratigraphy of the Imbrium-Procellarum Region. Merging the 5 deg. Th data with the shaded relief map, we observe that the highest Th concentrations are not related to pre-Imbrium upper crustal materials. The Apennines, Alpes, and Caucasus Mountains represent the pre-Imbrian highlands material and do not express concentrations of Th, FeO, and TiO2 as high as the most Th-fich materials exposed within the Procellarum KREEP Terrane. We observe that, in general, these massifs contain 10-14 wt% FeO and 4-7 ppm Th. Determining whether the Th signal is from KREEP basalts or KREEPy impact-melt breccias cannot be done with the Clementine data because the two rock types are compositionally and mineralogically too similar (e.g., the Th-rich, mafic impact-melt breccias in the Apollo sample collection are dominated by a KREEP-basalt like component. Mapping-the distribution and sizes of craters and whether they display elevated Th concentrations or not, should reveal the depth and thickness of the KREEP-rich materials, and whether they are ubiquitous (i.e., impact-melt breccia) or more randomly distributed; this might be taken as an indicator of localized KREEP-basalt flows. Within the southeastern region of the Imbrium basin, there are two Th hot spots. The first is associated with the crater Aristillus, and the latter with the Apennine Bench Formation. Adjacent to these two hot spots are craters with a lower Th signature: Archimedes and Autolycus. We observe in the ejecta of Aristillus, a region of significantly lower FeO (10-14 wt%) relative to the surrounding mare basalt. The crater Autolycus, 50 km to the south, did not excavate similar low-FeO material. We suggest that the lower-FeO material in the ejecta of Aristillus corresponds to Th-rich material; the FeO content observed in Aristillus ejecta is comparable to that of KREEP basalt or mafic impact melt breccia (10-12 wt% FeO). We determine that this low FeO, Th-rich material is volcanic KREEP, as opposed to Imbrium impact melt, on the basis that the low-Fe material is exposed more prominently in ejecta in the northern portion of Aristillus. Our assumption is that if the layer underlying Aristillus was continuous, a more widespread and uniform low-Fe signature wouldbe observed in the ejecta deposit. Archimedes, 110 km southwest of Aristillus, impacted the northern portion of the Apennine Bench prior to the eruption of KREEP basalt. Archimedes rim material is not as enriched in Th as the Apennine Bench, and there are differences between the two in FeO concentration and in their continuum slope. Archimedes exhibits a much steeperor "redder" continuum slope than the Apennine Bench. This steepslope suggests the presence of glassy material. The glassy material is concentrated around an unnamed crater on the southern rim of Archimedes (4.5W, 28.2N) and along the northern rim of Archimedes. We suggest two possibilities, or a combination of the two, to explain the low-Th signal from Archimedes: (1) The Apennine Bench prior to KREEP basalt eruption was lower in Th (4-7 ppm, e.g., similar to the Apennine massifs) and KREEP basalts are absent in the rim of Archimedes; or (2) the glassy (possibly pyroclastic) material layering the rim of the Archimedes, dilutes any high-Th material present with low-Th material. (Additional information is contained in the original)

Gillis, J. J.↗

NASA NTRS2010

Mud Volcanoes - Analogs to Martian Cones and Domes (by the Thousands!)

Mud volcanoes are mounds formed by low temperature slurries of gas, liquid, sediments and rock that erupt to the surface from depths of meters to kilometers. They are common on Earth, with estimates of thousands onshore and tens of thousands offshore. Mud volcanoes occur in basins with rapidly-deposited accumulations of fine-grained sediments. Such settings are ideal for concentration and preservation of organic materials, and mud volcanoes typically occur in sedimentary basins that are rich in organic biosignatures. Domes and cones, cited as possible mud volcanoes by previous authors, are common on the northern plains of Mars. Our analysis of selected regions in southern Acidalia Planitia has revealed over 18,000 such features, and we estimate that more than 40,000 occur across the area. These domes and cones strongly resemble terrestrial mud volcanoes in size, shape, morphology, associated flow structures and geologic setting. Geologic and mineralogic arguments rule out alternative formation mechanisms involving lava, ice and impacts. We are studying terrestrial mud volcanoes from onshore and submarine locations. The largest concentration of onshore features is in Azerbaijan, near the western edge of the Caspian Sea. These features are typically hundreds of meters to several kilometers in diameter, and tens to hundreds of meters in height. Satellite images show spatial densities of 20 to 40 eruptive centers per 1000 square km. Many of the features remain active, and fresh mud flows as long as several kilometers are common. A large field of submarine mud volcanoes is located in the Gulf of Cadiz, off the Atlantic coasts of Morocco and Spain. High-resolution sonar bathymetry reveals numerous km-scale mud volcanoes, hundreds of meters in height. Seismic profiles demonstrate that the mud erupts from depths of several hundred meters. These submarine mud volcanoes are the closest morphologic analogs yet found to the features in Acidalia Planitia. We are also conducting laboratory analyses of surface samples collected from mud volcanoes in Azerbaijan, Taiwan and Japan. X-ray diffraction, visible / near infrared reflectance spectroscopy and Raman spectroscopy show that the samples are dominated by mixed-layer smectite clays, along with quartz, calcite and pyrite. Thin section analysis by optical and scanning electron microscopy confirms the mineral identifications. These samples also contain chemical and morphological biosignatures, including common microfossils, with evidence of partial replacement by pyrite. The bulk samples contain approximately 1 wt% total organic carbon and 0.4 mg / gm volatile hydrocarbons. The thousands of features in Acidalia Planitia cited as analogous to terrestrial mud volcanoes clearly represent an important element in the sedimentary record of Mars. Their location, in the distal depocenter for massive Hesperian-age floods, suggests that they contain fine-grained sediments from a large catchment area in the martian highlands. We have proposed these features as a new class of exploration target that can provide access to minimally-altered material from significant depth. By analogy to terrestrial mud volcanoes, these features may also be excellent sites for the sampling martian organics and subsurface microbial life, if such exist or ever existed.

Allen, Carlton C.↗

NASA NTRSDate not supplied

Large Vehicle Lunar Landing Surface Interaction and In-Situ Resource Based Risk Mitigation: Landing & Launch Pads

A key capability required for the exploration of planetary bodies is the ability to land on the surface. Previous work performed by NASA and other institutions has primarily focused on landing small spacecraft on planetary surfaces and the associated small-to-medium thrusters required for the soft landing. In the case of human exploration—particularly the establishment of long duration exploration and habitation outposts—the ability to land large landers, such as the SpaceX Starship, is necessary. These larger landing systems require the use of more powerful engines, with higher engine exhaust temperatures and higher landing loads. Understanding the excavation of material by the engines, as well as the potential for the landing legs to sink into the subsurface, is key in ensuring reliable and safe landings. A further improvement in landing reliability can be achieved by constructing landing / launch pads, especially with in-situ resources. Some material excavation by the plume is inevitable, leaving at least a portion of the surface scoured and uneven under the lander and ejecting regolith particles and rocks at very high velocities. One possible solution would be to robotically build landing / launch pads (ideally autonomously) at the destination using in-situ materials. In this case, the first one or few landers will need to land on unimproved surfaces at higher risk; however, they would bring the required equipment to build the landing pads with mostly local resources, thus increasing the reliability of safe landing for subsequent larger landers. A number of methods to build in-situ landing and launch pads have already been developed. These methods include no, or some, addition of required binder additives to the local regolith material, different processing approaches and result in varying landing pad strengths. A sub-scale rocket engine plume, was used to simulate some of the conditions of a landing on the Moon to assess the effectiveness of various materials for an in-situ built landing pad, The GO2/GCH4rocket engine fired on a 1m2area coupons of representative pad materials. The results will allow continued development towards materials that satisfy the landing pad properties required for the effective risk reduction and increased reliability for landing people and equipment on the lunar surface. This work contained two parts: (1) computer modeling of a large rocket engine plume interacting with regolith on the Moon, using the Granular Gas Flow Solver (GGFS) provided by CFD Research Corporation as well as other computational fluid dynamics codes (CFD) such as Loci/CHEM. (2) Developing landing/launch pad materials that could be used for in-situ construction on the lunar surface in the future, to mitigate the calculated effects of a large vehicle rocket engine landing and launching on the Moon.

Lunar↗