Bernard A. Chouet

Source mechanisms of explosions at Stromboli Volcano, Italy, determined from moment‐tensor inversions of very‐long‐period data

[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. INDEX TERMS: 7215 Seismology: Earthquake parameters; 7280 Seismology: Volcano seismology (8419); 8414 Volcanology: Eruption mechanisms; KEYWORDS: very-long-period seismicity, moment tensor inversions, eruption mehanics

Earthquake classification, location, and error analysis in a volcanic environment: implications for the magmatic system of the 1989–1990 eruptions at redoubt volcano, Alaska

Abstract Determination of the precise locations of seismic events associated with the 1989–1990 eruptions of Redoubt Volcano posed a number of problems, including poorly known crustal velocities, a sparse station distribution, and an abundance of events with emergent phase onsets. In addition, the high relief of the volcano could not be incorporated into the hypoellipse earthquake location algorithm. This algorithm was modified to allow hypocenters to be located above the elevation of the seismic stations. The velocity model was calibrated on the basis of a posteruptive seismic survey, in which four chemical explosions were recorded by eight stations of the permanent network supplemented with 20 temporary seismographs deployed on and around the volcanic edifice. The model consists of a stack of homogeneous horizontal layers; setting the top of the model at the summit allows events to be located anywhere within the volcanic edifice. Detailed analysis of hypocentral errors shows that the long-period (LP) events constituting the vigorous 23-hour swarm that preceded the initial eruption on December 14 could have originated from a point 1.4 km below the crater floor. A similar analysis of LP events in the swarm preceding the major eruption on January 2 shows they also could have originated from a point, the location of which is shifted 0.8 km northwest and 0.7 km deeper than the source of the initial swarm. We suggest this shift in LP activity reflects a northward jump in the pathway for magmatic gases caused by the sealing of the initial pathway by magma extrusion during the last half of December. Volcano-tectonic (VT) earthquakes did not occur until after the initial 23-hour-long swarm. They began slowly just below the LP source and their rate of occurrence increased after the eruption of 01:52 AST on December 15, when they shifted to depths of 6 to 10 km. After January 2 the VT activity migrated gradually northward; this migration suggests northward propagating withdrawal of magma from a plexus of dikes and/or sills located in the 6 to 10 km depth range. Precise relocations of selected events prior to January 2 clearly resolve a narrow, steeply dipping, pencil-shaped concentration of activity in the depth range of 1–7 km, which illuminates the conduit along which magma was transported to the surface. A third event type, named hybrid, which blends the characteristics of both VT and LP events, originates just below the LP source, and may reflect brittle failure along a zone intersecting a fluid-filled crack. The distribution of hybrid events is elongated 0.2–0.4 km in an east-west direction. This distribution may offer constraints on the orientation and size of the fluid-filled crack inferred to be the source of the LP events.

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

Abstract 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.

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.

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 1014 N m and 3000 m3, respectively. Our results also point to the existence of a single force component with amplitude of 109 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 m3 in good agreement with the integrated volume flow rate of magma estimated near the eruptive site.

Source and path effects in the wave fields of tremor and explosions at Stromboli volcano, Italy

The wave fields generated by Strombolian activity are investigated using data from small-aperture seismic arrays deployed on the north flank of Stromboli and data from seismic and pressure transducers set up near the summit crater. Measurements of slowness and azimuth as a function of time clearly indicate that the sources of tremor and explosions are located beneath the summit crater at depths shallower than 200 m with occasional bursts of energy originating from sources extending to a depth of 3 km. Slowness, azimuth, and particle motion measurements reveal a complex composition of body and surface waves associated with topography, structure, and source properties. Body waves originating at depths shallower than 200 m dominate the wave field at frequencies of 0.5–2.5 Hz, and surface waves generated by the surficial part of the source and by scattering sources distributed around the island dominate at frequencies above 2.5 Hz. The records of tremor and explosions are both dominated by SH motion. Far-field records from explosions start with radial motion, and near-field records from those events show dominantly horizontal motion and often start with a low-frequency (1–2 Hz) precursor characterized by elliptical particle motion, followed within a few seconds by a high-frequency radial phase (1–10 Hz) accompanying the eruption of pyroclastics. The dominant component of the near- and far-field particle motions from explosions, and the timing of air and body wave phases observed in the near field, are consistent with a gas-piston mechanism operating on a shallow (<200 m deep), vertical crack-like conduit. Models of a degassing fluid column suggest that noise emissions originating in the collective oscillations of bubbles ascending in the magma conduit may provide an adequate self-excitation mechanism for sustained tremor generation at Stromboli.

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 ps are the densities of the fluid and rock matrix, respectively. The quality factor due to radiation loss (Qr) 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 (Qi) appropriate for these types of fluids, in which Q is given by Q−l = Qr−l + Qi−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.

Imaging the crustal magma sources beneath Mauna Loa and Kilauea volcanoes, Hawaii

Three-dimensional seismic P-wave traveltime tomography is used to image the magma sources beneath Mauna Loa and Kilauea volcanoes, Hawaii. High-velocity bodies (>6.4 km/s) in the upper 9 km of the crust beneath the summits and rift zones of the volcanoes correlate with zones of high magnetic intensities and are interpreted as solidified gabbro-ultramafic cumulates from which the surface volcanism is derived. The proximity of these high-velocity features to the rift zones is consistent with a ridge-spreading model of the volcanic flank. Southeast of the Hilina fault zone, along the south flank of Kilauea, low-velocity material (

Long-period events, the most characteristic seismicity accompanying the emplacement and extrusion of a lava dome in Galeras Volcano, Colombia, in 1991

Abstract Since its reactivation in 1988 the principal eruptions of Galeras Volcano occurred on May 4–9, 1989, July 16, 1992, and January 14, March 23, April 3, April 14 and June 7, 1993. The initial eruption was a phreatic event which clearly marked a new period of activity. A lava dome was extruded within the main crater in October 1991 and subsequently destroyed in an explosive eruption on July 16, 1992. The eruptions that followed were all vulcanian-type explosions. The seismicity accompanying the emplacement, extrusion, and destruction of the lava dome was dominated by a mix of long-period (LP) events and tremor displaying a variety of waveforms. Repetitive LP events with dominant periods in the range 0.2–1 s were observed in October and November 1991 and visually correlated with short energetic pulses of gas venting through a crack bisecting the dome surface. Each LP event was characterized by a weak precursory signal with dominant periods in the range 0.05–0.1 s lasting roughly 7 s. Using the fluid-driven crack model of Chouet (1988, 1992), we infer that two distinct cracks may have acted as sources for the LP and precursor signals. Spectral analyses of the data yield the following parameters for the LP source: crack length, 240–360 m; crack width, 130–150 m; crack aperture, 0.5–3.4 mm; crack stiffness, 100–500; sound speed of fluid, 880 m/s; and excess pressure, 0.01–0.19 MPa. Similar analyses yield the parameters of the precursor source: crack length, 20–30 m; crack width, 15–25 m; crack aperture, 2.3–8.7 mm; crack stiffness, 5–15; sound speed of fluid, 140 m/s; and excess pressure, 0.06–0.15 MPa. Combined with geologic and thermodynamic constraints obtained from field observations, these seismic parameters suggest a gas-release mechanism in which the episodic collapse of a foam layer trapped at the top of the magma column subjacent to the dome releases a slug of pressurized gas which escapes to the surface while dilating a preexisting system of cracks in the dome structure. Accordingly, the fracture observed on the crystallized dome body is the surface extension of the LP-source crack, where LP activity is induced by the rapid emission and expansion of gas flowing through this conduit. The width and aperture of the crack estimated in the model are in good agreement with the length and aperture of the fracture estimated from visual observations. The source parameters of the precursor signal are suggestive of a nozzle-like conduit connecting the LP-source crack to the underlying magma reservoir. Excitation of this conduit segment is attributed to the rapid emission and acceleration of the frothy fluid resulting from the collapse of the foam layer at the top of the reservoir. The calculated periodicity of foam collapse events is in agreement with the observed average rate of thirteen LP events per hour.

Characteristics of seismic waves composing Hawaiian volcanic tremor and gas‐piston events observed by a near‐source array

A correlation method, specifically designed for describing the characteristics of a complex wave field, is applied to volcanic tremor and gas-piston events recorded by a semicircular array of GEOS instruments set at the foot of the Puu Oo crater on the east rift of Kilauea volcano, Hawaii. The spatial patterns of correlation coefficients obtained as functions of frequency for the three components of motion over the entire array are similar for gas-piston events and tremor, and clearly depict dispersive waves propagating across the array from the direction of Puu Oo. The wave fields are composed of comparable amounts of Rayleigh and Love waves propagating with similar and extremely slow phase velocities ranging from 700 m/s at 2 Hz to 300 m/s at 8 Hz. The highly cracked solidified lava flow on which the array was deployed, and subjacent structure of alternating lava and ash layers formed during repeated eruptions of Puu Oo since 1983, appear to be responsible for the low velocities observed. The results from Puu Oo stand in sharp contrast to those obtained in an experiment conducted in 1976 on the partially solidified lava lake of Kilauea Iki. Rayleigh waves were not observed in Kilauea Iki, but well-developed trains of Love waves were seen to propagate there with velocities twice as high as those observed near Puu Oo. These differences in the propagation characteristics of surface waves at the two sites may be attributed to the presence of a soft horizontal layer of molten rock in Kilauea Iki, which may have lowered the phase velocity of Rayleigh waves more drastically than that of Love waves, resulting in severe scattering of the Rayleigh wave mode. On the other hand, the thin superficial pahoehoe flow under our array at Puu Oo may have favored the development of vertical columnar joints more extensively at this location than at Kilauea Iki, which may have reduced the shear moduli controlling the Love wave mode. The average phase velocities in the frequency band from 2 to 5 Hz found at Puu Oo are roughly similar to those determined for tremor at Klyuchevskoy volcano, Kamchatka, but the frequency dependence appears much stronger at Puu Oo.