Time History of Events and Macroscale Interactions during Substorms (THEMIS) spacecraft have been sampling the subsolar magnetosheath since the first dayside science phase in 2008, and we finally have observations over a solar cycle. However, we show that the solar wind coverage during these magnetosheath intervals is not always consistent with the solar wind conditions throughout the same year. This has implications for studying phenomena whose occurrence depends strongly on solar wind parameters. We demonstrate this with magnetosheath jets-flows of enhanced earthward dynamic pressure in the magnetosheath. Jets emerge from the bow shock, and some of them can go on and collide into the magnetopause. Their occurrence is highly linked to solar wind conditions, particularly the orientation of the interplanetary magnetic field, as jets are mostly observed downstream of the quasi-parallel shock. We study the yearly occurrence rates of jets recorded by THEMIS over solar cycle 24 (2008-2019) and find that they are biased due to differences in spacecraft orbits and uneven sampling of solar wind conditions during the different years. Thus, we instead use the THEMIS observations and their corresponding solar wind conditions to develop a model of how jet occurrence varies as a function of solar wind conditions. We then use OMNI data of the whole solar cycle to estimate the unbiased yearly jet occurrence rates. For comparison, we also estimate jet occurrence rates during solar cycle 23 (1996-2008). Our results suggest that there is no strong solar cycle dependency in jet formation.

Magnetosheath jets are localized plasma structures with high dynamic pressure which are frequently observed downstream of the Earth's bow shock. In this work we analyze Magnetospheric MultiScale magnetic field and plasma data and show that jets can be found in the quasi-perpendicular magnetosheath in regions permeated by Mirror mode waves (MMWs). We show that structures identified as jets by their enhanced dynamic pressure can have very different internal structure, with variable signatures in magnetic field magnitude and components, velocity, and density and can be associated to ion distribution functions of various types. This suggests that jets observed in the quasi-perpendicular magnetosheath are generated by different mechanisms. We find that jets can be related to traveling foreshocks, flux transfer events, and some have MMWs inside them. Our results suggest that some jets have a local source and their formation does not depend on upstream structures. We find that different types of ion distributions can exist inside the jets, while in some cases anisotropic distributions are present, in others counterstreaming distributions exist. We also show that for jets with MMWs inside them, ion distributions can be modulated. This highlights the importance of using ion distributions to identify and classify different types of jets.

The plasma environment at a comet can be divided into different regions with distinct plasma characteristics. Two such regions are the solar wind ion cavity, which refers to the part of the outer coma that does not contain any solar wind ions anymore; and the diamagnetic cavity, which is the region of unmagnetized plasma in the innermost coma. From theory and previous observations, it was thought that under usual circumstances no solar wind ion should be observable near or inside of the diamagnetic cavity. For the first time, we report on five observations that show that protons near solar wind energies can also be found inside the diamagnetic cavity. We characterize these proton signatures, where and when they occur, and discuss possible mechanisms that could lead to protons penetrating the inner coma and traversing the diamagnetic cavity boundary. By understanding these observations, we hope to better understand the interaction region of the comet with the solar wind under nonstandard conditions. The protons detected inside the diamagnetic cavity have directions and energies consistent with protons of solar wind origin. The five events occur only at intermediate gas production rates and low cometocentric distances. Charge transfer reactions, high solar wind dynamic pressure and a neutral gas outburst can be ruled out as causes. We suggest that the anomalous appearance of protons in the diamagnetic cavity is due to a specific solar wind configuration where the solar wind velocity is parallel to the interplanetary magnetic field, thus inhibiting mass-loading and deflection.

Magnetosheath jets are dynamic pressure enhancements observed in the terrestrial magnetosheath. Their generation mechanisms are currently debated but the majority of jets can be linked to foreshock processes. Recent results showed that jets are less numerous when coronal mass ejections (CMEs) cross the magnetosheath and more numerous when stream interaction regions (SIRs) cross it. Here, we show for the first time how the pronounced substructures of CMEs and SIRs are related to jet production. We distinguish between compression and magnetic ejecta (ME) regions for the CME as well as compression region associated with the stream interface and high-speed streams (HSSs) for the SIR. Based on THEMIS and OMNI data covering 2008-2021, we show the 2D probability distribution of jet occurrence using the cone angle and Alfvén Mach number. We compare this distribution with the values within each solar wind (SW) structure. We find that both high cone angles and low Alfvén Mach numbers within CME-MEs are unfavorable for jet production as they may inhibit a well-defined foreshock region. 1D histograms of all parameters show, which SW parameters govern jet occurrence in each SW structure. In terms of the considered parameters the most favorable conditions for jet generation are found for HSSs due to their associated low cone angles, low densities, and low magnetic field strengths.

