USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Volcanology and Geothermal Research”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17Linked to original sources
USGS2022

Revisiting the depth distribution of seismicity before and after the 2004–2008 eruption of Mount St. Helens

Hypocenter estimation at active volcanoes improves our understanding of their magmatic systems and indicates changing conditions at depth for continuously monitored volcanoes. The most active volcano in the Cascades Range, Mount St. Helens, has a multi-decadal earthquake catalog and it shows an apparent change in the depth distribution of seismicity before and after the 2004–2008 dome-building eruption and unrest sequence. We use two new resources to evaluate the accuracy of hypocenters and consequently the change in depth distribution of seismicity before and after the 2004–2008 eruption. First, we deployed a dense array of 136 three-component nodal seismographs for one month in 2017, including sub-arrays on the newly extruded dome and crater floor. Second, for events recorded during this month, we located their hypocenters using a three-dimensional (3D) wavefront-tracking location solver and a recently developed tomography model derived from active and passive source data. The relocated hypocenters are generally shallower and more concentrated beneath the crater compared to their catalog locations. The mean hypocenter movement from catalog locations is 2.86 km, with an averaged depth shift of 2.53 km upward. Comparison between 2017 hypocenters located using all of the available phase picks and those located using only the catalog picks from the permanent network suggests the improved hypocenters mostly resulted from the location solver and the 3D velocity model, with smaller changes due to the dense three-component array. Applying the same hypocenter estimation method to phase picks for all events from 1997 to 2004 and 2008–2021, we found a concentration of seismicity between sea level and ~ 1 km above it before and after the 2004–2008 eruption. The new results suggest that the seismogenic structure in the shallow magmatic system quickly re-equilibrated to its earlier state after the dome-building eruption.

Washington↗
USGS2022

Development, structure, and behavior of a perched lava channel at Kīlauea Volcano, Hawaiʻi, during 2007

Channelized lava flows are commonly produced during the early stages of basaltic eruptions. These channels usually maintain their morphology until the eruption ends or discharge is diverted. In some instances, narrower channels can roof over, developing into lava tubes. We report here on a channelized flow erupted at Kīlauea volcano in 2007 that evolved into a “perched lava channel” composed of a string of interconnected, elongate lava pools, forming a lava channel/lava pool hybrid. The lava channel, which had a time-averaged discharge rate of ~3–9 m 3 /s, initially fed a series of flow branches that exhibited cooling-limited and volume-limited controls on flow length, sometimes with each process controlling a different morphological aspect of a single flow branch. The perched lava channel grew vertically primarily by overplating of the channel levees from frequent overflows, forming a compound flow field. This vertical growth only occurred when the distal end of the channel was blocked. When levee failure at the distal end of the channel caused the lava level in the channel to drop below the levee rim, no vertical growth occurred. Seeps of spiny lava and slabby pāhoehoe were common, erupting from uplift scarps on the channel levees, apparently fed by sills from denser, relatively crystal-rich material filling the bottom of the channel. We infer that lava in the channel was stratified in vesicularity and velocity, with foamy, vesicular, faster-moving lava at the top of the lava stream and denser, relatively outgassed, slower-moving lava filling the bottom of the channel. The channel levees were unstable, failing on several occasions, perhaps triggered by the levee seeps. The appearance of seeps, therefore, is one way of assessing the collapse potential of similar perched lava structures.

Hawaii↗
USGS2022

Growth of complex volcanic ash aggregates in the Tierra Blanca Joven eruption of Ilopango Caldera, El Salvador

Aggregation processes control both the residence time and dispersal of volcanic ash during eruptions yet remain incompletely understood. The products of aggregation vary from simple ash clusters to large, complexly layered accretionary lapilli. Here we detail the micro-stratigraphy of a single population of accretionary lapilli that grew during the ∼431 CE Tierra Blanca Joven eruption from Ilopango Caldera, El Salvador. The accretionary lapilli were sampled 10 km from the caldera source within a sequence of ash-rich pyroclastic density current deposits and intercalated fall material, known as unit D, which is traceable >40 km from Ilopango. Scanning electron microscopy and image analysis reveal common facies that form distinct layers within the accretionary lapilli. Each facies is distinguished by quantitative and qualitative variations in particle size distribution , porosity, and particle fabric. We infer that these textures resulted from aggregation conditions that differed in terms of liquid water availability, particle concentration and grain size distributions. In our proposed model, a characteristic sequence of facies accreted from core to rim in the accretionary lapilli during passage through ash clouds generated by vent-derived plumes and pyroclastic density currents. The accretionary lapilli are mostly composed of smaller aggregates (ash clusters, ash pellets) and grew predominantly by accretion of already-formed aggregates, rather than by grain-by-grain accretion of individual particles. This finding is consistent with observations of rapid aggregate growth in volcanic plumes, suggesting a common evolutionary pathway for accretionary lapilli formation across diverse eruptions.

