Geology topics

Paul J. Wallace

Publications and source records attributed to Paul J. Wallace.

10 recordsLinked to original sources

USGS2026

Magmatic source of the opening phase of Kīlauea’s 2018 Lower East Rift Zone eruption

The 2018 eruption of Kīlauea volcano in its Lower East Rift Zone began with the discharge of evolved high-Ti basalt as weak lava fountains and short, slow-moving lava flows. The lavas were quickly geochemically recognized as being derived from magmas stored within the rift zone and remobilized by a new intrusion, a sequence that is common at Kīlauea. This initial phase of the 2018 eruption, referred to as phase 1a, lasted for 6 days and was followed by extrusion of mixed magma after a 3-day pause. Even though remobilization of older rift zone magmas is common within Kīlauea’s rift zones, it is difficult to determine which past intrusion(s) may have initially emplaced those stored magmas. This difficulty stems from the tendency for Kīlauea magmas to follow very similar differentiation paths without significant variations in major, minor, or even trace element chemistry. We investigate possible magma sources for the lavas erupted during phase 1a of the 2018 eruption using whole-rock, mineral, and glass major and trace element compositions from historical East Rift Zone eruptions with adjacent fissures. We consider two primary hypotheses for the phase 1a source: magmas associated with the 1955 Lower East Rift Zone eruption or the nine eruptions in the Middle and Upper East Rift Zone during the 1960s. Our results suggest that magma associated with the earliest phases of Kīlauea’s 1955 eruption was the most likely source of the 2018 phase 1a remobilized magma. We determine volatile saturation pressures from melt inclusion chemistry and find similar storage depths for the 2018 phase 1a and early 1955 magmas. The phase 1a and early 1955 lavas are nearly indistinguishable in all of the compositional criteria considered, implying that the leftover 1955 magma body barely cooled and differentiated in the 63 years between eruptions (cooling rates of ~0.1 °C/year). This study sheds light on the potential for protracted storage of eruptible magmas in rift zones at Kīlauea, and highlights some of the challenges and solutions to identifying genetic relationships between magmas at Kīlauea.

Hawaii↗

USGS2024

Insights into magma storage depths and eruption controls at Kīlauea Volcano during explosive and effusive periods of the past 500 years based on melt and fluid inclusions

Kīlauea Volcano experiences centuries-long cycles of explosive and effusive eruptive behavior, but the relation, if any, between these eruptive styles and changing conditions in the magma plumbing system remains poorly known. We analyze olivine-hosted melt and fluid inclusions to determine magma storage depths during the explosive-era Keanakākoʻi Tephra eruptions (∼1500–1840 CE) and compare these results to modern effusive-era Kīlauea eruptions (1959 Kīlauea Iki, 1960 Kapoho, 2018 lower East Rift Zone). We find that shallow (1–3 km) magma storage has persisted for centuries at Kīlauea, spanning both explosive and effusive periods. In contrast, mid-crustal zones of magma storage shallowed over time, from 5 to 8 km during the Keanakākoʻi sequence to 3–5 km during the modern effusive period. Melt and fluid inclusions in high-forsterite olivine (Fo 86–89 ) trapped at shallow depths indicate that high-temperature magmas (1200 to ∼1300 °C) commonly reach depths of ≤3 km. CO 2 -rich fluid inclusions are present in olivine from all investigated Kīlauea eruptions but are larger and much more abundant in Keanakākoʻi units, which we interpret as indicating that a greater volume fraction of exsolved CO 2 -rich fluid was present in pre-eruptive Keanakākoʻi melts. Increased amounts of CO 2 -rich fluids in the Keanakākoʻi-era magmas would have increased magma buoyancy and driven rapid magma ascent, thereby increasing eruption energy and enhancing near-surface magma-water interactions compared to the current effusive period.

Hawaii↗

USGS2021

The petrologic and degassing behavior of sulfur and other magmatic volatiles from the 2018 eruption of Kīlauea, Hawaiʻi: Melt concentrations, magma storage depths, and magma recycling

