Two methods to analyse radial diffusion ensembles: the peril of space- and time- dependent diffusion
Authors:
Sarah N. Bentley,
J. Stout,
Rhys Thompson,
Daniel J. Ratliff,
Clare E. J. Watt
Abstract:
Particle dynamics in Earth's outer radiation belt can be modelled using a diffusion framework, where large-scale electron movements are captured by a diffusion equation across a single adiabatic invariant, $L^{*}$ $``(L)"$. While ensemble models are promoted to represent physical uncertainty, as yet there is no validated method to analyse radiation belt ensembles. Comparisons are complicated by th…
▽ More
Particle dynamics in Earth's outer radiation belt can be modelled using a diffusion framework, where large-scale electron movements are captured by a diffusion equation across a single adiabatic invariant, $L^{*}$ $``(L)"$. While ensemble models are promoted to represent physical uncertainty, as yet there is no validated method to analyse radiation belt ensembles. Comparisons are complicated by the domain dependent diffusion, since diffusion coefficient $D_{LL}$ is dependent on $L$. We derive two tools to analyse ensemble members: time to monotonicity $t_m$ and mass/energy moment quantities $\mathcal{N}, \mathcal{E}$. We find that the Jacobian ($1/L^2$) is necessary for radiation belt error metrics. Components of $\partial\mathcal{E}/\partial t$ are explicitly calculated to compare the effects of outer and inner boundary conditions, and loss, on the ongoing diffusion. Using $t_m$, $\mathcal{N}$ and $\mathcal{E}$, we find that: (a) different physically motivated choices of outer boundary condition and location result in different final states and different rates of evolution; (b) the gradients of the particle distribution affect evolution more significantly than $D_{LL}$; (c) the enhancement location, and the amount of initial background particles, are both significant factors determining system evolution; (d) loss from pitch-angle scattering is generally dominant; it mitigates but does not remove the influence of both initial conditions and outer boundary settings, which are due to the $L$-dependence of $D_{LL}$. We anticipate this study will promote renewed focus on the distribution gradients, on the location and nature of the outer boundary in radiation belt modelling, and provide a foundation for systematic ensemble modelling.
△ Less
Submitted 5 July, 2024;
originally announced July 2024.
Statistics of Solar Wind Electron Breakpoint Energies Using Machine Learning Techniques
Authors:
Mayur R. Bakrania,
I. Jonathan Rae,
Andrew P. Walsh,
Daniel Verscharen,
Andy W. Smith,
Téo Bloch,
Clare E. J. Watt
Abstract:
Solar wind electron velocity distributions at 1 au consist of a thermal "core" population and two suprathermal populations: "halo" and "strahl". The core and halo are quasi-isotropic, whereas the strahl typically travels radially outwards along the parallel and/or anti-parallel direction with respect to the interplanetary magnetic field. With Cluster-PEACE data, we analyse energy and pitch angle d…
▽ More
Solar wind electron velocity distributions at 1 au consist of a thermal "core" population and two suprathermal populations: "halo" and "strahl". The core and halo are quasi-isotropic, whereas the strahl typically travels radially outwards along the parallel and/or anti-parallel direction with respect to the interplanetary magnetic field. With Cluster-PEACE data, we analyse energy and pitch angle distributions and use machine learning techniques to provide robust classifications of these solar wind populations. Initially, we use unsupervised algorithms to classify halo and strahl differential energy flux distributions to allow us to calculate relative number densities, which are of the same order as previous results. Subsequently, we apply unsupervised algorithms to phase space density distributions over ten years to study the variation of halo and strahl breakpoint energies with solar wind parameters. In our statistical study, we find both halo and strahl suprathermal breakpoint energies display a significant increase with core temperature, with the halo exhibiting a more positive correlation than the strahl. We conclude low energy strahl electrons are scattering into the core at perpendicular pitch angles. This increases the number of Coulomb collisions and extends the perpendicular core population to higher energies, resulting in a larger difference between halo and strahl breakpoint energies at higher core temperatures. Statistically, the locations of both suprathermal breakpoint energies decrease with increasing solar wind speed. In the case of halo breakpoint energy, we observe two distinct profiles above and below 500 km/s. We relate this to the difference in origin of fast and slow solar wind.
△ Less
Submitted 7 July, 2020; v1 submitted 26 May, 2020;
originally announced May 2020.