Influence of pre-eruptive storage conditions and volatile contents on explosive Plinian style eruptions of basic magma

Sub-Plinian to Plinian eruptions of basic magma present a challenge to modeling volcanic behavior because many models rely on magma becoming viscous enough during ascent to behave brittlely and cause fragmentation. Such models are unable, however, to strain low viscosity magma fast enough for it to behave brittlely. That assumes that such magmas actually have low viscosities, but the rare Plinian eruptions of basic magma may in fact result from them being anomalously viscous. Here, we examine two such eruptions, the 122 B.C. eruption of hawaiitic basalt from Mt. Etna and the late Pleistocene eruption of basaltic andesite from Masaya Caldera, to test whether they were anomalously viscous. We carried out hydrothermal experiments on both magmas and analyzed glass inclusions in plagioclase phenocrysts from each to determine their most likely pre-eruptive temperatures and water contents. We find that the hawaiite was last stored at 1,000–1,020°C, whereas the basaltic andesite was last stored at 1,010–1,060°C, and that both were water saturated with ∼3.0 wt.% water dissolved in them. Such water contents are not high enough to trigger Plinian explosive behavior, as much more hydrous basic magmas erupt less violently. In addition, despite being relatively cool, the viscosities of both magmas would range from ∼10^2.2–2.5 Pa s before erupting to ∼10⁴ Pa s when essentially degassed, all of which are too fluid to cause brittle disruption. Without invoking special external forces to explain all such eruptions, one of the more plausible explanations is that when the bubble content reaches some critical value the fragile foam-like magma disrupts. The rarity of Plinian eruptions of basic magma may be because such magmas must ascend fast enough to retain their bubbles.     

Holocene tephrochronology record of large explosive eruptions in the southernmost Patagonian Andes

For regionally widespread Holocene tephra layers in southernmost Patagonia, correlations based on both chemical and chronological data indicate their derivation from five large-volume (>1 km³) explosive eruptions of four different volcanoes in the southernmost Andes. Bulk-tephra and tephra-glass major and trace-element chemistry and Sr isotopic ratios unambiguously distinguish different source volcanoes, and imply that two of the regionally widespread tephra (MB₁ and MB₂) were derived from Mt. Burney (52°S), one (R₁) from Reclus (51°S), one (A₁) from Aguilera (50°S) and one (H₁) from Hudson volcano (46°S). The H₁ tephra derived from the Hudson volcano, which is located at the southern end of the Andean Southern Volcanic Zone (SVZ; 33–46°S), contains distinctive greenish andesitic glass with FeO > 4.5 wt.% and TiO₂ > 1.2 wt.%. In contrast, rhyolitic glass in tephra derived from the eruptions of Mt. Burney, Reclus and Aguilera volcanoes, which are located in the Andean Austral Volcanic Zone (AVZ; 49–55°S), is clear and transparent and has significantly lower FeO and TiO₂. Tephra derived from these three AVZ volcanoes all contain plagioclase, orthopyroxene, minor clinopyroxene and amphibole. Biotite occurs only in the Aguilera A₁ tephra, which also has the highest bulk-tephra and tephra-glass K₂O and Rb contents. Averages of new and published ¹⁴C ages determined on organic material in soil and sediment samples above and below these tephra constrain the uncalibrated ¹⁴C age of the R₁ eruption of Reclus volcano to 12,685 ± 260 years BP, the MB₁ and MB₂ eruptions of Mt. Burney to 8,425 ± 500 and 3,830 ± 390 years BP, the Hudson H₁ eruption to 6,850 ± 160 years BP, and the A₁ eruption of Aguilera volcano to 3,000 ± 100 years BP. The volume of the largest of these eruptions, H₁ of the Hudson volcano, is estimated as >18 km³. The volume of the Reclus R₁ eruption is estimated at >10 km³, the Aguilera A₁ eruption at between 4 and 9 km³, and the younger Mt. Burney MB₂ eruption at ≥2.8 km³. The volume of the older MB₁ Mt. Burney eruption is the least well constrained, but must have been larger than the younger MB₂ eruption. The data indicate that the frequency of explosive activity of volcanic centers in the AVZ is lower than in the southern SVZ.     

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     

The impact of moderate-scale explosive eruptions on stratospheric gas injections

