Search NASA⌕ Search

SEARCH · Search NASA

Results for “subsurface data”

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 361 records · Page 20

FracML: A Machine Learning Based Tool to Quantify Reservoir Scale Fracture Network for CO2 Storage

Poster on “FRACML: A Machine Learning Based Tool to Quantify Reservoir Scale Fracture Network for CO2 Storage” for the CCUS 2025 conference held in Houston, Texas March 3-5, 2025. The accurate characterization of subsurface fracture networks is essential for the secure operation of carbon capture, utilization, and storage (CCUS) projects. A thorough understanding of the spatial distribution of subsurface faults and fractures is crucial for predicting CO2 plume evolution and minimizing risks such as potential leakage into overlying formations or induced seismicity. In this context, robust fracture network quantification plays a pivotal role in reservoir management, providing the data necessary to fine-tune operational parameters, and ensure the environmental and economic viability of CCUS projects. As part of the U.S. Department of Energy’s SMART (Science-informed Machine Learning for Accelerating Real-time Decisions in Subsurface Applications) initiative, we focused on the development and application of a machine learning-based tool (FRACML) designed to quantify and map fracture networks using real-world (non-synthetic) data from an active CO2 injection site. Our objective is to demonstrate the utility of this tool in improving operational efficiency and safety across CCUS sites.

artifical intelligence / machine learning (AI/ML)↗

Thermo-hydro-mechanical analysis of subsurface ice-based thermal energy storage

Ice-based thermal energy storage systems are widely utilized for cooling and managing peak electrical demand globally, offering daily or weekly storage capabilities for both individual homes and larger office buildings. However, scaling these systems for district-level cooling or integrating them with renewable energy sources presents challenges, especially in accommodating larger volumes and addressing seasonal storage requirements in densely populated urban areas. This paper proposes a novel solution by evaluating subsurface ice-based thermal energy storage, in which the underground is subjected to seasonal freeze/thaw cycles. However, these cycles may influence ground behavior, affecting pore pressure and inducing ground movement. To systematically investigate these challenges, we enhance the TOUGH-FLAC simulator by integrating water/ice phase change capabilities and updating the effective stress–strain constitutive relation. Both modifications are validated against analytical solutions or experimental data. Through numerical simulations spanning a decade with ten seasonal freeze/thaw cycles, we evaluate the performance and long-term stability of a generic subsurface ice-based thermal energy storage system, considering factors such as ground permeability, freezing pipe spacing, freeze/thaw damage, and glycol solution temperature. The simulations indicate that ice formation induces pore pressure variations that drive seasonal surface heave and settlement, controlled by ground permeability, pipe spacing, and glycol solution temperature, along with tensile and localized shear deformation around freeze pipes. This highlights the need for accurate ground property characterization and geomechanical analysis for subsurface ice-based thermal energy storage.

Ground freezing↗

A cool anomaly off northern California - An investigation using IR imagery and in situ data

The OPTOMA (Ocean Prediction Through Observation, Modeling and Analysis) program is developing an ocean descriptive-predictive system for four-dimensional data assimilation. It is presently concerned with the mesoscale variability in the California Current System (CCS). The present paper has the objective to assess the relationship of surface temperature structure to subsurface temperature structure and flow fields. Surface temperature anomalies are related to the mesoscale horizontal advective and subsurface mass fields. The cool anomaly off northern California in summer 1982 is discussed, taking into account the temperature structure along transects, T-S variations and inferred water masses, variation of horizontal temperature patterns with depth, and synoptic information from IR data.

Rienecker, M. M.↗

The Project “Analogs for VENus’s GEologically Recent Surfaces” (AVENGERS): Unveiling Active Volcanism on Venus

