Temporal change of EIA asymmetry revealed by a beacon receiver network in Southeast Asia

Abstract

To reveal the temporal change of the equatorial ionization anomaly (EIA) asymmetry, a multipoint satellite-ground beacon experiment was conducted along the meridional plane of the Thailand–Indonesia sector. The observation includes one station near the magnetic equator and four stations at off-equator latitudes. This is the first EIA asymmetry study with high spatial resolution using GNU Radio Beacon Receiver (GRBR) observations in Southeast Asia. GRBR-total electron contents (TECs) from 97 polar-orbit satellite passes in March 2012 were analyzed in this study. Successive passes captured rapid evolution of EIA asymmetry, especially during geomagnetic disturbances. The penetrating electric fields that occur during geomagnetic disturbed days are not the cause of the asymmetry. Instead, high background TEC associated with an intense electric field empowers the neutral wind to produce severe asymmetry of the EIA. Such rapid evolution of EIA asymmetry was not seen during nighttime, when meridional wind mainly controlled the asymmetric structures. Additional data are necessary to identify the source of the variations, i.e., atmospheric waves. Precisely capturing the locations of the crests and the evolution of the asymmetry enhances understanding of the temporal change of EIA asymmetry at the local scale and leads to a future local modeling for TEC prediction in Southeast Asia.

Background

The equatorial ionization anomaly (EIA) is a significant combined ionospheric feature. During the daytime, the E × B plasma drift around the geomagnetic equator is upward due to the eastward electric field, whereas it reverses at night. The plasma is lifted upward during daytime and then it diffuses away from the magnetic equator along geomagnetic field lines due to the Earth’s gravitational and pressure gradient forces. The combination of electrodynamic drift (Martyn 1955) and diffusion (Mitra 1946) produces a fountain-like motion of the plasma known as the equatorial fountain effect. Consequently, ionization peaks, commonly known as EIA, form in both hemispheres. Namba and Maeda (1939) first discovered EIA before World War II (Rishbeth 2000). It was the major outcome of a study of the midday F2-layer critical frequency (foF2) from stations at different latitudes (Bailey 1948; Nishida 2009).

Interactions between transequatorial neutral winds and the strength of the equatorial fountain effect play important roles in the asymmetric development of the EIA crests. Aydogdu (1988) studied the thermospheric wind effect on EIA asymmetry for equinox seasons in African and West Asian regions and found that the NmF2 in the northern crest was greater than that in the southern crest during daytime. Maruyama (1996) proposed that it is possible to determine the direction of the meridional neutral wind by using the ionospheric height variations at the magnetic equator and at the magnetic conjugate points in the same magnetic meridian.

The first low-latitude ionospheric tomography network (LITN) covering 14.6° N to 31.3° N (or 3.3° N to 19.7° N geomagnetic latitude) along 120° E longitude was initiated in 1991 (Yeh et al. 1994). These researchers measured the total electron content (TEC) from satellite dual-band beacons and subsequently explored the unique science associated with the low-latitude ionosphere, e.g., EIA morphology and dynamics. Yeh et al. (2001) employed more than 350 ionospheric tomography images acquired in 1994. They found that the equatorial dynamo electric field is the main cause of the crest motion, and the correlation between the integrated strength of the electrojet and the strength of the equatorial fountain effect was high. They also showed that day-to-day variability of the EIA was controlled dominantly by the electrodynamics. These findings were deduced from observation in a low-solar-activity year from only the northern hemisphere. The LITN improved the study of the ionospheric structures at low latitudes on the Taiwan meridian.

Global-scale observations of ionospheric structures were conducted with six micro-satellites termed Formosa Satellite 3 and Constellation Observing System for Meteorology, Ionosphere, and Climate (FORMOSAT-3/COSMIC or F3/C for short) (Cheng et al. 2006). Lin et al. (2007a) investigated EIA crests during July–August 2006 by employing the F3/C GPS occultation experiment. The EIA peak altitude and density from F3/C revealed clear longitudinal variability that is probably due to differences in magnetic declination, E × B drift, and the neutral winds at different longitudes. However, clear seasonal asymmetry was not seen in the TEC map (integrated from a 100- to 500-km altitude), except in the 50°–150° E longitude regions. Although it provided valuable ionospheric information at the global scale, including over the oceans, the F3/C could not provide high spatial resolution of the ionosphere at the same longitude.

