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.     

Cross-hemisphere comparison of mid-latitude ionospheric variability during 1996–2009: Juliusruh vs. Hobart

Analysis of a long-time series of hourly median characteristics of the ionospheric plasma at two mid-latitude locations in the Northern and Southern hemisphere, Juliusruh (54.6N; 13.4E) and Hobart (42.9S; 147.3E), reveals patterns of their synchronous and independent variability. We studied timelines of GPS vTEC, ionogram-derived F2-layer peak electron density NmF2, ionospheric equivalent slab thickness τ, and their ratios at two locations during the complete 23rd solar cycle and its following period of the extremely low solar activity in 2008–2009. This study has also involved the comparative analysis of the observed data versus the model predictions by IRI-2012. During the high solar activity in 2000–2002, seasonal variations show a complicated cross-hemisphere behavior influenced by the winter and semi-annual anomalies, with the largest noon-time values of TEC and NmF2 observed around equinoxes. Strength of the winter anomaly in NmF2 was significantly greater at Juliusruh in comparison with Hobart. The winter anomaly in GPS vTEC values was much weaker than in NmF2 for the Northern hemisphere mid-latitudes and was entirely absent at the Southern hemisphere. Cross-hemisphere analysis of the equivalent slab thickness shows its clear seasonal dependence for all levels of solar activity: the day-time maximum τ_max is observed during local summer, whereas the day-time minimum τ_min is observed during local winter. The night-time values of τ were higher compared to the day-time values during the winter and equinox seasons. Comparative model-data study shows rather good IRI performance of the day-time NmF2 for mid-latitudes of both hemispheres and rather noticeable overestimations for the mid-night NmF2 values during high solar activity. Analysis of IRI vTEC demonstrates the model limitations, related with the absence of the plasmaspheric part, and actual demand in a reliable and standard ionosphere–plasmasphere model for analysis of GPS vTEC.     

Comparison of standard TEC models with a Neural Network based TEC model using multistation GPS TEC around the northern crest of Equatorial Ionization Anomaly in the Indian longitude sector during the low and moderate solar activity levels of the 24th solar cycle

The highest Total Electron Content (TEC) values in the world normally occur at Equatorial Ionization Anomaly (EIA) region resulting in largest ionospheric range delay values observed for any potential Space Based Augmentation System (SBAS). Reliable forecasting of TEC is crucial for satellite based navigation systems. The day to day variability of the location of the anomaly peak and its intensity is very large. This imposes severe limitations on the applicability of commonly used ionospheric models to the low latitude regions. It is necessary to generate a mathematical ionospheric forecasting and mapping model for TEC based on physical ionospheric influencing parameters. A model, IRPE-TEC, has been developed based on real time low latitude total electron content data using GPS measurements from a number of stations situated around the northern crest of the EIA during 2007 through 2011 to predict the vertical TEC values during the low and moderate solar activity levels of the 24th solar cycle. This model is compared with standard ionospheric models like International Reference Ionosphere (IRI) and Parameterized Ionospheric Model (PIM) to establish its applicability in the equatorial region for accurate predictions.     

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.     

Accuracy of IRI profiles of ionospheric density and temperatures derived from comparisons to Kharkov incoherent scatter radar measurements

The incoherent scatter radar (ISR) facility in Kharkov, Ukraine (49.6°N, 36.3°E) measures vertical profiles of electron density, electron and ion temperature, and ion composition of the ionospheric plasma up to 1100 km altitude. Acquired measurements constitute an accurate ionospheric reference dataset for validation of the variety of models and alternative measurement techniques. We describe preliminary results of comparing the Kharkov ISR profiles to the international reference ionosphere (IRI), an empirical model recognized for its reliable representation of the monthly-median climatology of the density and temperature profiles during quiet-time conditions, with certain extensions to the storm times. We limited our comparison to only quiet geomagnetic conditions during the autumnal equinoxes of 2007 and 2008. Overall, we observe good qualitative agreement between model and data both in time and with altitude. Magnitude-wise, the measured and modeled electron density and plasma temperatures profiles appear different. We discovered that representation accuracy improves significantly when IRI is driven by observed-averaged values of the solar activity index rather than their predictions. This result motivated us to study IRI performance throughout protracted solar minimum of the 24th cycle. The paper summarizes our observations and recommendations for optimal use of the IRI.     

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.     

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

