Thermal-image-derived dynamics of vertical ash plumes at Santiaguito volcano,Guatemala

Vertical ash plumes were imaged at Santiaguito (Guatemala) using a thermal camera to capture plume ascent dynamics. The plumes comprised a convecting plume front fed by a steady feeder plume. Of the 25 plumes imaged, 24 had a gas thrust region within which ascent velocities were 15–50 m s⁻¹. A transition to buoyant ascent occurred 20 to 50 m above the vent, where ascent velocities declined to 4–15 m s⁻¹. Plumes that attained greater heights had higher heat contents, wider feeder plumes and higher buoyant ascent velocities.     

Co-ignimbrite ash-fall deposits of the 1991 eruptions of Fugen-dake, Unzen Volcano, Japan

Fugen-dake, the main peak of Unzen Volcano, began a new eruption sequence on November 17, 1990. On May 20, 1991, a new lava dome appeared near the eastern edge of the Fugen-dake summit. Small-scale, 10⁴–10⁶ m³ in volume, Merapi-type block and ash flows were frequently generated from the growing lava dome during May–June, 1991. These pyroclastic flows were accompanied by co-ignimbrite ash plumes that deposited ash-fall deposits downwind of the volcano. Three examples of co-ignimbrite ash-fall deposits from Unzen pyroclastic flows are described. The volume of fall deposits was estimated to be about 30% by volume of the collapsed portions of the dome that formed pyroclastic flows. This proportion is smaller than that described for other larger co-ignimbrite ash-fall deposits from other volcanoes. Grain size distributions of the Unzen co-ignimbrite ash-fall deposits are bi-modal or tri-modal. Most ashes are finer than 4 phi and two modes were observed at around 4–7 phi and 9 phi. They are composed mainly of groundmass fragments. Fractions of another mode at around 2 phi are rich in crystals derived from dome lava. Some of the fine ash component fell as accretionary lapilli from the co-ignimbrite ash cloud indicating either moisture or electrostatic aggregation. We believe that the co-ignimbrite ash of Unzen block and ash flows were formed by the mechanical fracturing of the cooling lava blocks as they collapsed and moved down the slope. These ashes were entrained into the convective plumes generated off the tops of the moving flows.     

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     

Plinian and co-ignimbrite tephra fall from the

A study of pyroclastic deposits from the 1815 Tambora eruption reveals two distinct phases of activity, i.e., four initial tephra falls followed by generation of pyroclastic flows and the production of major co-ignimbrite ash fall. The first explosive event produced minor ash fall from phreatomagmatic explosions (F-1 layer). The second event was a Plinian eruption (F-2) correlated to the large explosion of 5 April 1815, which produced a column height of 33 km with an eruption rate of 1.1 × 10⁸ kg/s. The third event occurred during the lull in major activity from 5 to 10 April and produced minor ash fall (F-3). The fourth event produced a 43-km-high Plinian eruption column with an eruption rate of 2.8 × 10⁸ kg/s during the climax of activity on 10 April. Although very energetic, the Plinian events were of short duration (2.8 h each) and total erupted volume of the early (F-1 to F-4) fall deposits is only 1.8 km³ (DRE, dense rock equivalent). An abrupt change in style of activity occurred at end of the second Plinian event with onset of pyroclastic flow and surge generation. At least seven pyroclastic flows were generated, which spread over most of the volcano and Sanggar peninsula and entered the ocean. The volume of pyroclastic flow deposits on land is 2.6 km³ DRE. Coastal exposures show that pyroclastic flows entering the sea became highly fines depleted, resulting in mass loss of about 32%, in addition to 8% glass elutriation, as indicated by component fractionation. The subaqueous pyroclastic flows have thus lost about 40% of mass compared to the original erupted mixture. Pyroclastic flows and surges from this phase of the eruption are stratigraphically equivalent to a major ash fall deposit (F-5) present beyond the flow and surge zone at 40 km from the source and in distal areas. The F-5 fall deposit forms a larger proportion of the total tephra fall with increasing distance from source and represents about 80% of the total at a distance of 90 km and 92% of the total tephra fall from the 1815 eruption. The field relations indicate that the 20-km³ (DRE) F-5 deposit is a co-ignimbrite ash fall, generated largely during entrance of pyroclastic flows into the ocean. Based on the observed 40% fines depletion and component fractionation from the flows, the large volume of the F-5 co-ignimbrite ash requires eruption of 50 km³ (DRE, 1.4 × 10¹⁴ kg) pyroclastic flows.     

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 gas content and buoyancy of strombolian ash plumes

