Sakurajima volcano: a physico-chemical study of the health consequences of long-term exposure to volcanic ash

Regular eruptions from Sakurajima volcano, Japan, repeatedly cover local urban areas with volcanic ash. The frequency of exposure of local populations to the ash led to substantial concerns about possible respiratory health hazards, resulting in many epidemiological and toxicological studies being carried out in the 1980s. However, very few mineralogical data were available for determination of whether the ash was sufficiently fine to present a respiratory hazard. In this study, we review the existing studies and carry out mineralogical, geochemical and toxicological analyses to address whether the ash from Sakurajima has the potential to cause respiratory health problems. The results show that the amount of respirable (<4 μm) material produced by the volcano is highly variable in different eruptions (1.1–18.8 vol.%). The finest samples derive from historical, plinian eruptions but considerable amounts of respirable material were also produced from the most recent vulcanian eruptive phase (since 1955). The amount of cristobalite, a crystalline silica polymorph which has the potential to cause chronic respiratory diseases, is ~3–5 wt.% in the bulk ash. Scanning electron microscope and transmission electron microscope imaging showed no fibrous particles similar to asbestos particles. Surface reactivity tests showed that the ash did not produce significant amounts of highly reactive hydroxyl radicals (0.09–1.35 μmol m⁻² at 30 min.) in comparison to other volcanic ash types. A basic toxicology assay to assess the ability of ash to rupture the membrane of red blood cells showed low propensity for haemolysis. The findings suggest that the potential health hazard of the ash is low, but exposure and respiratory conditions should still be monitored given the high frequency and durations of exposure.     

Mineralogical analyses and in vitro screening tests for the rapid evaluation of the health hazard of volcanic ash at Rabaul volcano, Papua New Guinea

The continuous ash and gas emissions from the Tavurvur cone in Rabaul caldera, Papua New Guinea, during 2007–08, raised concerns regarding how exposure would affect the respiratory health of nearby populations and impact on the environment. As part of a formal evaluation of the effects of volcanic emissions on public health, we investigated the potential health hazard of the ash using a suite of selected mineralogical analyses and in vitro toxicity screening tests. The trachy-andesitic ash comprised 2.1–6.7 vol.% respirable (sub-4 μm diameter) particles. The crystalline silica content was 1.9–5.0 wt.% cristobalite (in the bulk sample) with trace amounts of quartz and/or tridymite. Scanning electron microscopy showed that the ash particles were angular with sparse, fibre-like particles (∼3–60 μm max. diameter) observed in some samples, which we confirmed to be CaSO₄ (gypsum, at <6 wt.% in the bulk samples) and not asbestiform fibres. The ash specific surface area was low (0.1–2.7 m² g⁻¹). The leached solution from one of the ash samples was slightly acidic (pH 5.6), but did not contain high levels of toxic metals (such as F, Cu, Zn, Mn, As, Ni and Cd) when compared to previously tested volcanic ash leachates. Ash samples generated potentially-harmful hydroxyl radicals through an iron-mediated catalytic reaction, in the range of 0.15–2.47 μmol m⁻² (after 30 min of reaction). However, measurement of particle oxidative capacity (potential oxidative stress reaction using ascorbic acid) and silica-like injury to red blood cells (erythrolysis assay, i.e. measurement of cell death) nevertheless revealed low biological reactivity. The findings suggest that acute exposure to the ash would have a limited potential to exacerbate pre-existing conditions such as asthma or chronic bronchitis, and the potential for chronic exposure leading to silicosis was low.     

