A MATHEMATICAL STATISTICS FORMULATION OF THE TELESEISMIC EXPLOSION IDENTIFICATION PROBLEM WITH MULTIPLE DISCRIMINANTS
Seismic monitoring for underground nuclear explosions answers three questions for all global seismic activity: Where is the seismic event located? What is the event source type (event identification)? If the event is an explosion, what is the yield? The answers to these questions involve processing seismometer waveforms with propagation paths predominately in the mantle. Four discriminants commonly used to identify teleseismic events are depth from travel time, presence of long-period surface energy (mb versus Ms ), depth from reflective phases, and polarity of first motion. The seismic theory for these discriminants is well established in the literature (see for example Pomeroy et al. (1982) and Blandford (1982)). However the physical basis of each has not been formally integrated into probability models to account for statistical error and provide discriminant calculations generally appropriate for multi-dimensional event identification. This paper develops a mathematical statistics formulation of these discriminants and offers a novel approach to multi-dimensional discrimination that is readily extensible to other discriminants. For each discriminant a probability model is formulated under a general null hypothesis of H0 : Explosion Characteristics. The veracity of the hypothesized model is measured with a p-value calculation (see Stuart et al. (1994) and Freedman et al. (1991)) that is filtered to be approximately normally distributed and is in the range [0, 1]. A value near zero rejects H0 and a moderate to large value indicates consistency with H0 . The hypothesis test formulation ensures that seismic phenomenology is tied to the interpretation of the p-value. These p-values are then embedded into a multi-discriminant algorithm that is developed from regularized discrimination methods proposed by Smidt and McDonald (1976), DiPillo (1976) and Friedman (1989). Performance of the methods is demonstrated with 102 teleseismic events with magnitudes (mb ) ranging from 5 to 6.5. Example p-value calculations are given for two of these events. Statistical properties of p-values are also presented.
- Research Organization:
- Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- AC05-76RL01830
- OSTI ID:
- 1327203
- Report Number(s):
- PNNL-SA-56382; NN2003000
- Resource Relation:
- Journal Volume: 97; Journal Issue: 5; Conference: Proceedings of the 29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies, September 25-27, 2007, Denver, Colorado, 526-530
- Country of Publication:
- United States
- Language:
- English
Similar Records
Regional discrimination of NTS (Nevada Test Site) explosions and western US earthquakes: Multivariate discriminants
Seismic data analysis center. Final techical report