Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic analysis”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Time-Lapse Vertical Seismic Profiling (VSP) for CO 2 Storage in a Depleted Oil Field in Northern Michigan

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. This report describes the Dover 33 VSP study to test the effectiveness of time-lapse Vertical Seismic Profile (VSP) for detecting and delineating a plume of more than 271,000 tonnes of CO2 injected into the Brown Niagaran and A-1 Carbonate formations within the Dover 33 reef between March 2013 and September 2016. Five 2D walkaway VSP (WVSP) source lines were acquired by SIGMA3 in September 2016 to investigate the possible time-lapse response in both P-wave and PS-wave seismic data. The data was compared to the same survey geometry acquired in March 2013 by SR2020. The seismic energy was recorded into an 80-level, three-component geophone array deployed into the Dover 1-33 well and placed just above the reef. P-wave and PS-wave reflection images were produced for each of the three source lines from both the 2013 baseline and 2016 monitor surveys. The images were then compared to look for changes in the reflectivity at and around the injection location that might indicate how the CO 2 has moved over this time period. This technique revealed several localized areas with sizable impedance differences inside the reef where CO 2 would be expected; however, a large number of similar impedance “hotspots” were also detected outside the reef in areas where injected CO 2 would not be expected. Therefore, a second analysis that involved calculating P-wave and S-wave travel time differences between the 2013 and 2016 VSPs was conducted to look for a change that could be caused by the CO 2 plume.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS1993

An analysis of nearfield normal mode amplitude anomalies of the Landers earthquake

The 1992 Landers earthquake (M(sub w) = 7.3) occurred in the middle of the TERRAscope network. Long-period Rayleigh waves recorded at the TERRAscope stations (delta less than or = 3 deg) after traveling around the Earth show large amplitude anomalies, one order of magnitude larger than spherical Earth predictions up to a period of about 600 s. The ground motions over the epicentral region at and after the arrival of R4-5 are in phase at all stations. These observations are inconsistent with the nearly vertical strike slip mechanism of the Landers earthquake. Synthetic seismograms for a rotating, elliptic, and laterally heterogeneous Earth model calculated by the variational method agree well with the observed waveforms. Calculations for various 3D Earth models demonstrate that the amplitudes are very sensitive to the large scale aspherical structure in the crust and the mantle. The anomalies for modes shorter than 300 s period can be explained by lateral heterogeneity shallower than the upper mantle. Rotation of the Earth and lower mantle heterogeneity are required to explain mode amplitudes at longer periods. Current whole mantle seismic tomographic models can fully explain the observed amplitudes longer than 300 s. To assess the effect of the high order lateral heterogeneity in the mantle more precise estimate of the crustal correction is required.

Watada, Shingo↗
DOE OSTI2020

Internal structure of the San Jacinto fault zone at the Ramona Reservation, north of Anza, California, from dense array seismic data

SUMMARY We image the internal structure of the San Jacinto fault zone (SJFZ) near Anza, California, with seismic data recorded by two dense arrays (RA and RR) from ∼42 000 local and ∼180 teleseismic events occurring between 2012 and 2017. The RA linear array has short aperture (∼470 m long with 12 strong motion sensors) and recorded for the entire analysed time window, whereas the RR is a large three-component nodal array (97 geophones across a ∼2.4 km × 1.4 km area) that operated for about a month in September–October 2016. The SJFZ at the site contains three near-parallel surface traces F1, F2 and F3 from SW to NE that have accommodated several Mw > 6 earthquakes in the past 15 000 yr. Waveform changes in the fault normal direction indicate structural discontinuities that are consistent with the three fault surface traces. Relative slowness from local events and delay time analysis of teleseismic arrivals in the fault normal direction suggest a slower SW side than the NE with a core damage zone between F1 and F2. This core damage zone causes ∼0.05 s delay at stations RR26–31 in the teleseismic P arrivals compared with the SW-most station, and generates both P- and S-type fault zone trapped waves. Inversion of S trapped waves indicates the core damaged structure is ∼100 m wide, ∼4 km deep with a Q value of ∼20 and 40 per cent S-wave velocity reduction compared with bounding rocks. Fault zone head waves observed at stations SW of F3 indicate a local bimaterial interface that separates the locally faster NE block from the broad damage zone in the SW at shallow depth and merges with a deep interface that separates the regionally faster NE block from rocks to the SW with slower velocities at greater depth. The multiscale structural components observed at the site are related to the geological and earthquake rupture history at the site, and provide important information on the preferred NW propagation of earthquake ruptures on the San Jacinto fault.

