-
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.
-
Space climate and space weather over the past 400 years: 2. Proxy indicators of geomagnetic storm and substorm
Authors:
Mike Lockwood,
Mathew J. Owens,
Luke A. Barnard,
Chris J. Scott,
Clare E. Watt,
Sarah Bentley
Abstract:
Using the reconstruction of power input to the magnetosphere given in Paper 1 (arXiv:1708.04904), we reconstruct annual means of geomagnetic indices over the past 400 years to within a 1-sigma error of +/-20 pc. In addition, we study the behaviour of the lognormal distribution of daily and hourly values about these annual means and show that we can also reconstruct the fraction of geomagnetically-…
▽ More
Using the reconstruction of power input to the magnetosphere given in Paper 1 (arXiv:1708.04904), we reconstruct annual means of geomagnetic indices over the past 400 years to within a 1-sigma error of +/-20 pc. In addition, we study the behaviour of the lognormal distribution of daily and hourly values about these annual means and show that we can also reconstruct the fraction of geomagnetically-active (storm-like) days and (substorm-like) hours in each year to accuracies of 50-60 pc. The results are the first physics-based quantification of the space weather conditions in both the Dalton and Maunder minima. We predict terrestrial disturbance levels in future repeats of these minima, allowing for the weakening of Earth's dipole moment.
△ Less
Submitted 8 January, 2018;
originally announced January 2018.
-
Space Climate and Space Weather over the past 400 years: 1. The Power input to the Magnetosphere
Authors:
Mike Lockwood,
Mathew J. Owens,
Luke A. Barnard,
Chris J. Scott,
Clare E. Watt
Abstract:
Using information on geomagnetic activity, sunspot numbers and cosmogenic isotopes, supported by historic eclipse images and in conjunction with models, it has been possible to reconstruct annual means of solar wind speed and number density and heliospheric magnetic field (HMF) intensity since 1611, when telescopic observations of sunspots began. These models are developed and tuned using data rec…
▽ More
Using information on geomagnetic activity, sunspot numbers and cosmogenic isotopes, supported by historic eclipse images and in conjunction with models, it has been possible to reconstruct annual means of solar wind speed and number density and heliospheric magnetic field (HMF) intensity since 1611, when telescopic observations of sunspots began. These models are developed and tuned using data recorded by near-Earth interplanetary spacecraft and by solar magnetograms over the past 53 years. In this paper, we use these reconstructions to quantify power input into the magnetosphere over the past 400 years. For each year, both the annual mean power input is computed and its distribution in daily means. This is possible because the distribution of daily values divided by the annual mean is shown to maintain the same lognormal form with a constant variance. This study is another important step towards the development of a physics-based, long-term climatology of space weather conditions.
△ Less
Submitted 17 August, 2017; v1 submitted 16 August, 2017;
originally announced August 2017.