Source Physics Experiment Research Articles

Seismic Characterization of Subsurface Structures at Rock Valley, Nevada

ABSTRACT The Source Physics Experiment (SPE) aims at improving explosion monitoring techniques by investigating source characteristics of chemical explosions in geologic formations. One of the critical tasks in Rock Valley Direct Comparison (RV/DC), SPE phase III, is to prepare for the main experiment by characterizing the subsurface structures at the test site. Based on the seismic data acquired during an accelerated-weight-drop (AWD) seismic survey at Rock Valley, we first pick the P-wave first-arrival travel times and derive a P-wave velocity model using the adjoint-state first-arrival travel-time tomography. We then apply reverse-time migration to the processed seismic data and obtain high-resolution images of the subsurface structures along the two main survey lines. Our migration results show several reflectors corresponding to major geologic formation boundaries. We employ a multitask machine learning model to enhance the reverse-time migration images and identify faults from these images. We find that our automatically picked faults correlate well with the locations of known faults in the region in addition to many geologically undetected faults. Our subsurface characterization results refine our understanding of the geology in this region and provide valuable velocity and structural information for RV/DC geologic model building and fault identification.

Spatially Dependent Seismic Wavefield Scattering from an Underground Chemical Explosion: Analysis of the Source Physics Experiment Dry Alluvium Geology Large-N Array

ABSTRACT Explosion sources have been observed to generate significant shear-wave energy despite their isotropic nature. To investigate this phenomenon, we conduct an analysis of the seismic data collected as part of the Source Physics Experiment (SPE): Dry Alluvium Geology (DAG) and investigate the generation of shear-wave energy via scattering. The data were produced by three underground chemical explosions and consist of three-component seismograms, which were recorded by the DAG Large-N array. Synthetic tests suggest that for the DAG experiments, small-scale stochastic heterogeneities, defined as features with correlation lengths of 10–100s of meters, are more effective than large-scale geologic structure (scales >1–10 km) at reproducing the scattering of explosion generated wavefields observed at DAG. We analyze the seismic data for spatially variable ratios between transversely and radially polarized seismic energy, and then estimate the mean free path of P and S waves. All analyses are conducted within a frequency band of 5–50 Hz. The ratio of transversely to radially polarized energy is the highest in the east and west portion of the Large-N array. In addition, the magnitude of the estimated S-wave mean free path is shorter in the eastern portion of the Large-N array. This variation indicates that the eastern area of the DAG array is where more scattering is occurring, suggesting azimuthal dependence of P-to-P and P-to-S scattering. This azimuthal dependence of P-to-S scattering can have implications for explosion discrimination based on spectral ratios of seismic wave types, because the general assumption is that explosions do not generate shear-wave energy. Synthetic tests modeling only larger-scale geologic structure had lower transversely polarized energy (only four stations showing a transversely to radially polarized energy ratio greater than 1) and fewer stations (<10) displaying shorter (<300 m) mean free paths than what was observed in the DAG data results.

How Good Is Your Location? Comparing and Understanding the Uncertainties in Location for the 1993 Rock Valley Sequence

AbstractAccurate event locations are important for many endeavors in seismology, and understanding the factors that contribute to uncertainties in those locations is complex. In this article, we present a case study that takes an in-depth look at the accuracy and precision possible for locating nine shallow earthquakes in the Rock Valley fault zone in southern Nevada. These events are targeted by the Rock Valley Direct Comparison phase of the Source Physics Experiment, as candidates for the colocation of a chemical explosion with an earthquake hypocenter to directly compare earthquake and explosion sources. For this comparison, it is necessary to determine earthquake hypocenters as accurately as possible so that different source types have nearly identical locations. Our investigations include uncertainty analysis from different sets of phase arrivals, stations, velocity models, and location algorithms. For a common set of phase arrivals and stations, we find that epicentral locations from different combinations of velocity models and algorithms are within 600 m of one another in most cases. Event depths exhibit greater uncertainties, but focusing on the S-P times at the nearest station allows for estimates within approximately 500 m.

Time-Varying Moment-Tensor Analysis with Application to Buried Chemical Explosions

