The magma budget of Volcan Arenal, Costa Rica from 1968 to 1980

The volume of magma emitted by Volcan Arenal from July 1968 to March 1980 has been calculated to be 304 × 10⁶ m³ (dense rock equivalent). Most of this magma has been emplaced as block lava flows on the western flanks of the volcano following the initial explosive eruptions in 1968. From 1968 to 1973 the volumetric discharge rate of magma decreased from about 3-2 m³ s⁻¹ to about 1 m³ s⁻¹. During a break in activity in late 1973 the site of effusion moved from Crater A to Crater C about 400 m higher. Subsequent effusion was at a lower rate (0.3 m³ s⁻¹) which remained constant for the next six years. Comparison of dry-tilt measurements during this latter period of steady-state effusion with numerical finite-element models of Arenal's elastic response to the evacuation of magma from an underlying reservoir favor a very shallow reservoir (< 2 km depth) to explain the data. However, the constraints imposed by the measured volumes of magma are not compatible with such a reservoir. Instead, it is argued that the steady downward tilting of the volcano's summit was caused by the loading of the western side of the volcano by about 19 × 10⁶ m³ of lava. Surface loading by lava flows may be an important deformational effect at other volcanoes. A system of magma supply involving open conduits (pipes) for the uppermost one kilometer and transitory conduits (cracks) to a crustal reservoir is proposed. This crustal reservoir initially contained a compositionally graded magma which was evacuated from 1968 to 1973. The subsequent abrupt decrease in effusion rate is compatible with the increased magmatic head required to reach Crater C. The constancy of magma composition and effusion rate from 1974 to 1980 implies a homogeneous magma reservoir.     

Eruptive patterns for the volcano piton de la fournaise, Reunion Island

A study of the historic record of activity of Piton de la Fournaise has revealed a cyclic pattern of eruption involving effusion of oceanite lava from major-flank centers every 20–40 years. Calculated volumes of the recent lava flows and pyroclastic ejecta have established an effusion rate of 3.9 m³ s⁻¹ since 1931 and 6.2 m³ s⁻¹ since 1951. Flank eruptions outside the present caldera define a distribution maximum which is expected to correlate with the depth range of a high-level magma reservoir.A model has been constructed which requires replenishment of a high-level magma chamber at a constant rate and regular eruption from summit and minor-flank centers, acting as “safety valves” to the magma chamber; when the magma chamber reaches its maximum expansion, a major-flank outburst of oceanitic lava occurs.The fact that calculated effusion rates are not consistent with radiometric dates implies an increase in effusion volume with time for the volcano.     

Stratigraphy and eruptive dynamics of a pulsating Plinian eruption of Somma-Vesuvius: the Pomici di Mercato (8900 years B.P.)

New volcanological studies allow reconstruction of the eruption dynamics of the Pomici di Mercato eruption (ca 8,900 cal. yr B.P.) of Somma-Vesuvius. Three main Eruptive Phases are distinguished based on two distinct erosion surfaces that interrupt stratigraphic continuity of the deposits, indicating that time breaks occurred during the eruption. Absence of reworked volcaniclastic deposits on top of the erosion surfaces suggests that quiescent periods between eruptive phases were short perhaps lasting only days to weeks. Each of the Eruptive Phases was characterised by deposition of alternating fall and pyroclastic density current (PDC) deposits. The fallout deposits blanketed a wide area toward the east, while the more restricted PDC deposits inundated the volcano slopes. Eruptive dynamics were driven by brittle magmatic fragmentation of a phonolitic magma, which, because of its mechanical fragility, produced a significant amount of fine ash. External water did not significantly contribute either to fragmentation dynamics or to mechanical energy release during the eruption. Column heights were between 18 and 22 km, corresponding to mass discharge rates between 1.4 and 6 × 10⁷ kg s⁻¹. The estimated on land volume of fall deposits ranges from a minimum of 2.3 km³ to a maximum of 7.4 km³. Calculation of physical parameters of the dilute pyroclastic density currents indicates speeds of a few tens of m s⁻¹ and densities of a few kg m⁻³ (average of the lowermost 10 m of the currents), resulting in dynamic pressures lower than 3 kPa. These data suggest that the potential impact of pyroclastic density currents of the Pomici di Mercato eruption was smaller than those of other Plinian and sub-Plinian eruptions of Somma-Vesuvius, especially those of 1631 AD and 472 AD (4–14 kPa), which represent reference values for the Vesuvian emergency plan. The pulsating and long-lasting behaviour of the Pomici di Mercato eruption is unique in the history of large explosive eruptions of Somma-Vesuvius. We suggest an eruptive scheme in which discrete magma batches rose from the magma chamber through a network of fractures. The injection and rise of the different magma batches was controlled by the interplay between magma chamber overpressure and local stress. The intermittent discharge of magma during a large explosive eruption is unusual for Somma-Vesuvius, as well as for other volcanoes worldwide, and yields new insights for improving our knowledge of the dynamics of explosive eruptions.     

