Sulfur,chlorine, and fluorine degassing and atmospheric loading by the 1783–1784 AD Laki (Skaftár Fires) eruption in Iceland

 The 1783–1784 Laki tholeiitic basalt fissure eruption in Iceland was one of the greatest atmospheric pollution events of the past 250 years, with widespread effects in the northern hemisphere. The degassing history and volatile budget of this event are determined by measurements of pre-eruption and residual contents of sulfur, chlorine, and fluorine in the products of all phases of the eruption. In fissure eruptions such as Laki, degassing occurs in two stages: by explosive activity or lava fountaining at the vents, and from the lava as it flows away from the vents. Using the measured sulfur concentrations in glass inclusions in phenocrysts and in groundmass glasses of quenched eruption products, we calculate that the total accumulative atmospheric mass loading of sulfur dioxide was 122 Mt over a period of 8 months. This volatile release is sufficient to have generated ∼250 Mt of H₂SO₄ aerosols, an amount which agrees with an independent estimate of the Laki aerosol yield based on atmospheric turbidity measurements. Most of this volatile mass (∼60 wt.%) was released during the first 1.5 months of activity. The measured chlorine and fluorine concentrations in the samples indicate that the atmospheric loading of hydrochloric acid and hydrofluoric acid was ∼7.0 and 15.0 Mt, respectively. Furthermore, ∼75% of the volatile mass dissolved by the Laki magma was released at the vents and carried by eruption columns to altitudes between 6 and 13 km. The high degree of degassing at the vents is attributed to development of a separated two-phase flow in the upper magma conduit, and implies that high-discharge basaltic eruptions such as Laki are able to loft huge quantities of gas to altitudes where the resulting aerosols can reside for months or even 1–2 years. The atmospheric volatile contribution due to subsequent degassing of the Laki lava flow is only 18 wt.% of the total dissolved in the magma, and these emissions were confined to the lowest regions of the troposphere and therefore important only over Iceland. This study indicates that determination of the amount of sulfur degassed from the Laki magma batch by measurements of sulfur in the volcanic products (the petrologic method) yields a result which is sufficient to account for the mass of aerosols estimated by other methods. Received: 30 May 1995 / Accepted: 19 April 1996     

An experimental facility for the investigation of magma fragmentation by rapid decompression

 A special experimental facility has been developed to investigate the fragmentation of vesicular magma undergoing rapid decompression. The facility operates in a regime similar to that of shock tubes and at temperatures up to 950  °C and pressures up to 200 bar. Cylindrical samples (diameter ca. 17 mm, length ca. 50 mm) undergo rapid decompression in a high-temperature, high-pressure section of the facility following the disruption of a diaphragm separating that section from a low-pressure, low-temperature section. Actual vesicular magma samples have been experimentally fragmented at elevated temperatures and pressures corresponding to those observed during explosive volcanic eruptions and the resulting pyroclastics have been photographically resolved in flight and collected for physical characterization. The results of these experiments show that the rapid decompression of highly viscous vesicular magma can generate pyroclastic ejecta via rapid and complete fragmentation of magma at high temperature. This new fragmentation facility is presented as a tool for experimental volcanology under well-constrained conditions. Received: 19 March 1996 / Accepted: 25 August 1996     

