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 217 records · Page 12Linked to original sources
USGS2015

Magma fracturing and degassing associated with obsidian formation: The explosive–effusive transition

This paper explores the role of melt fracturing in degassing rhyolitic volcanic systems. The Monte Pilato-Rocche Rosse eruptions in Italy evolved from explosive to effusive in style, and H 2 O content in quenched glasses changed over time from relatively H 2 O-rich (~ 0.90 wt.%) to H 2 O-poor dense obsidian (~ 0.10–0.20 wt.%). In addition, healed fractures have been recorded in all different eruptive materials, from the glass of early-erupted tube pumice and rinds of breadcrusted obsidian pyroclasts, to the glass of late-erupted dense obsidian pyroclasts, and throughout the final effusive Rocche Rosse lava flow. These rocks show multiple fault sets, some with crenulated fault planes indicating resumption of viscous flow after faulting, complex obsidian breccias with evidence for post-brecciation folding and stretching, and centimetre- to metre-thick tuffisite preserved in pyroclasts and lava, representing collapsed foam due to fracturing of vesicle walls. These microstructural observations indicate that multiple fracturing and healing events occurred during both explosive and effusive eruptions. H 2 O content in glass decreases by as much as 0.14 wt.% towards healed fractures/faults and decreases in stretched obsidian breccias towards regions of intense brecciation. A drop in pressure and/or increase in temperature along fractures caused diffusive H 2 O migration through melt towards fracture surfaces. Repetitive and pervasive fracturing and healing thereby create conditions for diffusive H 2 O loss into fractures and subsequent escape through permeable paths. This type of progressive magma degassing provides a potential mechanism to explain the formation of dense obsidian and the evolution from explosive to effusive eruption style.

Monte Pilato-Rocche Rosse↗
USGS2015

Dynamics within geyser conduits, and sensitivity to environmental perturbations: insights from a periodic geyser in the El Tatio Geyser Field, Atacama Desert, Chile

Despite more than 200 years of scientific study, the internal dynamics of geyser systems remain poorly characterized. As a consequence, there remain fundamental questions about what processes initiate and terminate eruptions, and where eruptions begin. Over a one-week period in October 2012, we collected down-hole measurements of pressure and temperature in the conduit of an exceptionally regular geyser (132 s/cycle) located in the Chilean desert. We identified four stages in the geyser cycle: (1) recharge of water into the conduit after an eruption, driven by the pressure difference between water in the conduit and in a deeper reservoir; (2) a pre-eruptive stage that follows the recharge and is dominated by addition of steam from below; (3) the eruption, which occurs by rapid boiling of a large mass of water at the top of the water column, and decompression that propagates boiling conditions downward; (4) a relaxation stage during which pressure and temperature decrease until conditions preceding the recharge stage are restored. Eruptions are triggered by the episodic addition of steam coming from depth, suggesting that the dynamics of the eruptions are dominated by geometrical and thermodynamic complexities in the conduit and reservoir. Further evidence favoring the dominance of internal processes in controlling periodicity is also provided by the absence of responses of the geyser to environmental perturbations (air pressure, temperature and probably also Earth tides).

Atacama Desert↗
USGS2015

Magmatic gas emissions at Holocene volcanic features near Mono Lake, California, and their relation to regional magmatism

