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.

 

Acoustic source inversions estimate the mass flow rate of volcanic explosions or yield of chemical explosions and provide insight into potential source directionality. However, the limitations of applying these methods to complex sources and their ability to resolve a stable solution have not been investigated in detail. We perform synthetic infrasound waveform inversions that use 3‐D Green’s functions for a variety of idealized and realistic deployment scenarios using both a flat plane and Yasur volcano, Vanuatu as examples. We investigate the ability of various scenarios to retrieve the input source functions and relative amplitudes for monopole and multipole (monopole and dipole) inversions. Infrasound waveform inversions appear to be a robust method to quantify mass flow rates from simple sources (monopole) using deployments of infrasound sensors placed around a source, but care should be taken when analyzing and interpreting results from more complex acoustic sources (multipole) that have significant directional components. In the examples we consider the solution is stable for monopole inversions with a signal‐to‐noise ratio greater than five and the dipole component is small. For most scenarios investigated, the vertical dipole component of the multipole explosion source is poorly constrained and can impact the ability to recover the other source term components. Because multipole inversions are ill‐posed for many deployments, a low residual does not necessarily mean the proper source vector has been recovered. Synthetic studies can help investigate the limitations and place bounds on information that may be missing using monopole and multipole inversions for potentially directional sources.

 

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.

 

Abstract Infrasound data from arrays can be used to detect, locate, and quantify a variety of natural and anthropogenic sources from local to remote distances. However, many array processing methods use a single broad frequency range to process the data, which can lead to signals of interest being missed due to the choice of frequency limits or simultaneous clutter sources. We introduce a new open-source Python code that processes infrasound array data in multiple sequential narrow frequency bands using the least-squares approach. We test our algorithm on a few examples of natural sources (volcanic eruptions, mass movements, and bolides) for a variety of array configurations. Our method reduces the need to choose frequency limits for processing, which may result in missed signals, and it is parallelized to decrease the computational burden. Improvements of our narrow-band least-squares algorithm over broad-band least-squares processing include the ability to distinguish between multiple simultaneous sources if distinct in their frequency content (e.g., microbarom or surf vs. volcanic eruption), the ability to track changes in frequency content of a signal through time, and a decreased need to fine-tune frequency limits for processing. We incorporate a measure of planarity of the wavefield across the array (sigma tau, στ) as well as the ability to utilize the robust least trimmed squares algorithm to improve signal processing and insight into array performance. Our implementation allows for more detailed characterization of infrasound signals recorded at arrays that can improve monitoring and enhance research capabilities. 

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. 

Infrasound data are routinely used to detect and locate volcanic and other explosions, using both arrays and single sensor networks. However, at local distances (<15 km) topography often complicates acoustic propagation, resulting in inaccurate acoustic travel times leading to biased source locations when assuming straight-line propagation. Here we present a new method, termed Reverse Time Migration-Finite-Difference Time Domain (RTM-FDTD), that integrates numerical modeling into the standard RTM back-projection process. Travel time information is computed across the entire potential source grid via FDTD modeling to incorporate the effects of topography. The waveforms are then back-projected and stacked at each grid point, with the stack maximum corresponding to the likely source. We apply our method to three volcanoes with different network configurations, source-receiver distances, and topography. At Yasur Volcano, Vanuatu, RTM-FDTD locates explosions within ∼20 m of the source and differentiates between multiple vents. RTM-FDTD produces a more accurate location for the two Yasur subcraters than standard RTM and doubles the number of detected events. At Sakurajima Volcano, Japan, RTM-FDTD locates the source within 50 m of the active vent despite notable topographic blocking. The RTM-FDTD location is similar to that from the Time Reversal Mirror method, but is more computationally efficient. Lastly, at Shishaldin Volcano, Alaska, RTM and RTM-FDTD both produce realistic source locations (<50 m) for ground-coupled airwaves recorded on a four-station seismic network. We show that RTM is an effective method to detect and locate infrasonic sources across a variety of scenarios, and by integrating numerical modeling, RTM-FDTD produces more accurate source locations and increases the detection capability.