Plinian plumes erupt with a bulk density greater than that of air, and depend upon air entrainment during their gas-thrust phase to become buoyant; if entrainment is insufficient, the column collapses into a potentially deadly pyroclastic flow. This study shows that strombolian ash plumes can be erupted in an initially buoyant state due to their extremely high initial gas content, and in such cases are thus impervious to column collapse. The high gas content is a consequence of decoupled gas rise in the conduit, in which particles are ultimately incidental. The relations between conduit gas flow, eruption style and plume density are explored here for strombolian scenarios and contrasted with conventional wisdom derived from plinian eruptions. Considering the inherent relation between gas content and initial plume density together with detailed measurements of plume velocities can help unravel ambiguities surrounding conduit processes, eruption styles and hazards at poorly understood volcanoes. Analysis of plume dynamics at Santiaguito volcano, Guatemala adds further support for a model involving decoupled gas rise in the conduit.     

The significance of vitric-enriched air-fall ashes associated with crystal-enriched ignimbrites

A distinctive type of fine-grained air-fall ash is found intimately associated with many ignimbrites. They have crystal/glass ratios systematically lower than artificially crushed pumice from the same ignimbrites. The crystal enrichment found in crystal-bearing ignimbrites indicates substantial losses of the vitric component, amounting to an average of at least 35% by weight of the original juvenile material, and this lost material is believed to occur in the ash-fall deposits. These ashes thus complement ignimbrite, and are here called “co-ignimbrite ashes”. The loss is believed to take place during ignimbrite eruptions as a result of: (1) the escape of fine ash and gas above a collapsing eruptive column; (2) the preferential entry of fine vitric ash into an upper turbulent cloud when (immediately following column collapse) the segregation of a dense pyroclastic flow from an initially highly turbulent, low-concentation density flow takes place; (3) the elutriation of fine vitric ash (generated in part within the pyroclastic flow) from the fluidised flow. Ash from all three mechanisms would be expected to rise to a great height in convective plumes and be dispersed by winds to produce extensive, vitric-enriched ash-fall deposits.The data indicate that the co-ignimbrite ashes must have volumes comparable with those of ignimbrites, and examples are given of particularly large ash-fall deposits (including some found in deep-sea cores) associated with large ignimbrites which may be of this type rather than fall-out from a preceding plinian phase as hitherto thought.     

Lateral variations within coarse co-ignimbrite lithic breccias of the Kos Plateau Tuff, Greece

 Coarse, co-ignimbrite lithic breccia, Ebx, occurs at the base of ignimbrite E, the most voluminous and widespread unit of the Kos Plateau Tuff (KPT) in Greece. Similar but generally less coarse-grained basal lithic breccias (Dbx) are also associated with the ignimbrites in the underlying D unit. Ebx shows considerable lateral variations in texture, geometry and contact relationships but is generally less than a few metres thick and comprises lithic clasts that are centimetres to a few metres in diameter in a matrix ranging from fines bearing (F2: 10 wt.%) to fines poor (F2: 0.1 wt.%). Lithic clasts are predominantly vent-derived andesite, although clasts derived locally from the underlying sedimentary formations are also present. There are no proximal exposures of KPT. There is a highly irregular lower erosional contact at the base of ignimbrite E at the closest exposures to the inferred vent, 10–14 km from the centre of the inferred source, but no Ebx was deposited. From 14 to <20 km from source, Ebx is present over a planar erosional contact. At 16 km Ebx is a 3-m-thick, coarse, fines-poor lithic breccia separated from the overlying fines-bearing, pumiceous ignimbrite by a sharp contact. This grades downcurrent into a lithic breccia that comprises a mixture of coarse lithic clasts, pumice and ash, or into a thinner one-clast-thick lithic breccia that grades upward into relatively lithic-poor, pumiceous ignimbrite. Distally, 27 to <36 km from source Ebx is a finer one-clast-thick lithic breccia that overlies a non-erosional base. A downcurrent change from strongly erosional to depositional basal contacts of Ebx dominantly reflects a depletive pyroclastic density current. Initially, the front of the flow was highly energetic and scoured tens of metres into the underlying deposits. Once deposition of the lithic clasts began, local topography influenced the geometry and distribution of Ebx, and in some cases Ebx was deposited only on topographic crests and slopes on the lee-side of ridges. The KPT ignimbrites also contain discontinuous lithic-rich layers within texturally uniform pumiceous ignimbrite. These intra-ignimbrite lithic breccias are finer grained and thinner than the basal lithic breccias and overlie non-erosional basal contacts. The proportion of fine ash within the KPT lithic breccias is heterogeneous and is attributed to a combination of fluidisation within the leading part of the flow, turbulence induced locally by interaction with topography, flushing by steam generated by passage of pyroclastic density currents over and deposition onto wet mud, and to self-fluidisation accompanying the settling of coarse, dense lithic clasts. There are problems in interpreting the KPT lithic breccias as conventional co-ignimbrite lithic breccias. These problems arise in part from the inherent assumption in conventional models that pyroclastic flows are highly concentrated, non-turbulent systems that deposit en masse. The KPT coarse basal lithic breccias are more readily interpreted in terms of aggradation from stratified, waning pyroclastic density currents and from variations in lithic clast supply from source. Received: 21 April 1997 / Accepted: 4 October 1997     

