Search NASA⌕ Search

SEARCH · Search NASA

Results for “SEISMICITY”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Small-Magnitude Seismic Swarms in Central Utah (US): Interactions of Regional Tectonics, Local Structures and Hydrothermal Systems

Swarms in Central Utah are situated in the complex transition between the Basin and Range (BR) province and the Colorado Plateau. Transecting transverse structures, volcanic deposits, and hydrothermal systems complicate the extensional BR horst and graben structures and provide a multitude of plausible triggering mechanisms. Revisiting the catalog of the University of Utah Seismograph Stations (1981–2022), we analyze spatio-temporal patterns and characteristic features of seismic sequences. Swarms with alternating seismicity rates, bursts, and longer swarms with persistent moment release exhibit a remarkable diversity in temporal evolution. Swarm durations do not scale with cumulative seismic moment: swarms lasting less than 1 day can have similar cumulative seismic moments as month-long swarms. Here, we observe stationary swarms re-occurring for years (e.g., Mineral Mountains), as well as singular swarms in low-seismicity areas (e.g., activating a local structure). The swarms show a pronounced heterogeneity in triggering and driving mechanisms, observed in the detailed analysis of exemplary sequences (detections, relocations, moment tensors, waveform-based clustering, and repeater analysis). The 2022 Sevier Valley sequence activated a BR-related normal fault, the first resolved fault plane in the valley since 1983. The 2011 Circleville sequence is interpreted as a swarm triggered by mainshock-aftershock activity characterized by increasing magnitudes, changing rupture mechanisms, and a concentration of highly similar events in the second part of the sequence. By jointly discussing exemplary sequences and catalog statistics, we draw a comprehensive picture of swarm activity and its relation to geothermal and tectonic activity.

58 GEOSCIENCES↗
DOE OSTI2024

A Frequency Domain Methodology for Quantitative Evaluation of Diffuse Wavefield With Applications to Seismic Imaging

Abstract Ambient Noise Imaging (ANI) of subsurface structures relies on seismic interferometry of diffuse seismic wavefields. However, the lack of effective methods to quantify and identify highly diffuse waves hampers applications of ANI, particularly in evaluating seismic attenuation and monitoring structural changes with high temporal resolution. Conventional ANI approaches require data normalization, which effectively suppresses the non‐diffuse component with large amplitude but also results in significant loss of amplitude and phase information in the continuous seismic records. In this study, we propose a frequency domain method to quantitatively evaluate the degree of diffuseness of seismic wavefields by analyzing their statistical characteristics of modal amplitudes for stationarity and randomness. Tests on synthetic waveform and field nodal records show that the proposed method can effectively distinguish between diffuse and non‐diffuse waveforms for either single‐ or three‐component data. As an application, we identify a 60‐s‐long diffuse coda of a local M 2.2 earthquake recorded by a dense nodal array on the San Jacinto Fault Zone, and successfully extract high‐quality dispersion curve andQ‐value without performing data normalization. These results are consistent with those obtained by conventional methods that assess the correlation between coherency and the Green's function, and by modeling ballistic waves generated by road traffic. Our proposed method can advance the imaging of subsurface velocity and attenuation structures as well as monitoring temporal changes for scientific studies and engineering applications.

Geochemistry & Geophysics↗
DOE OSTI2019

Energy of injection-induced seismicity predicted from in-situ experiments

The ability to predict the magnitude of an earthquake caused by deep fluid injections is an important factor for assessing the safety of the reservoir storage and the seismic hazard. Here, we propose a new approach to evaluate the seismic energy released during fluid injection by integrating injection parameters, induced aseismic deformation, and the distance of earthquake sources from injection. We use data from ten injection experiments performed at a decameter scale into fault zones in limestone and shale formations. We observe that the seismic energy and the hydraulic energy similarly depend on the injected fluid volume (V), as they both scale as V 3/2 . They show, however, a large discrepancy, partly related to a large aseismic deformation. Therefore, to accurately predict the released seismic energy, aseismic deformation should be considered in the budget through the residual deformation measured at the injection. Alternatively, the minimal hypocentral distance from injection points and the critical fluid pressure for fault reactivation can be used for a better prediction of the seismic moment in the total compilation of earthquakes observed during these experiments. Complementary to the prediction based only on the injected fluid volume, our approach opens the possibility of using alternative monitoring parameters to improve traffic-light protocols for induced earthquakes and the regulation of operational injection activities.

