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.     

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.     

Comparison of model and experimental ionospheric parameters at high latitudes

Modern use and study of the auroral region needs to attract a wider class of models for describing conditions of radio wave propagation in the ionosphere. In this paper the possibilities of the International Reference Ionosphere model, well-proven and widespread in the mid-latitudes, are investigated in the high latitude zone. Model and measured values of the critical frequency foF2 for two mid-latitude stations (Juliusruh and Goosebay) and four high-latitude ones (Loparskaya, Sodankyla, Sondrestrom, Thule) are compared. Deviations of medians, variations from day to day and solar activity trends seemed to be similar for both areas. This similarity is irrespective of the RZ12 index. Special attention is paid to the TEC parameter and its determination using 6 versions of models, a new version of the model IRI2010 (IRI-Plas) among them. It is shown that the IRI-Plas model significantly improves the definition of TEC in contrast to the versions of IRI2007 and the new model NeQuick. The use of the median of the experimental equivalent slab thickness, together with the current values of the TEC, increases by a factor of two the agreement between calculated and measured foF2 values as compared with the variations from day to day. This allows foF2 to be defined in near-real time.     

Midday bottomside electron density profiles during moderate solar activity and comparison with IRI-2001

Bottomside electron density (Ne-h) profiles during midday (10–14 h) are analyzed using modern digital ionosonde observations at a low-middle latitude station, New Delhi (28.6N, 77.2E, dip 42.4N), for the period from January 2003 to December 2003, pertaining to moderate solar activity (MSA). Each individual profile is normalized with respect to the peak height and density (hmF2, NmF2) of the F2-region. These profiles are compared with those obtained from the International Reference Ionosphere (IRI-2001) model. Bilitza [Bilitza, D. International Reference Ionosphere 2000. Radio Sci. 36 (2), 261–275, 2001] using both the options namely: Gulyaeva’s model [Gulyaeva, T.L. Progress in ionospheric informatics based on electron density profile analysis of ionograms. Adv. Space Res. 7 (6) 39–48, 1987] and B0 Tab. option [Bilitza, D., Radicella, S.M., Reinisch, B.W., Adeniyi, J.O., Mosert Gonzalez, M.E., Zhang, S.R., Obrou, O. New B0 and B1 models for IRI. Adv. Space Res. 25 (1), 89–95, 2000]. The study reveals that during summer and equinox, the IRI model with B0 Tab. option in general, produces better agreement with the observed median profiles, while the IRI predictions using Gulyaeva’s option, overestimate the electron density distribution at all the heights below the F2-peak. However, during winter, in general, the IRI model, using both the options, reveals shows fairly good agreement with the observations.     

Dynamics of global electron content in 1998–2005 derived from global GPS data and IRI modeling

We analyzed the dynamics of global electron content (GEC) for the period 1998–2005 and compared the estimated GEC with variations of the 10.7-cm solar radio emission and with and with GEC values obtained with IRI-2001. We found a strong resemblance between the curves’ shapes for the experimental and modeled GEC: strong semiannual variations are discernible in these series and both curves tend to increase the absolute GEC value during the period of maximum of solar activity. However, there are some significant distinctions, such as absence of 27-day fluctuations in the series of GEC computed by the IRI-2001. On the contrary, observational GEC reflects well dynamics of solar activity: 27-day variations of GEC are very similar to the ones of the index F10.7, but GEC undergoes a lagging of about of 30–60 h as compared to value of the F10.7 index. The relative amplitude of 27-day variations decreases from 8% at the rising and falling solar activity to 2% at the period of its maximum.     

Low solar activity variability and IRI 2007 predictability of equatorial Africa GPS TEC

Diurnal, seasonal and latitudinal variations of Vertical Total Electron Content (VTEC) over the equatorial region of the African continent and a comparison with IRI-2007 derived TEC (IRI-TEC), using all three options (namely; NeQuick, IRI01-corr and IRI-2001), are presented in this paper. The variability and comparison are presented for 2009, a year of low solar activity, using data from thirteen Global Positioning System (GPS) receivers. VTEC values were grouped into four seasons namely March Equinox (February, March, April), June Solstice (May, June, July), September Equinox (August, September, October), and December Solstice (November, December, January). VTEC generally increases from 06h00 LT and reaches its maximum value at approximately 15h00–17h00 LT during all seasons and at all locations. The NeQuick and IRI01-corr options of the IRI model predict reasonably well the observed diurnal and seasonal variation patterns of VTEC values. However, the IRI-2001 option gave a relatively poor prediction when compared with the other options. The post-midnight and post-sunset deviations between modeled and observed VTEC could arise because NmF2 or the shape of the electron density profile, or both, are not well predicted by the model; hence some improvements are still required in order to obtain improved predictions of TEC over the equatorial region of the Africa sector.     

A morphological study of GPS-TEC data at Agra and their comparison with the IRI model

