Skip to main content

Day-to-day variation of pre-reversal enhancement in the equatorial ionosphere based on GAIA model simulations

Abstract

Using a whole atmosphere–ionosphere coupled model GAIA (ground-to-topside model of atmosphere and ionosphere), we have investigated which parameters mainly control day-to-day variation of vertical plasma drift at the evening terminator over magnetic equator, so-called pre-reversal enhancement (PRE). Day-to-day variations of the peak PRE are compared with those of electron density, eastward current density and eastward neutral wind in the E- and F-region over Chumphon (10.7° N, 99.4° E; 0.86° N magnetic latitude), Thailand during equinoctial months in 2011–2013. Eastward neutral wind in the F-region shows positive correlation with peak PRE, indicating that the F-region winds control the peak PRE through the mechanisms of the F-region dynamo (including E- and F-region coupling processes). Daytime eastward electric current at an altitude of 110 km, corresponding to equatorial electro jet (EEJ), is also positively correlated with the peak PRE. Correlation between the EEJ and PRE is the largest at 1700 LT, approximately 1.5 h prior to the peak PRE.

Introduction

Equatorial plasma bubble (EPB) is characterized by localized plasma density depletion in the ionosphere, and are well-known due to their unique nature and adverse effects on communication and navigation systems. The Rayleigh–Taylor instability (RTI) is known to be the mechanism of EPB generation (Sultan 1996). Pre-reversal enhancement (PRE) is an enhancement of the vertical E × B drift due to the eastward electric field at the evening terminator at the magnetic equator. The eastward electric field strengthens before its reversal to westward direction post-sunset. PRE is believed to be one of the main controlling factors for the generation of EPB (Abdu 2001). There have been many studies in the past to show the connection of the PRE on the generation of EPB using ground-based and satellite observations (Abadi et al. 2015; Dabas et al. 2003; Huang and Hairston 2015; Tulasi Ram et al. 2006). However, mechanisms controlling day-to-day variation of the EPB occurrence are still unknown although it is important for the ionospheric community. Abadi et al. (2015) established that the magnitude of PRE governs the latitudinal extension of EPB along with the peak value of virtual height (h′F) during the developmental phase of EPB. Previous studies revealed that the daytime E × B drift and Equatorial Electro Jet (EEJ) have good positive correlation with the EPB generation (Dabas et al. 2003; Tulasi Ram et al. 2006).

Abdu et al. (1983) found that the onset of Equatorial Spread-F (ESF) is directly correlated with PRE peak during the southern summer solstice and equinoctial months in the equatorial ionosphere over Fortaleza. They found strong dependence of PRE and ESF intensity with solar flux index (F10.7) during the equinoctial months of solar maximum (mainly for strong ESF cases). Recently, Huang (2018) reported that PRE and ESF occurrence are linearly correlated. They found that the probability of occurrence of ESF is very less when peak PRE is negative, increases largely when peak PRE is 0–40 ms−1 and becomes 80% or higher above the threshold value of 40 ms−1. These studies reveal that plasma bubble occurrence largely depends on magnitude of PRE although various factors affect it (Abdu 2001). Uemoto et al. (2010) carried out the ionosonde and magnetometer observations over Chumphon and Phuket, respectively, during November 2007 to October 2008 to study the link among PRE, EEJ and ESF onset. They found that the increase in the virtual height (h’F) of bottomside F-layer and EEJ strength are interconnected with the generation of ESF.

