Ionospheric differential error determination using ray tracing for a short baseline

Since the United States government discontinued Selective Availability (SA) on 1 May 2000, ionospheric effects have been responsible for the largest errors in GPS systems. The standard Differential GPS (DGPS) method is incapable of completely eliminating the ionospheric error. This paper describes a new approach to determine the differential ionospheric error between geographically distributed receiver stations. The ray paths of GPS signals were simulated using a modified Jones 3D ray tracing programme that includes the effect of the geomagnetic field. A Nelder–Mead optimisation algorithm was embedded in the program to precisely determine the satellite-to-station path. A realistic ionospheric model is essential for accurate ray tracing results and for estimates of differential error that are accurate on sub-centimetre scales. Here, the ionospheric model used in the ray tracing programme was developed by fitting realistic ionosphere profiles with a number of exponential functions. Results were compared to the theoretical approach. Results show that the differential delay is about 1–5 cm at low elevation angles for a short baseline of 10 km, as reported in other literature. This delay is often neglected in DGPS application. The differential delay also shows a pattern similar to that predicted by the Klobuchar model. The method proposed here can be used to improve future GPS applications.     

Spatial distribution of TEC across India in 2005: Seasonal asymmetries and IRI prediction

Total electron content measured simultaneously at 10 locations over India during the low solar activity year 2005 is used to examine the temporal and spatial asymmetries and also to assess the predictability of the International Reference Ionosphere in respect of the observed asymmetrical distribution. The stations are distributed in latitude along 77°E and in longitude along 23°N forming a meridional and a zonal chain respectively. A longitudinal gradient positive towards east was observed in the daytime hours of equinox and summer. Equinoctial asymmetry was prevalent across India during this year. Within the crest and equator, winter anomaly has been observed. It is found that IRI 2012 (with Ne Quick option, URSI coefficients) is unable to fully capture the temporal variation and spatial gradients of the ionization density in the Indian sector during 2005. The amount of offset between the model and measurement varies with local time and location.     

Prediction of ionospheric total electron content using adaptive neural network with in-situ learning algorithm

The Ionospheric Total Electron Content is responsible for the group delay of the signals from the Navigation satellites. This delay causes ranging error, which in turn degrades the accuracy of position estimated by the receivers. For critical applications, single frequency receivers resort to Satellite Based Augmentation Systems in order to have improved accuracy and integrity. The performance of these systems in terms of accuracy can be improved if predictions of the delays are available simultaneously with real measurements. This paper attempts to predict the Total Electron Content using adaptive recurrent Neural Network at three different locations of India. These locations are selected at the magnetic equator, at the equatorial anomaly crest and outside the anomaly range, respectively. In-situ Learning Algorithm has been used for tracking the non-stationary nature of the variation. Prediction is done for different prediction intervals. It is observed that, for each case, the mean and root mean square values of prediction errors remain small enough for all practical applications. Analysis of Variance is also done on the results.     

Change in refractivity of the atmosphere and large variation in TEC associated with some earthquakes,observed from GPS receiver

The present study reports the analysis of GPS based TEC for our station Surat (21.16°N, 72.78°E) located at the northern crest of equatorial anomaly region in India at times close to some earthquake events (M ? 5) during the year 2009 in India and its neighbouring regions. The TEC data used in the study are obtained from GPS Ionospheric Scintillation and TEC Monitoring (GISTM) system. The TEC data has been analysed corresponding to 11 earthquakes in low solar activity period and quiet geomagnetic condition. We found that, out of 11 cases of earthquakes (M > 5) there were seven cases in which enhancement in TEC occurred on earthquake day and in other four cases there was depletion in TEC on earthquake day. The variation in refractivity prior to earthquake was significant for the cases in which the epicentre lied within a distance of 600 km from the receiving station. By looking into the features on temporal enhancement and depletion of TEC a prediction was made 3–2 days prior to an earthquake (on 28 October 2009 in Bhuj – India). The paper includes a brief discussion on the method of potentially identifying an impending earthquake from ionospheric data.     

Estimation of equatorial electrojet from total electron content at geomagnetic equator using Kalman filter

