Pyroclastic flows from the 1991 eruption of Unzen volcano,Japan

Pyroclastic flows from Unzen were generated by gravitational collapse of the growing lava dome. As soon as the parental lobe failed at the edge of the dome, spontaneous shattering of lava occurred and induced a gravity flow of blocks and finer debris. The flows had a overhanging, tongue-like head and cone- or rollershaped vortices expanding outward and upward. Most of the flows traveled from 1 to 3 km, but some flows reached more than 4 km, burning houses and killing people in the evacuated zone of Kita-kamikoba on the eastern foot of the volcano. The velocities of the flows ranged from 15 to 25 m/s on the gentle middle flank. Observations of the flows and their deposits suggest that they consisted of a dense basal avalanche and an overlying turbulent ash cloud. The basal avalanche swept down a topographic low and formed to tongue-like lobe having well-defined levees; it is presumed to have moved as a non-Newtonian fluid. The measured velocities and runout distances of the flows can be matched to a Bingham model for the basal avalanche by the addition of turbulent resistance. The rheologic model parameters for the 29 May flow are as follows: the density is 1300 kg/m³, the yield strength is 850 Pa, the viscosity is 90 Pa s, and the thickness of the avalanche is 2 m. The ash cloud is interpreted as a turbulent mixing layer above the basal avalanche. The buoyant portions of the cloud produced ash-fall deposits, whereas the dense portions moved as a surge separated from the parental avalanche. The ash-cloud surges formed a wide devastated zone covered by very thin debris. The initial velocities of the 3 June surges, when they detached from avalanches, are determined by the runout distance and the angle of the energy-line slope. A comparison between the estimated velocities of the 3 June avalanches and the surges indicates that the surges that extended steep slopes along the avalanche path, detached directly from the turbulent heads of the avalanches. The over-running surge that reached Kita-Kamikoba had an estimated velocity higher than that of the avalanche; this farther-travelled surge is presumed to have been generated by collapse of a rising ash-cloud plume.     

Pyroclastic flows and surges over water: an example from the 1883 Krakatau eruption

Pyroclastic deposits from the 1883 eruption of Krakatau are described from areas northeast of the volcano on the islands of Sebesi, Sebuku, and Lagoendi, and the southeast coast of Sumatra. Massive and poorly stratified units formed predominantly from pyroclastic flows and surges that traveled over the sea for distances up to 80 km. Granulometric and lithologic characteristics of the deposits indicate that they represent the complement of proximal subaerial and submarine pyroclastic flow deposits laid down on and close to the Krakatau islands. The distal deposits exhibit a decrease in sorting coefficient, median grain size, and thickness with increasing distance from Krakatau. Crystal fractionation is consistent with the distal facies being derived from the upper part of gravitationally segregated pyroclastic flows in which the relative amount of crystal enrichment and abundance of dense lithic clasts diminished upwards. The deposits are correlated to a major pyroclastic flow phase that occurred on the morning of 27 August at approximately 10 a.m. Energetic flows spread out away from the volcano at speeds in excess of 100 km/h and traveled up to 80 km from source. The flows retained temperatures high enough to burn victims on the SW coast of Sumatra. Historical accounts from ships in the Sunda Straits constrain the area affected by the flows to a minimum of 4x10³ km². At the distal edge of this area the flows were relatively dilute and turbulent, yet carried enough material to deposit several tens of centimeters of tephra. The great mobility of the Krakatau flows from the 10 a.m. activity may be the result of enhanced runout over the sea. It is proposed that the generation of steam at the flow/seawater interface may have led to a reduction in the sedimentation of particles and consequently a delay in the time before the flows ceased lateral motion and became buoyantly convective. The buoyant distal edge of these ash-and steam-laden clouds lifted off into the atmosphere, leading to cooling, condensation, and mud rain.     

Pyroclastic surges and flows from the 8–10 May 1997 explosive eruption of Bezymianny volcano, Kamchatka, Russia