Geochemistry & Geophysics↗
DOE OSTI2022

GraphAlign: Graph-Enabled Machine Learning for Seismic Event Filtering

This report summarizes results from a 2 year effort to improve the current automated seismic event processing system by leveraging machine learning models that can operated over the inherent graph data structure of a seismic sensor network. Specifically, the GraphAlign project seeks to utilize prior information on which stations are more likely to detect signals originating from particular geographic regions to inform event filtering. To date, the GraphAlign team has developed a Graphical Neural Network (GNN) model to filter out false events generated by the Global Associator (GA) algorithm. The algorithm operates directly on waveform data that has been associated to an event by building a variable sized graph of station waveforms nodes with edge relations to an event location node. This builds off of previous work where random forest models were used to do the same task using hand crafted features. The GNN model performance was analyzed using an 8 week IMS/IDC dataset, and it was demonstrated that the GNN outperforms the random forest baseline. We provide additional error analysis of which events the GNN model performs well and poorly against concluded by future directions for improvements.

58 GEOSCIENCES↗
DOE OSTI2022

Observations of Mantle Seismic Anisotropy Using Array Techniques: Shear-Wave Splitting of Beamformed SmKS Phases

Shear-wave splitting measurements are commonly used to resolve seismic anisotropy in both the upper and lowermost mantle. Typically, such techniques are applied to SmKS phases that have reflected (m-1) times off the underside of the core-mantle boundary before being recorded. Practical constraints for shear-wave splitting studies include the limited number of suitable phases as well as the large fraction of available data discarded because of poor signal-to-noise ratios (SNRs) or large measurement uncertainties. Array techniques such as beamforming are commonly used in observational seismology to enhance SNRs, but have not been applied before to improve SmKS signal strength and coherency for shear wave splitting studies. Here, we investigate how a beamforming methodology, based on slowness and backazimuth vespagrams to determine the most coherent incoming wave direction, can improve shear-wave splitting measurement confidence intervals. Through the analysis of real and synthetic seismograms, we show that (a) the splitting measurements obtained from the beamformed seismograms (beams) reflect an average of the single-station splitting parameters that contribute to the beam; (b) the beams have (on average) more than twice as large SNRs than the single-station seismograms that contribute to the beam; (c) the increased SNRs allow the reliable measurement of shear wave splitting parameters from beams down to average single-station SNRs of 1.3. Beamforming may thus be helpful to more reliably measure splitting due to upper mantle anisotropy. Moreover, here we show that beamforming holds potential to greatly improve detection of lowermost mantle anisotropy by demonstrating differential SKS–SKKS splitting analysis using beamformed USArray data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Stress Field Dynamics and Fault Slip Potential in the Paradox Basin

Abstract The Paradox Basin, straddling Utah, Colorado, Arizona, and New Mexico is characterized by an intricate amalgamation of evaporites and clastic layers and is dominated by prominent salt walls and related subsurface structures. Our research offers a new examination of the stress distribution across the basin, deriving from continuous and discrete stress measurements conducted in boreholes in the region and focal mechanism analysis, emphasizing variations over salt structures. Integrating Coulomb failure criteria with probabilistic methods, we assess potential fault movements resulting from fluid pressure alterations. Our approach provides a comprehensive understanding of the Paradox Basin's state of stress, showing a continuous change of the maximum horizontal stress orientation from N‐S at the Wasatch Fault Zone to WNW‐ESE in the northern part of the Paradox Basin and to WSW‐ENE in the southern part of the basin. Further East, into the Colorado Plateau and the Uncompahgre Uplift, the S H max orientation becomes E‐W. Decoding stress orientation dynamics has enabled critical insights into fault slip potential, especially in the basin's northern region. The salt wall faults are less likely to slip, and the Paradox Formation's evaporite and clastic rock sequence can serve as a potential low seismic risk target for carbon storage and hydrocarbon extraction.

