Global Detection of Infrasonic Signals from Three Large Bolides

Transcription

1 Earth Moon Planet (2008) 102: DOI /s z Global Detection of Infrasonic Signals from Three Large Bolides Stephen J. Arrowsmith Æ Doug ReVelle Æ Wayne Edwards Æ Peter Brown Received: 28 June 2007 / Accepted: 29 November 2007 / Published online: 13 December 2007 Ó Springer Science+Business Media B.V Abstract We present the infrasonic observations of three large bolides that were observed at numerous International Monitoring System (IMS) infrasound arrays on a global scale. First, a simple procedure for the global association of infrasound detections from large infrasound events is outlined. Infrasound signals are associated with large events based on arrival time, backazimuth and uniqueness at a given IMS array. Next, we apply the algorithm to three bolides and investigate some of the factors affecting the detectability of infrasound from large events. Our findings suggest that site-noise effects significantly degrade the capability of the IMS infrasound network, suggesting that more effort is required to reduce ambient site noise. These results have implications for the use of infrasound measurements (in particular those from IMS stations) as a tool for evaluating the global flux of near-earth objects. Keywords Bolide infrasound Meteor detection International Monitoring System 1 Introduction Infrasound is acoustic energy that propagates at frequencies below the 20 Hz hearing threshold of the human ear. Unlike audible sound, infrasound can travel for long distances with relatively little attenuation. Large meteors generate infrasound as they enter Earths atmosphere (ReVelle 1976), which can propagate across the globe. The purpose of this paper is to outline a simple method for identifying infrasonic signals from large events on a global scale, and to apply the method to three superbolides. Following Docobo et al. (1998), a superbolide has a peak optical luminosity that exceeds -17 stellar magnitudes (referenced to an elevation of 100 km in the zenith). We also provide an assessment of the S. J. Arrowsmith (&) D. ReVelle Atmospheric, Climate and Environmental Dynamics Group, Los Alamos National Laboratory, Los Alamos, NM 87545, USA sarrowsmith@gmail.com W. Edwards P. Brown Department of Physics and Astronomy, University of Western Ontario, London, ON, Canada

