Position error bound calculation for GNSS using measurement residuals

In safety-of-life applications of satellite navigation, the protection level (PL) equation translates what is known about the pseudo-range errors into a reliable limit on the positioning error. The current PL equations for satellite-based augmentation systems (SBAS) rely on Gaussian statistics. This approach is very practical: the calculations are simple and the receiver computation load is small. However, when the true distributions are far from Gaussian, such a characterization forces an inflation of the PLs that degrades performance. This happens in particular with errors with heavy tail distributions or for which there is not enough data to evaluate the distribution density up to small quantiles. We present a way of computing the optimal protection level when the pseudo-range errors are characterized by a mixture of Gaussian modes. First, we show that this error characterization adds a new flexibility and helps account for heavy tails without losing the benefit of tight core distributions. Then, we state the positioning problem using a Bayesian approach. Finally, we apply this method to PL calculations for the wide area augmentation system (WAAS) using real data from WAAS receivers. The results are very promising: vertical PLs are reduced by 50% without degrading integrity.