A new coarse-time GPS positioning algorithm using combined Doppler and code-phase measurements

A new coarse-time Global Positioning System (GPS) positioning algorithm based on the use of Doppler and code-phase measurements is proposed and described. The proposed method was demonstrated to be essential for reducing the time to first fix and the power consumption in a GPS receiver. Only 1 ms of data is required to obtain a positioning fix with accuracy comparable to that of the traditional GPS navigation algorithm. The algorithm is divided into two parts. In the first part, the Doppler measurement of the GPS signal is used to determine the coarse user position. With proper constraints, the required time accuracy for the Doppler measurements can be relaxed to be as long as 12 h. In the second part of the algorithm, the accurate user position is calculated by means of the 1 ms code-phase data. The traditional tracking process is no longer necessary in the proposed algorithm. Using the acquired 1-ms code-phase measurement, the positioning accuracy was determined to be approximately a few tens of meters in our experimental results. However, if the data length is extended to 10 ms, the positioning accuracy can be improved to within 10–20 m, which is similar to that of the traditional GPS positioning method. Various experiments were conducted to verify the usefulness of the proposed algorithm.     

A note on observation decorrelation,variances of residuals,and redundancy numbers

A technique for processing correlated observations suitable for large, sparse, least-squares adjustments is reviewed. Correlated coordinate differences derived from the Global Positioning System are used as illustrative examples. However, the methods examined are suitable for all types of correlated observations. The computation of variances of residuals, redundancy numbers, and marginally detectable errors is considered for sparse systems.     

Combination of GPS and VLBI on the observation level during CONT11—common parameters,ties and inter-technique biases

Multi-technique space geodetic analysis software has been developed which allows to combine data on the observation level. In addition to local tie information, site-wise common parameters, i.e., troposphere and clocks, can be estimated with this software. Thus, it will be discussed how common parameters have to be estimated and where biases/offsets need to be taken into account. To test such a novel concept, Global Positioning System (GPS) and Very Long Baseline Interferometry (VLBI) data from the CONT11 campaign are being utilized. Since the VLBI baselines of this campaign extend over several thousands of kilometers, GPS data are processed in precise-point positioning mode and satellite orbits and clocks are kept fixed to the IGS final products. From the obtained results, it can be shown that the combination of space geodetic data on the observation level leads to a consistent improvement of station position repeatability as well as nuisance parameters like troposphere estimates. Furthermore, estimation of common parameters (troposphere or clocks) at co-located sites helps to improve the solution further and derive an utmost physically consistent model of the concerned parameters.     

Non-linear motions of Australian geodetic stations induced by non-tidal ocean loading and the passage of tropical cyclones

We investigate daily and sub-daily non-tidal oceanic and atmospheric loading (NTOAL) in the Australian region and put an upper bound on potential site motion examining the effects of tropical cyclone Yasi that crossed the Australian coast in January/February 2011. The dynamic nature of the ocean is important, particularly for northern Australia where the long-term scatter due to daily and sub-daily oceanic changes increases by 20–55 % compared to that estimated using the inverted barometer (IB) assumption. Correcting the daily Global Positioning System (GPS) time series for NTOAL employing either a dynamic ocean model or the IB assumption leads to a reduction of up to 52 % in the weighted scatter of daily coordinate estimates. Differences between the approaches are obscured by seasonal variations in the GPS precision along the northern coast. Two compensating signals during the cyclone require modelling at high spatial and temporal resolution: uplift induced by the atmospheric depression, and subsidence induced by storm surge. The latter dominates ( \(>\) 135 %) the combined net effect that reaches a maximum of 14 mm, and 10 mm near the closest GPS site TOW2. Here, 96 % of the displacement is reached within 15 h due to the rapid transit of cyclones and the quasi-linear nature of the coastline. Consequently, estimating sub-daily NTOAL is necessary to properly account for such a signal that can be 3.5 times larger than its daily-averaged value. We were unable to detect the deformation signal in 2-hourly GPS processing and show that seasonal noise in the Austral summer dominates and precludes GPS detection of the cyclone-related subsidence.     

Improvement of the GPS time comparisons by simultaneous relative positioning of the receiver antennas

The Global Positioning System,GPS, is widely used for time comparisons between distant laboratories. Over distances of the order of 1000km or less, the system has the capability of 1 to 2ns accuracy. However this requires a relative positioning with errors lower than 30cm. We show that this positioning can be derived from theGPS time comparisons themselves. An example for European laboratories is given.     

A demonstration of 1–2 parts in 107 accuracy using GPS