Silicic lavas have erupted repeatedly in the Mono Basin over the past few thousand years, forming the massive domes and coulees of the Mono Craters chain and the smaller island vents in Mono Lake. We report here on the first systematic study of magmatic CO 2 emissions from these features, conducted during 2007–2010. Most notably, a known locus of weak steam venting on the summit of North Coulee is actually enclosed in a large area (~ 0.25 km 2 ) of diffuse gas discharge that emits 10–14 t/d of CO 2 , mostly at ambient temperature. Subsurface gases sampled here are heavily air-contaminated, but after standard corrections are applied, show average δ 13 C-CO 2 of − 4.72‰, 3 He/ 4 He of 5.89R A , and CO 2 / 3 He of 0.77 × 10 10 , very similar to the values in fumarolic gas from Mammoth Mountain and the Long Valley Caldera immediately to the south of the basin. If these values also characterize the magmatic gas source at Mono Lake, where CO 2 is captured by the alkaline lake water, a magmatic CO 2 upflow beneath the lake of ~ 4 t/d can be inferred. Groundwater discharge from the Mono Craters area transports ~ 13 t/d of 14 C-dead CO 2 as free gas and dissolved carbonate species, and adding in this component brings the estimated total magmatic CO 2 output to 29 t/d for the two silicic systems in the Mono Basin. If these emissions reflect intrusion and degassing of underlying basalt with 0.5 wt.% CO 2 , a modest intrusion rate of 0.00075 km 3 /yr is indicated. Much higher intrusion rates are required to account for CO 2 emissions from Mammoth Mountain and the West Moat of the Long Valley Caldera.

California↗
USGS2015

Coupling geophysical investigation with hydrothermal modeling to constrain the enthalpy classification of a potential geothermal resource.

An appreciable challenge in volcanology and geothermal resource development is to understand the relationships between volcanic systems and low-enthalpy geothermal resources. The enthalpy of an undeveloped geothermal resource in the Karckar region of Armenia is investigated by coupling geophysical and hydrothermal modeling. The results of 3-dimensional inversion of gravity data provide key inputs into a hydrothermal circulation model of the system and associated hot springs, which is used to evaluate possible geothermal system configurations. Hydraulic and thermal properties are specified using maximum a priori estimates. Limited constraints provided by temperature data collected from an existing down-gradient borehole indicate that the geothermal system can most likely be classified as low-enthalpy and liquid dominated. We find the heat source for the system is likely cooling quartz monzonite intrusions in the shallow subsurface and that meteoric recharge in the pull-apart basin circulates to depth, rises along basin-bounding faults and discharges at the hot springs. While other combinations of subsurface properties and geothermal system configurations may fit the temperature distribution equally well, we demonstrate that the low-enthalpy system is reasonably explained based largely on interpretation of surface geophysical data and relatively simple models.

Karckar Region↗
USGS2015

Origins of geothermal gases at Yellowstone

Gas emissions at the Yellowstone Plateau Volcanic Field (YPVF) reflect open-system mixing of gas species originating from diverse rock types, magmas, and crustal fluids, all combined in varying proportions at different thermal areas. Gases are not necessarily in chemical equilibrium with the waters through which they vent, especially in acid sulfate terrain where bubbles stream through stagnant acid water. Gases in adjacent thermal areas often can be differentiated by isotopic and gas ratios, and cannot be tied to one another solely by shallow processes such as boiling-induced fractionation of a parent liquid. Instead, they inherit unique gas ratios (e.g., CH 4 /He) from the dominant rock reservoirs where they originate, some of which underlie the Quaternary volcanic rocks. Steam/gas ratios (essentially H 2 O/CO 2 ) of Yellowstone fumaroles correlate with Ar/He and N 2 /CO 2 , strongly suggesting that H 2 O/CO 2 is controlled by addition of steam boiled from water rich in atmospheric gases. Moreover, H 2 O/CO 2 varies systematically with geographic location, such that boiling is more enhanced in some areas than others. The δ 13 C and 3 He/CO 2 of gases reflect a dominant mantle origin for CO 2 in Yellowstone gas. The mantle signature is most evident at Mud Volcano, which hosts gases with the lowest H 2 O/CO 2 , lowest CH 4 concentrations and highest He isotope ratios (~16Ra), consistent with either a young subsurface intrusion or less input of crustal and meteoric gas than any other location at Yellowstone. Across the YPVF, He isotope ratios ( 3 He/ 4 He) inversely vary with He concentrations, and reflect varied amounts of long- stored, radiogenic He added to the magmatic endmember within the crust. Similarly, addition of CH4 from organic-rich sediments is common in the eastern thermal areas at Yellowstone. Overall, Yellowstone gases reflect addition of deep, high-temperature magmatic gas (CO 2 -rich), lower-temperatures crustal gases ( 4 He- and CH 4 -bearing), and those gases (N 2 , Ne, Ar) added principally through boiling of the meteoric-water-derived geothermal liquid found in the upper few kilometers. We also briefly explore the pathways by which Cl, F, and S, move through the crust.

