A morphological study of GPS-TEC data at Agra and their comparison with the IRI model

The temporal and seasonal variations of Total Electron Content (TEC) are studied at Agra (Geographic Lat. 27.17°N, Long. 78.89°E, Dip: 41.4°), India, which is in the equatorial anomaly region, for a period of 12 months from 01 January to 31 December, 2007 using a Global Positioning System (GPS) receiver. The mean TEC values show a minimum at 0500 h LT (LT = UT + 5.5 h) and a peak value at about 1400 h LT. The lowest TEC values are observed in winter whereas largest values are observed in equinox and summer. Anomalous variations are found during the period of magnetic disturbances. These results are compared with the TEC derived from IRI-2007 using three different options of topside electron density, NeQuick, IRI01-corr, and IRI-2001. A good agreement is found between the TEC obtained at Agra and those derived from IRI models.     

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.     

Low solar activity variability and IRI 2007 predictability of equatorial Africa GPS TEC

Diurnal, seasonal and latitudinal variations of Vertical Total Electron Content (VTEC) over the equatorial region of the African continent and a comparison with IRI-2007 derived TEC (IRI-TEC), using all three options (namely; NeQuick, IRI01-corr and IRI-2001), are presented in this paper. The variability and comparison are presented for 2009, a year of low solar activity, using data from thirteen Global Positioning System (GPS) receivers. VTEC values were grouped into four seasons namely March Equinox (February, March, April), June Solstice (May, June, July), September Equinox (August, September, October), and December Solstice (November, December, January). VTEC generally increases from 06h00 LT and reaches its maximum value at approximately 15h00–17h00 LT during all seasons and at all locations. The NeQuick and IRI01-corr options of the IRI model predict reasonably well the observed diurnal and seasonal variation patterns of VTEC values. However, the IRI-2001 option gave a relatively poor prediction when compared with the other options. The post-midnight and post-sunset deviations between modeled and observed VTEC could arise because NmF2 or the shape of the electron density profile, or both, are not well predicted by the model; hence some improvements are still required in order to obtain improved predictions of TEC over the equatorial region of the Africa sector.     

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%.     

Comparison of GPS-TEC measurements with NeQuick2 and IRI model predictions in the low latitude East African region during varying solar activity period (1998 and 2008–2015)

This paper examines the performances of NeQuick2, the latest available IRI-2016, IRI-2012 and IRI-2007 models in describing the monthly and seasonal mean total electron content (TEC) over the East African region. This is to gain insight into the success of the various model types and versions at characterizing the ionosphere within the equatorial ionization anomaly. TEC derived from five Global Positioning System (GPS) receivers installed at Addis Ababa (ADD, 5.33°N, 111.99°E Geog.), Asab (ASAB, 8.67°N, 116.44°E Geog.), Ambo (ABOO, 5.43°N, 111.05°E Geog.), Nairobi (RCMN, ?4.48°N, 108.46°E Geog.) and Nazret (NAZR, 4.78°N, 112.43°E Geog.), are compared with the corresponding values computed using those models during varying solar activity period (1998 and 2008–2015). We found that different models describe the equatorial and anomaly region ionosphere best depending on solar cycle, season and geomagnetic activity levels. Our results show that IRI-2016 is the best model (compared to others in terms of discrepancy range) in estimating the monthly mean GPS-TEC at NAZR, ADD and RCMN stations except at ADD during 2008 and 2012. It is also found that IRI-2012 is the best model in estimating the monthly mean TEC at ABOO station in 2014. IRI show better agreement with observations during June solstice for all the years studied at ADD except in 2012 where NeQuick2 better performs. At NAZR, NeQuick2 better performs in estimating seasonal mean GPS-TEC during 2011, while IRI models are best during 2008–2009. Both NeQuick2 and IRI models underestimate measured TEC for all the seasons at ADD in 2010 but overestimate at NAZR in 2009 and RCMN in 2008. The periodic variations of experimental and modeled TEC have been compared with solar and geomagnetic indices at ABOO and ASAB in 2014 and results indicate that the F10.7 and sunspot number as indices of solar activity seriously affects the TEC variations with periods of 16–32?days followed by the geomagnetic activity on shorter timescales (roughly periods of less than 16?days). In this case, NeQuick2 derived TEC shows better agreement with a long term period variations of GPS-TEC, while IRI-2016 and IRI-2007 show better agreement with observations during short term periodic variations. This indicates that the dependence of NeQuick2 derived TEC on F10.7 is seasonal. Hence, we suggest that representation of geomagnetic activity indices is required for better performance over the low latitude region.     

