Energetic Electron Precipitation Research Articles

Global Survey of Energetic Electron Precipitation at Low Earth Orbit Observed by ELFIN

AbstractWe statistically evaluate the global distribution and energy spectrum of electron precipitation at low‐Earth‐orbit, using unprecedented pitch‐angle and energy resolved data from the Electron Losses and Fields INvestigation CubeSats. Our statistical results indicate that during active conditions, the ∼63 keV electron precipitation ratio peaks at L > 6 at midnight, whereas the spatial distribution of precipitating energy flux peaks between the dawn and noon sectors. ∼1 MeV electron precipitation ratio peaks near midnight at L > ∼6 but is enhanced near dusk during active times. The energy spectrum of the precipitation ratio shows reversal points indicating energy dispersion as a function of L shell in both the slot region and at L > ∼6, consistent with hiss‐driven precipitation and current sheet scattering, respectively. Our findings provide accurate quantification of electron precipitation at various energies in a broad region of the Earth's magnetosphere, which is critical for magnetosphere‐ionosphere coupling.

Satellite Observations of the Influence of Energetic Electron Precipitation on the Mesosphere and Stratosphere in the Northern Hemisphere

AbstractThe Earth’s atmosphere is influenced by energetic electrons coming from the magnetosphere. This energetic electron precipitation (EEP) is energized by the solar wind and directly affects in the high‐latitude mesosphere and lower thermosphere (MLT). EEP forms odd nitrogen (NOx) and hydrogen oxides (HOx) which destroy ozone. During winter EEP‐NOx descends to the stratosphere, establishing the indirect EEP effect. Several studies have found that EEP is related to changes in temperature and winds in the northern winter stratosphere. One of the most prominent effects of EEP is the influence on the northern polar vortex, a westerly wind system surrounding the winter pole in the middle atmosphere. Most studies of the EEP effect on dynamical features of the middle atmosphere have relied on either model simulations or reanalysis datasets which are mainly limited to stratospheric heights. We study here EEP effects on chemical and dynamical properties of the stratosphere and mesosphere in the northern hemisphere by using EOS Aura satellite’s measurements of atmospheric properties and POES satellites' measurements of precipitating electrons. We confirm earlier results showing that EEP decreases ozone and affects the temperature in the polar middle atmosphere and strengthens the stratospheric polar vortex. We show that EEP weakens the mesospheric polar vortex in late winter. This effect on polar vortex is partly due to changes in propagation and convergence of planetary waves. Accordingly, the EEP effect on the northern polar vortex depends on planetary waves not only in the stratosphere, as found in earlier studies, but also in the mesosphere.

The high-energy tail of energetic electron precipitation: solar wind drivers and geomagnetic responses

Compositional NOx changes caused by energetic electron precipitation (EEP) at a specific altitude and those co-dependent on vertical transport are referred to as the EEP direct and indirect effect, respectively. The direct effect of EEP at lower mesospheric and upper stratospheric altitudes is linked to the high-energy tail of EEP (≳ 300 keV). The relative importance of this direct effect on NOx, ozone, and atmospheric dynamics remains unresolved due to inadequate particle measurements and scarcity of polar mesospheric NOx observations. An accurate parameterization of the high-energy tail of EEP is, therefore, crucial. This study utilizes EEP flux data from MEPED aboard the POES/Metop satellites from 2004–2014. Data from both hemispheres (55–70° N/S) are combined in daily flux estimates. 164 peaks above the 90th percentile of the ≳ 30 keV flux are identified. These peaks are categorized into absolute E1 and E3 events representing weak and strong ≳ 300 keV responses, respectively. A subset of absolute E1 and E3 events with similar ≳ 30 keV responses is termed overlapping events. Additionally, relative E1 and E3 events are determined by the relative strength of the ≳ 300 keV response, scaled by the initial ≳ 30 keV flux. A comparison between E1 and E3 events aims to identify solar wind and geomagnetic conditions leading to high-energy EEP responses and to gain insight into the conditions that generate a high-energy tail, independent of the initial ≳ 30 keV flux level. Superposed epoch analysis of mesospheric NO density from SOFIE confirms an observable direct impact on lower mesospheric chemistry associated with the absolute E3 events. A probability assessment based on absolute events identifies specific thresholds in the solar wind-magnetosphere coupling function (epsilon) and the geomagnetic indices Kp*10 and Dst, capable of determining the occurrence or exclusion of absolute E1 and E3 events. Elevated solar wind speeds persisting in the recovery phase of a deep Dst trough appear characteristic of overlapping and relative E3 events. This study provides insight into which parameters are important for accurately modeling the high-energy tail of EEP.

