[1] Seismic data recorded in the 2-30 s band at Stromboli Volcano, Italy, are analyzed to quantify the source mechanisms of Strombolian explosions during September 1997. To determine the source-centroid location and source mechanism, we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous elastic medium that takes topography into account. Two source centroids are identified, each representative of the distinct event types associated with explosive eruptions from two different vents. The observed waveforms are well reproduced by our inversion, and the two source centroids that best fit the data are offset 220 and 260 m beneath and $160 m northwest of the active vents. The source mechanisms include both moment-tensor and single-force components. The principal axes of the moment tensor have amplitude ratios 1:1:2, which can be interpreted as representative of a crack, if one assumes the rock matrix at the source to have a Poisson ratio n = 1/3, a value appropriate for hot rock. Both imaged cracks dip $60°to the northwest and strike northeast-southwest along a direction parallel to the elongation of the volcanic edifice and a prominent zone of structural weakness, as expressed by lineaments, dikes, and brittle structures. For our data set, the volume changes estimated from the moments are $200 m 3 for the largest explosion from each vent. Together with the volumetric source is a dominantly vertical force with a magnitude of 10 8 N, consistent with the inferred movement of the magma column perched above the source centroid in response to the piston-like rise of a slug of gas in the conduit.

Abstract. We use data from broadband seismometers deployed around the summit of Kilauea Volcano to quantify the mechanism associated with a transient in the flow of magma feeding the east rift eruption of the volcano. The transient is marked by rapid inflation of the Kilauea summit peaking at 22 ttrad 4.5 hours after the event onset, followed by slow deflation over a period of 3 days. Superimposed on the summit inflation is a series of sawtooth displacement pulses, each characterized by a sudden drop in amplitude lasting 5-10 s followed by an exponential recovery lasting 1-3 min. The sawtooth waveforms display almost identical shapes, suggesting a process involving the repeated activation of a fixed source. The particle motion associated with each sawtooth is almost linear, and its major swing shows compressional motion at all stations. Analyses of semblance and particle motion are consistent with a point source located I km beneath the northeast edge of the Halemaumau pit crater. To estimate the source mechanism, we apply a moment tensor inversion to the waveform data, assuming a point source embedded in a homogeneous half-space with compressional and shear wave velocities representative of the average medium properties at shallow depth under Kilauea.

[1] The current eruption at Mount St. Helens is characterized by dome building and shallow, repetitive, long-period (LP) earthquakes. Waveform cross-correlation reveals remarkable similarity for a majority of the earthquakes over periods of several weeks. Stacked spectra of these events display multiple peaks between 0.5 and 2 Hz that are common to most stations. Lower-amplitude very-long-period (VLP) events commonly accompany the LP events. We model the source mechanisms of LP and VLP events in the 0.5-4 s and 8-40 s bands, respectively, using data recorded in July 2005 with a 19-station temporary broadband network. The source mechanism of the LP events includes: 1) a volumetric component modeled as resonance of a gently NNW-dipping, steam-filled crack located directly beneath the actively extruding part of the new dome and within 100 m of the crater floor and 2) a vertical single force attributed to movement of the overlying dome. The VLP source, which also includes volumetric and single-force components, is 250 m deeper and NNW of the LP source, at the SW edge of the 1980s lava dome. The volumetric component points to the compression and expansion of a shallow, magmafilled sill, which is subparallel to the hydrothermal crack imaged at the LP source, coupled with a smaller component of expansion and compression of a dike. The single-force components are due to mass advection in the magma conduit. The location, geometry and timing of the sources suggest the VLP and LP events are caused by perturbations of a common crack system.

) summit eruptive activity at Kilauea Volcano is characterized by explosive degassing bursts accompanied by very long period (VLP) seismic signals. We model the source mechanisms of VLP signals in the 10-50 s band using data recorded for 15 bursts with a 10-station broadband network deployed in the summit caldera. To determine the source centroid location and source mechanism, we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account. The VLP signals associated with the bursts originate in a source region ∼1 km below the eastern perimeter of Halemaumau pit crater. The observed waveforms are well explained by the combination of a volumetric component and a vertical single force component. For the volumetric component, several source geometries are obtained which equally explain the observed waveforms. These geometries include (1) a pipe dipping 64°t o the northeast; (2) two intersecting cracks including an east striking crack (dike) dipping 80°to the north, intersecting a north striking crack (another dike) dipping 65°to the east; (3) a pipe dipping 58°to the northeast, intersecting a crack dipping 48°to the westsouthwest; and (4) a pipe dipping 57°to the northeast, intersecting a pipe dipping 58°to the west-southwest. Using the dual-crack model as reference, the largest volume change obtained among the 15 bursts is ∼24,400 m 3 , and the maximum amplitude (peak to peak) of the force is ∼20 GN. Each burst is marked by a similar sequence of deflation and inflation, trailed by decaying oscillations of the volumetric source. The vertical force is initially upward, synchronous with source deflation, then downward, synchronous with source reinflation, followed by oscillations with polarity opposite to the volumetric oscillations. This combination of force and volume change is attributed to pressure and momentum changes induced during a fluid dynamic source mechanism involving the ascent, expansion, and burst of a large slug of gas within the upper ∼150 m of the magma conduit. As the slug expands upon approach to the surface and more liquid becomes wall supported by viscous shear forces, the pressure below the slug decreases, inducing conduit deflation and an upward force on the Earth. The final rapid slug expansion and burst stimulate VLP and LP oscillations of the conduit system, which slowly decay due to viscous dissipation and elastic radiation. Consideration of the fluid dynamic arguments leads us to prefer the dual-crack VLP source model as it is the only candidate model capable of producing plausible values of length scales and pressure changes. The magnitudes of the vertical forces observed in the 15 bursts appear consistent with slug masses of 10 4 to 10 6 kg.Citation: Chouet, B. A., P. B. Dawson, M. R. James, and S. J. Lane (2010), Seismic source mechanism of degassing bursts at Kilauea Volcano, Hawaii: Results from waveform inversion in the 10-50 s band,

