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 253 records · Page 14Linked to original sources
USGS2017

Monitoring gas and heat emissions at Norris Geyser Basin, Yellowstone National Park, USA based on a combined eddy covariance and Multi-GAS approach

We quantified gas and heat emissions in an acid-sulfate, vapor-dominated area (0.04-km 2 ) of Norris Geyser Basin, located just north of the 0.63 Ma Yellowstone Caldera and near an area of anomalous uplift. From 14 May to 3 October 2016, an eddy covariance system measured half-hourly CO 2 , H 2 O and sensible ( H ) and latent ( LE ) heat fluxes and a Multi-GAS instrument measured (1 Hz frequency) atmospheric H 2 O, CO 2 and H 2 S volumetric mixing ratios. We also measured soil CO 2 fluxes using the accumulation chamber method and temperature profiles on a grid and collected fumarole gas samples for geochemical analysis. Eddy covariance CO 2 fluxes ranged from − 56 to 885 g m − 2 d − 1 . Using wavelet analysis, average daily eddy covariance CO 2 fluxes were locally correlated with average daily environmental parameters on several-day to monthly time scales. Estimates of CO 2 emission rate from the study area ranged from 8.6 t d − 1 based on eddy covariance measurements to 9.8 t d − 1 based on accumulation chamber measurements. Eddy covariance water vapor fluxes ranged from 1178 to 24,600 g m − 2 d − 1 . Nighttime H and LE were considered representative of hydrothermal heat fluxes and ranged from 4 to 183 and 38 to 504 W m − 2 , respectively. The total hydrothermal heat emission rate ( H + LE + radiant) estimated for the study area was 11.6 MW and LE contributed 69% of the output. The mean ± standard deviation of H 2 O, CO 2 and H 2 S mixing ratios measured by the Multi-GAS system were 9.3 ± 3.1 parts per thousand, 467 ± 61 ppmv, and 0.5 ± 0.6 ppmv, respectively, and variations in the gas compositions were strongly correlated with diurnal variations in environmental parameters (wind speed and direction, atmospheric temperature). After removing ambient H 2 O and CO 2 , the observed variations in the Multi-GAS data could be explained by the mixing of relatively H 2 O-CO 2 -H 2 S-rich fumarole gases with CO 2 -rich and H 2 O-H 2 S-poor soil gases. The fumarole H 2 O/CO 2 and CO 2 /H 2 S end member ratios (101.7 and 27.1, respectively, on average) were invariant during the measurement period and fell within the range of values measured in direct fumarole gas samples. The soil gas H 2 O/CO 2 end member ratios (~ 15–30) were variable and low relative to the fumarole end member, likely resulting from water vapor loss during cooling and condensation in the shallow subsurface, whereas the CO 2 /H 2 S end member ratio was high (~ 160), presumably related to transport of CO 2 -dominated soil gas emissions mixed with trace fumarolic emissions to the Multi-GAS station. Nighttime eddy covariance ratios of H 2 O to CO 2 flux were typically between the soil gas and fumarole end member H 2 O/CO 2 ratios defined by Multi-GAS measurements. Overall, the combined eddy covariance and Multi-GAS approach provides a powerful tool for quasi-continuous measurements of gas and heat emissions for improved volcano-hydrothermal monitoring.

Wyoming↗
USGS2018

Insights into the emplacement of upper-crustal plutons and their relationship to large silicic calderas, from field relationships, geochronology, and zircon trace element geochemistry in the Stillwater – Clan Alpine caldera complex, western Nevada, USA

