Search for: All records

Abstract The 2018 eruption of Sierra Negra volcano, Galápagos, Ecuador has provided new insights into the mechanisms of caldera resurgence, subsidence, and fissuring at basaltic shield volcanoes. Here, we integrate local (∼0.4 km) seismo‐acoustic records and regional (∼85 km) infrasound array data to present new observations of the 2018 Sierra Negra eruption with improved time and spatial resolutions. These observations include: air‐to‐ground coupling ∼2 hr before the time of the eruption onset, migration of the infrasound tremor from 22:54 June 26 to 12:31 June 27 UT (all times in UT), and persistent infrasound detections during the weeks between 5 July and 18 August from an area that does not coincide with the previously documented eruptive fissures. We interpret air‐to‐ground coupling as infrasound tremor generated in the nearby fissures before the main eruptive phase started, although ambiguity remains in interpreting a single seismic‐infrasonic sensor pair. The progressive location change of the infrasound tremor agrees with the migration of the eruption down the north flank of Sierra Negra at a rate of ∼0.15 ± 0.04 m/s. The weeks‐long persistent detections coincide with a region that has thermal anomalies, co‐eruptive deformation, lava fields, and geological features that could be interpreted as multiple lava tube skylights. Our observations and interpretations provide constraints on the mechanisms underlying fissure formation and magma emplacement at Sierra Negra. 

Abstract Popocatépetl is a highly active stratovolcano in central Mexico with recurrent activity of Vulcanian-type explosions and frequent degassing. The proximity of Popocatépetl volcano to Mexico City, one of the most populated cities in the world, demands continuous monitoring to achieve an adequate volcano risk assessment. We present an overview of the first high-dynamic-range and high-broadband (0.01–200 Hz; 400 Hz sampling rate) seismoacoustic network (PoPiNet), which we operated around Popocatépetl volcano from August 2021 to May 2022. Here, we show preliminary results of the explosions recorded in September 2021. We deployed five seismoacoustic stations within 4–25 km horizontal distance (range) from the vent. We identify infrasonic waveforms associated with tremor and explosions, with pressures ranging from 16 to 134 Pa and dominant frequencies between 0.2 and 5.0 Hz. The frequency content of the recorded signals at the closest stations to the volcano spans the sub-bass (20–60 Hz) and bass (60–250 Hz) ranges. The associated seismic signals of moderate explosions exhibit air-to-ground coupled waves with maximum coherence values at frequencies up to 5 and 25 Hz for the farthest and closest stations to the volcano, respectively. Conversely, we observe infrasound signal amplitudes from relatively small explosions reaching maximum pressures of 10 Pa that do not couple into the ground, even at the closest stations. These infrasound signals are associated with type-I long-period events as reported in previous investigations. The waveform consistency suggests repetitive and nondestructive sources beneath the volcano. 

Abstract 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. 

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. 

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. 

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 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. 

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 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.