Analysis of GPS radio occultation data from the FORMOSAT-3/COSMIC and Metop/GRAS missions at CDAAC

This study investigates the noise level and mission-to-mission stability of Global Positioning System (GPS) radio occultation (RO) neutral atmospheric bending angle data at the UCAR COSMIC Data Analysis and Archive Center (CDAAC). Data are used from two independently developed RO instruments currently flying in orbit on the FORMOSAT-3/COSMIC (F3C) and Metop/GRAS (GNSS Receiver for Atmospheric Sounding) missions. The F3C 50 Hz RO data are post-processed with a single-difference excess atmospheric phase algorithm, and the Metop/GRAS 50 Hz closed loop and raw sampling (down-sampled from 1000 Hz to 50 Hz) data are processed with a zero-difference algorithm. The standard deviations of the F3C and Metop/GRAS bending angles from climatology between 60 and 80 km altitude from June–December 2009 are approximately 1.78 and 1.13 μrad, respectively. The F3C standard deviation reduces significantly to 1.44 μrad when single-difference processing uses GPS satellites on the same side of the spacecraft. The higher noise level for F3C bending angles can be explained by additional noise from the reference link phase data that are required with single-difference processing. The F3C and Metop/GRAS mean bending angles differences relative to climatology during the same six month period are statistically significant and have values of −0.05 and −0.02 μrad, respectively. A comparison of ~13 500 collocated F3C and Metop/GRAS bending angle profiles over this six month period shows a similar mean difference of ~0.02 ± 0.02 μrad between 30 and 60 km impact heights that is marginally significant. The observed mean difference between the F3C and Metop/GRAS bending angles of ~0.02–0.03 μrad is quite small and illustrates the high degree of re-produceability and mission independence of the GPS RO data at high altitudes. Collocated bending angles between two F3C satellites from early in the mission differ on average by up to 0.5% near the surface due to systematically lower signal-to-noise ratio for one of the satellites. Results from F3C and Metop/GRAS differences in the lower troposphere suggest the Metop/GRAS bending angles are negatively biased compared to F3C with a maximum of several percents near the surface in tropical regions. This bias is related to different tracking depths (deeper in F3C) and data gaps in Metop/GRAS which make it impossible to process the data from both missions in exactly the same way

Quantification of structural uncertainty in climate data records from GPS radio occultation

Global Positioning System (GPS) radio occultation (RO) has provided continuous observations of the Earth's atmosphere since 2001 with global coverage, all-weather capability, and high accuracy and vertical resolution in the upper troposphere and lower stratosphere (UTLS). Precise time measurements enable long-term stability but careful processing is needed. Here we provide climate-oriented atmospheric scientists with multicenter-based results on the long-term stability of RO climatological fields for trend studies. We quantify the structural uncertainty of atmospheric trends estimated from the RO record, which arises from current processing schemes of six international RO processing centers, DMI Copenhagen, EUM Darmstadt, GFZ Potsdam, JPL Pasadena, UCAR Boulder, and WEGC Graz. Monthly-mean zonal-mean fields of bending angle, refractivity, dry pressure, dry geopotential height, and dry temperature from the CHAMP mission are compared for September 2001 to September 2008. We find that structural uncertainty is lowest in the tropics and mid-latitudes (50° S to 50° N) from 8 km to 25 km for all inspected RO variables. In this region, the structural uncertainty in trends over 7 yr is <0.03% for bending angle, refractivity, and pressure, <3 m for geopotential height of pressure levels, and <0.06 K for temperature; low enough for detecting a climate change signal within about a decade. Larger structural uncertainty above about 25 km and at high latitudes is attributable to differences in the processing schemes, which undergo continuous improvements. Though current use of RO for reliable climate trend assessment is bound to 50° S to 50° N, our results show that quality, consistency, and reproducibility are favorable in the UTLS for the establishment of a climate benchmark record