The 8-10 May 1997 eruption of Bezymianny volcano began with extrusion of a crystallized plug from the vent in the upper part of the dome. Progressive gravitational collapses of the plug caused decompression of highly crystalline magma in the upper conduit, leading at 13:12 local time on 9 May to a powerful, vertical Vulcanian explosion. The dense pyroclastic mixture collapsed in boil-over style to generate a pyroclastic surge which was focused toward the southeast by the steep-walled, 1956 horseshoe-shaped crater. This surge, with a temperature <200 °C, covered an elliptical area >30 km² with deposits as much as 30 cm thick and extending 7 km from the vent. The surge deposits comprised massive to vaguely laminated, gravelly sand (Md -1.2 to 3.7J sorting 1.2 to 3J) of poorly vesiculated andesite (mean density 1.82 g cm^-3; vesicularity 30 vol%; SiO₂ content ~58.0 wt%). The deposits, with a volume of 5-15᎒⁶ m³, became finer grained and better sorted with distance; the maximal diameter of juvenile clasts decreased from 46 to 4 cm. The transport and deposition of the surge over a snowy landscape generated extensive lahars which traveled >30 km. Immediately following the surge, semi-vesiculated block-and-ash flows were emplaced as far as 4.7 km from the vent. Over time the juvenile lava in clasts of these flows became progressively less crystallized, apparently more silicic (59.0 to 59.9 wt% SiO₂) and more vesiculated (density 1.64 to 1.12 g cm^-3; vesicularity 37 to 57 vol%). At this stage the eruption showed transitional behavior, with mass divided between collapsing fountain and buoyant column. The youngest pumice-and-ash flows were accompanied by a sustained sub-Plinian eruption column ~14 km high, from which platy fallout clasts were deposited (~59.7% SiO₂; density 1.09 g cm^-3; vesicularity 58 vol%). The explosive activity lasted about 37 min and produced a total of ~0.026 km³ dense rock equivalent of magma, with an average discharge of ~1.2᎒⁴ m³ s^-1. A lava flow ~200 m long terminated the eruption. The evolutionary succession of different eruptive styles during the explosive eruption was caused by vertical gradients in crystallization and volatile content of the conduit magma, which produced significant changes in the properties of the erupting mixture.     

长白山天池火山碎屑流灾害区划

在回顾总结了国外火山碎屑流灾害分析模型研究历史的基础上,本文选取了Flow3D模型对我国东北地区长白山天池火山未来大喷发可能产生的火山碎屑流进行了灾害区域划分。以长白山天池火山现代地形为依据,设定了11条未来爆炸式火山喷发时产生的火山碎屑流的可能流动线路。模拟结果表明,在喷发柱高度为10km的情况下,灾害区划最大半径为13.7km;在喷发柱高度为20km的情况下,灾害区划最大半径为35.4km;在喷发柱高度为30km的情况下,灾害区划最大半径为57.8km。在此基础上,得出了长白山天池火山未来发生中规模、大规模和超大规模火山喷发时火山碎屑流的覆盖范围,完成了我国第一幅长白山天池火山碎屑流灾害区划图。     

Pyroclastic deposits from Sarychev Peak,Matua Island discharged in June 2009

We report the first results from a study of pyroclastic material discharged by Sarychev Peak in June 2009 (Matua I., Kuril Is.). We describe the deposits of pyroclastic flows and pyroclastic surges, the ash deposited from the plume of pyroclastic flows, and tephra.     

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.     

Effects of forest harvesting on the occurrence of landslides and debris flows in steep terrain of central Japan

Landslides and debris flows associated with forest harvesting can cause much destruction and the influence of the timing of harvesting on these mass wasting processes therefore needs to be assessed in order to protect aquatic ecosystems and develop improved strategies for disaster prevention. We examined the effects of forest harvesting on the frequency of landslides and debris flows in the Sanko catchment (central Japan) using nine aerial photo periods covering 1964 to 2003. These photographs showed a mosaic of different forest ages attributable to the rotational management in this area since 1912. Geology and slope gradient are rather uniformly distributed in the Sanko catchment, facilitating assessment of forest harvesting effects on mass wasting without complication of other factors. Trends of new landslides and debris flows correspond to changes in slope stability explained by root strength decay and recovery; the direct impact of clearcutting on landslide occurrence was greatest in forest stands that were clearcut 1 to 10 yr earlier with progressively lesser impacts continuing up to 25 yr after harvesting. Sediment supply rate from landslides in forests clearcut 1 to 10 yr earlier was about 10‐fold higher than in control sites. Total landslide volume in forest stands clearcut 0 to 25 yr earlier was 5·8 × 10³ m³ km^?2 compared with 1·3 × 10³ m³ km^?2 in clearcuts >25 yr, indicating a fourfold increase compared with control sites during the period when harvesting affected slope stability. Because landslide scars continue to produce sediment after initial failure, sediment supply from landslides continues for 45 yr in the Sanko catchment. To estimate the effect of forest harvesting and subsequent regeneration on the occurrence of mass wasting in other regions, changes in root strength caused by decay and recovery of roots should be investigated for various species and environmental conditions. Copyright © 2007 John Wiley & Sons, Ltd.     

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     

