Search for: All records

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. 

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. 

Abstract 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.