Dependence of the ionospheric response on the solar flare parameters based on the theoretical modeling and GPS data

The measurements of an increase in the total electron content (TEC) of the ionosphere during solar flares, obtained based on the GPS data, indicated that up to 30% of TEC increments corresponded to the ionospheric regions above 300 km altitude in some cases, and TEC increased mainly below altitudes of 300 km in other cases. The theoretical model of the ionosphere and plasmasphere was used to study the obtained effects. The altitude-time variations in the charged particle density in the ionospheric region from 100 to 1000 km were used depending on the solar flare spectrum. An analysis of the modeling results indicated that an intensification of the flare UV emission in the 55–65 and 85–95 nm spectral ranges results in a pronounced increase in the electron density in the topside ionosphere (above 300 km). The experimental dependences of the ionospheric TEC response amplitude on the localization and peak power of flares on the Sun in the X-ray range, obtained based on the GPS data, are also presented in the work.     

Variability of equatorial ionospheric electron density at fixed heights below the F2 peak

A study on variability of the equatorial ionosphere was carried out at fixed heights below the F2 peak for two different levels of solar activity. The study covered height range of 100 km up to the peak of F2 layer using a real height step increase of 10 km. The variability index used is the percentage ratio of standard deviation over the average value for the month. Daytime minimum variability of between 3% and 10% was observed at height range of about 150–210 km during low solar activity and between 2% and 7% at height range of 160–220 km during high solar activity. The nighttime maximum of between 70% and 187% was observed at height range of about 210–250 km during low solar activity and between 42% and 127% at height range of 210–250 km during high solar activity. The height range at which daytime minimum was observed falls within the F1 height of the ionosphere. The result obtained is consistent with previous works carried out in the low latitude locations for American sector.     

Prompt derivation of TEC from GEONET data for space weather monitoring in Japan

Continuous monitoring of ionospheric conditions is essential to monitoring and forecasting space weather. The worldwide use of global navigation satellite systems like the Gobal Positioning System (GPS) makes it possible to continuously monitor the total electron content (TEC) of the ionosphere and plasmasphere up to a height of about 20,000 km. We have developed a system for deriving the TEC from GEONET data rapidly and we use the TEC distribution over Japan in the daily operations of the Space Weather Forecast Center at NICT (RWC Tokyo of ISES). Using instrumental biases from a few days before enables us to drastically shorten the processing time for deriving TEC. The latest TEC values (with a delay of about 1 h) are obtained every 3 h, and most of the values are within 2 TEC units of the actual TEC. We have found our system for deriving TEC rapidly to be useful for continuously monitoring the progress of ionospheric storms under any ionospheric conditions, even those under which the usual ionosonde observations are unable to obtain F-region profiles.     

Ionospheric and geomagnetic disturbances during the 2005 Sumatran earthquakes

This paper investigates the ionospheric and geomagnetic responses during the 28 March 2005 and 14 May 2005 Sumatran earthquakes using GPS and magnetometer stations located in the near zone of the epicenters. These events occurred during low solar and geomagnetic activity. TEC oscillations with periods of 5–10 min were observed about 10–24 min after the earthquakes and have horizontal propagation velocities of 922–1259 m/s. Ionospheric disturbances were observed at GPS stations located to the northeast of the epicenters, while no significant disturbances were seen relatively east and south of the epicenters. The magnetic field measurements show rapid fluctuations of 4–5 s shortly after the earthquake, followed by a Pc5 pulsation of 4.8 min about 11 min after the event. The correlation between the ionospheric and geomagnetic responses shows a good agreement in the period and time lag of the peak disturbance arrival, i.e. about 11–13 min after the earthquake.     

Unusual sudden ionospheric disturbance from solar flare of 4 November 2003

