scholarly journals Wave-optics uncertainty propagation and regression-based bias model in GNSS radio occultation bending angle retrievals

2018 ◽  
Vol 11 (1) ◽  
pp. 111-125 ◽  
Author(s):  
Michael E. Gorbunov ◽  
Gottfried Kirchengast
Keyword(s):  
Regression Analysis ◽  
Systematic Uncertainty ◽  
Radio Occultation ◽  
Uncertainty Propagation ◽  
Lower Troposphere ◽  
Wave Optics ◽  
Bending Angle ◽  
Excess Phase ◽  
Excess Phases ◽  
Random Uncertainties

Abstract. A new reference occultation processing system (rOPS) will include a Global Navigation Satellite System (GNSS) radio occultation (RO) retrieval chain with integrated uncertainty propagation. In this paper, we focus on wave-optics bending angle (BA) retrieval in the lower troposphere and introduce (1) an empirically estimated boundary layer bias (BLB) model then employed to reduce the systematic uncertainty of excess phases and bending angles in about the lowest 2 km of the troposphere and (2) the estimation of (residual) systematic uncertainties and their propagation together with random uncertainties from excess phase to bending angle profiles. Our BLB model describes the estimated bias of the excess phase transferred from the estimated bias of the bending angle, for which the model is built, informed by analyzing refractivity fluctuation statistics shown to induce such biases. The model is derived from regression analysis using a large ensemble of Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) RO observations and concurrent European Centre for Medium-Range Weather Forecasts (ECMWF) analysis fields. It is formulated in terms of predictors and adaptive functions (powers and cross products of predictors), where we use six main predictors derived from observations: impact altitude, latitude, bending angle and its standard deviation, canonical transform (CT) amplitude, and its fluctuation index. Based on an ensemble of test days, independent of the days of data used for the regression analysis to establish the BLB model, we find the model very effective for bias reduction and capable of reducing bending angle and corresponding refractivity biases by about a factor of 5. The estimated residual systematic uncertainty, after the BLB profile subtraction, is lower bounded by the uncertainty from the (indirect) use of ECMWF analysis fields but is significantly lower than the systematic uncertainty without BLB correction. The systematic and random uncertainties are propagated from excess phase to bending angle profiles, using a perturbation approach and the wave-optical method recently introduced by Gorbunov and Kirchengast (2015), starting with estimated excess phase uncertainties. The results are encouraging and this uncertainty propagation approach combined with BLB correction enables a robust reduction and quantification of the uncertainties of excess phases and bending angles in the lower troposphere.

scholarly journals Wave-optics uncertainty propagation and regression-based bias model in GNSS radio occultation bending angle retrievals

2017 ◽  
Author(s):  
Michael Gorbunov ◽  
Gottfried Kirchengast
Keyword(s):  
Regression Analysis ◽  
Systematic Uncertainty ◽  
Radio Occultation ◽  
Uncertainty Propagation ◽  
Lower Troposphere ◽  
Wave Optics ◽  
Bending Angle ◽  
Excess Phase ◽  
Excess Phases ◽  
Random Uncertainties

Abstract. A new reference occultation processing system (rOPS) will include a Global Navigation Satellite System (GNSS) radio occultation (RO) retrieval chain with integrated uncertainty propagation. In this paper, we focus on wave-optics bending angle retrieval in the lower troposphere and introduce 1. an empirically estimated boundary layer bias (BLB) model then employed to reduce the systematic uncertainty of excess phases and bending angles in the lowest about two kilometers of the troposphere, and 2. the estimation of (residual) systematic uncertainties and their propagation together with random uncertainties from excess phase to bending angle profiles. Our BLB model describes the estimated bias of the excess phase transferred from the estimated bias of the bending angle, for which the model is built, informed by analyzing refractivity fluctuation statistics shown to induce such biases. The model is derived from regression analysis using a large ensemble of Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) RO observations and concurrent European Centre for Medium-Range Weather Forecasts (ECMWF) analysis fields. It is formulated in terms predictors and adaptive functions (powers and cross-products of predictors), where we use six main predictors derived from observations: impact altitude, latitude, bending angle and its standard deviation, canonical transform amplitude and its fluctuation index. Based on an ensemble of test days, independent of the days of data used for the regression analysis to establish the BLB model, we find the model very effective for bias reduction, capable of reducing bending angle and corresponding refractivity biases by about a factor of five. The estimated residual systematic uncertainty, after the BLB profile subtraction, is lower-bounded by the uncertainty from (indirect) use of ECMWF analysis fields but is significantly lower than the systematic uncertainty without BLB correction. The systematic and random uncertainties are propagated from excess phase to bending angle profiles, using a perturbation approach and the wave-optical method recently introduced by Gorbunov and Kirchengast (2015), starting with estimated excess phase uncertainties. The results are encouraging that this uncertainty propagation approach combined with BLB correction enables a robust reduction and quantification of the uncertainties of excess phases and bending angles in the lower troposphere.

