Michael E. Gorbunov

Analysis and validation of GPS/MET data in the neutral atmosphere

The Global Positioning System/Meteorology ( GPS/MET) Program was established in 1993 by the University Corporation for Atmospheric Research ( UCAR) to demonstrate active limb sounding of the Earths atmosphere using the radio occultation technique. The demonstration system observes occulted GPS satellite signals received by a low Earth orbiting ( LEO) satellite, MicroLab-1, launched April 3,1995. The system can profile ionospheric electron density and neutral atmospheric properties. Neutral atmospheric refractivity, density, pressure, and temperature are derived at altitudes where the amount of water vapor is low. At lower altitudes, vertical profiles of density, pressure, and water vapor pressure can be derived from the GPS/MET refractivity profiles if temperature data from an independent source are available. This paper describes the GPS/MET data analysis procedures and validates GPS/MET data with statistics and illustrative case studies. We compare more than 1200 GPS/MET neutral atmosphere soundings to correlative data from operational global weather analyses, radiosondes, and the GOES, TOVS, UARS/MLS and HALOE orbiting atmospheric sensors. Even though many GPS/MET soundings currently fail to penetrate the lowest 5 km of the troposphere in the presence of significant water vapor, our results demonstrate 1°C mean temperature agreement with the best correlative data sets between 1 and 40 km. This and the fact that GPS/MET observations are all-weather and self-calibrating suggests that radio occultation technology has the potential to make a strong contribution to a global observing system supporting weather prediction and weather and climate research.

GPS Sounding of the Atmosphere from Low Earth Orbit: Preliminary Results

Abstract This paper provides an overview of the methodology of and describes preliminary results from an experiment called GPS/MET (Global Positioning System/Meteorology), in which temperature soundings are obtained from a low Earth-orbiting satellite using the radio occultation technique. Launched into a circular orbit of about 750-km altitude and 70° inclination on 3 April 1995, a small research satellite, MicroLab 1, carried a laptop-sized radio receiver. Each time this receiver rises and sets relative to the 24 operational GPS satellites, the GPS radio waves transect successive layers of the atmosphere and are bent (refracted) by the atmosphere before they reach the receiver, causing a delay in the dual-frequency carrier phase observations sensed by the receiver. During this occultation, GPS limb sounding measurements are obtained from which vertical profiles of atmospheric refractivity can be computed. The refractivity is a function of pressure, temperature, and water vapor and thus provides informat...

COSMIC Radio Occultation Processing: Cross-Center Comparison and Validation

AbstractA radio occultation data processing system (OCC) was developed for numerical weather prediction and climate benchmarking. The data processing algorithms use the well-established Fourier integral operator–based methods, which ensure a high accuracy of retrievals. The system as a whole, or in its parts, is currently used at the Global Navigation Satellite System Receiver for Atmospheric Sounding (GRAS) Satellite Application Facility at the Danish Meteorological Institute, German Weather Service, and Wegener Center for Climate and Global Change. A statistical comparison of the inversions of the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) data by the system herein, University Corporation for Atmospheric Research (UCAR) data products, and ECMWF analyses is presented. Forty days of 2007 and 2008 were processed (from 5 days in the middle of each season) for the comparison of OCC and ECMWF, and 20 days of April 2009 were processed for the comparison of OCC, UCAR, and EC...

Spatiotemporal structure of a laser beam over 144 km in a Canary Islands experiment

We analyzed the observations of scintillations in a laser beam (532 nm, ~200 mW power) traveling along a 144 km path at an altitude of 2.2-2.4 km above sea level, just above the atmospheric boundary layer, between the islands of La Palma and Tenerife. The observations were performed during nighttime on 18 July 2011, by means of a telescope with an aperture diameter of 1 m. Strong scintillations were observed. The estimates of spatial spectra and correlation functions indicated that the observed intensity fields possess, statistically, a locally isotropic structure, which agrees with the idea of a locally isotropic turbulence. The estimates of spatial autospectra and autocorrelation functions of the intensity field indicated that the characteristic scale of the internal structure of the observed clusters is 6.5-8 mm, while the characteristic size of the clusters is 4-5 cm. The major contribution to the observed scintillations comes from the inhomogeneities of the intensity field with scales from 1-2 cm up to 10-12 cm. The analysis of the cross-spectra indicated that the hypothesis of frozen turbulence introduced by Taylor can be used for the description of spatiotemporal structure of intensity fluctuations of laser beams traveling through long paths in the atmosphere.

