Using amphibole phenocrysts to track vapor transfer during magma crystallization and transport: An example from Mount St. Helens,Washington

In order to evaluate and further constrain models for volatile movement and vapor enrichment of magma stored at shallow levels, amphibole phenocrysts from 2004–2005 Mount St. Helens dacite were analyzed for major and selected trace elements (Li, Cu, Zn, Mn, and REE) and Li isotopes. Several recent studies have examined fluid-mobile trace element abundances in phencryst phases and melt inclusions as a means of tracking volatile movement within subvolcanic magmatic systems, and high Li contents in plagioclase phenocrysts from 1980 and 2004 Mount St. Helens dacites have been interpreted as evidence that shallow magma was fluxed by a Li-bearing vapor phase prior to eruption.     

Fumaroles in ice caves on the summit of Mount Rainier—preliminary stable isotope, gas, and geochemical studies

The edifice of Mount Rainier, an active stratovolcano, has episodically collapsed leading to major debris flows. The largest debris flows are related to argillically altered rock which leave areas of the edifice prone to failure. The argillic alteration results from the neutralization of acidic magmatic gases that condense in a meteoric water hydrothermal system fed by the melting of a thick mantle of glacial ice. Two craters atop a 2000-year-old cone on the summit of the volcano contain the world's largest volcanic ice-cave system. In the spring of 1997 two active fumaroles (T=62°C) in the caves were sampled for stable isotopic, gas, and geochemical studies.Stable isotope data on fumarole condensates show significant excess deuterium with calculated δD and δ¹⁸O values (−234 and −33.2‰, respectively) for the vapor that are consistent with an origin as secondary steam from a shallow water table which has been heated by underlying magmatic–hydrothermal steam. Between 1982 and 1997, δD of the fumarole vapor may have decreased by 30‰.The compositions of fumarole gases vary in time and space but typically consist of air components slightly modified by their solubilities in water and additions of CO₂ and CH₄. The elevated CO₂ contents (δ¹³C_CO2=−11.8±0.7‰), with spikes of over 10,000 ppm, require the episodic addition of magmatic components into the underlying hydrothermal system. Although only traces of H₂S were detected in the fumaroles, most notably in a sample which had an air δ¹³C_CO2 signature (−8.8‰), incrustations around a dormant vent containing small amounts of acid sulfate minerals (natroalunite, minamiite, and woodhouseite) indicate higher H₂S (or possibly SO₂) concentrations in past fumarolic gases.Condensate samples from fumaroles are very dilute, slightly acidic, and enriched in elements observed in the much higher temperature fumaroles at Mount St. Helens (K and Na up to the ppm level; metals such as Al, Pb, Zn Fe and Mn up to the ppb level and volatiles such as Cl, S, and F up to the ppb level).The data indicate that the hydrothermal system in the edifice at Mount Rainier consists of meteoric water reservoirs, which receive gas and steam from an underlying magmatic system. At present the magmatic system is largely flooded by the meteoric water system. However, magmatic components have episodically vented at the surface as witnessed by the mineralogy of incrustations around inactive vents and gas compositions in the active fumaroles. The composition of fumarole gases during magmatic degassing is distinct and, if sustained, could be lethal. The extent to which hydrothermal alteration is currently occurring at depth, and its possible influence on future edifice collapse, may be determined with the aid of on site analyses of fumarole gases and seismic monitoring in the ice caves.     

Crystallization of microlites during magma ascent: the fluid mechanics of 1980–1986 eruptions at Mount St Helens

