An official website of the United States government
Abstract Steamboat Geyser in Yellowstone National Park is the tallest active geyser on Earth and is believed to have hydrologic connection to Cistern Spring, a hydrothermal pool ∼100 m southwest from the geyser vent. Despite broad scientific interest, rare episodic Steamboat eruptions have made it difficult to study its eruption dynamics and underground plumbing architecture. In response to the recent reactivation of Steamboat, which has produced more than 130 eruptions since March 2018, we deployed a dense seismic nodal array surrounding the enigmatic geyser in the summer of 2019. The array recorded abundant 1–5 Hz hydrothermal tremor originating from phase‐transition events within both Steamboat Geyser and Cistern Spring. To constrain the spatiotemporal distribution of the tremor sources, an interferometric‐based polarization analysis was developed. The observed tremor locations indicate that the conduit beneath Steamboat is vertical and extends down to ∼120 m depth and the plumbing of Cistern includes a shallow vertical conduit connecting with a deep, large, and laterally offset reservoir ∼60 m southeast of the surface pool. No direct connection between Steamboat and Cistern plumbing structures is found. The temporal variation of tremor combined within situtemperature and water depth measurements of Cistern reveals interaction between Steamboat and Cistern throughout the eruption/recharge cycles. The observed delayed responses of Cistern Spring in reaction to Steamboat eruptions and recharge suggest that the two plumbing structures may be connected through a fractured/porous medium instead of a direct open channel, consistent with our inferred plumbing structure.
more »
« less
Imaging the Deep Subsurface Plumbing of Old Faithful Geyser From Low‐Frequency Hydrothermal Tremor Migration
Abstract Old Faithful Geyser in Yellowstone is one of the most well‐known hydrothermal features in the world. Despite abundant geophysical studies, the structure of Old Faithful's plumbing system beneath ~20‐m depth remained largely elusive. By deploying a temporary dense three‐component geophone array, we observe 1–5 Hz low‐frequency hydrothermal tremor originating from Old Faithful's deeper conduit. By applying seismic interferometry and polarization analyses, we track seismic tremor source migration throughout the eruption/recharge cycle. The tremor source drops rapidly to ~80‐m depth right after the eruption and gradually ascends vertically back to ~20‐m depth, coinciding with the previously inferred bubble trap location. Likely excited by the liquid/steam phase transition, the observed tremor source migration can provide new constraints on the recharge process and deeper conduit geometry. Combined with the shallow conduit structure from previous studies, these results provide constraints on the major fluid pathway down to 80‐m depth.
more »
« less
- Award ID(s):
- 1753362
- PAR ID:
- 10360179
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 46
- Issue:
- 13
- ISSN:
- 0094-8276
- Page Range / eLocation ID:
- p. 7315-7322
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Geysers are rare geologic features that intermittently discharge liquid water and steam driven by heating and decompression boiling. The cause of variability in eruptive styles and the associated seismic signals are not well understood. Data collected from five broadband seismometers at Lone Star Geyser, Yellowstone National Park are used to determine the properties, location, and temporal patterns of hydrothermal tremor. The tremor is harmonic at some stages of the eruption cycle and is caused by near‐periodic repetition of discrete seismic events. Using the polarization of ground motion, we identify the location of tremor sources throughout several eruption cycles. During preplay episodes (smaller eruptions preceding the more vigorous major eruption), tremor occurs at depths of 7–10 m and is laterally offset from the geyser's cone by ~5 m. At the onset of the main eruption, tremor sources migrate laterally and become shallower. As the eruption progresses, tremor sources migrate along the same path but in the opposite direction, ending where preplay tremor originates. The upward and then downward migration of tremor sources during eruptions are consistent with warming of the conduit followed by evacuation of water during the main eruption. We identify systematic relations among the two types of preplays, discharge, and the main eruption. A point‐source moment tensor fit to low‐frequency waveforms of an individual tremor event using half‐space velocity models indicates averageVS ≳ 0.8 km/s, source depths ~4–20 m, and moment tensors with primarily positive isotropic and compensated linear vector dipole moments.more » « less
-
The intensity of explosive volcanic eruptions is correlated with the amplitude of eruption tremor, a ubiquitously observed seismic signal during eruptions. Here we expand upon a recently introduced theoretical model that attributes eruption tremor to particle impacts and dynamic pressure changes in the turbulent flow above fragmentation (Gestrich et al., 2020). We replace their point source model with Rayleigh wave Green's functions with full Green's functions and account for depth variation of input fields using conduit flow models. The latter self-consistently capture covariation of input fields like particle velocity, particle volume fraction, and density. Body wave contributions become significant above 2-3 Hz, bringing the power spectral density (PSD) closer to observations. Conditions at the vent are not representative of flow throughout the tremor source region and using these values overestimates tremor amplitude. Particle size and its depth distribution alter the PSD and where dominant source contributions arise within the conduit. Solutions with decreasing mass eruption rate, representing a waning eruption, reveal a shift in the dominant tremor contribution from turbulence to particle impacts. Our work demonstrates the ability to integrate conduit flow modeling with volcano seismology studies of eruption tremor, providing an opportunity to link observations to eruptive processes.more » « less
-
Abstract Fragmentation plays a critical role in eruption explosivity by influencing the eruptive jet and plume dynamics that may initiate hazards such as pyroclastic flows. The mechanics and progression of fragmentation during an eruption are challenging to constrain observationally, limiting our understanding of this important process. In this work, we explore seismic radiation associated with unsteady fragmentation. Seismic force and moment tensor fluctuations from unsteady fragmentation arise from fluctuations in fragmentation depth and wall shear stress (e.g., from viscosity variations). We use unsteady conduit flow models to simulate perturbations to a steady‐state eruption from injections of heterogeneous magma (specifically, variable magma viscosity due to crystal volume fraction variations). Changes in wall shear stress and pressure determine the seismic force and moment histories, which are used to calculate synthetic seismograms. We consider three heterogeneity profiles: Gaussian pulse, sinusoidal, and stochastic. Fragmentation of a high‐crystallinity Gaussian pulse produces a distinct very‐long‐period seismic signature and associated reduction in mass eruption rate, suggesting joint use of seismic, infrasound, and plume monitoring data to identify this process. Simulations of sinusoidal injections quantify the relation between the frequency or length scale of heterogeneities passing through fragmentation and spectral peaks in seismograms, with velocity seismogram amplitudes increasing with frequency. Stochastic composition variations produce stochastic seismic signals similar to observed eruption tremor, though computational limitations restrict our study to frequencies less than 0.25 Hz. We suggest that stochastic fragmentation fluctuations could be a plausible eruption tremor source.more » « less
-
Abstract The 2018 eruption of Sierra Negra volcano, Galápagos, Ecuador has provided new insights into the mechanisms of caldera resurgence, subsidence, and fissuring at basaltic shield volcanoes. Here, we integrate local (∼0.4 km) seismo‐acoustic records and regional (∼85 km) infrasound array data to present new observations of the 2018 Sierra Negra eruption with improved time and spatial resolutions. These observations include: air‐to‐ground coupling ∼2 hr before the time of the eruption onset, migration of the infrasound tremor from 22:54 June 26 to 12:31 June 27 UT (all times in UT), and persistent infrasound detections during the weeks between 5 July and 18 August from an area that does not coincide with the previously documented eruptive fissures. We interpret air‐to‐ground coupling as infrasound tremor generated in the nearby fissures before the main eruptive phase started, although ambiguity remains in interpreting a single seismic‐infrasonic sensor pair. The progressive location change of the infrasound tremor agrees with the migration of the eruption down the north flank of Sierra Negra at a rate of ∼0.15 ± 0.04 m/s. The weeks‐long persistent detections coincide with a region that has thermal anomalies, co‐eruptive deformation, lava fields, and geological features that could be interpreted as multiple lava tube skylights. Our observations and interpretations provide constraints on the mechanisms underlying fissure formation and magma emplacement at Sierra Negra.more » « less