Geologic mapping, new U-Pb zircon ages, and new and published 40 Ar/ 39 Ar sanidine ages document the timing and extent of Oligocene magmatism in the southern Stillwater Range and Clan Alpine Mountains of western Nevada, where Miocene extension has exposed at least six nested silicic calderas and underlying granitic plutons to crustal depths locally ≥ 9 km. Both caldera-forming rhyolitic tuffs and underlying plutons were emplaced in two episodes, one from about 30.4–28.2 Ma that included the Deep Canyon, Job Canyon, and Campbell Creek calderas and underlying plutons, and one from about 25.3–24.8 Ma that included the Louderback Mountains, Poco Canyon, and Elevenmile Canyon calderas and underlying plutons. In these two 1–2 m.y. periods, almost the entire Mesozoic upper crust was replaced by Oligocene intrusive and extrusive rocks to depths ≥ 9 km over an estimated total area of ~ 1500 km 2 (pre-extension). Zircon trace element geochemistry indicates that some plutonic rock can be solidified residual magma from the tuff eruptions. Most plutons are not solidified residual magma, although they directly underlie calderas and were emplaced along the same structures shortly after to as much as one million years after caldera formation. Magma chambers and plutons grew by floor subsidence accommodated by downward transfer of country rocks. If other Great Basin calderas are similar, the dense concentration of shallowly exposed calderas in central Nevada is underlain by a complexly zoned mid-Cenozoic batholith assembled in discrete pulses that coincided with formation of large silicic calderas up to 2500–5000 km 3 .

Nevada↗
USGS2019

Gas and ash emissions associated with the 2010–present activity of Sinabung Volcano, Indonesia

Sinabung Volcano (Sumatra, Indonesia) awoke from over 1200 years of dormancy with multiple phreatic explosions in 2010. After a period of quiescence, Sinabung activity resumed in 2013, producing frequent explosions, lava dome extrusion, and pyroclastic flows from dome collapses, becoming one of the world's most active volcanoes and displacing over 20,000 citizens. This study presents a compilation of the geochemical datasets collected by the Indonesian Center for Volcanology and Geological Hazard Mitigation (CVGHM) from 2010 - current (2016), which provides insights into the evolution of the eruption. Based on observations of SO2 emissions, ash componentry, leachate chemistry, and bulk ash geochemistry, the eruption can be split into five distinct phases. The initial stage of phreatic summit explosions occurred from August - October 2010, during which background SO2 emissions averaged ~550 ± 180 t/d (1 s.d.). An eruptive pause (phase two) starting in October 2010 abruptly ended in September 2013 with a resumption of conduit-clearing eruptions. This third phase had a relatively modest background SO2 emission rate (avg. ~410 ± 275 t/d) and produced ash consisting entirely of accidental ejecta with high S/Cl leachate ratios (up to 30), suggestive of deep-sourced magma and the incorporation of hydrothermal sulfur-bearing phases. The most intense phase of the eruption (phase four) occurred from December 2013 to February 2014, when juvenile magma first reached the surface. This period included dozens of large eruptions per day, high SO2 emission rates (average: 1,120 ± 1,030 t/d, peak: ~3,800 t/d), the onset of lava dome extrusion, and a dramatic drop in S/Cl ash leachates to ratios < 5, all reflecting increased degassing from shallow magma and the clearing out of sulfurous phases from the old hydrothermal system. From late February 2014 through the time of writing (September 2016), Sinabung settled into a relatively steady state of lower activity (phase five). Ash emissions now consist of dominantly juvenile material, and background SO2 emission rates have been progressively decreasing to an average of ~250 - 300 t/d. Starting August 2016, SO2 emissions started being measured in a continuous manner using a network of permanent scanning DOAS instruments. We find that long-term SO2 emission rates have been gradually declining at Sinabung since early 2014, consistent with an apparent decrease in magma supply. Our degassing model suggests that large explosions and pyroclastic flows could continue in the near-term owing to conduit plugging and dome collapses, remaining a major threat until the magma supply rate decreases further and the eruption ends.

Sumatra↗
USGS2018

The timing and origin of pre- and post-caldera volcanism associated with the Mesa Falls Tuff, Yellowstone Plateau volcanic field

