A statistical analysis of flank eruptions on Etna volcano

A singularly complete record exists for the eruptive activity of Etna volcano. The time series of occurrence of flank eruptions in the period 1600–1980, in which the record is presumably complete, is found to follow a stationary Poisson process. A revision of the available data shows that eruption durations are rather well correlated with the estimates of the volume of lava flows. This implies that the magnitude of an eruption can be defined directly by its duration. Extreme value statistics are then applied to the time series, using duration as a dependent variable. The probability of occurrence of a very long (300 days) eruption is greater than 50% only in time intervals of the order of 50 years. The correlation found between duration and total output also allows estimation of the probability of occurrence of a major event which exceeds a given duration and total flow of lava. The composite probabilities do not differ considerably from the pure ones. Paralleling a well established application to seismic events, extreme value theory can be profitably used in volcanic risk estimates, provided that appropriate account is also taken of all other variables.     

Temporal variation in the lavas of Mayon volcano, Philippines

Chemical and petrographic analyses of 51 sequential lava flows from the central vent of Mayon volcano show cyclical variation. In the two most recent cycles, from 1800 to 1876 and from 1881 to the present, one to three basaltic flows are followed by six to ten andesitic flows. Modal and whole-rock chemical parameters show the most regular cyclical variation; calculated groundmass chemical parameters vary less regularly. There is also a long-term trend, over approximately 1700 years of exposed section, toward more basic compositions.The cyclical variation in modes and the chemical composition of the lavas apparently results from periodic influxes of basaltic magma from depth into a shallow magma system. Fractional crystallization of olivine, augite, hypersthene, calcic plagioclase, magnetite and pargasitic hornblende produces successively more andesitic lavas until the next influx of basaltic magma. Differentiation in a deep zone of magma generation is not excluded by the data, but is more likely responsible for the overall change toward more basic compositions than for the cyclical variation.Three points in a cycle — the beginning of basaltic lavas, the beginning of andesitic lavas and a leveling-off of SiO₂, K₂ O and K₂O/Na₂O values — correspond roughly to the beginning of frequent effusive eruptions (with or without an early Plinian eruption), frequent weak to moderately explosive (Strombolian) eruptions, and less frequent explosive (Vulcanian) eruptions, respectively. Recognition of the current stage in a cycle can give a qualitative indication of the nature of forthcoming eruptions. Changes in several specific parameters may precede basaltic lavas and allow early detection of basaltic influxes. These include minima in the glass inclusion/plagioclase phenocryst and phenocryst/groundmass ratios, vesicularity and groundmass TiO₂, a decrease in hypersthene phenocrysts, and constant values for the whole-rock K₂O/Na₂O ratio. The Mayon area is densely populated, making prediction of eruption type important for safety and land-use planning.     

Audiomagnetotelluric soundings on White Island volcano

Audiomagnetotelluric soundings in the frequency range 100-1 Hz have been made at three locations within the crater region of White Island volcano, New Zealand. Exceptionally low apparent resistivities are prevalent. Despite the possibility of terrain effects and static shift, and the three-dimensional nature of the island, the measured values appear to be indicative of the true electrical structure. One-dimensional modelling of the E polarisation responses has been used to obtain approximations to the electrical structure beneath each site. The results suggest a variable depth to the acid hydrothermal system which is believed to underly the crater region, with low resistivities coming closest to the surface near the active fumaroles of Donald Mound.     

Volume,energy and cyclicity of eruptions of Arenal volcano,Costa Rica