Xiong et al. (2013) studied EIA morphology, inter-hemisphere asymmetry, and solar activity levels based on 9 years of CHAMP and GRACE observations. From their observations, crests were shown, based on global averages, to always be farther from the dip equator in the summer hemisphere at both CHAMP and GRACE altitudes (~400 and ~480 km, respectively). The electron density and magnetic latitudes of EIA crests were almost symmetric about the dip equator during the equinox season. The study using CHAMP and GRACE satellites globally revealed EIA information from a different altitude. However, it took 130 days for CHAMP (Reigber et al. 2002) and 160.5 days for GRACE (Tapley et al. 2004) to cover all local times in the same longitudinal sector.

Pancheva and Lysenko (1988) first reported a signature of the quasi-2-day oscillation in the peak electron concentration of the F-region (foF2). The wave period is in the range of 34 to 68 h (Pancheva et al. 2004). Chen (1992) found that there was a 2-day oscillation in the EIA. Forbes et al. (1997) suggested a possible connection between the 2-day oscillation in foF2 and the 2-day wave in the mesosphere and lower-thermosphere (MLT) region. Pancheva et al. (2006) observed a 2-day variability in geomagnetic fields due to E-region dynamo currents as well as modulation of tidal fields at E-region height. They suggested that the 2-day wave in the MLT region and the tidal field modulation could be attributed to the observed 2-day variation in the ionosphere. Takahashi et al. (2012) suggested that the 2-day-wave-induced thermospheric transequatorial wind could generate latitudinal transport of F-region plasma and therefore could cause day-to-day variability of EIA.

Although a number of studies have been carried out on EIA during the past seven decades, these studies have some limitations. For example, even though it provided ionospheric information with high spatial resolution, the LITN observation was limited to over only the northern hemisphere. Moreover, satellite observations are capable of capturing a global-scale distribution, but they are not good at capturing detailed structures. In addition, the F3/C did not allow consecutive observation at the same longitude, and the observations from the CHAMP and GRACE satellites took long collection times (130 and 160.5 days, respectively) to cover all local times at the same longitude. Finally, there has been no significant study of the precise structures of the EIA asymmetry across the geomagnetic equator, e.g. the locations of its crests (Chen et al. 2008). In this paper, we examine precise latitudinal structures of EIA over the Thailand–Indonesia region from a satellite-ground beacon experiment. The observation features are as follows: 1) EIA is observable in a narrow longitudinal coverage (~100 km). 2) The latitudinal extent at the ionospheric pierce point (IPP) covers ± 25 geomagnetic latitude. 3) Each observation has a quick scan within 20 min on average. 4) The availability of the SouthEast Asia Low-latitude IOnospheric Network (SEALION) provides information regarding the ionosphere’s condition and meridional wind behavior.

Methods

The GNU Radio Beacon Receiver (GRBR) was first designed and developed by Yamamoto (2008) to receive beacon signals in VHF and UHF bands from low earth orbit (LEO) satellites. A multipoint beacon experiment along 100° E began collecting observations in March 2012. The ionosonde chain under the SEALION project (Uemoto et al. 2007), which is a joint project among several countries including Japan, Thailand, and Indonesia, is also along this same meridian. Both GRBR and the ionosonde networks now cover the geomagnetic northern and southern hemispheres, including the geomagnetic equator, in Southeast Asia. The locations of these networks are summarized in Table 1. Data from the GRBRs and ionosondes from March 2012 are mainly used in this paper.

Table 1 GRBR and ionosonde network information

Generally, EIA is expected to be symmetric during the March equinox. To see whether or not the GRBR can reveal the precise EIA structure, we precisely investigated the asymmetry of EIA. In Southeast Asia, a plasma bubble influence is also active in March. Small asymmetry is seen during plasma bubble appearance. We thus selected the March equinox for this study to highlight the sensitivity of the GRBR (Yamamoto 2008) and the robustness of the TEC estimation technique (Watthanasangmechai et al. 2014).

