Subaqueous pyroclastic flows and ignimbrites: an assessment

An assessment of the literature on subaqueous pyroclastic flows and their deposits shows that the term pyroclastic flow is frequently used loosely to describe primary, hot gas-rich pyroclastic flows, mass-flows which resulted from the transformation of gassupported flows into water-supported ones, and secondary mass-flows carrying redeposited pyroclastic debris. Based on subaerial pyroclastic flows, the term pyroclastic flow should be restricted to demonstrably hot, gas-rich mass-flows of pyroclastic debris. Using this definition, very few examples of subaqueous pyroclastic deposits with evidence for hot emplacement and of having been wholly submerged have been described. In the majority of these cases, the evidence for a hot state of emplacement and for the subaqueous nature of the host depositional environment is inadequate. The only unequivocal cases of hot pyroclastic flow deposits with adequate supporting evidence are the Ordovician nearshore, shallow marine ignimbrites of Ireland and Wales, and Miocene ignimbrites of southwest Japan, resulting from the passage of subaerially erupted pyroclastic flows into shallow water. Other possible examples are near-vent dense clast deposits in the Donzurobo Formation of Japan, possible submarine intra-caldera ponded ignimbrite successions in California and Wales, and near-vent pumiceous deposits of Ramsay Island, Wales. All other purported cases are either clearly the result of water-supported mass-flow transportation and deposition (debris avalanches, debris flows, turbidity currents), or lack adequate supporting evidence regarding the heat state or the palaeoenvironment. Only the shallow marine ignimbrites of Ireland and Wales show adequate evidence of welding, but even these could have been nearly wholly exposed above sea-level when welding occurred. We conclude that when pyroclastic flows enter water they are generally disrupted explosively and/or ingest water and transform into water-supported mass-flows, and we suggest the various scenarios in which this occurs. There is no evidence to suggest that welding in wholly subaqueous environments is common.     

The role of fluidization in the emplacement of pyroclastic claws: An experimental approach

A series of experiments was carried out to review the process of fluidization for a number of particulate materials having various sorting and grain shape characteristics. The up (increasing gas velocity) curve on a gas velocity/bed-pressure drop plot for a poorly sorted mixture of irregularly shaped particles is divided into three sections; non-expanded, expanded, and segregating. These sections are used to define a threefold genetic classification of pyroclastic flows which can be directly linked to conditions within a semifluidized parent flow. Type 1 flows are ungraded and mostly result from hot-avalanche flows and pyroclastic flows formed of relatively dense material. Type 2 and type 3 flows are mainly pumiceous ignimbrites and are distinguished by expansion induced coarse-tail grading, coupled with segregation structures in the latter. The implications of this classification are discussed with reference to the flow regimes, deposition slope angles, crystal concentration and “fossil fumaroles” variously developed in the flows. The relevance of the source of fluidizing gas is discussed in relation to the zoning of flow types within a single flow unit and it is shown how this, and its attendant structures, can be used to estimate the relative importance of each gas source during and after flow emplacement.     

The entrance of pyroclastic flows into the sea, II. theoretical considerations on subaqueous emplacement and welding

The submarine counterparts of late Quaternary subaerial pyroclastic flow deposits off the western flanks of Dominica, Lesser Antilles, have been investigated by 3.5 kHz seismic profiling and dredging (cruise EN20 of R/V “Endeavor”). Block-and-ash flow deposits formed by dome collapse and a welded ignimbrite from a prominent fan at Grande Savanne, Dominica. This fan can be traced underwater as a major constructional ridge (2–4 km wide and 200–400 m thick) to over 13 km offshore at a water depth of 1800 m. The submarine ridge has a volume of 14 km³ and has the characteristic morphology of a debris flow apron composed of several individual units. The evidence suggests that pyroclastic flows can move underwater without losing their essential character.     

Modeling of pyroclastic flows of Colima Volcano,Mexico: implications for hazard assessment