The recently selected missions to Venus opened a new era for the exploration of the Earth’s hellish twin sister. One of the key questions that the selected missions to Venus will address is whether the planet is volcanically active. The identification of ongoing volcanic activity on Venus can indeed unveil important clues on the geologic evolution of the planet, also helping us to better understand how volcanic and atmospheric processes work on Earth, especially those related to the climate change produced by greenhouse gases. To fully contextualize the higher resolution datasets to be provided by the future missions to Venus, we need to be able to exploit the wealth of data to be produced in the best possible way. Different types of instruments will indeed perform unprecedent orbiting and in-situ investigations of the surface and subsurface of Venus. To this regard, the Project “Analogs for VENus’s GEologically Recent Surfaces” (AVENGERS) aims to build a comprehensive database of Terrestrial active volcanoes which may act as suitable analogs for identifying ongoing volcanic activity on Venus. To do this, the project AVENGERS will combine the implementation of terrestrial orbiter radar data analysis techniques (i.e., radar interferometry) with infrared and laboratory analyses of lava flow samples from Terrestrial active volcanoes. As a first suitable analog of the Project AVENGERS, we selected Mount Etna in Sicily for a direct comparison with potentially active volcanic structures of Venus, and with Idunn Mons, one of the most promising possibly active volcanoes of Venus. Given its composite nature, Mount Etna offers the opportunity to identify and study multiple types of eruptive products (both effusive and explosive) from the analysis of high-resolution radar images. The frequent eruptions characterizing this volcano will also allow to test techniques of radar interferometry. Moreover, its ease of access facilitates geologic field trips and in-situ investigations with related sample retrieving. Ongoing infrared laboratory analysis on several altered and unaltered samples retrieved by Mount Etna will help to distinguish spectral signatures of fresh vs weathered volcanic deposits.

P. D’Incecco↗

Subsurface Emission Effects in AMSR-E Measurements: Implications for Land Surface Microwave Emissivity Retrieval

An analysis of land surface microwave emission time series shows that the characteristic diurnal signature associated with subsurface emission in sandy deserts carry over to arid and semi-arid region worldwide. Prior work found that diurnal variation of Special Sensor Microwave/Imager (SSM/I) brightness temperatures in deserts was small relative to International Satellite Cloud Climatology Project land surface temperature (LST) variation and that the difference varied with surface type and was largest in sand sea regions. Here we find more widespread subsurface emission effects in Advanced Microwave Scanning Radiometer-EOS (AMSR-E) measurements. The AMSR-E orbit has equator crossing times near 01:30 and 13 :30 local time, resulting in sampling when near-surface temperature gradients are likely to be large and amplifying the influence of emission depth on effective emitting temperature relative to other factors. AMSR-E measurements are also temporally coincident with Moderate Resolution Imaging Spectroradiometer (MODIS) LST measurements, eliminating time lag as a source of LST uncertainty and reducing LST errors due to undetected clouds. This paper presents monthly global emissivity and emission depth index retrievals for 2003 at 11, 19, 37, and 89 GHz from AMSR-E, MODIS, and SSM/I time series data. Retrieval model fit error, stability, self-consistency, and land surface modeling results provide evidence for the validity of the subsurface emission hypothesis and the retrieval approach. An analysis of emission depth index, emissivity, precipitation, and vegetation index seasonal trends in northern and southern Africa suggests that changes in the emission depth index may be tied to changes in land surface moisture and vegetation conditions

Galantowicz, John F.↗

Estimates of the zonal slope and seasonal transport of the Atlantic North Equatorial Countercurrent

Data from six inverted echo sounder moorings and the Geosat satellite altimeter are used to examine the seasonal variability of sea surface elevation. Monthly sea level maps are constructed using a contemporaneous subsurface temperature survey to provide a reference sea level field. The maps are then used to describe the origin and structure of the western tropical Atlantic North Equatorial Countercurrent (NECC) during a two-year period beginning in November 1987. The data reveal a zonal current which is confined between 3 deg N and 9 deg N with a typical width of 300 km. The NECC flows strongly eastward during November and December 1986 and May 1987 through January 1988. The reappearance of the current is then delayed until August, but the current flows strongly from August until the end of the record in October 1988. Volume transport is estimated for the two-year period from surface elevation by approximating the vertical structure of the ocean as a two-layer fluid. It is found that the NECC has a maximum transport of 40 x 10 to the 6th cu m/s at 38 deg W.

Carton, James A.↗

Interference patterns of a horizontal electric dipole over layered dielectric media.