To capture EIA asymmetry, the signals from polar LEO satellites, for which orbital heights are 1000 km on average, are used to estimate GRBR-TEC. In this work, GRBR-TEC is the vertical TEC between the satellite altitude and ground, and the GRBR-TECs from 97 passes of polar-orbit satellites during March 2012 are analyzed. The units of GRBR-TEC are TEC units, 10¹⁶ el/m², and the satellite information is summarized in Table 2. The estimation technique of the absolute TEC from the latitudinal network of GRBR is explained in Watthanasangmechai et al. (2014). We assume that the EIA structure is stable over one satellite pass. Each satellite pass over the field of view of the five GRBRs covers about 20 min in time and 6000 km in distance at IPP altitude. IPP altitude is determined by the ionospheric height obtained from ionosondes. By connecting the TEC data from the meridional network of five GRBR sites, we can estimate the precise latitudinal TEC distribution around at most ± 25 geomagnetic latitude. This quick scan allows the EIA structures from both hemispheres to be captured clearly with high spatial resolution and without interpolation. Successive passes enable us to capture the temporal change of the EIA asymmetry.

Table 2 Information about the polar LEO satellites used in this work

The SEALION ionosonde is operated to continuously transmit radio waves from 2 to 30 MHz and receive echoes from the ionosphere in order to provide a bottom-side plasma density profile every 15 min (Kenpankho et al. 2011). The data from SEALION are used to estimate a proxy of the transequatorial wind by following the method proposed by Maruyama et al. (2008).

Results and discussion

EIA formation and daily variation

Figure 1 shows the daily variation of the 2-h binned GRBR-TEC of EIA for the whole period beginning from the 06–07 LT bin (top left panel) and ending at the 04–05 LT bin (bottom right panel). The K _p index for March 2012 is shown in Fig. 2. The data represented in Fig. 1 are during K _p <3−. The yellow line represents the latitudinal TEC profile. The numbers of latitudinal profiles used in each panel are as follows: two profiles for the 06–07 LT bin, one profile for the 08–09 LT bin, one profile for the 10–11 LT bin, two profiles for the 12–13 LT bin, five profiles for the 14–15 LT bin, eight profiles for the 16–17 LT bin, four profiles for the 18–19 LT bin, ten profiles for the 20–21 LT bin, five profiles for the 22–23 LT bin, four profiles for the 00–01 LT bin, five profiles for the 02–03 LT bin, and two profiles for the 04–05 LT bin. The horizontal axis on each plot presents geographic latitude. The geomagnetic equator for 100° E is located at 7.8° N geographic latitude. The inclination of the geomagnetic field is close to zero. The black error bars represent the standard deviation (σ) of the TEC at every degree. In the morning (06–07 LT), the plasma density is extremely low when the plasma is that surviving through the night. TEC is dominant from the sub-solar latitude (0° geographic latitude) to the south at post-sunrise (06–07 LT). TEC increases rapidly from 06 to 09 LT, which is an ordinal behavior after sunrise. TEC is dominant and has an arch shape around the dip equator. EIA crests are seen within the 10–11 LT bin, suggesting that EIA starts forming during this time. Whereas the northern crest is clearly dominant at the 10–11 LT and 12–13 LT bins, the southern crest becomes dominant after the 16–17 LT bin. The TECs at both hemispheres are quite symmetric during the 14–15 LT bin. The TECs at both hemispheres reach their peaks, as high as 70 TEC units (TECU) on average, during the 14–17 LT bin. The EIA then decays gradually from the 18–19 LT bin. The σ at the dip equator is larger at 18–19 and 20–21 LT compared with daytime. The EIA crests decay rapidly at 22–23 LT. The crest-to-crest width becomes smallest before the crests disappear at 00–01 LT. After midnight (02–05 LT), the TEC in both hemispheres and at the dip equator decreases further until it reaches the daily minimum just before sunrise (04–05 LT).

Fig. 1

Daily EIA asymmetry variations. Plotted in linear scale of the 2-h binned GRBR-TEC during K _p <3− in March 2012 along the meridian 100° E. The local time (LT) is indicated at the top of each panel. The green circle represents the average TEC value, the black bar represents the standard deviation (σ) at every 1°, and the yellow line represents the latitudinal TEC profile. The number of latitudinal profiles used in each panel is indicated in parentheses

