Protracted episodes of 0.5–7 Hz pre-eruptive volcanic tremor (PVT) are common at active stratovolcanoes. Reliable links to processes related to magma movement consequently enable a potential to use properties of PVT as diagnostic eruptive precursors. A challenging feature of PVT is that generic spectral and amplitude properties of the signal evolve similarly, independent of widely varying volcano structures and conduit geometries on which most physical models rely. The ‘magma wagging’ model introduced in Jellinek & Bercovici (2011) and extended by Bercovici et al. (2013), Liao et al. and Liao & Bercovici (2018) makes progress because it depends on magma dynamics that are only weakly sensitive to volcano architecture: The flow of gas through a permeable foamy annulus of gas bubbles excites, modulates and maintains a wagging oscillation of a central magma column rising in an erupting conduit. ‘Magma wagging’ and resulting PVT are driven through an energy transfer from a ‘Bernoulli mode’ related to azimuthal variations in annular gas flow speeds. Consistent with observations, spectral and amplitude properties of PVT are predicted to evolve before an eruption as the width of the annulus decreases with increased gas fluxes. To confirm this critical Bernoulli-to-wagging energy transfer we use extensive experiments and restricted numerical simulations on wagging oscillations excited on analogue viscoelastic columns by annular air flows. We also explore sensitivities of the spatial and temporal characters of wagging to asymmetric annular air flows that are intractable in the existing magma wagging model and expected to occur in nature with spatial variations in annulus permeability. From high-resolution time-series of linear and orbital displacements of analogue column tops and time-series of axial deflections and accelerations of the column centre line, we characterize the excitation, evolution, and steady-state oscillations in unprecedented detail over a broad range of conditions. We show that the Bernoulli mode corresponds to the timescale for the buildup of axial elastic bending stresses in response to pressure variations related to air flows over the heights of columns. We identify three distinct wagging modes: (i) rotational (cf. Liao et al. 2018); (ii) mixed-mode and (iii) chaotic. Rotational modes are favoured for symmetric, high intensity forcing and a maximal delivery of mechanical energy to the fundamental magma wagging mode. Mixed-mode oscillations regimes are favoured for a symmetric, intermediate intensity forcing. Chaotic modes, involving the least efficient delivery of energy to the fundamental mode, occur for asymmetric forcing and where the intensity of imposed airflow is low. Numerical simulations also show that where forcing frequencies are comparable to a natural mode of free oscillation, power delivered by peripheral air flows is concentrated at the lowest frequency fundamental mode generally and spread among higher frequency natural modes where air pressure and column elastic forces are comparable. Our combined experimental and numerical results make qualitative predictions for the evolution of the character of volcanic tremor and its expression in seismic or infrasound arrays during natural events that is testable in field-based studies of PVT and syn-eruptive volcanic tremor.
- Award ID(s):
- 2137467
- NSF-PAR ID:
- 10384560
- Date Published:
- Journal Name:
- Scientific Reports
- Volume:
- 12
- Issue:
- 1
- ISSN:
- 2045-2322
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
SUMMARY -
Subduction zones represent the interface between Earth’s interior (crust and mantle) and exterior (atmosphere and oceans), where carbon and other volatile elements are actively cycled between Earth reservoirs by plate tectonics. Helium is a sensitive tracer of volatile sources and can be used to deconvolute mantle and crustal sources in arcs; however it is not thought to be recycled into the mantle by subduction processes. In contrast, carbon is readily recycled, mostly in the form of carbon-rich sediments, and can thus be used to understand volatile delivery via subduction. Further, carbon is chemically-reactive and isotope fractionation can be used to determine the main processes controlling volatile movements within arc systems. Here, we report helium isotope and abundance data for 42 deeply-sourced fluid and gas samples from the Central Volcanic Zone (CVZ) and Southern Volcanic Zone (SVZ) of the Andean Convergent Margin (ACM). Data are used to assess the influence of subduction parameters (e.g., crustal thickness, subduction inputs, and convergence rate) on the composition of volatiles in surface volcanic fluid and gas emissions. He isotopes from the CVZ backarc range from 0.1 to 2.6 R A ( n = 23), with the highest values in the Puna and the lowest in the Sub-Andean foreland fold-and-thrust belt. Atmosphere-corrected He isotopes from the SVZ range from 0.7 to 5.0 R A ( n = 19). Taken together, these data reveal a clear southeastward increase in 3 He/ 4 He, with the highest values (in the SVZ) falling below the nominal range associated with pure upper mantle helium (8 ± 1 R A ), approaching the mean He isotope value for arc gases of (5.4 ± 1.9 R A ). Notably, the lowest values are found in the CVZ, suggesting more significant crustal inputs (i.e., assimilation of 4 He) to the helium budget. The crustal thickness in the CVZ (up to 70 km) is significantly larger than in the SVZ, where it is just ∼40 km. We suggest that crustal thickness exerts a primary control on the extent of fluid-crust interaction, as helium and other volatiles rise through the upper plate in the ACM. We also report carbon isotopes from ( n = 11) sites in the CVZ, where δ 13 C varies between −15.3‰ and −1.2‰ [vs. Vienna Pee Dee Belemnite (VPDB)] and CO 2 / 3 He values that vary by over two orders of magnitude (6.9 × 10 8 –1.7 × 10 11 ). In the SVZ, carbon isotope ratios are also reported from ( n = 13) sites and vary between −17.2‰ and −4.1‰. CO 2 / 3 He values vary by over four orders of magnitude (4.7 × 10 7 –1.7 × 10 12 ). Low δ 13 C and CO 2 / 3 He values are consistent with CO 2 removal (e.g., calcite precipitation and gas dissolution) in shallow hydrothermal systems. Carbon isotope fractionation modeling suggests that calcite precipitation occurs at temperatures coincident with the upper temperature limit for life (122°C), suggesting that biology may play a role in C-He systematics of arc-related volcanic fluid and gas emissions.more » « less
-
null (Ed.)Abstract Volcán Quizapu, Chile, is an under-monitored volcano that was the site of two historical eruptions: an effusive eruption in 1846–1847 and a Plinian eruption in 1932, both of which discharged ∼5 km3 (dense rock equivalent) of lava and/or tephra. The majority of material erupted in both cases is trachydacite, nearly identical for each event. We present H2O-saturated, phase equilibrium experiments on this end-member dacite magma, using a pumice sample from the 1932 eruption as the main starting material. At an oxygen fugacity (fO2) of ∼NNO + 0·2 (where NNO is the nickel–nickel oxide buffer), the phase assemblage of An25–30 plagioclase + amphibole + orthopyroxene, without biotite, is stable at 865 ± 10 °C and 110 ± 20 MPa H2O pressure (PH2O), corresponding to ∼4 km depth. At these conditions, experiments also reproduce the quenched glass composition of the starting pumice. At slightly higher PH2O and below 860 °C, biotite joins the equilibrium assemblage. Because biotite is not part of the observed Quizapu phase assemblage, its presence places an upper limit on PH2O. At the determined storage PH2O of ∼110 MPa, H2O undersaturation of the magma with XH2Ofluid = 0·87 would align Ptotal to mineral-based geobarometry estimates of ∼130 MPa. However, XH2Ofluid < 1 is not required to reproduce the Quizapu dacite phase assemblage and compositions. A second suite of experiments at lower fO2 shows that the stability fields of the hydrous silicates (amphibole and biotite) are significantly restricted at NNO – 2 relative to NNO + 0·2. Additional observations of Quizapu lava and pumice samples support the existing hypothesis that rapid pre-eruptive heating drove the effusive 1846–1847 eruption, with important refinements. We demonstrate that microlites in the end-member dacite lavas are consistent with in situ crystallization (during ascent), rather than transfer from an andesite. In one end-member dacite lava, newly identified reverse zoning in orthopyroxene and incipient destabilization of amphibole are consistent with small degrees of heating. Our work articulates a clear direction for future Quizapu studies, which are warranted given the active nature of the Cerro Azul–Descabezado Grande volcanic axis.more » « less
-
more » « less
Abstract On Dec. 22, 2018, at approximately 20:55–57 local time, Anak Krakatau volcano, located in the Sunda Straits of Indonesia, experienced a major lateral collapse during a period of eruptive activity that began in June. The collapse discharged volcaniclastic material into the 250 m deep caldera southwest of the volcano, which generated a tsunami with runups of up to 13 m on the adjacent coasts of Sumatra and Java. The tsunami caused at least 437 fatalities, the greatest number from a volcanically-induced tsunami since the catastrophic explosive eruption of Krakatau in 1883 and the sector collapse of Ritter Island in 1888. For the first time in over 100 years, the 2018 Anak Krakatau event provides an opportunity to study a major volcanically-generated tsunami that caused widespread loss of life and significant damage. Here, we present numerical simulations of the tsunami, with state-of the-art numerical models, based on a combined landslide-source and bathymetric dataset. We constrain the geometry and magnitude of the landslide source through analyses of pre- and post-event satellite images and aerial photography, which demonstrate that the primary landslide scar bisected the Anak Krakatau volcano, cutting behind the central vent and removing 50% of its subaerial extent. Estimated submarine collapse geometries result in a primary landslide volume range of 0.22–0.30 km3, which is used to initialize a tsunami generation and propagation model with two different landslide rheologies (granular and fluid). Observations of a single tsunami, with no subsequent waves, are consistent with our interpretation of landslide failure in a rapid, single phase of movement rather than a more piecemeal process, generating a tsunami which reached nearby coastlines within ~30 minutes. Both modelled rheologies successfully reproduce observed tsunami characteristics from post-event field survey results, tide gauge records, and eyewitness reports, suggesting our estimated landslide volume range is appropriate. This event highlights the significant hazard posed by relatively small-scale lateral volcanic collapses, which can occur
en-masse , without any precursory signals, and are an efficient and unpredictable tsunami source. Our successful simulations demonstrate that current numerical models can accurately forecast tsunami hazards from these events. In cases such as Anak Krakatau’s, the absence of precursory warning signals together with the short travel time following tsunami initiation present a major challenge for mitigating tsunami coastal impact. -
more » « less
Abstract The noble gas signature of incoming Pacific Bottom Water (PBW), when compared to North Atlantic Deep Water, indicates the addition of 450 ± 70 GT a−1glacial melt water to form AABW and subsequently PBW. The downstream evolution of this signature between the southern (20°S to equator) and northern (25°–45°N) bottom waters indicates a decrease in sea level pressure around Antarctica over the past two millennia. Vertical profiles of noble gases in the deep Pacific show exponential relationships with depth with scale heights identical to temperature and salinity. Unlike the other noble gases, helium isotopes show evidence of mid‐depth injection of non‐atmospheric helium. Using observed deviations from exponential behavior, we quantify its magnitude and isotope ratio. There is a clear latitude trend in the isotope ratio of this added helium that decreases from a high exceeding 9 RA(atmospheric3He/4He ratio) in the south to around 8 RAnear the equator. North of 30–40°N, it systematically decreases northward to a low of ∼2 RAnorth of 50°N. This decline results from a combination of northward decline in seafloor spreading, release of radiogenic helium from increased sediment thickness, and the possible emission of radiogenic helium through cold seeps along the Alaskan and North American margins. Finally, we derive an improved method of computing the excess helium isotope concentrations and that the distributions of bottom water3HeXS/4HeXSare consistent with what is known about bottom water flow patterns and the input of low3He/4He sedimentary helium.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1038/s41598-022-23991-3
