Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Volcanology and Geothermal Research”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2Linked to original sources

Assessment of leachable elements in volcanic ashfall: A review and evaluation of a standardized protocol for ash hazard characterization

Volcanic ash presents a widespread and common hazard during and after eruptions. Complex interactions between solid ash surfaces and volcanic gases lead to the formation of soluble salts that may be mobilized in aqueous environments. A variety of stakeholders may be concerned about the effects of ash on human and animal health, drinking water supplies, crops, soils and surface runoff. As part of the immediate emergency response, rapid dissemination of information regarding potentially hazardous concentrations of soluble species is critical. However, substantial variability in the methods used to characterize leachable elements makes it challenging to compare datasets and eruption impacts. To address these challenges, the International Volcanic Health Hazard Network ( www.ivhhn.org ) organized a two-day workshop to define appropriate methods for hazard assessment. The outcome of this workshop was a ‘consensus protocol’ for analysis of volcanic ash samples for rapid assessment of hazards from leachable elements, which was subsequently ratified by leading volcanological organizations. The purpose of this protocol is to recommend clear, standard and reliable methods applicable to a range of purposes during eruption response, such as assessing impacts on drinking-water supplies and ingestion hazards to livestock, and also applicable to research purposes. Where possible, it is intended that the methods make use of commonly available equipment and require little training. To evaluate method transferability, an interlaboratory comparison exercise was organized among six laboratories worldwide. Each laboratory received a split of pristine ash, and independently analyzed it according to the protocol for a wide range of elements. Collated results indicate good repeatability and reproducibility for most elements, thus indicating that the development of this protocol is a useful step towards providing standardized and reliable methods for ash hazard characterization. In this article, we review recent ash leachate studies, report the outcomes of the comparison exercise and present a revised and updated protocol based on the experiences and recommendations of the exercise participants. The adoption of standardized methods will improve and facilitate the comparability of results among studies and enable the ongoing development of a global database of leachate information relevant for informing volcanic health hazards assessment.

Journal of Volcanology and Geothermal Research↗

The thermal regime in the resurgent dome of Long Valley Caldera, California: Inferences from precision temperature logs in deep wells

Long Valley Caldera in eastern California formed 0.76Ma ago in a cataclysmic eruption that resulted in the deposition of 600km3 of Bishop Tuff. The total current heat flow from the caldera floor is estimated to be ~290MW, and a geothermal power plant in Casa Diablo on the flanks of the resurgent dome (RD) generates ~40MWe. The RD in the center of the caldera was uplifted by ~80cm between 1980 and 1999 and was explained by most models as a response to magma intrusion into the shallow crust. This unrest has led to extensive research on geothermal resources and volcanic hazards in the caldera. Here we present results from precise, high-resolution, temperature-depth profiles in five deep boreholes (327-1,158m) on the RD to assess its thermal state, and more specifically 1) to provide bounds on the advective heat transport as a guide for future geothermal exploration, 2) to provide constraints on the occurrence of magma at shallow crustal depths, and 3) to provide a baseline for future transient thermal phenomena in response to large earthquakes, volcanic activity, or geothermal production. The temperature profiles display substantial non-linearity within each profile and variability between the different profiles. All profiles display significant temperature reversals with depth and temperature gradients <50??C/km at their bottom. The maximum temperature in the individual boreholes ranges between 124.7??C and 129.5??C and bottom hole temperatures range between 99.4??C and 129.5??C. The high-temperature units in the three Fumarole Valley boreholes are at the approximate same elevation as the high-temperature unit in borehole M-1 in Casa Diablo indicating lateral or sub-lateral hydrothermal flow through the resurgent dome. Small differences in temperature between measurements in consecutive years in three of the wells suggest slow cooling of the shallow hydrothermal flow system. By matching theoretical curves to segments of the measured temperature profiles, we calculate horizontal groundwater velocities in the hydrothermal flow unit under the RD that range from 1.9 to 2.8m/yr, which corresponds to a maximum power flowing through the RD of 3-4MW. The relatively low temperatures and large isothermal segments at the bottom of the temperature profiles are inconsistent with the presence of magma at shallow crustal levels. ?? 2010.

Journal of Volcanology and Geothermal Research↗

Submarine geology of Hana Ridge and Haleakala Volcano's northeast flank, Maui

