On the uncertainties in the measurement of absolute (true) TEC over Indian equatorial and low latitude sectors

The Indian sector encompasses the equatorial and low latitude regions where the ionosphere is highly dynamic and is characterized by the equatorial ionization anomaly (EIA) resulting in large latitudinal electron density gradients causing errors and uncertainties in the estimation of range delays in satellite based navigation systems. The diurnal and seasonal variations of standard deviations in the TEC data measured during the low sunspot period 2004–2005 at 10 different Indian stations located from equator to the anomaly crest region and beyond are examined and presented. The day-to-day variability in TEC is found to be lowest at the equatorial station and increases with latitude up to the crest region of EIA and decreases beyond.     

GPS scintillation and TEC depletion near the northern crest of equatorial anomaly over South China

This study presents a statistical analysis of GPS L-band scintillation with data observed from July 2008 to March 2012 at the northern crest of equatorial anomaly stations in Guangzhou and Shenzhen of South China. The variations of the scintillation with local time, season, solar activity and duration of scintillation patches were investigated. The relationship between the scintillation and TEC depletion was also reported. Our results revealed that GPS scintillation occurred from 19:30 LT (pre-midnight) to 03:00 LT (post-midnight). During quiet solar activity years, the scintillation was only observed in pre-midnight hours of equinox months and patches durations were mostly less than 60 min. During high solar activity years, more scintillation occurred in the pre-midnight hours of equinox and winter months; and GPS scintillation started to occur in the post-midnight hours of summer and winter. The duration of scintillation patches extended to 180 min in high solar activity years. Solar activity had a larger effect to strong scintillations (S₄ > 0.6) than to weak scintillations (0.6 ? S₄ > 0.2). Strong scintillations were accompanied by TEC depletion especially in equinox months. We also discussed the relationship between TEC depletion and plasma bubble.     

Vertical TEC variations and model during low solar activity at a low latitude station,Xiamen

GPS-derived vertical TEC recorded at Xiamen (24.5°N, 118.1°E, geomagnetic latitude 13.2°N), China, during year 2006 is analyzed for the first time and compared to that predicted by ionosphere model SPIM recommend by ISO. A manifest seasonal anomaly is found with the high value during equinoctial season and low value during summer and winter season. Relative standard deviation for VTEC shows high value at around midnight and before sunrise. The correlation analysis exhibits that the variation of VTEC has a very weak relation with geomagnetic and solar activities (Dst, AP, SSN and F10.7). Comparative results reveal that the SPIM overestimates the observed VTEC at most of the time.     

Numerical simulation of equatorial plasma bubbles over Cachimbo: COPEX campaign

The problem of day-to-day variability in onset of equatorial spread F (ESF) is addressed using data from the 2002 COPEX observational campaign in Brazil and numerical modeling. The observational results show that for values of virtual height of the F layer base less than 355 km at around 18:35 LT, and for the prereversal peak enhancement of the vertical plasma drift (Vp) less than 30 m/s, the spread-F (ESF) was absent on four nights over Cachimbo (9.5°S, 54.8°W, dip latitude = −2.1°). In this work we analyze the geophysical conditions for the generation of the irregularities by comparing the nights with and without the ESF. In the comparison a numerical code is used to simulate plasma irregularity development in an extended altitude range from the bottom of the equatorial F   layer. The code uses the flux corrected transport method with Boris–Book’s flux limiter for the spatial integration and a predictor–corrector method for the direct time integration of the continuity equation for O⁺${\mathrm{O}}^{+}$ and the SOR (Successive-Over-Relaxation) method for electric potential equation. The code is tested with different evening eastward electric fields (or vertical drifts Vp < 30 m/s and Vp > 30 m/s) in order to study the influence of the prereversal enhancement in the zonal electric field on plasma bubble formation and development. The code also takes into account the zonal wind, the vertical electric field and the collision frequency of ions with neutrals and the amplitude of initial perturbation. The simulation shows a good agreement with the observational results of the ESF. The results of the code suggest that the instability can grow at the F layer bottomside by the Rayleigh–Taylor mechanism only when the Vp > 30 m/s. In the analyzed cases we have considered the competition of other geophysical parameters in the generation of plasma structures.     

Vertical TEC in the equatorial anomaly region from satellite altimetry data

The high-precision demands imposed by the ocean altimetry community of the late 1980 resulted in the TOPEX/Poseidon mission. This mission was the first to carry as its main instrument a dual-frequency sea-altimeter on board a satellite. This instrument together with other state-of-the-art technologies involved in the mission, led to sea-height determinations with precision better than 2 cm. As a by-product, the TOPEX/Poseidon mission provided vertical TEC determinations that since they became available, have demonstrated to be a powerful tool for ionospheric studies.     

