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.
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:
| System | Country / Organization | Satellites | Orbital Planes | Altitude / Inclination | Core Frequencies |
|---|---|---|---|---|---|
| GPS | United States (USSF) | ~31 | 6 MEO planes | $20,200 \text{ km} / 55^\circ$ | L1 ($1575.42 \text{ MHz}$), L2 ($1227.60 \text{ MHz}$), L5 ($1176.45 \text{ MHz}$) |
| GLONASS | Russia (Roscosmos) | ~24 | 3 MEO planes | $19,100 \text{ km} / 64.8^\circ$ | G1 ($1602 \text{ MHz}$ FDMA/CDMA), G2 ($1246 \text{ MHz}$) |
| Galileo | European Union (ESA) | ~24-30 | 3 MEO planes | $23,222 \text{ km} / 56^\circ$ | E1 ($1575.42 \text{ MHz}$), E5a/b, E6 |
| BeiDou (BDS) | China (CNSA) | ~35 | GEO, 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:
- Space Segment: The constellation of orbiting satellites transmitting PRN codes, carrier signals, and navigation messages (ephemerides, clock bias parameters).
- Control Segment: Master Control Station (MCS), ground antennas, and tracking stations monitoring orbits, calculating ephemerides, and uploading clock/orbital corrections.
- 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:
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}$):
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 Source | Magnitude | Physical Mechanism | Mitigation 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 Bias | Up to $100 \text{ km}$ equivalent | Internal 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:
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 Mode | Occupation Time | Baseline Length | Accuracy Level | Common 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:
- Single Difference: Formed between two receivers tracking the same satellite. Cancels satellite clock error ($dT$).
- 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.
- 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:
Misclosure vector magnitude $W$:
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).
In GNSS carrier-phase positioning, what is the main purpose of Double Differencing between two receivers and two satellites?
A geodetic survey team operates 5 GNSS receivers simultaneously in a single static session. How many independent baselines are generated from this session?
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)?
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)?