Equatorial M(3000)F2 estimation of F2-layer models peak heights during low solar activity

M(3000)F2 estimation of h_mF2 based on four different formulated models viz: (1) Shimazaki (1955) (2) Bradley and Dudeney (1973), (3) Dudeney (1974) and (4) Bilitza et al. (1979) at an equatorial station in West Africa during low solar activity period (1995) are used to validate its conformity with observed and International Reference Ionosphere (IRI) model. Local time analyses of data from fifteen (15) selected days during the January and July solstices and April and October equinoxes are used. The results obtained show that the M(3000)F2 estimation of h_mF2 from the ionosonde-measured values using the Ionospheric Prediction Service (IPS-42) sounder compared to the observed values which were deduced using an algorithm from scaled virtual heights of quiet day ionograms are highly correlated with Bilitza model. International Reference Ionosphere (IRI 2007) model for the equatorial region also agrees with the formulation developed by Bilitza et al. (1979) for the four different seasons of the year. h_mF2 is highest (425 km) in summer (June solstice) season and lowest (386 km) in autumn (September equinox) season with daytimes peaks occurring at 1100–1200 LT during the solstices and at 1000 LT during the equinoxes respectively. Also, the post-sunset peaks are highest (362 km) at the spring (March equinox) and lowest (308 km) at the summer (June solstice) both occurring between 1800 and 2000 LT.     

Nighttime thermospheric meridional neutral winds inferred from ionospheric h′F and hpF2 data

Nighttime thermospheric meridional winds aligned to the magnetic meridian have been inferred using h′F and hpF2 ionosonde data taken from two equatorial stations, Manaus (2.9°S, 60.0°W, dip latitude 6.0°N) and Palmas (10.17°S, 48.2°W, dip latitude 6.2°S), and one low-latitude station, Sao Jose dos Campos (23.21°S, 45.86°W, dip latitude 17.26°S), during geomagnetic quiet days of August and September, 2002. Using an extension of the ionospheric servo model and a simple formulation of the diffusive vertical drift velocity, the magnetic meridional component of the thermospheric neutral winds is inferred, respectively, at the peak (hpF2) and at the base (h′F) heights of the F region over Sao Jose dos Campos. An approach has been included in the models to derive the effects of the electrodynamic drift over Sao Jose dos Campos from the time derivative of hpF2 and h′F observed at the equatorial stations. The magnetic meridional winds inferred from the two methods, for the months of August and September, are compared with winds calculated using the HWM-90 model and with measurements from Fabry–Perot technique. The results show varying agreements and disagreements. Meridional winds calculated from hpF2 ionospheric data (servo model) may produce errors of about 59 m/s, whereas the method calculated from the F-region base height (h′F) ionospheric data gives errors of about 69 m/s during the occurrence of equatorial spread-F.     

Short-term forecasting regional model to predict M(3000)F2 over the European sector: Comparisons with the IRI model during moderate,disturbed, and very disturbed geomagnetic conditions

The hourly measurements of M(3000)F2 (M(3000)F2_meas) and the hourly quiet-time values of M(3000)F2 (M(3000)F2_QT) relative to the ionospheric observatories of Poitiers, Lannion, Dourbes, Slough, Rome, Juliusruh, Kaliningrad, Uppsala, Lyckesele, Sodankyla, and Kiruna as well as the hourly time-weighted accumulation series derived from the geomagnetic planetary index ap (ap(τ)), were considered during the period January 1957–December 2003 and used for the development of 11 short-term forecasting local models (STFLM) of M(3000)F2.     

The accuracy of data from ionosondes for the estimation of hmF2 and the identification of global change in the ionosphere

A long temporal series of simulated ionograms was generated with a superimposed secular variation of −14 km/century on the h_mF2 parameter. These ionograms were interpreted by the automatic scaling program Autoscala. By applying four different empirical formulas, four artificial series of h_mF2 were generated and then processed with the same methods used by other authors for real data sets. Data analysis of the simulated ionograms revealed the artificially imposed long-term trend. These results lead to the conclusion, that regardless of the empirical formula used, the accuracy of h_mF2 from ionosonde measurements would be adequate to observe a long-term trend of −14 km/century.     

