Reconciling pyroclastic flow and surge: the multiphase physics of pyroclastic density currents

Two end-member types of pyroclastic density current are commonly recognized: pyroclastic surges are dilute currents in which particles are carried in turbulent suspension and pyroclastic flows are highly concentrated flows. We provide scaling relations that unify these end-members and derive a segregation mechanism into basal concentrated flow and overriding dilute cloud based on the Stokes number (S_T), the stability factor (Σ_T) and the dense-dilute condition (D_D). We recognize five types of particle behaviors within a fluid eddy as a function of S_T and Σ_T: (1) particles sediment from the eddy, (2) particles are preferentially settled out during the downward motion of the eddy, but can be carried during its upward motion, (3) particles concentrate on the periphery of the eddy, (4) particles settling can be delayed or ‘fast-tracked’ as a function of the eddy spatial distribution, and (5) particles remain homogeneously distributed within the eddy. We extend these concepts to a fully turbulent flow by using a prototype of kinetic energy distribution within a full eddy spectrum and demonstrate that the presence of different particle sizes leads to the density stratification of the current. This stratification may favor particle interactions in the basal part of the flow and D_D determines whether the flow is dense or dilute. Using only intrinsic characteristics of the current, our model explains the discontinuous features between pyroclastic flows and surges while conserving the concept of a continuous spectrum of density currents.     

Rhyolitic fallout and pyroclastic density current deposits from a phreatoplinian eruption in the eastern Aegean Sea, Greece

The Kos Plateau Tuff consists of pyroclastic deposits from a major Quaternary explosive rhyolitic eruption, centred about 10 km south of the island of Kos in the eastern Aegean, Greece. Five main units are present, the first two (units A and B) were the product of a phreatoplinian eruption. The eruption style then changed to `dry' explosive style as the eruption intensity increased forming a sequence of ignimbrites and initiating caldera collapse. The final waning phase returned to phreatomagmatic eruptive conditions (unit F). The phreatomagmatic units are fine grained, poorly sorted, and dominated by blocky vitric ash, thickly ash-coated lapilli and accretionary lapilli. They are non-welded and were probably deposited at temperatures below 100°C. All existing exposures occur at distances between 10 km and 40 km from the inferred source. Unit A is a widespread (>42 km from source), thin (upwind on Kos) to very thick (downwind), internally laminated, dominantly ash bed with mantling, sheet-like form. Upwind unit A and the lower and middle part of downwind unit A are ash-rich (ash-rich facies) whereas the upper part of downwind unit A includes thin beds of well sorted fine pumice lapilli (pumice-rich facies). Unit A is interpreted to be a phreatoplinian fall deposit. Although locally the bedforms were influenced by wind, surface water and topography. The nature and position of the pumice-rich facies suggests that the eruption style alternated between `wet' phreatoplinian and `dry' plinian during the final stages of unit A deposition.Unit B is exposed 10–19 km north of the inferred source on Kos, overlying unit A. It is a thick to very thick, internally stratified bed, dominated by ash-coated, medium and fine pumice lapilli in an ash matrix. Unit B shows a decrease in thickness and grain size and variations in bedforms downcurrent that allow definition of several different facies and laterally equivalent facies associations. Unit B ranges from being very thick, coarse and massive or wavy bedded in the closest outcrops to source, to being partly massive and partly diffusely stratified or cross-bedded in medial locations. Pinch and swell, clast-supported pumice layers are also present in medial locations. In the most distal sections, unit B is stratified or massive, and thinner and finer grained than elsewhere and dominated by thickly armoured lapilli. Unit B is interpreted to have been deposited from an unsteady, density stratified, pyroclastic density current which decelerated and progressively decreased its particle load with distance from source. Condensation of steam during outflow of the current promoted the early deposition of ash and resulted in the coarser pyroclasts being thickly ash-coated. The distribution, texture and stratigraphic position of unit B suggest that the pyroclastic density current was generated from collapse of the phreatoplinian column following a period of fluctuating discharge when the eruptive activity alternated between `wet' and `dry'. The pyroclastic density current was transitional in particle concentration between a dilute pyroclastic surge and a high particle concentration pyroclastic flow. Unidirectional bedforms in unit B suggest that the depositional boundary was commonly turbulent and in this respect did not resemble conventional pyroclastic flows. However, unit B is relatively thick and poorly sorted, and was deposited more than 19 km from source, implying that the current comprised a relatively high particle concentration and in this respect, did not resemble a typical pyroclastic surge.     

