Geology topics

Shaul Hurwitz

Publications and source records attributed to Shaul Hurwitz.

At least 55 records · Page 3Linked to original sources

USGS2018

Combining multiphase groundwater flow and slope stability models to assess stratovolcano flank collapse in the Cascade Range

Hydrothermal alteration can create low‐permeability zones, potentially resulting in elevated pore‐fluid pressures, within a volcanic edifice. Strength reduction by rock alteration and high pore‐fluid pressures have been suggested as a mechanism for edifice flank instability. Here we combine numerical models of multiphase heat transport and groundwater flow with a slope‐stability code that incorporates three‐dimensional distributions of strength and pore‐water pressure to address the following questions: (1) What permeability distributions and contrasts produce elevated pore‐fluid pressures in a stratovolcano? (2) What are the effects of these elevated pressures on flank stability? (3) Finally, what are the effects of magma intrusion on potential flank failure in an edifice? Simulation results show that under a range of plausible parameters, water tables in a stratovolcano can be elevated or perched. These elevated water tables result in universally lower stability (lower factor of safety) compared with equivalent dry edifices, indicating a higher likelihood of flank collapse. Low‐permeability (<1 × 10 −17 m 2 ) layers such as altered pyroclastic deposits or breccias can result in locally saturated regions (perched water) and lower factors of safety near the ground surface but may actually reduce liquid water saturation and pore pressures in the core of the edifice and thus may favor small, shallow collapses over larger, deeper collapses. Magma intrusion into the base of the edifice increases pore‐fluid pressures and decreases the factor of safety. However, the shear strength of edifice rocks also exerts a significant control on stability, so both mechanical properties and pore‐fluid pressures are important for stability assessments.

Washington↗

USGS2018

Hot water in the Long Valley Caldera—The benefits and hazards of this large natural resource

The volcanic processes that have shaped the Long Valley Caldera in eastern California have also created an abundant supply of natural hot water. This natural resource provides benefits to many users, including power generation at the Casa Diablo Geothermal Plant, warm water for a state fish hatchery, and beautiful scenic areas such as Hot Creek gorge for visitors. However, some features can be dangerous because of sudden and unpredictable changes in the location and flow rate of boiling water. The U.S. Geological Survey monitors several aspects of the hydrothermal system in the Long Valley Caldera including temperature, flow rate, and water chemistry.

California↗

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↗

USGS2017

Probing magma reservoirs to improve volcano forecasts

When it comes to forecasting eruptions, volcano observatories rely mostly on real-time signals from earthquakes, ground deformation, and gas discharge, combined with probabilistic assessments based on past behavior [Sparks and Cashman, 2017] . There is comparatively less reliance on geophysical and petrological understanding of subsurface magma reservoirs.

Eos, Earth and Space Science News↗

USGS2017

The fascinating and complex dynamics of geyser eruptions

Geysers episodically erupt liquid and vapor. Despite two centuries of scientific study, basic questions persist—why do geysers exist? What determines eruption intervals, durations, and heights? What initiates eruptions? Through monitoring eruption intervals, analyzing geophysical data, taking measurements within geyser conduits, performing numerical simulations, and constructing laboratory models, some of these questions have been addressed. Geysers are uncommon because they require a combination of abundant water recharge, magmatism, and rhyolite flows to supply heat and silica, and large fractures and cavities overlain by low-permeability materials to trap rising multiphase and multicomponent fluids. Eruptions are driven by the conversion of thermal to kinetic energy during decompression. Larger and deeper cavities permit larger eruptions and promote regularity by isolating water from weather variations. The ejection velocity may be limited by the speed of sound of the liquid + vapor mixture.

Annual Review of Earth and Planetary Sciences↗

USGS2016

An experimental study of the role of subsurface plumbing on geothermal discharge

In order to better understand the diverse discharge styles and eruption intervals observed at geothermal features, we performed three series of laboratory experiments with differing plumbing geometries. A single, straight conduit that connects a hot water bath (flask) to a vent (funnel) can originate geyser-like periodic eruptions, continuous discharge like a boiling spring, and fumarole-like steam discharge, depending on the conduit length and radius. The balance between the heat loss from the conduit walls and the heat supplied from the bottom determines whether and where water can condense which in turn controls discharge style. Next, we connected the conduit to a cold water reservoir through a branch, simulating the inflow from an external water source. Colder water located at a higher elevation than a branching point can flow into the conduit to stop the boiling in the flask, controlling the periodicity of the eruption. When an additional branch is connected to a second cold water reservoir, the two cold reservoirs can interact. Our experiments show that branching allows new processes to occur, such as recharge of colder water and escape of steam from side channels, leading to greater variation in discharge styles and eruption intervals. This model is consistent with the fact that eruption duration is not controlled by emptying reservoirs. We show how differences in plumbing geometries can explain various discharge styles and eruption intervals observed in El Tatio, Chile, and Yellowstone, USA.