Large-scale magnetic field oscillations and their effects on modulating energetic electron precipitation

In this study, we present simultaneous multi-point observations of magnetospheric oscillations on a time scale of tens of minutes (forced-breathing mode) and modulated whistler-mode chorus waves, associated with concurrent energetic electron precipitation observed through enhanced BARREL X-rays. Similar fluctuations are observed in X-ray signatures and the compressional component of magnetic oscillations, spanning from ∼9 to 12 h in MLT and 5 to 11 in L shell. Such magnetospheric oscillations covering an extensive region in the pre-noon sector have been suggested to play a potential role in precipitating energetic electrons by either wave scattering or loss cone modulation, showing a high correlation with the enhancement in X-rays. In this event, the correlation coefficients between chorus waves (smoothed over 8 min), ambient magnetic field oscillations and X-rays are high. We perform an in-depth quasi-linear modeling analysis to evaluate the role of magnetic field oscillations in modulating energetic electron precipitation in the Earth’s magnetosphere through modulating whistler-mode chorus wave amplitude, resonance condition between chorus waves and electrons, as well as loss cone size. Model results further show that the modulation of chorus wave amplitude plays a dominant role in modulating the electron precipitation. However, the effect of the modulation in the resonant energy between chorus waves and energetic electrons due to the background magnetic field oscillations cannot be neglected. The bounce loss cone modulation, affected by the magnetic oscillations, has little influence on the electron precipitation modulation. Our results show that the low frequency magnetospheric oscillations could play a significant role in modulating the electron precipitation through modulating chorus wave intensity and the resonant energy between chorus waves and electron.

Multi-parametric study of seismogenic anomalies during the 2021 Crete earthquake (M=6.0)

It is well established that pre- and co-seismic irregularities in the earth’s atmosphere highly depend on a set of parameters. According to the Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) mechanism, these parameters are associated with various channels (thermal, chemical, acoustic, electromagnetic, etc.) through which an earthquake (EQ) preparation process can have its anomalous phenomena. In this research, we perform a multi-parametric observation using various channels during a strong EQ (M = 6.0) on September 27, 2021, on Crete Island in Greece. We investigate the acoustic and electromagnetic channel using ground and satellite-based observation. We present the Atmospheric Gravity Wave (AGW) in the acoustic channel using the temperature profile computed from the SABER/TIMED instrument. In the electromagnetic channel, ionospheric Total Electron Content (TEC) is recorded by the GNSS-IGS station DYNG in Greece. This TEC information is also used in the acoustic channel anomaly by computing the wave-like structures in the small-scale fluctuation of the TEC profile. We compute energetic (30 to 100 keV) electron precipitation in the inner radiation belt from the NOAA satellite. We also investigate the outcomes of the SWARM satellite to compute the magnetic field and electron density profile. All the parameters show significant seismogenic anomalies (mostly enhancement) before the EQ. To understand each parameter’s temporal and spatial sensitivity, we present a comparison using the anomalies in each parameter.

Electron precipitation caused by intense whistler-mode waves: combined effects of anomalous scattering and phase bunching

Resonant interactions with whistler-mode waves are a crucial mechanism that drives the precipitation of energetic electrons. Using test particle simulations, we investigated the impact of nonlinear interactions of whistler-mode waves on electron precipitation across a broad energy range (10 keV- 1 MeV). Specifically, we focused on the combined effects of conventional phase bunching and anomalous scattering, which includes anomalous trapping and positive bunching. It is shown that anomalous scattering transports electrons away from the loss cone and the only process directly causing precipitation in the nonlinear regime is the phase bunching. We further show that their combined effects result in a precipitation-to-trapped flux ratio lower than the quasilinear expectations in a quasi-equilibrium state. Additionally, we calculated the diffusion and advection coefficients associated with the nonlinear trapping and bunching processes, which are vital for understanding the underlying mechanisms of the precipitation. Based on these coefficients, we characterized the phase bunching boundary, representing the innermost pitch angle boundary where phase bunching can occur. A further analysis revealed that electrons just outside this boundary, rather than near the loss cone, are directly precipitated, while electrons within the boundary are prevented from precipitation due to anomalous scattering. Moreover, we demonstrated that the regime of dominant nonlinear precipitation is determined by the combination of the phase bunching boundary and the inhomogeneity ratio. This comprehensive analysis provides insights into the nonlinear effects of whistler-mode waves on electron precipitation, which are essential for understanding physical processes related to precipitation, such as microbursts, characterized by intense and bursty electron precipitation.