Interference patterns for electromagnetic fields due to a subsurface reflector below a layered lossy dielectric are calculated with the geometrical optics approximation for use in interpreting data to be collected on the moon by Apollo 17 as well as data currently being obtained on terrestrial glaciers. The radiating antenna lies on the surface. All six field components are calculated and studied. For the endfire solutions, the peak of the first reflected wave is found to be different from that of the broadside ones. To facilitate a physical discussion, we plotted the radiation patterns due to the antenna on the surface.

Tsang, L.↗

Elemental analysis of a comet nucleus by passive gamma ray spectrometry from a penetrator

The Comet Rendezvous Asteroid Flyby (CRAF) spacecraft, to be launched in 1991, was designed to study the physical and chemical properties of cometary and asteroid bodies. It is proposed that these properties can be determined by utilizing a penetrator experiment delivery system aboard the CRAF, which would deliver a passive gamma ray spectrometer to the comet and determine the composition of subsurface materials. Results of calculations for a model comet (50 percent CI carbonaceous chrondrite, 50 percent ices) provide data on the expected gamma ray fluxes, the minimum detectable limits, and the detection uncertainties for 35 isotopes. These results confirm that this technique can determine the composition of both the rocky and ice components, and the rock-to-ice ratio in the nucleus. In addition, useful information regarding the formation of the solar system is provided.

Evans, L. G.↗

The Evolution of Oblique Impact Flow Fields Using Maxwell's Z Model

Oblique impacts are the norm rather than the exception for impact craters on planetary surfaces. This work focuses on the excavation of experimental oblique impact craters using the NASA Ames Vertical Gun Range (AVGR). Three-dimensional particle image velocimetry (3D PIV) is used to obtain quantitative data on ejection positions, three-dimensional velocities and angles. These data are then used to test the applicability and limitations of Maxwell's Z Model in representing the subsurface evolution of the excavation-stage flow-field center during vertical and oblique impacts.

Anderson, J. L. B.↗

High-Power Radar Sounders for the Investigation of Jupiter Icy Moons

The high power and high data rate capability made available by a Prometheus class spacecraft could significantly enhance our ability to probe the subsurface of the planets/moons and asteroid/comets. The main technology development driver for our radar is the proposed Jupiter Icy Moon Orbiter (or JIMO) mission due to its harsh radiation environment. We plan to develop a dual-band radar at 5 and 50 MHz in response to the two major science requirements identified by the JIMO Science Definition Team: studying the near subsurface (less than 2 km) at high resolution and detection of the ice/ocean interface for Europa (depth up to 30 km). The 50-MHz band is necessary to provide high spatial resolution (footprint and depth) as required by the JIMO mission science requirements as currently defined. Our preliminary assessment indicates that the 50-MHz system is not required to be as high-power as the 5-MHz system since it will be more limited by the surface clutter than the Jupiter or galactic background noise. The low frequency band (e.g. 5 MHz), which is the focus of this effort, would be necessary to mitigate the performance risks posed by the unknown subsurface structure both in terms of unknown attenuation due to volumetric scattering and also the detection of the interface through the attenuative transition region at the ice/ocean interface. Additionally, the 5-MHz band is less affected by the surface roughness that can cause loss of coherence and clutter noise. However, since the Signal-to-Noise-Ratio (SNR) of the 5-MHz radar band is reduced due to Jupiter noise when operating in the Jupiter side of the moon, it is necessary to increase the radiated power. Our challenge is to design a high-power HF radar that can hnction in Jupiter's high radiation environment, yet be able to fit into spacecraft resource constraints such as mass and thermal limits. Our effort to develop the JIMO radar sounder will rely on our team's experience with planetary radar sounder design gained during our participation in the MARSIS radar sounder implementation.

Safaeinili, A.↗

Warming of the Global Ocean: Spatial Structure and Water-Mass Trends

