Geology topics

Thomas Shea

Publications and source records attributed to Thomas Shea.

9 recordsLinked to original sources

USGS2026

Timescales of cumulate mobilization and mixing for the 1868 A.D. eruption of Mauna Loa, Island of Hawai‘i

The deadly 1868 A.D. eruption of Mauna Loa’s lower Southwest Rift Zone (Island of Hawai‘i) included a M7.9 earthquake and associated tsunami and landslides, demonstrating the severe hazards posed by Earth’s largest active subaerial volcano. To better understand the relationship between intense seismic activity, dike emplacement, magma storage, transport histories, and mobilization of olivine cumulates at Mauna Loa, we examine compositional zoning of olivine in the 1868 lava flows. Samples range from basalt (< 10% olivine) to picrite (30–40% olivine). The olivine cargo is heterogeneous (Fo 78.2–89.2 ; forsterite = [Mg/(Mg + Fe) × 100]) but dominated by ~ Fo 89 cores that lie above the Fe-Mg equilibrium field of host glasses. Crystal rims < Fo 80 are due to post-eruptive modification in slow cooling lava flows. Minor element compositions fall within the range of other Mauna Loa olivine erupted in the past 200 years. Olivine crystals exhibit both normal and complex Fo zoning patterns that yield timescales of diffusive re-equilibration that range from 3 to 258 days, with 72% of crystals recording 71 days or less. These timescales correspond to magmatic priming of the summit reservoir system ~ 2 months prior to the eruption and the M7.9 earthquake likely facilitated the transport of the crystal-rich summit-derived magmas downrift shortly prior to eruption. If the recently proposed faster Fe-Mg diffusion coefficient is used, timescales instead range from < 1 day to 25 days, with most recording 1 week or less. In this scenario, most of the olivine zoning would have to have been generated after the M7.9 earthquake perturbed the system.

Hawaii↗

USGS2025

Rapid emplacement of the Keaiwa Lava Flow of 1823 from the Great Crack in the Southwest Rift Zone of Kilauea volcano

The Keaīwa Lava Flow of 1823 in the Southwest Rift Zone of Kīlauea volcano is unusual for its expansive pāhoehoe sheet flow morphology and lack of constructive vent topography, despite having a similar tholeiitic basalt composition to other lavas erupted from Kīlauea. This lava flow issued from a ∼10-km-long continuous fissure now known as the Great Crack, and has an unusually thin sheet flow morphology with margin thicknesses of ∼15–110 cm (average of 42 cm). Based on field observations of the lava flow at its fissure vent (e.g., drain-back features), we propose that the Great Crack formed, or at least significantly widened, just prior to and syn-eruptively with this 1823 eruption. The absence of pyroclastic cones or spatter ramparts indicates that the eruption consisted of a rapid outpouring of relatively degassed lava as the fissure unzipped. The rapidly moving lava flow overtopped pre-existing tumuli and scoria cones (e.g., Lava Plastered Cones) up to ∼10 m tall. Glass and whole-rock chemistry yield homogeneous compositions for the lavas erupted from the Great Crack, with glass compositions of 6.40 ± 0.10 wt% MgO and whole-rock compositions of 7.39 ± 0.07 wt% MgO. Lava pads erupted from a short western fissure system are richer in mafic minerals (e.g., olivine and clinopyroxene), and show slightly more MgO-rich whole-rock compositions (7.79 ± 0.05 wt%). MgO-in-glass thermometry on juvenile spatter yield eruption temperatures of 1153 ± 13°C that are typical of Kīlauea lavas. Thus, the extensive sheet-like lava flow morphology is not a direct consequence of unusual magmatic or rheological conditions (i.e., low viscosity). Instead, the flow morphology is associated with high effusion rates caused by sudden drainage of uprift magma as it erupted from the Great Crack. Lava flow modeling on a 2-m-resolution digital elevation model indicates that a minimum bulk effusion rate of ∼5800 m 3 /s (∼3500 m 3 /s dense rock equivalent) and a minimum flow velocity of ∼11 m/s are required for the lava flow to overcome the topography of the Lava Plastered Cones. This effusion rate is among the highest inferred for eruptions in Hawaiʻi and around the world. This study highlights a less frequent eruption style at Hawaiian volcanoes characterized by a sudden outpouring of lava from an unusual fissure system. Local eyewitness accounts indicate that the 1823 eruption was preceded by seismicity. Given the complex magmatic-volcanic-tectonic relations across Kīlauea, we speculate that the south flank could have slipped over one or more events that ultimately triggered unzipping of the Great Crack and passive release of briefly stored uprift magma. An eruption similar to 1823 at Kīlauea or Mauna Loa, with an eruptive timeframe that could be as short as an hour, with high effusion rates and rapid flow front velocities, would not easily allow for a timely response.