The 1968–73 (and continuing) eruption of Arenal Volcano, Costa Rica, a small 1633 m strato-volcano with long periods of repose, defines an eruptive cycle which is typical of Arenal’s pre-historic eruptions. An intense, short explosive phase (July 29–31, 1968) grades into an effusive phase, and is followed by a block lava flow. The eruptive rocks become increasingly less differentiated with time in a given cycle, ranging from andesite to basaltic andesite. Nuées ardentes are a characteristic of the initial explosions, and are caused by fall-back ejecta on slopes around the main crater — an explosion crater in the 1968 eruption — which coalesce into hot avalanches and descend major drainage channels. Total volume of pyroclastic flows was small, about 1.8 ± 0.5 × 10ⁿ m³, in the July 29–31 explosions, and are block and ash flows, with much accidental material. Overpressures, ranging up to perhaps 5 kilobars just prior to major explosions, were estimated from velocities of large ejected blocks, which had velocities of up to 600 m/sec. Total kinetic energy and volume of ejecta of all explosions are an estimated 3 × 10²² ergs and 0.03 km³, respectively. The block lava flow, emitted from Sept., 1968 to 1973 (and continuing) has a volume greater than 0.06 km³, and covers 2.7 km² at thicknesses ranging from 15 to over 100 m. The total volumes of the explosive and effusive phases for the 1968–73 eruption are about 0.05 km³ and 0.06 km³, respectively. The last eruption of Arenal occurred about 1500 AD. based on radiocarbon dating and archaeological means, and was about twice as voluminous as the current one (0.17 km³ versus 0.09 km³). The total thermal energies for this pre-historic eruption and the current one are 8 × 10²³ and 18 × 10²³, respectively. The total volume of Arenal’s cone is about 6 km³ from 1633 m (summit) to 500 m, and, estimates of age based on the average rate of cone growth from these two eruptions, suggest an age between 20,000 to 200,000 years.     

Video and seismic observations of Strombolian eruptions at Erebus volcano,Antarctica

Between 1986 and 1990 the eruptive activity of Erebus volcano was monitored by a video camera with on-screen time code and recorded on video tape. Corresponding seismic and acoustic signals were recorded from a network of 6 geophones and 2 infrasonic microphones. Two hundred Strombolian explosions and three lava flows which were erupted from 7 vents were captured on video. In December 1986 the Strombolian eruptions ejected bombs and ash. In November 1987 large bubble-bursting Strombolian eruptions were observed. The bubbles burst when the bubble walls thinned to ∼ 20 cm. Explosions with bomb flight-times up to 14.5 s were accompanied by seismic signals with our local size estimate, “unified magnitudes” (m_u), up to 2.3. Explosions in pools of lava formed by flows in the Inner Crater were comparatively weak.     

Electrical resistivity tomography study of Taal volcano hydrothermal system, Philippines

Taal volcano (311?m in altitude) is located in The Philippines (14°N, 121°E) and since 1572 has erupted 33 times, causing more than 2,000 casualties during the most violent eruptions. In March 2010, the shallow structures in areas where present-day surface activity takes place were investigated by DC resistivity surveys. Electrical resistivity tomography (ERT) lines were performed above the two identified hydrothermal areas located on the northern flank of the volcano and in the Main Crater, respectively. Due to rough topography, deep valleys, and dense vegetation, most measurements were collected using a remote method based on a laboratory-made equipment. This allowed retrieval of information down to a depth of 250?m. ERTs results detail the outlines of the two geothermal fields defined by previous self-potential, CO₂ soil degassing, ground temperature, and magnetic mapping (Harada et al. Japan Acad Sci 81:261–266, 2005; Zlotnicki et al. Bull Volcanol 71:29–49, 2009a, Phys Chem Earth 34:294–408, 2009b). Hydrothermal fluids originate mainly from inside the northern part of the Main Crater at a depth greater than the bottom of the Crater Lake, and flow upward to the ground surface. Furthermore, water from the Main Crater Lake infiltrates inside the surrounding geological formations. The hydrothermal fluids, outlined by gas releases and high temperatures, cross the crater rim and interact with the northern geothermal field located outside the Main Crater.     

Field survey of the 1994 Mindoro Island,Philippines tsunami

