Search NASA⌕ Search

SEARCH · Search NASA

Results for “subsurface imaging”

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

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

At least 559 records · Page 31

Understanding Fracture Aperture and Permeability Evolution due to Carbonate Mineralization Utilizing 3D Printing

Caprock formations are a crucial part of subsurface-engineered systems. Composed largely of shale, caprocks act as natural barriers that prevent the upward migration of fluids, thereby ensuring the containment of stored substances in subsurface formations. Fractures in these formations are potential leakage pathways for stored fluids. Mineral precipitation reactions in these fractures, particularly calcite, can significantly restrict the fluid permeability, reducing leakage potential. However, predictive capabilities of mineral precipitation in fractures and associated permeability evolution are limited due to a lack of fundamental understanding of such reactions in natural samples, complicated by mineral heterogeneity and the complexity of the fracture structure. In this study, 3D-printed fracture samples are used to understand the impact of carbonate mineralization on fracture aperture and permeability evolution. Samples were printed using a digital light processing (DLP) 3D printer and commercial liquid resin. Calcite precipitation was first tested on printed 2D films before conducting plug flow column experiments aimed to understand fracture permeability changes due to mineral precipitation. Contact angle measurement and Fourier transform infrared (FTIR) spectroscopy on printed 2D films show evidence of a substantial amount of surface energy for calcite nucleation and precipitation. Surface topography analysis of printed fractured surfaces reveals comparable values, highlighting the high replicability of the printed samples. During the column experiments, the permeability reduces exponentially due to a decrease in fracture aperture. Reductions in fracture aperture estimated from effluent concentration and 3D X-ray computed tomography (CT) show comparable results. Moreover, 3D X-ray CT images suggest the impact of local flow velocities on precipitation. The insights gained from this research contribute to a deeper understanding of the permeability evolution due to carbonate mineralization in caprock formations.

3D printing↗

Estimating biofuel contaminant concentration from 4D ERT with mixing models

Here, we present the results of a lab-scaled feasibility study to assess the performance of electrical resistivity tomography for detection, characterization, and monitoring of fuel grade ethanol releases to the subsurface. Further, we attempt to determine the concentration distribution of the ethanol from the electrical resistivity tomography data using mixing-models. Ethanol is a renewable fuel source as well as an oxygenate fuel additive currently used to replace the known carcinogen methyl tert-butyl ether; however, ethanol is preferentially biodegraded and a cosolvent. When introduced to areas previously impacted by nonethanol-based fuels, it will facilitate the persistence of carcinogenic fuel compounds like benzene and ethylbenzene, as well as remobilize them to the ground water. These compounds would otherwise be retained in the soil column undergoing active or passive remediation processes such as soil vapor extraction or natural attenuation. Here, we introduce ethanol to a saturated Ottawa sand in a tank instrumented for four-dimensional geoelectrical measurements. Forward model results suggest pure phase ethanol released into a water saturated silica sand should present a detectable target for electrical resistivity tomography relative to a saturated silica sand only. We observe the introduction of ethanol to the closed hydraulic system and subsequent migration over the duration of the experiment. One-dimensional and three–dimensional temporal data are assessed for the detection, characterization, and monitoring of the ethanol release. Results suggest one-dimensional geoelectrical measurements may be useful for monitoring a release, while three-dimensional geoelectrical field imaging would be useful to characterize, monitor, and design effective remediation approaches for an ethanol release, assuming field conditions do not preclude the application of geoelectrical methods. We then attempt to use predictive mixing models to calculate the distribution of ethanol concentration within the measurement domain. For this study we examine four different models: a nested parallel mixing model, a nested cubic mixing model, the complex refractive index model (CRIM), and the Lichtenecker-Rother (L-R) model. The L-R model, modified to include an electrical formation factor geometry term, provided the best agreement with expected EtOH concentrations.

09 BIOMASS FUELS↗

Striving to translate shale physics across ten orders of magnitude: What have we learned?

