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     

Eruptive history and petrology of Mount Drum volcano,Wrangell Mountains,Alaska

Mount Drum is one of the youngest volcanoes in the subduction-related Wrangell volcanic field (80×200 km) of southcentral Alaska. It lies at the northwest end of a series of large, andesite-dominated shield volcanoes that show a northwesterly progression of age from 26 Ma near the Alaska-Yukon border to about 0.2 Ma at Mount Drum. The volcano was constructed between 750 and 250 ka during at least two cycles of cone building and ring-dome emplacement and was partially destroyed by violent explosive activity probably after 250 ka. Cone lavas range from basaltic andesite to dacite in composition; ring-domes are dacite to rhyolite. The last constructional activity occurred in the vicinity of Snider Peak, on the south flank of the volcano, where extensive dacite flows and a dacite dome erupted at about 250 ka. The climactic explosive eruption, that destroyed the top and a part of the south flank of the volcano, produced more than 7 km³ of proximal hot and cold avalanche deposits and distal mudflows. The Mount Drum rocks have medium-K, calc-alkaline affinities and are generally plagioclase phyric. Silica contents range from 55.8 to 74.0 wt%, with a compositional gap between 66.8 and 72.8 wt%. All the rocks are enriched in alkali elements and depleted in Ta relative to the LREE, typical of volcanic arc rocks, but have higher MgO contents at a given SiO₂, than typical orogenic medium-K andesites. Strontium-isotope ratios vary from 0.70292 to 0.70353. The compositional range of Mount Drum lavas is best explained by a combination of diverse parental magmas, magma mixing, and fractionation. The small, but significant, range in ⁸⁷Sr/⁸⁶Sr ratios in the basaltic andesites and the wide range of incompatible-element ratios exhibited by the basaltic andesites and andesites suggests the presence of compositionally diverse parent magmas. The lavas show abundant petrographic evidence of magma mixing, such as bimodal phenocryst size, resorbed phenocrysts, reaction rims, and disequilibrium mineral assemblages. In addition, some dacites and andesites contain Mg and Ni-rich olivines and/or have high MgO, Cr, Ni, Co, and Sc contents that are not in equilibrium with the host rock and indicate mixing between basalt or cumulate material and more evolved magmas. Incompatible element variations suggest that fractionation is responsible for some of the compositional range between basaltic andesite and dacite, but the rhyolites have K, Ba, Th, and Rb contents that are too low for the magmas to be generated by fractionation of the intermediate rocks. Limited Sr-isotope data support the possibility that the rhyolites may be partial melts of underlying volcanic rocks. Received March 13, 1993/Accepted September 10, 1993     

Petrology of gabbroic xenoliths in 1960 Kilauea basalt: crystalline remnants of prior (1955) magmatism

The 1960 Kapoho lavas of Kilauea’s east rift zone contain 1–10 cm xenoliths of olivine gabbro, olivine gabbro-norite, and gabbro norite. Textures are poikilitic (ol+sp+cpx in pl) and intergranular (cpx+pl±ol±opx). Poikilitic xenoliths, which we interpret as cumulates, have the most primitive mineral compositions, Fo_82.5, cpx Mg# 86.5, and An_80.5. Many granular xenoliths (ol and noritic gabbro) contain abundant vesicular glass that gives them intersertal, hyaloophitic, and overall ‘open’ textures to suggest that they represent ‘mush’ and ‘crust’ of a magma crystallization environment. Their phase compositions are more evolved (Fo_80–70, cpx Mg# 82–75, and An_73–63) than those of the poikilitic xenoliths. Associated glass is basaltic, but evolved (MgO 5 wt%; TiO₂ 3.7–5.8 wt%). The gabbroic xenolith mineral compositions fit existing fractional crystallization models that relate the origins of various Kilauea lavas to one another. FeO/MgO crystal–liquid partitioning is consistent with the poikilitic ol-gabbro assemblage forming as a crystallization product from Kilauea summit magma with ∼8 wt% MgO that was parental to evolved lavas on the east rift zone. For example, least squares calculations link summit magmas to early 1955 rift-zone lavas (∼5 wt% MgO) through ∼28–34% crystallization of the ol+sp+cpx+pl that comprise the poikilitic ol-gabbros. The other ol-gabbro assemblages and the olivine gabbro-norite assemblages crystallized from evolved liquids, such as represented by the early 1955 and late 1955 lavas (∼6.5 wt% MgO) of the east rift zone. The eruption of 1960 Kapoho magmas, then, scoured the rift-zone reservoir system to entrain portions of cumulate and solidification zones that had coated reservoir margins during crystallization of prior east rift-zone magmas. Received: January 7, 1993/Accepted: November 23, 1993     

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     