Fig. 2

K _p index in March 2012. Plotted every 3 h with blue line. The red “o” mark indicates the K _p <3−

Figure 3 shows the EIA asymmetry parameters during K _p <3− in March 2012 with 2-h resolution. The TECs at the EIA crests and center are plotted in Fig. 3a. The TEC ratio, which is the ratio of the TEC of the northern crest to the TEC of the southern crest, is shown in Fig. 3b. The TEC ratio is larger than 1 when the TEC value of the northern crest is higher than that of the southern crest. The locations of the EIA crests and center are plotted in Fig. 3c. The crest peak latitude ratio, which is the ratio of the latitudinal width between the crest and the center in the northern hemisphere and that in the southern hemisphere, is shown in Fig. 3d. The crest peak latitude ratio is smaller than 1 when the crest peak latitude value in the northern hemisphere is smaller than that in the southern hemisphere. The latitudinal crest-to-crest width in degrees is illustrated in Fig. 3e. Figure 3f shows the expansion speed of the crests in degrees per hour. A positive speed means that the crests are expanding in the poleward direction, whereas a negative speed indicates equatorward motion of the crests. The TECs around the dip equator appear stable at about 50 TECU during 1030–1630 LT. This is probably due to a balance between ion production and the up-drift effect at the dip equator. As the fountain effect occurs repeatedly, the EIA crests at the off-equator continue developing until afternoon. The EIA crest-to-crest latitudinal width, which depends on the zonal electric field strength, continues to expand from 1030 LT to 1430 LT, as shown in Fig. 3e. The TEC is dominant in the northern hemisphere during 1000–1430 LT, when EIA expands rapidly.

Fig. 3

EIA asymmetry parameters. a The TEC values at the EIA crests and dip. b The TEC ratio of the northern crest and the southern crest. The TEC ratio is larger (smaller) than 1 when the TEC value of the northern crest is higher (lower) than that of the southern crest. c The EIA crest and dip locations. d The crest peak latitude ratio of the northern hemisphere to the southern hemisphere. The crest peak latitude ratio is smaller (larger) than 1 when the crest to latitude value in the northern hemisphere is smaller (larger) than that in the southern hemisphere. e The EIA crest-to-crest latitudinal width in degrees. f The crest expansion speed in degrees per hour. Positive (negative) sign indicates poleward (equatorward) motion of the EIA crests. All are plotted for March 2012 during K _p <3− with 2-h resolution

The transition of the daytime north–south asymmetry is at 15 LT, as shown in Fig. 3b. Before the transition hour, the expansion speed of the EIA decreases from 2° to 0°/h. The expansion speed is quite constant at around 0°/h after 15 LT. The average apex height of the EIA crests reaches 860 km. The TEC at the dip equator reaches its peak at 52.42 TECU at 1630 LT.

Near the dusk terminator, the equatorial F-layer rises due to the action of the dynamo electric fields; subsequently, it descends, and this descent is known as pre-reversal enhancement (PRE). The expansion speed increases suddenly around 1930 LT, probably due to the PRE. The crest-to-crest latitudinal width reaches its peak of 30° when the expansion speed increases from 0° to 1°/h. After 1930 LT, the expansion speed turns negative, which indicates a decrease of the EIA crest-to-crest width. The crest-to-crest latitude decreases rapidly from 30° in the 20–21 LT bin to 20° in the 22–23 LT bin with an expansion speed of −5°/h. The TEC continues to decrease after sunset, while the southern hemisphere continues to dominate until pre-midnight. The TEC at the dip equator increases slightly at pre-midnight, probably due to the midnight temperature maximum (MTM) (Colerico et al. 2006; Maruyama et al. 2008).

As a result of daily variation, the EIA is centered in the range of 7° N–9° N, except at 0030 LT (24.5 LT in Fig. 3) when the center appears at 4° N; the northern crest appears in the range of 15° N–21° N; and the southern crest appears in the range of 1° S–9° S. The average location of the EIA center is near the dip equator, and the average locations of northern crest and southern crest are 19.4° N and 6.5° S, respectively. By ignoring the crest and center information at 0030 LT, the crest is 2° farther from the dip equator in the southern hemisphere compared to the northern hemisphere. The asymmetry is also shown in Fig. 3d, as the average of the crest peak latitude ratio is 0.93. The EIA crests are fairly asymmetric after 1230 LT in this period average. The maximum TEC at both hemispheres and at the dip equator appears in the afternoon hours (1630 LT).

