Higher order ionospheric propagation effects on GPS radio occultation signals

With the increasing number of remote sensing satellites using the GPS radio occultation technique for atmospheric sounding, the estimation of higher order ionospheric effects and their mitigation have become relevant and important. Due to long ionospheric limb paths, GPS signals are strongly affected by ionospheric refraction during radio occultation. Standard dual-frequency GPS measurements may be used to estimate the first order term of the refractive index. However, non-linear terms such as the second and third order ionospheric terms and ray path bending effects are not considered in occultation measurements so far. Analysing selected CHAMP–GPS occultation events different higher order ionospheric terms are estimated and their effects on dual-frequency range estimation and total electron content (TEC) estimation are discussed. We have found that the separation between the GPS L1 and L2 ray paths exceeds the kilometer level during occultation for a vertical TEC level of more than 160 TEC units. Corresponding errors in the GPS dual-frequency range estimation and TEC estimation are found to exceed the meter and 10 TEC units level, respectively.     

Electron density retrieval from occulting GNSS signals using a gradient-aided inversion technique

In the coming years, opportunities for remote sensing of electron density in the Earth’s ionosphere will expand with the advent of Galileo, which will become part of the global navigation satellite system (GNSS). Methods for accurate electron density retrieval from radio occultation data continue to improve. We describe a new method of electron density retrieval using total electron content measurements obtained in low Earth orbit. This method can be applied to data from dual-frequency receivers tracking the GPS or Galileo transmitters. This simulation study demonstrates that the method significantly improves retrieval accuracy compared to the standard Abel inversion approach that assumes a spherically symmetric ionosphere. Our method incorporates horizontal gradient information available from global maps of Total Electron Content (TEC), which are available from the International GNSS Service (IGS) on a routine basis. The combination of ground and space measurements allows us to improve the accuracy of electron density profiles near the occultation tangent point in the E and F regions of the ionosphere.     

Validation of FORMOSAT-3/COSMIC radio occultation electron density profiles by incoherent scatter radar data

Low Earth Orbiting satellites carrying a dual frequency GPS receiver onboard offer a unique opportunity to remote sensing of the global ionosphere on a continuous basis. No other profiling technique unifies profiling through the entire F2-layer with global coverage. The FORMOSAT-3/COSMIC data can make a positive impact on the global ionosphere study providing essential information about the height electron density distribution and particularly over regions that are not accessible with ground-based measuring instruments such as ionosondes and GPS dual frequency receivers. Therefore, it is important to verify occultation profiles with other techniques and to obtain experience in the reliability of their derivation. In the given study we present results of comparison of the electron density profiles derived from radio occultation measurements on-board FS-3/COSMIC and from the Kharkov incoherent scatter radar sounding.     

Volcanism on Venus: Connecting link?

New observational data obtained from Venera-11 through Venera-14 and Pioneer-Venus data (1979–1982) call for searching for the interdependences of the phenomena occurring on this planet. It is proposed that the presence of electrical discharges in the near-surface layers of the atmosphere, the variable density of the above-cloud submicron haze, and the peculiarities of the cloud microphysics can be tied by a hypothesis that volcanic eruptions occur on Venus.     

An introduction to China FY3 radio occultation mission and its measurement simulation

GNSS (Global Navigation Satellite System) radio occultation mission for remote sensing of the Earth’s atmosphere will be performed by GNOS (GNSS Occultation Sounder) instrument on China FengYun-3 (FY3) 02 series satellites, the first of which FY3-C will be launched in the year 2013. This paper describes the FY3 GNOS mission and presents some results of measurement simulation. The key designed specifications of GNOS are also shown. The main objective of simulation is to provide scientific support for GNOS occultation mission on the FY3-C satellites. We used EGOPS software to simulate occultation measurements according to GNOS designed parameters. We analyzed the accuracy of retrieval profiles based on two typical occultation events occurring in China South–East area among total simulated events. Comparisons between the retrieval atmospheric profiles and background profiles show that GNOS occultation has high accuracy in the troposphere and lower stratosphere. The sensitivities of refractivity to three types of instrumental error, i.e. Doppler biases, clock stability and local multipath, were analyzed. The results indicated that the Doppler biases introduced by along-ray velocity error and GNOS clock error were the primary error sources for FY3-C occultation mission.     