We present new sanidine 40 Ar/ 39 Ar ages and paleomagnetic data for pre- and post-caldera rhyolites from the second volcanic cycle of the Yellowstone Plateau volcanic field, which culminated in the caldera-forming eruption of the Mesa Falls Tuff at ca. 1.3 Ma. These data allow for a detailed reconstruction of the eruptive history of the second volcanic cycle and provide new insights into the petrogenesis of rhyolite domes and flows erupted during this time period. 40 Ar/ 39 Ar age data for the biotite-bearing Bishop Mountain flow demonstrate that it erupted approximately 150 kyr prior to the Mesa Falls Tuff. Integrating 40 Ar/ 39 Ar ages and paleomagnetic data for the post-caldera Island Park rhyolite domes suggests that these five crystal-rich rhyolites erupted over a centuries-long time interval at 1.2905 ± 0.0020 Ma (2σ). The biotite-bearing Moonshine Mountain rhyolite dome was originally thought to be the downfaulted vent dome for the pre-caldera Bishop Mountain flow due to their similar petrographic and oxygen isotope characteristics, but new 40 Ar/ 39 Ar dating suggest that it erupted near contemporaneously with the Island Park rhyolite domes at 1.2931 ± 0.0018 Ma (2σ) and is a post-caldera eruption. Despite their similar eruption ages, the Island Park rhyolite domes and the Moonshine Mountain dome are chemically and petrographically distinct and are not derived from the same source. Integrating these new data with field relations and existing geochemical data, we present a petrogenetic model for the formation of the post-Mesa Falls Tuff rhyolites. Renewed influx of basaltic and/or silicic recharge magma into the crust at 1.2905 ± 0.0020 Ma led to [1] the formation of the Island Park rhyolite domes from the source region that earlier produced the Mesa Falls Tuff and [2] the formation of Moonshine Mountain dome from the source region that earlier produced the biotite-bearing Bishop Mountain flow. These magmas were stored in the crust for less than a few thousand years before being erupted contemporaneously along a 30 km long, structurally controlled vent zone related to extracaldera Basin and Range faults. These data highlight the rapidity with which magma can be generated and erupted over large distances at Yellowstone.

Idaho, Wyoming↗
USGS2018

Monogenetic origin of Ubehebe Crater maar volcano, Death Valley, California: Paleomagnetic and stratigraphic evidence

Paleomagnetic data for samples collected from outcrops of basaltic spatter at the Ubehebe Crater cluster, Death Valley National Park, California, record a single direction of remanent magnetization indicating that these materials were emplaced during a short duration, monogenetic eruption sequence ~ 2100 years ago. This conclusion is supported by geochemical data encompassing a narrow range of oxide variation, by detailed stratigraphic studies of conformable phreatomagmatic tephra deposits showing no evidence of erosion between layers, by draping of sharp rimmed craters by later tephra falls, and by oxidation of later tephra layers by the remaining heat of earlier spatter. This model is also supported through a reinterpretation and recalculation of the published age results (Sasnett et al., 2012) from an innovative and bold exposure-age study on very young materials. Their conclusion of multiple and protracted eruptions at Ubehebe Crater cluster is here modified through the understanding that some of their quartz-bearing clasts inherited from previous exposure on the fan surface (too old), and that other clasts were only exposed at the surface by wind and/or water erosion centuries after their eruption (too young). Ubehebe Crater cluster is a well preserved example of young monogenetic maar type volcanism protected within a National Park, and it represents neither a protracted eruption sequence as previously thought, nor a continuing volcanic hazard near its location.

California↗
USGS2018

Focused seismicity triggered by flank instability on Kīlauea's Southwest Rift Zone

Swarms of earthquakes at the head of the Southwest Rift Zone on Kīlauea Volcano, Hawaiʻi, reveal an interaction of normal and strike-slip faulting associated with movement of Kīlauea's south flank. A relocated subset of earthquakes between January 2012 and August 2014 are highly focused in space and time at depths that are coincident with the south caldera magma reservoir beneath the southern margin of Kīlauea Caldera. Newly calculated focal mechanisms are dominantly dextral shear with a north-south preferred fault orientation. Two earthquakes within this focused area of seismicity have normal faulting mechanisms, indicating two mechanisms of failure in very close proximity (10's of meters to 100 m). We suggest a model where opening along the Southwest Rift Zone caused by seaward motion of the south flank permits injection of magma and subsequent freezing of a plug, which then fails in a right-lateral strike-slip sense, consistent with the direction of movement of the south flank. The seismicity is concentrated in an area where a constriction occurs between a normal fault and the deeper magma transport system into the Southwest Rift Zone. Although in many ways the Southwest Rift Zone appears analogous to the more active East Rift Zone, the localization of the largest seismicity (>M2.5) within the swarms to a small volume necessitates a different model than has been proposed to explain the lineament outlined by earthquakes along the East Rift Zone.

