The solar eclipse and its associated ionospheric TEC variations over Indian stations on January 15, 2010

An annular solar eclipse occurred over the Indian subcontinent during the afternoon hours of January 15, 2010. This event was unique in the sense that solar activity was minimum and the eclipse period coincides with the peak ionization time at the Indian equatorial and low latitudes. The number of GPS receivers situated along the path of solar eclipse were used to investigate the response of total electron content (TEC) under the influence of this solar eclipse. These GPS receivers are part of the Indian Satellite Based Augmentation System (SBAS) named as ‘GAGAN’ (GPS Aided Geo Augmented Navigation) program. The eight GPS stations located over the wide range of longitudes allows us to differentiate between the various factors induced due to solar eclipse over the equatorial and low latitude ionosphere. The effect of the eclipse was detected in diurnal variations of TEC at all the stations along the eclipse path. The solar eclipse has altered the ionospheric behavior along its path by inducing atmospheric gravity waves, localized counter-electrojet and attenuation of solar radiation intensity. These three factors primarily control the production, loss and transport of plasma over the equatorial and low latitudes. The localized counter-electrojet had inhibited the equatorial ionization anomaly (EIA) in the longitude belt of 72°E–85°E. Thus, there was a negative deviation of the order of 20–40% at the equatorial anomaly stations lying in this ‘inhibited EIA region’. The negative deviation of only 10–20% is observed for the stations lying outside the ‘inhibited EIA region’. The pre-eclipse effect in the form of early morning enhancement of TEC associated with atmospheric gravity waves was also observed during this solar eclipse. More clear and distinctive spatial and temporal variations of TEC were detected along the individual satellite passes. It is also observed that TEC starts responding to the eclipse after 30 min from start of eclipse and the delay of the maximum TEC deviation from normal trend with respect to the maximum phase of the eclipse was close to one hour in the solar eclipse path.     

2009年7月22日日全食期间电离层总电子含量变化研究

利用中国区域内五个GPS台站(一个台站处于日全食区域、四个台站处于日偏食区域)观测数据, 研究2009年7月22日日全食期间电离层总电子含量(TEC)的变化, 结果表明, 日全食期间, 电离层TEC值经历了下降和恢复的过程, 最小TEC相对于最大食偏的时间延迟约为1～10min; 台站测得最小TEC的星下点(IPP)越靠近日全食带TEC下降量越大, 在日食期间武汉站(114.35°E, 30.53°N) TEC相对于各参考日期的TEC, 其平均下降量最大, 达到4.58TECU.      

The effect of total solar eclipse of October 3, 2005, on the total electron content over Europe

GPS observations from EUREF permanent GPS network were used to observe the response of TEC (Total Electron Content) to the total solar eclipse on October 3, 2005, under quiet geomagnetic conditions of the daytime ionosphere. The effect of the eclipse was detected in diurnal variations and more distinctly in the variations of TEC along individual satellite passes. The trough-like variations with a gradual decrease and followed by an increase of TEC at the time of the eclipse were observed over a large region. The depression of TEC amounted to 3–4 TECU. The maximum depression was observed over all stations located at the maximum path of the solar eclipse. The delay of a minimum level of TEC with respect to the maximum phase of the eclipse was about 20–30 min.     

Ionospheric total electron content response to annular solar eclipse on June 21, 2020

The effects of physical events on the ionosphere structure is an important field of study, especially for navigation and radio communication. The paper presents the spatio-temporal ionospheric TEC response to the recent annular solar eclipse on June 21, 2020, which spans across two continents, Africa and Asia, and 14 countries. This eclipse took place on the same day as the June Solstice. The Global Navigation Satellite System (GNSS) based TEC data of the Global Ionosphere Maps (GIMs), 9 International GNSS Service (IGS) stations and FORMOSAT-7/COSMIC-2 (F7/C2) were utilized to analyze TEC response during the eclipse. The phases of the TEC time series were determined by taking the difference of the observed TEC values on eclipse day from the previous 5-day median TEC values. The results showed clear depletions in the TEC time series on June 21. These decreases were between 1 and 9 TECU (15–60%) depending on the location of IGS stations. The depletions are relatively higher at the stations close to the path of annular eclipse than those farther away. Furthermore, a reduction of about ?10 TECU in the form of an equatorial plasma bubble (EPB) was observed in GIMs at ~20° away from the equator towards northpole, between 08:00–11:00 UT where its maximum phase is located in southeast Japan. Additionally, an overall depletion of ~10% was observed in F7/C2 derived TEC at an altitude of 240 km (hmF2) in all regions affected by the solar eclipse, whereas, significant TEC fluctuations between the altitudes of 100 km ? 140 km were analyzed using the Savitzky-Golay smoothing filter. To prove TEC depletions are not caused by space weather, the variation of the sunspot number (SSN), solar wind (V_SW), disturbance storm-time (Dst), and Kp indices were investigated from 16th to 22nd June. The quiet space weather before and during the solar eclipse proved that the observed depletions in the TEC time series and profiles were caused by the annular solar eclipse.     

