Seismic monitoring of explosions: a method to extract information on the isotropic component of the seismic source

The design of a monitoring system for detecting explosions is a very topical problem, both for routine data processing at seismological observatories as well as for the monitoring of a Comprehensive Test Ban Treaty. In this framework it is desirable to have the possibility to quantify the presence of the isotropic component in the seismic source. For this purpose a method is presented, which is based on waveform inversion for the full moment tensor retrieval. The method inverts either full waveforms or separate seismic phases and returns the mechanism and time history of a point source. Moreover, it allows to redefine the hypocentral depth of the event and, in a simplistic way, to optimize the structural model as well. In order to model strong laterally heterogeneous structures, different pairs of structural models can be used for each source-receiver path. The source is decomposed into a volumetric part (V), representing an explosive or implosive component, and into a deviatoric part, containing both the double couple (DC) and the compensated linear vector dipole (CLVD) components. The method is applied to an area in central Switzerland and to the network of the Swiss Seismological Service. The events of interest include both earthquakes and explosions. Despite some modelling inadequacies of the source-time function, the explosions can be well identified with the inverted isotropic component in the source, as long as the number of stations used for the inversion is larger than three. The results of the inversion are better for large epicenter-station distances of the order of 40–90 km.

[1]  Charles A. Langston,et al.  Source inversion of seismic waveforms: The Koyna, India, earthquakes of 13 September 1967 , 1981 .

[2]  John H. Woodhouse,et al.  Determination of earthquake source parameters from waveform data for studies of global and regional seismicity , 1981 .

[3]  Nicholas Deichmann,et al.  Structural and rheological implications of lower-crustal earthquakes below northern Switzerland , 1992 .

[4]  G. F. Panza,et al.  CHAPTER 4 – Complete Strong Motion Synthetics , 1987 .

[5]  Donat Fäh,et al.  Complete synthetic seismograms for high-frequency multimode SH -waves , 1991 .

[6]  NábÄ lek,et al.  Determination of earthquake source parameters from inversion of body waves , 1984 .

[7]  E. Kissling,et al.  Crustal structure beneath the eastern Swiss Alps derived from seismic refraction data , 1995 .

[8]  Giuliano F. Panza,et al.  Seismic moment tensor resolution by waveform inversion of a few local noisy records—I. Synthetic tests , 1996 .

[9]  K. Koch,et al.  Moment tensor inversion of local earthquake data—I. Investigation of the method and its numerical stability with model calculations , 1991 .

[10]  H. Kanamori,et al.  A moment magnitude scale , 1979 .

[11]  W. H. Lee,et al.  HYPO71 (revised; a computer program for determining hypocenter, magnitude, and first motion pattern of local earthquakes , 1972 .

[12]  Giuliano F. Panza,et al.  Waveform inversion for point source moment tensor retrieval with variable hypocentral depth and structural model , 1992 .

[13]  H. Maurer,et al.  Crustal structure beneath the northern margin of the Swiss Alps , 1992 .

[14]  F. Gilbert,et al.  An application of normal mode theory to the retrieval of structural parameters and source mechanisms from seismic spectra , 1975, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[15]  René Egloff,et al.  Sprengseismische Untersuchungen der Erdkruste in der Schweiz , 1979 .

[16]  G. Cespuglio,et al.  Seismic moment tensor resolution by waveform inversion of a few local noisy records—II. Application to the Phlegraean Fields (Southern Italy) volcanic tremors , 1996 .

[17]  G. F. Panza,et al.  Complete moment tensor retrieval for weak events: application to orogenic and volcanic areas , 1996 .

[18]  S. Sipkin,et al.  Estimation of earthquake source parameters by the inversion of waveform data: synthetic waveforms , 1982 .

[19]  P. Valasek,et al.  The tectonic structure of the Swiss Alpine crust interpreted from a 2D network of deep crustal seismic profiles and an evaluation of 3D effects , 1992 .