A quantitative look at the demise of a basaltic vent: the death of Kupaianaha,Kilauea Volcano,Hawai'i

 The Kupaianaha vent, the source of the 48th episode of the 1983-to-present Pu'u 'O'o–Kupaianaha eruption, erupted nearly continuously from July 1986 until February 1992. This investigation documents the geophysical and geologic monitoring of the final 10 months of activity at the Kupaianaha vent. Detailed very low frequency (VLF) electromagnetic profiles across the single lava tube transporting lava from the vent were used to determine the cross-sectional area of the molten lava within the tube. Combined with measurements of lava velocity, these data provide an estimate of the lava output of Kupaianaha. In addition, lava temperatures (calculated from analysis of quenched glass) and bulk-rock chemistry were obtained for samples taken from the tube at the same site. The combined data set shows the lava flux from Kupaianaha vent declining linearly from 250 000 m³/day in April 1991 to 54 000 m³/day by November 1991. During that time surface breakouts of lava from weak points along the tube occurred progressively closer to the vent, consistent with declining efficiency in lava transport. There were no significant changes in lava temperature or in bulk MgO content during this period. Another eruptive episode (the 49th) began uprift of Kupaianaha on 8 November 1991 and erupted lava concurrently with Kupaianaha for 18 days. Lava flux from Kupaianaha decreased in response to this new episode, but the response was delayed by approximately 1 day. After 14 November 1991, lava velocities were no longer measurable in the tube because the lava stream beneath the skylight had crusted over; however, the VLF-derived electrical conductances documented the decreasing flux of molten lava through the tube. Kupaianaha remained active, but output continued to decrease until early February 1992 when the last active surface flows were seen. In November 1991 we used the linearly decreasing effusion rate to accurately predict the date for the death of the Kupaianaha vent. The linear nature of the decline in lava tube conductance and the delayed and slow response of the Waha'ula tube conductances to the 49th eruptive episode led us to speculate that (a) the Kupaianaha vent shut down because of a decrease in driving pressure and not because of a freeze-up of the vent, and (b) that Pu'u 'O'o, episode 49, and Kupaianaha were fed nearly vertically from a source deep within the rift zone. Received: 29 September 1995 / Accepted: 21 November 1995     

Thick lava flows of Karisimbi Volcano, Rwanda: insights from SIR-C interferometric topography

 We use a digital elevation model (DEM) derived from interferometrically processed SIR-C radar data to estimate the thickness of massive trachyte lava flows on the east flank of Karisimbi Volcano, Rwanda. The flows are as long as 12 km and average 40–60 m (up to >140 m) in thickness. By calculating and subtracting a reference surface from the DEM, we derived a map of flow thickness, which we used to calculate the volume (up to 1 km³ for an individual flow, and 1.8 km³ for all the identified flows) and yield strength of several flows (23–124 kPa). Using the DEM we estimated apparent viscosity based on the spacing of large folds (1.2×10¹² to 5.5×10¹² Pa s for surface viscosity, and 7.5×10¹⁰ to 5.2×10¹¹ Pa s for interior viscosity, for a strain interval of 24 h). We use shaded-relief images of the DEM to map basic flow structures such as channels, shear zones, and surface folds, as well as flow boundaries. The flow thickness map also proves invaluable in mapping flows where flow boundaries are indistinct and poorly expressed in the radar backscatter and shaded-relief images. Received: 6 September 1997 / Accepted: 15 May 1998     

