Rayleigh Scattering: Phase Function

Definitions of Angles

See Figure 11. The angle between the incident flux and the local normal is given by alpha and µ0 ident cos(alpha). The angle between the local normal and the observer is given by epsilon and µ ident cos(epsilon). The phase angle, theta, is the angle between the incident and emission through the origin. The azimuthal angle, Deltapsi, is the projection of the solar incidence and emission directions onto the local horizon. Using spherical trigonometry cosine laws (e.g., see CRC pg. 146) we can easily obtain the value of Deltapsi) given µ0, µ, and theta. For solar scattering to an observer above the atmosphere the angles are related by:

cos(theta) = sqrt[(1-µ02)(1-µ2)] · cos(Deltapsi) -/+ µ µ0   (S = +, T = -)

This is the method employed by Chandrasekar, Liou, and others. We could also write the equations in terms of the scattering angles. The scattering angle between incidence and emission, Theta, which is related to the phase function, is given by Theta = pi - theta. The scattering angle in the horizontal plane, Deltaphi is related to the azimuthal angle, and is given by Deltaphi= pi - Deltapsi. This definition is used by Hansen, Tomasko, and Danielson. It is the method employed within the VIAMP programs. Note the sign change between the two methods.

cos(Theta) = sqrt[ (1-µ02)(1-µ2)] · cos(Deltaphi) ± µ µ0   (T = +, S = -)

The Rayleigh phase function is given in terms of scattering angles or phase angles by

Pray(Theta) = (3/4) (1 + cos2(Theta)) = (3/4)(1 + cos2(theta)) = Pray(theta)

The cos2 term can be expressed in terms of µ, µ0, and Deltaphi as

cos2(Theta) = (1-µ02)(1-µ2) · cos2(Deltaphi) ± 2µµ0 sqrt[ (1-µ02)(1-µ2) ] · cos(Deltaphi) + (µ2 µ02)

cos(2Deltaphi) = 2 cos2(phi) - 1   or,

cos2(phi) = 1/2 cos(2Deltaphi) + 1/2

cos2(Theta) = 1/2 (1-µ02)(1-µ2) · cos(2Deltaphi) ± 2 µ µ0 sqrt[ (1-µ02)(1-µ2) ] · cos(Deltaphi) + µ2 µ02 + 1/2 (1-µ02)(1-µ2)

So the Rayleigh phase function can be written in terms of 3 Fourier moments:

Pray(Theta) = P1 + P2 cos(Deltaphi) + P3 cos(2Deltaphi)

P0 = 3/4 [ 1 + µ2 µ02 + 1/2 (1-µ02)(1-µ2) ] = 3/8 ( 3 + 3 µ2 µ02022 )

P1 = ± 3/2 ( µ µ0 sqrt[(1-µ02)(1-µ2)] )   (T = +, S = -)

P2 = 3/8 (1-µ02)(1-µ2)