Ilopango caldera↗
USGS2022

A reappraisal of explosive–effusive silicic eruption dynamics: Syn-eruptive assembly of lava from the products of cryptic fragmentation

Silicic volcanic eruptions range in style from gently effusive to highly explosive, and may switch style unpredictably during a single eruption. Direct observations of subaerial rhyolitic eruptions (Chaiten 2008, Cordón Caulle 2011–2012, Chile) challenged long-standing paradigms of explosive and effusive eruptive styles and led to the formulation of new models of hybrid activity. However, the processes that govern such hybrid explosive–effusive activity remain poorly understood. Here, we bring together observations of the well-studied 2011–2012 Cordón Caulle eruption with new textural and petrologic data on erupted products, and video and still imagery of the eruption. We infer that all of the activity – explosive, effusive, and hybrid – was fed by explosive fragmentation at depth, and that effusive behaviour arose from sticking and sintering, in the shallow vent region, of the clastic products of deeper, cryptic fragmentation. We use a scaling approach to determine that there is sufficient time available, during emplacement, for diffusive pyroclast degassing and sintering to produce a degassed plug that occludes the shallow conduit, feeding clastogenic, apparently effusive, lava-like deposits. Based on evidence from Cordón Caulle, and from other similar eruptions, we further argue that hybrid explosive–effusive activity is driven by episodic gas-fracking of the occluding lava plug, fed by the underlying pressurized ash- and pyroclast-laden region. The presence of a pressurized pocket of ash-laden gas within the conduit provides a mechanism for generation of harmonic tremor, and for syn-eruptive laccolith intrusion, both of which were features of the Cordón Caulle eruption. We conclude that the cryptic fragmentation models is more consistent with available evidence than the prevailing model for effusion of silicic lava that assume coherent non-fragmental rise of magma from depth to the surface without wholesale explosive fragmentation.

Journal of Volcanology and Geothermal Research↗
USGS2022

Rock alteration mapping in and around fossil shallow intrusions at Mt. Ruapehu New Zealand with laboratory and aerial hyperspectral imaging

Diagnostic absorption features in hyperspectral data can be used to identify a specific mineral or mineral associations. However, it is unknown how accurate hyperspectral mapping can be for identifying alteration mineral compositions at the resolution required to describe structures such as fossil intrusions, or whether it can accurately quantify the alteration present. This study compared petrographic observation with visible, near-infrared (VNIR), and shortwave infrared (SWIR) hyperspectral remote sensing at laboratory- (centimetre-scale) and aerial- (metre-scale) scales to characterise the abundance of surface hydrothermal rock alteration in and around a shallow fossil intrusion on Pinnacle Ridge, Mt. Ruapehu, New Zealand. We classified a high-resolution aerial hyperspectral image to develop a new surface alteration map using Spectral Angle Mapper (SAM) algorithm. The petrographic thin-section and the laboratory and aerial hyperspectral imaging revealed a spectrum of hydrous alteration phases as indicated by the presence of an absorption feature at 2207 nm. Moderate correlation exists between the depth of the absorption feature at 2207 nm and the point counting-derived alteration percent values, indicating reliability of laboratory-based hyperspectral analytical methods. In contrast, aerial hyperspectral data failed to provide any clear correlations to field-mapped alteration using a band-depth approach, and we interpret this due to ‘oversampling’ of surface (supergene) alteration, spectral mixing, and sensor limitations (e.g., bandwidth, signal-to-noise ratio). The hyperspectral image-derived alteration map, created using supervised image classification, can loosely be translated to a geotechnical map where porosity and permeability play a major role in localizing hydrothermal fluid flow and the formation of alteration mineral associations.

Mt. Ruapehu↗
USGS2022

The source, fate, and transport of arsenic in the Yellowstone hydrothermal system - An overview