Hawaii↗

USGS2023

The presence of silicate melt may enhance rates of cation diffusion in olivine

Olivine is commonly used as a ‘crystal clock’ to extract timescales relevant to pre-eruptive perturbations within mafic magmatic systems. Diffusion chronometry applications require accurate calibrations for the rates at which Fe-Mg or other commonly measured elements like Ni, Mn, and Ca diffuse through the crystal lattice. In the past, these rates have been mainly characterized using solid-solid diffusion couple experiments involving olivine single crystals, thin films, or powder sources. Despite the presence of melt surrounding olivine in natural magmatic systems, very few experiments involving magma have been performed, largely because controlling interface reactions is difficult. For this study, we carried out olivine-melt diffusion experiments as a test of the diffusion chronometry method, and to determine whether the presence of melt influences the calculated timescales. To approximate a natural system, we incorporated small natural Kīlauea and San Carlos olivine seeds within a natural Kīlauea basalt and tracked diffusive re-equilibration through time. To better control interface reactions, after some equilibration period at an initial superliquidus temperature of 1290 °C, the runs were rapidly cooled to form a rim and left to dwell at various final temperatures (1200, 1220, 1240, 1255 °C) for 6–84 h. Concentration gradients for Fe-Mg, Mn, Ni, Ca were measured, and the step-wise nature of the core-rim transition was ascertained using slow diffusing elements like P or Al. When these gradients are modeled using published diffusivities, the timescales retrieved are typically 10 times longer than the actual experiment durations. Thus, measured diffusivities are an order of magnitude faster than those previously obtained in olivine-solid source experiments, but they are in excellent agreement with the only two other melt-olivine datasets. We explore reasons for why melt-bearing olivine diffusion experiments tend to yield faster rates. The possible effects of (1) growth during diffusion, (2) diffusion during any initial dissolution step, and (3) extended tube or planar defects at the interface on calculated diffusivities are all considered but found to be inconsequential. Instead, we argue that additional point defects (vacancies) are likely created at the interface by higher concentrations in elements like Al or H in the basalt melt compared to other solid couple diffusant sources. Future applications of diffusion chronometry in olivine may require a complete re-evaluation of published diffusivities using melt-bearing experimental configurations.

Earth and Planetary Science Letters↗

USGS2022

Trace elements in olivine fingerprint the source of 2018 magmas and shed light on explosive-effusive eruption cycles at Kīlauea Volcano