By interferometric analysis ofGPS phase observations made at Owens Valley, Mojave, and Mammoth Lakes, California, we determined the coordinate components of the71–245–313 km triangle of baselines connecting these sites. A separate determination was made on each of four days, April 1–4, 1985. The satellite ephemerides used in these determinations had been derived from observations on other baselines. Therms scatters of the four daily determinations of baseline vector components about their respective means ranged from a minimum of6 mm for the north component of the71-km baseline to a maximum of34 mm for the vertical component of the245-km baseline. To test accuracy, we compared the mean of ourGPS determinations of the245-km baseline between Owens Valley and Mojave with independent determinations by others using very-long-baseline interferometry(VLBI) and satellite laser ranging(SLR). TheGPS-VLBI difference was within 2 parts in10 ⁷ for every vector component. TheGPS-SLR difference was within6 parts in10 ⁸ in the horizontal coordinates, but83 mm in height.     

GPS inferred geocentric reference frame for satellite positioning and navigation

Accurate geocentric three dimensional positioning is of great importance for various geodetic and oceanographic applications. While relative positioning accuracy of a few centimeters has become a reality using Very Long Baseline Interferometry (VLBI), the uncertainty in the offset of the adopted coordinate system origin from the geocenter is still believed to be of the order of one meter. Satellite Laser Ranging (SLR) is capable of determining this offset to better than10 cm, though, because of the limited number of satellites, this requires a long arc of data. The Global Positioning System (GPS) measurements provide a powerful alternative for an accurate determination of this origin offset in relatively short period of time. Two strategies are discussed, the first utilizes the precise relative positions predetermined byVLBI, where as the second establishes a reference frame by holding only one of the tracking sites longitude fixed. Covariance analysis studies indicate that geocentric positioning to an accuracy of a few centimeters can be achieved with just one day of preciseGPS pseudorange and carrier phase data.     

Improved method for estimating the inter-frequency satellite clock bias of triple-frequency GPS

Considering the contribution of the hardware biases to the estimated clock errors, an improved method for estimating the satellite inter-frequency clock bias (IFCB) is presented, i.e., the difference in the satellite clock error as computed from ionospheric-free pseudorange and carrier phase observations using L1/L2 and P1/P2 versus L1/L5 and P1/P5. The IFCB is composed of a constant and a variable part. The constant part is the inter-frequency hardware bias (IFHB). It contains the satellite and receiver hardware delays and can be expressed as a function of the DCBs [DCB (P1 ? P2) and DCB (P1 ? P5)]. When a reference satellite is selected, the satellite IFHB can be computed but is biased by a reference satellite IFHB. This bias will not affect the utilization of IFCB in positioning since it can be absorbed by the receiver clock error. Triple-frequency observations of 30 IGS stations between June 1, 2013, and May 31, 2014, were processed to show the variations of the IFHB. The IFHB values show a long-term variation with time. When a linear and a fourth-order harmonic function are used to model the estimated IFCB, which contains contributions of the hardware delays and clock errors, the results show that 89 % of the IFCB can be corrected given the current five triple-frequency GPS satellites with the averaged fitting RMS of 1.35 cm. Five days of data are processed to test the estimated satellite clock errors using the strategy presented. The residuals of P1/P5 and L1/L5 have a STD of <0.27 m and 0.97 cm, respectively. In addition, most predicted satellite IFCBs reach an accuracy of centimeter level and its mean accuracy of 5 days is better than 7 cm.     

Correlation between ROTI and Ionospheric Scintillation Indices using Hong Kong low-latitude GPS data

The correlation between the rate of TEC index (ROTI) and scintillation indices S ₄ and σ _Φ for low-latitude region is analyzed in this study, using data collected from a Global Positioning System (GPS) scintillation monitoring receiver installed at the south of Hong Kong for the periods June–August of 2012 and May 2013 and July–December of 2013. The analysis indicates that the correlation coefficient between ROTI and S ₄/σ _Φ is about 0.6 if data from all GPS satellites are used together. If each individual satellite is considered, the correlation coefficients are above 0.6 on average and sometimes above 0.8. The analysis also shows that the ratio of ROTI and S ₄ varies between 1 and 4. The ratio ROTI/σ _Φ, varies between 2 and 9. In addition, it is also found that there is a good consistency between the temporal variations of ROTI with scintillation activity under different ionospheric conditions. ROTI has a high correlation relationship with scintillation indices on geomagnetically disturbed days or in solar active months. Moreover, the data observed at low elevation angles have weak correlation between ROTI and scintillation indices. These results demonstrate the feasibility of using ROTI derived from GPS observations recorded by common non-scintillation GPS receivers to characterize ionospheric scintillations.     

