Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic networks”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Seismic source characterization of the Arabian Peninsula and Zagros Mountains from regional moment tensor and coda envelopes

Reliable estimates of moment magnitude and source mechanism for seismic events in the Middle East can be challenging due to a small number of openly available stations, the complex tectonic setting, and regions of high attenuation. Access to high-quality waveform data from well-calibrated regional seismic stations is fundamental in producing robust and stable estimates of earthquake source parameters, particularly when measurements of absolute ground motion amplitudes are required. Earthquake source mechanism and moment magnitude are invaluable information in the assessment of seismic hazard, plate motions, and the characterization of faults and regional stress field. The expansion of the Kingdom of Saudi Arabia’s national seismic network in the past several decades provides an opportunity to develop the capabilities of routine focal mechanism and moment magnitude estimations. In this study, we use time-domain full waveform moment tensor inversion and coda envelope–derived amplitude measurements to solve for earthquake source mechanism, moment magnitudes, and their source type. Furthermore, we compared the moment magnitudes calculated from the two methods and publicly available earthquake catalogs and discuss the implications of the obtained source parameters.

58 GEOSCIENCES↗
DOE OSTI2021

Rayleigh Wave Propagation in the Bighorn Mountains Region, Wyoming

Short-period Rayleigh waves, Rg , provide strong constraints on the depth of shallow seismic events and are of interest for monitoring small explosions. Characterizing the seismic sources that generate Rg requires an understanding of how shallow crustal structure affects Rayleigh wave propagation. Here, in support of these efforts, this study utilizes observed waveforms from small shallow explosions recorded on temporary seismic network deployments in the Bighorn region, Wyoming. We study regional near-surface geology by measuring changes in surface-wave amplitude and polarization during propagation through basins, foothills, and mountains. We develop additional insight by carrying out surface-wave eigenfunction analyses and numerical-wave simulations, which together reproduce many characteristics seen in the observed waveforms. Our results show how sedimentary basins in the Bighorn region allow for amplified prograde-polarized higher-mode and retrograde-polarized fundamental-mode Rayleigh waves, whereas adjacent mountains only support retrograde motion. These different modes provide distinct constraints on the Earth structure and source characteristics, potentially enabling targeted inversions in future studies. Finally, our findings provide insight into Rg propagation through complex near-surface geology, improving our understanding of shallow propagation and source effects that are relevant to explosion monitoring efforts.

58 GEOSCIENCES↗
DOE OSTI2024

Inferring the Focal Depths of Small Earthquakes in Southern California Using Physics-Based Waveform Features

Determining the depths of small crustal earthquakes is challenging in many regions of the world, because most seismic networks are too sparse to resolve trade-offs between depth and origin time with conventional arrival-time methods. Precise and accurate depth estimation is important, because it can help seismologists discriminate between earthquakes and explosions, which is relevant to monitoring nuclear test ban treaties and producing earthquake catalogs that are uncontaminated by mining blasts. Here, we examine the depth sensitivity of several physics-based waveform features for ~8000 earthquakes in southern California that have well-resolved depths from arrival-time inversion. We focus on small earthquakes (2 < M L < 4) recorded at local distances (<150 km), for which depth estimation is especially challenging. We find that differential magnitudes (M w /M L –M c ) are positively correlated with focal depth, implying that coda wave excitation decreases with focal depth. We analyze a simple proxy for relative frequency content, Φ≡log 10 (M 0 )+3log 10 (f c ), and find that source spectra are preferentially enriched in high frequencies, or “blue-shifted,” as focal depth increases. Here, we also find that two spectral amplitude ratios Rg 0.5–2 Hz/Sg 0.5–8 Hz and Pg/Sg at 3–8 Hz decrease as focal depth increases. Using multilinear regression with these features as predictor variables, we develop models that can explain 11%–59% of the variance in depths within 10 subregions and 25% of the depth variance across southern California as a whole. We suggest that incorporating these features into a machine learning workflow could help resolve focal depths in regions that are poorly instrumented and lack large databases of well-located events. Some of the waveform features we evaluate in this study have previously been used as source discriminants, and our results imply that their effectiveness in discrimination is partially because explosions generally occur at shallower depths than earthquakes.

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↗
DOE OSTI2022

A Project Lifetime Approach to the Management of Induced Seismicity Risk at Geologic Carbon Storage Sites

