y the speed of light c , yields the range between the satellite and the receiver antenna. Since therange has several errors and biases lumped into it, it is called pseudorange (e.g. Langley, 1993). Thepseudorange observati<strong>on</strong> equati<strong>on</strong> can be expressed as (Komjathy, 1997):Heresi, L1( ) i<strong>on</strong> L trop rjL Pε1 P1P = ρ + c ⋅ dT − dt + d + d + b + b + mp + , (3.1)1 , 1 , 1si, L2( ) i<strong>on</strong> L trop rjL P2 P2P = ρ + c ⋅ dT − dt + γ ⋅ d + d + b + b + mp + ε . (3.2)2 , 1 , 2P1 , P2are pseudorange measurements <strong>on</strong> L1 and L2 in units of distance,ρ is the geometric range between the satellite and the receiver,dT,dt the receiver and satellite clock offsets from GPS time,di<strong>on</strong>, L1is the i<strong>on</strong>ospheric delay <strong>on</strong> L1,2 2⎛ f ⎞ ⎛1L ⎞1γ = ⎜ ⎟ = ⎜ ⎟ is the frequency ratio squared,f L⎝ 2 ⎠ ⎝ 2 ⎠dtropis tropospheric delay,si, L1 si, L2b , b are the satellite instrumental delays <strong>on</strong>1rj, L1 rj, L245P and P2respectively for satellite si,b , b are the receiver instrumental delays <strong>on</strong> P1and P2respectively for satellite rj,mp , mp are the multipaths <strong>on</strong> P1and P2measurements, andP1 P2ε , ε are the receiver noise <strong>on</strong> P1and P2respectively.P1 P2By subtracting equati<strong>on</strong> (3.2) from equati<strong>on</strong> (3.1) and rearranging the terms we obtain:where( γ )P − P = − ⋅ d + ⎡b − b + b − b ⎤ +⎣⎦si , L1 si, L21 21i<strong>on</strong>, L1 rj, L1 rj, L2νpp P1 P P2 1 P2, (3.3)ν = mp − mp + ε − ε .(3.4)It is crucially important to note that the i<strong>on</strong>ospheric term in equati<strong>on</strong> (3.3) cannot be measureddirectly from pseudorange measurements due to the term in the square brackets <strong>on</strong> the right hand
side of the equati<strong>on</strong>, which is slowly varying in time. However, since the i<strong>on</strong>ospheric term cannot bemeasured directly, it needs to be estimated al<strong>on</strong>g with the term in the square brackets in equati<strong>on</strong>(3.3).The first order equati<strong>on</strong> of the Applet<strong>on</strong>-Hartree formula (2.12) also gives an expressi<strong>on</strong> for thei<strong>on</strong>ospheric delay in terms of TEC as:di<strong>on</strong>TEC= d = 40.3 ⋅ .(3.5)i<strong>on</strong>, L1 2f1This equati<strong>on</strong> gives centimeter-level accuracy for i<strong>on</strong>ospheric delay (Komjathy, 1997). The equati<strong>on</strong>(3.5) neglects terms higher than the 2 nd order in the Applet<strong>on</strong>-Hartree formula as described insecti<strong>on</strong> 2.6.1.After substituting equati<strong>on</strong> (3.5) into equati<strong>on</strong> (3.3) and doing some algebra, the expressi<strong>on</strong> for theTEC using pseudorange (P) observati<strong>on</strong>s in TECU is given by:( )TEC = 9.52 P − P . (3.6)P2 1Since electromagnetic waves such as GNSS signals also experience phase advances when traversingthe i<strong>on</strong>osphere, the observati<strong>on</strong> equati<strong>on</strong>s for carrier phase ( Φ ) measurements <strong>on</strong> L1 and L2 are:where( )Φ = + c ⋅ dT − dt + N − d + d + b + b + mp + (3.7)ϕ , si, L11ρ λ ,1 1 i<strong>on</strong>, L1 trop ϕ , rj, L1 ϕ1 εϕ1 ( )Φ = ρ + c ⋅ dT − dt + λ N − γ . d + d + b + b + mp + ε ,(3.8)ϕ , si, L22 2 2 i<strong>on</strong>, L1 trop ϕ , rj, L2 ϕ 2 ϕ 2Φ1,Φ2are the carrier phase observati<strong>on</strong>s <strong>on</strong> L1 and L2 respectively in distance units,λ1 , λ2are the wavelengths of the L1 and L2 carriers respectively,N , N are the unknown L1 and L2 integer carrier phase ambiguities,b1 2, bϕ , si, L1 ϕ , si, L2, rj, L1 , rj, L2are the satellite instrumental delays <strong>on</strong> L1 and L2 carrier phase for satellite si,bϕ , bϕ are the receiver instrumental delays <strong>on</strong> L1 and L2 carrier phase for receiver rj,46
- Page 1 and 2: SOLAR CYCLE EFFECTS ON GNSS-DERIVED
- Page 3 and 4: ACKNOWLEDGEMENTSI would like to ack
- Page 5 and 6: SOHOSXITECVLBIVTECVTMUNB-IMTULMCARU
- Page 7 and 8: 2.6.3 Ionospheric group delay and p
- Page 9 and 10: 7.2 Recommendation for future work
- Page 11 and 12: Figure 4.4: The scatter plot of mid
- Page 13 and 14: Figure 6.5: Panel (a) depicts the c
- Page 15 and 16: Chapter 1INTRODUCTION1.1 Background
- Page 17 and 18: station using stochastic parameters
- Page 19 and 20: (GTEC) computed using UNB-IMT with
- Page 21 and 22: Chapter 2THE SUN, EARTH`S IONOSPHER
- Page 23 and 24: that the solar mat
- Page 25 and 26: charged ion is created, and the neg
- Page 27 and 28: 2.3.1 The D layerThe D layer is the
- Page 29 and 30: Figure 2.4: Major geographic region
- Page 31 and 32: latitude of 78 to 80 near local n
- Page 33 and 34: with different propagation velociti
- Page 35 and 36: 2.6 Ionospheric effects</st
- Page 37 and 38: where ω (radian.cωCZ = ,(2.5)ω1s
- Page 39 and 40: nwhich is the group refractive inde
- Page 41 and 42: (e.g. McKinnell, 2002). The ionogra
- Page 43 and 44: GPS is a satellite-based system for
- Page 45 and 46: solar flux between
- Page 47 and 48: Chapter 3THE IONOSPHERIC TOTAL ELEC
- Page 49: although efforts are continuously i
- Page 53 and 54: TECcombi= TEC −ϕin2∑nj=−2p
- Page 55 and 56: this reference frame compared to th
- Page 57 and 58: 1997). For more information about t
- Page 59 and 60: As first validation, Figure 3.2 ill
- Page 61 and 62: Chapter 4VALIDATION OF UNB-IMT USIN
- Page 63 and 64: esults using ITEC measurements as o
- Page 65 and 66: code described in Chapter 3 to comp
- Page 67 and 68: Figure 4.2: Comparison of daily mid
- Page 69 and 70: (b) Main phase: the TEC is expected
- Page 71 and 72: (Fedrizzi, et al., 2005). Further i
- Page 73 and 74: Figure 4.5: Computed difference ∆
- Page 75 and 76: The seasonal variation is difficult
- Page 77 and 78: Figure 4.7 (a)-(c) depicts the midd
- Page 79 and 80: 5. Summary and ConclusionsThis work
- Page 81 and 82: Chapter 5MAPPING GNSS-DERIVED TEC O
- Page 83 and 84: A detailed description of CELIAS/SE
- Page 85 and 86: 30 ° W 0 ° 30 ° E 60 ° E 90 °
- Page 87 and 88: latitudes is not associated with ge
- Page 89 and 90: However, a further investigation wa
- Page 91 and 92: Figure 5.6: The day 301, 2003 <stro
- Page 93 and 94: who reported that the largest TEC e
- Page 95 and 96: Figures 5.7 and 5.8 depict examples
- Page 97 and 98: dramatically to the X - ray and Ext
- Page 99 and 100: the X9 flare, (e) shows the TEC map
- Page 101 and 102:
solar cycl
- Page 103 and 104:
chapter demonstrated the capabiliti
- Page 105 and 106:
(Hobiger et al., 2006). The main pu
- Page 107 and 108:
(d) Based on (a), (b) and (c), VTEC
- Page 109 and 110:
Figure 6.3: Comparison of VTEC (dar
- Page 111 and 112:
However, for days in which data wer
- Page 113 and 114:
1994) and the follow-up to CONT95 (
- Page 115 and 116:
6.4.2 CONT05 CampaignCONT05 was a t
- Page 117 and 118:
To further investigate the capabili
- Page 119 and 120:
space geodetic techniques during th
- Page 121 and 122:
(2) Short-term variations of TEC (p
- Page 123 and 124:
conducted for the year 2002 (near <
- Page 125 and 126:
REFERENCESBartels, J., Heck, N. H.,
- Page 127 and 128:
Precision Geodesy using the Mark-II
- Page 129 and 130:
Fuller-Rowell, T. J., Codrescu, M.
- Page 131 and 132:
Jakowski, N., Heise, S., Wehrenpfen
- Page 133 and 134:
Langley, R. B., The GPS observables
- Page 135 and 136:
McNamara, L. F. Radio amateur guide
- Page 137 and 138:
Prölss, G. W. On explaining the lo
- Page 139 and 140:
Stubbe, P. The ionosphere as a Plas
- Page 141:
Zhang, D. H., Xiao, Z. Study of ion