# Tracking the Ionospheric Response to the Solar Eclipse of November 03, 2013

Authors: Emirant Bertillas Amabayo; Simon Katrini Anguma; Edward Jurua
Journal: International Journal of Atmospheric Sciences
 (2014)

---

## Abstract

The ionospheric dynamics is highly influenced by the solar radiation. During a solar eclipse, the moon occults the solar radiation from reaching the ionosphere, which may drastically affect the variability of the ionosphere. The variability of total electron content (TEC) observed by dual frequency Global Positioning System (GPS) receivers has made it possible to study effects of solar eclipse on the ionosphere. Total eclipse occurred on November 03, 2013, and the maximum amplitude was visible at Owiny in northern Uganda. Ionospheric behavior during this eclipse was analysed by using TEC data archived at Mbarara (MBAR), Malindi (MAL2), Eldoret (MOIU), and Kigali University (NURK) International GPS Satellite (IGS) stations. TEC variations of four consecutive days were used to study instantaneous changes of TEC during the eclipse event. The results generally show TEC decrease at the four stations. However, a maximum perturbation amplitude of ≥20 TECU was observed at MAL2 (18:00–20:00 UT) which is further south of the equator than the other stations. TEC enhancement and depletion were observed during the totality of the eclipse at MOIU, MBAR, NURK, and MAL2 (13:00–15:00 UT). This study found out that the ionospheric TEC over East Africa was modified by wave-like energy and momentum transport and obscuration of the solar disc due to the total solar eclipse.


---

## Body

## 1. Introduction

