scholarly journals Smoothing radio occultation bending angles above 40 km

2001 ◽  
Vol 19 (4) ◽  
pp. 459-468 ◽  
Author(s):  
S. B. Healy

Abstract. The 'statistically optimal' approach to smoothing bending angles derived from radio occultation (RO) measurements is outlined. This combines a measured bending angle profile with an a priori or background estimate derived from climatology, in order to obtain the most probable bending angle profile. However, the method is only optimal if the error statistics of both the measured and background profiles are known and applied accurately. In this work it is shown that correlations in the background estimate have a significant role in determining the degree of smoothing in the solution. We find that smooth profiles, consistent with the measured values, can be derived if the correlations are approximated analytically with a Gaussian, assuming a scale length of 6km. In regions where the observed and background error levels are comparable, the solutions take the general shape from the background estimate, centred on the observation data. The effects of correlated observation errors are also considered. It is shown that the quality of the temperature retrievals can be significantly affected by the choice of climatology used for background estimate.Key words. Atmosphere composition and structure (pressure, density and temperature) – Radio science (remote sensing)

2014 ◽  
Vol 7 (8) ◽  
pp. 8193-8231 ◽  
Author(s):  
B. Scherllin-Pirscher ◽  
S. Syndergaard ◽  
U. Foelsche ◽  
K. B. Lauritsen

Abstract. In this paper, we introduce a bending angle radio occultation climatology (BAROCLIM) based on Formosat-3/COSMIC (F3C) data. This climatology represents the monthly-mean atmospheric state from 2006 to 2012. Bending angles from radio occultation (RO) measurements are obtained from the accumulation of the change in the raypath direction of Global Positioning System (GPS) signals. Best quality of these near-vertical profiles is found from the middle troposphere up to the mesosphere. Beside RO bending angles we also use data from the Mass Spectrometer and Incoherent Scatter Radar (MSIS) model to expand BAROCLIM in a spectral model, which (theoretically) reaches from the surface up to infinity. Due to the very high quality of BAROCLIM up to the mesosphere, it can be used to detect deficiencies in current state-of-the-art analysis and reanalysis products from numerical weather prediction (NWP) centers. For bending angles derived from European Centre for Medium-Range Weather Forecasts (ECMWF) analysis fields from 2006 to 2012, e.g., we find a positive bias of 0.5% to % at 40 km, which increases to more than 2% at 50 km. BAROCLIM can also be used as a priori information in RO profile retrievals. In contrast to other a priori information (i.e., MSIS) we find that the use of BAROCLIM better preserves the mean of raw RO measurements. Global statistics of statistically optimized bending angle and refractivity profiles also confirm that BAROCLIM outperforms MSIS. These results clearly demonstrate the utility of BAROCLIM.


2015 ◽  
Vol 8 (1) ◽  
pp. 109-124 ◽  
Author(s):  
B. Scherllin-Pirscher ◽  
S. Syndergaard ◽  
U. Foelsche ◽  
K. B. Lauritsen

Abstract. In this paper, we introduce a bending angle radio occultation climatology (BAROCLIM) based on Formosat-3/COSMIC (F3C) data. This climatology represents the monthly-mean atmospheric state from 2006 to 2012. Bending angles from radio occultation (RO) measurements are obtained from the accumulation of the change in the raypath direction of Global Positioning System (GPS) signals. Best quality of these near-vertical profiles is found from the middle troposphere up to the mesosphere. Beside RO bending angles we also use data from the Mass Spectrometer and Incoherent Scatter Radar (MSIS) model (modified for RO purposes) to expand BAROCLIM in a spectral model, which (theoretically) reaches from the surface up to infinity. Due to the very high quality of BAROCLIM up to the mesosphere, it can be used to detect deficiencies in current state-of-the-art analysis and reanalysis products from numerical weather prediction (NWP) centers. For bending angles derived from European Centre for Medium-Range Weather Forecasts (ECMWF) analysis fields from 2006 to 2012, e.g., we find a positive bias of 0.5 to 1% at 40 km, which increases to more than 2% at 50 km. BAROCLIM can also be used as a priori information in RO profile retrievals. In contrast to other a priori information (i.e., MSIS) we find that the use of BAROCLIM better preserves the mean of raw RO measurements. Global statistics of statistically optimized bending angle and refractivity profiles also confirm that BAROCLIM outperforms MSIS. These results clearly demonstrate the utility of BAROCLIM.


2005 ◽  
Vol 5 (6) ◽  
pp. 1665-1677 ◽  
Author(s):  
A. von Engeln ◽  
G. Nedoluha

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.


2018 ◽  
Vol 11 (4) ◽  
pp. 2427-2440 ◽  
Author(s):  
Congliang Liu ◽  
Gottfried Kirchengast ◽  
Yueqiang Sun ◽  
Kefei Zhang ◽  
Robert Norman ◽  
...  

