Comparison between ionospheric peak parameters retrieved from COSMIC measurement and ionosonde observation over Sanya

In this paper we compared the ionospheric peak parameters (peak electron density of the F region, NmF2, and peak height of the F region, hmF2) retrieved from the FORMASAT-3/COSMIC (COSMIC for short) satellite measurement with those from ionosonde observation at Sanya (18.3°N, 109.6°E) during the period of 2008–2013. Although COSMIC NmF2 (hmF2) tends to be lower (higher) than ionosonde NmF2 (hmF2), the results show that the ionospheric peak parameters retrieved from COSMIC measurement generally agree well with ionosonde observation. For NmF2 the correlation between the COSMIC measurement and the ionosonde observation is higher than 0.89, and for hmF2 the correlation is higher than 0.80. The correlation of the ionospheric peak parameters decreases when solar activity increases. The performance of COSMIC measurement is acceptable under geomagnetic disturbed condition. The correlation of NmF2 between COSMIC and ionosonde measurements is higher (lower) during the nighttime (daytime), while the correlation of hmF2 is lower (higher) during the nighttime (daytime).     

Ionospheric response to the storm-time disturbance of 29 May, 2010

The Ionospheric F2-layer peak parameters response to a magnetic storm had been investigated over Ilorin, Nigeria (Lat. 8:53°N, Long. 4.5°E, dip angle, −2.96°), Jicamarca, Peru (11.95°S, 76.87°W, dip angle, 0.8°) and Hermanus, South Africa (34.42°S, 19.22°E, dip angle, −60.77°), using percentage enhancement/depletion values. Our results showed an enhancement in NmF2 at all of these stations. Averagely, pre-noon and post-noon peaks are highest at Ilorin during quiet time. The similar pattern observed for quiet condition between Ilorin and Jicamarca was due to their latitudinal positions. For disturbed NmF2 condition, Jicamarca and Ilorin recorded higher peaks at nighttime than during the daytime for the storms main phase, and the reverse over Hermanus. The nighttime and daytime increases were observed respectively at Ilorin and Hermanus during the recovery period. The hmF2 variation recorded higher enhancement at Jicamarca during the daytime and at Hermanus at nighttime during the main phase. During the recovery phase, the highest enhancement was recorded during the daytime at Jicamarca, and over Hermanus at nighttime. These observations find their explanation in the magnetospheric current, solar wind and E × B drift.     

Multi-station observation of ionospheric irregularities over South Africa during strong geomagnetic storms

This paper presents results pertaining to the response of the mid-latitude ionosphere to strong geomagnetic storms that occurred from 31 March to 02 April 2001 and 07–09 September 2002. The results are based on (i) Global Positioning Systems (GPSs) derived total electron content (TEC) variations accompanying the storm, (ii) ionosonde measurements of the ionospheric electrodynamic response towards the storms and (iii) effect of storm induced travelling ionospheric disturbances (TIDs) on GPS derived TEC. Ionospheric data comprising of ionospheric TEC obtained from GPS measurements, ionograms, solar wind data obtained from Advanced Composition Explorer (ACE) and magnetic data from ground based magnetometers were used in this study. Storm induced features in vertical TEC (VTEC) have been obtained and compared with the mean VTEC of quiet days. The response of the mid-latitude ionosphere during the two storm periods examined may be characterised in terms of increased or decreased level of VTEC, wave-like structures in VTEC perturbation and sudden enhancement in hmF2 and h′F. The study reveals both positive and negative ionospheric storm effects on the ionosphere over South Africa during the two strong storm conditions. These ionospheric features have been mainly attributed to the travelling ionospheric disturbances (TIDs) as the driving mechanism for the irregularities causing the perturbations observed. TEC perturbations due to the irregularities encountered by the satellites were observed on satellites with pseudo random numbers (PRNs) 15, 17, 18 and 23 between 17:00 and 23:00 UT on 07 September 2002.     

Variation of hmF2 and NmF2 deduced from DPS-4 over Multan (Pakistan) and their comparisons with IRI-2012 & IRI-2016 during the deep solar minimum between cycles 23 & 24