Changes in total electron content (TEC) during the annular solar eclipse of 15 January 2010

The solar eclipse of 15 January 2010 was an annular eclipse of the Sun with a maximum magnitude of 0.96 at 1.62°N, 69.29°E. To study the effect of this solar eclipse on the ionosphere the GPS data recorded at three different Indian stations Varanasi (Geographic latitude 25°, 16′N, longitude 82°, 59′E), Hyderabad (Geographic latitude 17°, 20′N, longitude 78°, 30′E) and Bengaluru (Geographic latitude 12°, 58′N, longitude 77°, 33′E) have been used to retrieve ionospheric total electron content (TEC). The ionospheric response to this rare event has been studied in terms of GPS-derived TEC observed at all the three Indian stations. A significant reduction in TEC reflected by all PRNs at all the three stations has been observed. The magnitude of the reduction in VTEC compared to quiet mean VTEC depends on latitude as well as longitude. The amount of reduction observed from different satellites (PRN) is different and depends on the location of the satellite from the solar eclipse path.     

2009年7月22日日全食电离层效应观测

首次利用廊坊中频雷达和武汉、嘉兴、廊坊等三站GPS对2009年7月22日日全食电离层效应进行了观测．日食期间,中频雷达D层78km高度上电子密度减小了约67％,电子密度为200cm^-3的高度上升了近10km,GPS／TEC减小了1TECU左右,其变化的最大相位与日食最大相位几乎同步;日食后,观测到周期为2个多小时的电离层扰动现象．     

Ionospheric responses to the 21 August 2017 great American solar eclipse – A multi-instrument study

Herein, we report on the ionospheric responses to a total solar eclipse that occurred on 21 August 2017 over the US region. Ground-based GPS total electron content (TEC) data along with ground-based measurements (Millstone Hill Observatory (MHO) and digital ionosondes) and space-based measurements (COSMIC radio occultation (RO) technique) allowed us to identify eclipse-associated ionospheric responses. TEC data at ~20°, ~30°, and ~40°N latitudes from the west to east longitudes show not only considerable depression but also wave-like characteristics in TEC both in the path of totality and away from it, exclusively on the day of eclipse. Interestingly, the observed depressions are associated with lesser (higher) magnitudes at stations over which the solar obscuration percentage was meager (significant), a clear indication of bow-wave-like features. The MHO observes a 30% reduction in F₂-layer electron densities between 180 and 220 km on eclipse day. Ionosonde-scaled parameters over Boulder (40.4°N, 100°E) and Austin (30.4°N, 94.4°E) show a significant decrease in critical frequencies while an altitude elevation is seen in the virtual heights of the F-layer only during the eclipse day and that decreases are associated with wave-like signatures, which could be attributed to eclipse-generated waves. The estimated vertical electron density profile from the COSMIC RO-based technique shows a maximum depletion of 40%. Relatively intense and moderate depths of TEC depression, considerable reductions in the F₂-layer electron densities measured by the MHO and COSMIC RO-measured densities at the F₂-layer peak, and elevations in virtual heights and reduction in the critical frequencies measured by ionosondes during the eclipse day could be due to the eclipse-induced dynamical effects such as gravity waves (GWs) and their associated electro-dynamical effects (modification of ionospheric electric fields due to GWs).     

Dayside ionospheric uplift during strong geomagnetic storms as detected by the CHAMP,SAC-C,TOPEX and Jason-1 satellites

