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 235 records · Page 13Linked to original sources
USGS2016

Shallow and deep controls on lava lake surface motion at Kīlauea Volcano

Lava lakes provide a rare window into magmatic behavior, and lake surface motion has been used to infer deeper properties of the magmatic system. At Halema'uma'u Crater, at the summit of Kīlauea Volcano, multidisciplinary observations for the past several years indicate that lava lake surface motion can be broadly divided into two regimes: 1) stable and 2) unstable. Stable behavior is driven by lava upwelling from deeper in the lake (presumably directly from the conduit) and is an intrinsic process that drives lava lake surface motion most of the time. This stable behavior can be interrupted by periods of unstable flow (often reversals) driven by spattering – a shallowly-rooted process often extrinsically triggered by small rockfalls from the crater wall. The bursting bubbles at spatter sources create void spaces and a localized surface depression which draws and consumes surrounding surface crust. Spattering is therefore a location of lava downwelling, not upwelling. Stable (i.e. deep, upwelling-driven) and unstable (i.e. shallow, spattering-driven) behavior often alternate through time, have characteristic surface velocities, flow directions and surface temperature regimes, and also correspond to changes in spattering intensity, outgassing rates, lava level and seismic tremor. These results highlight that several processes, originating at different depths, can control the motion of the lava lake surface, and long-term interdisciplinary monitoring is required to separate these influences. These observations indicate that lake surface motion is not always a reliable proxy for deeper lake or magmatic processes. From these observations, we suggest that shallow outgassing (spattering), not lake convection, drives the variations in lake motion reported at Erta 'Ale lava lake.

Hawaii↗
USGS2017

Olivine-melt relationships and syneruptive redox variations in the 1959 eruption of Kīlauea Volcano as revealed by XANES

The 1959 summit eruption of Kīlauea Volcano exhibited high lava fountains of gas-rich, primitive magma, containing olivine + chromian spinel in highly vesicular brown glass. Microprobe analysis of these samples shows that euhedral rims on olivine phenocrysts, in direct contact with glass, vary significantly in forsterite (Fo) content, at constant major-element melt composition, as do unzoned groundmass olivine crystals. Ferric/total iron (Fe + 3 /Fe T )ratios for matrix and interstitial glasses, plus olivine-hosted glass inclusions in eight 1959 scoria samples have been determined by micro X-ray absorption near-edge structure spectroscopy (μ-XANES). These data show that much of the variation in Fo content reflects variation in oxidation state of iron in the melt, which varies with sulfur concentration in the glass and (locally) with proximity to scoria edges in contact with air. Data for 24 olivine-melt pairs in the better-equilibrated samples from later in the eruption show K D averaging 0.280 ± 0.03 for the exchange of Fe and Mg between olivine and melt, somewhat displaced from the value of 0.30 ± 0.03 given by Roeder and Emslie (1970). This may reflect the low SiO 2 content of the 1959 magmas, which is lower than that in most Kīlauea tholeiites. More broadly, we show the potential of μ-XANES and electron microprobe to revisit and refine the value of K D in natural systems. The observed variations of Fe + 3 /Fe T ratios in the glasses reflect two distinct processes. The main process, sulfur degassing, produces steady decrease of the Fe + 3 /Fe T ratio. Melt inclusions in olivine are high in sulfur (1060–1500 ppm S), with Fe + 3 /Fe T = 0.160–0.175. Matrix glasses are degassed (mostly S < 200 ppm) with generally lower Fe + 3 /Fe T (0.114–0.135). Interstitial glasses within clumps of olivine crystals locally show intermediate levels of sulfur and Fe + 3 /Fe T ratio. The correlation suggests that (1) the 1959 magma was significantly reduced by sulfur degassing during the eruption and (2) the melts originally had Fe + 3 /Fe T ≥ 0.175, consistent with oxygen fugacity ( f O 2 ) at least 0.4 log units above the fayalite-magnetite-quartz (FMQ) buffer at 1 atm and magmatic temperature of 1200 °C. The second process is interaction between the melts and atmospheric oxygen, which results in higher Fe + 3 /Fe T ratios. Detailed μ-XANES traverses show gradients in Fe + 3 /Fe T of 0.145 to 0.628 over distances of 100–150 μm in thin, visibly reddened matrix glass bordering some scoriae, presumably caused by contact with air. This process was extremely rapid, giving insight into how fast the Fe + 3 /Fe T ratio can change in response to changes in external conditions.

Hawaii↗
USGS2016

Thermal mapping of a pahoehoe lava flow, Kilauea Volcano

