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 289 records · Page 16Linked to original sources
USGS2020

Combining ash analyses with remote sensing to identify juvenile magma involvement and fragmentation mechanisms during the 2018/19 small eruption of Peteroa volcano (Southern Andes)

The Planchón Peteroa Volcanic Complex (PPVC) is located on the border of Chile and Argentina, and is one of the most active volcanic systems in the Andes. Holocene activity has included magma-water interaction with an evolving series of crater lakes, mainly sourced from Peteroa volcano. This study examines data from the 2018/19 eruption, together with the volcanic history of the PPVC, to elucidate the complex interplay between magmatic activity and summit water and ice. From February 2016 to mid-2019, three seismic swarms occurred in the PPVC, preceding the explosive eruption from September 2018 to April 2019. The activity originated from a small vent nested within the easternmost crater, the most active portion of the complex (Peteroa). The explosions interacted with a crater lake, producing ash plumes up to 2 km above the crater and building a small tephra cone. To investigate the eruption mechanisms, we performed remote sensing analysis of plume dispersal, thermal anomalies and ground deformation, and characterized the volcanic products, including grain size, componentry, morphology, internal textures, composition and mineralogy. Our results suggest that the precursory seismicity beginning in 2016 was related to the intrusion of a new magma batch that reached the surface during the 2018/19 eruption. The eruption was also preceded by thermal anomalies, geomorphic changes and increased hydrothermal activity at the surface, though without any ground deformation recognized through radar interferometry (InSAR). The eruption initially produced predominantly recycled ash (phreatic activity), then evolved to increasing proportions of juvenile magma (phreatomagmatic) by April 2019. The juvenile clasts had a trachyandesite composition (~59 wt% SiO 2 ), with vesicular and dense scoria containing plagioclase and pyroxene. The ash surfaces show external quenching cracks and step fractures consistent with phreatomagmatic fragmentation within the active crater lake. Textural characteristics also point to a slowly ascending batch of magma that was relatively viscous by the time it interacted with water in the crater lake. Notably, these juvenile particles are distinctive from the pre-2018 products. Ash erupted from 2010/11 did not contain recognizable juvenile material, and is inferred to have been a mainly phreatic eruption. Our findings suggest that the interplay between phreatic and phreatomagmatic eruptions fed by small magma batches intruding at shallow levels characterize much of the eruptive behavior of the PPVC during the last three decades. Multi-parametric assessment is a powerful tool to discriminate between phreatic and phreatomagmatic eruptions.

Southern Andes↗
USGS2020

Dating silica sinter (geyserite): A cautionary tale

We describe a new effort to date hydrothermal silica sinter deposits (geyserite) from the Upper Geyser Basin of Yellowstone National Park using 14C of co-deposited organic matter, U-series and cosmogenic 10Be methods. A majority of the samples were collected from stratigraphic sections, mainly at Riverside, Giant, and Castle Geysers. Ages obtained from 41 14C analyses range from modern to 12.1 cal ka BP. Nearly all the 14C ages show inconsistencies with their stratigraphic positions, and several replicate 14C analyses from the same sample result in significantly different ages. The δ13C values of the organic material in the sinter range from -26.6‰ to -12.7‰. The more enriched values are attributed to microbial fixation of dissolved inorganic carbon (DIC), which has heavier δ13C values and is 14C-depleted relative to atmospheric CO2, leading to apparent older ages. U-series analyses on 4 samples yielded ages between 2.2 and 7.4 ka. Large 230Th/U age uncertainties in the sinter, due to low uranium concentrations along with elevated 232Th and associated initial 230Th, make these ages imprecise for use on Holocene deposits. A single cosmogenic 10Be exposure age of 596±18 ka is considerably older than the age of underlying rhyolite and is thus unreliable. This apparent old age results from contamination by meteoric 10Be trapped in the opal that overprints the very small amount of cosmogenic 10Be. By presenting the problems we encountered and discussing their probable cause, this paper highlights the difficulty in obtaining reliable, high-precision geochronological data necessary to use sinter deposits as paleoenvironmental and paleo-hydrothermal archives.

Wyoming↗
USGS2020

Sources, fate, and flux of riverine solutes in the Southwest Yellowstone Plateau Volcanic Field, USA

