Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Explosives”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Measurements of Gamow-Teller transitions from 59 Co via the 59 Co ⁢(𝑡, 3 He +𝛾) charge-exchange reaction and its application to the stellar electron-capture rates

Electron-capture reactions on iron-group nuclei play a crucial role in the late stages of massive star evolution. Since stellar evolution simulations depend on accurate electron-capture rates—which are highly sensitive to the detailed Gamow-Teller (GT) strength distributions—reliable theoretical models are essential. However, experimental data on GT strength distributions are scarce. High-resolution measurements are therefore vital for benchmarking and improving these theoretical calculations. To provide high-resolution data on Gamow-Teller strength distributions of iron-group nuclei and to compare these results with theoretical calculations within this mass region. Differential cross sections for the 59 Co ⁢(𝑡, 3 He)⁢ 59 Fe charge-exchange reaction at 115 MeV/u were measured using the S800 spectrometer. Furthermore, to resolve individual levels that are not distinguishable in the S800 particle singles data, coincident 𝛾 rays from the 59 Fe residual nucleus were detected by using the Gamma-Ray Energy Tracking In-beam Nuclear Array 𝛾-ray tracking array. Here, the Gamow-Teller transition strength distribution from the ground state of 59 Co to 59 Fe was extracted up to an excitation energy of 10 MeV. Additionally, transition strengths for several low-lying states were determined from coincident 𝛾-ray measurements. Electron-capture rates calculated using the present data indicate that these low-lying states contribute significantly to the overall rates in relevant stellar environments. The experimental results show reasonable agreement with theoretical predictions based on both shell-model and projected shell-model calculations. High-resolution data on Gamow-Teller strength distributions—particularly for individual low-lying states—are essential for accurately determining electron-capture rates in iron-group nuclei. Coincident 𝛾-ray measurements provide a powerful tool for obtaining such detailed information. While the present work demonstrates that shell-model calculations successfully reproduce the experimental results, such comparisons are scarce and more experimental data are desirable.

59 ≤ A ≤ 89↗
DOE OSTI2026

Enhanced antineutrino emission from 𝛽 decay in core-collapse supernovae with self-consistent weak decay rates

Nuclear weak-interaction rates are known to exert a prominent effect in the late stages of stellar collapse. Despite their importance, most studies to date on core-collapse supernovae (CCSNe) have focused primarily on the effects of electron captures, neglecting 𝛽 decay contributions. In this work, we present the first CCSNe simulation incorporating global 𝛽 decay rates from a microscopic theory. These are enabled by a large-scale evaluation of both electron capture and 𝛽 decay rates, obtained self-consistently utilizing the relativistic energy density functional theory and finite-temperature quasiparticle random-phase approximation. Including 𝛽 decay leads to a dramatic enhancement of the pre-bounce antineutrino signal as the antineutrino emissivity increases by more than two orders of magnitude and the luminosity by a factor of 50 relative to thermal emission alone, while the average antineutrino energy increases by over 1 MeV. Furthermore, it is expected that these new rates could help us constrain the model uncertainties related to weak-interaction processes, improving the prediction of antineutrino signal during the final stages of stellar death.

Beta decay↗
DOE OSTI2020

Influence of 73 Rb on the ashes of accreting neutron stars

We find that the proton separation energy, S(p), of 73 Rb is –640(40) keV, deduced from the observation of β-delayed ground-state protons following the decay of 73 Sr. This lower-limit determination of the proton separation energy of 73 Rb coupled with previous upper limits from nonobservation, provides a full constraint on the mass excess with ΔM ( 73 Rb) = –46.01 ± 0.04 MeV. With this new mass excess and the excitation energy of the J π = 5/2 – isobaric-analog state (T = 3/2) in 73 Rb, an improved constraint can be put on the mass excess of 73Sr using the isobaric-multiplet mass equation (IMME), and we find ΔM( 73 Sr) = –31.98 ± 0.37 MeV. Furthermore, these new data were then used to study the composition of ashes on accreting neutron stars following Type I x-ray bursts. Counterintuitively, we find that there should be an enhanced fraction of A > 102 nuclei with more negative proton separation energies at the 72 Kr rp-process waiting point. Larger impurities of heavier nuclei in the ashes of accreting neutron stars will impact the cooling models for such astrophysical scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Calculation of the C 12 + C 12 sub-barrier fusion cross section in an imaginary-time-dependent mean field theory