We present in-depth analysis of three southward-moving meso-scale (ion-to magnetohydrodynamic-scale) flux transfer events (FTEs) and subsequent crossing of a reconnecting magnetopause current sheet (MPCS), which were observed on 8 December 2015 by the Magnetospheric Multiscale spacecraft in the subsolar region under southward and duskward magnetosheath magnetic field conditions. We aim to understand the generation mechanism of ion-scale magnetic flux ropes (ISFRs) and to reveal causal relationship among magnetic field structures, electromagnetic energy conversion, and kinetic processes in magnetic reconnection layers. Results from magnetic field reconstruction methods are consistent with a flux rope with a length of about one ion inertial length growing from an electron-scale current sheet (ECS) in the MPCS, supporting the idea that ISFRs can be generated through secondary reconnection in an ECS. Grad-Shafranov reconstruction applied to the three FTEs shows that the FTEs had axial orientations similar to that of the ISFR. This suggests that these FTEs also formed through the same secondary reconnection process, rather than multiple X-line reconnection at spatially separated locations. Four-spacecraft observations of electron pitch-angle distributions and energy conversion rate suggest that the ISFR had three-dimensional magnetic topology and secondary reconnection was patchy or bursty. Previously reported positive and negative values of , with magnitudes much larger than expected for typical MP reconnection, were seen in both magnetosheath and magnetospheric separatrix regions of the ISFR. Many of them coexisted with bi-directional electron beams and intense electric field fluctuations around the electron gyrofrequency, consistent with their origin in separatrix activities.

The Martian crustal magnetic anomalies present a varied, asymmetric obstacle to the imposing draped interplanetary magnetic field (IMF) and solar wind plasma. Magnetic reconnection, a ubiquitous plasma phenomenon responsible for transferring energy and changing magnetic field topology, has been observed throughout the Martian magnetosphere. More specifically, reconnection can occur as a result of the interaction between crustal fields and the IMF, however, the global implications and changes to the overall magnetospheric structure of Mars have yet to be fully understood. Here, we present an analysis to determine these global implications by investigating external conditions that favor reconnection with the underlying crustal anomalies at Mars. To do so, we plot a map of the crustal anomalies' strength and orientation compiled from magnetic field data collected throughout the Mars Atmosphere and Volatile EvolutioN (MAVEN) mission. Then, we create "shear maps" which calculate and plot the angle of shear between the crustal fields and a chosen external field orientation. From there we define a "shear index" to quantify the susceptibility of a region to undergo reconnection based on a given overlaid, external field orientation and the resulting shear map for that region. We demonstrate that the shear analysis technique augments analysis of local reconnection events and suggests southward IMF conditions should favor dayside magnetic reconnection on a more global scale at Mars.

Current inference techniques for processing multi-needle Langmuir probe (m-NLP) data are often based on adaptations of the Orbital Motion-Limited (OML) theory which relies on several simplifying assumptions. Some of these assumptions, however, are typically not well satisfied in actual experimental conditions, thus leading to uncontrolled uncertainties in inferred plasma parameters. In order to remedy this difficulty, three-dimensional kinetic particle in cell simulations are used to construct a synthetic data set, which is used to compare and assess different m-NLP inference techniques. Using a synthetic data set, regression-based models capable of inferring electron density and satellite potentials from 4-tuples of currents collected with fixed-bias needle probes similar to those on the NorSat-1 satellite, are trained and validated. The regression techniques presented show promising results for plasma density inferences with RMS relative errors less than 20%, and satellite potential inferences with RMS errors less than 0.2 V for potentials ranging from -6 to -1 V. The new inference approaches presented are applied to NorSat-1 data, and compared with existing state-of-the-art inference techniques.

