SEARCH · Search USGS

Results for “Developments in Volcanology”

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

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4Linked to original sources

USGS2026

Lava flow emplacement dynamics during the Mauna Loa 2022 eruption

On November 27, 2022, Mauna Loa erupted for the first time in 38 years, threatening a major highway on the Island of Hawaiʻi. The opening phases of the eruption were focused at the volcano’s summit, before magma propagated into the upper Northeast Rift Zone to erupt at four fissure segments that localized at fissure 3 by December 2. This fissure produced a braided ʻaʻā flow that traversed the steep north flank at >1 km/hr, before slowing as slopes shallowed. Effusion rates peaked in the early hours of the eruption and then steadied at around 100 m3/s at fissure 3. Lava effusion became variable and dropped abruptly on December 7. The eruption ended on December 10 with a final flow length of 18 km. We tracked lava effusion rates, flow advance, channel dynamics, and 3-D flow field development with field and remote sensing observations throughout this eruption to inform lava flow forecasts and investigate ʻaʻā flow dynamics. Flows were mapped using optical and infrared imagery, uncrewed aircraft systems, and airborne single-pass InSAR. Video- and time-lapse cameras were deployed to record lava fountaining and channelized flow at frequent intervals, supplemented by field measurements. Analyses of these data built a detailed record of effusion rates and associated flow dynamics. We combine these high frequency observations of the proximal flow field with broader thermal and 3-D mapping of the evolution in the medial and distal flow field to examine ʻaʻā flow dynamics from source to toe. Within the steep, medial section of the flow field, channel blockages caused back-ups and overflows as the eruption progressed. Disrupted lava supply to the distal flow field and low slopes resulted in stalling, widening, and thickening at the toe. Overall, 142 million m3 of lava covered 36 km2 over 13 days. Our results inform the controls on development and evolution of lava flow fields, with implications for forecasting primary flow advance and secondary flow field growth.

Hawaii↗

USGS2023

Lessons learned from the 2022 CONVERSE Monogenetic Volcanism Response Scenario exercise

When volcanic unrest occurs, the scientific community can advance fundamental understanding of volcanic systems, but only with coordination before, during, and after the event across academic and governmental agencies. To develop a coordinated response plan, the Community Network for Volcanic Eruption Response (CONVERSE) orchestrated a scenario exercise centered around a hypothetical volcanic crisis in Arizona’s San Francisco Volcanic Field (SFVF). The exercise ran virtually from February 4 to March 4, 2022. Over 60 scientists from both academic and governmental spheres participated. The scenario exercise was assessed for its effectiveness in supporting collaborative production of knowledge, catalyzing transdisciplinary collaboration, supporting researcher confidence, and fostering a culture of inclusion within the volcanology community. This identified a need to support early career researchers through community and allyship. Overall, the 2022 CONVERSE exercise demonstrated how a fully remote, extended scenario can be authentically implemented and help broaden participation within the volcano science community.

Arizona↗

USGS2006

Carbon dioxide emissions from vegetation-kill zones around the resurgent dome of Long Valley caldera, eastern California, USA

A survey of diffuse CO 2 efflux, soil temperature and soil-gas chemistry over areas of localized vegetation-kill on and around the resurgent dome of Long Valley caldera California was performed to evaluate the premise that gaseous and thermal anomalies are related to renewed intrusion of magma. Some kill sites are long-lived features and others have developed in the past few years. Total anomalous CO 2 emissions from the thirteen areas average around 8.7 t per day; but the majority of the emissions come from four sites west of the Casa Diablo geothermal power plant. Geochemical analyses of the soil-gases from locations west and east of the plant revealed the presence of isobutane related to plant operations. The δ 13 C values of diffuse CO 2 range from − 5.7‰ to − 3.4‰, similar to values previously reported for CO 2 from hot springs and thermal wells around Long Valley. At many of the vegetation-kill sites soil temperatures reach boiling at depths ≤ 20 cm. Soil temperature/depth profiles at two of the high-emissions areas indicate that the conductive thermal gradient in the center of the areas is around 320 °C m − 1 . We estimate total heat loss from the two areas to be about 6.1 and 2.3 MW. Given current thinking on the rate of hydrothermal fluid flow across the caldera and using the CO 2 concentration in the thermal fluids, the heat and CO 2 loss from the kill areas is easily provided by the shallow hydrothermal system, which is sourced to the west of the resurgent dome. We find no evidence that the development of new areas of vegetation kill across the resurgent dome are related to new input of magma or magmatic fluids from beneath the resurgent dome. Our findings indicate that the areas have developed as a response to changes in the shallow hydrologic system. Some of the changes are likely related to fluid production at the power plant, but at distal sites the changes are more likely related to seismicity and uplift of the dome.