Fesen et al. (2000) simulated PRE using National Center for Atmospheric Research thermosphere/ionosphere/electrodynamic general circulation model (TIEGCM) and reported that electron density in the E-region is one of the most vital parameter to produce PRE. The E-region electron densities ≥ 104 cm−3 inhibit PRE development by ‘short-circuiting’ the F-region dynamo. Based on NCAR Whole Atmosphere Community Climate Model with thermosphere/ionosphere extension (WACCM‐X) simulations, Liu (2020) reported that day-to-day variation of PRE is affected by the large-scale day-to-day variability of the lower atmosphere, through E-region in the summer hemisphere, at the mid-latitudes. Recently, Shinagawa et al. (2018) have shown that daily and seasonal variations of the RTI growth rate simulated by the Ground-to-topside model of Atmosphere–ionosphere for Aeronomy (GAIA) are consistent with those of EPB observed by GPS scintillation and VHF radar. However, mechanisms controlling the day-to-day variation of PRE have not been investigated well. Using GAIA, in the present study, an attempt is made to investigate mechanisms controlling the day-to-day variation of PRE. This study could contribute to a goal of forecasting EPB before its actual incidence.

Data and instrumentation

In the present study, a whole atmosphere–ionosphere coupled model, GAIA consisting of three models combining an ionosphere model, a neutral atmosphere model, and an ionospheric electrodynamic model (Jin et al. 2011; Shinagawa et al. 2018), is used to investigate day-to-day variation of plasma density, neutral winds, and electric current in the ionosphere over magnetic equator at a location of (10° N, 100° E), which is close to Chumphon (10.7° N, 99.4° E; 0.86° N magnetic latitude) in Thailand. The longitudinal and latitudinal grid spacing in the ionospheric model is 2.5 and 2.0 degrees, respectively. The vertical spacing is 10 km below an altitude of 600 km. Temporal resolution of the output data used in this study is 30 min. The present GAIA simulation does not include the geomagnetically disturbed condition.

The vertical component of plasma drift due to the zonal component of the electric field at the magnetic equator, corresponding to PRE, is derived from the electric field and magnetic field data (around 1600–0000 LT) (local time) during the equinoctial months March, April, September and October of 2011–2013. The period 2011–2013 is chosen as it represents the high solar activity period when the occurrence rate of equatorial plasma bubble will be high. It is to be noted that day number 60–90 (61–91) corresponds to March, 91–120 (92–121) corresponds to April, 244–273 (245–274) corresponds to September, 274–304 (275–305) corresponds to October for the year 2011 and 2013 (2012). Figure 1 shows the day-to-day variations of upward E × B drift velocity with respect to local time at an altitude of 300 km at magnetic equator (10° N, 100° E) near Chumphon in Thailand during (a) spring equinox (March–April) and (b) autumn equinox (September–October) of 2011. It is clearly visible that upward drift velocity enhances during 1800–1900 LT and reaches a peak around 1830 LT. This enhancement of the upward drift velocity corresponds to PRE. In this study, a peak value of PRE between 1830 and 1900 LT is mentioned as peak PRE during both the autumn and spring months irrespective of the sunset time. In the figure, day-to-day variation can be seen in magnitude of the peak PRE. The peak PRE exceeds 30 ms−1 during a period of day-of-year between 62 and 72, but were intermittently smaller than 10 ms−1 during a period of day-of-year between 74 and 81.

Fig. 1
figure 1

Day-to-day variation of upward E × B drift velocity with respect to local time at an altitude of 300 km during a spring equinox (March–April) and b autumn equinox (September–October) of 2011 over (10° N, 100° E; GAIA co-ordinates) near Chumphon in Thailand

Results

