The Predictive Capabilities of the Auroral Electrojet Index for Medium Energy Electron Precipitation

The chemical imprint of the energetic electron precipitation on the atmosphere is now acknowledged as a part of the natural forcing of the climate system. It has, however, been questioned to which degree current proxies are able to quantify the medium energy electron (MEE) (≳30 keV) precipitation and the associated daily and decadal variability. It is particularly challenging to model the high energy tail (≳300 keV) of MEE, both in terms of the intensity as well as the timing. This study explores the predictive capabilities of the AE index for the MEE precipitation. MEE measurements from the NOAA/POES over a full solar cycle from 2004 to 2014 are applied. We combine observations from the MEPED 0° and 90° detectors together with theory of pitch angle diffusion by wave-particle interaction to estimate the precipitating fluxes. To explore the energy dependent time scales, each of the MEPED energy channels, > 43, >114, and >292 keV are evaluated independently. While there is a strong correlation between the daily resolved AE index and >43 keV fluxes, it is a poor predictor for the >292 keV fluxes. We create new AE based MEE proxies by accumulating the AE activity over multiple days, including terms counting for the associated lifetimes. The results indicate that AE based proxies can predict at least 70% of the observed MEE precipitation variance at all energies. The potential link between the AE index, substorms and the MEE precipitation is discussed.