Kīlauea Volcano’s 2018 lower East Rift Zone (LERZ) eruption produced exceptionally high lava effusion rates and record-setting SO 2 emissions. The eruption involved a diverse range of magmas, including primitive basalts sourced from Kīlauea’s summit reservoirs. We analyzed LERZ matrix glasses, melt inclusions, and host minerals to identify melt volatile contents and magma storage depths. The LERZ glasses and melt inclusions span nearly the entire compositional range previously recognized at Kīlauea. Melt inclusions in Fo 86-89 olivine from the main eruptive vent (fissure 8) underwent 70–170 °C cooling during transport in LERZ carrier melts, causing extensive post-entrapment crystallization and sulfide precipitation. Many of these melt inclusions have low sulfur (400–900 ppm) even after correction for sulfide formation. CO 2 and H 2 O vapor saturation pressures indicate shallow melt inclusion trapping depths (1–5 km), consistent with formation within Kīlauea’s Halemaʻumaʻu and South Caldera reservoirs. Many of these inclusions also have degassed δ 34 S values (− 1.5 to − 0.5‰). Collectively, these results indicate that some primitive melts experienced near-surface degassing before being trapped into melt inclusions. We propose that decades-to-centuries of repeated lava lake activity and lava drain-back during eruptions (e.g., 1959 Kīlauea Iki) recycled substantial volumes of degassed magma into Kīlauea’s shallow reservoir system. Degassing and magma recycling from the 2008–2018 Halemaʻumaʻu lava lake likely reduced the volatile contents of LERZ fissure 8 magmas, resulting in lower fountain heights compared to many prior Kīlauea eruptions. The eruption’s extreme SO 2 emissions were due to high lava effusion rates rather than particularly volatile-rich melts.

Hawaii↗

USGS2020

Vapor-bubble growth in olivine-hosted melt inclusions

Melt inclusions record the depth of magmatic processes, magma degassing paths, and volatile budgets of magmas. Extracting this information is a major challenge. It requires determining melt volatile contents at the time of entrapment when working with melt inclusions that have suffered post-entrapment modifications. Several processes decrease internal melt inclusion pressure, resulting in nucleation and growth of a vapor bubble and, time permitting, diffusion of volatiles (especially CO2) into the vapor bubble. Methods exist that attempt to reconstruct the entrapped CO2 contents, but they are difficult to apply and yield inconsistent results. Here, we explore bubble growth, evaluate CO2 reconstruction approaches, and develop improved experimental and computational approaches. Piston-cylinder experiments were conducted on olivine-hosted melt inclusions from Seguam (Alaska, USA) and Fuego (Guatemala) volcanoes at the following conditions: 500-800 MPa, 1140-1200 °C for Seguam and 1110-1140 °C for Fuego, 4-8 wt% H2O in the KBr brine, and run durations of 10-120 minutes. Heated melt inclusions form well-defined S-CO2 trends that can be described by degassing models. CO2 contents are enriched by a factor of ~2.5, on average, relative to those of the glasses within unheated melt inclusions, whereas S contents of heated and unheated melt inclusion glasses overlap, indicating insignificant amounts of S partition into the vapor bubble. Low closure temperatures enable CO2 diffusion into vapor bubbles during quench upon eruption, while a higher closure temperature for S limits its loss to vapor bubbles. We evaluate the timescales of post-entrapment processes and use the results to develop a new computational model to restore entrapped CO2 contents: MIMiC (Melt Inclusion Modification Corrections). Heated melt inclusion data are used as a benchmark to evaluate of the results from MIMiC and other published methods of CO2 reconstruction. The methods perform variably well. Key advantages to our experimental rehomogenization technique are that it enables accurate measurements of CO2 contents and allows for large quantities of melt inclusions to be rehomogenized efficiently. Our new computational model produces more accurate results than other computational methods, has similar accuracy to the Raman method of CO2 reconstruction in cases where Raman can be applied (i.e., no C-bearing phases in bubble), and can be applied to the vast body of published melt inclusion data. To obtain the most robust data on bubble-bearing melt inclusions, we recommend taking both experimental- and MIMiC-based approaches.

American Mineralogist↗

USGS2019

Direct measurements of copper speciation in basaltic glasses: Understanding the relative roles of sulfur and oxygen in copper complexation in melts

Micro-analytical determination of copper (Cu) speciation in natural magmatic glasses, equilibrated below the nickel – nickel oxide (NNO) buffer, reveals that two copper species are commonly stabilized in such basaltic melts. X-ray absorption fine structure (XAFS) spectroscopic analysis of basaltic matrix glasses and melt inclusions (MI) from samples of mid-ocean ridge basalt (MORB), and from Nyamuragira, Etna and Kīlauea volcanoes show that both Cu(I)-sulfide and Cu(I)-oxide species are stabilized. The proportion of each species correlates with the measured sulfur (S) abundance of the glass. In glasses with S abundances greater than ~1000 ppm, Cu(I)-sulfide species are dominant, whereas in glasses with S abundances between 500 and 1000 ppm, both species are found to coexist. The Cu(I)-oxide species dominate at S concentrations below 500 ppm. In 1 atm S-free experimental glasses of basaltic composition that we analyzed, only Cu(I)-oxide species are detectable, regardless of the oxygen fugacity (fO2), even at relatively high fO2 values well above the NNO buffer. Our results demonstrate that XAFS techniques are highly sensitive in measuring Cu speciation in reduced (below NNO) basaltic glasses and that both oxide and sulfide complexes can be stabilized. The relative proportion of these two species is highly dependent on the concentration of S in the melt, and thus the Cu speciation in natural melts changes as S is lost from the melt by low pressure degassing.