Geochemistry, Geophysics, Geosystems↗

USGS2016

Dissolved gases in hydrothermal (phreatic) and geyser eruptions at Yellowstone National Park, USA

Multiphase and multicomponent fluid flow in the shallow continental crust plays a significant role in a variety of processes over a broad range of temperatures and pressures. The presence of dissolved gases in aqueous fluids reduces the liquid stability field toward lower temperatures and enhances the explosivity potential with respect to pure water. Therefore, in areas where magma is actively degassing into a hydrothermal system, gas-rich aqueous fluids can exert a major control on geothermal energy production, can be propellants in hazardous hydrothermal (phreatic) eruptions, and can modulate the dynamics of geyser eruptions. We collected pressurized samples of thermal water that preserved dissolved gases in conjunction with precise temperature measurements with depth in research well Y-7 (maximum depth of 70.1 m; casing to 31 m) and five thermal pools (maximum depth of 11.3 m) in the Upper Geyser Basin of Yellowstone National Park, USA. Based on the dissolved gas concentrations, we demonstrate that CO 2 mainly derived from magma and N 2 from air-saturated meteoric water reduce the near-surface saturation temperature, consistent with some previous observations in geyser conduits. Thermodynamic calculations suggest that the dissolved CO 2 and N 2 modulate the dynamics of geyser eruptions and are likely triggers of hydrothermal eruptions when recharged into shallow reservoirs at high concentrations. Therefore, monitoring changes in gas emission rate and composition in areas with neutral and alkaline chlorine thermal features could provide important information on the natural resources (geysers) and hazards (eruptions) in these areas.

Idaho, Montana, Wyoming↗

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

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

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↗

USGS2015

Tracing chlorine sources of thermal and mineral springs along and across the Cascade Range using halogen and chlorine isotope compositions

In order to provide constraints on the sources of chlorine in spring waters associated with arc volcanism, the major/minor element concentrations and stable isotope compositions of chlorine, oxygen, and hydrogen were measured in 28 thermal and mineral springs along the Cascade Range in northwestern USA. Chloride concentrations in the springs range from 64 to 19,000 mg/L and View the MathML source values range from +0.2&permil; to +1.9&permil; (average=+1.0±0.4&permil;), with no systematic variation along or across the arc, nor correlations with their presumed underlying basement lithologies. Additionally, nine geochemically well-characterized lavas from across the Mt. St. Helens/Mt. Adams region of the Cascade Range (Leeman et al., 2004 and Leeman et al., 2005) were analyzed for their halogen concentrations and Cl isotope compositions. In the arc lavas, Cl and Br concentrations from the volcanic front are higher than in lavas from the forearc and backarc. F and I concentrations progressively decrease from forearc to backarc, similar to the trend documented for B in most arcs. View the MathML source values of the lavas range from −0.1 to +0.8&permil; (average = +0.4±0.3&permil;). Our results suggest that the predominantly positive View the MathML source values observed in the springs are consistent with water interaction with underlying 37Cl-enriched basalt and/or altered oceanic crust, thereby making thermal spring waters a reasonable proxy for the Cl isotope compositions of associated volcanic rocks in the Cascades. However, waters with View the MathML source values >+1.0&permil; also suggest additional contributions of chlorine degassed from cooling magmas due to subsurface vapor–liquid HCl fractionation in which Cl is lost to the aqueous fluid phase and 37Cl is concentrated in the ascending magmatic HCl vapor. Future work is necessary to better constrain Cl isotope behavior during volcanic degassing and fluid–rock interaction in order to improve volatile flux estimates through subduction zones.

Cascade Range↗

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↗

USGS2014

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

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.

Wyoming↗

USGS2014

Dynamics of the Yellowstone hydrothermal system

