Multicomponent geothermometry applied to a medium-low enthalpy carbonate-evaporite geothermal reservoir

Abstract To improve knowledge of the thermal state of medium to low-enthalpy thermal systems hosted in carbonate-evaporite rocks, a mineral-solution equilibrium model was compared to other theoretical geothermometers. We use the GeoT code, which uses as input the chemical composition of water and saturation indices of minerals to calculate water-rock equilibrium over a temperature range of interest. The calculations were applied to the medium and low enthalpy geothermal systems in the Tyrrhenian-Apennine area (central Italy). The lithology consists of a Paleozoic metamorphic basement, overlain by Mesozoic carbonate–evaporite- and Oligocene–Middle Miocene flysch formations, and Quaternary volcanic complexes associated with crustal extension. A regional aquifer is hosted in the carbonate-evaporite formations, and smaller aquifers are hosted in the volcanic rocks. Reservoir temperatures were calculated based on the chemical composition of springs and wells in Central Italy (sampled previously), and in the Cimino-Vicano hydrothermal system (sampled in 2012). Chalcedony and quartz geothermometers provide realistic temperatures. The sensitivity of the model is tested for CO 2 degassing and input minerals. The results of optimized GeoT simulations show that all the samples are affected by degassing during their rise to the surface and that for computing a realistic reservoir temperature it is necessary to consider the principal minerals of the geothermal reservoir (particularly gypsum, quartz, dolomite, aragonite and calcite). The equilibrium temperatures range from 48-115 °C. The statistical approach of "best clustering minerals" solves the problems related to cation or single component geothermometers. Multicomponent geothermometry coupled with optimization provides a reliable approach to reconstruct fluid composition at depth and estimate reservoir temperatures

Triggering and modulation of geyser eruptions in Yellowstone National Park by earthquakes, earth tides, and weather

Author Posting. © American Geophysical Union, 2014. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 119 (2014): 1718–1737, doi:10.1002/2013JB010803.We analyze intervals between eruptions (IBEs) data acquired between 2001 and 2011 at Daisy and Old Faithful geysers in Yellowstone National Park. We focus our statistical analysis on the response of these geysers to stress perturbations from within the solid earth (earthquakes and earth tides) and from weather (air pressure and temperature, precipitation, and wind). We conclude that (1) the IBEs of these geysers are insensitive to periodic stresses induced by solid earth tides and barometric pressure variations; (2) Daisy (pool geyser) IBEs lengthen by evaporation and heat loss in response to large wind storms and cold air; and (3) Old Faithful (cone geyser) IBEs are not modulated by air temperature and pressure variations, wind, and precipitation, suggesting that the subsurface water column is decoupled from the atmosphere. Dynamic stress changes of 0.1−0.2 MPa resulting from the 2002 M-7.9 Denali, Alaska, earthquake surface waves caused a statistically significant shortening of Daisy geyser's IBEs. Stresses induced by other large global earthquakes during the study period were at least an order of magnitude smaller. In contrast, dynamic stresses of >0.5 MPa from three large regional earthquakes in 1959, 1975, and 1983 caused lengthening of Old Faithful's IBEs. We infer that most subannual geyser IBE variability is dominated by internal processes and interaction with other geysers. The results of this study provide quantitative bounds on the sensitivity of hydrothermal systems to external stress perturbations and have implications for studying the triggering and modulation of volcanic eruptions by external forces.K. Luttrell and S. Hurwitz were supported by the USGS Volcano Hazards Program, and Michael Manga was supported by NSF grant EAR1114184.2014-09-0

Triggering and modulation of geyser eruptions in Yellowstone National Park by earthquakes, earth tides, and weather

Author Posting. © American Geophysical Union, 2014. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 119 (2014): 1718–1737, doi:10.1002/2013JB010803.We analyze intervals between eruptions (IBEs) data acquired between 2001 and 2011 at Daisy and Old Faithful geysers in Yellowstone National Park. We focus our statistical analysis on the response of these geysers to stress perturbations from within the solid earth (earthquakes and earth tides) and from weather (air pressure and temperature, precipitation, and wind). We conclude that (1) the IBEs of these geysers are insensitive to periodic stresses induced by solid earth tides and barometric pressure variations; (2) Daisy (pool geyser) IBEs lengthen by evaporation and heat loss in response to large wind storms and cold air; and (3) Old Faithful (cone geyser) IBEs are not modulated by air temperature and pressure variations, wind, and precipitation, suggesting that the subsurface water column is decoupled from the atmosphere. Dynamic stress changes of 0.1−0.2 MPa resulting from the 2002 M-7.9 Denali, Alaska, earthquake surface waves caused a statistically significant shortening of Daisy geyser's IBEs. Stresses induced by other large global earthquakes during the study period were at least an order of magnitude smaller. In contrast, dynamic stresses of >0.5 MPa from three large regional earthquakes in 1959, 1975, and 1983 caused lengthening of Old Faithful's IBEs. We infer that most subannual geyser IBE variability is dominated by internal processes and interaction with other geysers. The results of this study provide quantitative bounds on the sensitivity of hydrothermal systems to external stress perturbations and have implications for studying the triggering and modulation of volcanic eruptions by external forces.K. Luttrell and S. Hurwitz were supported by the USGS Volcano Hazards Program, and Michael Manga was supported by NSF grant EAR1114184.2014-09-0