SAMI2 model results for the quiet time low latitude ionosphere over India

Efficacy of SAMI2 model for the Indian low latitude region around 75°E longitudes has been tested for different levels of solar flux. With a slight modification of the plasma drift velocity the SAMI2 model has been successful in reproducing quiet time ionospheric low latitude features like Equatorial Ionization Anomaly. We have also showed the formation of electron hole in the topside equatorial ionosphere in the Indian sector. Simulation results show the formation of electron hole in the altitude range 800–2500?km over the magnetic equator. Indian zone results reveal marked differences with regard to the time of occurrence, seasonal appearances and strength of the electron hole vis-a-vis those reported for the American equatorial region.     

Variability study of the crest-to-trough TEC ratio of the equatorial ionization anomaly around 120°E longitude

In this paper, latitudinal profiles of the vertical total electron content (TEC) deduced from the dual-frequency GPS measurements obtained at ground stations around 120°E longitude were used to study the variability of the equatorial ionization anomaly (EIA). The present study mainly focuses on the analysis of the crest-to-trough TEC ratio (TEC-CTR) which is an important parameter representing the strength of EIA. Data used for the present study covered the time period from 01 January, 1998 to 31 December, 2004. An empirical orthogonal function analysis method is used to obtain the main features of the TEC-CTR’s diurnal and seasonal variations as well as its solar activity level dependency. Our results showed that: (1) The diurnal variation pattern of the TEC-CTR at 120°E longitude is characterized by two remarkable peaks, one occurring in the post-noon hours around 13–14 LT, and the other occurring in the post-sunset hours around 20–21 LT, and the post-sunset peak has a much higher value than the post-noon one. (2) Both for the north and south crests, the TEC-CTR at 120°E longitude showed a semi-annual variation with maximum peak values occurring in the equinoctial months. (3) TEC-CTR for the north crest has lower values in summer than in winter, whereas TEC-CTR for the south crest does not show this ‘winter anomaly’ effect. In other words, TEC-CTR for both the north and south crests has higher values in the northern hemispheric winter than in the northern hemispheric summer. (4) TEC-CTR in the daytime post-noon hours (12–14 LT) does not vary much with the solar activity, however, TEC-CTR in the post-sunset hours (19–21 LT) shows a clear dependence on the solar activity, its values increasing with solar activity.     

Investigation of ionospheric irregularities and scintillation using TEC at high latitude

The occurrence of ionospheric irregularities at high latitudes, with dimensions of several kms down to decameter scale size shows strong correlation with geomagnetic disturbance, season and solar activity. Transionospheric radio waves propagating through these irregularities experience rapid random fluctuations in phase and/or amplitude of the signal at the receiver, termed scintillation, which can degrade GNSS services. Thus, investigation and prediction of this scintillation effect is very important. To investigate such scintillation effects, a GISTM (GPS Ionospheric Scintillation and TEC Monitoring) NovAtel dual frequency (L1/L2) GPS receiver has been installed at Trondheim, Norway (63.41°$63.41°$ N, 10.4°$10.4°$ E), capable of collecting scintillation indices at a 1 min rate as well as the raw data (phase and intensity) of the satellite signals at a 50 Hz sampling rate and TEC (Total Electron Content) at a 1 Hz rate. Many researchers have reported that both phase and amplitude scintillation is closely associated with TEC fluctuations or associated with a significant developing enhancement or depletion in the TEC. In this study, a novel analogous phase index is developed which provides samples at a 1 min rate. Generally the scintillation indices can help in estimating the irregularity scintillation effect at a one minute rate, but such procedures are time consuming if DFTs of the phase and/or amplitude at a 50 Hz data are required. In this study, instead, this analogous phase index is estimated from 1 Hz rate TEC values obtained from the raw signals and is then compared for weak, moderate and strong scintillation at Trondheim for one year of data collected from the installed GPS receiver. The spectral index of the irregularities (that is the inverse power law of their spatial spectrum) is determined from the resultant phase scintillation psd. The correlations of the scintillation indices and spectral indices with the analogous phase index have been investigated under different geomagnetic conditions (represented by the Kp index) and an approximate linear correlation of phase scintillation with the analogous phase index was found. Then a principal advantage of this index is that it achieves this correlation without requiring a high sampling data rate and the need for DFTs. Thus, the index seems a good candidate for developing a simple means of ionospheric scintillation prediction which could also be utilized in the development of alerts using regional mappings.     

Development of a campaign to study equatorial ionospheric phenomena over Guam