California↗

USGS2002

Structure and physical characteristics of pumice from the climactic eruption of Mount Mazama (Crater Lake), Oregon

The vesicularity, permeability, and structure of pumice clasts provide insight into conditions of vesiculation and fragmentation during Plinian fall and pyroclastic flow-producing phases of the ~7,700 cal. year B.P. climactic eruption of Mount Mazama (Crater Lake), Oregon. We show that bulk properties (vesicularity and permeability) can be correlated with internal textures and that the clast structure can be related to inferred changes in eruption conditions. The vesicularity of all pumice clasts is 75–88%, with >90% interconnected pore volume. However, pumice clasts from the Plinian fall deposits exhibit a wider vesicularity range and higher volume percentage of interconnected vesicles than do clasts from pyroclastic-flow deposits. Pumice permeabilities also differ between the two clast types, with pumice from the fall deposit having higher minimum permeabilities (~5×10 –13 m 2 ) and a narrower permeability range (5–50×10 –13 m 2 ) than clasts from pyroclastic-flow deposits (0.2–330×10 –13 m 2 ). The observed permeability can be modeled to estimate average vesicle aperture radii of 1–5 µm for the fall deposit clasts and 0.25–1 µm for clasts from the pyroclastic flows. High vesicle number densities (~10 9 cm –3 ) in all clasts suggest that bubble nucleation occurred rapidly and at high supersaturations. Post-nucleation modifications to bubble populations include both bubble growth and coalescence. A single stage of bubble nucleation and growth can account for 35–60% of the vesicle population in clasts from the fall deposits, and 65–80% in pumice from pyroclastic flows. Large vesicles form a separate population which defines a power law distribution with fractal dimension D=3.3 (range 3.0–3.5). The large D value, coupled with textural evidence, suggests that the large vesicles formed primarily by coalescence. When viewed together, the bulk properties (vesicularity, permeability) and textural characteristics of all clasts indicate rapid bubble nucleation followed by bubble growth, coalescence and permeability development. This sequence of events is best explained by nucleation in response to a downward-propagating decompression wave, followed by rapid bubble growth and coalescence prior to magma disruption by fragmentation. The heterogeneity of vesicle sizes and shapes, and the absence of differential expansion across individual clasts, suggest that post-fragmentation expansion played a limited role in the development of pumice structure. The higher vesicle number densities and lower permeabilities of pyroclastic-flow clasts indicate limited coalescence and suggest that fragmentation occurred shortly after decompression. Either increased eruption velocities or increased depth of fragmentation accompanying caldera collapse could explain compression of the pre-fragmentation vesiculation interval.

Oregon↗

USGS1995

Graphical and PC-software analysis of volcano eruption precursors according to the Materials Failure Forecast Method (FFM)

The Materials Failure Forecasting Method for volcanic eruptions (FFM) analyses the rate of precursory phenomena. Time of eruption onset is derived from the time of “failure” implied by accelerating rate of deformation. The approach attempts to fit data, Ω, to the differential relationship Ω Ω ¨ = A Ω ˙ "> ¨ = AΩ˙ , where the dot superscript represents the time derivative, and the data Ω may be any of several parameters describing the accelerating deformation or energy release of the volcanic system. Rate coefficients, A and α, may be derived from appropriate data sets to provide an estimate of time to “failure”. As the method is still an experimental technique, it should be used with appropriate judgment during times of volcanic crisis. Limitations of the approach are identified and discussed. Several kinds of eruption precursory phenomena, all simulating accelerating creep during the mechanical deformation of the system, can be used with FFM. Among these are tilt data, slope-distance measurements, crater fault movements and seismicity. The use of seismic coda, seismic amplitude-derived energy release and time-integrated amplitudes or coda lengths are examined. Usage of cumulative coda length directly has some practical advantages to more rigorously derived parameters, and RSAM and SSAM technologies appear to be well suited to real-time applications. One graphical and four numerical techniques of applying FFM are discussed. The graphical technique is based on an inverse representation of rate versus time. For α = 2, the inverse rate plot is linear; it is concave upward for α < 2 and concave downward for α > 2. The eruption time is found by simple extrapolation of the data set toward the time axis. Three numerical techniques are based on linear least-squares fits to linearized data sets. The “linearized least-squares technique” is most robust and is expected to be the most practical numerical technique. This technique is based on an iterative linearization of the given rate-time series. The hindsight technique is disadvantaged by a bias favouring a too early eruption time in foresight applications. The “log rate versus log acceleration technique”, utilizing a logarithmic representation of the fundamental differential equation, is disadvantaged by large data scatter after interpolation of accelerations. One further numerical technique, a nonlinear least-squares fit to rate data, requires special and more complex software. PC-oriented computer codes were developed for data manipulation, application of the three linearizing numerical methods, and curve fitting. Separate software is required for graphing purposes. All three linearizing techniques facilitate an eruption window based on a data envelope according to the linear least-squares fit, at a specific level of confidence, and an estimated rate at time of failure.