Pāhoehoe lava flows are a major component of Hawaiian eruptive activity, and an important part of basaltic volcanism worldwide. In recent years, pāhoehoe lava has destroyed homes and threatened parts of Hawai‘i with inundation and disruption. In this study, we use oblique helicopter-borne thermal images to create high spatial resolution (~1 m) georeferenced thermal maps of the active pāhoehoe flow on Kīlauea Volcano’s East Rift Zone. Thermal maps were created on 27 days during 2014–2016 in the course of operational monitoring, encompassing a phase of activity that threatened the town of Pāhoa. Our results illustrate and reinforce how pāhoehoe flows are multicomponent systems consisting of the vent, master tube, distributary tubes and surface breakouts. The thermal maps accurately depict the distribution and character of pāhoehoe breakouts through time, and also delineate the subsurface lava tube. Surface breakouts were distributed widely across the pāhoehoe flow, with significant portions concurrently active well upslope of the flow front, often concentrated in clusters of activity that evolved through time. Gradual changes to surface breakout distribution and migration relate to intrinsic processes in the flow, including the slow evolution of the distributary tube system. Abrupt disruptions to this system, and the creation of new breakouts (and associated hazards), were triggered by extrinsic forcing—namely fluctuations in lava supply rate at the vent which disrupted the master lava tube. Although the total area of a pāhoehoe flow has been suggested to relate to effusion rate, our results show that changes in the proportion of expansion vs. overplating can complicate this relationship. By modifying existing techniques, we estimate time-averaged discharge rates for the flow during 2014–2016 generally in the range of 1–2 m3 s-1 (mean: 1.3±0.4 m3 s-1) – less than half of Kīlauea’s typical eruption rate on the East Rift Zone and suggestive of a weak eruptive regime during 2014–2016. We caution, however, that this discharge rate approach requires further independent corroboration. The thermal maps provide the first synoptic characterization of pāhoehoe flow activity at high spatial resolution, essential both for operational hazard assessment and fundamental understanding of pāhoehoe behavior.

Hawaii↗
USGS2017

Early postcaldera rhyolite and structural resurgence at Long Valley Caldera, California

After the 767-ka caldera-forming eruption of 650 km 3 of rhyolite magma as the Bishop Tuff, 90–100 km 3 of similar rhyolite erupted in the west-central part of Long Valley caldera in as many as 40 batches spread over the 110,000-year interval from ~ 750 ka to ~ 640 ka. Centrally, this Early Rhyolite (ER) is as thick as 622 m, but it spread radially to cover much of the caldera floor, where half its area is now concealed by post-ER sediments and lavas. At least 75% of the ER is aphyric rhyolite tuff. Drillholes encountered ~ 22 (altered) ER lava flows intercalated in the pyroclastic pile, and another 11 units of (largely fresh) ER lava are exposed on the caldera's resurgent dome and at Lookout Mountain. Exposed units have been distinguished, mapped, studied petrographically and chemically, and radioisotopically dated; each is described in detail. Their phenocryst contents range from 0 to 2.5 wt%. All the phyric units have plagioclase, orthopyroxene, and ilmenite; most have biotite and rare tiny magnetite, and a few contain rare zircon. The compositional range of fresh obsidians is narrow—74.3–75.0% SiO 2 , 1.21–1.37% FeO*, and 5.12–5.26% K 2 O, but wider variations in Ti, Ba, Sr, and Zr permit distinction of individual units and eruptive groups. The limited chemical and petrographic variability shown by so many ER batches released episodically for ~ 110,000 years suggests a thermally buffered and well-stirred reservoir. The ER central area, where ER eruptions had taken place, was uplifted ~ 400 m to form a structural dome ~ 10 km in diameter. Most of the inflation is attributable to 10 sills of ER that intrude the Bishop Tuff beneath the uplift, but other processes potentially contributing to resurgence are also considered. As shown by erratics of Mesozoic rocks ice-rafted from the Sierra Nevada and dropped on ER lavas, much of the ER had erupted early enough and at low enough elevation to be inundated by the intracaldera lake and was only later lifted by the resurgence that also raised clusters of the erratics hundreds of meters higher than any shoreline. Most of the uplift was over by ~ 570 ka, but dome-crossing faults that exhibit normal throw of 10–30 m cut lavas as young as 175–125 ka. For most elements, chemical ranges of the ER lie within those of the zoned Bishop Tuff, which had erupted earlier from the same place. Only Ba, Zr, Hf, and Eu/Eu* extend to ranges outside those of the Bishop Tuff, nominally to less evolved compositions. Initial 87 Sr/ 86 Sr values of ER are likewise within the range of the Bishop Tuff, but ER ratios of 143 Nd/ 144 Nd and 206 Pb/ 204 Pb extend beyond those of the Bishop Tuff to values slightly more influenced by upper-crustal contributions. FeTi-oxide geothermometry yields 752°–844 °C for ER, compared to 700°–820 °C for the Bishop Tuff. ER fO 2 values are 0.5–1.0 log units more reduced than those of the T–fO 2 array of the Bishop Tuff. The postcaldera reduction may reflect reaction with graphite from the black lithics of Paleozoic graphitic metapelite so abundant in the Bishop Tuff. Much of the pumice emplaced during the later half of the Bishop Tuff eruption has 10–25 wt% phenocrysts, dominantly quartz and sanidine, but the 100 km 3 of ER has only 0–2.5 wt% and completely lacks quartz and sanidine. Postcaldera processes, including mixing, volatile ascent, and crystal resorption, as well as potential contaminants and magmatic inputs, are all considered.

