Engineering topics

Walter, William R.

Publications and source records attributed to Walter, William R..

At least 19 records

DOE OSTI2024

Evaluating Physics-Informed Neural Network Performance for Seismic Discrimination between Earthquakes and Explosions

In this article, we evaluate adding a weak physics constraint, that is, a physics‐based empirical relationship, to the loss function with a physics‐informed manner in local distance explosion discrimination in the hope of improving the generalization capability of the machine learning (ML) model. We compare the proposed model with the two‐branch model we previously developed, as well as with a pure data‐driven model. Unexpectedly, the proposed model did not consistently outperform the pure data‐driven model. By varying the level of inconsistency in the training data, we find this approach is modulated by the strength of the physics relationship. In conclusion, this result has important implications for how to best incorporate physical constraints in ML models.

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↗

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 OSTI2022

The Source Physics Experiment (SPE) Science Plan

The Source Physics Experiment (SPE) series is a long-term NNSA research and development effort designed to improve U.S. arms control and nuclear nonproliferation verification and monitoring capabilities. The findings from the SPE will advance the United States’ nuclear explosion monitoring capabilities, particularly with respect to detection, discrimination and determination of yields associated with small nuclear explosions that can be lost amid the noisy seismo-acoustic background from other sources. The data generated from the SPE, a series of well-designed and recorded chemical explosions, will contribute to the development and validation of first-principles explosive source generated seismo-acoustic modeling codes. These codes will then facilitate the update of semi-empirical methods, currently based on historic test site data, such that key explosion observables can be reproduced, thus improving confidence in nuclear test monitoring in new areas and/or under novel emplacement conditions. The overall SPE project is comprised of both the development of the new explosion simulation codes and the chemical explosion test series. The chemical explosion test series will generate the empirical data required to both develop and validate the new simulation codes.

58 GEOSCIENCES↗

DOE OSTI2022

Analysis of the Source Physics Experiment Dry Alluvium Geology Ball Drop Sources

Two balls were dropped approximately 150 m down an open borehole as part of an investigation into the aftershocks from a 51 ton TNT-equivalent chemical explosion. The ball drop sources, a concrete-filled basketball and a bowling ball, were intended to mimic falling debris within an underground explosion borehole and cavity. Falling debris or sloughing of the borehole and cavity walls is thought to be one of the sources of explosion aftershocks. The signals from the ball drop sources were recorded by the same seismic network that monitored the aftershocks of the chemical explosion, so the ball drop recordings were used as templates to search through the continuous stream of aftershocks and detect correlated signals as a proportion of all aftershocks. The search resulted in only four similar event recordings out of the thousands of detected aftershocks. A synthetic test showed that the small, high-frequency ball drop templates were limited to producing highly correlated detections for events occurring in a small region around the impact location at 150 m depth. The high-frequency templates were necessary since the small ball drop signals were only above the noise at frequencies greater than 10 Hz. Unfortunately, the impact location is where the borehole was filled with sand in preparation for the subsequent chemical explosion in the test series, making it 150 m above the explosion created cavity centered at 300 m depth. This puts the ball drop templates 150 m away from the cavity, where many more events from falling debris would be expected. Further synthetic tests confirmed events emanating from the cavity 150 m below the impact site would have correlated with the lower frequency template and been detected had the template been usable below 10 Hz. Therefore, we conclude that the four detected events were most probably falling debris within the open borehole that impacted the floor at 150 m depth just as the ball drop sources had, and that this type of explosion aftershock mechanism comprises less than 1% of all explosion aftershocks.

58 GEOSCIENCES↗

DOE OSTI2022

Combining Deep Learning With Physics Based Features in Explosion-Earthquake Discrimination

This paper combines the power of deep-learning with the generalizability of physics-based features, to present an advanced method for seismic discrimination between earthquakes and explosions. The proposed method contains two branches: a deep learning branch operating directly on seismic waveforms or spectrograms, and a second branch operating on physics-based parametric features. These features are high-frequency P/S amplitude ratios and the difference between local magnitude (ML) and coda duration magnitude (MC). The combination achieves better generalization performance when applied to new regions than models that are developed solely with deep learning. Further, we also examined which parts of the waveform data dominate deep learning decisions (i.e., via Grad-CAM). Such visualization provides a window into the black-box nature of the machine-learning models and offers new insight into how the deep learning derived models use data to make decisions.

58 GEOSCIENCES↗

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 OSTI2021

A Big Problem for Small Earthquakes: Benchmarking Routine Magnitudes and Conversion Relationships with Coda Envelope-Derived Mw in Southern Kansas and Northern Oklahoma