Dynamics of carbon dioxide emissions,crystallization, and magma ascent: hypotheses,theory, and applications to volcano monitoring at Mount St. Helens

 Measurements of CO₂ fluxes from open-vent volcanos are rare, yet may offer special capabilities for monitoring volcanos and forecasting activity. The measured fluxes of CO₂ and SO₂ from Mount St. Helens decreased from July through November 1980, but the record includes variations of CO₂/SO₂ in the emitted gas and episodes of greatly increased fluxes of CO₂. We propose that the CO₂ flux variations reflect two gas components: (a) a component whose flux decreased in proportion to 1/ √t with a CO₂/SO₂ mass ratio of 1.7, and (b) a residual flux of CO₂ consisting of short-lived, large peaks with a CO₂/SO₂ mass ratio of 15. We propose two hypotheses: (a) the 1/ √t dependence was generated by crystallization in a deep magma body at rates governed by diffusion-limited heat transfer, and (b) the gas component with the higher CO₂/SO₂ was released from ascending magma, which replenished the same magma body. The separation of the total CO₂ flux into contributions from known processes permits quantitative inferences about the replenishment and crystallization rates of open-system magma bodies beneath volcanos. The flux separations obtained by using two gas sources with distinct CO₂/SO₂ ratios and a peak minus background approach to obtain the CO₂ contributions from an intermittent source and a continuously emitting source are similar. The flux separation results support the hypothesis that the second component was generated by episodic magma ascent and replenishment of the magma body. The diffusion-limited crystallization hypothesis is supported by the decay of minimum CO₂ and SO₂ fluxes with 1/ √t after 1 July 1980. We infer that the magma body at Mount St. Helens was replenished at an average rate (2.8×10⁶m³d^–1) which varied by less than 5% during July, August, and September 1980. The magma body volume (2.4–3.0 km³) in early 1982 was estimated by integrating a crystallization rate function inferred from CO₂ fluxes to maximum times (20±4 years) estimated from the increase of sample crystallinity with time. These new volcanic gas flux separation methods and the existence of relations among the CO₂ flux, crystallization rates, and magma body replenishment rates yield new information about the dynamics of an open-vent, replenished magma body. Received: 15 February 1995 / Accepted: 30 March 1996     

CO2 output and δ13C(CO2) from Mount Etna as indicators of degassing of shallow asthenosphere

 An estimated average CO₂ output from Etna's summit craters in the range of 13±3 Mt/a has recently been determined from the measured SO₂ output and measured CO₂/SO₂ molar ratios. To this amount the CO₂ output emitted diffusely from the soil (≈ 1 Mt/a) and the amount of CO₂ dissolved in Etna's aquifers (≈ 0.25 Mt/a) must be added. Data on the solubility of CO₂ in Etnean magmas at high temperature and pressure allow the volume of magma involved in the release of such an amount of this gas to be estimated. This volume of magma (≈ 0.7 km³/a) is approximately 20 times greater than the volume of magma erupted annually during the period 1971–1995. On the basis of C-isotopic data of CO₂ collected in the Etna area and of new hypotheses on the source of Mediterranean magmas, significant contributions of CO₂ from non-magmatic sources to the total output from Etna are unlikely. Such large outputs of CO₂ and also of SO₂ from Etna could be due to an anomalously shallow asthenosphere beneath the volcano that allows a continuous escape of gases toward the surface, even without migration of magma. Received: 7 August 1996 / Accepted: 9 November 1996     

Fragmentation of magma during Plinian volcanic eruptions

 The ratio of the volume of vesicles (gas) to that of glass (liquid) in pumice clasts (V _G /V _L ) reflects the degassing and dynamic history experienced by a magma during an explosive eruption. V _G /V _L in pumices from a large number of Plinian eruption deposits is shown here to vary by two orders of magnitude, even between pumices at a given level in a deposit. These variations in V _G /V _L do not correlate with crystallinity or initial water content of the magmas or their eruptive intensities, despite large ranges in these variables. Gas volume ratios of pumices do, however, vary systematically with magma viscosity estimated at the point of fragmentation, and we infer that pumices do not quench at the level of fragmentation but undergo some post-fragmentary evolution. On the timescale of Plinian eruptions, pumices with viscosities <10⁹ Pa s can expand after fragmentation, as long as their bubbles retain gas, at a rate inversely proportional to their viscosity. Once the bubbles connect to form a permeable network and lose their gas, expansion halts and pumices with viscosities <10⁵ Pa s can collapse under the action of surface tension. Textural evidence from bubble sizes and shapes in pumices indicates that both expansion and collapse have taken place. The magnitudes of expansion and collapse, therefore, depend critically on the timing of bubble connectivity relative to the final moment of quenching. We propose that bubbles in different pumices become connected at different times throughout the time span between fragmentation and quenching. After accounting for these effects, we derive new information on the fragmentation process from two characteristics of pumices. The most important is a relatively constant minimum value of V _G /V _L of ∼1.78 (64 vol.% vesicularity) in all samples with viscosities >10⁵ Pa s. This value is independent of magma composition and thus reflects a property of the eruptive mechanism. The other characteristic is that highly expanded pumices (>85 vol.% vesicularities) are common, which argues against overpressure in bubbles as a mechanism for fragmenting magma. We suggest that magma fragments when it reaches a vesicularity of ∼64 vol.%, but only if sheared sufficiently strongly. The intensity of shear varies as a function of velocity in the conduit, which is related to overpressure in the chamber, so that changes in overpressure with time are important in controlling the common progression from explosive to effusive activity at volcanoes. Received: 19 April 1995 / Accepted: 3 April 1996     

