Comparison of a Bottom-Up GNSS Radio Occultation Method to Measure D- and E-Region Electron Densities with Ionosondes and FIRI

Abstract

High-frequency skywave propagation can be heavily impacted by D- and E-region dynamics requiring accurate global measurements to optimize performance. A standard measurement technique is to use ionosondes, but they are unable to measure below 1 MHz and are only available at a limited number of land-locked sites around the globe. In contrast, the Global Navigation Satellite System radio occultation (GNSS-RO) bottom-up method is a new approach specifically designed to generate electron density profiles in the D- and E- region ionosphere. It takes advantage of satellite constellations that currently provide over 20,000 daily measurements and global coverage. In this paper, GNSS-RO profiles were compared against ionosonde profiles at four sites covering a wide latitudinal range, and FIRI modeled profiles corresponding to the same latitude and local solar time. This comparison was completed using daytime profiles when sporadic-E (Es) was not present. The average GNSS-RO profile is found to be a few kilometers higher in altitude than the ionosonde profiles at the minimum frequency, fmin. When the ionosonde profiles are shifted so that the altitudes match at fmin, they are in good agreement up to the E-region peak altitude, hmE. Below fmin, the GNSS-RO profile is in good agreement with the FIRI profile, indicating that the profiles can measure the D- to E- transition region. The frequency of the E-region peak, foE, showed general agreement between the GNSS-RO and ionosonde measurements; however, the hmE agreement was weaker and the GNSS-RO profiles tend to have an hmE in a narrow altitude range for all profiles. Virtual heights were simulated for the GNSS-RO profiles using a numerical ray tracer for direct comparison with ionosonde observations, which showed agreement for many of the virtual heights near fmin, but also indicated a positive bias in the GNSS-RO virtual heights that may be due to low foE or elevated hmE estimates. For a quiet ionosphere, the shifted GNSS-RO electron density profiles were a good match for both measured ionosonde profiles and modeled FIRI profiles and the method is capable of providing global coverage of the D- and E-regions. Future work will require more data for seasonal and morning–afternoon comparisons as well as comparisons for the disturbed ionosphere when the sporadic-E layer is present.