NO Radiative Cooling and Ionospheric Response to the HILDCAA Events Following Geomagnetic Storms

AbstractNitric oxide infrared (5.3 µm) radiative cooling plays an important role in Earth's upper atmospheric energy balance during space weather events. Its behavior during a HILDCAA event following a geomagnetic storm can be more complicated than what is observed in an isolated geomagnetic storm (without HILDCAA). In order to understand this contrast, two moderate (SYM‐Hmin −80 nT; April 2005 (event‐1) and −65 nT; December 2006 (event‐2)), and an intense (Dstmin −150 nT; May 2000 (event‐3)) geomagnetic storms followed by HILDCAA events are considered for this study. It is found that, despite the larger solar wind‐magnetospheric energy coupling during the main phases of event‐1 and event‐2, the NO radiative cooling (NORC) and energetic electron precipitation (EEP) rates are larger during the HILDCAA phase. TIE‐GCM simulations indicate that the NORC pattern can be mainly attributed to variation in temperature, [O], and [NO] during HILDCAA. A larger rates of EEP due to the continuous substorm activity contributed to an additional ionization in D‐region and E‐region of the northern polar ionosphere during HILDCAA period of event‐2. This work is the first attempt to understand the NORC and EEP fluctuations during HILDCAA events preceded by geomagnetic events like event‐1 and event‐2. This study presents a comprehensive investigation of NORC and energetic electron (>27.93 keV) precipitation rate during HILDCAA event preceded by a moderate geomagnetic storm. An atmospheric chemistry‐based NORC emission model is used to understand the variation of NORC during intense geomagnetic storm (event‐3).

Intense Energetic Electron Precipitation Caused by the Self‐Limiting of Space Radiation

