Abstract

The influence of the 15th January 2010 annular solar eclipse on traveling ionospheric disturbances (TIDs) and equatorial plasma bubbles (EPBs) is studied using data from six global navigation and satellite system (GNSS) receivers spread across the path of annularity over the low latitude region of East Africa. The GNSS receivers are stationed at Nairobi (RCMN), Malindi (MAL2), and Eldoret (MOIU) in Kenya; Mbarara (MBAR) in Uganda; Kigali (NURK) in Rwanda; and Mtwara (MTWA) in Tanzania. The study period ranges from 12th to 18th January 2010, three days before and after the 15th January 2010 annular solar eclipse. The year 2010 marked the beginning phase of solar cycle 24, evidently observed in low total electron content (TEC) values and the disturbed storm time index (Dst). The eclipse started at 7 : 06 LT and ended at 10 : 14 LT, with MOIU and RCMN experiencing eclipse magnitudes of 0.946 and 0.93, respectively. The maximum obscuration occurred between 8 : 21 LT and 8 : 34 LT across most of the stations. A detrending on vertical TEC (VTEC) derived from GNSS receivers across or close to the path of totality revealed a reduction of ∼2-3 TECU during the maximum phase of the eclipse. The level of reduction was highly close to the totality path and decreased smoothly away from the totality path. Using a background polynomial fitting technique on diurnal TEC, we analyzed TIDs along NURK-MBAR-MOIU and MOIU-RCMN-MAL2 GPS arrays. The results revealed a wavelike perturbation with a virtual horizontal velocity of 830m/s and ∼1 TECU amplitude propagating eastward along the MOIU-RCMN-MAL2 GPS array. The study reports a moderate scintillation activity of 0.5 ≤ ROTI ≤ 0.9 values, demonstrating the presence of few EPBs over the region. The results show a latitudinal variation in GPS-TEC scintillation activities and suggest a possible influence of the eclipse on the observed increase in average scintillation levels across East Africa.

1. Introduction

A solar eclipse happens when the Moon moves in front of the Sun and blocks its light entirely or partially on a portion of the Earth's surface, causing the solar illuminance changes for a short time [1], which also offers a rare chance to study the ionospheric reaction to the sudden cessation of solar radiation [2, 3]. The total solar eclipse only lasts for a few minutes in any given location due to the supersonic speed of the Moon's shadow traveling eastward across the Earth. Due to various Sun-Moon-Earth alignment geometries and varied lunar orbital properties, it does not repeat in any discernible pattern [4]. The ionosphere's response to a solar eclipse is influenced by a number of variables, including the sun's activity, the geomagnetic environment, the location, latitude, and longitude, and the time of day [3].

The rapid fluctuations in solar electromagnetic flux and ionizing ultraviolet (UV) radiation that occur during eclipses have a variety of impacts on the thermosphere and ionosphere, including a change in the temperature balance and the generation and loss of ionization in the lower ionosphere [5]. The eclipse alters the ionization in the E and F regions of the ionosphere where dynamic processes are crucial. As a result, the ionosphere is likely to be considerably affected both during and after a solar eclipse. Ionospheric response to solar eclipses at low latitudes and the underlying physical mechanisms is still an open question [6]. However, inconclusive results of solar eclipse observations arise from the fact that different solar eclipses produce different plasma motions [7].

Due to recent developments in radio remote sensing technology, the effects of solar eclipses on the ionosphere may now be continuously monitored over a wide area. Studying ionospheric dynamics across East Africa is now possible thanks to the installation of global navigation and satellite systems (GNSS) receivers throughout the region. During a solar eclipse, the total electron content (TEC), produced by the Global Positioning System (GPS) receivers, is a reliable indicator of the condition and dynamics of the ionospheric F-region [8, 9].

