Search USGS⌕ Search

Geology topics

B. A. Chouet

Publications and source records attributed to B. A. Chouet.

At least 19 recordsLinked to original sources
USGS2009

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.

Washington↗
USGS2008

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.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Gas slug ascent through changes in conduit diameter: Laboratory insights into a volcano-seismic source process in low-viscosity magmas

Seismic signals generated during the flow and degassing of low-viscosity magmas include long-period (LP) and very-long-period (VLP) events, whose sources are often attributed to dynamic fluid processes within the conduit. We present the results of laboratory experiments designed to investigate whether the passage of a gas slug through regions of changing conduit diameter could act as a suitable source mechanism. A vertical, liquid-filled glass tube featuring a concentric diameter change was used to provide canonical insights into potentially deep or shallow seismic sources. As gas slugs ascend the tube, we observe systematic pressure changes varying with slug size, liquid depth, tube diameter, and liquid viscosity. Gas slugs undergoing an abrupt flow pattern change upon entering a section of significantly increased tube diameter induce a transient pressure decrease in and above the flare and an associated pressure increase below it, which stimulates acoustic and inertial resonant oscillations. When the liquid flow is not dominantly controlled by viscosity, net vertical forces on the apparatus are also detected. The net force is a function of the magnitude of the pressure transients generated and the tube geometry, which dictates where, and hence when, the traveling pressure pulses can couple into the tube. In contrast to interpretations of related volcano-seismic data, where a single downward force is assumed to result from an upward acceleration of the center of mass in the conduit, our experiments suggest that significant downward forces can result from the rapid deceleration of relatively small volumes of downward-moving liquid. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Source process of a long-period event at Kilauea volcano, Hawaii

We analyse a long-period (LP) event observed by a dense seismic network temporarily operated at Kilauea volcano, Hawaii, in 1996. We systematically perform spectral analyses, waveform inversions and forward modeling of the LP event to quantify its source process. Spectral analyses identify two dominant spectral frequencies at 0.6 and 1.3 Hz with associated Q values in the range 10-20. Results from waveform inversions assuming six moment-tensor and three single-force components point to the resonance of a horizontal crack located at a depth of approximately 150 m near the northeastern rim of the Halemaumau pit crater. Waveform simulations based on a fluid-filled crack model suggest that the observed frequencies and Q values can be explained by a crack filled with a hydrothermal fluid in the form of either bubbly water or steam. The shallow hydrothermal crack located directly above the magma conduit may have been heated by volcanic gases leaking from the conduit. The enhanced flux of heat raised the overall pressure of the hydrothermal fluid in the crack and induced a rapid discharge of fluid from the crack, which triggered the acoustic vibrations of the resonator generating the LP waveform. The present study provides further support to the idea that LP events originate in the resonance of a crack. 2005 RAS.

Hawaii↗
USGS2003

Source mechanism of long-period events at Kusatsu-Shirane Volcano, Japan, inferred from waveform inversion of the effective excitation functions

We investigate the source mechanism of long-period (LP) events observed at Kusatsu-Shirane Volcano, Japan, based on waveform inversions of their effective excitation functions. The effective excitation function, which represents the apparent excitation observed at individual receivers, is estimated by applying an autoregressive filter to the LP waveform. Assuming a point source, we apply this method to seven LP events the waveforms of which are characterized by simple decaying and nearly monochromatic oscillations with frequency in the range 1-3 Hz. The results of the waveform inversions show dominant volumetric change components accompanied by single force components, common to all the events analyzed, and suggesting a repeated activation of a sub-horizontal crack located 300 m beneath the summit crater lakes. Based on these results, we propose a model of the source process of LP seismicity, in which a gradual buildup of steam pressure in a hydrothermal crack in response to magmatic heat causes repeated discharges of steam from the crack. The rapid discharge of fluid causes the collapse of the fluid-filled crack and excites acoustic oscillations of the crack, which produce the characteristic waveforms observed in the LP events. The presence of a single force synchronous with the collapse of the crack is interpreted as the release of gravitational energy that occurs as the slug of steam ejected from the crack ascends toward the surface and is replaced by cooler water flowing downward in a fluid-filled conduit linking the crack and the base of the crater lake. ?? 2003 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS2002

Waveform inversion of oscillatory signatures in long-period events beneath volcanoes

