The global response of terrestrial ionosphere to the December 2015 space weather event

This paper investigates the ionospheric storm of December 19–21, 2015, which was initiated by two successive CME eruptions that caused a G3 space weather event. We used the in situ electron density (N_e) and electron temperature (T_e) and the Total Electron Content (TEC) measurements from SWARM-A satellite, as well as the O/N₂ observations from TIMED/GUVI to study the ionospheric impact. The observations reveal the longitudinal and hemispherical differences in the ionospheric response to the storm event. A positive ionospheric storm was observed over the American, African and Asian regions on 20 December, and the next day showed a negative storm. Both these exhibited hemispheric differences. A positive storm was observed over the East Pacific region on 21 December. It is seen that the net effect of both the disturbance dynamo electric field and composition differences become important in explaining the observed variability in topside ionospheric densities. In addition, we also discuss the T_e variations that occurred as a consequence of the space weather event.     

Large-scale traveling ionospheric disturbances of auroral origin according to the data of the GPS network and ionosondes

The intensity of large-scale traveling ionospheric disturbances (LS TIDs), registered using measurements of total electron content (TEC) during the magnetic storms on October 29–31, 2003, and on November 7–11, 2004, had been compared with that of local electron density disturbances. The data of TEC measurements at ground-based GPS receivers located near the ionospheric stations and the corresponding values of the critical frequency of the ionospheric F region f_oF2 were used for this purpose. The variations of TEC and f_oF2 were similar for all events mentioned above. The previous assumption that the ionospheric region with vertical extension from 150 to 200 km located near the F-layer maximum mainly contributes to the TEC variations was confirmed for the cases when the electron density disturbance at the F region maximum was not more than 50%. However, this region probably becomes vertically more extended when the electron density disturbance in the ionospheric F region is about 85%.     

Empirical model of the TEC response to the geomagnetic activity over the North American region

The paper presents an empirical model of the total electron content (TEC) response to the geomagnetic activity described by the K_p-index. The model is built on the basis of TEC measurements covering the region of North America (50°W–150°W, 10°N–60°N) for the period of time between October 2004 and December 2009. By using a 2D (latitude-time) cross-correlation analysis it is found that the ionospheric response to the geomagnetic activity over the considered geographic region and at low solar activity revealed both positive and negative phases of response. The both phases of the ionospheric response have different duration and time delay with respect to the geomagnetic storm. It was found that these two parameters of the ionospheric response depend on the season and geographical latitude. The presence of two phases, positive and negative, of the ionospheric response imposed the implementation of two different time delay constants in order to properly describe the two different delayed reactions. The seasonal dependence of the TEC response to geomagnetic storms is characterized by predominantly positive response in winter with a short (usually ∼5–6 h) time delay as well as mainly negative response in summer with a long (larger than 15 h) time delay. While the TEC response in March and October is more close to the winter one the response in April and September is similar to the summer one.     

Ionospheric response to solar and magnetospheric protons during January 15–22, 2005: EAGLE whole atmosphere model results

We present an analysis of the ionosphere and thermosphere response to Solar Proton Events (SPE) and magnetospheric proton precipitation in January 2005, which was carried out using the model of the entire atmosphere EAGLE. The ionization rates for the considered period were acquired from the AIMOS (Atmospheric Ionization Module Osnabrück) dataset. For numerical experiments, we applied only the proton-induced ionization rates of that period, while all the other model input parameters, including the electron precipitations, corresponded to the quiet conditions. In January 2005, two major solar proton events with different energy spectra and proton fluxes occurred on January 17 and January 20. Since two geomagnetic storms and several sub-storms took place during the considered period, not only solar protons but also less energetic magnetospheric protons contributed to the calculated ionization rates. Despite the relative transparency of the thermosphere for high-energy protons, an ionospheric response to the SPE and proton precipitation from the magnetotail was obtained in numerical experiments. In the ionospheric E layer, the maximum increase in the electron concentration is localized at high latitudes, and at heights of the ionospheric F2 layer, the positive perturbations were formed in the near-equatorial region. An analysis of the model-derived results showed that changes in the ionospheric F2 layer were caused by a change in the neutral composition of the thermosphere. We found that in the recovery phase after both solar proton events and the enhancement of magnetospheric proton precipitations associated with geomagnetic disturbances, the TEC and electron density in the F region and in topside ionosphere/plasmasphere increase at low- and mid-latitudes due to an enhancement of atomic oxygen concentration. Our results demonstrate an important role of magnetospheric protons in the formation of negative F-region ionospheric storms. According to our results, the topside ionosphere/plasmasphere and bottom-side ionosphere can react to solar and magnetospheric protons both with the same sign of disturbances or in different way. The same statement is true for TEC and foF2 disturbances. Different disturbances of foF2 and TEC at high and low latitudes can be explained by topside electron temperature disturbances.     

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.     