Volcanic gases such as SO ₂, H ₂S, HCl and COS emitted during explosive eruptions significantly affect atmospheric chemistry and therefore the Earth's climate. We have evaluated the dependence of volcanic gas emission into the atmosphere on altitude, latitude, and tectonic setting of volcanoes and on the season in which eruptions occurred. These parameters markedly influence final stratospheric gas loading. The latitudes and altitudes of 360 active volcanoes were compared to the height of the tropopause to calculate the potential quantity of volcanic gases injected into the stratosphere. We calculated a possible stratospheric gas loading based on different volcanic plume heights (6, 10, and 15 km) generated by moderate-scale explosive eruptions to show the importance of the actual plume height and volcano location. At a plume height of 15 km for moderate-scale explosive eruptions, a volcano at sea level can cause stratospheric gas loading because the maximum distance to the tropopause is 15–16 km in the equatorial region (0–30°). Eruptions in the tropics have to be more powerful to inject gas into the stratosphere than eruptions at high latitudes because the tropopause rises from ca. 9–11 km at the poles to 15–16 km in the equatorial region (0–30°N and S). The equatorial region is important for stratospheric gas injection because it is the area with the highest frequency of eruptions. Gas injected into the stratosphere in equatorial areas may spread globally into both hemispheres.     

New radiometric and petrological constraints on the evolution of the Pichincha volcanic complex (Ecuador)

Fieldwork, radiometric (⁴⁰Ar/³⁹Ar and ¹⁴C) ages and whole-rock geochemistry allow a reconstruction of eruptive stages at the active, mainly dacitic, Pichincha Volcanic Complex (PVC), whose eruptions have repeatedly threatened Quito, most recently from 1999 to 2001. After the emplacement of basal lavas dated at ∼1100 to 900 ka, the eruptive activity of the old Rucu Pichincha volcano lasted from ∼850 ka to ∼150 ka before present (BP) and resulted in a 15 × 20 km-wide edifice, which comprises three main building stages: (1) A lower stratocone (Lower Rucu, ∼160 km³ in volume) developed from ∼850 to 600 ka; (2) This edifice was capped by a steeper-sided and less voluminous cone (the Upper Rucu, 40–50 km³), the history of which started 450–430 ka ago and ended around 250 ka with a sector collapse; (3) A smaller (8–10 km³) but more explosive edifice grew in the avalanche amphitheatre and ended Rucu Pichincha's history about 150 ka ago. The Guagua Pichincha volcano (GGP) was developed from 60 ka on the western flank of Rucu with four growth stages separated by major catastrophic events. (1) From ∼60 to 47 ka, a basal effusive stratocone developed, terminating with a large ash-and-pumice flow event. (2) This basal volcano was followed by a long-lasting dome building stage and related explosive episodes, the latter occurring between 28–30 and 22–23 ka. These first two stages formed the main GGP (∼30 km³), a large part of which was removed by a major collapse 11 ka BP. (3) Sustained explosive activity and viscous lava extrusions gave rise to a new edifice, Toaza (4–5 km³ in volume), which in turn collapsed around 4 ka BP. (4) The ensuing amphitheatre was partly filled by the ∼1-km³ Cristal dome, which is the historically active centre of the Pichincha complex. The average output rate for the whole PVC is 0.29 km³/ka. Nevertheless, the chronostratigraphic resolution we obtained for Lower Rucu Pichincha and for the two main edifices of Guagua Pichincha (main GGP and Toaza), leads to eruptive rates of 0.60–0.65 km³/ka during these construction stages. These output rates are compared to those of other mainly dacitic volcanoes from continental arcs. Our study also supports an overall SiO₂ and large-ion lithophile elements enrichment as the PVC develops. In particular, distinctive geochemical signatures indicate the involvement of a new magma batch at the transition between Rucu and Guagua. At the GGP, the same phenomenon occurs at each major collapse event marking the onset of the ensuing magmatic stage. Since the 11-ka-BP collapse event, this magmatic behaviour has led to increasingly explosive activity. Four explosive cycles of between 100 and 200 years long have taken place at the Cristal dome in the past 3.7 ka, and repose intervals between these cycles have tended to decrease with time. As a consequence, we suggest that the 1999–2001 eruptive period may have initiated a new eruptive cycle that might pose a future hazard to Quito (∼2 million inhabitants).     

