@article{HaasShpritsAllisonetal.2022, author = {Haas, Bernhard and Shprits, Yuri Y. and Allison, Hayley and Wutzig, Michael and Wang, Dedong}, title = {Which parameter controls ring current electron dynamics}, series = {Frontiers in astronomy and space sciences}, volume = {9}, journal = {Frontiers in astronomy and space sciences}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-987X}, doi = {10.3389/fspas.2022.911002}, pages = {11}, year = {2022}, abstract = {Predicting the electron population of Earth's ring current during geomagnetic storms still remains a challenging task. In this work, we investigate the sensitivity of 10 keV ring current electrons to different driving processes, parameterised by the Kp index, during several moderate and intense storms. Results are validated against measurements from the Van Allen Probes satellites. Perturbing the Kp index allows us to identify the most dominant processes for moderate and intense storms respectively. We find that during moderate storms (Kp < 6) the drift velocities mostly control the behaviour of low energy electrons, while loss from wave-particle interactions is the most critical parameter for quantifying the evolution of intense storms (Kp > 6). Perturbations of the Kp index used to drive the boundary conditions at GEO and set the plasmapause location only show a minimal effect on simulation results over a limited L range. It is further shown that the flux at L \& SIM; 3 is more sensitive to changes in the Kp index compared to higher L shells, making it a good proxy for validating the source-loss balance of a ring current model.}, language = {en} } @article{ShpritsDrozdovSpasojevicetal.2016, author = {Shprits, Yuri Y. and Drozdov, Alexander and Spasojevic, Maria and Kellerman, Adam C. and Usanova, Maria E. and Engebretson, Mark J. and Agapitov, Oleksiy V. and Zhelavskaya, Irina and Raita, Tero J. and Spence, Harlan E. and Baker, Daniel N. and Zhu, Hui and Aseev, Nikita}, title = {Wave-induced loss of ultra-relativistic electrons in the Van Allen radiation belts}, series = {Nature Communications}, volume = {7}, journal = {Nature Communications}, publisher = {Nature Publ. Group}, address = {London}, issn = {2041-1723}, doi = {10.1038/ncomms12883}, pages = {7}, year = {2016}, language = {en} } @article{AseevShpritsWangetal.2019, author = {Aseev, Nikita and Shprits, Yuri Y. and Wang, Dedong and Wygant, John and Drozdov, Alexander and Kellerman, Adam C. and Reeves, Geoffrey D.}, title = {Transport and loss of ring current electrons inside geosynchronous orbit during the 17 March 2013 storm}, series = {Journal of geophysical research : Space physics}, volume = {124}, journal = {Journal of geophysical research : Space physics}, number = {2}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2018JA026031}, pages = {915 -- 933}, year = {2019}, abstract = {Ring current electrons (1-100 keV) have received significant attention in recent decades, but many questions regarding their major transport and loss mechanisms remain open. In this study, we use the four-dimensional Versatile Electron Radiation Belt code to model the enhancement of phase space density that occurred during the 17 March 2013 storm. Our model includes global convection, radial diffusion, and scattering into the Earth's atmosphere driven by whistler-mode hiss and chorus waves. We study the sensitivity of the model to the boundary conditions, global electric field, the electric field associated with subauroral polarization streams, electron loss rates, and radial diffusion coefficients. The results of the code are almost insensitive to the model parameters above 4.5 RERE, which indicates that the general dynamics of the electrons between 4.5 RE and the geostationary orbit can be explained by global convection. We found that the major discrepancies between the model and data can stem from the inaccurate electric field model and uncertainties in lifetimes. We show that additional mechanisms that are responsible for radial transport are required to explain the dynamics of ≥40-keV electrons, and the inclusion of the radial diffusion rates that are typically assumed in radiation belt studies leads to a better agreement with the data. The overall effect of subauroral polarization streams on the electron phase space density profiles seems to be smaller than the uncertainties in other input parameters. This study is an initial step toward understanding the dynamics of these particles inside the geostationary orbit.}, language = {en} } @misc{ShpritsHorneKellermanetal.2018, author = {Shprits, Yuri Y. and Horne, Richard B. and Kellerman, Adam C. and Drozdov, Alexander}, title = {The dynamics of Van Allen belts revisited}, series = {Nature physics}, volume = {14}, journal = {Nature physics}, number = {2}, publisher = {Nature Publ. Group}, address = {London}, issn = {1745-2473}, doi = {10.1038/nphys4350}, pages = {102 -- 103}, year = {2018}, abstract = {In an effort to explain the formation of a narrow third radiation belt at ultra-relativistic energies detected during a solar storm in September 20121, Mann et al.2 present simulations from which they conclude it arises from a process of outward radial diffusion alone, without the need for additional loss processes from higher frequency waves. The comparison of observations with the model in Figs 2 and 3 of their Article clearly shows that even with strong radial diffusion rates, the model predicts a third belt near L* = 3 that is twice as wide as observed and approximately an order of magnitude more intense. We therefore disagree with their interpretation that "the agreement between the absolute fluxes from the model and those observed by REPT [the Relativistic Electron Proton Telescope] shown on Figs 2 and 3 is excellent." Previous studies3 have shown that outward radial diffusion plays a very important role in the dynamics of the outer belt and is capable of explaining rapid reductions in the electron flux. It has also been shown that it can produce remnant belts (Fig. 2 of a long-term simulation study4). However, radial diffusion alone cannot explain the formation of the narrow third belt at multi-MeV during September 2012. An additional loss mechanism is required. Higher radial diffusion rates cannot improve the comparison of model presented by Mann et al. with observations. A further increase in the radial diffusion rates (reported in Fig. 4 of the Supplementary Information of ref. 2) results in the overestimation of the outer belt fluxes by up to three orders of magnitude at energy of 3.4 MeV. Observations at 2 MeV, where belts show only a two-zone structure, were not presented by Mann et al. Moreover, simulations of electrons with energies below 2 MeV with the same diffusion rates and boundary conditions used by the authors would probably produce very strong depletions down to L = 3-3.5, where L is radial distance from the centre of the Earth to the given field line in the equatorial plane. Observations do not show a non-adiabatic loss below L ∼ 4.5 for 2 MeV. Such different dynamics between 2 MeV and above 4 MeV at around L = 3.5 are another indication that particles are scattered by electromagnetic ion cyclotron (EMIC) waves that affect only energies above a certain threshold. Observations of the phase space density (PSD) provide additional evidence for the local loss of electrons. Around L* = 3.5-4 PSD shows significant decrease by an order of magnitude starting in the afternoon of 3 September (Fig. 1a), while PSD above L* = 4 is increasing. The minimum in PSD between L* = 3.5-4 continues to decrease until 4 September. This evolution demonstrates that the loss is not produced by outward diffusion. Radial diffusion cannot produce deepening minima, as it works to smooth gradients. Just as growing peaks in PSD show the presence of localized acceleration5, deepening minima show the presence of localized loss. Figure 1: Time evolution of radiation profiles in electron PSD at relativistic and ultra-relativistic energies. figure 1 a, Similar to Supplementary Fig. 3 of ref. 2, but using TS07D model10 and for μ = 2,500 MeV G-1, K = 0.05 RE G0.5 (where RE is the radius of the Earth). b, Similar to Supplementary Fig. 3 of ref. 2, but using TS07D model and for μ = 700 MeV G-1, corresponding to MeV energies in the heart of the belt. Minimum in PSD in the heart of the multi-MeV electron radiation belt between 3.5 and 4 RE deepening between the afternoon of 3 September and 5 September clearly show that the narrow remnant belt at multi-MeV below 3.5 RE is produced by the local loss. Full size image The minimum in the outer boundary is reached on the evening of 2 September. After that, the outer boundary moves up, while the minimum decreases by approximately an order of magnitude, clearly showing that this main decrease cannot be explained by outward diffusion, and requires additional loss processes. The analysis of profiles of PSD is a standard tool used, for example, in the study about electron acceleration5 and routinely used by the entire Van Allen Probes team. In the Supplementary Information, we show that this analysis is validated by using different magnetic field models. The Supplementary Information also shows that measurements are above background noise. Deepening minima at multi-MeV during the times when the boundary flux increases are clearly seen in Fig. 1a. They show that there must be localized loss, as radial diffusion cannot produce a minimum that becomes lower with time. At lower energies of 1-2 MeV, which corresponds to lower values of the first adiabatic invariant μ (Fig. 1b), the profiles are monotonic between L* = 3-3.5, consistent with the absence of scattering by EMIC waves that affect only electrons above a certain energy threshold6,7,8,9. In summary, the results of the modelling and observations presented by Mann et al. do not lend support to the claim of explaining the dynamics of the ultra-relativistic third Van Allen radiation belt in terms of an outward radial diffusion process alone. While the outward radial diffusion driven by the loss to the magnetopause2 is certainly operating during this storm, there is compelling observational and modelling2,6 evidence that shows that very efficient localized electron loss operates during this storm at multi-MeV energies, consistent with localized loss produced by EMIC waves.}, language = {en} } @article{ZhelayskayaVasileShpritsetal.2019, author = {Zhelayskaya, Irina S. and Vasile, Ruggero and Shprits, Yuri Y. and Stolle, Claudia and Matzka, J{\"u}rgen}, title = {Systematic Analysis of Machine Learning and Feature Selection Techniques for Prediction of the Kp Index}, series = {Space Weather: The International Journal of Research and Applications}, volume = {17}, journal = {Space Weather: The International Journal of Research and Applications}, number = {10}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1542-7390}, doi = {10.1029/2019SW002271}, pages = {1461 -- 1486}, year = {2019}, abstract = {The Kp index is a measure of the midlatitude global geomagnetic activity and represents short-term magnetic variations driven by solar wind plasma and interplanetary magnetic field. The Kp index is one of the most widely used indicators for space weather alerts and serves as input to various models, such as for the thermosphere and the radiation belts. It is therefore crucial to predict the Kp index accurately. Previous work in this area has mostly employed artificial neural networks to nowcast Kp, based their inferences on the recent history of Kp and on solar wind measurements at L1. In this study, we systematically test how different machine learning techniques perform on the task of nowcasting and forecasting Kp for prediction horizons of up to 12 hr. Additionally, we investigate different methods of machine learning and information theory for selecting the optimal inputs to a predictive model. We illustrate how these methods can be applied to select the most important inputs to a predictive model of Kp and to significantly reduce input dimensionality. We compare our best performing models based on a reduced set of optimal inputs with the existing models of Kp, using different test intervals, and show how this selection can affect model performance.}, language = {en} } @article{KimShprits2018, author = {Kim, Kyung-Chan and Shprits, Yuri Y.}, title = {Survey of the Favorable Conditions for Magnetosonic Wave Excitation}, series = {Journal of geophysical research : Space physics}, volume = {123}, journal = {Journal of geophysical research : Space physics}, number = {1}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1002/2017JA024865}, pages = {400 -- 413}, year = {2018}, abstract = {The ratio of the proton ring velocity (VR) to the local Alfven speed (VA), in addition to proton ring distributions, plays a key factor in the excitation of magnetosonic waves at frequencies between the proton cyclotron frequency fcp and the lower hybrid resonance frequency fLHR in the Earth's magnetosphere. Here we investigate whether there is a statistically significant relationship between occurrences of proton rings and magnetosonic waves both outside and inside the plasmapause using particle and wave data from Van Allen Probe-A during the time period of October 2012 to December 2015. We also perform a statistical survey of the ratio of the ring energy (ER, corresponding to VR) to the Alfven energy (EA, corresponding to VA) to determine the favorable conditions under which magnetosonic waves in each of two frequency bands (fcp < f ≤ 0.5 fLHR and 0.5 fLHR < f < fLHR) can be excited. The results show that the magnetosonic waves in both frequency bands occur around the postnoon (12-18 magnetic local time, MLT) sector outside the plasmapause when ER is comparable to or lower than EA, and those in lower-frequency bands (fcp < f ≤ 0.5 fLHR) occur around the postnoon sector inside the plasmapause when ER/EA > ~9. However, there is one discrepancy between occurrences of proton rings and magnetosonic waves in low-frequency bands around the prenoon sector (6-12 MLT) outside the plasmapause, which suggests either that the waves may have propagated during active time from the postnoon sector after being excited during quiet time, or they may have locally excited in the prenoon sector during active time.}, language = {en} } @article{DelCorpoVellanteZhelavskayaetal.2022, author = {Del Corpo, Alfredo and Vellante, Massimo and Zhelavskaya, Irina and Shprits, Yuri Y. and Heilig, Balazs and Reda, Jan and Pietropaolo, Ermanno and Lichtenberger, Janos}, title = {Study of the average ion mass of the dayside magnetospheric plasma}, series = {Journal of geophysical research : Space physics}, volume = {127}, journal = {Journal of geophysical research : Space physics}, number = {10}, publisher = {American Geophysical Union}, address = {Washington, DC}, issn = {2169-9380}, doi = {10.1029/2022JA030605}, pages = {20}, year = {2022}, abstract = {The investigation of heavy ions dynamics and properties in the Earth's magnetosphere is still an important field of research as they play an important role in several space weather aspects. We present a statistical survey of the average ion mass in the dayside magnetosphere made comparing plasma mass density with electron number density measurements and focusing on both spatial and geomagnetic activity dependence. Field line resonance frequency observations across the European quasi-Meridional Magnetometer Array, are used to infer the equatorial plasma mass density in the range of magnetic L-shells 1.6-6.2. The electron number density is derived from local electric field measurements made on Van Allen Probes using the Neural-network-based Upper-hybrid Resonance Determination algorithm. The analysis is conducted separately for the plasmasphere and the plasmatrough during favorable periods for which both the plasma parameters are observed simultaneously. We found that throughout the plasmasphere the average ion mass is similar or equal to 1 amu for a wide range of geomagnetic activity conditions, suggesting that the plasma mainly consist of hydrogen ions, without regard to the level of geomagnetic activity. Conversely, the plasmatrough is characterized by a variable composition, highlighting a heavy ion mass loading that increases with increasing levels of geomagnetic disturbance. During the most disturbed conditions, the average radial structure shows a broad maximum around 3-4 Earth radii, probably correlated with the accumulation of oxygen ions near the plasmapause. Those ions are mostly observed in the post-dawn and pre-dusk longitudinal sectors.}, language = {en} } @article{ShpritsMeniettiDrozdovetal.2018, author = {Shprits, Yuri Y. and Menietti, J. D. and Drozdov, Alexander and Horne, Richard B. and Woodfield, Emma E. and Groene, J. B. and de Soria-Santacruz, M. and Averkamp, T. F. and Garrett, H. and Paranicas, C. and Gurnett, Don A.}, title = {Strong whistler mode waves observed in the vicinity of Jupiter's moons}, series = {Nature Communications}, volume = {9}, journal = {Nature Communications}, publisher = {Nature Publ. Group}, address = {London}, issn = {2041-1723}, doi = {10.1038/s41467-018-05431-x}, pages = {6}, year = {2018}, abstract = {Understanding of wave environments is critical for the understanding of how particles are accelerated and lost in space. This study shows that in the vicinity of Europa and Ganymede, that respectively have induced and internal magnetic fields, chorus wave power is significantly increased. The observed enhancements are persistent and exceed median values of wave activity by up to 6 orders of magnitude for Ganymede. Produced waves may have a pronounced effect on the acceleration and loss of particles in the Jovian magnetosphere and other astrophysical objects. The generated waves are capable of significantly modifying the energetic particle environment, accelerating particles to very high energies, or producing depletions in phase space density. Observations of Jupiter's magnetosphere provide a unique opportunity to observe how objects with an internal magnetic field can interact with particles trapped in magnetic fields of larger scale objects.}, language = {en} } @misc{ShpritsMeniettiDrozdovetal.2018, author = {Shprits, Yuri Y. and Menietti, J. D. and Drozdov, Alexander and Horne, Richard B. and Woodfield, Emma E. and Groene, J. B. and de Soria-Santacruz, M. and Averkamp, T. F. and Garrett, H. and Paranicas, C. and Gurnett, Don A.}, title = {Strong whistler mode waves observed in the vicinity of Jupiter's moons}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {695}, issn = {1866-8372}, doi = {10.25932/publishup-42627}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-426278}, pages = {6}, year = {2018}, abstract = {Understanding of wave environments is critical for the understanding of how particles are accelerated and lost in space. This study shows that in the vicinity of Europa and Ganymede, that respectively have induced and internal magnetic fields, chorus wave power is significantly increased. The observed enhancements are persistent and exceed median values of wave activity by up to 6 orders of magnitude for Ganymede. Produced waves may have a pronounced effect on the acceleration and loss of particles in the Jovian magnetosphere and other astrophysical objects. The generated waves are capable of significantly modifying the energetic particle environment, accelerating particles to very high energies, or producing depletions in phase space density. Observations of Jupiter's magnetosphere provide a unique opportunity to observe how objects with an internal magnetic field can interact with particles trapped in magnetic fields of larger scale objects.}, language = {en} } @article{AdolfsHoqueShprits2022, author = {Adolfs, Marjolijn and Hoque, Mohammed Mainul and Shprits, Yuri Y.}, title = {Storm-time relative total electron content modelling using machine learning techniques}, series = {Remote sensing}, volume = {14}, journal = {Remote sensing}, number = {23}, publisher = {MDPI}, address = {Basel}, issn = {2072-4292}, doi = {10.3390/rs14236155}, pages = {17}, year = {2022}, abstract = {Accurately predicting total electron content (TEC) during geomagnetic storms is still a challenging task for ionospheric models. In this work, a neural-network (NN)-based model is proposed which predicts relative TEC with respect to the preceding 27-day median TEC, during storm time for the European region (with longitudes 30 degrees W-50 degrees E and latitudes 32.5 degrees N-70 degrees N). The 27-day median TEC (referred to as median TEC), latitude, longitude, universal time, storm time, solar radio flux index F10.7, global storm index SYM-H and geomagnetic activity index Hp30 are used as inputs and the output of the network is the relative TEC. The relative TEC can be converted to the actual TEC knowing the median TEC. The median TEC is calculated at each grid point over the European region considering data from the last 27 days before the storm using global ionosphere maps (GIMs) from international GNSS service (IGS) sources. A storm event is defined when the storm time disturbance index Dst drops below 50 nanotesla. The model was trained with storm-time relative TEC data from the time period of 1998 until 2019 (2015 is excluded) and contains 365 storms. Unseen storm data from 33 storm events during 2015 and 2020 were used to test the model. The UQRG GIMs were used because of their high temporal resolution (15 min) compared to other products from different analysis centers. The NN-based model predictions show the seasonal behavior of the storms including positive and negative storm phases during winter and summer, respectively, and show a mixture of both phases during equinoxes. The model's performance was also compared with the Neustrelitz TEC model (NTCM) and the NN-based quiet-time TEC model, both developed at the German Aerospace Agency (DLR). The storm model has a root mean squared error (RMSE) of 3.38 TEC units (TECU), which is an improvement by 1.87 TECU compared to the NTCM, where an RMSE of 5.25 TECU was found. This improvement corresponds to a performance increase by 35.6\%. The storm-time model outperforms the quiet-time model by 1.34 TECU, which corresponds to a performance increase by 28.4\% from 4.72 to 3.38 TECU. The quiet-time model was trained with Carrington averaged TEC and, therefore, is ideal to be used as an input instead of the GIM derived 27-day median. We found an improvement by 0.8 TECU which corresponds to a performance increase by 17\% from 4.72 to 3.92 TECU for the storm-time model using the quiet-time-model predicted TEC as an input compared to solely using the quiet-time model.}, language = {en} } @article{SmirnovShpritsAllisonetal.2022, author = {Smirnov, Artem and Shprits, Yuri Y. and Allison, Hayley and Aseev, Nikita and Drozdov, Alexander and Kollmann, Peter and Wang, Dedong and Saikin, Anthony}, title = {Storm-Time evolution of the Equatorial Electron Pitch Angle Distributions in Earth's Outer Radiation Belt}, series = {Frontiers in astronomy and space sciences}, volume = {9}, journal = {Frontiers in astronomy and space sciences}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-987X}, doi = {10.3389/fspas.2022.836811}, pages = {15}, year = {2022}, abstract = {In this study we analyze the storm-time evolution of equatorial electron pitch angle distributions (PADs) in the outer radiation belt region using observations from the Magnetic Electron Ion Spectrometer (MagEIS) instrument aboard the Van Allen Probes in 2012-2019. The PADs are approximated using a sum of the first, third and fifth sine harmonics. Different combinations of the respective coefficients refer to the main PAD shapes within the outer radiation belt, namely the pancake, flat-top, butterfly and cap PADs. We conduct a superposed epoch analysis of 129 geomagnetic storms and analyze the PAD evolution for day and night MLT sectors. PAD shapes exhibit a strong energy-dependent response. At energies of tens of keV, the PADs exhibit little variation throughout geomagnetic storms. Cap PADs are mainly observed at energies < 300 keV, and their extent in L shrinks with increasing energy. The cap distributions transform into the pancake PADs around the main phase of the storm on the nightside, and then come back to their original shapes during the recovery phase. At higher energies on the dayside, the PADs are mainly pancake during pre-storm conditions and become more anisotropic during the main phase. The quiet-time butterfly PADs can be observed on the nightside at L> 5.6. During the main phase, butterfly PADs have stronger 90 degrees-minima and can be observed at lower L-shells (down to L = 5), then transitioning into flat-top PADs at L similar to 4.5 - 5 and pancake PADs at L < 4.5. The resulting PAD coefficients for different energies, locations and storm epochs can be used to test the wave models and physics-based radiation belt codes in terms of pitch angle distributions.}, language = {en} } @article{DrozdovAseevEffenbergeretal.2019, author = {Drozdov, Alexander and Aseev, Nikita and Effenberger, Frederic and Turner, Drew L. and Saikin, Anthony and Shprits, Yuri Y.}, title = {Storm Time Depletions of Multi-MeV Radiation Belt Electrons Observed at Different Pitch Angles}, series = {Journal of geophysical research : Space physics}, volume = {124}, journal = {Journal of geophysical research : Space physics}, number = {11}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2019JA027332}, pages = {8943 -- 8953}, year = {2019}, abstract = {During geomagnetic storms, the rapid depletion of the high-energy (several MeV) outer radiation belt electrons is the result of loss to the interplanetary medium through the magnetopause, outward radial diffusion, and loss to the atmosphere due to wave-particle interactions. We have performed a statistical study of 110 storms using pitch angle resolved electron flux measurements from the Van Allen Probes mission and found that inside of the radiation belt (L* = 3 - 5) the number of storms that result in depletion of electrons with equatorial pitch angle alpha(eq) = 30 degrees is higher than number of storms that result in depletion of electrons with equatorial pitch angle alpha(eq) = 75 degrees. We conclude that this result is consistent with electron scattering by whistler and electromagnetic ion cyclotron waves. At the outer edge of the radiation belt (L* >= 5.2) the number of storms that result in depletion is also large (similar to 40-50\%), emphasizing the significance of the magnetopause shadowing effect and outward radial transport.}, language = {en} } @article{BoydSpenceHuangetal.2016, author = {Boyd, A. J. and Spence, Harlan E. and Huang, Chia-Lin and Reeves, Geoffrey D. and Baker, Daniel N. and Turner, D. L. and Claudepierre, Seth G. and Fennell, Joseph F. and Blake, J. Bernard and Shprits, Yuri Y.}, title = {Statistical properties of the radiation belt seed population}, series = {Journal of geophysical research : Space physics}, volume = {121}, journal = {Journal of geophysical research : Space physics}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1002/2016JA022652}, pages = {7636 -- 7646}, year = {2016}, abstract = {We present a statistical analysis of phase space density data from the first 26 months of the Van Allen Probes mission. In particular, we investigate the relationship between the tens and hundreds of keV seed electrons and >1 MeV core radiation belt electron population. Using a cross-correlation analysis, we find that the seed and core populations are well correlated with a coefficient of approximate to 0.73 with a time lag of 10-15 h. We present evidence of a seed population threshold that is necessary for subsequent acceleration. The depth of penetration of the seed population determines the inner boundary of the acceleration process. However, we show that an enhanced seed population alone is not enough to produce acceleration in the higher energies, implying that the seed population of hundreds of keV electrons is only one of several conditions required for MeV electron radiation belt acceleration.}, language = {en} } @article{KimShprits2019, author = {Kim, Kyung-Chan and Shprits, Yuri Y.}, title = {Statistical Analysis of Hiss Waves in Plasmaspheric Plumes Using Van Allen Probe Observations}, series = {Journal of geophysical research : Space physics}, volume = {124}, journal = {Journal of geophysical research : Space physics}, number = {3}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2018JA026458}, pages = {1904 -- 1915}, year = {2019}, abstract = {Plasmaspheric hiss waves commonly observed in high-density regions in the Earth's magnetosphere are known to be one of the main contributors to the loss of radiation belt electrons. There has been a lot of effort to investigate the distributions of hiss waves in the plasmasphere, while relatively little attention has been given to those in the plasmaspheric plume. In this study, we present for the first time a statistical analysis of the occurrence and the spatial distribution of wave amplitudes and wave normal angles for hiss waves in plumes using Van Allen Probes observations during the period of October 2012 to December 2016. Statistical results show that a wide range of hiss wave amplitudes in plumes from a few picotesla to >100 pT is observed, but a modest (<20 pT) wave amplitude is more commonly observed regardless of geomagnetic activity in both the midnight-to-dawn and dusk sector. By contrast, stronger amplitude hiss occurs preferentially during geomagnetically active times in the dusk sector. The wave normal angles are distributed over a broad range from 0° to 90° with a bimodal distribution: a quasi-field-aligned population (<20°) with an occurrence rate of <60\% and an oblique one (>50°) with a relative low occurrence rate of ≲20\%. Therefore, from a statistical point of view, we confirm that the hiss intensity (a few tens of picotesla) and field-aligned hiss wave adopted in previous simulation studies are a reasonable assumption but stress that the activity dependence of the wave amplitude should be considered.}, language = {en} } @article{CastilloShpritsGanushkinaetal.2019, author = {Castillo, Angelica M. and Shprits, Yuri Y. and Ganushkina, Natalia and Drozdov, Alexander and Aseev, Nikita and Wang, Dedong and Dubyagin, Stepan}, title = {Simulations of the inner magnetospheric energetic electrons using the IMPTAM-VERB coupled model}, series = {Journal of Atmospheric and Solar-Terrestrial Physics}, volume = {191}, journal = {Journal of Atmospheric and Solar-Terrestrial Physics}, publisher = {Elsevier}, address = {Oxford}, issn = {1364-6826}, doi = {10.1016/j.jastp.2019.05.014}, pages = {17}, year = {2019}, abstract = {In this study, we present initial results of the coupling between the Inner Magnetospheric Particle Transport and Acceleration Model (IMPTAM) and the Versatile Electron Radiation Belt (VERB-3D) code. IMPTAM traces electrons of 10-100 keV energies from the plasma sheet (L = 9 Re) to inner L-shell regions. The flux evolution modeled by IMPTAM is used at the low energy and outer L* computational boundaries of the VERB code (assuming a dipole approximation) to perform radiation belt simulations of energetic electrons. The model was tested on the March 17th, 2013 storm, for a six-day period. Four different simulations were performed and their results compared to satellites observations from Van Allen probes and GOES. The coupled IMPTAM-VERB model reproduces evolution and storm-time features of electron fluxes throughout the studied storm in agreement with the satellite data (within similar to 0.5 orders of magnitude). Including dynamics of the low energy population at L* = 6.6 increases fluxes closer to the heart of the belt and has a strong impact in the VERB simulations at all energies. However, inclusion of magnetopause losses leads to drastic flux decreases even below L* = 3. The dynamics of low energy electrons (max. 10s of keV) do not affect electron fluxes at energies >= 900 keV. Since the IMPTAM-VERB coupled model is only driven by solar wind parameters and the Dst and Kp indexes, it is suitable as a forecasting tool. In this study, we demonstrate that the estimation of electron dynamics with satellite-data-independent models is possible and very accurate.}, language = {en} } @article{AseevShpritsDrozdovetal.2017, author = {Aseev, Nikita and Shprits, Yuri Y. and Drozdov, Alexander and Kellerman, Adam C. and Usanova, Maria E. and Wang, D. and Zhelavskaya, Irina}, title = {Signatures of Ultrarelativistic Electron Loss in the Heart of the Outer Radiation Belt Measured by Van Allen Probes}, series = {Journal of geophysical research : Space physics}, volume = {122}, journal = {Journal of geophysical research : Space physics}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1002/2017JA024485}, pages = {10102 -- 10111}, year = {2017}, abstract = {Up until recently, signatures of the ultrarelativistic electron loss driven by electromagnetic ion cyclotron (EMIC) waves in the Earth's outer radiation belt have been limited to direct or indirect measurements of electron precipitation or the narrowing of normalized pitch angle distributions in the heart of the belt. In this study, we demonstrate additional observational evidence of ultrarelativistic electron loss that can be driven by resonant interaction with EMIC waves. We analyzed the profiles derived from Van Allen Probe particle data as a function of time and three adiabatic invariants between 9 October and 29 November 2012. New local minimums in the profiles are accompanied by the narrowing of normalized pitch angle distributions and ground\&\#8208;based detection of EMIC waves. Such a correlation may be indicative of ultrarelativistic electron precipitation into the Earth's atmosphere caused by resonance with EMIC waves.}, language = {en} } @article{CaoNiSummersetal.2019, author = {Cao, Xing and Ni, Binbin and Summers, Danny and Shprits, Yuri Y. and Gu, Xudong and Fu, Song and Lou, Yuequn and Zhang, Yang and Ma, Xin and Zhang, Wenxun and Huang, He and Yi, Juan}, title = {Sensitivity of EMIC wave-driven scattering loss of ring current protons to wave normal angle distribution}, series = {Geophysical research letters}, volume = {46}, journal = {Geophysical research letters}, number = {2}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1029/2018GL081550}, pages = {590 -- 598}, year = {2019}, abstract = {Electromagnetic ion cyclotron waves have long been recognized to play a crucial role in the dynamic loss of ring current protons. While the field-aligned propagation approximation of electromagnetic ion cyclotron waves was widely used to quantify the scattering loss of ring current protons, in this study, we find that the wave normal distribution strongly affects the pitch angle scattering efficiency of protons. Increase of peak normal angle or angular width can considerably reduce the scattering rates of <= 10 keV protons. For >10 keV protons, the field-aligned propagation approximation results in a pronounced underestimate of the scattering of intermediate equatorial pitch angle protons and overestimates the scattering of high equatorial pitch angle protons by orders of magnitude. Our results suggest that the wave normal distribution of electromagnetic ion cyclotron waves plays an important role in the pitch angle evolution and scattering loss of ring current protons and should be incorporated in future global modeling of ring current dynamics.}, language = {en} } @misc{ShpritsAngelopoulosRusselletal.2017, author = {Shprits, Yuri Y. and Angelopoulos, V. and Russell, C. T. and Strangeway, R. J. and Runov, A. and Turner, D. and Caron, R. and Cruce, P. and Leneman, D. and Michaelis, I. and Petrov, V. and Panasyuk, M. and Yashin, I. and Drozdov, Alexander and Russell, C. L. and Kalegaev, V. and Nazarkov, I. and Clemmons, J. H.}, title = {Scientific Objectives of Electron Losses and Fields INvestigation Onboard Lomonosov Satellite}, series = {Space science reviews}, volume = {214}, journal = {Space science reviews}, number = {1}, publisher = {Springer}, address = {Dordrecht}, issn = {0038-6308}, doi = {10.1007/s11214-017-0455-4}, pages = {19}, year = {2017}, abstract = {The objective of the Electron Losses and Fields INvestigation on board the Lomonosov satellite ( ELFIN-L) project is to determine the energy spectrum of precipitating energetic electrons and ions and, together with other polar-orbiting and equatorial missions, to better understand the mechanisms responsible for scattering these particles into the atmosphere. This mission will provide detailed measurements of the radiation environment at low altitudes. The 400-500 km sun-synchronous orbit of Lomonosov is ideal for observing electrons and ions precipitating into the atmosphere. This mission provides a unique opportunity to test the instruments. Similar suite of instruments will be flown in the future NSF-and NASA-supported spinning CubeSat ELFIN satellites which will augment current measurements by providing detailed information on pitch-angle distributions of precipitating and trapped particles.}, language = {en} } @article{CaoShpritsNietal.2017, author = {Cao, Xing and Shprits, Yuri Y. and Ni, Binbin and Zhelavskaya, Irina}, title = {Scattering of Ultra-relativistic Electrons in the Van Allen Radiation Belts Accounting for Hot Plasma Effects}, series = {Scientific reports}, volume = {7}, journal = {Scientific reports}, publisher = {Nature Publ. Group}, address = {London}, issn = {2045-2322}, doi = {10.1038/s41598-017-17739-7}, pages = {7}, year = {2017}, abstract = {Electron flux in the Earth's outer radiation belt is highly variable due to a delicate balance between competing acceleration and loss processes. It has been long recognized that Electromagnetic Ion Cyclotron (EMIC) waves may play a crucial role in the loss of radiation belt electrons. Previous theoretical studies proposed that EMIC waves may account for the loss of the relativistic electron population. However, recent observations showed that while EMIC waves are responsible for the significant loss of ultra-relativistic electrons, the relativistic electron population is almost unaffected. In this study, we provide a theoretical explanation for this discrepancy between previous theoretical studies and recent observations. We demonstrate that EMIC waves mainly contribute to the loss of ultra-relativistic electrons. This study significantly improves the current understanding of the electron dynamics in the Earth's radiation belt and also can help us understand the radiation environments of the exoplanets and outer planets.}, language = {en} } @article{SaikinShpritsDrozdovetal.2021, author = {Saikin, Anthony and Shprits, Yuri Y. and Drozdov, Alexander and Landis, Daji August and Zhelavskaya, Irina and Cervantes Villa, Juan Sebastian}, title = {Reconstruction of the radiation belts for solar cycles 17-24 (1933-2017)}, series = {Space weather : the international journal of research and applications}, volume = {19}, journal = {Space weather : the international journal of research and applications}, number = {3}, publisher = {Wiley}, address = {New York}, issn = {1542-7390}, doi = {10.1029/2020SW002524}, pages = {24}, year = {2021}, abstract = {We present a reconstruction of the dynamics of the radiation belts from solar cycles 17 to 24 which allows us to study how radiation belt activity has varied between the different solar cycles. The radiation belt simulations are produced using the Versatile Electron Radiation Belt (VERB)-3D code. The VERB-3D code simulations incorporate radial, energy, and pitch angle diffusion to reproduce the radiation belts. Our simulations use the historical measurements of Kp (available since solar cycle 17, i.e., 1933) to model the evolution radiation belt dynamics between L* = 1-6.6. A nonlinear auto regressive network with exogenous inputs (NARX) neural network was trained off GOES 15 measurements (January 2011-March 2014) and used to supply the upper boundary condition (L* = 6.6) over the course of solar cycles 17-24 (i.e., 1933-2017). Comparison of the model with long term observations of the Van Allen Probes and CRRES demonstrates that our model, driven by the NARX boundary, can reconstruct the general evolution of the radiation belt fluxes. Solar cycle 24 (January 2008-2017) has been the least active of the considered solar cycles which resulted in unusually low electron fluxes. Our results show that solar cycle 24 should not be used as a representative solar cycle for developing long term environment models. The developed reconstruction of fluxes can be used to develop or improve empirical models of the radiation belts.}, language = {en} } @article{AseevShprits2019, author = {Aseev, Nikita and Shprits, Yuri Y.}, title = {Reanalysis of ring current electron phase space densities using van allen probe observations, convection model, and log-normal kalman filter}, series = {Space weather : the international journal of research and applications}, volume = {17}, journal = {Space weather : the international journal of research and applications}, number = {4}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1542-7390}, doi = {10.1029/2018SW002110}, pages = {619 -- 638}, year = {2019}, abstract = {Models of ring current electron dynamics unavoidably contain uncertainties in boundary conditions, electric and magnetic fields, electron scattering rates, and plasmapause location. Model errors can accumulate with time and result in significant deviations of model predictions from observations. Data assimilation offers useful tools which can combine physics-based models and measurements to improve model predictions. In this study, we systematically analyze performance of the Kalman filter applied to a log-transformed convection model of ring current electrons and Van Allen Probe data. We consider long-term dynamics of mu = 2.3 MeV/G and K = 0.3 G(1/2) R-E electrons from 1 February 2013 to 16 June 2013. By using synthetic data, we show that the Kalman filter is capable of correcting errors in model predictions associated with uncertainties in electron lifetimes, boundary conditions, and convection electric fields. We demonstrate that reanalysis retains features which cannot be fully reproduced by the convection model such as storm-time earthward propagation of the electrons down to 2.5 R-E. The Kalman filter can adjust model predictions to satellite measurements even in regions where data are not available. We show that the Kalman filter can adjust model predictions in accordance with observations for mu = 0.1, 2.3, and 9.9 MeV/G and constant K = 0.3 G(1/2) R-E electrons. The results of this study demonstrate that data assimilation can improve performance of ring current models, better quantify model uncertainties, and help deeper understand the physics of the ring current particles.}, language = {en} } @article{DobyndeEffenbergerKartashovetal.2019, author = {Dobynde, M. I. and Effenberger, Frederic and Kartashov, D. A. and Shprits, Yuri Y. and Shurshakov, V. A.}, title = {Ray-tracing simulation of the radiation dose distribution on the surface of the spherical phantom of the MATROSHKA-R experiment onboard the ISS}, series = {Life sciences in space research}, volume = {21}, journal = {Life sciences in space research}, publisher = {Elsevier}, address = {Amsterdam}, issn = {2214-5524}, doi = {10.1016/j.lssr.2019.04.001}, pages = {65 -- 72}, year = {2019}, abstract = {Space radiation is one of the main concerns for human space flights. The prediction of the radiation dose for the actual spacecraft geometry is very important for the planning of long-duration missions. We present a numerical method for the fast calculation of the radiation dose rate during a space flight. We demonstrate its application for dose calculations during the first and the second sessions of the MATROSHKA-R space experiment with a spherical tissue-equivalent phantom. The main advantage of the method is the short simulation time, so it can be applied for urgent radiation dose calculations for low-Earth orbit space missions. The method uses depth-dose curve and shield-and-composition distribution functions to calculate a radiation dose at the point of interest. The spacecraft geometry is processed into a shield-and-composition distribution function using a ray-tracing method. Depth-dose curves are calculated using the GEANT4 Monte-Carlo code (version 10.00.P02) for a double-layer aluminum-water shielding. Aluminum-water shielding is a good approximation of the real geometry, as water is a good equivalent for biological tissues, and aluminum is the major material of spacecraft bodies.}, language = {en} } @article{WoodfieldGlauertMeniettietal.2019, author = {Woodfield, Emma E. and Glauert, Saraha A. and Menietti, J. Douglas and Averkamp, Terrance F. and Horne, Richard B. and Shprits, Yuri Y.}, title = {Rapid Electron Acceleration in Low-Density Regions of Saturn's Radiation Belt by Whistler Mode Chorus Waves}, series = {Geophysical research letters}, volume = {46}, journal = {Geophysical research letters}, number = {13}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1029/2019GL083071}, pages = {7191 -- 7198}, year = {2019}, abstract = {Electron acceleration at Saturn due to whistler mode chorus waves has previously been assumed to be ineffective; new data closer to the planet show it can be very rapid (factor of 104 flux increase at 1 MeV in 10 days compared to factor of 2). A full survey of chorus waves at Saturn is combined with an improved plasma density model to show that where the plasma frequency falls below the gyrofrequency additional strong resonances are observed favoring electron acceleration. This results in strong chorus acceleration between approximately 2.5 R-S and 5.5 R-S outside which adiabatic transport may dominate. Strong pitch angle dependence results in butterfly pitch angle distributions that flatten over a few days at 100s keV, tens of days at MeV energies which may explain observations of butterfly distributions of MeV electrons near L = 3. Including cross terms in the simulations increases the tendency toward butterfly distributions. Plain Language Summary Radiation belts are hazardous regions found around several of the planets in our Solar System. They consist of very hot, electrically charged particles trapped in the magnetic field of the planet. At Saturn the most important way to heat these particles has for many years been thought to involve the particles drifting closer toward the planet. This paper adds to the emerging idea at Saturn that a different way to heat the particles is also possible where the heating is done by waves, in a similar way to what we find at the Earth. We use recent information from the Cassini spacecraft on the number and location of particles and also of the waves strength and location combined with computer simulations to show that a particular wave called chorus is excellent at heating the particles where the surrounding number of cold particles is low.}, language = {en} } @article{CervantesVillaShpritsAseevetal.2020, author = {Cervantes Villa, Juan Sebastian and Shprits, Yuri Y. and Aseev, Nikita and Allison, Hayley J.}, title = {Quantifying the effects of EMIC wave scattering and magnetopause shadowing in the outer electron radiation belt by means of data assimilation}, series = {Journal of geophysical research : Space physics}, volume = {125}, journal = {Journal of geophysical research : Space physics}, number = {8}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2020JA028208}, pages = {23}, year = {2020}, abstract = {In this study we investigate two distinct loss mechanisms responsible for the rapid dropouts of radiation belt electrons by assimilating data from Van Allen Probes A and B and Geostationary Operational Environmental Satellites (GOES) 13 and 15 into a 3-D diffusion model. In particular, we examine the respective contribution of electromagnetic ion cyclotron (EMIC) wave scattering and magnetopause shadowing for values of the first adiabatic invariant mu ranging from 300 to 3,000 MeV G(-1). We inspect the innovation vector and perform a statistical analysis to quantitatively assess the effect of both processes as a function of various geomagnetic indices, solar wind parameters, and radial distance from the Earth. Our results are in agreement with previous studies that demonstrated the energy dependence of these two mechanisms. We show that EMIC wave scattering tends to dominate loss at lower L shells, and it may amount to between 10\%/hr and 30\%/hr of the maximum value of phase space density (PSD) over all L shells for fixed first and second adiabatic invariants. On the other hand, magnetopause shadowing is found to deplete electrons across all energies, mostly at higher L shells, resulting in loss from 50\%/hr to 70\%/hr of the maximum PSD. Nevertheless, during times of enhanced geomagnetic activity, both processes can operate beyond such location and encompass the entire outer radiation belt.}, language = {en} } @article{DentonOfmanShpritsetal.2019, author = {Denton, Richard E. and Ofman, L. and Shprits, Yuri Y. and Bortnik, J. and Millan, R. M. and Rodger, C. J. and da Silva, C. L. and Rogers, B. N. and Hudson, M. K. and Liu, K. and Min, K. and Glocer, A. and Komar, C.}, title = {Pitch Angle Scattering of Sub-MeV Relativistic Electrons by Electromagnetic Ion Cyclotron Waves}, series = {Journal of geophysical research : Space physics}, volume = {124}, journal = {Journal of geophysical research : Space physics}, number = {7}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9402}, doi = {10.1029/2018JA026384}, pages = {5610 -- 5626}, year = {2019}, abstract = {Electromagnetic ion cyclotron (EMIC) waves have long been considered to be a significant loss mechanism for relativistic electrons. This has most often been attributed to resonant interactions with the highest amplitude waves. But recent observations have suggested that the dominant energy of electrons precipitated to the atmosphere may often be relatively low, less than 1 MeV, whereas the minimum resonant energy of the highest amplitude waves is often greater than 2 MeV. Here we use relativistic electron test particle simulations in the wavefields of a hybrid code simulation of EMIC waves in dipole geometry in order to show that significant pitch angle scattering can occur due to interaction with low-amplitude short-wavelength EMIC waves. In the case we examined, these waves are in the H band (at frequencies above the He+ gyrofrequency), even though the highest amplitude waves were in the He band frequency range (below the He+ gyrofrequency). We also present wave power distributions for 29 EMIC simulations in straight magnetic field line geometry that show that the high wave number portion of the spectrum is in every case mostly due to the H band waves. Though He band waves are often associated with relativistic electron precipitation, it is possible that the He band waves do not directly scatter the sub-megaelectron volts (sub-MeV) electrons, but that the presence of He band waves is associated with high plasma density which lowers the minimum resonant energy so that these electrons can more easily resonate with the H band waves.}, language = {en} } @article{RipollLoridanCunninghametal.2016, author = {Ripoll, Jean-Fran{\c{c}}ois and Loridan, Vivien and Cunningham, G. S. and Reeves, Geoffrey D. and Shprits, Yuri Y.}, title = {On the time needed to reach an equilibrium structure of the radiation belts}, series = {Journal of geophysical research : Space physics}, volume = {121}, journal = {Journal of geophysical research : Space physics}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1002/2015JA022207}, pages = {7684 -- 7698}, year = {2016}, abstract = {In this study, we complement the notion of equilibrium states of the radiation belts with a discussion on the dynamics and time needed to reach equilibrium. We solve for the equilibrium states obtained using 1-D radial diffusion with recently developed hiss and chorus lifetimes at constant values of Kp = 1, 3, and 6. We find that the equilibrium states at moderately low Kp, when plotted versus L shell (L) and energy (E), display the same interesting S shape for the inner edge of the outer belt as recently observed by the Van Allen Probes. The S shape is also produced as the radiation belts dynamically evolve toward the equilibrium state when initialized to simulate the buildup after a massive dropout or to simulate loss due to outward diffusion from a saturated state. Physically, this shape, intimately linked with the slot structure, is due to the dependence of electron loss rate (originating from wave-particle interactions) on both energy and L shell. Equilibrium electron flux profiles are governed by the Biot number (tau(Diffusion)/tau(loss)), with large Biot number corresponding to low fluxes and low Biot number to large fluxes. The time it takes for the flux at a specific (L, E) to reach the value associated with the equilibrium state, starting from these different initial states, is governed by the initial state of the belts, the property of the dynamics (diffusion coefficients), and the size of the domain of computation. Its structure shows a rather complex scissor form in the (L, E) plane. The equilibrium value (phase space density or flux) is practically reachable only for selected regions in (L, E) and geomagnetic activity. Convergence to equilibrium requires hundreds of days in the inner belt for E>300 keV and moderate Kp (<= 3). It takes less time to reach equilibrium during disturbed geomagnetic conditions (Kp = 3), when the system evolves faster. Restricting our interest to the slot region, below L = 4, we find that only small regions in (L, E) space can reach the equilibrium value: E similar to [200, 300] keV for L= [3.7, 4] at Kp= 1, E similar to[0.6, 1] MeV for L = [3, 4] at Kp = 3, and E similar to 300 keV for L = [3.5, 4] at Kp = 6 assuming no new incoming electrons.}, language = {en} } @article{WangShprits2019, author = {Wang, Dedong and Shprits, Yuri Y.}, title = {On How High-Latitude Chorus Waves Tip the Balance Between Acceleration and Loss of Relativistic Electrons}, series = {Geophysical research letters}, volume = {46}, journal = {Geophysical research letters}, number = {14}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1029/2019GL082681}, pages = {7945 -- 7954}, year = {2019}, abstract = {Modeling and observations have shown that energy diffusion by chorus waves is an important source of acceleration of electrons to relativistic energies. By performing long-term simulations using the three-dimensional Versatile Electron Radiation Belt code, in this study, we test how the latitudinal dependence of chorus waves can affect the dynamics of the radiation belt electrons. Results show that the variability of chorus waves at high latitudes is critical for modeling of megaelectron volt (MeV) electrons. We show that, depending on the latitudinal distribution of chorus waves under different geomagnetic conditions, they cannot only produce a net acceleration but also a net loss of MeV electrons. Decrease in high-latitude chorus waves can tip the balance between acceleration and loss toward acceleration, or alternatively, the increase in high-latitude waves can result in a net loss of MeV electrons. Variations in high-latitude chorus may account for some of the variability of MeV electrons.}, language = {en} } @article{RipollLoridanDentonetal.2019, author = {Ripoll, Jean-Francois and Loridan, Vivien and Denton, Michael H. and Cunningham, Gregory and Reeves, G. and Santolik, O. and Fennell, Joseph and Turner, Drew L. and Drozdov, Alexander and Cervantes Villa, Juan Sebastian and Shprits, Yuri Y. and Thaller, Scott A. and Kurth, William S. and Kletzing, Craig A. and Henderson, Michael G. and Ukhorskiy, Aleksandr Y.}, title = {Observations and Fokker-Planck Simulations of the L-Shell, Energy, and Times}, series = {Journal of geophysical research : Space physics}, volume = {124}, journal = {Journal of geophysical research : Space physics}, number = {2}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2018JA026111}, pages = {1125 -- 1142}, year = {2019}, abstract = {The evolution of the radiation belts in L-shell (L), energy (E), and equatorial pitch angle (alpha(0)) is analyzed during the calm 11-day interval (4-15 March) following the 1 March 2013 storm. Magnetic Electron and Ion Spectrometer (MagEIS) observations from Van Allen Probes are interpreted alongside 1D and 3D Fokker-Planck simulations combined with consistent event-driven scattering modeling from whistler mode hiss waves. Three (L, E, alpha(0)) regions persist through 11 days of hiss wave scattering; the pitch angle-dependent inner belt core (L similar to <2.2 and E < 700 keV), pitch angle homogeneous outer belt low-energy core (L > similar to 5 and E similar to < 100 keV), and a distinct pocket of electrons (L similar to [4.5, 5.5] and E similar to [0.7, 2] MeV). The pitch angle homogeneous outer belt is explained by the diffusion coefficients that are roughly constant for alpha(0) similar to <60 degrees, E > 100 keV, 3.5 < L < L-pp similar to 6. Thus, observed unidirectional flux decays can be used to estimate local pitch angle diffusion rates in that region. Top-hat distributions are computed and observed at L similar to 3-3.5 and E = 100-300 keV.}, language = {en} } @article{AseevShpritsDrozdovetal.2016, author = {Aseev, Nikita and Shprits, Yuri Y. and Drozdov, Alexander and Kellerman, Adam C.}, title = {Numerical applications of the advective-diffusive codes for the inner magnetosphere}, series = {Space Weather: The International Journal of Research and Applications}, volume = {14}, journal = {Space Weather: The International Journal of Research and Applications}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1542-7390}, doi = {10.1002/2016SW001484}, pages = {993 -- 1010}, year = {2016}, abstract = {In this study we present analytical solutions for convection and diffusion equations. We gather here the analytical solutions for the one-dimensional convection equation, the two-dimensional convection problem, and the one- and two-dimensional diffusion equations. Using obtained analytical solutions, we test the four-dimensional Versatile Electron Radiation Belt code (the VERB-4D code), which solves the modified Fokker-Planck equation with additional convection terms. The ninth-order upwind numerical scheme for the one-dimensional convection equation shows much more accurate results than the results obtained with the third-order scheme. The universal limiter eliminates unphysical oscillations generated by high-order linear upwind schemes. Decrease in the space step leads to convergence of a numerical solution of the two-dimensional diffusion equation with mixed terms to the analytical solution. We compare the results of the third- and ninth-order schemes applied to magnetospheric convection modeling. The results show significant differences in electron fluxes near geostationary orbit when different numerical schemes are used.}, language = {en} } @article{ShpritsVasileZhelayskaya2019, author = {Shprits, Yuri Y. and Vasile, Ruggero and Zhelayskaya, Irina S.}, title = {Nowcasting and Predicting the Kp Index Using Historical Values and Real-Time Observations}, series = {Space Weather: The International Journal of Research and Applications}, volume = {17}, journal = {Space Weather: The International Journal of Research and Applications}, number = {8}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1542-7390}, doi = {10.1029/2018SW002141}, pages = {1219 -- 1229}, year = {2019}, abstract = {Current algorithms for the real-time prediction of the Kp index use a combination of models empirically driven by solar wind measurements at the L1 Lagrange point and historical values of the index. In this study, we explore the limitations of this approach, examining the forecast for short and long lead times using measurements at L1 and Kp time series as input to artificial neural networks. We explore the relative efficiency of the solar wind-based predictions, predictions based on recurrence, and predictions based on persistence. Our modeling results show that for short-term forecasts of approximately half a day, the addition of the historical values of Kp to the measured solar wind values provides a barely noticeable improvement. For a longer-term forecast of more than 2 days, predictions can be made using recurrence only, while solar wind measurements provide very little improvement for a forecast with long horizon times. We also examine predictions for disturbed and quiet geomagnetic activity conditions. Our results show that the paucity of historical measurements of the solar wind for high Kp results in a lower accuracy of predictions during disturbed conditions. Rebalancing of input data can help tailor the predictions for more disturbed conditions.}, language = {en} } @article{ZhuShpritsSpasojevicetal.2019, author = {Zhu, Hui and Shprits, Yuri Y. and Spasojevic, M. and Drozdov, Alexander}, title = {New hiss and chorus waves diffusion coefficient parameterizations from the Van Allen Probes and their effect on long-term relativistic electron radiation-belt VERB simulations}, series = {Journal of Atmospheric and Solar-Terrestrial Physics}, volume = {193}, journal = {Journal of Atmospheric and Solar-Terrestrial Physics}, publisher = {Elsevier}, address = {Oxford}, issn = {1364-6826}, doi = {10.1016/j.jastp.2019.105090}, pages = {13}, year = {2019}, abstract = {New wave frequency and amplitude models for the nightside and dayside chorus waves are built based on measurements from the Electric and Magnetic Field Instrument Suite and Integrated Science (EMFISIS) instrument onboard the Van Allen Probes. The corresponding 3D diffusion coefficients are systematically obtained. Compared with previous commonly-used (typical) parameterizations, the new parameterizations result in differences in diffusion rates that depend on the energy and pitch angle. Furthermore, one-year 3D diffusive simulations are performed using the Versatile Electron Radiation Belt (VERB) code. Both typical and new wave parameterizations simulation results are in a good agreement with observations at 0.9 MeV. However, the new parameterizations for nightside chorus better reproduce the observed electron fluxes. These parameterizations will be incorporated into future modeling efforts.}, language = {en} } @article{LandisSaikinZhelavskayaetal.2022, author = {Landis, Daji August and Saikin, Anthony and Zhelavskaya, Irina and Drozdov, Alexander and Aseev, Nikita and Shprits, Yuri Y. and Pfitzer, Maximilian F. and Smirnov, Artem G.}, title = {NARX Neural Network Derivations of the Outer Boundary Radiation Belt Electron Flux}, series = {Space Weather: the international journal of research and applications}, volume = {20}, journal = {Space Weather: the international journal of research and applications}, number = {5}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1542-7390}, doi = {10.1029/2021SW002774}, pages = {18}, year = {2022}, abstract = {We present two new empirical models of radiation belt electron flux at geostationary orbit. GOES-15 measurements of 0.8 MeV electrons were used to train a Nonlinear Autoregressive with Exogenous input (NARX) neural network for both modeling GOES-15 flux values and an upper boundary condition scaling factor (BF). The GOES-15 flux model utilizes an input and feedback delay of 2 and 2 time steps (i.e., 5 min time steps) with the most efficient number of hidden layers set to 10. Magnetic local time, Dst, Kp, solar wind dynamic pressure, AE, and solar wind velocity were found to perform as predicative indicators of GOES-15 flux and therefore were used as the exogenous inputs. The NARX-derived upper boundary condition scaling factor was used in conjunction with the Versatile Electron Radiation Belt (VERB) code to produce reconstructions of the radiation belts during the period of July-November 1990, independent of in-situ observations. Here, Kp was chosen as the sole exogenous input to be more compatible with the VERB code. This Combined Release and Radiation Effects Satellite-era reconstruction showcases the potential to use these neural network-derived boundary conditions as a method of hindcasting the historical radiation belts. This study serves as a companion paper to another recently published study on reconstructing the radiation belts during Solar Cycles 17-24 (Saikin et al., 2021, ), for which the results featured in this paper were used.}, language = {en} } @article{ShpritsKellermanAseevetal.2017, author = {Shprits, Yuri Y. and Kellerman, Adam C . and Aseev, Nikita and Drozdov, Alexander and Michaelis, Ingo}, title = {Multi-MeV electron loss in the heart of the radiation belts}, series = {Geophysical research letters}, volume = {44}, journal = {Geophysical research letters}, number = {3}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1002/2016GL072258}, pages = {1204 -- 1209}, year = {2017}, abstract = {Significant progress has been made in recent years in understanding acceleration mechanisms in the Earth's radiation belts. In particular, a number of studies demonstrated the importance of the local acceleration by analyzing the radial profiles of phase space density (PSD) and observing building up peaks in PSD. In this study, we focus on understanding of the local loss using very similar tools. The profiles of PSD for various values of the first adiabatic invariants during the previously studied 17 January 2013 storm are presented and discussed. The profiles of PSD show clear deepening minimums consistent with the scattering by electromagnetic ion cyclotron waves. Long-term evolution shows that local minimums in PSD can persist for relatively long times. During considered interval of time the deepening minimums were observed around L* = 4 during 17 January 2013 storm and around L* = 3.5 during 1 March 2013 storm. This study shows a new method that can help identify the location, magnitude, and time of the local loss and will help quantify local loss in the future. This study also provides additional clear and definitive evidence that local loss plays a major role for the dynamics of the multi-MeV electrons.}, language = {en} } @article{ZhuChenLiuetal.2019, author = {Zhu, Hui and Chen, Lunjin and Liu, Xu and Shprits, Yuri Y.}, title = {Modulation of locally generated equatorial noise by ULF wave}, series = {Journal of geophysical research : Space physics}, volume = {124}, journal = {Journal of geophysical research : Space physics}, number = {4}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2018JA026199}, pages = {2779 -- 2787}, year = {2019}, abstract = {In this paper we report a rare and fortunate event of fast magnetosonic (MS, also called equatorial noise) waves modulated by compressional ultralow frequency (ULF) waves measured by Van Allen Probes. The characteristics of MS waves, ULF waves, proton distribution, and their potential correlations are analyzed. The results show that ULF waves can modulate the energetic ring proton distribution and in turn modulate the MS generation. Furthermore, the variation of MS intensities is attributed to not only ULF wave activities but also the variation of background parameters, for example, number density. The results confirm the opinion that MS waves are generated by proton ring distribution and propose a new modulation phenomenon.}, language = {en} } @article{LiemohnMcColloughJordanovaetal.2018, author = {Liemohn, Michael W. and McCollough, James P. and Jordanova, Vania K. and Ngwira, Chigomezyo M. and Morley, Steven K. and Cid, Consuelo and Tobiska, W. Kent and Wintoft, Peter and Ganushkina, Natalia Yu and Welling, Daniel T. and Bingham, Suzy and Balikhin, Michael A. and Opgenoorth, Hermann J. and Engel, Miles A. and Weigel, Robert S. and Singer, Howard J. and Buresova, Dalia and Bruinsma, Sean and Zhelavskaya, Irina and Shprits, Yuri Y. and Vasile, Ruggero}, title = {Model Evaluation Guidelines for Geomagnetic Index Predictions}, series = {Space Weather: The International Journal of Research and Applications}, volume = {16}, journal = {Space Weather: The International Journal of Research and Applications}, number = {12}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1542-7390}, doi = {10.1029/2018SW002067}, pages = {2079 -- 2102}, year = {2018}, abstract = {Geomagnetic indices are convenient quantities that distill the complicated physics of some region or aspect of near-Earth space into a single parameter. Most of the best-known indices are calculated from ground-based magnetometer data sets, such as Dst, SYM-H, Kp, AE, AL, and PC. Many models have been created that predict the values of these indices, often using solar wind measurements upstream from Earth as the input variables to the calculation. This document reviews the current state of models that predict geomagnetic indices and the methods used to assess their ability to reproduce the target index time series. These existing methods are synthesized into a baseline collection of metrics for benchmarking a new or updated geomagnetic index prediction model. These methods fall into two categories: (1) fit performance metrics such as root-mean-square error and mean absolute error that are applied to a time series comparison of model output and observations and (2) event detection performance metrics such as Heidke Skill Score and probability of detection that are derived from a contingency table that compares model and observation values exceeding (or not) a threshold value. A few examples of codes being used with this set of metrics are presented, and other aspects of metrics assessment best practices, limitations, and uncertainties are discussed, including several caveats to consider when using geomagnetic indices. Plain Language Summary One aspect of space weather is a magnetic signature across the surface of the Earth. The creation of this signal involves nonlinear interactions of electromagnetic forces on charged particles and can therefore be difficult to predict. The perturbations that space storms and other activity causes in some observation sets, however, are fairly regular in their pattern. Some of these measurements have been compiled together into a single value, a geomagnetic index. Several such indices exist, providing a global estimate of the activity in different parts of geospace. Models have been developed to predict the time series of these indices, and various statistical methods are used to assess their performance at reproducing the original index. Existing studies of geomagnetic indices, however, use different approaches to quantify the performance of the model. This document defines a standardized set of statistical analyses as a baseline set of comparison tools that are recommended to assess geomagnetic index prediction models. It also discusses best practices, limitations, uncertainties, and caveats to consider when conducting a model assessment.}, language = {en} } @article{SmirnovBerrendorfShpritsetal.2020, author = {Smirnov, Artem and Berrendorf, Max and Shprits, Yuri Y. and Kronberg, Elena A. and Allison, Hayley J. and Aseev, Nikita and Zhelavskaya, Irina and Morley, Steven K. and Reeves, Geoffrey D. and Carver, Matthew R. and Effenberger, Frederic}, title = {Medium energy electron flux in earth's outer radiation belt (MERLIN)}, series = {Space weather : the international journal of research and applications}, volume = {18}, journal = {Space weather : the international journal of research and applications}, number = {11}, publisher = {American geophysical union, AGU}, address = {Washington}, issn = {1542-7390}, doi = {10.1029/2020SW002532}, pages = {20}, year = {2020}, abstract = {The radiation belts of the Earth, filled with energetic electrons, comprise complex and dynamic systems that pose a significant threat to satellite operation. While various models of electron flux both for low and relativistic energies have been developed, the behavior of medium energy (120-600 keV) electrons, especially in the MEO region, remains poorly quantified. At these energies, electrons are driven by both convective and diffusive transport, and their prediction usually requires sophisticated 4D modeling codes. In this paper, we present an alternative approach using the Light Gradient Boosting (LightGBM) machine learning algorithm. The Medium Energy electRon fLux In Earth's outer radiatioN belt (MERLIN) model takes as input the satellite position, a combination of geomagnetic indices and solar wind parameters including the time history of velocity, and does not use persistence. MERLIN is trained on >15 years of the GPS electron flux data and tested on more than 1.5 years of measurements. Tenfold cross validation yields that the model predicts the MEO radiation environment well, both in terms of dynamics and amplitudes o f flux. Evaluation on the test set shows high correlation between the predicted and observed electron flux (0.8) and low values of absolute error. The MERLIN model can have wide space weather applications, providing information for the scientific community in the form of radiation belts reconstructions, as well as industry for satellite mission design, nowcast of the MEO environment, and surface charging analysis.}, language = {en} } @article{LongNiCaoetal.2022, author = {Long, Minyi and Ni, Binbin and Cao, Xing and Gu, Xudong and Kollmann, Peter and Luo, Qiong and Zhou, Ruoxian and Guo, Yingjie and Guo, Deyu and Shprits, Yuri Y.}, title = {Losses of radiation belt energetic particles by encounters with four of the inner Moons of Jupiter}, series = {Journal of geophysical research, Planets}, volume = {127}, journal = {Journal of geophysical research, Planets}, number = {2}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9097}, doi = {10.1029/2021JE007050}, pages = {13}, year = {2022}, abstract = {Based on an improved model of the moon absorption of Jovian radiation belt particles, we investigate quantitatively and comprehensively the absorption probabilities and particle lifetimes due to encounters with four of the inner moons of Jupiter (Amalthea, Thebe, Io, and Europa) inside L < 10. Our results demonstrate that the resultant average lifetimes of energetic protons and electrons vary dramatically between similar to 0.1 days and well above 1,000 days, showing a strong dependence on the particle equatorial pitch angle, kinetic energy and moon orbit. The average lifetimes of energetic protons and electrons against moon absorption are shortest for Io (i.e., similar to 0.1-10 days) and longest for Thebe (i.e., up to thousands of days), with the lifetimes in between for Europa and Amalthea. Due to the diploe tilt angle absorption effect, the average lifetimes of energetic protons and electrons vary markedly below and above alpha eq \${\alpha }_{\mathrm{e}\mathrm{q}}\$ = 67 degrees. Overall, the average electron lifetimes exhibit weak pitch angle dependence, but the average proton lifetimes are strongly dependent on equatorial pitch angle. The average lifetimes of energetic protons decrease monotonically and substantially with the kinetic energy, but the average lifetimes of energetic electrons are roughly constant at energies 5MeV electrons. The most crucial differences introduced by the hot plasma effects occur for >3MeV electron scattering rates by He+ band EMIC waves. Mainly due to the changes of resonant frequency and wave group velocity when the hot protons are included, the difference in scattering rates can be up to an order of magnitude, showing a strong dependence on both electron energy and equatorial pitch angle. Our study confirms the importance of including hot plasma effects in modeling the scattering of ultra-relativistic radiation belt electrons by EMIC waves.}, language = {en} } @article{YamazakiMatzkaStolleetal.2022, author = {Yamazaki, Yosuke and Matzka, J{\"u}rgen and Stolle, Claudia and Kervalishvili, Guram N. and Rauberg, Jan and Bronkalla, Oliver and Morschhauser, Achim and Bruinsma, Sean L. and Shprits, Yuri Y. and Jackson, David R.}, title = {Geomagnetic activity index Hpo}, series = {Geophysical research letters}, volume = {49}, journal = {Geophysical research letters}, number = {10}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1029/2022GL098860}, pages = {9}, year = {2022}, abstract = {The geomagnetic activity index Kp is widely used but is restricted by low time resolution (3-hourly) and an upper limit. To address this, new geomagnetic activity indices, Hpo, are introduced. Similar to Kp, Hpo expresses the level of planetary geomagnetic activity in units of thirds (0o, 0+, 1-, 1o, 1+, 2-, horizontal ellipsis ) based on the magnitude of geomagnetic disturbances observed at subauroral observatories. Hpo has a higher time resolution than Kp. 30-min (Hp30) and 60-min (Hp60) indices are produced. The frequency distribution of Hpo is designed to be similar to that of Kp so that Hpo may be used as a higher time-resolution alternative to Kp. Unlike Kp, which is capped at 9o, Hpo is an open-ended index and thus can characterize severe geomagnetic storms more accurately. Hp30, Hp60 and corresponding linearly scaled ap30 and ap60 are available, in near real time, at the GFZ website (https://www.gfz-potsdam.de/en/hpo-index).}, language = {en} } @misc{WoodfieldHorneGlauertetal.2018, author = {Woodfield, Emma E. and Horne, Richard B. and Glauert, Sarah A. and Menietti, John D. and Shprits, Yuri Y. and Kurth, William S.}, title = {Formation of electron radiation belts at Saturn by Z-mode wave acceleration}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1032}, issn = {1866-8372}, doi = {10.25932/publishup-46834}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-468342}, pages = {9}, year = {2018}, abstract = {At Saturn electrons are trapped in the planet's magnetic field and accelerated to relativistic energies to form the radiation belts, but how this dramatic increase in electron energy occurs is still unknown. Until now the mechanism of radial diffusion has been assumed but we show here that in-situ acceleration through wave particle interactions, which initial studies dismissed as ineffectual at Saturn, is in fact a vital part of the energetic particle dynamics there. We present evidence from numerical simulations based on Cassini spacecraft data that a particular plasma wave, known as Z-mode, accelerates electrons to MeV energies inside 4 R-S (1 R-S = 60,330 km) through a Doppler shifted cyclotron resonant interaction. Our results show that the Z-mode waves observed are not oblique as previously assumed and are much better accelerators than O-mode waves, resulting in an electron energy spectrum that closely approaches observed values without any transport effects included.}, language = {en} } @article{WoodfieldHorneGlauertetal.2018, author = {Woodfield, Emma E. and Horne, Richard B. and Glauert, S. A. and Menietti, J. D. and Shprits, Yuri Y. and Kurth, William S.}, title = {Formation of electron radiation belts at Saturn by Z-mode wave acceleration}, series = {Nature Communications}, volume = {9}, journal = {Nature Communications}, publisher = {Nature Publ. Group}, address = {London}, issn = {2041-1723}, doi = {10.1038/s41467-018-07549-4}, pages = {7}, year = {2018}, abstract = {At Saturn electrons are trapped in the planet's magnetic field and accelerated to relativistic energies to form the radiation belts, but how this dramatic increase in electron energy occurs is still unknown. Until now the mechanism of radial diffusion has been assumed but we show here that in-situ acceleration through wave particle interactions, which initial studies dismissed as ineffectual at Saturn, is in fact a vital part of the energetic particle dynamics there. We present evidence from numerical simulations based on Cassini spacecraft data that a particular plasma wave, known as Z-mode, accelerates electrons to MeV energies inside 4 RS (1 RS = 60,330 km) through a Doppler shifted cyclotron resonant interaction. Our results show that the Z-mode waves observed are not oblique as previously assumed and are much better accelerators than O-mode waves, resulting in an electron energy spectrum that closely approaches observed values without any transport effects included.}, language = {en} } @article{WalkerBoyntonShpritsetal.2022, author = {Walker, Simon N. and Boynton, Richard J. and Shprits, Yuri Y. and Balikhin, Michael A. and Drozdov, Alexander}, title = {Forecast of the energetic electron environment of the radiation belts}, series = {Space Weather: The International Journal of Research and Applications}, volume = {20}, journal = {Space Weather: The International Journal of Research and Applications}, number = {12}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1542-7390}, doi = {10.1029/2022SW003124}, pages = {21}, year = {2022}, abstract = {Different modeling methodologies possess different strengths and weakness. For instance, data based models may provide superior accuracy but have a limited spatial coverage while physics based models may provide lower accuracy but provide greater spatial coverage. This study investigates the coupling of a data based model of the electron fluxes at geostationary orbit (GEO) with a numerical model of the radiation belt region to improve the resulting forecasts/pastcasts of electron fluxes over the whole radiation belt region. In particular, two coupling methods are investigated. The first assumes an average value for L* for GEO, namely LGEO* L-GEO* = 6.2. The second uses a value of L* that varies with geomagnetic activity, quantified using the Kp index. As the terrestrial magnetic field responds to variations in geomagnetic activity, the value of L* will vary for a specific location. In this coupling method, the value of L* is calculated using the Kp driven Tsyganenko 89c magnetic field model for field line tracing. It is shown that this addition can result in changes in the initialization of the parameters at the Versatile Electron Radiation Belt model outer boundary. Model outputs are compared to Van Allen Probes MagEIS measurements of the electron fluxes in the inner magnetosphere for the March 2015 geomagnetic storm. It is found that the fixed LGEO* L-GEO* coupling method produces a more realistic forecast.}, language = {en} } @article{KimShpritsBlake2016, author = {Kim, Kyung-Chan and Shprits, Yuri Y. and Blake, J. Bernard}, title = {Fast injection of the relativistic electrons into the inner zone and the formation of the split-zone structure during the Bastille Day storm in July 2000}, series = {Journal of geophysical research : Space physics}, volume = {121}, journal = {Journal of geophysical research : Space physics}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1002/2015JA022072}, pages = {8329 -- 8342}, year = {2016}, abstract = {During the July 2000 geomagnetic storm, known as the Bastille Day storm, Solar, Anomalous, and Magnetospheric Particle Explorer (SAMPEX)/Heavy Ion Large Telescope (HILT) observed a strong injection of similar to 1MeV electrons into the slot region (L similar to 2.5) during the storm main phase. Then, during the following month, electrons were clearly seen diffusing inward down to L=2 and forming a pronounced split structure encompassing a narrow, newly formed slot region around L=3. SAMPEX observations are first compared with electron and proton observations on HEO-3 and NOAA-15 to validate that the observed unusual dynamics was not caused by proton contamination of the SAMPEX instrument. The time-dependent 3-D Versatile Electron Radiation Belt (VERB) simulation of 1MeV electron flux evolution is compared with the SAMPEX/HILT observations. The results show that the VERB code predicts overall time evolution of the observed split structure. The simulated split structure is produced by pitch angle scattering into the Earth atmosphere of similar to 1MeV electrons by plasmaspheric hiss.}, language = {en} } @article{ZhelavskayaShpritsSpasojevic2017, author = {Zhelavskaya, Irina and Shprits, Yuri Y. and Spasojevic, Maria}, title = {Empirical Modeling of the Plasmasphere Dynamics Using Neural Networks}, series = {Journal of geophysical research : Space physics}, volume = {122}, journal = {Journal of geophysical research : Space physics}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1002/2017JA024406}, pages = {11227 -- 11244}, year = {2017}, abstract = {We present the PINE (Plasma density in the Inner magnetosphere Neural network\&\#8208;based Empirical) model \&\#8208; a new empirical model for reconstructing the global dynamics of the cold plasma density distribution based only on solar wind data and geomagnetic indices. Utilizing the density database obtained using the NURD (Neural\&\#8208;network\&\#8208;based Upper hybrid Resonance Determination) algorithm for the period of 1 October 2012 to 1 July 2016, in conjunction with solar wind data and geomagnetic indices, we develop a neural network model that is capable of globally reconstructing the dynamics of the cold plasma density distribution for 2\&\#8804;L\&\#8804;6 and all local times. We validate and test the model by measuring its performance on independent data sets withheld from the training set and by comparing the model\&\#8208;predicted global evolution with global images of He+ distribution in the Earth's plasmasphere from the IMAGE Extreme UltraViolet (EUV) instrument. We identify the parameters that best quantify the plasmasphere dynamics by training and comparing multiple neural networks with different combinations of input parameters (geomagnetic indices, solar wind data, and different durations of their time history). The optimal model is based on the 96\&\#8201;h time history of Kp, AE, SYM\&\#8208;H, and F10.7 indices. The model successfully reproduces erosion of the plasmasphere on the nightside and plume formation and evolution. We demonstrate results of both local and global plasma density reconstruction. This study illustrates how global dynamics can be reconstructed from local in situ observations by using machine learning techniques.}, language = {en} } @article{DrozdovShpritsUsanovaetal.2017, author = {Drozdov, Alexander and Shprits, Yuri Y. and Usanova, Maria E. and Aseev, Nikita and Kellerman, Adam C. and Zhu, H.}, title = {EMIC wave parameterization in the long-term VERB code simulation}, series = {Journal of geophysical research : Space physics}, volume = {122}, journal = {Journal of geophysical research : Space physics}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1002/2017JA024389}, pages = {8488 -- 8501}, year = {2017}, abstract = {Electromagnetic ion cyclotron (EMIC) waves play an important role in the dynamics of ultrarelativistic electron population in the radiation belts. However, as EMIC waves are very sporadic, developing a parameterization of such wave properties is a challenging task. Currently, there are no dynamic, activity-dependent models of EMIC waves that can be used in the long-term (several months) simulations, which makes the quantitative modeling of the radiation belt dynamics incomplete. In this study, we investigate Kp, Dst, and AE indices, solar wind speed, and dynamic pressure as possible parameters of EMIC wave presence. The EMIC waves are included in the long-term simulations (1year, including different geomagnetic activity) performed with the Versatile Electron Radiation Belt code, and we compare results of the simulation with the Van Allen Probes observations. The comparison shows that modeling with EMIC waves, parameterized by solar wind dynamic pressure, provides a better agreement with the observations among considered parameterizations. The simulation with EMIC waves improves the dynamics of ultrarelativistic fluxes and reproduces the formation of the local minimum in the phase space density profiles.}, language = {en} } @article{SmirnovKronbergLatallerieetal.2019, author = {Smirnov, Artem G. and Kronberg, Elena A. and Latallerie, F. and Daly, Patrick W. and Aseev, Nikita and Shprits, Yuri Y. and Kellerman, Adam C. and Kasahara, Satoshi and Turner, Drew L. and Taylor, M. G. G. T.}, title = {Electron Intensity Measurements by the Cluster/RAPID/IES Instrument in Earth's Radiation Belts and Ring Current}, series = {Space Weather: The International Journal of Research and Applications}, volume = {17}, journal = {Space Weather: The International Journal of Research and Applications}, number = {4}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1542-7390}, doi = {10.1029/2018SW001989}, pages = {553 -- 566}, year = {2019}, abstract = {Plain Language Summary Radiation belts of the Earth, which are the zones of charged energetic particles trapped by the geomagnetic field, comprise enormous and dynamic systems. While the inner radiation belt, composed mainly of high-energy protons, is relatively stable, the outer belt, filled with energetic electrons, is highly variable and depends substantially on solar activity. Hence, extended reliable observations and the improved models of the electron intensities in the outer belt depending on solar wind parameters are necessary for prediction of their dynamics. The Cluster mission has been measuring electron flux intensities in the radiation belts since its launch in 2000, thus providing a huge dataset that can be used for radiation belts analysis. Using 16 years of electron measurements by the Cluster mission corrected for background contamination, we derived a uniform linear-logarithmic dependence of electron fluxes in the outer belt on the solar wind dynamic pressure.}, language = {en} } @misc{Shprits2017, author = {Shprits, Yuri Y.}, title = {Editorial: Topical Collection on the Lomonosov Mission}, series = {Space science reviews}, volume = {212}, journal = {Space science reviews}, publisher = {Springer}, address = {Dordrecht}, issn = {0038-6308}, doi = {10.1007/s11214-017-0393-1}, pages = {1685 -- 1686}, year = {2017}, language = {en} } @article{LaiMartinezSanchezCahoyetal.2017, author = {Lai, Shu T. and Martinez-Sanchez, Manuel and Cahoy, Kerri and Thomsen, Michelle F. and Shprits, Yuri Y. and Lohmeyer, Whitney and Wong, Frankie K.}, title = {Does spacecraft potential depend on the ambient electron density?}, series = {IEEE Transactions on Plasma Science}, volume = {45}, journal = {IEEE Transactions on Plasma Science}, publisher = {IEEE}, address = {New York}, issn = {0093-3813}, doi = {10.1109/TPS.2017.2751002}, pages = {2875 -- 2884}, year = {2017}, abstract = {In this paper, we address the question of whether spacecraft potential depends on the ambient electron density. In Maxwellian space plasmas, the onset of spacecraft charging does not depend on the ambient electron density. The balance of electron currents causes the incoming electrons to balance with the outgoing secondary electrons. The onset is controlled by the critical or anticritical temperature of the ambient electrons, but not the electron density. Above the critical temperature, charging to negative potential occurs. If the energy of the incoming electrons increases to well beyond the second crossing point of the secondary electron yield (SEY), the value of SEY decreases to well below unity. When the secondary electron current is negligible compared with the primary electron current, the spacecraft potential is governed solely by the balance of the incoming electrons and the sum of the currents of the repelled electrons and the attracted ions. In neutral space plasma, the electron and ion charges cancel each other. But if the space plasma deviates from being neutral, then the densities can have effect on the spacecraft potential. If the ambient plasma deviates significantly from equilibrium, a non-Maxwellian electron distribution may result. For a kappa distribution, one can show that the spacecraft charging level is independent of the ambient electron density. For a double Maxwellian distribution, the spacecraft charging level depends on the electron densities. For a conducting spacecraft charging in sunlight, the charging level is low and positive. It also depends on the ambient electron density. For a dielectric spacecraft in sunlight, the high-level negative-voltage charging on the shadowed side may extend to the sunlit side and block the photoelectrons trying to escape from the sunlit side. In this case, the charging level does not depend on ambient electron density. Using coordinated environmental and spacecraft charging data obtained from the Los Alamos National Laboratory geosynchronous satellites, we showed some results confirming that spacecraft potential is indeed often independent of the ambient electron density.}, language = {en} } @misc{ShpritsZhelavskayaGreenetal.2018, author = {Shprits, Yuri Y. and Zhelavskaya, Irina and Green, Janet C. and Pulkkinen, Antti A. and Horne, Richard B. and Pitchford, David and Glover, Alexi}, title = {Discussions on Stakeholder Requirements for Space Weather-Related Models}, series = {Space Weather: The International Journal of Research and Applications}, volume = {16}, journal = {Space Weather: The International Journal of Research and Applications}, number = {4}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1542-7390}, doi = {10.1002/2018SW001864}, pages = {341 -- 342}, year = {2018}, abstract = {Participants of the 2017 European Space Weather Week in Ostend, Belgium, discussed the stakeholder requirements for space weather-related models. It was emphasized that stakeholders show an increased interest in space weather-related models. Participants of the meeting discussed particular prediction indicators that can provide first-order estimates of the impact of space weather on engineering systems.}, language = {en} } @article{DrozdovAllisonShpritsetal.2022, author = {Drozdov, Alexander and Allison, Hayley J. and Shprits, Yuri Y. and Usanova, Maria E. and Saikin, Anthony and Wang, Dedong}, title = {Depletions of Multi-MeV Electrons and their association to Minima in Phase Space Density}, series = {Geophysical research letters}, volume = {49}, journal = {Geophysical research letters}, number = {8}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1029/2021GL097620}, pages = {11}, year = {2022}, abstract = {Fast-localized electron loss, resulting from interactions with electromagnetic ion cyclotron (EMIC) waves, can produce deepening minima in phase space density (PSD) radial profiles. Here, we perform a statistical analysis of local PSD minima to quantify how readily these are associated with radiation belt depletions. The statistics of PSD minima observed over a year are compared to the Versatile Electron Radiation Belts (VERB) simulations, both including and excluding EMIC waves. The observed minima distribution can only be achieved in the simulation including EMIC waves, indicating their importance in the dynamics of the radiation belts. By analyzing electron flux depletions in conjunction with the observed PSD minima, we show that, in the heart of the outer radiation belt (L* < 5), on average, 53\% of multi-MeV electron depletions are associated with PSD minima, demonstrating that fast localized loss by interactions with EMIC waves are a common and crucial process for ultra-relativistic electron populations.}, language = {en} } @article{KimShprits2017, author = {Kim, Kyung-Chan and Shprits, Yuri Y.}, title = {Dependence of the amplitude of magnetosonic waves on the solar wind and AE index using Van Allen Probes}, series = {Journal of geophysical research : Space physics}, volume = {122}, journal = {Journal of geophysical research : Space physics}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1002/2017JA024094}, pages = {6022 -- 6034}, year = {2017}, abstract = {We present the dependence of the magnetosonic wave amplitudes both outside and inside the plasmapause on the solar wind and AE index using Van Allen Probe-A spacecraft during the time period of 1 October 2012 to 31 December 2015, based on a correlation and regression analysis. Solar wind parameters considered are the southward interplanetary magnetic field (IMF B-S), solar wind number density (N-SW), and bulk speed (V-SW). We find that the wave amplitudes outside (inside) the plasmapause are well correlated with the preceding AE, IMF B-S, and N-SW with time delays, each corresponding to 2-3 h (3-4 h), 4-5 h (3-4 h), and 2-3 h (8-9 h), while the correlation with V-SW is ambiguous both inside and outside the plasmapause. As measured by the correlation coefficient, the IMF B-S is the most influential solar wind parameter that affects the dayside wave amplitudes both outside and inside the plasmapause, while N-SW contributes to enhancing the duskside waves outside the plasmapause. The AE effect on wave amplitudes is comparable to that of IMF B-S. More interestingly, regression with time histories of the solar wind parameters and the AE index preceding the wave measurements outside the plasmapause shows significant dependence on the IMF B-S, N-SW, and AE: the region of peak coefficients is changed with time delay for IMF B-S and AE, while isolated peaks around duskside remain gradually decrease with time for N-SW. In addition, the regression with magnetosonic waves inside the plasmapause shows high coefficients around prenoon sector with preceding IMF B-S and V-SW.}, language = {en} } @article{KronbergRashevDalyetal.2016, author = {Kronberg, Elena A. and Rashev, M. V. and Daly, P. W. and Shprits, Yuri Y. and Turner, D. L. and Drozdov, Alexander and Dobynde, M. and Kellerman, Adam C. and Fritz, T. A. and Pierrard, V. and Borremans, K. and Klecker, B. and Friedel, R.}, title = {Contamination in electron observations of the silicon detector on board}, series = {Space Weather: The International Journal of Research and Applications}, volume = {14}, journal = {Space Weather: The International Journal of Research and Applications}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1542-7390}, doi = {10.1002/2016SW001369}, pages = {449 -- 462}, year = {2016}, abstract = {Since more than 15 years, the Cluster mission passes through Earth's radiation belts at least once every 2 days for several hours, measuring the electron intensity at energies from 30 to 400 keV. These data have previously been considered not usable due to contamination caused by penetrating energetic particles (protons at >100 keV and electrons at >400 keV). In this study, we assess the level of distortion of energetic electron spectra from the Research with Adaptive Particle Imaging Detector (RAPID)/Imaging Electron Spectrometer (IES) detector, determining the efficiency of its shielding. We base our assessment on the analysis of experimental data and a radiation transport code (Geant4). In simulations, we use the incident particle energy distribution of the AE9/AP9 radiation belt models. We identify the Roederer L values, L\&\#8902;, and energy channels that should be used with caution: at 3\&\#8804;L\&\#8902;\&\#8804;4, all energy channels (40-400 keV) are contaminated by protons (\&\#8771;230 to 630 keV and >600 MeV); at L\&\#8902;\&\#8771;1 and 4-6, the energy channels at 95-400 keV are contaminated by high-energy electrons (>400 keV). Comparison of the data with electron and proton observations from RBSP/MagEIS indicates that the subtraction of proton fluxes at energies \&\#8771; 230-630 keV from the IES electron data adequately removes the proton contamination. We demonstrate the usefulness of the corrected data for scientific applications.}, language = {en} } @article{SaikinJordanovaZhangetal.2018, author = {Saikin, Anthony and Jordanova, Vania K. and Zhang, J. C. and Smith, C. W. and Spence, H. E. and Larsen, B. A. and Reeves, G. D. and Torbert, R. B. and Kletzing, C. A. and Zhelayskaya, I. S. and Shprits, Yuri Y.}, title = {Comparing simulated and observed EMIC wave amplitudes using in situ Van}, series = {Journal of Atmospheric and Solar-Terrestrial Physics}, volume = {177}, journal = {Journal of Atmospheric and Solar-Terrestrial Physics}, publisher = {Elsevier}, address = {Oxford}, issn = {1364-6826}, doi = {10.1016/j.jastp.2018.01.024}, pages = {190 -- 201}, year = {2018}, abstract = {We perform a statistical study calculating electromagnetic ion cyclotron (EMIC) wave amplitudes based off in situ plasma measurements taken by the Van Allen Probes' (1.1-5.8 Re) Helium, Oxygen, Proton, Electron (HOPE) instrument. Calculated wave amplitudes are compared to EMIC waves observed by the Electric and Magnetic Field Instrument Suite and Integrated Science on board the Van Allen Probes during the same period. The survey covers a 22-month period (1 November 2012 to 31 August 2014), a full Van Allen Probe magnetic local time (MLT) precession. The linear theory proxy was used to identify EMIC wave events with plasma conditions favorable for EMIC wave excitation. Two hundred and thirty-two EMIC wave events (103 H+-band and 129 He+-band) were selected for this comparison. Nearly all events selected are observed beyond L = 4. Results show that calculated wave amplitudes exclusively using the in situ HOPE measurements produce amplitudes too low compared to the observed EMIC wave amplitudes. Hot proton anisotropy (Ahp) distributions are asymmetric in MLT within the inner (L < 7) magnetosphere with peak (minimum) Ahp, ∼0.81 to 1.00 (∼0.62), observed in the dawn (dusk), 0000 < MLT ≤ 1200 (1200 < MLT ≤ 2400), sectors. Measurements of Ahp are found to decrease in the presence of EMIC wave activity. Ahp amplification factors are determined and vary with respect to EMIC wave-band and MLT. He+-band events generally require double (quadruple) the measured Ahp for the dawn (dusk) sector to reproduce the observed EMIC wave amplitudes.}, language = {en} } @article{ProlSmirnovHoqueetal.2022, author = {Prol, Fabricio S. and Smirnov, Artem G. and Hoque, M. Mainul and Shprits, Yuri Y.}, title = {Combined model of topside ionosphere and plasmasphere derived from radio-occultation and Van Allen Probes data}, series = {Scientific reports}, volume = {12}, journal = {Scientific reports}, number = {1}, publisher = {Macmillan Publishers Limited, part of Springer Nature}, address = {London}, issn = {2045-2322}, doi = {10.1038/s41598-022-13302-1}, pages = {11}, year = {2022}, abstract = {In the last years, electron density profile functions characterized by a linear dependence on the scale height showed good results when approximating the topside ionosphere. The performance above 800 km, however, is not yet well investigated. This study investigates the capability of the semi-Epstein functions to represent electron density profiles from the peak height up to 20,000 km. Electron density observations recorded by the Van Allen Probes were used to resolve the scale height dependence in the plasmasphere. It was found that the linear dependence of the scale height in the topside ionosphere cannot be directly used to extrapolate profiles above 800 km. We find that the dependence of scale heights on altitude is quadratic in the plasmasphere. A statistical model of the scale heights is therefore proposed. After combining the topside ionosphere and plasmasphere by a unified model, we have obtained good estimations not only in the profile shapes, but also in the Total Electron Content magnitude and distributions when compared to actual measurements from 2013, 2014, 2016 and 2017. Our investigation shows that Van Allen Probes can be merged to radio-occultation data to properly represent the upper ionosphere and plasmasphere by means of a semi-Epstein function.}, language = {en} } @article{ZhelavskayaSpasojevicShpritsetal.2016, author = {Zhelavskaya, Irina and Spasojevic, M. and Shprits, Yuri Y. and Kurth, William S.}, title = {Automated determination of electron density from electric field measurements on the Van Allen Probes spacecraft}, series = {Journal of geophysical research : Space physics}, volume = {121}, journal = {Journal of geophysical research : Space physics}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1002/2015JA022132}, pages = {4611 -- 4625}, year = {2016}, abstract = {We present the Neural-network-based Upper hybrid Resonance Determination (NURD) algorithm for automatic inference of the electron number density from plasma wave measurements made on board NASA's Van Allen Probes mission. A feedforward neural network is developed to determine the upper hybrid resonance frequency, fuhr, from electric field measurements, which is then used to calculate the electron number density. In previous missions, the plasma resonance bands were manually identified, and there have been few attempts to do robust, routine automated detections. We describe the design and implementation of the algorithm and perform an initial analysis of the resulting electron number density distribution obtained by applying NURD to 2.5 years of data collected with the Electric and Magnetic Field Instrument Suite and Integrated Science (EMFISIS) instrumentation suite of the Van Allen Probes mission. Densities obtained by NURD are compared to those obtained by another recently developed automated technique and also to an existing empirical plasmasphere and trough density model.}, language = {en} } @article{WangShpritsZhelayskayaetal.2019, author = {Wang, Dedong and Shprits, Yuri Y. and Zhelayskaya, Irina S. and Agapitov, Oleksiy and Drozdov, Alexander and Aseev, Nikita}, title = {Analytical chorus wave model derived from van Allen Probe Observations}, series = {Journal of geophysical research : Space physics}, volume = {124}, journal = {Journal of geophysical research : Space physics}, number = {2}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2018JA026183}, pages = {1063 -- 1084}, year = {2019}, abstract = {Chorus waves play an important role in the dynamic evolution of energetic electrons in the Earth's radiation belts and ring current. Using more than 5 years of Van Allen Probe data, we developed a new analytical model for upper-band chorus (UBC; 0.5fce < f < fce) and lower-band chorus (LBC; 0.05fce < f < 0.5fce) waves, where fce is the equatorial electron gyrofrequency. By applying polynomial fits to chorus wave root mean square amplitudes, we developed regression models for LBC and UBC as a function of geomagnetic activity (Kp), L, magnetic latitude (λ), and magnetic local time (MLT). Dependence on Kp is separated from the dependence on λ, L, and MLT as Kp-scaling law to simplify the calculation of diffusion coefficients and inclusion into particle tracing codes. Frequency models for UBC and LBC are also developed, which depends on MLT and magnetic latitude. This empirical model is valid in all MLTs, magnetic latitude up to 20°, Kp ≤ 6, L-shell range from 3.5 to 6 for LBC and from 4 to 6 for UBC. The dependence of root mean square amplitudes on L are different for different bands, which implies different energy sources for different wave bands. This analytical chorus wave model is convenient for inclusion in quasi-linear diffusion calculations of electron scattering rates and particle simulations in the inner magnetosphere, especially for the newly developed four-dimensional codes, which require significantly improved wave parameterizations.}, language = {en} } @article{ZhuShpritsChenetal.2018, author = {Zhu, Hui and Shprits, Yuri Y. and Chen, Lunjin and Liu, Xu and Kellerman, Adam C.}, title = {An event on simultaneous amplification of exohiss and chorus waves associated with electron density enhancements}, series = {Journal of geophysical research : Space physics}, volume = {123}, journal = {Journal of geophysical research : Space physics}, number = {11}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2017JA025023}, pages = {8958 -- 8968}, year = {2018}, abstract = {Whistler mode exohiss are the structureless hiss waves observed outside the plasma pause with featured equatorward Poynting flux. An event of the amplification of exohiss as well as chorus waves was recorded by Van Allen Probes during the recovery phase of a weak geomagnetic storm. Amplitudes of both types of the waves showed a significant increase at the regions of electron density enhancements. It is found that the electrons resonant with exohiss and chorus showed moderate pitch angle anisotropies. The ratio of the number of electrons resonating with exohiss to total electron number presented in-phase correlation with density variations, which suggests that exohiss can be amplified due to electron density enhancement in terms of cyclotron instability. The calculation of linear growth rates further supports above conclusion. We suggest that exohiss waves have potential to become more significant due to the background plasma fluctuation.}, language = {en} } @article{SmirnovShpritsAllisonetal.2022, author = {Smirnov, Artem and Shprits, Yuri Y. and Allison, Hayley and Aseev, Nikita and Drozdov, Alexander and Kollmann, Peter and Wang, Dedong and Saikin, Anthony}, title = {An empirical model of the equatorial electron pitch angle distributions in earth's outer radiation belt}, series = {Space Weather: the International Journal of Research and Applications}, volume = {20}, journal = {Space Weather: the International Journal of Research and Applications}, number = {9}, publisher = {American Geophysical Union}, address = {Washington, DC}, issn = {1542-7390}, doi = {10.1029/2022SW003053}, pages = {17}, year = {2022}, abstract = {In this study, we present an empirical model of the equatorial electron pitch angle distributions (PADs) in the outer radiation belt based on the full data set collected by the Magnetic Electron Ion Spectrometer (MagEIS) instrument onboard the Van Allen Probes in 2012-2019. The PADs are fitted with a combination of the first, third and fifth sine harmonics. The resulting equation resolves all PAD types found in the outer radiation belt (pancake, flat-top, butterfly and cap PADs) and can be analytically integrated to derive omnidirectional flux. We introduce a two-step modeling procedure that for the first time ensures a continuous dependence on L, magnetic local time and activity, parametrized by the solar wind dynamic pressure. We propose two methods to reconstruct equatorial electron flux using the model. The first approach requires two uni-directional flux observations and is applicable to low-PA data. The second method can be used to reconstruct the full equatorial PADs from a single uni- or omnidirectional measurement at off-equatorial latitudes. The model can be used for converting the long-term data sets of electron fluxes to phase space density in terms of adiabatic invariants, for physics-based modeling in the form of boundary conditions, and for data assimilation purposes.}, language = {en} } @misc{SmirnovKronbergDalyetal.2020, author = {Smirnov, Artem G. and Kronberg, Elena A. and Daly, Patrick W. and Aseev, Nikita and Shprits, Yuri Y. and Kellerman, Adam C.}, title = {Adiabatic Invariants Calculations for Cluster Mission: A Long-Term Product for Radiation Belts Studies}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {2}, issn = {1866-8372}, doi = {10.25932/publishup-52391}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-523915}, pages = {14}, year = {2020}, abstract = {The Cluster mission has produced a large data set of electron flux measurements in the Earth's magnetosphere since its launch in late 2000. Electron fluxes are measured using Research with Adaptive Particle Imaging Detector (RAPID)/Imaging Electron Spectrometer (IES) detector as a function of energy, pitch angle, spacecraft position, and time. However, no adiabatic invariants have been calculated for Cluster so far. In this paper we present a step-by-step guide to calculations of adiabatic invariants and conversion of the electron flux to phase space density (PSD) in these coordinates. The electron flux is measured in two RAPID/IES energy channels providing pitch angle distribution at energies 39.2-50.5 and 68.1-94.5 keV in nominal mode since 2004. A fitting method allows to expand the conversion of the differential fluxes to the range from 40 to 150 keV. Best data coverage for phase space density in adiabatic invariant coordinates can be obtained for values of second adiabatic invariant, K, similar to 10(2), and values of the first adiabatic invariant mu in the range approximate to 5-20 MeV/G. Furthermore, we describe the production of a new data product "LSTAR," equivalent to the third adiabatic invariant, available through the Cluster Science Archive for years 2001-2018 with 1-min resolution. The produced data set adds to the availability of observations in Earth's radiation belts region and can be used for long-term statistical purposes.}, language = {en} } @article{SmirnovKronbergDalyetal.2020, author = {Smirnov, Artem G. and Kronberg, Elena A. and Daly, Patrick W. and Aseev, Nikita and Shprits, Yuri Y. and Kellerman, Adam C.}, title = {Adiabatic Invariants Calculations for Cluster Mission: A Long-Term Product for Radiation Belts Studies}, series = {Journal of Geophysical Research: Space Physics}, volume = {125}, journal = {Journal of Geophysical Research: Space Physics}, number = {2}, publisher = {John Wiley \& Sons, Inc.}, address = {New Jersey}, pages = {12}, year = {2020}, abstract = {The Cluster mission has produced a large data set of electron flux measurements in the Earth's magnetosphere since its launch in late 2000. Electron fluxes are measured using Research with Adaptive Particle Imaging Detector (RAPID)/Imaging Electron Spectrometer (IES) detector as a function of energy, pitch angle, spacecraft position, and time. However, no adiabatic invariants have been calculated for Cluster so far. In this paper we present a step-by-step guide to calculations of adiabatic invariants and conversion of the electron flux to phase space density (PSD) in these coordinates. The electron flux is measured in two RAPID/IES energy channels providing pitch angle distribution at energies 39.2-50.5 and 68.1-94.5 keV in nominal mode since 2004. A fitting method allows to expand the conversion of the differential fluxes to the range from 40 to 150 keV. Best data coverage for phase space density in adiabatic invariant coordinates can be obtained for values of second adiabatic invariant, K, similar to 10(2), and values of the first adiabatic invariant mu in the range approximate to 5-20 MeV/G. Furthermore, we describe the production of a new data product "LSTAR," equivalent to the third adiabatic invariant, available through the Cluster Science Archive for years 2001-2018 with 1-min resolution. The produced data set adds to the availability of observations in Earth's radiation belts region and can be used for long-term statistical purposes.}, language = {en} } @article{ShpritsAllisonWangetal.2022, author = {Shprits, Yuri Y. and Allison, Hayley J. and Wang, Dedong and Drozdov, Alexander and Szabo-Roberts, Matyas and Zhelavskaya, Irina and Vasile, Ruggero}, title = {A new population of ultra-relativistic electrons in the outer radiation zone}, series = {Journal of geophysical research : Space physics}, volume = {127}, journal = {Journal of geophysical research : Space physics}, number = {5}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2021JA030214}, pages = {34}, year = {2022}, abstract = {Van Allen Probes measurements revealed the presence of the most unusual structures in the ultra-relativistic radiation belts. Detailed modeling, analysis of pitch angle distributions, analysis of the difference between relativistic and ultra-realistic electron evolution, along with theoretical studies of the scattering and wave growth, all indicate that electromagnetic ion cyclotron (EMIC) waves can produce a very efficient loss of the ultra-relativistic electrons in the heart of the radiation belts. Moreover, a detailed analysis of the profiles of phase space densities provides direct evidence for localized loss by EMIC waves. The evolution of multi-MeV fluxes shows dramatic and very sudden enhancements of electrons for selected storms. Analysis of phase space density profiles reveals that growing peaks at different values of the first invariant are formed at approximately the same radial distance from the Earth and show the sequential formation of the peaks from lower to higher energies, indicating that local energy diffusion is the dominant source of the acceleration from MeV to multi-MeV energies. Further simultaneous analysis of the background density and ultra-relativistic electron fluxes shows that the acceleration to multi-MeV energies only occurs when plasma density is significantly depleted outside of the plasmasphere, which is consistent with the modeling of acceleration due to chorus waves.}, language = {en} } @misc{SadovnichiiPanasyukAmelyushkinetal.2017, author = {Sadovnichii, V. A. and Panasyuk, M. I. and Amelyushkin, A. M. and Benghin, V. V. and Garipov, G. K. and Kalegaev, V. V. and Klimov, P. A. and Khrenov, B. A. and Petrov, V. L. and Sharakin, S. A. and Shirokov, A. V. and Svertilov, S. I. and Zotov, M. Y. and Yashin, I. V. and Gorbovskoy, E. S. and Lipunov, V. M. and Park, I. H. and Lee, J. and Jeong, S. and Kim, M. B. and Jeong, H. M. and Shprits, Yuri Y. and Angelopoulos, V. and Russell, C. T. and Runov, A. and Turner, D. and Strangeway, R. J. and Caron, R. and Biktemerova, S. and Grinyuk, A. and Lavrova, M. and Tkachev, L. and Tkachenko, A. and Martinez, O. and Salazar, H. and Ponce, E.}, title = {"Lomonosov" satellite-space observatory to study extreme phenomena in space}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {959}, issn = {1866-8372}, doi = {10.25932/publishup-42818}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-428185}, pages = {1705 -- 1738}, year = {2017}, abstract = {The "Lomonosov" space project is lead by Lomonosov Moscow State University in collaboration with the following key partners: Joint Institute for Nuclear Research, Russia, University of California, Los Angeles (USA), University of Pueblo (Mexico), Sungkyunkwan University (Republic of Korea) and with Russian space industry organi-zations to study some of extreme phenomena in space related to astrophysics, astroparticle physics, space physics, and space biology. The primary goals of this experiment are to study: -Ultra-high energy cosmic rays (UHECR) in the energy range of the Greizen-ZatsepinKuzmin (GZK) cutoff; -Ultraviolet (UV) transient luminous events in the upper atmosphere; -Multi-wavelength study of gamma-ray bursts in visible, UV, gamma, and X-rays; -Energetic trapped and precipitated radiation (electrons and protons) at low-Earth orbit (LEO) in connection with global geomagnetic disturbances; -Multicomponent radiation doses along the orbit of spacecraft under different geomagnetic conditions and testing of space segments of optical observations of space-debris and other space objects; -Instrumental vestibular-sensor conflict of zero-gravity phenomena during space flight. This paper is directed towards the general description of both scientific goals of the project and scientific equipment on board the satellite. The following papers of this issue are devoted to detailed descriptions of scientific instruments.}, language = {en} } @misc{SadovnichiiPanasyukAmelyushkinetal.2017, author = {Sadovnichii, V. A. and Panasyuk, M. I. and Amelyushkin, A. M. and Bogomolov, V. V. and Benghin, V. V. and Garipov, G. K. and Kalegaev, V. V. and Klimov, P. A. and Khrenov, B. A. and Petrov, V. L. and Sharakin, S. A. and Shirokov, A. V. and Svertilov, S. I. and Zotov, M. Y. and Yashin, I. V. and Gorbovskoy, E. S. and Lipunov, V. M. and Park, I. H. and Lee, J. and Jeong, S. and Kim, M. B. and Jeong, H. M. and Shprits, Yuri Y. and Angelopoulos, V. and Russell, C. T. and Runov, A. and Turner, D. and Strangeway, R. J. and Caron, R. and Biktemerova, S. and Grinyuk, A. and Lavrova, M. and Tkachev, L. and Tkachenko, A. and Martinez, O. and Salazar, H. and Ponce, E.}, title = {"Lomonosov" Satellite-Space Observatory to Study Extreme Phenomena in Space}, series = {Space science reviews}, volume = {212}, journal = {Space science reviews}, publisher = {Springer}, address = {Dordrecht}, issn = {0038-6308}, doi = {10.1007/s11214-017-0425-x}, pages = {1705 -- 1738}, year = {2017}, abstract = {The "Lomonosov" space project is lead by Lomonosov Moscow State University in collaboration with the following key partners: Joint Institute for Nuclear Research, Russia, University of California, Los Angeles (USA), University of Pueblo (Mexico), Sungkyunkwan University (Republic of Korea) and with Russian space industry organi-zations to study some of extreme phenomena in space related to astrophysics, astroparticle physics, space physics, and space biology. The primary goals of this experiment are to study: -Ultra-high energy cosmic rays (UHECR) in the energy range of the Greizen-ZatsepinKuzmin (GZK) cutoff; -Ultraviolet (UV) transient luminous events in the upper atmosphere; -Multi-wavelength study of gamma-ray bursts in visible, UV, gamma, and X-rays; -Energetic trapped and precipitated radiation (electrons and protons) at low-Earth orbit (LEO) in connection with global geomagnetic disturbances; -Multicomponent radiation doses along the orbit of spacecraft under different geomagnetic conditions and testing of space segments of optical observations of space-debris and other space objects; -Instrumental vestibular-sensor conflict of zero-gravity phenomena during space flight. This paper is directed towards the general description of both scientific goals of the project and scientific equipment on board the satellite. The following papers of this issue are devoted to detailed descriptions of scientific instruments.}, language = {en} }