Search for: All records

Seismicity at restless volcanoes commonly features a variety of signal types reflecting both volcanotectonic and fluid-driven source processes. However, traditional catalogs of seismicity are often incomplete, especially concerning events with emergent onsets such as those driven by the dynamics of magmatic and hydrothermal fluids. The detection of all discrete events and continuous seismic tremors, regardless of the underlying source processes, would therefore improve the ability of monitoring agencies to forecast eruptions and mitigate their associated hazards. We present a workflow for generalized detection of seismic events based on the network covariance matrix (Seydoux et al., 2016). Our contributions enable the method to simultaneously detect continuous and short-duration (<∼10 s) events, provide information about the frequency content of the signals, and to refine the initial detection times by an order of magnitude (from window lengths of 75 to 7.5 s). We test the workflow on a 15-month record of seismicity with 23 stations at Mammoth Mountain, California (July 2012–October 2013) and detect 62% of long-period events and 94% of volcanotectonic events in the existing Northern California Earthquake Data Center catalog. In addition, ∼3000 events are not included in the catalog, and thousands of tremor signals are found. The method is suitable for near-real-time analysis of continuous waveforms and can provide a valuable supplement to existing algorithms to improve the completeness of catalogs used for monitoring volcanoes.

 

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.

 

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.

 

Abstract Records of pressure variations on seismographs were historically considered unwanted noise; however, increased deployments of collocated seismic and acoustic instrumentation have driven recent efforts to use this effect induced by both wind and anthropogenic explosions to invert for near-surface Earth structure. These studies have been limited to shallow structure because the pressure signals have relatively short wavelengths (<∼300 m). However, the 2022 eruption of Hunga Tonga–Hunga Ha’apai (also called “Hunga”) volcano in Tonga generated rare, globally observed, high-amplitude infrasound signals with acoustic wavelengths of tens of kilometers. In this study, we examine the acoustic-to-seismic coupling generated by the Hunga eruption across 82 Global Seismographic Network (GSN) stations and show that ground motion amplitudes are related to upper (0 to ∼5 km) crust material properties. We find high (>0.8) correlations between pressure and vertical component ground motion at 83% of the stations, but only 30% of stations show this on the radial component, likely due to complex tilt effects. We use average elastic properties in the upper 5.2 km from the CRUST1.0 model to estimate vertical seismic/acoustic coupling coefficients (SV/A) across the GSN network and compare these to recorded observations. We exclude many island stations from these comparisons because the 1° resolution of the CRUST1.0 model places a water layer below these stations. Our simple modeling can predict observed SV/A within a factor of 2 for 94% of the 51 non-island GSN stations with high correlations between pressure and ground motion. These results indicate that analysis of acoustic-to-seismic coupling from the eruption could be used to place additional constraints on crustal structure models at stations with collocated seismic and pressure sensors. Ultimately, this could improve tomographic imaging models, which rely on methods that are sensitive to local structure. 

Abstract Laterally directed explosive eruptions are responsible for multiple fatalities over the past decade and are an increasingly important volcanology problem. To understand the energy dynamics for these events, we collected field-scale explosion data from nine acoustic sensors surrounding a tiltable cannon as part of an exploratory experimental design. For each cannon discharge, the blast direction was varied systematically at 0°, 12°, and 24° from vertical, capturing acoustic wavefield directivity related to the tilt angle. While each event was similar in energy discharge potential, the resulting acoustic signal features were variable event-to-event, producing non-repetitious waveforms and spectra. Systematic features were observed in a subset of individual events for vertical and lateral discharges. For vertical discharges, the acoustic energy had a uniform radiation pattern. The lateral discharges showed an asymmetric radiation pattern with higher frequencies in the direction of the blast and depletion of those frequencies behind the cannon. Results suggest that, in natural volcanic systems, near-field blast directionality may be elucidated from acoustic sensors in absence of visual data, with implications for volcano monitoring and hazard assessment. Graphical Abstract 

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. 

Real-time monitoring is crucial to assess hazards and mitigate risks of sustained volcanic eruptions that last hours to months or more. Sustained eruptions have been shown to produce a low frequency (infrasonic) form of jet noise. We analyze the lava fountaining at fissure 8 during the 2018 Lower East Rift Zone eruption of Kīlauea volcano, Hawaii, and connect changes in fountain properties with recorded infrasound signals from an array about 500 m from the fountain using jet noise scaling laws and visual imagery. Video footage from the eruption reveals a change in lava fountain dynamics from a tall, distinct fountain at the beginning of June to a low fountain with a turbulent, out-pouring lava pond surrounded by a tephra cone by mid-June. During mid-June, the sound pressure level reaches a maximum, and peak frequency drops. We develop a model that uses jet noise scaling relationships to estimate changes in volcanic jet diameter and jet velocity from infrasound sound pressure levels and peak frequencies. The results of this model indicate a decrease in velocity in mid-June which coincides with the decrease in fountain height. Furthermore, the model results suggest an increase in jet diameter, which can be explained by the larger width of the fountain that resembles a turbulent lava pond compared to the distinct fountain at the beginning of June. The agreement between the infrasound-derived and visually observed changes in fountain dynamics suggests that jet noise scaling relationships can be used to monitor lava fountain dynamics using infrasound recordings. 

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. 

The 21–22 June 2019 eruption of Raikoke volcano, Russia, provided an opportunity to explore how spatial trends in volcanic lightning locations provide insights into pulsatory eruption dynamics. Using satellite-derived plume heights, we examine the development of lightning detected by Vaisala’s Global Lightning Dataset (GLD360) from eleven, closely spaced eruptive pulses. Results from one-dimensional plume modeling show that the eruptive pulses with maximum heights 9–16.5 km above sea level were capable of producing ice in the upper troposphere, which contributed variably to electrification and volcanic lightning. A key finding is that lightning locations not only followed the main dispersal direction of these ash plumes, but also tracked a lower-level cloud derived from pyroclastic density currents. We show a positive relationship between umbrella cloud expansion and the area over which lightning occurs (the ‘lightning footprint’). These observations suggest useful metrics to characterize ongoing eruptive activity in near real-time.