Dynamics of global electron content in 1998–2005 derived from global GPS data and IRI modeling

We analyzed the dynamics of global electron content (GEC) for the period 1998–2005 and compared the estimated GEC with variations of the 10.7-cm solar radio emission and with and with GEC values obtained with IRI-2001. We found a strong resemblance between the curves’ shapes for the experimental and modeled GEC: strong semiannual variations are discernible in these series and both curves tend to increase the absolute GEC value during the period of maximum of solar activity. However, there are some significant distinctions, such as absence of 27-day fluctuations in the series of GEC computed by the IRI-2001. On the contrary, observational GEC reflects well dynamics of solar activity: 27-day variations of GEC are very similar to the ones of the index F10.7, but GEC undergoes a lagging of about of 30–60 h as compared to value of the F10.7 index. The relative amplitude of 27-day variations decreases from 8% at the rising and falling solar activity to 2% at the period of its maximum.     

Improvement of the real time total electron content based on the International Reference Ionosphere model

One of the methods to adapt the International Reference Ionosphere model to real time conditions is to use instantaneous values of the critical frequency of the ionosphere foF2. It is shown that there can be large discrepancies between model and experimental values of the total electron content TEC after this adaptation. Improvement can be provided by the use of an empirical model of the ionospheric slab thickness τ. This conclusion is based on analysis of contributions of various correction factors connected with foF2 and τ into discrepancies between model and experimental values of TEC.     

NeQuick 2 and IRI Plas VTEC predictions for low latitude and South American sector

Using vertical total electron content (VTEC) measurements obtained from GPS satellite signals the capability of the NeQuick 2 and IRI Plas models to predict VTEC over the low latitude and South American sector is analyzed. In the present work both models were used to calculate VTEC up to the height of GPS satellites. Also, comparisons between the performance of IRI Plas and IRI 2007 have been done. The data correspond to June solstice and September equinox 1999 (high solar activity) and they were obtained at nine stations. The considered latitude range extends from 18.4°N to ?64.7°N and the longitude ranges from 281.3°E to 295.9°E in the South American sector. The greatest discrepancies among model predictions and the measured VTEC are obtained at low latitudes stations placed in the equatorial anomaly region. Underestimations as strong as 40?TECU [1?TECU?=?10¹⁶?m^?2] can be observed at BOGT station for September equinox, when NeQuick2 model is used. The obtained results also show that: (a) for June solstice, in general the performance of IRI Plas for low latitude stations is better than that of NeQuick2 and, vice versa, for highest latitudes the performance of NeQuick2 is better than that of IRI Plas. For the stations TUCU and SANT both models have good performance; (b) for September equinox the performances of the models do not follow a clearly defined pattern as in the other season. However, it can be seen that for the region placed between the Northern peak and the valley of the equatorial anomaly, in general, the performance of IRI Plas is better than that of NeQuick2 for hours of maximum ionization. From TUCU to the South, the best TEC predictions are given by NeQuick2.The source of the observed deviations of the models has been explored in terms of CCIR foF2 determination in the available ionosonde stations in the region. Discrepancies can be also related to an unrealistic shape of the vertical electron density profile and or an erroneous prediction of the plasmaspheric contribution to the vertical total electron content. Moreover, the results of this study could be suggesting that in the case of NeQuick, the underestimation trend could be due to the lack of a proper plasmaspheric model in its topside representation. In contrast, the plasmaspheric model included in IRI, leads to clear overestimations of GPS derived TEC.     

Midday bottomside electron density profiles during moderate solar activity and comparison with IRI-2001

Bottomside electron density (Ne-h) profiles during midday (10–14 h) are analyzed using modern digital ionosonde observations at a low-middle latitude station, New Delhi (28.6N, 77.2E, dip 42.4N), for the period from January 2003 to December 2003, pertaining to moderate solar activity (MSA). Each individual profile is normalized with respect to the peak height and density (hmF2, NmF2) of the F2-region. These profiles are compared with those obtained from the International Reference Ionosphere (IRI-2001) model. Bilitza [Bilitza, D. International Reference Ionosphere 2000. Radio Sci. 36 (2), 261–275, 2001] using both the options namely: Gulyaeva’s model [Gulyaeva, T.L. Progress in ionospheric informatics based on electron density profile analysis of ionograms. Adv. Space Res. 7 (6) 39–48, 1987] and B0 Tab. option [Bilitza, D., Radicella, S.M., Reinisch, B.W., Adeniyi, J.O., Mosert Gonzalez, M.E., Zhang, S.R., Obrou, O. New B0 and B1 models for IRI. Adv. Space Res. 25 (1), 89–95, 2000]. The study reveals that during summer and equinox, the IRI model with B0 Tab. option in general, produces better agreement with the observed median profiles, while the IRI predictions using Gulyaeva’s option, overestimate the electron density distribution at all the heights below the F2-peak. However, during winter, in general, the IRI model, using both the options, reveals shows fairly good agreement with the observations.     

Spatial distribution of TEC across India in 2005: Seasonal asymmetries and IRI prediction

