Comparison of model and experimental ionospheric parameters at high latitudes

Modern use and study of the auroral region needs to attract a wider class of models for describing conditions of radio wave propagation in the ionosphere. In this paper the possibilities of the International Reference Ionosphere model, well-proven and widespread in the mid-latitudes, are investigated in the high latitude zone. Model and measured values of the critical frequency foF2 for two mid-latitude stations (Juliusruh and Goosebay) and four high-latitude ones (Loparskaya, Sodankyla, Sondrestrom, Thule) are compared. Deviations of medians, variations from day to day and solar activity trends seemed to be similar for both areas. This similarity is irrespective of the RZ12 index. Special attention is paid to the TEC parameter and its determination using 6 versions of models, a new version of the model IRI2010 (IRI-Plas) among them. It is shown that the IRI-Plas model significantly improves the definition of TEC in contrast to the versions of IRI2007 and the new model NeQuick. The use of the median of the experimental equivalent slab thickness, together with the current values of the TEC, increases by a factor of two the agreement between calculated and measured foF2 values as compared with the variations from day to day. This allows foF2 to be defined in near-real time.     

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

The impact of large solar events on the total electron content of the ionosphere at mid latitudes

Ionospheric disturbances associated with solar activity may occur via two basic mechanisms. The first is related to the direct impact on the ionosphere of EUV photons from a flare, and the second by prompt electric field penetration into the magnetosphere during geomagnetic storms. In this paper we examine the possibility that these two mechanisms may have an impact at mid latitudes by calculating the total electron content (TEC) from GPS stations in Mexico during several large X-ray flares. We have found that indeed large, complex flares, which are well located, may affect the mid latitude ionosphere. In fact, in the solar events of July 14, 2000 and April 2001 storms, ionospheric disturbances were observed to increase up to 138 and 150 TECu, respectively, due to the influence of EUV photons. Also, during the solar events of July 2000, April 2001, Halloween 2003, January 2005 and December 2006, there are large ionospheric disturbances (up to 393 TECu in the Halloween Storms), due to prompt penetration electric field, associated with CME producing geomagnetic storm.     

Computerized ionospheric tomography with the IRI model

Computerized ionospheric tomography (CIT) is a method to estimate ionospheric electron density distribution by using the global positioning system (GPS) signals recorded by the GPS receivers. Ionospheric electron density is a function of latitude, longitude, height and time. A general approach in CIT is to represent the ionosphere as a linear combination of basis functions. In this study, the model of the ionosphere is obtained from the IRI in latitude and height only. The goal is to determine the best representing basis function from the set of Squeezed Legendre polynomials, truncated Legendre polynomials, Haar Wavelets and singular value decomposition (SVD). The reconstruction algorithms used in this study can be listed as total least squares (TLS), regularized least squares, algebraic reconstruction technique (ART) and a hybrid algorithm where the reconstruction from the TLS algorithm is used as the initial estimate for the ART. The error performance of the reconstruction algorithms are compared with respect to the electron density generated by the IRI-2001 model. In the investigated scenario, the measurements are obtained from the IRI-2001 as the line integral of the electron density profiles, imitating the total electron content estimated from GPS measurements. It has been observed that the minimum error between the reconstructed and model ionospheres depends on both the reconstruction algorithm and the basis functions where the best results have been obtained for the basis functions from the model itself through SVD.     

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

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.     

新乡上空电离层总电子含量的经验模式

基于电离层总电子含量与太阳活动性之间的线性相关关系,利用在新乡长期测量所得的电离层总电子含量的资料,得出了总电子含量的一种经验模式。     

The effect of total solar eclipse of October 3, 2005, on the total electron content over Europe

GPS observations from EUREF permanent GPS network were used to observe the response of TEC (Total Electron Content) to the total solar eclipse on October 3, 2005, under quiet geomagnetic conditions of the daytime ionosphere. The effect of the eclipse was detected in diurnal variations and more distinctly in the variations of TEC along individual satellite passes. The trough-like variations with a gradual decrease and followed by an increase of TEC at the time of the eclipse were observed over a large region. The depression of TEC amounted to 3–4 TECU. The maximum depression was observed over all stations located at the maximum path of the solar eclipse. The delay of a minimum level of TEC with respect to the maximum phase of the eclipse was about 20–30 min.     

