ALBERT

Description: Knowledge of the plasma fluxes at geosynchronous orbit is important to both scientific and operational investigations. We present a new empirical model of the ion flux and the electron flux at geosynchronous orbit (GEO) in the energy range ~1 eV to ~40 keV. The model is based on a total of 82 satellite-years of observations from the Magnetospheric Plasma Analyzer instruments on Los Alamos National Laboratory satellites at GEO. These data are assigned to a fixed grid of 24 local-times and 40 energies, at all possible values of Kp. Bi-linear interpolation is used between grid points to provide the ion flux and the electron flux values at any energy and local-time, and for given values of geomagnetic activity (proxied by the 3-hour Kp index), and also for given values of solar activity (proxied by the daily F10.7 index). Initial comparison of the electron flux from the model with data from a Compact Environmental Anomaly Sensor II (CEASE-II), also located at geosynchronous orbit, indicate a good match during both quiet and disturbed periods. The model is available for distribution as a FORTRAN code that can be modified to suit user-requirements.

Description: A number of relativistic electron loss processes exist in the inner magnetosphere, and the extent to which MeV electron precipitation into Earth's atmosphere plays a role in radiation belt dynamics is a topic of much debate. In this work, we investigate the contribution of electron precipitation to radiation belt losses, looking at what times and locations precipitation is important. Through high-cadence low-altitude measurements from the SAMPEX satellite, we examine the distributions of millisecond (microburst) as well as longer duration (band-type) precipitation and the relative contributions of these two precipitation types to radiation belt dynamics during high speed stream (HSS) driven storms. Different local time and radial distributions between microbursts and precipitation bands suggest different scattering mechanisms as the causes of the two precipitation types. In a superposed epoch study of 42 HSS-driven storms, enhanced main and recovery phase losses to the atmosphere are observed. Microburst occurrence rates peak in the recovery phase of the storms, while their magnitudes remain fairly constant over storm phase. Precipitation bands show an increase in both occurrence and magnitude at storm onset, particularly at the inner edge of the outer radiation belt. The observations, enabled by the high time resolution and large geometric factor and field of view of the SAMPEX/HILT instrument, reveal when and where microburst and band-type precipitation are contributing to radiation belt dynamics during HSS-driven storms.

Description: The distribution of mass density along the field lines affects the ratios of toroidal (azimuthally oscillating) Alfvén frequencies, and given the ratios of these frequencies we can get information about that distribution. Here we assume the commonlyused power law form for the field line distribution, ρ m  =  ρ m , eq ( LR E / R ) α , where ρ m , eq is the value of the mass density ρ m at the magnetic equator, L is the L shell, R E is the Earth's radius, R is the geocentric distance to a point on the field line, and α is the power law coefficient. Positive values of α indicate that ρ m increases away fromthe magnetic equator, zero value indicates that ρ m is constant along the magnetic field line, and negative α indicates that there is a local peak in ρ m at the magnetic equator. Using 12 years of observations of toroidal Alfvén frequencies by the Geostationary Operational Environmental Satellites (GOES), we study the typical dependence of inferred values of α on the magnetic local time (MLT), the phase of the solar cycle as specified bythe F10.7 extreme ultraviolet solar flux, and geomagnetic activity as specified by the auroral electrojet (AE) index. Over the mostly dayside range of the observations, we find that α decreases with respect to increasing MLT and F10.7, but increases with respect to increasing AE. We develop a formula that depends on all three parameters, α 3Dmodel  = 2.2 + 1.3 ⋅  cos (MLT ⋅ 15 ∘ ) + 0.0026 ⋅ AE ⋅  cos ((MLT − 0.8) ⋅ 15 ∘ ) + 2.1 ⋅ 10 − 5  ⋅ AE ⋅ F10.7 − 0.010 ⋅ F10.7, that models the binned values of α within a standard deviation of 0.3. While we do not yet have a complete theoretical understanding of why α should depend on these parameters in such a way, we do make some observations and speculations about the causes. At least part of the dependence is related to that of ρ m , eq ; higher α , corresponding to steeper variation with respect to MLAT, occurs when ρ m , eq is lower.

Description: We investigate structural and thermodynamic properties of aqueous dispersions of ionic microgels—soft colloidal gel particles that exhibit unusual phase behavior. Starting from a coarse-grained model of microgel macroions as charged spheres that are permeable to microions, we perform simulations and theoretical calculations using two complementary implementations of Poisson-Boltzmann (PB) theory. Within a one-component model, based on a linear-screening approximation for effective electrostatic pair interactions, we perform molecular dynamics simulations to compute macroion-macroion radial distribution functions, static structure factors, and macroion contributions to the osmotic pressure. For the same model, using a variational approximation for the free energy, we compute both macroion and microion contributions to the osmotic pressure. Within a spherical cell model, which neglects macroion correlations, we solve the nonlinear PB equation to compute microion distributions and osmotic pressures. By comparing the one-component and cell model implementations of PB theory, we demonstrate that the linear-screening approximation is valid for moderately charged microgels. By further comparing cell model predictions with simulation data for osmotic pressure, we chart the cell model’s limits in predicting osmotic pressures of salty dispersions.

