Equatorial ionosphere–thermosphere system: Electrodynamics and irregularities

The equatorial ionosphere and thermosphere constitute a coupled system, with its electro dynamical and plasma physical processes being responsible for a variety of ionospheric phenomena peculiar to the equatorial region. The most important of these phenomena are: the equatorial electrojet (EEJ) current system and its instabilities, the equatorial ionization anomaly (EIA), and the plasma instabilities/irregularities of the night ionosphere (associated with the plasma bubble events – ESF). They constitute the major topics of investigations having both scientific and practical objectives. The tidal wind interaction with the geomagnetic field is responsible for the atmospheric dynamo electric fields, that together with the wind system, drives the major phenomena, under quiet conditions. Drastic modifications of these phenomena can occur due to magnetospheric forcing under solar-, interplanetary- and magnetospheric disturbances. They can also undergo significant modifications due to forcing by atmospheric waves (such as planetary- and atmospheric gravity waves) propagating upward or from extra tropics. This article will focus on the ambient conditions of the ionosphere–thermosphere system and the electro dynamics and plasma instability processes that govern the plasma irregularity generation. Major emphasis is given to problems related to the structuring of the equatorial night ionosphere through plasma bubble/ESF irregularity processes. Specific topics to be covered will include: equatorial electric fields, thermospheric winds, sunset electrodynamic processes, plasma drifts, EEJ plasma instability/irregularity generation, nighttime/post sunset plasma bubble irregularity generation, and very briefly, disturbance electric fields and winds and their effect on the ionization anomaly, the TEC and ESF/plasma bubble irregularities.     

On the uncertainties in the measurement of absolute (true) TEC over Indian equatorial and low latitude sectors

The Indian sector encompasses the equatorial and low latitude regions where the ionosphere is highly dynamic and is characterized by the equatorial ionization anomaly (EIA) resulting in large latitudinal electron density gradients causing errors and uncertainties in the estimation of range delays in satellite based navigation systems. The diurnal and seasonal variations of standard deviations in the TEC data measured during the low sunspot period 2004–2005 at 10 different Indian stations located from equator to the anomaly crest region and beyond are examined and presented. The day-to-day variability in TEC is found to be lowest at the equatorial station and increases with latitude up to the crest region of EIA and decreases beyond.     

Equatorial Ionospheric Anomaly (EIA) and comparison with IRI model during descending phase of solar activity (2005–2009)

The ionospheric variability at equatorial and low latitude region is known to be extreme as compared to mid latitude region. In this study the ionospheric total electron content (TEC), is derived by analyzing dual frequency Global Positioning System (GPS) data recorded at two stations separated by 325 km near the Indian equatorial anomaly region, Varanasi (Geog latitude 25°, 16^/ N, longitude 82°, 59^/ E, Geomagnetic latitude 16°, 08^/ N) and Kanpur (Geog latitude 26°, 18^/ N, longitude 80°, 12^/ E, Geomagnetic latitude 17°, 18^/ N). Specifically, we studied monthly, seasonal and annual variations as well as solar and geomagnetic effects on the equatorial ionospheric anomaly (EIA) during the descending phase of solar activity from 2005 to 2009. It is found that the maximum TEC (EIA) near equatorial anomaly crest yield their maximum values during the equinox months and their minimum values during the summer. Using monthly averaged peak magnitude of TEC, a clear semi-annual variation is seen with two maxima occurring in both spring and autumn. Results also showed the presence of winter anomaly or seasonal anomaly in the EIA crest throughout the period 2005–2009 only except during the deep solar minimum year 2007–2008. The correlation analysis indicate that the variation of EIA crest is more affected by solar activity compared to geomagnetic activity with maximum dependence on the solar EUV flux, which is attributed to direct link of EUV flux on the formation of ionosphere and main agent of the ionization. The statistical mean occurrence of EIA crest in TEC during the year from 2005 to 2009 is found to around 12:54 LT hour and at 21.12° N geographic latitude. The crest of EIA shifts towards lower latitudes and the rate of shift of the crest latitude during this period is found to be 0.87° N/per year. The comparison between IRI models with observation during this period has been made and comparison is poor with increasing solar activity with maximum difference during the year 2005.     