Shales will play an important role in the successful transition of energy from fossil-based resources to renewables in the coming decades. Aside from being a significant source of low-carbon intensity fuels, like natural gas, they also serve as geologic seals of subsurface formations that may be used to isolate nuclear waste, sequester CO 2 , or store intermittent energy (e.g., solar hydrogen). Despite their importance, shales pose significant engineering and environmental challenges due to their nanoporous structure and extreme heterogeneity that spans at least ~10 orders of magnitude in spatial scale. Two challenges inhibit a system-level understanding: (1) the physics of fluid flow and phase behavior in shales are poorly understood due to the dominant molecular interactions between minerals and fluids under confinement, and (2) the apparent lack of scale separation that prevents a reliable (closed) description of the physics at any single scale of observation. In this review, we focus on the latter issue and discuss scale translation, which in its broadest sense is transforming data or simulations from one spatiotemporal scale to another. While effective scale translation is not exclusive to shales, but all geologic porous media, the need for it is especially acute in shales given their high degree of heterogeneity. Classical theories like homogenization, while indispensable, fail when scales are not separated. Other methods, like numerical upscaling, scale-translate in only one direction: small to large, but not the reverse, called downscaling. However, the confluence of advances in three areas are bringing challenging problems such as shales within reach: increased computational power and scalable algorithms; high-resolution imaging and multi-modal data acquisition; and machine learning to process massive amounts of data. While these advances equip geoscientists with a wide array of experimental and computational tools, no individual tool can probe the entire gamut of heterogeneity in shales. Their effective use, therefore, requires an ability to bridge between various data types obtained at different scales. The aim of this review is to present a coherent account of computational and experimental methods that may be used to achieve just that, i.e., to perform scale translation. We provide a broader definition of scale translation, one that transcends classical homogenization and upscaling methods, but is consistent with them and accommodates notions like downscaling and data translation. After a brief introduction to homogenization, we review hybrid methods, numerical upscaling and its recent extensions, multiscale computing, high-resolution imaging, and machine learning. We place particular emphasis on multiscale computing and propose an algorithmic framework to bridge between the pore (micro) and Darcy (macro) scales. Throughout the paper, we draw comparisons between the various methods and highlight their (often hidden) similarities, differences, benefits, and pitfalls. We finally conclude with two case studies on shales that exemplify some of the methods presented.

58 GEOSCIENCES↗

Multichannel Analysis of Surface Waves Accelerated (MASWAccelerated): Software for efficient surface wave inversion using MPI and GPUs

Multichannel Analysis of Surface Waves (MASW) is a technique frequently used in geotechnical engineering and engineering geophysics to infer 1D layered models of seismic shear wave velocities in the top tens to hundreds of meters of the subsurface. We aim to accelerate MASW calculations by capitalizing on modern computer hardware available in the workstations of most engineers: multiple cores and graphics processing units (GPUs). We propose new parallel and GPU accelerated algorithms for computing 1D MASW inversion, and provide software implementations in C using Message Passing Interface (MPI) and CUDA. These algorithms take advantage of sparsity that arises in the problem, and the work balance between processes considers typical data trends. We compare our methods to an existing open source Matlab MASW tool. Our serial C implementation achieves a 2x speedup over the Matlab software, and we continue to see improvements by parallelizing the problem with MPI. Here we see nearly perfect strong and weak scaling for uniform data, and improve strong scaling for realistic data by repartitioning the problem to process mapping. By utilizing GPUs available on most modern workstations, we observe an additional 1.3x speedup over the serial C implementation on the first use of the method. We typically repeatedly evaluate theoretical dispersion curves as part of an optimization procedure, and on the GPU the kernel can be cached for faster reuse on later runs. We observe a 3.2x speedup on the cached GPU runs compared to the serial C runs. This work is the first open-source parallel or GPU-accelerated software tool for MASW imaging, and should enable geotechnical engineers to fully utilize all computer hardware at their disposal.

58 GEOSCIENCES↗

Physics-Informed and Data-Driven Prediction of Residual Stress in Three-Dimensional Machining

Efficient and reliable prediction of machining-induced residual stress (RS) is a key requirement for truly integrated computational materials engineering (ICME). Currently available process modeling approaches, including empirical, analytical, and numerical methodologies lack predictive power and require substantial calibration and validation data. Moreover, most model-based approaches consider only two-dimensional (2D) (i.e., orthogonal), cutting processes. Meanwhile, industrial processes such as milling, turning, and drilling are inherently three-dimensional (3D). The present work attempts to bridge the gap between 2D and 3D through careful consideration of the process physics, including geometric, kinematic, and size-effect constraints to realize robust prediction of how RS develops in 3D machining. Using a novel in-situ experimental technique and digital image correlation (DIC) to determine equivalent Hertzian contact widths, contact pressures, and friction coefficients, the proposed methodology leverages a discretized conversion algorithm that includes multi-pass shakedown effects. This paper presents a semi-analytical model to predict machining-induced RS in 3D turning operations, which are used representatively for 3D processes more generally. Rather than follow a ‘brute force’ 3D FEM approach or conduct countless experiments to train a purely data-driven machine learning algorithm, the proposed approach builds on previous 2D modeling work. Through careful consideration of the process physics, including complex geometry/kinematic considerations of 3D turning, the authors demonstrated an experimentally calibrated approach, as well as validation based on published RS data. Model predictions and previously published measurement data of RS depth profiles for turning of Inconel 718 were compared for a range of process parameters. Correlation between the proposed 3D model and validation data was found to be within the margin of experimental error for most conditions. The proposed model appears to capture the overall behavior of 3D RS depth profiles with acceptable accuracy, particularly the key metrics of near-surface stress, peak stress magnitude and location, as well as overall stress profile depth. This report presents a physics-informed, data-driven approach for efficient calibration of a 2D model for machining-induced RS through DIC analysis of in-situ characterized subsurface displacement fields.

