Using GPS-TEC data to calibrate VTEC computed with the IRI model over Nigeria

This paper presents the development of a Total Electron Content (TEC) map for the Nigerian ionosphere. In this work, TEC measurements obtained from the AFRL-SCINDA GPS (Air Force Research Laboratory-Scintillation Network Decision Aid, Global Positioning System) equipment installed at Nsukka (6.87°N, 7.38°E) are used to adapt the International Reference Ionosphere (IRI) model for the Nigerian Ionosphere. The map is being developed as a computer program (implemented in the MATLAB programming language) that shows spatial and temporal representations of TEC for the Nigerian ionosphere. The method is aimed at showing how the IRI model can be used to estimate VTEC over wide areas by incorporating GPS measurements. This method is validated by using GPS VTEC data collected from a station in Ilorin (8.50°N, 4.55°E).     

Equatorial total electron content (TEC) at low and high solar activity

Total electron content (TEC) derived from ionosonde data recorded at the station of Korhogo (Lat = 9.33°N, Long = 5.43°W, Dip = 0.67°S) are compared to the International Reference Ionosphere (IRI) model predicted TEC for high (1999) and low (1994) solar activity conditions. The results show that the model represents the diurnal variation of the TEC as well as a solar activity and seasonal dependence. This variation is closer to that of the ionosonde-inferred TEC at high solar activity. However, at low solar activity the IRI overestimates the ionosonde-inferred TEC. The relative deviation ΔTEC is more prominent in the equinoctial seasons during nighttime hours where it is as high as 70%. At daytime hours, the relative deviation is estimated to 0–30%.     

Estimation of GPS instrumental biases from small scale network

With 4 GPS receivers located in the equatorial anomaly region in southeast China, this paper proposes a grid-based algorithm to determine the GPS satellites and receivers biases, and at the same time to derive the total electron content (TEC) with time resolution of 15 min and spatial resolution of 1° by 3.5° in latitude and longitude. By assuming that the TEC is identical at any point within a given grid block and the biases do not vary within a day, the algorithm arranges unknown biases and TECs with slant path TEC from the 4 receivers’ observations into a set of equations. Then the instrumental biases and the TECs are determined by using the least squares fitting technique. The performance of the method is examined by applying it to the GPS receiver chain observations selected from 16 geomagnetically quiet days in four seasons of 2006. It is found that the fitting agrees with the data very well, with goodness of fit ranging from 0.452 TECU to 1.914 TECU. Having a mean of 0.9 ns, the standard deviations for most of the GPS satellite biases are less than 1.0 ns for the 16 days. The GPS receiver biases are more stable than that of the GPS satellites. The standard deviation in the 4 receiver bias is from 0.370 ns to 0.855 ns, with a mean of 0.5 ns. Moreover, the instrumental biases are highly correlated with those derived from CODE and JPL with independent methods. The typical precision of the derived TEC is 5 TECU by a conservative estimation. These results indicate that the proposed algorithm is valid and qualified for small scale GPS network.     

Relative amplitude of the total electron content variations depending on geomagnetic activity

We developed a method of estimation of a relative amplitude dI/I of the total electron content (TEC) variations in the ionosphere as deduced from the data of the global GPS receivers network. To obtain statistically significant results we picked out three latitudinal belts provided in the Internet by the maximum number of GPS sites. They are high-latitudinal belt (50–80°N, 200–300°E; 59 sites), mid latitude belt (20–50°N, 200–300°E; 817 sites), and equatorial belt (±20°N, 0–360°E; 76 sites). The results of the analysis of the diurnal and latitudinal dependencies of dI/I and dI/I distribution probability for 52 days with different levels of geomagnetic activity are presented. It was found that on average the relative amplitude of the TEC variations varies within the range 0–10% proportionally to the value of the K_p geomagnetic index. In quiet conditions the relative amplitude dI/I of the TEC variations at night significantly exceeds the daytime relative amplitude. At high levels of magnetic field disturbances, the geomagnetic control of the amplitude of TEC variations at high and middle latitudes is much more significant than the regular diurnal variations. At the equatorial belt, on average, the amplitude of TEC variations in quiet and disturbed periods almost does not differ. The obtained results may be useful for development of the theory of ionospheric irregularities.     

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.     

