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 451 records · Page 25

Synthesis and Assimilation Systems - Essential Adjuncts to the Global Ocean Observing System

Ocean assimilation systems synthesize diverse in situ and satellite data streams into four-dimensional state estimates by combining the various observations with the model. Assimilation is particularly important for the ocean where subsurface observations, even today, are sparse and intermittent compared with the scales needed to represent ocean variability and where satellites only sense the surface. Developments in assimilation and in the observing system have advanced our understanding and prediction of ocean variations at mesoscale and climate scales. Use of these systems for assessing the observing system helps identify the strengths of each observation type. Results indicate that the ocean remains under-sampled and that further improvements in the observing system are needed. Prospects for future advances lie in improved models and better estimates of error statistics for both models and observations. Future developments will be increasingly towards consistent analyses across components of the Earth system. However, even today ocean synthesis and assimilation systems are providing products that are useful for many applications and should be considered an integral part of the global ocean observing and information system.

Rienecker, Michele M.↗

Determining stress states using dike swarms: The Lauma Dorsa example

Initial examination of the Magellan coverage of Venus has revealed between 150 and 300 large, radially lineated landforms distributed across the planet's surface. Where the lineaments have been examined in detail, the majority fail to exhibit signatures indicative of relief at or above the resolution of the radar; however, when the sense of topographic relief may be ascertained, the lineaments commonly appear as fissures or flat-floored trenches interpreted as graben. Individual lineaments can display graben, fissure, and zero relief behavior along their length, suggesting either that these differences are a function of the resolution of the radar, or that the morphological distinctions are real but somehow genetically linked. In many instances, radial lineaments exhibiting these characteristics are directly associated with surface volcanism, including flanking and terminal flows, superimposed shield domes and pit chains, and central, calderalike topographic lows. These observable characteristics, as well as theoretical studies and comparison with similar terrestrial features, have led to the working hypothesis that many of the radial fracture systems on Venus are the surface manifestation of subsurface dikes propagating laterally from a central magma source. If this interpretation is correct, studies of terrestrial dikes suggest that the lineament directions, with localized exceptions and barring subsequent deformation, should be perpendicular to the orientation of the least compressive stress at the time of their formation. To test this hypothesis, we briefly examine a radial fracture system (63.7 degrees N, 195 degrees E) located between two deformation belts in Vinmara Planitia, and verify that the lineaments to the east behave in the expected manner. We have also chosen this feature, however, because of its proximity to Lauma Dorsa to the west. On the basis of Venera 15/16 data, both compressional and extensional origins for this deformation belt have been proposed. By examining the stratigraphy and applying our interpretation that the fracture system is linked to the presence of subsurface dikes, we present an independent evaluation of the stress state associated with Lauma Dorsa, and thus contribute to the assessment of its origin.

Grosfils, Eric B.↗

Microwave scattering and emission properties of large impact craters on the surface of Venus

Many of the impact craters on Venus imaged by the Magellan synthetic aperture radar (SAR) have interior floors with oblique incidence angle backscatter cross sections 2 to 16 times (3 dB to 12 dB) greater than the average scattering properties of the planet's surface. Such high backscatter cross sections are indicative of a high degree of wavelength-scale surface roughness and/or a high intrinsic reflectivity of the material forming the crater floors. Fifty-three of these (radar) bright floored craters are associated with 93 percent of the parabolic-shaped radar-dark features found in the Magellan SAR and emissivity data, features that are thought to be among the youngest on the surface of Venus. It was suggested by Campbell et al. that either the bright floors of the parabolic feature parent craters are indicative of a young impact and the floor properties are modified with time to a lower backscatter cross section or that they result from some property of the surface or subsurface material at the point of impact or from the properties of the impacting object. As a continuation of earlier work we have examined all craters with diameters greater than 30 km (except 6 that were outside the available data) so both the backscatter cross section and emissivity of the crater floors could be estimated from the Magellan data.

Stacy, N. J. S.↗

STILGAR End-of-Project Report