The temporal and seasonal variations of Total Electron Content (TEC) are studied at Agra (Geographic Lat. 27.17°N, Long. 78.89°E, Dip: 41.4°), India, which is in the equatorial anomaly region, for a period of 12 months from 01 January to 31 December, 2007 using a Global Positioning System (GPS) receiver. The mean TEC values show a minimum at 0500 h LT (LT = UT + 5.5 h) and a peak value at about 1400 h LT. The lowest TEC values are observed in winter whereas largest values are observed in equinox and summer. Anomalous variations are found during the period of magnetic disturbances. These results are compared with the TEC derived from IRI-2007 using three different options of topside electron density, NeQuick, IRI01-corr, and IRI-2001. A good agreement is found between the TEC obtained at Agra and those derived from IRI models.     

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

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

Modelling ionospheric vertical drifts over Africa low latitudes using Empirical Orthogonal functions and comparison with climatological model

For the first time, empirical model of daytime vertical $\mathbf{E}\times \mathbf{B}$ drift based on Empirical Orthogonal functions (EOF) decomposition technique is presented. Day-to-day variability of $\mathbf{E}\times \mathbf{B}$ drift inferred from horizontal (H) geomagnetic field data around dip latitude for the period of 2008–2013 is used to both develop and validate the model. Results show that the EOF technique is promising with modelled values and data giving correlation coefficient values of at least 0.90 for geomagnetic conditions of both $K_p?3$ and $K_p>3$ within 2008–2013. Independent model validation shows that in situ $\mathbf{E}\times \mathbf{B}$ values from ion velocity meter (IVM) instrument on-board C/NOFS satellite are closer to model $\mathbf{E}\times \mathbf{B}$ estimates than the climatological Scherliess-Fejer (SF) model incorporated within the International Reference Ionosphere (IRI).     

磁暴期间中低纬电离层暴形态和物理机制分析

本文利用东亚地区１２个低纬电离层台站的测高仪观测数据，对１９７８年８月２７日发生的一次曲型磁暴期间电离层峰值高度和密度的变化进行了分析。采用滑动平均区分开电离层中不同时间尺度的扰动，分析了影响中低纬度电离层暴的几种扰动形态特征，并对其物理机制进行了讨论。结果表明：伴随磁暴急始的磁层压缩，电离层中表现出峰值密度增加和峰值高度下降；磁暴主相期间热层大气暴环流及其所引起的中性大气成分变化控制着电离层的大     

A study on Nighttime Winter Anomaly (NWA) and other related Mid-latitude Summer Nighttime Anomaly (MSNA) in the light of International Reference Ionosphere (IRI) – Model

The ionospheric Nighttime Winter Anomaly (NWA) is a feature observed in the Northern Hemisphere at the American and in the Southern Hemisphere at the Asian longitude sector under low solar activity conditions. Jakowski et al. (2015) analyzed ground-based GPS derived TEC and peak electron density data from radio occultation measurements on Formosat-3/COSMIC satellites and confirmed the persistence of the phenomenon. Further, they assumed that Mid-latitude Summer Nighttime Anomaly (MSNA) and related special anomalies such as the Weddell Sea Anomaly (WSA) and the Okhotsk Sea Anomaly (OSA) are closely related to the NWA via enhanced wind-induced uplifting of the ionosphere. The aim of this paper is to study the factors causing these anomalies and also to investigate if these anomalies are re-produced by IRI. The results show that IRI model does include the NWA effect, though at a different longitude and could be improved for better predictions. The IRI-2016 model does show WSA in TEC but not in NmF2. Further, the IRI-2016 model could clearly predict the OSA both in NmF2 and TEC.     

Characteristics of nighttime E-region over Arecibo: Dependence on solar flux and geomagnetic variations

Electron concentration (Ne) inferred from Incoherent Scatter Radar (ISR) measurements has been used to determine the influence of solar flux and geomagnetic activity in the ionospheric E-region over Arecibo Observatory (AO). The approach is based on the determination of column integrated Ne, referred to as E-region total electron content (ErTEC) between 80 and 150 km altitude regions. The results discussed in this work are for the AO nighttime period. The study reveals higher ErTEC values during the low solar flux periods for all the seasons except for summer period. It is found that the E-region column abundance is higher in equinox periods than in the winter for low solar activity conditions. The column integrated Ne during the post-sunset/pre-sunrise periods always exceeds the midnight minima, independent of season or solar activity. This behavior has been attributed to the variations in the coupling processes from the F-region. The response of ErTEC to the geomagnetic variability is also examined for different solar flux conditions and seasons. During high solar flux periods, changes in Kp cause an ErTEC increase in summer and equinox, while producing a negative storm-like effect during the winter. Variations in ErTEC due to geomagnetic activity during low solar flux periods produce maximum variability in the E-region during equinox periods, while resulting in an increase/decrease in ErTEC before local midnight during the winter/summer periods, respectively.     

