IGS SSR Format Overview v1.00
IGS SSR Format Overview v1.00
00 2020-10-05
Acknowledgement: RTCM SC-104 SSR Working Group, IGS Real-Time Working Group
IGS SSR v1.00 2020-10-05
Table of Contents
0. REVISION HISTORY................................................................................................................ 1
1. THE PHILOSOPHY AND HISTORY OF IGS SSR ................................................................. 2
2. GENERAL FORMAT DESCRIPTION ..................................................................................... 3
3. BASIC DEFINITIONS ............................................................................................................... 3
4. SSR REPRESENTATION MODELS ........................................................................................ 4
4.1 SSR Satellite Orbit Correction .............................................................................................. 6
4.2 SSR Satellite Clock Correction ........................................................................................... 10
4.3 SSR Satellite Bias Corrections ............................................................................................ 11
4.4 SSR URA ............................................................................................................................ 15
4.5 SSR Ionosphere ................................................................................................................... 15
4.5.1 Ionosphere Vertical TEC .................................................................................................. 15
5. DATA TYPE SUMMARY ....................................................................................................... 18
6. SUB-TYPE MESSAGE SUMMARY ...................................................................................... 18
7. DATA FIELDS ......................................................................................................................... 22
7.1 Data Field Summary............................................................................................................ 22
7.2 Additional GNSS Specific Notes ........................................................................................ 33
7.2.1 GPS ............................................................................................................................... 33
7.2.2 GLONASS .................................................................................................................... 33
7.2.3 Galileo .......................................................................................................................... 33
7.2.4 QZSS ............................................................................................................................ 34
7.2.5 BDS .............................................................................................................................. 34
7.2.6 SBAS ............................................................................................................................ 34
8. SSR CORRECTION MESSAGES ........................................................................................... 34
8.1 GNSS SSR Orbit and Clock Corrections ............................................................................ 34
8.2 GNSS SSR Bias Corrections ............................................................................................... 38
8.3 GNSS SSR URA Messages ................................................................................................ 41
8.4 SSR IONOSPHERE CORRECTION MESSAGE .............................................................. 42
8.4.1 SSR Ionosphere VTEC Spherical Harmonics Message ............................................... 42
9 REFERENCES .......................................................................................................................... 45
ii
IGS SSR v1.00 2020-10-05
ACRONYMS
iii
IGS SSR v1.00 2020-10-05
iv
IGS SSR v1.00 2020-10-05
0. REVISION HISTORY
Version 1.00
05 Oct 2020 Release Version
1
IGS SSR v1.00 2020-10-05
The RTCM-SSR WG defined the development of SSR concepts and real-time capable messages
for all types of accuracies, including ambiguity resolution (RTK quality). The proposed schedule
or work plan consisted of the following major steps:
• stage 1 with satellite orbit, satellite clock and satellite code bias messages
to enable code-based Real-Time (RT)-PPP for dual-frequency receivers
• stage 2 with vertical Total Electron Content (VTEC) ionospheric message
to enable code-based RT-PPP for single frequency receivers and satellite phase bias
messages to enable phase-based RT-PPP
• stage 3 with ionospheric slant TEC (STEC) and tropospheric messages
to enable PPP-Real-time-Kinematic (RTK)
• stage 4 with compression to reduce bandwidth
The concept of SSR requires the transmission of a bias (code or phase bias) for every supported
signal. Signals with no bias information are unsupported and should not be used by receiving
applications.
Standardized RTCM-SSR Messages for GPS and GLONASS were first published in ”RTCM
STANDARD 10403.1 with Amendments 1-5, July 1, 2011”. Only RTCM stage 1 messages have
so far been standardized and out of those, only the ones that cover GPS and GLONASS.
The format was developed to be general and flexible to support different GNSS SSR estimation
and modeling approaches. Accordingly, the initial focus was put on the content of the messages,
shifting bandwidth optimization (compression) to the last development stage. In the meantime, the
general validity of the concept and the mapping of known modeling approaches have been
demonstrated in theory and practice.
The IGS has developed and published this standard with stage 1 and 2 messages for multi-GNSS
SSR correction dissemination. This is an open format that can be used by the wider community,
but in particular to support the IGS real-time service as well as broader research and scientific
application. Although a different message format is used, the basic contents of the IGS-SSR
representation are compatible with RTCM-SSR contents.
2
IGS SSR v1.00 2020-10-05
The IGS SSR format is an open standard for dissemination of real-time products to support the
IGS Real-Time Service and the wider community. The messages support multi-GNSS and include
corrections for orbits, clocks, code-biases resp. Differential Code Biases (DCBs), phase-biases and
ionospheric delays. Extensions to also cover satellite attitude, phase center offsets and variations
and group delay variations are planned in the near future. The goal is to create a self-contained and
scalable standard for a wide range of real-time applications.
RTCM proprietary message types are assigned to specific companies and organizations in the
RTCM 10403.3 Standard for the broadcast of proprietary information. The format of the
proprietary messages corresponds to the RTCM messages, in that the transport layer is defined in
the same way, and the first data field is a 12-bit message number for an organization. The
organization is free to define several sub-types of messages, but they must all utilize the assigned
message type.
The IGS RTCM Proprietary SSR messages uses the IGS assigned RTCM proprietary message
type:
Within the IGS Proprietary SSR Messages several sub-types of messages are defined, and they all
utilize the IGS assigned RTCM3 message type “4076”.
Within this document the term “shall” expresses a required specification. The term “should”
is dealt with as a recommendation.
3. BASIC DEFINITIONS
The principle of the state space concept is to provide information on individual error sources acting
on GNSS to enable an improved positioning for a user. Therefore, the term “State Space
Representation” is used.
3
IGS SSR v1.00 2020-10-05
The GNSS state vector of SSR consists of the following basic parameters:
Satellite orbit and clock errors are corrections to the broadcast information provided by a GNSS.
Satellite signal biases are delays to codes and carrier phases observations within satellite hard- and
software. The GNSS signals are affected by atmospheric propagation delays due to the ionosphere
and the troposphere.
The software system used to generate SSR corrections must apply different corrections due to
effects like reference station site displacements, relativity, phase wind-up, satellite and receiver
antenna phase center variations (PCVs) and group delay variations (GDVs), atmospheric delays
and others to the used GNSS observations. The estimated SSR parameters are then free of such
effects. Depending on the quality of the service, a client application (rover) must apply
corresponding corrections when using SSR parameters in order to determine its position in a
conventional coordinate system.
Reference for correction and correction models are the latest IERS Conventions and Wu et al.
(1991) for satellite phase wind-up.
4
IGS SSR v1.00 2020-10-05
parameters from different Update Intervals in order to obtain consistent sets of corrections. The
meaning of SSR Update Interval is different from the transmission interval/rate of SSR messages.
The validity interval of SSR parameters is at least the SSR Update Interval.
A state space parameter may consist of different constituents disseminated in different SSR
correction messages. The use of different SSR messages is intentional to support different
applications, update rates and accuracy requirements. Additional SSR correction messages will
consequently add additional resolution and positioning accuracy. This creates the need to know
the consistency of the SSR parameters for an application. Only a consistent set of SSR parameters
can build up a complete and accurate correction. The consistency of SSR data becomes more
important with increasing resolution provided by additional SSR correction messages.
Generally, the continuous chronology of messages can be used to check the consistency, but in
real-time applications messages may be lost or delayed. A consistency parameter is also the GNSS
specific Issue of Data (IOD), which is used to obtain consistency in the computation of SSR orbit
and clock corrections. A similar requirement exists for state space parameters distributed over
several SSR messages.
There are several SSR Update Intervals defined, which all start at the start of a day in the SSR time
scale (time: 00:00:00). The SSR message contains the SSR Epoch Time 1 s data field (IDF003).
The consistency of the SSR parameters can be verified from this information. The rover is then
capable to combine all relevant SSR parameters consistently and it secures the combination of
state parameters from different messages.
The SSR Update Interval serves to uniquely identify all consistent parameters from different
messages, which can be used to compute consistent corrections at one epoch.
The reference time for the correction terms is computed from the SSR Epoch Time 1 s data field
(IDF003) plus half the SSR Update Interval. Exception is SSR Update Interval “0”, which uses
the SSR Epoch Time as reference time.
GNSS IOD information of broadcast clocks is redundant and is not required in addition to the
broadcast orbit IOD. I.e., the SSR Orbit and Clock corrections refer to the same IOD. The service
provider should refer to the latest set of broadcast messages, which are generally also received in
real-time by a GNSS rover. However, it is recommended to delay the use of the latest broadcast
message for a period of 60 seconds, measured from the time of complete reception of ephemeris
and clock parameters, in order to accommodate rover applications to obtain the same set of
broadcast orbital and clock parameters. It is left to the service provider to send the RTCM Satellite
Ephemeris Data (e.g. Message Types (MT)1019 for GPS and MT1020 for GLONASS) before
switching to a new IOD. It is not allowed to send corrections for more than one IOD for the same
orbit Update Interval. The consistency of SSR Orbit and Clock messages is maintained through
the SSR Update Interval.
5
IGS SSR v1.00 2020-10-05
The orbit and clock messages contain data to be combined with the corresponding values obtained
from the satellites broadcast message. See Table 1 for the type of broadcast navigation messages
used per individual GNSS.
𝑿𝑜𝑟𝑏𝑖𝑡 = 𝑿𝑏𝑟𝑜𝑎𝑑𝑐𝑎𝑠𝑡 − 𝛿𝑿
6
IGS SSR v1.00 2020-10-05
with
𝑿𝑜𝑟𝑏𝑖𝑡 satellite position corrected by SSR Orbit Correction message
𝑿𝑏𝑟𝑜𝑎𝑑𝑐𝑎𝑠𝑡 satellite position computed according to corresponding
GNSS Interface Control Document (ICD) from broadcast ephemeris
parameter set identified by the GNSS IOD in SSR Orbit Correction
message
𝛿𝑿 satellite position correction
𝒓̇
𝒆𝑎𝑙𝑜𝑛𝑔 = |𝒓̇ |
𝒓 × 𝒓̇
𝒆𝑐𝑟𝑜𝑠𝑠 = |𝒓
× 𝒓̇ |
𝒆𝑟𝑎𝑑𝑖𝑎𝑙 = 𝒆𝑎𝑙𝑜𝑛𝑔 × 𝒆𝑐𝑟𝑜𝑠𝑠
𝛿𝑿 = [𝒆𝑟𝑎𝑑𝑖𝑎𝑙 𝒆𝑎𝑙𝑜𝑛𝑔 𝒆𝑐𝑟𝑜𝑠𝑠 ] 𝛿𝑶
with
𝒓 = 𝑿𝒃𝒓𝒐𝒂𝒅𝒄𝒂𝒔𝒕 Earth Centered Earth Fixed (ECEF) satellite broadcast position
vector
𝒓̇ = 𝑿̇𝒃𝒓𝒐𝒂𝒅𝒄𝒂𝒔𝒕 ECEF satellite broadcast velocity vector
𝒆𝑖 direction unit vector, i = {radial, along, cross}
𝛿𝑶 orbit correction vector
Note: The radial vector 𝒆𝑟𝑎𝑑𝑖𝑎𝑙 is used in the computation and should not be confused with a radial
vector of a circular orbit.
The complete orbit correction vector δO is computed from the individual correction terms and their
velocities:
𝛿𝑂𝑟𝑎𝑑𝑖𝑎𝑙 𝛿𝑂̇𝑟𝑎𝑑𝑖𝑎𝑙
𝛿𝑶 = [ 𝛿𝑂𝑎𝑙𝑜𝑛𝑔 ] + [ 𝛿𝑂̇𝑎𝑙𝑜𝑛𝑔 ] (𝑡 − 𝑡0 )
𝛿𝑂𝑐𝑟𝑜𝑠𝑠 𝛿𝑂̇𝑐𝑟𝑜𝑠𝑠
with
𝑡 time
𝑡0 reference time obtained from SSR Orbit Correction message
𝛿𝑂𝑖 , 𝛿𝑂̇𝑖 orbit correction terms from SSR Orbit message, i = {radial, along, cross}
The ECEF satellite broadcast velocity vector for Geostationary Earth Orbit satellites (GEO) can
be null, which gives no SSR orbit corrections using the above set of equations. The velocity vector
for GEO satellites is therefore defined in a non-rotating system parallel to the International
Terrestrial Reference Frame (ITRF) at SSR epoch time. The GPS value of the angular velocity of
the Earth around the Z-axis is used.
7
IGS SSR v1.00 2020-10-05
−Ω̇𝑒 ∙ 𝑌𝑏𝑟𝑜𝑎𝑑𝑐𝑎𝑠𝑡
𝑟̇ = 𝑟̇𝑏𝑟𝑜𝑎𝑑𝑐𝑎𝑠𝑡 + [ Ω̇𝑒 ∙ 𝑋𝑏𝑟𝑜𝑎𝑑𝑐𝑎𝑠𝑡 ]
0
with
Orbit representation requires the definition of a Coordinate Reference System (CRS). For global
services, the CRS should be related to the International Terrestrial Reference System (ITRS). For
regional services, a CRS related to the tectonic plate of the region can be used. Such continental
or regional CRS (e.g. European Terrestrial Reference System (ETRF), Canadian Spatial Reference
System (CSRS), National Spatial Reference System (NSRS), etc) adopt coordinate reference
frames close to the ITRF at a certain epoch. They drift away over time due to the tectonic
movement of the continental plates. The GNSS SSR orbit correction messages allow transmission
of orbit corrections in continental/regional CRS. In this case, it is not necessary for the rover to
perform the corresponding transformation. A regional datum is indicated by the Global/Regional
CRS Indicator, while the actual CRS will be identified by the configured stream of the service
provider or the RTCM Service CRS message.
Note: Using a regional coordinate system may result in degraded performance if the CRS has a
large scale factor (relative to ITRF).
The satellite orbit corrections must refer to a Satellite Reference Point (SRP). The satellite orbit
may be provided with respect to the satellite Center of Mass (CoM) or the Antenna Phase Center
(APC) of a specific carrier frequency.
When using APC, the SRP should be the average APC of the corresponding carrier frequency. In
the past, the APCs and PCVs for GPS and GLONASS L1, L2 frequencies have been treated as
being identical to the APC and PCV of the L1/L2 ionospheric free linear combination (L0), i.e.
APC(L1) = APC(L2) = APC(L0).
With the availability of new constellations, frequencies, and calibration procedures, a different
definition is necessary. For IGS-SSR the SRP is the APC of a primary frequency per GNSS (see
below Table 2). This approach avoids the definition of a specific linear combination and is
compatible with the old definition for GPS and GLONASS (L1 and L2 identical to L0).
If the SRP is not the CoM, the SRP for orbit corrections is the antenna phase center (APC) of the
reference frequency as listed in Table 2.
8
IGS SSR v1.00 2020-10-05
The SSR generation system shall correct all observations for receiver and satellite antenna PCVs
and GDVs.
APC of
GNSS
Reference Frequency
GPS L1
GLONASS G1
Galileo E1
SBAS L1
BDS B1-2
QZSS L1
NavIC/IRNSS L5
Services using the IGS SSR format can refer to either the satellite CoM or the average APC of a
GNSS specific reference frequency (see Table 2).
For an APC reference a user needs to correct for PCO, PCV and GDV as a function of nadir angle
and azimuth (in the satellite reference frame) by utilizing frequency dependent correction tables.
To do that, the PCOs as given in the satellite antenna correction for all frequencies have to be
reduced by the PCO of the reference frequency. Fortunately, the azimuth dependent satellite PCVs
are expected to be small (in the order of 1 mm) and can thus be neglected for most applications. The
correction of PCV and GDV therefore requires only the nadir distance of the ray to the receiving
antenna. This angle can be computed without knowledge of the correct satellite attitude model.
For APC related services the IGS Generic SSR Messages (IGM001-IGM007) are sufficient for
successful and interoperable operation. It is recommended to use antenna corrections from the
most recent IGS antenna correction model (e.g. igs14.atx file).
In case of services utilizing the CoM as SRP, additional information and computation is required
on the client side. The APC referenced observations must be corrected by the CoM-APC offset
vector. Since these offsets can amount to meters for the Z-offset and several decimeters for the
X/Y offsets for some satellite types, an azimuth dependency is present. This requires the
9
IGS SSR v1.00 2020-10-05
knowledge of the correct attitude of the satellite. It should be noted that the correct yaw angle of
the satellite might not be accessible at all times.
In order to allow consistent and interoperable services utilizing CoM as SRP, IGS SSR messages
to transport the full satellite attitude with reliability flags as well as the CoM/APC offsets are
required. Such messages are under development.
The SSR clock correction message contains the parameters to compute the clock correction 𝛿𝐶
applied to the broadcast satellite clock. The polynomial representation describes the clock
differences for a certain time period. The sign definition of the corrections is
𝛿𝐶
𝑡𝑠𝑎𝑡𝑒𝑙𝑙𝑖𝑡𝑒 = 𝑡𝑏𝑟𝑜𝑎𝑑𝑐𝑎𝑠𝑡 −
𝑆𝑝𝑒𝑒𝑑 𝑜𝑓 𝑙𝑖𝑔ℎ𝑡
with
𝛿𝐶 = 𝐶0 + 𝐶1 (𝑡 − 𝑡0 ) + 𝐶2 (𝑡 − 𝑡0 )2
with
𝑡 time
𝑡0 reference time obtained from SSR Clock Correction message
𝐶𝑖 polynomial coefficients from SSR Clock Correction message, i = {0, 1, 2}
10
IGS SSR v1.00 2020-10-05
A High Rate Clock SSR message supports higher resolution of the clock state and higher update
rates. Both constituents, Clock Message and High Rate Clock Message, define the complete state
of the satellite clock. The High Rate Clock correction is added to the corresponding clock
correction.
The complete SSR clock correction term 𝛿𝐶 (refer to RTCM 10403.3) with the High Rate Clock
correction term reads
𝛿𝐶 = 𝐶0 + 𝐶1 (𝑡 − 𝑡0 ) + 𝐶2 (𝑡 − 𝑡0 )2 + 𝛿𝐶𝐻𝑖𝑔ℎ 𝑅𝑎𝑡𝑒
with
Satellite clocks are determined from ionospheric free signals derived from observations used by
the service provider. Such observations are affected by delays introduced in the satellite hardware
(code biases). For example, GPS broadcast clocks are referenced to the ionospheric free linear
combination of the P codes on L1 and L2, ignoring any code biases of these signals. For SSR, the
selection of signals used to generate the satellite clock corrections and the treatment of code biases
are left to the service provider. The service provider shall ensure a consistent transmission of clock
and code bias parameters. A rover must then consistently apply the code biases and clock
corrections.
11
IGS SSR v1.00 2020-10-05
In order to avoid Bias Datum conventions for the definition of the SSR standards, the SSR satellite
bias correction messages contain absolute signal biases. I.e., for each supported signal component,
a bias parameter needs to be transmitted. This increases the number of bits required for one
additional parameter but provides maximum flexibility for a service provider to select supported
signals and Bias Datum definitions. The code and phase bias messages contain absolute values,
but also enable the alternative use of DCBs by setting one of the biases to zero.
In the literature it is sometimes stated that more than one bias can be omitted or set to zero because
it cannot be distinguished from the ionospheric effect. However, this is only true if the ionosphere
is treated or modeled as an individual satellite dependent line of sight bias. In case of more
sophisticated ionospheric modelling this assumption no longer holds. With respect to this, note the
Dispersive Bias Consistency Indicator in the SSR phase bias messages.
It is further stated sometimes that phase biases can be omitted because they cannot be distinguished
from phase ambiguities. This however only holds for applications where no ambiguity resolution
is intended and only for short periods of time because a carrier phase ambiguity is constant over
time, but the phase bias is not (varies with temperature).
An ideal SSR generation process would rigorously estimate all SSR parameters consistently from
all observations of a continuously operating reference stations (CORS) network utilizing an
adequate filtering or adjustment tool. All pseudorange and phaserange observations would be
processed simultaneously to estimate all geometric (mainly orbits, clocks, troposphere, receiver
clocks and site coordinates) as well as dispersive (ionosphere) parameters plus all satellite and
receiver dependent signal biases and carrier phase ambiguities. Practically such a system is
difficult to implement, due to the very high number of unknown parameters as well as usually very
high number of observations. Furthermore, the formulation of functional and stochastic models is
a challenging task.
To overcome these problems, different strategies to reduce the number of parameters and/or
observations to be processed are in use. One strategy is the processing of linear combinations of
observations. For example, the observation equation for the ionosphere free (IF) linear
combinations (LC) of observables does not contain ionospheric parameters and their usage can
thus drastically reduce the number of unknowns. On the other hand, geometry free (GF) linear
combinations do not depend on geometric parameters and allow the estimation of ionospheric
models. Linear combinations with both IF and GF property (IF+GF) only depend on signal biases
and carrier phase ambiguities. One important application for the IF+GF linear combinations is the
ambiguity resolution for wide lanes, often called the Melbourne-Wübbena (MW) linear
combination.
Any linear combination of observations contains the corresponding linear combination of satellite
and receiver signal biases. The SSR generating process must take care that the satellite signal
12
IGS SSR v1.00 2020-10-05
biases are consistent with all the other SSR parameters, especially if the SSR stream is a
combination of the results from different partial filters or adjustments.
For example, after correcting the observations of a client application (rover) with the received SSR
signal biases (and phase wind-up and antenna PCO/PCV), different IF LCs of PR and/or PHR
observations can be computed for dual- or multi-frequency observations. Except for a common
bias for all satellites with the same LCs, remaining carrier phase ambiguities and estimation and
observation inaccuracies, all such IF LCs shall provide the same geometric (IF) PR/PHR. Hence,
a plot of between satellites single differenced and time differenced values of the difference of two
such IF LCs shall not contain any bias or ambiguity, but just observation and estimation noise and
remaining multipath effects.
Note that in both examples above uncorrected higher order ionospheric effects may also cause
some differences.
The usage of single or double differences is another example for the reduction of the number of
parameters in the system. The consistency requirement holds here if the results of corresponding
filters are combined with non-differenced results within the same SSR stream.
Depending on the processing model, SSR services may differ concerning the properties of the
estimates for satellite signal biases.
The phase bias messages contain two consistency indicators (MW Consistency Indicator and
Dispersive Bias Consistency Indicator).
The MW Consistency Indicator indicates the coherent use of code and phase biases for MW LCs.
If this indicator is set to 1, the MW combination of observations corrected by the transmitted code
and phase biases are bias free (except for a common bias for all satellites) and thus can be used to
determine the integer wide-lane ambiguity directly.
The Dispersive Bias Consistency Indicator shall be set to zero (0) if the corrected client PR and
PHR may contain arbitrary ionosphere because the service providers Satellite Bias Datum
constrains not only one but two biases. A client application cannot assume to obtain true
13
IGS SSR v1.00 2020-10-05
ionosphere after applying bias corrections to its observations. True slant TEC (STEC) is e.g.
required to compute higher order ionosphere corrections.
The satellite code and phase bias message use a GNSS Signal and Tracking Mode Identifier to
describe the actual signal properties. The GNSS Signal and Tracking Mode Identifier maps the
RINEX 3.x observation types into a more compact storage scheme of integer indices. The RINEX
observation types use type (code, carrier, etc.), band (L1, L2, etc.) and attribute (tracking mode).
The band and tracking mode attribute are required to account for the variety of signal tracking.
The yaw angle is used for computation of phase wind-up correction (Wu et al. 1993). The yaw
angle describes the angle between the yaw origin and the satellites x-axis as a positive rotation
around the nadir-pointing z-axis.
The user shall not assume that the yaw angle describes the true orientation of the satellite. Yaw is
provided only for consistent computation of phase wind-up corrections. The service provider must
choose an appropriate satellite yaw model to ensure consistent and continuous provision of satellite
clock correction, code and phase biases and yaw angle. A yaw of 0 (zero) is a valid value and shall
be applied for the wind-up correction by users receiving it. Since the provided yaw angle may not
describe the true satellite yaw, the yaw cannot safely be used to convert CoM to SRP orbits or vice
versa.
For the effective use of satellite phase biases any discontinuities and the integer property of the
phase bias are of interest. The discontinuity and integer property of the phase biases is reported for
every satellite phase signal.
In addition to the satellite phase signals, the integer property of wide-lane linear combinations is
reported for up to three groups. The wide-lane groups simplify the bookkeeping of resolved wide-
lanes. Each one of the pair of signals of a resolved wide-lane is assigned to the same group. With
every satellite frequency the number of wide-lanes combinations increases. For two or three
frequencies one group, for four or five frequencies two groups and for six frequencies at most three
groups are required to track the possible different wide-lane combination. For details refer to
IDF030.
It should be noted that due to the linear dependency of satellite and receiver biases, the SSR biases
of a specific signal for all satellites may have a common time varying offset (bias), which affects
the estimate of the corresponding rover receiver bias but does not harm the positioning solution.
Such common offsets/biases may behave different from physical receiver biases. Especially the
time variations may be larger than physical bias variations and may not be continuous. A rover
application should thus not make any assumption about the variation of common biases or
correlation of common biases for different signals.
14
IGS SSR v1.00 2020-10-05
A special consideration is necessary for GLONASS frequency division multiple access (FDMA)
signals. The GLONASS receiver dependent Code Phase Biases (CPB) are important in conjunction
with carrier phase ambiguity resolution. The service provider shall ensure that for all processed
GLONASS FDMA observations the CPBs are corrected before generating GLONASS satellite
phase bias estimates. Furthermore, the FDMA code biases may contain significant satellite or
channel dependent individual contribution. The magnitude of such individual biases is highly
dependent on the signal tracking hard- and software, i.e. the antenna, receiver type and firmware.
In case of a network with mixed CORS types, the service provider is responsible for the correct
handling of different individual satellite/channel biases from different tracking hard- and software.
If the remaining “service” biases can be leveled to a certain tracking hardware and software, the
service provider may indicate this through a standard RTCM MT1033 “Receiver and Antenna
Descriptor”.
The provider shall support as many signals as possible and must report biases for all supported
signals even if they are zero and/or used as Bias Datum. A rover can consistently use signals for
which a bias is transmitted. It is not reliable for a rover to use a signal without retrieving a
corresponding bias from the data stream. The coherent use of code and phase signals relies on the
availability of information for both phase and code biases.
The sign definition for the biases reported in the SSR satellite bias messages is that the biases
should be added to the PR or PHR measurements of the corresponding signal to get corrected
ranges.
15
IGS SSR v1.00 2020-10-05
Note: A spherical harmonic expansion with high degree and order may suffer a weak
representation of the true VTEC in areas with weak observation coverage. It is the responsibility
of the service provider to take countermeasures to avoid unreasonable VTEC values. Negative
VTEC values must be ignored by the rover and shall be replaced by 0.0.
In order to simplify the algorithm and to avoid heavy computational load, the spherical harmonics
are defined for a spherical Earth model, with a mean radius of 6370 km. Latitude, longitude, height
as well as azimuth and elevation of the satellites are defined with respect to this model instead of
an ellipsoidal representation.
Note: For the computation of the ECEF satellite position’s azimuth and elevation, the ECEF
satellite coordinates computed at the time of signal transmission at the satellite shall be rotated to
the epoch of signal reception at the receiver to account for the Sagnac effect caused by Earth
rotation and signal propagation.
The cosine coefficients C and sine coefficients S are represented for a specific degree 𝑁 and order
𝑀 ≤ 𝑁 of a series of spherical harmonic functions to describe the VTEC in total electron content
𝑒𝑙𝑒𝑐𝑡𝑟𝑜𝑛𝑠
units (TECU) (1 TECU = 1016 ).
𝑚2
min (𝑛,𝑀)
𝑉𝑇𝐸𝐶(𝜑𝑃𝑃 , 𝜆𝑃𝑃 ) = ∑𝑁
𝑛=0 ∑𝑚=0 (𝐶𝑛𝑚 cos 𝑚𝜆𝑆 + 𝑆𝑛𝑚 sin 𝑚𝜆𝑆 )𝑃𝑛𝑚 (sin 𝜑𝑃𝑃 )
with
𝑁 degree of spherical expansion
𝑀 order of spherical expansion
𝑛, 𝑚 indices
𝐶𝑛𝑚 cosine coefficients for the layer [TECU]
𝑆𝑛𝑚 sine coefficients for the layer [TECU]
𝜆𝑆 mean sun fixed and phase shifted longitude of ionospheric pierce point
for the layer modulo 2𝜋 [radians]
𝜆𝑃𝑃 longitude of ionospheric pierce point for the layer [radians]
t SSR epoch time of computation epoch modulo 86400 [s]
𝜑𝑃𝑃 geocentric latitude of ionospheric pierce point for the layer [radians]
𝑃𝑛𝑚 () fully normalized associated Legendre functions
The mean sun fixed longitude phase shifted by 2 h to the approximate TEC maximum at 14:00
local time (resp. 50400 s) is computed by:
16
IGS SSR v1.00 2020-10-05
Longitude 𝜆𝑃𝑃 and latitude 𝜑𝑃𝑃 refer to the position of the ionospheric pierce point at the height
of the ionospheric layer. The ionospheric pierce point is the intersection of a straight line from the
rover position to the satellite and a sphere with the height of the ionospheric layer above the
spherical Earth model.
The position of the pierce point in the spherical Earth model is computed in radians, for example,
by:
with
The angle 𝜓𝑃𝑃 is the spherical Earth’s central angle between rover position and the projection of
the pierce point to the spherical Earth’s surface. It is computed in radians by:
𝑅 +ℎ
𝜓𝑃𝑃 = 𝜋⁄2 − 𝐸 − 𝑎𝑟𝑐𝑠𝑖𝑛 ( 𝑅𝑒 +ℎ𝑅 cos 𝐸)
𝑒 𝐼
Under the following conditions, which apply for rover position and satellites visible at rover
position
17
IGS SSR v1.00 2020-10-05
In all other cases the pierce point longitude is computed in radians according to
𝑉𝑇𝐸𝐶
𝑆𝑇𝐸𝐶𝑖 = sin(𝐸+𝜓𝑃𝑃 )
The total STEC is computed by the sum of the individual STECi for each layer i.
The influence of the ionosphere on the observations for a specific frequency f [Hz] is a
pseudorange delay in m:
or a phaserange advance in m:
All messages may be split over several messages using the SSR Multiple Message Indicator. It
increases the number of e.g. satellites to be transmitted, but it also enables the receiver to begin
immediately processing of data after decoding the first message of a multiple message type.
18
IGS SSR v1.00 2020-10-05
Sub-Type Notes
IGM Type GNSS IGM Type Name
Messages
IM000– IM020 Reserved
IGM01 GPS SSR Orbit Correction IM021
IGM02 GPS SSR Clock Correction IM022
SSR Combined Orbit and
IGM03 GPS IM023
Clock Correction
SSR High Rate Clock
IGM04 GPS IM024
Correction
IGM05 GPS SSR Code Bias IM025
IGM06 GPS SSR Phase Bias IM026
IGM07 GPS SSR URA IM027
IM028 – IM040 Reserved for GPS
IGM01 GLONASS SSR Orbit Correction IM041
IGM02 GLONASS SSR Clock Correction IM042
19
IGS SSR v1.00 2020-10-05
Sub-Type Notes
IGM Type GNSS IGM Type Name
Messages
SSR Combined Orbit and
IGM03 GLONASS IM043
Clock Correction
SSR High Rate Clock
IGM04 GLONASS IM044
Correction
IGM05 GLONASS SSR Code Bias IM045
IGM06 GLONASS SSR Phase Bias IM046
IGM07 GLONASS SSR URA IM047
IM048 – IM060 Reserved for GLONASS
IGM01 Galileo SSR Orbit Correction IM061
IGM02 Galileo SSR Clock Correction IM062
SSR Combined Orbit and
IGM03 Galileo IM063
Clock Correction
SSR High Rate Clock
IGM04 Galileo IM064
Correction
IGM05 Galileo SSR Code Bias IM065
IGM06 Galileo SSR Phase Bias IM066
IGM07 Galileo SSR URA IM067
IM068 – IM080 Reserved for Galileo
IGM01 QZSS SSR Orbit Correction IM081
IGM02 QZSS SSR Clock Correction IM082
SSR Combined Orbit and
IGM03 QZSS IM083
Clock Correction
SSR High Rate Clock
IGM04 QZSS IM084
Correction
IGM05 QZSS SSR Code Bias IM085
IGM06 QZSS SSR Phase Bias IM086
IGM07 QZSS SSR URA IM087
IM088 – IM100 Reserved for QZSS
IGM01 BDS SSR Orbit Correction IM101
IGM02 BDS SSR Clock Correction IM102
SSR Combined Orbit and
IGM03 BDS IM103
Clock Correction
SSR High Rate Clock
IGM04 BDS IM104
Correction
20
IGS SSR v1.00 2020-10-05
Sub-Type Notes
IGM Type GNSS IGM Type Name
Messages
IGM05 BDS SSR Code Bias IM105
IGM06 BDS SSR Phase Bias IM106
IGM07 BDS SSR URA IM107
IM108 – IM120 Reserved for BDS
IGM01 SBAS SSR Orbit Correction IM121
IGM02 SBAS SSR Clock Correction IM122
SSR Combined Orbit and
IGM03 SBAS IM123
Clock Correction
SSR High Rate Clock
IGM04 SBAS IM124
Correction
IGM05 SBAS SSR Code Bias IM125
IGM06 SBAS SSR Phase Bias IM126
IGM07 SBAS SSR URA IM127
IM128 – IM140 Reserved for SBAS
IM141 – IM160 Reserved for NavIC/IRNSS
IM161 – IM200 Reserved
SSR Ionosphere VTEC
- GNSS IM201
Spherical Harmonics
IM202 – IM254 Reserved
Table 6 is a list of general IGS SSR Sub-Type Messages, which do not depend on an IGS Generic
SSR Message type.
21
IGS SSR v1.00 2020-10-05
7. DATA FIELDS
7.1 Data Field Summary
Table 7 describes the additional data fields required to support the IGS SSR messages. The data
field number is prefixed with “IDF” to indicate the IGS SSR relevant data fields. RTCM3 data
field numbers are prefixed with “DF”.
22
IGS SSR v1.00 2020-10-05
23
IGS SSR v1.00 2020-10-05
24
IGS SSR v1.00 2020-10-05
25
IGS SSR v1.00 2020-10-05
26
IGS SSR v1.00 2020-10-05
27
IGS SSR v1.00 2020-10-05
28
IGS SSR v1.00 2020-10-05
29
IGS SSR v1.00 2020-10-05
30
IGS SSR v1.00 2020-10-05
31
IGS SSR v1.00 2020-10-05
32
IGS SSR v1.00 2020-10-05
7.2.1 GPS
The relativity correction has to be applied for GPS according to the GPS ICD to compute
𝑡𝑏𝑟𝑜𝑎𝑑𝑐𝑎𝑠𝑡 . The relativistic correction term ∆𝒕 𝑇 to be used is based on the orbital elements.
Note: The GPS ICD contains two different formulas for the relativistic correction. A second
formula is based on ECEF satellite broadcast position and velocity vector. The ICD states identity
between them, but there is actually a numerical difference.
7.2.2 GLONASS
The computation of satellite positions and satellite clock offsets from broadcast data must be
consistent. The GLONASS ICD refers to two different algorithms. For SSR, the GLONASS
broadcast ephemeris computation “Simplify of algorithm for re-calculation of ephemeris to current
time” must be used.
The relativistic effects are already accounted for in the broadcast clock parameters for GLONASS.
Note: The numerical integration of satellite coordinates from broadcast ephemeris should be
performed with sufficient accuracy. For instance, the Runge-Kutta method of 4th order with the
step of numerical integration that does not exceed 30 seconds can be used to obtain a precision of
0.1 mm.
7.2.3 Galileo
According to the Galileo OS SIS ICD, Issue 1.1, September 2010, there are two different clock
representations for Galileo satellite clock corrections, namely the F/NAV and the I/NAV clock.
The clocks are transmitted for different services and signals. F/NAV is for Dual-frequency (E1,
E5a) or Single-frequency E5a services. I/NAV is available for Dual-Frequency (E1, E5b), Single-
frequency E5b and Single-frequency E1. The clocks will be derived from dual frequency
ionosphere free linear combinations of observables on E1/E5A or E1/E5b respectively.
For both, F/NAV and I/NAV, clock polynomial coefficients and a reference time are provided by
Galileo. However, there is no information on consistency of both clocks in conjunction with group
delay parameters transmitted.
33
IGS SSR v1.00 2020-10-05
Clock corrections in RTCM-SSR are related to a broadcast reference clock. The I/NAV clock has
been chosen as the reference clock for RTCM Galileo SSR correction.
The relativity correction has to be applied for Galileo according to Galileo ICD to compute
𝑡𝑏𝑟𝑜𝑎𝑑𝑐𝑎𝑠𝑡 . The relativistic correction term ∆𝒕 𝑇 to be used is based on the orbital elements.
7.2.4 QZSS
The relativity correction has to be applied for QZSS according to QZSS ICD to compute 𝑡𝑏𝑟𝑜𝑎𝑑𝑐𝑎𝑠𝑡 .
The relativistic correction term ∆𝒕 𝑇 to be used is based on the orbital elements.
7.2.5 BDS
The relativity correction has to be applied for BDS according to BDS ICD to compute 𝑡𝑏𝑟𝑜𝑎𝑑𝑐𝑎𝑠𝑡 .
The relativistic correction term ∆𝒕 𝑇 to be used is based on the orbital elements.
7.2.6 SBAS
The Satellite-based Augmentation System (SBAS) broadcast ephemeris is identified by the Issue
of Data Navigation (IODN). The IODN is defined as the first 8 bits after the message type 9, called
IODN in RTCA DO229, Annex A and Annex B and called spare in Annex C (RINEX).
Although, the IODN is called spare in the RTCA DO229, it is used for the IGS SSR correction
messages as SBAS broadcast ephemeris indicate support and increasing count of the IODN.
There is no relativity correction for SBAS. Any relativity effects are removed by the control center
of the SBAS signal according to RTCA DO229.
34
IGS SSR v1.00 2020-10-05
Table 8 : Header Part of the GNSS SSR Orbit Correction Message IGM01
Table 9 : Satellite Specific Part of the GNSS SSR Orbit Correction Message IGM01
The GNSS satellite clock correction message content is defined in Table 10 and Table 11,
consisting of a header and satellite specific part.
35
IGS SSR v1.00 2020-10-05
Table 10 : Header Part of the GNSS SSR Clock Correction Message IGM02
Table 11 : Satellite Specific Part of the of the GNSS SSR Clock Correction Message IGM02
The GNSS combined satellite orbit and clock correction message content is defined in Table 12
and Table 13, consisting of a header and satellite specific part.
36
IGS SSR v1.00 2020-10-05
Table 12 : Header Part of the GNSS SSR Combined Orbit and Clock Correction Message
IGM03
Table 13 : Satellite Specific Part of the GNSS SSR Combined Orbit and Clock Correction
Message IGM03
The GNSS satellite high rate clock correction message content is defined in Table 14 and Table
15, consisting of a header and satellite specific part.
37
IGS SSR v1.00 2020-10-05
Table 14 : Header Part of the GNSS SSR High Rate Clock Correction Message IGM04
Table 15 : Satellite Specific Part of the GNSS SSR High Rate Clock Correction Message
IGM04
38
IGS SSR v1.00 2020-10-05
Table 16 : Header Part of the GNSS SSR Code Bias Message IGM05
Each SV Block consists of a satellite specific part immediately followed by a code specific part
per processed code bias for the corresponding satellite.
Table 17 : Satellite Specific Part of the of the GNSS SSR Code Bias Message IGM05
Table 18 : Code Specific Part of the of the GNSS SSR Code Bias Message IGM05
39
IGS SSR v1.00 2020-10-05
The GNSS satellite phase bias correction message content is defined in Table 19, Table 20 and
Table 21, consisting of a header, satellite and bias specific part.
Table 19 : Header Part of the GNSS SSR Phase Bias Message IGM06
Each SV Block consists of a satellite specific part immediately followed by a phase specific part
per processed phase bias for the corresponding satellite.
Table 20 : Satellite Specific Part of the of the GNSS SSR Phase Bias Message IGM06
40
IGS SSR v1.00 2020-10-05
Table 21 : Phase Specific Part of the of the GNSS SSR Phase Bias Message IGM06
41
IGS SSR v1.00 2020-10-05
Table 23 : Satellite Specific Part of the GNSS SSR URA Message IGM07
Table 24 : Header Part of the SSR Ionosphere VTEC Spherical Harmonics IM201
Every ionospheric layer model has an Ionosphere Layer Header (Table 25) and two subsequent
model parts with Cosine Coefficients (Table 26) and Sine Coefficients (Table 27).
42
IGS SSR v1.00 2020-10-05
Table 26 : Model Part of the SSR Ionosphere VTEC Spherical Harmonic Cosine Coefficients
IM201
The number of sine terms of the Sine Coefficients Part depends on the order and degree of the
spherical harmonic expansion. Refer to the corresponding Data Field Note for IDF040.
Table 27 : Model Part of the SSR Ionosphere VTEC Spherical Harmonic Sine Coefficients
IM201
The number of harmonic coefficients for one layer is defined by the corresponding degree 𝑁 and
order 𝑀. The overall number of cosine and sine coefficients for each layer is computed according
to:
The sequence of the harmonic coefficients in the message is first the cosine terms C followed by
the sine terms S within each group:
43
IGS SSR v1.00 2020-10-05
I.e. with degree 𝑁 = 3, order 𝑀 = 2 the sequence of the harmonic coefficients within one layer
with corresponding index n, m reads
44
IGS SSR v1.00 2020-10-05
9 REFERENCES
IERS Conventions (2010). Gérard Petit and Brian Luzum (eds.). “IERS Technical Note 36“,
Frankfurt am Main, Verlag des Bundesamts für Kartographie und Geodäsie, 2010. 179 pp., ISBN
3-89888-989-6
Schmitz, M. (2012). “RTCM State Space Representation Messages, Status and Plans.”
PPP-RTK & Open Standards Symposium, March 12-13, 2012, Frankfurt, Germany
Wu, J.T., S.C. Wu, G.A. Hajj, W.I. Bertiger, S.M. Lichten (1993). „Effects of Antenna Orientation
on GPS Carrier Phase”. Manuscripta Geodaetica, 18, 91-98
Wübbena, G. (2012): “RTCM State Space Representation (SSR) Overall Concepts Towards PPP-
RTK.” PPP-RTK & Open Standards Symposium, March 12-13, 2012, Frankfurt, Germany
Document: RINEX, The Receiver Independent Exchange Format, Version 3.04, November 23,
2018
Document: RTCM Standard 10403.3, Differential GNSS (Global Navigation Satellite Systems)
Services – Version 3, October 7, 2016
45