The 18–24 January 1913 eruption of Colima Volcano consisted of three eruptive phases that produced a complex sequence of tephra fall, pyroclastic surges and pyroclastic flows, with a total volume of 1.1 km³ (0.31 km³ DRE). Among these events, the pyroclastic flows are most interesting because their generation mechanisms changed with time. They started with gravitanional dome collapse (block-and-ash flow deposits, Merapi-type), changed to dome collapse triggered by a Vulcanian explosion (block-and-ash flow deposits, Soufrière-type), then ended with the partial collapse of a Plinian column (ash-flow deposits rich in pumice or scoria,). The best exposures of these deposits occur in the southern gullies of the volcano where Heim Coefficients (H/L) were obtained for the various types of flows. Average H/L values of these deposits varied from 0.40 for the Merapi-type (similar to the block-and-ash flow deposits produced during the 1991 and 1994 eruptions), 0.26 for the Soufrière-type events, and 0.17–0.26 for the column collapse ash flows. Additionally, the information of 1991, 1994 and 1998–1999 pyroclastic flow events was used to delimit hazard zones. In order to reconstruct the paths, velocities, and extents of the 20th Century pyroclastic flows, a series of computer simulations were conducted using the program FLOW3D with appropriate Heim coefficients and apparent viscosities. The model results provide a basis for estimating the areas and levels of hazard that could be associated with the next probable worst-case scenario eruption of the volcano. Three areas were traced according to the degree of hazard and pyroclastic flow type recurrence through time. Zone 1 has the largest probability to be reached by short runout (<5 km) Merapi and Soufrière pyroclastic flows, that have occurred every 3 years during the last decade. Zone 2 might be affected by Soufriere-type pyroclastic flows (∼9 km long) similar to those produced during phase II of the 1913 eruption. Zone 3 will only be affected by pyroclastic flows (∼15 km long) formed by the collapse of a Plinian eruptive column, like that of the 1913 climactic eruption. Today, an eruption of the same magnitude as that of 1913 would affect about 15,000 inhabitants of small villages, ranches and towns located within 15 km south of the volcano. Such towns include Yerbabuena, and Becerrera in the State of Colima, and Tonila, San Marcos, Cofradia, and Juan Barragán in the State of Jalisco.     

Basaltic pyroclastic flows of Fuji volcano, Japan: characteristics of the deposits and their origin

Fuji volcano is the largest active volcano in Japan, and consists of Ko-Fuji and Shin-Fuji volcanoes. Although basaltic in composition, small-volume pyroclastic flows have been repeatedly generated during the Younger stage of Shin-Fuji volcano. Deposits of those pyroclastic flows have been identified along multiple drainage valleys on the western flanks between 1,300 and 2,000 m a.s.l., and have been stratigraphically divided into the Shin-Fuji Younger pyroclastic flows (SYP) 1 to 4. Downstream debris flow deposits are found which contain abundant material derived from the pyroclastic flow deposits. The new¹⁴C ages for SYP1 to SYP4 are 3.2, 3.0, 2.9, and 2.5 ka, respectively, and correspond to a period where explosive summit eruptions generated many scoria fall deposits mostly toward the east. The SYP1 to SYP4 deposits consist of two facies: the massive facies is about 2 m thick and contains basaltic bombs of less than 50 cm in size, scoria lapilli, and fresh lithic basalt fragments supported in an ash matrix; the surge facies is represented by beds 1 to 15 cm thick, consisting mainly of ash with minor amount of fine lapilli. The bombs and scoria are 15 to 30% in volume within the massive facies. The ashes within the SYP deposits consist largely of comminuted basalt lithics and crystals that are derived from the Middle-stage lava flows exposed at the western flanks. SYP1 to SYP4 were only dispersed down the western flanks. The reason for this one-sided distribution is the asymmetric topography of the edifice; the western slopes of the volcano are the steepest (over 34 degrees). Most pyroclastic materials cannot rest stably on the slopes steeper than 33 degrees. Therefore, ejecta from the explosive summit eruptions that fell on the steep slopes tumbled down the slopes and were remobilized as high-temperature granular flows. These flows consisted of large pyroclastics and moved as granular avalanches along the valley bottom. Furthermore, the avalanching flows increased in volume by abrasion from the edifice and generated abundant ashes by the collision of clasts. The large amount of the fine material was presumably available within the transport system as the basal avalanches propagated below the angle of repose. Taking the typical kinetic friction coefficient of small pyroclastic flows, such flows could descend the western flanks where scattered houses are below 1,000 m a.s.l. A similar type of pyroclastic flow could result if explosive summit eruptions occur in the future.Editorial responsibility: R Cioni     

