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 is the surface-guided Lamb wave ( ≲ 0.01 Hz), which we observed propagating for four (+three antipodal) passages around the Earth over six days. Based on 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–20 Hz), long-range (~10,000 km) 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. We highlight exceptional observations of the atmospheric waves.

[1] The source mechanism of very long period (VLP) signals accompanying volcanic degassing bursts at Popocatépetl is analyzed in the 15-70 s band by minimizing the residual error between data and synthetics calculated for a point source embedded in a homogeneous medium. The waveforms of two eruptions (23 April and 23 May 2000) representative of mild Vulcanian activity are well reproduced by our inversion, which takes into account volcano topography. The source centroid is positioned 1500 m below the western perimeter of the summit crater, and the modeled source is composed of a shallow dipping crack (sill with easterly dip of 10°) intersecting a steeply dipping crack (northeast striking dike dipping 83°northwest), whose surface extension bisects the vent. Both cracks undergo a similar sequence of inflation, deflation, and reinflation, reflecting a cycle of pressurization, depressurization, and repressurization within a time interval of 3-5 min. The largest moment release occurs in the sill, showing a maximum volume change of 500-1000 m 3 , pressure drop of 3-5 MPa, and amplitude of recovered pressure equal to 1.2 times the amplitude of the pressure drop. In contrast, the maximum volume change in the dike is less (200-300 m 3 ), with a corresponding pressure drop of 1-2 MPa and pressure recovery equal to the pressure drop. Accompanying these volumetric sources are single-force components with magnitudes of 10 8 N, consistent with melt advection in response to pressure transients. The source time histories of the volumetric components of the source indicate that significant mass movement starts within the sill and triggers a mass movement response in the dike within a few seconds. Such source behavior is consistent with the opening of a pathway for escape of pent-up gases from slow pressurization of the sill driven by magma crystallization. The opening of this pathway and associated rapid evacuation of volcanic gases induces the pressure drop. Pressure recovery in the magma filling the sill is driven by diffusion of gases from the resulting supersaturated melt into bubbles. Assuming a penny-shaped crack at ambient pressure of 40 MPa, the observed pressure and volume variations can be modeled with the following attributes: crack radius (100 m), crack aperture (5 m), bubble number density (10 10 -10 12 m À3 ), initial bubble radius (10 À6 m), final bubble radius ($10 À5 m), and net decrease of gas concentration in the melt (0.01 wt %).Citation: Chouet, B., P. Dawson, and A. Arciniega-Ceballos (2005), Source mechanism of Vulcanian degassing at Popocatépetl Volcano, Mexico, determined from waveform inversions of very long period signals,

The Oaxaca subduction zone is an ideal area for detailed studies of plate boundary deformation as rapid convergent rates, shallow subduction, and short trench‐to‐coast distances bring the thermally defined seismogenic and transition zones of the plate interface over 100 km inland. Previous analysis of slow slip events in southern Mexico suggests that they may represent motion in the transition zone, defining the downdip edge of future megathrust earthquakes. A new deployment consisting of broadband seismometers distributed inland along the Oaxaca segment provide the means to examine whether nonvolcanic tremor (NVT) signals can also be used to characterize the boundary between the seismogenic and transition zones. In this study, we established that NVT exists in the Oaxaca region based on waxing and waning of seismic energy on filtered day‐long seismograms that were correlated across neighboring stations and were further supported by appropriate relative time moveouts in record sections and spectrograms with narrow frequency bands. Eighteen prominent NVT episodes that lasted upwards of a week were identified during the 15 months analyzed (June 2006 to September 2007), recurring as frequently as every 2–3 months in a given region. We analyze NVT envelope waveforms with a semiautomated process for identifying prominent energy bursts, and analyst‐refined relative arrival times are inverted for source locations. NVT burst epicenters primarily occur between the 40–50 km contours for depth of the plate interface, except in eastern Oaxaca where they shift toward the 30 km contour as the slab steepens. NVT hypocenters correlate well with a high conductivity zone that is interpreted to be due to slab fluids. NVT is more frequent, shorter in duration, and located further inland than GPS‐detected slow slip, while the latter is associated with a zone of ultra‐slow velocity interpreted to represent high pore fluid pressure. This zone of slow slip corresponds to approximately 350°C–450°C, with megathrust earthquakes, microseismicity, and strong long‐term coupling occurring immediately updip from it. This leaves NVT primarily in a region further inland from the thermally defined transition zone, suggesting that transition from locking to free slip may occur in more than one phase.

Discrete VLP signals accompanying long-period (LP) events also share similar signatures and have dominant periods that are nearly identical to those observed in the VLP waveforms of explosions. The VLP particle motions for eruption onsets consistently point to a source located a few km beneath the crater. The VLP ground displacement response to each explosion is marked by a compression, followed by a dilatation and terminating with another compression, suggesting a sequence of pressurization-depressurization-repressurization of the conduit. The repetitive nature of the waveforms points to a non-destructive source process which has remained active in the magmatic system of Popocatepetl at least since April 1997.