EIA day-to-day variability

The day-to-day variability of EIA can be interpreted using the σ in Fig. 1. Note that the σ at both ends of each TEC variation might be affected by the satellite pass coverage; therefore, these will be disregarded in the discussion. At post-sunrise (06–07 LT), the σ is small and quite constant because the background plasma is low and the dynamo electric field cannot greatly affect the day-to-day variability. In the morning (08–09 LT), the EIA day-to-day variability is controlled mainly by the day-to-day variability of electric field strength responsible for the background current (Stolle et al. 2008). In the 10–11 LT bin, the σ cannot be interpreted because only one overpass is available. In the 12–13 LT bin, the σ is small, suggesting that day-to-day variability of the EIA is small during its developing phase. In the 14–15 LT bin, the σ is still small at the dip equator, and it becomes slightly large around the EIA crests. This means that there is almost no day-to-day variability at the dip equator but that it exists around the crest latitudes. A large day-to-day variability around the crest latitudes is possibly contributed from the day-to-day variability of the crest locations and the EIA asymmetry during daytime. The σ at the dip equator is extremely low during daytime (10–15 LT), probably due to the balance between the ion production from solar radiation and the upward E × B drift. At 16–17 LT, the σ becomes large compared to the σ in the previous hours in all latitudes. At night, the σ is also large in all latitudes.

At night (18–01 LT), the EIA asymmetry is affected mainly by the temporal and day-to-day variability of the neutral wind (see below). In Southeast Asia, the March equinox is a favorable season for post-sunset plasma bubbles. In addition to the day-to-day variability itself, the plasma bubbles probably contribute to the large σ at night. In particular, relative variability of the TEC at the EIA crests is enhanced during 20–23 LT, when an equatorial spread F (ESF) and/or plasma bubbles are expected to occur. Shallow ripples considered as a plasma bubble characteristic on latitudinal GRBR-TEC (Watthanasangmechai et al. 2014) can be seen in the background yellow lines during 20–01 LT. In addition to the ripples in the GRBR-TEC profiles, the large variability of crest location contributes to the large σ at night. Note that the GRBR-TEC estimation error is about 3 % at night, whereas it reaches 10 % during ESF appearances (Watthanasangmechai et al. 2014).

Nighttime plasma bubbles are expected to occur around midnight (00–01 LT). Nevertheless, the March equinox during a high solar activity year is not a suitable condition for the event, and we did not find such an event from continuous observations using the Equatorial Atmosphere Radar. However, as shown in Fig. 1, while the averaged TEC is as large as 20 TECU, the relative σ becomes large for all latitudes. This enhanced relative σ might be caused by the day-to-day variability of the location of the EIA crests and the TEC value at decay time, and it may suggest a potential unstable situation. After the period, the relative σ decreases until predawn.

Figure 4 shows the TEC variations from March 5 to 9 (from top to bottom) from 11 LT to 01 LT. The time variation is separated into five bins: the 11–13 LT bin (noon), the 14–16 LT bin (afternoon), the 17–19 LT bin (evening), the 20–22 LT bin (night), and the 23–01 LT bin (midnight). These time bins are arranged from left to right, respectively, in Fig. 4. The daytime sector starts from noon to afternoon (11–16 LT), and the nighttime sector starts from evening to midnight (17–01 LT). The blank indicates lack of data. Well-developed EIA crests and the deep valley at the dip equator, which indicate intense electric fields, can be seen on March 8 at noon and on March 7 and 8 in the afternoon. High TECs at the EIA crests that are more than 80 TECU can be seen on March 8 at noon and on March 7–9 in the afternoon. The K _p index, considered as a proxy of geomagnetic activity, is shown in Fig. 2. The geomagnetic field condition was no-storm (K _p <4) on March 5–6, and there was a major storm (5 ≤ K _p ≤ 6) on March 7–8 and a severe storm (K _p >6) on March 9. During the morning hours (01–09 LT) of March 9, the Dst indices decreased to less than −100 nT, and the AE indices reached more than 1000 nT; this period is therefore classified as an intense geomagnetic storm. It is likely that the intensive EIA during the daytime on March 7–9 and the positive storm are correlated. These characteristics, however, were not noticeable in the nighttime sector. This suggests complex temporal variability of the background neutral wind. The TEC variations at night during March 7–9 were as normal as those on the other days.

