Search NASA⌕ Search

SEARCH · Search NASA

Results for “UNDERGROUND EXPLOSIONS”

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 91 records · Page 5

Dense Seismic Array Study of a Legacy Underground Nuclear Test at the Nevada National Security Site

The complex postdetonation geologic structures that form after an underground nuclear explosion are difficult to constrain because increased heterogeneity around the damage zone affects seismic waves that propagate through the explosion site. Generally, a vertical rubble-filled structure known as a chimney is formed after an underground nuclear explosion that is composed of debris that falls into the subsurface cavity generated by the explosion. Compared with chimneys that collapse fully, leaving a surface crater, partially collapsed chimneys can have remnant subsurface cavities left in place above collapsed rubble. The 1964 nuclear test HADDOCK, conducted at the Nevada test site (now the Nevada National Security Site), formed a partially collapsed chimney with no surface crater. Understanding the subsurface structure of these features has significant national security applications, such as aiding the study of suspected underground nuclear explosions under a treaty verification. In this study, we investigated the subsurface architecture of the HADDOCK legacy nuclear test using hybrid 2D–3D active source seismic reflection and refraction data. The seismic data were acquired using 275 survey shots from the Seismic Hammer (a 13,000 kg weight drop) and 65 survey shots from a smaller accelerated weight drop, both recorded by ~1000 three-component 5 Hz geophones. First-arrival, P-wave tomographic modeling shows a low-velocity anomaly at ~200 m depth, likely an air-filled cavity caused by partial collapse of the rock column into the temporary postdetonation cavity. A high-velocity anomaly between 20 and 60 m depth represents spall-related compaction of the shallow alluvium. Hints of low velocities are also present near the burial depth (~364 m). The reflection seismic data show a prominent subhorizontal reflector at ~300 m depth, a short-curved reflector at ~200 m, and a high-amplitude reflector at ~50 m depth. Comparisons of the reflection sections to synthetic data and borehole stratigraphy suggest that these features correspond to the alluvium–tuff contact, the partial collapse cavity, and the spalled layer, respectively.

58 GEOSCIENCES↗

Legacy surface change analyses of the 1993 Rock Valley earthquake sequence for direct comparison to planned NA-22 underground conventional high-explosive experiments (Source Physics Experiment 3 (SPE3) - RV/DC Task 1.3 FY22 Final Report)

Recent work under two previous phases of the NNSA NA-22 Source Physics Experiment (SPE) have shown that underground chemical high-explosive experiments can produce detectable surface changes that differ in spatial extent and vertical magnitude depending on the geologic media at the site (Schultz-Fellenz et al., 2018; 2020; Crawford et al., 2021). Neither of these two prior phases of SPE identified natural earthquake-related surface effects in the same region that occurred at a similar depth as the explosive experiments for direct comparison. While earthquakes can also produce surface changes that are detectable using remote sensing data analyses, it is expected that the pattern and spatial extent of surface changes would vary between earthquakes and explosions. However, no direct-observed surface-change signature comparison between earthquakes and explosions has ever been performed. The SPE Phase 3 Rock Valley/Direct Comparison (RV/DC) program presents a unique opportunity to investigate and characterize co-occurrence of both earthquakes and explosions. In this project, we worked to address five tasks in a workflow, as follows: 1. Identify and obtain existing high-resolution legacy satellite and aerial imagery as close in time before and after the 1993 Rock Valley earthquake sequence to temporally constrain the analyses. 2. Transform these pre-earthquake and post-earthquake datasets into digital elevation models (DEMs) using geospatial analysis software packages (e.g., ArcGIS, Agisoft Metashape, and Google Earth Engine). 3. Perform DEM differencing analyses to assess and quantify earthquake-related changes from the 1993 sequence, and develop map products that visualize these analyses. 4. Use the analyses from (3) to: (a) assess spatial distribution and magnitude of surface changes due to the 1993 Rock Valley earthquake sequence, and (b) determine parameters of forthcoming, planned explosion-related surface change data collection from sensors mounted on unmanned aerial vehicles (UAVs) (e.g., spatial extent of collection, design and density of survey control, sensors to deploy, forward speed and line spacing of UAV flight lines, flight altitude). 5. Develop a summary report on the analyses, including how the analyses define parameters and identify focus areas for any future surface change analytical field campaigns related to the explosive experiment. Analyzing these legacy data and identifying whether they can detect any surface changes related to the earthquake sequence facilitates opportunities for direct signature comparison of surface change from explosions at one location, which has never previously been performed. Comparing the surface change signatures from a co-located and depth-equivalent earthquake and an explosion could help to advance remote sensing event discrimination techniques. This report summarizes the work completed toward this ambitious goal.

