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.     

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.     

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.     

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.     

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.     

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.     

A simple algorithm for regional mid-latitude ionospheric modelling

A simplified model of the key vertical-incidence ionospheric characteristics such as f₀ F2, M(3000) F2, h'F, f₀F1, and f₀E that are used for prediction of operational parameters of HF telecommunication system in restricted area has been developed. The model is based upon the Fourier analysis of the monthly median values of these characteristics from seven ionospheric stations in Europe. It is shown that only 12 dominant Fourier coefficeints, expressed as function of local time, solar activity and geographic latitude are sufficient to reproduce the main features of the diurnal, seasonal and solar cycle behaviour of the mid-latitude ionosphere under median conditions.     

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.     

Evaluation of global modeling of M(3000)F2 and hmF2 based on alternative empirical orthogonal function expansions

Global modeling of M(3000)F2 and hmF2 based on three alternative EOF (empirical orthogonal function) expansion methods is described briefly. Data used for the model construction is the monthly median hourly values of M(3000)F2 from the ionosonde/digisonde stations distributed around the world for the period of 1975–1985 and the hmF2 data of the same period converted from the measured M(3000)F2 based on the strong anti-correlation existing between them. Independent data of a low (1965) and a high (1970) solar activity year are used to validate the three alternative models based on different EOF expansion methods. Comparisons between the modeled results and observed data for both the low (1965) and high (1970) solar activity years showed good agreement for both M(3000)F2 and hmF2 parameters. Statistical analysis on the differences between model values and observed data showed that all the three alternative models (Model A, B and C) based on the different EOF expansion methods have better agreement with the observed data than the models currently used in IRI. All three alternative EOF based models have almost the same accuracy. Discussion on the preference of the three alternative EOF based models is given.     

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.     

Comparative study of foF2 measurements with IRI-2007 model predictions during extended solar minimum

The unusually deep and extended solar minimum of cycle 23/24 made it very difficult to predict the solar indices 1 or 2 years into the future. Most of the predictions were proven wrong by the actual observed indices. IRI gets its solar, magnetic, and ionospheric indices from an indices file that is updated twice a year. In recent years, due to the unusual solar minimum, predictions had to be corrected downward with every new indices update. In this paper we analyse how much the uncertainties in the predictability of solar activity indices affect the IRI outcome and how the IRI values calculated with predicted and observed indices compared to the actual measurements. Monthly median values of F2 layer critical frequency (foF2) derived from the ionosonde measurements at the mid-latitude ionospheric station Juliusruh were compared with the International Reference Ionosphere (IRI-2007) model predictions. The analysis found that IRI provides reliable results that compare well with actual measurements, when the definite (observed and adjusted) indices of solar activity are used, while IRI values based on earlier predictions of these indices noticeably overestimated the measurements during the solar minimum. One of the principal objectives of this paper is to direct attention of IRI users to update their solar activity indices files regularly. Use of an older index file can lead to serious IRI overestimations of F-region electron density during the recent extended solar minimum.     

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.     

Equatorial predictions from a new neural network based global foF2 model

A new neural network (NN) based global empirical model for the foF2 parameter, which represents the peak ionospheric electron density, has been developed using extended temporal and spatial geophysical relevant inputs. It has been proposed that this new model be considered as a suitable replacement for the International Union of Radio Science (URSI) and International Radio Consultative Committee (CCIR) model options currently used within the International Reference Ionosphere (IRI) model for the purpose of F2 peak electron density predictions. The most recent version of the model has incorporated data from 135 global ionospheric stations including a number of equatorial stations.     

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.     

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

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

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.     

Parameterized Regional Ionospheric Model and a comparison of its results with experimental data and IRI representations

We describe a Parameterized Regional Ionospheric Model (PARIM) to calculate the spatial and temporal variations of the ionospheric electron density/plasma frequency over the Brazilian sector. The ionospheric plasma frequency values as calculated from an enhanced Sheffield University Plasmasphere–Ionosphere Model (SUPIM) were used to construct the model. PARIM is a time-independent 3D regional model (altitude, longitude/local time, latitude) used to reproduce SUPIM plasma frequencies for geomagnetic quiet condition, for any day of the year and for low to moderately high solar activity. The procedure to obtain the modeled representation uses finite Fourier series so that all plasma frequency dependencies can be represented by Fourier coefficients. PARIM presents very good results, except for the F region peak height (hmF2) near the geomagnetic equator during times of occurrence of the F₃ layer. The plasma frequency calculated by IRI from E region to bottomside of the F region present latitudinal discontinuities during morning and evening times for both solar minimum and solar maximum conditions. Both the results of PARIM and the IRI for the E region peak density show excellent agreement with the observational values obtained during the conjugate point equatorial experiment (COPEX) campaign. The IRI representations significantly underestimate the foF2 and hmF2 compared to the observational results over the COPEX sites, mainly during the evening–nighttime period.     

Controlling role of maximum usable frequencies in ionospheric informatics

Considering peak parameters of the ionospheric F region which might be observed/modelled/predicted a few relations connecting such data are pointed out. It is pointed out that the sub-peak semi-thickness of the profile depends on the peak height of the F2-layer. The dispersion range of monthly measurement of the MUFs (maximum usable frequency) inferred from M(3000)F2 is compared with that of the critical frequencies foF2 and found to be larger. The inverted Shimazaki formula yields straight-forward values of MUF from the F2 layer critical frequency and the peak height ZmF2. When comparing MUF-3000 values scaled routinely from ionograms with those obtained from foF2 and ZmF2 of the profile analysis, good agreement is obtained. It is felt that the parameter MUF-3000 is a valuable means for checking data obtained by measurements or by models or by predictions.     

Validate the IRI2007 model by the COSMIC slant TEC data during the extremely solar minimum of 2008

During 2008, the solar activity is extremely low. The satellite observations show that the ionospheric height and electron density is much lower than the predictions by the international reference ionosphere (IRI) model. In this paper, we compared the slant total electron content (TEC) observed by the COSMIC satellites during 2008 with the IRI model results. It is found that the IRI model with IRI2001 and IRI2001 Cor. topside options will always overestimate the electron density in both lower and higher altitudes. But the rest two topside options (NeQuick, and TTS) tend to overestimate the electron density in the F layer and underestimate it in the topside altitudes. The switch altitude between overestimation and underestimation and the latitude-local time distribution of the model deviation depend on the topside option. The current investigation might be useful for the model improvement as well as data assimilation work based on the IRI model and the LEO TEC data.     

Ionospheric heating due to solar flares as measured by SROSS-C2 satellite

The data on thermal fluctuations of the topside ionosphere have been measured by Retarding Potential Analyser (RPA) payload aboard the SROSS-C2 satellite over the Indian region for half of the solar cycle (1995–2000). The data on solar flare has been obtained from National Geophysical Data Center (NGDC) Boulder, Colorado (USA) and other solar indices (solar radio flux and sunspot number) were download from NGDC website. The ionospheric electron and ion temperatures show a consistent enhancement during the solar flares. The enhancement in the electron temperature is 28–92% and for ion temperature it is 18–39% compared to the normal day’s average temperature. The enhancement of ionospheric temperatures due to solar flares is correlated with the variation of sunspot and solar radio flux (F10.7cm). All the events studied in the present paper fall in the category of subflare with almost same intensity. The ionospheric electron and ion temperatures enhancement have been compared with the IRI model values.