Local transport dynamics of cold pulses in tokamak plasmas

For over two decades, our fundamental understanding of energy transport dynamics in the core of tokamak plasmas had been challenged by the striking observation of temperature perturbation reversals following the injection of cold pulses at the plasma edge. These phenomena were first discovered by Gentle et al. (Phys. Rev. Lett. 74(18):3620–3623, 1995) in 1995 and had long been suggested to be evidence of nonlocal transport effects. In recent years, a new explanation to these phenomena has emerged, fully consistent with the theory of turbulent transport in magnetized plasmas and in remarkable agreement with experiment. This article reviews the experimental observation of temperature reversals in tokamak plasmas and presents the explanation based on local transport physics.


Introduction
Understanding transport of energy, particle and momentum in the core of tokamak plasmas is key to the development of fusion energy. Performance of tokamaks is highly dependent on the temperature and density gradients attained in the plasma core, thus revealing the need of developing reliable and accurate models of neoclassical and turbulent transport to design new fusion devices and plan operational scenarios. Burning plasmas will soon be achieved in machines such as SPARC ) and ITER (Shimada et al. 2007), and increasing efforts in the area of integrated modeling are being devoted to "predictfirst" (Staebler et al. 2022) scenarios that are both safe and efficient enough to achieve performance goals.
Key to integrated modeling of the plasma core is the understanding of turbulent transport dynamics, known to be the dominant contributor to energy losses in modern tokamaks. Turbulent energy transport has come a long way since it was first identified as the main source of anomalous transport, dramatically increasing transport levels from collisional, neoclassical theory predictions. The gyrokinetic theory of plasma turbulence is now widely considered to provide the appropriate theoretical framework to simulate turbulence in magnetically confined plasmas (Brizard and Hahm 2007;Garbet et al. 2010), and has successfully explained many experimental observations of tokamak plasmas thanks to extensive validation exercises (White 2019, and references therein).
Nonetheless, there are still some fundamental phenomena that remain unexplained. Until recently, one prominent example was the seemingly "nonlocal" temperature response of the plasma core to edge perturbations, originally described by Gentle et al. (1995a) and that had led to dozens of journal publications for over two decades (Rodriguez-Fernandez 2019, and references therein). An example of such experimental observation is depicted in Fig. 1a for a lowdensity Alcator C-Mod plasma, in contrast to the standard core temperature response at high density (Fig. 1b). The observation of a fast (faster than energy confinement times) increase in the temperature of the core of the plasma caused by the cooling of the edge region (i.e. during a "cold-pulse" experiment) had long remained a mystery to be explained by modern turbulent transport theory. The expected locality of turbulent transport, owning to the small spatial scales of turbulent structures compared to the device size in modern tokamaks (radial correlation lengths of only a few ion gyroradii Fonck' et al. 1993;McKee et al. 2001;Rhodes et al. 2002;Conway 2008;Cabrera et al. 2021), could not allow such reversal of the temperature perturbation, and therefore nonlocal effects seemed to be required to explain the full dynamics. In cold-pulse experiments, the magnitude of the core temperature rise is significant and with associated time scales shorter than an energy confinement time. Were nonlocal effects needed to explain it, modern turbulent transport theory and models would need to be revisited from the ground up, thus casting doubt on the accuracy of the predictions of burning plasmas in SPARC Howard et al. 2021) and ITER ) with theory-based transport models.
Recent experimental and computational work resulting from an international and multi-institution collaboration between the Alcator C-Mod (Rodriguez-Fernandez et al. 2018a, 2019a, DIII-D (Rodriguez-Fernandez et al. 2019b) and ASDEX Upgrade (Angioni et al. 2019) teams has provided an explanation to the cold-pulse mystery without the need of nonlocal transport effects. Purely local transport models used nowadays to study and predict tokamak behavior are demonstrated to be enough to explain these effects, as long as multi-channel transport interactions are accurately captured within the model assumptions.
This paper provides a review of cold-pulse experiments and the local explanation for the temperature reversal effect. Section 2 provides a background of coldpulse experiments and associated phenomenology. In Sect. 3, the local paradigm is presented, including past work, simulation results and quality of agreement with experiment. In Sect. 4, discussion and thoughts on the topic are provided. Lastly, in Sect. 5 we present a few final remarks.

Background
Impurity injections at the edge of magnetically confined plasmas have long been used for impurity transport studies (Friichtenicht 1974;Marmar et al. 1975;Howard Fig. 1 a Example of temperature reversal effect in Alcator C-Mod in a low-density plasma, as measured by different channels of ECE diagnostic. b Example of standard cold-pulse propagation at high density et al. 2011; . When trace amounts of impurity neutral atoms are deposited at the edge of the plasma, they get ionized and transported inwards, reaching the central part of the plasma in a fast time scale. Once these impurity ions reach the plasma core, the time evolution of the concentration of impurity charge states as they decay back to background levels can be used to characterize impurity transport via inference frameworks ) and forward particle transport models, such as STRAHL (Dux 2004) and Aurora .
When larger amounts of impurities are introduced, the local increase of electron energy losses via ionization and radiation at the plasma edge can drive a sharp drop in temperature (Kissick et al. 1994). This "cold pulse" in temperature can propagate inwards, thus serving as a useful tool to survey local heat transport dynamics throughout the entire plasma core. Characteristic elements of energy transport such as electron heat stiffness and pinch radial profiles can be inferred by analyzing the time traces of the plasma response to temperature perturbations (Gao et al. 2014).
Not long after the first utilization of perturbative impurity injections to examine electron heat transport dynamics, the phenomenon of "core temperature reversal" was identified in the TEXT tokamak (Gentle et al. 1995a, b), as shown in Fig. 2.

Evidence as a transport phenomenon
Soon after their discovery, certain explanations for the temperature reversal effect were discarded based on experimental observations and basic intuition on spatial and time scales of plasma phenomena. Increased Ohmic heating in the plasma core as a consequence of current redistribution could not be possible, as well as changes in the local q-profile. Experiments in TEXT (Gentle et al. 1995a, b) showed no evidence of an increase in central current density or changes in Joule heating. Not only the current redistribution time scales are much longer than the observed temperature effect, but also such redistribution could cause the onset of sawteeth or significant variations in the plasma internal inductance, which were not observed experimentally. Furthermore, the later evidence of such effects in stellarators means that the phenomenology cannot be mediated by the current density (Gentle et al. 1997). Changes in the effective ion charge required to account for the change in resistivity needed for the core temperature heating would need to be much larger than those observed and inferred (through loop voltage measurements) experimentally, and certainly the amount of neutral impurities injected at the edge of the plasma would need to be much larger (Gentle et al. 1995a, b). The lack of heating power redistribution as a possible cause of these effects was also studied in Tore Supra (Zou et al. 2000), which again showed no evidence of such in discharges heated by lower hybrid waves.
Changes in plasma position and magnetic flux surfaces could compromise diagnostic measurements, and therefore they were studied early on (Gentle et al. 1995a, b;Kissick et al. 1996). Measured changes in plasma position and flux surfaces were negligible for the cold-pulse experiments studied. The effect of suprathermal electrons compromising temperature measurements from electron cyclotron emission 1 3 10 Page 6 of 29 (ECE) diagnostics was also discarded after analysis of Tore Supra experiments (Zou et al. 2000). Presence of MHD modes (already existing or triggered by the cold pulse) could affect transport properties and were also considered as candidates. Experiments in Tore Supra (Zou et al. 2000) and LHD (Tamura et al. 2007) showed no significant changes in the amplitude of low-m MHD modes during the temperature reversal phenomenon. Work in Alcator C-Mod (Rice et al. 2013) indicated that no tearing modes were seen during the cold-pulse experiments, likely discarding the effect of magnetic islands on transport as a mechanism for the temperature reversal.
In conclusion, all experimental evidence gathered over the years seemed to suggest that the temperature reversal had to be a transport effect through fast changes to core electron heat diffusivity, e . The question of which parameter (and where in the plasma) caused a significant change to electron heat transport properties with the right time scales was the big conundrum that had largely remained unanswered for over two decades.

Effect of density and plasma current
One of the most agreed upon observations of the phenomenology of the temperature reversal phenomenon is its strong dependency on density. Either through a smooth decrease of the temperature reversal magnitude, or by a sharper transition, all machines have observed the disappearance of the temperature reversal effect when increasing the average plasma density. Oftentimes, the linear scaling of the temperature reversal with density was better characterized when a T −1∕2 e factor was included. Such scaling was required to include discharges at different plasma current and heating schemes (Kissick et al. 1996;Ryter et al. 2000). Experimental studies in Alcator C-Mod (Rice et al. 2013;Gao et al. 2014) proposed a normalized collisionality dependence to explain the density and plasma current scalings. The ratio of the electron-ion collision frequency and the bounce frequency of trapped electrons simplified (Gao et al. 2014) as: was able to bring together observations from different tokamaks. Here, n e is the lineaveraged density, q 95 is the safety factor at the 95% poloidal flux surface, and R 0 is the device major radius. In Eq. 1, it was assumed that the factor Z eff ∕T 2 e is roughly constant, and that the aspect ratio is similar for the tokamaks included in the study. Figure 3 collects observations from different tokamaks (Gao et al. 2014), revealing that the transition to standard temperature drop depends on n e q 95 R 0 .
Importantly, experiments in RTP Hogeweij et al. 2000), ASDEX Upgrade (Ryter et al 2000) and Alcator C-Mod (Rice et al. 2013) demonstrated that the radial location of temperature profile "flex point" was affected by the plasma current. In particular, the flex point moved outwards (higher radii) with increased plasma current, suggesting some connection with the structure of the q-profile.

Effect of external heating and collisional equilibration
While temperature reversals are mostly characteristic of low-collisionality Ohmic plasmas, they have been observed with different levels of auxiliary heating. TEXT experiments (Gentle et al. 1997) showed that the temperature reversals got weaker with the addition of electron cyclotron heating (ECH) power, and experiments with neutral beam injections (NBI) in TFTR (Kissick et al. 1998) and JET (Mantica et al. 2002) showed no temperature reversals in such ion-heating dominated plasmas. Experiments in Tore Supra (Zou et al. 2000) observed, however, that temperature reversals were more present in discharges heated by lower hybrid (LH) waves, and studies in ASDEX Upgrade (Ryter et al. 2000;, RTP (Hogeweij et al. 2000) and HL-2A (Hong-Juan et al. 2010;Sun et al. 2011) showed that ECH had actually a positive effect in increasing the magnitude (and range of existence) of temperature reversals. It has long been suggested (Callen and Kissick 1997) that the temperature reversal phenomenon is closely related to the electron-ion power transfer. Consistent with the observation of the collisionality dependence from Eq. 1, past studies (Callen and Kissick 1997;Mantica et al. 2002), have noted that the effects take place when the ratio of collisional electron-ion power transfer to electron heat conduction is sufficiently small (i.e. thermally decoupled electrons and ions). Experiments in ASDEX Upgrade (Angioni et al. 2019) with Ohmic, ECH-and NBI-heated plasmas at different densities found that the core temperature perturbation was strongly dependent on the electron to ion heat flux ratio ( Q e ∕Q i ) at mid-radius.

Connection with other transport phenomena
The appearance of these effects in low-collisionality, non ion-heated plasmas suggested somehow a relationship with linear Ohmic confinement (LOC) characteristics, and the disappearance of the reversals at high density seemed to be correlated with the transition to saturated Ohmic confinement (SOC) (Rice et al. 2013. In fact, the intrinsic rotation reversal phenomenon (Camenen and Angioni 2017;Rice et al. 2011) has been hypothesized to play a role in the transport dynamics that could modify the core electron temperature response. Experiments with Ohmic plasmas in Alcator C-Mod (Rice et al. 2013;Gao et al. 2014) and with ECH plasmas in KSTAR (Shi et al. 2017) showed that the disappearance of the temperature reversal was concomitant to the reversal of intrinsic toroidal rotation.
More recent work in Alcator C-Mod (Rodriguez-Fernandez et al. 2017) demonstrated, using a large database study, that temperature reversals occur regardless of the direction of the intrinsic toroidal rotation (including in locked-mode, zero-rotation plasmas), which challenged the idea that these energy transport effects could be connected somehow to momentum transport phenomena. Figure 4 depicts a database study of core temperature perturbations and intrinsic rotation reversals. At high current ( I p = 1.1MA in Alcator C-Mod), the intrinsic rotation reversal had no effect on the behavior of the core temperature perturbations: temperature reversals occurred with both co-and counter-current toroidal rotation.

Nonlocal modeling of cold-pulse propagation
In parallel to the study of cold-pulse propagation with local models, which will be the focus of the remaining of this review, several theories based on nonlocal models have been suggested, and the reader is referred to a companion publication for details on this topic (Ida 2022).
In the context of nonlocal transport, the experimental observation of the fast temperature reversal phenomenon suggests that prompt changes in thermal diffusivity leading to the increase in core temperature are caused by changes in the outer region (Ida et al. 2015). It is proposed that mechanisms such as mesoand macro-scale fluctuations, energetic particles effects and turbulence spreading mechanisms could cause the connection between interior and outer plasma regions, even if they are separated by many micro-turbulence scale lengths. From these, the turbulence spreading model has been prominent in the literature of nonlocal explanations to the cold-pulse phenomenon (Naulin et al. 2015;Hariri et al. 2016) and has been compared qualitatively to experimental observations Liu et al. 2019;Shi et al. 2020). By this mechanism, coupling of distant regions can be possible due to fast spreading in the turbulence intensity field, which in principle can be independent from local plasma parameters and gradients. Another notable model discussed in the literature is the fractional diffusion model (Del-Castillo-Negrete et al. 2008), which utilizes transport operators that can make the heat flux at a given point depend on nonlocal contributions from distant regions or on global properties of the temperature profile.
While we recognize that such models can provide the basis for the fast propagation and reversal of cold pulses in tokamak plasmas, in this review paper we focus our attention on the local explanation to the cold-pulse phenomenology, as it has shown agreement with many experimental observations of the effect. And more importantly, it was able to reconcile the widely used and validated local drift-wave turbulent transport theory with the fast cold-pulse propagation, without the need of additional or different transport mechanisms. The reader is invited to read more details about nonlocal transport explations in Ida (2022).

The local paradigm
As described in Sect. 2, temperature reversals have broadly been considered as a consequence of transient changes in background energy transport. Most models to explain the temperature inversions were phenomenological (Gentle et al. 1995a, b;Kissick et al. 1996;Gentle et al. 1997;Hogeweij et al. 2000;Gao et al. 2014). In this context, a phenomenological model consists of ad-hoc modifications to the local diffusivity (or, in some cases, the local electron heat pinch (Gao et al. 2014)) so as to cause the temperature reversal. However, such models generally shift the problem from "What causes the fast temperature reversal?" to "What causes the fast electron heat transport reduction?" In other words, the changes in transport coefficients that would be required to produce a temperature reversal of the right magnitude and with the right time scale may be inconsistent with the underlying physics that dictates neoclassical and turbulent transport properties in a self-consistent manner.

First studies of local transport effects
Before the full explanation that was recently provided by the Alcator C-Mod other studies that pioneered the use of local transport models to explain the temperature reversal effect. However, the lack of accurate physics models to simulate turbulent transport dynamics at the time prevented the achievement of a full explanation for the experimental observations. Furthermore, in these first studies, multi-channel interactions (those related to the effect of the transient density gradient flattening on electron heat transport) were not properly accounted for. These were found to be key in later work, as it will be discussed in following sections. Kinsey et al. (1998) first proposed the use of local models (in particular, stiff critical gradient models) to simulate computationally the core temperature response to edge cold pulses. Simulations with the IFS-PPPL model (Kotschenreuther and Dorland 1995) reproduced a qualitative reversal of core temperature, mostly as a consequence of the effect of T i ∕T e on ion temperature gradient (ITG) stability that followed the cold-pulse injection. Quantitative match of steady-state kinetic profiles and the magnitude and time scale of the cold pulse propagation and reversal magnitude still remained largely elusive. Similar modeling of cold pulses in ASDEX Upgrade by Ryter et al. (2000) also yielded qualitative agreement, but the edge perturbation required for a moderate reversal was much larger than in experiment. Finally, simulation work by Mantica et al. (2002) of JET discharges also showed the possibility of temperature reversals with local physics, although the ion transport dynamics had to be modified for the effect in the electron channel to appear.

Improved quasilinear modeling
Building from previous work, simulations with local models within flux-matching (i.e. integrated modeling) frameworks were revisited after two decades of developments in quasilinear transport theory and computational techniques. The TGLF model (Staebler et al. 2005(Staebler et al. , 2007 solves an electromagnetic, linear, gyro-fluid set of equations for multiple plasma species to find the unstable drift-wave-type modes in the core of tokamak plasmas. It is suited to capture trapped ion and electron modes (TIM, TEM), ion and electron temperature gradient modes (ITG, ETG) and kinetic ballooning modes (KBM). One of the cornerstones of TGLF is its ability to capture trapped particle dynamics more accurately compared to the original gyro-Landau-fluid (GLF23) transport model (Waltz et al. 1997). This made TGLF a great candidate to explore cold-pulse dynamics. Given the low collisionality that characterizes experiments that exhibit temperature reversals, it is expected that TEM-like turbulence is one of the main contributors to electron heat transport in the plasmas of interest. In TGLF, a saturation rule is used to find the saturated level of potential fluctuations that may result from the linearly unstable modes. Saturation rules are constructed using databases of fully nonlinear gyrokinetic simulations, and one of the latest saturation rule (SAT-1) (Staebler et al. 2016(Staebler et al. , 2017 was demonstrated to successfully model the physics of the nonlinear upshift of the critical gradient (Dimits shift) (Dimits et al. 1996) and multiscale coupling effects . Validation studies of TGLF SAT-1 (e.g. Rodriguez-Fernandez et al. 2018b) have shown improved accuracy with respect to previous gyro-fluid solvers and previous saturation rules.
Integrated transport simulations with the TRANSP (Breslau et al. 2018) and ASTRA (Pereverzev and Yushmanov 2002) power balance frameworks were used to study the cold pulse behavior. TGLF SAT-1 was used to calculate the self-consistent time evolution of the turbulent electron and ion heat diffusivities following the cold pulse injections. In steady-state, conditions that exhibited temperature reversals were dominated by TEM modes at long wavelengths, which emphasized the need to include accurate trapped particle physics in the turbulence simulations to explain cold-pulse dynamics. Figure 5 depicts characteristic turbulence spectra for low and high density conditions that exhibit temperature reversals and diffusive drops respectively. At low density, the presence of TEMs causes the electron heat flux to be sensitive to the variation of density gradients. Linear ionscale gyrokinetic simulations confirmed the same results as the gyro-fluid modeling (Gao et al. 2014): at low density, the plasma core sits in a TEM-dominated regime, sensitive to density gradients, while the high-density plasma becomes more sensitive to ion temperature gradients.

The effect of density gradients
To introduce cold pulses at the edge of otherwise stationary plasmas, actuators that inject particles into the plasma are used. The most typical actuator used for such experiments is the laser blow-off (LBO) technique, but others such as transient gas puffs (e.g. Macwan et al. 2021), pellets (e.g.  or super-sonic molecular beam injections (SMBI) (e.g. Sun et al. 2011) have been used. The injection of neutral particles come with an associated increase in radiation and ionization losses, as well as a possible decrease in plasma average temperature as a consequence of equilibration of the plasma with the "colder" impurity ions (i.e. via an isobaric On the left, ion-scale a linear growth rate and b real frequency spectra for characteristic low and high density Ohmic plasmas in Alcator C-Mod, as calculated with the gyro-fluid TGLF model. Modes with < 0.01c s ∕a have been removed for plotting purposes. On the right, scans of driving gradients for the c low-and d high-density cases. Heat fluxes are calculated using SAT-1. In low-density Ohmic plasmas, the main microinstability in the plasma core is TEM, which depends strongly on electron temperature and density gradients. At high density, ITG becomes dominant at low-k ( k s < 0.6 ), and density gradients do not have a strong effect on electron heat transport process). LBO injections have the advantage of the localization of the impurity source at the edge, if a proper impurity species is used. One main consequence of the injection of cold pulses in tokamak plasmas is the unavoidable perturbation of background density profiles. Such perturbations have been reported in many devices. For example, pellet injections in RTP ) that displayed temperature reversals had a large associated increase in electron density (reaching up to 40% the background density), as inferred from the inversion of interferometer measurements. SMBI experiments in J-TEXT  showed changes in density of about 10% , and LBO injections in Alcator C-Mod (Gao et al. 2014) showcased changes on the order of 30% in average electron density. Recent work with gas puffs in ADITYA-U (Macwan et al. 2021) observed increases in electron density of about 15-20% , and reported a flattening of the core density profile. Because the change in electron density does not happen in a self-similar way, density gradient scale lengths can widely vary throughout the propagation of the cold pulse from the edge to the core. Because it is recognized that drift-wave-type turbulent transport is affected by density gradients (especially in low-collisionality TEM-dominated plasmas), the observation of changes in density profiles during the cold pulse propagation was a clear indication that multi-channel transport effects were required if local physics was to be successful in explaining the intriguing temperature reversals.

Heat transport simulations
Simulations with only energy transport prediction required the evolution of a density perturbation for the temperature reversal effect to appear in simulations. In Alcator C-Mod (Rodriguez-Fernandez et al. 2018a), fast interferometer data indicated the possible propagation of a density pulse that followed the cold pulse injection at the edge, but the spatial profile of the perturbation was unclear at the time. Dedicated experiments in DIII-D (Rodriguez-Fernandez et al. 2019a) demonstrated, thanks to high time and spatial resolution density profile reflectometry, that a density pulse originated at the edge and propagated inwards, reaching the core in a very fast time scale. The evolution of the electron density pulse was consistent with the evolution of the impurity profile as given by the STRAHL code (Dux 2004) and with experimentally inferred diffusion and convection impurity transport coefficients. ASDEX Upgrade experiments also measured a significant perturbation to the electron density, determined via the integrated data analysis suite (IDA) (Fischer et al. 2010).
When a density pulse propagates from edge to core in a monotonically increasing density profile (as it is often the case in the core of tokamaks), a reduction of density gradients at the front of the density pulse is always present. A reduction of density gradients in plasmas with density-gradient driven turbulence comes with an associated reduction of transport levels, and can be understood as a stabilization (or a transient transport barrier) effect. This is precisely what was first observed in simulation work of Alcator C-Mod cold-pulse experiments (Rodriguez-Fernandez et al. 2018a), when, for the first time, the impact of the density evolution was taken into account in the modeling of the temperature response.
The TGLF model successfully captured the stabilization of linear, mid-k TEM modes in the plasma core, as a manually imposed density pulse propagated from edge to core. From a quasilinear point of view, a reduction of the linear growth rate will cause a decrease of potential fluctuations and thus lower transport levels. This dynamics can be observed in Fig. 6, which depicts the core plasma behavior following the edge perturbations. At low density (left part of the figure), the density pulse causes the reduction of linear growth rate of electron modes, thus reducing heat flux, allowing the temperature to increase (evoking what can be understood as a transient heat transport barrier).
This phenomenology was observed in the simulations, and was used to predict-first experiments in the DIII-D tokamak (Rodriguez-Fernandez et al. 2019a). Near quantitative match of the magnitude of the temperature reversal in DIII-D (a machine that had never reported temperature reversals in the past) was achieved, evidencing the superb predictive power of the TGLF model. In the DIII-D study, it was determined that not only the reduction of the electron density gradient caused the reduction of electron heat transport levels, but also the increase in impurity content and reduction of impurity density gradient also contributed to TEM stabilization. Increased collisionality and effective ion charge also contributed, but their effect was not very strong given the small variation of these parameters during the cold-pulse propagation.
High time resolution density profile reflectometry in DIII-D (Zeng et al. 2006) was used to confirm the existence of a density pulse that propagates from edge to core, as depicted in Fig. 7. When such a measured pulse was fitted to a skewed Gaussian pulse and introduced in simulation, excellent match of the time traces of the temperature perturbation were obtained. This was further confirmation that the density pulse propagation from edge to core, as measured, was capable of producing temperature reversals with the right magnitude and time scale, and conclusively confirmed that a local transport model of heat transport was enough to reproduce the temperature reversal phenomenon. In fact, this work reconciled the steady-state prediction of kinetic profiles based on standard heat diffusion equations, with the much faster propagation of the cold pulse from edge to core.

Particle transport dynamics
Experimental and simulation work in Alcator C-Mod and DIII-D had proven that a local transport model of heat transport could explain the temperature reversal. Nonetheless, it still remained unclear whether the propagation of the measured density pulse could be explained by local physics as well. Simulation work of ASDEX Upgrade LBO-injected cold-pulse experiments  with the ASTRA transport code (Pereverzev and Yushmanov 2002) coupled to STRAHL (Dux 2004) for impurity transport calculations was used precisely to answer this question.
The carbon impurity ions present in the plasma following the LBO injection were included as kinetic species for the turbulent transport calculations with TGLF SAT-1, and both deuteron and carbon particle fluxes were used to self-consistently evolve the density profiles. Electron density was calculated by quasineutrality. Electron and ion temperatures, and deuterons and impurity densities were predicted all the way from the separatrix to the magnetic axis, and the arrival of carbon from the LBO was modeled as a local increase in carbon separatrix density. Simulations showed the exceptionally fast propagation of a density pulse from edge to core, quickly flattening the density gradients at mid-radius. Electron heat conductivity is also observed to be reduced as a consequence of the flattened density gradients in the TEM dominated plasmas, thus leading to the increase of the electron temperature peaking. Figure 8 shows the dynamics during these coupled simulations. The fast propagation of the density pulse was explained by the destabilization of an impurity density gradient driven mode as a consequence of the extremely large reversed impurity gradients that form locally after the LBO injection.
The work in ASDEX Upgrade was convincing proof that local transport physics not only explained the heat transport channel, but also the particle transport channel dynamics was fully explained by predictive transport simulations with the TGLF SAT-1 model.

Match with experimental phenomenology
As described in Sect. 2, a number of experiments have showed consistent trends with plasma parameters such as average density and plasma current. Increase of average density has universally led to a reduction of temperature reversal magnitude until their disappearance at high enough densities. Extensive simulation work via numerical experiments of Alcator C-Mod plasmas (Rodriguez-Fernandez et al. 2019b) has addressed this. Density scans at constant plasma current show the clear decrease of temperature reversal magnitude at higher plasma density (Fig. 6). The increased collisionality led to de-trapping of electrons and the consequent reduction of TEM growth rates compared to ITG at long wavelengths. When the relative strength of estimated mixing length transport of TEM becomes lower than ∼ 0.7× of ITG, temperature reversals disappeared. This is concomitant to the equalization of ion-electron collisional exchange and the conducted power density, which is consistent with the observation of stronger temperature reversals in plasmas with decoupled ions and electrons, widely reported in the literature of cold-pulse experiments, and as discussed in Sect. 2.3. The loss of dominance of TEM in relation to ITG meant that the turbulence stabilization effect of the density pulse propagation is weak and not strong enough to compensate the increase of temperature gradients (destabilization of temperature gradient modes) at the pulse front, which causes the ubiquitous drop of temperature in the high density cases. Modeling of ASDEX Upgrade plasmas (Angioni et al. 2019) with prescribed experimental electron density evolution also confirmed the absence of electron temperature reversals in conditions of intermediate to high density, as ITG turbulence becomes more dominant in the ion-scale range.
The effect of plasma current was also explored in numerical experiments of Alcator C-Mod plasmas (Rodriguez-Fernandez et al. 2019b). At fixed density, the increase in plasma current led to increase electron heating, which further destabilized background TEM turbulence. Presence of more TEM turbulence throughout the plasma led to increased sensitivity to the density pulse propagation, fully consistent with past work that showed that very high densities were required to make the temperature reversals disappear at high plasma current (Rodriguez-Fernandez et al. 2017) (Fig. 4). Simulation results of the effect of increasing plasma current are depicted in Fig. 9. Furthermore, the electron temperature flex point also moved outwards with increasing plasma current, following the motion of rational surfaces observed in previous studies (Rice et al. 2013;. Simulation work of ASDEX Upgrade (Angioni et al. 2019) cold pulses also explored the predictive capabilities of TGLF on cold-pulse propagation in low-density NBI-heated plasmas. In the case studied, a reduction of ion-scale electron thermal diffusivity was observed following the density gradient flattening, but the destabilization of electron temperature gradient modes at short wavelengths prevented the appearance of the temperature reversal, consistent with experimental observations (see Sect. 2.3). It was determined that conditions with comparable values of T e and T i are affected by ETG transport and therefore cannot benefit from the reduction of density gradients that result from the density pulse propagation.

Discussion and open questions
The local paradigm presented in Sect. 3 has been extraordinarily successful in explaining cold-pulse phenomenology, while at the same time being truthful to the widely known and validated physics of drift-wave-type turbulent transport in the core of tokamak plasmas, and capable of reproducing the fast propagation of heat and particle transport events along with the slower dynamics which takes place on confinement time scales. The ubiquity of density perturbations in cold-pulse experiments suggests that the dynamics of particle transport and associated time variation of gradients are the key to the understanding of temperature reversals and cold-pulse propagation. This is often referred to as multi-channel transport effects, and the understanding of their importance is critical for the development of reliable models for particle, momentum and heat transport in tokamak plasmas. The literature of temperature reversal experiments is very extensive, and several other phenomenological observations still remain to be attempted from the perspective of local transport modeling. Nevertheless, there is no evidence that local transport physics as described in Rodriguez-Fernandez et al. (2018a, 2019a; Angioni et al. (2019) is unable to explain the rest of experimental observations.

On the magnitude of the density perturbation
Only LBO-induced cold pulses were studied, as the localization of the particle sources facilitates the interpretation of the experimental results and the subsequent modeling. However, temperature reversals from cold pulses introduced via gas puffs, pellets and SMBI could be equally explained by the dynamics of the density profile and associated perturbations. As indicated in Sect. 3.3, density pulses are not exclusive from laser blow-off injections and therefore there is no basis as to assume that local transport physics does not explain the temperature reversals in such experiments.
Remarkably small density perturbations have been demonstrated (Rodriguez-Fernandez et al. 2019a;Rodriguez-Fernandez 2019), in both experiment and simulation, to be enough to cause the reversal of the pulse in temperature, as shown in Fig. 10. Significant changes to the density gradient are possible even with small changes in the absolute value of density (we hypothesize this as the reason why density pulses have long been disregarded in previous cold-pulse modeling studies). Transient changes in gradient scale lengths in situations where heat transport is stiff can cause the appearance of temperature reversals of significant magnitude and with fast time scales. Figure 11 shows that a 10 times smaller density perturbation can produce temperature increases of significant magnitude in Alcator C-Mod simulations.

Ion temperature evolution
Cold-pulse experiments throughout the years have focused on the electron temperature perturbations, and the ion temperature evolution has only been reported in a few cases, owning to the difficulty to measure ion temperature with a time resolution enough to resolve the transient perturbation. Measurements of the ion temperature profile using crystal spectroscopy in Alcator C-Mod (Rice et al. 2013;Gao et al. 2014) showed an increase in central ion temperature when reversals in the electron temperature occurred, although such evolution typically happened on a longer time scale. In the experiments reported, the edge ion temperature dropped promptly, and a moderate increase above pre-injection levels occurred later on. At high densities, similar to the electron temperature behavior, the increase in central ion temperature disappeared. Recent work in J-TEXT (Shi et al. 2020) further explored the ion temperature dynamics, and confirmed the Alcator C-Mod experimental results, with remarkable resemblance. The dynamics of the ion temperature during the density scan of Fig. 6 is in qualitative agreement with the experimental observations, as depicted in Fig. 12. Drops in edge ion temperature (followed by a later increase) are observed at low density, as well as an increase in core ion temperature. As the background density is increased, the core reversals also disappear.
However, the dynamics of the ion temperature during cold-pulse experiments is not as clear and robust as with the electron temperature. Experiments in Alcator C-Mod were reported where the edge ion temperature drop is not present (Rodriguez-Fernandez et al. 2018a). Instead, an increase (on a longer time scale) of ion temperature was observed throughout the entire plasma. Furthermore, it has also been observed in Alcator C-Mod that the rise of central ion temperature may appear even when there is no corresponding electron temperature reversal (Cao et al. 2016).
Quantitative and qualitative match of the diverse ion temperature responses during cold-pulse experiments still remains an open question. However, given the Fig. 12 Dynamics of the core and edge electron and ion temperatures during the density scan of Fig. 6 experimental evidence provided and the modeling results, it can be determined that electron and ion temperature dynamics are not strongly coupled. The evolution of the ion temperature played a negligible role in driving the central electron temperature response in the simulations described here (as evidenced, for example, by the negligible effect of a∕L T i on electron heat flux in the low-density TEM-dominated plasma in Fig. 5c), and the fact that experimentally the ion temperature "reversal" is observed regardless of the behavior of the electron temperature channel suggests that the two are de-coupled. The full ion temperature dynamics and the interplay with the turbulence evolution that causes the electron temperature response still requires further investigation. Nonetheless, as indicated at the beginning of Sect. 4, there is no evidence that it could not be explained by local transport physics.

Effect of ECH deposition location
Experiments in RTP (Hogeweij et al. 2000) showed that the ECH deposition location has one of the strongest effects on the magnitude of the temperature reversals, reaching values of up to 80% the background temperature when ECH is applied nearby the q = 1 flux surface. While this trend remains unexplained at the time of writing this review article, it could still be reconciled within the local transport physics explanation. As ECH heats up the approximate region of the temperature flex point, stronger TEM turbulence activity is expected at low collisionality. This means that electron heat transport may become more sensitive to density gradients, thus driving a stronger temperature reversal for the same density pulse perturbation. This is similar to the observation that at higher plasma current stronger temperature reversal responses at the same average density, as discussed in Sect. 3.6. The link to the location of rational surfaces could be due to the changes in what regions of the plasma are dominated by one or other type of turbulent modes, as discussed in Rodriguez-Fernandez et al. (2019b).
Nonetheless, these past experiments should be reconsidered in detail in light of the local modeling results presented in this review article. Modeling should be performed in the specific conditions of those experiments before one can conclude that the local paradigm can or cannot explain also those observations.

Final remarks
This review article has covered the experimental and simulation evidence that coldpulse propagation in the core of tokamak plasmas can be fully explained by local physics, including the seemingly "nonlocal" temperature reversal and speed of propagation. This experimental observation had been puzzling plasma physicists for over two decades, but thanks to experimental and computational work in Alcator C-Mod, DIII-D and ASDEX Upgrade, there is convincing evidence that the widelyvalidated theory of drift-wave turbulent transport can indeed fully explain the experimental behavior in such perturbative experiments while, at the same time, capture the steady-state confinement features. In fact, local or flux-tube transport models are sufficient to cause the fast reversal of core electron temperature. The variation of density gradients during the perturbative transport response has been demonstrated to both cause the temperature reversal at low collisionality (Rodriguez-Fernandez et al. 2018a) and also be in alignment with the local paradigm .
While this work provides further evidence that local transport modeling within the framework of the conventional quasilinear turbulent transport description and the usual diffusive heat transport equation has been successful in explaining the propagation of fast LBO-induced pulses, this study has not covered other perturbative transport phenomena (such as the propagation of heat pulses using modulation of ECH heating power), explored the validity of standard local transport theory in stellarators, or claimed perfect and conclusive understanding of all transport problems. Nonetheless, the coldpulse propagation in the core of tokamaks seems to be fully dominated by local turbulent transport effects. Truly nonlocal mechanisms were not needed to reconcile experiment and simulation, and while they may be important to study certain transport phenomena, they certainly were not the dominant effect in the cold-pulse experiments in tokamaks.
As it has been shown in this work, the fact that the plasma is able to react on time scales that are extremely fast (much faster than the energy confinement time or time scales derived from average heat diffusivities) cannot be considered evidence that nonlocal effects have to be present. This work has demonstrated that plasma dynamics can produce very fast responses (which at first sight may seem nonlocal) that are still fully compatible with conventional transport physics, as described by secondorder diffusive equations, and simultaneously also predict slower dynamics, on time scales which are comparable to experimentally-relevant confinement times. The compatibility of the conventional local transport paradigm with very fast phenomena is allowed by the rich nonlinearity present in the parametric dependencies of the transport coefficients. Such dependencies directly connect different plasma parameters, leading to interesting multichannel transport effects that must be accounted for in the predictions of transport in tokamak plasmas.