We report the results of ionospheric measurements from DPS-4 installed at Multan (Geog coord. 30.18°N, 71.48°E, dip 47.4°). The variations in F2-layer maximum electron density NmF2 and its peak height hmF2 are studied during the deep solar minimum between cycles 23 & 24 i.e 2008–2009 with comparisons conducted with the International Reference Ionosphere (IRI) versions 2012 & 2016. We find that the hmF2 observations peak around the pre-sunrise and sunrise hours depending on the month. Seasonally, the daytime variation of NmF2 is higher in the Equinox and Summer, while daytime hmF2 are slightly higher in the Equinox and Winter. High values of hmF2 around midnight are caused by an increase of upward drifts produced by meridional winds. The ionosphere over Multan, which lies at the verge of low and mid latitude, is affected by both $E\times B$ drifts and thermospheric winds as evident from mid-night peaks and near-sunrise dips in hmF2. The results of the comparison of the observed NmF2 and hmF2 for the year 2008–2009 with the IRI-2012 (both NmF2 and hmF2) and IRI-2016 (only hmF2) estimates indicate that for NmF2, IRI-2012 with Consultative Committee International Radio (CCIR) option produces values in better agreement with observed data. Whereas, for hmF2, IRI-2016 with both International Union of Radio Science (URSI) and CCIR SHU-2015 options, predicts well for nighttime hours throughout the year. However, the IRI-2012 with CCIR option produces better agreement with data during daytime hours. Furthermore, IRI-2012 with CCIR option gives better results during Equinox months, whereas, IRI-2016 with both URSI and CCIR SHU-2015 options predict well for Winter and Summer.     

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.     

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.     

Evaluation of the improvement of IRI-2016 over IRI-2012 at the India low-latitude region during the ascending phase of cycle 24

This paper investigated the performance of the latest International Reference Ionosphere model (IRI-2016) over that of IRI-2012 in predicting total electron content (TEC) at three different stations in the Indian region. The data used were Global Positioning System (GPS) data collected during the ascending phase of solar cycle 24 over three low-latitude stations in India namely; Bangalore (13.02°N Geographic latitude, 77.57°E Geographic longitude), Hyderabad (17.25°N Geographic latitude, 78.30°E Geographic longitude) and Surat (21.16°N Geographic latitude, 72.78°E Geographic longitude). Monthly, the seasonal and annual variability of GPS-TEC have been compared with those derived from International Reference Ionosphere IRI-2016 and IRI-2012 with two different options of topside electron density: NeQuick and IRI01-corr. It is observed that both versions of IRI (i.e., IRI-2012 and IRI-2016) predict the GPS-TEC with some deviations, the latest version of IRI (IRI-2016) predicted the TEC similar to those predicted by IRI-2012 for all the seasons at all stations except for morning hours (0500 LT to 1000?LT). This shows that the effect of the updated version is seen only during morning hours and also that there is no change in TEC values by IRI-2016 from those predicted by IRI-2012 for the rest of the time of the day in the Indian low latitude region. The semiannual variations in the daytime maximum values of TEC are clearly observed from both GPS and model-derived TEC values with two peaks around March-April and September-October months of each year. Further, the Correlation of TEC derived by IRI-2016 and IRI-2012 with EUV and F10.7 shows similar results. This shows that the solar input to the IRI-2016 is similar to IRI 2012. There is no significant difference observed in TEC, bottom-side thickness (B0) and shape (B1) parameter predictions by both the versions of the IRI model. However, a clear improvement is visible in hmF2 and NmF2 predictions by IRI-2016 to that by IRI-2012. The SHU-2015 option of the IRI-2016 gives a better prediction of NmF2 for all the months at low latitude station Ahmedabad compare to AMTB 2013.     

The variability and IRI2007-predictability of hmF2 over South Africa