The 12 C + 12 C sub-barrier fusion cross section is calculated within the framework of a time-dependent Hartree-Fock-based classical model using the Feynman path-integral method. The modified astrophysical S* factor is compared to direct and indirect experimental results. A good agreement with the direct data is found. In the lower-energy region where recent analyses of experimental data obtained with the Trojan horse method (THM) lead to contrasting results, the model predicts a nonresonant S* factor half-way between those results. Low-energy resonances revealed in the THM data are added to the calculation, and the relative reaction rate in the Gamow region is calculated. In particular, including 0 + resonances result in some agreement with the THM data. Here, the role of different resonances is discussed in detail, and their influence on the reaction rate at temperatures relevant to stellar evolution is investigated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Modified Brink-Axel hypothesis for astrophysical Gamow-Teller transitions

Weak interaction charged current transition strengths from highly excited nuclear states are fundamental ingredients for accurate modeling of compact object composition and dynamics, but are difficult to obtain either from experiment or theory. For lack of alternatives, calculations have often fallen back upon a generalized Brink-Axel hypothesis, that is, assuming the strength function (transition probability) is independent of the initial nuclear state but depends only upon the transition energy and the weak interaction properties of the parent nucleus ground state. Here we present numerical evidence for a modified “local” Brink-Axel hypothesis for Gamow-Teller transitions for pf-shell nuclei relevant to astrophysical applications. Specifically, while the original Brink-Axel hypothesis does not hold globally, strength functions from initial states nearby in energy are similar within statistical fluctuations. This agrees with previous work on strength function moments. Using this modified hypothesis, we can tackle strength functions at previously intractable initial energies, using semiconverged initial states at arbitrary excitation energy. Our work provides a well-founded method for computing accurate thermal weak transition rates for medium-mass nuclei at temperatures occurring in stellar cores near collapse. We finish by comparing results to previous calculations of astrophysical rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Finite-temperature electron-capture rates for neutron-rich nuclei near N = 50 and effects on core-collapse supernova simulations

In this paper the temperature dependence of stellar electron-capture (EC) rates is investigated, with a focus on nuclei near N = 50, just above Z = 28, which play an important role during the collapse phase of core-collapse supernovae (CCSN). Two new microscopic calculations of stellar EC rates are obtained from relativistic and nonrelativistic finite-temperature quasiparticle random-phase approximation approaches, for a conventional grid of temperatures and densities. In both approaches, EC rates due to Gamow-Teller transitions are included. In the relativistic calculation, contributions from first-forbidden transitions are also included and add strongly to the EC rates. The new EC rates are compared with large-scale shell-model calculations for the specific case of 86 Kr, providing insight into the finite-temperature effects on the EC rates. At relevant thermodynamic conditions for core collapse, the discrepancies between the different calculations of this paper are within about one order of magnitude. Numerical simulations of CCSN are performed with the spherically symmetric gr1d simulation code to quantify the impact of such differences on the dynamics of the collapse. These simulations also include EC rates based on two parametrized approximations. A comparison of the neutrino luminosities and enclosed mass at core bounce shows that differences between simulations with different sets of EC rates are relatively small (≈ 5 %), suggesting that the EC rates used as inputs for these simulations have become well constrained.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Bayesian estimation of the $S$ factor and thermonuclear reaction rate for 16 O(p, γ) 17 F

The 16 O(p, γ) 17 F reaction is the slowest hydrogen-burning process in the CNO mass region. Its thermonuclear rate sensitively impacts predictions of oxygen isotopic ratios in a number of astrophysical sites, including AGB stars. The reaction has been measured several times at low bombarding energies using a variety of techniques. The most recent evaluated experimental rates have a reported uncertainty of about 7.5% below 1 GK. However, the previous rate estimate represents a best guess only and was not based on rigorous statistical methods. We apply a Bayesian model to fit all reliable 16 O(p, γ) 17 F cross section data, and take into account independent contributions of statistical and systematic uncertainties. The nuclear reaction model employed is a single-particle potential model involving a Woods-Saxon potential for generating the radial bound state wave function. The model has three physical parameters, the radius and diffuseness of the Woods-Saxon potential, and the asymptotic normalization coefficients (ANCs) of the final bound state in 17 F. Here, we find that performing the Bayesian S -factor fit using ANCs as scaling parameters has a distinct advantage over adopting spectroscopic factors instead. Based on these results, we present the first statistically rigorous estimation of experimental 16 O(p, γ) 17 F reaction rates, with uncertainties (±4.2%) of about half the previously reported values.