Since the 1970s, temporal variations of hydrothermal discharge and thermal output from the numerous hydrothermal features in the Yellowstone Plateau Volcanic Field (YPVF) have been studied by measuring the chloride flux in the major rivers. In this study, the sources, fate, and flux of solutes in the Fall River and its major tributaries, in southwest Yellowstone National Park, were determined. The considerable precipitation in southwest YPVF and high groundwater flow through Quaternary rhyolites results in river solute fluxes that originate from shallow non-thermal groundwater and deep-thermal water. Specific conductance serves as a surrogate measure for thirteen riverine solute concentrations. Combining continuous 15-minute specific conductance and discharge data, the annual chloride, arsenic, fluoride, and silica fluxes from the Fall River were determined to be 11%, 5%, 25%, and 19% of the total flux exiting YPVF. Approximately 11% of the Fall River chloride flux is from non-thermal waters, which is larger than the previous estimate of 4 to 6%. Furthermore, a large proportion of fluoride and silica in the Fall River are derived from water-rock interaction in the shallow non-thermal groundwater system and the non-thermal weathering rate (30 ± 2 t/yr·km 2 ) is higher than other rivers draining the Yellowstone caldera. Consequently, 73 ± 3% of the annual total dissolved solid flux in the Fall River is from thermal sources. Synoptic sampling of river water and discharge measurements was performed during low-flow conditions that allowed for the determination of solute sources and their downstream fate. It was determined that chloride, sodium, arsenic, rubidium, lithium, and boron are primarily (>89%) associated with thermal waters and the Bechler River is the primary source of most hydrothermal solutes in the Fall River, but the major source of arsenic is Boundary Creek. Using the chloride inventory method, the thermal water discharge from several thermal areas was also determined.

Wyoming↗
USGS2020

Lava–water interaction and hydrothermal activity within the 2014–2015 Holuhraun Lava Flow Field, Iceland

Lava that erupted during the 2014–2015 Holuhraun eruption in Iceland flowed into a proglacial river system, resulting in aqueous cooling of the lava and an ephemeral hydrothermal system. We carried out a monitoring study of this system from 2015 to 2018 to document the cooling of the lava over this time, using thermocouple measurements and data-logging sensors. The heat loss rate from advection through this hydrothermal system in August 2015 was ~5.5 × 10 8 W; since eruption, aqueous cooling likely accounted for ~1% of the total heat loss from the lava. This estimate excludes steam losses from fumaroles as well as any groundwater that was not released to the surface, and thus is a lower bound. Near the terminus of the flow, advection of heat by flowing water may have locally accounted for tens of percent of the total cooling of that part of the flow. Our data quantify the importance of water cooling for this lava flow and can be compared with models to better understand lava–water interactions more generally. We also provide detailed methods for simple, low-cost monitoring of similar instances in the future.

Holuhraun Lava Flow Field↗
USGS2020

Understanding the storage conditions and fluctuating eruption style of a young monogenetic volcano: Blue Lake crater (<3 ka), High Cascades, Oregon

Blue Lake crater (<3 ka) is monogenetic volcano that produced one of the youngest eruptions in the central Oregon Cascades. Understanding monogenetic volcano behavior – from storage through eruption – is imperative in planning for future eruptions. Here we combine physical volcanology and geochemistry to determine the pre-eruptive storage conditions, ascent rate, eruption style, and deposit distribution of this young eruption. We find that the eruption of Blue Lake was initially phreatomagmatic, producing lithic-rich fall deposits and thin surge deposits and excavating the maar crater, before transitioning rapidly to a final voluminous magmatic-volatile driven explosive eruption. The mapped fall deposit has an estimated volume of 3.9 × 10 7 m 3 (2.2 × 10 7 m 3 DRE) which suggests a VEI of 3. Although similar in magnitude (as measured by fall deposit volume) to many other recent cinder cone eruptions in the Cascades, the Blue Lake crater eruption lacks an effusive phase. The absence of lava flows may reflect the lack of evidence for syn-eruptive magma storage at shallow levels. Indeed, corrected volatile contents of olivine-hosted melt inclusions (2.9–4.2 wt% H 2 O, 910–1330 ppm CO 2 ) are strikingly uniform and indicate storage and crystallization at a restricted pressure range (average ~ 235 MPa), equating to a depth of ~8.6 km. Melt inclusion geochemistry indicates that the basaltic andesite magma cooled and crystallized ~25% during storage at this pressure. Crystals in the Blue Lake magma show evidence of mixing with, or entrainment in, a more evolved magma. Feldspar crystals have large An-rich cores (An 80–85 ) and abrupt An-poor rims (An 60–70 ); olivine crystals have large, broad cores (~Fo 82–84 ) and thin rims with lower Fo and NiO contents. Diffusion modeling of olivine zoning suggests that an intrusion event occurred ~10–60 days prior to eruption. Diffusive loss of H + from melt inclusions was minimal (<1.3 wt% H 2 O) during magma ascent, from which we calculate minimum ascent times from 235 MPa of <1 day. Many inclusions indicate ascent times of <3 h, corresponding to ascent rates of ~1 to >13 m/s. This study illustrates the pre-eruptive and eruptive complexities of monogenetic volcanoes and highlights the minimal warning that may precede future eruptions.