The Subsurface Tunnel Imaging LeveraGed by Analysis of Rayleigh wave ellipticity (STILGAR) project demonstrated an integrated geophysical approach for detecting, locating, and characterizing underground structural changes using dense seismic arrays and advanced inversion techniques. Field campaigns were conducted at two operational mines—the Redmond salt mine (Utah) and Graymont Pleasant Gap limestone mine (Pennsylvania)—providing real-world testbeds for monitoring anthropogenic subsurface activity. At the Redmond salt mine, seismic interferometry combined with back-projection inversion successfully identified continuous, low-amplitude signals from mining operations. The approach differentiated stationary from migrating anthropogenic sources, captured daily operational cycles, and validated the potential of passive seismic monitoring for remote detection of underground activity. At the Graymont Pleasant Gap mine, two dense seismic deployments in the spring and fall of 2023 generated over 4 TB of high-resolution data. Key outcomes included the relocation of 199 underground and 8 surface explosions with accuracies within tens of meters and the development of a 3D P-wave velocity model using the triple-difference tomography algorithm (tomoTD) that resolved major structural features such as the mine entrance, low-velocity tunnels, and roof-collapse areas. Ambient noise cross-correlation and back-projection analyses revealed persistent sources linked to ongoing mining activity, whereas horizontal-to-vertical spectral ratio (HVSR) and ellipticity studies confirmed stable site responses across seasons and identified soil thickness trends consistent with regional erosional and depositional processes. Checkerboard and sensitivity tests further validated the robustness of the tomographic results. Overall, the findings emphasize that although significant progress has been made in subsurface imaging, further work is needed to enhance the detection and localization of underground structures. Accurate imaging requires higher frequencies, yet anthropogenic sources tend to dominate the seismic record at those frequencies, and high-frequency surface waves are affected by higher modes that complicate interpretation. The improved detection and localization of human-induced signals enabled detailed temporal and spatial mapping of daily mine operations, demonstrating the feasibility of continuous anthropogenic source monitoring. Sensitivity to signals from nontraditional sources, such as fan operations, highlights the broader applicability of this approach to other industrial environments in which continuous and impulsive signals are present. The field campaigns produced a substantial volume of high-quality seismic data, supporting the development and testing of new methods for seismic source characterization and subsurface imaging. Future deployments should include sensors capable of recording lower frequencies to probe deeper structures, increase bandwidth to enhance resolution and sensitivity to both shallow and deep targets, and collect additional large-scale datasets to refine imaging and source characterization techniques. Moreover, conducting 3D modeling studies of seismic wavefields at higher frequencies will provide a better understanding of wave scattering and cavity–wavefield interactions in complex underground environments. In conclusion, the STILGAR project demonstrated that integrated seismic monitoring can effectively characterize underground operations, capturing both natural and anthropogenic signals. The approaches developed provide a foundation for improved detection, localization, and imaging of subsurface structures and are directly transferable to broader industrial monitoring applications.

58 GEOSCIENCES↗

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

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

54 ENVIRONMENTAL SCIENCES↗

Geophysical Impacts and Spectroscopic Identification of a Hydrous Iron Sulfate on Icy Worlds

Over geologic time-scales, large volumes of exogenic sulfur ions from Io's plasma torus have been supplied to the surface of Europa and Ganymede, which, combined with recent interpretations of orbiter images, dynamical modeling, and surface-subsurface exchange, suggests further sulfur transport into the interior of the icy worlds. These observations motivate mixed-phase spectral modeling for interpreting orbiter spectroscopy data and determination of hydration states of candidate surface materials including hydrous sulfates. In this work, we present a combined experimental and theoretical study of the low temperature and high pressure vibrational spectral signature of the iron-sulfate monohydrate endmember, szomolnokite (FeSO 4 ·H 2 O). By employing synchrotron Fourier-transform infrared spectroscopy (FTIR) in the diamond anvil cell up to 23 GPa and down to 20 K, we explore the extreme range of pressure-temperature domains relevant to icy environments throughout our solar system and beyond. Combined with our density-functional theory quantum-mechanics molecular dynamics results, we demonstrate that experimentally observed infrared features in the O-H stretching region commonly associated with nH 2 O (n > 1) hydration states can be attributed to a pure monohydrate without the need for pressure-induced exsolved ice, other coexisting hydrous iron sulfates, or strong overtone and combination modes. We further discuss the possibility of lateral variations in density and shear properties on icy worlds associated with temperature variations and the high-pressure phases of kieserite group monohydrated sulfates.

Geosciences↗

Experimental remote sensing of subsurface temperature in natural ocean water