Hawaii↗
USGS2019

Construction of probabilistic event trees for eruption forecasting at Sinabung volcano, Indonesia 2013–14

Eruptions of Sinabung volcano , Indonesia have been ongoing since 2013. Since that time, the character of eruptions has changed, from phreatic to phreatomagmatic to magmatic explosive eruptions, and from production of a lava dome that collapsed to a subsequent thick lava flow that slowly ceased to be active, and later, to a new lava dome. As the eruption progressed, event trees were constructed to forecast eruptive behavior six times, with forecast windows that ranged from 2 weeks to 1 year: November 7–10, December 12–14, and December 27, 2013; and January 9–10, May 13, and October 7, 2014. These event trees were successful in helping to frame the forecast scenarios, to collate current monitoring information, and to document outstanding questions and unknowns. The highest probability forecasts closely matched outcomes of eruption size (including extrusion of the first dome), production of pyroclastic density currents, and pyroclastic density current runout distances. Events assigned low probabilities also occurred, including total collapse of the lava dome in January 2014 and production of a small blast pyroclastic density current in February 2014.

Sinabung volcano↗
USGS2018

Multi-year high-frequency hydrothermal monitoring of selected high-threat Cascade Range volcanoes

From 2009 to 2015 the U.S. Geological Survey (USGS) systematically monitored hydrothermal behavior at selected Cascade Range volcanoes in order to define baseline hydrothermal and geochemical conditions. Gas and water data were collected regularly at 25 sites on 10 of the highest-risk volcanoes in the Cascade Range. These sites include near-summit fumarole groups and springs/streams that show clear evidence of magmatic influence (high 3 He/ 4 He ratios and/or large fluxes of magmatic CO 2 or heat). Site records consist mainly of hourly temperature and hydrothermal-flux data. Having established baseline conditions during a multiyear quiescent period, the USGS reduced monitoring frequency from 2015 to present. The archived monitoring data are housed at (doi:10.5066/F72N5088). These data (1) are suitable for retrospective comparison with other continuous geophysical monitoring data and (2) will provide context during future episodes of volcanic unrest, such that unrest-related variations at these thoroughly characterized sites will be more clearly recognizable. Relatively high-frequency year-round data are essential to achieve these objectives, because many of the time series reveal significant diurnal, seasonal, and inter-annual variability that would tend to mask unrest signals in the absence of baseline data. Here we characterize normal variability for each site, suggest strategies to detect future volcanic unrest, and explore deviations from background associated with recent unrest.

California, Oregon, Washington↗
USGS2019

Application of an updated atmospheric model to explore volcano infrasound propagation and detection in Alaska

Winds and temperature gradients greatly affect the long-range propagation of infrasound . The spatio-temporal variability of these parameters must therefore be accurately characterized to correctly interpret recorded infrasound at long distances, specifically to differentiate between source and propagation effects. Here we present the first results of an open source reanalysis model, termed Alaska Volcano Observatory Ground-to-Space (AVO-G2S), constructed to accurately characterize the atmosphere and model long-range infrasound propagation from volcanic eruptions in Alaska. We select a number of case studies to examine recent eruptions of Alaskan volcanoes whose ash emissions posed a threat to air traffic, including the two most recent eruptions of Pavlof Volcano and two typical explosions from Cleveland Volcano. Strong tropospheric ducting and low noise at the station during the 21 July 2015 explosion of Cleveland Volcano led to an automated detection of the explosion at an infrasound array 992 km away, whereas low signal-to-noise ratio for the 6 November 2014 Cleveland Volcano explosion helps explain the non-detection in real-time of a predicted strong stratospheric arrival. For the November 2014 Pavlof eruption, discrepancies between local seismic data and a distal infrasound array 460 km away cannot be solely explained by changes in atmospheric conditions, though some features of the complex propagation predictions follow the trends in long-range infrasound signals. The most recent eruption of Pavlof Volcano in March 2016 shows minimal changes in propagation conditions throughout the eruption and therefore indicates that the signals detected at long-range primarily reflect source processes. These results show how detailed examination of the acoustic propagation conditions provides insight into detection capability and eruption dynamics. Future work will implement AVO-G2S and high-resolution long-range infrasound propagation modeling in real-time for Alaskan volcanoes of interest.

