Search NASA⌕ Search

SEARCH · Search NASA

Results for “geologic structure”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Molecular Origin of Wettability Alteration of Subsurface Porous Media upon Gas Pressure Variations

Upon extraction/injection of a large quantity of gas from/into a subsurface system in shale gas production or carbon sequestration, the gas pressure varies remarkably, which may significantly change the wettability of porous media involved. Mechanistic understanding of such changes is critical for designing and optimizing a related subsurface engineering process. Using molecular dynamics simulations, in this work we have calculated the contact angle of a water droplet on various solid surfaces (kerogen, pyrophyllite, calcite, gibbsite, and montmorillonite) as a function of CO 2 or CH 4 gas pressure up to 200 atm at a temperature of 300 K. The calculation reveals a complex behavior of surface wettability alteration by gas pressure variation depending on surface chemistry and structure, and molecular interactions of fluid molecules with surfaces. As the CO 2 gas pressure increases, a partially hydrophilic kerogen surface becomes highly hydrophobic, while a calcite surface becomes more hydrophilic. Considering kerogen and calcite being the major components of a shale formation, we postulate that the wettability alteration of a solid surface induced by a gas pressure change may play an important role in fluid flows in shale gas production and geological carbon sequestration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Rock Valley Accelerated Weight Drop Preliminary P-wave Tomographic Model

An active source experiment using an accelerated weight drop was conducted in Rock Valley, Nevada National Security Site, during the spring of 2021 in order to characterize the shallow seismic structure of the region. P-wave first arrival travel times picked from this experiment were used to construct a preliminary 3-D compressional wave speed model over an area that is roughly 4 km wide east-west and 8 km north-south to a depth of about 500-600 m below the surface, but with primary data concentration along the transects of the experimental lines. The preliminary model shows good correlation with basic geology and surface features, but geological interpretation is not the focus of this report. We describe the methods used in the tomographic inversion of the data and show results from this preliminary P-wave model.

42 ENGINEERING↗
DOE OSTI2020

GBCGE Subsurface Database Explorer and APIs

This submission defines a DOI for the Great Basin Center for Geothermal Energy's (GBCGE) Subsurface Database Explorer web application and underlying data services, and acknowledges the INGENIOUS project as a major source of funding for data compilation and quality assurance. The GBCGE Subsurface Database Explorer is an interactive web mapping application that provides public access to the GBCGE Subsurface Database, and its collection of datasets pertinent to geothermal exploration, oil and gas exploration, critical mineral exploration, and other subsurface characterization for the Great Basin Region, western US. This is a living database, and will be continuously updated with new data and datasets as funding and motivations allow. The underlying database views that populate the web application are on an automated refresh schedule. Data sources and acknowledgements: We thank our partners with the Nevada Division of Minerals (NDOM), the Southern Methodist University (SMU), and Great Basin State Geological Surveys for their active efforts in data curation, schema design, and quality assurance. We also thank contributors among the USGS, Oregon Institute of Technology, State Divisions of Water Resources, State Divisions of Oil, Gas, and Minerals, and State Geological Surveys for open data availability and direct contributions made under the National Geothermal Data System (NGDS).

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Quantifying particle-scale 3D granular dynamics during rapid compaction from time-resolved in situ 2D x-ray images

Understanding the particle-scale dynamics of granular materials during rapid compaction and flow is of fundamental importance for manufacturing, planetary science, geology, and defense applications. Time-resolved 2D radiography and static 3D x-ray tomography are powerful in situ tools for studying particle-scale dynamics but provide detail only in 2D or with significant time-scale limitations, respectively. Here, we introduce a new method that uses 2D in situ x-ray imaging for determining time-resolved 3D particle-scale dynamics in rapidly compressed granular materials. The method employs initial particle packing structures obtained from x-ray tomography, a 2D x-ray image generation algorithm, and an optimization algorithm. We first describe and validate the method using finite element simulations. Furthermore, we then apply the technique to x-ray phase-contrast images obtained during rapid compaction of granular materials with varying particle morphology and sample thickness. The depth-resolved particle-scale dynamics reveal complex velocity and porosity fields evolving heterogeneously along and perpendicular to the compaction direction. We characterize these features, their fluctuations near the compaction front, and the compaction front thickness. Our technique can be applied to understanding granular dynamics during rapid compaction events, and rearrangements during slower, but non-quasi-static, flows.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Novel ductile wellbore cementitious composite for geologic CO 2 storage