The three-dimensional P and S wave structure of Redoubt Volcano, Alaska, and the underlying crust to depths of 7-8 km is determined from 6219 P wave and 4008 S wave first-arrival times recorded by a 30-station seismograph network deployed on and around the volcano. First-arrival times are calculated using a finite-difference technique, which allows for flexible parameterization of the slowness model and easy inclusion of topography and source-receiver geometry. The three-dimensional P wave velocity structure and hypocenters are determined simultaneously, while the three-dimensional S wave velocity model is determined using the relocated seismicity and an initial S wave velocity model derived from the P wave velocity model assuming an average Vp/Vs ratio of 1.78. Convergence is steady with approximately 73% and 52% reduction in P and S wavearrival time RMS, respectively, after 10 iterations. The most prominent feature observed in the three-dimensional velocity models derived for both P and S waves is a relative lowvelocity, near-vertical, pipelike structure approximately 1 km in diameter that extends from 1 to 6 km beneath sea level. This feature aligns axially with the bulk of seismicity and is interpreted as a highly fractured and altered zone encompassing a magma conduit. The velocity structure beneath the north flank of the volcano between depths of 1 and 6 km is characterized by large lateral velocity variations. High velocities within this region are interpreted as remnant dikes and sills and low velocities as regions along which magma migrates. No large low-velocity body suggestive of a magma chamber is resolved inthe upper 7-8 km of the crust. Introduction Redoubt Volcano, Alaska, one of several active Quaternary volcanoes lying in south-central Alaska at the eastern end of the Aleutian volcanic arc (Figure la), erupted more than 20 times between December 13, 1989, and April 21, 1990. As part of its investigation of this eruption sequence the United States Geological Survey (USGS) deployed a network of 21 portable three-component seismographs in the vicinity of Redoubt Volcano for 3 weeks during July 1991 (Figure lb). The portable seismic network, complemented by a nine-station regional seismic network operated by the Alaskan Volcano Observatory (AVO), recorded thousands of volcano-tectonic (VT) earthquakes that are used to determine the three-dimensional P and S wave velocity structure of Redoubt Volcano and to gain a better understanding of the dynamics of magma migration, emplacement, and eruption.

No abstract

[1] The source mechanism of very long period (VLP) signals accompanying volcanic degassing bursts at Popocatépetl is analyzed in the 15-70 s band by minimizing the residual error between data and synthetics calculated for a point source embedded in a homogeneous medium. The waveforms of two eruptions (23 April and 23 May 2000) representative of mild Vulcanian activity are well reproduced by our inversion, which takes into account volcano topography. The source centroid is positioned 1500 m below the western perimeter of the summit crater, and the modeled source is composed of a shallow dipping crack (sill with easterly dip of 10°) intersecting a steeply dipping crack (northeast striking dike dipping 83°northwest), whose surface extension bisects the vent. Both cracks undergo a similar sequence of inflation, deflation, and reinflation, reflecting a cycle of pressurization, depressurization, and repressurization within a time interval of 3-5 min. The largest moment release occurs in the sill, showing a maximum volume change of 500-1000 m 3 , pressure drop of 3-5 MPa, and amplitude of recovered pressure equal to 1.2 times the amplitude of the pressure drop. In contrast, the maximum volume change in the dike is less (200-300 m 3 ), with a corresponding pressure drop of 1-2 MPa and pressure recovery equal to the pressure drop. Accompanying these volumetric sources are single-force components with magnitudes of 10 8 N, consistent with melt advection in response to pressure transients. The source time histories of the volumetric components of the source indicate that significant mass movement starts within the sill and triggers a mass movement response in the dike within a few seconds. Such source behavior is consistent with the opening of a pathway for escape of pent-up gases from slow pressurization of the sill driven by magma crystallization. The opening of this pathway and associated rapid evacuation of volcanic gases induces the pressure drop. Pressure recovery in the magma filling the sill is driven by diffusion of gases from the resulting supersaturated melt into bubbles. Assuming a penny-shaped crack at ambient pressure of 40 MPa, the observed pressure and volume variations can be modeled with the following attributes: crack radius (100 m), crack aperture (5 m), bubble number density (10 10 -10 12 m À3 ), initial bubble radius (10 À6 m), final bubble radius ($10 À5 m), and net decrease of gas concentration in the melt (0.01 wt %).Citation: Chouet, B., P. Dawson, and A. Arciniega-Ceballos (2005), Source mechanism of Vulcanian degassing at Popocatépetl Volcano, Mexico, determined from waveform inversions of very long period signals,

Abstract.High-resolution velocity models (0.5 km resolution) of the Kilauea caldera region are obtained by the tomographic inversion of both P-and S-wave arrival times. Data are from the permanent Hawaiian Volcano Observatory (HVO) seismic network, a broadband seismic network, and a temporary array of stations centered on the southern boundary of the caldera. A low-velocity P-wave anomaly is imaged centered on the southeastern edge of the caldera, with a velocity contrast of about 10% and a volume of 27 km 3. The