Correlation of IRTAM and FPMU data confirming the application of IRTAM to support ISS Program safety

A “Real-Time” plasma hazard assessment process was developed to support International Space Station (ISS) Program real-time decision-making providing solar array constraint relief information for Extravehicular Activities (EVAs) planning and operations. This process incorporates real-time ionospheric conditions, ISS solar arrays’ orientation, ISS flight attitude, and where the EVA will be performed on the ISS. This assessment requires real-time data that is presently provided by the Floating Potential Measurement Unit (FPMU) which measures the ISS floating potential (FP), along with ionospheric electron number density (Ne) and electron temperature (Te), in order to determine the present ISS environment. Once the present environment conditions are correlated with International Reference Ionosphere (IRI) values, IRI is used to forecast what the environment could become in the event of a severe geomagnetic storm. If the FPMU should fail, the Space Environments team needs another source of data which is utilized to support a short-term forecast for EVAs. The IRI Real-Time Assimilative Mapping (IRTAM) model is an ionospheric model that uses real-time measurements from a large network of digisondes to produce foF2 and hmF2 global maps in 15?min cadence. The Boeing Space Environments team has used the IRI coefficients produced in IRTAM to calculate the Ne along the ISS orbital track. The results of the IRTAM model have been compared to FPMU measurements and show excellent agreement. IRTAM has been identified as the FPMU back-up system that will be used to support the ISS Program if the FPMU should fail.     

A new global version of M(3000)F2 prediction model based on artificial neural networks

A new version of global empirical model for the ionospheric propagation factor, M(3000)F2 prediction is presented. Artificial neural network (ANN) technique was employed by considering the relevant geophysical input parameters which are known to influence the M(3000)F2 parameter. This new version is an update to the previous neural network based M(3000)F2 global model developed by Oyeyemi et al. (2007), and aims to address the inadequacy of the International Reference Ionosphere (IRI) M(3000)F2 model (the International Radio Consultative Committee (CCIR) M(3000)F2 model). The M(3000)F2 has been found to be relatively inaccurate in representing the diurnal structure of the low latitude region and the equatorial ionosphere. In particular, the existing hmF2 IRI model is unable to reproduce the sharp post-sunset drop in M(3000)F2 values, which correspond to a sharp post-sunset peak in the peak height of the F2 layer, hmF2. Data from 80 ionospheric stations globally, including a good number of stations in the low latitude region were considered for this work. M(3000)F2 hourly values from 1987 to 2008, spanning all periods of low and high solar activity were used for model development and verification process. The ability of the new model to predict the M(3000)F2 parameter especially in the low latitude and equatorial regions, which is known to be problematic for the existing IRI model is demonstrated.     

Dependence of the F-region peak electron density (foF2) on solar activity observed in the equatorial ionospheric anomaly region in the Brazilian sector

The long-term (solar cycle) changes in the Sun and how it affects the ionospheric F-region observed at São José dos Campos (23.2° S, 45.9° W), Brazil, a location under the southern crest of the equatorial ionospheric anomaly, have been investigated in this paper. The dependence of the F-region peak electron density (foF2) on solar activity during the descending phase of the 23rd solar cycle for the periods of high, medium, and low solar activity has been studied. The ionospheric F-region peak electron densities observed during high and medium solar activity show seasonal variations with maxima close to the equinox periods, whereas during the low solar activity the maxima during the equinox periods is absent. However, during the low solar activity only change observed is a large decrease from summer to winter months. We have further investigated changes in the different ionospheric F-region parameters (minimum virtual height of the F-region (h′F), virtual height at 0.834foF2 (hpF2), and foF2) during summer to winter months in low solar activity periods, 2006–2007 and 2007–2008. Large changes in the two ionospheric parameters (hpF2 and foF2) are observed during summer to winter months in the two low solar activity periods investigated.     

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.     

A prediction model of short-term ionospheric foF2 based on AdaBoost