An interpretation of the nature of the sudden ionospheric disturbance in terms of response to X-ray flux enhancement in the band 1–20 Å has been made by many authors. Last decades investigations revealed presence of important qualitative distinctions in spatially temporal pattern of geomagnetic response to solar flares featuring harder radiation spectra (with quanta energies above 100 keV). These distinctions can not be adequately described by classical theory implying ionization growing on E and D ionosperic layers and intensification of Sq-current system. In this respect, solar flare on 4 November 2003 characterizing by existence of two separate (time lag ~45 min) spectral maximums in X-rays range (average quantum energy <100 keV) and in γ-rays range (average quantum energy >100 keV), represents convenient proving ground for study of specifics the geomagnetic response to bursts marked by different hardness. In current article, we show that this flare has a number of unusual features including specific variation of accompanying current system and magnetospheric manifestation that is observed in trapped radiation fluxes and magnetic field on geosynchronous orbit. Possible physical mechanism leading to intensification of magnetospheric–ionospheric current system is discussed.     

Climatology of total electron content near the dip equator under geomagnetic quiet-conditions

This study analyzes the TEC data during 1998–2007, observed by the AREQ (16.5°S, 71.5°W) GPS station to investigate the equatorial ionospheric variations under geomagnetic quiet-conditions. The diurnal TEC values generally have a maximum value between 1330 and 1500 LT and a minimum around 0500 LT. For the seasonal variation, the semi-annual variation apparently exists in the daytime TEC with two peaks in equinoctial months. In contrast, this semi-annual variation is not found in the nighttime. Furthermore, the results of the annual variation show that the correlation between the daytime TEC value and the solar activity factor is highly positive.     

Ionospheric activity and possible connection with seismicity: Contribution from the analysis of long time series of GNSS signals

The modifications of some atmospheric physical properties prior to a high magnitude earthquake were debated in the frame of the Lithosphere Atmosphere Ionosphere Coupling (LAIC) model. In this work, among the variety of involved phenomena, the ionisation of air at the ionospheric levels triggered by the leaking of gases from the Earth’s crust was investigated through the analysis of GNSS (Global Navigation Satellite System) signals. In particular, the authors analysed a 5 year (2008–2012) long series of GNSS based ionospheric TEC to produce maps over an area surrounding the epicentre of the L’Aquila (Italy, M_w = 6.3) earthquake of April 6th, 2009. The series was used to detect and quantify amplitude and duration of episodes of ionospheric disturbances by a statistical approach and to discriminate local and global effects on the ionosphere comparing these series with TEC values provided by the analysis of GNSS data from international permanent trackers distributed over a wider region. The study found that during this time interval only three statistically meaningful episodes of ionospheric disturbances were observed. One of them, occurring during the night of 16th of March 2009, anticipated the main shock by 3 weeks and could be connected with the strong earthquake of 6th of April. The other two significant episodes were detected within periods that were not close to the main seismic events and are more likely due to various and global reasons.     

Ionospheric slab thickness and its seasonal variations observed by GPS

The ionospheric slab thickness, the ratio of the total electron content (TEC) to the F2-layer peak electron density (NmF2), is closely related to the shape of the ionospheric electron density profile Ne (h) and the TEC. Therefore, the ionospheric slab thickness is a significant parameter representative of the ionosphere. In this paper, the continuous GPS observations in South Korea are firstly used to study the equivalent slab thickness (EST) and its seasonal variability. The averaged diurnal medians of December–January–February (DJF), March–April–May (MAM), June–July–August (JJA) and September–October–November (SON) in 2003 have been considered to represent the winter, spring, summer and autumn seasons, respectively. The results show that the systematic diurnal changes of TEC, NmF2 and EST significantly appeared in each season and the higher values of TEC and NmF2 are observed during the equinoxes (semiannual anomaly) as well as in the mid-daytime of each season. The EST is significantly smaller in winter than in summer, but with a consistent variation pattern. During 14–16 LT in daytime, the larger EST values are observed in spring and autumn, while the smaller ones are in summer and winter. The peaks of EST diurnal variation are around 10–18 LT which are probably caused by the action of the thermospheric wind and the plasmapheric flow into the F2-region.     

Variations in ionospheric parameters during solar flares

