SEARCH · Engineering Papers

Results for “Earthquakes”

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.

SourceResearch connectionFrom yearThrough year

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

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

NASA NTRS2002

Charge Generation and Propagation in Igneous Rocks

Various electrical phenomena have been reported prior to or concurrent with earthquakes such as resistivity changes, ground potentials, electromagnetic (EM), and luminous signals. Doubts have been raised as to whether some of these phenomena are real and indeed precursory. One of the reasons for uncertainty is that, despite decades of intense work, there is still no physically coherent model. Using low- to medium-velocity impacts to measure electrical signals with microsecond time resolution, it has now been observed that when dry gabbro and diorite cores are impacted at relatively low velocities, approximately 100 m/s, highly mobile charge carriers are generated in a small volume near the impact point. They spread through the rocks, causing electric potentials exceeding +400 mV, EM, and light emission. As the charge cloud spreads, the rock becomes momentarily conductive. When a dry granite block is impacted at higher velocity, approximately 1.5 km/s, the propagation of the P and S waves is registered through the transient piezoelectric response of quartz. After the sound waves have passed, the surface of the granite block becomes positively charged, suggesting the same charge carriers as observed during the low-velocity impact experiments, expanding from within the bulk. During the next 2-3 ms the surface potential oscillates, indicating pulses of electrons injected from ground and contact electrodes. The observations are consistent with positive holes, e.g., defect electrons in the O(2-) sublattice, traveling via the O 2p-dominated valence band of the silicate minerals. Before activation, the positive holes lay dormant in the form of electrically inactive positive hole pairs (PHP), chemically equivalent to peroxy links, O3X/OO\XO3, with X=Si(4+), Al(3+), etc. PHPs are introduced into the minerals by way of hydroxyl,O3X-OH, which all nominally anhydrous minerals incorporate when crystallizing in H2O-laden environments. The fact that positive holes can be activated by low-energy impacts, and their attendant sound waves, suggests that they can also be activated by microfracturing. Depending on where in the stressed rock volume the charge carriers are activated, they will form rapidly moving or fluctuating charge clouds that may account for earthquake-related electrical signals and EM emission. Wherever such charge clouds intersect the surface, high fields are expected, causing electric discharges and earthquake lights.

Freund, Friedemann↗

NASA NTRS2000

Charge Generation and Propagation in Igneous Rocks

Resistivity changes, ground potentials, electromagnetic (EM) and luminous signals prior to or during earthquakes have been reported, in addition to ground uplift and tilt, and to changes in the seismic wave propagation parameters. However, no physical model exists that ties these diverse phenomena together. Through time-resolved impacts experiments it has been observed that, when igneous rocks (gabbro, diorite, granite) are impacted at low velocities (approx. 100 m/sec), highly mobile electronic charge carriers are generated, spreading from a small volume near the impact point, causing electric potentials, EM and light emission. The rock becomes momentarily conductive. When impacted at higher velocities (approx. 1.5 km/sec), the propagation of the P and S waves is registered through the transient piezoelectric response of quartz. At the same time, the rock volume is filled with mobile charge carriers, and a positive surface potential is registered. During the next 1-2 msec the surface potential oscillates, due to electron injection from ground. These observations are consistent with positive holes, e.g. defect electrons in the O(2-) sublattice, that can travel via the O 2p-dominated valence band of the silicate minerals at the speed of a phonon-mediated charge transfer. Before activation, the positive hole charge carriers lay dormant in form of positive hole pairs, PHP, electrically inactive, chemically equivalent to peroxy links in the structures of constituent minerals. PHPs are introduced by way of hydroxyl (O3Si-OH) incorporated into nominally anhydrous minerals when they crystallize in water-laden environments. Given that sound waves of even relatively low intensity appear to cause PHPs dissociation, thus generating mobile positive holes, it is proposed that microfracturing during rock deformation cause PHP dissociation. Depending on where and how much the rock volume is stressed, the positive holes are expected to form fluctuating charge clouds in the earthquake source region that may account for earthquake-related electrical phenomena and the reported low frequency EM signals.