The COMPLEIK subroutine of the IONORT-ISP system for calculating the non-deviative absorption: A comparison with the ICEPAC formula

The present paper proposes to discuss the ionospheric absorption, assuming a quasi-flat layered ionospheric medium, with small horizontal gradients. A recent complex eikonal model (Settimi et al., 2013b) is applied, useful to calculate the absorption due to the ionospheric D-layer, which can be approximately characterized by a linearized analytical profile of complex refractive index, covering a short range of heights between h₁ = 50 km and h₂ = 90 km. Moreover, Settimi et al. (2013c) have already compared the complex eikonal model for the D-layer with the analytical Chapman’s profile of ionospheric electron density; the corresponding absorption coefficient is more accurate than Rawer’s theory (1976) in the range of middle critical frequencies. Finally, in this paper, the simple complex eikonal equations, in quasi-longitudinal (QL) approximation, for calculating the non-deviative absorption coefficient due to the propagation across the D-layer are encoded into a so called COMPLEIK (COMPLex EIKonal) subroutine of the IONORT (IONOspheric Ray-Tracing) program ( Azzarone et al., 2012). The IONORT program, which simulates the three-dimensional (3-D) ray-tracing for high frequencies (HF) waves in the ionosphere, runs on the assimilative ISP (IRI-SIRMUP-P) discrete model over the Mediterranean area ( Pezzopane et al., 2011). As main outcome of the paper, the simple COMPLEIK algorithm is compared to the more elaborate semi-empirical ICEPAC formula (Stewart, undated), which refers to various phenomenological parameters such as the critical frequency of E-layer. COMPLEIK is reliable just like the ICEPAC, with the advantage of being implemented more directly. Indeed, the complex eikonal model depends just on some parameters of the electron density profile, which are numerically calculable, such as the maximum height.     

Improvement of the real time total electron content based on the International Reference Ionosphere model

One of the methods to adapt the International Reference Ionosphere model to real time conditions is to use instantaneous values of the critical frequency of the ionosphere foF2. It is shown that there can be large discrepancies between model and experimental values of the total electron content TEC after this adaptation. Improvement can be provided by the use of an empirical model of the ionospheric slab thickness τ. This conclusion is based on analysis of contributions of various correction factors connected with foF2 and τ into discrepancies between model and experimental values of TEC.     

Solar activity variations of ionosonde measurements and modeling results

The time series of hourly electron density profiles N(h) obtained at several mid-latitude stations in Europe have been used to obtain N(h) profiles on a monthly basis and to extract both the expected bottomside parameters and a proxy of the ionospheric variability as functions of time and height. With these data we present advances on a “Local Model” technique for the parameters B0 and B1, its applicability to other ionospheric stations, to other bottomside ionospheric parameters, and to modeling the time/height variability of the profile. The Local Model (LM) is an empirical model based on the experimental results of the solar activity dependence of the daily and seasonal behavior of the above parameters. The LM improves the IRI-2001 prediction of the B0 and B1 by factor of two at mid-latitudes. Moreover, the LM can be used to simulate other ionospheric parameters and to build mean N(h) profiles and the deviations from them. The modeling of both the average N(h) profiles and their deviations is an useful tool for ionospheric model users who want to know both the expected patterns and their deviations.     