Results of studying the ionospheric response to solar flares, obtained based on the incoherent scatter radar observations of the GPS signals and as a result of the model simulations, are presented. The method, based on the effect of partial “shadowing” of the atmosphere by the globe, has been used to analyze the GPS data. This method made it possible to estimate the value of a change in the electron content in the upper ionosphere during the solar flare of July 14, 2000. It has been shown that a flare can cause a decrease in the electron content at heights of the upper ionosphere (h > 300 km) according to the GPS data. Similar effects in the formation of a negative disturbance in the ionospheric F region were also observed during the solar flares of May 21 and 23, 1967, at the Arecibo incoherent scatter radar. The mechanism by which negative disturbances are formed in the upper ionosphere during solar flares has been studied based on the theoretical model of the ionosphere-plasmasphere coupling. It has been shown that an intense ejection of O⁺ ions into the above located plasmasphere under the action of a sharp increase in the ion production rate and the thermal expansion of the ionospheric plasma cause the formation of a negative disturbance in the electron concentration in the upper ionosphere.     

Dynamics of North American sector ionospheric and thermospheric response during the November 2004 superstorm

We present a study of ionospheric and thermospheric response during a November 9–10, 2004 major geomagnetic storm event (DsT ～?300 nT). We utilize the North American sector longitude chain of incoherent scatter radars at Arecibo, Millstone Hill, and Sondrestrom, operating as part of a coordinated international mesosphere/lower thermosphere coupling study experiment. Total electron content (TEC) determinations from global positioning system (GPS) ground receivers, ground magnetometer traces from the Canadian CANOPUS array, Defense Meteorological Satellite Platform (DMSP) topside data, and global convection patterns from the SuperDARN radar network are analyzed to place the detailed radar data in proper mesoscale context. The plasmaspheric boundary layer (PBL) expanded greatly in the dusk sector during ring current intensification to span more than 25° of magnetic latitude, reaching as far south as 30° invariant latitude. Strong sub-auroral polarization stream velocities of more than 1 km/s were accompanied by large upwards thermal O+ fluxes to the overlying magnetosphere. The large PBL expansion subsequently exposed both Millstone Hill and Sondrestrom to the auroral convection pattern, which developed a complex multicell and reverse convection response under strongly northward IMF conditions during a period of global interplanetary electric field penetration. Large traveling atmospheric and ionospheric disturbances caused significant neutral wind and ion velocity surges in the mid-latitude and tropical ionosphere and thermosphere, with substorm activity launching equatorward neutral wind enhancements and subsequent mid-latitude dynamo responses at Millstone Hill. However, ionosphere and thermosphere observations at Arecibo point to significant disturbance propagation modification in the post-dusk sector PBL region.     

Variability of total electron content near the crest of the equatorial anomaly during moderate geomagnetic storms

The ionospheric responses to a large number (116) of moderate (?50≥D_st>?100 nT) geomagnetic storms distributed over the period (1980–1990) are investigated using total electron content (TEC) data recorded at Calcutta (88.38°E, 22.58°N geographic, dip: 32°N). TEC perturbations exhibit a prominent dependence on the local times of main phase occurrence (MPO). The storms with MPO during daytime hours are more effective in producing larger deviations and smaller time delays for maximum positive deviations compared to those with nighttime MPO. Though the perturbations in the equinoctial and winter solstitial months more or less follow the reported climatology, remarkable deviations are detected for the summer solstitial storms. Depending on the local times of MPO, the sunrise enhancement in TEC is greatly perturbed. The TEC variability patterns are interpreted in terms of the storm time modifications of equatorial electric field, wind system and neutral composition.     

Disturbances in the topside ionosphere during solar flares

The results of studying the ionospheric response to solar flares, obtained from the data of the GPS signal observations and incoherent scatter radars and as a result of the model calculations, are presented. It is shown that, according to the GPS data, a flare can cause a decrease in the electron content at altitudes of the topside ionosphere (h > 300 km). Similar effects of formation of a negative disturbance in the ionospheric F region were also observed during the solar flares of May 21 and 23, 1967, with the Arecibo incoherent scatter radar. The mechanism by which negative disturbances appear in the topside ionosphere during solar flares has been studied in this work based on the theoretical model of the ionosphere-plasmasphere coupling. It has been indicated that the formation of the electron density negative disturbance in the topside ionosphere is caused by an intense removal of O⁺ ions into the overlying plasmasphere under the action of an abrupt increase in the ion production rate and thermal expansion of the ionospheric plasma.     