Freund, Friedemann↗

NASA NTRS2005

Glacial isostatic stress shadowing by the Antarctic ice sheet

Numerous examples of fault slip that offset late Quaternary glacial deposits and bedrock polish support the idea that the glacial loading cycle causes earthquakes in the upper crust. A semianalytical scheme is presented for quantifying glacial and postglacial lithospheric fault reactivation using contemporary rock fracture prediction methods. It extends previous studies by considering differential Mogi-von Mises stresses, in addition to those resulting from a Coulomb analysis. The approach utilizes gravitational viscoelastodynamic theory and explores the relationships between ice mass history and regional seismicity and faulting in a segment of East Antarctica containing the great Antarctic Plate (Balleny Island) earthquake of 25 March 1998 (Mw 8.1). Predictions of the failure stress fields within the seismogenic crust are generated for differing assumptions about background stress orientation, mantle viscosity, lithospheric thickness, and possible late Holocene deglaciation for the D91 Antarctic ice sheet history. Similar stress fracture fields are predicted by Mogi-von Mises and Coulomb theory, thus validating previous rebound Coulomb analysis. A thick lithosphere, of the order of 150-240 km, augments stress shadowing by a late melting (middle-late Holocene) coastal East Antarctic ice complex and could cause present-day earthquakes many hundreds of kilometers seaward of the former Last Glacial Maximum grounding line.

Coulomb stress↗

NASA NTRS2004

The QuakeSim Project: Numerical Simulations for Active Tectonic Processes

In order to develop a solid earth science framework for understanding and studying of active tectonic and earthquake processes, this task develops simulation and analysis tools to study the physics of earthquakes using state-of-the art modeling, data manipulation, and pattern recognition technologies. We develop clearly defined accessible data formats and code protocols as inputs to the simulations. these are adapted to high-performance computers because the solid earth system is extremely complex and nonlinear resulting in computationally intensive problems with millions of unknowns. With these tools it will be possible to construct the more complex models and simulations necessary to develop hazard assessment systems critical for reducing future losses from major earthquakes.

tectonic processes↗

NASA NTRS2011

High-Quality Seismic Observations of Sonic Booms

The SonicBREWS project (Sonic Boom Resistant Earthquake Warning Systems) is a collaborative effort between Seismic Warning Systems, Inc. and NASA Dryden Flight Research Center. This project aims to evaluate the effects of sonic booms on Earthquake Warning Systems in order to prevent such systems from experiencing false alarms due to sonic booms. The airspace above the Antelope Valley, California includes the High Altitude Supersonic Corridor and the Black Mountain Supersonic Corridor. These corridors are among the few places in the US where supersonic flight is permitted, and sonic booms are commonplace in the Antelope Valley. One result of this project is a rich dataset of high-quality accelerometer records of sonic booms which can shed light on the interaction between these atmospheric phenomena and the solid earth. Nearly 100 sonic booms were recorded with low-noise triaxial MEMS accelerometers recording 1000 samples per second. The sonic booms had peak overpressures ranging up to approximately 10 psf and were recorded in three flight series in 2010 and 2011. Each boom was recorded with up to four accelerometers in various array configurations up to 100 meter baseline lengths, both in the built environment and the free field. All sonic booms were also recorded by nearby microphones. We present the results of the project in terms of the potential for sonic-boom-induced false alarms in Earthquake Warning Systems, and highlight some of the interesting features of the dataset.

Wurman, Gilead↗

NASA NTRS2012

Plate Margin Deformation and Active Tectonics Along the Northern Edge of the Yakutat Terrane in the Saint Elias Orogen, Alaska and Yukon, Canada

