Comparing the improved Di Giovanni/Radicella model with sounding-based electron density profiles and with the IRI model

The Di Giovanni/Radicella model (DGR) /1/ determines a bottom side electron densty profile alone from the set of routinely scaled ionogram parameters foE, foF1, foF2 and M(3000)F2 and the total electron content; the smoothed sunspot number R12 appears in the calculation. Present designations are DGR2/2/ and DRR3 /3/ [see Appendix]; they are valid in the northern hemisphere. DGR is compared with electron density profiles derived from ionograms obtained at Juliusruh (54.6°N, 13.4°E), and with the (URSI-based) IRI90 at different conditiones. Experimental total electron content (TEC) data are compared to both models. At the considered station, the profiles obtained by both models are reasonably in agreement amongst themselves and with the experimental data.

The TEC derived from the DGR3 model is in good agreement with experimental TEC, whereas, at high solar activity, IRI90 gives too high TEC values, especially during daytime.     

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.     

Ionospheric parameters in the European sector during the magnetic storm of August 25–26, 2018

Variations of ionospheric parameters Total Electron Content (TEC) by GNSS, critical frequency (foF2) by vertical sounding and electron density (Ne) by low-altitude satellite were studied at high, mid and low latitudes of the European sector during the magnetic storm of August 25–26, 2018. During the main phase of the storm the ionospheric F2-layer was under the positive disturbance at mid and low latitudes. Then the transition from the positive to negative ΔfoF2 values occurred at all latitudes. The recovery phase was characterized by negative ionospheric disturbance at all latitudes. This is due to the decrease of thermospheric O/N2 ratio during the recovery phase of the storm. The intense Es layers screened the reflections from the F2-layer on August 26th at high and at low latitudes but at different times. Some blackouts occurred due to the high absorption level at high latitudes. In general, foF2 and TEC data were highly correlated. The major Ne changes were at the low latitudes. In general, Ne data confirmed the ionospheric dynamics revealed with foF2 and TEC.     

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.     

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.     

Response of the ionosphere to super geomagnetic storms: Observations and modeling

The relative importance of the main drivers of positive ionospheric storms at low-mid latitudes is studied using observations and modeling for the first time. In response to a rare super double geomagnetic storm during 07–11 November 2004, the low-mid latitude (17°–48°N geomag. lat.) ionosphere produced positive ionospheric storms in peak electron density (NmF2) in Japan longitudes (≈125°–145°E) on the day of main phase (MP1) onset (06:30 LT) and negative ionospheric storms in American longitudes (≈65°–120°W) on the following day of MP1 onset (13:00–16:00 LT). The relative effects of the main drivers of the positive ionospheric storms (penetrating daytime eastward electric field, and direct and indirect effects of equatorward neutral wind) are studied using the Sheffield University Plasmasphere Ionosphere Model (SUPIM). The model results show that the penetrating daytime (morning–noon) eastward electric field shifts the equatorial ionisation anomaly crests in NmF2 and TEC (total electron content) to higher than normal latitudes and reduces their values at latitudes at and within the anomaly crests while the direct effects of the equatorward wind (that reduce poleward plasma flow and raise the ionosphere to high altitudes of reduced chemical loss) combined with daytime production of ionisation increase NmF2 and TEC at latitudes poleward of the equatorial region; the later effects can be major causes of positive ionospheric storms at mid latitudes. The downwelling (indirect) effect of the wind increases NmF2 and TEC at low latitudes while its upwelling (indirect) effect reduces NmF2 and TEC at mid latitudes. The net effect of all main drivers is positive ionospheric storms at low-mid latitudes in Japan longitude, which qualitatively agrees with the observations.     

国际参考电离层IRI－90与中国实测电离层月平均电子密度剖面形状的比较

将太阳活动峰年期间中国4个电离层站垂直探测得到的月平均电子浓度剖面与国际参考电离层IKI-90进行了系统的比较。结果表明国际参考电离层所计算的峰下电子含量(或峰下半厚)总体来说偏大。一天中白天符合较好,晚上较差,对中纬台站较好,对低纬台站较差。      

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.     

Comparison of observed TEC values with IRI-2007 TEC and IRI-2007 TEC with optional foF2 measurements predictions at an equatorial region,Chumphon, Thailand