This is a report of the field survey of the November 15, 1994 Mindoro Island, Philippines, tsunami generated by an earthquake (M=7.0) with a strike-slip motion. We will report runup heights from 54 locations on Luzon, Mindoro and other smaller islands in the Cape Verde passage between Mindoro and Luzon. Most of the damage was concentrated along the northern coast of Mindoro. Runup height distribution ranged 3–4 m at the most severely damaged areas and 2–4 in neighboring areas. The tsunami-affected area was limited to within 10 km of the epicenter. The largest recorded runup value of 7.3 m was measured on the southwestern coast of Baco Island while a runup of 6.1 m was detected on its northern coastline. The earthquake and tsunami killed 62 people, injured 248 and destroyed 800 houses. As observed in other recent tsunami disasters, most of the casualties were children. Nearly all eyewitnesses interviewed described the first wave as a leading-depression wave. Eyewitnesses reported that the main direction of tsunami propagation was SW in Subaang Bay, SE in Wawa and Calapan, NE on Baco Island and N on Verde Island, suggesting that the tsunami source area was in the southern Pass of Verde Island and that the wave propagated rapidly in all directions. The fault plane extended offshore to the N of Mindoro Island, with its rupture originating S of Verde Island and propagating almost directly south to the inland of Mindoro, thereby accounting for the relatively limited damage area observed on the N of Mindoro.     

Geological and seismological evidence of increased explosivity during the 1986 eruptions of Pavlof volcano,Alaska

We present results of study of the best-documented eruptions of Pavlof volcano in historic time. The 1986 eruptions were mostly Strombolian in character; a strong initial phase may have been Vulcanian. The 1986 activity erupted at least 8×10⁶ m³ of feldspar-phyric basaltic andesite lava (SiO₂=53–54%), and a comparable volume of wind-borne tephra. During the course of the eruption, 5300 explosion earthquakes occurred, the largest of which was equivalent to an M _L =2.5 earthquake. Volcanic tremor was recorded for 2600 hours, and the strongest tremor was recorded out to a distance of 160 km and had an amplitude of at least 54 cm² reduced displacement. The 1986 eruptions modified the structure of the vent area for the first time in over two decades. A possible pyroclastic flow was observed on 19 June 1986, the first time such a phenomenon has been observed at the volcano. Overall, the 1986 eruptions were the strongest and longest duration eruptions in historic time, and changed a temporal pattern of activity that had persisted from 1973–1984.     

Vesiculation path of ascending magma in the 1983 and the 2000 eruptions of Miyakejima volcano, Japan

The vesiculation of magma during the 1983 eruption of Miyakejima Volcano, Japan, is discussed based on systematic investigations of water content, vesicularity, and bubble size distribution for the products. The eruption is characterized by simultaneous lava effusion and explosive sub-plinian (‘dry’) eruptions with phreatomagmatic (‘wet’) explosions. The magmas are homogeneous in composition (basaltic andesite) and in initial water content (H₂O = 3.9±0.9 wt%), and residual groundmass water contents for all eruption styles are low (H₂O <0.4 wt%) suggestive of extensive dehydration of magma. For the scoria erupted during simultaneous ‘dry’ and ‘wet’ explosive eruptions, inverse correlation was observed between vesicularity and residual water content. This relation can be explained by equilibrium exsolution and expansion of ca. 0.3 wt% H₂O at shallow level with different times of quenching, and suggests that each scoria with different vesicularity, which was quenched at a different time, provides a snapshot of the vesiculation process near the point of fragmentation. The bubble size distribution (BSD) varies systematically with vesicularity, and total bubble number density reaches a maximum value at vesicularity Φ ∼ 0.5. At Φ  ∼ 0.5, a large number of bubbles are connected with each other, and the average thickness of bubble walls reaches the minimum value below which they would rupture. These facts suggest that vesiculation advanced by nucleation and growth of bubbles when Φ < 0.5, and then by expansion of large bubbles with coalescence of small ones for Φ > 0.5, when bubble connection becomes effective. Low vesicularity and low residual water content of lava and spatter (Φ  < 0.1, H₂O  < 0.1 wt%), and systematic decrease in bubble number density from scoria through spatter to lava with decrease in vesicularity suggest that effusive eruption is a consequence of complete degassing by bubble coalescence and separation from magma at shallow levels when magma ascent rate is slow.

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

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

Major ash eruptions of Barren Island volcano (Andaman Sea) during the past 72 kyr: clues from a sediment core record