2 358 S. J. Arrowsmith et al. effect of site noise on the detectability of these large events. There have been relatively few studies of multi-station infrasound signals from large events on a global scale. Long before the advent of the present-day IMS, Wexler and Hass (1962) documented the global detection of infrasound signals from a large Soviet nuclear test. More recently, Brown et al. (2002) present the observations of infrasound at numerous IMS stations from two large bolides. Finally, Garces et al. (2005) and Le Pichon et al. (2006) document multistation infrasonic signals from the 2004 Sumatra earthquake and tsunami and the 2005 Chilean earthquake respectively. In contrast with earlier studies, this study focuses on the more general problem of associating infrasound signals with global events, and on an assessment of factors that affect the global detectability of large events. We report on the global detection of three superbolides which entered Earths atmosphere on September 3rd, 2004 over Antarctica; on October 7th, 2004 over the Indian Ocean; and on December 9th, 2006 over North Africa. 2 Global Infrasound Event Association The problem of associating an infrasound signal with a known event is non-trivial. The challenge can be phrased as follows: How can we confidently associate a large event (with Fig. 1 Flowchart illustrating the event-based global association procedure (schematic diagram for illustration purposes only)

3 Global Detection of Infrasonic Signals from Three Large Bolides 359 known event location and origin time) with an infrasound signal at a range of 100s to 10,000s of kilometers? Typically, event-based association utilizes the consistency of arrival-time and backazimuth at a given array with predicted values based on the known event location and origin time. However, standard detection algorithms can generate large numbers of coherent infrasonic signals with low signal-to-noise ratios (e.g., Brachet and Coyne 2006). This results in a significant risk of incorrectly associating an infrasonic signal with a given event. To mitigate against the possibility of mis-association, we introduce an additional requirement that a signal associated with a large bolide must be unique for a given array in a time window before and after the signal. Subsequently, we have developed a three-stage procedure for associating infrasound signals with large bolides (Fig. 1). First, for a known event location and origin time, we compute predicted arrival times and backazimuths at all IMS arrays located around the globe. Second, we apply an array-based signal detection algorithm to the data the PMCC algorithm (Cansi 1995). The basis of the PMCC method, in common with other array-based detectors, is that signals are spatially coherent arriving as plane waves, whereas noise is incoherent between the different elements in an array. By time-aligning arrivals at the different elements in an array, PMCC provides backazimuths and trace velocities for the incoming wave useful information for event association. We search for PMCC detections in a time-window encompassing the range of all possible infrasonic arrivals (where the earliest possible signal propagates with a horizontal group velocity of 0.34 km/s, and the latest signal propagates with a group velocity of 0.22 km/s, Ceplecha et al. 1998). We then consider detections within an allowed backazimuth deviation of the true great-circle path, occuring within the predicted time window, to be potential associations. The allowed backazimuth deviations are based on an empirical relation calculated by the U.S. Air Force for known explosions. Finally, we remove potential associations that are: (a) short-duration and narrow-band (and Fig. 2 Map showing stations that detected the September 3rd, 2004 Antarctic bolide (black triangles) and stations that did not (white circles). The event location is at the center of the map

4 360 S. J. Arrowsmith et al. Fig. 3 Example signals from the September 3rd, 2004 Antarctic superbolide at (a) the I55US array at Windless Bight, Antarctica; and (b) the I17CI array in the Ivory Coast. In (b), the signal from the superbolide is outlined by the dashed lines

5 Global Detection of Infrasonic Signals from Three Large Bolides 361 Table 1 Summary of associated detections at IMS infrasound arrays for the three superbolide events Event Array Range (km) Arrival time (UT) Duration (s) Peak amplitude (Pa) Back-azimuth ( ) 09/03/04 I27DE :15: /03/04 I55US :50: /03/04 I35NA :20: /03/04 I05AU :47: /03/04 I17CI :26: /03/04 I26DE 12,918 00:27: (on 09/04) 10/07/04 I52GB :50: /07/04 I32KE :38: /07/04 I55US :01: /07/04 I17CI :45: /07/04 I26DE :02: /07/04 I10CA :17: (on 10/08) 12/09/06 I26DE :08: /09/06 I35NA :40: /09/06 I30JP :13: /09/06 I41PY :05: /09/06 I56US :11: /09/06 I04AU :57: Note: The signal durations observed for these three superbolides do not increase as a function of range, as would be predicted theoretically. We speculate that the reason for this discrepancy is due to site noise therefore unlikely to be associated with a large, distant source), and (b) similar to other detections obtained before and/or after the signal time window (Fig. 1). Test (b) reduces the likelihood of incorrectly associating an unrelated signal with the event of interest, since detections from coherent noise at a given station (which are typically highly repetitive) will be removed. The criteria used for measuring the similarity with other detections include backazimuth, frequency and trace velocity. The threshold constraints for matching two detections are set as follows: Maximum deviation in backazimuth = 5, Maximum deviation in frequency = 1.2 Hz, Maximum deviation in trace velocity = 0.1 km/s. These constraints increase the likelihood that the signal is unique to a given event, i.e. it does not occur before or after the possible arrival window for an event of interest (Fig. 1). 3 Results and Conclusions The event-based association algorithm described above is applied to three superbolides, which occurred from 2004 to Figure 2 shows the locations of arrays that detected one of the three events, a superbolide that entered over Antarctica on September 3rd, 2004 (for further details on this event, see Klekociuk et al. 2005). The event was detected by six IMS arrays. As shown in Fig. 2, the distribution of observations is strongly asymmetrical, suggesting that detectability is not a simple function of range. An example signal, recorded at the I55US infrasound array (which is located at Windless Bight in

6 362 S. J. Arrowsmith et al. Fig. 4 Noise spectra for each of the three bolides. Solid black lines denote spectra for stations that detected the events and gray dashed lines denote spectra for stations that did not (for stations at ranges less than *10,000 km). Regions shaded gray denote the frequency band over which detections were observed for each bolide Antarctica), is shown in Fig. 3a. A long duration (*10 min long) signal is observed, with a backazimuth of 209 ± 4.8. The relatively fixed backazimuth and trace velocity indicate the signal is from a single event. The clear signal correlation between the separate array elements can be observed. Signals at greater ranges are typically lower in amplitude and harder to pick out from the ambient noise, requiring array-processing techniques to extract their full extent. An example of such a signal from the September 3rd, 2004 Antarctic bolide is shown in Fig. 3b, which was recorded at the I17CI infrasound array in the Ivory Coast (see map in Fig. 2). Similar patterns of observations were observed for the two other events, with six detections, and an asymmetric distribution of observations observed in each case. A summary of the observations from all three bolides is provided in Table 1. Whether or not a given station detects an event appears to be strongly governed by factors other than range. In particular, two factors that influence the detectability of infrasound from an event are: (1) ambient noise levels at a given station, and (2) propagation effects. In order to assess the importance of ambient noise levels on the signal detection, we have computed pre-event noise spectra for each station. Each power spectrum was computed for a beamformed waveform at each array using a time window of * 1 h (i.e., 2 16 data points), and smoothed by averaging the power spectrum over