Formation of high-grade ignimbrites Part II. A pyroclastic suspension current model with implications also for low-grade ignimbrites

 Analogue experiments in part I led to the conclusion that pyroclastic flows depositing very high-grade ignimbrite move as dilute suspension currents. In the thermo–fluid–dynamical model developed, the degree of cooling of expanded turbulent pyroclastic flows dynamically evolves in response to entrainment of air and mass loss to sedimentation. Initial conditions of the currents are derived from column-collapse modeling for magmas with an initial H₂O content of 1–3 wt.% erupting through circular vents and caldera ring-fissures. The flows spread either longitudinally or radially from source up to a runout distance that increases with higher mass flux but decreases with higher gas content, temperature, bottom slope and coarser initial grain size. Progressive dilution by entrainment and sedimentation causes pyroclastic currents to transform into buoyant ash plumes at the runout distance. The ash plumes reach stratospheric heights and distribute 30–80% of the erupted material as widespread co-ignimbrite ash. Pyroclastic suspension currents with initial mass fluxes of 10⁷-10¹² kg/s can spread for tens of kilometers with only limited cooling, although they move as supercritical, strongly entraining currents for the eruption conditions considered here. With increasing eruption mass flux, cooling during passage through the fountain diminishes while cooling during flow transport increases. The net effect is that eruption temperature exerts the prime control on emplacement temperature. Pyroclastic suspension currents can form welded ignimbrite across their entire extent if eruption temperature is T_o>1.3^.T_mw, the minimum welding temperature. High eruption rates, a large fraction of fine ash, and a ring-fissure vent favor the formation of extensive high-grade ignimbrite. For very hot eruptions producing sticky, partially molten pyroclasts, analysis of particle aggregation systematics shows that factors favoring longer runout also favor more efficient aggregation, which reduces runout. As a result, very high-grade ignimbrites cannot spread more than a few tens of kilometers from their source. In cooler pyroclastic currents, particles do not aggregate, and the sedimentation process may involve re-entrainment of particles, which potentially leads to more extensive cooling and longer runout; such effects, however, are only significant when net erosion of substrate occurs. Model results can be employed to estimate mass flux and duration of ignimbrite eruptions from measured ignimbrite masses and aspect ratios. The model also provides an alternative explanation of the observed decrease in H/Lratios with ignimbrite mass. Received: 10 May 1998 / Accepted: 21 October 1998     

Characterization of respirable volcanic ash from the Soufrière Hills volcano,Montserrat, with implications for human health hazards

Volcanic ash, generated in the long-lived eruption of the Soufrière Hills volcano, Montserrat, is shown to contain respirable (sub-4 μm) particles and cristobalite, a crystalline silica polymorph. Respirable particles of cristobalite can cause silicosis, raising the possibility that volcanic ash is a respiratory health hazard. This study considers some of the main factors which affect human exposure to respirable volcanic ash, namely, the composition and proportions of respirable ash, and the composition and concentrations of airborne suspended particulates. The composition, size distribution and proportion (by weight) of respirable particles in representative samples of the Soufrière Hills tephra (dome-collapse ash-fall deposits, dome-collapse pyroclastic-flow matrix, Vulcanian explosion ash and mixed ash) have been characterized. Dome-collapse ash-fall deposits are significantly richer in respirable particles (12 wt%) than the other tephra samples, in particular the matrices of dome-collapse pyroclastic-flow deposits (3 wt%). Within the respirable fraction, dome-collapse ash contains the highest proportion of crystalline silica particles (20–27 number%, of which 97 wt% is cristobalite), compared with other primary tephra types (0.4–5.6 number%). This enrichment of crystalline silica in the dome-collapse ash is most pronounced in the very fine particle fraction (sub-2 μm). The results are explained as being due to significant size fractionation during fragmentation of pyroclastic flows, resulting in a fines-depleted dome-collapse matrix and a fines-rich dome-collapse ash deposit. For all sample types, the sub-4 μm fraction comprises 45–55 wt% of the sub-10 μm fraction. Aeolian deposit, lahar deposit and airborne samples of suspended ash, collected on filters, were characterized. These samples show enrichment of crystalline silica in the respirable fraction (10–18 number%). The results are consistent with ash in the environment having a mixed origin but originating predominantly from dome-collapse eruptions. The reworked ash, however, contains low proportions of respirable ash (～3 wt%) compared to primary ash samples. The concentration of ash particles re-suspended by road vehicles on Montserrat is found to decrease exponentially with height above the ground, indicating higher exposure for children compared with adults: PM₄ concentration at 0.9 m (height of two-year-old child) is 3 times that at 1.8 m (adult height). The composition of the re-suspended road particles is similar to that re-suspended by the wind.     