We present a morphostructural analysis of the submarine portions of Haleakala Volcano and environs, based upon a 4-year program of geophysical surveys and submersible explorations of the underwater flanks of Hawaiian volcanoes that was conducted by numerous academic and governmental research organizations in Japan and the U.S. and funded primarily by the Japan Agency for Marine–Earth Science and Technology. A resulting reconnaissance geologic map features the 135-km-long Hana Ridge, the 3000 km 2 Hana slump on the volcano's northeast flank, and island-surrounding terraces that are the submerged parts of volcanic shields. Hana Ridge below 2000 m water depth exhibits the lobate morphology typical of the subaqueously erupted parts of Hawaiian rift zones, with some important distinctions: namely, subparallel crestlines, which we propose result from the down-rift migration of offsets in the dike intrusion zone, and an amphitheater at its distal toe, where a submarine landslide has embayed the ridge tip. Deformation of Haleakala's northeast flank is limited to that part identified as the Hana slump, which lies downslope from the volcano's submerged shield, indicating that flank mobility is also limited in plan, inconsistent with hypothesized volcanic spreading driven by rift-zone dilation. The leading edge of the slump has transverse basins and ridges that resemble the thrust ramps of accretionary prisms, and we present a model to describe the slump's development that emphasizes the role of coastally generated fragmental basalt on gravitational instability of Haleakala's northeast flank and that may be broadly applicable to other ocean-island slumps.

Hawai'i↗

Water volumes, heat flow, and solute discharge from Old Faithful Geyser eruptions, Yellowstone National Park, USA

The iconic Old Faithful Geyser in Yellowstone National Park, USA, has attracted a significant amount of research because of the relative regularity and impressive size of its eruptions. Numerous studies have included observations, measurements, and analyses that informed models of geyser eruptions. However, fundamental quantities, including the associated mass and heat discharged, remain poorly constrained. In April 2025 we measured the volume of water from 45 Old Faithful Geyser eruptions using a portable flume in an outflow channel and specific conductance measurements in the Firehole River. We used high-speed video to perform velocimetry, measured changes in water chemistry to calculate the volume of water evaporated along the outflow channels, and used thermodynamic calculations to estimate the volume of water erupted as steam and to quantify the geyser's heat output. The calculated average volume of water discharged by Old Faithful Geyser in each eruption is 27.9 ± 9.4 m 3 , with no relation between eruption volume and the length of the preceding eruption interval. Video analysis of the eruption's liquid-dominated phase yields similar volumes of 21–30 m 3 . The calculated heat flow from the geyser is 2.2–2.4 MW and the average annual discharge of chloride, fluoride, and arsenic are 63 tons, 3.9 tons, and 241 kg, respectively. Average annual silica deposition rate on the geyser cone and along the outflow channels is 7 tons. This study provides a methodology for future studies at geysers worldwide and a baseline for monitoring future activity changes at Old Faithful.

Wyoming↗

Historical eruptions of Merapi Volcano, Central Java, Indonesia, 1768-1998

Information on Merapi eruptive activity is scattered and much is remotely located. A concise and well-documented summary of this activity has been long needed to assist researchers and hazard-mitigation efforts, and the aim of this paper is to synthesize information from the mid-1700s to the present. A descriptive chronology is given, with an abbreviated chronology in a table that summarizes events by year, assigns preliminary Volcanic Explosivity Index (VEI) ratings and Hartmann classifications, and provides key references. The history of volcano monitoring is also outlined. The study reveals that a major difference in eruption style exists between the twentieth and nineteenth centuries, although the periodicity between larger events seems about the same. During the twentieth century, activity has comprised mainly the effusive growth of viscous lava domes and lava tongues, with occasional gravitational collapses of parts of oversteepened domes to produce the nue??es ardentes - commonly defined as "Merapi-type". In the 1800s, however, explosive eruptions of relatively large size occurred (to VEI 4), and some associated "fountain-collapse" nue??es ardentes were larger and farther reaching than any produced in the twentieth century. These events may also be regarded as typical eruptions for Merapi. The nineteenth century activity is consistent with the long-term pattern of one relatively large event every one or two centuries, based on the long-term eruptive record deduced by others from volcanic stratigraphy. It is uncertain whether or not a "recurrence-time" model continues to apply to Merapi, but if so, Merapi could soon be due for another large event and its occurrence with only modest (or inadequately appreciated) precursors could lead to a disaster unprecedented in Merapi's history because the area around the volcano is now much more densely populated.

Java, Merapi Volcano↗

Surface morphology of caldera-forming eruption deposits revealed by lidar mapping of Crater Lake National Park, Oregon- Implications for emplacement and surface modification