The first successful depth-resolved remote sensing measurements of subsurface ocean water temperature were obtained by spectral analysis of the 3400 per cm O-H stretching Raman band of liquid water. Raman spectral data were obtained from a research vessel at various depths from the surface to 10 meters below the surface in a tidal estuary. The temperature inferred from the spectra was consistent with ground truth temperature to within the shot noise limited accuracy of plus or minus 2 C. The performance of a future fully developed airborne laser Raman water temperature measurement system is estimated on the basis of these first tests.

Leonard, D. A.↗

Locating subsurface gravel with thermal imagery

A method was discussed for using 6 band thermal imagery to locate subsurface gravel deposits in vegetated areas. Geologic history is reviewed to select potential areas of study. An overflight was made using a thermal scanner. The data were processed with a computerized system to delineate areas showing a quartz signature radiated by a gravel deposit. The method was developed during a search for gravel on National Forest land in Louisiana. Processed data from thermal imagery was compared with known gravel deposits and exploratory drill hole logs. A high correlation was noted for a wide range of deposits, from commercial pits to trace deposits only a foot thick. Overburden at these sites varied from zero to sixty feet, near the maximum annual penetration by the thermal wave. It was concluded that the method can be used to locate buried gravel deposits and that more time and effort are needed to verify the usefulness for developing gravel pits adjacent to proposed construction sites.

Scholen, Douglas E.↗

Passive Microwave Remote Sensing of Soil Moisture

Microwave remote sensing provides a unique capability for direct observation of soil moisture. Remote measurements from space afford the possibility of obtaining frequent, global sampling of soil moisture over a large fraction of the Earth's land surface. Microwave measurements have the benefit of being largely unaffected by cloud cover and variable surface solar illumination, but accurate soil moisture estimates are limited to regions that have either bare soil or low to moderate amounts of vegetation cover. A particular advantage of passive microwave sensors is that in the absence of significant vegetation cover soil moisture is the dominant effect on the received signal. The spatial resolutions of passive Microwave soil moisture sensors currently considered for space operation are in the range 10-20 km. The most useful frequency range for soil moisture sensing is 1-5 GHz. System design considerations include optimum choice of frequencies, polarizations, and scanning configurations, based on trade-offs between requirements for high vegetation penetration capability, freedom from electromagnetic interference, manageable antenna size and complexity, and the requirement that a sufficient number of information channels be available to correct for perturbing geophysical effects. This paper outlines the basic principles of the passive microwave technique for soil moisture sensing, and reviews briefly the status of current retrieval methods. Particularly promising are methods for optimally assimilating passive microwave data into hydrologic models. Further studies are needed to investigate the effects on microwave observations of within-footprint spatial heterogeneity of vegetation cover and subsurface soil characteristics, and to assess the limitations imposed by heterogeneity on the retrievability of large-scale soil moisture information from remote observations.

Njoku, Eni G.↗

Advancing Geothermal Research: Fiscal Year 2025 Accomplishments Report

This is a summary of geothermal work done at the National Renewable Energy Laboratory (NREL) in Fiscal Year 2025. This year brought increased attention to the geothermal industry and NREL's geothermal research portfolio. With more than 70 active projects, NREL research spanned the areas of resource exploration and characterization; conventional and next-generation geothermal technologies; subsurface thermal energy storage; heating and cooling; co-production of geothermal with critical minerals and oil and gas; modeling and analysis leveraging expertise in data science and machine learning; and more.

15 GEOTHERMAL ENERGY↗

Advancing Geothermal Research: Fiscal Year 2025 Accomplishments Report

This is a summary of geothermal work done at the National Laboratory of the Rockies in Fiscal Year 2025. This year brought increased attention to the geothermal industry and NLR's geothermal research portfolio. With more than 70 active projects, NLR research spanned the areas of resource exploration and characterization; conventional and next-generation geothermal technologies; subsurface thermal energy storage; heating and cooling; co-production of geothermal with critical minerals and oil and gas; modeling and analysis leveraging expertise in data science and machine learning; and more.

15 GEOTHERMAL ENERGY↗

Machine Learning-Driven Quantification of CO2 Plume Dynamics at Illinois Basin Decatur Project Sites Using Microseismic Data