Vent area and depositional mechanisms of the Upper Member of the Neapolitan Yellow Tuff (Campi Flegrei, Italy): new insights from directional fabric through image analysis

In order to provide new information about the source area and depositional mechanisms of the Upper Member of the Neapolitan Yellow Tuff (NYT), a prominent pyroclastic deposit of the Campi Flegrei Volcanic District (southern Italy), statistics on directional fabric, by means of computer-assisted image analysis on 32 rock samples, were compiled. Seventeen samples were collected along vertical direction on two selected exposures and fifteen were taken from outcrops widely distributed all around the Campi Flegrei Volcanic District. Fabric measurements within the investigated successions reveal a vertically homogeneous direction of the mean particle iso-orientation, with considerable variability in the strength of particle iso-orientation even at cm-scale. The existence of particle iso-orientation can be related to continuous sedimentation from a concentrated bedload region beneath suspension currents, producing massive or inversely graded beds by traction carpet sedimentation. The considerable vertical variability in the strength of iso-orientation is the result of very unstable flow regimes, up to the extreme condition of discrete depositional events, with a variable combination of traction carpet and/or direct suspension sedimentation. The vertical homogeneity in the mean orientation values, found in the investigated sections, may derive from the sequential deposition of laminae to thin beds, whose relatively flat upper surfaces were unable to significantly deflect the depositional system of the following currents. According to the observed homogeneous mean particle orientation values along the investigated vertical profiles, samples collected through areal distribution are considered representative of the local paleo-flow directions of the whole deposit. The mean directions of the samples collected areally show two different coherent patterns which point to the existence of two different source areas. The first, which includes all samples from the northern outcrops, appears to converge in a narrow area about 2 km NE of the town of Pozzuoli, largely in coincidence with the inferred area on the basis of the pumice fall distribution. The second, which includes samples from Capo Miseno and Posillipo areas, points to the central part of the Pozzuoli Bay, about 4 km offshore the town of Pozzuoli.     

The 26.5 ka Oruanui eruption, New Zealand: an introduction and overview

The 26.5 ka Oruanui eruption, from Taupo volcano in the central North Island of New Zealand, is the largest known ‘wet’ eruption, generating 430 km³ of fall deposits, 320 km³ of pyroclastic density–current (PDC) deposits (mostly ignimbrite) and 420 km³ of primary intracaldera material, equivalent to 530 km³ of magma. Erupted magma is >99% rhyolite and <1% relatively mafic compositions (52.3–63.3% SiO₂). The latter vary in abundance at different stratigraphic levels from 0.1 to 4 wt%, defining three ‘spikes’ that are used to correlate fall and coeval PDC activity. The eruption is divided into 10 phases on the basis of nine mappable fall units and a tenth, poorly preserved but volumetrically dominant fall unit. Fall units 1–9 individually range from 0.8 to 85 km³ and unit 10, by subtraction, is 265 km³; all fall deposits are of wide (plinian) to extremely wide dispersal. Fall deposits show a wide range of depositional states, from dry to water saturated, reflecting varied pyroclast:water ratios. Multiple bedding and normal grading in the fall deposits show the first third of the eruption was very spasmodic; short-lived but intense bursts of activity were separated by time breaks from zero up to several weeks to months. PDC activity occurred throughout the eruption. Both dilute and concentrated currents are inferred to have been present from deposit characteristics, with the latter being volumetrically dominant (>90%). PDC deposits range from mm- to cm-thick ultra-thin veneers enclosed within fall material to >200 m-thick ignimbrite in proximal areas. The farthest travelled (90 km), most energetic PDCs (velocities >100 m s⁻¹) occurred during phase 8, but the most voluminous PDC deposits were emplaced during phase 10. Grain size variations in the PDC deposits are complex, with changes seen vertically in thick, proximal accumulations being greater than those seen laterally from near-source to most-distal deposits. Modern Lake Taupo partly infills the caldera generated during this eruption; a 140 km² structural collapse area is concealed beneath the lake, while the lake outline reflects coeval peripheral and volcano–tectonic collapse. Early eruption phases saw shifting vent positions; development of the caldera to its maximum extent (indicated by lithic lag breccias) occurred during phase 10. The Oruanui eruption shows many unusual features; its episodic nature, wide range of depositional conditions in fall deposits of very wide dispersal, and complex interplay of fall and PDC activity.     

Voluminous lava-like precursor to a major ash-flow tuff: low-column pyroclastic eruption of the Pagosa Peak Dacite, San Juan volcanic field, Colorado