The Yellowstone Plateau Volcanic Field (YPVF) contains >10,000 thermal features including hot springs, pools, geysers, mud pots, and fumaroles with diverse chemical compositions. Arsenic (As) concentrations in YPVF thermal waters typically range from 0.005 to 4 mg/L, but an As concentration of 17 mg/L has been reported. Arsenic data from thermal springs, outflow drainages, rivers, and from volcanic rocks and silica sinter were used to identify the sources, characterize geochemical and microbial processes affecting As, and quantify As fluvial transport. Arsenic in YPVF thermal waters is mainly derived from high temperature leaching of rhyolites. Arsenic concentrations in thermal waters primarily depend on water type, which is controlled by boiling, evaporation, mixing, and mineral precipitation and dissolution. Springs with low As concentrations include acid-SO 4 (0.1 ± 0.1 mg/L), NH 4 -SO 4 rich (0.003 ± 0.007 mg/L), and dilute thermal waters (0.1 ± 0.1 mg/L); travertine-forming waters have moderate As concentrations (0.4 ± 0.2 mg/L); and neutral- Cl waters (1.2 ± 0.8 mg/L) common in the western portion of the Yellowstone Caldera and Cl-rich waters (1.9 ± 1.2 mg/L) primarily from Basins near the Caldera boundary have elevated As concentrations. Reduced As species (arsenite and thiolated-As species) are most prevalent near the orifice of hot springs, and then As rapidly oxidizes to arsenate along drainages. Previously published cultivation-based studies and metagenomic data from microbial communities inhabiting a variety of hot springs indicate a widespread distribution of arsenite oxidation and arsenate reduction capabilities among the hot springs. Widespread use and transformation of As by thermophilic microorganisms promotes more soluble and toxic forms. Most of the water discharged from thermal springs eventually ends up in a nearby river where As remains soluble and exhibits little attenuation during downstream transport. Since 2010, 183 ± 10 metric tons/year of As were transported from Yellowstone National Park (YNP) via rivers. The discharge from YPVF thermal features impairs river water quality whereby As concentrations exceed 10 μg/L for many rivers reaches within and downstream from YNP.

Idaho, Montana, Wyoming↗
USGS2023

Precaldera mafic magmatism at Long Valley, California: Magma-tectonic siting and incubation of the Great Rhyolite System

The iconic volcanic center at Long Valley has released ∼820 km 3 of rhyolite in at least 110 eruptions. From 2.2 Ma until 0.23 Ma, products were exclusively rhyolitic, and ∼ 700 km 3 were high-silica rhyolite severely depleted in Sr, Ba, and Eu. The rhyolitic interval was preceded by an interval from 3.9 to 2.6 Ma with numerous basalt-andesite-dacite eruptions accompanied by no rhyolite at all. We have now mapped the circumcaldera products of this interval, defined 107 eruptive units, characterized them all chemically and petrographically, and dated many by 40 Ar/ 39 Ar. Here we display and describe them by sector around the caldera, interpret the nature of the transcrustal magma system that eventuated in the 35-km-wide Long Valley granite-rhyolite pluton, and analyze regional tectonic factors that did or did not contribute to siting the system. Nine Miocene (12–6 Ma) eruptive units close to Long Valley were followed by a Pliocene flare-up that released >300 mafic eruptions in a SW–NE swath 170 km long, centered across the later site of Long Valley. The basalts and their fractionates are intraplate alkalic products dominated by a continental lithosphere that had long been fluxed by Mesozoic subduction. Tertiary arc volcanism had not impinged on the area of Long Valley. Volumes estimated for the 107 Neogene precaldera eruptive units (only 40 of which exceeded 0.1 km 3 ) total ∼ 27 km 3 ± 50%—only ∼3% of the subsequent volume of rhyolite erupted. Such a volume of high-silica rhyolite with ultra-low Sr and Eu is not a product of partial melting but requires as proximate parent a leucogranitic crystal mush that is itself the upper level of a long-lived plutonic reservoir that extends to the lower crust. The 27 km 3 of Neogene magma that erupted was a small contingent of the mantle-derived basaltic flux needed to energize (and contribute its fractionated melt to) a 30-km-deep compositionally graded crustal column, which culminated in ∼10,000 km 3 of granitoid mush from which 820 km 3 of Quaternary high-silica rhyolitic melt escaped and erupted. Pliocene basaltic eruptions ceased at ∼2.6 Ma, probably because the basaltic flux intensified sufficiently to render the mushy upper crust impenetrable. The 2.6–2.2 Ma quiescent interval represented culmination of thermal activation of the plutonic column and refinement of its leucogranitic mushy upper layer, from which extreme melts escaped for the next 2 Myr. The Pliocene mafic swath crosses the Sierran rangefront fault zone coincident with a left-stepping extensional reentrant that also began developing at ∼3 Ma. Moreover, Long Valley overlies a dextral offset in the initial Sr-isotope 0.706 line, which may reflect the rifted or attenuated edge of Proterozoic crust and mantle lithosphere. Concatenation of these three influences may account for siting of intensified edge-focused magmatism that produced the great Quaternary pluton.

