Christian Rocken

GPS meteorology: Remote sensing of atmospheric water vapor using the global positioning system

We present a new approach to remote sensing of water vapor based on the global positioning system (GPS). Geodesists and geophysicists have devised methods for estimating the extent to which signals propagating from GPS satellites to ground-based GPS receivers are delayed by atmospheric water vapor. This delay is parameterized in terms of a time-varying zenith wet delay (ZWD) which is retrieved by stochastic filtering of the GPS data. Given surface temperature and pressure readings at the GPS receiver, the retrieved ZWD can be transformed with very little additional uncertainty into an estimate of the integrated water vapor (IWV) overlying that receiver. Networks of continuously operating GPS receivers are being constructed by geodesists, geophysicists, government and military agencies, and others in order to implement a wide range of positioning capabilities. These emerging GPS networks offer the possibility of observing the horizontal distribution of IWV or, equivalently, precipitable water with unprecedented coverage and a temporal resolution of the order of 10 min. These measurements could be utilized in operational weather forecasting and in fundamental research into atmospheric storm systems, the hydrologic cycle, atmospheric chemistry, and global climate change. Specially designed, dense GPS networks could be used to sense the vertical distribution of water vapor in their immediate vicinity. Data from ground-based GPS networks could be analyzed in concert with observations of GPS satellite occultations by GPS receivers in low Earth orbit to characterize the atmosphere at planetary scale.

GPS meteorology: mapping zenith wet delays onto precipitable water

Abstract Emerging networks of Global Positioning System (GPS) receivers can be used in the remote sensing of atmospheric water vapor. The time-varying zenith wet delay observed at each GPS receiver in a network can be transformed into an estimate of the precipitable water overlying that receiver. This transformation is achieved by multiplying the zenith wet delay by a factor whose magnitude is a function of certain constants related to the refractivity of moist air and of the weighted mean temperature of the atmosphere. The mean temperature varies in space and time and must be estimated a priori in order to transform an observed zenith wet delay into an estimate of precipitable water. We show that the relative error introduced during this transformation closely approximates the relative error in the predicted mean temperature. Numerical weather models can be used to predict the mean temperature with an rms relative error of less than 1%.

THE COSMIC/FORMOSAT-3 MISSION : Early Results

The radio occultation (RO) technique, which makes use of radio signals transmitted by the global positioning system (GPS) satellites, has emerged as a powerful and relatively inexpensive approach for sounding the global atmosphere with high precision, accuracy, and vertical resolution in all weather and over both land and ocean. On 15 April 2006, the joint Taiwan-U.S. Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC)/Formosa Satellite Mission 3 (COSMIC/FORMOSAT-3, hereafter COSMIC) mission, a constellation of six microsatellites, was launched into a 512-km orbit. After launch the satellites were gradually deployed to their final orbits at 800 km, a process that took about 17 months. During the early weeks of the deployment, the satellites were spaced closely, offering a unique opportunity to verify the high precision of RO measurements. As of September 2007, COSMIC is providing about 2000 RO soundings per day to support the research and operational communities. COSMIC RO dat...

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...

GPS Meteorology: Direct Estimation of the Absolute Value of Precipitable Water

Abstract A simple approach to estimating vertically integrated atmospheric water vapor, or precipitable water, from Global Positioning System (GPS) radio signals collected by a regional network of ground-based geodetic GPS receiver is illustrated and validated. Standard space geodetic methods are used to estimate the zenith delay caused by the neutral atmosphere, and surface pressure measurements are used to compute the hydrostatic (or “dry”) component of this delay. The zenith hydrostatic delay is subtracted from the zenith neutral delay to determine the zenith wet delay, which is then transformed into an estimate of precipitable water. By incorporating a few remote global tracking stations (and thus long baselines) into the geodetic analysis of a regional GPS network, it is possible to resolve the absolute (not merely the relative) value of the zenith neutral delay at each station in the augmented network. This approach eliminates any need for external comparisons with water vapor radiometer observation...

A Global Morphology of Gravity Wave Activity in the Stratosphere Revealed by the GPS Occultation Data (GPS/MET)

