Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic prediction”

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 19 records

Predicting seismic amplitudes with machine learning

The accurate estimation of seismic wave amplitude is vital to precisely determine the yield, magnitude, and event discrimination possible for a given network – a critical element in nuclear explosion monitoring. This task is complicated by several factors, including but not limited to radiation pattern, scattering effects, and crustal variations, which can lead to the attenuation or amplification of amplitude along a given raypath. In this report, we explore the novel application of machine learning to the task of seismic amplitude estimation by training a simple Artificial Neural Network (ANN) on an S-wave amplitude dataset from Lai et al. (2019). Attributes from this dataset used as input to the ANN included event-station distances, station locations (latitude, longitude), event locations (latitude, longitude), event depths, event magnitudes, radiation patterns, signal-to noise ratio (SNR) measurements (average-amplitude, peak-to-trough, maximum peak), and signal periods. We find that the trained ANN predicts S-wave amplitudes with a modest tendency toward underestimating the actual values, as indicated by a linear regression between predicted and actual data (slope: 0.892, intercept: -0.651). These results suggest that an ANN can perform this task, with potential for significant improvements through improved datasets, architectures, and parameter tuning.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Predicting Seafloor Seismic Noise from Sea Surface Wave Heights and Vice Versa

Seismic hazard monitoring and global seismic tomography efforts are improved by recording seismic signals close to sources and at a variety of azimuths surrounding the source. While seismic networks have been deployed at many locations on land, seafloor stations are still sparse. Inclusion of ocean-based data would greatly improve the global coverage of seismic networks. However, the use of seafloor seismic data to complement land-based detection and characterization of events is complicated by the fact that ambient noise on the seafloor is in general much higher than noise levels on land. This noise is driven primarily by sea surface waves and tides. In this project, I define an empirical relationship between wave height and seafloor seismic displacement, and tie this relationship to the microseism noise theory of Longuet-Higgins and Ursell (1948). I find that while there appears to be a linear relationship between wave height and seismic displacement, wave height values can only be used with low confidence to predict seismic displacement and vice versa, and furthermore, factors other than wave height are significant in determining the level of seismic noise at a particular site.

58 GEOSCIENCES↗

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↗

Seismic Features Predict Ground Motions During Repeating Caldera Collapse Sequence

Abstract Applying machine learning to continuous acoustic emissions, signals previously deemed noise, from laboratory faults and slowly slipping subduction‐zone faults, demonstrates hidden signatures are emitted that describe physical details, including fault displacement and friction. However, no evidence currently exists to demonstrate that similar hidden signals occur during seismogenic stick‐slip on earthquake faults—the damaging earthquakes of most societal interest. We show that continuous seismic emissions emitted during the 2018 multi‐month caldera collapse sequence at the Kı̄lauea volcano in Hawai'i contain hidden signatures characterizing the earthquake cycle. Multi‐spectral data features extracted from 30 s intervals of the continuous seismic emission are used to train a gradient boosted tree regression model to predict the GNSS‐derived contemporaneous surface displacement and time‐to‐failure of the upcoming collapse event. This striking result suggests that at least some faults emit such signals and provide a potential path to characterizing the instantaneous and future behavior of earthquake faults.

58 GEOSCIENCES↗

Projected Seismic Activity at the Tiger Stripe Fractures on Enceladus, Saturn, From an Analog Study of Tidally Modulated Icequakes Within the Ross Ice Shelf, Antarctica

