Ionospheric variation at Thailand equatorial latitude station: Comparison between observations and IRI-2001 model predictions

In this paper, the F2-layer critical frequency (foF2) and peak height (hmF2) measured by the FM/CW ionosonde at Thailand equatorial latitude station, namely Chumphon (10.72°N, 99.37°E, dip 3.22) are presented. The measurement data during low solar activity from January 2004 to December 2006 are analyzed based on the diurnal, seasonal variation. The results are then compared with IRI-2001 model predictions. Our study shows that: (1) In general, both the URSI and CCIR options of the IRI model give foF2 close to the measured ones, but the CCIR option produces a smaller range of deviation than the URSI option. The agreement during daytime is generally better than during nighttime. Overestimation mostly occurs in 2004 and 2006, while underestimation is during pre-sunrise hours in June solstice in 2005. The peak foF2 around sunset is higher during March equinox and September equinox than the other seasons, with longer duration of maximum levels in March equinox than September equinox. Large coefficients of variability foF2 occur during pre-sunrise hours. Meanwhile, the best agreement between the observed foF2 and the IRI model is obtained in June solstice. (2) In general, The IRI (CCIR) model predicts the observed hmF2 well during daytime in June solstice from 2004–2006, but it overestimates during March equinox, September equinox and December solstice. For nighttime, the model overestimates hmF2 values for all seasons especially during March equinox and September equinox. However, the model underestimates hmF2 values during September equinox and for some cases during June solstice and December solstice at pre-sunrise. The agreement between the IRI model and the hmF2(M3000_OBS) is worst around noontime, post-sunset and pre-sunrise hours. All comparative studies give feedback for new improvements of CCIR and URSI IRI models.     

Validation of the STORM model in IRI-2001 at a high latitude station

The International Reference Ionosphere IRI-2001 model contains geomagnetic activity dependence based on an empirical storm time ionospheric correction (STORM model). An extensive validation of the STORM model for the middle latitude region has been performed. In this paper the ability of the STORM model to predict foF2 values at high latitudes is analyzed. For this, ionosonde data obtained at Base Gral. San Martin (68.1°S, 293°E) are compared with those obtained by the IRI-2001 model with or without storm correction during four geomagnetic storms that occurred in 2000 (Rz₁₂ = 117) and 2001 (Rz₁₂ = 111). The results show that predicted values with the STORM model follow the behaviour of foF2 experimental data better than without the STORM model. The relative deviation between measured and predicted foF2 reaches values of up to 24% and 43% with and without the STORM model in IRI-2001, during the main phase of the storms. In order to explain increases of electron density that occurred prior to the storm onset and also decreases of electron density observed during the first part of the recovery of the storm, possible physical mechanisms are discussed.     

An evaluation of the IRI-2007 storm time model at low latitude stations

This paper discusses the ability of the International Reference Ionosphere IRI-2007 storm time model to predict foF2 ionospheric parameter during geomagnetic storm periods. Experimental data (based on availability) from two low latitude stations: Vanimo (geographic coordinates, 2.7 °S, 141.3 °E, magnetic coordinates, 12.3 °S, 212.50 °E) and Darwin (geographic coordinates, 12.45 °S, 130.95 °E, magnetic coordinates, 22.9 °S, 202.7 °E) during nine storms that occurred in 2000 (Rz₁₂ = 119), 2001(Rz₁₂ = 111) and 2003 (Rz₁₂ = 64) are compared with those obtained by the IRI-2007 storm model. The results obtained show that the percentage deviation between the experimental and IRI predicted foF2 values during these storm periods is as high as 100% during the main and recovery phases. Based on the values of “relative deviation module mean” (RDMM) obtained (i.e. between 0.08 and 0.60), it is observed that there is a reasonable to poor agreement between measured foF2 values and the IRI-storm model prediction values during main and recovery phases of the storms under investigation. As a result, in addition to other studies that have been carried out from different sectors, more studies are required to be carried out. This will enable IRI community to improve on the present performance of the model. In general the IRI-storm model predictions follow normal trend of the foF2 measured values but does not reproduce well the measured 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.     

The response of African equatorial foF2 to geomagnetic storms: Comparison between observations and IRI-2007 predictions