Tephra stratigraphy and eruptive volume of the May, 2008, Chaitén eruption,Chile

On May 1st 2008 Mount Chaitén (southern Chile) interrupted a long period of quiescence, generating a sequence of explosive eruptions and causing the evacuation of Chaitén town located a few kilometers south of the volcano. The activity was characterized by several explosive events each associated with plumes which reached up to about 19 km above sea level. The products were dispersed across a wide area, with the finest ash reaching the Atlantic coast of Argentina. Our field observations in the proximal-medial area (3–25 km from the vent) indicate that the May 2008 tephra deposit consists of numerous layers, most of which can be correlated with individual eruptive events. These layers vary from extremely fine-grained ash to layers of lapilli and blocks, composed of both juvenile and lithic material. Here we describe the stratigraphy and physical characteristics of the May 2008 deposits, and propose a reconstruction of the timing of the May 2008 events. The deposits are mainly associated with the three main explosive phases which occurred on 1st–2nd May, 3rd–5th May and 6th May, with an estimated bulk tephra volume of 0.5–1.0 km³ (integration of both exponential and power-law fitting). For the 6th May event, represented by a layer composed mainly of lithic lapilli and blocks (>2 mm), an isopleth map was compiled from which a 19 km plume height was determined, which is in good agreement with satellite observations.     

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.     

Strain-induced magma degassing: insights from simple-shear experiments on bubble bearing melts

Experiments have been performed to determine the effect of deformation on degassing of bubble-bearing melts. Cylindrical specimens of phonolitic composition, initial water content of 1.5 wt.% and 2 vol.% bubbles, have been deformed in simple-shear (torsional configuration) in an internally heated Paterson-type pressure vessel at temperatures of 798–848 K, 100–180 MPa confining pressure and different final strains. Micro-structural analyses of the samples before and after deformation have been performed in two and three dimensions using optical microscopy, a nanotomography machine and synchrotron tomography. The water content of the glasses before and after deformation has been measured using Fourier Transform Infrared Spectroscopy (FTIR). In samples strained up to a total of γ ∼ 2 the bubbles record accurately the total strain, whereas at higher strains (γ ∼ 10) the bubbles become very flattened and elongate in the direction of shear. The residual water content of the glasses remains constant up to a strain of γ ∼ 2 and then decreases to about 0.2 wt.% at γ ∼ 10. Results show that strain enhances bubble coalescence and degassing even at low bubble volume-fractions. Noticeably, deformation produced a strongly water under-saturated melt. This suggests that degassing may occur at great depths in the volcanic conduit and may force the magma to become super-cooled early during ascent to the Earth’s surface potentially contributing to the genesis of obsidian.     

Eruptive style of the young high-Mg basaltic-andesite Pelagatos scoria cone,southeast of México City

The eruption of the Pelagatos scoria cone in the Sierra Chichinautzin monogenetic field near the southern suburbs of Mexico City occurred less than 14,000 years ago. The eruption initiated at a fissure with an effusive phase that formed a 7-km-long lava flow, and continued with a phase of alternating and/or simultaneous explosive and effusive activity that built a 50-m-high scoria cone on the western end of the fissure and formed a compound lava flow-field near the vent. The eruption ended with the emplacement of a short lava flow that breached the cone and was accompanied by weak explosions at the crater. Products consist of a microlite-rich high-Mg basaltic andesite. Samples were analyzed to determine the magma’s initial properties as well as the effects of degassing-induced crystallization on eruptive style. Although distal ash fallout deposits from this eruption are not preserved, a recent quarry exposes a large section of the scoria cone. Detailed study of exposed layers allows us to elucidate the mode of cone-building activity. Petrological and textural data, combined with models calibrated by experimental work and melt-inclusion analyses of similar magmas elsewhere, indicate that the magma was initially hot (>1,200°C), gas-rich (up to 5 wt.% H₂O), crystal-poor (~10 vol.% Fo₉₀ olivine phenocrysts) and thus poorly viscous (40–80 Pa s). During the early phase, low magma ascent velocity at the fissure vent allowed low-viscosity magma to degas and crystallize during ascent, producing lava flows with elevated crystal contents at T < 1,100°C, and blocky surfaces. Later, the closure of the fissure by cooling dikes focused the magma flow at a narrow section of the fissure. This led to an increased magma ascent velocity. Rapid and shallow degassing (<3 km deep) triggered ~40 vol.% microlite crystallization. Limited times for gas-escape and higher magma viscosity (6 × 10⁵–4 × 10⁶ Pa s) drove strong explosions of highly (60–80 vol.%) and finely vesicular magma. Coarse clasts broke on landing, which implies brittle behavior due to complete solidification. This requires sufficient time to cool and in turn implies ejection heights of over 1 km, which is much higher than “normal” Strombolian activity. Hence, magma viscosity significantly impacts eruption style at monogenetic volcanoes because it affects the kinetics of shallow degassing. The long-lasting eruptions of Jorullo and Paricutin, which produced similar magmas in western México, were more explosive. This can be related to higher magma fluxes and total erupted volumes. Implications of this study are important because basaltic andesites are commonly erupted to form monogenetic scoria cones of the Trans-Mexican Volcanic Belt.     

