Equatorial anomaly according to the Interkosmos-19 data and IRI model: A comparison

The comparison of the IRI model with the foF2 distribution in the equatorial anomaly region obtained by topside sounding onboard the Interkosmos-19 satellite has been carried out. The global distribution of foF2 in terms of LT-maps was constructed by averaging Intercosmos-19 data for summer, winter, and equinox. These maps, in fact, represent an empirical model of the equatorial anomaly for high solar activity F10.7 ~ 200. The comparison is carried out for the latitudinal foF2 profiles in the characteristic longitudinal sectors of 30, 90, 210, 270, and 330°, as well as for the longitudinal variations in foF2 over the equator. The largest difference between the models (up to 60%) for any season was found in the Pacific longitudinal sector of 210°, where there are a few ground-based sounding stations. Considerable discrepancies, however, are sometimes observed in the longitudinal sectors, where there are many ground-based stations, for example, in the European or Indian sector. The discrepancies reach their maximum at 00 LT, since a decay of the equatorial anomaly begins before midnight in the IRI model and after midnight according to the Interkosmos-19 data. The discrepancies are also large in the morning at 06 LT, since in the IRI model, the foF2 growth begins long before sunrise. In the longitudinal variations in foF2 over the equator at noon, according to the satellite data, four harmonics are distinguished in the June solstice and at the equinox, and three harmonics in the December solstice, while in the IRI model only two and one harmonics respectively are revealed. In diurnal variations in foF2 and, accordingly, in the equatorial anomaly intensity, the IRI model does not adequately reproduce even the main, evening extremum.     

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.     

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.     

Variation of F-region critical frequency (foF2) over equatorial and low-latitude region of the Indian zone during 19th and 20th solar cycle

The effect of solar cycle and seasons on the daytime and nighttime F-layer ionization has been investigated over the equatorial and low-latitude region during 19th (1954–1964) and 20th (1965–1976) solar cycle. The F-layer critical frequency (foF2) data observed from the three Indian Ionosonde stations has been used for the present study. The dependence of foF2 on solar cycle has been examined by performing regression analysis between the foF2 values and R₁₂ (twelve month running average sunspot number). The result shows that the magnitude of the cycle, seasons and the location of station has considerable effects on foF2. There is a significant nonlinear relationship between the foF2 values and R₁₂ during 19th solar cycle as compared to 20th solar cycle. Further, the nighttime saturation effect is prominently seen during the 19th solar cycle and summer season. It is also observed that the most profound saturation effect appears at the equatorial ionization anomaly crest region. Seasonally, it is seen that all the stations exhibits semiannual anomaly. The phenomenon of winter anomaly decays as we move higher along the latitude and is prominently seen during the intense solar activity.     

Stability of the seasonal variations in diurnal and semidiurnal components of mid-latitude F2 layer parameters

We report work utilizing 15-min resolution ionospheric data obtained with DPS-4 digisonde in 2003–2011 to study the seasonal variations in amplitudes and phases of the most powerful spectral components of the F2 layer critical frequency (foF2) and peak height (hmF2) fluctuations over Irkutsk (52.5°N, 104.0°E). We show that fluctuations of both parameters contain quasi-harmonic components with periods of T_n = 24/n h (n = 1–7). The number of distinct spectral peaks varies from 3 in summer to 7 in winter. Amplitude and phase characteristics of the diurnal (n = 1) and semidiurnal (n = 2) components is studied using the data sets extracted from the original data sets with band-pass filter. It has been found that the amplitudes of diurnal/semidiurnal foF2 and diurnal hmF2 components are maximum in winter and minimum in summer. Amplitudes of the diurnal components vary gradually; those of the foF2 semidiurnal one, abruptly, thus forming a narrow winter maximum in November–January. The phase (local time of maximum) of the diurnal foF2 component increases gradually by 4–6 h from winter to summer. The phase of the semidiurnal foF2 component is nearly stable in winter/summer and sharply decreases (increases) by 2–3 h near the spring (autumn) equinox. The phase of the diurnal component of hmF2 (local time of minimum) varies slightly between 1130 and 1300 LT; that of the semidiurnal one decreases (increases) by 4–6 h from January to March (from September to November). The results obtained show that the main features of seasonal variations in the diurnal and semidiurnal components of the mid-latitude F2 layer parameters recur consistently during the solar activity growth and decline phases.     

Responses of ionospheric foF2 to geomagnetic activities in Hainan

