Search NASA⌕ Search

SEARCH · Search NASA

Results for “hypocenters”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Periodicity Analysis of Earthquake Occurrence and Hypocenter Depth Near Parkfield, California, 1994–2002 Versus 2006–2014

In central California, periodic earthquake occurrence suggests a relationship with annual hydrological, atmospheric, thermal, and tidal loadings. Here, in catalogs of declustered earthquakes within 100 km from Parkfield, CA, we study a semiannual periodicity for the monthly number of ≤7.2-km-deep earthquakes and the monthly median hypocenter depth over 1994–2002 and 2006–2014. Peak-trough months in fitted periodic components differ between time spans and event populations. Deeper earthquakes present no semiannual and a weak annual periodicity. Although modeled pore pressure shows a Spring peak, when added to the elastic Coulomb stress from surface hydrospheric loads, it fails to predict a 6-month periodicity for 2006–2014. In 1994–2002, the pore-pressure amplitude appears to be of same order as the elastic stress and may have had a stronger effect. In 2006–2014, we did not find load model parameters explaining the observed pattern or apparent changes following the 2003 San Simeon and 2004 Parkfield earthquakes.

58 GEOSCIENCES↗
DOE OSTI2021

Three-dimensional paganica fault morphology obtained from hypocenter clustering (L'Aquila 2009 seismic sequence, Central Italy)

In seismic modelling, fault planes are normally assumed to be flat due to the lack of data which can constrain fault morphology. However, incorporating 3D fault morphology is important for modelling several phenomena, for example calculating mainshock induced stress changes. Here we utilize a data-analytical method to unveil the 3D rupture morphology of faults using unsupervised clustering techniques applied to earthquake hypocenters in seismic sequences. We apply this method to the 2009 L'Aquila seismic sequence which involved a M W 6.1 mainshock on April 6th. We use a dataset of about 50,000 relocated events, mostly microearthquakes, reaching magnitude of completeness equal to 0.7. Clustering distinguishes the earthquakes as occurring in three main clusters along with other minor fault segments. We then represent the morphology of the main Paganica fault system (responsible for the largest mainshock) using splines. This method shows promise as a step toward robustly and quickly obtaining 3D rupture morphologies where earthquake sequences have been monitored. The 3D model is presented interactively online, and the processing is presented in an interactive Jupyter Notebook (https://bit.ly/2MnCFdj).

58 GEOSCIENCES↗
DOE OSTI2019

Optimal design of 3D borehole seismic arrays for microearthquake monitoring in anisotropic media during stimulations in the EGS collab project

Multiple U.S. national laboratories, universities and industrial collaborators are conducting collaborative research under the EGS Collab project supported by the U.S. Department of Energy, to understand the fracture creation and imaging during fracturing in enhanced geothermal systems. Microearthquake hypocenter locations and focal mechanisms are used to monitor hydraulic fracturing growth and characterization at the EGS Collab experimental site at the Sanford Underground Research Facility using seismic receivers in multiple monitoring wells. In this study, we develop a methodology for optimal design a 3D borehole seismic array for cost-effective seismic monitoring in anisotropic media using not only the relationship between receiver distributions and standard deviation errors of microearthquake hypocenter locations, but also that between receiver distributions and focal-mechanism inversion errors. Our results indicate that microearthquake hypocenter locations and focal mechanisms can be reasonably well reconstructed for the EGS Collab Experiment I using six monitoring wells, including four fracture-parallel monitoring wells and two orthogonal wells. Lastly, eight seismic receivers evenly distributed in four parallel monitoring wells or twelve receivers in all six monitoring wells are required for hypocenter location, and twelve receivers evenly distributed in six wells or sixteen receivers in four wells are needed for focal-mechanism inversion.

58 GEOSCIENCES↗
DOE OSTI2019

Effects of Fault Roughness on Coseismic Slip and Earthquake Locations

We report that fault zone structure is well known to exert strong controls on earthquake properties including coseismic slip distribution, rupture propagation direction, and hypocenter location. It has also been well established that the principal slip surface, which accommodates the majority of earthquake displacement, exhibits roughness at all scales following self-affine fractal distributions. Here we explore the relationship between fault roughness and specific earthquake properties including coseismic slip distribution and hypocenter location based on long-term simulations of earthquake catalogs on fractally rough faults. We begin by using the von Kármán autocorrelation function to procedurally generate single faults with different fractal roughness properties, which we place in a homogeneous elastic solid and apply pure right-lateral shear at a constant back slip rate with the earthquake simulator RSQSim. Running the simulations for 10,000 years each, we generate millions of earthquakes including thousands of events with Mw > 6.0, which rupture the surface. We show that the patterns of surface rupture in these large events follow self-affine fractal distributions with consistent fractal dimension related distinctly to the fractal dimension of the fault. In addition, the hypocenters of these large events occur in very specific predictable locations where the longest wavelength structure produces a stress asperity (i.e., restraining bend). The resulting patterns can explain many features observed on real fault systems, including clustered hypocenter locations, spatially variable coseismic slip distributions, and characteristic slip recurrence behavior. Lastly, these results demonstrate a quantitative link between a directly measurable fault property—roughness—and the properties of future earthquakes.

58 GEOSCIENCES↗
DOE OSTI2022

Unsupervised learning from three-component accelerometer data to monitor the spatiotemporal evolution of meso-scale hydraulic fractures

Enhanced geothermal systems can provide a substantial share of the global energy demand. There exist several hurdles in the engineering implementations of such geothermal systems. One such hurdle is the accurate monitoring of the fracture networks created in subsurface through hydraulic stimulation of these systems. Micro seismicity associated with the stimulation is the primary means to locate the event hypocenters for estimating the stimulated rock volume. Existing methods for location the hypocenters are restricted to only the highest amplitude impulsive signals that are simultaneously detected on several sensors. Consequently, a large portion (usually ~99%) of the measurements are left unused. In this paper, an unsupervised manifold-approximation followed by clustering of 3-component accelerometer data is used to analyze the seismicity recorded on a monitoring well. With this method, a larger portion of the measured signal is used for the monitoring of the hydraulic fracture network. We analyze the EGS Collab experiment 1 microseismic data, recorded at the Sanford Underground Research Facility, South Dakota. Using the data from a single three-component accelerometer, the polarization features viz. Azimuth, incidence, rectilinearity, and planarity are used as inputs for the unsupervised manifold approximation followed by clustering. Our study shows that density-based clusters in the projected 3D space correspond to distinct types of hydraulically fractured zones around the injection point. Finally, we show that the temporal evolution of these clusters can be used to track fracture creation and propagation.

58 GEOSCIENCES↗
DOE OSTI2018

Validation of Deterministic Broadband Ground Motion and Variability from Dynamic Rupture Simulations of Buried Thrust Earthquakes

According to the article, we numerically model broadband ground motion (up to 5–7.5 Hz⁠) from blind‐thrust scenario earthquakes matching the fault geometry of the 1994 Mw 6.7 Northridge earthquake. Several realizations are modeled (by varying the hypocenter location in the dynamic rupture simulation) in a 1D‐layered velocity profile. In addition, we include Q(f )⁠, nonlinear effects from Drucker–Prager plasticity, and superimpose small‐scale medium complexity in both a 1D‐layered and 3D velocity model within the subsequent wave propagation. We investigate characteristics of the ground motion and its variability up to 50 km from the fault by comparing them with ground‐motion prediction equations (GMPEs), simple proxy metrics, as well as strong ground motion records from the Northridge event. We find that median ground motion closely follows the trend predicted by GMPEs and that the intraevent standard deviation, although varying with hypocenter location, lies near that of GMPE models. Plasticity affects ground‐motion amplitudes in regions near the source, reducing intraevent variability above ∼0.5 Hz⁠. Heterogeneity in the velocity structure on both the regional and small scales is needed for the simulated data to match two proxy metrics: the period‐to‐period correlation of spectral acceleration (SA) and the ratio of maximum‐to‐median SA. Although small‐scale heterogeneity has a negligible effect on median SA for this style of rupture, it serves to significantly increase the cumulative absolute velocity, better agreeing with observations. When compared with strong‐motion data, we find that long‐wavelength velocity structure within our deterministic simulations reduces bias at both short and long periods. Finally, synthetic ground motion at both footwall and hanging‐wall sites has no clear dependence on the distance to rupture (at both short and long periods); directivity is likely overpowering any hanging‐wall effect.

Geochemistry & Geophysics↗
DOE OSTI2021

Complex Fault Geometry of the 2020 M ww 6.5 Monte Cristo Range, Nevada, Earthquake Sequence

On 15 May 2020 an M ww 6.5 earthquake occurred beneath the Monte Cristo Range in the Mina Deflection region of western Nevada. Rapid deployment of eight temporary seismic stations enabled detailed analysis of its productive and slowly decaying aftershock sequence (p=0.8), which included ~18,000 autodetected events in 3.5 months. Double-difference, waveform-based relative relocation of 16,714 earthquakes reveals a complex network of faults, many of which cross the inferred 35-km-long east–northeast-striking, left-lateral mainshock rupture. Seismicity aligns with left-lateral, right-lateral, and normal mechanism moment tensors of 128 of the largest earthquakes. The mainshock occurred near the middle of the aftershock zone at the intersection of two distinct zones of seismicity. In the western section, numerous subparallel, shallow, north-northeast-striking faults form a broad flower-structure-like fault mesh that coalesces at depth into a near-vertical, left-lateral fault. We infer the near-vertical fault to be a region of significant slip in the mainshock and an eastward extension of the left-lateral Candelaria fault. Near the mainshock hypocenter, seismicity occurs on a northeast-striking, west-dipping structure that extends north from the eastern Columbus Salt Marsh normal fault. Together, these two intersecting structures bound the Columbus Salt Marsh tectonic basin. East of this intersection and the mainshock hypocenter, seismicity occurs in a narrow, near-vertical, east-northeast-striking fault zone through to its eastern terminus. At the eastern end, the aftershock zone broadens and extends northwest toward the southern extension of the northwest-striking, right-lateral Petrified Springs fault system. The eastern section hosts significantly fewer aftershocks than the western section, but has more moment release. We infer that shallow aftershocks throughout the system highlight fault-fracture meshes that connect mapped fault systems at depth. Here, comparing earthquake data with surface ruptures and a simple geodetic fault model sheds light on the complexity of this recent M 6.5 Walker Lane earthquake.

58 GEOSCIENCES↗
DOE OSTI2023

How Good Is Your Location? Comparing and Understanding the Uncertainties in Location for the 1993 Rock Valley Sequence

Accurate event locations are important for many endeavors in seismology, and understanding the factors that contribute to uncertainties in those locations is complex. In this article, we present a case study that takes an in-depth look at the accuracy and precision possible for locating nine shallow earthquakes in the Rock Valley fault zone in southern Nevada. These events are targeted by the Rock Valley Direct Comparison phase of the Source Physics Experiment, as candidates for the colocation of a chemical explosion with an earthquake hypocenter to directly compare earthquake and explosion sources. For this comparison, it is necessary to determine earthquake hypocenters as accurately as possible so that different source types have nearly identical locations. Our investigations include uncertainty analysis from different sets of phase arrivals, stations, velocity models, and location algorithms. For a common set of phase arrivals and stations, we find that epicentral locations from different combinations of velocity models and algorithms are within 600 m of one another in most cases. Event depths exhibit greater uncertainties, but focusing on the S-P times at the nearest station allows for estimates within approximately 500 m.

58 GEOSCIENCES↗
DOE OSTI2019

Earthquake Locations

The location of an earthquake is a classic inverse problem. Arrival times of various seismic phases are combined with the location of the instrument to form the data vector. The inverse operator is applied in an iterative manor until the solution converges and the resultant solution is the hypocentral location of the earthquake in space and time. In this presentation, I will present how seismologists obtain local earthquake hypocenters. I will present how the generalized inverse operator is built using the simple distance = rate x time equation and how the inverse operator is used in an iterative manor to converge on the hypocenter of the earthquake. Once an earthquake location is determined, seismologists can then use the location to infer other information about the medium through which the waves traveled. One possible application of the earthquake locations is a branch of seismology known as tomography. In tomography, source and station locations are used to compute the travel times of seismic phases and these phases are then inverted to obtain a velocity model of the medium. Seismic tomography can be applied in both 2-D and 3-D settings and in this presentation, I will present examples from both. The first example comes from a technique known as ambient noise tomography. Rayleigh waves are obtained by cross correlating the background noise between two stations. Because we know the interstation distance between the two stations and the arrival times of the Rayleigh waves, we can use tomography to build a group velocity map of the study region. I will present results from a study of the Eastern Tennessee Seismic Zone using the EarthScope’s Transportable Array together with different seismic networks in the region, figure 1. The second example uses P and S-waves from 412 earthquakes which occurred during the 1994 Draney Peak swarm to calculate a 3-D velocity model of the area. Preliminary results show a strong velocity contrast at a depth of about 4 km where the majority of the earthquakes occur, figure 2.

58 GEOSCIENCES↗
DOE OSTI2018

Seismic Analysis of Spatio-Temporal Fracture Generation During EGS Resource Development - Deviatoric MT, Fracture Network, and Final Report

This submission contains 167 deviatoric moment tensor (MT) solutions for the seismicity observed two years prior and three years post start of injection activities at The Geysers Prati 32 EGS Demonstration. Also included is a statistical representation of the properties of 751 fractures derived from the analysis of seismicity observed two years prior and three years post start of injection activities at The Geysers Prati 32 EGS Demonstration Project. The locations of the fractures are taken from microseismic hypocenters, the fracture areas are derived from moment magnitudes via scaling relationships, and the azimuths (sigma 1) and dips (sigma 3) are derived from the results of stress analyses.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Simulation-based characterization of the variability of earthquake risk to buildings in the near-field

Recent advancements in high performance computing platforms and computational workflow for regional-scale simulations are enabling unprecedented modeling of fault-to-structure earthquake processes. Regional simulations resolving ground motions at frequencies relevant to engineered systems are becoming computationally viable and provide a new capability to improve understanding of the geographical distribution and intensity of risk to buildings and critical infrastructure. As computational capabilities advance, it is essential to move beyond illustrative single rupture realizations for scenario earthquake events towards the development of a full suite of rupture realizations that appropriately characterize the range of risk to building systems. The work described in this article investigates the application of a suite of fault rupture realizations with the objective of assessing near-fault, site-specific seismic demand variability for building structures. A representative high-performance regional-scale computational model is utilized to execute ground motion and building response simulations based on 18 kinematic rupture realizations of an M7 strike-slip scenario earthquake. The fault rupture models for the scenario earthquake are created by systematically perturbing the hypocenter location and stochastically generating rupture parameters (slip, rise time, rake angle) to represent a breadth of ground motion intensities resulting from the spatial and temporal variabilities of an earthquake rupture process. The resulting seismic demand variability for three-story (short period) and forty-story (long period) steel moment-resisting frame buildings is characterized in terms of the median and distribution of peak inter-story drift ratio for a range of near-fault sites. The full suite of 18 fault rupture realizations and approximately 280,000 nonlinear dynamic building simulations indicate that the three-story building undergoes higher median seismic demand and significantly greater variability of demand at a given site than the forty-story building, which has important implications for the level of certainty in predicting building performance during an earthquake. The simulations performed provide deeper insight into the relationship between fault rupture parameterization and building response, which is essential information for developing a representative suite of rupture realizations for specific earthquake scenarios.

58 GEOSCIENCES↗
DOE OSTI2022

Sequences of seismic and aseismic slip on bimaterial faults show dominant rupture asymmetry and potential for elevated seismic hazard

Here we perform numerical simulations of sequences of earthquake and aseismic slip on planar rate and state faults separating dissimilar material within the 2-D plane strain approximation. We resolve all stages of the earthquake cycle from aseismic slip to fast ruptures while incorporating full inertia effects during seismic event propagation. We show that bimaterial coupling results in favorable nucleation site and subsequent asymmetric rupture propagation. We demonstrate that increasing the material contrast enhances this asymmetry leading to higher slip rates and normal stress drops in the preferred rupture propagation direction. The normal stress drop, induced by the bimaterial effect, leads to strong dynamic weakening of the fault and may destabilize the creeping region on a heterogeneous rate and state fault, resulting in extended rupture propagation. Such rupture penetration into creeping patches may lead to more frequent opening of earthquake gates, causing increased seismic hazard. Furthermore, bimaterial coupling may lead to irregular seismicity pattern in terms of event length, peak slip rates, and hypocenter location, depending on the properties of the creeping patches bordering the seismogenically active part of the fault. Our results highlight robust characteristics of bimaterial interfaces that persist over long sequence of events and suggest the need for further exploration of the role of material contrast in earthquake physics and models of seismic hazard.

58 GEOSCIENCES↗
DOE OSTI2018

Microseismicity cloud can be substantially larger than the associated stimulated fracture volume: the case of the Paralana Enhanced Geothermal System

Here, the goal of hydraulic stimulation is to increase formation permeability in the near vicinity of a well. However, there remain technical challenges around measuring the outcome of this operation. During two Enhanced Geothermal System (EGS) stimulations in South Australia, Paralana in 2011 and Habanero in 2003, extensive catalogs of microseismicity were recovered. It is often assumed that shear failure of existing fractures is the main mechanism behind both the induced earthquakes and any permeability enhancement. This underpins a common notion, that the seismically active volume is also the stimulated reservoir. In this paper, we compute the density of earthquake hypocenters and provide evidence that, under certain conditions, this spatiotemporal quantity is a reasonable proxy for pore pressure increase. We then apply an inverse modeling approach that uses the earthquake observations and a modified reservoir simulator to estimate the parameters of a permeability evolution relation. The regime implied by the data indicates that most permeability enhancement occurred very near to the wellbore and was not coincident with the bulk of the seismicity, whose volume was about two orders of magnitude larger. Thus, we conclude that, in some cases, it is possible for permeability enhancement and induced seismicity to be decoupled, in which case the seismically active volume is a poor indicator of the stimulated reservoir. Our results raise serious questions about the effectiveness of hydroshearing as a stimulation mechanism in EGS. Finally, this study extends our understanding of the complex processes linking earthquakes, fluid pressure, and permeability in both natural and engineered settings.

58 GEOSCIENCES↗
DOE OSTI2020

Machine-Learning-Based Analysis of the Guy-Greenbrier, Arkansas Earthquakes: A Tale of Two Sequences

Here, we revisited the June 2010 to October 2011 Guy-Greenbrier earthquake sequence in central Arkansas using PhaseNet, a deep neural network trained to pick P and S arrival times. We applied PhaseNet to continuous waveform data and used phase association and hypocenter relocation to locate nearly 90,000 events. Our catalog suggests that the sequence consists of two adjacent earthquake sequences on the same fault and that the second sequence may be associated with the wastewater disposal well to the west of the Guy-Greenbrier Fault, rather than the wells to the north and the east that were previously implicated. We find that each sequence is composed of many small clusters that exhibit diffusion along the fault at shorter timescales. Our study demonstrates that machine-learning-based earthquake catalog development is now feasible and will yield new insights into earthquake behavior.

58 GEOSCIENCES↗
DOE OSTI2020

Fracture Caging to Limit Induced Seismicity

Geothermal resources offer a stable low-carbon energy source. However, geothermal sites can collocate with the hypocenters of large-magnitude seismic events. Large seismic events pose a risk to public safety and are therefore a liability for efforts to develop geothermal resources. In this work, we propose “fracture caging” to limit induced seismic event magnitudes and present evidence from numerical model predictions, laboratory experiments, and field observations. Fracture caging involves drilling tactical production wells around a geothermal injection zone to contain fluids in fracture-dominated flow systems. Prior to our work, the effect of small wells on the growth of large fractures and on flow through fractures was subject to debate. Our work shows that production wells can impede fracture growth and contain high-pressure fluids in fracture-dominated rocks. This containment offers a mechanism to limit induced seismicity.

58 GEOSCIENCES↗
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↗