6 ≤ A ≤ 19↗
DOE OSTI2020

First Direct Measurement of Mg (α,p) 25 Al and Implications for X-Ray Burst Model-Observation Comparisons

Type-I x-ray bursts can reveal the properties of an accreting neutron star system when compared with astrophysics model calculations. However, model results are sensitive to a handful of uncertain nuclear reaction rates, such as 22 Mg(α , p). We report the first direct measurement of 22 Mg(α , p), performed with the Active Target Time Projection Chamber. The corresponding astrophysical reaction rate is orders of magnitude larger than determined from a previous indirect measurement in a broad temperature range. Overall, our new measurement suggests a less-compact neutron star in the source GS1826-24.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Urca Cooling in Neutron Star Crusts and Oceans: Effects of Nuclear Excitations

We report the excited-state structure of atomic nuclei can modify nuclear processes in stellar environments. In this Letter, we study the influence of nuclear excitations on Urca cooling (repeated back-and-forth β decay and electron capture in a pair of nuclear isotopes) in the crust and ocean of neutron stars. We provide for the first time an expression for Urca process neutrino luminosity which accounts for a thermal Boltzmann distribution of excited states in both members of an Urca pair. We use our new formula with state-of-the-art nuclear structure inputs to compute neutrino luminosities of candidate Urca cooling pairs. Our nuclear inputs consist of the latest experimental data supplemented with calculations using the projected shell model. We show that, in contrast to previous results that only consider the ground states of both nuclei in the pair, our calculated neutrino luminosities for different Urca pairs vary sensitively with the environment temperature and can be radically different from those obtained in the one-transition approximation. We find that nuclear excitations can lead to an enhancement in total Urca neutrino luminosities in the accreted neutron star crust by about 5 times as compared with the previous Urca results, which is expected to cause significant observational effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Raman Spectroscopy Characterization of Nuclear Materials and Soot from High Explosives

Vibrational and luminescence spectroscopy were used in the investigation of U and Pu-bearing compounds and soot from high explosives testing. Key spectroscopic signatures of UF4 and its hydrates, PuO 2 , and soot from different explosives will be presented. The United States Government retains and the publisher, by accepting this article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this work, or allow others to do so, for United States Government purposes.

Villa-Aleman, Eliel↗
DOE OSTI2021

Exploring the Effects of Emplacement Conditions on Explosion P/S Ratios across Local to Regional Distances

High-frequency (~> 2 Hz) seismic P/S amplitude ratios are well-established as a discriminant to distinguish between natural earthquakes and underground explosions at regional distances (~200–1500 km). As research shifts toward identifying lower-yield events, work has begun to investigate the potential of this discriminant for use at local distances (<200 km), in which initial results raise questions about its effectiveness. Here, we utilize data from several chemical explosion experiment series at the Nevada National Security Site in southern Nevada in the United States to study explosion Pg/Lg ratios across the range of local to regional distances. Additionall, the experiments are conducted over differing emplacement conditions, with contrasting geologies and a variety of yields and depths of burial, including surface explosions. We first establish the similarities of Pg/Lg ratios from chemical explosions to those from historic nuclear tests and conclude that, as previous data have suggested, chemical explosion ratios are good proxies for nuclear tests. We then examine Pg/Lg ratios from the new experiment series as functions of distance, yield, depth of burial, and scaled depth of burial (SDOB). At far-local and regional distances, we observe consistently higher ratios from hard-rock explosions compared to ones in a weaker dry alluvium medium, consistent with prior regional distance results. No other trends with yield, depth of burial, or SDOB are strongly evident. Scatter in the observed ratios is very high, particularly at the shortest event-to-station distances, suggesting that small-scale path effects play a significant role. On average, the local distance explosion Pg/Lg ratios show remarkable consistency across all the variations in emplacement. Explosion source models will need to reproduce these results.

