Scroll Top

Precise Point Positioning with Ambiguity Resolution

Precise Point Positioning System Overview

Precise Point Positioning with Ambiguity Resolution (PPP-AR) Working Group

Established: 2018

Chair: Jianghui Geng

Mailing List: IGS PPP-AR Working Group Mailing List

The precise point positioning with ambiguity resolution (PPP-AR) working group (WG) investigates the inter-operability of PPP-AR products generated by various analysis centers. The aim of the WG is to analyze the feasibility and benefits of having the IGS adopt a modernized combination process considering the consistency of the satellite clock and bias products.

Charter

Precise point positioning (PPP) is characterized by its undifferenced data processing with predetermined precise satellite products. After resolving undifferenced ambiguities with the aid of code and phase bias products, PPP has the potential to achieve a millimeter-level positioning precision on a global scale, which is consistent with conventional network analysis techniques. Over the past decade, PPP ambiguity resolution (PPP-AR) has been well studied by the GNSS community, and the IGS established the PPP-AR Working Group (WG) in 2018.

The goal of the PPP-AR WG at the first stage is to encourage IGS analysis centers (ACs) to generate phase bias products and further combine them to study their interoperability. At the moment, CODE (Center for Orbit Determination in Europe), GFZ (German Research Center for Geosciences), CNES/CLS (Centre National d’Etudes Spatiales/Collecte Localisation Satellites) and Wuhan University are routinely providing GPS/Galileo phase bias products in the form of observable-specific signal biases (OSBs). Since such code/phase biases are highly correlated with clocks, their cross-validation and combination entail both signal biases and clocks, and an interoperability test has been performed throughout IGS ACs (Banville et al. 2020). In addition, satellite attitudes are also recommended as quaternion products to cope with eclipsing seasons and thus improve the consistency of satellite clocks/biases among ACs (Loyer et al. 2021). In the IGS Repro3, four ACs or institutions (i.e., CODE, CNES/CLS, NRCan, TUGraz) provided OSBs starting from the year of 2000. Their interoperability has been examined and a combined product has been released (Lin et al. 2023). A multi-GNSS satellite phase clock/bias combination software package has also been transferred to the IGS for routine combination tasks.

Based on the achievements above, the goal of the IGS PPP-AR WG for the next stage is to routinely cross-validate each contributing AC’s phase clock/bias products and encourage their continuing efforts in improving phase bias qualities and developing multi-frequency phase clock/bias products.

 

Goals

Phase 1
  • Survey the technical practices and problems from current phase bias estimates.
  • Encourage more ACs to produce GPS/Galileo/BDS phase clock/bias products.
  • Start a cross-validation task to assess the phase clock/bias qualities and improve the consistency among the contributing products.
  • Work with the Bias WG to ensure continuity of clocks/biases at day boundaries by using the integer properties of the combined products.
  • Initiate a project to expose the combined clock/bias solution to open testing for the public.
  • Develop a phase clock/bias combination software open-sourced across the IGS community.
  • Recommend the GB and ACC to introduce a combined clock/bias solution as an official IGS satellite product.
Phase 2
  • Work in partnership with the Clock WG on time scale issues associated with continuous “phase” clocks.
  • Encourage the RT WG to compute and broadcast real-time phase biases by sharing knowledge on phase bias generation.
  • Survey proper models, algorithms, and requirements for multi-frequency phase bias products, and start a pilot project to generate such products.
  • Assess the versatility and scalability of multi-frequency phase biases for all conceivable observable combinations and frequency choices, and start a pilot combination task.

Members

The following table contains a list of people who are supporting the PPP-AR WG initiative and willing to contribute products and/or expertise. This list is likely to grow as more people are aware of the WG.