Detection of large scale TIDs associated with the dayside cusp using SuperDARN data

Variations in the dayside ionosphere parameters as a result of a large-scale acoustic gravity wave (LS AGW) were studied for the 17 February 1998 substorm using the super dual auroral radar network (SuperDARN) measurements. This event was characterised by a sharp rise in the AE index with a maximum of ～900 nT. The source of the disturbance responsible for the LS AGW appears to have been located within the plasma convection throat and in the dayside cusp region. The location of the source was obtained from studies of a number of datasets including high-latitude convection maps, data from 4 DMSP satellites and networks of ground-based magnetometers. The propagation of the LS AGWs caused quasi-periodic variations in the skip distance (with an amplitude up to 220–260 km) of the ground backscatter measured by up to 6 SuperDARN radars, including Goose Bay and Kapuskasing, resulting in two large-scale travelling ionospheric disturbances (LS TIDs). The LS TIDs had wave periods of 1.5 and 2 h, a velocity of ～400 m/s for both, and wavelengths of 2200 and 2900 km, respectively. These quasi-periodic variations were also present in the peak electron density and height of the F2 layer measured by the Goose Bay ionosonde. The numerical simulation of the inverse problem show good agreement between Goose Bay radar and Goose Bay ionosonde measurements. But these LS TIDs would be difficult to deduce from the ground based ionospheric station data alone, because h_mF2 variations were 10–40 km only and f_OF2 variations between 10% and 20%. The results demonstrate how important SuperDARN radars can be, and that this is a more powerful technique than routine ground-based sounding for studies of weak quasi-periodic variations in the dayside subauroral ionosphere related to LS AGW.     

TEC variations over the Mediterranean before and during the strong earthquake (M = 6.5) of 12th October 2013 in Crete,Greece

In this paper, the total electron content (TEC) data from eight global positioning system (GPS) stations of the EUREF network, provided by IONOLAB (Turkey), were analyzed using discrete Fourier analysis to investigate the TEC variations over the Mediterranean before and during the strong earthquake of 12th October 2013, which occurred west of Crete, Greece. In accordance with the results of similar analyses in the area, the main conclusions of this study are the following: (a) TEC oscillations in a broad range of frequencies occur randomly over an area of several hundred km from the earthquake and (b) high frequency oscillations (f ⩾ 0.0003 Hz, periods T ⩽ 60 m) may point to the location of the earthquake with questionable accuracy. The fractal characteristics of the frequency distribution may point to the locus of the earthquake with higher accuracy. We conclude that the lithosphere–atmosphere–ionosphere coupling (LAIC) mechanism through acoustic or gravity waves could explain this phenomenology.     

Asymmetric response of the topside ionosphere to large-scale IGW generated during the November 30, 1979, substorm

We used bottomside ground observations and topside sounding data from the Intercosmos-19 satellite to study a Travelling Ionospheric Disturbance (TID) that occurred in response to Large-Scale Internal Gravity Wave (LSIGW) propagation during a substorm on November 30, 1979. We built a global scheme for the wavelike ionospheric variations during this medium substorm (AE_max ～800 nT). The area where the TID was observed looks like a wedge since it covers the nighttime hours at subauroral latitudes but contracts to a ～02 h local sector at low latitudes. The ionospheric response is strongly asymmetric because the wedge area and the TID amplitude are larger in the winter hemisphere than in the summer hemisphere. Clear evidence was obtained indicating that the more powerful TID from the Northern (winter) hemisphere propagated across the equator into the low latitude Southern (summer) hemisphere. Intercosmos-19 observations show that the disturbance covers the entire thickness of the topside ionosphere, from h_mF2 up to at least the 1000 km satellite altitude at post-midnight local times. F-layer lifting reached ～200 km, N_e increases in the topside ionosphere by up to a factor of ～1.9 and variations in N_mF2 of both signs were observed. Assumptions are made concerning the reason for the IGW effect at high altitudes in the topside ionosphere. The relationship between TID parameters and source characteristics determined from a global network of magnetometers are studied. The role of the dayside cusp in the generation of the TID in the daytime ionosphere is discussed. The magnetospheric electric field effects are distinguished from IGW effects.     

