Haselgrove equations [2]. Furthermore, recastingEquation (16) in terms of q [2], we obtain⎡ 2( ) ( ) ( )2 2 2q 1− X − Y + q Y p 2 1 X Y⎤⎢• + − −⎥+ 1− X = 0⎣⎦(23)from which q can be derived (note that there are twopossible values of q, corresponding to two of the possiblesolutions to Equation (11)). The corresponding current,and hence the electric field through Equation (5), can beobtained from the homogeneous Equations (11) using thecalculated values of p on the ray. However, Equations (11)do not determine the magnitude of the field. This must becalculated by tracking the power out from the source. Thiscan be done by tracing out ray tubes, and noting that within20 2 0Ap E η , must remainconstant (A is the cross-sectional area of the tube). Themagnitude of the electric field can then be calculated fromthe power.these tubes, the power, ( )There remains the issue of solving the above raytracingequations. Starting with some initial value, w () 0 ,for w , the system of ordinary differential equations,dw dt = W( w), can be solved using Runge-Kuttatechniques. This was suggested in the Haselgrove andHaselgrove paper [2]. If we start a ray at the Earth’s surface,the initial values of p are given by the direction cosines atthis starting point. However, the extreme variations that canexist in the ionosphere often make it necessary to vary theincrement in t at each stage of the solution process, in orderto maintain accuracy. Haselgrove and Haselgrove [2]suggested the introduction of a scaling to produce thiseffect, but the Runge-Kutta-Fehlberg (RKF) method (see[8] for example) can also achieve this.Consider the solution at time t and prospectiveincrement ∆ t . Form the vectors k1, k2, k3, k4, k5, andk through the process6k1( )=∆ tW w , (24)( )k =∆ tW w+ k , (25)2 1 4(k =∆ tW w+ 439k 216 − 8k + 3680k5135 1 2 3(4)− 845k4104 , (28)k =∆tW w− 8k 27 + 2k −3544k25656 1 2 3+ 1859k4104 − 11k40 . (29)4 5The approximate solution, wt () +∆ w at t+∆ t , is obtainedusing∆ w = 25k 216 + 1408k25654 1 3+ 2197k4104 − k 5(29)4 5for the fourth-order Runge-Kutta method, and by∆ w = 16k 135 + 6656k128255 1 3+ 28561k 56430 − 9k 50 + 2k55 (30)4 5 6for the fifth-order method. The global truncation error will45be O( ∆ t ) for the fourth-order method and O( ∆ t ) forthe fifth-order method. The magnitude of the differencebetween the fourth- ( w 4 ) and fifth- ( w 5 ) order estimates,∆ w= w5 − w4, gives an estimate of the global truncationerror in the fourth-order method.The Runge-Kutta-Fehlberg method proceeds byadjusting the step at each stage so that this global errorremains close to a pre assigned value, ∆ E . That is, at eachstage we adjust ∆ t so that1⎛∆E⎞4∆ tnew=∆told⎜ ⎟⎝∆w⎠). (31)( 3 32 9 32)k =∆ tW w+ k + k , (26)3 1 2(k =∆ tW w+ 1932k 2197 −7200k21974 1 23)+ 7296k2197 , (27)This approach has been used in [4] to provide an efficientnumerical implementation of the Haselgrove equations.3. The Weak-Background-Field LimitFor the collision-free case ( U = 1), Equation (3) canbe inverted to yield20The<strong>Radio</strong> <strong>Science</strong> <strong>Bulletin</strong> No <strong>325</strong> (<strong>June</strong> <strong>2008</strong>)
− jωε0XJ = E− jY× E− Y•EY21−Y( ). (32)Up until now, we have assumed that both the plasmafrequency, ω p , and the gyro frequency, ω H = eB0m ,are of the same order as the wave frequency, ω . We willnow assume that ω >> ωpand, hence, we can ignore the2term of second order, Y , in Equation (32). Substitutingthe trial solution E = E 0 exp ( − jβϕ ) into Equation (1),2the β and β terms yieldand( ) ( )−∇ϕ∇ ϕ• E + ∇ ϕ•∇ϕE − 1− X E = 0 (33)0 0 02( ) 2−∇ϕ∇• E − ∇ ∇ ϕ• E + ∇ ϕ•∇ E + ∇ ϕE0 0 0 0= β XY× E , (34)0respectively. By taking the divergence of Equation (1), we∇• E = j ωε ∇• J , from which and ∇ ϕ • E0 = 0obtain ( )0( 1− X) ∇• E =∇ X • E −βX∇ ϕ• ( Y×E )0 0 0for the two leading orders. Consequently, Equations (33)and (34) reduce to⎡ 1 dx dx= β X ⎢Y× E − • Y×E⎣ 1−X dt dt( )0 0⎤ .⎥⎦(37)Taking the dot product of Equation (37) with E 0 , we find2 2 2that dE0 dt+∇ ϕ E0 = 0 . Combining this withEquation (37), we obtain an equation for the polarizationvector, P = E0 E0, of the formdP dx 22 + ∇ lnn• Pdt dt⎡ 1 dx dx ⎤= β X ⎢Y× P− •( Y×P)1−X dt dt⎥⎣⎦.(38)Essentially, the background magnetic field manifests itselfas a rotation of the polarization vector about the ray direction.4. The Negligible-Background-Field LimitWhen the magnetic field of the Earth can be ignored( B 0 = 0 ), the Haselgrove equations can be reduced tod ⎛ dx⎞ n =∇ n , (39)⎜ ⎟ds ⎝ ds ⎠and( ϕ ϕ X) E 0∇ •∇ − 1+ = 0 (35)−∇ϕβX∇ X • E + ∇ ϕ ∇ ϕ • Y × E1−X 1−X( )0 0where s is the distance along the ray path. The equationsare themselves the Euler-Lagrange (EL) equations for thevariational principleR∫T() 0δ n x ds=, (40)20 0 0+∇ 2 ϕ•∇ E +∇ ϕE = βXY× E (36)From Equation (35), we obtain the ray-tracing equationsdp dt n 2 2=∇ 2 and dx dt = p , where dt = dϕp (i.e.,2t is the group distance) and n = 1− X . These equations donot contain the background magnetic field, and hence canbe solved far more efficiently than the original Haselgroveequations. Once again, the Runge-Kutta-Fehlberg methodcan be used with great effect.Equation (36) can now be recast as a differentialequation for E0along a ray:dE0dx2 + ∇ lnn • E +∇ Edt dt2 20 ϕ 0i.e., Fermat’s principle. If we now consider the case for a raypath that remains in a plane, we can use a polar coordinatesystem ( r,θ ), and the variational Euler-Lagrange principlebecomes12R ⎡ 2⎛ dr ⎞ ⎤δ n( r, θ)⎢ + 1 ⎥ dθ= 0 . (41)∫ ⎜ ⎟dθT⎢⎣⎝⎠ ⎥⎦The above ray-tracing equations can be used when anydeviations from the great-circle path (certainly the case fora uniform ionosphere) are negligible. Coordinate r is thedistance from the center of the Earth, and θ is an angularcoordinate along the great-circle path. The Euler-Lagrangeequations for the above variational principle can be reduced[9] toThe<strong>Radio</strong> <strong>Science</strong> <strong>Bulletin</strong> No <strong>325</strong> (<strong>June</strong> <strong>2008</strong>) 21
- Page 4 and 5: URSI Forum on Radio Scienceand Tele
- Page 6 and 7: URSI Accounts 2007Income in 2007 wa
- Page 8 and 9: EURO EUROA2) Routine Meetings 5,372
- Page 10 and 11: URSI Awards 2008The URSI Board of O
- Page 12 and 13: • Compatibility refers to the abi
- Page 14 and 15: Introduction to Special Sections Ho
- Page 16 and 17: In the December 2008 issue, Part 2
- Page 18 and 19: mwÆ = eE + ew× B −mνw , (2)whe
- Page 22 and 23: anddQ 1 r ∂n= + n −Qdθ2 2 2n
- Page 24 and 25: Ray Tracing ofMagnetohydrodynamic W
- Page 26 and 27: Figure 1a. The refractive-index sur
- Page 28 and 29: 2 ⎡ ∂VAj, ∂VS⎤ω ⎢VAj, +
- Page 30 and 31: ( ω kV j j)dx ∂i 0 + ∂ω0,= =
- Page 32 and 33: compared with the wavelength, i.e.,
- Page 34 and 35: decreasing to half their value in a
- Page 36 and 37: Practical Applicationsof Haselgrove
- Page 38 and 39: **.modeled TID ionospheres using Ha
- Page 40 and 41: ∂ϕ∂ϕ1 kϕ∂r=− +, (12)∂
- Page 42 and 43: screen (which is implemented in a s
- Page 44 and 45: directed boresight and, similar to
- Page 46 and 47: Typical OTH radar ionospheric sound
- Page 48 and 49: ay-homing calculations for analyzin
- Page 50 and 51: Figure 1. A Cartesian geometry cons
- Page 52 and 53: orc dϕh′ = , (5a)4πdfincreases
- Page 54 and 55: in four runs by five rays, giving a
- Page 56 and 57: GPS : A Powerfull Tool forTime Tran
- Page 58 and 59: altitude of ~26562 km will tick mor
- Page 60 and 61: Each component of the errors is ass
- Page 62 and 63: Figure 3. The effect of GPS timeerr
- Page 64 and 65: Figure 6. The frequencystability of
- Page 66 and 67: Figure 9. The participating laborat
- Page 68 and 69: unlike the possibility of a gap in
- Page 70 and 71:
Table 5. The capabilities of GPS ti
- Page 72 and 73:
12. References1. C. Audoin and J. V
- Page 74 and 75:
ConferencesCONFERENCE REPORT12TH IN
- Page 76 and 77:
URSI CONFERENCE CALENDARURSI cannot
- Page 78 and 79:
The Journal of Atmospheric and Sola
- Page 80:
APPLICATION FOR AN URSI RADIOSCIENT