Search NASA⌕ Search

Engineering topics

Phillips, William Scott

Publications and source records attributed to Phillips, William Scott.

A Regional Phase Amplitude Model of 2-D Attenuation for North America

We analyzed seismic attenuation patterns across the North America using over 70,000,000 Lg wave amplitudes recorded by the various networks at frequencies of 0.1-32 Hz. Our inversion solved for laterally varying attenuation, site terms, moments, and apparent stress following Phillips et al., (2016). The inversion was anchored by independently constrained: corner frequencies (via coda spectral ratios) to control the attenuation-stress tradeoff and moment measurements of teleseismic (GCMT, USGS) and regional (St. Louis University and UC Berkeley) earthquakes to provided absolute scaling. The quality factor (Q) shows clear regional patterns: low values in coastal, volcanic, and tectonically active regions, and high values in stable areas like the Great Plains and major plateaus throughout North America. These 2-D Q models enable improved regional source characterization, magnitude estimation, and yield determination.

58 GEOSCIENCES↗

Probabilistic Programming for Transportable Source Characterization and Uncertainty Quantification of the North Korean Nuclear Tests 2006–2017

Here, we introduce a transportable technique to determine the yield and depth of burial (DOB) from seismic source spectra of underground nuclear explosions. We demonstrate this technique on the six declared North Korean nuclear tests. This approach derives source spectra in absolute units from regional phase (Pg) amplitudes by correcting the observations for geometric spreading, attenuation, and site amplification. We couple the source spectra and explosion source models with a probabilistic programming framework that integrates deep learning techniques and Bayesian modeling. This approach permits the exchange of information across various data categories to quantify both the data and model uncertainty. This technique stands out as an innovative use of broad‐area propagation models, making it transportable across various geologic settings. This method proves to be effective in scenarios with diverse and/or limited observational data, even when the source depth is unknown. We present new independent estimates of absolute yield and DOB that are consistent with the prior assessments, underscoring the potential of this method in enhancing transportable nuclear explosion monitoring capabilities.

58 GEOSCIENCES↗

Transportable absolute yields of underground nuclear explosions. Application to the North Korean nuclear tests

The method developed and presented here improves our nuclear explosion monitoring capability by providing a new independent and transportable method to determine the explosive yields of seismically recorded events associated with foreign weapons development. Here we report new independent absolute yield and depth of burial estimates of the North Korean nuclear tests which are commensurate with previously determined source parameters. The key advantage of this method over other methods which estimate absolute explosive yields is that this method does not require any a priori calibration, and can be immediately applied to any region of interest. This method uses the source information obtained from the spectral ratios of envelopes of measured seismic coda waves to invert for model-derived yield and depth of burial. The only requirements of this method are that the region of interest must contain several (> 3) seismic sources which have been well recorded (SNR > 2) by a set of shared stations. For regions with only a few events (~ 3) and/or extremely band-limited observations for which the depth of burial is unknown, an earthquake with a known magnitude may be required to obtain an accurate and robust yield estimate.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Identifying Suspect Instrument Intervals Using Midnight Noise Time Histories

With an intense work ethic, and high levels of devotion to their craft, seismologists continue to battle the elements and wrestle with complex logistics in their efforts to field instrumentation over the entire globe. We are fortunate to see data volumes accelerating in size, and quality continuing to improve; for example, universal timing issues can now be considered rare. However, instrument response issues are still common, and can hinder research that seeks to understand amplitudes of seismic wavefields. We are interested in using the accumulations of seismic data to develop models that will predict high frequency (0.2-20 Hz, and higher) signal amplitudes for explosion monitoring purposes over broad areas, which requires large data sets and extensive quality control. To identify response and station health issues, we have collected noise time histories for global seismic data, focusing on measurements near (but not restricted to) midnight to eliminate diurnal variations, and have manually determined time intervals that appear inconsistent with background behavior. We assign descriptive labels, but do not attempt to diagnose causes. We use these intervals to discard data. To date, we have examined 39,260 channels from 11,105 stations, heavily weighted toward IRIS holdings, dates through 2017 (depending on station, roughly the date we started our manual review), bands between 1 and 8 Hz, finding 24,733 anomalous time intervals. The great majority (90%) of these intervals appear to be shifts of constant offset, often bounded by times of known instrument changes, likely the result of poor documentation of response parameters at one of many stages between the field and the plotter. We hope these results can be of use to our colleagues, and would encourage community efforts to diagnose anomalous behavior, and fix poor responses. We also hope these results will support automation efforts, including application of supervised learning techniques.

47 OTHER INSTRUMENTATION↗

Explosion Source Models and the Scattering Origin of Regional Phases from SPE Phase 1 Coda Spectral Ratios: Monitoring Implications –Discrimination and Yield Estimation