Long term exposure to respirable volcanic ash on Montserrat: a time series simulation

Frequent ash fallout from long-lived eruptions (with active phases greater than 5 years) may lead to local populations experiencing unacceptably high cumulative exposures to respirable particulate matter. Ash from Montserrat has been shown to contain significant levels of cristobalite and other reactive agents that are associated with an increased risk of developing pneumoconiosis (including silicosis) and other long-term health problems. There are a number of difficulties associated with estimating risks in populations due to uncertain and wide ranging individual exposures, change in behaviour with time and the natural variation in individual response. Present estimates of risk in workers and other population groups are simplifications based on a limited number of exposure measurements taken on Montserrat (1996–1999), and exposure−response curves from epidemiological studies of coal workers exposed to siliceous dust. In this paper we present a method for calculating the long-term cumulative exposure to cristobalite from volcanic ash by Monte Carlo simulation. Code has been written to generate synthetic time series for volcanic activity, rainfall, ash deposition and erosion to give daily ash deposit values and cristobalite fraction at a range of locations. The daily mean personal exposure for PM₁₀ and cristobalite is obtained by sampling from a probability distribution, with distribution parameters dependent on occupation, ground deposit depth and daily weather conditions. Output from multiple runs is processed to calculate the exceedance probability for cumulative exposure over a range of occupation types, locations and exposure periods. Results are interpreted in terms of current occupational standards, and epidemiological exposure−response functions for silicosis are applied to quantify the long-term health risk. Assuming continuing volcanic activity, median risk of silicosis (profusion 1/0 or higher) for an average adult after 20 years continuous exposure is estimated to be approximately 0.5% in northern Montserrat to 1.6% in Cork Hill. The occupational group with the highest exposure to ash are gardeners, with a corresponding 2% to 4% risk of silicosis. In situations where opportunities for in-depth exposure studies are limited, computer simulations provide a good indication of risk based on current expert knowledge. By running the code for a range of input scenarios, the cost-benefit of mitigation measures (such as a programme of active ash clearance) can be estimated. Results also may be used to identify situations where full exposure studies or fieldwork would be beneficial. Editorial responsibility: J Stix     

Erratum to Compositional variations and magma mixing in the 1991 eruptions of Hudson volcano,Chile