Multiple levels of magma storage during the 1980 summer eruptions of Mount St. Helens, WA

Transitions in eruptive style—explosive to effusive, sustained to pulsatory—are a common aspect of volcanic activity and present a major challenge to volcano monitoring efforts. A classic example of such transitions is provided by the activity of Mount St. Helens, WA, during 1980, where a climactic Plinian event on May 18 was followed by subplinian and vulcanian eruptions that became increasing pulsatory with time throughout the summer, finally progressing to episodic growth of a lava dome. Here we use variations in the textures, glass compositions and volatile contents of melt inclusions preserved in pyroclasts produced by the summer 1980 eruptions to determine conditions of magma ascent and storage that may have led to observed changes in eruptive activity. Five different pyroclast types identified in pyroclastic flow and fall deposits produced by eruptions in June 12, July 22 and August 7, 1980, provide evidence for multiple levels of magma storage prior to each event. Highly vesicular clasts have H₂O-rich (4.5–5.5 wt%) melt inclusions and lack groundmass microlites or hornblende reaction rims, characteristics that require magma storage at P≥160 MPa until shortly prior to eruption. All other clast types have groundmass microlites; P_H20 estimated from both H₂O-bearing melt inclusions and textural constraints provided by decompression experiments suggest pre-eruptive storage pressures of ∼75, 40, and 10 MPa. The distribution of pyroclast types within and between eruptive deposits can be used to place important constraints on eruption mechanisms. Fall and flow deposits from June 12, 1980, lack highly vesicular, microlite-free pyroclasts. This eruption was also preceded by a shallow intrusion on June 3, as evidenced by a seismic crisis and enhanced SO₂ emissions. Our constraints suggest that magma intruded to a depth of ≤4 km beneath the crater floor fed the June eruption. In contrast, eruptions of July and August, although shorter in duration and smaller in volume, erupted deep volatile-rich magma. If modeled as a simple cylinder, these data require a step-wise decrease in effective conduit diameter from 40–50 m in May and June to 8–12 m in July and August. The abundance of vesicular (intermediate to deep) clast types in July and August further suggests that this change was effected by narrowing the shallower part of the conduit, perhaps in response to solidification of intruded magma remaining in the shallow system after the June eruption. Eruptions from July to October were distinctly pulsatory, transitioning between subplinian and vulcanian in character. As originally suggested by Scandone and Malone (1985), a growing mismatch between the rate of magma ascent and magma disruption explains the increasingly pulsatory nature of the eruptions through time. Recent fragmentation experiments Spieler et al. (2004) suggest this mismatch may have been aided by the multiple levels at which magma was stored (and degassed) prior to these events.Editorial responsibility: J Stix     

Analytical models for bubble growth during decompression of high viscosity magmas