scholarly journals Integrating uncertainty propagation in GNSS radio occultation retrieval: from excess phase to atmospheric bending angle profiles

2017 ◽  
Author(s):  
Jakob Schwarz ◽  
Gottfried Kirchengast ◽  
Marc Schwaerz
Keyword(s):  
Monte Carlo ◽  
Radio Occultation ◽  
Uncertainty Propagation ◽  
Processing System ◽  
Background Information ◽  
Climate Variables ◽  
Free Atmosphere ◽  
Bending Angle ◽  
Covariance Propagation ◽  
Excess Phase

Abstract. Global Navigation Satellite System (GNSS) radio occultation (RO) observations are highly accurate, long-term stable data sets, and are globally available as a continuous record since 2001. Essential climate variables for the thermodynamic state of the free atmosphere, such as pressure, temperature and tropospheric water vapor profiles (involving background information), can be derived from these records, which therefore have the potential to serve as climate benchmark data. However, to exploit this potential, atmospheric profile retrievals need to be very accurate and the remaining uncertainties quantified and traced throughout the retrieval chain from raw observations to essential climate variables. The new Reference Occultation Processing System (rOPS) at the Wegener Center aims to deliver such an accurate RO retrieval chain with integrated uncertainty propagation. Here we introduce and demonstrate the algorithms implemented in the rOPS for uncertainty propagation from excess phase to atmospheric bending angle profiles, for estimated systematic and random uncertainties, including vertical error correlations and resolution estimates. We estimated systematic uncertainty profiles with the same operators as used for the basic state profiles retrieval. The random uncertainty is traced through covariance propagation and validated using Monte-Carlo ensemble methods. The algorithm performance is demonstrated using test-day ensembles of simulated data as well as real RO event data from the satellite missions CHAMP and COSMIC. The results of the Monte-Carlo validation show that our covariance propagation delivers correct uncertainty quantification from excess phase to bending angle profiles. The results from the real RO event ensembles demonstrate that the new uncertainty estimation chain performs robustly. Together with the other parts of the rOPS processing chain this part is thus ready to provide integrated uncertainty propagation through the whole RO retrieval chain for the benefit of climate monitoring and other applications

scholarly journals Integrating uncertainty propagation in GNSS radio occultation retrieval: from excess phase to atmospheric bending angle profiles

2018 ◽  
Vol 11 (5) ◽  
pp. 2601-2631 ◽  
Author(s):  
Jakob Schwarz ◽  
Gottfried Kirchengast ◽  
Marc Schwaerz
Keyword(s):  
Monte Carlo ◽  
Radio Occultation ◽  
Uncertainty Propagation ◽  
Processing System ◽  
Background Information ◽  
Climate Variables ◽  
Free Atmosphere ◽  
Bending Angle ◽  
Covariance Propagation ◽  
Excess Phase

Abstract. Global Navigation Satellite System (GNSS) radio occultation (RO) observations are highly accurate, long-term stable data sets and are globally available as a continuous record from 2001. Essential climate variables for the thermodynamic state of the free atmosphere – such as pressure, temperature, and tropospheric water vapor profiles (involving background information) – can be derived from these records, which therefore have the potential to serve as climate benchmark data. However, to exploit this potential, atmospheric profile retrievals need to be very accurate and the remaining uncertainties quantified and traced throughout the retrieval chain from raw observations to essential climate variables. The new Reference Occultation Processing System (rOPS) at the Wegener Center aims to deliver such an accurate RO retrieval chain with integrated uncertainty propagation. Here we introduce and demonstrate the algorithms implemented in the rOPS for uncertainty propagation from excess phase to atmospheric bending angle profiles, for estimated systematic and random uncertainties, including vertical error correlations and resolution estimates. We estimated systematic uncertainty profiles with the same operators as used for the basic state profiles retrieval. The random uncertainty is traced through covariance propagation and validated using Monte Carlo ensemble methods. The algorithm performance is demonstrated using test day ensembles of simulated data as well as real RO event data from the satellite missions CHAllenging Minisatellite Payload (CHAMP); Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC); and Meteorological Operational Satellite A (MetOp). The results of the Monte Carlo validation show that our covariance propagation delivers correct uncertainty quantification from excess phase to bending angle profiles. The results from the real RO event ensembles demonstrate that the new uncertainty estimation chain performs robustly. Together with the other parts of the rOPS processing chain this part is thus ready to provide integrated uncertainty propagation through the whole RO retrieval chain for the benefit of climate monitoring and other applications.

