Ray Tracing ofMagnetohydrodynamic Wavesin GeospaceA.D.M. WalkerAbstractThe method of ray tracing is reviewed for MHDwaves in stationary media and for media in the steady state.The method is generalized to allow for media that changearbitrarily but slowly in space and time. This requires anadditional equation representing the change in frequencyalong the ray as a result of Doppler shifts. Explicit expressionsfor ray tracing of transverse Alfvén waves and magnetosonicwaves are presented. Some illustrative examples arepresented.1. IntroductionFor more than fifty years, since its introduction byHaselgrove [1], ray tracing has been a powerful tool forstudying radio propagation in the magnetosphere [2]. Thishas not been the case for lower-frequencymagnetohydrodynamic waves, the chief reason being thatoften the wavelength is comparable to the size of themagnetosphere, and in no sense can the medium be regardedas slowly varying – a sine qua non for the validity of the raytracingapproximation. For shorter-period waves in the Pc3range (periods between 10 and 45 seconds), ray tracing maybe valid in the magnetosheath and solar wind, while in thesolar wind it may be valid for much longer periods. Anotherdifficulty arises in these regions: the medium is moving.Ray tracing in a moving medium has been discussed forsound waves by Lighthill [3], for example, and his discussionwas generalized for MHD waves by Walker [4, pp. 429-430.]. The treatment in this paper follows these treatmentsfor steady-state motions, but goes beyond them in allowingfor non-steady flow and the resulting Doppler shifts. We donot apply the methods to specific problems of ULF wavesin geospace, which will be the subject of a later paper, butinstead provide some illustrative numerical examples.As an example of the issues that are likely to arise insuch geospace applications, consider a transverse Alfvénwave. It is well known that the group velocity of such awave in a stationary medium is directed precisely along themagnetic field, and, since the medium is nondispersive, thegroup velocity is equal to the Alfvén velocity. The problemof ray tracing in such a case is trivial: the ray paths coincidewith the magnetic field lines. However, in a moving mediuma wave packet travels with a velocity that is the resultant ofthe Alfvén velocity and the velocity of the medium. Thesolar-wind velocity is generally much larger than the Alfvénvelocity, so that the ray paths are more nearly radial thanalong the magnetic-field direction. The case of themagnetosonic waves is even more complicated.In this paper, we first discuss the geometrical optics ofMHD waves in uniform media. We then review thederivation of the ray-tracing equations in moving media inthe steady state. Finally, we show how these equations canbe generalized to include arbitrary flows, provided changesin the background flow with time are slow compared to theperiod of the wave. We present explicit equations to allowthe numerical computations of MHD rays, and performsome computations to illustrate the techniques.2. Geometrical Optics of MHDWaves in Uniform Media2.1 Stationary MediaIn a uniform stationary compressible plasma, atfrequencies much less than the lowest ion gyrofrequency,there are three characteristic magnetohydrodynamic wavesthat can be propagated. These can be identified as thetransverse Alfvén wave, with the dispersion relationA. D. M. Walker is with the School of Physics, UniversityA. D. M. Walker is with the School of Physics, Universityof KwaZulu-Natal, Durban 4000, South Africa;e-mail: walker@ukzn.ac.za.24The<strong>Radio</strong> <strong>Science</strong> <strong>Bulletin</strong> No <strong>325</strong> (<strong>June</strong> <strong>2008</strong>)
orω =± k•V (1)A2 2 2 2ω = kV A cos θ , (2 )and the fast and slow magnetosonic waves, with thedispersion relationor( ) ( ) 2A S S Aω 4 − k 2 V 2 + V 2 + k 2 V2 k• V = 0 (3)( )4 2 2 2 2 4 2 2 2k VA VS k VAVSω − + ω + cos θ = 0 , (4)where V A is the Alfvén speed, V S is the sound speed, andθ is the angle between the direction of the magnetic fieldand the wave vector, k [4, Equations (7.5) and (7.6)].These dispersion relations have a special property.They show that the phase velocity, VP≡ ω k , is a functionof θ but not of frequency, ω . The propagation of the wavesis therefore anisotropic (the phase speed depends ondirection) but nondispersive (the phase speed does notdepend on frequency). Anisotropic waves have a groupvelocity (and hence a direction of energy propagation) thatin general is not parallel to the wave vector, which is normalto the wavefronts.As described, for example, by Walker [2], propagationof waves in uniform media can be understood in terms of theproperties of three characteristic surfaces: (i) The wavevectorsurface, given by the magnitude of the wave vectoras a function of its direction as represented by the polarangles, θ and φ ; (ii) the group-velocity surface, given bythe magnitude of the group velocity as a function of itsdirection; (iii) the ray surface, given by the magnitude of theray velocity as a function of the direction of energypropagation. The ray velocity is defined as having amagnitude equal to the speed at which wavefronts movealong the direction of the group velocity. In the special caseof MHD waves, because they are nondispersive, thecomponent of the group velocity in the direction of the wavenormal is equal to the phase velocity, and the ray- andgroup-velocity surfaces are identical.In radio propagation theory, it is useful to normalizethe k vector in the form of a refractive index vector,n = ck ω . In MHD, the phase speed, ω k , is severalorders of magnitude less than the speed of light, c, so that3 4this leads to refractive indices of the order of 10 or 10 .It is better to express the normalization in terms of adifferent characteristic speed. The Alfvén speed,V ≡ B µρ, (5)A0is suitable. We define a “refractive-index” vector, ornormalized wave vector,n ≡ VAk ω . (6)The dispersion relations of Equations (2 ) and (4) thenbecomencosθ =± 1 , (7)( S)4 2 2 2 2SnU cos θ − n 1+ U + 1 = 0 ,where US ≡ VS VA. The refractive index is a function ofthe polar angle, θ , and independent of the azimuthal angle,φ . The refractive-index surface in a stationary medium istherefore a surface of rotation about the magnetic-fielddirection. In particular, the surface for the Alfvén wave,represented by the first equation, is simply a pair of planeswhere the component of the refractive-index vector parallelto the field is unity. The refractive-index surfaces for thethree characteristic waves are shown in Figure 1. The fastwave has the smallest refractive index of the three waves forany direction of propagation. It is a closed surface with anovid shape. The slow wave has the largest refractive index.Propagation is not possible for directions approximatelyperpendicular to the magnetic-field direction, so there aretwo separate surfaces of rotation that are asymptotic to thepair of cones defined by−1cosθres U S V A V S=± ≡ . (8)The two planes representing the transverse Alfvén wave liebetween these two surfaces, with n || = 1 .The group velocity is given by Equation (2),V G =∇ k ω , (9)and thus it is normal to the refractive-index surfaces, whichare surfaces of constant ω in k space. We immediately seethat for the transverse Alfvén wave, the group velocity, andhence the direction of energy propagation, is exactly alongthe magnetic-field direction. For the fast wave, although itis not parallel to the wave normal, it can make any anglewith the magnetic field, so that energy can be propagated inany direction through the medium of the fast wave. For theslow wave, the direction of energy propagation generallymakes a small angle with the magnetic field, so that energyis always propagated approximately in the direction of thefield and cannot be propagated transverse to it.The ray surface for the transverse Alfvén wave isdegenerate. No matter what the direction of the wavenormal, the ray velocity is exactly equal in magnitude to theAlfvén speed, and is directed along the magnetic field. Theray- and group-velocity surfaces are therefore a pair ofThe<strong>Radio</strong> <strong>Science</strong> <strong>Bulletin</strong> No <strong>325</strong> (<strong>June</strong> <strong>2008</strong>) 25
- 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 20 and 21: Haselgrove equations [2]. Furthermo
- Page 22 and 23: anddQ 1 r ∂n= + n −Qdθ2 2 2n
- 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