1. Introduction
The ionosphere is the medium of radio wave propagation. In addition to the regular variations in the background ionosphere, there are some special phenomena and structures in the ionosphere, such as sporadic E (Es), equatorial plasma bubble, spread F, etc. They can cause scintillation and loss of locks of the radio wave signals propagating through the ionosphere and seriously affect the reliability of communication and navigation, positioning, and other high-tech systems relying on the trans-ionospheric propagation services [
1,
2,
3,
4]. Es layers are thin and highly dense plasma layers appearing in the lower ionosphere. They can become denser than the normal E layer or even the peak F layer [
5]. Es layers were studied extensively for decades using rocket and satellite in situ measurements, as well as incoherent scatter radar (ISR) and ionosonde observations [
6,
7,
8,
9].
Rocket measurements showed that the Es layer was mainly comprised of metallic ions. The dominant E-region ion O
+ was not detected within the Es layer during the flight of the rocket [
6,
7]. Ground-based observations could trace the vertical motion of the Es layer. ISR is a powerful instrument to monitor the ionospheric structure with superb sensitivity and good range and time resolution. Statistical results from 14 years of Arecibo ISR observations showed that the Es layer generally descended in a diurnal and semidiurnal periodicity [
8]. Based on ionosonde observation at Chung-li, an equatorial anomaly station, the occurrence and vertical motion of higher Es layers were studied. Statistical results indicated that Es layers appeared primarily during the daytime in the spring/winter months. The layers descended from high to low altitudes and different tidal motions controlled the layers in different seasons [
10]. At mid-latitude, a method of height-time-frequency analysis was applied on the ionogram recorded at Milos (36.7°N, 24.5°E) to study the vertical motion of the Es layers with high temporal resolution ionograms (5 min/2 min during the daytime/nighttime) [
9]. The Es layer was characterized by a daytime layer at 120 km and a nighttime layer at above 125 km, of which both moved downward but with a higher descending rate during nighttime. This method was recently applied to study the daily and seasonal variability of Es layers over three European stations. A 12 h periodicity prevailed throughout the year, while 6 h and 8 h periodicities were also found in winter and summer, respectively [
11]. In the Brazil sector, the occurrence frequency of the descending Es layer is >60% over São Luís and >90% over Cachoeira Paulista. In most cases, the layers occur during the day at altitudes 130–180 km and descend to lower altitudes [
12].
It is now generally accepted that the Es layers are majorly composed of metallic ions through a wind shear mechanism [
5,
13]. The long-lived metallic ions move vertically and converge into thin and dense plasma layers when meridional winds are poleward (upside) and equatorward (downside) and zonal winds are eastward (downside) and westward (upside) in the northern hemisphere through the combined action of ion-neutral collisional coupling and geomagnetic Lorentz forcing [
5]. The atmospheric tides primarily produce the vertical shears of horizontal winds and cause the Es layers to descend with the downward propagations of the tide phases. The ions composing Es layers are suggested to be metallic (monoatomic) ions of a meteoric origin undergoing slow radiative recombination, since Es layers live a long time while ambient molecular ions live a short time due to dissociative recombination [
5,
13].
For the observed descending Es layers, their initial heights were much higher than the general heights of meteor ablation, 90–110 km. The results raise a question: what is the source of these high altitude Es layers and how do they appear at such high altitudes? Hanson and Sanitani [
14] suggested a possibility that the E × B drift driven by the equatorial dynamo electric field was the active agent, which lifts the E region metallic ions to higher altitudes around the magnetic equator. For the lifted metallic ions around the magnetic equator, the fountain effect takes over and transports the ions, via gravity and diffusion, to higher latitudes along the magnetic field lines. The deposit of metallic ions from the equatorial F region to low latitudes via the fountain effect provides the metallic ions in the thermosphere. Later numerical simulations confirmed the electric field control in the formation of high altitude metallic ions [
15,
16,
17,
18].
Whereas the previous theoretical and model simulation studies have shown that the high altitude metallic ions are sourced from the equatorial E region, there is almost no direct observational evidence of the upward transport of metallic ions from the E region to F region. To our best knowledge, only two cases of ascending thin layer events were briefly reported, at Jicamarca [
19] and São Luís [
12]. However, the ascending rates of thin layers were quite fast, and the whole evolving process of the thin layers was not very clear due to the low time resolution of the ionograms in these two cases (30 min/10 min for the Jicamarca case and São Luís case, respectively). In this context, an ionosonde with the capability to detect ionospheric variations with high time/height resolution could provide more insight to the vertical motion of the thin layer.
In this paper, a large database of rapid-run ionosonde observations at low latitude Sanya (18.3°N, 109.6°E; dip lat. 12.2°N) during the periods December 2017–June 2018 and April 2019–June 2020 are investigated. Together with the equatorial ionosonde measurements at Chumphon (10.7°N, 99.4°E; dip lat. 3.8°N), we report thin layers ascending from the E region to the F region, simultaneously observed at low latitude and the magnetic equator for the first time. A statistical analysis shows that the low latitude ascending ion layer is not a rare phenomenon. In addition to the equatorial fountain effect, our results suggest a new possibility that the low latitude F region metallic ions could be directly from the local E region. This may have important implications for understanding the sources of the high altitude metallic layers reported previously.
3. Results and Discussion
Figure 1 presents a sequence of ionograms extracted from the 1 min resolution ionograms on 4 April 2018 (LT = UT + 7 h 18 min), a geomagnetic quiet day with a sum of Kp 8.3. Similar to the ionogram scaling in previous studies [
23], here the ordinary wave traces were used to determine the electron density profile. As shown in the figure, two types of Es layers, the low Es in the virtual height of ~110 km (indicated by L) and the high Es in the virtual height of 145–190 km (indicated by H) occurred in the ionogram at 0230 UT (0948LT). We focus on the high Es layer, which shows obvious upward drift in later hours and is, thus, named the ascending ion layer. During 0245–0314 UT, the shape of the ascending layer showed obvious changes, with an evident cusp structure. The minimum virtual height of the structure increased to 200 km at 0314 UT, being connected with the F1 layer trace. In order to present the movement, we marked the position of the high Es layer with a black arrow pointing to the turning point in F1 layer trace. The true height of the ascending ion layer can be calculated by the electron density profile inversion algorithm (NHPC) embedded in the SAOExplorer software [
24]. As indicated by the arrows in panel g, the virtual height and true height of the ion layer are about 217 km and 141 km, respectively. The ion layer continued ascending from a virtual height of 217 km at 0330 UT to 292 km at 0400 UT and finally disappeared at 0432 UT (1150 LT). The frequency at the turning point labeled with the black arrow also increased from 3.9 MHz to 4.4 MHz during 0307–0400 UT. This frequency, however, is actually the combined effect of background electron density and the Es layer density. During 0300–0400UT, it appeared that the F1 layer disappeared between 0314–0330UT. Actually, the fact is that the upper part of the F1 layer trace, which is higher than the high Es layer, became steeper as the high Es layer ascended. This made the F1 layer trace appear to break off due to the limited frequency resolution of the ionogram. An animation of ionograms (S1) obtained during the period 2200 UT, 3 April-0500UT, 4 April is provided in the data repository of this paper at
https://doi.org/10.12197/2022GA006 (last accessed 25 October 2022). It clearly shows the detailed evolution of the ascending ion layer. At the nearby site Fuke (19.5°N, 109.1°E; dip lat. 13.4°N), the ascending ion layer was also observed in the 15 min resolution ionograms (
Figure S2 in data repository).
In order to obtain the whole evolution process of the ascending ion layer, the ionograms were manually scaled with the SAOExplorer software [
23]. The ordinary wave traces and ionospheric parameters (e.g., virtual height, true height) were retrieved. As shown in
Figure 2a, the ion layer (indicated by the red circles) started to ascend at 0210 UT and reached the peak virtual height at 0425 UT. During this period, the F2 layer height also showed an ascending trend.
Figure 2b shows the summary plot of the manually scaled ordinary wave trace with a time resolution of 1 min (different plasma frequencies are color-coded) during the period 2200 UT, 3 April–0500 UT, 4 April (0518–1218 LT, 4 April). The virtual height of the ascending ion layer is marked by red dots. A low Es layer appeared steadily in the virtual heights of 105–120 km with plasma frequencies of 1.5–3 MHz. The high Es layer initially appeared at 2240 UT, 3 April (0558 LT of 4 April), with a critical frequency of 1.6 MHz. It ascended from the upper E region to the F region, shown as a rising structure in the plot. Before ~0200 UT, the virtual height of the ion layer did not change much, except around 0100 UT when a slight enhancement was observed. The variation of the F layer virtual height around the plasma frequencies 3–4 MHz also showed a slight enhancement, similar to that of the high Es layer. After 0200 UT, the ion layer started to move upward, connected to the F1 layer trace at 0310 UT, and reached the maximum virtual height of 510 km at 0422 UT when it moved to the F1 ledge. The virtual height of the ion layer decreased as it continued moving to the F2 layer. Finally, the ion layer disappeared in the F2 layer trace at 0432 UT. During the ascending of the ion layer, the F layer virtual height (around the plasma frequencies 4–6 MHz) also increased slowly.
The ascending ion layer at low latitude Sanya, from the upper E to F region, is similar to the rising structure observed by the Jicamarca radar at the magnetic equator [
25,
26]. It was suggested that the equatorial rising structure was due to the fringe electric fields associated with the initial development of equatorial spread-F (ESF). This, however, cannot explain the present observations because the ascending ion layer was observed during daytime, when generally no ESF was generated. Based on the Arecibo ISR observations, high Es layers with downward movement were frequently observed in daytime and suggested to be controlled by semi-diurnal tides [
8]. Recent simulations supported that the day to day variations in the daytime high Es layer were generated mainly by the vertical shears of zonal winds driven by the semi-diurnal tides [
27]. In the present study, the high Es layer (ascending ion layer) initially appeared at 0558 LT at the virtual (true) height of 125 km (115–120 km) over Sanya. As pointed out by Haldoupis [
5], at the altitudes below 125 km, the zonal wind shear mechanism is dominant on the Es layer formation. Therefore, the generation of the high Es layer observed over Sanya was most likely driven by the zonal wind shear. However, the uplift of the ion layer in later hours was unlikely to be controlled by the tidal wind, because the zonal wind shear node driven by tidal winds normally shows a downward propagation phase [
8,
27,
28].
Based on the feature that the uplift of the ion layer was in a pattern similar to that of the F region, we surmise that the uplift could be mainly driven by a background electric field. In this regard, the measurements by an equatorial ionosonde may provide evidence to support this argument. Fortunately, the ionosonde at Chumphon, an equatorial station in this longitude sector, recorded continuous ionograms on that day (part of the ionograms are shown in
Figure S3 in data repository).
Figure 2c displays the summary plot of manually scaled ordinary wave trace with a time resolution of 5 min at Chumphon. Two interesting features can be noted. One is that an ascending ion layer was also observed at the magnetic equator around the same period, which started from ~0200 UT, connected to the F region bottomside at ~0300 UT, and then disappeared in the F2 layer trace. This result indicates that the ascending of the ion layer occurred simultaneously over a large latitude region from the magnetic equator to the low latitude (>12° dip latitude) far away from the magnetic equator. The other is that at the magnetic equator, the uplift of the ion layer was also in a pattern similar to the vertical drift of the F region, providing a solid evidence for the driving of the ion layer uplift by the background eastward electric field.
Figure 3 shows the true height variations of the ascending ion layer and the F layer peak height (hmF2) at Sanya and Chumphon during 0000–0600 UT on 4 April 2018. Note that for the ion layer, only when its trace connected with the F layer trace could the true height be retrieved by the NHPC algorithm. It can be seen from
Figure 3 that over Chumphon, the F layer peak height (hmF2) increased roughly from 263 km to 312 km during 0215–0255 UT, showing an average upward drifting velocity of ~20 m/s. During 0240–0250 UT, the ion layer height varied from 140 km to 150 km, indicating an average ascending velocity of ~17 m/s. Over Sanya, the vertical drifts of hmF2 and the ion layer during 0310–0415 UT were estimated to be ~10 m/s and ~8 m/s, respectively. It is evident that the ion layer and F region show similar upward drift at both equatorial and low latitudes, which was suggested to be driven by daytime background eastward electric field. Furthermore, there are two notable differences. One is that the ascending velocity of the ion layer is slightly lower than that of hmF2. This could be attributed to higher ion-neutral collision frequency at a lower height, which reduces the vertical drift velocities of the ion layer in the upper E region to the lower F region [
14]. The other is that the upward drifts of hmF2 and the ion layer at low latitude are approximately half of those at the magnetic equator. This may be due to the fact that the vertical drifts driven by electric fields (
E ×
B/B
2) are strongly dependent on latitude. Due to the decrease in electric field strength and the increase in magnetic field inclination with latitude, the vertical drift velocity decreases correspondingly in the poleward direction [
16].
Lee et al. [
19] reported an ascending thin layer event during the morning hour in the equatorial F region at Jicamarca. The thin layer, with topside and bottomside scale heights of 2.8–5.0 and 1.2–2.4 km, was lifted upward from the F1 layer to the F2 layer with a velocity of 15.5 m/s. The ascending velocity of the thin layer observed at Chumphon is similar to that at Jicamarca. The Jicamarca ascending thin layer case provided evidence for the vertical transport of metallic ions in the equatorial F region. Santos et al. [
12] reported another case of an ascending thin layer at São Luís in the American sector. The thin layer ascended from 150 km to 420 km (both virtual height) in 1.5 h. Our results on ascending velocity of ion layer at Chumphon is in line with those from the American equatorial region (Jicamarca and São Luís). All the ascending layer events occurred in morning hours in the spring equinox season (October in Southern American and April in East Asia). Moreover, the present observations showed the clear ascending process of the thin layer from the E region to the F region. However, due to the low signal strength of Chumphon ionosonde data, the uplift of the thin ion layer in the F2 region could not be clearly observed.
By using the OGO 6 satellite measurements, Hanson and Sanatani [
14] have reported the presence of metallic ions above the F2 layer peak. They surmised that the E × B drift driven by the equatorial dynamo electric field was the active agent lifting the E region metallic ions to the F region topside around the magnetic equator. The metallic ions in the equatorial F region topside could then be transported to higher latitudes through the effect of gravity and diffusion (i.e., the equatorial fountain effect). This process was simulated by Carter and Forbes [
16]. The present observations at Chumphon provide direct evidence for the upward moving of the ion layer from the equatorial E region to the F region. Furthermore, the simultaneous observations at low latitude Sanya show that in addition to the equatorial fountain effect, the ion layers could be lifted directly from the local E to F region at the latitudes far away from the magnetic equator, providing another source for the F region metallic ions. This result offers another possibility for the explanation of a metallic ion source responsible for the high altitude metallic layers observed at off-equatorial latitudes. In this regard, one important question is how often this type of ascending ion layer occurs at low latitudes. Based on the daily summary plots (see, for example,
Figure 2b) obtained at Sanya, a statistical analysis was performed. If an ascending ion layer was seen from the summary plot, we then manually checked each ionogram during the interval of the ascending ion layer and obtained its onset time and duration. Note that at Chumphon, some ionograms are not good enough to identify the ascending ion layer, and thus, the statistical feature cannot be well obtained.
Figure 4 displays the statistical results of ascending ion layer occurrence over Sanya during the periods December 2017–June 2018 and April 2019–June 2020. Out of a total of 511 days with data, 38 ascending ion layer events were observed. As shown in
Figure 4a, the occurrence rate of the ascending ion layer shows two peaks in spring (March and April) and quasi-autumn months (August, September). No ascending ion layer event was observed in January and December.
Figure 4b,d show that all the ion layer events started to ascend in the morning hours and ended before 1300 LT, mostly with durations of 0.5–2.5 h. The results show that the ascending ion layer is entirely a morning phenomenon. These results are also in line with the descending upper Es layer trace observed over Arecibo, which mostly descended from 130 km starting at 9–12LT [
8]. On the days with an ascending ion layer, both high and low Es layers appeared, but it was always the high type Es layer that moved upward to the F region (similar to that shown in
Figure 2).
Next, we consider what may control the seasonal variation of the ascending ion layer. As presented above, the background eastward electric field is suggested to be the most likely driving mechanism for the uplift of the thin layer simultaneously observed at Chumphon and Sanya. However, for the statistical results over Sanya, we cannot rule out the possibility of a neutral wind effect. Theoretical simulations have shown that neutral wind can also transport metallic ions upward at low latitudes away from the magnetic equator and at the polar region [
27,
29]. For example, Andoh et al. [
27] reported two cases of an ascending metallic ion layer due to the Lorentz force driven by the eastward wind over Arecibo.
To investigate if there is any statistical relationship between the seasonal variations of the ascending ion layer and of the background eastward electric field, we employed the magnetometer data at the equatorial station Da Lat (11.9°N, 108.5°E, dip lat. 5.3°N) and low latitude station Sanya to calculate the difference in the magnetic field horizontal component (∆H). In previous studies, it was proven that the ∆H provides a direct measure of the daytime electrojet current and, in turn, the magnitude of the vertical E × B drift velocity in the F region ionosphere [
30]. A positive ∆H indicates an eastward electric field. During the periods December 2017–June 2018 and April 2019–June 2020, there were a total of 416 days with simultaneous ∆H and ionosonde data.
Figure 4c shows the variation in average ∆H in high, moderate, and low occurrence months for the ascending ion layer. It is evident that in the months when the occurrence rate of the ascending ion layer is high (M, A, A, S, O months), the background eastward electric field (characterized by the ∆H) is the strongest. In the J, F, N, D months when the occurrence rate of the ascending ion layer is low, the background eastward electric field is the weakest. On the other hand, it is found that the average ∆H on the days with (without) ascending ion layer, the background eastward electric field in the morning hours is relatively strong (weak) (
figures in Supplementary Materials). The close correlation between the strength of the background eastward electric field and the occurrence of the ascending ion layer suggests that a stronger background eastward electric field is more likely to drive the uplift of the ion layer.
Another interesting aspect is that the seasonal variation of the ascending ion layer observed at Sanya is quite consistent with that of the Mg
+ intensity reported previously. Based on the ESRO TD1 satellite observations, Gérard [
31] investigated the annual variation of Mg
+ intensity during sunset hours at the F region topside in 125–160°E. Their results showed two occurrence peaks of Mg
+ in the equinoctial month, with the major and minor peaks in spring and autumn, respectively. This good agreement suggests that the uplifted ion layer from the local E to F regions at low latitudes could contribute significantly to the F region metallic ions observed by satellite-based in situ payloads.
It is worth noting that though the ascending thin layer is likely to be driven by the background electric field, the reason why the ions could be concentrated to maintain a thin layer during their ascending is unknown. Further investigation with new observations and simulations is required. In the observation aspect, lidar observation can provide fine temporal variation in metallic species at given locations [
27,
28,
29], which would shed new light on the study of the ascending high Es layer. However, this is beyond the scope of this study and will be addressed in future studies.