On the sensitivity of total electron content (TEC) to upper atmospheric/ionospheric parameters

The total electron content (TEC) is a key ionospheric parameter for various space weather applications. Over the last decade an extensive database of TEC measurements has become available from both space- and ground-based observations, and these measurements have established the general morphology of the global TEC distributions. In particular, the TOPEX TEC measurements have shown strong longitudinal variations of TEC in addition to the observed day-to-day variabilities. To better understand the observed TEC variations and to better guide its modeling, we have studied the sensitivity of quiet-time TEC to the following key atmospheric and ionospheric parameters: neutral density, neutral wind, plasma temperatures, plasmaspheric flux, and the O⁺–O collision frequency. These parameters are often only roughly known and can cause large uncertainties in model results. For this study, we have developed a numerical mid-latitude ionospheric model, which solves the momentum and continuity equations for the O⁺ density and a simplified set of equations for the H⁺ density. To obtain TEC, the calculated ion densities have been integrated from the bottom altitude (100 km) to the altitude of the TOPEX satellite (1336 km). Our study shows that during the day the neutral wind and the neutral composition have the most important effect on TEC. In particular, the zonal component of the neutral wind can have a large effect on TEC in the southern hemisphere where the magnetic declination angle is large. During the night, most of the above-mentioned parameters can play a significant role in the TEC morphology, except for the plasma temperature, which has only a small effect on TEC. Finally, the TEC varies roughly linearly with respect to all of the parameters except for the neutral wind.     

On the connection between the atmospheric electric field measured at the surface and the ionospheric electric field in the Central Antarctica

Regular measurements of the atmospheric electric field made at Vostok Station (φ=78.45°S; λ=106.87°E, elevation 3500 m) in Antarctica demonstrate that extremely intense electric fields (1000–5000 V/m) can be observed during snow storms. Usually the measured value of the atmospheric electric field at Vostok is about 100–250 V/m during periods with “fair weather” conditions. Actual relation between near-surface electric fields and ionospheric electric fields remain to be a controversial problem. Some people claimed that these intense electric fields produced by snowstorms or appearing before strong earthquakes can re-distribute electric potential in the ionosphere at the heights up to 300 km. We investigated interrelation between the atmospheric and ionospheric electric fields by both experimental and theoretical methods. Our conclusion is that increased near-surface atmospheric electric fields do not contribute notably to distribution of ionospheric electric potential.     

First incoherent scatter radar observations of ionospheric heating on the second electron gyro-harmonic

We report first results from a unique experiment performed at the HIPAS ionospheric modification facility in conjunction with the Poker Flat incoherent scatter radar in Alaska. High-power radio waves at 2.85 MHz, which corresponds to the second electron gyro-harmonic at ~245 km altitude, were transmitted into the nighttime ionosphere. Clear evidence of F-region ionospheric electron temperature enhancements were found, for the first time at this pump frequency, maximizing when the pump frequency is close to the second gyro-harmonic and double resonance. This is consistent with previous pump-enhanced artificial optical observations. We estimate the plasma heating efficiency to be approximately double that for higher pump frequencies.     

The polar wind: Recent observations

