On the linear relation between mb and Ms for discrimination between explosions and earthquakes

Summary. The statistical capability of the m _b: M _s discriminant for the discrimination of earthquake and explosion populations is examined by application of discriminant functions to a group of 83 explosions and 72 earthquakes in Eurasia. Equations are derived for the probability that an event is an earthquake or an explosion. The positive sign of DIS in the decision index equation, DIS _i = 34.3383 – 11.9569 mb _t + 7.1161 M _si , indicates that the event i is an earthquake. Its negative sign indicates that event i is an explosion. The probability of correct classification for an event, P _i , is related to its DIS _i value, by P _i = [1-exp (DIS _i )]⁻¹, where a large, positive DIS indicates a high probability that an event is an earthquake and a large, negative DIS indicates a high probability that an event is an explosion. The discrimination line M _s = 1.680 m _b– 4.825, or m _b= 0.595 M _s+ 2.872 very successfully separates the explosion population from the earthquake population. The points on this line have an equal chance of being an earthquake or an explosion; moreover, for any event, the distance parallel to the M _s-axis from the point representing that event in the m _b: M _s plane to this line is a measure of the probability for the correct classification of that event.     

Dynamics of explosive volcanism at Unzen volcano: an experimental contribution

Knowledge of the dynamics of magma fragmentation is necessary for a better understanding of the explosive behaviour of silicic volcanoes. Here we have measured the fragmentation speed and the fragmentation threshold of five dacitic samples (6.7–53.5 vol% open porosity) from Unzen volcano, Kyushu, Japan. The measurements were carried out using a shock-tube-based fragmentation apparatus modified after Alidibirov and Dingwell (1996a,b). The results of the experimental work confirm the dominant influence of porosity on fragmentation dynamics. The velocity of the fragmentation front increases and the value of the fragmentation threshold decreases with increasing porosity. Further, we observe that the fragmentation speed is strongly influenced by the initial pressure difference and the texture of the dacite. At an initial pressure difference of 30 MPa, the fragmentation speed varies from 34 m/s for the least porous sample to 100 m/s for the most porous sample. These results are evaluated by applying them to the 1990–1995 eruptive activity of Unzen volcano. Emplacements of layered lava dome lobes, Merapi-type pyroclastic flows and minor explosive events dominated this eruption. The influence of the fragmentation dynamics on dome collapse and Vulcanian events is discussed.     

Explosive activity and generation mechanisms of pyroclastic flows at Arenal volcano, Costa Rica between 1987 and 2001

Explosive activity at Arenal and associated tephra fall that has occurred over the 14-year period from 1987–2001 is described. Explosions have been notably variable in both frequency and size. A marked decrease in both frequency and quantity of tephra fallout occurred in early 1998 until the end of 2001. Grainsize distributions of cumulative tephra samples collected once a month are typically bimodal. Aggregation causing premature fallout of fine ash and possibly fallout from ash plumes produced by pyroclastic flows are considered responsible for the bimodality of fallout. Scanning electron microscopy of the glass component of tephra from single explosions show predominantly blocky and blocky/fluidal clast types, interpreted as being the product of vulcanian type explosions. Fragmentation of a mainly rigid, degassed magma body, and a minor molten component is inferred for these explosions. Pyroclastic flows were produced either associated with the larger explosions by a mechanism of column collapse (1987–1990), or unrelated to explosions by partial collapse of the crater wall (1993, 1998, 2000, 2001). Pyroclastic flow activity has migrated from west to north during the period reported. Pyroclastic flow deposits are variable in the quantity of juvenile material and any associated surge component. Large juvenile blocks were partially molten on emplacement and many have a typical cauliform texture. Blocks with both juvenile and lithic textures indicate that at the summit magma was in intimate contact with the pre-existing edifice, rather than as a simple open crater or lava pool. Crater wall collapse may have been promoted by the reduction in explosive activity, which has increased the lava accumulation at the summit and in turn increased instability of the summit region. Thus although explosive activity has waned, if the lava output is maintained, the hazard of pyroclastic flows is likely to continue.Editorial responsibility: R. Cioni     

2D model of the crust and uppermost mantle along rift profile, Siberian craton