Journal of Volcanology and Geothermal Research↗

USGS2008

Characteristics, extent and origin of hydrothermal alteration at Mount Rainier Volcano, Cascades Arc, USA: Implications for debris-flow hazards and mineral deposits

Hydrothermal alteration at Mount Rainier waxed and waned over the 500,000-year episodic growth of the edifice. Hydrothermal minerals and their stable-isotope compositions in samples collected from outcrop and as clasts from Holocene debris-flow deposits identify three distinct hypogene argillic/advanced argillic hydrothermal environments: magmatic-hydrothermal, steam-heated, and magmatic steam (fumarolic), with minor superimposed supergene alteration. The 3.8 km 3 Osceola Mudflow (5600 y BP) and coeval phreatomagmatic F tephra contain the highest temperature and most deeply formed hydrothermal minerals. Relatively deeply formed magmatic-hydrothermal alteration minerals and associations in clasts include quartz (residual silica), quartz–alunite, quartz–topaz, quartz–pyrophyllite, quartz–dickite/kaolinite, and quartz–illite (all with pyrite). Clasts of smectite–pyrite and steam-heated opal–alunite–kaolinite are also common in the Osceola Mudflow. In contrast, the Paradise lahar, formed by collapse of the summit or near-summit of the edifice at about the same time, contains only smectite–pyrite and near-surface steam-heated and fumarolic alteration minerals. Younger debris-flow deposits on the west side of the volcano (Round Pass and distal Electron Mudflows) contain only low-temperature smectite–pyrite assemblages, whereas the proximal Electron Mudflow and a < 100 y BP rock avalanche on Tahoma Glacier also contain magmatic-hydrothermal alteration minerals that are exposed in the avalanche headwall of Sunset Amphitheater, reflecting progressive incision into deeper near-conduit alteration products that formed at higher temperatures. The pre-Osceola Mudflow alteration geometry is inferred to have consisted of a narrow feeder zone of intense magmatic-hydrothermal alteration limited to near the conduit of the volcano, which graded outward to more widely distributed, but weak, smectite–pyrite alteration within 1 km of the edifice axis, developed chiefly in porous breccias. The edifice was capped by a steam-heated alteration zone, most of which resulted from condensation of fumarolic vapor and oxidation of H 2 S in the unsaturated zone above the water table. Weakly developed smectite–pyrite alteration extended into the west and east flanks of the edifice, spatially associated with dikes that are localized in those sectors; other edifice flanks lack dikes and associated alteration. The Osceola collapse removed most of the altered core and upper east flank of the volcano, but intensely altered rocks remain on the uppermost west flank. Major conclusions of this study are that: (1) Hydrothermal–mineral assemblages and distributions at Mount Rainier can be understood in the framework of hydrothermal processes and environments developed from studies of ore deposits formed in analogous settings. (2) Frequent eruptions supplied sufficient hot magmatic fluid to alter the upper interior of the volcano hydrothermally, despite the consistently deep (≥ 8 km) magma reservoir which may have precluded formation of economic mineral deposits within or at shallow depths beneath Mount Rainier. The absence of indicator equilibrium alteration-mineral assemblages in the debris flows that effectively expose the volcano to a depth of 1–1.5 km also suggests a low potential for significant high-sulfidation epithermal or porphyry-type mineral deposits at depth. (3) Despite the long and complex history of the volcano, intensely altered collapse-prone rocks were spatially restricted to near the volcano's conduit system and summit, and short distances onto the upper east and west flanks, due to the necessary supply of reactive components carried by ascending magmatic fluids. (4) Intensely altered rocks were removed from the summit, east flank, and edifice interior by the Osceola collapse, but remain on the upper west flank in the Sunset Amphitheater area and present a continuing collapse hazard. (5) Visually conspicuous rocks on the lower east and mid-to-lower west flanks are not intensely altered and probably have not significantly weakened the rock, and thus do not present significant collapse hazards. (6) Alteration developed most intensely within breccia units, because of their high permeability and porosity. Volcanoes with abundant near-conduit upper-edifice breccias are prone to alteration increasing the possibility of collapse, whereas those that are breccia-poor (e.g., massive domes) are less prone to alteration.

