Assessment of the NeQuick model at mid-latitudes using GNSS TEC and ionosonde data

The modelling of the total electron content (TEC) plays an important role in global navigation satellite systems (GNSS) accuracy, especially for single-frequency receivers, the most common ones constituting the mass market. For the latter and in the framework of Galileo, the NeQuick model has been chosen for correcting the ionospheric error contribution and will be integrated into a global algorithm providing the users with daily updated information.     

Electric currents from thunderstorms to the ionosphere during a solar cycle: Quasi-static modeling of the coupling mechanism

Conduction and displacement currents, and their sum the Maxwell current, generated over a thunderstorm (TS) with recurrent lightning discharges are investigated theoretically. The aim is to study better the influence of different factors on these currents, which form the link between thunderstorms and the ionosphere in the global atmospheric electrical circuit. The factors studied concern the thunderstorm characteristics (the charge separation current, and the lightning discharge parameters), as well as the atmospheric and cloud conductivity. Some of these factors may show long-term changes with the 11-year solar cycle, possibly realized through an inverse dependence of the cosmic ray flux on solar activity. Earlier investigations have suggested that the lightning-related charge redistribution and subsequent relaxation, rather than the high intensity current, is mainly the source of the energy coupled to the ionosphere. With respect to this, a quasi-electrostatic analytical model is proposed, based on Maxwell’s equations. The currents are generated by a TS modeled as a positive vertical dipole with charges which are first accumulated and then destroyed by lightning. Our computations show that the mean and peak values of the conduction and total Maxwell currents to the ionosphere depend significantly on the charge moment change. The mean currents are also sensitive to the reduction of the conductivity in thunderclouds. Small variations of the stratospheric conductivity (20% at geomagnetic latitude 40° and 40–50% at 55°) with the solar activity do not influence the currents to the ionosphere very much.     

Deformation of the ionosphere structure during the space weather events on October–November 2003

We present the spatial maps of the ionosphere–plasmasphere slab thickness τ (ratio of the vertical total electron content, TEC, to the F-region peak electron density, NmF2) during the intense ionospheric storms of October–November 2003. The model-assisted technology for estimate of the upper boundary of the ionosphere, hup, from the slab thickness components in the bottomside and topside ionosphere – eliminating the plasmasphere contribution of τ – is applied at latitudes 35° to 70°N and longitudes −10° to 40°E, from the data of 20 observatories of GPS-TEC and ionosonde networks, for selected days and hours of October and November 2003. The daily–hourly values of NmF2, hmF2 and TECgps are used as the constrained parameters for the International Reference Ionosphere extended to the plasmasphere, IRI-Plas, during the ionospheric quiet days, positive and negative storm phases for estimate of τ and hup. Good correlation has been found between the slab thickness and the upper boundary of the ionosphere for the intense ionospheric storms at October–November 2003. During the negative phase of the ionospheric storm, when the ionospheric plasma density is exhausted, the nighttime upper boundary of the ionosphere is greatly uplifted towards the magnetosphere tail, while the daytime upper boundary of the ionosphere is reduced below 500 km over the Earth.     

An improved approach to model regional ionosphere and accelerate convergence for precise point positioning

Given the severe effects of the ionosphere on global navigation satellite system (GNSS) signals, single-frequency (SF) precise point positioning (PPP) users can only achieve decimeter-level positioning results. Ionosphere-free combinations can eliminate the majority of ionospheric delay, but increase observation noise and slow down dual-frequency (DF) PPP convergence. In this paper, we develop a regional ionosphere modeling and rapid convergence approach to improve SF PPP (SFPPP) accuracy and accelerate DF PPP (DFPPP) convergence speed. Instead of area model, ionospheric delay is modeled for each satellite to be used as a priori correction. With the ionospheric, wide-lane uncalibrated phase delay (UPD) and residuals satellite DCBs product, the wide-lane observations for DF users change to be high-precision pseudorange observations. The validation of a continuously operating reference station (CORS) network was analyzed. The experimental results confirm that the approach considerably improves the accuracy of SFPPP. For DF users, convergence time is substantially reduced.     

国际参考电离层IRI－90与中国实测电离层月平均电子密度剖面形状的比较