INTRODUCTION
Precipitating auroral electrons (≲ 30 keV) and protons (≲ 1 MeV) from the plasma sheet will ionize the lower thermosphere and upper mesosphere. Medium energy electrons (MEE) (≳ 30 keV) from the radiation belts have sufficient energy to penetrate deep into the mesosphere (Turunen et al., 2009). The associated ionization enhances the production rate of NOx and HOx gasses, which in turn can reduce ozone in catalytic processes (e.g., Saetre et al., 2004;Andersson et al., 2012;Sinnhuber et al., 2016;Zawedde et al., 2016). During polar winter, the chemical impact can be long lasting and influence temperatures, winds, and wave propagation. This chain of reactions can impact the strength of the stratospheric polar vortex, causing a dynamical signal that may propagate all the way down to the surface (Seppälä et al., 2013;Maliniemi et al., 2016). To account for this natural solar forcing of the atmosphere, a parametrization of energetic electron precipitation isfor the first timepart of the official input to the Coupled Model Intercomparison Project phase 6 (CMIP 6) going into the Intergovernmental Panel on Climate Change sixth assessment report (Matthes et al., 2017). The MEE ionization rate dataset therein is based on observations from the Medium Energy Proton and Electron Detector (MEPED) instrument on board the NOAA/Polar Orbiting Environmental Satellites (POES), and the geomagnetic Ap index is used as a proxy to provide an extended time series beyond the satellite observation period (van de Kamp et al., 2016). There is, however, an active discussion to what extent this approach gives a representative flux and ionization rate level (Mironova et al., 2019;Nesse Tyssøy et al., 2019;Pettit et al., 2019;Clilverd et al., 2020). The CMIP6 flux is a general underestimate, largely ascribed to the use of the vertical (0°) detector on MEPED which only covers a small fraction of the loss cone (Nesse Tyssøy et al., 2019). Recent studies, however, also point out that the role of substorms in driving MEE precipitation is not readily explained by a single magnetic index value. Hence, substorms is one of the main unknowns in the existing proxies when considering the MEE precipitation and especially its high energy tail (≳ 300 keV) (Partamies et al., 2021).
Several processes in the magnetosphere contribute to the driving of the MEE precipitation. A globally induced electric field during southward interplanetary magnetic field accelerates and transports electrons from the magnetotail into the inner magnetosphere where they become a part of the radiation belts. In parallel, localized transient-induced electric fields, created from the magnetotail collapse during substorm activity, also energize and increase the MEE population in the radiation belts. So-called seed particles of 10-100 s keV can be directly injected into the inner magnetosphere during a substorm (Li et al., 2009;Jaynes et al., 2015). Further acceleration also occurs as the injection of source particles (tens of keV) gives rise to Very Low Frequency (VLF) wave growth, which may resonantly interact and accelerate radiation belt electrons to MEE (Borovsky & Yakymenko, 2017). Repeated substorm activity has been shown to be particularly important for MEE fluxes Partamies et al., 2021), possibly due the induction electric field directly energizing the trapped electrons as they drift across the nightside during substorm expansion phases (Dai et al., 2014). Partamies et al. (2021) identified the substorm occurrence by a regional AE index and used cosmic noise absorption enhancement as a measure of the MEE precipitation. They showed that for multi-night substorm events, the first night was rarely associated with the most intense absorption. Instead, the high-energy electron population, needed to cause the strongest absorption, was built up over one to two additional nights of substorm activity. This was further confirmed by MEPED in situ particle spectra. In the expansion phases the bulk of the spectra showed a local maximum flux in the range of a few keV to 10 keV, while in the recovery phases higher fluxes were seen in the range of tens of keV to hundreds of keV. Based on the SuperMAG substorm event list, Rodger et al. (2016) constructed a superposed epoch analysis differing between isolated and recurrent events. Even though their main focus was the trapped MEE fluxes, they also showed higher precipitating MEE fluxes in the epoch analysis for the recurrent events compared to isolated events, in particular with respect to the high energy tail (>300 keV). Although not commented on, Figure A1 in Rodger et al. (2016) showed an instant effect as well as a gradual build-up over the consecutive days. Seppälä et al. (2015) used the substorm model developed by Beharrell et al. (2015) to investigate the production of HOx and NOx and the subsequent depletion of mesospheric ozone with the Sondakylä Ion and Neutral Chemistry model (Turunen et al., 2009). The peak loss of mesospheric ozone was observed during the third and fourth day in the period of the repetitive substorm activity. This implies that the creation of MEE precipitation involves processes with various time constants and potentially accumulated effects, where there are increasing delays in flux buildup with energy (Boynton et al., 2016;Ødegaard et al., 2017;Stepanov et al., 2021).
Out of the International Association of Geomagnetism and Aeronomy (IAGA) recognized magnetic indices, the Auroral Electrojet (AE) index best corresponds to substorm activity. The AE index goes back to the 1960s and is constructed from the horizontal magnetic field component recorded with 1-min time resolution at 10-13 magnetic observatories located under the average auroral oval in the Northern Hemisphere (NH) (geomagnetic latitudes 60°-70°) (Davis & Sugiura, 1966;Kauristie et al., 2017, and references therein). The upper and lower envelope curves are defined to be the AU and AL indices which characterize the intensity of eastward and westward electrojets, respectively. The difference, AU-AL, defines the AE-index (Davis & Sugiura, 1966). The magnetotail energy release associated with the substorm expansion phase affects the intensity and spatial distribution of electric currents in the auroral oval region. In particular, the substorm current wedge causes sudden enhancements in the westward electrojet (Kauristie et al., 2017).
In this study, we explore the prediction capability of AE in regard to MEE precipitating fluxes. We combine observations from both the MEPED 0°and 90°detectors together with the theory of pitch angle diffusion by wave-particle interaction to quantify the MEE flux in the bounce loss cone. We treat the MEPED energy channels, > 43, >114, and >292 keV, independently to explore their associated time delays in respect to the geomagnetic activity. The objective is to determine the potential of AE as a proxy for MEE precipitating fluxes with a particular focus on the high energy tail. The paper is organised as follows: Section 2 describes the MEPED detectors and the methods applied to estimate the loss cone fluxes. It provides a short introduction to the AE index followed by a comparison between the AE index and the SuperMAG substorms list from the period 2004-2014. Section 3 starts with a simple inspection of the correlation between the AE index and the MEE fluxes, where the investigation of time delays and build-up effects points toward a more advanced model. Section 4 provides a discussion of the implication of the results and the potential role of substorms as the physical process linking the AE based proxies and the MEE fluxes.
Frontiers in Astronomy and Space Sciences | www.frontiersin.org October 2021 | Volume 8 | Article 714146 2 2 DATA AND METHODS