Oregon↗
USGS2021

Correcting the historical record for Kīlauea Volcano's 1832, 1868, and 1877 summit eruptions

Three fissure eruptions are known to have occurred along the northeastern edge of Kīlauea's summit caldera in the 19th century—in the years 1832, 1868, and 1877. Modern portrayal of these eruptions on maps and in written sources indicates that the 1832 eruption was from a fissure on the side of the Poliokeawe scarp south of Byron Ledge, the 1868 eruption was from a fissure on the southern wall of Kīlauea Iki Crater and fed a lava flow that covered the bottom of that crater, and the eruption in 1877 occurred on the floor of Keanakākoʻi Crater, as well as from a fissure of uncertain location on the east wall of the caldera below Byron Ledge. New geologic mapping and a review of historical documents and maps contradict these views. We find, instead, that: (1) the 1832 eruption discharged from a fissure on Byron Ledge (not Poliokeawe scarp), from another fissure on the southwestern wall of Kīlauea Iki Crater, and from at least one fissure along the east side of Kīlauea caldera below Byron Ledge; (2) the 1868 lava erupted through the floor of Kīlauea Iki Crater, not from a fissure in its southwestern wall; and (3) the 1877 lava erupted from Kīlauea Iki Crater's mid-wall fissure (until now believed to have opened in 1868), from the fissure previously assigned an 1832 date on Poliokeawe escarpment, and from a precisely relocated vent on the northeastern wall of the caldera. Finally, no conclusive first-hand accounts of the late 19th century eruption in Keanakākoʻi Crater were identified, leaving in doubt the often-inferred 1877 date for this event. Possible alternative dates include 1868, 1879, and 1881.

Hawaii↗
USGS2021

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

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.

Alaska↗
USGS2021

Insight into the May 2015 summit inflation event at Kīlauea Volcano, Hawai‘i

We use ground and space geodetic data to study surface deformation at Kīlauea Volcano from January to September 2015. This period includes an episode of heightened activity in April and May 2015 that culminated in a magmatic intrusion beneath the volcano's summit. The data set consists of Global Navigation Satellite System (GNSS), tilt, visual and seismic time series along with 25 descending and 15 ascending acquisitions of the Sentinel-1 satellite. We identify four different stages of surface deformation and volcanic activity, which we attribute to pressure changes and the movement of magma in response to an imbalance between magma supply and withdrawal in the shallow plumbing system, eventually leading to an intrusion beneath the summit area. In particular, we model the deformation as due to pressure changes in two subsurface magma bodies: the Halema‘uma‘u Reservoir (HMMR) and South Caldera Reservoir (SCR). The SCR was best described by an ellipsoidal source at 2.8 (2.65–3.07 at 95% confidence) km depth below the south caldera region. The HMMR was modeled as a point source located just east of Halema‘uma‘u crater at 1.5 (0.95–2.62) km depth. We suggest that a short-term increase in the magma supply rate to the volcano is a potential mechanisms for the intrusion, although other factors, like the filling of available void space or a reduced efficiency of magma transport through the volcano's East Rift Zone, may also play a role.

Hawaii↗
USGS2021

Zircon geochronology and geochemistry of Quaternary rhyolite domes of the Coso volcanic field, Inyo County, California