58 GEOSCIENCES↗
DOE OSTI2016

Inverting seismic data for rock physical properties; Mathematical background and application

The basic concept behind seismic inversion is that mathematical assumptions can be established to relate seismic to geological formation properties that caused their seismic responses. In this presentation we address some widely used seismic inversion method in hydrocarbon reservoirs identification and characterization. A successful use of the inversion in real example from gas sand reservoir in Boonsville field, Noth Central Texas is presented. Seismic data was not unambiguous indicator of reservoir facies distribution. The use of the inversion led to remove the ambiguity and reveal clear information about the target.

58 GEOSCIENCES↗
DOE OSTI2012

TIDAL DISSIPATION COMPARED TO SEISMIC DISSIPATION: IN SMALL BODIES, EARTHS, AND SUPER-EARTHS

While the seismic quality factor and phase lag are defined solely by the bulk properties of the mantle, their tidal counterparts are determined by both the bulk properties and the size effect (self-gravitation of a body as a whole). For a qualitative estimate, we model the body with a homogeneous sphere, and express the tidal phase lag through the lag in a sample of material. Although simplistic, our model is sufficient to understand that the lags are not identical. The difference emerges because self-gravitation pulls the tidal bulge down. At low frequencies, this reduces strain and the damping rate, making tidal damping less efficient in larger objects. At higher frequencies, competition between self-gravitation and rheology becomes more complex, though for sufficiently large super-Earths the same rule applies: the larger the planet, the weaker the tidal dissipation in it. Being negligible for small terrestrial planets and moons, the difference between the seismic and tidal lagging (and likewise between the seismic and tidal damping) becomes very considerable for large exoplanets (super-Earths). In those, it is much lower than what one might expect from using a seismic quality factor. The tidal damping rate deviates from the seismic damping rate, especially in the zero-frequency limit, and this difference takes place for bodies of any size. So the equal in magnitude but opposite in sign tidal torques, exerted on one another by the primary and the secondary, have their orbital averages going smoothly through zero as the secondary crosses the synchronous orbit. We describe the mantle rheology with the Andrade model, allowing it to lean toward the Maxwell model at the lowest frequencies. To implement this additional flexibility, we reformulate the Andrade model by endowing it with a free parameter {zeta} which is the ratio of the anelastic timescale to the viscoelastic Maxwell time of the mantle. Some uncertainty in this parameter's frequency dependence does not influence our principal conclusions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A viscoelastic model for seismic attenuation using fractal mechanical networks

SUMMARY Seismic attenuation (quantified by the quality factor Q) has a significant impact on the seismic waveforms, especially in the fluid-saturated rocks. This dissipative process can be phenomenologically represented by viscoelastic models. Previous seismological studies show that the Q value of Earth media exhibits a nearly frequency-independent behaviour (often referred to as constant-Q in literature) in the seismic frequency range. Such attenuation can be described by the mathematical Kjartansson constant-Q model, which lacks of a physical representation in the viscoelastic sense. Inspired by the fractal nature of the pore fluid distribution in patchy-saturated rocks, here we propose two fractal mechanical network (FMN) models, that is, a fractal tree model and a quasi-fractal ladder model, to phenomenologically represent the frequency-independent Q behaviour. As with the classic viscoelastic models, the FMN models are composed of mechanical elements (spring and dashpots) arranged in different hierarchical patterns. A particular parametrization of each model can produce the same complex modulus as in the Kjartansson model, which leads to the constant-Q. Applying the theory to several typical rock samples, we find that the seismic attenuation signature of these rocks can be accurately represented by either one of the FMN models. Besides, we demonstrate that the ladder model in particular exhibits the realistic multiscale fractal structure of the saturated rocks. Therefore, the FMN models as a proxy could provide a new way to estimate the microscopic rock structure property from macroscopic seismic attenuation observation.