California↗
USGS2019

Character of community response to volcanic crises at Sinabung and Kelud volcanoes

Volcano disaster mitigation is a priority in Indonesia's national development plan due to the large number of active volcanoes, high frequency of eruptions and dense population within volcanic hazards zones. As a consequence of these factors and the many decades of experience in mitigation of volcanic risk, Indonesia demonstrates that effective community response is critical to avoid disasters. Understanding of hazards, risks and early-warning systems are the main factors in building an effective community response. The level of understanding of these factors is proportional to a community's independence in taking appropriate actions and is reflected in the community's preparedness. Further, the experience gained from responses to the eruptions of Sinabung and Kelud as well as from other Indonesian volcanoes shows that a local leader plays an important role in mobilizing a community. Such a leader utilizes knowledge of local culture (local wisdom) and an understanding of the community's character to encourage community members to participate and empower themselves. Sharing of information and the formulation of a disaster mitigation plan by government and community also results in a shared commitment to increase participation and community empowerment. Indonesia uses four volcano alert levels to activate community mitigation plans. In order of increasing criticality, these are: Normal, Advisory (Waspada), Watch (Siaga) and Warning (Awas). Implementation of policies, strategies and actions are tied to each level. In addition, mandatory actions by certain stakeholders are required at certain alert levels. These mandatory actions increase the effectiveness of disaster mitigation by both government and community. A comparison of the crisis responses of the Sinabung and Kelud communities is carried out here in order to better understand problems, learn lessons and improve the process of Indonesia's community mitigation policies. Local culture, the start time of community involvement (in normal versus crisis time), degree and level of government involvement, roles of local leaders, local perception of hazards, and political intervention are all factors that influence a community's response at the time of eruption. We stress that an understanding of the “ community character ,” which includes a number of cultural, social and knowledge parameters is essential for effective crisis management. The responses of Sinabung and Kelud illustrate both the key role of community in disaster mitigation and the synergy that can result from close coordination and collaboration between government and community.

Journal of Volcanology and Geothermal Research↗
USGS2017

Eruptive history of the Ubehebe Crater Cluster, Death Valley, California

A sequence of late Holocene eruptions from the Ubehebe Crater cluster in Death Valley was short-lived, emplacing several phreatomagmatic and magmatic deposits. Seven craters form the main group, which erupted along a north-south alignment 1.5 km long. At least five more make a 500-m east-west alignment west of the main crater group. One more is an isolated shallow crater ~ 400 m south of that alignment. All erupted through Miocene fanglomerate and sandstone, which are now distributed as comminuted matrix and lithic clasts in all Ubehebe deposits. Stratigraphic evidence showing that all Ubehebe strata were emplaced within a short time interval includes: (1) deposits from the many Ubehebe vents make a multi-package sequence that conformably drapes paleo-basement topography with no erosive gullying between emplacement units; (2) several crater rims that formed early in the eruptive sequence are draped smoothly by subsequent deposits; and (3) tack-welded to agglutinated spatter and bombs that erupted at various times through the sequence remained hot enough to oxidize the overlying youngest emplacement package. In addition, all deposits sufficiently consolidated to be drilled yield reliable paleomagnetic directions, with site mean directions showing no evidence of geomagnetic secular variation. Chemical analyses of juvenile components representing every eruptive package yield a narrow range in major elements [SiO 2 (48.65–50.11); MgO (4.98–6.23); K 2 O (2.24–2.39)] and trace elements [Rb (28–33); Sr (1513–1588); Zr (373–404)]. Despite lithologic similarities, individual fall units can be traced outward from vent by recording layer thicknesses, maximum scoria and lithic sizes, and juvenile clast textural variations. This permits reconstruction of the eruptive sequence, which produced a variety of eruptive styles. The largest and northernmost of the craters, Ubehebe Crater, is the youngest of the group. Its largely phreatomagmatic deposits drape all of the others, thicken in paleogullies and thin over several newly created crater rims. Evidence in-hand virtually requires that the Ubehebe cluster of craters erupted over a brief time interval, not protracted over centuries.

