USGS Science⌕ Search

Geology topics

David Fee

Publications and source records attributed to David Fee.

At least 55 records · Page 3Linked to original sources

Constraints on eruption processes and event masses for the 2016–2017 eruption of Bogoslof volcano, Alaska, through evaluation of IASI satellite SO2 masses and complementary datasets

Bogoslof volcano, Alaska, experienced at least 70 explosive eruptions between 12 December 2016 and 31 August 2017. Due to its remote location and limited local monitoring network, this eruption was monitored and characterized primarily using remote geophysical and satellite techniques. SO 2 emissions from Bogoslof were persistently detected by the Infrared Atmospheric Sounding Interferometer (IASI) satellite sensors. Of Bogoslof’s 70 explosive events, 50% produced measurable SO 2 masses ranging from 0.1 to 21.5 kt, with a median and standard deviation of 0.7 ± 4.0 kt SO 2 , respectively. Here, we compare IASI-derived SO 2 masses from Bogoslof events to complementary geophysical datasets to provide insights into eruption source processes, namely the degree of seawater scrubbing of water-soluble SO 2 and variations in magma flux. Correlations with the number of lightning strokes and infrasound energy are expected to indicate magma-flux as a controlling process, while correlations with infrasound frequency index are expected to indicate variations in vent-water content as a controlling factor. These comparisons suggest that the measured SO 2 masses are primarily a function of eruption magnitude (degassed magma mass) and that scrubbing of SO 2 emissions by vent seawater may have exerted a minor effect on the observed SO 2 masses. SO 2 masses were combined with petrologic constraints on melt inclusion and matrix glass S concentrations to calculate degassed magma masses and volumes. The cumulative SO 2 -derived degassed magma mass and estimated volume (dense-rock equivalent) for the full Bogoslof eruption were found to be 2.8 × 10 10 kg and 9.3 × 10 6 m 3 , respectively. When individual event masses are compared against event masses calculated using an empirical plume-height method, a strong correlation is found ( R 2 = 0.83), with better than order-of-magnitude agreement in most cases. These estimates of eruption masses provide useful information on the magnitude, behavior, and associated hazards of the 2016–2017 eruption, and potentially future unrest at Bogoslof volcano.

Alaska↗

Characteristics of thunder and electromagnetic pulses from volcanic lightning at Bogoslof volcano, Alaska

We combine global detections of volcanic lightning with acoustic and hydroacoustic data to investigate novel indications of plume electrification in ground-based, geophysical data streams during the 2016–2017 eruption of Bogoslof volcano, Alaska. Such signals offer additional ways to diagnose the occurrence of volcanic lightning and confirm whether eruptive activity is producing significant amounts of ash. We discuss three signatures of lightning activity: volcanic thunder, electromagnetic pulses arising from lightning-induced voltages in cabling, and hydroacoustic signals associated with volcanic lightning. Observations of these phenomena provide additional insights into volcanic lightning activity and reveal several periods of electrical activity that were not otherwise detected during the Bogoslof eruption.

Alaska↗

Goals and development of the Alaska Volcano Observatory Seismic Network and application to forecasting and detecting volcanic eruptions

The Alaska Volcano Observatory (AVO) seismic network has been in operation since 1988 and during this time has grown from 29 to 217 seismic stations providing real-time monitoring of 32 active volcanoes in Alaska, as well as useful data for regional earthquake monitoring. Since 1988, AVO has detected 59 volcanic eruptions at Aleutian arc volcanoes, and 31 of these have been captured by local seismic instrumentation. As part of this monitoring effort, AVO has cataloged more than 120,000 earthquake hypocenters and magnitudes associated with volcanic processes throughout the arc. This high rate of volcanic activity provides an excellent opportunity to study seismicity associated with magmatic and eruptive processes and develop and refine analytical techniques to track volcanic seismicity and warn of hazardous eruptions. The network is currently undergoing an extensive upgrade, replacing aging short-period analog seismometers with digital broadband instruments. These are expected to improve AVO’s seismic capability and further facilitate other geophysical instrumentation such as continuous Global Positioning System receivers, infrasound sensors, and web cams.

