More Like this

1. Infrasonic directivity of monopole, dipole and bipole ground-surface reflected sources

   https://doi.org/10.1093/gji/ggaf180  

   Iezzi, Alexandra M; Matoza, Robin S; Opper, Emma V; Kim, Keehoon (June 2025, Geophysical Journal International)

   SUMMARY Infrasound (acoustic waves below 20 Hz) can be used to detect, locate and quantify activity in the atmosphere such as volcanic eruptions and anthropogenic explosions. Attempts to quantify volcanic eruption parameters such as exit velocity, plume height and mass flow rate using infrasound data depend strongly on assumptions of the acoustic source type. Infrasonic sources may produce omnidirectional or directional wavefields, while propagation effects, such as interaction with topography, can induce further wavefield directivity that is measured by field instrumentation. Limited sampling of these wavefields can hinder our ability to infer the underlying source, and thus our understanding of the eruption characteristics. Equivalent sources are often used to represent acoustic source mechanisms and resultant wavefields. In this study, we review equivalent acoustic sources as they pertain to infrasonic scale and wavelengths commonly encountered in very local ($<$5 km range) geophysical field deployments. We highlight the equivalent infrasonic bipole source that can be induced by ground-reflection of an elevated monopole; we are not aware of any prior infrasound studies that use the bipole source concept. We use analytical and numerical methods to explore source directivity of monopole, dipole and bipole ground-reflected sources at infrasonic frequencies as well as the additional directivity complications introduced by interactions with topography. We illustrate that for typical volcano-infrasound wavelengths, increasing height above the ground as well as increasing source frequency leads to increased wavefield directivity. Numerical modelling using a simple omnidirectional monopole source embedded in topography further illustrates that both horizontal and vertical infrasound directionality can be induced by topography at the distance scales appropriate for local volcano infrasound monitoring. Information summarized in this analytical and numerical exploration of infrasound directivity may be used to help guide future volcano-infrasound field deployments intended to estimate source parameters or quantify wavefield directivity. Analytic solutions for simple whole-space or half-space atmospheres provide useful formulations for planning or initially analysing geophysical field-scale experimental data; however, especially at very local distances from the source ($<$5 km), 3-D simulations are necessary to account for complex topography commonly encountered in volcano-infrasound applications. 

   more » « less
2. Evidence for near-source nonlinear propagation of volcano infrasound from Strombolian explosions at Yasur Volcano, Vanuatu

   https://doi.org/10.1007/s00445-022-01552-w  

   Maher, Sean P.; Matoza, Robin S.; Jolly, Arthur; de Groot-Hedlin, Catherine; Gee, Kent L.; Fee, David; Iezzi, Alexandra M. (April 2022, Bulletin of Volcanology)

   Abstract 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$$ ≤ 10 −3  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. 

   more » « less
3. High-rate very-long-period seismicity at Yasur volcano, Vanuatu: Source mechanism and decoupling from surficial explosions and infrasound

   https://doi.org/10.1093/gji/ggab533  

   Matoza, Robin S; Chouet, Bernard A; Jolly, Arthur D; Dawson, Phillip B; Fitzgerald, Rebecca H; Kennedy, Ben M; Fee, David; Iezzi, Alexandra M; Kilgour, Geoff N; Garaebiti, Esline; et al (January 2022, Geophysical Journal International)

   Summary 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 broadband seismic and infrasonic network coverage from 27 July to 3 August 2016. 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 to 10 sec. The VLP seismic events are “high-rate”, reoccurring near-continuously throughout the dataset with short interevent times (∼20 to 60 sec). 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 time-scale 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–1,000 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–1,000 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. 

   more » « less
4. Seismo-acoustic characterisation of the 2018 Ambae (Manaro Voui) eruption, Vanuatu

   https://doi.org/10.1007/s00445-021-01474-z  

   Park, Iseul; Jolly, Arthur; Matoza, Robin S.; Kennedy, Ben; Kilgour, Geoff; Johnson, Richard; Garaebiti, Esline; Cevuard, Sandrine (September 2021, Bulletin of Volcanology)

   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. 

   more » « less
5. Evaluating the temporal capability of empirical climatologies for rapid long-range volcanic infrasound propagation estimates using a multidecadal data set of persistent Vanuatu volcanic eruptions

   https://doi.org/10.1093/gji/ggaf027  

   De Negri, Rodrigo S; Matoza, Robin S; Hupe, Patrick; Le Pichon, Alexis; Rose, Kaelynn M; Cevuard, Sandrine; Niroa, John J (February 2025, Geophysical Journal International)

   SUMMARY Powerful infrasound (acoustic waves $$\lt $$20 Hz) can be produced by explosive volcanic eruptions. The long-range propagation capability, over hundreds to thousands of kilometers, of atmospheric infrasound motivates the development of regional or even global scale volcano-infrasound monitoring systems. Infrasound propagation paths are subject to spatiotemporal atmospheric dynamics, which lead to deviations in the direction-of-arrival (back-azimuth) observed at sensor arrays and contribute to source location uncertainty. Here, we further investigate the utility of empirical climatologies combined with 3-D ray tracing for providing first-order estimates of infrasound propagation paths and back-azimuth deviation corrections. The intended application is in scenarios requiring rapid or pre-computed infrasound propagation calculations, such as for a volcano-infrasound monitoring system. Empirical climatologies are global observationally based function fitting models of the atmosphere, representing robust predictors of the bulk diurnal to seasonal atmospheric variability. Infrasound propagation characteristics have previously been shown to have strong seasonal and diurnal components. At the International Monitoring System infrasound station IS22, New Caledonia, quasi-continuous multiyear infrasound array detections show oscillating azimuthal variations for arrivals from volcanoes in Vanuatu, including Yasur ($$\sim$$400 km range), Ambrym ($$\sim$$670 km range) and Lopevi ($$\sim$$650 km range). We perform 3-D ray tracing to model infrasound propagation from the Ambrym and Yasur volcano locations to IS22 every six hours (00:00, 06:00, 12:00 and 18:00 UTC) for every day of 2004 and 2019 for Ambrym and Yasur, respectively and evaluate the results as compared to the multiyear observations. We assess a variety of models and parametrizations, including both empirical climatologies and hybrid descriptions; range-independent and range-dependent atmospheric discretizations; and unperturbed and perturbed range-independent empirical climatologies. The hybrid atmospheric descriptions are composed of fifth generation reanalysis descriptions (ERA 5) from the European Centre for Medium-Range Weather Forecasts below $$\sim$$80 km altitude combined with empirical climatologies above. We propose and employ simple parametric perturbations to the empirical climatologies, which are designed to enhance the stratospheric duct and compensate for missing gravity wave perturbations not included in the climatologies, and thereby better match observations. We build year-long back-azimuth deviation interpolations from the simulations and compare them with three different multiyear array detection data sets from IS22 covering from 2003 up to 2022. Through a systematic comparison, we find that the range-independent empirical climatologies can capture bulk azimuth deviation variability and could thus be useful for rapid infrasound propagation calculation scenarios, particularly during favourable sustained propagation ducting conditions. We show that the hybrid models better describe infrasound propagation during periods of weak stratospheric ducting and during transient atmospheric changes such as stratospheric wind reversals. Overall, our results support the notion that climatologies, if perturbed to compensate for missing gravity wave structure, can improve rapid low-latency and pre-computed infrasound source discrimination and location procedures. 

   more » « less