Abstract. The Global Navigation Satellite System (GNSS) radio occultation (RO) technique is widely used to observe the atmosphere for applications such as numerical weather prediction and global climate monitoring. The ionosphere is a major error source to RO at upper stratospheric altitudes, and a linear dual-frequency bending angle correction is commonly used to remove the first-order ionospheric effect. However, the higher-order residual ionospheric error (RIE) can still be significant, so it needs to be further mitigated for high-accuracy applications, especially from 35 km altitude upward, where the RIE is most relevant compared to the decreasing magnitude of the atmospheric bending angle. In a previous study we quantified RIEs using an ensemble of about 700 quasi-realistic end-to-end simulated RO events, finding typical RIEs at the 0.1 to 0.5 µrad noise level, but were left with 26 exceptional events with anomalous RIEs at the 1 to 10 µrad level that remained unexplained. In this study, we focused on investigating the causes of the high RIE of these exceptional events, employing detailed along-ray-path analyses of atmospheric and ionospheric refractivities, impact parameter changes, and bending angles and RIEs under asymmetric and symmetric ionospheric structures. We found that the main causes of the high RIEs are a combination of physics-based effects – where asymmetric ionospheric conditions play the primary role, more than the ionization level driven by solar activity – and technical ray tracer effects due to occasions of imperfect smoothness in ionospheric refractivity model derivatives. We also found that along-ray impact parameter variations of more than 10 to 20 m are possible due to ionospheric asymmetries and, depending on prevailing horizontal refractivity gradients, are positive or negative relative to the initial impact parameter at the GNSS transmitter. Furthermore, mesospheric RIEs are found generally higher than upper-stratospheric ones, likely due to being closer in tangent point heights to the ionospheric E layer peaking near 105 km, which increases RIE vulnerability. In the future we will further improve the along-ray modeling system to fully isolate technical from physics-based effects and to use it beyond this work for additional GNSS RO signal propagation studies.


2021 ◽  
Vol 13 (18) ◽  
pp. 3644
Author(s):  
Yong Chen ◽  
Xi Shao ◽  
Changyong Cao ◽  
Shu-peng Ho

The Global Navigation Satellite System (GNSS) radio occultation (RO) is a remote sensing technique that uses International System of Units (SI) traceable GNSS signals for atmospheric limb soundings. The RO bending angle/sounding profiles are needed for assimilation in Numerical Weather Prediction (NWP) models, weather, climate, and space weather applications. Evaluating these RO data to ensure the high data quality for these applications is becoming more and more critical. This study presents a method for predicting radio occultation events, from which simultaneous radio occultation (SRO) for a pair of low-Earth-orbit (LEO) satellites on the limb to the same GNSS satellite can be obtained. The SRO method complements the Simultaneous Nadir Overpass (SNO) method (for nadir viewing satellite instruments), which has been widely used to inter-calibrate LEO to LEO and LEO to geosynchronous-equatorial-orbit (GEO) satellites. Unlike the SNO method, the SRO method involves three satellites: a GNSS and two LEO satellites with RO receivers. The SRO method allows for the direct comparison of bending angles when the simultaneous RO measurements for two LEO satellites receiving the same GNSS signal pass through approximately the same atmosphere within minutes in time and within less than 200 km of distance from each other. The prediction method can also be applied to radiosonde overpass prediction, and coordinate radiosonde launches for inter-comparisons between RO and radiosonde profiles. The main advantage of the SRO comparisons of bending angles is the significantly reduced uncertainties due to the much shorter time and smaller atmospheric path differences than traditional RO comparisons. To demonstrate the usefulness of this method, we present a comparison of the Constellation Observing System for Meteorology, Ionosphere, and Climate-2 (COSMIC-2) and GeoOpitcs RO profiles using SRO data for two time periods: Commercial Weather Data (CWD) data delivery order-1 (DO-1): 15 December 2020–15 January 2021 and CWD DO-2: 17 March 2021–31 August 2021. The results show good agreement in the bending angles between the COSMIC-2 RO measurements and those from GeoOptics, although systematic biases are also found in the inter-comparisons. Instrument and processing algorithm performances for the signal-to-noise ratio (SNR), penetration height, and bending angle retrieval uncertainty are also characterized. Given the efficiency of this method and the many RO measurements that are publicly and commercially available as well as the expansion of receiver capabilities to all GNSS systems, it is expected that this method can be used to validate/inter-calibrate GNSS RO measurements from different missions.


2016 ◽  
Vol 9 (2) ◽  
pp. 335-346 ◽  
Author(s):  
Z. Zeng ◽  
S. Sokolovskiy ◽  
W. Schreiner ◽  
D. Hunt ◽  
J. Lin ◽  
...  