Solar eclipses are known to generate atmospheric gravity waves (AGWs) [10, 11]. AGWs induce traveling ionospheric disturbances (TIDs) via neutral-ion collisions [12]. TIDs are evident as fluctuations in TEC derived from GPS receivers [14]. Although these TIDs appear to be self-propagating, they are not. If the underlying AGW dissipates, the TID will disappear because of the lack of neutral-ion collisions. Therefore, TIDs induced by AGWs are not exact tracers of the AGWs. However, their presence does indicate the existence of underlying AGWs, which can be quite useful for diagnosing AGWs and their sources [14]. The source of the AGWs could be located by combining observations from three stations using a reverse ray-tracing technique [3]. A differential total electron content (ΔTEC) approach detects tiny spatially coherent disturbances caused by numerous physical mechanisms [15]. Local conditions such as local time, background wind, and topography might affect AGW propagation, causing significant variation in the ionospheric electron density during an eclipse [16].

A study by Emirant [17] on the ionospheric response to the 3rd November 2013 solar eclipse showed that TEC over East Africa was modified by wave-like energy and momentum transport. When studying detection and description of the different ionospheric disturbances that appeared during the solar eclipse of 21st August 2017 over the United States, Heng [18] observed medium-scale TIDs (MSTIDs) with wavelength and velocity in the range of 250–300 km and 150–300 m/s, respectively, propagating in the same direction as that of the lunar disc during the early stages of the eclipse transit. The MSTIDs were detected less than 30 minutes after the eclipse's onset. Earlier observations during the same eclipse event over North America by Zhang [12] revealed strong solar eclipse-induced ionospheric bow waves with a wavelength ranging between 270 and 350 km, a 25-minute period and a phase velocity of 280 m/s in a direction aligned with totality. In a study to investigate the effects of the 15th January 2010 solar eclipse event on equatorial and low latitude regions of India, Panda [19] observed that the faster downward movement of plasma was the primary cause of increased electron density over the equator and substantially decreased TEC at all latitudes. When investigating the ionospheric response over South America during the 2nd July 2019 total solar eclipse using the Chilean GPS eyeball, [16] observed AGWs with a period of 30 to 60 minutes at the station located north of totality.

The intensity of ionospheric scintillations is frequently described quantitatively using the rate of change of TEC index (ROTI) [20, 21]. Equatorial plasma bubbles (EPBs) are believed to be the dominant source of irregularities that cause scintillation of the GPS satellite signals at L band frequencies [22, 23]. Values of ROTI ≥ 0.5 TECU/min are considered moderate scintillation activities [24], which could potentially reveal the presence of EPBs. Strong scintillation in the signal amplitude and phase contributes to the loss of lock on GPS signals, decreasing the number of available satellites and consequently reducing the possibility to get good satellite geometry, hence degrading navigation and positioning measurements [25]. A study by Norsuzila [26] on the effect of the 15th January 2010 solar eclipse over Malaysia also reported a reduction in TEC and amplitude scintillation. When examining the effects of the 21st August 2017 total solar eclipse on variations of TEC over the North American region, the author found that the solar eclipse decreased ionospheric scintillation occurrence at the totality region by up to 28 percent and TEC values decreased by up to 7 TECU [27].

We report the effects of the 15th January 2010 annular solar eclipse on low-latitude TIDs and EPBs over the East Africa region using data from GNSS receivers stationed at Nairobi (RCMN), Malindi (MAL2), and Eldoret (MOIU) in Kenya; Mbarara (MBAR) in Uganda; Kigali (NURK) in Rwanda; and Mtwara (MTWA) in Tanzania.

2. Materials and Methods

2.1. Preparation and Processing of Data

The 15th January 2010 annular solar eclipse was visible across most parts of the East Africa region, with RCMN and MOIU receiver stations being along the umbra path, while MAL2, NURK, MBAR, and MTWA receiver stations were along the penumbral path of the eclipse. Generally, the eclipse onset was at 07 : 06 LT, reaching its maximum phase at about 08 : 28 LT and ending at 10 : 14 LT. The period of annularity lasted between 6 and 8 minutes. Further details are available at https://eclipse.gsfc.nasa.gov/SEmono/ASE2010/ASE2010.html. The geographic locations of the six GNSS receivers, namely, RCMN, MOIU, MAL2, NURK, MBAR, and MTWA are given in Figure 1.

