Evaluation of the IRI-2007 model options for the topside electron density

The International Reference Ionosphere (IRI) 2007 provides two new options for the topside electron density profile: (a) a correction of the IRI-2001 model, and (b) the NeQuick topside formula. We use the large volume of Alouette 1, 2 and ISIS 1, 2 topside sounder data to evaluate these two new options with special emphasis on the uppermost topside where IRI-2001 showed the largest discrepancies. We will also study the accurate representation of profiles in the equatorial anomaly region where the profile function has to accommodate two latitudinal maxima (crests) at lower altitudes but only a single maximum (at the equator) higher up. In addition to IRI-2001 and the two new IRI-2007 options we also include the Intercosmos-based topside model of Triskova, Truhlik, and Smilauer [Triskova, L., Truhlik, V., Smilauer, J. An empirical topside electron density model for calculation of absolute ion densities in IRI. Adv. Space Res. 37 (5), 928–934, 2006] (TTS model) in our analysis. We find that overall IRI-2007-NeQ gives the best results but IRI-2007-corrected provides a more realistic representation of the altitudinal–latitudinal structure in the equatorial anomaly region. The applicability of the TTS model is limited by the fact that it is not normalized to the F2 peak density and height.     

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.     

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.     

On the use of NeQuick topside option in IRI-2007

The latest version of IRI includes various options for the computation of the topside electron density profile. One of the possible choices is based on NeQuick model. Its inclusion in IRI has been made transferring all the formulations used in NeQuick model. In details, an Epstein layer function is used to describe the electron density profile and the topside shape is controlled by an empirical parameter, connected to the NeQuick F2 bottomside thickness parameter, B2bot. It is computed also in this IRI topside option in order to maintain self-consistency with its original formulation. This paper analyses the possibility of using the IRI bottomside parameters for this option and its impact on the profile and TEC. The case of experimental peak values given as input is also analysed.     

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.     

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.     

Validate the IRI2007 model by the COSMIC slant TEC data during the extremely solar minimum of 2008

During 2008, the solar activity is extremely low. The satellite observations show that the ionospheric height and electron density is much lower than the predictions by the international reference ionosphere (IRI) model. In this paper, we compared the slant total electron content (TEC) observed by the COSMIC satellites during 2008 with the IRI model results. It is found that the IRI model with IRI2001 and IRI2001 Cor. topside options will always overestimate the electron density in both lower and higher altitudes. But the rest two topside options (NeQuick, and TTS) tend to overestimate the electron density in the F layer and underestimate it in the topside altitudes. The switch altitude between overestimation and underestimation and the latitude-local time distribution of the model deviation depend on the topside option. The current investigation might be useful for the model improvement as well as data assimilation work based on the IRI model and the LEO TEC data.     

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

Prediction of total electron content using the International Reference Ionosphere

The International Reference Ionosphere (IRI) is a model of the ionosphere, based on experimental data, which has been proposed as a standard ionospheric model. As such, it should be tested extensively to determine its range of validity. One of the ways in which the electron denisty profile given by the IRI, especially above the peak of the F layer, can be tested is to compare calculated and observed values of total electron content (TEC). We have therefore studied the discrepancies between calculated and observed values of TEC recorded at 15 stations covering a wide range of longitudes and latitudes, mainly in the northern hemisphere, and mainly for high levels of solar activity. W have found that the IRI produces reasonably accurate values of TEC at mid and high latitudes, but that it greatly underestimates the daytime values of TEC at low latitudes. We conclude therefore that the daytime electron density profile given by the IRI is reasonably accurate at mid and high latitudes, at least above the peak of the F2 layer. The situation at low latitudes clearly requires more work, and we have suggested two possible lines of study. The generally low discrepancies at night indicate that the night-time electron density profiles given by the IRI correspond fairly closely to the actual profiles.     

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 topside electron density computed by ionospheric models and plasma density observed by DMSP satellites

Electron density obtained by IRI (topside options NeQuick and IRI-Corr) and NeQuick models in their standard versions have been compared with plasma density values measured by F13 and F15 DMSP satellites for years of different solar activities. A statistical study of the differences between modeled and experimental data has been carried out to investigate each model performance.     

Comparing TOPEX TEC measurements with IRI predictions

TEC values obtained from TOPEX satellite were compared with the International Reference Ionosphere (IRI) 2001 model estimates. The present work also shows results of the IRI model with the option of a new topside electron density distribution (NeQuick model). TOPEX TEC measurements, which include years of high and middle to low solar activity (2000 and 2004), were analyzed by binning the region covered by the satellite (±66°) every five degrees of modip. In general, there is good agreement between IRI predictions and Topex measurements. Cases with large disagreements are observed at low and high latitudes during high solar activity. Comparing the model predictions using the default IRI2001 model and the NeQuick topside option show that the default IRI 2001 version represents the observed data in a more realistic way, but appears to be less reliable at high and low latitudes in some cases.     

Performance of the IRI-2007 model for equatorial topside ion density in the African sector for low and extremely low solar activity