The geologic storage of carbon dioxide (CO 2 ) is one method that can help reduce atmospheric CO 2 by sequestering it into the subsurface. Large-scale deployment of geologic carbon storage, however, may be accompanied by induced seismicity. We present a project lifetime approach to address the induced seismicity risk at these geologic storage sites. This approach encompasses both technical and nontechnical stakeholder issues related to induced seismicity and spans the time period from the initial consideration phase to postclosure. These recommendations are envisioned to serve as general guidelines, setting expectations for operators, regulators, and the public. They contain a set of seven actionable focus areas, the purpose of which are to deal proactively with induced seismicity issues. Although each geologic carbon storage site will be unique and will require a custom approach, these general best practice recommendations can be used as a starting point to any site-specific plan for how to systematically evaluate, communicate about, and mitigate induced seismicity at a particular reservoir.

58 GEOSCIENCES↗
DOE OSTI2025

Cardinal: Seismic and Geoacoustic Array Processing

Data collected via seismic and infrasound array deployments are leveraged in the geosciences to detect and characterize a myriad of natural and anthropogenic sources. These deployments consist of numerous sensors placed in a predetermined configuration to amplify signal strength and improve the efficacy of array processing techniques used to measure signal directionality and waveform coherence. High‐fidelity feature extraction is often predicated on interstation distance as well as the frequency content and wavelength of an incident signal. Numerous array processing softwares analyze data in sequential frequency bands to obtain a more detailed characterization of a signal. However, current algorithms are limited in their ability to determine optimal array configuration for each band. We introduce an open‐source Python code, called Cardinal, to process seismic and infrasound array data in discretized time–frequency space with the option of applying an adaptive array design to determine optimal subarray configuration for each frequency band. To reduce computational time, the array processing step can be run in parallel using multithreading. Furthermore, the software has the capability to aggregate array processing results from different time–frequency pixels to produce separate sets of detections, or families, with added utility via the application of an adaptive semblance threshold, which aids in isolating signals‐of‐interest from coherent background noise. Upon appropriate configuration, Cardinal exhibits the potential to combine distinct seismic and infrasound phases into separate families.

Adaptive Array↗
DOE OSTI2025

PySolate : A Python‐Based Thresholding Tool to Denoise or Designal Seismic Waveforms Based on the Continuous Wavelet Transform

PySolate is a Python‐based toolset that implements the continuous wavelet transform and nonlinear thresholding operations to denoise or designal seismic data, following Langston and Mousavi (2019). This filtering approach can remove microseismic noise to isolate intermediate‐period seismic signals that are key to enabling full‐waveform modeling and analysis of smaller‐magnitude regional events. This approach is best for the application to signals with frequency or time separation of signal and noise, in contrast to Fourier analysis, which is effective when signal and noise are separated in frequency. We demonstrate the Python toolset using the six announced Democratic People’s Republic of Korea declared nuclear tests, showing the effectiveness of isolating the seismic signal compared to standard bandpass filtering. In conclusion, we also demonstrate the ease of using the toolset with any Python processing tools.

Asia↗
DOE OSTI2020

Report on the recovery of legacy 1997-1998 Depth-of-Burial / Tunnel-and-Hole-Closure explosion data from the Kazakhstan former Soviet nuclear test site, Semipalatinsk.

In 1996, negotiations between the National Nuclear Center of the Republic of Kazakhstan (NNC), the U.S. Defense Special Weapons Administration (DSWA) and the U.S. Department of Energy led to an agreement for a joint series of experiments at the former Soviet nuclear test site Semipalatinsk (STS). These experiments had a dual purpose: to close holes and tunnels originally created for the purpose of testing nuclear weapons, and to characterize explosive sources at a former nuclear test site. This report covers the recovery of legacy data from the experiments conduct in 1997 and 1998. Two teams collected data. Los Alamos National Laboratory (LANL), in collaboration with the NNC, collected the near-source (local) data out to approximately 20 km. Lamont-Doherty Earth Observatory (LDEO), with funding from Lawrence Livermore Laboratory (LLNL), also in collaboration with NNC, collected regional data, 200 – 1000 km distant, leveraging the KZ regional seismic network that LDEO had recently installed in collaboration with NNC. The two teams also worked together to some extent, regarding instrumentation and planning. This report primarily addresses the near-source data collected by the LANL/NNC team.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2020