In Fourier time-frequency power spectrograms of satellite magnetic field data, electromagnetic ion cyclotron (EMIC) waves may feature discrete, rising tone structures that rapidly increase in frequency. Using data from the Van Allen Probes Electric and Magnetic Field Instrument Suite and Integrated Science (EMFISIS) fluxgate magnetometer, we conducted a statistical study of EMIC waves from September 2012 through June 2016. We compared the occurrence rates and spatial distributions for all EMIC waves with those for rising tone EMIC waves as a function of magnetic local time (MLT) and shell, as well as a function of and in solar-magnetic (SM) coordinates. Overall, EMIC waves occurred during 2.4% of the time period considered, but rising tone EMIC waves were only found during 0.2% of the time period considered. About 7%-8% of the minutes of orbital coverage with H+ or He+ band EMIC waves had rising tones. The regions of peak occurrence rates for H+ and He+ band waves, as well as waves with rising tones, were found in the noon and dusk sectors for 4 <  < 6. The preferred regions for H+ waves as a function of and SM suggest an association with magnetospheric compressions near noon and interactions between plumes and the ring current near dusk. Peak occurrence rates for O+ band waves were found between 2 <  < 4 at all MLT, and over a wide range of shells near dusk. No rising tones were found in the O+ band.

The Thermosphere-Ionosphere-Mesosphere Energetics and Dynamics (TIMED) satellite has been making observations of the mesosphere and lower thermosphere (MLT) region for two decades. The TIMED Doppler Interferometer (TIDI) measures the neutral winds using four orthogonal telescopes. In this study, the line of sight (LOS) winds from individual telescopes are compared to the measurements from the Ionospheric Connection Explorer's (ICON's) Michelson Interferometer for Global High-resolution Thermospheric Imaging (MIGHTI) instrument from 90 to 100 km altitude during 2020. With the MIGHTI vector winds projected onto the LOS direction of each TIDI telescope, coincidences of the two data sets are found. The four telescopes perform differently and the performance depends on the satellite configuration and local solar zenith angle. Measurements from the coldside telescopes, Telescope 1 (Tel1) and Telescope 2 (Tel2), are better correlated with the MIGHTI winds in general with Tel2 having higher correlation coefficients across all conditions. The performance of Tel1 is comparable to that of Tel2 during backward flight while showing systematic errors larger than the average wind speeds during forward flight. The warmside LOS winds from Telescope 3 (Tel3) and Telescope 4 (Tel4) vary widely in magnitude, especially on the nightside. Compared with MIGHTI winds, the Tel4 measurements have the weakest correlation, while the Tel3 performance is comparable to that of the coldside telescopes during the ascending phase but deteriorates during the descending phase. Based on the TIDI/MIGHTI comparisons, figures of merit are generated to quantify the quality of measurements from individual telescopes in different configurations.

Thin current sheets (TCSs) have been postulated to be a necessary precondition for reconnection onset. Magnetic reconnection X-lines in the magnetotail have been observed to be more common duskward of midnight. We take advantage of the MMS tetrahedral formation during the 2017-2020 MMS tail seasons to calculate the thickness of the cross-tail neutral sheet relative to ion gyroradius. While a similar technique was applied to Cluster data, current sheet thickness over a broader range of radial distances has not been robustly explored before this study. We compare our analysis to recent theories regarding mechanisms of tail current sheet thinning and to recent simulations. We find MMS spent more than twice as long in ion-scale TCSs in the pre-midnight sector than post-midnight, despite nearly even plasma sheet dwell time. The dawn-dusk asymmetry in the distribution of Ion Diffusion Regions, as previously reported in relation to regions of TCSs, is also analyzed.

We analyze observations of a solar energetic particle (SEP) event at Rosetta's target comet 67P/Churyumov-Gerasimenko during 6-10 March 2015. The comet was 2.15 AU from the Sun, with the Rosetta spacecraft approximately 70 km from the nucleus placing it deep inside the comet's coma and allowing us to study its response. The Eastern flank of an interplanetary coronal mass ejection (ICME) also encountered Rosetta on 6 and 7 March. Rosetta Plasma Consortium data indicate increases in ionization rates, and cometary water group pickup ions exceeding 1 keV. Increased charge exchange reactions between solar wind ions and cometary neutrals also indicate increased upstream neutral populations consistent with enhanced SEP induced surface activity. In addition, the most intense parts of the event coincide with observations interpreted as an infant cometary bow shock, indicating that the SEPs may have enhanced the formation and/or intensified the observations. These solar transient events may also have pushed the cometopause closer to the nucleus. We track and discuss characteristics of the SEP event using remote observations by SOHO, WIND, and GOES at the Sun, in situ measurements at Solar Terrestrial Relations Observatory Ahead, Mars and Rosetta, and ENLIL modeling. Based on its relatively prolonged duration, gradual and anisotropic nature, and broad angular spread in the heliosphere, we determine the main particle acceleration source to be a distant ICME which emerged from the Sun on 6 March 2015 and was detected locally in the Martian ionosphere but was never encountered by 67P directly. The ICME's shock produced SEPs for several days which traveled to the in situ observation sites via magnetic field line connections.