Observation of wave structures in the upper atmosphere by means of radio holographic analysis of the radio occultation data

A radio holographic approach, developed by Pavelyev (1998), Hocke (1999), Igarashi (2000), is applied to observation of wave phenomena in the upper atmosphere using Global Positioning System — “Microlab-1” satellite (GPS/MET) radio occultation data. In the current state the radio holography approach uses the radar focused synthetic aperture principle to obtain high spatial resolution, and to remove the interference part corresponding to scattering from the upper ionosphere. High spatial resolution and accuracy of the radio halographic method is validated by means of revealing the weak signal reflected from the sea in the GPS/MET radio occultation data. The radio holographic method gives a new possibility to measure directly the vertical gradient of the electron density altitude profile in the D-layer using the radio occultation signal. The results of the application of radio holographic analysis to two GPS/MET occultation events (07 February 1997, No. 0447, 0158), in the D-region of the ionosphere, are discussed. Wave structures in the electron density concentration with a vertical spatial period of 1.4–6 km, and variations in the electron density gradient from ±5·10⁹ to ±8·10⁹ [1/(m³km)], have been retrieved from the D-layer data. The features observed in the vertical electron density profiles may be connected with breaking of gravity waves in the D-layer of the ionosphere.     

A global survey of COSMIC ionospheric peak electron density and its height: A comparison with ground-based ionosonde measurements

With a network of ground-based ionosondes distributed around the world, the ionospheric peak electron density and its height measured by FORMOSAT-3/COSMIC satellites in terms of GPS radio occultation technique are extensively examined in this article. It is found that, in spite of the latitude, the mean values of the peak electron density measured by COSMIC satellites are systematically smaller than those observed by ground-based ionosondes. The discrepancy between them is dependent on the latitude, namely, it is small in low and mid-latitudes and large in high-latitude region. Moreover, statistical analysis shows that the slopes of the regression line that is best fitted to the scatter diagram of occultation-retrieved peak electron density (ordinate axis) versus ionosonde-observed peak density (abscissa axis) are universally less than one. This feature is believed to be the result of path average effect of non-uniform distribution of the electron density along the GSP ray during the occultation. A comparison between COSMIC-measured peak height and ionosonde-derived peak height hmF2 indicates that the former is systematically higher than the latter. The difference in the two can be as large as 20% or more in equatorial and low-latitude regions. This result implies that the peak height hmF2 derived from the virtual height through true height analysis based on Titheridge method seems to underestimate the true peak height. The correlation between COSMIC and ionosonde peak electron densities is analyzed and the result reveals that correlation coefficient seems to be dependent on the fluctuation of the occultation-retrieved electron density profile. The correlation will be higher (lower) for the electron density profiles with smaller (larger) fluctuations. This feature suggests that the inhomogeneous distribution of the electron density along the GPS ray path during the occultation plays an important role affecting the correlation between COSMIC and ionosonde measurements.     

GRAS radio occultation on-board of Metop

The GRAS radio occultation instrument is flying on Metop-A and belongs to the EPS (EUMETSAT Polar System). GRAS observes GPS satellites in occultation. Within this work, validation of GRAS closed-loop bending angle data against co-located ECMWF profiles extracted from model fields and occultations from the COSMIC constellation of radio occultation instruments is shown. Results confirm the high data quality and robustness, where GRAS shows lower bending angle noise against ECMWF than COSMIC and in terms of occultations per day, one GRAS ≈ two COSMIC satellites. This is partly due to the operational setup of EPS. For the investigation we focus on two observation periods where updates in the ECMWF (March 2009) and COSMIC processing (October 2009) have improved the statistics further. Bending angles biases agree to within 0.5% against ECMWF and to within 0.1% against COSMIC after the updates for altitudes between 8 and 40 km. In addition, we also analyze the impact of the Metop orbit processing on the derived GRAS bending angle data, where different GPS and Metop orbit solutions are analyzed. Results show that a batch based orbit processing would improve in particular the bending angle bias behavior at higher altitudes. Requirements for the operational processing of GRAS data are briefly outlined, options to ease the use of other positioning system satellites in the near future are discussed. A simplified analysis on the observation of several of these systems, e.g. GPS and Galileo, from one platform shows that about 16% of occultations are found within 300 km, ±3 h, thus providing similar information. A constellation of 2 GRAS like instruments would have only about 10% close-by.     