Hawaii↗

USGS2016

Slab melting and magma formation beneath the southern Cascade arc

The processes that drive magma formation beneath the Cascade arc and other warm-slab subduction zones have been debated because young oceanic crust is predicted to largely dehydrate beneath the forearc during subduction. In addition, geochemical variability along strike in the Cascades has led to contrasting interpretations about the role of volatiles in magma generation. Here, we focus on the Lassen segment of the Cascade arc, where previous work has demonstrated across-arc geochemical variations related to subduction enrichment, and H-isotope data suggest that H 2 O in basaltic magmas is derived from the final breakdown of chlorite in the mantle portion of the slab. We use naturally glassy, olivine-hosted melt inclusions (MI) from the tephra deposits of eight primitive ( MgO>7 wt% ) basaltic cinder cones to quantify the pre-eruptive volatile contents of mantle-derived melts in this region. The melt inclusions have B concentrations and isotope ratios that are similar to mid-ocean ridge basalt (MORB), suggesting extensive dehydration of the downgoing plate prior to reaching sub-arc depths and little input of slab-derived B into the mantle wedge. However, correlations of volatile and trace element ratios (H 2 O/Ce, Cl/Nb, Sr/Nd) in the melt inclusions demonstrate that geochemical variability is the result of variable addition of a hydrous subduction component to the mantle wedge. Furthermore, correlations between subduction component tracers and radiogenic isotope ratios show that the subduction component has less radiogenic Sr and Pb than the Lassen sub-arc mantle, which can be explained by melting of subducted Gorda MORB beneath the arc. Agreement between pMELTS melting models and melt inclusion volatile, major, and trace element data suggests that hydrous slab melt addition to the mantle wedge can produce the range in primitive compositions erupted in the Lassen region. Our results provide further evidence that chlorite-derived fluids from the mantle portion of the slab ( ∼7–9 km below the slab top) cause flux melting of the subducted oceanic crust, producing hydrous slab melts that migrate into the overlying mantle, where they react with peridotite to induce further melting.

Cascades↗

USGS2016

Magma transport and olivine crystallization depths in Kīlauea’s East Rift Zone inferred from experimentally rehomogenized melt inclusions

Concentrations of H 2 O and CO 2 in olivine-hosted melt inclusions can be used to estimate crystallization depths for the olivine host. However, the original dissolved CO 2 concentration of melt inclusions at the time of trapping can be difficult to measure directly because in many cases substantial CO 2 is transferred to shrinkage bubbles that form during post-entrapment cooling and crystallization. To investigate this problem, we heated olivine from the 1959 Kīlauea Iki and 1960 Kapoho (Hawai‘i) eruptions in a 1-atm furnace to temperatures above the melt inclusion trapping temperature to redissolve the CO 2 in shrinkage bubbles. The measured CO 2 concentrations of the experimentally rehomogenized inclusions (⩽590 ppm for Kīlauea Iki [n=10]; ⩽880 ppm for Kapoho, with one inclusion at 1863 ppm [n=38]) overlap with values for naturally quenched inclusions from the same samples, but experimentally rehomogenized inclusions have higher within-sample median CO 2 values than naturally quenched inclusions, indicating at least partial dissolution of CO 2 from the vapor bubble during heating. Comparison of our data with predictions from modeling of vapor bubble formation and published Raman data on the density of CO 2 in the vapor bubbles suggests that 55-85% of the dissolved CO 2 in the melt inclusions at the time of trapping was lost to post-entrapment shrinkage bubbles. Our results combined with the Raman data demonstrate that olivine from the early part of the Kīlauea Iki eruption crystallized at <6 km depth, with the majority of olivine in the 1-3 km depth range. These depths are consistent with the interpretation that the Kīlauea Iki magma was supplied from Kīlauea’s summit magma reservoir (&sim;2-5 km depth). In contrast, olivine from Kapoho, which was the rift zone extension of the Kīlauea Iki eruption, crystallized over a much wider range of depths (&sim;1-16 km). The wider depth range requires magma transport during the Kapoho eruption from deep beneath the summit region and/or from deep beneath Kīlauea’s east rift zone. The deeply derived olivine crystals and their host magma mixed with stored, more evolved magma in the rift zone, and the mixture was later erupted at Kapoho.

