Search NASA⌕ Search

SEARCH · Search NASA

Results for “Surface acoustic wave”

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

Search for Infrasound Signals in InSight Data Using Coupled Pressure/Ground Deformation Methods

The unprecedented quality and sampling rate of seismometer and pressure sensors of the InSight Mars mission allow us to investigate infrasound through its pressure and ground deformation signals. This study focuses on compliance effects induced by acoustic waves propagating almost horizontally close to the surface. The compliance of acoustic waves is first estimated using the compliance estimates from pressure perturbations moving at wind speed. Then, a marker of compliance events is used to detect events of ground deformation induced by pressure variations, in three frequency bands from 0.4 to 3.2 Hz, from InSight sol 180 to 690. Additional selection criteria are imposed on the detected events to focus on acoustic waves and to remove various noise sources (e.g., wind effects or seismometer artifacts). After an automated selection, the visual inspection of the records allows us to validate two infrasound candidates that cannot be related to pressure perturbations moving at wind speed nor to known noise sources. For our highest quality infrasound candidate, the relation between this event and a convective vortex occuring 10 s later is tested. The azimuth of the vortex position at the time of infrasound detection is not consistent with the arrival azimuth of the suspected infrasound inferred from the polarization of seismometer records, thus the link between these two phenomena cannot be demonstrated. Further investigations would require a better understanding of wind-related noise impacting InSight sensors and of the effects of lateral variations of subsurface mechanical properties on the ground deformations induced by atmospheric pressure variations.

58 GEOSCIENCES↗
DOE OSTI2023

Seismoacoustic Analysis of the Large Surface Explosion Coupling Experiment Using a Large- N Seismic Array

Seismoacoustic signals at local distance (<~10 km) are widely used as important constraints on source parameters for near-surface events, yet the seismoacoustic wave generation and energy partitioning are not fully understood. Spatially dense sensors could provide observations in high resolution to capture the full wavefield for better understanding wave propagation and improving source estimation. Recently, spatially dense observations of the local seismoacoustic wavefield produced by a pair of 1-ton surface explosions have been recorded using a large-N seismic array. This large-N array consists of 446 geophones and covers an area of approximately 2×2.5 km 2 . The two surface explosions occur at the same location but at different times with different atmospheric conditions. Both seismic and air–ground coupled acoustic waves from the two surface explosions are well observed. Analyses of signals recorded by the large-N seismic array show different acoustic wave speed and amplitude for the two explosions. A strong spatial variability in acoustic wave speed and amplitude for each explosion is also observed. Finally, the observations suggest the important role of local atmosphere state on wave propagation and source estimation and demonstrate how the use of a large-N capability can improve characterization of the propagation medium and source.

58 GEOSCIENCES↗
DOE OSTI2022

Identification of Acoustic Wave Signatures in the Ionosphere From Conventional Surface Explosions Using MF/HF Doppler Sounding

Here, we present an experiment to detect one ton TNT-equivalent chemical explosions using pulsed Doppler radar observations of isodensity layers in the ionospheric E region during two campaigns. The first campaign, conducted on 15 October 2019, produced potential detections of all three shots. The detections closely resemble the temporal and spectral properties predicted using the InfraGA ray tracing and weakly nonlinear waveform propagation model. Here the model predicts that within 6.5–7.25 min of each shot a waveform peaking between 0.9 and 0.4 Hz will impact the ionosphere at 100 km. As the waves pass through this region, they will imprint their signal on an isodensity layer, which is detectable using a Doppler radar operating at the plasma frequency of the isodensity. Within the time windows of each of the three shots in the first campaign, we detect enhanced wave activity peaking near 0.5 Hz. These waves were imprinted on the Doppler signal probing an isodensity layer at 2.785 MHz near 100 km altitude. Despite these detections, the method appears to be unreliable as none of the six shots from the second campaign, conducted on 10 July 2020 were detected. The observations from this campaign were characterized by an increased acoustic noise environment in the microbarom band and persistent scintillation on the radar returns. These effects obscured any detectable signal from these shots and the baseline noise was well above the detection levels of the first campaign.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Infrasonic directivity of monopole, dipole and bipole ground-surface reflected sources