Comparing IRI and a regional model with ionosonde measurements in Pakistan

We have used the technique of expansion in Empirical Orthogonal Functions (EOFs) to develop regional models of the critical frequencies of E and F₂ layers (f_oE, f_oF₂), peak height (h_mF₂), and semi-thickness of F₂ layer (Y_mF₂) over Pakistan. In the present study levels of solar activity specified by Smoothed Sunspot Number (R) from 10 to 200 are taken into account. The magnetic dip angle for the model ranges from 30° to 60°. We have compared the regional model and the International Reference Ionosphere (IRI) with measurements of three ionosondes in Pakistan. The model parameters f_oE and f_oF₂ are found overall comparable to the observed hourly median values during daytime at Karachi (geographic latitude = 24.95°N, longitude = 67.13°E, magnetic inclination = 37°), Multan (30.18°N, 71.48°E, 45°) and Islamabad (33.75°N, 73.13°E, 51.5°) during the years 1988, 1996 and 2000. For h_mF₂ the computed values by regional and IRI model for the year 1995 are found close to each other. However, for Y_mF₂the results are better during daytime as compared to nighttime.     

A different method to update monthly median hmF2 values

The height, h_mF2, and the electron density, N_mF2, of the F2 peak are key model parameters to characterize the actual state of the ionosphere. These parameters, or alternatively the propagation factor, M3000F2, and the critical frequency, f_oF2, of the F2 peak, which are related to h_mF2 and N_mF2, are used to anchor the electron density vertical profile computed with different models such as the International Reference Ionosphere ( Bilitza, 2002), as well as for radio propagation forecast purposes. Long time series of these parameters only exist in an inhomogeneous distribution of points over the surface of Earth, where dedicated instruments (typically ionosondes) have been working for many years. A commonly used procedure for representing median values of the aforementioned parameters all over the globe is the one recommended by the ITU-R ( ITU-R, 1997). This procedure, known as the Jones and Gallet mapping technique, was based on ionosondes measurements gathered from 1954 to 1958 by a global network of around 150 ionospheric stations (  and ). Even though several decades have passed since the development of that innovative work, only few efforts have been dedicated to establish a new mapping technique for computing h_mF2 and N_mF2 median values at global scale or to improve the old method using the increased observational database. Therefore, in this work three different procedures to describe the daily and global behavior of the height of the F2 peak are presented. All of them represent a different and simplified method to estimate h_mF2 and are based on different mathematical expressions. The advantages and disadvantages of these three techniques are analyzed, leading to the conclusion that the recommended procedure to represent h_mF2 is best characterized by a Spherical Harmonics expansion of degree and order equal to 15, since the differences between the h_mF2 values obtained with the Jones and Gallet technique and those obtained using the abovementioned procedure are of only 1%.     

Proposal of new models of the bottom-side B0 and B1 parameters for IRI

The time series of hourly electron density profiles N(h) obtained from 27 ionosonde stations distributed world-wide have been used to obtain N(h) average profiles on a monthly basis and to extract the expected bottom-side parameters that define the IRI profile under quiet conditions. The time series embrace the time interval from 1998 to 2006, which practically contains the entire solar cycle 23. The Spherical Harmonic Analysis (SHA) has been used as an analytical technique for modeling globally the B0 and B1 parameters as general functions on a spherical surface. Due to the irregular longitudinal distribution of the stations over the globe, it has been assumed that the ionosphere remains approximately constant in form for a given day under quiet conditions for a particular coordinate system. Since the Earth rotates under a Sun-fixed system, the time differences have been considered to be equivalent to longitude differences. The time dependence has been represented by a two-degree Fourier expansion to model the annual and semiannual variations and the year-by-year analyses of the B0 and B1 have furnished nine sets of spherical harmonic coefficients for each parameter. The spatial–temporal yearly coefficients have been further expressed as linear functions of Rz12 to model the solar cycle dependence. The resultant analytical model provides a tool to predict B0 and B1 at any location distributed among the used range of latitudes (70°N–50°S) and at any time that improves the fit to the observed data with respect to IRI prediction.     