Abstract. For inversions of the GPS radio occultation (RO) data in the neutral atmosphere, this study investigates an optimal transition height for replacing the standard ionospheric correction using the linear combination of the L1 and L2 bending angles with the correction of the L1 bending angle by the L1–L2 bending angle extrapolated from above. The optimal transition height depends on the RO mission (i.e., the receiver and firmware) and is different between rising and setting occultations and between L2P and L2C GPS signals. This height is within the range of approximately 10–20 km. One fixed transition height, which can be used for the processing of currently available GPS RO data, can be set to 20 km. Analysis of the L1CA and the L2C bending angles shows that in some occultations the errors of standard ionospheric correction substantially increase around the strong inversion layers (such as the top of the boundary layer). This error increase is modeled and explained by the horizontal inhomogeneity of the ionosphere.


2015 ◽  
Vol 8 (7) ◽  
pp. 2999-3019 ◽  
Author(s):  
C. L. Liu ◽  
G. Kirchengast ◽  
K. Zhang ◽  
R. Norman ◽  
Y. Li ◽  
...  

Abstract. The radio occultation (RO) technique using signals from the Global Navigation Satellite System (GNSS), in particular from the Global Positioning System (GPS) so far, is currently widely used to observe the atmosphere for applications such as numerical weather prediction and global climate monitoring. The ionosphere is a major error source in RO measurements at stratospheric altitudes, and a linear ionospheric correction of dual-frequency RO bending angles is commonly used to remove the first-order ionospheric effect. However, the residual ionospheric error (RIE) can still be significant so that it needs to be further mitigated for high-accuracy applications, especially above about 30 km altitude where the RIE is most relevant compared to the magnitude of the neutral atmospheric bending angle. Quantification and careful analyses for better understanding of the RIE is therefore important for enabling benchmark-quality stratospheric RO retrievals. Here we present such an analysis of bending angle RIEs covering the stratosphere and mesosphere, using quasi-realistic end-to-end simulations for a full-day ensemble of RO events. Based on the ensemble simulations we assessed the variation of bending angle RIEs, both biases and standard deviations, with solar activity, latitudinal region and with or without the assumption of ionospheric spherical symmetry and co-existing observing system errors. We find that the bending angle RIE biases in the upper stratosphere and mesosphere, and in all latitudinal zones from low to high latitudes, have a clear negative tendency and a magnitude increasing with solar activity, which is in line with recent empirical studies based on real RO data although we find smaller bias magnitudes, deserving further study in the future. The maximum RIE biases are found at low latitudes during daytime, where they amount to within −0.03 to −0.05 μrad, the smallest at high latitudes (0 to −0.01 μrad; quiet space weather and winter conditions). Ionospheric spherical symmetry or asymmetries about the RO event location have only a minor influence on RIE biases. The RIE standard deviations are markedly increased both by ionospheric asymmetries and increasing solar activity and amount to about 0.3 to 0.7 μrad in the upper stratosphere and mesosphere. Taking also into account the realistic observation errors of a modern RO receiving system, amounting globally to about 0.4 μrad (unbiased; standard deviation), shows that the random RIEs are typically comparable to the total observing system error. The results help to inform future RIE mitigation schemes that will improve upon the use of the linear ionospheric correction of bending angles and also provide explicit uncertainty estimates.


2015 ◽  
Vol 8 (1) ◽  
pp. 759-809 ◽  
Author(s):  
C. L. Liu ◽  
G. Kirchengast ◽  
K. Zhang ◽  
R. Norman ◽  
Y. Li ◽  
...  

Abstract. The radio occultation (RO) technique using signals from the Global Navigation Satellite System (GNSS), in particular from the Global Positioning System (GPS) so far, is meanwhile widely used to observe the atmosphere for applications such as numerical weather prediction and global climate monitoring. The ionosphere is a major error source in RO measurements at stratospheric altitudes and a linear ionospheric correction of dual-frequency RO bending angles is commonly used to remove the first-order ionospheric effect. However, the residual ionopheric error (RIE) can still be significant so that it needs to be further mitigated for high accuracy applications, especially above about 30 km altitude where the RIE is most relevant compared to the magnitude of the neutral atmospheric bending angle. Quantification and careful analyses for better understanding of the RIE is therefore important towards enabling benchmark-quality stratospheric RO retrievals. Here we present such an analysis of bending angle RIEs covering the stratosphere and mesosphere, using quasi-realistic end-to-end simulations for a full-day ensemble of RO events. Based on the ensemble simulations we assessed the variation of bending angle RIEs, both biases and SDs, with solar activity, latitudinal region, and with or without the assumption of ionospheric spherical symmetry and of co-existing observing system errors. We find that the bending angle RIE biases in the upper stratosphere and mesosphere, and in all latitudinal zones from low- to high-latitudes, have a clear negative tendency and a magnitude increasing with solar activity, in line with recent empirical studies based on real RO data. The maximum RIE biases are found at low latitudes during daytime, where they amount to with in −0.03 to −0.05 μrad, the smallest at high latitudes (0 to −0.01 μrad; quiet space weather and winter conditions). Ionospheric spherical symmetry or asymmetries about the RO event location have only a minor influence on RIE biases. The RIE SDs are markedly increased both by ionospheric asymmetries and increasing solar activity and amount to about 0.3 to 0.7 μrad in the upper stratosphere and mesosphere. Taking into account also realistic observation errors of a modern RO receiving system, amounting globally to about 0.4 μrad (un-biased; SD), shows that the random RIEs are typically comparable to the total observing system error. The results help to inform future RIE mitigation schemes that will improve upon the use of the linear ionospheric correction of bending angles and that will also provide explicit uncertainty estimates.


