Higher order ionospheric error correction in BDS precise orbit determination

The ionospheric error affects the accuracy of the Global Navigation Satellite Systems observation and precise orbit determination. Usually, only the first order ionospheric error is considered, which can be eliminated by the ionospheric-free linear combination observation. But the remaining higher order ionospheric error will affect the accuracy of observations and their applications. In this paper, the influence of the higher order ionospheric error have been studied by using the International Geomagnetic Reference Field 13 and the Global Ionosphere Maps model produced by the Center for Orbit Determination in Europe. Focus on ionospheric error, the experiment of paper at doy 302 in 2019, which show that the second order ionospheric error impacting BeiDou Navigation Satellite System (BDS) B1I and B3I observation is 6.3569 mm and 11.8484 mm, respectively. Whereas, the third order ionospheric error impacting BDS B1I and B3I observation is 0.1734 mm and 0.3977 mm, respectively. Due to the current measurement accuracy of BDS carrier-phase observation can reach 2 mm, the influence of high order ionospheric error on observation should be considered. For BDS precise orbit determination, the orbit overlapping results are indicated that its orbit accuracy can be improved approximately 5 mm with the higher order ionospheric error correction, which is also in agreement with the results of Satellite Laser Ranging in this work.     

An update global model of hmF2 from values estimated from ionosonde and COSMIC/FORMOSAT-3 radio occultation

In this paper, we present our recent work on developing an updated global model of the ionospheric F2 peak height hmF2 parameter by combining data from the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC/FORMOSAT-3) radio occultation (RO) measurements and from the extended global ionosonde stations. In particular, 10 Chinese ionosonde stations’ data are newly introduced into this study. The modeling technique used is based on a two-layer empirical orthogonal function (EOF) expansion. Global distributions of hmF2 maps calculated using the newly constructed global model and the one provided by the International Reference Ionosphere model (IRI-ITU-R) are compared with the global distributions of hmF2 obtained by the COSMIC RO measurements and quantitative statistical analysis of the differences between the model results and those of the COSMIC RO measurements is made for the low (2008) and high (2012) solar activity years. The obtained average root-mean-square differences (RMSEs) for our model are 27.7 km (11.1%) and 31.0 km (9.8%), respectively for the years 2008 and 2012, whereas those for the IRI-ITU-R model are 39.9 km (16.9%) and 35.0 km (11.6%), respectively. Comparison of the results calculated both by our model and the IRI-ITU-R model with the digisonde observation is also made. The comparisons show that the newly constructed global hmF2 model can reproduce reasonably well the observations and perform better than IRI-ITU-R model.     

Global characteristics of ionospheric electric fields and disturbances during the first hours of magnetic storms

The ionospheric plasma density can be significantly disturbed during magnetic storms. In the conventional scenario of ionospheric storms, the negative storm phases with plasma density decreases are caused by neutral composition changes, and the positive storm phases with plasma density increases are often related to atmospheric gravity waves. However, recent studies show that the global redistribution of the ionospheric plasma is dominated primarily by electric fields during the first hours of magnetic storms. In this paper, we present the measurements of ionospheric disturbances by the DMSP satellites and GPS network during the magnetic storm on 6 April 2000. The DMSP measurements include the F region ion velocity and density at the altitude of ∼840 km, and the GPS receiver network provides total electron content (TEC) measurements. The storm-time ionospheric disturbances show the following characteristics. The plasma density is deeply depleted in a latitudinal range of ∼20° over the equatorial region in the evening sector, and the depletions represent plasma bubbles. The ionospheric plasma density at middle latitudes (20°–40° magnetic latitudes) is significantly increased. The dayside TEC is increased simultaneously over a large latitudinal range. An enhanced TEC band forms in the afternoon sector, goes through the cusp region, and enters the polar cap. All the observed ionospheric disturbances occur within 1–5 h from the storm sudden commencement. The observations suggest that penetration electric fields play a major role in the rapid generation of equatorial plasma bubbles and the simultaneous increases of the dayside TEC within the first 2 h during the storm main phase. The ionospheric disturbances at later times may be caused by the combination of penetration electric fields and neutral wind dynamo process.     

The occurrence of the mid-latitude ionospheric trough in GPS-TEC measurements

Simultaneous GPS observations from about 150 stations of European Permanent Network (EPN) have been used for studying dynamics of latitudinal profiles and structure of mid-latitude ionospheric trough (MIT). For the analyses, the TEC maps over Europe were created with high spatial and temporal resolution. The latitudinal profiles were produced from TEC maps with one-hour interval for geographic latitude range from 35N to 75N. The structure of latitudinal profiles relates to the occurrence of the ionospheric trough. The location of the trough depends on season, local time, and both geophysical and geomagnetic conditions. The trough structure in GPS-TEC demonstrates a smooth shape. The trough occurrence as a distinguished structure is more distinct during winter. The relation of TEC in the trough minimum to the equator and polar walls amounted to a factor of 2–4.     

A study of long-term climatology of ionospheric irregularities by using GPS phase fluctuations at the Brazilian longitudes

We have examined the ionospheric plasma irregularities that were recorded by using three ground-based receivers of the global positioning system (GPS) located at Brazilian longitudes during the period of a complete solar cycle, 1995–2005. The statistic results show that ionospheric irregularities are very easy to occur in December solstice months but rare to occur in June solstice months. Besides, the occurrence rates of irregularities in both December and June solstice months are little dependent on solar activity. However, in equinoctial months, the development of irregularities is obviously dependent on solar activity. There is a new finding in this study that if strong irregularities are distinguished from moderate ones, their occurrence rates would increase with solar activity during the December solstice months.     