Paleomagnetic determination of emplacement temperatures of pyroclastic deposits: an under-utilized tool

Paleomagnetic data from lithic clasts collected from Mt. St. Helens, USA, Volcán Láscar, Chile, Volcán de Colima, Mexico and Vesuvius, Italy have been used to determine the emplacement temperature of pyroclastic deposits at these localities and to highlight the usefulness of the paleomagnetic method for determining emplacement temperatures. At Mt. St. Helens, the temperature of the deposits (T _dep ) at three sites from the June 12, 1980 eruption was found to be ≥532°C, ≥509°C, and 510–570°C, respectively. One site emplaced on July 22, 1980 was emplaced at ≥577°C. These new paleomagnetic temperatures are in good agreement with previously published direct temperature measurements and paleomagnetic estimates. Lithic clasts from pyroclastic deposits from the 1993 eruption of Láscar were fully remagnetized above the respective Curie temperatures, which yielded a minimum T _dep of 397°C. Samples were also collected from deposits thought to be pyroclastics from the 1913, 2004 and 2005 eruptions of Colima. At Colima, the sampled clasts were emplaced cold. This is consistent with the sampled clasts being from lahar deposits, which are common in the area, and illustrates the usefulness of the paleomagnetic method for distinguishing different types of deposit. T _dep of the lower section of the lithic rich pyroclastic flow (LRPF) from the 472 A.D. deposits of Vesuvius was ~280–340°C. This is in agreement with other, recently published paleomagnetic measurements. In contrast, the upper section of the LRPF was emplaced at higher temperatures, with T _dep ~520°C. This temperature difference is inferred to be the result of different sources of lithic clasts between the upper and lower sections, with the upper section containing a greater proportion of vent-derived material that was initially hot. Our studies of four historical pyroclastic deposits demonstrates the usefulness of paleomagnetism for emplacement temperature estimation.     

Morphology and emplacement of flows from the Deccan Volcanic Province,India

The present study is probably the first of its kind in the Deccan Volcanic Province (DVP) that deals in detail with the morphology and emplacement of the Deccan Trap flows, and employs modern terminology and concepts of flow emplacement. We describe in detail the two major types of flows that occur in this province. Compound pahoehoe flows, similar to those in Hawaii and the Columbia River Basalts (CRB) constitute the older stratigraphic Formations. These are thick flows, displaying the entire range of pahoehoe morphology including inflated sheets, hummocky flows, and tumuli. In general, they show the same three-part structure associated with pahoehoe flows from other provinces. However, in contrast to the CRB, pahoehoe lobes in the DVP are smaller, and hummocky flows are quite common. 'Simple' flows occur in the younger Formations and form extensive sheets capped by highly vesicular, weathered crusts, or flow-top breccias. These flows have few analogues in other provinces. Although considered to be a'a flows by previous workers, the present study clearly reveals that the simple flows differ considerably from typical a'a flows, especially those of the proximal variety. This is very significant in the context of models of flood basalt emplacement. At the same time, they do not display direct evidence of endogenous growth. Understanding the emplacement of these flows will go a long way in determining whether all extensive flows are indeed inflated flows, as has recently been postulated.Most of the studies relating to the emplacement of Continental Flood Basalt (CFB) lavas have relied on observations of flows from the CRB. Much of the current controversy surrounding the emplacement of CFB flows centers around the comparison of Hawaiian lava flows to those from the CRB. We demonstrate that the DVP displays a variety of lava features that are similar to those from the CRB as well as those from Hawaii. This suggests that there may have been more than one mechanism or style for the emplacement of CFB flows. These need to be taken into account before arriving at any general model for flood basalt emplacement.Editorial responsibility: T. Druitt     