The Medium Energy Electrons Loss Cone Fluxes
The MEPED instrument is mounted on the NOAA/POES series and three EUMETSAT/MetOp spacecraft (Evans & Greer, 2004). The satellites are Sun-synchronous, low-altitude (∼ 850 km), polar orbiting spacecrafts. Their orbital period is about 100 min, resulting in 14-15 orbits for each satellite each day. The combined measurements offer a long, near continuous observation of MEE from 1979 until today. During the latest decades a constellation of up to six operating satellites has allowed for a more global magnetic local time coverage.
The MEPED instrument consists of two directional electron telescopes and two directional proton telescopes, as well as an omni-directional detector for very energetic protons measured over a wide range of angles (Evans & Greer, 2004). The field of view of both the 0°and 90°telescopes is 30°full width. The nominal energy limits of the MEPED telescopes are given as > 30, >100, and >300 keV. The true detector efficiency, however, will depend on the incoming energy spectrum (Yando et al., 2011). Ødegaard et al. (2017) utilize the geometric factors given in Yando et al. (2011) to determine new optimized effective integral energy limits >43, >114, and >292 keV and associated geometric factors based on a series of realistic power law and exponential spectra. Furthermore, the spurious response to contaminating protons is accounted for. The observed proton fluxes are first corrected for degradation due to radiation damage by applying correction factors derived by Sandanger et al. (2015) and Ødegaard et al. (2016). Subsequently, a monotonic piecewise cubic Hermite interpolating polynomial is applied to the corrected proton fluxes. The proton flux in the energy ranges known to impact the respective electron channels (Evans & Greer, 2004), are then retrieved and subtracted from the original measured electron fluxes.
A detailed discussion on which radiation belt populations the 0°and 90°telescopes measure has been presented in Appendix A in Rodger et al. (2010). In general, it shows that at middle and high latitudes the 0°telescopes measure particle fluxes that will be lost to the atmosphere, whereas the 90°telescopes detect precipitating particle fluxes and/or trapped particles in the radiation belts. This implies that in the frequent case of pitch angle anisotropy the 0°detector will underestimate, while the 90°d etector will overestimate the flux of precipitating electrons Nesse Tyssøy et al., 2019). A more realistic estimate can be achieved by combining fluxes from both the 0°and 90°telescopes together with electron pitch angle distributions from theory of wave-particle interactions in the magnetosphere. We solve the Focker-Planck equation for particle diffusion (Kennel & Petschek, 1966;Theodoridis & Paolini, 1967) for a wide range of diffusion coefficients (The specific equations are also given in Nesse Tyssøy et al. (2016)). The solutions are then transformed to the satellite altitude and saved in a look-up table. When comparing the theoretical pitch angle distributions with the measured particle fluxes the procedure is as follows: • Determine the pitch angles of center look directions of the 0°a nd 90°telescopes. • Calculate the ratio between of the fluxes detected by the 0°a nd 90°detector. • Calculate the ratio between the fluxes for the theoretical pitch angle distributions taking into account the look directions. • Determine which of the theoretical pitch angle distributions best corresponds to the observed ratio.
Finally, the size of the loss cone, predicted based on the International Geomagnetic Reference Field (IGRF) model, are applied to estimate the precipitating fluxes. The loss cone flux estimate is done separately for each energy channel as the level of particle diffusion will vary with energy. A detail explanation of the method can be found in Nesse Tyssøy et al. (2016).
We use MEE precipitation estimates from a full solar cycle from 2004 to 2014. This includes the active years in the declining phase of cycle 23, the extreme minimum, and the moderate maximum of cycle 24. Figure 1 illustrates the number of NOAA and MetOp satellites applied throughout the 11 years alongside the MLT coverage in the NH. Although, the combination of satellites has a near global coverage on a daily scale, it is not equally distributed in terms of MLT. In particular, the evening sector and midnight sector are poorly represented in the NH. The daily average of the fluxes is therefore first calculated for four separate MLT sectors, 0-6 MLT, 6-12 MLT, 12-18 MLT, and 18-24 MLT over the CGM latitude band 55°-70°. Thereafter, a daily MEE flux is achieved as the average of the four MLT regions to ensure that the MEE flux is weighted equally in respect to MLT. Figure 2 shows the resulting daily integral fluxes for the energies >43 keV (blue line), >114 keV (red line), and >292 keV (black line) over the full solar cycle. The values are given as the logarithm of the flux value, illustrating that the flux of the high energy tail >292 keV are typically two order of magnitude less than the >43 keV fluxes.