Table 1 shows the geographic latitude (GLAT), geographic longitude (GLON), geomagnetic latitude (MLAT), and geomagnetic longitude (MLON) of the GNSS receiver stations, obtained from http://wdc.kugi.kyoto-u.ac.jp/igrf/gggm/index.html.

The data used for this research for the period between 12th and 18th January 2010 were obtained from the GNSS receivers stations at six locations over the East African region, as shown in Figure 1. The data for each GNSS receiver stations were downloaded from https://cddis.gsfc.nasa.gov/pub/gps/data/daily/ in receiver independent exchange (RINEX) format. The TEC values for each station were derived using the GPS-TEC algorithm developed by [28]. The application first calculates the relative values of slant TEC (STEC) by eliminating the effect of tropospheric water vapor and clock errors using the phase and code values of the two GPS frequencies [29], then later the absolute value after compensating for the satellite and receiver biases. The absolute STEC along the oblique signal paths are also mapped to vertical TEC (VTEC) values at the ionospheric pierce point (IPP) assumed to be 350 km by means of the single layer model (SLM) mapping function [30, 31]. The VTEC values were sampled at 30-second intervals.

2.2. Traveling Ionospheric Disturbances (TIDs)

To determine the background or unperturbed ionospheric TEC, each satellite's time series of vertical TEC data were fitted with a fourth-order polynomial across all stations. The equation for the fourth-order polynomial fitting iswhere is fitted TEC at time t, i = 1, 2, …, 6 (no. of receiver stations), j = 1, 2, …, 31 (no. of satellites), coefficients a, b, c, and d are obtained through the least squares method, and ε is the residual error of the fitting process [32]. The dominant periods of the structures observed, were obtained using detrended TEC (ΔTEC), expressed aswhere is the actual TEC at time, t, and implies for all i and j receiver stations and satellites, respectively. By monitoring , the amplitude of the wavelike structures (TIDs) that modulate TEC were inferred [32, 33].

The virtual horizontal velocities of the TIDs were calculated by monitoring the time delays between consecutive peaks and troughs of the wavelike structures and determining the horizontal distances between the GNSS receiver stations [34]. The potential wavelike structures were monitored across three stations along similar or close latitudes. The horizontal distances between stations along the NURK-MBAR-MOIU and MOIU-RCMN-MAL2 arrays, as shown in Table 2, were calculated using the following equation:where RE is the radius of the Earth and ∆σ is the central angle between two receivers given aswhere , with the subscript s for the first latitude/longitude and subscript f for the second latitude/longitude [35].

2.3. Equatorial Plasma Bubbles (EPBs)

The rate of change of TEC (ROT) was calculated for every considered GNSS station and then detrended for all individual satellite tracks using the equation:where i is the visible satellite, k is the time of epoch, k−1 is the subsequent time of epoch, and ∆t is the time difference (60 seconds). ROTI, defined as the standard deviation of the detrended ROT values [36], was calculated using the following equation:

The presence of EPBs is inferred by monitoring the intensity of ionospheric scintillations. A value of ROTI ≥ 0.5 TECU/min is considered moderate scintillation activity [24] and could potentially reveal the presence of EPBs.

3. Results and Discussion

3.1. Effect of the Eclipse on TEC

Figure 2 shows the daily mean VTEC over all stations for days 12th to 18th January 2010, with a minimum of 14.19 TECU and a maximum of 19.78 TECU. Only NURK (Figure 2(f)) shows considerable decline in daily mean VTEC values on the eclipse day (15th day) in comparison with the 14th and 16th day. Generally, the eclipse impact on daily mean TEC variation was not so pronounced since the year 2010 marked the beginning phase of the solar minimum of solar cycle 24.

To monitor TEC variations, we obtained the detrended TEC over MOIU, RCM, MAL2, MBAR, NURK, and MTWA receiver stations on 15th January 2010, during the onset, maximum, and offset phases of the eclipse as shown in Figure 3.