Heat flow in vapor dominated areas of the Yellowstone Plateau Volcanic Field: Implications for the thermal budget of the Yellowstone Caldera

Characterizing the vigor of magmatic activity in Yellowstone requires knowledge of the mechanisms and rates of heat transport between magma and the ground surface. We present results from a heat flow study in two vapor dominated, acid-sulfate thermal areas in the Yellowstone Caldera, the 0.11 km² Obsidian Pool Thermal Area (OPTA) and the 0.25 km² Solfatara Plateau Thermal Area (SPTA). Conductive heat flux through a low permeability layer capping large vapor reservoirs is calculated from soil temperature measurements at >600 locations and from laboratory measurements of soil properties. The conductive heat output is 3.6 ± 0.4 MW and 7.5 ± 0.4 MW from the OPTA and the SPTA, respectively. The advective heat output from soils is 1.3 ± 0.3 MW and 1.2 ± 0.3 MW from the OPTA and the SPTA, respectively and the heat output from thermal pools in the OPTA is 6.8 ± 1.4 MW. These estimates result in a total heat output of 11.8 ± 1.4 MW and 8.8 ± 0.4 MW from OPTA and SPTA, respectively. Focused zones of high heat flux in both thermal areas are roughly aligned with regional faults suggesting that faults in both areas serve as conduits for the rising acid vapor. Extrapolation of the average heat flux from the OPTA (103 ± 2 W·m⁻²) and SPTA (35 ± 3 W·m⁻²) to the ∼35 km² of vapor dominated areas in Yellowstone yields 3.6 and 1.2 GW, respectively, which is less than the total heat output transported by steam from the Yellowstone Caldera as estimated by the chloride inventory method (4.0 to 8.0 GW)

Eruptions at Lone Star Geyser, Yellowstone National Park, USA: 1. Energetics and eruption dynamics

Author Posting. © American Geophysical Union, 2013. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 118 (2013): 4048–4062, doi:10.1002/jgrb.50251.Geysers provide a natural laboratory to study multiphase eruptive processes. We present results from a 4 day experiment at Lone Star Geyser in Yellowstone National Park, USA. We simultaneously measured water discharge, acoustic emissions, infrared intensity, and visible and infrared video to quantify the energetics and dynamics of eruptions, occurring approximately every 3 h. We define four phases in the eruption cycle (1) a 28±3 min phase with liquid and steam fountaining, with maximum jet velocities of 16–28 m s−1, steam mass fraction of less than ∼0.01. Intermittently choked flow and flow oscillations with periods increasing from 20 to 40 s are coincident with a decrease in jet velocity and an increase of steam fraction; (2) a 26±8 min posteruption relaxation phase with no discharge from the vent, infrared (IR), and acoustic power oscillations gliding between 30 and 40 s; (3) a 59±13 min recharge period during which the geyser is quiescent and progressively refills, and (4) a 69±14 min preplay period characterized by a series of 5–10 min long pulses of steam, small volumes of liquid water discharge, and 50–70 s flow oscillations. The erupted waters ascend from a 160–170°C reservoir, and the volume discharged during the entire eruptive cycle is 20.8±4.1 m3. Assuming isentropic expansion, we calculate a heat output from the geyser of 1.4–1.5 MW, which is <0.1% of the total heat output from Yellowstone Caldera.Support comes from NSF (L. Karlstrom, M. Manga), the USGS Volcano Hazards program (S. Hurwitz, F. Murphy, M.J.S. Johnston, and R.B. McCleskey), and WHOI (R. Sohn).2014-02-1

Eruptions at Lone Star geyser, Yellowstone National Park, USA: 2. Constraints on subsurface dynamics