Analytical models for decompressional bubble growth in a viscous magma are developed to establish the influence of high magma viscosity on vesiculation and to assess the time-scales on which bubbles respond to decompression. Instantaneous decompression of individual bubbles, analogous to a sudden release of pressure (e.g. sector collapse), is considered for two end-member cases. The infinite melt model considers the growth of an isolated bubble before significant bubble interaction occurs. The shell model considers the growth of a bubble surrounded by a thin shell and is analogous to bubble growth in a highly vesicular magmatic foam. Results from the shell model show that magmas less viscous than 10⁹ Pa s can freely expand without developing strong overpressures. The timescales for pressure re-equilibration are shortened by increased ratios of bubble radius to shell thickness and by larger decompression. Time-scales for isolated bubbles in rhyolitic melts (infinite melt model) are significantly longer, implying that such bubbles could experience internal pressures greater than the ambient pressure for at least a few hours following a sudden release of pressure. The shell model is developed to assess bubble growth during the linear decompression of a magma body of constant viscosity. For the range of decompression rates and viscosities associated with actual volcanic eruptions, bubble growth continues at approximately the equilibrium rate, with no attendant excess of internal pressure. The results imply that viscosity does not have any significant role in preventing the explosive expansion of high viscosity foams. However, for viscosities of >10⁹ Pa s there is the potential for a viscosity quench under the extreme decompression rates of an explosive eruption. It is proposed that the typical vesicularities of pumice of 0.7–0.8 are a consequence of the viscosity of the degassing magmas becoming sufficiently high to inhibit bubble expansion over the characteristic time-scale of eruption. For fully degassed silicic lavas with viscosities in the range 10¹⁰ to 10¹² Pa s time-scales for decompression of isolated bubbles can be hours to many months.     

Shallow seismicity beneath Ruapehu Crater Lake: results of a 1994 seismometer deployment

 Virtually all the seismicity within Ruapehu Volcano recorded during a 2-month deployment in early 1994, with 14 broadband seismographs around the Tongariro National Park volcanoes in the North Island of New Zealand, was associated with the active vent and occurred within approximately 1 km of Ruapehu Crater Lake. High-frequency volcano-tectonic earthquakes and low-frequency events (similar to bursts of 2 Hz volcanic tremor) were both found to have sources in this region. The high-frequency events, which often consisted of a smaller precursor event followed approximately 2 s later by the main event, had sharp onsets and were locatable using standard techniques. The depth of these events ranged from the surface down to approximately 1500 m below Crater Lake. The low-frequency events did not have sharp onsets and were located by phase-correlation methods. Nearly all occurred under a small region on the east side of Crater Lake, at depths from 200 to 1000 m below the surface. This low-frequency earthquake source region, in which no high-frequency events occurred, may be the steam zone within the actual vent of Ruapehu Volcano. Received: 30 June 1996 / Accepted: 16 February 1998     

Impacts of volcanism on pre-European inhabitants of Taveuni, Fiji

 Taveuni is a Fijian ocean-island volcano that sporadically erupted throughout the Holocene. The 437-km² island is an active monogenetic volcanic field with a constantly shifting locus of activity along a single apparent rift axis. Although the eruptions were not large ( ≤VEI 2), unexpected shifts in Taveuni volcanism had the potential to affect habitation sites. Since known human settlement of the Fiji Group (ca. 950–750 BC), there have been at least 58 eruptions on Taveuni. Up to 25 of these eruptions potentially affected pre-European inhabitants of the island and at least four former occupation sites are known to have been affected by volcanic products. Despite apparent earliest settlement of Taveuni post-dating other nearby islands by up to 600 years, volcanism probably did not hinder or stall settlement of Taveuni compared with neighbouring islands. However, a period of voluminous eruptions between 300 and 500 AD covered much of south Taveuni with lava and/or thick tephra, apparently causing abandonment of at least this portion of Taveuni until approximately 1100 AD. Most eruptions were not of catastrophic proportions and, due to their localised effects, re-settlement was rapid in marginal unaffected areas. Localised stories and a relict place name survive to describe former eruption locations and effects since approximately 120–320 AD. Knowledge of the impacts on Taveuni's past inhabitants forms the basis of volcanic disaster-mitigation strategies to minimise future effects on the current 14,500 residents. Received: 9 September 1999 / Accepted: 21 February 1999     

