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   
                    
                            
                            Origin and Properties of Hydrothermal Tremor at Lone Star Geyser, Yellowstone National Park, USA
                        
                    
    
            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   
        
    
                            - Award ID(s):
- 1724986
- PAR ID:
- 10447741
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 125
- Issue:
- 12
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            Abstract With substantial postglacial rhyolite eruptions and ongoing rapid uplift, the Laguna del Maule volcanic field in the southern Andes provides an exceptional opportunity to study the dynamics of an active silicic magmatic system. Using 4,093Parrivals from 137 distant earthquakes recorded by 44 local stations over∼2.25 years, we conduct teleseismic tomography to image the crustal structure down to 40 km below the volcanic field. A prominent low‐velocity body at depths between∼0 and 12 km below sea level (b.s.l.), characterized by a volume of∼500 km3and a peak anomaly of−400 m/s (∼9%), overlaps the location of the upper‐crustal magma reservoir detected in recent gravity and surface wave tomography studies. Its estimated averagePwave velocity of∼4.6 km/s corresponds to an average melt fraction of about 14% and a melt volume of∼70 km3. Petrologic observations are also consistent with generation and storage of rhyolitic melts at depths corresponding to the anomalous zone. Moreover, the tomographic results support a lower crust zone of MASH (melting, assimilation, storage, and homogenization) from a depth of∼25 km to the base of the model, which likely reflects a deep crustal source of magma that contributes to and incubates the shallow silicic reservoir.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 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
- 
            Abstract A new episode of unrest and phreatic/phreatomagmatic/magmatic eruptions occurred at Ambae volcano, Vanuatu, in 2017–2018. We installed a multi-station seismo-acoustic network consisting of seven 3-component broadband seismic stations and four 3-element (26–62 m maximum inter-element separation) infrasound arrays during the last phase of the 2018 eruption episode, capturing at least six reported major explosions towards the end of the eruption episode. The observed volcanic seismic signals are generally in the passband 0.5–10 Hz during the eruptive activity, but the corresponding acoustic signals have relatively low frequencies (< 1 Hz). Apparent very-long-period (< 0.2 Hz) seismic signals are also observed during the eruptive episode, but we show that they are generated as ground-coupled airwaves and propagate with atmospheric acoustic velocity. We observe strongly coherent infrasound waves at all acoustic arrays during the eruptions. Using waveform similarity of the acoustic signals, we detect previously unreported volcanic explosions at the summit vent region based on constant-celerity reverse-time-migration (RTM) analysis. The detected acoustic bursts are temporally related to shallow seismic volcanic tremor (frequency content of 5–10 Hz), which we characterise using a simplified amplitude ratio method at a seismic station pair with different distances from the vent. The amplitude ratio increased at the onset of large explosions and then decreased, which is interpreted as the seismic source ascent and descent. The ratio change is potentially useful to recognise volcanic unrest using only two seismic stations quickly. This study reiterates the value of joint seismo-acoustic data for improving interpretation of volcanic activity and reducing ambiguity in geophysical monitoring.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