Geochemistry & Geophysics↗
DOE OSTI2024

Characteristics of Vertical Ground Motions and Their Effect on the Seismic Response of Bridges in the Near-Field: A State-of-the-Art Review

Despite the evidence from past earthquakes and several numerical investigations demonstrating the detrimental impact of vertical ground motions (VGMs) on the integrity of bridge structures, incorporating their effects into seismic assessment and design procedures has traditionally been given limited consideration. Current codes utilize rather simplistic approaches to account for the concurrent effects of vertical and horizontal motions in structural performance evaluations, potentially leading to unconservative estimates of structural demands. This paper reviews the main features of VGMs and their effect on the seismic response of bridges. The methods and empirical models available to estimate vertical motions for design purposes are discussed, and research gaps and related research needs are identified. Finally, the emerging role of physics-based ground-motion simulations, as well as their limitations, in supporting future research and informing the development of simplified design procedures is examined. The main areas of interest for future research are identified in the need to carry out systematic sensitivity studies to gain insight into the main earthquake parameters that influence key VGMs features, understand the influence of soil nonlinearities on VGMs amplitude and frequency content, inform the development of empirical models that cover a range of site conditions and source-to-site distances where current models are poorly constrained, generate arrays of motions to update coherency models to properly inform the analysis of distributed infrastructure, investigate the impulsive character of VGMs, and assess the approximations made in estimating VGMs with 1D site response analyses. Furthermore, specific focus is laid on large-magnitude earthquakes in the near-field.

Asynchronous ground motion↗
DOE OSTI2021

Oregon Rivers Characterization Tool

DHS’s Regional Resiliency Assessment Program (RRAP) undertook the Oregon State Transportation Systems RRAP project between 2018-2021 to assess the seismic vulnerabilities of the state’s transportation system to a Cascadia Subduction Zone (CSZ) earthquake (CISA 2021). The project primarily focused on the effects of a CSZ earthquake on the mobility of emergency response efforts into the region. In the event that bridges located at crucial river crossings fail during a CSZ earthquake, rivers will become natural barriers affecting the ability to support response and recovery efforts. The characteristics of rivers (e.g., width, depth, and flow) are critical data to inform planning for the recovery of transportation systems across river barriers and determine reopening times following a CSZ earthquake. As part of the transportation system assessment, an Argonne National Laboratory (Argonne) project team developed a River Characterization Tool (RCT) to quantify the physical characteristics of the river system in Oregon, which can be used to inform a system-level analysis of transportation recovery and reopening times. This report mainly discusses approaches, methods, development, and applications of the RCT. The use of the river data generated from this tool for a system-level analysis of transportation recovery is discussed in the Oregon Transportation Systems RRAP Resiliency Assessment report (CISA 2021). The primary goal of the RCT is to generate statewide dataset of widths, depths, and flows for all rivers across Oregon to assess the effects of river barriers on emergency response and recovery efforts following a CSZ earthquake. The RCT leverages remote sensing data that multispectral sensors collect from satellites, employs a deep learning model to characterize river predictions, and computes river channel geometric dimensions with hydraulic relationships and river predictions. The entire process of river characterization is automated and implemented through the RCT.

58 GEOSCIENCES↗
DOE OSTI2023

Analysis of the 2021 Milford, Utah earthquake swarm: Enhanced earthquake catalog and migration patterns