Characteristics of submarine pumice-rich density current deposits sourced from turbulent mixing of subaerial pyroclastic flows at the shoreline: field and experimental assessment

This study investigates the types of subaqueous deposits that occur when hot pyroclastic flows turbulently mix with water at the shoreline through field studies of the Znp marine tephra in Japan and flume experiments where hot tephra sample interacted with water. The Znp is a very thick, pumice-rich density current deposit that was sourced from subaerial pyroclastic flows entering the Japan Sea in the Pliocene. Notable characteristics are well-developed grain size and density grading (lithic-rich base, pumice-rich middle, and ash-rich top), preponderance of sedimentary lithic clasts picked up from the seafloor during transport, fine ash depletion in coarse facies, and presence of curviplanar pumice clasts. Flume experiments provide a framework for interpreting the origin and proximity to source of the Znp tephra. On contact of hot tephra sample with water, steam explosions produced a gas-supported pyroclastic density current that advanced over the water while a water-supported density current was produced on the tank floor from the base of a turbulent mixing zone. Experimental deposits comprise proximal lithic breccia, medial pumice breccia, and distal fine ash. Experiments undertaken with cold, water-saturated slurries of tephra sample and water did not produce proximal lithic breccias but a medial basal lithic breccia beneath an upper pumice breccia. Results suggest the characteristics and variations in Znp facies were strongly controlled by turbulent mixing and quenching, proximity to the shoreline, and depositional setting within the basin. Presence of abundant curviplanar pumice clasts in submarine breccias reflects brittle fracture and dismembering that can occur during fragmentation at the vent or during quenching. Subsequent transport in water-supported pumiceous density currents preserves the fragmental textures. Careful study is needed to distinguish the products of subaerial versus subaqueous eruptions.     

Flow behavior of large-scale pyroclastic flows — Evidence obtained from petrofabric analysis

The grain orientations within the matrix of two large-scale welded, two small-scale nonwelded and two nonwelded low-aspect ratio pyroclastic flow deposits are measured to analyze flow behavior. Preferred grain alignments are especially apparent in the middle part of layer 2 of each deposit. Preferred grain alignments do not vary laterally within a 10 m interval. The grain alignments obtained are radial from the source caldera, especially in proximal to medial and plateau-forming facies of pyroclastic flow deposits. Grain alignments are controlled by valley-channel directions for the valley-ponded facies of pyroclastic flow deposits, especially at medial to distal locations. Such local topographic factors strongly affect the data for high-aspect ratio and smallscale deposits, and weakly affect the data for widespread low-aspect ratio pyroclastic flow deposits. This work suggests that grain alignment analysis should be used with care when attempting to determine the location of an unknown source.     

Inviscid behaviour of fines-rich pyroclastic flows inferred from experiments on gas-particle mixtures