The August 1991 eruptions of Hudson volcano produced ~2.7 km³ (dense rock equivalent, DRE) of basaltic to trachyandesitic pyroclastic deposits, making it one of the largest historical eruptions in South America. Phase 1 of the eruption (P1, April 8) involved both lava flows and a phreatomagmatic eruption from a fissure located in the NW corner of the caldera. The paroxysmal phase (P2) began several days later (April 12) with a Plinian-style eruption from a different vent 4 km to the south-southeast. Tephra from the 1991 eruption ranges in composition from basalt (phase 1) to trachyandesite (phase 2), with a distinct gap between the two erupted phases from 54–60 wt% SiO₂. A trend of decreasing SiO₂ is evident from the earliest part of the phase 2 eruption (unit A, 63–65 wt% SiO₂) to the end (unit D, 60–63 wt% SiO₂). Melt inclusion data and textures suggest that mixing occurred in magmas from both eruptive phases. The basaltic and trachyandesitic magmas can be genetically related through both magma mixing and fractional crystallization processes. A combination of observed phase assemblages, inferred water content, crystallinity, and geothermometry estimates suggest pre-eruptive storage of the phase 2 trachyandesite at pressures between ~50–100 megapascal (MPa) at 972 ± 26°C under water-saturated conditions (log fO₂ –10.33 (±0.2)). It is proposed that rising P1 basaltic magma intersected the lower part of the P2 magma storage region between 2 and 3 km depth. Subsequent mixing between the two magmas preferentially hybridized the lower part of the chamber. Basaltic magma continued advancing towards the surface as a dyke to eventually be erupted in the northwestern part of the Hudson caldera. The presence of tachylite in the P1 products suggests that some of the magma was stalled close to the surface (<0.5 km) prior to eruption. Seismicity related to magma movement and the P1 eruption, combined with chamber overpressure associated with basalt injection, may have created a pathway to the surface for the trachyandesite magma and subsequent P2 eruption at a different vent 4 km to the south-southeast. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Mafic Plinian volcanism and ignimbrite emplacement at Tofua volcano,Tonga

Tofua Island is the largest emergent mafic volcano within the Tofua arc, Tonga, southwest Pacific. The volcano is dominated by a distinctive caldera averaging 4 km in diameter, containing a freshwater lake in the south and east. The latest paroxysmal (VEI 5–6) explosive volcanism includes two phases of activity, each emplacing a high-grade ignimbrite. The products are basaltic andesites with between 52 wt.% and 57 wt.% SiO₂. The first and largest eruption caused the inward collapse of a stratovolcano and produced the ‘Tofua’ ignimbrite and a sub-circular caldera located slightly northwest of the island’s centre. This ignimbrite was deposited in a radial fashion over the entire island, with associated Plinian fall deposits up to 0.5 m thick on islands >40 km away. Common sub-rounded and frequently cauliform scoria bombs throughout the ignimbrite attest to a small degree of marginal magma–water interaction. The common intense welding of the coarse-grained eruptive products, however, suggests that the majority of the erupted magma was hot, water-undersaturated and supplied at high rates with moderately low fragmentation efficiency and low levels of interaction with external water. We propose that the development of a water-saturated dacite body at shallow (<6 km) depth resulted in failure of the chamber roof to cause sudden evacuation of material, producing a Plinian eruption column. Following a brief period of quiescence, large-scale faulting in the southeast of the island produced a second explosive phase believed to result from recharge of a chemically distinct magma depleted in incompatible elements. This similar, but smaller eruption, emplaced the ‘Hokula’ Ignimbrite sheet in the northeast of the island. A maximum total volume of 8 km³ of juvenile material was erupted by these events. The main eruption column is estimated to have reached a height of ∼12 km, and to have produced a major atmospheric injection of gas, and tephra recorded in the widespread series of fall deposits found on coral islands 40–80 km to the east (in the direction of regional upper-tropospheric winds). Radiocarbon dating of charcoal below the Tofua ignimbrite and organic material below the related fall units imply this eruption sequence occurred post 1,000 years BP. We estimate an eruption magnitude of 2.24 × 10¹³ kg, sulphur release of 12 Tg and tentatively assign this eruption to the AD 1030 volcanic sulphate spike recorded in Antarctic ice sheet records.     

Multiple origins of obsidian pyroclasts and implications for changes in the dynamics of the 1300 B.P. eruption of Newberry Volcano,USA