Event-stratigraphy of a caldera-forming ignimbrite eruption on Tenerife: the 273 ka Poris Formation

The 273 ka Poris Formation in the Bandas del Sur Group records a complex, compositionally zoned explosive eruption at Las Cañadas caldera on Tenerife, Canary Islands. The eruption produced widespread pyroclastic density currents that devastated much of the SE of Tenerife, and deposited one of the most extensive ignimbrite sheets on the island. The sheet reaches ~ 40-m thick, and includes Plinian pumice fall layers, massive and diffuse-stratified pumiceous ignimbrite, widespread lithic breccias, and co-ignimbrite ashfall deposits. Several facies are fines-rich, and contain ash pellets and accretionary lapilli. Eight brief eruptive phases are represented within its lithostratigraphy. Phase 1 comprised a fluctuating Plinian eruption, in which column height increased and then stabilized with time and dispersed tephra over much of the southeastern part of the island. Phase 2 emplaced three geographically restricted ignimbrite flow-units and associated extensive thin co-ignimbrite ashfall layers, which contain abundant accretionary lapilli from moist co-ignimbrite ash plumes. A brief Plinian phase (Phase 3), again dispersing pumice lapilli over southeastern Tenerife, marked the onset of a large sustained pyroclastic density current (Phase 4), which then waxed (Phase 5), covering increasingly larger areas of the island, as vents widened and/or migrated along opening caldera faults. The climax of the Poris eruption (Phase 6) was marked by widespread emplacement of coarse lithic breccias, thought to record caldera subsidence. This is inferred to have disturbed the magma chamber, causing mingling and eruption of tephriphonolite magma, and it changed the proximal topography diverting the pyroclastic density current(s) down the Güimar valley (Phase 7). Phase 8 involved post-eruption erosion and sedimentary reworking, accompanied by minor down-slope sliding of ignimbrite. This was followed by slope stabilization and pedogenesis. The fines-rich lithofacies with abundant ash pellets and accretionary lapilli record agglomeration of ash in moist ash plumes. They resemble phreatomagmatic deposits, but a phreatomagmatic origin is difficult to establish because shards are of bubble-wall type, and the moisture may have arisen by condensation within ascending thermal co-ignimbrite ash plumes that contained atmospheric moisture enhanced by that derived from the evaporation of seawater where the hot pyroclastic currents crossed the coast. Ash pellets formed in co-ignimbrite ash-clouds and then fell through turbulent pyroclastic density currents where they accreted rims and evolved into accretionary lapilli.Editorial Responsibility: J. Stix     

Characterisation of the 1997 Vulcanian explosions of Soufrière Hills Volcano,Montserrat, by video analysis