CO 2 leakage through damaged wellbore cement sheaths is a major risk of geologic CO 2 storage (GCS), as conventional wellbore cement is brittle and can be damaged due to acid attack and downhole stresses during CO 2 injection. Here we examine a novel fiber-reinforced engineered cementitious composite (ECC) proposed as a substitute to conventional wellbore cement due to its superior ductility and intrinsic crack width control. ECC and conventional wellbore cement coupons were exposed to water in equilibrium with CO 2 at 50 °C and 10 MPa. The samples were retrieved after several days and their mechanical performance was evaluated using a four-point bending test, microhardness, and compressive strength analyses. Optical microscopy and mercury intrusion porosimetry were used to characterize the progression of the carbonation front and pore structures of the specimens. Control experiments were conducted under the same temperature and pressure conditions but with a N 2 headspace to isolate the impact of CO 2 . It was found that carbonation increased the ultimate flexural strength of ECC but decreased its ductility. However, the ductility of carbonated ECC remained higher than that of conventional wellbore cements that exhibited brittle failure under all test conditions. Additionally, ECC exhibited minimal material loss and continued resistance to deformation in comparison to conventional wellbore cements. This suggests that while the exposure of ECC to CO 2 will alter its mechanical properties, altered ECC will continue to exhibit mechanical performance superior to conventional wellbore cement, and therefore shows promise as a highly durable wellbore cementing material for GCS applications.

58 GEOSCIENCES↗
DOE OSTI2021

Soil and Water Chemistry and Trace Metal Extractability and Speciation in Wetland Soils from Illinois and South Carolina and Stream Sediments from Tennessee

Dataset revised on October 15, 2021. This revision adds sulfur and iron X-ray absorption near-edge structure spectra for the wetland soils and stream sediments from the field areas. It also renames the sample locations in a way that is more intuitive to readers of the companion paper that is under review. Finally, the data filenames and organization have been updated in their labeling to parallel the data sources in the associated paper. The abstract text and methods were also revised to reflect the data that was added to the dataset.Trace metals are essential for microbially-mediated biogeochemical processes occurring in anoxic wetland soils and stream bed sediments, such as denitrification, methanogenesis, and mercury methylation. Low availability of these elements may potentially inhibit key components of anaerobic carbon and nitrogen cycling and contaminant transformation. The solid-phase speciation of trace metals likely plays an important role in controlling their bioavailability. Metal speciation is well studied in contaminated soils and sediments as well as those naturally elevated in trace metals. However, less is known regarding the chemical forms of trace metals in systems having concentrations similar to geological background levels, the very settings where metal limitations may be most prevalent. We have investigated trace metal concentrations, extractability, and solid-phase speciation in three freshwater subsurface aquatic systems: marsh wetland soils, riparian wetland soils, and the sediments of a streambed.Data are provided for marsh wetland soils at Argonne National Laboratory, riparian wetland soils in the Tims Branch watershed at Savannah River National Laboratory, and stream bed sediments from East Fork Poplar Creek near Oak Ridge National Laboratory. Soil and sediment elemental abundances, mineralogy, and extractable nutrients as well as dissolved major elements, anions, trace metals, and nutrients in the overlying surface waters are provided. In addition, the results of sequential chemical extraction for the trace metals cobalt, nickel, copper, and zinc from the soils and sediment are reported as well as X-ray absorption near-edge structure (XANES) spectra in these materials are reported. To aid interpretation of these data, XANES spectra of sulfur in the soils and sediments as well as both XANES and extended X-ray absorption fine structure (EXAFS) spectra of iron in these materials are reported. The data package also includes the XANES spectra of reference standards and a potential interferent in the measurements. All data are provided in text-based CSV format with header sections indicating the data contained in each file and the corresponding units. Note that "u" is used in place of Greek lower case mu to indicate the micro prefix on units.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Learning the factors controlling mineral dissolution in three-dimensional fracture networks: applications in geologic carbon sequestration