Montana, Wyoming↗
USGS2015

Temperature profile around a basaltic sill intruded into wet sediments

The transfer of heat into wet sediments from magmatic intrusions or lava flows is not well constrained from field data. Such field constraints on numerical models of heat transfer could significantly improve our understanding of water–lava interactions. We use experimentally calibrated pollen darkening to measure the temperature profile around a basaltic sill emplaced into wet lakebed sediments. It is well known that, upon heating, initially transparent palynomorphs darken progressively through golden, brown, and black shades before being destroyed; however, this approach to measuring temperature has not been applied to volcanological questions. We collected sediment samples from established Miocene fossil localities at Clarkia, Idaho. Fossils in the sediments include pollen from numerous tree and shrub species. We experimentally calibrated changes in the color of Clarkia sediment pollen and used this calibration to determine sediment temperatures around a Miocene basaltic sill emplaced in the sediments. Results indicated a flat temperature profile above and below the sill, with T > 325 °C within 1 cm of the basalt-sediment contact, near 300 °C at 1–2 cm from the contact, and ~ 250 °C at 1 m from the sill contact. This profile suggests that heat transport in the sediments was hydrothermally rather than conductively controlled. This information will be used to test numerical models of heat transfer in wet sediments on Earth and Mars.

Journal of Volcanology and Geothermal Research↗
USGS2015

Heat flux from magmatic hydrothermal systems related to availability of fluid recharge

Magmatic hydrothermal systems are of increasing interest as a renewable energy source. Surface heat flux indicates system resource potential, and can be inferred from soil CO 2 flux measurements and fumarole gas chemistry. Here we compile and reanalyze results from previous CO 2 flux surveys worldwide to compare heat flux from a variety of magma-hydrothermal areas. We infer that availability of water to recharge magmatic hydrothermal systems is correlated with heat flux. Recharge availability is in turn governed by permeability, structure, lithology, rainfall, topography, and perhaps unsurprisingly, proximity to a large supply of water such as the ocean. The relationship between recharge and heat flux interpreted by this study is consistent with recent numerical modeling that relates hydrothermal system heat output to rainfall catchment area. This result highlights the importance of recharge as a consideration when evaluating hydrothermal systems for electricity generation, and the utility of CO 2 flux as a resource evaluation tool.

Journal of Volcanology and Geothermal Research↗
USGS2015

Magma extrusion during the Ubinas 2013–2014 eruptive crisis based on satellite thermal imaging (MIROVA) and ground-based monitoring

After 3 years of mild gases emissions, the Ubinas volcano entered in a new eruptive phase on September 2nd, 2013. The MIROVA system (a space-based volcanic hot-spot detection system), allowed us to detect in near real time the thermal emissions associated with the eruption and provided early evidence of magma extrusion within the deep summit crater. By combining IR data with plume height, sulfur emissions, hot spring temperatures and seismic activity, we interpret the thermal output detected over Ubinas in terms of extrusion rates associated to the eruption. We suggest that the 2013–2014 eruptive crisis can be subdivided into three main phases: ( i ) shallow magma intrusion inside the edifice, (ii) extrusion and growing of a lava plug at the bottom of the summit crater coupled with increasing explosive activity and finally, (iii) disruption of the lava plug and gradual decline of the explosive activity. The occurrence of the 8.2 Mw Iquique (Chile) earthquake (365 km away from Ubinas) on April 1st, 2014, may have perturbed most of the analyzed parameters, suggesting a prompt interaction with the ongoing volcanic activity. In particular, the analysis of thermal and seismic datasets shows that the earthquake may have promoted the most intense thermal and explosive phase that culminated in a major explosion on April 19th, 2014. These results reveal the efficiency of space-based thermal observations in detecting the extrusion of hot magma within deep volcanic craters and in tracking its evolution. We emphasize that, in combination with other geophysical and geochemical datasets, MIROVA is an essential tool for monitoring remote volcanoes with rather difficult accessibility, like those of the Andes that reach remarkably high altitudes.