The Auroral Electrojet Index
The AE-index is designed by Davis and Sugiura (1966) to monitor the electrojet activity. The link to substorm activity was assumed from the very beginning as the substorm current system lies in the ionosphere. However, its deficiency with respect to coarse geographic latitude and longitude coverage was quickly pointed out. The observations, limited to a geomagnetic latitude band of 60°-70°N, could not always detect the dynamic auroral oval, both expanding equatorward and contracting poleward of the 12-13 stations. Short-term and localized events in the midnight sector, such as substorms, can easily be missed by the coarse network of stations.
Besides substorms, pseudo-breakups, steady magnetospheric convection events, sawtooth injections, poleward boundary intensifications, or a mixture of these modes have been recognized in the AE data (McPherron, 2015). To overcome some of these challenges, regional electrojet-indices have been created (Tanskanen, 2009), and multiple regional magnetometer chains have been combined to compile globally denser network of stations (Gjerloev, 2012). Analysis have been performed to identify individual substorms. However, no generally accepted Frontiers in Astronomy and Space Sciences | www.frontiersin.org October 2021 | Volume 8 | Article 714146 method exists to identify substorm events from auroral electrojet indices, and different sets of criteria are used in different studies (e.g., Tanskanen et al., 2002;Newell & Gjerloev, 2011). Figure 3 shows the daily averaged AE index alongside the daily number of substorms identified by Newell and Gjerloev (2011) from 2004 to 2014. The Pearson correlation coefficient, r, between the two datasets over the entire solar cycle is 0.89, giving r 2 0.79. It should, however, be emphasized that using the AE index on a daily scale makes identification of substorm/non-substorm modes problematic in particular for long active periods. Therefore, Figure 3 shows a potential (but not proven) link between the daily averaged AE index and substorms.
We note that the daily AE index has a pronounced seasonal bias where the AU and AL indices maximize during summer and equinoctial months, respectively (Ahn et al., 2000). The equinox bias is due to the seasonal bias of solar wind driving which also applies to MEE. While the summer maxima, clearly evident in Figure 3, is due to increased background ionization from UV. This implies that the relation between the AE index and other parameters such as the number of substorms or MEE flux will vary with season. It also means that if AE, based on observations only from NH, is to be used as a global proxy for both hemispheres the seasonal bias needs to be addressed. In the MEE proxy developed in this study, we remove the seasonal trend by subtracting the minimum daily AE value found in a moving window of ±14 days from the daily resolved AE index.

RESULTS
The daily AE index vary by three orders of magnitude over the 11 years. The daily fluxes of >43, >114, and >292 keV loss cone    Table 1 lists the Pearson correlation coefficient between the daily resolved logarithmic value of the AE index and the logarithmic value of >43, >114, and >292 keV The square value of the correlation coefficients of 0.88 and 0.73 implies that the AE regression models fit with 77% and 53% of the variability of the >43 keV and >114 keV fluxes, respectively. The correlation coefficient decreases with energy, where a value of 0.47 implies that the AE regression model can only fit 22% of the variability of the >292 keV fluxes. This might reflect that a larger fraction of the >43 keV and >114 keV electron fluxes are directly injected during the substorm activity, while relatively fewer electrons >292 keV are part of the initial seed population. Table 1 also shows that the AE index is best correlated with >43 keV and >114 keV electron fluxes in the post-midnight MLT sector which supports the link to substorm nightside injection and the subsequent eastward electron drift around the Earth and westward electrojet enhancements. The MLT bias is, however, not as prominent for >292 keV electron fluxes. Figure 5A shows the Pearson correlation coefficients between the logarithmic value of the AE index 0-8 days prior to the logarithmic value of the observed MEE fluxes. The highest correlation is found on day zero, zero, and two for >43, >114, and >292 keV fluxes, respectively. The increasing offset as a function of energy substantiates that time is a prerequisite for the MEE high energy tail. The broad correlation peak and gradual decay further suggest that the MEE fluxes are influenced by the geomagnetic activity level of several days. Figure 5B shows the correlation between the logarithmic value of the AE index accumulated over consecutive longer periods prior to the logarithmic value of the observed MEE fluxes. The correlation coefficients between the accumulated AE index and the >43 keV fluxes peaks at 0.88 when both the zero and first preceding day are taken into account, which based on r-squared corresponds to approximately 77% of the flux variability. The correlation coefficient of the AE index for the >114 keV fluxes reaches about 0.83 when 3-4 days are accumulated, increasing the predictive capability of AE index from about 53% to 69%. The correlation coefficient between the accumulated AE index and the >292 keV fluxes increases rapidly the first few days and as many as seven preceding days are required for it to reach its peak of 0.81. Now, the predictive capability of the AE index has increased from about 22% to 66%.
Based on the correlation coefficients shown in Figure 5A is unlikely that each of the preceding days are equally important as assumed in this simple model shown in Figure 5B. As time evolves the seed population injected during a specific substorm is not only accelerated, but it is likely to be lost to the atmosphere or to the magnetopause. A more realistic model, where the AE index is weighted in respect to lifetimes of the particle population, could increase the predictive capability of the AE index even more. We assume the following weighting of geomagnetic activity:

MEE proxy
where i refers to the preceding days and τ is the assumed lifetime. Due to the broad peak found in Figure 5 for the >292 keV fluxes, the equation includes ten offset days. Furthermore, Figure 3 shows a clear seasonal trend in the daily AE index which could impact the correlation with the MEE fluxes. As described in Section 2.2, we remove the seasonal trend in the AE index by subtracting the minimum daily AE value found in a moving window of ±14 days from all AE values. Now, the maximum correlation coefficients increases to 0.91, 0.89, and 0.84, using the optimized lifetime, τ, of one, three, and 9 days for >43 keV, > 114 keV, and >292 keV, respectively. For the >43 keV fluxes a correlation coefficient of 0.91 is a small improvement from the 0.88 found in Figure 5. In addition, to better model low and high extremes, we fit two separate linear equations based on the logarithmic value of the MEE proxy . The optimized separating boundary and the associated model equations are achieved by stepwise moving the limit over the     Frontiers in Astronomy and Space Sciences | www.frontiersin.org October 2021 | Volume 8 | Article 714146 7 entire MEE proxy interval while performing two independent linear regression fits for the fluxes associated with the MEE proxy values below and above the limit. The boundary giving the highest correlation coefficient and smallest RMSE for the two model fits combined are selected. An overview of the AE based MEE proxy , including the mathematical expressions, lifetime τ, and linear equations are given in Table 2 (Note that the limits are based on the MEE proxy and not the daily AE values.) These simple measures increase the correlation coefficients to 0.93, 0.90, and 0.85 for >43 keV, > 114 keV, and >292 keV, respectively. This implies that about 86%, 81% and 72% of the daily flux variability, considering all MLTs and the full solar cycle, can be accounted for by the MEE proxy models.
The upper panel in Figure 6 shows the resulting linear fits to the AE-based MEE proxy for >43 keV electron fluxes. The model captures the day-to-day fluctuations. The lower panel shows the differences between the observed and modelled fluxes. Similarly to the flux, the error are given as log [J (s −1 cm 2 sr −1 )], which implies that values larger than one would correspond to one order of magnitude difference. The 5/95 percentile black, dashed lines demonstrate that the typical error is less than 0.4 log [J (s −1 cm 2 sr −1 )], corresponding to the value 10 0.40 ∼ 2.5. This means that for 90% of the days the model predicts fluxes that deviates from the observed fluxes by less than a factor of 2.5. The largest error found is 0.99 which correspond to a factor of 9.8. Hence, all values are within one order of magnitude of the observed fluxes. We note, however, specific periods in time where the model appears to have a bias. In the declining phase the model underestimates the flux level, while the opposite occurs for the extreme solar minimum year of 2009. There also appear to be a seasonal bias where the model overestimate/underestimate the flux values during summer/ winter. Figure 7 shows the scatter plot of the MEE proxy model vs the observed >43 keV electron flux. It illustrates that the largest errors are found during low to moderate activity. Compared to Figure 4 the improvement using the MEE proxy model on a detrended AE index is readily evident. Figure 8 shows the resulting linear fits to the AE-based MEE proxy for >114 keV electron fluxes. The accuracy and weaknesses of the modelled >114 keV fluxes are similar to the modelled >43 keV fluxes as shown in Figure 8. However, Figure 9 shows a tendency of the MEE proxy model to overestimate the fluxes during high activity. Figure 10 and Figure 11 confirm a similar trend for the modelled >292 keV fluxes. All correlation coefficients, including the MLT dependence are listed in Table 3.