Name
Institution
Role
Jianghui Geng Wuhan University Chair
Yoaz Bar-Sever JPL Analysis center contribution (future)
Larry Romans JPL Analysis center contribution (future)
Ant Sibthorpe JPL Analysis center contribution (future)
Michael Coleman NRL Time alignment advisor
Qiyuan Zhang Wuhan University Analysis center contribution (MGEX rapid)
Jing Guo Wuhan University Analysis center contribution (MGEX final)
Guo Chen Wuhan University Orbit combination
Zhe Yan Wuhan University Clock/bias combination and web developer
Denis Laurichesse CNES Analysis center contribution
Sylvain Loyer CLS Analysis center contribution
Benjamin Maennel GFZ Analysis center contribution
Zhiguo Deng GFZ Analysis center contribution
Giulio Tagliaferro BIPM Clock analysis
Francesco Gini ESA Analysis center contribution (future)
Erik Schoenemann ESA Analysis center contribution (future)
Nacho Romero CSC (@ESA/ESOC) Analysis center contribution (future)
Tim Springer PosiTim (@ESA/ESOC) Analysis center contribution (future)
Stefan Schaer CODE Analysis center contribution
Baocheng Zhang CAS Technical analysis
Salim Masoumi GA Analysis center coordinator (ACC)
Patrick Dumitraschkewitz TU Graz Analysis center contribution and reference attitude

Recommended software for PPP-AR product validation

PRIDE PPP-AR (Precise Point Positioning with Ambiguity Resolution): https://github.com/PrideLab/PRIDE-PPPAR/

Formats

SINEX Bias Solution INdependent EXchange format for GNSS biases
ORBEX GNSS attitude quaternions

Publications related to IGS PPP-AR WG activities

Banville S, Geng J, Loyer S, Schaer S, Springer T, Strasser S (2020) On the interoperability of IGS products for precise point positioning with ambiguity resolution. Journal of Geodesy, 94, 1-15. https://doi.org/10.1007/s00190-019-01335-w

Lin J, Geng J, Yan Z, Masoumi S, Zhang Q (2023) Correcting antenna phase center effects to reconcile the code/phase bias products from the third IGS reprocessing campaign. GPS solutions, 27, 70. https://doi.org/10.1007/s10291-023-01405-9

Loyer S, Banville S, Geng J, Strasser S (2021) Exchanging satellite attitude quaternions for improved GNSS data processing consistency. Advances in Space Research, 68, 2441-2452. https://doi.org/10.1016/j.asr.2021.04.049

Strasser S, Banville S, Kvas A, Loyer S, Mayer-Gürr T (2021) Comparison and generalization of GNSS satellite attitude models. EGU General Assembly. https://doi.org/10.5194/egusphere-egu21-7825

IGSMail related to IGS PPP-AR WG activities

IGS Repro3 orbit. clock, bias and attitude solutions available (IGSMAIL-8248)

A new keyword in Bias-SINEX for antenna phase center corrections on geometry-free biases (IGSMAIL-8279)

Special Issue “GNSS Precise Point Positioning: Towards Global Instantaneous cm-Level Accuracy”

A special issue of Remote Sensing (ISSN 2072-4292).
Deadline for manuscript submissions: 31 December 2021

Precise point positioning (PPP) using global navigation satellite systems (GNSS) enables accurate positioning worldwide. Recent advances, including improved error source modeling and the modernization of GNSS constellations, have reduced the time required to achieve cm-level accuracies from hours to seconds, creating new possibilities for several applications.

The objective of this Special Issue is to address the remaining technical challenges associated with global instantaneous high-accuracy GNSS positioning. Timely and reliable positioning can only be achieved through careful attention to detail in all system components, including space segment, signal propagation medium, and receiver design. Research topics of interest include but are not limited to satellite orbit dynamics (solar radiation pressure, attitude), atmospheric augmentation (functional and stochastic models, contribution of low-Earth orbit (LEO) satellites, correction dissemination), and end-user algorithms (stochastic models, time correlation, ambiguity validation). Since integrity plays a critical role in many emerging applications, novel network and end-user techniques for improved error source characterization are also of interest.

I hope that by focusing our efforts on these technical challenges, end-users in a wide range of applications including autonomous vehicles, marine navigation, and many more, can experience fast and accurate positioning globally

For more information: https://www.mdpi.com/journal/remotesensing/special_issues/GNSS_PPP

Remote Sensing special issue PPP

Center for Orbit Determination in Europe (CODE) [COD]

