ISSI Meeting - Bern, Switzerland, 23-28 March 214 3D simulations of atmospheric response caused by precipitating electrons and solar protons at both polar regions induced by geomagnetic storms Alexei Krivolutsky alexei.krivolutsky@rambler.ru Laboratory for Atmospheric Chemistry and Dynamiics Central Aerological Observatory (CAO), Dolgoprudny, Russia
Outline 1. Particles-Ionization-Chemistry mechanism. 2. Ionization rates caused by the particles (October-November, 23; HEPPA). 3. Models description 4. Particle effects in composition 5. Particle effects in temperature and dynamics 6. Conclusions
Daily averaged energy values of solar protons during 1994-23 (1-1 MeV, GOES data) Averaged energy of protons (MeV) 18 17 16 15 14 13 12 11 1 9 8 7 6 5 4 3 2 1 July 14 Nov. 4 Oct. 28 1994 1995 1996 1997 1998 1999 2 21 22 23
N and OH production caused by cosmic rays (Jackman et al., 1976; Heaps, 1978; Solomon and Crutzen, 1981) 7 km 5 N OH 3 1. 1.2 1.4 1.6 1.8 2. molecular / ion pair
Ionization rates 75 N (66 km) caused by energetic particles during October-November 23 (calculated by Maik Wissing) 1E+11 1E+1 [ p ] [ e ] [ alfa ] Ionization raites 1E+9 1E+8 1E+7 75 N 66 km 1E+6 1E+5 1 2 3 5 7 9 Hours after 26 October 23 (: UT)
Ionizatin raites Ionization rates caused at 75 S (66 km) by energetic particles during October-November 23 (calculated by Maik Wissing) 1E+11 1E+1 1E+9 1E+8 1E+7 [ p ] [ e ] [ alfa ] 75 S 66 km 1E+6 1E+5 2 1 Hours after 26 October 23 (: UT)
Calculated Ionization rates caused by solar protons at 66 km (6 AM 28 October 23) 6H 28oct p 35 longitude 3 25 2 15 1 5 - - - -2 2 latitude 1.5E+1 9.E+9 7.E+9 5.E+9 3.E+9 1.E+9.E+
longitude Calculated Ionization rates caused by solar protons at 66 km (6 PM 28 October 23) 35 3 25 2 15 1 5 18H 28oct p - - - -2 2 latitude 1.4E+9 1.2E+9 1.E+9 8.E+8 6.E+8 4.E+8 2.E+8.E+
Calculated Ionization rates caused by relativistic electrons at 66 km (6 AM 28 October 23) 35 6H 28oct E longitude 3 25 2 15 1 5 - - - -2 2 latitude 1.4E+1 1.2E+1 1.E+1 8.E+9 6.E+9 4.E+9 2.E+9.E+
Calculated Ionization rates caused by electrons at 66 km (6 PM 28 October 23) 35 18H 28oct E longitude 3 25 2 15 1 5 - - - -2 2 latitude 1.7E+1 1.5E+1 1.3E+1 1.1E+1 9.E+9 7.E+9 5.E+9 3.E+9 1.E+9.E+
Ionization rates induced by solar protons at 75 N km 7 5 3 2 1 1 2 3 5 7 hours after 27 October 23 1.8E+1 1.6E+1 1.4E+1 1.2E+1 1.E+1 8.E+9 6.E+9 4.E+9 2.E+9.E+
Ionization rates caused relativistic electrons at 75 N km 7 5 3 2 1 4.E+9 3.6E+9 3.2E+9 2.8E+9 2.4E+9 2.E+9 1.6E+9 1.2E+9 8.E+8 4.E+8.E+ 1 2 3 5 7 hours after 27 October 23
Ionization rates caused by relativistic electrons at 85 S IO_e_85S km 7 5 3 2 1 1 2 3 5 7 hours after 27 October 23 9.E+9 8.E+9 7.E+9 6.E+9 5.E+9 4.E+9 3.E+9 2.E+9 1.E+9.E+
CHARM CHemical Atmospheric Research Model equations (-9 km): Model µ µ µ µ + U + V + W = P L AD µ t a cosθ λ a θ z P, L photochemical Sources and Losses U, V, W wind components (from ARM) Number of gas-phase reactions 73; number of photolysis reactions - 38 Chemical families method - Turco, Whitten, 1974 Advection scheme: Prather, 1986
Ozone mixing ratio (ppmv) for 1st January (3D simulation with CAO model) 7 5 3 2 1 - - - -2 2 1 8 6 4 2
NO mixing ratio (ppbv) for 1st January (3D simulations with CAO model) 7 5 3 2 1 - - - -2 2 2 15 1 7 5 2 1
HNO3 mixing ratio (ppbv) for 1st January (3D simulations with CAO model) 7 5 3 2 1 - - - -2 2 1. 8. 6. 4. 2..
N2O mixing ratio (ppbv) for 1st January (3D photochemical simulations) 7 5 3 2 1 - - - -2 2 3 2 22 1 1 1 2 5
ClNO3 mixing ratio (ppbv) for 1st January (3D photochemical simulations) 7 5 3 2 1 - - - -2 2 2.4 2.2 2. 1.5 1..5.1.
ARM -Atmospheric Research Model (GCM) Altitudes: -135 км Resolutions: vertical 1 km; longitudinal 1 ; latitudinal 5 time step 5 min. Paramaterizations: Heating - О2, О3, Н2О (Strobel, 1978; Chou et al., 22); IR cooling- СО2, О3, H2O, NО ( Chou et al., 22; Fomichev, 23; Kockarts, 19), GWs (Lindzen, 1981) Planetary waves at lower boundary
Temperature of the atmosphere for July (ARM model runs)
Zonal wind (m/s) for January (ARM model runs) Height 12 1 2 5 3 2 1 5-5 -1-2 -3 - -5 - -7 - - - - -2 2 Latitude
Tidal component (D) in zonal wind for 18 July 12 11 1 9 7 5 3 2 - -7 - -5 - -3-2 -1 1 2 3 5 7 32 3 28 26 24 22 2 18 16 14 12 1 8 6 4 2
Tidal component (SD) in zonal wind for 18 July 12 11 1 9 7 5 3 2 - -7 - -5 - -3-2 -1 1 2 3 5 7 42 38 36 34 32 3 28 26 24 22 2 18 16 14 12 1 8 6 4 2
Longitude 3D Simulation of ozone response (%) at 66 km induced by joint effect of energetic particles (54 hours after - UT 26 October 23) 35 3 25 2 15 1 5 - - - -2 2 Latitude -1-1 -2-3 - -5 - -65
3D Simulation of ozone response (%) at 66 km induced by joint effect of solar protons and precipitating electrons (66 hours after - UT 26 October 23) Longitude 35 3 25 2 15 1 5 - - - -2 2 latitude -1-1 -2-3 - -5 - -7 - -9
3D Simulation of NOy response (ppbv) at 66 km induced by joint effect of solar protons and precipitating electrons (54 hours after - UT 26 October 23) km 35 3 25 2 15 1 5 - - - -2 2 latitude 12 1 2
Longitude 3D Simulation of NOy response (ppbv) at 66 km induced by joint effect of solar protons and precipitating electrons (66 hours after - UT 26 October 23) 35 3 25 2 15 1 5 - - - -2 2 latitude 2 1 1 1 12 1 2 1 1
3D Simulation of ozone response induced by joint effect of solar protons and precipitating electrons (28 October 23 12- UT) Height 2 - - - -2 2 Latitude 15 1 5-5 -1-15 -2-25 -3-35 - -45-5 -55 - -65-7 -75
3D Simulation of ozone response induced by joint effect of solar protons and precipitating electrons (2 November 23 12- UT) 1-5 Heigh(km) 2-1 -15-2 -25-3 -35 - -45 - - - -2 2 Latitude -5
3D Simulation of ozone response induced by joint effect of solar protons and precipitating electrons (4 November 23 12- UT) Height(km) 2 - - - -2 2 Latitude 2 15 1 5-5 -1-15 -2-25 -3-35 - -45-5
3D Simulation of ozone response (%) at 75 N induced by joint effect of solar protons and precipitating electrons after 26 October 23 Height(km) 7 5 3 2-1 -5-1 -15-2 -25-3 -35 - -45 1 1 2 3 5 7 9 1 11 12 Hours -5-55
3D Simulation of ozone response (%) at 75 S induced by joint effect of solar protons and precipitating electrons after 26 October 23 Height(km) 7 5 3 2 1 1 2 3 5 7 9 1 11 12 Hours 2 15 1 5-1 -5-1 -15-2 -25-3 -35 - -45-5 -55 - -65-7 -75 -
Ozone response (%) at 75 N induced by solar protons after 27 October 23 (3D simulation) Height(km) 7 5 3 2 1 2 1 1 1 Hours -5-1 -2-3 - -45
Ozone response (%) at 75 S induced by solar protons after 27 October 23 (3D simulation) Height(km) 7 5 3 2 1 2 1121112-5 -1-2 -3 - -5 - -7 - -9 Hours
Ozone response (%) at 75 N induced by electrons after 27 October 23 (3D simulation) Height(km) 7 5 3 2 1 2 1 1 1 22 Hours -5-1 -2-3 -
Ozone response (%) at 75 S induced by electrons after 27 October 23 (3D simulation) Height(km) 7 5 3 2 1 2 1 1 1 22 Hours -5-1 -2-3 - -5 - -7
Ozone response (%) at 75 N induced by alfa-particles after 27 October 23 (3D simulation) Height(km) 7 5 3 2 1 2 1 12 1 1 1 2 22-5 -1-15 -2-25 -3-32 Hours
Ozone response (%) at 75 S induced by alfa particles after 27 October 23 (3D simulation) Height(km) 7 5 3 2 1 2 1 1 1 22-5 -1-2 -3 - -5 - -7-75 Hours
OH response at 66 km induced by joint effect of energetic particles for 11- UT 28 October 23 (3D model simulations) 35 Longitude 3 25 2 15 1 5 - - - -2 2 Latitude 1.E+7 8.E+6 6.E+6 4.E+6 2.E+6.E+
NOy response (%) induced by energetic particles in 23-24 at 75 N
NOy response (%) induced by energetic particles in 23-24 at 75 S
Simulated O3 response (%) induced by energetic particles in 23-24 at 75 N
O3 response (%) induced by energetic particles in 23-24 at 75 S
Temperature changes induced by SPE 28.1.23 (ARM model runs) 12 1 1 1 Height(km) 2 - - - - -1-1 -1 - - - -2 2 Latitude
Changes in zonal wind (m/s) induced by SPE of 28.1.28 (ARM model runs) Height (km) 12 1 2.6.5.4.3.2.1 -.1 -.2 -.3 -.4 -.5 -.6 -.7 -.8 -.9-1. -1.1-1.2-1.3 - - - -2 2 Latitude
Changes in zonal wind magnitudes (absolute values) induced by SPE for 4 November 12- UT 23 (simulations with ARM) Height(km) 12 1 2 2.2 2. 1.8 1.6 1.4 1.2 1..8.6.4.2 -.2 -.4 -.6 -.8-1. -1.2-1.4-1.6 - - - -2 2 Latitude
Changes in temperature caused by SPE 28.1.23 at 72 km (GCM model runs) Latitude 2-2 - - - 5 1 15 2 25 3 Hours after 26.1.23 [12:].5.4.3.2.1. -.1 -.2 -.3 -.4 -.5 -.6 -.7 -.8 -.9-1. -1.1-1.2
Changes in UV heating (K/day) at 72 km caused by SPE on 1st of November 12- UT 23 (GCM simulations) Latitude 2-2 - - - 5 1 15 2 25 3 Hours after 26.1.23 [12:].4.3.2.1. -.1 -.2 -.3 -.4 -.5 -.6 -.7 -.8 -.9-1. -1.1-1.2-1.3-1.4-1.5
Changes in zonal wind induced by SPE 28.1.23 at 72 km (GCM model runs) Latitude 2-2 - - - 5 1 15 2 25 3 Hours after 26.1.23 [12:] 2. 1.8 1.6 1.4 1.2 1..8.6.4.2. -.2 -.4 -.6 -.8-1. -1.2
GWs input to SPE-induced temperature changes (K/day) on 1.11.23 (GCM model runs) Height [km] 12 1 2 - - - -2 2 1.6 1.5 1.4 1.3 1.2 1.1 1..9.8.7.6.5.4.3.2.1. -.1 -.2 -.3 -.4 -.5 Latitude
Changes in zonal wind (m/s/day) induced by GWs on 1.11.23 as simulated by GCM Height [km] 12 1 2 - - - -2 2 3 3 2 2 1 1-1 -1-2 -2-3 -3-4 -4-5 -5-6 -6 Latitude
Changes in absolute values of zonal wind (m/s) initiated by particle-caused ozone depletion at 75 N (simulations with GCM)
Changes in absolute values of zonal wind (m/s) initiated by particle-caused ozone depletion at 75 S (simulations with GCM)
Changes in absolute values of zonal wind (m/s) initiated by particle-caused ozone depletion at 75 N (simulations with GCM)
Changes in temperature (K) initiated by particle-caused ozone depletion at 75 N (simulations with GCM)
Changes in temperature (K) initiated by particle-caused ozone depletion at 75 N (simulations with GCM)
Changes in temperature (K) on global scale initiated by particle-caused ozone depletion at 37 km (simulations with GCM)
Conclusions 1) Energetic particles change chemical composition and transform wind and temperature fields via ozone depletion. 2) Such effects may induce long-term consequences due to downward transport during polar night. 3) The interaction between chemistry and dynamics leads to temperature effects in the troposphere.
Thank you for your attention!