In Figure 3, RCMN, MAL2, MBAR, and MOIU stations that were either within or close to the path of totality during the maximum phase of the eclipse showed a considerable decline in TEC of 1-2 TEC units. The moon occulted the sun’s radiation over the East African region, thus briefly preventing photoionization. This was evident, especially during the maximum phase of the eclipse at about 08 : 27 LT. Similar results were observed by [37] when studying the effect of the 15th January 2010 solar eclipse over the Indian region, in which a decline in TEC during the eclipse period was dependent on the station’s proximity to the eclipse path.

3.2. Effect of the Eclipse on TIDs

TIDs normally manifest themselves as plasma density fluctuations which propagate as a wave away from the region of totality. Since TIDs are known to modulate TEC, we investigated TEC enhancements and depletions observed across receiver stations along similar or close latitudes. The GNSS receiver station arrays investigated were NURK-MBAR-MOIU and MOIU-RCMN-MAL2, thus giving longitudinal variation regarding the TIDs propagation. Figure 4 shows the graph of VTEC and fitted VTEC against local time (LT) across NURK-MBAR-MOIU and MOIU-RCMN-MAL2 station arrays on 15th January 2010.

In Figure 4, the TEC perturbations manifesting as a wave-like structure along the arrays are presented as a double peak, indicated by the arrows before and during the eclipse period. The first peak is observed between 4 : 00 LT and 5 : 00 LT, while the second peak is observed between 7 : 00 LT and 8 : 30 LT. Although the first peak is not quite visible along the NURK-MBAR-MOIU array, it is significantly noticeable along the MOIU-RCMN-MAL2 array. For the MOIU-RCMN-MAL2 GPS array, the first peak is observed at MOIU, then RCMN, and finally MAL2 station, with the time of contact being about 06 : 50 LT, 07 : 05 LT, and 07 : 10 LT, respectively. The virtual horizontal velocity of the observed wave-like perturbation was calculated using time delays between the peaks of structures observed and the horizontal distances presented in Table 2. The average virtual horizontal velocity of the wave-like perturbation was found to be 830 m/s in the eastward direction. Since the wavelike structure is observed at the same time as the eclipse and propagates in the same direction as the lunar disc, it could potentially be influenced by the eclipse event. The detrended TEC plots revealing the amplitude of the perturbation along NURK-MBAR-MOIU and MOIU-RCMN-MAL2 receiver array on 15th January 2010 are shown in Figure 5.

In Figure 5, a TEC depletion of about 1.5 TECU was observed along NURK-MBAR-MOIU and MOIU-RCMN-MAL2 receiver array on 15th January 2010, the eclipse day. This is because these stations were either along or close to the path of totality of the eclipse.

Figures 611 show plots of actual and fitted VTEC against local time (LT) across NURK-MBAR-MOIU and MOIU-RCMN-MAL2 station arrays on 12th, 13th, 14th, 16th, 17th, and 18th January 2010, respectively.

We note in Figures 611 that the control days, 12th, 13th, 14th, 16th, 17th, and 18th January 2010, respectively, did not reveal similar wave-like structures. The VTEC perturbations observed across different arrays in Figures 611 did not reveal continuous propagation of any wavelike structures. Thus, there was minimal perturbation from the background VTEC on the control days, 12th, 13th, 14th, 16th, 17th, and 18th January 2010 resulting from an underlying propagating disturbance at ionospheric heights.

3.3. Effects of the Solar Eclipse on Ionospheric Scintillations

The ROTI derived from TEC observations is a good proxy for ionospheric scintillations. The plots of ROTI against LT derived from equation (6) are presented in Figures 1217 for RCMN, MAL2, MOIU, MBAR, MTWA, and NURK stations, respectively.

In Figure 12, the highest ROTI value of 0.9 TECU/min between 20 : 00 LT and 23 : 00 LT was observed over the RCMN station on 15th January 2010. Clearly, there was more intense ionospheric scintillation over RCMN station on the eclipsed day than on the control days. It is important to note that the RCMN station was within the path of totality.

In Figure 13, the MOIU receiver station experienced weak scintillation activities both on the eclipsed day and control days as observed even though it was along the path of totality.

