USGS Science⌕ Search

USGS · 70282855

Alaska upper crustal velocities revealed by air-to-ground coupled waves from the 2022 Hunga Tonga-Hunga Ha’apai eruption

Abstract

Pressure changes in the atmosphere couple to the solid earth, producing ground motions that contain information about local crustal elastic parameters. This type of air-to-ground coupled wave was observed globally following the largest explosion of the instrumental age, the climactic eruption of the Hunga Tonga-Hunga Ha’apai volcano on 15 th January, 2022. We utilize this unprecedented source, along with the presence of colocated seismometers, infrasound sensors, and barometers in Alaska, to examine coupling and reveal elastic parameters beneath the stations. We derive coupling spectra by forming seismic--to--pressure amplitude ratios as a function of frequency, and identify passbands of high coherence between the pressure and seismic records. By relating coupling spectra in high-coherence bands to elastic parameters, we estimate mean shear wave velocities under stations to a depth encompassing much of the upper crust. Our velocity estimates from low-frequency coupling exhibit good agreement with a previously existing tomographic velocity model from Berg et al. (2020), while estimates from high-frequency coupling show considerable scatter when compared to proxy V s30 , even though the overall values are reasonable. In addition to providing velocity estimates, our results also indicate that, for the broadband pressure signals from the Hunga Tonga-Hunga Ha’apai eruption, microseismic noise exerts a strong effect on the frequency bands where coupling is observed, and that the air-to-ground coupled waves exhibit significant complexity not necessarily described by theory. Our results show that coupling observations provide a simple forward observation of mean seismic velocities beneath seismoacoustic stations, without the need to resort to complex inversion schemes. It is remarkable that pressure waves generated thousands of kilometers away are able to reveal the seismic velocity structure of Alaska to several kilometers depth.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 50° to 72° latitude; -175° to -130° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Kenneth Macpherson, David Fee, Stefan Awender, Bryant Chow, Juliann Colwell, Sam Delamere, Matthew M. Haney. 2025-07-08. Alaska upper crustal velocities revealed by air-to-ground coupled waves from the 2022 Hunga Tonga-Hunga Ha’apai eruption. https://doi.org/10.26443/seismica.v4i2.1557

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

CRESCENT earthquake dynamic rupture, earthquake cycle, and tsunami code verification platform

Physics-based simulations are critical for understanding natural hazards. The increasing complexity of numerical codes requires benchmark exercises to verify that different computational methods yield consistent results when solving the same governing equations. Here, we present an open-access web platform designed for the verification of earthquake dynamic rupture, seismic cycle, and tsunami simulations. The platform architecture utilizes a modular, serverless backend on Amazon Web Services (AWS) to provide scalable file processing and visualization. A lightweight static web application provides a secure interface for uploading and managing results, while the browser-based data visualization enables interactive analysis of time series and surface grid data. By using structured JavaScript Object Notation (JSON) text files to define benchmark structures, the system remains fully extensible, allowing the addition of new scenarios without modifying the underlying software logic. The platform hosts the "The Tsunami Problem Versions" (TTPV) 1 & 2, two benchmarks for 3D fully coupled earthquake dynamic rupture and tsunami generation, and provides a framework for earthquake cycle models. This community resource aims to build trust in numerical simulations and facilitate long-term collaborative code verification as modeling software continues to evolve.

Seismica↗

Small earthquake moment magnitude and implications for frequency–magnitude scaling of injection induced earthquakes of the Raton Basin

Accurate estimation of earthquake source parameters—such as moment magnitudes, corner frequencies, and stress drops—is essential for improving seismic hazard assessments and understanding earthquake physics. In this study, moment magnitudes ( M W ) are calculated for 31,581 earthquakes associated with wastewater injection in the Raton Basin (located along the border between northern New Mexico and southern Colorado) between 2016 and 2024 using radiative transfer theory to fit coda decay envelopes. Our results show that it is feasible to estimate moment magnitudes down to M W ~1 with coda envelopes from a small local monitoring network. Significant differences were found between M W and local magnitudes ( M L ) for small earthquakes ( M < 3.0). A linear relationship was optimized to convert M L to M W : M W = 0.7 M L + 0.96 and M W = 0.73 M L + 0.99 (for the events reported by the U.S. Geological Survey), which can be applied in future studies of Raton Basin seismicity. We find that b -values calculated employing different methods and using M L are approximately 1.0, while those using M W range from 1.2 to 1.4. A larger estimate of the b -value could influence interpretations of the statistical behavior of earthquakes associated with injection and consequently seismic hazard assessments based on a magnitude–frequency distribution. The potential differences between local versus moment magnitude-based earthquake statistics should be considered in other seismically active regions.

Colorado, New Mexico↗

3D Converted wave reverse time migration imaging

We describe a newly developed method for recovering high-resolution images of seismic discontinuities, such as subducting slabs, in 3D. Our method makes use of converted S → P or P → S waves observed by dense arrays of seismometers to infer the locations and relative strengths of seismic discontinuities at depth in a target region. Observed direct and converted waves are backpropagated to their times of origin. The time-reversed wavefield is then separated into its constituent P and S components via the Helmholtz decomposition, and those separated wavefields are used to compute imaging functions that characterize the locations and relative strengths of seismic discontinuities. Imaging functions may be designed to use either S → P or P → S waves, so that users can target those arrivals expected to be most dominant in a given dataset. We have previously demonstrated the efficacy of our method in two dimensions, and we now present a 3D implementation of our technique which addresses the significant computational challenges posed by the size of volumetric wavefield data in three dimensions. Through a series of synthetic examples, we demonstrate that our method is capable of recovering the fine scale structure of a subducting slab given realistic station coverage and earthquake sources. We investigate optimal seismic station geometries for our technique and explore image interpretability in regions with poor data coverage. We find that linear station geometries yield more optimal, interpretable imaging functions than collections of small arrays can. We also show that our method can successfully recover bothS → P or P → S images when realistic shear earthquake sources are used, and we explore the additional computational challenges presented by the high frequency content of S waves. Our results demonstrate the potential for our technique to recover high-resolution information about subducting slabs in real-world regions, given that relatively sparse seismic arrays with only approximately 100 stations are capable of recovering interpretable imaging functions from just a few realistic earthquake sources for multiple discontinuities at significant depth in an area of approximately 400~sq~km.

Seismica↗