The source mechanism of long-period (LP) events is examined using synthetic waveforms generated by the acoustic resonance of a fluid-filled crack. We perform a series of numerical tests in which the oscillatory signatures of synthetic LP waveforms are used to determine the source time functions of the six moment tensor components from waveform inversions assuming a point source. The results indicate that the moment tensor representation is valid for the odd modes of crack resonance with wavelengths 2 L / n , 2 W / n , n = 3, 5, 7, …, where L and W are the crack length and width, respectively. For the even modes with wavelengths 2 L / n , 2 W / n , n = 2, 4, 6, …, a generalized source representation using higher-order tensors is required, although the efficiency of seismic waves radiated by the even modes is expected to be small. We apply the moment tensor inversion to the oscillatory signatures of an LP event observed at Kusatsu-Shirane Volcano, central Japan. Our results point to the resonance of a subhorizontal crack located a few hundred meters beneath the summit crater lakes. The present approach may be useful to quantify the source location, geometry, and force system of LP events, and opens the way for moment tensor inversions of tremor.

Journal of Geophysical Research B: Solid Earth↗
USGS2001

Trends in long-period seismicity related to magmatic fluid compositions

Sound speeds and densities are calculated for three different types of fluids: gas-gas mixture; ash-gas mixture; and bubbly liquid. These fluid properties are used to calculate the impedance contrast (Z) and crack stiffness (C) in the fluid-driven crack model (Chouet: J. Geophys. Res., 91 (1986) 13,967; 101 (1988) 4375; A seismic model for the source of long-period events and harmonic tremor. In: Gasparini, P., Scarpa, R., Aki, K. (Eds.), Volcanic Seismology, IAVCEI Proceedings in Volcanology, Springer, Berlin, 3133). The fluid-driven crack model describes the far-field spectra of long-period (LP) events as modes of resonance of the crack. Results from our calculations demonstrate that ash-laden gas mixtures have fluid to solid density ratios comparable to, and fluid to solid velocity ratios lower than bubbly liquids (gas-volume fractions <10%). This difference results in synthetic far-field spectra with higher impedance contrasts and narrower spectral bandwidths for ash-laden gas mixture than spectra for bubbly liquids. Spectral characteristics are described in terms of the quality factor Q-1. Q-1 is measured by the ratio of the frequency of the dominant spectral peak to the bandwidth of the peak measured at one half of its amplitude. This factor expresses the losses of energy due to elastic radiation Q-1r and other dissipative mechanisms Q-1i at the source, Q-1 = Q-1r + Q-1i. Spectra for LP events recorded at active volcanoes such as Galeras in Colombia and Kilauea in Hawaii, have Q-1 factors in the range of 0.1-0.002. The Q-1r factors due to radiation loss calculated for a sphere filled with a H2O-CO2 or H2O-SO2 gas mixture, vary between 0.0015 and 0.0040 with a change in wt% H2O at 800-1600 K and 10-50 MPa. For gas-rich mixtures, Q-1r has a strong dependence on resonator geometry (spherical versus rectangular). The spectra from a resonating sphere filled with gas-rich mixture yields values of Q-1r an order of magnitude smaller than those from a rectangular crack. For a resonating crack filled with an ash-gas mixture (or pseudogas), Q-1r varies parabolically from ???0.006 for an ash-rich mixture, to 0.0015 or 0.0023 for a H2O-rich or CO2-rich mixture at 800 K and 25 MPa. For low (<20%) gas-volume fraction fluids (foams, bubbly fluids and ash-rich pseudogases), the magnitudes for Q-1r are independent of crack geometry. Spectra associated with a foam (gas-volume fractions 10-90%) or bubbly basalt (gas-volume fractions <10%) may have a dominant spectral peak with values of Q-1r on the order of 0.01 and 0.1, respectively. The spectra from a resonating sphere filled with a foam containing >20% gas-volume fraction yields values of Q-1r similar to those for a rectangular crack. As with gas-gas and ash-gas mixtures, an increase in mass fraction narrows the bandwidth of the dominant mode and shifts the spectra to lower frequencies. Including energy losses due to dissipative processes in a bubbly liquid increases attenuation. Attenuation may also be higher in ash-gas mixtures and foams if the effects of momentum and mass transfer between the phases were considered in the calculations. ?? 2001 Elsevier Science B. V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS2001

Evolution of the December 14, 1989 precursory long-period event swarm at Redoubt volcano, Alaska

