Probing mid-mantle heterogeneity using PKP coda waves

We investigate the utility of PKP coda waves for studying weak scattering from small-scale heterogeneity in the mid-mantle. Coda waves are potentially a useful probe of heterogeneity in the mid-mantle because they are not preferentially scattered near the CMB, as PKP precursors are, but are sensitive to scattering at all depths. PKP coda waves have not been used for this purpose historically because of interference with other late-arriving energy due to near-surface resonance and scattering. Any study of deep mantle scattering using coda waves requires the removal of near-surface effects from the data. We have analyzed 3624 recordings of PKP precursors and coda made by stations in the Incorporated Research Institutions for Seismology (IRIS) Global Seismographic Network (GSN). To study the range and time dependence of the scattered waves, we binned and stacked envelopes of the recordings. We have considered precursors that arrive within a 20 s window before PKP and coda waves in a 60 s window after PKP. The PKP scattered waves increase in amplitude rapidly with range as predicted by scattering theory. At ranges below 125°, we predict and observe essentially no scattered energy preceding PKP. Coda amplitudes at these ranges are independent of range and provide an estimate of energy due to near-surface effects that we can expect at all ranges. We use the average coda amplitude at ranges from 120 to 125° to correct coda amplitudes at other ranges. PKP coda waves show a strong dependence on time and range and are clearly influenced by scattering in the lower mantle. PKP coda waves, however, do not provide a tighter constraint on the vertical distribution of mantle heterogeneity than is provided by precursors. This is due, in part, to relatively large scatter in coda amplitudes as revealed by a resampling analysis. Modeling using Rayleigh–Born scattering theory and an exponential autocorrelation function shows that PKP coda amplitudes are not highly sensitive to the vertical distribution of heterogeneity in the mantle. To illustrate this we consider single-scattering in two extreme models of mantle heterogeneity. One allows heterogeneity just at the CMB; the other includes heterogeneity throughout the mantle. The amplitudes of precursors are tightly constrained by our stack and support our earlier conclusion that small-scale heterogeneity is uniformly distributed throughout the lower mantle. The best-fit model includes 8 km scale length heterogeneity with an rms velocity contrast throughout the mantle of 1%.     

Identification of Mining Blasts at Mid- to Far-regional Distances Using Low Frequency Seismic Signals

—?This paper reports results from two recent monitoring experiments in Wyoming. Broadband seismic recordings of kiloton class delay-fired cast blasts and instantaneous calibration shots in the Black Thunder coal mine were made at four azimuths at ranges from 1° to 2°. The primary focus of this experiment was to observe and to explain low-frequency signals that can be seen at all azimuths and should routinely propagate above noise to mid-regional distances where most events will be recorded by International Monitoring System (IMS) stations.¶The recordings clearly demonstrate that large millisecond delay-fired cast blasts routinely produce seismic signals that have significant spectral modulations below 10?Hz. These modulations are independent of time, the azimuth from the source and the orientation of the sensor. Low-frequency modulations below 5?Hz are seen beyond 9°. The modulations are not due to resonance as they are not produced by the calibration shots. Linear elastic modeling of the blasts that is guided by mine-blast reports fails to reproduce the fine detail of these modulations but clearly indicates that the enhanced “spectral roughness” is due to long interrow delays and source finiteness. The mismatch between the data and the synthetics is likely due to source processes, such as nonlinear interactions between shots, that are poorly understood and to other effects, such as variations of shot time and yield from planned values, that are known to be omnipresent but cannot be described accurately. A variant of the Automated Time-Frequency Discriminant (Hedlin, 1998b), which uses low-frequency spectral modulations, effectively separates these events from the calibration shots.¶The experiment also provided evidence that kiloton class cast blasts consistently yield energetic 2–10 second surface waves. The surface waves are strongly dependent on azimuth but are seen beyond 9°. Physical modeling of these events indicates that the surface waves are due mainly to the extended source duration and to a lesser extent to the slap-down of spalled material. The directionality is largely a path effect. A discriminant that is based on the partitioning of energy between surface and body waves routinely separates these events from the calibration shots.¶The Powder River Basin has essentially no natural seismic activity. How these mining events compare to earthquake observations remains to be determined.     