Evaluation of the STORM model storm-time corrections for middle latitude

This paper presents results from the Storm-Time Ionospheric Correction Model (STORM) validation for selected Northern and Southern Hemisphere middle latitude locations. The created database incorporated 65 strong-to-severe geomagnetic storms, which occurred within the period 1995–2007. This validation included data from some ionospheric stations (e.g., Pruhonice, El Arenosillo) that were not considered in the development or previous validations of the model. Hourly values of the F2 layer critical frequency, foF2, measured for 5–7 days during the main and recovery phases of each selected storm were compared with the predicted IRI 2007 foF2 with the STORM model option activated. To perform a detailed comparison between observed values, medians and predicted foF2 values the correlation coefficient, the root-mean-square error (RMSE), and the percentage improvement were calculated. Results of the comparative analysis show that the STORM model captures more effectively the negative phases of the summer ionospheric storms, while electron density enhancement during winter storms and the changeover of the different storm phases is reproduced with less accuracy. The STORM model corrections are less efficient for lower-middle latitudes and severe geomagnetic storms.     

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.     

Near real-time PPP-based monitoring of the ionosphere using dual-frequency GPS/BDS/Galileo data

Ionosphere delay is very important to GNSS observations, since it is one of the main error sources which have to be mitigated even eliminated in order to determine reliable and precise positions. The ionosphere is a dispersive medium to radio signal, so the value of the group delay or phase advance of GNSS radio signal depends on the signal frequency. Ground-based GNSS stations have been used for ionosphere monitoring and modeling for a long time. In this paper we will introduce a novel approach suitable for single-receiver operation based on the precise point positioning (PPP) technique. One of the main characteristic is that only carrier-phase observations are used to avoid particular effects of pseudorange observations. The technique consists of introducing ionosphere ambiguity parameters obtained from PPP filter into the geometry-free combination of observations to estimate ionospheric delays. Observational data from stations that are capable of tracking the GPS/BDS/GALILEO from the International GNSS Service (IGS) Multi-GNSS Experiments (MGEX) network are processed. For the purpose of performance validation, ionospheric delays series derived from the novel approach are compared with the global ionospheric map (GIM) from Ionospheric Associate Analysis Centers (IAACs). The results are encouraging and offer potential solutions to the near real-time ionosphere monitoring.     

GPS scintillation and TEC depletion near the northern crest of equatorial anomaly over South China

This study presents a statistical analysis of GPS L-band scintillation with data observed from July 2008 to March 2012 at the northern crest of equatorial anomaly stations in Guangzhou and Shenzhen of South China. The variations of the scintillation with local time, season, solar activity and duration of scintillation patches were investigated. The relationship between the scintillation and TEC depletion was also reported. Our results revealed that GPS scintillation occurred from 19:30 LT (pre-midnight) to 03:00 LT (post-midnight). During quiet solar activity years, the scintillation was only observed in pre-midnight hours of equinox months and patches durations were mostly less than 60 min. During high solar activity years, more scintillation occurred in the pre-midnight hours of equinox and winter months; and GPS scintillation started to occur in the post-midnight hours of summer and winter. The duration of scintillation patches extended to 180 min in high solar activity years. Solar activity had a larger effect to strong scintillations (S₄ > 0.6) than to weak scintillations (0.6 ? S₄ > 0.2). Strong scintillations were accompanied by TEC depletion especially in equinox months. We also discussed the relationship between TEC depletion and plasma bubble.     

Response of the EIA ionosphere to the 7-8 May 2005 geomagnetic storm

In this paper, response of low latitude ionosphere to a moderate geomagnetic storm of 7–8 May 2005 (SSC: 1920 UT on 7 May with Sym-H minimum, ∼−112 nT around 1600 UT on 8 May) has been investigated using the GPS measurements from a near EIA crest region, Rajkot (Geog. 22.29°N, 70.74°E, Geomag.14°), India. We found a decrease in total electron content (TEC) in 12 h after the onset of the storm, an increase during and after 6 h of Sym-H deep minimum with a decrease below its usual-day level on the second day during the recovery phase of the storm. On 8 May, an increase of TEC is observed after sunset and during post-midnight hours (maximum up to 170%) with the formation of ionospheric plasma bubbles followed by a nearly simultaneous onset of scintillations at L-band frequencies following the time of rapid decrease in Sym-H index (−30 nT/h around 1300 UT).     