This paper presents an investigation into the variability and predictability of the maximum height of the ionospheric F2 layer, hmF2 over the South African region. Data from three South African stations, namely Madimbo (22.4°S, 26.5°E, dip angle: −61.47°), Grahamstown (33.3°S, 26.5°E, dip angle: −64.08°) and Louisvale (28.5°S, 21.2°E, dip angle: −65.44°) were used in this study. The results indicate that hmF2 shows a larger variability around midnight than during the daytime for all seasons. Monthly median hmF2 values were used in all cases and were compared with predictions from the IRI-2007 model, using the URSI (Union Radio-Scientifique Internationale) coefficient option. The analysis covers the diurnal and seasonal hourly hmF2 values for the selected months and time sectors e.g. January, July, April and October for 2003 and 2005. The time ranges between (03h00–23h00 UT; LT = UT + 2h) representing the local sunrise, midday, sunset and midnight hours. The time covers sunrise, midday, sunrise, and midnight hours (03–06h00 UT, 07–11h00 UT, sunrise 16–18h00 UT and 22–23h00 UT; LT = UT + 2h). The dependence of the results on solar activity levels was also investigated. The IRI-2007 predictions follow fairly well the diurnal and seasonal variation patterns of the observed hmF2 values at all the stations. However, the IRI-2007 model overestimates and underestimates the hmF2 value during different months for all the solar activity periods.     

Time delay of ionospheric TEC storms to geomagnetic storms and pre-storm disturbance events in East Asia

Using the TEC data at Beijing (39.61°N, 115.89°E)/Yakutsk (62.03°N, 129.68°E) stations of East Asia regions and relevant geomagnetic data from 2010 to 2017, we have studied the time delay of ionospheric storms to geomagnetic storms and compare it with our previous results of Taoyuan (25.02°N, 121.21°E) station (Zhang et al., 2020). The data shows a well-known local time dependence of the time delay, and seasonal dependences are different at these stations. In addition, there is no correlation between the time delay and the magnetic storm intensity /solar activity, except the time delay of negative storms has weakly negative dependence on the solar activity. Comparing with the results of Taoyuan station which is located at EIA region in East Asia, we find that the time delay increases nonlinearly as the latitude decreases due to different ionospheric backgrounds at these places. Moreover, the pre-storm disturbance events are found to have similar statistical characteristics as the pre-storm enhancement in Europe middle latitudes (Bure?ová and La?tovi?ka, 2007). By subtracting the common features of the pre-storm disturbance events, we preliminarily infer that auroral activity might be main driver of the pre-storm disturbance events.     

A 2-D comparison of ionospheric convection derived from SuperDARN and DMSP measurements

Ion drift vectors measured by the DMSP satellites are compared with plasma convection vectors obtained by the SuperDARN HF radars through the standard Map Potential algorithm of Ruohoniemi and Baker [Ruohoniemi, M., Baker, K.B. Large-scale imaging of high-latitude convection with super dual auroral radar network HF radar observations. J. Geophys. Res. 103, 20797–20811, 1998]. Despite significant data spread, the agreement can be qualified as reasonable for a data set comprising of 149 satellite passes over the Northern Hemisphere at high latitudes. The slope of the best-fit line relating SuperDARN and DMSP velocity magnitudes is of the order of 0.3 with a tendency for the SuperDARN velocities to be smaller. The agreement between the azimuths of the ion drift and convection is better with the slope of the best-fit line being close to 1. It is shown that consistency between the radar and satellite measurements is much better if the SuperDARN line-of-sight velocities are compared with the DMSP cross-track ion drifts for events showing slow spatial and temporal variations of the convection. If areas of strong convection changes are included into comparison, the degree of agreement deteriorates drastically. This result implies that differences in the spatial and temporal resolutions of DMSP and SuperDARN measurements are crucial factors contributing to the observed discrepancies. In addition, some differences are introduced when the SuperDARN line-of-sight velocities are filtered and reprocessed into vectors with the application of a background convection model.     

A global survey of COSMIC ionospheric peak electron density and its height: A comparison with ground-based ionosonde measurements