Infrasound (acoustic waves below 20 Hz) can be used to detect, locate and quantify activity in the atmosphere such as volcanic eruptions and anthropogenic explosions. Attempts to quantify volcanic eruption parameters such as exit velocity, plume height and mass flow rate using infrasound data depend strongly on assumptions of the acoustic source type. Infrasonic sources may produce omnidirectional or directional wavefields, while propagation effects, such as interaction with topography, can induce further wavefield directivity that is measured by field instrumentation. Limited sampling of these wavefields can hinder our ability to infer the underlying source, and thus our understanding of the eruption characteristics. Equivalent sources are often used to represent acoustic source mechanisms and resultant wavefields. In this study, we review equivalent acoustic sources as they pertain to infrasonic scale and wavelengths commonly encountered in very local (⁠<5 km range) geophysical field deployments. We highlight the equivalent infrasonic bipole source that can be induced by ground-reflection of an elevated monopole; we are not aware of any prior infrasound studies that use the bipole source concept. We use analytical and numerical methods to explore source directivity of monopole, dipole and bipole ground-reflected sources at infrasonic frequencies as well as the additional directivity complications introduced by interactions with topography. We illustrate that for typical volcano-infrasound wavelengths, increasing height above the ground as well as increasing source frequency leads to increased wavefield directivity. Numerical modelling using a simple omnidirectional monopole source embedded in topography further illustrates that both horizontal and vertical infrasound directionality can be induced by topography at the distance scales appropriate for local volcano infrasound monitoring. Information summarized in this analytical and numerical exploration of infrasound directivity may be used to help guide future volcano-infrasound field deployments intended to estimate source parameters or quantify wavefield directivity. Analytic solutions for simple whole-space or half-space atmospheres provide useful formulations for planning or initially analysing geophysical field-scale experimental data; however, especially at very local distances from the source (⁠<5 km), 3-D simulations are necessary to account for complex topography commonly encountered in volcano-infrasound applications.

Infrasound↗
DOE OSTI2020

Surface Josephson plasma waves in a high-temperature superconductor

Abstract Electron density oscillations with acoustic dispersions and sustained at boundaries between different media provide information about surface and interface properties of heterostructures. In ultrathin metallic films these plasmonic excitations are heavily damped. Superconductivity is predicted to reduce dissipation allowing detection of these resonances. Emerging low-loss interface Cooper-pair waves have been studied before, however, the observation of surface-confined Josephson plasmons in highly anisotropic superconductors has remained elusive. Here, we report on generation and coupling to these excitations in an ultrathin single-crystal film of high-temperature superconductor La 1.85 Sr 0.15 CuO 4 . The film becomes brighter than Au below the critical temperature when probed with sub-gap THz photons. We show that the enhanced signal in the superconducting state, which can be visualized with a spatial resolution better than λ/3000, originates from near-field coupling of light to surface Josephson plasmons. Our results open a path towards non-invasive investigation of enhanced superconductivity in artificial multilayers, buried interface states in topological heterostructures, and non-linear phenomena in Josephson devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Nonlinear ultrasonic technique for the characterization of microstructure in additive materials

This study employs nonlinear ultrasonic techniques to track microstructural changes in additively manufactured metals. The second harmonic generation technique based on the transmission of Rayleigh surface waves is used to measure the acoustic nonlinearity parameter, β. Stainless steel specimens are made through three procedures: traditional wrought manufacturing, laser-powder bed fusion, and laser engineered net shaping. The β parameter is measured through successive steps of an annealing heat treatment intended to decrease dislocation density. Dislocation density is known to be sensitive to manufacturing variables. In agreement with fundamental material models for the dislocation-acoustic nonlinearity relationship in the second harmonic generation, β drops in each specimen throughout the heat treatment before recrystallization. Geometrically necessary dislocations (GNDs) are measured from electron back-scatter diffraction as a quantitative indicator of dislocations; average GND density and β are found to have a statistical correlation coefficient of 0.852 showing the sensitivity of β to dislocations in additively manufactured metals. Moreover, β shows an excellent correlation with hardness, which is a measure of the macroscopic effect of dislocations.

42 ENGINEERING↗
DOE OSTI2025

Advancing Tunnel Detection Via Vertical Acoustic Profiling

