Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic source”

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 109 records · Page 6

Imaging the Shallow Structure of the Yucca Flat at the Source Physics Experiment Phase II Site with Horizontal-to-Vertical Spectral Ratio Inversion and a Large- N Seismic Array

The Source Physics Experiment (SPE) is a series of chemical explosions at the Nevada National Security Site (NNSS) with the goal of understanding seismic-wave generation and propagation of underground explosions. To understand explosion source physics, accurate geophysical models of the SPE site are needed. Here, we utilize a large-N seismic array deployed at the SPE phase II site to generate a shallow subsurface model of shear-wave velocity. The deployment consists of 500 geophones and covers an area of, approximately, 2.5 × 2 km. The array is located in the Yucca Flat in the northeast corner of the NNSS, Nye County, Nevada. Using ambient-noise recordings throughout the large-N seismic array, we calculate horizontal-to-vertical spectral ratios (HVSRs) across the array. We obtain 2D seismic images of shear-wave velocities across the SPE phase II site for the shallow structure of the basin. In this work the results clearly image two significant seismic impedance interfaces at ~150–500 and ~350–600 m depth. The shallower interface relates to the contrast between Quaternary alluvium and Tertiary volcanic rocks. The deeper interface relates to the contrast between Tertiary volcanic rocks and the Paleozoic bedrock. The 2D subsurface models support and extend previous understanding of the structure of the SPE phase II site. This study shows that the HVSR method in conjunction with a large-N seismic array is a quick and effective method for investigating shallow structures.

58 GEOSCIENCES↗

Humming Trains in Seismology: An Opportune Source for Probing the Shallow Crust

Abstract Seismologists are eagerly seeking new and preferably low-cost ways to map and track changes in the complex structure of the top few kilometers of the crust. By understanding it better, they can build on what is known regarding important, practical issues. These include telling us whether imminent earthquakes and volcanic eruptions are generating telltale underground signs of hazard, about mitigation of induced seismicity such as from deep injection of wastewater, how the Earth and its atmosphere couple, and where accessible natural resources are. Passive seismic imaging usually relies on blind correlations within extended recordings of Earth’s ceaseless “hum” or coda of well-mixed, small vibrations. In this article, we propose a complementary approach. It is seismic interferometry using opportune sources—specifically ones not stationary in time and moving in a well-understood configuration. Its interpretation relies on an accurate understanding of how these sources radiate seismic waves, precise timing, careful placement of pairs of listening stations, and seismic phase differentiation (surface and body waves). Massive freight trains were only recently recognized as such a persistent, powerful cultural (human activity-caused) seismic source. One train passage may generate a tremor with an energy output of a magnitude 1 earthquake and be detectable for up to 100 km from the track. We discuss the source mechanisms of train tremors and review the basic theory on sources. Finally, we present case studies of body- and surface-wave retrieval as an aid to mineral exploration in Canada and to monitoring of a southern California fault zone. We believe noise recovery from this new signal source, together with dense data acquisition technologies such as nodes or distributed acoustic sensing, will deeply transform our ability to monitor activity in the shallow crust at sharpened resolution in time and space.

Geochemistry & Geophysics↗

Seismic Event Characterization Using Full Moment Tensors on the Hypersphere

Moment tensor solutions provide insights into the deformation that has occurred in the source region of a seismic event and are therefore of great value in identifying different types of seismic sources, such as when monitoring for underground nuclear tests. Despite this utility, inversion of waveforms recorded by seismometers for their full seismic moment tensor is not yet routine, and development of robust methods to classify events based on this information is in its infancy. Here, we assemble an inventory of 1405 full moment tensor solutions that include explosive, earthquake, and collapse events, and investigate the use of anisotropic probability distribution functions on the 5D hypersphere to discriminate between these sources. Using a Bayesian classifier, we obtain optimal success rates of 98.4% across all events and demonstrate that modification of the prior probabilities provides a natural way to alter the balance between not missing desirable events (such as explosions) versus misclassifying large numbers of undesired events (such as earthquakes). The approach is specifically designed to progress from traditional, bipolar event screening metrics to more generalized event identification across multiple types of seismic sources. Despite current databases containing insufficient numbers of events to definitively demonstrate at present, we also find intriguing evidence of subgroupings within individual source populations on the hypersphere, for example, between chemical and nuclear explosions, raising the potential possibility of discriminating between these event types in the future.

Geosciences↗

New Architecture to Support Integration and Processing of Seismic Data from Heterogeneous Sources

