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.     

Surface latent heat flux anomalies quasi-synchronous with ionospheric disturbances before the 2007 Pu’er earthquake in China

One of various mechanisms of pre-earthquake lithosphere–atmosphere–ionosphere coupling as possible explanation of the seismo-ionospheric effects is connected with the release of latent heat. Abnormal variations of ionospheric electromagnetic parameters possibly related to the 2007 Ms 6.4 Pu’er earthquake in China were reported. This paper attempts to examine whether there were abnormal changes of surface latent heat flux (SLHF) linked with this pre-earthquake ionospheric disturbances. The spatio-temporal statistical analyzes of multi-years SLHF data from USA NCEP/NCAR Reanalysis Project reveal that local SLHF enhancements appeared 11, 10 and 7 days before the Pu’er earthquake, respectively. As contrasted to the formerly reported local ionospheric Ne enhancement 9 days before the shock observed by DEMETER satellite, it is discovered that the SLHF and Ne anomalies are quasi-synchronous and have good spatial correspondence with the epicentre and the local active faults. This is valuable for understanding the seismogenic coupling processes and for recognizing earthquake anomaly with multiple parameters from integrated Earth observation system.     

Statistical analysis of the ionospheric ion density recorded by DEMETER in the epicenter areas of earthquakes as well as in their magnetically conjugate point areas

Results of a statistical variation of total ion density observed in the vicinity of epicenters as well as around magnetically conjugated points of earthquakes are presented in this paper. Two data sets are used: the ion density measured by DEMETER during about 6.5?years and the list of strong earthquakes (M_W?≥?4.8) occurring globally during this period (14,764 earthquakes in total). First of all, ionospheric perturbations with 23–120?s observation time corresponding to spatial scales of 160–840?km are automatically detected by a software (64,287 anomalies in total). Second, it is checked if a perturbation could be associated either with the epicenter of an earthquake or with its magnetically conjugated point (distance?<?1500?km and time?<?15?days before the earthquake). The index Kp?<?3 is also considered in order to reduce the effect of the geomagnetic activity on the ionosphere during this period. The results show that it is possible to detect variations of the ionospheric parameters above the epicenter areas as well as above their conjugated points. About one third of the earthquakes are detected with ionospheric influence on both sides of the Earth. There is a trend showing that the perturbation length increases as the magnitude of the detected EQs but it is more obvious for large magnitude. The probability that a perturbation appears is higher on the day of the earthquake and then gradually decreases when the time before the earthquake increases. The spatial distribution of perturbations shows that the probability of perturbations appearing southeast of the epicenter before an earthquake is a little bit higher and that there is an obvious trend because perturbations appear west of the conjugated point of an earthquake.     

Study of the 2013 Lushan M7.0 earthquake coseismic ionospheric disturbances

On April 20, 2013, an earthquake of M7.0 occurred in Lushan, Sichuan province, China. This paper investigates the coseismic ionospheric anomalies using GPS (Global Positioning System) data from 23 reference stations in Sichuan province that are a part of the Crustal Movement Observation Network of China (CMONOC). The recorded results show that a clear ionospheric anomaly occurred within 15 min after the earthquake near the epicenter, and the occurrence time of the anomalies recorded by various stations is related to the distance from the epicenter. The maximum anomaly is 0.25 TECu, with a 2 min duration and the distance of the recording station to the epicenter is 83 km. Acoustic waves generated by the crustal vertical movement during the earthquake propagate up to the height of the ionosphere lead to the ionospheric anomaly, and the propagation speed of the acoustic wave is calculated as 0.72 ± 0.04 km/s based on the propagation time and propagation distance, consistent with the average speed of sound waves within a 0–450 km atmospheric height.     

Ionospheric disturbances on 8 September, 2010: Was it connected with the incoming moderate Chongqing earthquake?

Earthquake prediction stimulates the searches for a correlation between seismic activity and ionospheric anomalies. Contrary to common focuses on strong earthquakes, we report the ionospheric disturbances, 2 days before a moderate Ms = 4.7 Chongqing earthquake (29.4°N, 105.5°E, depth = 7.0 km, occurred at 21:21 LT, 10 September, 2010) with the data of ground-based ionosondes and IGS receivers. The data covering the period under the quiet geomagnetic conditions and a geomagnetic storm was analyzed with upper and lower bounds. It is found that there were significant enhancements of foF2 and total electron content (TEC) on the afternoon of 8 September, 2010, with a limited area close to the epicentre, which was different from the feature of ionospheric perturbations triggered by the geomagnetic storm on 15 September. Taking into account the heliogeomagnetical condition, we conclude that the observed ionospheric enhancements were very likely associated with the forthcoming moderate Chongqing earthquake, which implies that the relationship between the amplitudes of ionospheric disturbances and earthquakes is very complicated.     