Alaska↗
USGS2018

Three-dimensional geophysical mapping of shallow water saturated altered rocks at Mount Baker, Washington: Implications for slope stability

Water-saturated hydrothermal alteration reduces the strength of volcanic edifices, increasing the potential for catastrophic sector collapses that can lead to far traveled and destructive debris flows. Intense hydrothermal alteration significantly lowers the resistivity and magnetization of volcanic rock and therefore hydrothermally altered rocks can be identified with helicopter electromagnetic and magnetic measurements. Geophysical models constrained by rock properties and geologic mapping show that intensely altered rock is restricted to two small (500 m diameter), >150 m thick regions around Sherman Crater and Dorr Fumarole Field at Mount Baker, Washington. This distribution of alteration contrasts with much thicker and widespread alteration encompassing the summits of Mounts Adams and Rainier prior to the 5600 year old Osceola collapse, which is most likely due to extreme erosion and the limited duration of summit magmatism at Mount Baker. In addition, the models suggest that the upper ~300 m of rock contains water which could help to lubricate potential debris flows. Slope stability modeling incorporating the geophysically modeled distribution of alteration and water indicates that the most likely and largest (~0.1 km 3 ) collapses are from the east side of Sherman Crater. Alteration at Dorr Fumarole Field raises the collapse hazard there, but not significantly because of its lower slope angles. Geochemistry and analogs from other volcanoes suggest a model for the edifice hydrothermal system.

Washington↗
USGS2019

Monitoring, forecasting collapse events, and mapping pyroclastic deposits at Sinabung volcano with satellite imagery

During the ongoing (2013–present) eruption of Sinabung volcano, north Sumatra, we have routinely used a variety of satellite remote sensing data to observe and forecast lava dome and lava flow collapse events, to map the resulting pyroclastic deposits, and to estimate effusion rates. In this paper, we focus on the first two years of the current eruption (September 2013–December 2015), and we summarize major events in 2016. We divide the eruption into 5 major phases: 1) phreatomagmatic (July 2013–18 December 2013), 2) first dome growth and collapse (18 December 2013–10 January 2014), 3) lava-flow (10 January 2014–mid-September 2014), 4) second lava dome and collapse (mid-September 2014–July 2015), 5) lava dome collapse and ash explosion phase (August 2015–present). Throughout the eruption, remotely sensed information has been instrumental in assessing the stability of the lava dome and flow and to forecast collapse events that produce pyroclastic density currents (PDCs: block-and-ash flows, co-ignimbrite surges, and blasts). Forecasts based on remote sensing data in combination with seismic, geodetic and gas-monitoring data have also helped inform decisions related to alert levels and evacuations. Relatively unusual aspects of the Sinabung eruption include the transition from dome to flow morphology (phase 2 to phase 3 transition) and the frequent occurrence during phase 3 of collapses from the lava flow-front and flow-margins—collapses that produced extensive pyroclastic density currents. By analogy to the well-known “Merapi type” collapses and pyroclastic deposits, we propose that lava flow-front and flow-margin collapses with associated PDCs be known as “Sinabung type.” Although detailed study of deposits has not been possible due to continuing hazards, our observations suggest that the transition from lava dome to lava flow and the occurrence of flow-front and flow-margin collapses reflect a particular combination of lava viscosity and steepness of slope. Our observations also show clear evidence of at least one slope-parallel high-velocity and dilute PDC (a “blast”) that emanated from a lava-margin collapse site 500 m downslope from the vent. This 1 February 2014 blast downed and singed a forest out to at least 3.9 km from the collapse site and killed 16 people. We also use a combination of field and remotely sensed data to map the distribution of Sinabung deposits. We estimate eruptive volumes and extrusion rates by combining sequential measurements of lava surface and pyroclastic flow areas with thickness estimates derived from simple geometric assumptions, oblique photographs and Digital Elevation Models (DEMs) derived from remotely sensed data. Our estimates of short-term effusion rates vary widely on a daily to weekly basis, from <1 to >20 m 3 s −1 . In a few cases, periods of increased extrusion precede lava flow-front collapses by a few days to a week, suggesting delays in transmittance of effusion pulses as lava moves from vent to flow front. We find that, as of 1 January 2016, the total area of deposits is 10 7 m 2 , and their approximate deposit volume is about 0.3 km 3 , equivalent to 0.2 km 3 Dense Rock Equivalent (DRE). We anticipate that our deposit maps will be valuable in the future as a framework for the study of the magmatic and textural evolution of eruptive products through time.