Aftershocks of a Chemical Explosion in Granite from the Source Physics Experiment Phase I

The Source Physics Experiment (SPE) is a multiphase experiment to better understand explosion source physics and thereby improve explosion monitoring. The rst phase (Phase I) examined the explosion source in a hard-rock medium and took place in a granite outcrop called the Climax Stock section of the Nevada National Security Site (NNSS). This series of chemical explosions are called the Wet Granite Geology (WGG) events, but are usually referred to simply as the SPE events (i.e., SPE-1, SPE-2, etc.). Event information for the chemical explosions is given in Table 1. The study of explosion aftershocks can aid in understanding of in-situ stresses and be used as a discriminant for earthquakes [Ford and Walter, 2010]. Sweeney and Harben reported no observed aftershocks for SPE-1 and -2, respectively. They calculated a minimum magnitude of observation based on the noise oor of the seismic array and predicted at least one or two events in the week after the events using the hard-rock aftershock model from Ford and Labak. Possible explanations for the null observation ranged from instrumental - poor seismometer coupling, to physical - the shallow depth of burial was unable to access the deep tectonic stress that causes some explosion aftershocks. We extend the analysis of Sweeney and Harben to the largest of the SPE Phase I chemical explosions, SPE-5, which was approximately a 5 tonne TNT-equivalent shot at 76.5 m depth. This event registered as an M L 2 on the University of Nevada, Reno seismic network and should produce a measurable aftershock sequence.

58 GEOSCIENCES↗
DOE OSTI2021

The Modeled Performance of the ORNL Matched Filter Detector for Seismic Signatures of Multi-Axle Tractor Trailers

The MINOS program seeks methods to passively monitor the transfer of special materials, which are reproducible and quantitatively defensible. To date, we have demonstrated that short-period seismometers deployed near vehicle tracks can detect and monitor signatures of multi-axle tractor trailers impacting small road obstacles as means to detect and identify some transfer activities. In this report, we model the seismic velocity signal output from a tractor-trailer impacting a small obstacle (a template waveform) to parameterize a matched filter that identifies signals with similar shapes buried in noise, and that indicate repeating events (target waveforms). We then use this model to perform thousands of synthetic experiments, in which we apply our detector to identify noisy waveforms output from a similarly modeled vehicle data set; we process all data in the same manner as we process real data collected from the Oak Ridge National Laboratory (ORNL) seismic network. Our experiments, while limited, reveal that the matched filter shows a high discrimination capability, and at very low SNR values (at visibility thresholds). This means that template waveforms can match target waveforms that are sourced by models that share the same vehicle speed (target vehicles), and have a low probability of mis-identification. In particular, a matched filter designed to detect a five-axle tractor-trailer will reliably identify that target vehicle down to an SNR of 11.25 dB, and within ±0.1 km/hr, while successfully screening target vehicles with three axles.

42 ENGINEERING↗
DOE OSTI2021

Central Asia Seismic Hazard Assessment (Casha) "Complete" Catalog of Earthquakes

Probabilistic Seismic Hazard Assessments (PSHA) greatly benefit from reliable earthquake catalogs. Central Asia, while rich in seismic monitoring history, reports a fraction of their earthquake data to international data collection centres. Furthermore, since the collapse of Soviet Union until now, monitoring and data collection efforts in Central Asia have not been centralized, making monitoring and reporting uneven and disconnected in this time period. Lawrence Livermore National Laboratory (LLNL) has undertaken a multi-year effort to systematically work with national seismic monitoring and data centres and seismic networks in Kazakhstan, Kyrgyzstan, and Tajikistan to compile a unified seismological bulletin for Central Asia. In tandem with this effort, an earthquake catalog that is processed to be PSHA ready was also generated. The processing for this effort included further duplicate identification and removal, magnitude harmonization and catalog completeness analyses. This report describes the “complete” PSHA-ready earthquake catalog, including the completeness intervals that were applied.

58 GEOSCIENCES↗
DOE OSTI2022

Analysis of the Source Physics Experiment Dry Alluvium Geology Ball Drop Sources