Analysis of electron content variations over Japan during solar minimum: Observations and modeling

Comparative analysis of GPS TEC data and FORMOSAT-3/COSMIC radio occultation measurements was carried out for Japan region during period of the extremely prolonged solar minimum of cycle 23/24. COSMIC data for different seasons corresponded to equinox and solstices of the years 2007–2009 were analyzed. All selected electron density profiles were integrated up to the height of 700 km (altitude of COSMIC satellites), the monthly median estimates of Ionospheric Electron Content (IEC) were retrieved with use of spherical harmonics expansion. Monthly medians of TEC values were calculated from diurnal variations of GPS TEC estimates during considered month. Joint analysis of GPS TEC and COSMIC data allows us to extract and estimate electron content corresponded to the ionosphere (its bottom and topside parts) and the plasmasphere (h > 700 km) for different seasons of 2007–2009. Percentage contribution of ECpl to GPS TEC indicates the clear dependence from the time and varies from a minimum of about 25–50% during day-time to the value of 50–75% at night-time. Contribution of both bottom-side and topside IEC has minimal values during winter season in compare with summer season (for both day- and night-time). On average bottom-side IEC contributes about 5–10% of GPS TEC during night and about 20–27% during day-time. Topside IEC contributes about 15–20% of GPS TEC during night and about 35–40% during day-time. The obtained results were compared with TEC, IEC and ECpl estimates retrieved by Standard Plasmasphere–Ionosphere Model that has the plasmasphere extension up to 20,000 km (GPS orbit).     

TEC variability near northern EIA crest and comparison with IRI model

Monthly median values of hourly total electron content (TEC) is obtained with GPS at a station near northern anomaly crest, Rajkot (geog. 22.29°N, 70.74°E; geomag. 14.21°N, 144.9°E) to study the variability of low latitude ionospheric behavior during low solar activity period (April 2005 to March 2006). The TEC exhibit characteristic features like day-to-day variability, semiannual anomaly and noon bite out. The observed TEC is compared with latest International Reference Ionosphere (IRI) – 2007 model using options of topside electron density, NeQuick, IRI01-corr and IRI-2001 by using both URSI and CCIR coefficients. A good agreement of observed and predicted TEC is found during the daytime with underestimation at other times. The predicted TEC by NeQuick and IRI01-corr is closer to the observed TEC during the daytime whereas during nighttime and morning hours, IRI-2001 shows lesser discrepancy in all seasons by both URSI and CCIR coefficients.     

Highlights about the performances of storm-time TEC modelling techniques for low/equatorial and mid-latitude locations

A statistical evaluation of storm-time total electron content (TEC) modelling techniques over various latitudes of the African sector and surrounding areas is presented. The source of observational TEC data used in this study is the Global Navigation Satellite Systems (GNSS), specifically the Global Positioning Systems (GPS) receiver networks. For each selected receiver station, three different storm-time models based on empirical orthogonal functions (EOF) analysis, non-linear regression analysis (NLRA) and Artificial neural networks (ANN), were implemented. Storm-time GPS TEC data used for both development and validation of the models was selected based on the storm criterion of $\mathit{Dst}?-50$ nT or $K_p?4$ to take into account both coronal mass ejections (CMEs) and co-rotating interaction regions (CIRs) driven storms, respectively. To make an independent test of the models, storm periods considered for validation were excluded from datasets used during the implementation of the models and results are compared with observations, monthly median values, and International Reference Ionosphere (IRI-2016) predictions. Considering GPS TEC as reference, a statistical analysis performed over six storm periods reserved for validation revealed that ANN model is about 10%, 26%, and 58% more accurate than EOF, NLRA, and IRI models, respectively. It was further found that, EOF model performs 15%, and 44% better than NLRA, and IRI models, respectively, while NLRA is 25% better than IRI. On the other hand, results are also discussed referring to the background ionosphere represented by monthly median TEC (MM TEC) and statistics are provided. Moreover, strengths and weaknesses of each model are highlighted.     

