An official website of the United States government
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
Ultra and Very Long Period Seismic Signatures of Unsteady Eruptions Predicted From Conduit Flow Models
Abstract Explosive volcanic eruptions radiate seismic waves as a consequence of pressure and shear traction changes within the conduit/chamber system. Kinematic source inversions utilize these waves to determine equivalent seismic force and moment tensor sources, but relation to eruptive processes is often ambiguous and nonunique. In this work, we provide an alternative, forward modeling approach to calculate moment tensor and force equivalents of a model of eruptive conduit flow and chamber depressurization. We explain the equivalence of two seismic force descriptions, the first in terms of traction changes on conduit/chamber walls, and the second in terms of changes in magma momentum, weight, and momentum transfer to the atmosphere. Eruption onset is marked by a downward seismic force, associated with loss of restraining shear tractions from fragmentation. This is followed by a much larger upward seismic force from upward drag of ascending magma and reduction of magma weight remaining in the conduit/chamber system. The static force is upward, arising from weight reduction. We calculate synthetic seismograms to examine the expression of eruptive processes at different receiver distances. Filtering these synthetics to the frequency band typically resolved by broadband seismometers produces waveforms similar to very long period seismic events observed in strombolian and vulcanian eruptions. However, filtering heavily distorts waveforms, accentuating processes in early, unsteady parts of eruptions and eliminating information about longer (ultra long period time scale depressurization and weight changes that dominate unfiltered seismograms. Our workflow can be utilized to directly and quantitatively connect eruption models with seismic observations.
more »
« less
- Award ID(s):
- 1930979
- PAR ID:
- 10445076
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 127
- Issue:
- 6
- ISSN:
- 2169-9313
- 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
-
Abstract Water and carbon dioxide are the most abundant volatile components in terrestrial magmas. As they exsolve into magmatic vapour, they promote magma buoyancy, accelerating ascent and modulating eruptive dynamics. It is commonly thought that an increase in pre-eruptive volatile content produces an increase in eruption intensity. Using a conduit model for basaltic eruptions, covering the upper 6 km of conduit, we show that for the same chamber conditions mass eruption rate is not affected by CO2content, whereas an increase in H2O up to 10 wt.% produces an increase in eruption rate of an order of magnitude. It is only when CO2is injected in the magma reservoir from an external source that the resulting pressurisation will generate a strong increase in eruption rate. Results also show that ascent velocity and fragmentation depth are strongly affected by pre-eruptive volatile contents demonstrating a link between volatile content and eruptive style.more » « less
-
Abstract Volcanic summit craters are typically noted to form by roof collapse into a depressurized magma chamber or by explosive excavation. Recent examples of effusive activity (e.g., Kilauea Volcano, Hawai'i) allowed specifically for quantification of the collapse process. However, small spatiotemporal morphologic change related to background mass wasting and low‐level explosive activity has not been well quantified in volcanic craters. Telica volcano, Nicaragua, is a persistently restless basaltic‐andesite stratovolcano. Telica's persistent restlessness is caused by a long‐lived magmatic‐hydrothermal system with high‐temperature crater fumaroles and low‐frequency seismicity, punctuated by subdecadal, low‐explosivity (VEI 1–2) phreatic eruptions. We use photographic observations (1994 to 2017) and structure‐from‐motion point cloud construction and differencing (2011 to 2017) to analyze changes at Telica in the context of summit crater formation and eruptive precursors. Crater wall retreat (up to 40 m) spatially correlates with long‐lived high‐temperature fumaroles in the crater walls, whereas eruptions eject material (>5 m) from the crater floor through vent formation and/or clearing. These processes sustain a morphology similar to that of pit craters but without a shallow depressurized magma chamber. Our observations indicate system‐wide sealing prior to eruption by viscous magma in the conduit and eruption of a dome in 2017 and hydrothermal mineralization, not from vent covering talus; though, vent covering talus can redirect the shallow conduit. This study shows promise for photogrammetric techniques in correlating morphologic change with summit crater formation and volcanic activity and the power of long‐term visual observations in understanding active volcanic processes.more » « less
-
Key Points The 15 January 2022 Hunga Tonga‐Hunga Ha'apai eruption had four episodic seismic subevents with similar waveforms within ∼300 s An unusual upward force jump‐started each subevent A magma hammer explains the force and estimates the subsurface magma mass flux which fits the vent discharge rate based on satellite datamore » « less