Using temperature profiles obtained by the GPS/MET (GPS Meteorology) experiment from April 1995 to February 1997, we have extracted mesoscale temperature perturbations with vertical wavelengths ranging from 2 to 10 km and background Brunt-Vaisala frequency squared, N2. For each occultation event, we can evaluate a potential energy Ep which is assumed to be caused by atmospheric gravity waves. The monthly mean values of Ep at 15–20 km around Japan showed an annual variation with an enhancement in winter, which is consistent with the climatological behavior of the kinetic energy of gravity waves observed with the MU (middle and upper atmosphere) radar (34.9°N, 136.0°E) from 1985 to 1989. We have then derived the global distribution of Ep at 20–30 km during Northern Hemisphere winter (from November to February). Our analysis shows that the largest Ep values are generally centered around the equator between 25°N and 25°S with considerable longitude variations. Longitudinal variations of Ep at 20–30 km in a latitude range of 30°–60°N are also analyzed, resulting in larger Ep values over the continents than over the Pacific Ocean. Using GPS/MET data without antispoofing, latitudinal variations of Ep are determined in 15–45 km. Although large Ep values are concentrated near the equator at 20–30 km, Ep tends to become larger at midlatitudes at 30–40 km and higher-altitude regions. At midlatitudes, Ep is found to be larger in winter months in both hemispheres. Height variations of Ep indicate a decrease at 25–30 km and a monotonic increase above 30 km.

Analysis and validation of GPS/MET radio occultation data in the ionosphere

Global Positioning System (GPS) radio occultation signals received by a low Earth orbit (LEO) satellite provide information about the global distribution of electron den- sity in the ionosphere. We examine two radio occultation inversion algorithms. The first algo- rithm utilizes the Abel integral transform, which assumes spherical symmetry of the electron density field. We test this algorithm with two approaches: through the computation of bend- ing angles and through the computation of total electron content (TEC) assuming straight line propagation. We demonstrate that for GPS frequencies and for observations in LEO, the as- sumption of straight-line propagation (neglecting bending) introduces small errors when monitoring the F2 layer. The second algorithm, which also assumes straight-line propagation, is a three-dimensional (3-D) inversion constrained with the horizontal structure of a priori electron density fields. As a priori fields we use tomographic solutions and the parameterized real-time ionospheric specification model (PRISM) when adjusted with ionosonde data or ground-based GPS vertical TEC maps. For both algorithms we calibrate the occultation data by utilizing observations from the part of the LEO that is closer to the GPS satellite. For in- versions we use dual-frequency observational data (the difference of L1 and L2 phase ob- servables) which cancel orbit errors (without applying precise orbit determination) and clock errors (without requiring synchronous ground data) and thus may allow inversions to be computed close to real time in the future. The Abel and 3-D constrained algorithms are vali- dated by statistically comparing 4 days of inversions with critical frequency (foF2) data from a network of 45 ionosonde stations and with vertical TEC data from the global network of GPS ground receivers. Globally, the Abel inversion approach agrees with the foF2 correlative data at the 13% rms level, with a negligible mean difference. All tested 3-D constrained in- version approaches possess a statistically significant mean difference when compared with the ionosonde data. The vertical TEC correlative comparisons for both the Abel and 3-D con- strained inversions are significantly biased (-30%) by the electrons above the 735-km LEO altitude.

Sensing atmospheric water vapor with the global positioning system

Global Positioning System (GPS) receivers, water vapor radiometers (WVRs), and surface meteorological equip- ment were operated at both ends of a 50-kin baseline in Colorado to measure the precipitable water vapor (PWV) and wet delay in the line-of-sight to GPS satellites. Using high pre- cision orbits, WVR-measured and GPS-inferred PWV differences between the two sites usually agreed to better than 1 min. Using less precise on-line broadcast orbits increased the discrepancy by 30%. Data simulations show that GPS mea- surements can provide ram-level separate PWV estimates for the two sites, as opposed to just their difference, if baselines exceed 500 km and the highest accuracy GPS orbits are used.

GPS/STORM—GPS Sensing of Atmospheric Water Vapor for Meteorology

Abstract Atmospheric water vapor was measured with six Global Positioning System (GPS) receivers for 1 month at sites in Colorado, Kansas, and Oklahoma. During the time of the experiment from 7 May to 2 June 1993, the area experienced severe weather. The experiment, called “GPS/STORM,” used GPS signals to sense water vapor and tested the accuracy of the method for meteorological applications. Zenith wet delay and precipitable water (PW) were estimated, relative to Platteville, Colorado, every 30 min at five sites. At three of these five sites the authors compared GPS estimates of PW to water vapor radiometer (WVR) measurements. GPS and WVR estimates agree to 1–2 mm rms. For GPS/STORM site spacing of 500–900 km, high-accuracy GPS satellite orbits are required to estimate 1–2-mm-level PW. Broadcast orbits do not have sufficient accuracy. It is possible, however, to estimate orbit improvements simultaneously with PW. Therefore, it is feasible that future meteorological GPS networks provide near-real-time hig...