Barren Island (Andaman Sea) is the northernmost active volcano of the Indonesian Arc. To construct the eruptive history of this little studied volcano, we measured ¹⁴C dates of inorganic carbon in sediment beds, and Sr and Nd isotopic ratios of seven discrete ash layers, in a marine sediment core collected from 32 km southeast of the volcano. The study reveals that the volcano had seven major ash eruptions at ~70, 69, 61, 24, 19, 15, and 10 ka. The ash layers erupted from 70 ka through 19 ka have highly uniform Nd isotopic composition, and since the ~15 ka eruption to the present the isotopic composition has been highly variable. Between ~24 ka and ~10 ka, the volcano had large ash eruptions spaced at 4,500 year intervals. Isotopically correlating the precaldera lavas and ash exposed on the volcano to the uppermost ash layer in the core, we infer that the caldera of Barren Island volcano is younger than 10 ka.     

The submarine eruption of the Bombarda volcano,Milos Island,Cyclades, Greece

The Milos volcanic field includes a well-exposed volcaniclastic succession which records a long history of submarine explosive volcanism. The Bombarda volcano, a rhyolitic monogenetic center, erupted ∼1.7 Ma at a depth <200 m below sea level. The aphyric products are represented by a volcaniclastic apron (up to 50 m thick) and a lava dome. The apron is composed of pale gray juvenile fragments and accessory lithic clasts ranging from ash to blocks. The juvenile clasts are highly vesicular to non-vesicular; the vesicles are dominantly tube vesicles. The volcaniclastic apron is made up of three fades: massive to normally graded pumice-lithic breccia, stratified pumice-lithic breccia, and laminated ash with pumice blocks. We interpret the apron beds to be the result of water-supported, volcaniclastic mass-How emplacement, derived directly from the collapse of a small-volume, subaqueous eruption column and from syn-eruptive, down-slope resedimentation of volcaniclastic debris. During this eruptive phase, the activity could have involved a complex combination of phreatomagmatic explosions and minor submarine effusion. The lava dome, emplaced later in the source area, is made up of flow-banded lava and separated from the apron by an obsidian carapace a few meters thick. The near-vertical orientation of the carapace suggests that the dome was intruded within the apron. Remobilization of pyroclastic debris could have been triggered by seismic activity and the lava dome emplacement. Published online: 30 January 2003 Editorial responsibility: J. McPhie     

Seismic pattern recognition techniques to predict large eruptions at the Popocatépetl,Mexico, volcano

Using pattern recognition techniques, we formulate a simple prediction rule for a retrospective prediction of the three last largest eruptions of the Popocatépetl, Mexico, volcano that occurred on 23 April–30 June 1997 (Eruption 1; VEI ~ 2–3); 11 December 2000–23 January 2001 (Eruption 2; VEI ~ 3–4) and 7 June–4 September 2002 (Eruption 3; explosive dome extrusion and destruction phase). Times of Increased Probability (TIP) were estimated from the seismicity recorded by the local seismic network from 1 January 1995 to 31 December 2005. A TIP is issued when a cluster of seismic events occurs under our algorithm considerations in a temporal window several days (or weeks) prior to large volcanic activity providing sufficient time to organize an effective alert strategy. The best predictions of the three analyzed eruptions were obtained when averaging seismicity rate over a 5-day window with a threshold value of 12 events and declaring an alarm for 45 days. A TIP was issued about six weeks before Eruption 1. TIPs were detected about one and four weeks before Eruptions 2 and 3, respectively. According to our objectives, in all cases, the observed TIPs would have allowed the development of an effective civil protection strategy. Although, under our model considerations the three eruptive events were successfully predicted, one false alarm was also issued by our algorithm. An analysis of the epicentral and depth distribution of the local seismicity used by our prediction rule reveals that successful TIPs were issued from microearthquakes that took place below and towards SE of the crater. On the contrary, the seismicity that issued the observed false alarm was concentrated below the summit of the volcano. We conclude that recording of precursory seismicity below and SE of the crater together with detection of TIPs as described here, could become an important tool to predict future large eruptions at Popocatépetl. Although our model worked well for events that occurred in the past, it is necessary to verify the real capability of the model for future eruptive events.     

西马塔海底火山喷发的超远程水声监测

