Correct atmospheric optics modelling: Theory and Experiment Irina Melnikova Observatory of Environmental Safety Resource Center, Research Park St.Petersburg State University St.Petersburg. Russia. irina.melnikova@pobox.spbu.ru 1
Objectives: 1. Constructing the simple optical model of homogeneous atmosphere 2. Solution of the direct problem of atmospheric optics with operative varying optical parameters for elucidating the interaction between key atmospheric parameters and radiative characteristics 3. The solution of the direct problem is calculation of radiation characteristics (radiant heat. radiation balance at the tropopause level) 4. Comparison with the results of airborne measurements 2
Optical parameters Optical thickness of the clear atmosphere : τ = τ a.sc + τ a.ab ab + τ R + τ m.ab. τ a.sc the optical thickness of aerosol scattering. τ a.ab ab - the optical thickness of aerosol absorption. τ R - the optical thickness of molecular scattering. τ m.ab - the optical thickness of molecular absorption; Cloud optical thickness τ Single scattering albedo: for clear atmosphere ω = (τ a.sc + τ R )/τ ; For cloudy atmosphere ω = (τ + τ a.sc + τ R )/( τ + τ ) ; The phase function asymmetry parameter: g=. for clear atmosphere and g=.85 for cloud atmosphere; The ground albedo A s 3
1..9 y = 7E-5x +.889.8.7 Ground albedo.6.5.4.3 y = -9E-7x 2 +.1x -.456 WATER SAND SNOW.2.1. 3 5 7 9 Wavelength. nm y = 3E-7x 2 -.x +.216 Spectral dependence of the ground albedo of different surfaces from observational data processing (C.A. Varotsos. I.N. Melnikova. A.P. Cracknell. C. Tzanis. A.V. Vasilyev. New spectral functions of the near-ground albedo derived from aircraft diffraction spectrometer observations. Atmospheric Chemistry and Physics. v. 13. pp. 16211-16245) 4
Data of the SPSU RC lidar is used for aerosol optical thickness modelling Elastic channels Aerosol lidar Raman channels Multiwavelength lidar Polarization filter 355 nm Doppler (wind) lidar Tunable titan-sapphier lasers (at the mobile complex) 5
Stationary lidar system: A Doppler lidar for measuring the wind speed and direction up to 12 km height An aerosol lidar for measuring the atmospheric aerosol parameters up to 25 km height 1. Provide regular monitoring the dynamics of an atmospheric pollution above the big city center. 2. Retrieving atmospheric dust parameters: size, extinction coefficient, backscattering coefficient, real and imagine parts of the refractive index, and content
Stationary lidar system
Aerosol lidar 164 nm - 4 mj 532 nm - 16 mj 355 nm - 1 mj Doppler (wind) lidar for wind velocity and direction profile Pulse repetition rate 1kHz 8
signal 78 76 74 72 7 68 66 64 62 6 58 56 54 52 5 48 46 44 42 4 38 36 34 32 3 28 26 24 22 2 18 16 14 12 1 8 6 4 2 Lidar 38 Screenshot of the received Lidar signal 31.1.213 1 355_Anlg 2 355_PhCt 3 355_De_Anlg 4 355_De_PhCt 5 532_Anlg 6 532_PhCt 7 387_Anlg 8 387_PhCt 9 68_Anlg 1 68_PhCt 11 164_Anlg 12 48_PhCt E+ 2E+3 4E+3 6E+3 Distance [m] 8E+3 1E+4 1.2E+4 1.4E+4 9
Lidar sounding above St. Petersburg city The extinction coefficient above St. Petersburg 25 March 213. The vertical profile till 4 km and 25 km during 1 hour λ=532 nm The maximum of hat pollution at.7 km, disappeared during 45 min (15:3-16:15) α(z).56km -1 The stratosphere aerosol Yunge layer at 17-22km α(z).15km -1 Total optical thickness is τ =.735 1
Mol scatt VALUES OF OPTICAL PARAMETERS OF THE CLEAR ATMOSPHERE λ. µm.3.3.4.5.6.7.8.9 τ Rel(z=) 1.222.364.14.67.36.21.13 Aer τ a scatt. τ a abs..25 Aer I τ a scatt τ a abs.1 Aer III τ a scatt.7 τ a abs.4.26184 I ω.29876 III.33328 I III τ...18...16...14...13...12...12.1.1.1.1.1.1.1 4.667 4.927 5.767.5.4 1..98195.68354.364.554 1.264.5.4 1..96774.61539.14.31 1.4.5.4.91781.93243.58757.72.222.965.4.4 1..94318.52153.36.176.836.3.4 1..93378.44522.21.151.721.3.4 1..937.43899.13.143 11.713
Clear atmosphere 2. Optical thickness.6 Single scattering co-albedo Optical thickness 1.6 1.2.8.4 Sing gle scattering co-albedo.5.4.3.2.1..3.5.7.9 Wavelength. µm A-1 A-2 A-3 Thin lines optical thickness of scattering Thick lines optical thickness of absorption..3.5.7.9 Wavelength. µm A A1 A3 A2 Spectral dependence of the single scattering co-albedo for 4 aerosol models 12
Clear atmosphere. Lidar sounding Dynamics of the variation of aerosol extinction from lidar observation in SPSU (Donchenko V.K.. Samulenkov D.A.. Melnikova I.N.. Boreysho A.S.. Chugreev A.V. Laser systems of the St.Petersburg State University Resources Center. Possibilities. Problems Statement and the First Results. The contemporary problems of the Earth remote sensing form the Space. Moscow. 213. Том 1. 3. p 122-134) 13
Experimental values of the aerosol optical thickness in St.Petersburg and suburbs Wavelength. nm Experiment St. Petersburg city. Lidar sounding The Ladoga Lake Airborne observations Peterhoff city. Ground observations 355 532 7.136.842.1.6.15.3.7.9 14
Calculation of radiative characteristics Clear atmosphere. Radiative divergence Radiative div vergence. Wcm -2 µm -1 6 5 4 3 2 1 Aerosol Аs=..5..9.2.3.4.5.6.7.8.9 1. Wavelength. µm A= A=.5 Simulation and airborne observation of radiative divergence for models of Aerosol and 1 15
Clear atmosphere. Radiative divergence (continuation) Radiat tive divergence. Wcm -2 µm -1 3 2 1 Aerosol 3.2.4.6.8 1 Wavelength. µm А=, 8 Kara-Kum 12.1.1983, 7 Simulation (Aerosol 3) and airborne observation of radiative divergence after the sand storm (Melnikova I.. Vasilyev A. Short-wave solar radiation in the Earth atmosphere. Calculation. Observation. Interpretation. Springer- Verlag. Heidelberg. 24. 35 p.) 16
CLOUD 1 τ = 5 and 1 for all wavelength. added to the scattering optical thickness of the clear atmosphere CLOUD 2 2-layer atmosphere : cloud 1 (in layer -1 km) + clear layer above the cloud The partly scattered light falls to cloud top and cloud spherical albedo is assumed as ground albedo for above - cloud layer CLOUD 3 τ (λ) - Spectral dependent optical thickness 17
OPTICAL PARAMETERS OF THE CLOUD-1 MODEL I II III I II III I II III I II III λ. µm.3.4.5.6.7.8.9 τ Σ (τ =1) 14.667 1.364 1.14 1.117 1.36 1.21 1.13 14.927 1.554 1.31 1.287 1.176 1.151 1.143 CLOUD 1 15.78 1.68 1.43 1.427 1.296 1.261 1.253 15.767 11.264 11.4 11.37 1.836 1.721 1.713 ω (τ =1).76512 1 1.9965 1 1 1.76854.9995.9993.99222.9992.9992.9991 CLOUD 1.76887.99626.99617.98945.99612.9961.9961 τ Σ (τ =5) CLOUD 1 ω (τ =5) CLOUD 1.75614 9.667 9.927 1.78 1.767.64363.65196.6542.64289.96449 5.364 5.554 5.68 6.264 1.9982.99296.93614.96377 5.14 5.31 5.43 6.4 1.99812.99263.93378.95742 5.117 5.287 5.427 6.37.99218.98487.97973.92215.9639 5.36 5.176 5.296 5.836 1.9987.99245.93146.96269 5.21 5.151 5.261 5.721 1.9986.9924.938.96266 5.13 5.143 5.253 5.713 1.9986.9924.938 18
I II III OPTICAL PARAMETERS IN THE CLEAR ABOVE-CLOUD LAYER (P Z=1 λ, µm.3.4.5.6.7.8.9 τ Σ 4.545 4.85 4.956 5.645.328.518.644 1.228.13.3.42 1.3.11.215.365.965.32.172.292.832.19.149.259.719 =1KM).12.142.252.712 III 5.645 1.228 1.3.965.832.719.712 I II III ω.2422.2897.29681.31887 1.987.93789.67427 1.96667.9476.61165.54546.9323.84932.5831 1.94186.8631.51923 1 1.93289.92958.84556.84127.44367.4382 19
OPTICAL PARAMETERS OF THE CLOUD-3 MODEL λ, µm.3.4.5.6.7.8.9 τ scattτ 58 25 16 12 1 1 1 Rel τ (z>) 1.222.364.14.67.36.21.13 scatt τ Σ τ Σ I II III 62.667 25.364 16.14 12.72 1.36 1.21 1.13 62.85 25.518 16.3 12.21 1.172 1.149 1.142 62.956 25.644 16.42 12.35 1.292 1.259 1.252 63.645 26.228 17.3 12.96 1.832 1.719 1.712 I II III ω.94527.944988.944644.93959 1.99968.99844.984749 1.99939.99756.97651.99959.99918.99676.96914 1.9992.99611.9637 1.9992.9961.96268 1.9991.9961.96266 2
14 12 Cloudy atmosphere 12 Optical thickness 1 Optical thickness for Cloud 1 (red) Cloud 3 (lila) Optical thickness 1 8 6 4 2 Op ptical thickness 8 6 4 2.2.4.6.8 1 Wavelength. µm А А1 А2 А3.2.3.4.5.6.7.8.9 1 Wavelength, µm Optical thickness of Cloud-1 model (τ =1) (upper group of curves) and above-cloud atmosphere (lower group of curves) for 4 Aerosol models А Cl1 А3 Cl1 А3-CL3 Optical thickness of Cloud-1 and Cloud 3 models 21
Single scattering co-albedo.14.12.1.8.6.4 Cloudy atmosphere (continuation) Single scattering co-albedo Single scattering co-albedo for 4 aerosol models and Cloud-1 model with optical thickness 5 and 1 and.2 experimental data from..3.5.7.9 Wavelength. µm TAU=1,A TAU=1,A1 TAU=1,A2 TAU=1,A3 3: 1 Apr.1971 4: 5 Oct.1972 5: 5 Dec.1972 9: 1 Oct.1972 11: 3 May 1976 6: 24.9.1972 7: 2.4.1985 8: 12.4.1996 1: 12 July 1974 2: 4 Aug 1974 NASA 13.9.2 (T) TAU=5,A TAU=5,A1 TAU=5,A2 (Melnikova I. Vasilyev A. Short-wave solar radiation in the Earth atmosphere. Calculation. Observation. Interpretation. Springer-Verlag. Heidelberg. 24. 35 p.) 22
Cloudy atmosphere. Radiative divergence Radiative di ivergence. Wcm -2 µm -1 25 2 15 1 5 CLOUD 1. 2; Aerosol 1 A=.25.45.65.85 Wavelength. µm Aerosol 2 Cloud 1, 6 Ladoga, 24,9,1972, 63 Cloud 2, A=, 6 Simulation for Aerosol 1. Cloud 1 (red) and 2 (green) models and results of airborne observation of the radiative divergence above the Ladoga Lake (Melnikova I.. Vasilyev A. Short-wave solar radiation in the Earth atmosphere. Calculation. Observation. Interpretation. Springer-Verlag. Heidelberg. 24. 35 p.) 23
Cloudy atmosphere. Radiative divergence (continuation) Radiative divergence. W cm -2 mm -1 6 5 4 3 2 1 CLOUD 1. Aerosol 1. 3; A=..9 Simulation for Aerosol 1 and 3). Cloud 1 model and results of airborne observations of radiative divergence after the sand storm (Sakhara dust) above the Atlantic Ocean and in clean atmosphere above the Ladoga Lake.3.5.7.9 1.1 Wavelength. µm Aerosol 2, A=,9, 6 GATE, 4,8,1974, 15 GATE, 12,7,1974, 15 Ladoga, 2.4.1985, snow, 5 Aerosol 3, A=,9, 6 (Melnikova I.. Vasilyev A. Short-wave solar radiation in the Earth atmosphere. Calculation Observation. Interpretation. Springer-Verlag. Heidelberg. 24. 35 p.) 24
1..9 Cloudy atmosphere. Relative radiative divergence. Cloud - 3 model Relat tive radiative divergence. r.u...8.7.6.5.4.3.2.1..3.5.7.9 Wavelength. µm GATE 12.7.1974 GATE 4.8.1974 Black Sea, 1.4.1971 Azov Sea 5.1.1972 Aer 1, A=, Div, 2 Aer 3, A=, Div, 2 Aer 2, A=, Div, 2 Simulation (Aerosol 1.2 and 3) for Cloud 3 model and airborne observation of relative radiative divergence in cloudy atmosphere (Melnikova I.. Vasilyev A. Short-wave solar radiation in the Earth atmosphere. Calculation. Observation. Interpretation. Springer-Verlag. Heidelberg. 24. 35 p.) 25
Local instantaneous radiative forcing (variations of the net flux at the troposphere top) Clear atmosphere (1-F ) forsing Aerosol Cloud 1 (Smoothed cloud) (1-F ) forsing forsing Aerosol cloud Cloud 2 (2-layer atmosphere) (1-F ) forsing forsing Aerosol cloud Aer=.761..4843. -.2767.7743. +.133 W/m 2-5.152 +2.41 Aer=1.758 -.24.4831 -.12 -.2755.534 -.243 -.2246 W/m 2 6 -.435 -.2175-49.93-43.55-4.71 Aer=2.768 +.65.492 +.77 -.2755.584 -.1949 -.1835 W/m 2 +1.178 +1.396-49.93-35.14-33.26 Aer=3.863.121.6788 +.1945 -.1843.8219 +.476 -.412 W/m 2 +18.5 +35.25-33.4 +8.628-7.468 f aeros =[(1-F ) Aer - (1-F )]F f = cloud [(1-F ) clear -(1-F )]F 26
Estimating the heating rate of the atmospheric layer in the shortwave range S λ = 1 J/(s m 2 ) - the solar constant in shortwave range (.3 1. µm); r = 1 kg/m 3 - the air density at the level 8 mb; C p = 15 J/(kg deg) - the specific heat of the dry air in clear atmosphere C p = 1952 J/(kg deg) - the specific heat of water vapor at constant pressure; dt dt S rc dr dz C p = 4218 J/(kg deg) - the specific heat of liquid water at C; The average value C p = 2392 J/(kg deg) in clouds. Model CLEAR p = λ p dt/dt. dt. degree / day CLOUD A S.9.9 Aerosol 1 2.7 3. 2.7 2.9 Aerosol 3 12.3 24.5 5.6 8.3 27
CONCLUSIONS: 1.Lidar sounding in the Research Park of SPSU provides the construction of adequate optical models of the atmosphere 2.The simplest optical model provides suitable results of radiative characteristics calculation that shows an acceptable accordance with airborne observation 3.These model allows clearly demonstrate influence of chosen optical parameter on radiative characteristics. 4.The presence of aerosols in the atmosphere greatly affects the optical and radiative properties of clouds 5.Even the simple models confirm that simulation of the atmosphere optical and radiative characteristics should accurate account for atmospheric pollution and correct 28 forecast of global environmental changes
Research Park (1 MEuros, 21-212)
Research park portal
Unified electronic system for applications One time registration Possibility to work with all the centers "in one window The ability to track the status of requests and keep track of time of its execution Remote getting results List of publications using equipment Research park System of assessments and reviews to improve the work of the Research park http://researchpark.spbu.ru
Research Park of SPSU READY FOR JOINT RESEARCH PROJECTS! http://researchpark.spbu. ru/en
Thanks for your attention 33