Effects of eruption and lava drainback on the H2O contents of basaltic magmas at Kilauea Volcano

 Lava drainback has been observed during many eruptions at Kilauea Volcano: magma erupts, degasses in lava fountains, collects in surface ponds, and then drains back beneath the surface. Time series data for melt inclusions from the 1959 Kilauea Iki picrite provide important evidence concerning the effects of drainback on the H₂O contents of basaltic magmas at Kilauea. Melt inclusions in olivine from the first eruptive episode, before any drainback occurred, have an average H₂O content of 0.7±0.2 wt.%. In contrast, many inclusions from the later episodes, erupted after substantial amounts of surface degassed lava had drained back down the vent, have H₂O contents that are much lower (≥0.24 wt.% H₂O). Water contents in melt inclusions from magmas erupted at Pu'u 'O'o on the east rift zone vary from 0.39–0.51 wt.% H₂O in tephra from high fountains to 0.10–0.28 wt.% H₂O in spatter from low fountains. The low H₂O contents of many melt inclusions from Pu'u 'O'o and post-drainback episodes of Kilauea Iki reveal that prior to crystallization of the enclosing olivine host, the melts must have exsolved H₂O at pressures substantially less than those in Kilauea's summit magma reservoir. Such low-pressure H₂O exsolution probably occurred as surface degassed magma was recycled by drainback and mixing with less degassed magma at depth. Recognition of the effects of low-pressure degassing and drainback leads to an estimate of 0.7 wt.% H₂O for differentiated tholeiitic magma in Kilauea's summit magma storage reservoir. Data for MgO-rich submarine glasses (Clague et al. 1995) and melt inclusions from Kilauea Iki demonstrate that primary Kilauean tholeiitic magma has an H₂O/K₂O mass ratio of ∼1.3. At transition zone and upper mantle depths in the Hawaiian plume source, H₂O probably resides partly in a small amount of hydrous silicate melt. Received: 31 March 1997 / Accepted: 17 November 1997     

Temporal variations in gravity at Mt. Etna (Italy) associated with the 1989 and 1991 eruptions

 Results are presented from 11 microgravity surveys on Mt. Etna between 1987 and 1993, a period including the major 1989 and 1991–1993 flank eruptions and subordinate 1990 activity. Measurements were made with LaCoste and Romberg D-62 and D-157 gravity meters along a network around the volcano between 1000 and 1900 m a.s.l. and, since 1992, a N–S summit profile. Gravity changes of as much as 200 μGal were observed at scales from the size of the summit region to that of the volcano. None was associated with significant changes in ground elevation. The data show an increase in gravity for 2 years before the 1989 eruption. The increase is attributed to the accumulation of magma (0.25–1.7×10⁹m³) in an elongate zone, oriented NNW–SSE, between 2.5 and 6 km below sea level. Part of this magma was injected into the volcanic pile to supply the 1989 and 1990 eruptions. It also probably fed the start of the 1991–1993 eruption, since this event was not preceded by significant gravity changes. A large gravity increase (up to 140 μGal) detected across the volcano between June and September 1992 is consistent with the arrival in the accumulation zone of 0.32–2.2×10⁹m³ of new magma, thus favoring continued flank effusion until 1993. A large gravity decrease (200 μGal) in the summit region marked the closing stages of the 1991–1993 event and is associated with magma drainage from the upper levels of Etna's central feeding system. Received: 15 July 1995 / Accepted: 27 October 1997     

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     

