Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasmasphere modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

Simulating Long–Term Dynamics of Radiation Belt Electrons Using DREAM3D Model

Here, we compared the performance of DREAM3D simulations in reproducing the long-term radiation belt dynamics observed by Van Allen Probes over the entire year of 2017 with various boundary conditions (BCs) and model inputs. Specifically, we investigated the effects of three different outer boundary conditions, two different low-energy boundary conditions for seed electrons, four different radial diffusion (RD) coefficients (DLL), four hiss wave models, and two chorus wave models from the literature. Using the outer boundary condition driven by GOES data, our benchmark simulation generally well reproduces the observed radiation belt dynamics inside L* = 6, with a better model performance at lower μ than higher μ, where μ is the first adiabatic invariant. By varying the boundary conditions and inputs, we find that: (a) The data-driven outer boundary condition is critical to the model performance, while adding in the data-driven seed population doesn't further improve the performance. (b) The model shows comparable performance with DLL from Brautigam and Albert (2000, https://doi.org/10.1029/1999ja900344), Ozeke et al. (2014, https://doi.org/10.1002/2013ja019204), and Liu et al. (2016, https://doi.org/10.1002/2015gl067398), while with DLL from Ali et al. (2016, https://doi.org/10.1002/2016ja023002) the model shows less RD compared to data. (c) The model performance is similar with data-based hiss models, but the results show faster loss is still needed inside the plasmasphere. (d) The model performs similarly with the two different chorus models, but better capturing the electron enhancement at higher μ using the Wang et al. (2019, https://doi.org/10.1029/2018ja026183) model due to its stronger wave power, since local heating for higher energy electrons is under-reproduced in the current model.

79 ASTRONOMY AND ASTROPHYSICS↗

Comparison of electron loss models in the inner magnetosphere during the 2013 St. Patrick's Day geomagnetic storm

Electrons with energies in the keV range play an important role in the dynamics of the inner magnetosphere. Therefore, accurately modeling electron fluxes in this region is of great interest. However, these calculations constitute a challenging task since the lifetimes of electrons that are available have limitations. In this study, we simulate electron fluxes in the energy range of 20 eV-100 keV to assess how well different electron loss models can account for the observed electron fluxes during the GEM Challenge Event of the 2013 St. Patrick's Day storm. Three models (Case 1, Case 2, and Case 3) of electron lifetimes due to wave-induced pitch angle scattering are used to compute the fluxes, which are compared with measurements from the Van Allen Probes. The three models consider electron losses due to interactions with whistler mode hiss waves inside the plasmasphere and with whistler mode chorus waves outside the plasmasphere. The Case 1 (historical) model produces excessive loss at low L shells before and after the storm, suggesting that it overestimates losses due to hiss during quiet times. During the storm main phase and early recovery all three models show good agreement with the observations, indicating that losses due to chorus during disturbed times are in general well accounted for by the models. Furthermore, the more recent Case 2 and Case 3 models show overall better agreement with the observed fluxes.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulations of Van Allen Probes Plasmaspheric Electron Density Observations

Here we simulate equatorial plasmaspheric electron densities using a physics-based model (Cold PLasma, CPL; used in the ring current-atmosphere interactions model) of the source and loss processes of refilling and erosion driven by empirical inputs. The performance of CPL is evaluated against in situ measurements by the Van Allen Probes (Radiation Belt Storm Probes) for two events: the 31 May to 5 June and 15 to 20 January 2013 geomagnetic storms observed in the premidnight and postmidnight magnetic local time (MLT) sectors, respectively. Overall, CPL reproduces the radial extent of the plasmasphere to within a mean absolute difference of |ΔL avg |=0.5 ± 0.1 L. The model electric field responsible for E × B convection and the parameterization of geomagnetic conditions (under the Kp-index and solar wind properties) implemented by CPL did not account for localized enhancements in the duskward electric field during increased activity. Rather, it was found to be largely dependent on the measure of the quiet time background. Finally, this property indicates that the agreement between these simulations and observations does not account for the complete set of physical processes during extreme (strong or weak) geomagnetic conditions impacting the plasmasphere. Nevertheless, at the presented resolution of the model CPL does provide good agreement in reproducing Radiation Belt Storm Probes observations of plasmaspheric density and plasmapause location.

79 ASTRONOMY AND ASTROPHYSICS↗

Assessment of Variability of the TEC in the Equatorial Anomaly Region with a Focus over Africa Using Rz and F10.7 as Input Drivers