The Geophysical Monitoring Program (GMP) at Lawrence Livermore National Lab (LLNL) maintains a database and supporting infrastructure for geophysical data used in support of the Nuclear Detonation Detection mission. This database includes data from multiple sources, many of which do not distribute data to the public or for which there is no automated means of access. For example, Figure 1 shows (left) the distribution of waveform data in our database by source. The Incorporated Research Institutions for Seismology Data Management Center (IRISDMC) is our major source of waveform data and those data may be retrieved at will using the Federated Digital Seismograph Networks FDSN web Application Programming Interface (API). However, the next 6 most important sources of waveform data have no or only limited automated access to waveforms. As Figure 1 (right) shows, it is very common for waveform records associate with an event in our database to come from two or more sources, and in some cases data come from 10 sources. This diversity of data sources drives our need for efficient and correct integration of metadata, parametric data, and waveform data.

58 GEOSCIENCES↗

Unraveling the Wrinkle in Time-Variable Sources with Lunes and Synthetic Seismic Data

In this report, we describe how to estimate the time-variable components of the seismic moment tensor and compare these estimates to the more conventional analysis that incorporates an assumption of the source time function (STF) across all components of the seismic moment tensor. The advantage of our method is that we are able to independently estimate the time-evolution of each component of the seismic moment tensor, which may help to resolve the complex source phenomena associated with buried explosions. By performing an eigen decomposition of the time-evolving seismic moment tensor components, we are able to plot the seismic mechanism as a trajectory on a lune diagram. This technique enables interpretation of the seismic mechanism as a function of time, as opposed to the more conventional analysis which assumes that the seismic mechanism is time invariant. Finally, we describe the differences between the seismic moment and the seismic moment rate STFs, how to implement each one in inversion schemes, and the relative strengths/weaknesses of each. Our key take-away is that we are able to distinguish nearly-overlapping sources with highly different mechanisms, such as an explosion immediately following an earthquake, by estimating moment rate from seismic data through a STF-invariant inversion for the full time-variable moment tensor.

58 GEOSCIENCES↗

Estimating an airborne dipole source using 3-D wavefield simulations and seismic receivers on the ground

Acoustic signals can couple to the ground, providing an opportunity to use seismic stations to investigate airborne sources. The study of Bishop et al. used wavefield simulations in a fluid-solid medium to quantify the role of topography on the seismic (ground) recordings of a monopole source in the air. We build upon this study by linking wavefield forward modelling with the source estimation code MTUQ, which can accommodate point forces or moment tensors in a solid medium, as well as sources in the air (or water) if they are enabled by the forward-modelling solver. We perform a series of synthetic numerical experiments to demonstrate that a dipole airborne source can be estimated using ground-based receivers, even within the presence of realistic topography. We investigate the influence of receiver coverage, topography, and assumed source location on the estimated results. The established capabilities raise the prospects for future efforts to estimate dipole sources in 3-D models that include heterogeneity in the air and the earth in addition to topography.

58 GEOSCIENCES↗

The 2017 Pohang, South Korea, Mw 5.4 main shock was either natural or triggered, but not induced

Understanding the cause of the November 2017 Pohang main shock is of considerable importance to the future of the geothermal industry because of its large magnitude compared to prior expectations based on case histories of other projects involving underground fluid injection. Of the three possibilities – induced, triggered or natural, “induced” can be ruled out based on the disproportionately large seismic moment of the main shock. Furthermore, whether natural or triggered, the source of seismic hazard at Pohang was tectonic strain accumulation, not fluid injection. Arguably, the most timely indicator of seismic hazard and risk in the environs of Pohang was the September 2016 Mw 5.4 Gyeongju earthquake, which was natural and located about 40 km south of Pohang along the same active fault system.

58 GEOSCIENCES↗

FY 2025 End of Year Report: Seismic Monitoring of Underground Vibration Sources using Distributed Acoustic Sensing (DAS) and Seismometers

This end-of-year report summarizes progress on using seismic monitoring to detect, associate, and locate anomalous vibration signals that may indicate potential containment breaches. The work focused on four key tasks: 1. Developing a database of continuous waveforms and ground-truth event data from multiple sensing modalities. 2. Refining and implementing detection and association algorithms to generate a catalog of anomalous underground activities. 3. Testing and improving distributed acoustic sensing amplitude-based geolocation methods to build an event location catalog. 4. Testing and refining seismic array polarization-based geolocation methods to build an event location catalog. This report provides a brief recap of results from the FY25 midyear report (Tasks 1 and 2) and presents new findings from geolocation methods (Tasks 3 and 4).