The activity of Convention at Montserrat Soufrière Hills Volcano, Montserrat, during the period 1995–1999 included numerous violent explosions. Two major cycles of Vulcanian explosions occurred in 1997: a first of 13 explosions between 4 and 12 August and a second of 75 between 22 September and 21 October. The explosions were short-lived events lasting a few tens of seconds during which partial fountain collapse generated pyroclastic surges and pyroclastic flows, and buoyant plumes ascended 3–15 km into the atmosphere. Each explosion discharged on average 3×10⁵ m³ (dense-rock equivalent, DRE) of magma, draining the conduit to depths of 1–2 km. The paper focuses on the first few seconds of three explosions of the 75 that occurred in September/October 1997: 6 October 1997 at 17:50, 7 October 1997 at 16:02 and 9 October 1997 at 12:32. Physical parameters such as exit velocities, magmatic water contents and magma pressures at fragmentation are estimated by following and modelling the ascent of individual momentum-dominated finger jets visible on videos during the initial stages of each explosion. The model treats each finger jet as an incompressible flow sustained by a steady flux of gas and particles during the few seconds of ascent, and produces results that compare favourably with those using a multiphase compressible code run using similar eruptive parameters. Each explosion reveals a progressive increase in eruptive intensity with time, jet exit velocities increasing from 40 m s^–1 at the beginning of the explosion up to 140 m s^–1 after a few seconds. Modelling suggests that the first magma to exit was largely degassed, whereas that discharged after a few seconds contained up to 2 wt% water. Magma overpressures up to ~10 MPa are estimated to have existed in the conduit immediately prior to each explosion. Progressive increases in jet exit velocity with time over the first few seconds of each explosion provide direct evidence for strong pre-eruptive gradients in water content and magma pressure in the upper reaches (probably 100–500 m) of the conduit. Fountain collapse occurred during the first 10–20 s of each explosion because the discharging jets had bulk densities up to 100 times that of the atmosphere and were unable to entrain enough air to become buoyant. Such high eruptive densities were due to the presence of partially degassed magma in the conduit.Editorial responsibility: A. Woods     

Strombolian explosive styles and source conditions: insights from thermal (FLIR) video

Forward Looking Infrared Radiometer (FLIR) cameras offer a unique view of explosive volcanism by providing an image of calibrated temperatures. In this study, 344 eruptive events at Stromboli volcano, Italy, were imaged in 2001–2004 with a FLIR camera operating at up to 30 Hz. The FLIR was effective at revealing both ash plumes and coarse ballistic scoria, and a wide range of eruption styles was recorded. Eruptions at Stromboli can generally be classified into two groups: Type 1 eruptions, which are dominated by coarse ballistic particles, and Type 2 eruptions, which consist of an optically-thick, ash-rich plume, with (Type 2a) or without (Type 2b) large numbers of ballistic particles. Furthermore, Type 2a plumes exhibited gas thrust velocities (>15 m s⁻¹) while Type 2b plumes were limited to buoyant velocities (<15 m s⁻¹) above the crater rim. A given vent would normally maintain a particular gross eruption style (Type 1 vs. 2) for days to weeks, indicating stability of the uppermost conduit on these timescales. Velocities at the crater rim had a range of 3–101 m s⁻¹, with an overall mean value of 24 m s⁻¹. Mean crater rim velocities by eruption style were: Type 1 = 34 m s⁻¹, Type 2a = 31 m s⁻¹, Type 2b = 7 m s⁻¹. Eruption durations had a range of 6–41 s, with a mean of 15 s, similar among eruption styles. The ash in Type 2 eruptions originates from either backfilled material (crater wall slumping or ejecta rollback) or rheological changes in the uppermost magma column. Type 2a and 2b behaviors are shown to be a function of the overpressure of the bursting slug. In general, our imaging data support a broadening of the current paradigm for strombolian behavior, incorporating an uppermost conduit that can be more variable than is commonly considered.     

Numerical study of pyroclastic surges

