Evaluation of the IRI-2016 and NeQuick electron content specification by COSMIC GPS radio occultation,ground-based GPS and Jason-2 joint altimeter/GPS observations

We examined performance of two empirical profile-based ionospheric models, namely IRI-2016 and NeQuick-2, in electron content (EC) and total electron content (TEC) representation for different seasons and levels of solar activity. We derived and analyzed EC estimates in several representative altitudinal intervals for the ionosphere and the plasmasphere from the COSMIC GPS radio occultation, ground-based GPS and Jason-2 joint altimeter/GPS observations. It allows us to estimate a quantitative impact of the ionospheric electron density profiles formulation in several altitudinal intervals and to examine the source of the model-data discrepancies of the EC specification from the bottom-side ionosphere towards the GPS orbit altitudes. The most pronounced model-data differences were found at the low latitude region as related to the equatorial ionization anomaly appearance. Both the IRI-2016 and NeQuick-2 models tend to overestimate the daytime ionospheric EC and TEC at low latitudes during all seasons of low solar activity. On the contrary, during high solar activity the model results underestimated the EC/TEC observations at low latitudes. We found that both models underestimated the EC for the topside ionosphere and plasmasphere regions for all levels of solar activity. For low solar activity, the underestimated EC from the topside ionosphere and plasmasphere can compensate the overestimation of the ionospheric EC and, consequently, can slightly decrease the resulted model overestimation of the ground-based TEC. For high solar activity, the underestimated EC from the topside ionosphere and plasmasphere leads to a strengthening of the model underestimation of the ground-based TEC values. We demonstrated that the major source of the model-data discrepancies in the EC/TEC domain comes from the topside ionosphere/plasmasphere system.     

Thermosphere densities derived from Swarm GPS observations

After the detection of many anomalies in the Swarm accelerometer data, an alternative method has been developed to determine thermospheric densities for the three-satellite mission. Using a precise orbit determination approach, non-gravitational and aerodynamic-only accelerations are estimated from the high-quality Swarm GPS data. The GPS-derived non-gravitational accelerations serve as a baseline for the correction of the Swarm-C along-track accelerometer data. The aerodynamic accelerations are converted directly into thermospheric densities for all Swarm satellites, albeit at a much lower temporal resolution than the accelerometers would have been able to deliver. The resulting density and acceleration data sets are part of the European Space Agency Level 2 Swarm products.To improve the Swarm densities, two modifications have recently been added to our original processing scheme. They consist of a more refined handling of radiation pressure accelerations and the use of a high-fidelity satellite geometry and improved aerodynamic model. These modifications lead to a better agreement between estimated Swarm densities and NRLMSISE-00 model densities. The GPS-derived Swarm densities show variations due to solar and geomagnetic activity, as well as seasonal, latitudinal and diurnal variations. For low solar activity, however, the aerodynamic signal experienced by the Swarm satellites is very small, and therefore it is more difficult to accurately resolve latitudinal density variability using GPS data, especially for the higher-flying Swarm-B satellite. Therefore, mean orbit densities are also included in the Swarm density product.     

A method of estimating GPS instrumental biases with a convolution algorithm

This paper presents a method of deriving the instrumental differential code biases (DCBs) of GPS satellites and dual frequency receivers. Considering that the total electron content (TEC) varies smoothly over a small area, one ionospheric pierce point (IPP) and four more nearby IPPs were selected to build an equation with a convolution algorithm. In addition, unknown DCB parameters were arranged into a set of equations with GPS observations in a day unit by assuming that DCBs do not vary within a day. Then, the DCBs of satellites and receivers were determined by solving the equation set with the least-squares fitting technique. The performance of this method is examined by applying it to 361?days in 2014 using the observation data from 1311 GPS Earth Observation Network (GEONET) receivers. The result was crosswise-compared with the DCB estimated by the mesh method and the IONEX products from the Center for Orbit Determination in Europe (CODE). The DCB values derived by this method agree with those of the mesh method and the CODE products, with biases of 0.091?ns and 0.321?ns, respectively. The convolution method's accuracy and stability were quite good and showed improvements over the mesh method.     

Real-time clock comparison and monitoring with multi-GNSS precise point positioning: GPS,GLONASS and Galileo

The precise point positioning (PPP) technique is widely used in time and frequency applications. Because of the real-time service (RTS) project of the International GNSS Service, we can use the PPP technique for real-time clock comparison and monitoring. As a participant in the RTS, the Centre National d’Etudes Spatiales (CNES) implements the PPPWIZARD (Precise Point Positioning with Integer and Zero-difference Ambiguity Resolution Demonstrator) project to validate carrier phase ambiguity resolution. Unlike the Integer-PPP (IPPP) of the CNES, fixing ambiguities in the post-processing mode, the PPPWIZARD operates in the real-time mode, which is also called real-time IPPP (RT-IPPP). This paper focuses on applying the RT-IPPP for real-time clock comparison and monitoring. We review the principle of real-time clock comparison and monitoring, and introduce the methodology of the RT-IPPP technique. The observations of GPS, GLONASS and Galileo were processed for the experiments. Five processing modes were provided in the experiment to analyze the benefits of ambiguity resolution and multi-GNSS. In the clock comparison experiment, the average reduction ratios of standard deviations with respect to the G PPP mode range from 9.7% to 35.0%. In the clock monitoring experiment, G PPP mode can detect clock jumps whose magnitudes are larger than 0.9 ns. The RT-IPPP technique with GRE PPP AR (G) mode allows for the detection of any clock jumps larger than 0.6 ns. For frequency monitoring, G PPP mode allows detection of frequency changes larger than 1.1 × 10⁻¹⁴. When the RT-IPPP technique is applied, monitoring with GRE PPP AR (G) mode can detect frequency changes larger than 6.1 × 10⁻¹⁵.     