Abstract Dissipation of tidal energy is expected to generate seismicity on icy‐ocean worlds; however, the distribution and timing of this seismic activity throughout an orbital cycle is not known. We used new observations from an icy‐ocean‐world analog environment on Earth to examine the relationship between tidally driven tensile stress and seismic activity within an ice shell. We investigated a pair of rifts within Antarctica's Ross Ice Shelf which are tidally stressed in a manner analogous to the orbital cycle of tidal stress experienced by Enceladus' Tiger Stripe Fractures. We found that seismic activity at the Antarctic rifts is sensitive to both the amplitude and the rate of tensile stress across the rifts. We combined these findings with calculated stress values along Enceladus' Tiger Stripe Fractures to predict seismic‐activity levels expected along the ice‐shell fractures. We predict a peak in seismicity along the four Tiger Stripe Fractures when Enceladus is 90°–120° past pericenter in its orbit around Saturn, at which point tensile stresses would reach ∼2/3 of their maximum value. We also used the magnitude distribution of icequakes along Antarctic rifts to investigate implications for the likely size of stick‐slip rupture patches along icy faults on Enceladus. Our findings predict that the Tiger Stripe Fractures should produce sustained, low‐magnitude seismic events that involve rupture along discrete portions of each fracture's total length. We predict that seismicity would fall to 50% of peak levels when stresses across the Tiger Stripe Fractures are dominantly compressional.

Olsen, Kira G.↗

Updates to the Regional Seismic Travel Time (RSTT) Model: 1. Tomography

Abstract A function of global monitoring of nuclear explosions is the development of Earth models for predicting seismic travel times for more accurate calculation of event locations. Most monitoring agencies rely on fast, distance-dependent one-dimensional (1D) Earth models to calculate seismic event locations quickly and in near real-time. RSTT (Regional Seismic Travel Time) is a seismic velocity model and computer software package that captures the major effects of three-dimensional crust and upper mantle structure on regional seismic travel times, while still allowing for fast prediction speed (milliseconds). We describe updates to the RSTT model using a refined data set of regional phases (i.e., Pn, Pg, Sn, Lg) using the Bayesloc relative relocation algorithm. The tomographic inversion shown here acts to refine the previous RSTT public model ( rstt201404um ) and displays significant features related to areas of global tectonic complexity as well as further reduction in arrival residual values. Validation of the updated RSTT model demonstrates significant reduction in median epicenter mislocation (15.3 km) using all regional phases compared to the iasp91 1D model (22.1 km) as well as to the current station correction approach used at the Comprehensive Nuclear-Test-Ban Treaty Organization International Data Centre (18.9 km).

58 GEOSCIENCES↗

Crustal Structure of the Iceland Region from Spectrally Correlated Free-air and Terrain Gravity Data

Seismic refraction studies have provided critical, but spatially restricted constraints on the structure of the Icelandic crust. To obtain a more comprehensive regional view of this tectonically complicated area, we spectrally correlated free-air gravity anomalies against computed gravity effects of the terrain for a crustal thickness model that also conforms to regional seismic and thermal constraints. Our regional crustal thickness estimates suggest thickened crust extends up to 500 km on either side of the Greenland-Scotland Ridge with the Iceland-Faeroe Ridge crust being less extended and on average 3-5 km thinner than the crust of the Greenland-Iceland Ridge. Crustal thickness estimates for Iceland range from 25-35 km in conformity with seismic predictions of a cooler, thicker crust. However, the deepening of our gravity-inferred Moho relative to seismic estimates at the thermal plume and rift zones of Iceland suggests partial melting. The amount of partial melting may range from about 8% beneath the rift zones to perhaps 20% above the plume core where mantle temperatures may be 200-400 C above normal. Beneath Iceland, areally limited regions of partial melting may also be compositionally and mechanically layered

Leftwich, T. E.↗

Crustal Structure of the Iceland Region from Spectrally Correlated Free-air and Terrain Gravity Data

Seismic refraction studies have provided critical, but spatially restricted constraints on the structure of the Icelandic crust. To obtain a more comprehensive regional view of this tectonically complicated area, we spectrally correlated free-air gravity anomalies against computed gravity effects of the terrain for a crustal thickness model that also conforms to regional seismic and thermal constraints. Our regional crustal thickness estimates suggest thickened crust extends up to 500 km on either side of the Greenland-Scotland Ridge with the Iceland-Faeroe Ridge crust being less extended and on average 3-5 km thinner than the crust of the Greenland-Iceland Ridge. Crustal thickness estimates for Iceland range from 25-35 km in conformity with seismic predictions of a cooler, thicker crust. However, the deepening of our gravity-inferred Moho relative to seismic estimates at the thermal plume and rift zones of Iceland suggests partial melting. The amount of partial melting may range from about 8% beneath the rift zones to perhaps 20% above the plume core where mantle temperatures may be 200-400 C above normal. Beneath Iceland, areally limited regions of partial melting may also be compositionally and mechanically layered and intruded. The mantle plume appears to be centered at (64.6 deg N, 17.4 deg W) near the Vatnajokull Glacier and the central Icelandic neovolcanic zones.