This paper mainly intends to address the assessment of altitudinal variability of the vertical total electron content (VTEC) extracted from the Global Positioning System (GPS), the new versions of the International Reference Ionosphere Extended to the Plasmasphere (IRI-Plas 2017), and a quick-run ionospheric electron density (NeQuick 2) model in the equatorial anomaly region with a focus over Africa. This is manipulated employing the daily sunspot number (Rz) and 10.7 cm solar radio flux (F10.7) indices as the driver for the models during 2013–2016. The results show that the GPS-derived VTEC values are generally smaller than the modeled VTEC (IRI-Plas 2017 and NeQuick 2 VTEC) values, especially while utilizing the models with the Rz index. The modeled VTEC values are closer to the GPS VTEC values while using the models with the F10.7 index than Rz. The VTEC values obtained while using both models with the F10.7 index tend to increase while traversing from the high to the low solar activity years (especially in the March equinox and December solstice). The top-side ionospheric and plasmaspheric layers contribute the highest and lowest values to the total VTEC. Hence, the signal propagation through the ionosphere is largely affected when the signal crosses the top-side layer. Unlike other layers, due to the limited impact of the EUV on the ionization of neutral particles in the plasmasphere, the VTEC in the plasmaspheric layer (ECpl) shows smooth pattern with similar hourly values. Moreover, the IRI-Plas 2017 model does not effectively respond to the geomagnetic storm time variability of the VTEC with altitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling of the cold electron plasma density for radiation belt physics

This review focusses strictly on existing plasma density models, including ionospheric source models, empirical density models, physics-based and machine-learning density models. This review is framed in the context of radiation belt physics and space weather codes. The review is limited to the most commonly used models or to models recently developed and promising. A great variety of conditions is considered such as the magnetic local time variation, geomagnetic conditions, ionospheric source regions, radial and latitudinal dependence, and collisional vs. collisionless conditions. These models can serve to complement satellite observations of the electron plasma density when data are lacking, are for most of them commonly used in radiation belt physics simulations and can improve our understanding of the plasmasphere dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗

Scattering by whistler-mode waves during a quiet period perturbed by substorm activity

We study the dynamics of radiation belt electrons during a 10-day quiet period perturbed by substorm activity and preceding a high-speed stream (HSS), aiming at a global description of the radiation belts in L-shell, L in [2, 6], and energy [0.1, 10] MeV. We combine Van Allen Probes observations and Fokker-Planck numerical simulations of pitch-angle diffusion. The Fokker-Planck model uses event-driven pitch angle diffusion coefficients from whistler-mode waves, built from the wave properties and the ambient plasma density measurements from the Van Allen Probes. We first find this event has some similar characteristics to regular quiet times previously studied; a widely extended plasmasphere within which we observe strong and varying whistler-mode waves. These ambient conditions lead to strong pitch-angle scattering, which contributes to the creation of a wide slot region as well as a significant decay of the outer radiation belts, which are observed and qualitatively well simulated. In addition, we find the substorm activity causes short duration (within ± 4h) decay of the plasma density and a lowering amplitude of the whistler-mode waves within the plasmasphere, both causing opposite effects in terms of pitch angle diffusion. This leads to a diminution of pitch-angle diffusion at the time of the main substorm activity. Conversely, whistler-mode waves become enhanced in the time periods between the substorm injections. All effects cumulated, we find an enhancement of pitch angle diffusion by whistler-mode waves above L~4.7 during the 10-day period. This directly relates to the combination of quietness and substorm activity which allows pitch angle diffusing of up to 1 MeV electrons in the outer belt. Relativistic electrons of 1–2 MeV remain trapped in the outer belt, from L~4.7 to L~5.2, forming, in both the observations and the simulations, a distinct pocket of remnant electrons.

79 ASTRONOMY AND ASTROPHYSICS↗

Generation of Subauroral Longitudinally Extended Emissions Following Intensifications of the Poleward Boundary of the Substorm Bulge and Streamer Production

Observations showing the development of Subauroral Longitudinally Extended Emigerssions on a global scale are presented. It is demonstrated that they occur as a separation of an east-west arc-like band of luminosity that detaches away from the equatorward edge of the auroral distribution following episodes of auroral streamer production. They persist for time intervals on the order of ~30 min and devolve into patchy disjointed segments before they fade. Emissions are seen in both 557.7 nm OI and 391.4 nm N + 2 1NG lines, and 630.0 nm emissions are often observed equatorward. The 391.4 nm emissions are typically weaker and fade away more quickly than the 557.7 nm emissions, suggesting that an auroral precipitation source is initially present but is rapidly depleted as the forms age. All cases are associated with enhanced Subauroral Polarization Stream (SAPS) flows, and one event shows clear association with large-scale Giant Undulations (GUs) and the formation of spur-like forms that fold equatorward under the main oval. A model is proposed for the production of subauroral longitudinally extended forms in which nonlinear growth of SAPS-induced surface waves on the plasmapause results in a disruption of the boundary. It is suggested that hot plasma of plasma sheet origin becomes entrained in the plasmasphere to produce transient precipitation-associated auroral emissions that may decay into STEVE emissions and that cold dense plasma from the plasmasphere becomes entrained onto open drift paths to feed long-lived drainage plumes. This process may occur quasi periodically during intense substorms and storm-time conditions.

79 ASTRONOMY AND ASTROPHYSICS↗

Key Elements of Auroral Substorm Development and their Relationship to Recent Observations of Detached Sub-Auroral Phenomena Including STEVE-like Emissions