California↗
USGS2017

Magmatic degassing, lava dome extrusion, and explosions from Mount Cleveland volcano, Alaska, 2011–2015: Insight into the continuous nature of volcanic activity over multi-year timescales

Mount Cleveland volcano (1730 m) is one of the most active volcanoes in the Aleutian arc, Alaska, but heightened activity is rarely accompanied by geophysical signals, which makes interpretation of the activity difficult. In this study, we combine volcanic gas emissions measured for the first time in August 2015 with longer-term measurements of thermal output and lava extrusion rates between 2011 and 2015 calculated from MODIS satellite data with the aim to develop a better understanding of the nature of volcanic activity at Mount Cleveland. Degassing measurements were made in the month following two explosive events (21 July and 7 August 2015) and during a period of new dome growth in the summit crater. SO 2 emission rates ranged from 400 to 860 t d − 1 and CO 2 /SO 2 ratios were < 3, consistent with the presence of shallow magma in the conduit and the observed growth of a new lava dome. Thermal anomalies derived from MODIS data from 2011 to 2015 had an average repose time of only 4 days, pointing to the continuous nature of volcanic activity at this volcano. Rapid increases in the cumulative thermal output were often coincident with visual confirmation of dome growth or accumulations of tephra in the crater. The average rate of lava extrusion calculated for 9 periods of rapid increase in thermal output was 0.28 m 3 s − 1 , and the total volume extruded from 2011 to 2015 was 1.9–5.8 Mm 3 . The thermal output from the lava extrusion events only accounts for roughly half of the thermal budget, suggesting a continued presence of shallow magma in the upper conduit, likely driven by convection. Axisymmetric dome morphology and occasional drain back of lava into the conduit suggests low-viscosity magmas drive volcanism at Mount Cleveland. It follows also that only small overpressures can be maintained given the small domes and fluid magmas, which is consistent with the low explosivity of most of Mount Cleveland's eruptions. Changes between phases of dome growth and explosive activity are somewhat unpredictable and likely result from plugs that are related to the dome obtaining a critical dimension, or from small variations in the magma ascent rate that lead to crystallization-induced blockages in the upper conduit, thereby reducing the ability of magma to degas. We suggest the small magma volumes, slow ascent rates, and low magma viscosity lead to the overall lack of anomalous geophysical signals prior to eruptions, and that more continuous volcanic degassing measurements might lead to more successful eruption forecasting at this continuously-active open-vent volcano.

Alaska↗
USGS2017

Plateaus and sinuous ridges as the fingerprints of lava flow inflation in the Eastern Tharsis Plains of Mars

The Tharsis Montes rift aprons are composed of outpourings of lava from chaotic terrains to the northeast and southwest flank of each volcano. Sinuous and branching channel networks that are present on the rift aprons suggest the possibility of fluvial processes in their development, or erosion by rapidly emplaced lavas, but the style of lava flow emplacement throughout rift apron development is not clearly understood. To better characterize the style of lava emplacement and role of fluvial processes in rift apron development, we conducted morphological mapping of the Pavonis Mons southwest rift apron and the eastern Tharsis plains using images from the High Resolution Imaging Science Experiment (HiRISE), Mars Orbiter Camera (MOC), Context Camera (CTX), Thermal Emission Imaging System (THEMIS), and High Resolution Stereo Camera (HRSC) along with the Mars Orbiter Laser Altimeter (MOLA) Precision Experiment Data Records (PEDRs) and gridded data. Our approach was to: (1) search for depositional fans at the slope break between the rift apron and adjacent low slope plains; (2) determine if there is evidence that previously formed deposits might have been buried by plains units; (3) characterize the Tharsis plains morphologies east of Pavonis Mons; and (4) assess their relationship to the rift apron units. We have not identified topographically significant depositional fans, nor did we observe evidence to suggest that plains units have buried older rift apron units. Flow features associated with the rift apron are observed to continue across the slope break onto the plains. In this area, the plains are composed of a variety of small fissures and low shield vents around which broad channel-fed and tube-fed flows have been identified. We also find broad, flat-topped plateaus and sinuous ridges mixed among the channels, tubes and vents. Flat-topped plateaus and sinuous ridges are morphologies that are analogous to those observed on the coastal plain of Hawai‘i, where lava flows have advanced from the volcano's several degree flank onto the nearly zero degree coastal plain. When local volumetric flow rates are low, flow fronts tend to spread laterally and often thicken via endogenous growth, or inflation, of the sheet-like flow units. If flow advance is restricted by existing topography into narrow pathways, inflation can be focused into sinuous, elongate ridges. The presence of plateaus and ridges—emplaced from the rift zones, across the plains to the east of Pavonis Mons—and a lack of fan-like features, or evidence for their burial, are consistent with rift apron lavas crossing a slope break with low local volumetric flow rates that led to inflation of sheet-like and tube-fed lava flows.