Abstract The Source Physics Experiment (SPE) Phase I consisted of a series of over-buried, horizontally colocated chemical explosions at the Nevada National Security Site. Seismic waveforms from these explosions recorded at near-source accelerometers, local geophone arrays, and regional seismic stations provided a rich suite of observations suitable for resolving fine source details. To investigate the time-varying source history of the explosions, we used the frequency-domain moment-tensor inversion method described in Yang et al. (2018) with added regularization and reconstruction to suppress the nonuniqueness evident in unconstrained inversion results. The inverted moment-rate spectra are accurate within the response band of the local geophones and, in all cases, display predominately isotropic characteristics. For SPE-4Prime, SPE-5, and SPE-6, we resolve predominately isotropic moment release followed by double couple and compensated linear vector dipole (CLVD) release later in the time-varying source history. We interpret these results both in terms of absolute depth and scaled depth of burial. The apparent non-isotropic release from SPE-4Prime and SPE-5 may simply reflect increased resolving power related to improved Earth model accuracy at greater absolute depths, whereas the non-isotropic release from SPE-6 likely reflects the larger damage associated with an event at a shallower scaled depth. These results provide insight into the time-varying source characteristics of shallow explosions and motivation to study shear-wave generation by inverting for fracture, spallation, induced slip, and other temporally delayed source processes through time-varying methods.

Global Nuclear Explosion Discrimination Using a Convolutional Neural Network

AbstractUsing P‐wave seismograms, we trained a seismic source classifier using a Convolutional Neural Network. We trained for three classes: earthquake P‐wave, underground nuclear explosion (UNE) P‐wave, and noise. With the current absence of nuclear testing by countries that have signed the Comprehensive Test Ban Treaty, high quality seismic data from UNEs is limited. Even with limited training data, our model can accurately characterize most events recorded at regional and teleseismic distances, finding over 95% signals in the validation set. We applied the model on holdout datasets of the North Korean test explosions to evaluate the performance on unique region and station‐source pairs, with promising results. Additionally, we tested on the Source Physics Experiment events to investigate the potential for chemical explosions to act as a surrogate for nuclear explosions. We anticipate that machine‐learning models like our classifier system can have broad application for other seismic signals including volcanic and non‐volcanic tremor, anomalous earthquakes, ice‐quakes or landslide‐quakes.

Seismoacoustic Signatures Observed During a Long-Term Deployment of Infrasound Sensors at the Nevada National Security Site

ABSTRACT Earthquakes have repeatedly been shown to produce inaudible acoustic signals (<20 Hz), otherwise known as infrasound. These signals can propagate hundreds to thousands of kilometers and still be detected by ground-based infrasound arrays depending on the source strength, distance between source and receiver, and atmospheric conditions. Another type of signal arrival at infrasound arrays is the seismic induced motion of the sensor itself, or ground-motion-induced sensor noise. Measured acoustic and seismic waves produced by earthquakes can provide insight into properties of the earthquake such as magnitude, depth, and focal mechanism, as well as information about the local lithology and atmospheric conditions. Large earthquakes that produce strong acoustic signals detected at distances greater than 100 km are the most commonly studied; however, more recent studies have found that smaller magnitude earthquakes (Mw<2.0) can be detected at short ranges. In that vein, this study will investigate the ability for a long-term deployment of infrasound sensors (deployed as part of the Source Physics Experiments [SPE] from 2014 to 2020) to detect both seismic and infrasonic signals from earthquakes at local ranges (<50 km). Methods used include a combination of spectral analysis and automated array processing, supported by U.S. Geological Survey earthquake bulletins. This investigation revealed no clear acoustic detections for short range earthquakes. However, secondary infrasound from an Mw 7.1 earthquake over 200 km away was detected. Important insights were also made regarding the performance of the SPE networks including detections of other acoustic sources such as bolides and rocket launches. Finally, evaluation of the infrasound arrays is performed to provide insight into optimal deployments for targeting earthquake infrasound.

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