Statistical ionospheric E layer properties measured with the Cyprus Digisonde and comparisons with IRI predictions

This paper reports the diurnal, seasonal, and long term variability of the E layer critical frequency (foE) and peak height (hmE) derived from Digisonde measurements from 2009 to 2016 at the low-middle latitude European station of Nicosia, Cyprus (geographical coordinates: 35°N, 33°E, geomagnetic lat. 29.38°N, I = 51.7°). Manually scaled monthly median values of foE and hmE are compared with IRI-2012 predictions with a view to assess the predictability of IRI. Results show that in general, IRI slightly overestimates foE values both at low and high solar activity. At low solar activity, overestimations are mostly limited to 0.25?MHz (equivalent electron density, 0.775?×?10³?el/m^?3) but can go as high as 0.5?MHz (equivalent electron density, 3.1?×?10³?el/m^?3, during noon) around equinox. In some months, underestimations, though sporadic in nature, up to 0.25?MHz are noted (mostly during sunrise and sunset). At high solar activity, a similar pattern of over-/underestimation is evident. During the entire period of study, over-/under estimations are mostly limited to 0.25?MHz. In very few cases, these exceed 0.25?MHz but are limited to 0.5?MHz. Analysis of hmE reveals that: (1) hmE remains almost constant during ±2 to ±4?h around local noon, (2) hmE values are higher in winter than in spring, summer and autumn, (3) there are two maxima near sunrise and sunset with a noontime minimum in between. During the entire period of study, significant differences between observed hmE and the IRI predictions have been noted. IRI fails to predict hmE and outputs a constant value of 110?km, which is higher than most of the observed values. Over- and under estimations range from 3 to 13?km and from 0 to 3?km respectively.     

Total electron content variations in equatorial anomaly region

Diurnal variations in the total electron content (TEC) at Makerere University (00°19′N, 32°40′E, Geo Dip −22°), Uganda, have been investigated using a NovAtel GSV400B GPS receiver for the year 2010. The highest TEC values occurred from 13h00 to 17h00 local time (LT) throughout the year, with the highest values being exhibited during equinoctial months. In addition, there was some correlation between this high TEC and the moderate storms that occurred in 2010. These high TEC values have been attributed to the solar EUV ionization coupled with the upward vertical ExB drift. Nighttime enhancements were also found to be seasonally dependant, attaining maximum values during equinoctial months. These results were also compared with modeled TEC values by the IRI-2007 model. The modeled values were in good agreement with the measured values except for these two points: (1) the model had a short-fall in predicting the nighttime enhancements; and (2) the model’s minimum TEC did not coincide with the measured minimum in most of the months. Observed TEC depletions were found to correlate with an increase in the S₄ index and have been identified as a manifestation of the plasma density depletions of the equatorial origin.     

Deformation of the ionosphere structure during the space weather events on October–November 2003

We present the spatial maps of the ionosphere–plasmasphere slab thickness τ (ratio of the vertical total electron content, TEC, to the F-region peak electron density, NmF2) during the intense ionospheric storms of October–November 2003. The model-assisted technology for estimate of the upper boundary of the ionosphere, hup, from the slab thickness components in the bottomside and topside ionosphere – eliminating the plasmasphere contribution of τ – is applied at latitudes 35° to 70°N and longitudes −10° to 40°E, from the data of 20 observatories of GPS-TEC and ionosonde networks, for selected days and hours of October and November 2003. The daily–hourly values of NmF2, hmF2 and TECgps are used as the constrained parameters for the International Reference Ionosphere extended to the plasmasphere, IRI-Plas, during the ionospheric quiet days, positive and negative storm phases for estimate of τ and hup. Good correlation has been found between the slab thickness and the upper boundary of the ionosphere for the intense ionospheric storms at October–November 2003. During the negative phase of the ionospheric storm, when the ionospheric plasma density is exhausted, the nighttime upper boundary of the ionosphere is greatly uplifted towards the magnetosphere tail, while the daytime upper boundary of the ionosphere is reduced below 500 km over the Earth.     

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.