A COSMIC-1/FORMOSAT-3 (Constellation Observing System for Meteorology, Ionosphere, and Climate-1 and Formosa Satellite Mission 3) follow-on mission, COSMIC-2/FORMOSAT-7, had been successfully launched into low-inclination orbits on 25 June 2019. COSMIC-2 has a significantly increased Signal-to-Noise ratio (SNR) compared to other Radio Occultation (RO) missions. This study summarized the initial assessment of COSMIC-2 data quality conducted by the NOAA (National Oceanic and Atmospheric Administration) Center for Satellite Applications and Research (STAR). We use validated data from other RO missions to quantify the stability of COSMIC-2. In addition, we use the Vaisala RS41 radiosonde observations to assess the accuracy and uncertainty of the COSMIC-2 neutral atmospheric profiles. RS41 is currently the most accurate radiosonde observation system. The COSMIC-2 SNR ranges from 200 v/v to about 2800 v/v. To see if the high SNR COSMIC-2 signals lead to better retrieval results, we separate the COSMIC-2–RS41 comparisons into different SNR groups (i.e., 0–500 v/v group, 500–1000 v/v group, 1000–1500 v/v group, 1500–2000 v/v group, and >2000 v/v group). In general, the COSMIC-2 data quality in terms of stability, precision, accuracy, and uncertainty of the accuracy is very compatible with those from COSMIC-1. Results show that the mean COSMIC-2–RS41 water vapor difference from surface to 5 km altitude for each SNR groups are equal to −1.34 g/kg (0–500 v/v), −1.17 g/kg (500–1000 v/v), −1.33 g/kg (1000–1500 v/v), −0.93 g/kg (1500–2000 v/v), and −1.52 g/kg (>2000 v/v). Except for the >2000 v/v group, the high SNR measurements from COSMIC-2 seem to improve the mean water vapor difference for the higher SNR group slightly (especially for the 1500–2000 v/v group) comparing with those from lower SNR groups.
Recently, the 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 its numerical weather prediction systems, 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 CWDP 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 RO data products for weather and climate applications. Spire RO data collected from commercial CubeSats are carefully compared with data from Formosa Satellite Mission 7–Constellation Observing System for Meteorology, Ionosphere, and Climate-2 (COSMIC-2), the fifth-generation European Centre for Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis (ERA5), and high-quality radiosonde data. The results demonstrate that, despite their generally lower Signal-Noise-Ratio (SNR), Spire RO data show a pattern of lowest penetration height similar to that of COSMIC-2. The Spire and COSMIC-2 penetration heights are between 0.6 and 0.8 km altitude over tropical oceans. Although using different GNSS RO receivers, the precision of Spire STRATOS receivers is of the same quality as those of the COSMIC-2 TriG (Global Positioning System—GPS, GALILEO, and GLObal NAvigation Satellite System—GLONASS) RO Receiver System (TGRS) receivers. Furthermore, the Spire and COSMIC-2 retrieval accuracies are quite comparable. We validate the Spire temperature and water vapor profiles by comparing them with collocated radiosonde observation (RAOB) 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 of <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 estimate the error covariance matrix for Spire, COSMIC-2, and KOMPSAT-5. The results show that the COSMIC-2 estimated error covariance values are slightly more significant than those from Spire over the oceans at the mid-latitudes (45°N–30°N and 30°S–45°S), which may be owing to COSMIC-2 SNR being relatively lower at those latitudinal zones.
In response to its users’ needs, the National Oceanic and Atmospheric Administration (NOAA) initiated reanalysis (RAN) of the Advanced Very High Resolution Radiometer (AVHRR) Global Area Coverage (GAC; 4 km) sea surface temperature (SST) data employing its Advanced Clear Sky Processor for Oceans (ACSPO) retrieval system. Initially, AVHRR/3 data from five NOAA and two Metop satellites from 2002 to 2015 have been reprocessed. The derived SSTs have been matched up with two reference SSTs—the quality controlled in situ SSTs from the NOAA in situ Quality Monitor (iQuam) and the Canadian Meteorological Centre (CMC) L4 SST analysis—and analyzed in the NOAA SST Quality Monitor (SQUAM) online system. The corresponding clear-sky ocean brightness temperatures (BT) in AVHRR bands 3b, 4 and 5 (centered at 3.7, 11, and 12 µm, respectively) have been compared with the Community Radiative Transfer Model simulations in another NOAA online system, Monitoring of Infrared Clear-sky Radiances over Ocean for SST (MICROS). For some AVHRRs, the time series of “AVHRR minus reference” SSTs and “observed minus model” BTs are unstable and inconsistent, with artifacts in the SSTs and BTs strongly correlated. In the official “Reanalysis version 1” (RAN1), data from only five platforms—two midmorning (NOAA-17 and Metop-A) and three afternoon (NOAA-16, -18 and -19)—were included during the most stable periods of their operations. The stability of the SST time series was further improved using variable regression SST coefficients, similarly to how it was done in the NOAA/NASA Pathfinder version 5.2 (PFV5.2) dataset. For data assimilation applications, especially those blending satellite and in situ SSTs, we recommend bias-correcting the RAN1 SSTs using the newly developed sensor-specific error statistics (SSES), which are reported in the product files. Relative performance of RAN1 and PFV5.2 SSTs is discussed. Work is underway to improve the calibration of AVHRR/3s and extend RAN time series, initially back to the mid-1990s and later to the early 1980s.
The radio occultation technique provides stable atmospheric measurements that can work as a benchmark for calibrating and validating satellite-sounding data. Launched on 25 June 2019, the Constellation Observing System for Meteorology, Ionosphere, and Climate 2 and Formosa Satellite Mission 7 (COSMIC-2/FORMOSAT-7) are expected to produce about 5000 high-quality RO observations daily over the tropics and subtropics. COSMIC-2 constellation consists of 6 Low Earth Orbit (LEO) satellites in 24° inclination orbits at 720 km altitude and distributed mainly between 45°N to 45°S. The COSMIC-2 observations have uniform temporal coverage between 30°N to 30°S. This paper presents an independent inversion algorithm to invert COSMIC-2 geometry and phase data to bending angle and refractivity. We also investigate the quality of Global Navigation Satellite System (GNSS) and LEO position vectors derived from the UCAR COSMIC Data Analysis and Archive Center (CDAAC). The GNSS and LEO position vectors are stable with LEO position variations < 1.4 mm/s. The signal-to-noise ratio (SNR) on the L1 band ranges from 300–2600 v/v with a mean of 1600 v/v. The inversion algorithm developed at NOAA Center for Satellite Applications and Research (STAR) uses the Full Spectrum Inversion (FSI) method to invert COSMIC-2 geometry and phase data to bending angle and refractivity profiles. The STAR COSMIC-2 bending angle and refractivity profiles are compared with in situ radiosonde, the current COSMIC-2 products derived from CDAAC, and the collocated European Center for Medium-Range Weather Forecasts (ECMWF) climate reanalysis data ERA5. The mean bias at 8–40 km altitude among the UCAR, ERA5, and STAR is <0.1% for both bending and refractivity, with a standard deviation in the range of 1.4–2.3 and 0.9–1.1% for bending angles refractivity, respectively. In the lowest 2 km, the RO bias relative to ERA-5 shows a strong latitudinal and SNR dependence.
Global Navigation Satellite System (GNSS) Radio Occultation (RO) data play an essential role in improving numerical weather prediction (NWP) and monitoring climate change. The NOAA Commercial RO Purchase Program (CDP) purchased RO data provided by Spire Global Inc. To ensure the data quality from Spire Global Inc. is consistent with other RO missions, we need to quantify their accuracy and retrieval uncertainty carefully. In this work, Spire Wet Profile (wet temperature profile) data from 7 September 2021 to 31 October 2022, processed by the University Corporation for Atmospheric Research (UCAR), and COSMIC-2 (Constellation Observing System for Meteorology, Ionosphere, and Climate-2/Formosa Satellite Mission 7) data are evaluated through comparison with NOAA-20 Advanced Technology Microwave Sounder (ATMS) microwave sounder measurements and collocated RS41 radiosonde measurements. Through the Community Radiative Transfer Model (CRTM) simulation, we convert the Spire and COSMIC-2 RO retrievals to ATMS brightness temperature (BT) at sounding channels CH07 to CH14 (temperature channels), with weighting function peak heights from 8 km to 35 km, and CH19 to CH22 (water vapor channels), with weighting function peak heights ranging from 3.2 km to 6.7 km, and compare the simulations with the collocated NOAA-20 ATMS measurements over ocean. Using ATMS observations as references, Spire and COSMIC-2 BTs agree well with ATMS within 0.07 K for CH07-14 and 0.20 K for CH19-22. The trends between Spire and COSMIC-2 are consistent within 0.07 K/year over the oceans for ATMS CH07-CH13 and CH19-22, indicating that Spire/COSMIC-2 wet profiles are, in general, compatible with each other over oceans. The RO retrievals and RS41 radiosonde observation (RAOB) comparison shows that above 0.2 km altitude, RS41 RAOB matches Spire/COSMIC-2 temperature profiles well with a temperature difference of <0.13 K, and the trends between Spire and COSMIC-2 are consistent within 0.08 K/year over land, indicating that Spire/COSMIC-2 wet profiles are overall compatible with each other through RS41 RAOB measurements over land. In addition, the consistency of Spire and COSMIC-2 based on different latitude intervals, local times, and signal-to-noise ratios (SNRs) through ATMS was evaluated. The results show that the performance of Spire is comparable to COSMIC-2, even though COSMIC-2 has a higher SNR. The high quality of RO profiles from Spire is expected to improve short- and medium-range global numerical weather predictions and help construct consistent climate temperature records.
The global navigation satellite system (GNSS) radio occultation (RO) is becoming an essential component of National Oceanic and Atmospheric Administration (NOAA) observation systems. The constellation observing system for meteorology, ionosphere, and climate (COSMIC) 2 mission and the Formosa satellite mission 7, a COSMIC follow-on mission, is now the NOAA’s backbone RO mission. The NOAA’s dedicated GNSS RO SAtellite processing and science Application Center (RO-SAAC) was established at the Center for Satellite Applications and Research (STAR). To better quantify how the observation uncertainty from clock error and geometry determination may propagate to bending angle and refractivity profiles, STAR has developed the GNSS RO data processing and validation system. This study describes the COSMIC-2 neutral atmospheric temperature and moisture profile inversion algorithms at STAR. We used RS41 and ERA5, and UCAR 1D-Var products (wetPrf2) to validate the accuracy and uncertainty of the STAR 1D-Var thermal profiles. The STAR-RS41 temperature differences are less than a few tenths of 1 K from 8 km to 30 km altitude with a standard deviation (std) of 1.5–2 K. The mean STAR-RS41 water vapor specific humidity difference and the standard deviation are −0.35 g/kg and 1.2 g/kg, respectively. We also used the 1D-Var-derived temperature and water vapor profiles to compute the simulated brightness temperature (BTs) for advanced technology microwave sounder (ATMS) and cross-track infrared sounder (CrIS) channels and compared them to the collocated ATMS and CrIS measurements. The BT differences of STAR COSMIC-2-simulated BTs relative to SNPP ATMS are less than 0.1 K over all ATMS channels.