This study demonstrates that placing an acoustic source at depth significantly enhances the detection of tunnels. By adapting the methodology of vertical seismic profiling to acoustic-range frequencies, we introduce a novel approach for tunnel detection. To our knowledge, this report documents the first application of this technique. The field study was conducted on a cleared site at the Oak Ridge National Laboratory; the subsurface was clay. A series of vertical shafts were bored immediately before and after the placement of a 40 ft long, 3 ft diameter steel pipe located 10 ft below the surface. Acoustic signatures were produced at various depths up to 30 ft and an array of surface geophones was used to collect the resulting signals. The innovation lies in using subsurface acoustic signatures to better enable the collection of wide-angle diffracted waves by surface geophones. These diffracted waves appear as subharmonic frequencies of the source signal and are recorded by a surface geophone array. While the concept draws from vertical seismic profiling used in oil and gas exploration, a key distinguishing feature is the use of acoustic frequencies, which are necessary for detecting small subsurface features. The findings have important implications. Critical information pertaining to diffracted source signals are not measured by current state-of-the-art methods. Although this study confirms enhanced detectability, accurate localization and imaging will require additional research. Future work should focus on leveraging diffracted-wave strength and time delays to enable full imaging and characterization.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Numerical Modeling of Air-Blast Suppression as a Function of Explosive-Charge Burial Depth

As a chemical explosion is buried, the mechanism for acoustic wave generation transitions from fully gas-generated at the surface to completely spall-induced at full containment depth. The fully gas-generated and completely spall-induced signals in the acoustic waveform are well described; however, the transition between these two end-members eludes numerical modeling because of the complex phenomena that are involved. The phenomena of crater formation and explosive cloud evolution are simulated using an Eulerian hydrocode that incorporates geomaterials with strength and porosity. Having accurately modeled these phenomena, we can confidently predict the propagation and relative strength of the gas-generated and spall-induced pulses in the recorded acoustic waveform. The numerical predictions agree with observations from the historical Stagecoach experiment as well as modern recordings from the Source Physics Experiment. In particular, the peak pressure p generated by an explosion is initially due to the gas-generated mechanism and decays with scaled depth of burial d s (depth d scaled by the cube-root of explosive yield w 1/3 ) as exp(-d s ) but then transitions near a scaled depth of 6 m/ton 1/3 to the spall-generated mechanism in which the decay is d$_{s}^{-7/4}$. This decay form is related to the strong ground-motion attenuation relationship that affects spall strength. So these results can improve seismoacoustic inverse models for the explosive source that need to account for the gas-generated and spall-induced signals and their effect on peak pressures and other acoustic signal features.

58 GEOSCIENCES↗
DOE OSTI2021

Strain wave pathway to semiconductor-to-metal transition revealed by time-resolved X-ray powder diffraction

One of the main challenges in ultrafast material science is to trigger phase transitions with short pulses of light. Here we show how strain waves, launched by electronic and structural precursor phenomena, determine a coherent macroscopic transformation pathway for the semiconducting-to-metal transition in bistable Ti 3 O 5 nanocrystals. Employing femtosecond powder X-ray diffraction, we measure the lattice deformation in the phase transition as a function of time. We monitor the early intra-cell distortion around the light absorbing metal dimer and the long range deformations governed by acoustic waves propagating from the laser-exposed Ti 3 O 5 surface. We developed a simplified elastic model demonstrating that picosecond switching in nanocrystals happens concomitantly with the propagating acoustic wavefront, several decades faster than thermal processes governed by heat diffusion.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Monitoring urban hydrological environment monitoring using fiber optical sensing

The goal of this research is to develop an innovative approach to monitoring urban dynamic hydrological environment using optical-fiber-based distributed acoustic sensing (DAS). We propose to employ surface-wave dispersion inversion and coda-wave interferometry inversion based on natural and anthropogenic seismic ambient noises acquired using a DAS array to accurately invert for time-lapse changes of groundwater in urban areas. Through this work, we aim to provide the first cost-effective, city-scale, high-resolution approach to monitoring urban groundwater changes that is otherwise difficult or expensive to achieve using conventional techniques. This work will create LANL’s DAS-based near-surface and groundwater monitoring capability, which can also be leveraged in solving relevant applied energy and national security problems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

High-Resolution Near-Surface Imaging at the Basin Scale Using Dark Fiber and Distributed Acoustic Sensing: Toward Site Effect Estimation in Urban Environments