The northwest directed motion of the Pacific plate is accompanied by migration and collision of the Yakutat terrane into the cusp of southern Alaska. The nature and magnitude of accretion and translation on upper crustal faults and folds is poorly constrained, however, due to pervasive glaciation. In this study we used high-resolution topography, geodetic imaging, seismic, and geologic data to advance understanding of the transition from strike-slip motion on the Fairweather fault to plate margin deformation on the Bagley fault, which cuts through the upper plate of the collisional suture above the subduction megathrust. The Fairweather fault terminates by oblique-extensional splay faulting within a structural syntaxis, allowing rapid tectonic upwelling of rocks driven by thrust faulting and crustal contraction. Plate motion is partly transferred from the Fairweather to the Bagley fault, which extends 125 km farther west as a dextral shear zone that is partly reactivated by reverse faulting. The Bagley fault dips steeply through the upper plate to intersect the subduction megathrust at depth, forming a narrow fault-bounded crustal sliver in the obliquely convergent plate margin. Since . 20 Ma the Bagley fault has accommodated more than 50 km of dextral shearing and several kilometers of reverse motion along its southern flank during terrane accretion. The fault is considered capable of generating earthquakes because it is linked to faults that generated large historic earthquakes, suitably oriented for reactivation in the contemporary stress field, and locally marked by seismicity. The fault may generate earthquakes of Mw <= 7.5.

Bruhn, Ronald L.↗

DOE OSTI2023

Moment tensor event identification for collapses

SUMMARY We introduce a seismic identification method for collapse events using moment tensors (MTs). We start by computing full (six-element) MT solutions for 43 identified collapse events from around the world, and statistically characterizing the population on the MT hypersphere. We then test a large data set of over 1000 full MTs for the western U.S. against the distribution of collapses using a MT-based identification method similarly as used for testing explosions. Known collapses and explosions are readily identified, along with other anomalous events in the Geysers and central California coast. Misidentification rates are determined for various screening angles with optimal misidentification rates between earthquakes and collapses on the order of 3 per cent. The method is demonstrated to be very effective at identifying non-earthquake sources with a 97–98 per cent accuracy. It is likely to be transportable to other regions, and can be used for event identification anywhere full MT solutions are routinely calculated.

58 GEOSCIENCES↗

DOE OSTI2024

Precise relative magnitude measurement improves fracture characterization during hydraulic fracturing

SUMMARY Microseismic monitoring is an important technique to obtain detailed knowledge of in-situ fracture size and orientation during stimulation to maximize fluid flow throughout the rock volume and optimize production. Furthermore, considering that the frequency of earthquake magnitudes empirically follows a power law (i.e. Gutenberg–Richter), the accuracy of microseismic event magnitude distributions is potentially crucial for seismic risk management. In this study, we analyse microseismicity observed during four hydraulic fracture treatments of the legacy Cotton Valley experiment in 1997 at the Carthage gas field of East Texas, where fractures were activated at the base of the sand-shale Upper Cotton Valley formation. We perform waveform cross-correlation to detect similar event clusters, measure relative amplitude from aligned waveform pairs with a principal component analysis, then measure precise relative magnitudes. The new magnitudes significantly reduce the deviations between magnitude differences and relative amplitudes of event pairs. This subsequently reduces the magnitude differences between clusters located at different depths. Reduction in magnitude differences between clusters suggests that some attenuation-related biases could be effectively mitigated with relative magnitude measurements. The maximum likelihood method is applied to understand the magnitude frequency distributions and quantify the seismogenic index of the clusters. Statistical analyses with new magnitudes suggest that fractures that are more favourably oriented for shear failure have lower b-value and higher seismogenic index, suggesting higher potential for relatively larger earthquakes, rather than fractures subparallel to maximum horizontal principal stress orientation.

58 GEOSCIENCES↗

DOE OSTI2022

High-Precision Characterization of Seismicity from the 2022 Hunga Tonga-Hunga Ha'apai Volcanic Eruption