California↗
USGS2023

Geochemistry and fluxes of gases from hydrothermal features at Newberry Volcano, Oregon, USA

We present the chemical and isotopic compositions of gases and fluxes of CO 2 from the hydrothermal features of Newberry Volcano, a large composite volcano located in Oregon's Cascade Range with a summit caldera that hosts two lakes, Paulina and East Lakes. Gas samples were collected from 1982 to 2021 from Paulina Hot Springs (PHS) on the shore of Paulina Lake, East Lake Hot Springs (ELHS) on the shore of East Lake, and Obsidian Flow Gas Seep (OFGS), an area of diffuse gas emissions. Surveys of CO 2 flux were conducted in 2020 at OFGS (1400 m 2 ) and East Lake (4.1 km 2 ). Gases from all three sites were CO 2 -rich (≥79 mol% in dry gas) but showed considerable compositional variability over time due to interaction with ground and surface water. An increase in H 2 S concentrations and decline in CO 2 /H 2 S ratios in ELHS gases coincided with a drop in East Lake water level from 1999 to 2021. ELHS and OFGS gases were high in CH 4 relative to PHS and the δ 13 C of CH 4 values for ELHS gases (−72.2 and − 63.6 ‰) reflected a predominantly biogenic origin. The dominant source of N 2 and Ar in PHS, ELHS, and OFGS samples was likely groundwater. Helium isotopic ratios (6.47 to 8.02 R c /R a ) support a persistent source of magmatic He beneath Newberry caldera and consistently high values measured at OFGS and PHS relative to ELHS suggest distinct fluid flow paths from depth to the surface features. The δ 13 C of CO 2 and CO 2 / 3 He values (−8.9 to −5.35 ‰ and 1.3 × 10 9 to 4.6 × 10 10 , respectively) measured in gases reflect contributions of CO 2 from both mantle and crustal sources. Measured CO 2 fluxes at OFGS and East Lake ranged from 1 to 8808 and < 1 to 364 g m −2 d −1 , respectively. A CO 2 emission rate of 0.5 t d −1 was calculated for OFGS. The CO 2 emission rate estimated for East Lake was 30 t d −1 and when compared to prior estimates, reflects steady-state lake degassing. An enhanced geochemical monitoring plan, including annual sampling of gases at ELHS, OFGS, and PHS for geochemical analysis, installation of a continuous lake-level monitoring station at East Lake, and annual CO 2 flux surveys at OFGS, would provide valuable background data and insights into any precursor volcanic activity. Integrating geochemical data with data from the real-time seismic and GPS network at Newberry Volcano could better resolve and interpret potential changes in its magma-hydrothermal system.

Oregon↗
USGS2023

Paleomagnetically defined brief lifespans for two large shield volcanoes in the Cascades Arc

Mafic to intermediate shield volcanoes with multi-cubic-kilometer eruptive volumes are common in the Cascades Volcanic Arc, but little is known about their eruptive histories as either singular or sustained episodes, or the total time required for their construction. Paleomagnetic data were collected from the lava flows of Ash Creek Butte (17 sites) and Crater Mountain (14 sites) in northern California; both volcanoes are large shields with total volumes of ∼11 km 3 each. Tightly clustered paleomagnetic results at both volcanoes, when coupled with analysis of geomagnetic secular variation, suggest that each edifice was built in only a few centuries, possibly in as little time as 50–90 years, indicating sub-century to century scale eruptive durations for two sizeable regional shield volcanoes within the Cascades Arc. These rapidly built shield volcanoes are substantially larger than typically defined monogenetic volcanoes, yet both are wholly formed within a single ‘episode’ at the limits of temporal resolution. Paleomagnetic methods provide a high-resolution tool that can be applied to understanding the tempo of regional volcanism in arcs.

California, Oregon, Washington↗
USGS2023

Dating individual zones in phenocrysts from the 2016–2017 eruption of Bogoslof volcano provides constraints on timescales of magmatic processes