Eruptions of Mount St Helens (Washington, USA) decreased in intensity and explosivity after the main May 18, 1980 eruption. As the post-May 18 eruptions progressed, albitic plagioclase microlites began to appear in the matrix glass, although the bulk composition of erupted products, the phenocryst compositions and magmatic temperatures remained fairly constant. Equilibrium experiments on a Mount St Helens white pumice show that at 160 MPa water pressure and 900°C, conditions deduced for the 8 km deep magma storage zone, the stable plagioclase is An₄₇. The microlites in the natural samples, which are more albitic, had to grow at lower water pressures during ascent. Isothermal decompression experiments reported here demonstrate that a decrease in water pressure from 160 to 2 MPa over four to eight days is capable of producing the albitic groundmass plagioclase and evolved melt compositions observed in post-May 18 1980 dacites. Because groundmass crystallization occurs over a period of days during and after decreases in pressure, microlite crystallization in the Mount St Helens dacites must have occurred during the ascent of each magma batch from a deep reservoir rather than continuously in a shallow holding chamber. This is consistent with data on the kinetics of amphibole breakdown, which require that a significant portion of magma vented in each eruption ascended from a depth of at least 6.5 km (160 MPa water pressure) in a matter of days. The size and shape of the microlite population have not been studied because of the small size of the experimental samples; it is possible that the texture continues to mature long after chemical equilibrium is approached. As the temperature, composition, crystal content and water content of magma in the deep reservoir remained approximately constant from May 1980 to at least March 1982, the spectacular decrease in eruption intensity during this period cannot be attributed to changes in viscosity or density of the magma. Simple fluld mechanical considerations indicate, however, that the observed changes in mass flux of magma can be modelled by a five-fold decrease in conduit radius from 35 to 7 m, produced perhaps by plating of magma along the conduit walls. The decreased ascent rates which accompanied the decrease in conduit radius can explain the change from closed-system to open-system degassing and the shift from explosive to effusive eruptions during 1980.     

Water contents, temperatures and diversity of the magmas of the catastrophic eruption of Nevado del Ruiz, Colombia, November 13, 1985

The petrology of the highly phyric two-pyroxene andesitic to dacitic pyroclastic rocks of the November 13, 1985 eruption of Nevado del Ruiz, Colombia, reveals evidence of: (1) increasingly fractionated bulk compositions with time; (2) tapping of a small magma chamber marginally zoned in regard to H₂O contents (1 to 4%), temperature (960–1090°C), and amount of residual melt (35 to 65%); (3) partial melting and assimilation of degassed zones in the hotter less dense interior of the magma chamber; (4) probable heating, thermal disruption and mineralogic and compositional contamination of the magma body by basaltic magma “underplating”; and (5) crustal contamination of the magmas during ascent and within the magma chamber. Near-crater fall-back or “spill-over” emitted in the middle of the eruptive sequence produced a small pyroclastic flow that became welded in its central and basal portions because of ponding and thus heat conservation on the flat glaciated summit near the Arenas crater. The heterogeneity of Ruiz magmas may be related to the comparatively small volume (0.03 km³) of the eruption, nearly ten times less than the 0.2 km³ of the Plinian phase of Mount St. Helens, and probable steep thermal and P_H2O gradients of a small source magma chamber, estimated at 300 m long and 100 m wide for an assumed ellipsoidal shape.     

Modelling the dynamics and thermodynamics of volcanic degassing

 The rates of passive degassing from volcanoes are investigated by modelling the convective overturn of dense degassed and less dense gas-rich magmas in a vertical conduit linking a shallow degassing zone with a deep magma chamber. Laboratory experiments are used to constrain our theoretical model of the overturn rate and to elaborate on the model of this process presented by Kazahaya et al. (1994). We also introduce the effects of a CO₂–saturated deep chamber and adiabatic cooling of ascending magma. We find that overturn occurs by concentric flow of the magmas along the conduit, although the details of the flow depend on the magmas' viscosity ratio. Where convective overturn limits the supply of gas-rich magma, then the gas emission rate is proportional to the flow rate of the overturning magmas (proportional to the density difference driving convection, the conduit radius to the fourth power, and inversely proportional to the degassed magma viscosity) and the mass fraction of water that is degassed. Efficient degassing enhances the density difference but increases the magma viscosity, and this dampens convection. Two degassing volcanoes were modelled. At Stromboli, assuming a 2 km deep, 30% crystalline basaltic chamber, containing 0.5 wt.% dissolved water, the ∼700 kg s^–1 magmatic water flux can be modelled with a 4–10 m radius conduit, degassing 20–100% of the available water and all of the 1 to 4 vol.% CO₂ chamber gas. At Mount St. Helens in June 1980, assuming a 7 km deep, 39% crystalline dacitic chamber, containing 4.6 wt.% dissolved water, the ∼500 kg s^–1 magmatic water flux can be modelled with a 22–60 m radius conduit, degassing ∼2–90% of the available water and all of the 0.1 to 3 vol.% CO₂ chamber gas. The range of these results is consistent with previous models and observations. Convection driven by degassing provides a plausible mechanism for transferring volatiles from deep magma chambers to the atmosphere, and it can explain the gas fluxes measured at many persistently active volcanoes. Received: 26 September 1997 / Accepted: 11 July 1998     

