{"id":2415,"date":"2026-06-09T13:37:40","date_gmt":"2026-06-09T12:37:40","guid":{"rendered":"https:\/\/research.northumbria.ac.uk\/drbentley\/?p=2415"},"modified":"2026-06-09T13:37:40","modified_gmt":"2026-06-09T12:37:40","slug":"what-did-we-learn-from-radial-diffusion-ensembles","status":"publish","type":"post","link":"https:\/\/research.northumbria.ac.uk\/drbentley\/?p=2415","title":{"rendered":"What did we learn from radial diffusion ensembles?"},"content":{"rendered":"\n<p class=\"wp-block-paragraph\"><em>I wrote a post about this paper a year ago &#8211; embarrassingly, it appears I never pressed &#8220;Publish&#8221;. I&#8217;ve updated it for you too.<\/em><\/p>\n\n\n\n<p class=\"has-medium-font-size wp-block-paragraph\">One of the current challenges in space weather modelling is to capture the variety of observations in our physical models. In meteorology, one such technique is ensemble modelling &#8211; we applied this to simple radiation belt models, to build a foundation for future ensembles in forecasting. Our paper \u201c<em><a href=\"https:\/\/doi.org\/10.1063\/5.0217268\" data-type=\"link\" data-id=\"https:\/\/doi.org\/10.1063\/5.0217268\">Two methods to analyze radial diffusion ensembles: The perils of space- and time-dependent diffusion<\/a>\u201d<\/em> was published in <em>Physics of Plasmas<\/em> as a Featured Article. This 4+ year journey led us to new techniques and to several unexpected and intriguing results that throw up some questions about everyday choices in radiation belt diffusion modelling. These mainly arise because we are modelling in the <math data-latex=\"L^*\"><semantics><msup><mi>L<\/mi><mo>\u2217<\/mo><\/msup><annotation encoding=\"application\/x-tex\">L^*<\/annotation><\/semantics><\/math> domain\u2026 and the amount of diffusion varies as <math data-latex=\"L^*\"><semantics><msup><mi>L<\/mi><mo>\u2217<\/mo><\/msup><annotation encoding=\"application\/x-tex\">L^*<\/annotation><\/semantics><\/math> varies.<\/p>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"has-medium-font-size wp-block-paragraph\">Our two key results explained here:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li class=\"has-medium-font-size\">Both the location \/ distance of the outer boundary, and the outer boundary condition change the way in which the system evolves &#8211; but there is no clear \u201cmost physical\u201d choice. Observations can help limit this for hindcasting, but&#8230; <strong>what should we use?<\/strong><\/li>\n\n\n\n<li class=\"has-medium-font-size\">In our experiments, the gradient of the phase space density (<math data-latex=\"PSD\"><semantics><mrow><mi>P<\/mi><mi>S<\/mi><mi>D<\/mi><\/mrow><annotation encoding=\"application\/x-tex\">PSD<\/annotation><\/semantics><\/math>) contributed more to the evolution of the system, rather than the diffusion coefficient  <math data-latex=\"D_{LL}\"><semantics><msub><mi>D<\/mi><mrow><mi>L<\/mi><mi>L<\/mi><\/mrow><\/msub><annotation encoding=\"application\/x-tex\">D_{LL}<\/annotation><\/semantics><\/math>. <strong>If true, should we focus on determining the <\/strong><math data-latex=\"PSD\"><semantics><mrow><mi>P<\/mi><mi>S<\/mi><mi>D<\/mi><\/mrow><annotation encoding=\"application\/x-tex\">PSD<\/annotation><\/semantics><\/math><strong> than on <math data-latex=\"D_{LL}\"><semantics><msub><mi>D<\/mi><mrow><mi>L<\/mi><mi>L<\/mi><\/mrow><\/msub><annotation encoding=\"application\/x-tex\">D_{LL}<\/annotation><\/semantics><\/math>? If it\u2019s not true, are 1d-radial diffusion models inappropriate radiation belt toy models?<\/strong><\/li>\n<\/ul>\n<\/blockquote>\n\n\n\n<h1 class=\"wp-block-heading\">The model<\/h1>\n\n\n\n<p class=\"has-medium-font-size wp-block-paragraph\">In radiation belt modelling, we can approximate belt evolution by considering diffusion across three adiabatic invariants of magnetically trapped particles. The third invariant (magnetic flux enclosed on a drift path around the Earth) is much longer than the other timescales and can be modelled by itself as an approximation of the belts. We use a Fokker-Planck model,<\/p>\n\n\n\n<div class=\"wp-block-math\"><math display=\"block\"><semantics><mrow><mfrac><mrow><mi>\u2202<\/mi><mi>f<\/mi><\/mrow><mrow><mi>\u2202<\/mi><mi>t<\/mi><\/mrow><\/mfrac><mo>=<\/mo><msup><mi>L<\/mi><mn>2<\/mn><\/msup><mrow><mo fence=\"true\" form=\"prefix\">(<\/mo><mfrac><msub><mi>D<\/mi><mrow><mi>L<\/mi><mi>L<\/mi><\/mrow><\/msub><msup><mi>L<\/mi><mn>2<\/mn><\/msup><\/mfrac><mfrac><mrow><mi>\u2202<\/mi><mi>f<\/mi><\/mrow><mrow><mi>\u2202<\/mi><mi>L<\/mi><\/mrow><\/mfrac><mo fence=\"true\" form=\"postfix\">)<\/mo><\/mrow><mo>\u2212<\/mo><mi class=\"mathcal\">\u2112<\/mi><mi>f<\/mi><\/mrow><annotation encoding=\"application\/x-tex\">\\frac{\\partial f}{\\partial t} = L^2 \\left( \\frac{D_{LL}}{L^2}\\frac{\\partial f}{\\partial L} \\right) &#8211; \\mathcal{L} f<\/annotation><\/semantics><\/math><\/div>\n\n\n\n<p class=\"has-medium-font-size wp-block-paragraph\">where f is the phase space density,  <math data-latex=\"L^*\"><semantics><msup><mi>L<\/mi><mo>\u2217<\/mo><\/msup><annotation encoding=\"application\/x-tex\">L^*<\/annotation><\/semantics><\/math> is related to the third adiabatic invariant and as a rough guide can be considered &#8220;how far the drift path is from the Earth, in Earth radii&#8221;, <math data-latex=\"D_{LL} = D_0L^*\"><semantics><mrow><msub><mi>D<\/mi><mrow><mi>L<\/mi><mi>L<\/mi><\/mrow><\/msub><mo>=<\/mo><msub><mi>D<\/mi><mn>0<\/mn><\/msub><msup><mi>L<\/mi><mo>\u2217<\/mo><\/msup><\/mrow><annotation encoding=\"application\/x-tex\">D_{LL} = D_0L^*<\/annotation><\/semantics><\/math> is the diffusion coefficient and <math data-latex=\"\\mathcal{L}\"><semantics><mi class=\"mathcal\">\u2112<\/mi><annotation encoding=\"application\/x-tex\">\\mathcal{L}<\/annotation><\/semantics><\/math> incorporates loss.<\/p>\n\n\n\n<h1 class=\"wp-block-heading\">Major results<\/h1>\n\n\n\n<h2 class=\"wp-block-heading\">New tools to measure the difference between ensemble members<\/h2>\n\n\n\n<p class=\"has-medium-font-size wp-block-paragraph\">We can&#8217;t compare model results to any ground truth; nevertheless, in any radiation belt model a variety of metrics are used to compare a variety of outputs. We found that:<\/p>\n\n\n\n<h3 class=\"wp-block-heading\"><em>The shape of the distribution is a useful and representative measure<\/em><\/h3>\n\n\n\n<p class=\"has-medium-font-size wp-block-paragraph\">We found that any threshold-based methods to find when the simulation was \u201cdone\u201d (e.g. these two distributions are the same when mean square error is less than X) were completely dependent on the (arbitrary!) choice of threshold. Instead, we chose a measure of the diffusion that represents what we\u2019re often interested in &#8211; how long for radial diffusion to reduce an enhancement by diffusing it away (&#8220;<strong>time to monotonicity<\/strong>&#8220;)<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"937\" height=\"245\" src=\"https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-3.png\" alt=\"Doodle showing that once an enhancement in the phase space density has diffused away, we are left with a monotonic function.\" class=\"wp-image-2426\" srcset=\"https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-3.png 937w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-3-300x78.png 300w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-3-768x201.png 768w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-3-230x60.png 230w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-3-350x92.png 350w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-3-480x126.png 480w\" sizes=\"auto, (max-width: 937px) 100vw, 937px\" \/><figcaption class=\"wp-element-caption\">Doodle showing that once an enhancement in the phase space density has diffused away, we are left with a monotonic function.<\/figcaption><\/figure>\n\n\n\n<h3 class=\"wp-block-heading\"><em>There are extremely useful analytic approaches we have never used<\/em><\/h3>\n\n\n\n<p class=\"has-medium-font-size wp-block-paragraph\">The second metrics we used were analytic: the mass- and energy-like density<\/p>\n\n\n\n<div class=\"wp-block-math\"><math display=\"block\"><semantics><mtable displaystyle=\"true\" columnalign=\"right left\" class=\"tml-jot\"><mtr><mtd class=\"tml-right\" style=\"padding-left:0em;padding-right:0em;\"><mi class=\"mathcal\">\ud835\udca9<\/mi><\/mtd><mtd class=\"tml-left\" style=\"padding-left:0em;padding-right:0em;\"><mrow><mo>=<\/mo><mo movablelimits=\"false\">\u222b<\/mo><mi>f<\/mi><mfrac><mn>1<\/mn><msup><mi>L<\/mi><mn>2<\/mn><\/msup><\/mfrac><mi>d<\/mi><mi>L<\/mi><\/mrow><\/mtd><\/mtr><mtr><mtd class=\"tml-right\" style=\"padding-left:0em;padding-right:0em;\"><mi class=\"mathcal\">\u2130<\/mi><\/mtd><mtd class=\"tml-left\" style=\"padding-left:0em;padding-right:0em;\"><mrow><mo>=<\/mo><mfrac><mn>1<\/mn><mrow><mn>2<\/mn><msub><mi>c<\/mi><mi class=\"mathcal\">\u2130<\/mi><\/msub><\/mrow><\/mfrac><mo movablelimits=\"false\">\u222b<\/mo><msup><mrow><mo fence=\"true\" form=\"prefix\">|<\/mo><mi>f<\/mi><mo fence=\"true\" form=\"postfix\">|<\/mo><\/mrow><mn>2<\/mn><\/msup><mfrac><mn>1<\/mn><msup><mi>L<\/mi><mn>2<\/mn><\/msup><\/mfrac><mi>d<\/mi><mi>L<\/mi><\/mrow><\/mtd><\/mtr><\/mtable><annotation encoding=\"application\/x-tex\">\\begin{align*} \\mathcal{N} &amp;= \\int f \\frac{1}{L^2} dL \\\\ \\mathcal{E} &amp;= \\frac{1}{2c_\\mathcal{E}}\\int \\left| f \\right|^2 \\frac{1}{L^2} dL \\end{align*}<\/annotation><\/semantics><\/math><\/div>\n\n\n\n<p class=\"has-medium-font-size wp-block-paragraph\"><math data-latex=\"\\mathcal{N}\"><semantics><mi class=\"mathcal\">\ud835\udca9<\/mi><annotation encoding=\"application\/x-tex\">\\mathcal{N}<\/annotation><\/semantics><\/math> clearly represents the total mass in our simulation at a given time step. <math data-latex=\"\\mathcal{E}\"><semantics><mi class=\"mathcal\">\u2130<\/mi><annotation encoding=\"application\/x-tex\">\\mathcal{E}<\/annotation><\/semantics><\/math> is more interesting; it represents how far we are from a low-energy or equilibrium state. We use it to quantify how the system evolves towards equilibrium.<\/p>\n\n\n\n<h2 class=\"wp-block-heading\">New insights about the boundary<\/h2>\n\n\n\n<p class=\"has-medium-font-size wp-block-paragraph\">It&#8217;s long been understood that the outer boundary of radiation belt models is important but tricky: we only have data far inside the true boundary, which is difficult to define and difficult to find in observations. An additional question we had going in was &#8220;What about the  <math data-latex=\"L^*\"><semantics><msup><mi>L<\/mi><mo>\u2217<\/mo><\/msup><annotation encoding=\"application\/x-tex\">L^*<\/annotation><\/semantics><\/math> dependence of diffusion?&#8221; Many models in practice are cut short to a domain where data can constrain them. But if the amount of diffusion depends on where we are in <math data-latex=\"L^*\"><semantics><msup><mi>L<\/mi><mo>\u2217<\/mo><\/msup><annotation encoding=\"application\/x-tex\">L^*<\/annotation><\/semantics><\/math>, what does that imply about our abbreviated models?<\/p>\n\n\n\n<h3 class=\"wp-block-heading\"><em>Outer boundary conditions used currently are not physical<\/em><\/h3>\n\n\n\n<p class=\"has-medium-font-size wp-block-paragraph\">We tested both Neumann (fixed gradient; the radiation belt continues past the model domain) and Dirichlet (fixed value; say that data is available to constrain it) outer boundaries. Both essentially impose unphysical results.<\/p>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"396\" src=\"https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-2-1024x396.png\" alt=\"Figure adapted from Figure 4 of the paper: starting with different boundary conditions changes _how_ the enhancement diffuses away, not just how long it takes; both have unphysical implications.\" class=\"wp-image-2424\" srcset=\"https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-2-1024x396.png 1024w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-2-300x116.png 300w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-2-768x297.png 768w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-2-1000x387.png 1000w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-2-230x89.png 230w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-2-350x135.png 350w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-2-480x186.png 480w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-2.png 1148w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><figcaption class=\"wp-element-caption\">Figure adapted from Figure 4 of the paper: starting with different boundary conditions changes _how_ the enhancement diffuses away, not just how long it takes; both have unphysical implications.<\/figcaption><\/figure>\n\n\n\n<p class=\"has-medium-font-size wp-block-paragraph\">With a Neumann outer boundary the injection diffuses away quickly because the particles can move outwards and the phase space density at the outer limit can rise &#8211; but this represents a large increase outside the domain. With a Dirichlet outer boundary, this cannot happen, but nor can material flow outwards. In our experiments we found that far more ensemble members with a Neumann outer boundary reached monotonicity (where the injection diffused away) Recent work has shown that data-derived Dirichlet models which are most common are indeed the best current compromise (<a href=\"https:\/\/agupubs.onlinelibrary.wiley.com\/doi\/10.1029\/2023JA032286\" data-type=\"link\" data-id=\"https:\/\/agupubs.onlinelibrary.wiley.com\/doi\/10.1029\/2023JA032286\">Lee et al., 2024<\/a>) but these will not work when forecasting, as such data is obviously not yet available. Mixed boundary conditions will be more physical to represent the continuation of material outside the modelled domain.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\"><em>Cutting short the domain changes the physics<\/em><\/h3>\n\n\n\n<p class=\"has-medium-font-size wp-block-paragraph\">The amount of diffusion that occurs is dependent on where we are within the domain; it depends on <math data-latex=\"L^*\"><semantics><msup><mi>L<\/mi><mo>\u2217<\/mo><\/msup><annotation encoding=\"application\/x-tex\">L^*<\/annotation><\/semantics><\/math>. The rate of diffusion at <math data-latex=\"L^*=6\"><semantics><mrow><msup><mi>L<\/mi><mo>\u2217<\/mo><\/msup><mo>=<\/mo><mn>6<\/mn><\/mrow><annotation encoding=\"application\/x-tex\">L^*=6<\/annotation><\/semantics><\/math> will be higher than at <math data-latex=\"L^*=4\"><semantics><mrow><msup><mi>L<\/mi><mo>\u2217<\/mo><\/msup><mo>=<\/mo><mn>4<\/mn><\/mrow><annotation encoding=\"application\/x-tex\">L^*=4<\/annotation><\/semantics><\/math>; more material is lost when the domain is longer. However, domain length also interacts with the type of boundary condition.<\/p>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"275\" src=\"https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-1-1024x275.png\" alt=\"Doodle showing the main finding of Figure 12 from the paper: that longer domains (more representative of the real outer radiation belt) reduce the differences we find between two choices of outer boundary condition.\" class=\"wp-image-2422\" srcset=\"https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-1-1024x275.png 1024w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-1-300x81.png 300w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-1-768x206.png 768w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-1-1000x269.png 1000w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-1-230x62.png 230w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-1-350x94.png 350w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-1-480x129.png 480w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-1.png 1083w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><figcaption class=\"wp-element-caption\">Doodle showing the main finding of Figure 12 from the paper: that longer domains (more representative of the real outer radiation belt) reduce the differences we find between two choices of outer boundary condition.<\/figcaption><\/figure>\n\n\n\n<p class=\"has-medium-font-size wp-block-paragraph\">The doodles above summarise this: the energy-like quantity <math data-latex=\"\\mathcal{E}\"><semantics><mi class=\"mathcal\">\u2130<\/mi><annotation encoding=\"application\/x-tex\">\\mathcal{E}<\/annotation><\/semantics><\/math> which represents how far we are from equilibrium is different for Neumann vs Dirichlet conditions, but the longer the domain, the more similar they become.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><\/p>\n\n\n\n<h2 class=\"wp-block-heading\">Should we spend time improving the diffusion coefficient  <math data-latex=\"D_{LL}\"><semantics><msub><mi>D<\/mi><mrow><mi>L<\/mi><mi>L<\/mi><\/mrow><\/msub><annotation encoding=\"application\/x-tex\">D_{LL}<\/annotation><\/semantics><\/math> or the phase space density <math data-latex=\"PSD\"><semantics><mrow><mi>P<\/mi><mi>S<\/mi><mi>D<\/mi><\/mrow><annotation encoding=\"application\/x-tex\">PSD<\/annotation><\/semantics><\/math>?<\/h2>\n\n\n\n<p class=\"has-medium-font-size wp-block-paragraph\">Using an analytic tool such as <math data-latex=\"\\mathcal{E}\"><semantics><mi class=\"mathcal\">\u2130<\/mi><annotation encoding=\"application\/x-tex\">\\mathcal{E}<\/annotation><\/semantics><\/math> allows us to compare the relative contribution of different model components to the resulting diffusion by considering <math data-latex=\"\\mathcal{E}_t\"><semantics><msub><mi class=\"mathcal\">\u2130<\/mi><mi>t<\/mi><\/msub><annotation encoding=\"application\/x-tex\">\\mathcal{E}_t<\/annotation><\/semantics><\/math>. We find that:<\/p>\n\n\n\n<figure class=\"wp-block-image aligncenter size-medium\"><img loading=\"lazy\" decoding=\"async\" width=\"300\" height=\"227\" src=\"https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-4-300x227.png\" alt=\"Part of Figure 9 from the paper, showing that in our numerical experiment the PSD gradient dominates the evolution to an equilibrium state - not the diffusion coefficient. \" class=\"wp-image-2432\" srcset=\"https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-4-300x227.png 300w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-4-230x174.png 230w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-4-350x265.png 350w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-4-480x363.png 480w, https:\/\/research.northumbria.ac.uk\/drbentley\/wp-content\/uploads\/2026\/06\/image-4.png 725w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><figcaption class=\"wp-element-caption\">Part of Figure 9 from the paper, showing that in our numerical experiment<math data-latex=\"(f_L)^2 &gt; \\frac{D_{LL}}{L^2}\"><semantics><mrow><mo form=\"prefix\" stretchy=\"false\">(<\/mo><msub><mi>f<\/mi><mi>L<\/mi><\/msub><msup><mo form=\"postfix\" stretchy=\"false\">)<\/mo><mn>2<\/mn><\/msup><mo>&gt;<\/mo><mfrac><msub><mi>D<\/mi><mrow><mi>L<\/mi><mi>L<\/mi><\/mrow><\/msub><msup><mi>L<\/mi><mn>2<\/mn><\/msup><\/mfrac><\/mrow><annotation encoding=\"application\/x-tex\">(f_L)^2 &gt; \\frac{D_{LL}}{L^2}<\/annotation><\/semantics><\/math> and hence the PSD gradient dominates the evolution to an equilibrium state &#8211; not the diffusion coefficient. <\/figcaption><\/figure>\n\n\n\n<p class=\"has-text-align-left has-medium-font-size wp-block-paragraph\">You&#8217;ll notice that this is <em>all red<\/em>. This means that across the model domain, the gradient of the phase space density contributes more to the evolution of the particle distribution than the diffusion coefficient <math data-latex=\"D_{LL}\"><semantics><msub><mi>D<\/mi><mrow><mi>L<\/mi><mi>L<\/mi><\/mrow><\/msub><annotation encoding=\"application\/x-tex\">D_{LL}<\/annotation><\/semantics><\/math>. This directly contradicts current practice in radiation belt modelling, where the focus is on developing appropriate  <math data-latex=\"D_{LL}\"><semantics><msub><mi>D<\/mi><mrow><mi>L<\/mi><mi>L<\/mi><\/mrow><\/msub><annotation encoding=\"application\/x-tex\">D_{LL}<\/annotation><\/semantics><\/math>.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><\/p>\n\n\n\n<h1 class=\"wp-block-heading\">How significant would this be in a less idealised scenario?<\/h1>\n\n\n\n<p class=\"has-medium-font-size wp-block-paragraph\">Obviously this is a very simple model, but we need to understand how model changes affect the final ensemble. Changes in the initial phase space density and in the outer boundary location are actually (in general) dominated by loss processes, once we add them in &#8211; although if the plasmasphere causes loss at high <math data-latex=\"L^*\"><semantics><msup><mi>L<\/mi><mo>\u2217<\/mo><\/msup><annotation encoding=\"application\/x-tex\">L^*<\/annotation><\/semantics><\/math> than a simple monotonic distribution is never reached. Nevertheless, boundary conditions and location still need to be more physical when we don&#8217;t have data to constrain them &#8211; i.e. for forecasting and nowcasting.<\/p>\n\n\n\n<p class=\"has-medium-font-size wp-block-paragraph\">The results comparing relative contributions from <math data-latex=\"D_{LL}\"><semantics><msub><mi>D<\/mi><mrow><mi>L<\/mi><mi>L<\/mi><\/mrow><\/msub><annotation encoding=\"application\/x-tex\">D_{LL}<\/annotation><\/semantics><\/math> and <math data-latex=\"PSD\"><semantics><mrow><mi>P<\/mi><mi>S<\/mi><mi>D<\/mi><\/mrow><annotation encoding=\"application\/x-tex\">PSD<\/annotation><\/semantics><\/math>  is concerning. Is this an indicator suggesting a shift of focus? Conversations with modellers of the full radiation belts suggest that this is not what they see in practice. Perhaps this model is too idealised? We need to test this further.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><em>And to answer my original question, &#8220;What are realistic radial diffusion timescales for an injection?&#8221;: From hours to a week. Very helpful.<\/em><\/p>\n\n\n\n<h1 class=\"wp-block-heading\">Links:<\/h1>\n\n\n\n<ul class=\"wp-block-list\">\n<li>The paper can be found at <a href=\"https:\/\/doi.org\/10.1063\/5.0217268\">https:\/\/doi.org\/10.1063\/5.0217268<\/a><\/li>\n\n\n\n<li>A Python notebook to make figures shown in the paper can be found at <a href=\"https:\/\/zenodo.org\/doi\/10.5281\/zenodo.13336681\">https:\/\/zenodo.org\/doi\/10.5281\/zenodo.13336681<\/a><\/li>\n\n\n\n<li>I presented this at a <a href=\"https:\/\/www.youtube.com\/watch?v=UW4q09vIQK8\" data-type=\"link\" data-id=\"https:\/\/www.youtube.com\/watch?v=UW4q09vIQK8\">MIST seminar<\/a> and at a <a href=\"https:\/\/mediacentral.princeton.edu\/id\/1_5uyr4rfp\" data-type=\"link\" data-id=\"https:\/\/mediacentral.princeton.edu\/id\/1_5uyr4rfp\">Physics of Plasmas webinar<\/a>, both of which can be watched online.<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\"><\/p>\n","protected":false},"excerpt":{"rendered":"<p>I wrote a post about this paper a year ago &#8211; embarrassingly, it appears I never pressed &#8220;Publish&#8221;. I&#8217;ve updated it for [&hellip;]<\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[22],"tags":[20,17,16,14],"class_list":["post-2415","post","type-post","status-publish","format-standard","hentry","category-paper-summaries","tag-mathematics","tag-radial-diffusion","tag-radiation-belts","tag-research"],"_links":{"self":[{"href":"https:\/\/research.northumbria.ac.uk\/drbentley\/index.php?rest_route=\/wp\/v2\/posts\/2415","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/research.northumbria.ac.uk\/drbentley\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/research.northumbria.ac.uk\/drbentley\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/research.northumbria.ac.uk\/drbentley\/index.php?rest_route=\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/research.northumbria.ac.uk\/drbentley\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=2415"}],"version-history":[{"count":6,"href":"https:\/\/research.northumbria.ac.uk\/drbentley\/index.php?rest_route=\/wp\/v2\/posts\/2415\/revisions"}],"predecessor-version":[{"id":2437,"href":"https:\/\/research.northumbria.ac.uk\/drbentley\/index.php?rest_route=\/wp\/v2\/posts\/2415\/revisions\/2437"}],"wp:attachment":[{"href":"https:\/\/research.northumbria.ac.uk\/drbentley\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=2415"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/research.northumbria.ac.uk\/drbentley\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=2415"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/research.northumbria.ac.uk\/drbentley\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=2415"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}