The Pagosa Peak Dacite is an unusual pyroclastic deposit that immediately predated eruption of the enormous Fish Canyon Tuff (5000 km³) from the La Garita caldera at 28 Ma. The Pagosa Peak Dacite is thick (to 1 km), voluminous (>200 km³), and has a high aspect ratio (1:50) similar to those of silicic lava flows. It contains a high proportion (40–60%) of juvenile clasts (to 3–4 m) emplaced as viscous magma that was less vesiculated than typical pumice. Accidental lithic fragments are absent above the basal 5–10% of the unit. Thick densely welded proximal deposits flowed rheomorphically due to gravitational spreading, despite the very high viscosity of the crystal-rich magma, resulting in a macroscopic appearance similar to flow-layered silicic lava. Although it is a separate depositional unit, the Pagosa Peak Dacite is indistinguishable from the overlying Fish Canyon Tuff in bulk-rock chemistry, phenocryst compositions, and ⁴⁰Ar/³⁹Ar age.The unusual characteristics of this deposit are interpreted as consequences of eruption by low-column pyroclastic fountaining and lateral transport as dense, poorly inflated pyroclastic flows. The inferred eruptive style may be in part related to synchronous disruption of the southern margin of the Fish Canyon magma chamber by block faulting. The Pagosa Peak eruptive sources are apparently buried in the southern La Garita caldera, where northerly extensions of observed syneruptive faults served as fissure vents. Cumulative vent cross-sections were large, leading to relatively low emission velocities for a given discharge rate. Many successive pyroclastic flows accumulated sufficiently rapidly to weld densely as a cooling unit up to 1000 m thick and to retain heat adequately to permit rheomorphic flow. Explosive potential of the magma may have been reduced by degassing during ascent through fissure conduits, leading to fracture-dominated magma fragmentation at low vesicularity. Subsequent collapse of the 75×35 km² La Garita caldera and eruption of the Fish Canyon Tuff were probably triggered by destabilization of the chamber roof as magma was withdrawn during the Pagosa Peak eruption.     

The Sarikavak Tephra, Galatea, north central Turkey: a case study of a Miocene complex plinian eruption deposit

The Sarikavak Tephra from the central Galatean Volcanic Province (Turkey) represents the deposit of a complex multiple phase plinian eruption of Miocene age. The eruptive sequence is subdivided into the Lower-, Middle-, and Upper Sarikavak Tephra (LSKT, MSKT, USKT) which differ in type of deposits, lithology and eruptive mechanisms.The Lower Sarikavak Tephra is characterised by pumice fall deposits with minor interbedded fine-grained ash beds in the lower LSKT-A. Deposits are well stratified and enriched in lithic fragments up to >50 wt% in some layers. The upper LSKT-B is mainly reversely graded pumice fall with minor amounts of lithics. It represents the main plinian phase of the eruption. The LSKT-A and B units are separated from each other by a fine-grained ash fall deposit. The Middle Sarikavak Tephra is predominantly composed of cross-bedded ash-and-pumice surge deposits with minor pumice fall deposits in the lower MSKT-A and major pyroclastic flow deposits in the upper MSKT-B unit. The Upper Sarikavak Tephra shows subaerial laminated surge deposits in USKT-A and subaqueous tephra beds in USKT-B.Isopach maps of the LSKT pumice fall deposits as well as the fine ash at the LSKT-A/B boundary indicate NNE–SSW extending depositional fans with the source area in the western part of the Ovaçik caldera. The MSKT pyroclastic flow and surge deposits form a SW-extending main lobe related to paleotopography where the deposits are thickest.Internal bedding and lithic distribution of the LSKT-A result from intermittent activity due to significant vent wall instabilities. Reductions in eruption power from (partial) plugging of the vent produced fine ash deposits in near-vent locations and subsequent explosive expulsion of wall rock debris was responsible for the high lithic contents of the lapilli fall deposits. A period of vent closure promoted fine ash fall deposition at the end of LSKT-A. The subsequent main plinian phase of the LSKT-B evolved from stable vent conditions after some initial gravitational column collapses during the early ascent of the re-established eruption plume. The ash-and-pumice surges of the MSKT-A are interpreted as deposits from phreatomagmatic activity prior to the main pyroclastic flow formation of the MSKT-B.     

Fluidization and de-aeration of pyroclastic mixtures: The influence of fines content,polydispersity and shear flow

