Quantitative evaluation of multipath rejection capabilities of GNSS antennas

Multipath error remains the largest error source in many high precision GPS applications. To counteract this problem, solutions at both software and hardware level have been studied. Software processing by means of measurement redundancy or error predictability can be used in order to mitigate the multipath effects. In general, these techniques work properly only when the length of a reflection path exceeds that of the direct path by more than 10–20 m. Unfortunately, in most cases, reflections are generated in the area near the receiving antenna. For this reason, multipath rejection actuated at the antenna level is one of the most valid means to improve the accuracy of GPS systems. The scope of this work is twofold. First, a review of low-multipath reception requirements will be proposed for comparing different classes of high precision GNSS antennas. Based on this discussion, we introduce a quantitative evaluation of multipath rejection capabilities of a GNSS antenna. The proposed assessment technique is focused on the antenna pattern, but, in contrast to other parameters evaluating the antenna radiation characteristics, it is specifically conceived to capture the effects of multipath signals.     

双极性天线在多路径效应探测中的应用

多路径效应是影响卫星定位精度和稳定性的重要因素。针对多路径环境复杂多变难以进行分析的问题,该文提出了一种利用双极性天线探测多路径的方法。详细介绍了反射信号的极性特征;通过判断信号的极性可以有效地区分直射信号和反射信号。分析了通过信号极性实现多路径探测的理论:组成双极性天线的右螺旋圆极化天线和左螺旋圆极化天线能够分别输出卫星信号中的右螺旋圆极化分量和左螺旋圆极化分量,如果两个分量的载噪比差值大于一定的门槛值,就可以推断卫星信号在传输过程中经历了多路径效应。最后通过仿真实验验证了该方法的可行性和有效性。     

Multipath mitigation by conventional antennas with ground planes and passive vertical structures

An antenna ground plane (GP) plays a major role in the mitigation of multipath coming from underneath the antenna. A distinct trend in user antenna size and weight reduction has been observed during recent years. With this focus in mind, an overview of conventional flat conductive and impedance GP is provided. Results of modern developments of passive vertically stacked structures are discussed. The possibility of obtaining a reasonable level of multipath mitigation with vertical structures of several centimeters is shown.     

Carrier phase multipath error characterization and reduction in single aircraft relative positioning

This paper describes GPS carrier phase (CP) multipath characterization and error reduction techniques and their application in single aircraft relative positioning. In particular, the single aircraft relative positioning scenario has applications for high-accuracy multi-sensor stabilization where CP multipath is a major error source that limits system performance. We will briefly review the requisite multipath theory and discuss models for quantifying the error characteristics. Field-test data will be used to validate the multipath models considering the underlying assumptions. A basic geometric reflection point theory will be presented to demonstrate the physical environmental sensitivity of CP multipath to parameters such as surface flatness and antenna height. Several different quantities will be described as multipath indicators for time-domain detection and will be compared with a frequency-domain technique. A new multipath detection approach will be introduced that is suitable to track multipath from time-varying reflection/diffraction points. Finally, the narrow-lane antenna baseline processing technique will be presented as a real-time approach to mitigate CP noise and multipath errors that is well suited for a very short baseline single aircraft high accuracy relative positioning system. Field-test results and analysis will be used to illustrate the key concepts in this paper and to help characterize the total navigation system performance.     

GPS snow depth meter with geometry-free linear combinations of carrier phases

Multipath in global positioning system (GPS) is the interference of the microwave signals directly from satellites and those reflected before reaching the antenna, typically by the ground. Because reflected signals cause positioning errors, GPS antennas are designed to reduce such interference. Recent studies show that multipath could be utilized to infer the properties of the ground around the antenna. Here, we report one such application, i.e. a fixed GPS station used as a snow depth meter. Because the satellite moves in the sky, the excess path length of reflected waves changes at rates dependent on the antenna height. This causes quasi-periodic variations of the amplitude and phase of the received signals. Accumulation of snow reduces effective antenna heights, and we can see it by analyzing multipath signatures. Signal-to-noise ratios (SNR) are often used to analyze multipath, but they are not always available in raw GPS data files. Here, we demonstrate that the geometry-free linear combination (L4), normally used to study the ionosphere, can also be used to analyze multipath signatures. We obtained snow depth time series at a GPS station in Hokkaido, Japan, from January to April in 2009 using L4 and SNR. Then, we compared their precisions. We also discuss mechanisms responsible for the possible underestimation of the snow depth by GPS. Finally, we investigate the possibility of inferring physical conditions of the snow surface using amplitudes of multipath signatures.     

