The Tonga volcano eruption at 04:14:45 UT on 2022-01-15 released enormous amounts of energy into the atmosphere, triggering very significant geophysical variations not only in the immediate proximity of the epicenter but also globally across the whole atmosphere. This study provides a global picture of ionospheric disturbances over an extended period for at least 4 days. We find traveling ionospheric disturbances (TIDs) radially outbound and inbound along entire Great-Circle loci at primary speeds of ∼300–350 m/s (depending on the propagation direction) and 500–1,000 km horizontal wavelength for front shocks, going around the globe for three times, passing six times over the continental US in 100 h since the eruption. TIDs following the shock fronts developed for ∼8 h with 10–30 min predominant periods in near- and far- fields. TID global propagation is consistent with the effect of Lamb waves which travel at the speed of sound. Although these oscillations are often confined to the troposphere, Lamb wave energy is known to leak into the thermosphere through channels such as atmospheric resonance at acoustic and gravity wave frequencies, carrying substantial wave amplitudes at high altitudes. Prevailing Lamb waves have been reported in the literature as atmospheric responses to the gigantic Krakatoa eruption in 1883 and other geohazards. This study provides substantial first evidence of their long-duration imprints up in the global ionosphere. This study was enabled by ionospheric measurements from 5,000+ world-wide Global Navigation Satellite System (GNSS) ground receivers, demonstrating the broad implication of the ionosphere measurement as a sensitive detector for atmospheric waves and geophysical disturbances.

This paper presents observations of postsunset super plasma bubbles over China and adjacent areas during the second main phase of a storm on 8 September 2017. The signatures of the plasma bubbles can be seen or deduced from (1) deep field‐aligned total electron content depletions embedded in regional ionospheric maps derived from dense Global Navigation Satellite System networks, (2) significant equatorial and midlatitudinal plasma bite‐outs in electron density measurements on board Swarm satellites, and (3) enhancements of ionosonde virtual height and scintillation in local evening associated with strong southward interplanetary magnetic field. The bubbles/depletions covered a broad area mainly within 20°–45°N and 80°–110°E with bifurcated structures and persisted for nearly 5 hr (∼13–18 UT). One prominent feature is that the bubbles extended remarkably along the magnetic field lines in the form of depleted flux tubes, reaching up to midlatitude of around 50°N (magnetic latitude: 45.5°N) that maps to an altitude of 6,600 km over the magnetic equator. The maximum upward drift speed of the bubbles over the magnetic equator was about 700 m/s and gradually decreased with altitude and time. The possible triggering mechanism of the plasma bubbles was estimated to be storm time eastward prompt penetration electric field, while the traveling ionospheric disturbance could play a role in facilitating the latitudinal extension of the depletions.

Postsunset midlatitude traveling ionospheric disturbances (TIDs) and equatorial plasma bubbles (EPBs) were simultaneously observed over American sector during the geomagnetic storm on 8 September 2017. The characteristics of TIDs are analyzed by using a combination of the Millstone Hill incoherent scatter radar data and 2‐D detrended total electron content (TEC) from ground‐based Global Navigation Satellite System receivers. The main results associated with EPBs are as follows: (1) stream‐like structures of TEC depletion occurred simultaneously at geomagnetically conjugate points, (2) poleward extension of the TEC irregularities/depletions along the magnetic field lines, (3) severe equatorial and midlatitude electron density (Ne) bite outs observed by Defense Meteorological Satellite Program and Swarm satellites, and (4) enhancements of ionosphere F layer virtual height and vertical drifts observed by equatorial ionosondes near the EPBs initiation region. The stream‐like TEC depletions reached 46° magnetic latitudes that map to an apex altitude of 6,800 km over the magnetic equator using International Geomagnetic Reference Field. The formation of this extended density depletion structure is suggested to be due to the merging between the altitudinal/latitudinal extension of EPBs driven by strong prompt penetration electric field and midlatitude TIDs. Moreover, the poleward portion of the depletion/irregularity drifted westward and reached the equatorward boundary of the ionospheric main trough. This westward drift occurred at the same time as the sudden expansion of the convection pattern and could be attributed to the strong returning westward flow near the subauroral polarization stream region. Other possible mechanisms for the westward tilt are also discussed.