Sinabung volcano↗
USGS2018

A retrospective look at the February 1993 east rift zone intrusion at Kīlauea volcano, Hawaii

The February 1993 dike intrusion in the East Rift Zone (ERZ) of Kīlauea Volcano, Hawai'i, was recognized from tilt and seismic data, but ground-based geodetic data were too sparse to constrain the characteristics of the intrusion. Analysis of Interferometric Synthetic Aperture Radar (InSAR) from the Japan Aerospace Exploration Agency (JAXA) JERS-1 satellite reveals a maximum of ~30 cm of line-of-sight (LOS) displacement occurring near Makaopuhi Crater in the middle ERZ of Kīlauea. We model this deformation signal as a subvertical dike using a 3D-Mixed Boundary Element Method (3D-MBEM) paired with a nonlinear inversion algorithm to find the best-fit model. The best-fit dike is located just to the west of Makaopuhi Crater striking N50°W, extends to within 100 m of the surface, is ~1.3 km in length by ~4.2 km in width along strike, and has a total volume of ~7.4 × 10 6 m 3 . In addition, a post-intrusion interferogram from JERS-1 spanning 1993–1997 was analyzed. Guided by previous results, our model for the 4-year period consists of opening of the deep rift zones by about 0.5 m at 3–8.5 km depth beneath the Southwest Rift Zone (SWRZ), ERZ and the summit. A sub-horizontal detachment fault is connected to the seaward side of the vertical dike-like source to mimic the décollement known to exist beneath the volcano. We classify the 1993 dike intrusion as a passive intrusion similar to those that occurred in 1997 and 1999. Passive intrusions lack precursory inflation at Kīlauea's summit, and the likely triggering mechanism is persistent deep rift opening combined with seaward motion of the south flank along the basal décollement. Passive intrusions make forecasting and hazard assessment difficult since they are not preceded by inflation nor by large increases in seismicity.

Hawaii↗
USGS2018

Linking the Ukinrek 1977 maar-eruption observations to the tephra deposits: New insights into maar depositional processes

The Ukinrek Maars erupted 30 March to 9 April 1977, forming two maars, a line of small pit craters and a tephra blanket extending to ~2 km from the vents. We combine photographic and written observations with stratigraphic analysis to reconstruct the eruption. The eruption began with very low (a few meters high) fountaining from small craters above an inferred east-west-trending dike, creating local scoria/spatter agglomerate ramparts with a sandy matrix. The eruption very quickly (in minutes to hours) centered on the West Maar. The West Maar eruption lasted 1–2 days, starting and ending with phreatomagmatic explosions with weak phreato-Strombolian activity in between. Initial explosions formed a 30-m-wide crater, enlarged by crater-wall collapse, and columns as high as 6500 m. Phreato-Strombolian activity produced ~72% of the erupted volume, including a small spatter cone and a scoria blanket around the vent. A final explosion series emplaced a lithic-rich breccia as ballistic blocks, possibly as the northern half of the final crater collapsed into the southern vent area. The East Maar formed over the last nine days of the eruption and represents ~93% of the total volume (4.6 × 10 6 m 3 ) of the Ukinrek eruption. Initial explosions were probably shallower than 10–20 m but most of the eruption occurred from explosions at 50–60 m below the pre-eruptive surface, with evidence of explosions to 90 m depth only at the very end of the eruption. The East Maar eruption mostly produced columns of lapilli, ash, and steam and the deposits are mostly fallout. Winds blew fallout mostly to the north for the first 5–6 days and to the south for the last three days of the eruption. Wind-directed pyroclastic density currents collapsed from the column, producing fines-rich layers within the coarser fallout. Sporadic explosions produced weak density currents in the first few days and lithic-and juvenile-block-rich breccias in the last few days of the eruption. We interpret that collapse of the crater walls made a slurry that in part provided the water for phreatomagmatic interaction. Explosions came from depths <90 m below the pre-eruptive surface except for a few explosions at the end of the eruption, with most occurring at <70 m depth. The East Maar crater was open to 40–60 m depth throughout most of the eruption, so the explosions were rarely, if ever, deeper than 30 m below the crater floor. Thus, we infer there is no classic, well-formed diatreme structure below the maar. Collapse of the East Maar crater walls provided a supply of water-saturated sediment for much of the phreatomagmatic activity, which came from two vents that did not migrate much, if at all, during the eruption. The Ukinrek Maars deposits were nearly entirely emplaced by fallout, rather than density currents, from explosions and low columns.