Effects of azimuthal multipath asymmetry on long GPS coordinate time series

Carrier phase multipath is currently one source of unmodeled signals that may bias GPS coordinate time series significantly. We investigate the effect of simulated carrier phase multipath on time series of several sites covering the period 2002.0–2008.0 and spanning a range of observation geometries. High-, mid-, and low-latitude IGS sites are investigated as well as sites with significant signal obstructions. We examine the effect of multipath in different sectors of the sky, considering time-constant, horizontal reflectors at each of 0.1, 0.2, and 1.5 m below the antenna. The differences between a horizontally uniform multipath source are analyzed, and it is shown that positioning errors are generally larger when unmodeled carrier phase multipath is azimuthally heterogeneous. Using the adopted multipath model, height biases reach ±1 mm in case of the symmetric multipath and ±5 mm for the asymmetric multipath but this increases to being ±10 mm in the worst case. In addition to mean bias, low-frequency variations in the bias also exist, including periodic signals and leading to velocity biases of up to ±0.1 mm/year in the symmetric case and ±1 mm/year in the asymmetric case over the considered period. In contrast to the generally slowly varying observation geometry that is typically experienced, we show the effects of an abrupt change in geometry due to receiver/antenna hardware changes; in the case considered, we see changed pattern of temporal variation in the bias in addition to an instantaneous offset.     

Evaluating pseudorange multipath effects at stations in the National CORS Network

A preliminary study was conducted to evaluate the amount of pseudorange multipath at 390+ sites in the National Continuously Operating Reference Station (CORS) Network. The National CORS Network is a cooperative effort involving over 110 different agencies, universities, and private companies who seek to make GPS data from dual-frequency receivers located throughout the United States and its territories available to the general public. For CORS users, pseudorange multipath can seriously degrade the accuracy of any application that relies on precise measurements of the pseudorange observable over a short period of time, including differential pseudorange navigation, kinematic and rapid-static surveying, and ionospheric monitoring. The main objectives of this study were to identify the most affected and least affected sites in the network, to closely investigate problematic sites, and to compare various receiver/antenna combinations. Dual-frequency carrier phase and pseudorange measurements were used to estimate the amount of L1 and L2 pseudorange multipath at each site over a one-year period. Some of the most severely affected sites were maritime Differential GPS and Nationwide Differential GPS (DGPS/NDGPS) sites. Photographs obtained for these sites verified the presence of transmission towers and other reflectors in close proximity to the GPS antennas. Plotting the variations of the L1 and L2 pseudorange multipath with respect to azimuth and elevation further verified that even above a 60° elevation angle there was still as much as five meters of pseudorange multipath at some sites. The least affected sites were the state networks installed in Ohio and Michigan; these sites used excellent antenna mounts, choke ring antennas, and new receiver technology. A comparison of the 12 most commonly used receiver/antenna combinations in the CORS Network indicated that newer receivers such as the Ashtech UZ-12, Leica RS-500, and Trimble 5700 help to significantly mitigate pseudorange multipath, while the receivers/antennas at some DGPS/NDGPS sites, and the antennas formerly used at the Wide Area Augmentation System (WAAS) sites, are among those most affected by pseudorange multipath. The receiver/antenna comparison did not take into account the potential presence of reflectors at the sites (i.e., it is possible that a well-performing receiver/antenna combination could have been consistently placed at very poor site locations, and vice-versa).Product Disclaimer: Mention of a commercial company or product does not constitute an endorsement by the National Oceanic and Atmospheric Administration. Use for publicity or advertisement purposes of information from this paper concerning proprietary products or the comparison of such products is not authorized.     

Development and testing of a new ray-tracing approach to GNSS carrier-phase multipath modelling