Large explosive eruptions of silicic magma can produce widespread pumice fall, extensive ignimbrite sheets, and collapse calderas. The surfaces of voluminous ignimbrites are rarely preserved or documented because most terrestrial examples are heavily vegetated, or severely modified by post-depositional processes. Much research addresses the internal sedimentary characteristics, flow processes, and depositional mechanisms of ignimbrites, however, surface features of ignimbrites are less well documented and understood, except for comparatively small-volume deposits of historical eruptions. The ~7,700 calendar year B.P. climactic eruption of Mount Mazama, USA vented ~50 km3 of magma, deposited first as rhyodacite pumice fall and then as a zoned rhyodacite-to-andesite ignimbrite as Crater Lake caldera collapsed. Lidar collected during summer 2010 reveals the remarkably well-preserved surface of the Mazama ignimbrite and related deposits surrounding Crater Lake caldera in unprecedented detail despite forest cover. The ±1 m lateral and ±4 cm vertical resolution lidar allows surface morphologies to be classified. Surface morphologies are created by internal depositional processes and can point to the processes at work when pyroclastic flows come to rest. We describe nine surface features including furrow-ridge sets and wedge-shaped mounds in pumice fall eroded by high-energy pyroclastic surges, flow- parallel ridges that record the passage of multiple pyroclastic flows, perched benches of marginal deposits stranded by more-mobile pyroclastic-flow cores, hummocks of dense clasts interpreted as lag deposit, transverse ridges that mark the compression and imbrication of flows as they came to rest, scarps indicating ignimbrite remobilization, fields of pit craters caused by phreatic explosions, fractures and cracks caused by extensional processes resulting from ignimbrite volume loss, and stream channels eroded in the newly formed surface. The nine morphologies presented here illustrate a dynamic depositional environment that varied spatially and with time during the eruption, and show that multiple processes modified the ignimbrite after deposition, both during and after the eruption.

Oregon↗

Constraints on the mechanism of long-term, steady subsidence at Medicine Lake volcano, northern California, from GPS, leveling, and InSAR

Leveling surveys across Medicine Lake volcano (MLV) have documented subsidence that is centered on the summit caldera and decays symmetrically on the flanks of the edifice. Possible mechanisms for this deformation include fluid withdrawal from a subsurface reservoir, cooling/crystallization of subsurface magma, loading by the volcano and dense intrusions, and crustal thinning due to tectonic extension (Dzurisin et al., 1991 [Dzurisin, D., Donnelly-Nolan, J.M., Evans, J.R., Walter, S.R., 1991. Crustal subsidence, seismicity, and structure near Medicine Lake Volcano, California. Journal of Geophysical Research 96, 16, 319-16, 333.]; Dzurisin et al., 2002 [Dzurisin, D., Poland, M.P., Bürgmann, R., 2002. Steady subsidence of Medicine Lake Volcano, Northern California, revealed by repeated leveling surveys. Journal of Geophysical Research 107, 2372, doi:10.1029/2001JB000893.]). InSAR data that approximate vertical displacements are similar to the leveling results; however, vertical deformation data alone are not sufficient to distinguish between source mechanisms. Horizontal displacements from GPS were collected in the Mt. Shasta/MLV region in 1996, 1999, 2000, 2003, and 2004. These results suggest that the region is part of the western Oregon block that is rotating about an Euler pole in eastern Oregon. With this rotation removed, most sites in the network have negligible velocities except for those near MLV caldera. There, measured horizontal velocities are less than predicted from ∼10 km deep point and dislocation sources of volume loss based on the leveling data; therefore volumetric losses simulated by these sources are probably not causing the observed subsidence at MLV. This result demonstrates that elastic models of subsurface volume change can provide misleading results where additional geophysical and geological constraints are unavailable, or if only vertical deformation is known. The deformation source must be capable of causing broad vertical deformation with comparatively smaller horizontal displacements. Thermoelastic contraction of a column of hot rock beneath the volcano cannot reproduce the observed ratio of vertical to horizontal surface displacements. Models that determine deformation due to loading by the volcano and dense intrusions can be made to fit the pattern of vertical displacements by assuming a weak upper crust beneath MLV, though the subsidence rates due to surface loading must be lower than the observed displacements. Tectonic extension is almost certainly occurring based on fault orientations and focal mechanisms, but does not appear to be a major contributor to the observed deformation. We favor a model that includes a combination of sources, including extension and loading of a hot weak crust with thermal contraction of a cooling mass of rock beneath MLV, which are processes that are probably occurring at MLV. Future microgravity surveys and the planned deployment of an array of continuous GPS stations as part of a Plate Boundary Observatory volcano cluster will help to refine this model.

California↗

Petrologic and geochemical characterization of the first year of the 2024–ongoing episodic Halemaumau eruption of Kilauea (Episodes 1–39)

