The May 2024 super storm is one of the strongest geomagnetic storms during the past 20 years. One of the most remarkable ionospheric responses to this event over East and Southeast Asia is the complex ionospheric fluctuations following the storm commencement. The fluctuations created multiple oscillations of total electron content (TEC) embedded in the diurnal variation, with amplitudes up to 10 TECu. Along the same latitude, the fluctuations were nearly synchronized over a wide longitude span up to 35°. In the meridional direction, the fluctuations over low latitudes were the most significant and complex, which contained two main components, the poleward extending oscillations originated from the magnetic equator, and the equatorward propagating traveling ionospheric disturbances (TIDs) from high latitudes. The TIDs likely occurred around the globe. The storm‐time interplanetary electric fields penetrating into equatorial latitudes of the ionosphere and the auroral energy input were suggested to drive the poleward extending oscillations and the equatorward TIDs, respectively.
Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher.
Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?
Some links on this page may take you to non-federal websites. Their policies may differ from this site.
-
Abstract Free, publicly-accessible full text available December 1, 2025 -
Abstract Plasma blob is generally a low‐latitude phenomenon occurring at the poleward edge of equatorial plasma bubble (EPB) during post‐sunset periods. Here we report a case of midlatitude ionospheric plasma blob‐like structures occurring along with super EPBs over East Asia around sunrise during the May 2024 great geomagnetic storm. Interestingly, the blob‐like structures appeared at both the poleward and westward edges of EPBs, reached up to 40°N magnetic latitudes, and migrated westward several thousand kilometers together with the bubble. The total electron content (TEC) inside the blob‐like structures was enhanced by ∼50 TEC units relative to the ambient ionosphere. The blob‐like structure at the EPB poleward edge could be partly linked with field‐aligned plasma accumulation due to poleward development of bubble. For the blob‐like structure at the EPB west side, one possible mechanism is that it was formed and enhanced accompanying the bubble evolution and westward drift.
Free, publicly-accessible full text available November 16, 2025 -
Abstract The global 3‐dimensional structure of the concentric traveling ionospheric disturbances (CTIDs) triggered by 2022 Tonga volcano was reconstructed by using the 3‐dimensional computerized ionospheric tomography (3DCIT) technique and extensive global navigation satellite system (GNSS) observations. This study provides the first estimation of the CTIDs vertical wavelengths, ∼736 km, which was much larger than the gravity wave (GW) vertical wavelength, 240–400 km, estimated using ICON neutral wind observations. Notable trend with the variation of azimuth was also found in horizontal speeds at 200 and 500 km altitudes and differences between them. These results imply that (a) the global propagation of Lamb waves determined the arrival time of local ionospheric disturbances, and (b) the arriving Lamb waves caused vertical atmospheric perturbations that are not typical of GWs, resulting in local thermospheric horizontal wave propagation which is faster than the Lamb wave propagation at lower altitudes.
-
Abstract This paper conducts a multi‐instrument and data assimilation analysis of the three‐dimensional ionospheric electron density responses to the total solar eclipse on 08 April 2024. The altitude‐resolved electron density variations over the continental US and adjacent regions are analyzed using the Millstone Hill incoherent scatter radar data, ionosonde observations, Swarm in situ measurements, and a novel TEC‐based ionospheric data assimilation system (TIDAS) with SAMI3 model as the background. The principal findings are summarized as follows: (a) The ionospheric hmF2 exhibited a slight enhancement in the initial phase of the eclipse, followed by a distinct reduction of 20–30 km in the recovery phase of the eclipse. The hmF2 in the umbra region showed a post‐eclipse fluctuation, characterized by wavelike perturbations of 10–25 km in magnitude and a period of 30 min. (b) There was a substantial reduction in ionospheric electron density of 20%–50% during the eclipse, with the maximum depletion observed in the F‐region around 200–250 km. The ionospheric electron density variation exhibited a significant altitude‐dependent feature, wherein the response time gradually delayed with increasing altitude. (c) The bottomside ionospheric electron density displayed an immediate reduction after local eclipse began, reaching maximum depletion 5–10 min after the maximum obscuration. In contrast, the topside ionospheric electron density showed a significantly delayed response, with maximum depletion occurring 1–2.5 hr after the peak obscuration.
Free, publicly-accessible full text available September 1, 2025 -
Abstract The storm‐enhanced density (SED) is a large‐scale midlatitude ionospheric electron density enhancement in the local afternoon sector, which exhibits substantial spatial gradients and thus can impose detrimental effects on modern navigation and communication systems, causing potential space weather hazards. This study has identified a comprehensive list of 49 SED events over the continental US and adjacent regions, by examining strong geomagnetic storms occurring between 2000 and 2023. The ground‐based Global Navigation Satellite System (GNSS) total electron content and data from a new TEC‐based ionospheric data assimilation system were used to analyze the characteristics of SED. For each derived SED events, we have quantified its morphology by employing a Gaussian function to parameterize key characteristics of the SED, such as the plume intensity, central longitude, and half‐width. A statistical analysis of SEDs was conducted for the first time to characterize their climatological features. We found that the SED distribution exhibits a higher peak intensity and a narrower width as geomagnetic activity strengthens. The peak intensity of SED has maximum values around the equinoxes in their seasonal distribution. Additionally, we observed a solar cycle dependence in the SED distribution, with more events occurring during the solar maximum and declining phases compared to the solar minimum. SED plumes exhibit a sub‐corotation feature with respect to the Earth, characterized by a westward drift speed between 50 and 400 m/s and a duration of 3–10 hr. These information advanced the current understanding of the spatial‐temporal variation of SED characteristics.
-
Abstract This study investigates the ionospheric total electron content (TEC) responses in the 2‐D spatial domain and electron density variations in the 3‐D spatial domain during the annular solar eclipse on 14 October 2023, using ground‐based Global Navigation Satellite System (GNSS) observations, a novel TEC‐based ionospheric data assimilation system (TIDAS), ionosonde measurements, and satellite in situ data. The main results are summarized as follows: (a) The 2‐D TEC responses exhibited distinct latitudinal differences. The mid‐latitude ionosphere exhibited a more substantial TEC decrease of 25%–40% along with an extended recovery time of 3–4 hr. In contrast, the equatorial and low‐latitude ionosphere experienced a smaller TEC reduction of 10%–25% and a faster recovery time of 20–50 min. The minimal eclipse effect was observed near the northern equatorial ionization anomaly crest region. (b) The ionospheric electron density variations during the eclipse were effectively reconstructed by TIDAS data assimilation in the 3‐D domain, providing important altitude information with validity. (c) The ionospheric electron density variations showed a notable altitude‐dependent feature. The eclipse led to a substantial electron density reduction of 30%–50%, with the maximum depletion occurring around the ionospheric F2‐layer peak height (hmF2) of 250–350 km. The post‐eclipse recovery of electron density exhibited a relatively slower pace near the F2‐layer peak height than that at lower and higher altitudes.
Free, publicly-accessible full text available March 1, 2025 -
Abstract This paper investigates the midlatitude ionospheric disturbances over the American/Atlantic longitude sector during an intense geomagnetic storm on 23 April 2023. The study utilized a combination of ground‐based observations (Global Navigation Satellite System total electron content and ionosonde) along with measurements from multiple satellite missions (GOLD, Swarm, Defense Meteorological Satellite Program, and TIMED/GUVI) to analyze storm‐time electrodynamics and neutral dynamics. We found that the storm main phase was characterized by distinct midlatitude ionospheric density gradient structures as follows: (a) In the European‐Atlantic longitude sector, a significant midlatitude bubble‐like ionospheric super‐depletion structure (BLISS) was observed after sunset. This BLISS appeared as a low‐density channel extending poleward/westward and reached ∼40° geomagnetic latitude, corresponding to an APEX height of ∼5,000 km. (b) Coincident with the BLISS, a dynamic storm‐enhanced density plume rapidly formed and decayed at local afternoon in the North American sector, with the plume intensity being doubled and halved in just a few hours. (c) The simultaneous occurrence of these strong yet opposite midlatitude gradient structures could be mainly attributed to common key drivers of prompt penetration electric fields and subauroral polarization stream electric fields. This shed light on the important role of storm‐time electrodynamic processes in shaping global ionospheric disturbances.
Free, publicly-accessible full text available March 1, 2025 -
Abstract Volcanic eruptions provide broad spectral forcing to the atmosphere and understanding the primary mechanisms that are relevant to explain the variety in waveform characteristics in the Ionosphere‐Thermosphere (IT) is still an important open question for the community. In this study, Global Navigation Satellite System (GNSS) Total Electron Content (TEC) data are analyzed and compared to simulations performed by the Global Ionosphere‐Thermosphere Model with Local Mesh Refinement (GITM‐R) for the first phase of the 2015 Calbuco eruption that occurred on 22 April. A simplified source representation and spectral acoustic‐gravity wave (AGW) propagation model are used to specify the perturbation at the lower boundary of GITM‐R at 100 km altitude. Two assumptions on the propagation structure, Direct Spherical (DS) and Ground Coupled (GC), are compared to the GNSS data and these modeling specifications show good agreement with different aspects of the observations for some waveform characteristics. Most notably, GITM‐R is able to reproduce the relative wave amplitude of AGWs as a function of radial distance from the vent, showing acoustic dominant forcing in the near field (<500 km) and gravity dominant forcing in the far‐field (>500 km). The estimated apparent phase speeds from GITM‐R simulations are consistent with observations with ∼10% difference from observation for both acoustic wave packets and a trailing gravity mode. The relevance of the simplifications made in the lower atmosphere to the simulated IT response is then discussed.
-
Abstract The generation of medium‐scale traveling ionospheric disturbances (MSTIDs) in the mid‐latitude F region ionosphere, particularly in the presence of sporadic E (Es) layers or geomagnetically conjugate features, has been the focus of extensive investigation using both observational and numerical modeling approaches. Recent observations have revealed the occurrence of nighttime MSTIDs over the continental US during storm conditions even without invoking the Es instability. While this phenomenon is considered to be electrified and likely associated with the Perkins instability, the influences of storm‐enhanced density (SED), electric fields, and winds on the excitation of nighttime MSTIDs remain a complicated issue and require further quantitative analysis. In this study, we develop a two‐dimensional numerical model of the nighttime ionospheric electrodynamics at midlatitudes using the ionospheric ion continuity equation and the electric field Poisson equation to investigate the characteristics of MSTIDs in the SED base region during storm conditions. We demonstrate that the magnetic inclination effect can explain the lower latitude preference of the MSTIDs during magnetic storms, while the development of MSTIDs is primarily influenced by intense storm electric fields under the background ionospheric condition of large density gradients associated with SED. However, the impact of neutral winds on the MSTIDs growth varies, depending on their specific direction determined by the strongly dynamic spatiotemporal variation of the thermosphere and ionosphere during storms. Therefore, the MSTIDs stormtime scenario results from a combination of multiple important factors.
-
Abstract This study investigates the global distribution of electron temperature enhancement observed by Defense Meteorological Satellite Program F16 satellite and its dependence on the season and solar activity for the solar maximum (2014) and minimum (2018) years during geomagnetic quiet times (maximum per day ap <10). Electron temperature enhancements occurred mainly over the North American‐Atlantic (260°–360°E) and Eurasia (0°–160°E) (Southern Oceania (80°–280°E)) sector in the Northern (Southern) Hemisphere and are prominent in the winter hemispheres and solar maximum year. They have obvious longitude characteristics. Interestingly, they could extend to geomagnetic equatorial regions in the North American‐Atlantic sector from high to low latitudes in the December Solstice, further crossed the magnetic equator, and merged into the Southern Hemisphere in 2014, where the maximum temperature reached ∼3500 K. Our analysis indicates that low‐energy electrons (<100 eV) associated with photoelectron from the conjugate sunlit hemisphere, can contribute to these enhancements. Furthermore, the local geomagnetic declination, magnetic equator position, and terminator position at magnetic conjugate points together can impact the global distribution of photoelectrons of different energies and therefore the electron temperature enhancement distribution. Other processes (including local electron density variation) may play certain roles as well.