The Quaternary Coso volcanic field (CVF) is a compositionally bimodal volcanic field located within a releasing bend along the eastern range-front Sierra Nevada fault zone in California's southern Owens Valley. The erupted products of CVF silicic magmatism since ~1 Ma comprise 38 high-silica rhyolite domes, with the volumetric majority (~99%) of rhyolite emplaced within the past ~300 ka. The CVF hosts an economically important geothermal field driven by heat associated with a shallow (~5 km) igneous intrusion. The CVF is potentially an immature analog to the nearby Long Valley system, which culminated in generation and eruption of the voluminous and widespread Bishop Tuff. As such, the CVF represents a considerable volcanic hazard, making a detailed understanding of the eruptive history and pre-eruptive conditions of the system critically important. We present uranium-series isochron dates from zircon ± allanite crystal surfaces and zircon trace element geochemical data on the youngest 17 rhyolite domes at Coso, which represent ~60% (by volume) of the silicic magma erupted by the system. These data suggest: (1) a shorter emplacement duration than previously recognized for these domes, with a duration of 20 ± 5 ka; (2) 4 shorter-duration eruption pulses within this interval, all of which occur during the marine isotope stage (MIS) 5 interglacial period; (3) an uptick in the volume of CVF magma erupted between ~200 ka and ~ 78 ka relative to that emplaced over the lifetime of the system; (4) near-coeval eruption of geochemically distinct magma in close geographic proximity, either sourced from different portions of the same magma system at depth or from discrete, uncommunicating bodies; (5) ambiguity with respect to whether or not CVF magmatism is time-predictable, as previously suggested, or erupted as a series of punctuated episodes; (6) no rhyolite volcanism in the past ~78 kyr.

California↗
USGS2021

The Independent Volcanic Eruption Source Parameter Archive (IVESPA, version 1.0): A new observational database to support explosive eruptive column model validation and development

Eruptive column models are powerful tools for investigating the transport of volcanic gas and ash, reconstructing past explosive eruptions, and simulating future hazards. However, the evaluation of these models is challenging as it requires independent estimates of the main model inputs (e.g. mass eruption rate) and outputs (e.g. column height). There exists no database of independently estimated eruption source parameters (ESPs) that is extensive, standardized, maintained, and consensus-based. This paper introduces the Independent Volcanic Eruption Source Parameter Archive (IVESPA, ivespa.co.uk), a community effort endorsed by the International Association of Volcanology and Chemistry of the Earth’s Interior (IAVCEI) Commission on Tephra Hazard Modelling. We compiled data for 134 explosive eruptive events, spanning the 1902-2016 period, with independent estimates of: i) total erupted mass of fall deposits; ii) duration; iii) eruption column height; and iv) atmospheric conditions. Crucially, we distinguish plume top versus umbrella spreading height, and the height of ash versus sulphur dioxide injection. All parameter values provided have been vetted independently by at least two experts. Uncertainties are quantified systematically, including flags to describe the degree of interpretation of the literature required for each estimate. IVESPA also includes a range of additional parameters such as total grain size distribution, eruption style, morphology of the plume (weak versus strong), and mass contribution from pyroclastic density currents, where available. We discuss the future developments and potential applications of IVESPA and make recommendations for reporting ESPs to maximize their usability across different applications. IVESPA covers an unprecedented range of ESPs and can therefore be used to evaluate and develop eruptive column models across a wide range of conditions using a standardized dataset.

Journal of Volcanology and Geothermal Research↗
USGS2021

Long-term year-round observations of magmatic CO2 emissions on Mammoth Mountain, California, USA

Diffuse emission of magmatic CO 2 is one of the main indicators of volcanic unrest at Mammoth Mountain, but the presence of deep seasonal snowpack at the site has hindered year-round CO 2 flux observations. A permanent eddy covariance station was established at the largest area of diffuse CO 2 degassing on Mammoth Mountain (Horseshoe Lake tree kill) that measured CO 2 fluxes ( F c ) and meteorological parameters on a half-hourly basis. From July 22, 2014 to May 24, 2020, F c ranged from −35 to 10,546 g m −2 d −1 . F c decreased on average by 53% over the study period, tracking the long-term decline in CO 2 emissions following the last major increase that occurred at the Horseshoe Lake tree kill area from 2009 to 2011. Statistical and spectral analyses were applied to the F c and ancillary meteorological parameter time series to understand (1) relationships between these parameters, (2) their dominant periodicities, and (3) changes in F c that may be unexplained by meteorological forcing. Variations in detrended F c ( F cdt ) were most strongly correlated with wind direction and atmospheric temperature, followed by atmospheric pressure on diurnal to annual time scales, but wind direction likely exerted the most direct control on F cdt . Comparison of the smoothed (180-d span) F cdt time series to the time series of average-daily snow water equivalent measured ~1 km away suggested that snowpack may have suppressed CO 2 emissions. No evidence of a change in CO 2 emissions related to the last major seismic swarm beneath Mammoth Mountain on February 2–18, 2014 was observed.

