SEARCH · Engineering Papers
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.
This project is in support of the Source Physics Experiment SPE (Snelson et al. 2013), which aims to develop new seismic source models of explosions. One priority of this program is first principle numerical modeling to validate and extend current empirical models.
The Source Physics Experiments (SPE) were designed to improve our physics-based understanding of explosion sources for the purposes of nuclear test monitoring. Phase I consisted of 6 chemical explosions in the Climax Stock Granite of the Nevada National Security Site (NNSS), while Phase II consisted of 4 explosions in a contrasting dry alluvium geology (DAG) in Yucca Flat, providing essential data in various media and emplacement conditions to further modeling efforts. For Phase III, the Rock Valley Direct Comparison (RVDC) seeks to directly compare earthquake and explosion source types. An unusually shallow series of events in 1993 along the Rock Valley Fault Zone in the southeastern portion of the NNSS has been targeted for this direct comparison. Depth ranges for the events, previously estimated to be less than 3 km, is achievable by modern drilling techniques and accessibility to the epicentral locations would require minimal improvements to the infrastructure. The events providing this unique opportunity for direct comparison are the focus of this report.
The Large Surface Explosion Coupling Experiment (LSECE) project will include two surface energetic events (chemical explosions) conducted as close as possible to ground zero (GZ) of the Source Physics Experiment (SPE) Phase II Dry Alluvium (DAG) site in Yucca Flat at the Nevada National Security Site (NNSS). Seismic and acoustic data from these events will complement DAG buried source data and improve interpretation of speed-of-sound signals from explosions for national security programs. A satellite image of the region near GZ along with some of the recording stations that we will use for the far- eld simulation is shown in Figure 1. The working design of the HE package is shown in the Figure 2 and is roughly a cube of C-4 bricks.
Acoustic signals recorded for explosive events are sensitive to pressure changes in the source region and can provide valuable information about source characteristics. However, acoustic wave propagation in the atmosphere is significantly affected by weather conditions (e.g., temperature, pressure, and wind), and the atmospheric propagation effects should be considered properly to obtain accurate source information. Acoustic waveform inversion method has been increasing useful to determine quantitative source parameters by incorporating accurate Green’s function into the inversion and by effectively removing propagation path effects. This paper reviews recent progress in the acoustic waveform inversion and its application to explosion yield/emplacement (e.g., height-of-burst or depth-of-burial) estimation. The inversion method is applied to acoustic signals recorded during Source Physics Experiment Phase I and Forensics Surface Event conducted at the Nevada National Security Site (NNSS). Inverted acoustic source strengths are related to the explosion yield and emplacement, and their estimation uncertainties are explored by the comparison with the ground-truth data.
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.
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.
The Rock Valley Direct Comparison (RV/DC) is the third phase of the Source Physics Experiment (SPE; Snelson et al., 2014), a project aimed at improving seismic discrimination between explosive sources and natural earthquakes. Phases I and II of SPE focused on studying the generation of shear waves in endmember geologies (granite and alluvium, respectively), while Phase III aims to detonate a chemical explosion co-located with the source region of the 1993 Rock Valley earthquake sequence, located in the Rock Valley Fault Zone (RVFZ), in the southern portion of the Nevada Test Site (now the Nevada National Security Site, or NNSS).
Here, we develop high–resolution, laterally varying attenuation models for the regional crustal phases of Pg and Lg in the area surrounding the Basin and Range Province in the western United States. The models are part of the characterization effort for the Source Physics Experiment (SPE), a series of chemical explosions at the Nevada National Security Site designed to improve our understanding of explosion source phenomenology. To aid in SPE modeling efforts, we focus on improving our ability to accurately predict amplitudes in a set of narrow frequency bands ranging from 0.5 to 16.0 Hz. To explore constraints at higher frequencies where data become more sparse, we test the robustness of the empirically observed power–law relationship between quality factor Q and frequency (Q=Q 0 f γ ). Our methodology uses a staged approach to consider attenuation, physics–based source terms, site terms, and geometrical spreading contributions to amplitude measurements. Tomographic inversion results indicate that the frequency dependence is a reasonable assumption as attenuation varies laterally for this region through all frequency bands considered. Our 2D Pg and Lg attenuation models correlate with underlying physiographic provinces, with the highest Q located in the Sierra Nevada Mountains and the Colorado plateau. Compared to a best–fitting 1D model for the region, the 2D model provides an 81% variance reduction overall for Lg residuals and a 75% reduction for Pg . These detailed attenuation maps at high frequencies will facilitate further study of local and regional distance P/S amplitude discriminants that are typically used to distinguish between earthquakes and underground explosions.
One challenging task in explosion seismology is the development of physical models for explaining the generation of S waves during underground explosions. Recent analysis of ground motion from chemical explosions during the Source Physics Experiment (Pitarka et al., 2015) suggests that, although a large component of shear motion was generated directly at the source, additional scattering from heterogeneous velocity structure and topography is necessary to better match the recorded data. In our paper, we used a stochastic representation of small-scale velocity variability to produce high-frequency scattering and to analyze its implication on shear-motion generation during underground explosions. In our stochastic velocity model, the key parameters that affect scattering are the correlation length and the relative amplitude of velocity perturbations. Finally, based on finite-difference simulations of elastic wave propagation from an isotropic explosion source, we find that higher velocity perturbations result in larger shear motion, whereas the correlation length, which controls the scatterers size, affects the frequency range at which relative transverse motion is larger.
In Rock Valley, Nevada, the USA, past earthquake sequences and shallow faulting in Paleozoic carbonates remain poorly understood. Carbonates are typically more ductile rocks that do not experience large stress drops or fracture coalescence, which contradicts previous observations in the region. As part of the Source Physics Experiment, an effort has been made to experimentally characterize the petrophysical and geomechanical properties of carbonates from Rock Valley. Well core and outcrop samples were used to determine the difference between shallow-buried Tertiary limestones and deeply buried Paleozoic limestones and dolostones. The carbonates possessed porosities between 3 and 9%, with V P and V S in the Tertiary carbonates around 2700 and 1800 m/s, respectively, and 6000 and 3100 m/s, respectively, in the Paleozoic carbonates. Microstructural characterization revealed that the Paleozoic carbonates contained significantly more pre-existing damage and alteration than the Tertiary carbonates, though the deformation varies from localized to diffuse. Unconfined Brazilian tests and triaxial tests with confining pressures between 0 and 50 MPa showed that, although samples all experienced failure, the dolostones typically failed at greater stresses and possessed greater E/ν ratios than the limestones. Furthermore, the Hoek–Brown failure criterion was used to construct failure envelopes with the compressive and tensile failure tests. Velocity, density, and porosity measurements were used to construct a hypothetical velocity-depth profile and compare the results with the theoretical petrophysical measurements. Using brittleness analysis, the likely failure conditions were determined to be low-porosity dolomitic rock for fault nucleation in the Paleozoic basement.
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. Further, 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.
Accurate 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.
We visually examined surface seismic data from the Source Physics Experiment (SPE) Phase-1 chemical explosions in preparation for waveform analysis and interpretation. Recording stations included far-field (near source) geophone lines, and both permanent network, reoccupied USArray, and temporary stations at local and regional distances. We found that near-source waveforms are dominated by broad band compressional (P), and band limited surface wave (Rg) (< 10 Hz) phases. Distinct, consistent shear (S) phases are difficult to observe visually even with a comprehensive bank of frequency bands. Secondary phases are indistinct at near-local distances, unlike for earthquake sources, especially at the high frequencies of interest for discrimination (e.g. 4-8 Hz), likely due to intense near surface scattering. S phases begin to emerge from the scattered background at and beyond 100 km. We identified regional Pn, Pg and Lg phases to 446 km for SPE-5. The data set is, unfortunately, plagued by poorly recorded instrument responses, or, possibly, time varying site amplification due to human interaction with sensors. Instrument response issues should be fixable, and this should be attempted to make this a truly, high-quality data set for use by the monitoring community.
In this, reduced displacement potentials (RDPs) for chemical explosions of the Source Physics Experiments (SPE) in granite at the Nevada Nuclear Security Site are estimated from free-field ground motion recordings. Far-field P wave source functions are proportional to the time derivative of RDPs. Frequency domain comparisons between measured source functions and model predictions show that high-frequency amplitudes roll off as ω -2 , but models fail to predict the observed seismic moment, corner frequency, and spectral overshoot. All three features are fit satisfactorily for the SPE-2 test after cavity radius R c is reduced by 12%, elastic radius is reduced by 58%, and peak-to-static pressure ratio on the elastic radius is increased by 100%, all with respect to the Mueller-Murphy model modified with the Denny-Johnson R c scaling law. A large discrepancy is found between the cavity volume inferred from RDPs and the volume estimated from laser scans of the emplacement hole. The measurements imply a scaled R c of ~5 m/kt 1/3 , more than a factor of 2 smaller than nuclear explosions. Less than 25% of the seismic moment can be attributed to cavity formation. A breakdown of the incompressibility assumption due to shear dilatancy of the source medium around the cavity is the likely explanation. New formulas are developed for volume changes due to medium bulking (or compaction). A 0.04% decrease of average density inside the elastic radius accounts for the missing volumetric moment. Assuming incompressibility, established R c scaling laws predicted the moment reasonable well, but it was only fortuitous because dilation of the source medium compensated for the small cavity volume.
Direct Shear (DS) and Triaxial Shear (FCT) tests from Core holes U-15n and U-15n#10 are part of a larger material characterization effort for the Source Physics Experiment (SPE) project. This larger effort encompasses characterizing a granite body from Nevada both before and after each SPE shot. Core hole U-15n is the vertically oriented source hole for all SPE shots; pre shot core was taken from this hole for DS and FCT testing. After two SPE shots were executed, an inclined core hole (U-15n#10) was drilled; both DS and FCT tests were conducted from this core hole. The first shot (SPE-1) conducted on May 3, 2011 was a calibration shot. SPE-1 was an order of magnitude smaller than the second shot (SPE-2). After SPE-2 was conducted on October 25, 2011 the aforementioned inclined core hole (U-15n#10) was drilled. At its bottom, the inclined core hole intersects the source hole. The third shot (SPE-3) occurred on July 24, 2012. Vertical and inclined core holes were drilled post SPE-3 and specimens will soon be selected for geomechanical characterization. At the time of this writing, work is ongoing at Nevada in preparation for the fourth SPE shot (SPE-4).
As 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. Here, 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.
Space-based and airplane-based synthetic aperture RADAR (SAR) can monitor ground height using interferometric SAR (InSAR) collections. However, fielding the airplane-based SAR is expensive and coordinating the frequency and timing of ground experiments with space-based SAR is challenging. This research explored the possibility of using a small, mobile unmanned aerial vehicle- base (UAV) SAR to see if it could provide a quick and inexpensive InSAR option for the Source Physics Experiment (SPE) Phase III project. Firstly, a local feasibility collection using a UAV-based SAR showed that InSAR products and height measurements were possible, but that in-scene fiducials were needed to assist in digital elevation model (DEM) construction. Secondly, an InSAR collection was planned and executed over the SPE Phase III site using the same platform configuration. We found that the image formation by the SAR manufacturer creates discontinuities, and that noise impacted the generation and accuracy of height maps. These processing artifacts need to be overcome to generate an accurate height map.
The DTRA sponsored 2020-2022 Large Surface Explosion Coupling Experiment (LSECE) consists of two large ground surface chemical explosions, data collection, analysis and modeling carried out in 2020-2022. The LSECE chemical explosions were carried out at the site of prior NNSA sponsored buried chemical explosions that were part of the Source Physics Experiment (SPE) Phase II in Dry Alluvium Geology (DAG). This report described the data collected under LSECE that is being made publicly available. The prior buried explosion DAG data are described in a separate report (https://www.osti.gov/biblio/1825534) and are also publicly available. The LSECE explosions, data analysis and numerical modeling work sponsored by DTRA were intended to address three objectives: 1) Generate seismo-acoustic data to test and improve numerical models of explosion energy coupling in dry alluvium geology; 2) To improve seismo-acoustic yield estimation techniques as a function of depth and medium properties; 3) To study and improve acoustic propagation modeling under two different atmospheric conditions. The LSECE data that were collected have expanded the prior DAG buried explosion dataset to include surface chemical explosions recorded on a common set of seismo-acoustic stations. They have allowed a more detailed study of above and below surface seismo-acoustic energy coupling with applications to explosion monitoring and assessment. The two approximately 1-ton TNT equivalent yield LSECE explosions were detonated at different times of the day to explore the effects of the different atmospheric conditions. The first chemical explosion “Artemis” was conducted before dawn when temperature inversions were present. The second chemical explosion “Apollo” was conducted on a sunny afternoon when the temperature gradient was more linear. The LSECE chemical explosion data were also collected across a variety of different sensor types to allow evaluation of their effectiveness in recovering useful information. The LSECE instrumentation included fiber optic or Distributed Acoustic Sensing (DAS), a dense array or Large-N array of seismometers, borehole accelerometers, DAS and a velocity meter, airborne acoustic instruments and a variety of visual and remote sensing data.