Fig. 4

Day-to-day TEC variations on March 5–9. The TEC profiles during March 5–9, 2012, are arranged from top to bottom. The time variation is separated into five bins (from left to right): the 11–13 LT bin (noon), the 14–16 LT bin (afternoon), the 17–19 LT bin (evening), the 20–22 LT bin (night), and the 23–01 LT bin (midnight). The daytime sector starts from noon to afternoon (11–16 LT), and the nighttime sector starts from evening to midnight (17–01 LT). The blank indicates a lack of the TEC variation

During intense geomagnetic storms, the enhanced neutral wind and composition changes can penetrate all the way to the equatorial region. During the early stages of a magnetic storm, strong electric fields can penetrate to the low-latitude ionosphere (Gulyaeva et al. 2011) and create a region of strong plasma drift (Kelley 2009). Zonally eastward electric field penetration causes an uplift of the daytime equatorial ionosphere, removes ions from the recombination zone, and allows for further production of plasma in sunlight at higher altitudes. This uplift eventually leads to an enhanced fountain effect and an expansion of the EIA crests to higher latitudes. The electron density also increases as a result of storm dynamics (Cherniak et al. 2014). This is known as a positive ionospheric storm (Schunk and Nagy 2009; Simi et al. 2013).

Figure 5 shows the EIA asymmetry variation in March 2012 for quiet days (K _p <3−) and disturbed days (K _p ≥3−). The 14–15 LT bin is used to represent the daytime variation, and the 22–23 LT bin is used to represent the nighttime variation. The effect of the positive ionospheric storm (Liu et al. 2004; Klimenko et al. 2012) on the TEC variation during the daytime can be seen. As shown in Fig. 5a, c, the σ and average TEC in the daytime during the disturbed days are large compared to those during the quiet days, and similar behaviors have been described in other reports (Kutiev et al. 2006; Lin et al. 2007b; Kumar et al. 2014). The high σ around the EIA crests during the disturbed days could be caused by the highly variable penetration of the electric field that controls the strength of EIA. On the other hand, there is no significant difference in the σ during the quiet days and on the disturbed days at night. This suggests that high/low K _p affects the behavior of daytime but not nighttime EIA.

Fig. 5

EIA asymmetry variations in March 2012 along the 100° E meridian. Plotted for K _p <3− in a the 14–15 LT bin and b the 22–23 LT bin and for K _p ≥3− in c the 14–15 LT bin and d the 22–23 LT bin. The local time (LT) is indicated at the top of each panel. The green circle represents the average TEC value, the black bar represents the standard deviation (σ) at every 1°, and the yellow line represents the latitudinal TEC profile. The number of latitudinal profiles used in each panel is indicated in parentheses

As mentioned above, strong day-to-day variability around the EIA crests was most dominant from 20 to 23 LT. We thus investigated TEC in the 20–23 LT bin for nighttime day-to-day variation. The nighttime TEC variation from March 3 to 9 is illustrated in Fig. 6. Figure 6d is picked from the previous bin, where the pass on March 6 at 19:52 LT is available to fill the data gap. North–south asymmetry is clearly shown. The TEC was larger in the northern hemisphere on March 3, larger in the southern hemisphere on March 4, and larger in the northern hemisphere on March 5. Additionally, the TEC was larger in the northern hemisphere on March 6, larger in the southern hemisphere on March 7, and larger in the northern hemisphere on March 8. The TEC variations on March 3–5 and on March 6–8 have the same trend. The TEC on March 9 was larger in the northern hemisphere as expected. There are no data on March 10. The TEC on March 11 was expected to be dominant in the northern hemisphere; however, as shown in Fig. 7a, the opposite occurred. The consecutive observations during March 3–9, 2012, however, cover both quiet and disturbed days. The neutral wind inducing the EIA asymmetry during a geomagnetic storm could also take part in producing the TEC variation. In this work, we discussed the effect of the meridional wind on the nighttime EIA asymmetry.