Geochemistry & Geophysics↗
DOE OSTI2018

Sensitivity of Seismically Cued Antineutrino Detectors to Nuclear Explosions

Here, we evaluate the sensitivity of large gadolinium-doped water detectors to antineutrinos released by nuclear-fission explosions, using updated signal and background models and taking advantage of the capacity for seismic observations to provide an analysis trigger. Under certain realistic conditions, the antineutrino signature of a 250-kton pure-fission explosion could be identified several hundred kilometers away, in a detector about the size of the largest module currently proposed for a basic physics experiment. In principle, such an observation could provide rapid confirmation that the seismic signal coincided with a fission event, possibly useful for international monitoring of nuclear-weapon tests. We discuss the limited potential for seismically cued antineutrino observations to constrain fission yield, differentiate pure fission from fusion-enhanced weapon tests, indicate that the seismic evidence of an explosion had been intentionally masked, or verify the absence of explosive testing in a targeted area. We conclude that advances in seismic monitoring and neutrino physics have made the detection of explosion-derived antineutrinos more conceivable than previously asserted, but the size and cost of sufficiently sensitive detectors continue to limit applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Operational Forecasting of Induced Seismicity

Operational Forecasting of Induced Seismicity (ORION) is an open-source toolkit developed as part of the US Department of Energy SMART and NRAP programs. ORION is designed to help both expert and non-expert users understand the seismic hazard at a site. This tool analyzes how the seismic hazard evolves during injection of fluid, and can be used to identify potential mitigation strategies for observed or forecasted seismic activity. ORION is built using Python, and has an optional GUI where users can easily configure model parameters, evaluate temporal and spatial forecasts, and inspect key data (seismicity, fluid injection, pressure models, etc.).

Kroll, Kayla↗
DOE OSTI2025

Seismic H2: Version 1.0

Seismic-H2 is an integrated software package for geological hydrogen reservoir simulation, optimization, and leakage monitoring. The package includes multiple components: (1) code used for modeling seismic wave propagation in 3D heterogeneous elastic media based on finite-difference method to support detection of geological hydrogen storage reservoir leakage; (2) 3D reservoir simulations of leaks from an underground reservoir and 3D simulations of saline aquifers and depleted gas reservoirs; (3) seismic monitoring costs of passive and active seismic monitoring required for UHS; (4) rock physics calculations and interpolations for converting the reservoir simulations from part (2) into the elastic media models in part (1); (5) pre-processing seismic data; and lastly (6), a GUI interface that combines these different components.

Creasy, Neala↗
DOE OSTI2017

In situ measurement of velocity-stress sensitivity using crosswell continuous active-source seismic monitoring

The ability to characterize time-varying reservoir properties, such as the state of stress, has fundamental implications in subsurface engineering, relevant to geologic sequestration of CO 2 . Stress variation, here in the form of changes in pore fluid pressure, is one factor known to affect seismic velocity. Induced variations in velocity have been used in seismic studies to determine and monitor changes in the stress state. Previous studies conducted to determine velocity-stress sensitivity at reservoir conditions rely primarily on laboratory measurements of core samples or theoretical relationships. In this paper, we have developed a novel field-scale experiment designed to study the in situ relationship between pore-fluid pressure and seismic velocity using a crosswell continuous active-source seismic monitoring (CASSM) system. At the Cranfield, Mississippi, CO 2 sequestration field site, we actively monitored seismic response for five days with a temporal resolution of 5 min; the target was a 26 m thick injection zone at approximately 3.2 km depth in a fluvial sandstone formation (lower Tuscaloosa Formation). The variation of pore fluid pressure was obtained during discrete events of fluid withdrawal from one of the two wells and monitored with downhole pressure sensors. The results indicate a correlation between decreasing CASSM time delay (i.e., velocity change for a raypath in the reservoir) and periods of reduced fluid pore pressure. The correlation is interpreted as the velocity-stress sensitivity measured in the reservoir. This observation is consistent with published laboratory studies documenting a velocity ( V ) increase with an effective stress increase. A traveltime change ( d t ) of 0.036 ms is measured as the consequence of a change in pressure of approximately 2.55 MPa ( d P e ). For T = 13 ms total traveltime, the velocity-stress sensitivity is d V / V / d P e = d t / T / d P e = 10.9 × 10 - 4 / MPa . Finally, the overall results suggest that CASSM measurements represent a valid technique for in situ determination of velocity-stress sensitivity in field-scale monitoring studies.