We investigate the rates of magmatic processes using clinopyroxene diffusion chronometry on volcanic products erupted in August 2017 at the end of the 9-month eruption of Bogoslof volcano. The eruptive products contain plagioclase, clinopyroxene, and amphibole, all of which exhibit sharp chemical boundaries and are occasionally observed in multi-phase crystal clots with shared zoning boundaries across different mineral phases. At the shared boundaries in crystal clots, clinopyroxene and plagioclase continued to grow but abruptly changed composition from Mg# 81.7 ± 5.8 to 72.9 ± 3.0 and An 82.5±1.4 to An 61.3±5.7 , respectively. Additionally, the sharp boundary marks where amphibole became unstable and began forming a reaction rim. Synthesizing these observations, we were able to determine that the shared boundaries formed as a result of rapid decompression during magma ascent, followed by storage in a shallow cryptodome, where magma accumulated prior to erupting. In order to determine the timescales of magma ascent and subsequent crystal residence times, we applied diffusion chronometry on zoned clinopyroxene phenocrysts using Mg# concentrations at 1056 °C determined from Fe Ti oxide pairs. Our diffusion modeling results show that diffusion began at the stepwise boundaries in clinopyroxenes no more than180 days before the final explosive event. These results were then used to calculate crystal growth rates for shared plagioclase and amphibole rims, as shared zones in crystal clots indicate that the boundaries in all three phases formed contemporaneously. We calculate growth rates of plagioclase crystals (1.7 ± 0.99 × 10 −6 um/s) and amphibole reaction rims (2.8 ± 0.47 × 10 −6 um/s). The calculated natural growth rate of plagioclase was then used to constrain additional magmatic timescales from growth rate chronometry, results of which support our diffusion timescales. Our results indicate that the distinct boundaries in all three mineral phases formed due to ascent-driven decompression followed by shallow emplacement of mafic magma that occurred continually throughout the course of the eruption. By subtracting diffusion timescales from the date that the samples were erupted, the oldest crystal boundaries correspond to March 2017, seemingly correlating with increases in both seismicity and SO 2 emissions. These observations may suggest that our petrochronometric results are supported by interdisciplinary observations.

Alaska↗
USGS2023

The 2013−2020 seismic activity at Sabancaya Volcano (Peru): Long lasting unrest and eruption

Sabancaya volcano is the youngest and second most active volcano in Peru. It is part of the Ampato-Sabancaya volcanic complex which sits to the south of the ancient Hualca Hualca volcano and several frequently active faults, thus resulting in complex volcano-tectonic interactions. After 15 years of repose, in 2013, a series of 4 earthquakes with magnitude >4.5 occurred within 24 h, marking the beginning of a new episode of unrest. Several additional swarms of earthquakes occurred in the following years until magmatic eruptive activity started on 6 November 2016. This activity is ongoing as of this writing, with an average of 50 explosions per day. In this study, we present results of multiparametric monitoring of Sabancaya's activity observed during 2013–2020. Seismic data are used to create a one-dimensional seismic velocity model, to catalog, locate, and characterize earthquakes, to detect repeating earthquake families, and to monitor seismic velocity variations by ambient noise cross-correlation. These analyses are complemented by visual and remote sensing observations and ground deformation measurements. All monitored parameters showed significant changes on 6 November 2016, the day of eruption onset, thus dividing the eruptive activity into pre-eruptive and eruptive stages. The unrest is characterized by high levels of seismic activity with hundreds of events detected per day. Volcano-tectonic (VT) earthquakes were dominant during the pre-eruptive period while long-period (LP) events and explosions have been most numerous since the eruption onset. Earthquake locations highlight long-lasting seismogenic zones along multiple previously active regional faults, as well as along newly identified faults. This VT seismicity is mainly distributed in a sector from the northwest to the east of the volcanic complex at distances of up to 30 km from the crater. We focus our analysis on two eruptive episodes: the eruption onset and subsequent crater migration from south to north, and the increase of lava dome extrusion rate in 2019. Both episodes are accompanied by seismic velocity decreases of up to 0.2% and are preceded by a few weeks by bursts of distal VT activity, including numerous repeating earthquakes. These repeated events were located on several remote tectonic faults (5–25 km from the vent). We suggest that these phenomena could be due to the injection of a batch of magma in the deep reservoir and/or conduit, which would generate 1) a pressure wave propagating in the hydrothermal system, triggering the bursts of seismic activity and 2) slow rising of magma by melting old material filling the conduit that eventually produced the eruptive and dome growth acceleration events.