58 GEOSCIENCES↗
DOE OSTI2025

PySolate : A Python‐Based Thresholding Tool to Denoise or Designal Seismic Waveforms Based on the Continuous Wavelet Transform

PySolate is a Python‐based toolset that implements the continuous wavelet transform and nonlinear thresholding operations to denoise or designal seismic data, following Langston and Mousavi (2019). This filtering approach can remove microseismic noise to isolate intermediate‐period seismic signals that are key to enabling full‐waveform modeling and analysis of smaller‐magnitude regional events. This approach is best for the application to signals with frequency or time separation of signal and noise, in contrast to Fourier analysis, which is effective when signal and noise are separated in frequency. We demonstrate the Python toolset using the six announced Democratic People’s Republic of Korea declared nuclear tests, showing the effectiveness of isolating the seismic signal compared to standard bandpass filtering. In conclusion, we also demonstrate the ease of using the toolset with any Python processing tools.

Asia↗
DOE OSTI2022

Developing the S 32 ( p , d ) S * 31 ( p ) ( γ ) reaction to probe the P 30 ( p , γ ) S 31 reaction rate in classical novae

Background: The 30 P(p,γ) 31 S reaction rate is one of the largest remaining sources of uncertainty in the final abundances of nuclei created in a classical nova involving a ONe white dwarf. The reaction rate directly influences silicon isotopic ratios, which are used as identifiers of presolar grains with nova origins. Additionally, the uncertainty in the 30 P(p,γ) 31 S reaction rate has been found to limit the use of nova nuclear thermometers based on observations of elemental ratios in nova ejecta. Purpose: Reduce uncertainties in the nuclear data for proton-unbound states in 31 S, which act as resonances for the 30 P(p,γ) 31 S reaction at classical nova temperatures, and develop a technique for high efficiency, high-resolution reaction-decay coincidence measurements. Methods: Here, the 32 S(p,d) 31 S reaction was used to populate the states of interest in 31 S. The experiment was performed at the Texas A&M Cyclotron Institute using the LLNL Hyperion array for the detection of charged particles and γ rays. A downstream silicon telescope was used to select reaction deuterons, and a single upstream silicon detector was used to measure protons emitted in the decay of unbound 31 S levels. Results: Several states in 31 S above the proton separation energy were observed to have been populated. Decay protons from the resonant states in 31 S were identified as events in the upstream silicon detectors that came in coincidence with deuterons in the downstream telescope. Protons emitted from these states were measured and branching ratios extracted. Conclusions: While no new reaction rate is derived, spin-parity assignments for several higher-lying proton unbound states have been confirmed. Measured p 0 branching ratios for these levels have been compared to previous measurements with good agreement, and in some cases provided a reduction in uncertainty. The previously identified T = 3/2 state may have been incorrectly assigned a large p 0 branching ratio in a previous measurement. The technique of measuring reaction-decay coincidences with a particle-gamma setup appears promising.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The Low-Yield Nuclear Monitoring (LYNM) Experimental Science Plan

The Low-Yield Nuclear Monitoring (LYNM) Program is a long-term NNSA research and development effort designed to improve the United States’ explosion monitoring capabilities, particularly with respect to low-yield and potentially evasive underground nuclear testing. The LYNM Program focuses on researching, discovering, and exploiting unique and useful signatures, from all available technologies and sensors (e.g., seismic, acoustic, electromagnetic, gases, and particulates (both stable and radioactive)). Four Department of Energy laboratories participate in LYNM and together are referred to as the ‘quad-lab’. The R&D program execution is performed under NNSA defined structures known as “ventures”. Four science ventures were established: 1) Explosion Source Functions; 2) Containment of Low-Yield Underground Tests; 3) Local Signatures; 4) Dynamic Monitoring Networks. To organize and execute the large field scale experiment a fifth venture was established: 5) Physics Experiment One (PE-1). As part of the LYNM Program, a series of experiments are planned at various scales and levels of venture involvement. These vary from those that involve a subset of labs and/or LYNM ventures (e.g., small experiments), to full quad-lab LYNM Program field-scale integrated experiments. Such experiments may involve chemical explosions with tracer materials or other means of simulating the expected signals from a nuclear explosion. The LYNM Program does not conduct actual nuclear explosions. Since 1992, the U.S. has observed a moratorium on underground nuclear explosions. This document is intended to provide the underlying scientific basis for the LYNM planned experimental work. Each specific LYNM experiment will develop a goals, objectives, and requirements (GOR) plan following the guidance in this document. The LYNM technical staff will define the numbers and types of experiments required over the course of the Program based on technical needs and within funding constraints. As with any scientific experiment series, the number and types of experiments may change based upon the experimental results obtained. An experiment that agrees with models/codes/software signature predictions may need fewer repetitions/variations, depending upon the level of statistical rigor desired, as compared to one in which the predictions and experimental data do not match. The large LYNM field-scale integrated experiments require the longest lead-time for planning, and these are discussed in more detail near the end of this document.