Responses of low-latitude ionospheric critical frequency of F2 layer to geomagnetic activities in different seasons and under different levels of solar activity are investigated by analyzing the ionospheric foF2 data from DPS-4 Digisonde in Hainan Observatory during 2002–2005. The results are as follows: (1) the response of foF2 to geomagnetic activity in Hainan shows obvious diurnal variation except for the summer in low solar activity period. Generally, geomagnetic activity will cause foF2 to increase at daytime and decrease at nighttime. The intensity of response of foF2 is stronger at nighttime than that at daytime; (2) seasonal dependence of the response of foF2 to geomagnetic activity is very obvious. The negative ionospheric storm effect is the strongest in summer and the positive ionospheric storm effect is the strongest in winter; (3) the solar cycle has important effect on the response of foF2 to geomagnetic activity in Hainan. In high solar activity period, the diurnal variation of the response of foF2 is very pronounced in each season, and the strong ionospheric response can last several days. In low solar activity period, ionospheric response has very pronounced diurnal variation in winter only; (4) the local time of geomagnetic activities occurring also has important effect on the responses of foF2 in Hainan. Generally, geomagnetic activities occurred at nighttime can cause stronger and longer responses of foF2 in Hainan.     

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.     

A new global F2 peak electron density model for the International Reference Ionosphere (IRI)

A new neural network (NN) based global empirical model for the F2 peak electron density (NmF2) has been developed using extended temporal and spatial geophysical relevant inputs. Measured ground based ionosonde data, from 84 global stations, spanning the period 1995 to 2005 and, for a few stations from 1976 to 1986, obtained from various resources of the World Data Centre (WDC) archives (Space Physics Interactive Data Resource SPIDR, the Digital Ionogram Database, DIDBase, and IPS Radio and Space Services) have been used for training a NN. The training data set includes all periods of quiet and disturbed magnetic activity. A comprehensive comparison for all conditions (e.g., magnetic storms, levels of solar activity, season, different regions of latitudes, etc.) between foF2 value predictions using the NN based model and International Reference Ionosphere (IRI) model (including both the International Union of Radio Science (URSI) and International Radio Consultative Committee (CCIR) coefficients) with observed values was investigated. The root-mean-square (RMS) error differences for a few selected stations are presented in this paper. The results of the foF2 NN model presented in this work successfully demonstrate that this new model can be used as a replacement option for the URSI and CCIR maps within the IRI model for the purpose of F2 peak electron density predictions.     

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.     

Temporal and spatial variation of equatorial ionization anomaly by using multistation ionosonde data for the 19th solar cycle over the Indian region

The paper deals with the study of temporal and spatial variation of equatorial ionization anomaly (EIA) phenomenon along with its dependence on solar activity and season during the 19th solar cycle by using seven Indian ionosonde stations. Present study is an attempt to carry out the comprehensive study of EIA by using the limited number of ground based instruments. This has been achieved by performing the Gaussian fitting over the latitudinal distribution of F2-region critical frequency (foF2) data. Results reveal that the phenomenon of EIA has a strong dependence on solar activity and seasons. The EIA crest exhibits the feature of latitudinal shifting and expansion with increasing solar activity. It is found out that the effect of solar cycle and seasons on EIA is local time dependent. The observations were also compared with the IRI-2001 model predictions and results reveal that the model values are in general agreement with the observed values with some discrepancies, particularly during the high solar activity period and morning sector. The results have been discussed in the light of relative contribution from transequatorial interhemispheric neutral wind and strength of equatorial fountain process during different local time, season and solar activity levels. Furthermore, an attempt is made to parameterize the location and foF2 of the EIA crest by using the regression analysis. These results can be used to predict the latitudinal position and foF2 of the EIA crest for any given 12-month running average sunspot number (R12).     

Latitudinal variation of F2-region response to geomagnetic disturbance

The effect of geomagnetic storms on the F2 region was studied by calculating the deviation, ΔfoF2, of foF2 during 40 magnetic storms, ranging from moderate (Dst < −50 nT) to very intense (Dst < −200 nT) of the 21st solar cycle. In order to study the variation of storm-time foF2 with latitude, season and storm strength, ionosonde data were obtained from eight stations spanning a latitudinal range of +60–−60°. The stations chosen lay in a narrow longitudinal range of 140–151°, so that local time difference between the stations is practically negligible. The features exhibited by positive and negative phases were essentially different. The storm time ΔfoF2 clearly exhibited a latitudinal variation and this variation were found to be coupled with the seasonal variation. As for the variation with storm intensity, though ΔfoF2 was found to vary even between two storms of almost equal intensity, the amplitude of a positive or negative phase, |ΔfoF2max| showed a distinct upper limit for each intensity category of storms.     