58 GEOSCIENCES↗

FY 2026 Midyear Report: Seismic Monitoring of Underground Vibration Sources Using Distributed Acoustic Sensing and Seismometers

Safeguards-relevant temporal changes in underground facilities can be observed using geophysical monitoring techniques. Seismic waves, in particular, provide valuable insights into subsurface activities and can serve as an important tool for detecting anomalous events that may indicate containment breaches at geological repositories. This midyear report summarizes ongoing efforts to automatically and rapidly detect and locate anomalous vibration signals that could be indicative of potential containment breaches. Previous work during FY25 focused on compiling continuous seismic datasets from two underground sites and developing a database of continuous waveforms and ground-truth event data derived from multiple sensing modalities. Building on this foundation, we are adapting anomaly detection and geolocation algorithms to explore methods for monitoring underground activities using two relatively low-maintenance sensing technologies: a dense surface geophone array deployed at the Pleasant Gap mine in Pennsylvania, and a three-dimensional fiber-optic cable array for distributed acoustic sensing (DAS) installed in the subsurface at the Sanford Underground Research Facility (SURF) in South Dakota. This report summarizes work conducted during the first two quarters of FY26, during which we refined a dynamic power spectral density (PSD)-based detector, applied it independently to each geophone station, and then combined the per‑station detections with density-based spatial clustering of applications with noise (DBSCAN) to cluster events and produce spatial maps over a nine‑day interval. In addition, we outline plans for a field trial at the Waste Isolation Pilot Plant (WIPP) in New Mexico to compare traditional seismic monitoring approaches with DAS techniques and to evaluate the benefits of combined data analysis. Activities during the past two quarters have included the preparation and submission of a Field Test Plan to WIPP for approval, as well as submission to headquarters for review and feedback.

58 GEOSCIENCES↗

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↗

Screening IMS Seismic Detections using Dynamic Correlation Processor: Progress report for 2020

LLNL completed the FY20 workplan entailing application of the Dynamic Correlation Processor (DCP) to International Monitoring System (IMS) primary seismic data in 2019 and 2020. The effective analysis period was curtailed to January 1 through November 12, 2019 due to a data compression change that rendered later data unusable. Fixing the data compression error is a priority for our data source, and we plan to extend the period of data analysis when the fix is implemented. We found that the percentage overlap between detections generated by the International Data Centre (IDC) DFX software and detections generated by DCP varied widely by station. Overlap is 30-40% at some stations, suggesting that a large percentage of DFX detections could be screened (removed) from the detection-event association algorithm. However, detection overlap was much lower, a few percent, at many stations and the percentage overlap was low at stations known to have many mines and other local sources of seismicity. Unexpectedly low detection overlap at many stations prompted us to examine results for the ARCES station in detail. There is considerable mine activity near ARCES, but our initial FY20 analysis found only ~2% DCP/DFX detection overlap. Differences in the pre-processing (e.g. beam recipes) used by DCP and DFX appear to be the cause of low detection overlap at ARCES. DCP uses wideband filters to improve detector robustness and the IDC uses narrow-band pre-filtering to improve sensitivity. After reprocessing ARCES using the IDC beam recipe, DCP/DFX overlap increased to ~74%. This result shows that we must reprocess the IMS network using the IDC station-specific beam recipes if we are to effectively screen DFX detections and improve automatic event building performance at the IDC.

58 GEOSCIENCES↗

Assessing the Accuracy of a Borehole-Controlled Velocity Model of Yucca Flat, Nevada, Using Large-N Seismic Data

The Source Physics Experiment (SPE) aims to improve nuclear test monitoring capabilities through a series of chemical explosions in Yucca Flat within the Nevada National Security Site (NNSS) (Snelson et al., 2013). With geologic data from over 950 boreholes, Yucca Flat has excellent borehole control on stratigraphy. These data were used to create a Geologic Framework Model (GFM) of the basin to aid in the study of seismic wave propagation through the Earth. However, the seismic velocities of the geologic units in the GFM are less well-defined. With the acquisition of five new active source Large-N datasets, we can now quantitatively assess the accuracy of the seismic model previously controlled by borehole data alone.

58 GEOSCIENCES↗

Utah FORGE 3-2535: Joint EM-Seismic-InSAR Imaging of Fracture Properties - Workshop Presentation