Seismological Research Letters↗

Seismo-acoustic evidence for vent drying during shallow submarine eruptions at Bogoslof volcano, Alaska

Characterizing the state of the volcanic vent is key for interpreting observational datasets and accurately assessing volcanic hazards. This is particularly true for remote, complex eruptions such as the 2016–2017 Bogoslof volcano, Alaska eruption sequence. Bogoslof’s eruptions in this period were either shallow submarine or subaerial, or some combination of both. Our results demonstrate how low-frequency sound waves (infrasound), integrated with seismic and satellite data, can provide unique insight into shallow vent processes, otherwise not available. We use simple metrics, such as the infrasound frequency index (FI), event duration, and acoustic-seismic amplitude ratio, to look at changes in the elastic energy radiation and infer changes in seawater access to the vent. Satellite imagery before and after selected eruptions is used to ground-truth inferences on vent conditions. High FI and gradual increases in infrasound frequency content at Bogoslof correspond with transitions from submarine to subaerial vent conditions and a diminished or absent role of water, likely resulting in a drying out of the vent region. Event durations generally correlate with high FI and the range of FI values for each event, suggesting long duration events were more effective at drying out the vent region. A trend from low to high acoustic-seismic amplitude ratios for some long duration events also suggests an increase in acoustic efficiency as the vent dried out. We demonstrate that infrasound can serve as a robust indicator of seawater involvement for Bogoslof and other shallow submarine eruptions that may not be inferable from other datasets, particularly in near-real-time.

Alaska↗

Infrasound from giant bubbles during explosive submarine eruptions

Shallow submarine volcanoes pose unique scientific and monitoring challenges. The interaction between water and magma can create violent explosions just below the surface, but the inaccessibility of submerged volcanoes means they are typically not instrumented. This both increases the risk to marine and aviation traffic and leaves the underlying eruption physics poorly understood. Here we use low-frequency sound in the atmosphere (infrasound) to examine the source mechanics of shallow submarine explosions from Bogoslof volcano, Alaska. We show that the infrasound originates from the oscillation and rupture of magmatic gas bubbles that initially formed from submerged vents, but that grew and burst above sea level. We model the low-frequency signals as overpressurized gas bubbles that grow near the water–air interface, which require bubble radii of 50–220 m. Bubbles of this size and larger have been described in explosive subaqueous eruptions for more than a century, but we present a unique geophysical record of this phenomenon. We propose that the dominant role of seawater during the effusion of gas-rich magma into shallow water is to repeatedly produce a gas-tight seal near the vent. This resealing mechanism leads to sequences of violent explosions and the release of large, bubble-forming volumes of gas—activity we describe as hydrovulcanian.

Alaska↗

AVO-G2S: A modified, open-source Ground-to-Space atmospheric specification for infrasound modeling

To facilitate infrasound propagation studies, we present AVO-G2S, an open-source, Ground-to-Space model which provides temperature and wind specifications from the surface to an altitude of 225 km. This model provides a means of smoothly characterizing atmospheric conditions using multiple numerical weather prediction forecast and reanalysis products, along with upper-atmospheric empirical models. Regional atmospheric reconstructions only require a limited domain and can utilize high-resolution numerical weather prediction forecasts typically provided on a projected grid. The use of a projected grid allows for faster spectral transform libraries to be employed. The AVO-G2S software can also provide global reconstructions that rely on global numerical weather prediction products and spherical harmonic decompositions. AVO-G2S is inspired by a global Ground-to-Space model developed by the Naval Research Laboratory, and relies on their empirical descriptions of upper-atmospheric conditions. Alaska Volcano Observatory has implemented this model for near-real-time infrasound monitoring of volcanic eruptions and historical research projects.

Computers & Geosciences↗

Application of an updated atmospheric model to explore volcano infrasound propagation and detection in Alaska