7 Global Detection of Infrasonic Signals from Three Large Bolides 363 1/8 octave intervals. The resultant power spectra for stations at ranges \10,000 km are plotted in Fig. 4. The plot demonstrates that for two events there is a clear separation between noise-levels at stations that detected the event and stations that did not detect the events. In fact, for these two events the difference in noise levels between stations that detected the events, and stations that did not, is approximately 20 db on average. This suggests that site-noise effects play a dominant role on the detectability of infrasound from such large events. We note that such a separation is not observed for more distant stations ([10,000 km), suggesting that propagation effects may dominate noise effects at larger ranges. We have not performed detailed propagation modeling in this study but plan to do so for a future publication. As discussed previously, both site noise and propagation effects will affect the detectability of infrasound. However, our present results clearly suggest that site noise effects are a significant factor on the global detectability of the three superbolides. Acknowledgements We gratefully appreciate financial support from the GNEM Program of the US DOE-HQ in NA-22. Waveform data were obtained from the SMDC monitoring website. References N. Brachet, J. Coyne, The Current Status of Infrasound Data Processing at the International Data Centre, Proceedings of the 28th Seismic Research Review, Orlando, Florida (2006) P.G. Brown, R.W. Whitaker, D.O. ReVelle, Multi-station infrasonic observations of two large bolides: signal interpretation and implications for monitoring of atmospheric explosions. Geophys. Res. Lett. 29 (2002). doi: /2001gl Y. Cansi, An automated seismic event processing for detection and location: The P.M.C.C. method. Geophys. Res. Lett. 22, (1995) Z. Ceplecha, J. Borovicka, W.G. Elford, D.O. ReVelle, R.L. Hawkes, V. Porubcan, M. Simel, Meteor phenomena and bodies. Space Sci. Rev. 84, (1998) J.A. Docobo, R.E. Spalding, Z. Ceplecha, F. Diaz-Fierros, V. Tamazian, Y. Onda, Investigation of a bright flying object over Northwest Spain, 1994 January 18. Meteorit Planet. Sci. 33, (1998) M. Garces, P. Caron, C. Hetzer, A. Le Pichon, H. Bass, D. Drob, J. Bhattacharyya, Deep infrasound radiated by the Sumatra earthquake and Tsunami. Eos 86, (2005) A.R. Klekociuk, P.G. Brown, D.W. Pack, D.O. ReVelle, W.N. Edwards, R.E. Spalding, E. Tagliaferri, B.B. Yoo, J. Zagari, Meteoritic dust from the atmospheric disintegration of a large meteoroid. Nature 436, (2005) A. Le Pichon, P. Mialle, J. Guilbert, J. Vergoz, Multistation infrasonic observations of the Chilean earthquake of 2005 June 13. Geophys. J. Int. 167, (2006) D. ReVelle, On meteor-generated infrasound, J. Geophys. Res. 81, (1796) H. Wexler, W.A. Hass, Global atmospheric pressure effects of the October 30, Explosion. J. Geophys. Res. 67, (1962)