The intermittency pattern and evolution in waveforms of long-period (LP) seismic events during the intense, 23-h swarm that preceded the December 14, 1989 eruption of Redoubt volcano are investigated. Utilizing cross correlation to exploit the high degree of similarity among waveforms, a substantially more complete event catalog is generated than was available from near realtime detection based on short-term/long-term amplitude ratios, which was saturated by the high rate of activity. The temporal magnitude distribution of the predominant LP events is found to have an unusual banded structure in which the average magnitude of each band slowly increases and then decreases through time. A bifurcation that appears in the uppermost band shortly after the peak in magnitudes is characterized by a quasi-periodicity in intermittency and magnitude that is reminiscent of one of the classic routes to chaotic behavior in some non-linear systems. The waveforms of the predominant events evolve slowly but unsteadily through time. These gradual changes appear to result from variations in the relative amplitudes of spectral peaks that remain stable in frequency, which suggests that they are due to differential excitation of a single, resonant source. Two other previously unrecognized, repetitive waveforms are also identified, but the signals from these secondary events are not clearly recorded at distances beyond the closest station. Similarities among the spectra of the predominant and secondary events suggest that the signals from these events also could represent different modes of exciting the same source. Significant changes in the rates and the sizes of the largest of these secondary events appear to coincide with the peak in the size distribution of the predominant LPs. At least some of the non-repetitive LP waveforms in the swarm appear to be the result of the superposition of signals from the rapid repetition of predominant LP source, thus placing a constraint on the repeat time of the triggering mechanism for this source. A lone hybrid event, which has a waveform character intermediate between the predominant LP events and high-frequency volcano-tectonic events, was also identified in the swarm; the occurrence of this event provides important evidence that the low-frequency character of the LP events is a source rather than a path or site effect. ?? 2001 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS2001

The dependence of acoustic properties of a crack on the resonance mode and geometry

We examine the dependence of the acoustic properties of a crack containing magmatic or hydrothermal fluids on the resonance mode and geometry to quantify the source properties of long-period (LP) events observed in volcanic areas. Our results, based on spectral analyses of synthetic waveforms generated with a fluid-driven crack model, indicate that the basic features of the dimensionless frequency (??) and quality factor (Qr) for a crack containing various types of fluids are not strongly affected by the choice of mode, although the actual ranges of Q?? and ?? both depend on the mode. The dimensionless complex frequency systematically varies with changes in the crack geometry, showing increases in both Qr and ?? as the crack length to aperture ratio decreases. The present results may be useful for the interpretation of spatial and temporal variations in the observed complex frequencies of LP events.

Geophysical Research Letters↗
USGS2000

Acoustic properties of a crack containing magmatic or hydrothermal fluids

We estimate the acoustic properties of a crack containing magmatic or hydrothermal fluids to quantify the source properties of long-period (LP) events observed in volcanic areas assuming that a crack-like structure is the source of LP events. The tails of synthetic waveforms obtained from a model of a fluid-driven crack are analyzed by the Sompi method to determine the complex frequencies of one of the modes of crack resonance over a wide range of the model parameters α/ a and ρ f /ρ s , where αis the P wave velocity of the rock matrix, a is the sound speed of the fluid, and ρ f and p s are the densities of the fluid and rock matrix, respectively. The quality factor due to radiation loss ( Q r ) for the selected mode almost monotonically increases with increasing α/ a , while the dimensionless frequency ( v ) of the mode decreases with increasing α/ a and ρ f /ρ s . These results are used to estimate Q and v for a crack containing various types of fluids (gas-gas mixtures, liquid-gas mixtures, and dusty and misty gases) for values of a , ρ f , and quality factor due to intrinsic losses ( Q i ) appropriate for these types of fluids, in which Q is given by Q −l = Q r −l + Q i −1 . For a crack containing such fluids, we obtain Q ranging from almost unity to several hundred, which consistently explains the wide variety of quality factors measured in LP events observed at various volcanoes. We underscore the importance of dusty and misty gases containing small-size particles with radii around 1 μm to explain long-lasting oscillations with Q significantly larger than 100. Our results may provide a basis for the interpretation of spatial and temporal variations in the observed complex frequencies of LP events in terms of fluid compositions beneath volcanoes.

Journal of Geophysical Research B: Solid Earth↗
USGS1999

Three-dimensional velocity structure of the Kilauea caldera, Hawaii

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³. The V P /V S model mimics the spatial extent of the P-wave anomaly, but is partitioned into two discrete anomalous volumes centered on the southern boundary of the caldera and on the upper east rift of the volcano. The corresponding Poisson's ratio in these zones is high (ν=0.25–0.32) which is consistent with a densely-cracked, hot volume which may contain partial melt. The large-scale features of the models are consistent with results obtained from an earlier, larger-scale (2 km resolution) tomographic image of Kilauea Volcano based on HVO network data.

