Abstract:The Earth's inner magnetosphere is a very dynamic system, mostly driven by the external solar wind forcing exerted upon the magnetic field of our planet. Disturbances in the solar wind, such as coronal mass ejections and co-rotating interaction regions, cause geomagnetic storms, which lead to prominent changes in charged particle populations of the inner magnetosphere - the plasmasphere, ring current, and radiation belts. Satellites operating in the regions of elevated energetic and relativistic electron fluxes can be damaged by deep dielectric or surface charging during severe space weather events. Predicting the dynamics of the charged particles and mitigating their effects on the infrastructure is of particular importance, due to our increasing reliance on space technologies. The dynamics of particles in the plasmasphere, ring current, and radiation belts are strongly coupled by means of collisions and collisionless interactions with electromagnetic fields induced by the motion of charged particles. Multidimensional numerical models simplify the treatment of transport, acceleration, and loss processes of these particles, and allow us to predict how the near-Earth space environment responds to solar storms. The models inevitably rely on a number of simplifications and assumptions that affect model accuracy and complicate the interpretation of the results. In this dissertation, we quantify the processes that control electron dynamics in the inner magnetosphere, paying particular attention to the uncertainties of the employed numerical codes and tools. We use a set of convenient analytical solutions for advection and diffusion equations to test the accuracy and stability of the four-dimensional Versatile Electron Radiation Belt (VERB-4D) code. We show that numerical schemes implemented in the code converge to the analytical solutions and that the VERB-4D code demonstrates stable behavior independent of the assumed time step. The order of the numerical scheme for the convection equation is demonstrated to affect results of ring current and radiation belt simulations, and it is crucially important to use high-order numerical schemes to decrease numerical errors in the model. Using the thoroughly tested VERB-4D code, we model the dynamics of the ring current electrons during the 17 March 2013 storm. The discrepancies between the model and observations above 4.5 Earth's radii can be explained by uncertainties in the outer boundary conditions. Simulation results indicate that the electrons were transported from the geostationary orbit towards the Earth by the global-scale electric and magnetic fields. We investigate how simulation results depend on the input models and parameters. The model is shown to be particularly sensitive to the global electric field and electron lifetimes below 4.5 Earth's radii. The effects of radial diffusion and subauroral polarization streams are also quantified. We developed a data-assimilative code that blends together a convection model of energetic electron transport and loss and Van Allen Probes satellite data by means of the Kalman filter. We show that the Kalman filter can correct model uncertainties in the convection electric field, electron lifetimes, and boundary conditions. It is also demonstrated how the innovation vector - the difference between observations and model prediction - can be used to identify physical processes missing in the model of energetic electron dynamics. We computed radial profiles of phase space density of ultrarelativistic electrons, using Van Allen Probes measurements. We analyze the shape of the profiles during geomagnetically quiet and disturbed times and show that the formation of new local minimums in the radial profiles coincides with the ground observations of electromagnetic ion-cyclotron (EMIC) waves. This correlation indicates that EMIC waves are responsible for the loss of ultrarelativistic electrons from the heart of the outer radiation belt into the Earth's atmosphere.

Numerical Modeling of Energetic Charged-particle Transport with SPECTRUM Software: General Approach and Artificial Effects due to Field Discretization

Application of the Monte Carlo Method in Modeling Transport and Acceleration of Solar Energetic Particles

Segmental Interpolating Spectra for Solar Particle Events and in Situ Validation

Implementing accelerated particle beams in a 3D simulation of the quiet Sun

Spectral element method for 3-D controlled-source electromagnetic forward modelling using unstructured hexahedral meshes

Hybrid particle-spectral method for kinetic plasma simulations

Efficient Modeling of Heavy Cosmic Rays Propagation in Evolving Astrophysical Environments

Geometric Particle-In-Cell discretizations of a plasma hybrid model with kinetic ions and mass-less fluid electrons

Effects of dispersive wave modes on charged particles transport

Modeling and understanding dynamics of charged particles in the Earth's inner magnetosphere

Monte Carlo simulations of charged particle transport in plasmas: A fast single-scattering model

Mitigation of Numerical Cerenkov Radiation and Instability Using a Hybrid Finite Difference-Fft Maxwell Solver and a Local Charge Conserving Current Deposit

Prompt: Probability-Conserved Cross Section Biasing Monte Carlo Particle Transport System

Hybrid Numerical Simulations of Pulsar Magnetospheres

Towards Modelling AR Sco: Generalised Particle Dynamics and Strong Radiation-Reaction Regimes

Coupled MHD-Focused Transport Simulations for Modeling Solar Particle Events

Multi-physics and Multi-scale Electromagnetic Modeling and High Performance Algorithms

A Kinetic-magnetohydrodynamic Model with Adaptive Mesh Refinement for Modeling Heliosphere Neutral-plasma Interaction

Numerical integration of particle orbits in discontinuous fields using VENUS-LEVIS and SPEC

Numerical Investigation of Non-equilibrium Electron Effects on the Collisional Ionization Rate in the Collisional-Radiative Model