An earthquake swarm occurred in Spring 2021 in south-central Utah near the town of Milford. The University of Utah Seismograph Stations located 125 earthquakes between March 19 and May 10 with magnitudes ranging from 0.5 to 3.2. We implement a matched-filter technique in order to identify additional earthquakes that went undetected during the routine network location. The 125 network-located earthquakes are used as templates and are cross-correlated with continuous data for the dates Feb 17—June 10. This time period corresponds to approximately 1 month before the earthquake swarm began through 1 month after it ended. For the matched-filter analysis, we rely heavily on station FOR1, which is located within 5 km of most template events. Four other stations within 20–30 km of most template events provide a supplement to the closest station. The matched-filter implementation results in the detection of over 600 earthquakes in addition to the original 125 catalog events. This is one of the largest swarms ever recorded in Utah, and no previous large swarms have been recorded in this location. We use HypoDD to obtain relative double-difference locations of the catalog events. Both routine locations and HypoDD relocations of the catalog events suggest a fault dipping west beneath the mountain range, opposite of typical Basin and Range normal faults that dip beneath alluvial valleys. Moment tensors for the largest five events show normal faulting consistent with the west-dipping fault seen in the seismicity. Hydrothermal features in the area, including a geothermal power plant, suggest that fluids may be a contributing factor to the earthquake swarm triggering. We examine the role of fluids by exploring bounds on diffusion parameters and investigating spatial migration characteristics of the swarm seismicity. We conclude that this swarm is the result of heterogenous stress conditions in a prefractured region.

58 GEOSCIENCES↗
NASA NTRS1974

Mineral resources, geological structure, and landform surveys

Diagnostic ERTS imagery has been used to pinpoint surface conditions associated with known mining districts. These include enhancements which depict hitherto unrecognized surface alteration and allow analysis of ore-controlling fractures distribution in a regional context. ERTS has likewise provided observational data containing previously unrecognized surface anomalies in large oil-producing basins which correlate closely with known oil fields. These observational data offer promise of providing new and powerful techniques for oil exploration, especially if further work using more sophisticated enhancement-processing proves capable of emphasizing the anomalies. ERTS is showing a better-than-anticipated potential for producing accurate small-scale (large-area) geologic maps, often containing details that were previously not recorded on similar regional maps. The maps produced from ERTS imagery can be prepared more effectively than previously possible, mainly because of the synoptic, multispectral, and repetitive character of ERTS data. ERTS has also provided extensive information on possible geologic hazards. Many new fractures have been identified in several regions of the Pacific Coast seismic belt that have histories of recent earthquakes. This has obvious implications for engineering projects such as dams, aqueducts, and transportation routes. In the mid-continent area, ERTS data have been used to predict zones of rooffall danger in a working coal mine from newly discovered lineations (probably fractures) used as indicators of hazards.

Short, N. M.↗
NASA NTRS2014

Non-Cooled Power System for Venus Lander

The Planetary Science Decadal Survey of 2013-2022 stated that the exploration of Venus is of significant interest. Studying the seismic activity of the planet is of particular importance because the findings can be compared to the seismic activity of Earth. Further, the geological and atmospheric properties of Venus will shed light into the past and future of Earth. This paper presents a radioisotope power system (RPS) design for a small low-power Venus lander. The feasibility of the new power system is then compared to that of primary batteries. A requirement for the power source system is to avoid moving parts in order to not interfere with the primary objective of the mission - to collect data about the seismic activity of Venus using a seismometer. The target mission duration of the lander is 117 days, a significant leap from Venera 13, the longest-lived lander on the surface of Venus, which survived for 2 hours. One major assumption for this mission design is that the power source system will not provide cooling to the other components of the lander. This assumption is based on high-temperature electronics technology that will enable the electronics and components of the lander to operate at Venus surface temperature. For the proposed RPS, a customized General Purpose Heat Source Radioisotope Thermoelectric Generator (GPHSRTG) is designed and analyzed. The GPHS-RTG is chosen primarily because it has no moving parts and it is capable of operating for long duration missions on the order of years. This power system is modeled as a spherical structure for a fundamental thermal analysis. The total mass and electrical output of the system are calculated to be 24 kilograms and 26 Watts, respectively. An alternative design for a battery-based power system uses Sodium Sulfur batteries. To deliver a similar electrical output for 117 days, the battery mass is calculated to be 234 kilograms. Reducing mission duration or power required will reduce the required battery mass. Finally, the advantages and disadvantages of both power systems with regard to science return, risk, and cost are briefly compared. The design of the radioisotope power system is considerably riskier because it is novel and would require additional years of further refinement, manufacturing, safety analysis, and testing that the primary batteries do not need. However, the lifetime of the radioisotope power system makes its science return more promising.