Pre-eruptive and syn-eruptive conditions in the Black Butte, California dacite: Insight into crystallization kinetics in a silicic magma system

A series of experiments and petrographic analyses have been run to determine the pre-eruption phase equilibria and ascent dynamics of dacitic lavas composing Black Butte, a dome complex on the flank of Mount Shasta, California. Major and trace element analyses indicate that the Black Butte magma shared a common parent with contemporaneously erupted magmas at the Shasta summit. The Black Butte lava phenocryst phase assemblage (20 v.%) consists of amphibole, plagioclase (core An_77.5), and Fe–Ti oxides in a fine-grained (< 0.5 mm) groundmass of plagioclase, pyroxene, Fe–Ti oxides, amphibole, and cristobalite. The phenocryst assemblage and crystal compositions are reproduced experimentally between 890 °C and 910 °C, ≥ 300 MPa, X_H2O = 1, and oxygen fugacity = NNO + 1. This study has quantified the extent of three crystallization processes occurring in the Black Butte dacite that can be used to discern ascent processes. Magma ascent rate was quantified using the widths of amphibole breakdown rims in natural samples, using an experimental calibration of rim development in a similar magma at relevant conditions. The majority of rims are 34 ± 10 μm thick, suggesting a time-integrated magma ascent rate of 0.004–0.006 m/s among all four dome lobes. This is comparable to values for effusive samples from the 1980 Mount St. Helens eruption and slightly faster than those estimated at Montserrat. A gap between the compositions of plagioclase phenocryst cores and microlites suggests that while phenocryst growth was continuous throughout ascent, microlite formation did not occur until significantly into ascent. The duration of crystallization is estimated using the magma reservoir depth and ascent rate, as determined from phase equilibria and amphibole rim widths, respectively. Textural analysis of the natural plagioclase crystals yields maximum growth rates of plagioclase phenocryst rims and groundmass microlites of 8.7 × 10^− 8 and 2.5 × 10^− 8 mm/s, respectively. These rates are comparable to values determined from time-sequenced samples of dacite erupted effusively from Mount St. Helens during 1980 and indicate that syneruptive crystallization processes were important during the Black Butte eruptive cycle.     

Tritium and stable isotopes of magmatic waters

To investigate the isotopic composition and age of water in volcanic gases and magmas, we analyzed samples from 11 active volcanoes ranging in composition from tholeiitic basalt to rhyolite: Mount St. Helens (USA), Kilauea (USA), Pacaya (Guatemala), Galeras (Colombia), Satsuma Iwo-Jima (Japan), Sierra Negra and Alcedo (Ecuador), Vulcano (Italy), Parícutin (Mexico), Kudryavy (Russia), and White Island (New Zealand). Tritium at relatively low levels (0.1–5 T.U.) is found in most emissions from high-temperature volcanic fumaroles sampled, even at discharge temperatures >700°C. Although magmatic fluids sampled from these emissions usually contain high CO₂, S_total, HCl, HF, B, Br, ³He R/R_A, and low contents of air components, stable isotope and tritium relations of nearly all such fluids show mixing of magmatic volatiles with relatively young meteoric water (model ages≤75 y). Linear δD/δ¹⁸O and ³H/δ¹⁸O mixing trends of these two end-members are invariably detected at arc volcanoes. Tritium is also detected in fumarole condensates at hot spot basalt volcanoes, but collecting samples approaching the composition of end-member magmatic fluid is exceedingly difficult. In situ production of ³H, mostly from spontaneous fission of ²³⁸U in magmas is calculated to be <0.001 T.U., except for the most evolved compositions (high U, Th, and Li and low H₂O contents). These values are below the detection limit of ³H by conventional analytical techniques (about 0.01 T.U. at best). We found no conclusive evidence that natural fusion in the Earth produces anomalous amounts of detectable ³H (>0.05 T.U.).     

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     

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     

