Ionospheric perturbations due to earthquakes as determined from VLF and GPS-TEC data analysis at Agra,India

Employing SoftPAL receiver, amplitude variations of VLF transmitter signals NWC (19.8?kHz) and NPM (21.4?kHz) are analyzed at Agra station in India (Geograph. lat. 27.2°N, long. 78°E) ±15?days from five major earthquakes of magnitude M?=?6.9–8.5 occurred in Indian subcontinent during the years 2011–2013. We apply nighttime fluctuation (NF) method and show that in almost all cases the trend decreases and dispersion and NF increase on the same days corresponding to each earthquake about 11–15?days prior to the main shock. Assuming that the ionospheric perturbations are caused by atmospheric gravity waves (AGW), we also calculate AGW modulation index for each case and find its values increased on the days amplitude fluctuations take place. Its value is decreased in one case only where the perturbations may be attributed to penetration of seismogenic electric field.In order to support the above results we also present GPS-TEC data analyzed by us corresponding to three of the above earthquakes. We study the TEC anomalies (unusual enhancements) and find that in one case the precursory period is almost the same as that found in NF method.     

Assessment of improvement of the IRI model over Ethiopia for the modeling of the variability of TEC during the period 2013–2016

This paper discusses the monthly and seasonal variation of the total electron content (TEC) and the improvement of performance of the IRI model in estimating TEC over Ethiopia during the solar maximum (2013–2016) phase employing as reference the GPS derived TEC data inferred from four GPS receivers installed in different regions of Ethiopia; Assosa (geog 10.05°N, 34.55°E, Geom. 7.01°N), Ambo (8.97°N, 37.86°E, Geom. 5.42°N), Nazret (8.57°N, 39.29°E, Geom. 4.81°N) and Arba Minch (6.06°N, 37.56°E, Geom. 2.62°N). The results reveal that, in the years 2013–2016, the highest peak GPS-derived diurnal VTEC is observed in the March equinox in 2015 over Arba Minch station. Moreover, both the arithmetic mean GPS-derived and modelled VTEC values, generally, show maximum and minimum values in the equinoctial and June solstice months, respectively in 2014–2015. However, in 2013, the minimum and maximum arithmetic mean GPS-derived values are observed in the March equinox and December solstice, respectively. The results also show that, even though overestimation of the modelled VTEC has been observed on most of the hours, all versions of the model are generally good to estimate both the monthly and seasonal diurnal hourly VTEC values, especially in the early morning hours (00:00–03:00?UT or 03:00–06:00?LT). However, it has also been shown that the IRI 2007 and IRI 2012 versions generally perform best in matching the diurnal GPS derived TEC values as compared to that of the IRI 2016 version. In addition, the IRI 2012 version with IRI2001 option for the topside electron density shows the highest overestimation of the VTEC as compared to the other options. None of the versions of the IRI model are proved to be able to capture the effects of geomagnetic storms.     

Comparison of GPS TEC variations with IRI-2007 TEC prediction at equatorial latitudes during a low solar activity (2009–2011) phase over the Kenyan region

This work presents an analysis of the Total Electron Content (TEC) derived from the International GNSS Service (IGS) receivers at Malindi (mal2: 2.9^oS, 40.1^oE, dip −26.813^o), Kasarani (rcmn: 36.89^oE, 1.2^oS, dip −23.970^o), Eldoret (moiu: 35.3^oE, 0.3^oN, dip −21.037^o) and GPS-SCINDA (36.8^oE, 1.3^oS, dip −24.117^o) receiver located in Nairobi for the period 2009–2011. The diurnal, monthly and seasonal variations of the GPS derived TEC (GPS-TEC) and effects of space weather on TEC are compared with TEC from the 2007 International Reference Ionosphere model (IRI-TEC) using the NeQuick option for the topside electron density. The diurnal peaks in GPS-TEC is maximum during equinoctial months (March, April, October) and in December and minimum in June solstice months (May, June, July). The variability in GPS-TEC is minimal in all seasons between 0:00 and 04:00 UT and maximum near noon between 10:00 and 14:00 UT. Significant variability in TEC at post sunset hours after 16:00 UT (19:00 LT) has been noted in all the seasons except in June solstice. The TEC variability of the post sunset hours is associated with the occurrence of the ionization anomaly crest which enhances nighttime TEC over this region. A comparison between the GPS-TEC and IRI-TEC indicates that both the model and observation depicts a similar trend in the monthly and seasonal variations. However seasonal averages show that IRI-TEC values are higher than the GPS-TEC. The IRI-TEC also depicts a double peak in diurnal values unlike the GPS-TEC. This overestimation which is primarily during daytime hours could be due to the model overestimation of the equatorial anomaly effect on levels of ionospheric ionization over the low latitude regions. The IRI-TEC also does not show any response to geomagnetic activity, despite the STORM option being selected in the model; the IRI model generally remains smooth and underestimates TEC during a storm. The GPS-TEC variability indicated by standard deviation seasonal averages has been presented as a basis for extending the IRI-model to accommodate TEC-variability.     

