9.3 Global Navigation Satellite Systems (GNSS/GPS), Carrier-Phase Positioning, RTK, and Static Baselines

Key Takeaways

  • GNSS constellations (GPS, GLONASS, Galileo, BeiDou) operate across core frequencies (L1, L2, L5) with three functional segments: Space, Control, and User.
  • Carrier-phase positioning achieves millimeter-to-centimeter accuracy by resolving integer phase ambiguity (N), compared to code-based pseudorange positioning which yields meter-level accuracy.
  • Ionospheric delay is dispersive and eliminated via dual-frequency combination (L_IF), whereas tropospheric delay is non-dispersive and mitigated through zenith delay modeling.
  • Dilution of Precision (DOP) quantifies geometric satellite constellation strength, where lower values (PDOP < 3) reflect optimal geometry and higher precision.
  • High-precision survey modes include Static GNSS (primary control), Real-Time Kinematic (RTK/NTRIP for real-time boundaries), Precise Point Positioning (PPP), and rigorous baseline loop closure analysis.
Last updated: July 2026

9.3 Global Navigation Satellite Systems (GNSS/GPS), Carrier-Phase Positioning, RTK, and Static Baselines

Global Navigation Satellite Systems (GNSS) have revolutionized geodetic surveying, replacing traditional optical triangulation networks with satellite-based baseline vectors referencing global geodetic datums (such as WGS84 and ITRF/PRS92).


1. GNSS Constellations & System Segments

Modern multi-constellation GNSS integrates signals from four major global satellite networks:

SystemCountry / OrganizationSatellitesOrbital PlanesAltitude / InclinationCore Frequencies
GPSUnited States (USSF)~316 MEO planes$20,200 \text{ km} / 55^\circ$L1 ($1575.42 \text{ MHz}$), L2 ($1227.60 \text{ MHz}$), L5 ($1176.45 \text{ MHz}$)
GLONASSRussia (Roscosmos)~243 MEO planes$19,100 \text{ km} / 64.8^\circ$G1 ($1602 \text{ MHz}$ FDMA/CDMA), G2 ($1246 \text{ MHz}$)
GalileoEuropean Union (ESA)~24-303 MEO planes$23,222 \text{ km} / 56^\circ$E1 ($1575.42 \text{ MHz}$), E5a/b, E6
BeiDou (BDS)China (CNSA)~35GEO, IGSO, MEO$21,150 \text{ km} \text{ (MEO)} / 55^\circ$B1 ($1561.098 \text{ MHz}$), B2, B3

Functional Segments

Every GNSS operates via three interconnected segments:

  1. Space Segment: The constellation of orbiting satellites transmitting PRN codes, carrier signals, and navigation messages (ephemerides, clock bias parameters).
  2. Control Segment: Master Control Station (MCS), ground antennas, and tracking stations monitoring orbits, calculating ephemerides, and uploading clock/orbital corrections.
  3. User Segment: GNSS antennas, field survey receivers, and post-processing software used by geodetic engineers.

2. Pseudorange vs. Carrier-Phase Measurements

GNSS receivers determine distance to satellites using two fundamental measurement types:

Pseudorange (Code) Measurement

Pseudorange ($P$) measures time of arrival ($\Delta t$) of the satellite's Pseudo-Random Noise (PRN) code:

P=ρ+c(dtdT)+I+T+ϵPP = \rho + c \cdot (dt - dT) + I + T + \epsilon_P

where $\rho$ is true geometric range, $c$ is speed of light, $dt$ is receiver clock error, $dT$ is satellite clock error, $I$ is ionospheric delay, $T$ is tropospheric delay, and $\epsilon_P$ is code noise.

  • Accuracy: Code-based positioning (standalone GNSS) yields accuracies of $0.5 \text{ m}$ to $3.0 \text{ m}$.

Carrier-Phase Measurement

Carrier-phase ($\Phi$) measures the phase of the continuous sinusoidal carrier wave (e.g. L1 wavelength $\lambda_1 \approx 19.05 \text{ cm}$):

Φ=1λρ+N+cλ(dtdT)Iλ+Tλ+ϵΦ\Phi = \frac{1}{\lambda} \rho + N + \frac{c}{\lambda}(dt - dT) - \frac{I}{\lambda} + \frac{T}{\lambda} + \epsilon_\Phi

where $N$ is the Integer Phase Ambiguity—the unknown whole number of carrier cycles between satellite and receiver at signal lock-on.

  • Accuracy: Resolving $N$ to correct integers (Fixed Solution) yields millimeter-to-centimeter precision ($1 \text{ mm} - 10 \text{ mm}$).

3. Atmospheric Delays and Error Sources

Error SourceMagnitudePhysical MechanismMitigation Strategy
Ionospheric Refraction ($I$)$1 \text{ m} - 50 \text{ m}$Dispersive delay caused by free electrons in ionosphere ($I \propto 1/f^2$).Dual-Frequency Combination ($L_{\text{IF}}$) cancels 99% of 1st-order delay. Differential processing.
Tropospheric Refraction ($T$)$2 \text{ m} - 25 \text{ m}$Non-dispersive refraction in lower atmosphere (90% hydrostatic/dry, 10% wet).Standard models (Saastamoinen, Hopfield) + zenith troposphere parameter estimation.
Multipath Effect$0.1 \text{ m} - 5 \text{ m}$Signal reflections from buildings, ground, or structures arriving out of phase.Choke ring antennas, elevation masks ($15^\circ$), avoiding reflective surroundings.
Satellite Clock & Orbit$1 \text{ m} - 3 \text{ m}$Orbit ephemeris drift and satellite clock instability.Differential double differencing, precise IGS orbits.
Receiver Clock BiasUp to $100 \text{ km}$ equivalentInternal quartz crystal clock offset.Eliminated by double-differencing or estimating as receiver parameter.

