Hemisphere-coupled modeling of nighttime medium-scale traveling ionospheric disturbances

Nighttime medium-scale traveling ionospheric disturbances (MSTIDs), which have tilted frontal structures in the midlatitude ionosphere, are investigated by the midlatitude ionosphere electrodynamics coupling (MIECO) model in this study. It has been proposed that the electrodynamic coupling between the E and F regions plays an important role in generating MSTIDs within a few hours. An intriguing aspect of MSTIDs is that they were simultaneously observed at magnetic conjugate locations in the Northern and Southern Hemispheres. In order to study the hemisphere-coupled electrodynamics, the MIECO model has been upgraded to consist of two simulation domains for both hemispheres in which the electrostatic potential is solved by considering electrodynamics in both hemispheres. The simultaneous occurrence of MSTIDs at the magnetic conjugate stations has clearly been reproduced when the F-region neutral wind satisfies the unstable condition in both hemispheres and a sporadic-E layer is given only at the Northern (summer) Hemisphere. Even if the unstable condition is satisfied in the summer hemisphere, an unfavorable F-region neutral wind in the winter hemisphere largely suppresses the growth of MSTIDs in both hemispheres.     

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.     

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.     

The glow of the night ionosphere due to energetic ion beams

This paper discusses photometric measurements made of the ionospheric excitation of the line $\text{λ = 5577}\text{A}\text{?}$ at the time of electron beam injection from a rocket into the Earth's ionosphere. The gradual increase of the glow intensity per impulse occurs due to accumulation of the energy of excited states of N₂(A³Σ⁺_u) and O(′S) during their lifetimes. The large disturbed zone in the near-rocket environment (size >500 m) is connected via the interaction of ions accelerated in the rocket potential field with ionospheric components. The glow intensity modulation is observed at a height of ～98 km during the electron beam injection simultaneously with the ignition of the beam-plasma discharge (BPD). The intensity minima are explained by a decrease of the energy of accelerated ions due to effective neutralization of the rocket body by the BPD plasma. The height profile of the glow intensity revealed two maxima at heights of ～103 km and ～115 km. The second maximum (at ～115 km) indicates that, at these heights, both collision and collision-free mechanisms of accelerated ion energy transport to ionospheric components exist.     

Collocated ionosonde and dense GPS/GLONASS network measurements of midlatitude MSTIDs

To analyze midlatitude medium-scale travelling ionospheric disturbances (MSTIDs) over Kazan (55.5°N, 49°E), Russia, the sufficiently dense network of GNSS receivers (more than 150 ground-based stations) were used. For the first time, daytime MSTIDs in the form of their main signature (band structure) on high-resolution two-dimensional maps of the total electron content perturbation (TEC maps) are compared with ionosonde data with a high temporal resolution. For a pair of events, a relationship between southwestward TEC perturbations and evolution of F2 layer traces was established. So F2 peak frequency varied in antiphase to TEC perturbations. The ionograms show that during the movement of plasma depletion band (overhead ionosonde) the F2 peak frequency is the highest, and vice versa, for the plasma enhancement band, the F2 peak frequency is the lowest. One possible explanation may be a greater inclination of the radio beam from the vertical during the placement of a plasma enhancement band above the ionosonde, as evidenced by the absence of multiple reflections and the increased occurrence rate of additional cusp trace. Another possible explanation may be the redistribution of the electron content in the topside ionosphere with a small decrease in the F peak concentration of the layer with a small increase in TEC along the line-of-sight. Analysis of F2 peak frequency variation shows that observed peak-to-peak values of TEC perturbation equal to 0.4 and 1 TECU correspond to the values of ΔN/N equal to 13% and 28%. The need for further research is evident.     

Global comparison of the model results of GSM TIP with IRI for summer conditions