58 GEOSCIENCES↗
DOE OSTI2016

Advances in Rotational Seismic Measurements

Rotational motion is increasingly understood to be a significant part of seismic wave motion. Rotations can be important in earthquake strong motion and in Induced Seismicity Monitoring. Rotational seismic data can also enable shear selectivity and improve wavefield sampling for vertical geophones in 3D surveys, among other applications. However, sensor technology has been a limiting factor to date. The US Department of Energy (DOE) and Applied Technology Associates (ATA) are funding a multi-year project that is now entering Phase 2 to develop and deploy a new generation of rotational sensors for validation of rotational seismic applications. Initial focus is on induced seismicity monitoring, particularly for Enhanced Geothermal Systems (EGS) with fracturing. The sensors employ Magnetohydrodynamic (MHD) principles with broadband response, improved noise floors, robustness, and repeatability. This paper presents a summary of Phase 1 results and Phase 2 status.

58 GEOSCIENCES↗
DOE OSTI2025

STILGAR Seismic Array Data

Two three-component dense seismic arrays were deployed above an active limestone mine in Pleasant Gap, Pennsylvania; one in the fall and one in the spring of 2023. Each campaign lasted approximately four weeks and included 80 to 100 seismic stations with stations spaced between 100 and 400 meters apart. Each campaign also included several small dense seismic arrays. Two types of seismic stations were deployed including FairfieldNodal ZLand 3C All-in-One seismometers and Geospace 3C GS-ONE geophone paired with GSX3 and GSX4 dataloggers. This dataset includes ~ 4 TB of three-component continuously recorded seismic data in the miniseed format data from both campaigns as well as a readme file for the data structure.

58 GEOSCIENCES↗
DOE OSTI2018

Surface Seismic Monitoring of Hydraulic Fracturing Test Site (HFTS) in the Midland Basin, Texas

Low-frequency seismic events of long duration were recorded by a surface array of 48, 3-component geophones during the fracture stimulation of two offset horizontal wells in the Middle Wolfcamp Formation (8000-8100 ft. below the geophone array) at the Hydraulic Fracturing Test Site (HFTS) in Reagan County, TX. Each recorded event was approximately 5-10 seconds long, with no clear first wave arrival; the seismic energy was largely visible on the vertical component suggesting that the events are dominated by P waves. A finite temporal moveout was observed across the seismic network in event arrival times indicating a causal relationship with the local source of deformation. The low-frequency seismic events observed in this study can be distinguished from the more commonly known microseismic events by their: 1) emergent waveform, 2) lack of clear P wave and S wave arrivals, 3) longer duration (few seconds), and 4) lower-frequency energy content (10-60 Hz). Previous researchers have observed similar low-frequency, long duration events in conjunction with the hydraulic fracturing of hydrocarbon reservoirs and have referred to them as long-period, long-duration (LPLD) events in the literature. LPLD events may be important to a comprehensive understanding of fracture stimulation because these events are thought to represent: 1) movement along pre-existing discontinuities such as bedding planes and fractures/faults, 2) deformation of ductile rocks with high clay content, and 3) the jerky tensile opening of hydraulic fractures. LPLD events resulting from slow shear-slip on faults and fractures of relatively larger dimension or rupture in the clay rich rock are collectively referred to as slow-slip seismicity, which was first observed in tectonic and volcanic settings, areas with high pore pressure like hydraulically fractured reservoirs. Because LPLDs are thought to represent deformation mechanisms not captured by microseismic analysis (which dominantly indicate shear failure of brittle rock), we believe that incorporating LPLD events into a new methodology for delimiting stimulated reservoir volume might improve predictions of well production.

