Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic sources”

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 145 records · Page 8

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↗

Reverse Time Migration of Explosion-generated Infrasound and Ground-Coupled Airwaves

Accurately locating sources of seismic and infrasonic energy is integral to global monitoring of earthquakes and explosions. Infrasound arrivals times can be used to calculate the origins of events that generate acoustic energy. Picking times of emergent infrasound arrivals, however, can be difficult and prone to uncertainty. Reverse time migration (RTM) is a waveform based location method that does not rely on picked arrival times. Here we use RTM to locate a known chemical explosion that generated acoustic and acoustic-to-seismic signals on 26 and 108 receivers, respectively. All location predictions are less than 24 km from the known location with time errors of less than three minutes. We find strong overall agreement between our results and those of existing RTM and arrival time based methods. Our initial results suggest that RTM is a promising method of event location using acoustic arrivals recorded on both infrasound and seismic instrumentation.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Designing Yield Grids for Airborne Nuclear Explosions to Maintain Decrements in Rayleigh Waveform Amplitudes

This document details a method for selecting grid points that sample the yields of airborne, low yield nuclear explosions that constrains the relative size of seismic waveforms sourced by explosions that are parameterized by this grid. In detail, we impose this constraint so that the vertical component of ground displacement that is sourced by explosions with a fixed height of burst (HoB) at grid point k change by a fixed factor of a when compared to the waveform amplitudes sourced by explosions at grid point k + 1, as observed at a fixed standoff distance. This gridding thereby preserves amplitude ratios along the grid so that they are agnostic to k. We assert that this method is useful in tests that apply algorithms against synthetic or augmented datasets for which the performance of the algorithm must be measured against source size, and waveform amplitudes must change uniformly (e.g., Berg et al, 2025).

58 GEOSCIENCES↗

Three-Dimensional Magnetic Restructuring in Two Homologous Solar Flares in the Seismically Active NOAA AR 11283

We carry out a comprehensive investigation comparing the three-dimensional magnetic field restructuring, flare energy release, and the helioseismic response of two homologous flares, the 2011 September 6 X2.1 (FL1) and September 7 X1.8 (FL2) flares in NOAA AR 11283. In our analysis, (1) a twisted flux rope (FR) collapses onto the surface at a speed of 1.5 km s(exp−1) after a partial eruption in FL1. The FR then gradually grows to reach a higher altitude and collapses again at 3 km s(exp−1) after a fuller eruption in FL2. Also, FL2 shows a larger decrease of the flux-weighted centroid separation of opposite magnetic polarities and a greater change of the horizontal field on the surface. These imply a more violent coronal implosion with corresponding more intense surface signatures in FL2. (2) The FR is inclined northward and together with the ambient fields, it undergoes a southward turning after both events. This agrees with the asymmetric decay of the penumbra observed in the peripheral regions. (3) The amounts of free magnetic energy and nonthermal electron energy released during FL1 are comparable to those of FL2 within the uncertainties of the measurements. (4) No sunquake was detected in FL1; in contrast, FL2 produced two seismic emission sources S1 and S2 both lying in the penumbral regions. Interestingly, S1 and S2 are connected by magnetic loops, and the stronger source S2 has a weaker vertical magnetic field. We discuss these results in relation to the implosion process in the low corona and the sunquake generation.

Sun: X-rays↗

Combining translational and rotational seismic motions to invert local-scale seismic data for time-variable moment tensors: do rotational motions help for high-frequency seismic data produced by underground explosions?

SUMMARY We present an analysis of combining translational and rotational seismic data in an inversion for the time-variable source time functions corresponding to the components of the seismic moment tensor. We conduct a series of numerical experiments where the data are simulated by a combination of an underground explosion and a co-located double couple shear source and recorded on surface-mounted seismometers within 1–2 km of the source. The experiments are designed to mimic explosion seismology experiments, and thus the data are in the 1–10 Hz frequency range and contain very few surface waves. We use a Monte Carlo method to propagate Earth model uncertainty into the estimates of seismic source parameters. In our experiments, we find that the uncertainty of the estimated seismic source parameters increases when we add rotational seismic motions to the inversion when using a constant number of data channels. In this case, the increased degree of uncertainty in the final results is most likely due to the near-surface Earth model uncertainty that we introduce in our simulations. However, for a fixed number of seismic stations, adding rotational seismic motions to the inversion acts to decrease the uncertainty of the estimated seismic source parameters, most likely due to the increase in the number of data channels used in the inversion.

Poppeliers, Christian (ORCID:0000000159526849)↗

HFT events - Shallow moonquakes

A few large distant seismic events of distinctly high signal frequency, designated HFT (high-frequency teleseismic) events, are observed yearly by the Apollo lunar seismic network. Their sources are located on or near the surface of the moon, leaving a large gap in seismic activity between the zones of HFT sources and deep moonquakes. No strong regularities are found in either their spatial or temporal distributions. Several working hypotheses for the identity of these sources have advanced, but many characteristics of the events seem to favor a hypothesis that they are shallow moonquakes. Simultaneous observations of other lunar phenomena may eventually enable the determination of their true identity.

Nakamura, Y.↗

Transportable absolute yields of underground nuclear explosions. Application to the North Korean nuclear tests

