The wavelet method as an alternative for reducing ionospheric effects from single-frequency GPS receivers

The ionospheric effect is one of the major errors in GPS data processing over long baselines. As a dispersive medium, it is possible to compute its influence on the GPS signal with the ionosphere-free linear combination of L1 and L2 observables, requiring dual-frequency receivers. In the case of single-frequency receivers, ionospheric effects are either neglected or reduced by using a model. In this paper, an alternative for single-frequency users is proposed. It involves multiresolution analysis (MRA) using a wavelet analysis of the double-difference observations to remove the short- and medium-scale ionosphere variations and disturbances, as well as some minor tropospheric effects. Experiments were carried out over three baseline lengths from 50 to 450 km, and the results provided by the proposed method were better than those from dual-frequency receivers. The horizontal root mean square was of about 0.28 m (1 sigma).

Integration of a regional GPS network within ITRF using precise point positioning

GPS precise point positioning (PPP) can provide high precision 3-D coordinates. Combined pseudorange and carrier phase observables, precise ephemeris and satellite clock corrections, together with data from dual frequency receivers, are the key factors for providing such levels of precision (few centimeters). In general, results obtained from PPP are referenced to an arbitrary reference frame, realized from a previous free network adjustment, in which satellite state vectors, station coordinates and other biases are estimated together. In order to obtain consistent results, the coordinates have to be transformed to the relevant reference frame and the appropriate daily transformation parameters must be available. Furthermore, the coordinates have to be mapped to a chosen reference epoch. If a velocity field is not available, an appropriated model, such as NNR-NUVEL-IA, has to be used. The quality of the results provided by this approach was evaluated using data from the Brazilian Network for Continuous Monitoring of the Global Positioning System (RBMC), which was processed using GIPSY-OASIS 11 software. The results obtained were compared to SIRGAS 1995.4 and ITRF2000, and reached precision better than 2cm. A description of the fundamentals of the PPP approach and its application in the integration of regional GPS networks with ITRF is the main purpose of this paper.

Integridade no posicionamento RTK e RTK em rede

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

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.

Transformação de Helmert generalizada no posicionamento de alta precisão: fundamentação teórica e exemplificações

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.

Spacecraft's attitude prediction: solar radiation torque and the Earth's shadow

A semi-analytical approach is proposed to study the rotational motion of an artificial satellite, under the influence of torque due to the solar radiation pressure, and taking into account the influence of Earth's shadow. Using Andoyer variables the equations for the rotational motion are presented in extended Hamiltonian form. In order to get a solution for the state variables close to an actual motion, the considered model for the shadow function takes into account physical and geometric factors and three specific regions: shadow, penumbra and full light. A mapping for the shadow function is proposed and a semi-analytical process is applied. When the satellite is totally illuminated or it is inside the penumbra, a known analytical solution is used to compute the satellite's attitude. A numerical simulation shows, when the penumbra region is included, the attenuation of the rotational motion during the transition from the shadow to the illuminate region and vice versa. (c) 2005 Published by Elsevier Ltd on behalf of COSPAR.

Analytical attitude prediction of spin stabilized spacecrafts perturbed by magnetic residual torque

An analytical approach for spin-stabilized spacecraft attitude prediction is presented for the influence of the residual magnetic torques. Assuming an inclined dipole model for the Earth's magnetic field, an analytical averaging method is applied to obtain the mean residual torque every orbital period. The orbit mean anomaly is utilized to compute the average components of residual torque in the spacecraft body frame reference system. The theory is developed for time variations in the orbital elements, and non-circular orbits, giving rise to many curvature integrals. It is observed that the residual magnetic torque does not have component along the spin axis. The inclusion of this torque on the rotational motion differential equations of a spin stabilized spacecraft yields conditions to derive an analytical solution. The solution shows that residual torque does not affect the spin velocity magnitude, contributing only for the precession and the drift of the spin axis of the spacecraft. (c) 2005 COSPAR. Published by Elsevier Ltd. All rights reserved.

Inner satellites of Neptune: I the disturbing function

In this work, we focus our attention to the expansion of the disturbing function (R) which governs the dynamics of a satellite (natural or artificial) in the Neptune-Triton system. What makes this problem quite unusual, is the fact that a small inner satellite can be strongly disturbed by Triton which is moving in a highly inclined and retrograde orbit. These features are unique in our solar system. Although a lot of retrograde satellites are currently known, all of them have negligible mass and the), do not offer almost any perturbation on the others satellites. However, in the case of the inner satellites of Neptune, Triton is an interesting exception. In a highly inclined orbit, the perturbation it exerts on the neighbouring satellites of Neptune cannot be ignored even for the present scenario. However, in the future, this perturbation will be much more important because due to the tides, the orbit of Triton is contracting, whereas the semi major axes of the remaining inner satellites of Neptune will remain almost unaffected by the tides. In this work we first obtain the disturbing function in the retrograde case. After that, we generalize R for arbitrary inclination. Several numerical tests are presented and a possible future case of resonant configuration is briefly discussed as well. (c) 2005 COSPAR. Published by Elsevier Ltd. All rights reserved.