Journal of Volcanology and Geothermal Research↗
USGS2019

Twenty years (1990–2010) of geodetic monitoring of Galeras volcano (Colombia) from continuous tilt measurements.

Galeras - an andesitic stratovolcano part of the Galeras Volcanic Complex - is one of the most active volcanoes in Colombia. Historic activity is centered on a small-volume cone inside the youngest amphitheater, which breaches the west flank of the volcano. At least 30 confirmed eruption periods have been recorded in the past 480 years, with episodes of unrest ranging from weak fumarolic activity and ash emissions to larger explosive events. The most recent eruption periods, recorded instrumentally since 1988, have been characterized by minor explosive eruptions, and the emplacement of three crater domes and small pyroclastic flow deposits. In this paper, we discuss the evolution of volcanic activity using a 20-year-long record of tilt measurements. In particular, we focus on three episodes of unrest occurred in 1991, 2006 and 2008, when the deformation was clearly associated with shallow magma intrusions, and the emplacement and destruction of crater domes. The depth of the intrusions varied from a few hundred meters (August 2005) to two kilometers (January 2009), while the volume change ranged from 104 m3 (May–October 2009) to 106 m3 (January 2009). A comparison with seismic data indicates that the deformation sources were located within the cloud of hypocenters of the volcano-tectonic events. The lack of a clear correlation between the volume change (and depth) of the sources and the total SO2 flux could indicate that the unrest at Galeras was related to a larger intrusive event with only a small part of the magma erupted in the form of tephra and lava domes.

Galeras volcano↗
USGS2019

Using a process-based model of pre-eruptive seismic patterns to forecast evolving eruptive styles at Sinabung Volcano, Indonesia

Most volcanoes worldwide are not monitored in real-time; for those that are, patterns of pre-eruptive earthquakes coupled with conceptual models of magma ascent enable short-term forecasting of eruption onset. Basic event locations, characterization of background seismicity, and recognition of changes in earthquake types and energy release are most important to successful eruption forecasting. During renewed activity at Sinabung volcano, Indonesia, this approach was used by the Center for Volcanology and Geological Hazards Mitigation (CVGHM) and the USGS Volcano Disaster Assistance Program to forecast eruption onset, identify changes in eruptive styles and raise or lower alert levels and extend or contract evacuation zones. After > 400 years of quiescence, Sinabung began erupting in August 2010. The volcano was unmonitored at the onset of these eruptions, which were phreatic, but soon after a monitoring network was installed by CVGHM. Increasing swarms of high-frequency volcano tectonic (VT) earthquakes were used to forecast continuing phreatic eruptions. Volcanic activity decreased in mid-September 2010, while additional intrusions at depth (inferred from continued distal VT swarms) continued through September 2013, when explosive phreatic eruptions recurred. Explosive eruptions were forecast based on increases in real-time seismic amplitude measurement (RSAM) and VT seismicity. Seismicity changed markedly in late November and early December 2013 with the occurrence of deep earthquakes and an overall transition from low-frequency (LF) dominated and irregular (in time and magnitude) earthquakes to more regular LF and hybrid seismicity – a transition that accompanied the continued rise, eventual emergence and growth of a lava dome in the summit crater. This lava dome was first observed on 18 December. In late December 2013 to early January 2014, the eruptive style changed again as additional ascending magma deformed the summit and the dome grew beyond the capacity of the summit crater, resulting in the en masse collapse of the lava dome (2 Mm3) on 11 January and the largest pyroclastic flow to date. The collapse was forecast on the basis of a several order of magnitude increase in RSAM, continued strong distal VT seismicity, an increase in proximal seismicity, and large-scale observed deformation of the summit area. Similarly, a later collapse of a second summit lava dome on 1 February 2014 was forecast on the basis of increased distal seismicity. Here, we demonstrate how a process-based volcano seismicity model was used in combination with real-time data to forecast the time and magnitude of eruptions, as well as changes in eruption style.

Sinabung Volcano↗
USGS2017

Revised tephra volumes for Cascade Range volcanoes

