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.

Over the past two decades, the amount of available seismic data has increased significantly, fueling the need for automatic processing to use the vast amount of information contained in such data sets. Detecting seismicity in temporary aftershock networks is one important example that has become a huge challenge because of the high seismicity rate and dense station coverage. Additionally, the need for highly accurate earthquake locations to distinguish between different competing physical processes during the postseismic period demands even more accurate arrival‐time estimates of seismic phase. Here, we present a convolutional neural network (CNN) for classifying seismic phase onsets for local seismic networks. The CNN is trained on a small dataset for deep‐learning purposes (411 events) detected throughout northern Chile, typical for a temporary aftershock network. In the absence of extensive training data, we demonstrate that a CNN‐based automatic phase picker can still improve performance in classifying seismic phases, which matches or exceeds that of historic methods. The trained network is tested against an optimized short‐term average/long‐term average (STA/LTA) based method (Rietbrock et al., 2012) in classifying phase onsets for a separate dataset of 3878 events throughout the same region. Based on station travel‐time residuals, the CNN outperforms the STA/LTA approach and achieves location residual distribution close to the ones obtained by manual inspection.

Summary. The purpose of this work is to gain insights into the 2011–2012 eruption of El Hierro (Canary Islands) by mapping the evolution of the seismic b-value