Solar eclipse occurred on November 03, 2013, and this was observed in different parts of the world. The maximum amplitude of this total eclipse was observed at Owiny in northern Uganda. During solar eclipse, the moon passes between the sun and the earth, thereby occulting solar radiation from reaching the earth when the three celestial bodies are aligned in a syzygy. This may drastically affect the dynamics of the ionosphere. Observation of solar eclipse events has been applied in the field of ophthalmology and optometry to study the risk of solar retinopathy. Psychiatrists have investigated the impact of total solar eclipses on the incidences of suicide. In addition, sociologists investigate the impacts of media hypercoverage, greater social cohesion, tourism, and public disruptions related to the eclipse events. Scientists have investigated the response of the earth’s environment to the abrupt and brief disturbance of solar radiation during the eclipse events [1] (references herein). These environmental effects cover photochemistry, meteorology, gravity waves, boundary layer physics, and the ionosphere [1–4] (references herein). The ionospheric response to solar eclipse events has mainly been investigated using techniques such as rockets, incoherent radar systems, ionosondes, satellite measurements, GPS, and theoretical modeling [5] (references herein).The ionization of gas molecules and/or atoms in the ionosphere by the ultraviolet radiation and X-rays from the sun depends on solar intensity, solar cycle, and seasons [6]. These external factors affect ionization in the ionosphere simultaneously with the internal processes such as chemical reactions amongst the gas molecules and atoms, electrodynamic drifts, and neutral wind effects among others [6]. Since solar eclipse progressively reduces the solar spectral flux, the layers of the ionosphere are equally disturbed. Furthermore, at the lower edges of the ionosphere, recombination coefficient is higher than the recombination rate at higher altitudes. This leads to the lowering of ionization at the lower edges of the ionosphere. Again, the loss rate of ionization at different layers of the ionosphere depends on the chemical composition of each layer [6] (references herein). Further asserted that, in the F1 layer, the NO+ and O+ ions become dominant. A decrease in electron density is observed during an eclipse event in the F1 layer. The decrease in electron density occurs without any time delay as the loss rate of NO+ and O2+ is higher in the F1 region. In the F2 region, a time delay occurs due to slow recombination rate of O+, which is the dominant chemical species in the F2 region.Reference [7] showed that, during an eclipse event, the electron density profile experiences a trough in the eclipse zone. The trough is generated by an increase in the charge order processes expressed as
(1)O++N2⟶NO++NO++O2⟶O2++O
The two charge processes in which O+ ions, in the F2 region, are driven through neutral atmosphere produce NO+ and O2+ as the dominant species in the F1 region. Thus, the above processes are known to cause a faster recovery of foF1 (critical frequency of F1 layer) and delayed recovery of foF2 (critical frequency of F2 layer) after the end of the eclipse. Furthermore, during an eclipse event, the stratification of the F layer into F1 and F2 sublayers becomes distinct. According to [8], this is due to a large upward drift of the F layer to higher altitudes during the solar eclipse. This drift may also slow down the diffusion and hence reduce the ionization of the F2 region and increase the accumulation of ionization in the F1 region, making the stratification more prominent. Besides, the G factor is given by
(2)G=B2aq,
where B is the linear loss coefficient, a is the recombination coefficient, and q is the rate of production of ions and electrons. The G factor quantifies the stratification of the F1 and F2 regions [9]. During solar eclipse, q decreases as the solar disc is obscured and a also decreases as the temperature decreases due to the eclipse shadow region. Hence, the value of G becomes larger around the solar eclipse time, which is attributed to increased stratification of the F layer [10].There is rich literature on the ionospheric variability due to total solar eclipse events. For example, [11] used GPS receivers located at Wuchang and Guangzhou in China to determine the effect of October 24, 1995, solar eclipse on the ionospheric TEC. The study observed a negative deviation in TEC during the eclipse, which is proportional to the degree of obscuration by the lunar disc. Furthermore, the observed maximum negative deviation from the reference levels recorded during the pre- and posteclipse days ranged from 20% to 50%. The absolute maximum negative deviation was found to be 1.0×1017 electrons per square metre at Wuchang, while it was twice at Guangzhou. The maximum deviation was more at Guangzhou due to the fact that it was closer to the center of totality than Wuchang. Therefore, a total eclipse negatively affects the TEC in the ionosphere. In another study, [12] studied the ionospheric response to the total solar eclipse of July 22, 2009, by analyzing the data from four low-latitude GPS stations in India. The ionospheric vertical TEC (VTEC) showed a significant reduction at all the four stations. The authors further showed that the amount of reduction in VTEC decreased as the obscuration decreased. This is a clear manifestation of the direct dependence of electron production on photoionization process. The solar eclipse of July 22, 2009, was again studied by [13]. The authors investigated the behaviour of the ionosphere by measuring the variation in TEC for three (pre-, post-, and eclipse) consecutive days. This eclipse event was complicated by the simultaneous occurrence of a geomagnetic storm which lasted for 10 hours. The results of the GPS data collected from Chongqing and Wuhan and IGS stations at Wuhan and Shao indicated that the instantaneous TEC values declined by about 1–4 TECU at the measurement points depending on the location and time during which the measurement was taken.Furthermore, [5] studied the effect of the solar eclipse of January 15, 2010, and observed that the decrease in the electron density during the eclipse event occurred throughout the E and F1 layers simultaneously. In the F2 layer, the decrease in electron density began at the lower end and extended progressively towards the peak of the F2 region. The experimental measurements and theoretical simulation of results indicate that the eclipse events have more pronounced effects during the midday than in the morning and afternoon hours [5]. They attributed these effects to the fact that a decrease in solar radiation due to the eclipse event resulted in a decrease in the electron production rate and, consequently, the electron concentration. However, the change in the electron concentration in the F2 layer was affected by photochemical processes as well as electrodynamical and neutral forcing processes.Reference [2] described observations of the sporadic E (Es) behaviour during the total solar eclipse of July 22, 2009. Es is a thin layer with dense patches of ionization around the E region. It is generally observed between the altitudes of 95 km and 120 km. This study observed the variations in Es using high frequency ground-based radio systems and found that the critical frequency of the Es (foEs) obtained using ionosonde observations at Wuhan was enhanced. Again, they found that the hop number and the wave power as recorded by oblique incidence rays from Wuhan to Suzhou (both in China) were increased. Furthermore, the wave Doppler frequency shift (DFS) value of Es was increased. Finally, they observed that the cooling effect of the moon’s shadow was believed to induce a powerful meridional air flow in the atmosphere. This was thought to accelerate the ionized clouds in the Es layer and formed the wind shear which raised the observed DFS and foEs values, respectively. During an eclipse event, foEs increases because of the decrease in solar radiation and does not result in a reduction in the electron concentration in Es; instead, the electron concentration increases. This is due to the fact that the meridional air flow accelerates the ionized clouds in Es layer to form the wind shear. It is this wind shear that induces the intensification in the Es layer [3, 4, 14] and increases the foEs during solar eclipse events [2]. In this current study, the ionospheric response to the solar eclipse of November 3, 2013, was investigated using TEC data as measured at four GPS receiver stations in Uganda (MBAR), Kenya (MAL2 and MOIU), and Rwanda (NURK). The next subsection briefly describes this particular solar eclipse event.The solar eclipse of November 3, 2013, was widely referred to as a hybrid solar eclipse due to the occurrence of both annular and total solar eclipse phenomena. The path of lunar umbral shadow moved from the northern Atlantic Ocean through Gabon, Republic of Congo, Democratic Republic of Congo, Uganda, northwestern Kenya, and Ethiopia and ended in Somalia as shown in Figure1.Figure 1
Map showing the path of the solar eclipse event of November 03, 2013 (Courtesy of NASA).The totality of this eclipse lasted for a maximum period of about 1 minute and 39 seconds. In Uganda, the lunar disc was observed to touch the solar disc at around 13:06 UT (16:06 LT), culminating into a total eclipse at 14:22 UT (17:22 LT). The total eclipse had a magnitude of 1.03. The duration of the total solar eclipse at Owiny (2.56°N, 31.42°E) was 22.3 seconds.