Journal of Volcanology and Geothermal Research↗
USGS2016

Volcano-tectonic earthquakes: A new tool for estimating intrusive volumes and forecasting eruptions

We present data on 136 high-frequency earthquakes and swarms, termed volcano-tectonic (VT) seismicity, which preceded 111 eruptions at 83 volcanoes, plus data on VT swarms that preceded intrusions at 21 other volcanoes. We find that VT seismicity is usually the earliest reported seismic precursor for eruptions at volcanoes that have been dormant for decades or more, and precedes eruptions of all magma types from basaltic to rhyolitic and all explosivities from VEI 0 to ultraplinian VEI 6 at such previously long-dormant volcanoes. Because large eruptions occur most commonly during resumption of activity at long-dormant volcanoes, VT seismicity is an important precursor for the Earth's most dangerous eruptions. VT seismicity precedes all explosive eruptions of VEI ≥ 5 and most if not all VEI 4 eruptions in our data set. Surprisingly we find that the VT seismicity originates at distal locations on tectonic fault structures at distances of one or two to tens of kilometers laterally from the site of the eventual eruption, and rarely if ever starts beneath the eruption site itself. The distal VT swarms generally occur at depths almost equal to the horizontal distance of the swarm from the summit out to about 15 km distance, beyond which hypocenter depths level out. We summarize several important characteristics of this distal VT seismicity including: swarm-like nature, onset days to years prior to the beginning of magmatic eruptions, peaking of activity at the time of the initial eruption whether phreatic or magmatic, and large non-double couple component to focal mechanisms. Most importantly we show that the intruded magma volume can be simply estimated from the cumulative seismic moment of the VT seismicity from: Log10 V = 0.77 Log ΣMoment − 5.32, with volume, V, in cubic meters and seismic moment in Newton meters. Because the cumulative seismic moment can be approximated from the size of just the few largest events, and is quite insensitive to precise locations, the intruded magma volume can be quickly and easily estimated with few short-period seismic stations. Notable cases in which distal VT events preceded eruptions at long-dormant volcanoes include: Nevado del Ruiz (1984–1985), Pinatubo (1991), Unzen (1989–1995), Soufriere Hills (1995), Shishaldin (1989–1999), Tacana' (1985–1986), Pacaya (1980–1984), Rabaul (1994), and Cotopaxi (2001). Additional cases are recognized at frequently active volcanoes including Popocateptl (2001–2003) and Mauna Loa (1984). We present four case studies (Pinatubo, Soufriere Hills, Unzen, and Tacana') in which we demonstrate the above mentioned VT characteristics prior to eruptions. Using regional data recorded by NEIC, we recognized in near-real time that a huge distal VT swarm was occurring, deduced that a proportionately huge magmatic intrusion was taking place beneath the long dormant Sulu Range, New Britain Island, Papua New Guinea, that it was likely to lead to eruptive activity, and warned Rabaul Volcano Observatory days before a phreatic eruption occurred. This confirms the value of this technique for eruption forecasting. We also present a counter-example where we deduced that a VT swarm at Volcan Cosiguina, Nicaragua, indicated a small intrusion, insufficient to reach the surface and erupt. Finally, we discuss limitations of the method and propose a mechanism by which this distal VT seismicity is triggered by magmatic intrusion.

Journal of Volcanology and Geothermal Research↗
USGS2016

Radiocarbon dating of silica sinter deposits in shallow drill cores from the Upper Geyser Basin, Yellowstone National Park

To explore the timing of hydrothermal activity at the Upper Geyser Basin (UGB) in Yellowstone National Park, we obtained seven new accelerator mass spectrometry (AMS) radiocarbon 14 C ages of carbonaceous material trapped within siliceous sinter. Five samples came from depths of 15–152 cm within the Y-1 well, and two samples were from well Y-7 (depths of 24 cm and 122 cm). These two wells, at Black Sand and Biscuit Basins, respectively, were drilled in 1967 as part of a scientific drilling program by the U.S. Geological Survey ( White et al., 1975 ). Even with samples as small as 15 g, we obtained sufficient carbonaceous material (a mixture of thermophilic mats, pollen, and charcoal) for the 14 C analyses. Apparent time of deposition ranged from 3775 ± 25 and 2910 ± 30 14 C years BP at the top of the cores to about 8000 years BP at the bottom. The dates are consistent with variable rates of sinter formation at individual sites within the UGB over the Holocene. On a basin-wide scale, though, these and other existing 14 C dates hint that hydrothermal activity at the UGB may have been continuous throughout the Holocene.

Wyoming↗
USGS2016

Timing and composition of continental volcanism at Harrat Hutaymah, western Saudi Arabia

Harrat Hutaymah is an alkali basalt volcanic field in north-central Saudi Arabia, at the eastern margin of a large Neogene continental, intraplate magmatic province. Lava flow, tephra and spatter cone compositions in the field include alkali olivine basalts and basanites. These compositions contrast with the predominantly tholeiitic, fissure-fed basalts found along the eastern margin of the Red Sea. The Hutaymah lava flows were erupted through Proterozoic arc-associated plutonic and meta-sedimentary rocks of the Arabian shield, and commonly contain a range of sub-continental lithospheric xenoliths, although the lavas themselves show little indication of crustal contamination. Previous radiometric dating of this volcanic field (a single published K–Ar age; 1.8 Ma) is suspiciously old given the field measurement of normal magnetic polarity only (i.e. Brunhes interval, ≤ 780 Ka). We report new age determinations on 14 lava flows by the 40Ar–39Ar laser step heating method, all younger than ~ 850 Ka, to better constrain the time frame of volcanism, and major, trace and rare earth element compositions to describe the chemical variation of volcanic activity at Harrat Hutaymah. Crystal fractionation was dominated by olivine ± clinopyroxene at a range of upper mantle and crustal pressures. Rapid ascent and eruption of magma is indicated by the array of lower crustal and lithospheric xenoliths observed in lava flows and tephra. Modeling suggests 1–7% melting of an enriched asthenospheric mantle source occurred beneath Harrat Hutaymah under a relatively thick lithospheric cap (60–80 km).

Harrat Hutaymah↗
USGS2016

Results of the eruptive column model inter-comparison study

This study compares and evaluates one-dimensional (1D) and three-dimensional (3D) numerical models of volcanic eruption columns in a set of different inter-comparison exercises. The exercises were designed as a blind test in which a set of common input parameters was given for two reference eruptions, representing a strong and a weak eruption column under different meteorological conditions. Comparing the results of the different models allows us to evaluate their capabilities and target areas for future improvement. Despite their different formulations, the 1D and 3D models provide reasonably consistent predictions of some of the key global descriptors of the volcanic plumes. Variability in plume height, estimated from the standard deviation of model predictions, is within ~ 20% for the weak plume and ~ 10% for the strong plume. Predictions of neutral buoyancy level are also in reasonably good agreement among the different models, with a standard deviation ranging from 9 to 19% (the latter for the weak plume in a windy atmosphere). Overall, these discrepancies are in the range of observational uncertainty of column height. However, there are important differences amongst models in terms of local properties along the plume axis, particularly for the strong plume. Our analysis suggests that the simplified treatment of entrainment in 1D models is adequate to resolve the general behaviour of the weak plume. However, it is inadequate to capture complex features of the strong plume, such as large vortices, partial column collapse, or gravitational fountaining that strongly enhance entrainment in the lower atmosphere. We conclude that there is a need to more accurately quantify entrainment rates, improve the representation of plume radius, and incorporate the effects of column instability in future versions of 1D volcanic plume models.

Journal of Volcanology and Geothermal Research↗
USGS2016

Inter-comparison of three-dimensional models of volcanic plumes

We performed an inter-comparison study of three-dimensional models of volcanic plumes. A set of common volcanological input parameters and meteorological conditions were provided for two kinds of eruptions, representing a weak and a strong eruption column. From the different models, we compared the maximum plume height, neutral buoyancy level (where plume density equals that of the atmosphere), and level of maximum radial spreading of the umbrella cloud. We also compared the vertical profiles of eruption column properties, integrated across cross-sections of the plume (integral variables). Although the models use different numerical procedures and treatments of subgrid turbulence and particle dynamics, the inter-comparison shows qualitatively consistent results. In the weak plume case (mass eruption rate 1.5 × 10 6 kg s − 1 ), the vertical profiles of plume properties (e.g., vertical velocity, temperature) are similar among models, especially in the buoyant plume region. Variability among the simulated maximum heights is ~ 20%, whereas neutral buoyancy level and level of maximum radial spreading vary by ~ 10%. Time-averaging of the three-dimensional (3D) flow fields indicates an effective entrainment coefficient around 0.1 in the buoyant plume region, with much lower values in the jet region, which is consistent with findings of small-scale laboratory experiments. On the other hand, the strong plume case (mass eruption rate 1.5 × 10 9 kg s − 1 ) shows greater variability in the vertical plume profiles predicted by the different models. Our analysis suggests that the unstable flow dynamics in the strong plume enhances differences in the formulation and numerical solution of the models. This is especially evident in the overshooting top of the plume, which extends a significant portion (~ 1/8) of the maximum plume height. Nonetheless, overall variability in the spreading level and neutral buoyancy level is ~ 20%, whereas that of maximum height is ~ 10%. This inter-comparison study has highlighted the different capabilities of 3D volcanic plume models, and identified key features of weak and strong plumes, including the roles of jet stability, entrainment efficiency, and particle non-equilibrium, which deserve future investigation in field, laboratory, and numerical studies.

Journal of Volcanology and Geothermal Research↗
USGS2016

Mihi Breccia: A stack of lacustrine sediments and subaqueous pyroclastic flows within the Taupo Volcanic Zone, New Zealand

The Taupo Volcanic Zone (TVZ), New Zealand, encompasses a wide variety of arc-related strata, although most of its small-volume (non-caldera-forming) eruptions are poorly-exposed and extensively hydrothermally altered. The Mihi Breccia is a stratigraphic sequence consisting of interbedded rhyolitic pyroclastic flows and lacustrine sediments with eruption ages of 281 ± 18 to at least 239 ± 6 ka (uncertainties at 2σ). In contrast to other small-volume rhyolitic eruptions within the TVZ, Mihi Breccia is relatively well-exposed within the Paeroa fault block, and contains minimal hydrothermal alteration. Pyroclastic flow characteristics and textures including: 1) breadcrusted juvenile clasts, 2) lack of welding, 3) abundant ash-rich matrix, 4) lack of fiamme and eutaxitic textures, 5) lack of thermal oxidation colors, 6) lack of cooling joints, 7) exclusive lacustrine sediment lithic clasts, and 8) interbedding with lacustrine sediments, all indicating that Mihi Breccia strata originated in a paleo-lake system. This ephemeral paleo-lake system is inferred to have lasted for > 50 kyr (based on Mihi Breccia age constraints), and referred to as Huka Lake. Mihi Breccia pyroclastic flow juvenile clast geochemistry and petrography correspond with similar-aged (264 ± 8, 263 ± 10, and 247 ± 4 ka) intra-caldera rhyolite domes filling the Reporoa caldera (source of the 281 ± 81 Kaingaroa Formation ignimbrite). These exposed intra-caldera rhyolite domes (as well as geophysically inferred subsurface domes) are proposed to be source vents for the Mihi Breccia pyroclastic flows. Soft-sediment deformation associated with Mihi Breccia strata indicate either seismic shock, rapid sediment loading during pyroclastic flow emplacement, or both. Thus, the Mihi Breccia reflects a prolonged series of subaqueous rhyolite dome building and associated pyroclastic flows, accompanied by seismic activity, emplaced into a large paleo-lake system within the TVZ.

