Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic methods”

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 55 records · Page 3

Seismicity, shear failure and modes of deformation in deep subduction zones

The joint hypocentral determination method is used to relocate deep seismicity reported in the International Seismological Center catalog for earthquakes deeper than 400 km in the Honshu, Bonin, Mariannas, Java, Banda, and South America subduction zones. Each deep seismic zone is found to display planar features of seismicity parallel to the Harvard centroid-moment tensor nodal planes, which are identified as planes of shear failure. The sense of displacement on these planes is one of resistance to deeper penetration.

Lundgren, Paul R.↗

Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties

This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.

15 GEOTHERMAL ENERGY↗

On the relative scattering of P- and S-waves

Using a single-scattering approximation, equations for the scattering attenuation coefficients of P-body and S-body waves are derived. The results are discussed in the light of the energy-renormalization approaches of Wu (1980, 1982) and Sato (1982) to seismic wave scattering. Practical methods for calculating the scattering attenuation coefficients for various earth models are emphasized. The conversions of P-waves to S-waves and S-waves to P-waves are included in the theory. The earth models are assumed to be randomly inhomogeneous, with their properties known only through their average-wavenumber power spectra. The power spectra are approximated with piecewise constant functions, each segment of which contributes to the net frequency-dependent scattering attenuation coefficient. The smallest and largest wavenumbers of a segment can be plotted along with the wavevectors of the incident and scattered waves on a wavenumber diagram. This diagram gives a geometric interpretation for the frequency behavior associated with each spectral segment, including a transition peak that is due entirely to the wavenumber limits of the segment. For regions of the earth where the inhomogeneity spectra are concentrated in a band of wavenumbers, it should be possible to observe such a peak in the apparent attenuation of seismic waves. Both the frequency and distance limits on the accuracy of the theoretical results are given.

Malin, P. E.↗

Regularization by denoising diffusion models for solving inverse PDE problems with application to full waveform inversion

Partial differential equation (PDE)-governed inverse problems are fundamental across various scientific and engineering applications; yet they face significant challenges due to nonlinearity, ill-posedness, and sensitivity to noise. Here, we introduce a computational framework, regularization by denoising using diffusion models for partial differential equations (RED-DiffEq), by integrating physics-driven inversion and data-driven learning. RED-DiffEq leverages pretrained diffusion models as a regularization mechanism for PDE-governed inverse problems. We apply RED-DiffEq to solve the full waveform inversion problem in geophysics, a challenging seismic imaging technique that seeks to reconstruct high-resolution subsurface velocity models from seismic measurement data. Our method shows enhanced accuracy and robustness compared to benchmark methods. Additionally, it exhibits strong generalization and domain decomposition capacity, enabling the inversion of more complex velocity models with larger domains than those used in training the diffusion model. Our framework can also be directly applied to diverse PDE-governed inverse problems.