California↗
USGS2021

Evaluating the state-of-the-art in remote volcanic eruption characterization Part I: Raikoke volcano, Kuril Islands

Raikoke, a small, unmonitored volcano in the Kuril Islands, erupted in June 2019. We integrate data from satellites (including Sentinel-2, TROPOMI, MODIS, Himawari-8), the International Monitoring System (IMS) infrasound network, and global lightning detection network (GLD360) with information from local authorities and social media to retrospectively characterize the eruptive sequence and improve understanding of the pre-, syn- and post- eruptive behavior. We observe six infrasound pulses beginning on 21 June at 17:49:55 UTC as well as the main Plinian phase on 21 June at 22:29 UTC. Each pulse is tracked in space and time using lightning and satellite imagery as the plumes drift eastward. Post-eruption visible satellite imagery shows expansion of the island's surface area, an increase in crater size, and a possibly-linked algal bloom south of the island. We use thermal satellite imagery and plume modeling to estimate plume height at 10–12 km asl and 1.5–2 × 106 kg/s mass eruption rate. Remote infrasound data provide insight into syn-eruptive changes in eruption intensity. Our analysis illustrates the value of interdisciplinary analyses of remote data to illuminate eruptive processes. However, our inability to identify deformation, pre-eruptive outgassing, and thermal signals, which may reflect the relatively short duration (~12 h) of the eruption and minimal land area around the volcano and/or the character of closed-system eruptions, highlights current limitations in the application of remote sensing for eruption detection and characterization.

Hokkaido, Sakhalin Oblast↗
USGS2021

The structure and volume of large geysers in Yellowstone National Park, USA and the mineralogy and chemistry of their silica sinter deposits

Siliceous sinter is formed by biogenic and abiogenic opal deposition around hot springs and geysers. Using Structure-from-Motion photogrammetry we generated three-dimensional models of Giant and Castle Geysers from the Upper Geyser Basin of Yellowstone National Park. We use these models to calculate an approximate mass of sinter for each (~2 and ~ 5 kton, respectively) and estimate a range of plausible long-term deposition rates for Castle Geyser (470 to 940 kg·yr −1 ). We estimate ~2% of the silica discharged from Castle Geyser is deposited as sinter in the cone and proximal terraces. We collected 15 sinter samples following the stratigraphy of each geyser from an older terrace to a younger cone and examined them using a variety of analytical methods. We find that young opaline sinter with a water content of <12 wt% (from loss on ignition) contains higher concentrations of major and trace elements, notably As, Sb, Rb, Ga and Cs, relative to older dehydrated sinter. Rare earth element (REE) concentrations in sinter are 2–3 orders of magnitude higher than in the thermal water from which they are deposited. Sinter deposits are enriched in light REE, Gd and Yb when normalized to concentrations in thermal water and enriched in Eu, Tm, and Yb when normalized to the underlying rhyolite. Sinter samples with the highest REE concentrations are also enriched in organic material, implying either microbial uptake of REE, or that organic molecules are efficient ligands that form metal complexes.

Wyoming↗
USGS2021

New insights on faulting and intrusion processes during the June 2007, East Rift Zone eruption of Kilauea volcano, Hawai'i

The East Rift Zone (ERZ) of Kīlauea Volcano, Hawai'i, represents one of the most volcanically active regions in the world. The 2007 Father's Day (FD) dike intrusion, eruption, and accompanying slow-slip event (SSE) has been previously modeled using geodetic data to constrain the geometry of the intrusion and the timing and magnitude of the SSE. Here, we perform inversions of three interferometric synthetic aperture radar (InSAR) datasets and a new intensity offset tracking dataset to assess the effect of integrating intensity cross-correlation offsets into inversion problems and explore additional potential models for the intrusion geometry of the FD event based on this additional data. The overall lowest misfit single Okada model for all datasets opens 2.3 m, strikes 73 degrees while dipping sub-vertically at 83 degrees, and extends approximately 2.9 km to the ENE and 2.4 km downdip. The differences are minor between complex en-echelon distributed Okada and decollement model of (Montgomery-Brown et al., 2010) or 3D-MBEM breaching models including multiple surface breaches and free-slipping decollement movement. Finally, we examine the static Coulomb stress changes for the proposed decollement fault created by our preferred model and a representative model of deep rift opening and find that deep rift zones dilation, not shallow ERZ intrusions, are likely modulating slip on the decollement.