Kīlauea has been erupting episodically since December 23, 2024 (and is ongoing at the time of this writing in 2026) with two consistently active vents in Halema‘uma‘u producing lava fountains ranging in height from 40 to 540 m. In the first year of the eruption (through December 24, 2025, episode 39), the U.S. Geological Survey's Hawaiian Volcano Observatory collected 178 samples for eruption monitoring purposes and used a subset of 80 for near-real-time geochemical analyses. Tephra, lava, and recycled lava lake ooze-out samples have been studied on inter-episode timescales using diverse microanalytical techniques, where research-scale investigations have proven to be advantageous pseudo-near-real-time techniques. Here we present a geochemical and petrologic overview of episodes 1–39 to better understand eruption dynamics and identify changes at Kīlauea following the 2018 summit caldera collapse. There are no significant intra-episode geochemical changes that might be related to variable fountain heights or episode durations, suggesting that eruption dynamics are instead primarily controlled by viscosity and conduit geometry changes throughout a single episode and between episodes. There is variable but clear influence of the still-molten lava lake interior on compositions of this eruption, including in tephra samples from lava fountains, indicating that the vent conduits are connected in some way at depth with the still molten passive lava lake core. Tephra glasses (7.0–7.8 wt.% MgO) show cyclic changes in MgO that were decoupled from lava fountain heights and behavior in episodes 1–28 but positively correlated after episode 29. Additionally, whole-rock and glass compositions are enriched in minor elements such as K 2 O and TiO 2 compared to 2008–2023 summit eruptions, tracking the arrival of a new magma composition that may be related to the prolonged episodic eruptive behavior. Olivine-hosted melt inclusions from five episodes have ~0.35–0.43 wt.% H2O, < 300 ppm CO 2 and < 1500 ppm S, with no significant variations between episodes. Melt and fluid inclusion barometry show that magmas driving these episodes are stored and sourced from the shallow 1–2 km Halemaʻumaʻu reservoir, as were olivine from the 2023 summit eruption. The ongoing episodic eruption also has higher bulk-rock MgO compared to the five previous summit eruptions (2020−2023) and the forsterite (Fo) content of olivine has varied systematically from 2015 to 2025 with the equilibrium olivine composition of the Halema‘uma‘u reservoir increasing from Fo 80.7 in 2015 to Fo 82.7 by 2025. We argue that the Halemaʻumaʻu reservoir is heating up because the rate of magma supply is outpacing effusion and cooling following caldera collapse in 2018.

Hawaii↗

Globally detected volcanic lightning and umbrella dynamics during the 2014 eruption of Kelud, Indonesia

Volcanic lightning shows considerable promise as a monitoring and research tool to characterize explosive eruptions. Its key strengths are rapid and remote detection, because the radio signals produced by lightning can propagate thousands of km at the speed of light. Despite these tantalizing properties, the scientific work on volcanic lightning has only recently started gaining momentum. Much more is needed to understand what lightning reveals about the evolution of an eruption in near-real time. Here we examine the timing and energy release of lightning generated by the eruption of Kelud volcano in Indonesia on 13 February 2014, as detected by the World Wide Lightning Location Network (WWLLN). The eruption column reached at least 26 km above sea level, representing the highest plume since the advent of global lightning networks in the last decade. Therefore, it provides valuable constraints on the electrification of end-member, sustained Plinian columns. We investigate the lightning in context with satellite images, photographs, and other published studies. Results show that the earliest satellite-detected activity was a thermal anomaly at ~15:46 UTC, corresponding to a directed blast at the onset of eruption (and only a few lightning strokes). Following a brief pause, the eruption produced a sustained column and umbrella cloud that spread outward into the tropical stratosphere. Rates of umbrella expansion provide an average mass eruption rate (MER) in the range of 8 × 10 7 –1 × 10 8 kg s −1 . A more nuanced picture emerges from the time-varying MERs (determined between each satellite pass), which show rapid intensification during the first hour of eruption, followed by constant MER for about an hour, and waning toward the end (after ~17:50 UTC). At this stage, decreasing flux into the umbrella cloud coincides with column instability and formation of pyroclastic density currents, as recorded by photos from the ground ~17:45 UTC. We infer that some of the erupted mass partitioned into ground-hugging currents, leading to a lower apparent MER. Interestingly, there is not a 1:1 correlation between lightning intensity and MER over the course of eruption. Stroke rates increase sharply within the first 30–40 min (during rapid intensification of the plume), and then drop below 2 strokes per min once the MER remains constant. This suggests that electrification was controlled by the rate of increase in MER—in other words, the acceleration of particles out of the vent. We also show that lightning stroke-rates and energies are greatest within 50 km of the vent, even when the ash cloud extends >200 km downwind, indicating that lightning was focused in the regions of highest particle concentration and turbulence. Overall, we conclude that abrupt changes in lightning rates are clearly linked to changes in eruption behavior, and that rapid detection could aid monitoring efforts to characterize eruption rates or styles.

Journal of Volcanology and Geothermal Research↗

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

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

California↗

Transition of basaltic lava from pahoehoe to aa, Kilauea Volcano, Hawaii: Field observations and key factors