Magnetospheric disturbances associated with the 13 December 2006 solar flare and their ionospheric effects over North-East Asia

We present an observational study of magnetospheric and ionospheric disturbances during the December 2006 intense magnetic storm associated with the 4В/Х3.4 class solar flare. To perform the study we utilize the ground data from North–East Asian ionospheric and magnetic observatories (60–72°N, 88–152°E) and in situ measurements from LANL, GOES, Geotail and ACE satellites. The comparative analysis of ionospheric, magnetospheric and heliospheric disturbances shows that the interaction of the magnetosphere with heavily compressed solar wind and interplanetary magnetic field caused the initial phase of the magnetic storm. It was accompanied by the intense sporadic E and F2 layers and the total black-out in the nocturnal subauroral ionosphere. During the storm main phase, LANL-97A, LANL 1994_084, LANL 1989-046 and GOES_11 satellites registered a compression of the dayside magnetosphere up to their orbits. In the morning–noon sector the compression was accompanied by an absence of reflections from ionosphere over subauroral ionospheric station Zhigansk (66.8°N, 123.3°E), and a drastic decrease in the F2 layer critical frequency (foF2) up to 54% of the quite one over subauroral Yakutsk station (62°N, 129.7°E). At the end of the main phase, these stations registered a sharp foF2 increase in the afternoon sector. At Yakutsk the peak foF2 was 1.9 time higher than the undisturbed one. The mentioned ionospheric disturbances occurred simultaneously with changes in the temperature, density and temperature anisotropy of particles at geosynchronous orbit, registered by the LANL-97A satellite nearby the meridian of ionospheric and magnetic measurements. The whole complex of disturbances may be caused by radial displacement of the main magnetospheric domains (magnetopause, cusp/cleft, plasma sheet) with respect to the observation points, caused by changes in the solar wind dynamic pressure, the field of magnetospheric convection, and rotation of the Earth.     

Ground-based and satellite DC-ULF electric field anomalies around Wenchuan M8.0 earthquake

The ground-based and satellite DC-ULF electric field data were analyzed around Wenchuan M8.0 earthquake on May 12, 2008 in China. The results show that ground electric field anomalies occurred at 3 stations located to the north and south of the epicenter with the amplitude of 3–100 mV/km. The change shapes and their amplitude of ground electric field anomalies are different largely due to their individual underground layer conductivity, water level and so on. The analysis of long time series illustrates that the abnormal geoelectric field started since March 2008. Onboard the DEMETER satellite, the ULF waveforms of electric field were collected and processed by wavelet transform method. The disturbances in the ionosphere were about 3–5 mV/m at a frequency band lower than 0.5 Hz. When the ground and space electric field anomalies were compared, their occurrence time and spatial distribution points are consistent with each other, including the long time anomalies from March 2008 and the short term ones 1–2 days before the Wenchuan earthquake. Finally, the coupling mechanism was discussed.     

Statistical seismo-ionospheric precursors of M7.0+ earthquakes in Circum-Pacific seismic belt by GPS TEC measurements

The Circum-Pacific seismic belt is the region heavily affected by earthquakes in the world. The relationship between earthquake (e.g., the geographic location, occurrence time, magnitude, and focal depth) and ionospheric anomalies in the belt was investigated using 100 M7.0+ earthquakes during 2006–2015. The ground-based GPS measurements and global ionosphere map (GIM) data were used for the analyses of the ionospheric variations preceding the earthquakes. The results indicated that the occurrence rate of total electron content (TEC) anomalies was proportional to the magnitude and inversely proportional to the focal depth to a certain degree, and the occurrence frequency of anomalies had a rising trend with the days getting close to the main shock. The occurrence rate of TEC anomalies in the Southern hemisphere was larger than that in the Northern hemisphere. Besides, the spatial characteristics of TEC anomalies showed that the anomalies in low-middle latitudes did not coincide with the epicenter, sometimes the anomalies were also observed in the corresponding conjugated region. However, the TEC anomalies in the high latitude usually appeared around the epicenter and within the seismogenic zone while no TEC anomalies appeared in the conjugated area. These results may have potential applications to the earthquake prediction in the Circum-Pacific seismic belt.     

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.     