The caldera-forming eruption of Ksudach volcano about cal. A.D. 240: the greatest explosive event of our era in Kamchatka, Russia

The largest Plinian eruption of our era and the latest caldera-forming eruption in the Kuril-Kamchatka region occurred about cal. A.D. 240 from the Ksudach volcano. This catastrophic explosive eruption was similar in type and characteristics to the 1883 Krakatau event. The volume of material ejected was 18–19 km³ (8 km³ DRE), including 15 km³ of tephra fall and 3–4 km³ of pyroclastic flows. The estimated height of eruptive column is 22–30 km. A collapse caldera resulting from this eruption was 4 × 6.5 km in size with a cavity volume of 6.5–7 km³. Tephra fall was deposited to the north of the volcano and reached more than 1000 km. Pyroclastic flows accompanied by ash-cloud pyroclastic surges extended out to 20 km. The eruption was initially phreatomagmatic and then became rhythmic, with each pulse evolving from pumice falls to pyroclastic flows. Erupted products were dominantly rhyodacite throughout the eruption. During the post-caldera stage, when the Shtyubel cone started to form within the caldera, basaltic-andesite and andesite magma began to effuse. The trigger for the eruption may have been an intrusion of mafic magma into the rhyodacite reservoir. The eruption had substantial environmental impact and may have produced a large acidity peak in the Greenland ice sheet.     

Simple stochastic modelling of the eruption history of a basaltic volcano: Nyamuragira,Zaire

Three simple models of the behaviour of a series of basaltic eruptions have been tested against the eruptive history of Nyamuragira. The data set contains the repose periods and the volumes of lava emitted in 22 eruptions since 1901. Model 1 is fully stochastic and eruptions of any volume with random repose intervals are possible. Models 2 and 3 are constrained by deterministic limits on the maximum capacity of the magma reservoir and on the lowest drainage level of the reservoir respectively. The method of testing these models involves (1) seeking change points in the time series to determine regimes of uniform magma supply rate, and (2) applying linear regression to these regimes, which for models 2 and 3 are the determinsstic limits to those models. Two change points in the time series for Nyamuragira, in 1958 and 1980, were determined using a Kolmogorov-Smirnov technique. The latter change involved an increase in the magma supply rate by a factor of 2.5, from 0.55 to 1.37 m³s^-1. Model 2 provides the best fit to the behavior of Nyamuragira with the ratio of variation explained by the model to total variation. R², being greater than 0.9 for all three regimes. This fit can be interpreted to mean that there is a determinstic limit to the elastic strength of the magma reservoir 4–8 km below the summit of the volcano.     

The October 1902 plinian eruption of Santa Maria volcano, Guatemala

The October, 1902, eruption of Santa Maria Volcano, Guatemala, was one of the largest this century. It was preceded by a great earthquake on April 19 centered at the volcano, as well as numerous other major earthquakes. The 18–20 hour-long plinian eruption on October 25 produced a column at least 28 km high, reaching well into the stratosphere.The airfall pumice deposit covered more than 1.2 million km² with a trace of ash and was only two meters thick at the vent. White dacitic pumice, dark gray scoriaceous basalt (with physically and chemically mixed intermediate pumice) and loose crystals of plagioclase, hornblende, hypersthene, biotite and magnetite make up the juvenile components of the deposit. Lithic fragments are of volcanic, plutonic, and metamorphic origin. The plinian deposit is a fine-grained, crystal-rich, single pumice fall unit and shows inverse grading. Mapping of the deposit gives a volume of 8.3 km³ within the one mm isopach. Crystal concentration studies show that the true volume erupted was at least 20 km³ (equivalent to 8.5 km³ of dense dacite) and that 90% of the ejecta was less than 2 mm in diameter.The plinian volume eruption rate averaged 1.2 × 10⁵ m³s⁻¹ and the average gas muzzle velocity of the column exceeded 270 ms⁻¹. A total of 8.3 × 10¹⁸ J of energy were released by the eruption. A knowledge of both theoretically derived eruption parameters and contemporary information allows a detailed analysis of eruption mechanisms.This eruption was the major stratospheric aerosol injection in the 1902–1903 period. However, mid- to low- latitude northern hemisphere temperature deviation data for the years following the eruption show no significant temperature decrease. This may be explained by the sulfur-poor nature of dacite magmas, suggesting that volatile composition, rather than mass of volatiles, is the controlling parameter in climatic response to explosive eruptions.     