Abstract 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 ds (depth d scaled by the cube-root of explosive yield w1/3) as exp(−ds) but then transitions near a scaled depth of 6 m/ton1/3 to the spall-generated mechanism in which the decay is ds−7/4. This decay form is related to the strong ground-motion attenuation relationship that affects spall strength. 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.

Finite-Difference Simulation for Infrasound Generated by Finite-Extent Ground Motions

Abstract Underground explosions can produce infrasound in the atmosphere, and the wavefield characteristics are often governed by the ground surface motions. Finite-difference methods are popular for infrasound simulation as their generality and robustness allow for complex atmospheric structures and surface topography. A simple point-source approximation is often used because infrasound wavelengths tend to be large relative to the source dimensions. However, this assumption may not be able to capture the complexity of explosion-induced ground motions if the surface area is not compact, and appropriate source models must be incorporated into the finite-difference simulations for accurate infrasound prediction. In this study, we develop a point source representation of the complex ground motions for infrasound sources. Instead of a single point source, we use a series of point sources distributed over the source area. These distributed point sources can be equivalent to air volume changes produced by the ground motions in the atmosphere. We apply the distributed point-source method to a series of buried chemical explosions conducted during the Source Physics Experiment Phase I. Epicentral ground-motion measurements during the experiments provide a way to calculate accurate distributed point sources. We validate and evaluate the accuracy of distributed point source approach for infrasound simulations by direct comparison with acoustic observations in the field experiment.

Inversion of Infrasound Time Series for Seismoacoustic Source Parameters Produced by a Buried Chemical Explosion at the Source Physics Experiment Phase II: Dry Alluvium Geology

ABSTRACT We invert infrasonic timeseries produced by a single buried chemical explosion detonated as part of the Source Physics Experiment phase II Dry Alluvium Geology (DAG). The key to our inversion is that we use a three dimensional, fully coupled, linear, elastic-to-acoustic forward model to predict the recorded acoustic wavefield. We assume that the fourth buried chemical explosion in the series, DAG-4, produces an atmospheric acoustic wavefield due to a combination of direct elastic-to-acoustic coupling at the air–earth interface and ground upheaval, or spall, at the Earth’s surface. In our linear model we approximate these two phenomena with distinct source terms separated in space and time: (1) a buried point source consisting of six source time functions, each corresponding to a single component of a generalized moment tensor and (2) a vertically directed, time-variable force applied to the Earth’s surface. Inverting the infrasonic data with this linear forward model results in estimated source parameters that accurately predict the observed infrasound signals. To test our estimated buried source, we use the Rayleigh integral (RI) to model the infrasound signal from the measured and predicted acceleration surfaces. First, we simulate the infrasound signal waveforms using vertical accelerometer data and find that these closely match the observed infrasound. Next, using the estimated buried seismic source model without spall, we estimate the acceleration of the Earth’s surface using a linear approximation. When using the predicted surface acceleration in the RI, we simulate infrasound signals that reproduce pulse shape characteristics but not the amplitude of the observed infrasound, indicating our linear model does not completely account for nonlinear spall effects. Based on these results, we argue that for the scale of this experiment, purely linear models can estimate seismoacoustic source attributes and predict the far-field infrasonic signal, but the primary contribution to the infrasonic signal is the spall source.

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

AbstractThe 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. 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.

Source Separation and Medium Change of Contained Chemical Explosions from Coda Wave Interferometry

Abstract Differences in the seismic coda of neighboring events can be used to investigate source location offsets and medium change with coda wave interferometry (CWI). We employ CWI to infer the known relative location between two chemical explosions in Phase I of the Source Physics Experiment (SPE). The inferred displacement between the first, SPE-1, and second, SPE-2, chemical explosion is between 6 and 18 m, with an expectation of 9.2 m, where the known separation is close to 9.4 m. We also employ CWI to find any velocity perturbation due to damage from SPE-2, by comparing its coda with the collocated third SPE chemical explosion, SPE-3. We find that damage due to SPE-2 must be confined to a spherical region with radius less than 10 m and velocity perturbation less than 25%.

Investigating the detection of underground nuclear explosions by radon displacement

