“…Important questions remain regarding: the relative role of large-scale versus mesoscale convection; the statistical distribution of bubble lengthscales and the physical mechanisms that govern the size of bubbles; and the mechanisms of ion energization. Sparse spacecraft coverage of the transition region makes it difficult to spatially localize injections, although there has been progress using multi-mission conjunctions (Turner et al, 2017) and remote imaging (Keesee et al, 2021). Modeling multiscale transport in the transition requires: energydependent drifts, not included in global magnetohydrodynamic (MHD) models which approximate the bulk flow with the E × B drift; fast flows outside the quasi-static slow-flow approximation of typical ring current models (see e.g., Toffoletto, 2020, and references therein); global models that produce mesoscale magnetotail structures (Wiltberger et al, 2015;Cramer et al, 2017;Merkin et al, 2019), requiring high spatial resolution and low dissipation algorithms; a representation of wave-particle interactions that contribute to ion heating (e.g., Chaston et al, 2014;Cheng et al, 2020); and finite gyroradius effects, as H+ in the transition region is known to be quasi-adiabatic (Runov et al, 2017), while heavier ion species like O+, likely exhibit highly nonadiabatic behavior (e.g., Moebius et al, 1987;Delcourt et al, 1997;Nosé et al, 2000;Keika et al, 2013;Bingham et al, 2020).…”