Taupo Volcanic Zone↗
USGS2016

Quantification of the intrusion process at Kīlauea volcano, Hawai'i

The characteristic size of two types of intrusions identified beneath Kīlauea's East Rift zone are uniquely estimated by combining time constraints from fractional crystallization and the rates of magma solidification during cooling. Some intrusions were rapidly emplaced as dikes, but stalled before reaching the surface, and cooled and crystallized to feed later fractionated eruptions. More specifically, using the observed time interval between initial emplacement and eruption of fractionated lava, whose degree of fractionation is estimated from petrologic mixing calculations, the extent of solidification or cooling needed to produce this amount of fractionation can be directly inferred. And from the known erupted volumes the spatial extent or size of this fractionated volume can be analytically related to the full size of the source body itself. Two examples yield dike widths of 82 and 68 m. Other intrusions remain close to the east rift magma transport path and are observed to last for decades or longer as viable magma bodies that may participate in feeding later eruptions. The thickness of semi-permanent reservoirs near the East Rift Zone magma transport path can be estimated by assuming a resupply rate that is sufficiently frequent to restrict cooling to < 10 °C. It is inferred that both types of intrusions likely began as dike offshoots from the East Rift Zone magma transport path, but the frequently resupplied bodies may have later been converted to sills or laccoliths of heights estimated at 43–62 m. Our modeled intrusions contrast with models of rapidly emplaced thinner dikes feeding shallow intrusions, which are accompanied by intense rift earthquake swarms and are often associated with eruptions. These calculations show that long-term heating of the wallrock of the magma transport paths serves to slow conduit cooling, which may be partly responsible for sustaining long East Rift Zone eruptions. Adjacent to the vertical transport path beneath Kīlauea's summit, the combined effects of heating and ever-increasing magma supply rate may have forced a commensurate enlarging of the conduit, perhaps explaining the occurrence of a temporary burst of deep (5–15 km) long-period earthquake swarms between 1987 and 1992.