Alaska↗
USGS2018

Hydrothermal discharge from the El Tatio basin, Atacama, Chile

El Tatio in northern Chile is one of the best-studied geothermal fields in South America. However, there remain open questions about the mass and energy budgets, water recharge rates and residence time in the subsurface, origin of dissolved solutes, and processes affecting the phase and chemical composition of groundwater and surface water. We measured and sampled surface manifestations of the geothermal system (geysers perpetual spouters, mud pools/volcanoes, and non-eruptive hot springs) and meteoric water. From the isotopic composition we infer that the thermal water has a meteoric origin that is different from the composition of local meteoric water. The absence of detectable tritium in thermal waters indicates that most of the recharge occurred pre-1950. Boiling and steam separation from the deep reservoir appear to be the main subsurface processes affecting the thermal fluids. A large amount of heat is lost to the atmosphere by evaporation from surface water and by steam emitted from erupting geysers. Using the chloride inventory method, we estimate thermal water discharge to be 218 to 234 L/s, and the advective heat flow to be 120 to 170 MW.

El Tatio Basin↗
USGS2019

Seismic and acoustic signatures of surficial mass movements at volcanoes

Surficial mass movements, such as debris avalanches, rock falls, lahars, pyroclastic flows, and outburst floods, are a dominant hazard at many volcanoes worldwide. Understanding these processes, cataloging their spatio-temporal occurrence, and detecting, tracking, and characterizing these events would advance the science of volcano monitoring and help mitigate hazards. Seismic and acoustic methods show promise for achieving these objectives: many surficial mass movements generate observable seismic and acoustic signals, and many volcanoes are already monitored. Significant progress has been made toward understanding, modeling, and extracting quantitative information from seismic and infrasonic signals generated by surficial mass movements. However, much work remains. In this paper, we review the state of the art of the topic, covering a range of scales and event types from individual rock falls to sector collapses. We consider a full variety of volcanic settings, from submarine to subaerial, shield volcano to stratovolcano. Finally, we discuss future directions toward operational seismo-acoustic monitoring of surficial mass movements at volcanoes.

Journal of Volcanology and Geothermal Research↗
USGS2019

Eruption and fountaining dynamics of selected 1985–1986 high fountaining episodes at Kīlauea volcano, Hawai'i, from quantitative vesicle microtexture analysis

Tephra from the early Hawaiian fountaining episodes of the ongoing eruption of Pu'u 'Ō'ō in the East Rift Zone (ERZ) of Kīlauea provides an opportunity to study the vesicle microtextures of pyroclasts erupted from a single vent over a prolonged period of time. We report the results of microtextural analysis of pyroclasts from five of Pu'u 'Ō'ō's high (>200 m) Hawaiian fountaining episodes (episodes 32, 37, 40, 44 and 45) erupted during 1985–1986. This analysis was carried out to constrain the parameters that led to large variations in fountain height at Pu'u 'Ō'o, and the extent to which pyroclast residence times in the fountain modified microtextures. Our results confirm the finding of Stovall et al., 2011 , Stovall et al., 2012 that pyroclasts from a single Hawaiian fountain can vary greatly in texture (from bubbly to foamy), and have vesicle volume densities (N m v ) and vesicle to melt ratios (V G /V L ) that vary by an order of magnitude. This range in vesicle texture and population is due to extensive growth and coalescence of vesicles within the fountain after fragmentation. Only one pyroclast from four of five episodes was found to have textures interpreted as indicative of the vesicle population near the moment of fragmentation: bubbly texture, high density (typically >500 kg m −3 ), high N m v (2.2 × 10 6 to 4.4 × 10 6 ), and low V G /V L of 2.06 to 4.65. We demonstrate a linear correlation between Δ(V G /V L ) and peak fountain height across a range of Hawaiian fountains from Kilauea. This correlation could be used to infer peak heights of unobserved Hawaiian fountaining eruptions after further testing using well-recorded events.