Winds and temperature gradients greatly affect the long-range propagation of infrasound . The spatio-temporal variability of these parameters must therefore be accurately characterized to correctly interpret recorded infrasound at long distances, specifically to differentiate between source and propagation effects. Here we present the first results of an open source reanalysis model, termed Alaska Volcano Observatory Ground-to-Space (AVO-G2S), constructed to accurately characterize the atmosphere and model long-range infrasound propagation from volcanic eruptions in Alaska. We select a number of case studies to examine recent eruptions of Alaskan volcanoes whose ash emissions posed a threat to air traffic, including the two most recent eruptions of Pavlof Volcano and two typical explosions from Cleveland Volcano. Strong tropospheric ducting and low noise at the station during the 21 July 2015 explosion of Cleveland Volcano led to an automated detection of the explosion at an infrasound array 992 km away, whereas low signal-to-noise ratio for the 6 November 2014 Cleveland Volcano explosion helps explain the non-detection in real-time of a predicted strong stratospheric arrival. For the November 2014 Pavlof eruption, discrepancies between local seismic data and a distal infrasound array 460 km away cannot be solely explained by changes in atmospheric conditions, though some features of the complex propagation predictions follow the trends in long-range infrasound signals. The most recent eruption of Pavlof Volcano in March 2016 shows minimal changes in propagation conditions throughout the eruption and therefore indicates that the signals detected at long-range primarily reflect source processes. These results show how detailed examination of the acoustic propagation conditions provides insight into detection capability and eruption dynamics. Future work will implement AVO-G2S and high-resolution long-range infrasound propagation modeling in real-time for Alaskan volcanoes of interest.

Alaska↗

Short-term forecasting and detection of explosions during the 2016–2017 eruption of Bogoslof volcano, Alaska

We describe a multidisciplinary approach to forecast, rapidly detect, and characterize explosive events during the 2016–2017 eruption of Bogoslof volcano, a back-arc shallow submarine volcano in Alaska’s Aleutian arc. The eruptive sequence began in December 2016 and included about 70 discrete explosive events. Because the volcano has no local monitoring stations, we used distant stations on the nearest volcanoes, Okmok (54 km) and Makushin (72 km), combined with regional infrasound sensors and lightning detection from the Worldwide Lightning Location Network (WWLLN). Pre-eruptive seismicity was detected for 12 events during the first half of the eruption; for all other events co-eruptive signals allowed for detection only. Monitoring of activity used a combination of scheduled checks combined with automated alarms. Alarms triggered on real-time data included real-time seismic amplitude measurement (RSAM); infrasound from several arrays, the closest being on Okmok; and lightning strokes detected from WWLLN within a 20-km radius of the volcano. During periods of unrest, a multidisciplinary response team of four people fulfilled specific roles to evaluate geophysical and remote-sensing data, run event-specific ash-cloud dispersion models, ensure interagency coordination, and develop and distribute of formalized warning products. Using this approach, for events that produced ash clouds ≥7.5 km above sea level, Alaska Volcano Observatory (AVO) called emergency response partners 15 min, and issued written notices 30 min, after event onset (mean times). Factors that affect timeliness of written warnings include event size and number of data streams available; bigger events and more data both decrease uncertainty and allow for faster warnings. In remote areas where airborne ash is the primary hazard, the approach used at Bogoslof is an effective strategy for hazard mitigation.

Alaska↗

Infrasound signal detection and back azimuth estimation using ground-coupled airwaves on a seismo-acoustic sensor pair

We present a new infrasonic signal detection and back azimuth determination technique that requires just one microphone and one three-component seismometer. Ground-coupled airwaves (GCAs) occur when an incident atmospheric acoustic wave impinges on the ground surface and is partially transmitted as a seismic wave. GCAs are commonly detected hundreds of kilometers away on seismic networks and are observed to have retrograde particle motion. Horizontally propagating acoustic waves and GCAs have previously been observed on collocated infrasound and seismic sensor pairs as coherent with a 90° phase difference. If the sensors are spatially separated, an additional propagation-induced phase shift is present. The additional phase shift depends on the direction from which the acoustic wave arrives, as each back azimuth has a different apparent distance between the sensors. We use the additional phase shift, the coherence, and the characteristic particle motion on the three-component seismometer to determine GCA arrivals and their unique back azimuth. We test this technique with synthetic seismo-acoustic data generated by a coupled Earth-atmosphere 3-D finite difference code, as well as three seismo-acoustic data sets from Mount St. Helens, Mount Cleveland, and Mount Pagan volcanoes. Results from our technique compare favorably with traditional infrasound array processing and provide robust GCA detection and back azimuth determination. Assuming adequate station spacing and sampling, our technique provides a new and robust method to detect infrasonic signals and determine their back azimuth, and may be of practical benefit where resources are limited and large sensor networks or arrays are not feasible.

