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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Enhanced Microseismicity During Production Pumping Cessation at the San Emidio Geothermal Field (Nevada, USA) in December 2016

Abstract Tectonic activity, geothermal fluids, and microseismic events (MSEs) tend to occur in similar locations as a result of spatiotemporal changes in the subsurface stress state. To quantify this association, we analyze data from a dense seismic array deployed at the San Emidio geothermal field, Nevada for 1 week in December 2016 to coincide with a 19.45‐hr shutdown of all injection and production pumping operations. 123 MSEs were detected, of which 101 occurred during the shutdown. The spatial association of the MSEs with the production wells suggests a causal relationship between the production cessation and the MSEs. Here we performed a detailed analysis to investigate reservoir material properties, distribution of seismically activated faults, and local stress state. We determined the hypocenters, magnitudes, and focal mechanisms for the MSEs, P‐wave tomographic velocity model, and local stress tensor. The results show that most MSEs occurred near the production wells. Magnitudes fall between −2.2 and 0.0 with larger events located closer to the production wells. Most MSEs occurred within a westward‐dipping normal fault zone in the reservoir associated with anomalously low P‐wave velocity values. The focal mechanism and stress inversion results show predominantly normal faulting with the maximum horizontal stress oriented north‐south. We suggest that the MSEs during shutdown were triggered on pre‐existing, small‐scale, critically stressed fault patches in the reservoir as the pore pressure increased around the production wells when the production pumping ceased. We interpret the larger MSE magnitudes closer to the production wells as a result of higher pore pressure increase.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Fracture Network Localization Preceding Catastrophic Failure in Triaxial Compression Experiments on Rocks

We quantify the spatial distribution of fracture networks throughout six in situ X-ray tomography triaxial compression experiments on crystalline rocks at confining stresses of 5–35 MPa in order to quantify how fracture development controls the final macroscopic failure of the rock, a process analogous to those that control geohazards such as earthquakes and landslides. Tracking the proportion of the cumulative volume of fractures with volumes >90th percentile to the total fracture volume, ∑ v 90 / v t o t indicates that the fracture networks tend to increase in localization toward these largest fractures for up to 80% of the applied differential stress. The evolution of this metric also matches the evolution of the Gini coefficient, which measures the deviation of a population from uniformity. These results are consistent with observations of localizing low magnitude seismicity before large earthquakes in southern California. In both this analysis and the present work, phases of delocalization interrupt the general increase in localization preceding catastrophic failure, indicating that delocalization does not necessarily indicate a reduction of seismic hazard. However, the proportion of the maximum fracture volume to the total fracture volume does not increase monotonically. Experiments with higher confining stress tend to experience greater localization. To further quantify localization, we compare the geometry of the largest fractures, with volumes >90th percentile, to the best fit plane through these fractures immediately preceding failure. The r 2 scores and the mean distance of the fractures to the plane indicate greater localization in monzonite than in granite. The smaller mean mineral diameter and lower confining stress in the granite experiments may contribute to this result. Tracking these various metrics of localization reveals a close association between macroscopic yielding and the acceleration of fracture network localization. Near yielding, ∑ v 90 / v t o t and the Gini coefficient increase while the mean distance to the final failure plane decreases. Macroscopic yielding thus occurs when the rate of fracture network localization increases.

58 GEOSCIENCES↗
NASA NTRS1985

Evolution of basin and range structure in the Ruby Mountains and vicinity, Nevada

Results from various age dating techniques, seismic reflection profiling hydrocarbon maturation studies, and structural analysis were used to evaluate the Cenozoic deformation in the Ruby Mountains and adjoining ranges (pinyon Range and Cortez Range) in Elko and Eureka Counties, Nevada. Age dating techniques used include potassium-argon ages of biotites from granites published by Kistler et al. (1981) and fission track ages from apatite and zircon. Fission track ages from apatite reflect a closing temperature of 100 plus or minus 20 deg C. Zircon fission track ages reflect a closing temperature of 175 plus or minus 25 deg C and potassium-argon ages from brotite reflect a closing temperature of 250 plus or minus 30 deg C. Thus these results allow a reasonably precise tracking of the evolution of the ranges during the Cenozoic. Seismic reflection data are available from Huntington Valley. Access to seismic reflection data directly to the west of the Harrison Pass Pluton in the central Ruby Mountains was obtained. In addition results are available from several deep exploration holes in Huntington Valley.