In this paper, the AdaBoost-BP algorithm is used to construct a new model to predict the critical frequency of the ionospheric F2-layer (foF2) one hour ahead. Different indices were used to characterize ionospheric diurnal and seasonal variations and their dependence on solar and geomagnetic activity. These indices, together with the current observed foF2 value, were input into the prediction model and the foF2 value at one hour ahead was output. We analyzed twenty-two years’ foF2 data from nine ionosonde stations in the East-Asian sector in this work. The first eleven years’ data were used as a training dataset and the second eleven years’ data were used as a testing dataset. The results show that the performance of AdaBoost-BP is better than those of BP Neural Network (BPNN), Support Vector Regression (SVR) and the IRI model. For example, the AdaBoost-BP prediction absolute error of foF2 at Irkutsk station (a middle latitude station) is 0.32 MHz, which is better than 0.34 MHz from BPNN, 0.35 MHz from SVR and also significantly outperforms the IRI model whose absolute error is 0.64 MHz. Meanwhile, AdaBoost-BP prediction absolute error at Taipei station from the low latitude is 0.78 MHz, which is better than 0.81 MHz from BPNN, 0.81 MHz from SVR and 1.37 MHz from the IRI model. Finally, the variety characteristics of the AdaBoost-BP prediction error along with seasonal variation, solar activity and latitude variation were also discussed in the paper.     

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

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

Modeling M(3000)F2 based on extreme learning machine

This paper presents an novel extreme learning machine (ELM)-based prediction model for the ionospheric propagation factor M(3000)F2 at Darwin station (12.4°S, 131.5°E; −44.5°dip) in Australia. The proposed ELM model is trained with hourly daily values of M(3000)F2 from the period 1998–2014 except 2001 and 2009. The hourly daily values of 2001 (high solar activity) and 2009 (low solar activity) are used for validating the prediction accuracy. The proposed ELM for modeling M(3000)F2 can achieve faster training process and similar testing accuracy compared with backward propagation neural network (BPNN). In addition, the performance of the ELM is verified by comparing the predicted values of M(3000)F2 with observed values and the international reference ionosphere (IRI −2016) model predicted values. Based on the error differences (the root mean square error (RMSE) and the M(3000)F2 percentage improvement values M(3000)F2_IMP(%)), the result demonstrates the effectiveness of the ELM model compared with the IRI-2016 model at hourly, daily, monthly, and yearly in high (2001) and low (2009) solar activity years. The ELM also shows good agreement with observations compared with the IRI during disturbed magnetic activity.     

Initial results of the evaluation of IRI hmF2 performance for minima 22–23 and 23–24

The performance of the International Reference Ionosphere (IRI) in predicting the height of the maximum of electron density (hmF2) has been evaluated for similar geomagnetic latitudes stations in the northern hemisphere (NH) and southern hemisphere (SH), and for the last two minima. As truth-sites, the digisonde stations of Millstone Hill (42.6°N, 288.5°E), USA, and Grahamstown (33.3°S, 26.5°E), South Africa, were considered. A monthly averaged diurnal variation was obtained from all the observations and model output in the months studied, and the corresponding difference was also calculated. For this initial study data from summer and winter in the NH and SH were selected for the solstice comparison, and October data for both stations were used to represent equinox conditions. The choice of these periods depended on data availability and quality. The results show that for the earlier minimum in 1996, in general IRI hmF2 values are in reasonable agreement with the observations. The exceptions are October and December in the SH, where IRI hmF2 tends to high, particularly on the dayside, and also July for which the daytime measured values tend to be larger than the IRI ones. For the recent minimum in 2008, IRI tends to over-estimate the hmF2 in most of the observations. The results support the general assertion that thermospheric temperatures were cooler during the last solar minimum as a consequence of an unusually low, and extended, minimum in solar extreme-ultraviolet flux, and in response to continually increasing long-term trend in anthropogenic carbon dioxide. The cooler temperatures not only decrease density at a fixed height, but also make the corresponding contraction of the atmosphere lower the height of the F-region peak.     

A 12 year comparison of MIDAS and IRI 2007 ionospheric Total Electron Content