Geophysical Research Letters↗
USGS1999

The complex frequencies of long-period seismic events as probes of fluid composition beneath volcanoes

Long-period (LP) events have been widely observed in relation to magmatic and hydrothermal activities in volcanic areas. LP waveforms characterized by their harmonic signature have been interpreted as oscillations of a fluid-filled resonator, and mixtures of liquid and gas in the form of bubbly liquids have been mainly assumed for the fluid. To investigate the characteristic properties of the resonator system, we analyse waveforms of LP events observed at four different volcanoes in Hawaii, Alaska, Colombia and Japan using a newly developed spectral method. This method allows an estimation of the complex frequencies of decaying sinusoids based on an autoregressive model. The results of our analysis show a wide variety of Q factors, ranging from tens to several hundred. We compare these complex frequencies with those predicted by the fluid-filled crack model for various mixtures of liquid, gas and ash. Although the oscillations of LP events with Q smaller than 50 can be explained by various combinations of liquids and gases, we find that ash-laden gases are required to explain long-lasting oscillations with Q larger than 100. The complex frequencies of LP events yield useful information on the types of fluids. Temporal and spatial variations of the complex frequencies can be used as probes of fluid composition beneath volcanoes.

Geophysical Journal International↗
USGS1998

Waveform inversion of very long period impulsive signals associated with magmatic injection beneath Kilauea Volcano, Hawaii

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 μrad 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 1 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. Synthetic waveforms are constructed by a superposition of impulse responses for six moment tensor components and three single force components. The origin times of individual impulses are distributed along the time axis at appropriately small, equal intervals, and their amplitudes are determined by least squares. In this inversion, the source time functions of the six tensor and three force components are determined simultaneously. We confirm the accuracy of the inversion method through a series of numerical tests. The results from the inversion show that the waveform data are well explained by a pulsating transport mechanism operating on a subhorizontal crack linking the summit reservoir to the east rift of Kilauea. The crack acts like a buffer in which a batch of fluid (magma and/or gas) accumulates over a period of 1–3 min before being rapidly injected into a larger reservoir (possibly the east rift) over a timescale of 5–10 s. The seismic moment and volume change associated with a typical batch of fluid are approximately 10 14 N m and 3000 m 3 , respectively. Our results also point to the existence of a single force component with amplitude of 10 9 N, which may be explained as the drag force generated by the flow of viscous magma through a narrow constriction in the flow path. The total volume of magma associated with the 4.5-hour-long activation of the pulsating source is roughly 500,000 m 3 in good agreement with the integrated volume flow rate of magma estimated near the eruptive site.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

A numerical investigation of choked flow dynamics and its application to the triggering mechanism of long-period events at Redoubt Volcano, Alaska

We use numerical simulations of transonic flow through a crack to study the dynamics of the formation of shock waves downstream from a nozzle-like constriction inside the crack. The model solves the full set of Navier-Stokes equations in two dimensions via an explicit multifield finite difference representation. The crack walls are assumed to be perfectly rigid, and elastic coupling to the solid is not considered. The simulations demonstrate how the behavior of unsteady shock waves near the walls can produce recurring step-like pressure transients in the flow, which in turn induce resonance of the fluid-filled crack. The motion of the shock waves is governed primarily by smooth, low-amplitude pressure fluctuations at the outlet of the crack. The force induced on the walls scales with the amplitude of the shock, which is a function of the magnitude of the inlet pressure, aperture of the constriction, and thickness of the boundary layer. The applied force also scales in proportion to the spatial extent of the shock excursion, which depends on the fluctuation rate of outlet pressure. Using the source parameters of long-period (LP) events at Redoubt Volcano, Alaska, as a guide for our simulations, we infer that coupling of the shock to the walls occurs for crack inlet to outlet pressure ratios p i /p o >2.31 and that the position of the shock front becomes most sensitive to outlet pressure fluctuations for flow regimes with p i /p o >2.48. For such regimes, fluctuations of outlet pressure of up to ±0.5 MPa at rates up to 3 MPa/s are sufficient to induce pressure transients with magnitudes up to 12.5 MPa over 0.1–2.5 m of the walls within ∼0.5 s. These flow parameters may be adequate for triggering the LP events in the precursory swarm to the December 14, 1989, eruption of Redoubt. According to the flow model the recurrence rate and amplitudes of L.P events are inferred to be a manifestation of the response of a shallow hydrothermal reservoir to the sustained injection of superheated steam from a magma column roofing below this reservoir.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Three-dimensional P and S wave velocity structure of Redoubt Volcano, Alaska

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 wave arrival 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 low‐velocity, 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 in the upper 7–8 km of the crust.