4. Dilution of Precision (DOP)

Dilution of Precision (DOP) is a dimensionless scalar factor describing the impact of satellite constellation geometry on positioning accuracy:

Positioning Error σpos=DOP×σrange\text{Positioning Error } \sigma_{\text{pos}} = \text{DOP} \times \sigma_{\text{range}}

Types of DOP

  • GDOP (Geometric DOP): Overall geometry combining 3D position and time ($\text{GDOP} = \sqrt{\text{PDOP}^2 + \text{TDOP}^2}$).
  • PDOP (Position DOP): 3D spatial geometry ($\text{PDOP} = \sqrt{\text{HDOP}^2 + \text{VDOP}^2}$).
  • HDOP (Horizontal DOP): 2D horizontal positional strength.
  • VDOP (Vertical DOP): Vertical height determination strength (typically $1.5\times$ to $2\times$ higher than HDOP due to satellites only being visible above horizon).

Rule of Thumb: Ideal surveys require $\text{PDOP} < 3$. If $\text{PDOP} > 6$, satellite geometry is weak, and precision survey observations must be suspended.


5. GNSS Survey Positioning Modes

Positioning ModeOccupation TimeBaseline LengthAccuracy LevelCommon Geodetic Applications
Static GNSS$30 \text{ min} - 4+ \text{ hours}$Long ($> 20 \text{ km}$)$\pm (3 \text{ mm} + 0.5 \text{ ppm})$Primary geodetic control networks (PRS92/PGP), tectonic monitoring.
Rapid Static$5 \text{ min} - 20 \text{ min}$Short ($< 10-15 \text{ km}$)$\pm (5 \text{ mm} + 1 \text{ ppm})$Secondary control, DENR boundary surveys.
Real-Time Kinematic (RTK)Instantaneous ($1 - 5 \text{ sec}$)Short ($< 10 \text{ km}$)$\pm (10 \text{ mm} + 1 \text{ ppm})$Boundary stakeout, topographic surveys, engineering construction.
Network RTK (VRS/MAC)Instantaneous ($1 - 5 \text{ sec}$)Extended ($< 50 \text{ km}$)$\pm (10 \text{ mm} + 0.5 \text{ ppm})$Continuous CORS networks (e.g. NAMRIA Pag-ASA CORS).
Precise Point Positioning (PPP)$15 \text{ min} - 1 \text{ hour}$ (convergence)Unlimited (Global)$\pm (1 \text{ cm} - 5 \text{ cm})$Remote regional mapping, aerial survey ground control.

6. Baseline Processing, Differencing, and Loop Closure

Differencing Techniques

In differential GNSS processing, carrier-phase observations are differenced to cancel systematic errors:

  1. Single Difference: Formed between two receivers tracking the same satellite. Cancels satellite clock error ($dT$).
  2. Double Difference: Formed by differencing two single-differences between two receivers and two satellites. Cancels both receiver and satellite clock biases, leaving integer phase ambiguities ($N$) as solveable integers.
  3. Triple Difference: Formed by differencing double-differences between consecutive time epochs. Cancels integer phase ambiguities; used primarily for cycle slip detection.

Independent vs. Dependent Baselines

For a session with $n$ GNSS receivers observing simultaneously:

  • Total Baselines Formed: $N_{\text{total}} = \frac{n(n - 1)}{2}$
  • Independent (Non-Trivial) Baselines: $N_{\text{ind}} = n - 1$

Critical Exam Rule: Processing software must only include independent baselines ($n-1$) in network adjustments. Including trivial (dependent) baselines artificially inflates statistical degrees of freedom.

Loop Closure Analysis

A closed loop formed by GNSS baseline vectors must sum to zero:

ΔX=0,ΔY=0,ΔZ=0\sum \Delta X = 0, \quad \sum \Delta Y = 0, \quad \sum \Delta Z = 0

Misclosure vector magnitude $W$:

W=(ΔX)2+(ΔY)2+(ΔZ)2W = \sqrt{(\sum \Delta X)^2 + (\sum \Delta Y)^2 + (\sum \Delta Z)^2}

According to DENR standards (DAO 2007-29), primary control GNSS baseline loops must satisfy relative error tolerances $\le 1 : 100,000$ or maximum misclosure limit $W \le 10 \text{ mm} + 2 \text{ ppm} \times L$ (where $L$ is total loop length).

Test Your Knowledge

In GNSS carrier-phase positioning, what is the main purpose of Double Differencing between two receivers and two satellites?

A
B
C
D
Test Your Knowledge

A geodetic survey team operates 5 GNSS receivers simultaneously in a single static session. How many independent baselines are generated from this session?

A
B
C
D
Test Your Knowledge

If the Horizontal Dilution of Precision is HDOP = 1.2 and the Vertical Dilution of Precision is VDOP = 1.6, what is the Position Dilution of Precision (PDOP)?

A
B
C
D
Test Your Knowledge

Which GNSS positioning mode uses a single receiver without a local base station by applying precise satellite orbit and clock products provided by global ground networks (such as IGS)?

A
B
C
D