Performance evaluation of IRI-2007 at equatorial latitudes and its Matlab version for GNSS applications

International Reference Ionosphere (IRI) model is the widely used empirical model for ionospheric predictions, especially TEC which is an important parameter for radio navigation and communication. The Fortran based IRI-2007 does not support real-time interactive visualization and debugging. Therefore, the source code is converted into Matlab and is validated for the purposes of this study. This facilitates easy representation of results and for near real-time implementation of IRI in the applications including spacecraft launching, now casting, pseudolite based navigation systems etc. In addition, the vertical delay results over the equatorial region derived from IRI and GPS data of three IGS stations namely Libreville (Garbon, Africa), Brasilia (Brazil, South America) and Hyderabad (India, Asia) are compared. As the IRI model does not account for plasmasphere TEC, the vertical delays are underestimated compared to vertical delays of GPS signals. Therefore, the model should be modified accordingly for precise TEC estimation.     

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.     

Evaluation of the improvement of IRI-2016 over IRI-2012 at the India low-latitude region during the ascending phase of cycle 24

This paper investigated the performance of the latest International Reference Ionosphere model (IRI-2016) over that of IRI-2012 in predicting total electron content (TEC) at three different stations in the Indian region. The data used were Global Positioning System (GPS) data collected during the ascending phase of solar cycle 24 over three low-latitude stations in India namely; Bangalore (13.02°N Geographic latitude, 77.57°E Geographic longitude), Hyderabad (17.25°N Geographic latitude, 78.30°E Geographic longitude) and Surat (21.16°N Geographic latitude, 72.78°E Geographic longitude). Monthly, the seasonal and annual variability of GPS-TEC have been compared with those derived from International Reference Ionosphere IRI-2016 and IRI-2012 with two different options of topside electron density: NeQuick and IRI01-corr. It is observed that both versions of IRI (i.e., IRI-2012 and IRI-2016) predict the GPS-TEC with some deviations, the latest version of IRI (IRI-2016) predicted the TEC similar to those predicted by IRI-2012 for all the seasons at all stations except for morning hours (0500 LT to 1000?LT). This shows that the effect of the updated version is seen only during morning hours and also that there is no change in TEC values by IRI-2016 from those predicted by IRI-2012 for the rest of the time of the day in the Indian low latitude region. The semiannual variations in the daytime maximum values of TEC are clearly observed from both GPS and model-derived TEC values with two peaks around March-April and September-October months of each year. Further, the Correlation of TEC derived by IRI-2016 and IRI-2012 with EUV and F10.7 shows similar results. This shows that the solar input to the IRI-2016 is similar to IRI 2012. There is no significant difference observed in TEC, bottom-side thickness (B0) and shape (B1) parameter predictions by both the versions of the IRI model. However, a clear improvement is visible in hmF2 and NmF2 predictions by IRI-2016 to that by IRI-2012. The SHU-2015 option of the IRI-2016 gives a better prediction of NmF2 for all the months at low latitude station Ahmedabad compare to AMTB 2013.     

Comparison of standard TEC models with a Neural Network based TEC model using multistation GPS TEC around the northern crest of Equatorial Ionization Anomaly in the Indian longitude sector during the low and moderate solar activity levels of the 24th solar cycle

The highest Total Electron Content (TEC) values in the world normally occur at Equatorial Ionization Anomaly (EIA) region resulting in largest ionospheric range delay values observed for any potential Space Based Augmentation System (SBAS). Reliable forecasting of TEC is crucial for satellite based navigation systems. The day to day variability of the location of the anomaly peak and its intensity is very large. This imposes severe limitations on the applicability of commonly used ionospheric models to the low latitude regions. It is necessary to generate a mathematical ionospheric forecasting and mapping model for TEC based on physical ionospheric influencing parameters. A model, IRPE-TEC, has been developed based on real time low latitude total electron content data using GPS measurements from a number of stations situated around the northern crest of the EIA during 2007 through 2011 to predict the vertical TEC values during the low and moderate solar activity levels of the 24th solar cycle. This model is compared with standard ionospheric models like International Reference Ionosphere (IRI) and Parameterized Ionospheric Model (PIM) to establish its applicability in the equatorial region for accurate predictions.     

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.     

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.     

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.     

A study on Nighttime Winter Anomaly (NWA) and other related Mid-latitude Summer Nighttime Anomaly (MSNA) in the light of International Reference Ionosphere (IRI) – Model