Blackwell, D. D.↗
DOE OSTI2024

Integrating Atmospheric Specifications into Seismoacoustic Event Localization

This report investigates the integration of infrasound and seismic data to improve event localization accuracy, specifically focusing on a surface explosion at the Utah Training and Testing Range (UTTR). Utilizing the Seismoacoustic Bayesian Event Locator (SABEL) framework, we incorporated atmospheric specifications derived from Ground to Space (G2S) profiles to enhance celerity-range priors. Our analysis revealed that while the combination of infrasound and seismic observations significantly reduced localization uncertainty, challenges remained, particularly with returns at distances less than 200 km from the source and the influence of specific observations on location estimates. The results indicate that broader celerity distributions, such as those from Blom et al. (2020), facilitate better alignment with ground truth locations compared to narrower models. Overall, this work demonstrates the promise of seismoacoustic approaches in refining event localization and highlights the need for further exploration of celerity-range models to ensure reliable outcomes.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2023

Characteristics and selection of near-fault simulated earthquake ground motions for nonlinear analysis of buildings

Earthquake-induced ground shaking near rupturing faults is highly sensitive to the rupture characteristics, seismic wave propagation patterns and site conditions, and field recordings of near-fault shaking are relatively sparse. These challenges complicate the assessment of the seismic performance of near-fault structures. A common approach to representing near-fault ground motion in engineering analysis is to explicitly consider and select records with strong directivity pulses (pulse records). We use three-dimensional high-resolution physics-based earthquake simulations to test this approach in the context of scenario-based ground motion record selection, and to study the important characteristics of near-fault ground shaking. We highlight the deficiencies associated with classifying near-fault simulated records as “pulse” or “non-pulse,” based on the presence of a single dominating pulse in the velocity time history. We show that this approach is inadequate for characterizing near-fault shaking on soft soils which can be dominated by both forward rupture directivity and basin amplification effects. We conduct ground motion selection experiments for the analysis of near-fault structures with and without explicit classification of the pulse features in the records, and evaluate the bias in the predicted structural demands. We find that the maximum interstory drift demands on building structures imposed by unscaled site-specific simulated ground motion records selected based on relevant spectral shape features are not sensitive to the classification of records as pulse/non-pulse. Therefore, with regard to predicting the maximum interstory drifts in near-fault buildings, we do not find justification for the binary pulse classification of near-fault records.

42 ENGINEERING↗
DOE OSTI2024

Accelerating uncertainty quantification in incremental dynamic analysis using dimension reduction-based surrogate modeling

We propose a surrogate modeling framework based on dimension reduction to facilitate the quantification of seismic risk of structural systems in performance-based earthquake engineering. The framework adopts incremental dynamic analysis (IDA) for addressing hazard variability, and promotes significant computational efficiency improvement for propagating epistemic uncertainties associated with the structural models. It utilizes both linear and nonlinear dimension reduction approaches, equipped with inverse mappings, to learn a functional between the input parameter space (e.g., the epistemic uncertainties of the structure) to the high-dimensional output space created through the IDA implementation across different ground motions and seismic intensity levels. Polynomial chaos expansion is adopted as the surrogate model to learn this functional in the reduced space. A nine-story steel moment-resisting frame with uncertain structural properties is used as a testbed. Furthermore, we select the seismic fragility curves as a measure of the structure’s seismic performance, since it provides an estimate of the probability of entering specified damage states for given levels of ground shaking.

42 ENGINEERING↗
DOE OSTI2020

Clustering earthquake signals and background noises in continuous seismic data with unsupervised deep learning

The continuously growing amount of seismic data collected worldwide is outpacing our abilities for analysis, since to date, such datasets have been analyzed in a human-expert intensive, supervised fashion. Moreover, analyses that are conducted can be strongly biased by the standard models employed by seismologists. In response to both of these challenges, we develop a new unsupervised machine learning framework for detecting and clustering seismic signals in continuous seismic records. Our approach combines a deep scattering network and a Gaussian mixture model to cluster seismic signal segments and detect novel structures. To illustrate the power of the framework, we analyze seismic data acquired during the June 2017 Nuugaatsiaq, Greenland landslide. We demonstrate the blind detection and recovery of the repeating precursory seismicity that was recorded before the main landslide rupture, which suggests that our approach could lead to more informative forecasting of the seismic activity in seismogenic areas.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Multiscale Data-Driven Seismic Full-Waveform Inversion With Field Data Study