Leftwich, T. E.↗

Evaluation of the PhaseNet Model Applied to the IMS Seismic Network

Producing a complete and accurate set of signal detections is essential for automatically building and characterizing seismic events of interest for nuclear explosion monitoring. Signal detection algorithms have been an area of research for decades, but still produce large quantities of false detections and misidentify real signals that must be detected to produce a complete global catalog of events of interest. Deep learning methods have shown promising capabilities in effectively characterizing seismic signals for complex tasks such as identifying phase arrival times. We use the PhaseNet model, a UNet-based Neural Network, trained on local distance data from northern California to predict seismic arrivals on data from the International Monitoring System (IMS) global network. We use an analyst-curated bulletin generated from this data set to compare the performance of PhaseNet to that of the Short-Term Average/Long-Term Average (STA/LTA) algorithm. We find that PhaseNet has the potential of outperforming traditional processing methods and recommend the training of a new model with the IMS data to achieve optimal performance.

58 GEOSCIENCES↗

Integrating deep neural networks with full-waveform inversion: Reparameterization, regularization, and uncertainty quantification

Full-waveform inversion (FWI) is an accurate imaging approach for modeling the velocity structure by minimizing the misfit between recorded and predicted seismic waveforms. However, the strong nonlinearity of FWI resulting from fitting oscillatory waveforms can trap the optimization in local minima. We have adopted a neural-network-based full-waveform inversion (NNFWI) method that integrates deep neural networks with FWI by representing the velocity model with a generative neural network. Neural networks can naturally introduce spatial correlations as regularization to the generated velocity model, which suppresses noise in the gradients and mitigates local minima. Furthermore, the velocity model generated by neural networks is input to the same partial differential equation (PDE) solvers used in conventional FWI. The gradients of the neural networks and PDEs are calculated using automatic differentiation, which back propagates gradients through the acoustic PDEs and neural network layers to update the weights of the generative neural network. Experiments on 1D velocity models, the Marmousi model, and the 2004 BP model determine that NNFWI can mitigate local minima, especially for imaging high-contrast features such as salt bodies, and it significantly improves the inversion in the presence of noise. Adding dropout layers to the neural network model also allows analyzing the uncertainty of the inversion results through Monte Carlo dropout. NNFWI opens a new pathway to combine deep learning and FWI for exploiting the characteristics of deep neural networks and the high accuracy of PDE solvers. Because NNFWI does not require extra training data and optimization loops, it provides an attractive and straightforward alternative to conventional FWI.

58 GEOSCIENCES↗

Seismicity as a guide to global tectonics and earthquake prediction.

From seismicity studies, evidence is presented for several aspects of plate-tectonic theory, including ideas of sea-floor spreading, transform faulting and underthrusting of the lithosphere in island arcs. Recent advances in seismic instrumentation, the use of computers in earthquake location, and the installation of local networks of instruments are shown to have vastly increased the data available for seismicity studies. It is pointed out that most of the world's earthquakes are located in very narrow zones along active plate margins and are intimately related to global processes in an extremely coherent manner. Important areas of uncertainty calling for further studies are also pointed out.

Sykes, L. R.↗

A prediction of Mars seismicity from surface faulting

The shallow seismicity of Mars has been estimated by measurement of the total slip on faults visible on the surface of the planet throughout geologic time. Seismicity was calibrated with estimates based on surface structures on the moon and measured lunar seismicity that includes the entire seismogenic lithosphere. Results indicate that Mars is seismically active today, with a sufficient number of detectable marsquakes to allow seismic investigations of its interior.

Golombek, Matthew P.↗