42 ENGINEERING↗

Beagle 2: The Next Exobiology Mission to Mars

Beagle 2 is a 60 kg probe (with a 30 kg lander) developed in the United Kingdom for inclusion on the European Space Agency's 2003 Mars Express. Beagle 2 will deliver to the Martian surface a payload which consists of a high percentage of science instruments to landed spacecraft mass. Beagle 2 will be launched in June, 2003 with Mars Express on a Soyuz-Fregat rocket from the Baikonur Cosmodrome in Kazakhstan. Beagle 2 will land on Mars on December 26, 2003 in the Isidis Planitia basin (approximately 10 degrees N and 275 degrees W), a large sedimentary basin that overlies the boundary between ancient highlands and northern plains. Isidis Planitia, the third largest basin on Mars, which is possibly filled with sediment deposited at the bottom of long-standing lakes or seas, offers an ideal environment for preserving traces of life. Beagle 2 was developed to search for organic material and other volatiles on and below the surface of Mars in addition to the study of the inorganic chemistry and mineralogy. Beagle 2 will utilize a mechanical mole and grinder to obtain samples from below the surface, under rocks and inside rocks. A pair of stereo cameras will image the landing site along with a microscope for examination of surface and rock samples. Analyses will include both rock and soil samples at various wavelengths, X-ray spectrometer and Mossbauer spectrometer as well as a search for organics and other light element species (e.g. carbonates and water) and measurement of their isotopic compositions. Beagle 2 has as its focus the goal of establishing whether evidence for life existed in the past on Mars at the Isidis Planitia site or at least establishing if the conditions were ever suitable. Carbonates and organic components were first recognized as existing on Mars when they were found in the Martian meteorite Nakhla. Romanek et al showed the carbonates in ALH84001 were formed at low temperatures. McKay et al noted possible evidence of early life on Mars within the ALH84001 meteorite. Thomas-Keprta et al showed the magnetite biomarkers in ALH84001's carbonates are indistinguishable from those formed by magnetotactic bacteria found on Earth. Gibson et al showed there was significant evidence for liquid water and biogenic products present on Mars across a 3.9 billion year period. A mechanical arm (PAW) operates from the lander and is used for science operations along with sample acquisition. Instruments attached to the PAW include: stereo cameras, Mossbauer instrument, X-ray fluorescence instrument, microscope, environmental sensors, rock corer/grinder, a spoon, mirror, brushes, a mole attachment for acquisition of subsurface to depths of 1 to 2 meters and an illumination device. Each camera has 14 filters which have been optimized for mineralogy composition, dust and water vapor detection. The microscope's camera is designed for viewing the size and shape of dust particles, rock surfaces, microfossils, and characteristics of the samples prior to introduction into the gas analysis package (GAP). The camera has a resolution of 4 microns/pixel, features four color capability (red, green, blue and UV (ultraviolet) fluoresence), a depth of focus of 40 micrometers and translation stage of +3 millimeters. The heart of the Beagle 2's life detection package is the gas analysis package (GAP), which consists of a mass spectrometer with collectors at fixed masses for precise isotopic ratio measurements and voltage scanning for spectral analysis. Primary aim of the GAP is to search for the presence of bulk constituents, individual species, and isotopic fractionations for both extinct and extent life along with studying the low-temperature geochemistry of the hydrogen, carbon, nitrogen and oxygen components on Mars from both the surface and atmosphere. GAP is a magnetic sector mass spectrometer with the range of 1 to 140 amu which can be operated in both the static and dynamic modes. A triple Faraday collector array will be used for C, N and O ratios along with a double Faraday array for H/D. Pulse counting electron multiplier will be utilized for noble gases and selected organics. Anticipated detection limits are at the picomole level for operation in the static mode of operation and high precision isotopic measurements will be made in the dynamic mode. Sample processing and preparation system consists of reaction vessels along with references. Sample ovens capable of being heated are attached to the manifold for sample combustion. Surface, subsurface materials and interior rock specimens will be combusted in pure oxygen gas at various temperature intervals to release organic matter and volatiles. Combustion process will permit detection of all forms and all atoms of carbon present in the samples. A chemical processing system is capable of a variety of conversion reactions. Gases are manipulated either by cryogenic or chemical reactions and passed through the gas handling portion of the vacuum system. There are two modes of operation: quantitative analysis and precise isotopic measurements. Additional information is contained in the original extended abstract.

