Global Ionospheric Electron Density from GNSS-POD Limb Measurements

Abstract

GNSS-LEO radio links from Precise Orbital Determination (POD) and Radio Occultation (RO) antennas have been used increasingly in studying and monitoring the global ionospheric electron density (Nₑ). In this study we developed an optimal estimation (OE) method to retrieve Ne profiles from the slant total electron content (hTEC) measurements acquired by the GNSS-POD links at negative elevation angles. The hmF2 and NmF2 from the OE retrieval are validated against ground-based ionosondes and radar observations, showing generally good agreements in NmF2 from all sites. Nighttime hmF2 measurements tend to agree better than the daytime when the ionosonde heights tend to be slightly lower. The OE algorithm has been applied to all GNSS-POD data acquired from the COSMIC-1 (2006-2019), COSMIC-2 (2019- present), and Spire (2019-present) constellations, showing a consistent ionospheric Ne morphology. A detailed analysis of the frequency-wavenumber spectra is made for the Ne variability at different heights. In the lower ionosphere (~150 km) we found significant spectral power in DE1, DW6, DW4, SW5, and SE4 wave components, in addition to well-known DW1, SW2 and DE3 waves. In the upper ionosphere (~450 km), additional wave components are still present, including DE4, DW4, DW6, SE4 and SW4. The co-exist of eastward and westward propagating WN4 components implies the presence of a stationary planetary wave4 (sPW4), as suggested by other earlier studies. GNSS-LEO radio links from Precise Orbital Determination (POD) and Radio Occultation (RO) antennas have been used increasingly in studying and monitoring the global ionospheric electron density (Ne). In this study we developed an optimal estimation (OE) method to retrieve Ne profiles from the slant total electron content (hTEC) measurements acquired by the GNSS-POD links at negative elevation angles. The hmF2 and NmF2 from the OE retrieval are validated against ground-based ionosondes and radar observations, showing generally good agreements in NmF2 from all sites. Nighttime hmF2 measurements tend to agree better than the daytime when the ionosonde heights tend to be slightly lower. The OE algorithm has been applied to all GNSS-POD data acquired from the COSMIC-1 (2006-2019), COSMIC-2 (2019- present), and Spire (2019-present) constellations, showing a consistent ionospheric Ne morphology. A detailed analysis of the frequency-wavenumber spectra is made for the Ne variability at different heights. In the lower ionosphere (~150 km) we found significant spectral power in DE1, DW6, DW4, SW5, and SE4 wave components, in addition to well-known DW1, SW2 and DE3 waves. In the upper ionosphere (~450 km), additional wave components are still present, including DE4, DW4, DW6, SE4 and SW4. The co-exist of eastward and westward propagating WN4 components implies the presence of a stationary planetary wave4 (sPW4), as suggested by other earlier studies.