Validation of B2bot in the NeQuick model during high solar activity at Hainan station

NeQuick ionospheric electron density model, which has been developed to version 2, produces the full electron density profile in the ionosphere. Each part of the profile is modeled using Epstein layer formalism. Simple empirical relations are used to compute the thicknesses of each layer. In order to validate the B2_bot parameter in the NeQuick model during high solar activity, we use the data at Hainan, China (109.1°E, 19.5°N; Geomagnetic coordinates: 178.95°E, 8.1°N), measured with DPS-4, and study the diurnal and seasonal variations of B2_bot, ΔB2 (B2_best − B2_NeQuick 2) and the seasonal median values of B2_best/B2_NeQuick 2 at that region. The results show that, (1) The differences between B2_best and B2_NeQuick 2 have diurnal and seasonal variations. (2) The diurnal variations of B2_NeQuick 2 are smaller than those of B2_best. (3) Generally, except for early morning the experimental values are properly reproduced. (4) Generally, during morning the NeQuick model has an underestimation. The magnitude of underestimation varies with LT and season.     

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.     

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.     

Correction of B2bot for NeQuick during low solar activity at Hainan station

This work is a continuation of the previous article and it focuses on low solar activity and modeling effort. NeQuick model uses Epstein layer formalism to model each part of the profile. We study the diurnal and seasonal variations of B2_bot, ΔB2 (B2_best − B2_NeQuick2) and R (B2_best/B2_{NeQuick 2}) at Hainan station during low solar activity. The results show it is possible to improve the B2_bot parameter of the NeQuick model at that region during low solar activity. Then, we use a function ?(t) with LT in different seasons to correct the B2_bot formula of NeQuick 2. The correction shows that (1) By the correction formula, the B2_bot of NeQuick is improved. The maximum standard deviation is improved for 9 km. (2) The correction formula is more effective in summer than in equinox and winter and performs better during early morning hours than during the rest of the day.     

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.     

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.     

Behavior of the equivalent slab thickness over three European stations

The total electron content (TEC) derived from the global positioning system (GPS) and the F2-layer peak electron density obtained from Digisonde data have been used to study the diurnal, seasonal and solar activity variations of the ionospheric equivalent slab thickness (τ) over three European stations located at Pruhonice (50.0°N, 15.0°E), Ebro (40.8°N, 0.5°E) and El Arenosillo (37.1°N, 353.3°E). The diurnal variation of the τ is characterized by daytime values lower than nighttime ones for all seasons at low solar activity while daytime values larger than nighttime characterizes the diurnal variation for summer at high solar activity. A double peak is noticeable at dusk and at dawn, better expressed for winter at low solar activity. The seasonal variations of τ depend on local time and solar activity, the daytime values of τ increases from winter to summer whereas nighttime values of τ show the opposite. The effect of the solar activity on τ depends on local time and season, there being very sensitive for winter nighttime values of τ. The results of this study are compared with those presented by other authors.     

Comparison of GPS TEC variations with IRI-2007 TEC prediction at equatorial latitudes during a low solar activity (2009–2011) phase over the Kenyan region

