1. Introduction
Manifested as irregular electron density or total electron content (TEC) variations, the equatorial plasma bubbles (EPBs) are a phenomenon related to the collisional Rayleigh-Taylor instability (RTI) mechanisms. Although the onset condition remains unclear, the RTI is believed to occur after sunset at low latitudes near magnetic equator and generates plasma depletions which grow up from the bottom side to the upper F layer. While moving upward, EPBs extend along the magnetic flux tube to higher latitudes and often reach the equatorial ionization anomaly crest [
1]. These plasma irregularities disturb radio communication and navigation by causing phase and amplitude scintillations or even signal loss, which is an important issue for space weather [
2]. Furthermore, the study of the RTI and EPBs is relevant to progress in understanding and modeling of ionospheric plasma physics.
The study of EPBs has gone a long way since it was first recorded as equatorial spread F (ESF) with ionosonde [
3]. The term EPBs was proposed when ground-based incoherent scatter radar detected low-density irregularities rising from the bottom side of F layer to an altitude far beyond the ionospheric peak [
4]. These irregularities were directly confirmed by AE-C satellite in-situ measurements to be regions with an abrupt decrease of electron density by two orders of magnitude with tens of kilometers scale [
5]. The optical imaging techniques probed a two-dimensional structure of the plasma bubble [
6]. Effects of scintillations on radio waves from the stars or satellites have been also used to investigate the behavior of EPBs. TEC fluctuation measured with ground-based Global Navigation Satellite System (GNSS) networks can capture occurrence and evolution of EPBs and reveal their longitudinal occurrence characteristics as well [
2,
7]. Spaceborne GNSS radio occultation of Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) complements the ground-based observations of EPBs with its global coverage. Moreover, it is crucial to obtain the vertical structure of the ionosphere [
8]. For this reason, basing on GNSS and radio occultation measurements, a three-dimensional morphology of EPBs was deduced and a tomography method was attempted to reconstruct EPBs in a Brazilian sector [
9,
10].
Recent researches on RTI and EPBs have been driven mostly by space weather concerns to identify the key environmental parameters leading to the presence of EPBs. Comprehensive observation can reveal the structure and evolution of EPBs that are connected with their initial conditions. This paper aims to investigate the characteristics of EPBs in space- and ground-based GNSS data. The scintillation index S4 and electron density of the ionosphere, profiled by spaceborne GNSS radio occultation of COSMIC, can provide new information on their vertical structure, occurrence, and evolution process of EPBs. Together with ground-based observation, three-dimensional aspects of EPBs can be inspected. We focus on low latitude regions at longitude ~110° E in south China in March 2014 when EPBs occurred every day of the month, and therefore a common underlying cause to seed the RTI can be examined. We use GNSS radio occultation of COSMIC to derive the vertical structure of EPBs and relate their characteristics with ground-based GNSS measurements in China. Global distribution of irregularities occurrence is also presented.
Section 2 introduces data and analysis methods. Our results are described in
Section 3, where TEC fluctuation and occurrence from ground-based GNSS observation are displayed in
Section 3.1, as well as temporal variations of electron density and EPBs occurrence profiles, global occurrence of scintillations caused by EPBs, and sporadic E (E
s) from COSMIC are shown in
Section 3.2.
Section 4 presents a discussion of a comprehensive view of EPBs by combining both ground-based and space-based observations and referencing previous research findings. The conclusions are drawn in
Section 5.
2. Data and Methods
Ground-based GNSS observations are provided by the Crustal Movement Observation Network of China (CMONOC). Eighteen GNSS receivers are selected to form a chain-like configuration along 110° E of longitude, covering a range in latitude between 16.8° N and 35.1° N (corresponding to a magnetic latitude between 7.0° N and 25.2° N), as shown in
Figure 1.
With a thin layer ionosphere model at 400 km, the tracks of the ionospheric piercing point (IPP) have latitude coverage from ~12° N to ~40° N. The slant total electron content (TEC) was calculated with carrier phase measurement performed every 30 s for each receiver. To mitigate multipath effects, only data with satellite elevation larger than 30 degrees was used. The rate of change of TEC (ROT) was determined by taking the difference between the slant TECs at two consecutive times. Calculated on a 5 min time window with 11 successive measurements, the rate of TEC index (ROTI) is defined as the standard deviation of ROT and it is used to quantify TEC fluctuation in time and space [
7]. An irregularity encounter is reckoned if more than 20 consecutive ROTIs are larger than a threshold, which is around 0.2 TECU/min; this threshold, that could be different for different receivers, is described as the sum of the average plus 10 times root mean square (RMS) of all ROTIs collected during daytime 6:00 LT to 18:00 LT for a receiver. In this study, to show vertical TEC fluctuation, we use the vROTI, which can be converted from ROTI, analogously as TEC can be converted from slant TEC. An image of vROTI in latitudinal-temporal variation is obtained by binning vROTI into 0.3° × 30 min latitude-local time grid cells. Occurrence of irregularities is defined by the ratio of ROTI larger than the threshold to all ROTI recorded within the grid cell.
The space-based observations we used are electron density (Ne) and scintillation index S4 from COSMIC radio occultation of GNSS. Defined by the normalized RMS of GPS L1 signal intensity fluctuation, S4 is derived from the 50 Hz measurements of L1 amplitude and the output of the S4 profiles is in 1 Hz sampling rate. A large electron density fluctuation tends to cause a large signal intensity fluctuation. S4 > 0.2 is taken as a threshold to identify the ionospheric irregularities. The mean of Ne profiles at different time is obtained by binning the Ne data into 10 km × 1 h grid cell and taking the arithmetic average within ±20° magnetic latitudes. The vertical profiles of S4 in the altitude range of 90–710 km are used to display height variations of irregularities. To examine the global distribution of S4 at different heights, the data was binned into 5° × 10° latitude-longitude grid cells and occurrence of irregularities is defined by the ratio of samples with S4 > 0.2 to all S4 measured within the grid cell.
The ionosphere at low latitude can be related to that at the magnetic equator by magnetic lines. On the basis of the International Geomagnetic Reference Field (IGRF),
Figure 2 shows a sketch of magnetic lines at the longitude of 110° E which can relate the apex height at the magnetic equator (7° N) to the latitudinal coordinate. The magnetic line crossing the ionosphere (at 400 km) at 23° N has an apex height of ~1000 km at the magnetic equator, while the magnetic line crossing the ionosphere at 30° N has an apex height of ~1500 km at the magnetic equator.
4. Discussion
It is understood that the pre-reversal eastward electric field enhancement (PRE) near sunset uplifts the F layer and contributes to the growth rate of the generalized RTI which is responsible for EPBs. During March 2014, EPBs occurred every day up to the latitude of 20° N at ~110° E longitude, as shown in
Figure 3. It is well known that the occurrence of EPBs tends to peak at equinoxes especially at high solar activity. The year 2014 was a year of maximum activity of solar cycle 24, even though the geomagnetic activity was generally low in March (the largest −Dst was 52 nT on March 1 and the largest K-index was 4 on March 25). Under magnetically quiet conditions, higher solar activity implies greater pre-reversal eastward electric field, as well as earlier occurrence and earlier decay of EPBs [
11,
12]. Nishioka et al. reported a 70% occurrence rate in the Asian region in previous solar maximum 2002 (a year during maximum phase of solar cycle 23) when the solar radio flux F10.7 had a similar level to that of 2014 [
13]. Such a 100% monthly occurrence rate of EPBs is very special.
The monthly mean vROTI and EPB-induced occurrence in the 0.3° × 30 min latitude-local time map in
Figure 4 showed average properties of EBPs in March 2014. In the latitude range between 12° N and 19° N, EPBs started to appear at ~19:00 LT. vROTI had high value from ~20:00 LT to ~23:00 LT. The strongest value was observed at ~20:00–21:00 LT around 16° N; while the largest occurrence (90%) was observed around 13° N at ~22:00 LT. From ~23:00 LT, vROTI decreased with time until EPBs disappeared at ~2:00 LT after midnight. From 19° N to 30° N, the occurrence of the irregularities delayed with latitude and existed for a shorter time with increasing latitude. The highest latitude of 10% occurrence reached ~28° N. Referencing
Figure 2, the magnetic field line at 400 km at ~19–28° N has an apex height of ~750–1000 km at the magnetic equator. This suggests that every day near the magnetic equator EPBs reached the altitude of 750 km about 1 h after their occurrence and some days came to the altitude of 1000 km.
The space-based GNSS observations agreed well with the ground-based ones. The vertical structure of the ionosphere background and the occurrence of EPBs at magnetic equator in
Figure 5 showed that after sunset, the occurrence and development of EPBs is accompanied by the uplift of the F layer, giving direct evidence of eastward electric field. EPBs mainly occurred at and above the F peak. The observed S4 scintillation occurrence was generally smaller than those shown by ROTI. The timing and the latitude dependence of the EPBs occurrence rate agreed between the ground-based and space-based GNSS data, however, the evolution of the EPBs seen from ground-based and space-based GNSS receivers was different. Together with the declining of EPBs, the F peak declined from 450 km at 20:00 LT to 300 km at 24:00 LT in the equatorial ionosphere as seen by COSMIC in
Figure 5. This can be attributed to the strong post-reversal westward background electric field in high solar activity [
11,
12]. As GNSS observation reflects line integral effects, the height variation of the irregularities would not change vROTI much. There is no declination of vROTI occurrence at ~20:00 LT in
Figure 4, actually the occurrence for 17–28° N continued to increase after 20:00 LT. The occurrence from 12° N to 19° N started to decrease at ~21:00–22:00 LT and for higher latitudes it started to decrease after midnight. The irregularities in latitudes lower than ~19° N persisted longer than those in latitudes larger than ~19° N.
Furthermore, the spaceborne GNSS also detected scintillation at 100 km to 130 km which indicates the existence of E
s layer at low magnetic latitudes. Looking at lower heights in
Figure 5, it appears that E
s occurred about 20 min earlier than EPBs. E
s occurrence increased with time and reached the largest at ~22:00 LT, 1 h later than EPBs. They disappeared after midnight at about the same time as that of EPBs. Although there are many researches on low to mid latitudes E
s and their interaction with bottom side of F layer to seed RTI [
14], the nighttime E
s at magnetic equator has not often been reported, especially the E
s appeared in the evening and continued to post midnight. Actually, E
s near magnetic equator could be disrupted while the F layer was uplifted to higher altitude, which is associated with competing conditions of zonal wind shear and upward vertical electric field [
15]. With data of COSMIC for 2006 to 2011, E
s was very rare in the magnetic equatorial zone. Model simulation showed E
s layer occurrences in low-latitude and equatorial regions do not correlate well with the zonal wind shear [
16]. The downward vertical electric field resulted from post-reversal westward electric field probably play a role in the E
s formation at the magnetic equator, however, this cannot explain the earlier occurrence of E
s in the evening.
The global distribution of EPBs was confined to equatorial ionization anomaly regions as shown in
Figure 7. Strongest scintillation occurred at magnetic equator and both equatorial ionization anomaly regions, as reported in previous studies [
1]. EPBs were stronger and had higher occurrence over Atlantic, Africa, and Asia. This is in agreement with the result derived from global ground-based GPS receiver networks [
13]. The longitudinal distribution of E
s was similar to that of EPBs. However, E
s distributed in a larger latitudinal range of ±20° magnetic latitudes, and it was stronger at magnetic equator and ~20° away from the magnetic equator. The magnetic field line at ~100 km at 20° magnetic latitude had an apex height of ~400 km, excluding the possibility of this latitude E
s for seeding EPBs.
Although the mechanisms for seeding EPBs and producing E
s are different, their correlation of occurrence in time implies coupling between the E and F layers, and a common factor probably underlies the formation of both irregularities. Gravity wave seeding from below might be a cause for generation of E
s and EPBs in the equatorial post-sunset ionosphere [
17,
18]. In addition, there is a decrease of both EPBs and sporadic E in the equatorial post-sunset ionosphere over the Pacific. These features suggest a connection between EPBs and sporadic E, at the magnetic equator.