Monthly mean foF2 model for an African low-latitude station and comparison with IRI

We have employed the hourly values of the ionospheric F-region critical frequency (foF2) obtained from Ouagadougou ionosonde, Burkina Faso (geographic coordinates 12° N, 1.8° W) during the interval of 1985–1995 (solar cycle 22) and solar radio flux of 10 cm wavelength (F10.7) to develop a local model (LM) for the African low-latitude station. The model was developed from regression analysis method, using the two-segmented regression analysis. We validated LM with foF2 data from Korhogo observatory, Cote d’Ivorie (geographical coordinates 9.3° N, 5.4° W). LM as well as the International Reference Ionosphere (IRI) agrees well with observations. LM gave some improvement on the IRI-predicted foF2 values at the sunrise (06 LT) at all solar flux levels and in all seasons except June solstice. The performance of the models at the representing the salient features of the equatorial foF2 was presented. Considering daytime and nighttime performances, LM and IRI are comparable in low solar activity (LSA), LM performed better than IRI in moderate solar activity (MSA), while IRI performed better than LM in high solar activity (HSA). CCIR has a root mean square error (r.m.s.e), which is only 0.10 MHz lower than that of LM while LM has r.m.s.e, which is about 0.05 MHz lower than that of URSI. In general, our result shows that performance of IRI, especially the CCIR option of the IRI, is quite comparable with the LM. The improved performance of IRI is a reflection of the numerous contributions of ionospheric physicists in the African region, larger volume of data for the IRI and the diversity of data sources, as well as the successes of the IRI task force activities.     

Peculiarities of the nighttime winter foF2 increase over Irkutsk

An analysis of properties and peculiarities of the nighttime winter foF2 increases (NWI) in the East Siberia is made on data of ionospheric station Irkutsk in the periods 1958–1992 and 2002–2009 and the empirical model of the F2 layer critical frequency under the geomagnetic quiet conditions deduced from these data (model Q-F2). It is revealed, that the NWI is the stable regularity of the quiet ionosphere over Irkutsk. The amplitude of the NWI (the difference between maximum and minimum foF2 values at night hours) is the greatest in December–January and nearly the same at low and middle solar activity. It is a peculiarity of the quiet ionosphere in the East Siberia. Maximum in night foF2 under quiet geomagnetic conditions is observed mainly after midnight (02-04 LT) and is shifted to predawn hours as solar activity increases. At low solar activity the quiet ionosphere at ∼02–04 LT shows the following properties: (a) the fluctuations of foF2 and hmF2 are in the reverse correlation but this dependence is weak; (b) very strong fluctuations of foF2 (|δfoF2| > 30%) occur seldom (∼4% of events) and almost all of them are positive; an example of very strong fluctuations of foF2 up to 60% can be an extreme increase in the foF2 on 19.12.2008; (c) the very strong enhancements of foF2 in the NWI maximum can be observed at the low geomagnetic activity, they occur more often during substorms but very seldom during geomagnetic storms. Possible reasons of these properties of NWI are discussed.     

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.     

Wave-like perturbations in the ionospheric F2-layer observed after the Ms8.1 Samoa earthquake of September 29, 2009

A study of the critical frequency foF2 variations after the large earthquake (M_s = 8.1) which occurred on 29 September, 2009 in the region of Samoa Islands in the Pacific Ocean is carried out using data of the ionospheric station of Kwajalein. The epicenter of the earthquake was located at about 184 km southwest from Apia (the capital of West Samoa). It was found that wave-like perturbations of foF2 were observed for ∼3 h above the station (located approximately 3560 km northwest from the epicenter). The amplitude of the disturbance was as large as ∼20% of the average magnetic quiet day foF2 values. A comparison of the observed perturbations of foF2 with the ones detected at Stanford ionospheric station after the Alaska earthquake of 28 March 1964 (M_s = 8.4) showed a close similarity of the wave-like perturbations of foF2 in both cases.     

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.     

Longitudinal variation of the equatorial ionosphere: Modeling and experimental results