Role of carbon dioxide in the dynamics of magma ascent in explosive eruptions

 The role of carbon dioxide in the dynamics of magma ascent in explosive eruptions is investigated by means of numerical modeling. The model is steady, one-dimensional, and isothermal; it calculates the separated flow of gas and a homogeneous mixture of liquid magma and crystals. The magma properties are calculated on the basis of magma composition and crystal content and are allowed to change along the conduit due to pressure decrease and gas exsolution. The effect of the presence of a two-component (water + carbon dioxide) exsolving gas phase is investigated by performing a parametric study on the CO₂/(H₂O+CO₂) ratio, which is allowed to vary from 0 to 0.5 at either constant total volatile or constant water content. The relatively insoluble carbon dioxide component plays an important role in the location of the volatile-saturation and magma-fragmentation levels and in the distribution of the flow variables in the volcanic conduit. In detail, the results show that an increase of the proportion of carbon dioxide produces a decrease of the mass flow rate, pressure, and exit mixture density, and an increase of the exit gas volume fraction and depth of the fragmentation level. A relevant result is the different role played by water and carbon dioxide in the eruption dynamics; an increasing amount of water produces an increase of the mass flow rate, and an increasing amount of carbon dioxide produces a decrease. Even small amounts of carbon dioxide have major consequences on the eruption dynamics, implying that the multicomponent nature of the volcanic gas must be taken into account in the prediction of the eruption scenario and the forecasting of volcanic hazard. Received: 6 March 1998 / Accepted: 28 October 1998     

Statistical analyses and characteristics of volcanic tremor on Stromboli Volcano (Italy)

 A study of volcanic tremor on Stromboli is carried out on the basis of data recorded daily between 1993 and 1995 by a permanent seismic station (STR) located 1.8 km away from the active craters. We also consider the signal of a second station (TF1), which operated for a shorter time span. Changes in the spectral tremor characteristics can be related to modifications in volcanic activity, particularly to lava effusions and explosive sequences. Statistical analyses were carried out on a set of spectra calculated daily from seismic signals where explosion quakes were present or excluded. Principal component analysis and cluster analysis were applied to identify different classes of spectra. Three clusters of spectra are associated with two different states of volcanic activity. One cluster corresponds to a state of low to moderate activity, whereas the two other clusters are present during phases with a high magma column as inferred from the occurrence of lava fountains or effusions. We therefore conclude that variations in volcanic activity at Stromboli are usually linked to changes in the spectral characteristics of volcanic tremor. Site effects are evident when comparing the spectra calculated from signals synchronously recorded at STR and TF1. However, some major spectral peaks at both stations may reflect source properties. Statistical considerations and polarization analysis are in favor of a prevailing presence of P-waves in the tremor signal along with a position of the source northwest of the craters and at shallow depth. Received: 15 December 1996 / Accepted: 31 March 1998     

Premixing of magma and water in MFCI experiments

 Experimental studies have been performed to evaluate pre-explosive water–melt mixes with respect to explosive volcanic molten–fuel–coolant interaction (MFCI), i.e., phreatomagmatic explosion. Remolten ultrabasic volcanic rock was used as a magma simulant. Measurement of the explosion intensity was used to determine optimal premixing conditions. A well-defined optimal range was found for the hydrodynamic mixing energy (differential flow speed of 4.2 m/s), as well as for the water/melt mass ratio (0.03 to 0.04) under experimental conditions. The mass flux of water had a minor influence on the explosion intensity. Additionally, transparent mixing experiments with silicon oil and inked water were carried out. They indicate a direct dependence of the pre-explosive water-melt interface area on the explosion intensity. The experimental results show that the contact conditions of water and melt required for explosive MFCI may easily be established in natural volcanic systems. Thus, explosive MFCI is a probable mechanism of explosive volcanism. Received: 23 July 1996 / Accepted: 16 December 1996     