This study examines the response of the African equatorial ionospheric foF2 to different levels of geomagnetic storms, using the foF2 hourly data for the year 1989 from Ouagadougou (12.4°N, 1.5°W, dip: 2.8°N). The study also compares the observed data for the selected storm periods with the latest IRI model (IRI-2007). The foF2 values (both observed and predicted) show typical features of daytime peak and post-midnight minimum peak. The response of the ionospheric foF2 over Ouagadougou to storms events, during the night-time and post-midnight hours indicates negative responses of the ionospheric foF2, while that of the daytime hours indicates positive responses. For the investigation on the variability of the observed foF2 with respect to IRI-2007 model, with the exception of the analysis of the 20–22, October, 1989 data, where a midday peak was also observed on the first day, this study reveals two characteristic daily foF2 variability peaks: post-midnight and evening peaks. In addition, for all the geomagnetic storms considered, the URSI and CCIR coefficients of the IRI-2007 model show reasonable correspondence with each other, except for some few discrepancies. For instance, the event of 28–30 August, 1989 shows comparatively higher variability for the URSI coefficient, and at the foF2 peak values, the event of 20–22 October, 1989 shows that the CCIR coefficient is more susceptible to foF2 variability than the URSI coefficient. This study is aimed at providing African inputs for the future improvement of the IRI model.     

Bottom side profiles for two close stations at the southern crest of the EIA: Differences and comparison with IRI-2012 and NeQuick2 for low and high solar activity

Bottom side electron density profiles for two stations at the southern crest of the Equatorial Ionization Anomaly (EIA), São José dos Campos (23.1°S, 314.5°E, dip latitude 19.8°S; Brazil) and Tucumán (26.9°S, 294.6°E, dip latitude 14.0°S; Argentina), located at similar latitude and separated by only 20° in longitude, have been compared during equinoctial, winter and summer months under low (year 2008, minimum of the solar cycle 23/24) and high solar activity (years 2013–2014, maximum of the solar cycle 24) conditions. An analysis of parameters describing the bottom side part of the electron density profile, namely the peak electron density NmF2, the height hmF2 at which it is reached, the thickness parameter B₀ and the shape parameter B₁, is carried out. Further, a comparison of bottom side profiles and F-layer parameters with the corresponding outputs of IRI-2012 and NeQuick2 models is also reported. The variations of NmF2 at both stations reveal the absence of semi-annual anomaly for low solar activity (LSA), evidencing the anomalous activity of the last solar minimum, while those related to hmF2 show an uplift of the ionosphere for high solar activity (HSA). As expected, the EIA is particularly visible at both stations during equinox for HSA, when its strength is at maximum in the South American sector. Despite the similar latitude of the two stations upon the southern crest of the EIA, the anomaly effect is more pronounced at Tucumán than at São José dos Campos. The differences encountered between these very close stations suggest that in this sector relevant longitudinal-dependent variations could occur, with the longitudinal gradient of the Equatorial Electrojet that plays a key role to explain such differences together with the 5.8° separation in dip latitude between the two ionosondes. Furthermore at Tucumán, the daily peak value of NmF2 around 21:00 LT during equinox for HSA is in temporal coincidence with an impulsive enhancement of hmF2, showing a kind of “elastic rebound” under the action of the EIA. IRI-2012 and NeQuick2 bottom side profiles show significant deviations from ionosonde observations. In particular, both models provide a clear underestimation of the EIA strength at both stations, with more pronounced differences for Tucumán. Large discrepancies are obtained for the parameter hmF2 for HSA during daytime at São José dos Campos, where clear underestimations made by both models are observed. The shape parameter B₀ is quite well described by the IRI-2012 model, with very good agreement in particular during equinox for both stations for both LSA and HSA. On the contrary, the two models show poor agreements with ionosonde data concerning the shape parameter B₁.     

Comparison of ionosphere characteristic parameters obtained by ionosonde with IRI-2007 model over Southeast Asia

In this work, the foF2 and hmF2 parameters at the conjugate points near the magnetic equator of Southeast Asia are studied and compared with the International Reference Ionosphere (IRI) model. Three ionosondes are installed nearly along the magnetic meridian of 100°E; one at the magnetic equator, namely Chumphon (10.72°N, 99.37°E, dip angle 3.0°N), and the other two at the magnetic conjugate points, namely Chiang Mai (18.76°N, 98.93°E, dip angle 12.7°N) and Kototabang (0.2°S, 100.30°E, dip angle 10.1°S). The monthly hourly medians of the foF2 and hmF2 parameters are calculated and compared with the predictions obtained from the IRI-2007 model from January 2004 to February 2007. Our results show that: the variations of foF2 and hmF2 predicted by the IRI-2007 model generally show the similar feature to the observed data. Both parameters generally show better agreement with the IRI predictions during daytime than during nighttime. For foF2, most of the results show that the IRI model overestimates the observed foF2 at the magnetic equator (Chumphon), underestimates at the northern crest (Chiang Mai) and is close to the measured ones at the southern crest of the EIA (Kototabang). For hmF2, the predicted hmF2 values are close to the hmF2(M3000F2_OBS) during daytime. During nighttime, the IRI model gives the underestimation at the magnetic equator and the overestimation at both EIA crests. The results are important for the future improvements of the IRI model for foF2 and hmF2 over Southeast Asia region.     