Seasonal,annual, and interannual variability in MLT quasi-two-day waves over the low-latitude region Kolhapur (16.8°N; 74.2°E)

This study presents the quasi-two-day wave (Q2DW) characteristics of the mesosphere and lower thermosphere (MLT) region obtained by taking hourly mean values of horizontal wind velocities for 4? years (August 2013–July 2017) through continuous measurements using a medium-frequency (MF) radar (operating frequency – 1.98?MHz) located at the low-latitude Indian station Kolhapur (16.8°N; 74.2°E). The MF radar located at Kolhapur was upgraded in 2013, and these results of Q2DW have been reported for the first time after upgrading. The present study investigated variability in seasonal, annual, interannual, and solar indices of Q2DWs traveling in zonal (EW) and meridional (NS) components in the MLT region. The Q2DW activity is observed to be stronger during austral summer (January–February) (EW?=?～5?m/s and NS?=?～8–10?m/s) than during boreal summer (June–July) (EW = ～5 m/s and NS = ～6–8?m/s). The Q2DW amplitudes are larger in the meridional component than in the zonal one. A strong semiannual oscillation (SAO) has been observed in Q2DWs, with peak during January–February and June–July. In addition, small enhancement is seen in meridional Q2DW in October (～5–6?m/s). It is observed that the entire spectrum (40–60?h) measured between 86 and 94?km contributes to the SAO amplitudes during January–February and June–July, whereas the waves measured between 42?h and 52?h contribute to enhancement in October similar to that reported elsewhere. In general, the Q2DW amplitude shows large interannual variability. The easterlies developed in the global circulation model in Northern hemisphere during May intensify up to around summer solstice. Q2DW activity peaks during westerly shear zone and intensifies with time at a lower thermospheric altitude (above 90?km). Small positive correlations (r?=?0.2 for sunspot number and r?=?0.1 for 10.7?cm solar flux) have been observed between Q2DW amplitudes and solar activity.     

Nighttime ionosphere–thermosphere coupling observed during an intense geomagnetic storm

The electrodynamics of the ionosphere in the tropical region presents various scientific aspects, which remain subject of intensive investigations and debates by the scientific community. During the year 2002, in a joint project between the Universidade do Vale do Paraíba (UNIVAP) and Universidade Luterana do Brasil (ULBRA), a chain of three Canadian Advanced Digital Ionosondes (CADIs) was established nearly along the geomagnetic meridian direction, for tropical ionospheric studies, such as, changes and response due to geomagnetic disturbances and thermosphere–ionosphere coupling and the generation and dynamics of ionospheric irregularities, in the Brazilian sector. The locations of the three ionosondes stations are São José dos Campos (23.2°S, 45.9°W, dip latitude 17.6°S – under the southern crest of equatorial ionospheric anomaly), Palmas (10.2°S, 48.2°W, dip latitude 5.5°S – near the magnetic equator) and Manaus (2.9°S, 60.0°W, dip latitude 6.4°N – between the geographic and geomagnetic dip equators). It should be pointed out that Palmas and Manaus are located on the opposite sides of the magnetic equator but both are south of the geographic equator. The three CADIs work in time-synchronized mode and obtain ionograms every 5 min. This configuration of the ionospheric sounding stations allowed us to study the F-region dynamics during geomagnetically disturbed period in the meridional direction. Just after the installation and testing of the three CADIs, on September 05, 2002 a coronal mass ejection (CME) left the Sun and about 2 days after the CME left the Sun, it reached the Earth’s magnetosphere and complex and multi step events took place during the period September 07–09. In the study we note that the equatorial stations located north (Manaus, dip latitude 6.4°N) and south (Palmas, dip latitude 5.5°S) of the dip equator presented significant F-layer height asymmetries during the storm main phase. In addition, the low-latitude station SJC (dip latitude 17.6°S) presented decrease in the F-layer densities (negative phase), whereas Palmas presented increase in the F-layer densities (positive phase) during the main phase. This was followed by positive phase at both the stations. During the first night of the recovery phase a strong formation and evolution of large-scale ionospheric irregularities (equatorial spread-F (ESF)) was observed, but on the second night of the recovery phase, there was strong and almost simultaneous sporadic E (Es) formation at all three stations. During the presence of Es, spread-F formation is not observed, indicating the suppression of spread-F, possibly by sporadic E.     

TEC derived from some GPS stations in Nigeria and comparison with the IRI and NeQuick models

Total electron content (TEC) measured simultaneously using Global Positioning System (GPS) ionospheric monitors installed at some locations in Nigeria during the year 2011 (Rz = 55.7) was used to study the diurnal, seasonal, and annual TEC variations. The TEC exhibits daytime maximum, seasonal variation and semiannual variations. Measured TEC were compared with those predicted by the improved versions of the International Reference Ionosphere (IRI) and NeQuick models. The models followed the diurnal and seasonal variation patterns of the observed values of TEC. However, IRI model produced better estimates of TEC than NeQuick at all locations.