This paper presents the results of the numerical calculations thermosphere/ionosphere parameters which were executed with using of the Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP)and comparison of these results with empirically-based model IRI-2001. Model GSM TIP was developed in West Department of IZMIRAN and solves self-consistently the time-dependent, 3-D coupled equations of the momentum, energy and continuity for neutral particles (O₂, N₂, O), ions (O⁺, H⁺), molecular ions (M⁺) and electrons and largescale eletric field of the dynamo and magnetospheric origin in the range of height from 80 km to 15 Earth’s radii. The empirically derived IRI model describes the E and F regions of the ionosphere in terms of location, time, solar activity and season. Its output provides a global specification not only of N_e but also on the ion and electron temperatures and the ion composition. These two models represent a unique set of capabilities that reflect major differences in along with a substantial approaches of the first-principles model and global database model for the mapping ionosphere parameters. We focus on global distribution of the N_e, T_i, T_e and TEC for the one moment UT and fixed altitudes: 110 km, h_mF2, 300 km and 1000 km. The calculations were executed with using of GSM TIP and IRI models for August 1999, moderate solar activity and quiet geomagnetic conditions. Results present as the global differences between the IRI and GSM TIP models predictions. The discrepancies between model results are discussed.     

Planetary wave coupling of the atmosphere–ionosphere system during the Northern winter of 2008/2009

This paper presents the global spatial (latitude and altitude) structure and temporal variability of the ∼23-day ionospheric zonally symmetric (s = 0) planetary wave (PW) seen in the Northern winter of 2008/2009 (October 2008–March 2009). It is shown that these ∼23-day ionospheric oscillations are forced from PWs propagating from below. The COSMIC ionospheric parameters f_oF2 and h_mF2 and electron density at fixed altitudes and the SABER temperatures were utilized in order to define the waves which are present simultaneously in the atmosphere and ionosphere. The long-period PWs from the two data sets have been extracted through the same data analysis method. The similarity between the lower thermospheric ∼23-day (s = 0) temperature PW and its ionospheric electron density response provides valuable and strong experimental evidence for confirming the paradigm of atmosphere–ionosphere coupling.     

Physical implications of discrepancy between summer and winter PC indices observed in the course of magnetospheric substorms

The PC index based on a statistically justified relationship between the polar cap magnetic activity and the interplanetary electric field E_KL has been derived as a value standardized for the E_KL intensity regardless of season, UT and hemisphere. As a result, the summer and winter PC indices are consistent with one another under ordinary conditions. Discrepancies between the summer and winter PC indices arising in the course of magnetospheric substorms are analyzed in this paper. It is argued that the channel of enhanced conductivity, formed in the auroral oval owing to intense auroral particle precipitation, strongly improves the conditions for closure of the Region 1 field-aligned currents in the winter dark polar region but only trivially affects the conditions of the Region 1 FAC closure in the summer sunlit ionosphere. Since the coefficients describing the relationship between E_KL and the polar cap magnetic activity were derived for statistically justified (i.e., mean) conditions, their application to such abnormal situation, as intense field-aligned currents in the winter dark polar region, leads to overestimation of the winter PC index. The summer and winter PC indices level off as soon as the intense auroral particle precipitation terminates and the auroral ionosphere in the winter and summer polar caps returns to the ordinary (statistically justified) state.     

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.     

Modeling of ionospheric scintillation at low-latitude

The presence and movement of plasma density fluctuations in the F-region of the ionosphere are studied by monitoring phase and amplitude of radio waves propagating through the region. In this paper, we have used weak scattering theory and assumed the plasma density fluctuations to behave like phase changing diffraction screen. Appropriate relations for scintillation index S₄, and phase variance δ? are derived and computed for different parameters of the plasma density irregularities of the ionosphere. SROSS-C2 satellite in situ measurements of plasma density fluctuations, which provide direct information about the structure and morphology of irregularities that are responsible for scintillation of radio waves, were used first time to develop a scintillation model for low latitude. It is observed that the scintillation index S₄ and phase variance δ? depends on the strength of the plasma turbulence. Finally, the results obtained from modeling are compared and discussed with the available recent results.     

