Evaluation of the IRI-2016 and NeQuick electron content specification by COSMIC GPS radio occultation,ground-based GPS and Jason-2 joint altimeter/GPS observations

We examined performance of two empirical profile-based ionospheric models, namely IRI-2016 and NeQuick-2, in electron content (EC) and total electron content (TEC) representation for different seasons and levels of solar activity. We derived and analyzed EC estimates in several representative altitudinal intervals for the ionosphere and the plasmasphere from the COSMIC GPS radio occultation, ground-based GPS and Jason-2 joint altimeter/GPS observations. It allows us to estimate a quantitative impact of the ionospheric electron density profiles formulation in several altitudinal intervals and to examine the source of the model-data discrepancies of the EC specification from the bottom-side ionosphere towards the GPS orbit altitudes. The most pronounced model-data differences were found at the low latitude region as related to the equatorial ionization anomaly appearance. Both the IRI-2016 and NeQuick-2 models tend to overestimate the daytime ionospheric EC and TEC at low latitudes during all seasons of low solar activity. On the contrary, during high solar activity the model results underestimated the EC/TEC observations at low latitudes. We found that both models underestimated the EC for the topside ionosphere and plasmasphere regions for all levels of solar activity. For low solar activity, the underestimated EC from the topside ionosphere and plasmasphere can compensate the overestimation of the ionospheric EC and, consequently, can slightly decrease the resulted model overestimation of the ground-based TEC. For high solar activity, the underestimated EC from the topside ionosphere and plasmasphere leads to a strengthening of the model underestimation of the ground-based TEC values. We demonstrated that the major source of the model-data discrepancies in the EC/TEC domain comes from the topside ionosphere/plasmasphere system.     

Numerical study of traveling ionosphere disturbances with vertical incidence data

High frequency ionosphere vertical sounding, as an important and representative application for detecting the ionosphere and studying the characteristics of radio propagation, can be utilized to monitor the ionosphere continuously variation and to acquire the ionosphere asymmetrical features of diverse scale above the ionosphere vertical sounding stations. This is a first article on real time application of numerical methods to obtain the parameters of traveling ionosphere disturbances (TIDs) using vertical incident ionograms. In this paper, the distribution of ionosphere electron density with TIDs is constructed using a background ionosphere model superimposed a perturbation theory model. The background ionosphere electron density is modelled by the inversion of vertical incident ionograms which are observed before the appearance of the disturbance. Based on the fourth order Adams-Bashforth-Moulton (the so-called ABM) predictor corrector method, instead of Runge-Kutta method, the fast digital ray tracing method is established. According to process of the disturbed trace simulation and parameters inversion, the characteristic parameters of ionosphere disturbance at different detection time can be obtained in real time. The numerical analysis of TIDs is then captured completely.     

Case study of inclined sporadic E layers in the Earth’s ionosphere observed by CHAMP/GPS radio occultations: Coupling between the tilted plasma layers and internal waves

We have used the radio occultation (RO) satellite data CHAMP/GPS (Challenging Minisatellite Payload/Global Positioning System) for studying the ionosphere of the Earth. A method for deriving the parameters of ionospheric structures is based upon an analysis of the RO signal variations in the phase path and intensity. This method allows one to estimate the spatial displacement of a plasma layer with respect to the ray perigee, and to determine the layer inclination and height correction values. In this paper, we focus on the case study of inclined sporadic E (E_s) layers in the high-latitude ionosphere based on available CHAMP RO data. Assuming that the internal gravity waves (IGWs) with the phase-fronts parallel to the ionization layer surfaces are responsible for the tilt angles of sporadic plasma layers, we have developed a new technique for determining the parameters of IGWs linked with the inclined E_s structures. A small-scale internal wave may be modulating initially horizontal E_s layer in height and causing a direction of the plasma density gradient to be rotated and aligned with that of the wave propagation vector k. The results of determination of the intrinsic wave frequency and period, vertical and horizontal wavelengths, intrinsic vertical and horizontal phase speeds, and other characteristics of IGWs under study are presented and discussed.     