Hawaii↗

USGS2015

Slab melting beneath the Cascades Arc driven by dehydration of altered oceanic peridotite

Water is returned to Earth’s interior at subduction zones. However, the processes and pathways by which water leaves the subducting plate and causes melting beneath volcanic arcs are complex; the source of the water—subducting sediment, altered oceanic crust, or hydrated mantle in the downgoing plate—is debated; and the role of slab temperature is unclear. Here we analyse the hydrogen-isotope and trace-element signature of melt inclusions in ash samples from the Cascade Arc, where young, hot lithosphere subducts. Comparing these data with published analyses, we find that fluids in the Cascade magmas are sourced from deeper parts of the subducting slab—hydrated mantle peridotite in the slab interior—compared with fluids in magmas from the Marianas Arc, where older, colder lithosphere subducts. We use geodynamic modelling to show that, in the hotter subduction zone, the upper crust of the subducting slab rapidly dehydrates at shallow depths. With continued subduction, fluids released from the deeper plate interior migrate into the dehydrated parts, causing those to melt. These melts in turn migrate into the overlying mantle wedge, where they trigger further melting. Our results provide a physical model to explain melting of the subducted plate and mass transfer from the slab to the mantle beneath arcs where relatively young oceanic lithosphere is subducted.

Cascade Arc↗

USGS2013

Magma-ice-sediment interactions and the origin of lava/hyaloclastite sequences in the Síða formation, South Iceland

Products of subglacial volcanism can illuminate reconstructions of paleo-environmental conditions on both local and regional scales. Competing interpretations of Pleistocene conditions in south Iceland have been proposed based on an extensive sequence of repeating lava-and-hyaloclastite deposits in the Síða district. We propose here a new eruptive model and refine the glacial environment during eruption based on field research and analytical data for the Síða district lava/hyaloclastite units. Field observations from this and previous studies reveal a repeating sequence of cogenetic lava and hyaloclastite deposits extending many kilometers from their presumed eruptive source. Glasses from lava selvages and unaltered hyaloclastites have very low H 2 O, S, and CO 2 concentrations, indicating significant degassing at or close to atmospheric pressure prior to quenching. We also present a scenario that demonstrates virtual co-emplacement of the two eruptive products. Our data and model results suggest repeated eruptions under thin ice or partially subaerial conditions, rather than eruption under a thick ice sheet or subglacial conditions as previously proposed.

Bulletin of Volcanology↗

USGS2005

Storage and interaction of compositionally heterogeneous magmas from the 1986 eruption of Augustine Volcano, Alaska

Compositional heterogeneity (56–64 wt% SiO 2 whole-rock) in samples of tephra and lava from the 1986 eruption of Augustine Volcano, Alaska, raises questions about the physical nature of magma storage and interaction beneath this young and frequently active volcano. To determine conditions of magma storage and evolutionary histories of compositionally distinct magmas, we investigate physical and chemical characteristics of andesitic and dacitic magmas feeding the 1986 eruption. We calculate equilibrium temperatures and oxygen fugacities from Fe-Ti oxide compositions and find a continuous range in temperature from 877 to 947°C and high oxygen fugacities (ΔNNO=1–2) for all magmas. Melt inclusions in pyroxene phenocrysts analyzed by Fourier-transform infrared spectroscopy and electron probe microanalysis are dacitic to rhyolitic and have water contents ranging from <1 to ∼7 wt%. Matrix glass compositions are rhyolitic and remarkably similar (∼75.9–76.6 wt% SiO 2 ) in all samples. All samples have ∼25% phenocrysts, but lower-silica samples have much higher microlite contents than higher-silica samples. Continuous ranges in temperature and whole-rock composition, as well as linear trends in Harker diagrams and disequilibrium mineral textures, indicate that the 1986 magmas are the product of mixing between dacitic magma and a hotter, more mafic magma. The dacitic endmember is probably residual magma from the previous (1976) eruption of Augustine, and we interpret the mafic endmember to have been intruded from depth. Mixing appears to have continued as magmas ascended towards the vent. We suggest that the physical structure of the magma storage system beneath Augustine contributed to the sustained compositional heterogeneity of this eruption, which is best explained by magma storage and interaction in a vertically extensive system of interconnected dikes rather than a single coherent magma chamber and/or conduit. The typically short repose period (∼10 years) between Augustine's recent eruptive pulses may also inhibit homogenization, as short repose periods and chemically heterogeneous magmas are observed at several volcanoes in the Cook Inlet region of Alaska.

Alaska↗