Victims from volcanic eruptions: a revised database

 The number of victims from volcanism and the primary cause(s) of death reported in the literature show considerable uncertainty. We present the results of investigations carried out either in contemporary accounts or in specific studies of eruptions that occurred since A.D. 1783. More than 220 000 people died because of volcanic activity during this period, which includes approximately 90% of the recorded deaths throughout history. Most of the fatalities resulted from post-eruption famine and epidemic disease (30.3%), nuées ardentes or pyroclastic flows and surges (26.8%), mudflows or lahars (17.1%), and volcanogenic tsunamis (16.9%). At present, however, international relief efforts might reduce the effects of post-eruption crop failure and disease, and at least some of the lahars could be anticipated in time by adequate scientific and social response. Thus, mitigation of hazards from pyroclastic flows and tsunamis will become of paramount importance to volcanologists and civil authorities. Received: 3 August 1997 / Accepted: 10 April 1998     

Permeability and degassing of dome lavas undergoing rapid decompression: An experimental determination

The gas permeability of volcanic rocks may influence various eruptive processes. The transition from a quiescent degassing dome to rock failure (fragmentation) may, for example, be controlled by the rocks permeability, in as much as it affects the speed by which a gas overpressure in vesicles is reduced in response to decompression. Using a modified shock-tube-based fragmentation bomb (Alidibirov and Dingwell 1996a,b; Spieler et al. 2003a), we have measured unsteady-state permeability at a high initial pressure differential. Following sudden decompression above the rock cylinder, pressurized gas flows through the sample. Two pressure transducers record the pressure signals above and below the sample. A transient 1D filtration code has been developed to calculate permeability using the experimental decay curve of the lower pressure transducer. Additionally an analytical steady-state method to achieve permeability is presented as an alternative to swiftly predict the sample permeability in a sufficiently precise manner. Over 100 permeability measurements have been performed on samples covering a wide range of porosity. The results show a general positive relationship between porosity and permeability with a high data scatter. Our preferred interpretation of the results is a combination of two different, but overlapping effects. We propose that at low porosities, gas escape occurs predominantly through microcracks or elongated micropores and therefore could be described by simplified forms of Kozeny–Carman relations (Carman 1956) and fracture flow models. At higher porosities, the influence of vesicles becomes progressively stronger as they form an increasingly connected network. Therefore, a model based on the percolation theory of fully penetrable spheres is used, as a first approximation, to describe the permeability-porosity trend. In the data acquired to date it is evident, that in addition to the porosity control, the samples bubble size, shape and distribution strongly influence the permeability. This leads to a range of permeability values up to 2.5 orders of magnitude at a given porosity.     

Stabilization of volcanic flanks by dike intrusion: an example from Kilauea

 Dike propagation and dilation increases the compression of adjacent rocks. On volcanoes, especially oceanic shields, dikes are accordingly thought to be structurally destabilizing. As compression is incremented, volcanic flanks are driven outward or downslope and thus increase their susceptibility to destructive earthquakes and giant landslides. We show, however, that the 2-m-thick dike emplaced along the east rift zone of Kilauea in 1983 actually stabilized that volcano's flank. Specifically, production of flank earthquakes dropped more than twofold after 1983 as maximum downslope motion slowed to 6 cm·year^–1 from approximately 40 cm·year^–1 during 1980–1982. As much as 65 cm of deflationary subsidence above Kilauea's summit and upper rift zones accompanied the dike intrusion. According to recent estimates, this deflation corresponds to a reduction in magma-reservoir pressure of approximately 4 MPa, probably about as much as the driving pressure of the 1983 dike. The volume of the dike, approximately 0.10–0.15 km³, is orders of magnitude less than the estimated 200- to 250-km³ volume of Kilauea's reservoir of magma and nearby hot, mushy rock. Thus, deflation of that reservoir reduces the compressional load on the flank over a much larger area than intrusion of the dike adds to it, particularly at the dominant depth of seismicity, 8–9 km. A Coulomb block model for flank motion during intervals between major earthquakes requires the low-angle fault beneath Kilauea's flank to exhibit slip weakening, conducive to earthquake instability. Accordingly, the triggering mechanism of destructive earthquakes, several of which have struck Hawaii during the past 150 years, need not require stresses accumulated by dike intrusions. Received: 27 October 1998 / Accepted: 24 May 1999