GPS derived TEC and foF2 variability at an equatorial station and the performance of IRI-model

The ionosphere induces a time delay in transionospheric radio signals such as the Global Positioning System (GPS) signal. The Total Electron Content (TEC) is a key parameter in the mitigation of ionospheric effects on transionospheric signals. The delay in GPS signal induced by the ionosphere is proportional to TEC along the path from the GPS satellite to a receiver. The diurnal monthly and seasonal variations of ionospheric electron content were studied during the year 2010, a year of extreme solar minimum (F10.7 = 81 solar flux unit), with data from the GPS receiver and the Digisonde Portable Sounder (DPS) collocated at Ilorin (Geog. Lat. 8.50°N, Long. 4.50°E, dip −7.9°). The diurnal monthly variation shows steady increases in TEC and F2-layer critical frequency (foF2) from pre-dawn minimum to afternoon maximum and then decreases after sunset. TEC show significant seasonal variation during the daytime between 0900 and 1900 UT (LT = UT + 1 h) with a maximum during the March equinox (about 35 TECU) and minimum during the June solstice (about 24 TECU). The GPS-TEC and foF2 values reveal a weak seasonal anomaly and equinoctial asymmetry during the daytime. The variations observed find their explanations in the amount of solar radiation and neutral gas composition. The measured TEC and foF2 values were compared with last two versions of the International Reference Ionosphere (IRI-2007 and IRI-2012) model predictions using the NeQuick and CCIR (International Radio Consultative Committee) options respectively in the model. In general, the two models give foF2 close to the experimental values, whereas significant discrepancies are found in the predictions of TEC from the models especially during the daytime. The error in height dependent thickness parameter, daytime underestimation of equatorial drift and contributions of electrons from altitudes above 2000 km have been suggested as the possible causes.     

Comparison of GPS TEC measurements with IRI-2007 TEC prediction over the Kenyan region during the descending phase of solar cycle 23

This paper presents an analysis of the Total Electron Content (TEC) derived from the International GNSS Service receiver (formerly IGS) at Malindi (2.9°S, 40.1°E), Kenya for the periods 2004–2006 during the declining phase of solar cycle 23. The diurnal, monthly and seasonal variations of the TEC are compared with TEC from the latest International Reference Ionosphere model (IRI-2007). The GPS–TEC exhibits features such as an equatorial noon time dip, semi-annual variations, Equatorial Ionization Anomaly and day-to-day variability. The lowest GPS–TEC values are observed near the June solstice and September equinox whereas largest values are observed near the March equinox and December solstice. The mean GPS–TEC values show a minimum at 03:00 UT and a peak value at about 10:00 UT. These results are compared with the TEC derived from IRI-2007 using the NeQuick option for the topside electron density (IRI–TEC). Seasonal mean hourly averages show that IRI-2007 model TEC values are too high for all the seasons. The high prediction primarily occur during daytime hours till around midnight hours local time for all the seasons, with the highest percentage deviation in TEC of more 90% seen in September equinox and lowest percentage deviation in TEC of less than 20% seen in March equinox. Unlike the GPS–TEC, the IRI–TEC does not respond to geomagnetic storms and does overestimate TEC during the recovery phase of the storm. While the modeled and observed data do correlate so well, we note that IRI-2007 model is strongly overestimating the equatorial ion fountain effect during the descending phase of solar cycle, and this could be the reason for the very high TEC estimations.     