The temporal variation of the equatorial electrojet is estimated utilising a suitably designed Kalman filter and using the established empirical relations between the anomaly component of equatorial TEC and the modified electrojet. TEC data obtained from dual frequency GPS receivers are used for the purpose. Estimation requires the a-priori knowledge of the peak electrojet value of the day and hence can be made in post temporal scenario only. Initial results obtained during a low solar activity time in an equinoctial month shows acceptable accuracy of the proposed algorithm. Limited analysis is done by segregating the results into temporal sessions of pre-attainment and post-attainment of the electrojet peak.     

A simulation analysis of systematic errors in ion arrival angle and drift velocity from the satellite borne ion drift meter

As one payload of a Chinese seismic satellite program, an ion drift meter (IDM) will measure drift velocity of thermal ions at an altitude of 500 km. Previous works have shown that such instruments use biased grids to create nonuniform potential in the grid planes, which brings systematic errors to the inferred parameters. A commercial finite element analysis software is used to simulate this instrument in the exact size. The error sources from thermal velocity, nonuniform transparency of real grids and potential depression in the grid planes are explained. The simulation results show that the arrival angle and drift velocity will be underestimated in all the conditions and the maximal error will be about −0.87° and −121 m/s, respectively. Furthermore, the relative error of the inferred arrival angle and the drift velocity will be inversely correlated with the arrival angle because of the lensing effect of the potential depression. This simulation provides a quantificational method of evaluating and correcting the data during in situ operation.     

On the application of Exploratory Data Analysis for characterization of space weather data sets

This paper presents an overview of the mathematical foundations for techniques in Exploratory Data Analysis (EDA) for the purpose of investigating the relationships among the numerous variables in large sets of multivariate space weather data. Specifically, we cover techniques in Principal Components Analysis (PCA) and Common Factor Analysis (CFA). These techniques are illustrated using space weather activity indices collected during the year 2002 and the corresponding noon-time hmF2 data from the International Reference Ionosphere (IRI). A CFA is used to categorize the activity indices, and a PCA is used to derive two macro-indices of activity to ascertain the strength of solar and geomagnetic activity. These macro-indices are then used to compare and contrast IRI’s noon-time hmF2 values at six different geographic stations. It was found that the correlation between hmF2 and the macro-indices more accurately represented the variation of this correlation with latitude found in previous studies than if we used an isolated conventional index, such as SSN and AE. We also found that the daily maximum value of the Polar Cap Index was dependent on both solar and geomagnetic activity, but the closely-related cross-Polar Cap Potential was solely associated with elevated levels of geomagnetic activity, which is a unique result compared to previous studies. We argue that the discrepancy can be explained by the difference in experiment designs between the two studies. This paper demonstrates the usefulness of EDA in space weather studies of large multivariate data sets.     

The lower ionosphere response to its disturbances by powerful radio waves

The paper presents data from some campaigns at Sura heating facility in 2011–1016. The experiments on probing of the artificial disturbed region of the lower ionosphere were carried out at two observation sites. One of them was located near Vasil’sursk 1 km from Sura facility (56.1°N; 46.1°E) and the other site was located at the Observatory (55.85°N; 48.8°E) of Kazan State University, 170 km to the East. Investigation of the features of the disturbed region of the lower ionosphere based on its diagnostics by the methods of the vertical sounding and oblique backscattering is the main goal of this paper. Ionosphere disturbance was fulfilled by the effect of the powerful radio wave of the ordinary or extraordinary polarization emitted by transmitters of the Sura facility with effective radiated power ERP = 50–120 MW at the frequency of 4.3, 4.7 and 5.6 MHz. Pumping waves were emitted with period from 30 s to 15 min. The disturbed region of the ionosphere in Vasil’sursk was probed by the vertical sounding technique using the partial reflexion radar at the frequency of 2.95 and 4.7 MHz. For the oblique sounding of the disturbed region the modified ionosonde Cyclon-M, operating at ten frequencies from 2.01 to 6.51 MHz was used at the Observatory site. On many heating sessions simultaneous variations of the probing partial reflection signals in Vasil’sursk and backscattered signals in Observatory were observed at the height at 40–100 km below the reflection height of the pumping wave. These observations were correlated with the pumping periods of the Sura facility. Possible mechanisms of the appearance of the disturbance in the lower ionosphere and its effect on the probing radio waves are discussed.