Download Rayleigh scattering. Scattering and absorption by aerosols and

Lecture 9
Rayleigh scattering.
Scattering and absorption by aerosols and clouds.
Objectives:
1. Principles of scattering.
2. Rayleigh scattering.
3. Scattering and absorption by aerosol and cloud particles.
Required Reading:
G: 5.1-5.4, 5.6, 5.7
Additional/advanced Reading:
G: 5.5,
1. Principles of scattering.
Incident radiation
Scattered radiation
Obstacle
Consider a single arbitrary particle consisted of many individual dipoles. The incident
electromagnetic field induces dipole oscillations. The dipoles oscillate at the frequency of
the incident field and therefore scatter radiation in all directions. In a given direction of
observation, the total scattered field is a superposition of the scattered wavelets of these
dipoles, accounting for their phase difference: scattering by the dipoles is coherent (i.e.,
there is a definite relation between phases).
1
Scattering phase function P(cosΘ
Θ) is defined as a non-dimensional parameter to
describe the angular distribution of the scattered radiation as
1
4π
∫ P (cos
Θ )dΩ = 1
[9.1]
Ω
where Θ is called the scattering angle between the direction of incidence and observation.
NOTE: Another form of [9.1]
1
4π
2π
π
∫∫
0
P(cos Θ) sinθdθdϕ = 1
[9.2]
0
NOTE: The phase function is often expressed as
P (cos Θ ) = P (θ ′, ϕ ′, θ , ϕ ′)
where (θ
θ', ϕ') and (θ
θ, ϕ) are the spherical coordinates of incident beam and direction of
observation. Using the spherical geometry it can be shown (see G: Appendix 1) that
cos(Θ)
Θ) = cos(θ
θ'))cos(θ)
θ) + sin(θ
θ'))sin(θ)
θ) cos(ϕ
(ϕ'-ϕ
(ϕ ϕ)
The asymmetry factor g is defined as
1
g =
2
1
∫ P (cos Θ ) cos Θ d (cos Θ )
[9.3]
−1
g=0 for isotropic scattering
Forward scattering refers to the observations directions for which Θ < π/2:
g>0 scattering in the forward direction
Backward scattering refers to the observations directions for which Θ > π/2:
g<0 scattering in the backward direction
NOTE: Θ = 00 is often called forward scattering and Θ = 1800 is called backscattering.
2
Scattering domains:
Rayleigh scattering: 2π
πr/λ
λ <<1, and the refractive index m is arbitrary (applies to
scattering by molecules and small aerosol particles)
Rayleigh-Gans scattering: (m –1) <<1 (not useful for atmospheric application)
Mie-Debye scattering: 2π
πr/λ
λ and m are both arbitrary but for spheres only (applies
to scattering by aerosol and cloud particles)
Geometrical optics: 2π
πr/λ
λ is very large and m is real (applies to scattering by large
cloud droplets and ice crystals).
2. Rayleigh scattering
•
Because the sizes of atmospheric molecules are much smaller than the
wavelengths of solar and IR radiation, the scattering by atmospheric gases is the
Rayleigh scattering.
•
In the Rayleigh scattering approximation, a molecule (or a small particle) is
considered as an individual dipole.
Consider a small homogeneous spherical particle (e.g., a molecule) with size smaller than
r
r
the wavelength of incident radiation E 0 . Let p 0 be the induced dipole moment, then
from the classical electromagnetic theory we have
r
r
p 0 = αE0
[9.4]
where α is the polarizability of the particle.
NOTE: Do not confuse the polarization of the medium with polarization associated with
the EM wave.
The scattered electric field at the large distance r (called far field scattering) from the
dipole is given (in cgs units) by
r
r
1 1 ∂p
E = 2
sin( γ )
c r ∂t
3
[9.5]
r
where γ is the angle between the scattered dipole moment p and the direction of
observation. In oscillating periodic field, the dipole moment is given in terms of induced
dipole moment by
r
r
p = p 0 exp( − ik ( r − ct ))
[9.6]
and thus the electrical field is
r
r exp( − ik ( r − ct ) 2
k α sin( γ )
E = −E0
r
[9.7]
Decomposing the electrical vector on two orthogonal components perpendicular and
parallel to the plane of scattering (a plane containing the incident and scattering beams),
We have
E
r
= −E
E l = − E 0l
0r
exp( − ik ( r − ct ) 2
k α sin( γ 1 )
r
exp( − ik ( r − ct ) 2
k α sin( γ 2 )
r
E0r
[9.9]
pr
γ1
Dipole
E0l
[9.8]
Direction of
incident radiation
γ1=π/2; γ2=π/2-Θ
pl
γ2
Θ
Direction of
scattering
(out of page)
Using that
I=
1 c
E
∆Ω 4π
2
[9.10]
the perpendicular and parallel intensities (or linear polarized intensities) are
I r = I 0 r k 4α 2 / r 2
[9.11]
I l = I 0l k 4α 2 cos 2 (Θ) / r 2
[9.12]
4
Using that the natural light (incident beam) in not polarized (I0r =I0l=I0/2) and that
k=2π/λ
π/λ,
π/λ we have
2
I
 2 π  1 + cos ( Θ )
I = I r + I l = 02 α 2 

r
2
 λ 
4
[9.13]
Eq.[ 9.13] gives the intensity scattered by molecules (Rayleigh scattering) for
unpolarized incident light.
Rayleigh scattering phase function for incident unpolarized radiation is
P (cos( Θ )) =
3
(1 + cos 2 ( Θ ))
4
[9.14]
Eq.[ 9.13] may be rewritten in the form
I (cos( Θ )) =
I 0 2 128 π 5 P ( Θ )
α
r2
3λ 4
4π
[9.15]
Eq. [9.15] may be rewritten in the terms of the scattering cross section
I (cos( Θ )) =
I0
P (Θ )
σs
2
r
4π
[9.16]
Here the scattering cross section (in units or area) by a single molecule is
σs =α2
128 π
3λ 4
5
[9.17]
The polarizability is given by the Lorentz-Lorenz formula
α =
3
4π N s
 m2 −1
 2
m + 2



[9.18]
where Ns in the number of molecules per unit volume and m=n – ik in the refractive
index.
For air molecules in solar spectrum: n is about 1 but depends on λ, and k =0.
Thus the polarizability can be approximated as
α ≈
1
( n 2 − 1)
4π N s
5
[9.19]
Therefore, the scattering cross section of an air molecule becomes
σ
s
=
8 π 3 ( n 2 − 1) 2
3 λ 4 N s2
f (δ )
[9.20]
where f(δ
δ) is the correction factor for the anisotropic properties of air molecules, defined
as f(δ
δ) =(6+3δ
δ)/(6-7δ
δ) and δ=0.035
Using this scattering cross section, one can calculate the optical depth of the entire
atmosphere due to molecular scattering as
τ (λ ) = σ s (λ )
top
∫ N ( z ) dz
[9.21]
0
NOTE: The Rayleigh scattering cross section and hence optical depth are inversely
proportional to the fourth power of the wavelength (=> blue color of the sky)
Approximation of molecular Rayleigh optical depth (i.e., optical depth due to molecular
scattering) down to pressure level p in the Earth’s atmosphere:
p

 − 4 .15 + 0 .2 λ
τ ( λ ) ≈ 0 . 0088 
λ
 1013 mb 
•
[9.22]
Rayleigh scattering results in the sky polarization. The degree of linear
polarization is
LP ( Θ ) = −
I − Ir
Q
cos 2 Θ − 1
sin 2 Θ
= − l
=
=
I
Il + Ir
cos 2 Θ + 1 cos 2 Θ + 1
6
[9.23]
3. Scattering and absorption by aerosol and cloud particles: Mie-Debye
theory.
NOTE: Mie-Debye theory is often called Mie theory or Lorentz-Mie theory.
Mie theory outline:
Assumptions:
i) Particle is a sphere;
ii) Particle is homogeneous (therefore it is characterized by a single refractive index
m=n - ik at a given wavelength);
NOTE: Mie theory requires the relative refractive index = refractive index of a
particle/refractive index of a medium. But for air m is about 1, so one needs to know the
refractive index of the particle (i.e., refractive index of the material of which the particle
is composed).
NOTE: If a particle has complex chemical composition, the effective refractive index
must be calculated at a given wavelength.
•
Mie theory calculates the scattered electromagnetic field at all points in the
particle (called internal field) and at all points of the homogeneous medium in
which the particle is embedded. For all practical applications in the atmosphere,
light scattering observations are carried out in the far-field zone (i.e., at the large
distances from a particle).
In the far-field zone (i.e., at the large distances r from a sphere), the solution of the vector
wave equation can be obtained as
 E ls  exp( − ikr + ikz )  S 2
 s =
S
ikr
 4
 E r 
S 3   E li 
 
S 1   E ri 
i
[9.24]
i
here k =2π/λ, El and E r are the parallel and perpendicular components of incident
s
s
electrical field, and El and E r are the parallel and perpendicular components of
scattered electrical field,
7
 S 2 (Θ )
 S (Θ )
 4
S 3 (Θ ) 
is the amplitude scattering matrix (unitless)
S 1 ( Θ ) 
Θ) = S4(Θ
Θ) = 0, and thus Eq.[9.24] gives
For spheres: S3(Θ
 E ls  exp( − ikr + ikz )  S 2 ( Θ )
 s =
 0
ikr

 E r 
0   E li 
 
S 1 ( Θ )   E ri 
[9.25]
Eq.[9.25] is a fundamental equation of scattered radiation by a sphere including
polarization.
Mie theory scattering amplitudes:
S1 (Θ ) =
S 2 (Θ ) =
∞
∑
n =1
∞
∑
n =1
2n + 1
[a n π n (cos Θ ) + b nτ n (cos Θ ) ]
n ( n + 1)
[9.26]
2n + 1
[b n π n (cos Θ ) + a nτ n (cos Θ ) ]
n ( n + 1)
[9.27]
where πn and τn are Mie angular functions
π n (cos Θ ) =
1
Pn1 (cos Θ )
sin( Θ )
[9.28]
d
Pn1 (cos Θ )
dΘ
[9.29]
τ n (cos Θ ) =
where
Pn1 are the associated Legendre polynomials.
Mie theory also gives the scattering phase matrix P(Θ) that relates the Stokes
parameters {I0, Q0, U0 and V0 } of incident radiation field and the Stokes parameters {I,
Q, U and V}of scattered radiation
I o 
I 
Q 
Q 
  = σ s P o 
U  4πr 2 U 0 
 
 
V 
Vo 
8
[9.30]
0 
 P11 P12 0
P P
0
0 
12
22

P=
 0 0 P33 − P34 


 0 0 P34 P44 
[9.31]
For spheres: P22 = P11 and P 44 = P33
NOTE: In general, for a particle of any shape, the scattering phase matrix consists of 16
independent elements, but for a sphere this number reduces to four.
Thus for spheres, Eq.[9.30] reduces to
I 
Q 
  = σs
U  4π r 2
 
V 
 P11
P
 12
 0

 0
P12
P11
0
0
0
0
P33
P34
0 I o 
 
0   Q o 
− P34  U 0 
 
P33  V o 
[9.32]
where each element of the scattering phase matrix is expressed via the scattering
amplitudes S1(Θ) and S2(Θ)
P11(Θ) =P(Θ) is the scattering phase function of a particle.
From Mie theory it follows that the extinction cross-section of a particle is
σe =
4π
Re[ S (0 0 )]
2
k
[9.33]
But for the forward direction (i.e. Θ =00) from Eqs.[9.26]-9.27], we have
1 ∞
S1 (0 ) = S2 (0 ) = ∑(2n +1)(an + bn )
2 n=1
0
0
Efficiencies (or efficiency factors) for extinction, scattering and absorption are defined
as
Q
e
=
σ e
πr 2
Qs =
σs
πr 2
Q
a
=
σ a
πr 2
[9.34]
where πr2 is the particle area projected onto the plane perpendicular to the incident beam.
9
Mie efficiency factors are derived from the Mie scattering amplitude
∞
Qe =
2
x2
∑
Qs =
2
x2
∑
n =1
∞
n =1
( 2 n + 1 ) Re[ a n + b n ]
( 2 n + 1 )[ a n
2
+ bn
2
[9.35]
]
[9.36]
and the absorption efficiency can be calculated as
Qa = Qe − Qs
•
[9.37]
Extinction efficiency vs. size parameter x (assuming NO ABSORPTION):
1) small in Rayleigh limit: Q e ∝ x 4
2) largest Qe when particles and wavelength have similar size
3) Qe --> 2 in geometric limit ( x → ∞ )
Figure 9.1 Examples of Qe calculated with Mie theory for several refractive indexes.
----------------------------------------END of Mie theory outline-----------------------------------
Thus, for a single spherical particle Mie theory gives the extinction, scattering and
absorption cross-sections and efficiency factors, the scattering amplitudes and phase
matrix.
10
Strategy to compute optics of an ensemble of spherical particles:
Particle size distribution, N(r);
Refractive index, m(λ)
Mie theory
Scattering cross section, σs
Absorption cross section, σa
Extinction cross section, σe
Phase function, P11(Θ)
(as a function of a particle size
and wavelength)
Integration over size
distribution N(r)
Scattering coefficient, κs
Absorption coefficient, κa
Extinction coefficient, κe
Phase function, P11(Θ)
(as a function of the wavelength)
Integration over the size distribution:
For a given type of particles characterized by the size distribution N(r), the extinction,
scattering and absorption coefficients (in units LENGTH-1) are determined as
r2
r2
r1
r1
k e = ∫ σ e ( r ) N (r ) dr = ∫ πr 2 Qe N ( r ) dr
r2
r2
r1
r1
k s = ∫ σ s (r ) N (r )dr = ∫ πr 2 Q s N ( r )dr
r2
r2
r1
r1
k a = ∫ σ a (r ) N ( r )dr = ∫ πr 2 Qa N (r )dr
11
[9.38]
[9.39]
[9.40]
The single scattering albedo gives the percentage of light which will be scattered in a
single scattered event and it is defined as
ω0 =
ks
ke
[9.41]
NOTE: No absorption (conservative scattering): ω0 = 1
No scattering: ω0 = 0
Scattering phase function of particles characterized by the size distribution N(r)
r2
P ( Θ) =
∫ P(Θ, r )σ
s
N (r )dr
r1
[9.42]
r2
∫σ
s
N ( r ) dr
r1
Figure 9.2 Examples of scattering phase functions at a wavelength of 0.5 µm.
12
Optical properties of the external mixture (i.e., the mixture of several types of
particles)
ke = ∑ k ei
k s = ∑ k si
i
where
i
kei , k si and kai
ka = ∑ kai
[9.43]
i
are calculated for each particle type characterized by its particle
size distribution Ni(r) and a refractive index (or effective refractive index) mi.
NOTE: Do not sum the single scattering albedo and scattering phase functions !!!
•
In general, an atmospheric layer has molecules, aerosols and/or cloud particles.
Thus, one needs to calculate the effective optical properties of this layer as an
external mixture of the optical properties of these constituents.
How to calculate the effective optical properties of an atmospheric layer consisting of
gas and aerosols (or clouds):
Effective (also called total) optical depth:
τ λ = τ aM,λ + τ sM,λ + τ aA,λ + τ sA,λ
where
τ aM,λ
and
τ sM,λ
[9.44]
are optical depth due to absorption by gases and molecular
(Rayleigh) scattering, respectively;
τ aA,λ
and τ s ,λ are optical depth due to absorption and scattering by aerosol particles,
A
respectively.
Effective single scattering albedo:
τ sM,λ + τ sA,λ
τλ
[9.45]
τ sM,λ PλM ( Θ ) + τ sA,λ PλA ( Θ )
Pλ ( Θ ) =
τ sM,λ + τ sA,λ
[9.46]
ω 0 ,λ =
Effective scattering phase function:
13
Effective asymmetry parameter:
gλ =
τ sA,λ g λA
τ sM,λ + τ sA,λ
[9.47]
Optics of cloud particles
•
In contrast to atmospheric aerosols, the cloud droplets consist of a single species –
water. Thus, to calculate the absorption and scattering cross sections of cloud
droplets, one needs to know the size of droplets and the refractive index of water
versus wavelength.
NOTE: Cloud droplet sizes vary from a few micrometers to 100 micrometers with
average diameter in 10 to 20 µm range.
Using the effective radius (see lecture 8)
πr N (r )dr
=∫
∫ πr N (r )dr
3
re
2
and the liquid water content (LWC) (see Lecture 8):
LWC = ρ wV =
4
ρ w ∫ π r 3 N (r )dr
3
the extinction coefficient of cloud droplets is
k e = ∫ σ e ( r ) N ( r )dr = ∫ Qeπr 2 N ( r )dr
and that Qe ≈ 2 for water droplets at solar wavelengths, we have
ke ≈
3 LWC
2 re ρ w
14
[9.48]
Example of cloud optical properties
Figure 9.3 Example of optical properties of typical cumulus and stratus clouds (for a
cloud droplet size distribution reff = 6 µm) . Here the normalized extinction coefficient is
k e ( λ ) / k e ( 0 . 5 µ m ) and k e ( 0 . 5 µ m ) = 42 . 8 km − 1 .
15
Differences in optics of water drops and ice particles are due to
i)
Particle shapes;
ii)
Differences between the refractive indices of water and ice.
Approaches to predict ice crystal optics:
i)
Ray-optics approach (also called ray tracing approach) based on geometrical
optics for the size parameters (2πr/λ) of 15-20 and larger;
ii)
Finite differences time domain (FDTD) method for the size parameters
smaller than about 20;
iii)
Laboratory measurements of optical properties (e.g., phase function) of the ice
crystals of known shapes and sizes.
16