JGR Solid Earth↗

Local, regional, and remote seismo-acoustic observations of the April 2015 VEI 4 eruption of Calbuco volcano, Chile

The two major explosive phases of the 22-23 April 2015 eruption of Calbuco volcano, Chile produced powerful seismicity and infrasound. The eruption was recorded on seismo-acoustic stations out to 1,540 km and on 5 stations (IS02, IS08, IS09, IS27, and IS49) of the International Monitoring System (IMS) infrasound network at distances from 1,525 to 5,122 km. The remote IMS infrasound stations provide an accurate explosion chronology consistent with the regional and local seismo-acoustic data, and with previous studies of lightning and plume observations. We use the IMS network to detect and locate the eruption signals using a brute-force, grid search, cross-bearings approach. After incorporating azimuth deviation corrections from stratospheric cross-winds using 3D ray-tracing, the estimated source location is 172 km from true. This case study highlights the significant capability of the IMS infrasound network to provide automated detection, characterization, and timing estimates of global explosive volcanic activity. Augmenting the IMS with regional seismo-acoustic networks will dramatically enhance volcanic signal detection, reduce latency, and improve discrimination capability.

Calbuco volcano↗

Volcanic thunder from explosive eruptions at Bogoslof volcano, Alaska

Lightning often occurs during ash-producing eruptive activity, and its detection is now being used in volcano monitoring for rapid alerts. We report on infrasonic and sonic recordings of the related, but previously undocumented, phenomenon of volcanic thunder. We observe volcanic thunder during the waning stages of two explosive eruptions at Bogoslof volcano, Alaska, on a microphone array located 60 km away. Thunder signals arrive from a different direction than coeruptive infrasound generated at the vent following an eruption on 10 June 2017, consistent with locations from lightning networks. For the 8 March 2017 eruption, arrival times and amplitudes of high-frequency thunder signals correlate well with the timing and strength of lightning detections. In both cases, the thunder is associated with lightning that continues after significant eruptive activity has ended. Infrasonic and sonic observations of volcanic thunder offer a new avenue for studying electrification processes in volcanic plumes.

Alaska↗

Seismic equivalents of volcanic jet scaling laws and multipoles in acoustics

We establish analogies between equivalent source theory in seismology (moment-tensor and single-force sources) and acoustics (monopoles, dipoles and quadrupoles) in the context of volcanic eruption signals. Although infrasound (acoustic waves < 20 Hz) from volcanic eruptions may be more complex than a simple monopole, dipole or quadrupole assumption, these elementary acoustic sources are a logical place to begin exploring relations with seismic sources. By considering the radiated power of a harmonic force source at the surface of an elastic half-space, we show that a volcanic jet or plume modelled as a seismic force has similar scaling with respect to eruption parameters (e.g. exit velocity and vent area) as an acoustic dipole. We support this by demonstrating, from first principles, a fundamental relationship that ties together explosion, torque and force sources in seismology and highlights the underlying dipole nature of seismic forces. This forges a connection between the multipole expansion of equivalent sources in acoustics and the use of forces and moments as equivalent sources in seismology. We further show that volcanic infrasound monopole and quadrupole sources exhibit scalings similar to seismicity radiated by volume injection and moment sources, respectively. We describe a scaling theory for seismic tremor during volcanic eruptions that agrees with observations showing a linear relation between radiated power of tremor and eruption rate. Volcanic tremor over the first 17 hr of the 2016 eruption at Pavlof Volcano, Alaska, obeyed the linear relation. Subsequent tremor during the main phase of the eruption did not obey the linear relation and demonstrates that volcanic eruption tremor can exhibit other scalings even during the same eruption.

Geophysical Journal International↗

Volcanic tremor and plume height hysteresis from Pavlof Volcano, Alaska