Understanding magma genesis and the evolution of intensive parameters (temperature, pressure, composition, degree of melting) in the mantle source of highly active volcanic systems is crucial for interpreting magma supply changes over time and recognizing cyclic behavior to anticipate future volcanic behavior. Major and trace elements in olivine are commonly used to study variations in mantle lithologies and melting conditions (e.g., temperature, pressure, oxygen fugacity) affecting the mantle over time. Here, we track the temporal evolution of primary melts through the most recent cycle of explosive and effusive eruptions at Kīlauea (Hawai‘i), which spans the last ∼500 years. We report major and trace elements in olivine from the last explosive period (∼1500 – early 1820’s Keanakāko‘i Tephra) and the most recent decade of the current effusive period (2018 LERZ, 2015–2018 Pu‘u‘ō‘ō, 2008–2018 lava lake and 2020 eruption in Halema‘uma‘u). Scandium concentrations in olivine allow characterizing changes in mantle source between 1500 and 2018, and suggest that the recent (2015–2018) magma feeding the Pu‘u‘ō‘ō cone did not significantly interact with the magma that erupted in the LERZ in 2018. The evolution of olivine and melt compositions over the past 500 years is not easily reconcilable with variations in mantle potential temperature, pressure of mantle melt pooling and storage, or oxygen fugacity. Instead, Sc, Mn, and Co concentrations and Ni/Mg ratio in high forsterite (Fo >87) olivine advocate for an increase in the proportion of clinopyroxene in the mantle source associated with a slightly higher degree of partial melting from 1500 to 2018. Changes in primitive melt compositions and degrees of mantle melting may well modulate magma supply to the crust and formation-replenishment of steady or ephemeral summit reservoirs, and thereby control transitions between explosive and effusive periods at Kīlauea. Analyzing trace elements in olivine at Kīlauea and elsewhere could therefore provide important clues on subtle changes occurring at the mantle level that might herald changes in volcanic behavior.

Hawaii↗

USGS2021

Evaluating lava flow propagation models with a case study from the 2018 eruption of Kīlauea Volcano, Hawai'i

The 2018 lower East Rift Zone (LERZ) eruption of Kīlauea, Hawai’i, provides an excellent natural laboratory with which to test models of lava flow propagation. During early stages of eruption crises, the most useful lava flow propagation equations utilize readily determined parameters and require fewer a priori assumptions about future behavior of the flow. Here, we leverage the numerous observations of lava flows collected over the duration of the eruption crisis at Kīlauea in 2018 to test simple lava flow propagation models. These models track the one-dimensional propagation of the flows according to three main rheological restraining forces: bulk viscosity, yield strength, and growth of a surface crust. We calculate the predicted changes in length through time of three flows that vary in bulk composition, crystal content, and total flow length. Cooler flows that are more crystal-rich tend to be more dominated by crust growth, though early stages of propagation can be controlled by bulk viscosity. We find that variations in effusion rate significantly impact flows that are short-lived; flows that are produced during steady-state effusion are readily approximated by average values for the entire flow. Thus, accurate knowledge of variations in effusion rate are critical to accurate lava flow propagation forecasting.

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

Capturing, preserving and digitizing legacy seismic data from the Montserrat Volcano Observatory analog seismic network, July 1995 – December 2004

An eruption of the Soufrière Hills Volcano (SHV) on the eastern Caribbean island of Montserrat began on 18 July 1995 and continued until February 2010. Within nine days of the eruption onset, an existing four‐station analog seismic network (ASN) was expanded to 10 sites. Telemetered data from this network were recorded, processed, and archived locally using a system developed by scientists from the U.S. Geological Survey (USGS) Volcano Disaster Assistance Program (VDAP). In October 1996, a digital seismic network (DSN) was deployed with the ability to capture larger amplitude signals across a broader frequency range. These two networks operated in parallel until December 2004, with separate telemetry and acquisition systems (analysis systems were merged in March 2001). Although the DSN provided better quality data for research, the ASN featured superior real‐time monitoring tools and captured valuable data including the only seismic data from the first 15 months of the eruption. These successes of the ASN have been rather overlooked. This article documents the evolution of the ASN, the VDAP system, the original data captured, and the recovery and conversion of more than 230,000 seismic events from legacy SUDS, Hypo71, and Seislog formats into Seisan database with waveform data in miniSEED format. No digital catalog existed for these events, but students at the University of South Florida have classified two-thirds of the 40,000 events that were captured between July 1995 and October 1996. Locations and magnitudes were recovered for ~10,000 of these events. Real-time seismic amplitude measurement, seismic spectral amplitude measurement, and tiltmeter data were also captured. The result is that the ASN seismic dataset is now more discoverable, accessible, and reusable, in accordance with FAIR data principles. These efforts could catalyze new research on the 1995–2010 SHV eruption. Furthermore, many observatories have data in these same legacy data formats and might benefit from procedures and codes documented here.