The Source Physics Experiment (SPE) explosions provide comprehensive ground truth information, and near-source recordings that are rarely available in typical monitoring scenarios. Here, we analyze these data together with local and regional distance recordings to isolate source and propagation effects, and to potentially improve monitoring at distance. We examine source spectral ratios, which can be measured to high precision by taking advantage of many SPE and independently operated stations, as well as the redundancy available from coda waves. Classical source models (Mueller and Murphy, 1971, MM71; Denny and Johnson, 1991, DJ91) predict the source ratios poorly; however, a hybrid model (MM71 with DJ91 cavity radii) performs better. We observe a distinct spectral modulation at 6-9 Hz that is not predicted by classical models, most likely caused by short period surface waves (Rg) interfering with those produced by spallation of near surface layers. The same modulation is observed for compressional (P) and shear (S) waves at distance, indicating that local and regional phases originate as near-source Rg that is scattered into body waves. The scattering process must be accounted for quantitatively to best use remotely recorded signals to monitor treaty compliance, discriminate event types, and estimate explosion yields over broad regions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Traditional Look at SPE Phase 1 Surface Seismic Data: Far-Field Instruments, Measurements, and Analysis

We visually examined surface seismic data from the Source Physics Experiment (SPE) Phase-1 chemical explosions in preparation for waveform analysis and interpretation. Recording stations included far-field (near source) geophone lines, and both permanent network, reoccupied USArray, and temporary stations at local and regional distances. We found that near-source waveforms are dominated by broad band compressional (P), and band limited surface wave (Rg) (< 10 Hz) phases. Distinct, consistent shear (S) phases are difficult to observe visually even with a comprehensive bank of frequency bands. Secondary phases are indistinct at near-local distances, unlike for earthquake sources, especially at the high frequencies of interest for discrimination (e.g. 4-8 Hz), likely due to intense near surface scattering. S phases begin to emerge from the scattered background at and beyond 100 km. We identified regional Pn, Pg and Lg phases to 446 km for SPE-5. The data set is, unfortunately, plagued by poorly recorded instrument responses, or, possibly, time varying site amplification due to human interaction with sensors. Instrument response issues should be fixable, and this should be attempted to make this a truly, high-quality data set for use by the monitoring community.

58 GEOSCIENCES↗

Explosion Source Models and the Scattering Origin of Regional Phases from SPE Phase 1 Coda Spectral Ratios

The Source Physics Experiment(SPE)explosions provide comprehensive ground truth information,and nearby recordings that are rarely available in typical monitoring scenarios. Here, we look for relationships with local and regional distance data to understand source and propagation effects in an effort to improve monitoring at a distance. We examine source spectral ratios, which can be measured to high precision by taking advantage of many SPE and independently operated stations, as well as the redundancy available from coda waves. Classical source models (Mueller and Murphy, 1971, MM71; Denny and Johnson, 1991, DJ91)predict the source ratios poorly; however, a hybrid model(MM71 with DJ91 cavity radii)performs better, and by varying cavity and elastic radii we obtain models that fit to within measurement precision. Based on this, the best prediction model for the SPE suite is the hybrid model with cavity radius decreased by 95%, and elastic radius decreased by 75%. This works best for mid-range scaled depths of burial (SPE-1, 2, 3, and 5). We do not guarantee transportability to other geologies,depth and scaled depth ranges. Finally, we observe a distinct spectral modulation at 6-9Hzthat is not predicted by classical models, and is most likely the result of interference of short period surface waves(Rg)with those produced by spallation of near surface layers.The same modulationis observed for compressional (P) and shear(S) waves at distance, indicating the origin of these waves as scattered Rg.The scattering process must be accounted for quantitatively when using distant signals to monitor treaty compliance over broad regions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Identifying and Correcting Timing Errors at Seismic Stations in and around Iran

A fundamental component of seismic research is the use of phase arrival times, which are central to event location, Earth model development, and phase identification, as well as derived products. Hence, the accuracy of arrival times is crucial. However, errors in the timing of seismic waveforms and the arrival times based on them may go unidentified by the end user, particularly when seismic data are shared between different organizations. Here, we present a method used to analyze travel-time residuals for stations in and around Iran to identify time periods that are likely to contain station timing problems. For the 14 stations with the strongest evidence of timing errors lasting one month or longer, timing corrections are proposed to address the problematic time periods. Finally, two additional stations are identified with incorrect locations in the International Registry of Seismograph Stations, and one is found to have erroneously reported arrival times in 2011.

58 GEOSCIENCES↗

Pn Geometric Spreading

Pn geometric-spreading model that accounts for the sphericity of the Earth.

97 MATHEMATICS AND COMPUTING↗