Total electron content measured simultaneously at 10 locations over India during the low solar activity year 2005 is used to examine the temporal and spatial asymmetries and also to assess the predictability of the International Reference Ionosphere in respect of the observed asymmetrical distribution. The stations are distributed in latitude along 77°E and in longitude along 23°N forming a meridional and a zonal chain respectively. A longitudinal gradient positive towards east was observed in the daytime hours of equinox and summer. Equinoctial asymmetry was prevalent across India during this year. Within the crest and equator, winter anomaly has been observed. It is found that IRI 2012 (with Ne Quick option, URSI coefficients) is unable to fully capture the temporal variation and spatial gradients of the ionization density in the Indian sector during 2005. The amount of offset between the model and measurement varies with local time and location.     

Ionospheric variability at taiwan low latitude station: Comparison between observations and IRI-2001 model

This paper presents the observed ionospheric F-region critical frequency, foF2, and peak height, hmF2, at northern crest of equatorial ionization anomaly (EIA) area station, namely Chung-Li (24.9°N, 121.1°E, dip 35°), and to be compared with International Reference Ionosphere model (IRI-2001) predictions for the period from 1994 to 1999, corresponding to half of the 23rd solar cycle. The diurnal and seasonal variation of foF2 and hmF2 are analyzed for different solar phases, respectively. The result shows the largest discrepancies were observed during nighttime for foF2 and hmF2, respectively. The value of foF2 both CCIR and URSI selected in the IRI model produced a good agreement during the daytime and underestimated during the noon time for high solar activities. The underestimation at noon time is mainly caused by the fountain effect from equator. Further, the peak height hmF2 shows a larger variability around the midnight than daytime in the equinox and winter seasons and reserved in summer, respectively. The study shows that the monthly median values of observed hmF2 is somewhat lower than those predicated by the IRI model, at night time in all the seasons except the period of 04:00–06:00 LT and reverse at daytime in summer. In general the IRI model predictions with respect to the observed in hmF2 is much better than foF2. The percentage deviation of the observed foF2 (hmF2) values with respect to the IRI model varies from 5% to 80% (0–25%) during nighttime and 2–17% (0–20%) at daytime, respectively. In general, the model generates good results, although some improvements are still necessary to implement in order to obtain better simulations for ionospheric low-latitudes region.     

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.     

Adaptation of the bottomside electron density profile of the IRI model to data of topside radiosounding

A method is proposed for reconstructing the electron density profiles N(h) of the IRI model from ionograms of topside satellite sounding of the ionosphere. An ionograms feature is the presence of traces of signal reflection from the Earth's surface. The profile reconstruction is carried out in two stages. At the first stage, the N(h) –profile is calculated from the lower boundary of the ionosphere to the satellite height (total profile) by the method presented in this paper using the ionogram. In this case, the monotonic profile of the topside ionosphere is calculated by the classical method. The profile of the inner ionosphere is represented by analytical functions, the parameters of which are calculated by optimization methods using traces of signal reflection, both from the topside ionosphere and from the Earth. At the second stage, the profile calculated from the ionogram is used to obtain the key parameters: the height of the maximum hmF2 of the F2 layer, the critical frequency foF2, the values of B0 and B1, which determine the profile shape in the F region in the IRI model. The input of key parameters, time of observation, and coordinates of sounding into the IRI model allows obtaining the IRI-profile corrected to real experimental conditions. The results of using the data of the ISIS-2 satellite show that the profiles calculated from the ionograms and the IRI profiles corrected from them are close to each other in the inner ionosphere and can differ significantly in the topside ionosphere. This indicates the possibility of obtaining a profile in the inner ionosphere close to the real distribution, which can significantly expand the information database useful for the IRTAM (IRI Realmax Assimilative Modeling) model. The calculated profiles can be used independently for local ionospheric research.     

Predicting indices of the ionosphere response to solar activity for the ascending phase of the 25th solar cycle

The international reference ionosphere, IRI, and its extension to plasmasphere, IRI-Plas, models require reliable prediction of solar and ionospheric proxy indices of solar activity for nowcasting and forecasting of the ionosphere parameters. It is shown that IRI prediction errors could increase for the F2 layer critical frequency foF2 and the peak height hmF2 due to erroneous predictions of the ionospheric global IG index and the international sunspot number SSN1 index on which IRI and IRI-Plas models are built. Regression relation is introduced to produce daily SSN1 proxy index from new time series SSN2 index provided from June 2015, after recalibration of sunspots data. To avoid extra errors of the ionosphere model a new solar activity prediction (SAP) model for the ascending part of the solar cycle SC25 is proposed which expresses analytically the SSN1 proxy index and the 10.7-cm radio flux F10.7 index in terms of the phase of the solar cycle, Φ. SAP model is based on monthly indices observed during the descending part of SC24 complemented with forecast of time and amplitude for SC25 peak. The strength of SC25 is predicted to be less than that of SC24 as shown with their amplitudes for eight types of indices driving IRI-Plas model.     

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.