Petrology and sulfur and chlorine emissions of the 1963 eruption of Gunung Agung,Bali, Indonesia

 The 1963 eruption of Gunung Agung produced 0.95 km³ dense rock equivalent (DRE) of olivine±hornblende-bearing, weakly phyric, basaltic andesite tephra and lava. Evidence for magma mixing in the eruptive products includes whole-rock compatible and incompatible trace element trends, reverse and complex compositional zoning of mineral phases, disequilibrium mineral assemblages, sieve-textured plagioclase phenocrysts, and augite rims on reversely zoned orthopyroxene. Basalt magma mixed with pre-existing andesite magma shortly before eruption to yield basaltic andesite with a temperature of 1040–1100  °C at an assumed pressure of 2 kb, f O₂>NNO, and an average melt volatile content (H₂O±CO₂) of 4.3 wt.%. Magma-mixing end members may have provided some of the S and Cl emitted in the eruption. Glass inclusions in phenocrysts contain an average of 650 ppm S and 3130 ppm Cl as compared with 70 ppm and 2220 ppm, respectively, in the matrix glass. Maximum S and Cl contents of glass inclusions approach 1800 and 5000 ppm, respectively. Application of the petrologic method to products of the 1963 eruption for estimating volatile release yields of 2.5×10¹²g (Mt) of SO₂ and 3.4 Mt of Cl released from the 0.65 km³ of juvenile tephra which contributed to stratospheric injection of H₂SO₄ aerosols on 17 March and 16 May, when eruption column heights exceeded 20 km above sea level. An independent estimate of SO₂ release from atmospheric aerosol loading (11–12 Mt) suggests that approximately 7 Mt of SO₂ was injected into the stratosphere. The difference between the two estimates can be most readily accounted for by the partitioning of S, as well as some Cl, from the magma into a water-rich vapor phase which was released upon eruption. For other recent high-S-release eruptions of more evolved and oxidized magmas (El Chichón, Pinatubo), the petrologic method gives values two orders of magnitude less than independent estimates of SO₂ emissions. Results from this study of the Agung 1963 magma and its volatile emissions, and from related studies on eruptions of more mafic magmas, suggest that SO₂ emissions from eruptions of higher-S-solubility magma may be more reliably estimated by the petrologic method than may those from more-evolved magma eruptions. Received: 29 June 1994 / Accepted: 25 April 1996     

Degassing during magma ascent in the Mule Creek vent (USA)

 The structures and textures of the rhyolite in the Mule Creek vent (New Mexico, USA) indicate mechanisms by which volatiles escape from silicic magma during eruption. The vent outcrop is a 300-m-high canyon wall comprising a section through the top of a feeder conduit, vent and the base of an extrusive lava dome. Field relations show that eruption began with an explosive phase and ended with lava extrusion. Analyses of glass inclusions in quartz phenocrysts from the lava indicate that the magma had a pre-eruptive dissolved water content of 2.5–3.0 wt% and, during eruption, the magma would have been water-saturated over the vertical extent of the present outcrop. However, the vesicularity of the rhyolite is substantially lower than that predicted from closed-system models of vesiculation under equilibrium conditions. At a given elevation in the vent, the volume fraction of primary vesicles in the rhyolite increases from zero close to the vent margin to values of 20–40 vol.% in the central part. In the centre the vesicularity increases upward from approximately 20 vol.% at 300 m below the canyon rim to approximately 40 vol.% at 200 m, above which it shows little increase. To account for the discrepancy between observed vesicularity and measured water content, we conclude that gas escaped during ascent, probably beginning at depths greater than exposed, by flow through the vesicular magma. Gas escape was most efficient near the vent margin, and we postulate that this is due both to the slow ascent of magma there, giving the most time for gas to escape, and to shear, favouring bubble coalescence. Such shear-related permeability in erupting magma is supported by the preserved distribution of textures and vesicularity in the rhyolite: Vesicles are flattened and overlapping near the dense margins and become progressively more isolated and less deformed toward the porous centre. Local zones have textures which suggest the coalescence of bubbles to form permeable, collapsing foams, implying the former existence of channels for gas migration. Local channelling of gas into the country rocks is suggested by the presence of sub-horizontal syn-eruptive rhyolitic tuffisite veins which depart from the vent margin and invade the adjacent country rock. In the central part of the vent, similar local channelling of gas is indicated by steep syn-eruption tuffisite veins which cut the rhyolite itself. We conclude that the suppression of explosive eruption resulted from gas separation from the ascending magma and vent structure by shear-related porous flow and channelling of gas through tuffisite veins. These mechanisms of gas loss may be responsible for the commonly observed transition from explosive to effusive behaviour during the eruption of silicic magma. Received: 24 May 1995 / Accepted: 13 March 1996     