The March 2016 eruption of Pavlof Volcano, Alaska produced ash plumes that canceled over 100 flights in North America. The eruption produced strong tremor recorded by seismic and remote low-frequency acoustic (infrasound) stations, including the EarthScope Transportable Array. The relationship between the tremor amplitudes and plume height changes considerably between the waxing and waning portions of the eruption. Similar hysteresis has been observed between seismic river noise and discharge during storms, suggesting flow and erosional processes in both rivers and volcanoes can produce irreversible structural changes detectable in geophysical data. We propose that the time-varying relationship at Pavlof arose from changes in the tremor source related to volcanic vent erosion. This relationship may improve estimates of volcanic emissions and characterization of eruption size and intensity.

Alaska↗

Seismic envelope-based detection and location of ground-coupled airwaves from volcanoes in Alaska

Volcanic explosions and other infrasonic sources frequently produce acoustic waves that are recorded by seismometers. Here we explore multiple techniques to detect, locate, and characterize ground‐coupled airwaves (GCA) on volcano seismic networks in Alaska. GCA waveforms are typically incoherent between stations, thus we use envelope‐based techniques in our analyses. For distant sources and planar waves, we use f ‐ k beamforming to estimate back azimuth and trace velocity parameters. For spherical waves originating within the network, we use two related time difference of arrival (TDOA) methods to detect and localize the source. We investigate a modified envelope function to enhance the signal‐to‐noise ratio and emphasize both high energies and energy contrasts within a spectrogram. We apply these methods to recent eruptions from Cleveland, Veniaminof, and Pavlof Volcanoes, Alaska. Array processing of GCA from Cleveland Volcano on 4 May 2013 produces robust detection and wave characterization. Our modified envelopes substantially improve the short‐term average/long‐term average ratios, enhancing explosion detection. We detect GCA within both the Veniaminof and Pavlof networks from the 2007 and 2013–2014 activity, indicating repeated volcanic explosions. Event clustering and forward modeling suggests that high‐resolution localization is possible for GCA on typical volcano seismic networks. These results indicate that GCA can be used to help detect, locate, characterize, and monitor volcanic eruptions, particularly in difficult‐to‐monitor regions. We have implemented these GCA detection algorithms into our operational volcano‐monitoring algorithms at the Alaska Volcano Observatory.

Alaska↗

The 2013 eruption of Pavlof Volcano, Alaska: a spatter eruption at an ice- and snow-clad volcano

The 2013 eruption of Pavlof Volcano, Alaska began on 13 May and ended 49 days later on 1 July. The eruption was characterized by persistent lava fountaining from a vent just north of the summit, intermittent strombolian explosions, and ash, gas, and aerosol plumes that reached as high as 8 km above sea level and on several occasions extended as much as 500 km downwind of the volcano. During the first several days of the eruption, accumulations of spatter near the vent periodically collapsed to form small pyroclastic avalanches that eroded and melted snow and ice to form lahars on the lower north flank of the volcano. Continued lava fountaining led to the production of agglutinate lava flows that extended to the base of the volcano, about 3&ndash;4 km beyond the vent. The generation of fountain-fed lava flows was a dominant process during the 2013 eruption; however, episodic collapse of spatter accumulations and formation of hot spatter-rich granular avalanches was a more efficient process for melting snow and ice and initiating lahars. The lahars and ash plumes generated during the eruption did not pose any serious hazards for the area. However, numerous local airline flights were cancelled or rerouted, and trace amounts of ash fall occurred at all of the local communities surrounding the volcano, including Cold Bay, Nelson Lagoon, Sand Point, and King Cove.

Alaska↗

Distinguishing high surf from volcanic long-period earthquakes

Repeating long-period (LP) earthquakes are observed at active volcanoes worldwide and are typically attributed to unsteady pressure fluctuations associated with fluid migration through the volcanic plumbing system. Nonvolcanic sources of LP signals include ice movement and glacial outburst floods, and the waveform characteristics and frequency content of these events often make them difficult to distinguish from volcanic LP events. We analyze seismic and infrasound data from an LP swarm recorded at Pagan volcano on 12–14 October 2013 and compare the results to ocean wave data from a nearby buoy. We demonstrate that although the events show strong similarity to volcanic LP signals, the events are not volcanic but due to intense surf generated by a passing typhoon. Seismo-acoustic methods allow for rapid distinction of volcanic LP signals from those generated by large surf and other sources, a critical task for volcano monitoring.