Hawaii↗
USGS2016

Multireaction equilibrium geothermometry: A sensitivity analysis using data from the Lower Geyser Basin, Yellowstone National Park, USA

A multireaction chemical equilibria geothermometry (MEG) model applicable to high-temperature geothermal systems has been developed over the past three decades. Given sufficient data, this model provides more constraint on calculated reservoir temperatures than classical chemical geothermometers that are based on either the concentration of silica (SiO 2 ), or the ratios of cation concentrations. A set of 23 chemical analyses from Ojo Caliente Spring and 22 analyses from other thermal features in the Lower Geyser Basin of Yellowstone National Park are used to examine the sensitivity of calculated reservoir temperatures using the GeoT MEG code (Spycher et al. 2013, 2014) to quantify the effects of solute concentrations, degassing, and mineral assemblages on calculated reservoir temperatures. Results of our analysis demonstrate that the MEG model can resolve reservoir temperatures within approximately ±15°C, and that natural variation in fluid compositions represents a greater source of variance in calculated reservoir temperatures than variations caused by analytical uncertainty (assuming ~5% for major elements). The analysis also suggests that MEG calculations are particularly sensitive to variations in silica concentration, the concentrations of the redox species Fe(II) and H 2 S, and that the parameters defining steam separation and CO 2 degassing from the liquid may be adequately determined by numerical optimization. Results from this study can provide guidance for future applications of MEG models, and thus provide more reliable information on geothermal energy resources during exploration.

