Radio Occultation Observations Research Articles

PlanetiQ Radio Occultation: Preliminary Comparative Analysis of Neutral Profiles vs. COSMIC and NWP Models

Radio Occultation (RO) is pivotal for profiling the neutral and ionized atmosphere, with the PlanetiQ mission, via its GNOMES satellites, striving to establish an advanced atmospheric observing system. However, an assessment of the spatiotemporal distributions of PlanetiQ observations and comparisons with reliable datasets are lacking. This study addresses this gap by examining the temporal and spatial distribution of RO observations from PlanetiQ during its initial 198 operational days in 2023, alongside comparisons with COSMIC and Numerical Weather Prediction (NWP) models. Data from GN02, GN03, and GN04 satellites, yielding 1099, 1313, and 1843 RO events per day, respectively, were analyzed. The satellite constellation’s observations demonstrate a generally well-distributed pattern, albeit minor deficiencies in equatorial and polar regions. Single-profile comparisons with COSMIC data reveal strong correlations for pressure, temperature, Water Vapor Pressure (WVP), and refractivity profiles, with temperature exhibiting larger variations (RMSE = 1.24 °C). Statistical analyses confirm statistically insignificant differences between the PlanetiQ and COSMIC profiles at the same spatio-temporal coordinates. Comparisons with NWP models show slight differences with GFS, with overall RMSE values of 0.23 mb (WVP), 0.6 mb (pressure), 1.3 (refractivity), and 1.5 °C (temperature). However, assessments against GFS/ECMWF models indicate overall compatibility, with insignificant differences between PlanetiQ profiles and model observations.

The Study of Internal Gravity Waves in the Earth’s Atmosphere by Radio Occultations: A Review

Internal gravity waves (IGWs) in the middle atmosphere are the main source of mesoscale fluctuations of wind and temperature. The parameterization of IGWs and study of their climatology is necessary for the development of global atmospheric circulation models. In this review, we focus on the application of Radio Occultation (RO) observations for the retrieval of IGW parameters. (1) The simplest approach employs the retrieved temperature profiles. It is based on the fact that IGWs are highly anisotropic structures and can be accurately retrieved by RO. The basic assumption is that all the temperature fluctuations are caused by IGWs. The smoothed background temperature profile defines the the Brunt–Väisälä frequency, which, together with the temperature fluctuations, defines the IGW specific potential energy. Many studies have derived the distribution and climatology of potential energy, which is one of the most important characteristics of IGWs. (2) More detailed analysis of the temperature profiles is based on the derivation of the temperature fluctuation spectra. For saturated IGWs, the spectra must obey the power law with an exponent of −3. Such spectra are obtained by using Wave Optical (WO) processing. (3) More advanced analysis employs space–frequency analysis. It is based on phase-sensitive techniques like cross S- or wavelet transforms in order to identify propagating IGWs. (4) Another direction is the IGW parameter estimate from separate temperature profiles applying the stability condition in terms of the Richardson number. In this framework, a necessary condition is formulated that defines whether or not the temperature fluctuations can be related to IGW events. The temperature profile retrieval involves integral transforms and filtering that constitute the observation filter. (5) A simpler filter is implemented by the analysis of the RO amplitude fluctuation spectra, based on the diffraction theory in the framework of the phase screen and weak fluctuation approximations. The two spectral parameters, the external scale and the structural characteristic, define the specific potential energy. This approach allows the derivation of the spacial and seasonal distributions of IGW activity. We conclude that the success of IGW study by RO is stimulated by a large number of RO observations and advanced techniques based on Fourier and space–time analysis, physical equations describing IGWs, and diffraction theory.

Monitoring urban heat island intensity based on GNSS tomography technique