Earthquake magnitudes are widely relied upon measures of earthquake size. Although moment magnitude (M w ) has become the established standard for moderate and large earthquakes, difficulty in reliably measuring seismic moments for small (generally M w <4) earthquakes has meant that magnitudes for these events remain plagued by a patchwork of inconsistent measurement scales. Because of this, magnitudes of small earthquakes and statistics derived from them can be biased. Furthermore, because small earthquakes are much more numerous than large ones, many applications, such as seismic hazard modeling, depend critically on analysis of events characterized by magnitudes other than M w . Therefore, to assess this problem, we apply coda envelope analysis to reliably determine moment magnitudes for a case study of small earthquakes from northern Oklahoma and southern Kansas. Not surprisingly, we find significant differences among M L , m bLg , and M w for M ~2–4 earthquakes examined here. More troublingly, we find that relations designed to convert other magnitudes to M w , which are relied upon for important applications such as seismic hazard analysis, often increase rather than decrease this bias for our dataset. In our case study, we find that converted magnitudes can result in a systematic bias sometimes exceeding 0.5 magnitude units, a difference that typically corresponds to a factor of ~3 in seismicity rate. Moreover, we find a correspondingly large bias in Gutenberg–Richter b-values, controlled primarily by inaccurate magnitude scaling in the conversion relationships. Although this study focuses on a relatively small geographic area, we can expect that similar issues exist with varying severity in other regions. Therefore, magnitudes of small earthquakes and their associated statistics, including seismicity rates and b-values, should be treated with caution.

58 GEOSCIENCES↗

DOE OSTI2021

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

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

58 GEOSCIENCES↗

DOE OSTI2021

Toward Robust and Routine Determination of Mw for Small Earthquakes: Application to the 2020 Mw 5.7 Magna, Utah, Seismic Sequence

To better characterize seismic hazard, particularly, for induced seismicity, there is an increasing interest in methods to estimate moment magnitude (M w ) for small earthquakes. M w is generally preferred over other magnitude types, but, it is difficult to estimate M w for earthquakes with local magnitude (M L ) <3-3.5, using conventional moment tensor (MT) inversion. The 2020 M ww 5.7 Magna, Utah, seismic sequence provides an opportunity to illustrate and evaluate the value of spectral methods for this purpose. Starting with a high-quality seismic catalog of 2103 earthquakes (M L <5.6), we estimate M w using two independent spectral methods—one based on direct waves, yielding M w,direct , and the other based on coda waves, yielding M w,coda . For the direct-wave method, we present a non-parametric (NP) inversion scheme that solves for apparent geometrical spreading, G(R), and site effects (S), similar to other NP procedures that have been used to calibrate regional M L scales. The NP inversion is constrained using M ws derived from MTs for nine events in the Magna sequence. We recover statistically robust and physically reasonable G(R) and S and compute M w,direct for 635 Magna earthquakes down to M L 0.7. For the coda-wave method, we consider two separate calibration schemes involving previous MT solutions and compute M w,coda for 311 earthquakes down to M L 1.0. For 280 of the events that were processed with both methods—M w,direct and M w,coda —are strongly correlated (r = 0.98), with a mean difference of only 0.05. We compare M w,direct and M w,coda with M L and find reasonably good agreement for M L <3.6 with the theoretically predicted relationship of M w =(2/3)M L +C, in which C is a regional constant. Our results imply that seismic network operators can use spectral-based M w estimates to replace M L estimates for events with M L ≥1.0, and possibly smaller. The main requirement is the existence of a small number of MT solutions for calibration purposes.

58 GEOSCIENCES↗

DOE OSTI2020

New Opportunities to Study Earthquake Precursors

The topic of earthquake prediction has a long history, littered with failed attempts. Part of the challenge is that possible precursory signals are usually reported after the event, and the systematic relationships between potential precursors and main events, should they exist, are unclear. Furthermore, several recent studies have shown the potential of new approaches to simultaneously detect earthquake foreshocks and slow-slip phenomena through ground deformation, seismic, and gravitational transients—weeks to months before large subduction zone earthquakes. The entire international community of earthquake researchers should be engaged in deploying instrumentation, sharing data in real time, and improving physical models to resolve the extent to which slow slip events and earthquake swarms enhance the likelihood (or not) for later, larger earthquakes.

58 GEOSCIENCES↗

DOE OSTI2020

Rescuing Legacy Seismic Data FAIR’ly

The Earth’s interconnected and dynamical systems operate on a spectrum of time scales from millions of years to fractions of a second. While the evolution of the Earth and its deep interior are beyond the time span of human observations, understanding of many natural phenomena operating on human times scales have benefited from direct scientific observation. Continuous processes and those that are repeated over time shape the environment we live in. Furthermore, as we are faced with unprecedented changes to climate, understanding the deeper patterns and trends in natural systems through time have taken on new importance (Research Data Alliance, 2019). The call to reuse data is driven not only by economics but also by the recognition of their complete uniqueness (observations of natural systems are not repeatable) and scientific value in enhancing current understandings as well as potential new discoveries especially in the era of big data. These data are part of the historical record and our scientific heritage (American Geophysical Union, 2019) not only in explicitly recording earth observations but implicitly recording, and thus providing the evidence that addresses the manner in which science was conducted.