Kīlauea summit overflows: their ages and distribution in the Puna District, Hawai'i

 The tube-fed pāhoehoe lava flows covering much of the northeast flank of Kīlauea Volcano are named the 'Ailā'au flows. Their eruption age, based on published and six new radiocarbon dates, is approximately AD 1445. The flows have distinctive paleomagnetic directions with steep inclinations (40°–50°) and easterly declinations (0°–10°E). The lava was transported ∼40 km from the vent to the coast in long, large-diameter lava tubes; the longest tube (Kazumura Cave) reaches from near the summit to within several kilometers of the coast near Kaloli Point. The estimated volume of the 'Ailā'au flow field is 5.2±0.8 km³, and the eruption that formed it probably lasted for approximately 50 years. Summit overflows from Kīlauea may have been nearly continuous between approximately AD 1290 and 1470, during which time a series of shields formed at and around the summit. The 'Ailā'au shield was either the youngest or the next to youngest in this series of shields. Site-mean paleomagnetic directions for lava flows underlying the 'Ailā'au flows form only six groups. These older pāhoehoe flows range in age from 2750 to <18,000 BP, and the region was inundated by lava flows only three times in the past 5000 years. The known intervals between eruptive events average ∼1600 years and range from ∼1250 years to >2200 years. Lava flows from most of these summit eruptions also reached the coast, but none appears as extensive as the 'Ailā'au flow field. The chemistry of the melts erupted during each of these summit overflow events is remarkably similar, averaging approximately 6.3 wt.% MgO near the coast and 6.8 wt.% MgO near the summit. The present-day caldera probably formed more recently than the eruption that formed the 'Ailā'au flows (estimated termination ca. AD 1470). The earliest explosive eruptions that formed the Keanakāko'i Ash, which is stratigraphically above the 'Ailā'au flows, cannot be older than this age. Received: 10 October 1998 / Accepted: 12 May 1999     

The Taurewa Eruptive Episode: evidence for climactic eruptions at Ruapehu volcano, New Zealand

 The ca. 10,500 years B.P. eruptions at Ruapehu volcano deposited 0.2–0.3 km³ of tephra on the flanks of Ruapehu and the surrounding ring plain and generated the only known pyroclastic flows from this volcano in the late Quaternary. Evidence of the eruptions is recorded in the stratigraphy of the volcanic ring plain and cone, where pyroclastic flow deposits and several lithologically similar tephra deposits are identified. These deposits are grouped into the newly defined Taurewa Formation and two members, Okupata Member (tephra-fall deposits) and Pourahu Member (pyroclastic flow deposits). These eruptions identify a brief (<ca. 2000-year) but explosive period of volcanism at Ruapehu, which we define as the Taurewa Eruptive Episode. This Episode represents the largest event within Ruapehu's ca. 22,500-year eruptive history and also marks its culmination in activity ca. 10,000 years B.P. Following this episode, Ruapehu volcano entered a ca. 8000-year period of relative quiescence. We propose that the episode began with the eruption of small-volume pyroclastic flows triggered by a magma-mingling event. Flows from this event travelled down valleys east and west of Ruapehu onto the upper volcanic ring plain, where their distal remnants are preserved. The genesis of these deposits is inferred from the remanent magnetisation of pumice and lithic clasts. We envisage contemporaneous eruption and emplacement of distal pumice-rich tephras and proximal welded tuff deposits. The potential for generation of pyroclastic flows during plinian eruptions at Ruapehu has not been previously considered in hazard assessments at this volcano. Recognition of these events in the volcanological record is thus an important new factor in future risk assessments and mitigation of volcanic risk at Tongariro Volcanic Centre. Received: 5 July 1998 / Accepted: 12 March 1999     

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     

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     

Major-ion bulk deposition around an active volcano (Mt. Etna, Italy)

Bulk atmospheric deposition of major cations (Na, K, Ca, Mg) and anions (Cl, F, SO₄) were measured at 15 sites around an active volcano, Mount Etna, from 2001 to 2003. Their composition indicates several natural sources, among which deposition of plume-derived volcanogenic gas compounds is prevalent for F, Cl and S. Plume-derived acidic compounds are also responsible for the prevailing acidic composition of the samples collected on the summit of the volcano (pH in the 2.45–5.57 range). Cation species have complex origin, including deposition of plume volcanogenic ash and aerosols and soil-dust wind re-suspension of either volcanic or carbonate sedimentary rocks. Variation of the deposition rates during the March 2001–March 2003 period, coupled with previous measurements from 1997 to 2000 (Appl Geochem 16:985–1000, 2001), were compared with the variation of SO₂ flux, volcanic activity and rainfall. The deposition rate was mainly controlled by rainfall. Commonly, about 0.1–0.9% of HF, HCl and SO₂ emitted by the summit crater's plume were deposited around the volcano. We estimate that ∼2 Gg of volcanogenic sulphur were deposited over the Etnean area during the 2002–2003 flank eruption, at an average rate of ∼24 Mg day⁻¹ which is two orders of magnitude higher than that typical of quiescent degassing phases.     