We perform a set of high-fidelity simulations of geochemical reactions within three-dimensional discrete fracture networks (DFN) and use various machine learning techniques to determine the primary factors controlling mineral dissolution. The DFN are partially filled with quartz that gradually dissolves until quasi-steady state conditions are reached. At this point, we measure the quartz remaining in each fracture within the domain as our primary quantity of interest. We observe that a primary sub-network of fractures exists, where the quartz has been fully dissolved out. This reduction in resistance to flow leads to increased flow channelization and reduced solute travel times. However, depending on the DFN topology and the rate of dissolution, we observe substantial variability in the volume of quartz remaining within fractures outside of the primary subnetwork. This variability indicates an interplay between the fracture network structure and geochemical reactions. We characterize the features controlling these processes by developing a machine learning framework to extract their relevant impact. Specifically, we use a combination of high-fidelity simulations with a graph-based approach to study geochemical reactive transport in a complex fracture network to determine the key features that control dissolution. We consider topological, geometric and hydrological features of the fracture network to predict the remaining quartz in quasi-steady state. We found that the dissolution reaction rate constant of quartz and the distance to the primary sub-network in the fracture network are the two most important features controlling the amount of quartz remaining. This study is a first step towards characterizing the parameters that control carbon mineralization using an approach with integrates computational physics and machine learning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Groundwater Circulation Within the Mountain Block: Combining Flow and Transport Models With Magnetotelluric Observations to Untangle Its Nested Nature

Abstract Mountains are vital water sources for humans and ecosystems, continuously replenishing lowland aquifers through surface runoff and mountain recharge. Quantifying these fluxes and their relative importance is essential for sustainable water resource management. However, our mechanistic understanding of the flow and transport processes determining the connection between the mountain block and the basin aquifer remains limited. Traditional conceptualizations assume groundwater circulation within the mountain block is predominantly shallow. This view neglects the role of deep groundwater flowpaths significantly contributing to the water, solute, and energy budgets. Overcoming these limitations requires a holistic characterization of the multiscale nature of groundwater flow along the mountain‐to‐valley continuum. As a proof‐of‐concept, we use a coupled groundwater flow and transport model to design a series of numerical experiments that explore the role of geology, topography, and weathering rates in groundwater circulation and their resulting resistivity patterns. Our results show that accumulating solutes near stagnation zones create contrasting electrical resistivity patterns that separate local, intermediate, and regional flow cells, presenting a target for magnetotelluric observations. To demonstrate the sensitivity of magnetotelluric data to features in our resistivity models, we use the MARE2DEM electromagnetic modeling code to perform forward and inverse simulations. This study highlights the potential of magnetotelluric surveys to image the resistivity structure resulting from multiscale groundwater circulation through relatively impervious crystalline basement rocks in mountainous terrains. This capability could change our understanding of the critical zone, offering a holistic perspective that includes deep groundwater circulation and its role in conveying solutes and energy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Locating Seismic Events with Local-Distance Data

As the seismic monitoring community advances toward detecting, identifying, and locating ever-smaller natural and anthropogenic events, the need is constantly increasing for higher resolution, higher fidelity data, models, and methods for accurately characterizing events. Local-distance seismic data provide robust constraints on event locations, but also introduce complexity due to the significant geologic heterogeneity of the Earth’s crust and upper mantle, and the relative sparsity of data that often occurs with small events recorded on regional seismic networks. Identifying the critical characteristics for improving local-scale event locations and the factors that impact location accuracy and reliability is an ongoing challenge for the seismic community. Using Utah as a test case, we examine three data sets of varying duration, finesse, and magnitude to investigate the effects of local earth structure and modeling parameters on local-distance event location precision and accuracy. We observe that the most critical elements controlling relocation precision are azimuthal coverage and local-scale velocity structure, with tradeoffs based on event depth, type, location, and range.

42 ENGINEERING↗
DOE OSTI2021