Nearly all Hawaiian basaltic lava erupts as pahoehoe, and some changes to aa during flowage and cooling; factors governing the transition involve certain critical relations between viscosity and rate of shear strain. If the lava slows, cools, and stops in direct response to concomitant increase in viscosity before these critical relations are reached, it remains pahoehoe. But, if flow mechanics (flow rate, flow dimensions, slope, momentum, etc.) impel the lava to continue to move and deform even after it has become highly viscous, the critical relations may be reached and the lava changes to aa. Typical modes of transition from pahoehoe to aa include: (1) spontaneous formation of relatively stiff clots in parts of the flowing lava where shear rate is highest; these clots grow into discrete, rough, sticky masses to which the remaining fluid lava incrementally adheres; (2) fragmentation and immersion of solid or semi-solid surface crusts of pahoehoe by roiling movements of the flow, forming cores of discrete, tacky masses; (3) sudden renewed movement of lava stored and cooled within surface reservoirs to form clots. The masses, fragments, and clots in these transition modes are characterized by spinose, granulated surfaces; as flow movement continues, the masses and fragments aggregate, fracture, and grind together, completing the transition to aa. Observations show that the critical relation between viscosity and rate of shear strain is inverse: if viscosity is low, a high rate of shear is required to begin the transition to aa; conversely, if viscosity is high, a much lower rate of shear will induce the transition. These relations can be demonstrated qualitatively with simple graphs, which can be used to examine the flow history of any selected finite lava element by tracing the path represented by its changing viscosity and shear rate. A broad, diffuse “transition threshold zone” in these graphs portrays the inverse critical relation between viscosity and shear rate; the transition to aa is represented by the path of the lava element crossing this zone. Moving lava flows can be regarded as natural viscometers, by which shear stress and rate of shear strain at selected points can be determined and viscosity can be computed. By making such determinations under a wide range of conditions on pahoehoe, aa, and transitional flow types, the critical relations that control the pahoehoe-aa transition can be quantified.

Hawaii↗

Huge landslide blocks in the growth of piton de la fournaise, La réunion, and Kilauea volcano, Hawaii

Piton de la Fournaise, on the island of La R&eacute;union, and Kilauea volcano, on the island of Hawaii, are active, basaltic shield volcanoes growing on the flanks of much larger shield volcanoes in intraplate tectonic environments. Past studies have shown that the average rate of magma production and the chemistry of lavas are quite similar for both volcanoes. We propose a structural similarity &mdash; specifically, that periodic displacement of parts of the shields as huge landslide blocks is a common mode of growth. In each instance, the unstable blocks are within a rift-zone-bounded, unbuttressed flank of the shield. At Kilauea, well-documented landslide blocks form relatively surficial parts of a much larger rift-zone-bounded block; scarps of the Hilina fault system mark the headwalls of the active blocks. At Fournaise, Hilina-like slump blocks are also present along the unbuttressed east coast of the volcano. In addition, however, the existence of a set of faults nested around the present caldera and northeast and southeast rift zones suggests that past chapters in the history of Fournaise included the slumping of entire rift-zone-bounded blocks themselves. These nested faults become younger to the east southeast and apparently record one of the effects of a migration of the focus of volcanism in that direction. Repeated dilation along the present set of northeast and southeast rift zones, most recently exemplified by an eruption in 1977, suggests that the past history of rift-zone-bounded slumping will eventually be repeated. The record provided by the succession of slump blocks on Fournaise is apparently at a relatively detailed part of a migration of magmatic focus that has advanced at least 30 km to the east-southeast from neighboring Piton des Neiges, an extinct Pliocene to Pleistocene volcano. ?? 1982.

Journal of Volcanology and Geothermal Research↗

Leachability of uranium and other elements from freshly erupted volcanic ash

A study of leaching of freshly erupted basaltic and dacitic air-fall ash and bomb fragment samples, unaffected by rain, shows that glass dissolution is the dominant process by which uranium is initially mobilized from air-fall volcanic ash. Si, Li, and V are also preferentially mobilized by glass dissolution. Gaseous transfer followed by fixation of soluble uranium species on volcanic-ash particles is not an important process affecting uranium mobility. Gaseous transfer, however, may be important in forming water-soluble phases, adsorbed to ash surfaces, enriched in the economically and environmentally important elements Zn, Cu, Cd, Pb, B, F, and Ba. Quick removal of these adsorbed elements by the first exposure of freshly erupted ash to rain and surface water may pose short-term hazards to certain forms of aquatic and terrestrial life. Such rapid release of material may also represent the first step in transportation of economically important elements to environments favorable for precipitation into deposits of commercial interest. Ash samples collected from the active Guatemalan volcanoes Fuego and Pacaya (high-Al basalts) and Santiaguito (hornblende-hypersthene dacite); bomb fragments from Augustine volcano (andesite-dacite), Alaska, and Heimaey (basalt), Vestmann Islands, Iceland; and fragments of "rhyolitic" pumice from various historic eruptions were subjected to three successive leaches with a constant water-to-ash weight ratio of 4:1. The volcanic material was successively leached by: (1) distilled-deionized water (pH = 5.0-5.5) at room temperature for 24 h, which removes water-soluble gases and salts adsorbed on ash surfaces during eruption; (2) dilute HCl solution (pH = 3.5-4.0) at room temperature for 24 h, which continues the attack initiated by the water and also attacks acid-soluble sulfides and oxides; (3) a solution 0.05 M in both Na,CO, and NaHCO, (pH = 9.9) at 80&deg;C for one week, which preferentially dissolves volcanic glass. The first two leaches mimic interaction of ash with rain produced in the vicinity of an active eruption. The third leach accelerates the effect of prolonged contact of volcanic ash with alkaline ground water present during ash diagenesis.

