Shaul Hurwitz

Hydrothermal explosions are a significant geological hazard in some active volcanic systems; however, the timing and triggering mechanisms of these explosions are poorly constrained. This study applies luminescence dating techniques to hydrothermal explosion deposits in the Yellowstone Plateau volcanic field to constrain explosion chronologies and evaluate potential triggering mechanisms. We tested four luminescence dating techniques: K-feldspar post-infrared infrared stimulated luminescence (pIRIR ₂₂₅ ), quartz blue light optically stimulated luminescence (BLOSL), quartz blue thermoluminescence (BTL), and quartz red thermoluminescence (RTL). The pIRIR ₂₂₅ and RTL protocols produce consistent age estimates that agree with independent radiocarbon ages and with the timing of the Pinedale deglaciation. This study focuses on two craters, Mary Bay, along the northern shore of Yellowstone Lake, and Pocket Basin in Lower Geyser Basin. The mean pIRIR ₂₂₅ ages from Mary Bay deposits (11.99 ± 0.68 ka) agree with previous radiocarbon constraints. The mean pIRIR ₂₂₅ results from Pocket Basin deposits (13.44 ± 1.06 ka) suggest a history of explosion following Pinedale deglaciation, followed by recent hydrothermal alteration. Luminescence dating techniques are a promising tool for reconstructing the timing of hydrothermal explosions in the Late Pleistocene and Holocene, helping to constrain recurrence intervals of the largest hydrothermal systems, informing risk, and improving hazard assessments.

To improve our understanding of hydrothermal activity on the Yellowstone Plateau volcanic field, we collected and analyzed a large data set of δ²H, δ¹⁸O, and the ³H concentrations of circum‐neutral and alkaline waters. We find that (a) hot springs are fed by recharge throughout the volcanic plateau, likely focused through fractured, permeable tuff units. Previous work had stressed the need for light δ²H water recharge restricted to the northern part of the plateau or recharge during past cold periods. However, new data from the Y‐7 drill hole suggests that recharge is not restricted to a certain area or a cold period. (b) δ¹⁸O values of thermal waters in the geyser basins are shifted from the global meteoric water line by temperature‐dependent water‐rock reactions with higher subsurface temperatures resulting in a greater shift. (c) Large temporal variations in the isotopic composition of meteoric water recharge and small temporal variability in the isotopic composition of hot spring discharge implies that the volume of groundwater in, and around the Yellowstone caldera is substantially larger than the volume of annual water recharge. (d) Hot springs discharged through different rhyolitic units correlate with identifiable differences in δ²H and δ¹⁸O compositions, ³H concentrations, and water chemistry that imply equilibration at different temperatures and travel along different flow paths. (e) Based on measured ³H concentrations, we calculate that hot spring waters in the central part of the geyser basins mostly contain <2% post‐1950 meteoric water, whereas waters discharged at the basin margins contain larger fractions of post‐1950s meteoric water.