Study of maximum electron density NmF2 at Karachi and Islamabad during solar minimum (1996) and solar maximum (2000) and its comparison with IRI

The monthly hourly medians of maximum electron density, N_mF₂, at two Pakistani ionospheric stations, Karachi and Islamabad, have been determined for solar minimum (1996) and solar maximum (2000) and compared with IRI predictions using the URSI coefficients. At night and pre-noon period the N_mF₂ values at both stations are almost equal during the 2 years. However, at post-noon the values at Karachi are considerably larger than those at Islamabad due to the equatorial or geomagnetic anomaly. Karachi (geomag. coord. 16.44°N, 139.08°E) lies near the region of the equatorial anomaly (+20 and −20 geomagnetic latitude), so most of the N_mF₂ values at Karachi are larger than those at Islamabad (geomag. coord. 24.46°N, 145.67°E). The maximum monthly values of N_mF₂ show a semi-annual variation at Karachi and Islamabad both during 1996 and 2000 as predicted by IRI.     

Comparison of the CHAMP radio occultation data with the Canadian advanced digital ionosonde in the Polar Regions

In this paper we present the results of the comparison of the retrieved electron density profiles of the Ionospheric Radio Occultation (IRO) experiment on board CHAMP (CHAllenging Minisatellite Payload), with the ground ionosonde profiles for the Polar Regions. IRO retrieved electron density profiles from CHAMP are compared with Canadian Advanced Digital Ionosonde (CADI) measurements at two vertical sounding stations well within the Polar Cap, Eureka (geog. 80°13′ N; 86°11′ W) and Resolute Bay (geog. 74°41′ N; 94°54′ W). We compared the ionospheric parameters such as the peak electron density of the F-layer (NmF₂) and the peak height of the F-layer (h_mF₂) for a 3-year period, 2004–2006. CHAMP derived NmF₂ shows reasonable agreement with the ionosonde retrieved NmF₂ for both the stations (0.76 and 0.71 correlation coefficient, for Eureka and Resolute Bay, respectively) whereas the h_mF₂ agreement is not that acceptable (0.25 and 0.37 correlation coefficient, respectively). The h_mF₂ from vertical sounding showed less spread than the CHAMP h_mF₂.     

Towards a new reference model of hmF2 for IRI

A numerical model of the peak height of the F2 layer, hmF2_top, is derived from the topside sounding database of 90,000 electron density profiles for a representative set of conditions provided by ISIS1, ISIS2, IK19 and Cosmos-1809 satellites for the period of 1969–1987. The model of regular hmF2 variations is produced in terms of local time, season, geomagnetic latitude, geodetic longitude and solar radio flux. No geomagnetic activity trends were discernible in the topside sounding data. The nighttime peak of hmF2_top evident for mid-latitudes disappears near the geomagnetic equator where a maximum of hmF2_top occurs at sunset hours when it can exceed 500 km at solar maximum. The hmF2 given by the IRI exceeds hmF2_top at the low solar activities. The hmF2_top, obtained by extrapolation of the first derivative of the topside profile to zero shows saturation similar to foF2 the greater the solar activity. The proposed model differs from hmF2 given by IRI based on M(3000)F2 to hmF2 conversion by empirical relationships in terms of foF2, foE and R12 with these quantities mapped globally by the ITU-R (former CCIR) from ground-based ionosonde data. The differences can be attributed to the different techniques of the peak height derivation, different epochs and different global distribution of the source data as well as the different mathematical functions involved in the maps and the model presentation.     

Equatorial vertical E × B drift velocities inferred from ionosonde measurements over Ouagadougou and the IRI-2007 vertical ion drift model