Longitudinal behaviors of the IRI-B parameters of the equatorial electron density profiles retrieved from FORMOSAT-3/COSMIC radio occultation measurements

The electron density profiles in the bottomside F₂-layer ionosphere are described by the thickness parameter B0 and the shape parameter B1 in the International Reference Ionosphere (IRI) model. We collected the ionospheric electron density (N_e) profiles from the FORMOSAT-3/COSMIC (F3/C) radio occultation measurements from DoY (day number of year) 194, 2006 to DoY 293, 2008 to investigate the daytime behaviors of IRI-B parameters (B0 and B1) in the equatorial regions. Our fittings confirm that the IRI bottomside profile function can well describe the averaged profiles in the bottomside ionosphere. Analysis of the equatorial electron density profile datasets provides unprecedented detail of the behaviors of B0 and B1 parameters in equatorial regions at low solar activity. The longitudinal averaged B1 has values comparable with IRI-2007 while it shows little seasonal variation. In contrast, the observed B0 presents semiannual variation with maxima in solstice months and minima in equinox months, which is not reproduced by IRI-2007. Moreover, there are complicated longitudinal variations of B0 with patterns varying with seasons. Peaks are distinct in the wave-like longitudinal structure of B0 in equinox months. An outstanding feature is that a stable peak appears around 100°E in four seasons. The significant longitudinal variation of B0 provides challenges for further improving the presentations of the bottomside ionosphere in IRI.     

Variation of ionospheric total electron content in Indian low latitude region of the equatorial anomaly during May 2007–April 2008

The ionospheric total electron content (TEC), derived by analyzing dual frequency signals from the Global Positioning System (GPS) recorded near the Indian equatorial anomaly region, Varanasi (geomagnetic latitude 14°, 55′N, geomagnetic longitude 154°E) is studied. Specifically, we studied monthly, seasonal and annual variations as well as solar and geomagnetic effects on the equatorial ionospheric anomaly (EIA) during the solar minimum period from May 2007 to April 2008. It is found that the daily maximum TEC near equatorial anomaly crest yield their maximum values during the equinox months and their minimum values during the summer. Using monthly averaged peak magnitude of TEC, a clear semiannual variation is seen with two maxima occurring in both spring and autumn. Statistical studies indicate that the variation of EIA crest in TEC is poorly correlated with Dst-index (r = −0.03) but correlated well with Kp-index (r = 0.82). The EIA crest in TEC is found to be more developed around 12:30 LT.     

Diagnostics of equatorial and low latitude ionosphere by TEC mapping over Brazil

The total electron content (TEC) in the equatorial and low-latitude ionosphere over Brazil was monitored in two dimensions by using 2011 data from the ground-based global navigation satellite system (GNSS) receiver network operated by the Brazilian Institute for Geography and Statistics. It was possible to monitor the spatial and temporal variations in TEC over Brazil continuously during both day and night with a temporal interval of 10 min and a spatial resolution of about 400 km. The daytime equatorial ionization anomaly (EIA) and post-sunset plasma enhancement (PS-EIA) were monitored over an area corresponding to a longitudinal extension of 4000 km in South America. Considerable day-to-day variation was observed in EIA and PS-EIA. A large latitudinal and longitudinal gradient of TEC indicated a significant ionospheric range error in application of the GNSS positioning system. Large-scale plasma bubbles after sunset were also mapped over a wide range. Depletions with longitudinally separated by more than 800 km were observed. They were extended by more than 2000 km along the magnetic field lines and drifted eastward. It is expected that 2-dimensional TEC mapping can serve as a useful tool for diagnosing ionospheric weather, such as temporal and spatial variation in the equatorial plasma trough and crest, and particularly for monitoring the dynamics of plasma bubbles.     

IMF polarity effects on the equatorial ionospheric F-region