Journal of Volcanology and Geothermal Research↗

USGS1989

The campi flegrei (Italy) geothermal system: A fluid inclusion study of the mofete and San Vito fields

A fluid inclusion study of core from the Mofete 1, Mofete 2, Mofete 5, San Vito 1, and San Vito 3 geothermal wells (Campi Flegrei, Campania, Italy) indicates that the hydrothermal minerals were precipitated from aqueous fluids (??CO2) that were moderately saline (3-4 wt.% NaCl equiv.) to hypersaline (> 26 wt.% NaCl equiv.) and at least in part, boiling. Three types of primary fluid inclusions were found in authigenic K-feldspar, quartz, calcite, and epidote: (A) two-phase [liquid (L) + vapor (V)], liquid-rich inclusions with a range of salinity; (B) two-phase (L + V), vaporrich inclusions with low salinity; and (C) three-phase [L + V + crystals (NaCL)], liquid-rich inclusions with hypersalinity. Results of microthermometric and crushing studies are reported for twenty drill core samples taken from the lower portions of the five vertical wells. Data presented for selected core samples reveal a general decrease in porosity and increase in bulk density with increasing depth and temperature. Hydrothermal minerals commonly fill fractures and pore-spaces and define a zonation pattern, similar in all five wells studied, in response to increasing depth (pressure) and temperature. A greenschist facies assemblage, defined by albite + actinolite, gives way to an amphibolite facies, defined by plagioclase (andesine) + hornblende, in the San Vito 1 well at about 380??C. The fluid inclusion salinity values mimic the saline and hypersaline fluids found by drilling. Fluid inclusion V/L homogenization temperatures increase with depth and generally correspond to the extrapolated down-hole temperatures. However, fluid inclusion data for Mofete 5 and mineral assemblage data for San Vito 3, indicate fossil, higher-temperature regimes. A limited 87Sr/86Sr study of leachate (carbonate) and the leached cores shows that for most samples (except San Vito 3) the carbonate deposition has been from slightly 87Sr-enriched fluids and that Sr isotopic exchange has been incomplete. However, San Vito 3 cores show an approach to fluid/rock Sr equilibrium with a fluid similar to modern ocean water in 87Sr/86Sr ratio. The Campi Flegrei volcanic system has evolved undersaturated products, mostly trachyte, and defines a large (??? 12 km) caldera. The hydrothermal system developed in this location can be used as an analog for fossil systems in similar trachytic environments. The potential for ore mineralization is expressed by the recognition, from fluid inclusion and drilling data, of ore-forming environments such as boiling and brine stratification. ?? 1989.

Journal of Volcanology and Geothermal Research↗

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

Understanding the drivers of volcano deformation through geodetic model verification and validation

Volcano geodesy often involves the use of models to explain observed surface deformation. A variety of forward models are used, from analytical point sources to numerical simulations that consider complex magma system geometries, topography, and material properties. Various inversion methods can then be used to relate observed volcano data to models. Ideally, forward models should be verified through intercomparison, to check for implementation errors and quantify the error induced by any approximations used. Additionally, forward models and inversion methods should be validated through tests with synthetic and/or real data, to determine their ability to match data and estimate parameter values within uncertainty. However, to date, there have not been comprehensive verification and validation efforts in volcano geodesy. Here, we report on the first phase of the Drivers of Volcano Deformation (DVD) exercises, which were designed to build community involvement through web-based exercises involving calculations of static elastic displacement around pressurized magma reservoirs. The forward model exercises begin with a spherical reservoir in a homogeneous half space, then introduce topography, heterogeneous elastic properties, and spheroidal geometries. The inversion exercises provide synthetic noisy surface displacement data for a spherical reservoir in a homogeneous half space and assess consistency in estimates of reservoir location and volume/pressure change. There is variability in the results from both forward modeling and inversions, which highlights the strengths and limitations of different forward models, as well as the importance of inversion method choice and uncertainty quantification. This first phase of the DVD exercises serves as a community resource and will facilitate further efforts to develop standards of reproducibility.

