Investigating Terrestrial Weather Impacts on the Ionosphere‐Thermosphere‐Mesosphere System Using Satellite Observations, Assimilative Modeling, and Climatological Tidal Modeling

TIMED Doppler interferometer (TIDI) observations of migrating diurnal and semidiurnal tides

The global ionospheric response to a stratospheric sudden warming (SSW) is studied using three‐dimensional electron density maps derived from radio occultation observations of FORMOSAT‐3/COSMIC during the 2009 SSW periods. Results show that the ionospheric electron density at EIA crests exhibit a morning/early afternoon increase followed by an afternoon decrease and an evening increase, indicative of a semidiurnal component during the SSW period, which is consistent with recent studies. The latitude‐altitude electron density slice maps show that the SSW related modifications of the equatorial plasma fountain interact with the existing summer‐to‐winter neutral winds and resulting in a north–south asymmetry. The global ionospheric response shows a clear longitudinal dependence in the equatorial plasma fountain enhancement during morning/early afternoon, inferred from the duration of the equatorial ionization anomaly (EIA) enhancement. Following the enhancement, prominent global EIA reductions resulting from the equatorial plasma fountain weakening in the afternoon sector are seen. The ionospheric response to the 2009 SSW event is also compared with the usual seasonal variation during January–February 2007. Instead of showing the electron density increase in the northern hemisphere and decrease in the southern hemisphere as the usual seasonal variation does, the SSW period ionosphere shows prominent global electron density reductions in the afternoon period during the 2009 SSW event.

Variability of the equatorial ionization anomaly over the South American sector: Effects of electric field and effective meridional wind

Ground-based GNSS (Global Navigation Satellite System) reference stations lack the capacity to provide data for ocean regions with sufficient spatial-temporal resolution, limiting the detailed study of the equatorial ionization anomaly (EIA). Thus, this study collected kinematic multi-GNSS data on the ionospheric Total Electron Content (TEC) during two research cruises across the equator in the Western Pacific Ocean in 2014 (31 October–8 November) and 2015 (29 March–6 April). The purpose of the study was to use sufficient spatial–temporal resolution data to conduct a detailed analysis of the diurnal variation of the equatorial ionization anomaly in different seasons. The two-year data collected were used to draw the following conclusions. During the test in 2014, the EIA phenomenon in the Northern and Southern Hemispheres was relatively obvious. The maximum values occurred at local time (LT) 15:00 (~136TECu) and LT22:00 (~107TECu) in the Northern Hemisphere and at LT14:00 (100TECu) and LT22:00 (80TECu) in the Southern Hemisphere. During the test in 2015, the EIA in the Southern Hemisphere reached its maximum level at LT14:00 (~115TECu) and LT20:00 (~85TECu). However, the EIA phenomenon in the Northern Hemisphere was weakened, and a maximum value occurred only at LT 15:00 (~85TECu). The intensity contrast was reversed. The EIA phenomenon manifests a strong hemisphere asymmetry in this region.

Interannual variability of the nonmigrating semidiurnal tide at high latitudes and stationary planetary wave in the opposite hemispheres

This paper presents an investigation of ionospheric response to great (Dst ≤−350 nT) geomagnetic storms that occurred during solar cycle 23. The storm periods analyzed are 29 March to 2 April 2001, 27–31 October 2003, 18–23 November 2003, and 6–11 November 2004. Global Navigation Satellite System, total electron content (TEC), and ionosonde critical frequency of F2 layer (foF2) data over Southern Hemisphere (African sector) and Northern Hemisphere (European sector) midlatitudes were used to study the ionospheric responses within 15°E–40°E longitude and ±31° to ±46° geomagnetic latitude. Midlatitude regions within the same longitude sector in both hemispheres were selected in order to assess the contribution of the low‐latitude changes especially the expansion of equatorial ionization anomaly (EIA) also called the dayside ionospheric superfountain effect during these storms. In all storm periods, both negative and positive ionospheric responses were observed in both hemispheres. Negative ionospheric responses were mainly due to changes in neutral composition, while the expansion of the EIA led to pronounced positive storm effects at midlatitudes for some storm periods. In other cases (e.g., 29 October 2003), penetration electric fields, EIA expansion, and large‐scale traveling ionospheric disturbances were found to be present during the positive storm effect at midlatitudes in both hemispheres. An increase in TEC on the 28 October 2003 was because of the large solar flare with previously determined intensity of X45 ±5.

Temporal evolution of the EIA along 95°E as obtained from GNSS TEC measurements and SAMI3 model