We examine the along-strike transition from flat to steeper subduction in Oaxaca, Mexico, to provide a better understanding of what controls the slab morphology. Prior studies have suggested the slab tends to tear along the transitions in dip as the slab rolls back. We determine the slab geometry based on local seismicity, nonvolcanic tremor (NVT), and slow slip utilizing a deployment of broadband seismometers and continuous GPS receivers distributed in and around Oaxaca. We construct depth contours of the subducting slab surface down to 100 km, which illustrate that the transition from flat to steeper subduction occurs rapidly via a sharper flexure than previously recognized. The prior catalog of NVT in Oaxaca is extended using the same method and additional stations that extend further west. The band of NVT follows the new slab contours, widening toward the west with the downdip extent gradually moving inland. The amount of NVT also correlates with the strength of an ultraslow-velocity layer. There are no gaps in seismicity, NVT, or slow slip across the rapid transition in slab dip, further supporting the notion that the slab is not currently torn in the updip region. We propose that the sharp flexure is possible in this region due to bending moment saturation that leads to greater curvature in both the downdip and along-strike directions. A similar set of observations in southern Peru suggests this is a viable alternative to tearing that accommodates the large strains from variable rates of slab rollback.FASOLA ET AL.

[1] The seismicity of Popocatépetl is dominated by longperiod and very-long period signals associated with hydrothermal processes and magmatic degassing. We model the source mechanism of repetitive long-period signals in the 0.4-2 s band from a 15-station broadband network by stacking long-period events with similar waveforms to improve the signal-to-noise ratio. The data are well fitted by a point source located within the summit crater $250 m below the crater floor and $200 m from the inferred magma conduit. The inferred source includes a volumetric component that can be modeled as resonance of a horizontal steam-filled crack and a vertical single force component. The long-period events are thought to be related to the interaction between the magmatic system and a perched hydrothermal system. Repetitive injection of fluid into the horizontal fracture and subsequent sudden discharge when a critical pressure threshold is met provides a non-destructive source process.

Following an initial phreatic eruption on 21 December 1994, activity at Popocatepetl has been dominated by fumarolic emissions interspersed with more energetic emissions of ashes and gases. A phase of repetitive dome-building and dome-destroying episodes began in March 1996 and is still ongoing at present. We describe the long-period (LP) seismicity accompanying eruptive activity at Popocatepetl from December 1994 through May 2000, using data from a three-component broadband seismometer located 5 km from the summit crater. The broadband records display a variety of signals, with periods ranging in the band 0.04-90 s. Long-period events and tremor with typical dominant periods in the range 0.3-2.0 s are the most characteristic signals observed at Popocatepetl. These signals appear to reflect volumetric sources driven by pressure fluctuations associated with the unsteady transport of gases beneath the crater. Very-long-period (VLP) signals are also observed in association with LP events and tremor. The VLP signals which accompany LP events display Ricker-like wavelets with periods near 36 s, whereas VLP signals associated with tremor waveforms typically show sustained oscillations at periods ranging up to 90 s. The spectra and particle motion patterns remain similar from event to event for the majority of LP and tremor signals analyzed during the time span of this study, suggesting a repeated, non-destructive activation of a common source. Hypocenters determined by phase pick analyses of selected LP events recorded by the seven-station, permanent Popocatepetl short-period network suggest that the majority of these events are confined to a source region in the top 1.5 km below the crater floor. The repetitive occurrences of VLP signals with closely matched waveform characteristics are consistent with a non-destructive reactivation of at least two sources. One source appears to coincide with the main source region of LP seismicity, whereas the second is a deeper source whose activity appears to be intimately linked with episodes of monochromatic tremor.

We simulate gas-burst and volcanic explosions under controlled laboratory conditions inducing fragmentation of volcanic rocks by rapid depressurization. A series of experiments were performed in a shocktube apparatus at room temperature and a pressure range of 4 to 20 MPa using Argon (Ar) gas and, particles or pumice samples of different porosities. The instrumentation of this system with high-precision piezoelectric sensors enabled us to capture elastic waves and to recognize their characteristic signatures. By relating these signals to physical processes in the wave field, we have been able to characterize the conduit mechanism and the source dynamics. We compare and discuss conspicuous features of the waveforms and frequency spectra of these experimental signals with those of volcanic origin. Despite the fact that these signals are different in amplitude (resulting from different scale conditions); our observations indicate that the physical processes that occur during simulated explosions and those that occur during volcanic eruptions yield comparable signatures in their respective records. The effects of the source-receiver configuration and resonance also have significant implications. All this suggests that the physical processes (e.g. pressurization and depressurization of a system) involve a system response that causes similar distinctive effects independent of the system size, reflecting its intrinsic dynamics. These similarities imply that powerful constraints on the source mechanisms of volcanic seismicity can emerge from seismic investigations of experimental simulations of volumetric sources. Such constraints may yield significant advances in the understanding of volcanic conduit dynamics and in the interpretation of seismic unrest at volcanic centers.