Idaho, Montana, Wyoming↗
USGS2016

Juvenile magma recognition and eruptive dynamics inferred from the analysis of ash time series: The 2015 reawakening of Cotopaxi volcano

Forecasting future activity and performing hazard assessments during the reactivation of volcanoes remain great challenges for the volcanological community. On August 14, 2015 Cotopaxi volcano erupted for the first time in 73 years after approximately four months of precursory activity, which included an increase in seismicity, gas emissions, and minor ground deformation. Here we discuss the use of near real-time petrological monitoring of ash samples as a complementary aid to geophysical monitoring, in order to infer eruption dynamics and evaluate possible future eruptive activity at Cotopaxi. Twenty ash samples were collected between August 14 and November 23, 2015 from a monitoring site on the west flank of the volcano. These samples contain a range of grain types that we classified as: hydrothermal/altered, lithic, juvenile, and free crystals. The relative proportions of theses grains evolved as the eruption progressed, with increasing amounts of juvenile material and a decrease in hydrothermally altered material. In samples from the initial explosion, juvenile grains are glassy, microlite-poor and contain hydrothermal minerals (opal and alunite). The rising magma came in contact with the hydrothermal system under confinement, causing hydro-magmatic explosions that cleared the upper part of the plumbing system. Subsequently, the magmatic column produced a thermal aureole in the conduit and dried out the hydrothermal system, allowing for dry eruptions. Magma ascent rates were low enough to allow for efficient outgassing and microlite growth. Constant supply of magma from below caused quasi-continuous disruption of the uppermost magma volume through a combination of shear-deformation and gas expansion. The combination of increasing crystallinity of juvenile grains, and high measured SO 2 flux indicate decreasing integrated magma ascent rates and clearing of the hydrothermal system along transport pathways in a system open to gas loss. The near real-time monitoring of ash samples combined with traditional geophysical monitoring techniques during the reawakening of Cotopaxi allowed us to gain a much clearer understanding of events than when using traditional geophysical monitoring alone.