The United States Air Force Academy (USAFA) is in the process of developing a series of ground-based and space-based experiments to investigate the equatorial ionosphere over Guam and the southern crest of the Equatorial Appleton Anomaly over New Guinea. On the ground the Digital Ionospheric Sounder (University of Massachusetts, Lowell DPS-4 unit) and a dual-frequency GPS TEC/scintillation monitor will be used to investigate ionospheric phenomena in both campaign and long-term survey modes. In campaign mode, we will combine these observations with those collected from space during USAFA’s FalconSAT-3 and FalconSAT-5 low Earth orbit satellite missions, which will be active over a period of several years beginning in the first quarter of the 2007 calendar year. Additionally, we will investigate the long-term morphology of key ionospheric characteristics useful for driving the International Reference Ionosphere, such as critical frequencies (f_oE, f_oF1, f_oF2, etc.), the M(3000) F2 parameter (the maximum useable frequency for a signal refracted within the F2 layer and received on the ground at a distance of 3000 km away), and a variety of other characteristics. Specific targets of investigation include: (a) a comparison of TEC observed by the GPS receiver with those calculated by IRI driven by DPS-4 observations, (b) a comparison of plasma turbulence observed on-orbit with ionospheric conditions as measured from the ground, and (c) a comparison between topside ionospheric satellite in situ measurements of plasma density during an overpass of a Digisonde versus the calculated value based on extrapolation of the electron density profiles using Digisonde data and a topside α-Chapman function. This last area of investigation is discussed in detail in this paper.     

Signatures of strong geomagnetic storms in the equatorial latitude

Ionosonde data from two equatorial stations in the African sector have been used to study the signatures of four strong geomagnetic storms on the height – electron density profiles of the equatorial ionosphere with the objective of investigating the effects and extent of the effects on the three layers of the equatorial ionosphere. The results showed that strong geomagnetic storms produced effects of varying degrees on the three layers of the ionosphere. Effect of strong geomagnetic storms on the lower layers of the equatorial ionosphere can be significant when compared with effect at the F2-layer. Fluctuations in the height of ionization within the E-layer were as much as 0% to +20.7% compared to −12.5% to +8.3% for the F2-layer. The 2007 version of the International Reference Ionosphere, IRI-07 storm-time model reproduced responses at the E-layer but overestimated the observed storm profiles for the F1- and F2-layers.     

Climatology of ionospheric scintillations and TEC trend over the Ugandan region

This study presents results on the investigation of the diurnal, monthly and seasonal variability of Total Electron Content (TEC), phase (_σΦ${\sigma }_{\mathrm{\Phi }}$) and amplitude (S4) scintillation indices over Ugandan (Low latitude) region. Scintillation Network Decision Aid (SCINDA) data was obtained from Makerere (0.34°N, 32.57°E) station, Uganda for two years (2011 and 2012). Data from two dual frequency GPS receivers at Mbarara (0.60°S, 30.74°E) and Entebbe (0.04°N, 32.44°E) was used to study TEC climatology during the same period of scintillation study. The results show that peak TEC values were recorded during the months of October–November, and the lowest values during the months of July–August. The diurnal peak of TEC occurs between 10:00 and 14:00 UT hours. Seasonally, the ascending and descending phases of TEC were observed during the equinoxes (March and September) and solstice (June and December), respectively. The scintillations observed during the study were classified as weak (0.1≤$\le$S4,_σΦ≤${\sigma }_{\mathrm{\Phi }}\le$0.3) and strong (0.3<$<$S4,_σΦ≤${\sigma }_{\mathrm{\Phi }}\le$1.0). The diurnal scintillation pattern showed peaks between 17:00 and 22:00 UT hour, while the seasonal pattern follows the TEC pattern mentioned above. Amplitude scintillation was more dominant than phase scintillation during the two years of the study. Scintillation peaks occur during the months of March–April and September–October, while the least scintillations occur during the months of June–July. Therefore, the contribution of this study is filling the gap in the current documentation of amplitude scintillation without phase scintillation over the Ugandan region. The scintillations observed have been attributed to wave-like structures which have periods of about 2–3 h, in the range of that of large scale travelling ionospheric disturbances (LSTIDs).     

Spatial distribution of TEC across India in 2005: Seasonal asymmetries and IRI prediction