This study investigates the multidecadal warming and interannual-to-decadal heat content changes in the upper ocean (0-700 m), focusing on vertical and horizontal patterns of variability. These results support a nearly monotonic warming over much of the World Ocean, with a shift toward Southern Hemisphere warming during the well-observed past decade. This is based on objectively analyzed gridded observational datasets and on a modeled state estimate. Besides the surface warming, a warming climate also has a subsurface effect manifesting as a strong deepening of the midthermocline isopycnals, which can be diagnosed directly from hydrographic data. This deepening appears to be a result of heat entering via subduction and spreading laterally from the high-latitude ventilation regions of subtropical mode waters. The basin-average multidecadal warming mainly expands the subtropical mode water volume, with weak changes in the temperature-salinity (u-S) relationship (known as ''spice'' variability). However, the spice contribution to the heat content can be locally large, for example in Southern Hemisphere. Multidecadal isopycnal sinking has been strongest over the southern basins and weaker elsewhere with the exception of the Gulf Stream/North Atlantic Current/subtropical recirculation gyre. At interannual to decadal time scales, wind-driven sinking and shoaling of density surfaces still dominate ocean heat content changes, while the contribution from temperature changes along density surfaces tends to decrease as time scales shorten.

Global Ocean↗

3D Modeling of Solar-Type Stars to Characterize Stellar Jitter

Detection of Earth-mass planets requires measurements of radial velocity with extreme precision. To capture the tiny disturbances caused by a planet's motion, it is necessary to understand and characterize the host star’s turbulent dynamics in order to apply proper filtering to the observational data. We take advantage of current computational and technological capabilities to develop 3D realistic models of the stellar subsurface convection and atmospheres and thereby estimate the photospheric jitter. We have identified an initial set of target stars, obtained initial conditions using the MESA code, and obtained initial 3D radiative models of the stellar surfaces and atmospheres with a spatial resolution of 50km. We present initial 3D radiative hydrodynamic model results of the planet-hosting star HD209458.

SMD↗

Role of Oxygen Fugacity on the Melting Properties of Enstatite Chondrites and Implications for Mercury’s Magmatic Evolution

Mercury, the innermost terrestrial planet, is the most reduced planet in our solar system. Insights from MESSENGER mission data revealed the presence of several distinct geochemical terranes, evidence of complex magmatic processes and collisional processes exposing subsurface materials. Surface chemical analysis indicated elevated S (~2-3 wt%) and low FeO (~1.5 wt%) concentrations. The high sulfur concentration indicates reduced conditions with an average of 5.4 log units below the Iron-Wüstite (IW) oxygen fugacity ( O 2 ) buffer (IW-5.4). Surface compositions also show a range in redox conditions during mantle melting and eruption, with inferred log f O 2 , ranging from IW-6.5 to IW-3.5. The effects of oxygen fugacity on magmatic differentiation are poorly constrained, despite their significance in our understanding of mantle-crust differentiation in the solar system. Enstatite High-Fe (EH) chondrites are very reduced undifferentiated meteorites with elevated concentrations of Fe and volatiles like S, Cl, Na, and K as compared to other chondrites. These characteristics suggest that EH chondrites are a potential analog for Mercury’s building blocks. However, the comparison of Mercury’s surface composition with melting products of EH chondrites is necessary to determine whether Mercury surface materials may be derived from EH chondrite- like materials. Here, we investigate the role of f O 2 on EH melting properties and its implications for Mercury’s accretion and differentiation.

Mercury↗

Satellite microwave observations of the Utah Great Salt Lake Desert

Microwave data acquired over the Great Salt Lake Desert area by sensors aboard Skylab and Nimbus 5 indicate that the microwave emission and backscatter were strongly influenced by contributions from subsurface layers of sediment saturated with brine. This phenomenon was observed by Skylab's S-194 radiometer operating at 1.4 GHz, S-193 RADSCAT (Radiometer-Scatterometer) operating at 13.9 GHz, and the Nimbus 5 ESMR (Electrically Scanning Microwave Radiometer) operating at 19.35 GHz. The availability of ESMR data over an 18-month period allowed an investigation of temporal variations.

Ulaby, F. T.↗

Satellite ocean color measurements

The application of pattern recognition to ocean color data analysis is considered. Due to weight, cost, and data transmission rate limitations, any mapping type remote sensor of ocean color must necessarily collect light from the sea in a finite number of channels. The optical properties of the sea are discussed together with the remote sensing of ocean color, an optical model of natural water, the microscopic optical model, the macroscopic optical model, multiple scattering theory, measurements of subsurface oceanographic parameters, measurements of bulk absorption and scattering properties, measurements of the up- and down-welling light field, and the techniques for ocean color data analysis.

Mccluney, W. R.↗

