On the equivalence of ionospheric-free, differenced, and undifferenced, uncombined GNSS network processing

Global Navigation Satellite System (GNSS) observations are commonly processed using various established methods, including the processing of ionospheric-free (IF) linear combinations (LC), differences of observations as well as the processing of undifferenced, uncombined (UDUC) observations. The most general and flexible approach to GNSS processing is widely regarded to be the UDUC approach, as this approach is based on the raw observation equations of potentially all available observations. The IF approach uses IF-LC obtained by observations on different frequencies but the same receiver-satellite link to eliminate the ionospheric slant delay. The differencing approach also uses LC of observations, with the distinction that observations of the same signal but different receivers and/or satellites are being used. The purpose of differencing is the elimination of satellite/receiver clocks and biases and, in some cases, reducing or even eliminating atmospheric delays. We aim to uncover the implicit model assumptions made when using various IF and differencing approaches and in what case they are equivalent to the processing of the UDUC observations. This is achieved by introducing a reformulation method, which is then applied to the UDUC observations of code division multiple access GNSS to obtain the functional models of various IF and differencing approaches. The underlying assumptions in this reformulation can then be identified. The results of this theoretical contribution will provide insight into the most appropriate method for processing GNSS observations in different cases and what implicit assumptions are being made when the respective method is being used.


Introduction
Precise point positioning (PPP) (Malys and Jensen 1990;Zumberge et al. 1997;Kouba and Héroux 2001;Kouba et al. 2017) is a technique used to process GNSS observations to achieve highly accurate point positions.PPP can be divided into the network and the user side.The network side consists of globally distributed receivers observing one or multiple GNSS.With this network, precise satellite orbits and clocks can be computed and provided to the user for highly accurate point solutions.In order to generate accurate satellite orbits and clocks on the network side and accurate point positions on the user side, code and phase observations are required.PPP-real-time kinematic (PPP-RTK) is a variation of PPP that includes the satellite phase bias parameters in addition to the satellite orbits and clocks.This enables the network and the user to estimate integer ambiguities through proper satellite phase bias and ambiguity datum handling (Geng et al. 2009;Laurichesse et al. 2009;Teunissen et al. 2010).PPP-RTK is a more general formulation of PPP, as it considers existing phase biases of receivers and satellites.For a critical review of different PPP-RTK methods, see Teunissen and Khodabandeh (2015).Here, we focus on the network side of PPP-RTK, as this is the most general case that includes the user side as well, considering that a user can be interpreted as another receiver in the network.
One commonly used method for processing dual-frequency observations is the use of IF-LC instead of the original observations.In the IF-LC, the ionospheric slant delay has been eliminated and does not have to be accounted for in the functional model.Two methods have emerged regarding dual-frequency observations using the IF-LC.The first method, often used in the conventional PPP approach where 173 Page 2 of 10 no ambiguity fixing is desired, uses one IF-LC using code observations and another IF-LC only using phase observations.For the second method, an additional third IF-LC is introduced, the Hatch-Melbourne-Wübbena LC (HMW-LC) (Hatch 1982;Melbourne 1985;Wübbena 1985) and is more commonly used in PPP-RTK see Ge et al. (2008), Collins et al. (2010), and Geng et al. (2019).
Another method of processing GNSS observations is double differencing (DD), which is commonly used in relative positioning but also used to process global networks and generate satellite orbits (Dach et al. 2021).In the DD approach, observations of the same signal are differenced between receivers and satellites, so satellite and receiver clocks and biases cancel out.When this approach is used for short baselines, it can also be assumed that atmospheric delays and orbit errors are eliminated or significantly reduced.Single differencing (SD) is also a well-known approach to GNSS processing.In the SD approach, observations are differenced between receivers or between satellites to eliminate the satellite clocks and biases or the receiver clocks and biases, respectively.
Finally, the approach of processing undifferenced and uncombined (UDUC) observations at the network and user levels (Teunissen et al. 2010;Schönemann et al. 2011;Zhang et al. 2011) has emerged.This is widely regarded as the most general approach, as the original observations are being processed, and no parameters are eliminated by LC or differencing of observations.It is regarded as the most flexible approach, offering a general formulation of the functional model while maintaining the original physical and stochastic characteristics of the observations.Schaffrin and Grafarend (1986) introduced a concept of nuisance parameter elimination at an observational level and applied this theory to single-frequency phase observations.It was found that SD and DD are equivalent reformulations of the system of equations of single-frequency phase observations in which the satellite and/or receiver clock are eliminated.In contrast to the previously mentioned work, the reformulation method we present is applied to the general UDUC network equations, including multi-signal code and phase observations and considers existing rank deficiencies.
First, an overview of the processing of UDUC observations in the case of GNSS with CDMA is provided.This is done by briefly introducing the code and phase observation equations of a GNSS network, the resulting design matrix, and its rank deficiencies.For this, the relevant aspects of the S-system theory (Teunissen 1985) and the contribution of rank deficiencies in the UDUC processing of GNSS networks by Odijk et al. (2015) are used in this work.Then a reformulation of the system of network equations is presented.In subsequent chapters, we demonstrate that the two IF-LC and the differencing approaches are specific applications of the presented reformulation methods.By examining the implicit assumptions made in these approaches, we can show their mathematical equivalence to the respective results of the processing of UDUC observations.