Analysis of the data processing strategies of spherical harmonic expansion model on global ionosphere mapping for moderate solar activity

Spherical harmonic (SH) expansion is widely used to model the global ionosphere map (GIM) of vertical total electron content (VTEC). According to the impact of different data processing methods of the SH expansion model on the VTEC maps, we specifically performed comprehensive analysis in terms of the data sampling rate, the time resolution, the spherical harmonic degree, and the relative constraint. One month of GPS data (January in 2016) from the International GNSS (Global Navigation Satellite System) Service (IGS) network in a moderate ionospheric activity period at the descending phase of Solar Cycle 24 was processed. To improve the computational efficiency of the daily GIM generation, the data sampling rate of 5?min was recommended allowing the GIM precision loss within 0.10 TECU (total electron content unit). The global VTEC map could be better represented in temporal and spatial domains with higher time resolution and higher spherical harmonic degree, especially at low latitude bands and in the southern hemisphere. The GIM precision improvement was about 10.91% for 1-h and about 15.15% for 0.5-h compared with the commonly used 2-h time resolution. The use of spherical harmonic degree 17 or 20 instead of 15 could improve the precision by 3.19% or 6.06%. We also found that an optimal relative constraint had to be found experimentally considering both the GIM precision and the GIM root mean square (RMS) map.     

Impact of geomagnetic variation over sub-auroral ionospheric region during high solar activity year 2014

The present work is an attempt to evaluate the impact of changing space weather condition over sub-auroral ionosphere during high solar activity year 2014. In view of this, the GPS based TEC along with Ionosonde data over Indian permanent scientific base “Maitri”, Antarctica (70°46′00″S, 11°43′56″E) has been utilized. The results suggested that the nature of ionospheric responses to the geomagnetic disturbances not only depended upon the status of high latitudinal electro-dynamic processes but also influenced by the seasonal variations. The results revel both negative and positive type of ionospheric response in a single year but during different seasons. The study suggested that the combination of equator-ward plasma transportation along with ionospheric compositional changes causes a negative ionospheric impact especially during summer and equinox seasons. However, the combination of pole-ward contraction of the oval region along with particle precipitation may lead to exhibit positive ionospheric response during the winter season. The plasma transportation direction has been validated with the help of convection boundary (HM boundary) deduced with the help of SuperDARN observations. The ground based ionosonde observations clearly provided the evidence of deep penetration of high energetic particles up to the E-layer heights which results a sudden and strong appearance of E-layer. The strengthening of E-layer is responsible for modification of auroral electrojet and field-aligned current system. Also, the sudden appearance of E-layer along with a decrease in F-layer electron density suggested the dominance of NO⁺ over O⁺ in a considered region under geomagnetic disturbed condition.     

Satellite launch vehicle effect on the Earth’s lower ionosphere: A case study

The effect of the rocket exhaust products on the D-region of the ionosphere is investigated with the help of Very low frequency (VLF) electromagnetic wave propagation characteristics within the Earth-ionosphere waveguide. The changes in the electron density profile are computed from the observed VLF signal amplitude perturbations about 3 dB during the rocket launch. We find a localized electron depletion in the lower ionosphere at an altitude of around 58 km, that is thought to be originated by the attachment of ionospheric ion and molecular hydrogen along with water molecule in the exhaust product of first stage burn of Geosynchronous Launch Vehicle (GSLV) rocket at the time of GSLV launched from Sriharikota, India, on 27 August 2015 at 11:22 UT (16:52 IST). The ionospheric depletion perturbed the navigational VLF signal (VTX = 17 kHz) 134 s after the launch of the GSLV rocket.     

Characteristics of field line resonance type geomagnetic pulsations and variations of plasmaspheric plasma density distribution

