Stochastic modelling considering ionospheric scintillation effects on GNSS relative and point positioning

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

RINEX_HO: second- and third-order ionospheric corrections for RINEX observation files

When GNSS receivers capable of collecting dual-frequency data are available, it is possible to eliminate the first-order ionospheric effect in the data processing through the ionosphere-free linear combination. However, the second- and third-order ionospheric effects still remain. The first-, second- and third-order ionospheric effects are directly proportional to the total electron content (TEC), although the second- and third-order effects are influenced, respectively, by the geomagnetic field and the maximum electron density. In recent years, the international scientific community has given more attention to these kinds of effects and some works have shown that for high precision GNSS positioning these effects have to be taken into consideration. We present a software tool called RINEX_HO that was developed to correct GPS observables for second- and third-order ionosphere effects. RINEX_HO requires as input a RINEX observation file, then computes the second- and third-order ionospheric effects, and applies the corrections to the original GPS observables, creating a corrected RINEX file. The mathematical models implemented to compute these effects are presented, as well as the transformations involving the earth's magnetic field. The use of TEC from global ionospheric maps and TEC calculated from raw pseudorange measurements or pseudoranges smoothed by phase is also investigated.

Improving the GNSS positioning stochastic model in the presence of ionospheric scintillation

Ionospheric scintillations are caused by time-varying electron density irregularities in the ionosphere, occurring more often at equatorial and high latitudes. This paper focuses exclusively on experiments undertaken in Europe, at geographic latitudes between similar to 50 degrees N and similar to 80 degrees N, where a network of GPS receivers capable of monitoring Total Electron Content and ionospheric scintillation parameters was deployed. The widely used ionospheric scintillation indices S4 and sigma(phi) represent a practical measure of the intensity of amplitude and phase scintillation affecting GNSS receivers. However, they do not provide sufficient information regarding the actual tracking errors that degrade GNSS receiver performance. Suitable receiver tracking models, sensitive to ionospheric scintillation, allow the computation of the variance of the output error of the receiver PLL (Phase Locked Loop) and DLL (Delay Locked Loop), which expresses the quality of the range measurements used by the receiver to calculate user position. The ability of such models of incorporating phase and amplitude scintillation effects into the variance of these tracking errors underpins our proposed method of applying relative weights to measurements from different satellites. That gives the least squares stochastic model used for position computation a more realistic representation, vis-a-vis the otherwise 'equal weights' model. For pseudorange processing, relative weights were computed, so that a 'scintillation-mitigated' solution could be performed and compared to the (non-mitigated) 'equal weights' solution. An improvement between 17 and 38% in height accuracy was achieved when an epoch by epoch differential solution was computed over baselines ranging from 1 to 750 km. The method was then compared with alternative approaches that can be used to improve the least squares stochastic model such as weighting according to satellite elevation angle and by the inverse of the square of the standard deviation of the code/carrier divergence (sigma CCDiv). The influence of multipath effects on the proposed mitigation approach is also discussed. With the use of high rate scintillation data in addition to the scintillation indices a carrier phase based mitigated solution was also implemented and compared with the conventional solution. During a period of occurrence of high phase scintillation it was observed that problems related to ambiguity resolution can be reduced by the use of the proposed mitigated solution.

Second and Third Order Ionospheric Effects on GNSS Positioning: A Case Study in Brazil

The Global Positioning System (GPS) transmits signals in two frequencies. It allows the correction of the first order ionospheric effect by using the ionosphere free combination. However, the second and third order ionospheric effects, which combined may cause errors of the order of centimeters in the GPS measurements, still remain. In this paper the second and third order ionospheric effects, which were taken into account in the GPS data processing in the Brazilian region, were investigated. The corrected and not corrected GPS data from these effects were processed in the relative and precise point positioning (PPP) approaches, respectively, using Bernese V5.0 software and the PPP software (GPSPPP) from NRCAN (Natural Resources Canada). The second and third order corrections were applied in the GPS data using an in-house software that is capable of reading a RINEX file and applying the corrections to the GPS observables, creating a corrected RINEX file. For the relative processing case, a Brazilian network with long baselines was processed in a daily solution considering a period of approximately one year. For the PPP case, the processing was accomplished using data collected by the IGS FORT station considering the period from 2001 to 2006 and a seasonal analysis was carried out, showing a semi-annual and an annual variation in the vertical component. In addition, a geographical variation analysis in the PPP for the Brazilian region has confirmed that the equatorial regions are more affected by the second and third order ionospheric effects than other regions.

