Science of the Van Allen Probes Science Operations Centers

The Van Allen Probes mission operations materialized through a distributed model in which operational responsibility was divided between the Mission Operations Center (MOC) and separate instrument specific SOCs. The sole MOC handled all aspects of telemetering and receiving tasks as well as certain scientifically relevant ancillary tasks. Each instrument science team developed individual instrument specific SOCs proficient in unique capabilities in support of science data acquisition, data processing, instrument performance, and tools for the instrument team scientists. In parallel activities, project scientists took on the task of providing a significant modeling tool base usable by the instrument science teams and the larger scientific community. With a mission as complex as Van Allen Probes, scientific inquiry occurred due to constant and significant collaboration between the SOCs and in concert with the project science team. Planned cross-instrument coordinated observations resulted in critical discoveries during the seven-year mission. Instrument cross-calibration activities elucidated a more seamless set of data products. Specific topics include post-launch changes and enhancements to the SOCs, discussion of coordination activities between the SOCs, SOC specific analysis software, modeling software provided by the Van Allen Probes project, and a section on lessons learned. One of the most significant lessons learned was the importance of the original decision to implement individual team SOCs providing timely and well-documented instrument data for the NASA Van Allen Probes Mission scientists and the larger magnetospheric and radiation belt scientific community.

GeV Proton Detection in the 8 November 2000 Solar Event

In this study, we analyze the L3 precision muon spectrometer data from November 2000. The results showed that a 4.7σmuon excess appeared at a time coincident with the solar flare of 8 November 2000. This muon excess corresponded to primary protons above 40 GeV, coming from a sky cell of solid angle 0.048 sr. The probability of being a background fluctuation was estimated to be about 0.1%. It is interesting and noteworthy that an M-class solar flare may also accelerate solar protons to such high energies.

Van Allen Probes Observations of Symmetric Stormtime Compressional ULF Waves

Previous spacecraft studies showed that stormtime poloidal ultralow-frequency (ULF) waves in the ring current region have an antisymmetric (second harmonic) mode structure about the magnetic equator. This paper reports Van Allen Probes observations of symmetric ULF waves in the postnoon sector during a moderate geomagnetic storm. The mode structure is determined from the presence of purely compressional magnetic field oscillations at the equator accompanied by strong transverse electric field perturbations. Antisymmetric waves were also detected but only very late in the recovery phase. The symmetric waves were detected outside the plasmasphere at L = 3.0-5.5 and had peak power at 4-10 mHz, lower than the frequency of the local fundamental toroidal standing Alfvén wave. During the wave events, the flux of protons was enhanced at energies below ∼5 keV, which appears to be a prerequisite for the waves. The protons may provide free energies to waves through drift resonance instability or drift compressional instability, which occur in the presence of radial gradients of plasma parameters.

Unsteady Magnetopause Reconnection Under Quasi-Steady Solar Wind Driving

The intrinsic temporal nature of magnetic reconnection at the magnetopause has been an active area of research. Both temporally steady and intermittent reconnection have been reported. We examine the steadiness of reconnection using space-ground conjunctions under quasi-steady solar wind driving. The spacecraft suggests that reconnection is first inactive, and then activates. The radar further suggests that after activation, reconnection proceeds continuously but unsteadily. The reconnection electric field shows variations at frequencies below 10 mHz with peaks at 3 and 5 mHz. The variation amplitudes are ∼10-30 mV/m in the ionosphere, and 0.3-0.8 mV/m at the equatorial magnetopause. Such amplitudes represent 30%-60% of the peak reconnection electric field. The unsteadiness of reconnection can be plausibly explained by the fluctuating magnetic field in the turbulent magnetosheath. A comparison with a previous global hybrid simulation suggests that it is the foreshock waves that drive the magnetosheath fluctuations, and hence modulate the reconnection.

The August 2018 Geomagnetic Storm Observed by the High-Energy Particle Detector on Board the CSES-01 Satellite

On 25 August 2018, a G3-class geomagnetic storm reached the Earth’s magnetosphere, causing a transient rearrangement of the charged particle environment around the planet, which was detected by the High-Energy Particle Detector (HEPD) on board the China Seismo-Electromagnetic Satellite (CSES-01). We found that the count rates of electrons in the MeV range were characterized by a depletion during the storm’s main phase and a clear enhancement during the recovery caused by large substorm activity, with the key role played by auroral processes mapped into the outer belt. A post-storm rate increase was localized at L-shells immediately above ∼3 and mostly driven by non-adiabatic local acceleration caused by possible resonant interaction with low-frequency magnetospheric waves.