Fluidization of pyroclastic solids has long been indicated as one key to explain the enhanced mobility of dense pyroclastic gravity currents and their associated hazard. However there is a lack of characterization of the actual pattern and extent of fluidization establishing in real pyroclastic flows and some authors still raise arguments about the relevance of fluidization to the mobility of dense pyroclastic gravity currents. The present paper addresses the fluidization of pyroclastic granular solids with a specific focus on the analysis of factors that may promote homogeneous fluidization and retard solids de-aeration and consolidation. These factors include fines content, particle polydispersity and the establishment of shear flow.     

Toward a revised hazard assessment at Merapi volcano, Central Java

Of 1.1 million people living on the flanks of the active Merapi volcano, 440,000 are at relatively high risk in areas prone to pyroclastic flows, surges, and lahars. For the last two centuries, the activity of Merapi has alternated regularly between long periods of viscous lava dome extrusion, and brief explosive episodes at 8–15 year intervals, which generated dome-collapse pyroclastic flows and destroyed part of the pre-existing domes. Violent explosive episodes on an average recurrence of 26–54 years have generated pyroclastic flows, surges, tephra-falls, and subsequent lahars. The 61 reported eruptions since the mid-1500s killed about 7000 people. The current hazard-zone map of Merapi (Pardyanto et al., 1978) portrays three areas, termed ‘forbidden zone’, ‘first danger zone’ and ‘second danger zone’, based on successively declining hazards. Revision of the hazard map is desirable, because it lacks details necessary to outline hazard zones with accuracy, in particular the valleys likely to be swept by lahars, and excludes some areas likely to be devastated by pyroclastic gravity-currents such as the 22 November 1994 surge. In addition, risk maps should be developed to incorporate social, technical, and economic factors of vulnerability.Eruptive hazard assessment at Merapi is based on reconstructed eruptive history, on eruptive behavior and scenarios, and on existing models and preliminary numerical modeling. Firstly, the reconstructed eruptive activity, in particular for the past 7000 years and from historical accounts of eruptions, helps to define the extent and recurrence frequency of the most hazardous phenomena (Newhall et al., 2000; Camus et al., 2000). Pyroclastic flows traveled as far as 9–15 km from the source, pyroclastic surges swept the flanks as far as 9–20 km away from the vent, thick tephra fall buried temples in the vicinity of Yogyakarta 25 km to the south, and subsequent lahars spilled down the radial valleys as far as 30 km to the west and south. At least one large edifice collapse has occurred in the past 7000 years (Newhall et al., 2000; Camus et al., 2000). Secondly, four eruption scenarios are portrayed as hazardous zones on two maps and derived from the past eruptive behavior of Merapi and from the most affected areas in the past. Thirdly, simple numerical simulation, based on a Digital Elevation Model, a stereo-pair of SPOT satellite images, and one 2D-orthoimage helps to simulate pyroclastic and lahar flowage on the flanks and in radial valley channels, and to outline areas likely to be devastated.Three major threats are identified: (1) a collapse of the summit dome in the short-to mid-term, that can release large-volume pyroclastic flows and high-energy surges towards the south–southwest sector of the volcano; (2) an explosive eruption, much larger than any since 1930, may sweep all the flanks of Merapi at least once every century; (3) a potential collapse of the summit area, involving the fumarolic field of Gendol and part of the southern flank, which can contribute to moderate-scale debris avalanches and debris flows.     

The Holocene Secche di Lazzaro phreatomagmatic succession (Stromboli,Italy): evidence of pyroclastic density current origin deduced by facies analysis and AMS flow directions

The edifice of Stromboli volcano gravitationally collapsed several times during its volcanic history (>100 ka–present). The largest Holocene event occurred during the final stage of the Neostromboli activity (∼13–5 ka), and was accompanied by the emplacement of phreatomagmatic and lahar deposits, known as the Secche di Lazzaro succession. A stratigraphic and paleomagnetic study of the Secche di Lazzaro deposits allows the interpretation of the emplacement and the eruptive processes. We identify three main units within the succession that correspond to changing eruption conditions. The lower unit (UA) consists of accretionary lapilli-rich, thinly bedded, parallel- to cross-stratified ash deposits, interpreted to indicate the early stages of the eruption and emplacement of dilute pyroclastic density currents. Upward, the second unit (UB) of the deposit is more massive and the beds thicker, indicating an increase in the sedimentation rate from pyroclastic density currents. The upper unit (UC) caps the succession with thick, immediately post-eruptive lahars, which reworked ash deposited on the volcano’s slope. Flow directions obtained by Anisotropy of Magnetic Susceptibility (AMS) analysis of the basal bed of UA at the type locality suggest a provenance of pyroclastic currents from the sea. This is interpreted to be related to the initial base-surges associated with water–magma interaction that occurred immediately after the lateral collapse, which wrapped around the shoulder of the sector collapse scar. Upward in the stratigraphy (upper beds of UA and UB) paleoflow directions change and show a provenance from the summit vent, probably related to the multiple collapses of a vertical, pulsatory eruptive column.     