Mount St. Helens a decade after the 1980 eruptions: magmatic models,chemical cycles,and a revised hazards assessment

Available geophysical and geologic data provide a simplified model of the current magmatic plumbing system of Mount St. Helens (MSH). This model and new geochemical data are the basis for the revised hazards assessment presented here. The assessment is weighted by the style of eruptions and the chemistry of magmas erupted during the past 500 years, the interval for which the most detailed stratigraphic and geochemical data are available. This interval includes the Kalama (A. D. 1480–1770s?), Goat Rocks (A.D. 1800–1857), and current eruptive periods. In each of these periods, silica content decreased, then increased. The Kalama is a large amplitude chemical cycle (SiO₂: 57%–67%), produced by mixing of arc dacite, which is depleted in high field-strength and incompatible elements, with enriched (OIB-like) basalt. The Goat Rocks and current cycles are of small amplitude (SiO₂: 61%–64% and 62%–65%) and are related to the fluid dynamics of magma withdrawal from a zoned reservoir. The cyclic behavior is used to forecast future activity. The 1980–1986 chemical cycle, and consequently the current eruptive period, appears to be virtually complete. This inference is supported by the progressively decreasing volumes and volatile contents of magma erupted since 1980, both changes that suggest a decreasing potential for a major explosive eruption in the near future. However, recent changes in seismicity and a series of small gas-release explosions (beginning in late 1989 and accompanied by eruption of a minor fraction of relatively low-silica tephra on 6 January and 5 November 1990) suggest that the current eruptive period may continue to produce small explosions and that a small amount of magma may still be present within the conduit. The gas-release explosions occur without warning and pose a continuing hazard, especially in the crater area. An eruption as large or larger than that of 18 May 1980 (0.5 km³ dense-rock equivalent) probably will occur only if magma rises from an inferred deep (7 km), relative large (5–7 km³) reservoir. A conservative approach to hazard assessment is to assume that this deep magma is rich in volatiles and capable of erupting explosively to produce voluminous fall deposits and pyroclastic flows. Warning of such an eruption is expectable, however, because magma ascent would probably be accompanied by shallow seismicity that could be detected by the existing seismic-monitoring system. A future large-volume eruption (0.1 km³) is virtually certain; the eruptive history of the past 500 years indicates the probability of a large explosive eruption is at least 1% annually. Intervals between large eruptions at Mount St. Helens have varied widely; consequently, we cannot confidently forecast whether the next large eruption will be years decades, or farther in the future. However, we can forecast the types of hazards, and the areas that will be most affected by future large-volume eruptions, as well as hazards associated with the approaching end of the current eruptive period.     

Eruptive activity at Mount St Helens,Washington, USA, 1984–1988: a gas geochemistry perspective

The results from two different types of gas measurement, telemetered in situ monitoring of reducing gases on the dome and airborne measurements of sulfur dioxide emission rates in the plume by correlation spectrometry, suggest that the combination of these two methods is particularly effective in detecting periods of enhanced degassing that intermittently punctuate the normal background leakage of gaseous effluent from Mount St Helens to the atmosphere. Gas events were recorded before lava extrusion for each of the four dome-building episodes at Mount St Helens since mid-1984. For two of the episodes, precursory reducing gas peaks were detected, whereas during three of the episodes, COSPEC measurements recorded precursory degassing of sulfur dioxide. During one episode (October 1986), both reducing gas monitoring and SO₂ emission rate measurements simultaneously detected a large gas release several hours before lava extrusion. Had both types of gas measurements been operational during each of the dome-building episodes, it is thought that both would have recorded precursory signals for all four episodes. Evidence from the data presented herein suggests that increased degassing at Mount St Helens becomes detectable when fresh upward-moving magma is between 2 km and a few hundred meters below the base of the dome and between about 60 and 12 hours before the surface extrusion of lava.     

Forecasts and predictions of eruptive activity at Mount St. Helens, USA: 1975–1984