This study characterizes the African equatorial ionization anomaly (EIA) during the maximum phase of Solar Cycle 24 (2012–2015). Total electron content (TEC) data were obtained from a chain of 13 African Global Navigation Satellite System (GNSS) receivers (within 36–42°E geographic longitude; 29°N to 10°S geographic latitude; ±20°N magnetic latitude) to study quiet time variations of TEC and thereafter construct the EIA profiles. The correlations of TEC for pairs of conjugate stations within the African EIA were determined. TEC station‐to‐trough ratio (TEC‐STR) was used as an index to determine the strength of the EIA crests. Furthermore, the variability of the EIA crests on a 3‐hourly basis over each month was investigated. Overall, equinoxes recorded the highest values of TEC, with the widest latitudinal location of the EIA crests, and June solstice recorded the least, with the characteristic equatorward collapse of the EIA crests. TEC and the location of the crests of the EIA also increased with solar activity. TEC showed good correlations for conjugate pairs of stations, especially for stations closer to the trough. The African EIA morphology showed asymmetry in the locations of the Northern Hemisphere (NH) and Southern Hemisphere (SH) crests. The formation of EIA starts around 0600–0800 LT and decays around 2100–2300 LT. Generally, SH crest forms first before NH crest and also decays first; we attributed this to the effects of the meridional neutral winds. Comparing TEC at African EIA crests with other longitudinal sectors, the Asian sector recorded the highest TEC, while the American sector recorded the least.

Assessment of IRI-2012, NeQuick-2 and IRI-Plas 2015 models with observed equatorial ionization anomaly in Africa during 2009 sudden stratospheric warming event

In this study we present for the first time the nonmigrating tidal spectrum for the electron density in the topside ionosphere on global scale for different seasons at both solar maximum and minimum conditions. The electron density observations from the CHAMP (CHAllenging Minisatellite Payload) satellite provide evidence for prominent nonmigrating tides at different latitude regions. At middle and high latitudes the most prominent diurnal tides are DE1, D0, and DW2, as well as semidiurnal tides SW1 and SW3 preferably during equinox seasons. DE1 is only found in the northern middle and high latitudes, while D0 and DW2 are much stronger in the Southern Hemisphere. These tides are believed to be driven by the thermospheric winds through ion‐neutral interactions. At equatorial and low‐latitude regions the most prominent diurnal tides are DE2 and DE3. DE2 is found to be present at low‐latitude regions throughout the whole year with larger amplitudes in the Southern Hemisphere, while DE3 shows largest amplitudes (symmetric above the dip equator) at the equatorial ionization anomaly (EIA) crest region during September equinox. A general antiphase behavior between the EIA crest and trough is observed for the tides DE3, DE2, DW2, and SW3. We consider this as strong evidence for the modulation of the EIA electron density by the E layer zonal electric field via the ion fountain effect. An exception is the tidal component D0, which exhibits antiphase variations between the two hemispheres. The phase value at the equatorial trough is halfway between those of the two hemispheres. Presently, we cannot give an explanation for it, and study concerning special model effort is further needed.

Semidiurnal tides from the extended Canadian Middle Atmosphere Model (CMAM) and comparisons with TIMED Doppler interferometer (TIDI) and meteor radar observations

This study presents longitudinal structures of the mid- and low-latitude ionosphere using the GPS radio occultation observation on board the COSMIC satellite mission. The longitudinal structure seen in the equatorial and low-latitude ionospheric regions results from modification of the daily dynamo electric field by upward propagating atmospheric tides that are generated by latent heat release of the tropical rainstorms. Changes of the dynamo electric field modify the equatorial plasma fountain and thereby enhance the equatorial ionization anomaly (EIA). With capability of three-dimensional global ionospheric observation, altitudinal, local time, and monthly variations of this recent discovered fascinating feature are obtained for further understanding of the underlying physical mechanism. Through comparison between electron densities at various altitudes, the longitudinal structure is prominently seen at upper part of the ionosphere. Additionally, COSMIC observations provide three-dimensional structure of the Weddell Sea anomaly which is featured by the greater nighttime electron density than daytime. Not only occurring at the southern hemisphere near the Weddell Sea region of the Antarctica, a similar nighttime density enhancement feature is also found in the northern hemisphere during local summer by COSMIC observations. The anomalous signatures in both hemispheres share very similar characteristics in electron density structure, latitudinal distribution, and appearance time. They are, therefore, categorized as the mid-latitude summer nighttime anomaly (MSNA).

 Knowing the locations of the north and south Equatorial Ionization Anomaly (EIA) crests and their corresponding widths is essential for characterizing the spatiotemporal and solar activity variations of the low latitude ionosphere. The crest region is characterized by strong electron density gradients that significantly affect GNSS applications. However, there is still a lack of complete characterization and modeling of the spatiotemporal and solar activity variations of the EIA crest positions and widths. Therefore, the purpose of this study is to characterize and model the spatiotemporal and solar activity variations of EIA crest widths and positions. These characteristics are described and modeled using 13 years (2009-2017 and 2020-2023) NmF2 data, which are obtained from radio occultation electron density profiles of GRACE, COSMIC-1, and COSMIC-2 satellites over the globe. Crest positions and widths exhibit diurnal, semi-diurnal, and annual variations. There is a slight linear correlation between the solar activity and the crest positions and widths. Furthermore, longitudinal variations in the geomagnetic field declination and the interplay between wave number 3 diurnal tides and planetary waves may be related to longitudinal variations in the crest positions and widths. The results of the analysis are used to create semi-empirical EIA crest position and width models. These models provide a good description of the crest widths and positions verified by vTEC altimeter measurements. These models can serve as subroutines for the improvement of current ionospheric TEC and F2-layer maximum electron density models, applicable for improving terrestrial communications and GNSS positioning and navigation. 