Origin of reverse-graded bedding in air-fall pumice, Coso Range, California

The origin of reverse grading in air-fall pyroclastic deposits has been ascribed to: (1) changing conditions at an erupting vent; (2) deposition in water; or (3) rolling of large clasts over smaller clasts on the surface of a steep slope. Structural features in a deposit of air-fall pumice lapilli in the Coso Range, California, indicate that reverse grading there formed by a fourth mechanism during flow of pumice. Reverse-graded beds in this deposit occur where pumice lapilli fell on slopes at or near the angle of repose and formed as parts of the blanket of accumulating pumice became unstable and flowed downslope. The process of size sorting during such flow is probably analogous to that which sorts sand grains in a reverse fashion during avalanching on the slip faces of sand dunes, attributed by Bagnold (1954a) to a grain-dispersive pressure acting on particles subjected to a shear stress. In view of the several ways in which air-fall pyroclastic debris may become reverse graded, caution is advised in interpretation of the origin of this structure both in modern and in ancient deposits.     

The 1984 nuée-ardente deposits of Merapi volcano,Central Java,Indonesia: stratigraphy,textural characteristics,and transport mechanisms

For many centuries Merapi volcano has generated hot avalanches of blocks, lapilli and ashes, derived from the destruction of partially solidified, viscous lava domes (Merapi-type nuées ardentes). On 15 June 1984, at least four nuées ardentes came down the southwest slope of the Merapi, the first and the last being responsible for more than 99% of the deposits which are now exposed. The first nuée ardente, a Merapi-type nuée ardente, was produced by the destruction of the dome, travelled 7 km from the crater, leaving a measured deposit, 2.7 m thick, 4 km from the crater, near its upper depositional limit, regularly increasing to a maximum measured thickness of 12 m at the front of the deposit. The lower contact is sharp, non-erosive, with pines still rooted in the underlying paleosol. The deposit consists of 50% ash, 33% lapilli, and 17% blocks, with two subpopulations (one Rosin and one normal), and is finespoor, with less than 4% of fine ash (d finer than 4 ). The deposit displays reverse population grading of both vesiculated and massive clasts, and of the maximum grain size. The maximum size significantly increases regularly down-current over most of the exposed length of unit 1, and bed thickness increases for the entire length of the deposit. The deposit of the second nuée ardente is only 6–21 cm thick, and of very limited lateral extent. It is a normally graded, coarse to fine ash, with a finespoor base. The third unit consists of fines-poor, normally graded coarse ash, exposed in low-amplitude (20–40 cm), 12-m-wavelength dunes. The deposit of the fourth nuée ardente rests in sharp erosive contact on the underlying unit, increasing in thickness down-flow. It consists of transitional coarse and fine-grained strata, 6–130 c cm thick, dipping 5–10° down-flow. The deposit, made up of two subpopulations (one Rosin and one normal), is normally graded over the entire bed, but coarsegrained strata are reversely graded. The relative content of vesiculated clasts increases up-bed in both strata types, from 12% at the base to 40% at the top. The characteristics of unit 1 suggest that it accumulated from a concentrated suspension of cohesionless solids exhibiting non-Newtonian behavior, where dispersive pressure played an important role in the suspension of the clasts. Units 2 and 3 were probably deposited from dilute turbulent suspensions, whereas the upper unit (4) is a classic example of deposition from a high-density turbulent suspension leading to the formation of multiple traction carpets driven by the overlying, lower-density, surge. The horizontal distance travelled by a hot rock avalanche may be influenced by its transport mechanism. Debris flows are mobile on very low slopes-as low as 1°-whereas grain flows, even density-modified grain flows, require relatively high slopes-more than 6° at Merapi-to maintain their mobility. If the present Merapi dome were to collapse and produce a debris flow, its present volume coupled with the minimal 1.5 km vertical drop could travel a distance ranging between 15 and 30 km. However, if transport were by grain flow mechanisms, the mass could come to rest as it reaches a 5–10° slope.     

Interrelations among pyroclastic surge,pyroclastic flow,and lahars in Smith Creek valley during first minutes of 18 May 1980 eruption of Mount St. Helens,USA