Monitoring urban heat island (UHI) effect is critical because it causes health problems and excessive energy consumption more energy when cooling buildings. In this study, we propose an approach for UHI monitoring by fusing data from ground-based global navigation satellite system (GNSS), space-based GNSS radio occultation (RO), and radiosonde. The idea of the approach is as follows: First, the first and second grid tops are defined based on historical RO and radiosonde observations. Next, the wet refractivities between the first and second grid tops are fitted to higher-order spherical harmonics and they are used as the inputs of GNSS tomography. Then, the temperature and water vapor partial pressure are estimated by using best search method based on the tomography-derived wet refractivity. In the end, the UHI intensity is evaluated by calculating the temperature difference between the urban regions and nearby rural regions. Feasibility of the UHI intensity monitoring approach was evaluated with GNSS RO and radiosonde data in 2010–2019, as well as ground-based GNSS data in 2020 in Hong Kong, China, by taking synoptic temperature data as reference. The result shows that the proposed approach achieved an accuracy of 1.2 K at a 95% confidence level.

Comparison of the Heights of Sporadic E Layers and Vertical Ion Convergence Parameters

Sporadic E (Es) layers are thin layers of enhanced electron density that commonly appear at altitudes of 90–130 km, often impacting radio communications and navigation systems. The wind shear theory posits that the vertical ion drift, influenced by atmospheric neutral winds and the magnetic field, serves as a significant dynamic driver for the formation and movement of Es layers. In current studies, both the heights of ion vertical velocity null (IVN) and the maximum vertical ion convergence (VICmax) have been proposed as the potential height of Es layer occurrence. In this study, utilizing the neutral atmospheric wind data derived from the WACCM-X (The Whole Atmosphere Community Climate Model with thermosphere and ionosphere extension), we computed and compared these two parameters with the observed Es layer heights recorded by the FORMOSAT-3/COSMIC (FORMOsa SATellite-3/Constellation Observing System for Meteorology, Ionosphere, and Climate) radio occultation (RO) observations. The comparative analysis suggests that IVN is a more likely node for Es layer occurrence than VICmax. Subsequently, we examined the height–time distributions of IVN and Es layers, as well as their respective descent rates at different latitudes. These results demonstrated a notable agreement in height variations between IVN and Es layers. The collective results presented in this paper provide strong support that the ion vertical velocity null plays a crucial role in determining the height of Es layers.

Variations in the Tropopause Temperature and Height over the Vu Gia-Thu Bon River Basin in Vietnam: Insights from GNSS Radio Occultation Observations

Variations in the Tropopause Temperature and Height over the Vu Gia-Thu Bon River Basin in Vietnam: Insights from GNSS Radio Occultation Observations

Analysis of Two Selected Solar Events in 2011 and 2015 With Mars Express Radio Occultation Data

AbstractThe temporal behavior of the Martian ionosphere is highly variable due to various dynamic processes including space weather events. Here, we study the effect of solar flares and coronal mass ejections (CMEs) on the Martian ionosphere for two selected solar events in 2011 and 2015, using the publicly available Mars EXpress (MEX) radio occultation (RO) data (MaRS). We developed a data processing software that converts the calibrated Radio Occultation (RO) Doppler data to scientifically valuable atmospheric profiles. Using this software and previously unexplored MaRS observations, the variations in ionospheric parameters (electron density profiles and total electron content (TEC)) are calculated in order to evaluate the ionospheric changes due to solar flares and CMEs. The RO measurements mostly available 1–4 days apart from the peak events, showed no evident change in the TEC nor in the shape of electron density profiles except for a possible gradual decrease in altitude of M2 (main layer) peak density following the arrival of CMEs. To better quantify the effect of solar events on electron density profiles, RO observations near the time of arrival of solar flares and CMEs are crucial. This can be achieved by frequent RO measurements by various Mars orbiters including spacecraft‐to‐spacecraft measurements assisted by multi‐instrument monitoring of the ionosphere.

Preliminary Estimations of Mars Atmospheric and Ionospheric Profiles from Tianwen-1 Radio Occultation One-Way, Two-Way, and Three-Way Observations