Journal of Volcanology and Geothermal Research↗

The compositionally zoned eruption of 1912 in the Valley of Ten Thousand Smokes, Katmai National Park, Alaska

On June 6–8, 1912, ∼ 15 km 3 of magma erupted from the Novarupta caldera at the head of the Valley of Ten Thousand Smokes (VTTS), producing ∼ 20 km 3 of air-fall tephra and 11–15 km 3 of ash-flow tuff within ∼ 60 hours. Three discrete periods of ash-fall at Kodiak correlate, respectively, with Plinian tephra layers designated A, CD, and FG by Curtis (1968) in the VTTS. The ash-flow sequence overlapped with but outlasted pumice fall A, terminating within 20 hours of the initial outbreak and prior to pumice fall C. Layers E and H consist mostly of vitric dust that settled during lulls, and Layer B is the feather edge of the ash flow. The fall units filled and obscured the caldera, but arcuate and radial fissures outline a 6-km 2 depression. The Novarupta lava dome and its ejecta ring were emplaced later within the depression. At Mt. Katmai, 10 km east of the 1912 vent, a 600-m-deep caldera of similar area also collapsed at about this time, probably owing to hydraulic connection with the venting magma system; but all known ejecta are thought to have erupted at Novarupta. Mingling of three distinctive magmas during the eruption produced an abundance of banded pumice, and mechanical mixing of chilled ejecta resulted in deposits with a wide range of bulk composition. Pumice in the initial fall unit (A) is 100% rhyolite, but fall units atop the ash flow are > 98% dacite; black andesitic scoria is common only in the ash flows and in near-vent air-fall tephra. Pumice counts show the first half of the ash-flow deposit to be 91–98% rhyolite, but progressive increases of dacite and andesite eventually reduced the rhyolitic component to < 2%. The later, rhyolite-poor flows were hotter, less mobile, and widely produced partially welded tuff and vapor-indurated sillar. The main ash flow was too deflated and sluggish 16 km from the vent to surmount a 25-m-high moraine in its path but was diverted around it and continued 5 km down-valley, engulfing and charring trees but not toppling all of them. Thin ash-flow veneers feather 30–40 m up the enclosing valley walls but only where a constriction in the central VTTS locally raised the flow level. In the upper VTTS, the “high sand mark” is not a veneer but a marginal bench formed in thick tuff by differential compaction. Flooding from adjacent glaciers led to phreatic explosions that ejected blocks of tuff more welded than any yet exposed. A cluster of phreatic craters dammed a lake atop the tuff, the breaching of which caused a flood that scoured the ash-flow surface in the central VTTS, transported 50-cm blocks of welded tuff > 20 km to the lowermost VTTS, and deposited 1–8 m of debris there. Rhyolitic ejecta contain only 1–2% phenocrysts but andesite and dacite have 30–45%. Quartz is present and augite absent only in the rhyolite, but all ejecta contain plagioclase, orthopyroxene, titanomagnetite, ilmenite, apatite, and pyrrhotite; rare olivine occurs in the andesite. The zoning ranges of phenocrysts in the rhyolitic and intermediate ejecta do not overlap. New chemical data show the bulk SiO 2 range to be: rhyolite 77 ± 0.6, dacite 66-64.5, and andesite 61.5–58.5%. The dacitic and andesitic ejecta contrast in color and density, and it is not certain whether they form a compositional continuum. Analyses reported by Fenner within the 66–76% SiO 2 range were of banded pumice and lava and of bulk tephra that mechanically fractionated and mixed during flight. Despite the gap of 10% SiO 2 , Fe-Ti-oxide temperatures show a continuous range from rhyolite (805–850°C) through dacite (855–955°C) to andesite (955–990°C). Thermal continuity and isotopic and trace-element data suggest that all were derived from a single magmatic system, whether or not they were physically contiguous before eruption. If the rhyolitic liquid separated from dacitic magma, extraction was so efficient that no dacitic phenocrysts were retained and no bulk compositions in the range 66–76% SiO 2 were created; if it were a partial melt of roof rocks atop an intermediate magma body, then such rocks had no O- or Sr-isotopic contrast with the andesite-dacite magma and clearly did not include the Jurassic arkosic or granitic basement. The presence of Holocene domes of pre-1912 glassy dacite adjacent to the 1912 vent suggest that the 7 km 3 (or more) of high-silica rhyolitic magma (a composition rare in the Aleutian arc) was generated in less than a few thousand years. The 1912 vent is semi-encircled by several andesitic stratocones and is as close to Mageik, Trident, and Griggs volcanoes as it is to Mt. Katmai. The erupted magma probably occupied only shallow levels of an extensive system of injection and storage under a cluster of several stratovolcanoes. Although Quaternary basalt is not known to have erupted here, the intrusion of basaltic magma probably sustains the greater-VTTS magmatic system.

