Search for: All records

Volcanic eruption source parameters may be estimated from acoustic pressure recordings dominant at infrasonic frequencies (< 20 Hz), yet uncertainties may be high due in part to poorly understood propagation dynamics. Linear acoustic propagation of volcano infrasound is commonly assumed, but nonlinear processes such as wave steepening may distort waveforms and obscure the sourcing process in recorded waveforms. Here we use a previously developed frequency-domain nonlinearity indicator to quantify spectral changes due to nonlinear propagation primarily in 80 signals from explosions at Yasur Volcano, Vanuatu. We find evidence for$$\le$$$\le$10⁻³ dB/m spectral energy transfer in the band 3–9 Hz for signals with amplitude on the order of several hundred Pa at 200–400 m range. The clarity of the nonlinear spectral signature increases with waveform amplitude, suggesting stronger nonlinear changes for greater source pressures. We observe similar results in application to synthetics generated through finite-difference wavefield simulations of nonlinear propagation, although limitations of the model complicate direct comparison to the observations. Our results provide quantitative evidence for nonlinear propagation that confirms previous interpretations made on the basis of qualitative observations of asymmetric waveforms.

 

Sound waves generated by erupting volcanoes can be used to infer important source dynamics, yet acoustic source‐time functions may be distorted during propagation, even at local recording distances (15 km). The resulting uncertainty in source estimates can be reduced by improving constraints on propagation effects. We aim to quantify potential distortions caused by wave steepening during nonlinear propagation, with the aim of improving the accuracy of volcano‐acoustic source predictions. We hypothesize that wave steepening causes spectral energy transfer away from the dominant source frequency. To test this, we apply a previously developed single‐point, frequency domain, quadspectral density‐based nonlinearity indicator to 30 acoustic signals from Vulcanian explosion events at Sakurajima Volcano, Japan, in an 8‐day data set collected by five infrasound stations in 2013 with 2.3‐ to 6.2‐km range. We model these results with a 2‐D axisymmetric finite‐difference method that includes rigid topography, wind, and nonlinear propagation. Simulation results with flat ground indicate that wave steepening causes up to2 dB (1% of source level) of cumulative upward spectral energy transfer for Sakurajima amplitudes. Correction for nonlinear propagation may therefore provide a valuable second‐order improvement in accuracy for source parameter estimates. However, simulations with wind and topography introduce variations in the indicator spectra on order of a few decibels. Nonrandom phase relationships generated during propagation or at the source may be misinterpreted as nonlinear spectral energy transfer. The nonlinearity indicator is therefore best suited to small source‐receiver distances (e.g.,2 km) and volcanoes with simple sources (e.g., gas‐rich strombolian explosions) and topography.

 

Yasur volcano, Vanuatu is a continuously active open-vent basaltic-andesite stratocone with persistent and long-lived eruptive activity. We present results from a seismo-acoustic field experiment at Yasur, providing locally dense broad-band seismic and infrasonic network coverage from 2016 July 27 to August 3. We corroborate our seismo-acoustic observations with coincident video data from cameras deployed at the crater and on an unoccupied aircraft system (UAS). The waveforms contain a profusion of signals reflecting Yasur’s rapidly occurring and persistent explosive activity. The typical infrasonic signature of Yasur explosions is a classic short-duration and often asymmetric explosion waveform characterized by a sharp compressive onset and wideband frequency content. The dominant seismic signals are numerous repetitive very-long-period (VLP) signals with periods of ∼2–10 s. The VLP seismic events are ‘high-rate’, reoccurring near-continuously throughout the data set with short interevent times (∼20–60 s). We observe variability in the synchronization of seismic VLP and acoustic sources. Explosion events clearly delineated by infrasonic waveforms are underlain by seismic VLPs. However, strong seismic VLPs also occur with only a weak infrasonic expression. Multiplet analysis of the seismic VLPs reveals a systematic progression in the seismo-acoustic source decoupling. The same dominant seismic VLP multiplet occurs with and without surficial explosions and infrasound, and these transitions occur over a timescale of a few days during our field campaign. We subsequently employ template matching, stacking, and full-waveform inversion to image the source mechanism of the dominant VLP multiplet. Inversion of the dominant VLP multiplet stack points to a composite source consisting of either a dual-crack (plus forces) or pipe-crack (plus forces) mechanism. The derived mechanisms correspond to a point-source directly beneath the summit vents with centroid depths in the range ∼900–1000 m below topography. All mechanisms suggest a northeast trending crack dipping relatively shallowly to the northwest and indicate a VLP source centroid and mechanism controlled by a stable structural geologic feature beneath Yasur. We interpret the results in the framework of gas slug ascent through the conduit responsible for Yasur explosions. The VLP mechanism and timing with infrasound (when present) are explained by a shallow-buffered top-down model in which slug ascent is relatively aseismic until reaching the base of a shallow section. Slug disruption in this shallow zone triggers a pressure disturbance that propagates downward and couples at the conduit base (VLP centroid). If the shallow section is open, an explosion propagates to the surface, producing infrasound. In the case of (the same multiplet) VLPs occurring without surficial explosions and weak or no infrasound, the decoupling of the dominant VLPs at ∼900–1000 m depth from surficial explosions and infrasound strongly indicates buffering of the terminal slug ascent. This buffering could be achieved by a variety of conditions at or directly beneath the vents, such as a high-viscosity layer of crystal-rich magma, a debris cap from backfill, a foam layer, or a combination of these. The dominant VLP at Yasur captured by our experiment has a source depth and mechanism separated from surface processes and is stable over time.

 

