Atomium is a software developped by the ROB for geodetic positionning and time transfer. Based on GNSS code and carrier phase measurements analysis, it was originally dedicated to perform Precise Point Positioning (PPP), but it can also perform GNSS based time and frequency transfer in Single Difference mode. The PPP procedure is fully described in (Kouba-Heroux 2001). It uses a least-square adjustment to estimate several parameters :

  • Station position
  • Station clock synchronistation error (with the respect to the satellite clock correction time scale)
  • Tropospheric zenith wet delay
  • Phase ambiguities (float)

Atomium is based on the ionosphere-free combinations of L1 and L2 carrier phase measurements and P1 and P2 pseudorange measurements, named L3 and P3 respectively. The observation can be used at 5-minute or 30-second sampling rate.

Main Features

  • Performs PPP and single difference Time Transfer at the 0.1 ns precision level
  • Compued station position in PPP at the cm level
  • Allows to combine GPS and GLONASS observables in PPP and in single differences
  • Possibility to correct the measurements for second and third order ionospheric delays
  • Possibility to change the sampling rate of the tropospheric zenith wet delay estimation and to constrain its evolution
  • Several troposphere mapping function implemented (NMF and GMF so far)
  • Adjustable threshold for the data screening
  • Pseudorange smoothing


  • Earth Tide
    • Station position is corrected for its time variation due to degrees 2 and 3 solid Earth tides as given in the IERS convention (McCarthy-Petit 2003). The rotational deformation correction on station position due to Earth pole mean motion are also taken into account in the station motion.
  • Ocean loading
    • The FES 2004 model (Lettier et al 2004) is implemented in order to correct ocean loading displacements.
  • Troposphere
    • The hydrostatic part of the tropospheric delay is computed with the Saastamoinan model of dry atmosphere (Saastamoinen 1972).
    • The tropospheric mapping function implemented in Atomium are the Niell mapping function (Niell 1996) and Global Mapping Function (GMF, Boehm et al 2006).
    • The wet part of the zenith path delay is estimated in the least-square adjustment.
  • Ionosphere
    • The ionosphere free combination is used to eliminate first order perturbations. Higher order ionosphere effects (second and third order delays on code and phase) are implemented as described in the latest IERS 2009 conventions.
  • Satellite orbits and clocks
    • The satellite positions are obtained with a Neville interpolation on 12 points of the IGS orbit files.
    • The satellite clock correction are taken from the IGS CLK files in which the sampling rate is 5 minutes.
    • For receiver measuring only C1 and P2 pseudorange measurements, the satellite biases C1-P1 provided by CODE are removed from the C1 pseudorange measurements.
  • Antenna Phase center variations
    • Both satellite and receiver antennas are corrected for respectively nadir-dependent and elevation-dependent phase center variations. Since GPS week 1400, absolute corrections provided by the IGS (file igs05.atx) are applied (Gendt 2006).
    • The carrier-phase measurements are corrected for phase windup (Wu et al 1993) taking into account the satellite attitude and eclipse events.