## 2. Data Sources and Processing

The data used in this study was archived at four GPS receiver stations shown in Figure2. The geophysical parameters of these stations are shown in Table 1.Table 1
Geophysical parameters of the dual frequency GPS receiver stations.
Station
Code
Geogra. coordinates
Geomag. coordinates
Long. (°E)
Lat. (°S)
Long. (°E)
Lat. (°S)
Malindi
MAL2
40.19
2.30
111.63
6.73
Eldoret
MOIU
35.29
-
0.29
107.64
2.69
Mbarara
MBAR
30.74
0.60
102.99
2.82
Kigali
NURK
30.09
1.94
102.13
4.03Figure 2
Map showing the GPS receiver stations from which data was obtained.The GPS_TEC algorithm developed by [15] was used to derive TEC values at the equatorial stations with an elevation mask of 20°. The GPS_TEC algorithm uses the phase and code values associated with the L1 and L2 GPS frequencies to remove the tropospheric water vapour and clock error effects to calculate relative slant TEC. Absolute TEC values are obtained by taking into account the differential satellite biases and the receiver biases [16]. The equivalent VTEC at altitude of 350 km was calculated by using the thin shell model and averaging the TEC for individual visible satellites. The VTEC derived from the STEC was used to derive the disturbance (VTEC perturbation) in the ionospheric electron density. The data for each GPS station for 2-3 days before and after the eclipse was used as background (quiet day) data to calculate the background mean VTEC (VTEC¯) of the data points for every 30 seconds. The mean VTEC was then subtracted from the disturbed day VTEC to obtain the perturbation (ΔVTEC), similar to [12]:
(3)ΔTEC=VTEC-VTEC¯.
Wavelet analysis was done using the TEC perturbation data and the results are shown in Figure 6. Wavelet analysis provides information on the temporal characteristics of the signal. Spectral analysis performed in this study implemented wavelet transform presented in [17]. The Morlet wavelet was chosen as the mother wavelet because it consists of a plane wave modulated by a Gaussian.


## 3. Results