The pyroclastic deposits of the 1300 B.P. eruption of Newberry Volcano, OR, USA, contain minor amounts of obsidian (1–6 wt.%). The volatile (H₂O and CO₂) contents and textures of these clasts vary considerably. FTIR measurements of H₂O in obsidian pyroclasts range from 0.1 to 1.5 wt.% indicating equilibration pressures ≤20 MPa. CO₂ contents are low (<10 ppm) except in clasts that also contain xenolith powder that provided a local CO₂ source. Obsidian clasts exhibit a range of color and textural types that changed in relative proportion as the eruption progressed. Together these data indicate that there were multiple origins of obsidian and that the dominant source changed during the eruption. Early in the eruption, obsidian was almost entirely black or grey (microlite-bearing) and probably derived from dikes or wall rock fractures filled with vanguard magma or tuffisite that, together with wall rocks, were eroded and incorporated into the eruption column as the vent widened. Later in the eruption, following a brief cessation of activity, the proportion of obsidian to wallrock lithic clasts increased and new types of obsidian dominated, types that represent remnants of a shallow conduit plug, welded fallback material from within the conduit, and sheared and degassed magma from near the conduit walls. Analysis of bubble shapes preserved within obsidian indicates that shear stresses and shear rates varied by over two orders of magnitude, with maxima of 88 kPa and 10^−2.3 s⁻¹, respectively, based on an assumed magma temperature of 850°C. Furthermore, the highest shear rates and stresses, and the shortest flow times (10–20 min), are preserved in clasts that also contain wall rock. The longest deformation times (5 and 8 h) correspond to two microlite-rich clasts, suggesting that the higher microlite content results from slower ascent rates and/or longer magma residence times at shallow levels. Differences between obsidian pyroclasts from the Newberry eruption and those of the Mono Craters may relate to the nature of the conduit feeding the two events. From this comparison, we conclude that obsidian can provide information on time scales and mechanisms of pre-fragmentation magma ascent.     

Eruptive history and magmatic evolution of the 1.9 kyr Plinian dacitic Chiltepe Tephra from Apoyeque volcano in west-central Nicaragua

The youngest dacitic Plinian eruption in west-central Nicaragua, forming the 18 km³ Chiltepe Tephra (CT), occurred about nineteen hundred years ago at Apoyeque stratovolcano, which dominates the Chiltepe volcanic complex 15 km north of the capital Managua, where the CT is 2 m thick. We have traced the CT from its proximal facies at the crater rim, through the medial facies in the lowlands around Apoyeque, and to the distal facies up to 550 km offshore in the Pacific. While medial and distal facies consist of widespread Plinian fall deposits, the proximal facies reveals the complexity of this eruption, which we divide into four phases (I–IV). Interaction of rising magma with a pre-existing crater lake generated the phreatomagmatic opening phase I of the eruption, which produced ash fall with accretionary lapilli. Phase II marked a rapid change to persistent magmatic activity that yielded several large Plinian eruptions, declining through a period of unstable eruption conditions, followed by a short hiatus. Phase III began with unstable conditions, probably as a result of eastward migration and widening of the vent, leading to a second period of Plinian eruptions with three major events reaching magma discharge rates five times larger than those of phase II. Phase III again declined through unstable eruption conditions before magmatic activity terminated. Numerous explosions in the shallow hydrothermal system during the final phase IV resulted in the formation of a phreatic tuff ring on the rim of Apoyeque crater. The white, highly-vesicular, dacitic CT pumice contains plagioclase (An_45–68), orthopyroxene, clinopyroxene, and minor hornblende, apatite and titanomagnetite phenocrysts. A very subordinate fraction of gray pumice has the highest crystal content, the least evolved bulk-rock, but the most evolved matrix-glass composition. The CT dacite has two unusual compositional features: (1) all white dacite has the same melt (matrix-glass) composition such that variations in bulk-rock compositions (64–68 wt% SiO₂) simply reflect different phenocryst contents of 10–35%, interpreted as the result of gradual phenocryst settling in the magma chamber. (2) Abundant olivine crystals with a bimodal distribution in Mg# (modes at Mg# = 0.75 and Mg# = 0.8) are dispersed throughout the erupted dacite. These are clearly out of equilibrium with the dacitic melt and are interpreted as xenocrysts derived from the basaltic Nejapa-Miraflores volcanic lineament that intersects the Chiltepe volcanic complex and was contemporaneously active. Thermobarometric estimates place the dacitic CT magma reservoir in the upper crust (<250 MPa), with a temperature of about 890°C and about 5 wt% water dissolved in the melt. Comparing water and chlorine contents with respective solubility models suggests that volatile degassing began in the magma reservoir and triggered the CT eruption. From the vertical compositional variation pattern of the CT we deduce that the conduit tapped the magma chamber not at the top but from the side, at some deeper level, and that subsequent magma withdrawal was governed by both variations in discharge rate and possible upward migration and/or widening of the conduit entrance.     