A devastating pyroclastic surge and resultant lahars at Mount St. Helens on 18 May 1980 produced several catastrophic flowages into tributaries on the northeast volcano flank. The tributaries channeled the flows to Smith Creek valley, which lies within the area devastated by the surge but was unaffected by the great debris avalanche on the north flank. Stratigraphy shows that the pyroclastic surge preceded the lahars; there is no notable “wet” character to the surge deposits. Therefore the lahars must have originated as snowmelt, not as ejected water-saturated debris that segregated from the pyroclastic surge as has been inferred for other flanks of the volcano. In stratigraphic order the Smith Creek valley-floor materials comprise (1) a complex valley-bottom facies of the pyroclastic surge and a related pyroclastic flow, (2) an unusual hummocky diamict caused by complex mixing of lahars with the dry pyroclastic debris, and (3) deposits of secondary pyroclastic flows. These units are capped by silt containing accretionary lapilli, which began falling from a rapidly expanding mushroom-shaped cloud 20 minutes after the eruption's onset. The Smith Creek valley-bottom pyroclastic facies consists of (a) a weakly graded basal bed of fines-poor granular sand, the deposit of a low-concentration lithic pyroclastic surge, and (b) a bed of very poorly sorted pebble to cobble gravel inversely graded near its base, the deposit of a high-concentration lithic pyroclastic flow. The surge apparently segregated while crossing the steep headwater tributaries of Smith Creek; large fragments that settled from the turbulent surge formed a dense pyroclastic flow along the valley floor that lagged behind the front of the overland surge. The unusual hummocky diamict as thick as 15 m contains large lithic clasts supported by a tough, brown muddy sand matrix like that of lahar deposits upvalley. This unit contains irregular friable lenses and pods meters in diameter, blocks incorporated from the underlying dry and hot pyroclastic material that had been deposited only moments earlier. The hummocky unit is the deposit of a high-viscosity debris flow which formed when lahars mingled with the pyroclastic materials on Smith Creek valley floor. Overlying the debris flow are voluminous pyroclastic deposits of pebbly sand cut by fines-poor gas-escape pipes and containing charred wood. The deposits are thickest in topographic lows along margins of the hummocky diamict. Emplaced several minutes after the hot surge had passed, this is the deposit of numerous secondary pyroclastic flows derived from surge material deposited unstably on steep valley sides.     

Pyroclastic flow deposits of the 1991 eruption of Volcán de Colima,Mexico

The April 16, 1991, eruption of Volcán de Colima represents a classical example of partial dome collapse with the generation of progressively longer-runout, Merapi-type pyroclastic flows that traveled up to 4 km along the El Cordoban gullies (East, Central and West). The flows filled the gullies with block-and-ash flow deposits up to 10 m thick, of which, after 7 years of erosion, only remnants remained in the El Cordoban West and East gullies. The El Cordoban Central gully, however, provided a well-preserved and incised longitudinal section of the 1991 deposits. The deposits were emplaced as proximal and distal facies, separated by a change in slope angle from >30° to <20°. The proximal facies consists of massive, clast-supported flow units (up to 1 m thick) with andesite blocks locally supported by a matrix of coarse ash and devoid of segregation structures or grading. The distal facies consists of a massive, matrix-supported deposit up to 8 m thick, which contains dispersed andesite blocks in a fine ash matrix. In the distal facies, a train of blocks marks flow-unit upper boundaries and, although sorting is poor, some grading is present. Thin, finely stratified, or dune-bedded layers of fine ash material are locally present above or below units of both facies. Sedimentologic parameters show that the size or fraction of large pyroclasts (larger than –1 ) decreases from proximal to distal facies, as the percentage of matrix (0 to 4 ) increases, especially immediately beyond the break in slope. We propose that the propagation of the Colima pyroclastic flows is critically dependent on local slope angle, the presence of erodible slope debris, and the decrease in grain size with distance from the vent. The progressive fining is probably caused by some combination of erosion, clast breakup and deposition of larger pyroclasts, and is itself influenced by the slope angle. In the proximal region, the flows moved as granular avalanches, in which interacting grains ground each other and erosion occurred to produce an overriding dilute ash cloud. The maximum runout distance of the avalanches was controlled by the angle of repose of the material, and the volume and grain size of source and eroded material. Because the slope angle is close to the repose angle for this debris, granular avalanches were not able to propagate far beyond the change in slope. If, however, an avalanche had enough mass in finer grain size fractions, at least part of the flow continued beyond the break in slope and across the volcano apron, propagating in a turbulent state and depositing surge layers, or in an otherwise settling-modified state and depositing block-and-ash flow layers.Editorial responsibility: T Druitt     