In Figure 14, the MAL2 station experienced two moderate ionospheric scintillation events on 12th and 15thJanuary 2010 with ROTI values of 0.45 TECU/min and 0.5 TECU/min, respectively, occurring between 20 : 00 and 23 : 00 LT.

In Figure 15, the MBAR station had the highest ROTI value of 0.4 TECU/min, occurring between 20 : 00 and 23 : 00 LT on 12thJanuary 2010. The eclipsed day, however, did not reveal significant scintillation activities.

Figure 16 shows that most of the days over MTWA station, including the eclipsed day, had an average ROTI value of 0.1 TECU/min. Therefore, the ionosphere over MTWA station had no major plasma density depletions that could potentially cause strong GPS signal scintillation. There was no observed eclipse effect on EPBs over MTWA. This could be due to its proximity to the path of annularity as compared with other stations.

In Figure 17, NURK station had the highest ROTI value of about 0.5 TECU/min between 20 : 00 and 22 : 00 LT on 12thJanuary 2010. The eclipsed day, however, did not reveal significant scintillation activities.

Generally, as shown in Figures 1217 for RCMN, MAL2, MOIU, MBAR, MTWA, and NURK stations, respectively, there were minimal scintillation events with moderate ROTI values ranging from 0.5 to 0.9 TECU/min observed across all the stations during the period of study. This is in agreement with [38, 39] who reported fewer and moderate GPS scintillations over the same region in January 2010. This implies that there were few plasma density depletions over the East Africa region. Stations closer to the geographic equator, such as RCMN, MBAR, MAL2, and NURK, recorded considerable scintillation events in contrast to MTWA station, which had none. On the eclipse day, 15th January 2010, considerable scintillation events were only observed over RCMN and MAL2 stations, with ROTI values of about 0.9 TECU/min and 0.5 TECU/min, as seen in Figures 12 and 13, respectively. This is because low latitude ionospheric scintillation is known to be a postsunset phenomenon. The influence of the eclipse on TEC earlier in the day was more pronounced over these two stations. However, [27] reported a reduction in scintillation occurrence of up to 28% during the 21st August 2017 solar eclipse over the North American region.

4. Conclusion and Summary

This study presents the first observation of the 15th January 2010 solar eclipse effect on TIDs and EPBs over six stations across the low-latitude regions of East Africa during the ascending phase of the solar minimum of solar cycle 24. The control days were three days before and after the eclipse. TEC depletion was 1.5 TECU along MBAR, MOIU, RCMN, and MAL2 stations. These four stations were either along or close to the totality path of the eclipse. Further analysis revealed the passage of an eastward propagating TID with an average virtual horizontal velocity of 830 m/s along the MOIU-RCMN-MAL2 GPS array during the eclipse period. No TIDs were observed on control days along both GPS arrays. Equatorial plasma bubbles are evident as signal scintillation with values of ROTI ≥0.5. RCMN and MAL2 stations recorded moderate postsunset scintillation levels on the eclipsed day with ROTI values of 0.9 TECU/min and 0.5 TECU/min, respectively, suggesting a possible influence of the eclipse on the occurrence of EPBs, especially over RCMN and MAL2 stations. However, it is not clear why MOIU station being along the totality path showed contrasting observations. Nonetheless, the study reveals evidence of solar eclipse-induced TID propagating eastward in the direction of totality and an increase in scintillation levels over MAL2 and RCMN stations on 15thJanuary 2010. More observations and simulations are necessary to further understand the complex dynamic processes in both the ionospheric and space environment during and after the solar eclipse.

Data Availability

The data that are archived in receiver independent format were downloaded from https://cddis.gsfc.nasa.gov/pub/gps/data/daily/ and used to derive TEC necessary for further analysis.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

The authors sincerely would like to thank the International GNSS Service for the data used in this study, https://cdaw.gsfc.nasa.gov/CME_list/index.html and http://wdc.kugi.kyoto-u.ac.jp/, Prof. Gopi Seemala of Indian Institute of Geomagnetism for the GPS-TEC analysis software, https://eclipse.gsfc.nasa.gov/SEmono/ASE2010/ASE2010.html, and Masinde Muliro University of Science and Technology for library and office space.