Isopach maps from tephra eruptions from Mount St. Helens were reported in Carey et al. (1995) and for tephra eruptions from Glacier Peak in Gardner et al. (1998). For exponential thinning, the isopach data only define a single slope on a log thickness versus square root of area plot. Carey et al. (1995) proposed a model that was used to estimate a second slope, and volumes were presented in both studies using this model. A study by Sulpizio (2005) for estimating the second slope and square root of area where the lines intersect involves a systematic analysis of many eruptions to provide correlation equations. The purpose of this paper is to recalculate the volumes of Cascades eruptions and compare results from the two methods. In order to gain some perspective on the methods for estimating the second slope, we use data for thickness versus distance beyond the last isopach that are available for some of the larger eruptions in the Cascades. The thickness versus square root of area method is extended to thickness versus distance by developing an approximate relation between the two assuming elliptical isopachs with the source at one of the foci. Based on the comparisons made between the Carey et al. (1995) and Sulpizio (2005) methods, it is felt that the later method provides a better estimate of the second slope. For Mount St. Helens, the estimates of total volume using the Sulpizio (2005) method are generally smaller than those using the Carey et al. (1995) method. For the volume estimates of Carey et al. (1995), the volume of the May 18, 1980, eruption of Mount St. Helens is smaller than six of the eight previous eruptions. With the new volumes using the Sulpizio (2005) method, the 1980 eruption is smaller in volume than the upper end of the range for only three of the layers (Wn, Ye, and Yn) and is the same size as layer We. Thus the 1980 eruption becomes representative of the mid-range of volumes rather than being in the lower range.

Cascade Range volcanoes↗
USGS2017

The spectrum of persistent volcanic flank instability: A review and proposed framework based on Kīlauea, Piton de la Fournaise, and Etna

Persistent motion of the south flank of Kīlauea Volcano, Hawai'i, has been known for several decades, but has only recently been identified at other large basaltic volcanoes—namely Piton de la Fournaise (La Réunion) and Etna (Sicily)—thanks to the advent of space geodetic techniques. Nevertheless, understanding of long-term flank instability is based largely on the example of Kīlauea, despite the large differences in the manifestations and mechanisms of the process when viewed through a comparative lens. For example, the rate of flank motion at Kīlauea is several times that of Etna and Piton de la Fournaise and is accommodated on a slip plane several km deeper than is probably present at the other two volcanoes. Gravitational spreading also appears to be the dominant driving force at Kīlauea, given the long-term steady motion of the volcano's south flank regardless of eruptive/intrusive activity, whereas magmatic activity plays a larger role in flank deformation at Etna and Piton de la Fournaise. Kīlauea and Etna, however, are both characterized by heavily faulted flanks, while Piton de la Fournaise shows little evidence for flank faulting. A helpful means of understanding the spectrum of persistent flank motion at large basaltic edifices may be through a framework defined on one hand by magmatic activity (which encompasses both magma supply and edifice size), and on the other hand by the structural setting of the volcano (especially the characteristics of the subvolcanic basement or subhorizontal intravolcanic weak zones). A volcano's size and magmatic activity will dictate the extent to which gravitational and magmatic forces can drive motion of an unstable flank (and possibly the level of faulting of that flank), while the volcano's structural setting governs whether or not a plane of weakness exists beneath or within the edifice and can facilitate flank slip. Considering persistent flank instability using this conceptual model is an alternative to using a single volcano as a “type example”—especially given that the example is usually Kīlauea, which defines an extreme end of the spectrum—and can provide a basis for understanding why flank motion may or may not exist on other large basaltic volcanoes worldwide.

Etna, Kīlauea Volcano, Piton de la Fournaise↗
USGS2017

Fluid-driven uplift at Long Valley Caldera, California: Geologic perspectives

