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⁻¹⁵.     

Identification of seismic activity sources on the subsatellite track by ionospheric plasma disturbances detected with the Sich-2 onboard probes

Using a cylindrical Langmuir probe and the authors’ proprietary two-channel pressure transducer, ionospheric plasma parameter distributions along the orbit of the Sich-2 satellite (Ukraine, 2011–2012) were measured. This paper is concerned with identifying the space–time location of ionospheric plasma disturbance sources, including the epicenters of actual earthquakes (before or during the satellite flyover) and incipient earthquakes on the subsatellite track, from the measured distributions of the electron density and temperature and the neutral particle temperature along the satellite orbit. To do this, the measured ionospheric plasma parameter distributions are connected to the coordinates on the subsatellite track.It is shown that local disturbances in the electron density and temperature and neutral particle temperature distributions in the satellite orbit in the ionosphere may serve as indicators of seismic activity on the subsatellite track. The epicenters of incipient earthquakes may be set off from other plasma parameter disturbance sources associated with seismic activity using information provided by special monitoring and survey centers that monitor the current seismic situation.     

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.     

Statistical study of the time delay of ionospheric TEC storms to geomagnetic storms in Taoyuan,Taiwan

We have studied the time delay of ionospheric storms to geomagnetic storms at a low latitude station Taoyuan (25.02°N, 121.21°E), Taiwan using the Dst and TEC data during 126 geomagnetic storms from the year 2002 to 2014. In addition to the known local time dependence of the time delay, the statistics show that the time delay has significant seasonal characteristics, which can be explained within the framework of the seasonal characteristics of the ionospheric TEC. The data also show that there is no correlation between the time delay and the intensity of magnetic storms. As for the solar activity dependence of the time delay, the results show that there is no relationship between the time delay of positive storms and the solar activity, whereas the time delay of negative storms has weakly negative dependence on the solar activity, with correlation coefficient −0.41. Especially, there are two kinds of extreme events: pre-storm response events and long-time delay events. All of the pre-storm response events occurred during 15–20 LT, manifesting the Equator Ionospheric Anomaly (EIA) feature at Taoyuan. Moreover, the common features of the pre-storm response events suggest the storm sudden commencement (SSC) and weak geomagnetic disturbance before the main phase onset (MPO) of magnetic storms are two main possible causes of the pre-storm response events. By analyzing the geomagnetic indices during the events with long-time delay, we infer that this kind of events may not be caused by magnetic storms, and they might belong to ionospheric Q-disturbances.     

A comparative study of ionospheric IRIEup and ISP assimilative models during some intense and severe geomagnetic storms

Three-dimensional (3-D) electron density matrices, computed in the Mediterranean area by the IRI climatological model and IRIEup and ISP nowcasting models, during some intense and severe geomagnetic-ionospheric storms, were ingested by the ray tracing software tool IONORT, to synthesize quasi-vertical ionograms. IRIEup model was run in different operational modes: (1) assimilating validated autoscaled electron density profiles only from a limited area which, in our case, is the Mediterranean sector (IRIEup_re(V) mode); (2) assimilating electron density profiles from a larger region including several stations spread across Europe: (a) without taking care of validating the autoscaled data in the assimilation process (IRIEup(NV)); (b) validating carefully the autoscaled electron density profiles before their assimilation (IRIEup(V)).The comparative analysis was carried out comparing IRI, IRIEup_re(V), ISP, IRIEup(NV), and IRIEup(V) foF2 synthesized values, with corresponding foF2 measurements autoscaled by ARTIST, and then validated, at the truth sites of Roquetes (40.80°N, 0.50°E, Spain), San Vito (40.60°N, 17.80°E, Italy), Athens (38.00°N, 23.50°E, Greece), and Nicosia, (35.03°N, 33.16°E, Cyprus). The outcomes demonstrate that: (1) IRIEup_re(V), performs better than ISP in the western Mediterranean (around Roquetes); (2) ISP performs slightly better than IRIEup_re(V) in the central part of Mediterranean (around Athens and San Vito); (3) ISP performance is better than the IRIEup_re(V) one in the eastern Mediterranean (around Nicosia); (4) IRIEup(NV) performance is worse than the IRIEup(V) one; (5) in the central Mediterranean area, IRIEup(V) performance is better than the IRIEup_re(V) one, and it is practically the same for the western and eastern sectors.Concerning the overall performance, nowcasting models proved to be considerably more reliable than the climatological IRI model to represent the ionosphere behaviour during geomagnetic-ionospheric storm conditions; ISP and IRIEup(V) provided the best performance, but neither of them has clearly prevailed over the other one.     