A novel approach is proposed to detect underground nuclear explosions (UNEs) through the displacement of natural radon isotopes (222Rn and 220Rn). Following an explosion, it is hypothesized that the disturbance and pressurization of the sub-surface would facilitate the movement of radon from the depth of the UNE towards the surface resulting in increased soil gas activity. The resulting signal may be magnified by a factor of 2.0–4.9 by the decay of radon to its short-lived progeny. Increases in background activity may be useful for identifying locations to perform additional measurements, or as a detectable signal at monitoring stations. To validate this hypothesis, radon detection instrumentation was deployed at the Dry Alluvium Geology (DAG) site of the Source Physics Experiment (SPE) at the Nevada National Security Site (NNSS). Natural fluctuations in the soil gas activity due to barometric pumping, and the lower yield of the chemical explosions (1–50 t) made it difficult to confirm a displacement of radon from the explosions, and further study to validate the proposed hypothesis is recommended.

Using Dense Array Waveform Correlations to Build a Velocity Model with Stochastic Variability

ABSTRACTIn an ongoing effort to improve 3D seismic-wave propagation modeling for frequencies up to 10 Hz, we used cross correlations between vertical-component waveforms from an underground chemical explosion to estimate the statistical properties of small-scale velocity heterogeneities. The waveforms were recorded by a dense 2D seismic array deployed during the Source Physics Experiments for event number 5 (SPE-5) in a series of six underground chemical explosions, conducted at the Nevada National Security Site. The array consisted of 996 geophones with a 50–100 m grid spacing, deployed at the SPE site at the north end of the Yucca Flat basin. The SPE were conducted to investigate the generation and propagation of seismic and acoustic waves from underground explosions.Comparisons of decay rates of waveform cross correlations as function of interstation distance, computed for observed and synthetic seismograms from the SPE-5 chemical explosion, were used to constrain statistical properties of correlated stochastic velocity perturbations representing small-scale heterogeneities added to a geology-based velocity model of the Yucca Flat basin. Using comparisons between recorded and simulated waveform cross correlations, we were able to recover sets of statistical properties of small-scale velocity perturbations in the velocity model that produce the best-fit between the recorded and simulated ground motion. The stochastic velocity fluctuations in the velocity model that produced the smallest misfits have a horizontal correlation distance of between 400 and 800 m, a vertical correlation distance between 100 and 200 m, and a standard deviation of 10% from the nominal model velocity in the alluvium basin layers. They also have a horizontal correlation distance of 1000 m, a vertical correlation distance of 250 m, and a standard deviation of 6% in the underlying and consolidated sedimentary layers, up to a depth of 4 km.Comparisons between observed and simulated wavefields were used to assess the proposed small-scale heterogeneity enhancements to the Yucca Flat basin model. We found that adding a depth-resolved stochastic variability to the geology-based velocity model improves the overall performance of ground-motion simulations of an SPE-5 explosion in the modeled frequency range up to 10 Hz. The results may be applicable to other similar basins.

Comparison of Surface Phenomena Created by Underground Chemical Explosions in Dry Alluvium and Granite Geology From Fully Polarimetric VideoSAR Data

Phase I of the Source Physics Experiment (SPE) series involved six underground chemical explosions, all of which were conducted at the same experimental pad. Research from the sixth explosion of the series (SPE-6) demonstrated that polarimetric synthetic aperture radar (PolSAR) is a viable technology for monitoring an underground chemical explosion when the geologic structure is Cretaceous granitic intrusive. It was shown that a durable signal is measurable by the H/A/α polarimetric decomposition parameters. After the SPE-6 experiment, the SPE program moved to the Phase II location, which is composed of dry alluvium geology (DAG). The loss of wavefront energy is greater through dry alluvium than through granite. In this article, we compare the SPE-6 analysis to the second DAG (DAG-2) experiment. We hypothesize that despite the geology at the DAG site being more challenging than at the Phase I location, combined with the DAG-2 experiment having a 3.37 times deeper scaled depth of burial than the SPE-6, a durable nonprompt signal is still measurable by a PolSAR sensor. We compare the PolSAR time-series measures from videoSAR frames, from the SPE-6 and DAG-2 experiments, with accelerometer data. We show which PolSAR measures are invariant to the two types of geology and which are geology dependent. We compare a coherent change detection (CCD) map from the DAG-2 experiment with the data from a fiber-optic distributed acoustic sensor to show the connection between the spatial extent of coherence loss in CCD maps and spallation caused by the explosion. Finally, we also analyze the spatial extent of the PolSAR measures from both explosions.