Alaska↗

Eruptive history of Mount Mazama and Crater Lake Caldera, Cascade Range, U.S.A.

New investigations of the geology of Crater Lake National Park necessitate a reinterpretation of the eruptive history of Mount Mazama and of the formation of Crater Lake caldera. Mount Mazama consisted of a glaciated complex of overlapping shields and stratovolcanoes, each of which was probably active for a comparatively short interval. All the Mazama magmas apparently evolved within thermally and compositionally zoned crustal magma reservoirs, which reached their maximum volume and degree of differentiation in the climactic magma chamber ∼ 7000 yr B.P. The history displayed in the caldera walls begins with construction of the andesitic Phantom Cone ∼ 400,000 yr B.P. Subsequently, at least 6 major centers erupted combinations of mafic andesite, andesite, or dacite before initiation of the Wisconsin Glaciation ∼ 75,000 yr B.P. Eruption of andesitic and dacitic lavas from 5 or more discrete centers, as well as an episode of dacitic pyroclastic activity, occurred until ∼ 50,000 yr B.P.; by that time, intermediate lava had been erupted at several short-lived vents. Concurrently, and probably during much of the Pleistocene, basaltic to mafic andesitic monogenetic vents built cinder cones and erupted local lava flows low on the flanks of Mount Mazama. Basaltic magma from one of these vents, Forgotten Crater, intercepted the margin of the zoned intermediate to silicic magmatic system and caused eruption of commingled andesitic and dacitic lava along a radial trend sometime between ∼ 22,000 and ∼ 30,000 yr B.P. Dacitic deposits between 22,000 and 50,000 yr old appear to record emplacement of domes high on the south slope. A line of silicic domes that may be between 22,000 and 30,000 yr old, northeast of and radial to the caldera, and a single dome on the north wall were probably fed by the same developing magma chamber as the dacitic lavas of the Forgotten Crater complex. The dacitic Palisade flow on the northeast wall is ∼ 25,000 yr old. These relatively silicic lavas commonly contain traces of hornblende and record early stages in the development of the climatic magma chamber. Some 15,000 to 40,000 yr were apparently needed for development of the climactic magma chamber, which had begun to leak rhyodacitic magma by 7015 ± 45 yr B.P. Four rhyodacitic lava flows and associated tephras were emplaced from an arcuate array of vents north of the summit of Mount Mazama, during a period of ∼ 200 yr before the climactic eruption. The climactic eruption began 6845 ± 50 yr B.P. with voluminous airfall deposition from a high column, perhaps because ejection of ∼ 4−12 km 3 of magma to form the lava flows and tephras depressurized the top of the system to the point where vesiculation at depth could sustain a Plinian column. Ejecta of this phase issued from a single vent north of the main Mazama edifice but within the area in which the caldera later formed. The Wineglass Welded Tuff of Williams (1942) is the proximal featheredge of thicker ash-flow deposits downslope to the north, northeast, and east of Mount Mazama and was deposited during the single-vent phase, after collapse of the high column, by ash flows that followed topographic depressions. Approximately 30 km 3 of rhyodacitic magma were expelled before collapse of the roof of the magma chamber and inception of caldera formation ended the single-vent phase. Ash flows of the ensuing ring-vent phase erupted from multiple vents as the caldera collapsed. These ash flows surmounted virtually all topographic barriers, caused significant erosion, and produced voluminous deposits zoned from rhyodacite to mafic andesite. The entire climactic eruption and caldera formation were over before the youngest rhyodacitic lava flow had cooled completely, because all the climactic deposits are cut by fumaroles that originated within the underlying lava, and part of the flow oozed down the caldera wall. A total of ∼ 51−59 km 3 of magma was ejected in the precursory and climactic eruptions, and ∼ 40−52 km 3 of Mount Mazama was lost by caldera formation. The spectacular compositional zonation shown by the climactic ejecta — rhyodacite followed by subordinate andesite and mafic andesite — reflects partial emptying of a zoned system, halted when the crystal-rich magma became too viscous for explosive fragmentation. This zonation was probably brought about by convective separation of low-density, evolved magma from underlying mafic magma. Confinement of postclimactic eruptive activity to the caldera attests to continuing existence of the Mazama magmatic system.

