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.     

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.     

Comparison of GPS derived TEC with the TEC predicted by IRI 2012 model in the southern Equatorial Ionization Anomaly crest within the Eastern Africa region

We have compared the TEC obtained from the IRI-2012 model with the GPS derived TEC data recorded within southern crest of the EIA in the Eastern Africa region using the monthly means of the 5 international quiet days for equinoxes and solstices months for the period of 2012 – 2013. GPS-derived TEC data have been obtained from the Africa array and IGS network of ground based dual-frequency GPS receivers from four stations (Kigali (1.95°S, 30.09°E; Geom. Lat. 11.63°S), Malindi (2.99°S, 40.19°E; Geom. Lat. 12.42°S), Mbarara (0.60°S, 30.74°E; Geom. Lat. 10.22°S) and Nairobi (1.22°S, 36.89°E; Geom. Lat. 10.69°S)) located within the EIA crest in this region. All the three options for topside Ne of IRI-2012 model and ABT-2009 for bottomside thickness have been used to compute the IRI TEC. Also URSI coefficients were considered in this study. These results are compared with the TEC estimated from GPS measurements. Correlation Coefficients between the two sets of data, the Root-Mean Square Errors (RMSE) of the IRI-TEC from the GPS-TEC, and the percentage RMSE of the IRI-TEC from the GPS-TEC have been computed. Our general results show that IRI-2012 model with all three options overestimates the GPS-TEC for all seasons and at all stations, and IRI-2001 overestimates GPS-TEC more compared with other options. IRI-Neq and IRI-01-corr are closely matching in most of the time. The observation also shows that, GPS TEC are underestimated by TEC from IRI model during noon hours, especially during equinoctial months. Further, GPS-TEC values and IRI-TEC values using all the three topside Ne options show very good correlation (above 0.8). On the other hand, the TEC using IRI-Neq and IRI-01- corr had smaller deviations from the GPS-TEC compared to the IRI-2001.     

Statistical study of the time delay of ionospheric TEC storms to geomagnetic storms in Taoyuan,Taiwan

We have studied the time delay of ionospheric storms to geomagnetic storms at a low latitude station Taoyuan (25.02°N, 121.21°E), Taiwan using the Dst and TEC data during 126 geomagnetic storms from the year 2002 to 2014. In addition to the known local time dependence of the time delay, the statistics show that the time delay has significant seasonal characteristics, which can be explained within the framework of the seasonal characteristics of the ionospheric TEC. The data also show that there is no correlation between the time delay and the intensity of magnetic storms. As for the solar activity dependence of the time delay, the results show that there is no relationship between the time delay of positive storms and the solar activity, whereas the time delay of negative storms has weakly negative dependence on the solar activity, with correlation coefficient −0.41. Especially, there are two kinds of extreme events: pre-storm response events and long-time delay events. All of the pre-storm response events occurred during 15–20 LT, manifesting the Equator Ionospheric Anomaly (EIA) feature at Taoyuan. Moreover, the common features of the pre-storm response events suggest the storm sudden commencement (SSC) and weak geomagnetic disturbance before the main phase onset (MPO) of magnetic storms are two main possible causes of the pre-storm response events. By analyzing the geomagnetic indices during the events with long-time delay, we infer that this kind of events may not be caused by magnetic storms, and they might belong to ionospheric Q-disturbances.     

TEC differences for the mid-latitude ionosphere in both sides of the longitudes with zero declination

Based on measurements of ground-based GPS station network, differences of the mid-latitude ionospheric TEC in the east and west sides of North America, South America and Oceania have been analyzed in this paper. Results show that for nearly all seasons from 2001 to 2010 and in both sides of the longitudes with zero declination, there exist systematic differences for the mid-latitude ionospheric TEC in the regions mentioned above and the features of these differences markedly depend upon the local time but less depend upon seasons and the level of solar activity. Theory analysis shows that the longitude variations of both declination and zonal thermospheric winds are one of important factors to cause differences of the mid-latitude ionospheric TEC in both sides of the longitudes with zero declination.     

L-band nighttime scintillations at the northern edge of the EIA along 95°E during the ascending half of the solar cycle 24

The characteristics of nighttime ionospheric scintillations measured at the L-band frequency of 1.575 GHz over Dibrugarh (27.5°N, 95°E, MLAT ～ 17°N, 43° dip) during the ascending half of the solar cycle 24 from 2010 to 2014 have been investigated and the results are presented in this paper. The measurement location is within or outside the zone of influence of the equatorial ionization anomaly depending on solar and geomagnetic activity. Maximum scintillation is observed in the equinoxes irrespective of solar activity with clear asymmetry between March and September. The occurrence frequency in the solstices shifts from minimum in the June solstice in low solar activity to a minimum in the December solstice in high solar activity years. A significant positive correlation of occurrence of scintillations in the June solstice with solar activity has been observed. However, earlier reports from the Indian zone (～75°E) indicate negative or no correlation of scintillation in June solstice with solar activity. Scintillations activity/occurrence in solstices indicates a clear positive correlation with Es recorded simultaneously by a collocated Ionosonde. In equinoxes, maximum scintillations occur in the pre-midnight hours while in solstices the occurrence frequency peaks just after sunset. The incidence of strong scintillations (S₄ ≥ 0.4) increases with increase in solar activity. Strong (S₄ ≥ 0.4) ionospheric scintillations accompanied by TEC depletions in the pre-midnight period is attributed to equatorial irregularities whereas the dusk period scintillations are related to the sporadic-E activity. Present results thus indicate that the current location at the northern edge of the EIA behaves as low as well as mid-latitude location.     