Experiments were carried out on granular flows generated by instantaneous release of gas-fluidised, bidisperse mixtures and propagating into a horizontal channel. The mixture consists of fine (< 100 μm) and coarse (> 100 μm) particles of same density, with corresponding grain size ratios of ∼ 2 to 9. Initial fluidisation of the mixture destroys the interparticle frictional contacts, and the flow behaviour then depends on the initial bed packing and on the timescale required to re-establish strong frictional contacts. At a fines mass fraction (α) below that of optimal packing (∼ 40%), the initial mixtures consist of a continuous network of coarse particles with fines in interstitial voids. Strong frictional contacts between the coarse particles are probably rapidly re-established and the flows steadily decelerate. Some internal friction reduction appears to occur as α and the grain size ratio increases, possibly due to particle rolling and the lower roughness of internal shear surfaces. Segregation only occurs at large grain size ratio due to dynamical sieving with fines concentrated at the flow base. In contrast, at α above that for optimal packing, the initial mixtures consist of coarse particles embedded in a matrix of fines. Flow velocities and run-outs are similar to that of the monodisperse fine end-member, thus showing that the coarse particles are transported passively within the matrix whatever their amount and grain size are. These flows propagate at constant height and velocity as inviscid fluid gravity currents, thus suggesting negligible interparticle friction. We have determined a Froude number of 2.61 ± 0.08 consistent with the dam-break model for fluid flows, and with no significant variation as a function of α, the grain size ratio, and the initial bed expansion. Very little segregation occurs, which suggests low intensity particle interactions during flow propagation and that active fluidisation is not taking place. Strong frictional contacts are only re-established in the final stages of emplacement and stop the flow motion. We infer that fines-rich (i.e. matrix-supported) pyroclastic flows propagate as inviscid fluid gravity currents for most of their emplacement, and this is consistent with some field data.     

Objective rapid delineation of areas at risk from block-and-ash pyroclastic flows and surges

Assessments of pyroclastic flow (PF) hazards are commonly based on mapping of PF and surge deposits and estimations of inundation limits, and/or computer models of varying degrees of sophistication. In volcanic crises a PF hazard map may be sorely needed, but limited time, exposures, or safety aspects may preclude fieldwork, and insufficient time or baseline data may be available for reliable dynamic simulations. We have developed a statistically constrained simulation model for block-and-ash type PFs to estimate potential areas of inundation by adapting methodology from Iverson et al. (Geol Soc America Bull 110:972–984, 1998) for lahars. The predictive equations for block-and-ash PFs are calibrated with data from several volcanoes and given by A = (0.05 to 0.1)V ^2/3, B = (35 to 40)V ^2/3, where A is cross-sectional area of inundation, B is planimetric area and V is deposit volume. The proportionality coefficients were obtained from regression analyses and comparison of simulations to mapped deposits. The method embeds the predictive equations in a GIS program coupled with DEM topography, using the LAHARZ program of Schilling (1998). Although the method is objective and reproducible, any PF hazard zone so computed should be considered as an approximate guide only, due to uncertainties on the coefficients applicable to individual PFs, the authenticity of DEM details, and the volume of future collapses. The statistical uncertainty of the predictive equations, which imply a factor of two or more in predicting A or B for a specified V, is superposed on the uncertainty of forecasting V for the next PF to descend a particular valley. Multiple inundation zones, produced by simulations using a selected range of volumes, partly accommodate these uncertainties. The resulting maps show graphically that PF inundation potentials are highest nearest volcano sources and along valley thalwegs, and diminish with distance from source and lateral distance from thalweg. The model does not explicitly consider dynamic behavior, which can be important. Ash-cloud surge impact limits must be extended beyond PF hazard zones and we provide several approaches to do this. The method has been used to supply PF and surge hazard maps in two crises: Merapi 2006; and Montserrat 2006–2007.     

Depositional processes and gas pore pressure in pyroclastic flows: an experimental perspective