Case studies of height structure of TID propagation characteristics using cross-correlation analysis of incoherent scatter radar and DPS-4 ionosonde data

The height structure of TID characteristics is studied on the base of the electron density profiles measured by two beams of the incoherent scatter radar and DPS-4 ionosonde. The height profiles of the TID propagation characteristics are obtained by means of cross-correlation and spectrum analysis of the radar and ionosonde data. The noticeable height variability of the TID parameters is observed. The variability is explained by interference of several TIDs. The obtained TID propagation characteristics are compared with known results of the TID studies.     

On the clustering of foF2 time series corresponding to disturbed ionospheric periods

This paper presents an analysis of a set of time series that represent foF2 disturbances during storm conditions, using clustering tools. The time series under study have been drawn from ionospheric observations obtained from eight European middle latitude ionosondes during a significant number of storm-time intervals and they are divided into eight groups according to the latitude (middle to low and middle to high) and the local time of the observation point at storm onset (afternoon, evening, morning, prenoon). The time series in each group have been analyzed using clustering-based methods. Specifically, each time series has been represented using two different ways of representation: the first is the raw representation while the second is through the parameters of the autoregressive (AR) model that best represents it. For each representation a hierarchy of clusterings is produced via the complete link algorithm. The two produced hierarchies are combined to a single one and the final clustering results are extracted from the produced hierarchy. The obtained results are in close agreement with the theoretical formulations concerning ionospheric storm effects at middle latitudes. In addition, they may be proved useful in the development of more accurate ionospheric forecasting methods.     

Modeling medium-scale TEC structures,observed by Belgian GPS receivers network

GALOCAD project “Development of a Galileo Local Component for the nowcasting and forecasting of atmospheric disturbances affecting the integrity of high precision Galileo applications” aims to perform a detailed study on ionospheric small- and medium-scale structures and to assess the influence of these structures on the reliability of Galileo precise positioning applications. GPS-derived TEC (total electron content) is obtained from the Belgium Dense Network (BDN), consisting of 67 permanent GPS stations. An empirical 3-D model is developed for studying these ionospheric structures. The model, named LLT model, described temporal variations of TEC in latitude/longitude frame (46°, 52°)N and (−1°, 11°)E. The spatial variations of TEC are modeled by Tchebishev base functions, while the temporal variations are described by a trigonometric basis. To fit the model to the data, the observed area is divided into bins with (1° × 1°) geographic scale and 6 min on time axis. LLT model is made flexible, with varying number of coefficients along each axis. This allows different degree of smoothing, which is the key element of the present approach. Model runs with higher number of coefficients, capturing in details medium-scale TEC structures are subtracted from results obtained with smaller number of coefficients; the latter represent the background ionosphere. The residual structures are localized and followed as they travel across the observed area. In this way, the size, velocity, and direction of the irregular structures are obtained.     

Seasonal variability of GPS-VTEC and model during low solar activity period (2006–2007) near the equatorial ionization anomaly crest location in Chinese zone

Variability of vertical TEC recorded at Fuzhou (26.1°N, 119.3°E, geomagnetic latitude 14.4°N), Xiamen (24.5°N, 118.1°E, geomagnetic latitude 13.2°N), Nanning (22.8°N, 108.3°E, geomagnetic latitude 11.4°N), China, during the low solar activity in 2006–2007 have been analyzed and discussed. Remarkable seasonal anomaly was found over three stations with the highest value during spring and the lowest value during summer. The relative standard deviation of VTEC is over 20% all the time, with steady and smooth variation during daytime while it has a large fluctuation during nighttime. The biggest correlation coefficient was found in the VTEC-sunspot pair with a value of over 0.5. It seems that solar activity has a better correlation ship than geomagnetic activity with the variation of VTEC and better correlations are found with more long-term data when comparing our previous study. The results of comparing observation with model prediction in three sites reveal again that the SPIM model overestimates the measured VTEC in the low latitude area.     

Global empirical models of the density peak height and of the equivalent scale height for quiet conditions

Monthly average electron density profiles have been calculated from hourly electron density N(h) recorded in 26 digisonde stations distributed worldwide encompassing the time interval 1998–2006. The ionospheric electron density peak height of the F2 region, hmF2, and the effective scale height at the hmF2, Hm, deduced from average profiles have been analyzed to obtain the quiet-time behavior and have been analytically modeled by the spherical harmonic analysis (SH) technique using the modip latitude as the coordinate of the reference system. The coefficients of the SH models of hmF2 and Hm are bounded to the solar activity, and the temporal and seasonal variations are considered by Fourier expansion of the coefficients. The SH models provide a tool to predict hmF2 and Hm located anywhere in the range of latitudes between of 70°N and 70°S and at any time. The SH analytical model for hmF2 improves the fit to the observations by 10% in average compared to the IRI prediction, and it might improve the IRI prediction of hmF2 by more than 30% at high and low latitudes. The analytical model for Hm predicts the quiet behavior of the effective scale height with accuracy better than 15% in average which enables to obtain a good estimation of vertical profiles. These results could be useful to estimate information for the topside profile formulation.