Geochemistry of waters in the Valley of Ten Thousand Smokes region, Alaska

Meteoric waters from cold springs and streams outside of the 1912 eruptive deposits filling the Valley of Ten Thousand Smokes (VTTS) and in the upper parts of the two major rivers draining the 1912 deposits have similar chemical trends. Thermal springs issue in the mid-valley area along a 300-m lateral section of ash-flow tuff, and range in temperature from 21 to 29.8°C in early summer and from 15 to 17°C in mid-summer. Concentrations of major and minor chemical constituents in the thermal waters are nearly identical regardless of temperature. Waters in the downvalley parts of the rivers draining the 1912 deposits are mainly mixtures of cold meteoric waters and thermal waters of which the mid-valley thermal spring waters are representative. The weathering reactions of cold waters with the 1912 deposits appear to have stabilized and add only subordinate amounts of chemical constituents to the rivers relative to those contributed by the thermal waters. Isotopic data indicate that the mid-valley thermal spring waters are meteoric, but data is inconclusive regarding the heat source. The thermal waters could be either from a shallow part of a hydrothermal system beneath the 1912 vent region or from an incompletely cooled, welded tuff lens deep in the 1912 ash-flow sheet of the upper River Lethe area.Bicarbonate-sulfate waters resulting from interaction of near-surface waters and the cooling 1953–1968 southwest Trident plug issue from thermal springs south of Katmai Pass and near Mageik Creek, although the Mageik Creek spring waters are from a well-established, more deeply circulating hydrothermal system. Katmai caldera lake waters are a result of acid gases from vigorous drowned fumaroles dissolving in lake waters composed of snowmelt and precipitation.     

Complex proximal deposition during the Plinian eruptions of 1912 at Novarupta,Alaska

Proximal (<3 km) deposits from episodes II and III of the 60-h-long Novarupta 1912 eruption exhibit a very complex stratigraphy, the result of at least four transport regimes and diverse depositional mechanisms. They contrast with the relatively simple stratigraphy (and inferred emplacement mechanisms) for the previously documented, better known, medial–distal fall deposits and the Valley of Ten Thousand Smokes ignimbrite. The proximal products include alternations and mixtures of both locally and regionally dispersed fall ejecta, and numerous thin complex deposits of pyroclastic density currents (PDCs) with no regional analogs. The locally dispersed component of the fall deposits forms sector-confined wedges of material whose thicknesses halve radially from and concentrically about the vent over distances of 100–300 m (cf. several kilometers for the medial–distal fall deposits). This locally dispersed fall material (and many of the associated PDC deposits) is rich in andesitic and banded pumices and richer in shallow-derived wall-rock lithics in comparison with the coeval medial fall units of almost entirely dacitic composition. There are no marked contrasts in grain size in the near-vent deposits, however, between locally and widely dispersed beds, and all samples of the proximal fall deposits plot as a simple continuation of grain size trends for medial–distal samples. Associated PDC deposits form a spectrum of facies from fines-poor, avalanched beds through thin-bedded, landscape-mantling beds to channelized lobes of pumice-block-rich ignimbrite. The origins of the Novarupta near-vent deposits are considered within a spectrum of four transport regimes: (1) sustained buoyant plume, (2) fountaining with co-current flow, (3) fountaining with counter-current flow, and (4) direct lateral ejection. The Novarupta deposits suggest a model where buoyant, stable, regime-1 plumes characterized most of episodes II and III, but were accompanied by transient and variable partitioning of clasts into the other three regimes. Only one short period of vent blockage and cessation of the Plinian plume occurred, separating episodes II and III, which was followed by a single PDC interpreted as an overpressured "blast" involving direct lateral ejection. In contrast, regimes 2 and 3 were reflected by spasmodic sedimentation from the margins of the jet and perhaps lower plume, which were being strongly affected by short-lived instabilities. These instabilities in turn are inferred to be associated with heterogeneities in the mixture of gas and pyroclasts emerging from the vent. Of the parameters that control explosive eruptive behavior, only such sudden and asymmetrical changes in the particle concentration could operate on time scales sufficiently short to explain the rapid changes in the proximal 1912 products.Editorial responsibility: R. Cioni     

Geochemistry and isotopes of fluids from sulphur springs, Valles Caldera, New Mexico