A new inversion algorithm for HF sky-wave backscatter ionograms

HF sky-wave backscatter sounding system is capable of measuring the large-scale, two-dimensional (2-D) distributions of ionospheric electron density. The leading edge (LE) of a backscatter ionogram (BSI) is widely used for ionospheric inversion since it is hardly affected by any factors other than ionospheric electron density. Traditional BSI inversion methods have failed to distinguish LEs associated with different ionospheric layers, and simply utilize the minimum group path of each operating frequency, which generally corresponds to the LE associated with the F₂ layer. Consequently, while the inversion results can provide accurate profiles of the F region below the F₂ peak, the diagnostics may not be so effective for other ionospheric layers. In order to resolve this issue, we present a new BSI inversion method using LEs associated with different layers, which can further improve the accuracy of electron density distribution, especially the profile of the ionospheric layers below the F₂ region. The efficiency of the algorithm is evaluated by computing the mean and the standard deviation of the differences between inverted parameter values and true values obtained from both vertical and oblique incidence sounding. Test results clearly manifest that the method we have developed outputs more accurate electron density profiles due to improvements to acquire the profiles of the layers below the F₂ region. Our study can further improve the current BSI inversion methods on the reconstruction of 2-D electron density distribution in a vertical plane aligned with the direction of sounding.     

The variability and predictability of the IRI B0, B1 parameters over Grahamstown,South Africa

The International Reference Ionosphere (IRI) parameters B₀ and B₁ provide a representation of the thickness and shape, respectively, of the F2 layer of the bottomside ionosphere. These parameters can be derived from electron density profiles that are determined from vertical incidence ionograms. This paper aims to illustrate the variability of these parameters for a single mid latitude station and demonstrate the ability of the Neural Network (NN) modeling technique for developing a predictive model for these parameters. Grahamstown, South Africa (33.3°S, 26.5°E) was chosen as the mid latitude station used in this study and the B₀ and B₁ parameters for an 11 year period were determined from electron density profiles recorded at that station with a University of Massachusetts Lowell Center for Atmospheric Research (UMLCAR) Digisonde. A preliminary single station NN model was then developed using the Grahamstown data from 1996 to 2005 as a training database, and input parameters known to affect the behaviour of the F2 layer, such as day number, hour, solar and magnetic indices. An analysis of the diurnal, seasonal and solar variations of these parameters was undertaken for the years 2000, 2005 and 2006 using hourly monthly median values. Comparisons between the values derived from measured data and those predicted using the two available IRI-2001 methods (IRI tables and Gulyaeva, T. Progress in ionospheric informatics based on electron density profile analysis of ionograms. Adv. Space Res. 7(6), 39–48, 1987.) and the newly developed NN model are also shown in this paper. The preliminary NN model showed that it is feasible to use the NN technique to develop a prediction tool for the IRI thickness and shape parameters and first results from this model reveal that for the mid latitude location used in this study the NN model provides a more accurate prediction than the current IRI model options.     

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.     

Parameterized Regional Ionospheric Model and a comparison of its results with experimental data and IRI representations

We describe a Parameterized Regional Ionospheric Model (PARIM) to calculate the spatial and temporal variations of the ionospheric electron density/plasma frequency over the Brazilian sector. The ionospheric plasma frequency values as calculated from an enhanced Sheffield University Plasmasphere–Ionosphere Model (SUPIM) were used to construct the model. PARIM is a time-independent 3D regional model (altitude, longitude/local time, latitude) used to reproduce SUPIM plasma frequencies for geomagnetic quiet condition, for any day of the year and for low to moderately high solar activity. The procedure to obtain the modeled representation uses finite Fourier series so that all plasma frequency dependencies can be represented by Fourier coefficients. PARIM presents very good results, except for the F region peak height (hmF2) near the geomagnetic equator during times of occurrence of the F₃ layer. The plasma frequency calculated by IRI from E region to bottomside of the F region present latitudinal discontinuities during morning and evening times for both solar minimum and solar maximum conditions. Both the results of PARIM and the IRI for the E region peak density show excellent agreement with the observational values obtained during the conjugate point equatorial experiment (COPEX) campaign. The IRI representations significantly underestimate the foF2 and hmF2 compared to the observational results over the COPEX sites, mainly during the evening–nighttime period.     