Near-surface seismic structure, particularly the shear wave velocity (V s ), can strongly affect local site response, and should be accurately estimated for ground motion prediction during seismic hazard assessment. The Imperial Valley (California), occupying the southern end of the Salton Trough, is a seismically active basin with thick surficial lacustrine sedimentary deposits. In this study, we utilize ambient noise records and local earthquake events for high-resolution near-surface characterization and site effect estimation with an unlit fiber-optic telecommunication infrastructure (dark fiber) in Imperial Valley by using the distributed acoustic sensing (DAS) technique. We apply ambient noise interferometry to retrieve coherent surface waves from DAS records, and evaluate performances of three different surface wave methods on DAS ambient noise dispersion imaging. We develop a quality control workflow to improve the dispersion measurement of noisy portions of the DAS data set by using a data selection strategy. Using the joint inversion of both the fundamental mode and higher overtones of Rayleigh waves, a high resolution two-dimensional (2D) V s structure down to 70 m depth is obtained. We successfully achieve an improved V s 30 (the time-averaged shear-wave velocity in the top 30 m) model with higher spatial-resolution and reliability compared to the existing community model for the area. We also explore the potential for utilizing DAS earthquake events for site amplification estimation. The preliminary results reveal a clear anti-correlation between the approximated site response and the V s 30 profile. In conclusion, our results indicate the potential utility of DAS deployed on dark fiber for near-surface characterization in appropriate contexts.

58 GEOSCIENCES↗
DOE OSTI2023

Sound Assisted Low Temperature Spalling: Upscaling and Throughput

In evaluating the state-of-the-art spalling techniques, the key factor that is missing is the control of crack propagation, which is responsible for low quality wafers due to high surface roughness and Total Thickness Variation, TTV. This project has shown that acoustic waves can alter the conditions at the crack tip during thermal spalling, modifying the propagation speed and the stress field around it (stress intensity factor, K). Both factors are directly related to the roughness of the spalled surface, so as we control speed and stress, we can fine-tune the surface roughness through the entire wafer as suggested by Arakawa.

42 ENGINEERING↗
DOE OSTI2022

Multiscale and multiphysics FEA simulation and materials optimization for laser ultrasound transducers

In this study, the relationship between the nanocomposite design and the laser ultrasound transducer (LUT) characteristics was investigated through simulations in multiple scale levels for material behavior, device response, and acoustic wave propagation in media. First, the effects of the nanoparticle size and concentration on the effective properties of composites were quantitatively investigated with the finite element analysis (FEA) method. Second, the effective properties of the nanocomposite were assigned to the layer, which is modeled as a homogeneous material, in the FEA for the LUT simulating the energy conversion from the incident laser to the acoustic wave. Finally, the ultrasound propagation in the water was calculated by a theoretical wave propagation model. The FEA-based prediction was compared with the experimental data in the literature and a theoretical analysis for LUT based on Thermal-Acoustic coupling. As a result, the ultrasound waves on the transducer surface and at a distance in the water could be predicted. Based on the hierarchically integrated prediction procedure, the optimal conditions of the photoacoustic nanocomposites were investigated through the parametric study with the particle size and concentration as variables. The results guide the material designs optimized for different device characteristics, such as high pressure and broad bandwidth.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Kimberlina 1.2 CCUS Geophysical Models and Synthetic Data Sets