DISCUSSION AND CONCLUSION
There is increasing evidence that substorms are key in driving MEE precipitation Partamies et al., 2021). Nevertheless, the substorms remain partly unresolved in the existing proxies when considering the MEE precipitation (van de Kamp et al., 2016). In this study, we have explored the prediction capability of AE in regard to MEE precipitation on a daily scale over a full solar cycle. The potential link between the AE index, substorms, and MEE precipitation will be discussed, alongside the progressive time delay between the geomagnetic activity and the MEE precipitation.
The initial correlation study, summarized in Figure 5, reveals a high coherence between the daily AE index with zero lag and >43 keV precipitating electron fluxes. Similarly, Figure 3 suggests a strong coherence with the daily AE index and the daily number of substorms. The magnetotail dipolarization during substorms will directly inject source and seed particles in this energy range, some of which will precipitate into the atmosphere. It is therefore a realistic physical link between the AE intensity and >43 keV fluxes, despite the fact that the >43 keV fluxes deposit their energy below 100 km and do not directly contribute to the intensity of the electrojets. The strong correlation between the AE index and >43 keV fluxes is in agreement with recent studies of the trapped radiation belt electrons. Based on an extensive database of 16 years of corrected MEE flux observations (40-400 keV) from the Research with Adaptive Particle Imaging Detector (RAPID)/Imaging Electron Spectrometer (IES) instrument on board the Cluster mission, Smirnov et al. (2019) reveal that the variability of the outer belt electrons (L-shell 4-6) exhibits a pattern very close to the AE index. Furthermore, Katsavrias et al. (2021) confirm, based on 9 years of electron measurements from GOES-13, 14 and 15, that the trapped electron fluxes at energies in the interval 10-100 keV are well correlated with the AE index consistent with substorm injected source particles.
In the case of the high energy tail of the MEE precipitation, Figure 5 suggests that only 22% of the >292 keV flux variability are described by the daily AE variability. Ødegaard et al. (2017) shows that these higher energy electrons typically peak 1-2 days after the onset of a geomagnetic storm. The progressive time delays of relativistic electrons has also been identified by e.g., Boynton et al. (2016) and Mourenas et al. (2019). The delay implies that it takes time to accelerate the injected seed electrons into several 100 s keV. It will also require a continuous acceleration mechanism. This could be provided during repeated substorm activity, as the source particles will fuel VLF waves and the injected seed particles can be energized as they drift across the substorm induced electric field on the nightside. This means that the precipitating >292 keV fluxes are a product of both the ongoing activity, as well as the substorm activity during the previous days. This potential link is made viable by the simple accumulation of the AE values in Figure 5B. However, the identification of substorm/nonsubstorm modes might be problematic using the daily resolved AE index, in particular for long active geomagnetic periods. Moreover, VLF acceleration and scattering are not limited to substorm periods. Katsavrias et al. (2021) also find a reduced correlation between 100-350 keV trapped electron flux and the AE index compared to the 10-100 keV electrons. These higher energies display, however, a strong dependence on the solar wind speed. They suggest that it implies that the acceleration and loss of the seed energies are not purely substorm driven but rather depend on convection driven by the fast solar wind and/or ULF driven inward diffusion. Smirnov et al. (2019) demonstrate a high positive correlation between the 40-400 keV radiation belt electrons and the AE index and solar wind dynamic pressure. Furthermore, Stepanov et al. (2021) confirm the role of solar wind speed as one of the most important predictors for transporting electron fluxes from the plasmasheet to the radiation belt region. Similarly, Ødegaard et al. (2017) show how the >292 keV flux variability depends on the Akasufos's coupling function. Boynton et al. (2016) developed forecast models for MEE and highly relativistic electrons using the coupling function proposed by (Boynton et al., 2011) and the Dst index. The latter confirms the solar wind dependence while simultaneously accounting for the direction of the Northward interplanetary field direction. Neither the solar wind dependence, nor the coupling function dependence does, however, exclude substorms and the associated VLF wave generation as the working mechanism. In  Frontiers in Astronomy and Space Sciences | www.frontiersin.org October 2021 | Volume 8 | Article 714146 9 fact, Miyoshi et al. (2013) find that High Speed Solar Wind Streams (HSSWS) alone is not sufficient to cause relativistic electron flux enhancement in the outer radiation belt, but strongly depend on IMF-Bz. They state that this would not be the case if radial transport via ULF waves is the primary mechanism. Miyoshi and Kataoka (2008) suggested that the IMF-Bz dependence could be attributed to substorm occurrence during HSSWS events where the electrons are accelerated via VLF waves. This is in line with the relativistic electron flux enhancements found during intervals of prolonged substorm activity operating on a timescale of the order of days (Bühler & Desorgher, 2002;Meredith et al., 2003).
The lifetimes applied in the MEE proxy are impacted by the lifespan of the MEE particles in the outer radiation belt. The lifespan depends on the energy, radial distance from Earth, and the level of geomagnetic activity. The lifetimes for 100 keV electrons at a radial distance of about four are approximately 3.6 days and 13 h for quiet and active geomagnetic conditions, respectively, and are increasing rapidly with energy to about 131 and 17 days for 1 MeV electrons (Orlova et al., 2016). The expected lifetime can, however, be rapidly reduced in the case of magnetopause shadowing where trapped particles over several radial distances are lost to the magnetopause. This is a consequence of a sudden dynamic pressure increase in the solar wind, alongside convection driven outward radial transport (Turner et al., 2012). Applying a fixed lifetime as suggested in the MEE proxy is therefore unrealistic. The optimized lifetime applied in the MEE proxy is thus only the empirical average based on the applied data. It does, however, demonstrate the potential of an AE based MEE proxy , where the precipitation fluxes are an accumulated effect of both the current and previous geomagnetic activity.
This study reveals the high predictive capabilities of the AE index for MEE precipitation, and how to account for the delayed response of the high energy tail (> 300 keV). Hence, it offers a potential improvement to the current MEE parameterization included in CMIP6 recommendation (Matthes et al., 2017). Nesse Tyssøy et al. (2019) compares the loss cone estimate to  the CMIP recommendation and shows an overall underestimation of basic flux strength about one order of magnitude arises from utilizing 0°detector electron fluxes. As this is the same data used to develop the MEE proxy model it is likely that the same assessment will apply if compared to the CMIP recommendation. Furthermore, Nesse Tyssøy et al. (2019) showed that the CMIP recommendation generally captured the initial phase of the storm fluxes, but fell short in respect to reproducing elevated flux levels during the recovery phase of CIRdriven storms. As such, the energy dependent lifetimes applied in the MEE proxy , taking into account the accumulated geomagnetic activity, is likely to avoid this pitfall. However, for AE to be used as a proxy for MEE precipitation in a more advanced MEE precipitation model, the seasonal bias needs further examination as it could be a source of unequal distribution between the two hemispheres. In addition, Figures 6, 8, 10 show a potential solar cycle bias of a general underestimation in the declining phase. To which extent this is due to the dynamical expansion of the auroral oval and the equatorward shift of the electrojets, where the AE-stations cannot reliably monitor their intensity, needs to be explored. The small number of magnetometer stations and their uneven spatial distribution implies that small perturbations (e.g., isolated substorms and pseudo-breakups) can be undetected and large deflections underestimated if they are constrained in longitude or are located at latitudes poleward or equatorward of the AE station network (Gjerloev et al., 2004). Similarly, the MEE precipitation region does not cover a fixed latitude band, and its dependence with geomagnetic activity is explored in a parallel study. Alternatively, the solar cycle bias might reflect the type of solar wind driver responsible for the geomagnetic disturbances as Corotating Interaction Region (CIR)/ High Speed Solar Wind Streams (HSSWS) and Coronal Mass Ejection (CME) driven geomagnetic storms dominate different phases of the solar cycle. For example, the lifetime applied in the MEE proxy can vary in weak but long lasting CIR/HSSWS compared to a short but powerful CME event.
In summary, this study demonstrates that simple AE based MEEproxies have the capability of explaining at 72-86% of the detected MEE precipitation variance on a daily scale. The model shows, however, caveats in respect to the solar cycle and extreme events that summon further investigations. Nevertheless, 90% of the modelled flux values deviate less than a factor of 2.5 from the observed NOAA/ POES MEPED fluxes throughout a full solar cycle. By evaluating the different energy channels, > 43, >114, and >292 keV, independently, the model enables a realistic description of the time dependent energy spectrum. This finding will form the base of a new MEE model to be used for future studies of the energetic electron precipitation impact on the atmosphere.

AUTHOR CONTRIBUTIONS
HNT conceived idea of this paper, in charge of the analysis, and writing as lead author. NP contributed to the idea of this paper and editing of the manuscript. EMB is responsible for constructing the MLT dependent daily resolved electron flux and editing of the manuscript. CS-J contributed to editing of the manuscript. JAS contributed to editing of the manuscript.

FUNDING
The study is supported by the Norwegian Research Council (NRC) under contract 223252, 302040 and 287427. HNT further acknowledge the Young CAS (Centre for Advanced Studies) fellow program.