* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Download Rayleigh scattering. Scattering and absorption by aerosols and
Nonimaging optics wikipedia , lookup
Ellipsometry wikipedia , lookup
Atomic absorption spectroscopy wikipedia , lookup
Harold Hopkins (physicist) wikipedia , lookup
Retroreflector wikipedia , lookup
X-ray fluorescence wikipedia , lookup
Vibrational analysis with scanning probe microscopy wikipedia , lookup
Birefringence wikipedia , lookup
Magnetic circular dichroism wikipedia , lookup
Phase-contrast X-ray imaging wikipedia , lookup
Ultraviolet–visible spectroscopy wikipedia , lookup
Anti-reflective coating wikipedia , lookup
Resonance Raman spectroscopy wikipedia , lookup
Nonlinear optics wikipedia , lookup
Optical tweezers wikipedia , lookup
Dispersion staining wikipedia , lookup
Sol–gel process wikipedia , lookup
Atmospheric optics wikipedia , lookup
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