An exploratory study is made of the influence, during the equinoxes, of the interplanetary magnetic field (IMF) sector structure on the ionospheric F-region using ionosonde data from several equatorial stations for a 3-yr period around the 19th sunspot cycle maximum. It is found that, compared with days having positive IMF polarity, the post-sunset increase of h'F near the dip equator and the depth of the equatorial ionization anomaly (EIA) are reduced during the vernal equinox and enhanced during the antumnal equinox on days with negative IMF polarity. Similar trends are also noted in the data for the 20th sunspot cycle maximum, but with reduced amplitude. The systematic changes in the F-region characteristics suggest a modification of the equatorial zonal electric fields in association with the IMF polarity-related changes in the semi-annual variation of geomagnetic activity.     

Effects of CME and CIR induced geomagnetic storms on low-latitude ionization over Indian longitudes in terms of neutral dynamics

This paper presents the response of the ionosphere during the intense geomagnetic storms of October 12–20, 2016 and May 26–31, 2017 which occurred during the declining phase of the solar cycle 24. Total Electron Content (TEC) from GPS measured at Indore, Calcutta and Siliguri having geomagnetic dips varying from 32.23°N, 32°N and 39.49°N respectively and at the International GNSS Service (IGS) stations at Lucknow (beyond anomaly crest), Hyderabad (between geomagnetic equator and northern crest of EIA) and Bangalore (near magnetic equator) in the Indian longitude zone have been used for the storms. Prominent peaks in diurnal maximum in excess of 20–45 TECU over the quiet time values were observed during the October 2016 storm at Lucknow, Indore, Hyderabad, Bangalore and 10–20 TECU for the May 2017 storm at Siliguri, Indore, Calcutta and Hyderabad. The GUVI images onboard TIMED spacecraft that measures the thermospheric O/N₂ ratio, showed high values (O/N₂ ratio of about 0.7) on October 16 when positive storm effects were observed compared to the other days during the storm period. The observed features have been explained in terms of the O/N₂ ratio increase in the equatorial thermosphere, CIR-induced High Speed Solar Wind (HSSW) event for the October 2016 storm. The TEC enhancement has also been explained in terms of the Auroral Electrojet (AE), neutral wind values obtained from the Horizontal Wind Model (HWM14) and equatorial electrojet strength from magnetometer data for both October 2016 and May 2017 storms. These results are one of the first to be reported from the Indian longitude sector on influence of CME- and CIR-driven geomagnetic storms on TEC during the declining phase of solar cycle 24.     

Variability of foE in the equatorial ionosphere with solar activity

This research examined the variability of f_oE in the equatorial ionosphere with solar activity within the equatorial ionospheric anomaly region. Ionosonde data recorded at Ouagadougou (lat. 12.4°N, long. 1.5°W and magnetic dip 1.43°N) were engaged to study the transient variations of the critical frequency of the E-layer (f_oE) and its dependence on solar activity. The study revealed that f_oE increases with the increase in solar intensity of the sun. The variability of the f_oE decreases with increases in the solar activity. The maximum value of the f_oE is at local noon when the ionosphere is stable; the variability at this local time is minimal. The minimum value of the f_oE is at sunrise and sunset, at this period on local time the equatorial ionosphere recorded its maxima variability. Irrespective of the degree of solar activity, f_oE is observed to be maximum in June solstice, followed by the equinoxes and minimum in December solstice. Equinoctial asymmetry occurred in the variation of the relative standard deviation of f_oE with maximum in September/March equinox for low/high solar activity.     

Variations of equatorial electrojet as possible seismo-ionospheric precursor at the occurrence of TEC anomalies before strong earthquake

The measurements of GPS signal delays show that the local areas of increased/decreased Total Electron Content (TEC) of the ionosphere can be observed before strong earthquakes. The main possible cause of these TEC disturbances is the vertical plasma drift under the action of zonal electric field. The spatial pattern of electric potentials for such electric field was proposed. The model calculations were done to investigate the efficiency of the proposed mechanism. The calculation results revealed the agreement with TEC variations observed before strong earthquakes and showed that the equatorial electrojet variations can be considered as precursors of earthquakes.     

