A new slow transition process for energy transport in magnetically confined plasmas is reported. The slow transition is characterized by the change between two metastable transport conditions characterized by a weak and a strong electron temperature (T e ) dependence of normalized heat flux. These two branches are found to merge at the critical T e gradient. In metastable transport, the derivative of normalized heat flux to the T e gradient, @ Q e =n e =@ ÿrT e , is positive, while it becomes negative during the transition phase. The time for the transition increases as the normalized T e gradient is increased and exceeds the transport time scale characterized by the global energy confinement time. High-temperature plasmas, which are seen in nature or confined in laboratory experimental devices, are very often far from equilibrium. One of the noticeable features is that the observed temperature profiles are realized as a balance between the flux of energy and turbulence-driven transport. A sudden and distinctive jump in the profiles has been found, e.g., the H-mode phenomena in tokamak plasmas [1]. The concept of a profile transition has been employed in order to understand the formation and dynamics of a self-sustained radial profile in inhomogeneous plasmas [2]. The transitions in the profiles in H-mode plasmas are caused by the S-curve property (cusp-type catastrophe) in the gradient-flux relation, showing a feature analogous to the first-order phase transition [3,4]. As the second-order phase transition exists in terrestrial matter, it is possible that the other type of transition exists in the process of transport barrier formation in plasmas [4]. We here report the discovery of a new, slowly evolving transition between two transport branches that have different electron temperature (T e ) dependences, in toroidal plasmas.Heat transport in a magnetically confined toroidal plasma is strongly governed by turbulence. The mechanism of the H-mode transition has been studied both experimentally [5,6] and theoretically [7][8][9][10]. The bifurcation of the radial electric field and the associated suppression of turbulence are the key parameters to explain the fast transition between L mode and H mode [2,[7][8][9]. This transition mechanism appears as various phenomena such as the electric pulsation observed in the toroidal helical plasma and the internal transport barrier [11]. In parallel with these rapidly evolving phenomena, slow evolutions into the state of improved plasma confinement have also been observed in a slowly developing edge transport barrier [12] and in other confinement improved modes. Improved Ohmic confinement mode [13], counter-neutral beam injection improved mode [14], pellet enhanced performance mode [15], radiatively improved mode [16] associated with the change of density profiles without a clear fast transition. Although these improved modes associated with the slow transition are as important as the improved mode with a fast transition, the mechanism of the slow transition has been studie...