Detailed geochemistry supported by geologic mapping has been used to investigate Sulphur Springs, an acid-sulfate hot spring system that issues from the western flank of the resurgent dome inside Valles Caldera. The most intense activity occurs at the intersection of faults offsetting caldera-fill deposits and post-caldera rhyolites. Three geothermal wells in the area have encountered pressures <1 MPa and temperatures of 200°C at depths of 600 to 1000 m. Hot spring and fumarole fluids may discharge at boiling temperatures with pH 1.0 and SO₄ 8000 mg/l. These conditions cause argillic alterations throughout a large area.Non-condensible gases consist of roughly 99% CO₂ with minor amounts of H₂S, H₂, and CH₄. Empirical gas geothermometry suggests a deep reservoir temperature of 215 to 280°C. Comparison of ¹³C and ¹⁸O between CaCO₃ from well cuttings and CO₂ from fumarole steam indicates a fractionation temperature between 200 and 300°C by decarbonation of hydrothermally altered Paleozoic limestone and vein calcite in the reservoir rocks. Tritium concentrations obtained from steam condensed in a mudpot and deep reservoir fluids (Baca #13, 278°C) are 2.1 and 1.0 T.U. respectively, suggesting the steam originates from a reservoir whose water is mostly >50 yrs old. Deuterium contents of fumarole steam, deep reservoir fluid, and local meteoric water are practically identical even though ¹⁸O contents range through 4‰, thus, precipitation on the resurgent dome of the caldera could recharge the hydrothermal system by slow percolation. From analysis of D and ¹⁸O values between fumarol steam and deep reservoir fluid, steam reaches the surface either (1) by vaporizing relatively shallow groundwater at 200°C or (2) by means of a two-stage boiling process through an intermediate level reservoir at roughly 200°C.Although many characteristics of known vapor-dominated geothermal systems are found at Sulphur Springs, fundamental differences exist in temperature and pressure of our postulated vapor-zone. We propose that the reservoir beneath Sulphur Springs is too small or too poorly confined to sustain a “true” vapor-dominated system and that the Sulphur Springs system may be a “dying” vapor-dominated system that has practically boiled itself dry.     

New evidence on the hydrothermal system in Long Valley caldera, California, from wells, fluid sampling, electrical geophysics, and age determinations of hot-spring deposits

Data collected since 1985 from test drilling, fluid sampling, and geologic and geophysical investigations provide a clearer definition of the hydrothermal system in Long Valley caldera than was previously available. This information confirms the existence of high-temperature (> 200°C) reservoirs within the volcanic fill in parts of the west moat. These reservoirs contain fluids which are chemically similar to thermal fluids encountered in the central and eastern parts of the caldera. The roots of the present-day hydrothermal system (the source reservoir, principal zones of upflow, and the magmatic heat source) most likely occur within metamorphic basement rocks beneath the western part of the caldera. Geothermometer-temperature estimates for the source reservoir range from 214 to 248°C. Zones of upflow of hot water could exist beneath the plateau of moat rhyolite located west of the resurgent dome or beneath Mammoth Mountain. Lateral flow of thermal water away from such upflow zones through reservoirs in the Bishop Tuff and early rhyolite accounts for temperature reversals encountered in most existing wells. Dating of hot-spring deposits from active and inactive thermal areas confirms previous interpretations of the evolution of hydrothermal activity that suggest two periods of extensive hot-spring discharge, one peaking about 300 ka and another extending from about 40 ka to the present. The onset of hydrothermal activity around 40 ka coincides with the initiation of rhyolitic volcanism along the Mono-Inyo Craters volcanic chain that extends beneath the caldera's west moat.     

New perspectives on the eruption of 1912 in the valley of ten thousand smokes,Katmai National Park,Alaska

New data extend our understanding of the 1912 eruption, its backfilled vent complex at Novarupta, and magma-storage systems beneath adjacent stratovolcanoes. Initial Plinian rhyolite fallout is confined to a narrow downwind sector, and its maximum thickness may occur as far as 13 km from source. In contrast, the partly contemporaneous rhyolite-rich ash flows underwent relatively low-energy emplacement, their generation evidently being decoupled from the high column. Flow veneers 1–13 m thick on near-vent ridge crests exhibit a general rhyolite-to-andesite sequence like that of the much thicker valley-confined ignimbrite into which they merge downslope. Lithics in both the initial Plinian and the ignimbrite are predominantly fragments of the Jurassic Naknek Formation, which extends from the surface to a depth of ca. 1500 m. Absence of lithics from the underlying sedimentary section limits to < 1.5 km the fragmentation level and the structural depth of the vent, which is thought to be funnel-shaped, flaring shallowly to a surface diameter of 2 km. Overlying the ignimbrite are layers of Plinian dacite fallout, > 100 m thick near source and 10 m thick 3 km away, which dip back into an inner vent <0.5 km wide, nested inside the earlier vent funnel of the ignimbrite. The dacite fallout is poor in Naknek lithics but contains abundant fragments of vitrophyre, most of which was vent-filling, densely welded tuff reejected during later phases of the 3-day eruption. Adjacent to the inner vent, a 225-m-high asymmetrical accumulation of coarse near-vent ejecta is stratigraphically continuous with the regional dacite fallout. Distensional faulting of its crest may reflect spreading related to compaction and welding. Nearby andesite-dacite stratovolcanoes, i.e., Martin, Mageik, Trident, and Katmai, display at least 12 vents that define a linear volcanic front trending N65°E. The 1912 vent and adjacent dacite domes are disposed parallel to the front and ca. 4 km behind it. Mount Griggs, 10 km behind the front, is more potassic than other centers, taps isotopically more depleted source materials, and reflects a wholly independent magmatic plumbing system. Geochemical differences among the stratovolcanoes, characteristically small eruptive volumes ( < 0.1 to 0.4 km³), and the dominance of andesite and low-SiO₂ dacite suggest complex crustal reservoirs, not large integrated magma chambers. Linear fractures just outside the 1912 vent strike nearly normal to the volcanic front and may reflect dike transport of magma previously stored beneath Trident 3–5 km away. Caldera collapse at Mount Katmai may have taken place in response to hydraulic transfer of Katmai magma toward Novarupta via reservoir components beneath Trident. The voluminous 1912 eruption (12–15 km³ DRE) was also unusual in producing high-silica rhyolite (6–9 km³ DRE), a composition rare in this arc and on volcanic fronts in general. Isotopic data indicate that rhyolite genesis involved little assimilation of sedimentary rocks, pre-Tertiary plutonic rocks, or hydrothermally altered rocks of any age. Trace-element data suggest nonetheless that the rhyolite contains a nontrivial crustal contribution, most likely partial melts of Late Cenozoic arc-intrusive rocks. Because the three compositions (77%, 66–64.5%, and 61.5–58.5% SiO₂) that intermingled in 1912 vented both concurrently and repeatedly (after eruptive pauses hours in duration), the compositional gaps between them must have been intrinsic to the reservoir, not merely effects of withdrawal dynamics.     

Abrupt transitions during sustained explosive eruptions: examples from the 1912 eruption of Novarupta, Alaska

Plinian/ignimbrite activity stopped briefly and abruptly 16 and 45 h after commencement of the 1912 Novarupta eruption defining three episodes of explosive volcanism before finally giving way after 60 h to effusion of lava domes. We focus here on the processes leading to the termination of the second and third of these three episodes. Early erupted pumice from both episodes show a very similar range in bulk vesicularity, but the modal values markedly decrease and the vesicularity range widens toward the end of Episode III. Clasts erupted at the end of each episode represent textural extremes; at the end of Episode II, clasts have very thin glass walls and a predominance of large bubbles, whereas at the end of Episode III, clasts have thick interstices and more small bubbles. Quantitatively, all clasts have very similar vesicle size distributions which show a division in the bubble population at 30 μm vesicle diameter and cumulative number densities ranging from 10⁷–10⁹ cm^–3. Patterns seen in histograms of volume fraction and the trends in the vesicle size data can be explained by coalescence signatures superimposed on an interval of prolonged nucleation and free growth of bubbles. Compared to experimental data for bubble growth in silicic melts, the high 1912 number densities suggest homogeneous nucleation was a significant if not dominant mechanism of bubble nucleation in the dacitic magma. The most distinct clast populations occurred toward the end of Plinian activity preceding effusive dome growth. Distributions skewed toward small sizes, thick walls, and teardrop vesicle shapes are indicative of bubble wall collapse marking maturation of the melt and onset of processes of outgassing. The data suggest that the superficially similar pauses in the 1912 eruption which marked the ends of episodes II and III had very different causes. Through Episode III, the trend in vesicle size data reflects a progressive shift in the degassing process from rapid magma ascent and coupled gas exsolution to slower ascent with partial open-system outgassing as a precursor to effusive dome growth. No such trend is visible in the Episode II clast assemblages; we suggest that external changes involving failure of the conduit/vent walls are more likely to have effected the break in explosive activity at 45 h.     

Reactivation of basement faults beneath volcanoes: a new model of flank collapse

In order to explain the presence of voluminous volcanic debris avalanche deposits around a stratovolcano, reactivation of vertical faults beneath a volcanic cone has been tested using analogue models. Reactivation of a single vertical fault beneath a cone generates a normal fault and an upturning of the layers creating a bulge on the flank. The upturning induces a flank collapse characterized by a typical horseshoe-shaped scar called an avalanche caldera. Reactivation of two vertical faults beneath a cone also generates a normal fault and a summit bulge. This bulge may result from the movement along a reverse fault. A large collapse is generated within the angle created by the two vertical faults. The angle of the collapse can be up to 140° whereas this angle is typically 120° for a dome intrusion. Collapse is instantaneous and is favoured by the presence of ductile layers (ash-and-pumice formations in the example considered) in a stratovolcano complex. The model may be applicable to volcanoes in a state of dormancy (or extinction) in regions with active regional tectonism. We suggest this mechanism of collapse in the case of the Cantal stratovolcano (Massif Central, France) to explain the presence of voluminous volcanic debris avalanche deposits around this volcano.     

Compositional heterogeneity of distal tephra deposits from the 1912 eruption of Novarupta, Alaska

Physical and chemical analyses of distal tephra from the 1912 eruption of Novarupta, Alaska, show considerable variations in glass and mineral compositions. A combination of a 150°C range in temperature deduced from iron-titanium oxide geothermometry, and curved patterns in bivariant element plots of glass compositions indicate that a chamber of compositionally zoned magma existed prior to the eruption. Magma-mixing cannot explain these features. The magma chamber may have resembled the model recently proposed by McBirney (1980). A highly silicic, quartz-phyric magma with mean phenocryst compositions of An₂₅ plagioclase, Fs₄₂ orthopyroxene, at a temperature of 880°C and a water pressure of 1.4 kbar, was located above a more mafic, hotter magma, bearing phenocrysts of An₄₅ plagioclase and Fs₃₅, orthopyroxene.Our results on distal tephras compare favorably with those from a recently completed study at source by Hildreth (1983), suggesting that useful petrologic information about distant volcanoes can be obtained from both types of deposits. Compositionally heterogeneous abyssal tephra layers are common in the Gulf of Alaska. Eruptions from chambers of zoned magma may account for many of these layers.     

New fumarolic activity on Mt. Baker: observations during April through July, 1975

On 10 March, 1975, large amounts of steam were first sighted coming from the Sherman Crater of Mt. Baker, long the site of mild fumarole activity. To document and evaluate the new activity, we conducted a series of aerial and ground-based observations and experiments including the collection of samples within the crater, installation of 35-mm automatic sequence cameras at Park and Anderson Buttes, and overflights with a mapping camera and infrared scanner. Significant changes characterizing the new activity are: development of a new, energetic fumarole in an area of earlier mild activity; growth of crevasses within and concentric to the crater walls; collapse of a central 70 m diameter plug of ice to form a warm lake; and ejection of tephra rich in analcite and sulfides from the new fumarole. Although the tephra contains minor amounts of unaltered glass, it does not differ in this or other respects from material collected from the west crater wall, where fumaroles have long been active. This suggests that the glass is older tephra derived from the fumarole vent walls.Weather and lighting conditions were such that the Park Butte sequence camera provided acceptable images of the crater area on 27% of the frames for the period 1825 Pacific Daylight Time June 5 to 1825 PDT August 2. Of these, 8% show a significant plume (>300 m above the vent). Large plumes occur during favorable weather conditions of low wind and high relative humidity, and all appear white.While there are abundant manifestations of increased heat emission from Sherman Crater, there is as yet no sign of a fundamental change in the style of activity, or that the rise of fresh magma is the cause.     

Fumarole distribution,morphology, and encrustation mineralogy associated with the 1986 eruptive deposits of mount St. Augustine,Alaska

Numerous rootless fumaroles were developed on pyroclastic flows and a lava flow generated during the March 1986 eruptive cycle of Mount St. Augustine. Gases issued from fumarole vents with four different shapes: fissure, phreatic explosion crater, single/multiple ovoid opening, and diffuse, multiple opening. Fumarole distribution and morphology were controlled by preeruption drainage and topography, as well as by the thickness, compaction, and settling of the flow deposits. Fumarole temperatures measured in June and July 1986 ranged from 75°–394°C. Varying amounts of colorful and often roughly zoned encrustations are associated with all fumarole vent shapes. Only six types of crystalline phases were detected by X-ray diffraction, with gypsum the most abundant mineral, followed by anhydrite, sulfur, tridymite, halite, and soda alum. Scanning electron microscopy and energy dispersive X-ray analysis revealed a number of amorphous phases, mainly halogen-rich, as well as other minor crystalline phases. The mineral assemblages in the encrustations suggest formation conditions for these deposits within a general range of 25°–250°C in an oxidizing environment. Many of the amorphous phases are metastable and upon cooling of the fumarole lose nonstructural water and crystallize to more stable forms. The high halogen contents of the fumarole condensates and the mineralogy, chemistry, and morphology of the encrustations support leaching of the andesitic ash and lava flow by condensed acid vapors as the primary source for the chemical components contained in the encrustations. Comparison of traceelement (Sr, Ba, V, Co, Ni, and Cr) contents in unaltered and altered ash suggests that trace-element distribution follows a pattern of isomorphic substitution in the encrustation phases.     

Fumaroles in ice caves on the summit of Mount Rainier—preliminary stable isotope, gas, and geochemical studies

The edifice of Mount Rainier, an active stratovolcano, has episodically collapsed leading to major debris flows. The largest debris flows are related to argillically altered rock which leave areas of the edifice prone to failure. The argillic alteration results from the neutralization of acidic magmatic gases that condense in a meteoric water hydrothermal system fed by the melting of a thick mantle of glacial ice. Two craters atop a 2000-year-old cone on the summit of the volcano contain the world's largest volcanic ice-cave system. In the spring of 1997 two active fumaroles (T=62°C) in the caves were sampled for stable isotopic, gas, and geochemical studies.Stable isotope data on fumarole condensates show significant excess deuterium with calculated δD and δ¹⁸O values (−234 and −33.2‰, respectively) for the vapor that are consistent with an origin as secondary steam from a shallow water table which has been heated by underlying magmatic–hydrothermal steam. Between 1982 and 1997, δD of the fumarole vapor may have decreased by 30‰.The compositions of fumarole gases vary in time and space but typically consist of air components slightly modified by their solubilities in water and additions of CO₂ and CH₄. The elevated CO₂ contents (δ¹³C_CO2=−11.8±0.7‰), with spikes of over 10,000 ppm, require the episodic addition of magmatic components into the underlying hydrothermal system. Although only traces of H₂S were detected in the fumaroles, most notably in a sample which had an air δ¹³C_CO2 signature (−8.8‰), incrustations around a dormant vent containing small amounts of acid sulfate minerals (natroalunite, minamiite, and woodhouseite) indicate higher H₂S (or possibly SO₂) concentrations in past fumarolic gases.Condensate samples from fumaroles are very dilute, slightly acidic, and enriched in elements observed in the much higher temperature fumaroles at Mount St. Helens (K and Na up to the ppm level; metals such as Al, Pb, Zn Fe and Mn up to the ppb level and volatiles such as Cl, S, and F up to the ppb level).The data indicate that the hydrothermal system in the edifice at Mount Rainier consists of meteoric water reservoirs, which receive gas and steam from an underlying magmatic system. At present the magmatic system is largely flooded by the meteoric water system. However, magmatic components have episodically vented at the surface as witnessed by the mineralogy of incrustations around inactive vents and gas compositions in the active fumaroles. The composition of fumarole gases during magmatic degassing is distinct and, if sustained, could be lethal. The extent to which hydrothermal alteration is currently occurring at depth, and its possible influence on future edifice collapse, may be determined with the aid of on site analyses of fumarole gases and seismic monitoring in the ice caves.     

Correlation of seismic activity and fumarole temperature at the Mt. Merapi volcano (Indonesia) in 2000

In this paper we present densely sampled fumarole temperature data, recorded continuously at a high-temperature fumarole of Mt. Merapi volcano (Indonesia). These temperature time series are correlated with continuous records of rainfall and seismic waveform data collected at the Indonesian–German multi-parameter monitoring network. The correlation analysis of fumarole temperature and precipitation data shows a clear influence of tropical rain events on fumarole temperature. In addition, there is some evidence that rainfall may influence seismicity rates, indicating interaction of meteoric water with the volcanic system. Knowledge about such interactions is important, as lava dome instabilities caused by heavy-precipitation events may result in pyroclastic flows. Apart from the strong external influences on fumarole temperature and seismicity rate, which may conceal smaller signals caused by volcanic degassing processes, the analysis of fumarole temperature and seismic data indicates a statistically significant correlation between a certain type of seismic activity and an increase in fumarole temperature. This certain type of seismic activity consists of a seismic cluster of several high-frequency transients and an ultra-long-period signal (<0.002 Hz), which are best observed using a broadband seismometer deployed at a distance of 600 m from the active lava dome. The corresponding change in fumarole temperature starts a few minutes after the ultra-long-period signal and simultaneously with the high-frequency seismic cluster. The change in fumarole temperature, an increase of 5 °C on average, resembles a smoothed step. Fifty-four occurrences of simultaneous high-frequency seismic cluster, ultra-long period signal and increase of fumarole temperature have been identified in the data set from August 2000 to January 2001. The observed signals appear to correspond to degassing processes in the summit region of Mt. Merapi.     

Two-dimensional kinematic and rheological modeling of the 1912 pyroclastic flow,Katmai, Alaska

Pyroclastic flow emplacement is strongly influenced by eruption column height. A surface along which kinetic energy is zero theoretically connects the loci of eruption column collapse with all coeval ignimbrite termini. This surface is reconstructed as a two-dimensional energy line for the 1912 Katmai pyroclastic flow in the Valley of Ten Thousand Smokes from mapped flow termini and the runup of the ignimbrite onto obstructions and through passes. Extrapolation of the energy line to the vicinity of the source vent at Novarupta suggests the eruption column which generated the ignimbrite eruption was approximately 425 m high. The 1912 pyroclastic flow travelled about 25 km downvalley. Empirical velocity data calculated from runup elevations and surveyed centrifugal superelevations indicate initial velocities near Novarupta were greater than 79–88 m s^–1. The flow progressively decelerated and was travelling only 2–8 m s^–1 when it crossed a moraine 16 km downvalley. The constant slope of the energy line away from Novarupta suggests the flow was systematically slowed by internal and basal friction. Using a simple physical model to calculate flow velocities and a constant kinetic friction coefficient (Heim coefficient) of 0.04 derived from the reconstructed energy line, the flow is estimated to have decelerated at an average rate of –0.16 m s^–2 and to have taken approximately 9.5 minutes to travel 25 km down the Valley of Ten Thousand Smokes. The shear strength of the flowing ignimbrite at the moraine was approximately 0.5 kPa, and its Bingham viscosity when it crossed the moraine was 3.5 × 10³ P. If the flow was Newtonian, its viscosity was 4.2 × 10³ P. Reynolds and Froude numbers at the moraine were only 41–62 and 0.84–1.04, respectively, indicating laminar, subcritical flow.     

Three-dimensional P-wave velocity structure derived from local earthquakes at the Katmai group of volcanoes,Alaska

The three-dimensional P-wave velocity structure beneath the Katmai group of volcanoes is determined by inversion of more than 10,000 rays from over 1000 earthquakes recorded on a local 18 station short-period network between September 1996 and May 2001. The inversion is well constrained from sea level to about 6 km below sea level and encompasses all of the Katmai volcanoes; Martin, Mageik, Trident, Griggs, Novarupta, Snowy, and Katmai caldera. The inversion reduced the average RMS travel-time error from 0.22 s for locations from the standard one-dimensional model to 0.13 s for the best three-dimensional model. The final model, from the 6th inversion step, reveals a prominent low velocity zone (3.6–5.0 km/s) centered at Katmai Pass and extending from Mageik to Trident volcanoes. The anomaly has values about 20–25% slower than velocities outboard of the region (5.0–6.5 km/s). Moderately low velocities (4.5–6.0 km/s) are observed along the volcanic axis between Martin and Katmai Caldera. Griggs volcano, located about 10 km behind (northwest of) the volcanic axis, has unremarkable velocities (5.0–5.7 km/s) compared to non-volcanic regions. The highest velocities are observed between Snowy and Griggs volcanoes (5.5–6.5 km/s). Relocated hypocenters for the best 3-D model are shifted significantly relative to the standard model with clusters of seismicity at Martin volcano shifting systematically deeper by about 1 km to depths of 0 to 4 km below sea level. Hypocenters for the Katmai Caldera are more tightly clustered, relocating beneath the 1912 scarp walls. The relocated hypocenters allow us to compare spatial frequency-size distributions (b-values) using one-dimensional and three-dimensional models. We find that the distribution of b is significantly changed for Martin volcano, which was characterized by variable values (0.8 < b < 2.0) with standard locations and more uniform values (0.8 < b < 1.2) after relocation. Other seismic clusters at Mageik (1.2 < b < 2.2), Trident (0.5 < b < 1.5) and Katmai Caldera (0.8 < b < 1.8) had stable b-values indicating the robustness of the observations. The strong high b-value region at Mageik volcano is mainly associated with an earthquake swarm in October, 1996 that possibly indicates a shallow intrusion or influx of gas. The new velocity and spatial b-value results, in conjunction with prior gravity (Bouguer anomalies up to − 40 mgal) and interferometry (several cm uplift) data, provide strong evidence in favor of partially molten rock at shallow depths beneath the Mageik–Katmai–Novarupta region. Moderately low velocities beneath Martin and Katmai suggest that old, mostly solidified intrusions exist beneath these volcanoes. Higher relative velocities beneath the Griggs and Snowy vents suggest that no magma is resident in the shallow crust beneath these volcanoes.     

Chemical evolution and volcanic activity of the active crater lake of Poás volcano, Costa Rica, 1993–1997

Concentrations of chloride and sulfate and pH in the hot crater lake (Laguna Caliente) at Poás volcano and in acid rain varied over the period 1993–1997. These parameters are related to changes in lake volume and temperature, and changes in summit seismicity and fumarole activity beneath the active crater. During this period, lake level increased from near zero to its highest level since 1953, lake temperature declined from a maximum value of 70°C to a minimum value of 25°C, and pH of the lake water increased from near zero to 1.8. In May 1993 when the lake was nearly dry, chloride and sulfate concentrations in the lake water reached 85,400 and 91,000 mg l⁻¹, respectively. Minimum concentrations of chloride and sulfate after the lake refilled to its maximum volume were 2630 and 4060 mg l⁻¹, respectively. Between January 1993 and May 1995, most fumarolic activity was focused through the bottom of the lake. After May 1995, fumarolic discharge through the bottom of the lake declined and reappeared outside the lake within the main crater area. The appearance of new fumaroles on the composite pyroclastic cone coincided with a dramatic decrease in type B seismicity after January 1996. Between May 1995 and December 1997, enhanced periods of type A seismicity and episodes of harmonic tremor were associated with an increase in the number of fumaroles and the intensity of degassing on the composite pyroclastic cone adjacent to the crater lake. Increases in summit seismic activity (type A, B and harmonic tremor) and in the height of eruption plumes through the lake bottom are associated with a period of enhanced volcanic activity during April–September 1994. At this time, visual observations and remote fumarole temperature measurements suggest an increase in the flux of heat and gases discharged through the bottom of the crater lake, possibly related to renewed magma ascent beneath the active crater. A similar period of enhanced seismic activity that occurred between August 1995 and January 1996, apparently caused fracturing of sealed fumarole conduits beneath the composite pyroclastic cone allowing the focus of fumarolic degassing to migrate from beneath the lake back to the 1953–1955 cone. Changes in the chemistry of summit acid rain are correlated changes in volcanic activity regardless of whether fumaroles are discharging into the lake or are discharging directly into the atmosphere.     

Alvin-SeaBeam studies of the Sumisu Rift, Izu-Bonin arc

Bimodal volcanism, normal faulting, rapid sedimentation, and hydrothermal circulation characterize the rifting of the Izu-Bonin arc at 31°N. Analysis of the zigzag pattern, in plan view, of the normal faults that bound Sumisu Rift indicates that the extension direction (080° ± 10°) is orthogonal to the regional trend of the volcanic front. Normal faults divide the rift into an inner rift on the arc side, which is the locus for maximum subsidence and sedimentation, and an outer rift further west. Transfer zones that link opposing master faults and/or rift flank uplifts further subdivide the rift into three segments along strike. Volcanism is concentrated along the ENE-trending transfer zone which separates the northern and central rift segments. The differential motion across the zone is accommodated by interdigitating north-trending normal faults rather than by ENE-trending oblique-slip faults. Volcanism in the outer rift has built 50–700 m high edifices without summit craters whereas in the inner rift it has formed two multi-vent en echelon ridges (the largest is 600 m high and 16 km long). The volcanism is dominantly basaltic, with compositions reflecting mantle sources little influenced by arc components. An elongate rhyolite dome and low-temperature hydrothermal deposits occur at the en echelon step in the larger ridge, which is located at the intersection of the transfer zone with the inner rift. The chimneys, veins, and crusts are composed of silica, barite and iron oxide, and are of similar composition to the ferruginous chert that mantles the Kuroko deposits. A 1.2-km transect of seven alvin heat flow measurements at 30°48.5′N showed that the inner-rift-bounding faults may serve as water recharge zones, but that they are not necessarily areas of focussed hydrothermal outflow, which instead occurs through the thick basin sediments. The rift basin and arc margin sediments are probably dominated by permeable rhyolitic pumice and ash erupted from submarine arc calderas such as Sumisu and South Sumisu volcanoes.     

Pyroclastic phases of a rhyolitic dome-building eruption: Puketarata tuff ring,Taupo Volcanic Zone,New Zealand

The 14 ka Puketarata eruption of Maroa caldera in Taupo Volcanic Zone was a dome-related event in which the bulk of the 0.25 km³ of eruption products were emplaced as phreatomagmatic fall and surge deposits. A rhyolitic dike encountered shallow groundwater during emplacement along a NE-trending normal fault, leading to shallow-seated explosions characterised by low to moderate water/magma ratios. The eruption products consist of two lava domes, a proximal tuff ring, three phreatic collapse craters, and a widespread fall deposit. The pyroclastic deposits contain dominantly dense juvenile clasts and few foreign lithics, and relate to very shallow-level disruption of the growing dome and its feeder dike with relatively little involvement of country rock. The distal fall deposit, representing 88% of the eruption products is, despite its uniform appearance and apparently subplinian dispersal, a composite feature equivalent to numerous discrete proximal phreatomagmatic lapilli fall layers, each deposited from a short-lived eruption column. The Puketarata products are subdivided into four units related to successive phases of:(A) shallow lava intrusion and initial dome growth; (B) rapid growth and destruction of dome lobes; (C) slower, sustained dome growth and restriction of explosive disruption to the dome margins; and (D) post-dome withdrawal of magma and crater-collapse. Phase D was phreatic, phases A and C had moderate water: magma ratios, and phase B a low water: magma ratio. Dome extrusion was most rapid during phase B, but so was destruction, and hence dome growth was largely accomplished during phase C. The Puketarata eruption illustrates how vent geometry and the presence of groundwater may control the style of silicic volcanism. Early activity was dominated by these external influences and sustained dome growth only followed after effective exclusion of external water from newly emplaced magma.     

Chemistry and mineralogy of high-temperature gas discharges from Colima volcano, Mexico. Implications for magmatic gas–atmosphere interaction

Gases, condensates and silica tube precipitates were collected from 400°C (Z2) and 800°C (Z3) fumaroles at Colima volcano, Mexico, in 1996–1998. Volcanic gases at Colima were very oxidized and contain up to 98% air due to mixing with air inside the dome interior, close to the hot magmatic body. An alkaline trap method was used to collect gas samples, therefore only acidic species were analysed. Colima volcanic gases are water-rich (95–98 mol%) and have typical S/C/Cl/F ratios for a subduction type volcano. δD-values for the high-temperature Z3 fumarolic vapour vary from −26 to −57‰. A negative δD–Cl correlation for the Z3 high-temperature fumarole may result from magma degassing: enrichment in D and decrease in the Cl concentration in condensates are likely a consequence of input of “fresh” batches of magma and an increasing of volcanic activity, respectively.The trace element composition of Colima condensates generally does not differ from that of other volcanoes (e.g. Merapi, Kudryavy) except for some enrichment in V, Cu and Zn. Variations in chemical composition of precipitates along the silica tube from the high-temperature fumarole (Colima 1, fumarole Z3), in contrast to other volcanoes, are characterized by high concentrations of Ca and V, low concentration of Mo and a lack of Cd. Mineralogy of precipitates differs significantly from that described for silica tube experiments at other volcanoes with reduced volcanic gas. Thermochemical modelling was used to explain why very oxidized gas at Colima does not precipitate halite, sylvite, and Mo- and Cd-minerals, but does precipitate V-minerals and native gold, which have not been observed before in mineral precipitates from reduced volcanic gases.     

GPR investigation of tephra fallout, Cerro Negro volcano, Nicaragua: a method for constraining parameters used in tephra sedimentation models

Most tephra fallout models rely on the advection–diffusion equation to forecast sedimentation and hence volcanic hazards. Here, we test the application of the advection–diffusion equation to tephra sedimentation using data collected on the proximal (350 to ~1,200?m from the vent) to medial (greater than ~1,200?m from the vent) tephra blanket of a basaltic cinder cone, Cerro Negro volcano, located in Nicaragua. Our understanding of tephra depositional processes at this volcano is significantly improved by combination of sample pit data in the medial zone and high-resolution ground-penetrating radar (GPR) data collected in the near vent and proximal zones. If the advection–diffusion equation applies, then the thickness of individual tephra deposits should have Gaussian crosswind profiles and exponential decay with distance away from the vent. At Cerro Negro, steady trade winds coupled with brief eruptions of relatively low energy (VEI 2–3) create relatively simple deposits. GPR data were collected along three crosswind profiles at distances of 700–1,600?m from the vent; sample pits were used to estimate thickness of the 1992 tephra deposit up to 13?km from the vent. Horizons identified in proximal GPR profiles exhibit Gaussian distributions with a high degree of statistical confidence, with diffusion coefficients of ~500?m²?s^?1 estimated for the deposits, confirming that the advection–diffusion equation is capable of modeling sedimentation in the proximal zone. The thinning trend downwind of the vent decreases exponentially from the cone base (350?m) to ~1,200?m from the vent. Beyond this distance, deposit overthickening occurs, identified in both GPR and sample pit datasets. The combined data reveal three depositional regimes: (1) a near-vent region on the cone itself, where fallout remobilizes in granular flows upon deposition; (2) a proximal zone in which particles fall from a height of less than ~2?km; and (3) a medial zone, in which particles fall from ~4 to 7?km and the deposit is thicker than expected based on thinning trends observed in the proximal zone of the deposit. This overthickening of the tephra blanket, defining the transition from proximal to medial depositional facies, is indicative of transition from sedimentation dominated by fallout from plume margins to that dominated by fallout from the buoyant eruption cloud—a feature of deposits previously identified in larger-volume eruptions. We interpret this change to represent a change in diffusion law, occurring at total particle fall times (the fall time threshold of numerical models) of ~400?s. Thus, the detailed GPR profiles and pit data collected at Cerro Negro help to validate current numerical models of tephra sedimentation.     

Basal layered deposits of the Peach Springs Tuff,northwestern Arizona,USA

Basal layered deposits of the large-volume Peach Springs Tuff occur beneath the main pyroclastic flow deposit over a minimum lateral distance of 70 km in northwestern Arizona (USA). The basal deposits are interpreted to record initial blasting and pyroclastic surge events at the beginning of the eruption; the pyroclastic surges traveled a minimum of 100 km from the (as yet unknown) source. Changes in bedding structures with increasing flow distance are related to the decreasing sediment load of the surges. Some bed forms in the most proximal part of the study area (Kingman, Arizona) can be interpreted as being shock induced, reflecting a blast origin for the surges. Component analyses support a hydrovolcanic origin for some of the blasting and subsequent pyroclastic surges. The eruption apparently began with magmatic blasts, which were replaced by hydrovolcanic blasts. Hydrovolcanic activity may be partially related to failure of the conduit walls that temporarily plugged the vent. A single large-volume pyroclastic flow immediately followed the blast phase, and no evidence has been observed for a Plinian eruption column. The stratigraphic sequence indicates that powerful hydrovolcanic blasting rapidly widened the vent, thus bypassing a Plinian fallout phase and causing rapid evolution to a collapsing eruption column. Similar processes may occur in other large-volume ignimbrite eruptions, which commonly lack significant Plinian fallout deposits.