Multipath is one of the most important error sources in Global Navigation Satellite System (GNSS) carrier-phase-based precise relative positioning. Its theoretical maximum is a quarter of the carrier wavelength (about 4.8 cm for the Global Positioning System (GPS) L1 carrier) and, although it rarely reaches this size, it must clearly be mitigated if millimetre-accuracy positioning is to be achieved. In most static applications, this may be accomplished by averaging over a sufficiently long period of observation, but in kinematic applications, a modelling approach must be used. This paper is concerned with one such approach: the use of ray-tracing to reconstruct the error and therefore remove it. In order to apply such an approach, it is necessary to have a detailed understanding of the signal transmitted from the satellite, the reflection process, the antenna characteristics and the way that the reflected and direct signal are processed within the receiver. This paper reviews all of these and introduces a formal ray-tracing method for multipath estimation based on precise knowledge of the satellite–reflector–antenna geometry and of the reflector material and antenna characteristics. It is validated experimentally using GPS signals reflected from metal, water and a brick building, and is shown to be able to model most of the main multipath characteristics. The method will have important practical applications for correcting for multipath in well-constrained environments (such as at base stations for local area GPS networks, at International GNSS Service (IGS) reference stations, and on spacecraft), and it can be used to simulate realistic multipath errors for various performance analyses in high-precision positioning.     

An open source GPS multipath simulator in Matlab/Octave

Multipath is detrimental for both GPS positioning and timing applications. However, the benefits of GPS multipath for reflectometry have become increasingly clear for monitoring soil moisture, snow depth, and vegetation growth. In positioning applications, a simulator can support multipath mitigation efforts in terms of, e.g., site selection, antenna design, receiver performance assessment, and in relating different observations to a common parameterization. For reflectometry, in order to convert observed multipath parameters into useable environmental products, it is important to be able to explicitly link the GPS observables to known characteristics of the GPS receiver/antenna and the reflecting environment. Existing GPS multipath software simulators are generally not readily available for the general scientific community to use and/or modify. Here, a simulator has been implemented in Matlab/Octave and is made available as open source code. It can produce signal-to-noise ratio, carrier phase, and code pseudorange observables, based on L1 and L2 carrier frequencies and C/A, P(Y), and L2C modulations. It couples different surface and antenna types with due consideration for polarization and coherence. In addition to offering predefined material types (water, concrete, soil, etc.), it allows certain dimensional properties to be varied, such as soil moisture and snow density.     

Characterization of multipath phase rates in different multipath environments

A well-known effect of multipath propagation is multipath fading that typically causes periodic signal variations. Such signal variations may become visible in some basic GNSS observables such as the code minus carrier observable, single or double differences or in C/N ₀ time series. The frequency of these variations—also called fading frequency or multipath phase rate—strongly depends on the multipath environment, i.e. on the actual geometric conditions which can be described by the location of the satellite causing the multipath signal and the reflector location with respect to the receiving antenna. This paper gives a detailed insight on the expected multipath phase rates in different multipath environments. Different geometric conditions are analyzed, from arbitrary reflector positions to the point of dealing with the special case of ground multipath. Fading frequencies are determined by means of an empirical approach using the characteristics of real satellite passes. The approach results in distributions of multipath phase rates which are computed for a multitude of possible reflector locations and from which minimum, mean and maximum multipath phase rates can be derived.     

A Methodology for Monitoring Vertical Dynamic Sub-Centimeter Displacments with GPS

This article describes a methodology to monitor dynamic vertical sub-centimeter displacement, of a GPS antenna. The dynamic movement of an antenna is determined by choosing the appropriate reference satellite for measurement differencing and by applying a FFT filter on the double-difference phase residuals. The validity of the method depends on the time variations of the GPS residuals and errors, such as, receiver noise, atmospheric contribution, multipath effects, and the antenna movement. This research is under development and results for simulated motion are presented here. ? 2002 Wiley Periodicals, Inc.     

Forward modeling of GPS multipath for near-surface reflectometry and positioning applications