The depositional processes and gas pore pressure in pyroclastic flows are investigated through scaled experiments on transient, initially fluidized granular flows. The flow structure consists of a sliding head whose basal velocity decreases backwards from the front velocity (U _f) until onset of deposition occurs, which marks transition to the flow body where the basal deposit grows continuously. The flows propagate in a fluid-inertial regime despite formation of the deposit. Their head generates underpressure proportional to U _f ² whereas their body generates overpressure whose values suggest that pore pressure diffuses during emplacement. Complementary experiments on defluidizing static columns prove that the concept of pore pressure diffusion is relevant for gas-particle mixtures and allow characterization of the diffusion timescale (t _d) as a function of the material properties. Initial material expansion increases the diffusion time compared with the nonexpanded state, suggesting that pore pressure is self-generated during compaction. Application to pyroclastic flows gives minimum diffusion timescales of seconds to tens of minutes, depending principally on the flow height and permeability. This study also helps to reconcile the concepts of en masse and progressive deposition of pyroclastic flow units or discrete pulses. Onset of deposition, whose causes deserve further investigation, is the most critical parameter for determining the structure of the deposits. Even if sedimentation is fundamentally continuous, it is proposed that late onset of deposition and rapid aggradation in relatively thin flows can generate deposits that are almost snapshots of the flow structure. In this context, deposition can be considered as occurring en masse, though not strictly instantaneously.     

Flow and deposition of pyroclastic granular flows: A type example from the 1975 Ngauruhoe eruption,New Zealand

Small-volume pyroclastic density currents (PDCs) are generated frequently during explosive eruptions with little warning. Assessing their hazard requires a physical understanding of their transport and sedimentation processes which is best achieved by the testing of experimental and numerical models of geophysical mass flows against natural flows and/or deposits. To this end we report on one of the most detailed sedimentological studies ever carried out on a series of pristine small-volume PDC deposits from the 1975 eruption of Ngauruhoe volcano, whose emplacement were also witnessed during eruption. Using high-resolution GPS surveys, a series of lateral excavations across the deposits, and bulk sedimentological analysis we constrained the geomorphology, internal structure and texture of the deposits with respect to laterally varying modes of deposition.     

Morphology of rubbly pahoehoe (simple) flows from the Deccan Volcanic Province: Implications for style of emplacement

Lava flows with preserved bases and brecciated upper crusts constitute a morphological type that differs in character from typical pahoehoe and a'a: such flows have been reported from many provinces around the world. Previous studies had referred to these flows informally as ‘pahoehoe flows with rubbly tops’, ‘broken-top pahoehoe’ and ‘rubbly pahoehoe’. Recent studies have formally applied the latter term to describe parts of the well-studied Laki flow in Iceland as well as flows from the Columbia River Basalt province. Rubbly pahoehoe flows are abundant in the upper stratigraphic formations of the Deccan Volcanic Province (DVP), and are more commonly known as simple flows. This study presents detailed observations of such flows from various parts of the DVP and discusses their implications for understanding flow emplacement. These flows, which appear to be single units at the outcrop-scale, are generally much thicker and significantly more extensive than individual pahoehoe lobes that dominate the lower formations of the Deccan stratigraphy. They are characterised by preserved, gently undulating tachylitic bases but variably disrupted crustal zones that grade into flow-top breccias. The breccias are constituted of highly vesicular and oxidised fragments of varying sizes that appear to have been derived from previously formed pahoehoe crusts. Previous work has indicated that the morphology of these flows might be related to initial inflation, accompanied by rapid volatile exsolution and an increase in effusion rate and/or viscosity with time. This agrees reasonably well with the qualitative and quantitative models of emplacement developed for the Laki flow. The abundance of such flows in the upper formations of the Deccan stratigraphy clearly hints at a significant shift in the nature of the Deccan eruptions; this could be indicative of higher eruption rates during this period. This, in turn, raises the possibility of hazardous impact on the climate during the eruption of these flows, which is also discussed in the paper.     

A methodology for the evaluation of long-term volcanic risk from pyroclastic flows in Campi Flegrei (Italy)