Author Posting. © American Geophysical Union, 2014. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 119 (2014): 8688–8707, doi:10.1002/2014JB011526.We use seismic, tilt, lidar, thermal, and gravity data from 32 consecutive eruption cycles of Lone Star geyser in Yellowstone National Park to identify key subsurface processes throughout the geyser's eruption cycle. Previously, we described measurements and analyses associated with the geyser's erupting jet dynamics. Here we show that seismicity is dominated by hydrothermal tremor (~5–40 Hz) attributed to the nucleation and/or collapse of vapor bubbles. Water discharge during eruption preplay triggers high-amplitude tremor pulses from a back azimuth aligned with the geyser cone, but during the rest of the eruption cycle it is shifted to the east-northeast. Moreover, ~4 min period ground surface displacements recur every 26 ± 8 min and are uncorrelated with the eruption cycle. Based on these observations, we conclude that (1) the dynamical behavior of the geyser is controlled by the thermo-mechanical coupling between the geyser conduit and a laterally offset reservoir periodically filled with a highly compressible two-phase mixture, (2) liquid and steam slugs periodically ascend into the shallow crust near the geyser system inducing detectable deformation, (3) eruptions occur when the pressure decrease associated with overflow from geyser conduit during preplay triggers an unstable feedback between vapor generation (cavitation) and mass discharge, and (4) flow choking at a constriction in the conduit arrests the runaway process and increases the saturated vapor pressure in the reservoir by a factor of ~10 during eruptions.Funding for USGS team members was provided by the USGS Volcano Hazards Program. R. Sohn's participation was supported by the WHOI Green Technology Program. M. Manga, L. Karlstrom and M. Rudolph did receive salary from the National Science Foundation to spend time on this project.2015-06-0

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, δ18O, 87Sr/86Sr and δ11B isoto- pic 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,HCO3-rich fresh waters (b24.2 °C) flowing in shallow volcanic rocks and SO4-rich thermal waters (25.3 °C to 62.2 °C) ascending from deep, high permeability Mesozoic limestones. The (equivalent) SO4/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 δ18O (−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 (87Sr/86Sr) 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 (b10%) of thermal water as well. The boron isotopic composition (δ11B) of fresh waters (−5.00 and 6.12‰) is similar to that of the volcanic cover, but the δ11B of thermal waters (−8.37‰ to − 4.12‰) is a mismatch for the Mesozoic reservoir rocks and instead reflects dissolution of secondary boron min- erals during fluid ascent through flysch units that overlie the reservoir. A slow and tortuous ascent enhances ex- traction 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

Heat flux from a vapor-dominated hydrothermal field beneath Yellowstone Lake

Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 126(5),(2021): e2020JB021098, https://doi.org/10.1029/2020JB021098.We report results from 149 heat flux measurements made over an ∼2-year interval at sites in and around a vapor-dominated geothermal field located at water depths of ∼100–120 m in Yellowstone Lake, Wyoming. Measurements of both in situ temperature and thermal conductivity as a function of depth were made with a 1 m probe via a remotely operated vehicle, and are combined to compute the vertical conductive heat flux. Inside the ∼55.5 × 103 m2 bathymetric depression demarcating the vapor-dominated field, the median conductive flux is 13 W m−2, with a conductive output of 0.72 MW. Outside the thermal field, the median conductive flux is 3.5 W m−2. We observed 49 active vents inside the thermal field, with an estimated mass discharge rate of 56 kg s−1, a median exit-fluid temperature of 132°C, and a total heat output of 29 MW. We find evidence for relatively weak secondary convection with a total output of 0.09 MW in thermal area lake floor sediments. Our data indicate that vapor beneath the thermal field is trapped by a low-permeability cap at a temperature of ∼189°C and a depth of ∼15 m below the lake floor. The thermal output of the Deep Hole is among the highest of any vapor-dominated field in Yellowstone, due in part to the high boiling temperatures associated with the elevated lake floor pressures.This work was funded by U.S. National Science Foundation (NSF) grants EAR-1515283 to R. N. Harris and J. E. Favorito, EAR-1516361 to R. A. Sohn, and EAR-1514865 to K. M. Luttrell All work in Yellowstone National Park was completed under research permit (YELL-2018-SCI-7018) and the authors thank Annie Carlson from the Yellowstone Center for Resources for logistical help.2021-11-1

Ge/Si and Si isotopes in thermal waters and rivers draining the Yellowstone Plateau Volcanic Field, USA

Yellowstone Plateau Volcanic Field, USA: Ge concentrations measured by inductively coupled plasma mass spectrometry (ICP-MS), Si concentrations measured by inductively coupled plasma optical emission spectrometry (ICP-OES), Ge/Si ratio, Si isotope compositions (d30Si and standard deviation SD) measured by Multicollector ICP-MS, concentrations in Ca, Na, Mg, K measured by ICP-OES, and concentrations in SO4 and Cl measured by ion chromatography in thermal waters, major rivers draining the Yellowstone Plateau Volcanic Field, and creeks flowing into Yellowstone Lake