Gibson, Everett K., Jr.↗

NASAs Evolvable Mars Campaign: Mars Moons Robotic Precursor

Human exploration missions to the moons of Mars are being considered within NASA's Evolvable Mars Campaign (EMC) as an intermediate step for eventual human exploration and pioneering of the surface of Mars. A range of mission architectures is being evaluated in which human crews would explore one or both moons for as little as 14 days or for as long as 500 days with a variety of orbital and surface habitation and mobility options being considered. Relatively little is known about the orbital, surface, or subsurface characteristics of either moon. This makes them interesting but challenging destinations for human exploration missions during which crewmembers must be able to effectively conduct scientific exploration without being exposed to undue risks due to radiation, dust, micrometeoroids, or other hazards. A robotic precursor mission to one or both moons will be required to provide data necessary for the design and operation of subsequent human systems and for the identification and prioritization of scientific exploration objectives. This paper identifies and discusses considerations for the design of such a precursor mission based on current human mission architectures. Objectives of a Mars' moon precursor in support of human missions are expected to include: 1) identifying hazards on the surface and the orbital environment at up to 50-km distant retrograde orbits; 2) collecting data on physical characteristics for planning of detailed human proximity and surface operations; 3) performing remote sensing and in situ science investigations to refine and focus future human scientific activities; and 4) prospecting for in situ resource utilization. These precursor objectives can be met through a combination or remote sensing (orbital) and in-situ (surface) measurements. Analysis of spacecraft downlink signals using radio science techniques would measure the moon's mass, mass distribution, and gravity field, which will be necessary to enable trajectory planning. Laser altimetry would precisely measure the moon's shape and improve the accuracy of radio science measurements. A telescopic imaging camera would map the moon at submeter resolution and photograph selected areas of interest at subcentimeter resolution and a visible and near-infrared (0.4-3.0 mm) imaging spectrograph would produce a global map of mineral composition variations at a resolution of tens of meters and maps of selected areas of interest at meter resolution. Additional remote sensing capabilities could include a thermal infrared imager (heat flow, thermal inertia, and grain size distributions), a gamma-ray and neutron detector (atomic composition), a ground-penetrating radar (internal structure), and a magnetometer and Langmuir probe (magnetic properties and plasma field). Once on the surface of Phobos or Deimos, necessary instrumentation would include a penetrometer (regolith compressive strength), a motion-imagery camera (to observe the penetrometer tests before, during, and after contact), a dust-adhesion witness plate and camera (dust levitation), a microimager (dust particle sizes and shapes), and an alpha-proton-X-ray, X-ray fluorescence, Mossbauer, or Raman spectrometer (atomic and mineral composition of surface materials) and an optional temperature probe (regolith thermal properties). A variety of robotic mission design options to enable both orbital and surface measurements are being considered that include fully integrated and modular approaches. In-situ measurements from at least one surface location would be required, with additional measurement locations possible through use of multiple landers, through propulsive relocation of a single lander, or through electromechanical surface translation by a walking or hopping lander vehicle, which could also serve to evaluate such mobility capabilities for subsequent human missions. Preliminary orbital analysis suggests that remote sensing would likely be performed while in a distant retrograde orbit around the target moon. Mission design options to enable characterization of both Mars’ moons in a single mission are also being studied.

Gernhardt, Michael L.↗

Data for: Subsurface Interface Structure Controlling Local Electronic Properties of Epitaxial Graphene on SiC(0001)