58 GEOSCIENCES↗

DOE OSTI2019

The Source Physics Experiment (SPE) Science Plan (Version 2.4c)

58 GEOSCIENCES↗

DOE OSTI2019

Investigating the Effectiveness of P/S Amplitude Ratios for Local Distance Event Discrimination

The ability to discriminate between underground explosions and earthquakes is an important operation for global security and test-ban-treaty monitoring, and seismology plays a key role in this endeavor. The use of seismic P- to S-wave amplitude ratios has been demonstrated to be an effective discriminant at regional distances for a variety of emplacement conditions and test sites around the world. Increasing interest in detecting and discriminating smaller events, both for purposes of nuclear monitoring and seismic hazard estimation, requires the use of data with shorter event-to-station paths. At these local distances (≤ 200 km), the use of P/S ratios for event discrimination is poorly studied compared with regional distances (200-1500 km). The path averaging that allows for relatively simple corrections to P/S ratios to remove bias from along-path structural differences at regional distances may be insufficient to remove enough scatter from the ratios to allow for accurate event identification at local distances. We use data from a series of chemical explosions detonated as part of the Source Physics Experiment at the Nevada National Security Site in southern Nevada and surrounding seismicity to investigate the discrimination performance of P/S ratios at local distances. We find that discrimination at individual stations is highly variable, regardless of frequency. Network averaging of ratios improves discrimination performance above 6 Hz, as a result of the reduction of scatter among the earthquake ratios. Furthermore, azimuthal and distance–limited averaging suggests that explosion ratios are affected by a complicated source and path relationship that may lead to substantial variability in local P/S ratios at different locations.

58 GEOSCIENCES↗

DOE OSTI2018

Preface to the Focus Section on North Korea’s September 2017 Nuclear Test and Its Aftermath

On September 3, 2017 a large seismic event was detected near the Punggye-ri nuclear test site in North Korea (officially known as the Democratic People Republic of Korea or DPRK) by several geophysical monitoring institutions (for example, the USGS announced it as mb 6.3 event, the CTBTO announced it as a mb 6.1 event). Within a few hours North Korea declared it had carried out its sixth underground nuclear explosion, a test of a hydrogen bomb.

58 GEOSCIENCES↗

DOE OSTI2018

Body-Wave Methods of Distinguishing between Explosions, Collapses, and Earthquakes: Application to Recent Events in North Korea

Examination of regional distance seismic data from historic nuclear test sites has led to a variety of very effective discriminants between explosions, earthquakes, and collapses. We focus on the body–wave methods. We show that ratios between P– and S–wave amplitudes (P/S ratios) above about ~2 Hz very effectively separate the six Democratic People’s Republic of Korea (DPRK) declared nuclear tests between 2006 and 2017 from natural earthquakes in the region. Similarly, P/S ratios separate historic Nevada Test Site (NTS) nuclear explosions from western U.S. earthquakes. We show that combining P/S ratios with ratios of low–frequency to high–frequency S–wave amplitudes can effectively separate postexplosion collapse events, such as the 1982 NTS Atrisco collapse, and the apparent collapse that followed about eight and a half minutes after the 3 September 2017 DPRK explosion. Explosions often produce fewer and smaller aftershocks than comparably sized earthquakes, which has been proposed as a potential discriminant. We apply the body–wave techniques to the recent seismicity following the largest DPRK event, after first using correlation methods to build a more complete catalog of these events. Despite the empirical effectiveness of the regional body–wave discriminants, the physical basis for the generation of explosion S waves, and therefore the predictability of P/S and low/high frequency techniques, as a function of path, frequency, and event properties such as size, depth, and geology, remains incompletely understood. Furthermore, a goal of current research, such as the Source Physics Experiments (SPE), is to improve our physical understanding of the mechanisms of explosion S–wave generation and advance our ability to numerically model and predict them.

58 GEOSCIENCES↗

DOE OSTI2018

A Preliminary Explosion Seismic Spectral Model for Saturated/Hard Rock

This short technical report describes a new seismic spectral model for chemical and nuclear explosions in saturated and/or hard rock. Results from the Source Physics Experiments (SPE) indicate that existing explosion spectral models can be poor fits to small and/or over-buried explosions. Here we develop a preliminary improved explosion model to fit both SPE and historic nuclear explosion data. The model is based on revising the Mueller and Murphy (1972) and Denny and Johnson (1991) seismic spectral models. A FORTRAN code was developed to provide the spectral models in SAC format for use in explosion monitoring applications such as detection and yield estimation. This model should be considered preliminary and work is underway to better generalize and expand the model to include explosions in dry porous and/or soft rock media and to better fit explosion S-waves.

58 GEOSCIENCES↗

DOE OSTI2017

High-Resolution 2D Lg and Pg Attenuation Models in the Basin and Range Region with Implications for Frequency-Dependent Q

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.

58 GEOSCIENCES↗