Pahoehoe toe dimensions, morphology, and branching relationships at Mauna Ulu, Kilauea Volcano, Hawai'i

 Pahoehoe toe dimensions, morphology, and branching relationships were analyzed in flows emplaced during 1972 at Mauna Ulu, a satellitic shield on the east rift zone of Kilauea Volcano, Hawai'i. In order to characterize regions within flow fields dominated by networks of pahoehoe toes, measurements of toe length, width, thickness, and orientation were completed for 445 toes at 13 sites. Variations in site characteristics, including slope, substrate, and position in the flow field allow an evaluation of the effects of such parameters on toe dimensions. Toe surface morphology (ropy or smooth), local flow lobe position (interior or margin), and connective relationships between toes were documented in the form of detailed toe maps. These maps show the number of branches connecting a given toe to other toes in its local pahoehoe network and illustrate branching patterns. Statistical analyses of toe dimensions and comparisons of pahoehoe toe study sites and sub-populations combined with field observations, evaluation of toe maps, and qualitative examination of toe dimension size distributions show the following: (a) Although there are significant variations at a given site, toes typically have mean lengths (101 cm) greater than mean widths (74 cm) and mean widths greater than mean thicknesses (19 cm); sites that have mean widths greater than mean lengths are those with lower slopes. (b) Where significant site-to-site variations in mean values of a given toe dimension were apparent, these differences could not be directly related to site characteristics. (c) Ropy toes have significantly larger mean values of length, width, and number of branches than smooth toes, and toes with three or more branches have greater lengths, widths, and thicknesses than toes with two or fewer branches, suggesting concentration of flow in these toe types. (d) The skewness of all size distributions of toe length and width to larger values suggests that toes are transitional to larger sheets and channels, consistent with field observations; and (e) Two distinct types of branching patterns (called monolayer and centrally ridged) were observed in preserved pahoehoe flow lobes. The significant variability in measured toe dimensions at Mauna Ulu suggests that toe dimensions are influenced by numerous locally defined, random factors, and that an approach based on stochastic methods can be used to model pahoehoe flow emplacement. Received: 19 January 1998 / Accepted: 24 March 1999     

Statistical analyses and characteristics of volcanic tremor on Stromboli Volcano (Italy)

 A study of volcanic tremor on Stromboli is carried out on the basis of data recorded daily between 1993 and 1995 by a permanent seismic station (STR) located 1.8 km away from the active craters. We also consider the signal of a second station (TF1), which operated for a shorter time span. Changes in the spectral tremor characteristics can be related to modifications in volcanic activity, particularly to lava effusions and explosive sequences. Statistical analyses were carried out on a set of spectra calculated daily from seismic signals where explosion quakes were present or excluded. Principal component analysis and cluster analysis were applied to identify different classes of spectra. Three clusters of spectra are associated with two different states of volcanic activity. One cluster corresponds to a state of low to moderate activity, whereas the two other clusters are present during phases with a high magma column as inferred from the occurrence of lava fountains or effusions. We therefore conclude that variations in volcanic activity at Stromboli are usually linked to changes in the spectral characteristics of volcanic tremor. Site effects are evident when comparing the spectra calculated from signals synchronously recorded at STR and TF1. However, some major spectral peaks at both stations may reflect source properties. Statistical considerations and polarization analysis are in favor of a prevailing presence of P-waves in the tremor signal along with a position of the source northwest of the craters and at shallow depth. Received: 15 December 1996 / Accepted: 31 March 1998     

Continuous monitoring of volcanic eruption dynamics: a review of various techniques and new results from a frequency-modulated radar Doppler system

 In situ measurement of volcanic eruption velocities is one of the great challenges left in geophysical volcanology. In this paper we report on a new radar Doppler technique for monitoring volcanic eruption velocities. In comparison with techniques employed previously (e.g., photographic methods or acoustic Doppler measurements), this method allows continuous recordings of volcanic eruptions even during poor visibility. Also, radar Doppler instruments are usually light weight and energy efficient, which makes them superior to other Doppler techniques based on laser light or sound. The proposed new technique was successfully tested at Stromboli Volcano in late 1996 during a period of low activity. The recorded data allow a clear distinction between particles rising from the vent and particles falling back towards the vent. The mean eruption velocity was approximately 10 m/s. Most of the eruptions recorded by radar were correlated to seismic recordings. The correlation between the magnitude of the volcanic shocks and the eruption force index defined in the paper may provide new insights into magma transport in the conduit. Received: 15 May 1998 / Accepted: 15 December 1998     