Andes Mountains, Sabancaya Volcano↗
USGS2023

Ash aggregate-rich pyroclastic density currents of the 431 CE Tierra Blanca Joven eruption, Ilopango caldera, El Salvador

The VEI 6, Tierra Blanca Joven pyroclastic sequence (30–90 km 3 DRE volume), erupted from Ilopango caldera, El Salvador, in 431 CE, is the product of one of the largest eruptions of the last two millennia. The eruption devastated Central America's Mayan civilization. The eruption began with a short-lived phase of ash and pumice fall deposition and transitioned to a ‘wet’ explosive phase during which pyroclastic density currents flowed >40 km from the caldera. Detailed field and sedimentological analyses are provided for the deposits of ash-aggregate-rich pyroclastic density currents generated during early phases of the eruption. The first phase of pyroclastic density current inundation incinerated forests and deposited up to 30 m of, non-welded, ash-rich ignimbrite in proximal regions, along with ash fall layers of co-ignimbrite origin. Following fallout of a thin layer of pumice and lithic lapilli, a second phase of pyroclastic density current inundation and co-ignimbrite ash fall commenced. A range of ash aggregate types is present in the pyroclastic density current deposits and interbedded co-ignimbrite ash fall layers. Whole and broken concentrically layered ash aggregates (accretionary lapilli) reach >50 vol% in some horizons within some beds in the pyroclastic density current deposits. The evidence indicates that the ash aggregates grew within overriding co-ignimbrite ash plumes and subsequently fell into ground-hugging currents. Our findings suggest that the aggregate-rich nature of the pyroclastic density current deposits originated through incorporation of lake water into eruptive plumes, which in turn triggered rapid, pervasive aggregation within ash clouds and co-ignimbrite plumes.

Ilopango caldera↗
USGS2023

Microgravity as a tool for eruption forecasting

Detection of gravity change over time has been used to better understand magmatic activity at volcanoes for decades, but the technique is not commonly applied to forecasting eruptions. In contrast, other tools, notably seismic, deformation, and gas monitoring have made exceptional strides in the past several decades and form the foundation for eruption forecasting, especially during the final buildup to an eruption. Reasons for this gap include the high cost and fragile nature of gravity instruments, and the ambiguous nature of many results. But this is changing. Instrumentation is becoming more robust and accurate, expenses may soon diminish thanks to technological advances, and the record of success in tracking subsurface mass change (either from magma or hydrothermal fluids) in volcanic areas is growing. Here we review how gravity change can be applied to forecasting volcanic eruptions across a variety of spatial and temporal scales. We argue that microgravity has untapped potential as a forecasting tool in three specific ways: constraining probabilistic assessments, detecting long-term mass change that may occur prior to the onset of vigorous seismicity and deformation, and identifying transient activity that indicates magma ascent or other changes that immediately precede new eruptions or changes in ongoing eruptions. As with any volcano-monitoring method, microgravity has strengths and weaknesses, but the varied forms of data collection—for instance, campaign versus continuous, and relative versus absolute—offer the potential to record a broad range of signals at volcanoes with a diversity of magmatic systems. The technique is currently underutilized; additional attention, investment, and application at more volcanoes could help to realize its promise.

Hawai'i↗
USGS2023

Eruption of stagnant lava from an inactive perched lava lake

Lava flow hazards are usually thought to end when the erupting vent becomes inactive, but this is not always the case. At Kīlauea in August 2014, a spiny ʻaʻā flow erupted from the levee of a crusted perched lava lake that had been inactive for a month, and the surface of the lava lake subsided as the flow advanced downslope over the following few days. Topography constructed from oblique aerial photographs using structure-from-motion (SfM) software shows that the volume of the flow (∼68,000 m 3 ) closely matches the volume of subsidence of the crusted lava lake (∼64,000 m 3 ). The similarity of these volumes, along with the textural characteristics of the lava, shows that the lava that fed the flow had been stored beneath the surface of the perched lava lake, and that the flow was not generated by reactivation of the vent. This extends the duration of the local lava flow hazard presented by perched lava lakes and similar flow field structures that store lava, such as rootless shields. The flow probably occurred because the density of the lava beneath the crusted surface of the perched lava lake increased through loss of gas bubbles until it was able to penetrate the less-dense levee, which was composed of relatively vesicular overflows. The flow is thus equivalent to the lava seeps described previously at Kīlauea and elsewhere. We present a simple physical model for the pressure change at the base of a densifying body of lava, which we apply to this case study, and which could be applied to similar scenarios elsewhere.

