Volcanic seismicity beneath Chuginadak Island, Alaska (Cleveland and Tana volcanoes): Implications for magma dynamics and eruption forecasting

John A. Power; Diana Roman; John J. Lyons; Matthew M. Haney; Daniel J. Rasmussen; Terry A. Plank; K. P. Nicolaysen; Pavel Izbekov; C. Werner; A Kaufman

doi:10.1016/j.jvolgeores.2021.107182

USGS · 70249479

Volcanic seismicity beneath Chuginadak Island, Alaska (Cleveland and Tana volcanoes): Implications for magma dynamics and eruption forecasting

John A. Power·Diana Roman·John J. Lyons·Matthew M. Haney·Daniel J. Rasmussen·Terry A. Plank·K. P. Nicolaysen·Pavel Izbekov·C. Werner·A Kaufman

View original source ↗Export metadata ↓

Abstract

Cleveland and Tana are remote volcanoes located in the central Aleutian volcanic arc on the eastern end of the Islands of Four Mountains (IFM). The persistently active Mount Cleveland volcano, on the western side of Chuginadak Island, is surrounded by several closely spaced Quaternary volcanic centers including Carlisle, Herbert, Kagamil, Tana, and Uliaga, and numerous small satellite vents on Chiginadak between Cleveland and Tana. The Alaska Volcano Observatory (AVO) installed two permanent broadband seismometers on Chuginadak Island in 2014, and we operated a temporary broadband network focused on the western side of the island in 2015–2016. Collectively, these stations provided the first seismic observations of this frequently active volcano and the surrounding Holocene-aged volcanic vents. During the study period (July 2014–January 2019), eruptive activity at Cleveland was characterized by small explosions separated by periods of lava effusion that formed small domes in the volcano's summit crater. We characterize seismicity beneath Chuginadak Island through automated analysis of event waveform frequency content, development of a one-dimensional P-wave velocity model, calculation of earthquake hypocenters, magnitudes, focal mechanisms, and identification of earthquake families. This analysis reveals the full range of seismic event types expected in a highly active volcanic environment and includes Volcano-Tectonic (VT) earthquakes, Long-Period (LP) events, and explosion signals. LP events appear to cluster at shallow depth beneath the active crater of Mount Cleveland and almost all of the explosions occur without identifiable short-term (hours to days) seismic precursors. VT earthquakes beneath Mount Cleveland occur at depths of 2 to 8 km below sea level (BSL) and range in magnitude from −0.2 to 1.8. VT focal mechanisms have horizontal P-axes that align with the regional axis of maximum stress. These observations, and a relatively slow one-dimensional seismic velocity model, are consistent with a shallow body of magma that is fed through a deeper conduit system. The time-history of VT earthquakes and shallow LP events suggest their occurrence may track the transfer of magma and fluids from the mid-crust to the shallow portions of the conduit system and may provide a means to anticipate future explosions and periods of dome growth. VT hypocenters also extend ~7 km northeast of Cleveland's summit at depths of 5 to 10 km BSL, under a group of Holocene-aged vents between Mount Cleveland and Tana. These earthquakes have vertically-oriented P-axes and a greater percentage occur in families. These observations, combined with observations of vent orientation and morphology and gas flux, suggest the area between Cleveland and Tana represents a zone of complicated volcano-tectonic interaction, similar to calderas elsewhere in the Aleutian arc. The presence of a larger volcanic system in the eastern IFM could influence magmatism and account for the multiple closely spaced volcanic centers in this region.

Explore related subjects

Journal of Volcanology and Geothermal Research Article

Source-reported bounding extent: 52.76373370379605° to 52.904805932105404° latitude; -170.01669419707966° to -169.65098413866693° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

John A. Power, Diana Roman, John J. Lyons, Matthew M. Haney, Daniel J. Rasmussen, Terry A. Plank, K. P. Nicolaysen, Pavel Izbekov, C. Werner, A Kaufman. 2021. Volcanic seismicity beneath Chuginadak Island, Alaska (Cleveland and Tana volcanoes): Implications for magma dynamics and eruption forecasting. https://doi.org/10.1016/j.jvolgeores.2021.107182

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Tree-ring isotope and satellite records of vegetation changes reveal transition from cold to hot degassing in the Yellowstone Plateau volcanic field

Magmatic degassing often provides the earliest clues to changes in magmatic systems, yet detection and characterization of the escaping gases is challenging. Recent work suggests that satellite-based normalized difference vegetation index (NDVI) calculations may provide an indirect method for detecting changes in degassing. Here, we use tree ring δ 13 ⁢C , δ 18 ⁢O , and radiocarbon ( 14 C) data to assess NDVI as a proxy for diffuse magmatic degassing and to better understand the development of the Tern Lake thermal area (TLTA) in the Yellowstone Plateau volcanic field, USA. We measured δ 13 ⁢C , δ 18⁢ O , and 14 C in wood from four TLTA trees and two control trees. The 14 C data indicate that TLTA trees were exposed to significant magmatic CO 2 degassing during the 1990s, corroborating satellite-based NDVI changes observed at a time when no other indications of degassing existed. Magmatic degassing also increased water availability for trees closest to the degassing center, as indicated by decreasing δ 13 ⁢C and δ 18⁢ O values relative to background throughout the 1980s and 1990s. Increased H 2 S degassing in the early 2000s impaired stomatal functioning, increasing δ 13 ⁢C and δ 18⁢ O values. These isotope values were strongly correlated with NDVI, demonstrating that ground- and satellite-based approaches yield similar conclusions. Together, the satellite and isotope datasets suggest that TLTA degassing was characterized by cold CO 2 emissions with minimal H 2 S until the early 2000s, when it transitioned to a high-temperature system with significant CO 2 and H 2 S emissions. These results add to a growing body of evidence that vegetation responses can be an effective proxy for magmatic degassing.