Total electron content measured simultaneously at 10 locations over India during the low solar activity year 2005 is used to examine the temporal and spatial asymmetries and also to assess the predictability of the International Reference Ionosphere in respect of the observed asymmetrical distribution. The stations are distributed in latitude along 77°E and in longitude along 23°N forming a meridional and a zonal chain respectively. A longitudinal gradient positive towards east was observed in the daytime hours of equinox and summer. Equinoctial asymmetry was prevalent across India during this year. Within the crest and equator, winter anomaly has been observed. It is found that IRI 2012 (with Ne Quick option, URSI coefficients) is unable to fully capture the temporal variation and spatial gradients of the ionization density in the Indian sector during 2005. The amount of offset between the model and measurement varies with local time and location.     

Investigation of the characteristics of wavelike oscillations of post-sunset equatorial ionospheric irregularity by decomposing fluctuating TEC

Post-sunset ionospheric irregularities are common features of the equatorial ionosphere that affect radio communication and navigation systems; their triggering physical mechanism is not yet fully understood. Atmospheric gravity wave is considered as a seeding mechanism for the occurrence of ionospheric irregularities (Abdu et al., 2009). To understand the effects of atmospheric waves, characteristics of wavelike oscillation from ionospheric total electron content (TEC) fluctuation that can be obtained from superposition of different oscillation modes have been investigated. Decomposing fluctuating TEC into different oscillation modes and investigating oscillation characteristics of each component is also important to get insight about the characteristics of individual atmospheric waves that may cause TEC fluctuation. In this paper we have investigated characteristics of components of fluctuating TEC obtained from SCINDA GPS receiver installed at Bahir Dar, (geographic coordinate, 11.5°N, 37.6° E, and dip latitude of 2.5°N) Ethiopia during April 2012. First Empirical Mode Decomposition (EMD) has been applied to decompose TEC fluctuation into different oscillation modes that are known as Intrinsic Mode Function (IMF). Hilbert-Huang Transform (HHT) and Continuous Wavelet Transform (CWT) have been applied to investigate the characteristics of wave-like oscillations. Applying EMD on fluctuating vTEC corresponding to a GPS satellite, five components are found. Results from HHT and CWT have shown excellent agreement. In addition, it is found out that the median periods of oscillation of those five components are 9, 17, 47, 78, and 118 min. Of these periods, 17 and 47 min respectively are oscillation periods of components of TEC fluctuation with occurrence frequency of 92% and 91% that may be interpreted as the manifestation of two frequently occurring components of atmospheric gravity waves that are likely generated by the motion of solar terminator.     

Role of Equatorial Anomaly in Earthquake time precursive features: A few strong events over West Pacific zone

The earthquake (EQ) time coupling processes between equator-low-mid latitude ionosphere are complex due to inherent dynamical status of each latitudinal zone and qualified geomagnetic roles working in the system. In an attempt to identify such process, the paper presents temporal and latitudinal variations of ionization density (foF₂) covering 45°N to 35°S, during a number of earthquake events (M?>?5.5). The approaches adopted for extraction of features by the earthquake induced preparatory processes are discussed in the paper through identification of parameters like the ‘EQ time modification in density gradient’ defined by δ?=?(foF_{2 max}???foF_{2 min})∕τ_mm, where τ_mm – time span (in days) between EQ modified density maximum and minimum, and the Earthquake time Equatorial Anomaly, i.e. EEA, one of the most significant phenomenon which develops even during night time irrespective of epicenter position. Based on the observations, the paper presents the seismic time coupling dynamics through anomaly like manifestations between equator, low and mid latitude ionosphere bringing in the global Total Electron Content (TEC) features as supporting indices.     

Short term variations of total electron content (TEC) fitting to a regional GPS network over the Kingdom of Saudi Arabia (KSA)

The ionosphere is a dispersive medium for radio waves with the refractive index which is a function of frequency and total electron content (TEC). TEC has a strong diurnal variation in addition to monthly, seasonal and solar cycle variations and small and large scale irregularities. Dual frequency GPS observations can be utilized to obtain TEC and investigate its spatial and temporal variations. We here studied short term TEC variations over the Kingdom of Saudi Arabia (KSA). A regional GPS network is formed consisting of 16 sites in and around KSA. GPS observations, acquired between 1st and 11th February 2009, were processed on a daily basis by using the Bernese v5.0 software and IGS final products. The geometry-free zero difference smoothed code observables were used to obtain two hour interval snapshots of TEC and their RMS errors at 0.5 × 0.5 degree grid nodes and regional ionosphere models in a spherical harmonics expansion to degree and order six. The equatorial ionized anomaly (EIA) is recovered in the south of 20°N from 08:00 to 12:00 UT. We found that day-by-day TEC variation is more stable than the night time variation.     

The equatorial ionization anomaly at the topside F region of the ionosphere along 75°E