Investigation on the MLT tidal variability during September 2019 minor sudden stratospheric warming

Tidal variability in the mesosphere and lower thermosphere (MLT) during September 2019 Southern hemisphere minor sudden stratospheric warming (SSW) is investigated utilizing ground-based meteor radar wind observations from the equatorial, extratropical, middle, and high latitude stations and global reanalysis dataset. The polar warming is found to move from the mesosphere to the stratosphere until the peak warming day (PWD) of the SSW. The diurnal and semidiurnal tides at individual observational sites do not exhibit any consistent response during the observational interval, but a notable and consistent variability in some specific zonal wavenumber components, i. e., DW1 (migrating diurnal tide), DE3 (nonmigrating eastward wavenumber 3 diurnal tide), and SW2 (migrating semidiurnal tide) is found in the global reanalysis dataset. Incidentally, the warming event occurs during Spring equinox when a dominant seasonal change in the tidal activities generally takes place and hence seasonal variability is also looked into while identifying the SSW impact during the observational interval. It is found that the seasonal broad changes in the DW1, DE3, and SW2 amplitudes can be explained by the variability in the tidal sources, i.e., water vapor, convective activity, ozone, etc during the observational period. However, the extracted short-term variability in the global tidal modes on removing seasonal trend reveals noticeable response in connection with the warming event. The deseasoned amplitude of the DW1 significantly enhances around the PWD at most of the present latitudes. The deseasoned DE3 amplitude responds significantly in the middle atmosphere at low latitudes during the warming phase. The deseasoned SW2 exhibit clear enhancement around the PWD at all the latitudes. However, the deseasoned tidal features do not seem to correlate well with that of the source species unlike the seasonal ones that imply involvement of complex processes during the warming event, seeking further future investigations in this regard.     

Dayside ionospheric uplift during strong geomagnetic storms as detected by the CHAMP,SAC-C,TOPEX and Jason-1 satellites

Ionosphere response to severe geomagnetic storms that occurred in 2001–2003 was analyzed using data of global ionosphere maps (GIM), altimeter data from the Jason-1 and TOPEX satellites, and data of GPS receivers on-board CHAMP and SAC-C satellites. This allowed us to study in detail ionosphere redistribution due to geomagnetic storms, dayside ionospheric uplift and overall dayside TEC increase. It is shown that after the interplanetary magnetic field turns southward and intensifies, the crests of the equatorial ionization anomaly (EIA) travel poleward and the TEC value within the EIA area increases significantly (up to ∼50%). GPS data from the SAC-C satellite show that during the main phase of geomagnetic storms TEC values above the altitude of 715 km are 2–3 times higher than during undisturbed conditions. These effects of dayside ionospheric uplift occur owing to the “super-fountain effect” and last few hours while the enhanced interplanetary electric field impinged on the magnetopause.     

The equatorial ionization anomaly at the topside F region of the ionosphere along 75°E

Electron density measured by the Indian satellite SROSS C2 at the altitude of ∼500 km in the 75°E longitude sector for the ascending half of the solar cycle 22 from 1995 to 1999 are used to study the position and density of the equatorial ionization anomaly (EIA). Results show that the latitudinal position and peak electron density of the EIA crest and crest to trough ratios of the anomaly during the 10:00–14:00 LT period vary with season and from one year to another. Both EIA crest position and density are found to be asymmetric about the magnetic equator and the asymmetry depends on season as well as the year of observation, i.e., solar activity. The latitudinal position of the crest of the EIA and the crest density bears good positive correlation with F_10.7 and the strength of the equatorial electrojet (EEJ).     

赤道地区热层纬向风反转时间和风速的经度分布差异