Public statements about volcanic activity at Mount St. Helens include factual statements, forecasts, and predictions. A factual statement describes current conditions but does not anticipate future events. A forecast is a comparatively imprecise statement of the time, place, and nature of expected activity. A prediction is a comparatively precise statement of the time, place, and ideally, the nature and size of impending activity. A prediction usually covers a shorter time period than a forecast and is generally based dominantly on interpretations and measurements of ongoing processes and secondarily on a projection of past history. The three types of statements grade from one to another, and distinctions are sometimes arbitrary.Forecasts and predictions at Mount St. Helens became increasingly precise from 1975 to 1982. Stratigraphic studies led to a long-range forecast in 1975 of renewed eruptive activity at Mount St. Helens, possibly before the end of the century. On the basis of seismic, geodetic and geologic data, general forecasts for a landslide and eruption were issued in April 1980, before the catastrophic blast and landslide on 18 May 1980. All extrusions except two from June 1980 to the end of 1984 were predicted on the basis of integrated geophysical, geochemical, and geologic monitoring. The two extrusions that were not predicted were preceded by explosions that removed a substantial part of the dome, reducing confining pressure and essentially short-circuiting the normal precursors.     

An analysis of the earthquake time series at Mount St. Helens, Washington, in the framework of recent volcanic activity

The analysis of the earthquake time distribution at Mount St. Helens reveals a good correlation between the physical state of the volcano and statistical parameters of earthquake sequence. There are three main seismic phases in the whole 1980–1986 period. The first one precedes the main eruption of May 18, 1980. It begins with a sudden increase of the seismicity level in late March and continues with an Utsu (1961) type decay of the seismic occurrence rate, characterized by a small value of the decay coefficient, β. The second phase lasts from the cataclysmic eruption on May 18, 1980 until the continuous dome building episode in 1983 and is characterized by a very slow exponential increase of the background level of seismicity. The third phase covers the remaining part of the sample and is characterized by a stationary earthquake clustering process episodically interrupted by peaks of activity related to eruptions. The trends in seismic occurrence rate within each phase, as well as the statistical parameter variations at each transition, are analyzed and discussed in the framework of volcanic activity. This leads to the conclusion that statistical techniques may give a significant contribution in understanding changes in volcanic processes such as those at Mount St. Helens.     

Radionuclides and sulfur content in Mount Etna plume in 1983–1995: new constraints on the magma feeding system

Radionuclide activities (²¹⁰Pb, ²¹⁰Bi, ²¹⁰Po) were investigated in Mount Etna plume from 1983 to 1995. At SE crater the long-term observation (12 years) of the ²¹⁰Po/²¹⁰Pb ratio shows that it can behave as a degassing vent not directly related to the main magma reservoir depending on the magma level inside the volcano. Since 1992, the simultaneous determination of radionuclides and sulfur in the main plume results in new constraints on the degassing model of Lambert et al. (Earth Planet. Sci. Lett., 76 (1986) 185). The ²¹⁰Po/SO₂ and ²¹⁰Pb/SO₂ ratios enable us to identify two sources of ²¹⁰Po in the plume: one is magmatic, correlated with SO₂, the other one is an additional component issued from the decay of ²¹⁰Pb in the shallow degassing cell, and depends directly on the residence time of the gases before their emission. Estimations of the volume of degassing magma, the residence time of the gases and the proportion of undegassed magma renewing the shallow degassing cell are given for the period 1992–1995. During the 1992 eruption, the rate of degassing magma volume is estimated to have been as high as 5×10⁶ m³/day, and the volume of the shallow degassing magma reservoir about 0.5 km³. In 1994 and 1995 the rate of non-erupted degassing magma volume was estimated to have been about 0.18 km³/year. During the entire 1983–1995 period, only 15–20% of the degassed magma was erupted.     

Earthquake-induced variations in the composition of the water in the geothermal reservoir at Vulcano Island, Italy