Product Line
Bias SINEX
Phase Biases
Start Date
End Date
Sampling Rate
Link
Rapid Yes GPS L1/L2 2019-001 1 D http://ftp.aiub.unibe.ch/CODE/yyyy_M/
https://cddis.nasa.gov/archive/gnss/products/wwww/
Galileo L1/L5 2019-272 1 D
Final Yes GPS L1/L2 2019-001 1 D http://ftp.aiub.unibe.ch/CODE/yyyy/
https://cddis.nasa.gov/archive/gnss/products/wwww/
Galileo L1/L5 2019-272 1 D
MGEX Yes GPS L1/L2 2018-182 1 D http://ftp.aiub.unibe.ch/CODE_MGEX/CODE/yyyy/
https://cddis.nasa.gov/archive/gnss/products/wwww/
Galileo L1/L5
Repro3 Yes GPS L1/L2 2000-124 2022-330 1 D http://ftp.aiub.unibe.ch/REPRO_2020/CODE/yyyy/
https://cddis.nasa.gov/archive/gnss/products/wwww/repro3/
Galileo L1/L5 2014-001 2022-330 1 D

Natural Resources Canada (NRCan) [EMR]

Product Line
Bias SINEX
Phase Biases
Start Date
End Date
Sampling Rate
Link
Repro3 Yes GPS L1/L2 2000-001 2020-366 1 D https://cddis.nasa.gov/archive/gnss/products/wwww/repro3/

German Research Center for Geosciences (GFZ) [GFZ]

Product Line
Bias SINEX
Phase Biases
Start Date
End Date
Sampling Rate
Link
MGEX No GPS WL L1/L2 2021-167 2021-189 1 D ftp://ftp.gfz-potsdam.de/pub/GNSS/products/mgex/wwww/
Galileo WL L1/L5
BDS-2 WL L2/L6
BDS-3 WL L2/L6
Yes GPS L1/L2 2021-190 1 D https://cddis.nasa.gov/archive/gnss/products/wwww/
Galileo L1/L5
BDS-2 L2/L6
BDS-3 L2/L6

Centre National d’Etudes Spatiales (CNES) / Collecte Localisation Satellites (CLS) [GRG]

Product Line
Bias SINEX
Phase Biases
Start Date
End Date
Sampling Rate
Link
Rapid No GPS WL L1/L2 2023-003 2023-330 1 D ftp://ftpsedr.cls.fr/pub/igsac/Wide_lane_GPS_satellite_biais.wsb
ftp://ftpsedr.cls.fr/pub/igsac/Wide_lane_GAL_satellite_biais.wsb
https://cddis.nasa.gov/archive/gnss/products/wwww/
Galileo WL L1/L5
Yes GPS L1/L2 2023-331 1 D https://cddis.nasa.gov/archive/gnss/products/wwww/
Galileo L1/L5
BDS-3 L2/L6
Final Yes GPS L1/L2 2022-331 1 D https://cddis.nasa.gov/archive/gnss/products/wwww/
Galileo L1/L5
MGEX No GPS WL L1/L2 2009-305 2021-135 1 D ftp://ftpsedr.cls.fr/pub/igsac/Wide_lane_GPS_satellite_biais.wsb
ftp://ftpsedr.cls.fr/pub/igsac/Wide_lane_GAL_satellite_biais.wsb
https://cddis.nasa.gov/archive/gnss/products/wwww/
Galileo WL L1/L5 2018-214 2021-135 1 D
Yes GPS L1/L2 2021-136 1 D https://cddis.nasa.gov/archive/gnss/products/wwww/
Galileo L1/L5 2021-136 1 D
BDS-2 L2/L6 2023-316 1 D
BDS-3 L2/L6 2023-316 1 D
Repro3 No GPS WL L1/L2 2000-124 2020-366 1 D ftp://ftpsedr.cls.fr/pub/igsac/Wide_lane_GPS_satellite_biais.wsb
ftp://ftpsedr.cls.fr/pub/igsac/Wide_lane_GAL_satellite_biais.wsb
https://cddis.nasa.gov/archive/gnss/products/wwww/repro3/
Galileo WL L1/L5 2017-001 2020-366 1 D

Graz University of Technology [TUG]

Product Line
Bias SINEX
Phase Biases
Start Date
End Date
Sampling Rate
Link
Repro3 Yes GPS L1/L2 1994-001 2020-366 1 D https://cddis.nasa.gov/archive/gnss/products/wwww/repro3/
GPS L5 2010-240 2020-366 1 D
Galileo L1/L5/L7/L8 2013-001 2020-366 1 D

Centre National d’Etudes Spatiales (CNES) / Navigation Team