Shallow and deep reservoirs involved in magma supply of the 1944 eruption of Vesuvius

 During the 1944 eruption of Vesuvius a sudden change occurred in the dynamics of the eruptive events, linked to variations in magma composition. K-phonotephritic magmas were erupted during the effusive phase and the first lava fountain, whereas the emission of strongly porphyritic K-tephrites took place during the more intense fountain. Melt inclusion compositions (major and volatile elements) highlight that the magmas feeding the eruption underwent differentiation at different pressures. The K-tephritic volatile-rich melts (up to 3 wt.% H₂O, 3000 ppm CO₂, and 0.55 wt.% Cl) evolved to reach K-phonotephritic compositions by crystallization of diopside and forsteritic olivine at total fluid pressure higher than 300 MPa. These magmas fed a very shallow reservoir. The low-pressure differentiation of the volatile-poor K-phonotephritic magmas (H₂O<1 wt.%) involved mixing, open-system degassing, and crystallization of leucite, salite, and plagioclase. The eruption was triggered by intrusion of a volatile-rich magma batch that rose from a depth of 11–22 km into the shallow magma chamber. The first phase of the eruption represents the partial emptying of the shallow reservoir, the top of which is within the volcanic edifice. The newly arrived magma mixed with that resident in the shallow reservoir and forced the transition from the effusive to the lava fountain phase of the eruption. Received: 14 September 1998 / Accepted: 10 January 1999     

The contemporaneous emission of low-K and high-K trachybasalts and the role of the NE Rift during the 2002 eruptive event,Mt. Etna,Italy

Mount Etna volcano erupted almost simultaneously on its northeastern and southern flanks between October 27 and November 3, 2002. The eruption on the northeastern flank lasted for 8 days, while on the southern flank it continued for 3 months. The northeastern flank eruption was characterized by the opening of a long eruptive fracture system between 2,900 and 1,900 m.a.s.l. A detailed survey indicates that the fractures’ direction shifted during the opening from N10W (at the NE Crater, 2,900 m) to N45E (at its lowest portion, 1,900 m) and that distinct magma groups were erupted at distinct fracture segments. Based on their petrological features, three distinct groups of rocks have been identified. The first group, high-potassium porphyritic (HKP), is made up of porphyritic lavas with a Porphyritic Index (P.I.) of 20–32 and K₂O content higher than 2 wt%. The second group is represented by lavas and tephra with low modal phenocryst abundance (P.I. < 20) named here oligo-phyric (low-phyric), and K₂O content higher than 2 wt% (HKO, high-potassium oligophyric). The third group, low-potassium oligophyric (LKO), consists of tephra with oligophyric texture (P.I. < 20) but K₂O content < 2 wt%. K-rich magmas (HKP and HKO) are similar to the magma erupted on the southern flank, and geochemical variations within these groups can be accounted for by a variable degree of fractionation from a single parent magma. The K-poor magma (LKO), erupted only in the upper segment of the fracture, cannot be placed on the same liquid line of descent of the HK groups, and it is similar to the magmas that fed the activity of Etna volcano prior to the eruption of 1971. This is the first time since then that a magma of this composition has been documented at Mt. Etna, thus providing a strong indication for the existence of distinct batches of magma whose rise and differentiation are independent from the main conduit system. The evolution of this eruption provides evidence that the NE Rift plays a very active role in the activity of Mt. Etna volcano, and that its extensional tectonics allows the intrusion and residence of magma bodies at various depths, which can therefore differentiate independently from the main open conduit system.     

Volcanological and petrological evolution of Nyiragongo volcano,Virunga volcanic field,Zaire