Use of varying shell heights derived from ionosonde data in calculating vertical total electron content (TEC) using GPS – New method

The dispersive nature of the ionosphere makes it possible to measure its total electron content (TEC). Thus Global Positioning System, which uses dual-frequency radio signals, is an ideal system to measure TEC. When data from an ionosonde situated in polar region was observed, the height of an approximated thin shell of electrons (shell height) used in GPS studies was seen not to be fixed but rather changing with time. Here we introduce a new method in which we included the varying shell heights derived from the ionosonde to map the slant total electron content from GPS to obtain a more precise vertical total electron content of the ionosphere contrary to some previous methods which used fixed shell heights. In this paper we also compared the ionosonde derived TEC with the GPS derived vertical TEC (vTEC) values. These GPS vTEC values were obtained from GPS slant TEC (sTEC) measurements using both fixed shell height and varying shell heights (from ionosonde measurements). For the polar regions, the varying shell height approach produced better results than the fixed shell height and compared to exponential function, Chapman function seems to be a better function to model the topside ionosphere.     

The solar eclipse and its associated ionospheric TEC variations over Indian stations on January 15, 2010

An annular solar eclipse occurred over the Indian subcontinent during the afternoon hours of January 15, 2010. This event was unique in the sense that solar activity was minimum and the eclipse period coincides with the peak ionization time at the Indian equatorial and low latitudes. The number of GPS receivers situated along the path of solar eclipse were used to investigate the response of total electron content (TEC) under the influence of this solar eclipse. These GPS receivers are part of the Indian Satellite Based Augmentation System (SBAS) named as ‘GAGAN’ (GPS Aided Geo Augmented Navigation) program. The eight GPS stations located over the wide range of longitudes allows us to differentiate between the various factors induced due to solar eclipse over the equatorial and low latitude ionosphere. The effect of the eclipse was detected in diurnal variations of TEC at all the stations along the eclipse path. The solar eclipse has altered the ionospheric behavior along its path by inducing atmospheric gravity waves, localized counter-electrojet and attenuation of solar radiation intensity. These three factors primarily control the production, loss and transport of plasma over the equatorial and low latitudes. The localized counter-electrojet had inhibited the equatorial ionization anomaly (EIA) in the longitude belt of 72°E–85°E. Thus, there was a negative deviation of the order of 20–40% at the equatorial anomaly stations lying in this ‘inhibited EIA region’. The negative deviation of only 10–20% is observed for the stations lying outside the ‘inhibited EIA region’. The pre-eclipse effect in the form of early morning enhancement of TEC associated with atmospheric gravity waves was also observed during this solar eclipse. More clear and distinctive spatial and temporal variations of TEC were detected along the individual satellite passes. It is also observed that TEC starts responding to the eclipse after 30 min from start of eclipse and the delay of the maximum TEC deviation from normal trend with respect to the maximum phase of the eclipse was close to one hour in the solar eclipse path.     

Dayside ionospheric uplift during strong geomagnetic storms as detected by the CHAMP,SAC-C,TOPEX and Jason-1 satellites

Ionosphere response to severe geomagnetic storms that occurred in 2001–2003 was analyzed using data of global ionosphere maps (GIM), altimeter data from the Jason-1 and TOPEX satellites, and data of GPS receivers on-board CHAMP and SAC-C satellites. This allowed us to study in detail ionosphere redistribution due to geomagnetic storms, dayside ionospheric uplift and overall dayside TEC increase. It is shown that after the interplanetary magnetic field turns southward and intensifies, the crests of the equatorial ionization anomaly (EIA) travel poleward and the TEC value within the EIA area increases significantly (up to ∼50%). GPS data from the SAC-C satellite show that during the main phase of geomagnetic storms TEC values above the altitude of 715 km are 2–3 times higher than during undisturbed conditions. These effects of dayside ionospheric uplift occur owing to the “super-fountain effect” and last few hours while the enhanced interplanetary electric field impinged on the magnetopause.     

