Geology topics

Katharine V. Cashman

Publications and source records attributed to Katharine V. Cashman.

At least 19 recordsLinked to original sources

USGS2026

Syn-emplacement crystallization of Mauna Loa 2022 lava flows, Hawaiʻi

Advances in near-real-time forecasts of lava flow advance rates, final travel distance, and areal coverage are transforming hazard response, particularly in locations such as Hawaiʻi that experience frequent lava flow activity. The most severe threats are posed by rapidly advancing channelized ʻaʻā flows, such as those that characterized the November–December 2022 eruption of Mauna Loa volcano. To constrain rheological inputs to flow forecasting models during the eruption, samples were collected from the two most persistent lava flows and rapidly assessed for downflow changes in crystallinity and vesicularity. As observed in other channelized lava flows, the lava lost bubbles and cooled during transport, with consequent increases in groundmass crystallinity and bulk viscosity. Lava rheology, however, is controlled by more than bulk bubble and crystal contents; also important are the melt viscosity and the size and shape distributions of bubbles and crystals. Here we further interrogate the sample suite to document a down flow increase of almost two orders of magnitude in melt viscosity alone (a function of cooling combined with compositional change) and assess the extreme anisotropy of the plagioclase crystals as a function of size and flow type (pāhoehoe or ʻaʻā). When combined with previous studies, these data constrain observed relations between surface morphology, flow temperature and crystal content and indicate that remote assessments of morphological transitions and core lava temperature could be used for channel-wide monitoring of rheological evolution for assimilation into flow models.

Hawaii↗

USGS2025

New insights into gas-driven phase segregation in andesitic enclaves from Mt. Mazama (Crater Lake), USA

A key process in active magmatic systems is the “recharge” of deep-sourced mafic magma into cooler, more evolved, and crystal-rich shallow reservoirs; recharge may be the cause of, or response to, eruptive activity. Although compositional evidence for recharge has been extensively documented, physical models of recharge are limited, particularly processes that separate exsolving volatiles and melts from rapidly growing crystals. To improve constraints on phase separation behaviors, we re-examine andesitic enclaves in silicic andesite lava flows of Mt. Mazama (Crater Lake), USA, that provided early evidence of gas-driven filter pressing (Bacon, 1986). 2D and 3D imaging shows that enclaves have a sample-spanning crystal framework that is disrupted by melt patches, indicating that initially deformable crystal networks were subject to early phase reorganization. Small enclaves are poorly vesicular and require early gas loss. Large enclaves have porous cores with angular (diktytaxitic) voids that are well-connected in 3D and denser rinds with isolated pores. Large enclave rinds have similar bulk compositions to small enclaves but their less evolved cores require ~ 20% melt removal. In the large enclave, diktytaxitic core textures and gas fingering structures at the core–rind boundary suggest relatively slow late-stage outward gas migration. Both scaling arguments and evidence of outward gas/melt migration require a resistant rind. Rind formation is best explained by differential cooling and demonstrates the importance of thermal gradients for gas-driven filter pressing. A corollary is a limited time scale of recharge, enclave formation, and vesiculation to produce diktytaxitic textures, suggesting that recharge was (near) synchronous with eruption.

Oregon↗

USGS2023

Decompression and degassing, repressurization, and regassing during cyclic eruptions at Guagua Pichincha volcano, Ecuador, 1999–2001

In 1999–2001, Guagua Pichincha volcano, Ecuador, produced a series of cyclic explosive and effusive eruptions. Rock samples, including dense blocks and pumiceous clasts collected during the eruption sequence, and ballistic bombs later collected from the crater floor, provide information about magma storage, ascent, decompression, degassing, repressurization, and regassing prior to eruption. Pairs of Fe-Ti oxides indicate equilibrium within 1.2–1.5 log units above the NNO oxidation buffer and equilibrium temperatures from 805 to 905 °C. Melt inclusions record H 2 O contents of 2.7–4.6 wt% and CO 2 contents (uncorrected for CO 2 segregation into bubbles) from 19 to 310 ppm. Minimum melt inclusion saturation pressures fall between 69 and 168 MPa, or equilibration depths of 2.8 and 6.8 km, the lower end of which is coincident with the maximum inferred equilibration depths for the most vesicular breadcrust bombs sampled. Amphibole phenocrysts lack breakdown rims (except for one sample) and plagioclase phenocrysts have abundant oscillatory compositional zones. Plagioclase areal microlite number densities ( N a ) range over less than one order of magnitude (8.9×10 3 –8.7×10 4 mm -2 ) among all samples, with the exception of a dense, low crystallinity sample ( N a = 3.0×10 3 mm −2 ) and a pumiceous sample erupted on 17 December 1999 ( N a = 1.7×10 3 mm −2 ). Plagioclase microlite shapes include tabular, hopper, and swallowtail forms. Taken together, the relatively high plagioclase microlite number densities, the high number of oscillatory zones in plagioclase phenocrysts, the presence of CO 2 in groundmass glass, seismicity, and time-varying tilt cycles provide a picture of sudden evacuation of magma residing at different levels in the shallow conduit. Explosive eruptions punctuate inter-eruptive repose periods marked by time-varying rates of degassing (volatile fluxing) and re-pressurization. Shallow residence time in the conduit was sufficient to allow precipitation of silica-phase in the groundmass, but insufficient to allow breakdown of hornblende phenocrysts, with the one exception of the final dome sample from 2000, which has the longest preceding repose time. These results support a model of cyclic pressure cycling, volatile exsolution and regassing, and magma decompression decoupled from ascent.