将太阳活动峰年期间中国4个电离层站垂直探测得到的月平均电子浓度剖面与国际参考电离层IKI-90进行了系统的比较。结果表明国际参考电离层所计算的峰下电子含量(或峰下半厚)总体来说偏大。一天中白天符合较好,晚上较差,对中纬台站较好,对低纬台站较差。      

On the detection of anomalous seismo-ionospheric behavior in the presence of space weather stimuli for large earthquakes

Anomalous behavior of ionospheric total electron content (TEC) prior to earthquake has been observed in many studies. Evidence of such seismo-ionospheric coupling effects suggests that it is plausible to rely on TEC signatures for early earthquake warning. However, the detection of pre-earthquake TEC anomalies (PETA) has not been adopted in practice due to two pertinent issues. Firstly, the effects of space weather activity can affect TEC levels and cause anomalous behavior in the TEC. Usually arbitrary thresholds are set for space weather indices to eliminate TEC anomaly due to space weather effects. Secondly, the choice regarding moving time-window length used to characterise background variation of TEC within the statistical envelope approach has an effect on detection of PETA. While the rule-of-thumb in selecting the moving window length is to have a time window capable of capturing background variability and short-term fluctuations, the length of the time window used in the literature varies with little justification. In this study, a critical examination is conducted on the statistical envelope approach and in particular, to eliminate the effect of space weather activity without the use of arbitrary space indices to detect PETA. A two-part PETA identification procedure is proposed, consisting of wavelet analyses isolating non-earthquake TEC contributions, followed by the statistical envelope method using a moving window length standardized based on observed periodicities and statistical implications is suggested. The approach is tested on a database of 30 large earthquakes (M?≥?7.0). The proposed procedure shows that PETA can be detected prior to earthquakes at higher confidence levels without the need to separately check for space weather activity. More importantly, the procedure was able to detect PETA for studies where it was previously reported that PETA could not be detected or detected convincingly.     

Kalman filter based estimation of differential hardware biases with triangular interpolation technique for IRNSS

Ionosphere delay correction is the main error correction to the computation of single frequency user position using satellite navigation. However ionosphere delay consists of not only delay but also frequency dependent differential hardware biases from satellite and receiver ends. For ionosphere point of view, Indian Regional Navigation Satellite System (IRNSS) service area comes in equatorial anomaly region. It is a unique satellite navigation system which operates at L5 and S frequencies and consists of Geostationary Earth Orbit (GEO) and Geo Synchronous Orbit (GSO) satellite constellation. With IRNSS measurements availability, there is a good opportunity to estimate and analyse differential hardware biases with GEO/GSO combination and with equatorial ionosphere variation. In this paper, Kalman filter based estimation with triangular interpolation technique is used to estimate differential hardware biases for all IRNSS satellites and reference receivers at L5 frequency. The standard deviation of the 15?days of daily estimation of satellite differential hardware biases is in the range of 0.32 to 1.17 TECU for all IRNSS satellites. Similarly, the standard deviation of the 15?days of daily estimation varies up to 2.85 and 6.0 TECU for receiver differential hardware biases during calm and stormy period respectively. The ionosphere delay computed using estimated differential hardware biases is compared with Global Ionosphere Map (GIM) data. A rigorous analysis is carried out to study the error in the estimation in terms of input data noise level, satellite constellation and effect of latitude. Our result reveals that over IRNSS service area, there is an exponential increase in the error in the estimation of receiver differential hardware biases with respect to latitude.     

A realistic simulation framework to evaluate ionospheric tomography