In this research, as part of working towards improving the IRI over equatorial region, the total electron content (TEC) derived from GPS measurements and IRI-2007 TEC predictions at Chumphon station (10.72°N, 99.37°E), Thailand, during 2004–2006 is analyzed. The seasonal variation of the IRI-2007 TEC predictions is compared with the TEC from the IRI-2007 TEC model with the option of the actual F2 plasma frequency (foF2) measurements as well as the TEC from the GPS and International GNSS service (IGS). The Chumphon station is located at the equatorial region and the low latitude of 3.22°N. For a declining phase of the solar cycle (2004–2006), the study shows that the IRI-2007 TEC underestimates the IRI-2007 TEC with the foF2 observation at the nighttime by about 5 TECU. The maximum differences are about 15 TECU during daytime and 5 TECU during nighttime. The overestimation is more evident at daytime than at nighttime. When compared in terms of the root-mean square error (RMSE), we find that the highest RMSE between GPS TEC and IRI 2007 TEC is 14.840 TECU at 1230 LT in 2004 and the lowest average between them is 1.318 TECU at 0630 LT in 2006. The noon bite-out phenomena are clearly seen in the IRI-2007 TEC with and without optional foF2 measurements, but not on the GPS TEC and IGS TEC. The IRI TEC with optional foF2 measurements gives the lowest RMSE values between IRI TEC predicted and TEC measurement. However, the TEC measurements (GPS TEC and IGS TEC) are more correct to use at Chumphon station.     

TEC variability near northern EIA crest and comparison with IRI model

Monthly median values of hourly total electron content (TEC) is obtained with GPS at a station near northern anomaly crest, Rajkot (geog. 22.29°N, 70.74°E; geomag. 14.21°N, 144.9°E) to study the variability of low latitude ionospheric behavior during low solar activity period (April 2005 to March 2006). The TEC exhibit characteristic features like day-to-day variability, semiannual anomaly and noon bite out. The observed TEC is compared with latest International Reference Ionosphere (IRI) – 2007 model using options of topside electron density, NeQuick, IRI01-corr and IRI-2001 by using both URSI and CCIR coefficients. A good agreement of observed and predicted TEC is found during the daytime with underestimation at other times. The predicted TEC by NeQuick and IRI01-corr is closer to the observed TEC during the daytime whereas during nighttime and morning hours, IRI-2001 shows lesser discrepancy in all seasons by both URSI and CCIR coefficients.     

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

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

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.     

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.     

Modelling bottom and topside electron density and TEC with profile data from topside ionograms

The paper describes the technique that has been implemented to model the electron density distribution above and below the F2 peak making use of only the profiles obtained from the INTERCOSMOS-19 topside ionograms. Each single profile from the satellite height to the ionosphere peak has been fitted by a semi-Epstein layer function of the type used in the DGR model with shape factor variable with altitude. The topside above the satellite height has been extrapolated to match given values of plasmaspheric electron densities to obtain the full topside profile. The bottomside electron density has been calculated by using the maximum electron density and its altitude estimated from the topside ionogram as input for a modified version of the DGR derived profiler that uses model values for the foF1 and foE layers of the ionosphere. Total electron content has also been calculated. Longitudinal cross sections of vertical profiles from latitudes 50° N to 50° S latitude are shown for low and high geomagnetic activity. These cross sections indicate the equatorial anomaly effect and the changes of the shape of low latitude topside ionosphere during geomagnetic active periods. These results and the potentiality of the technique are discussed.     

Global empirical models of the density peak height and of the equivalent scale height for quiet conditions

Monthly average electron density profiles have been calculated from hourly electron density N(h) recorded in 26 digisonde stations distributed worldwide encompassing the time interval 1998–2006. The ionospheric electron density peak height of the F2 region, hmF2, and the effective scale height at the hmF2, Hm, deduced from average profiles have been analyzed to obtain the quiet-time behavior and have been analytically modeled by the spherical harmonic analysis (SH) technique using the modip latitude as the coordinate of the reference system. The coefficients of the SH models of hmF2 and Hm are bounded to the solar activity, and the temporal and seasonal variations are considered by Fourier expansion of the coefficients. The SH models provide a tool to predict hmF2 and Hm located anywhere in the range of latitudes between of 70°N and 70°S and at any time. The SH analytical model for hmF2 improves the fit to the observations by 10% in average compared to the IRI prediction, and it might improve the IRI prediction of hmF2 by more than 30% at high and low latitudes. The analytical model for Hm predicts the quiet behavior of the effective scale height with accuracy better than 15% in average which enables to obtain a good estimation of vertical profiles. These results could be useful to estimate information for the topside profile formulation.     