During 1979–1989, variations were observed in the oxygen composition of the water contained in the geothermal reservoir at Vulcano Island, Italy.The reservoir water, that has a magmatic origin, showed an oxygen composition of +1.0±0.5‰ δ¹⁸O during periods without local tectonic earthquakes, and an oxygen composition of +3.4±0.5‰ δ¹⁸O after the highest-energy seismic activity that occurred recently near the island. A slight increase of the δ¹⁸O value in the reservoir water was also observed after a low-energy sequence of tectonic earthquakes that occurred at very shallow depth just beneath Vulcano Island. These ¹⁸O variations in the reservoir water are consistent with earthquake-induced increases in the contribution from high-temperature δ¹⁸O-rich magmatic condensate to the geothermal reservoir, and with subsequent decreases in the δ¹⁸O value due to ¹⁸O exchanges at the temporarily increased reservoir temperature during reactions between the highly reactive magmatic condensate and the local rocks.Only minor changes in the deuterium composition of the reservoir water occurred with time, as the δD value in the magmatic condensate released from the magma after major local earthquakes quickly approached the δD value of the water contained in the geothermal reservoir.Also the chloride concentration in the reservoir water appears to be linked to the contribution from the magmatic fluid. This chloride content seems not to have undergone major changes with time, as it may be buffered by temporary increases in the reservoir temperature up to values >300°C induced by major local earthquakes. This mechanism may possibly occur also in other magmatic–hydrothermal systems.     

Crystallization of Mount St. Helens 1980–1986 dacite: A quantitative textural approach

Quantitative measurements of crystal size distributions (CSDs) have been used to obtain kinetic information on crystallization of industrial compounds (Randolph and Larson 1971) and more recently on Hawaiian basalts (Cashman and Marsh 1988). The technique is based on a population balance resulting in a differential equation relating the population densityn of crystals to crystal sizeL, i.e., at steady staten =n ^o exp(–L/itG), wheren ^o is nucleation density,G is the average crystal growth rate, is the average growth time, and the nucleation rateJ =n ^o G. CSD (Inn vsL) plots of plagioclase phenocrysts in 12 samples of Mount St. Helens blast dacite and 14 samples of dacite from the 1980–1986 Mount St. Helens dome are similar and averageG = 9.6 (± 1.1) × 10^–3 cm andn ^o = 1–2 × 10⁶ cm^–4. Reproducibility of the measurements was tested by measuring CSDs of 12 sections cut from a single sample in three mutually perpendicular directions; precision of the size distributions is good in terms of relative, but not necessarily absolute values (± 10%). Growth and nucleation rates for plagioclase have been calculated from these measurements using time brackets of = 30–150 years; growth ratesG are 3–10 × 10^–12cm/s, and nucleation ratesJ are 5–21 × 10^–6/cm³ s.G andJ for Fe-Ti oxides calculated from CSD data areG = 2–13 ± 10^–13 cm/sec andJ = 7–33 × 10^–5/cm³ s, respectively. The higher nucleation rate and lower growth rate of oxides resulted in a smaller average crystal size than for plagioclase. Sizes of plagioclase microlites (<0.01 mm) in the blast dacite groundmass have been measured from backscatter SEM photographs. Nucleation of these microlites was probably triggered by intrusion of material into the cone of Mount St. Helens in spring 1980. This residence time of 52 days gives minimum crystallization estimates ofG = 1–3 × 10^–11 cm/s andJ = 9–16 × 1O³/cm³ s. The skeletal form of the microlites provides evidence for nucleation and growth at high values of undercooling (T) relative to the phenocrysts. A comparison of nucleation and growth rates for the two crystal populations (phenocrysts vs microlites) suggests that while growth rate seems to be only slightly affected by changes inT, nucleation rate is a very strong function of undercooling. A comparison of plagioclase nucleation and growth rates measured in the Mount St. Helens dacite and in basalt from Makaopuhi lava lake in Hawaii suggests that plagioclase nucleation rates are not as dependent on composition. Groundmass textures suggest that plagioclase microphenocrysts crystallized at depth rather than in the conduit, in the dome, or after extrusion onto the surface. Most of this crystallization appears to be in the form of crystal growth (coarsening) of groundmass microphenocrysts at small degrees of undercooling rather than extensive nucleation of new crystals. This continuous crystallization in a shallow magmatic reservoir may provide the overpressurization needed to drive the continuing periodic domebuilding extrusions, which have been the pattern of activity at Mount St. Helens since December 1980.     

Evaluation of gases,condensates, and SO2 emissions from Augustine volcano,Alaska: the degassing of a Cl-rich volcanic system