Observation of the 3-dimensional (3-D) electron density of the ionosphere is useful to study large-scale physical processes in space weather events. Ionospheric data assimilation and ionospheric tomography are methods that can create an image of the 3-D electron density distribution. While multiple techniques have been developed over the past 30 years, there are relatively few studies that show the accuracy of the algorithms. This paper outlines a novel simulation approach to test the quality of an ionospheric tomographic inversion. The approach uses observations from incoherent scatter radar (ISR) scans and extrapolates them spatially to create a realistic ionospheric representation. A set of total electron content (TEC) measurements can then be simulated using real geometries from satellites and ground receivers. This data set, for which the ‘truth’ ionosphere is known, is used as input for a tomographic inversion algorithm to estimate the spatial distribution of electron density. The reconstructed ionospheric maps are compared with the truth ionosphere to calculate the difference between the images and the truth.To demonstrate the effectiveness of this simulation framework, an inversion algorithm called MIDAS (Multi-Instrument Data Analysis Software) is evaluated for three geographic regions with differing receiver networks. The results show the importance of the distribution and density of GPS receivers and the use of a realistic prior conditioning of the vertical electron density profile. This paper demonstrates that when these requirements are met, MIDAS can reliably estimate the ionospheric electron density. When the region under study is well covered by GPS receivers, as in mainland Europe or North America, the errors in vertical total electron content (vTEC) are smaller than 1 TECu (2–4%) . In regions with fewer and more sparsely distributed receivers, the errors can be as high as 20–40%. This is caused by poor data coverage and poor spatial resolution of the reconstruction, which has an important effect on the calibration process of the algorithm.     

Modeling M(3000)F2 based on extreme learning machine

This paper presents an novel extreme learning machine (ELM)-based prediction model for the ionospheric propagation factor M(3000)F2 at Darwin station (12.4°S, 131.5°E; −44.5°dip) in Australia. The proposed ELM model is trained with hourly daily values of M(3000)F2 from the period 1998–2014 except 2001 and 2009. The hourly daily values of 2001 (high solar activity) and 2009 (low solar activity) are used for validating the prediction accuracy. The proposed ELM for modeling M(3000)F2 can achieve faster training process and similar testing accuracy compared with backward propagation neural network (BPNN). In addition, the performance of the ELM is verified by comparing the predicted values of M(3000)F2 with observed values and the international reference ionosphere (IRI −2016) model predicted values. Based on the error differences (the root mean square error (RMSE) and the M(3000)F2 percentage improvement values M(3000)F2_IMP(%)), the result demonstrates the effectiveness of the ELM model compared with the IRI-2016 model at hourly, daily, monthly, and yearly in high (2001) and low (2009) solar activity years. The ELM also shows good agreement with observations compared with the IRI during disturbed magnetic activity.     

Comparative analysis of four different single-frequency PPP models on positioning performance and atmosphere delay retrieval

Single-frequency precise point positioning (SF-PPP) has attracted increasing attention due to its high precision and cost effectiveness. With various strategies to handle the dominant error, i.e., ionosphere delay, the ionosphere-float (IF), ionosphere-free-half (IFH), ionosphere-corrected (IC), and ionosphere-weighted (IW) SF-PPP models are certain to possess different characteristics and performance levels. This study is dedicated to assessing and comparing the four models from model characteristics, positioning performance, and atmosphere delay retrieval. The model comparison shows that IC and IW models are full-rank while IF and IFH models have a rank deficiency of size one that will result in biased estimations, which means the better solvability of IC and IW models. The experiments are carried out based on the 7-day Global Positioning System (GPS) observations collected at 57 global Multi-GNSS Experiment (MGEX) stations and Global Ionosphere Map (GIM) products. The results indicate that the IW model can accelerate SF-PPP convergence and achieve higher positioning accuracy compared to the other three SF-PPP models, especially in kinematic mode. With convergence criteria of 0.25 m in horizontal and 0.5 m in vertical, the east/north/up convergence times of IW model are 0.5/15.0/25.0 min and 0.5/16.0/36.5 min for static and kinematic modes, respectively. The IW model is able to achieve an instantaneous positioning accuracy of 0.28/0.35/0.75 m. In addition, a real kinematic test also demonstrates the best positioning solutions of IW model. Regarding troposphere delay retrieval, the IF, IFH, and IW models obtain a comparable daily accuracy of 3.0 cm on average, while the IC model achieves the worst accuracy of 8.0 cm. For precise ionosphere delay estimation, IW model only needs an average initialization time of 34.3 min, but a longer initialization time of 51.6 min is required for IF model. The daily precision of ionosphere delay estimation for IW model can reach up to 10.8 cm. At the present accuracy of GIM products, it is suggested that the IW model should be adopted for SF-PPP first due to its superior performance in positioning and atmosphere delay retrieval.