The radio occultation method, one of the methods used to provide planetary atmospheric profiles with high vertical resolution, was applied to China’s first Mars mission, Tianwen-1. We carried out observations based on the Chinese Deep Space Network, and one-way, two-way, and three-way modes were used for Doppler observations from the Tianwen-1 spacecraft. We successfully obtained effective observations from Tianwen-1 on 22 and 25 March 2022. An inversion system developed for Tianwen-1 radio occultation observations enabled the derivation of neutral atmospheric density, pressure, temperature, and electron density profiles of Mars. Utilizing one-way tracking data, Martian ionospheric electron density profiles were retrieved at latitudes between 68.7 and 70.7 degrees (N). However, the presence of strong random walk noise in one-way tracking data led to poor inversion results. Meanwhile, Martian ionospheric electron density and neutral atmosphere profiles were extracted from two-way and three-way tracking data at latitudes between 55.1 and 57.0 degrees (S) on 22 March and at latitudes between 62.8 and 63.4 degrees (S) on 25 March. Importantly, our inversion results from Tianwen-1 maintained consistency with results from the Mars Express and the Chapman theory (mainly in the M2 layer). Through two days’ observation experiments, we established and verified the occultation solution system and prepared for the follow-up occultation plans.

Limits on the Density of the Lunar Ionosphere: ARTEMIS Observations

The Moon of our Earth has a tenuous atmosphere, known as an exosphere. The ionization of this exosphere is speculated to possibly form a weak ionosphere. Some radio occultation (RO) experiments have suggested the presence of a dense ionosphere with an electron density on the order of hundreds of cm−3 near the surface. Using in situ measurements from the ARTEMIS mission during 2012–2021, we conduct statistical analyses and case studies to investigate the plasma density at near-surface altitudes. ARTEMIS measurements reveal no plasma densities at altitudes between 10 and 50 km that exceed 35 cm−3, and therefore they provide no evidence for a steady-state or global lunar ionosphere at the level suggested by some RO observations. Density profiles with local time and altitude show higher density in the sunlit sector than in the shadowed sector. These observations suggest that the natural variation of solar wind plasma flux with solar zenith angle plays a critical role in controlling the plasma population near the surface. This research provides a reference for a comparison with RO observations and a statistical view of the low-altitude plasma environment near the lunar surface.

Quality Assessment of the Atmospheric Radio Occultation Profiles from FY-3E/GNOS-II BDS and GPS Measurements

The Fengyun-3E (FY-3E) satellite carrying the advanced Global Navigation Satellite System (GNSS) Radio Occultation Sounder-II (GNOS-II) is already in operation for radio occultation (RO) observation, with the BeiDou Navigation Satellite System (BDS-2 and BDS-3) and Global Positioning System (GPS) signals tracking capability. FY-3E BDS and GPS RO signals tracking capability were first evaluated by comparing their penetration depths, and then the quality of the refractivity, temperature, and specific humidity profiles was analyzed with the fifth-generation European Centre for Medium-Range Weather Forecasts reanalysis (ERA5) data. Results show the mean penetration depth of BDS occultations was 1.65 km compared to 1.09 km of GPS occultations. Between 5 and 25 km, the mean refractivity bias of the BDS (GPS) was −0.14% (0.01%) with the mean standard deviation (SD) being 1.11% (1.52%); the mean temperature biases of both were within ±0.1 K, and the mean SD of BDS was 1.1 K compared to 1.2 K for the GPS; BDS/GPS specific humidity bias was within ±0.3 g/kg with corresponding SD being less than 1.3 g/kg. Seasonal deviations of specific humidities were largest in summer and smallest in winter. Latitudinal deviations over the tropics were generally higher than in other areas. Enriched quantity and high accuracy and precision after careful calibration will promote the FY-3E RO profiles as a reliable data source for the RO community.

Estimating Individual Radio Occultation Uncertainties Using the Observations and Environmental Parameters