Late Holocene eruptive activity at Nevado Cayambe Volcano, Ecuador

 Four Late Holocene pyroclastic units composed of block and ash flows, surges, ashfalls of silicic andesite and dacite composition, and associated lahar deposits represent the recent products emitted by domes on the upper part of Nevado Cayambe, a large ice-capped volcano 60 km northeast of Quito. These units are correlated stratigraphically with fallout deposits (ash and lapilli) exposed in a peat bog. Based on ¹⁴C dating of the peat and charcoal, the following ages were obtained: ∼910 years BP for the oldest unit, 680–650 years BP for the second, and 400–360 years BP for the two youngest units. Moreover, the detailed tephrochronology observed in the peat bog and in other sections implies at least 21 volcanic events during the last 4000 years, comprising three principal eruptive phases of activity that are ∼300, 800, and 900 years in duration and separated by repose intervals of 600–1000 years. The last phase, to which the four pyroclastic units belong, has probably not ended, as suggested by an eruption in 1785–1786. Thus, Cayambe, previously thought to have been dormant for a long time, should be considered active and potentially dangerous to the nearby population of the Interandean Valley. Received: 5 July 1997 / Accepted: 21 October 1997     

Simulating bubble number density of rhyolitic pumices from Plinian eruptions: constraints from fast decompression experiments

Decompression experiments of a crystal-free rhyolitic liquid with ≈ 6.6 wt. % H₂O were carried out at a pressure range from 250 MPa to 30–75 MPa in order to characterize effects of magma ascent rate and temperature on bubble nucleation kinetics, especially on the bubble number density (BND, the number of bubbles produced per unit volume of liquid). A first series of experiments at 800°C and fast decompression rates (10–90 MPa/s) produced huge BNDs (≈ 2 × 10¹⁴ m⁻³ at 10 MPa/s ; ≈ 2 × 10¹⁵ m⁻³ at 90 MPa/s), comparable to those in natural silicic pumices from Plinian eruptions (10¹⁵–10¹⁶ m⁻³). A second series of experiments at 700°C and 1 MPa/s produced BNDs (≈ 9×10¹² m⁻³) close to those observed at 800°C and 1 MPa/s (≈ 6 × 10¹² m⁻³), showing that temperature has an insignificant effect on BNDs at a given decompression rate. Our study strengthens the theory that the BNDs are good markers of the decompression rate of magmas in volcanic conduits, irrespective of temperature. Huge number densities of small bubbles in natural silicic pumices from Plinian eruptions imply that a major nucleation event occurs just below the fragmentation level, at which the decompression rate of ascending magmas is a maximum (≥ 1 MPa/s).     

Geochemistry of the volcano-hydrothermal system of El Chichón Volcano, Chiapas, Mexico

 The 1982 eruption of El Chichón volcano ejected more than 1 km³ of anhydrite-bearing trachyandesite pyroclastic material to form a new 1-km-wide and 300-m-deep crater and uncovered the upper 500 m of an active volcano-hydrothermal system. Instead of the weak boiling-point temperature fumaroles of the former lava dome, a vigorously boiling crater spring now discharges  / 20 kg/s of Cl-rich (∼15 000 mg/kg) and sulphur-poor ( / 200 mg/kg of SO₄), almost neutral (pH up to 6.7) water with an isotopic composition close to that of subduction-type magmatic water (δD=–15‰, δ¹⁸O=+6.5‰). This spring, as well as numerous Cl-free boiling springs discharging a mixture of meteoric water with fumarolic condensates, feed the crater lake, which, compared with values in 1983, is now much more diluted (∼3000 mg/kg of Cl vs 24 030 mg/kg), less acidic (pH=2.6 vs 0.56) and contains much lower amounts of S ( / 200 mg/kg of SO₄, vs 3550 mg/kg) with δ³⁴S=0.5–4.2‰ (+17‰ in 1983). Agua Caliente thermal waters, on the southeast slope of the volcano, have an outflow rate of approximately 100 kg/s of 71  °C Na–Ca–Cl water and are five times more concentrated than before the eruption (B. R. Molina, unpublished data). Relative N₂, Ar and He gas concentrations suggest extensional tectonics for the El Chichón volcanic centre. The ³He/⁴He and ⁴He/²⁰Ne ratios in gases from the crater fumaroles (7.3R_a, 2560) and Agua Caliente hot springs (5.3R_a, 44) indicate a strong magmatic contribution. However, relative concentrations of reactive species are typical of equilibrium in a two-phase boiling aquifer. Sulphur and C isotopic data indicate highly reducing conditions within the system, probably associated with the presence of buried vegetation resulting from the 1982 eruption. All Cl-rich waters at El Chichón have a common source. This water has the appearence of a "partially matured" magmatic fluid: condensed magmatic vapour neutralized by interaction with fresh volcaniclastic deposits and depleted in S due to anhydrite precipitation. Shallow ground waters emerging around the volcano from the thick cover of fresh pumice deposits (Red waters) are Ca–SO₄–rich and have a negative oxygen isotopic shift, probably due to ongoing formation of clay at low temperatures. Received: 21 July 1997 / Accepted: 4 December 1997     