Since persistent seismicity began in the Sierra Nevada adjacent to Long Valley caldera in 1978–1980, intracaldera unrest has been marked by (1) episodes of uplift totaling ~ 83 cm, centered on the middle Pleistocene resurgent dome, and (2) recurrent earthquake swarms along a 12-km-long segment of the caldera's ring-fault zone that is contiguous with both the dome and the Sierran seismogenic domain. Others have attributed the recent unrest to magmatic intrusion(s), but it is argued here that evidence for new magma is lacking and that ongoing uplift and ring-fault-zone seismicity are both promoted by ascent of aqueous fluid released by second boiling of the residue of the enormous Pleistocene rhyolitic reservoir terminally crystallizing at depths ≥ 10 km. For 2 Myr, eruptive vent clusters migrated southwestward from Glass Mountain to Mammoth Mountain. There has been no eruption on the resurgent dome since 500 ka, and since 230 ka volcanism has been restricted to the caldera's west moat and contiguous Sierran terrain, both outside the structural caldera. High-temperature hydrothermal activity in the central caldera waned after ~ 300 ka, cooling the Pleistocene rhyolitic focus to the extent that drilling on the resurgent dome found mid-caldera temperature to be only 100 °C and isothermal at depths of 2–3 km. Beneath most of the resurgent dome, there is little seismicity at any depth, no emission of magmatic CO 2 or other magmatic gases, no elevated 3 He/ 4 He ratios, and only normal to below-normal heat flow. Most of the 75-km-long ring-fault zone is likewise aseismic, excepting only the 12-km segment contiguous with the extracaldera seismogenic domain in the Sierra. Since 1980, the Sierran seismicity has released 3.6 times more cumulative seismic energy than have intracaldera earthquakes. The caldera seismicity is not driven by stresses associated with the adjacent uplift but, instead, by the extracaldera tectonic stressfield. Sierran seismicity activated the directly contiguous south-moat segment of the ring-fault zone, which had originated in the caldera-forming eruption at 767 ka and everywhere else remains sealed. Hypocenter relocation studies of 1000s of earthquakes along the seismic segment have resolved recurrent upward-migrating swarms within networks of cryptic faults, apparently triggered by rapidly ascending pulses of high-pressure low-viscosity aqueous fluid. Entering the brittle crust at depths of 8–10 km, such fluid is just what should be expected from second boiling of the late-stage CO 2 -poor rhyolitic residue. The fluid provides the pressure source above the apex of the crystallizing caldera-wide pluton and then escapes laterally to the newly reactivated southern segment of the ring-fault zone, its only available permeable pathway, where it mediates the ongoing south-moat seismicity.

California↗
USGS2017

Evidence for degassing of fresh magma during the 2004-2008 eruption of Mount St. Helens: Subtle signals from the hydrothermal system

Results from chemical and isotopic analyses of water and gas collected between 2002 and 2016 from sites on and around Mount St. Helens are used to assess magmatic degassing related to the 2004-2008 eruption. During 2005 the chemistry of hot springs in The Breach of Mount St. Helens showed no obvious response to the eruption, and over the next few years, changes were subtle, giving only slight indications of perturbations in the system. By 2010 however, water chemistry, temperatures, and isotope compositions (δD and δ 18 O) clearly indicated some inputs of volatiles and heat associated with the eruption, but the changes were such that they could be attributed to a pre-existing, gas depleted magma. An increase of ~ 1.5‰ in the δ 13 C values of dissolved carbon in the springs was noted in 2006 and continued through 2009, a change that was mirrored by a similar shift in δ 13 C-CO 2 in bubble gas emissions. These changes require input of a new source of carbon to the hydrothermal system and provide clear evidence of CO 2 from an undegassed body of magma. Rising trends in 3 He/ 4 He ratios in gas also accompanied the increases in δ 13 C. Since 2011 maximum R C /R A values are ≥ 6.4 and are distinctly higher than 5 samples collected between 1986 and 2002, and provide additional evidence for some involvement of new magma as early as 2006, and possibly earlier, given the unknown time needed for CO 2 and He to traverse the system and arrive at the springs.

Washington↗
USGS2017

Causes of distal volcano-tectonic seismicity inferred from hydrothermal modeling

Distal volcano-tectonic (dVT) seismicity typically precedes eruption at long-dormant volcanoes by days to years. Precursory dVT seismicity may reflect magma-induced fluid-pressure pulses that intersect critically stressed faults. We explored this hypothesis using an open-source magmatic-hydrothermal code that simulates multiphase fluid and heat transport over the temperature range 0 to 1200 °C. We calculated fluid-pressure changes caused by a small (0.04 km 3 ) intrusion and explored the effects of flow geometry (channelized vs. radial flow), magma devolatilization rates (0–15 kg/s), and intrusion depths (5 and 7.5 km, above and below the brittle-ductile transition). Magma and host-rock permeabilities were key controlling parameters and we tested a wide range of permeability ( k ) and permeability anisotropies ( k h / k v ), including k constant, k ( z ), k ( T ), and k ( z , T , P ) distributions, examining a total of ~ 1600 realizations to explore the relevant parameter space. Propagation of potentially causal pressure changes (Δ P ≥ 0.1 bars) to the mean dVT location (6 km lateral distance, 6 km depth) was favored by channelized fluid flow, high devolatilization rates, and permeabilities similar to those found in geothermal reservoirs ( k ~ 10 − 16 to 10 − 13 m 2 ). For channelized flow, magma-induced thermal pressurization alone can generate cases of ∆ P ≥ 0.1 bars for all permeabilities in the range 10 − 16 to 10 − 13 m 2 , whereas in radial flow regimes thermal pressurization causes ∆ P < 0.1 bars for all permeabilities. Changes in distal fluid pressure occurred before proximal pressure changes given modest anisotropies ( k h / k v ~ 10–100). Invoking k ( z , T , P ) and high, sustained devolatilization rates caused large dynamic fluctuations in k and P in the near-magma environment but had little effect on pressure changes at the distal dVT location. Intrusion below the brittle-ductile transition damps but does not prevent pressure transmission to the dVT site.