Following the 2022 Tonga Volcano eruption, dramatic suppression and deformation of the equatorial ionization anomaly (EIA) crests occurred in the American sector ∼14,000 km away from the epicenter. The EIA crests variations and associated ionosphere‐thermosphere disturbances were investigated using Global Navigation Satellite System total electron content data, Global‐scale Observations of the Limb and Disk ultraviolet images, Ionospheric Connection Explorer wind data, and ionosonde observations. The main results are as follows: (a) Following the eastward passage of expected eruption‐induced atmospheric disturbances, daytime EIA crests, especially the southern one, showed severe suppression of more than 10 TEC Unit and collapsed equatorward over 10° latitudes, forming a single band of enhanced density near the geomagnetic equator around 14–17 UT, (b) Evening EIA crests experienced a drastic deformation around 22 UT, forming a unique X‐pattern in a limited longitudinal area between 20 and 40°W. (c) Thermospheric horizontal winds, especially the zonal winds, showed long‐lasting quasi‐periodic fluctuations between ±200 m/s for 7–8 hr after the passage of volcano‐induced Lamb waves. The EIA suppression and X‐pattern merging was consistent with a westward equatorial zonal dynamo electric field induced by the strong zonal wind oscillation with a westward reversal.

This paper investigates the local and global ionospheric responses to the 2022 Tonga volcano eruption, using ground‐based Global Navigation Satellite System total electron content (TEC), Swarm in situ plasma density measurements, the Ionospheric Connection Explorer (ICON) Ion Velocity Meter (IVM) data, and ionosonde measurements. The main results are as follows: (a) A significant local ionospheric hole of more than 10 TECU depletion was observed near the epicenter ∼45 min after the eruption, comprising of several cascading TEC decreases and quasi‐periodic oscillations. Such a deep local plasma hole was also observed by space‐borne in situ measurements, with an estimated horizontal radius of 10–15° and persisted for more than 10 hr in ICON‐IVM ion density profiles until local sunrise. (b) Pronounced post‐volcanic evening equatorial plasma bubbles (EPBs) were continuously observed across the wide Asia‐Oceania area after the arrival of volcano‐induced waves; these caused a N_e decrease of 2–3 orders of magnitude at Swarm/ICON altitude between 450 and 575 km, covered wide longitudinal ranges of more than 140°, and lasted around 12 hr. (c) Various acoustic‐gravity wave modes due to volcano eruption were observed by accurate Beidou geostationary orbit (GEO) TEC, and the huge ionospheric hole was mainly caused by intense shock‐acoustic impulses. TEC rate of change index revealed globally propagating ionospheric disturbances at a prevailing Lamb‐wave mode of ∼315 m/s; the large‐scale EPBs could be seeded by acoustic‐gravity resonance and coupling to less‐damped Lamb waves, under a favorable condition of volcano‐induced enhancement of dusktime plasma upward E×B drift and postsunset rise of the equatorial ionospheric F‐layer.

This study investigates midlatitude ionospheric variations during the super geomagnetic storm on 10–11 May 2024, utilizing multi‐instrument data from ground‐based sources (Global Navigation Satellite Systems receivers and a Fabry–Perot Interferometer) and space‐based measurements (Swarm and DMSP). We observed several distinct density gradient structures in the midlatitude ionosphere, with the main findings summarized as follows: (a) Significant zonal plasma density enhancements developed continuously in local dusk across the American‐Pacific‐Asian longitude sectors around geomagnetic latitude. These midlatitude peaks exhibited a wide longitudinal extension exceeding 150 and a prolonged duration of 12–15 hr during the late main phase and early recovery phase of the storm. (b) Strong storm‐enhanced density (SED) was observed in both hemispheres yet with different longitudinal and universal time preferences. In the Northern Hemisphere, significant SED occurred over the American longitude sector during 20:30–22:30 UT on May 10. In the Southern Hemisphere, pronounced SED was observed not only in the American longitudes during 20:30–22:30 UT on May 10 but also in the Australian longitude sector during 02:00–04:00 UT on May 11.