Electron density measured by the Indian satellite SROSS C2 at the altitude of ∼500 km in the 75°E longitude sector for the ascending half of the solar cycle 22 from 1995 to 1999 are used to study the position and density of the equatorial ionization anomaly (EIA). Results show that the latitudinal position and peak electron density of the EIA crest and crest to trough ratios of the anomaly during the 10:00–14:00 LT period vary with season and from one year to another. Both EIA crest position and density are found to be asymmetric about the magnetic equator and the asymmetry depends on season as well as the year of observation, i.e., solar activity. The latitudinal position of the crest of the EIA and the crest density bears good positive correlation with F_10.7 and the strength of the equatorial electrojet (EEJ).     

Variability study of ionospheric total electron content at crest of equatorial anomaly in China from 1997 to 2007

The variation of TEC data at Wuhan station (geographic coordinate: 30.5°N, 114.4°E; geomagnetic coordinate: 19.2°N, 183.8°E) at crest of equatorial anomaly in China from January 1997 to December 2007 were analyzed. Variability with solar activity, annual, semiannual, diurnal and seasonal variation were also analyzed. The MSIS00 model and ISR model were used to analyze the possible mechanisms of the variabilities found in the results. The TEC data in 1997 and 2001 deduced from another crest station Xiamen (geographic coordinate: 24.4°N, 118.1°E; geomagnetic coordinate: 13.2°N, 187.4°E) were used to contrast. Analysis results show that long-term variations of TEC at Xiamen station are mainly controlled by the variations of solar activities. Typical diurnal variation behaves as a minimum of the TEC in the pre-dawn hours around 05:00–06:00LT and a maximum on the afternoon hours around 13:00–15:00LT. Some features like the semiannual anomaly and winter anomaly in TEC have been reported. The anomaly may be the result of common action of the electric field over the magnetic equatorial and the [O/N₂] at the crest station.     

Highlights about the performances of storm-time TEC modelling techniques for low/equatorial and mid-latitude locations

A statistical evaluation of storm-time total electron content (TEC) modelling techniques over various latitudes of the African sector and surrounding areas is presented. The source of observational TEC data used in this study is the Global Navigation Satellite Systems (GNSS), specifically the Global Positioning Systems (GPS) receiver networks. For each selected receiver station, three different storm-time models based on empirical orthogonal functions (EOF) analysis, non-linear regression analysis (NLRA) and Artificial neural networks (ANN), were implemented. Storm-time GPS TEC data used for both development and validation of the models was selected based on the storm criterion of $\mathit{Dst}?-50$ nT or $K_p?4$ to take into account both coronal mass ejections (CMEs) and co-rotating interaction regions (CIRs) driven storms, respectively. To make an independent test of the models, storm periods considered for validation were excluded from datasets used during the implementation of the models and results are compared with observations, monthly median values, and International Reference Ionosphere (IRI-2016) predictions. Considering GPS TEC as reference, a statistical analysis performed over six storm periods reserved for validation revealed that ANN model is about 10%, 26%, and 58% more accurate than EOF, NLRA, and IRI models, respectively. It was further found that, EOF model performs 15%, and 44% better than NLRA, and IRI models, respectively, while NLRA is 25% better than IRI. On the other hand, results are also discussed referring to the background ionosphere represented by monthly median TEC (MM TEC) and statistics are provided. Moreover, strengths and weaknesses of each model are highlighted.     

Using GPS-SCINDA observations to study the correlation between scintillation,total electron content enhancement and depletions over the Kenyan region

This paper presents the first results of total electron content (TEC) depletions and enhancement associated with ionospheric irregularities in the low latitude region over Kenya. At the low latitude ionosphere the diurnal behavior of scintillation is driven by the formation of large scale equatorial depletions which are formed by post-sunset plasma instabilities via the Rayleigh–Taylor instability near the magnetic equator. Data from the GPS scintillation receiver (GPS-SCINDA) located at the University of Nairobi (36.8°E, 1.27°S) for March 2011 was used in this study. The TEC depletions have been detected from satellite passes along the line of sight of the signal and the detected depletions have good correspondence with the occurrence of scintillation patches. TEC enhancement has been observed and is not correlated with increases in S₄ index and consecutive enhancements and depletions in TEC have also been observed which results into scintillation patches related to TEC depletions. The TEC depletions have been interpreted as plasma irregularities and inhomogeneities in the F region caused by plasma instabilities, while TEC enhancement have been interpreted as the manifestation of plasma density enhancements mainly associated with the equatorial ionization anomaly crest over this region. Occurrence of scintillation does happen at and around the ionization anomaly crest over Kenyan region. The presence of high ambient electron densities and large electron density gradients associated with small scale irregularities in the ionization anomaly regions have been linked to the occurrence of scintillation.