A ray-tracing operator and its adjoint for the use of GPS/MET refraction angle measurements

The development of small, high-performance instruments to receive Global Positioning System (GPS) signals has created an opportunity for active remote sounding of the Earth's atmosphere by radio occultation techniques. A prototype demonstration of this capability has been provided by the GPS Meteorology (GPS/MET) experiment. Although it was shown that high vertical resolution profiles of atmospheric refractivity, temperature, and geopotential height of constant pressure levels can be derived from the GPS measurements, with high accuracy under many circumstances, many issues remain. These include the existence of multipath propagation, the ambiguity between water vapor and temperature in moist regions of the atmosphere, and the difficulty in retrieving an accurate refractivity profile from the GPS refraction angle measurements over regions where the horizontal gradient of the refractivity is large. The aim of this paper is to begin the development of a methodology for incorporating the GPS "raw" measurements (refraction angles) directly into numerical weather analysis and/or prediction systems in order to alleviate the above mentioned problems. First, a ray-tracing observation operator that links the atmospheric state to the GPS refraction angle measurements is developed, the physics and numerics involved are described, and the simulated refraction angles, based on the NOAA National Centers for Environmental Prediction (NCEP) global analysis, are compared with the observed GPS/MET refraction angle measurements. Second, the tangent linear and adjoint of the ray-tracing operator are developed. These three operators are required for the direct use of GPS refraction angle measurements in a variational data analysis system. A single observation experiment reveals that the direct use of GPS refraction angles in a variational analysis causes changes in the temperature and specific humidity fields that are not limited to the occultation location but in an elongated band of 6300 km in its occultation plane. On s-levels, changes from the use of one GPS occultation occur in an area of about 600 km 3 600 km large which is centered around the ray perigee point. Finally, the advantages and disadvantages are discussed for the use of the GPS refractivities versus refraction angles. Errors made by using local estimates of refractivity are also assessed.

[1]  E. Kursinski.  Initial results of radio occultation observations of Earth's atmophere: using the Global Positioning , 1996 .

[2]  Lennart Bengtsson,et al.  Advanced algorithms of inversion of GPS/MET satellite data and their application to reconstruction of temperature and humidity , 1996 .

[3]  X. Zou,et al.  Introduction to Adjoint Techniques and the MM5 Adjoint Modeling System. , 1997 .

[4]  J. Barnett,et al.  COSPAR International reference atmosphere grand mean , 1990 .

[5]  Joseph Sela,et al.  Spectral Modeling at the National Meteorological Center , 1980 .

[6]  V. V. Vorob’ev,et al.  Estimation of the accuracy of the atmospheric refractive index recovery from Doppler shift measurements at frequencies used in the NAVSTAR system , 1994 .

[7]  X. Zou,et al.  Analysis and validation of GPS/MET data in the neutral atmosphere , 1997 .

[8]  H. H. Rachford,et al.  The Numerical Solution of Parabolic and Elliptic Differential Equations , 1955 .

[9]  N. N. Yanenko,et al.  The Method of Fractional Steps , 1971 .

[10]  N. McFarlane,et al.  The mineral dust aerosol cycle during the Last Glacial Maximum , 1999 .

[11]  Steven Businger,et al.  GPS Sounding of the Atmosphere from Low Earth Orbit: Preliminary Results , 1996 .

[12]  M. E. Gorbunov,et al.  Microlab‐1 experiment: Multipath effects in the lower troposphere , 1998 .

[13]  P. Courtier,et al.  The ECMWF implementation of three‐dimensional variational assimilation (3D‐Var). I: Formulation , 1998 .

[14]  Christian Rocken,et al.  A GPS/MET Sounding through an Intense Upper-Level Front. , 1998 .

[15]  Stephen S. Leroy,et al.  Measurement of geopotential heights by GPS radio occultation , 1997 .

[16]  Ying-Hwa Kuo,et al.  Assimilation of Atmospheric Radio Refractivity Using a Nonhydrostatic Adjoint Model , 1995 .

[17]  Michael E. Gorbunov,et al.  Remote sensing of refractivity from space for global observations of atmospheric parameters , 1993 .

[18]  S. K. Runcorn,et al.  International Dictionary of Geophysics , 1967 .

[19]  J. Schofield,et al.  Observing Earth's atmosphere with radio occultation measurements using the Global Positioning System , 1997 .

[20]  Larry J. Romans,et al.  Observing tropospheric water vapor by radio occultation using the Global Positioning System , 1995 .

[21]  Y. Kravtsov,et al.  Geometrical optics of inhomogeneous media , 2019, Geometrical Optics of Weakly Anisotropic Media.

[22]  R. J. List.  Smithsonian Meteorological Tables , 2018, Nature.

[23]  A. Kliore,et al.  The neutral atmosphere of Venus as studied with the Mariner V radio occultation experiments , 1971 .

[24]  W. G. Melbourne,et al.  Initial Results of Radio Occultation Observations of Earth's Atmosphere Using the Global Positioning System , 1996, Science.