Acoustic waveform inversions can provide estimates of volume flow rate and erupted mass, enhancing the ability to estimate volcanic emissions. Previous studies have generally assumed a simple acoustic source (monopole); however, more complex and accurate source reconstructions are possible with a combination of equivalent sources (multipole). We deployed a high‐density acoustic network around Yasur volcano, Vanuatu, including acoustic sensors on a tethered aerostat that was moved every ∼15–60 min. Using this unique data set we invert for the acoustic multipole source mechanism using a grid search approach for 80 events to examine volume flow rates and dipole strengths. Our method utilizes finite‐difference time‐domain modeling to obtain the full 3‐D Green's functions that account for topography. Inversion results are compared using a monopole‐only, multipole (monopole and dipole), simulations that do not include topography, and those that use a subset of sensors. We find that the monopole source is a good approximation when topography is considered. However, initial compression amplitude is not fully captured by a monopole source so source directionality cannot be ruled out. The monopole solution is stable regardless of whether a monopole‐only or multipole inversion is performed. Inversions for the dipole components produce estimates consistent with observed source directionality, though these inversions are somewhat unstable given station configurations of typical deployments. Our results suggest that infrasound waveform inversion shows promise for realistic quantitative source estimates, but additional work is necessary to fully explore inversion stability, uncertainty, and robustness.

 

Abstract Over the past two decades (2000–2020), volcano infrasound (acoustic waves with frequencies less than 20 Hz propagating in the atmosphere) has evolved from an area of academic research to a useful monitoring tool. As a result, infrasound is routinely used by volcano observatories around the world to detect, locate, and characterize volcanic activity. It is particularly useful in confirming subaerial activity and monitoring remote eruptions, and it has shown promise in forecasting paroxysmal activity at open-vent systems. Fundamental research on volcano infrasound is providing substantial new insights on eruption dynamics and volcanic processes and will continue to do so over the next decade. The increased availability of infrasound sensors will expand observations of varied eruption styles, and the associated increase in data volume will make machine learning workflows more feasible. More sophisticated modeling will be applied to examine infrasound source and propagation effects from local to global distances, leading to improved infrasound-derived estimates of eruption properties. Future work will use infrasound to detect, locate, and characterize moving flows, such as pyroclastic density currents, lahars, rockfalls, lava flows, and avalanches. Infrasound observations will be further integrated with other data streams, such as seismic, ground- and satellite-based thermal and visual imagery, geodetic, lightning, and gas data. The volcano infrasound community should continue efforts to make data and codes accessible and to improve diversity, equity, and inclusion in the field. In summary, the next decade of volcano infrasound research will continue to advance our understanding of complex volcano processes through increased data availability, sensor technologies, enhanced modeling capabilities, and novel data analysis methods that will improve hazard detection and mitigation. 

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. 

Atmospheric acoustic waves from volcanoes at infrasonic frequencies (0.01–20 Hz) can be used to estimate source parameters for hazard modeling, but signals are often distorted by wavefield interactions with topography, even at local recording distances (<15 km). We present new developments toward a simple empirical approach to estimate attenuation by topographic diffraction at reduced computational cost. We investigate the applicability of a thin screen diffraction relationship developed by Maekawa [1968, doi: https://doi.org/10.1016/0003-682X(68)90020- 0]. We use a 2D axisymmetric finite-difference method to show that this relationship accurately predicts power losses for infrasound diffraction over an idealized kilometer-scale screen; thus validating the scaling to infrasonic wavelengths. However, the Maekawa relationship overestimates attenuation for realistic volcano topography (using Sakurajima Volcano as an example). The attenuating effect of diffraction may be counteracted by constructive interference of multiple reflections along concave volcano slopes. We conclude that the Maekawa relationship is insufficient as formulated for volcano infrasound, and suggest modifications that may improve the prediction capability.