58 GEOSCIENCES↗
DOE OSTI2023

Secure Storage: Historical Documentation of TA-08-0032, TA-11-0036, TA-22-0016, TA-22-0023, TA-22-0025, TA-22-0035, TA-37-0006, TA-37-0009, and TA-37-0020

The U.S. Department of Energy, National Nuclear Security Administration, Los Alamos Field Office (NA-LA), is pursuing the decommissioning and demolition (D&D) of facilities contaminated with high-explosives residues at the Los Alamos National Laboratory (Laboratory or LANL). This effort affects nine facilities associated with high-explosives and detonator research, development, and storage: Technical Area (TA) 8 Facility 32 (TA-08-0032), TA-11-0036, TA-22-0016, TA-22-0023, TA-22-0025, TA-22-0035, TA-37-0006, TA-37-0009, and TA-37-0020. All nine facilities proposed for D&D have been evaluated for listing in the National Register of Historic Places (NRHP) and determined eligible. NA-LA previously requested the State Historic Preservation Officer (SHPO) to concur with the NRHP-eligibility determinations of these nine properties presented in four reports: (1) TA-08-0032 was determined eligible for listing in the NRHP in the report, From Ranching to Radiography: An Assessment of Historic Buildings at Anchor West Site (TA-8), Vol. 1 (McGehee et al. 2008a). The SHPO concurred with this eligibility determination on November 26, 2008. (2) TA-11-0036 was determined eligible for listing in the NRHP in the report, ESA Division’s Five-Year Plan: Consolidation and Revitalization at Technical Areas 3, 8, 11, and 16, Vol. 1) (McGehee et al. 2003). The SHPO concurred with this eligibility assessment on June 22, 2003. (3) TA-22-0016, TA-22-0023, TA-22-0025, and TA-22-0035 were determined eligible for listing in the NRHP in the report, DX Division’s Facility Strategic Plan: Consolidation and Revitalization at Technical Areas 6, 8, 9, 14, 15, 22, 36, 39, 40, 60, and 69, Vol. 1 (McGehee et al. 2005a). The SHPO concurred with these eligibility determinations on April 18, 2006. (4) TA-37-0006, TA-37-0009, and TA-37-0020 were determined eligible for listing in the NRHP in the report, High Explosives and the Nuclear Stockpile: An Assessment of Historic Buildings at Magazine Area C (TA-37), Vol. 1 (McGehee et al. 2008b). The SHPO concurred with these eligibility determinations on April 17, 2008. In a letter dated January 24, 2020, NA-LA acknowledged that the D&D of these nine NRHP-eligible facilities was an adverse effect that requires resolution through mitigation. NA-LA proposed the use of standard mitigation practices as defined in the Programmatic Agreement (PA) among the U.S. Department of Energy, National Nuclear Security Administration, Los Alamos Field Office, the New Mexico State Historic Preservation Office, and the Advisory Council on Historic Preservation Concerning Management of the Historic Properties at Los Alamos National Laboratory, Los Alamos, New Mexico. The PA states that adverse effects to NRHP-eligible buildings and structures will be resolved according to the standard practices defined in Part II, Section 10, of the Laboratory’s Cultural Resources Management Plan, A Plan for the Management of the Cultural Heritage at Los Alamos National Laboratory, New Mexico (Purtzer et al. 2019), and Section 2.B of Appendix D of the PA itself. The standard practice documentation package includes the following components: (1) Interior and exterior photography and production of archival-quality digital photographs; (2) Documentation and curation of historically significant equipment and artifacts; (3) A list of all known drawings for the property; (4) Reduced-scale reproductions of selected drawings for the property; (5) A location map that shows the location of the property relative to the entire Laboratory property; (6) Reproduction of historical TA maps; (7) A TA map that depicts the footprint of each eligible and non-eligible facility; and (8) An expanded historic context that uses oral-history interviews, if available. On March 3, 2020, the SHPO concurred with the adverse effect determination and the mitigation plan. The documentation package, as previously described, is provided in Volumes 1 and 2 of this report.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