Orbital propagation for Brazilian satellites using NORAD models

The performance of the NORAD models for near Earth satellites (SGP, SGP4, SGP8) using two Brazilian flying satellites, SCD-1 and SCD-2, and the Chinese-Brazilian satellite CBERS-1 is compared. The achievable accuracy of such models is compared against the predicted 2-lines elements set for the satellites. Every week an updated fresh set of 2-lines elements for these satellites is made available through the Internet. About ten years of observations of the SCD-1 satellite are available and therefore solar activity influences on the orbit can be analyzed. Data are selected considering also orbital (for CBERS-1) and attitude (for SCD-2) maneuvers - (C) 2002 COSPAR. Published by Elsevier B.V. Ltd. All rights reserved.

Optimal space manoeuvres in a non-Keplerian force field

In this work the problem of a spacecraft bi-impulsive transfer between two given non coplanar elliptical orbits, with minimum fuel consumption, is solved considering a non-Keplerian force field (the perturbing forces include Earth gravity harmonics and atmospheric drag). The problem is transformed in the Two Point Boundary Value Problem. It is developed and implemented a new algorithm, that uses the analytical expressions developed here. A dynamics that considered a Keplerian force field was used to produce an initial guess to solve the Two Point Boundary Value Problem. Several simulations were performed to observe the spacecraft orbital behaviour by different kind of perturbations and constraints, on a fuel consumption optimization point of view. (C) 2002 COSPAR. Published by Elsevier B.V. Ltd. All rights reserved.

Numerical study of low-cost alternative orbits around the Moon

In this paper, we have investigated a region of direct stable orbits around the Moon, whose stability is related to the H2 Family of periodic orbits and to the quasi-periodic orbits that oscillate around them. The stability criteria adopted was that the path did not escape from the Moon during an integration period of 1000 days (remaining with negative two-body Moon-probe orbital energy during this period). Considering the three-dimensional four-body Sun-Earth-Moon-probe problem, we investigated the evolution of the size of the stability region, taking into account the eccentricity of the Earth's orbit, the eccentricity and inclination of the Moon's orbit, and the solar radiation pressure on the probe. We also investigated the evolution of the region's size and its location by varying the inclination of the probe's initial osculating orbit relative to the Moon's orbital plane between 0 degrees and 180 degrees. The size of the stability region diminishes; nevertheless, it remains significant for 0 <= i <= 25 degrees and 35 degrees <= i <= 45 degrees. The orbits of this region could be useful for missions by space vehicles that must remain in orbit around the Moon for periods of up to 1000 days, requiring low maintenance costs. (c) 2005 COSPAR. Published by Elsevier Ltd. All rights reserved.

Multipath reduction from GPS double differences using wavelets: How far can we go?

Wavelets are being extensively used in Geodetic applications. In this paper, the Multi-Resolution Analysis (MRA) using wavelets is applied to pseudorange and carrier phase GPS double differences (DDs) in order to reduce multipath effects. The wavelets were already applied to GPS carrier phase DDs, but some questions remain: How good can be the results, and are all multipath effects reduced? The answers to these questions are discussed in this paper. Thus, the wavelet transform is used to decompose the DD signals, splitting them in lower resolution components. After the decomposition process, the wavelet shrinkage is performed by thresholding to eliminate the components relative to multipath effects. Then, the DD observation can be reconstructed. This new DD signal is used to perform the baseline processing. The daily multipath repeatability was verified. With the application of the proposed approach, the results showed that the reliability of the ambiguity resolution and accuracy of the results improved when compared with the standard procedure. Furthermore, the method showed to be very efficient computationally, because, it is not noticed, at practical level, difference in the time span between the processing with and without application of the proposed method. However, only the high frequency multipath was eliminated.

Conseqüências de uma tempestade geomagnética no posicionamento relativo com receptores GPS de simples freqüência

The error associated with the ionosphere depends on Total Electron Content (TEC) of the ionosphere. The geomagnetic field exerts strong influence in the TEC variation, because it controls the movement of the electrons. After solar events the magnetic lines of force can be compressed, characterizing the geomagnetic storm. The aim of this paper is to present to geodesic community the effects of a geomagnetic storm in the relative positioning. The processing of the data was accomplished with an interval of two hours, with a 430 km baseline. The analyze of the obtained results have been carried out from the discrepancies between the true coordinates and corresponding ones obtained in the processing of the baseline. The used data in this paper include the period of 30/03/2001 up to 02/04/2001. In March 31 a strong geomagnetic storm happened. One day after, that it corresponds to main phase of the storm, the values of the discrepancies decreased significantly. For instance, in 01:00-03:00 UT period, the value of the planimetric discrepancy reached 20 m in the storm day. However, in the main phase of the storm, the planimetric discrepancy decreased to 0.1 m.