This section presents results based on TEC measurements from four equatorial GPS stations. The results entail variability of VTEC, VTEC perturbation (ΔTEC), and TEC along signal path (STEC). The results also include wavelet analysis for determining the presence of wave-like structures and traveling ionospheric disturbances (TIDs). The period from 01 to 04 November 2013 was geomagnetically quiet; thus, no results based on geomagnetic data have been presented. The derived VTEC showed a significant reduction with a maximum of ≤35 TECU in Figure 3(b) between 17:00 and 19:00 UT. The reduction in VTEC depends on latitude as well as longitude, which also depends on the location of the satellite with respect to the solar eclipse path [12]. The amount of reduction in VTEC is proportional to the obscuration of the lunar disc, which is directly related to the electron production by the photoionization process [12]. Figures 3(a), 3(b), 3(c), and 3(d) show that the VTEC decreases gradually from around 14:00 UT to minimum values at 17:17 UT, 17:37 UT, 17:37 UT, and 17:24 UT, respectively. The maximum VTEC values are at ~70 TECU in Figures 3(a) and 3(d) and 75 TECU in Figures 3(b) and 3(c). Thus, the decrease in VTEC closely follows the evolutionary pattern of the eclipse. As the obscuration decreased, VTEC started recovering (after 18:00 UT), possibly due to photoionization by partial solar radiation falling into the atmosphere [12]. After 18:00 UT (21:00 LT), solar eclipse ceases and normal photoionization resumed.Variation of VTEC recorded at four equatorial GPS stations from 01 to 04 November 2013. East African standard time, LT = UT + 3 hours.
(a)
(b)
(c)
(d)Several studies of TEC variability [18, 19] have shown higher daytime TEC values during the equinoctial months, moderate values during the summer, and the least values in the winter months. This particular eclipse of November 03, 2013, lies within winter season. During sunrise, the magnetic flux tubes get filled up because of their small volume resulting in sudden increase in ionization due to increasing thermospheric temperatures [18, 19]. The gradual increase in TEC to a maximum value at peak hours of the day at these equatorial latitude stations could be attributed to the solar extreme ultraviolet (EUV) ionization coupled with the upward verticalE ×B drift. The observed nighttime TEC enhancements may be attributed to tidal winds which blow the ionization across geomagnetic fields. Meanwhile, the nighttime decrease could be due to the small size of the magnetic flux tubes that makes the electron content in these tubes collapse rapidly after sunset in response to the low temperatures in the thermosphere at night, leading to low TEC values [18, 19].Generally, the results showed that TEC enhancement was observed from 13:00 to 17:45 UT (refer to Figure4(a)), and later a gradual decrease in TEC. The peak of the TEC enhancement was observed at ~17:43 UT (NURK), during which opposite polarity of VTEC perturbation was observed at MAL2, MOIU, and MBAR. The perturbation levels were low (≤10 TECU) in Figure 4(b), except between 17:30 and 21:30 UT. The TEC depletions observed after the eclipse event in Figure 4(d) were higher than those before (Figure 4(b)) the eclipse event. The study observed a negative deviation in TEC (≥20 TECU) during the eclipse proportional to the degree of obscuration by the lunar disc as shown in Figure 4(c). Furthermore, the observed maximum negative deviation recorded during the pre- and posteclipse days ranges from 5 to 15 TECU. Measurements from the four stations showed reduction in the amount of incoming solar radiation after the eclipse onset. The size of the reduction is closely related to the eclipse magnitude [20]. Changes induced in the spectral solar irradiance during the eclipse resulted in decreased photoionization activity in the ionosphere. The ionospheric response to solar eclipse varies with latitude variations. Figure 4(c) shows that maximum perturbation amplitude of ≥20 TECU (18:00–20:00 UT) was observed at MAL2 which is further south of the equator than the other stations. The same figure shows that TEC enhancement and depletion occurred during the totality of the eclipse at MOIU, MBAR, NURK, and MAL2 (13:00 to 15:00 UT). The different behavior of F2 layer may be due to the fact that electron density is more controlled by the transport process than by the production rate [12]. However, TEC decrease during solar eclipse is dependent on latitude and local time [12] (references herein). The latitudinal variation may arise due to the latitudinal variation of magnetic inclination angle, which influences plasma diffusion. The dependence on local time is attributed to the rate of electron loss, which is related to the local-time dependent background density [12]. However, processes such as equatorial electrojet (EEJ) and the equatorial ionospheric anomaly (EIA) may complicate detection of the impact of solar eclipse in the equatorial region. The oscillatory nature of VTEC perturbation at all stations may be attributed to gravity wave effects generated in the lower atmosphere during the total solar eclipse [12].Variability of VTEC perturbation for two days before eclipse event ((a) and (b)), during eclipse event (c), and one day after the eclipse event (d).
(a)
(b)
(c)
(d)The TEC observed along the ray path between pseudorandom number (PRN) 20 and the GPS receivers at MOIU, MBAR, and NURK shows significant depletions between 17:00 and 20:00 UT (Figure5(a)). A similar TEC pattern (but of lower magnitude) was observed along the ray path between PRN 23 and the receivers at MAL2, MBAR, and NURK between 16:00 and 20:00 UT (Figure 5(b)). TEC enhancements of low amplitude were observed along the ray paths satellite 30 and the four receiver stations between 08:00 and 14:00 UT (Figure 5(c)). The TEC along the ray path between PRN 32 and the four receiver stations revealed depletions between 15:00 and 20:00 UT (Figure 5(d)). These depletions were observed when the satellite was at least above 20° elevation. The amplitude of the TEC depletions is at least 40 TECU for all the four satellites as observed from the four stations. The start of the eclipse was observed by PRN 20 as seen from MBAR (12:56 UT) and NURK (13:21 UT) as shown in Figure 5(a). The eclipse was observed at the low latitude station (MBAR) before that of the slightly higher latitude station (NURK). The effect of the eclipse therefore seems to show latitudinal dependence. This could be attributed to theE ×B vertical plasma drift that may influence the equatorial ionosphere during a solar eclipse. Most likely, the large plasma depletions at the low latitudes get transported to high latitudes through the plasma vertical drift [12] (references herein).TEC observed along the ray paths between the receiver stations and PRN 20 (panel (a)), PRN 23 (panel (b)), PRN 30 (panel (c)), and PRN 32 (panel (d)) on 2013-11-03.
(a)
(b)
(c)
(d)Wavelet transforms of diurnal TEC for one day before ((a) and (b)), after ((e) and (f)), and during ((c) and (d)) the eclipse event.
(a)
(b)
(c)
(d)
(e)
(f)Spectral analysis of the diurnal TEC shows the presence of gravity waves of slightly different periods before, during, and after the eclipse event. The red line shows the cone of influence, which describes the region within which information can be considered significant. Wave-like structures of periods ranging from 1.0 to 1.5 hours were observed at MBAR (06–12 UT and 16–19 UT), MAL2 (06–08 UT and 14–20 UT), MOIU (16–20 UT), and MOIU (10–12 UT), on the day of eclipse. Similar wave-like structures with periods between 1.0 hour and 2.5 hours were also observed at the four stations a day after the eclipse event. A day before the eclipse, high intensity gravity waves having periods of 1-2 hours were observed at MBAR (16–22 UT), MAL2 (20–22 UT), and MOIU (10–12 UT). According to [21], these waves have periods corresponding to large scale TIDs (LSTIDs). TIDs are known to be generated due to the interaction between thermospheric gravity waves and the ionosphere.The TEC perturbations in Figure4(c) consist of contributions from wave-like structures (or TIDs) and obscuration of the solar radiation during the total eclipse event. To isolate contributions to the TEC perturbation due to the solar eclipse, TEC perturbations for November 01, 02, and 04, 2013, were averaged and then subtracted from that of November 03, 2013. The results are shown in Figure 7.Residual diurnal TEC perturbations for one day before (a) and after (b) the eclipse event compared with that of November 03, 2013. Wavelet transforms of the residual TEC perturbations during the total eclipse event ((d), (e), (f), and (g)).
(a)
(b)
(c)
(d)
(e)
(f)
(g)Figures7(b) and 7(c) show that the daily residual TEC perturbations for November 01, 02, and 04, 2013, are approximately zero. Meanwhile, significant residual TEC perturbations were observed on the day of eclipse (November 03, 2013) as shown in Figure 7(a). TEC depletions of ~3 TECU (4–6 UT) and ~5 TECU (16–19 UT) were observed at MAL2, after which TEC started recovering. The other three stations mostly revealed small TEC enhancements of varying amplitude. To further confirm the presence of wave-like structures (TIDs) that can modulate TEC, wavelet analysis was applied on the residual TEC perturbations. The results in Figures 7(d), 7(e), 7(f), and 7(g) confirm that TID contributions to the TEC perturbations in Figure 4(c) have been minimized. It can be assumed here that the TEC perturbations in Figure 7(a) are entirely due to obscuration of solar UV radiation, which led to varying ionization levels in the ionosphere. Hence, ionospheric TEC on November 03, 2013, was modified by the passage of wave-like structures and the obscuration of the solar disc during the solar eclipse event.


## 4. Conclusion

The study investigated the response of the equatorial ionosphere towards the solar eclipse of November 03, 2013. GPS L-band TEC data over the East African region was used for tracking the effect of the eclipse on the ionospheric electron density concentration. The results show that TEC enhancement and depletion were observed during the time of totality of the eclipse at MOIU, MBAR, NURK, and MAL2 (13:00 to 15:00 UT). Maximum TEC perturbation amplitude of ≥20 TECU (18:00–20:00 UT) was observed at MAL2, further south of the equator than the other three stations. TEC depletions (≥40 TECU) were observed along the ray path between PRN 20 (17:00–20:00 UT) and the GPS receivers at MOIU, MBAR, and NURK. TEC enhancements of low amplitude were observed along the ray paths between PRN 30 and the four receiver stations between 08:00 and 14:00 UT. The eclipse started at 12:56 UT and 13:21 UT as observed by PRN 20 from MBAR and NURK, respectively. Generally, the ionospheric TEC on November 03, 2013, was observed to be modified by wave-like energy and momentum transfer and obscuration of the solar disc during the total eclipse.



---

*Source: 127859-2014-10-23.xml*