42 ENGINEERING↗

Discrimination of Small Earthquakes and Buried Single-Fired Chemical Explosions at Local Distances (<150 km) in the Western United States from Comparison of Local Magnitude ( M L ) and Coda Duration Magnitude ( M C )

Seismologists distinguish underground nuclear explosions from more commonly occurring earthquakes using moment tensor inversion, high-frequency P/S amplitude ratios, m b :M s comparisons, and P–pP differential travel times. These methods are generally successful for large seismic events (M >3-4) well recorded at regional-to-teleseismic distances (>150 km); however, it is unclear whether they can be modified to work for small events (M <3) well recorded only at local distances (<150 km). Here, we evaluate a recently proposed, local-distance seismic source discriminant—the difference between local magnitude (M L ) and coda duration magnitude (M C )—using seismograms of earthquakes and buried, single-fired chemical explosions recorded in three regions of the western United States. The quantity M L –M C was previously found to be sensitive to source depth, effectively discriminating mine blasts, induced earthquakes, and very shallow tectonic earthquakes from deeper crustal earthquakes. In this study, we report the first evaluation of M L –M C as a depth discriminant using data from buried, single-fired explosions that, unlike the seismic sources studied earlier, are good analogs for underground nuclear explosions. We find that even when using generic, uncalibrated methods of assigning magnitudes, M L –M C separates single-fired explosions and earthquakes. The area under the receiver operating characteristic curve is 0.92 for 19 explosions and 14 earthquakes in Washington, 0.90 for 22 explosions and 90 earthquakes in Wyoming, and 0.99 for three explosions and 149 earthquakes in Nevada. Finally, M L :M C comparisons have the potential to enhance discrimination based on high-frequency P/S amplitudes ratios—which perform less well at local than regional distances—because the two metrics have complementary sensitivities.

58 GEOSCIENCES↗

Modeling the use of mobile modular gas samplers in near-field detection

The Wireless Independent Noble Gas Sampler (WINGS) is a mobile, modular gas sampling system designed for use in low infrastructure environments. In the case of a suspected underground nuclear explosion (UNE), WINGS units could be deployed to detect noble gases emanating from the suspected explosion site or identify gases originating from another source entering the local area from offsite. In conclusion, this work uses the atmospheric transport modeling tool inline WRF-HYSPLIT to test the effectiveness of WINGS units around a suspected UNE site as a function of sampler density and distance from the emission point of gases from a UNE.

Air sampling↗

UNESE Argon-39 Measurement Techniques: Developing an above-ground Argon-39 Measurement Capability

The Underground Nuclear Explosion Signatures Experiment (UNESE) sought to use 37 Ar as a tracer for measuring noble-gas migration in the soil surrounding historic Underground Nuclear Explosions (UNE). One unexpected observation was the presence of the much longer-lived isotope 39 Ar from historic UNEs. Quantifying the activity of 39 Ar proved difficult due to the lack of capability to measure significantly-above-background levels of 39 Ar and a general discomfort to repeatedly expose Ultra-Low-Background Proportional Counters (ULBPCs) to significant radioactivity. Because the whole-air samples collected for the 37 Ar tracer measurement were already being measured on the above-ground argon capability, it was decided to expand that capability to include 39 Ar. This document describes the efforts required to achieve quantitative reporting of the 39 Ar backgrounds measured during experiments at the sites of the historic Barnwell and Disko Elm UNEs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