The improvement of IRI profiles of O2+, NO+ and O+ by means of airglow measurements

Methods were suggested in /1,2,3/ to find the maximum electron density N_mF (or N_mF2) and the) corresponding height h_mF (h_mF2) in the ionosphere through a combination of measurements on two airglow lines: the UV oxygen line at λ135.6 nm and the red oxygen line near λ630 nm.     

Moderate geomagnetic storms of January 22–25, 2012 and their influences on the wave components in ionosphere and upper stratosphere-mesosphere regions

Moderate geomagnetic storms occurred during January 22–25, 2012 period. The geomagnetic storms are characterized by different indices and parameters. The SYM-H value on January 22 increased abruptly to 67 nT at sudden storm commencement (SSC), followed by a sharp decrease to −87 nT. A second SSC on January 24 followed by a shock on January 25 was also observed. These SSCs before the main storms and the short recovery periods imply the geomagnetic storms are CME  -driven. The sudden jump of solar wind dynamic pressure and IMF _Bz$B_z$ are also consistent with occurrence of CMEs. This is also reflected in the change in total electron content (TEC) during the storm relative to quiet days globally. The response of the ionospheric to geomagnetic storms can also be detected from wave components that account for the majority of TEC variance during the period. The dominant coherent modes of TEC variability are diurnal and semidiurnal signals which account upto 83% and 30% of the total TEC variance over fairly exclusive ionospheric regions respectively. Comparison of TEC anomalies attributed to diurnal (DW1) and semidiurnal (SW2) tides, as well as stationary planetary waves (SPW1) at 12 UTC shows enhancement in the positive anomalies following the storm. Moreover, the impact of the geomagnetic storms are distinctly marked in the daily time series of amplitudes of DW1, SW2 and SPW1. The abrupt changes in amplitudes of DW1 (5 TECU) and SW2 (2 TECU) are observed within 20°S–20°N latitude band and along 20°N respectively while that of SPW1 is about 3 TECU. Coherent oscillation with a period of 2.4 days between interplanetary magnetic field and TEC was detected during the storm. This oscillation is also detected in the amplitudes of DW1 over EIA regions in both hemispheres. Eventhough upward coupling of quasi two day wave (QTDWs) of the same periodicity, known to have caused such oscillation, are detected in both ionosphere and upper stratosphere, this one can likely be attributed to the geomagnetic storm as it happens after the storm commencement. Moreover, further analysis has indicated that QTDWs in the ionosphere are strengthened as a result of coherent oscillation of interplanetary magnetic field with the same frequency as QTDWs. On the otherhand, occurrences of minor SSW and geomagnetic storms in quick succession complicated clear demarcation of attribution of the respective events to variability of QTDWs amplitudes over upper stratosphere.     

Response of low latitude D-region ionosphere to the total solar eclipse of 22 July 2009 deduced from ELF/VLF analysis

Response of the D-region of the ionosphere to the total solar eclipse of 22 July 2009 at low latitude, Varanasi (Geog. lat., 25.27° N; Geog. long., 82.98° E; Geomag. lat. = 14° 55’ N) was investigated using ELF/VLF radio signal. Tweeks, a naturally occurring VLF signal and radio signals from various VLF navigational transmitters are first time used simultaneously to study the effect of total solar eclipse (TSE). Tweeks occurrence is a nighttime phenomena but the obscuration of solar disc during TSE in early morning leads to tweek occurrence. The changes in D-region ionospheric VLF reflection heights (h) and electron density (n_e: 22.6–24.6 cm⁻³) during eclipse have been estimated from tweek analysis. The reflection height increased from ∼89 km from the first occurrence of tweek to about ∼93 km at the totality and then decreased to ∼88 km at the end of the eclipse, suggesting significant increase in tweek reflection height of about 5.5 km during the eclipse. The reflection heights at the time of totality during TSE are found to be less by 2–3 km as compared to the usual nighttime tweek reflection heights. This is due to partial nighttime condition created by TSE. A significant increase of 3 dB in the strength of the amplitude of VLF signal of 22.2 kHz transmitted from JJI-Japan is observed around the time of the total solar eclipse (TSE) as compared to a normal day. The modeled electron density height profile of the lower ionosphere depicts linear variation in the electron density with respect to solar radiation as observed by tweek analysis also. These low latitude ionospheric perturbations on the eclipse day are discussed and compared with other normal days.     