Description: We simulate whistler mode waves using a hybrid code. There are four species in the simulations, hot electrons initialized with a bi-Maxwellian distribution with temperature in the direction perpendicular to background magnetic field greater than that in the parallel direction, warm isotropic electrons, cold inertialess fluid electrons and protons as an immobile background. The density of the hot population is a small fraction of the total plasma density. Comparison between the dispersion relation of our model and other dispersion relations shows that our model is more accurate for lower frequency whistlers than for higher frequency whistlers. Simulations in 2-D Cartesian coordinates agree very well with those using a full dynamics code. In the 1-D simulations along the dipole magnetic field, the predicted frequency and wave number are observed. Rising tones are observed in the 1/14 scale simulations that have larger than realistic magnetic field spatial inhomogeneity. However, in the full scale 1-D simulation in a dipole field, the waves are more broadband and do not exhibit rising tones. In the 2-D simulations in a meridional plane, the waves are generated with propagation approximately parallel to the background magnetic field. However, the wave fronts become oblique as they propagate to higher latitudes. Simulations with different plasma density profiles across L -shell are performed to study the effect of the background density on whistler propagation.

Description: Relativistic electrons have been thought to more easily resonate with electromagnetic ion cyclotron (EMIC) waves if the total density is large. We show that, for a particular EMIC mode, this dependence is weak due to the dependence of the wave frequency and wave vector on the density. A significant increase in relativistic electron minimum resonant energy might occur for the H band EMIC mode only for small density, but no changes in parameters significantly decrease the minimum resonant energy from a nominal value. The minimum resonant energy depends most strongly on the thermal velocity associated with the field line motion of the hot ring current protons that drive the instability. High density due to a plasmasphere or plasmaspheric plume could possibly lead to lower minimum resonance energy by causing the He band EMIC mode to be dominant. We demonstrate these points using parameters from a ring current simulation.

Description: We present observations of the radiation belts from the HOPE and MagEIS particle detectors on the Van Allen Probes satellites that illustrate the energy-dependence and L-shell dependence of radiation belt enhancements and decays. We survey events in 2013 and analyze an event on March 1 in more detail. The observations show: (a) At all L-shells, lower-energy electrons are enhanced more often than higher energies; (b) Events that fill the slot region are more common at lower energies; (c) Enhancements of electrons in the inner zone are more common at lower energies; and (d) Even when events do not fully fill the slot region, enhancements at lower-energies tend to extend to lower L-shells than higher energies. During enhancement events the outer zone extends to lower L-shells at lower energies while being confined to higher L-shells at higher energies. The inner zone shows the opposite with an outer boundary at higher L-shells for lower energies. Both boundaries are nearly straight in log(energy) vs. L-shell space. At energies below a few hundred keV radiation belt electron penetration through the slot region into the inner zone is commonplace but the number and frequency of “slot filling” events decreases with increasing energy. The inner zone is enhanced only at energies that penetrate through the slot. Energy- and L-shell dependent losses (that are consistent with whistler hiss interactions) return the belts to more quiescent conditions.

Description: Most theoretical wave models require the power in the wave magnetic field in order to determine the effect of chorus waves on radiation belt electrons. However, researchers typically use the cold plasma dispersion relation to approximate the magneticwave power when only electric field data are available. In this study, the validity of using the the cold plasma dispersion relation in this context is tested using EMFISIS observations of both the electric and magnetic spectral intensities in the chorus wave band (0.1-0.9 f ce ). Results from this study indicate that the calculated wave intensity is least accurate during periods of enhanced wave activity. For observed wave intensities 〉10¯ 3 nT 2 , using the cold plasma dispersion relation results in an underestimate of the wave intensity by a factor of 2 or greater, 56% of the time over the full chorus wave band, 60% of the time for lower band chorus, and 59% of the time for upper band chorus. Hence during active periods, empirical chorus wave models that are reliant on the cold plasma dispersion relation will underestimate chorus wave intensities to a significant degree, thus causing questionable calculation of wave-particle resonance effects on MeVelectrons.

Description: With data from Van Allen Probes, we investigate EMIC wave excitation using simultaneously observed ion distributions. Strong He-band waves occurred while the spacecraft was moving through an enhanced density region. We extract from Helium, Oxygen, Proton, and Electron (HOPE) Mass Spectrometer measurement the velocity distributions of warm heavy ions as well as anisotropic energetic protons that drive wave growth through the ion cyclotron instability. Fitting the measured ion fluxes to multiple sin m -type distribution functions, we find that the observed ions make up about 15% of the total ions, but about 85% of them are still missing. By making legitimate estimates of the unseen cold (below ~2 eV) ion composition from cutoff frequencies suggested by the observed wave spectrum, a series of linear instability analyses and hybrid simulations are carried out. The simulated waves generally vary as predicted by linear theory. They are more sensitive to the cold O+ concentration than the cold He+ concentration. Increasing the cold O+ concentration weakens the He-band waves but enhances the O-band waves. Finally, the exact cold ion composition is suggested to be in a range when the simulated wave spectrum best matches the observed one.