Diurnal and seasonal variation of F2-layer ionospheric parameters at equatorial ionization anomaly crest region and their comparison with IRI-2001

Diurnal and seasonal variations of critical frequency of ionospheric F2-region ‘foF2’ and the height of peak density ‘hmF2’ are studied using modern digital ionosonde observations of equatorial ionization anomaly (EIA) crest region, Bhopal (23.2°N, 77.6°E, dip 18.5°N), during solar minimum period 2007. Median values of these parameters are obtained at each hour using manually scaled data during different seasons and compared with the International Reference Ionosphere-2001 model predictions. The observations suggest that on seasonal basis, the highest values of foF2 are observed during equinox months, whereas highest values of hmF2 are obtained in summer and lowest values of both foF2 and hmF2 are observed during winter. The observed median and IRI predicted values of foF2 and hmF2 are analyzed with upper and lower bound of inter-quartile range (IQR) and it is find out that the observed median values are well inside the inter-quartile range during the period of 2007. Comparison of the recorded foF2 and hmF2 values with the IRI-2001 output reveals that IRI predicted values exhibit better agreement with hmF2 as compared to foF2. In general, the IRI model predictions show some agreement with the observations during the year 2007. Therefore it is still necessary to implement improvements in order to obtain better predictions for EIA regions.     

Comparison between observed ionospheric foF2 and IRI-2001 predictions over periods of severe geomagnetic activities at Grahamstown,South Africa

The observed ionospheric F2 critical frequency (foF2) values over a South Africa mid-latitude station, Grahamstown, (geographic coordinates: 33.3°S, 26.5°E), were analysed and compared with International Reference Ionosphere (IRI) model, using the CCIR (Comite´ Consultatif International des Radio communications) and URSI (Union Radio-Scientifique Internationale) coefficients, during four geomagnetically disturbed days in the year 2000. These days are April 5, May 23, August 10 and September 15. The data were analysed for five days around the storm day. Comparisons between the IRI-2001 predicted foF2 values, using both CCIR and URSI coefficients and the observed values are shown with their root-mean-square error (RMSE) and the relative deviation module mean (rdmm) for the various storm periods. The CCIR option performed more accurately than the URSI option.     

The subauroral electron density trough: Comparison between satellite observations and IRI-2007 model estimates

We compare electron density predictions of the International Reference Ionosphere (IRI-2007) model with in situ measurements of the satellites CHAMP and GRACE for the years from 2005 to 2010 over the subauroral regions. The electron density between 58° and 68° Mlat are considered. The trough region Ne peaks during local summers and attain the valley during local winter. Around −100°E and 60°E, two larger electron density sectors features can be seen in both hemispheres during all three seasons, which attributed to the electron extending from middle latitude to trough region. From 2005 to the beginning of 2010, the model overestimates the trough region Ne by 20% on average and the decrease of Ne in this region can also be seen during the last solar minimum. In the southern hemisphere, the model prediction shows quite well consistence with the observation during all three seasons while the huge difference between observations and model estimation implies that the IRI-2007 model needs significant improvement to predict better the trough region in northern hemisphere.     

Comparison of peak characteristics of F2 ionospheric layer over Tehran region at a low solar activity period with IRI-2001 and IRI-2007 models predictions

In the present work values of peak electron density (NmF2) and height of F2 ionospheric layer (hmF2) over Tehran region at a low solar activity period are compared with the predictions of the International Reference Ionosphere models (IRI-2001 and IRI-2007). Data measured by a digital ionosonde at the ionospheric station of the Institute of Geophysics, University of Tehran from July 2006 to June 2007 are used to perform the calculations. Formulations proposed by  and  are utilized to calculate the hmF2. The International Union of Radio Science (URSI) and International Radio Consultative Committee (CCIR) options are employed to run the IRI-2001 and IRI-2007 models. Results show that both IRI-2007 and IRI-2001 can successfully predict the NmF2 and hmF2 over Tehran region. In addition, the study shows that predictions of IRI-2007 model with CCIR coefficient has closer values to the observations. Furthermore, it is found that the monthly average of the percentage deviation between the IRI models predictions and the values of hmF2 and NmF2 parameters are less than 10% and 21%, respectively.     

