Abstract

A moderate geomagnetic storm occurred on April 20, 2018. Using vertical total electron content (VTEC) maps provided by the Center for Orbit Determination in Europe, ionospheric responses to the geomagnetic storm could be identified in generated two-dimensional differential VTEC maps. During the day of the storm the enhancement of the equatorial ionization anomaly (EIA), product of the super-fountain effect was identified. A localized TEC enhancement (LTE) was also observed to the south of the EIA on April 20, 2018. It was also possible to visualize this LTE in a longitudinal section of the EIA as a third crest. The maximum increment of VTEC for the LTE was 204%. This LTE is quite unique because it happened during the expected solar cycle 24 and 25 minimum, and according to a previous study no LTE observation could be done for the last solar two-cycle minimum. The origin of the observed LTE is suggested to be partly product of the super-fountain effect. Finally, a localized TEC decrement (LTD) was observed towards the end of the day, April 20, 2018. Because this LTD consisted in the disappearance of the northern and southern crests of the EIA and this occurred during the recovery phase of the geomagnetic storm, it can be suggested that the LTD origin is due to the westward disturbance electric field. This mechanism was put forward by a past study that also analyzed the responses to a geomagnetic storm (the 2015 St. Patrick’s day storm), being one of the responses the inhibition of both crests of the EIA.

1. Introduction

Geomagnetic storms can be very disruptive. For instance, they can induce currents in main electricity grids, can adversely affect global positioning system by considerably reducing its accuracy and also can cause the loss of high-frequency radio communications (e.g., [1, 2]). The ionosphere, which is a layer of the atmosphere located between 60 and 1000 km in altitude, responds to geomagnetic storms in the form of ionospheric storms. Ionospheric storms can be classified as positive or negative (e.g., [35]), i.e., when there is an increment or decrement in the total electron density (TEC), respectively.

At low latitudes, a noticeable feature in the layer of the ionosphere is the equatorial ionization anomaly (EIA; [6]). This anomaly consists of two regions of enhanced plasma located at ±15° of the magnetic equator. This phenomenon is mainly due to eastward zonal electric field which in combination with the Earth’s magnetic field at equatorial latitudes produces an upward x drift of the plasma. Once high altitudes are reached, the plasma diffuses along the geomagnetic field lines due to pressure gradient forces in conjunction with the action of gravity [7]. This whole process is also known as the equatorial fountain effect. On the other hand when a geomagnetic storm sets in, it has been previously observed in several studies (e.g., [812]) that the EIA is also considerably enhanced, which is due to the penetration of the interplanetary electric field from high to low latitudes and which lasts for several hours [13]. As a consequence of the electric field penetration, there is an upward displacement of plasma to higher altitudes and also an expansion of the EIA in the direction of the poles (e.g., [12]). In this case, this phenomenon is known as the super-fountain effect.

Besides the expected diurnal response of the EIA to geomagnetic storms in the low and midlatitude ionosphere, there are still not well-defined physical mechanisms that can account for the anomalous increase or decrease of TEC outside the extent of the EIA. Very recently, Edemskiy et al. [14] have observed that the ionospheric response to the August 15, 2015, geomagnetic storm presented not only the regular increment of TEC in the EIA, but also two regions of localized TEC enhancements (LTEs or localized positive ionospheric storms) at central-middle (~35°S) and lower-high latitudes (~62°S) were identified. The aforementioned study suggested that the excursion of the vertical interplanetary magnetic field () was necessary but not a definite requirement for the appearance of the observed LTEs.

In this paper, publicly available vertical TEC (VTEC) data from the Center for Orbit Determination in Europe (CODE) were used to investigate the ionospheric responses to the April 20, 2018, geomagnetic storm that happened during the expected time of the solar cycle 24 and 25 minimum [15].

2. Data and Method

2.1. Ionospheric Data

In order to detect ionospheric anomalies due to the geomagnetic storm of April 20, 2018, global maps of VTEC provided by CODE were used. As indicated by previous studies [16, 17], global ionospheric maps (GIMs) are good indicators of space weather events.