Full Moment Tensor Solutions of U.S. Underground Nuclear Tests for Event Screening and Yield Estimation

Moment tensor (MT) solutions are proving increasingly valuable in explosion monitoring, especially now that they are more routinely calculated for the unconstrained, full (six component) MT. In this study, we have calculated MTs for U.S. underground nuclear tests conducted at the Nevada National Security Site using seismic recordings primarily from the Livermore Nevada Network. We are able to determine them for 130 nuclear explosions from 1970 to 1996 for a range of yields and under a variety of material conditions, which we have supplemented with 10 additional chemical explosions at the test site. The result is an extensive database of MTs that can be used to assess the performance of important monitoring tasks such as event identification and yield determination. We test the explosion event screening on the fundamental lune of the MT eigensphere and find MT screening to be a robust discriminant between earthquakes and explosions. We then explore the estimation of moment-derived yield, in which we find that material properties are the largest contributor to differences in the estimated moment-to-yield ratio. Additionally, further research conducted on this dataset can be used to develop, test, and improve various explosion monitoring methodologies.

58 GEOSCIENCES↗
DOE OSTI2025

Shear-Wave Radiation Patterns from Explosive and Earthquake Sources in Scattering, Heterogeneous Media

Distinguishing whether a seismic event is an earthquake or an explosion is a core problem in explosion monitoring. For simplistic models, earthquakes produce a predictable S-wavefield based on their radiation pattern, while an explosion produces no S-wavefield. However, observations from nuclear tests show that explosions can produce significant S-wave energy on both horizontal components. We perform numerical experiments using SPECFEM to constrain when an S-wavefield generated from an explosion might differ from one produced by an earthquake. We generate the S-wavefield for our explosions by placing the source location within a small region consisting of large velocity heterogeneities. For the 2D cases investigated here, a crucial condition is that the size of the region of heterogeneities must be much smaller (approximately 1/12th) than the minimum wavelength of the wavefield for our simulation setups. We investigate both isotropic heterogeneities and anisotropic heterogeneities. Initial results demonstrate that the heterogeneous region lowers the P/S amplitude ratios of an explosion. The presence of strong anisotropic heterogeneities near the source can also produce an S-wavefield with similar amplitudes to the P-wavefield. Lastly, we show that S-wave polarization angles vary as a function of azimuth and distance for explosive sources in the presence of small-scale heterogeneities.

58 GEOSCIENCES↗
DOE OSTI2023

The Observed Inefficiency of Explosions to Produce Large Aftershocks: Båth’s Law for Explosions is 2.5

Underground explosions are observed to produce fewer and smaller aftershocks than similar size earthquakes. The seismic magnitude difference $Δm_x$ between an explosion and its largest aftershock is an expression of Båth’s law for explosions. Based on an analysis of a compilation of aftershock studies from Soviet testing at the Semipalatinsk test site in Kazakhstan and observations from American testing at the Nevada National Security Site (NNSS), here we find that the average magnitude difference for explosions $\overline{Δm_x}$ is about 2.5. Based on the NNSS data, two standard deviations of $Δm_x$ is about 1.5. In all the cases studied, from ton to megaton yield, from shallow to overburied depth, and chemical or nuclear source, no explosion aftershock has been larger than the explosion that preceded it. In fact, the two events at the NNSS with the largest aftershock magnitudes relative to the explosion are associated with the collapse of the cavity created by the explosion. This is similar to observations from North Korean testing at the Punggye-ri Test Site, where the largest seismic event following the test is attributed to the collapse after the 2017 explosion and is from 0.8 to 2 magnitude units less than the mainshock.

58 GEOSCIENCES↗