F-region vertical plasma drift velocities were deduced from the hourly hmF2 values acquired from ionogram data over a near dip equatorial station Ouagadougou (12.4°N, 358.5°E, dip angle 5.9°N) in Africa. Our results are compared against the global empirical model of Scherliess and Fejer (1999) incorporated in the IRI model (IRI-2007) for 1600 to 0800 LT from 1 year of data during sunspot maximum year of 1989 (yearly average solar flux intensity, F10.7 = 192) corresponding to the peak phase of solar cycle 22, under magnetically quiet conditions. The drifts are entirely downward between 2000 and 0500 LT bin for both techniques and the root mean square error (RMSE) between the modeled and the ionosonde vertical plasma drifts during these periods is 3.80, 4.37, and 4.74 m/s for June solstice, December solstice and equinox, respectively. Ouagadougou average vertical drifts show evening prereversal enhancement (PRE) velocity peaks (V_ZP) of about 16, 14, and 17 m/s in June solstice, December solstice, and equinox, respectively, at 1900–2000 LT; whereas global empirical model average drifts indicate V_ZP of approximately 33 m/s (June solstice), 29 m/s (December solstice), and 50 m/s (equinox) at 1800 LT. We find very weak and positive correlation (+0.10376) between modeled V_ZP versus F10.7, while ionosonde V_ZP against F10.7 gives worst and opposite correlation (−0.05799). The results also show that modeled V_ZP–Ap indicates good and positive correlation (+0.64289), but ionosonde V_ZP–Ap exhibits poor and negative correlation (−0.22477).     

Estimation of the elastic Earth parameters (h2, l2) using SLR data

We present results for the global elastic parameters h₂ and l₂ derived from the analysis of Satellite Laser Ranging (SLR) data. SLR data for the two satellites LAGEOS 1 and LAGEOS 2 observed during 2.5 years from January 3, 2005 until July 1, 2007 with 18 globally distributed ground stations were analysed using different approaches. The analysis was done separately for the two satellites and approaches to estimate the two elastic parameters independently and together were performed. We do a sequential analysis and study the stability of the estimates as a function of length of the data set used. The adjusted final values for h₂ equal to 0.6151 ± 0.0008 and 0.6152 ± 0.0008, and those for l₂ equal to 0.0886 ± 0.0003 and 0.0881 ± 0.0003 for LAGEOS 1 and LAGEOS 2 tracking data are compared to other independently derived estimates. These parameters and their errors achieve stability at about the 24 and 27 month time interval for h₂ and l₂, respectively.     

Effect of vertical plasma transport on ionospheric F2-region parameters at equatorial latitude

The variability of the F2-layer even during magnetically quiet times are fairly complex owing to the effects of plasma transport. The vertical E × B drift velocities (estimated from simplified electron density continuity equation) were used to investigate the seasonal effects of the vertical ion drifts on the bottomside daytime ionospheric parameters over an equatorial latitude in West Africa, Ibadan, Nigeria (Geographic: 7.4°N, 3.9°E, dip angle: 6°S) using 1 year of ionsonde data during International Geophysical Year (IGY) of 1958, that correspond to a period of high solar activity for quiet conditions. The variation patterns between the changes of the vertical ion drifts and the ionospheric F2-layer parameters, especially; f_oF₂ and h_mF₂ are seen remarkable. On the other hand, we observed strong anti-correlation between vertical drift velocities and h′F in all the seasons. We found no clear trend between N_mF₂ and h_mF₂ variations. The yearly average value of upward daytime drift at 300 km altitude was a little less than the generally reported magnitude of 20 ms⁻¹ for equatorial F-region in published literature, and the largest upward velocity was roughly 32 ms⁻¹. Our results indicate that vertical plasma drifts; ionospheric F2-layer peak height, and the critical frequency of F2-layer appear to be somewhat interconnected.     

Variation in the ionospheric propagating factor M(3000)F2 at Ouagadougou,Burkina Faso