Plume dynamics, thermal energy and long-distance transport of vulcanian eruption clouds from Augustine Volcano, Alaska

Augustine, an island volcano in Lower Cook Inlet, southern Alaska, erupted in January, 1976, after 12 years of dormancy. By April, when the eruptions ended, a new lava dome had been extruded into the summit crater and about 0.1 km³ of pyroclastics had been deposited on the island, mainly as pyroclastic debris avalanches and pumice flows. The ventclearing phase in January was highly explosive and we have been able to document 13 major vulcanian eruptions.The timing, thermal energy, mass loading of fine particles and the horizontal dispersion of these eruption clouds were determined from radar measurements of cloud height, reports of pilots flying in plumes, satellite photography, seismic records and infrasonic detection of air waves. A lower estimate of the mass of fine (r < 68 μm) particles injected into the troposphere from the 13 main eruptions in January is 5.5–18 × 10¹² g. The corresponding mass loading of fine particles within individual eruption clouds is 0.3–1 g m⁻³. We calculated thermal energies of 4 × 10¹⁴ to 35 × 10¹⁴ J for individual eruptions by applying convective plume rise theory to observed cloud heights and seismically determined eruption durations. This energy range compares favorably with the 4–16 × 10¹⁴ J of thermal energy, calculated from the cooling of juvenile material contained in a typical eruption cloud.The vulcanian eruption clouds stayed intact for at least 700 km downwind. Satellite images in both visible and infrared wavebands, showing the Gulf of Alaska just after sunrise on January 23, reveal a series of puffs strung out downwind from the volcano, 20–30 km in diameter and with their tops at altitudes of about 8 km, overlying a continuous plume at altitude 4 km. Each puff corresponded to a seismically and infrasonically timed eruption. A substantial portion of the material injected into the atmosphere between January 22 and 25 was rapidly transported by the subpolar jet stream through southwestern Canada and the western United States, then northeast across the States into the Atlantic. The clouds were observed passing over Tucson, Arizona, on January 25 at an elevation of 7 km.Several of the eruptions penetrated into the stratosphere. Sun photometer measurements, taken at Mauna Loa, Hawaii, six weeks after the eruption, showed an increased stratospheric optical thickness of 0.01 (wavelength 0.5 μm), which decayed in about 5 months. The maximum column mass loading of the veil was 4–10 × 10⁻⁷ g cm⁻². The mass of the veil, spread-ever a fourth of the earth's surface, is 10 to 100 times larger than can be accounted for by assuming that injected ash and converted sulfate particles from the 13 main Augustine eruptions are the only components contributing to the stratospheric turbidity observed at Mauna Loa.     

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.     

Explosive volcanism (VEI 6) without caldera formation: insight from Huaynaputina volcano, southern Peru