Bulletin of Volcanology↗

USGS2021

Repeating earthquakes during multiple phases of unrest and eruption at Mount Agung, Bali, Indonesia, 2017

In 2017, Mount Agung produced a small (VEI 2) eruption that was preceded by an energetic volcano-tectonic (VT) swarm (>800 earthquakes per day up to M4.9) and two months of declining activity. The period of decreased seismic activity complicated forecasting efforts for scientists monitoring the volcano. We examine the time history of earthquake families at Mount Agung in search of additional insight into the temporal changes in the shallow crust prior to eruption. Specifically, we analyze the period of declining seismic activity about five weeks prior to the eruption when forecasting uncertainty was greatest. We use REDPy (Hotovec-Ellis and Jeffries, 2016) to build a catalog of 6,508 earthquakes from 18 October 2017–15 February 2018 and group them into families of repeating earthquakes based on waveform similarity using a cross-correlation coefficient threshold of 0.8. We show that the evolution of earthquake families provides evidence that Mount Agung was progressing toward eruption even though overall earthquake rates and seismic-energy-release declined. We find that earthquake families that dominated seismicity during the beginning of the crisis ceased near the onset of tremor on 12 November 2017. Then, earthquake families took on characteristics commonly observed during effusive phases of eruptions on 15 November—a full six days before the first phreatomagmatic eruption on 21 November 2017 and a full ten days before the actual onset of lava effusion on 25 November 2017. We interpret the transitions in seismicity as the manifestation of a three-phase physical model including an Intrusion Phase, a Transition Phase, and a Eruptive Phase. During the Intrusion Phase, seismicity was dominated by VT earthquakes with a relatively high percentage of repeaters (59%) grouped into numerous (65) simultaneous families. During the Eruptive Phase, seismicity included both VT and low frequency earthquakes that grouped into relatively long-lived families despite a low overall percentage of repeaters (14%). The Transition Phase exhibited characteristics of earthquake families between the Intrusion Phase and Eruptive Phase. We conclude that the time history of earthquake families provides insight into the evolution of the stress distribution in the volcanic edifice, the development of the volcanic conduit, and seismogenesis of magma effusion. Finally, we discuss the role that repeating earthquakes could play in real-time monitoring at restless volcanoes. Our work suggests eruption forecasts can be improved by incorporating automatic processing codes to assist seismologists during sustained periods of high earthquake rates, even at sparsely monitored volcanoes.

Bali, Mount Agung↗

USGS1990

High-precision 40Ar/39Ar sanidine geochronology of ignimbrites in the Mogollon-Datil volcanic field, southwestern New Mexico

40Ar/39Ar age spectra have been obtained from 85 sanidine separates from 36 ignimbrites and one rhyolitic lava in the latest Eocene-Oligocene Mogollon-Datil volcanic field of southwestern New Mexico. Of the 97 measured age spectra, 94 yield weighted-mean plateau ages each giving single-spectrum 1?? precision of??0.25%-0.4% (??0.07-0.14 Ma). Replicate plateau age determinations for eight different samples show within-sample 1?? precisions averaging ??0.25%. Plateau ages from multiple (n=3-8) samples of individual ignimbrites show 1?? within-unit precision of ??0.1%-0.4% (??0.04-0.13 Ma). This within-unit precision represents a several-fold improvement over published K-Ar data for the same ignimbrites, and is similar to the range of precisions reported from single-crystal laser fusion studies. A further indication of the high precision of unit-mean 40Ar/30Ar ages is their close agreement with independently established stratigraphic order. Two samples failed to meet plateau criteria, apparently due to geologic contamination by older feldspars. Effects of minor contamination are shown by six other samples, which yielded slightly anomalous plateau ages. 40Ar/39Ar plateau ages permit resolution of units differing in age by 0.5% (0.15 Ma) or less. This high resolution, combined with paleomagnetic studies, has helped to correlate ignimbrites among isolated ranges and has allowed development of an integrated timestratigraphic framework for the volcanic field. Mogollon-Datil ignimbrites range in age from 36.2 to 24.3 Ma. Ignimbrite activity was strongly episodic, being confined to four brief (<2.6 m.y.) eruptive episodes separated by 1-3 m.y. gaps. Ignimbrite activity generally tended to migrate from the southeast toward the north and west. ?? 1990 Springer-Verlag.

