The Bayesian detection of discontinuities in a polynomial regression and its application to the cycle-slip problem

This paper deals with the problem of detecting and correcting cycle-slips in Global Navigation Satellite System (GNSS) phase data by exploiting the Bayesian theory. The method is here applied to undifferenced observations, because repairing cycle-slips already at this stage could be a useful pre-processing tool, especially for a network of permanent GNSS stations. If a dual frequency receiver is available, the cycle-slips can be easily detected by combining two phase observations or phase and range observations from a single satellite to a single receiver. These combinations, expressed in a distance unit form, are completely free from the geometry and depend only on the ionospheric effect, on the electronic biases and on the initial integer ambiguities; since these terms are expected to be smooth in time, at least in a short period, a cycle-slip in one or both the two carriers can be modelled as a discontinuity in a polynomial regression. The proposed method consists in applying the Bayesian theory to compute the marginal posterior distribution of the discontinuity epoch and to detect it as a maximum a posteriori (MAP) in a very accurate way. Concerning the cycle-slip correction, a couple of simultaneous integer slips in the two carriers is chosen by maximazing the conditional posterior distribution of the discontinuity amplitude given the detected epoch. Numerical experiments on simulated and real data show that the discontinuities with an amplitude 2 or 3 times larger than the noise standard deviation are successfully identified. This means that the Bayesian approach is able to detect and correct cycle-slips using undifferenced GNSS observations even if the slip occurs by one cycle. A comparison with the scientific software BERNESE 5.0 confirms the good performance of the proposed method, especially when data sampled at high frequency (e.g. every 1 s or every 5 s) are available.     

Geodetic intersection on the ellipsoid

Through each of two known points on the ellipsoid a geodesic is passing in a known azimuth. We solve the problem of intersection of the two geodesics. The solution for the latitude is obtained as a closed formula for the sphere plus a small correction, of the order of the eccentricity of the ellipsoid, which is determined by numerical integration. The solution is iterative. Once the latitude is obtained, the longitude is determined without iteration.     

Prediction of UT1–UTC,LOD and AAM χ<Subscript>3</Subscript> by combination of least-squares and multivariate stochastic methods

This article presents the application of a multivariate prediction technique for predicting universal time (UT1–UTC), length of day (LOD) and the axial component of atmospheric angular momentum (AAM χ ₃). The multivariate predictions of LOD and UT1–UTC are generated by means of the combination of (1) least-squares (LS) extrapolation of models for annual, semiannual, 18.6-year, 9.3-year oscillations and for the linear trend, and (2) multivariate autoregressive (MAR) stochastic prediction of LS residuals (LS + MAR). The MAR technique enables the use of the AAM χ ₃ time-series as the explanatory variable for the computation of LOD or UT1–UTC predictions. In order to evaluate the performance of this approach, two other prediction schemes are also applied: (1) LS extrapolation, (2) combination of LS extrapolation and univariate autoregressive (AR) prediction of LS residuals (LS + AR). The multivariate predictions of AAM χ ₃ data, however, are computed as a combination of the extrapolation of the LS model for annual and semiannual oscillations and the LS + MAR. The AAM χ ₃ predictions are also compared with LS extrapolation and LS + AR prediction. It is shown that the predictions of LOD and UT1–UTC based on LS + MAR taking into account the axial component of AAM are more accurate than the predictions of LOD and UT1–UTC based on LS extrapolation or on LS + AR. In particular, the UT1–UTC predictions based on LS + MAR during El Niño/La Niña events exhibit considerably smaller prediction errors than those calculated by means of LS or LS + AR. The AAM χ ₃ time-series is predicted using LS + MAR with higher accuracy than applying LS extrapolation itself in the case of medium-term predictions (up to 100 days in the future). However, the predictions of AAM χ ₃ reveal the best accuracy for LS + AR.     

Bringing GPS into harsh environments for fully automated deformation monitoring