This synthetic multi-scale and multi-physics data set was produced in collaboration with teams at the Lawrence Berkeley National Laboratory, National Energy Technology Laboratory, Los Alamos National Laboratory, and Colorado School of Mines through the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative. Data are associated with the following publication: Alumbaugh, D., Gasperikova, E., Crandall, D., Commer, M., Feng, S., Harbert, W., Li, Y., Lin, Y., and Samarasinghe, S., “The Kimberlina Synthetic Geophysical Model and Data Set for CO2 Monitoring Investigations”, The Geoscience Data Journal, 2023, DOI: 10.1002/gdj3.191. The dataset uses the Kimberlina 1.2 CO2 reservoir flow model simulations based on a hypothetical CO2 storage site in California (Birkholzer et al., 2011; Wainwright et al., 2013). Geophysical properties models (P- and S-wave seismic velocities, saturated density, and electrical resistivity) were produced with an approach similar to that of Yang et al. (2019) and Gasperikova et al. (2022) for 100 Kimberlina 1.2 reservoir models. Links to individual resources are provided below: [CO2 Saturation Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-co2-saturation-models); Resistivity Models – [part 1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-1), [part 2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-2), and [part 3](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-3); [Vp Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vp-velocity-models); [Vs Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vs-velocity-models); [Density Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-density-models). The 3D distributions of geophysical properties for the 33 time stamps of the SIM001 model were used to generate synthetic seismic, gravity, and electromagnetic (EM) responses for 33 times between zero and 200 years. Synthetic surface seismic data were generated using 2D and 3D finite-difference codes that simulate the acoustic wave equation (Moczo et al., 2007). 2D data were simulated for six point-pressure sources along a 2D line with 10 m receiver spacing and a time spacing of 0.0005 s. 3D simulations were completed for 25 surface pressure sources using a source separation of 1 km in both the x and y directions and a time spacing of 0.001 s. Links to individual resources are provided below: [2D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-velocity-models) and [2D surface seismic data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-surface-seismic-data). [3D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-velocity-models), and 3D seismic data [year0](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year0), [year1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year1), [year2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year2), [year5](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year5), [year10](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year10), [year15](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year15), [year20](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year20), [year25](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year25), [year30](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year30), [year35](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year35), [year40](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year40), [year45](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year45), [year49](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year49), [year50](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year50), [year51](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year51), [year52](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year52), [year55](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year55), [year60](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year60), [year65](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year65), [year70](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year70), [year75](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year75), [year80](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year80), [year85](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year85), [year90](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year90), [year95](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year95), [year100](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year100), [year110](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year110), [year120](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year120), [year130](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year130), [year140](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year140), [year150](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year150), [year175](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year175), [year200](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year200). The Python scripts to read these models and data are provided [here](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-python-scripts). EM simulations used a borehole-to-surface survey configuration, with the source located near the reservoir level and receivers on the surface using the code developed by Commer and Newman (2008). Pseudo-2D data for the source at [2500 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz2500m) and [3025 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz3025m), used a 2D inline receiver configuration to simulate a response over 3D resistivity models. The [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-csem-data) contain electric fields generated by borehole sources at monitoring well locations and measured over a surface receiver grid. Vector gravity data, both on the surface and in boreholes, were simulated using a modeling code developed by Rim and Li (2015). The simulation scenarios were parallel to those used for the EM: [pseudo-2D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were calculated along the same lines and within the same boreholes, and [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were simulated over 3D models on the surface and in three monitoring wells. A series of [synthetic well logs](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-well-logs) of CO2 saturation, acoustic velocity, density, and induction resistivity in the injection well and three monitoring wells are also provided at 0, 1, 2, 5, 10, 15, and 20 years after the initiation of injection. These were constructed by combining the low-frequency trend of the geophysical models with the high-frequency variations of actual well logs collected in the Kimberlina 1 well that was drilled at the proposed site. Measurements of permeability and pore connectivity were made on cores of Vedder Sandstone, which forms the primary reservoir unit: [CT micro scans](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-ct-micro-scans-of-vedder-formation) and [Industrial CT Images](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-industrial-ct-images-vedder-formation). These measurements provide the range of scales in the otherwise synthetic data set to be as close to a real-world situation as possible. References: Birkholzer, J.T., Zhou, Q., Cortis, A. and Finsterle, S., 2011. A sensitivity study on regional pressure buildup from large-scale CO2 storage projects. Energy Procedia, 4, 4371-4378. Commer, M., and Newman, G.A., 2008. New advances in three-dimensional controlled-source electromagnetic inversion, Geophysical Journal International, 172, 513-535. Gasperikova, E., Appriou, D., Bonneville, A., Feng, Z., Huang, L., Gao, K., Yang, X., Daley, T., 2022, Sensitivity of geophysical techniques for monitoring secondary CO2 storage plumes, Int. J. Greenh. Gas Control, Volume 114, 103585, ISSN 1750-5836, https://doi.org/10.1016/j.ijggc.2022.103585. Moczo, P., J.O. Robertsson and L. Eisner, 2007, The finite-difference time-domain method for modeling of seismic wave propagation: Advances in geophysics, 48, 421-516. Rim, H., and Y. Li, 2015, Advantages of borehole vector gravity in density imaging, Geophysics, 80, G1-G13. Wainwright, H. M.; Finsterle, S.; Zhou, Q.; Birkholzer, J. T., 2013. Modeling the Performance of Large-Scale CO2 Storage Systems: A Comparison of Different Sensitivity Analysis Methods. International Journal of Greenhouse Gas Control, 17, 189205. https://doi.org/10.1016/j.ijggc.2013.05.007, DOI: 10.18141/1603331. Yang, X., Buscheck, T.A., Mansoor, K., Wang, Z., Gao, K., Huang, L., Appriou, D., and Carroll, S.A., 2019. Assessment of geophysical monitoring methods for detection of brine and CO2 leakage in drinking water aquifers, International Journal of Greenhouse Gas Control, 90, 102803, https://doi.org/10.1016/j.ijggc.2019.102803.

CCUS↗
DOE OSTI2021

Joint Bayesian Inference for Near-Surface Explosion Yield and Height-of-Burst

Forensic capabilities to understand chemical and nuclear explosions are greatly aided by an accurate estimate of explosive yield with uncertainty. The relationship between explosive size and geophysical observations of seismic, acoustic, and optical waves can be exploited to provide an estimate of yield. Any near-surface yield estimate is complicated by the surface interaction, so an estimate for the explosion height-of-burst is necessarily included in the relationship. Additionally, the relationship dictates a trade-off between estimates of yield and height-of-burst. Fortunately, the surface interaction for each type of observation is different, which breaks the trade-off, and the inclusion of height-of-burst with multiple data types improves yield estimation. We define simple parametric forward models to relate seismoäcoustoöptic observations from a data set of known explosive yields and height-of-bursts. The parameters of the models and a prediction for the yield and height-of-burst of a new event can then be estimated given new observations via Bayesian inference. We report posterior distribution estimates of the parametric models using a Markov chain Monte Carlo sampling technique. These models are then used to predict the yield and height-of-burst of SUGAR, a historical near-surface nuclear explosion, using its reported historical observations. The reported yield of 1.2 ktonne Trinitrotoluene (TNT)-equivalent (Department of Energy, 2015) is within the estimated posterior. Yield uncertainty can be estimated from the spread of the posterior, which is between 0.9 and 2.1 ktonne TNT-equivalent. The posterior for height-of-burst has a wider range between 10 m below and 8 m above ground that includes the true height-of-burst of 1 m.

58 GEOSCIENCES↗
DOE OSTI2023

Seismoacoustic Explosion Yield and Depth Estimation: Insights from the Large Surface Explosion Coupling Experiment

The Large Surface Explosion Coupling Experiment (LSECE) is a chemical explosion experiment conducted in Yucca Flat at the Nevada National Security Site in 2020. The experiment included two surface detonations of ~1000 kg trinitrotoluene equivalent. The main goal of this experiment was to provide the ground-truth data for seismoacoustic wave excitation by large chemical explosions near the ground surface. The seismic and acoustic energy partitioning between the surface is strongly governed by the depth or height of explosions, and either seismic or acoustic-only analysis may have inherent ambiguity in determining explosion yield and depth simultaneously. Previous studies suggested that joint seismoacoustic analysis can resolve the trade-off and reduce the uncertainty of yield and depth estimation dramatically. Here we demonstrate the capability of seismoacoustic analysis to improve the accuracy of explosion yield and depth estimation with the LSECE data. Local acoustic wave propagation in the atmosphere can be substantially affected by constantly varying weather conditions. Consisting of two detonations before dawn and in the afternoon, LSECE provides unique data to evaluate the model accuracy of acoustic wave propagation and seismoacoustic energy partitioning depending on local atmospheric conditions. We quantitatively evaluate the accuracy of yield and depth estimation depending on atmospheric variability and the improvement achieved by the joint seismoacoustic approach.

58 GEOSCIENCES↗
DOE OSTI2025

Acoustic waves from the 2023 April 20 SpaceX Starship rocket explosion travelling in the elevated ‘AtmoSOFAR’ channel

SUMMARY The ability to detect low-frequency sounds from distant energetic events depends on the temperature and wind structure of the atmosphere. Thus, from time to time surface-based acoustic detectors may not be able to capture sounds arriving from certain directions. However, the temperature minimum at the tropopause may create an acoustic duct called the ‘AtmoSOFAR’ channel that could transmit acoustic waves laterally—but perhaps not to the ground. If true, then elevated sensors such as those borne aloft by balloons would record the signatures even in regions where ground-based sensors cannot. This has been difficult to prove because high altitude acoustic sources are rare and balloon deployments are sporadic. This work describes the detection and characterization of powerful acoustic waves generated during the launch and terminal explosion of the SpaceX Starship rocket on 2023 April 20 using a pair of microbarometers on a stratospheric balloon. The signals travelled through the AtmoSOFAR channel, carrying information about the size and nature of their source. This channel also appears to leak some acoustic energy to the ground, in agreement with previous studies. The acoustic yield of the Starship terminal explosion was on the order of $10^3$ tons trinitrotoluene equivalent, which agrees with the estimated fuel load of the vehicle to about a factor of 2, but is two orders of magnitude larger than optical estimates. These results support an earlier study that claimed lateral transmission of sound from a smaller rocket through the AtmoSOFAR channel. The transmission of source information through the AtmoSOFAR channel motivates its use for monitoring other natural and anthropogenic events using balloon-borne sensors. This may become increasingly important as more and more private and government entities conduct spacecraft launches and reentries. It may also provide a means of monitoring clear air turbulence and other sound-generating atmospheric phenomena at a distance.

58 GEOSCIENCES↗