Characterization of ionospheric amplitude scintillations using wavelet entropy detrended GNSS data

The extensive monitoring networks of Global Navigation Satellite System (GNSS) ionospheric scintillation have been established to continuously log observation data. Further, the amplitude scintillation index and the phase scintillation index, which are derived from scintillation observations, are anticipated to accommodate the accuracy requirement of both the user level and the monitoring station level. However, raw scintillation observations essentially measure superposed waveform impairments of GNSS signals propagating through ionosphere and troposphere. It implies that fluctuations of raw scintillation observations are caused by multiple factors from the entire radio propagation environment. Hence, it is crucial to characterize ionospheric scintillations from GNSS observation data. And the characterization is implemented through extracting fluctuations of raw observations merely induced by ionospheric scintillations. Designed to address this problem by means of Fourier filtering detrending, the present work investigates the influence of varying detrending cutoff frequencies on wavelet statistical energy and wavelet entropy distributions of scintillation data. It consequently derives criteria on the optimum detrending cutoff frequency for three types of raw amplitude scintillation data, which are classified by their wavelet energy distributions. Results of the present work verify that detrending with specific optimum cutoff frequencies rather than the fixed and universally applicable one renders the validity and credibility of characterizing ionospheric scintillations as the part of GNSS observation fluctuations purely induced by ionosphere electron density irregularities whose scale sizes are comparable with or smaller than the Fresnel scale.     

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.     

L-band nighttime scintillations at the northern edge of the EIA along 95°E during the ascending half of the solar cycle 24

The characteristics of nighttime ionospheric scintillations measured at the L-band frequency of 1.575 GHz over Dibrugarh (27.5°N, 95°E, MLAT ～ 17°N, 43° dip) during the ascending half of the solar cycle 24 from 2010 to 2014 have been investigated and the results are presented in this paper. The measurement location is within or outside the zone of influence of the equatorial ionization anomaly depending on solar and geomagnetic activity. Maximum scintillation is observed in the equinoxes irrespective of solar activity with clear asymmetry between March and September. The occurrence frequency in the solstices shifts from minimum in the June solstice in low solar activity to a minimum in the December solstice in high solar activity years. A significant positive correlation of occurrence of scintillations in the June solstice with solar activity has been observed. However, earlier reports from the Indian zone (～75°E) indicate negative or no correlation of scintillation in June solstice with solar activity. Scintillations activity/occurrence in solstices indicates a clear positive correlation with Es recorded simultaneously by a collocated Ionosonde. In equinoxes, maximum scintillations occur in the pre-midnight hours while in solstices the occurrence frequency peaks just after sunset. The incidence of strong scintillations (S₄ ≥ 0.4) increases with increase in solar activity. Strong (S₄ ≥ 0.4) ionospheric scintillations accompanied by TEC depletions in the pre-midnight period is attributed to equatorial irregularities whereas the dusk period scintillations are related to the sporadic-E activity. Present results thus indicate that the current location at the northern edge of the EIA behaves as low as well as mid-latitude location.     

中纬电离层暴时形态的理论模式研究

结合观测结果对中纬电离层暴时形态进行理论模式研究.分析了两次电离层暴变事件中影响其基本形态的原子－分子含量比、上部输运通量和中性风等因素的行为,发现在这两次具有不同特点的事件中,中性风较为平稳,而原子－分子含量比和耦合输运通量的相对作用则各不相同.