With a network of ground-based ionosondes distributed around the world, the ionospheric peak electron density and its height measured by FORMOSAT-3/COSMIC satellites in terms of GPS radio occultation technique are extensively examined in this article. It is found that, in spite of the latitude, the mean values of the peak electron density measured by COSMIC satellites are systematically smaller than those observed by ground-based ionosondes. The discrepancy between them is dependent on the latitude, namely, it is small in low and mid-latitudes and large in high-latitude region. Moreover, statistical analysis shows that the slopes of the regression line that is best fitted to the scatter diagram of occultation-retrieved peak electron density (ordinate axis) versus ionosonde-observed peak density (abscissa axis) are universally less than one. This feature is believed to be the result of path average effect of non-uniform distribution of the electron density along the GSP ray during the occultation. A comparison between COSMIC-measured peak height and ionosonde-derived peak height hmF2 indicates that the former is systematically higher than the latter. The difference in the two can be as large as 20% or more in equatorial and low-latitude regions. This result implies that the peak height hmF2 derived from the virtual height through true height analysis based on Titheridge method seems to underestimate the true peak height. The correlation between COSMIC and ionosonde peak electron densities is analyzed and the result reveals that correlation coefficient seems to be dependent on the fluctuation of the occultation-retrieved electron density profile. The correlation will be higher (lower) for the electron density profiles with smaller (larger) fluctuations. This feature suggests that the inhomogeneous distribution of the electron density along the GPS ray path during the occultation plays an important role affecting the correlation between COSMIC and ionosonde measurements.     

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.     

Bottom side profiles for two close stations at the southern crest of the EIA: Differences and comparison with IRI-2012 and NeQuick2 for low and high solar activity

Bottom side electron density profiles for two stations at the southern crest of the Equatorial Ionization Anomaly (EIA), São José dos Campos (23.1°S, 314.5°E, dip latitude 19.8°S; Brazil) and Tucumán (26.9°S, 294.6°E, dip latitude 14.0°S; Argentina), located at similar latitude and separated by only 20° in longitude, have been compared during equinoctial, winter and summer months under low (year 2008, minimum of the solar cycle 23/24) and high solar activity (years 2013–2014, maximum of the solar cycle 24) conditions. An analysis of parameters describing the bottom side part of the electron density profile, namely the peak electron density NmF2, the height hmF2 at which it is reached, the thickness parameter B₀ and the shape parameter B₁, is carried out. Further, a comparison of bottom side profiles and F-layer parameters with the corresponding outputs of IRI-2012 and NeQuick2 models is also reported. The variations of NmF2 at both stations reveal the absence of semi-annual anomaly for low solar activity (LSA), evidencing the anomalous activity of the last solar minimum, while those related to hmF2 show an uplift of the ionosphere for high solar activity (HSA). As expected, the EIA is particularly visible at both stations during equinox for HSA, when its strength is at maximum in the South American sector. Despite the similar latitude of the two stations upon the southern crest of the EIA, the anomaly effect is more pronounced at Tucumán than at São José dos Campos. The differences encountered between these very close stations suggest that in this sector relevant longitudinal-dependent variations could occur, with the longitudinal gradient of the Equatorial Electrojet that plays a key role to explain such differences together with the 5.8° separation in dip latitude between the two ionosondes. Furthermore at Tucumán, the daily peak value of NmF2 around 21:00 LT during equinox for HSA is in temporal coincidence with an impulsive enhancement of hmF2, showing a kind of “elastic rebound” under the action of the EIA. IRI-2012 and NeQuick2 bottom side profiles show significant deviations from ionosonde observations. In particular, both models provide a clear underestimation of the EIA strength at both stations, with more pronounced differences for Tucumán. Large discrepancies are obtained for the parameter hmF2 for HSA during daytime at São José dos Campos, where clear underestimations made by both models are observed. The shape parameter B₀ is quite well described by the IRI-2012 model, with very good agreement in particular during equinox for both stations for both LSA and HSA. On the contrary, the two models show poor agreements with ionosonde data concerning the shape parameter B₁.     

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.     

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.