Through examination of the vent region of Volcán Huaynaputina, Peru, we address why some major explosive eruptions do not produce an equivalent caldera at the eruption site. Here, in 1600, more than 11 km³ DRE (VEI 6) were erupted in three stages without developing a volumetrically equivalent caldera. Fieldwork and analysis of aerial photographs reveal evidence for cryptic collapse in the form of two small subsidence structures. The first is a small non-coherent collapse that is superimposed on a cored-out vent. This structure is delimited by a partial ring of steep faults estimated at 0.85 by 0.95 km. Collapse was non-coherent with an inwardly tilted terrace in the north and a southern sector broken up along a pre-existing local fault. Displacement was variable along this fault, but subsidence of approximately 70 m was found and caused the formation of restricted extensional gashes in the periphery. The second subsidence structure developed at the margin of a dome; the structure has a diameter of 0.56 km and crosscuts the non-coherent collapse structure. Subsidence of the dome occurred along a series of up to seven concentric listric faults that together accommodate approximately 14 m of subsidence. Both subsidence structures total 0.043 km³ in volume, and are much smaller than the 11 km³ of erupted magma. Crosscutting relationships show that subsidence occurred during stages II and III when ∼2 km³ was erupted and not during the main plinian eruption of stage I (8.8 km³). The mismatch in erupted volume vs. subsidence volume is the result of a complex plumbing system. The stage I magma that constitutes the bulk of the erupted volume is thought to originate from a ∼20-km-deep regional reservoir based on petrological constraints supported by seismic data. The underpressure resulting from the extraction of a relatively small fraction of magma from the deep reservoir was not sufficient enough to trigger collapse at the surface, but the eruption left a 0.56-km diameter cored-out vent in which a dome was emplaced at the end of stage II. Petrologic evidence suggests that the stage I magma interacted with and remobilized a shallow crystal mush (∼4–6 km) that erupted during stage II and III. As the crystal mush erupted from the shallow reservoir, depressurization led to incremental subsidence of the non-coherent collapse structure. As the stage III eruption waned, local pressure release caused subsidence of the dome. Our findings highlight the importance of a connected magma reservoir, the complexity of the plumbing system, and the pattern of underpressure in controlling the nature of collapse during explosive eruptions. Huaynaputina shows that some major explosive eruptions are not always associated with caldera collapse. Editorial responsibility: J Stix     

The December 4, 1983 to February 18, 1984 eruption of Piton de la Fournaise (La Reunion, Indian Ocean): Description and interpretation

On December 4, 1983 an eruption started at vents located 1.5 km southwest of the summit of Piton de la Fournaise at the base of the central cone. After 31 months of quiescence this was one of the longest repose period in the last fifty years. The eruption had two phases: December 4 to January 18 and January 18 to February 18. Phase 1 produced about 8 × 10⁶ m³ of lava and Phase II about 9 × 10⁶ m³. The erupted lava is an aphyric basalt whose mineralogical and geochemical composition is close to that of other lavas emitted since 1977.The precursors of the December 4 outbreak were limited to two-week shallow (1.5–3 km) seismic crisis of fewer than 50 events. No long-term increase was noted in the local seismicity which is very quiet during repose periods and no long-term ground inflation preceded the eruption. Outbreaks of Phases I and II were preceded by short (2.5 hours and 1.5 hours) seismic swarms corresponding to the rise of magma toward the surface from a shallow reservoir. Large ground deformation explained by the emplacement of the shallow intrusions, was recorded during the seismic swarms. A summit inflation was observed in early January, before the phase II outbreak, while the phase I eruption was still continuing.Piton de la Fournaise volcanological observatory was installed in 1980. Seismic and ground deformation data now available for a period of 4 years including the 1981 and the 1983–1984 eruptions, allow us to describe the physical behavior of the volcano during this period. These observations lead us to propose that the magma transfer from deep levels to the shallow magma reservoir is not a continuous process but a periodic one and that the shallow magma reservoir was not resupplied before the 1981 and 1983–1984 eruptions. Considerations on the eruptive history and the composition of recent lavas indicate that the reservoir was refilled in 1977.     

The Middle Scoria sequence: A Holocene violent strombolian, subplinian and phreatomagmatic eruption of Okmok volcano, Alaska

The Middle Scoria deposit represents an explosive eruption of basaltic andesite magma (54 wt. % SiO₂) from Okmok volcano during mid-Holocene time. The pattern of dispersal and characteristics of the ejecta indicate that the eruption opened explosively, with ash textural evidence for a limited degree of phreatomagmatism. The second phase of the eruption produced thick vesicular scoria deposits with grain texture, size and dispersal characteristics that indicate it was violent strombolian to subplinian in style. The third eruptive phase produced deposits with a shift towards grain shapes that are dense, blocky, and poorly vesicular, and intermittent surge layers, indicating later transitions between magmatic (violent strombolian) to phreatomagmatic (vulcanian) eruptive styles. Isopach maps yield bulk volume estimates that range from 0.06 to 0.43 km³, with ~ 0.04 to 0.25 km³ total DRE. The associated column heights and mass discharge values calculated from isopleth maps of individual Middle Scoria layers are 8.5 – 14 km and 0.4 to 45 × 10⁶ kg/s. The Middle Scoria tephras are enriched in plagioclase microlites that have the textural characteristics of rapid magma ascent and relatively high degrees of effective undercooling. Those textures probably reflect the rapid magma ascent accompanying the violent strombolian and subplinian phases of the eruption. In the later stages of the eruption, the plagioclase microlite number densities decrease and textures include more tabular plagioclase, indicating a slowing of the ascent rate. The findings on the Middle Scoria are consistent with other explosive mafic eruptions, and show that outside of the two large caldera-forming eruptions, Okmok is also capable of producing violent mafic eruptions, marked by varying degrees of phreatomagmatism.     