Comparison of observed TEC values with IRI-2007 TEC and IRI-2007 TEC with optional foF2 measurements predictions at an equatorial region,Chumphon, Thailand

In this research, as part of working towards improving the IRI over equatorial region, the total electron content (TEC) derived from GPS measurements and IRI-2007 TEC predictions at Chumphon station (10.72°N, 99.37°E), Thailand, during 2004–2006 is analyzed. The seasonal variation of the IRI-2007 TEC predictions is compared with the TEC from the IRI-2007 TEC model with the option of the actual F2 plasma frequency (foF2) measurements as well as the TEC from the GPS and International GNSS service (IGS). The Chumphon station is located at the equatorial region and the low latitude of 3.22°N. For a declining phase of the solar cycle (2004–2006), the study shows that the IRI-2007 TEC underestimates the IRI-2007 TEC with the foF2 observation at the nighttime by about 5 TECU. The maximum differences are about 15 TECU during daytime and 5 TECU during nighttime. The overestimation is more evident at daytime than at nighttime. When compared in terms of the root-mean square error (RMSE), we find that the highest RMSE between GPS TEC and IRI 2007 TEC is 14.840 TECU at 1230 LT in 2004 and the lowest average between them is 1.318 TECU at 0630 LT in 2006. The noon bite-out phenomena are clearly seen in the IRI-2007 TEC with and without optional foF2 measurements, but not on the GPS TEC and IGS TEC. The IRI TEC with optional foF2 measurements gives the lowest RMSE values between IRI TEC predicted and TEC measurement. However, the TEC measurements (GPS TEC and IGS TEC) are more correct to use at Chumphon station.     

Ionospheric scintillations at Guilin detected by GPS ground-based and radio occultation observations

The occurrence of ionospheric scintillations with S4 ? 0.2 was studied using GPS measurements at Guilin, China (25.29°N, 110.33°E; geomagnetic: 15.04°N, 181.98°E), a station located near the northern crest of the equatorial anomaly. The results are presented for data collected from January 2009 to March 2010. The results show that nighttime amplitude scintillations only took place in February and March of the considered years, while daytime amplitude scintillations occurred in August and December of 2009. Nighttime amplitude scintillations, observed in the south of Guilin, always occurred with phase scintillations, TEC (Total Electron Content) depletions, and ROT (Rate Of change of TEC) fluctuations. However, TEC depletions and ROT fluctuations were weak during daytime amplitude scintillations, and daytime amplitude scintillations always took place simultaneously for most of the GPS satellites which appeared over Guilin in different azimuth directions. Ground-based GPS scintillation/TEC observations recorded at Guilin and signal-to-noise-ratio (SNR) measurements obtained from GPS-COSMIC radio occultation indicate that nighttime and daytime scintillations are very likely caused by ionospheric F region irregularities and sporadic E, respectively. Moreover, strong daytime amplitude scintillations may be associated with the plasma density enhancements in ionospheric E region caused by the Perseid and Geminid meteor shower activities.     

Seasonal variation of total electron content at crest of equatorial anomaly station during low solar activity conditions

The variability of total electron content (TEC) over the crest of equatorial anomaly station Bhopal has been studied during the low solar activity period (2005–2006) using global positioning system (GPS) data. Diurnal variation of TEC is studied for different seasons. Interesting features like the winter anomaly, semiannual anomaly and noon bite out in TEC have been reported. GPS derived TEC is then compared with International Reference Ionosphere (IRI) 2001 model and the difference between predictions and observation is being studied. Using the variability index we have also studied the TEC variability for different seasons and also during quiet and disturbed conditions. A higher variability is observed on quiet days as compared to disturbed days during daytime and nighttime hours.     

TEC derived from some GPS stations in Nigeria and comparison with the IRI and NeQuick models