The volcanological history of Campi Flegrei suggests that the most frequent eruptions are characterized by the emplacement of pyroclastic flow and surge deposits erupted from different vents scattered over a 150-km² caldera. The evaluation of volcanic risk in volcanic fields is complex because of the lack of a central vent. To approach this problem, we subdivided the entire area of Campi Flegrei into a regular grid and evaluated the relative spatial probability of opening of vents based on geological, geophysical and geochemical data. We evaluated the volcanic risk caused by pyroclastic flows based on the formula proposed by UNESCO (1972), R=H×V×Va, where H is the hazard, V is the vulnerability and Va is the value of the elements at risk. The product H×V was obtained by performing simulations of type eruptions centered in each cell of the grid. The simulation is based on the energy cone scheme proposed by Sheridan and Malin [J. Volcanol. Geotherm. Res. 17 (1983) 187–202], hypothesizing a column collapse height of 100 m for eruptions of VEI=3 and 300 m for eruptions of VEI=4 with a slope angle of 6°. Each simulation has been given the relative probability value associated with the corresponding cell. We made use of the GIS software ArcView 3.2 to evaluate the intersection between the energy cone and the topography. The superposition of the areas invaded by pyroclastic flows (124 simulations for VEI=3 and 37 for VEI=4) was used to obtain the relative hazard map of the area. The relative volcanic risk map is obtained by superimposing the urbanization maps.     

Three types of pyroclastic currents and their deposits: examples from the Vulsini Volcanoes, Italy

This paper deals with ground-hugging, gas–pyroclast currents from explosive volcanic eruptions and their deposits. Key field observations and laboratory determinations are proposed to relate specific deposit types with flow regimes and particle concentration in the transport and depositional systems. Three relevant flow scenarios and corresponding deposit types have been recognized from a survey of pyroclastic successions of the Vulsini Volcanic District (central Italy): (1) dilute, turbulent, pyroclastic currents producing normally or multiply graded beds by direct suspension sedimentation; (2) concentrated bedload regions beneath suspension currents, depositing inversely graded beds by traction carpet sedimentation; (3) self-sustained, high particle concentration, laminar, mass flows developing massive, poorly sorted bodies, with opposite grading of coarse lithic and pumice clasts, overlying fine-grained, inversely graded, basal layers. Main distinguishing criteria include the occurrence and pattern of clast grading, clast–thickness relationships, grain size, ash matrix componentry and pyroclast size–density relationships. Downcurrent and temporal transitions among identified flow scenarios are likely to occur for changing energy conditions and gas–pyroclast ratio both on regional and local scales. The nature and efficiency of magma fragmentation, volatile content, conduit geometry (which determine the characteristics of the erupted mixture and possible lateral blast component at the vent), and the angle of incidence of the column collapse, are suggested as the main factors controlling the generation of one type over the other at flow inception. Dilute, fine-grained, overpressured eruption clouds are thought to favor the formation of low particle concentration turbulent currents. Column collapse over slightly inclined volcano slopes, causing a high degree of compression of the collapsing mixture and of gas expulsion, would favor the generation of high particle concentration pyroclastic currents.     

Influence of fluctuating supply on the emplacement dynamics of channelized lava flows

The evolution of lava flows emplaced on Mount Etna (Italy) in September 2004 is examined in detail through the analysis of morphometric measurements of flow units. The growth of the main channelized flow is consistent with a layering of lava blankets, which maintains the initial geometry of the channel (although levees are widened and raised), and is here explicitly related to the repeated overflow of lava pulses. A simple analytical model is introduced describing the evolution of the flow level in a channelized flow unit fed by a fluctuating supply. The model, named FLOWPULSE, shows that a fluctuation in the velocity of lava extrusion at the vent triggers the formation of pulses, which become increasingly high the farther they are from the vent, and are invariably destined to overflow within a given distance. The FLOWPULSE simulations are in accordance with the observed morphology, characterized by a very flat initial profile followed by a massive increase in flow unit cross-section area between 600 and 700 m downflow. The modeled emplacement dynamics provides also an explanation for the observed substantial “loss” of the original flowing mass with increasing distance from the vent.