An empirical model of electron density in low latitude at 600 km obtained by Hinotori satellite

An empirical model of electron density (Ne) was constructed by using the data obtained with an impedance probe on board Japanese Hinotori satellite. The satellite was in circular orbit of the height of 600 km with the inclination of 31 degrees from February 1981 to June 1982. The constructed model gives Ne at any local time with the time resolution of 90 min and between −25 and 25 degrees in magnetic latitude with its resolution of 5 degrees in the range of F_10.7 from 150 to 250 under the condition of K_p < 4. Spline interpolations are applied to the functions of day of year, geomagnetic latitude and solar local time, and linear interpolation is applied to the function of F_10.7. Longitude dependence of Ne is not taken into account. Our density model can reproduce solar local time variation of electron density at 600 km altitude better than current International Reference Ionosphere (IRI2001) model which overestimates Ne in night time and underestimates Ne in day time. Our density model together with electron temperature model which has been constructed before will enable more understanding of upper ionospheric phenomenon in the equatorial region.     

Vertical TEC variations and model during low solar activity at a low latitude station,Xiamen

GPS-derived vertical TEC recorded at Xiamen (24.5°N, 118.1°E, geomagnetic latitude 13.2°N), China, during year 2006 is analyzed for the first time and compared to that predicted by ionosphere model SPIM recommend by ISO. A manifest seasonal anomaly is found with the high value during equinoctial season and low value during summer and winter season. Relative standard deviation for VTEC shows high value at around midnight and before sunrise. The correlation analysis exhibits that the variation of VTEC has a very weak relation with geomagnetic and solar activities (Dst, AP, SSN and F10.7). Comparative results reveal that the SPIM overestimates the observed VTEC at most of the time.     

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.     

Ionospheric behavior of the F2 peak parameters foF2 and hmF2 at Hainan and comparisons with IRI model predictions

Monthly median values of foF2, hmF2 and M(3000)F2 parameters, with hourly time interval resolution for the diurnal variation, obtained with DPS-4 digisonde observations at Hainan (19.4°N, 109.0°E) are used to study the low latitude ionospheric variation behavior. The observational results are compared with the International Reference Ionospheric Model (IRI) predictions. The time period coverage of the data used for the present study is from March 2002 to February 2005. Our present study showed that: (1) In general, IRI predictions using CCIR and URSI coefficients follow well the diurnal and seasonal variation patterns of the experimental values of foF2. However, CCIR foF2 and URSI foF2 IRI predictions systematically underestimate the observed results during most time period of the day, with the percentage difference ΔfoF2 (%) values changing between about −5% and −25%, whereas for a few hours around pre-sunrise, the IRI predictions generally overestimate the observational ones with ΔfoF2 (%) sometimes reaching as large as ∼30%. The agreement between the IRI results and the observational ones is better for the year 2002 than for the other years. The best agreement between the IRI results and the observational ones is obtained in summer when using URSI coefficients, with the seasonal average values of ΔfoF2 (%) being within the limits of ±10%. (2) In general, the IRI predicted hmF2 values using CCIR M(3000)F2 option shows a poor agreement with the observational results. However, when using the measured M(3000)F2 as input, the diurnal variation pattern of hmF2 given by IRI2001 has a much better agreement with the observational one with the detailed fine structures including the pre-sunrise and post-sunset peaks reproduced reasonably well. The agreement between the IRI predicted hmF2 values using CCIR M(30,000)F2 option and the observational ones is worst for the afternoon to post-midnight hours for the high solar activity year 2002. During daytime hours the agreement between the hmF2 values obtained with CCIR M(30,000)F2 option and the observational ones is best for summer season. The discrepancy between the observational hmF2 and that obtained with CCIR M(30,000)F2 option stem from the CCIR M(3000)F2 model, which does not produce the small scale structures observed in the measured M(3000)F2.     

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.     

Variability of the behavior of the bottomside (B0, B1) parameters obtained from the ground-based ionograms at China’s low latitude station