The earthquake swarm accompanying the January 2022 Hunga Tonga-Hunga Ha'apai (HTHH) volcanic eruption includes a large number of posteruptive moderate-magnitude seismic events and presents a unique opportunity to use remote monitoring methods to characterize and compare seismic activity with other historical caldera-forming eruptions. We compute improved epicentroid locations, magnitudes, and regional moment tensors of seismic events from this earthquake swarm using regional to teleseismic surface-wave cross correlation and waveform modeling. Precise relative locations of 91 seismic events derived from 59,047 intermediate-period Rayleigh- and Love-wave cross-correlation measurements collapse into a small area surrounding the volcano and exhibit a southeastern time-dependent migration. Regional moment tensors and observed waveforms indicate that these events have a similar mechanism and exhibit a strong positive compensated linear vector dipole component. Precise relative magnitudes agree with regional moment tensor moment magnitude ($M_w$) estimates while also showing that event sizes and frequency increase during the days after the eruption followed by a period of several weeks of less frequent seismicity of a similar size. The combined information from visual observation and early geologic models indicate that the observed seismicity may be the result of a complex series of events that occurred after the explosive eruption on 15 January, possibly involving rapid resupply of the magma chamber shortly after the eruption and additional faulting and instability in the following weeks. In addition, we identify and characterize an $M_w$ 4.5 event five days before the paroxysmal explosion on 15 January, indicating that additional seismic events preceding the main eruption could have been identified with improved local monitoring. As a result, our analysis of the HTHH eruption sequence demonstrates the value of potentially utilizing teleseismic surface-wave cross correlation and waveform modeling methods to assist in the detailed analysis of remote volcanic eruption sequences.

58 GEOSCIENCES↗

DOE OSTI2022

A phase-field model for quasi-dynamic nucleation, growth, and propagation of rate-and-state faults

Despite its critical role in the study of earthquake processes, numerical simulation of the entire stages of fault rupture remains a formidable task. The main challenges in simulating a fault rupture process include the complex evolution of fault geometry, frictional contact, and off-fault damage over a wide range of spatial and temporal scales. Here, in this paper, we develop a phase-field model for quasi-dynamic fault nucleation, growth, and propagation, which features two standout advantages: (i) it does not require any sophisticated algorithms to represent fault geometry and its evolution; and (ii) it allows for modeling fault nucleation, propagation, and off-fault damage processes with a single formulation. Built on a recently developed phase-field framework for shear fractures with frictional contact, the proposed formulation incorporates rate- and state-dependent friction, radiation damping, and their impacts on fault mechanics and off-fault damage. We show that the numerical results of the phase-field model are consistent with those obtained from well-verified approaches that model the fault as a surface of discontinuity, without suffering from the mesh convergence issue in the existing continuous approaches to fault rupture (e.g., the stress glut method). Further, through numerical examples of fault propagation in various settings, we demonstrate that the phase-field approach may open new opportunities for investigating complex earthquake processes that have remained overly challenging for the existing numerical methods.

58 GEOSCIENCES↗

DOE OSTI2020

Geodetic Measurements of Slow-Slip Events Southeast of Parkfield, CA

Tremor and low-frequency earthquakes are presumed to be indicative of surrounding slow, aseismic slip that is often below geodetic detection thresholds. This study uses data from borehole seismometers and long-baseline laser strainmeters to observe both the seismic and geodetic signatures of episodic tremor and slip on the Parkfield region of the San Andreas Fault near Cholame, CA. We note the observed occurrence rates of both the tremors and co-located families of low-frequency earthquakes are not steady but instead exhibit quasiperiodic bursts of increased activity.We show that these periods of elevated seismic activity correlate with statistically significant stacked strain signals consisting of 44 slow-slip events. Modeled individual slow-slip events and their total summed moment, which are constrained by seismic signals and stacked strain, respectively, indicate that the individual moment magnitudes of these events range from $M_w$ 4.6–5.2.We find that the measured geodetic signal likely precedes the seismic signal by several hours, consistent with the aseismic slip preceding and driving the observed seismic tremor activity.We confirm that strike-slip faults, in addition to subduction zones, are capable of producing episodic tremor and slip.

58 GEOSCIENCES↗

DOE OSTI2025

Prediction of Complex Observed Shear Wave Splitting Patterns at Ryukyu Subduction Zone Using a Strong Intra-Slab Anisotropy Model