Cone-building block-and-ash flows: the Senyama volcanic products of O’e Takayama volcano,SW Japan

The Senyama volcanic products of the late Pliocene to early Pleistocene O’e Takayama volcano overlie a 100-m-thick, late Pliocene coastal quartz-sandstone and are intruded by an early Pleistocene dacite dome. The Senyama volcanic products are the remains of a cone that retains a basal part 1.5 km across and 150–250 m high from the substrate. The cone comprises dacite block-and-ash flow deposits and minor base-surge deposits occur at the base. Single beds of the block-and-ash flow deposits are 1–16 m thick and dip inward 20–40° at the base of the cone and inward or outward 10–20° at the summit. Juvenile fragments in the block-and-ash flow deposits are non- to poorly vesicular and commonly have curviplanar surfaces and prismatic joints extending inward from the surfaces, which imply quenching and brittle fracturing of dacite lava. They are variably hydrothermally altered. Nevertheless, juvenile blocks appear to retain a uniform direction of the magnetization vector residual during thermal demagnetization between 280°C and 625°C. At the time of the eruption, the well-sorted sand of the substrate was at the coast and a good aquifer that facilitated explosive interaction of water and the ascending dacite lava. The mechanism of the explosion perhaps involved thermal contraction cracking of the dacite lava, water-inflow into the interior of the lava, and explosive expansion of the water. Initial phreatomagmatic explosions opened the vent. Succeeding phreatomagmatic or phreatomagmatic–vulcanian explosions produced block-and-ash flow deposits around the vent. Hydrothermal silver-ore deposits and manganese-oxide deposits occur in the Senyama volcanic products and the underlying sandstone, respectively. They could represent post-eruptive activity of the hydrothermal system developed in and around the cone.     

Pyroclastic deposits as a guide for reconstructing the multi-stage evolution of the Somma-Vesuvius Caldera

 The evolution of the Somma-Vesuvius caldera has been reconstructed based on geomorphic observations, detailed stratigraphic studies, and the distribution and facies variations of pyroclastic and epiclastic deposits produced by the past 20,000 years of volcanic activity. The present caldera is a multicyclic, nested structure related to the emptying of large, shallow reservoirs during Plinian eruptions. The caldera cuts a stratovolcano whose original summit was at 1600–1900 m elevation, approximately 500 m north of the present crater. Four caldera-forming events have been recognized, each occurring during major Plinian eruptions (18,300 BP "Pomici di Base", 8000 BP "Mercato Pumice", 3400 BP "Avellino Pumice" and AD 79 "Pompeii Pumice"). The timing of each caldera collapse is defined by peculiar "collapse-marking" deposits, characterized by large amounts of lithic clasts from the outer margins of the magma chamber and its apophysis as well as from the shallow volcanic and sedimentary units. In proximal sites the deposits consist of coarse breccias resulting from emplacement of either dense pyroclastic flows (Pomici di Base and Pompeii eruptions) or fall layers (Avellino eruption). During each caldera collapse, the destabilization of the shallow magmatic system induced decompression of hydrothermal–magmatic and hydrothermal fluids hosted in the wall rocks. This process, and the magma–ground water interaction triggered by the fracturing of the thick Mesozoic carbonate basement hosting the aquifer system, strongly enhanced the explosivity of the eruptions. Received: 24 November 1997 / Accepted: 23 March 1999     

Generation of block and ash flows during the 1990–1995 eruption of Unzen Volcano, Japan

Processes generating block and ash flows by gravitational dome collapse (Merapi-type pyroclastic flow) were observed in detail during the 1990–1995 eruption of Unzen volcano, Japan. Two different types were identified by analysis of video records and observations during helicopter flights. Most of the block and ash flows erupted during the 1991–1993 exogenous dome growth stage initially involved crack propagation due to cooling and flowage of the dome lava lobes. The mass around the crack became unstable, locally decreasing in tensile strength. Finally, a slab separated from the lobe front, fragmented progressively from the base to the top within a few seconds, and became a block and ash flow. Rock falls immediately followed, in response to local instability of the lobe front. Clasts in these rock falls fragmented and merged with the preceding flow. In contrast, block and ash flows during the endogenous dome growth stage in 1994 were generated due to local bulge of the dome. Unstable lava blocks collapsed and subsequently fragmented to produce block and ash flows.     

