skip to main content

Title: Infrasound single-channel noise reduction: application to detection and localization of explosive volcanism in Alaska using backprojection and array processing

Infrasound sensors are deployed in a variety of spatial configurations and scales for geophysical monitoring, including networks of single sensors and networks of multisensor infrasound arrays. Infrasound signal detection strategies exploiting these data commonly make use of intersensor correlation and coherence (array processing, multichannel correlation); network-based tracking of signal features (e.g. reverse time migration); or a combination of these such as backazimuth cross-bearings for multiple arrays. Single-sensor trace-based denoising techniques offer significant potential to improve all of these various infrasound data processing strategies, but have not previously been investigated in detail. Single-sensor denoising represents a pre-processing step that could reduce the effects of ambient infrasound and wind noise in infrasound signal association and location workflows. We systematically investigate the utility of a range of single-sensor denoising methods for infrasound data processing, including noise gating, non-negative matrix factorization, and data-adaptive Wiener filtering. For the data testbed, we use the relatively dense regional infrasound network in Alaska, which records a high rate of volcanic eruptions with signals varying in power, duration, and waveform and spectral character. We primarily use data from the 2016–2017 Bogoslof volcanic eruption, which included multiple explosions, and synthetics. The Bogoslof volcanic sequence provides an opportunity to investigate regional infrasound detection, association, and location for a set of real sources with varying source spectra subject to anisotropic atmospheric propagation and varying noise levels (both incoherent wind noise and coherent ambient infrasound, primarily microbaroms). We illustrate the advantages and disadvantages of the different denoising methods in categories such as event detection, waveform distortion, the need for manual data labelling, and computational cost. For all approaches, denoising generally performs better for signals with higher signal-to-noise ratios and with less spectral and temporal overlap between signals and noise. Microbaroms are the most globally pervasive and repetitive coherent ambient infrasound noise source, with such noise often referred to as clutter or interference. We find that denoising offers significant potential for microbarom clutter reduction. Single-channel denoising of microbaroms prior to standard array processing enhances both the quantity and bandwidth of detectable volcanic events. We find that reduction of incoherent wind noise is more challenging using the denoising methods we investigate; thus, station hardware (wind noise reduction systems) and site selection remain critical and cannot be replaced by currently available digital denoising methodologies. Overall, we find that adding single-channel denoising as a component in the processing workflow can benefit a variety of infrasound signal detection, association, and location schemes. The denoising methods can also isolate the noise itself, with utility in statistically characterizing ambient infrasound noise.