使用CHAMP卫星观测资料和TIEGCM模拟数据,研究了太阳活动低年春秋分季节赤道地区热层纬向风的反转时间和风速在经度分布上的差异,重点分析了离子拖曳力在不同地磁场构型下对纬向风速地理经度分布的影响.研究发现,纬向风一般在清晨时段由东向转为西向,下午时段由西向转为东向,但反转时间存在明显的经度差异,最大经度差异可达1....     

Excitations of nonmigrating diurnal tides in the mesosphere and lower thermosphere simulated by the Kyushu-GCM

Excitation mechanisms of nonmigrating diurnal tides in the MLT region simulated by the Kyushu-GCM are examined. It is shown that the westward propagating diurnal tide with zonal wavenumber s = 2 is mainly excited by nonlinear interactions between the migrating diurnal tide and the stationary planetary wave with zonal wavenumber s = 1, while the nonlinear excitation of the standing diurnal tide with zonal wavenumber s = 0 is less important than the excitation by tropospheric heating. Nonlinear interactions between the migrating diurnal tide and the stationary planetary wave with zonal wavenumber s = 2 are not dominant to excite the westward propagating diurnal tide with zonal wavenumber s = 3, and it is shown that the excitation by tropospheric heating is comparable to the nonlinear excitation. It is also shown that other nonmigrating diurnal tides are excited by tropospheric heating.     

The solar wind control of the equatorial ionosphere dynamics during geomagnetic storms

The interplanetary magnetic field, geomagnetic variations, virtual ionosphere height h′F, and the critical frequency foF2 data during the geomagnetic storms are studied to demonstrate relationships between these phenomena. We study 5-min ionospheric variations using the first Western Pacific Ionosphere Campaign (1998–1999) observations, 5-min interplanetary magnetic field (IMF) and 5-min auroral electrojets data during a moderate geomagnetic storm. These data allowed us to demonstrate that the auroral and the equatorial ionospheric phenomena are developed practically simultaneously. Hourly average of the ionospheric foF2 and h′F variations at near equatorial stations during a similar storm show the same behavior. We suppose this is due to interaction between electric fields of the auroral and the equatorial ionosphere during geomagnetic storms. It is shown that the low-latitude ionosphere dynamics during these moderate storms was defined by the southward direction of the B_z-component of the interplanetary magnetic field. A southward IMF produces the Region I and Region II field-aligned currents (FAC) and polar electrojet current systems. We assume that the short-term ionospheric variations during geomagnetic storms can be explained mainly by the electric field of the FAC. The electric fields of the field-aligned currents can penetrate throughout the mid-latitude ionosphere to the equator and may serve as a coupling agent between the auroral and the equatorial ionosphere.     

Variation of F-region critical frequency (foF2) over equatorial and low-latitude region of the Indian zone during 19th and 20th solar cycle

The effect of solar cycle and seasons on the daytime and nighttime F-layer ionization has been investigated over the equatorial and low-latitude region during 19th (1954–1964) and 20th (1965–1976) solar cycle. The F-layer critical frequency (foF2) data observed from the three Indian Ionosonde stations has been used for the present study. The dependence of foF2 on solar cycle has been examined by performing regression analysis between the foF2 values and R₁₂ (twelve month running average sunspot number). The result shows that the magnitude of the cycle, seasons and the location of station has considerable effects on foF2. There is a significant nonlinear relationship between the foF2 values and R₁₂ during 19th solar cycle as compared to 20th solar cycle. Further, the nighttime saturation effect is prominently seen during the 19th solar cycle and summer season. It is also observed that the most profound saturation effect appears at the equatorial ionization anomaly crest region. Seasonally, it is seen that all the stations exhibits semiannual anomaly. The phenomenon of winter anomaly decays as we move higher along the latitude and is prominently seen during the intense solar activity.     

Day-to-day variability of equatorial anomaly in GPS-TEC during low solar activity period