Three major phases are distinguished during the growth of Nyiragongo, an active volcano at the western limit of the Virunga Range, Zaire. Lavas erupted during phase 1 are strongly undersaturated melilitites characterized by the presence of kalsilite phenocrysts, perovskite, and the abundance of calcite in the matrix. Such lavas crop out mainly on the inner crater wall and progressively evolve toward more aphyric melilite nephelinites well represented on the flanks of the volcano. Adventive vents lying at the base of the cone developed along radial fracture systems and erupted olivine and/or clinopyroxene – rich melilitites or nephelinites. Stage 2 lavas are melilite-free nephelinites. Clinopyroxene is the main phenocryst and feldspathoids are abundant in the lavas exposed on the crater wall. These flows result from periodic overflowing of a magma column from an open crater. Extensive fissure flows which erupted from the base of the cone at the end of this stage are related to widespread draining out of magma which in turn induces the formation of the summit pit crater. Magmas erupted during stage 3 are relatively aphyric melilite nephelinites and the main volcanological characteristic is the permanent lava lake observed into the pit crater until the 1977 eruption. Fluctuations of the level of the lava lake was responsible for the development of the inner terraces. Periodic overflowing of the lava lake from the central pit formed the nepheline aggregate lava flows. Petrography and major element geochemistry allow the determination of the principal petrogenetic processes. Melilitites and nephelinites erupted from the summit crater are lavas derived, via clinopyroxene fractionation, from a more primitive melt. The abundance of feldspathoids in these lavas is in keeping with nepheline flotation. Aphyric melilite nephelinites covering the flanks and the extensive fissure flows have a homogeneous chemical composition; rocks from the historical lava lake are slightly more evolved. All these lavas differentiated in a shallow reservoir. Lavas erupted from the parasitic vents are mainly olivine and/or clinopyroxene-phyric rocks. Rushayite and picrites from Muja cone are peculiar high-magnesium lavas resulting from the addition of olivine xenocrysts to melilitic or nephelinitic melts. Fluid and melt inclusions in olivine and clinopyroxene phenocrysts indicate a crystallization depth of 10–14 km. A model involving two reservoirs located at different depths and periodically connected is proposed to explain the petrography of the lavas; this hypothesis is in accordance with geophysical data. Received: July 8, 1993/Accepted: September 10, 1993     

Petrological characteristics and volatile content of magma from the 2000 eruption of Miyakejima Volcano, Japan

Among the series of eruptions at Miyakejima volcano in 2000, the largest summit explosion occurred on 18 August 2000. During this explosion, vesiculated bombs and lapilli having cauliflower-like shapes were ejected as essential products. Petrological observation and chemical analyses of the essential ejecta and melt inclusions were carried out in order to investigate magma ascent and eruption processes. SEM images indicate that the essential bombs and lapilli have similar textures, which have many tiny bubbles, crystal-rich and glass-poor groundmass and microphenocrysts of plagioclase, augite and olivine. Black ash particles, which compose 40% of the air-fall ash from the explosion, also have similar textures to the essential bombs. Whole-rock analyses show that the chemical composition of all essential ejecta is basaltic (SiO₂=51–52 wt%). Chemical analyses of melt inclusions in plagioclase and olivine phenocrysts indicate that melt in the magma had 0.9–1.9 wt% H₂O, <0.011 wt% CO₂, 0.04–0.17 wt% S and 0.06–0.1 wt% Cl. The variation in volatile content suggests degassing of the magma during ascent up to a depth of about 1 km. The ratio of H₂O and S content of melt inclusions is similar to that of volcanic gas, which has been intensely and continuously emitted from the summit since the end of August 2000, indicating that the 18 August magma is the source of the gas emission. Based on the volatile content of the melt inclusions and the volcanic gas composition, the initial bulk volatile content of the magma was estimated to be 1.6–1.9 wt% H₂O, 0.08–0.1 wt% CO₂, 0.11–0.17 wt% S and 0.06–0.07 wt% Cl. The basaltic magma ascended from a deeper chamber (10 km) due to decrease in magma density caused by volatile exsolution with pressure decrease. The highly vesiculated magma, which had at least 30 vol% bubbles, may have come into contact with ground water at sea level causing the large explosion of 18 August 2000.Editorial responsibility: S. Nakada, T. DuittAn erratum to this article can be found at     