Fundamental changes in the activity of the natrocarbonatite volcano Oldoinyo Lengai, Tanzania

On September 4, 2007, after 25 years of effusive natrocarbonatite eruptions, the eruptive activity of Oldoinyo Lengai (OL), N Tanzania, changed abruptly to episodic explosive eruptions. This transition was preceded by a voluminous lava eruption in March 2006, a year of quiescence, resumption of natrocarbonatite eruptions in June 2007, and a volcano-tectonic earthquake swarm in July 2007. Despite the lack of ground-based monitoring, the evolution in OL eruption dynamics is documented based on the available field observations, ASTER and MODIS satellite images, and almost-daily photos provided by local pilots. Satellite data enabled identification of a phase of voluminous lava effusion in the 2 weeks prior to the onset of explosive eruptions. After the onset, the activity varied from 100 m high ash jets to 2–15 km high violent, steady or unsteady, eruption columns dispersing ash to 100 km distance. The explosive eruptions built up a ∼400 m wide, ∼75 m high intra-crater pyroclastic cone. Time series data for eruption column height show distinct peaks at the end of September 2007 and February 2008, the latter being associated with the first pyroclastic flows to be documented at OL. Chemical analyses of the erupted products, presented in a companion paper (Keller et al. 2010), show that the 2007–2008 explosive eruptions are associated with an undersaturated carbonated silicate melt. This new phase of explosive eruptions provides constraints on the factors causing the transition from natrocarbonatite effusive eruptions to explosive eruptions of carbonated nephelinite magma, observed repetitively in the last 100 years at OL.     

Horizontal ground deformation patterns and magma storage during the Puu Oo eruption of Kilauea volcano,Hawaii: episodes 22–42

Horizontal ground deformation measurements were made repeatedly with an electronic distance meter near the Puu Oo eruption site approximately perpendicular to Kilauea's east rift zone (ERZ) before and after eruptive episodes 22–42. Line lengths gradually extended during repose periods and rapidly contracted about the same amount following eruptions. The repeated extension and contraction of the measured lines are best explained by the elastic response of the country rock to the addition and subsequent eruption of magma from a local reservoir. The deformation patterns are modeled to constrain the geometry and location of the local reservoir near Puu Oo. The observed deformation is consistent with deformation patterns that would be produced by the expansion of a shallow, steeply dipping dike just uprift of Puu Oo striking parallel to the trend of the ERZ. The modeled dike is centered about 800 m uprift of Puu Oo. Its top is at a depth of 0.4 km, its bottom at about 2.9 km, and the length is about 1.6 km; the dike strikes N65° E and dips at about 87°SE. The model indicates that the dike expanded by 11 cm during repose periods, for an average volumetric expansion of nearly 500 000 m³. The volume of magma added to the dike during repose periods was variable but correlates positively with the volume of erupted lava of the subsequent eruption and represents about 8% of the new lava extruded. Dike geometry and expansion values are used to estimate the pressure increase near the eruption site due to the accumulation of magma during repose periods. On average, vent pressures increased by about 0.38 MPa during the repose periods, one-third of the pressure increase at the summit. The model indicates that the dikelike body below Puu Oo grew in volume from 3 million cubic meters (Mm³) to about 10–12 Mm³ during the series of eruptions. The width of this body was probably about 2.5–3.0 m. No net long-term deformation was detected along the measured deformation lines.     

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     

The storage and release of magma on Mount Etna