Shan, Siming [Yale University, New Haven, CT (Unit↗

A statistical technique for processing radio interferometer data

The output of a radio interferometer is the Fourier transform of the object under investigation. Due to the limited coverage of the Fourier plane, the reconstruction of the image of the source is blurred by the beam of the synthesized array. A maximum-likelihood processing technique is described which uses the statistical properties of the received noise-like signals. This technique has been used extensively in the processing of large-aperture seismic arrays. This inversion method results in a synthesized beam that is more uniform, has lower sidelobes, and higher resolution than the normal Fourier transform methods. The maximum-likelihood method algorithm was applied successfully to very long baseline and short baseline interferometric data.

Papadopoulos, G. D.↗

Using GPS to Detect Imminent Tsunamis

A promising method of detecting imminent tsunamis and estimating their destructive potential involves the use of Global Positioning System (GPS) data in addition to seismic data. Application of the method is expected to increase the reliability of global tsunami-warning systems, making it possible to save lives while reducing the incidence of false alarms. Tsunamis kill people every year. The 2004 Indian Ocean tsunami killed about 230,000 people. The magnitude of an earthquake is not always a reliable indication of the destructive potential of a tsunami. The 2004 Indian Ocean quake generated a huge tsunami, while the 2005 Nias (Indonesia) quake did not, even though both were initially estimated to be of the similar magnitude. Between 2005 and 2007, five false tsunami alarms were issued worldwide. Such alarms result in negative societal and economic effects. GPS stations can detect ground motions of earthquakes in real time, as frequently as every few seconds. In the present method, the epicenter of an earthquake is located by use of data from seismometers, then data from coastal GPS stations near the epicenter are used to infer sea-floor displacements that precede a tsunami. The displacement data are used in conjunction with local topographical data and an advanced theory to quantify the destructive potential of a tsunami on a new tsunami scale, based on the GPS-derived tsunami energy, much like the Richter Scale used for earthquakes. An important element of the derivation of the advanced theory was recognition that horizontal sea-floor motions contribute much more to generation of tsunamis than previously believed. The method produces a reliable estimate of the destructive potential of a tsunami within minutes typically, well before the tsunami reaches coastal areas. The viability of the method was demonstrated in computational tests in which the method yielded accurate representations of three historical tsunamis for which well-documented ground-motion measurements were available. Development of a global tsunami-warning system utilizing an expanded network of coastal GPS stations was under consideration at the time of reporting the information for this article.

Song, Y. Tony↗

Seismic response of vertical dry storage casks under three-dimensional earthquake motions

Ensuring the long-term seismic safety of dry storage casks (DSCs) is becoming increasingly critical as these systems evolve from temporary to de facto permanent repositories for spent nuclear fuels. Traditional seismic soil–structure interaction (SSI) assessment methods use one-dimensional deconvolution or simplified boundary conditions to model incident waves. Although computationally appealing, simplifying assumptions may alter the seismic risk by neglecting the full complexity of three-dimensional (3D) wave propagation effects. To address this challenge, this paper introduces a novel high-fidelity computational framework that leverages the Domain Reduction Method (DRM) with perfectly matched layers (PML) to accurately transfer complex, 3D seismic wavefields from regional-scale fault-rupture simulations into local-scale finite element models of DSCs. Using broadband, physics-based ground motions from a generic M w 7.0 strike-slip event, both single-cask and multi-cask configurations were investigated under near- and far-field conditions. Emphasis is placed on capturing complex SSI, spatial variability in the ground motion, and nonlinear phenomena such as cask rocking and sliding. Numerical results demonstrate that near-field conditions, where forward directivity and fling-step effects dominate, lead to significantly higher DSC rocking and sliding. Far-field cases, by contrast, generally exhibit modest responses. Incorporating SSI tends to amplify or alter DSC response spectra and introduce response variability, which underscores the need for site-specific evaluations and robust modeling approaches to ensure the seismic integrity of DSCs in interim spent fuel storage installations.

Das, Tonmoy↗

Seismicity-constrained fault detection and characterization with a multitask machine learning model

Geological fault detection and characterization are crucial for understanding subsurface dynamics across scales. While methods for fault delineation based on either seismicity location analysis or seismic image reflector discontinuity are well-established, a systematic approach that integrates both data types remains absent. We develop a novel machine learning model that unifies seismic reflector images and seismicity location information to automatically identify geological faults and characterize their geometrical properties. The model encodes a seismic image and a seismicity location image separately, and fuses the encoded features with a spatial-channel attention fusion module to improve the learning of important features in both inputs. We design an automated strategy to generate high-quality synthetic training data and labels. To improve the realism of the seismicity location image, we include random seismicity noise and missing seismicity location associated with some of the faults. We validate the model’s efficacy and accuracy using synthetic data examples and two field data examples. Moreover, we show that fine-tuning the trained model with a small, domain-specific dataset enhances its fidelity for field data applications. The results demonstrate that integrating seismicity location and seismic images into a unified framework allows the end-to-end neural network to achieve higher fidelity and accuracy in delineating subsurface faults and their geometrical properties compared with image-only fault detection methods. Our approach offers an adaptive data-driven tool for geological fault characterization and seismic hazard mitigation, bridging the gap between seismicity location and image-based fault detection methods.

58 GEOSCIENCES↗

Strong earthquakes, novae and cosmic ray environment

Observations about the relationship between seismic activity and astronomical phenomena are discussed. First, after investigating the seismic data (magnitude 7.0 and over) with the method of superposed epochs it is found that world seismicity evidently increased after the occurring of novae with apparent magnitude brighter than 2.2. Second, a great many earthquakes of magnitude 7.0 and over occurred in the 13th month after two of the largest ground level solar cosmic ray events (GLEs). The causes of three high level phenomena of global seismic activity in 1918-1965 can be related to these, and it is suggested that according to the information of large GLE or bright nova predictions of the times of global intense seismic activity can be made.

Yu, Z. D.↗

Underground hydrogen storage leakage detection and characterization based on machine learning of sparse seismic data

Underground hydrogen storage (UHS) is considered as a scalable approach for massive storage and seasonal extraction of hydrogen (H 2 ). Although conventional leakage detection and characterization methods based on time-lapse seismic imaging and inversion generally apply to H 2 leakage detection problem, a high-fidelity yet cost effective geophysics approach is still missing to reliably inform leakage location and properties based on very sparse data. In response, we develop a novel supervised machine learning method to detect and characterize H 2 leakage from UHS. The input to our neural network are sparse time-lapse seismic waveforms, while the output from the neural network includes the spatial location and physical properties of a H 2 leakage. Here, we generate high-quality time-lapse waveforms using the elastic-wave equations to train the neural network. We train and validate our machine learning model and find that it attains high accuracy in using extremely sparse time-lapse seismic data to detect and characterize H 2 leakage. Our investigation is the first systematic study that focuses on applying machine learning to subsurface H 2 leakage detection and characterization and could potentially serve as a cost-effective geophysical tool for underground hydrogen leakage detection and characterization with high fidelity.

08 HYDROGEN↗

Exploring Data Set Bias and Decision Support with Predictive Uncertainty Through Bayesian Approximations and Convolutional Neural Networks

Individual seismic catalogs can contain multiscale observations from fault level to global scales and associated waveforms from discrete events reflect crustal structure across many different scales and locations. Seismic network aperture, geographic location, and observation distance may not provide informative guidance or intuition on how different catalogs will behave across models trained under different conditions. We rely on uncertainty to provide guardrails for when to trust model decisions, but understanding when our uncertainty is trustworthy is an open challenge. Here, in this work, we explore Bayesian approximation methods for assigning predictive uncertainty in seismic event classification problems. We find that computationally expensive Bayesian approximations do not outperform simple ensemble methods. We also find that when exploiting multiple seismic event catalogs, joint training with data from all the catalogs combined with Bayesian approximations and supervised training for classification can obscure bias and result in less robust uncertainty while also not providing substantial performance benefits compared to training individual models for each catalog.

58 GEOSCIENCES↗

SIIOS in Alaska: Testing an "In-Vault" Option for a Europa Lander Seismometer Experiment

The icy moons of Europa and Enceladus are thought to have global subsurface oceans in contact with mineral-rich silicate interiors, likely providing the three ingredients needed for life as we know it: liquid water, essential chemicals, and a source of energy. The possibility of life forming in their subsurface oceans relies in part on transfer of oxidants from the irradiated ice surface to the sheltered ocean below. Constraining the mechanisms and location of material exchange between the ice surface, the ice shell, and the subsurface ocean, however, is not possible without knowledge of ice thickness and liquid water depths. In a future lander-based experiment seismic measurements will be a key geophysical tool for obtaining this critical knowledge. The Seismometer to Investigate Ice and Ocean Structure (SIIOS) field-tests flight-ready technologies and develops the analytical methods necessary to make a seismic study of Europa and Enceladus a reality. We have been performing small-array seismology with a flight-candidate sensor in analog environments that exploit passive sources. Determining the depth to a subsurface ocean and any intermediate bodies of water is a priority for Ocean Worlds missions as it allows assessment of the habitability of these worlds and provides vital information for evaluating the spacecraft technologies required to access their oceans.

Bray, Veronica J.↗

Determination of design and operation parameters for upper atmospheric research instrumentation to yield optimum resolution with deconvolution, appendix 2

This thesis reviews the technique established to clear channels in the Power Spectral Estimate by applying linear combinations of well known window functions to the autocorrelation function. The need for windowing the auto correlation function is due to the fact that the true auto correlation is not generally used to obtain the Power Spectral Estimate. When applied, the windows serve to reduce the effect that modifies the auto correlation by truncating the data and possibly the autocorrelation has on the Power Spectral Estimate. It has been shown in previous work that a single channel has been cleared, allowing for the detection of a small peak in the presence of a large peak in the Power Spectral Estimate. The utility of this method is dependent on the robustness of it on different input situations. We extend the analysis in this paper, to include clearing up to three channels. We examine the relative positions of the spikes to each other and also the effect of taking different percentages of lags of the auto correlation in the Power Spectral Estimate. This method could have application wherever the Power Spectrum is used. An example of this is beam forming for source location, where a small target can be located next to a large target. Other possibilities extend into seismic data processing. As the method becomes more automated other applications may present themselves.

Ioup, George E.↗

Modeling of Seismic Waves Through Geologic Metamaterials

This project conducted a modeling study on seismic invisibility cloaks that render geologic targets invisible to seismic waves, using the concept of seismic metamaterials. We present a parametric numerical study on the behaviors of seismic waves through cloaks with different design parameters as well as degrees of geologic heterogeneity. In addition, a seismic cloaking strategy is proposed for a future field-scale experiment at a real-world test bed. This feasibility study will guide future field experiment designs and ultimately allow us to conduct systematic field-scale tests employing Sandia’s existing resources and field expertise. The ultimate goal is to develop methods and design parameters of seismic invisibility cloaks to protect against natural and man-made seismic waves. Seismic cloaking has potential applications in several areas of national security, energy, and natural hazard reduction.

58 GEOSCIENCES↗

Relocation of the 8 September 2023 High Atlas, Morocco, Earthquake Aftershock Sequence

The earthquake that occurred on 8 September 2023, with a magnitude of 6.8, was the most destructive earthquake event in Morocco in the past decade. This earthquake took place in the Al Haouz region, located in the western part of the High Atlas Mountain range. To better understand what caused and triggered this earthquake, the earthquake catalogs including P and S arrival times were collected from the Moroccan seismic network and combined with regional data from the International Seismological Centre. The mainshock and aftershocks were relocated by using iLoc, a state-of-the-art single-event location algorithm, and then by the multiple event location double-difference algorithm, hypoDD. The improved earthquake relocations using iLoc and the double-difference methods provide sharper lineation of seismicity and agree well with tomographic images of the earthquake zone. Finally, the seismicity distribution and the focal mechanism of the mainshock indicate that the earthquake sequence has occurred along the South Atlas fault system.

58 GEOSCIENCES↗

Technical Report on Subsurface Monitoring of the Brady Hot Spring Geothermal Site, Nevada, based upon Full Waveform Inversion

Abilities to accurately characterize the subsurface in a geothermal setting is key to assess and support production. An important element of geothermal reservoir monitoring is also the ability to investigate fluid transport within fracture network. This report focuses on improving subsurface imaging and monitoring in geothermal settings using full waveform inversion based on the adjoint method and time-lapse imaging. To assess our method, we rely on a dense seismic dataset collected in 2016 at the Brady Hot Springs geothermal site in Nevada for the DOE-funded project Poroelastic Tomography by Adjoint Inverse Modeling of Data from Seismology, Geodesy, and Hydrology. This dataset captures subsurface changes across four stages of geothermal power plant operations, which involve varying rates of fluid injection and extraction. Two velocity models were previously derived from this dataset using different methods: one based on travel times and another on sweep interferometry. Our first step is to refine these models using adjoint tomography, which has been applied successfully at global and regional-scales but is less common at the reservoir-scale. Two approaches are then explored for time-lapse analysis: directly comparing refined tomographic models from different stages or backpropagating waveform differences relative to a baseline tomographic model. The main take away is that both approaches highlight similar reservoir behaviors, but the latter approach is more computationally effective in capturing small-scale changes in subsurface properties. For this work, we leverage the use of Salvus (www.mondaic.com), an end-to-end seismic imaging solution, relying on the spectral element method to compute forward and adjoint simulations, and developed by Mondaic Ltd. It includes integrated workflow management that handles waveform and metadata, launches simulations, computes waveform misfits and adjoint sources, and iterates for model updates by nonlinear optimization.

15 GEOTHERMAL ENERGY↗

Optimum seismic structural design based on random vibration and fuzzy graded damages

This paper presents the fuzzy dynamical reliability and failure probability as well as the basic principles and the analytical method of loss assessment for nonlinear seismic steel structures. Also presented is the optimization formulation and a numerical example for double objectives, initial construction cost and expected failure loss, and dynamical reliability constraints. The earthquake ground motion is based on a stationary filtered non-white noise and the fuzzy damage grade is described by damage index.

Cheng, Franklin Y.↗