The analytical description of the topside ionosphere included in the NeQuick model is studied in detail. First, the modeled scale height behavior is analyzed at infinity and for the lowest part of the topside region; in the latter case, the analysis is done through an expansion in Taylor series near the F2-layer peak. Moreover, the significant influence of the three NeQuick topside parameters in the modeling of the topside profile is carefully investigated and, in particular, it is shown that for the lowest part of the topside the model assumes a linearly increasing trend of the topside scale height. Second, the topside formulation is inverted to derive a fully analytical expression of the topside scale height as a function of the electron density and F2-layer peak parameters. This expression has been applied to a selected and very reliable dataset of COSMIC/FORMOSAT-3 (Constellation Observing System for Meteorology, Ionosphere, and Climate) Radio Occultation profiles. Statistical analyses strongly support the hypothesis made by NeQuick regarding the linear trend of the topside scale height for the lowest part of the topside.

This paper presents a review on the PECASUS service, which provides advisories on enhanced space weather activity for civil aviation. The advisories are tailored according to the Standards and Recommended Practices of the International Civil Aviation Organization (ICAO). Advisories are disseminated in three impact areas: radiation levels at flight altitudes, GNSS-based navigation and positioning, and HF communication. The review, which is based on the experiences of the authors from two years of running pilot ICAO services, describes empirical models behind PECASUS products and lists ground- and space-based sensors, providing inputs for the models and 24/7 manual monitoring activities. As a concrete example of PECASUS performance, its products for a post-storm ionospheric F2-layer depression event are analyzed in more detail. As PECASUS models are particularly tailored to describe F2-layer thinning, they reproduce observations more accurately than the International Reference Ionosphere model (IRI(STORM)), but, on the other hand, it is recognized that the service performance is much affected by the coverage of its input data. Therefore, more efforts will be directed toward systematic measuring of the availability, timeliness and quality of the data provision in the next steps of the service development.

The series of X and M class flares and associated coronal mass ejections that occurred on the first days of September 2017 induced significant perturbations on the low‐latitude ionospheric electrodynamics. On 8 September in the Indian sector, the storm caused a severe modification of the equatorial electrojet (EEJ) with a consequent variation of the ionospheric structuring and dynamics. In our analysis, we propose an original method to isolate and identify EEJ variations from geomagnetic data and we detect the presence of equatorial plasma bubbles (EPB) from L‐band total electron content (TEC) data in order to understand their movement. Our results provide evidence of independent EPBs appearance freshly generated and inherited from a migrating plasma structure. The EPB (or EPBs) occurring in the south of India is/are freshly generated just above the magnetic equator, and is/are likely triggered by the sudden increase of EEJ just before the local sunset, acting as a pre‐reversal enhancement. The EPB appearing in the North‐East Indian region is associated with a migrating structure, resulting in a northward movement with a velocity of about 650 m/s, possibly testifying the passage of a large‐scale traveling ionospheric disturbance. The occurrence of severe post‐sunset scintillations in the northeastern sector suggests a possible cascade process forming small‐scale irregularities from the migrating EPB.

In this work, we investigate various types of ionospheric disturbances observed over Europe following the earthquake that occurred in Turkey on 6 February 2023. By combining observations from Doppler sounding systems, ionosondes, and GNSS receivers, we are able to discern different types of disturbances, propagating with different velocities and through different mechanisms. We can detect co‐seismic ionospheric disturbances close to the epicenter, as well as ionospheric signatures of acoustic waves propagating as a consequence of propagating seismic waves. Unlike the vast majority of past ionospheric co‐seismic disturbance studies that are primarily based on Total Electron Content variations, reflecting disturbances propagating around the F‐region peak, the focus of the present study is the manifestation of disturbances at different ionospheric altitudes by exploiting complementary ionospheric remote sensing techniques. This is particularly highlighted through ionospheric earthquake‐related signatures established as specific ionogram deformations known as multiple‐cusp signatures which appear as additional cusps at the base of the F‐region attributed to electron density irregularities generated by Rayleigh surface waves that generate acoustic waves propagating up to the ionosphere. Therefore this study underlines the advantage that multi‐instrument investigations offer in identifying the propagation of earthquake‐related ionospheric disturbances at different ionospheric altitudes and distances from the earthquake epicenter.

The bottomside thickness and shape (B0 and B1) are the critical key elements for depicting a realistic electron density profile in the International Reference Ionosphere (IRI) model. We investigated their longitudinal variability using a large database of FORMOSAT‐3/COSMIC (Constellation Observing System for Meteorology, Ionosphere, and Climate) electron density profiles covering the period 2006–2015 from which these parameters are derived using a simple analytic Chapman least squares fitting. The parameters are then compared with the IRI‐2012 submodels, that is, Bil‐2000, Gul‐1987, and ABT‐2009 during different seasons and solar activity. Besides apparent adherence of B0 to the equatorial region, the opposite occurrence pattern of maximum B0 to that of NmF2 is realized, manifesting semiannual distribution and hemispheric asymmetry. Unlike NmF2, the presence of single anomaly crest in B0 in all seasons except March equinox and their leading local time stamps are among the key observations. In spite of regular representation of IRI submodels during equinox, the performance of Gul‐1987 seems to be inferior during solstice with respect to geomagnetic field and equatorial electrodynamics, thereby depicting an erroneous distribution of B0. The striking wave‐like feature in B0 longitudinal variability is also revealed indicating its possible connection with the nonmigrating diurnal/semidiurnal tides. Moreover, the solar activity effect on B0 appears to be relatively higher over the Asian and Pacific longitudes (60–150°E). Although speculations are made over B1 variability, its complete longitudinal characteristics are indistinguishable in our results, except the subordinate magnitude as compared to IRI estimations. The important findings from this study may reinforce the understanding of bottomside variability and future improvements in IRI.

The present study analyzes the global ionospheric response to the intense geomagnetic storm of 10–11 October 2024 (SYM—H minimum of −346 nT), using observations from COSMIC—2 and Swarm satellites, GNSS TEC, and Digisondes. Significant uplift of the F-region was observed across both Hemispheres on the dayside, primarily driven by equatorward thermospheric winds and prompt penetration electric fields (PPEFs). However, this uplift did not correspond with increases in foF2 due to enhanced molecular nitrogen-promoting recombination in sunlit regions and the F2 peak rising beyond the COSMIC—2 detection range. In contrast, in the Southern Hemisphere nightside ionosphere exhibited pronounced Ne depletion and low hmF2 values, attributed to G-conditions and thermospheric composition changes caused by storm-time circulation. Strong vertical plasma drifts exceeding 100 m/s were observed during both the main and recovery phases, particularly over Ascension Island, driven initially by southward IMF—Bz-induced PPEFs and later by disturbance dynamo electric fields (DDEFs) as IMF—Bz turned northward. Swarm data revealed a poleward expansion of the Equatorial Ionization Anomaly (EIA), with more pronounced effects in the Southern Hemisphere due to seasonal and longitudinal variations in ionospheric conductivity. Additionally, the storm excited Large-Scale Travelling Ionospheric Disturbances (LSTIDs), triggered by thermospheric perturbations and electrodynamic drivers, including PPEFs and DDEFs. These disturbances, along with enhanced westward thermospheric wind and altered zonal electric fields, modulated ionospheric irregularity intensity and distribution. The emergence of anti-Sq current systems further disrupted quiet-time electrodynamics, promoting global LSTID activity. Furthermore, storm-induced equatorial plasma bubbles (EPBs) were observed over Southeast Asia, initiated by enhanced PPEFs during the main phase and suppressed during recovery, consistent with super EPB development mechanisms.