This study utilizes machine learning to quantify CO2 plume extents by analyzing microseismic data from the Illinois Basin Decatur Project (IBDP). Leveraging a unique dataset of well logs, microseismic records, and CO2 injection metrics, this work aims to predict the temporal evolution of subsurface CO2 saturation plumes. The findings illustrate that machine learning can predict plume dynamics, revealing vertical clustering of microseismic events over distinct time periods within certain proximities to the injection well, consistent with an invasion percolation model. The buoyant CO2 plume partially trapped within sandstone intervals periodically breaches localized barriers or baffles, which act as leaky seals and impede vertical migration until buoyancy overcomes gravity and capillary forces, leading to breakthroughs along vertical zones of weakness. Between different unsupervised clustering techniques, K-Means and DBSCAN were applied and analyzed in detail, where K-means outperformed DBSCAN in this specific study by indicating the combination of the highest Silhouette Score and the lowest Davies–Bouldin Index. The predictive capability of machine learning models in quantifying CO2 saturation plume extension is significant for real-time monitoring and management of CO2 sequestration sites. The models exhibit high accuracy, validated against physical models and injection data from the IBDP, reinforcing the viability of CO2 geological sequestration as a climate change mitigation strategy and enhancing advanced tools for safe management of these operations.

Iyegbekedo, Ikponmwosa↗

The nature of the subsurface of Mercury from microwave observations at several wavelengths

Interferometric observations of Mercury at 3.71 cm are used to determine the dielectric constant of the subsurface material at a depth of 50 cm to be 2.4 plus or minus 0.3. The surface emissivity in the radio range is determined by the value of the dielectric constant which is thus necessary in order to interpret the thermal properties of the regolith. New integrated-disk data at 6- and 18-cm wavelengths are found to be consistent with a regolith of particulate geological material in which radiative transfer of heat is important. The results are very similar to the lunar case where regolith density, dielectric constant, and loss tangent are observed to increase smoothly with depth in the subsurface.

Cuzzi, J. N.↗

Lunar and Planetary Science XXXV: Mars: Radar, Gamma Ray Spectrometer, and Cratering Mineralogy

The session "Mars: Radar, Gamma Ray Spectrometer, and Cratering Mineralogy" contained the following reports:The Dielectric Properties of Martian Soil Simulant JSC Mars-1 in the Range from 20Hz to 10kHz; Eastern Sahara Geology from Orbital Radar: Potential Analog to Mars On the Dielectric Properties of the Martian-like Surface Sediments; Radar Observations of Recent Mars Landing Sites; Sounding of Subsurface Water Through Conductive Media in Mars Analog Environments Using Transient Electromagnetics and Low Frequency Ground-penetrating Radar; Burial Depth of the Reservoirs of Hydrogen at the Equatorial Latitudes on Mars; Elemental Composition Variations for Large Dusty and Rocky Regions on Mars Using Gamma-Ray Data from the Mars Odyssey Gamma-Ray Spectrometer; The Distribution of Non-Volatile Elements on Mars: Mars Odyssey GRS Results; Using Mars Odyssey GRS Data to Assess Models for the Bulk Composition of Mars; Mars: The Terra Arabia Low Epithermal Neutron Flux Anomaly; The Isidis Basin of Mars: New Results from MOLA, MOC, and THEMIS; Remote Sensing of the Haughton Impact Structure (HIS): A Terrestrial Proof of Concept for Using the Remote Sensing of Martian Craters as a Probe of Subsurface Composition;and Thermal Emission Spectra of Impact Glass and Shocked Deccan Basalt from Lonar Crater, India and Implications for Remote Sensing of Mars.

Source record↗

The Cretaceous-Tertiary Impact Crater and the Cosmic Projectile that Produced it