UDUC network equations
The linearized code and phase observation equations for a receiver r, satellite s, and signal j can be formulated with (Hofmann-Wellenhof et al. 2007) in case of CDMA GNSS.Here, E(⋅) denotes the expected value of the inherently stochastic observations.It is assumed that the observations have been corrected a priori with appropriate models (Kouba et al. 2017).Thus, terms usually not estimated as parameters in the context of PPP-RTK, such as phase wind-up, phase center offsets, and variations, are excluded in (1) and ( 2).
Given a network with r = 1, … , n receivers, s = 1, … , m satellites observing j = 1, … , f signals and i = 1, … , p epochs, we can formulate a system of equations with L + v = A ⋅ x , with L containing the linearized code and phase observations, v their residuals, A the design matrix and x the unknown parameters.The general design matrix A of an all-in view, i.e., all receivers observe all satellites and signals, network can be formulated with (Odijk et al. 2015) where here e n denotes the n-vector of ones, I n ∈ ℝ n×n the identity matrix, ⊗ the Kronecker matrix product, Λ is a diagonal matrix containing the wavelength of the observed signals, i.e., Λ = diag 1 , … , f and = 1 2 … f T .The design matrix A geo references the geometric parameters (receiver/satellite coordinates and tropospheric delays) and is not further specified, as it is not of further relevance hereafter.All other submatrices have been created assuming (1) (2) (5) Then, all code observations (i) are ordered in the same man- ner and finally The parameter vector x is divided into subvectors of specific where x geo contains all geometric parameters, i.e., receiver coordinates and tropospheric delays, x rec all receiver clocks and biases, x sat all satellite clocks and biases, x ion all iono- spheric slant delays, and finally x amb all ambiguities.The design matrix is formulated similarly, with the submatrices referencing the respective subvector.
The design matrix A is of rank(A) = dim (R(A)) = r ≤ k with R(A) denoting the range or column space.The rank deficiency k − r of matrix A is equal to the dimension of its null space, i.e., dim(N(A)) = def(A) .Here, N denotes the null space of a matrix.We can find a basis matrix of N(A) with V ∈ ℝ k×k−r .Then, R(V) = N(A) and A ⋅ V = 0 .With a basis matrix S that has a complementary range space to V, i.e., R(S) ⊗ R(V) = ℝ k , we obtain the parameter vector x as a unique linear combination of the form x = S ⋅ ã + V ⋅  .We can then reformulate the system of equations to The basis matrix S is not unique, and its choice defines what estimable parameters one solves.The matrix S ⊥ is defined as a basis matrix of the orthogonal complement of R(S), i.e., S ⊥ T S = 0 .Its transpose S ⊥ T ∈ ℝ (k−r)×k can be interpreted as containing the (k − r) restrictions imposed on the param- eter vector, i.e., S ⊥ T x = 0 (Khodabandeh and Teunissen 2019), making A T S ⊥ T a full rank matrix.
The rank deficiencies in A can be described by the iden- tification of the rank deficiencies in the submatrices and the additional rank deficiencies occurring between the submatrices, as explained in the following Corollary.