We present a statistical study of Jupiter's disk X-ray emissions using 19 years of Chandra X-Ray Observatory (CXO) observations. Previous work has suggested that these emissions are consistent with solar X-rays elastically scattered from Jupiter's upper atmosphere. We showcase a new pulse invariant (PI) filtering method that minimizes instrumental effects which may produce unphysical trends in photon counts across the nearly two-decade span of the observations. We compare the CXO results with solar X-ray flux data from the Geostationary Operational Environmental Satellites X-ray Sensor for the wavelength band 1-8 Å (long channel), to quantify the correlation between solar activity and Jovian disk counts. We find a statistically significant Pearson's Correlation Coefficient of 0.9, which confirms that emitted Jovian disk X-rays are predominantly governed by solar activity. We also utilize the high spatial resolution of the High Resolution Camera Instrument on-board the CXO to map the disk photons to their positions on Jupiter's surface. Voronoi tessellation diagrams were constructed with the Juno Reference Model through Perijove 9 internal field model overlaid to identify any spatial preference of equatorial photons. After accounting for area and scattering across the curved surface of the planet, we find a preference of Jovian disk emission at 2-3.5 Gauss surface magnetic field strength. This suggests that a portion of the disk X-rays may be linked to processes other than solar scattering: the spatial preference associated with magnetic field strength may imply increased precipitation from the radiation belts, as previously postulated.

The properties and acceleration mechanisms of electrons (<200 keV) associated with a pair of tailward traveling flux ropes and accompanied reconnection X-lines in Earth's plasma sheet are investigated with MMS measurements. Energetic electrons are enhanced on both boundaries and core of the flux ropes. The power-law spectra of energetic electrons near the X-lines and in flux ropes are harder than those on flux rope boundaries. Theoretical calculations show that the highest energy of adiabatic electrons is a few keV around the X-lines, tens of keV immediately downstream of the X-lines, hundreds of keV on the flux rope boundaries, and a few MeV in the flux rope cores. The X-lines cause strong energy dissipation, which may generate the energetic electron beams around them. The enhanced electron parallel temperature can be caused by the curvature-driven Fermi acceleration and the parallel electric potential. Betatron acceleration due to the magnetic field compression is strong on flux rope boundaries, which enhances energetic electrons in the perpendicular direction. Electrons can be trapped between the flux rope pair due to mirror force and parallel electric potential. Electrostatic structures in the flux rope cores correspond to potential drops up to half of the electron temperature. The energetic electrons and the electron distribution functions in the flux rope cores are suggested to be transported from other dawn-dusk directions, which is a 3-dimensional effect. The acceleration and deceleration of the Betatron and Fermi processes appear alternately indicating that the magnetic field and plasma are turbulent around the flux ropes.

Observations of magnetic clouds, within interplanetary coronal mass ejections (ICMEs), are often well described by flux rope models. Most of these assume either a cylindrical or toroidal geometry. In some cases, these models are also capable of accounting for non-axisymmetric cross-sections but they generally all assume axial invariance. It can be expected that any ICME, and its flux rope, will be deformed along its axis due to influences such as the solar wind. In this work, we aim to develop a writhed analytical magnetic flux rope model which would allow us to analytically describe a flux rope structure with varying curvature and torsion so that we are no longer constrained to a cylindrical or toroidal geometry. In this first iteration of our model we will solely focus on a circular cross-section of constant size. We describe our flux rope geometry in terms of a parametrized flux rope axis and a parallel transport frame. We derive expressions for the axial and poloidal magnetic field components under the assumption that the total axial magnetic flux is conserved. We find an entire class of possible solutions, which differ by the choice of integration constants, and present the results for a specific example. In general, we find that the twist of the magnetic field locally changes when the geometry deviates from a cylinder or torus. This new approach also allows us to generate completely new types of in situ magnetic field profiles which strongly deviate from those generated by cylindrical or toroidal models.