Effect of Random 3D Correlated Velocity Perturbations on Numerical Modeling of Ground Motion from the Source Physics Experiment

ABSTRACTExplosions are traditionally discriminated from earthquakes, using the relative amplitude of compressional and shear waves at regional and teleseismic distances known as the P/S discriminant. Pyle and Walter (2019) showed this technique to be less robust at shorter distances, in detecting small-magnitude earthquakes and low-yield explosions. The disparity is largely due to ground motion from small, shallow sources being significantly impacted by near-surface structural complexities. To understand the implications of wave propagation effects in generation of shear motion and P/S ratio during underground chemical explosions, we performed simulations of the Source Physics Experiment (SPE) chemical explosions using 1D and 3D velocity models of the Yucca Flat basin. All simulations used isotropic point sources in the frequency range 0–5 Hz. We isolate the effect of large-scale geological structure and small-scale variability at shallow depth (<5 km), using a regional 3D geologic framework model (GFM) and the GFM-R model derived from the GFM, by adding correlated stochastic velocity perturbations. A parametric study of effects of small-scale velocity variations on wave propagation, computed using a reference 1D velocity model with stochastic perturbations, shows that the correlation length and depth of stochastic perturbations significantly impact wave scattering, near-surface wave conversions, and shear-wave generation. Comparisons of recorded and simulated waveforms for the SPE-5 explosion, using 3D velocity models, demonstrate that the shallow structure of the Yucca Flat basin contributes to generation of observed shear motion. The inclusion of 3D wave scattering, simulated by small-scale velocity perturbations in the 3D model, improves the fit between the simulated and recorded waveforms. In addition, a relatively low intrinsic attenuation, combined with small-scale velocity variations in our models, can confirm the observed wave trapping and its effect on duration of coda waves and the spatial variation of P/S ratio at basin sites.

Comparison of PolSAR Surface Measurements From Underground Chemical Explosions to Recorded and Predicted Surface Ground Motion

The Source Physics Experiment (SPE) Phase I conducted six underground chemical explosions at the same experimental pad with the goal of characterizing underground explosions to enhance the United States (U.S.) ability to detect and discriminate underground nuclear explosions (UNEs). A fully polarimetric synthetic aperture RADAR (PolSAR) collected imagery in VideoSAR mode during the fifth and sixth explosions in the series (SPE-5 and SPE-6). Previously, we reported the prompt PolSAR surface changes cause by SPE-5 and SPE-6 explosions within seconds or minutes of the underground chemical explosions, including a drop of spatial coherence and polarimetric scattering changes. Therein it was hypothesized that surface changes occurred when surface particles experienced upward acceleration greater than 1 g. Because the SPE site was instrumented with surface accelerometers, we explore that hypothesis and report our findings in this article. We equate explosion-caused prompt surface expressions measured by PolSAR to the prompt surface movement measured by accelerometers. We tie these findings to UNE detection by comparing the PolSAR and accelerometer results to empirical ground motion predictions derived from accelerometer recordings of UNEs collected prior to cessation of U.S. nuclear testing. We find the single threshold greater than 1 g hypothesis is not correct for it does not explain the PolSAR results. Our findings show PolSAR surface coherence spatial extent is highly correlated with surface velocity, both measured and predicted, and the resulting surface deformation extent is corroborated by accelerometer records and the predicted lateral spall extent. PolSAR scattering changes measured during SPE-6 are created by the prompt surface displacement being larger than the spall gap.

Seismoacoustic analysis of underground explosions using the Rayleigh integral