Recently realized high-mobility semiconducting epitaxial graphene on silicon carbide, provided an important step towards integration of the graphene-based system into active components in post-silicon micro- and nano-electronics. However, the exact atomic-scale structure and the complex bonding configurations of the first epitaxial graphene carbon layer remain an open problem. Our recent report has shed new light on understanding this interface, where the external transverse electric field-dependent dynamic switching behavior of the Cbuffer-SiC bonds was observed. Here, using scanning tunneling microscopy and spectroscopy (STM and STS), we present the direct evidence of silicon (Si) vacancies at the interface and provide their distribution at the topmost reconstructed SiC(0001) layer. Experimental STM and density functional theory modeling data were used in the preparation of figures in a published article in the Journal of Physical Chemistry Letters. Files related to the figures and supplementary materials in the article are present in this dataset in .txt format.

Condensed matter imaging↗

Data for: Subsurface Interface Structure Controlling Local Electronic Properties of Epitaxial Graphene on SiC(0001)

Recently realized high-mobility semiconducting epitaxial graphene on silicon carbide, provided an important step towards integration of the graphene-based system into active components in post-silicon micro- and nano-electronics. However, the exact atomic-scale structure and the complex bonding configurations of the first epitaxial graphene carbon layer remain an open problem. Our recent report has shed new light on understanding this interface, where the external transverse electric field-dependent dynamic switching behavior of the Cbuffer-SiC bonds was observed. Here, using scanning tunneling microscopy and spectroscopy (STM and STS), we present the direct evidence of silicon (Si) vacancies at the interface and provide their distribution at the topmost reconstructed SiC(0001) layer. Experimental STM and density functional theory modeling data were used in the preparation of figures in a published article in the Journal of Physical Chemistry Letters. Files related to the figures and supplementary materials in the article are present in this dataset in .txt format.

Condensed matter imaging↗

Exploring for Superhot Geothermal Targets in Magmatic Settings: 2022 Field Campaign at Newberry Volcano

This paper presents preliminary results from a subset of work carried out as part of a multinational research project entitled DErisking Exploration for multiple geothermal Plays in magmatic ENvironments (DEEPEN), supported by the U.S. Department of Energy (DOE) and Geothermica, a joint effort by EU member states and associated countries. The DEEPEN project will develop a customized approach to exploration for supercritical and superhot geothermal plays in magmatic systems, which will be applied to two demonstration sites. This paper summarizes field activities carried out at the U.S. demonstration site, Newberry Volcano in central Oregon. The objective of this work effort is to refine the subsurface model of Newberry Volcano, with special focus on deeper zones including the magmatic plumbing system and other key geologic elements. New data collection included gravity and wideband magnetotelluric (MT) surveys, as well as reinstallation of a seismic network. The National Renewable Energy Laboratory (NREL) and Enthalpion Energy LLC (Enthalpion) worked with the Deschutes National Forest Fort Rock District to use a low ground disturbance method of MT deployment to collect MT data inside the caldera and other restricted areas inside the National Volcanic Monument. This opened these areas to geophysical exploration for the first time in decades. Sites along and adjacent to the south rim of the volcano constituted the primary survey objectives. A team from Lawrence Berkeley National Laboratory (LBNL), the U.S. Geological Survey (USGS), and AltaRock also began the process of reinstalling the seismic network from the AltaRock enhanced geothermal system (EGS) demonstration in anticipation of further development activities at the site. The data ingestion, reduction, and analysis phase of the project is ongoing. We are currently processing the MT and gravity data and are developing a new, highly GPU-accelerated, 3D joint MT and gravity inversion to better localize the south rim/south flank conductive target and better understand its relationship to deep heat, fluid sources, and surface extrusive features. Joint inversions, which have not yet been undertaken at Newberry, will allow us to obtain constraints on the geologic model that cannot be determined from each method in isolation, improving our ability to image key geologic features at depth.

geophysics↗

Remote Sensing of Icy Galilean Moon Surface and Atmospheric Composition Using Low Energy (1 eV-4 keV) Neutral Atom Imaging