Hawaii↗
USGS2022

Characterizing unrest: A retrospective look at 20 years of gas emissions and seismicity at Iliamna Volcano, Alaska

Episodes of unrest are not as well documented as eruptions at most volcanoes globally. Iliamna is an andesitic stratovolcano in the Cook Inlet of Alaska that has experienced several episodes of unrest. Unrest in 1996 was previously studied. Here we present data from a minor period of unrest between 2002 and 2006, and a more significant period in 2012. None of the episodes led to an eruption. A dike intrusion was suggested for the 1996 unrest based on increases in gas emissions and seismic analysis. The 2002–2006 period was characterized by a slight increase in the rate of seismicity to 13 events per day and was particularly notable due to an increase in deep long period (DLP) seismic events between 15 and 37 km that were not observed at other times. This period also included one airborne gas measurement with and elevated CO 2 /SO 2 molar ratio (17). In 2012, Iliamna unrest was characterized by significantly elevated gas emissions (up to 582 t/d SO 2 and 1385 t/d CO 2 ) and up to 49 located earthquakes per day (M > 0), and was remarkably similar to the 1996 unrest. Differences in the observed evolution of the CO 2 /SO 2 gas ratio in 2012 (2.2–4) compared to that in 1996 (up to 18) suggests that no new deep magma was involved in 2012, however this does not preclude the movement of a previously intruded magma. A months-long increase in the SO 2 /H 2 S molar ratio from 8 to 17 during the peak of the activity could reflect a temperature increase on the order of 10–30 °C of the emitted gas. Compared to pre-eruptive unrest at other Cook Inlet volcanoes, Iliamna unrest in 2012 differed in that gas emissions were < 1500 t/d and seismicity lacked a rapidly escalating sequence of earthquakes and volcanic tremor, which is normally observed in the hours to days before eruption. The observation of DLPs, the fact that Iliamna produces moderately elevated degassing over decadal timeframes, and the persistent dominance of SO 2 over H 2 S, suggests that periodic input of fresh magma from the lower crust sustains the shallower magmatic system over time, which sets it apart from neighboring volcanoes in the Cook Inlet that show minimal activity between eruptions. Various scenarios could explain why Iliamna did not proceed to eruption in 2012. Finally, we present criteria by which monitoring data may suggest an increased likelihood of eruption at Iliamna in the future.

Alaska↗
USGS2022

Olivine and glass chemistry record cycles of plumbing system recovery after summit collapse events at Kīlauea Volcano, Hawai‘i