Soufrière Hills Volcano↗

USGS2020

The historic events at Kilauea Volcano in 2018: Summit collapse, rift zone eruption, and Mw 6.9 earthquake: Preface to the special issue

Kīlauea Volcano, on the Island of Hawaiʻi, has had a prominent role in the science of volcanology, and a long history of generating new insights into how volcanoes operate (Tilling et al. 2014; Garcia 2015). Native Hawaiians shared ideas on the behavior of the volcano with early Western visitors to Kīlauea, addressing the basic geometry of magma supply and transport (Ellis 1825; Bishop 1827). The recognition that magma originated at the summit and was transferred at shallow levels to the flanks implied that these ideas were rooted in centuries of observation preceding Western contact. The lava lake activity at Kīlauea’s summit in the 1800s and early 1900s fascinated early geologists, such as James Dana (1890), who published one of the first inquiries into the fundamental processes of Hawaiian volcanoes. The sustained activity led to the 1912 founding of the Hawaiian Volcano Observatory, one of the world’s first volcano observatories, by Thomas Jaggar (Tilling et al. 2014). Kīlauea’s activity in the 20th century contributed to the development of many modern volcano monitoring techniques (Tilling et al. 2014), which helped refine conceptual models of how volcanoes behave (Eaton and Murata, 1960).

Hawaii↗

USGS2017

Timescales of mixing and storage for Keanakāko‘i Tephra magmas (1500-1823 C.E.), Kīlauea Volcano, Hawai‘i

The last 2500 years of activity at Kīlauea Volcano (Hawai‘i) have been characterized by centuries-long periods dominated by either effusive or explosive eruptions. The most recent period of explosive activity produced the Keanakāko‘i Tephra (KT; ca. 1500–1820 C.E.) and occurred after the collapse of the summit caldera (1470–1510 C.E.). Previous studies suggest that KT magmas may have ascended rapidly to the surface, bypassing storage in crustal reservoirs. The storage conditions and rapid ascent hypothesis are tested here using chemical zoning in olivine crystals and thermodynamic modeling. Forsterite contents (Fo; [Mg/(Mg + Fe) × 100]) of olivine core and rim populations are used to identify melt components in Kīlauea’s prehistoric (i.e., pre-1823) plumbing system. Primitive (≥Fo 88 ) cores occur throughout the 300+ years of the KT period; they originated from mantle-derived magmas that were first mixed and stored in a deep crustal reservoir. Bimodal olivine populations (≥Fo 88 and Fo 83–84 ) record repeated mixing of primitive magmas and more differentiated reservoir components shallower in the system, producing a hybrid composition (Fo 85–87 ). Phase equilibria modeling using MELTS shows that liquidus olivine is not stable at depths >17 km. Thus, calculated timescales likely record mixing and storage within the crust. Modeling of Fe–Mg and Ni zoning patterns (normal, reverse, complex) reveal that KT magmas were mixed and stored for a few weeks to several years before eruption, illustrating a more complex storage history than direct and rapid ascent from the mantle as previously inferred for KT magmas. Complexly zoned crystals also have smoothed compositional reversals in the outer 5–20 µm rims that are out of Fe–Mg equilibrium with surrounding glasses. Diffusion models suggest that these rims formed within a few hours to a few days, indicating that at least one additional, late-stage mixing event may have occurred shortly prior to eruption. Our study illustrates that the lifetimes of KT magmas are more complex than previously proposed, and that most KT magmas did not rise rapidly from the mantle without modification during shallow crustal storage.

Hawaii↗