X-ray Computed Tomography on UNESE Core: FY2020 Data Report to Support Fracture and Multiphase Fluid Flow Studies

Natural and induced fractures are potential preferential pathways for migration of radioactive gases to earths surface from underground nuclear explosions (UNEs). This report documents X-ray computed tomography (XRCT) imaging on 26 samples of rock core that was collected to support the Underground Nuclear Explosion Signatures Experiment (UNESE) program. The XRCT datasets are intended to help fill a data gap on the three-dimensional (3D) characteristics of natural and/or induced fractures at the centimeter and smaller scale, which may strongly influence multiphase fluid flow and transport properties of preferential flow paths and interaction with the matrix of the surrounding host rock. Pre- and post-UNE rock samples were carefully chosen to enable comparison of fractures as a function of lithologic and petrophysical properties, as well as distance to the past UNEs. This report serves as documentation for the data, including an introduction with the research motivation, a methods and materials section, descriptions of the XRCT datasets without post-processing, and recommendations for 3D quantification via image analysis and digital rock physics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Comparison of near-background concentrations of Argon-37 and Xenon-133 in the Atmosphere

Radioisotopes of noble gases xenon and argon can be important indicators of underground nuclear explosions. The Comprehensive Nuclear-Test-Ban Treaty (CTBT) includes monitoring capabilities to identify potential nuclear explosions conducted in violation of the CTBT. This monitoring currently focuses on measurement of the xenon isotopes 133Xe, 135Xe, 131mXe and 133mXe. However, it is predicted that within 100 days of an underground nuclear explosion (UNE) 37Ar would be released to the atmosphere at higher concentrations than xenon isotopes (Haas et al. 2010) and with a higher signal to background ratio, depending on the radioxenon background levels. Therefore, inclusion of 37Ar measurement capabilities at atmospheric International Monitoring System (IMS) stations may represent an improvement in the capability to detect a nuclear explosion. At an IMS station, it could be difficult to determine what constitutes an elevated concentration as a result of an UNE without first understanding the expected range of background concentrations. This work describes our analysis of atmospheric samples for 37Ar to evaluate the range of background concentrations. Samples were collected at multiple locations, with approximately half coming from a sampler co-located with an IMS xenon monitoring station (RN75). The range of 37Ar concentrations measured in atmospheric air samples was relatively narrow; for samples considered detectable, the minimum and maximum measured concentrations were 0.56 and 2.3 mBq/m3, respectively. Comparison of 37Ar and 133Xe concentrations measured at the IMS station indicated some correlation between the measured concentrations. The results presented here demonstrate the capability to detect background concentrations of 37Ar in atmospheric air and provide a basis for potential implementation of 37Ar monitoring at IMS stations.

Fritz, Brad G.↗

Measurement of Argon-37 Background Concentrations in the Atmosphere

Radioisotopes of noble gases, primarily xenon and argon, are important indicators of underground nuclear explosions. The Comprehensive Nuclear-Test-Ban Treaty (CTBT) includes monitoring capabilities to identify potential nuclear explosions conducted in violation of the CTBT. This monitoring currently focuses on measurement of the xenon isotopes 133 Xe, 135 Xe, 131 mXe and 133 mXe. However, it is predicted that within 100 days of an underground nuclear explosion (UNE) 37 Ar would be released to the atmosphere at higher concentrations than xenon isotopes (Haas et al. 2010). Therefore, inclusion of 37 Ar measurement capabilities at atmospheric International Monitoring System (IMS) stations could improve the capability to detect a nuclear explosion. At an IMS station, it could be difficult to determine what constitutes an elevated concentration as a result of an UNE without first understanding the expected range of background concentrations. This work describes our analysis of atmospheric samples for 37 Ar to evaluate the range of background concentrations. The range of 37 Ar concentrations measured in atmospheric air samples was relatively narrow ranging between less than 1 mBq/m 3 to greater than 2 mBq/m 3 . The average atmospheric concentration of all samples was 0.98 mBq/m 3 , with a standard deviation of 0.5 mBq/m 3 . This is consistent with previous published results. The results presented here demonstrate the capability to detect background concentrations of 37 Ar in atmospheric air.