We describe a low energy neutral atom imager suitable for composition measurements Europa and other icy Galilean moons in the Jovian magnetosphere. This instrument employs conversion surface technology and is sensitive to either neutrals converted to negative ions, neutrals converted to positive ions and the positive ions themselves depending on the power supply. On a mission such as the Jupiter Icy Moons Orbiter (JIMO), two back-to-back sensors would be flown with separate power supplies fitted to the neutral atom and iodneutral atom sides. This will allow both remote imaging of 1 eV < E < 4 keV neutrals from icy moon surfaces and atmospheres, and in situ measurements of ions at similar energies in the moon ionospheres and Jovian magnetospheric plasma. The instrument provides composition measurements of the neutrals and ions that enter the spectrometer with a mass resolution dependent on the time-of-flight subsystem and capable of resolving molecules. The lower energy neutrals, up to tens of eV, arise from atoms and molecules sputtered off the moon surfaces and out of the moon atmospheres by impacts of more energetic (keV to MeV) ions from the magnetosphere. Direct Simulation Monte Carlo (DSMC) models are used to convert measured neutral abundances to compositional distributions of primary and trace species in the sputtered surfaces and atmospheres. The escaping neutrals can also be detected as ions after photo- or plasma-ionization and pickup. Higher energy, keV neutrals come from charge exchange of magnetospheric ions in the moon atmospheres and provide information on atmospheric structure. At the jovicentric orbits of the icy moons the presence of toroidal gas clouds, as detected at Europa's orbit, provide M e r opportunities to analyze both the composition of neutrals and ions originating from the moon surfaces, and the characteristics of magnetospheric ions interacting with neutral cloud material. Charge exchange of low energy ions near the moons, and directional distributions of the resultant neutrals, allow indirect global mapping of magnetic field structures around the moons. Temporal variation of the magnetic structures can be linked to induced magnetic fields associated with subsurface oceans.

Collier, M. R.↗

Electrical Resistivity Tomography based monitoring of stress perturbations to optimize placement of high-precision strain meters

The Center for Understanding Subsurface Signals and Permeability is a new U.S. Department of Energy Earthshot Center focused on understanding and predicting the long-term evolution of permeability in enhanced geothermal systems. The center will use a highly instrumented testbed within the Sanford Underground Research Facility to conduct field scale experiments that elucidate and test capabilities to simulate geochemical-geomechanical interactions and permeability evolution. Here we demonstrate initial developments using previously collected electrical resistivity tomography (ERT) monitoring data with high-performance multi-physics modelling advancements to inform the optimal location of two new monitoring boreholes. Specifically, ERT monitoring data collected during shear stimulation testing shows marked responses to changes in stress during borehole pressurization. We demonstrate how the same response is being simulated, ultimately to train a machine-learning algorithm to estimate rock properties and enable enhanced prediction of stress and strain responses anticipated during future testing campaigns.

Stress, EGS, CUSSP, 3D Electrical Imaging↗

Fracture Intersections under Stress: Laboratory Data and Code [Data set]

The connectivity of natural and induced fractures governs the injection and withdrawal of fluids from subsurface reservoirs. Connectivity depends on intersections that control how fluids mix and move through the entire system. Here, we present data sets from 3D X-ray microscopy measurements of simple fracture networks under stress. 3D printing was used to create prismatic blocks that formed fracture networks composed of 2 orthogonal fractures. The network orientation was either "x" or "+" relative to an applied vertical stress. 3D data sets were collected for normal loads of 25, 100 and 200 Newtons for samples with fracture surfaces with either correlated or uncorrelated asperity distributions. The file contains data from the 12 samples analyzed along with an example code used to extract the intersection geometry. Additional experimental details can be found in the manuscript "Geologic Stress Modulates Fluid Mixing at Fracture Intersections" (10.1038/s43247-026-03525-9)and supplemental information to appear in Communications Earth & Environment in 2026.

02 PETROLEUM↗

Mapping Rock and Soil Units in the MPF IMP SuperPan Using a Kohonen Self Organizing Map

The 1997 Mars Pathfinder mission provided information on a site in the Ares Vallis floodplain. Initial analysis of multispectral data from the Imager for Mars Pathfinder (IMP) indicated the presence of only a single rock type, the 'gray rock' spectral class and various coated variants thereof (e.g., 'maroon rock'). Continued analysis of the IMP 'SuperPan' mosaic has confirmed multiple examples of a second 'black rock' spectral class existing as small cobbles in the near field and as boulders in the far field. These results are consistent with recent analysis of MGS Thermal Emission Spectrometer (TES) data which indicates that there is likely a mix of both 'Surface Type 1' (ST1) and 'Surface Type 2' (ST2) spectral classes at the MPF landing site. Nominally, the black rock spectral class would correspond to ST1 (basalts) and 'gray rock' would correspond to ST2 (andesites). Orbital remote sensing has also revealed the pervasive presence of layering on Mars. Recently it was suggested that there are extensive outcrops of the black rock spectral class in the SuperPan far field on the flanks of the Twin Peaks and on the rim of Big Crater. These authors suggested that these exposures represented outcrops of black rock from beneath a surficial, flood deposited layer. In this work, we have reexamined the MPF IMP SuperPan mosaic using an artificial neural network self organizing map (SOM) processing architecture in order to classify the distribution of spectral classes within the SuperPan. In this paper, we present initial results from that work and draw specific attention to a subset of the identified spectral classes in order to address questions relating to whether there are extensive exposures of black rock in the IMP far field, what other materials might be exposed in the far field, and what evidence there is for subsurface layering at the MPF landing site.