The dynamics of pyroclastic surges accompanied by co-ignimbrite plumes is investigated numerically. The numerical simulations are performed with a newly developed numerical model, which is based on the Navier–Stokes equations for time-dependent flows of a compressible fluid in two-dimensional Cartesian coordinates. We regard pyroclastic surges as dilute turbulent suspensions in which hot gases and fine solid particles are homogeneously mixed owing to vigorous turbulence. In other words, the gas–particle mixture is treated as a single-phase fluid whose bulk density is represented by averaging the density of each component in the numerical model. We focus on the effect of buoyancy forces generated by the thermal expansion of the air mixed into pyroclastic surges from the calm surroundings. For our purpose, the numerical model is designed to simulate relatively simple flows spreading over a horizontal flat surface. Topographic irregularity and the sedimentation process of solid particles are neglected in the present simulations. The motion of pyroclastic surges is generated by the instantaneous release of a gas–particle mixture whose density is initially larger than the ambient air density and changes nonlinearly with the temperature and concentration of suspended solid particles. Turbulent mixing is evaluated by adopting the Smagorinsky model. By employing cubic interpolated pseudo-particle (CIP) method and C-CUP method, we obtain the fine structure of flows. The behavior of calculated flows agrees fairly well with observed pyroclastic surges in nature. The current head, which remains hot and dense, keeps spreading over a horizontal surface at a speed of about 20 m s⁻¹. The spreading speed is of the order of the speed of a gravity current that excludes the influence of thermal expansion. Besides, turbulent mixing between the basal dense layer and the ambient air is enhanced by the successive development of an interfacial less-dense layer. This results in the formation of a number of buoyant plumes rising above a horizontally spreading current. Consequently, the tails of the current thickens as time progresses. A parametric study shows that the initial temperature of a gas–particle mixture should be higher than about 600 K when buoyant plumes occur owing to the thermal expansion of mixed air. The result is quantitatively interpreted by introducing a diagram that describes the relationship among the bulk density, temperature and concentration of solid particles suspended in pyroclastic surges.     

Internal structure and occurrence of accretionary lapilli — a case study at Laacher See Volcano

Accretionary lapilli are common in fine-grained pyroclastic flow and surge deposits and related co-ignimbrite/co-surge ash layers of Laacher See volcano. Two morphologically different types are distin-guished: (1) Rim-type lapilli are composed of a coarse-grained core surrounded by a fine-grained rim. Rims are internally graded or made up of several layers of alternating fine and very-fine grained ash. (2) Core-type lapilli lack fine-grained rims. Field relationships, internal, and grain-size characteristics are specific to accretionary lapilli from different types of tephra deposits. Accretionary lapilli may therefore be a helpful tool to infer the origin of tephra of different origin. In co-ignimbrite ashfall, accretionary lapilli are generally concentrated at the base, whereas pyroclastic flow and surge deposits contain lapilli in the upper parts of individual, thin-bedded layers. Rim-type lapilli are found in pyroclastic flow and surge deposits up to 4 km from the source. Core-type lapilli occur at greater distances or are associated with vesiculated tuffs where they are within 1 km from the vent. Accretionary lapilli from co-ignimbrite/co-surge ash show open framework textures and edge-to-face contacts of individual ash particles. Vesicularity is generally low but the overall porosity of 40% to 50% results in an average density of 1200 kg/m³. Accretionary lapilli in pyroclastic flow and surge deposits are more densely packed and platy particles are often in face-to-face contacts. Vesicularity of those from pyroclastic flow deposits is significantly higher; the overall porosity is about 30% to 40% and the average density 1600 kg/m³. Grain-size analyses show that the accretionary lapilli in co-ignimbrite/co-surge ashfall deposits are the most fine-grained with a median (Md) of 20 to 30 m and a maximum grain size of 250 to 350 m. Accretionary lapilli from pyroclastic flow deposits have intermediate Md-values of 30 to 50 m and a maximum grain size of 350 to 500 m. Those of surge deposits are the coarsest grained with Md-values of 30 to >63 m and a maximum grain size up to 2 mm.     

Incursion of a large-volume,spatter-bearing pyroclastic density current into a caldera lake: Pavey Ark ignimbrite,Scafell caldera,England