61 RADIATION PROTECTION AND DOSIMETRY↗

Source Characterization of the Declared North Korean Nuclear Tests From Regional Distance Coda Wave Spectral Ratios

Abstract Seismic observations of underground nuclear explosions provide crucial data on source yield and depth that cannot easily be estimated from other geophysical methods. However, it is difficult to obtain reliable yield estimates for test sites for which we do not have direct seismic calibration experiments. To obtain source information from uncalibrated sites and paths, local and regional seismic records of six, proximal, declared underground nuclear explosions in North Korea are used to compute spectral ratios of narrow‐band waveform envelopes of body‐wave coda that remove path and site effects to reveal precise, relative source moment. The yields of these explosions are obtained from the observed source ratios by simultaneously fitting the classical source model of Mueller and Murphy (1971), https://doi.org/10.1785/bssa0610061675 to all event pairs. The source model provides an impressive fit to the observations considering that the P phase coda derived source spectral ratios did not require a prior knowledge of the source or regionally calibrated corrections to be applied to the data. However, the observed corner frequencies from S wave coda spectral ratios are lower than the source model predictions, but are well fit by models calculated using the corner frequency consistent with the Fisk conjecture. The results presented here provide novel constraints on the spectral distributions of the energy radiated by the sources of the DPRK test series, and allows for an independent evaluation of existing estimated source model parameters.

58 GEOSCIENCES↗

UNESE Phase 2: Injection and measurement of gaseous tracers at U-12p Tunnel

In June 2018 radioactive ( 37 Ar and 127 Xe) and stable (SF 6 ) tracers were injected into the chimney of the Disko Elm underground nuclear explosion, located in the U-12p tunnel complex. After the injection sampling was performed from locations within the tunnel, a borehole drilled from the surface, and from the surface of the mesa. Additionally, the 127 Xe activity within the chimney was continuously monitored for two months after the injection. While none of the injected radiotracers were detected in the monitoring borehole, elevated levels of 39 Ar (an activation product produced by underground nuclear explosions) were detected and quantified in each sample.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Reconstruction of Pre- and Post-Test Surfaces at the U-20az BARNWELL Site, Nevada National Security Site

Technical report on results of photogrammetric reconstruction task at NNSS. Modern photogrammetric methods were applied to historic aerial photographs taken six days before and four days after the BARNWELL underground nuclear explosion (UNE) to recreate the pre-UNE topography and the fresh post-UNE surface. The surfaces were compared to each other to characterize surface effects of the BARNWELL UNE. Topographic changes between the pre-UNE and post-UNE surfaces were limited to minor rock fall and intentional changes related to construction activities. Methods developed in this study could potentially be used across a much wider area to discover and catalogue as yet undescribed subtle underground nuclear explosion surface effects.

58 GEOSCIENCES↗

UNESE Data Analysis-Disko Elm Gas Transport Characteristics (NA-22 Quarterly Report)

