A real-time ionospheric model based on GNSS Precise Point Positioning

This paper proposes a method of real-time monitoring and modeling the ionospheric Total Electron Content (TEC) by Precise Point Positioning (PPP). Firstly, the ionospheric TEC and receiver’s Differential Code Biases (DCB) are estimated with the undifferenced raw observation in real-time, then the ionospheric TEC model is established based on the Single Layer Model (SLM) assumption and the recovered ionospheric TEC. In this study, phase observations with high precision are directly used instead of phase smoothed code observations. In addition, the DCB estimation is separated from the establishment of the ionospheric model which will limit the impacts of the SLM assumption impacts. The ionospheric model is established at every epoch for real time application. The method is validated with three different GNSS networks on a local, regional, and global basis. The results show that the method is feasible and effective, the real-time ionosphere and DCB results are very consistent with the IGS final products, with a bias of 1–2 TECU and 0.4 ns respectively.     

Analysis of ionospheric TEC response to solar and geomagnetic activities at different solar activity stages

Studying the relationship of total electron content (TEC) to solar or geomagnetic activities at different solar activity stages can provide a reference for ionospheric modeling and prediction. On the basis of solar activity indices, geomagnetic activity parameters, and ionospheric TEC data at different solar activity stages, this study analyzes the overall variation relationships of solar and geomagnetic activities with ionospheric TEC, the characteristics of the quasi-27-day periodic oscillations of the three variables at different stages, and the delayed TEC response of solar activity by conducting correlation analysis, Butterworth band-pass filtering, Fourier transform, and time lag analysis. The following results are obtained. (1) TEC exhibits a significant linear relationship with solar activity at different solar activity stages. The correlation coefficients |R| are arranged as follows: |R|_EUV > |R|_F10.7 > |R|sunspot number. No significant linear relationship exists between TEC and geomagnetic activity parameters (|R| < 0.35). (2) TEC, solar activity indices, and geomagnetic activity parameters have a period of 10.5 years. The maximum amplitudes of the Fourier spectrum for TEC and solar activity indices are nearly 27 days and those of geomagnetic activity parameters are nearly 27 and 13.5 days. (3) The deviations of the quasi-27-day significant periodic oscillation of TEC and solar activity indices are consistent. (4) No evident relationship exists between the quasi-27-day periodic oscillation of TEC and geomagnetic activity parameters. (5) The delay time of TEC for the 10.7 cm solar radio flux and extreme ultraviolet is always consistent, whereas that for sunspot number varies at each stage.     

SHAKING: Adjusted spherical harmonics adding KrigING method for near real-time ionospheric modeling with multi-GNSS observations

Precise positioning based on Global Navigation Satellite System (GNSS) technique requires high accuracy ionospheric total electron content (TEC) correction models to account for the ionospheric path delay errors. We present an adjusted Spherical Harmonics Adding KrigING method (SHAKING) approach for regional ionospheric vertical TEC (VTEC) modeling in real time. In the proposed SHAKING method, the VTEC information over the sparse observation data area is extrapolated by the Adjusted Spherical Harmonic (ASH) function, and the boundary distortion in regional VTEC modeling is corrected by the stochastic VTEC estimated using Kriging interpolation. Using real-time GPS, GLONASS and BDS-2/3 data streams of the Crust Movement Observation Network of China (CMONOC), the SHAKING-based regional ionospheric VTEC maps are re-constructed over China and its boundary regions. Compared to GNSS VTECs derived from the independent stations, the quality of SHAKING solution improves by 13–31% and 6–33% with respect to the ASH-only solution during high and low geomagnetic periods, respectively. Compared to the inverse distance weighting (IDW) generated result, significant quality improved of SHAKING-based VTEC maps is also observed, especially over the edge areas with an improvement of 60–80%. Overall, the proposed SHAKING method exhibits notable advantage over the existing regional VTEC modeling techniques, which can be used for regional TEC modeling and associated high-precision positioning applications.     

Aspects related to variability of radiative cooling by NO in lower thermosphere,TEC and O/N2 correlation,and diffusion of NO into mesosphere during the Halloween storms

Nitric Oxide is a very important trace species which plays a significant role acting as a natural thermostat in Earth’s thermosphere during strong geomagnetic activity. In this paper, we present various aspects related to the variation in the NO Infrared radiative flux (IRF) exiting the thermosphere by utilizing the TIMED/SABER (Thermosphere Ionosphere Mesosphere Energetics and Dynamics/ Sounding of the Atmosphere using Broadband Emission Radiometry) observational data during the Halloween storm which occurred in late October 2003. The Halloween storm comprised of three intense-geomagnetic storms. The variability of NO infrared flux during these storm events and its connection to the strength of the geomagnetic storms were found to be different in contrast to similar super storms. The connection between the quantum of energy outflux from the upper atmosphere into space in terms of NO IRF and the duration of storms is established. The NO radiative cooling, and the closely correlated depletion in O/N₂ ratio are controlled by the Joule heating intensity (proxied by AE-index). The collisional excitation rate of NO, calculated using the modelled datasets of WACCM-X (Whole Atmosphere Community Climate Model with thermosphere and ionosphere extension), correlates well with the observed pattern of radiative emission by NO. Observational datasets from TIMED/GUVI (Global Ultra-Violet Imager) and MIT Haystack observatory madrigal GNSS (Global navigation satellite system) total electron content (TEC) database shows that the TEC and O/N₂ enhancement in low-mid northern hemispheric latitudes are mainly controlled by the z-component of Interplanetary magnetic field (IMF-B_z). The penetration of eastward electric field during the storm events is found to be responsible for the overall enhancement of TEC. The contribution of enhanced day-side TEC in observed variation of O/N₂ ratio by GUVI is also reported. It is also seen that during substorms related events the night-time polar region experiences more cooling due to NO than the daytime polar region. The connections between the mid- and low-latitude enhancement in NO IRF with the propagation of LSTIDs (Large-scale traveling ionospheric disturbances) in combination with the O/N₂ variability, and the altitudinal variation in NO flux with the progression of the storm is also investigated. This study presents the evidence on the role of diffusion processes in the large scale enhancement of NO in the mesospheric altitudes.     

Ionospheric storm due to solar Coronal mass ejection in September 2017 over the Brazilian and African longitudes

Coronal mass ejection (CME) occurs when there is an abrupt release of a large amount of solar plasma, and this cloud of plasma released by the Sun has an intrinsic magnetic field. In addition, CMEs often follow solar flares (SF). The CME cloud travels outward from the Sun to the interplanetary medium and eventually hits the Earth’s system. One of the most significant aspects of space weather is the ionospheric response due to SF or CME. The direction of the interplanetary magnetic field, solar wind speed, and the number of particles are relevant parameters of the CME when it hits the Earth’s system. A geomagnetic storm is most geo-efficient when the plasma cloud has an interplanetary magnetic field southward and it is accompanied by an increase in the solar wind speed and particle number density. We investigated the ionospheric response (F-region) in the Brazilian and African sectors during a geomagnetic storm event on September 07–10, 2017, using magnetometer and GPS-TEC networks data. Positive ionospheric disturbances are observed in the VTEC during the disturbed period (September 07–08, 2017) over the Brazilian and African sectors. Also, two latitudinal chains of GPS-TEC stations from the equatorial region to low latitudes in the East and West Brazilian sectors and another chain in the East African sector are used to investigate the storm time behavior of the equatorial ionization anomaly (EIA). We noted that the EIA was disturbed in the American and African sectors during the main phase of the geomagnetic storm. Also, the Brazilian sector was more disturbed than the African sector.