The ionospheric total electron content (TEC) in the northern hemispheric equatorial ionization anomaly (EIA) crest region is investigated by using dual-frequency signals of the Global Positioning System (GPS) acquired from Rajkot (Geog. Lat. 22.29°N, Geog. Long. 70.74°E; Geom. Lat. 14.21°N, Geom. Long. 144.90°E), India. The day-to-day variability of EIA characteristics is examined during low solar activity period (F10.7∼83 sfu). It is found that the daily maximum TEC at EIA crest exhibits a day-to-day and strong semi-annual variability. The seasonal anomaly and equinoctial asymmetry in TEC at EIA is found non-existent and weaker, respectively. We found a moderate and positive correlation of daily magnitude of crest, Ic with daily F10.7 and EUV fluxes with a correlation coefficient of 0.43 and 0.33, respectively indicating an existence of a short-term relation between TEC at EIA and the solar radiation even during low solar activity period. The correlation of daily Ic with Dst index is also moderate (r = −0.35), whereas no correlation is found with the daily Kp index (r = 0.14) respectively. We found that the magnitude of EIA crest is moderately correlated with solar flux in all seasons except winter where it is weakly related (0.27). The magnitude of EIA crest is also found highly related with EEJ strength in spring (r = 0.69) and summer (r = 0.65) than autumn (0.5) and winter (r = 0.47), though EEJ is stronger in autumn than spring.     

Validation of the Neustrelitz Global Model according to the low latitude ionosphere

Empirical modeling including empirical model for the total electron content (TEC) is important for the study of the ionosphere and practical applications. In this paper goodness of new Neustrelitz Global Model (then NGM) at low latitudes is studied. The NGM model includes such parameters as the maximal electron density (NmF2) and altitude of the maximum (hmF2). As of today, besides NGM there are several empirical models for NmF2 and hmF2. Therefore, a comparison of these parameters of the NGM model, not only with the experimental data, but also with two versions of the International Reference Ionosphere (the IRI model): IRI2001 and IRI-Plas would be instructive. Because the NGM model incorporates special factor describing the equatorial anomaly, the comparison in lower latitude areas is particularly interesting. As one can see from the presented example of the data from low latitude stations located in the northern and southern hemispheres near the Greenwich meridian, the NGM model may have certain advantages over the IRI model versions. In particular, NGM TEC is preferable regardless of solar activity level while NGM NmF2 is only preferable under high solar activity conditions. Next, NGM equivalent slab thickness of the ionosphere: τ(NGM) = TEC(NGM)/NmF2(NGM) has been calculated and tested to answer the question whether τ(NGM) can be used as a proxy of the slab thickness of the ionosphere for an empirical modeling. The answer is positive for the near equatorial stations and periods of high solar activity, and under such conditions predicted τ(NGM) can be used for deriving NmF2 from the experimental values of TEC(CODE) in real time.     

Variation of ionospheric total electron content in Taiwan region of the equatorial anomaly from 1994 to 2003

The ionospheric total electron content (TEC) in the northern hemispheric equatorial ionospheric anomaly (EIA) region is studied by analyzing dual-frequency signals of the Global Position System (GPS) acquired from a chain of nine observational sites clustered around Taiwan (21.9–26.2°N, 118.4–112.6°E). In this study, we present results from a statistical study of seasonal and geomagnetic effects on the EIA during solar cycle 23: 1994–2003. It is found that TEC at equatorial anomaly crests yield their maximum values during the vernal and autumnal months and their minimum values during the summer (except 1998). Using monthly averaged I_c (magnitude of TEC at the northern anomaly crest), semi-annual variations is seen clearly with two maxima occurring in both spring and autumn. In addition, I_c is found to be greater in winter than in summer. Statistically monthly values of I_c were poorly correlated with the monthly Dst index (r = −0.22) but were well correlated with the solar emission F_10.7 index (r = 0.87) for the entire database for the period during 1994–2003. In contrast, monthly values of I_c were correlated better with Dst (r ? 0.72) than with F_10.7 (r ? 0.56) in every year during the low solar activity period (1994–1997). It suggests that the effect of solar activity on I_c is a longer term (years), whereas the effect of geomagnetic activity on I_c is a shorter term (months).