Eddy covariance imaging of diffuse volcanic CO<Subscript>2</Subscript> emissions at Mammoth Mountain,CA, USA

Use of eddy covariance (EC) techniques to map the spatial distribution of diffuse volcanic CO₂ fluxes and quantify CO₂ emission rate was tested at the Horseshoe Lake tree-kill area on Mammoth Mountain, California, USA. EC measurements of CO₂ flux were made during September–October 2010 and ranged from 85 to 1,766 g m⁻² day⁻¹. Comparative maps of soil CO₂ flux were simulated and CO₂ emission rates estimated from three accumulation chamber (AC) CO₂ flux surveys. Least-squares inversion of measured eddy covariance CO₂ fluxes and corresponding modeled source weight functions recovered 58–77% of the CO₂ emission rates estimated based on simulated AC soil CO₂ fluxes. Spatial distributions of modeled surface CO₂ fluxes based on EC and AC observations showed moderate to good correspondence (R ² = 0.36 to 0.70). Results provide a framework for automated monitoring of volcanic CO₂ emissions over relatively large areas.     

Nature and significance of small volume fall deposits at composite volcanoes: Insights from the October 14, 1974 Fuego eruption,Guatemala

The first of four successive pulses of the 1974 explosive eruption of Fuego volcano, Guatemala, produced a small volume (∼0.02 km³ DRE) basaltic sub-plinian tephra fall and flow deposit. Samples collected within 48 h after deposition over much of the dispersal area (7–80 km from the volcano) have been size analyzed down to 8 φ (4 μm). Tephra along the dispersal axis were all well-sorted (σ _φ = 0.25–1.00), and sorting increased whereas thickness and median grain size decreased systematically downwind. Skewness varied from slightly positive near the vent to slightly negative in distal regions and is consistent with decoupling between coarse ejecta falling off the rising eruption column and fine ash falling off the windblown volcanic cloud advecting at the final level of rise. Less dense, vesicular coarse particles form a log normal sub-population when separated from the smaller (Md_φ < 3φ or < 0.125 mm), denser shard and crystal sub-population. A unimodal, relatively coarse (Md_φ = 0.58φ or 0.7 mm σ _φ = 1.2) initial grain size population is estimated for the whole (fall and flow) deposit. Only a small part of the fine-grained, thin 1974 Fuego tephra deposit has survived erosion to the present day. The initial October 14 pulse, with an estimated column height of 15 km above sea level, was a primary cause of a detectable perturbation in the northern hemisphere stratospheric aerosol layer in late 1974 to early 1975. Such small, sulfur-rich, explosive eruptions may substantially contribute to the overall stratospheric sulfur budget, yet leave only transient deposits, which have little chance of survival even in the recent geologic record. The fraction of finest particles (Md_φ = 4–8φ or 4–63 μm) in the Fuego tephra makes up a separate but minor size mode in the size distribution of samples around the margin of the deposit. A previously undocumented bimodal–unimodal–bimodal change in grain size distribution across the dispersal axis at 20 km downwind from the vent is best accounted for as the result of fallout dispersal of ash from a higher subplinian column and a lower “co-pf” cloud resulting from pyroclastic flows. In addition, there is a degree of asymmetry in the documented grain-size fallout pattern which is attributed to vertically veering wind direction and changing windspeeds, especially across the tropopause. The distribution of fine particles (<8 μm diameter) in the tephra deposit is asymmetrical, mainly along the N edge, with a small enrichment along the S edge. This pattern has hazard significance.     