Radio occultation experiments with Comet Halley

Previous radio occultation investigations on cometary comae and tails have included refraction measurements and intensity scintillations of natural radio sources used to derive the density and structure of the cometary plasma. Significant improvements in the coverage and sensitivity of these measurements will be achieved during the present apparition of Comet Halley. The comet missions GIOTTO and VEGA will also feature passive radio science expeirments designed to measure comet-induced Doppler shifts of the dual-frequency spacecraft signals during Halley flyby. A brief survey of these radio occultation measurement techniques and their application in the specific case of Comet Halley are presented.     

Assessment of precision in ionospheric electron density profiles retrieved by GPS radio occultations

The Constellation Observing System for Meteorology Ionosphere and Climate (COSMIC) is a six satellite radio occultation mission that was launched in April 2006. The close proximity of these satellites during some months after launch provides a unique opportunity to evaluate the precision of Global Positioning System (GPS) radio occultation (RO) retrievals of ionospheric electron density from nearly collocated and simultaneous observations. RO data from 30 consecutive days during July and August 2006 are divided into ten groups in terms of daytime or nighttime and latitude. In all cases, the best precision values (about 1%) are found at the F peak height and they slightly degrade upwards. For all daytime groups, it is seen that electron density profiles above about 120 km height exhibit a substantial improvement in precision. Nighttime groups are rather diverse: in particular, the precision becomes better than 10% above different levels between 120 and 200 km height. Our overall results show that up to 100–200 km (depending on each group), the uncertainty associated with the precision is in the order of the measured electron density values. Even worse, the retrieved values tend sometimes to be negative. Although we cannot rely directly on electron density values at these altitudes, the shape of the profiles could be indicative of some ionospheric features (e.g. waves and sporadic E layers). Above 200 km, the profiles of precision are qualitatively quite independent from daytime or latitude. From all the nearly collocated pairs studied, only 49 exhibited a difference between line of sight angles of both RO at the F peak height larger than 10°. After analyzing them we find no clear indications of a significant representativeness error in electron density profiles due to the spherical assumption above 120 km height. Differences in precision between setting and rising GPS RO may be attributed to the modification of the processing algorithms applied to rising cases during the initial period of the COSMIC mission.     

月球电离层掩星探测研究

日本SELENE/KAGUYA探测任务提供了研究月球电离层的机会。采用无线电掩星探测技术和趋势外推算法,消除地球电离层和行星际等离子体的干扰影响,残余的信号相位信息的变化反映了月球电离层的信息,估算出月球周围附近近似对称分布的稀薄电离层中电子总含量(TEC)约为每立方米10-14个。     

Total electron content measurements in ionospheric physics

With the advent of modern global networks of dual-frequency Global Positioning System (GPS), total electron content (TEC) measurements along slant paths connecting GPS receivers and satellites at 22,000 km have become the largest data set available to ionospheric scientists. The TEC can be calculated from the time and phase delay in the GPS signal using the GPS Toolkit, but an unknown bias will remain. In addition, UHF/VHF radio beacons on board low-Earth-orbiting satellites can also be used to measure the electron content. However, the TEC measurements are obtained by integrating TEC differences between slant paths, but also contain biases. It is often necessary to use data assimilative algorithms like the Ionospheric Data Assimilation Three-Dimensional (IDA3D), and to treat both GPS- and LEO-beacon TEC measurements as relative data in order to conduct ionospheric studies.     

A topside investigation over a mid-latitude digisonde station in Cyprus