Ionosphere response to severe geomagnetic storms that occurred in 2001–2003 was analyzed using data of global ionosphere maps (GIM), altimeter data from the Jason-1 and TOPEX satellites, and data of GPS receivers on-board CHAMP and SAC-C satellites. This allowed us to study in detail ionosphere redistribution due to geomagnetic storms, dayside ionospheric uplift and overall dayside TEC increase. It is shown that after the interplanetary magnetic field turns southward and intensifies, the crests of the equatorial ionization anomaly (EIA) travel poleward and the TEC value within the EIA area increases significantly (up to ∼50%). GPS data from the SAC-C satellite show that during the main phase of geomagnetic storms TEC values above the altitude of 715 km are 2–3 times higher than during undisturbed conditions. These effects of dayside ionospheric uplift occur owing to the “super-fountain effect” and last few hours while the enhanced interplanetary electric field impinged on the magnetopause.     

Response of the mid-latitude D-region ionosphere to the total solar eclipse of 22 July 2009 studied using VLF signals in South Korean peninsula

In this paper, we analyze VLF signals received at Busan to study the the D-region changes linked with the solar eclipse event of 22 July 2009 for very short (∼390 km) transmitter–receiver great circle path (TRGCP) during local noon time 00:36–03:13 UT (09:36–12:13 KST). The eclipse crossed south of Busan with a maximum obscuration of ∼84%. Observations clearly show a reduction of ∼6.2 dB in the VLF signal strength at the time of maximum solar obscuration (84% at 01:53 UT) as compared to those observed on the control days. Estimated values of change in Wait ionospheric parameters: reflection height (h′) in km and inverse scale height parameter (β) in km⁻¹ from Long Wave Propagation Capability (LWPC) model during the maximum eclipse phase as compared to unperturbed ionosphere are 7 km and 0.055 km⁻¹, respectively. Moreover, the D-region electron density estimated from model computation shows 95% depletion in electron density at the height of ∼71 km. The reflection height is found to increase by ∼7 km in the D-region during the eclipse as compared to those on the control days, implying a depletion in the Lyman-α flux by a factor of ∼7. The present observations are discussed in the light of current understanding on the solar eclipse induced D-region dynamics.     

2009年7月22日日全食期间电离层参量的变化

利用多个电离层垂测站的数据和IGS-TEC数据资料, 结合日地空间环境指数, 分析了2009年7月22日日全食期间中国地区电离层参量(反射回波最低频率fmin及f0F2和TEC)的变化特征. 结果表明, 日食发生后fmin迅速降低, 日食结束后fmin迅速恢复到正常水平; 在食甚时刻附近, f0F2和TEC出现明显的降低, 显示了明显的光食效应. 日食结束后5~6 h, f0F2和TEC出现不同程度的正扰动, 在驼峰区更明显; 日食结束后9~10 h, f0F2和TEC出现较显著的负扰动. 由于此次日食发生时伴随着中等强度的磁暴和低纬电场穿透等空间天气事件, 给此次日食电离层效应的深入分析带来很大困难.      

Collocated ionosonde and dense GPS/GLONASS network measurements of midlatitude MSTIDs

To analyze midlatitude medium-scale travelling ionospheric disturbances (MSTIDs) over Kazan (55.5°N, 49°E), Russia, the sufficiently dense network of GNSS receivers (more than 150 ground-based stations) were used. For the first time, daytime MSTIDs in the form of their main signature (band structure) on high-resolution two-dimensional maps of the total electron content perturbation (TEC maps) are compared with ionosonde data with a high temporal resolution. For a pair of events, a relationship between southwestward TEC perturbations and evolution of F2 layer traces was established. So F2 peak frequency varied in antiphase to TEC perturbations. The ionograms show that during the movement of plasma depletion band (overhead ionosonde) the F2 peak frequency is the highest, and vice versa, for the plasma enhancement band, the F2 peak frequency is the lowest. One possible explanation may be a greater inclination of the radio beam from the vertical during the placement of a plasma enhancement band above the ionosonde, as evidenced by the absence of multiple reflections and the increased occurrence rate of additional cusp trace. Another possible explanation may be the redistribution of the electron content in the topside ionosphere with a small decrease in the F peak concentration of the layer with a small increase in TEC along the line-of-sight. Analysis of F2 peak frequency variation shows that observed peak-to-peak values of TEC perturbation equal to 0.4 and 1 TECU correspond to the values of ΔN/N equal to 13% and 28%. The need for further research is evident.     

