1 / 46

ISSI Meeting, Bern 19-23 January 2014

ISSI Meeting, Bern 19-23 January 2014. New 3D photochemical global model with ions in D-region: The instrument for solar-atmospheric relations study Alexei Krivolutsky Lidiya Cherepanova, Tatyana Vyushkova, and Alexander Repnev Laboratory for Atmospheric Chemistry and Dynamics

Download Presentation

ISSI Meeting, Bern 19-23 January 2014

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. ISSI Meeting, Bern 19-23 January 2014 New 3D photochemical global model with ions in D-region: The instrument for solar-atmospheric relations study Alexei Krivolutsky Lidiya Cherepanova, Tatyana Vyushkova, and Alexander Repnev Laboratory for Atmospheric Chemistry and Dynamics Central Aerological Observatory, Dolgoprudny, Moscow Region Russia

  2. Outline 1. Model description 2. Results of simulations: - neutral compounds; - electrons; - other ions. 3. Effects of solar cycle 4. Conclusions

  3. MODELCHARM – I (CHemical Atmospheric Research Model with Ions)(Krivolutsky et al., 2015) Heights: 0-90 km P, L – photochemical sources and losses U, V, W – wind components, µ – mixing ratio number of species: neutrals – 41; ions – 23 number of chemical reactions (total): 194 Photodissociation and ionization rates (total) : 48 Methods:“chemical families” for neutrals(Turco, Whitten, 1974) “electroneutrality” for ions Prather’s scheme for advection ( Prather, 1986) Resolution: 2 km Х 5 Х 5 deg., Time step: 100 s

  4. List of neutral species • “Families” • Ox = O3 + O(3P) + O(1D); • NOy =N + NO + NO2 + NO3 + 2N2O5 + HNO3 + HO2NO2 + ClNO3+N(2D); • Cly=Cl + ClO + OClO + ClOO + HOCl + HCl; • HOx=H + OH + HO2 + 2H2O2 ; others • CH3, CH2O, CH3O2, CH3O2H, CH3O, CHO, CO. • О2(1g) • Source-gases • CH4, CO2, N2O, СF2Cl2, CFCl3, H2, Cl4, Cl2, СН3Cl, CH2Cl, • О2, N2 (fixed profiles), H2O(fixed global field/HALOE).

  5. PHOTODISSOCIATION RATES(CHARM-I) O2+h O+O(1D) N2O5+h NO2+NO3 H2O+h H+OH O2+h O+O HNO3+h OH+NO2 CF2Cl2+h products O3+h O+O2 CLONO2+h Cl+NO3 CFCl3+h products O3+h O(1D)+O2 HCl+h H+Cl CH4+h CH3+H H2O2+h OH+OH ClO+h Cl+O CH4+h CH2+H2 NO2+h NO+O(1 D) NO3+h NO+O2 CCl4+h products HNO3+h H+NO3 H2O2+h OH+OH CH3Cl+h CH3+Cl HOCl+h Cl+HO CO2+h CO+O N2O+h N2+O(1D) N2O5+h 2NO2+O HO2NO2+h HO2+NO2 Cl2+h Cl+Cl NO+h N+O NO2+h NO+O NO3+h NO2+O

  6. List of ionized compounds • Positive: • O2+ O4+ O2+(H2O) H+(H2O) H+(H2O)3 • H+(H2O)4 H+ (H2O)2 • NO+N2 NO+CO2 NO+(H2O) NO+(H2O)2 NO+(H2O)3 NO+ • Negative: [e] • O2 - O3- O4- CO4- O- OH- CO3- O2 -(H2O)HCO3-

  7. Ionization • 1-10 nm ( X-Rays) • 102,7-111,8 nm О2(1g) • q(z)=n(O2(1g))0,54910-9exp(-2,406 10-20 • N(O2)+2,61410-9exp(-8,50810-20N(O2)) • 121,6 nm ( Lα ) NO • GCRs (Heaps, 1978)

  8. ARM -Atmospheric Research Model (GCM)(Krivolutsky et al., 2012) Altitudes:0-135 км Resolutions:vertical– 1 km; longitudinal – 100; latitudinal– 50 time step – 5 min. Paramaterizations: Heating -О2, О3, Н2О (Strobel, 1978; Chou et al., 2002); IR cooling-СО2, О3, H2O, NО (Chou et al., 2002; Fomichev, 2003; Kockarts, 1980), IGWs (Lindzen, 1981) Planetary waves at lower boundary(S=1,2.3)

  9. Global temperature field for July (К)(Krivolutsky et al, 2012)

  10. Zonal wind structure (m/s)for July (Krivolutsky et. al, 2012) Height (km) latitude

  11. Amplitude of D tidal component in zonal wind (m/s)July (Krivolutsky et al., 2012)Height (km) latitude

  12. Amplitude of SD tidal component in zonal wind (m/s) July(Krivolutsky et al., 2012)Height (km) latitude

  13. Neutral compounds: O3 ( ppmv) January

  14. Neutral compounds: NOy ( ppbv) January(CHARM-I)

  15. Neutral compounds: HNO3 ( ppbv) January(CHARM-I)

  16. Neutral compounds: N2O ( ppbv) January

  17. Calculated ionization rate byLα

  18. Electron density,80 km (number/cm**3) 1st January (00:00 UT)LatitudeLongitude 50 0 -50 0 50 100 150 200 250 300 350

  19. IONS: electron density,60 km (number/cm**3) 1st January (00:00 UT)

  20. NO+ (number/cm**3) at 80 km 1st January (00:00 UT)

  21. NO+ (number/cm**3) at 70 km 1st January (00:00 UT)

  22. O2+ (number/cm**3) at 80 km 1st January (00:00 UT)

  23. O2+ (number/cm**3) at 70 km 1st January (00:00 UT)

  24. O-2 (number/cm**3) at 80 km1stJanuary (00:00 UT)

  25. O-2 (number/cm**3) at 60 km1stJanuary (00:00 UT)

  26. 1st January (noon)

  27. 1st January (noon)

  28. 1st January (noon)

  29. electrons (number/m**3) 45 N (noon)

  30. POSITIVE IONS: O2+ (number/m**3) 45 N (noon)

  31. POSITIVE IONS: H+(H2O)4(number/m**3) 45 N (noon)

  32. POSITIVE IONS: O2+H2O45 N (noon)

  33. Effects of solar cycle simulated • with CHARM-I

  34. Solar cycle in UV radiation(Matthew et al., 2012)

  35. Solar UV variations (164,5 нм)

  36. Solar UV spectrum variations used in model runs

  37. Ozone change (%) max-min

  38. Ox change (%) max-min

  39. HOx change (%) max-min

  40. NOy change (%) max-min

  41. Simulated changes in electron density (%)between max. and min. of solar cycle(January)

  42. Simulated changes in NO+ (%)between max. and min. of solar cycle(January)

  43. Simulated changes in NO+ (H2O) (%)between max. and min. of solar cycle (January)

  44. Simulated changes in O2+(%)between max. and min. of solar cycle(January)

  45. Concluding remarks 1. It seems that CHARM-I reproduces ion and neutral composition well. 2. UV variations disturb neutrals ( ozone etc) and ion composition due to its interactions. 3. Solar cycle in ionization was included only by Lα , Lβ 4. Effect of particles will be included.

  46. Thank you for your attention!

More Related