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.
Multi-Phenomenology Explosion Monitoring (MultiPEM) Toolbox provides a capability for post-detonation characterization focusing on yield estimation; Synthesize information from multiple sensor types; Leverage historical databases to inform statistical models; Allows for rapid or comprehensive analysis; Implemented in R code.
Distributed Acoustic Sensing (DAS) is an emerging technology capable of recording the acoustic wavefield at unprecedented spatial resolution. However, this new tool requires significant refinements before it becomes operational for explosion monitoring objectives. Recent studies have shown significant development of array processing with DAS data. In this report we explore three such array processing methods including DAS strain-rate data versus geophone measured ground motion, beamforming for event parameters, and machine learning based denoising. Prior to applying these algorithms to the Source Physics Experiment Phase II and Phase III data, we validate these methods through replication analysis.
The accurate estimation of seismic wave amplitude is vital to precisely determine the yield, magnitude, and event discrimination possible for a given network – a critical element in nuclear explosion monitoring. This task is complicated by several factors, including but not limited to radiation pattern, scattering effects, and crustal variations, which can lead to the attenuation or amplification of amplitude along a given raypath. In this report, we explore the novel application of machine learning to the task of seismic amplitude estimation by training a simple Artificial Neural Network (ANN) on an S-wave amplitude dataset from Lai et al. (2019). Attributes from this dataset used as input to the ANN included event-station distances, station locations (latitude, longitude), event locations (latitude, longitude), event depths, event magnitudes, radiation patterns, signal-to noise ratio (SNR) measurements (average-amplitude, peak-to-trough, maximum peak), and signal periods. We find that the trained ANN predicts S-wave amplitudes with a modest tendency toward underestimating the actual values, as indicated by a linear regression between predicted and actual data (slope: 0.892, intercept: -0.651). These results suggest that an ANN can perform this task, with potential for significant improvements through improved datasets, architectures, and parameter tuning.
Co-Author Voon Hui Lai, Australian National University, will be attending IUGG and presenting the work. The session of interested is "S05 - Advances in Earthquake and Explosion Monitoring Using Distributed Acoustic Sensing"
This document explains use of the Multi-Phenomenology Explosion Monitoring (Multi PEM) Toolbox, a collection of R scripts for estimating the unknown device parameters of a new event with uncertainty quantification. The methodology and application used for il lustration in this user manual are published in Geophysical Journal International and fully documented in a Los Alamos National Laboratory technical report hereafter designated “WPA” for reference. Additional details on the application are found in a recent journal article. Two assessment types are available: rapid and complete.
To be presented at the 2021 NEM. The March 15-19, 2021 virtual Nuclear Explosion Monitoring (NEM) Program Review combines reviews of several related programs within the Office of Defense Nuclear Nonproliferation Research and Development (DNN R&D) of the National Nuclear Security Administration.
The INL Noble Gas Laboratory provides intercomparison samples for the noble gas analysis laboratories as part of the CTBTO PrepCom IMS. Xe-133m is one of the four relevant radionuclides in nuclear explosion monitoring. Without commercially available Xe-133m calibration standards laboratories must create and improve calibration methods. Improvements in calibration methods at the INL NGL benefit the CTBTO PrepCom through better certified values for Xe-133m intercomparison samples. Calibration of High Purity Germanium detectors for Xe-133m quantification is complicated by the coexistence of Xe-133 in samples under analysis. Xe-133 is typically produced in larger quantities, has higher gamma emission probabilities, and its gammas are detected more efficiently than Xe-133m. Xe-133m activity of samples can be indirectly inferred through the 133:133m activity ratio of a batch of material, and the Xe-133 counts in the assay of a small aliquot of the same material. This indirect quantification method can be leveraged to perform detector calibrations for quantification of Xe-133m. Calibrations can be performed by inferring the Xe-133m to certify the sample, and direct counting to determine detector efficiency. A comparison of method results will be shown.
The Comprehensive Nuclear-Test-Ban Treaty (CTBT) bans the testing of nuclear weapons anywhere on the earth (atmospheric, surface, underwater and subsurface). Identification of nuclear explosions in the atmosphere, surface, and underwater is relatively straightforward considering a wide range of signatures resulting from such an event. However, for a subsurface explosion, most of the signatures traditionally associated with a nuclear explosion are not readily available. Therefore, the international community has increasingly relied on the atmospheric measurement of noble gases to identify subsurface nuclear weapon explosions. This chapter initially covers the basic principles of subsurface nuclear explosion identification and the importance of detecting radioxenon. This is followed by reviewing some of the early radioxenon detection systems that were developed by research groups around the world in the late 1990s and early 2000s. The detection media employed, results from laboratory and field testing, and some challenges/drawbacks for these systems are detailed. The next section of the chapter is dedicated to innovative detector concepts that have emerged in the past ten to fifteen years using novel detection material, algorithms, and signal readout techniques. The advances achieved in terms of energy resolution, coincidence detection efficiencies, system performance, and the minimum detectable concentration are covered. The final section goes over some of the potential improvements that can be incorporated in the design to enhance detector sensitivity and new detection material that can be explored in the field of radioxenon detection.
The noble gas collection and measurement stations in the International Monitoring System (IMS) are heavily influenced by releases from medical isotope production facilities. The ability to reliably model the movement of radioxenon from the points of release to these IMS samplers has improved enough that a routine aspect of the analysis of IMS radioxenon data should be the prediction of the effect of releases from civilian nuclear facilities on the sample concentrations. In this work, predicted concentrations at IMS noble gas systems in Germany and Sweden based on measured releases from Institute for Radioelements (IRE) in Belgium and atmospheric transport modeling for a four-month period are presented and discussed.
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.
The ability to differentiate between atmospheric radionuclide signatures from underground nuclear explosions (UNEs) and signals from other sources, such as medical isotope-production facilities and nuclear reactors, can be critical to the detection and monitoring of unannounced, low-yield nuclear events. Signatures having anomalously high amplitudes, compared to background levels, remain the best indicator in screening for a UNE. However, isotopic composition can further validate a suspected UNE signature, but separation from any atmospheric background composition is first necessary. To date, evaluating the challenges of performing this separation has typically involved comparing an observed background with a highly idealized deterministic model of radioxenon signature production by a UNE that does not consider the influence of post-detonation chemical/physical processes in the detonation cavity or the subsequent gas transport mechanisms that can also affect the isotopic composition of the detected gas signature. In addition, purely deterministic models, as previously employed, overlook the uncertainty inherent in estimating critical parameters characterizing the UNE and its detonation environment. In this paper, we create detailed, multi-parameter models of radionuclide evolution using the widely accepted England and Rider post-detonation radionuclide decay-chain network coupled to detailed models simulating physical production and transport processes affecting the gas signature. Because these models are governed by uncertain parameters including barometric fluctuations, realistic ranges of variation for each of the parameters influencing isotopic composition are then defined. A Latin-Hypercube sampling approach is used to obtain a random distribution of isotopic production and gas transport results associated with a given value of each parameter. We apply these results to background histories of two stations, one providing 4-isotope background measurements and the other providing two-isotope measurements associated with the 2013 DPRK announced UNE.
As part of a scientific research and development project, the radionuclide fingerprint of an operating advanced gas-cooled reactor 25 (AGR) has been studied across several facets (Goodwin et al., 2024). One part of this project was to deploy an array of 26 radioxenon samplers to the region for a period of around 1 year, to measure any radioxenon emissions from the reactors of the 27 Hartlepool nuclear power station at a range of tens of kilometres away. The array of 3 sensors was operational for around 12 28 months from March 2022 and detected many occurrences of isotopes of radioxenon. Here we provide a detailed analysis and 29 interpretation of the data and where possible, attribution of detections to a source or region. Whilst a key part of this work is to 30 measure any emissions of radioxenon from Hartlepool, this work presents one of the most comprehensive efforts to determine the 31 source of a great number of (mostly) 133Xe detections. A combination of different types of atmospheric dispersion modelling 32 techniques, including the use of stack monitoring data from nearby civil radioxenon-emitting nuclear facilities has enabled the 33 majority of detections to be attributed to one or more possible sources. Whilst emissions from Hartlepool have been detected on 34 the systems, the majority of detections are associated with a medical isotope production facility in Fleurus, Belgium (IRE).
Current underground nuclear explosion (UNE) detection strategies rely heavily on atmospheric noble gas sampling of radioxenon. However, discriminating nuclear weapons testing programs from civilian sources is difficult due to highly variable atmospheric radioxenon backgrounds and processes affecting subsurface transport of parent radionuclides. Here, we aim to study the transport of gases produced by subsurface explosions as novel stable signatures for underground nuclear explosion (UNE) monitoring. These gases may be produced in large quantities with distinct molecular ratios, which will be impacted by subsurface transport processes. To demonstrate how ratios of gases produced by explosions can change during transport in geomaterials, we conducted laboratory benchtop experiments on the transport of carbon dioxide (CO 2 ) and hydrogen (H 2 ) gases through variably saturated zeolitic tuff, which is abundant at the historic US testing site. We observed that zeolitic tuff sorbs substantial quantities of CO 2 while allowing H 2 to transport more freely, leading to changes in the molecular ratios of the two gases along the transport pathway. Gas uptake in the dry zeolitic tuff core was 72.3% for CO 2 , compared with 53.4% for xenon and 7.6% for H 2 . The presence of 20% water saturation disrupted the CO 2 sorption process, though to a lesser extent than observed for noble gases, with a 36.7% drop in xenon sorption compared with a 21.9% drop for CO 2 . These results represent the first observations of zeolite sorption altering explosive gas ratios during transport through geomedia relevant to nuclear proliferation monitoring.
The PIPS–SrI 2 (Eu) is a prototype atmospheric radioxenon detection system designed at Oregon State University in support of international efforts towards monitoring clandestine nuclear weapon testing activities. This detector aims to address some shortcomings found in currently deployed beta–gamma atmospheric radioxenon detection systems, such as lackluster energy resolution and memory effect, by employing modern detection materials and readout. The system uses a PIPSBox, a silicon-based gas cell, for electron detection, and a pair of ultrabright, D-shaped SrI 2 (Eu) scintillators coupled to silicon photomultipliers for photon detection. A custom eight-channel digital pulse processor equipped with a field programmable gate-array (FPGA) identifies electron–photon coincidences between the volumes in near real-time. Gas samples of the four radioxenon isotopes of interest were independently measured with the PIPS–SrI 2 (Eu) detection system to determine energy resolution and efficiency. Application of FPGA-based coincidence discrimination in near real-time reduced the ambient background count rate by 95.85 ± 0.04%. Using parameters from the Xenon International gas processing unit and assuming a blank sample and zero memory effect the minimum detectable concentrations (MDCs) for the isotopes were calculated to be 0.12 ± 0.03, 0.27 ± 0.05, 0.15 ± 0.02, and 1.00 ± 0.08 mBq/m 3 air for 131m Xe, 133 Xe, 133m Xe, and 135 Xe, respectively. These MDC estimates compare well with other radioxenon detection systems employed in the International Monitoring System (IMS) and indicate that the PIPS–SrI 2 (Eu) is in compliance with the Comprehensive Nuclear Test-Ban-Treaty Organization (CTBTO) sensitivity requirement of ≤ 1 mBq/m 3 for 133 Xe.
Gamma-gamma coincidence spectrometers have the potential to significantly enhance detection sensitivity for ultra-trace radionuclide measurements. The implementation of these spectrometers, however, is limited by the complexity of acquisition hardware, data processing and quantification. This work reports development of a novel radionuclide quantification software for γ-γ coincidence measurements. For any radionuclide, the software parses the Evaluated Nuclear Structure Data File (ENSDF) database, recursively simulating all possible γ-γ coincidence signatures and their respective emission and detection probabilities. Implemented using Python programming language, the software employs several strategies to boost overall computational performance. Since coincidence-based spectrometers are of notable interest in monitoring compliance for the Comprehensive Nuclear-Test-Ban Treaty (CTBT), the software’s execution was tested for 84 CTBT-relevant radionuclides. To date, the software has been experimentally validated for 15 radionuclides using the Advanced Radionuclide Gamma spectrOmeter (ARGO) at Pacific Northwest National Laboratory, USA (PNNL). Notably, the software can be operated in convergence mode, whereby coincidence detection efficiency’s convergence behavior can help avoid unreliable radionuclide activity estimates. With growing number of coincidence spectrometers worldwide, this paper aims to assist the radiation metrology community in developing similar software for their system.
Gamma attenuation attributable to xenon and carbon dioxide carrier gas in radioactive xenon assay samples has been calculated using Monte Carlo n-Particle (MCNP) transport code for the 0.65 mL flame sealed ampoule (bean), and 15 mL and 50 mL Schlenk tube geometries. The attenuation was correlated as a % loss per standard mL of carrier gas. In this work, high Purity Ge gamma spectroscopy experiments were performed to measure against the MCNP results. The experimental results from the 133 Xe in xenon carrier gas in a bean geometry showed 0.97 ± 0.4% loss per standard mL xenon, which was in agreement with the MCNP predictions.
For this report a method was developed to measure trace noble gas element adsorption to the surfaces of geologic materials in the presence of a background gas that could potentially compete for surface adsorption sites. Adsorption of four noble gas elements (Ne, Ar, Kr, and Xe) at a concentration of 100 ppm in helium and nitrogen were measured on a sample of crushed tuff at 0, 15, 30, and 45°C. In addition, Ne, Ar, Kr, and Xe at 250 ppm and 500 ppm in nitrogen at 15°C were measured. Noble gas adsorption was found to increase with increasing atomic mass and decreasing temperature. It was also observed that the relative increase in noble gas element adsorption with decreasing temperature tends to increase with increasing atomic mass. Noble gas adsorption from nitrogen was approximately an order of magnitude greater than for the noble gases in helium. As the noble gas concentrations in nitrogen increased, adsorption increased in a slightly non-linear fashion which could be modeled using a Freundlich isotherm. For noble gas concentrations that were ≤100 ppm Henry's Law constant were calculated.
Fiber optic distributed acoustic sensors (DAS) are becoming a widely used tool for seismic sensing. Here we examine recordings of two subsurface chemical explosions, DAG-1 and DAG-3, each of which was about one metric ton (TNT equivalent), that were recorded from a helical fiber installed in two boreholes 80 m away from the source location. Several clear phases including the initial P wave, a weak S wave, and a surface reflected P wave are observed on the helical DAS data. We estimate a velocity model using arrival times measured from the fiber. The DAS waveform data were compared with colocated accelerometers at specific depths in both frequency and time domains. Furthermore, the spectra of the DAS data matched spectra estimated from the accelerometer records. Comparisons of observed waveform shape between the accelerometer records and the fiber measurements (strain-rate) show reasonable agreement except for the data near the event depth. The DAS data and the accelerometer agreed in relative amplitudes but we had difficulties in matching absolute amplitudes, possibly due to errors in metadata. Synthetic strain-rate waveforms were calculated using a 2D wavenumber algorithm and matched the waveform shape and relative amplitudes. In general, DAS is effective at recording strong ground motions at high spatial density. Comparison of the synthetic seismograms with observed data indicate that the waveforms are not consistent with a pure isotropic explosion source and that the observed S waves originate from very near the source region.