Kumar, Abhash↗
DOE OSTI2018

Seismic Models for Near-Surface Explosion Yield Estimation in Alluvium and Sedimentary Rock

Seismic ground-motion data provide valuable constraints on explosion characteristics, such as yield and height-of-burst/depth-of-burial (HOB/DOB). In this study, we investigate a range of seismic amplitude features and their efficacy in minimizing errors associated with yield estimation. Using a set of explosion recordings from experiments conducted in alluvium and sedimentary rock geologies, we investigated the effectiveness of three different seismic feature types over a range of different frequencies. Using both velocity and displacement data, we investigated the zero-to-peak (ZTP) amplitude of the first arriving P wave, the peak-to-peak (PTP) amplitude of the first arriving P wave, and the P rms amplitude over various time windows starting with the P wave. These three basic features were measured on both vertical component only data and the vertical–radial vector sum. In total, there were 56 different combinations investigated. Our seismic models combine the effects of scaled range and scaled HOB/DOB on the observed seismic features. The results show that the vertical–radial vector sum of the ZTP amplitudes measured on displacement seismograms most consistently produce models with the smallest difference between predicted amplitude values and observed amplitude values when accounting for both alluvium and sedimentary rock lithology. The improvement in fit to the data when incorporating the difference in source lithology is significant. Lastly, if an alluvium lithology is assumed for a sedimentary rock lithology, the difference in yield could be up to ~3.0 times the appropriate yield in the far field.

58 GEOSCIENCES↗
DOE OSTI2020

Dense Seismic Array Study of a Legacy Underground Nuclear Test at the Nevada National Security Site

The complex postdetonation geologic structures that form after an underground nuclear explosion are difficult to constrain because increased heterogeneity around the damage zone affects seismic waves that propagate through the explosion site. Generally, a vertical rubble-filled structure known as a chimney is formed after an underground nuclear explosion that is composed of debris that falls into the subsurface cavity generated by the explosion. Compared with chimneys that collapse fully, leaving a surface crater, partially collapsed chimneys can have remnant subsurface cavities left in place above collapsed rubble. The 1964 nuclear test HADDOCK, conducted at the Nevada test site (now the Nevada National Security Site), formed a partially collapsed chimney with no surface crater. Understanding the subsurface structure of these features has significant national security applications, such as aiding the study of suspected underground nuclear explosions under a treaty verification. In this study, we investigated the subsurface architecture of the HADDOCK legacy nuclear test using hybrid 2D–3D active source seismic reflection and refraction data. The seismic data were acquired using 275 survey shots from the Seismic Hammer (a 13,000 kg weight drop) and 65 survey shots from a smaller accelerated weight drop, both recorded by ~1000 three-component 5 Hz geophones. First-arrival, P-wave tomographic modeling shows a low-velocity anomaly at ~200 m depth, likely an air-filled cavity caused by partial collapse of the rock column into the temporary postdetonation cavity. A high-velocity anomaly between 20 and 60 m depth represents spall-related compaction of the shallow alluvium. Hints of low velocities are also present near the burial depth (~364 m). The reflection seismic data show a prominent subhorizontal reflector at ~300 m depth, a short-curved reflector at ~200 m, and a high-amplitude reflector at ~50 m depth. Comparisons of the reflection sections to synthetic data and borehole stratigraphy suggest that these features correspond to the alluvium–tuff contact, the partial collapse cavity, and the spalled layer, respectively.

58 GEOSCIENCES↗
DOE OSTI2022