Journal of Geophysical Research B: Solid Earth↗
USGS1994

The 1989-1990 eruptions of Redoubt Volcano: an introduction

Redoubt Volcano, located on the west side of Cook Inlet in south-central Alaska, erupted explosively on over 20 separate occasions between December 14, 1989 and April 21, 1990. Fourteen lava domes were emplaced in the summit area, thirteen of which were subsequently destroyed. The eruption caused economic losses estimated at over $160,000,000 making this the second most costly eruption in U.S. history. This economic impact provided the impetus for a integrated comprehensive account of an erupting volcano using both modern and classical research and modern techniques which in turn led to advances in eruption monitoring and interpretation. Research on such topics as dome formation and collapse and the resulting pyroclastic flows, elutriated ash, lightning, tephra, and flooding was blended with the rapid communication of associated hazards to a large user group. The seismology successes in predicting and monitoring eruption dynamics were due in part to (1) the recognition of long-period seismic events as indicators of the readiness of the volcano to erupt, and (2) to the development of new tools that allowed the seismicity to be assessed instantaneously. Integrated studies of the petrology of erupted products and volatile content over time gave clues as to the progress of the eruption towards completion.

Alaska↗
USGS1994

Precursory swarms of long-period events at Redoubt Volcano (1989-1990), Alaska: Their origin and use as a forecasting tool

During the eruption of Redoubt Volcano from December 1989 through April 1990, the Alaska Volcano Observatory issued advance warnings of several tephra eruptions based on changes in seismic activity related to the occurrence of precursory swarms of long-period (LP) seismic events (dominant period of about 0.5 s). The initial eruption on December 14 occurred after 23 years of quiescence and was heralded by a 23-hour swarm of LP events that ended abruptly with the eruption. After a series of vent-clearing explosions over the next few days, dome growth began on December 21. Another swarm, with LP events similar to those of the first, began on the 26th and ended in a major tephra eruption on January 2. Eruptions continued over the next two weeks and then ceased until February 15, when a large eruption initiated a long phase of repetitive dome-building and dome-destroying episodes that continued into April. Warnings were issued before the major events on December 14 and January 2, but as the eruptive sequence continued after January 2, the energy of the swarms decreased and forecasting became more difficult. A significant but less intense swarm preceded the February 15 eruption, which was not forecast. This eruption destroyed the only seismograph on the volcanic edifice and stymied forecasting until March 4, when the first of three new stations was installed within 3 km of the active vent. From March 4 to the end of the sequence on April 21, there were eight eruptions, six of which were preceded by detectable swarms of LP events. Although weak, these swarms provided the basis for warnings issued before the eruptions on March 23 and April 6. The initial swarm on December 13 had the following features: (1) short duration (23 hours); (2) a rapidly accelerating rate of seismic energy release over the first 18 hours of the swarm, followed by a decline of activity during the 5 hours preceding the eruption; (3) a magnitude range from −0.4 to 1.6; (4) nearly identical LP signatures with a dominant period near 0.5 s; (5) dilatational first motions everywhere; and (6) a stationary source location at a depth of 1.4 km beneath the crater. This occurrence of long-period events suggests a model involving the interaction of magma with groundwater in which magmatic gases, steam and water drive a fixed conduit at a stationary point throughout the swarm. The initiation of that sequence of events is analogous to the failure of a pressure-relief valve connecting a lower, supercharged magma-dominated reservoir to a shallow hydrothermal system. A three-dimensional model of a vibrating fluid-filled crack recently developed by Chouet is found to be compatible with the seismic data and yields the following parameters for the LP source: crack length, 280–380 m; crack width, 140–190 m; crack thickness, 0.05–0.20 m; crack stiffness, 100–200; sound speed of fluid, 0.8–1.3 km/s; compressional-wave speed of rock, 5.1 km/s; density ratio of fluid to rock, ≈0.4; and ratio of bulk modulus of fluid to rigidity of rock, 0.03–0.07. The fluid-filled crack is excited intermittently by an impulsive pressure drop that varies in magnitude within the range of 0.4 to 40 bar. Such disturbance appears to be consistent with a triggering mechanism associated with choked flow conditions in the crack.

Alaska↗