The Scafell caldera-lake volcaniclastic succession is exceptionally well exposed. At the eastern margin of the caldera, a large andesitic explosive eruption (>5 km³) generated a high-mass-flux pyroclastic density current that flowed into the caldera lake for several hours and deposited the extensive Pavey Ark ignimbrite. The ignimbrite comprises a thick (≤125 m), proximal, spatter- and scoria-rich breccia that grades laterally and upwards into massive lapilli-tuff, which, in turn, is gradationally overlain by massive and normal-graded tuff showing evidence of soft-state disruption. The subaqueous pyroclastic current carried juvenile clasts ranging from fine ash to metre-scale blocks and from dense andesite through variably vesicular scoria to pumice (<10³ kg m⁻³). Extreme ignimbrite lithofacies diversity resulted via particle segregation and selective deposition from the current. The lacustrine proximal ignimbrite breccia mainly comprises clast- to matrix-supported blocks and lapilli of vesicular andesite, but includes several layers rich in spatter (≤1.7 m diameter) that was emplaced in a ductile, hot state. In proximal locations, rapid deposition of the large and dense clasts caused displacement of interstitial fluid with elutriation of low-density lapilli and ash upwards, so that these particles were retained in the current and thus overpassed to medial and distal reaches. Medially, the lithofacies architecture records partial blocking, channelling and reflection of the depletive current by substantial basin-floor topography that included a lava dome and developing fault scarps. Diffuse layers reflect surging of the sustained current, and the overall normal grading reflects gradually waning flow with, finally, a transition to suspension sedimentation from an ash-choked water column. Fine to extremely fine tuff overlying the ignimbrite forms ∼25% of the whole and is the water-settled equivalent of co-ignimbrite ash; its great thickness (≤55 m) formed because the suspended ash was trapped within an enclosed basin and could not drift away. The ignimbrite architecture records widespread caldera subsidence during the eruption, involving volcanotectonic faulting of the lake floor. The eruption was partly driven by explosive disruption of a groundwater-hydrothermal system adjacent to the magma reservoir.     

Volcanic stratigraphy of large-volume silicic pyroclastic eruptions during Oligocene Afro-Arabian flood volcanism in Yemen

A new stratigraphy for bimodal Oligocene flood volcanism that forms the volcanic plateau of northern Yemen is presented based on detailed field observations, petrography and geochemical correlations. The >1 km thick volcanic pile is divided into three phases of volcanism: a main basaltic stage (31 to 29.7 Ma), a main silicic stage (29.7 to 29.5 Ma), and a stage of upper bimodal volcanism (29.5 to 27.7 Ma). Eight large-volume silicic pyroclastic eruptive units are traceable throughout northern Yemen, and some units can be correlated with silicic eruptive units in the Ethiopian Traps and to tephra layers in the Indian Ocean. The silicic units comprise pyroclastic density current and fall deposits and a caldera-collapse breccia, and they display textures that unequivocally identify them as primary pyroclastic deposits: basal vitrophyres, eutaxitic fabrics, glass shards, vitroclastic ash matrices and accretionary lapilli. Individual pyroclastic eruptions have preserved on-land volumes of up to ∼850 km³. The largest units have associated co-ignimbrite plume ash fall deposits with dispersal areas >1×10⁷ km² and estimated maximum total volumes of up to 5,000 km³, which provide accurate and precisely dated marker horizons that can be used to link litho-, bio- and magnetostratigraphy studies. There is a marked change in eruption style of silicic units with time, from initial large-volume explosive pyroclastic eruptions producing ignimbrites and near-globally distributed tuffs, to smaller volume (<50 km³) mixed effusive-explosive eruptions emplacing silicic lavas intercalated with tuffs and ignimbrites. Although eruption volumes decrease by an order of magnitude from the first stage to the last, eruption intervals within each phase remain broadly similar. These changes may reflect the initiation of continental rifting and the transition from pre-break-up thick, stable crust supporting large-volume magma chambers, to syn-rift actively thinning crust hosting small-volume magma chambers.Electronic Supplementary Material Supplementary material is available for this article at     

Causes and consequences of bimodal grain-size distribution of tephra fall deposited during the August 2006 Tungurahua eruption (Ecuador)

The violent August 16–17, 2006 Tungurahua eruption in Ecuador witnessed the emplacement of numerous scoria flows and the deposition of a widespread tephra layer west of the volcano. We assess the size of the eruption by determining a bulk tephra volume in the range 42–57 × 10⁶ m³, which supports a Volcanic Explosivity Index 3 event, consistent with calculated column height of 16–18 km above the vent and making it the strongest eruptive phase since the volcano’s magmatic reactivation in 1999. Isopachs west of the volcano are sub-bilobate in shape, while sieve and laser diffraction grain-size analyses of tephra samples reveal strongly bimodal distributions. Based on a new grain-size deconvolution algorithm and extended sampling area, we propose here a mechanism to account for the bimodal grain-size distribution. The deconvolution procedure allows us to identify two particle subpopulations in the deposit with distinct characteristics that indicate dissimilar transport-depositional processes. The log-normal coarse-grained subpopulation is typical of particles transported downwind by the main volcanic plume. The positively skewed, fine-grained subpopulation in the tephra fall layer shares close similarities with the elutriated co-pyroclastic flow ash cloud layers preserved on top of the scoria flow deposits. The area with the higher fine particle content in the tephra layer coincides with the downwind prolongation of the pyroclastic flow deposits. These results indicate that the bimodal distribution of grain size in the Tungurahua fall deposit results from synchronous deposition of lapilli from the main plume and fine ash elutriated from scoria flows emplaced on the western flank of the volcano. Our study also reveals that inappropriate grain-size data processing may produce misleading determination of eruptive type.     