Recently published papers have demonstrated the impact of sorption mechanisms on gas transport in the subsurface following an underground nuclear explosion. To include sorption processes in the multi-physics models of noble gas production and transport, researchers have conducted experiments of gas sorption in geologic materials and calibrated a linear (Henry’s law) sorption model for all gases in various geologic materials. As a part of Task #4, we requested the experimental data from Sandia Laboratories and derived a dual-mode (Henry’s absorption and Langmuir adsorption) sorption model. The derived model better describes the nonlinear sorption mechanism with much higher fitness with the experimental data. When gas concentrations are relatively low (in the conditions of field experiments or underground nuclear explosions) compared to those used sorption experiments, the nonlinear sorption model can be approximated as a linear sorption model. The equivalent Henry’s constant, which equals the product of Langmuir capacity and affinity, is systematically higher than the literature values. Further analyses and comparison of sorption models are planned in FY2020 Q4 and the sorption model will be coupled in models of nonisothermal multiphase transport for studying noble gas detectability.

58 GEOSCIENCES↗

Dynamic Characteristics of Multistory Buildings

In order to evaluate unit stresses or to estimate possible damage in multistory buildings from earthquake or ground motion from underground nuclear explosion it is essential to determine how much of the dynamic response at any level is due to various types of freedom. The effects of shear deformation between floors, joint rotation, over-all flexure, and ground compliance are considered. Joint rotation and over-all flexure are evaluated over a wide range of building characteristics with simply determined indices. Further, a period synthesis concept and a pseudo-stiffness procedure are proposed to enable the determination of periods, mode shapes and stiffnesses with consideration of joint rotation, over-all flexure and ground compliance while performing simple labor-saving analyses for assumed rigid-floor shear buildings. The first three modes of vibration are considered in elastic free vibration. The buildings and models are symmetric in plan without torsional coupling.

42 ENGINEERING↗

A multi-Physics Experiment for Low-Yield Nuclear Explosion Monitoring

A series of multi-physics experiments, referred to as Physics Experiment 1 (PE1) is underway at the United States’ Nevada National Security Site (NNSS). The PE1 series includes detonations of three underground chemical explosions in P-tunnel, with fully coupled (PE1 A), partially decoupled (PE1 D L ), and fully decoupled (PE1 B) emplacements. Canisters with gas tracers are imbedded in the explosives, and the tracers are released when the canister is destroyed by the detonation. A dedicated electromagnetic (EM) experiment (EMX) generates well-characterized EM signals at an underground location near the chemical explosive experiments. A series of atmospheric experiments (METEX, REACT, and METREX) release smoke and radioactive tracers around Aqueduct Mesa to test gas transport in complex topography. Each of the chemical explosive experiments includes a network of sensors to record seismic, acoustic, and electromagnetic waves, measurement of atmospheric conditions, and air sample collection for measurement of tracer concentration. EMX records EM signals underground and on the surface of Aqueduct Mesa. METEX, REACT, and METREX include measurement of atmospheric condition, as well as tracking smoke releases. REACT and METREX add low-level radioactive gas tracers to the atmospheric releases.

58 GEOSCIENCES↗

Multiphenomenology explosion monitoring (MultiPEM): a general framework for data interpretation and yield estimation

SUMMARY An underground nuclear explosion (UNE) couples mechanical energy into crustal rock, which propagates as seismic and acoustic waves. These different physical phenomena transport, by different pathways, to standoff detectors at varying distances. The transport pathways attenuate the original signal but in different ways. Enabled by correct statistical weighting, signal attenuation models can be used to combine these disparate sensor data to estimate the yield of an UNE. Contemporaneous statistical models, used in yield estimation, can be improved with an advanced partition of error for these physical signal propagation models. We present an advanced multivariate approach to error modelling of multiphenomenology physical signatures. In addition to measurement error, our error model represents physical model biases as random with a physics-based covariance structure. To illustrate this proposed framework, we demonstrate the estimation of explosion yield using openly available seismic and acoustic data from chemical single-point explosions.

Williams, Brian J.↗

Review of Cavity Radius and Chimney Height Information for Underground Nuclear Tests at Nevada National Security Site