The period of field line resonance (FLR) type geomagnetic pulsations depends on the length of the field line and on the plasma density in the inner magnetosphere (plasmasphere), where field lines are closed. Here as FLR period, the period belonging to the maximum occurrence frequency of the occurrence frequency spectrum (equivalent resonance curve) of pulsations has been considered. The resonance system may be replaced by an equivalent resonant circuit. The plasma density would correspond to the ohmic load. The plasma in the plasmasphere originates from the ionosphere, thus FLR period, occurrence frequency are also affected by the maximum electron density in the ionosphere. The FLR period has shown an enhancement with increasing F region electron density, while the occurrence frequency indicated diminishing trend (possible damping effect). Thus, the increased plasma density may be the cause of the decreased occurrence of FLR type pulsations in the winter months of solar activity maximum years (winter anomaly).     

Application of small-size antennas for estimation of angles of arrival of HF signals scattered by ionospheric irregularities

A modification of the Doppler Interferometry Technique is suggested to enable estimating angles of arrival of comparatively broadband HF signals scattered by random irregularities of the ionospheric plasma with the use of small-size weakly directional antennas. The technique is based on the measurements of cross-spectra phases of the probe radiation recorded at least in three spatially separated points. The developed algorithm has been used to investigate the angular and frequency-time characteristics of HF signals propagating at frequencies above the maximum usable one (MUF) for the direct radio path Moscow-Kharkiv. The received signal spectra show presence of three families of spatial components attributed, respectively, to scattering by plasma irregularities near the middle point of the radio path, ground backscatter signals and scattering of the sounding signals by the intense plasma turbulence associated with auroral activations. It has been shown that the regions responsible for the formation of the third family components are located well inside the auroral oval. The drift velocity and direction of the auroral ionosphere plasma have been determined. The obtained estimates are consistent with the classical conception of the ionospheric plasma convection at high latitudes and do not contradict the results of investigations of the auroral ionosphere dynamics using the SuperDARN network.     

TEC differences for the mid-latitude ionosphere in both sides of the longitudes with zero declination

Based on measurements of ground-based GPS station network, differences of the mid-latitude ionospheric TEC in the east and west sides of North America, South America and Oceania have been analyzed in this paper. Results show that for nearly all seasons from 2001 to 2010 and in both sides of the longitudes with zero declination, there exist systematic differences for the mid-latitude ionospheric TEC in the regions mentioned above and the features of these differences markedly depend upon the local time but less depend upon seasons and the level of solar activity. Theory analysis shows that the longitude variations of both declination and zonal thermospheric winds are one of important factors to cause differences of the mid-latitude ionospheric TEC in both sides of the longitudes with zero declination.     

Eliminating negative VTEC in global ionosphere maps using inequality-constrained least squares

Currently, ground-based Global Navigation Satellite System (GNSS) stations of the International GNSS Service (IGS) are distributed unevenly around the world. Most of them are located on the mainland, while only a small part of them are scattered on some islands in the oceans. As a consequence, many unreasonable zero values (in fact negative values) appear in Vertical Total Electron Content (VTEC) of European Space Agency (ESA) and Center for Orbit Determination in Europe (CODE) IONEX products, especially in 2008 and 2009 when the solar activities were rather quiet. To improve this situation, we directly implement non-negative physical constraints of ionosphere for global ionosphere maps (GIM) with spherical harmonic functions. Mathematically, we propose an inequality-constrained least squares method by imposing non-negative inequality constraints in the areas where negative VTEC values may occur to reconstruct GIM models. We then apply the new method to process the IGS data in 2008. The results have shown that the new algorithm efficiently eliminates the unwanted behavior of negative VTEC values, which could otherwise often be seen in the current CODE and ESA GIM products in both middle and high latitude areas of the Southern Hemisphere (45°S∼90°S) and the Northern Hemisphere (50°N∼90°N). About 64% of GPS receivers’ DCBs have been significantly improved. Finally, we compare the GIM results between with and without the inequality constraints, which has clearly shown that the GIM result with inequality constraints is significantly better than that without the inequality constraints. The inequality-constrained GIM result is also highly consistent with the final IGS products in terms of root mean squared (RMS) and mean VTEC.