Multi-Instrument Characterization of Magnetospheric Cold Plasma Dynamics in the June 22, 2015 Geomagnetic Storm

We present a comparison of magnetospheric plasma mass/electron density observations during an 11-day interval which includes the geomagnetic storm of June 22, 2015. For this study we used: Equatorial plasma mass density derived from geomagnetic field line resonances (FLRs) detected by Van Allen Probes and at the ground-based magnetometer networks EMMA and CARISMA; in situ electron density inferred by the Neural-network-based Upper hybrid Resonance Determination algorithm applied to plasma wave Van Allen Probes measurements. The combined observations at L ∼ 4, MLT ∼ 16 of the two longitudinally separated magnetometer networks show a temporal pattern very similar to that of the in situ observations: A density decrease by an order of magnitude about 1 day after the Dst minimum, a partial recovery a few hours later, and a new strong decrease soon after. The observations are consistent with the position of the measurement points with respect to the plasmasphere boundary as derived by a plasmapause test particle simulation. A comparison between plasma mass densities derived from ground and in situ FLR observations during favorable conjunctions shows a good agreement. We find however, for L < ∼3, the spacecraft measurements to be higher than the corresponding ground observations with increasing deviation with decreasing L, which might be related to the rapid outbound spacecraft motion in that region. A statistical analysis of the average ion mass using simultaneous spacecraft measurements of mass and electron density indicates values close to 1 amu in plasmasphere and higher values (∼2-3 amu) in plasmatrough.

Determining EMIC Wave Vector Properties Through Multi-Point Measurements: The Wave Curl Analysis

Electromagnetic ion cyclotron (EMIC) waves play important roles in particle loss processes in the magnetosphere. Determining the evolution of EMIC waves as they propagate and how this evolution affects wave-particle interactions requires accurate knowledge of the wave vector, k. We present a technique using the curl of the wave magnetic field to determine k observationally, enabled by the unique configuration and instrumentation of the Magnetospheric MultiScale (MMS) spacecraft. The wave curl analysis is demonstrated for synthetic arbitrary electromagnetic waves with varying properties typical of observed EMIC waves. The method is also applied to an EMIC wave interval observed by MMS on October 28, 2015. The derived wave properties and k from the wave curl analysis for the observed EMIC wave are compared with the Waves in Homogenous, Anisotropic, Multi-component Plasma (WHAMP) wave dispersion solution and with results from other single- and multi-spacecraft techniques. We find good agreement between k from the wave curl analysis, k determined from other observational techniques, and k determined from WHAMP. Additionally, the variation of k due to the time and frequency intervals used in the wave curl analysis is explored. This exploration demonstrates that the method is robust when applied to a wave containing at least 3-4 wave periods and over a rather wide frequency range encompassing the peak wave emission. These results provide confidence that we are able to directly determine the wave vector properties using this multi-spacecraft method implementation, enabling systematic studies of EMIC wave k properties with MMS.

Relative contributions of large-scale and wedgelet currents in the substorm current wedge

We examined how much large-scale and localized upward and downward currents contribute to the substorm current wedge (SCW), and how they evolve over time, using the THEMIS all-sky imagers (ASIs) and ground magnetometers. One type of events is dominated by a single large-scale wedge, with upward currents over the surge and broad downward currents poleward-eastward of the surge. The other type of events is a composite of large-scale wedge and wedgelets associated with streamers, with each wedgelet having comparable intensity to the large-scale wedge currents. Among 17 auroral substorms with wide ASI coverage, the composite current type is more frequent than the single large-scale wedge type. The dawn-dusk size of each wedgelet is ~ 600 km in the ionosphere (~ 3.2 R _E in the magnetotail, comparable to the flow channel size). We suggest that substorms have more than one type of SCW, and the composite current type is more frequent.

EMIC Waves in the Outer Magnetosphere: Observations of an Off-Equator Source Region

