Sanae-Inoue Itoh 1952–2019: a memorial note for a pioneer researcher of plasma bifurcation

This memorial note for Professor Sanae-I Itoh presents her specific achievements in physics research alongside her wider record of accomplishment in the field of magnetically confined plasmas. The topics include bifurcation phenomena (e.g., H-mode and improved confinement modes), turbulence-generated structures (e.g., zonal flows and streamers), and fundamental concepts and processes in plasma turbulence (e.g. nonlinear couplings and energy transfer. The note focuses initially on results obtained through her integration of theory, simulation, and experiment, particularly those arising from a low temperature plasma facility at Kyushu University. We then describe contemporary challenges in plasma turbulence which Sanae addressed with great interest, and consider some of the perspectives that were opened by her achievements.


Introduction
During the past half century, there has been substantial progress in obtaining and understanding the high-performance plasmas that are necessary for nuclear fusion (Gibney 2022; Clery 2021Clery , 2022;;Tollefson and Gibney 2022).The conditions 21 Page 2 of 27 required for a burning plasma have been largely achieved, and the next generation of experimental facilities will exploit this progress systematically.Plasma turbulence is among the key physics challenges on this journey, raising questions and intriguing phenomena that are also important elsewhere in contemporary physics: for example, structural bifurcation, non-local property of transport, turbulence generated structures.Prof. Sanae-I.Itoh was one of the great plasma physicists who challenged these mysteries in the research field of magnetic confinement plasma.Many of Sanae's results were obtained in collaboration with her husband and longterm research partner Prof Kimitaka Itoh [whose obituary of Sanae was published in PPCF (Itoh 2020)], and with an extensive network of colleagues internationally and within Japan.Amongst her physics achievements, the most famous one is the first proposal of a theory that reminds us of the significant importance of radial electric field, E r and its nature of bifurcation (Itoh and Itoh 1988; Ida et al. 1990; Fujisawa  1997) after being inspired by the discovery of H-mode transition by Prof. F. Wagner and colleagues (Wagner et al. 1982).Moreover, she pushed forward her research, always indicating new directions and advanced concepts in understandings of turbulent plasmas (Itoh and Itoh 1996), i.e., self-sustained nature of plasma turbulence, turbulence suppression related to electric field structure (Itoh et al. 1994a), abrupt events often accompanied with collapses (Itoh et al. 1991(Itoh et al. , 1998)), physics of zonal flows and streamers (Diamond et al. 2005;Fujisawa et al. 2004Fujisawa et al. , 2007;;Inagaki et al. 2011), and developing tools for plasma turbulence research (Nagashima et al. 2005;Itoh et al. 2005).Furthermore, she promoted her own integrated project, called Itoh project, in which theory has initiative on simulation and experiment.Particularly, she constructed her own experimental facility at Kyushu University, having noted the excellent potential of low temperature plasma for turbulence research, to produce various fundamental results (Yamada et al. 2008;Inagaki et al. 2016).This memorial note, focusing largely on the work of Sanae Itoh, provides a brief description of her above-mentioned achievements and associated work alongside her research progress.Section 2 introduces theories for H-mode and improved confinement, as well as their experimental confirmations, while Sect. 3 focuses on the physics achievements internationally in her own project (or Itoh project), with ones obtained in her constructed facility.Concluding remarks are placed in Sect.4, where we discuss the unsolved problems of contemporary physics, relevant to plasma turbulence, in which Sanae took great interest.

H-mode and improved confinement
The H mode was discovered in ASDEX tokamak for the first time in 1982 (Wagner et al. 1982).The significance of H mode was to demonstrate that toroidal plasma can bifurcate, and that the plasma confinement can be improved above a heating power threshold.The latter gave a strong impact to fusion research at the time when the problem was a phenomenon, so-called power degradation, as empirical scaling laws show the dependence i.e., E ∝ 1∕ √ P , where E and P are the energy confinement time and auxiliary heating power, respectively (Kaye and Goldston 1985;Sudo et al. 1990).The discovery of H-mode stimulated searching other favorable operational regimes (Wagner and Stroth 1993;Burrell 1997;Fujita 2002), and really found various kinds of improved confinement modes (ICMs) classified three categories; a family of edge transport barrier (ETB) like VH-mode in DIII-D (Jackson et al. 1991), ones with internal transport barriers (ITB) (Ida and Fujita 2018;Hugon 1992;Koide et al. 1994;Levinton et al. 1995;Strait et al. 1995), and ones without a transport barrier, such as supershot (Strachan 1987) and radiative improved (RI) mode (Messian 1996).On the other hand, similar kinds of ICMs were found in stellarators (Fujisawa 2003;Wagner et al. 2006), ETBs (Erckmann et al. 1993;Toi et al. 1993).ITBs (Fujisawa 1999;Stroth 2001;Alejaldre 2001;Ida 2003) and ICMs without a barrier (Morita et al. 1993;Ida 1996Ida , 1999;;Brakel 1997).Nowadays, utilizing these ICMs, reactor-relevant operational scenarios have been developed to pursue the compatibility of steady and high performance in various plasma confinement systems; i.e., current hole, in JT-60 (Fujita et al. 2001), I-mode in tokamaks (Whyte et al. 2010;Stroth 2022;Fenstermacher 2022), HDH-mode in W7-AS (McCormick 2002), super dense core (Ohyabu et al. 2006) and ITB with an impurity hole in LHD (Ida et al. 2009).The clarification of ICM formation mechanisms, particularly H-mode, was one of the research subjects that Sanae pursued throughout her research career with Kimitaka.

Theory of H-mode
Bifurcation theory Sanae and Kimitaka proposed the first theory to describe LHtransition based on the bifurcation of the radial electric field, E r (Itoh and Itoh  1988).These theories predicted that changes in the radial electric field could be a cause of the H-mode transition.An equation to describe E r -generation is introduced to explain the theories.The equation is given as (Itoh and Itoh 1996) where ⟂ is the dielectric constant of plasma representing the effect of the polarization current (Shaing et al. 1992).The radial current on the right-hand side is composed of a number of elements: j Neo , j loss , j vis , j pvis , j trs and j flc are the currents driven by neoclassical bipolarity (Kovrizhnykh 1984;Hastings 1984), orbit loss (Itoh et al. 1989(Itoh et al. , 1991;;Sanuki et al. 1993), viscosity (Rozhanski and Tendler 1992), parallel viscosity, turbulence Reynolds stress (Diamond and Kim 1991) and fluctuations (Solomon and Shats 2001), respectively.Some of these current elements show nonlinear characteristics to cause E r -bifurcation.In their theory, the two terms are considered to be dominants on the right-hand-side, i.e., anomalous and ion-orbit loss term.If one solves the equation j flc (E r ) + j loss (E r ) = 0 , three solutions, as is shown in Fig. 1, can be found to correspond to bifurcated branches of an unstable, L and H modes.The theory was followed by an analogous one proposed by Prof. K. C. Shaing independently (Shaing and Crume 1989).Although the model was based on poloidal flow dynamics, radial electric field is coupled with poloidal flow through radial force balance, thus, the poloidal flows can be regarded as being equivalent (1) to the radial electric field.In case of the Shaing's theory, the combination of the ion-orbit loss and the parallel viscosity term gives similar solutions.While these theories revealed stimulating and contributed to the development of other models, such as one based on Stringer's spin-up (Hassam et al. 1991) and the predator-prey model (Diamond et al. 1994), radial electric field became a long-standing topic in the research of magnetic confinement fusion.
Self-sustained turbulence The next challenge to explain the H-mode formation was why the turbulence transport was reduced after the transition of E r -structure.Essen- tial theories are needed to explain experimental observation of turbulence and transport (Liewer 1985;Wootton et al. 1990) and their reductions at a barrier including H-mode (Connor and Wilson 2000).Turbulence should be saturated to satisfy the balance between its growth and decorrelation (or damping) rate, as dI∕dt = L I − DΔI, where I is the turbulence intensity and the second term represents the damping due to diffusion whose constant is expressed by D. Setting Δ = k 2 ⟂ in a steady state, the above equation gives the well-known Kadomtsev's expression of the diffusion coefficient, . Under the presence of E r -shear (or velocity shear), the turbulence structure is deformed along the direction of velocity, which is characterized by a change of effective perpendicular wavenumber, , assuming that the perpendicular velocity varies as obeying V ⟂ = E×B x = (E r ∕B) � x , where x is the perpendicular coordinate (Itoh et al. 1994a).Then, the turbulence diffusion coefficient is obtained, as . The consideration is a guidance for similar expressions that most of theories expect for turbulence reduction due to E r -shear (Itoh et al. 1994a; Shaing et al. 1988;  Biglari et al. 1990).Similarly, thermal coefficient should obey the formula like , where and E are the thermal diffusivity of the L-mode and a numerical coefficient, respectively.At present this phenomenon is well known as E r -shearing effect that should regulate or suppress turbulence.
Sanae and Kimitaka advanced their thoughts on nonlinear growth of instabilities, and gave a birth to a concept of self-sustained turbulence.A mode can be destabilized through nonlinear interactions among fluctuations even if the mode is linearly stable.Thus, the linear growth rate L needs to be replaced with nonlinear growth rate NL in the above expression of diffusivity.Moreover, the nonlinear growth rate of a mode could be balanced with that of nonlinear damping in a stationary state.This leads to the concept that turbulence is self-sustained above a threshold in fluctuation amplitude (Itoh et al. 1992a, b, c), and that such a mode can grow abruptly with a property of subcritical onset.Actually, examples have been found that feature such nonlinear growth of instabilities in simulations (Scott 1990(Scott , 1992;;Carreras et al. 1992;Dimits et al. 2000).The concept was first applied on interchange mode turbulence (Itoh et al. 1994a;Yagi et al. 1994Yagi et al. , 1995)).Later, the theory was extended to the ballooning mode turbulence including both L and H modes (Itoh et al. 1994b).For the case of ballooning mode turbulence, the expression of thermal diffusivity is obtained as follows (Itoh et al. 1994b;Yagi et al. 1993;Fukuyama et al. 1994): where (= −q 2 R � ) , h 1 , h 2 and f are numerical coefficients, and E1 and E2 repre- sent the variables related to electric field shear and electric field curvature, with p , v A , c, q, s and R being the plasma frequency, the Alfven velocity, the light veloc- ity, the safety factor, the magnetic shear parameter and the major radius of plasma, respectively.The formula demonstrates the importance of E r -curvature in addition to E r -shear, as is shown in Sect. 2.4.
Predator-Prey model The nonlinear coupling of fluctuations can generate plasma flows, or radial electric field in toroidal systems, through a phenomenon known as inverse cascade (Hasegawa 1985).In other words, radial currents are induced through the turbulent Reynolds stress and result in the formation of the radial electric field.The turbulence Reynolds stress to generate poloidal flows has been paid attention as an important player in barrier formation, supported by experimental observation (Hidalgo et al. 1999).Taking into account the feedback effect of the generated E r -shear on tur- bulence (Diamond and Kim 1991;Hidalgo et al. 1999), the interplay between fluctuations and radial electric field is recognized as a mechanism for barrier formation and structural bifurcation in toroidal plasmas (Diamond et al. 1994).Based on the abovementioned concept, a model termed predator-prey model was proposed to explain the H mode transition.In the model, the interplay between turbulence and its generated flows is expressed as the following simultaneous equations, where F(= (ñ∕n) 2 ) and U = (V � E ) 2 are the normalized turbulence and the shear flow intensities, respectively.The steady state solutions of the above equations, i.e., dF∕dt = dU∕dt = 0 , correspond to L-and H-mode; that is, U = 0 , F = 0 ∕ 1 (L-mode) and U = ( 0 − 1 ∕ 3 )∕ 2 , F = ∕ 3 (H-mode).In addition to two sta- ble solutions corresponding to the bifurcated states, the equations also predict limit cycle oscillations as is the Lotka-Volterra equations.Experimental support was obtained (Xu et al. 2000), and more recent simulation with realistic configurations shows that the H-mode transition should occur as synergy of the Reynolds stress and ion-orbit loss effects (Chang 2017).

Experimental confirmation
Electric field bifurcation In response to the bifurcation model, the structural change of radial electric field before and after LH-transition was measured with spectroscopic measurements in DIII-D and JFT-2M tokamaks (Ida et al. 1990;Groebner et al. 1990).Figure 2 shows the measurement of E r -profiles in JFT-2M tokamak.The force balance equation, E r = v ϕ B θ − v θ B ϕ + ( p i ∕ r)∕Zen i , allows to evaluate radial electric field from flow velocities and impurity temperature measured with spectroscopy (Fonck 1984), where v , v , Z, n i , and p∕ r , repre- sent poloidal, toroidal velocities, ion charge, ion density, and ion pressure gradient, respectively.The experiments were performed, soon after the proposal of E r -bifurcation model, to confirm if the radial electric field really controlled plasma confinement or created transport barrier.In continuous current tokamak (CCT), ( 4) TUMAN 3 and TEXTOR tokamak, the hypothesis was experimentally confirmed using a bias electrode inserted to the plasma edge.As a result of employing the biasing electrode, plasmas were successfully forced to make a transition to a better confinement mode (Taylor et al. 1989;Askinazi et al. 1992;Weynants et al. 1992).The E r -structure generated with the bias electrode was discussed theoreti- cally (Rozhansky and Tendler 1992;Stringer 1993;Itoh et al. 1998).The result demonstrates that the change in E r -structure should be the cause of confinement improvement, although there may be a possibility that the mechanism can be different from the case of spontaneous transition.

Causality problem
The force balance equation raised a causality problem since it showed negative pressure gradient should result in negative radial electric field.Thus, one can consider that the change in radial electric field after H-mode transition can be a result of the barrier formation.The time scale of electric field change was the key to answer the question.The time scale of spontaneous E r -transition was first observed with HIBP in CHS plasma (Fujisawa 1997;Fujisawa et al. 2000) to show that the observed time scale, a few dozen microseconds, was much faster than the confinement time of the plasma, and the change of radial electric field after the transition is positive direction, which may lessen the pressure gradient in the force balance equation.Later in a sawtooth-triggered H-mode transition in JFT-2M, it was also shown that the E r -transition could occur in a few dozens of microseconds, much faster, compared to confinement time scale, in the JFT-2M tokamak (Miura et al. 2001;Ido et al. 2000;Ido 2001).These facts strongly support the hypothesis that the E r -shear should trigger the barrier forma- tion.The barrier formation in CHS is found to completely obey the E r -bifurcation scenario that Sanae and Kimitaka predicted (Itoh and Itoh 1988;Fujisawa et al. 2000;Fujisawa 2021).
Turbulence reduction The conditions for turbulence reduction have been tested using the expression of Eq. ( 2) at a barrier.The formula shows that the fluctuation reduction is effective when L ∼ ExB .Simulations and experiments demonstrated the validity of the condition,  L <  ExB , for E r -shear reduction (Waltz et al. 1994;  Hahm and Burrell 1995; Waltz et al. 1997); for reviews (Burrell 1997a, b).Moreover, intensive efforts have been made to observe E r -structure and turbulence.Pro- gress in diagnostic techniques for high temperature plasmas (Bretz 1997), such as HIBP (Jobes and Hickok 1970;Schoch et al. 1988;Crowley 1994;Fujisawa et al. 1996), spectroscopy using a diagnostic beam (Fonck 1984;Levinton et al. 1989;Fonck et al. 1993), phase contrast imaging (PCI) (Coda and Porkolab 1995), and reflectometry (Manso 1993;Mazzucato and Nazikian 1993;Mazzucato 1998;Hirsch et al. 2001), has made it possible to measure both fluctuations and radial electric field when transport reduction or transport barrier formation occurs.Experiments using advanced diagnostics have been carried out in DIII-D and W7-AS, reporting that the reduction of local fluctuations is accompanied with E r -shear for- mation around the ETBs (Doyle et al. 1991;Moyer 1995;Kislov and T-10 Team 2001;Burrel 2001;Wagner 1994).In the enhanced reversed shear (ERS) mode of the TFTR, reflectometry showed drastic suppression of fluctuations of rather long-wavelength inside the barrier (Mazzucato and Nazikian 1993).Moreover, a micro-wave-scattering system found that short wavelength fluctuations of electron 21 Page 8 of 27 skin-depth ( c∕ pe ∼ 1 mm) was suppressed, and that its temporal behavior was well correlated with the electron transport (Wong et al. 1997).Furthermore, development of diagnostics, such as gas-puff imaging technique, revealed successfully how turbulent eddies should be broken into small pieces in the bias-electrode experiments (Shesterikov et al. 2013).

Contemporary issues for improved confinement
Edge localized modes and intermediate states The physics of edge localized modes (ELMs) is a critical issue at present, because their control is extremely important for fusion application in two aspects; to remove confined impurities and ashes from the plasma, and to avoid damage to the surrounding wall due to the energy pulses associated with ELMs.In the early phase of their discoveries, Sanae and Kimitaka regarded the ELMs as repetitive transitions between bifurcates states.According to the concept, they first proposed a model of dithering ELMs, which occurs near the power threshold to H-mode transition, as the repetitive transitions between L-and H-mode (Itoh et al. 1991).Then the model was extended to describe type-I ELMs as repetitive transitions between M-, L-and H-mode; M-mode was introduced as the bifurcated state characterized by braided magnetic field in high-plasma (Itoh et al. 1996).The experiment confirmed that the model should well describe dithering ELMs (Zohm 1994).On the other hand, self-excited or limit cycle oscillations (LCO) were observed in stellarators (Stroth 2001;Fujisawa 1998;Maassberg et al. 2000), Particularly, CHS experiments reported the presence of many bifurcated states in potential with various transition patterns between them, called electric pulsation (Fujisawa et al. 2000).Later, in tokamaks, the finding of an intermediate regime (IM-mode) between L and H mode, characterized by periodic bursts (Colchin et al. 2002), extended the predator and prey model to explain the phenomena (Kim and Diamond 2003).Moreover, later in 2010s, the further findings of I-phase, LCOs between two states dominated by GAM and turbulence (Conway et al. 2011), simulated exploring transient phase between L-and H-mode.Many experiments followed to report similar transient behavior or intermediate states before the full H-mode transition (Estrada et al. 2010;Xu et al. 2011;Schmitz et al. 2012;Estrada et al. 2011;Cheng et al. 2013;Kobayashi et al. 2013;Yan et al. 2014), providing active research field of H-mode formation mechanisms and scenarios involving physics of zonal flows (Cziegler et al. 2014;Czieler et al. 2015;Schmitz 2017).
Shock structure The importance of two-dimensional flow structure or poloidal asymmetry in electric field for H-mode transition was pointed out (Shaing et al. 1992;Kasuya and Itoh 2005).As is shown in Fig. 3a, where a shock structure with poloidal asymmetry is expected to be induced during H-mode transition.The asymmetric potential on magnetic field surfaces produces inward convective flows that give a large impact on particle transport and a rapid rise of density inside.A mysterious phenomenon to suggest shock formation was really observed during electric pulsation in CHS, in which repetitive transitions occurs to form steady state oscillation.During the electric pulsation, abrupt density rises at the plasma edge were found to occur simultaneously with drops in potential at the plasma center.The time scales of the rises and drops are in a dozen microseconds, which is much faster than a diffusive one of millisecond (Fujisawa et al. 2000).Two-dimensional observation should be needed to investigate transient heat and particle fluxes in such phenomena as well as rapid electric field change.Particularly, such two dimensional measurement should be of significant importance also for considering dangerously large heat (and particle) fluxes accompanied with ELM crashes.
Curvature effects of radial electric field While the magnetic field structure is one of the most important ingredients for improved confinement (Shafer et al. 2004;Waltz et al. 2006), more detailed investigation showed that the strongest reduction of turbulence could occur at the position of minimum radial electric field where E r -shear is rather small (Kamiya et al. 2017).The observation indi- cated that only E r -shear should not be sufficient to explain turbulence reduction in barrier regions.Then, an important role of E r -curvature was revisited (Itoh et al.  1992a, b, c) to be another ingredient for barrier formation theoretically (Itoh and Itoh 2016).Sanae and Kimitaka introduced a term for a saturation mechanism, which represented nonlinear energy transfer from turbulence to larger scale structure like zonal flows through modulational couplings.Then the equation for turbulence saturation is written as where E r -structure parameter, Z = (E � r ) 2 − E r E �� , with E ′′ representing the E r -cur- vature.In equilibrium, the above equation gives a modified formula of Eq. ( 2), D = D 0 ∕(1 + Z) .Then the relation between transport and the parameter, Z, has been investigated and the hypothesis was confirmed experimentally.Figure 3b, c (5) show examples in JT-60U and LHD that gradient at barrier positions should linearly increase to be saturated with the parameter, Z, and that temporal behavior of turbulence amplitude should be correlated with the E r -curvature.The observations dem- onstrate significant impact of E r -curvature on turbulence reduction (Kamiya et al.  2016; Tokuzawa et al. 2014; Tokuzawa 2017).

Establishment of research center for plasma turbulence
It was 2004 that Sanae started her own project for implementing her integrated research of plasma turbulence in her original manner that unified theory, simulation and experiment.The project was supported by two large grants of Japanese Society for the Promotion of Science (JSPS) (KAKENHI in Japanese) that she acquired sequentially.Thanks to the project, she promoted domestic and international collaborations in addition to pushing her own works forward.Besides, Sanae had an opportunity to construct her own low temperature facility, and succeeded in obtaining a number of significant results on fundamental processes of turbulent plasmas.This section introduces the achievements mainly on zonal flow physics obtained in experiments strongly supported in Itoh project: CHS, JFT-2M and her own facility.

Physics of zonal flows
Zonal flow identification Sanae and Kimitaka were among the pioneers of the study of zonal flows in plasma physics (Itoh et al. 2006).In 1996, a simulation result based on gyrofluid model sensationally called into question whether ITER could reach ignition (Glanz 1996).A theoretical answer was that turbulence driven flows should not be strongly damped in collisionless processes (Landau damping), and that the undamped flows, called zonal flows, should play a significant role in turbulence saturation (Rosenbluth and Hinton 1998).The situation stimulated simulations and experiments on zonal flows (Fujisawa 2009;Tynan et al. 2009;Garbet 2010).
Simulation works demonstrated that short-wavelength turbulence could generate zonal flows (Hasagawa and Wakatani 1987;Lin et al. 1998) through nonlinear couplings, and that the zonal flows could regulate the short-wavelength turbulence again to control transport (Hahm et al. 1999).Then, many experimental trials were carried out to search for zonal flows (Fujisawa et al. 2004;Xu et al. 2003;Gupta et al. 2005), while Sanae published the well-known review for the physics of zonal flows with Prof. P. H. Diamond and colleagues (Diamond et al. 2005).
The zonal flows are fluctuating mesoscale structures, which are driven by microscale instabilities such as drift waves.The zonal flows are symmetric around the magnetic axis, and have no poloidal wave number with a finite radial one in mesoscale.During the period, the identification of zonal flows was completed with twin HIBPs in CHS plasma (Fujisawa et al. 2004).Figure 4 shows the spectrum of the measured electric field fluctuations.In the measurements, it was confirmed that the low frequency fluctuations less than ∼ 1 kHz showed a long-distance correlation.This is one of the characteristics for zonal flows to satisfy.The correlation analysis deduced a spatiotemporal structure of the low-frequency electric field fluctuations; note that Sanae called this method 'Correlation Hunting'.The result shows that it should have a sinusoidal structure of a finite radial wavelength of ∼ 1.5 cm in mes- oscale, and that the structure should decay in timescale of ∼ 1.5 ms.In the same experiments, zonal magnetic fields (Diamond et al. 2005;Gruzinov et al. 2002;Guzdar et al. 2001) were found to coexist with zonal flows with extended use of the HIBP (Fujisawa et al. 2007).The result demonstrates that turbulence has an ability to generate magnetic fields in mesoscale.
Geodesic acoustic modes Geodesic acoustic modes (GAM) form one of the physics issues that attracted the interests of Sanae and Kimitaka.The GAMs are oscillatory branches of zonal flows that can only emerge in toroidal plasmas.The oscillatory property makes it easier to detect GAMs than stationary zonal flows, thus, a huge number of experimental works have been reported, as is explicitly shown in a recent exhaustive review (Conway et al. 2022), both in tokamaks and stellarators (Nagashima et al. 2005;Shats 2002;Jakubowski et al. 2009;Nagashima 2006;Zhao . 2006;Schoch 2003;Hamada 2005;Ido 2006;Melnikov 2006;Mckee 2003;Krämer-Flecken et al. 2006;Fujisawa et al. 2007).These inter-experimental comparisons confirmed the essential characteristics of GAM frequencies, f GAM ∝ c s ∕R , where c s is sound speed, the spatial symmetry ( m = n = 0 ) in potential fluctua- tions, and m = 1 nature in density fluctuations (Zhao et al. 2006).On the other hand, experimental observations often reported the eigenmode characteristics of GAM; the frequency appeared constant over a certain radial extent of the identical GAM profile that exhibits discrete and intermittent changes, regardless of radial variation of electron temperature to give dependence of local GAM frequency.After paying attentions on the eigenmode characteristics (Itoh et al. 2006), Sanae and Kimitaka devised a novel method, called GAM spectroscopy, that determines the effective ion mass through the ion sound velocity (Itoh et al. 2007) in addition to the exact position of plasma boundary (Nagashima 2009;Itoh 2016).The GAM spectroscopy was further extended to infer the safety factor under co-existence of ion acoustic mode (Itoh et al. 2009).The possibility of ion mass detection for evaluating fuel ratio in fusion plasmas, and the other essential plasma parameters could be beneficial in future burning plasmas.In addition to the above works, they led many of experimental, theoretical and simulation works on GAMs, including their basic excitation mechanisms (Itoh et al. 2005;Sasaki 2009), turbulence enhancement and suppression due to its trapping of GAM structure (Sasaki 2017;Sasaki et al. 2018), experimental confirmation of the turbulent Reynolds stress drive of GAMs (Kobayashi 2018), and the mechanism of energetic particle excited GAMs (E-GAMs) that should work as transfer channels of momentum and energy to bulk ions (Sasaki et al. 2011;Sasaki 2017), Zonal flow effects on turbulence and transport Zonal flows affect turbulence and confinement through effects of their turbulence shearing, their quasi-energy conservation with turbulence (or energy partition), and turbulence trapping (Kaw 2002) that can contribute to turbulence spreading (Mattor and Diamond 1994;Garbet et al. 1994;Hahm et al. 2004;Kobayashi 2022).Zonal flow shearing was confirmed in several experiments of smaller plasma experiments in the early phase of the research, for example, H1 heliac (Shats et al. 2003) and TJ-K (Birkenmeier et al. 2013;Ullman et al. 2020).The effects of energy partition was demonstrated in transport barrier formation of CHS.After transport barrier is formed, the increase in electron temperature inside the barrier was observed as a puzzling fact, since the inside E r -shear was not large enough to reduce the turbulence.Detailed investigation of energy partition showed that the energy of zonal flows should increase after the barrier formation, hence the turbulence was reduced consequently to improve confinement without E r -shear inside.Note that the cause of the increase of zonal flow fraction is ascribed to the reduction of neoclassical viscosity (Fujisawa et al. 2008).The result clearly demonstrates the importance of the energy partition for confinement (Itoh 2007;Toda et al. 2007).The direct interaction between zonal flows and turbulence, probably turbulence trapping, was also observed in CHS (van Milligen et al. 1995;Fujisawa 2006).As is shown in Fig. 5a, increase (or decrease) in turbulence amplitude is synchronized with decrease (or increase) in zonal flow's amplitude.Moreover, the nonlinear couplings between elemental waves in turbulence should be stronger at minima of zonal flows (Fujisawa et al. 2007a, b).A simulation work also showed turbulence trapping of GAMs. Figure 5b shows the simulation results to demonstrate that GAMs could trap turbulence in their valley and convey the turbulence along their trajectories of propagation and could result in turbulence spreading that enhanced transport in a turbulence-free region (Sasaki 2017).
Anatomy of turbulence and findings Sanae and Kimitaka were ambitious in developing advanced methods to analyze turbulence, or visualize turbulence interaction, which they call anatomy of turbulence (Itoh et al. 2005(Itoh et al. , 2017)).Following early works (Moyer et al. 2001;Tynan et al. 2001), the nonlinear  (Sasaki et al. 2018).It is clearly shown that the trapped turbulence is transferred inward Fig. 6 Analyses of nonlinear couplings.a Bicoherence analysis on frequency domain to prove that GAM is generated by nonlinear couplings with background turbulence in JFT-2M (Nagashima et al. 2005).b Energy transfer analysis on wavenumber domain to prove that stationary zonal flows are generated by nonlinear couplings with background turbulence in TJ-K (Manz et al. 2009b).c The top and bottom panels show a reconstructed image of the giant fluctuating structure, and the result of bicoherence analysis to prove that a giant fluctuating structure is generated with background turbulence in LHD, respectively (Inagaki et al. 2011) couplings between GAMs and background turbulence was confirmed, as is shown in Fig. 6a, using bicoherence analysis on a reciprocating Langmuir probe data in ohmically heated plasma of JFT-2M (Nagashima et al. 2005).The analysis also succeeded in evaluating coherent biphase, i.e,, phase relationship between GAM and turbulence modulation (Nagashima 2006).Before the bicoherence analysis was widely provided, the power transfer function (PTF) analysis was invented to elucidate not only the presence of nonlinear couplings but energy transfer direction between different scale fluctuations or structures in plasma turbulence (Ritz et al. 1989;Kim et al. 1996Kim et al. , 1997;;Xia and Shats 2003;Nagashima et al. 2008).The analysis method, based on Hasegawa-Mima equation (Hasegawa and Mima 1977) evaluates the power transfer function T(k 1 , k 2 ) using the following formula, where k and P k represent linear growth rate and the power of a mode.After the first PFT analysis on frequency domain (Ritz et al. 1989;Xu et al. 2010), it was applied on wavenumber domain using two dimensional (2D) probe data (Manz et al. 2008(Manz et al. , 2009a)).In TJ-K torsatron, the energy transfer was evaluated the two-dimensional wavenumber space, and finally as is shown in Fig. 6b, succeeded in proving the nonlocal energy transfer from drift-waves to zonal flows (Manz et al. 2009b).On the other hand, the integrated analysis of correlation hunting and visualizing nonlinear processes in Itoh project allowed a discovery of a giant fluctuating structure in LHD (Inagaki et al. 2011).The structure fluctuates over the plasma cross-section, as is shown in Fig. 6c with a frequency in a few kHz range.The structure is considered to be a nonlinearly excited metastable mode, probably ones related to dissipative trapped ion modes (DTIM) (Kadomtsev and Pogutse 1971), by the background turbulence; the nonlinear couplings were confirmed using the cross-bicoherence analysis in Fig. 6c.The experimental result demonstrates that turbulence can generate a macroscopic structure as well as mesoscale ones like zonal flows.

Construction of low temperature facility and its physics achievements
Promotion of cylindrical plasmas Low temperature plasma experiments, whether linear or toroidal in geometry, can provide a suitable arena for studying the physics of plasma turbulence.The advantageous points of such experimental facilities are as follows; (i) high accessibility, superior flexibility, quick turn-around of experiments, (ii) excellent controllability of plasma, (iii) easy acquisition of a huge amount of data for statistics, and (iv) suitability to use traditional probes and to develop new advanced diagnostics.In Itoh project, a linear plasma system named large magnetic device (LMD) (Shinohara et al. 2009) was upgraded to LMD-U, and later was fully modified to a plasma system, named Plasma assembly for nonlinear turbulence analysis (PANTA), for plasma turbulence research.The physics target of the experiment is to investigate fundamental processes of plasma turbulence mainly using multi-channel probe arrays (Oldenbürger et al. 2012;Yamada et al. 2007;Nagashima et al. 2011).The system has produced versatile results in plasma physics as well as contributions to advanced diagnostics and turbulence analyzing techniques (Kobayashi et al. 2010(Kobayashi et al. , 2011(Kobayashi et al. , 2012(Kobayashi et al. , 2021)).The physics achievements include identification of fundamental processes and properties of turbulence (Nagashima et al. 2009(Nagashima et al. , 2011;;Ohyama et al. 2012), findings of solitary wave oscillations (Arakawa et al. 2009(Arakawa et al. , 2010(Arakawa et al. , 2011) ) and streamers (Yamada et al. 2008(Yamada et al. , 2010;;Kin et al. 2018;Yamada et al. 2018;Kin et al. 2019), visuallization of interaction between nonlinear structures (Arakawas et al. 2016), remote energy transfer between local turbulence through zonal flows (Nagashima et al. 2011), difference in local and global-averaged transport fluxes (Nagashima et al. 2011), while development of fundamental tools does bispectral analyses (Kobayashi et al. 2010(Kobayashi et al. , 2011(Kobayashi et al. , 2012) ) and tomography (Fujisawa et al. 2015(Fujisawa et al. , 2016;;Yamasaki et al. 2020;Moon et al. 2021).Three selected topics amongst them are explained below.
Nonlocal energy transfer In the earlier phase of the project co-existence of zonal flows and drift waves was confirmed in a plasma parameter region of LMD-U.In the plasma, energy transfer between zonal flows and drift waves was evaluated using the cross-bispectrum described as ⟨E r (f ZF ) * E r (f 1 )E (f 2 )⟩ , where ⟨⟩ represents statistical average.The turbulence Reynolds stress term (or so-called Reynolds force)  r ũr ũ to generate zonal flows is expressed as the above cross-bispectrum in the Fourier domain.As is shown in Fig. 7a, the result indicates that the zonal flows should gain the drift wave energy in the plasma interior where a steep density gradient exists, while in the outer region the energy transfer is opposite (Nagashima et al. 2009).The result suggests that the spatial energy transfer should occur through the zonal flows.On the other hand, a similar analysis applied to beam emission spectroscopy (BES) data in DIII-D revealed a role of zonal flows; GAM convection mediated the energy transfer from low ( < 40 kHz) to high ( > 150 kHz) frequency fluctuations (Holland et al. 2007).
Identification of streamers Existence of streamers were identified in LMD-U and PANTA (Yamada et al. 2008;Kin et al. 2018Kin et al. , 2019)).They were searched in torus plasma due to their importance on transport, however a few experiments Fig. 7 Example of fundamental experiments using multichannel probe arrays for plasma turbulence.a Energy transfer analysis based on cross-bicoherence that shows non-local energy exchange between zonal flows and local turbulence in LMD-U (Nagashima et al. 2009).b Signals of streamers in ion saturation current measured with 64-channel azimuthal probes (Yamada et al. 2008).The 64 channel probes are located at an outer radius of the plasma.The carrier wave (stripes in figure ) propagates in electron diamagnetic direction (positive direction in azimuthal coordinate), while the streamers (spots on stripes) rotate in the opposite (or ion diamagnetic) direction.c A reconstructed image of streamers obtained using conditional averaged technique (Kin et al. 2019).The image shows the structure that should be azimuthally localized and radially elongated reported only signature of the existence (Politzer 2000;Hamada et al. 2006).Some of the phenomena could be ascribed to so-called avalanches caused by a domino-like process as expected in the model of self-organized criticality (SOC) (Bak et al. 1987;Diamond and Hahm 1995;Carreras et al. 1996;Newman et al. 1996;Kosuga et al. 2013;Hahm and Diamond 2018).The streamers are shortlived localized structure, whose shapes are radially elongated and poloidally (or azimuthally) localized, They are expected to be formed from the nonlinear couplings of background turbulence as is similar to zonal flows.However, streamers, differently from zonal flows, enhance the transport (Diamond et al. 2001;Nozaki et al. 1979;Beyer et al. 2000).The anatomy using developed analysis tools confirmed that the streamers are generated through nonlinear couplings between the modes called carriers (drift-waves) and mediator, and that streamers are symmetric along the magnetic axis direction (Yamada et al. 2018).Figure 7b, c shows that their temporal behavior and reconstructed two-dimensional image of streamer (Yamada et al. 2008(Yamada et al. , 2010)).Moreover, later it was confirmed that they really induced ballistic transport to directly connect the plasma core to the edge in one order faster time scale than that of stationary turbulence, and that the envelope of streamers and mediator was found to have the feature of solitons; both become more sharpened as their amplitude increases, demonstrating the characteristics of solitons (Kin et al. 2018).Moreover, it was recently theoretical expected that the nonlinear evolution of streamers should lead to transient excitation of breaths (or Rogue waves in ocean) (Kosuga 2022).
Strongly correlated turbulence In a plasma parameter region of PANTA, coexistence of two kinds of instabilities, resistive drift wave and D'Angelo modes (or parallel velocity gradient (PVG) mode) (D'Angelo 1965;Catto 1973;Barbet 1999) was found: the PVG mode has been receiving a lot of attention recently owing to its capability of zonal flow generation due to potential vorticity mixing (Kosuga 2017), and three-dimensional helical flow generation (Kosuga et al. 2020).The two instabilities played a major role of sustaining the plasma structure (Inagaki et al. 2016).These two instabilities induce particle and momentum fluxes through their associated fluctuations (Kosuga et al. 2015;Sasaki et al. 2017).The D'Angelo mode drives the inward directing particle fluxes to sustain the density gradient to cause the drift mode, while the drift mode does momentum fluxes to sustain the velocity shear to cause the D'Angelo mode.Therefore, two different kinds of instabilities sustain each other, that is, one maintains the free energy source or structural inhomogeneity to cause the other.The present case of PANTA provides an example of spatial structure sustained through collaboration of different instabilities that are strongly correlated with each other.Note that a simulation work reproduced repetitive transition between L-and H-mode due to temporal competition between resistive drift wave and Kelvin-Helmholtz instablities (Numata 2007).This observation led Sanae and Kimitaka to acquire a new concept of strongly correlated turbulence, and to become aware of the advanced problems how structure and dynamics are selected in a nonequilibrium system where various kinds of inhomogeneity coexist and compete to drive different instabilities.

Unsolved questions of interest to Sanae
The preceding sections have briefly described Sanae's achievements in H-mode physics and the Itoh project.These are only a subset that can be introduced here among other works with so wide and multifaceted prospects: for instance, the series of statistical theories (Itoh and Itoh 1998, 1999, 2000a, b, c) and exhaustive monographs such as 'Transport and Structural Formation in Plasmas' (Itoh et al. 1999;Yoshizawa et al. 2003;Diamond et al. 2010).This section introduces some of remaining physics problems that Sanae challenged intensively.
Abruptness and symmetry-breaking in collapse phenomena Sanae published a review paper entitled 'Physics of events in toroidal plasmas' in 1998.She discussed of structure and dynamics in catastrophic phenomena, such as sawtooth oscillations, ELMs, and disruptions (Itoh et al. 1998).Sanae focused on their probabilistic occurrence and the explosive growth of catastrophic event, and notified that the key to clarify the physics of collapse should be the symmetrybreaking perturbation, which is distinguished from the precursors, to trigger the explosive growth of instabilities.In other words, abrupt structural change could be triggered by a local event to break symmetry of plasma, which could be subcritical nature of instabilities originating from nonlinear interactions between elemental waves and modes in turbulence (Itoh et al. 1998;Dimits et al. 2000).
Recently, several experiments have suggested the important roles of symmetrybreaking phenomena for abrupt changes of global structure.It has been recently found that the onsets of bursty phenomena should happen to begin at a local position; for example, ELMs (Lee et al. 2017), sawtooth instability (Yun et al. 2011), and tongue instabilities in LHD (Ida et al. 2016).The common feature of these phenomena is the occurrence of a local event to trigger the abrupt and global changes of structure.Moreover, the subcritical onset of GAMs was observed, which was caused by phase space distortion (Ido et al. 2016;Lesur et al. 2016;Itoh 2016;Wang et al. 2018).The concept of symmetry breaking should be key to understand the underlying physics of the abrupt changes in the entire plasma structure.
Nonlocal transport and its hysteresis Nonlocal transport is a long-standing problem.The famous examples should be the Gentle's cold pulse (Gentle 1995), and the Cordey's core-edge link on H-mode transition (Cordey 1995), while many others are available in both tokamak and stellarator experiments (Luce et al. 1992;Lopes Cardozo 1995;Kissick et al. 1996;Galli et al. 1998;Mantica et al. 1999;Galli et al. 1999;Milligen et al. 2002;Inagaki et al. 2004Inagaki et al. , 2006a, b;, b;Tamura et al. 2007).The mysterious points of these phenomena are that plasma transport at remote location from a position at which a perturbation is applied start to change in much shorter time scale than the confinement time or diffusive time.A possible candidate to explain the nonlocal transport or global linkage should be the nonlinear couplings between disparate scale of fluctuations.Communication between turbulence at two spatial points, even if located away more than the correlation length, is possible through the nonlinear couplings with long correlation structure, such as zonal flows, streamers and giant fluctuating structure.
On the other hand, it is commonly observed, for example, in ECH modulation experiment in LHD (Inagaki 2013), that the local transport can change immediately at the time a perturbation is applied without waiting for changes of local plasma parameters.Sanae and Kimitaka advocated a theory to explain the observations, with a simple phrase to express their theory 'Heating heats turbulence'.The phrase means that heating itself should cause a dynamics force through the distortion of phase space distribution (Itoh and Itoh 2012) to enhance additional turbulence, according to the formula (Itoh and Itoh 2013) where I 0 , 0 and h represent fluctuation intensity without heating, turbulence heat conductivity and growth rate parameter, respectively.This formula shows that fluctuations should increase (Itoh and Itoh 2012) if h ∼ 0 k 2 ⟂ ; the condition can be satisfied easier in smaller wave numbers.Nonlocal transport, which could be the result of the change in a macroscopic fluctuation structure, should occur with auxiliary heatings.Actually in an ECR-heating modulation experiment in LHD, the rise and drop of turbulence amplitude and resultant flux were really observed to induce transport flux changes, showing a transport hysteresis, simultaneously with increase and decrease in the heating power (Ida 2015).

Effects of hydrogen isotope and fueling
The hydrogen isotope effect is also a longstanding problem (Bessenrodt-Weberpals et al. 1993;Stroth 1998) giving underlying physics for the operational scenario and modeling for fusion performance strategy (Garcia et al. 2019).The essence of the hydrogen isotope effect is that deuterium plasma shows better confinement than that of hydrogen in contradiction to the prediction of gyro-Bohm dependence.Many experimental trials have been made to clarify the origin of the isotope effect both in tokamak (Xu et al. 2013) and stellarators, including LHD (Yamada et al. 2019;Kobayashi et al. 2019), in addition to theoretical works to investigate isotope effects (Hahm 2013;Itoh and Itoh 2016).
More recently, Sanae and Kimitaka proposed another theory in somewhat analogous concept of heating effects.They expressed the essence of the theory with a simple phrase 'Fuel feeds turbulence' (Itoh and Itoh 2012), suggesting the importance of atomic and molecular processes for confinement.The neutrals imprinted by the fluctuations in scrape-off-layer (SOL) penetrate into core plasmas, and enhance the core plasma fluctuations (Itoh et al. 2017).Then heavier ions should have less effect on the core fluctuations.They applied the concept on improved confinement associated with wall conditioning and gas puffing, such as supershot and improved ohmic confinement (Strachan 1987;Soldner 1988;Itoh et al. 1995).It should be mentioned that they derive the following formula, which indicates that density peaking should be restricted by neutral sources (Itoh and Itoh 2019) where K, L n , Δ n , p , C s , a represent the typical wave number of gradient-driven tur- bulence, gradient scale length, the penetration length of neutral particles, the global particle confinement time the ion sound velocity and the minor radius, respectively.The concepts need future experimental confirmation.
Advanced diagnostics and phase space structure and dynamics Sanae and Kimitaka provided several key concepts for further understandings of structural formation and dynamics of turbulent plasmas: bifurcation, global linkage, symmetry-breaking, abruptness, strongly correlated instabilities, effects of heating and fueling including atomic and molecular processes, etc.They also tried to promote experiments and simulations along these concepts.The research along these key concepts could need advanced experiments with new diagnostics that can cover the entire plasma region with locally fine resolution; it may need new plasma experimental system that enables such advanced measurements.A system of tomography, which has been used to detect MHD activities (Nagayama 1988;Park et al. 2006;Nagayama 1987), was developed to examine its feasibility of the measurement for the entire field of plasma turbulence a wide range from micro-scale ( ∼ ion Lamor radius) to macro- scale ( ∼ plasma radius) under a project of KAKENHI (specially promoted research) entitled 'Plasma Turbulence Observation System for Puzzling Out the Principles of Structural Formation and Functional Expression in Turbulent Plasmas' from 2017 to 2021.The feasibility for such turbulence measurement was confirmed (Fujisawa et al. 2015(Fujisawa et al. , 2016;;Yamasaki et al. 2020), and the system has been upgraded to measure three dimensional structure and dynamics of such turbulence field, and to reveal its unique features in PANTA (Moon et al. 2021).Synthetic diagnostics, which she started in Itoh project (Kasuya 2018;Sasaki 2018), have been developed to simulate actual turbulence signals presumed in measurements, and contributes to clarifying its own limitation of diagnostics by evaluating deformation of the measured signals from the real turbulence.
Finally, in contrast to neutral fluids, phase space structure and dynamics should play an important role in structural formation and dynamics of turbulent plasmas (Biglari 1988;Kosuga 2017), as vortex-like structure in phase space, called granulations, was suggested in 1970s (Dupree 1972;Dupress 1970;Kadomtsev and Pogutse 1970).Various kinds of phenomena related to phase space can be expected such as conversion of poloidal flows to toroidal flows (Kosuga 2013), frequency chirping instabilities (Berk 1999), explosive growth of instabilities (Lesur 2017), zonal flow generation (Kosuga and Diamond 2012), and the direct enhancement of turbulence due to auxiliary heatings (Itoh and Itoh 2012).Then the importance of the phase space measurement becomes crutial for further understanding of turbulent plasmas.At present, a new system of charge exchange recombination spectroscopy (CXRS) is being developed to aim at measuring phase space dynamics.Development of advanced diagnostics to access phase space information should be a challenging subject in future to clarify turbulent plasmas under a project of KAKENHI

Concluding remarks
This memorial note shows, by reviewing the achievements of Prof. Sanae-I.Itoh, that she made tremendous contributions to research in plasma physics and nuclear fusion, in experiment and simulation as well as in theory.At present, Wendelstein 7-X stellarator has started operation to aim at compatibility of steady-state operation and better confinement, in other words, to combine the advantages of the stellarator and tokamak confinement systems, in the sense of the symmetry of the magnetic field configurations (Klinger et al. 2019).Regarding the fusion performance, it was good news that the JET tokamak recently achieved doubling the record of fusion production energy (Gibney 2022), while recently the inertial fusion experiment in National Ignition Facility (NIF) demonstrated its great progress and promising future in achieving the fusion product Q = 0.67 (Clery 2021), and the latest one Q ∼ 1.5 (Clery 2022) , where Q is the ratio of fusion product to the input power.Moreover, using artificial intelligence (A.I.) to support the search for improved magnetic configurations has enabled excellent plasma performance in Tokamak à configuration variable (TCV) (Degrave et al. 2022).All these achievements tell us that really a new age has come to fusion research.However, a more precise and accurate prediction of plasma performance, based on physical first principles, is mandatory for efficient and early stage realization of a fusion reactor, particularly for avoiding uneconomical experimental costs for modern and future huge plasma confinement systems.Thus, research into the physics of plasma turbulence should be maintained to obtain solutions of many unsolved problems that Sanae continued to challenge.For obtaining sufficient understandings of plasma turbulence, it is essential to proceed the research under tight collaboration between theory, simulation and experiment.In order to improve the understanding of plasma turbulence to the level required, it is essential that research proceeds with close collaboration between between theory, simulation and experiment.In parallel with fusion-oriented research, specialized experiments in facilities oriented towards the fundamentals of plasma turbulence are called for.(Fujisawa 2021;Fujisawa and Conf 2018;Itoh et al. 2018).Finally, it is worthwhile to mention and her dreams should have been much more ambitious and so passionate; to mention that Sanae promoted an ambitious unified research plan across all fields of plasma physics and technology such as fusion, laser produced and processing plasmas, under the title of 'Research Network on Non-equilibrium and Extreme State Plasmas'.Sanae's phrase, 'Plasma is rich in bifurcation', expresses concisely the attractive nature of plasma as a subject of the modern physics, a system of ultimate states in non-equilibrium.

Fig. 1
Fig. 1 A theory for H-mode transition based on bifurcation of radial electric field.The theory assumes that the ion orbit loss current is balanced with anomalous electron current.a Ion and electron currents.The Gaussian-like curve and the lines represent the ion orbit loss and anomalous (or fluctuation-driven) current, respectively.b Obtained solutions of radial electric field as a function of a variable determining anomalous electron flux.The H-mode solution corresponds to lower balanced flux in the left graph (a), and the stronger electric field branch in the right graph (b)

Fig. 2
Fig.2Measurement of radial electric field in L-and H-mode in JFT-2M, which shows the clear difference of radial electric field profiles.a Radial electric field measured with charge exchange spectroscopy using force balance equation(Ida et al. 1990), and b that measured with a heavy ion beam probe(Ido et al. 2000)

Fig. 3 a
Fig. 3 a Top and bottom panels show poloidal shock structure and induced inward particle flux, predicted in a simulation.The particle fluxes induced by asymmetric electric field structure (or convective fluxes) can explain rapid increase in density observed during H-mode transition (Kasuya and Itoh 2005).b The relation between density gradient scale and E r -structure parameter, Z = (E � r ) 2 − E r E ��, where E ′′ represents the E r -curvature.The density gradient becomes larger as the parameter Z increases (JT60U)(Kamiya et al. 2017).c Temporal behavior of turbulence and the E r -curvature in LHD.The top and bottom panels show temporal evolution of the E r -curvature and density fluctuations.It is clear that the turbulence amplitude increases when E r -curvature is expected to enhance turbulence(Tokuzawa et al. 2014) The project was initiated by the first grant entitled 'Research on Structural Formation and Selection Rules in Turbulent Plasmas' from Japanese fiscal year 2004 to 2008, then it was extended by the second one 'Integrated Research on Dynamic Response and Transport in Turbulent Plasmas' from 2009 to 2013.Two notable enterprises were undertaken for the period.One is the creation of 'Kyushu University Itoh Project Prize' in 2005 in collaboration with Institute of Physics Publishing and Division of Plasma Physics of European Physics Society, which aimed at finding promising young scientists in her/his early career.The other one is the founding of her own research center, named 'Research Center for Plasma Turbulence' in 2009.Many worldwide collaborators, including Profs.R. Balescu, G. Bonhomme, R. Dendy, P. Diamond, O. Grulke, K. Hallaschek, T. Klinger, H. Kamitsubo, T. Ohkawa, R. D. Sagdeev, U. Stroth, H. Toki, M. Q. Tran, G. Tynan, F. Wagner, joined the project.

Fig. 4
Fig. 4 Identification of zonal flows in CHS (Fujisawa et al. 2004).a Spectra of electric field fluctuations.b Spatiotemporal correlation functions of electric field between two toroidal positions.The dashed line shows the position where the correlation value is zero.c The spatial correlation functions at = 0 , 1, and 2 ms.The shapes of the deduced correlation function correspond to the structures of zonal flows

Fig. 5
Fig. 5 Turbulence trapping and spreading due to zonal flows.a Experimental results to suggest turbulence trapping due to stationary zonal flows (Fujisawa et al. 2007a).b Simulation results of turbulence spreading caused by GAM trapping.The left panels show the wave number distribution as a function of radius at the initial and turbulence saturated state, while the right panel does the turbulence propagation in the radius/time domain (Sasaki et al. 2018).It is clearly shown that the trapped turbulence is transferred inward research) entitled 'Research for New Paradigm of Transport due to Phase Space Fluctuations in Fusion Plasma' from 2021 to 2025.