Bulletin of Volcanology↗

USGS2020

Evolution of the submarine–subaerial edifice of Bogoslof volcano, Alaska, during its 2016–2017 eruption based on analysis of satellite imagery

The 2016–2017 eruption of Bogoslof volcano involved at least 70 detected eruptive events between mid-December 2016 and August 30, 2017. Acquisition of high-resolution satellite imagery throughout the duration of the eruptive period allowed us to document and map the various morphologic changes that occurred on the subaerial part of Bogoslof Island. The emplacement of pyroclastic-flow and surge deposits caused the island to increase in area by about 1.5 km 2 . The dominant volcanic landforms of the eruption were a series of tuff rings emplaced around various submarine vents. Many of the tuff rings were mantled with surface dunes and impressive amounts of ballistic ejecta, likely derived from erupting magma bodies or previously emplaced submarine lava domes. Debris-flow deposits and surface channels extending over tuff ring surfaces apparent in multiple satellite images are evidence for explosive ejection of seawater. In most cases, erupting vents were initially submarine or began at subaerial lava domes and were largely flooded by seawater suggesting that water-magma ratios were likely high. Under such conditions where water is abundant, eruptive products typically reflect a high degree of water involvement and are dominated by the formation of wet tephra jets and flows and associated deposits typically consist of fine ash and lapilli, contain accretionary lapilli and ash aggregates, and usually form tuff cones and mounds. We observed none of these features in our analysis of satellite data or during our examination of eruptive deposits on Bogoslof Island in 2018. On the contrary, the dominant landform associated with the Bogoslof eruption was tuff rings. The development of tuff rings and surface dunes are commonly associated with the formation of pyroclastic base surges that are by comparison emplaced relatively dry. Dry base surge deposits can be generated from phreatomagmatic explosions involving superheated steam. It is possible that shallow submarine, magma–wet sediment interactions were a characteristic and possibly a dominant eruptive process of the 2016–2017 Bogoslof eruption.

Alaska↗

USGS2002

The oligocene Lund Tuff, Great Basin, USA: A very large volume monotonous intermediate

Unusual monotonous intermediate ignimbrites consist of phenocryst-rich dacite that occurs as very large volume (> 1000 km3) deposits that lack systematic compositional zonation, comagmatic rhyolite precursors, and underlying plinian beds. They are distinct from countless, usually smaller volume, zoned rhyolite-dacite-andesite deposits that are conventionally believed to have erupted from magma chambers in which thermal and compositional gradients were established because of sidewall crystallization and associated convective fractionation. Despite their great volume, or because of it, monotonous intermediates have received little attention. Documentation of the stratigraphy, composition, and geologic setting of the Lund Tuff - one of four monotonous intermediate tuffs in the middle-Tertiary Great Basin ignimbrite province - provides insight into its unusual origin and, by implication, the origin of other similar monotonous intermediates. The Lund Tuff is a single cooling unit with normal magnetic polarity whose volume likely exceeded 3000 km3. It was emplaced 29.02 ?? 0.04 Ma in and around the coeval White Rock caldera which has an unextended north-south diameter of about 50 km. The tuff is monotonous in that its phenocryst assemblage is virtually uniform throughout the deposit: plagioclase > quartz ??? hornblende > biotite > Fe-Ti oxides ??? sanidine > titanite, zircon, and apatite. However, ratios of phenocrysts vary by as much as an order of magnitude in a manner consistent with progressive crystallization in the pre-eruption chamber. A significant range in whole-rock chemical composition (e.g., 63-71 wt% SiO2) is poorly correlated with phenocryst abundance. These compositional attributes cannot have been caused wholly by winnowing of glass from phenocrysts during eruption, as has been suggested for the monotonous intermediate Fish Canyon Tuff. Pumice fragments are also crystal-rich, and chemically and mineralogically indistinguishable from bulk tuff. We postulate that convective mixing in a sill-like magma chamber precluded development of a zoned chamber with a rhyolitic top or of a zoned pyroclastic deposit. Chemical variations in the Lund Tuff are consistent with equilibrium crystallization of a parental dacitic magma followed by eruptive mixing of compositionally diverse crystals and high-silica rhyolite vitroclasts during evacuation and emplacement. This model contrasts with the more systematic withdrawal from a bottle-shaped chamber in which sidewall crystallization creates a marked vertical compositional gradient and a substantial volume of capping-evolved rhyolite magma. Eruption at exceptionally high discharge rates precluded development of an underlying plinian deposit. The generation of the monotonous intermediate Lund magma and others like it in the middle Tertiary of the western USA reflects an unusually high flux of mantle-derived mafic magma into unusually thick and warm crust above a subducting slab of oceanic lithosphere. ?? 2002 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