We describe a new version of the Parameterized Regional Ionospheric Model (PARIM) which has been modified to include the longitudinal dependences. This model has been reconstructed using multidimensional Fourier series. To validate PARIM results, the South America maps of critical frequencies for the E (foE) and F (foF2) regions were compared with the values calculated by Sheffield Plasmasphere-Ionosphere Model (SUPIM) and IRI representations. PARIM presents very good results, the general characteristics of both regions, mainly the presence of the equatorial ionization anomaly, were well reproduced for equinoctial conditions of solar minimum and maximum. The values of foF2 and hmF2 recorded over Jicamarca (12°S; 77°W; dip lat. 1°N; mag. declination 0.3°) and sites of the conjugate point equatorial experiment (COPEX) campaign Boa Vista (2.8°N; 60.7°W; dip lat. 11.4°; mag. declination −13.1°), Cachimbo (9.5°S; 54.8°W; dip lat. −1.8°; mag. declination −15.5°), and Campo Grande (20.4°S; 54.6°W; dip lat. −11.1°; mag. declination −14.0°) have been used in this work. foF2 calculated by PARIM show good agreement with the observations, except during morning over Boa Vista and midnight-morning over Campo Grande. Some discrepancies were also found for the F-region peak height (hmF2) near the geomagnetic equator during times of F₃ layer occurrences. IRI has underestimated both foF2 and hmF2 over equatorial and low latitude sectors during evening-nighttimes, except for Jicamarca where foF2 values were overestimated.     

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.     

F-region ionospheric parameters observed in the equatorial and low latitude regions during medium solar activity in the Brazilian sector and comparison with the IRI-2007 model results

The ionospheric sounding observations using the Canadian Advanced Digital Ionosondes (CADIs) operational at Palmas (PAL; 10.2°S, 48.2°W; dip latitude 6.6°S; a near-equatorial station), and São José dos Campos (SJC, 23.2°S, 45.9°W; dip latitude 17.6°S; a low-latitude station located under the southern crest of the equatorial ionospheric anomaly), Brazil, are analyzed during the different seasons viz., winter (June and July 2003), spring (September and October 2003), summer (December 2003 and January 2004), and fall (March and April 2004). The period used has medium solar activity (sunspot number between 77.4 and 39.3). The seasonal mean variations (using only geomagnetically quiet days) of the ionospheric parameters foF2 (critical frequency of the F-region), hpF2 (virtual height at 0.834 foF2; considered to be close to hmF2 (peak height of the F-region)), and h’F (minimum virtual height of the F-region) are calculated and compared between PAL and SJC. The prominent differences between PAL and SJC are as follows: h’F variations show strong post-sunset enhancement at PAL during the seasons of spring, summer, and fall; hpF2 variations show pre-sunrise uplifting of the F-layer at both stations during all the seasons and the hpF2 values during the daytime are lower at SJC compared with PAL during all the seasons; the foF2 variations show mid-day bite-out at PAL during all the seasons and SJC shows strong equatorial ionospheric anomaly during summer and fall seasons. Also, the seasonal variations of the ionospheric parameters foF2 and hpF2 (with ±1 standard deviation) observed at PAL and SJC are compared with the IRI-2007 model results of foF2 and hmF2. In addition, variations of the foF2 and hpF2 observed at SJC are compared with the IRI-2001 model results of foF2 and hmF2. It should be pointed out that the ionospheric parameter hpF2 is much easier to obtain using computer program developed at UNIVAP compared with hmF2 (using POLAN program). During the daytime due to underlying ionization hpF2 estimated is higher (approximately 50 km) than the true peak height hmF2. During the nighttime hpF2 is fairly close to hmF2. The comparison between the foF2 variations observed at PAL and SJC with the IRI-2007 model results shows a fairly good agreement during all the seasons. However, the comparison between the hpF2 variations observed at PAL and SJC with the hmF2 variations with the IRI-2007 model results shows: (1) a fairly good agreement during the nighttime in all the seasons; (2) the model results do not show the pre-sunrise uplifting of the F-layer at PAL and SJC in any season; (3) the model results do not show the post-sunset uplifting of the F-layer at PAL; (4) considering that, in general, hpF2 is higher than hmF2 during the daytime by about 50 km, the model results are in good agreement at PAL and SJC during all the seasons except summer at SJC, when large discrepancies in the observed hpF2 and modeled hmF2 are observed. Also, it has been observed that, in general, hmF2 values for SJC calculated using IRI-2001 are higher than IRI-2007 during the daytime in winter, summer, and fall. However, hmF2 values for SJC calculated using IRI-2001, are lower than IRI-2007 during the nighttime in spring.