We use hourly monthly median values of propagation factor M(3000)F2 data observed at Ouagadougou Ionospheric Observatory (geographic12.4°N, 1.5°W; 5.9^o dip), Burkina Faso (West Africa) during the years Januar1987–December1988 (average F10.7 < 130 × 10⁻²² W/m²/Hz, representative of low solar flux conditions) and for January 1989–December1990 (average F10.7 ? 130 × 10⁻²² W/m²/Hz, representative of high solar epoch) for magnetically quiet conditions to describe local time, seasonal and solar cycle variations of equatorial ionospheric propagation factor M(3000)F2 in the African region. We show that that seasonal trend between solar maximum and solar minimum curves display simple patterns for all seasons and exhibits reasonable disparity with root mean square error (RMSE) of about 0.31, 0.29 and 0.26 for December solstice, June solstice and equinox, respectively. Variability Σ defined by the percentage ratio of the absolute standard deviation to the mean indicates significant dissimilarity for the two solar flux levels. Solar maximum day (10–14 LT) and night (22–02 LT) values show considerable variations than the solar minimum day and night values. We compare our observations with those of the IRI 2007 to validate the prediction capacity of the empirical model. We find that the IRI model tends to underestimate and overestimate the observed values of M(3000)F2, in particular, during June solstice season. There are large discrepancies, mainly during high solar flux equinox and December solstice between dawn and local midnight. On the other hand, IRI provides a slightly better predictions for M(3000)F2 between 0900 and 1500 LT during equinox low and high solar activity and equinox high sunspot number. Our data are of great importance in the area of short-wave telecommunication and ionospheric modeling.     

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.     

Equatorial F-region vertical ion drifts during quiet solar maximum

Median values of ionosonde h′F data acquired at Ibadan (Geographic:7.4°N, 3.9°E, Magnetic: dip 6°S, and magnetic declination, 3°W), Nigeria, West Africa, have been used to determine vertical ion drift (electric field) characteristics in the postsunset ionosphere in the African region during a time of high solar activity (average F_10.7 −208). The database spans from January and December 1958 during the era of International Geophysical Year (IGY) for geomagnetic quiet conditions. Bimonthly averaged diurnal variations patterns are very similar, but differ significantly in magnitude and in the evening reversal times. Also, monthly variations of F-region vertical ion drift reversal times inferred from the time of h′F maximum indicates early reversal during equinoxes and December solstice months except for the month of April. Late reversal is observed during the June solstice months. The equatorial evening prereversal enhancement in vertical ion drift (V_zp) occurs largely near 1900 LT with typical values 20–45 m/s. Comparison of Ibadan ionosonde V_zp with the values of prereversal peak velocity reported for Jicamarca (South America), Kodaikanal (India), and Scherliess and Fejer global model show considerable disparity. The changes of postsunset peak in virtual height of F-layer (h′F_P) with prereversal velocity peak V_zp are anti-correlated. Investigation of solar effects on monthly values of V_zp and h′F_P revealed that these parameters are independent of monthly averaged solar flux intensity during quiet-time sunspot maximum conditions.     

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.     

Variation of monthly mean foF2 and hmF2 over a mid latitude station during the period 1997–2006

Ionosonde data of a mid latitude station Novosibirsk (Geog. Lat. 54.6°N, Geog. Long. 83.2°E) has been analyzed for the years 1997–2006 that covers the major part of solar cycle 23. Our results show the presence of winter anomaly in the daytime F2 layer critical frequency during different phases of solar activity. Results also reveal a semiannual variation of foF₂ with two maxima and a minimum that always appears in summer. While the first maximum is in the spring equinox, the second one is found to shift from autumn to winter with the increase of solar activity. The maximum height of F2 layer during the daytime shows variation with the solar activity. It is higher during the higher activity periods and lower during the periods of low activity. Results of ionosonde observations have been compared with those obtained from IRI-2007 model and it is found that model reproduces gross features of foF₂ variation. However, the modeled hmF₂ variations during equinoxes are significantly different from the ones derived using the ionosonde data. The model also underestimates the hmF₂ values.     

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.     

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.