Core Characterization of Bakken Shale from the Bedwell 33-52-1-1H Well, Sheridan County, Montana

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia were used to characterize core of the Bedwell 33-52-1-1H well (API 25-091-21920), a wildcat field in Sheridan County, Montana. The primary impetus of this work is to provide a unique dataset to researchers at universities, national laboratories, geological surveys, and other locations for additional analyses. The resultant datasets are presented in this report and can be accessed from NETL's Energy Data eXchange (EDX) online system using the following link: https://edx.netl.doe.gov/dataset/bedwell-33-52-1-1h-well. The equipment and techniques used to characterize the full core were non-destructive, enabling future examinations and analyses to be performed on these cores. However, none of the equipment used was suitable for direct visualization of the pore space in fine-grained structures such as those found in shales; fractures, discontinuities, and millimeter scale features were readily detectable with the methods tested. Imaging with the NETL medical CT scanner was performed on the entire core. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), P-wave, and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for more detailed analysis. A selection of samples (Table 1) were cored, sliced, and powdered for additional characterization analyses. These analyses did involve destructive subsampling of the full core at discrete sections, and the results are presented here in the context of the larger core description. Selected sections of the core were removed from the full core for detailed analyses using the following destructive techniques: bulk inductively coupled plasma-optical emission spectroscopy (ICP-OES), isolation of carbon and sulfur, Fourier transform infrared (FTIR) spectroscopy, scanning electron microscopy (SEM), and higher resolution CT scan with NETL’s micro-CT scanning systems. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provides a multi-scale analysis of the core; the resulting macro and micro descriptions are relevant to many subsurface energy related examinations traditionally performed at NETL.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2025

NE-COST plug-in: Expanding ACCERT's Capabilities for Life-Cycle Cost Modeling

The Algorithm for the Capital Cost Estimation of Reactor Technologies (ACCERT) is a structured methodology and software tool designed to simplify and standardize cost estimation for nuclear reactor technologies [1]. By utilizing a relational database structure and modular cost estimation algorithms, ACCERT delivers a robust, flexible, and scalable framework for evaluating costs across various reactor types and configurations [2]. The recent integration of the NE-COST plugin further expands ACCERT’s scope by introducing detailed life-cycle cost modeling and probabilistic analysis of uncertainties. This addition enables users to evaluate costs across front-end processes such as uranium enrichment and fabrication, as well as back-end activities including waste disposal and geologic storage. Through Monte Carlo statistical cost simulations, the plugin provides probabilistic insights into cost ranges, offering critical decision-making support for stakeholders including reactor developers, policymakers, and researchers.

Zhou, Jia↗
DOE OSTI2020

Surface and bulk characterization of reservoir and cap-rocks: Past, present, and future

Geologic formations provide potentially some of the largest volume capacities for CO 2 storage or sequestration. Potential storage sites can be deep saline aquifers, depleted oil reservoirs, and coal seams, surrounded by sealing layers to prevent CO 2 from leaking. It is therefore critical to understand mechanisms contributing to CO 2 trapping and CO 2 leaks. Both phenomena are governed by reactions at the interfaces of the reservoir and cap-rocks and are controlled by the complex chemistry and pore structures of rocks. Mechanisms at the macroscale are affected by the processes occurring at the nanoscale. This review highlights the necessity of multitechnique, multiscale characterization of rocks and points to the importance of surface analysis and surface science studies. Two shale rocks (seals) from Niobrara and Agardhfjellet formations with complex surface chemistry are used as examples throughout the paper. Typically, evaluation of rocks with x-ray diffraction, thermogravimetric analysis, Rock-Eval pyrolysis, gas adsorption, and electron microscopy combined with energy dispersive x-ray spectroscopy is conducted to provide valuable information about the bulk mineralogy, elemental composition, pore volume, and adsorbed species on the sample surface. These studies are necessary prior to designing surface sensitive experiments with x-ray photoelectron spectroscopy (XPS), guiding both sample preparation and sample analysis. XPS has been widely used to study the surface composition of rocks during the investigations of their fine-scale wettability, and the main findings are highlighted here. This paper also reviews the existing literature on ambient-pressure XPS, which provides new opportunities to study in situ chemical alteration due to interactions with CO 2 and offers recommendations for adapting this technique to study rock-fluid interactions, allowing for the identification of fundamental interactions during CO 2 sequestration and guide selection of formation sites for improved storage.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Mineral detection of neutrinos and dark matter. A whitepaper