Farrand, W.↗

InversionNet3D: Efficient and Scalable Learning for 3-D Full-Waveform Inversion

Seismic full-waveform inversion (FWI) techniques aim to find a high-resolution subsurface geophysical model provided with waveform data. Some recent effort in data-driven FWI has shown some encouraging results in obtaining 2-D velocity maps. However, due to high computational complexity and large memory consumption, the reconstruction of 3-D high-resolution velocity maps via deep networks is still a great challenge. Here, in this article, we present InversionNet3D (InvNet3D), an efficient and scalable encoder–decoder network for 3-D FWI. The proposed method employs group convolution in the encoder to establish an effective hierarchy for learning information from multiple sources while cutting down unnecessary parameters and operations at the same time. The introduction of invertible layers further reduces the memory consumption of intermediate features during training and, thus, enables the development of deeper networks with more layers and higher capacity as required by different application scenarios. Experiments on the 3-D Kimberlina dataset demonstrate that InvNet3D achieves state-of-the-art reconstruction performance with lower computational cost and lower memory footprint compared to the baseline.

58 GEOSCIENCES↗

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.↗

Unveiling Clues from Spacecraft Missions to Comets and Asteroids through Impact Experiments

The Deep Impact Spacecraft mission was the first to boldly face the challenge of impacting the surface of a comet, 9P/Tempel 1, to investigate surface and subsurface 'pristine' materials. The Stardust mission to Comet 81P/Wild 2 brought back an exciting surprise: shocked minerals which were likely altered during the comet's lifetime. Signatures of shock in meteorites also suggest that the violent past of the solar system has left our small bodies with signatures of impacts and collisions. These results have led to the question: How have impacts affected the evolutionary path taken by comets and asteroids, and what signatures can be observed? A future planetary mission to a near-Earth asteroid is proposing to take the next steps toward understanding small bodies through impacts. The mission would combine an ESA led AIM (Asteroid Impact Mission) with a JHU/APL led DART (Double Asteroid Redirect Mission) spacecraft to rendezvous with binary near-Earth asteroid 65803 Didymus (1996 G2). DART would impact the smaller asteroid, 'Didymoon' while AIM would characterize the impact and the larger Didymus asteroid. With these missions in mind, a suite of experiments have been conducted at the Experimental Impact Laboratory (EIL) at NASA Johnson Space Center to investigate the effects that collisions may have on comets and asteroids. With the new capability of the vertical gun to cool targets in the chamber through the use of a cold jacket fed by liquid nitrogen, the effects of target temperature have been the focus of recent studies. Mg-rich forsterite and enstatite (orthopyroxene), diopside (monoclinic pyroxene) and magnesite (Mg-rich carbonate) were impacted. Target temperatures ranged from 25 deg to -100 deg, monitored by connecting thermocouples to the target container. Impacted targets were analyzed with a Fourier Transform Infrared Spectrometer (FTIR) and Transmission Electron Microscope (TEM). Here we present the evidence for impact-induced shock in the minerals through both spectra and TEM imaging and compare with unshocked samples.

Lederer, Susan M.↗

Autonomous Ocean World Exploration: Advancement of a Virtual Testbed