Multipath is detrimental for both GPS positioning and timing applications. However, the benefits of GPS multipath for reflectometry have become increasingly clear for soil moisture, snow depth, and vegetation growth monitoring. Most multipath forward models focus on the code modulation, adopting arbitrary values for the reflection power, phase, and delay, or they calculate the reflection delay based on a given geometry and keep reflection power empirically defined. Here, a fully polarimetric forward model is presented, accounting for right- and left-handed circularly polarized components of the GPS broadcast signal and of the antenna and surface responses as well. Starting from the fundamental direct and reflected voltages, we have defined the interferometric and error voltages, which are of more interest in reflectometry and positioning applications. We examined the effect of varying coherence on signal-to-noise ratio, carrier phase, and code pseudorange observables. The main features of the forward model are subsequently illustrated as they relate to the broadcast signal, reflector height, random surface roughness, surface material, antenna pattern, and antenna orientation. We demonstrated how the antenna orientation—upright, tipped, or upside-down—involves a number of trade-offs regarding the neglect of the antenna gain pattern, the minimization of CDMA self-interference, and the maximization of the number of satellites visible. The forward model was also used to understand the multipath signature in GPS positioning applications. For example, we have shown how geodetic GPS antennas offer little impediment for the intake of near-grazing reflections off natural surfaces, in contrast to off metal, because of the lack of diversity with respect to the direct signal—small interferometric delay and Doppler, same sense of polarization, and similar direction of arrival.     

NLOS GPS signal detection using a dual-polarisation antenna

The reception of indirect signals, either in the form of non-line-of-sight (NLOS) reception or multipath interference, is a major cause of GNSS position errors in urban environments. We explore the potential of using dual-polarisation antenna technology for detecting and mitigating the reception of NLOS signals and severe multipath interference. The new technique computes the value of the carrier-power-to-noise-density (C/N ₀) measurements from left-hand circular polarised outputs subtracted from the right-hand circular polarised C/N ₀ counterpart. If this quality is negative, NLOS signal reception is assumed. If the C/N ₀ difference is positive, but falls below a threshold based on its lower bound in an open-sky environment, severe multipath interference is assumed. Results from two experiments are presented. Open-field testing was first performed to characterise the antenna behaviour and determine a suitable multipath detection threshold. The techniques were then tested in a dense urban area. Using the new method, two signals in the urban data were identified as NLOS-only reception during the occupation period at one station, while the majority of the remaining signals present were subject to a mixture of NLOS reception and severe multipath interference. The point positioning results were dramatically improved by excluding the detected NLOS measurements. The new technique is suited to a wide range of static ground applications based on our results.     

A multi-technique approach for characterizing the SVN49 signal anomaly, part 2: chip shape analysis

Due to a satellite internal reflection at the L5 test payload, the SVN49 (PRN1) GPS satellite exhibits a static multipath on the L1 and L2 signals, which results in elevation-dependent tracking errors for terrestrial receivers. Using a 30-m high-gain antenna, code and carrier phase measurements as well as raw in-phase and quadrature radio frequency samples have been collected during a series of zenith passes in mid-April 2010 to characterize the SVN49 multipath and its impact on common users. Following an analysis of the receiver tracking data and the IQ constellation provided in Part 1 of this study, the present Part 2 provides an in-depth investigation into chip shapes for the L1 and L2 signals. A single reflection model is found to be compatible with the observed chip shape distortions and key parameters for an elevation dependent multipath model are derived. A good agreement is found between multipath parameters derived independently from raw IQ-samples and measurements of a so-called Vision Correlator. The chip shapes and their observed variation with elevation can be used to predict the multipath response of different correlator types within a tracking receiver. The multipath model itself is suitable for implementation in a signal simulator and thus enables laboratory testing of actual receiver hardware.     

A multi-technique approach for characterizing the SVN49 signal anomaly, part 1: receiver tracking and IQ constellation

A characterization of the signal anomaly of SVN49 is presented. A mathematical model is developed to relate the observed multipath to the internal signal reflection. The analyses provided are based on measurements, which have been collected during a dedicated tracking campaign with a 30-m dish antenna. Data on the L1 and L2 frequency have been collected with four different receivers. In addition, IQ samples have been recorded directly with a spectrum analyzer. The multipath combination of the receiver measurements on L1 and L2 is analyzed to demonstrate the effect of the signal reflections on different correlator spacing. The capability to suppress the signal reflection with receiver multipath mitigation methods is demonstrated. Finally, preliminary estimates of the attenuation, delay, and phase shift over elevation are obtained from an IQ sample analysis.     

Analysis of high-frequency multipath in 1-Hz GPS kinematic solutions