The recent availability of new data sets during the recent extreme solar minimum provides an opportunity for testing the performance of the International Reference Ionosphere in historically undersampled regions. This study presents averages and variability of topside ionospheric densities over Africa as a function of season, local time, altitude, and magnetic dip latitude as measured by the Communication/Navigation Outage Forecast System (C/NOFS) satellite. The results are compared to the three topside model options available in IRI-2007. Overall, the NeQuick model is found to have the best performance, though during the deepest part of the solar minimum all three options significantly overestimate density.     

An analysis of the topside ionosphere parameters based on the long-duration Irkutsk incoherent scatter radar measurements

The topside ionosphere parameters are studied based on the long-duration Irkutsk incoherent scatter radar (52.9N, 103.3E) measurements conducted in September 2005, June and December 2007. As a topside ionosphere parameter we chose the vertical scale height (VSH) related to the gradient of the electron density logarithm above the peak height. For morphological studies we used median electron density profiles. Besides the median behavior we also studied VSH disturbances (deviations from median values) during the magnetic storm of September 11th 2005. We compared the Irkutsk incoherent scatter radar data with the Millstone Hill and Arecibo incoherent scatter radar observations, the IRI-2007 prediction (using the two topside options) and VSH derived from the Irkutsk DPS-4 Digisonde bottomside measurements.     

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.     

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.     

On the variabilities of the Total Electron Content (TEC) over the Indian low latitude sector

The Total Electron Content (TEC) from four locations in the Indian sector namely, Trivandrum (8.47°N, 76.91°E, Geomag.0.63°S, 0.3° dip), Waltair (17.7° N, 83.3°E, Geomag. 6.4°N, 20° dip), Bhopal (23.28°N, 77.34°E, Geomag.14.26°N, 33.2° dip), and Delhi (28.58°N, 77.21°E, Geomag.19.2°N, 43.4° dip) during a low sunspot year of 2004 are used to study the variabilities of the TEC. The day time TEC values are higher over Waltair and Bhopal compared to those at Trivandrum and Delhi. Considerable day-to-day variations in the diurnal values of TEC are observed at the anomaly crest locations. The observed GPS-TEC has been compared with the IRI-2007 model derived TEC considering three different options (IRI-2001, IRI-2001 corrected and Ne-Quick) available in the model for the topside electron density. The TEC derived with Ne-Quick and IRI-01 corrected options show better agreement with GPS-TEC while the TEC from IRI-01 method shows larger deviations. From the correlation analysis carried out between TEC value at 1300 h LT and solar indices parameters namely sunspot number (SSN), F10.7 and EUV, it is observed that the correlation is more during equinoctial months and less during summer months. The correlation coefficients observed over the anomaly locations, Bhopal and Delhi are lower compared to those at Trivandrum and Waltair.     

Comparison of ionospheric F2 peak parameters foF2 and hmF2 with IRI2001 at Hainan

Monthly median values of foF2, hmF2 and M(3000)F2 parameters, with quarter-hourly time interval resolution for the diurnal variation, obtained with DPS4 digisonde at Hainan (19.5°N, 109.1°E; Geomagnetic coordinates: 178.95°E, 8.1°N) are used to investigate the low-latitude ionospheric variations and comparisons with the International Reference Ionosphere (IRI) model predictions. The data used for the present study covers the period from February 2002 to April 2007, which is characterized by a wide range of solar activity, ranging from high solar activity (2002) to low solar activity (2007). The results show that (1) Generally, IRI predictions follow well the diurnal and seasonal variation patterns of the experimental values of foF2, especially in the summer of 2002. However, there are systematic deviation between experimental values and IRI predictions with either CCIR or URSI coefficients. Generally IRI model greatly underestimate the values of foF2 from about noon to sunrise of next day, especially in the afternoon, and slightly overestimate them from sunrise to about noon. It seems that there are bigger deviations between IRI Model predictions and the experimental observations for the moderate solar activity. (2) Generally the IRI-predicted hmF2 values using CCIR M(3000)F2 option shows a poor agreement with the experimental results, but there is a relatively good agreement in summer at low solar activity. The deviation between the IRI-predicted hmF2 using CCIR M(3000)F2 and observed hmF2 is bigger from noon to sunset and around sunrise especially at high solar activity. The occurrence time of hmF2 peak (about 1200 LT) of the IRI model predictions is earlier than that of observations (around 1500 LT). The agreement between the IRI hmF2 obtained with the measured M(3000)F2 and the observed hmF2 is very good except that IRI overestimates slightly hmF2 in the daytime in summer at high solar activity and underestimates it in the nighttime with lower values near sunrise at low solar activity.     

Mapping seasonal trends of electron temperature in the topside ionosphere based on DEMETER data

The diurnal, seasonal and latitudinal variations of the electron temperature in the Earth‘s topside ionosphere during relatively low solar activity period of 2005 – 2008 are investigated. In order to examine seasonal variations and morphology of the topside ionospheric plasma temperature, CNES micro-satellite DEMETER ISL data are used. Presented study is oriented on the dataset gathered in 2005 and 2008. Within conducted analysis, global maps of electron temperature for months of equinoxes and solstices have been developed. Furthermore, simultaneous studies on two-dimensional time series based on DEMETER measurements and predictions obtained with the IRI-2012 model supply examination of the topside ionosphere during recent deep solar minimum. Comparison with the IRI-2012 model reveals discrepancies between data and prediction, that are especially prominent during the periods of very low solar activity.