Journal of Volcanology and Geothermal Research↗
USGS2019

Overview of the eruptions of Sinabung Volcano, 2010 and 2013-present, and details of the 2013 phreatomagmatic phase

A small phreatic eruption of Sinabung Volcano, North Sumatra on August 28, 2010, at 18:30 local time marked the first eruption in the past ~1200 years. The eruption took place from two small vents in the south crater area. Explosions and ash emissions from these vents generated multiple ash plumes that reached altitudes of up to 5 km during early- to mid-September. By the end of September, only low level steam plumes were visible and the alert level was reduced from Level 4 (highest) to Level 3. The 2010 eruption effectively ended at this time. Beginning two days after the initial 2010 eruption, activity of the eruption has been monitoring continuously by a telemetered seismic network surrounding the volcano and by remotely sensed observations. This monitoring system was supplemented with a near-field continuous GPS network, beginning in February 2011. Persistent fumarolic emissions continued for almost 3 years following the 2010 eruption, before a new eruption began on 15 September 2013. This eruption continues to the present. The ongoing eruption is divided into 5 major phases: 1) phreatomagmatic phase (July 2013 - 18 December 2013); 2) first dome and collapse phase with pyroclastic density currents (PDCs; block-and-ash flows and related surges) to south (18 December 2013 - 10 January 2014); 3) lava-flow and collapse phase (10 January 2014 - mid-September 2014); 4) second lava dome and collapse phase with PDCs to south (mid-September 2014 - July 2015); 5) lava dome collapse and ash explosion phase with PDCs to southeast and east (August 2015 - present). The volcano erupted intermittently during the early phreatomagmatic phase with small vertical ash explosions. Then the eruption became increasingly vigorous with more repetitive and intense vertical ash explosions during late October through November. The first small pyroclastic density currents (PDCs) began on November 1. These pyroclastic flows descended the southeastern flank to a distance of 2 km.

Journal of Volcanology and Geothermal Research↗
USGS2017

Volcano Geodesy: Recent developments and future challenges

Ascent of magma through Earth's crust is normally associated with, among other effects, ground deformation and gravity changes. Geodesy is thus a valuable tool for monitoring and hazards assessment during volcanic unrest, and it provides valuable data for exploring the geometry and volume of magma plumbing systems. Recent decades have seen an explosion in the quality and quantity of volcano geodetic data. New datasets (some made possible by regional and global scientific initiatives), as well as new analysis methods and modeling practices, have resulted in important changes to our understanding of the geodetic characteristics of active volcanism and magmatic processes, from the scale of individual eruptive vents to global compilations of volcano deformation. Here, we describe some of the recent developments in volcano geodesy, both in terms of data and interpretive tools, and discuss the role of international initiatives in meeting future challenges for the field.

Journal of Volcanology and Geothermal Research↗
USGS2019

Joint 3-D tomographic imaging of Vp, Vs and Vp/Vs and hypocenter relocation at Sinabung volcano, Indonesia from November to December 2013

We conducted travel time tomography using P- and S-wave arrival times of volcanic-tectonic (VT) events that occurred between November and December 2013 to determine the three-dimensional (3D) seismic velocity structure (Vp, Vs, and Vp/Vs) beneath Sinabung volcano, Indonesia in order to delineate geological subsurface structure and to enhance our understanding of the volcanism itself. This was a time period when phreatic explosions became phreatomagmatic and then magma migrated to the surface forming a summit lava dome. We used 4846 VT events with 16,138 P- and 16,138 S-wave arrival time phases recorded by 6 stations for the tomographic inversion. The relocated VTs collapse into three clusters at depths from the surface to sea level, from 2 to 4 km below sea level, and from 5 to 8.5 km below sea level. The tomographic inversion results show three prominent regions of high Vp/Vs (~ 1.8) beneath Sinabung volcano at depths consistent with the relocated earthquake clusters. We interpret these anomalies as intrusives associated with previous eruptions and possibly surrounding the magma conduit, which we cannot resolve with this study. One anomalous region might contain partial melt, at sea level and below the eventual eruption site at the summit. Our results are important for the interpretation of a conceptual model of the “plumbing system” of this hazardous volcano.

Sinabung Volcano↗