Geology Reports⌕ Search

Geology topics

David Fee

Publications and source records attributed to David Fee.

At least 19 recordsLinked to original sources
USGS2026

Volcanic activity in Alaska and the Northern Mariana Islands in 2024—Summary of events and response of the Alaska Volcano Observatory

In 2024, the Alaska Volcano Observatory responded to eruptions, volcanic unrest or suspected unrest, increased seismicity, and other significant activity at 10 volcanic centers in Alaska and the Northern Mariana Islands. Eruptive activity in Alaska consisted of the eruption of a thick lava flow within the summit crater at Great Sitkin Volcano throughout the entire year and a small explosion at Atka volcano in March. Uplift, elevated seismicity, and the redevelopment of a summit meltwater lake was detected at Mount Spurr. A brief seismic swarm occurred at Mount Gareloi, while at Trident Volcano, Shishaldin Volcano, and Kanaga Volcano, seismicity declined to background levels, marking the end of unrest and volcanic activity for those volcanoes. Other activity that the Alaska Volcano Observatory responded to in 2024 included ash resuspension events at Mount Katmai, and a period of unrest at Mount Cleveland. Finally, at Ahyi seamount, in the Commonwealth of the Northern Marianas Islands, a plume of discolored ocean water observed in satellite data indicated underwater eruptive activity there.

Alaska↗
USGS2026

Rapid seismic and infrasound assessment of large landslides: A case study from Denali National Park and Preserve (Alaska)

Large, rapid landslides are a global hazard that can occur in remote, mountainous areas. Eyewitness reports of landslides and satellite imagery can often be limited or delayed, particularly during inclement weather. However, landslide-generated seismic and infrasound (low-frequency atmospheric sound) waves can be remotely detected in near real-time. This information can significantly expedite characterization and possible landslide response activities. Here, we highlight these capabilities using a > 4 million m³ ice–rock avalanche in Denali National Park and Preserve (Alaska). This event was detected via a landslide-specific seismic location and volume estimation algorithm deployed in Alaska, and — notably — by standard earthquake monitoring systems. Following rapid detection of this event, we combined its seismic and infrasound dataset with optical, synthetic aperture radar, and oblique aerial imagery, multitemporal digital elevation models, and a numerical flow model to reconstruct its failure timeline and dynamics. We apply array processing to infrasound signals traveling > 250 km and find that two precursory events occurred minutes prior to the main failure. We use long-period seismic signals to infer the force exerted by the landslide on the Earth and constrain the rheological parameters of our numerical flow simulation with this result and deposit morphology. The main failure produced a steeply-dipping impulsive initial downward force and reached speeds exceeding 60 m/s. This impulsive force generated relatively strong seismic body waves, which contributed to the earthquake system detection. This large, remote Alaska landslide underscores the key value of seismic and infrasound analysis for rapid landslide assessment and motivates efforts to further operationalize these approaches.

Alaska↗
USGS2025

The morphological response to the 13 May 2022 major explosive event at Stromboli volcano, Italy, characterised by high-resolution UAS surveys

Since 2019, the frequency of major explosive eruptions at Stromboli volcano (Italy) has increased, heightening the exposure of population and scientists to the hazards posed by ejecta. Morphological changes can directly alter the hazard potential associated with these phenomena. Here, we present a quantitative morphological analysis of changes of the crater terrace area linked to the 13 May 2022 major explosive event. High resolution (2.5 cm pixel -1 ) aerial imagery was acquired by unoccupied aircraft systems 2 hours before and 19 hours after the event. The 13 May 2022 major explosive event consisted of a minimum of seven explosions from four vents located in the south-central crater area. The opportune timing of this campaign enabled the quantification of morphological changes at Stromboli related to a single major explosive event at high temporal and spatial resolution. A total of 12.7 × 10 3 m 3 was excavated and 5.5 × 10 3 m 3 deposited. Via the mapping and classification of bomb distributions we observe that angular blocks make up the largest fraction of ballistics >0.2 m, from which we infer a strong interaction with wall rock and/or fragmentation of solidified plugs in the shallow plumbing system. The morphological changes observed provide valuable constraints on how much material is displaced, and the shift in location and the number of active vents during major explosive events at Stromboli.

Stromboli volcano↗
USGS2025

A generalized deep learning model to detect and classify volcano seismicity