Seismic full-waveform inversion (FWI), which uses iterative methods to estimate high-resolution subsurface models from seismograms, is a powerful imaging technique in exploration geophysics. In recent years, the computational cost of FWI has grown exponentially due to the increasing size and resolution of seismic data. Moreover, it is a nonconvex problem and can encounter local minima due to the limited accuracy of the initial velocity models or the absence of low frequencies in the measurements. To overcome these computational issues, we develop a multiscale data-driven FWI method based on fully convolutional networks (FCNs). In preparing the training data, we first develop a real-time style transform method to create a large set of synthetic subsurface velocity models from natural images. We then develop two convolutional neural networks with encoder-decoder structures to reconstruct the low- and high-frequency components of the subsurface velocity models, separately. To validate the performance of our data-driven inversion method and the effectiveness of the synthesized training set, we compare it with conventional physics-based waveform inversion approaches using both synthetic and field data. Finally, these numerical results demonstrate that, once our model is fully trained, it can significantly reduce the computation time and yield more accurate subsurface velocity models in comparison with conventional FWI.

58 GEOSCIENCES↗
DOE OSTI2020

Seismic monitoring and high-frequency noise using arrays in the Arabian Peninsula

Seismic arrays provide useful tools for regional seismic monitoring. Two small-aperture, regionally-oriented arrays, QWAR and HQAR, were deployed in Saudi Arabia and Oman in 2012 and 2016. We utilize a time-domain beampacking method, similar to frequency-wavenumber analysis, to examine the performance of the arrays in terms of slowness and azimuthal bias and event detection capabilities. Additionally, we investigate persistent ambient noise sources recorded by the arrays. We find that the arrays provide slowness vectors with biases comparable with similar-sized arrays of the International Monitoring System in other locations around the world. At QWAR, regional events of magnitude 3.0 and above are detected a majority of the time, and as magnitudes increase to 4.0 and above, the detection rate is greater than 82%. Finally, strong noise generation, primarily with slownesses characteristic of Lg waves, is found in the directions of the northern Arabian/Persian Gulf and southern Red Sea and may be a factor in event detection capabilities.

58 GEOSCIENCES↗
DOE OSTI2023

Characterizing the seismic response of a molten salt nuclear reactor

Abstract One Generation IV nuclear reactor, which uses a fluoride salt as a coolant, graphite reflector blocks as a moderator, and circulating buoyant TRISO pebbles as fuel is at an advanced stage of development. To characterize the seismic behavior of components of this reactor, validate numerical models for analysis, and develop recommendations for design, a set of earthquake‐simulator experiments on a scaled model of the reactor vessel and its internals was executed on a six‐degree‐of‐freedom earthquake simulator. The model was seismically isolated at its base using two types of spherical sliding bearings. The scaled model involved representations of the prototype reactor vessel, core barrel, reflector blocks, coolant, and spherical fuel pebbles. The material and geometric properties of different test components were selected based on a dynamic similitude scaling analysis and an approximate length scale of 0.4. Four sets of three‐component earthquake motions were used as inputs for testing. Instrumentation on the test specimen recorded the dynamic responses of the outer vessel, core barrel, and reflector‐block assembly, the hydrodynamic responses (sloshing and hydrodynamic pressure) of the liquid coolant, pebble consolidation under earthquake shaking, and the behavior of the isolation systems. This paper describes the design of the experiments and presents key results from the tests. The dynamic responses of the outer vessel, core barrel, and the reflector blocks revealed that the components responded as a unit for the intense shaking used in the experiments. The sloshing response of the fluid in a thin annulus near the perimeter of the vessel was heavily damped. The change in the packing fraction of the pebble bed under repeated, intense 3D earthquake shaking was less than 3%. Seismically isolating the vessel substantially reduced demands on its internal components.

Engineering↗
DOE OSTI2024

Empowering Machine Learning Forecasting of Labquake Using Event‐Based Features and Clustering Characteristics