直接观察海底火山是非常困难的, 但监测其喷发产生的声信号可以为长期及近实时的海底火山监测提供一种有效手段.本文利用位于SOFAR声道轴的两套水听器三联体, 对西马塔海底火山开展了为期14个月的超远程水声监测研究.提出了基于加权慢度和聚类分析的火山喷发声信号检测方法, 采用抛物方程模型和绝热简正波-抛物方程模型对火山喷发活动进行了二维和三维的超远程声传播仿真, 火山喷发声接收距离超过9000 km.喷发活动频次和强度估算结果表明, 2009年1月—2010年2月该火山至少发生了37次喷发事件, 喷发活动的强度差异可达10 dB.

     

Reconstruction of the most recent volcanic eruptions from the Valles caldera,New Mexico

Products of the latest eruptions from the Valles caldera, New Mexico, consist of the El Cajete Pyroclastic Beds and Battleship Rock Ignimbrite, a sequence of pyroclastic fall and density current deposits erupted at ~ 55 ka, capped by the later Banco Bonito Flow erupted at ~ 40 ka, and collectively named the East Fork Member of the Valles Rhyolite. The stratigraphy of the East Fork Member has been the subject of conflicting interpretations in the past; a long-running investigation of short-lived exposures over a period of many years enables us to present a more complete event stratigraphy for these eruptions than has hitherto been possible. The volume of rhyolitic magma erupted during the 55 ka event may have been more than 10 km³, and for the 40 ka event can be estimated with rather more confidence at 4 km³. During the earlier event, plinian eruptions dispersed fallout pumice over much of the Valles caldera, the southern Jemez Mountains, and the Rio Grande rift. We infer a fallout thickness of several decimeters at the site of the city of Santa Fe, and significant ash fall in eastern New Mexico. In contrast, pyroclastic density currents were channeled within the caldera moat and southwestward into the head of Cañon de San Diego, the principal drainage from the caldera. Simultaneous (or rapidly alternating) pyroclastic fallout and density current activity characterized the ~ 55 ka event, with density currents becoming more frequent as the eruption progressed through two distinct stages separated by a brief hiatus. One early pyroclastic surge razed a forest in the southern caldera moat, in a similar manner to the initial blast of the May 18, 1980 eruption of Mt. St. Helens. Ignimbrite outflow from the caldera through the drainage notch may have been restricted in runout distance due to steep, rugged topography in this vicinity promoting mixing between flows and air, and the formation of phoenix clouds. Lavas erupted during both the ~ 55 and ~ 40 ka events were largely confined to the caldera moat. Any future rhyolitic eruptions of similar magnitude in the southern or western parts of the Valles caldera will likely affect similar areas.     

Eruptive patterns for the volcano piton de la fournaise, Reunion Island

A study of the historic record of activity of Piton de la Fournaise has revealed a cyclic pattern of eruption involving effusion of oceanite lava from major-flank centers every 20–40 years. Calculated volumes of the recent lava flows and pyroclastic ejecta have established an effusion rate of 3.9 m³ s⁻¹ since 1931 and 6.2 m³ s⁻¹ since 1951. Flank eruptions outside the present caldera define a distribution maximum which is expected to correlate with the depth range of a high-level magma reservoir.A model has been constructed which requires replenishment of a high-level magma chamber at a constant rate and regular eruption from summit and minor-flank centers, acting as “safety valves” to the magma chamber; when the magma chamber reaches its maximum expansion, a major-flank outburst of oceanitic lava occurs.The fact that calculated effusion rates are not consistent with radiometric dates implies an increase in effusion volume with time for the volcano.     

Magma discharge variations during the 2011 eruptions of Shinmoe-dake volcano, Japan, revealed by geodetic and satellite observations

We present precise geodetic and satellite observation-based estimations of the erupted volume and discharge rate of magma during the 2011 eruptions of Kirishima-Shinmoe-dake volcano, Japan. During these events, the type and intensity of eruption drastically changed within a week, with three major sub-Plinian eruptions on January 26 and 27, and a continuous lava extrusion from January 29 to 31. In response to each eruptive event, borehole-type tiltmeters detected deflation of a magma chamber caused by migration of magma to the surface. These measurements enabled us to estimate the geodetic volume change in the magma chamber caused by each eruptive event. Erupted volumes and discharge rates were constrained during lava extrusion using synthetic aperture radar satellite imaging of lava accumulation inside the summit crater. Combining the geodetic volume change and the volume of lava extrusion enabled the determination of the erupted volume and discharge rate during each sub-Plinian event. These precise estimates provide important information about magma storage conditions in magma chambers and eruption column dynamics, and indicate that the Shinmoe-dake eruptions occurred in a critical state between explosive and effusive eruption.     