The polar wind is an ambipolar outflow of thermal plasma from the high-latitude ionosphere to the magnetosphere, and it primarily consists of H⁺, He⁺ and O⁺ ions and electrons. Statistical and episodic studies based primarily on ion composition observations on the ISIS-2, DE-1, Akebono and Polar satellites over the past four decades have confirmed the existence of the polar wind. These observations spanned the altitude range from 1000 to ∼50,500 km, and revealed several important features in the polar wind that are unexpected from “classical” polar wind theories. These include the day–night asymmetry in polar wind velocity, which is 1.5–2.0 times larger on the dayside; appreciable O⁺ flow at high altitudes, where the velocity at 5000–10,000 km is of 1–4 km/s; and significant electron temperature anisotropy in the sunlit polar wind, in which the upward-to-downward electron temperature ratio is 1.5–2. These features are attributable to a number of “non-classical” polar wind ion acceleration mechanisms resulting from strong ionospheric convection, enhanced electron and ion temperatures, and escaping atmospheric photoelectrons. The observed polar wind has an averaged ion temperature of ∼0.2–0.3 eV, and a rate of ion velocity increase with altitude that correlates strongly with electron temperature and is greatest at low altitudes (<4000 km for H⁺). The rate of velocity increase below 4000 km is larger at solar minimum than at solar maximum. Above 4000 km, the reverse is the case. This suggests that the dominant polar wind ion acceleration process may be different at low and high altitudes, respectively. At a given altitude, the polar wind velocity is highly variable, and is on average largest for H⁺ and smallest for O⁺. Near solar maximum, H⁺, He⁺, and O⁺ ions typically reach a velocity of 1 km/s near 2000, 3000, and 6000 km, respectively, and velocities of 12, 7, and 4 km/s, respectively, at 10,000 km altitude. Near solar minimum, the velocity of all three species is smaller at high altitudes. Observationally it is not always possible to unambiguously separate an energized “non-polar-wind” ion such as a low-energy “cleft ion fountain” ion that has convected into a polar wind flux tube from an energized “polar-wind” ion that is accelerated locally by “non-classical” polar-wind ion acceleration mechanisms. Significant questions remain on the relative contribution between the cleft ion fountain, auroral bulk upflow, and the topside polar-cap ionosphere to the O⁺ polar wind population at high altitudes, the effect of positive spacecraft charging on the lowest-energy component of the H⁺ polar wind population, and the relative importance of the various classical and non-classical ion acceleration mechanisms. These questions pose several challenges in future polar wind observations: These include measurement of the lowest-energy component in the presence of positive spacecraft potential, definitive determination and if possible active control of the spacecraft potential, definitive discrimination between polar wind and other inter-mixed thermal ion populations, measurement of the three-dimensional ion drift velocity vector and the parallel and perpendicular ion temperatures or the detailed three-dimensional velocity distribution function, and resolution of He⁺ and other minor ion species in the polar wind population.     

Estimation of the high-latitude topside electron heat flux using DMSP plasma density measurements

The high-latitude ionosphere interfaces with the hot, tenuous, magnetospheric plasma, and a heat flow into the ionosphere is expected, which has a large impact on the plasma densities and temperatures in the high-latitude ionosphere. The value of this magnetospheric heat flux is unknown. In an effort to estimate the value of the magnetospheric heat flux into the high-latitude ionosphere, and show its effect on the high-latitude ionospheric plasma densities, we ran an ensemble of model runs using the Ionosphere Forecast Model (IFM) with different values of the heat flux through the upper boundary. These model runs included heating from both auroral and solar sources. Then, for each heat flux value, the plasma densities obtained from the model runs, at 840 km, were compared to the corresponding values measured by the DMSP F13 satellite. The heat flux value that gave the best comparison between the measured and calculated plasma densities was considered to be the best estimate for the topside heat flux. The comparison was conducted for a 1-year data set of the DMSP F13 measured plasma densities (4300 consecutive orbits). Our systematic IFM/DMSP plasma density comparisons indicate that when a zero magnetospheric downward heat flux is assumed at the upper boundary of the IFM model, on the average, the IFM underestimates the measured plasma densities by a factor of 2. A good IFM/DMSP plasma density comparison was achieved for each month in 1998 when for each month a constant heat flux was assumed at the upper boundary of the model. For the 12-month period, the heat flux values that gave the best IFM/DMSP plasma density comparisons varied on the average from −0.5×10¹⁰ to −1.5×10¹⁰ eV cm⁻² s⁻¹.