Abstract Estimation of uncertainties (random error statistics) of radio occultation (RO) observations is important for their effective assimilation in numerical weather prediction (NWP) models. Average uncertainties can be estimated for large samples of RO observations and these statistics may be used for specifying the observation errors in NWP data assimilation. However, the uncertainties of individual RO observations vary, and so using average uncertainty estimates will overestimate the uncertainties of some observations and underestimate those of others, reducing their overall effectiveness in the assimilation. Several parameters associated with RO observations or their atmospheric environments have been proposed to estimate individual RO errors. These include the standard deviation of bending angle (BA) departures from either climatology in the upper stratosphere and lower mesosphere (STDV) or the sample mean between 40 and 60 km (STD4060), the local spectral width (LSW), and the magnitude of the horizontal gradient of refractivity (|∇HN|). In this paper we show how the uncertainties of two RO datasets, COSMIC-2 and Spire BA, as well as their combination, vary with these parameters. We find that the uncertainties are highly correlated with STDV and STD4060 in the stratosphere, and with LSW and |∇HN| in the lower troposphere. These results suggest a hybrid error model for individual BA observations that uses an average statistical model of RO errors modified by STDV or STD4060 above 30 km, and LSW or |∇HN| below 8 km. Significance Statement These results contribute to the understanding of the sources of uncertainties in radio occultation observations. They could be used to improve the effectiveness of these observations in their assimilation into numerical weather prediction and reanalysis models by improving the estimation of their observational errors.

Quality assessment of the ionospheric density profiles based on long-term COSMIC 1 and 2 radio occultation observations

Quality assessment of the ionospheric density profiles based on long-term COSMIC 1 and 2 radio occultation observations

Preliminary Assessment of BDS Radio Occultation Retrieval Quality and Coverage Using FY-3E GNOS II Measurements

The FengYun-3E Global Navigation Satellite System (GNSS) occultation sounder II (FY-3E GNOS II) was launched on 5 July 2021. For the first time, based on the new GNOS II sensor, this mission realizes radio occultation (RO) and reflectometry observations using the navigation signals from the third-generation BeiDou System (BDS-3), and it is hence important to assess and analyze the BDS-3 remote sensing performances relative to other systems. Here, we assessed FY-3E GNOS II RO atmospheric retrievals by inter-comparing with corresponding data from the NCEP FNL global atmospheric analysis and FY-3D GNOS mission. The GNOS RO data quality and consistency of the different FY-3 meteorological satellites, i.e., FY-3D and FY-3E, as well as different GNSS systems (GPS, BDS-2, BDS-3) were analyzed. We find that the FY-3E GNOS II RO data exhibit better quality than FY-3D GNOS, particularly in the number, penetration height toward surface, and global coverage by BDS RO profiles, due to the integration of BDS-2 and BDS-3. Additionally, comparing with co-located NCEP FNL analysis profiles, the mean difference (and standard deviation) of the FY-3E GNOS II RO atmospheric refractivity profile retrievals is found to be smaller than 0.2% (and 1%), in the upper troposphere and lower stratosphere, from 5 to 30 km, and remains consistent at this accuracy and precision level with the FY-3D GNOS RO data. These features provide clear evidence for a high utility of the new GNOS II RO data for weather and climate research and applications.

Recent Impact of COSMIC-2 with Improved Radio Occultation Data Assimilation Algorithms

Abstract A Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) follow-on constellation, COSMIC-2, was successfully launched into equatorial orbit on 24 June 2019. With an increased signal-to-noise ratio due to improved receivers and digital beam-steering antennas, COSMIC-2 is producing about 5000 high-quality radio occultation (RO) profiles daily over the tropics and subtropics. The initial evaluation of the impact of assimilating COSMIC-2 into NOAA’s Global Forecast System (GFS) showed mixed results, and adjustments to quality control procedures and observation error characteristics had to be made prior to the assimilation of this dataset in the operational configuration in May 2020. Additional changes in the GFS that followed this initial operational implementation resulted in a larger percentage of rejection (∼90%) of all RO observations, including COSMIC-2, in the mid- to lower troposphere. Since then, two software upgrades directly related to the assimilation of RO bending angle observations were developed. These improvements were aimed at optimizing the utilization of COSMIC-2 and other RO observations to improve global weather analyses and forecasts. The first upgrade was implemented operationally in September 2021 and the second one in November 2022. This study describes both RO software upgrades and evaluates the impact of COSMIC-2 with this most recently improved configuration. Specifically, we show that the assimilation of COSMIC-2 observations has a significant impact in improving temperature and winds in the tropics, though benefits also extend to the extratropical latitudes.