Hawaii↗
USGS2023

Chemistry and petrography of early 19th century basaltic andesites and basalts from the Kamakaiʻa Hills in the Southwest Rift Zone of Kīlauea volcano, Hawaiʻi

Kīlauea is a frequently active, open-system volcano on the Island of Hawaiʻi known for erupting olivine-dominated tholeiitic basalt compositions. On rare occasions it erupts more differentiated magmas (<1% of erupted volume), such as basaltic andesites and andesites, from its rift zones. These differentiated magmas offer an opportunity to understand better the petrology, magma storage, magma mixing, and eruptive triggers that occur in Kīlauea's rift zone reservoirs. This study focuses on an eruption from the Southwest Rift Zone of Kīlauea, which is dominantly basaltic andesite with subordinate basalt. This eruption originated at the Kamakaiʻa Hills during the early 19th century and has two eruptive phases: 1) an early ‘a‘ā phase that is primarily exposed in the eastern part of the flow field, with minor western lobes, and 2) a late pāhoehoe phase that makes up most of the western part of the flow field. The early ‘a‘ā phase covers at least 5.8 km 2 with an erupted volume of ∼150 × 10 6 m 3 and consists of uniform composition basaltic andesites with 3.72–4.15 wt% MgO over its ∼7 km flow length. The late pāhoehoe phase reached >10 km from its vent, covers an area of ∼7.1 km 2 , has a volume of ∼100 × 10 6 m 3 , and initially erupted basaltic andesite near its vent (4.50–5.64 wt% MgO extending to 3.8 km from vent) with channel and tube-fed basalt (6.21–12.38 wt% MgO sampled at >3.8 km from vent) emplaced during its waning stages. Most Kamakaiʻa Hills lavas are crystal-poor, containing ≤1.5% glomerocrysts and individual phenocrysts of plagioclase + clinopyroxene + Fe Ti oxides ± orthopyroxene, as well as olivine in lavas with >6 wt% MgO. Major-oxide and trace-element concentrations throughout the Kamakaiʻa Hills lavas demonstrate the involvement of three distinct magmatic processes. First, the basaltic andesites of the early ‘a‘ā phase are the products of fractionation of plagioclase + clinopyroxene + Fe Ti oxides ± orthopyroxene, indicative of magmas that have been stored in rift zone reservoirs for decades or longer. Second, the near-vent (within ∼400 m of vent) basaltic andesites of the late pāhoehoe phase yield chemical concentrations that indicate magma mixing with a more differentiated magma (of a similar evolved composition to basaltic andesites at ∼55–56 wt% SiO 2 and ∼3.4–4.1 wt% MgO that erupted in the lower East Rift Zone in 2018). Third, the progressively more mafic magma (containing olivine + plagioclase + clinopyroxene) that continued to erupt throughout the waning stages of activity suggests an eruptive triggering process whereby an intruding summit or uprift reservoir basalt overpressurized and forced out the stored, differentiated magma of the Kamakaiʻa Hills rift zone reservoir.

Hawaiʻi↗
USGS2024

Origins and nature of large explosive eruptions in the lower East Rift Zone of Kīlauea volcano, Hawaii: Insights from ash characterization and geochemistry

Several powerful explosive eruptions have taken place in the populated lower East Rift Zone of Kīlauea within the past ∼750 years. These have created distinctive landforms, including a tephra rim enclosing Puʻulena Crater immediately south of the Puna Geothermal Venture power station, a tuff cone at Kapoho Crater near the eastern cape of the Island of Hawaiʻi, and a set of littoral cones, the Sand Hill in Nānāwale, where the 1840 lava flow poured into the ocean. Kapoho Crater tuff cone is the largest of these recent pyroclastic features. Mineral, glass, and melt inclusion analyses of tuff cone ash and later fissure-related scoriaceous materials also found within the crater indicate slightly evolved basaltic magmas (1120–1130 °C) that are compositionally similar to parts of the effusive lower East Rift Zone eruptions in 1955 and 2018. Tuff cone magmas were stored at depths of ∼2.5–3.5 km and had pre-eruptive volatile contents (0.5–0.8 wt% H 2 O, 280–340 ppm CO 2 , 1400–1800 ppm S) similar to other Kīlauea eruptions (e.g., 1959, 1960), suggesting that internal magma properties were unlikely to account for the unusual explosiveness of this eruption. Tephra componentry, grain-size analyses, and field observations confirm that the cone grew during a phreatomagmatic eruption mostly of vitric ash, probably where a fissure opened across the coastline or shallow ocean floor nearby. Supporting this hypothesis is the identification of at least two genera of marine diatoms within tuff cone strata. Sand Hill littoral cone ash is also vitric like that of Kapoho Crater, but distinctly coarser with abundant fluidal ejecta represented. In contrast, the Puʻulena Crater eruption deposited lithic ash and related blocks with minor juvenile magmatic contribution; a phreatomagmatic eruption that was dominantly phreatic. Differences in eruption styles are related to unique mechanics that tephra analyses help us interpret. While powerful explosive eruptions in the lower East Rift Zone are rare, they present a definite future hazard for inhabitants in this part of Hawaii.