A New Approach to 3D Dense LiDAR Data Classification in Urban Environment

Classification of Mobile Mapping LiDAR (Light Detection and Ranging) data is a challenge in the research community since the day when laser scanner system were integrated and mounted on vehicles for collection of 3D data in urban environment. The approach proposed here for classifying LiDAR data is analogous to the process followed for classifying data from satellite images. Pixel based and segmentation based methods have been employed in past for classifying images obtained from satellites. These methods were based on spectral properties of objects present in the images. But for Mobile mapping LiDAR data this approach has been applied and tested for the first time. The properties of this data are completely different from that of satellite images. So even if the basic approach remains the same, many changes have to be made in the entire classification process. The paper here aims to propose the basic procedure of using pixel-wise classification on dense 3D LiDAR data.     

Demonstration of sub-meter GPS orbit determination and 1.5 parts in 108 three-dimensional baseline accuracy

Differential tracking of theGPS satellites in high-earth orbit provides a powerful relative positioning capability, even when a relatively small continental U.S. fiducial tracking network is used with less than one-third of the fullGPS constellation. To demonstrate this capability, we have determined baselines of up to2000 km in North America by estimating high-accuracyGPS orbits and ground receiver positions simultaneously. The2000 km baselines agree with very long baseline interferometry(VLBI) solutions at the level of1.5 parts in10 ⁸ and showrms daily repeatability of0.3–2 parts in10 ⁸. The orbits determined for the most thoroughly trackedGPS satellites are accurate to better than1 m. GPS orbit accuracy was assessed from orbit predictions, comparisons with independent data sets, and the accuracy of the continental baselines determined along with the orbits. The bestGPS orbit strategies included data arcs of at least one week, process noise models for tropospheric fluctuations, estimation ofGPS solar pressure coefficients, and combined processing ofGPS carrier phase and pseudorange data. For data arcs of two weeks, constrained process noise models forGPS dynamic parameters significantly improved the solutions.     

Monitoring selective availability dither frequencies and their effect on GPS data

Summary The signals transmitted by Block II satellites of the Global Positioning System (GPS) can be degraded to limit the highest accuracy of the system (10 m or better point positioning) to authorized users. This mode of degraded operation is called Selective Availability (S/A). S/A involves the degradation in the quality of broadcast orbits and satellite clock dithering. We monitored the dithered satellite oscillator and investigated the effect of this clock dithering on high accuracy relative positioning. The effect was studied over short 3-meter and zero-baselines with two GPS receivers. The equivalent S/A effects for baselines ranging from 0 to >10,000 km can be examined with short test baselines if the receiver clocks are deliberately mis-synchronized by a known and varying amount. Our results show that the maximum effect of satellite clock dithering on GPS double difference phase residuals grows as a function of the clock synchronization error according to: S/A _effect =0.04 cm/msec, and it increases as a function of baseline length like: S/A _effect =0.014 cm/100 km. These are equations for maximum observed values of post-fit residuals due to S/A. The effect on GPS baselines is likely to be smaller than the 0.14 mm for a baseline separation of 100 km. We therefore conclude, for our limited data set, and for the level of S/A during our tests, that S/A clock dithering has negligible effect on all terrestrial GPS baselines if double difference processing techniques are employed and if the GPS receivers remain synchronized to better than 10 msec. S/A may constitute a problem, however, if accurate point processing is required, or if GPS receivers are not synchronized. We suggest and test two different methods to monitor satellite frequency offsets due to S/A. S/A modulates GPS carrier frequencies in the range of-2 Hz to +2 Hz over time periods of several minutes. The methods used in this paper to measure the satellite clock dither could be applied by the civilian GPS community to continuously monitor S/A clock dithering. The monitored frequencies may aid high accuracy point positioning applications in a postprocessing mode (Malys and Ortiz 1989), and differential GPS with poorly synchronized receivers (Feigl et al. 1991).     

Very high-rate GPS for measuring dynamic seismic displacements without aliasing: performance evaluation of the variometric approach

