Source Mechanism of Small Long-Period Events at Mount St. Helens in July 2005 Using Template Matching, Phase-Weighted Stacking, and Full-Waveform Inversion

Long-period (LP, 0.5-5 Hz) seismicity, observed at volcanoes worldwide, is a recognized signature of unrest and eruption. Cyclic LP “drumbeating” was the characteristic seismicity accompanying the sustained dome-building phase of the 2004–2008 eruption of Mount St. Helens (MSH), WA. However, together with the LP drumbeating was a near-continuous, randomly occurring series of tiny LP seismic events (LP “subevents”), which may hold important additional information on the mechanism of seismogenesis at restless volcanoes. We employ template matching, phase-weighted stacking, and full-waveform inversion to image the source mechanism of one multiplet of these LP subevents at MSH in July 2005. The signal-to-noise ratios of the individual events are too low to produce reliable waveform inversion results, but the events are repetitive and can be stacked. We apply network-based template matching to 8 days of continuous velocity waveform data from 29 June to 7 July 2005 using a master event to detect 822 network triggers. We stack waveforms for 359 high-quality triggers at each station and component, using a combination of linear and phase-weighted stacking to produce clean stacks for use in waveform inversion. The derived source mechanism points to the volumetric oscillation (∼10 m3) of a subhorizontal crack located at shallow depth (∼30 m) in an area to the south of Crater Glacier in the southern portion of the breached MSH crater. A possible excitation mechanism is the sudden condensation of metastable steam from a shallow pressurized hydrothermal system as it encounters cool meteoric water in the outer parts of the edifice, perhaps supplied from snow melt

Seismic source mechanism of degassing bursts at Kilauea Volcano, Hawaii:results from waveform inversion in the 10–50 s band

The current (March 2008 to February 2009) 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° to 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 west–southwest; 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 m3, 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 104 to 106 kg

Eruption dynamics at Mount St. Helens imaged from broadband seismic waveforms: Interaction of the shallow magmatic and hydrothermal systems

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, magma-filled 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

Experimental observations of pressure oscillations and flow regimes in an analogue volcanic system.

Gas-liquid flows, designed to be analogous to those in volcanic conduits, are generated in the laboratory using organic gas-gum rosin mixtures expanding in a vertically mounted tube. The expanding fluid shows a range of both flow and pressure oscillation behaviors. Weakly supersaturated source liquids produce a low Reynolds number flow with foam expanding from the top surface of a liquid that exhibits zero fluid velocity at the tube wall; i.e., the conventional “no-slip” boundary condition. Pressure oscillations, often with strong long-period characteristics and consistent with longitudinal and radial resonant oscillation modes, are detected in these fluids. Strongly supersaturated source liquids generate more energetic flows that display a number of flow regimes. These regimes include a static liquid source, viscous flow, detached flow (comprising gas-pockets-at-wall and foam-in-gas annular flow, therefore demonstrating strong radial heterogeneity), and a fully turbulent transonic fragmented or mist flow. Each of these flow regimes displays characteristic pressure oscillations that can be related to resonance of flow features or wall impact phenomena. The pressure oscillations are produced by the degassing processes without the need of elastic coupling to the confining medium or flow restrictors and valvelike features. The oscillatory behavior of the experimental flows is compared to seismoacoustic data from a range of volcanoes where resonant oscillation of the fluid within the conduit is also often invoked as controlling the observed oscillation frequencies. On the basis of the experimental data we postulate on the nature of seismic signals that may be measured during large-scale explosive activity

The source of infrasound associated with long-period events at mount St. Helens

During the early stages of the 2004-2008 Mount St. Helens eruption, the source process that produced a sustained sequence of repetitive long-period (LP) seismic events also produced impulsive broadband infrasonic signals in the atmosphere. To assess whether the signals could be generated simply by seismic-acoustic coupling from the shallow LP events, we perform finite difference simulation of the seismo-acoustic wavefield using a single numerical scheme for the elastic ground and atmosphere. The effects of topography, velocity structure, wind, and source configuration are considered. The simulations show that a shallow source buried in a homogeneous elastic solid produces a complex wave train in the atmosphere consisting of P/SV and Rayleigh wave energy converted locally along the propagation path, and acoustic energy originating from , the source epicenter. Although the horizontal acoustic velocity of the latter is consistent with our data, the modeled amplitude ratios of pressure to vertical seismic velocity are too low in comparison with observations, and the characteristic differences in seismic and acoustic waveforms and spectra cannot be reproduced from a common point source. The observations therefore require a more complex source process in which the infrasonic signals are a record of only the broadband pressure excitation mechanism of the seismic LP events. The observations and numerical results can be explained by a model involving the repeated rapid pressure loss from a hydrothermal crack by venting into a shallow layer of loosely consolidated, highly permeable material. Heating by magmatic activity causes pressure to rise, periodically reaching the pressure threshold for rupture of the valve sealing the crack. Sudden opening of the valve generates the broadband infrasonic signal and simultaneously triggers the collapse of the crack, initiating resonance of the remaining fluid. Subtle waveform and amplitude variability of the infrasonic signals as recorded at an array 13.4 km to the NW of the volcano are attributed primarily to atmospheric boundary layer propagation effects, superimposed upon amplitude changes at the source. Copyright 2009 by the American Geophysical Union