Figure 2 shows the day-to-day variation of the peak PRE (magenta line) with (a, b) electron density, (c, d) eastward current density, and (e, f) eastward neutral wind at the time when the PRE reached a peak (around 1830–1900 LT) during (top) March–April and (bottom) September–October of the year 2012 over (10° N, 100° E; GAIA co-ordinates) near Chumphon in Thailand. The horizontal axis depicts the day number of the year. The day number 61–91 corresponds to March and day number 92–121 corresponds to April. The electron density peak (F2 peak) exists at an altitude between 350 and 550 km, and varies day by day. The F2 peak is elevated to higher altitudes by PRE around evening terminator. The electron density in the E-region is quite smaller than in the F-region because the E-region plasma disappears due to rapid recombination soon after the solar extreme ultraviolet (EUV) radiation ceases. The rapid recombination of the plasma also makes a steep gradient of the electron density at the bottomside of the F-layer. Day-to-day variation in the altitude of bottomside F-layer is also seen. From comparison with the peak PRE, it is found that the F2 peak and bottomside altitudes ascend and the electron density at the F2 peak altitude decreases when the peak PRE is intense (Fig. 2a, b). In Fig. 2c, d, eastward electric current is shown. The eastward current density is enhanced at an altitude of 110–120 km. This intense current density corresponds to EEJ. The intensity of EEJ varies day by day, and its day-to-day variation ranges between 0.05 and 0.12 Am−2 approximately. Intensification in the eastward current density coincides with enhancement of the peak in PRE, indicating positive correlation between peak PRE and current density on most of the days except for day numbers 68, 275, 284, 294 and 301. As shown in Fig. 2e, f, the neutral wind is eastward in the F-region and westward in the altitude range between 120 and 200 km. Day-to-day variation of the eastward wind in the F-region ranges from 50 to 100 ms−1, and displays distinct positive correlation with peak PRE. When the PRE is weak, the neutral wind weakens and correlation between the peak PRE and E-region neutral wind is not discernible. On few days, the positive correlation is not clearly seen (such as day numbers 69, 83, 88, 115, 248, 254, 266–271 and 281) which may lead to weak positive correlation.

Fig. 2
figure 2

Day-to-day variations of the peak PRE (magenta line) and a, b electron density, c, d eastward current density, and e, f eastward neutral wind at the time when the PRE reached a peak during (top) March–April and (bottom) September–October in 2012 over (10° N, 100° E; GAIA co-ordinates) near Chumphon in Thailand

We have investigated quantitatively correlation of day-to-day variation of the PRE with that of the electron density, eastward current density and eastward neutral wind at the time when the peak PRE occurs. Figure 3 displays the Pearson correlation coefficient of the peak PRE with (a) electron density, (b) eastward current density, and (c) eastward neutral wind during the equinoctial months of March, April, September and October at the altitudes of 100–600 km for the years 2011 (red line), 2012 (blue line) and 2013 (green line). In Fig. 3a, it is observed that the electron density in the E-region (at altitudes of 100–150 km) shows negative correlation with the peak PRE, and the correlation coefficients are approximately − 0.4 for 2011, − 0.2 to − 0.6 for 2012 and − 0.3 to − 0.4 for 2013. The electron density in the F-region (at altitudes of 300–450 km) also show negative correlation with the peak PRE with maximum negative cross-correlation coefficients of approximately − 0.5 for 2011 and 2013 and − 0.6 for 2012. The eastward current density correlation coefficient ranges between − 0.3 and 0.4 approximately in the E-region (100–150 km) for all the years (2011–2013) and positive correlation (0 to 0.6) above 250 km altitude (Fig. 3b). The cross-correlation coefficient shows two positive peaks. One is seen around 150–160 km altitudes, and the correlation coefficients are approximately 0.6 during 2011 and 2012, and 0.2 in 2013. The other is around 350–450 km altitudes, and the cross-correlation coefficients are approximately 0.6 for 2011 and 2012, and 0.4 for 2013. In the F-region, the Pederson current flows in the same direction as the electric field (leading to the positive correlation). Thus, increase in eastward electric field intensifies the Pederson current. The day-to-day variation of current density may be primarily controlled by the electric field. However, the less correlation coefficient values in both E- and F-region in 2013 is not known. As shown in Fig. 3c, the cross-correlation coefficients between the peak PRE and eastward neutral wind are between − 0.2 and 0.2 below 150 km altitude for all the years (2011–2013). Above 300 km altitude, the cross-correlation coefficient is approximately 0.4 and does not show altitude dependence (nearly constant). This could be because the neutral wind does not change largely with altitude due to high viscosity in the F-region (Yerg 1955).

Fig. 3
figure 3

Correlation between the day-to-day variation of the peak PRE and electron density, eastward current density and eastward neutral wind during 2011 (red line), 2012 (blue line) and 2013 (green line) over (10° N, 100° E; GAIA co-ordinates) near Chumphon in Thailand