planetary↗
NASA NTRS1994

Lower mantle thermal structure deduced from seismic tomography, mineral physics and numerical modelling

The long-wavelength thermal anomalies in the lower mantle have been mapped out using several seismic tomographic models in conjunction with thermodynamic parameters derived from high-pressure mineral physics experiments. These parameters are the depth variations of thermal expansivity and of the proportionality factor between changes in density and seismic velocity. The giant plume-like structures in the lower mantle under the Pacific Ocean and Africa have outer fringes with thermal anomalies around 300-400 K, but very high temperatures are found in the center of the plumes near the base of the core-mantle boundary. These extreme values can exceed +1500 K and may reflect large hot thermal anomalies in the lower mantle, which are supported by recent measurements of high melting temperatures of perovskite and iron. Extremely cold anomalies, around -1500 K, are found for anomalies in the deep mantle around the Pacific rim and under South America. Numerical simulations show that large negative thermal anomalies in the mid-lower mantle have modest magnitudes of around -500 K. correlation pattern exists between the present-day locations of cold masses in the lower mantle and the sites of past subduction since the Cretaceous. Results from correlation analysis show that the slab mass-flux in the lower mantle did not conform to a steady-state nature but exhibited time-dependent behavior.

Cadek, O.↗
DOE OSTI2020

Prediction of gas hydrate saturation using machine learning and optimal set of well-logs

We report resistivity and acoustic logs are widely used to estimate gas hydrate saturation in various sedimentary systems using one of the two popular methods ((1) acoustic velocity and (2) electrical resistivity), but the limitations of these two methods are often overlooked, which include (i) well-specific calibration of empirical exponents in the electrical resistivity method, (ii) assumption of known pore morphology for gas hydrates in the acoustic velocity method, and (iii) presence of unknown mineralogy and bulk modulus terms in the acoustic velocity method. NMR-density porosity-derived gas hydrate saturation based on the analysis of the transverse magnetization relaxation time (T2) is considered the most precise method, but acquisition of NMR-based logs is limited at relatively recent drilled sites; additionally, its use in conventional oil and gas reservoirs is not that common due to higher cost and operational deployment limitations associated with acquiring NMR well-logs. This study proposes a new method that predicts gas hydrate saturation (S h ) for any well using porosity, bulk density, and compressional wave (P wave) velocity well-logs with neural network (or stochastic gradient descent regression) without any well-specific calibration and/or other aforementioned shortcomings of the existing methods. The method is developed by examining the underlying dependency between S h and different combinations of well-logs, chosen from 6 routine logs, with 12 different machine learning (ML) algorithms. The accuracy of the proposed method in predicting S h is ~ 84%, which is better than the accuracy of seismic and electrical resistivity methods (≤ 75%) per the results reported by three different studies. The robustness of the method in the specific case of permafrost-associated gas hydrates is demonstrated with well-log data from two wells drilled on the Alaska North Slope.

58 GEOSCIENCES↗
NASA NTRS2022

Tidal Insights into Rocky and Icy Bodies: An Introduction and Overview