Collocated ionosonde and dense GPS/GLONASS network measurements of midlatitude MSTIDs

To analyze midlatitude medium-scale travelling ionospheric disturbances (MSTIDs) over Kazan (55.5°N, 49°E), Russia, the sufficiently dense network of GNSS receivers (more than 150 ground-based stations) were used. For the first time, daytime MSTIDs in the form of their main signature (band structure) on high-resolution two-dimensional maps of the total electron content perturbation (TEC maps) are compared with ionosonde data with a high temporal resolution. For a pair of events, a relationship between southwestward TEC perturbations and evolution of F2 layer traces was established. So F2 peak frequency varied in antiphase to TEC perturbations. The ionograms show that during the movement of plasma depletion band (overhead ionosonde) the F2 peak frequency is the highest, and vice versa, for the plasma enhancement band, the F2 peak frequency is the lowest. One possible explanation may be a greater inclination of the radio beam from the vertical during the placement of a plasma enhancement band above the ionosonde, as evidenced by the absence of multiple reflections and the increased occurrence rate of additional cusp trace. Another possible explanation may be the redistribution of the electron content in the topside ionosphere with a small decrease in the F peak concentration of the layer with a small increase in TEC along the line-of-sight. Analysis of F2 peak frequency variation shows that observed peak-to-peak values of TEC perturbation equal to 0.4 and 1 TECU correspond to the values of ΔN/N equal to 13% and 28%. The need for further research is evident.     

The variability and predictability of the IRI B0, B1 parameters over Grahamstown,South Africa

The International Reference Ionosphere (IRI) parameters B₀ and B₁ provide a representation of the thickness and shape, respectively, of the F2 layer of the bottomside ionosphere. These parameters can be derived from electron density profiles that are determined from vertical incidence ionograms. This paper aims to illustrate the variability of these parameters for a single mid latitude station and demonstrate the ability of the Neural Network (NN) modeling technique for developing a predictive model for these parameters. Grahamstown, South Africa (33.3°S, 26.5°E) was chosen as the mid latitude station used in this study and the B₀ and B₁ parameters for an 11 year period were determined from electron density profiles recorded at that station with a University of Massachusetts Lowell Center for Atmospheric Research (UMLCAR) Digisonde. A preliminary single station NN model was then developed using the Grahamstown data from 1996 to 2005 as a training database, and input parameters known to affect the behaviour of the F2 layer, such as day number, hour, solar and magnetic indices. An analysis of the diurnal, seasonal and solar variations of these parameters was undertaken for the years 2000, 2005 and 2006 using hourly monthly median values. Comparisons between the values derived from measured data and those predicted using the two available IRI-2001 methods (IRI tables and Gulyaeva, T. Progress in ionospheric informatics based on electron density profile analysis of ionograms. Adv. Space Res. 7(6), 39–48, 1987.) and the newly developed NN model are also shown in this paper. The preliminary NN model showed that it is feasible to use the NN technique to develop a prediction tool for the IRI thickness and shape parameters and first results from this model reveal that for the mid latitude location used in this study the NN model provides a more accurate prediction than the current IRI model options.     

Equatorial total electron content (TEC) at low and high solar activity

Total electron content (TEC) derived from ionosonde data recorded at the station of Korhogo (Lat = 9.33°N, Long = 5.43°W, Dip = 0.67°S) are compared to the International Reference Ionosphere (IRI) model predicted TEC for high (1999) and low (1994) solar activity conditions. The results show that the model represents the diurnal variation of the TEC as well as a solar activity and seasonal dependence. This variation is closer to that of the ionosonde-inferred TEC at high solar activity. However, at low solar activity the IRI overestimates the ionosonde-inferred TEC. The relative deviation ΔTEC is more prominent in the equinoctial seasons during nighttime hours where it is as high as 70%. At daytime hours, the relative deviation is estimated to 0–30%.     

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.