The eruptive activity of Kīlauea Volcano (Hawai‘i) in the past 2500 years has alternated between centuries-long periods dominated either by explosive or effusive eruptions. The onset of explosive periods appears to be marked by caldera collapse events at the volcano's summit accompanied by draining of Kīlauea's magmatic plumbing system. Here we leverage >1800 olivine forsterite (Fo) contents, >900 glass MgO contents, and estimated magma supply rates from the past six centuries to describe the relationships between summit collapse and the composition of erupted material. On a first order basis, the major element chemistry of the centuries-long eruptive periods largely originates from fundamental differences between fractional crystallization of shallowly stored magmas during high-supply effusive-dominated periods versus little evolution of mafic recharge magmas during low-supply explosive-dominated periods. The modern effusive period (1820s-present) is dominated by relatively evolved olivine forsterite contents (Fo81–82) for Kīlauea, which is interpreted to reflect a buffered crustal reservoir system in which shallow storage and fractional crystallization control the composition of magmas. In contrast, olivine crystals from the explosive Keanakāko‘i Tephra (1500 - early 1800s C.E.) are dominated by higher olivine forsterite contents (Fo89) which are interpreted to reflect more primitive compositions, are correlated with glass MgO compositions extending to high values (e.g.,11.0 wt%), and show signs of magma mixing (zoned olivine, bimodal Fo populations). These signatures reflect a disrupted reservoir system in which high-MgO recharge melts mix with melts left over from draining of the shallow (<5 km) magma plumbing. Superimposed on these explosive-effusive periods are three decades- to centuries long periods of progressively evolving olivine and glass compositions. Eruptions that occur after caldera collapse in ~1500C.E. and smaller scale crater collapse events in 1790 (inferred) and 1924 have heterogeneous olivine populations dominated by ≥Fo88 and typically high MgO glasses. These compositions reflect inefficient mixing of stored and primitive recharge magmas after the disruption of the shallow plumbing system. After these collapses, olivine Fo and glass MgO subsequently evolve to <Fo82 and <7.0 wt% compositions, reflecting the recovery of the crustal plumbing system to an end-member system state characterized by efficient mixing of recharge and stored magmas that serve to buffer the shallow magma reservoirs. These evolved signatures suggest that a mature and buffered reservoir system may be a key condition for significant disruptions of volcanic plumbing systems. Plumbing system recovery is slower following large-scale caldera collapse (hundreds of years) compared to recovery following smaller crater collapse (tens of years), which may be modulated by differences in magma supply rates. Following the 2018 crater collapse olivine populations have high-Fo but glasses are low MgO, suggesting that this collapse might have disrupted shallow magma pathways but not strongly impacted the reservoir(s). Ultimately, olivine and glass major element chemistry record the impacts of caldera and smaller but significant summit crater collapse events at Kīlauea and could be used to provide a framework for better characterizing long-term volcano evolution in Hawai‘i and shield volcanoes elsewhere.

Hawaii↗
USGS2022

Volcanic unrest at Nevados de Chillán (Southern Andean Volcanic Zone) from January 2019 to November 2020, imaged by DInSAR

The volcanic complex of Nevados de Chillán, located in the Southern Volcanic Zone (SVZ) of the Andes, has been active for the past 640 ± 20 ka. Its volcanic activity includes dome forming eruptions, explosive events, and lava flows. The most recent eruption cycle started in January 2016. We employ DInSAR time-series from Sentinel-1 data to investigate the unrest episode from January 2019 to November 2020. Two distinct periods of unrest are recognized in the time series. The first period (from January to October 2019) coincides with explosive events, dome growth inside the active crater, and a decrease in seismic activity but does not present a significant deformation. The second period (October 2019 to November 2020) is characterized by a displacement towards the sensor's line-of-sight of 100–120 mm. The observed surface deformation is compatible with an inflation source approximately 1.5 km south-southwest of the present active vent, at 5.5 ± 0.5 km depth from the surface, and with a volume change of 0.044 ± 0.014 km 3 . The most likely explanation for the observed inflation of Nevados de Chillan is the intrusion of magma in a reservoir feeding the current eruption cycle.

Nevados de Chillan↗
USGS2022

Outgassing through magmatic fractures enables effusive eruption of silicic magma

Several mechanisms have been proposed to allow highly viscous silicic magma to outgas efficiently enough to erupt effusively. There is increasing evidence that challenges the classic foam-collapse model in which gas escapes through permeable bubble networks, and instead suggests that magmatic fracturing and/or accompanying localized fragmentation and welding within the conduit play an important role in outgassing. The 2011–2012 eruption at Cordón Caulle volcano, Chile, provides direct observations of the role of magmatic fractures. This eruption exhibited a months-long hybrid phase, in which rhyolitic lava extrusion was accompanied by vigorous gas-and-tephra venting through fractures in the lava dome surface. Some of these fractures were preserved as tuffisites (tephra-filled veins) in erupted lava and bombs. We integrate constraints from petrologic analyses of erupted products and video analyses of gas-and-tephra venting to construct a model for magma ascent in a conduit. The one-dimensional, two-phase, steady-state model considers outgassing through deforming permeable bubble networks, magmatic fractures, and adjacent wall rock. Simulations for a range of plausible magma ascent conditions indicate that the eruption of low-porosity lava observed at Cordón Caulle volcano occurs because of significant gas flux through fracture networks in the upper conduit. This modeling emphasizes the important role that outgassing through magmatic fractures plays in sustaining effusive or hybrid eruptions of silicic magma and in facilitating explosive-effusive transitions.

Puyehue Cordon Caulle volcanic complex↗