New formulae for estimating lava flow volumes at Mt. Etna Volcano, Sicily

 A new data set of Etna lava flows erupted since 1868 has been compiled from eight topographic maps of the volcano published at intervals since then. Volumes of 59 flows or groups of flows were measured from topographic difference maps. Most of these volumes are likely to be considerably more accurate than those published previously. We cut the number of flow volumes down to 25 by selecting those examples for which the volume of an individual eruption could be derived with the highest accuracy. This refined data set was searched for high correlations between flow volume and more directly measurable parameters. Only two parameters showed a correlation coefficient of 70% or greater: planimetric flow area A (70%) and duration of the eruption D (79%). If only short duration (<18 days) flows were used, flow length cubed, L³, had a correlation coefficient of 98%. Using combinations of measured parameters, much more significant correlations with volume were found. Dh had a correlation coefficient of 90% (h is the hydrostatic head of magma above the vent), and  , 92% (where W is mean width and E is the degree of topographic enclosure), and a combination of the two , 97%. These latter formulae were used to derive volumes of all eruptions back to 1868 to compare with those from the complete data set. Values determined from the formulae were, on average, lower by 16% (Dh), 7% (, and 19% . Received: 30 November 1998 / Accepted: 20 June 1999     

Evolution of the seismic activity from the 1982 eruption of El Chichon Volcano, Chiapas, Mexico

 We analyzed more than 1700 earthquakes related to the 1982 eruption of El Chichon volcano in southern Mexico. The data were recorded at specific periods throughout the whole eruptive interval of March to April 1982, by three different networks. The seismic activity began several months before the first eruption on 28 March. During this period the seismicity consisted of hybrid and long-period shallow earthquakes most likely related to processes of faulting, fracturing, and fluid movement underneath the volcano. The foci of events occurring before the eruption circumscribe an aseismic zone from approximately 7 to 13 km below the volcano. After the eruption, the seismic activity consisted of tectonic-type earthquakes that peaked at 1200 events/h. This later activity occurred over a wide range of depths, mostly between 5 and 20 km, that includes the former aseismic zone and is roughly limited by the major tectonic faults in the area. Received: 19 May 1998 / Accepted: 13 June 1999     

The enormous Chillos Valley Lahar: an ash-flow-generated debris flow from Cotopaxi Volcano, Ecuador

 The Chillos Valley Lahar (CVL), the largest Holocene debris flow in area and volume as yet recognized in the northern Andes, formed on Cotopaxi volcano's north and northeast slopes and descended river systems that took it 326 km north–northwest to the Pacific Ocean and 130+ km east into the Amazon basin. In the Chillos Valley, 40 km downstream from the volcano, depths of 80–160 m and valley cross sections up to 337 000 m² are observed, implying peak flow discharges of 2.6–6.0 million m³/s. The overall volume of the CVL is estimated to be ≈3.8 km³. The CVL was generated approximately 4500 years BP by a rhyolitic ash flow that followed a small sector collapse on the north and northeast sides of Cotopaxi, which melted part of the volcano's icecap and transformed rapidly into the debris flow. The ash flow and resulting CVL have identical components, except for foreign fragments picked up along the flow path. Juvenile materials, including vitric ash, crystals, and pumice, comprise 80–90% of the lahar's deposit, whereas rhyolitic, dacitic, and andesitic lithics make up the remainder. The sand-size fraction and the 2- to 10-mm fraction together dominate the deposit, constituting ≈63 and ≈15 wt.% of the matrix, respectively, whereas the silt-size fraction averages less than ≈10 wt.% and the clay-size fraction less than 0.5 wt.%. Along the 326-km runout, these particle-size fractions vary little, as does the sorting coefficient (average=2.6). There is no tendency toward grading or improved sorting. Limited bulking is recognized. The CVL was an enormous non-cohesive debris flow, notable for its ash-flow origin and immense volume and peak discharge which gave it characteristics and a behavior akin to large cohesive mudflows. Significantly, then, ash-flow-generated debris flows can also achieve large volumes and cover great areas; thus, they can conceivably affect large populated regions far from their source. Especially dangerous, therefore, are snow-clad volcanoes with recent silicic ash-flow histories such as those found in the Andes and Alaska. Received: 28 April 1997 / Accepted: 19 August 1997     

Comment on " Volume of magma accumulation or withdrawal estimated from surface uplift or subsidence, with application to the 1960 collapse of Kīlauea volcano" by P. T. Delaney and D. F. McTigue

 In volcanoes that store a significant quantity of magma within a subsurface summit reservoir, such as Kīlauea, bulk compression of stored magma is an important mode of deformation. Accumulation of magma is also accompanied by crustal deformation, usually manifested at the surface as uplift. These two modes of deformation – bulk compression of resident magma and deformation of the volcanic edifice – act in concert to accommodate the volume of newly added magma. During deflation, the processes reverse and reservoir magma undergoes bulk decompression, the chamber contracts, and the ground surface subsides. Because magma compression plays a role in creating subsurface volume to accommodate magma, magma budget estimates that are derived from surface uplift observations without consideration of magma compression will underestimate actual magma volume changes. Received: 30 September 1998 / Accepted: 27 July 1999     

Postshield volcanism and catastrophic mass wasting of the Waianae Volcano, Oahu, Hawaii

 The 3.9- to 2.9-Ma Waianae Volcano is the older of two volcanoes making up the island of Oahu, Hawaii. Exposed on the volcanic edifice are tholeiitic shield lavas overlain by transitional and alkalic postshield lavas. The postshield "alkalic cap" consists of aphyric hawaiite of the Palehua Member of the Waianae Volcanics, overlain unconformably by a small volume of alkalic basalt of the Kolekole Volcanics. Kolekole Volcanics mantle erosional topography, including the uppermost slopes of the great Lualualei Valley on the lee side of the Waianae Range. Twenty new K–Ar dates, combined with magnetic polarity data and geologic relationships, constrain the ages of lavas of the Palehua member to 3.06–2.98 Ma and lavas of the Kolekole Volcanics to 2.97–2.90 Ma. The geochemical data and the nearly contemporaneous ages suggest that the Kolekole Volcanics do not represent a completely independent or separate volcanic event from earlier postshield activity; thus, the Kolekole Volcanics are reduced in rank, becoming the Kolekole Member of the Waianae Volcanics. Magmas of the Palehua and Kolekole Members have similar incompatible element ratios, and both suites show evidence for early crystallization of clinopyroxene consistent with evolution at high pressures below the edifice. However, lavas of the Kolekole Member are less fractionated and appear to have evolved at greater depths than the earlier Palehua hawaiites. Postshield primary magma compositions of the Palehua and Kolekole Members are consistent with formation by partial melting of mantle material of less than 5–10% relative to Waianae shield lavas. Within the section of Palehua Member lavas, an increase with respect to time of highly incompatible to moderately incompatible element ratios is consistent with a further decrease in partial melting by approximately 1–2%. This trend is reversed with the onset of eruption of Kolekole Member lavas, where an increase in extent of partial melting is indicated. The relatively short time interval between the eruption of Palehua and Kolekole Member lavas appears to date the initial formation of Lualualei Valley, which was accompanied by a marked change in magmatic conditions. We speculate that the mass-wasting event separating lavas of the Palehua and Kolekole Members may be related to the formation of a large submarine landslide west and southwest of Waianae Volcano. Enhanced decompression melting associated with removal of the equivalent volume of this landslide deposit from the edifice is more than sufficient to produce the modeled increase of 1–2% in extent of melting between the youngest Palehua magmas and the posterosional magmas of the Kolekole Member. The association between magmatic change and a giant landsliding event suggests that there may be a general relationship between large mass-wasting events and subsequent magmatism in Hawaiian volcano evolution. Received: 1 September 1996 / Accepted: 26 November 1996     

Holocene explosive activity of Hudson Volcano, southern Andes

 Fallout deposits in the vicinity of the southern Andean Hudson Volcano record at least 12 explosive Holocene eruptions, including that of August 1991 which produced ≥4 km³ of pyroclastic material. Medial isopachs of compacted fallout deposits for two of the prehistoric Hudson eruptions, dated at approximately 3600 and 6700 BP, enclose areas at least twice that of equivalent isopachs for both the 1991 Hudson and the 1932 Quizapu eruptions, the two largest in the Andes this century. However, lack of information for either the proximal or distal tephra deposits from these two prehistoric eruptions of Hudson precludes accurate volume estimates. Andesitic pyroclastic material produced by the 6700-BP event, including a  1 10-cm-thick layer of compacted tephra that constitutes a secondary thickness maximum over 900 km to the south in Tierra del Fuego, was dispersed in a more southerly direction than that of the 1991 Hudson eruption. The products of the 6700-BP event consist of a large proportion of fine pumiceous ash and accretionary lapilli, indicating a violent phreatomagmatic eruption. This eruption, which is considered to be the largest for Hudson and possibly for any volcano in the southern Andes during the Holocene, may have created Hudson's 10-km-diameter summit caldera, but the age of the caldera has not been dated independently. Received: 31 January 1997 / Accepted: 29 October 1997     

Reappraisal of the 1982 eruptions of El Chichón Volcano, Chiapas, Mexico: new data from proximal deposits

 Additional data from proximal areas enable a reconstruction of the stratigraphy and the eruptive chronology of phases III and IV of the 1982 eruption of El Chichón Volcano. Phase III began on 4 April at 0135 GMT with a powerful hydromagmatic explosion that generated radially fast-moving (∼100 ms^–1) pyroclastic clouds that produced a surge deposit (S1). Due to the sudden reduction in the confining pressure the process continued by tapping of magma from a deeper source, causing a new explosion. The ejected juvenile material mixed with large amounts of fragmented dome and wall rock, which were dispersed laterally in several pulses as lithic-rich block-and-ash flow (F1). Partial evacuation of juvenile material from the magmatic system prompted the entrance of external water to generate a series of hydromagmatic explosions that dispersed moisture-rich surge clouds and small-volume block-and-ash flows (IU) up to distances of 3 km from the crater. The eruption continued by further decompression of the magmatic system, with the ensuing emission of smaller amounts of gas-rich magma which, with the strong erosion of the volcanic conduit, formed a lithic-rich Plinian column that deposited fallout layer B. Associated with the widening of the vent, an increase in the effective density of the uprising column took place, causing its collapse. Block-and-ash flows arising from the column collapse traveled along valleys as a dense laminar flow (F2). In some places, flow regime changes due to topographic obstacles promoted transformation into a turbulent surge (S2) which attained minimum velocities of approximately 77 ms^–1 near the volcano. The process continued with the formation of a new column on 4 April at 1135 GMT (phase IV) that emplaced fall deposit C and was followed by hydromagmatic explosions which produced pyroclastic surges (S3). Received: 13 May 1996 / Accepted: 12 November 1996     

Stratigraphic framework of Holocene volcaniclastic deposits, Akutan Volcano, east-central Aleutian Islands, Alaska

 Akutan Volcano is one of the most active volcanoes in the Aleutian arc, but until recently little was known about its history and eruptive character. Following a brief but sustained period of intense seismic activity in March 1996, the Alaska Volcano Observatory began investigating the geology of the volcano and evaluating potential volcanic hazards that could affect residents of Akutan Island. During these studies new information was obtained about the Holocene eruptive history of the volcano on the basis of stratigraphic studies of volcaniclastic deposits and radiocarbon dating of associated buried soils and peat. A black, scoria-bearing, lapilli tephra, informally named the "Akutan tephra," is up to 2 m thick and is found over most of the island, primarily east of the volcano summit. Six radiocarbon ages on the humic fraction of soil A-horizons beneath the tephra indicate that the Akutan tephra was erupted approximately 1611 years B.P. At several locations the Akutan tephra is within a conformable stratigraphic sequence of pyroclastic-flow and lahar deposits that are all part of the same eruptive sequence. The thickness, widespread distribution, and conformable stratigraphic association with overlying pyroclastic-flow and lahar deposits indicate that the Akutan tephra likely records a major eruption of Akutan Volcano that may have formed the present summit caldera. Noncohesive lahar and pyroclastic-flow deposits that predate the Akutan tephra occur in the major valleys that head on the volcano and are evidence for six to eight earlier Holocene eruptions. These eruptions were strombolian to subplinian events that generated limited amounts of tephra and small pyroclastic flows that extended only a few kilometers from the vent. The pyroclastic flows melted snow and ice on the volcano flanks and formed lahars that traveled several kilometers down broad, formerly glaciated valleys, reaching the coast as thin, watery, hyperconcentrated flows or water floods. Slightly cohesive lahars in Hot Springs valley and Long valley could have formed from minor flank collapses of hydrothermally altered volcanic bedrock. These lahars may be unrelated to eruptive activity. Received: 31 August 1998 / Accepted: 30 January 1999     

Thermal monitoring of Lascar Volcano, Chile, using infrared data from the along-track scanning radiometer: a 1992–1995 time series

 Lascar Volcano (22°22'S, 67°44'W) is the most active volcano of the central Andes of northern Chile. Activity since 1984 has been characterised by periods of lava dome growth and decay within the active crater, punctuated by explosive eruptions. We present herein a technique for monitoring the high-temperature activity within the active crater using frequent measurements of emitted shortwave infrared (SWIR) radiation made by the spaceborne along-track scanning radiometer (ATSR). The ATSR is an instrument of low spatial resolution (pixels 1 km across) that shares certain characteristics with the MODIS instrument, planned for use as a volcano monitoring tool in the NASA EOS Volcanology Project. We present a comprehensive time series of over 60 cloud- and plume-free nighttime ATSR observations for 1992–1995, a period during which Lascar experienced its largest historical eruption. Variations in short wavelength infrared flux relate directly to changes in high-temperature surfaces within the active crater. From these data, interpretations can be made that supplement published field reports and that can document the presence and status of the lava dome during periods where direct, ground-based, observations are lacking. Our data agree with less frequent information collected from sensors with high spatial resolution, such as the Landsat thematic mapper (Oppenheimer et al. 1993) and are consistent with field observations and models that relate subsidence of the dome to subsequent explosive eruptions (Matthews et al., 1997). Most obviously, Lascar's major April 1993 eruption follows a period in which the magnitude of emitted shortwave infrared radiation fell by 90%. At this time subsidence of the 1991–1992 lava dome was reported by field observers and this subsidence is believed to have impeded the escape of hot volatiles and ultimately triggered the eruption (Smithsonian Institution 1993a). Extrapolating beyond the period for which field observations of the summit are available, our data show that the vulcanian eruption of 20 July 1995 occurred after a period of gradual increase in short wavelength infrared flux throughout 1994 and a more rapid flux decline during 1995. We attribute this additional, otherwise undocumented, cycle of increasing and decreasing SWIR radiance as most likely representing variations in degassing through fumaroles contained within the summit crater. Alternatively, it may reflect a cycle of dome growth and decay. The explosive eruption of 17 December 1993 appears to have followed a similar, but shorter, variation in SWIR flux, and we conclude that large explosive eruptions are more likely when the 1.6-μm signal has fallen from a high to a low level. The ATSR instrument offers low-cost data at high temporal resolution. Despite the low spatial detail of the measurements, ATSR-type instruments can provide data that relate directly to the status of Lascar's lava dome and other high-temperature surfaces. We suggest that such data can therefore assist with predictions of eruptive behaviour, deduced from application of physical models of lava dome development at this and similar volcanoes. Received: 1 October 1996 / Accepted: 13 January 1997     

The 1984 to 1996 cyclic activity of Lascar Volcano, northern Chile: cycles of dome growth, dome subsidence, degassing and explosive eruptions

 Lascar Volcano (5592 m; 23°22'S, 67°44'W) entered a new period of vigorous activity in 1984, culminating in a major explosive eruption in April 1993. Activity since 1984 has been characterised by cyclic behaviour with recognition of four cycles up to the end of 1993. In each cycle a lava dome is extruded in the active crater, accompanied by vigorous degassing through high-temperature, high-velocity fumaroles distributed on and around the dome. The fumaroles are the source of a sustained steam plume above the volcano. The dome then subsides back into the conduit. During the subsidence phase the velocity and gas output of the fumaroles decrease, and the cycle is completed by violent explosive activity. Subsidence of both the dome and the crater floor is accommodated by movement on concentric, cylindrical or inward-dipping conical fractures. The observations are consistent with a model in which gas loss from the dome is progressively inhibited during a cycle and gas pressure increases within and below the lava dome, triggering a large explosive eruption. Factors that can lead to a decrease in gas loss include a decrease in magma permeability by foam collapse, reduction in permeability due to precipitation of hydrothermal minerals in the pores and fractures within the dome and in country rock surrounding the conduit, and closure of open fractures during subsidence of the dome and crater floor. Dome subsidence may be a consequence of reduction in magma porosity (foam collapse) as degassing occurs and pressurisation develops as the permeability of the dome and conduit system decreases. Superimposed upon this activity are small explosive events of shallow origin. These we interpret as subsidence events on the concentric fractures leading to short-term pressure increases just below the crater floor. Received: 12 December 1996 / Accepted: 6 May 1997     

Geochemical and isotopic evidence for seawater contamination of the hydrothermal system of Taal Volcano, Luzon, the Philippines

 The hydrologic structure of Taal Volcano has favored development of an extensive hydrothermal system whose prominent feature is the acidic Main Crater Lake (pH<3) lying in the center of an active vent complex, which is surrounded by a slightly alkaline caldera lake (Lake Taal). This peculiar situation makes Taal prone to frequent, and sometimes catastrophic, hydrovolcanic eruptions. Fumaroles, hot springs, and lake waters were sampled in 1991, 1992, and 1995 in order to develop a geochemical model for the hydrothermal system. The low-temperature fumarole compositions indicate strong interaction of magmatic vapors with the hydrothermal system under relatively oxidizing conditions. The thermal waters consist of highly, moderately, and weakly mineralized solutions, but none of them corresponds to either water–rock equilibrium or rock dissolution. The concentrated discharges have high Na contents (>3500 mg/kg) and low SO₄/Cl ratios (<0.3). The Br/Cl ratio of most samples suggests incorporation of seawater into the hydrothermal system. Water and dissolved sulfate isotopic compositions reveal that the Main Crater Lake and spring discharges are derived from a deep parent fluid (T≈300  °C), which is a mixture of seawater, volcanic water, and Lake Taal water. The volcanic end member is probably produced in the magmatic-hydrothermal environment during absorption of high-temperature gases into groundwater. Boiling and mixing of the parent water give rise to the range of chemical and isotopic characteristics observed in the thermal discharges. Incursion of seawater from the coastal region to the central part of the volcano is supported by the low water levels of the lakes and by the fact that Lake Taal was directly connected to the China sea until the sixteenth century. The depth to the seawater-meteoric water interface is calculated to be 80 and 160 m for the Main Crater Lake and Lake Taal, respectively. Additional data are required to infer the hydrologic structure of Taal. Geochemical surveillance of the Main Crater Lake using the SO₄/Cl, Na/K, or Mg/Cl ratio cannot be applied straightforwardly due to the presence of seawater in the hydrothermal system. Received: 12 February 1997 / Accepted: 26 January 1998     

Tsunami generation by pyroclastic flow during the 3500-year B.P. caldera-forming eruption of Aniakchak Volcano, Alaska

 A discontinuous pumiceous sand, a few centimeters to tens of centimeters thick, is located up to 15 m above mean high tide within Holocene peat along the northern Bristol Bay coastline of Alaska. The bed consists of fine-to-coarse, poorly to moderately well-sorted, pumice-bearing sand near the top of a 2-m-thick peat sequence. The sand bed contains rip-up clasts of peat and tephra and is unique in the peat sequence. Major element compositions of juvenile glass from the deposit and radiocarbon dating of enclosing peat support correlation of the pumiceous sand with the caldera-forming eruption of Aniakchak Volcano. The distribution of the sand and its sedimentary characteristics are consistent with emplacement by tsunami. The pumiceous sand most likely represents redeposition by tsunami of climactic fallout tephra and beach sand during the approximately 3.5 ka Aniakchak caldera-forming eruption on the Alaska Peninsula. We propose that a tsunami was generated by the sudden entrance of a rapidly moving, voluminous pyroclastic flow from Aniakchak into Bristol Bay. A seismic trigger for the tsunami is unlikely, because tectonic structures suitable for tsunami generation are present only south of the Alaska Peninsula. The pumiceous sand in coastal peat of northern Bristol Bay is the first documented geologic evidence of a tsunami initiated by a volcanic eruption in Alaska. Received: 3 December 1997 / Accepted: 11 April 1998     

Kinematic significance of mingling-rolling structures in lava flows: a case study from Porri Volcano (Salina, Southern Tyrrhenian Sea)

 A basaltic andesite lava flow from Porri Volcano (Salina, Southern Tyrrhenian Sea) is composed of two different magmas. Magma A (51 vol.% of crystals) has a dacitic glass composition, and magma B (18 vol.% of crystals), a basaltic glass composition. Magma B is hosted in A and consists of sub-spherical enclaves and boudin-like, banding and rolling structures (RS). Four types of RS have been recognized: σ–type;δ–type; complex σ-δ–types and transitional structures between sub-spherical enclaves and rolling structures. An analysis of the RS has been performed in order to reconstruct the flow kinematics and the mechanism of flow emplacement. Rolling structures have been selected in three sites located at different distances from the vent. In all sites most RS show the same sense of shear. Kinematic analysis of RS allows the degree of flow non-coaxiality to be determined. The non-coaxiality is expressed by the kinematic vorticity number Wk, a measure of the ratio Sr between pure shear strain rate and simple shear strain rate. The values of Wk calculated from the measured shapes of microscopic RS increase with increasing distance from the vent, from approximately 0.5 to 0.9. Results of the structural analysis reveal that the RS formed during the early–intermediate stage of flow emplacement. They represent originally sub-spherical enclaves deformed at low shear strain. At higher strain, RS deformed to give boudin-like and stretched banding structures. Results of the kinematic analysis suggest that high viscosity lava flows are heterogeneous non-ideal shear flows in which the degree of non-coaxiality increases with the distance from the vent. In the vent area, deformation is intermediate between simple shear and pure shear. Farther from the vent, deformation approaches ideal simple shear. Lateral extension processes occur only in the near-vent zone, where they develop in response to the lateral push of magma extruded from the vent. Lateral shortening processes develop in the distal zone and record the gravity-driven movement of the lava. The lava flow advanced by two main mechanisms, lateral translation and rolling motion. Lateral translation equals rolling near the vent, while rolling motion prevailed in the distal zones. Received: 6 November 1997 / Accepted: 29 November 1997     

Petrology of lavas from the Puu Oo eruption of Kilauea Volcano: III. The Kupaianaha episode (1986–1992)

 The Puu Oo eruption has been remarkable in the historical record of Kilauea Volcano for its duration (over 13 years), volume (>1 km³) and compositional variation (5.7–10 wt.% MgO). During the summer of 1986, the main vent for lava production moved 3 km down the east rift zone and the eruption style changed from episodic geyser-like fountaining at Puu Oo to virtually continuous, relatively quiescent effusion at the Kupaianaha vent. This paper examines this next chapter in the Puu Oo eruption, episodes 48 and 49, and presents new ICP-MS trace element and Pb-, Sr-, and Nd-isotope data for the entire eruption (1983–1994). Nearly aphyric to weakly olivine-phyric lavas were erupted during episodes 48 and 49. The variation in MgO content of Kupaianaha lavas erupted before 1990 correlates with changes in tilt at the summit of Kilauea, both of which probably were controlled by variations in Kilauea's magma supply rate. These lavas contain euhedral olivines which generally are in equilibrium with whole-rock compositions, although some of the more mafic lavas which erupted during 1990, a period of frequent pauses in the eruption, accumulated 2–4 vol.% olivine. The highest forsterite content of olivines (∼85%) in Kupaianaha lavas indicates that the parental magmas for these lavas had MgO contents of ∼10 wt.%, which equals the highest observed value for lavas during this eruption. The composition of the Puu Oo lavas has progressively changed during the eruption. Since early 1985 (episode 30), when mixing between an evolved rift zone magma and a more mafic summit reservoir-derived magma ended, the normalized (to 10 wt.% MgO) abundances of highly incompatible elements and CaO have systematically decreased with time, whereas ratios of these trace elements and Pb, Sr, and Nd isotopes, and the abundances of Y and Yb, have remained relatively unchanged. These results indicate that the Hawaiian plume source for Puu Oo magmas must be relatively homogeneous on a scale of 10–20 km³ (assuming 5–10% partial melting), and that localized melting within the plume has apparently progressively depleted its incompatible elements and clinopyroxene component as the eruption continued. The rate of variation of highly incompatible elements in Puu Oo lavas is much greater than that observed for Kilauea historical summit lavas (e.g., Ba/Y 0.09 a^–1 vs ∼0.03 a^–1). This rapid change indicates that Puu Oo magmas did not mix thoroughly with magma in the summit reservoir. Thus, except for variable amounts of olivine fractionation, the geochemical variation in these lavas is predominantly controlled by mantle processes. Received: 8 March 1996 / Accepted: 30 April 1996