Gas and particle emissions from Soufrière Hills Volcano, Montserrat, West Indies: characterization and health hazard assessment

 The Soufrière Hills Volcano, Montserrat, erupting since 18 July 1995, intensified its degassing in early 1996 with the continuing growth of the lava dome inside the summit crater. During this period of increased activity, between 11 and 18 March 1996, we measured gases and particles within the visible plume to determine whether at that time it posed a health risk to the population of Plymouth, the capital town, which is 5 km southwest (downwind) and was then still occupied. Gravimetric measurements were made of total suspended particles (TSP) and particles having an aerodynamic diameter of less than 10 μm (PM₁₀). Measurements were made of sulphur dioxide (SO₂), hydrochloric acid (HCl), hydrofluoric acid (HF), nitric acid (HNO₃), acetic acid (CH₃COOH), formic acid (HCOOH), and particulate sulphate (SO₄ ^2–), chloride (Cl^–), nitrate (NO₃ ^–), fluoride (F^–), methanesulphonate (CH₃SO₃ ^–), acetate (CH₃COO^–), formate (HCOO^–), ammonium (NH₄ ⁺), sodium (Na⁺) and acidity (H⁺). Trace metals having human health implications [chromium (Cr), nickel (Ni), cobalt (Co), copper (Cu), zinc (Zn), arsenic (As), selenium (Se), cadmium (Cd), tin (Sn), mercury (Hg) and lead (Pb)] were also determined. Mean concentrations of HCl, SO₂ and HF obtained in the town of Plymouth were 14.0, 5.9 and 0.8 ppbv, respectively. Corresponding concentrations in the mixed plume on the crater edge were 533, 168 and 22 ppbv. There were no direct emissions of HNO₃, although nitrate was detected in coarse particles at the source. Higher concentrations of CH₃COOH and HCOOH were measured close to the crater. Mean TSP and PM₁₀ were 64 and 15 μg m^–3 in Plymouth, and 455 and 47 μg m^–3 on the upper volcano slope. Aerosols were highly acidic at the source but rapidly neutralised during transport. Trace metals were enriched in the aerosol relative to crater surface material. The concentrations of the acid gases, sulphur dioxide in particular, and particles were found to be too small to pose a health hazard at the time of these measurements, when relatively modest volcanic activity was occurring. Received: 9 September 1998 / Accepted: 29 August 1999     

Paleomagnetic dating of the most recent silicic eruptive activity at Pantelleria (Strait of Sicily)

We report on the paleomagnetism of ten sites in the products of the most recent silicic eruptive cycle of Pantelleria, Strait of Sicily. Previously radiometrically dated at 5–10 ka, our comparison with proxies for geomagnetic field directions allows us to narrow considerably the time window during which these eruptions occurred. The strongly peralkaline composition causes the magmas to have low viscosities, locally resulting in strong agglutination of proximal fall deposits. This allows successful extraction of paleomagnetic directions from the explosive phases of eruptions. One of our sites was located in the Serra della Fastuca fall deposit, produced by the first explosive event of the eruptive cycle. The other nine sites were located in the most recent explosive (pumice fall and agglutinate from Cuddia del Gallo and Cuddia Randazzo) and effusive (Khaggiar lava) products. The (very similar) paleomagnetic directions gathered from eight internally consistent sites were compared to reference geomagnetic field directions of the last 5–10 ka. Directions from Cuddia del Gallo agglutinate and Khaggiar flows translate into 5.9- to 6.2-ka ages, whereas the Fastuca pumices yield a slightly older age of 6.2–6.8 ka. Hence, the most recent silicic eruptive cycle lasted at most a millennium and as little as a few centuries around 6.0 ka. Paleomagnetically inferred ages are in good agreement with published (and calibrated by us) ¹⁴C dates from paleosols/charcoals sampled below the studied volcanic units, whereas K/Ar data are more scattered and yield ∼30% older ages. Our data show that the time elapsed since the most recent silicic eruptions at Pantelleria is comparable to the quiescence period separating the two latest volcanic cycles.     

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     