Two balls were dropped approximately 150 m down an open borehole as part of an investigation into the aftershocks from a 51 ton TNT-equivalent chemical explosion. The ball drop sources, a concrete-filled basketball and a bowling ball, were intended to mimic falling debris within an underground explosion borehole and cavity. Falling debris or sloughing of the borehole and cavity walls is thought to be one of the sources of explosion aftershocks. The signals from the ball drop sources were recorded by the same seismic network that monitored the aftershocks of the chemical explosion, so the ball drop recordings were used as templates to search through the continuous stream of aftershocks and detect correlated signals as a proportion of all aftershocks. The search resulted in only four similar event recordings out of the thousands of detected aftershocks. A synthetic test showed that the small, high-frequency ball drop templates were limited to producing highly correlated detections for events occurring in a small region around the impact location at 150 m depth. The high-frequency templates were necessary since the small ball drop signals were only above the noise at frequencies greater than 10 Hz. Unfortunately, the impact location is where the borehole was filled with sand in preparation for the subsequent chemical explosion in the test series, making it 150 m above the explosion created cavity centered at 300 m depth. This puts the ball drop templates 150 m away from the cavity, where many more events from falling debris would be expected. Further synthetic tests confirmed events emanating from the cavity 150 m below the impact site would have correlated with the lower frequency template and been detected had the template been usable below 10 Hz. Therefore, we conclude that the four detected events were most probably falling debris within the open borehole that impacted the floor at 150 m depth just as the ball drop sources had, and that this type of explosion aftershock mechanism comprises less than 1% of all explosion aftershocks.

58 GEOSCIENCES↗
DOE OSTI2024

Collaborative Research: Properties and Dynamics of the Shallow Crust (Final Report)

Ground motions recorded at one location are often extrapolated to nearby regions within a given radius. Here, we explore the appropriateness of spatial extrapolation using data from seven small aperture seismic network deployments in southern California. Six of these deployments are linear arrays of 4-13 stations, and one is a 2D array of 13 stations at Pinyon Flats Observatory. The spatial footprint array diameters are 3 km or less, and each array was operational for a year or more. From our base catalog (M2.5+ earthquakes; 4038 events; September 2010 - June 2023), automated methods remove temporally overprinted waveforms from nearby events (< 5 km) in quick succession (< 5 min) and data with nonviable waveforms. These 200 samples per second data are filtered at 0.5-25 Hz and must have signal-tonoise ratios (SNRs) of 2.5+. Peak ground acceleration (PGA) and peak ground velocity (PGV) are derived individually from the maximum absolute values of each of the 3-component waveforms (vertical, northsouth, and east-west). Five of the seven arrays traverse the San Jacinto fault, and two do not. Ground motion observations are compared with theoretical estimates from Abrahamson et al., 2014. On average, arrays deployed within and across fault zones consistently record ground motions above theoretical expectations, whereas off-fault arrays record ground motions at or slightly below theoretical expectations. We attribute these differences to site conditions because these trends prevail for the full data suites. For each network and each individual channel, the coefficient of variation indicates that the standard deviations are ~30±6% of the mean. Exploring relative ground motion contributions from all three channels (ternary plots), as expected, most data show that vertical ground motions are attenuated compared to horizontal ones. However, this is not always the case for RA array data, where vertical motions can be ~2-3 times larger than horizontal motions for select events near Cahuilla, CA. These anomalously high vertical motions are focal mechanism-related. These results suggest that ternary plots created using only a small amount of data can be used as a data quality metric and a tool to find anomalous features in three-component data.

58 GEOSCIENCES↗
DOE OSTI2020

Modulation of seismic noise near the San Jacinto fault in southern California: origin and observations of the cyclical time dependence and associated crustal properties