The Yellowstone Plateau Volcanic Field is characterized by extensive seismicity, episodes of uplift and subsidence, and a hydrothermal system that comprises more than 10,000 thermal features, including geysers, fumaroles, mud pots, thermal springs, and hydrothermal explosion craters. The diverse chemical and isotopic compositions of waters and gases derive from mantle, crustal, and meteoric sources and extensive water-gas-rock interaction at variable pressures and temperatures. The thermal features are host to all domains of life that utilize diverse inorganic sources of energy for metabolism. The unique and exceptional features of the hydrothermal system have attracted numerous researchers to Yellowstone beginning with the Washburn and Hayden expeditions in the 1870s. Since a seminal review published a quarter of a century ago, research in many fields has greatly advanced our understanding of the many coupled processes operating in and on the hydrothermal system. Specific advances include more refined geophysical images of the magmatic system, better constraints on the time scale of magmatic processes, characterization of fluid sources and water-rock interactions, quantitative estimates of heat and magmatic volatile fluxes, discovering and quantifying the role of thermophile microorganisms in the geochemical cycle, defining the chronology of hydrothermal explosions and their relation to glacial cycles, defining possible links between hydrothermal activity, deformation, and seismicity; quantifying geyser dynamics; and the discovery of extensive hydrothermal activity in Yellowstone Lake. Discussion of these many advances forms the basis of this review.

Idaho, Montana, Wyoming↗

USGS2014

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

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.

Wyoming↗

USGS2013

Eruptions at Lone Star Geyser, Yellowstone National Park, USA, part 1: energetics and eruption dynamics

Geysers provide a natural laboratory to study multiphase eruptive processes. We present results from a four–day experiment at Lone Star Geyser in Yellowstone National Park, USA. We simultaneously measured water discharge, acoustic emissions, infraredintensity, and visible and infrared video to quantify the energetics and dynamics of eruptions, occurring approximately every three hours. We define four phases in the eruption cycle: 1) a 28 ± 3 minute phase with liquid and steam fountaining, with maximum jet velocities of 16–28 m s − 1 , steam mass fraction of less than &sim; 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 minute post–eruption relaxation phase with no discharge from the vent, infrared (IR) and acoustic power oscillations gliding between 30 and 40 s; 3) a 59 ± 13 minute recharge period during which the geyser is quiescent and progressively refills, and 4) a 69 ± 14 minute pre–play period characterized by a series of 5–10 minute–long pulses of steam, small volumes of liquid water discharge and 50–70 s flow oscillations. The erupted waters ascend froma 160 − 170° C reservoir and the volume discharged during the entire eruptive cycle is 20.8 ± 4.1 m 3 . 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.

Wyoming↗

USGS2013

Constraints on the upper crustal magma reservoir beneath Yellowstone Caldera inferred from lake-seiche induced strain observations

Seiche waves in Yellowstone Lake with a ~78-minute period and heights <10 cm act as a load on the solid earth observed by borehole strainmeters with subnanostrain sensitivity throughout the Yellowstone Caldera. The far-field strain induced by the load of the seiche waves calculated with a homogeneous elastic model representing the upper crust is more than an order of magnitude smaller than the measured strain amplitude ~30 km from the lake shore. By contrast, the observed far field strain amplitudes are consistent with the seiche load on a two-layered viscoelastic model representing an elastic upper crust overlying a partially molten body deeper than 3–6 km with Maxwell viscosity less than 10 11 Pa s. These strain observations and models provide independent evidence for the presence of partially molten material in the upper crust, consistent with seismic tomography studies that inferred 10%–30% melt fraction in the upper crust.

Wyoming↗

USGS2013

Ca, Sr, O and D isotope approach to defining the chemical evolution of hydrothermal fluids: example from Long Valley, CA, USA

We present chemical and isotopic data for fluids, minerals and rocks from the Long Valley meteoric-hydrothermal system. The samples encompass the presumed hydrothermal upwelling zone in the west moat of the caldera, the Casa Diablo geothermal field, and a series of wells defining a nearly linear, ∼16 km long, west-to-east trend along the likely fluid flow path. Fluid samples were analyzed for the isotopes of water, Sr, and Ca, the concentrations of major cations and anions, alkalinity, and total CO 2 . Water isotope data conform to trends documented in earlier studies, interpreted as indicating a single hydrothermal fluid mixing with local groundwater. Sr isotopes show subtle changes along the flow path, which requires rapid fluid flow and minimal reaction between the channelized fluids and the wallrocks. Sr and O isotopes are used to calculate fracture spacing using a dual porosity model. Calculated fracture spacing and temperature data for hydrothermal fluids indicate the system is (approximately) at steady-state. Correlated variations among total CO 2 , and the concentration and isotopic composition of Ca suggest progressive fluid degassing (loss of CO 2 ), which drives calcite precipitation as the fluid flows west-to-east and cools. The shifts in Ca isotopes require that calcite precipitated at temperatures of 150–180 °C is fractionated by ca. −0.3‰ to −0.5‰ relative to aqueous species. Our data are the first evidence that Ca isotopes undergo kinetic fractionation at high temperatures (>100 °C) and can be used to trace calcite precipitation along hydrothermal fluid flow paths.

California↗