Abstract Following recent advances of machine learning (ML), we present a novel approach to extract spatiotemporal seismo‐mechanical features from Acoustic Emission (AE) catalogs to empower ML‐based forecasting. The AE data were recorded during laboratory stick‐slip experiments on granite samples cut by rough faults. Based on the features computed for a past time window, a random forest (RF) classifier is used to forecast the occurrence of a large magnitude event ( M AE > 3.5) in the next time window. Event‐based features allow us to associate informative time‐space characteristics to each feature and nearest‐neighbor clustering analysis enables us to separate background and clustered seismicity and train individual models. The results show that the separation of AEs enhances the forecasting accuracy from 73.2% for the entire catalog up to 82.1% and 89.0% if background and clustered events are used separately. The presented new approach may be upscaled for applications to forecast tectonic earthquakes.

Karimpouli, Sadegh↗
DOE OSTI2021

WHOLESCALE: Seismic Survey Data from San Emidio Nevada 2021

This dataset includes raw and processed seismic data from the 2021 seismic survey at the San Emidio geothermal field in Nevada. In April and May 2021, 37 tri-axial short period seismographs were deployed in a 1.8km diameter cluster centered on 40.367278 N, 119.409019 W. The first data record started at 2021-04-06T07:09:10Z UTC and the last record ended at 2021-05-11T02:58:52Z UTC. The 37 stations included 29 SmartSolo IGU-16HR 3C all-in-one 5 Hz seismographs and 8 DataCube seismographs with 4.5 Hz HGS HG-6(B coil) tri-axial geophones. The raw format (level 0) data includes 353 GB of 500 sps SmartSolo data in native DLD format, 113 GB of 400 sps DataCube data in native DataCube format, and 3.4 GB of GPS data collected during the RTK GPS survey. The SAC data (level 1) files include 564 GB of hourly SAC files. The experiment was designed to monitor seismic activity before, during, and after the planned three-day plant maintenance shutdown April 19-21, 2021. The pumping stopped at 2021/04/19 12:51:45 UTC and resumed about 2021/04/21 21:00:00 UTC. The dataset is hosted in an AWS data lake, with an associated GDR data set providing the metadata. Links to both of these resources are included below. Additionally, this collection features data and metadata from a 2016 seismic survey at the same site for comparative analysis.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Waveform Simulation Framework: User Manual with Tutorials

This manuscript describes the Waveform Simulation Framework (WSF), a Python-based framework that provides a unified, programmable interface for generating synthetic seismograms for applications such as seismic array design, method development, and special event analysis. WSF standardizes how users define sources, receivers, and velocity models while abstracting simulator-specific configuration details, enabling workflows that are largely independent of the underlying numerical engine. The document provides installation guidance and tutorial-driven examples for three WSF simulator wrappers—WSF PyFK, WSF SW4, and WSF SPECFEM2D—illustrating end-to-end workflows from forward waveform simulation to common post-processing tasks (e.g., visualization and backprojection) using consistent data products (e.g., ObsPy Stream objects and SAC files).

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

Structure of the lunar crust at highland site Apollo Station 16

The seismic crustal structure of the moon is known in the region of Oceanus Procellarum from the analysis of artificial impact data. To extend this knowledge we have used data from natural lunar seismic events to search for secondary seismic wave arrivals in the form of peg-leg multiples caused by reflections at crustal interfaces and converted arrivals caused by refractions at crustal interfaces. A polarization filter has been applied to the data to enhance the rectilinear particle motion expected for the onset of these body wave arrivals in the scattered coda of lunar seismograms. The results of this work tentatively indicate that the highland site at station 16 has a 75 km thick crust with an intermediate 20 km interface, compared to the 60 km crust and 20 km interface in Oceanus Procellarum, a mare region. Since the 20 km upper crust appears to exist at both highland and mare sites, it probably does not represent a mare basalt layer but rather a more general feature of the crust. Crustal thickness may partially control elevation by isostasy.

Goins, N. R.↗
NASA NTRS2003

Exploring Asteroid Interiors: The Deep Interior Mission Concept

Deep Interior is a mission to determine the geophysical properties of near-Earth objects, including the first volumetric image of the interior of an asteroid. Radio reflection tomography will image the 3D distribution of complex dielectric properties within the ~1 km rendezvous target and hence map structural, density or compositional variations. Laser altimetry and visible imaging will provide high-resolution surface topography. Smart surface pods culminating in blast experiments, imaged by the high frame rate camera and scanned by lidar, will characterize active mechanical behavior and structure of surface materials, expose unweathered surface for NIR analysis, and may enable some characterization of bulk seismic response. Multiple flybys en route to this target will characterize a diversity of asteroids, probing their interiors with non-tomographic radar reflectance experiments. Deep Interior is a natural follow-up to the NEARShoemaker mission and will provide essential guidance for future in situ asteroid and comet exploration. While our goal is to learn the interior geology of small bodies and how their surfaces behave, the resulting science will enable pragmatic technologies required of hazard mitigation and resource utilization.