scholarly journals Global 3D Features of Error Variances of GPS Radio Occultation and Radiosonde Observations

2020 ◽  
Vol 13 (1) ◽  
pp. 1
Author(s):  
Xu Xu ◽  
Xiaolei Zou
Keyword(s):  
Radio Occultation ◽  
Weather Prediction ◽  
Lower Troposphere ◽  
Error Variance ◽  
Specific Humidity ◽  
Temperature Error ◽  
Observation Error ◽  
Bending Angle ◽  
Gps Radio Occultation ◽  
Low Latitudes

Global Positioning System (GPS) radio occultation (RO) and radiosonde (RS) observations are two major types of observations assimilated in numerical weather prediction (NWP) systems. Observation error variances are required input that determines the weightings given to observations in data assimilation. This study estimates the error variances of global GPS RO refractivity and bending angle and RS temperature and humidity observations at 521 selected RS stations using the three-cornered hat method with additional ERA-Interim reanalysis and Global Forecast System forecast data available from 1 January 2016 to 31 August 2019. The global distributions, of both RO and RS observation error variances, are analyzed in terms of vertical and latitudinal variations. Error variances of RO refractivity and bending angle and RS specific humidity in the lower troposphere, such as at 850 hPa (3.5 km impact height for the bending angle), all increase with decreasing latitude. The error variances of RO refractivity and bending angle and RS specific humidity can reach about 30 N-unit2, 3 × 10−6 rad2, and 2 (g kg−1)2, respectively. There is also a good symmetry of the error variances of both RO refractivity and bending angle with respect to the equator between the Northern and Southern Hemispheres at all vertical levels. In this study, we provide the mean error variances of refractivity and bending angle in every 5°-latitude band between the equator and 60°N, as well as every interval of 10 hPa pressure or 0.2 km impact height. The RS temperature error variance distribution differs from those of refractivity, bending angle, and humidity, which, at low latitudes, are smaller (less than 1 K2) than those in the midlatitudes (more than 3 K2). In the midlatitudes, the RS temperature error variances in North America are larger than those in East Asia and Europe, which may arise from different radiosonde types among the above three regions.

scholarly journals Estimation and evaluation of COSMIC radio occultation excess phase using undifferenced measurements

2017 ◽  
Vol 10 (5) ◽  
pp. 1813-1821
Author(s):  
Pengfei Xia ◽  
Shirong Ye ◽  
Kecai Jiang ◽  
Dezhong Chen
Keyword(s):  
Radio Occultation ◽  
Weather Prediction ◽  
Lower Troposphere ◽  
Double Difference ◽  
Cosmic Radio ◽  
Processing Strategy ◽  
Single Difference ◽  
Excess Phase ◽  
Mean Square Errors ◽  
Better Than

Abstract. In the GPS radio occultation technique, the atmospheric excess phase (AEP) can be used to derive the refractivity, which is an important quantity in numerical weather prediction. The AEP is conventionally estimated based on GPS double-difference or single-difference techniques. These two techniques, however, rely on the reference data in the data processing, increasing the complexity of computation. In this study, an undifferenced (ND) processing strategy is proposed to estimate the AEP. To begin with, we use PANDA (Positioning and Navigation Data Analyst) software to perform the precise orbit determination (POD) for the purpose of acquiring the position and velocity of the mass centre of the COSMIC (The Constellation Observing System for Meteorology, Ionosphere and Climate) satellites and the corresponding receiver clock offset. The bending angles, refractivity and dry temperature profiles are derived from the estimated AEP using Radio Occultation Processing Package (ROPP) software. The ND method is validated by the COSMIC products in typical rising and setting occultation events. Results indicate that rms (root mean square) errors of relative refractivity differences between undifferenced and atmospheric profiles (atmPrf) provided by UCAR/CDAAC (University Corporation for Atmospheric Research/COSMIC Data Analysis and Archive Centre) are better than 4 and 3 % in rising and setting occultation events respectively. In addition, we also compare the relative refractivity bias between ND-derived methods and atmPrf profiles of globally distributed 200 COSMIC occultation events on 12 December 2013. The statistical results indicate that the average rms relative refractivity deviation between ND-derived and COSMIC profiles is better than 2 % in the rising occultation event and better than 1.7 % in the setting occultation event. Moreover, the observed COSMIC refractivity profiles from ND processing strategy are further validated using European Centre for Medium-Range Weather Forecasts (ECMWF) analysis data, and the results indicate that the undifferenced method reduces the noise level on the excess phase paths in the lower troposphere compared to the single-difference processing strategy.