The ionospheric Nighttime Winter Anomaly (NWA) is a feature observed in the Northern Hemisphere at the American and in the Southern Hemisphere at the Asian longitude sector under low solar activity conditions. Jakowski et al. (2015) analyzed ground-based GPS derived TEC and peak electron density data from radio occultation measurements on Formosat-3/COSMIC satellites and confirmed the persistence of the phenomenon. Further, they assumed that Mid-latitude Summer Nighttime Anomaly (MSNA) and related special anomalies such as the Weddell Sea Anomaly (WSA) and the Okhotsk Sea Anomaly (OSA) are closely related to the NWA via enhanced wind-induced uplifting of the ionosphere. The aim of this paper is to study the factors causing these anomalies and also to investigate if these anomalies are re-produced by IRI. The results show that IRI model does include the NWA effect, though at a different longitude and could be improved for better predictions. The IRI-2016 model does show WSA in TEC but not in NmF2. Further, the IRI-2016 model could clearly predict the OSA both in NmF2 and TEC.     

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.     

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.     

On the performance of IRI-2016 to predict the North-South asymmetry of the equatorial ionization anomaly around 73°E longitude

The ionospheric total electron content (TEC) in both northern and southern Equatorial anomaly regions are examined by using the Global Positioning System (GPS) based TEC measurements around 73°E Longitude in the Asian sector. The TEC contour charts obtained at SURAT (21.16°N; 72.78°E; 12.9°N Geomagnetic Lat.) and DGAR (7.27°S; 72.37°E; 15.3°S Geomagnetic Lat.) over 73°E longitude during a very low solar activity phase (2009) and a moderate solar activity (2012) phase are used in this study. The results show the existence of hemispheric asymmetry and the effects of solar activity on the EIA crest in occurrence time, location and strength. The results are also compared with the TEC derived by IRI-2016 Model and it is found that the North-South asymmetry at the EIA region is clearly depicted by IRI-2016 with some discrepancies (up to 20% in the northern hemisphere at SURAT and up to 40% in the southern hemisphere at DGAR station for June Solstice and up to 10% both for SURAT and DGAR for December Solstice). This discrepancy in the IRI-2016 model is found larger during the year 2012 than that during the solar minimum year 2009 at both the hemispheres. Further, an asymmetry index, (A_i) is determined to illustrate the North-South asymmetry observed in TEC at EIA crest. The seasonal, annual and solar flux dependence of this index are investigated during both solstices and compared with the TEC derived by IRI.     

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.     

Improved predictions of equatorial TEC for real time applications

It is important to use models developed specifically for the equatorial ionospheric estimation for real-time applications, particularly in Satellite Navigation. This work demonstrates a methodology for improved predictions of VTEC in real time using the model developed for the equatorial ionosphere by the authors. This work has been done using TEC data of the low solar activity period of 2005 obtained using dual frequency GPS receivers installed under the GAGAN project of ISRO. For the purpose, the model is first used in conjunction with Kriging technique. Improvement in accuracy is observed when compared with the estimations from the model alone using the measurements as true reference. Further improvement is obtained by Bayesian combination of these estimates with independent Neural Network based predictions. Statistical performance of improvement is provided. An improvement of ∼1 m in confidence level of estimation of VTEC is obtained.     

Global validation of IRI TEC for high and medium solar activity conditions

The IRI model offers a choice of options for the computation of the electron density profile and electron content (TEC). Recently new options for the topside electron density profile have been developed, which have a strong impact on TEC. Therefore it is important to test massively the IRI and the new options with experimental data. A large number of permanent stations record dual frequency GPS data from which it is possible to obtain TEC values. Thirty-one worldwide distributed stations have been selected to investigate the capabilities of the IRI to reproduce experimental TEC. Data for years 2000 (high solar activity) and 2004 (medium solar activity) have been analyzed computing modeled values with the IRI-2001 and the IRI-2007-NeQuick topside options. It is found that IRI-2007-NeQuick option generally improves the estimate of the slant TEC, especially in the case of high latitudes stations during high solar activity.     

国际参考电离层模型输出参数峰值高度性能

