Search USGS⌕ Search

Geology topics

Kathleen McKee

Publications and source records attributed to Kathleen McKee.

9 recordsLinked to original sources

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↗

Remote characterization of the 12 January 2020 eruption of Taal Volcano, Philippines, using seismo-acoustic, volcanic lightning, and satellite observations

On 12 January 2020, an eruption began on the shores of the Main Crater Lake (MCL) of Taal Volcano—a caldera system on the southern end of Luzon Island in the Philippines. Taal, one of the most active volcanoes in the Philippines, is located 30 km south of Manila—a major metropolitan area with a population of 13.5 million people. Eruptive activity intensified throughout the day on 12 January, producing prolific volcanic lightning, ashfall, and a sustained plume that reached 16–17 km altitude. The chronology of events was well documented by the Philippine Institute of Volcanology and Seismology and the Tokyo Volcanic Ash Advisory Center. The wealth of data collected during the eruption provides a unique opportunity to investigate how the combination of different remote sensing methods may complement local observations and monitoring. Remote systems tend to provide lower resolution data but are also less likely to be compromised by the eruptive activity, thus providing continuous records of eruptive processes. Here, we present a postevent analysis of the 12 January activity, including data from long‐range lightning, infrasound, and seismic arrays located at distances up to several thousands of kilometers from the volcano. By combining these datasets, we distinguish five phases of activity and infer a major shift in eruption behavior around 12:00 on 12 January (UTC). The remote observations suggest that the most of the water within the MCL ( ⁠∼42 million m 3⁠ ) was vaporized and incorporated into the volcanic plume within the first 12 hr of the eruption.

Taal Volcano↗

Spatial analysis of globally detected volcanic lightning from the June 2019 eruption of Raikoke volcano, Kuril Islands

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.

Kuril Islands, Raikoke volcano↗

Atmospheric waves and global seismoacoustic observations of the January 2022 Hunga eruption, Tonga

The 15 January 2022 climactic eruption of Hunga volcano, Tonga, produced an explosion in the atmosphere of a size that has not been documented in the modern geophysical record. The event generated a broad range of atmospheric waves observed globally by various ground-based and spaceborne instrumentation networks. Most prominent was the surface-guided Lamb wave (≲0.01 hertz), which we observed propagating for four (plus three antipodal) passages around Earth over 6 days. As measured by the Lamb wave amplitudes, the climactic Hunga explosion was comparable in size to that of the 1883 Krakatau eruption. The Hunga eruption produced remarkable globally detected infrasound (0.01 to 20 hertz), long-range (~10,000 kilometers) audible sound, and ionospheric perturbations. Seismometers worldwide recorded pure seismic and air-to-ground coupled waves. Air-to-sea coupling likely contributed to fast-arriving tsunamis. Here, we highlight exceptional observations of the atmospheric waves.

Hunga volcano↗

Evaluating the state-of-the-art in remote volcanic eruption characterization Part I: Raikoke volcano, Kuril Islands

Raikoke, a small, unmonitored volcano in the Kuril Islands, erupted in June 2019. We integrate data from satellites (including Sentinel-2, TROPOMI, MODIS, Himawari-8), the International Monitoring System (IMS) infrasound network, and global lightning detection network (GLD360) with information from local authorities and social media to retrospectively characterize the eruptive sequence and improve understanding of the pre-, syn- and post- eruptive behavior. We observe six infrasound pulses beginning on 21 June at 17:49:55 UTC as well as the main Plinian phase on 21 June at 22:29 UTC. Each pulse is tracked in space and time using lightning and satellite imagery as the plumes drift eastward. Post-eruption visible satellite imagery shows expansion of the island's surface area, an increase in crater size, and a possibly-linked algal bloom south of the island. We use thermal satellite imagery and plume modeling to estimate plume height at 10–12 km asl and 1.5–2 × 106 kg/s mass eruption rate. Remote infrasound data provide insight into syn-eruptive changes in eruption intensity. Our analysis illustrates the value of interdisciplinary analyses of remote data to illuminate eruptive processes. However, our inability to identify deformation, pre-eruptive outgassing, and thermal signals, which may reflect the relatively short duration (~12 h) of the eruption and minimal land area around the volcano and/or the character of closed-system eruptions, highlights current limitations in the application of remote sensing for eruption detection and characterization.

Hokkaido, Sakhalin Oblast↗

Evaluating the state-of-the-art in remote volcanic eruption characterization Part II: Ulawun volcano, Papua New Guinea

Retrospective eruption characterization is valuable for advancing our understanding of volcanic systems and evaluating our observational capabilities, especially with remote technologies (defined here as a space-borne system or non-local, ground-based instrumentation which include regional and remote infrasound sensors). In June 2019, the open-system Ulawun volcano, Papua New Guinea, produced a VEI 4 eruption. We combined data from satellites (including Sentinel-2, TROPOMI, MODIS, Himawari-8), the International Monitoring System infrasound network, and GLD360 globally detected lightning with information from the local authorities and social media to characterize the pre-, syn - and post-eruptive behaviour. The Rabaul Volcano Observatory recorded ~24 h of seismicity and detected SO 2 emissions ~16 h before the visually-documented start of the Plinian phase on 26 June at 04:20 UTC. Infrasound and SO 2 detections suggest the eruption started during the night on 24 June 2019 at 10:39 UTC ~38 h before ash detections with a gas-dominated jetting phase. Local reports and infrasound detections show that the second phase of the eruption started on 25 June 19:28 UTC with ~6 h of jetting. The first detected lightning occurred on 26 June 00:14 UTC, and ash emissions were first detected by Himawari-8 at 01:00 UTC. Post-eruptive satellite imagery indicates new flow deposits to the south and north of the edifice and ash fall to the west and southwest. In particular, regional infrasound data provided novel insight into eruption onset and syn -eruptive changes in intensity. We conclude that, while remote observations are sufficient for detection and tracking of syn-eruptive changes, key challenges in data latency, acquisition, and synthesis must be addressed to improve future near-real-time characterization of eruptions at minimally-monitored or unmonitored volcanoes.

Bismark volcanic arc, New Britain Island, Ulawun v↗

Co-eruptive tremor from Bogoslof volcano: Seismic wavefield composition at regional distances

We analyze seismic tremor recorded during eruptive activity over the course of the 2016–2017 eruption of Bogoslof volcano, Alaska. Only regional recordings of the tremor wavefield exist for Bogoslof, making it a challenge to place the recordings in context with other eruptions that are normally captured by local seismic data. We apply a technique of time-frequency polarization analysis to three-component seismic data to reveal the wavefield composition of Bogoslof eruption tremor. We find that at regional distances, the tremor is dominated by P-waves in the band from 1.5 to 10 Hz. Using this information, along with an enriched Bogoslof earthquake catalog, we obtain estimates of average reduced displacement ( D R ) for eruption tremor during 25 of the 70 Bogoslof events. D R reaches as high as approximately 40 cm 2 for two of the major events, similar to other VEI~3 eruptions in Alaska. Overall, average reduced displacement displays a weak correlation to plume height during the first half of the 9-month-long eruption sequence, with a few notable exceptions. The two events with the highest D R values also generated measurable eruption tremor at very-long-periods (VLP) between 0.05 and 0.15 Hz.

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↗