Fig. 6

Day-to-day asymmetry variations during 20–23 LT. Plotted for a March 3, b March 4, c March 5, d March 6, e March 7, f March 8, and g March 9

Fig. 7

Examples of the TEC response on the meridional wind. a TEC variation on March 11, 2012, during 2120–2140 LT from COSMOS2463. b TEC variation on March 27, 2012, during 0016–0039 LT from COSMOS2454

EIA asymmetry

The neutral wind is typically considered the primary source of EIA interhemispheric asymmetry, whereas the zonal electric field is secondary. The magnetic inclinations at the bottom F-layer altitude are 25.413° and −19.073° at Chiang Mai and Kototabang, respectively. The thermospheric meridional wind is linked with the parallel upward motion along the geomagnetic field lines. The vertical component of the parallel upward velocity of plasma can be derived from the change of the F-layer bottom height. Because zonal electric fields also affect the upward drift of plasma, the effect from zonal electric fields is removed in order to use the F-layer bottom height for the transequatorial meridional wind estimation. The application to link meridional wind with an actual F-layer bottom height is limited to nighttime (Maruyama et al. 2008).

In this study, we estimate meridional neutral wind during nighttime in a magnetic meridional plane by taking the double difference of the ionospheric height change from a no-wind reference at conjugate points. The wind estimation method is explained in Maruyama et al. (2008). The double differential height is used as a proxy of the transequatorial meridional wind at about 250-km altitude. Figure 7a, b shows examples of the double differential heights as a proxy of the meridional wind on March 11 and March 26, 2012, which are the quiet days. The height change in kilometers is read as the wind velocity in meters per second, and the positive (negative) sign corresponds to the northward (southward) direction (Maruyama et al. 2009).

The response to the northward (southward) wind is expected because the northward (southward) wind induces decreasing TEC in the northern (southern) hemisphere. This is because the northward (southward) wind drags the plasma away from the southern (northern) hemisphere along the trough across the dip equator to the other hemisphere. The plasma is then forced to lower altitude by the neutral wind and by the Earth’s gravitational force. Since there is no solar input to generate photoionization at night, the plasma that reaches lower altitude thus decreases in content due to the high recombination rate.

Figure 8 shows the estimated meridional winds corresponding to the latitudinal TEC distributions of Fig. 6. Figure 6a shows the TEC response to the northward wind on March 11, 2012, during 2120–2140 LT from COSMOS2463. The TEC was 32.1 TECU in the northern hemisphere and was 44.0 TECU in the southern hemisphere. The TEC in the northern hemisphere was 11.9 TECU lower than the TEC in the southern hemisphere when the northward wind appeared. The northward wind speed was 36 m/s. The ratio of the crest-to-crest TEC deviation to the meridional wind speed was about 0.3 TECU per 1 m/s. The meridional wind had a northward direction from 1950 LT. This means that the delay between the TEC and the meridional wind was probably within one and a half hours. Figure 6b shows the TEC response to the southward wind on March 27, 2012, during 0016–0039 LT from COSMOS2454. The TEC was 29.3 TECU in the southern hemisphere and was 37.5 TECU in the northern hemisphere. The TEC in the southern hemisphere was 8.2 TECU lower than the TEC in the northern hemisphere when the southward wind appeared. The southward wind speed was 27.5 m/s. The ratio of the crest-to-crest TEC deviation to the meridional wind speed was 0.3 TECU per 1 m/s. The meridional wind had a southward direction from 2300 LT on March 26, 2012. The time delay would probably be within one and quarter hours.

Fig. 8

Double differential of height as a proxy of the meridional wind. The double differential heights were derived from the no-wind reference at the conjugate points on a March 11, 2012, and b March 26, 2012. The height change in kilometers is read as the wind velocity in meters per second, and the positive (negative) sign corresponds to the northward (southward) direction