Cotopaxi volcano↗
USGS2016

The chemistry and isotopic composition of waters in the low-enthalpy geothermal system of Cimino-Vico Volcanic District, Italy

Geothermal energy exploration is based in part on interpretation of the chemistry, temperature, and discharge rate of thermal springs. Here we present the major element chemistry and the δD, δ 18 O, 87 Sr/ 86 Sr and δ 11 B isotopic ratio of groundwater from the low-enthalpy geothermal system near the city of Viterbo in the Cimino-Vico volcanic district of west-Central Italy. The geothermal system hosts many thermal springs and gas vents, but the resource is still unexploited. Water chemistry is controlled by mixing between low salinity,HCO 3 -rich fresh waters (<24.2°C) flowing in shallow volcanic rocks and SO 4 -rich thermal waters (25.3°C to 62.2°C) ascending from deep, high permeability Mesozoic limestones. The (equivalent) SO 4 /Cl (0.01–0.02), Na/Cl (2.82–5.83) and B/Cl ratios (0.02–0.38) of thermal waters differs from the ratios in other geothermal systems from Central Italy, probably implying a lack of hydraulic continuity across the region. The δ 18 O (−6.6‰ to −5.9‰) and δD (−40.60‰ to −36.30‰) isotopic composition of spring water suggest that the recharge area for the geothermal system is the summit region of Mount Cimino. The strontium isotope ratios ( 87 Sr/ 86 Sr) of thermal waters (0.70797–0.70805) are consistent with dissolution of the Mesozoic evaporite-carbonate units that constitute the reservoir, and the ratios of cold fresh waters mainly reflect shallow circulation through the volcanic cover and some minor admixture (<10%) of thermal water as well. The boron isotopic composition (δ 11 B) of fresh waters (−5.00 and 6.12‰) is similar to that of the volcanic cover, but the δ 11 B of thermal waters (−8.37‰ to −4.12‰) is a mismatch for the Mesozoic reservoir rocks and instead reflects dissolution of secondary boron minerals during fluid ascent through flysch units that overlie the reservoir. A slow and tortuous ascent enhances extraction of boron but also promotes conductive cooling, partially masking the heat present in the reservoir. Overall data from this study is consistent with previous studies that concluded that the geothermal system has a large energy potential.

Journal of Volcanology and Geothermal Research↗