Total electron content (TEC) measured simultaneously using Global Positioning System (GPS) ionospheric monitors installed at some locations in Nigeria during the year 2011 (Rz = 55.7) was used to study the diurnal, seasonal, and annual TEC variations. The TEC exhibits daytime maximum, seasonal variation and semiannual variations. Measured TEC were compared with those predicted by the improved versions of the International Reference Ionosphere (IRI) and NeQuick models. The models followed the diurnal and seasonal variation patterns of the observed values of TEC. However, IRI model produced better estimates of TEC than NeQuick at all locations.     

Towards a GPS-based TEC prediction model for Southern Africa with feed forward networks

In this paper, first results from a national Global Positioning System (GPS) based total electron content (TEC) prediction model over South Africa are presented. Data for 10 GPS receiver stations distributed through out the country were used to train a feed forward neural network (NN) over an interval of at most five years. In the NN training, validating and testing processes, five factors which are well known to influence TEC variability namely diurnal variation, seasonal variation, magnetic activity, solar activity and the geographic position of the GPS receivers were included in the NN model. The database consisted of 1-min data and therefore the NN model developed can be used to forecast TEC values 1 min in advance. Results from the NN national model (NM) were compared with hourly TEC values generated by the earlier developed NN single station models (SSMs) at Sutherland (32.38°S, 20.81°E) and Springbok (29.67°S, 17.88°E), to predict TEC variations over the Cape Town (33.95°S, 18.47°E) and Upington (28.41°S, 21.26°E) stations, respectively, during equinoxes and solstices. This revealed that, on average, the NM led to an improvement in TEC prediction accuracy compared to the SSMs for the considered testing periods.     

Application of neural networks to South African GPS TEC modelling

The propagation of radio signals in the Earth’s atmosphere is dominantly affected by the ionosphere due to its dispersive nature. Global Positioning System (GPS) data provides relevant information that leads to the derivation of total electron content (TEC) which can be considered as the ionosphere’s measure of ionisation. This paper presents part of a feasibility study for the development of a Neural Network (NN) based model for the prediction of South African GPS derived TEC. The South African GPS receiver network is operated and maintained by the Chief Directorate Surveys and Mapping (CDSM) in Cape Town, South Africa. Vertical total electron content (VTEC) was calculated for four GPS receiver stations using the Adjusted Spherical Harmonic (ASHA) model. Factors that influence TEC were then identified and used to derive input parameters for the NN. The well established factors used are seasonal variation, diurnal variation, solar activity and magnetic activity. Comparison of diurnal predicted TEC values from both the NN model and the International Reference Ionosphere (IRI-2001) with GPS TEC revealed that the IRI provides more accurate predictions than the NN model during the spring equinoxes. However, on average the NN model predicts GPS TEC more accurately than the IRI model over the GPS locations considered within South Africa.     

Short term variations of total electron content (TEC) fitting to a regional GPS network over the Kingdom of Saudi Arabia (KSA)

The ionosphere is a dispersive medium for radio waves with the refractive index which is a function of frequency and total electron content (TEC). TEC has a strong diurnal variation in addition to monthly, seasonal and solar cycle variations and small and large scale irregularities. Dual frequency GPS observations can be utilized to obtain TEC and investigate its spatial and temporal variations. We here studied short term TEC variations over the Kingdom of Saudi Arabia (KSA). A regional GPS network is formed consisting of 16 sites in and around KSA. GPS observations, acquired between 1st and 11th February 2009, were processed on a daily basis by using the Bernese v5.0 software and IGS final products. The geometry-free zero difference smoothed code observables were used to obtain two hour interval snapshots of TEC and their RMS errors at 0.5 × 0.5 degree grid nodes and regional ionosphere models in a spherical harmonics expansion to degree and order six. The equatorial ionized anomaly (EIA) is recovered in the south of 20°N from 08:00 to 12:00 UT. We found that day-by-day TEC variation is more stable than the night time variation.