Underground nuclear explosions release an immense amount of thermal and mechanical energy that alter the geologic materials surrounding the “working point” (WP = underground location) of the detonation. Explosion effects produce mechanically and thermally altered zones such as the “cavity,” “chimney,” and “melt glass puddle,”as well as concentric zones of crushed, compressed, and fractured rock extending away from the cavity wall.

58 GEOSCIENCES↗

Gas diffusion through variably-water-saturated zeolitic tuff: Implications for transport following a subsurface nuclear event

Noble gas transport through geologic media has important applications in the characterization of underground nuclear explosions (UNEs). Without accurate transport models, it is nearly impossible to distinguish between xenon signatures originating from civilian nuclear facilities and UNEs. Understanding xenon transport time through the earth is a key parameter for interpreting measured xenon isotopic ratios. One of the most challenging aspects of modeling gas transport time is accounting for the effect of variable water saturation of geological media. In this study, we utilize bench-scale laboratory experiments to characterize the diffusion of krypton, xenon, and sulfur hexafluoride (SF6) through intact zeolitic tuff under different saturations. Here, we demonstrate that the water in rock cores with low partial saturation dramatically affects xenon transport time compared to that of krypton and SF6 by blocking sites in zeolitic tuff that preferentially adsorb xenon. This leads to breakthrough trends that are strongly influenced by the degree of the rock saturation. Xenon is especially susceptible to this phenomenon, a finding that is crucial to incorporate in subsurface gas transport models used for nuclear event identification. We also find that the breakthrough of SF6 diverges significantly from that of noble gases within our system. When developing field scale models, it is important to understand how the behavior of xenon deviates from chemical tracers used in the field, such as SF6 (Carrigan et al., 1996). These new insights demonstrate the critical need to consider the interplay between rock saturation and fission product sorption during transport modeling, and the importance of evaluating specific interactions between geomedia and gases of interest, which may differ from geomedia interactions with chemical tracers.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Preferential adsorption of noble gases in zeolitic tuff with variable saturation: A modeling study of counter-intuitive diffusive-adsorptive behavior

Noble gas transport through geologic media has important applications in the prediction and characterization of measured gas signatures related to underground nuclear explosions (UNEs). Retarding processes such as adsorption can cause significant species fractionation of radionuclide gases, which has implications for measured and predicted signatures used to distinguish radioxenon originating from civilian nuclear facilities or from UNEs. Accounting for the effects of variable water saturation in geologic media on tracer transport is one of the most challenging aspects of modeling gas transport because there is no unifying relationship for the associated tortuosity changes between different rock types, and reactive transport processes such as adsorption that are affected by the presence of water likewise behave differently between gas species. In this study, we perform numerical diffusive-adsorptive transport simulations to estimate gas transport parameters associated with bench-scale laboratory diffusion cell experiments measuring breakthrough in zeolitic and non-zeolitic rocks for a gaseous mixture of xenon, krypton, and SF 6 at varying degrees of water saturation (S w ). Counter-intuitive transport behavior was observed in the zeolitic rock experiments whereby breakthrough concentrations were significantly higher when the core was partially saturated (S w = 17 %) than under dry (S w = 0 %) conditions. Breakthrough of xenon was especially retarded in the dry core – likely due to comparatively high affinity of xenon for zeolitic adsorption sites – and estimated effective diffusion coefficients for all gases were approximately an order of magnitude lower than what is predicted by porosity-tortuosity models. We propose the counter-intuitive behavior observed is because water infiltration of zeolite nanopores reduces both the adsorptive capacity of the rock and the tortuosity of connected flow paths. We developed a two-site competitive kinetic Langmuir adsorption reaction for the porous media transport simulator in order to constrain transport parameters within zeolitic tuff, where differential adsorption to zeolite and non-zeolite pores was observed. We determined that liquid saturation-dependent diffusive-adsorptive transport is affected by subtle and at times competing processes that are specific to different gases, which have a significant overall influence on effective transport parameters.

58 GEOSCIENCES↗