Comparison of GPS derived TEC with the TEC predicted by IRI 2012 model in the southern Equatorial Ionization Anomaly crest within the Eastern Africa region

We have compared the TEC obtained from the IRI-2012 model with the GPS derived TEC data recorded within southern crest of the EIA in the Eastern Africa region using the monthly means of the 5 international quiet days for equinoxes and solstices months for the period of 2012 – 2013. GPS-derived TEC data have been obtained from the Africa array and IGS network of ground based dual-frequency GPS receivers from four stations (Kigali (1.95°S, 30.09°E; Geom. Lat. 11.63°S), Malindi (2.99°S, 40.19°E; Geom. Lat. 12.42°S), Mbarara (0.60°S, 30.74°E; Geom. Lat. 10.22°S) and Nairobi (1.22°S, 36.89°E; Geom. Lat. 10.69°S)) located within the EIA crest in this region. All the three options for topside Ne of IRI-2012 model and ABT-2009 for bottomside thickness have been used to compute the IRI TEC. Also URSI coefficients were considered in this study. These results are compared with the TEC estimated from GPS measurements. Correlation Coefficients between the two sets of data, the Root-Mean Square Errors (RMSE) of the IRI-TEC from the GPS-TEC, and the percentage RMSE of the IRI-TEC from the GPS-TEC have been computed. Our general results show that IRI-2012 model with all three options overestimates the GPS-TEC for all seasons and at all stations, and IRI-2001 overestimates GPS-TEC more compared with other options. IRI-Neq and IRI-01-corr are closely matching in most of the time. The observation also shows that, GPS TEC are underestimated by TEC from IRI model during noon hours, especially during equinoctial months. Further, GPS-TEC values and IRI-TEC values using all the three topside Ne options show very good correlation (above 0.8). On the other hand, the TEC using IRI-Neq and IRI-01- corr had smaller deviations from the GPS-TEC compared to the IRI-2001.     

Analysis of the data processing strategies of spherical harmonic expansion model on global ionosphere mapping for moderate solar activity

Spherical harmonic (SH) expansion is widely used to model the global ionosphere map (GIM) of vertical total electron content (VTEC). According to the impact of different data processing methods of the SH expansion model on the VTEC maps, we specifically performed comprehensive analysis in terms of the data sampling rate, the time resolution, the spherical harmonic degree, and the relative constraint. One month of GPS data (January in 2016) from the International GNSS (Global Navigation Satellite System) Service (IGS) network in a moderate ionospheric activity period at the descending phase of Solar Cycle 24 was processed. To improve the computational efficiency of the daily GIM generation, the data sampling rate of 5?min was recommended allowing the GIM precision loss within 0.10 TECU (total electron content unit). The global VTEC map could be better represented in temporal and spatial domains with higher time resolution and higher spherical harmonic degree, especially at low latitude bands and in the southern hemisphere. The GIM precision improvement was about 10.91% for 1-h and about 15.15% for 0.5-h compared with the commonly used 2-h time resolution. The use of spherical harmonic degree 17 or 20 instead of 15 could improve the precision by 3.19% or 6.06%. We also found that an optimal relative constraint had to be found experimentally considering both the GIM precision and the GIM root mean square (RMS) map.     

TEC differences for the mid-latitude ionosphere in both sides of the longitudes with zero declination

Based on measurements of ground-based GPS station network, differences of the mid-latitude ionospheric TEC in the east and west sides of North America, South America and Oceania have been analyzed in this paper. Results show that for nearly all seasons from 2001 to 2010 and in both sides of the longitudes with zero declination, there exist systematic differences for the mid-latitude ionospheric TEC in the regions mentioned above and the features of these differences markedly depend upon the local time but less depend upon seasons and the level of solar activity. Theory analysis shows that the longitude variations of both declination and zonal thermospheric winds are one of important factors to cause differences of the mid-latitude ionospheric TEC in both sides of the longitudes with zero declination.     

利用GPS进行姿态估计的一种算法

首先建立了全球定位系统（ＧＰＳ）姿态确定的观测方程;然后给出了利用ＧＰＳ进行飞行器姿态估计的模型,并对该模型进行了线性处理;最后利用攻推广卡尔曼滤波技术,针对某飞行器进行了仿真计算。计算结果表明,对于不同的测量噪声和系统噪声,滤波器都有较好的估计,姿态估计的精度明显高于单纯ＧＰＳ姿态确定的精度,可以满足大多数飞行器对姿态确定的要求,证实了模型和算法可用性。     

INTEGRATED GPS/INS OF PSEUDO-RANGE DELTA-RANGE SEMI-PHYSICAL SIMULATION'S RESEARCH

In order to further study the performance of tightly integrated navigation system of GPS/INS, a semi-physical simulation of tightly coupled system has been done based on the data gathered from the experiment of integrated system of GPS and INS. The closed-loop Kalman Filter and U-D discompose algorithm have been used. The simulation results associated to four integrated models of pseudo-range, delta-range, pseudo-range and delta-range alternation, and pseudo-range and delta-range synthesis have been provided, and the actual effects of variously integrated models have been analyzed. The results show that the pseudo-range and delta-range synthesis coupled model is the most effective to improve the coupled system performance and the individual delta-range coupled model had better be avoided in application.     

GPS精密点定位的数据处理

在GPS精密点定位数据处理中,通常采用所谓无差分模型.然而,GPS单站观测量也可通过不同差分组合以减少未知参数数目,且在理论上与无差分模型等价.文中利用实际GPS观测数据,比较了不同点定位差分模型的定位精度和其他指标.