Very high-rate global positioning system (GPS) data has the capacity to quickly resolve seismically related ground displacements, thereby providing great potential for rapidly determining the magnitude and the nature of an earthquake’s rupture process and for providing timely warnings for earthquakes and tsunamis. The GPS variometric approach can measure ground displacements with comparable precision to relative positioning and precise point positioning (PPP) within a short period of time. The variometric approach is based on single-differencing over time of carrier phase observations using only the broadcast ephemeris and a single receiver to estimate velocities, which are then integrated to derive displacements. We evaluate the performance of the variometric approach to measure displacements using 50 Hz GPS data, which were recorded during the 2013 M_W 6.6 Lushan earthquake and the 2011 M_W 9.0 Tohoku-Oki earthquake. The comparison between 50 and 1 Hz seismic displacements demonstrates that 1 Hz solutions often fail to faithfully manifest the seismic waves containing high-frequency seismic signals due to aliasing, which is common for near-field stations of a moderate-magnitude earthquake. Results indicate that 10–50 Hz sampling GPS sites deployed close to the source or the ruptured fault are needed for measuring dynamic seismic displacements of moderate-magnitude events. Comparisons with post-processed PPP results reveal that the variometric approach can determine seismic displacements with accuracies of 0.3–4.1, 0.5–2.3 and 0.8–6.8 cm in the east, north and up components, respectively. Moreover, the power spectral density analysis demonstrates that high-frequency noises of seismic displacements, derived using the variometric approach, are smaller than those of PPP-derived displacements in these three components.     

A low-complexity algorithm for fast acquisition of weak DSSS signal in high dynamic environment

The initial acquisition of direct-sequence spread-spectrum (DSSS) signal transmitted in bursts by ground terminals at satellite-borne receiver poses an engineering challenge. We propose a low-complexity acquisition algorithm that is capable of capturing extremely weak DSSS signal in the presence of large Doppler dynamics. The algorithm uses fast Fourier transform (FFT)-based frequency-domain techniques to implement circular correlations between the received signal and the local pseudo-random noise (PRN) code, and it coherently accumulates the correlation results across multiple PRN code periods, to achieve a sufficient signal–noise ratio for reliable acquisition. We invoke another FFT procedure to perform the coherent accumulation and the fine compensation for the residual Doppler frequency offset. To highlight the advantage of the proposed algorithm, we make a complexity comparison among the proposed algorithm and two other benchmark strategies, namely the modified double block zero padding (MDBZP) and two-dimensional exhaustive search (2D-ES). It is shown that the proposed algorithm has the lowest complexity, which is particularly desirable for satellite-borne receivers where the computational resource is limited. The acquisition performance of the proposed algorithm is verified by theoretical analysis and Monte Carlo simulations and compared with that of MDBZP and 2D-ES. Moreover, we have carried out extensive tests on a hardware verification system, and we show the claimed tradeoff between performance and cost is indeed attainable with the suggested algorithm. Numerically, it is found the proposed algorithm can achieve a detection rate of 0.9 and a false alarm rate of \(10^{ - 5}\) at C/N ₀ = 29.5 dBHz over a Doppler frequency offset range of \(\left[ { - 7.5\,{\text{kHz}},7.5\,{\text{kHz}}} \right]\) in floating-point simulation, which coincides with the analytical results. The same performance is achieved at C/N ₀ = 31 dBHz in fixed-point simulation and at C/N ₀ = 31.5 dBHz on a hardware system.     

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.     

Stochastic estimation of tropospheric path delays in global positioning system geodetic measurements

Water vapor radiometric (WVR) and surface meteorological (SM) measurements taken during three Global Positioning System (GPS) geodetic experiments are used to calculate process noise levels for random walk and first-order Gauss-Markov temporal models of tropospheric path delays. Entire wet and combined wet and dry zenith delays at each network site then are estimated simultaneously with the geodetic parameters without prior calibration. The path delays and corresponding baseline estimates are compared to those obtained with calibrated data and stochastic residual delays. In this manner, the marginal utility of a priori tropospheric calibration is assessed given the ability to estimate the path delays directly using only theGPS data. Estimation of total zenith path delays with appropriate random walk or Gauss-Markov models yields baseline repeatabilities of a few parts in 10⁸. This level of geodetic precision, and accuracy as suggested by analyses on collocated baselines estimated independently by very long baseline interferometry, is comparable to or better than that obtained after path delay calibration usingWVR and/orSM measurements. Results suggest thatGPS data alone have sufficient strength to resolve centimeter-level zenith path delay fluctuations over periods of a few minutes.     

Cycle slip detection and repair for undifferenced GPS observations under high ionospheric activity

