An official website of the United States government
Abstract We present the transverse coherence minimization method (TCM)—an approach to estimate the back-azimuth of infrasound signals that are recorded on an infrasound microphone and a colocated three-component seismometer. Accurate back-azimuth information is important for a variety of monitoring efforts, but it is currently only available for infrasound arrays and for seismoacoustic sensor pairs separated by 10 s of meters. Our TCM method allows for the analysis of colocated sensor pairs, sensors located within a few meters of each other, which may extend the capabilities of existing seismoacoustic networks and supplement operating infrasound arrays. This approach minimizes the coherence of the transverse component of seismic displacement with the infrasound wave to estimate the infrasound back-azimuth. After developing an analytical model, we investigate seismoacoustic signals from the August 2012 Humming Roadrunner experiment and the 26 May 2021 eruption of Great Sitkin Volcano, Alaska, U.S.A., at the ranges of 6.5–185 km from the source. We discuss back-azimuth estimates and potential sources of deviation (1°–15°), such as local terrain effects or deviation from common analytical models. This practical method complements existing seismoacoustic tools and may be suitable for routine application to signals of interest.
more »
« less
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
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
- Award ID(s):
- 1847736
- PAR ID:
- 10572337
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 241
- Issue:
- 1
- ISSN:
- 0956-540X
- Format(s):
- Medium: X Size: p. 268-290
- Size(s):
- p. 268-290
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Mount Michael stratovolcano, South Sandwich Islands is extremely remote and challenging to observe, but eruptive activity has been sporadically observed since 1820 and captured by satellite methods since 1989. We identify long‐range infrasound signals recorded by the International Monitoring System attributable to episodes of persistent eruptive activity at Mount Michael. Analysis of multi‐year (2004–2020) infrasound array data at station IS27, Antarctica (range 1,672 km) reveals candidate signals especially from May 2005 to January 2008 and from May 2016 to April 2018. By combining ray‐tracing with empirical climatologies and atmospheric specifications, we show that systematic variations in the observed backazimuth of the signals (at IS27) are broadly consistent with annual variability in stratospheric propagation conditions for a source at Mount Michael. Observed signal amplitudes combined with transmission loss estimates are consistent with moderate explosive eruption. We highlight a selection of infrasound signals that correspond to satellite observation of eruptions.more » « less
-
null (Ed.)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.more » « less
-
Abstract We recorded a MWR3.6 earthquake in Idaho (USA) on 7 April 2020 with a network of six three‐element infrasound arrays and co‐located broadband seismometers situated within 25 km of the hypocenter. Infrasound array processing is used to identify the arrival of seismic‐to‐atmospheric coupled phases and as much as 90 s of infrasound coda. Apparent velocities ranging from seismic speeds to subhorizontal atmospheric sound speeds are attributed to a superposition of coincident waves arriving at the arrays. We find that the arriving infrasound originates from a broad range of back azimuths that deviates from epicentral back azimuth and indicates the ubiquity of secondary radiators for this relatively small earthquake. Secondary radiators, which often locate in regions of elevated topography, are identified using backprojections and earthquake initiation time. Analysis of infrasound sources from proximal earthquakes can be used to map ground shaking distributions, which are important for assessment of earthquake hazards.more » « less
-
Abstract Volcano infrasound data contain a wealth of information about eruptive patterns, for which machine learning (ML) is an emerging analysis tool. Although global catalogs of labeled infrasound events exist, the application of supervised ML to local (<15 km) volcano infrasound signals has been limited by a lack of robust labeled datasets. Here, we automatically generate a labeled dataset of >7500 explosions recorded by a five-station infrasound network at the highly active Yasur Volcano, Vanuatu. Explosions are located via backprojection and associated with one of Yasur’s two summit subcraters. We then apply a supervised ML approach to classify the subcrater of origin. When trained and tested on data from the same station, our chosen algorithm is >95% accurate; when training and testing on different stations, accuracy drops to about 75%. The choice of waveform features provided to the algorithm strongly influences classification performance.more » « less