High-frequency multipath would be problematic for studies at seismic or antenna dynamical frequencies as one could mistakenly interpret them as signals. A simple procedure to identify high-frequency multipath from global positioning system (GPS) time series records is presented. For this purpose, data from four GPS base stations are analyzed using spectral analyses techniques. Additional data, such as TEQC report files of L1 pseudorange multipath, are also used to analyze the high-frequency multipath and confirmation of the high-frequency multipath inferred from the phase records. Results show that this simple procedure is effective in identification of high-frequency multipath. The inferred information can aid interpretation of multipath at the GPS site, and is important for a number of reasons. For example, the information can be used to study GPS site selections and/or installations.

Multipath mitigation via component analysis methods for GPS dynamic deformation monitoring

Multipath is one of the main error sources in high-precision global positioning system (GPS) dynamic deformation monitoring, as it is difficult to be mitigated by differencing between observations. In addition, since a specific frequency threshold value between multipath and deformation signals may not exist, multipath is usually inseparable from the low-frequency vibration signal using conventional frequency-domain filter methods. However, the multipath repeats in two sidereal days when the surroundings of a GPS antenna remain unchanged. This characteristic can be exploited to model and thus mitigate multipath effectively in dynamic deformation monitoring. Unfortunately, a major issue is that the degree of repeatability decreases as the interval between first day and subsequent days increases. To overcome this problem, we develop a new sidereal filtering referred to as reference EMD-ICA (EMD-ICA-R), where empirical mode decomposition (EMD) and independent component analysis (ICA) are jointly used to model multipath and renew the reference multipath. For the successful implementation of the EMD-ICA-R, an a priori denoised multipath signal is needed as a reference. We further propose to use the principal component analysis (PCA) method to extract more accurate reference multipath signal and form a combined PCA-EMD-ICA-R approach. Simulation experiments with a motion simulation platform were conducted, and the testing results indicate that the proposed methods can mitigate the multipath by around 67 % when a reliable reference multipath signal is extracted from a static situation. Furthermore, simulation experiments with different deformation signals added into the coordinate time series of three consecutive days show that the two proposed methods are also effective in a dynamic situation. Since wavelet filtering is used to denoise the reference multipath signals in the new approaches, simulation experiments with several wavelet filters are tested, and the results indicate that the PCA-EMD-ICA-R approach can work well with various wavelet filters.     

Site-specific multipath characteristics of global IGS and CORS GPS sites

The site-specific multipath characteristics of 217 Global Positioning System (GPS) sites worldwide were analyzed using the variability of the post-fit phase residuals. Because the GPS satellite constellation returns to the same configuration in a sidereal day (23 h 56 min 4 s), the multipath repeats on that period. However, daily GPS position estimates are usually based on the solar day. When several days of GPS data are processed, this steady change in the orientation of the satellite constellation with respect to the station manifests itself in the form of patterns in the post-fit phase residuals which shift by 3 min 56 s per day. It was found that the mean root mean square of the time-shifted post-fit phase residuals is highly dependent on the GPS antenna type. The conclusions derived from the analysis of the time-shifted post-fit residuals were verified by performing a cross-correlation of the post-fit residuals across many days for selected sites.     

Tests of phase center variations of various GPS antennas, and some results

Phase variations of GPS receiving antennas are a significant error component in precise GPS applications. A calibration procedure has been developed by Geo++ and the Institut für Erdmessung, which directly determines absolute phase center variations (PCVs) without any multipath influence by field measurements. The precision and resolution of the procedure allows the determination of reliable azimuthal variations. PCV may affect long-term static GPS differently than real-time GPS, depending on the applications. At the same time, different antenna types are involved. Less investigations have been done on absolute PCV of rover antennas than on geodetic antennas which, however, becomes more important due to the mixed antenna situation in GPS reference networks and RTK networks. The concepts of the absolute PCV field calibration are summarized and emphasis is placed on a variety of absolute PCV patterns of geodetic and rover antennas. Electronic Publication     

Correction of apparent position shifts caused by GNSS antenna changes

Antenna changes at GNSS reference stations frequently produce discontinuities in the coordinate time series. These apparent position shifts are mainly caused by changes of carrier-phase multipath effects and different errors in the antenna phase center corrections. A monitoring method was developed and successfully tested, which requires additional GNSS observations from a local, temporary reference station. Changes of carrier-phase measurement errors due to the antenna change are determined and stored in L1 and L2 phase maps. These phase maps provide corrections to be applied either to the observation data obtained before the antenna change or to the observation data obtained after the antenna change. The observation corrections are able to remove coordinate discontinuities independent of the selected coordinate estimation algorithm.