Verification and validation of the COSMIC-2 excess phase and bending angle algorithms for data quality assurance at STAR Remote Sensing

In recent years, Global Navigation Satellite System (GNSS) radio occultation (RO) has become a critical observation system for global operational numerical weather prediction. Constellation Observing System for Meteorology, Ionosphere, Climate (COSMIC) 2 (COSMIC-2) has been a backbone RO mission for NOAA. NOAA also began to purchase RO data from commercial sources in 2020. To ensure the consistent quality of RO data from different sources, NOAA Center for Satellite Applications and Research (STAR) has developed capabilities to process all available RO data from different missions. This paper describes the STAR RO processing systems which convert the pseudorange and carrier phase observations to excess phases and bending angles (BAs). We compared our COSMIC-2 data products with those processed by the University Corporation for Atmospheric Research (UCAR) COSMIC Data Analysis and Archive Center (CDAAC).We processed more than twelve thousand COSMIC-2 occultation profiles. Our results show that the excess phase difference between UCAR and STAR is within a few centimeters at high altitudes, although the difference increases towards the lower atmosphere. The BA profiles derived from the excess phase are consistent with UCAR. The mean relative BA differences at impact height from 10 to 30 km are less than 0.1% for GLObal NAvigation Satellite System (GLONASS) L2C signals and Global Positioning System (GPS) L2C and L2P signals. The standard deviations are 1.15%, 1.15%, and 1.32% for GLONASS L2C signal and for GPS L2C and L2P signals, respectively. The BA profiles agree with those derived from European Center for Medium-rangeWeather Forecast (ECMWF) reanalysis version 5 (ERA5). The Signal-to-Noise-Ratio (SNR) plays an essential role in the processing. The STAR BA profiles with higher L1 SNRs (L1 at 80 km) tend to yield more consistent results than those from UCAR, with a negligible difference and a smaller deviation than lower SNR profiles. Profiles with lower SNR values tend to show a more significant standard deviation towards the surface during the open-loop stage in the lower troposphere than those of higher SNR. We also found that the different COSMIC-2 clock solutions could contribute to the significant relative BA difference at high altitudes; however, it has little effect on the lower troposphere comparisons given larger BA values. Grant no. NA19NES4320002