Figure 4 portrays day-to-day variations of peak PRE with (a, b) electron density, (c, d) eastward current density, and (e, f) eastward neutral wind at an altitudes at (top) 1700 LT and (bottom) 1500 LT in an altitude range of 100–600 km (left Y-axis; contour plot) over (10° N, 100° E) near Chumphon during the equinoctial months of March and April 2012. The electron density in the F-region is lower at 1500 LT than 1700 LT. The negative correlation between the peak PRE and electron density, which can be seen at the time of PRE, is not seen at 1500 LT. Regarding the day-to-day variation of eastward current density around 110 km, that is EEJ, intensity of EEJ enhances intermittently (Fig. 4c, d). Enhancement of EEJ intensity coincide well with increase of peak PRE. It is noted that the magnitude of the current density for EEJ increases up to 2.0 Am−2 at 1500 LT and up to 0.7 Am−2 at 1700 LT (Fig. 4c). The current density of EEJ decreases with time because the electron density in the E-region decreases rapidly at sunset. Figure 4e, f displays the day-to-day variation of eastward neutral wind at altitudes of 100–600 km and peak PRE over (10° N, 100° E). Correlation between the neutral wind and peak PRE at 1500 and 1700 LT is not discernible although the correlation is significant at the time of the peak PRE.

Fig. 4
figure 4

Day-to-day variation of electron density ranging from 100 to 600 km (contour, in m−3), eastward current density ranging from 100 to 150 km (contour, in A m−2) and eastward neutral wind ranging from 100 to 600 km (contour, in m s−1) along with peak PRE (magenta line, in m s−1) over (10° N, 100° E; GAIA co-ordinates) near Chumphon in Thailand during March and April of 2012 at 1700 LT and 1500 LT

Figure 5 represents the cross-correlation coefficient of the peak PRE with (a) electron density, (b) eastward current density and (c) eastward neutral wind at altitudes of 110 km, 250 km, 300 km and 350 km at different local time (1100 LT, 1300 LT, 1500 LT, 1700 LT and 1900 LT) over the magnetic equator during the equinoctial months of 2011–2013. In case of plasma density, the correlation coefficient goes on decreasing from 1100 LT (around 0.3) to 1900 LT (ranging from − 0.4 to − 0.6 for different altitudes). The cross-correlation coefficients with eastward current density increases with local time from approximately 0.2–0.35 at 1100 LT to 0.4 at 1700 LT. The maximum correlation coefficient value of 0.4 is obtained at 1700 LT for all the altitudes 110 km, 250 km, 300 km and 350 km. Therefore, it is observed that the pre-sunset (at 1700 LT) eastward current density can predict the generation of PRE in GAIA as the peak PRE is observed around 1830 to 1900 LT mostly. In case of neutral wind, the maximum correlation is seen in the F-region (250, 300 and 350 km) with maximum correlation coefficient of 0.5 at 1900 LT.

Fig. 5
figure 5

Correlation coefficients of the day-to-day variation of the peak PRE with a electron density, b eastward current density, and c eastward neutral wind at different local time at the altitudes of 110 km (red line), 250 km (blue line), 300 km (green line) and 350 km (magenta line) over (10° N, 100° E; GAIA co-ordinates) near Chumphon in Thailand

Figure 6 displays the scatter plot showing the correlation of eastward current density at 1700 LT versus the peak PRE (during the equinoctial months of 2011–2013) at the altitude of 110 km over (10° N, 100° E; GAIA co-ordinates) near Chumphon in Thailand. It is found that the peak PRE and eastward current density depicts a positive correlation with correlation coefficient (R) of 0.407.

Fig. 6
figure 6

Scatter plot showing the correlation of eastward current density around 1700 LT versus the peak PRE (during the equinoctial months of 2011–2013) at the altitude of 110 km over (10° N, 100° E; GAIA co-ordinates) near Chumphon in Thailand

Discussion