Product Line
Bias SINEX
Phase Biases
Start Date
End Date
Sampling Rate
Link
Rapid Yes GPS L1/L2/L5 2019-001 30 S http://www.ppp-wizard.net/products/POST_PROCESSED/
Galileo L1/L5/L6/L7
BDS-2 L2/L6/L7
BDS-3 L1/L2/L5/L6/L7
Real time Yes GPS L1/L2/L5 2019-001 30 S http://www.ppp-wizard.net/products/REAL_TIME/
Galileo L1/L5/L6/L7
BDS-2 L2/L6/L7
BDS-3 L2/L6/L7

Wuhan University [WUM]

Product Line
Bias SINEX
Phase Biases
Start Date
End Date
Sampling Rate
Link
MGEX (Rapid) Yes GPS L1/L2 2020-001 2022-365 1 D
ftp://igs.gnsswhu.cn/pub/whu/phasebias/
Galileo L1/L5
BDS-2 L2/L6
BDS-3 L2/L6
Yes GPS L1/L2/L5 2023-001 1 D
Galileo L1/L5/L6/L7/L8
BDS-2 L2/L6/L7
BDS-3 L1/L2/L5/L6/L7
MGEX (Final) Yes GPS L1/L2 2023-163 1 D https://cddis.nasa.gov/archive/gnss/products/wwww/
Galileo L1/L5
BDS-2 L2/L6
BDS-3 L2/L6

2021-09-16

Meeting to discuss the impact of phase center offsets on bias computation.

Summary: PPPAR-WG_20210916_Summary.pdf
Presentation: PPPAR-WG_20210916_Presentation1.pdf

2022-06-30—2022-08-16

Email communications to discuss a new keyword in Bias-SINEX denoting APC corrections.

IGS Mail announcement: IGSMAIL-8279

2023-11-14

Virtual meeting to discuss a new charter and the webpages showing the combination results, etc.

Summary: PPPAR-WG_20231114_Summary.pdf
Presentation: PPPAR-WG_20231114_Presentation.pdf

Combination and analysis of rapid products

The rapid clock and phase bias products provided by different ACs are combined and evaluated.