Thermosphere–ionosphere response to a severe magnetic storm: A case study

This paper reports the response of the ionosphere–thermosphere system to an intense geomagnetic storm. For that, data taken by instruments on board Dynamic Explorer 2 at heights of the F2-layer (molecular nitrogen N₂ and atomic oxygen O compositions, neutral temperature Tg and electron density Ne) were used. The ionospheric response is characterized by a negative storm effect expanding from mid–high to low latitude. It is observed during this severe geomagnetic storm that negative effects were caused mainly by an increase in molecular nitrogen composition N₂ and almost no changes in atomic oxygen composition O.     

An empirical model of electron density in low latitude at 600 km obtained by Hinotori satellite

An empirical model of electron density (Ne) was constructed by using the data obtained with an impedance probe on board Japanese Hinotori satellite. The satellite was in circular orbit of the height of 600 km with the inclination of 31 degrees from February 1981 to June 1982. The constructed model gives Ne at any local time with the time resolution of 90 min and between −25 and 25 degrees in magnetic latitude with its resolution of 5 degrees in the range of F_10.7 from 150 to 250 under the condition of K_p < 4. Spline interpolations are applied to the functions of day of year, geomagnetic latitude and solar local time, and linear interpolation is applied to the function of F_10.7. Longitude dependence of Ne is not taken into account. Our density model can reproduce solar local time variation of electron density at 600 km altitude better than current International Reference Ionosphere (IRI2001) model which overestimates Ne in night time and underestimates Ne in day time. Our density model together with electron temperature model which has been constructed before will enable more understanding of upper ionospheric phenomenon in the equatorial region.     

Using the radial basis function neural network to predict ionospheric critical frequency of F2 layer over Wuhan

Neural networks (NNs) have been applied to ionospheric predictions recently. This paper uses radial basis function neural network (RBF-NN) to forecast hourly values of the ionospheric F₂ layer critical frequency(f_oF₂), over Wuhan (30.5N, 114.3E), China. The false nearest neighbor method is used to determine the embedding dimension, and the principal component analysis (PCA) is used to reduce noise and dimension. The whole study is based on a sample of about 26,000 observations of f_oF₂ with 1-h time resolution, derived during the period from January 1981 to December 1983. The performance of RBF-NN is estimated by calculating the normalized root-mean-squared (NRMSE) error, and its results show that short-term predictions of f_oF₂ are improved.     

Evaluation of the STORM model storm-time corrections for middle latitude

This paper presents results from the Storm-Time Ionospheric Correction Model (STORM) validation for selected Northern and Southern Hemisphere middle latitude locations. The created database incorporated 65 strong-to-severe geomagnetic storms, which occurred within the period 1995–2007. This validation included data from some ionospheric stations (e.g., Pruhonice, El Arenosillo) that were not considered in the development or previous validations of the model. Hourly values of the F2 layer critical frequency, foF2, measured for 5–7 days during the main and recovery phases of each selected storm were compared with the predicted IRI 2007 foF2 with the STORM model option activated. To perform a detailed comparison between observed values, medians and predicted foF2 values the correlation coefficient, the root-mean-square error (RMSE), and the percentage improvement were calculated. Results of the comparative analysis show that the STORM model captures more effectively the negative phases of the summer ionospheric storms, while electron density enhancement during winter storms and the changeover of the different storm phases is reproduced with less accuracy. The STORM model corrections are less efficient for lower-middle latitudes and severe geomagnetic storms.     