We have investigated cross-correlations between day-to-day variation of the peak PRE and that of the electron density, eastward current density and eastward neutral wind, and found the followings: (1) F-region winds at the time of the peak PRE is positively correlated with the peak PRE. The cross-correlation coefficient is 0.5. (2) The eastward current density at an altitude of 110 km has positive correlation with peak PRE. Their cross-correlation coefficient is approximately 0.4 at 1700 LT. (3) The electron density in the F-region shows negative correlation with peak PRE. Below, the possible mechanisms of PRE and its relationship with the above parameters will be discussed.

The eastward neutral winds in the F-region (250, 300 and 350 km) show significant correlation with the peak PRE with correlation coefficient of 0.5 (approximately) at 1900 LT while the eastward electric current depicts correlation coefficient of 0.4 (approximately) at 1700 LT (pre-sunset). According to Farley et al. (1986), mechanism for generating PRE is explained below. The neutral wind in the F-region (U) blows eastward at the evening terminator because the neutral winds in the F-region is driven by the pressure gradient of the neutral atmosphere due to the solar heating. Due to the U, which blows across the geomagnetic field (\(\varvec{B}\)), downward polarization electric field (\(\varvec{E}_{\phi }\)), which is in the direction of \(- \varvec{U} \times \varvec{B}\), is generated through the F-region dynamo mechanism (Rishbeth 1971). \(\varvec{E}_{\phi }\) maps along \(\varvec{B}\) to the E-region. In the E-region, \(\varvec{E}_{\phi }\) has an equatorward component so that it induces a westward Hall current. The Hall current flows in the dayside because of high conductivity, but it does not flow in the nightside because the plasma density in the E-region disappears rapidly after the sunset. Due to the discontinuity of the Hall current at the evening terminator, negative charges are accumulated at the evening terminator. Consequently, eastward and westward electric fields are generated in the dayside and nightside, respectively. These electric fields are transmitted to the F-region. The eastward electric field in the dayside corresponds to PRE. Considering this mechanism, it is expected that the peak PRE could be intense when \(\varvec{U}\) is intense. The intense eastward electric fields make the growth rate of the RTI large, and thus EPB likely occurs. Otsuka et al. (2006) has pointed out close relationship between the eastward drift velocity and EPB occurrence. Their observations show that EPB occurrence rate is higher in March equinox than September equinox, and that the eastward plasma drift velocity (downward electric field) is larger in March equinox than September equinox. They have studied relationship between the downward electric field and EPB occurrence at the basis of seasonal variation. Our study suggest that the eastward neutral wind, which drive the downward electric field through the F-region dynamo mechanism, is closely related to the peak PRE, and that this scenario is applicable to their day-to-day variations. In this study, we have also investigated contribution of the zonal winds in the E-region on PRE because the E-region winds drive Hall current, which is not considered in the theory of Farley et al. (1986). Our results show that the E-region winds over magnetic equator have very less correlation with PRE. On the other hand, using WACCM-X model, Liu (2020) has shown that the day-to-day variation of PRE over the July period, when PRE is weak, is strongly affected by the large-scale lower atmosphere variability through the E-region dynamo at middle latitudes in the summer hemisphere. These results suggest that the PRE may be controlled by global distribution of electric fields.

GAIA model has shown negative correlation between the peak PRE and electron density in the F-region at magnetic equator. The eastward electric field of PRE elevates the F-region plasma at the magnetic equator to higher altitudes by E × B drift. The plasma moves down along the magnetic field lines under the influence of gravitational forces. This phenomenon is known as the “equatorial fountain”. This fountain effect causes plasma density depression or trough at the magnetic equator fringed by two nearly symmetric crests at conjugate low magnetic latitudes (± 15°), characterizing the equatorial ionization anomaly (EIA) (Moffett and Hanson 1965). Consequently, negative correlation between the PRE and F-region electron density over the magnetic equator represents the trough of plasma density caused by the equatorial fountain. When PRE is intense, the F-region electron density decreases at the magnetic equator.