Evaluation of new radio occultation observations among small satellites at Venus by data assimilation

We conducted observing system simulation experiments (OSSEs) for radio occultation measurements (RO) among small satellites, which are expected to be useful for future Venus missions. The effectiveness of the observations based on realistic orbit calculations was evaluated by reproduction of the “cold collar”, a unique thermal structure in the polar atmosphere of Venus. Pseudo-temperature observations for the OSSEs were provided from the Venus atmospheric GCM in which the cold collar was reproduced by the thermal forcing. The vertical temperature distributions between 40 and 90 km altitudes at observation points were assimilated. The result showed that the cold collar was most clearly reproduced in the case where the temperature field in high-latitudes was observed twice a day, suggesting that the proposed observation is quite effective to improve the polar atmospheric structure at least. Although the cold collar was also reproduced in the OSSEs for Longwave Infrared Camera (LIR) observations, the result seemed unrealistic and inefficient compared to that obtained in the RO OSSEs. The present study shows that the OSSEs can be used to evaluate observation plans and instruments in terms of reproducibility of specific atmospheric phenomena, and applied to future missions targeting planetary atmospheres.

An evaluation of GNSS radio occultation atmospheric profiles from Sentinel-6

The Global Navigation Satellite System (GNSS) radio occultation (RO) significantly impacts climate change monitoring. GNSS RO is based mainly on the refraction of the GNSS signals transmitted by satellites and received by a receiver on a Low Earth Orbit (LEO) satellite. Sentinel-6 (S6) is the first satellite from the Sentinel missions to retrieve atmospheric parameters using the RO technique. Temperature, pressure, and water vapor profiles can be retrieved using S6 based on the refraction of GNSS signals. This research evaluated the distribution and coverage of GNSS RO data from S6. This research also aims to evaluate the results of GNSS RO products, especially temperature from S6, from January 2022 to June 2022. The results of S6 are compared with the results from the Meteorological Operational Polar Satellite series (MetOp-B and MetOp-C) and the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis v5 (ERA5) dataset. Our results indicate that S6 allows for many RO observations in space by receiving signals from the Global Positioning System (GPS) and Russia's Global Navigation Satellite System (GLONASS). S6 provides monthly RO profiles from GPS that are greater than monthly RO profiles from GLONASS. Results show that S6 provides more GNSS RO data than MetOp-B and MetOp-C. The findings indicate that the distribution of GNSS RO events from S6 over the globe is uniform in longitude and non-uniform at latitude, with the minimum number at the poles. This paper makes a significant contribution to exploring the atmosphere with the S6 as a promising satellite for weather forecasting and monitoring climate change.

A New Method of Frequency Fluctuation Estimation and IPS Processing Results Based on the Downlink Signal of Tianwen-1

The radio-occultation observations taken by Tianwen-1 are aiming to study the properties of solar wind. A new method of frequency fluctuation (FF) estimation is presented for processing the down-link signals of Tianwen-1 during the occultation period to study the properties of the coronal plasma at the heliocentric distances of 4.48–19 R ⊙. Because of low S/N as well as the phase fluctuation phenomena caused by solar activity, a Kalman based on polynomial prediction methods is proposed to avoid the phase locked loop loss lock. A new detrend method based on multi-level iteration correction is proposed to estimate Doppler shift to get more accurate power density spectra of FF in the low frequency region. The data analyze procedure is used to get the properties of the solar corona during the occultation. The method was finally verified at the point when the solar offset is 5.7 R ⊙, frequency tracking was successfully performed on data with a carrier-to-noise ratio of about 28 dBHz. The density spectra obtained by the improved method are basically the same when the frequency is greater than 2 mHz, the uncertainty in the result of the rms of the FF obtained by removing the trend term with different order polynomials is less than 3.3%. The data without eliminating interference show a large error for different detrending orders, which justifies the need for an improved approach. Finally, the frequency fluctuation results combined with the information on intensity fluctuation obtained by the new method are compared with the results of the integrated Space Weather Analysis system and theoretical formula, which verifies that the processing results in this paper have a certain degree of credibility.