Engineering projects that require deformation monitoring frequently utilize geodetic sensors to measure displacements of target points located in the deformation zone. In situations where control stations and targets are separated by a kilometer or more, GPS can offer higher precision position updates at more frequent intervals than can normally be achieved using total station technology. For large-scale deformation projects requiring the highest precision, it is therefore advisable to use a combination of the two sensors. In response to the need for high precision, continuous GPS position updates in harsh deformation monitoring environments, a software has been developed that employs triple-differenced carrier-phase measurements in a delayed-state Kalman filter. Two data sets were analyzed to test the capabilities of the software. In the first test, a GPS antenna was displaced using a translation stage to mimic slow deformation. In the second test, data collected at a large open pit mine were processed. It was shown that the delayed-state Kalman filter developed could detect millimeter-level displacements of a GPS antenna. The actual precision attained depends upon the amount of process noise infused at each epoch to accommodate the antenna displacements. Higher process noise values result in quicker detection times, but at the same time increase the noise in the solutions. A slow, 25 mm displacement was detected within 30 min of the full displacement with sigma values in E, N and U of ±10 mm or better. The same displacement could also be detected in less than 5 h with sigma values in E, N and U of ±5 mm or better. The software works best for detecting long period deformations (e.g., 20 mm per day or less) for which sigma values of 1–2 mm are attained in all three solution components. It was also shown that the triple-differenced carrier-phase observation can be used to significantly reduce the effects of residual tropospheric delay that would normally plague double-differenced observations in harsh GPS environments.

Spin rate estimation of sounding rockets using GPS wind-up

Carrier phase wind-up is a well-known effect that arises from the relative rotation between a transmitting and receiving antenna. In GPS measurements at L1 frequency, this effect translates into an error of 19.029 cm per full relative rotation of antennas. Since this effect is independent of the satellite elevation for pure rotation about the antenna boresight axis, it is usually absorbed by the clock estimation in navigation algorithms. Therefore, the impact of wind-up is usually neglected for applications that do not require accuracies to the cm level like RTK. However, in receiving platforms with high rotation rate, the accumulated wind-up value can be important and actually be larger than receiver noise or even ionospheric variations. Therefore, in such scenarios, the wind-up contribution can be isolated and used as a source of information to compute the spin rate of such platforms using an appropriate combination of GPS observables. This work shows some results of a coarse, yet simple, approach to monitor the rotation angle and spin-rate of spin stabilized sounding rockets flown by DLR.     

An approach to aid INS using time-differenced GPS carrier phase (TDCP) measurements

This paper presents a simple and effective approach that incorporates single-frequency, L1 time-differenced GPS carrier phase (TDCP) measurements without the need of ambiguity resolution techniques and the complexity to accommodate the delayed-state terms. Static trial results are included to illustrate the stochastic characteristics and effectiveness of the TDCP measurements in controlling position error growth. The formulation of the TDCP observation model is also described in a 17-state tightly-coupled GPS/INS iterative, extended Kalman filter (IEKF) approach. Preliminary land vehicle trial results are also presented to illustrate the effectiveness of the TDCP which provides sub-meter positional accuracies when operating for more than 10 min.     

GPS signal acquisition using the repeatability of successive code phase measurements

This paper introduces a new method for GPS signal acquisition, which is based on the repeatability of successive code phase measurements and the M-of-N search algorithm. The performance of the proposed method in terms of probability of signal detection is similar to that of traditional methods, except that the calculation of the probability of detection does not rely on the noise distribution or the Carrier-to-Noise ratio (C/N₀). The code phase repeatability-based method is presented along with equations for probability of detection and probability of false detection. If the distribution of the noise is known, it also provides an estimate of the C/N₀. The proposed method is illustrated for coherent and non-coherent acquisition and C/N₀ estimation.     

On the Use of CloudSat and MODIS Data for Estimating Hurricane Intensity

This letter presents preliminary results concerning the use of new observations from the A-Train Constellation for testing a new technique of remotely sensing hurricane intensity from space based on modeling a hurricane as a balanced, convectively neutral vortex. The key observational requirements are simultaneous, accurate measurements of cloud-top height, cloudtop temperature, and cloud profiling information across the center of the storm, although there are ways to bypass the need for cloud-top temperature. In this letter, the Moderate Resolution Imaging Spectroradiometer onboard Aqua provides an estimation of the cloud-top temperature, and the near-simultaneous CloudSat observations provide the essential cloud-top height and cloud profiling information. Initial results indicate that the new technique is a promising method for estimating storm intensity when compared post facto to the best track database. Potential uncertainties and room for further refinement of the technique are discussed.