After the March–April 1986 explosive eruption a comprehensive gas study at Augustine was undertaken in the summers of 1986 and 1987. Airborne COSPEC measurements indicate that passive SO₂ emission rates declined exponentially during this period from 380±45 metric tons/day (T/D) on 7/24/86 to 27±6 T/D on 8/24/87. These data are consistent with the hypothesis that the Augustine magma reservoir has become more degassed as volcanic activity decreased after the spring 1986 eruption. Gas samples collected in 1987 from an 870°C fumarole on the andesitic lava dome show various degrees of disequilibrium due to oxidation of reduced gas species and condensation (and loss) of H₂O in the intake tube of the sampling apparatus. Thermochemical restoration of the data permits removal of these effects to infer an equilibrium composition of the gases. Although not conclusive, this restoration is consistent with the idea that the gases were in equilibrium at 870°C with an oxygen fugacity near the Ni–NiO buffer. These restored gas compositions show that, relative to other convergent plate volcanoes, the Augustine gases are very HCl rich (5.3–6.0 mol% HCl), S rich (7.1 mol% total S), and H₂O poor (83.9–84.8 mol% H₂O). Values of D and ¹⁸O suggest that the H₂O in the dome gases is a mixture of primary magmatic water (PMW) and local seawater. Part of the Cl in the Augustine volcanic gases probably comes from this shallow seawater source. Additional Cl may come from subducted oceanic crust because data by Johnston (1978) show that Cl-rich glass inclusions in olivine crystals contain hornblende, which is evidence for a deep source (>25km) for part of the Cl. Gas samples collected in 1986 from 390°–642°C fumaroles on a ramp surrounding the inner summit crater have been oxidized so severely that restoration to an equilibrium composition is not possible. H and O isotope data suggest that these gases are variable mixtures of seawater, FMW, and meteoric steam. These samples are much more H₂O-rich (92%–97% H₂O) than the dome gases, possibly due to a larger meteoric steam component. The 1986 samples also have higher Cl/S, S/C, and F/Cl ratios, which imply that the magmatic component in these gases is from the more degassed 1976 magma. Thus, the 1987 samples from the lava dome are better indicators than the 1986 samples of degassing within the Augustine magma reservoir, even though they were collected a year later and contain a significant seawater component. Future gas studies at Augustine should emphasize fumaroles on active lava domes. Condensates collected from the same lava-dome fumarole have enrichments ot 10⁷–10² in Cl, Br, F, B, Cd, As, S, Bi, Pb, Sb, Mo, Zn, Cu, K, Li, Na, Si, and Ni. Lower-temperature (200°–650°C) fumaroles around the volcano are generally less enriched in highly volatile elements. However, these lower-termperature fumaroles have higher concentration of rock-forming elements, probably derived from the wall rock.     

Some regularity in the process of re-awakening of andesite and dacite volcanoes: specific features of the 1982 El Chichón volcano,México reactivation

El Chichón volcano is an andesite stratovolcano in southern México. It erupted in March 1982, after about 550 years of quiescence. The 1982 eruption of El Chichón has not been followed by the growth of a lava dome within the newly formed crater. This is rather anomalous since the construction of a new dome after the destruction of an old one is a common process during the eruptions at andesite and dacite volcanoes. To discuss this anomalous aspect of the El Chichón eruption, some regularity in the process of re-awakening of dormant (here defined as a period of quiescence of more than 100 years) andesite and dacite volcanoes are studied based on the seismic activity recorded at the volcanoes Bezymianny, Mount St. Helens, El Chichón, Unzen, Pinatubo and Soufrière Hills. Three stages were identified in the re-awakening activity of these volcanoes: (1) preliminary seismic activity, leading up to the first phreatic explosion; (2) activity between the first and the largest explosions; (3) post-explosion dome-building process. The eruptions were divided into two groups: low-VEI (Volcanic Explosivity Index) and the long duration stage-1 events (Unzen, 1991 and Soufrière Hills volcano, 1995) and high-VEI and the short duration stage-1 events (Bezymianny, 1956; Mount St. Helens, 1980; El Chichón, 1982 and Pinatubo, 1992). The comparative analysis of the seismo-eruptive activity of two eruptions of the second group, the 1980 of Mt. St. Helens and the 1982 of El Chichón, produced an explanation the absence of new dome building during the 1982 eruption of El Chichón volcano. It may be explained in terms of the unusually rapid emission of gas and water from the magmatic and hydrothermal system beneath the volcano during a relatively short sequence of large explosions that could have sharply increased the viscosity of the magma making impossible its exit to the surface.     