The search for life (extinct or extant) and potentially habitable bodies in our solar system and beyond is one of the 12 priority science questions outlined in the National Acadamies’ 2022 decadal survey [5]. Extraterrestrial destinations containing liquid water present an opportunity to search for life as we know it, and in recent years an increasing number of such locations have been discovered within our solar system. Several Jovian moons—Europa, Ganymede, and Callisto [10]—and the Saturnian moons Enceladus [8] and Titan [9] are known or suspected to harbor massive subsurface oceans. Of these "ocean worlds", Europa is the focus of at least one planned NASA orbiter mission, Europa Clipper [4], and an early lander mission concept, the Europa Lander [2, 3]. Whereas most robotic missions to the Moon and Mars (e.g. orbiters, rovers, landers) to date have had ground controllers on Earth tightly involved in mission operations, missions to more distant worlds will require a high degree of onboard autonomy due to long communication lags and blackouts, harsh environments (radiation, cold), and more limited battery and hardware life. The past decade has seen great advances in both AI technologies and computing scalability and performance that offer promising solutions for spacecraft autonomy and motivate the software system and research programs described in this paper. The Ocean Worlds Autonomy Testbed for Exploration, Research, and Simulation (OceanWATERS) [1], which has been in development at the NASA Ames Research Center since 2018, is a virtual environment for testing lander autonomy solutions. It is built on the Robot Operating System (ROS), runs on consumer-grade Linux workstations, and was released as open source in 2020. OceanWATERS provides a physical and visual simulation of a prototypical lander in a Europa-like environment (Figure 1). The lander was modeled after requirements and specifications made in JPL’s Europa Lander Study of 2016 [3]. Simulated lander systems include stereo cameras and spotlights mounted on an antenna mast that pans and tilts, a 6 degrees of freedom (DoF) robotic arm with a force-torque sensor and two interchangeable end effectors, and a battery pack power system. The environment consists of multiple terrain models including a highly detailed model sourced from the FROST dataset [11], simulation of surrounding planetary bodies based on an ephemeris model, and lighting from the sun with associated surface illumination, reflectance, and shadows. Operations supported by OceanWATERS include panoramic and directed imaging of the environment and lander workspace, Cartesian and joint-level arm commanding, grinding of the terrain surface (e.g. digging a trench), and scooping of ground material (Figure 2) which can be discarded or collected as science samples in a receptacle that can be emptied (science operations themselves are not simulated). These operations are realized as ROS Actions and are complimented by a wide selection of telemetry that is continually produced by each lander subsystem. The power system model is driven by the open-source Generic Software Architecture for Prognostics (GSAP) [11] that predicts the battery’s remaining useful life and other characteristics. As a testbed for high-level autonomy, OceanWATERS provides an execution framework based on PLEXIL [12], an open-source plan specification language and execution engine developed largely at Ames. NASA's initial development of OceanWATERS, as well the Ocean Worlds Lander Autonomy Testbed (OWLAT) [6], a complimentary physical testbed developed at JPL, was the first step in a plan for realizing candidate onboard autonomy solutions for such planetary landers. In 2020 NASA solicited applications for its Autonomous Robotics Research for Ocean Worlds (ARROW) program, and in 2021 the similar Concepts for Ocean worlds Life Detection Technology (COLDTech) program. Collectively six research teams, based in universities and companies across the United States, were awarded grants to develop and demonstrate autonomy solutions on OceanWATERS and OWLAT. These 1–2-year projects have now finished or are nearing completion, and a wide variety of autonomy challenges in ocean world surface missions were addressed. Prototyped and demonstrated solutions have included autonomous discovery, response and adaptation to system faults and unexpected environmental events, world model synthesis through perception, plan synthesis using learned models, methods to optimize sample target selection and prioritize science data transmission, extension of PLEXIL for stochastic decision-making, and an integration of a model of JPL’s mission-ready COLDArm [7]. Technologies used in these projects include many forms of machine learning, causal reasoning, automated planning, Markov decision processes, formal methods, and other advanced techniques. A more detailed summary of the ARROW and COLDTech projects is given herein. OceanWATERS has had significant enhancements since its open-source release in 2020. Many of its new features were driven or shaped by feedback from the ARROW and COLDTech teams and requirements of their projects. In support of enabling autonomous adaptation to spacecraft faults (a specific capability solicited by both programs), a fault injection and detection framework was developed that supports a wide and growing range of fault types such as locked joints, image loss, and battery failures. The power system model was completed and integrated into the simulator, starting as a single-cell battery model and later upgraded to a multi-cell model with associated faults such as cell disconnection. Arm/terrain interaction was improved by adding a force-torque sensor and associated faults, and an analytic dig force model based on the Balovnev bucket force equations. Environment fidelity was increased by modeling terrain deformation resulting from digging and scooping; visual improvements were made in textures, lighting, and shadows. To facilitate interoperation with OWLAT, a unified command and telemetry interface between the testbeds was developed at the ROS level, along with a PLEXIL interface. The number of lander operations was greatly expanded (e.g. with Cartesian-based arm and antenna movement), and a framework was designed for users to build their own lander actions. A GUI for PLEXIL plan selection was created (Figure 3), and an expansive set of plans were added, such as those that illustrate patterns for fault handling. This paper provides a self-contained high-level description of OceanWATERS, focusing on more detailed coverage of the aforementioned enhancements. It provides a high-level summary of the projects undertaken by participants in the ARROW and COLDTech programs and how these efforts have helped shape OceanWATERS. Finally, potential future work and directions for the testbed are listed, as likely informed by the recent planetary science decadal survey [5].

K Michael Dalal↗