scholarly journals Retrieval of temperature and water vapor profiles from radio occultation refractivity and bending angle measurements using an Optimal Estimation approach: a simulation study

2005 ◽  
Vol 5 (6) ◽  
pp. 1665-1677 ◽  
Author(s):  
A. von Engeln ◽  
G. Nedoluha
Keyword(s):  
Water Vapor ◽  
Radio Occultation ◽  
A Priori ◽  
Estimation Method ◽  
Lower Troposphere ◽  
Optimal Estimation ◽  
Pronounced Difference ◽  
Bending Angle ◽  
Advantages And Disadvantages ◽  
The Impact

Abstract. The Optimal Estimation Method is used to retrieve temperature and water vapor profiles from simulated radio occultation measurements in order to assess how different retrieval schemes may affect the assimilation of this data. High resolution ECMWF global fields are used by a state-of-the-art radio occultation simulator to provide quasi-realistic bending angle and refractivity profiles. Both types of profiles are used in the retrieval process to assess their advantages and disadvantages. The impact of the GPS measurement is expressed as an improvement over the a priori knowledge (taken from a 24h old analysis). Large improvements are found for temperature in the upper troposphere and lower stratosphere. Only very small improvements are found in the lower troposphere, where water vapor is present. Water vapor improvements are only significant between about 1 km to 7 km. No pronounced difference is found between retrievals based upon bending angles or refractivity. Results are compared to idealized retrievals, where the atmosphere is spherically symmetric and instrument noise is not included. Comparing idealized to quasi-realistic calculations shows that the main impact of a ray tracing algorithm can be expected for low latitude water vapor, where the horizontal variability is high. We also address the effect of altitude correlations in the temperature and water vapor. Overall, we find that water vapor and temperature retrievals using bending angle profiles are more CPU intensive than refractivity profiles, but that they do not provide significantly better results.

scholarly journals A Quality Control Procedure Based on Bending Angle Measurement Uncertainty for Radio Occultation Data Assimilation in the Tropical Lower Troposphere

2018 ◽  
Vol 35 (10) ◽  
pp. 2117-2131 ◽  
Author(s):  
Hui Liu ◽  
Ying-Hwa Kuo ◽  
Sergey Sokolovskiy ◽  
Xiaolei Zou ◽  
Zhen Zeng ◽  
...  
Keyword(s):  
Quality Control ◽  
Data Assimilation ◽  
Radio Occultation ◽  
Lower Troposphere ◽  
Spectral Width ◽  
Angle Measurement ◽  
Control Procedure ◽  
Random Errors ◽  
Bending Angle ◽  
Occultation Data

AbstractThe fluctuation of radio occultation (RO) signals in the presence of refractivity irregularities in the moist lower troposphere results in uncertainties of retrieved bending angle and refractivity profiles. In this study the local spectral width (LSW) of RO signals, transformed to impact parameter representation, is used for the characterization of the uncertainty (random error) of retrieved bending angle and refractivity profiles. A large LSW has some correlation with the large mean difference (bias) of retrieved refractivity and bending angle from radiosondes and European Centre for Medium-Range Weather Forecasts analyses based on data from 2008 to 2014. An LSW-based quality control (QC) procedure is developed to eliminate low-quality (large random errors and biases) profiles from data assimilation. The LSW-based QC procedure is tested and evaluated in the assimilation of Constellation Observing System for Meteorology, Ionosphere and Climate RO data using the NCAR Data Assimilation Research Testbed and the Weather Research and Forecasting Model. Preliminary results, based on a 2-week data assimilation cycle, show that the LSW-based QC procedure improves water vapor analyses in the moist lower troposphere.

Sign in / Sign up

Export Citation Format

Share Document