skip to main content


Search for: All records

Creators/Authors contains: "Matoza, Robin S."

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. 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. 
    more » « less
  2. 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
  3. 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 
    more » « less
  4. 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
  5. 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. 
    more » « less
  6. 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
  7. 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. 
    more » « less
  8. 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