Estimating Atmospheric Boundary Layer Depth Using COSMIC Radio Occultation Data

This study presents an algorithm for estimating atmospheric boundary layer (ABL) depth from Global Positioning System (GPS) radio occultation (RO) data. The algorithm is applied to the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) RO data and validated using high-resolution radiosonde data from the island of St. Helena (16.0 degrees S, 5.7 degrees W), tropical (30 degrees S-30 degrees N) radiosondes collocated with RO, and European Centre for Medium-Range Weather Forecasts (ECMWF) high-resolution global analyses. Spatial and temporal variations of the ABL depth obtained from COSMIC RO data for a 1-yr period over tropical and subtropical oceans are analyzed. The results demonstrate the capability of RO data to resolve geographical and seasonal variations of ABL height. The spatial patterns of the variations are consistent with those derived from ECMWF global analysis. However, the ABL heights derived from ECMWF global analysis, on average, are negatively biased against those estimated from COSMIC GPS RO data. These results indicate that GPS RO data can provide useful information on ABL height, which is an important parameter for weather and climate studies

Assimilation of GPS radio occultation data for Numerical Weather

With the availability of approximately 4,000 radio occultation soundings per day within three hours of observation, COSMIC has the potential to contribute significantly to global and regional weather analysis and prediction. However, the basic radio occultation measurements (phase delays) are very different from traditional meteorological measurements (i.e., temperature, water vapor), and to effectively assimilate them into weather prediction models is a challenging task. Over the past five years, considerable progress has been made in the development of an effective strategy for the assimilation of GPS radio occultation data. In this paper, we (1) review the measurement and data reduction procedures, (2) discuss the error characteristics of the GPS radio occultation data, (3) discuss the various strategies for data assimilation, (4) review results from recent data assimilation research, and (5) provide suggestions for future research. Results from recent studies have led to the conclusion that the best strategy to assimilate GPS radio occultation data is a mixture of bending angles below 10 km and refractivity above 10 km using a variational approach. The assimilation of GPS radio occultation data is likely to have a significant positive impact on global and regional weather prediction through improved definition of wate

Ionospheric correction of GPS radio occultation data in the troposphere

For inversions of the GPS radio occultation (RO) data in the neutral atmosphere, this study investigates an optimal transition height for replacing the standard ionospheric correction using the linear combination of the L1 and L2 bending angles with the correction of the L1 bending angle by the L1–L2 bending angle extrapolated from above. The optimal transition height depends on the RO mission (i.e., the receiver and firmware) and is different between rising and setting occultations and between L2P and L2C GPS signals. This height is within the range of approximately 10–20 km. One fixed transition height, which can be used for the processing of currently available GPS RO data, can be set to 20 km. Analysis of the L1CA and the L2C bending angles shows that in some occultations the errors of standard ionospheric correction substantially increase around the strong inversion layers (such as the top of the boundary layer). This error increase is modeled and explained by the horizontal inhomogeneity of the ionosphere

Error analysis of Abel retrieved electron density profiles from radio occultation measurements

This letter reports for the first time the simulated error distribution of radio occultation (RO) electron density profiles (EDPs) from the Abel inversion in a systematic way. Occultation events observed by the COSMIC satellites are simulated during the spring equinox of 2008 by calculating the integrated total electron content (TEC) along the COSMIC occultation paths with the "true" electron density from an empirical model. The retrieval errors are computed by comparing the retrieved EDPs with the "true" EDPs. The results show that the retrieved NmF2 and hmF2 are generally in good agreement with the true values, but the reliability of the retrieved electron density degrades in low latitude regions and at low altitudes. Specifically, the Abel retrieval method overestimates electron density to the north and south of the crests of the equatorial ionization anomaly (EIA), and introduces artificial plasma caves underneath the EIA crests. At lower altitudes (E- and F1-regions), it results in three pseudo peaks in daytime electron densities along the magnetic latitude and a pseudo trough in nighttime equatorial electron densities