Complex shear wave splitting (SWS) patterns in subduction zones are often interpreted geodynamically as resulting from complex mantle flow; however, this may not always be necessary. We analyzed 7,093 high-quality SWS measurements from teleseismic S waves recorded by Hi-net stations across the Ryukyu arc in Japan. Our findings show a systematic rotation of the fast S polarization from trench-parallel to trench-perpendicular depending on the earthquake backazimuth. For the same earthquake, the measured splitting patterns also vary spatially across the southwest Japan. Using full-wave seismic modeling, we showed that a dipping slab with ~30% shear anisotropy of the tilted transverse isotropy (TTI) type, with a symmetry axis perpendicular to the slab interface, can predict the observed delay times and polarization rotation. Our results highlight the importance of considering dipping anisotropic slabs in interpreting SWS at subduction zones.

58 GEOSCIENCES↗

DOE OSTI2025

Precursory Locking Precedes Slip Events on Laboratory Fault

Earthquakes occur through repeated cycles of locking and slip. What is the relationship between transiently locked patches and the size, location, and timing of the subsequent slip events? Here we use an experimental fault made of transparent rubber to image the slip history of over 18,000 events. The soft elastic system allows for contained ruptures and thus a natural evolution of stress heterogeneity over multiple earthquake cycles. We observe a precursory locking phase preceding slip nucleation for events of all sizes. The integrated locking distance increases linearly with the moment of the eventual event. The most intense (largest slip deficit) locking events are spatially and temporally correlated with the hypocenters of the largest rupture events. These results show that the transient locking is statistically, but not deterministically, predictive of the location, timing, and size of the eventual slip event.

earthquakes↗

DOE OSTI2022

On various mechanisms of shear wave generation from underground chemical explosions in hard rocks

Significant shear waves observed during underground explosions can make them look similar to natural earthquakes. In the near field, they can be caused by the source geometry and emplacement conditions, heterogeneous and anisotropic properties of the rock mass, relaxation of a non-uniform in situ stress around the source after the departure of the shock wave as well as by activation of nearby faults. Here, the stress relaxation mechanism is important at depth when the difference in main principal stresses becomes significant. Another mechanism of shear motion generation is due to explosive gas flow from the source cavity into the cracks initiated by the shock wave and joints connected to the cavity. Also, the shock wave from the explosion can cause weakening of the rock either due to permanent damage or shock-induced liquefaction at saturated joints. Numerical modelling of underground chemical explosions in hard rock at various confinement stresses has been conducted to compare various mechanisms of shear wave generation and their seismic characteristics with that of a simple fault rupture process observed during an earthquake.

58 GEOSCIENCES↗

DOE OSTI2023

The effect of shear strain and shear localization on fault healing

The seismic cycle of repeated earthquake failure requires that faults regain frictional strength during the interseismic phase, when the fault is locked or undergoing quasi-static creep. Fault healing plays a central role in determining earthquake stress drop, recurrence interval, elastic radiation frequency and other source parameters. In particular, the longer a fault remains quasi-stationary, the stronger it becomes and the larger the potential stress drop can be for the next event. Here, we address the role of shear strain and strain localization on fault healing and healing rate. We performed slide-hold-slide friction experiments on quartz gouge in the double-direct shear configuration for shear strain up to 25 and hold times from 10 to 1000 s. The results show that both healing and healing rate increase nonlinearly with increasing shear strain. Frictional healing scales with volumetric strain within the laboratory fault zone. Using the volumetric strain upon reshear as a proxy for strain localization, we demonstrate that the capacity of a fault to heal is directly proportional to shear bandwidth and degree of strain localization. The more the deformation is localized, the higher are the healing and healing rate, and thus, the fault strength. Our data provide a framework for understanding variations in fault strength over the seismic cycle and the role of brecciation and strain localization on spatiotemporal variations in fault strength.

58 GEOSCIENCES↗

DOE OSTI2024

Volcanic arc rigidity variations illuminated by coseismic deformation of the 2011 Tohoku-oki M9