Since Akasofu introduced the concept of the substorm in 1964, numerous refinements have been added to the phenomenological model of how auroral substorms develop. Here, we review several of the most important of these including; polar cap patches, intensifications of the Polar Activated Band (PAB) at the poleward edge of the bulge, the development of new arcs poleward, the generation of auroral streamers, the evolution of streamers into auroral torches and omega bands, the association of streamers with particle injection and BBF activity in the tail, and wedgelet-like magnetic perturbations on the ground. We also review the concept of “contact breakups” which are breakups that appear to be triggered by the arrival of streamer-like forms in the equatorward regions of the oval. Finally, we remind the community that the addition of detached sub-auroral emissions was made to the phenomenological picture of substorm develop in the late 1960s along with the concept of auroral streamers. Newer observations of detached sub-auroral (STEVE-like emissions) from the Viking/UVI and POLAR/VIS imagers are presented which confirm that they are east-west aligned bands of emission that separate away from the equatorward edge of the auroral oval in response to intensifications of the poleward edge of the bulge and subsequent streamer production. We propose a new model for the formation of detached sub-auroral (STEVE-like) emissions in which a non-linear growth of a SAPS-driven instability (e.g. like the shear-flow/ballooning instability) at the plasmapause results in the disruption of the boundary separating hot plasma sheet particle from the cold plasmaspheric plasma. It is proposed that the resulting intermixing of plasma populations leads to both the observed STEVE-like emissions and also provides a source of cold plasma on open drift paths that can feed the long-lived drainage plumes that have recently been discovered.

58 GEOSCIENCES↗

Propagation of Chorus Waves Generated in Minimum‐B Pockets

Abstract We present the first ray tracing results of chorus waves generated in minimum‐B pockets. A new ray tracer is developed to accommodate an arbitrary magnetic field model, which enables us to construct minimum‐B pocket geometries using the Tsyganenko (1989) model. Rays with different frequencies and initial wave normal angles are launched from minimum‐B pockets and traced under different geomagnetic activity and hot electron density conditions. Our results indicate that chorus waves generated from minimum‐B pockets are usually highly localized and thus are unlikely to be detected unless the spacecraft is very close to the source region. These waves are also very unlikely to propagate into the plasmasphere. The propagation distances generally decrease with increasing geomagnetic activity.

Kang, Ning↗

The Composition of Plasma inside Geostationary Orbit Based on Van Allen Probes Observations

The composition of the inner magnetosphere is of great importance for determining the plasma pressure, and thus the currents and magnetic field configuration. In this study, we perform a statistical survey of equatorial plasma pressure distributions and investigate the relative contributions of ions and electron with different energies inside of geostationary orbit under two AE levels based on over sixty months of observations from the HOPE and RBSPICE mass spectrometers on board Van Allen Probes. We find that the total and partial pressures of different species increase significantly at high AE levels with Hydrogen (H + ) pressure being dominant in the plasmasphere. The pressures of the heavy ions and electrons increase outside the plasmapause and develop a strong dawn-dusk asymmetry with ion pressures peaking at dusk and electron pressure peaking at dawn. In addition, ring current H + with energies ranging from 50 keV up to several hundred keV is the dominant component of plasma pressure during both quiet (> 90%) and active times (> 60%), while Oxygen (O + ) with 10 < E < 50 keV and electrons with 0.1 < E < 40 keV become important during active times contributing more than 25% and 20% on the nightside, respectively, while the Helium (He + ) contribution is generally small. The results presented in this study provide a global picture of the equatorial plasma pressure distributions and the associated contributions from different species with different energy ranges, which advance our knowledge of wave generation and provide models with a systematic baseline of plasma composition.

79 ASTRONOMY AND ASTROPHYSICS↗

The Cold Ion Population at Geosynchronous Orbit and Transport to the Dayside Magnetopause: September 2015 to February 2016

During intervals of enhanced magnetospheric convection, a high-density plume of cold ions is eroded from the plasmasphere and can flow toward the dayside magnetopause where it has the potential to reduce the rate of magnetic reconnection. In any interval of long-duration enhanced convection, tons of ions may follow such a trajectory. The study here concerns cold ion observations from geosynchronous orbit (GEO) during both calm and active periods. Probability distributions of the cold ion density and cold ion flow speed are determined during the 6-month period from September 2015 to February 2016, inclusive. During low geomagnetic activity, the cold, dense ions are in corotation with the Earth, and flow speeds seldom exceed 10 km/s. During elevated geomagnetic activity, the cold ions between 12 and 18 magnetic local time are observed to flow toward the dayside magnetosphere with a speed of >10 km/s approximately 50% of the time. The Shue et al. (1998, https://doi.org/10.1029/98JA01103) model of the magnetopause location is used to derive the distribution of approximate minimum times for the cold ions to be transported from GEO to the dayside magnetopause. On average, during enhanced convection periods (Kp > 3) ions will take a mean time of ~4.5 hr to travel from GEO to the dayside magnetopause.

79 ASTRONOMY AND ASTROPHYSICS↗