Pagan Volcano↗

Evaluation of Redoubt Volcano's sulfur dioxide emissions by the Ozone Monitoring Instrument

The 2009 eruption of Redoubt Volcano, Alaska, provided a rare opportunity to compare satellite measurements of sulfur dioxide (SO 2 ) by the Ozone Monitoring Instrument (OMI) with airborne SO 2 measurements by the Alaska Volcano Observatory (AVO). Herein we: (1) compare OMI and airborne SO 2 column density values for Redoubt's tropospheric plume, (2) calculate daily SO 2 masses from Mount Redoubt for the first three months of the eruption, (3) develop simple methods to convert daily measured SO 2 masses into emission rates to allow satellite data to be directly integrated with the airborne SO 2 emissions dataset, (4) calculate cumulative SO 2 emissions from the eruption, and (5) evaluate OMI as a monitoring tool for high-latitude degassing volcanoes. A linear correlation (R 2 ~ 0.75) is observed between OMI and airborne SO 2 column densities. OMI daily SO 2 masses for the sample period ranged from ~ 60.1 kt on 24 March to below detection limit, with an average daily SO 2 mass of ~ 6.7 kt. The highest SO 2 emissions were observed during the initial part of the explosive phase and the emissions exhibited an overall decreasing trend with time. OMI SO 2 emission rates were derived using three methods and compared to airborne measurements. This comparison yields a linear correlation (R 2 ~ 0.82) with OMI-derived emission rates consistently lower than airborne measurements. The comparison results suggest that OMI's detection limit for high latitude, springtime conditions varies from ~ 2000 to 4000 t/d. Cumulative SO 2 masses calculated from daily OMI data for the sample period are estimated to range from 542 to 615 kt, with approximately half of this SO 2 produced during the explosive phase of the eruption. These cumulative masses are similar in magnitude to those estimated for the 1989–90 Redoubt eruption. Strong correlations between daily OMI SO 2 mass and both tephra mass and acoustic energy during the explosive phase of the eruption suggest that OMI data may be used to infer relative eruption size and explosivity. Further, when used in conjunction with complementary datasets, OMI daily SO 2 masses may be used to help distinguish explosive from effusive activity and identify changes in lava extrusion rates. The results of this study suggest that OMI is a useful volcano monitoring tool to complement airborne measurements, capture explosive SO 2 emissions, and provide high temporal resolution SO 2 emissions data that can be used with interdisciplinary datasets to illuminate volcanic processes.

Alaska↗

Detecting hidden volcanic explosions from Mt. Cleveland Volcano, Alaska with infrasound and ground-couples airwaves

In Alaska, where many active volcanoes exist without ground-based instrumentation, the use of techniques suitable for distant monitoring is pivotal. In this study we report regional-scale seismic and infrasound observations of volcanic activity at Mt. Cleveland between December 2011 and August 2012. During this period, twenty explosions were detected by infrasound sensors as far away as 1827 km from the active vent, and ground-coupled acoustic waves were recorded at seismic stations across the Aleutian Arc. Several events resulting from the explosive disruption of small lava domes within the summit crater were confirmed by analysis of satellite remote sensing data. However, many explosions eluded initial, automated, analyses of satellite data due to poor weather conditions. Infrasound and seismic monitoring provided effective means for detecting these hidden events. We present results from the implementation of automatic infrasound and seismo-acoustic eruption detection algorithms, and review the challenges of real-time volcano monitoring operations in remote regions. We also model acoustic propagation in the Northern Pacific, showing how tropospheric ducting effects allow infrasound to travel long distances across the Aleutian Arc. The successful results of our investigation provide motivation for expanded efforts in infrasound monitoring across the Aleutians and contributes to our knowledge of the number and style of vulcanian eruptions at Mt. Cleveland.

Alaska↗