地磁场与电离层异常现象及其与地震的关系

利用中国地磁台网与电离层台站资料,总结了大地震前出现的地磁低点位移、地磁日变异常及电离层f0F2(F2层临界频率)异常现象.对比研究了1997年11月8日玛尼7.5级与2001年11月14日昆仑山口西8.1级地震前磁场与电离层异常分布及特征.结果显示,两次巨大地震前磁场与电离层短临异常时空分布特征有较好的一致性,震中周围出现日变异常、拉萨台出现电离层f0F2明显异常;震前约1个月出现地磁低点位移,其突变分界线通过震中地区.      

An annual variation analysis of the ionospheric spatial gradient over a regional area for GNSS applications

An ionospheric spatial gradient represents the ionosphere delay difference between different locations, and its variation over a specific area is important for implementing differential GNSS systems. An estimation method for the ionospheric spatial gradient over a small regional area is proposed. A plate map model is implemented for the direct estimation of the gradients. Nine years of GPS data were processed to figure out the annual variation of the mean gradient at the mid-geomagnetic latitude of 30° N. Gradients along the north–south direction have a mean of 0.65 mm/km and follow solar-cycle variations.     

Ionospheric anomalies of local earthquakes detected by GPS TEC measurements using data from Tashkent and Kitab stations

This paper reports the ionospheric anomalies observed during strong local earthquakes (M?5.0) which occurred mostly in and around Uzbekistan in seismically active zones, during years 2006 to 2009 within approximately 1000 km distance from the observing GPS stations located in Tashkent and Kitab, Uzbekistan. The solar and geomagnetic conditions were quiet during occurrence of the selected strong earthquakes. We produce Total Electron Content (TEC) time series over both sites and apply them to detect anomalous TEC signals preceding or accompanying the local earthquakes. The results show anomalous increase or decrease of TEC before or during the earthquakes. In general the anomalies occurred 1–7 days before the earthquakes as ionospheric electromagnetic precursors. To identify the anomalous values of TEC we calculated differential TEC (dTEC). dTEC is obtained by subtracting monthly averaged diurnal vTEC from the values of observed vTEC at each epoch. This procedure removes normal diurnal variations of vTEC. The present results are in good agreement with the previous observations on ionospheric earthquake precursors reported by various researchers.     

Symmetry and asymmetry of ionospheric weather at magnetic conjugate points for two midlatitude observatories

Variations of the ionospheric weather W-index for two midlatitude observatories, namely, Grahamstown and Hermanus, and their conjugate counterpart locations in Africa are studied for a period from October 2010 to December 2011. The observatories are located in the longitude sector, which has consistent magnetic equator and geographic equator so that geomagnetic latitudes of the line of force are very close to the corresponding geographic latitudes providing opportunity to ignore the impact of the difference of the gravitational field and the geomagnetic field at the conjugate points on the ionosphere structure and dynamics. The ionosondes of Grahamstown and Hermanus provide data of the critical frequency (foF2), and Global Ionospheric Maps (GIM) provide the total electron content (TECgps) along the magnetic field line up to the conjugate point in the opposite hemisphere. The global model of the ionosphere, International Reference Ionosphere, extended to the plasmasphere altitude of 20,200 km (IRI-Plas) is used to deliver the F2 layer peak parameters from TECgps at the magnetic conjugate area. The evidence is obtained that the electron gas heated by day and cooled by night at the summer hemisphere as compared with the opposite features in the conjugate winter hemisphere testifies on a reversal of plasma fluxes along the magnetic field line by the solar terminator. The ionospheric weather W-index is derived from NmF2 (related with foF2) and TECgps data. It is found that symmetry of W-index behavior in the magnetic conjugate hemispheres is dominant for the equinoxes when plasma movement along the magnetic line of force is imposed on symmetrical background electron density and electron content. Asymmetry of the ionospheric storm effects is observed for solstices when the plasma diffuse down more slowly into the colder winter hemisphere than into the warmer summer hemisphere inducing either plasma increase (positive phase) or decrease (negative phase of W-index) in the ionospheric and plasmaspheric plasma density.     

Pre-earthquake ionospheric anomalies before three major earthquakes by GPS-TEC and GIM-TEC data during 2015–2017