GAIA model has succeeded in reproducing EEJ, which is an intense electro jet flowing at an altitude of approximately 100 km during daytime. GAIA model shows that EEJ is correlated with the peak PRE with cross-correlation coefficient of approximately 0.4 at 1700 LT and approximately 0.3 at 1100 LT. It should be noted that the PRE occurs between 1830 and 1900 LT. Therefore, EEJ can be utilized to predict the probable development of PRE observationally. Uemoto et al. (2010) have compared EEJ with the PRE strength and the onset of ESF over Chumphon during a period from November 2007 to October 2008, and shown that EEJ 1–2 h prior to sunset is positively correlated with PRE strength and spread-F occurrence. This led them to conclude that the pre-sunset E-region dynamo and/or electric field can be related to the F-region dynamics and onset of ESF around sunset. The present results show that GAIA reproduced well the relationship between the EEJ and PRE.

Summary and conclusions

Using an atmosphere–ionosphere coupled model (GAIA), we have studied what parameter mainly controls day-to-day variation of PRE. We have compared day-to-day variations of the peak PRE with electron density, eastward current density and eastward neutral wind over magnetic equator at (10° N, 100° E) near Chumphon, Thailand during the equinoctial months (March, April, September and October) in 2011–2013 during a high solar activity period. The important findings in this study are described as follows:

  1. 1.

    The F-region winds (at altitudes of 250–350 km) are positively correlated with the peak PRE, and their cross-correlation coefficient is approximately 0.4 when the PRE reaches a peak. This result indicates that the F-region winds control the peak PRE through the mechanisms of the F-region dynamo and PRE including E- and F-region coupling processes. On the other hand, correlation of E-region winds with the peak PRE is not seen.

  2. 2.

    The eastward current density at an altitude of 110 km during daytime, which is EEJ, has positive correlation with peak PRE. The cross-correlation coefficient between the current density and peak PRE increases with time until pre-sunset, and is approximately 0.4 at 1700 LT. This result is consistent with the previous observations. It should be noted that the EEJ strength can be used to predict magnitude of the PRE a few hours prior to the sunset.

  3. 3.

    The electron density in the F-region shows negative correlation with peak PRE. This could be due to transport of the plasma density to the low latitudes by ‘fountain effect’ which is responsible for the EIA.

The present study shows electro-dynamical forces coupled between E- and F-region could play an important role in day-to-day variation of the peak PRE.

Availability of data and materials

To access the GAIA simulation data, contact H. Jin (jin@nict.go.jp).

Abbreviations

EEJ:

Equatorial electro jet

EIA:

Equatorial ionization anomaly

EPB:

Equatorial plasma bubble

ESF:

Equatorial Spread-F

EUV:

Extreme ultraviolet

GAIA:

Ground-to-topside model of atmosphere–ionosphere for aeronomy

LT:

Local time

PRE:

Pre-reversal enhancement

RTI:

Rayleigh–Taylor instability

TIEGCM:

Thermosphere/ionosphere/electrodynamic general circulation model

References

Download references

Acknowledgements

Simulations and data analyses in this work were performed using Hitachi SR16000/M1 and the NICT Science Cloud System, Japan. The authors duly acknowledge the reviewers for their valuable comments.

Funding

The present study was supported by the Japan Society for the Promotion of Science, KAKENHI (15H05815), Project for Solar-Terrestrial Environment Prediction (PSTEP).

Author information

Authors and Affiliations

Authors

Contributions

PG carried out the data analysis and wrote the manuscript. YO conceived the idea, involved in elucidation of data analysis and organization of the paper. SM was engaged in the data interpretation. HS participated in reviewing the paper and is responsible person for GAIA model simulations. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Priyanka Ghosh.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ghosh, P., Otsuka, Y., Mani, S. et al. Day-to-day variation of pre-reversal enhancement in the equatorial ionosphere based on GAIA model simulations. Earth Planets Space 72, 93 (2020). https://doi.org/10.1186/s40623-020-01228-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40623-020-01228-9

Keywords