Electromagnetic ion cyclotron (EMIC) waves at large L shells were observed away from the magnetic equator by the Magnetospheric MultiScale (MMS) mission nearly continuously for over four hours on 28 October 2015. During this event, the wave Poynting vector direction systematically changed from parallel to the magnetic field (toward the equator), to bidirectional, to antiparallel (away from the equator). These changes coincide with the shift in the location of the minimum in the magnetic field in the southern hemisphere from poleward to equatorward of MMS. The local plasma conditions measured with the EMIC waves also suggest that the outer magnetospheric region sampled during this event was generally unstable to EMIC wave growth. Together, these observations indicate that the bidirectionally propagating wave packets were not a result of reflection at high latitudes but that MMS passed through an off-equator EMIC wave source region associated with the local minimum in the magnetic field.

Assessing the global Alfvén wave power flow into and out of the auroral acceleration region during geomagnetic storms

Geomagnetic storms are large space weather events with potentially tremendous societal implications. During these storms, the transfer of energy from the solar wind into geospace is largely increased, leading to enhanced energy flow and deposition within the magnetosphere and ionosphere. While various energy forms participate, the rate of total Alfvén wave energy flowing into the auroral acceleration region-where the magnetosphere and ionosphere couple-has not been quantified. Here, we report a fourfold increase in hemispherical Alfvénic power (from 2.59 to 10.05 GW) over a largely expanded oval band covering all longitudes and latitudes between 50° and 85° during the main storm phase compared with nonstorm periods. The Poynting flux associated with individual Alfvén waves reached values of up to about 0.5 W/m² (mapped to ionospheric altitude). These results demonstrate that Alfvén waves are an important component of geomagnetic storms and associated energy flow into the auroral acceleration region.

Transport and Loss of Ring Current Electrons Inside Geosynchronous Orbit During the 17 March 2013 Storm

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 R E R E, which indicates that the general dynamics of the electrons between 4.5 R E 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.

The Largest Electron Differential Energy Flux Observed at Mars by the Mars Express Spacecraft, 2004-2016

The goal of this paper is to understand the processes by which solar wind electrons are energized in the Martian magnetosphere and how this compares to processes at Venus and Earth. Each is unique in the source of its magnetic field topology and how this influences electron energization. To achieve this goal, 24 million spectra spanning 13 years have been examined using the electron spectrometer from the Mars Express spacecraft between about 12,000 km and about 250 km altitude, and from all latitudes and local times. The top 10 largest differential energy flux at energies above the differential energy flux peak have been found: seven spectra from the magnetosheath near noon, three from the dark tail (the largest two from the middle and ionospheric edge of the magnetosheath). Spectral comparisons show a decade range in the peak of the electron distributions; however, all distributions show a similar energy maximum dictated by solar wind/planet interaction. Similarly derived, the largest Venus spectrum occurred near the magnetosheath bow shock and had the same shape as the most intense Mars inner magnetosheath spectrum. The Mars and Venus dayside spectra compared to the Mars nightside spectrum that included an enhanced optical signal attributed to discrete "auroral" precipitation show a similar shape. These spectra are also compared to a selected auroral zone electron spectra from the Earth. The Mars and Venus results suggest that there is no more energy needed to generate electrons forming the nightside precipitation than is gained during the solar wind/planet interaction.

Modular model for Mercury's magnetospheric magnetic field confined within the average observed magnetopause

Accurate knowledge of Mercury's magnetospheric magnetic field is required to understand the sources of the planet's internal field. We present the first model of Mercury's magnetospheric magnetic field confined within a magnetopause shape derived from Magnetometer observations by the MErcury Surface, Space ENvironment, GEochemistry, and Ranging spacecraft. The field of internal origin is approximated by a dipole of magnitude 190 nT R_M³, where R_M is Mercury's radius, offset northward by 479 km along the spin axis. External field sources include currents flowing on the magnetopause boundary and in the cross-tail current sheet. The cross-tail current is described by a disk-shaped current near the planet and a sheet current at larger (≳ 5 R_M ) antisunward distances. The tail currents are constrained by minimizing the root-mean-square (RMS) residual between the model and the magnetic field observed within the magnetosphere. The magnetopause current contributions are derived by shielding the field of each module external to the magnetopause by minimizing the RMS normal component of the magnetic field at the magnetopause. The new model yields improvements over the previously developed paraboloid model in regions that are close to the magnetopause and the nightside magnetic equatorial plane. Magnetic field residuals remain that are distributed systematically over large areas and vary monotonically with magnetic activity. Further advances in empirical descriptions of Mercury's magnetospheric external field will need to account for the dependence of the tail and magnetopause currents on magnetic activity and additional sources within the magnetosphere associated with Birkeland currents and plasma distributions near the dayside magnetopause.