USGS ScienceSearch

USGS · 70175554

Volcano-tectonic earthquakes: A new tool for estimating intrusive volumes and forecasting eruptions

Abstract

We present data on 136 high-frequency earthquakes and swarms, termed volcano-tectonic (VT) seismicity, which preceded 111 eruptions at 83 volcanoes, plus data on VT swarms that preceded intrusions at 21 other volcanoes. We find that VT seismicity is usually the earliest reported seismic precursor for eruptions at volcanoes that have been dormant for decades or more, and precedes eruptions of all magma types from basaltic to rhyolitic and all explosivities from VEI 0 to ultraplinian VEI 6 at such previously long-dormant volcanoes. Because large eruptions occur most commonly during resumption of activity at long-dormant volcanoes, VT seismicity is an important precursor for the Earth's most dangerous eruptions. VT seismicity precedes all explosive eruptions of VEI ≥ 5 and most if not all VEI 4 eruptions in our data set. Surprisingly we find that the VT seismicity originates at distal locations on tectonic fault structures at distances of one or two to tens of kilometers laterally from the site of the eventual eruption, and rarely if ever starts beneath the eruption site itself. The distal VT swarms generally occur at depths almost equal to the horizontal distance of the swarm from the summit out to about 15 km distance, beyond which hypocenter depths level out. We summarize several important characteristics of this distal VT seismicity including: swarm-like nature, onset days to years prior to the beginning of magmatic eruptions, peaking of activity at the time of the initial eruption whether phreatic or magmatic, and large non-double couple component to focal mechanisms. Most importantly we show that the intruded magma volume can be simply estimated from the cumulative seismic moment of the VT seismicity from: Log10 V = 0.77 Log ΣMoment − 5.32, with volume, V, in cubic meters and seismic moment in Newton meters. Because the cumulative seismic moment can be approximated from the size of just the few largest events, and is quite insensitive to precise locations, the intruded magma volume can be quickly and easily estimated with few short-period seismic stations. Notable cases in which distal VT events preceded eruptions at long-dormant volcanoes include: Nevado del Ruiz (1984–1985), Pinatubo (1991), Unzen (1989–1995), Soufriere Hills (1995), Shishaldin (1989–1999), Tacana' (1985–1986), Pacaya (1980–1984), Rabaul (1994), and Cotopaxi (2001). Additional cases are recognized at frequently active volcanoes including Popocateptl (2001–2003) and Mauna Loa (1984). We present four case studies (Pinatubo, Soufriere Hills, Unzen, and Tacana') in which we demonstrate the above mentioned VT characteristics prior to eruptions. Using regional data recorded by NEIC, we recognized in near-real time that a huge distal VT swarm was occurring, deduced that a proportionately huge magmatic intrusion was taking place beneath the long dormant Sulu Range, New Britain Island, Papua New Guinea, that it was likely to lead to eruptive activity, and warned Rabaul Volcano Observatory days before a phreatic eruption occurred. This confirms the value of this technique for eruption forecasting. We also present a counter-example where we deduced that a VT swarm at Volcan Cosiguina, Nicaragua, indicated a small intrusion, insufficient to reach the surface and erupt. Finally, we discuss limitations of the method and propose a mechanism by which this distal VT seismicity is triggered by magmatic intrusion.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Randall A. White, Wendy McCausland. 2016. Volcano-tectonic earthquakes: A new tool for estimating intrusive volumes and forecasting eruptions. https://doi.org/10.1016/j.jvolgeores.2015.10.020

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

KEEP EXPLORING

Related USGS reports

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

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

Seismic characteristics of the April 2024 eruption of Ruang Volcano, North Sulawesi, Indonesia

The two most recent eruptions of Ruang volcano, in 2002 and 2024, have been large (VEI 4) with very short precursory periods, 3 and 6 days, respectively. Background seismicity at Ruang volcano is generally low with 0–2 volcano-tectonic (VT) earthquakes per day. The 2024 eruptions were preceded by a notable increase in VT earthquakes from 11 to 15 April with a rapid escalation from 16 to 17 April in the form of increased VT swarms, and the appearance of low frequency earthquakes and short bursts of VT drumbeats shortly before the eruption. Similar VT seismic unrest was recorded in March 2015, May 2016, and April 2022, although none of these episodes culminated in an eruption. A comparison of seismic energy release between April 2022 and April 2024 shows distinct differences. In 2024, cumulative seismic energy reached its peak within 5 days, just before the onset of eruption on 17 April. In contrast, during the 2022-episode, energy release followed a more gradual increase over 8 days peaking on 16 April, then declining gradually without leading to eruption. Notably, low frequency (LF) and drumbeat earthquakes were absent in non-eruptive unrest episodes. The transition to proximal LFs in 2024 likely signals pressure fluctuations, intensified magma degassing, and conduit formation, acting as an important precursor to eruption. Additionally, the occurrence of short bursts of drumbeat seismicity likely reflects the onset of the shallow magma ascent and was an important short-term precursor to the eruption. Finally, we examine the potential influence of regional tectonic earthquakes on the escalation of precursory unrest in 2024.

North Sulawesi, Ruang volcano