The dual-frequency satellite-based measurements from Global Positioning System (GPS) may provide feasible ways of studying and potentially detecting of earthquake (EQ) related anomalies in the ionosphere. In this paper, GPS based Total Electron Content (TEC) data are studied for three major M?>?7.0 EQs in Nepal and Iran-Iraq border during 2015–2017 by implementing statistical procedures on temporal and spatial scale. Previous studies presented different time interval of pre-seismic ionospheric anomalies, however, this study showed that EQs ionospheric precursors may occur within 10?days. Furthermore, the ionospheric anomalies on the suspected day occurred during UT?=?08:00–12:00?h before the main shock. The Global Ionospheric Map TEC (GIM-TEC) data retrieved over the epicenter of M7.8 (Nepal EQ) showed a significant increase of 6 TECU on April 24, 2015 (one day before the main shock), which is recorded by the ground GPS station data of Islamabad (station lies within the EQ preparation zone). Furthermore, the spatial GIM-TEC result imply significant anomalies over the epicenter during the time interval between UT?=?08:00–12:00?h (LT?=?13:00–17:00). For M7.3 (Nepal EQ), the TEC anomalies were detected on May 10, 2015 (2?days before the event) in the temporal data. The spatial TEC data imply the huge clouds over the epicenter at about UT?=?08:00–12:00?h on May 10, 2015, that may be associated with this EQ in the quiet geomagnetic storm conditions. Similarly, temporal and spatial TEC showed anomaly on November 3, 2017, during UT?=?08:00–12:00 (9?days before the Iran-Iraq border EQ) after implementing the statistical method on it. Conversely, there exists a short-term but low magnitude TEC anomaly synchronized with a geomagnetic storm on November 7–8, 2017 (4 to 5?days prior to M7.3 Iran-Iraq border EQ). The diurnal and hourly GIM-TEC and VTEC data also imply the execution of ionospheric anomalies within 10?days prior to all events. All these positive anomalies in TEC may be due to the existence of a huge energy from the epicenter during the EQ preparation period.     

GPS derived TEC and foF2 variability at an equatorial station and the performance of IRI-model

The ionosphere induces a time delay in transionospheric radio signals such as the Global Positioning System (GPS) signal. The Total Electron Content (TEC) is a key parameter in the mitigation of ionospheric effects on transionospheric signals. The delay in GPS signal induced by the ionosphere is proportional to TEC along the path from the GPS satellite to a receiver. The diurnal monthly and seasonal variations of ionospheric electron content were studied during the year 2010, a year of extreme solar minimum (F10.7 = 81 solar flux unit), with data from the GPS receiver and the Digisonde Portable Sounder (DPS) collocated at Ilorin (Geog. Lat. 8.50°N, Long. 4.50°E, dip −7.9°). The diurnal monthly variation shows steady increases in TEC and F2-layer critical frequency (foF2) from pre-dawn minimum to afternoon maximum and then decreases after sunset. TEC show significant seasonal variation during the daytime between 0900 and 1900 UT (LT = UT + 1 h) with a maximum during the March equinox (about 35 TECU) and minimum during the June solstice (about 24 TECU). The GPS-TEC and foF2 values reveal a weak seasonal anomaly and equinoctial asymmetry during the daytime. The variations observed find their explanations in the amount of solar radiation and neutral gas composition. The measured TEC and foF2 values were compared with last two versions of the International Reference Ionosphere (IRI-2007 and IRI-2012) model predictions using the NeQuick and CCIR (International Radio Consultative Committee) options respectively in the model. In general, the two models give foF2 close to the experimental values, whereas significant discrepancies are found in the predictions of TEC from the models especially during the daytime. The error in height dependent thickness parameter, daytime underestimation of equatorial drift and contributions of electrons from altitudes above 2000 km have been suggested as the possible causes.     

Deformation of the ionosphere structure during the space weather events on October–November 2003

We present the spatial maps of the ionosphere–plasmasphere slab thickness τ (ratio of the vertical total electron content, TEC, to the F-region peak electron density, NmF2) during the intense ionospheric storms of October–November 2003. The model-assisted technology for estimate of the upper boundary of the ionosphere, hup, from the slab thickness components in the bottomside and topside ionosphere – eliminating the plasmasphere contribution of τ – is applied at latitudes 35° to 70°N and longitudes −10° to 40°E, from the data of 20 observatories of GPS-TEC and ionosonde networks, for selected days and hours of October and November 2003. The daily–hourly values of NmF2, hmF2 and TECgps are used as the constrained parameters for the International Reference Ionosphere extended to the plasmasphere, IRI-Plas, during the ionospheric quiet days, positive and negative storm phases for estimate of τ and hup. Good correlation has been found between the slab thickness and the upper boundary of the ionosphere for the intense ionospheric storms at October–November 2003. During the negative phase of the ionospheric storm, when the ionospheric plasma density is exhausted, the nighttime upper boundary of the ionosphere is greatly uplifted towards the magnetosphere tail, while the daytime upper boundary of the ionosphere is reduced below 500 km over the Earth.     

