@article{EiblHainzlVeselyetal.2019, author = {Eibl, Eva P. S. and Hainzl, Sebastian and Vesely, Nele I. K. and Walter, Thomas R. and Jousset, Philippe and Hersir, Gylfi Pall and Dahm, Torsten}, title = {Eruption interval monitoring at strokkur Geyser, Iceland}, series = {Geophysical research letters}, volume = {47}, journal = {Geophysical research letters}, number = {1}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1029/2019GL085266}, pages = {10}, year = {2019}, abstract = {Geysers are hot springs whose frequency of water eruptions remain poorly understood. We set up a local broadband seismic network for 1 year at Strokkur geyser, Iceland, and developed an unprecedented catalog of 73,466 eruptions. We detected 50,135 single eruptions but find that the geyser is also characterized by sets of up to six eruptions in quick succession. The number of single to sextuple eruptions exponentially decreased, while the mean waiting time after an eruption linearly increased (3.7 to 16.4 min). While secondary eruptions within double to sextuple eruptions have a smaller mean seismic amplitude, the amplitude of the first eruption is comparable for all eruption types. We statistically model the eruption frequency assuming discharges proportional to the eruption multiplicity and a constant probability for subsequent events within a multituple eruption. The waiting time after an eruption is predictable but not the type or amplitude of the next one.
Plain Language Summary Geysers are springs that often erupt in hot water fountains. They erupt more often than volcanoes but are quite similar. Nevertheless, it is poorly understood how often volcanoes and also geysers erupt. We created a list of 73,466 eruption times of Strokkur geyser, Iceland, from 1 year of seismic data. The geyser erupted one to six times in quick succession. We found 50,135 single eruptions but only 1 sextuple eruption, while the mean waiting time increased from 3.7 min after single eruptions to 16.4 min after sextuple eruptions. Mean amplitudes of each eruption type were higher for single eruptions, but all first eruptions in a succession were similar in height. Assuming a constant heat inflow at depth, we can predict the waiting time after an eruption but not the type or amplitude of the next one.}, language = {en} } @article{EiblMuellerWalteretal.2021, author = {Eibl, Eva P. S. and M{\"u}ller, Daniel and Walter, Thomas R. and Allahbakhshi, Masoud and Jousset, Philippe and Hersir, Gylfi P{\´a}ll and Dahm, Torsten}, title = {Eruptive cycle and bubble trap of Strokkur Geyser, Iceland}, series = {Journal of geophysical research : JGR. B: Solid earth}, volume = {126}, journal = {Journal of geophysical research : JGR. B: Solid earth}, number = {4}, publisher = {Wiley}, address = {Hoboken, NJ}, issn = {2169-9313}, doi = {10.1029/2020JB020769}, pages = {20}, year = {2021}, abstract = {The eruption frequency of geysers can be studied easily on the surface. However, details of the internal structure including possible water and gas filled chambers feeding eruptions and the driving mechanisms often remain elusive. We used a multidisciplinary network of seismometers, video cameras, water pressure sensors and one tiltmeter to study the eruptive cycle, internal structure, and mechanisms driving the eruptive cycle of Strokkur geyser in June 2018. An eruptive cycle at Strokkur always consists of four phases: (1) Eruption, (2) post-eruptive conduit refilling, (3) gas filling of the bubble trap, and (4) regular bubble collapse at shallow depth in the conduit. For a typical single eruption 19 +/- 4 bubble collapses occur in Phase 3 and 8 +/- 2 collapses in Phase 4 at a mean spacing of 1.52 +/- 0.29 and 24.5 +/- 5.9 s, respectively. These collapses release latent heat to the fluid in the bubble trap (Phase 3) and later to the fluid in the conduit (Phase 4). The latter eventually reaches thermodynamic conditions for an eruption. Single to sextuple eruptions have similar spacings between bubble collapses and are likely fed from the same bubble trap at 23.7 +/- 4.4 m depth, 13-23 m west of the conduit. However, the duration of the eruption and recharging phase linearly increases likely due to a larger water, gas and heat loss from the system. Our tremor data provides documented evidence for a bubble trap beneath a pool geyser.}, language = {en} } @article{NamikiRivaltaWoithetal.2018, author = {Namiki, Atsuko and Rivalta, Eleonora and Woith, Heiko and Willey, Timothy and Parolai, Stefano and Walter, Thomas R.}, title = {Volcanic activities triggered or inhibited by resonance of volcanic edifices to large earthquakes}, series = {Geology}, volume = {47}, journal = {Geology}, number = {1}, publisher = {American Institute of Physics}, address = {Boulder}, issn = {0091-7613}, doi = {10.1130/G45323.1}, pages = {67 -- 70}, year = {2018}, abstract = {The existence of a causal link between large earthquakes and volcanic unrest is widely accepted. Recent observations have also revealed counterintuitive negative responses of volcanoes to large earthquakes, including decreased gas emissions and subsidence in volcanic areas. In order to explore the mechanisms that could simultaneously explain both the positive and negative responses of volcanic activity to earthquakes, we here focus on the role played by topography. In the laboratory, we shook a volcanic edifice analogue, made of gel, previously injected with a buoyant fluid. We find that shaking triggers rapid migration of the buoyant fluid upward, downward, or laterally, depending on the fluid's buoyancy and storage depth; bubbly fluids stored at shallow depth ascend, while low-buoyancy fluids descend or migrate laterally. The migration of fluids induced by shaking is two orders of magnitude faster than without shaking. Downward or lateral fluid migration may decrease volcanic gas emissions and cause subsidence as a negative response, while upward migration is consistent both with an increase in volcanic activity and immediate unrest (deformation and seismicity) after large earthquakes. The fluid migration is more efficient when the oscillation frequency is close to the resonance frequency of the edifice. The resonance frequency for a 30-km-wide volcanic mountain range, such as those where subsidence was observed, is ∼0.07 Hz. Only large earthquakes are able to cause oscillation at such low frequencies.}, language = {en} } @article{TuGeWangetal.2014, author = {Tu, Rui and Ge, Maorong and Wang, Rongjiang and Walter, Thomas R.}, title = {A new algorithm for tight integration of real-time GPS and strong-motion records, demonstrated on simulated, experimental, and real seismic data}, series = {Journal of seismology}, volume = {18}, journal = {Journal of seismology}, number = {1}, publisher = {Springer}, address = {Dordrecht}, issn = {1383-4649}, doi = {10.1007/s10950-013-9408-x}, pages = {151 -- 161}, year = {2014}, abstract = {The complementary advantages of GPS and seismic measurements are well recognized in seismotectonic monitoring studies. Therefore, integrated processing of the two data streams has been proposed recently in an attempt to obtain accurate and reliable information of surface displacements associated with earthquakes. A hitherto still critical issue in the integrated processing is real-time detection and precise estimation of the transient baseline error in the seismic records. Here, we report on a new approach by introducing the seismic acceleration corrected by baseline errors into the state equation system. The correction is performed and regularly updated in short epochs (with increments which may be as short as seconds), so that station position, velocity, and acceleration can be constrained very tightly and baseline error can be estimated as a random-walk process. With the adapted state equation system, our study highlights the use of a new approach developed for integrated processing of GPS and seismic data by means of sequential least-squares adjustment. The efficiency of our approach is demonstrated and validated using simulated, experimental, and real datasets. The latter were collected at collocated GPS and seismic stations around the 4 April 2010, E1 Mayor-Cucapah earthquake (Mw, 7.2). The results have shown that baseline errors of the strong-motion sensors are corrected precisely and high-precision seismic displacements are real-timely obtained by the new approach.}, language = {en} } @article{TuWangWalteretal.2014, author = {Tu, Rui and Wang, Rongjiang and Walter, Thomas R. and Diao, FaQi}, title = {Adaptive recognition and correction of baseline shifts from collocated GPS and accelerometer using two phases Kalman filter}, series = {Advances in space research}, volume = {54}, journal = {Advances in space research}, number = {9}, publisher = {Elsevier}, address = {Oxford}, issn = {0273-1177}, doi = {10.1016/j.asr.2014.07.008}, pages = {1924 -- 1932}, year = {2014}, abstract = {The real-time recognition and precise correction of baseline shifts in strong-motion records is a critical issue for GPS and accelerometer combined processing. This paper proposes a method to adaptively recognize and correct baseline shifts in strong-motion records by utilizing GPS measurements using two phases Kalman filter. By defining four kinds of learning statistics and criteria, the time series of estimated baseline shifts can be divided into four time intervals: initialization, static, transient and permanent. During the time interval in which the transient baseline shift is recognized, the dynamic noise of the Kalman filter system and the length of the baseline shifts estimation window are adaptively adjusted to yield a robust integration solution. The validations from an experimental and real datasets show that acceleration baseline shifts can be precisely recognized and corrected, thus, the combined system adaptively adjusted the estimation strategy to get a more robust solution. (C) 2014 COSPAR. Published by Elsevier Ltd. All rights reserved.}, language = {en} } @article{ZornLeCorvecVarleyetal.2019, author = {Zorn, Edgar Ulrich and Le Corvec, Nicolas and Varley, Nick R. and Salzer, Jacqueline T. and Walter, Thomas R. and Navarro-Ochoa, Carlos and Vargas-Bracamontes, Dulce M. and Thiele, Samuel T. and Ar{\´a}mbula Mendoza, Ra{\´u}l}, title = {Load stress controls on directional lava dome growth at Volcan de Colima, Mexico}, series = {Frontiers in Earth Science}, volume = {7}, journal = {Frontiers in Earth Science}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-6463}, doi = {10.3389/feart.2019.00084}, pages = {18}, year = {2019}, abstract = {During eruptive activity of andesitic stratovolcanoes, the extrusion of lava domes, their collapse and intermittent explosions are common volcanic hazards. Many lava domes grow in a preferred direction, in turn affecting the direction of lava flows and pyroclastic density currents. Access to active lava domes is difficult and hazardous, so detailed data characterizing lava dome growth are typically limited, keeping the processes controlling the directionality of extrusions unclear. Here we combine TerraSAR-X satellite radar observations with high-resolution airborne photogrammetry to assess morphological changes, and perform finite element modeling to investigate the impact of loading stress on shallow magma ascent directions associated with lava dome extrusion and crater formation at Volcan de Colima, Mexico. The TerraSAR-X data, acquired in similar to 1-m resolution spotlight mode, enable us to derive a chronology of the eruptive processes from intensity-based time-lapse observations of the general crater and dome evolution. The satellite images are complemented by close-range airborne photos, processed by the Structure-from-Motion workflow. This allows the derivation of high-resolution digital elevation models, providing insight into detailed loading and unloading features. During the observation period from Jan-2013 to Feb-2016, we identify a dominantly W-directed dome growth and lava flow production until Jan-2015. In Feb-2015, following the removal of the active summit dome, the surface crater widened and elongated along a NE-SW axis. Later in May-2015, a new dome grew toward the SW of the crater while a separate vent developed in the NE of the crater, reflecting a change in the direction of magma ascent and possible conduit bifurcation. Finite element models show a significant stress change in agreement with the observed magma ascent direction changes in response to the changing surface loads, both for loading (dome growth) and unloading (crater forming excavation) cases. These results allow insight into shallow dome growth dynamics and the migration of magma ascent in response to changing volcano summit morphology. They further highlight the importance of detailed volcano summit morphology surveillance, as changes in direction or location of dome extrusion may have major implications regarding the directions of potential volcanic hazards, such as pyroclastic density currents generated by dome collapse.}, language = {en} }