A study of melt inclusions at Vulcano (Aeolian Islands, Italy): insights on the primitive magmas and on the volcanic feeding system

 This work presents the results of a microthermometric and EPMA-SIMS study of melt inclusions in phenocrysts of rocks of the shoshonitic eruptive complex of Vulcano (Aeolian Islands, Italy). Different primitive magmas related to two different evolutionary series, an older one (50–25 ka) and a younger one (15 ka to 1890 A.D.), were identified as melt inclusions in olivine Fo_88–91 crystals. Both are characterized by high Ca/Al ratio and present very similar Rb/Sr, B/Be and patterns of trace elements, with Nb and Ti anomalies typical of a subduction zone. The two basalts present the same temperature of crystallization (1180±20  °C) and similar volatile abundances. The H₂O, S and Cl contents are relatively high, whereas magmatic CO₂ concentrations are very low, probably due to CO₂ loss before low-pressure crystallization and entrapment of melt inclusions. The mineral chemistry of the basaltic assemblages and the high Ca/Al ratio of melt inclusions indicate an origin from a depleted, metasomatized clinopyroxene-rich peridotitic mantle. The younger primitive melt is characterized with respect to the older one by higher K₂O and incompatible element abundances, by lower Zr/Nb and La/Nb, and by higher Ba/Rb and LREE enrichment. A different degree of partial melting of the same source can explain the chemical differences between the two magmas. However, some anomalies in Sr, Rb and K contents suggest either a slightly different source for the two magmas or differing extents of crustal contamination. Low-pressure degassing and cooling of the basaltic magmas produce shoshonitic liquids. The melt inclusions indicate evolutionary paths via fractional crystallization, leading to trachytic compositions during the older activity and to rhyolitic compositions during the recent one. The bulk-rock compositions record a more complex history than do the melt inclusions, due to the syneruptive mixing processes commonly affecting the magmas erupted at Vulcano. The composition and temperature data on melt inclusions suggest that in the older period of activity several shallow magmatic reservoirs existed; in the younger one a relatively homogeneous feeding system is active. The shallow magmatic reservoir feeding the recent eruptive activity probably has a vertical configuration, with basaltic magma in the deeper zones and differentiated magmas in shallower, low-volume, dike-like reservoirs. Received: 11 March 1998 / Accepted: 14 July 1998     

Grain-size variation of tephra derived from volcanic umbrella clouds

The relationships among the thickness and grain-size of tephra-fall deposits and the volumetric flow rate of their source umbrella clouds are analytically obtained. The logarithm of the ratio of the probability distribution function based on grain size (ln R _f) in fall deposits at two localities from the vent (r ₁ and r ₂, respectively) has a linear relationship with the particle-settling velocity, v, as: where Q is the volumetric flow rate of the umbrella cloud and A is a constant for a given pair of localities. The volumetric flow rate of the umbrella cloud can be estimated from granulometric data using this formula. Generally, the thickness–distance relationship of tephra-fall deposits depends on the initial grain-size distribution and the volumetric flow rate of the umbrella cloud. The empirical relationship of the exponential thinning behaviour can be extrapolated towards infinite distance only for a specific initial grain size which is similar to a log-normal distribution with σ _φ=2.5, otherwise it holds only in a limited range of distances. In applying these results to the 1991 eruption of Mt. Pinatubo, it is shown that the volumetric flow rate of the umbrella cloud during the climactic phase of 15 June was approximately 5×10¹⁰ m³/s, which is fairly consistent with the expansion rate of the umbrella cloud observed in the satellite images. Received March 20, 1993/Accepted September 11, 1993     

Magma degassing and basaltic eruption styles: a case study of 2000 year BP Xitle volcano in central Mexico

To investigate the relationship between volatile abundances and eruption style, we have analyzed major element and volatile (H₂O, CO₂, S) concentrations in olivine-hosted melt inclusions in tephra from the 2000 yr BP eruption of Xitle volcano in the central Trans-Mexican Volcanic Belt. The Xitle eruption was dominantly effusive, with fluid lava flows accounting for 95% of the total dense rock erupted material (1.1 km³). However, in addition to the initial, Strombolian, cinder cone-building phase, there was a later explosive phase that interrupted effusive activity and deposited three widespread ash fall layers. Major element compositions of olivine-hosted melt inclusions from these ash layers range from 52 to 58 wt.% SiO₂, and olivine host compositions are Fo_84–86. Water concentrations in the melt inclusions are variable (0.2–1.3 wt.% H₂O), with an average of 0.45±0.3 (1σ) wt.% H₂O. Sulfur concentrations vary from below detection (50 ppm) to 1000 ppm but are mostly ≤200 ppm and show little correlation with H₂O. Only the two inclusions with the highest H₂O have detectable CO₂ (310–340 ppm), indicating inclusion entrapment at higher pressures (700–900 bars) than for the other inclusions (≤80 bars). The low and variable H₂O and S contents of melt inclusions combined with the absence of less soluble CO₂ indicates shallow-level degassing before olivine crystallization and melt inclusion formation. Olivine morphologies are consistent with the interpretation that most crystallization occurred rapidly during near-surface H₂O loss. During cinder cone eruptions, the switch from initial explosive activity to effusive eruption probably occurs when the ascent velocity of magma becomes slow enough to allow near-complete degassing of magma at shallow depths within the cone as a result of buoyantly rising gas bubbles. This allows degassed lavas to flow laterally and exit near the base of the cone while gas escapes through bubbly magma in the uppermost part of the conduit just below the crater. The major element compositions of melt inclusions at Xitle show that the short-lived phase of renewed explosive activity was triggered by a magma recharge event, which could have increased overpressure in the storage reservoir beneath Xitle, leading to increased ascent velocities and decreased time available for degassing during ascent.     

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     