Solid body tides provide key information on the interior structure, evolution, and origin of the planetary bodies. Our Solar system harbours a very diverse population of planetary bodies, including those composed of rock, ice, gas, or a mixture of all. While a rich arsenal of geophysical methods has been developed over several years to infer knowledge about the interior of the Earth, the inventory of tools to investigate the interiors of other Solar-system bodies remains limited. With seismic data only available for the Earth, the Moon, and Mars, geodetic measurements, including the observation of the tidal response, have become especially valuable and therefore, has played an important role in understanding the interior and history of several Solar system bodies. To use tidal response measurements as a means to obtain constraints on the interior structure of planetary bodies, appropriate understanding of the viscoelastic reaction of the materials from which the planets are formed is needed. Here, we review the fundamental aspects of the tidal modeling and the information on the present-day interior properties and evolution of several planets and moons based on studying their tidal response. We begin with an outline of the theory of viscoelasticity and tidal response. Next, we proceed by discussing the information on the tidal response and the inferred structure of Mercury, Venus, Mars and its moons, the Moon, and the largest satellites of giant planets, obtained from the analysis of the data that has been provided by space missions. We also summarise the upcoming possibilities offered by the currently planned missions.

Tides↗
NASA NTRS1997

Using Laser Altimetry to Detect Topographic Change in Long Valley Caldera, California

Long Valley caldera California, is a site of extensive volcanism, persistent seismicity, and uplift of a resurgent dome, currently at a rate of about 3 cm/year. Airborne laser altimetry was used to determine the surface topography of the region in 1993. A repeat mission occurred in 1995. Three different laser altimeters were flown, dubbed ATLAS, SLICER and RASCAL. Data processing consists of the combination of the aircraft trajectory and attitude data with the laser range, the determination of an atmospheric delay, laser pulse timing errors, laser system biases, and data geolocation to obtain the position of the laser spot on the ground. Results showed that using the ATLAS and SLICER instruments, the elevation of an overflown lake is determined to precisions of 3.3 cm and 2.9 cm from altitudes of 500 m and 3 km above the ground, and about 10 cm using the RASCAL instrument from 500 m above ground. Comparison with tide gauge data showed the laser measurements are able to resolve centimeter-level changes in the lake elevation over time. Repeat pass analysis of tracks over flat surfaces indicate no systematic biases affect the measurement procedure of the ATLAS and SLICER instruments. Comparison of GPS and laser-derived elevations of easily-identifiable features in the caldera confirm the horizontal accuracy of the measurement is within the diameter of the laser footprint, and vertical accuracy is within the error inherent in the measurement. Crossover analysis shows that the standard error of the means at track intersection points within the caldera, and dome (i.e., where zero and close to the maximum amount of uplift is expected) are about I cm, indicating elevation change at the 3 cm/year level should be detectable. We demonstrate one of the powerful advantages of scanning laser altimetry over other remote sensing techniques; the straightforward creation of precise digital elevation maps of overflown terrain. Initial comparison of the 1993-1995 data indicates uplift occurred, but filtering is required to remove vegetation effects. Although research continues to utilize the full potential of laser altimetry data, the results constitute a successful demonstration that the technique may be used to perform geodetic monitoring of surface topographic change.

Hofton, M. A.↗
NASA NTRS1997

Using Laser Altimetry to Detect Topographic Change at Long Valley Caldera, California

Long Valley caldera, California, is a site of extensive volcanism, persistent seismicity, and uplift of a resurgent dome, currently at a rate of approximately 3 cm/year. Airborne laser altimetry was used to determine the surface topography of the region in 1993. A repeat mission occurred in 1995. Three different laser altimeters were flown, dubbed ATLAS, SLICER and RASCAL. Data processing consists of the combination of the aircraft trajectory and attitude data with the laser range, the determination of an atmospheric delay, laser pulse timing errors, laser system biases, and data geolocation to obtain the position of the laser spot on the ground. Results showed that using the ATLAS and SLICER instruments, the elevation of an overflown lake is determined to precisions of 3.3 cm and 2.9 cm from altitudes of 500 m and 3 km above the ground, and approximately 10 cm using the RASCAL instrument from 500 m above ground. Comparison with tide gauge data showed the laser measurements are able to resolve centimeter-level changes in the lake elevation over time. Repeat pass analysis of tracks over flat surfaces indicate no systematic biases affect the measurement procedure of the ATLAS and SLICER instruments. Comparison of GPS and laser-derived elevations of easily-identifiable features in the caldera confirm the horizontal accuracy of the measurement is within the diameter of the laser footprint, and vertical accuracy is within the error inherent in the measurement. Crossover analysis shows that the standard error of the means at track intersection points within the caldera and dome (i.e., where zero and close to the maximum amount of uplift is expected) are about 1 cm, indicating elevation change at the 3 cm/year level should be detectable. We demonstrate one of the powerful advantages of scanning laser altimetry over other remote sensing techniques; the straightforward creation of precise digital elevation maps of overflown terrain. Initial comparison of the 1993-1995 data indicates uplift occurred, but filtering is required to remove vegetation effects. Although research continues to utilize the full potential of laser altimetry data, the results constitute a successful demonstration that the technique may be used to perform geodetic monitoring of surface topographic change.