We examine the systematic differences between topside electron density measurements recorded by different techniques over the low-middle latitude operating European station in Nicosia, Cyprus (geographical coordinates: 35.14^oN, 33.2^oE), (magnetic coordinates 31.86^oN, 111.83 ^oE). These techniques include space-based in-situ data by Langmuir probes on board.European Space Agency (ESA) Swarm satellites, radio occultation measurements on board low Earth orbit (LEO) satellites from the COSMIC/FORMOSAT-3 mission and ground-based extrapolated topside electron density profiles from manually scaled ionograms. The measurements are also compared with International Reference Ionosphere Model (IRI-2016) topside estimations and IRI-corrected NeQuick topside formulation (method proposed by Pezzopane and Pignalberi (2019)). The comparison of Swarm and COSMIC observations with digisonde and IRI estimations verifies that in the majority of cases digisonde underestimates while IRI overestimates Swarm observations but in general, IRI provides a better topside representation than the digisonde. For COSMIC and digisonde profiles matched at the F layer peak the digisonde systematically underestimates topside COSMIC electron density values and the relative difference between COSMIC and digisonde increases with altitude (above hmF2), while IRI overestimates the topside COSMIC electron density but after a certain altitude (~150 km above hmF2) this overestimation starts to decrease with altitude. The IRI-corrected NeQuick underestimates the majority of topside COSMIC electron density profiles and relative difference is lower up to approximately 100 km (above the hmF2) and then it increases. The overall performance of IRI-corrected NeQuick improves with respect to IRI and digisonde.     

An introduction of mountain-based GPS radio occultation experiments in China

The mountain-based GPS radio occultation is a novel approach to lower atmospheric profiling. The experiments of the mountain-based GPS radio occultation were conducted on the top of Mt. Yaogu (29.38°N, 113.68°E, ∼1240 m) on December 17, 2003, and on the top of Mt. Jiugong (29.39°N, 114.65°E, ∼1550 m) on July 24, 2004. Based on these observation data, the scientific data processing software has been developed and is used to retrieve successfully the atmospheric refractivity profiles. The validation experiment was performed on the top of Mt. Wuling (40.60°N, 117.48°E, ∼2118 m) during August 1–29, 2005. Collocated automatic weather station and the radiosondes nearby were operated simultaneously for the comparison campaign. Results show that the radio occultation technique obtained about 40 profiles every day with the receiver antenna pointing to the south. Comparisons show that the refractivity measured by occultation agree well with those by the radiosondes, but not well with those by the automatic weather station due to their much different geographic locations of measurements. Results of these experiments suggest that the mountain-based GPS radio occultation is an economic reliable novel technique monitoring temporal and spatial variations of local lower atmospheric environments.     

An alternative Shell inversion technique - Analysis and validation based on COSMIC and ionosonde data

Multi-channel Global Positioning System (GPS) carrier phase signals, received by the six low Earth orbiting (LEO) satellites from the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) program, were used to undertake active limb sounding of the Earth’s atmosphere and ionosphere via radio occultation. In the ionospheric radio occultation (IRO) data processing, the standard Shell inversion technique (SIT), transformed from the traditional Abel inversion technique (AIT), is widely used, and can retrieve good electron density profiles. In this paper, an alternative SIT method is proposed. The comparison between different inversion techniques will be discussed, taking advantage of the availability of COSMIC datasets. Moreover, the occultation results obtained from the SIT and alternative SIT at 500 km and 800 km, are compared with ionosonde measurements. The electron densities from the alternative SIT show excellent consistency to those from the SIT, with strong correlations over 0.996 and 0.999 at altitudes of 500 km and 800 km, respectively, and the peak electron densities (N_mF₂) from the alternative SIT are equivalent to the SIT, with 0.839 vs. 0.844, and 0.907 vs. 0.909 correlation coefficients when comparing to those by the ionosondes. These results show that: (1) the N_mF₂ and h_mF₂ retrieved from the SIT and alternative SIT are highly consistent, and in a good agreement with those measured by ionosondes, (2) no matter which inversion technique is used, the occultation results at the higher orbits (∼800 km) are better than those at the lower orbits (∼500 km).     

利用COSMIC掩星数据监测电离层的异常变化

分析了COSMIC掩星数据反演电子密度的方法,利用实例研究反演方法的特点,并采用ISR非相干散射雷达获取的电子密度数据进行验证,进而反演了长三角区域SHAO（IGS）站上空在日全食和太阳风暴期间的电子密度廓线图. 通过与平静日期间电离层电子密度进行比较,发现日全食及太阳风暴导致电离层发生的异常变化,从而提出COSMIC掩星数据反演电子密度在监测电离层变化时所具有的优势.      