These products are accessed from the CDDIS archive (https://cddis.nasa.gov/archive/gnss/products/) and the corresponding ACs are listed in the following table. Currently this combination and analysis are only applied to the GPS, Galileo and BDS since the phase bias products provided by ACs are available for these constellations. And other systems are not included in the analysis at this stage.

ID
Rapid
name system orbit clock attitude code bias phase bias APC model
COD COD0OPSRAP GE 300 s 30 s 30 s GE (OSB) GE (OSB) applied
EMR EMR0OPSRAP G 900 s 300 s
ESA ESA0OPSRAP G 900 s 300 s
GFZ GFZ0MGXRAP GEC 300 s 30 s 30 s GEC (OSB) GEC (OSB) applied
GRG GRG0OPSRAP GEC 300 s 30 s 30 s GEC (OSB) GEC (OSB) applied
IGS IGS0OPSRAP G 900 s 300 s
JGX JGX0OPSRAP GE 300 s 30 s
JPL JPL0OPSRAP G 300 s 30 s 30 s
USN USN0OPSRAP G 900 s 300 s
WUM WUM0MGXRAP GEC 300 s 30 s 30 s GEC (OSB) GEC (OSB) applied
Note:
a. In the combination and analysis, rapid products from AC are represented using their ID + “r”, e.g. “CODr” is for the rapid products from COD.
b. GEC stand for GPS, Galileo and BDS-2/3 respectively.
c. OSB stands for observable specific bias, IRC stands for integer recovery clock.

Evaluation of combined results

For the methods used in the clock and bias combination, please refer to the tabs Resources.

All the products in the table are included in the combination except for the IGS products, which is used for comparison. The COD’s rapid orbit products are used as reference orbit, and the reference attitude is generated by the open source software GROOPS. The GNSS rapid clock products are combined with a sampling interval of 300 s. All the results are presented on a weekly basis, with a latency of one week.

The combined weight

The clock combination employs an iterative weighting method in which the weight assigned to a specific product depends on its residuals in relation to the combined clock.

The weekly clock/bias RMSE

The clock/bias RMSE refers to the daily RMSE of clock/bias for each AC with respect to the combined integer clock, which reflects the precision of clock and bias combination and the consistency between individual ACs. Each grid represents a satellite on a particular day. Blank grids mean unavailable products and a slash inside means an outlier excluded from the combination. The line chart below shows the satellite clock outlier rate per day for each AC, with gray block indicating that relevant satellite clocks do not participate in the comparison. The statistics at the bottom indicate the overall RMSE of AC clock/bias for this week.

The accumulated clock/bias RMSE

The accumulated clock/bias RMSE refers to the daily RMSE of clock/bias for all the ACs with respect to the combined integer clock, which reflects the change of the combined accuracy over time. The line chart below shows the overall RMSE for each constellation.

PPP-AR validation

The GPS/Galileo/BDS data with a sampling interval of 300 s from 10 globally distributed stations are processed for PPP-AR in a static mode with the PRIDE PPP-AR software. The fixing rate and position precision of each single constellation solution are presented in the figure below. The JAXA daily SINEX products are used as reference solution.  “CMBr” stands for the combined rapid products.

Combination and analysis of final products

The final clock and phase bias products provided by different ACs are combined and evaluated.

These products are accessed from the CDDIS archive (https://cddis.nasa.gov/archive/gnss/products/) and the corresponding ACs are listed in the following table. Currently this combination and analysis are only applied to the GPS, Galileo and BDS since the phase bias products provided by ACs are available for these constellations. And other systems are not included in the analysis at this stage.

ID
Final
name system orbit clock attitude code bias phase bias APC model
COD COD0OPSFIN GE 300 s 30 s 30 s GE (OSB) GE (OSB) applied
ESA ESA0OPSFIN GE 300 s 30 s
GFZ GFZ0OPSFIN GE 300 s 30 s
GRG GRG0MGXFIN GEC 300 s 30 s 30 s GEC (OSB) GEC (OSB) applied
IGS IGS0OPSFIN G 900 s 30 s
JAX JAX0MGXFIN GE 300 s 30 s
JGX JGX0OPSFIN GE 300 s 30 s
JPL JPL0OPSFIN G 300 s 30 s 30 s
WUM WUM0MGXFIN GEC 300 s 30 s 30 s GEC (OSB) GEC (OSB) applied
Note:
a. In the combination and analysis, final products from AC are represented using their ID + “f”, e.g., “CODf” is for the final products from COD.
b. GEC stand for GPS, Galileo and BDS-2/3 respectively.
c. OSB stands for observable specific bias, IRC stands for integer recovery clock.

Evaluation of combined results

For the methods used in the clock and bias combination, please refer to the tabs Resources.

All the products in the table are included in the combination except for the IGS products, which is used for comparison. The reference orbit products are generated from the Wuhan University and the reference attitude is generated by the open source software GROOPS. The GNSS final clock products are combined with a sampling interval of 30 s. All the results are presented on a weekly basis, with a latency of three weeks.

The combined weight

The clock combination employs an iterative weighting method in which the weight assigned to a specific product depends on its residuals in relation to the combined clock.

The weekly clock/bias RMSE

The clock/bias RMSE refers to the daily RMSE of clock/bias for each AC with respect to the combined integer clock, which reflects the precision of clock and bias combination and the consistency between individual ACs. Each grid represents a satellite on a particular day. Blank grids mean unavailable products and a slash inside means an outlier excluded from the combination. The line chart below shows the satellite clock outlier rate per day for each AC, with gray block indicating that relevant satellite clocks do not participate in the comparison. The statistics at the bottom indicate the overall RMSE of AC clock/bias for this week.

The accumulated clock/bias RMSE

The accumulated clock/bias RMSE refers to the daily RMSE of clock/bias series for all the ACs with respect to the combined integer clock, which reflects the change of the combined accuracy over time. The line chart below shows the overall RMSE for each constellation.

PPP-AR validation

The GPS/Galileo/BDS data with a sampling interval of 300 s from 10 globally distributed stations are processed for PPP-AR in a static mode with the PRIDE PPP-AR software. The fixing rate and position precision of each single constellation solution are presented in the figure below. IGS daily SINEX products are used as reference solution. “CMBf” stands for the combined final products.

Last Updated on 8 Apr 2024 09:20 UTC

Privacy Preferences
When you visit our website, it may store information through your browser from specific services, usually in form of cookies. Here you can change your privacy preferences. Please note that blocking some types of cookies may impact your experience on our website and the services we offer.