USGS2022

Volcano geodesy using InSAR in 2020: The past and next decades

The study of volcano deformation has grown significantly through they year 2020 since the development of interferometric synthetic aperture radar (InSAR) in the 1990s. This relatively new data source, which provides evidence of changes in subsurface magma storage and pressure without the need for ground-based equipment, has matured during the past decade. It now provides a means to address previously inaccessible questions and offers input to increasingly complex models of magmatic processes. Here, we review how technological advances in InSAR during 2010-2020 have facilitated our ability to monitor and interpret volcanic processes, primarily through rapid and accurate observations of the changing surfaces at active volcanoes worldwide. Specifically, we examine how current systems achieve excellent resolution in time and space, provide global coverage, and generate products that are easy to use by non-specialists—factors that have often limited the practical study of volcanoes using radar measurements. We also look to the future, offering our perspective about how advancements in technology and data management in the decade to come will increase the value and accessibility of InSAR applied to the geodetic study of volcanoes and monitoring of hazardous volcanic processes. New developments will include the launch of additional satellites by both public space agencies and private companies, as well as implementation of algorithms for exploiting the growing volumes of data. To meet their full potential, these efforts will require coordination between data users and data providers so that the relevant imagery is acquired, made available to volcanologists in a timely fashion, and utilized to assess and mitigate volcanic hazards.

Bulletin of Volcanology↗

USGS2016

Monogenetic volcanoes fed by interconnected dikes and sills in the Hopi Buttes volcanic field, Navajo Nation, USA

Although monogenetic volcanic fields pose hazards to major cities worldwide, their shallow magma feeders (<500 m depth) are rarely exposed and, therefore, poorly understood. Here, we investigate exposures of dikes and sills in the Hopi Buttes volcanic field, Arizona, to shed light on the nature of its magma feeder system. Shallow exposures reveal a transition zone between intrusion and eruption within 350 m of the syn-eruptive surface. Using a combination of field- and satellite-based observations, we have identified three types of shallow magma systems: (1) dike-dominated, (2) sill-dominated, and (3) interconnected dike-sill networks. Analysis of vent alignments using the pyroclastic massifs and other eruptive centers (e.g., maar-diatremes) shows a NW-SE trend, parallel to that of dikes in the region. We therefore infer that dikes fed many of the eruptions. Dikes are also observed in places transforming to transgressive (ramping) sills. Estimates of the observable volume of dikes (maximum volume of 1.90 × 10 6 m 3 ) and sills (minimum volume of 8.47 × 10 5 m 3 ) in this study reveal that sills at Hopi Buttes make up at least 30 % of the shallow intruded volume (∼2.75 × 10 6 m 3 total) within 350 m of the paeosurface. We have also identified saucer-shaped sills, which are not traditionally associated with monogenetic volcanic fields. Our study demonstrates that shallow feeders in monogenetic fields can form geometrically complex networks, particularly those intruding poorly consolidated sedimentary rocks. We conclude that the Hopi Buttes eruptions were primarily fed by NW-SE-striking dikes. However, saucer-shaped sills also played an important role in modulating eruptions by transporting magma toward and away from eruptive conduits. Sill development could have been accompanied by surface uplifts on the order of decimeters. We infer that the characteristic feeder systems described here for the Hopi Buttes may underlie monogenetic fields elsewhere, particularly where magma intersects shallow, and often weak, sedimentary rocks. Results from this study support growing evidence of the important role of shallow sills in active monogenetic fields.