2017 ◽  
Author(s):  
Congliang Liu ◽  
Gottfried Kirchengast ◽  
Yueqiang Sun ◽  
Kefei Zhang ◽  
Robert Norman ◽  
...  

Abstract. The Global Navigation Satellite System (GNSS) radio occultation (RO) technique is widely used to observe the atmosphere for applications such as numerical weather prediction and global climate monitoring. The ionosphere is a major error source to RO at upper stratospheric altitudes and a linear dual-frequency bending angle correction is commonly used to remove the first-order ionospheric effect. However, the residual higher-order ionospheric error (RIE) can still be significant so that it needs to be further mitigated for high accuracy applications, especially from 30 km altitude upward where the RIE is most relevant compared to the decreasing magnitude of the atmospheric bending angle. In a previous study we quantified RIEs using an ensemble of about 700 quasi-realistic end-to-end simulated RO events, finding typical RIEs at the 0.1 to 0.5 μrad noise level, but were left with 26 exceptional events with anomalous RIEs at the 1 to 10 μrad level that remained unexplained. In this study, we focused on investigating the causes of the high RIE of these exceptional events, employing detailed along-raypath analyses of atmospheric and ionospheric refractivities, impact parameter changes, and bending angles and RIEs under asymmetric and symmetric ionospheric structures. We found that the main causes of the high RIEs are a combination of physics-based effects, where asymmetric ionospheric conditions play the primary role, more than the ionization level driven by solar activity, and technical ray tracer effects due to occasions of imperfect smoothness in ionospheric refractivity model derivatives. We also found that along-ray impact parameter variations of more than 10 to 20 m are well possible due to ionospheric asymmetries, and depending on prevailing horizontal refractivity gradients are positive or negative relative to the initial impact parameter at the GNSS transmitter. Furthermore, mesospheric RIEs are found generally higher than upper stratospheric ones, likely due to being closer in tangent point heights to the ionospheric E layer peaking near 105 km, which increases RIE vulnerability. In future we will further improve the along-ray modeling system to fully isolate technical from physics-based effects and to use it beyond this work for additional GNSS RO signal propagation studies.


2018 ◽  
Author(s):  
Julia Danzer ◽  
Marc Schwärz ◽  
Veronika Proschek ◽  
Ulrich Foelsche ◽  
Hans Gleisner

Abstract. Recently a new approach for the production of GNSS radio occultation climatologies has been proposed. The idea is to perform the averaging of individual profiles already in bending angle space and propagating the mean bending angle profiles through the Abel transform. Climatological products of refractivity, density, pressure, and temperature are directly retrieved from the mean bending angles. The averaging of a large number of profiles suppresses noise in the data, enabling observed bending angle data to be used up to 80 km without the need of a priori information. Above that altitude some background information for the Abel integral is still necessary. This work is a follow up study, having the focus on the comparison of the average profile inversion climatologies (API) from the two processing centers WEGC and DMI, studying monthly COSMIC data from January to March 2011. The impact of different backgrounds above 80 km is tested, and different implementations of the Abel integral are investigated. Results are compared for the climatological products against ECMWF analysis, MIPAS, and SABER data. It is shown that different implementations of the Abel integral have only little impact on the average profile inversion climatologies. On the other hand, different expansions of the bending angle profile above 80 km play a key role on the resulting monthly mean refractivities above 35 km altitude. Below that respective altitude the API climatologies show a good agreement between the two processing centers WEGC and DMI. Due to the downward propagation within the retrieval, effects of the upper initialization lead to differences in dry temperature climatologies already at 20 km altitude. Applying at both centers an exponential extrapolation to the bending angles above 80 km, dry temperature climatologies agree between WEGC, DMI, ECMWF analysis, and MIPAS up to 35 km altitude within ± 0.5 K, and up to 40 km altitude within ± 1 K. We conclude that up to the lower stratosphere the average profile inversion is a valid – and in computation time much faster – alternative for the production of dry atmospheric RO climatologies.


Sign in / Sign up

Export Citation Format

Share Document