Small Seismic Events in Oklahoma Detected and Located by Machine Learning–Based Models

A complete earthquake catalog is essential to understand earthquake nucleation and fault stress. Following the Gutenberg–Richter law, smaller, unseen seismic events dominate the earthquake catalog and are invaluable for revealing the fault state. The published earthquake catalogs, however, typically miss a significant number of small earthquakes. Part of the reason is due to a limitation of conventional algorithms, which can hardly extract small signals from background noise in a reliable and efficient way. To address this challenge, we utilized a machine learning method and developed new models to detect and locate seismic events. These models are efficient in processing a large amount of seismic data and extracting small seismic events. We applied our method to seismic data in Oklahoma, United States, and detected ~14 times more earthquakes compared with the standard Oklahoma Geological Survey catalog. The rich information contained in the new catalog helps better understand the induced earthquakes in Oklahoma.

58 GEOSCIENCES↗
DOE OSTI2025

Bayesian Inference for the Seismic Moment Tensor Using Regional Waveforms and Teleseismic- P Polarities with a Data-Derived Distribution of Velocity Models and Source Locations

The largest source of uncertainty in any source inversion is the velocity model used in the transfer function that relates observed ground motion to the seismic moment tensor. However, standard inverse procedure often does not quantify uncertainty in the seismic moment tensor due to error in the Green’s functions from uncertain event location and Earth structure. Here, we incorporate this uncertainty into an estimation of the seismic moment tensor using a data-derived distribution of velocity models based on complementary geophysical data sets, including thickness constraints, velocity profiles, gravity data, surface-wave group velocities, and regional body-wave travel times. The data-derived distribution of velocity models is then used as a prior distribution of Green’s functions for use in Bayesian inference of an unknown seismic moment tensor using regional and teleseismic-P waveforms. The use of multiple data sets is important for gaining resolution to different components of the moment tensor. The combined likelihood is estimated using data-specific error models and the posterior of the seismic moment tensor is estimated and interpreted in terms of the most probable source type.

58 GEOSCIENCES↗
DOE OSTI2021

Toward Robust and Routine Determination of Mw for Small Earthquakes: Application to the 2020 Mw 5.7 Magna, Utah, Seismic Sequence

To better characterize seismic hazard, particularly, for induced seismicity, there is an increasing interest in methods to estimate moment magnitude (M w ) for small earthquakes. M w is generally preferred over other magnitude types, but, it is difficult to estimate M w for earthquakes with local magnitude (M L ) <3-3.5, using conventional moment tensor (MT) inversion. The 2020 M ww 5.7 Magna, Utah, seismic sequence provides an opportunity to illustrate and evaluate the value of spectral methods for this purpose. Starting with a high-quality seismic catalog of 2103 earthquakes (M L <5.6), we estimate M w using two independent spectral methods—one based on direct waves, yielding M w,direct , and the other based on coda waves, yielding M w,coda . For the direct-wave method, we present a non-parametric (NP) inversion scheme that solves for apparent geometrical spreading, G(R), and site effects (S), similar to other NP procedures that have been used to calibrate regional M L scales. The NP inversion is constrained using M ws derived from MTs for nine events in the Magna sequence. We recover statistically robust and physically reasonable G(R) and S and compute M w,direct for 635 Magna earthquakes down to M L 0.7. For the coda-wave method, we consider two separate calibration schemes involving previous MT solutions and compute M w,coda for 311 earthquakes down to M L 1.0. For 280 of the events that were processed with both methods—M w,direct and M w,coda —are strongly correlated (r = 0.98), with a mean difference of only 0.05. We compare M w,direct and M w,coda with M L and find reasonably good agreement for M L <3.6 with the theoretically predicted relationship of M w =(2/3)M L +C, in which C is a regional constant. Our results imply that seismic network operators can use spectral-based M w estimates to replace M L estimates for events with M L ≥1.0, and possibly smaller. The main requirement is the existence of a small number of MT solutions for calibration purposes.

58 GEOSCIENCES↗