Estimation of meridional wind, however, is successful for only limited cases. The availability of ionosonde data, which can be used to estimate meridional wind, is limited by the ionosonde operating condition and the absence of an intense ESF. Because an intense ESF appears as a spread trace in an ionogram, it is difficult to scale the ionogram precisely. To investigate the nighttime TEC response to the meridional wind in Southeast Asia, simultaneous availability of TEC data and ionosonde data is necessary. The number of satellite passes and the quality of the beacon signal, which is determined by its signal-to-noise-ratio (SNR) level, limit TEC availability. With these limitations, nine cases in March 2012 were found to be suitable for the investigation. The nine cases included four cases for the TEC response to the northward wind, three cases for the TEC response to the southward wind, one case during the quiet wind on the ESF active day, and one case for no response. By discarding the quiet wind case, seven of the remaining eight cases behaved as expected. The response delay was found to be within one and one quarter hours. The average ratio of the crest-to-crest TEC deviation to the meridional wind speed was 0.3 TECU per 1 m/s.

Daytime asymmetry is shown in Fig. 4. The locations of the crests varied drastically by time. Strong EIA strength indicates that intense electric fields were seen in the daytime of March 7–9, when geomagnetic activity was high. This is the result of the production of an enhanced penetrating electric field during the geomagnetic storm. Daytime EIA is quite symmetric on March 7, whereas it is asymmetric on March 8–9. The asymmetry on March 8–9 could be generated by the neutral wind. The changes in EIA asymmetry on March 8–9, 2012, from the 11–13 LT bin to the 14–16 LT bin indicate rapid evolution of the TEC asymmetry, suggesting a rapid change of the neutral wind during the times.

Considering Fig. 5, there is no strong north–south asymmetry in the 17–19 LT bin on March 6, 8, and 9, 2012. Drastic TEC changes from the 14–16 LT bin to the 17–19 LT bin are seen on both quiet (March 6, 2012) and disturbed days (March 8 and 9, 2012). The TEC values decrease about 50 % from the 14–16 LT bin to the 17–19 LT bin. On March 9, 2012, it is clear that the shapes of the north–south asymmetry and the TEC values in the 14–16 LT bin and in the 17–19 LT bin are drastically different. The severe asymmetry in the 14–16 LT bin becomes fairly asymmetric in the 17–19 LT bin with a reversed dominant crest. Note that the indicated geomagnetic activity on March 9, 2012, is severe storm. Such a rapid evolution of the north–south asymmetry on March 9, 2012, makes temporal prediction of the TEC around the EIA crests difficult.

As mentioned above, meridional wind is the main source of the EIA asymmetry. Neutral wind variation can be generated by forcing from the lower atmosphere, i.e., tides and waves, or forcing from the magnetosphere, i.e., the wind induced by the geomagnetic storm. However, we could not identify the source of the wind behavior in this research. Local operation of the GRBR network improves the dataset of temporal change of EIA asymmetry, especially when we compare it to in situ observation for which data are available at only certain altitudes.

Conclusions

The meridional GRBR network in Thailand–Indonesia was used to clarify the precise asymmetric structure of EIA. EIA asymmetry was successfully captured at the March equinox, when symmetric EIA was expected. From 1 month of data, 97 passes in March 2012, we successfully showed the day-to-day variability of EIA. Daily variation of the TEC values and the locations of the EIA crests along longitude ~100° E was precisely observed. North/south asymmetries were found at both daytime and nighttime. In daytime, the EIA strength variability was found to be proportional to the enhanced penetrating electric field associated with magnetospheric activity. Rapid evolution of EIA asymmetry during daytime was seen on geomagnetic disturbed days, although the penetrating electric field could not have formed the asymmetry. Instead, when neutral wind forces on the intense EIA exist during the penetrating electric field enhancement, it produces severe asymmetry of the EIA with respect to the high background TEC. Such rapid evolution of EIA asymmetry was not seen during nighttime. Meridional wind mainly controlled the asymmetric structures at nighttime. We found direct variation between the speed of the thermospheric northward (southward) wind and reduction of the TEC value in the northern (southern) hemisphere, which is evidence of an influence from the bottom-side ionosphere. Local operation of the GRBR network improves the dataset of precise temporal change of EIA asymmetry, which leads to an improvement of capability for ionospheric prediction (Watthanasangmechai et al. 2010, 2012) in Southeast Asia. We plan to expand the dataset to deepen knowledge of EIA behavior in future work.