Quality of TEC Estimated with Mod Ion Using GPS and GLONASS Data

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Higher order ionospheric effects in GNSS positioning in the European region

After removal of the Selective Availability in 2000, the ionosphere became the dominant error source for Global Navigation Satellite Systems (GNSS), especially for the high-accuracy (cm-mm) demanding applications like the Precise Point Positioning (PPP) and Real Time Kinematic (RTK) positioning.The common practice of eliminating the ionospheric error, e. g. by the ionosphere free (IF) observable, which is a linear combination of observables on two frequencies such as GPS L1 and L2, accounts for about 99% of the total ionospheric effect, known as the first order ionospheric effect (Ion1). The remaining 1% residual range errors (RREs) in the IF observable are due to the higher - second and third, order ionospheric effects, Ion2 and Ion3, respectively. Both terms are related with the electron content along the signal path; moreover Ion2 term is associated with the influence of the geomagnetic field on the ionospheric refractive index and Ion3 with the ray bending effect of the ionosphere, which can cause significant deviation in the ray trajectory (due to strong electron density gradients in the ionosphere) such that the error contribution of Ion3 can exceed that of Ion2 (Kim and Tinin, 2007).The higher order error terms do not cancel out in the (first order) ionospherically corrected observable and as such, when not accounted for, they can degrade the accuracy of GNSS positioning, depending on the level of the solar activity and geomagnetic and ionospheric conditions (Hoque and Jakowski, 2007). Simulation results from early 1990s show that Ion2 and Ion3 would contribute to the ionospheric error budget by less than 1% of the Ion1 term at GPS frequencies (Datta-Barua et al., 2008). Although the IF observable may provide sufficient accuracy for most GNSS applications, Ion2 and Ion3 need to be considered for higher accuracy demanding applications especially at times of higher solar activity.This paper investigates the higher order ionospheric effects (Ion2 and Ion3, however excluding the ray bending effects associated with Ion3) in the European region in the GNSS positioning considering the precise point positioning (PPP) method. For this purpose observations from four European stations were considered. These observations were taken in four time intervals corresponding to various geophysical conditions: the active and quiet periods of the solar cycle, 2001 and 2006, respectively, excluding the effects of disturbances in the geomagnetic field (i.e. geomagnetic storms), as well as the years of 2001 and 2003, this time including the impact of geomagnetic disturbances. The program RINEX_HO (Marques et al., 2011) was used to calculate the magnitudes of Ion2 and Ion3 on the range measurements as well as the total electron content (TEC) observed on each receiver-satellite link. The program also corrects the GPS observation files for Ion2 and Ion3; thereafter it is possible to perform PPP with both the original and corrected GPS observation files to analyze the impact of the higher order ionospheric error terms excluding the ray bending effect which may become significant especially at low elevation angles (Ioannides and Strangeways, 2002) on the estimated station coordinates.

Resolução espacial da Grade Ionosférica e do Give

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Efeitos de segunda e terceira ordem da ionosfera no posicionamento GNSS no Brasil