Experimental simulations of gas-driven eruptions: kinetics of bubble growth and effect of geometry

 Simulated gas-driven eruptions using CO₂–water-polymer systems are reported. Eruptions are initiated by rapidly decompressing CO₂–saturated water containing up to 1.0 wt.% CO₂. Both cylindrical test cells and a flask test cell were used to examine the effect of magma chamber/conduit geometry on eruption dynamics. Bubble-growth kinetics are examined quantitatively in experiments using cylindrical test cells. Uninhibited bubble growth can be roughly expressed as dr/dt≈λD(β-1)/(γt ^1/3) for a CO₂–water-polymer system at 0–22  °C and with viscosities up to 5 Pa·s, where r is the radius of bubbles, λ and D are the Ostwald solubility coefficient and diffusivity of the gas in the liquid, β is the degree of saturation (decompression ratio), and γ characterizes how the boundary layer thickness increases with time and is roughly 1.0×10^–5 m/s^1/3 in this system. Unlike the radius of cylindrical test cells, which does not affect the eruption threshold and dynamics, the shape of the test cells (flask vs cylindrical) affects the dynamics but not the threshold of eruptions. For cylindrical test cells, the front motion is characterized by constant acceleration with both Δh (the height increase) and ΔV (the volume increase) being proportional to t ²; for the flask test cell, however, neither Δh nor ΔV is proportional to t ² as the conduit radius varies. Test-cell geometry also affects foam stability. In the flask test cell, as it moves from the wider base chamber into the narrower conduit, the bubbly flow becomes fragmented, affecting the eruption dynamics. The fragmentation may be caused by a sudden increase in acceleration induced by conduit-shape change, or by the presence of obstacles to the bubbly flow. This result may help explain the range in vesicularities of pumice and reticulite. Received: 16 May 1997 / Accepted: 11 October 1997     

Ignimbrites of basaltic andesite and andesite compositions from Tanna,New Hebrides Arc

Tanna island is part of a large volcanic complex mainly subsided below sea-level. On-land, two series of hydroclastic deposits and ignimbrites overlie the subaerial remains of a basal, mainly effusive volcano. The ‘Older’ Tanna Ignimbrite series (OTI), Late Pliocene or Pleistocene in age, consists of ash flows and ash- and scoria-flow deposits associated with fallout tephra layers, overlain by indurated pumice-flow deposits. Phreatomagmatic features are a constant characteristic of these tuffs. The ‘younger’ Late Pleistocene pyroclastics, the Siwi sequence, show basal phreatomagmatic deposits overlain by two successive flow units, each comprising a densely welded layer and a nonwelded ash-flow deposit. Whole-rock analyses of 17 juvenile clasts from the two sequences (vitric blocks from the phreatomagmatic deposits, welded blocks, scoriaceous bombs and pumices from the ignimbrites) show basaltic andesite and andesite compositions (SiO₂=53–60%). In addition, 296 microprobe analyses of glasses in these clasts show a wide compositional range from 51 to 69% SiO₂. Dominant compositions at ∼54, 56, 58.5 and 61–62% SiO₂ characterize the glass from the OTI. Glass compositions in the lower – phreatomagmatic – deposits from the Siwi sequence also show multimodal distribution, with peaks at SiO₂=55, 57.5, 61–62 and 64% whereas the upper ignimbrite has a predominant composition at 61–62% SiO₂. In both cases, mineralogical data and crystal fractionation models suggest that these compositions represent the magmatic signature of a voluminous layered chamber, the compositional gradient of which is the result of fractional crystallization. During two major eruptive stages, probably related to two caldera collapses, the OTI and Siwi ignimbrites represent large outpourings from these magmatic reservoirs. The successive eruptive dynamics, from phreatomagmatic to Plinian, emphasize the role of water in initiating the eruptions, without which the mafic and intermediate magmas probably would not have erupted. Received: February 19, 1993/Accepted October 10, 1993     

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