Variability of the ionosphere over Irkutsk at low solar activity

Statistical and spectral analyses are performed to investigate variations of two ionosphere F2 layer key parameters, the critical frequency (foF2) and the peak height (hmF2), that were measured over Irkutsk (52.5°N, 104.0°E) from December 2006 to January 2008 under solar minimum. The analyses showed that both parameters contain quasi-harmonic oscillations with periods of T_n = 24/n hours (n = 1–7), among which the diurnal (n = 1) and semidiurnal (n = 2) ones are the strongest. Seasonal variations are explored of mean and median values, spectrum, amplitude, and phase of the diurnal and semidiurnal components of foF2 and hmF2.     

Comparison of topside ionosphere scale height modeled by the Topside Sounder Model and incoherent scatter radar ionospheric model

The topside ionosphere scale height extracted from two empirical models are compared in the paper. The Topside Sounder Model (TSM) provides directly the scale height (H_T), while the incoherent scatter radar ionospheric model (ISRIM) provides electron density profiles and its scale height (H_R) is determined by the lowest gradient in the topside part of the profile. H_T and H_R are presented for 7 ISR locations along with their dependences on season, local time, solar flux F10.7, and geomagnetic index ap. Comparison reveals that H_T values are systematically lower than respective H_R values as the average offset for all 7 stations is 55 km. For the midlatitude stations Arecibo, Shigaraki, and Millstone Hill this difference is reduced to 43 km. The range of variations of H_R is much larger than that of H_T, as the H_T range overlaps the lower part of the H_R range. Dependences on ap, DoY and LT are much stronger in the ISRIM than in TSM. This results in much larger values of H_R at higher ap. Diurnal amplitude of H_R is much larger than that of H_T, with large maximum of H_R at night. The present comparison yields the conclusion that the ISR measurements provide steeper topside Ne profiles than that provided by the topside sounders.     

Multi-station observation of ionospheric irregularities over South Africa during strong geomagnetic storms

This paper presents results pertaining to the response of the mid-latitude ionosphere to strong geomagnetic storms that occurred from 31 March to 02 April 2001 and 07–09 September 2002. The results are based on (i) Global Positioning Systems (GPSs) derived total electron content (TEC) variations accompanying the storm, (ii) ionosonde measurements of the ionospheric electrodynamic response towards the storms and (iii) effect of storm induced travelling ionospheric disturbances (TIDs) on GPS derived TEC. Ionospheric data comprising of ionospheric TEC obtained from GPS measurements, ionograms, solar wind data obtained from Advanced Composition Explorer (ACE) and magnetic data from ground based magnetometers were used in this study. Storm induced features in vertical TEC (VTEC) have been obtained and compared with the mean VTEC of quiet days. The response of the mid-latitude ionosphere during the two storm periods examined may be characterised in terms of increased or decreased level of VTEC, wave-like structures in VTEC perturbation and sudden enhancement in hmF2 and h′F. The study reveals both positive and negative ionospheric storm effects on the ionosphere over South Africa during the two strong storm conditions. These ionospheric features have been mainly attributed to the travelling ionospheric disturbances (TIDs) as the driving mechanism for the irregularities causing the perturbations observed. TEC perturbations due to the irregularities encountered by the satellites were observed on satellites with pseudo random numbers (PRNs) 15, 17, 18 and 23 between 17:00 and 23:00 UT on 07 September 2002.     