The Super Dual Auroral Radar Network (SuperDARN) is a collection of radars built to study ionospheric convection. We use a 7-year archive of SuperDARN convection maps, processed in 3 different ways, to build a statistical understanding of dusk-dawn asymmetries in the convection patterns. We find that the data set processing alone can introduce a bias which manifests itself in dusk-dawn asymmetries. We find that the solar wind clock angle affects the balance in the strength of the convection cells. We further find that the location of the positive potential foci is most likely observed at latitudes of 78° for long periods (>300 min) of southward interplanetary magnetic field (IMF), as opposed to 74° for short periods (<20 min) of steady IMF. For long steady dawnward IMF the median is also at 78°. For long steady periods of duskward IMF, the positive potential foci tends to be at lower latitudes than the negative potential and vice versa during dawnward IMF. For long periods of steady Northward IMF, the positive and negative cells can swap sides in the convection pattern. We find that they move from ∼0-9 MLT to 15 MLT or ∼15-23 MLT to 10 MLT, which reduces asymmetry in the average convection cell locations for Northward IMF. We also investigate the width of the region in which the convection returns to the dayside, the return flow width. Asymmetries in this are not obvious, until we select by solar wind conditions, when the return flow region is widest for the negative convection cell during Southward IMF.

Here, we extend the scope of the Gamayunov and Engebretson (2021, hereinafter Paper 1), https://doi.org/10.1029/2021JA029247 work by analyzing the low frequency ultra-low-frequency (ULF) wave power spectra in the Earth's inner magnetosphere during high speed stream (HSS) and quiet solar wind (QSW) driving conditions in the upstream solar wind (SW) and comparing our results to the results of Paper 1, where the statistics of ULF wave power spectra during coronal mass ejections (CMEs) are presented. The most important results of our statistical and comparative analyses are as follows. (a) During CMEs, HSSs, and QSW, the magnetic field power spectra of the transverse and compressional fluctuations are well approximated by power laws in the ∼mHz-Hz frequency range, where on average the parameters of power law fits during CMEs and HSSs are close, and those during QSW differ considerably from the respective parameters during CMEs and HSSs. (b) The dominance of the average compressional power over the average transverse power for the low frequency ULF waves during the 0 < SYM/H ≲ 25 nT geomagnetic conditions may serve as a proxy of HSSs in the upstream SW, whereas the opposite relation between the average powers is an indication of CMEs. (c) Independently of the SW driving conditions, a turbulent energy cascade from low frequencies in the ULF wave frequency range into the higher frequency range exists in the Earth's inner magnetosphere, supplying the nonthermal electromagnetic seed fluctuations needed for the growth of electromagnetic ion cyclotron waves (∼Hz) due to relaxation of unstable distributions of energetic magnetospheric ions.

Electron diffusion by whistler-mode chorus waves is one of the key processes controlling the dynamics of relativistic electron fluxes in the Earth's radiation belts. It is responsible for the acceleration of sub-relativistic electrons injected from the plasma sheet to relativistic energies as well as for their precipitation and loss into the atmosphere. Based on analytical estimates of chorus wave-driven quasi-linear electron energy and pitch-angle diffusion rates, we provide analytical steady-state solutions to the corresponding Fokker-Planck equation for the relativistic electron distribution and flux. The impact on these steady-state solutions of additional electromagnetic ion cyclotron waves, and of ultralow frequency waves are examined. Such steady-state solutions correspond to hard energy spectra at 1-4 MeV, dangerous for satellite electronics, and represent attractors for the system dynamics in the presence of sufficiently strong driving by continuous injections of 10-300 keV electrons. Therefore, these analytical steady-state solutions provide a simple means for estimating the most extreme electron energy spectra potentially encountered in the outer radiation belt, despite the great variability of injections and plasma conditions. These analytical steady-state solutions are compared with numerical simulations based on the full Fokker-Planck equation and with relativistic electron flux spectra measured by satellites during one extreme event and three strong events of high time-integrated geomagnetic activity, demonstrating a good agreement.