Data assimilation in conventional meteorological applications uses measurements in conjunction with a physical model. In the case of the ionised region of the upper atmosphere, the ionosphere, assimilation techniques are much less mature. The empirical model known as the International Reference Ionosphere (IRI) could be used to augment data-sparse regions in an ionospheric now-cast and forecast system. In doing so, it is important that it does not introduce systematic biases to the result. Here, the IRI model is compared to ionospheric observations from the Global Positioning System satellites over Europe and North America. Global Positioning System data are processed into hour-to-hour monthly averages of vertical Total Electron Content using a tomographic technique. A period of twelve years, from January 1998 to December 2009, is analysed in order to capture variations over the whole solar cycle. The study shows that the IRI model underestimates Total Electron Content in the daytime at solar maximum by up to 37% compared to the monthly average of GPS tomographic images, with the greatest differences occurring at the equinox. IRI shows good agreement at other times. Errors in TEC are likely due to peak height and density inaccuracies. IRI is therefore a suitable model for specification of monthly averages of Total Electron Content and can be used to initialise a data assimilation process at times away from solar maximum. It may be necessary to correct for systematic deviations from IRI at solar maximum, and to incorporate error estimation into a data assimilation scheme.     

Calibration of IRI–ITU-R peak density and height over the oceans with topside sounding data

Accuracy of IRI electron density profile depends on the F2 layer peak density and height converted by empirical formulae from the critical frequency and M3000F2 factor provided by the ITU-R (former CCIR). The CCIR/ITU-R maps generated from ground-based ionosonde measurements suffer from model assumptions, in particular, over the oceans where relatively few measurements are available due to a scarcity of ground-based ionosondes. In the present study a grid-point calibration of IRI/ITU-R maps for the foF2 and hmF2 over the oceans is proposed using modeling results based on the topside true-height profiles provided by ISIS1, ISIS2, IK-19 and Cosmos-1809 satellites for the period of 1969–1987. Topside soundings results are compared with IRI and the Russian standard model of ionosphere, SMI, and grouped to provide an empirical calibration coefficient to the peak density and height generated from ITU-R maps. The grid-point calibration coefficients maps are produced in terms of the solar activity, geodetic latitude and longitude, universal time and season allowing update of IRI–ITU-R predictions of the F2 layer peak parameters.     

The quiet nighttime low-latitude ionosphere as observed by TIMED/GUVI

In this paper, we examine the nighttime ionosphere climatology structure in the low latitude region and discrepancies between Global Ultraviolet Imager (GUVI) observations and the IRI model predictions using (1) the magnetic zonal mean of electron number density as a function of altitude and magnetic latitude, (2) vertical electron density profiles at various levels of F10.7 index, (3) nighttime descent and magnitude decrease of the ionosphere, (4) point-to-point comparisons of F-peak height (hmF2) and density (NmF2), and (5) the magnetic longitudinal variations of hmF2 and NmF2. The data collected from the Thermosphere, Ionosphere, Mesosphere, Energetics, and Dynamics (TIMED) mission since its launch in December 2001 have provided great opportunities for many scientific investigations of the ionosphere. In this analysis, we investigate the climatology of the nighttime low-latitude ionosphere under low geomagnetic activity (kp ? 4) using the electron density profiles inferred from the airglow measurements obtained by the GUVI aboard the TIMED spacecraft and compared with the results obtained from IRI (International Reference Ionosphere) model-2001. The observed climatology is an essential tool for further understanding the electrodynamics in the low-latitude region and improving the model’s prediction capability. The time range of the GUVI data used in this study is from 2002 (day 053) to 2006 (day 304), and the IRI model predictions were produced at every GUVI location. The ionosphere observed is generally of greater density than what IRI predicts throughout the night for all four seasons for low and moderate solar activity while the model over-predicts the electron density near the F-region peak at high solar activity before midnight. Observations show that the height of the F-region peak has a steep descent from dusk to midnight and near midnight the height of layer is insensitive to solar conditions, significantly different than what is predicted by IRI. Longitudinal features shown in GUVI data are present in the low-latitude ionosphere after sunset and continue through to midnight after which the low-latitude ionosphere is largely zonally symmetric.     

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.