Uncertainty propagation through wave optics retrieval of bending angles from GPS radio occultation: Theory and simulation results

The wave optical technique for bending angle retrieval in processing radio occultation observations is nowadays widely used by different data processing and assimilation groups and centers. This technique uses Fourier Integral Operators that map the observed records of the amplitude and phase into the impact parameter representation, which allows for the retrieval of bending angle as a function of impact parameter. We investigate the propagation of uncertainty in the observed amplitude and excess phase to the retrieved bending angle. We construct a simple linear approximation, where the excess phase uncertainty is mapped into the bending angle uncertainty. This results in a simple analytical expression for the final uncertainty. To verify our approximation, we perform numerical Monte Carlo simulations for three example occultation events (tropical, middle, and polar latitude profiles from an atmospheric analysis). We demonstrate that our approximation basically gives good results in all cases over the entire troposphere. Exception is the narrow area near the top of the sharp boundary layer, especially in tropics, where, due to nonlinear effects, a significant systematic error arises accompanied by increased uncertainty.

Fluctuations of refractivity as a systematic error source in radio occultations

The fact that fluctuations of refractivity may result in a systematic negative shift of the phase of waves propagating in a random medium is known for a long time. Tatarskii was the first to reveal it, and von Eshleman put this into the context of the radio occultation sounding of planetary atmospheres. In this paper, we show that this effect may also be one of the causes of the negative bias of refractivity retrieved for radio occultation observations of the Earths atmosphere. We perform theoretical estimates of this effect based on the Rytov approximation. These estimates, however, do not consider the regular refraction, which may significantly change the magnitude of this effect. We perform numerical simulations of radio occultations, based on the Kolmogorov-von Karman isotropic spectrum of refractivity fluctuations, with the internal and external scales and magnitude tuned so as to reproduce the realistic level of the variance of retrieved refractivity and the amplitude fluctuations of the modeled signals. The model of the regular atmosphere is based on analyses of the European Centre for Medium-Range Weather Forecasts. We show that it is possible to set up a vertical profile of the structural constant of the fluctuation spectrum such that it will result in a systematic shift and variances of the retrieved refractivity consistent with those observed for COSMIC measurements.

Canonical Transform Methods for Analysis of Radio Occultations

We discuss the application of canonical transform (CT) methods for analyzing radio occultations with multipath behavior. In the present work we discuss a mapping to the representation of approximate impact parameter. This method generalizes the full spectrum inversion (FSI) method for the case of noncircular orbits. The method is based on mapping the field by a Fourier integral operator that maps directly from the measured time-dependent field to the impact parameter representation without first doing a back propagation. Furthermore, our method allows for simple, asymptotic direct modelling of wave propagation. We perform processing of simulated radio occultations and show that our method provides the same accuracy and resolution as the standard CT method.

Assimilation of CHAMP and GRACE-A Radio Occultation Data in the GME Global Meteorological Model of the German Weather Service

The assimilation of GPS radio occultations within the three-dimensional variational data assimilation system of the German Weather Service requires GPS radio occultation bending angle forward operators. To optimize the forward operator setup, different one- and three-dimensional bending angle forward operators are evaluated. The innovation statistics for radio occultation data from the CHAMP, GRACE-A and FORMOSAT-3/COSMIC satellites are compared with estimates based on the background and observation errors specified in the assimilation scheme. Numerical experiments are performed to assess the impact of assimilated radio occultation data on the weather forecast scores.