Rock strength has long been linked to lithospheric deformation and seismicity. However, independent constraints on the related elastic heterogeneity are missing, yet could provide key information for solid Earth dynamics. Using coseismic Global Navigation Satellite Systems (GNSS) data for the 2011 M9 Tohoku-oki earthquake in Japan, we apply an inverse method to infer elastic structure and fault slip simultaneously. We find compliant material beneath the volcanic arc and in the mantle wedge within the partial melt generation zone inferred to lie above ~100 km slab depth. We also identify low-rigidity material closer to the trench matching seismicity patterns, likely associated with accretionary wedge structure. Along with traditional seismic and electromagnetic methods, our approach opens up avenues for multiphysics inversions. Those have the potential to advance earthquake and volcano science, and in particular once expanded to InSAR type constraints, may lead to a better understanding of transient lithospheric deformation across scales.

58 GEOSCIENCES↗

DOE OSTI2024

Performance evaluation of the USGS velocity model for the San Francisco Bay Area

In this study, we evaluated the performance of the United States Geological Survey velocity model developed for the San Francisco Bay Area (SFBA), version 21.1. The evaluation was performed through high-resolution three-dimensional physics-based ground motion simulations of seven small-magnitude earthquakes (ranging from magnitude 3.8 to 4.4) that occurred on the eastern side of the San Francisco Bay. The simulations were performed in the frequency range from 0 to 5 Hz with a minimum shear-wave velocity of 250 m/s, which allowed the capture of wave propagation effects of the near-surface soft materials that characterize local basins. Based on the direct comparison of Fourier amplitude spectra between recorded and simulated ground motions for more than 250 stations, we found that the velocity model generally performs well in the frequency range of 0.2–5 Hz. The median value of the Fourier amplitude residuals was found to be near zero for all seven earthquakes. The slight over-prediction of 0.2 log-natural units at frequencies above 3 Hz in our simulations was attributed to the potentially inaccurate representation of the source radiation pattern by a double-couple point source model, and simple representation of shallow small-scale underground structural complexity in the velocity model. Maps of spectral amplitude differences between the simulated and recorded data were used to identify areas responsible for systematic ground motion over-predictions or under-predictions. For example, while some sub-domains over soft sediments show over-prediction patterns, the block east of the Hayward fault is prone to exhibit patterns of under-prediction. These maps can be used to guide future refinements of the SFBA velocity model. Since our simulation methodology allows for the decoupling of the source and wave propagation effects, the ground motion data generated by our simulations can also be used to quantify the epistemic uncertainty due to the velocity model, in empirically based ground motion estimates for the SFBA.

58 GEOSCIENCES↗

DOE OSTI2021

High-Quality Revision of the Israeli Seismic Bulletin

Seismic bulletins, with trustworthy phase picks, origin times, and source locations are key for regional seismic studies, such as travel-time (TT) tomography, attenuation tomography, and anisotropy studies. To lay the groundwork for such studies in Israel, we revised the seismic bulletin of Israel and the surrounding area and obtained a trustworthy TT data set. From the earthquake and explosion bulletins of the Geophysical Institute of Israel, we compiled a starting data set of about 123,000 earthquakes and explosions that occurred during the past 40 yr. After screening out the poorly recorded events, we were left with a data set of ~38,000 well-recorded events. We then revised the remaining data set in two consecutive steps. In the first, we reviewed and updated station metadata, including changes in station metadata parameters over time. In the second step, we jointly relocated a list of selected seismic events, using the Bayesian hierarchical location software package (BayesLoc) of Myers et al. (2007) that performs joint relocation of multiple events. We observed striking dissimilarities between the spatial distributions of the newly relocated catalog and the initial locations. Although the depth distribution of the starting catalog is trimodal with peaks at 0, 5, and 10 km, the distribution in this study is unimodal, with a broad peak between 7.5 and 12.5 km. By differencing the observed arrival times and the origin times obtained through relocation with BayesLoc, we obtained a revised TT database that consists of 261,336 Pg, 132,876 Pn, 114,816 Sg, and 60,394 Sn arrivals, from a set of 30,458 jointly relocated seismic sources. In this work, we compared prerevision and postrevision TTs as a function of epicentral distance and concluded that the revised data set contains far fewer outliers and inconsistencies than the original data set. The revised TT data set may be used for seismic studies, such as TT tomography, attenuation tomography, and anisotropy studies.

58 GEOSCIENCES↗