Hofton, M. A.↗
DOE OSTI2026

Detecting and Characterizing Fracture Zones Using a Convolutional Neural Network

This project directly supports the Geothermal Technologies Office (GTO) objectives outlined in the Multi-Year Program Plan (MYPP) by advancing two key research areas: “Exploration and Characterization” and “Data, Modeling, and Analysis.” This project has successfully demonstrated a pre-drilling ability to image and characterize the distribution and connectivity of subsurface faults and fractures, key parameters for identifying permeable pathways that enable geothermal fluids to circulate and produce energy. Specifically, we developed and implemented innovative machine learning methodologies to enhance geothermal exploration. Large-scale faults were detected using a Convolutional Neural Network (CNN), while small-scale fractures were characterized using a novel Double-Beam Neural Network (DBNN). These tools have proven both technically effective and cost-efficient by reducing reliance on expensive exploratory drilling. Through collaboration with our geothermal industry partner, this research has significantly advanced techniques for identifying hidden geothermal systems and extending the productive lifespan of existing geothermal fields. We applied our methods to two geothermal fields—Soda Lake (Nevada) and Lightning Dock (New Mexico)—to identify shallow steam-charged fracture zones and characterize deep faults at depths of 1.5-2 km. The steam zone identified at the Soda Lake geothermal field showed excellent agreement with prior drilling data, validating the effectiveness of our approaches. In addition, the analysis revealed three new prospective drilling targets for further development and verification. The outcomes of this project improve our scientific understanding of geothermal reservoir behavior, enhance exploration efficiency, extend the economic life of existing geothermal plants. Ultimately, these advancements contribute to GTO’s goal of achieving more sustainable, affordable, and data-driven geothermal energy development across the United States.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Late Holocene history of Tokotan Lake (Kuril Archipelago, Russian Far East): The use of lacustrine records for paleoclimatic reconstructions from geologically dynamic settings

A high-resolution, multiproxy analysis of a core from Tokotan Lake, Urup Island, illustrates the complex paleoenvironmental signals that can be preserved in sedimentary deposits from geologically active areas, such as the Kuril Archipelago. Diatom and geochemical analyses of the Tokotan core suggest a shallowing and lessening of the lake area over the last ~3450 cal BP. Without other proxies, the record could be interpreted as showing a general decrease in precipitation over the Late Holocene. However, palynological data indicate no change in vegetation, thereby implying no significant climatic variation. Neither of these paleoclimatic interpretations agrees with Late Holocene reconstructions from the Kuril Archipelago or the Russian mainland. Paleomagnetic data imply variable sediment sources during the lake's history, indicating alterations within the Tokotan catchment (e.g., shifts in gradients or sediment loads in the inflowing streams). Such changes would be particularly important in the Shabalina River, which today has formed a partial dam across the lake's outlet. Here, the results of the Tokotan study indicate that other factors (e.g., volcanism, seismic activity, sea-level changes) rather than climate were more probable determinants of the Late Holocene lake record. Given the geologic dynamism of the region, it is perhaps not surprising that paleo-records from the Kuril Islands display an inconsistency in the numbers and/or duration of proposed Holocene climatic oscillations, underscoring the importance of multiproxy investigations in areas that are or have been geologically dynamic.

58 GEOSCIENCES↗