Rotational Raman scattering in the Earth's atmosphere ... - SRON
Rotational Raman scattering in the Earth's atmosphere ... - SRON
Rotational Raman scattering in the Earth's atmosphere ... - SRON
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
A doubl<strong>in</strong>g-add<strong>in</strong>g method to <strong>in</strong>clude multiple orders of rotational <strong>Raman</strong> <strong>scatter<strong>in</strong>g</strong> 37<br />
where σ ram,i (ν;ν ′ ) is <strong>the</strong> rotational <strong>Raman</strong> cross section for each l<strong>in</strong>e, b i (J → J ′ ) are <strong>the</strong> so-called<br />
Placzek-Teller coefficients, J and J ′ are <strong>the</strong> rotational quantum numbers of <strong>the</strong> <strong>in</strong>itial state and f<strong>in</strong>al<br />
state respectively, and f i (T,J) is fractional population of <strong>the</strong> <strong>in</strong>itial rotational state where T is <strong>the</strong><br />
temperature of <strong>the</strong> gas [Penney et al., 1974, Jo<strong>in</strong>er et al., 1995, Sioris, 2001]. For our calculations<br />
we use a maximum rotational number of <strong>the</strong> <strong>in</strong>itial state of J = 50. The wavenumber shifts can be<br />
calculated easily from <strong>the</strong> energy differences between <strong>the</strong> <strong>in</strong>itial and f<strong>in</strong>al state.<br />
The Rayleigh, Cabannes, and <strong>Raman</strong> s<strong>in</strong>gle <strong>scatter<strong>in</strong>g</strong> albedo can be determ<strong>in</strong>ed from <strong>the</strong> cross<br />
sections as,<br />
ω ray,air (ν ′ ) = σ ray,air (ν ′ )/σ ext (ν ′ ),<br />
ω cab,air (ν ′ ) = σ cab,air (ν ′ )/σ ext (ν ′ ),<br />
ω ram,air (ν;ν ′ ) = σ ram,air (ν;ν ′ )/σ ext (ν ′ ),<br />
where <strong>the</strong> total ext<strong>in</strong>ction cross section,<br />
(2.30)<br />
σ ext (ν ′ ) = σ ray,air (ν ′ ) + σ abs (ν ′ ) (2.31)<br />
is <strong>the</strong> sum of <strong>the</strong> <strong>scatter<strong>in</strong>g</strong> cross section of all scatterers plus <strong>the</strong> absorption cross section of all<br />
absorbers.<br />
Besides <strong>the</strong> s<strong>in</strong>gle <strong>scatter<strong>in</strong>g</strong> albedo, <strong>the</strong> phase function is needed to describe a <strong>scatter<strong>in</strong>g</strong> event.<br />
The phase functions of Rayleigh, Cabannes, and <strong>Raman</strong> <strong>scatter<strong>in</strong>g</strong> depend only on <strong>the</strong> wavenumber<br />
of <strong>the</strong> <strong>in</strong>com<strong>in</strong>g light and <strong>the</strong> <strong>scatter<strong>in</strong>g</strong> angle Θ, which is related to zenith and azimuthal angles µ,<br />
µ ′ , ϕ and ϕ ′ as<br />
cos Θ = µµ ′ + √ 1 − µ 2√ 1 − µ ′2 cos[ϕ − ϕ ′ ]. (2.32)<br />
The effective Rayleigh and Cabannes phase functions for air are given by<br />
P ray,air (ν ′ , Θ) =<br />
∑ [ ]<br />
σray,i (ν ′ )<br />
χ i P<br />
σ ray,air (ν ′ ray,i (ν ′ , Θ), (2.33)<br />
)<br />
i=N 2 ,O 2 ,Ar<br />
P cab,air (ν ′ , Θ) =<br />
∑ [ ]<br />
σcab,i (ν ′ )<br />
χ i P<br />
σ cab,air (ν ′ cab,i (ν ′ , Θ), (2.34)<br />
)<br />
i=N 2 ,O 2 ,Ar<br />
with<br />
[ ] [ 135 +<br />
P ray,i (ν ′ 39εi (ν ′ ) 135 + 3εi (ν ′ )<br />
, Θ) =<br />
+<br />
180 + 40ε i (ν ′ ) 180 + 40ε i (ν ′ )<br />
[ ] [ 540 +<br />
P cab,i (ν ′ 39εi (ν ′ ) 540 + 3εi (ν ′ )<br />
, Θ) =<br />
+<br />
720 + 40ε i (ν ′ ) 720 + 40ε i (ν ′ )<br />
and <strong>the</strong> phase function for <strong>Raman</strong> <strong>scatter<strong>in</strong>g</strong> is given by<br />
]<br />
cos 2 Θ, (2.35)<br />
]<br />
cos 2 Θ, (2.36)<br />
P ram,air (Θ) = 39<br />
40 + 3 40 cos2 Θ, (2.37)<br />
which clearly is <strong>in</strong>dependent of wavenumber and molecular species (e.g. Stam et al. [2002] and references<br />
<strong>the</strong>re<strong>in</strong>).