Solar flux effects on equatorial ionization anomaly and total electron content over Brazil: Observational results versus IRI representations

The F layer critical frequency (foF2) as measured by Digisondes in the equatorial and low latitude locations in Brazil is analyzed to investigate the seasonal and solar flux controls of the intensity of the equatorial ionization anomaly (EIA) in the equinoctial month of March. The analysis also included the total electron content (TEC) as measured by a GPS receiver operated at the EIA crest location. The foF2 data set covered a period of large solar flux variation from 1996 to 2003, while the GPS TEC data was for a period in 2002–2003 when the solar flux parameter F10.7 underwent large variations, permitting in both cases an examination of the solar flux effects on these parameters. The seasonal variation pattern in TEC shows a maximum in equinoctial months and a minimum in June solstice, with similar variations for foF2. The solar flux dependence of the TEC is a maximum during equinoxes, especially for post-sunset TEC values at times when the latitudinal distribution is controlled by the equatorial evening plasma fountain processes. Significant variations with local time are found in the degree of solar flux dependence for both the TEC and EIA. The EIA intensity shows large dependence on F10.7 during post-sunset to midnight hours. These results are discussed in comparison with their corresponding IRI representations.     

Geomagnetic storm effects on the occurrences of ionospheric irregularities over the African equatorial/low-latitude region

The study investigated the effects of intense geomagnetic storms of 2015 on the occurrences of large scale ionospheric irregularities over the African equatorial/low-latitude region. Four major/intense geomagnetic storms of 2015 were analyzed for this study. These storms occurred on 17th March 2015 (?229?nT), 22nd June 2015 (?204?nT), 7th October 2015 (?124?nT), and 20th December 2015 (?170?nT). Total Electron Content (TEC) data obtained from five African Global Navigation Satellite Systems (GNSS) stations, grouped into eastern and western sectors were used to derive the ionospheric irregularities proxy indices, e.g., rate of change of TEC (ROT), ROT index (ROTI) and ROTI daily average (ROTI_AVE). These indices were characterized alongside with the disturbance storm time (Dst), the Y component of the Interplanetary Electric Field (IEFy), polar cap (PC) index and the H component of the Earth’s magnetic field from ground-based magnetometers. Irregularities manifested in the form of fluctuations in TEC. Prompt penetration of electric field (PPEF) and disturbance dynamo electric field (DDEF) modulated the behaviour of irregularities during the main and recovery phases of the geomagnetic storms. The effect of electric field over both sectors depends on the local time of southward turning of IMF Bz. Consequently, westward electric field inhibited irregularities during the main phase of March and October 2015 geomagnetic storms, while for the June 2015 storm, eastward electric field triggered weak irregularities over the eastern sector. The effect of electric field on irregularities during December 2015 storm was insignificant. During the recovery phase of the storms, westward DDEF suppressed irregularities.     

Comparison of GPS-TEC measurements with NeQuick2 and IRI model predictions in the low latitude East African region during varying solar activity period (1998 and 2008–2015)

This paper examines the performances of NeQuick2, the latest available IRI-2016, IRI-2012 and IRI-2007 models in describing the monthly and seasonal mean total electron content (TEC) over the East African region. This is to gain insight into the success of the various model types and versions at characterizing the ionosphere within the equatorial ionization anomaly. TEC derived from five Global Positioning System (GPS) receivers installed at Addis Ababa (ADD, 5.33°N, 111.99°E Geog.), Asab (ASAB, 8.67°N, 116.44°E Geog.), Ambo (ABOO, 5.43°N, 111.05°E Geog.), Nairobi (RCMN, ?4.48°N, 108.46°E Geog.) and Nazret (NAZR, 4.78°N, 112.43°E Geog.), are compared with the corresponding values computed using those models during varying solar activity period (1998 and 2008–2015). We found that different models describe the equatorial and anomaly region ionosphere best depending on solar cycle, season and geomagnetic activity levels. Our results show that IRI-2016 is the best model (compared to others in terms of discrepancy range) in estimating the monthly mean GPS-TEC at NAZR, ADD and RCMN stations except at ADD during 2008 and 2012. It is also found that IRI-2012 is the best model in estimating the monthly mean TEC at ABOO station in 2014. IRI show better agreement with observations during June solstice for all the years studied at ADD except in 2012 where NeQuick2 better performs. At NAZR, NeQuick2 better performs in estimating seasonal mean GPS-TEC during 2011, while IRI models are best during 2008–2009. Both NeQuick2 and IRI models underestimate measured TEC for all the seasons at ADD in 2010 but overestimate at NAZR in 2009 and RCMN in 2008. The periodic variations of experimental and modeled TEC have been compared with solar and geomagnetic indices at ABOO and ASAB in 2014 and results indicate that the F10.7 and sunspot number as indices of solar activity seriously affects the TEC variations with periods of 16–32?days followed by the geomagnetic activity on shorter timescales (roughly periods of less than 16?days). In this case, NeQuick2 derived TEC shows better agreement with a long term period variations of GPS-TEC, while IRI-2016 and IRI-2007 show better agreement with observations during short term periodic variations. This indicates that the dependence of NeQuick2 derived TEC on F10.7 is seasonal. Hence, we suggest that representation of geomagnetic activity indices is required for better performance over the low latitude region.     