Volcanic eruption clouds and the thermal power output of explosive eruptions

The maximum height attained by a volcanic eruption cloud is principally determined by the convective buoyancy of the mixture of volcanic gas + entrained air + fine-sized pyroclasts within the cloud. The thermal energy supplied to convection processes within an eruption cloud is derived from the cooling of pyroclastic material and volcanic gases discharged by an explosive eruption. Observational data from six recent eruptions indicates that the maximum height attained by volcanic eruption clouds is positively correlated with the rate at which pyroclastic material is produced by an explosive eruption (correlation coefficient r = + 0.97). The ascent of industrial hot gas plumes is also governed by the thermal convection process. Empirical scaling relationships between plume height and thermal flux have been developed for industrial plumes. Applying these scaling relationships to volcanic eruption clouds suggests that the rate at which thermal energy is released into the atmosphere by an explosive eruption increases in an approximately linear manner as an eruption's pyroclastic production rate increases.     

Estimates of eruption velocity and plume height from infrasonic recordings of the 2006 eruption of Augustine Volcano,Alaska

The 2006 eruption of Augustine Volcano, Alaska, began with an explosive phase comprising 13 discrete Vulcanian blasts. These events generated ash plumes reaching heights of 3–14 km. The eruption was recorded by a dense geophysical network including a pressure sensor located 3.2 km from the vent. Infrasonic signals recorded in association with the eruptions have maximum pressures ranging from 13–111 Pa. Eruption durations are estimated to range from 55–350 s. Neither of these parameters, however, correlates with eruption plume height. The pressure record, however, can be used to estimate the velocity and flux of material erupting from the vent, assuming that the sound is generated as a dipole source. Eruptive flux, in turn, is used to estimate plume height, assuming that the plume rises as a buoyant thermal. Plume heights estimated in this way correlate well with observations. Events that exhibit strongly impulsive waveforms are underestimated by the model, suggesting that flow may have been supersonic.     

The Holocene Pucón eruption of Volcán Villarrica, Chile: deposit architecture and eruption chronology

The Pucón eruption was the largest Holocene explosive outburst of Volcán Villarrica, Chile. It discharged >1.0 km³ of basaltic-andesite magma and >0.8 km³ of pre-existing rock, forming a thin scoria-fall deposit overlain by voluminous ignimbrite intercalated with pyroclastic surge beds. The deposits are up to 70 m thick and are preserved up to 21 km from the present-day summit, post-eruptive lahar deposits extending farther. Two ignimbrite units are distinguished: a lower one (P1) in which all accidental lithic clasts are of volcanic origin and an upper unit (P2) in which basement granitoids also occur, both as free clasts and as xenoliths in scoria. P2 accounts for ∼80% of the erupted products. Following the initial scoria fallout phase, P1 pyroclastic flows swept down the northern and western flanks of the volcano, magma fragmentation during this phase being confined to within the volcanic edifice. Following a pause of at least a couple of days sufficient for wood devolatilization, eruption recommenced, the fragmentation level dropped to within the granitoid basement, and the pyroclastic flows of P2 were erupted. The first P2 flow had a highly turbulent front, laid down ignimbrite with large-scale cross-stratification and regressive bedforms, and sheared the ground; flow then waned and became confined to the southeastern flank. Following emplacement of pyroclastic surge deposits all across the volcano, the eruption terminated with pyroclastic flows down the northern flank. Multiple lahars were generated prior to the onset of a new eruptive cycle. Charcoal samples yield a probable eruption age of 3,510 ± 60 ¹⁴C years BP.