A method of estimating GPS instrumental biases with a convolution algorithm

This paper presents a method of deriving the instrumental differential code biases (DCBs) of GPS satellites and dual frequency receivers. Considering that the total electron content (TEC) varies smoothly over a small area, one ionospheric pierce point (IPP) and four more nearby IPPs were selected to build an equation with a convolution algorithm. In addition, unknown DCB parameters were arranged into a set of equations with GPS observations in a day unit by assuming that DCBs do not vary within a day. Then, the DCBs of satellites and receivers were determined by solving the equation set with the least-squares fitting technique. The performance of this method is examined by applying it to 361?days in 2014 using the observation data from 1311 GPS Earth Observation Network (GEONET) receivers. The result was crosswise-compared with the DCB estimated by the mesh method and the IONEX products from the Center for Orbit Determination in Europe (CODE). The DCB values derived by this method agree with those of the mesh method and the CODE products, with biases of 0.091?ns and 0.321?ns, respectively. The convolution method's accuracy and stability were quite good and showed improvements over the mesh method.     

Total electron content measurements in ionospheric physics

With the advent of modern global networks of dual-frequency Global Positioning System (GPS), total electron content (TEC) measurements along slant paths connecting GPS receivers and satellites at 22,000 km have become the largest data set available to ionospheric scientists. The TEC can be calculated from the time and phase delay in the GPS signal using the GPS Toolkit, but an unknown bias will remain. In addition, UHF/VHF radio beacons on board low-Earth-orbiting satellites can also be used to measure the electron content. However, the TEC measurements are obtained by integrating TEC differences between slant paths, but also contain biases. It is often necessary to use data assimilative algorithms like the Ionospheric Data Assimilation Three-Dimensional (IDA3D), and to treat both GPS- and LEO-beacon TEC measurements as relative data in order to conduct ionospheric studies.     

Comparative analysis of TEC disturbances over tropical cyclone zones in the North–West Pacific Ocean

In this paper, we study ionospheric total electron content (TEC) disturbances associated with tropical cyclones (TCs). The study relies on the statistical analysis of six cyclones of different intensity which occurred in the North–West Pacific Ocean in September–November 2005. We have used TEC data from the international network of two-frequency ground-based GPS receivers and NCEP/NCAR meteorological archive. TEC variations of different period ranges (02–20 and 20–60 min) are shown to be more intense during TC peaks under quiet geomagnetic conditions. The highest TEC variation amplitudes are registered when the wind speed in the cyclone and the TC area are maximum. The intensification of TEC disturbances is more pronounced when several cyclones occur simultaneously. We have revealed that the ionospheric response to TC can be observed only after the cyclone has reached typhoon intensity. The ionospheric response is more pronounced at low satellite elevation angles.     

Comparison of observed TEC values with IRI-2007 TEC and IRI-2007 TEC with optional foF2 measurements predictions at an equatorial region,Chumphon, Thailand

In this research, as part of working towards improving the IRI over equatorial region, the total electron content (TEC) derived from GPS measurements and IRI-2007 TEC predictions at Chumphon station (10.72°N, 99.37°E), Thailand, during 2004–2006 is analyzed. The seasonal variation of the IRI-2007 TEC predictions is compared with the TEC from the IRI-2007 TEC model with the option of the actual F2 plasma frequency (foF2) measurements as well as the TEC from the GPS and International GNSS service (IGS). The Chumphon station is located at the equatorial region and the low latitude of 3.22°N. For a declining phase of the solar cycle (2004–2006), the study shows that the IRI-2007 TEC underestimates the IRI-2007 TEC with the foF2 observation at the nighttime by about 5 TECU. The maximum differences are about 15 TECU during daytime and 5 TECU during nighttime. The overestimation is more evident at daytime than at nighttime. When compared in terms of the root-mean square error (RMSE), we find that the highest RMSE between GPS TEC and IRI 2007 TEC is 14.840 TECU at 1230 LT in 2004 and the lowest average between them is 1.318 TECU at 0630 LT in 2006. The noon bite-out phenomena are clearly seen in the IRI-2007 TEC with and without optional foF2 measurements, but not on the GPS TEC and IGS TEC. The IRI TEC with optional foF2 measurements gives the lowest RMSE values between IRI TEC predicted and TEC measurement. However, the TEC measurements (GPS TEC and IGS TEC) are more correct to use at Chumphon station.     