Anomalous TEC variations associated with the strong Pakistan-Iran border region earthquake of 16 April 2013 at a low latitude station Agra,India

Employing a dual frequency GPS-receiver, ionospheric total electron (TEC) measurements have been in progress at Agra station in India (Geograph. lat. 27.2°N, long. 78°E) since 1 April 2006. In this paper, the TEC data have been analyzed for a period of one month from 1 April-1 May 2013 to examine the effect of multiple earthquakes, some of which occurred on the same day of 16 April 2013, and others occurred in the same month of April, 2013 in India and neighboring countries. We process the data using quartile and epoch analysis based statistical techniques and show that out of all the earthquakes, the one of the largest magnitude (M = 7.8) that occurred on Pakistan-Iran border caused anomalous enhancements and depletions in TEC 1–9 days before the occurrence of main sock. The E × B drift mechanism is suggested for the anomalies to occur in which the seismogenic electric field E is generated in a process suggested by Pulinets (2004).     

A study on TEC reduction during the tail phase of the 21st June 2020 annular solar eclipse

Ionospheric response during the annular solar eclipse of June 21, 2020, has been examined in terms of the Total Electron Content (TEC) obtained from six Global Positioning System (GPS) receivers positioned in the Chinese-Taiwanese region. We have shown TEC variation from satellites designated by PRNs (Pseudo-Random Noise code) 2, 6, and 19. PRN wise TEC trend was observed to depend upon satellite-pass trajectory to the receiver's location during the eclipse period. A time lag of ~15–30 min is also observed in maximum TEC decrement after the phase of maximum eclipse. Instead of the percentage of eclipse magnitude, a reduction in TEC is seen more for the station for which the orbital track of respective satellites was in closer view relative to receivers for more hours of eclipse window. Additionally, the eclipse day diurnal variations are compared with the pre-eclipse day TEC trend, and observed results show a clear decrease in TEC values at all chosen stations after the eclipse onset then reached the lowest value a few minutes afterward the maximum eclipse phase.     

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.     

Spatio-temporal analysis of ionospheric disturbances for ground based augmentation systems over a midlatitude region

Forcings from above and below the ionosphere can cause disturbances that need to be detected and corrected for navigation systems. Ground Based Augmentation Systems (GBAS) are used to give corrections to aircraft navigation systems while landing. These systems use regional ionosphere monitoring algorithms to detect the anomalies in the ionosphere. The aim of this study is to understand occurrence of ionosphere anomalies and their trends over Turkey. A comprehensive analysis of spatio-temporal variability of ionosphere is carried out for a midlatitude GPS network using Slant Total Electron Content (STEC). Differential Rate Of TEC (DROT), which is a measure of the amount of deviation of temporal derivative of TEC from its trend, is used to detect and classify the level of such disturbances. The GPS satellite tracks are grouped into north, east, west and over directions. The 24 h is divided into six time intervals. The percentage occurrence of each DROT category and the deviation from STEC trend in magnitude are calculated and grouped into satellite track directions and time intervals for 2010 (low solar activity), 2011 and 2012 (medium solar activity). The highest level of disturbances is observed in north and west directions, and during sunrise and sunset hours. The dominant periods of percentage occurrences are diurnal (22–25 h), semidiurnal (12–13 h) and terdiurnal (8–9 h) followed by quasi two-day and quasi 16-day periods. Disturbances corresponding to 50% < DROT < 70% are mostly visible during low solar activity years with magnitudes from 1 to 2 TECU. Geomagnetic storms can cause aperiodic larger scale disturbances that are mostly correlated with DROT > 70%. In 2012, the magnitude of such disturbances can reach 5 TECU. The anisotropic and dynamic nature of midlatitude ionosphere is reflected in the spatio-temporal and spectral distributions of DROT, and their percentage occurrences. This study serves a basis for future studies about development of a regional ionosphere monitoring for Turkey.