GIMs are routinely produced every hour and they come in daily IONosphere Map EXchange files (IONEX; [18]). CODE makes IONEX files publicly available via their ftp site (ftp://ftp.aiub.unibe.ch/CODE/). GIMs provide VTEC in a global geographical grid with a resolution of 5° x 2.5° (longitude and latitude, respectively). The total number of cells in latitude is 71 and in longitude is 73. The VTEC values are given in TEC Units (TECU), where 1 TECU = 1016 electrons/m2. The accuracy of the GIMs is between 2 to 9 TECU. Due to the bit of structural complexity of the IONEX files, a program using the library NumPy from Python was written. NumPy has the ability to work efficiently with N-dimensional arrays [19]. The result after this preprocessing stage was to obtain 3D cubes of VTEC per each IONEX file.

Ionospheric disturbances can be properly identified from GIMs assuming that the VTEC within a certain range of days follows a Gaussian distribution (e.g., [2023]). Under this consideration, the upper bound () and lower bounds () can be defined as + 2 and - 2, respectively (where is the mean and the standard deviation). In this way, a 3D cube of differential VTEC (VTEC) at a confidence level of 95% was constructed based on the 3D VTEC cubes. Positive and negative ionospheric anomalies are identified in the 3D VTEC cube when

2.2. Geomagnetic Conditions

In Figure 1, the variability of two important geomagnetic indices, Dst and Kp, during the month of April 2018 can be observed. While the Kp index indicates the severity of a storm, the Dst index is a proxy for the strength of westward terrestrial ring current [24]. As seen in Figure 1, the Kp index for April 20, 2018, between 06:00 and 12:00 has a value of Kp = 6, which according to the National Oceanic and Atmospheric Administration (NOAA) space weather scales points to a moderate G2 storm. On April 19, 2018, at 22:00 UT the Dst index rises to a value of 11 nT, which marks the storm sudden commencement (SSC), and it also indicates the start of the initial phase of the storm [25, 26]. On April 20, 2018, at 05:00 UT the Dst decreased drastically to a negative value of -19 nT. The Dst index reached its minimum of -66 nT at 09:00 UT, time in which the main phase finished. Kp and Dst data were downloaded from the German Research Center for Geosciences (https://www.gfz-potsdam.de/en/kp-index/) and the World Data Center for Geomagnetism in Kyoto (http://wdc.kugi.kyoto-u.ac.jp/wdc/Sec3.html), respectively.

In addition to the geomagnetic indices, a very crucial parameter that drives the development of a geomagnetic storm is the vertical interplanetary magnetic field (). Gonzalez and Tsurutani [27] have indicated that, for strong storms, the observed southward is lower than -10 nT and it remains below this threshold for at least three hours. On the other hand, Gonzalez et al. [25] pointed that, for moderate storms (-100 nT < Dst ≤ -50 nT), ≤ -5 nT for at least 2 hours. It is the latter what is observed during the geomagnetic storm of April 20, 2018. In the bottom plot of Figure 1, it can be observed that -5 nT, most of the time ranging within April 20, 2018, at 01:00 until 16:00 UT. Data for in geomagnetic solar magnetospheric coordinate system was retrieved as hourly averages from the OMNI database (https://omniweb.gsfc.nasa.gov/form/dx1.html).

3. Results and Discussion

In order to better appreciate how maps of VTEC improve the detection of ionospheric anomalies due to geomagnetic storms, Figure 2 shows a GIM (left), a map containing the mean VTEC (center), and the VTEC map (right) for April 20, 2018, at 10:00 UT. It can be observed in the VTEC map that the two enhanced plasma regions that form the EIA are displaced towards the poles (30°N and ~25°S), confirming what it was noticed for other geomagnetic storms in past studies (e.g., [812, 28, 29]). Furthermore in the VTEC map of Figure 2, a well-defined LTE in the southern hemisphere at higher-middle latitudes (~44°S) can be clearly identified. The GIM is also showing this LTE; however, there is no doubt when looking at the VTEC map that the LTE is properly displayed showing clearly its complete spatial distribution.

3.1. Localized TEC Enhancement

After carefully inspecting the generated VTEC cube, as previously seen in Figure 2, a LTE was observed to appear on April 20, 2018, at ~06:00 UT (Figure 3). The LTE appeared about 8 hours after the SSC to the south of the EIA (~44°S) right over the Indian Ocean. At later times the LTE started to drift westward (same as the direction of the Sun). Edemskiy et al. [14] have also observed such displacement for the two LTEs they detected during the August 15, 2015, geomagnetic storm. At ~10:00 UT the LTE reached its maximum intensity with a VTEC of 4.07 TECU. At this hour the EIA reached also its peak VTEC intensity of 4.33 TECU in the northern plasma region. As seen in Figure 3, the enhancement of the EIA lasted only for about 6 hours (08:00 until 14:00 UT), whereas the LTE is almost visible until April 21, 2018, at 00:00 UT. At this time the LTE is located over the south Pacific to finally fade away. It is worth mentioning that approximately two hours after (16:00 UT) the enhanced EIA disappears, returns to its regular variability. A day before (April 19, 2018) and after (April 21, 2018) the occurrence of the geomagnetic storm, the VTEC maps for these days revealed that there was no EIA enhancement or appearance of LTEs.

In order to appreciate the shape of the EIA and the LTE, a meridional profile along the 58.5°E at 10:00 UT (14:00 local time for this longitudinal coordinate) is shown in Figure 4. At this time and longitudinal coordinate, which is pretty much the same for the capital city of Oman, Muscat, the LTE reached, as previously mentioned, its highest intensity. The EIA shape is also shown at the same UT for 2 days before and after the geomagnetic storm. From Figure 4 it can be noticed that, for April 18, 19, 21, and 22, 2018, at 10:00 UT, the EIA shape has a well-defined double-crest with a trough form [6], which is expected during daytime because of the equatorial fountain effect. During the day of the storm at 10:00 UT, Figure 4 also shows that the EIA is enhanced by showing a considerably increase of the VTEC in the northern hemisphere of ~133% and at the same time a poleward displacement of the crests, this being product of the super-fountain effect (e.g., [812, 29]). On the other hand, an unexpected response to the storm in the EIA shape is the additional presence of a third crest (~44°S), which corresponds to the LTE observed in the southern hemisphere in Figure 3. However, it can be seen that, during the days before the storm (April 18 and 19, 2018) and after (April 21 and 22, 2018), the meridional profiles show as well a slight increment of their TEC in the same position of the observed third crest, for the LTE the peak increment in VTEC is of ~204%. Here it could be suggested that the super-fountain effect, which is caused due to southern flip of , with some other physical mechanism may be playing an important role in the generation of this LTE.

Edemskiy et al. [14] have not only detected the two LTEs observed during the geomagnetic storm of August 15, 2015, but they also discovered further 26 LTEs for the period of 2002-2015. However, they indicated that, during the time of the solar activity minimum of 2006-2009, there was no detection of LTEs. Thus in the present study, the detected LTE during the April 20, 2018, storms is the first identification of such type of event (LTE) during a solar cycle minimum. According to Hathaway and Upton [15], between 2016 and 2019 the solar cycle 24 and 25 minimum is expected.

3.2. Localized TEC Decrement

In Figure 3 at 22:00 UT on April 20, 2018, and at 00:00 UT on April 21, 2018, a negative storm or a localized TEC decrement (LTD) could be already observed. However, for completeness, the whole event is presented in Figure 5. This LTD appears at 22:00 UT on April 20, 2018, right over the southern half of Mexico. At 00:00 UT on April 21, 2018, the LTD moved to the west and reached its minimum peak intensity with a VTEC of -3.21 TECU. This negative storm was quite transitory; by 04:00 UT it completely disappeared. On the other hand the GIMs for this period of hours, also presented in Figure 5, show apparently that the northern plasma region of the EIA disappears, although from the GIMs alone this is not completely clear.

In order to have a better view of what is occurring with the EIA, meridional profiles for 100°W at 00:00 UT between April 19 and 23, 2018, are presented in Figure 6. The 100°W longitudinal coordinate is about the same as the one for Mexico City (capital of Mexico), where the local time is 19:00 hours. On April 19, 20, 22, and 13, 2018, at 00:00 UT (19:00 local time) the EIA shape shows a double-crest with a trough, which is still expected after sunset. However, what was unforeseen, as to what is seen in Figure 5, was the disappearance of both crests on April 21, 2018, at 00:00 UT, time when the recovery phase is already under way. The disappearance of the crests corresponds to the negative ionospheric storm observed in the VTEC maps of Figure 5. In Figure 6 it can be seen that the peak of the meridional profile for April 21, 2018, at 00:00 UT falls nearly in the location of the magnetic equator. A similar situation, where the northern and southern crests faded away, happened during the 2015 St. Patrick’s Day storm on March 18, 2015, at the beginning of the recovery phase [30]. In the aforementioned study it was indicated that the westward disturbance electric field is the most likely mechanism that could inhibit the formation of both crests. Thus it can be suggested that the same process is causing the disappearance of both crests observed for the April 20, 2018 storm during the recovery phase.

4. Conclusions

Analysis of differential VTEC content maps during the G2 (moderate) geomagnetic storm that occurred on April 20, 2018, has revealed different ionospheric responses. The expected response to the super-fountain effect was observed for a period of 6 hours with a maximum increment of VTEC of 133% in the northern crest of the equatorial ionization anomaly.

A localized TEC enhancement has also been observed during the storm, located to the south of the EIA. The maximum increment of VTEC for the LTE is of 204%. This LTE appeared two hours before the EIA was enhanced and also lasted until nearly the end of the day. The LTE shows itself as a third crest in the longitudinal profile of the EIA. It is suggested that the super-fountain effect may be playing an important role in the origin of the LTE. Moreover, this LTE is quite unique because, according to a previous study which detected 26 LTEs during 2002-2015, no LTEs have been observed during the solar two-cycle minimum of 2006-2009, and the April 20, 2018, storm occurred during the solar cycle 24 and 25 expected minimum.

Finally, a localized TEC decrement (or negative ionospheric storm) was also identified. This LTD revealed itself in the longitudinal profile of the EIA as the depletion of both its northern and southern crest. The LTD was observed to occur during the recovery phase of the storm. This similar behavior was observed for the 2015 St. Patrick’s Day geomagnetic storm. The authors of that study attributed the inhibition of both crests to the westward disturbance dynamo electric field. Hence, the origin of the LTD during the April 20, 2018, storm is considered to be due to the same mechanism. However, more observations of such types of LTDs are needed in order to understand better their origin.

Data Availability

The ionospheric data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The author declares that he has no conflicts of interest.

Acknowledgments

The author would like to acknowledge the Center for Orbit Determination in Europe for making IONEX files available to the public.