Mount Mazama, Crater Lake, Cascade Range↗

Thermal areas on Kilauea and Mauna Loa Volcanoes, Hawaii

Active thermal areas are concentrated in three areas on Mauna Loa and three areas on Kilauea. High-temperature fumaroles (115-362°C) on Mauna Loa are restricted to the summit caldera, whereas high-temperature fumaroles on Kilauea are found in the upper East Rift Zone (Mauna Ulu summit fumaroles, 562°C), middle East Rift Zone (1977 eruptive fissure fumaroles), and in the summit caldera. Solfataric activity that has continued for several decades occurs along border faults of Kilauea caldera and at Sulphur Cone on the southwest rift zone of Mauna Loa. Solfataras that are only a few years old occur along recently active eruptive fissures in the summit caldera and along the rift zones of Kilauea. Steam vents and hot-air cracks also occur at the edges of cooling lava ponds, on the summits of lava shields, along faults and graben fractures, and in diffuse patches that may reflect shallow magmatic intrusions.

Hawaii↗

Sulfur and oxygen isotopic systematics of the 1982 eruptions of El Chichón Volcano, Chiapas, Mexico

Thermometers based on sulfur and oxygen isotopic compositions of anhydrite, pyrrhotite, titanomagnetite, and plagioclase crystals from fresh pumices of the 1982 eruptions of El Chich&oacute;n Volcano indicate a pre-eruption temperature of 810 &plusmn; 40&deg;C, confirming textural evidence that the anhydrite precipitated directly from the melt. The isotopic composition of sulfate leached from fresh ashfall samples shows it to be a mixture of anhydrite microphenocrysts and adsorbed sulfate derived from oxidized sulfur (SO 2 ) in the eruption plume. The leachate data show no evidence for rapid oxidation of significant amounts of H 2 S in the eruption cloud even though the fugacity ratio of H 2 S/SO 2 in the gas phase of the magma was >400. This may indicate kinetic inhibition of H 2 S to SO 2 conversion in the eruption cloud. Prior to eruption, the magma contained an estimated 2.6 wt. % sulfur (as SO 3 ). The estimated &delta; 34 S of the bulk magma is 5.8&permil;. Such a high value may reflect assimilation of 34 S-enriched evaporites or the prior loss of 34 S-depleted H 2 S to a fluid or gas phase during formation of a small prophyry-type hydrothermal system or ore deposit. In either case, the original magma must have been very sulfur rich. It is likely that the initial high sulfur content of the magma and at least some of its 34 S enrichment reflects involvement of subducted volcanogenic massive sulfides deposits during Benioff-zone partial melting. Isotopic data on mineralized, accidental lithic fragments support the possible development of a porphyry-type system at El Chich&oacute;n.

Journal of Volcanology and Geothermal Research↗

Crater lake and post-eruption hydrothermal activity, El Chichón Volcano, Mexico

Explosive eruptions of Volc&aacute;n El Chich&oacute;n in Chiapas, Mexico on March 28 and April 3&ndash;4, 1982 removed 0.2 km 3 of rock to form a 1-km-wide 300-m-deep summit crater. By late April 1982 a lake had begun to form on the crater floor, and by November 1982 it attained a maximum surface area of 1.4 &times; 10 5 m 2 and a volume of 5 &times; 10 6 m 3 . Accumulation of 4&ndash;5 m of rainfall between July and October 1982 largely formed the lake. In January 1983, temperatures of fumaroles on the crater floor and lower crater walls ranged from 98 to 115&deg;C; by October 1983 the maximum temperature of fumarole emissions was 99&deg;C. In January 1983 fumarole gas emissions were greater than 99 vol. % H 2 O with traces of CO 2 , SO 2 , and H 2 S. The water of the lake was a hot ( T = 52&ndash;58&deg;C), acidic (pH = 0.5), dilute solution (34,046 mg L &minus;1 dissolved solids; Cl/S = 20.5). Sediment from the lake contains the same silicate minerals as the rocks of the 1982 pyroclastic deposits, together with less than 1% of elemental sulfur. The composition and temperature of the lake water is attributed to: (1) solution of fumarole emissions; (2) reaction of lake water with hot rocks beneath the lake level; (3) sediments washed into the lake from the crater walls; (4) hydrothermal fluids leaching sediments and formational waters in sedimentary rocks of the basement; (5) evaporation; and (6) precipitation.

Journal of Volcanology and Geothermal Research↗