Harmonic analysis of the ionospheric electron densities retrieved from FORMOSAT-3/COSMIC radio occultation measurements

In order to investigate the regular variations of the ionosphere, the least-squares harmonic estimation is applied to the time series of ionospheric electron densities derived from about five years of Global Positioning System radio occultation observations by FORMOSAT-3/COSMIC satellites. The analysis is done for different latitudes and altitudes in the region of Iran. The least-squares harmonic estimation is found to be a powerful tool for the frequency analysis of the completely unevenly spaced time series of radio occultation measurements. Although the obtained results are slightly different from the exact expected cycles (i.e. annual and diurnal components with their Fourier decompositions, and the 27-day period) due to the low horizontal resolution of radio occultation measurements, high vertical resolution of the observations enables us to detect not only the total electron content variations but also periodic patterns of electron densities at different altitudes of the ionosphere. The dominant diurnal and annual signals together with their Fourier series decompositions are obtained, which are consistent with the previous analyses on the total electron content. In the equatorial anomaly band, the annual component is weaker than its Fourier decomposition periods. In particular, the semiannual period dominates the annual component, indicating the relationship between the semiannual variation of the electron densities and the ionospheric equatorial anomaly. From detection of the phases of the components, it is revealed that the annual signal generally has its maximum value in summers at high altitudes, and in the winters at low altitudes. This is probably due to the higher [O/N₂] ratios in winter than in the summer in the lower ionosphere. Furthermore, the semiannual component mostly peaks around solstices or about a month before/after them.     

风云三号C星GNOS北斗掩星电离层探测初步结果

利用风云三号卫星C星GNOS掩星探测仪电离层数据,分析了2013年10月FY-3C GNOS探测的北斗掩星电离层廓线分布,将2013年10月1日至2015年10月10日期间FY-3C GNOS观测的F₂层峰值电子密度（N_mF₂）与地面电离层测高仪观测结果进行对比,验证了FY-3C GNOS北斗电离层掩星的探测精度.结果表明,FY3-C GNOS北斗电离层掩星与电离层测高仪探测的N_mF₂数据相关系数为0.96,平均偏差为10.21%,标准差为19.61%.在不同情况下其数据精度有如下特征:白天精度高于夜晚;夏季精度高于分季,分季精度高于冬季;中纬地区精度高于低纬地区,低纬地区精度高于高纬地区; BDS倾斜同步轨道（IGSO）卫星精度高于同步轨道（GEO）卫星和中轨道（MEO）卫星.FY-3C GNOS北斗电离层掩星与国际上其他掩星电离层数据精度的一致性对GNSS掩星探测资料的综合利用具有重大意义.      

Measurements of the upper troposphere and lower stratosphere during tropical cyclones using the GPS radio occultation technique

Water vapour transport to the upper troposphere and lower stratosphere by deep convective storms affects the radiation balance of the atmosphere and has been proposed as an important component of climate change. The aim of the work presented here is to understand if the GPS radio occultation technique is useful for characterization of this process. Our assessment addresses the question if severe storms leave a significant signature in radio occultation profiles in the upper troposphere/lower stratosphere. Radio occultation data from the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) were analyzed, focusing on two particular tropical cyclones with completely different characteristics, the hurricane Bertha, which formed in the Atlantic Basin during July 2008 and reached a maximum intensity of Category 3, and the typhoon Hondo, which formed in the south Indian Ocean during 2008 reaching a maximum intensity of Category 4. The result is positive, suggesting that the bending angle of a GPS radio occultation signal contains interesting information on the atmosphere around the tropopause, but not any information regarding the water vapour. The maximum percentage anomaly of bending angle between 14 and 18 km of altitude during tropical cyclones is typically larger than the annual mean by 5–15% and it can reach 20% for extreme cases. The results are discussed in connection to the GPS radio occultation receiver which will be part of the Atomic Clock Ensemble in Space (ACES) payload on the International Space Station.