Asphaug, E.↗
NASA NTRS2014

Broadscale Postseismic Gravity Change Following the 2011 Tohoku-Oki Earthquake and Implication for Deformation by Viscoelastic Relaxation and Afterslip

The analysis of GRACE gravity data revealed post-seismic gravity increase by 6 micro-Gal over a 500 km scale within a couple of years after the 2011 Tohoku-Oki earthquake, which is nearly 40-50% of the co-seismic gravity change. It originates mostly from changes in the isotropic component corresponding to the M(sub rr) moment tensor element. The exponential decay with rapid change in a year and gradual change afterward is a characteristic temporal pattern. Both viscoelastic relaxation and afterslip models produce reasonable agreement with the GRACE free-air gravity observation, while their Bouguer gravity patterns and seafloor vertical deformations are distinctly different. The post-seismic gravity variation is best modeled by the bi-viscous relaxation with a transient and steady state viscosity of 10(exp 18) and 10(exp 19) Pa s, respectively, for the asthenosphere. Our calculated higher-resolution viscoelastic relaxation model, underlying the partially ruptured elastic lithosphere, yields the localized post-seismic subsidence above the hypocenter reported from the GPS-acoustic seafloor surveying.

GRACE↗
DOE OSTI2021

Intelligent Monitoring Systems and Advanced Well Integrity and Mitigation

Long-term seismic monitoring of carbon capture and storage projects is needed to verify that the injected gas is safely stored in the subsurface until permanence can be assured. Conventional surface seismic monitoring techniques are usually expensive, require highly invasive surface operations, and need significant time investments on the part of personnel for both the field effort and processing the acquired data. For these reasons, permanent reservoir monitoring technologies are preferred, as they can offer a cost-effective solution for long-term monitoring. As part of the monitoring program of the Archer Daniels Midland’s large-scale injection of CO 2 in Decatur, Illinois, USA, a continuous seismic monitoring array was installed using a combination of surface orbital vibrator (SOV) sources and fiber-optic cables for distributed acoustic sensing (DAS) acquisition with the objective to build a continuous monitoring array. The aim of the presented project was to build a monitoring array and platform that integrates real-time seismic data with conventional data streams and provides continuous data analysis using dynamic computational models to deliver a comprehensive real-time assessment of subsurface conditions. It is in this context that the Intelligent Monitoring Systems and Advanced Well Integrity and Mitigation project was proposed with the objective to develop an integrated architecture that utilizes a permanent seismic monitoring network, combines the real-time geophysical and process data with reservoir flow and geomechanical models to create a comprehensive monitoring, visualization, and control system that delivers critical information for process surveillance and optimization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Coupling geophysical, geological, geochemical and mineralogical assessments to examine preferential contaminant transport pathways in interbedded fractured bedrock

This study shows that a multi-faceted approach, combining borehole geophysical logging and surface seismic P-wave first-arrival tomography with confirmatory coring, well installation, and chemical and mineralogical analysis, is effective for identifying difficult-to-locate preferential contaminant transport pathways in deeper fractured bedrock. Seismic tomography detected porous 10–20 m wide elongated fractured conduits that allow acidic groundwater contaminated with uranium (U) and nitrate (NO 3 - ) to migrate within interbedded shale-limestone bedrock over 1000 m from a former disposal facility (S-3 Ponds site) located at the DOE Y-12 National Security Complex in Tennessee (USA). Conventional drilling techniques would easily miss these conduits because they are oriented parallel with fractured bedding planes. Synchrotron analysis of aquifer solids revealed that > 95 % of the U is hexavalent (U VI ). This uranyl (UO 2 2+ ) species is coordinated with carbonate, iron oxide, silicate and phosphate minerals within cemented white to yellow precipitates, which contain U concentrations as high as ∼21.6 % by weight fraction. Identifying the presence of these mineral phases, enables a further understanding of the potential effectiveness of remediation actions. The combination of methodologies presented here can also be applied to other explorations, such as the detection of water supply.

Fate and transport↗