The caldera-forming eruption of Volcán Ceboruco, Mexico

3 of magma erupted, ∼95% of which was deposited as fall layers. During most of the deposition of P1, eruptive intensity (mass flux) was almost constant at 4–8×10⁷ kg s⁻¹, producing a Plinian column 25–30 km in height. Size grading at the top of P1 indicates, however, that mass flux waned dramatically, and possibly that there was a brief pause in the eruption. During the post-P1 phase of the eruption, a much smaller volume of magma erupted, although mass flux varied by more than an order of magnitude. We suggest that caldera collapse began at the end of the P1 phase of the eruption, because along with the large differences in mass flux behavior between P1 and post-P1 layers, there were also dramatic changes in lithic content (P1 contains ∼8% lithics; post-P1 layers contain 30–60%) and magma composition (P1 is 98% rhyodacite; post-P1 layers are 60–90% rhyodacite). However, the total volume of magma erupted during the Jala pumice event is close to that estimated for the caldera. These observations appear to conflict with models which envision that, after an eruption is initiated by overpressure in the magma chamber, caldera collapse begins when the reservoir becomes underpressurized as a result of the removal of magma. The conflict arises because firstly, the P1 layer makes up too large a proportion (∼75%) of the total volume erupted to correspond to an overpressurized phase, and secondly, the caldera volume exceeds the post-P1 volume of magma by at least a factor of three. The mismatches between model and observations could be reconciled if collapse began near the beginning of the eruption, but no record of such early collapse is evident in the tephra sequence. The apparent inability to place the Jala pumice eruptive sequence into existing models of caldera collapse, which were constructed to explain the formation of calderas much greater in volume than that at Ceboruco, may indicate that differences in caldera mechanics exist that depend on size or that a more general model for caldera formation is needed. Received: 18 November 1998 / Accepted: 23 October 1999     

Characterization of shape and terminal velocity of tephra particles erupted during the 2002 eruption of Etna volcano,Italy

In this paper, we present a complete morphological characterization of the ash particles erupted on 18 December 2002 from Etna volcano, Italy. The work is based on the acquisition and processing of bidimensional digital images carried out by scanning electron microscopy (SEM) to obtain shape parameters by image analysis. We measure aspect ratio (AR), form factor (FF), compactness (CC), and rectangularity (RT) of 2065 ash particles with size between 0.026 and 1.122 mm. We evaluate the variation of these parameters as a function of the grain-size. Ash particles with a diameter of <0.125 mm vary from mostly equant to very equant, ash particles between 0.125 and 0.250 mm have an intermediate shape, and particles with diameters >0.250 mm are subelongate. We find that, on average, particles with a diameter of <0.250 mm are subrounded, particles between 0.250 and 0.50 mm are subangular, and particles >0.50 mm are angular. Using this morphological analysis and an empirical relation between the drag coefficient (C _D) and Reynolds number (R _e) of Wilson and Huang (Earth Planet Sci Lett 44:311–324, 1979), we calculate the terminal settling velocities (V _WH). The comparisons between these velocities and those calculated with the formula of Kunii and Levenspiel (Fluidization engineering. Wiley, New York, pp 97, 1969) (V _KL), which considers ash particles as spheres, show that V _KL are in average 1.28 greater than V _WH. Hence, we quantify the systematic error on the spatial distribution of the mass computed around the volcano carried out by tephra dispersal models when the particles are assumed to be spherical.     

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.     