O Sistema de Posicionamento Global (GPS) transmite seus sinais em duas freqüências, o que permite eliminar matematicamente os efeitos de primeira ordem da ionosfera através da combinação linear ionosphere free. Porém, restam os efeitos de segunda e terceira ordem, os quais podem provocar erros da ordem de centímetros nas medidas GPS. Esses efeitos, geralmente, são negligenciados no processamento dos dados GPS. Os efeitos ionosféricos de primeira, segunda e terceira ordem são diretamente proporcionais ao TEC presente na ionosfera, porém, no caso dos efeitos de segunda e terceira ordem, comparecem também o campo magnético da Terra e a máxima densidade de elétrons, respectivamente. Nesse artigo, os efeitos de segunda e terceira ordem da ionosfera são investigados, sendo que foram levados em consideração no processamento de dados GPS na região brasileira para fins de posicionamento. Serão apresentados os modelos matemáticos associados a esses efeitos, as transformações envolvendo o campo magnético da Terra e a utilização do TEC advindo dos Mapas Globais da Ionosfera ou calculados a partir das observações GPS de pseudodistância. O processamento dos dados GPS foi realizado considerando o método relativo estático e cinemático e o posicionamento por ponto preciso (PPP). Os efeitos de segunda e terceira ordem foram analisados considerando períodos de alta e baixa atividade ionosférica. Os resultados mostraram que a não consideração desses efeitos no posicionamento por ponto preciso e no posicionamento relativo para linhas de base longas pode introduzir variações da ordem de poucos milímetros nas coordenadas das estações, além de variações diurnas em altitude da ordem de centímetros.

Período de mínima atividade solar: melhora no desempenho do posicionamento relativo

In the relative positioning, even considering that part of the errors due to ionosphere is canceled with the double-difference observations, strong ionospheric effects can occur in maximum solar activity period. However, in minimum solar activity period, the ionospheric effects decrease significantly and therefore an improvement of the relative positioning performance takes place. In this paper we aim at showing that improvement for the scientific and GPS community users. So, have been experiments by using GPS data of two stations of the Brazilian Network for Continuous Monitoring of GPS, forming a baseline of 430 km. The processing were use accomplished with interval of two hours, and only L1 carrier data have been used. The analysis of the obtained results has been carried out from the discrepancies between the "true" coordinates and corresponding ones obtained in the processing. In maximum solar activity period the discrepancy value reached 25 m. on the other hand, in minimum solar activity period, the discrepancy value reached 5,5 m. It is important to emphasize that the majority of the discrepancy values didn't exceed 0,50 m, and in some cases only reached 0,10 m. This shows the increase of application possibilities of the relative positioning using single-frequency GPS receivers in minimum solar activity period.

Desempenho do modelo global da ionosfera do IGS: avaliação no posicionamento por ponto na região sul do Brasil em período de alta atividade solar

One of the main drawbacks of the GPS accuracy for L1 users is the error due to ionosphere. This error depends on the total electron content presents in the ionosphere, as well as of the carrier frequency. Some models have been developed to correct GPS observables of the systematic error due to the ionosphere. The model more known and used is the Klobuchar model, which corrected 50-60% of the ionospheric error approximately. Alternatively, IGS (International GNSS Service) also has developed a model called Global Ionospheric Map (GIM). These maps, in format IONEX, are available in the site of the IGS, and one of the applications of them is to correct the GPS observables of the error due to ionosphere. This work aims at evaluating the quality of GPS point positioning using the IGS ionospheric model in the southerm region of Brazil. Tests carried out had shown an average improvement in the horizontal and vertical determination of 44% and 77%, respectively, when GIM is used in the point positioning.

Análise e comparação do comportamento da ionosfera e do posicionamento por ponto em períodos de alta e baixa atividade solar

The ionosphere is a major source of systematic error in the GPS observables. As this error is directly proportional to the TEC (Total Electron Content), the quality of GPS positioning (especially with single frequency receivers) can be significantly affected by regular changes of TEC. The ionosphere factor is even more relevant in the Brazilian region, where ionospheric phenomena, such as the Equatorial Anomaly, intensify these variations. Taking the above mentioned factors into account, experiments were conducted in this research to evaluate the daily and seasonal behavior of the TEC and the point positioning with GPS (single frequency) in periods of high and low solar activity in the Brazilian region. The results showed a direct correlation between the decrease in electrons density in the ionosphere (period of low solar activity) and improvement in positioning accuracy, as well as a large influence of Equatorial Anomaly on the results of point positioning.