Hawaii↗
USGS2024

Event classification, seismicity, and eruption forecasting at Great Sitkin Volcano, Alaska: 1999–2023

The frequency content of volcanogenic seismicity is often used to classify events and their spatial and temporal progression is then used to map subsurface volcanic processes. The progression of volcano-seismic events and associated source processes also plays a critical role in eruption forecasting. Here we develop and evaluate a computerized methodology for characterizing volcano-seismic event types using Frequency Index and Average Peak Frequency. We apply and test this technique at Great Sitkin Volcano, Alaska, classifying over 9000 hypocenters between 1999 and 2023. This 24-year time span covers periods of seismic quiescence, earthquake activity on nearby tectonic (bookshelf) faults, precursory unrest from 2016 to 2021, and the explosive onset in May 2021 of the ongoing effusive eruption. We use the spatial and temporal evolution of classified event types to map the active volcanic and tectonic processes, develop a conceptual model of the subsurface magmatic system, and perform a retrospective analysis of eruption forecasts at Great Sitkin Volcano between 2016 and the present. The classification and progression of hypocenters suggests the subsurface Great Sitkin Volcano magmatic system consists of a mid- to lower- crustal source zone between 10 and 40 km depth and an upper crustal magma storage area between −1 and 10 km depth (hypocenter depth is referenced to sea level and negative depths reflect height above sea level). The earliest precursors occurred in July 2016 and consisted of deep long-period and volcano-tectonic earthquakes at mid-crustal depths suggesting the subsequent unrest and eruption were triggered by a deeper intrusion of magma. This mid-crustal seismic activity was immediately followed by the onset upper-crustal long-period events and volcano-tectonic earthquakes VTs suggesting a strong linkage between the shallow and deeper portions of the magmatic system. The upper crustal area was likely capped by the 1974 lava dome until the magmatic explosion on May 26, 2021.

Alaska↗
USGS2024

Unique seismic and eruption precursors to the 1996 and ongoing magmatic eruptions of Popocatépetl: Coupled and fluidized bed events

We describe three unique types of seismicity at Popocatépetl volcano that accompanied the initial vent-clearing eruptive activity in December 1994 through the eruption of the first two domes in 1996. We identify and describe two types of coupled events, 1) spasmodic burst coupled events, a burst of volcano tectonic (VT) events coupled with a large eruptive explosion (21 December 1994 and 5 March 1996), and 2) explosion-couplet coupled events, a pair of events with a deeper first event and whose second event correlates with a small gas emission or explosion. Explosion-couplets occurred with the onset of magmatic ash dominated eruptions and their properties are very useful for forecasting magmatic eruptions at Popocatépetl volcano. Measurable quantities including the time between the first and second phase, daily numbers of events, and the amplitude ratio between the second and first phase systematically changed as the first two domes approached the surface between March and June 1996. Interevent times decreased from many tens of seconds to a few seconds, while amplitude ratios increased from about 2 to 10 or more. Event numbers increased prior to and during initial dome extrusion. These changes were used to forecast the eruption of the first two domes. We also report on another unusual type of low frequency seismicity that accompanied emissions called fluidized bed events. Fluidized bed events occurred following the initial eruption on 21 December 1994, beginning mid-January 1995, when SO 2 emissions were initially elevated, and the conduit was open. Fluidized bed events began to wane in late March through the end of May 1995 as did SO 2 emissions. Consistent gas, visual and seismic observations by Centro Nacional de Prevención de Desastres (CENAPRED) enabled the direct correlation of the gas, magma, and seismic phenomena.

Popocatépetl↗