Guagua Pichincha Volcano↗

USGS2022

Can lava flow like water? Assessing applications of critical flow theory to channelized basaltic lava flows

Flowing lava and water have dramatically different physical properties but can form similar hydraulic structures, including undular hydraulic jumps, or standing wave trains. In water flows, undular hydraulic jumps are evidence of critical flow (Froude number ∼1) and open-channel hydraulic theory provides a powerful tool for estimating flow depth and velocity. Monitoring these parameters in an active lava channel is inherently challenging, but essential for calculating lava discharge (effusion rate), a primary control on the rate of flow front advance and ultimate flow runout distance. We analyze undular hydraulic jumps in both water and lava flows to assess the conditions under which they form and, by extension, the potential use of critical flow theory to estimate, in real time, lava flow velocity, depth, and discharge. Experimental data for water flows show that these structures mark the transition from supercritical to subcritical flow. Undular hydraulic jumps in the near-vent lava channel of the 2018 lower East Rift Zone eruption of Kīlauea, Hawaiʻi also reflect critical flow conditions; their wavelengths scale with flow depth and velocity, consistent with hydraulic theory. Calculated lava effusion rates are similar to estimates made using more traditional approaches (Jeffreys', 1925, https://doi.org/10.1080/14786442508634662 , equation based on lava viscosity, density, and channel slope) and with lava volumes derived from topographic-change mapping. From this we conclude that critical flow phenomena show great potential to track flow dynamics and inform hazard assessment for a wide range of geophysical fluids.

Journal of Geophysical Research - Earth Surface↗

USGS2022

Modelling the transport and deposition of ash following a magnitude 7 eruption: The distal Mazama tephra

Volcanic ash transport and dispersion models (VATDMs) are necessary for forecasting tephra dispersal during volcanic eruptions and are a useful tool for estimating the eruption source parameters (ESPs) of prehistoric eruptions. Here we use Ash3D, an Eulerian VATDM, to simulate the tephra deposition from the ~ 7.7 ka climactic eruption of Mount Mazama. We investigate how best to apply a VATDM using the ESPs characteristic of a large magnitude eruption (M ≥ 7). We simplify the approach to focus on the distal deposit as if it were formed by a single phase of Plinian activity. Our results demonstrate that it is possible to use modern wind profiles to simulate the tephra dispersal from a prehistoric eruption; however, this introduces an inherent uncertainty to the subsequent simulations where we explore different ESPs. We show, using the well-documented distal Mazama tephra, that lateral umbrella cloud spreading, rather than advection–diffusion alone, must be included in the VATDM to reproduce the width of the isopachs. In addition, the Ash3D particle size distribution must be modified to simulate the transport and deposition of distal fine-grained (< 125 µm) Mazama ash. With these modifications, the Ash3D simulations reproduce the thickness and grain size of the Mazama tephra deposit. Based on our simulations, however, we conclude that the exact relationship between mass eruption rate and the scale of umbrella cloud spreading remains unresolved. Furthermore, for ground-based grain size distributions to be input directly into Ash3D, further research is required into the atmospheric and particle processes that control the settling behaviour of fine volcanic ash.

Bulletin of Volcanology↗

USGS2020

Understanding the storage conditions and fluctuating eruption style of a young monogenetic volcano: Blue Lake crater (<3 ka), High Cascades, Oregon

Blue Lake crater (<3 ka) is monogenetic volcano that produced one of the youngest eruptions in the central Oregon Cascades. Understanding monogenetic volcano behavior – from storage through eruption – is imperative in planning for future eruptions. Here we combine physical volcanology and geochemistry to determine the pre-eruptive storage conditions, ascent rate, eruption style, and deposit distribution of this young eruption. We find that the eruption of Blue Lake was initially phreatomagmatic, producing lithic-rich fall deposits and thin surge deposits and excavating the maar crater, before transitioning rapidly to a final voluminous magmatic-volatile driven explosive eruption. The mapped fall deposit has an estimated volume of 3.9 × 10 7 m 3 (2.2 × 10 7 m 3 DRE) which suggests a VEI of 3. Although similar in magnitude (as measured by fall deposit volume) to many other recent cinder cone eruptions in the Cascades, the Blue Lake crater eruption lacks an effusive phase. The absence of lava flows may reflect the lack of evidence for syn-eruptive magma storage at shallow levels. Indeed, corrected volatile contents of olivine-hosted melt inclusions (2.9–4.2 wt% H 2 O, 910–1330 ppm CO 2 ) are strikingly uniform and indicate storage and crystallization at a restricted pressure range (average ~ 235 MPa), equating to a depth of ~8.6 km. Melt inclusion geochemistry indicates that the basaltic andesite magma cooled and crystallized ~25% during storage at this pressure. Crystals in the Blue Lake magma show evidence of mixing with, or entrainment in, a more evolved magma. Feldspar crystals have large An-rich cores (An 80–85 ) and abrupt An-poor rims (An 60–70 ); olivine crystals have large, broad cores (~Fo 82–84 ) and thin rims with lower Fo and NiO contents. Diffusion modeling of olivine zoning suggests that an intrusion event occurred ~10–60 days prior to eruption. Diffusive loss of H + from melt inclusions was minimal (<1.3 wt% H 2 O) during magma ascent, from which we calculate minimum ascent times from 235 MPa of <1 day. Many inclusions indicate ascent times of <3 h, corresponding to ascent rates of ~1 to >13 m/s. This study illustrates the pre-eruptive and eruptive complexities of monogenetic volcanoes and highlights the minimal warning that may precede future eruptions.

Oregon↗

USGS2017

Benchmarking computational fluid dynamics models of lava flow simulation for hazard assessment, forecasting, and risk management

Numerical simulations of lava flow emplacement are valuable for assessing lava flow hazards, forecasting active flows, designing flow mitigation measures, interpreting past eruptions, and understanding the controls on lava flow behavior. Existing lava flow models vary in simplifying assumptions, physics, dimensionality, and the degree to which they have been validated against analytical solutions, experiments, and natural observations. In order to assess existing models and guide the development of new codes, we conduct a benchmarking study of computational fluid dynamics (CFD) models for lava flow emplacement, including VolcFlow, OpenFOAM, FLOW-3D, COMSOL, and MOLASSES. We model viscous, cooling, and solidifying flows over horizontal planes, sloping surfaces, and into topographic obstacles. We compare model results to physical observations made during well-controlled analogue and molten basalt experiments, and to analytical theory when available. Overall, the models accurately simulate viscous flow with some variability in flow thickness where flows intersect obstacles. OpenFOAM, COMSOL, and FLOW-3D can each reproduce experimental measurements of cooling viscous flows, and OpenFOAM and FLOW-3D simulations with temperature-dependent rheology match results from molten basalt experiments. We assess the goodness-of-fit of the simulation results and the computational cost. Our results guide the selection of numerical simulation codes for different applications, including inferring emplacement conditions of past lava flows, modeling the temporal evolution of ongoing flows during eruption, and probabilistic assessment of lava flow hazard prior to eruption. Finally, we outline potential experiments and desired key observational data from future flows that would extend existing benchmarking data sets.

Journal of Applied Volcanology↗

USGS2016

Holocene volcanism of the upper McKenzie River catchment, central Oregon Cascades, USA

To assess the complexity of eruptive activity within mafic volcanic fields, we present a detailed geologic investigation of Holocene volcanism in the upper McKenzie River catchment in the central Oregon Cascades, United States. We focus on the Sand Mountain volcanic field, which covers 76 km 2 and consists of 23 vents, associated tephra deposits, and lava fields. We find that the Sand Mountain volcanic field was active for a few decades around 3 ka and involved at least 13 eruptive units. Despite the small total volume erupted (∼1 km 3 dense rock equivalent [DRE]), Sand Mountain volcanic field lava geochemistry indicates that erupted magmas were derived from at least two, and likely three, different magma sources. Single units erupted from one or more vents, and field data provide evidence of both vent migration and reoccupation. Overall, our study shows that mafic volcanism was clustered in space and time, involved both explosive and effusive behavior, and tapped several magma sources. These observations provide important insights on possible future hazards from mafic volcanism in the central Oregon Cascades.

Geological Society of America Bulletin↗

USGS2015

Diverting lava flows in the lab

Recent volcanic eruptions in Hawai'i, Iceland and Cape Verde highlight the challenges of mitigating hazards when lava flows threaten infrastructure. Diversion barriers are the most common form of intervention, but historical attempts to divert lava flows have met with mixed success and there has been little systematic analysis of optimal barrier design. We examine the interaction of viscous flows of syrup and molten basalt with barriers in the laboratory. We find that flows thicken immediately upslope of an obstacle, forming a localized bow wave that can overtop barriers. Larger bow waves are generated by faster flows and by obstacles oriented at a high angle to the flow direction. The geometry of barriers also influences flow behaviour. Barriers designed to split or dam flows will slow flow advance, but cause the flow to widen, whereas oblique barriers can effectively divert flows, but may also accelerate flow advance. We argue that to be successful, mitigation of lava-flow hazards must incorporate the dynamics of lava flow–obstacle interactions into barrier design. The same generalizations apply to the effect of natural topographic features on flow geometry and advance rates .

Nature Geoscience↗

USGS2014

The dynamics of Hawaiian-style eruptions: A century of study

This chapter, prepared in celebration of the Hawaiian Volcano Observatoryʼs centennial, provides a historical lens through which to view modern paradigms of Hawaiian-style eruption dynamics. The models presented here draw heavily from observations, monitoring, and experiments conducted on Kīlauea Volcano, which, as the site of frequent and accessible eruptions, has attracted scientists from around the globe. Long-lived eruptions in particular—Halema‘uma‘u 1907–24, Kīlauea Iki 1959, Mauna Ulu 1969–74, Pu‘u ‘Ō‘ō-Kupaianaha 1983–present, and Halema‘uma‘u 2008–present—have offered incomparable opportunities to conceptualize and constrain theoretical models with multidisciplinary data and to field-test model results. The central theme in our retrospective is the interplay of magmatic gas and near-liquidus basaltic melt. A century of study has shown that gas exsolution facilitates basaltic dike propagation; volatile solubility and vesiculation kinetics influence magma-rise rates and fragmentation depths; bubble interactions and gas-melt decoupling modulate magma rheology, eruption intensity, and plume dynamics; and pyroclast outgassing controls characteristics of eruption deposits. Looking to the future, we anticipate research leading to a better understanding of how eruptive activity is influenced by volatiles, including the physics of mixed CO 2 -H 2 O degassing, gas segregation in nonuniform conduits, and vaporization of external H 2 O during magma ascent.

Hawaii↗

USGS2013

Compaction and gas loss in welded pyroclastic deposits as revealed by porosity, permeability, and electrical conductivity measurements of the Shevlin Park Tuff

Pyroclastic flows produced by large volcanic eruptions commonly densify after emplacement. Processes of gas escape, compaction, and welding in pyroclastic-flow deposits are controlled by the physical and thermal properties of constituent material. Through measurements of matrix porosity, permeability, and electrical conductivity, we provide a framework for understanding the evolution of pore structure during these processes. Using data from the Shevlin Park Tuff in central Oregon, United States, and from the literature, we find that over a porosity range of 0%–70%, matrix permeability varies by almost 10 orders of magnitude (from 10 –20 to 10 –11 m 2 ), with over three orders of magnitude variation at any given porosity. Part of the variation at a given porosity is due to permeability anisotropy, where oriented core samples indicate higher permeabilities parallel to foliation (horizontally) than perpendicular to foliation (vertically). This suggests that pore space is flattened during compaction, creating anisotropic crack-like networks, a geometry that is supported by electrical conductivity measurements. We find that the power law equation: k 1 = 1.3 × 10 –21 × ϕ 5.2 provides the best approximation of dominant horizontal gas loss, where k 1 = permeability, and ϕ = porosity. Application of Kozeny-Carman fluid-flow approximations suggests that permeability in the Shevlin Park Tuff is controlled by crack- or disk-like pore apertures with minimum widths of 0.3 and 7.5 μm. We find that matrix permeability limits compaction over short times, but deformation is then controlled by competition among cooling, compaction, water resorption, and permeable gas escape. These competing processes control the potential for development of overpressure (and secondary explosions) and the degree of welding in the deposit, processes that are applicable to viscous densification of volcanic deposits in general. Further, the general relationships among porosity, permeability, and pore geometry are relevant for flow of any fluid through an ignimbritic host.

Oregon↗

USGS2012

Faulting within the Mount St. Helens conduit and implications for volcanic earthquakes

The 2004–2008 eruption of Mount St. Helens produced seven dacite spines mantled by cataclastic fault rocks, comprising an outer fault core and an inner damage zone. These fault rocks provide remarkable insights into the mechanical processes that accompany extrusion of degassed magma, insights that are useful in forecasting dome-forming eruptions. The outermost part of the fault core consists of finely comminuted fault gouge that is host to 1- to 3-mm-thick layers of extremely fine-grained slickenside-bearing ultracataclasite. Interior to the fault core, there is an ∼2-m-thick damage zone composed of cataclastic breccia and sheared dacite, and interior to the damage zone, there is massive to flow-banded dacite lava of the spine interior. Structures and microtextures indicate entirely brittle deformation, including rock breakage, tensional dilation, shearing, grain flow, and microfaulting, as well as gas and fluid migration through intergranular pores and fractures in the damage zone. Slickenside lineations and consistent orientations of Riedel shears indicate upward shear of the extruding spines against adjacent conduit wall rocks. Paleomagnetic directions, demagnetization paths, oxide mineralogy, and petrology indicate that cataclasis took place within dacite in a solidified steeply dipping volcanic conduit at temperatures above 500 °C. Low water content of matrix glass is consistent with brittle behavior at these relatively high temperatures, and the presence of tridymite indicates solidification depths of <1 km. Cataclasis was coincident with the eruption’s seismogenic zone at <1.5 km. More than a million small and low-frequency “drumbeat” earthquakes with coda magnitudes (M d ) <2.0 and frequencies <5 Hz occurred during the 2004–2008 eruption. Our field data provide a means with which to estimate slip-patch dimensions for shear planes and to compare these with estimates of slip patches based on seismic moments and shear moduli for dacite rock and granular fault gouge. Based on these comparisons, we find that aseismic creep is achieved by micron-scale displacements on Riedel shears and by granular flow, whereas the drumbeat earthquakes require millimeter to centimeter displacements on relatively large (e.g., ∼1000 m 2 ) slip patches, possibly along observed extensive principal shear zones within the fault core but probably not along the smaller Riedel shears. Although our field and structural data are compatible with stick-slip models, they do not rule out seismic and infrasound models that call on resonance of steam-filled fractures to generate the drumbeat earthquakes. We suggest that stick-slip and gas release processes may be coupled, and that regardless of the source mechanism, the distinctive drumbeat earthquakes are proving to be an effective precursor for dome-forming eruptions. Our data document a continuous cycle of deformation along the conduit margins beginning with episodes of fracture in the damage zone and followed by transfer of motion to the fault core. We illustrate the cycle of deformation using a hypothetical cross section of the Mount St. Helens conduit, extending from the surface to the depth of magmatic solidification.

Washington↗

USGS2012

Estimating rates of decompression from textures of erupted ash particles produced by 1999-2006 eruptions of Tungurahua volcano, Ecuador

Persistent low- to moderate-level eruptive activity of andesitic volcanoes is difficult to monitor because small changes in magma supply rates may cause abrupt transitions in eruptive style. As direct measurement of magma supply is not possible, robust techniques for indirect measurements must be developed. Here we demonstrate that crystal textures of ash particles from 1999 to 2006 Vulcanian and Strombolian eruptions of Tungurahua volcano, Ecuador, provide quantitative information about the dynamics of magma ascent and eruption that is difficult to obtain from other monitoring approaches. We show that the crystallinity of erupted ash particles is controlled by the magma supply rate (MSR); ash erupted during periods of high magma supply is substantially less crystalline than during periods of low magma supply. This correlation is most easily explained by efficient degassing at very low pressures (<<50 MPa) and degassing-driven crystallization controlled by the time available prior to eruption. Our data also suggest that the observed transition from intermittent Vulcanian explosions at low MSR to more continuous periods of Strombolian eruptions and lava fountains at high MSR can be explained by the rise of bubbles through (Strombolian) or trapping of bubbles beneath (Vulcanian) vent-capping, variably viscous (and crystalline) magma.

Mt. Tungurahua↗

USGS2012

Tracking lava flow emplacement on the east rift zone of Kilauea, Hawai’i with synthetic aperture radar (SAR) coherence

Lava flow mapping is both an essential component of volcano monitoring and a valuable tool for investigating lava flow behavior. Although maps are traditionally created through field surveys, remote sensing allows an extraordinary view of active lava flows while avoiding the difficulties of mapping on location. Synthetic aperture radar (SAR) imagery, in particular, can detect changes in a flow field by comparing two images collected at different times with SAR coherence. New lava flows radically alter the scattering properties of the surface, making the radar signal decorrelated in SAR coherence images. We describe a new technique, SAR Coherence Mapping (SCM), to map lava flows automatically from coherence images independent of look angle or satellite path. We use this approach to map lava flow emplacement during the Pu‘u ‘Ō‘ō-Kupaianaha eruption at Kīlauea, Hawai‘i. The resulting flow maps correspond well with field mapping and better resolve the internal structure of surface flows, as well as the locations of active flow paths. However, the SCM technique is only moderately successful at mapping flows that enter vegetation, which is also often decorrelated between successive SAR images. Along with measurements of planform morphology, we are able to show that the length of time a flow stays decorrelated after initial emplacement is linearly related to the flow thickness. Finally, we use interferograms obtained after flow surfaces become correlated to show that persistent decorrelation is caused by post-emplacement flow subsidence.

Hawai'i↗

USGS2009

The length of channelized lava flows: Insight from the 1859 eruption of Mauna Loa Volcano, Hawai‘i

The 1859 eruption of Mauna Loa Volcano, Hawai'i, produced paired 'a'ā and pāhoehoe flows of exceptional length (51 km). The 'a'ā flow field is distinguished by a long (> 36 km) and well-defined pāhoehoe-lined channel, indicating that channelized lava remained fluid to great distances from the vent. The 1859 eruption was further unusual in initiating at a radial vent on the volcano's northwest flank, instead of along the well-defined rift zone that has been the source of most historic activity. As such, it presents an opportunity both to examine controls on the emplacement of long lava channels and to assess hazards posed by future flank eruptions of Mauna Loa. Here we combine evidence from historical chronicles with analysis of bulk compositions, glass geothermometry, and microlite textures of samples collected along the 1859 lava flows to constrain eruption and flow emplacement conditions. The bulk compositions of samples from the 'a'ā and pāhoehoe flow fields are bimodally distributed and indicate tapping of two discrete magma bodies during eruption. Samples from the pāhoehoe flow field have bulk compositions similar to those of historically-erupted lavas (< 8 wt.% MgO); lava that fed the 'a'ā channel is more primitive (> 8 wt.% MgO), nearly aphyric, and was erupted at high temperatures (1194–1216 °C). We suggest that the physical properties of proximal channel-fed lava (i.e., high-temperature, low crystallinity, and low bulk viscosity) promoted both rapid flow advance and development of long pāhoehoe-lined channels. Critical for the latter was the large temperature decrease (~ 50 °C) required to reach the point at which plagioclase and pyroxene started to crystallize; the importance of phase constraints are emphasized by our difficulty in replicating patterns of cooling and crystallization recorded by high-temperature field samples using common models of flow emplacement. Placement of the 1859 eruption within the context of historic activity at Mauna Loa suggests that the formation of radial vents and eruptions of high-temperature magma may not only be linked, but may also be a consequence of periods of high magma supply (e.g., 1843–1877). Flank eruptions could therefore warrant special consideration in models and hazard mitigation efforts.

Hawai'i↗

USGS2008

From dome to dust: shallow crystallization and fragmentation of conduit magma during the 2004-2006 dome extrusion of Mount St. Helens, Washington

An unusual feature of the 2004-6 eruptive activity of Mount St. Helens has been the continuous growth of successive spines that are mantled by thick fault gouge. Fault gouge formation requires, first, solidification of ascending magma within the conduit, then brittle fragmentation and cataclastic flow. We document these processes through field relations, hand samples, and thin-section textures. Field observations show that the gouge zone is typically 1-3 m thick and that it includes cataclasite and, locally, breccia in addition to unconsolidated (true) gouge. The gouge contains multiple slickenside sets oriented subparallel to each other and to the striation direction, as well as surface striations parallel to extrusion direction. Hand specimens show the cataclasite and gouge to be composed of a wide size range of broken dome and wallrock fragments. This grain-size heterogeneity is even more pronounced in thin section, where individual samples contain fragments that span more than four orders of magnitude in size (from more than 10 to less than 10 -3 mm). Textures preserved within the gouge zone provide evidence of different processes operating in time and space. Most individual fragments are holocrystalline, suggesting that crystallization of the ascending magma preceded gouge formation. Cataclasite samples preserve a wide range of clast sizes; pronounced rounding of many clasts indicates extensive abrasion during transport. Within the gouge, crystals and lava fragments adjacent to finely comminuted shear zones (slickensides) are shattered into small, angular fragments that are either preserved in place, with little disruption, or incorporated into shear trains, creating a well-developed foliation. Together, evidence of initial grain shattering, followed by shear, grinding, and wear, suggests extensive transport distances (large strains). Textural transitions are often abrupt, indicating extreme shear localization during transport. Comparison of groundmass textures from dome lavas and fault gouge further suggests that brittle fracture was confined to the upper 400-500 m of the conduit. Observed magma extrusion (ascent) rates of ~7 m/d (8×10 -5 m/s) permit several weeks for magma ascent from ~1,000 m (where groundmass crystallization becomes important) to ~500 m (where solidification nears completion). Brittle fracture, cataclastic flow, and shear localization (slickenside formation) probably dominated in the upper 500 m of the conduit. Comparison of eruptive conditions during the 2004-6 activity at Mount St. Helens with those of other spine-forming eruptions suggests that magma ascent rates of about 10 -4 m/s or less allow sufficient degassing and crystallization within the conduit to form large volcanic spines of intermediate composition (andesite to dacite). Solidification deep within the conduit, in turn, requires transport of the solid plug over long distances (hundreds of meters); resultant large strains are responsible for extensive brittle breakage and development of thick gouge zones. Moreover, similarities between gouge textures and those of ash emitted by explosions from spine margins indicate that fault gouge is the origin for the ash. As the comminution and generation of ash-sized particles was clearly a multistep process, this observation suggests that fragmentation preceded, rather than accompanied, these explosions.

Washington↗

USGS2008

Petrology of the 2004-2006 Mount St. Helens lava dome -- implications for magmatic plumbing and eruption triggering

Eighteen years after dome-forming eruptions ended in 1986, and with little warning, Mount St. Helens began to erupt again in October 2004. During the ensuing two years, the volcano extruded more than 80×10 6 m 3 of gas-poor, crystal-rich dacite lava. The 2004-6 dacite is remarkably uniform in bulk-rock composition and, at 65 percent SiO 2 , among the richest in silica and most depleted in incompatible elements of the magmas erupted at Mount St. Helens during the past 500 years. Since shortly after the first spine of lava appeared, samples have been collected using a steel box dredge (“Jaws”) suspended 20-35 m below a helicopter and, occasionally, by hand sampling. As of the spring of 2006, 25 age-controlled samples have been collected from the seven spines of the new lava dome. Samples were obtained from both the interiors of spines and from their carapaces, which are composed of fault gouge and cataclasite 1-2 m thick. The dacite lava is crystal rich, with 40-50 percent phenocrysts. The groundmass is extensively crystallized to a cotectic assemblage of quartz, tridymite, and Na- and K-rich feldspar microlites, raising the total crystal content to more than 80 percent on a vesicle-free basis in all but the earliest erupted samples. Early samples and those collected from near the spine margin are more glassy and vesicular that those collected later and from the interior of the spines. Oxide thermobarometer determinations for the earliest erupted samples we collected cluster at temperatures of approximately 850°C and at an oxygen fugacity one log unit above the nickel-nickel oxide (NNO) buffer curve. In contrast, samples from relatively glass-poor samples erupted in late 2004 and early 2005 have zoned oxides with apparent temperatures that range to greater than 950°C. The higher temperatures in these microlite-rich rocks are attributed to latent heat evolved during extensive and rapid groundmass crystallization. Low volatile contents of matrix glasses and presence of tridymite and quartz in the high-silica rhyolite matrix glass indicate extensive shallow (<1 km) crystallization of the matrix, driven by degassing of water and solidifying the magma below the level of the vent. The mode of eruption of the dacite as a series of fault-gouge-mantled spines is explained by this process of extensive subvent degassing and solidification. Although the dacite from this eruption is more silica rich than 1980-86 dome rocks, most major and trace element concentrations of the 1980-86 and 2004-6 magma batches are similar, and magmatic gas emissions have been low and have had similar ratios to those of the 1980s, raising the possibility that the magma might be residual from the 1980–86 reservoir. However, titanium and chromium are enriched slightly relative to the most recent 1980-86 and Goat Rocks (A.D. 1800-1857) eruptive cycles, and heavy rare-earth-element abundances are slightly depleted relative to those erupted during the past 500 years at Mount St. Helens. These data suggest either addition of new gas-poor dacite magma or tapping of a region of the preexisting reservoir that was not erupted previously. A relatively low pressure of last phenocryst growth suggests that the magma was derived from near the apex of the Mount St. Helens magma reservoir at a depth of about 5 km. Viewed in the context of seismic, deformation, and gas-emission data, the petrologic and geochemical data can be explained by ascent of a geochemically distinct batch of magma into the apex of the reservoir during the period 1987-97, followed by upward movement of magma into a new conduit beginning in late September 2004. The question of new versus residual magma has implications for the long-term eruptive behavior of Mount St. Helens, because arrival of a new batch of dacitic magma from the deep crust could herald the beginning of a new long-term cycle of eruptive activity. It is also important to our understanding of what triggered the eruption and its future course. Two hypotheses for triggering are considered: (1) top-down fracturing related to the shallow groundwater system and (2) an increase in reservoir pressure brought about by recent magmatic replenishment. With respect to the future course of the eruption, similarities between textures and character of eruption of the 2004-6 dome and the long-duration (greater than 100 years) pre-1980 summit dome, along with the low eruptive rate of the current eruption, suggest that the eruption could continue sluggishly or intermittently for years to come.

Washington↗

USGS2008

Evolving magma storage conditions beneath Mount St. Helens inferred from chemical variations in melt inclusions from the 1980-1986 and current (2004-2006) eruptions

Major element, trace element, and volatile concentrations in 187 glassy melt inclusions and 25 groundmass glasses from the 1980-86 eruption of Mount St. Helens are presented, together with 103 analyses of touching FE-Ti oxide pairs from the same samples. These data are used to evaluate the temporal evolution of the magmatic plumbing system beneath the volcano during 1980-86 and so provide a framework in which to interpret analyses of melt inclusions from the current (2004-2006) eruption. Major and trace element concentrations of all melt inclusions lie at the high SiO 2 end of the data array defined by eruptive products of the late Quaternary age from Mount St. Helens. For several major and trace elements, the glasses define a trend that is oblique to the whole-rock trend, indicating that different mineral assemblages were responsible for the two trends. The whole-rock trend can be ascribed to differentiation of hydrous basaltic parents in a deep-seated magma reservoir, probably at depths great enough to stabilize garnet. In contrast, the glass trends were generated by closed-system crystallization of the phenocryst and microlite mineral assemblages at low pressures. The dissolved H 2 O content of the melt inclusions from 1980-86, as measured by the ion microprobe, ranges from 0 to 6.7 wt. percent, with the highest values obtained from the plinian phase of May 18, 1980. Water contents decrease with increasing SiO 2 , consistent with decompression-driven crystallization. Preliminary data for dissolved CO 2 in melt inclusions from the May 18 plinian phase from August 7, 1980, indicate that X H 2 O in a vapor phase was approximately constant at 0.80, irrespective of H 2 O content, suggestive of closed-system degassing with a high bubble fraction or gas streaming through the subvolcanic system. Temperature and f O 2 estimates for touching Fe-Ti oxides show evidence for heating during crystallization owing to release of latent heat. Consequently, magmas with the highest microlite crystallinities record the highest temperatures. Magmas also become progressively reduced during ascent and degassing, probably as a result of redox equilibria between exsolving S-bearing gases and magmas. The lowest temperature oxides have f O 2 ≈ NNO, similar to high-temperature fumarole gases from the volcano. The temperature and f O 2 of the magma tapped by the plinian phase of May 18, 1980, are 870-875°C and NNO+0.8, respectively. The dissolved volatile contents of the melt inclusions have been used to calculate sealing pressures; that is, the pressure at which chemical exchange between inclusion and matrix melt ceased. These are greatest for the May 18 plinian magma (120 to 320 MPa); lower pressures are recorded by samples of the preplinian cryptodome and by all post-May 18 magmas. Magma crystallinity, calculated from melt-inclusion Rb contents, is negatively correlated with sealing pressure, consistent with decompression crystallization. Elevated contents of Li in melt inclusions from the cryptodome and post-May 18 samples are consistent with transfer of Li in a magmatic vapor phase from deeper parts of the magma system to magma stored at shallower levels. The Li enrichment attains its maximum extent at ~150 MPa, which is ascribed to separation of a single vapor phase into H 2 O-rich gas and dense Li-rich brine at the top of the magma column. There are striking correlations between melt-inclusion chemistry and monitoring data for the 1980-86 eruption. Dissolved SO 2 contents of melt inclusions from any given event, multiplied by the mass of magma erupted during that event, correlate with the measured flux of SO 2 at the surface, suggesting that magma degassing and melt-inclusion sealing are closely related in time and space. Textural and chemical evidence indicates that melt inclusions became effectively sealed (physically or kinetically) shortly before eruption. Thus by converting pressure to depth using a density model and edifice-loading algorithm for the volcano, changing depths of magma extraction with time can be tracked and compared to the seismic record. The plinian eruption of May 18, 1980, involved magma stored 5-11 km below sea level; this is inferred to be the subvolcanic magma chamber. The preceding eruptions, including the May 18, 1980, blast, involved magma withdrawal from the cryptodome and conduit down to 5 km below sea level. Subsequent 1980 eruptions tapped magma down to depths of ≤10 km below sea level. Tapping of magma stored deeper than 2 km below sea level stopped abruptly at the end of 1980, coincident with the onset of extensive shallow seismicity and a change from explosive to effusive eruption style from 1981 to 1986. Overall, the 1980-86 eruption is consistent with the evisceration of a thin, vertically extensive body of magma extending from 5 to at least 11 km below sea level and connected to the surface by a thin conduit. In the absence of sustained high magma-supply rates from depth, decompression crystallization of magma ascending through the system leads eventually to plugging of the conduit. The current eruption of Mount St. Helens shares some similarities with the 1981-86 dome-building phase of the previous eruption, in that there is extensive shallow seismicity and extrusion of highly crystalline material in the form of a sequence of flows and spines. Melt inclusions from the current eruption have low H 2 O contents, consistent with magma extraction from shallow depths. Highly enriched Li in melt inclusions suggests that vapor transport of Li is a characteristic feature of Mount St. Helens. Melt inclusions from the current eruption have subtly different trace-element chemistry from all but one of the 1980-86 melt inclusions, with steeper rareearth-element (REE) patterns and low U, Th, and high-fieldstrength elements (HFSE), indicating addition of a new melt component to the magma system. It is anticipated that increasing involvement of the new melt component will be evident as the current eruption proceeds.

Washington↗