Machine Learning Applications in Analyzing the Role of Shale Barriers and Baffles for CO2 Storage

This study uses machine learning to analyze microseismic data from the Illinois Basin Decatur Project (IBDP) and quantify CO₂ plume extents. By leveraging well logs, microseismic records, and CO₂ injection metrics, the research predicts subsurface CO₂ plume dynamics. Findings show vertical clustering of microseismic events near the injection well, with CO₂ periodically breaching barriers due to buoyancy. K-Means clustering performed best, achieving the highest Silhouette Score and lowest Davies-Bouldin Index. This capability is crucial for real-time monitoring and management of CO₂ sequestration sites, validated against physical models and IBDP data, reinforcing CO₂ geological sequestration's viability and enhancing management tools.

Carr, Timothy↗

A Rover Deployed Ground Penetrating Radar on Mars

Radar is a fundamental tool capable of addressing a variety of geological problems on Mars via collection of data suitable for interpreting variations in surface morphology and reflectivity. Surface-deployed ground penetrating radar (GPR) can help further constrain the geology and structure of the near surface of Mars by directly measuring the range and character of in situ radar properties. In recognition of this potential, a miniaturized, easily modified GPR is being developed for possible deployment on a future Mars rover and will enable definition of radar stratigraphy at high spatial resolution to depths of 10-20 meters. Ongoing development of a Mars impulse GPR with industry partners at Geophysical Survey Systems, Inc., focuses on design and testing of a prototype transducer array (with both high frequency bistatic and low frequency monostatic components) in parallel with fabrication of a low power, mass, and volume control unit. The operational depth of 10-20 meters is geared towards definition of stratigraphy, subsurface blocks, and structure at the decimeter to meter scale that is critical for establishing the geologic setting of the rover. GPR data can also be used to infer the degree of any post-depositional pedogenic alteration or weathering that has subsequently taken place, thereby enabling assessment of pristine versus secondary morphology at the landing site. As is the case for most remote sensing instruments, a GPR may not detect water unambiguously. Nevertheless, any local, near-surface occurrence of liquid water will lead to large, easily detected dielectric contrasts. Moreover, definition of stratigraphy and setting will help in evaluating the history of aqueous activity and where any water might occur and be accessible. Most importantly perhaps, GPR can provide critical context for other rover and orbital instruments/data sets. Hence, GPR deployment along well positioned transects in the vicinity of a lander should enable 3-D mapping of stratigraphy and could serve to guide direct subsurface sampling. Additional information is contained in the original extended abstract.

Grant, J. A.↗

Exploration of the Moon with Remote Sensing, Ground-Penetrating Radar, and the Regolith-Evolved Gas Analyzer (REGA)

There are two important reasons to explore the Moon. First, we would like to know more about the Moon itself: its history, its geology, its chemistry, and its diversity. Second, we would like to apply this knowledge to a useful purpose. namely finding and using lunar resources. As a result of the recent Clementine and Lunar Prospector missions, we now have global data on the regional surface mineralogy of the Moon, and we have good reason to believe that water exists in the lunar polar regions. However, there is still very little information about the subsurface. If we wish to go to the lunar polar regions to extract water, or if we wish to go anywhere else on the Moon and extract (or learn) anything at all, we need information in three dimensions an understanding of what lies below the surface, both shallow and deep. The terrestrial mining industry provides an example of the logical steps that lead to an understanding of where resources are located and their economic significance. Surface maps are examined to determine likely locations for detailed study. Geochemical soil sample surveys, using broad or narrow grid patterns, are then used to gather additional data. Next, a detailed surface map is developed for a selected area, along with an interpretation of the subsurface structure that would give rise to the observed features. After that, further sampling and geophysical exploration are used to validate and refine the original interpretation, as well as to make further exploration/ mining decisions. Integrating remotely sensed, geophysical, and sample datasets gives the maximum likelihood of a correct interpretation of the subsurface geology and surface morphology. Apollo-era geophysical and automated sampling experiments sought to look beyond the upper few microns of the lunar surface. These experiments, including ground-penetrating radar and spectrometry, proved the usefulness of these methods for determining the best sites for lunar bases and lunar mining operations.

Cooper, B. L.↗