Sawtooth substorms generated under conditions of the steadily high solar wind energy input into the magnetosphere: Relationship between PC, AL and ASYM indices

Periodicity in occurrence of magnetic disturbances in polar cap and auroral zone under conditions of steady and powerful solar wind influence on the magnetosphere is analyzed on the example of 9 storm events with distinctly expressed sawtooth substorms (N = 48). Relationships between the polar cap magnetic activity (PC-index), magnetic disturbances in the auroral zone (AL-index) and value of the ring current asymmetry (ASYM index) were examined within the intervals of the PC growth phase and the PC decline phase inherent to each substorm. It is shown that the substorm sudden onsets are always preceded by the PC growth and that the substorm development does not affect the PC growth rate. On achieving the disturbance maximum, the PC and AL indices are simultaneously fall down to the level preceding the substorm, so that the higher the substorm intensity, the larger is the AL and PC drop in the decline phase. The ASYM index increases and decreases in conformity with the PC and AL behavior, the correlation between ASYM and PC being better than between ASYM and AL. Level of the solar wind energy input into the magnetosphere determines periodicity and intensity of disturbances: the higher the coupling function E_KL, the higher is substorm intensity and shorter is substorm length. Taking into account the permanently high level of auroral activity and inconsistency of aurora behavior and magnetic onsets during sawtooth substorms, the conclusion is made that auroral ionosphere conductivity is typically high and ensures an extremely high intensity of field-aligned currents in R1 FAC system. The periodicity of sawtooth substorms is determined by recurrent depletions and restorations of R1 currents, which are responsible for coordinated variations of magnetic activity in the polar cap and auroral zone.     

Formation mechanism of great positive TEC disturbances prior to Wenchuan earthquake on May 12, 2008

The problems of physical explanation and possible mechanisms of the seismo-ionospheric effects formation are under discussion now. There are proposed different mechanisms of such effects, for example, large- and small-scale internal gravity waves (IGWs), atmospheric electric field, electromagnetic fields and emissions. However, the appearance of local large-scale seismo-ionospheric anomalies in Total Electron Content (TEC) is possible to explain only by two mechanisms: an atmospheric electric field and/or small-scale IGWs. In this paper, the simulation results for reproduction of the observed seismo-ionospheric great positive effects in TEC prior to strong Wenchuan earthquake are presented. The obtained results confirm the proposed mechanism of seismo-ionospheric effects formation by the penetration of the seismogenic electric field from the atmosphere into the ionosphere. It is suggested that so great TEC enhancement observed 3 days prior to Wenchuan earthquake could be explained by combined action of seismogenic vertical electric field and IGWs generated by the solar terminator.     

Imaging Global Electron Content backwards in time more than 160 years ago

The Global Electron Content, GEC, represents the total number of electrons in the spherical layer over the Earth restricted by orbit of Global Positioning Satellite system (20,200 km). GEC is produced from Global Ionospheric Map of Total Electron Content, GIM-TEC, transformed to the electron density varying with height using the International Reference Ionosphere and Plasmasphere model, IRI-Plas. The climatologic GEC model is developed from GIM-TEC maps for a period 1999–2012 including the solar activity, annual and semiannual cycles as the most important factors affecting daily GEC variation. The proxy R_zp of the international sunspot numbers, R_i, is used as a measure of solar activity composed of 3 day smoothed R_i, 7 day and 81 day backwards mean of R_i scaled to the range of 1–40 proxy units, p.u. The root mean square error of the GEC climatologic model is found to vary from 8% to 13% of GEC. Taking advantage of a long history of sunspot numbers, the climatologic GEC model is applied for GEC reconstruction backwards in time for more than 160 years ago since 1850. The extended set of GEC values provides the numerical representation of the ionosphere and plasmasphere electron content coherent with variations of solar activity as a potential proxy index driving the ionosphere models.     

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).