Occurrence conditions of hyperpycnal flows, and their significance for organic-matter sedimentation in a Holocene estuary, Niigata Plain, Central Japan

River floods influence sedimentary environments and ecosystems from the terrestrial to the deep-marine. This study documents the occurrence conditions of hyperpycnal flows generated by river floods and related organic-matter sedimentation for Holocene sediments of the Niigata Plain, Central Japan, based on detailed sedimentary facies, total sulfur and total organic carbon content, diatom assemblages and organic-matter composition. Holocene sediments of the Niigata Plain consist of sand, mud and gravel that were deposited in estuarine and fluvial systems during a sea-level rise (15 000–6800 years BP) and stillstand (after 6800 years BP) following the Last Glacial Maximum. Hyperpycnites are present in the upper part of the estuarine lagoon sediments. The depositional age is considered to be about 5000 years BP. The hyperpycnites comprise two successions of a top fining-up unit and a basal coarsening-up unit, and include abundant terrigenous organic matter and freshwater diatoms. A large volume of freshwater is inferred to have flowed into the lagoon during deposition of the upper part of the lagoon sediments. In consequence, hyperpycnal flows may have readily formed in the lagoon, because the halocline was weak. The hyperpycnal flows also produced a layer of concentrated terrigenous organic matter in the uppermost part of the hyperpycnites. The abundant organic matter on the estuarine floor is inferred to have produced anoxic bottom conditions owing to oxidative decomposition by benthic bacteria.     

Turbulence structures in non-uniform flows

This study investigates turbulence structures in steady and non-uniform flows. Equations of Reynolds shear stress and turbulent velocity fluctuations are derived and their physical interpretations are explained. The theoretical results show that, different from previous studies, the variation of water surface can generate the wall-normal velocity, resulting in deviations of Reynolds shear stress and turbulence intensities from those in uniform flows. A self-similarity relationship is found between the Reynolds shear stress and turbulence intensities in non-uniform flows. The existence of self-similarity indicates that the effect of non-uniformity does not influence the mixing length. An empirical equation has been proposed to express the relationship based on experimental data available in the literature. Good agreement is achieved between the measured and predicted turbulence intensities by applying the self-similarity relationship.     

The 15 September 1991 pyroclastic flows at Unzen Volcano (Japan): a flow model for associated ash-cloud surges

Large-scale collapse of a dacite dome in the late afternoon of 15 September 1991 generated a series of pyroclastic-flow events at Unzen Volcano. Pyroclastic flows with a volume of 1×10⁶ m³ (as DRE) descended the northeastern slope of the volcano, changing their courses to the southeast due to topographic control. After they exited a narrow gorge, an ash-cloud surge rushed straight ahead, detaching the main body of the flow that turned and followed the topographic lows to the east. The surge swept the Kita-Kamikoba area, which had been devastated by the previous pyroclastic-flow events, and transported a car as far as 120 m. Following detachment, the surge lost its force after it moved several hundred meters, but maintained a high temperature. The deposits consist of a bottom layer of better-sorted ash (unit 1), a thick layer of block and ash (unit 2), and a thin top layer of fall-out ash (unit 3). Unit 2 overlies unit 1 with an erosional contact. The upper part of unit 2 grades into better-sorted ash. At distal block-and-ash flow deposits, the bottom part of unit 2 also consists of better-sorted ash, and the contact with the unit 1 deposits becomes ambiguous. Video footage of cascading pyroclastic flows during the 1991–1995 eruption, traveling over surfaces without any topographic barriers, revealed that lobes of ash cloud protruded intermittently from the moving head and sides, and that these lobes surged ahead on the ground surface. This fact, together with the inspection by helicopter shortly after the events, suggests that the protruded lobes consisted of better-sorted ash, and resulted in the deposits of unit 1. The highest ash-cloud plume at the Oshigadani valley exit, and the thickest deposition of fall-out ash over Kita-Kamikoba and Ohnokoba, indicate that abundant ash was also produced when the flow passed through a narrow gorge. In the model presented here, the ash clouds from the pyroclastic flows were composed of a basal turbulent current of high concentration (main body), an overriding and intermediate fluidization zone, and an overlying dilute cloud. Release of pressurized gas in lava block pores, due to collisions among blocks and the resulting upward current, caused a zone of fluidization just above the main body. The mixture of gas and ash sorted in the fluidization zone moved ahead and to the side of the main body as a gravitational current, where the ash was deposited as surge deposits. The main body, which had high internal friction and shear near its base, then overran the surge deposits, partially eroding them. When the upward current of gas (fluidization) waned, better-sorted ash suspended in the fluidization zone was deposited on block-and-ash deposits. In the distal places of block-and-ash deposits, unit 2 probably was deposited in non-turbulent fashion without any erosion of the underlying layer (unit 1).     