The GLONASS-M satellite yaw-attitude model

The proper modeling of the satellites’ yaw-attitude is a prerequisite for high-precision Global Navigation Satellite System (GNSS) positioning and poses a particular challenge during periods when the satellite orbital planes are partially eclipsed. Whereas a lot of effort has been put in to examine the yaw-attitude control of GPS satellites that are in eclipsing orbits, hardly anything is known about the yaw-attitude behavior of eclipsing GLONASS-M satellites. However, systematic variations of the carrier phase observation residuals in the vicinity of the orbit’s noon and midnight points of up to ±27 cm indicate significant attitude-related modeling issues. In order to explore the GLONASS-M attitude laws during eclipse seasons, we have studied the evolution of the horizontal satellite antenna offset estimates during orbit noon and orbit midnight using a technique that we refer to as “reverse kinematic precise point positioning”. In this approach, we keep all relevant global geodetic parameters fixed and estimate the satellite clock and antenna phase center positions epoch-by-epoch using 30-second observation and clock data from a global multi-GNSS ground station network. The estimated horizontal antenna phase center offsets implicitly provide the spacecraft’s yaw-attitude. The insights gained from studying the yaw angle behavior have led to the development of the very first yaw-attitude model for eclipsing GLONASS-M satellites. The derived yaw-attitude model proves to be much better than the nominal yaw-attitude model commonly being used by today’s GLONASS-capable GNSS software packages as it reduces the observation residuals of eclipsing satellites down to the normal level of non-eclipsing satellites and thereby prevents a multitude of measurements from being incorrectly identified as outliers. It facilitates continuous satellite clock estimation during eclipse and improves in particular the results of kinematic precise point positioning of ground-based receivers.     

Variations of total electron content over high latitude region during the ascending phase of 24th solar cycle

The solar cycle variation and seasonal changes significantly affects the ionization process of earth’s ionosphere and required to be monitored in real time basis for regional level refinement of existing models. In view of this, the present study has been carried out by using the ionospheric Total Electron Content (TEC) data observed with the help of Global Ionospheric Scintillation and TEC monitoring (GISTM) system installed at Indian Antarctic Research Station, “Maitri” [70°46′00″S 11°43′56″E] during the ascending phase of 24th solar cycle. The daily values of solar extreme ultraviolet (EUV) flux (0.1–50?nm wavelength), 10.7?cm radio flux F_10.7 and Sunspot number (SSN) has been taken as a proxy to represent the solar cycle variation to correlate with TEC. The linear regression results revels better correlation of TEC with EUV flux rather than F_10.7 and SSN. Also, the EUV and TEC show better agreement during summer as compared to winter and equinox period. Correlation between TEC and EUV appears significantly noticeable during ten internationally defined quiet days of each month (stable background geophysical condition) as compared to the overall days (2010–2014). Further, saturation effect has been observed on TEC values during the solar maxima year 2014. The saturation effects are more prominent during the night hours of winter and equinox season due to transportation losses manifested by the equator-ward direction of meridional wind.     

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.     

Storm time spatial variations in TEC during moderate geomagnetic storms in extremely low solar activity conditions (2007–2009) over Indian region

The total electron content (TEC) measurements from a network of GPS receivers were analyzed to investigate the storm time spatial response of ionosphere over the Indian longitude sector. The GPS receivers from the GPS Aided Geo Augmented Navigation (GAGAN) network which are uniquely located around the ∼77°E longitude are used in the present study so as to get the complete latitudinal coverage from the magnetic equator to low mid-latitude region. We have selected the most intense storms but of moderate intensity (−100 nT < Dst < −50 nT) which occurred during the unusually extremely low solar activity conditions in 2007–2009. Though the storms are of moderate intensity, their effects on equatorial to low mid-latitude ionosphere are found to be very severe as TEC deviations are more than 100% during all the storms studied. Interesting results in terms of spatial distribution of positive/negative effects during the main/early recovery phase of storms are noticed. The maximum effect was observed at crest region during two storms whereas another two storms had maximum effect near the low mid-latitude region. The associated mechanisms like equatorial electrodynamics and neutral dynamics are segregated and explained using the TIMED/GUVI and EEJ data during these storms. The TEC maps are generated to investigate the storm time development/inhibition of equatorial ionization anomaly (EIA).