This work presents an analysis of the Total Electron Content (TEC) derived from the International GNSS Service (IGS) receivers at Malindi (mal2: 2.9^oS, 40.1^oE, dip −26.813^o), Kasarani (rcmn: 36.89^oE, 1.2^oS, dip −23.970^o), Eldoret (moiu: 35.3^oE, 0.3^oN, dip −21.037^o) and GPS-SCINDA (36.8^oE, 1.3^oS, dip −24.117^o) receiver located in Nairobi for the period 2009–2011. The diurnal, monthly and seasonal variations of the GPS derived TEC (GPS-TEC) and effects of space weather on TEC are compared with TEC from the 2007 International Reference Ionosphere model (IRI-TEC) using the NeQuick option for the topside electron density. The diurnal peaks in GPS-TEC is maximum during equinoctial months (March, April, October) and in December and minimum in June solstice months (May, June, July). The variability in GPS-TEC is minimal in all seasons between 0:00 and 04:00 UT and maximum near noon between 10:00 and 14:00 UT. Significant variability in TEC at post sunset hours after 16:00 UT (19:00 LT) has been noted in all the seasons except in June solstice. The TEC variability of the post sunset hours is associated with the occurrence of the ionization anomaly crest which enhances nighttime TEC over this region. A comparison between the GPS-TEC and IRI-TEC indicates that both the model and observation depicts a similar trend in the monthly and seasonal variations. However seasonal averages show that IRI-TEC values are higher than the GPS-TEC. The IRI-TEC also depicts a double peak in diurnal values unlike the GPS-TEC. This overestimation which is primarily during daytime hours could be due to the model overestimation of the equatorial anomaly effect on levels of ionospheric ionization over the low latitude regions. The IRI-TEC also does not show any response to geomagnetic activity, despite the STORM option being selected in the model; the IRI model generally remains smooth and underestimates TEC during a storm. The GPS-TEC variability indicated by standard deviation seasonal averages has been presented as a basis for extending the IRI-model to accommodate TEC-variability.     

Seasonal dependence of TEC and SLAB thickness

Several years of hourly daily GPS measurements of the vertical total electron content (TEC) and of the equivalent slab thickness made at different European locations are analysed by using the linear regression technique to demonstrate the response of these two ionospheric parameters to seasonal variations. It is found that both TEC and slab thickness are highly correlated with season. Analytical relationships are determined expressing the seasonal dependence of the vertical TEC and of the equivalent slab thickness as a function of the seasonal parameter cos χ at noon in each location.     

IRI-2007预测TEC在广州地区的适用性分析

利用广州站(23.2°N, 113.3°E) GPS双频接收机监测的电离层TEC数据和IRI-2007模型不同电离层输入参数计算得到的TEC预测值, 对比分析了太阳活动低年(2008年)广州地区TEC的变化特征. 结果表明, TEC观测值周日变化在16:00LT左右达到最大值, 而IRI-TEC最大值出现时间较GPS-TEC提前1h左右. TEC季节变化在春秋分较高, 两至季节较低, 表现出明显的半年特性和季节依赖性, 并出现冬季异常现象. IRI-TEC与GPS-TEC在白天具有较好的一致性, 夜间偏差较大. 不同电离层输入参数得到的TEC预测值也相差较大, 选用顶部电子密度参数NeQuick、底部厚度参数B0 Table并用URSI系数计算F₂层峰值参数时, 能较好地反映TEC观测值的变化特征. 在对磁暴的响应上, 预测值无明显变化, 观测值则有比较明显的表现. 通过对比, 初步分析了利用IRI-2007模型预测TEC在广州地区的适用性, 并给出了合理的参数选择方案.      

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.     

Solar flux effects on equatorial ionization anomaly and total electron content over Brazil: Observational results versus IRI representations

The F layer critical frequency (foF2) as measured by Digisondes in the equatorial and low latitude locations in Brazil is analyzed to investigate the seasonal and solar flux controls of the intensity of the equatorial ionization anomaly (EIA) in the equinoctial month of March. The analysis also included the total electron content (TEC) as measured by a GPS receiver operated at the EIA crest location. The foF2 data set covered a period of large solar flux variation from 1996 to 2003, while the GPS TEC data was for a period in 2002–2003 when the solar flux parameter F10.7 underwent large variations, permitting in both cases an examination of the solar flux effects on these parameters. The seasonal variation pattern in TEC shows a maximum in equinoctial months and a minimum in June solstice, with similar variations for foF2. The solar flux dependence of the TEC is a maximum during equinoxes, especially for post-sunset TEC values at times when the latitudinal distribution is controlled by the equatorial evening plasma fountain processes. Significant variations with local time are found in the degree of solar flux dependence for both the TEC and EIA. The EIA intensity shows large dependence on F10.7 during post-sunset to midnight hours. These results are discussed in comparison with their corresponding IRI representations.     