Evidence on the Koseda coast of Yakushima Island of a tsunami during the 7.3 ka Kikai caldera eruption

Previous research indicates that Yakushima Island, southwestern Japan, may have been struck by a huge tsunami before or soon after the arrival of the Koya pyroclastic flow during the 7.3 ka caldera‐forming Kikai eruption, but this has not yet been confirmed. This paper describes sedimentological and chronostratigraphic evidence showing that Unit TG, one of three gravel beds exposed on the Koseda coast of northeast Yakushima Island and investigated here, is a tsunami deposit. Unit TG is a poorly sorted, 30 cm thick gravel bed overlying a wave‐cut bench and underlying a Koya pyroclastic flow deposit. Sparse wood fragments in Unit TG were dated at 7 416–7 167 cal year BP. The constituent gravel clasts of Unit TG are similar in composition to those of modern beach and river deposits along the Koseda coast. Unit TG also contains pumice clasts whose chemistry is identical to that of pumice derived from the 7.3 ka eruption at Kikai caldera. The long‐axis orientations and composition of gravel clasts in Unit TG suggest that they were transported by a landward‐travelling high‐particle‐concentration flow, which suggests that Unit TG was deposited by a tsunami run‐up flow during the 7.3 ka Kikai caldera eruption, just before the arrival of the major Koya pyroclastic flow at the Koseda coast. Whether the 7.3 ka tsunami was caused by a volcanic eruption or an earthquake remains unclear, but Unit TG demonstrates that a tsunami arrived immediately before emplacement of a Koya pyroclastic flow.     

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     

Pyroclastic current dynamic pressure from aerodynamics of tree or pole blow-down

The common occurrence of tree and pole blow-down from pyroclastic currents provides an opportunity to estimate properties of the currents. Blow-down may occur by uprooting (root zone rupture), or flexure or shear at some point on the object. If trees are delimbed before blow-down, each tree or pole can be simulated by a cylinder perpendicular to the current. The force acting on a cylinder is a function of flow dynamic pressure, cylinder geometry, and drag coefficient. Treated as a cantilever of circular cross-section, the strength for the appropriate failure mode (rupture, uprooting or flexure) can then be used to estimate the minimum necessary current dynamic pressure. In some cases, larger or stronger standing objects can provide upper bounds on the dynamic pressure. This analysis was treated in two ways: (1) assuming that the current properties are vertically constant; and (2) allowing current velocity and density to vary vertically according to established models for turbulent boundary layers and stratified flow. The two methods produced similar results for dynamic pressure. The second, along with a method to approximate average whole-current density, offers a means to estimate average velocity and density over the height of the failed objects. The method is applied to several example cases, including Unzen, Mount St. Helens, Lamington, and Merapi volcanoes. Our results compare reasonably well with independent estimates. For several cases, we found that it is possible to use the dynamic pressure equations developed for vertically uniform flow, along with the average cloud density multiplied by a factor of 2–5, to determine average velocity over the height of the failed object.     

The Ottaviano eruption of Somma-Vesuvio (8000 y B.P.): a magmatic alternating fall and flow-forming eruption

The Ottaviano eruption occurred in the late neolithic (8000 y B.P.). 2.40 km³ of phonolitic pyroclastic material (0.61 km³ DRE) were emplaced as pyroclastic flow, surge and fall deposits. The eruption began with a fall phase, with a model column height of 14 km, producing a pumice fall deposit (LA). This phase ended with short-lived weak explosive activity, giving rise to a fine-grained deposit (L1), passing to pumice fall deposits as the result of an increasing column height and mass discharge rate. The subsequent two fall phases (producing LB and LC deposits), had model column heights of 20 and 22 km with eruption rates of 2.5 × 10⁷ and 2.81 × 10⁷ kg/s, respectively. These phases ended with the deposition of ash layers (L2 and L3), related to a decreasing, pulsing explosive activity. The values of dynamic parameters calculated for the eruption classify it as a sub-plinian event. Each fall phase was characterized by variations in the eruptive intensity, and several pyroclastic flows were emplaced (F1 to F3). Alternating pumice and ash fall beds record the waning of the eruption. Finally, owing to the collapse of a eruptive column of low gas content, the last pyroclastic flow (F4) was emplaced.     

An ignimbrite veneer deposit: The trail-marker of a pyroclastic flow