The method developed and presented here improves our nuclear explosion monitoring capability by providing a new independent and transportable method to determine the explosive yields of seismically recorded events associated with foreign weapons development. Here we report new independent absolute yield and depth of burial estimates of the North Korean nuclear tests which are commensurate with previously determined source parameters. The key advantage of this method over other methods which estimate absolute explosive yields is that this method does not require any a priori calibration, and can be immediately applied to any region of interest. This method uses the source information obtained from the spectral ratios of envelopes of measured seismic coda waves to invert for model-derived yield and depth of burial. The only requirements of this method are that the region of interest must contain several (> 3) seismic sources which have been well recorded (SNR > 2) by a set of shared stations. For regions with only a few events (~ 3) and/or extremely band-limited observations for which the depth of burial is unknown, an earthquake with a known magnitude may be required to obtain an accurate and robust yield estimate.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

The 26 December 2004 tsunami source estimated from satellite radar altimetry and seismic waves

The 26 December 2004 Indian Ocean tsunami was the first earthquake tsunami of its magnitude to occur since the advent of both digital seismometry and satellite radar altimetry. Both have independently recorded the event from different physical aspects. The seismic data has then been used to estimate the earthquake fault parameters, and a three-dimensional ocean-general-circulation-model (OGCM) coupled with the fault information has been used to simulate the satellite-observed tsunami waves. Here we show that these two datasets consistently provide the tsunami source using independent methodologies of seismic waveform inversion and ocean modeling. Cross-examining the two independent results confirms that the slip function is the most important condition controlling the tsunami strength, while the geometry and the rupture velocity of the tectonic plane determine the spatial patterns of the tsunami.

ocean general circulation model (OGCM)↗

Dynamic Rupture Simulations of the Mw7.2 1992 Landers,California, Earthquake

An important part of the simulation technique implemented in the Recipe is the earthquake rupture characterization based on established empirical scaling relationships between source parameters and seismic moment. The scaling relationships used in the rupture characterization are derived from recorded data and are tested against dynamic and kinematic rupture models of crustal earthquakes. Currently, a revised three-stage source scaling relationship between the source rupture area and seismic moment of inland crustal earthquakes, initially proposed by Irikura and Miyake, is being used. Recent analysis of results obtained with the Recipe for the 2016 M7.1 Kumamoto earthquake performed by several researchers, including the PI, Arben Pitarka, indicate that in order to reproduce the near-fault long-period ground motions and the permanent displacements, observed during surface-rupturing earthquakes, the recipe needs modifications to the kinematic rupture parametrization in the shallow crustal weak zone. In particular, improvements are required in characterizing the slip rate function at shallow depths, above the sesimogenic zone. The main objective of this research project was the development of a dynamic rupture model for the Mw7.2 1992 Landers, California earthquake. The dynamic rupture model was used to derive kinematic rupture parameters, such as slip, slip-rate, rupture time distributions, and slip rate functions in the weak (0-4km) and seismogenic crustal zones. In addition to the knowledge gained from the dynamic rupture modeling of Chi-chi, Ridgecrest, Kumamoto and other crustal earthquakes, studied by the PI, the rupture characterization of the Landers earthquake can be helpful in updating the conventional recipe for broadband ground motion simulations of crustal earthquakes with shallow rupture. The project emphasizes the need for performing physics-based earthquake rupture modeling which ultimately will inform further improvements of rupture characterization. In particular, the modeling of rupture dynamics for recorded earthquakes, and numerical simulations for different rupture scenarios, provides physical constrains on slip evolution in time and space, used in adjustments of empirical source scaling relationships for large crustal earthquakes.

58 GEOSCIENCES↗

BayesMT: A Probabilistic Bayesian Framework for the Seismic Moment Tensor

Moment tensors (MTs) have long been used in earthquake and explosion source analysis, and there has been a renewed interest in how they can inform us about the seismic source, particularly in the geophysical monitoring community due to its application in event identification and yield analysis. However, parameter uncertainties in seismic MT inversion are rarely available. The inverse procedure often does not quantify MT model errors such as event location, data noise and Earth model that are essential for estimating solution robustness. To address this need, we propose to adopt the Bayesian probabilistic framework to incorporate uncertainties in MT inversions. In this study, we present the theoretical background of a probabilistic Bayesian framework for MT inversion accounting for model and measurements errors and illustrate the implementation of the method using a synthetic example.

58 GEOSCIENCES↗

Programmatic Advantages of Linear Equivalent Seismic Models

Underground explosions nonlinearly deform the surrounding earth material and can interact with the free surface to produce spall. However, at typical seismological observation distances the seismic wavefield can be accurately modeled using linear approximations. Although nonlinear algorithms can accurately simulate very near field ground motions, they are computationally expensive and potentially unnecessary for far field wave simulations. Conversely, linearized seismic wave propagation codes are orders of magnitude faster computationally and can accurately simulate the wavefield out to typical observational distances. Thus, devising a means of approximating a nonlinear source in terms of a linear equivalent source would be advantageous both for scenario modeling and for interpretation of seismic source models that are based on linear, far-field approximations. This allows fast linear seismic modeling that still incorporates many features of the nonlinear source mechanics built into the simulation results so that one can have many of the advantages of both types of simulations without the computational cost of the nonlinear computation. In this report we first show the computational advantage of using linear equivalent models, and then discuss how the near-source (within the nonlinear wavefield regime) environment affects linear source equivalents and how well we can fit seismic wavefields derived from nonlinear sources.

58 GEOSCIENCES↗