We develop a new approach for cycle slip detection and repair under high ionospheric activity using undifferenced dual-frequency GPS carrier phase observations. A forward and backward moving window averaging (FBMWA) algorithm and a second-order, time-difference phase ionospheric residual (STPIR) algorithm are integrated to jointly detect and repair cycle slips. The FBMWA algorithm is proposed to detect cycle slips from the widelane ambiguity of Melbourne–Wübbena linear combination observable. The FBMWA algorithm has the advantage of reducing the noise level of widelane ambiguities, even if the GPS data are observed under rapid ionospheric variations. Thus, the detection of slips of one cycle becomes possible. The STPIR algorithm can better remove the trend component of ionospheric variations compared to the normally used first-order, time-difference phase ionospheric residual method. The combination of STPIR and FBMWA algorithms can uniquely determine the cycle slips at both GPS L ₁ and L ₂ frequencies. The proposed approach has been tested using data collected under different levels of ionospheric activities with simulated cycle slips. The results indicate that this approach is effective even under active ionospheric conditions.     

GPS scintillation over the European Arctic during the November 2004 storms

Small-scale irregularities in the background electron density of the ionosphere can cause rapid fluctuations in the amplitude and phase of radio signals passing through it. These rapid fluctuations are known as scintillation and can cause a Global Positioning System (GPS) receiver to lose lock on a signal. This could compromise the integrity of a safety of life system based on GPS, operating in auroral regions. In this paper, the relationship between the loss of lock on GPS signals and ionospheric scintillation in auroral regions is explored. The period from 8 to 14 November 2004 is selected for this study, as it includes both geomagnetically quiet and disturbed conditions. Phase and amplitude scintillation are measured by GPS receivers located at three sites in Northern Scandinavia, and correlated with losses of signal lock in receivers at varying distances from the scintillation receivers. Local multi-path effects are screened out by rejection of low-elevation data from the analysis. The results indicate that losses of lock are more closely related to rapid fluctuations in the phase rather than the amplitude of the received signal. This supports the idea, suggested by Humphreys et al. (2005) (performance of GPS carrier tracking loops during ionospheric scintillations. Proceedings Internationsl Ionospheric Effects Symposium 3–5 May 2005), that a wide loop bandwidth may be preferred for receivers operating at auroral latitudes. Evidence from the Imaging Riometer for Ionospheric Studies (IRIS) appears to suggest that, for this particular storm, precipitation of particles in the D/E regions may be the mechanism that drives the rapid phase fluctuations in the signal.

Building Extraction from LIDAR Point Cloud Data Using Marked Point Process

This paper presents a new algorithm for building extraction from LIDAR (Light Detection and Ranging) point cloud data on the basis of a marked point process based building model. In this building model, the positions and geometries of buildings are modeled by a point process and its marks, respectively. The geometric marks for buildings include their length, width, direction, height. By Bayesian paradigm, a posterior distribution for the marked point process conditional on the LIDAR point cloud data is obtained. The Reversible Jump Markov Chain Monte Carlo (RJMCMC) based scheme is designed to simulate the posterior distribution. Finally, Maximum A Posteriori (MAP) strategy is used to obtain the optimal building detection. The proposed algorithm is tested by a set of LIDAR point cloud data. The results show its efficiency in complex residential environments.     

Comparison of GOCE-GPS gravity fields derived by different approaches

Several techniques have been proposed to exploit GNSS-derived kinematic orbit information for the determination of long-wavelength gravity field features. These methods include the (i) celestial mechanics approach, (ii) short-arc approach, (iii) point-wise acceleration approach, (iv) averaged acceleration approach, and (v) energy balance approach. Although there is a general consensus that—except for energy balance—these methods theoretically provide equivalent results, real data gravity field solutions from kinematic orbit analysis have never been evaluated against each other within a consistent data processing environment. This contribution strives to close this gap. Target consistency criteria for our study are the input data sets, period of investigation, spherical harmonic resolution, a priori gravity field information, etc. We compare GOCE gravity field estimates based on the aforementioned approaches as computed at the Graz University of Technology, the University of Bern, the University of Stuttgart/Austrian Academy of Sciences, and by RHEA Systems for the European Space Agency. The involved research groups complied with most of the consistency criterions. Deviations only occur where technical unfeasibility exists. Performance measures include formal errors, differences with respect to a state-of-the-art GRACE gravity field, (cumulative) geoid height differences, and SLR residuals from precise orbit determination of geodetic satellites. We found that for the approaches (i) to (iv), the cumulative geoid height differences at spherical harmonic degree 100 differ by only \({\approx }10~\%\) ; in the absence of the polar data gap, SLR residuals agree by \({\approx }96~\%\) . From our investigations, we conclude that real data analysis results are in agreement with the theoretical considerations concerning the (relative) performance of the different approaches.