skip to main content


Title: Comparison of COSMIC and COSMIC-2 Radio Occultation Refractivity and Bending Angle Uncertainties in August 2006 and 2021
We compare the random error statistics (uncertainties) of COSMIC (Constellation Observing System for Meteorology, Ionosphere and Climate, C1) and COSMIC-2 (C2) radio occultation (RO) bending angles and refractivities for the months of August 2006 and 2021 over the tropics and subtropics using the three-cornered hat method. The uncertainty profiles are similar for the two RO missions in the troposphere. However, a higher percentage of C2 profiles reach close to the surface in the moisture-rich tropics, an advantage of the higher signal-to-noise ratio (SNR) in C2. C2 uses signals from both GPS (Global Positioning System) and GLONASS Global Navigation System Satellites (GNSS). The GPS occultations show smaller uncertainties in the stratosphere and lower mesosphere (30–60 km) than the GLONASS occultations, a result of more accurate GPS clocks. Therefore, C2 (GPS) uncertainties are smaller than C1 uncertainties between 30–60 km while the C2 (GLONASS) uncertainties are larger than those of C1. The uncertainty profiles vary with latitude at all levels. We find that horizontal gradients in temperature and water vapor, and therefore refractivity, are the major cause of uncertainties in the tropopause region and troposphere through the violation of the assumption of spherical symmetry in the retrieval of bending angles and refractivity.  more » « less
Award ID(s):
2054356
NSF-PAR ID:
10436586
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Atmosphere
Volume:
13
Issue:
5
ISSN:
2073-4433
Page Range / eLocation ID:
790
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Estimation of uncertainties (random error statistics) of radio occultation (RO) observations is important for their effective assimilation in numerical weather prediction (NWP) models. Average uncertainties can be estimated for large samples of RO observations and these statistics may be used for specifying the observation errors in NWP data assimilation. However, the uncertainties of individual RO observations vary, and so using average uncertainty estimates will overestimate the uncertainties of some observations and underestimate those of others, reducing their overall effectiveness in the assimilation. Several parameters associated with RO observations or their atmospheric environments have been proposed to estimate individual RO errors. These include the standard deviation of bending angle (BA) departures from either climatology in the upper stratosphere and lower mesosphere (STDV) or the sample mean between 40 and 60 km (STD4060), the local spectral width (LSW), and the magnitude of the horizontal gradient of refractivity (|∇HN|). In this paper we show how the uncertainties of two RO datasets, COSMIC-2 and Spire BA, as well as their combination, vary with these parameters. We find that the uncertainties are highly correlated with STDV and STD4060 in the stratosphere, and with LSW and |∇HN| in the lower troposphere. These results suggest a hybrid error model for individual BA observations that uses an average statistical model of RO errors modified by STDV or STD4060 above 30 km, and LSW or |∇HN| below 8 km.

    Significance Statement

    These results contribute to the understanding of the sources of uncertainties in radio occultation observations. They could be used to improve the effectiveness of these observations in their assimilation into numerical weather prediction and reanalysis models by improving the estimation of their observational errors.

     
    more » « less
  2. Abstract FORMOSAT-3/COSMIC (F3/C) constellation of six micro-satellites was launched into the circular low-earth orbit at 800 km altitude with a 72-degree inclination angle on 15 April 2006, uniformly monitoring the ionosphere by the GPS (Global Positioning System) Radio Occultation (RO). Each F3/C satellite is equipped with a TIP (Tiny Ionospheric Photometer) observing 135.6 nm emissions and a TBB (Tri-Band Beacon) for conducting ionospheric tomography. More than 2000 RO profiles per day for the first time allows us globally studying three-dimensional ionospheric electron density structures and formation mechanisms of the equatorial ionization anomaly, middle-latitude trough, Weddell/Okhotsk Sea anomaly, etc. In addition, several new findings, such as plasma caves, plasma depletion bays, etc., have been reported. F3/C electron density profiles together with ground-based GPS total electron contents can be used to monitor, nowcast, and forecast ionospheric space weather. The S4 index of GPS signal scintillations recorded by F3/C is useful for ionospheric irregularities monitoring as well as for positioning, navigation, and communication applications. F3/C was officially decommissioned on 1 May 2020 and replaced by FORMOSAT-7/COSMIC-2 (F7/C2). F7/C2 constellation of six small satellites was launched into the circular low-Earth orbit at 550 km altitude with a 24-degree inclination angle on 25 June 2019. F7/C2 carries an advanced TGRS (Tri Gnss (global navigation satellite system) Radio occultation System) instrument, which tracks more than 4000 RO profiles per day. Each F7/C2 satellite also has a RFB (Radio Reference Beacon) on board for ionospheric tomography and an IVM (Ion Velocity Meter) for measuring ion temperature, velocity, and density. F7/C2 TGRS, IVM, and RFB shall continue to expand the F3/C success in the ionospheric space weather forecasting. 
    more » « less
  3. Radio occultation (RO) can provide high-vertical-resolution thermodynamic soundings of the planetary boundary layer (PBL). However, sharp moisture gradients and strong temperature inversion lead to large gradients in refractivity N and often cause ducting. Ducting results in systematically negative RO N biases resulting from a nonunique Abel inversion problem. Using 8 years (2006–13) of Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) RO soundings and collocated European Centre for Medium-Range Weather Forecasts interim reanalysis (ERA-I) data, we confirm that the large lower-tropospheric negative N biases are mainly located in the subtropical eastern oceans and we quantify the contribution of ducting for the first time. The ducting-contributed N biases in the northeast Pacific Ocean (160°–110°W; 15°–45°N) are isolated from other sources of N biases using a two-step geometric-optics simulation. Negative bending angle biases in this region are also observed in COSMIC RO soundings. Both the negative refractivity and bending angle biases in COSMIC soundings mainly lie below ~2 km. Such bending angle biases introduce N biases that are in addition to those caused by ducting. Following the increasing PBL height from the southern California coast westward to Hawaii, centers of maxima bending angles and N biases tilt southwestward. In areas where ducting conditions prevail, ducting is the major cause of the RO N biases. Ducting-induced N biases with reference to ERA-I compose over 70% of the total negative N biases near the southern California coast, where strongest ducting conditions prevail, and decrease southwestward to less than 20% near Hawaii.

     
    more » « less
  4. Abstract. Current climate models have difficulty representing realistic wave–mean flow interactions, partly because the contribution from waves with fine vertical scales is poorly known. There are few direct observations of these waves, and most models have difficulty resolving them. This observational challenge cannot be addressed by satellite or sparse ground-based methods. The Strateole-2 long-duration stratospheric superpressure balloons that float with the horizontal wind on constant-density surfaces provide a unique platform for wave observations across a broad range of spatial and temporal scales. For the first time, balloon-borne Global Navigation Satellite System (GNSS) radio occultation (RO) is used to provide high-vertical-resolution equatorial wave observations. By tracking navigation signal refractive delays from GPS satellites near the horizon, 40–50 temperature profiles were retrieved daily, from balloon flight altitude (∼20 km) down to 6–8 km altitude, forming an orthogonal pattern of observations over a broad area (±400–500 km) surrounding the flight track. The refractivity profiles show an excellent agreement of better than 0.2 % with co-located radiosonde, spaceborne COSMIC-2 RO, and reanalysis products. The 200–500 m vertical resolution and the spatial and temporal continuity of sampling make it possible to extract properties of Kelvin waves and gravity waves with vertical wavelengths as short as 2–3 km. The results illustrate the difference in the Kelvin wave period (20 vs. 16 d) in the Lagrangian versus ground-fixed reference and as much as a 20 % difference in amplitude compared to COSMIC-2, both of which impact estimates of momentum flux. A small dataset from the extra Galileo, GLONASS, and BeiDou constellations demonstrates the feasibility of nearly doubling the sampling density in planned follow-on campaigns when data with full equatorial coverage will contribute to a better estimate of wave forcing on the quasi-biennial oscillation (QBO) and improved QBO representation in models. 
    more » « less
  5. Assimilation of radio occultation (RO) observations into numerical weather prediction (NWP) models has improved forecasts, where RO is typically one of the top five observational systems contributing to forecast accuracy. By measuring the phase delay of radio waves traversing Earth’s atmosphere between global positioning system (GPS) and low-Earth orbiting satellites, RO obtains quasi-vertical profiles of bending angles (BA) of the radio waves’ trajectories. BA are the RO observation most often assimilated into NWP models, but since they are not computed or analyzed in the models, they must be computed from model data using a forward model. First, model refractivity N is computed from variables specified on the model’s vertical levels, then using the Abel integral, BA are computed from N. The forward model requires vertical differentiation of N, and accurate differentiation requires vertical interpolation of N between model levels. The interpolation results in errors, which then propagate through the forward model to produce BA errors. In this study, we investigate the sensitivity of forward-modeled BA to five different methods of vertical interpolation of N between model levels to determine a method that minimizes interpolation errors. We use RO-observed N to isolate the interpolation errors and determine an accurate method that can be applied to any NWP model. Of the five methods investigated, the log-spline interpolation reduces N and BA errors the most, regardless of the vertical resolution of the model grid.

     
    more » « less