A simplified approach for generating GNSS radio occultation refractivity climatologies

Abstract. The possibility of simplifying the retrieval scheme required to produce GNSS radio occultation refractivity climatologies is investigated. In a new, simplified retrieval approach, the main statistical analysis is performed in bending angle space and an estimate of the average bending angle profile is then propagated through an Abel transform. The average is composed of means and medians of ionospheric corrected bending angles up to 80 km. Above that, the observed profile is exponentially extrapolated to infinity using a fixed a priori scale height. The new approach circumvents the need to introduce a "statistical optimisation" processing step in which individual bending angle profiles are merged with a priori data, often taken from a climatology. This processing step can be complex, difficult to interpret, and is generally recognised as a potential source of structural uncertainty. The new scheme is compared with the more conventional approach of averaging individual refractivity profiles, produced with the implementation of statistical optimisation used in the EUMETSAT Radio Occultation Meteorology Satellite Application Facility (ROM SAF) operational processing. It is shown that the two GNSS radio occultation climatologies agree to within 0.1% from 5 km up to 35–40 km, for the three months January, February, and March 2011. During this time period, the new approach also produces slightly better agreement with ECMWF analyses between 40–50 km, which is encouraging. The possible limitations of the new approach caused by mean residual ionospheric errors and low observation numbers are discussed briefly, and areas for future work are suggested.

[1]  Gottfried Kirchengast,et al.  Advancements of Global Navigation Satellite System radio occultation retrieval in the upper stratosphere for optimal climate monitoring utility , 2004 .

[2]  S. B. Healy,et al.  Monitoring twenty‐first century climate using GPS radio occultation bending angles , 2008 .

[3]  Stig Syndergaard,et al.  Processing and validation of refractivity from GRAS radio occultation data , 2011 .

[4]  John R. Christy,et al.  UNCERTAINTIES IN CLIMATE TRENDS Lessons from Upper-Air Temperature Records , 2005 .

[5]  M. E. Gorbunov,et al.  Ionospheric correction and statistical optimization of radio occultation data , 2002 .

[6]  M. Lohmann,et al.  Application of dynamical error estimation for statistical optimization of radio occultation bending angles , 2005 .

[7]  Ying-Hwa Kuo,et al.  Estimating the uncertainty of using GPS radio occultation data for climate monitoring: Intercomparison of CHAMP refractivity climate records from 2002 to 2006 from different data centers , 2009 .

[8]  Ying-Hwa Kuo,et al.  Quantification of structural uncertainty in climate data records from GPS radio occultation , 2012 .

[9]  William H. Press,et al.  Numerical Recipes in Fortran 77 , 1992 .

[10]  J. Wickert,et al.  Reproducibility of GPS radio occultation data for climate monitoring: Profile‐to‐profile inter‐comparison of CHAMP climate records 2002 to 2008 from six data centers , 2012 .

[11]  A. Engeln.  A first test of climate monitoring with radio occultation instruments: Comparing two processing centers , 2006 .

[12]  A. Hedin.  Extension of the MSIS Thermosphere Model into the middle and lower atmosphere , 1991 .

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

[14]  Anthony J. Mannucci,et al.  Improving GPS Radio occultation stratospheric refractivity retrievals for climate benchmarking , 2012 .