The term “ignimbrite veneer deposit” (IVD) is proposed for a new kind of pyroclastic deposit which is found associated with, and passes laterally into, Taupo ignimbrite of valley pond type in New Zealand. It forms a thin layer mantling the landscape over 15,000 km², and is regarded as the deposit from the trailing “tail” of a pyroclastic flow, where a relaxation of shear stress favoured the deposition of the basal part of the flow. The IVD differs little in grain-size from the associated ignimbrite, but it shows a crude internal stratification attributed to the deposition of a succession of layers, one after the passage of each pulse of the pyroclastic flow. It locally contains laterally-discontinuous lenses of coarse pumice (“lee-side lenses”) on the far-vent side of topographic obstacles. In nearvent exposures the Taupo IVD shows lensoid and cross-stratified bed-forms even where it stands on a planar surface, attributed to deposition from a flow travelling at an exceedingly high velocity.An IVD can be distinguished from a poorly sorted pyroclastic fall deposit because the beds in it show more rapid lateral variations in thickness, it may show a low-angle cross-stratification, and it contains carbonised wood from trees not in the position of growth; from the deposit of a wet base surge because it lacks vesicles and strong antidune-like structures and contains carbonised vegetation, and from a hot and dry pyroclastic surge deposit because it possesses a high content of pumice and “fines”.The significance of an IVD is that it records the passage of a pyroclastic flow, where the flow itself has moved farther on.     

The 1982 eruptions of El Chichon volcano,Mexico (3): Physical properties of pyroclastic surges

Two major pyroclastic surges generated during the 4 April 1982 eruption of El Chichon devastated an area of 153 km² with a quasi-radial distribution around the volcano. The hot surge clouds carbonized wood throughout their extent and were too hot to allow accretionary lapilli formation by vapor condensation. Field evidence indicates voidage fraction of 0.99 in the surge cloud with extensive entrainment of air. Thermal calculations indicate that heat content of pyroclasts can heat entrained air and maintain high temperatures in the surge cloud. The dominant bed form of the surge deposits are sand waves shaped in dune forms with vertical form index of 10–20, characterized by stoss-side erosion and lee-side deposition of 1–10 cm reversely graded laminae. A systematic decrease in maximum lithic diameter with distance from source is accompanied by decrease in wavelength and amplitude. Modal analysis indicates fractionation of glass and pumice from the surge cloud relative to crystals, resulting in loss of at least 10%–25% of the cloud mass due to winnowing out of fines during surge emplacement. Greatest fractionation from the –1.0–0.0– grain sizes reflects relatively lower pumice particle density in this range and segregation in the formative stages of the surge cloud. Extensive pumice rounding indicates abrasion during bed-load transport. Flow of pyroclastic debris in the turbulent surge cloud was by combination of bed-load and suspended-load transport. The surges are viewed as expanding pyroclastic gravity flows, which entrain and mix with air during transport. The balance between sedimentation at the base of the surge cloud and expansion due to entrainment of air contributed to low cloud density and internal turbulence, which persisted to the distal edge of the surge zone.     

Nuées ardentes of 22 November 1994 at Merapi volcano, Java, Indonesia

Nuées ardentes associated with dome collapse on 22 November 1994, at Merapi volcano traveled to the south–southwest as far as 6.5 km, and collectively accumulated roughly 2.5–3 million cubic meters of deposits. The damaged area comprises 9.5 km² and is covered by two nuée ardente facies, a conventional “Merapi-type”, valley-fill block-and-ash flow facies and a pyroclastic surge facies. The proximal deposits reflect the accumulation of dozens of nuées ardentes, with many subsidiary flow units. The distal deposits are more simply organized, as only a few individual events reached to distances >3.5 km. The stratigraphic relationships north of Turgo hill indicate that the surge deposits are a facies of particularly mobile nuées ardentes that also deposited channeled block-and-ash flow facies. They further suggest that the surge facies beyond the channel margins correlate laterally with a finer-grained sublayer locally developed at the base of the block-and-ash flow facies. Eyewitness reports suggest that the emplacement of the block-and-ash flow facies in the distal part of the Boyong river may have followed, by a short time interval, the destruction and deposition of the surge facies at Turgo village. The stratigraphy is in accord with the eyewitness reports. The surge facies was emplaced by a dilute surge current, detached from the same dome-collapse nuée ardente that, as a separate flow unit, subsequently emplaced the distal block-and-ash deposit in the Boyong valley. The detachment occurred at higher elevations, likely at or above the slope break at about 2000 m elevation. This flow separation enabled the surge current to shortcut over the landscape and to emplace its deposit even as the block-and-ash flow continued its tortuous southward movement in the Boyong channel. Dome-collapse nuée ardente activity formed the bulk of the eruption, which was accompanied by virtually no significant vertical summit explosive activity.