A two-dimensional model of the crust and uppermost mantle for the western Siberian craton and the adjoining areas of the Pur-Gedan basin to the north and Baikal Rift zone to the south is determined from travel time data from recordings of 30 chemical explosions and three nuclear explosions along the RIFT deep seismic sounding profile. This velocity model shows strong lateral variations in the crust and sub-Moho structure both within the craton and between the craton and the surrounding region. The Pur-Gedan basin has a 15-km thick, low-velocity sediment layer overlying a 25-km thick, high-velocity crystalline crustal layer. A paleo-rift zone with a graben-like structure in the basement and a high-velocity crustal intrusion or mantle upward exists beneath the southern part of the Pur-Gedan basin. The sedimentary layer is thin or non-existent and there is a velocity reversal in the upper crust beneath the Yenisey Zone. The Siberian craton has nearly uniform crustal thickness of 40–43 km but the average velocity in the lower crust in the north is higher (6.8–6.9 km/s) than in the south (6.6 km/s). The crust beneath the Baikal Rift zone is 35 km thick and has an average crustal velocity similar to that observed beneath the southern part of craton. The uppermost mantle velocity varies from 8.0 to 8.1 km/s beneath the young West Siberian platform and Baikal Rift zone to 8.1–8.5 km/s beneath the Siberian craton. Anomalous high Pn velocities (8.4–8.5 km/s) are observed beneath the western Tunguss basin in the northern part of the craton and beneath the southern part of the Siberian craton, but lower Pn velocities (8.1 km/s) are observed beneath the Low Angara basin in the central part of the craton. At about 100 km depth beneath the craton, there is a velocity inversion with a strong reflecting interface at its base. Some reflectors are also distinguished within the upper mantle at depth between 230 and 350 km.     

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.     

Observations of teleseismicP wave coda for underground explosions

The earlyP wave coda (5–15 sec after the first arrival) of underground explosions at the Nevada Test Site is studied in the time domain using 2082 teleseismic short-period recordings, with the intent of identifying near-source contributions to the signals in the frequency range 0.2–2.0 Hz. Smaller magnitude events tend to have relatively high coda levels in the 0.4–0.8 Hz frequency band for both Yucca Flat and Pahute Mesa explosions. Coda complexity in this low-frequency passband is negatively correlated with burial depth for Pahute Mesa events but is only weakly correlated with depth for Yucca Flat events. Enhanced excitation of relatively long-period scattered waves for smaller, less deeply buried events is required to explain this behavior. Coda complexity in the 0.8–1.1 Hz band is positively correlated with magnitude and depth for Pahute Mesa events, but has no such dependence for Yucca Flat events. This may result from systematic variations between the spectra of direct signals and coda arrivals caused bypP interference for the largest events, all of which were detonated at Pahute Mesa. Another possible explanation is a frequency-dependent propagation effect on the direct signals of the larger events, most of which were located in the center of the mesa overlying strong lateral velocity gradients in the crust and upper mantle. Event average complexity varies spatially for both test sites, particularly in the 0.8–1.1 Hz band, providing evidence for frequency-dependent focussing or scattering by near-source structure. Both the direct arrivals and the early coda have strong azimuthal amplitude patterns that are produced by defocussing by mantle heterogeneity. The direct arrivals have stronger coherent azimuthal patterns than the early coda for Pahute Mesa events, indicating more pronounced deep crustal and shallow mantle defocussing for the direct signals. However, for Yucca Flat events the direct arrivals have less coherent azimuthal patterns than the coda, suggesting that a highly variable component of near-source scattering preferentially affecting the downgoing energy is superimposed on a pattern produced by mantle heterogeneity that affects the entire signal. This complicated behavior of the direct arrivals may be the result of triplications and caustics produced by the complex basement structure known to underlie the Yucca Flat test site. The presence of strong azimuthal patterns in the early coda indicates that source strength estimates based on early coda are subject to biases similar to those affecting estimates based on direct arrivals.     

P/S震相幅值比判据对低震级地震事件的适用性检验c

基于中国新疆及中亚地区大量天然地震和地下爆炸，检验了较大震级情况下得到的区域P/S震相幅值比识别判据对低震级地震事件的适用性. 分析结果表明，对较大震级地震事件具有较好识别效果的P/S震相幅值比判据对低震级事件同样适用，而且对低震级事件也具有与较大震级事件大致相当的识别效果. 据WMQ、BLK、MUL和MAK台的25个识别效果较好的P/S震相幅值比判据统计， 大于ML3.5和ML3.5以下地震事件的误识率之差约为2个百分点.