The dynamics and thermodynamics of large ash flows

 Ash flow deposits, containing up to 1000 km³ of material, have been produced by some of the largest volcanic eruptions known. Ash flows propagate several tens of kilometres from their source vents, produce extensive blankets of ash and are able to surmount topographic barriers hundreds of metres high. We present and test a new model of the motion of such flows as they propagate over a near horizontal surface from a collapsing fountain above a volcanic vent. The model predicts that for a given eruption rate, either a slow (10–100 m/s) and deep (1000–3000 m) subcritical flow or a fast (100–200 m/s) and shallow (500–1000 m) supercritical flow may develop. Subcritical ash flows propagate with a nearly constant volume flux, whereas supercritical flows entrain air and become progressively more voluminous. The run-out distance of such ash flows is controlled largely by the mass of air mixed into the collapsing fountain, the degree of fragmentation and the associated rate of loss of material into an underlying concentrated depositional system, and the mass eruption rate. However, in supercritical flows, the continued entrainment of air exerts a further important control on the flow evolution. Model predictions show that the run-out distance decreases with the mass of air entrained into the flow. Also, the mass of ash which may ascend from the flow into a buoyant coignimbrite cloud increases as more air is entrained into the flow. As a result, supercritical ash flows typically have shorter runout distances and more ash is elutriated into the associated coignimbrite eruption columns. We also show that one-dimensional, channellized ash flows typically propagate further than their radially spreading counterparts. As a Plinian eruption proceeds, the erupted mass flux often increases, leading to column collapse and the formation of pumiceous ash flows. Near the critical conditions for eruption column collapse, the flows are shed from high fountains which entrain large quantities of air per unit mass. Our model suggests that this will lead to relatively short ash flows with much of the erupted material being elutriated into the coignimbrite column. However, if the mass flux subseqently increases, then less air per unit mass is entrained into the collapsing fountain, and progressively larger flows, which propagate further from the vent, will develop. Our model is consistent with observations of a number of pyroclastic flow deposits, including the 1912 eruption of Katmai and the 1991 eruption of Pinatubo. The model suggests that many extensive flow sheets were emplaced from eruptions with mass fluxes of 10⁹–10¹⁰ kg/s over periods of 10³–10⁵ s, and that some indicators of flow "mobility" may need to be reinterpreted. Furthermore, in accordance with observations, the model predicts that the coignimbrite eruption columns produced from such ash flows rose between 20 and 40 km. Received: 25 August 1995 / Accepted: 3 April 1996     

The size and frequency of the largest explosive eruptions on Earth

A compilation and analysis of the size and frequency of the largest known explosive eruptions on Earth are presented. The largest explosive events are defined to be those eruptions yielding greater than 10¹⁵ kg of products (>150 times the mass of the 1991 eruption of Mt. Pinatubo). This includes all known eruptions with a volcanic explosivity index (VEI) of 8. A total of 47 such events, ranging in age from Ordovician to Pleistocene, are identified, of which 42 eruptions are known from the past 36 Ma. A logarithmic magnitude scale of eruption size is applied, based on erupted mass, to these events. On this scale, 46 eruptions >10¹⁵ kg are defined to be of magnitude M8. There is one M9 event known so far, the Fish Canyon Tuff, with an erupted mass of >10¹⁶ kg and a magnitude of 9.2. Analysis of this dataset indicates that eruptions of size M8 and larger have occurred with a minimum frequency of 1.4 events/Ma in two pulses over the past 36 Ma. On the basis of the activity during the past 13.5 Ma, there is at least a 75% probability of a M8 eruption (>10¹⁵ kg) occurring within the next 1 Ma. There is a 1% chance of an eruption of this scale in the next 460–7,200 years. While the effect of any individual M8 or larger eruption is considerable, the time-averaged impact (i.e., erupted mass×frequency) of the very largest eruptions is small, due to their rarity. The long-term, time-averaged erupted mass flux from magnitude 8 and 9 eruptions is ~10–100 times less than for M7 eruptions; the time-averaged mass eruption rate from M7 eruptions is 9,500 kg s^–1, whereas for M8 and M9 eruptions it is ~70–1,000 kg s^–1. Comparison of the energy release by volcanic eruptions with that due to asteroid impacts suggests that on timescales of <100,000 years, explosive volcanic eruptions are considerably more frequent than impacts of similar energy yield. This has important implications for understanding the risk of extreme events.Editorial responsibility: R. Cioni