Natural radiogenic isotopes (primarily ⁸⁷Sr/⁸⁶Sr) from hot springs in the Upper Geyser Basin of the Yellowstone Plateau volcanic field and associated rocks were used to evaluate groundwater flow patterns, water‐rock reactions, and the extent of mixing between various groundwater sources. Thermal waters have very low uranium concentrations and ²³⁴U/²³⁸U activity ratios near 1.0, which limit their utility as tracers in this reducing setting. Thermal waters have higher Sr concentrations (87Sr/⁸⁶Sr values that vary both temporally at individual discharge sites and between adjacent springs, indicating that conduits tap different subsurface reservoirs to varying degrees. Sr from local rhyolites have ⁸⁷Sr/⁸⁶Sr compositions that bound the range of values observed in groundwater throughout the basin. Non‐boiling springs on the west flank of the basin discharge water with low ⁸⁷Sr/⁸⁶Sr consistent with flow through young volcanic rocks exposed at the surface. Boiling springs in the central basin have higher ⁸⁷Sr/⁸⁶Sr values reflecting interactions with older, more radiogenic volcanic rocks. Variability in upwelling thermal waters requires mixing with a low ⁸⁷Sr/⁸⁶Sr component derived from young lava or glacial sediments, or more likely, from deeper sources of hot groundwater circulating through buried Lava Creek Tuff having intermediate ⁸⁷Sr/⁸⁶Sr. Isotope data constrain basin‐wide output of thermal water to 110–140 kg·s⁻¹. Results underscore the utility of radiogenic Sr isotopes as valuable tracers of hydrothermal flow patterns and improve the understanding of temperature‐dependent water‐rock reactions in one of the largest continental hydrothermal systems on Earth.

Less than a year after the 2018 Kīlauea caldera collapse and eruption, water appeared in newly deepened Halemaʻumaʻu crater. The lake—unprecedented in the written record—grew to a depth of ∼50 m before lava from the December 2020 eruption boiled it away. Surface water heightened concerns of potential phreatic or phreatomagmatic explosions but also offered a new means of possibly identifying eruption precursors. The U.S. Geological Survey Hawaiian Volcano Observatory (HVO) monitored the lake via direct visual observation, webcams, thermal imaging, colorimetry, and laser rangefinders. HVO also employed uncrewed aircraft systems to sample the water and measure near‐lake gas composition. The lake's δD and δ¹⁸O indicate a groundwater source with substantial evaporation. The initial sample had a salinity (total dissolved solids concentration) of 71,000 mg/L and was rich in sulfate (∼53,000 mg/L), iron (∼500 mg/L), and magnesium (∼10,000 mg/L). Subsequent samples were slightly more dilute. The water's pH (∼4), δ³⁴S (+4.3‰), and surface temperatures (up to 85°C) suggest, rather than significant scrubbing of magmatic volatiles, leaching of basalt and reactions with sulfate minerals resulted in high concentrations of sulfate and other solutes. Thermodynamic modeling and precipitate mineralogy indicate that water composition was controlled by iron oxidation and sulfate dissolution. Although the lake exhibited no detectable precursors before the next eruption, and phreatic or phreatomagmatic explosions did not materialize, our multi‐parameter approach to monitoring yielded an enhanced understanding of the hydrologic, geologic, and magmatic conditions that led to the formation of the unique and short‐lived lake.

Volcanic unrest can trigger appreciable change to surface waters such as streams, springs, and volcanic lakes. Magma degassing produces gases and soluble salts that are absorbed into groundwater that feeds streams and lakes. As magma ascends, the amount of heat and degassing will increase, and so will any related geochemical and thermal signal. Subsurface magma movement can cause pressurization that alters hydrostatic head and may induce groundwater discharge. Fluid-pressure changes have been linked to distal volcano-tectonic earthquakes (White and McCausland, 2016; Coulon and others, 2017) and phreatic eruptions (for example, Yamaoka and others, 2016). Clearly, changes in groundwater and surface waters are both indicators of unrest and clues to how and where magma is rising toward the surface. Where possible, it is prudent to incorporate real-time hydrologic data into multiparameter monitoring of restless volcanoes. Hydrologic dynamics can also be tracked by changes in groundwater levels that are commonly measured in shallow boreholes (see chapter K, this volume, on boreholes; Hurwitz and Lowenstern, 2024).Although inferred to be common, relatively few volcano-hydrology anomalies are well documented, and many are essentially anecdotal (Newhall and others, 2001), reflecting the fact that high-resolution time series remain rare. Extreme examples include the 2008 eruption of Nevado del Huila, Colombia, where relatively minor phreatomagmatic eruptions were accompanied by expulsion of as much as 300 million cubic meters of groundwater from fissures high on the volcano (Worni and others, 2011), generating large lahars. Substantial decreases in flow rate from springs about 8 kilometers from the summit of Mayon Volcano, Philippines, have been noted before most eruptions in the 20th century (Newhall and others, 2001). Stream monitoring at Redoubt Volcano in 2009 allowed Werner and others (2012) to recognize that groundwater was unable to absorb (or scrub) the high flux of volcanic gas and that a high CO2/SO2 precursor signal had been evident for 5 months prior to the eruption. A key to better interpreting hydrologic anomalies—or even identifying them—is therefore obtaining adequate baseline data.Most hydrologic monitoring at U.S. volcanoes has been accomplished by intermittent sampling surveys with annual or less frequent sampling (for example, https://hotspringchem.wr.usgs.gov/index.php). More frequent sampling, however, generally is needed to establish reliable baselines. A recent hydrologic and hydrothermal monitoring experiment at 25 sites and 10 of the 12 level 4 (very high threat) volcanoes in the U.S. portion of the Cascade Range demonstrated that there is sufficient temporal variability in hydrothermal fluxes, even during quiescent periods, that one-time measurements will commonly have limited interpretive value (Crankshaw and others, 2018). Thus, surveys are best augmented with data from streamgages (for example, Evans and others, 2004; Bergfeld and others, 2008). Streamflow (water discharge) data allow measured temperature and specific conductance to be converted to heat and solute mass fluxes, which could be insightful parameters for detecting anomalous activity (McCleskey and others, 2012). At the Yellowstone Caldera, long-term monitoring of river solutes has allowed calculation of the chloride flux, a proxy for heat discharge (Hurwitz and others, 2007; McCleskey and others, 2016) from the subsurface magma. This is readily accomplished because data from streamgages are continuously recorded and archived by the U.S. Geological Survey (USGS) National Water Information System (NWIS) (USGS, 2024).Similar studies on stratovolcanoes or shield volcanoes would be scientifically useful, and yet are logistically challenging, requiring streamgages on numerous radial drainages complemented by either frequent manual sampling or numerous deployments of equipment to measure water temperature and specific conductance as a proxy for water chemistry. Another challenge is that some volcanic areas, especially shield volcanoes, are characterized by near-surface porous rocks and soils, such that surface streams are rare and replaced by distant, dilute large-volume springs with only a trace of any original volcanically sourced water (Manga, 2001; Hurwitz and others, 2021).Volcanic lakes are worthy of special attention for monitoring efforts, as their temperature and composition can provide evidence of increased flux of volatile-rich fluids from below. Quantifying changes in volatile and heat release from magma can be simpler in lakes than for volcanoes with radial drainages and no major lakes. Moreover, volcanic lakes pose a range of hazards themselves, including phreatomagmatic eruptions, debris flows, flank collapse, tsunamis, and toxic gas release (Mastin and Witter, 2000; Delmelle and others, 2015; Manville, 2015; Rouwet and others, 2015)—hazards that have historically been responsible for substantial loss of life at many volcanoes worldwide (Manville, 2015). Catastrophic CO2 release at Lake Nyos, Cameroon, in 1986 suffocated about 1,750 people and about 3,500 livestock and was probably triggered by a large landslide into the gas-saturated lake (Kling and others, 1987; Evans and others, 1993). Gas-charged springs in Soda Bay within Clear Lake (California) have caused almost a dozen deaths to bathers in the past hundred years (ABC News, 2000). A 2005 example of lake overturn and abundant gas release was documented at Mount Chiginagak in Alaska (Schaefer and others, 2008) but did not result in any human casualties. Although thermally stratified lakes, which promote trapping of exsolved magmatic gas, tend to develop in tropical regions, the phenomenon can also arise where salinity creates meromixis (a condition in which a lake does not mix completely), as occurs in Mono Lake, California (Jellison and Melack, 1993; Jellison and others, 1998).If magma erupts or flows into a lake, the...

Installation of instrument packages in deep (several hundred to several thousand meters) boreholes near volcanoes is relatively expensive (a few million to tens of millions of U.S. dollars), but can provide a low-noise, high-quality source of geophysical (seismic, strain, tilt, and pore pressure), physical (temperature and water level), and geochemical data. Observations from instruments at depth have the potential to provide insights into processes associated with magma intrusion, unrest, and eruption that would not otherwise be possible (Lowenstern and others, 2017; Eichelberger, 2020). Examples of instrumented boreholes in volcanic areas include the 3-kilometer (km)-deep Long Valley Exploratory Well (LVEW) in California (for example, Priest and others, 1998; Prejean and Ellsworth, 2001; Fischer and others, 2003; Roeloffs and others, 2003; Sorey and others, 2003), the 1,262 meter-deep NSF Well (commonly referred to as the “Keller Well”) within the summit caldera of Kīlauea, Hawaiʻi (Keller and others, 1979; Myren and others, 2006), and the Caribbean Andesite Lava Island-volcano Precision Seismo-geodetic Observatory (CALIPSO) project at Soufrière Hills, Montserrat, which includes a series of four 200-meter (m)-deep holes (for example, Mattioli and others, 2004; Voight and others, 2006). The Plate Boundary Observatory (PBO) of the National Science Foundation’s Earthscope project placed seismometers, tiltmeters, strainmeters, and pore-pressure sensors at depths of 100 to 250 m in more than 100 boreholes scattered in western North America, including at Mount St. Helens, Washington, and Yellowstone Caldera, Wyoming. The total cost for an instrumented PBO borehole ranged from $250,000 to $270,000 U.S. dollars (USD) and a few thousand USD are required annually for maintenance (David Mencin, UNAVCO, written commun., October 2020).

In the past century, most eruptions of Steamboat Geyser in Yellowstone National Park's Norris Geyser Basin were mainly clustered in three episodes: 1961–1969, 1982–1984, and ongoing since 2018. These eruptive episodes resulted in extensive disturbance to surrounding trees. To characterize tree response over time as an indicator of geyser activity adjustments to climate variability, aerial and ground images were analyzed to document changes in tree coverage around the geyser since 1954. Radiocarbon dating of silicified tree remnants from within 14 m of the geyser vent was used to examine geyser response to possible variations in decadal to centennial precipitation patterns. We searched for atypical or absent growth rings in cores from live trees in years associated with large geyser eruptions. Photographs indicate that active eruptive phases have adversely affected trees up to 30 m from the vent, primarily in the dominant downwind direction. Radiocarbon dates indicate that the geyser formed before 1878, in contrast to the birthdate reported in historical documents. Further, the ages of the silicified trees cluster within three episodes that are temporally correlated with periods of relative drought in the Yellowstone region during the 15th–17th centuries. The discontinuous growth of trees around the geyser suggests that changes in eruptive patterns occur in response to decadal to multidecadal droughts. This inference is supported by the lack of silicified specimens with more than 20 annual rings and by the existence of atypical or missing rings in live trees during periods of extended geyser activity.

On 20 December 2020, after more than 2 years of quiescence at Kīlauea Volcano, Hawaiʻi, renewed volcanic activity in the summit crater caused boiling of the water lake over a period of ∼90 min. The resulting water‐rich, electrified plume rose to 11–13 km above sea level, which is among the highest plumes on record for Kīlauea. Although conventional models would infer a high mass flux from explosive magma‐water interaction, the plume was not associated with an infrasound signal indicative of “explosive” activity, nor did it produce a measurable ash‐fall deposit. We use multisensor data to characterize lava‐water interaction and plume generation during this opening phase of the 2020–21 eruption. Satellite, weather radar, and eyewitness observations revealed that the plume was rich in water vapor and hydrometeors but transported less ash than expected from its maximum height. Volcanic lightning flashes detected by ground‐based cameras were confined to freezing altitudes of the upper cloud, suggesting that the ice formation drove the electrification of this plume. The low acoustic energy from lava‐water interaction points to a weakly explosive style of hydrovolcanism. Heat transfer calculations show that the lava to water heat flux was sufficient to boil the lake within 90 min. Limited mixing of lava and water inhibited major steam explosions and fine fragmentation. Results from one‐dimensional plume modeling suggest that the models may underpredict plume height due to overestimation of crosswind air‐entrainment. Our findings shed light on an unusual style of volcanism in which weakly explosive lava‐water interaction generated an outsized plume.

The Upper Geyser Basin at Yellowstone National Park (Wyoming, USA) harbors the greatest concentration of geysers worldwide. Research suggests that individual geysers are not isolated but rather are hydraulically connected in the subsurface with other geysers and thermal springs. To quantify such connections, we combined techniques from machine learning, causal inference, and dynamical systems to characterize the collective eruptive behavior of a set of 10 geysers over 18 months (April 2007 – September 2008) focusing on geyser‐geyser interactions. Model predictions were up to 15 times more accurate when we sought to predict a geyser's eruption time series based on outflow channel temperatures from the network than based on its own time series alone, suggesting the existence of a complex interconnected subsurface groundwater system. On average, cone‐type geysers had larger impacts on other geysers than did fountain‐type geysers. Similarly, cone‐type geysers were on average more insulated from other geysers. However, substantial unexplained variation remained after considering the cone versus fountain dichotomy. Distance between geysers also affected interactions: nearby geysers had stronger effects on focal geysers than did geysers located farther away. Collectively, results support the hypothesis of geyser interdependence at timescales of 5 min–10 days. Our analyses highlight the existence of quantifiable geyser‐to‐geyser interactions that can be resolved through pairwise and system‐level analyses. These findings emphasize the subsurface interconnectedness of thermal features, provide information relevant to visitor experiences in Yellowstone National Park, and suggest strategies for exploring patterns of interdependence that may exist among other episodic geological phenomena.