电离层峰值高度H_mF₂是描述电离层形态的重要参数之一，国际参考电离层模型IRI-2016中融入了大量电离层测高仪和无线电掩星探测数据，用以提升H_mF₂的预测精度.本文利用太阳活动低年（2007—2010年）气象、电离层和气候卫星联合观测系统COSMIC探测数据描述全球范围内COSMIC H_mF₂的三维形态变化，对比分析了IRI-2016与IRI-2012模型的预测结果，同时分析了IRI-2016模型输出H_mF₂的性能.结果表明，IRI模型在中高纬度地区的输出结果高于COSMIC反演结果，而赤道及低纬地区则大都偏低.与IRI-2012模型相比，IRI-2016模型的输出结果在夜间至凌晨时段呈现较为明显的纬向梯度变化且大部分区域输出值偏高，但在白天时段赤道附近区域的输出值大都偏低.上述结果为电离层IRI模型的完善提供了一定参考.      

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.     

IRI 2001/90 TEC predictions over a low latitude station

Total electron content (TEC) over Tucumán (26.9°S, 294.6°W) measured with Faraday technique during the high solar activity year 1982, is used to check IRI 2001 TEC predictions at the southern crest of the equatorial anomaly region. Comparisons with IRI 90 are also made. The results show that in general IRI overestimates TEC values around the daily minimum and underestimates it the remaining hours. Better predictions are obtained using ground ionosonde measurements as input coefficients in the IRI model. The results suggest that for hours of maximum TEC values the electron density profile is broader than that assumed by the model. The main reason for the disagreement would be the IRI shape of the electron density profile.     

Doppler utilised Kalman estimation (DUKE) of ionospheric delay for satellite navigation

The ionospheric delay experienced by the satellite navigation signals depends upon the Total Electron Content (TEC) and needs to be corrected. While the single frequency receivers always use parametric models to correct this delay, dual frequency receivers, when suffers a loss of lock of one of its signal, also has to resort to these models. Here, an alternative method, based on Doppler, surrogated by range rate variation, has been attempted to estimate the ionospheric delay using a Kalman filter. GPS data have been used for all visible satellites over four days selected around the equinox and solstice with nominal geomagnetic conditions and estimations done in continuous and calibrated modes. Results of continuous estimation, obtained for a mid latitude station, showed moderate accuracy while it was significantly better for the calibrated mode with no seasonal dependence. Estimations done for station within the extent of equatorial anomaly, has not only resulted in relative deterioration in performance, but also shown seasonal dependence. Compared with estimates of Klobuchar model, the Calibrated estimation showed superior performance, conspicuously in the mid latitude station. However, for the continuous mode, performance was at par with the model at higher latitudes but inferior to it in regions within the extent of the equatorial anomaly.     

Assessment of improvement of the IRI model over Ethiopia for the modeling of the variability of TEC during the period 2013–2016

This paper discusses the monthly and seasonal variation of the total electron content (TEC) and the improvement of performance of the IRI model in estimating TEC over Ethiopia during the solar maximum (2013–2016) phase employing as reference the GPS derived TEC data inferred from four GPS receivers installed in different regions of Ethiopia; Assosa (geog 10.05°N, 34.55°E, Geom. 7.01°N), Ambo (8.97°N, 37.86°E, Geom. 5.42°N), Nazret (8.57°N, 39.29°E, Geom. 4.81°N) and Arba Minch (6.06°N, 37.56°E, Geom. 2.62°N). The results reveal that, in the years 2013–2016, the highest peak GPS-derived diurnal VTEC is observed in the March equinox in 2015 over Arba Minch station. Moreover, both the arithmetic mean GPS-derived and modelled VTEC values, generally, show maximum and minimum values in the equinoctial and June solstice months, respectively in 2014–2015. However, in 2013, the minimum and maximum arithmetic mean GPS-derived values are observed in the March equinox and December solstice, respectively. The results also show that, even though overestimation of the modelled VTEC has been observed on most of the hours, all versions of the model are generally good to estimate both the monthly and seasonal diurnal hourly VTEC values, especially in the early morning hours (00:00–03:00?UT or 03:00–06:00?LT). However, it has also been shown that the IRI 2007 and IRI 2012 versions generally perform best in matching the diurnal GPS derived TEC values as compared to that of the IRI 2016 version. In addition, the IRI 2012 version with IRI2001 option for the topside electron density shows the highest overestimation of the VTEC as compared to the other options. None of the versions of the IRI model are proved to be able to capture the effects of geomagnetic storms.