Assessment of GRAS Ionospheric Measurements for Ionospheric Model Assimilation

We conducted a study to assess the GNSS (Global Navigation Satellite System) Receiver for Atmospheric Sounding (GRAS) ionospheric data quality by processing Radio Occultation (RO) observations of ionospheric products. The main objective of the study is to validate ionospheric data generated at EUMETSAT, such as ionospheric bending angle profiles, amplitude and phase scintillations, topside Total Electron Content (TEC) from MetOp-A GRAS instrument as well as generating and validating new ionospheric products derived from GRAS RO observations such as the TEC, rate of TEC and vertical electron density profiles. The assessment is conducted by comparing and evaluating the systematic differences between similar products from other Low Earth Orbit (LEO) satellite missions or from ground-based ionospheric measurements. The study confirms that the GNSS topside and RO observations recorded by the GRAS instrument onboard MetOp satellites are of good quality and are a valuable source of data for ionospheric research.

The Four‐Dimensional Variational Neustrelitz Electron Density Assimilation Model: NEDAM

AbstractWith the development of the Global Navigation Satellite System (GNSS) and regional augmentation systems, the influence of the ionosphere on GNSS signals is becoming increasingly important. It is also of particular interest to retrieve the electron density distribution from GNSS observations. We have successfully developed a four‐dimensional variational assimilation scheme, the Neustrelitz Electron Density Assimilation Model (NEDAM), and verified NEDAM by a simulation study using a European ground‐based GNSS network. The performance of NEDAM is validated using two ionosondes and COSMIC‐1 radio occultation observations during the September 2017 geomagnetic storm period. The critical frequency of the F2 layer in NEDAM is much more accurate than that of a physics‐based model driven by observed geophysical indices or the Neustrelitz Electron Density Model (NEDM), when compared to data from two ionosondes. During the storm, the root mean square error of the F2‐layer critical frequency with respect to the two ionosondes is improved by 0.54 and 0.42 MHz, respectively. We also compare two co‐located electron density profiles from the COSMIC‐1 mission with NEDAM and NEDM. It is found that NEDAM is able to reconstruct well not only the peak density but also the peak density height, which is missing in the previous research.

An Improved One-Dimensional Bending Angle Forward Operator for the Assimilation of Radio Occultation Profiles in the Lower Troposphere

Abstract Under very large vertical gradients of atmospheric refractivity, which are typical at the height of the planetary boundary layer, the assimilation of radio occultation (RO) observations into numerical weather prediction (NWP) models presents several serious challenges. In such conditions, the assimilation of RO bending angle profiles is an ill-posed problem, the uncertainty associated with the RO observations is higher, and the one-dimensional forward operator used to assimilate these observations has several theoretical deficiencies. As a result, a larger percentage of these RO observations are rejected at the NWP centers by existing quality control procedures, potentially limiting the benefits of this data type to improve weather forecasting in the lower troposphere. To address these problems, a new methodology that enables the assimilation of RO data to be extended to the lower moist troposphere has been developed. Challenges associated with larger atmospheric gradients of refractivity are partially overcome by a reformulation that has minimal effect at higher altitudes. As a first step toward this effort, this study presents both the theoretical development of this new methodology and a forecast impact assessment of it using the NCEP NWP system. Though using a conservative approach, benefits in the lower tropical troposphere are already noticeable. The encouraging results of this work open the potential for further exploitation and optimization of RO assimilation.

Effect of GNSS radio occultation observations on the prediction of the 2021 Henan rainstorm

Effect of GNSS radio occultation observations on the prediction of the 2021 Henan rainstorm