This is a presentation on the Joint Electromagnetic/Seismic/InSAR Imaging of Spatial-Temporal Fracture Growth and Estimation of Physical Fracture Properties During EGS Resource Development project by Lawrence Berkeley National Laboratory, presented by Dr. David Alumbaugh, Staff Scientist. The project's objective was to develop a set of technologies and workflow to image induced fracture generation and growth for an Enhanced Geothermal System (EGS). The project anticipates imaging of the fracture generation and growth at FORGE using a combination of passive seismic, active source borehole EM, and INSAR technology. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗

Multiwell Fiber Optic Sensing Reveals Effects of CO2 Flow on Triggered Seismicity

Induced seismicity is one of the main risks for gigaton-scale geological storage of carbon dioxide (⁠CO 2 ⁠). Thus, passive seismic monitoring is often recommended as a necessary component of the monitoring systems for CO 2 storage projects, with a particular forcus on risk mitigation. We present the first field study, CO2CRC Otway Project Stage 3 (Victoria, Australia), where distributed acoustic sensing (DAS) enabled high-precision tracking of the induced seismicity triggered by a small CO 2 injection and also informed the reservoir models. In 610 days of passive seismic monitoring of the Stage 3 injection, we detected 17 microseismic events (maximum moment magnitude Mw 0.1) using five deep boreholes equipped with enhanced-sensitivity optical fiber. The DAS array has sensitivity sufficient for detection and location of induced events with M w ~ ₋2 in a monitoring borehole located up to 1500 m away. Thanks to the dense spatial sampling by the DAS, we were able to estimate the focal mechanisms for events with M w > ₋1.5; although the monitoring boreholes provided very limited angular coverage. The main cluster of the events has the same location and source mechanism as the one triggered by the previous CO 2 injection at the Otway Project site, Stage 2C. Surprizingly, the Stage 2C and Stage 3 events closely followed the actual movement of the CO 2 saturation plume front (not the pressure front), as observed using controlled-source reflection seismic images. The nature of the plume-fault interaction remains unclear, but some alteration of the fault gouge by CO 2 might be responsible for the faults’ reactivation by the pressure perturbation. Importantly, the seismogenic fault could not be identified in the seismic images and was only revealed by DAS observations, which also demonstrated the signature of fluid–rock interaction, that may control the CO 2 flow.

58 GEOSCIENCES↗

3D characterisation and quantification of an offshore freshened groundwater system in the Canterbury Bight

Although offshore freshened groundwater (OFG) systems have been documented in numerous continental margins worldwide, their geometry, controls and emplacement dynamics remain poorly constrained. Here we integrate controlled-source electromagnetic, seismic reflection and borehole data with hydrological modelling to quantitatively characterise a previously unknown OFG system near Canterbury, New Zealand. The OFG system consists of one main, and two smaller, low salinity groundwater bodies. The main body extends up to 60 km from the coast and a seawater depth of 110 m. We attribute along-shelf variability in salinity to permeability heterogeneity due to permeable conduits and normal faults, and to recharge from rivers during sea level lowstands. A meteoric origin of the OFG and active groundwater migration from onshore are inferred. However, modelling results suggest that the majority of the OFG was emplaced via topographically-driven flow during sea level lowstands in the last 300 ka. Global volumetric estimates of OFG will be significantly revised if active margins, with steep coastal topographies like the Canterbury margin, are considered.

54 ENVIRONMENTAL SCIENCES↗

Seismic constraints from a Mars impact experiment using InSight and Perseverance

NASA’s InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) mission has operated a sophisticated suite of seismology and geophysics instruments on the surface of Mars since its arrival in 2018. On 18 February 2021, we attempted to detect the seismic and acoustic waves produced by the entry, descent and landing of the Perseverance rover using the sensors onboard the InSight lander. Similar observations have been made on Earth using data from both crewed and uncrewed spacecraft, and on the Moon during the Apollo era, but never before on Mars or another planet. This was the only seismic event to occur on Mars since InSight began operations that had an a priori known and independently constrained timing and location. It therefore had the potential to be used as a calibration for other marsquakes recorded by InSight. Here we report that no signal from Perseverance’s entry, descent and landing is identifiable in the InSight data. Nonetheless, measurements made during the landing window enable us to place constraints on the distance–amplitude relationships used to predict the amplitude of seismic waves produced by planetary impacts and place in situ constraints on Martian impact seismic efficiency (the fraction of the impactor kinetic energy converted into seismic energy).

79 ASTRONOMY AND ASTROPHYSICS↗

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