Evidence gathered to date from topographic data, geophysical data, well logs, and drill-core samples indicates that the buried Chicxulub basin, the source crater for the approximately 65 Ma Cretaceous-Tertiary (K/T) boundary deposits, is approximately 300 km in diameter. A prominent topographic ridge and a ring of gravity anomalies mark the position of the basin rim at approximately 150 km from the center. Wells in this region recovered thick sequences of impact-generated breccias at 200-300 m below present sea level. Inside the rim, which has been severely modified by erosion following impact, the subsurface basin continues to deepen until near the center it is approximately 1 km deep. The best planetary analog for this crater appears to be the 270 km-diameter Mead basin on Venus. Seismic reflection data indicate that the central zone of downward displacement and excavation (the transient crater is approximately 130 km in diameter, consistent with previous studies of gravity anomaly data). Our analysis of projectile characteristics utilizes this information, coupled with conventional scaling relationships, and geochemical constraints on the mass of extraterrestrial material deposited within the K/T boundary layer. Results indicate that the Chicxulub crater would most likely be formed by a long-period comet composed primarily of nonsilicate materials (ice, hydrocarbons, etc.) and subordinate amounts (less than or equal to 50 percent) primitive chondritic material. This collision would have released the energy equivalent to between 4 x 10(exp 8) and 4 x 10(exp 9) megatons of TNT. Studies of terrestrial impact rates suggest that such an event would have a mean production rate of approximately 1.25 x 10(exp -9) y(exp -1). This rate is considerably lower than that of the major mass extinctions over the last 250 million years (approximately 5 x 10(exp -7) y(exp -1). Consequently, while there is substantial circumstantial evidence establishing the cause-effect link between the Chicxulub basin forming event and the K/T biological extinctions, the results of our analysis do not support models of impact as a common or singular causative agent of mass extinctions on Earth.

Sharpton, Virgil L.↗

A detailed model for simulation of catchment scale subsurface hydrologic processes

A catchment scale numerical model is developed based on the three-dimensional transient Richards equation describing fluid flow in variably saturated porous media. The model is designed to take advantage of digital elevation data bases and of information extracted from these data bases by topographic analysis. The practical application of the model is demonstrated in simulations of a small subcatchment of the Konza Prairie reserve near Manhattan, Kansas. In a preliminary investigation of computational issues related to model resolution, we obtain satisfactory numerical results using large aspect ratios, suggesting that horizontal grid dimensions may not be unreasonably constrained by the typically much smaller vertical length scale of a catchment and by vertical discretization requirements. Additional tests are needed to examine the effects of numerical constraints and parameter heterogeneity in determining acceptable grid aspect ratios. In other simulations we attempt to match the observed streamflow response of the catchment, and we point out the small contribution of the streamflow component to the overall water balance of the catchment.

Paniconi, Claudio↗

Multi-resolution Arctic Shrub Cover Dataset Derived from UAS and Airborne SfM and LiDAR (2013-2025)

We synthesized 177 unoccupied aerial system flights and 77 airborne flights across the Arctic and created a multi-resolution benchmark data of low-to-tall shrub fractional cover leveraging Structure-from-Motion and Light Detection and Ranging. The resulting dataset covered a total of 1899 km2 across Alaska, Western Canada, Sweden, and Siberian Arctic, including key sites from the Oro Arctic to the High Arctic. The dataset is organized into 6 primary data collection directories (“Abisko,” “AWI,” “ERE,” “Fairbanks,” “NGEE,” “Toolik”), each containing site and flight subdirectories. Flight directories include shrub cover rasters (*.tifs) at 1 m, 5 m, and 30 m resolution, the canopy height model at 1 m resolution (*.tifs), and a bounding box *.kml file. For the AWI, Abisko, NGEE, and Fairbanks collections, we also include the GCC raster at 1 m resolution (*.tif). Files are organized by Collection > Site > Flight Name > Data Files. Flight rasters are in the local UTM zone and the .kml files are in the geographic coordinate system EPSG 4326. We also include a .csv file that details the source datasets for every flight. The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research. Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

canopy height model↗

Rover mounted ground penetrating radar as a tool for investigating the near-surface of Mars and beyong

In spite of the highly successful nature of recent planetary missions to the terrestrial planets and outer satellites a number of questions concerning the evolution of their surfaces remain unresolved. For example, knowledge of many characteristics of the stratigraphy and soils comprising the near-surface on Mars remains largely unknown, but is crucial in order to accurately define the history of surface processes and near-surface sedimentary record. Similar statements can be made regarding our understanding of near-surface stratigraphy and processes on other extraterrestrial planetary bodies. Ground penetrating radar (GPR) is a proven and standard instrument capable of imaging the subsurface at high resolution to 10's of meters depth in a variety of terrestrial environments. Moreover, GPR is portable and easily modified for rover deployment. Data collected with a rover mounted GPR could resolve a number of issues related to planetary surface evolution by defining shallow stratigraphic records and would provide context for interpreting results of other surface analyses (e.g. elemental or mineralogical). A discussion of existing GPR capabilities is followed first by examples of how GPR might be used to better define surface evolution on Mars and then by a brief description of possible GPR applications to the Moon and other planetary surfaces.

Grant, J. A.↗