Minerals are solid state nuclear track detectors — nuclear recoils in a mineral leave latent damage to the crystal structure. Depending on the mineral and its temperature, the damage features are retained in the material from minutes (in low-melting point materials such as salts at a few hundred °C) to timescales much larger than the 4.5 Gyr-age of the Solar System (in refractory materials at room temperature). The damage features from the O (50) MeV fission fragments left by spontaneous fission of 238U and other heavy unstable isotopes have long been used for fission track dating of geological samples. Laboratory studies have demonstrated the readout of defects caused by nuclear recoils with energies as small as O (1) keV. This whitepaper discusses a wide range of possible applications of minerals as detectors for E R ≳ O (1) keV nuclear recoils: Using natural minerals, one could use the damage features accumulated over O (10) Myr–O (1) Gyr to measure astrophysical neutrino fluxes (from the Sun, supernovae, or cosmic rays interacting with the atmosphere) as well as search for Dark Matter. Using signals accumulated over months to few-years timescales in laboratory-manufactured minerals, one could measure reactor neutrinos or use them as Dark Matter detectors, potentially with directional sensitivity. Research groups in Europe, Asia, and America have started developing microscopy techniques to read out the O (1)–O (100) damage features in crystals left by O (0.1)–O (100) keV nuclear recoils. We report on the status and plans of these programs. Furthermore, the research program towards the realization of such detectors is highly interdisciplinary, combining geoscience, material science, applied and fundamental physics with techniques from quantum information and Artificial Intelligence.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Acoustic Impedance Inversions for Offshore CO 2 Storage: South Georgia Embayment, United States

This study builds and expands on previous CO 2 storage resource assessment studies of the southeastern offshore Atlantic margin by providing a detailed evaluation on how rock porosities and permeabilities are distributed across the Upper Cretaceous strata restricted to the South Georgia Embayment (SGE). Using legacy industry two-dimensional seismic reflection and well data, this assessment is the first application of multiple seismic inversion techniques in this area. This workflow provides a reliable and repeatable model-based inversion which gives an improved image to discriminate lithology and predict porosity. The workflow is applicable to future CO 2 storage resource assessment studies elsewhere. The inversion results indicate that distinct porosity and permeability regimes are present and distributed in the Upper Cretaceous strata within the SGE. The impedance and porosity relationships show well-founded and reliable correlation. These relationships reveal low impedance coincident to the high porosity intervals which are proposed as potential reservoir intervals for CO 2 storage. In addition, the result shows that the Upper Cretaceous strata have two main potential reservoirs in the lower part. These are overlain by a thick impermeable interval, mostly shale, which has high impedance, low porosity, and low permeability and extends within the SGE. This result is in agreement with a previous study that also proposed two significant storage reservoirs for CO 2 in the Upper Cretaceous strata. Since porosity distribution is estimated using multiple methods, it follows the trends of seismic signature and structures of the Upper Cretaceous strata. The extracted values of porosity, ranging from 15 to 36%, and permeability, ranging from 1 to 100 mD, are close to the measured values from the well core data at the Upper Cretaceous strata interval.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Multi-Scale Seismic Measurements for Site Characterization and CO2 Monitoring in an Enhanced Oil Recovery/Carbon Capture, Utilization, and Sequestration Project, Farnsworth Field, Texas

To address the challenges of climate change, significantly more geologic carbon sequestration projects are beginning. The characterization of the subsurface and the migration of the plume of supercritical carbon dioxide are two elements of carbon sequestration that can be addressed through the use of the available seismic methods in the oil and gas industry. In an enhanced oil recovery site in Farnsworth, TX, we employed three separate seismic techniques. The three-dimensional (3D) surface seismic survey required significant planning, design, and processing, but produces both a better understanding of the subsurface structure and a three-dimensional velocity model, which is essential for the second technique, a timelapse vertical seismic profile, and the third technique, cross-well seismic tomography. The timelapse 3D Vertical Seismic Profile (3D VSP) revealed both significant changes in the reservoir between the second and third surveys and geo-bodies that may represent the extent of the underground carbon dioxide. The asymmetry of the primary geo-body may indicate the preferential migration of the carbon dioxide. The third technique, cross-well seismic tomography, suggested a strong correlation between the well logs and the tomographic velocities, but did not observe changes in the injection interval.

