We have used Matlab and Google Earth to construct a prototype application for modeling the performance of local seismic networks for monitoring small, contained explosions. Published equations based on refraction experiments provide estimates of peak ground velocities as a function of event distance and charge weight. Matlab routines implement these relations to calculate the amplitudes across a network of stations from sources distributed over a geographic grid. The amplitudes are then compared to ambient noise levels at the stations, and scaled to determine the smallest yield that could be detected at each source location by a specified minimum number of stations. We use Google Earth as the primary user interface, both for positioning the stations of a hypothetical local network, and for displaying the resulting detection threshold contours.
To make use of some portions of the National Nuclear Security Administration (NNSA) Knowledge Base (KB) for which no current operational monitoring applications were available, Sandia National Laboratories have developed a set of prototype regional analysis tools (MatSeis, EventID Tool, CodaMag Tool, PhaseMatch Tool, Dendro Tool, Infra Tool, etc.), and we continue to maintain and improve these. Individually, these tools have proven effective in addressing specific monitoring tasks, but collectively their number and variety tend to overwhelm KB users, so we developed another application - the KB Navigator - to launch the tools and facilitate their use for real monitoring tasks. The KB Navigator is a flexible, extensible java application that includes a browser for KB data content, as well as support to launch any of the regional analysis tools. In this paper, we will discuss the latest versions of KB Navigator and the regional analysis tools, with special emphasis on the new overarching inter-tool communication methodology that we have developed to make the KB Navigator and the tools function together seamlessly. We use a peer-to-peer communication model, which allows any tool to communicate with any other. The messages themselves are passed as serialized XML, and the conversion from Java to XML (and vice versa) is done using Java Architecture for XML Binding (JAXB).
To improve the nuclear event monitoring capability of the U.S., the NNSA Ground-based Nuclear Explosion Monitoring Research & Engineering (GNEM R&E) program has been developing a collection of products known as the Knowledge Base (KB). Though much of the focus for the KB has been on the development of calibration data, we have also developed numerous software tools for various purposes. The Matlab-based MatSeis package and the associated suite of regional seismic analysis tools were developed to aid in the testing and evaluation of some Knowledge Base products for which existing applications were either not available or ill-suited. This presentation will provide brief overviews of MatSeis and each of the tools, emphasizing features added in the last year. MatSeis was begun in 1996 and is now a fairly mature product. It is a highly flexible seismic analysis package that provides interfaces to read data from either flatfiles or an Oracle database. All of the standard seismic analysis tasks are supported (e.g. filtering, 3 component rotation, phase picking, event location, magnitude calculation), as well as a variety of array processing algorithms (beaming, FK, coherency analysis, vespagrams). The simplicity of Matlab coding and the tremendous number of available functions make MatSeis/Matlab an ideal environment for developing new monitoring research tools (see the regional seismic analysis tools below). New MatSeis features include: addition of evid information to events in MatSeis, options to screen picks by author, input and output of origerr information, improved performance in reading flatfiles, improved speed in FK calculations, and significant improvements to Measure Tool (filtering, multiple phase display), Free Plot (filtering, phase display and alignment), Mag Tool (maximum likelihood options), and Infra Tool (improved calculation speed, display of an F statistic stream). Work on the regional seismic analysis tools (CodaMag, EventID, PhaseMatch, and Dendro) began in 1999 and the tools vary in their level of maturity. All rely on MatSeis to provide necessary data (waveforms, arrivals, origins, and travel time curves). CodaMag Tool implements magnitude calculation by scaling to fit the envelope shape of the coda for a selected phase type (Mayeda, 1993; Mayeda and Walter, 1996). New tool features include: calculation of a yield estimate based on the source spectrum, display of a filtered version of the seismogram based on the selected band, and the output of codamag data records for processed events. EventID Tool implements event discrimination using phase ratios of regional arrivals (Hartse et al., 1997; Walter et al., 1999). New features include: bandpass filtering of displayed waveforms, screening of reference events based on SNR, multivariate discriminants, use of libcgi to access correction surfaces, and the output of discrim{_}data records for processed events. PhaseMatch Tool implements match filtering to isolate surface waves (Herrin and Goforth, 1977). New features include: display of the signal's observed dispersion and an option to use a station-based dispersion surface. Dendro Tool implements agglomerative hierarchical clustering using dendrograms to identify similar events based on waveform correlation (Everitt, 1993). New features include: modifications to include arrival information within the tool, and the capability to automatically add/re-pick arrivals based on the picked arrivals for similar events.
The International Monitoring System (IMS) proposed for verifying compliance with the Comprehensive Nuclear-Test-Ban Treaty will include an infrasound network for detecting and identifying explosions in the atmosphere. As is the case with seismic monitoring, data collected from historic events of interest are vital for improving infrasonic monitoring capabilities. Unfortunately, however, infrasonic recordings of such events are rare and thus any additional data sets that might be available should be pursued. Towards that end, we will digitize, as a result of the ROA01-38 award, paper records and extract from 9-track tapes several unique data sets from Sandia National Laboratories and Los Alamos National Laboratory that have not been available to the monitoring community. These data sets include recordings of surface and atmospheric explosions representing different yields, altitudes and weather conditions, as well as bolides and other natural phenomena that may be detected by the international infrasound monitoring network. Once the data are all in digital form, we will convert them to the standard CSS format, including event and station information. The complete set of database tables and binary waveform files will be the ultimate product of our work.
In order to exploit the information on surface wave propagation that is stored in large seismic event datasets, Sandia and Lawrence Livermore National Laboratories have developed a MatSeis interface for performing phase-matched filtering of Rayleigh arrivals. MatSeis is a Matlab-based seismic processing toolkit which provides graphical tools for analyzing seismic data from a network of stations. Tools are available for spectral and polarization measurements, as well as beam forming and f-k analysis with array data, to name just a few. Additionally, one has full access to the Matlab environment and any functions available there. Previously the authors reported the development of new MatSeis tools for calculating regional discrimination measurements. The first of these performs Lg coda analysis as developed by Mayeda and coworkers at Lawrence Livermore National Laboratory. A second tool measures regional phase amplitude ratios for an event and compares the results to ratios from known earthquakes and explosions. Release 1.5 of MatSeis includes the new interface for the analysis of surface wave arrivals. This effort involves the use of regionalized dispersion models from a repository of surface wave data and the construction of phase-matched filters to improve surface wave identification, detection, and magnitude calculation. The tool works as follows. First, a ray is traced from source to receiver through a user-defined grid containing different group velocity versus period values to determine the composite group velocity curve for the path. This curve is shown along with the upper and lower group velocity bounds for reference. Next, the curve is used to create a phase-matched filter, apply the filter, and show the resultant waveform. The application of the filter allows obscured Rayleigh arrivals to be more easily identified. Finally, after screening information outside the range of the phase-matched filter, an inverse version of the filter is applied to obtain a cleaned raw waveform which can be used for amplitude measurements. Because all the MatSeis tools have been written as Matlab functions, they can be easily modified to experiment with different processing details. The performance of the propagation models can be evaluated using any event available in the repository of surface wave events.
The United States conducted over 100 atmospheric nuclear tests at the Nevada Test Site from 1951 through 1962. Some of the earliest tests caused unexpected damage, primarily broken glass and cracked plaster, in Las Vegas and other surrounding communities. To address this problem, Sandia initiated a program to monitor and predict the pressure waves around NTS. Infrasound recording systems were developed, then field for all tests beginning with Operation Buster in October 1951. Investigators soon discovered that near-surface temperature inversions and wind profiles caused the damaging pressures in Las Vegas. A typical test was recorded at about a dozen stations from the Control Point on NTS to as far away as Pasadena, CA. In addition, some tests in the South Pacific were monitored, as well as numerous chemical explosions. Strip charts recorded signals in the frequency band from 0.05 to 30 Hz, and the paper tapes were achieved at Sandia in the early 1970s. The NTS events ranged in yield from below 1 ton to 74 kilotons; source altitudes varied from near ground level (including some cratering experiments) to as high as 11 km. The resulting data contain a wealth of information on the source function, yield scaling and regional propagation of infrasound signals from atmospheric explosions. The renewed interest in infrasonic monitoring for CTBT verification has prompted the authors to exhume some of the archived records. The authors plan to digitize the signals from several tests and evaluate their applicability to CTBT issues. In addition, they will collect any existing parametric measurements for these records (arrival times, amplitudes, etc.). All data will be converted to CSS database format and made available to the research community. If appropriate, the resulting information could also be included in the Knowledge Base under development for CTBT monitoring.
The Seismic Verification Program at Sandia focuses on designing and building seismic monitoring systems which could be deployed within the Soviet Union. To support this effort, Sandia also is involved in developing and testing seismic components, and in evaluating system performance and new monitoring techniques. Seismic studies conducted under the latter task concentrate on analyzing regional seismic signals, recorded within 2000 km of the source. In-country monitoring stations would be able to exploit these regional signals to improve on the capabilities of external stations at teleseismic ranges. The principal advantages which regional signals offer are greater amplitudes and broader frequency bands. However, such signals are more complex and less well understood than teleseismic records. Ongoing studies at Sandia will help define the performance which can be expected from an in-country network. Other studies are testing new monitoring concepts which may further improve in-country capabilities. This presentation gives a brief overview of some of these projects.