The present article aims at a consistent understanding of observation, theoretical model, and simulation with the geomagnetic sudden commencement (SC) observed in the morning and afternoon at high and middle latitudes in the northern and southern hemispheres and at the noontime equator on 12 May 2021. The SC in Bx- and By-components of the geomagnetic field, SCx,y, was composed of the positive/negative preliminary (PI) and main impulses (MI) as SCx (+ -) and SCy (- +) in the morning and SCx (- +) and SCy (+ -) in the afternoon at middle latitudes in the northern hemisphere. SCx in the southern hemisphere is in the same polarity as those in the northern hemisphere, except for SCx (+ +) in the morning. SCy in the southern hemisphere has reverse polarity to those in the northern hemisphere. The PIx in the northern hemisphere matches the well-established two-cell Hall current vortices with anti-clockwise and clockwise directions in the morning and afternoon, respectively, and the MIx matches reverse Hall current vortices. The PIx and MIx in the southern hemisphere meet the Hall currents that are mirror images of those in the northern hemisphere with respect to the equator except for the positive MI in the morning. The PIy in the northern hemisphere is shown to meet the northward and southward Pedersen currents in the morning and afternoon, respectively, and the MIy meets reverse Pedersen currents. The PIy and MIy in the southern hemisphere are found to meet the Pedersen currents that are mirror images of those in the northern hemisphere. At the equator, typical SCx (- +) is observed, meeting the Cowling currents that should be supplied by the Pedersen currents responsible for the observed midlatitude SCy in the northern and southern hemispheres. The electric fields of the PI and MI observed by the HF Doppler sounders at the middle latitudes in the northern hemisphere are westward and eastward, respectively, in both the morning and afternoon, meeting the conventional dusk-to-dawn PI and dawn-to-dusk MI electric fields. The onset of the PI is found to be simultaneous with the resolution of a few seconds from high latitude to the equator in both the northern and southern hemispheres, indicating instantaneous achievement of the Pedersen–Cowling currents from high latitude to the equator. The instantaneous achievement of the energy-consuming Pedersen–Cowling currents is explained by the TM0/TEM mode wave in the Earth-ionosphere waveguide/transmission line rather than the compressional waves in the magnetosphere and F-region ionosphere. REPPU (REProduce Plasma Universe) global simulation model equipped with a potential solver at the inner boundary of the model magnetosphere reproduces the PI and MI electric fields at middle latitudes and SCx (- +) at the dayside equator. The simulation results are found to be consistent with most features of observations, such as the time scale of PI and MI, direction of the midlatitude electric field and generation of the Cowling currents. The simulation proves that the electric fields and FACs are generated in the outer magnetosphere, transmitted to the polar ionosphere and then to the equator in the Pedersen–Cowling current circuit.

We have retrieved the latitudinal and vertical structures of the 11 year solar cycle modulation on ionospheric electron density using 14 years of satellite‐based radio occultation measurements utilizing the Global Positioning System. The densities at the crests of the equatorial ionization anomaly (EIA) in the subtropics near 300 km in 2003 and 2014 (high solar activity with solar 10.7 cm flux, F10.7 ≈ 140 solar flux unit (sfu)) were 3 times higher than that in 2009 (low solar activity F10.7 ≈ 70 sfu). The higher density is attributed to the elevated solar extreme ultraviolet and geomagnetic activity during high solar activity periods. The location of the EIA crests moved ~50 km upward and ~10° poleward, because of the enhanced E × B force. The EIA in the northern hemisphere was more pronounced than that in the southern hemisphere. This interhemispheric asymmetry is consistent with the effect of enhanced transequatorial neutral wind. The above observations were reproduced qualitatively by the two benchmark runs of the Thermosphere‐Ionosphere‐Electrodynamics General Circulation Model. In addition, we have studied the impact of the 11 year solar cycle on the 27 day solar cycle response of the ionospheric electron density. Beside the expected modulation on the amplitude of the 27 day solar variation due to the 11 year solar cycle, we find that the altitude of the maximal 27 day solar response is unexpectedly ~50 km higher than that of the 11 year solar response. This is the first time that a vertical dependence of the solar responses on different time scales is reported.