Hawai'i↗
USGS2019

The Miocene Atastra Creek sinter (Bodie Hills volcanic field, California and Nevada): 4D evolution of a geomorphically intact siliceous hot spring deposit

The Atastra Creek siliceous hot spring deposit, or sinter, occurs in the Paramount-Bald Peak alteration zone, due north of the Bodie precious metals mining district in the Miocene Bodie Hills volcanic field, California and Nevada, U.S.A. Distinctive features include its geomorphically intact geyser vent mounds, the presence of growth-fault-stepped sinter terraces that developed westward along a NNW trending fault, and atypical Hg, Sb and As concentrations. High-temperature (>75 °C), subaerial, proximal apron sinter occurs in two settings – (1) radiating botryoidal geyserite in vent mounds and (2) bedded geyserite (columnar, nodular, size-sorted geyser eggs) intercalated with mid- to low-temperature (<60 °C), middle to distal apron terrace sinter. The active geyser vent mounds migrated from south to north across the apron terrace through time. A newly recognized, visually striking type of sinter – termed silica infiltrate herein – formed directly beneath the subaerial geyserite deposits. It is characterized by mainly parallel, thin, irregular beds of white silica with narrow borders of dark blue-gray silica, as well as swirly to globular white silica, all chaotically linked together by silica ‘necks’ and ‘bridges’. The silica infiltrate permeates the pre-existing sinter, and is interpreted to represent super-heated fluid injections into the immediately overlying geyser vent mounds and pools. The adjacent sinter sheets preserve fluid-flow directions of warm-water channels that traversed the discharge apron, including microbial streamer fabric and wavy laminated sinter with flattened bubbles. More distal sinter textures are recorded by rimmed terracettes with mammillated and palisade microbial fabrics and rare, warm pool-related stromatolites. However, plants are absent, probably due to the high metalloids, in particular As. Post-depositional events include an inferred hydrothermal eruption breccia and an interpreted drop in the water table accompanied by incursion of steam acid condensate, as evidenced by pervasive hydrothermal alteration of the Atastra Creek deposits. Hence, the well-exposed sedimentary facies distributions, well-preserved geomorphology, and stratigraphic and structural relationships together allow an unusually complete, paleoenvironmental reconstruction of the Atastra Creek paleohydrothermal deposit.

Journal of Volcanology and Geothermal Research↗
USGS2019

Lava flow morphology at an erupting andesitic stratovolcano: A satellite perspective on El Reventador, Ecuador

Lava flows pose a significant hazard to infrastructure and property located close to volcanoes, and understanding how flows advance is necessary to manage volcanic hazard during eruptions. Compared to low-silica basaltic flows, flows of andesitic composition are infrequently erupted and so relatively few studies of their characteristics and behaviour exist. We use El Reventador, Ecuador as a target to investigate andesitic lava flow properties during a 4.5 year period of extrusive eruption between February 2012 and August 2016. We use satellite radar to map the dimensions of 43 lava flows and look at variations in their emplacement behaviour over time. We find that flows descended the north and south flanks of El Reventador, and were mostly emplaced over durations shorter than the satellite repeat interval of 24 days. Flows ranged in length from 0.3 to 1.7 km, and the length of flows generally decreased over the observation period. We measure a decrease in flow volume with time that is correlated with a long-term exponential decrease in effusion rate, and propose that this behaviour is caused by temporary magma storage in the conduit acting as magma capacitor between the magma reservoir and the surface. We use the dimensions of flow levées and widths to estimate the flow yield strengths. We observe that some flows were diverted by topographic obstacles, and compare measurements of decreased channel width and increased flow thickness at the obstacles with observations from laboratory experiments. Radar observations, such as those presented here, could be used to map and measure properties of evolving lava flow fields at other remote or difficult to monitor volcanoes.

El Reventador↗