Arizona↗

USGS2024

Following the tug of the audience from complex to simplified hazards maps at Cascade Range volcanoes

Volcano-hazard maps are broadly recognized as important tools for forecasting and managing volcanic crises and for disseminating spatial information to authorities and people at risk. As scientists, we might presume that hazards maps can be developed at the time and with the methods of our discretion, yet the co-production of maps with stakeholder groups, who have programmatic needs of their own, can sway the timing, usability, and acceptance of map products. We examine two volcano hazard map-making efforts by staff at the U.S. Geological Survey. During the 1990s and early 2000s scientists developed a series of hazard assessments and maps with detailed zonations for volcanoes in Washington and Oregon. In 2009, the National Park Service expressed the need for simplified versions of the existing hazard maps for a high-profile visitor center exhibit. This request created an opportunity for scientists to rethink the objectives, scope, content, and map representations of hazards. The primary focus of this article is a discussion of processes used by scientists to distill the most critical information within the official parent maps into a series of simplified maps using criteria specified. We contextualize this project with information about development of the parent maps, public response to the simplified hazard maps, the value of user engagement in mapmaking, and with reference to the abundance of guidance available to the next generation of hazard-mapmakers. We argue that simplified versions of maps should be developed in tandem with any hazard maps that contain technical complexities, not as a replacement, but as a mechanism to broaden awareness of hazards. We found that when scientists endeavor to design vivid and easy-to-understand maps, people in many professions find uses for them within their organization’s information products, resulting in extensive distribution.

Washington↗

USGS2014

Pyroclast/snow interactions and thermally driven slurry formation. Part 2: Experiments and theoretical extension to polydisperse tephra

Erosion of snow by pyroclastic flows and surges presumably involves mechanical scour, but there may be thermally driven phenomena involved as well. To investigate this possibility, layers of hot (up to 400??C), uniformly sized, fine- to medium-grained sand were emplaced vertically onto finely shaved ice ('snow'); thus there was no relative shear motion between sand and snow and no purely mechanical scour. In some cases large vapor bubbles, commonly more than 10 mm across, rose through the sand layer, burst at the surface, and caused complete convective overturn of the sand, which then scoured and mixed with snow and transformed into a slurry. In other cases no bubbling occurred and the sand passively melted its way downward into the snow as a wetting front moved upward into the sand. A continuum of behaviors between these two cases was observed. Vigorous bubbling and convection were generally favored by high temperature, small grain size, and small layer thickness. A physically based theory of heat- and mass transfer at the pyroclast/snow interface, developed in Part 1 of this paper, does a good job of explaining the observations as a manifestation of unstable vapor-driven fluidization. The theory, when extrapolated to the behavior of actual, poorly sorted pyroclastic flow sediments, leads to the prediction that the observed 'thermal-scour' phenomenon should also occur for many real pyroclastic flows passing over snow. 'Thermal scour' is therefore likely to be involved in the generation of lahars.

Bulletin of Volcanology↗

USGS2020

Indonesia and the United States team up to reduce impacts from dangerous volcanoes

With 75 historically active volcanoes, Indonesia is the world’s most volcanically active nation. Its volcanoes are legendary throughout the world, with the notorious 19th-century eruptions at Mount Tambora (1815) and Krakatau (1883), and the eruption that created the giant Toba Caldera in Sumatra (75,000 years ago)—the Earth’s largest volcanic eruption in the past 100,000 years. Just in the past 20 years, more than 36 volcanoes have erupted, some multiple times, threatening millions of people. The Government of Indonesia responds to this daunting volcanic threat from Indonesia’s Center for Volcanology and Geologic Hazards Mitigation (CVGHM) offices in Bandung. From there, dozens of staff members coordinate operations at 77 small observatories spread out across the 5,000-kilometer-long archipelago that is Indonesia. For more than 50 years, the United States has worked in partnership with the Government of Indonesia to reduce the risk of volcanic eruptions and to save lives. Today, the two governments cooperate through an agreement between CVGHM and the U.S. Geological Survey’s Volcano Disaster Assistance Program (VDAP). Each year CVGHM and VDAP representatives convene to discuss priorities for the upcoming 12 months. They outline a series of training efforts, field campaigns, equipment deployments, and timely topics that require further exploration and development.

Fact Sheet↗