In this paper, data of (B0, B1) parameters deduced from the electron density profiles that are inverted from the ionograms recorded at Hainan (19.4°N, 109.0°E), China during a three year period from March 2002 to February 2005 are used to study the diurnal and seasonal variation of (B0, B1) parameters at low latitude. The observational results are compared with the IRI2001 model predictions. Variability study of (B0, B1) in terms of percentage ratio of the inter-quartiles to the median values and correlative analysis between (B0, B1) parameters and other ionospheric characteristics such as hmF2 and M(3000)F2 are also made. Our present study showed that: (1) for daytime hours, the IRI2001 model results with new table option (B0_Tab) is in a better agreement with the observational results (B0_Obs) than the IRI2001 model results with Gulyaeva option (B0_Gul) for summer season, whereas B0_Gul is in a better agreement with B0_Obs than B0_Tab for winter season. For nighttime, in general, B0_Gul is in a better agreement with B0_Obs than B0_Tab. For other occasions, both B0_Tab and B0_Gul showed some systematic deviations from the observational ones. Moreover, the deviations of B0_Tab and B0_Gul from B0_Obs showed opposite trends; (2) the monthly upper (lower) quartiles of (B0, B1) parameter showed a good linear relationship with the monthly median values, this makes it possible to do the regression analysis between the monthly upper (lower) quartiles and the monthly median values, which can give a measure of the variability of these parameters. In terms of the percentage ratio of the inter-quartiles to the median values, the variability of B0 showed a diurnal variation ranging between 22% and 36% with maximum value occurring at pre-sunrise hours, whereas the variability of B1 showed a diurnal variation ranging between 15% and 30% with higher value by daytime than at night; (3) B0 shows high linear correlative relationships with hmF2 and M(3000)F2 for most of the local time period of a day except for a few hours around midnight, whereas B1 showed high linear correlations with B0, hmF2 for daytime hours, but not for nighttime hours. This suggests that it maybe is possible to obtain the synthetic database of (B0, B1) parameter or to construct the model of (B0, B1) using database of hmF2 or M(3000)F2 which is much easier to obtain from experimental measurements.     

Validation of B0 and B1 in the IRI 2001 model at low solar activity for Ilorin an equatorial station

A new set of data obtained at low solar activity from Ilorin, Nigeria (geog. latitude 8.5°N, geog longitude, 4.6°E, dip 4.1°S) is used to validate the IRI 2001 model at low solar activity. The results show in general a good agreement between model and observed B0 at night but an over estimation during daytime. The overestimation is greatest during the morning period (0600LT–1000LT). The model prediction for B1 is fairly good at night and during the day. A dependence of B0 on solar zenith angle χ is observed during the daytime. A formulation of the form B0 = A[cos(χ)ⁿ] is therefore proposed. Values of the constants n and A were determined for the period of low solar activity for this station.     

Evaluation of the IRI-2016 and NeQuick electron content specification by COSMIC GPS radio occultation,ground-based GPS and Jason-2 joint altimeter/GPS observations

We examined performance of two empirical profile-based ionospheric models, namely IRI-2016 and NeQuick-2, in electron content (EC) and total electron content (TEC) representation for different seasons and levels of solar activity. We derived and analyzed EC estimates in several representative altitudinal intervals for the ionosphere and the plasmasphere from the COSMIC GPS radio occultation, ground-based GPS and Jason-2 joint altimeter/GPS observations. It allows us to estimate a quantitative impact of the ionospheric electron density profiles formulation in several altitudinal intervals and to examine the source of the model-data discrepancies of the EC specification from the bottom-side ionosphere towards the GPS orbit altitudes. The most pronounced model-data differences were found at the low latitude region as related to the equatorial ionization anomaly appearance. Both the IRI-2016 and NeQuick-2 models tend to overestimate the daytime ionospheric EC and TEC at low latitudes during all seasons of low solar activity. On the contrary, during high solar activity the model results underestimated the EC/TEC observations at low latitudes. We found that both models underestimated the EC for the topside ionosphere and plasmasphere regions for all levels of solar activity. For low solar activity, the underestimated EC from the topside ionosphere and plasmasphere can compensate the overestimation of the ionospheric EC and, consequently, can slightly decrease the resulted model overestimation of the ground-based TEC. For high solar activity, the underestimated EC from the topside ionosphere and plasmasphere leads to a strengthening of the model underestimation of the ground-based TEC values. We demonstrated that the major source of the model-data discrepancies in the EC/TEC domain comes from the topside ionosphere/plasmasphere system.