Ionospheric disturbances under low solar activity conditions

The paper is focused on ionospheric response to occasional magnetic disturbances above selected ionospheric stations located at middle latitudes of the Northern and Southern Hemisphere under extremely low solar activity conditions of 2007–2009. We analyzed changes in the F2 layer critical frequency foF2 and the F2 layer peak height hmF2 against 27-days running mean obtained for different longitudinal sectors of both hemispheres for the initial, main and recovery phases of selected magnetic disturbances. Our analysis showed that the effects on the middle latitude ionosphere of weak-to-moderate CIR-related magnetic storms, which mostly occur around solar minimum period, could be comparable with the effects of strong magnetic storms. In general, both positive and negative deviations of foF2 and hmF2 have been observed independent on season and location. However positive effects on foF2 prevailed and were more significant. Observations of stormy ionosphere also showed large departures from the climatology within storm recovery phase, which are comparable with those usually observed during the storm main phase. The IRI STORM model gave no reliable corrections of foF2 for analyzed events.     

Response of the ionosphere to super geomagnetic storms: Observations and modeling

The relative importance of the main drivers of positive ionospheric storms at low-mid latitudes is studied using observations and modeling for the first time. In response to a rare super double geomagnetic storm during 07–11 November 2004, the low-mid latitude (17°–48°N geomag. lat.) ionosphere produced positive ionospheric storms in peak electron density (NmF2) in Japan longitudes (≈125°–145°E) on the day of main phase (MP1) onset (06:30 LT) and negative ionospheric storms in American longitudes (≈65°–120°W) on the following day of MP1 onset (13:00–16:00 LT). The relative effects of the main drivers of the positive ionospheric storms (penetrating daytime eastward electric field, and direct and indirect effects of equatorward neutral wind) are studied using the Sheffield University Plasmasphere Ionosphere Model (SUPIM). The model results show that the penetrating daytime (morning–noon) eastward electric field shifts the equatorial ionisation anomaly crests in NmF2 and TEC (total electron content) to higher than normal latitudes and reduces their values at latitudes at and within the anomaly crests while the direct effects of the equatorward wind (that reduce poleward plasma flow and raise the ionosphere to high altitudes of reduced chemical loss) combined with daytime production of ionisation increase NmF2 and TEC at latitudes poleward of the equatorial region; the later effects can be major causes of positive ionospheric storms at mid latitudes. The downwelling (indirect) effect of the wind increases NmF2 and TEC at low latitudes while its upwelling (indirect) effect reduces NmF2 and TEC at mid latitudes. The net effect of all main drivers is positive ionospheric storms at low-mid latitudes in Japan longitude, which qualitatively agrees with the observations.     

Study the behavior of the ionospheric frequencies at Karachi during three consecutive solar minima of cycle 21, 22 & 23 and their comparison with IRI-2016

The behavior of critical frequencies of ionospheric E and F2 layers (foE & foF2) along with minimum ionospheric frequency (fmin) is studied for solar minima of cycle 21 (1986), 22 (1996) and 23 (2008) over Karachi (24.95$°$N, 67.13$°$E), Pakistan. The station is located at the crest of equatorial ionization anomaly region. Beside seasonal differences, pronounced change in the values of frequencies is noted from one solar minimum to another solar minimum. A strong and direct correlation of foF2 with Smoothed Sunspot Number (SSN) and F10.7?cm solar flux is observed. In the minimum of cycle 23, reduction in foF2 is noted due to reduction of solar EUV as compared to other minima. Also disappearance of semi-annual variations in foF2 is noted in cycle 23 minimum. Unexpectedly higher values of foE and fmin are observed in minimum of cycle 23 as compared to other minima. It is difficult to explain this unusual behavior of fmin and foE along with disappearance of semi-annual variation in foF2. It is possible that during very low solar activity, thermospheric conditions are changed which in turn altered the ionosphere. Further investigation of atmosphere-ionosphere coupling is required to understand this complex behavior. On comparison of observed values with IRI-2016, higher deviations are observed in foE before noon hours while in case of foF2, large deviations are noted during daytime. The absence of foF2 semi-annual variation in cycle 23 is not reproduced by IRI-2016. It is suggested that IRI-2016 need some modification for extremely low solar activity condition.     

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.