Many geologic materials and minerals are seismically anisotropic, with the most general anisotropic material having up to 21 independent elastic coefficients. This report outlines the development of a 3-D, generally anisotropic, linear elastic full waveform finite-difference solver. First, a mathematical description of the solution equations will be described. The finite-difference implementation of these equations will then be shown. Finally, a comparison of results from this new solver to other solutions will be provided as verification that the new algorithm can accurately replicate these solutions. ACKNOWLEDGEMENTS The author also acknowledges the National Nuclear Security Administration, Defense Nuclear Nonproliferation Research and Development (DNN R&D), and the Source Physics Experiment (SPE) working group, a multi-institutional and interdisciplinary group of scientists and engineers.
ParelastiFWl is a python-based frontend to the seismic full waveform inversion process using Sandia Geophysics Department's 3-D isotropic elastic full waveform simulation code, Parelasti. The arguments one provides to ParelastiFWl guide the full waveform inversion process, including resolution of the inversion grid and basic regularization. This report outlines the user flags and ParelastiFWI usage to control the full waveform inversion procedure. ACKNOWLEDGEMENTS The author also acknowledges the National Nuclear Security Administration, Defense Nuclear Nonproliferation Research and Development (DNN R&D), and the Source Physics Experiment (SPE) working group, a multi-institutional and interdisciplinary group of scientists and engineers.
We invert far-field infrasound data for the equivalent seismoacoustic timedomain moment tensor to assess the effects of variable atmospheric models and source phenomena. The infrasound data were produced by a series of underground chemical explosions that were conducted during the Source Physics Experiment (SPE), which was originally designed to study seismoacoustic signal phenomena. The first goal of this work is to investigate the sensitivity of the inversion to the variability of the estimated atmospheric model. The second goal is to determine the relative contribution of two presumed source mechanisms to the observed infrasonic wavefield. Rather than using actual atmospheric observations to estimate the necessary atmospheric Green’s functions, we build a series of atmospheric models that rely on publicly available, regional-scale atmospheric observations. The atmospheric observations are summarized and interpolated onto a 3D grid to produce a model of sound speed at the time of the experiment. For each of four SPE acoustic datasets that we invert, we produced a suite of three atmospheric models for each chemical explosion event, based on 10 yrs of meteorological data: an average model, which averages the atmospheric conditions for 10 yrs prior to each SPE event, as well as two extrema models. To parameterize the inversion, we assume that the source of infrasonic energy results from the linear combination of explosion-induced surface spall and linear seismic-to-elastic mode conversion at the Earth’s free surface. We find that the inversion yields relatively repeatable results for the estimated spall source. Conversely, the estimated isotropic explosion source is highly variable. This suggests that 1) the majority of the observed acoustic energy is produced by the spall and/or 2) our modeling of the elastic energy, and the subsequent conversion to acoustic energy, is too simplistic.
As part of the Source Physics Experiment (SPE) Phase I shallow chemical detonation series, multiple surface and borehole active-source seismic campaigns were executed to perform high- resolution imaging of seismic velocity changes in the granitic substrate. Cross-correlation data processing methods were implemented to efficiently and robustly perform semi-automated change detection of first-arrival times between campaigns. The change detection algorithm updates the arrival times, and consequently the velocity model, of each campaign. The resulting tomographic imagery reveals the evolution of the subsurface velocity structure as the detonations progressed. ACKNOWLEDGEMENTS The authors thank Dan Herold, Bob White, Kale Mc Lin, Ryan Emmit, Maggie Townsend, Curtis Obi, Fred Helsel, Rebekah Lee, Liam Toney, Matt Geuss, and Josh Feldman for their direct and invaluable contributions to this work. Sandia National Laboratories is a multi-mission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525. This paper describes objective technical results and analysis. Any subjective views or opinions that might be expressed in the paper do not necessarily represent the views of the U.S. Department of Energy or the United States Government. Note that a more detailed manuscript for this work is being prepared for publication in the Bulletin of the Seismological Society of America (BSSA).
This document serves to guide a researcher through the process of running the Weather Research and Forecasting (WRF) model and incorporating observations into coarse resolution reanalysis products to model atmospheric conditions at high (50 m) resolution. This documentation is specific to WRF and the WRF Preprocessing System (WPS) version 3.8.1 and the Objective Analysis (OBSGRID) code released on April 8, 2016. Output from WRF serves as an input into the Time-Domain Atmospheric Acoustic Propagation Suite (TDAAPS) which performs staggered-grid finite difference modeling of the acoustic velocity pressure system to produce Green's functions through these atmospheric models.
Full waveform inversion allows the seismologist to utilize an entire waveform and all the information it contains to help image the 3-D structure of the interior of the earth. This report summarizes the basic theory that has been developed in full waveform seismic inversion, primarily related to computation of sensitivity kernels. It then describes the implementation of this theory using Sandia Geophysics Department's Parelasti code, a 3-D full waveform elastic simulation algorithm. Finally, the code is validated using synthetics from simple homogeneous elastic earth models.