Energy & Fuels↗
DOE OSTI2021

Convolutional Neural Networks Based Remote Sensing Scene Classification under Clear and Cloudy Environments

Remote sensing (RS) scene classification has wide applications in the environmental monitoring and geological survey. In the real-world applications, the RS scene images taken by the satellite might have two scenarios: clear and cloudy environments. However, most of existing methods did not consider these two environments simultaneously. In this paper, we assume that the global and local features are iscriminative in either clear or cloudy environments. Many existing Convolution Neural Networks (CNN) based models have made excellent achievements in the image classification, however they somewhat ignored the global and local features in their network structure. In this paper, we propose a new CNN based network (named GLNet) with the Global Encoder and Local Encoder to extract the discriminative global and local features for the RS scene classification, where the constraints for inter-class dispersion and intra-class compactness are embedded in the GLNet training. The experimental results on two publicized RS scene classification datasets show that the proposed GLNet could achieve better performance based on many existing CNN backbones under both clear and cloudy environments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Validation of Fault Displacements from Dynamic Rupture Simulations against the Observations from the 1992 Landers Earthquake

ABSTRACT Coseismic fault displacements in large earthquakes have caused significant damage to structures and lifelines on and near fault lines. Coseismic displacements represent a real threat, especially to distributed infrastructure systems. For infrastructure systems that can not avoid active faults, engineering displacement demands are defined using probabilistic fault-displacement hazard analyses (PFDHA). However, PFDHA models are sparse and poorly constrained partly due to the scarcity of detailed fault-displacement observations. Advancements in dynamic rupture simulation methods make them an attractive approach to address this important issue. Because fault displacements can be simulated for various geologic conditions as constrained by current knowledge about earthquake processes, they can be used to supplement the observation datasets. In addition to providing on-fault displacements, when used with appropriate constitutive models for the bulk medium, they can capture off-fault distributed inelastic deformations as well. For viable extrapolation, simulations must first be validated against data. In this article, we summarize the calibration and validation of the dynamic rupture model against the observations of the well-documented 1992 Landers earthquake. We defined a preferred model that reproduces several first-order fault-displacement metrics such as the on-fault partition of the total displacement, the mean fault-zone width, and the location of the peak displacement. Simulated ground motions consistent with the observations ensure that all physics important to modeling have been properly parameterized. For the extrapolation, we generated a suite of dynamic rupture models to quantify expected fault-displacement metrics, their intercorrelations, and magnitude dependencies, which are in part supported by the Landers and other recent earthquakes. Our validation and extrapolation exercise paves the way for using dynamic rupture modeling to quantitatively address fault-displacement hazard on a broader scale. The results are promising and are expected to be useful to inform PFDHA model development.

Geochemistry & Geophysics↗
DOE OSTI2025

Geospatial modeling of near subsurface temperatures of the contiguous United States for assessment of materials degradation

Abstract Understanding subsurface temperature variations is crucial for assessing material degradation in underground structures. This study maps subsurface temperatures across the contiguous United States for depths from 50 to 3500 m, comparing linear interpolation, gradient boosting (LightGBM), neural networks, and a novel hybrid approach combining linear interpolation with LightGBM. Results reveal heterogeneous temperature patterns both horizontally and vertically. The hybrid model performed best achieving a root mean square error of 2.61 °C at shallow depths (50–350 m). Model performance generally decreased with depth, highlighting challenges in deep temperature prediction. State-level analyses emphasized the importance of considering local geological factors. This study provides valuable insights for designing efficient underground facilities and infrastructure, underscoring the need for depth-specific and region-specific modeling approaches in subsurface temperature assessment.

Science & Technology - Other Topics↗