Volcano seismicity is often detected and classified based on its spectral properties. However, the wide variety of volcano seismic signals and increasing amounts of data make accurate, consistent, and efficient detection and classification challenging. Machine learning (ML) has proven very effective at detecting and classifying tectonic seismicity, particularly using Convolutional Neural Networks (CNNs) and leveraging labeled datasets from regional seismic networks. Progress has been made applying ML to volcano seismicity, but efforts have typically been focused on a single volcano and are often hampered by the limited availability of training data. We build on the method of Tan et al. [2024] ( 10.1029/2024JB029194 ) to generalize a spectrogram-based CNN termed the VOlcano Infrasound and Seismic Spectrogram Neural Network ( VOISS-Net ) to detect and classify volcano seismicity at any volcano. We use a diverse training dataset of over 270,000 spectrograms from multiple volcanoes: Pavlof, Semisopochnoi, Tanaga, Takawangha, and Redoubt volcanoes\replaced (Alaska, USA); Mt. Etna (Italy); and Kīlauea, Hawai`i (USA). These volcanoes present a wide range of volcano seismic signals, source-receiver distances, and eruption styles. Our generalized VOISS-Net model achieves an accuracy of 87 % on the test set. We apply this model to continuous data from several volcanoes and eruptions included within and outside our training set, and find that multiple types of tremor, explosions, earthquakes, long-period events, and noise are successfully detected and classified. The model occasionally confuses transient signals such as earthquakes and explosions and misclassifies seismicity not included in the training dataset (e.g. teleseismic earthquakes). We envision the generalized VOISS-Net model to be applicable in both research and operational volcano monitoring settings.

Volcanica↗
USGS2025

Identifying precursors and tracking pulses of magma ascent in multidisciplinary data during the 2018–2023 phreatomagmatic eruption at Semisopochnoi Island, Alaska

The 2018–2023 phreatomagmatic eruptions at Semisopochnoi Island, Alaska produced abundant long-period (LP) seismicity, harmonic and broadband tremor, and explosion signals over several well-monitored periods of eruption and quiescence. The corresponding dataset provides an excellent opportunity to investigate precursory and syn -eruptive geophysical signals of long-lived phreatomagmatic eruptions using multiparameter observations. We generated explosion and LP event catalogs through novel implementations of the REDPy ( Hotovec-Ellis, 2024 ) repeating event detector in mid-2021 following a network upgrade and the onset of a new phase of the eruption. The hundreds of detected explosions show a high degree of infrasound waveform similarity over more than a year, indicating a repeating source mechanism likely associated with explosive magma-water interaction. The seismic LP catalog shows that events began over a month prior to renewed explosive activity at the beginning of August 2021, and that lower frequency index (FI) LPs were generated in the week prior to the onset of explosions. We applied a recently developed machine learning tool (VOISS-Net, Tan et al., 2024 ) to catalog abundant broadband and harmonic seismic tremor recorded before and during the renewed explosive activity, along with LPs and explosions. The tremor catalogs complement the LP and explosion catalogs by filling out the seismic sequence with the dominant signal types. Together, these catalogs reveal a seismic sequence of renewed unrest that started with several weeks of LP events, followed by LPs with lower FI values and harmonic tremor in the days prior to explosive activity, and finally the onset of discrete explosions and broadband eruption tremor. We interpret this sequence as the ascent of a new pulse of magma that first interacted with the hydrothermal/groundwater system to produce LPs, followed by harmonic tremor, and that ultimately drove explosive magma-water interactions and periods of continuous ash emissions. The 2021 seismic sequence, in combination with long-term records of satellite SO 2 emissions, deformation from interferometric synthetic aperture radar (InSAR) analysis, ash sample analysis, infrasound, and volcano tectonic seismicity, allows us to interpret the entire 9-year period of unrest and eruption that began with an intrusion and earthquake swarm in 2014.

Alaska↗
USGS2025

Vent-specific unrest at Pavlof Volcano, Alaska: Insights from multidisciplinary data

Pavlof Volcano, a frequently active volcano in the Aleutian arc, has erupted six times between 2007 and 2022. Many Pavlof Volcano eruptions lack detectable precursory seismicity or ground deformation, making them challenging to forecast. Since 2007, these eruptions have ranged from Strombolian to Subplinian (VEI 2–3), with their characteristics seemingly related to their active vent location. Summit vent eruptions tend to be abrupt and explosive, while southeast flank vent eruptions tend to be more effusive and longer-lived. To better understand Pavlof Volcano’s plumbing system and improve eruption forecasts, we integrate geophysical, petrologic, and satellite-based thermal infrared and gas data to conduct a multidisciplinary data synthesis for eruptions between 2007 and 2022. Notably, a recently developed volcano seismicity detection model reveals unique pre- and syn-eruptive seismic tremor regimes associated with each vent system. The summit vent eruptions show greater tremor diversity and resurgent eruptive behavior, whereas the southeast flank vent eruptions show a steady tremor transition over a single eruptive phase. We use local infrasound data and air-to-ground coupled waves to build explosion catalogs, revealing gas-rich explosions during the 2021–2022 eruption from the summit vent, which was previously thought to be inactive during that eruption. Whole-rock composition and ash analyses indicate a near-uniform basaltic andesite composition spanning multiple eruptions, and a consistent anomalous presence of altered olivines in erupted material. We use these findings to propose a new conceptual model for Pavlof Volcano’s plumbing system: a shallow T-junction outlet controlling seismic tremor diversity and the partitioning of gas-charged magma; and an elongated, heated conduit transporting magma from depth and storing partly molten, residual magma between eruptions. We suggest that magma ascent rate and vent sealing modulate eruption style at Pavlof Volcano, where the T-junction directs rapidly ascending, gas-rich magma to the summit vent, and diverts slower rising, gas-poor magma towards the southeast flank vent.

Alaska↗
USGS2024

Infrasound for volcano monitoring

Introduction Volcanic eruptions produce acoustic waves when volcanic gases and hot material rapidly expand in the atmosphere. Volcanic activity can produce acoustic signals with a wide range of frequencies, from very long period (>10 seconds) to audible (>20 hertz [Hz]), but the most energetic band is typically in the infrasound from 0.5 to 20 Hz. Studies of volcanic infrasound and the deployment of infrasound for volcano monitoring have increased rapidly in the past two decades as sensors have improved and as analytical tools have become more widely available. Improved sensors and tools have led to a growing diversity of eruptive activity being recorded and characterized, from Hawaiian to Plinian eruption styles at scales from local to global (Johnson and Ripepe, 2011; Fee and Matoza, 2013). Infrasound sensors on volcanoes are most commonly deployed locally with seismic stations, and the combination of co-located seismic and infrasound is more useful for characterizing unrest and detecting changes in activity than either data stream alone (for example, Lyons and others, 2016; Fee and others, 2017a; Matoza and others, 2018). At local (<15 kilometers [km]) to regional (15–250 km) distances from volcanoes, arrays of infrasound sensors are commonly deployed to detect coherent signals, constrain the direction to the source, and provide information on eruption dynamics; thus, infrasound is well suited to regional monitoring of volcanoes when local sensor networks are not feasible. A common usage of infrasound data in an observatory is to provide rapid confirmation that an explosion has occurred (for example, Coombs and others, 2018), although near-real-time eruption intensity quantification is also possible (Fee and others, 2010a; Ripepe and others, 2018; fig. C1). Infrasound is well suited to this task because it is not affected by clouds or precipitation and can propagate long distances with little attenuation. However, wind and ocean noise also produce infrasound, and spatiotemporal variability in the atmosphere can affect the propagation of infrasound, so care must be taken when deploying, analyzing, and interpreting the data. In addition to detecting and monitoring explosive activity, investigations of infrasound records from eruptions help constrain source processes, which in turn enhance syneruptive forecasting capabilities (for example, Fee and others, 2017b; Lyons and others, 2019). The following is a description of the capabilities recommended for real-time monitoring of eruptive phenomena with infrasound. Infrasound is also beginning to be used for tracking hazardous surface flows that occur on volcanoes, including pyroclastic density currents (Ripepe and others, 2010), lahars (Johnson and Palma, 2015), debris flows (Marchetti and others, 2019), snow avalanches (Havens and others, 2014), and lava flows (Patrick and others, 2019). Please refer to the chapter on lahars (this volume; Thelen and others, 2024a) for more information on this application.

Scientific Investigations Report↗
USGS2024

From field station to forecast: Managing data at the Alaska Volcano Observatory

The Alaska Volcano Observatory (AVO) uses multidisciplinary data to monitor and study dozens of active and potentially active volcanoes. Here, we provide an overview of internally and externally generated data types, tools and resources used in their management, and challenges faced. Data sources include the following: (1) a multiparameter (seismic, infrasound, GNSS, web cameras) ground-based monitoring network that spans 3000 km and transmits data in real time; (2) a variety of satellite-borne sensors that provide information about surface change and volcanic emissions; (3) geologic and gas field campaigns; and (4) other external data products that provide situation awareness. Each data type requires distinct acquisition, processing, storage, visualization, and archiving approaches. AVO uses a variety of externally and internally developed tools to handle individual data types as well as multidisciplinary volcanological data. A primary tool is the Geologic Database of Information on Volcanoes in Alaska (GeoDIVA), which stores detailed, searchable information on more than 140 volcanoes and over 1000 eruptions and unrest events, including images, eruption descriptions, and geologic station and sample data, metadata, and analyses. It interacts with other internal tools that store monitoring reports and other operational records. Additional data management resources used by AVO assist with alarms and alerts, state-of-health monitoring, and multiparameter visualization. Requirements for 24/7 accessibility, the ever-expanding portfolio of data, and transitioning new tools from development to operations are all challenges faced by AVO and other volcano observatories. AVO strives to meet FAIR data practices and ensure that data are available to national and international community efforts using external repositories as well as those hosted by AVO and its parent institutions.

Alaska↗
USGS2024

2021 Volcanic activity in Alaska and the Commonwealth of the Northern Mariana Islands—Summary of events and response of the Alaska Volcano Observatory

In 2021, the Alaska Volcano Observatory responded to eruptions, volcanic unrest or suspected unrest, increased seismicity, and other significant activity at 15 volcanic centers in Alaska and the Commonwealth of the Northern Mariana Islands. Eruptive activity in Alaska consisted of repeated small, ash-producing, phreatomagmatic explosions from Mount Young on Semisopochnoi Island; an explosion at Great Sitkin Volcano followed by the eruption of a thick lava flow that filled and overflowed the summit crater; weak explosive activity and the eruption of small, channelized flows at Pavlof Volcano; and a short-lived eruption at Mount Veniaminof that produced ash emissions from an intracaldera cone, as well as lava flows confined to a melt pit in the ice mantling the cone’s flank. Mount Cleveland had a period of unrest, but no eruptive activity took place there. Anomalous seismicity was also detected at Atka volcanic complex, Mount Gareloi, and Davidof volcano. New warm springs opened and deposited mud at the summit and north base of Shrub mud volcano. Other activity of note in Alaska consisted of large ice and rock avalanches at Iliamna Volcano and Mount Spurr, ash resuspension events at Mount Katmai and Aniakchak Crater, and anomalous deformation at Mount Okmok that was consistent with a shallow intrusion of magma. In the Commonwealth of the Northern Marianas Islands, a brief, ash-producing eruption occurred at Mount Pagan.

Alaska↗
USGS2023

Back-azimuth estimation of air-to-ground coupled infrasound from transverse coherence minimization

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.

The Seismic Record↗
USGS2023

Volcanic earthquake catalog enhancement using integrated detection, matched-filtering, and relocation tools

Volcanic earthquake catalogs are an essential data product used to interpret subsurface volcanic activity and forecast eruptions. Advances in detection techniques (e.g., matched-filtering, machine learning) and relative relocation tools have improved catalog completeness and refined event locations. However, most volcano observatories have yet to incorporate these techniques into their catalog-building workflows. This is due in part to complexities in operationalizing, automating, and calibrating these techniques in a satisfactory way for disparate volcano networks and their varied seismicity. In an effort to streamline the integration of catalog-enhancing tools at the Alaska Volcano Observatory (AVO), we have integrated four popular open-source tools: REDPy, EQcorrscan, HypoDD, and GrowClust. The combination of these tools offers the capability of adding seismic event detections and relocating events in a single workflow. The workflow relies on a combination of standard triggering and cross-correlation clustering (REDPy) to consolidate representative templates used in matched-filtering (EQcorrscan). The templates and their detections are then relocated using the differential time methods provided by HypoDD and/or GrowClust. Our workflow also provides codes to incorporate campaign data at appropriate junctures, and calculate magnitude and frequency index for valid events. We apply this workflow to three datasets: the 2012–2013 seismic swarm sequence at Mammoth Mountain (California), the 2009 eruption of Redoubt Volcano (Alaska), and the 2006 eruption of Augustine Volcano (Alaska); and compare our results with previous studies at each volcano. In general, our workflow provides a significant increase in the number of events and improved locations, and we relate the event clusters and temporal progressions to relevant volcanic activity. We also discuss workflow implementation best practices, particularly in applying these tools to sparse volcano seismic networks. We envision that our workflow and the datasets presented here will be useful for detailed volcano analyses in monitoring and research efforts.

Alaska↗
USGS2023

UAS-Based Observations of Infrasound Directionality at Stromboli Volcano, Italy

Infrasound (low frequency sound waves) can be used to monitor and characterize volcanic eruptions. However, infrasound sensors are usually placed on the ground, thus providing a limited sampling of the acoustic radiation pattern that can bias source size estimates. We present observations of explosive eruptions from a novel uncrewed aircraft system (UAS)-based infrasound sensor platform that was strategically hovered near the active vents of Stromboli volcano, Italy. We captured eruption infrasound from short-duration explosions and jetting events. While potential vertical directionality was inconclusive for the short-duration explosion, we find that jetting events exhibit vertical sound directionality that was observed with a UAS close to vertical. This directionality would not have been observed using only traditional deployments of ground-based infrasound sensors, but is consistent with jet noise theory. This proof-of-concept study provides unique information that can improve our ability to characterize and quantify the directionality of volcanic eruptions and their associated hazards.

Stromboli Volcano↗
USGS2022

Analog field-scale acoustic study of volcanic eruption directivity using a tiltable liquid nitrogen-charged water cannon

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.

Earth, Planets and Space↗
USGS2022

Lava fountain jet noise during the 2018 eruption of fissure 8 of Kīlauea volcano

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.

Hawaii↗
USGS2022

Return from dormancy: Rapid inflation and seismic unrest driven by transcrustal magma transfer at Mt. Edgecumbe (L’´ux Shaa) Volcano, Alaska

In April 2022, a seismic swarm near Mt. Edgecumbe in southeast Alaska suggested renewed activity at this transform fault volcano, which was last active ≈800 years ago. Previously, thin rhyolitic tephras were deposited 5 and 4 ka. Satellite radar data from 2014 to 2022 resolves line-of-sight rapid inflation up to 7.1 cm/yr beginning in August 2018. Bayesian modeling suggests a transcrustal system of a deflating (−0.528 km 3 ) dipping sill at 20 km depth recharging a magma chamber at 10 km (0.222 km 3 ). A near-vertical conduit could capture the volume difference without noticeable surface deformation. Reanalyzed seismicity, recorded 25 km away, shows increases since July 2019. Magma ascent through ductile material and brittle strain release in a stressed overburden could explain the time delay. Cloud-native open data and workflows enabled discovery and analysis of this signal within days after going unnoticed for >3 years.

Alaska↗
USGS2022

Infrasound observations and constraints on the 2018 eruption of Kīlauea Volcano, Hawaii

The 2018 eruption of Kīlauea Volcano was a dynamic event involving explosions, collapses, and fountaining at multiple vents spread over tens of kilometers. The permanent infrasound network operated by the USGS Hawaiian Volcano Observatory (HVO) was well prepared to observe the collapse of the summit, and additional deployments permitted infrasound observations during fissuring in the lower East Rift Zone (LERZ). We provide a summary of infrasound observations, including lava lake spattering, collapses, explosions, rockfall, and lava fountaining, using seismicity and tilt at times to help constrain our interpretations. At the summit of Kīlauea Volcano, we document the process of partial caldera collapse and examine a set of “proto-collapse” events that precede the widely observed events but share many of the same qualities as the larger collapses. For the initial twelve collapse events, we compare the timing of collapse onset to other observations and illustrate the repeatable characteristics of the recorded waveforms and infrasound characteristics associated with each episode of caldera collapse. In the LERZ, we match the acoustic signals with visual observations, including fissure migration, explosions near fissures, and littoral explosions. Lastly, we document and discuss the performance of infrasound alarms during the 2018 Kīlauea eruption. In general, alarming became successful in detecting collapse events at the summit of the volcano after tuning and became a key discriminant in the initial determination of collapse events, especially when visual observations were not available.

Hawaii↗
USGS2022

Evidence for near-source nonlinear propagation of volcano infrasound from Strombolian explosions at Yasur Volcano, Vanuatu

Volcanic eruption source parameters may be estimated from acoustic pressure recordings dominant at infrasonic frequencies (< 20 Hz), yet uncertainties may be high due in part to poorly understood propagation dynamics. Linear acoustic propagation of volcano infrasound is commonly assumed, but nonlinear processes such as wave steepening may distort waveforms and obscure the sourcing process in recorded waveforms. Here we use a previously developed frequency-domain nonlinearity indicator to quantify spectral changes due to nonlinear propagation primarily in 80 signals from explosions at Yasur Volcano, Vanuatu. We find evidence for &#x2264; "> ≤ ≤ 10 −3 dB/m spectral energy transfer in the band 3–9 Hz for signals with amplitude on the order of several hundred Pa at 200–400 m range. The clarity of the nonlinear spectral signature increases with waveform amplitude, suggesting stronger nonlinear changes for greater source pressures. We observe similar results in application to synthetics generated through finite-difference wavefield simulations of nonlinear propagation, although limitations of the model complicate direct comparison to the observations. Our results provide quantitative evidence for nonlinear propagation that confirms previous interpretations made on the basis of qualitative observations of asymmetric waveforms.

Bulletin of Volcanology↗