Recently, NOAA has included GNSS (Global Navigation Satellite System) Radio Occultation (RO) data as one of the crucial long-term observables for weather and climate applications. To include more GNSS RO data in the numerical weather prediction system, the NOAA Commercial Weather Data Pilot program (CWDP) started to explore the commercial RO data available on the market. After two rounds of pilot studies, the CDWP decided to award the first Indefinite Delivery Indefinite Quantity (IDIQ) contract to GeoOptics and Spire Incs. in 2020. This study examines the quality of Spire data products for weather and climate applications. Spire RO data are collected from commercial CubeSats through careful comparison with the data from Formosa Satellite Mission 7–Constellation Observing System for Meteorology, Ionosphere, and Climate-2 (COSMIC-2), ERA-5, and high-quality radiosonde data. The results demonstrated that although with lower Signal-Noise-Ratio (SNR) in general, the pattern of the lowest penetration height for Spire is similar to those for COSMIC-2. The Spire and COSMIC-2 penetrate heights are between 0.6 and 0.8 km altitude at the tropical oceans. Although using different GNSS RO receivers, the precision of Spire STRASP receivers is of the same quality as those of COSMIC-2 Global Positioning System - GPS, GALILEO, and GLObal NAvigation Satellite System – GLONASS (TGRS) receivers. The retrieval accuracy from Spire is very compatible with those from COSMIC-2. We validated Spire temperature and water vapor profiles by comparing them with collocated radiosonde data. Generally, over the height region between 8 km and 16.5 km, the Spire temperature profiles match those from RS41 RAOB very well with temperature biases < 0.02 K. Over the height range from 17.8 to 26.4 km, the temperature differences are ~-0.034 K with RS41 RAOB being warmer. We also estimated the error covariance matrix for Spire, COSMIC-2, and KOMPSAT-5. Results showed that the COSMIC-2 estimated error covariance values are slightly more significant over the oceans at the mid-latitudes (45oN-30oN and 30oS-45oS), which may also be owing to COSMIC-2 SNR being lower at those latitudinal zones.