The Temporal Morphology of Infrasound Propagation

Expert knowledge suggests that the performance of automated infrasound event association and source location algorithms could be greatly improved by the ability to continually update station travel-time curves to properly account for the hourly, daily, and seasonal changes of the atmospheric state. With the goal of reducing false alarm rates and improving network detection capability we endeavor to develop, validate, and integrate this capability into infrasound processing operations at the International Data Centre of the Comprehensive Nuclear Test-Ban Treaty Organization. Numerous studies have demonstrated that incorporation of hybrid ground-to-space (G2S) enviromental specifications in numerical calculations of infrasound signal travel time and azimuth deviation yields significantly improved results over that of climatological atmospheric specifications, specifically for tropospheric and stratospheric modes. A robust infrastructure currently exists to generate hybrid G2S vector spherical harmonic coefficients, based on existing operational and emperical models on a real-time basis (every 3- to 6-hours) (Drob et al., 2003). Thus the next requirement in this endeavor is to refine numerical procedures to calculate infrasound propagation characteristics for robust automatic infrasound arrival identification and network detection, location, and characterization algorithms. We present results from a new code that integrates the local (range-independent) τp ray equations to provide travel time, range, turning point, and azimuth deviation for any location on the globe given a G2S vector spherical harmonic coefficient set. The code employs an accurate numerical technique capable of handling square-root singularities. We investigate the seasonal variability of propagation characteristics over a five-year time series for two different stations within the International Monitoring System with the aim of understanding the capabilities of current working knowledge of the atmosphere and infrasound propagation models. The statistical behaviors or occurrence frequency of various propagation configurations are discussed. Representative examples of some of these propagation configuration states are also shown.     

Surveying Infrasonic Noise on Oceanic Islands

v--vAn essential step in the establishment of an International Monitoring System (IMS) infrasound station is the site survey. The survey seeks a location with relatively low infrasonic noise and the necessary logistical support. This paper reports results from our surveys of two of the oceanic sites in the IMS - the Azores and Cape Verde. Each survey sampled infrasonic noise, wind velocity, air temperature and humidity for ~3 weeks at 4 sites near the nominal IMS locations. The surveys were conducted on Sao Miguel (the main island in the Azores) and Maio (Cape Verde). Infrasonic noise was measured using the French MB2000 microbarometer.¶During our 3-week experiment in January the trade winds at Cape Verde varied little from an azimuth of 63°. Because of the unvarying wind azimuth, the experiment gave us an opportunity to examine the effectiveness of a forest at reducing both wind speed and infrasonic noise. We find that the thick Acacia forest on Maio reduces wind speeds at a 2 m elevation by more than 50% but does not reduce infrasonic noise at frequencies below 0.25 Hz. This forest serves as a high-frequency filter and clearly does not reduce long-period noise levels which are due to large-scale turbulence in the atmospheric boundary layer above the forest. This is consistent with our observations in the Azores where the relationship between infrasonic noise and wind speed is more complex due to frequent changes in wind azimuth.¶In Cape Verde, wind speed and infrasonic noise are relatively constant. The diurnal variations are clearly seen however the microbarom is only rarely sensed. In the Azores, during our 3-week experiment in November and December of 1998, wind speed and infrasonic noise change rapidly. At this location, daily noise level swings of 40 to 50 dB at 0.1 Hz are not uncommon in the early winter and are due to changes in wind speed and atmospheric turbulence. The effectiveness of an infrasound station in the Azores will be strongly dependent on time during the winter season.¶The two surveys illustrate some of the difficulties inherent in the selection of sites for 1 to 3 km aperture arrays on oceanic islands. Due to elevated noise levels at these sites, 8 element, 2 km aperture arrays are strongly preferred.