Wyoming↗

USGS2026

Dynamics and evolution of the Kı̄lauea lower East Rift Zone 2018 fissure 8 lava flow and implications for multiphase magma properties

The 2018 Kı̄lauea lower East Rift Zone (KLERZ) eruption was one of the most voluminous eruptions on the Island of Hawai’i in the past 200 years, leading to major disruption and destroying over 700 homes and structures. The majority of the erupted magma was emitted as a lava flow from Ahu’ailā’au (fissure 8), which was active from late May to early August. To better understand the evolution of long-lived channelized lava flows, we examined the evolution of velocity, texture, and inferred rheology of the fissure 8 lava in space and time. We quantified lava flow surface velocities using particle image velocimetry in more than 200 aerial videos that span the lava flow duration and length. Velocity measurements were analyzed together with vesicularity and crystallinity measurements from 9 co-located post-eruptive field samples to understand the textural evolution of this flow and its impact on lava rheology and flow velocity. The fissure 8 flow was highly vesicular, with 79%–88% vesicularity at the vent, decreasing to 16%–26% vesicularity 12.5 km from the vent. The volume fraction occupied by crystals 50 in size increased from 6% at the vent to about 18% at 12.5 km downstream. We find that the effective flow viscosity increased at a quadratic rate with distance. Using experimentally determined liquid viscosity and applying established models to account for the effect of crystals and bubbles, we attribute this increase primarily to textural evolution driven initially by the near-vent loss of deformable bubbles and later by cooling and crystal growth. We demonstrate the importance of accounting for the evolution of vesicularity and the role of vesicles by showing that utilizing this capability in the open-source thermo-rheological lava flow propagation model PyFLOWGO allows for more accurate predictions of the observed flow velocities. Our modeling results suggest that small bubbles behaving rigidly are required to simulate the observed flow length, speed, and viscosities. Flow velocities of the fissure 8 lava also varied with time, driven by near-daily collapse events of the summit caldera. Temporal velocity changes were characterized by a period of steep acceleration, with the volumetric flux peaking around 4 h after a caldera collapse, followed by a period of gradual deceleration lasting up to 40 h or until the next collapse event. We use this temporal behavior to estimate the compressibility of the magma inside the plumbing system between the summit reservoir and the lower East Rift Zone. Overall, quantifying the spatial and temporal evolution of the KLERZ eruption provides information about magma and lava properties that can inform predictive modeling and hazard assessment during an eruption.

Hawaii↗

USGS2026

Small-volume tephra deposits of the May 1924 explosions from Halemaʻumaʻu, Kīlauea volcano, and their origin

More than 50 explosive eruptions occurred from Halemaʻumaʻu at Kīlauea volcano over 17 days from May 11 to 27, 1924. Ballistics weighing as much as 14,000 kg were ejected and most landed within 2 km of the vent. Fine ash made up a major component of the tephra and was dispersed tens of kilometers downwind. Draining of the Halemaʻumaʻu lava lake occurred in late February 1924, with the crater floor eventually subsiding by a further ∼70 m (to ∼180 m below the crater rim) by the time the first explosions took place during the night of May 10–11. The largest explosions occurred on May 17–18 and smaller explosions continued until May 27, at which point Halemaʻumaʻu had more than doubled in width and depth. The explosions generated plumes reaching up to ∼10 km high with ballistics ejected up to 2 km from the crater. Almost 100 years later, we investigate and characterize the preserved tephra deposits within ∼3 km of the 1924 crater rim. Grain size and shape analyses were performed on 202 samples collected from 34 tephra profiles using dynamic image analysis, with a subset of layers from nine tephra profiles used for componentry (200 grains per layer in the 0.5–1 mm size fraction). Additionally, we characterize the average diameters (using the five largest clasts) at 216 locations and measure the average diameters of 2291 ballistics (largest per ∼100 m 2 area). Physical descriptions from fieldwork and grain size distributions were used to subdivide the tephra layers into five lithofacies: coarse homogeneous, fine homogenous, red ash, accretionary lapilli-bearing, and finely laminated. Grain size versus shape data show a range of values that demonstrate most grains are dense, smooth, and equant, in alignment with lithic clasts dominating the tephra componentry. The fine grained and accretionary lapilli-bearing nature of some of these lithofacies confirms that water influenced the style of the explosions. However, we also note juvenile clasts within many of the tephra layers, indicating that many of the layers were formed during phreatomagmatic explosions (sensu stricto), despite the eruptive mechanism being dominantly phreatic. Juvenile clasts are more abundant higher in the tephra profiles, suggesting that juvenile magma was more involved later in the explosive sequence. Thermal and hydrologic modeling indicate that groundwater inflow into a short-lived, small-diameter volcanic conduit (10-m to 120-m-diameter used for modeling) during the 78–85 days preceding the first explosion provides a physically plausible mechanism for this eruptive sequence.

Hawaii↗