AbstractUnderstanding intense electron precipitation is crucial for characterizing radiation belt loss and assessing related impacts on the atmosphere. We investigate the evolution of electron flux during an ensemble of 70 geomagnetic storms, focusing on equatorial and low‐Earth orbit observations of trapped and precipitating ∼30–100 keV energy electrons. We reveal that the most intense electron precipitation is associated with equatorial flux capping through self‐limiting processes, for example, as described theoretically by Kennel and Petschek (1966, https://doi.org/10.1029/jz071i001p00001). Our results indicate that the most intense electron precipitation is caused by electron injections associated with self‐limiting processes. Dawn side injections are observed to have fluxes that exceed the Kennel‐Petschek limit, consistent with the excitation of strong chorus waves and resulting in intense precipitation and return of the trapped flux to the Kennel‐Petschek limit. Our results clearly demonstrate the important role of self‐limiting processes in affecting the dynamics of newly injected electrons and driving intense electron precipitation.

Ionospheric Plasma Density Gradients Associated With Night‐Side Energetic Electron Precipitation

AbstractEnergetic electron precipitation from the equatorial magnetosphere into the atmosphere plays an important role in magnetosphere‐ionosphere coupling: precipitating electrons alter ionospheric properties, whereas ionospheric outflows modify equatorial plasma conditions affecting electromagnetic wave generation and energetic electron scattering. However, ionospheric measurements cannot be directly related to wave and energetic electron properties measured by high‐altitude, near‐equatorial spacecraft, due to large mapping uncertainties. We aim to resolve this by projecting low‐altitude measurements of energetic electron precipitation by ELFIN CubeSats onto total electron content (TEC) maps serving as a proxy for ionospheric density structures. We examine three types of precipitation on the nightside: precipitation of <200 keV electrons in the plasma sheet, bursty precipitation of <500 keV electrons by whistler‐mode waves, and relativistic (>500 keV) electron precipitation by EMIC waves. All three types of precipitation show distinct features in TEC horizontal gradients, and we discuss possible implications of these features.

Energetic Electron Precipitation During Slot Region Filling Events

AbstractThe slot region marks the equatorward boundary of the energetic electron precipitation (EEP). There are, however, numerous reports where energetic electrons cross these boundaries and fill the slot region. The ensuing EEP will occur long after the geomagnetic activity subsides. This is a missing energy input in current EEP estimates scaled by geomagnetic indices. This study explores the occurrence rate, duration, and local time dependence of slot region filling events using observations from the National Oceanic and Atmospheric Administration/Polar Orbiting Environmental Satellites over a full solar cycle from 2004 to 2014. The EEP flux estimates are based on the Medium Energy Proton Electron Detector 0° and 90° detectors and the theory of pitch angle diffusion by wave‐particle interaction. The occurrence rates of >43, >114, and >292 keV events are found to be strongly energy and solar cycle dependent. Higher energy events are more likely to be associated with Coronal Mass Ejections and stronger geomagnetic deflections compared to lower energy events. Solar wind speed, Bz, and Ey reveal a calm period before the events, potentially important for preconditioning the ensuing magnetospheric mass convection. The slot region reforms more efficiently closer to the plasmapause, which creates a double EEP band throughout the recovery period. The slot region EEP maximizes around noon throughout the afternoon/evening sector, consistent with pitch angle scattering from plasmaspheric hiss and lightning induced whistler mode waves. Concurrent with slot region filling events, the Michelson Interferometer for Passive Atmospheric Sounding/Envisat nitric oxide density show an increase at <55° corrected geomagnetic latitudes. This demonstrates the importance of including slot region EEP when assessing the EEP impact on the atmosphere.

Conjugate Observation of Whistler Mode Chorus, ECH Waves and Dayside Diffuse Aurora by MMS and Ground‐Based Yellow River Station

AbstractDayside diffuse auroras (DDAs) originate from the precipitation of energetic electrons in the central plasma sheet, induced by wave‐particle interactions. This study examines a specific case characterized by the detection of whistler mode chorus and electron cyclotron harmonic (ECH) waves using the Magnetospheric Multiscale (MMS) satellites, along with the conjugate observation of DDA at the Yellow River (YR) station. ECH waves are observed shortly after the satellite enters the magnetosphere from the magnetosheath and persist throughout the magnetosphere region. Chorus waves are also observed within the magnetosphere, and probably generated near the B minima, followed by the propagation along magnetic field lines. Analysis of MMS particle data reveals an increase in low‐energy electrons and scattering of high‐energy electrons during periods of enhanced whistler mode chorus waves. These precipitated electrons manifest as DDA phenomena. The observed evidence demonstrates the simultaneous occurrence of ECH, chorus waves, electron pitch angle scattering, and DDA. Moreover, the concurrent presence of these waves in the outer magnetosphere during the dayside raises important considerations regarding the relative roles of different plasma waves, particle energization, precipitation dynamics, and the overall development of DDA phenomena.

An Updated Geomagnetic Index‐Based Model for Determining the Latitudinal Extent of Energetic Electron Precipitation

AbstractEnergetic Electron Precipitation (EEP) from the Earth's plasma sheet and the radiation belts is an important feature of atmospheric dynamics through their destruction of ozone in the lower thermosphere and mesosphere. Therefore, understanding the magnitude of the atmospheric impact of the Sun‐Earth interaction requires a comprehensive understanding of the intensity and location of EEP. This study improves the accuracy of a previous pressure‐corrected Dst model that predicts the equatorward extent of >43, >114, and >292 keV EEP using the measurements from the Medium Energy Proton Electron Detector detector of six National Oceanic and Atmospheric Administration/Polar Orbiting Environmental Satellites and EUMETSAT/METOP satellites. The improvement is achieved through multiple linear regression of pressure‐corrected Dst and pressure‐corrected Ring Current (RC) indices. The RC index mitigates the baseline variation of the Dst index that created an inherent solar cycle bias in the previous model. The new model is then extended to the Southern Hemisphere (SH) after removing the South Atlantic Anomaly longitudes from the data. More than 80% of the residuals lie within ±1.8° Corrected Geomagnetic Latitude (CGMLat) in the Northern Hemisphere and within ±1.98° CGMLat in the SH.

Estimation of Energy Spectrum of Precipitating Magnetospheric Electrons Based on Bremsstrahlung X‐Ray Fluxes Recorded in the Atmosphere

AbstractNumerous energetic electron precipitation events were recorded since 1963 in the atmosphere at polar latitudes in the course of regular measurements of charged particle fluxes in the Earth's atmosphere being performed by the Lebedev Physical Institute. The experimental data obtained represent the only uniform database in the world on electron precipitation events recorded directly in the atmosphere well below satellite orbits. The precipitating electrons are absorbed in the upper atmosphere. However, they generate X‐rays that can penetrate deep into the atmosphere sometimes down to altitudes of 20–35 km accessible to balloon measurements. These experimental data allow studying energy, temporal and spatial characteristics of electron precipitation events. In particular, based on PLANETOCOSMICS/GEANT4 we developed a method for evaluating the energy spectra of the primary flux of precipitating electrons. This method was used for evaluation of primary spectra of precipitating electrons assuming exponential energy distribution of incident electrons. Now, for the development of the method, we present the possibility of determining the energy spectra of electrons in the power‐law form, using new software RUSCOSMICS code. This new method, based on the GEANT4 software, makes it possible to describe the transport of precipitating electrons in the atmosphere taking into account the evolution of the energy and pitch‐angle distributions of electrons and X‐ray photons.

External and Internal Causes of the Stormtime Intensification of the Dawnside Westward Auroral Electrojet

AbstractThe dawn‐dusk asymmetry of magnetic depression is a characteristic feature of the storm main phase. Recently Ohtani (2021, https://doi.org/10.1029/2021JA029643) reported that its magnitude is correlated with the dawnside westward auroral electrojet (AEJ) intensity, and suggested that the dawnside AEJ intensification is a fundamental process of the stormtime magnetosphere‐ionosphere coupling. In this study we observationally address the cause of the dawnside AEJ intensification in terms of four scenarios. That is, the dawnside AEJ intensifies because (a) the external driving of global convection strengthens, (b) solar wind compression enhances energetic electron precipitation, and therefore, ionospheric conductance, through wave‐particle interaction, (c) the substorm current wedge forms in the dawn sector, and (d) energetic electrons injected by nightside substorms drift dawnward, and the subsequent precipitation enhances ionospheric conductance. We find an event that fits each scenario, and therefore, none of these scenarios can be precluded. However, the result of a superposed epoch analysis shows that some causes are more prevalent than others. More specifically, (a) although the enhancement of external driving may precondition the dawnside AEJ intensification, it is rarely the direct cause; (b) external compression probably explains only a small fraction of the events; (c) prior to the dawnside AEJ intensification, the westward AEJ tends to intensify in the midnight sector along with mid‐latitude positive bays, which suggests that the substorm injection of energetic electrons is the most prevalent cause. This last result may also be explained by the dawnside expansion of the substorm current wedge, which, however, is arguably far less common.

Observed Loss of Polar Mesospheric Ozone Following Substorm‐Driven Electron Precipitation

AbstractSeveral drivers cause precipitation of energetic electrons into the atmosphere. While some of these drivers are accounted for in proxies of energetic electron precipitation (EEP) used in atmosphere and climate models, it is unclear to what extent the proxies capture substorm‐induced EEP. The energies of these electrons allow them to reach altitudes between 55 and 95 km. EEP‐driven enhanced ionization is known to result in production of HOx and NOx, which catalytically destroy ozone. Substorm‐driven ozone loss has previously been simulated, but has not been observed before. We use mesospheric ozone observations from the Microwave Limb Sounder and Global Ozone Monitoring by Occultation of Stars instruments, to investigate the loss of ozone during substorms. Following substorm onset, we find reductions of polar mesospheric (∼76 km) ozone by up to 21% on average. This is the first observational evidence demonstrating the importance of substorms on the ozone balance within the polar atmosphere.

The Global Variation of Low Ionosphere Under Action of Energetic Electron Precipitation

AbstractThe density in the D region is hard to observe on a large scale, because the density in this region is too low. In addition to solar radiation changes with a sinusoidal period of 24 hr, the energetic electron precipitation (EEP) is the major source for the ionospheric D region. In this study, we utilized observations from NOAA/MetOp satellites to derive the electron density profile from 60 to 150 km under the influence of EEP. It seems that the derived electron profile is roughly consistent with the observed one. Additionally, using above derived method of low ionospheric (D and E regions) electron profile based on the observations of NOAA/MetOp satellites, the actions of EEP on global ionosphere under different level of geomagnetic activity are analyzed. We find that EEP has a significant impact on the electron density in the altitude range from 60 to 90 km in the L ∼ 4.5–8 and MLT ∼ 9–14 regions, particularly during strong geomagnetic activity. A peak in Total Electron Content (TEC) was observed at MLT ∼ 11–12 and L ∼ 5.5–6, exhibiting an enhancement approximately 19.96 times higher than that observed during quiet periods. However, the effects of EEP decrease as altitude increases. The TEC enhancement at altitudes from 90 to 120 km during strong geomagnetic activity is like that at 60–90 km, but only about 2.08 times higher than that during quiet periods. In addition, only a weak increase is detected around midnight at altitudes from 120 to 150 km.

Surface‐Plasma Interactions at Europa in Draped Magnetospheric Fields: The Contribution of Energetic Electrons to Energy Deposition and Sputtering

AbstractWe calculate the time‐varying spatial distribution of energetic magnetospheric electron influx onto Europa's surface by combining a hybrid model of the moon's draped electromagnetic environment with a relativistic particle tracer. We generate maps of the energetic electron influx patterns at four distinct locations of Europa relative to the center of the Jovian magnetospheric current sheet. For a full synodic rotation of Jupiter, these results are applied to constrain the averaged number and energy influx patterns as well as the O2 sputtering rates associated with energetic electron precipitation. We also determine the relative contributions of magnetospheric ions and electrons to surface erosion and exospheric genesis at Europa. Our major results are: (a) Except for a small region near Europa's downstream apex, the moon's entire surface is exposed to heavy irradiation by magnetospheric electrons. (b) The spatial distribution of energetic electron influx onto Europa's surface is only slightly modified by field line draping and the induced magnetic field from the moon's subsurface ocean. (c) The contributions of magnetospheric electron and ion impacts to energy deposition onto Europa's surface are of the same order of magnitude. (d) Within uncertainties, impinging magnetospheric electrons and ions make similar contributions to O2 sputtering from Europa's surface. (e) The spatial distribution of electron energy influx and the observed concentrations of sulfuric acid (H2SO4) are only weakly correlated, suggesting that energy deposition by magnetospheric electron impacts is not a necessary agent for H2SO4 production within Europa's surface.

Investigation of the drivers and atmospheric impacts of energetic electron precipitation

Investigation of the drivers and atmospheric impacts of energetic electron precipitation

A Ground‐Based Instrument Suite for Integrated High‐Time Resolution Measurements of Pulsating Aurora With Arase

AbstractA specialized ground‐based system has been developed for simultaneous observations of pulsating aurora (PsA) and related magnetospheric phenomena with the Arase satellite. The instrument suite is composed of (a) six 100 Hz sampling high‐speed all‐sky imagers (ASIs), (b) two 10 Hz sampling monochromatic ASIs observing 427.8 and 844.6 nm auroral emissions, (c) a 20 Hz sampling fluxgate magnetometer. The 100 Hz ASIs were deployed in four stations in Scandinavia and two stations in Alaska, which have been used for capturing the main pulsations and quasi 3 Hz internal modulations of PsA at the same time. The 10 Hz sampling monochromatic ASIs have been operative in Tromsø, Norway with the 20 Hz sampling magnetometer. Combination of these multiple instruments with the European Incoherent SCATter (EISCAT) radar enables us to detect the low‐altitude ionization due to energetic electron precipitation during PsA and further to reveal the ionospheric electrodynamics behind PsA. Since the launch of the Arase satellite, the data from these instruments have been examined in comparison with the wave and particle data from the satellite in the magnetosphere. In the future, the system can be utilized not only for studies of PsA but also for other classes of aurora in close collaboration with the planned EISCAT_3D project.

Mesospheric Ozone Depletion Depending on Different Levels of Geomagnetic Disturbances and Seasons

Energetic electron precipitation (EEP) into the atmosphere are considered to play an important role in the natural forcing of the ozone variability and dynamics of the middle atmosphere during magnetospheric and geomagnetic disturbances. Energetic electrons from the radiation belt spill out into the atmosphere during geomagnetic disturbances and cause additional ionization rates in the polar middle atmosphere. These rates of induced atmospheric ionization lead to the formation of radicals in ion-molecular reactions at the heights of the mesosphere with the formation of reactive compounds of odd nitrogen groups NOy and odd hydrogen groups HOx. These compounds are involved in catalytic reactions that destroy ozone. The percentage of ozone destruction can depend not only intensity of EEP but also on season where it happens. In this work, we study mesospheric ozone depletion depending on seasons and precipitating energetic electrons with energies from keV up to relativistic energies about 1 MeV, based on the NOAA POES satellites observations in 2003. For estimation ozone deplation we use a one-dimensional radiative-convective model with ion chemistry. As one of the main results, we show that, despite the intensity of EEP-induced ionization rates, polar mesospheric ozone cannot be destroyed by EEP in summer in the presence of UV radiation. In winter time, the maximum ozone depletion, at altitude of about 80 km, can reach up to 80% during strong geomagnetic disturbances. In fall and spring, the maximum ozone depletion is less intense and can reach 20% during strong geomagnetic disturbances. Linear relation of EEP induced maximum mesospheric ozone depletion depending on geomagnetic disturbances and seasons have been obtained.