Corollary 1 ( Rank deficiencies in submatrices) Let matrix
(6) , in all other cases, we can extend V * to a basis matrix Note that all vectors in W i are linearly independent of V i .(Proof: Let us assume a vector w * 1 in W 1 can be expressed as a LC of vectors in V 1 .Then, we can find a matrix C , so that The rank deficiencies of A can therefore be interpreted as the r 1 and r 2 rank deficiencies of the submatrices A 1 and A 2 and the l between A 1 and A 2 rank deficiencies.We denote with the rank deficiency between matrices. In Odijk et al. (2015), the different types of rank deficiencies are specifically declared as.The size of these rank deficiencies in the scope of this research is understood to exist at each epoch, except for the rank deficiencies between ambiguities and biases, as ambiguities are introduced as constant in (5).Often used S-basis restrictions to eliminate these rank deficiencies are the "Common Clocks" (CC) S-basis.It is further classified into the CC-R and the CC-S S-basis.The CC-R S-basis uses the receiver clock and biases of one receiver as restrictions to eliminate rank deficiencies 1a and 1b, whereas the CC-S S-basis uses the mean satellite clock and biases of a satellite system.
Since none of the individual S-basis restrictions in the CC models span across multiple parameter groups, we can 173 Page 4 of 10 express S ⊥ T as a block diagonal matrix, with each matrix containing the restrictions imposed on a specific parameter group, as presented in the following Corollary.

Corollary 2 (Restrictions in GNSS processing
We can find a basis matrix Applying Corollary 2 to our specific GNSS parameter g r o u p s y i e l d s amb T , a n d c o n s e q u e n t l y, w e c a n f i n d S = diag S rec , S sat , S ion , S amb in a similar notation.The dimensions of the submatrices of S ⊥ T and S are defined as outlined in the corollary and vary depending on the number of S-basis restrictions imposed on each parameter group.

Transforming the network equations
This chapter is structured as follows.First, the theory required to achieve an equivalent reformulation of the UDUC network equations is presented.The presented reformulation takes existing rank deficiencies in the original design matrix and how they propagate into the transformed design matrix into account.This reformulation is then applied to obtain the two different IF, and the SD and DD processing approaches.

Theory
It is possible to transform the original system of equations as defined by the design matrix (3) of the UDUC observations with no loss of information.This is described in the following Lemma.

Lemma 1 (Reformulation of the Least Squares (LSQ) Problem)
The solution of a system of equations of the form is given with and is also known as the weighted LSQ solution with weight matrix P being the inverse of the variance covariance matrix (VCM) of the observations L.
An equal solution introducing a regular matrix V can be found with Lemma 1 can be summarized as follows: It is possible to use the LC L * = VL instead of the original observations L to obtain solution ã .The solution ã using observations L * is equal to the solution obtained by the original obser- vations L , as long as V is a regular matrix and the design and weight matrix of the observations are transformed as demonstrated.
The motivation to perform such a transformation is given in the following Lemma.

Lemma 2 (Partial LSQ solution)
The partial solution x 1 of the LSQ problem can be expressed with if A 3 is a regular matrix and P a symmetric and positive definite matrix.
Therefore, by applying an appropriate transformation, the matrix Ã * = V Ã can be transformed into the form of the design matrix presented in Lemma 2, allowing for the pre-elimination of certain parameters x 2 at an observational level.This enables us to obtain the partial solution x 1 using a reduced observation vector.This partial solution is mathematically equivalent to the partial solution obtained when all observations are considered.
A generalization of this approach, including possible complications arising due to rank deficiencies, is given in the following corollary, which serves as the foundation for the next chapters, as the following methods of GNSS processing can be shown to be an application thereof.Note that in Corollary 3, it is not assumed that rank A 1 A 2 = rank A 1 + rank A 2 like it has been assumed in Schaffrin and Grafarend (1986).
Corollary 3 ( Reformulation of rank-deficient LSQ) For a rank-deficient LSQ problem The admissible S-basis restrictions are chosen so and S ⊥ 2 ∈ ℝ r 2 ×k 2 .It follows that we can find We can then reformulate the LSQ problem with where matrix V 2 A 2 S 2 is a regular and V 1 A 1 S 1 a full rank matrix.The partial solution ã1 can be computed by solving the reduced LSQ problem, as shown in Lemma 2.
The condition regarding the S-basis restrictions can be interpreted as follows: The restrictions S ⊥ 2 T must exclusively eliminate the r rank deficiencies in A 2 , and the restrictions S ⊥ 1 T must eliminate the r 1 rank deficiencies occurring in A 1 , as well as the b rank deficiencies occurring between A 1 and A 2 (see Lemma 3).

Lemma 3 ( Null space after transformation) Let
be a basis matrix of N(A) , V i a basis matrix of A i , and F T 1 a basis matrix of N A T 2 .Then,

Ionospheric-free linear combinations
In the next chapters, the GNSS processing approach of formulating IF-LC is introduced and the equivalence to the processing of UDUC observations is demonstrated.First, the general case is presented, which provides an admissible and equivalent reformulation of the UDUC network equations as introduced with Eqs.(3-5) for any number of signals.Then, two common dual-frequency IF approaches are explored.
The first approach, commonly used in PPP, involves using one code and one phase IF-LC.The second approach, often used in PPP-RTK, involves using the HMW-LC in addition to the previously mentioned IF-LC.

General case
LC in the context of GNSS processing specifically refer to LC using only code and phase observations of a specific receiver-satellite link.They can be formulated with where l are the resulting LC and Important properties of LC might be: IF, geometric free (GF), ionospheric maintaining (IM), and geometry maintaining (GM).After introducing B = e f e f − ∈ ℝ 2×l , we can derive specific properties of a LC v ⋅ In the context of GNSS processing, LC are commonly used to eliminate the ionospheric slant delay present in the UDUC observations.For the given matrix B , we can find a matrix and applying the LC defined by V to the observations of all satellite-receiver links, we obtain the regular matrix W of the entire network: ion .The transformed design matrix Ã * then reads Since there exist no rank deficiencies in A ion , i.e., def A ion = 0 , we can find S ⊥ T = S ⊥ 1 T 0 and S = S 1 0 0 I (see Corollary 2).The full rank design matrix Ã can then be obtained with where W T IM ⋅ A ion ∈ ℝ pnm×pnm is a regular matrix.Thus, we can obtain all parameters, excluding the ionospheric slant delays, only by processing the IF observations (see Lemma 2).Both the CC-R the CC-S basis are valid choices for an IF-LC approach since ionospheric slant delays are not used as restrictions to eliminate the rank deficiencies occurring between A rec A sat and A ion .Those rank deficiencies are eliminated via restrictions on the receiver and satellite biases instead.
Note that the IF-LC combination matrix W T IF is a basis matrix of N A T ion when epoch-wise ionospheric slant delay parameters are assumed and l − 1 linear independent IF-LC for each receiver-satellite link are used.If additional temporal or spatial ionospheric constraints were introduced into the submatrix A ion , then W T IF would not be a basis matrix of the null space of the transposed partial design matrix that includes the temporal and/or spatial constraints.With the greater dimension of the null space, more observations would be required. (20) We can conclude that the LSQ solution of the IF-LC approach using l − 1 linear independent IF-LC for each receiver-satellite link is mathematically equivalent to the solution when evaluating the UDUC network equations when ionospheric slant delays are assumed to be epoch-wise parameters.We also found that both the CC-R and CC-S S-basis restrictions are viable options for the IF-LC approach.Given that the estimation of epoch-wise ionospheric slant delays is a common assumption when processing UDUC observations, the IF-LC approach presents a valid alternative that reduces the number of observations and parameters when compared to the UDUC approach.

Current practice in dual frequency IF processing
In PPP, it is common practice to process dual-frequency GNSS observations using one IF and GM code and phase LC.The LC matrix V can then be formulated with where As demonstrated in the previous chapter and explicitly mentioned in Teunissen (2020), it is evident that this is not a valid IF reformulation of the functional model described by (3) to (5).This is due to the requirement of three linearly independent IF LC for every receiver-satellite link when considering dual-frequency observations.
When reformulating the design matrix of the ionospheric delays to introduce different ionospheric slant delay parameters for code and phase observations, i.e., the transposed of the matrix becomes a basis matrix of N A T 2ion and a valid reformulation can be achieved according to Corollary 3. The conclusion is that the implicit model assumptions, when processing only two linear independent IF-LC, are epoch-wise ionospheric slant delay parameters that are different for code and phase observations.
For the purpose of PPP-RTK based on IF-LC, it is common practice to use the HMW-LC in addition to the previously mentioned LC (24).Then, a valid reformulation with pre-eliminated ionospheric slant delays as described in Corollary 3 is achieved, and integer ambiguity resolution can be performed analog to when the UDUC observations are ( 24) processed.Note that this approach does not require a reparameterization of the ambiguities.Integer estimable ambiguities can be achieved as described by Khodabandeh and Teunissen (2019), and an a priori wide-lane and narrow-lane reparameterization like it is often performed (Ge et al. 2008;Laurichesse et al. 2008) and is not required.The reparameterization of the ambiguities via an admissible ambiguity transformation (Teunissen 1995) can be done independently of the chosen IF-LC.The same applies to all other remaining parameters, including receiver and satellite phase biases.
They can be estimated in the same manner as in the general UDUC approach, with the only exception being that ionospheric slant delays cannot be used as S-basis restrictions.

Single differencing
SD refers to a differencing of observations between satellites (BSSD) or between receivers (BRSD) and is most commonly performed on the same signal.
Between receiver A and B , SD Δ * s AB,j can be formulated with Δ * s AB,j = * s A,j − * s B,j , where * is a wildcard for either a code or phase observations.Observations from different receivers but the same signal and satellite are differenced.In the resulting observation Δ * s AB,j , satellite clocks and biases of satellite s have been eliminated.BSSD can be formulated similarly with Δ * r,j , i.e., observa- tions of a single receiver and signal are differenced between the satellites s A and s B .Here, the receiver clocks and biases are eliminated in the resulting BSSD observation.In this chapter, both the BRSD and BSSD matrices are introduced, but for brevity, only the derivation of the implicit assumptions of BRSD is given as the derivation of the BSSD can be done similarly.
A general formulation of the BRSD and BSSD differencing matrix based on the observation order defined by ( 6) and ( 7) is given with: A well-known property of BRSD is that satellite clocks and biases are eliminated, i.e., and is a full rank matrix, it is a basis matrix of N A T sat .We can obtain a square and invertible matrix Analog to the presented method in Corollary 3, we now use F to reformulate the design matrix order of the submatrices in A has been changed for ease of notation) and obtain: with F r A sat S sat being a regular matrix if and only if rank defi- ciencies occurring between A 1 and A sat have been eliminated with restrictions imposed on parameters belonging to A 1 .
Neither the CC-R nor the CC-S S-basis restrictions can be used to obtain the partial solution since parameters in x sat are used as restrictions to eliminate rank deficiencies between A 1 and A sat , or more precisely: satellite biases are used to eliminate rank deficiencies of type 3b.A valid S-basis alternative can be found by using one ionospheric slant delay per satellite as a restriction instead.The corresponding S-basis is the CC-R-SD S-basis in the case of BRSD and the CC-S-SD basis in the case of BSSD given in Table 1.One can now solve for the remaining parameters only using the SD observations F SD ⋅ L .The solution is then equal to the parameters obtained if one had used all observations and the S-basis restrictions shown in Table 1.
With this, it has been established that the LSQ solution of the BRSD (BSSD) GNSS processing approach is mathematically equivalent to the solution of the general UDUC when satellite (receiver) clocks and biases are estimated as epoch-wise parameters.The SD approach pre-eliminates the respective parameters at an observational level, reducing the number of observations and parameters in the functional model.The pre-eliminated parameters can no longer be used to eliminate rank deficiencies, and the SD observation design matrix "inherits" the rank deficiencies between its parameters and the pre-eliminated parameters (Lemma 3).This has to be accounted for when choosing the S-basis restrictions.

Double differencing
DD can be expressed with ( 28) Table 1 Adjustments to the CC-R and CC-S basis Note that the adjustments for CC-R-DD and CC-S-DD include the previous adjustment above.The column rd# refers to the type of rank deficiency (see chapter on UDUC network equations) the new restriction eliminates.The corresponding restriction in the original CC-R or CC-S S-basis (Odijk et al. 2015)  and can be interpreted as the difference of BRSD/BSSD between different satellites/receivers.In the DD observation ΔΔ * s A s B AB,j , all satellite and receiver clocks and biases cancel out.
For the all-in-view network case of GNSS processing, the DD observation equations are given as, and eliminate all receiver and satellite clocks a n d b i a s e s , i .e ., F DD A rec A sat = 0 .We c a n extend the DD-matrix F DD to a full rank matrix F = F T DD F T DR T , and we transform the Design Matrix The same principles apply to the DD case as to the SD case: It can be shown that F DD is a basis matrix of N A rec A sat T .One can then solve the partial system of equations using only the DD observations and obtain an equal solution to the UDUC approach in which receiver and satellite clocks and biases are estimated as epoch-wise parameters when similar S-basis restrictions are chosen.Since the rank deficiencies 2b and 3b must be eliminated in the S-basis restrictions S ⊥ 1 T , ionospheric delays must be used as restrictions to eliminate rank deficiencies 2b and 3b instead of receiver/satellite biases, as is done in the CC models.Since biases have been assumed epoch-wise, the rank deficiencies between them and the ionospheric delays must also be eliminated at every epoch.An example S-basis is given in Table 1, denoted as CC-R-DD and CC-S-DD.The estimable ionospheric delays are then of double differenced nature.

On differencing with a constant bias assumption
As shown in the previous section, the classic SD and DD approaches are equivalent to the processing of UDUC observations when all biases and clocks are estimated as epochwise parameters.Since biases are usually estimated as constant daily parameters (Geng et al. 2019;Schaer et al. 2021) or as varying slowly in time, the functional model of SD and DD GNSS processing is expected to be weaker due to the greater number of implicit parameters.
The SD and DD approaches can be modified so that the assumption of constant biases is made, resulting in identical outcomes to undifferenced processing with constant biases.To (30) obtain a modification of the SD and DD approaches with the assumption of constant biases, one must look at the null space of the transposed partial design matrices under the assumption of constant biases.The design matrices are given with Note that the parameter order has been changed: The first column of both matrices contains the epoch-wise receiver/satellite clocks, and the second column the constant biases.The dimensions of the null spaces are given with The SD and DD observations ( 27) and ( 30) are also valid for the constant bias assumption, i.e., are linearly independent vectors in the null space of the respective matrix, but additional observations are required to achieve an equivalent reformulation under the assumption of constant biases.The number of additional observations increases with the number of epochs p , signals f , and receivers m/satellites n, respectively, as shown by ( 34) to (36) and will require between epoch and signal differencing in addition to the regular SD or DD.Given these considerations, the differencing approach, when assuming constant biases, does not seem to offer any significant benefits over an undifferenced approach.

Conclusions
We presented a method of reformulating the functional model of the UDUC observations to pre-eliminate parameters at an observational level.This reformulation addresses existing rank deficiencies and demonstrates how they propagate into the reformulated model.It was demonstrated that the IF and differencing approaches for processing GNSS observations are specific applications of this reformulation.By applying the reformulation methods to the UDUC observations of CDMA GNSS to obtain the various other approaches, we were able to identify the implicit assumptions made in the two IF-LC and the multiple differencing approaches presented here.All approaches, or more precisely the parameters and their VCM, are mathematically equivalent to the results obtained when the UDUC functional model is evaluated with the identified assumptions.Specifically, when using one IF code and phase observation, the dual frequency IF-LC approach assumes epochwise ionospheric slant delay parameters that are different for the code and phase observations.The latter assumption is discarded when using a third IF-LC, e.g., the HMW-LC, in addition to the IF code and phase observations.
The various differencing approaches make the assumptions that the eliminated receiver and or satellite clocks and biases are epoch-wise parameters in the functional model of the UDUC observations.With the general theory behind the reformulation of the original model, it would be possible to adapt the differencing approaches to make the implicit assumption of constant biases.
Understanding the implicit assumptions in the various methods of processing GNSS observations helps us to fairly compare the different approaches, their performance, and the estimated parameters.When combining satellite clocks and biases computed by different data centers into one consistent solution, it is essential to be aware of the underlying functional models and imposed S-basis restrictions since they affect the interpretation of those parameters.Furthermore, the presented reformulation can be applied to other methods of GNSS processing where LC of observations are used.One example is given with cycle slip detection in the preprocessing of GNSS observations.With the presented methods and ideas, implicit assumptions in various established cycle slip detection methods could be detected and applied to, e.g., UDUC GNSS observations.

1.
Between receiver and satellite parameters a. Between receiver and satellite clocks (size 1) b. Between receiver and satellite biases (size 2f ) 2. Between receiver biases and a. Receiver clocks (size n − 1) b.Ionospheric Slant Delays (size n − 1) 3. Between satellite biases and a. Satellite clocks (size m) b.Ionospheric Slant Delays (size m) 4. Between receiver biases and ambiguities (size f (n − 1)) 5. Between satellite biases and ambiguities (size fm).
a receiver-satellite link on f = l 2 different signals.Here, it is assumed that only observations of the same epoch are used.