Solar activity dependence of total electron content derived from GPS observations over Mbarara

Vertical total electron content (VTEC) observed at Mbarara (geographic co-ordinates: 0.60°S, 30.74°E; geomagnetic coordinates: 10.22°S, 102.36°E), Uganda, for the period 2001–2009 have been used to study the diurnal, seasonal and solar activity variations. The daily values of the 10.7 cm radio flux (F_10.7) and sunspot number (R) were used to represent Solar Extreme Ultraviolet Variability (EUV). VTEC is generally higher during high solar activity period for all the seasons and increases from 0600 h LT and reaches its maximum value within 1400 h–1500 h LT. All analysed linear and quadratic fits demonstrate positive VTEC-F_10.7 and positive VTEC-R correlation, with all fits at 0000 h and 1400 h LT being significant with a confidence level of 95% when both linear and quadratic models are used. All the fits at 0600 h LT are insignificant with a confidence level of 95%. Generally, over Mbarara, quadratic fit shows that VTEC saturates during all seasons for F_10.7 more than 200 units and R more than 150 units. The result of this study can be used to improve the International Reference Ionosphere (IRI) prediction of TEC around the equatorial region of the African sector.     

A prediction model of short-term ionospheric foF2 based on AdaBoost

In this paper, the AdaBoost-BP algorithm is used to construct a new model to predict the critical frequency of the ionospheric F2-layer (foF2) one hour ahead. Different indices were used to characterize ionospheric diurnal and seasonal variations and their dependence on solar and geomagnetic activity. These indices, together with the current observed foF2 value, were input into the prediction model and the foF2 value at one hour ahead was output. We analyzed twenty-two years’ foF2 data from nine ionosonde stations in the East-Asian sector in this work. The first eleven years’ data were used as a training dataset and the second eleven years’ data were used as a testing dataset. The results show that the performance of AdaBoost-BP is better than those of BP Neural Network (BPNN), Support Vector Regression (SVR) and the IRI model. For example, the AdaBoost-BP prediction absolute error of foF2 at Irkutsk station (a middle latitude station) is 0.32 MHz, which is better than 0.34 MHz from BPNN, 0.35 MHz from SVR and also significantly outperforms the IRI model whose absolute error is 0.64 MHz. Meanwhile, AdaBoost-BP prediction absolute error at Taipei station from the low latitude is 0.78 MHz, which is better than 0.81 MHz from BPNN, 0.81 MHz from SVR and 1.37 MHz from the IRI model. Finally, the variety characteristics of the AdaBoost-BP prediction error along with seasonal variation, solar activity and latitude variation were also discussed in the paper.     

Total electron content at low latitudes and its comparison with the IRI90

We used total electron content (TEC) data measured by Faraday rotation technique over Cachoeira Paulista (22.5°S, 45°W), in Brazil, to study the TEC variations with the solar flux at 10.7 cm (F10.7) and to compare the results with the IRI90 predictions. The data were divided into summer, equinox and winter. During the analysed period F10.7 varied from 66 up to 330. Our data shows that the observed TEC at 1600 LT (around the diurnal maximum) and at 0500 LT (around the diurnal minimum) increases with F10.7 until saturation is reached which occurs at F10.7≈210 to 220 for equinox and summer, and at F10.7≈180 for winter months. Comparison with the IRI90 predictions shows that IRI overestimates the TEC at 0500 LT for all solar flux values. At 1600 LT, IRI overestimates the observed TEC for low solar flux but underestimates it for high solar flux values.