Effect of vertical plasma transport on ionospheric F2-region parameters at equatorial latitude

The variability of the F2-layer even during magnetically quiet times are fairly complex owing to the effects of plasma transport. The vertical E × B drift velocities (estimated from simplified electron density continuity equation) were used to investigate the seasonal effects of the vertical ion drifts on the bottomside daytime ionospheric parameters over an equatorial latitude in West Africa, Ibadan, Nigeria (Geographic: 7.4°N, 3.9°E, dip angle: 6°S) using 1 year of ionsonde data during International Geophysical Year (IGY) of 1958, that correspond to a period of high solar activity for quiet conditions. The variation patterns between the changes of the vertical ion drifts and the ionospheric F2-layer parameters, especially; f_oF₂ and h_mF₂ are seen remarkable. On the other hand, we observed strong anti-correlation between vertical drift velocities and h′F in all the seasons. We found no clear trend between N_mF₂ and h_mF₂ variations. The yearly average value of upward daytime drift at 300 km altitude was a little less than the generally reported magnitude of 20 ms⁻¹ for equatorial F-region in published literature, and the largest upward velocity was roughly 32 ms⁻¹. Our results indicate that vertical plasma drifts; ionospheric F2-layer peak height, and the critical frequency of F2-layer appear to be somewhat interconnected.     

Development of a campaign to study equatorial ionospheric phenomena over Guam

The United States Air Force Academy (USAFA) is in the process of developing a series of ground-based and space-based experiments to investigate the equatorial ionosphere over Guam and the southern crest of the Equatorial Appleton Anomaly over New Guinea. On the ground the Digital Ionospheric Sounder (University of Massachusetts, Lowell DPS-4 unit) and a dual-frequency GPS TEC/scintillation monitor will be used to investigate ionospheric phenomena in both campaign and long-term survey modes. In campaign mode, we will combine these observations with those collected from space during USAFA’s FalconSAT-3 and FalconSAT-5 low Earth orbit satellite missions, which will be active over a period of several years beginning in the first quarter of the 2007 calendar year. Additionally, we will investigate the long-term morphology of key ionospheric characteristics useful for driving the International Reference Ionosphere, such as critical frequencies (f_oE, f_oF1, f_oF2, etc.), the M(3000) F2 parameter (the maximum useable frequency for a signal refracted within the F2 layer and received on the ground at a distance of 3000 km away), and a variety of other characteristics. Specific targets of investigation include: (a) a comparison of TEC observed by the GPS receiver with those calculated by IRI driven by DPS-4 observations, (b) a comparison of plasma turbulence observed on-orbit with ionospheric conditions as measured from the ground, and (c) a comparison between topside ionospheric satellite in situ measurements of plasma density during an overpass of a Digisonde versus the calculated value based on extrapolation of the electron density profiles using Digisonde data and a topside α-Chapman function. This last area of investigation is discussed in detail in this paper.     

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.     

Acceleration of energetic particles by whistler waves in active space experiment with charged particle beams injection

In this paper the investigation of wave-particle interaction during simultaneous injection of electron and xenon ion beams from the satellite Intercosmos-25 (IK-25) carried out using the data of the double satellite system with subsatellite Magion-3 (APEX). Results of active space experiment devoted to the beam-plasma instability are partially presented in the paper Baranets et al. (2007). A specific feature of the experiment carried out in orbits 201, 202 was that charged particle flows were injected in the same direction along the magnetic field lines B₀ so the oblique beam-into-beam injection have been produced. Results of the beam-plasma interaction for this configuration were registered by scientific instruments mounted on the station IK-25 and Magion-3 subsatellite. Main attention is paid to study the electromagnetic and longitudinal waves excitation in different frequency ranges and the energetic electron fluxes disturbed due to wave-particle interaction with whistler waves. The whistler wave excitation on the 1st electron cyclotron harmonic via normal Doppler effect during electron beam injection in ionospheric plasma are considered.