SUMMARYA coupled seismoacoustic model is developed for the analysis of acoustic signals produced by underground explosive events with an aim to develop a means of improving estimated depth and yield for explosion monitoring. A ground spall model is used to predict surface motion characteristics produced by an underground explosion and the Rayleigh integral is applied to relate the surface motion to the acoustic signal some distance from surface ground zero. The low-frequency component of the ground motion associated with the prolonged free fall of lofted material during spall is found to dominate the acoustic signal propagating away from surface ground zero at shallow angles. The model is applied to observed ground motion and acoustic signals recorded during the Source Physics Experiment (SPE) with promising results. In addition to accurately predicting characteristics of the observed acoustic signals from several of the SPE events, the model provides a means of explaining the lack of signals observed during several events in the SPE due to the directionality of the higher frequency acoustic signals associated with the uplift and closure components as well as the lack of a spatially localized, longer duration dwell.

The Effects of Atmospheric Models on the Estimation of Infrasonic Source Functions at the Source Physics Experiment

ABSTRACTWe invert infrasound signals for an equivalent seismoacoustic source function using different atmospheric models to produce the necessary Green’s functions. The infrasound signals were produced by a series of underground chemical explosions as part of the Source Physics Experiment (SPE). In a previous study, we inverted the infrasound data using so-called predictive atmospheric models, which were based on historic, regional-scaled, publicly available weather observations interpolated onto a 3D grid. For the work presented here, we invert the same infrasound data, but using atmospheric models based on weather data collected in a time window that includes the approximate time of the explosion experiments, which we term postdictive models. We build two versions of the postdictive models for each SPE event: one that is based solely on the regional scaled observations, and one that is based on both regional scaled observations combined with on-site observations obtained by a weather sonde released at the time of the SPE. We then invert the observed data set three times, once for each atmospheric model type. We find that the estimated seismoacoustic source functions are relatively similar in waveform shape regardless of which atmospheric model that we used to construct the Green’s functions. However, we find that the amplitude of the estimated source functions is systematically dependent on the atmospheric model type: using the predictive atmospheric models to invert the data generally yields estimated source functions that are larger in amplitude than those estimated using the postdictive models.

Characterization of Spall in Hard Rock from Observations and Simulations of the Source Physics Experiment Phase I

ABSTRACTSpall signals from the Source Physics Experiments are presented, analyzed, and modeled for insight to the explosion source. The observed signal is similar in nature to nearby historical nuclear explosions, and the surface force-time history or velocity can be interpreted with the same model. We use the models for peak spall velocity, spalled mass, and spall depth and radius derived from historical nuclear explosions to parameterize the physical force-time history model from Stump (1985) and show that this parameterized model can be used for spall prediction. The spall signal is also investigated with a numerical continuum model that incorporates gravity. Peak velocity and dwell time are well predicted, and the multiple slap-down phases are captured if one includes a weak near-surface layer similar to the geologic observation.

Seismic Characterization of the Nevada National Security Site Using Joint Body Wave, Surface Wave, and Gravity Inversion

ABSTRACTAs a part of the series of Source Physics Experiments (SPE) conducted on the Nevada National Security Site in southern Nevada, we have developed a local-to-regional scale seismic velocity model of the site and surrounding area. Accurate earth models are critical for modeling sources like the SPE to investigate the role of earth structure on the propagation and scattering of seismic waves. We combine seismic body waves, surface waves, and gravity data in a joint inversion procedure to solve for the optimal 3D seismic compressional and shear-wave velocity structures and earthquake locations subject to model smoothness constraints. Earthquakes, which are relocated as part of the inversion, provide P- and S-body-wave absolute and differential travel times. Active source experiments in the region augment this dataset with P-body-wave absolute times and surface-wave dispersion data. Dense ground-based gravity observations and surface-wave dispersion derived from ambient noise in the region fill in many areas where body-wave data are sparse. In general, the top 1–2 km of the surface is relatively poorly sampled by the body waves alone. However, the addition of gravity and surface waves to the body-wave dataset greatly enhances structural resolvability in the near surface. We discuss the methodology we developed for simultaneous inversion of these disparate data types and briefly describe results of the inversion in the context of previous work in the region.