Multiple edifice failures, debris avalanches and associated eruptions in the Holocene history of Shiveluch volcano, Kamchatka, Russia

 Investigation of well-exposed volcaniclastic deposits of Shiveluch volcano indicates that large-scale failures have occurred at least eight times in its history: approximately 10,000, 5700, 3700, 2600, 1600, 1000, 600 ¹⁴C BP and 1964 AD. The volcano was stable during the Late Pleistocene, when a large cone was formed (Old Shiveluch), and became unstable in the Holocene when repetitive collapses of a portion of the edifice (Young Shiveluch) generated debris avalanches. The transition in stability was connected with a change in composition of the erupting magma (increased SiO₂ from ca. 55–56% to 60–62%) that resulted in an abrupt increase of viscosity and the production of lava domes. Each failure was triggered by a disturbance of the volcanic edifice related to the ascent of a new batch of viscous magma. The failures occurred before magma intruded into the upper part of the edifice, suggesting that the trigger mechanism was indirectly associated with magma and involved shaking by a moderate to large volcanic earthquake and/or enhancement of edifice pore pressure due to pressurised juvenile gas. The failures typically included: (a) a retrogressive landslide involving backward rotation of slide blocks; (b) fragmentation of the leading blocks and their transformation into a debris avalanche, while the trailing slide blocks decelerate and soon come to rest; and (c) long-distance runout of the avalanche as a transient wave of debris with yield strength that glides on a thin weak layer of mixed facies developed at the avalanche base. All the failures of Young Shiveluch were immediately followed by explosive eruptions that developed along a similar pattern. The slope failure was the first event, followed by a plinian eruption accompanied by partial fountain collapse and the emplacement of pumice flows. In several cases the slope failure depressurised the hydrothermal system to cause phreatic explosions that preceded the magmatic eruption. The collapse-induced plinian eruptions were moderate-sized and ordinary events in the history of the volcano. No evidence for directed blasts was found associated with any of the slope failures. Received: 28 June 1998 / Accepted: 28 March 1999     

The Taurewa Eruptive Episode: evidence for climactic eruptions at Ruapehu volcano, New Zealand

 The ca. 10,500 years B.P. eruptions at Ruapehu volcano deposited 0.2–0.3 km³ of tephra on the flanks of Ruapehu and the surrounding ring plain and generated the only known pyroclastic flows from this volcano in the late Quaternary. Evidence of the eruptions is recorded in the stratigraphy of the volcanic ring plain and cone, where pyroclastic flow deposits and several lithologically similar tephra deposits are identified. These deposits are grouped into the newly defined Taurewa Formation and two members, Okupata Member (tephra-fall deposits) and Pourahu Member (pyroclastic flow deposits). These eruptions identify a brief (<ca. 2000-year) but explosive period of volcanism at Ruapehu, which we define as the Taurewa Eruptive Episode. This Episode represents the largest event within Ruapehu's ca. 22,500-year eruptive history and also marks its culmination in activity ca. 10,000 years B.P. Following this episode, Ruapehu volcano entered a ca. 8000-year period of relative quiescence. We propose that the episode began with the eruption of small-volume pyroclastic flows triggered by a magma-mingling event. Flows from this event travelled down valleys east and west of Ruapehu onto the upper volcanic ring plain, where their distal remnants are preserved. The genesis of these deposits is inferred from the remanent magnetisation of pumice and lithic clasts. We envisage contemporaneous eruption and emplacement of distal pumice-rich tephras and proximal welded tuff deposits. The potential for generation of pyroclastic flows during plinian eruptions at Ruapehu has not been previously considered in hazard assessments at this volcano. Recognition of these events in the volcanological record is thus an important new factor in future risk assessments and mitigation of volcanic risk at Tongariro Volcanic Centre. Received: 5 July 1998 / Accepted: 12 March 1999