more » « less
Award ID(s):
1847736 1614855
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Geophysical Journal International
Page Range / eLocation ID:
p. 1684-1712
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The International Monitoring System (IMS) infrasound network has been established to detect nuclear explosions and other signals of interest embedded in the station‐specific ambient noise. The ambient noise can be separated into coherent infrasound (e.g., real infrasonic signals) and incoherent noise (such as that caused by wind turbulence). Previous work statistically and systematically characterized coherent infrasound recorded by the IMS. This paper expands on this analysis of the coherent ambient infrasound by including updated IMS data sets with data up to the end of 2020 for all 53 of the currently certified IMS infrasound stations using an updated configuration of the Progressive Multi‐Channel Correlation (PMCC) method. This paper presents monthly station‐dependent reference curves for the back azimuth, trace velocity, and root mean squared amplitude, which provides a means to determine the deviation from the nominal monthly behavior. In addition, a daily Ambient Noise Stationarity (ANS) factor based on deviations from the reference curves is determined for a quick reference to the coherent signal quality compared to the nominal situations. Newly presented histograms provide a higher resolution spectrum, including the observations of the microbarom peak, as well as additional peaks reflecting station‐dependent environmental noise. The aim of these reference curves is to identify periods of suboptimal operation (e.g., nonoperational sensor) or instances of strong abnormal signals of interest.

    more » « less
  2. null (Ed.)
    Abstract Erosion, hydrothermal activity, and magmatism at volcanoes can cause large and unexpected mass wasting events. Large fluidized debris flows have occurred within the past 6000 yr at Mount Adams, Washington, and present a hazard to communities downstream. In August 2017, we began a pilot experiment to investigate the potential of infrasound arrays for detecting and tracking debris flows at Mount Adams. We deployed a telemetered four-element infrasound array (BEAR, 85 m aperture), ~11 km from a geologically unstable area where mass wasting has repeatedly originated. We present a preliminary analysis of BEAR data, representing a survey of the ambient infrasound and noise environment at this quiescent stratovolcano. Array processing reveals near continuous and persistent infrasound signals arriving from the direction of Mount Adams, which we hypothesize are fluvial sounds from the steep drainages on the southwest flank. We interpret observed fluctuations in the detectability of these signals as resulting from a combination of (1) wind-noise variations at the array, (2) changes in local infrasound propagation conditions associated with atmospheric boundary layer variability, and (3) changing water flow speeds and volumes in the channels due to freezing, thawing, and precipitation events. Suspected mass movement events during the study period are small (volumes <105  m3 and durations <2 min), with one of five visually confirmed events detected infrasonically at BEAR. We locate this small event, which satellite imagery suggests was a glacial avalanche, using three additional temporary arrays operating for five days in August 2018. Events large enough to threaten downstream communities would likely produce stronger infrasonic signals detectable at BEAR. In complement to recent literature demonstrating the potential for infrasonic detection of volcano mass movements (Allstadt et al., 2018), this study highlights the practical and computational challenges involved in identifying signals of interest in the expected noisy background environment of volcanic topography and drainages. 
    more » « less
  3. Abstract. The International Monitoring System (IMS) was established in the late 1990sfor verification of the Comprehensive Nuclear-Test-Ban Treaty (CTBT). Uponcompletion, 60 infrasound stations distributed over the globe will monitorthe Earth's atmosphere for low-frequency pressure waves. In this study, wepresent advanced infrasound data products of the 53 currently certified IMSinfrasound stations for atmospheric studies and civilian applications. Forthis purpose, 18 years of raw IMS infrasound waveform data (2003–2020) were reprocessed using the Progressive Multi-Channel Correlation (PMCC) method. A one-third octave frequency band configuration between 0.01 and 4 Hz was chosen to run this array-processing algorithm which detects coherent infrasound waves within the background noise. From the comprehensive detection lists, four products were derived for each of the certified 53 IMS infrasound stations. The four products cover different frequency ranges and are provided at the following different temporal resolutions: a very low-frequency set (0.02–0.07 Hz, 30 min;, Hupe et al., 2021a), two so-called microbarom frequency sets – covering both the lower (0.15–0.35 Hz, 15 min;, Hupe et al., 2021b) and a higher (0.45–0.65 Hz, 15 min;, Hupe et al., 2021c) part – named after the dominant ambient noise of interacting ocean waves that are quasi-continuously detected at IMS stations, and observations with center frequencies of 1 to 3 Hz (5 min), called the high-frequency product (, Hupe et al., 2021d). Within these frequency ranges and time windows, the dominant repetitive signal directions are summarized. Along with several detection parameters, calculated quantities for assessing the relative quality of the products are provided. The validity of the data products is demonstrated through example case studies of recent events that produced infrasound detected at IMS infrasound stations and through a global assessment and summary of the products. The four infrasound data products cover globally repeating infrasound sources such as ocean ambient noise or persistently active volcanoes, which have previously been suggested as sources for probing the winds in the middle atmosphere. Therefore, our infrasound data products open up the IMS observations also to user groups who do not have unconstrained access to IMS data or who are unfamiliar with infrasound data processing using the PMCC method. These types of data products could potentially serve as a basis for volcanic eruption early warning systems in the future. 
    more » « less
  4. 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. 
    more » « less
  5. Abstract

    Since the 1919 foundation of the International Association of Volcanology and Chemistry of the Earth’s Interior (IAVCEI), the fields of volcano seismology and acoustics have seen dramatic advances in instrumentation and techniques, and have undergone paradigm shifts in the understanding of volcanic seismo-acoustic source processes and internal volcanic structure. Some early twentieth-century volcanological studies gave equal emphasis to barograph (infrasound and acoustic-gravity wave) and seismograph observations, but volcano seismology rapidly outpaced volcano acoustics and became the standard geophysical volcano-monitoring tool. Permanent seismic networks were established on volcanoes (for example) in Japan, the Philippines, Russia, and Hawai‘i by the 1950s, and in Alaska by the 1970s. Large eruptions with societal consequences generally catalyzed the implementation of new seismic instrumentation and led to operationalization of research methodologies. Seismic data now form the backbone of most local ground-based volcano monitoring networks worldwide and play a critical role in understanding how volcanoes work. The computer revolution enabled increasingly sophisticated data processing and source modeling, and facilitated the transition to continuous digital waveform recording by about the 1990s. In the 1970s and 1980s, quantitative models emerged for long-period (LP) event and tremor sources in fluid-driven cracks and conduits. Beginning in the 1970s, early models for volcano-tectonic (VT) earthquake swarms invoking crack tip stresses expanded to involve stress transfer into the wall rocks of pressurized dikes. The first deployments of broadband seismic instrumentation and infrasound sensors on volcanoes in the 1990s led to discoveries of new signals and phenomena. Rapid advances in infrasound technology; signal processing, analysis, and inversion; and atmospheric propagation modeling have now established the role of regional (15–250 km) and remote (> 250 km) ground-based acoustic systems in volcano monitoring. Long-term records of volcano-seismic unrest through full eruptive cycles are providing insight into magma transport and eruption processes and increasingly sophisticated forecasts. Laboratory and numerical experiments are elucidating seismo-acoustic source processes in volcanic fluid systems, and are observationally constrained by increasingly dense geophysical field deployments taking advantage of low-power, compact broadband, and nodal technologies. In recent years, the fields of volcano geodesy, seismology, and acoustics (both atmospheric infrasound and ocean hydroacoustics) are increasingly merging. Despite vast progress over the past century, major questions remain regarding source processes, patterns of volcano-seismic unrest, internal volcanic structure, and the relationship between seismic unrest and volcanic processes.

    more » « less