Errors due to the tropospheric and ionospheric effects on the geographic location of data collection platforms

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Análise da ionosfera usando dados de receptores GPS durante um período de alta atividade solar e comparação com dados de Digissondas

Cada vez mais é crescente o uso do GPS (Global Positioning System ) em estudos da atmosfera terrestre. Neste artigo, a atmosfera superior da Terra, denominada ionosfera, foi estudada durante um período de alta atividade solar (ano de 2001) usando dados de receptores GPS de dupla freqüência localizados na região brasileira, pertencentes à RBMC (Rede Brasileira de Monitoramento Contínuo). A partir dos dados GPS foram calculados os valores de TEC (Total Electron Content - Conteúdo Total de Elétrons) da ionosfera. Adicionalmente foram incluídos no estudo dados de freqüência crítica da camada F2 advinda de Digissondas localizadas em São Luís/MA (3ºS; 44ºW) e Cachoeira Paulista/SP (22ºS; 45ºW), para fins de comparação. de uma forma geral, os resultados mostraram maiores valores de TEC durante os meses próximos aos equinócios e menores durante os meses de inverno. Quanto à variação diária do TEC, os menores valores foram verificados por volta das 4-6 HL (Hora Local) e os maiores durante o período da tarde, com valores um pouco maiores para São Luís. O segundo pico da anomalia equatorial foi verificado em Cachoeira Paulista nos meses próximos aos equinócios e verão. Na maioria dos meses, uma alta correlação linear foi verificada quando realizada a comparação entre os valores de freqüência crítica da camada F2 e os de TEC, principalmente para São Luís.

Unusual early morning development of the equatorial anomaly in the Brazilian sector during the Halloween magnetic storm

The solar events that occurred at the end of October 2003 gave rise to very strong geomagnetic disturbances that peaked twice with Dst values reaching -345 nT around 0000 UT on 30 October and -400 nT around 2300 UT, on the same day. Disturbances in several ionospheric parameters were observed over Brazil. This work will focus on the ionospheric response to the initial westward prompt penetration electric field and on the strong intensification of the equatorial ionization anomaly that occurred because of the electric field polarity reversal that followed in the early morning hours of 29 October. The F layer peak height over the equator first decreased under the strong prompt penetration westward electric field, which was followed by significant height increase under eastward electric field. We have used Sheffield University Plasmasphere Ionosphere Model (SUPIM) with an intensified westward disturbed electric field in the presunrise hours, presumably due to prompt penetration from the magnetosphere, in order to study the effect of such a field in the ionosphere. The simulation results showed that prompt penetration of magnetospheric electric fields of westward polarity to the nightside equatorial region seems to be the most probable cause of the initial F layer height decreases. The intensification of the equatorial ionization anomaly and the unusual enhancement on F layer peak density, which was not modeled by the SUPIM, are explained as caused by the strong eastward electric field that followed the initial phase in combination with a highly variable disturbed meridional/transequatorial wind system as inferred from the F2 layer peak height variations. The highly dynamic wind pattern, with a short-term response (2-4 hours), is compatible with the predictions of some previous theoretical model calculations reported in the literature.

Using cubic splines to mitigate systematic errors in GPS relative positioning

Systematic errors can have a significant effect on GPS observable. In medium and long baselines the major systematic error source are the ionosphere and troposphere refraction and the GPS satellites orbit errors. But, in short baselines, the multipath is more relevant. These errors degrade the accuracy of the positioning accomplished by GPS. So, this is a critical problem for high precision GPS positioning applications. Recently, a method has been suggested to mitigate these errors: the semiparametric model and the penalised least squares technique. It uses a natural cubic spline to model the errors as a function which varies smoothly in time. The systematic errors functions, ambiguities and station coordinates, are estimated simultaneously. As a result, the ambiguities and the station coordinates are estimated with better reliability and accuracy than the conventional least square method.