Pyroclastic deposits of the November 13, 1985 eruption of Nevado del Ruiz volcano, Colombia

The November 13, 1985 eruption of Nevado del Ruiz produced a series of pyroclastic flows and surges that eroded channels on the surface of the summit glacier and generated lahars which descended down most of the rivers that drain the volcano. The stratigraphy of the proximal pyroclastic deposits indicates that there were at least four episodes to the eruption. Episode I, deposited an unusual surge consisting of small pieces of ice mixed with ash and exhibiting planar stratification. Ballistically emplaced fragments are also intercalated with this unit. During Episode II, at least two pyroclastic flows were erupted. Their deposits contain the most evolved pumice of the entire eruption; SiO₂ content of matrix glass ranges between 74.5 and 74.9%. Episode III is marked by the emplacement of a welded tuff with an average SiO₂ content of about 66% in the matrix glass. The final Episode IV was characterized by the development of a high-altitude eruption column and the emplacement of several nonwelded pyroclastic flows. Banded pumice are common in the pyroclastic flow as well as in the pumice fall deposits. Co-existing dark and light pumice bands differ in SiO₂ content by 3.5% and in general are similar to the composition of the welded pumice from Episode III.The compositional zonation of the pyroclastic deposits from Episode I to IV suggests that a nearsurface compositionally-stratified portion of the magma body was tapped during Episode II. During Episodes III and IV the main body of magma was involved although the coexistence of the compositionally distinct pumice clasts at similar stratigraphic levels argues for mixing of magma from different levels in the chamber during the eruptive process.     

Stochastic overland flows

A theoretical solution framework to the nonlinear stochastic partial differential equations (SPDE) of the kinematic wave and diffusion wave models of overland flows under stochastic inflows/outflows, stochastic surface roughness field and stochastic state of flows was obtained. This development was realized by means of an eigenfunction representation of the time-space overland flow depths, and by transforming the problem into the phase space. By using Van Kampen's lemma and the cumulant expansion theory of Kubo-Van Kampen-Fox, the deterministic partial differential equation (PDE) for the evolutionary probability density function (pdf) of overland flow depths was finally obtained. Once this deterministic PDE is solved for the time-varying pdf of overland flow depths, then the time-space varying pdf of overland flow depths can be obtained by a transformation given in the text. In this solution framework it is possible to incorporate the stochastic dynamic behavior of the parameters and of the forcing functions of the overland flow process. For example, not only the individual rainfall duration and fluctuating rain intensity characteristics but also the sequential behavior of rainfall patterns is incorporated into the evolutionary probability density function of overland flow depths.     

Segregation structures in vapor-differentiated basaltic flows

 Vesicle cylinders represent a spectacular kind of segregation structure involving residual liquids formed in situ during the cooling of lava flows. These vertical pipes, commonly found within basalt flows typically 2–10 m thick, are interpreted as the product of a vapor-driven differentiation process. The olivine phenocrysts and the earliest generation of groundmass olivines found in cylinder-bearing basalts appear to have been generally affected by magmatic oxidation, resulting in high-temperature iddingsite (HTI) alteration. This feature is also observed within cylinder-free basalt flows which exhibit other kinds of vesicular segregation structures, such as vesicle-rich pegmatoid segregation sheets and/or segregation vesicles. Detailed textural, petrological, and geochemical characteristics of two types of cylinders, three types of vesicle sheets, and five types of segregation vesicles are described, based on the study of 12 occurrences of HTI-bearing basalt flows from oceanic shield volcanoes or continental basalt plateaus. We propose a general classification of these segregation structures likely to derive from vapor differentiation. Flow thickness is probably the main factor influencing their morphology. Finally, we suggest that the concomitant occurrence of olivine oxidation and vapor-differentiation effects results from the late persistence of water oversaturation after eruption, perhaps due to a high rate of magma ascent. Received: 27 March 1999 / Accepted: 15 February 2000