3.4.1. Potential Drivers of Equatorial and Mid-Latitude Ionospheric Irregularities
Changes in thermospheric composition are an important driver of positive or negative ionospheric irregularities. Therefore, the O/N
2 density ratio measured by the Global UltraViolet Imager (GUVI) on board the space-borne TIMED satellite (~625 km) was analyzed. It is important to note that the O/N
2 density ratio is a height-integrated quantity within the orbit altitudes of the GNSS constellation and the GUVI satellite. In addition, global topside TECs derived from Swarm-A and Swarm-B were investigated. TEC signifies the integrated electrons within the altitudes from the orbit of the Swarm microsatellite to the orbit of the GNSS constellation.
Figure 9 provides an overview of daily topside TEC and O/N
2 during the storm. It is also worth noting that the TEC distribution differs slightly among the data from the two satellites, most likely due to the ~80 km difference in altitude. Both the Swarm-A and Swarm-B profiles conclude that the topside TECs over the Asian–Australian and American sectors were quiet before 12 UT on 25 August, with an average TEC of under 6 TECU.
Figure 9b shows that the TEC profile over the American sector suddenly increased, and the expanded profiles covered the eastern Pacific. This phenomenon was also validated by the observation of Swarm-B. During 00–12 UT on 26 August, the TEC over the Asia–Australia sector increased remarkably, with a maximum value exceeding 12 TECU. After, the enhanced TEC profiles gradually decreased and returned to a normal level. The profiles derived from Swarm-B agreed well with those of Swarm-A in that large-scale positive TEC irregularities appeared over the Asian–Australian and American sectors.
Figure 9i illustrates that the density ratio O/N
2 exhibited an inverse relationship with geographical latitude, with a range of 0.3 to 0.6 on quiet days. However, the O/N
2 ratio underwent a sudden change with the onset of a storm. Specifically, the O/N2 in low to middle latitudes increased significantly, reaching a maximum value of 0.8, while in polar regions, it decreased to 0.2. Moreover, the enhanced O/N
2 had a southward excursion in the American sector. As shown in
Figure 9j, the O/N
2 in North America was approximately 0.4, while in South America, it increased to 0.8. In the following two days, the storm-induced O/N
2 gradually decreased as the geomagnetic field returned to normal levels. The change in the O/N
2 density ratio was consistent with the TEC irregularities, indicating that it could be a crucial driver in generating plasma disturbances. O/N
2 has a strong positive correlation with plasma density and is a reliable indicator of neutral composition disturbances for analyzing ionospheric storms [
28]. The loss rate of ionospheric ion density is proportional to the molecular concentration. An increase in the mean molecular mass results in a decrease in the electron density, whereas a decrease in the molecular concentration leads to a positive disturbance.
At altitudes of 90 to 130 km, many electrons move westward under the influence of the dayside electric field. This westward electron flow generates a dayside eastward electric current through the equatorial dynamo effect, which is known as the equatorial electrojet (EEJ). During strong geomagnetic storms, the EEJ can be severely disrupted by disturbed electric fields from the magnetosphere. The signature of the EEJ can be estimated by taking the difference between the horizontal components measured by a pair of magnetometers located off the equator and at the equator. In this study, the horizontal components of the PHU, DLT, SJG, and HUA magnetometers located in the Asian–Australian and American sectors were used to investigate the storm-induced changes in the EEJ. The observations were obtained from the International Real-Time Magnetic Observatory Network (
https://intermagnet.org/index-eng.php, Accessed on 20 May 2021).
Figure 10a,b illustrate the changes in the EEJ signatures I, Accessed on 20 May 2021 the Asian–Australian and American sectors in response to the geomagnetic storm that occurred on 25 August 2018, as well as the corresponding variations in the horizontal components of the PHU, DLT, SJG, and HUA magnetometers. It was observed that the Dst index and magnetometers’ horizontal components experienced a sudden decrease during the storm. In
Figure 10a, the equatorial magnetometer GUA recorded a minimum value of approximately −270 nT during UT 06–08 on 26 August, and the differential component between GUA and KAK showed a negative perturbation with a minimum value of −100 nT. This was accompanied by a maximum TEC enhancement of 18 TECU observed at the PIMO station during the most severe moment. The changes in the EEJ, estimated by the difference between HUA and KOU, as shown in
Figure 10b, were consistent with those observed in the Asian–Australian sector. The EEJ signature showed two distinct perturbations in the afternoons of 26–27 August, with an amplitude of approximately 50 nT. The differential TEC over the UNSA station also corresponded well with the changes in the EEJ signature. During the recovery phase of the geomagnetic storm, two TEC enhancements were observed, with the maximum delta TEC reaching 14 TECU. These results suggest that changes in the EEJ may be a significant driver in triggering storm-effect TEC disturbances. However, it should be noted that the slight fluctuations in the EEJ cannot fully account for the strong TEC enhancements, indicating the need for further analysis of other drivers responsible for the significant ionospheric disturbances.
The equatorial plasma fountain effect plays a dominant role in generating the equatorial ionization anomaly (EIA) by driving equatorial plasmas upward to higher altitudes via the vertical E×B drift. The accumulated plasmas then diffuse down to higher latitudes along the geomagnetic field lines, resulting in two highly concentrated plasma crests on either side of the magnetic equator. A stronger E×B drift in the fountain process can lift more plasmas to higher altitudes, resulting in stronger and more poleward EIA crests generated by the plasma diffusion process. Therefore, changes in the vertical E×B drift may be another driver of equatorial and mid-latitude plasma perturbations. The latitudinal changes of the vertical E×B drifts along the −70° and 110° longitudes were simulated using the TIE-GCM, as shown in
Figure 10c,d.
Figure 10c shows that the E×B drifts along the meridian −70° began to increase from UT 14 on 25 August, with a slight enhancement of 3 to 5 m/s. However, the nighttime E×B drifts suddenly weakened, with a maximum decrement of −15 m/s at dawn on 26 August (local time). Subsequently, the daytime differential E×B drifts turned positive from UT 8 on 26 August, with a maximum increment of 15 m/s. On the forenoon of 27 August (DOY 239), the differential E×B had a hemispheric asymmetrical structure, with the differential E×B on DOY 239 enhancing by ~5 m/s in the Southern Hemisphere. In the Asian–Australian sector, the E×B drift enhanced from 20 UT on 25 August, with a magnitude of 5 to 10 m/s (local time). A slight positive E×B irregularity was also observed on the following day.
The results of this study not only confirm some previous papers but also conclude new findings. For instance, Astafyeva et al. [
29] utilized ground-based instruments, such as GNSS, SuperDARN, and ionosondes, as well as space-borne tools, including Swarm, CSES, and DMSP, to detect significant storm-time enhancements in the Asian and American sectors on a large scale. The storm-enhanced equatorial E×B drift was the main driver for the most significant 100–300% large-scale enhancement occurred over the Asian region. In the American sector, the complex plasma irregularities were caused by the action of multiple drivers. This study confirms the enhanced vertical E×B drift successfully, which coincide with the Astafyeva et al.’s study. However,
Figure 9 and
Figure 10 also indicate that the thermospheric density ratio of O/N
2 and equatorial electrojet were severely impacted by the strong geomagnetic storm. Specifically, the equatorial electrojet was activated by the disturbed electric field that penetrated from the magnetosphere, and the changes in vertical E×B drifts may be associated with PPEFs and DDEFs. Therefore, it is believed that the irregularities in the equatorial and mid-latitude ionosphere are the result of a combination of multiple physical–chemical processes. The enhanced density ratio O/N
2, vertical E×B drift, and equatorial electrojet played a decisive role in inducing positive irregularities. In the recovery phase, the hemispheric asymmetrical O/N
2 and E×B drift on 27 August may be responsible for the asymmetrical TEC observed over the American sector, as shown in
Figure 3.
3.4.2. Potential Drivers of Auroral Ionospheric Irregularities
The results from the GNSS ROT, sounder density profiles, and TIE-GCM simulations reveal significant auroral ionospheric irregularities induced by the storm. To further investigate the vertical structures of these irregularities, the TIE-GCM simulated plasma irregularities between altitudes of 96 and 400 km along the −70°, 0°, and 110° longitudes, as shown in
Figure 11a–c. In
Figure 6, the TIE-GCM simulation finds that the largest plasma irregularities, with a changing of 250%, are located in the 70°S–80°S latitude span at a layer of 150 km. Therefore, the geographical latitude of 77.5°S was selected, as shown in
Figure 11a–c. It is important to note that the vertical scale of each panel is similar to
Figure 7 but represents the amplitude of the plasma irregularities. The results show significant plasma enhancements in the topside and bottom-side of the Antarctic ionosphere along all three longitudes. During the main phase, the increment along the meridian −70° was the largest, with a value up to 6 × 10
5 el/cm
3, followed by the meridian 0°, while the weakest increment of 4 × 10
5 el/cm
3 occurred along the meridian 110°. It is noted that the difference between the −70° and 110° meridians is 180°, and the occurrence times (LT) of largest increments over two sectors are opposite. In addition, plasma fluctuations were also observed in the bottom-side ionosphere along the meridians −70° and 0°, except for the meridian 110°.
Was the disturbance caused by the geomagnetic storm limited only to the bottom-side ionosphere in western Antarctica? To investigate, we analyzed the temporal variations in the differential bottom-side plasma (150 km) as a function of geographical latitude along the meridians −70°, 0°, and 110°, as shown in
Figure 11d–f. The results indicate that the storm-induced plasma increments were significant in the bottom-side ionosphere across all three longitudinal sectors, but the geographic latitudes of the increments varied. In
Figure 11d, two plasma crests were located at latitudes 60°N and 80°S. These crests moved northward in the Eastern Hemisphere, with the plasma crests in the Asian–Australian sector shifting to latitudes 80°N and 60°S, with weaker values of 3 × 10
5 el/cm
3. The latitudinal motion of bottom-side plasma enhancements in different longitudinal sectors was influenced by the asymmetrical structure of the geomagnetic field.
Figure 11g–i depict the longitudinal structures of the differential plasma (150 km) as a function of the day of year (DOY) for three geographic latitudes: 77.5°S, 67.5°S, and 57.5°S. In
Figure 11g, positive plasma irregularities occurred in the −120°~0° longitudinal span on 26 August, with a value of 4 × 10
5 el/cm
3. Along the 67.5°S latitude, the plasma irregularities had a double-peak structure in the −180°~60° longitudinal span. Along the 57.5°S latitude, double-peak plasma increments were observed in the 0°~90° and 150°~−120° longitudinal spans. The intensity of the bottom-side plasma irregularities decreased in the middle geographical latitudes compared to the auroral plasma irregularities. The results suggest that the strong storm not only induced topside plasma fluctuations but also triggered positive bottom-side plasma irregularities near the auroral zone (~>50°N/S), which is consistent with the sounder profiles of the radar PQ052, as shown in
Figure 5a.
To explain the development of the auroral ionospheric irregularities, the TIE-GCM was used to simulate the Joule heating and the O
+ and O
2+ ion densities within the altitudes of 100 to 400 km. The neutral mass density decreased exponentially with height, and the Joule heating per unit mass was much larger at higher altitudes than at lower altitudes. Neglecting the bottom-side change, the ratio of changed Joule heating during storms compared to quiet background values was also investigated. The study area in the Arctic was located at 110°E, 67.5°N.
Figure 12a–d show the vertical changes in the Joule heating, the ratio of enhanced Joule heating, and the O
+ and O
2+ ion density over the Arctic. The results show that the Joule heating was enhanced from the main phase with a magnitude of 1 × 10
4 erg/g/s, reached a maximum in the recovery phase with a value of 3 × 10
4 erg/g/s, and gradually recovered to background values thereafter.
Similar to the variations in the topside Joule heating, the positive Joule heating disturbances also appeared below 200 km on DOY 237–238, with a slight value of 1 × 10
4 erg/g/s.
Figure 12b shows that the Joule heating in the main and recovery phases enhanced more than 200 times compared to the background values, with the maximum Joule heating enhancements located at altitudes of 100–150 km. The changes in the O
+ ion density, as shown in
Figure 12c, agree well with the Joule heating, with the O
+ ion density increasing from DOY 237 and growing stronger on DOY 238 above the 200 km layer, with the maximum value reaching 1.5 × 10
5/cm
3. An O
2+ increment generated from UT 14, DOY 237, grew to 1 × 10
5 el/cm
3 on DOY 238. Unlike the O
+ ions, the O
2+ increments mainly occurred below the 200 km layer, consistent with the behavior of the bottom-side enhanced Joule heating, as shown in
Figure 12b.
In the Antarctic study area, located at 70°W, 77.5°S, the changes in the Joule heating were found to agree with those over the Arctic. However, the positive Joule heating disturbances only appeared in the main phase, which is consistent with the variation in the storm-effect on Antarctic ionospheric plasma, as shown in
Figure 11a. The magnitude of the enhanced Joule heating over the Antarctic was several times larger than that over the Arctic.
Figure 12e shows that the topside and bottom-side Joule heating were enhanced by approximately 8 × 10
4 erg/g/s and 2 × 10
4 erg/g/s, respectively. However,
Figure 12f indicates that the changed ratio of the bottom-side Joule heating in the main phase was significantly larger than that in the topside, and the maximum enhanced ratio exceeded 1 × 10
4. Similar to
Figure 12c–d, in the main phase, the O
+ ion density over the Antarctic was enhanced by approximately 5 × 10
5/cm
3 above the 200 km layer, and the bottom-side O
2+ ion density was enhanced by approximately 2 × 10
5 el/cm
3.
On DOY 237–238, a positive Joule heating disturbance also appeared under 200 km with a slight value of 1 × 10
4 erg/g/s, similar to the variation of the topside Joule heating.
Figure 12b shows that in the main and recovery phases, the Joule heating enhancements were more than 200 times that of the background values, and the maximum enhancements were located at altitudes of 100–150 km. The changes in the O
+ ion density, as shown in
Figure 12c, agreed well with the Joule heating, with the O
+ ion density enhancing from DOY 237 and growing stronger on DOY 238 above the 200 km layer, reaching a maximum value of 1.5 × 10
5/cm
3. On DOY 238, an O
2+ increment was generated at UT 14, growing to 1 × 10
5 el/cm
3, and mainly occurring under the 200 km layer, consistent with the behaviors of the bottom-side enhanced Joule heating, as shown in
Figure 12b, which is different from the behavior of the O
+ ions.
Similar to the variations in topside Joule heating, a positive Joule heating disturbance also appeared below 200 km on DOY 237–238, with a slight value of 1 × 10
4 erg/g/s.
Figure 12b shows that, different from the absolute change in the differential Joule heating in
Figure 12a, the Joule heating in the main and recovery phases was enhanced by more than 200 times that of the background values, with the maximum enhancements located at altitudes of 100–150 km. The changes in the O
+ ion density, as shown in
Figure 12c, showed good agreement with the Joule heating, with the O
+ ion density increasing from DOY 237 and growing stronger on DOY 238 above the 200 km layer, with a maximum value of 1.5 × 10
5/cm
3. An O
2+ increment was generated from UT 14, DOY 237, and grew to 1 × 10
5 el/cm
3 on DOY 238. Different from the O
+ ions, the O
2+ increments mainly occurred under the 200 km layer, which is consistent with the behavior of the bottom-side enhanced Joule heating, as shown in
Figure 12b.
In the Antarctic, the study area was selected at 70°W, 77.5°S. The changes in the Joule heating over the Antarctic were found to be in good agreement with those over the Arctic. However, positive Joule heating disturbances only appeared in the main phase, which is consistent with the variation in the storm-effect Antarctic ionospheric plasma, as shown in
Figure 11a. The magnitude of the enhanced Joule heating over the Antarctic was several times greater than that over the Arctic.
Figure 12e illustrates that the topside and bottom-side Joule heating was enhanced by approximately 8 × 10
4 erg/g/s and 2 × 10
4 erg/g/s, respectively. However,
Figure 12f indicates that the changed ratio of the bottom-side Joule heating in the main phase was significantly larger than that in the topside, with the maximum enhanced ratio exceeding 1 × 10
4. Similar to
Figure 12c,d, in the main phase, the O
+ ion density over the Antarctic was enhanced by approximately 5 × 10
5/cm
3 above the 200 km layer, while the bottom-side O
2+ ion density was enhanced by approximately 2 × 10
5 el/cm
3.
The amplitudes of the Joule heating and the O
+ and O
2+ ion density over the Antarctic were considerably stronger than those over the Arctic, which is consistent with the magnitude of the polar ionospheric plasma disturbance, as depicted in
Figure 11. During a space weather event, sudden enhanced energy could ionize the primary neutral gases, O
2 and N
2, resulting in an increase in the ion density [
30]. Our simulations have confirmed the theory that enhanced Joule heating can accelerate the Polar ionospheric ionization process, and the enhanced O
+ and O
2+ ion densities are responsible for the topside and bottom-side plasma increments, respectively.
Geomagnetic storms not only create storm-enhanced densities (SEDs) in low–middle latitudes and tongues-of-ionization at the polar cap but also alter the global magnetic field and strengthen the ionospheric–magnetospheric current systems, as noted in [
31]. Since the ionosphere acts as a conductor, electric fields in the polar ionosphere can generate ionospheric currents similar to that of close field-aligned currents flowing in the ionospheric–magnetospheric system. This phenomenon can lead to Joule heating in the upper atmosphere and even influence the circulation of ionospheric plasma, thereby altering the polar ionospheric electron density structure. The ionospheric electric potential contour maps, estimated using the Super Dual Auroral Radar Network’s (SuperDARN) Assimilative Mapping procedure, accessible at
http://vt.superdarn.org/tiki-index.php?page=ASCIIData, Accessed on 20 May 2021, were utilized to examine the spatial–temporal variations in the polar convection patterns during the main phase of the geomagnetic storm.
Figure 13 shows the red and blue contours indicating the positive and negative ionospheric electric potentials, respectively. Typically, the potential pattern reaches a maximum near dawn and a minimum near dusk, and the difference between these values is referred to as the cross-polar voltage. The spatial–temporal evolution of the electric potentials over the Arctic from UT 20, 25 August, to UT 18, 26 August, is presented in
Figure 13a–f. At UT 20, 25 August, the negative electric potential with a minimum magnitude of −32 kv occurred in a −60°~120° longitudinal span, whereas a positive electric potential with a maximum value of 41 kv was distributed in a −90°~−180° longitudinal span. Starting from 26 August, the polar electric potential remarkably intensified with the decreasing Dst index, and the strongest electric potential was observed in UT 04–08, 26 August, coinciding with a peak AE index of 1500 to 2000 nT. For instance, in
Figure 13d, the negative potential decreased to −53 kv, while the positive potential increased to 36 kv, resulting in a cross-polar voltage of 89 kv. During the main phase’s end, the convection zone in the Arctic extended to approximately 50°N, consistent with the geographical latitude of radar PQ052, as shown in
Figure 2a, the station with the minimum latitude capable of detecting bottom-side ionospheric regularities. Similarly, the scale and scope of the electric potential over Antarctica significantly expanded and intensified. The negative potential decreased from −32 kv at UT 20, 25 August, to −61 kv at UT 04, 26 August, while the corresponding positive potential reduced from 40 kv to 33 kv. The maximum cross-polar voltage on UT 04, 26 August, reached 94 kv. The scale of the Antarctic storm-effect electric potential was stronger than that over the Arctic, although the scope was smaller.
During active space weather events, sudden energy and momentum are deposited in the high-latitude ionosphere and thermosphere, primarily in the form of particle precipitation and Joule heating. The incident particles gradually transfer energy to different ionospheric layers, ionizing more charged particles as the energy deposition increases. The accelerated ionization process enhances the ionospheric current flowing in the medium, while particle precipitation and Joule heating control the variations in the short-scale structures of the ionosphere–thermosphere, resulting in increased electrical conductivity and heating of the ionosphere–thermosphere system. TIE-GCM simulations of the 5 April 2010 geomagnetic storm concluded that additional particle precipitation not only significantly increases ionospheric conductivity but also causes remarkable Joule heating enhancements [
32]. The enhanced conductivity, electric field, or a combination of both could intensify the ionospheric electric currents. The current density is directly proportional to ionospheric conductivity, which is itself proportional to the plasma density [
33]. Thus, there is a close connection among magnetosphere energy deposition, particle precipitation, ionospheric current intensification, Joule heating, and SEDs generation.
Figure 12 and
Figure 13 reveal that the plasma density, Joule heating, and ionospheric electric potential affected by the storm were all significantly enhanced, further confirming the charged particle diffusion process reported in the previous literature. Therefore, it is believed that the storm-induced polar plasma irregularities are associated with the additional energy input from the particle precipitation, Joule heating, and ionospheric current intensification.