The eruptive history of Etna during the past 450 years provides data on effusion rates, volumes of magma involved, and the nature of the eruptive conduits. These data are interpreted in terms of a two-part intravolcanic magma reservoir which feeds the flank eruptions through dike-like conduits. The structural framework of the volcano which controls the spatial distribution of eruptive sites is partly inherited from the basement and partly controlled by the central magma column and the surrounding caldera boundary faults. Hydraulic fracturing theory predicts that the central magma column will fail at depths below 1 km if the tensile strength of the conduit rocks is about 100 bars and that a peak fracturing capability will be reached between 1 and 2 km depth. This inference agrees well with the peak of flank eruptive activity at 1.4 km below the summit observed in the data on the loci of eruptions. The average flank-eruption feeding dike is defined and shown to be capable of the observed maximum effusion rates (20–100 m³ s⁻¹) from magmatic pressure differences of 30–150 bars     

Magma discharge variations during the 2011 eruptions of Shinmoe-dake volcano, Japan, revealed by geodetic and satellite observations

We present precise geodetic and satellite observation-based estimations of the erupted volume and discharge rate of magma during the 2011 eruptions of Kirishima-Shinmoe-dake volcano, Japan. During these events, the type and intensity of eruption drastically changed within a week, with three major sub-Plinian eruptions on January 26 and 27, and a continuous lava extrusion from January 29 to 31. In response to each eruptive event, borehole-type tiltmeters detected deflation of a magma chamber caused by migration of magma to the surface. These measurements enabled us to estimate the geodetic volume change in the magma chamber caused by each eruptive event. Erupted volumes and discharge rates were constrained during lava extrusion using synthetic aperture radar satellite imaging of lava accumulation inside the summit crater. Combining the geodetic volume change and the volume of lava extrusion enabled the determination of the erupted volume and discharge rate during each sub-Plinian event. These precise estimates provide important information about magma storage conditions in magma chambers and eruption column dynamics, and indicate that the Shinmoe-dake eruptions occurred in a critical state between explosive and effusive eruption.     

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.     

Volumetric characteristics of lava flows from interferometric radar and multispectral satellite data: the 1995 Fernandina and 1998 Cerro Azul eruptions in the western Galápagos

We have used a suite of remotely sensed data, numerical lava flow modeling, and field observations to determine quantitative characteristics of the 1995 Fernandina and 1998 Cerro Azul eruptions in the western Galápagos Islands. Flank lava flow areas, volumes, instantaneous effusion rates, and average effusion rates were all determined for these two eruptions, for which only limited syn-eruptive field observations are available. Using data from SPOT, TOPSAR, ERS-1, and ERS-2, we determined that the 1995 Fernandina flow covers a subaerial area of 6.5×10⁶ m² and has a subaerial dense rock equivalent (DRE) volume of 42×10⁶ m³. Field observations, ATSR satellite data, and the FLOWGO numerical model allow us to determine that the effusion rate declined exponentially from a high of ~60–200 m³ s^-1 during the first few hours to <5 m³ s^-1 prior to ceasing after 73 days, with a mean effusion rate of 4–16 m³ s^-1. Integrating the ATSR-derived, exponentially declining effusion rate over the eruption duration produces a total (subaerial + submarine) DRE volume of between 27 and 100×10⁶ m³, the range in values being due to differing assumptions about heat loss characteristics; only values in the higher part of this range are consistent with the independently derived subaerial volume. Using SPOT, TOPSAR, ERS-1, and ERS-2 data, we determine that the 1998 Cerro Azul flow is 16 km long, covers 16 km², and has a DRE volume of 54×10⁶ m³. FLOWGO produces at-vent velocity and effusion rate values of 11 m s^-1 and ~600 m³ s^-1, respectively. The velocity value agrees well with the 12 m s^-1 estimated in the field. The mean effusion rate (total DRE volume/duration) was 7–47 m³ s^-1. Dike dimensions, fissure lengths, and pressure gradients along the conduit based on magma chamber depth estimates of 3–5 km produce mean effusion rates for the two eruptions that range over nearly four orders of magnitude, the range being due to uncertainty in the magma viscosity, dike dimensions, and pressure gradient between magma chamber and vent. Although somewhat consistent with mean effusion rates from other techniques, their wide range makes them less useful. The exponentially declining effusion rates during both eruptions are consistent with release of elastic strain being the driving mechanism of the eruptions. Our results provide independent input parameters for previously published theoretical relationships between magma chamber pressurization and eruption rates that constrain chamber volumes and increases in volume prior to eruption, as well as time constants of exponential decay during the eruption. The results and theoretical relationships combine to indicate that at both volcanoes probably 25–30% of the volumetric increase in the magma chamber erupted as lava onto the surface. In both eruptions the lava flow volumes are less than 1% of the magma chamber volume.     