Under quasi-radial interplanetary magnetic fields (IMF), foreshock turbulence can have an impact on the magnetosheath and cusps depending on the location of the quasi-parallel shock. We perform three-dimensional simulations of Earth's dayside magnetosphere using the hybrid code HYPERS, and compare northward and southward quasi-radial IMF configurations. We study the magnetic field configuration, fluctuations in the magnetosheath and the plasma in the regions around the northern cusp. Under northward IMF with Earthward , there is a time-varying plasma depletion layer immediately outside the northern cusp. In the southward IMF case, the impact of foreshock turbulence and high-speed jets, together with magnetopause reconnection, can lead to strong density enhancements in the cusp.

Fast divergent flows of electrons and ions in the magnetotail plasma sheet are conventionally interpreted as a key reconnection signature caused by the magnetic topology change at the X-line. Therefore, reversals of the -component ( ) of the plasma flow perpendicular to the magnetic field must correlate with the sign changes in the north-south component of the magnetic field ( ). Here we present observations of the flow reversals that take place with no correlated reversals. We report six such events, which were measured with the high-resolution plasma and fields instruments of the Magnetospheric Multiscale mission. We found that electron flow reversals in the absence of reversals (a) have amplitudes of ∼1,000-2,000 km s and durations of a few seconds; (b) are embedded into larger-scale ion flow reversals with enhanced ion agyrotropy; and (c) compared with conventional reconnection outflows around the electron diffusion regions (EDRs), have less (if ever) pronounced electron agyrotropy, dawnward electron flow amplitude, and electric field strength toward the neutral sheet, although their energy conversion parameters, including the Joule heating rate, are quite substantial. These results suggest that such flow reversals develop in the ion-demagnetization regions away from electron-scale current sheets, in particular the EDRs, and yet they play an important role in the energy conversion. These divergent flows are interpreted as precursors of the flow-driven reconnection onsets provided by the ion tearing or the ballooning/interchange instability.

The dynamics of the Jovian magnetosphere is controlled by the interplay of the planet's fast rotation, its solar-wind interaction and its main plasma source at the Io torus, mediated by coupling processes involving its magnetosphere, ionosphere, and thermosphere. At the ionospheric level, these processes can be characterized by a set of parameters including conductances, field-aligned currents, horizontal currents, electric fields, transport of charged particles along field lines including the fluxes of electrons precipitating into the upper atmosphere which trigger auroral emissions, and the particle and Joule heating power dissipation rates into the upper atmosphere. Determination of these key parameters makes it possible to estimate the net transfer of momentum and energy between Jovian upper atmosphere and equatorial magnetosphere. A method based on a combined use of Juno multi-instrument data and three modeling tools was developed by Wang et al. (2021, https://doi.org/10.1029/2021ja029469) and applied to an analysis of the first nine orbits to retrieve these parameters along Juno's magnetic footprint. We extend this method to the first 30 Juno science orbits and to both hemispheres. Our results reveal a large variability of these parameters from orbit to orbit and between the two hemispheres. They also show dominant trends. Southern current systems are consistent with the generation of a region of sub-corotating ionospheric plasma flows, while both super-corotating and sub-corotating plasma flows are found in the north. These results are discussed in light of the previous space and ground-based observations and currently available models of plasma convection and current systems, and their implications are assessed.

We present a comprehensive analysis of the processes that lead to quasilinear pitch-angle-scattering loss of electrons from the  < 4 region of the Earth's inner magnetosphere during geomagnetically quiet times. We consider scattering via Coulomb collisions, hiss waves, lightning-generated whistler (LGW) waves, waves from ground-based very-low-frequency (VLF) transmitters, and electromagnetic ion cyclotron (EMIC) waves. The amplitude, frequency, and wave normal angle spectra of these waves are parameterized with empirical wave models, which are then used to compute pitch-angle diffusion coefficients. From these coefficients, we estimate the decay timescales, or lifetimes, of 30 keV to 4 MeV electrons and compare the results with timescales obtained from in-situ observations. We demonstrate good quantitative agreement between the two over most of the and energy range under investigation. Our analysis suggests that the electron decay timescales are very sensitive to the choice of plasmaspheric density model. At  < 2, where our theoretical lifetimes do not agree well with the observations, we show that including Coulomb energy drag (ionization energy loss) in our calculations significantly improves the quantitative agreement with the observed decay timescales. We also use an accurate model of the geomagnetic field to provide an estimate of the effect that the drift-loss cone has on the theoretically calculated electron lifetimes, which are usually obtained using an axisymmetric dipole field.