Petrology of lavas from the Puu Oo eruption of Kilauea Volcano: III. The Kupaianaha episode (1986–1992)

 The Puu Oo eruption has been remarkable in the historical record of Kilauea Volcano for its duration (over 13 years), volume (>1 km³) and compositional variation (5.7–10 wt.% MgO). During the summer of 1986, the main vent for lava production moved 3 km down the east rift zone and the eruption style changed from episodic geyser-like fountaining at Puu Oo to virtually continuous, relatively quiescent effusion at the Kupaianaha vent. This paper examines this next chapter in the Puu Oo eruption, episodes 48 and 49, and presents new ICP-MS trace element and Pb-, Sr-, and Nd-isotope data for the entire eruption (1983–1994). Nearly aphyric to weakly olivine-phyric lavas were erupted during episodes 48 and 49. The variation in MgO content of Kupaianaha lavas erupted before 1990 correlates with changes in tilt at the summit of Kilauea, both of which probably were controlled by variations in Kilauea's magma supply rate. These lavas contain euhedral olivines which generally are in equilibrium with whole-rock compositions, although some of the more mafic lavas which erupted during 1990, a period of frequent pauses in the eruption, accumulated 2–4 vol.% olivine. The highest forsterite content of olivines (∼85%) in Kupaianaha lavas indicates that the parental magmas for these lavas had MgO contents of ∼10 wt.%, which equals the highest observed value for lavas during this eruption. The composition of the Puu Oo lavas has progressively changed during the eruption. Since early 1985 (episode 30), when mixing between an evolved rift zone magma and a more mafic summit reservoir-derived magma ended, the normalized (to 10 wt.% MgO) abundances of highly incompatible elements and CaO have systematically decreased with time, whereas ratios of these trace elements and Pb, Sr, and Nd isotopes, and the abundances of Y and Yb, have remained relatively unchanged. These results indicate that the Hawaiian plume source for Puu Oo magmas must be relatively homogeneous on a scale of 10–20 km³ (assuming 5–10% partial melting), and that localized melting within the plume has apparently progressively depleted its incompatible elements and clinopyroxene component as the eruption continued. The rate of variation of highly incompatible elements in Puu Oo lavas is much greater than that observed for Kilauea historical summit lavas (e.g., Ba/Y 0.09 a^–1 vs ∼0.03 a^–1). This rapid change indicates that Puu Oo magmas did not mix thoroughly with magma in the summit reservoir. Thus, except for variable amounts of olivine fractionation, the geochemical variation in these lavas is predominantly controlled by mantle processes. Received: 8 March 1996 / Accepted: 30 April 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     

Simulation of the 1980 eruption of Mount St. Helens using the ash-tracking model PUFF

The dispersal of volcanic ash from the May 18, 1980 eruption of Mount St. Helens (MSH) has been simulated using the Lagrangian ash-tracking model PUFF. Previous applications of the model were limited to smaller, short-lived eruptions with ash dispersal occurring mainly within the troposphere. Two high-resolution atmospheric reanalysis datasets (ERA-40 and NCEP/NCAR-40) allowed MSH ash cloud dispersal to be simulated up to 30 km elevation. The 1980 eruption was divided into two distinct eruptive phases, (1) an initial, relatively short-lived blast/surge phase that injected ash up to 30 km and (2) a subsequent nine-hour plinian phase that maintained an average eruption column height of 16 km. Using PUFF, the two phases of the MSH eruption were modeled separately based on a range of individual input parameters and then combined to produce an integrated simulation of the entire eruption. The trajectory and areal extent of the modeled atmospheric ash cloud best match the actual distribution of MSH ash when input parameters are set to values inferred from satellite and radar data collected on May 18, 1980. The prevailing wind field exerts the strongest control on the advection and ultimate position of the modeled ash cloud, making the maximum column height and the vertical distribution of ash the most sensitive of the PUFF input parameters for this event. The results indicate that the PUFF model works well at simulating the dispersal of ash injected well into the lower stratosphere from a moderate, relatively long-lived eruption, such as MSH. However, attempts to use PUFF to recreate some granulometric aspects of the MSH fallout deposit, such as the maximum particle size as a function of distance from source, were not successful. PUFF consistently predicts much greater fallout distances for small ash particles (< 500 µm) than actually observed in the MSH deposit. The effective settling velocities used by the PUFF model appear to be too slow to accurately predict fallout distances of small ash particles. As a consequence the PUFF model may overestimate the duration of ash loading in the atmosphere associated with the distal fine ash component of explosive eruptions.