This study presents a statistical analysis of the ionospheric irregularities and topside ionospheric scintillation at low and middle latitudes by using in situ electron density and upward-looking total electron content data measured by the Swarm constellation during 2014–2021. The main purpose of this study is to determine whether the phase scintillation could present similar seasonal, longitudinal, latitudinal, local time, and solar activity features as the in situ ionospheric irregularities do at low and middle latitudes, and how the irregularities affect the phase scintillation. The results are summarized as follows: (1) At low latitudes, the occurrence rate of equatorial plasma irregularities (EPIs) at the equinoxes and December solstice peaks before midnight, but during the June solstice, the EPIs mainly occur after midnight. The occurrence rate of EPIs has a positive correlation with solar activity. The distribution of topside scintillation occurrence is relatively consistent with EPIs, but during the June solstice, the scintillation occurrence rate remains at a very low level. (2) The midlatitude irregularities mainly occur after midnight, and their occurrence rate is negatively correlated with solar activity. Midlatitude irregularities mainly occur during the solstices, concentrated over the Pacific region during the June solstice and over the Pacific American sector during the December solstice. Especially, the distribution of midlatitude irregularities has hemispheric asymmetry, with a higher occurrence rate in the winter hemisphere. However, the occurrence of midlatitude scintillation is comparable in both hemispheres during the June solstice, and it concentrates in the southern hemisphere during the December solstice. (3) The EPIs concentrate more at the altitudes of Swarm A, while the midlatitude irregularities mainly occur at the altitudes of Swarm B.

Equatorial plasma bubbles (EPBs) occur frequently in low-latitude areas and have a non-negligible impact on navigation satellite signals. To systematically analyze the effects of a single EPB event on multi-frequency signals of GPS, Galileo, GLONASS, and BDS, all-sky airglow images over South China are jointly used to visually determine the EPB structure and depletion degree. The results reveal that scintillations, or GNSS signal fluctuations, are directly linked to EPBs and that the intensity of scintillation is positively correlated with the airglow depletion intensity. The center of the airglow depletion often corresponds to stronger GNSS scintillation, while the edge of the bubble, which is considered to have the largest density gradient, corresponds to relatively smaller scintillation instead. This work also systematically analyzes the responses of multi-constellation and multi-frequency signals to EPBs. The results show that the L2 and L5 frequencies are more susceptible than the L1 frequency is. For different constellations, Galileo’s signal has the best tracking stability during an EPB event compared with GPS, GLONASS, and BDS. The results provide a reference for dual-frequency signal selection in precise positioning or TEC calculation, that is, L1C and L2L for GPS, L1C and L5Q for Galileo, L1P and L2C for GLONASS, and L1P and L5P for BDS. Notably, BDS-2 is significantly weaker than BDS-3. And inclined geosynchronous orbit (IGSO) satellites have abnormal data error rates, which should be related to the special signal path trajectory of the IGSO satellite.

Previous efforts to retrieve thermospheric density using satellite payloads have been limited to a small number of satellites equipped with GNSS (Global Navigation Satellite System) receivers and accelerometers. These satellites are confined to a few orbital planes, and analysis can only be conducted after the data are processed and updated, resulting in sparse and delayed thermospheric density datasets. In recent years, the Starlink constellation, developed and deployed by SpaceX, has emerged as the world’s largest low Earth orbit (LEO) satellite constellation, with over 6000 satellites in operations as of October 2024. Through the strategic use of multiple orbital shells featuring various inclinations and altitudes, Starlink ensures continuous near-global coverage. Due to extensive coverage and frequent maneuvers, SpaceX has publicly released predicted ephemeris data for all Starlink satellites since May 2021, with updates approximately every 8 h. With the ephemeris data of Starlink satellites, we first apply a maneuver detection algorithm based on mean orbital elements to analyze their maneuvering behavior. The results indicate that Starlink satellites exhibit more frequent maneuvers during thermospheric disturbances. Then, we calculate the mechanical energy loss caused by non-conservative forces (primarily atmospheric drag) through precise dynamical models. The results demonstrate that, despite certain limitations in Starlink ephemeris data, the calculated mechanical energy loss still effectively captures thermospheric density variations during both quiet and disturbed geomagnetic periods. This finding is supported by comparisons with Swarm-B data, revealing that SpaceX incorporates the latest space environment conditions into its orbit extrapolation models during each ephemeris update. With a maximum lag of only 8 h, this approach enables near-real-time monitoring of thermospheric density variations using Starlink ephemeris.