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.
- NSF-PAR ID:
- 10349819
- Date Published:
- Journal Name:
- Seismological Research Letters
- ISSN:
- 0895-0695
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
SUMMARY -
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; https://doi.org/10.25928/bgrseis_bblf-ifsd, Hupe et al., 2021a), two so-called microbarom frequency sets – covering both the lower (0.15–0.35 Hz, 15 min; https://doi.org/10.25928/bgrseis_mblf-ifsd, Hupe et al., 2021b) and a higher (0.45–0.65 Hz, 15 min; https://doi.org/10.25928/bgrseis_mbhf-ifsd, 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 (https://doi.org/10.25928/bgrseis_bbhf-ifsd, 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
-
ABSTRACT Earthquakes generate infrasound in multiple ways. Acoustic coupling at the surface from vertical seismic velocity, termed local infrasound, is often recorded by infrasound sensors but has seen relatively little study. Over 140 infrasound stations have recently been deployed in Alaska. Most of these stations have single sensors, rather than arrays, and were originally installed as part of the EarthScope Transportable Array. The single sensor nature, paucity of ground-truth signals, and remoteness makes evaluating their data quality and utility challenging. In addition, despite notable recent advances, infrasound calibration and frequency response evaluation remains challenging, particularly for large networks and retrospective analysis of sensors already installed. Here, we examine local seismoacoustic coupling on colocated seismic and infrasound stations in Alaska. Numerous large earthquakes across the region in recent years generated considerable vertical seismic velocity and local infrasound that were recorded on colocated sensors. We build on previous work and evaluate the full infrasound station frequency response using seismoacoustic coupled waves. By employing targeted signal processing techniques, we show that a single seismometer may be sufficient for characterizing the response of an entire nearby infrasound array. We find that good low frequency (<1 Hz) infrasound station response estimates can be derived from large (Mw>7) earthquakes out to at least 1500 km. High infrasound noise levels at some stations and seismic-wave energy focused at low frequencies limit our response estimates. The response of multiple stations in Alaska is found to differ considerably from their metadata and are related to improper installation and erroneous metadata. Our method provides a robust way to remotely examine infrasound station frequency response and examine seismoacoustic coupling, which is being increasingly used in airborne infrasound observations, earthquake magnitude estimation, and other applications.more » « less
-
more » « less
Abstract Infrasound (low‐frequency acoustic waves) has proven useful to detect and characterize subaerial volcanic activity, but understanding the infrasonic source during sustained eruptions is still an area of active research. Preliminary comparison between acoustic eruption spectra and the jet noise similarity spectra suggests that volcanoes can produce an infrasonic form of jet noise from turbulence. The jet noise similarity spectra, empirically derived from audible laboratory jets, consist of two noise sources: large‐scale turbulence (LST) and fine‐scale turbulence (FST). We fit the similarity spectra quantitatively to eruptions of Mount St. Helens in 2005, Tungurahua in 2006, and Kīlauea in 2018 using nonlinear least squares fitting. By fitting over a wide infrasonic frequency band (0.05–10 Hz) and restricting the peak frequency above 0.15 Hz, we observe a better fit during times of eruption versus non‐eruptive background noise. Fitting smaller overlapping frequency bands highlights changes in the fit of LST and FST spectra, which aligns with observed changes in eruption dynamics. Our results indicate that future quantitative spectral fitting of eruption data will help identify changes in eruption source parameters such as velocity, jet diameter, and ash content which are critical for effective hazard monitoring and response.
-
more » « less
SUMMARY Protracted episodes of 0.5–7 Hz pre-eruptive volcanic tremor (PVT) are common at active stratovolcanoes. Reliable links to processes related to magma movement consequently enable a potential to use properties of PVT as diagnostic eruptive precursors. A challenging feature of PVT is that generic spectral and amplitude properties of the signal evolve similarly, independent of widely varying volcano structures and conduit geometries on which most physical models rely. The ‘magma wagging’ model introduced in Jellinek & Bercovici (2011) and extended by Bercovici et al. (2013), Liao et al. and Liao & Bercovici (2018) makes progress because it depends on magma dynamics that are only weakly sensitive to volcano architecture: The flow of gas through a permeable foamy annulus of gas bubbles excites, modulates and maintains a wagging oscillation of a central magma column rising in an erupting conduit. ‘Magma wagging’ and resulting PVT are driven through an energy transfer from a ‘Bernoulli mode’ related to azimuthal variations in annular gas flow speeds. Consistent with observations, spectral and amplitude properties of PVT are predicted to evolve before an eruption as the width of the annulus decreases with increased gas fluxes. To confirm this critical Bernoulli-to-wagging energy transfer we use extensive experiments and restricted numerical simulations on wagging oscillations excited on analogue viscoelastic columns by annular air flows. We also explore sensitivities of the spatial and temporal characters of wagging to asymmetric annular air flows that are intractable in the existing magma wagging model and expected to occur in nature with spatial variations in annulus permeability. From high-resolution time-series of linear and orbital displacements of analogue column tops and time-series of axial deflections and accelerations of the column centre line, we characterize the excitation, evolution, and steady-state oscillations in unprecedented detail over a broad range of conditions. We show that the Bernoulli mode corresponds to the timescale for the buildup of axial elastic bending stresses in response to pressure variations related to air flows over the heights of columns. We identify three distinct wagging modes: (i) rotational (cf. Liao et al. 2018); (ii) mixed-mode and (iii) chaotic. Rotational modes are favoured for symmetric, high intensity forcing and a maximal delivery of mechanical energy to the fundamental magma wagging mode. Mixed-mode oscillations regimes are favoured for a symmetric, intermediate intensity forcing. Chaotic modes, involving the least efficient delivery of energy to the fundamental mode, occur for asymmetric forcing and where the intensity of imposed airflow is low. Numerical simulations also show that where forcing frequencies are comparable to a natural mode of free oscillation, power delivered by peripheral air flows is concentrated at the lowest frequency fundamental mode generally and spread among higher frequency natural modes where air pressure and column elastic forces are comparable. Our combined experimental and numerical results make qualitative predictions for the evolution of the character of volcanic tremor and its expression in seismic or infrasound arrays during natural events that is testable in field-based studies of PVT and syn-eruptive volcanic tremor.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1785/0220220042