2000年7月和2003年10月大磁暴期间东亚地区中低纬电离层的GPS TEC的响应研究

利用武汉电离层观象台研制的GPS TEC的现报方法及现报系统,对东亚地区GPS台网的观测数据进行处理分析,特别对2000年7月14-18日和2003年10月28日至11月1日两次特大磁暴期间的数据进行了对比考察,文中分析了两次磁暴间的电离层响应,得到对应不同磁暴时段电离层TEC的不同变化情况,着重揭示了TEC赤道异常峰的压缩和移动以及赤道异常随时间的压缩—反弹—恢复的过程,并结合高纬电离层的部分响应机制进行了说明,结果显示,两次磁暴期的电离层响应表现出了各自不同的特点,从而反映出因季节变化引起的高纬电离层暴时能量注入的不同而造成的全球性电离层扰动的不同形态,由此看出,磁暴期间电离层TEC的变化直接与太阳扰动发生的时间及其对高纬电离层的耦合有关,若短时期内连续发生多次磁暴,则电离层反应更加复杂,不能简单地当做单一磁暴叠加处理。     

玉树地震前的电离层异常现象分析

分析了玉树地震前地基电离层探测临界频率、GPS TEC和卫星探测原位等离子体参量等多个参数的扰动变化信息, 研究了不同高度异常变化的时空关联性. 分析发现, 在地震前一天的4月13日, 多个电离层参量出现同步扰动异常, 电离层临界频率f₀F₂异常相对滑动中值增大40%, 异常空间上存在从震中东南向西南漂移的特性; GPS TEC异常增强15TECU (1TECU=1016m-2)左右, 分布于震中南部经度15°范围内, 且有明显的磁共轭效应; DEMETER观测的原位氧离子密度N_i(O+) 4月13日为1-4月中最强的一天, 异常分布偏向赤道区, 但仅局限在30°-50°左右的经度范围内. 综合三个参量的异常特征发现, 无论是空间的局地性还是时间上的密切关联均反映这次电离层扰动可能与玉树地震孕育有关. 结合其他观测信息, 进一步探讨了这次地震孕育过程的地震电离层耦合机理.      

Investigation of GPS-TEC measurements using ANN method indicating seismo-ionospheric anomalies around the time of the Chile (Mw = 8.2) earthquake of 01 April 2014

In 94 km NW of Iquique in Chile (19.610°S, 70.776°W) a powerful earthquake of M_w = 8.2 took place at 23:46:47 UTC (20:46:47 LT) on April 01, 2014. Using GPS-TEC (Total Electron Content) measurements, potential unusual variations around the time and location of the Chile earthquake have been detected based on the median and Artificial Neural Network (ANN) methods. The indices D_st, K_p, A_p and F10.7 were used to distinguish pre-earthquake anomalies from the other anomalies related to the solar-geomagnetic activities. Using the median method, striking anomalies in time series of TEC data are observed 4 days before the earthquake at 14:00 and 16:00 UTC. The ANN method detected a number of anomalies, 4 (02:00 and 16:00 UTC) and 13 (24:00 UTC) days preceding the earthquake. The results indicate that the ANN method due to its capability of non linear learning is quite promising and deserves serious attention as a robust predictor tool for seismo-ionospheric anomalies detection.     

Traveling ionospheric disturbances observed at South African midlatitudes during the 29–31 October 2003 geomagnetically disturbed period

This paper presents traveling ionospheric disturbances (TIDs) observations from GPS measurements over the South African region during the geomagnetically disturbed period of 29–31 October 2003. Two receiver arrays, which were along two distinct longitudinal sectors of about 18°-20° and 27°-28° were used in order to investigate the amplitude, periods and virtual propagation characteristics of the storm induced ionospheric disturbances. The study revealed a large sudden TEC increase on 28 October 2003, the day before the first of the two major storms studied here, that was recorded simultaneously by all the receivers used. This pre-storm enhancement was linked to an X-class solar flare, auroral/magnetospheric activities and vertical plasma drift, based on the behaviour of the geomagnetic storm and auroral indices as well as strong equatorial electrojet. Diurnal trends of the TEC and foF2 measurements revealed that the geomagnetic storm caused a negative ionospheric storm; these parameters were depleted between 29 and 31 October 2003. Large scale traveling ionospheric disturbances were observed on the days of the geomagnetic storms (29 and 31 October 2003), using line-of-sight vertical TEC (vTEC) measurements from individual satellites. Amplitude and dominant periods of these structures varied between 0.08–2.16 TECU, and 1.07–2.13 h respectively. The wave structures were observed to propagate towards the equator with velocities between 587.04 and 1635.09 m/s.