SUMMARY We examine the cyclic amplitude variation of seismic noise recorded by continuous three-component broad-band seismic data with durations spanning 91–713 d (2008–2011) from three different networks: Anza seismic network, IDA network and the transportable seismic array. These stations surround the San Jacinto Fault Zone (SJFZ) in southern California. We find the seismic noise amplitudes exhibit a cyclical variation between 0.3 and 7.2 Hz. The high-frequency (≥0.9 Hz) noise variations can be linked to human activity and are not a concern. Our primary interest is signals in the low frequencies (0.3–0.9 Hz), where the seismic noise is modulated by semi-diurnal tidal mode M2. These long-period (low-frequency) variations of seismic noise can be attributed to a temporal change of the ocean waves breaking at the shoreline, driven by ocean tidal loading. We focus on the M2 variation of seismic noise at f = 0.6 Hz, travelling distances of ∼92 km through the crust from offshore California to the inland Anza, California, region. Relative to the shoreline station, data from the inland stations show a phase lag of ∼ –12°, which we attribute to the cyclic change in M2 that can alter crustal seismic attenuation. We also find that for mode M2 at 0.6 Hz, the amplitude variations of the seismic quality factor (Q) depend on azimuth and varies from 0.22 per cent (southeast to northwest) to 1.28 per cent (northeast to southwest) with Q = 25 for Rayleigh waves. We propose the direction dependence of the Q variation at 0.6 Hz reflects the preferred orientation of subfaults parallel to the main faulting defined by the primarily N45°W strike of the SJFZ.

Geochemistry & Geophysics↗
DOE OSTI2024

Evaluating Physics-Informed Neural Network Performance for Seismic Discrimination between Earthquakes and Explosions

In this article, we evaluate adding a weak physics constraint, that is, a physics‐based empirical relationship, to the loss function with a physics‐informed manner in local distance explosion discrimination in the hope of improving the generalization capability of the machine learning (ML) model. We compare the proposed model with the two‐branch model we previously developed, as well as with a pure data‐driven model. Unexpectedly, the proposed model did not consistently outperform the pure data‐driven model. By varying the level of inconsistency in the training data, we find this approach is modulated by the strength of the physics relationship. In conclusion, this result has important implications for how to best incorporate physical constraints in ML models.

58 GEOSCIENCES↗
DOE OSTI2022

Earthquake Phase Association Using a Bayesian Gaussian Mixture Model

Earthquake phase association algorithms aggregate picked seismic phases from a network of seismometers into individual seismic events and play an important role in earthquake monitoring and research. Dense seismic networks and improved phase picking methods produce massive seismic phase datasets, particularly for earthquake swarms and aftershocks occurring closely in time and space, making phase association a challenging problem. Here, we present a new association method, the Gaussian Mixture Model Association (GaMMA), that combines the Gaussian mixture model with earthquake location, origin time, and magnitude estimation. We treat earthquake phase association as an unsupervised clustering problem in a probabilistic framework, where each earthquake corresponds to a cluster of P and S phases with a hyperbolic moveout of arrival times and a decay of amplitude with distance. We use the multivariate Gaussian distribution to model the collection of phase picks of an event; and the mean of the multivariate Gaussian distribution is given by the predicted arrival time and amplitude from the causative event. We carry out the pick assignment to each earthquake and determine earthquake source parameters (i.e., earthquake location, origin time, and magnitude) under the maximum likelihood criterion using the Expectation-Maximization algorithm. The GaMMA method does not require typical association steps of other algorithms, such as grid-search or supervised training. The results for both synthetic tests and for the 2019 Ridgecrest earthquake sequence show that GaMMA effectively associates phases from a temporally and spatially dense earthquake sequence while producing useful estimates of earthquake location and magnitude.

58 GEOSCIENCES↗
DOE OSTI2022

Deep compressed seismic learning for fast location and moment tensor inferences with natural and induced seismicity

Fast detection and characterization of seismic sources is crucial for decision-making and warning systems that monitor natural and induced seismicity. However, besides the laying out of ever denser monitoring networks of seismic instruments, the incorporation of new sensor technologies such as Distributed Acoustic Sensing (DAS) further challenges our processing capabilities to deliver short turnaround answers from seismic monitoring. In response, this work describes a methodology for the learning of the seismological parameters: location and moment tensor from compressed seismic records. In this method, data dimensionality is reduced by applying a general encoding protocol derived from the principles of compressive sensing. The data in compressed form is then fed directly to a convolutional neural network that outputs fast predictions of the seismic source parameters. Thus, the proposed methodology can not only expedite data transmission from the field to the processing center, but also remove the decompression overhead that would be required for the application of traditional processing methods. An autoencoder is also explored as an equivalent alternative to perform the same job. We observe that the CS-based compression requires only a fraction of the computing power, time, data and expertise required to design and train an autoencoder to perform the same task. Implementation of the CS-method with a continuous flow of data together with generalization of the principles to other applications such as classification are also discussed.

54 ENVIRONMENTAL SCIENCES↗