Sulfur and chlorine degassing from primitive arc magmas: temporal changes during the 1982–1983 eruptions of Galunggung (West Java, Indonesia)

The 1982–1983 eruptions of Galunggung represent a nine-month period of intermittent volcanic activity with significant changes in explosivity and emission of volatiles. Eruptions started with Vulcanian explosions but changed gradually to Strombolian activity. Compositions of juvenile material changed from basaltic andesite to high-Mg basalt, which are among the most primitive rock types known in the Indonesian arc system. Although bulk compositions suggest a single evolution trend, we infer from the compositions of melt inclusions in olivine phenocrysts that the magmas represent derivatives of a complex spectrum of primary melts. Primitive inclusions in olivine phenocrysts from magma erupted during the Strombolian phase contain up to 2000 ppm sulfur, but concentrations decrease rapidly with increasing SiO₂ down to matrix glass values (50–100 ppm). ‘Vulcanian’ inclusions appear to be degassed before eruption (200 ppm S). Chlorine concentrations increase from 750 to 2200 ppm in Strombolian, and from 800 to 1500 in Vulcanian magmas, whereas matrix glass contains about 1000 ppm in both cases. Ash leachates show two cycles of decreasing S/Cl ratios: from 9.7 to 5.6 at the start of the activity, and from 12.2 to 2.0 after four months. As the second cycle follows upon increased seismic activity at shallow depth, it probably reflects degassing of fresh sulfur-rich magma arriving in the shallow Galunggung reservoir. In contrast to the degassed state of Vulcanian magma, the significant amounts of adsorbed sulfur on the ashes point to an excess source of sulfur, which was most likely derived from intruding Strombolian magma. Hence, the observed sulfur flux of 2 Mt is not in accordance with a petrologic estimate of 0.09 Mt. Using a published value of 550 Mt of erupted material about 0.34 km³ fresh undegassed magma is needed to account for the observed sulfur flux. This is close to the erupted volume of Vulcanian magma (0.26 km³), which presumably was replaced completely by Strombolian magma during the eruption. Using the petrologic method, we calculate a total release of 0.3 Mt chlorine, which agrees well with an output of 0.47 Mt estimated independently from S/Cl ratios of the ash leachates and TOMS sulfur yields. Ash leachates show that about 35% of the sulfur and 30% of the chlorine was scavenged from the eruption plumes. Our results suggest that sulfur and chlorine were largely decoupled during degassing, which resulted in considerable variations in S/Cl ratios during the Galunggung eruptions. We infer that sulfur degassing reflects the arrival of fresh magma at shallow depth, whereas chlorine is largely derived from simultaneously erupted material. As a consequence, the petrologic estimates are more consistent with observed emissions for chlorine than for sulfur.     

Physical properties of the 1980 Mount St. Helens cryptodome magma

 Physical properties of cryptodome and remelted samples of the Mount St. Helens grey dacite have been measured in the laboratory. The viscosity of cryptodome dacite measured by parallel–plate viscometry ranges from 10.82 to 9.94 log₁₀ η (Pa s) (T=900–982  °C), and shrinkage effects were dilatometrically observed at T>900  °C. The viscosity of remelted dacite samples measured by the micropenetration method is 10.60–9.25 log₁₀ η (Pa s) (T=736–802  °C) and viscosities measured by rotational viscometry are 3.22–1.66 log₁₀ η (Pa s) (T=1298–1594  °C). Comparison of the measured viscosity of cryptodome dacitic samples with the calculated viscosity of corresponding water-bearing melt demonstrates significant deviations between measured and calculated values. This difference reflects a combination of the effect of crystals and vesicles on the viscosity of dacite as well as the insufficient experimental basis for the calculation of crystal-bearing vesicular melt viscosities at low temperature. Assuming that the cryptodome magma of the 18 May 1980 Mount St. Helens eruption was residing at 900  °C with a phenocryst content of 30 vol.%, a vesicularity of 36 vol.% and a bulk water content of 0.6 wt.%, we estimate the magma viscosity to be 10^10.8 Pa s. Received: 25 August 1996 / Accepted: 19 July 1997