Volcanic hazards from Bezymianny- and Bandai-type eruptions

Major slope failures are a significant degradational process at volcanoes. Slope failures and associated explosive eruptions have resulted in more than 20 000 fatalities in the past 400 years; the historic record provides evidence for at least six of these events in the past century. Several historic debris avalanches exceed 1 km³ in volume. Holocene avalanches an order of magnitude larger have traveled 50–100 km from the source volcano and affected areas of 500–1500 km². Historic eruptions associated with major slope failures include those with a magmatic component (Bezymianny type) and those solely phreatic (Bandai type). The associated gravitational failures remove major segments of the volcanoes, creating massive horseshoe-shaped depressions commonly of caldera size. The paroxysmal phase of a Bezymianny-type eruption may include powerful lateral explosions and pumiceous pyroclastic flows; it is often followed by construction of lava dome or pyroclastic cone in the new crater. Bandai-type eruptions begin and end with the paroxysmal phase, during which slope failure removes a portion of the edifice. Massive volcanic landslides can also occur without related explosive eruptions, as at the Unzen volcano in 1792.The main potential hazards from these events derive from lateral blasts, the debris avalanche itself, and avalanche-induced tsunamis. Lateral blasts produced by sudden decompression of hydrothermal and/or magmatic systems can devastate areas in excess of 500km² at velocities exceeding 100 m s^–1. The ratio of area covered to distance traveled for the Mount St. Helens and Bezymianny lateral blasts exceeds that of many pyroclastic flows or surges of comparable volume. The potential for large-scale lateral blasts is likely related to the location of magma at the time of slope failure and appears highest when magma has intruded into the upper edifice, as at Mount St. Helens and Bezymianny.Debris avalanches can move faster than 100 ms^–1 and travel tens of kilometers. When not confined by valley walls, avalanches can affect wide areas beyond the volcano's flanks. Tsunamis from debris avalanches at coastal volcanoes have caused more fatalities than have the landslides themselves or associated eruptions. The probable travel distance (L) of avalanches can be estimated by considering the potential vertical drop (H). Data from a catalog of around 200 debris avalanches indicates that the H/L rations for avalanches with volumes of 0.1–1 km³ average 0.13 and range 0.09–0.18; for avalanches exceeding 1 km³, H/L ratios average 0.09 and range 0.5–0.13.Large-scale deformation of the volcanic edefice and intense local seismicity precede many slope failures and can indicate the likely failure direction and orientation of potential lateral blasts. The nature and duration of precursory activity vary widely, and the timing of slope faliure greatly affects the type of associated eruption. Bandai-type eruptions are particularly difficult to anticipate because they typically climax suddenly without precursory eruptions and may be preceded by only short periods of seismicity.     

Dynamics of ash-dominated eruptions at Vesuvius: the post-512 AD AS1a event

Recent stratigraphic studies at Vesuvius have revealed that, during the past 4,000 years, long lasting, moderate to low-intensity eruptions, associated with continuous or pulsating ash emission, have repeatedly occurred. The present work focuses on the AS1a eruption, the first of a series of ash-dominated explosive episodes which characterized the period between the two Subplinian eruptions of 472 AD and 1631 AD. The deposits of this eruption consist of an alternation of massive and thinly laminated ash layers and minor well sorted lapilli beds, reflecting the pulsatory injection into the atmosphere of variably concentrated ash-plumes alternating with Violent Strombolian stages. Despite its nearly constant chemical composition, the juvenile material shows variable external clast morphologies and groundmass textures, reflecting the fragmentation of a magma body with lateral and/or vertical gradients in both vesicularity and crystal content. Glass compositions and mineralogical assemblages indicate that the eruption was fed by rather homogeneous phonotephritic magma batches rising from a reservoir located at ~ 4 km (100 MPa) depth, with fluctuations between magma delivery and magma discharge. Using crystal size distribution (CSD) analyses of plagioclase and leucite microlites, we estimate that the transit time of the magma in the conduit was on the order of ~ 2 days, corresponding to an ascent rate of around 2 × 10⁻² ms⁻¹. Accordingly, assuming a typical conduit diameter for this type of eruption, the minimum duration of the AS1a event is between about 1.5 and 6 years. Magma fragmentation occurred in an inertially driven regime that, in a magma with low viscosity and surface tension, can act also under conditions of slow ascent.