Newly uncovered physics of MHD instabilities using 2-D electron cyclotron emission imaging system in toroidal plasmas

ABSTRACT Validation of physics models using the newly uncovered physics with a 2-D electron cyclotron emission imaging (ECEi) system for magnetic fusion plasmas has either enhanced the confidence or substantially improved the modeling capability. The discarded “full reconnection model” in sawtooth instability is vindicated and established that symmetry and magnetic shear of the 1/1 kink mode are critical parameters in sawtooth instability. For the 2/1 instability, it is demonstrated that the 2-D data can determine critical physics parameters with a high confidence and the measured anisotropic distribution of the turbulence and its flow in presence of the 2/1 island is validated by the modelled potential and gyro-kinetic calculation. The validation process of the measured reversed-shear Alfvéneigenmode (RSAE) structures has improved deficiencies of prior models. The 2-D images of internal structure of the ELMs and turbulence induced by the resonant magnetic perturbation (RMP) have provided an opportunity to establish firm physics basis of the ELM instability and role of RMPs. The importance of symmetry in determining the reconnection time scale and role of magnetic shear of the 1/1 kink mode in sawtooth instability may be relevant to the underlying physics of the violent kink instability of the filament ropes in a solar flare. Graphical abstract


Introduction
The Magneto-hydro-dynamic (MHD) and Alfvénic instabilities in solar, space and laboratory plasmas have attracted researchers for many decades. These instabilities are closely associated with explosive events such as disruptions in tokamak plasmas [1], coronal mass ejection (CME) in solar flares [2], and geomagnetic storms [3]. The physics of enhanced energy transport and disruptions induced by the MHD and Alfvénic instabilities in fusion plasma devices has been vigorously studied to achieve improved confinement and stability of the plasmas. The ultimate goals are to perfect theoretical models with a predictive capability and to achieve preemptive control of harmful instabilities. Such progress has been accelerated through a precise validation with multi-dimensional advanced imaging diagnostic systems in recent years. The synergy between theoretical modelings and multi-dimensional experimental constraints in the validation process is essential for physics understanding of complex plasma dynamics. For example, new findings in magnetic reconnection experiment with sophisticated multi-dimensional probe arrays in laboratories [4,5] have partly validated the physics of reconnection events in the solar flare and geomagnetic storms. Note that such probe array measurements are useful in low-temperature plasmas, and the physics can be extrapolated to different plasma regimes with a caution. In geomagnetic plasma studies, unresolved physics will be benefited by the magnetospheric multiscale (MMS) program [6]. Presently the performance of fusion plasmas has routinely reached to a temperature above~10keV at a high plasma density and International Thermonuclear Experimental Reactor (ITER) based on the tokamak concept is under construction to demonstrate a burning plasma physics creating significant fusion power (~500 MW). Effective control of the critical MHD and Alfvénic instabilities that affect the stability and confinement of the plasma [7] can be achieved with validated models. This will ensure the success of ITER and future fusion reactors will greatly benefit from this knowledge. Since the present theoretical modelings have been validated with conventional diagnostic information [8], deficiencies in modeling can be further clarified with advanced diagnostics. Therefore, multi-dimensional measurements with adequate spatial and time resolution in fusion plasma research are critical so that the physics validation can be performed with precision and tight constraints. The 2-D electron cyclotron emission imaging (ECEi) system [9] has been developed for such needs, and the newly-uncovered physics has provided an opportunity to improve the theoretical models for toroidal fusion devices. This paper is organized as follows: In Section 2, a brief review of toroidal fusion devices is given with a definition of the rotational transform (magnetic field helicity) which is essential for a stable equilibrium in these devices. Among the MHD instabilities driven by the current and pressure gradient in tokamak plasmas [7], the sawtooth, neoclassical tearing mode (NTM), Alfvén eigenmode and edge localized mode (ELM) instabilities are addressed [7]. In Section 3, advantages and shortfalls of diagnostic techniques employed for the study of these critical instabilities are discussed: They are Faraday rotation polarimetry and motional Stark effect (MSE) spectroscopy for the current density profile measurement and soft x-ray tomography for 2-D presentation of critical instabilities [8]. A brief review of the electron cyclotron emission (ECE) measurement of the electron temperature (T e ) profile and evolution toward the 2-D ECEi system for visualization of T e fluctuations follow. In Section 4, the newly-uncovered physics in the study of the critical instabilities with the 2-D ECEi diagnostic system are addressed together with the background of these instabilities. They are vindication of the 'full reconnection model' for the sawtooth instability, improved accuracy in determination of the critical parameters of the 2/1 mode and multiscale interaction between the 2/1 mode and micro-turbulence in NTM instability, deficiency in modeling of the Alfvénic eigenmode instability, and newly measured dynamics of the edge-localized modes (ELMs) and ELM-crash in 2-D will further improve understanding. Relevance of these new findings to the solar and space plasma physics is addressed in Section 5 and a summary section follows.

Review of fusion plasma devices and critical instabilities in tokamak plasmas
In the course of fusion plasma research, various concepts of toroidal fusion devices have been proposed. Among them, promising fusion plasma devices are merged into two groups; stellarator and tokamak families as illustrated in Figure 1. In both families, a stable equilibrium requires a helical magnetic field line (i.e. field line pitch) instead of straight solenoidal field in a closed torus. The field line pitch is defined as a 'rotational transform' (t/2π) in the stellarator and 'safety factor' (q) in the tokamak, and they are related by q = 2π/t [10]. In the stellarator family, the rotational transform is induced by sets of complex external coil systems as shown in Figure 1b, and it is non-trivial to change the rotational transform once the device is assembled. There are two ways to achieve the rotational transform in stellarators. The complex external coil structures can be arranged to generate rotation of elongated flux surfaces around the torus (e.g. LHD, Japan) or to make the magnetic axis non-planar (e.g. Wendelstein-7X, Germany). In the tokamak family (tokamaks and reverse field pinches), a helically axisymmetric nested magnetic surface is generated by superposition of an applied toroidal magnetic field (B T ) by external coils and poloidal magnetic field (B P ) induced by driven toroidal plasma current (I P ) as shown in Figure 1a. The degree of helicity proportional to (B p /B T ) at fixed geometry (R, r, and κ, where R is the major radius, r is the minor radius and κ is the elongation factor of the nested flux surface) is characterized by the safety factor q(r) as shown in Equation (1).
In the tokamak plasma, the safety factor profile, q(r), can easily be manipulated by an external current drive system to influence the stability and confinement of the plasma. Note that the current required for a stable equilibrium is also a source of instability in tokamak plasmas. The safety factor profile, q(r), contains a series of rational surfaces (q = m/n = 1, 2, 3, …, where m and n are poloidal and toroidal mode number, respectively) on a nested poloidal surface. On a rational surface, the field line returns to the same point after m and n of the poloidal and toroidal angles, respectively (e.g. on the q = 1 surface, the point fully wraps around poloidally for one full toroidal turn), and a resonant mode can grow on these rational surfaces. Since the safety factor profile, q(r), is the key parameter for the stability of critical MHD and Alfvénic modes, it is extremely important to measure the q(r) accurately with a good spatial and temporal resolution. The plasma pressure profile gradient in the toroidal plasma is mostly positive definite toward the center of the plasma and free energy from it is a potential source of instabilities. There are various modes of plasma operation and distinctive profile characteristics can be found either in the core or edge of the plasma pressure profile. Historically, tokamak plasma operation was started with relatively low energy confinement mode (L-mode) where the edge pedestal height was low (i.e. moderate pressure gradient in this region). As the shape of the plasma was evolved from a circular plasma with the plasma in contact with a material surface (a 'limiter' plasma) to an elongated plasma with magnetic divertors, a high confinement mode (H-mode) with high pedestal pressure [11] was discovered in the ASDEX device in Germany and the energy confinement time was almost twice higher compared to that of the L-mode for the same heating power. Note that there is a wide variation in confinement time within L and H mode plasmas. While high edge pedestal in the H-mode plasma has the advantage of higher confinement and bootstrap current, higher free energy due to steep pressure gradient and large bootstrap current at the edge can each be a source of edge instabilities. At the edge of H-mode plasmas, the edge-localized mode ('ELM'), characterized by a periodic burst of D α light due to a partial collapse of the edge pedestal pressure, has been observed. Another improved confinement mode under study is the internal transport barrier (ITB) mode [12], where the core confinement is significantly enhanced with the L-mode edge plasmas. This mode of operation has a relatively broad current profile due to the heating beam driven current. The ITB mode with high-pressure gradient in the core of the plasma can excite internal instabilities such as the NTM instability.
The free energy from the pressure gradient or energetic particles at rational surfaces of the q profile can also drive benign and/or harmful resonant MHD and Alfvénic instabilities. Ironically, the harmful modes driven by the plasma current which is required for a stable equilibrium, limit the increase of the confined energy due to enhanced transport and a rapid growth of these modes often leads to minor and/or major disruptions. Therefore, effort for an effective control of these harmful instabilities has been continued, but the progress has been mainly empirical. In order to develop a predictive capability in theoretical modelings for full control of these instabilities, fully validated theoretical models are needed. This capability is essential for future advanced fusion devices, where the plasma pressure and current density will be significantly higher than those of the present devices. Furthermore, any new underlying physics of these MHD models validated in fusion plasma devices can also be applied for interpretation of similar MHD phenomena in a nature such as explosive behavior of magnetic flux ropes in a solar flare. Among MHD instabilities at the resonant rational surfaces, the critical MHD instabilities are sawtooth oscillation (m/n = 1/1 mode in kink or tearing) within the q~1 surface in the core, neoclassical tearing modes (NTM) (2/1 tearing mode) near the q~2 layer and edge localized mode (ELM) (high m/n mode) at the pedestal of the H-mode plasma where the pressure gradient is extremely steep and the edge current profile is significantly modified by the bootstrap current. Alfvénic instabilities, such as the reverse shear Alfvénic eigenmode (RSAE) and toroidal Alfvénic eigenmode (TAE) are driven by the energetic particles.

Conventional and advanced diagnostics essential for the study of critical MHD and Alfvénic instabilities
In the course of studying the physics of MHD and Alfvénic instabilities, various diagnostic tools [8] have been employed to measure physical quantities that can influence the stability and confinement of the current-carrying plasma. To monitor the current flow of the plasma, various types of magnetic probes positioned at various places on the vacuum vessel wall are employed to measure the change of the magnetic flux induced by the plasma current flow and its fluctuations [8]. They are simple and easy to be deployed but a drawback is lack of internal plasma information. Therefore, a local measurement of the plasma current profile (or safety factor profile) with an adequate spatial and temporal resolution is essential to understand the dynamics of the current flow in these instabilities.
High-temperature plasmas in a magnetic field create a wide range of electromagnetic emissions from radio frequency to gamma ray. Besides the cyclotron emission from gyration motion of the charged particles, emissions originate from various mechanisms such as atomic transitions in the main and impurity particle species, acceleration in a transiently developed electric field, bremsstrahlung, etc. These emissions have been widely used to characterize the events in plasmas as well as to monitor instabilities. Among them, soft x-ray (SXR) emission [8] has been widely used for the study of MHD and Alfvénic instabilities, since SXR has a strong dependence on the electron temperature (T e ) and implementation is relatively easy compare to the active systems that require sophisticated instrumentation which may not be easily deployed in some cases. Since comprehensive measurement with high spatiotemporal resolution [9] is essential to the study the localized MHD and Alfvénic instabilities, tomographic techniques perfected in the medical field have been applied but precise physics information may be limited due to accessibility constraints (i.e. the diagnostic requires many views for a unique tomographical solution) and the dependence of the emission on more than one plasma parameter.
In order to address fast time dynamics of the localized MHD modes and Alfvén waves in linear and non-linear phases, both high spatial and temporal resolutions are essential. Therefore, a fast detection of the electron cyclotron emission (ECE) has been widely used for this application [8]. Note that the electron temperature is known to be a flux surface quantity, since electrons can rapidly follow changes of the local magnetic field due to its high thermal velocity. For example, the timescale of thermal electron motion with T e~1 keV at B T~2 T is a nano-second timescale whereas the upper bound timescale of the MHD instabilities is~μ-second timescale, and the particle collision frequency is in similar timescale.

Measurement of the current density profile in tokamak plasmas
As noted in the previous section, precise measurement of the plasma current density and pressure profile is critical for stability analysis of the MHD and Alfvénic instabilities in tokamak plasmas. The confidence in measurement of the plasma pressure profile (i.e. plasma density and temperature profiles) in tokamak plasmas is well established and will not be discussed in this paper. There are two principal methods for measurement of the current density (or equivalently, poloidal field) in tokamak plasmas and both methods are technically complex and many factors can limit the required accuracy for the stability analysis: polarimetric measurement based on Faraday rotation [13][14][15] and spectroscopy based on the motional Stark effect (MSE) [16][17][18][19][20][21]. Since the required accuracy level varies for different instabilities under study, it is important to be aware of limitations of these methods which are critical for the model validation process.
The polarimeter based on the Faraday rotation [13][14][15] requires a double inversion process, since the measured Faraday rotation angle is a convolution of the plasma density and magnetic field along the path of an electromagnetic wave propagating through the plasma. Also, the unfolding process relies on the accuracy of the equilibrium geometry. As an example, the most successful measurement of the current density profile was the polarimetry system on TEXTOR [13] where the plasma shape was circular. The unfolding process is easier in the plasma in a circular cross-section compared to a shaped plasma. Much attention was given to instrument calibration, and various tests were performed in order to minimize errors in measurement and the double inversion process as illustrated in Figure 2. Here, chordal arrangement of the probe beams on the poloidal cross-section is given. The increased data points of the line-integrated electron density and Faraday rotation angle are achieved by a radial 'jogging' of the plasma to enhance the spatial resolution. Then, double inversion process is applied to obtain the safety factor profile, q (r), as shown with q 0~0 .75 in this figure. However, an intrinsic weakness of the polarimetry is that the signal is typically small near the magnetic axis, where the poloidal field (B p ) approaches zero, and also near the edge of the plasma where the plasma density is close to zero. Since the q value at the edge of the tokamak plasma is usually a large number (>3), relative accuracy can be reasonable. On the other hand, the q value near the axis, including q 0 , is small (q~1.0) in a positive definite current profile; then, the relative error becomes large. Note that the q 0 is proportional to the derivative of the measured Faraday rotation angle profile across the axis. In addition, an application of this technique to shaped plasmas (plasma with finite elongation) further complicates the analysis due to the nature of chordal measurement and complexity of the instrumentation.
The spectroscopic method based on the Motional Stark Effect (MSE) [16][17][18][19][20][21] and Zeeman splitting [22] measures polarization angle of the MSE or Zeeman splitting and the resultant pitch angle measured depends on system geometry from equilibrium reconstruction. The nature of the MSE spectroscopy is the local measurement of the magnetic field line pitch angle by measuring the polarization of the emission of the injected neutral beam (which is usually the neutral beam used for heating the plasma) with known energy as illustrated in Figure 3. Here, the polarization of the emission is sensitive to the local electric field induced by the local magnetic field direction and velocity of the beam ions. The MSE system has been evolved technically to allow more reliable measurement and is now more routinely employed in contemporary tokamak plasmas. The poloidal field information is based on knowledge of the equilibrium Figure 2. Illustration of double inversion process of the measured Faraday rotation angles using 9 chords of probe beams. In order to increase the precision, radial jogging of the plasma was performed as shown in the plot of line densities and Faraday rotation angles (left side). Then Abel inversion is introduced to invert the line integral density into the local density. Then the local density profile is used to evaluate the poloidal field which is used to calculate safety factor profile (right side).
Source: partial figure from H. Soltwisch magnetic geometry as shown in Figure 3. Here, the measured poloidal field (B z ) is compared with the calculated value from equilibrium reconstruction and the agreement is excellent. The measured q profiles before (q 0 = 0.97) and after (q 0 = 1.05) the sawtooth crash are illustrated in this figure. Similar to the polarimetry, the MSE measurement faces the same sensitivity issue near the axis and edge, since the polarization angle of the π component emission is proportional to the field line pitch angle, tan (B z /B T ) and the B z is again zero on axis, and the signal level becomes poor at the edge due to low emission level (i.e. low edge density). There are a number of sources of error and corrections in this measurement. The first one could be the Zeeman effect and the estimated ratio of the Stark splitting to the Zeeman splitting is~4 for the signal that MSE measures from a hydrogen atom experiencing Balmer-alpha transition. Therefore, the Zeeman effect may not be a big issue but caution is needed for small signals. Also note that there is an additional intrinsic electric field due to ambipolar diffusion of the plasma particles [19] which can impact the measured polarization and the detected emission signal has to go through a careful subtraction of an unwanted background emission including reflected signal from metal walls [21]. At the same time, an accurate magnetic geometry from a reconstructed equilibrium, which needs the current profile information, is critical for confidence of the measured poloidal magnetic field. For instance, the pitch angle has a dependence on an elongation of the magnetic surface. Therefore, the measured q 0 value can be subject to significant error due to this and other effects. Viewing angles with respect to the plasma are shown with typical positions of the axis and separatrix are~1.55m and~2.13m, respectively. The MSE measured poloidal field (B z ) and B z from equilibrium reconstruction (solid line) are compared just before the sawtooth crash. The q profiles for before (q 0 = 0.97) and after (q 0 = 1.05) the sawtooth crash are ploted.
Source: B.W., Rice, Fusion Engineering and Design, [34][35]1997 3.2. Soft X-ray emission measurement for monitoring MHD instabilities In high-temperature plasmas in the tokamak, it is difficult to avoid high-Z impurities due influx of wall materials such as C, Fe, Ni, etc. which can diffuse into the core plasma. The emission from high Z elements at high temperature are from atomic transitions from the outer shell to hydrogenlike elements. The emission near the soft X-ray (SXR) region avoiding line radiation has been well understood and has a strong dependence on T e with relatively moderate dependence on n e and Z eff . The SXR detection system is capable of high time-resolution on the order of a μ-second, which is well-suited for the study of T e fluctuations in the core of tokamak plasmas. Note that the first observation of the sawtooth oscillation was by SXR emission [23]. Therefore, the SXR spectroscopy system has been widely used in studying the critical MHD instabilities such as the sawtooth oscillation and NTMs in tokamak plasmas [23][24][25][26][27][28]. Since MHD instabilities in tokamak plasmas are located near their corresponding rational surfaces as described in the previous section, analysis of the amplitude and phase information along the chords can provide information regarding the q-profile. Furthermore, one or more poloidal arrays of the SXR emission detection can be used to construct a 2-D image of MHD dynamics [23][24][25][26][27] using a tomographic technique. Such a reconstruction from the line-integrated measurement demands a unique solution of the inverse problem. However, often times, the solution to the relevant inverse problem is not unique with the limited views; so multiple views are needed in order to find an unambiguous solution. The methodology has been perfected in the context of the Magnetic Resonance Imaging (MRI) where multiple views of an object (i.e. 360°view) and a long integration time (~minutes) are used to determine an accurate image of the object. The most widely used diagnostic system for the study of critical MHD and Alfvénic instabilities as well as the electron temperature profile in tokamak plasmas is the ECE measurement [8,[29][30][31][32][33][34][35][36], which provides a spatial and temporal evolution of the local electron temperature and its fluctuations. In a magnetized toroidal plasma, the electron gyro-motion results in an emission at the electron cyclotron frequency, ω ¼ eB m and its harmonics [8], where B is the magnetic field strength, e is the electron charge, and m is the electron mass. The principle of ECE for T e measurement is well established, and this technique has been routinely used to measure the local electron temperature in fusion plasmas for over four decades, since the first measurement on the CLEO tokamak [29].
When the electron density and temperature are sufficiently high, the plasma is considered to be optically thick [8,37] with respect to the first (ordinary) and second harmonic (extraordinary) mode. Then, the emission intensity, I(ω), approaches to that of black body radiation which is given by Equation (2).
where ω is the radiation frequency, c is the speed of light and T e is the electron temperature. The optical thickness in ECE has been extensively studied and a guideline for a valid T e measurement for given plasma density and temperature is provided in Ref [8,37]. In the regime, where the optical thickness is reduced to the grey level, localized measurement is possible but the emission intensity is no longer proportional to the local T e alone (i.e. the emission depends on a combination of the electron density and temperature). When the optical thickness is lowered further and it becomes thin (i.e. at very edge of the plasma), the emission neither is localized nor represents the local T e and the measurement is contaminated by down-shifted emissions. Since the intensity of the emission is directly proportional to T e as long as the region is optically thick, the measurement of ECE at a fixed frequency enables a measurement of the local T e as well as its fluctuations. Here, the ECE frequency has a spatial dependence due to the radial dependence of the toroidal magnetic field, B T / 1=R, where B T and R are the toroidal field and major radius of the tokamak plasma, respectively. The advantage of ECE for T e measurement can be compromised somewhat by a number of factors: spatial variations introduced by the poloidal magnetic field, which can be a significant fraction of the toroidal magnetic field as the toroidal field is lowered and/or the plasma current is raised, by harmonic overlap when the same frequency is resonant at two different spatial locations in a low magnetic field device. A good example is an application of ECE measurement to a spherical tokamak where the toroidal field is comparable to the poloidal field. The other effect comes from a wave refraction effect which may introduce uncertainties in viewing directions when the plasma density gradient is high. Allowances need to be made for these potential sources of errors in interpretation of the observed emission. The local T e can be obtained from the measurement of ECE in a number of ways. These include usage of the heterodyne ECE receiver [30], grating polychromator [31] and Fourier transform spectrometer [32]. The Fourier transform spectrometer [32] has the advantage of ease of calibration due to a large throughput but has poor time resolution. In the grating polychromator system [31], the time resolution is significantly improved, but achieving an absolute calibration is difficult. In the heterodyne system [30], a single or two-stage down-converter is required to process the high-frequency signals to a convenient intermediate frequency (IF) for a final detection. Among the three techniques, the heterodyne system yields the highest sensitivity with absolute calibration. The principal limitation of this system is the spatial coverage due to bandwidth limitation and harmonic overlaps. In order to cover the entire plasma, one may need more than one detection system. A schematic of a typical arrangement of instruments required for a 1-D heterodyne ECE radiometer is shown in Figure 4a. Here, emissions from the plasma on the mid-plane, where the refraction effect is minimum, are collected by a single detector through optics. The detected emissions are splitted with adequate radial spacings and then they are down converted with local oscillator (LO) sources to measure the local T e . A typical resolution of the 1-D ECE system is 5 cm x 5 cm. Since the ECE system is a local measurement with high time resolution, the measured T e across the plasma provides flux surface position and corresponding T e profile. Note that T e is constant on the same flux surface (i.e. flux surface quantity), since the electron has a high thermal velocity. In a toroidally rotating plasma, toroidally axisymmetric MHD structures can be reconstructed to form a 2-D image using the time-dependent 1-D data with an assumption that the MHD dynamics is invariant during one full rotation.
In order to investigate non-axisymmetric and fast transient dynamics of MHD mode activity and related turbulence (i.e. growth of the 1/1 kink, 2/1 TM and NTM, reconnection events during sawtooth dynamics and ELM-crash, disruption, etc.), it is desirable to develop a 2-D ECE imaging (ECEi) system which can simultaneously measure 2-D images of the MHD activity with high spatial and time resolution (e.g. instead of~5 x 5 cm in 1-D ECE, it could bẽ 1.0 cm x 1.0 cm). This requires improved sensitivity of the detection system and integrated signal processing technology. As the sensitivity of the detection system has been improved and the size miniaturized with advances in submillimeter technology [38], the single detection element is replaced with a 1-D array detector where each detection element is vertically stacked as shown in Figure 4b. As shown in this figure, the vertically extended large optical elements (e.g. triplet lens system with zooming capability in the KSTAR system) are employed to have a one-to-one mapping of the vertically resolved beamlines of the emissions from the plasma within the focal depth to each detection element in the vertical array through a mini-lens to collect the emissions effectively and simultaneously. The detection element consists of a Schottky diode and dipole antenna structure as shown in the figure. The measured emission at the ECE frequency (~100 GHz corresponds to second harmonics emission) is downconverted using a LO to the desired intermediate frequency (IF) range (2-9 GHz). The IF signals from each detection element are divided into eight channels with 700-900 MHz bandwidth and processed simultaneously with state-of-the-art electronics to form a 2-D image [9,36]. The first prototype 2-D ECEi system with a vertical array of 16 detection elements was developed for the TEXTOR device [9]. Here, the ECEi system was combined with the microwave imaging reflectometer (MIR) for simultaneous measurement of the fluctuations of n e and T e . The sampling area at the focal plane was~20 cm (vertical) x~8 cm (radial), and the vertical resolution was determined by the antenna pattern of  resolution was~1 cm across the core of the tokamak plasma (total of 128 [16 x 8] channels). The time resolution was primarily limited by the digitizer, and the nominal time resolution was~5 μsec. The time-dependent signals from 128 channels were simultaneously processed to assemble a 2-D image. In principle, each channel can be absolutely calibrated so that it can measure the T e profile as well as the fluctuation. In practice, there are difficulties in calibrating all channels due to the variation in the local oscillator power at various frequencies.
So the processed image is based on a relatively calibrated value (δT e = T e /<T e >) for each channel, where <T e > is the time-averaged value for a certain time zone.
In recent years, the microwave technology has been further improved and an integrated system with a highly sensitive detection element, and primary electronics has been developed at the lower frequency where the technology is available [39]. This technology will be valuable for future application. When the first ECEi system was installed on the TEXTOR tokamak plasma, validation of the system performance was tested against the sawtooth instability, since dynamics of the 1/1 kink mode in a sawtooth oscillation were well known and the size of the cold island and 1/1 kink mode were large enough to crosscheck the spatial and temporal resolutions of the system as shown in Figure 5. The poloidal cross-section of the TEXTOR plasma with the normalized change of T e (δT e /<T e >) before and after the crash of the sawtooth across q~1 surface is shown in this figure. The monotonically peaked T e profile prior to the crash (yellow color) across the q~1 surface (double white line) and increased T e in the mixing zone and decreased T e (blue color) inside the double white line right after the crash (yellow color) are well represented by 2-D ECEi data. The 2-D images are consistent with the classical dynamics of the 1/1 kink mode with a time resolution of~5 μs. Following the first successful ECEi system on TEXTOR, the ECEi system has become a standard advanced diagnostic system on tokamaks as well as stellarators; ASDEX-U [40], DIII-D [41], EAST [42], HL-2A [43], LHD [44], KSTAR [45], and WEST [46]. On KSTAR, the first ECEi system with two poloidal views and 24 vertical array detectors was commissioned in 2008. The second system with a single view, toroidally separated by 22.9°, was added to form a semi 3-D system along with the first system. The nominal view window is~40 cm (vertical) x~12 cm (radial) with the radial resolution of~1.5 cm and the vertical zoom factor of~3 (total of 192 channels). The time resolution can be adjusted by the digitizing speed, and the nominal time speed is~2 μsec. Two simultaneously measured images of an ELM from two edge toroidal views are illustrated in Figure 6. An important information such as the magnetic pitch angle, velocities, and self-consistent mode numbers can be comprehensively obtained as illustrated in this figure. In each tokamak device, the ECEi system has produced many new physics observations that were not available with conventional diagnostic systems. The sawtooth instability (slow rise of the core plasma pressure followed by a sudden growth an m/n = 1/1 kink (or tearing) mode which leads to a fast 'crash' of the measured emission) is the most fundamental tokamak instability and occurs in the core (inside the q~1 surface) of the plasma. This instability is benign as long as the amplitude of the crash is small. As the amplitude of the crash is increased, a sudden burst of the plasma pressure and perhaps current density at the q~1 surface generates an outward ballistic perturbation that can be a 'seed' of harmful MHD like NTMs and ELMs. Despite of a long history of study, the stability criteria have been disputed until recently. This instability was first discovered in the ST tokamak at Princeton [23] in which the soft X-ray detection system revealed an oscillation pattern similar to a 'sawteeth-shaped' pattern as shown in Figure 5. Since the observed crash timescale was significantly faster than the transport timescale, other mechanisms such as magnetic field reconnection or an interchange mode instability was suspected. Figure 6. Arrangement of the KSTAR quasi 3-D ECEi system is illustrated. The 1st ECEi system equipped with two views for simultaneous measurement at two poloidal planes (e.g. core/ edge or high field side/low field side) is shown in red color box. The 2nd ECEi system with blue color is added at the toroidal plane separated by 22.9°and simultaneously measured ELM images from two edge views are shown. Here, the pitch angle, velocity and mode numbers can easily be calculated.

Source: H. Park
Following the first observation of this instability, Kadomstev proposed the 'full (or complete) reconnection model' [47]. In this model, an excess core current density and plasma pressure responsible for the 1/1 kink (or tearing type) instability are fully removed through a magnetic reconnection event and the removed excess current and pressure are piled up at the mixing zone in a characteristic timescale, , where τ r is the resistive time and τ A is the Alfvén time. The timescale estimation was based on the 'Sweet-Parker' reconnection model [48]. After the crash, the central current density slowly builds up with the rising central temperature, a displacement of the 1/1 kink mode occurs as the q 0 drops below 1.0 along the helically symmetric surface within the q~1 surface. Then, the reconnection event occurs at the region where the pressure gradient is steepest on the poloidal plane. Here, the reconnection zone is helically axisymmetric on the toroidal plane and the q 0 returns above~1.0 after the discharge of the excess current and pressure is completed. An example of 1-D T e profile during the sawooth oscillation and the detailed 2-D electron temperature (T e ) transport event occurring over a~100μs interval is shown in Figure 5. Increased T e is represented in a yellow color in frame 1 (increase in T e ) before crash event. The increased T e in the mixing zone is shown in the frame 3 as yellow color with bluish color in the core (decrease of T e ) after the crash event is illustrated in this figure. There is a discrepancy in the time of the dynamics based on theory: the measured crash time of~100μs does not agree with the expected τ c~a few ms. This is one of the reason why the 'full reconnection model' was abandoned.
After the proposed full reconnection model, it took more than a decade to measure the current density (or q) profile in tokamak plasmas. The measurements of the q 0 value during sawtooth oscillation by polarimeters (TEXTOR [13], JET [14] and MTX [15]) and Spectroscopic measurements (MSE on TFTR tokamak [17] and Zeeman effect on TEXT [22]) reported that the measured q 0 was q 0 ,0:75 AE 0:05 which was well below~1.0 even after the crash. An immediate dilemma was how the 1/1 kink mode can be stable while q 0 is well below~1.0 after the crash. Numerous new theoretical models called 'incomplete (or partial) reconnection models' were proposed [49,50]. Meantime, technical advances of MSE spectroscopy were made, and new measurements were conducted on DIII-D tokamak device [18,19]. Here, they reported the measured q 0 ,1:0 AE 0:05 and supported the full reconnection model. Note that there was an earlier MSE measurement of q 0~1 .0 on TEXT [16]. The MHD community was puzzled about the two vastly different measurements. One suspicion was that the flux surface elongation might explain the difference in q 0 , since the first two measurements (~0.75) were from circular plasmas (TFTR) and the others (~1.0) were from the shaped plasmas (DIII-D), where < > represents average value. Note that there were measurements of mixtures of~0.75 [21] and~1.0 [16] on TEXT (circular plasma) and the group who measured~0.75 in TFTR [17] reported~0.86 in PBX-M device where the core κ was~1.36 [51] and the dependence of q 0 on the plasma shape factor may not be convincing. The experiment conducted on DIII-D [52] to compare the sawtooth behavior in vastly different plasma shapes (bean and oval) could not find any further resolution. In recent experiments in plasmas with elongation of~1.2 on KSTAR, MSE spectroscopy reconfirmed that q 0 ,1:0 AE 0:03 in the plasma with sawtooth oscillation as shown in Figure 7 [53].
In summary, the measurement of q 0 in plasmas with sawtooth oscillations is divided into two groups. They are q 0 ,1:0 AE 0:03 and q 0 ,0:75 AE 0:05. Here, the difference in median value of two group is Δq~0.25 with error range of AE 0:03. However, the relative change of q 0 (i.e. δq~0.07) before and after the crash is common for both groups and this change is consistent with the expected small current diffusion after the crash. The difference, Δq~0.2, between two groups cannot be explained by the estimated experimental error. Even if the second group (~1.0) is correct, it is still difficult to validate the model, since the experimental error is too large or the required absolute error is too stringent to address that the q 0 is above~1.0 right after the crash. Therefore, it is imperative to introduce a comprehensive supplementary experiment which can identify, if q 0 is above or below~1.0, because a comprehensible validation of the model with a direct measurement is not possible with present-day instruments. This statement does not imply that the MSE and Faraday rotation measurements are inaccurate. In fact, they are extremely accurate in spite of the complexity of instrumentation except that the change in current density is relatively small in the sawtooth dynamic as discussed in the previous section. The observed crash time of the sawtooth instability has been unanimously faster than the proposed τ c by the full reconnection model even in a large tokamak [54]. This was the primary reason why the full reconnection model was abandoned as soon as the experimental observation was determined even before the q 0 issue became controversial. A decade later, a new model, the 'quasi interchange mode model' was proposed by Wesson [55] with an assumption that the q profile is flat at 1.0 inside the q~1 surface (no magnetic shear) to explain the observed fast crash timescale. Note that this model does not require field line reconnection and is based on the interchange mode instability developed at the vulnerable region of the q profile (i.e. dq/dr is maximum near the edge of the q = 1 surface). This model was supported by observation of the interchange mode with the SXR tomography system on the JET tokamak device [25]. A later experiment on the TCV device demonstrated that a unique solution with the limited views of tomography might not be feasible [26]. However, the probability might have been low, but there was a chance to have the q = 1 condition in the core so that the observed interchange mode could explain the observations. Certainly, this model can explain the fast reconnection time based on observed interchange mode, but the majority of cases are consistent with the full reconnection model where the field line reconnection appears to be due to the 1/1 kink (or tearing) mode rather than interchange mode. Recently, a plasmoid model that has been popular in explaining the fast reconnection process in astrophysical plasmas and solar flares has been used as an attempt to explain the fast reconnection process in sawtooth oscillations [56].
Detailed analysis of the sawtooth crash performed on TFTR reported many interesting observations through reconstructed 2-D image with 1-D ECE data using the plasma rotation [57]. While most of the observed sawtooth crash processes were consistent with the full reconnection model in which the reconnection happened at bad curvature side (i.e. low field side) except that the reconnection zone is localized along the magnetic pitch on the toroidal plane (i.e. 3-D nature) with a fast reconnection time compared to the τ c . Here, the hot spot (1/1 kink or tearing mode) had a feature of the ballooning mode which was discovered for the first time on TFTR [58]. In addition, examples of slow reconnection processes were observed in ohmic and ion cyclotron radio-frequency heated (ICRH) plasmas. Here, the reconstructed 2-D image revealed that the 1/1 kink mode was axisymmetric (i.e. 2-D nature) and reconstruction of 2-D image for these cases was difficult due to the slow reconnection time compared to the plasma rotation speed. The core plasma pressure was expelled and flattened but little change of the core current density was noted after the crash based on the measurement from Reference [17].
While the sawtooth pattern has a slow rise of T e in an MHD quiescent period after the crash and a rapid growth of the 1/1 kink mode before the crash (i.e. precursor), there have been many reports on variations of the sawtooth pattern such as middle cursor, post-cursor, compound and humpback sawteeth [59,60]. For instance, the middle cursor and compound sawtooth [57,61] has the 1/1 kink oscillation in the middle of the rising phase of T e followed by small crashes before the main crash. Another example is the post-cursor oscillation which has a long decay time of the 1/1 kink mode after the first crash [57,60]. The physical interpretation of these variations may be self-explanatory, once the basic physical mechanisms, such as the sensitivity of the 1/1 kink mode on q 0 and its shear are understood. This subject will be revisited in the following section. Prior to visualization of the sawtooth instability with 2-D ECEi imaging, the 2-D crash pattern reconstructed with the SXR and ECE tomography demonstrated that the crash was likely to be observed at the lower field side in an axisymmetric toroidal plane due to the ballooning nature of the 1/1 kink mode [26,58,61]. When the 2-D ECEi system was introduced in the study of the sawooth crash process, it was found that the crash could occur localized at the high field side as well as the low field side. As demonstrated in Figure 8a, a 'finger' with bulge (i.e. ballooning mode) with a highly distorted 1/1 kink mode develops before the crash event and the 'finger' punches through the q~1 surface at the high field side. The core temperature (T e ) is removed through a narrow reconnection zone on the poloidal plane and then the T e profile is flattened inside the q~1 surface after the crash [62]. A similar crash pattern was also observed at the low field side as shown in Figure 8b. The other characteristic of the ballooning mode is a finite reconnection zone on the toroidal plane, and this was confirmed by the observed 2-D image during the crash time. These localized crashes, which occur randomly along the q~1 surface, may have been difficult to be detect with the 1-D measurement on the horizontal mid-plane. An example of the ballooning mode observed with the conventional 1-D ECE system in TFTR is shown in Figure 9 [57,58,61]. Here, two toroidally separated polychrometer systems measured T e on the midplane and a 'finger' structure of the ballooning mode is demonstrated on the low field side of the plasma. Since the motion of the narrow reconnection zone moves along the helical magnetic field line, the observation of it on the mid-plane can have a very short duration. With the time resolved 2-D images, the toroidal extent of the 1/1 kink mode with ballooning signature (bulge of T e at the reconnection zone) in a toroidally rotating plasma could be estimated to have a length of~1/3 of the plasma circumference (2πR 0 ) [63]. In summary, the observed bulge or finger during the crash process of the sawtooth instability is 'ballooning like' event, and it is not necessarily the classical ballooning mode.
The advent of a direct measurement of the 2-D image of the T e perturbation finds an opportunity to compare with the 2-D image from simulation. A direct comparison between the measured partial images of the 1/1 kink mode/cold island and simulated images from the 'full reconnection model' and 'quasi-interchange mode model' is demonstrated in Figure 10 [64]. The ECEi results show a closer resemblance to the full reconnection model. Since the first measurement of 2-D images of the sawtooth dynamics in TEXTOR, the sawtooth crash process has been studied with the 2-D ECEi system on DIII-D [65], EAST [66], ASDEX-U [67] and HL-2A [68], HT-7 [69] and each reporting that the observations are consistent with the findings from the TEXTOR plasma. The measured 2-D images from the sawtooth crash process in EAST, ASDEX-U and HT-7 are illustrated in Figure. 11a-11c. The fact that the observed sawtooth crash processes appear closer to the full reconnection model rather than the quasi-interchange mode model does not exclude the possibility of the quasi-interchange mode model. As discussed in the previous section, if q = 1 is nearly constant inside the q~1 surface, this instability can grow as the model suggests [55] but the probability of a nearly flat q = 1 region is small. Also, the present form of the plasmoid model introduced to explain the fast reconnection time [56] may produce a similar conclusion as the quasi-interchange mode model. Note that the crash process does not have to be unique and variations can be rare but are always feasible due to anomalies of the core current density and pressure in each crash event. It is important to focus on the physics of the dominant explanation of the event and modify the understanding by allowing deviations from the dominant event.
As discussed in the previous section, if the measured~0.75~0.75 is valid, we definitely need new comprehensive theoretical models, but there have been no new convincing theoretical models brought forward so far. On the other hand, even if the measured~1.0~1.0 valid, direct measurement of the current density profile alone cannot definitively validate the q 0 value before and after the crash event due to the required precision of the measurement of a small change of the core current density before and after the crash (~5%). A new comprehensive experiment, that can validate whether q 0 is above or below~1.0 right after the crash with no other MHD instabilities present, should be conducted with today's measurement capabilities. A couple of ideas have been suggested: The first one is based on measurement of the Alfvén eigenmode dispersion that is sensitive to the central q value and its shear [70]. The second one is excitation of the double tearing mode inside the q~1 surface. The second idea is similar to the off-axis sawtooth crash of the double tearing mode near the q~2 surface [71]. Note that there was an independent theoretical study to explain double tearing modes inside the q~1 surface [72]. In order to design an experiment for excitation of the double tearing mode with the current blip by electron cyclotron current drive (ECCD) inside the q~1 surface, a simulation study [73] with the reduced MHD code M3D-C 1 [74] was performed to examine growth rates of resonant modes (m/n = 2/2, 3/3, 4/4, etc.) with the modeled q(r) modified by a current blip (i.e. dip in q profile at r/a = 0.2) within the q~1 surface as shown in Figure 12a. For the background q 0 < 1.0 (= 0.98, and 0.80), there is no double-crossing points, and the m/n = 1/1 mode is the only mode that has significant growth rate [73]. When q 0 = 1.04, there are double-crossing points across the q = 1 layer and growth rates for different resonant modes at different positions of the current blip are clearly demonstrated in Figure 12a. In order to establish the double-crossing dip in the q profile with the current blip in the core region of the KSTAR plasma, the ECCD system at 170 GHz was slowly scanned along the resonance layer (vertical direction) as illustrated in Figure 12b. Here, the trajectory and driven current of the launched ECCD was modeled with the TORAY code [75]. As illustrated in this figure, the observed tearing modes right after each crash event are well matched with the modes with the highest growth rate from the modeling for the q 0 = 1.04 case. These experimental results confirm that the q 0 after the crash is likely above 1.0 and this is consistent with the direct measurements from the second group (~1.0) as measured on KSTAR.
The time evolution of one sawtooth cycle without and with ECCD is shown in Figure 13. In the cycle without ECCD, the only observable instability is the rapid growing 1/1 kink mode before the crash. In the cycle with ECCD, the m/n = 3/3 tearing mode is excited right after the crash and transforms into the 2/2 mode. Then, the 2/2 mode transforms into the 1/1 mode which is followed by the crash. Based on simulation results, transformation of the 3/3 mode can be interpreted as follows: As T e builds up on the flattened T e profile after the crash (i.e. peaking of T e ), the increasing core current density induces the driven current blip to move toward the core (i.e. q 0 drops faster with shear) so that the 2/2 mode can grow. As the q value goes below 1.0, the growth rate of the 1/1 kink mode dominates. Note that it is common to have a discharge without sawtooth oscillation when q 0 is increased above~1.0. Among them, there is a 'hybrid mode' of operation where the q 0 is above slightly~1.0 due to broadened current profile in the core. When the same experiment was performed in a hybrid mode plasma, the excited double tearing modes (m/n = 3/3 and 5/5) did not change for~0.5 s as shown in Figure 5 of reference [73]. These results suggest that the double tearing modes was excited and sustained as Maximum growth rates of the resonant modes calculated using M3D-C 1 code are shown in the right side as the dip in the q profile with q 0 = 1.04 is scanned from the center to the edge of the q~1 surface (zone 1, 2, 3, and 4). (b) Excited resonant mode is shown when the current blip (dip in the q profile) is scanned from the center of the plasma to the edge of the q~1 using ECCD. In zone 1 (r/a < 0.15), a hot spot in the core is observed. In zone 2 (0.1< r/a < 0.22), the 2/2 mode is excited. In zone 3 (0.2< r/a < 0.27), the 3/3 mode is observed. Near the q = 1 region, higher order mode is excited. long as the q 0 remains above 1.0 and further support the full reconnection model.
In the 'full reconnection mode' model, the reconnection timescale (τ c ) is based on the helically symmetric 1/1 kink mode, and the reconnection zone is poloidally localized and helically axisymmetric along the toroidal plane (i.e. 2-D nature). As shown in Figures 8 and 11, the observed 1/1 kink mode near the crash time is highly distorted (i.e. not symmetric in either the poloidal or toroidal planes) and the reconnection zone is highly localized on both the poloidal and toroidal planes (i.e. non-axisymmetric in toroidal plane as in the ballooning mode case [57,[61][62][63]). This is consistent with the observation of the fast reconnection time in ballooning mode case [i.e. 3-D nature] on TFTR [57]. Distortion of the 1/1 kink mode and localization in toroidal plane can be interpreted as a higher harmonic Fourier components of the 1/1 kink mode and there are numerous experimental reports on observation of harmonic modes of the 1/1 kink mode [65][66][67] as shown in Figure 11a-c. In the experiment, sawtooth crashes are dominated by the reconnection event of the highly distorted 1/1 kink mode with non-axisymmetry in the poloidal and toroidal planes and the corresponding crash time is an order of magnitude faster than τ c deduced from the model in a helically axisymmetric system. So, it is important to find a case of crash events that is close to the model where symmetry of the 1/1 kink mode is preserved so that the reconnection timescale in the helically axisymmetric case can be compared to that of the modeling.
The 2-D nature of the crash cases with slow reconnection time in TFTR discussed in earlier section is consistent with the full reconnection model which is based on 2-D modeling. These cases are rare compared to the crash cease with fast reconnection time. The 'post-cursor oscillation' is also a rare Figure 13. One sawtooth period without and with the current blip is compared. The sawtooth without current blip only has the 1/1 kink mode just before crash (upper left). One sawtooth cycle with the current blip in zone 3 in Figure 12, has the 3/3 mode in early rising phase after the crash. The 3/3 mode transforms into the 2/2 and 1/1 mode before the crash. The 2-D images of the 3/3, 2/2 and 1/1 mode are shown in the right side. The transformation is likely due to the change of the background q value as the current builds up with the increasing T e in the core. event of sawtooth oscillations and has been observed in many devices [57,59,60]. This oscillation is distinguished from the precursor oscillation. In precursor phase, the 1/1 kink mode is rapidly growing without reconnected field line and then crashes through the reconnection process. After the first reconnection event, the 1/1 kink mode in a reconnected state is slowly decaying in the post-cursor oscillation. In a toroidally rotating system, the slowly decaying 1/1 kink mode in a reconnected state (post-cursor) appears as an oscillatory motion and one period represents one full toroidal rotation of the plasma. The measured 2-D image of the post-cursor oscillation shows new information that is not available from conventional diagnostics as shown in Figure 14. Here, the remnant of the clear and symmetric 1/1 kink mode remains in a reconnected state for a long time after the first crash as shown in Figure 14a. The 1/1 kink mode is relaxing slowly as the T e rises. The continuously rotating reconnection zone in the view window due to the toroidal rotation suggests that the reconnection zone is helically symmetric around the torus and the reconnection timescale (i.e. end of rotation) is significantly slower than that of the dominant fast reconnection cases that we have experienced. One clear 2-D image shows a flow of the heat from the remaining 1/1 kink mode to outside of the q~1 layer along the field line through a wide reconnection zone (~7 cm) compared to narrow reconnection zone in the fast reconnection case is illustrated in Figure 14b. Note that the timescale of the change of T e (~μs order) is much slower than that of the electron thermal speed (~ns order) as we have discussed before. A reconstructed image of the reconnection process with the reconnected 1/1 kink mode at the high field side and cold island formation at the low field side is illustrated in Figure 14c. Through an extensive analysis of the post-cursor oscillation with 2-D images, the post-cursor case is consistent with the axisymmetric reconnection model used in the full reconnection model. The conclusion is that the crash time of the symmetric 1/1 kink mode in a helically axisymmetric reconnection zone along the toroidal direction (post-cursor case) is an order of magnitude slower than that of the more probable nonaxisymmetric crash cases. The helically axisymmetric case is rare due to the dominance of the nonlinear growth of the 1/1 kink mode bounded by q~1 resonance surface. Two distinctively different reconnection patterns are compared in Figure 15. In the axisymmetric reconnection process of the 'post-cursor' case in Figures 15b and 14b, the reconnection zone is long (~7 cm) along the poloidal direction compared to that of the 'ballooning like' nonaxisymmetric case (<2 cm) (Figure 15a) where the reconnection is likely induced at the tip of the 'finger' as shown in this figure. This pattern is consistent with the previous data from other devices ( Figure. 11a-c) including the TEXTOR case ( Figure 8). The slow crash timescale observed in the 'post-cursor' case (~8 ms) is consistent with τ c of the full reconnection mode whereas the crash time of the dominant fast crash case is much shorter~100μs. The common feature of the cases with the slow crash time is 2-D nature (helically axisymmetric system) while that of the cases with the fast crash time is 3-D nature (helically non-axisymmetric system). The two distinctive cases can be compared with the two well-known reconnection models: Sweet Parker and Petscheck models [48]. Physically the current sheet of the Sweet Parker model much longer than that of the Petscheck model as shown in the middle of these figures. Note that the Petscheck model was introduced to explain a fast reconnection time which cannot be explained with the Sweet Parker model.

Brief review of NTM instability
The neoclassical tearing mode (NTM) instability with low m/n mode numbers (m/n = 2/1, 3/1, 3/2, 4/1, 4/2, 4/3, … ., where m and n are poloidal and toroidal mode numbers, respectively) was first discovered on TFTR [76]. These classical internal instabilities not only have a potential of limiting the energy confinement of the core plasma but also can lead to disruption as the plasma beta (β) is further increased, where the β is the ratio of the plasma energy to the magnetic energy [77]. This instability is excited due to lack of helical bootstrap current inside the island (O-point) where the pressure profile is flattened near the resonant rational surfaces (i.e. q = 2,3,4, ...). A schematic illustration of the positions and shapes of these NTMs for KSTAR is depicted in Figure 16. The physics-based model of NTMs [78] indicates the importance of the trigger (or seed) to form resonant modes with islands. The 'seed' that starts the instability often comes from bursts of other MHD activity such as the sawtooth instability. While higher order NTMs are typically less problematic, the 3/1 NTM near the q~3 surface arises at modest plasma β and can limit the increase of the core plasma energy [78]. As the plasma β is further increased, the m/n = 2/1 mode can be excited at the q~2 surface. Onset of this instability imposes a significant limit to the growth of the core plasma energy due to altered cross-field transport in the presence of the island. Often rapid growth of the island above the critical size leads to a slowing of mode rotation in the plasma, with the mode magnetic perturbation eventually locking to the conducting vessel, causing a disruption. Therefore, control mechanisms to maintain the island size below a critical size is essential for high β tokamak operation. NTMs can be linearly stable and nonlinearly unstable and can be controlled with the injection of the current inside the island to compensate for the loss of the helical current inside the flattened island. While control of the 2/1 mode has been empirically successful with the injection of the current in the island using an ECCD system [79,80], the stability physics of the 2/1 mode has complex dependence on various plasma parameters and the cross-field transport has not been comprehensively validated yet. Among them, a comprehensive validation of the stability criteria of the 2/1 instability as well as the modified turbulence critical for the cross-field transport across the q~2 surface in the presence of the 2/1 island, has been attempted in the course of modeling studies. However, it has been a non-trivial problem to solve due to the lack of accurate diagnostics that can measure the dynamics of the 2/1 mode, the modified turbulence distribution and flow shear with adequate spatial and temporal resolutions.
The temporal change of the island size is expressed in terms of the classical linear stability index (Δ 0 ), which can be positive or negative and the destabilizing term from loss of the helical bootstrap current inside the island. The dynamics of the stability of the 2/1 mode has been described in the framework of the modified Rutherford equation (MRE) as given in where ω is the island half width, a 1 and a 2 are flux geometry-related coefficients, r s is the minor radius, τ r is the current diffusion time, τ s is the resistive time, ε is the inverse aspect ratio, β θ is the plasma poloidal beta, and L q and L p are scale lengths of the safety factor and pressure profile, respectively. It is critical to determine the sign of Δ' and ω c with high confidence in this equation, since the first term containing Δ' is closely related to the equilibrium current profile, and the second term is the destabilizing contribution stemming from the lack of the bootstrap current inside the island due to the flattened island. The ω is the half width of the island and relatively easy to determine in experiment. When the second term is significant (high poloidal beta and large island), the mode is referred to as a neoclassical tearing mode (NTM) and control of the NTM is critical for high-performance plasmas. Since precise determination of the critical parameters like Δ' and ω c is key, many experiments have been performed previously with conventional 1-D diagnostics such as ECE and Thomson scattering systems to aid in computing these quantities. However, the results were not conclusive due in part to the lack of spatiotemporal resolution of the data set. To overcome this deficiency of the data for accurate analysis, plasma jogging in a slow timescale was introduced assuming that the jogging does not influence the stability of the 2/1 mode [81]. It is extremely important to determine these two critical parameters (Δ' and ω c ) with confidence so that any deficiencies in modeling can be challenged for various cases for more effective control of the NTM. It is well known that the presence of the 2/1 mode (i.e. large magnetic islands) degrades (or limits) the plasma stored energy through increased cross-field energy transport [82]. The cross-field transport physics is complex, since the presence of the large island modifies the internal helical magnetic structure and creates fluctuations in 3-D magnetic flux and plasma pressure profile. Recent studies [83,84] indicate that multi-scale interactions between the large scale 2/1 mode and micro-turbulence can be a seed of the island formation and can influence the growth rate of the island non-linearly as well as the enhanced cross-field transport. Here, modified gradients of T e and n e profiles in presence of the 2/1 mode will subsequently modify the distribution of micro-turbulence as well as its flow pattern in the regions surrounding the island such as in the vicinity of O and X points [85,86]. The reduction of the micro-turbulence was measured in the flattened 2/1 island [86] and the importance of the modified flow shear across the 2/1 island was emphasized through measurement of micro-turbulence in the vicinity of the 2/1 island. Also, simulation studies suggest various multiscale interactions between the magnetic island and micro-turbulence via pressure and flow profiles [85]. They are the localized micro-turbulence distribution, poloidal vortex flow around the magnetic island [87], strong flow shear outside the 2/1 island [87-90] and insignificant micro-turbulence level across the O-point region due to small or flattened pressure gradient inside the magnetic island. The turbulence-induced transport is predicted to be significant in close proximity to the X-point [86,87,91,92]. Simultaneous measurements of the micro-turbulence and its flow in two-dimensional (2-D) space are ideal to fully validate the findings of numerical simulations.

What are the new findings in study of the 2/1 instability with the 2-D ECEi system?
For the validation of theoretical models and smart control of the 2/1 mode, clear and speedy identification of key parameters such as the stability parameter (Δ') and the critical island width (ω c ) of the 2/1 mode. The difficulty in determination of these parameters with confidence may have been largely due to lack of spatial and time resolutions using conventional 1-D data. Therefore, it is important to test if the confidence level can be improved in determining these parameters with 2-D data and the first such attempt was made in a study of the 2/1 mode with 2-D ECEi data on KSTAR [88]. In Figure 17a, the ECEi window is positioned to capture the 2/1 mode and the helical coordinate system (r, ζ) used in analysis is also indicated on the KSTAR plasma in this figure. The  1-D data is demonstrated for analysis of T e distribution inside the island and this is mainly due to high spatial resolution in the vertical coordinate (ζ) not in the radial coordinate (r) owing to the helical geometry and plasma rotation.
The detailed analysis can be found in Ref [88]. and summary of the results is provided here. The analysis of Δ' is relatively straightforward and is deduced from the calculated magnetic flux function based on magnetic field diffusion.
Since T e inside the island has a strong dependence on the parallel and perpendicular heat transport coefficients as well as the geometric factors of the island with critical width, synthetic T e data are assembled from the modeled T e profile [93] with dimensionless geometric factors and ratio of the parallel and perpendicular transport coefficients. And then statistical analysis between the measured and synthetic 2-D T e data is performed for optimum fits for each parameter sets. It is emphasized that higher spatial resolution is critical in resolving the synthetic T e data with different sets of parameters. As expected, the 2-D data cannot increase the resolving ability in the radial direction, but it can significantly improve the resolving ability in the helical axis utilizing data from other vertical channels and/or time-dependent data. In order to demonstrate a clear advantage of 2-D over 1-D data, seven channels of T e data along the mid-plane was assembled, and the same comparative study was repeated. As illustrated in Figure 18, the fitted T e data set for geometric factors and transport coefficients with the minimum χ 2 (Figure 18a,b) are used to estimate Δ' and ω c and statistical significance of the fits for sets of parameters in the 2-D modeled T e is significantly better compared to the fits with 1-D data in determination of the stability parameter (Δ') and critical island width (ω c ) as shown in Figure 18c.
Empirical understanding of the role of macroscopic MHD islands on limiting the core energy confinement as well as leading to disruption has to be supported by the comprehensive physical process of the cross-field transport physics based on modified micro-turbulence in presence of the 2/1 island. In recent years, high-frequency turbulence spectral measurement of T e employing a correlation technique has been routinely performed using 2-D ECEi data on KSTAR. Together with the measured macroscopic MHD fluctuation (i.e. 2/1 island) induced by the resonant magnetic perturbation (RMP), measurement of small-scale turbulences and the 2/1 island has provided an opportunity to study multiscale interaction between the macrofluctuation and micro-turbulence [94]. The simultaneously measured macroscopic fluctuation of the 2/1 island and dynamics of micro-turbulences in the proximity of the island in a discharge with R 0 = 180 cm, a 0 = 50 cm, and q 954 .6 are illustrated in Figure 19  This observation suggests that there are other mechanisms for this strong asymmetry in turbulence level, since the turbulence level should be higher at the steepened gradient region. On the outer side, the observed spectrum is narrow, and the level of turbulence is higher near the X-point. On the other hand, the measured flow speed of the micro-turbulence is striking as shown in Figure 19c. Here, the cross phase of two adjacent ECEI channels is calculated to measure 2-D flow velocity in the laboratory frame as shown in this figure.
The observed poloidal flow of the higher frequency spectrum is in the electron  Global gyrokinetic simulations were carried out to understand the underlying physics of the measured anisotropy of the micro-turbulence and strong flow shear in presence of the stationary 2/1 island [95]. The global equilibrium E x B flow perturbed by the 2/1 island was calculated using XGC1 code [96] as shown in Figure 20a. Here, the E x B flow is significantly enhanced in the presence of the O-and X-points of the 2/1 island. It was found that the flow shear is maximum and minimum near the O-and X-point of the island, respectively. From global micro-instability analyses using the gKPSP code [97], it was found that the collisionless-trapped electron mode (CTEM) and the ion temperature gradient (ITG) mode could be excited around the island as shown in Figure 20b. This is quite consistent with the measured high   portion of the CTEM is suppressed in the region near the O-point as shown in Figure 20c. On the other hand, the flow shear around the X-point is not strong enough to stabilize ambient instabilities, which is consistent with the experimental measurements shown in Figure 19c The Alfvén wave instability is a widely known phenomenon in laboratory, solar and space plasmas [98,99]. Alfvén waves in the plasma can be destabilized by the presence of resonating energetic particles. The energetic particles comes from various sources: accelerated ions and electrons from coronal mass ejection (CME) through magnetic reconnection processes in space and fast ions by heating systems [radio frequency (RF) and neutral beam injection (NBI)] in tokamaks. The Alfvén eigenmode instability [99] driven by fast ions are routinely observed in present tokamaks. This instability is particularly important for ITER and future fusion reactors, since the substantially increased fusion products such as 3.5 MeV α-particles in the core of plasmas can drive these instabilities and the transport of the α-particles can be significantly enhanced, thereby leading to losses of the heating source. A variety of mode structures such as the reverse shear Alfvén eigenmode (RSAE) [98] and toroidal Alfvén eigenmode (TAE) [98] have been extensively studied, and theoretical modeling is well advanced for tokamak plasmas. Full understanding of the physics of Alfvénic instabilities is important for advanced modes of operation where the current profile is likely reversed and energetic ions are abundant. Therefore, a detailed validation process is essential to avoid the deleterious instabilities and/ or develop preemptive control of them. In the validation process, measurement of the spatial extent and related dynamics of the eigenmode structure in the core of the tokamak plasma are extremely important. The 2-D ECEi system can be an ideal validation tool for these instabilities.

What are the new findings in the study of Alfvén mode instabilities with the 2-D ECEi system?
Studies of Alfvén eigenmode instabilities using 2-D ECEi systems were performed in both ASDEX-U, Germany [100] and DIII-D, USA [101] tokamak plasmas. Here, the plasma condition ideal for the RSAE and TAE is in the early stage of the discharge where the current profile is reversed, since penetration of the ohmically driven current is relatively slow. In order to destabilize the Alfvén eigenmode, early injection of the heating beam was introduced producing abundant fast ions in the core of the plasma. The fast-ion dynamics, which induce measurable modifications to the 2-D structure of RSAE and TAE instabilities were studied with the 2-D ECEi system. Since the electron temperature fluctuation amplitude of these eigenmodes is extremely small (~2.5 eV where the background plasma is~500 eV), signal enhancing techniques such as singular value decomposition (SVD) and Fourier decomposition were employed to eliminate uncorrelated noise.
As illustrated in Figure 21,22a, the measured 2-D images of RSAE structures at two different times (α and β) show two different mode numbers with distinctive spatial extents on ASDEX-U. Similar measurement was obtained in DIII-D, and an initial validation was carried out with the NOVA code [102] which is an ideal MHD code. The characteristics of the RSAE, such as the location, amplitude, frequency and mode number (n = 3) were in good agreement with the measured RSAE as shown in Figure 21,22b. Note that the maximum amplitude of this instability is located at the minimum q, -q min -, position of the reversed-shear q profile. With improved signal level and data analysis, the measurement revealed spiraling (poloidal shearing) of the tail of the mode in the ion diamagnetic direction as shown in Figure 20c. The NOVA code failed to reproduce this spiraling tail which is the poloidal shearing effect due to fast ions. Since the effect of fast ions is not included in NOVA, an ideal MHD code utilizing a perturbation approach, the non-perturbative gyro fluid MHD code, TAE/FL [103] was introduced. The simulation of the instabilities in TAE/FL at the corresponding toroidal mode number, n = 4 is compared with the experimentally obtained mode structures in Figure 21,22c. As shown in this figure, the radial localization and extent of the mode observed with the 2-D ECEi system is in excellent agreement with the simulation results including the spiraling tail. More simulation was carried out and each RSAE structure is dominated by a single poloidal mode number, m, which is positively identified in each case. The most distinguishing feature of the mode shown in this figure is a distinct poloidal shearing or 'spiraling' that was missed by the NOVA code. This mode shearing occurs as the regions at greater minor radius (observed on the outboard mid-plane) are shifted in the ion diamagnetic drift direction. Also, a discrepancy in mode frequency between the simulation and measurement arises due to the omission of compressibility in the simulation model. The difference is found to be consistent with the geodesic acoustic shift of the Alfvén continuum. This is an excellent example of how improved modeling improves the validation process along with the critical 2-D imaging data. The non-linear 'edge-localized mode (ELM)' instability is known as a periodic bursting event observed routinely in edge region of the H-mode plasmas (higher confinement mode) with the instability driven by high pedestal pressure gradient due to the formation of an edge transport barrier as schematically illustrated in Figure 23a. The transient heat and particles from partial collapse of the edge pedestal pressure are rapidly transported to open field regions. As the edge pressure gradient is reduced,  The H-mode plasma with high edge pedestal formed by the transport barrier is shown with the cross-section of the plasma with divertor where the heat from 'ELM' follows the magnetic field lines. (b) Due to burst of 'ELM', the edge pressure/current gradients are reduced and loss of the core plasma energy is followed during 'ELM' event. After 'ELM' event, pressure/current profile recover.
Source: Conner, J.W. www.ccfe.ac.uk/assets/Documents/AIPCONFPROC103p174.pdf the edge pedestal height is lowered and loss of the core plasma energy follows as shown in Figure 23b. The expelled transient high heat and particle flux can damage the divertor plates and first wall materials. Therefore, the 'ELM' in the H-mode plasma has drawn great attention and there are many research papers on both experimental and theoretical studies [104,105]. Especially, the impact of the transported transient heat and particles in large devices like ITER and comparable reactor devices could be detrimental to the first wall and divertor. Development of adequate control tools to mitigate, suppress and/or eliminate the 'ELM' instability has been actively pursued.
The 'ELM' event regulates the edge plasma pressure through repetitive non-linear bursting events and study of 'ELMs' has relied on D α light measurements in the divertor region as a measure of the 'ELM' onset and certain mode characteristics. Note that the rapidly transported plasma particles interact with neutral particles near the divertor plates, and D α lights are subsequently emitted. Various 'types' (i.e. type I-V, grassy, etc.) of 'ELMs' have been empirically categorized based on the frequency and amplitude of measured D α lights. In general, the type of 'ELM' changes from grassy (i.e. tiny bursts) to type-I (i.e. large burst with a low frequency) as the plasma energy is increased with increased auxiliary heating power. Variations in 'ELM' types are observed as the edge condition and plasma shape (i.e. triangularity and elongation) are changed [105]. For instance, Figure 23. Schematic presentation of the peeling-ballooning mode model: as the plasma edge pressure is increased with the heating power, the point is moving along the blue arrow (α-axis) until it reaches the ballooning boundary. Then slowly increasing current moves the point along the ballooning boundary (red line). As it hits the peeling-ballooning limit, crash occurs and loses edge pressure and returns to the initial position (green line).
Source: Conner, J.W. www.ccfe.ac.uk/assets/Documents/AIPCONFPROC103p174.pdf the type-I 'ELM' is characterized as a ballooning mode type due to the fact that the observation has been made mainly near the ballooning limit with no clear precursor prior the crash. On the other hand, type-III 'ELM' is shown to have a precursor and is classified as a resistive mode, since it has been observed in highly collisional plasmas. This difference between type-I and III can be compared to the sawtooth crash with precursor (tearing type) and without precursor (kink type).
The underlying physics of the 'ELM' has been suspected as a crash event of a rapidly growing ideal MHD mode with intermediate to high toroidal mode number (with a coupled poloidal mode spectrum) at the pedestal region of the plasma. A strong collisionality dependence of the 'ELM' dynamics may suggest that the nature of the 'ELM' may be a resistive rather than an ideal MHD mode. These high m/n modes are suspected to be driven by the steep pressure gradient and/or strong bootstrap current at the pedestal region of the plasma. The physical understanding of this instability has been largely through empirical experiences, since the nature of the eigenmode structure (e.g. spatial and temporal extent of clear m/n mode numbers) and corresponding dynamics have not been sufficiently diagnosed until recent years. Meantime, the 'peeling and ballooning' mode model [106,107], based on two driving energy sources for this instability have been developed: the edge pressure gradient and current density due to ohmic and bootstrap current at the edge, respectively. This model qualitatively has explained the stability boundary of the 'ELM' in a space spanned by the parallel current density and pressure gradient near the plasma edge as shown in Figure 24. As the edge plasma pressure builds up by the applied heating power, the point moves along α-axis (pressure) close to the critical value (α critical ) (path 1 in this figure). Slow increase of the current density brings the point along the ballooning limit to the peeling-ballooning limit (path 2 in this figure). When the point crosses to the unstable region, a crash occurs, with a partial loss of the pressure and current, and returns to the stable region where it started (path 3 in this figure). Validation of this model has been attempted with the ballooning and peeling boundary determined by the experimental data and modeled mode numbers with various assumptions. The extensive validation results of this model can be found in references [106,107]. Note that many critical parameters such as bootstrap current and pressure gradient are non-trivial to be measured and validated dynamics of the internal MHD modes are limited.
If the 'ELM' is based on ideal or resistive MHD modes, the modeling should be developed based on linear and non-linear dynamics of these modes through a step-by-step validation procedure. Then, building physics model of the transient crash and transport mechanism of the heat and particles to divertor is the next step. The advent of fast camera systems has provided an opportunity to study the fast phenomena like the 'ELM'. Here, fast-moving global filamentary structures such as created by the non-linear time evolution of the high m/n modes with a hint of fast transport event such as magnetic reconnection process [108][109][110] are captured by the camera in MAST as shown in Figure 25. Study of the 'ELM' based on filamentary mode structure was a significant step forward, since prior studies had been largely based on D α light spikes observed in the divertor region as discussed in the previous section. These did not provide visualization of the filamentary structures formed in the non-linear evolution period. Note that the captured filamentary structure is generally aligned along the helical magnetic field line and is also based on D α lights from the interaction between the radially moving filaments and neutrals at the scrape-off layer (SOL) region (outside of the separatrix). It is interesting to note that the filamentary structures are present in the inter-'ELM' period and L-mode and 'ELM' phase with the intensity increasing as the mode number is lowered from the inter-'ELM' period, L-mode to the 'ELM' phase as shown in Figure 25 [108]. Note that there was early observation of the filament structure on COMPASS-D, UK [111]. The COMPASS-D had low error field comparable to that of KSTAR. The dynamics of the filamentary structure with well-defined mode numbers are better physical quantity than D α lights from divertor region to build the physics basis of this instability. Ultimately, it is important to identify the internal mode structure and its linear and nonlinear behavior to properly implement MHD physics models. Here, the 2-D ECEi system can be instrumental, if the internal edge mode structures can be identified through measurement.
As noted before, large 'ELM' energy bursts in present research devices (reaching 6% of the total plasma stored energy) may not be detrimental but the transient heat and particles in an ITER-size device should not be underestimated. To reduce the risk to ITER operation, active mitigation, suppression or elimination of the 'ELM' without significant loss of confinement is essential. Many methods of mitigation or suppression have been pursued with external tools that can regulate the pedestal pressure such as the resonant magnetic perturbation (RMP) [112][113][114], pellets [115,116], magnetic control [117] and supersonic gas injection [118] together with control of the plasma separatrix strike points on the divertor plates, plasma shaping factors, etc. Essentially, all these external tools have been tested to regulate the edge pressure gradient/current so that the growth of unstable modes is maintained below the threshold level. Among external control tools, the RMP method was applied for control of the 'ELM' on the DIII-D tokamak for the first time and was demonstrated to be very effective in suppression and/or mitigation of the 'ELM' in a narrow plasma operation space [119]. Here, the conjecture was that the stochastic magnetic islands, induced by RMP coils with n = 3 or n = 4 toroidal modes in the pedestal region of the plasma, might enhance radial transport of the plasma pressure effectively reducing the edge pedestal pressure gradient so that the unstable mode is maintained below the threshold level. However, details of the RMP configuration for optimum penetration and the nature of the stochastic magnetic islands in the presence of the plasma are not clearly understood yet.
As ITER operation is near and design of future fusion reactors is presently based on H-mode operation, active research on control of the 'ELM' has been carried out using RMPs in various devices and results have been quite promising. Most of the perturbation experiments have been carried out with high toroidal mode numbers (n > 2) due to fear that deeper penetration of the perturbed field at lower toroidal mode numbers (n = 1 or 2) could provide a seed for destabilizing internal MHD instabilities such as the sawtooth and NTMs in the core of the plasmas. The three rows of in-vessel control coils (IVCC) with n = 1 and/or n = 2 toroidal mode numbers are used for suppression of the 'ELM' on KSTAR [120]. 'ELM'-free H-mode operation of a plasma with a pulse length over~30 s was produced [121]. However, the operation window of the 'ELM' suppressed regime is quite limited and efforts to expand the operating window have been largely empirical. Recently a successful suppression experiment with a pre-designed polar diagram, calculated with IPEC code including the plasma response model, spanned by the amplitude and phase of RMP coils has been conducted [122]. Here, one can select the amplitude and phase of RMP coils to operate in the region where suppression of the 'ELM' is expected for given RMP coil configurations for a specific helicity and set of physical parameters.
The first RMP experiment [112] introduced stochastic islands that may be responsible for suppression of the 'ELM' and a magnetic perturbation model with only the vacuum response may not be sufficient. Inclusion of the plasma in the model is essential. Therefore, it is important to measure the response of applied RMPs. In a recent DIII-D experiment, a substantial increase of long wavelength density turbulence was observed by the beam emission spectroscopy (BES) diagnostic during the 'ELM'-suppressed period under an applied RMP and the radial extent of the turbulence increased to r/a = 0.6 [123]. In another experiment, a sudden bifurcation of the perpendicular electron flow to near zero was suggested as a physical mechanism for suppression of the 'ELM' from a mitigated state using the n = 2 RMP configuration [124]. The role of RMPs in control of the 'ELM' can be divided in two stages: the first stage is preparing the optimum condition to couple in the magnetic perturbation induced by the IVCC configuration as demonstrated in construction of the polar diagram [122]. The second stage is to explore the mechanism from the 'ELMing' state to the mitigated or suppressed state through interaction between the filamentary mode structure and magnetic perturbation induced by the RMP. Therefore, it is important to clarify the relationship between the growth/decay of the filamentary structure, turbulence induced by RMP and enhanced particle transport near the pedestal (i.e. density pump out).
Up to now, the term 'ELM' represents the bursting event rather than the mode itself, because the study of this instability started based on measurements of transient D α lights near the divertor region assuming that the explosive growth of the edge-localized mode (ELM) is responsible without knowing the internal ELM dynamics. Then, the fast camera revealed the filamentary mode with clear mode structures. However, the images captured by the camera are based on D α light from outside of the separatrix where the radially displaced mode interacts with neutrals at the SOL region. The 2-D ECEi system is capable to measure the dynamics of the ELM near the pedestal region, and the term ELM represents the mode itself and the 'ELM-crash' will be used in remainder of the paper to mean the 'ELM' in order to avoid the confusion. The KSTAR plasma is ideal for studying the ELM dynamics and role of RMPs in suppression/mitigation of the ELM-crash, since KSTAR is equipped with three unique features and capabilities [121]; 1) a higher degree of symmetry of the tokamak plasma compared to other devices due to small intrinsic error field (approximately an order of magnitude lower error field compared to the tokamaks) and low magnetic ripple [125], 2) three rows of IVCC coils with the n = 1 and n = 2 toroidal mode numbers to use for the application of magnetic perturbations [122], 3) and an advanced 2-D ECEi system to study response of the magnetic perturbation [45]. As H-mode operation has been achieved on KSTAR, the 2-D ECEi system has been devoted to study instabilities at the plasma edge. As we discussed in the previous diagnostic section, interpretation of the ECE signals from the edge plasma may need special attention, since the optical thickness at the edge rapidly changes due to the steep density gradient in that region. Prior to the KSTAR experiment, 2-D ECEi system was applied to the edge physics study on ASDEX-U [126] and DIII-D [127].
In KSTAR experiments, a coherent mode structure near the pedestal region was detected with self-consistently determined toroidal mode numbers during the inter-ELM-crash period. An example of the ELM with n = 8 toroidal mode number is shown in Figure 26 [128]. Here, the 2-D image of the ELM and position of the ECEi window are marked on the KSTAR geometry (Figure 26a,b) together with D α light from the outer divertor region (Figure 26c) and spectrogram of one of the ECEi channels ( Figure  26d). The eigenmode structure of the n = 8 toroidal mode number is calculated with the resistive MHD code (BOUT++) using a reconstructed plasma equilibrium at the time when the mode was measured as shown in Figure 27a [128]. It is striking that the shape of the observed ELM in Figure 26b is quite different from the simulated mode structure in Figure  27a. In order to understand the apparent discrepancy, a synthetic image of the calculated eigenmode is constructed accounting for the intrinsic spatial resolution of each pixel as shown in Figure 27b. Here, even the downshifted images of the mode are shown in the far-right corner of this figure. Then, the background noise of the detection system is added as shown in Figure 27c. The final synthetic image in Figure 27c is directly compared with the measured one in Figure 27d. The discrepancy between the measured ELM radial width (~3 cm) and the simulated one (~1 cm) is understood based on the instrumental broadening of the present KSTAR ECEi system when the objects measured are smaller than the system resolution. Note that there is no loss of resolution in the poloidal direction, since the mode size is much larger than the vertical resolution (~2.5 cm).
Based on the clear validation of the ELM structure at the edge of the KSTAR plasma, the time evolution of the ELM cycle from growth to the crash has been studied and an example is illustrated in Figure 28 [129]. As shown in this figure, multiple modes are growing at the initial stage and the dominant mode is immediately saturated before the crash as it is rotating in the counterclockwise direction as depicted in Figure 28a. At the time of the crash, the images of multiple bursting processes of the filament are shown in Figure 28b,c. The bursting behavior resembles a local magnetic reconnection event that we observed in the non-axisymmetric reconnection event of the sawtooth instability (earlier section of this paper). Two examples of a single filament structure with a size less than~3 cm is evolving in time. As they rotate in the counter-clockwise direction, both filaments are slightly enlarged and elongated. The changes in size and elongation can be real, since the radial and vertical resolutions are~1.2 cm and~1.0 cm, respectively. The images are constructed based on the interpolation of each pixel due to the antenna pattern overlap between channels (~30%). Also, the optical thickness may not be changed as the background plasma pressure would be almost constant during~25μs. However, it is difficult to resolve the local reconnection event as we have seen in the sawtooth crash case due to a limited size of the filament, but a clear bursting of the filament is illustrated in these figures. It is also common to observe the ELM-crash with a toroidal mode number in the range n = 4~20 in KSTAR. In general, the toroidal mode number of the ELM filament is changing from high to low (e.g. n = 11 to n = 4) as the amplitude of D α spike is increased. Qualitatively, the observed mode structure is consistent with the filamentary structure from the fast camera illustrated in Figure 25. Another important transient event is observation of the solitary perturbation [130] before the ELM-crash event. This mode is a localized filament with a low mode structure (i.e. n = 1). The solitary perturbation appears together with a quasi-steady mode (QSM) in the inter ELM-crash (n~7), and the rotation direction is opposite each other. This solitary perturbation has been rarely observed, and it is premature to connect this event to the crash physics of the ELM at this point. If the physics of the ELM and ELM-crash is based on the peeling-ballooning mode model, it would be expected to find that the mode strength as well as the localized bursting would occur on the low field side. While a majority of experiments focused on measuring the low field side of the plasma, two Figure 30. The ELM-crash suppression experiment with a slow RMP current ramp up and down time scale (green) comparable to the magnetic diffusion time scale (~2s) for the KSTAR edge plasma parameters. (a) The intensity of D α spikes is linearly reduced as the amplitude of the broadband turbulence is linearly increased with the RMP current ramp up while the perpendicular flow speed of the turbulence suddenly dropped to minimum at about the same time when the ELM-crash is suppressed. The decay of the turbulence level is significantly delayed as the RMP ramp down is started and D α spikes returned when the amplitude of turbulence dropped to the level where the suppression was started and the perpendicular rotation is increased suddenly. (b) The coherence spectra for poloidal grows as the RMP is ramped up and radial spread of the turbulence occurs when the turbulence level is fully saturated poloidally. Source: Park, H imaging stations were simultaneously deployed to both high and low field sides at the edge of the plasma to study the ELM dynamics on the same flux surface as illustrated in Figure 29. The simultaneously observed 2-D images of the ELM on the high and low field side during inter ELM-crash period are depicted in this figure and the timing of this measurement is~0.5 ms before the crash. There are three distinctive observations in this study: First, the observed mode intensity at both the high and low field side is comparable. This is inconsistent with the simulated results from MHD codes such as the BOUT++ code and M3D-C 1 in which the simulated intensity of the mode at high field is an order of magnitude weaker than that at the low field side due to the ballooning nature of the ELMs. However, the recent numerical simulation from the kinetic code, XGC, produced a comparable ELM mode strength at both high and low field side [131]. Further study is needed to clarify the different results between the MHD and kinetic simulation codes. On the other hand, this result is striking but it is not surprising, since the classical ballooning mode model for the sawtooth instability failed to explain the observation of the crash on the high field side as reported in sawtooth instability section. Note that difference of the magnetic pressure between high and low field side in the sawtooth case is much smaller than in the ELM-crash case due to the dependence of the magnetic field strength on major radius. Second, the observed poloidal wavenumber of the ELM at high field side is almost twice lower than the expected poloidal wavenumber on the same flux surface corresponding to that of the low field side as shown in white lines in this figure. This may suggest that the eigenmode at the low field side can be different from that of the high field side due to the different boundary conditions. Third, the observed rotation direction was opposite each other with different speed. The mode is rotating in clockwise direction on the low field side while it is rotating in counter-clockwise direction on the high field side. This non-intuitive observation may come from an asymmetry of E r at the high and low field sides, since the observed rotation is a combination of E r x B, drift velocity, drift velocity and phase velocity of the mode. So far, the physics basis of the ELM is still in an infant stage, and dynamics and global behaviors of the observed eigenmode in the ELM cycle will provide an opportunity to establish an improved basis of the ELM physics.
Without a fundamental physics basis of the non-linear ELM dynamics, understanding of RMP dynamics in suppression of the ELM-crash may be premature. Empirically the RMP has been very effective in suppression and/or mitigation of the ELM-crash as described in the previous section. For instance, predictive design of the operation window for the ELM-crash suppression based on the amplitude and phase of RMP coils [120] is quite encouraging but the recipe is not universal yet. The essence of this design suggests that the modeling has to account for the presence of the plasma in modeling of the RMP-induced suppression. In general, when the ELMcrash is mitigated by the RMP, edge density pump-out is followed and the pump-out is less from the mitigation to the suppression phase. At this stage, it is important to clarify the role of turbulence in the growth/decay of the ELM amplitude and enhanced transport of the edge density so that an appropriate transport physics model can be applied to minimize the loss of edge pressure due to density pump out while the suppression period is sustained. If the observed turbulence is in response to the RMP in presence of plasmas, non-linear interaction between the ELM and turbulence can slow down growth of the ELM below the threshold level. At the same time, turbulence can enhance the radial transport of plasma density. Note that a wide spread of the density turbulence, while the RMP is active, has been reported [123] as discussed in an earlier section. However, identification of the observed turbulence is challenging, since the characteristics of the measured turbulence of n e , and T e , and magnetic turbulence is quite similar to each other in the wavelength and frequency space.
In an early experiment on the ELM-crash suppression using RMPs on KSTAR [132], as the RMP was applied, the ELM structure was not changed immediately but the change in toroidal mode number (from n~15 to n2 0) is notable with marginally stable phase as shown in the last image in Figure 30a. The time history of the sum of the broadband turbulence (~30-70 kHz) and the coexisting eigenmode (ELM) (~5-30 kHz) is plotted as the RMP current is ramped up as shown in Figure 30b. It is interesting to note that the RMP is ramped up with a fast ramp-up timescale (~0.2 s) but the sum of the broadband turbulence in response to the RMP is slowly increasing while the ELM (not ELM-crash) amplitude decreases on a similar timescale as illustrated in this figure. This observation may suggest that penetration of the magnetic perturbation from the applied RMP is slower than the particle transport timescale. If the slowly increasing turbulence level drains the energy of the ELM eigenmode through cross interaction between the turbulence induced by RMP and the ELM, the timescale of the decreasing D a spike amplitude should occur in a similar timescale but the suppression occurs as the fast ramp-up is completed. The fact that the toroidal mode number of the inter-ELM period is changed from n~15 to n~20, implies that the ELM has been marginally unstable. The growth of the ELM may have been slightly above the threshold level and small increase of the turbulence could drain the energy of the ELM and suppress the ELM-crash.
In a recent experiment for further study on the role of the turbulence as a suppression mechanism of the ELM-crash, the RMP current is ramped up and down slowly in a timescale comparable to the magnetic diffusion timescale (τ md~2 s) for the edge plasma parameters as shown in Figure  31. The amplitude of the sum of the turbulence spectra from the inter-ELM-crash period is increased almost linearly with the RMP current while the amplitude of D α light spikes (i.e. ELM-crash) is decreasing linearly until the suppression time. The perpendicular electron flow speed is suddenly dropped when the ELM-crash is suppressed. The sudden bifurcation of the perpendicular electron flow is consistent with the previous observation in Reference [124] where magnetic signals have been used to address the electron flow. Note that direct measurement of turbulence is used to address these issues in KSTAR for the first time. Since the poloidal flow speed is likely dominated by the E r x B term, it is logical to conclude that the sudden change of the perpendicular flow speed has to do with the change of E r . It is important to separate this effect from the E r x B physics of the H-mode near the pedestal region which is extremely narrow [133]. As the RMP current is ramped down, amplitude of the turbulence has a long-time delay before it starts decreasing. At about the same amplitude of the turbulence when the ELM-crash was suppressed, the ELM-crash (i.e. D α light) and perpendicular flow speed return at much higher intensity of D α light. Does this implicate recovery of the E r x B? If the sudden bifurcation of the electron flow speed is associated with the E r change [133], it is reasonable to conclude that the time history of the turbulence is not correlated with that of E r . As shown in Figure 31b, the time evolution of the observed turbulence on poloidal plane shows that the level is increasing with broad spectra as it approaches to the flat top of the RMP current ramp-up phase. On the other hand, the radial spread of the turbulence is only shown as the poloidal turbulence is saturated. This observation suggests that the turbulence induced by the RMP is radially very narrow (~5 cm) and the spread only occurs when the penetration of the magnetic perturbation is fully saturated. More experiment should be conducted to clarify these issues such as the coupling efficiency of the magnetic perturbation, role of the turbulence in enhanced transport of the electron density and nonlinear interaction with the ELMs.
The introduction of the 2-D visualization of the ELM structure with a wide range of toroidal mode numbers (n = 4~n = 20) and its dynamics at the pedestal region, has provided abundant new information such as mode validation through synthetic imaging of the calculated eigenmode, observation of the fast reconnection crash process, and the role of turbulence in suppression of the ELM-crashes using RMPs. The most striking observation is from the RMP-induced suppression experiment with short and long timescale up to the magnetic diffusion timescale. While the suppression state can be clearly defined when there is no ELM-crash, the gradual decrease of the ELM-crash amplitude before the suppression raises an issue on defining the mitigation state. The fact that the timescale of the observed turbulence is quite different from the perpendicular flow speed suggests that the relationship between the broadband turbulence, E r physics and enhanced transport may need more experiments to be resolved. Any phenomenology based on 'a posteriori' measurement such as D α spikes and filamentary structure by the fast camera is subject to individual device conditions and the experimental set-up. It is desirable to build physics model of the ELMs, ELM-crashes and control dynamics based on universal physics parameters such as the internal mode dynamics. Eventually, the internal dynamics of the ELMs will be connected to the global filament structures by the fast camera and D α spikes in the divertor.

Relevancy of the new knowledge of MHD events to non-fusion plasmas
In general, significant progress has been made in understanding linear to non-linear evolution of MHD physics in the laboratory [4,5] and solar plasmas [2,3] in which boundary conditions and constraints are vastly different. The naturally occurring non-linear MHD events that we have observed in tokamak plasmas are more complex than present symmetrical physical model based. Similar events in solar and space plasmas may have more variations due to freer boundary conditions. Therefore, it is not surprising that the validation process using measurements from conventional diagnostic systems may have experienced difficulties and/or limitations, since the limited aspects of the measurement may overlook asymmetries and anisotropies while the theoretical models are developed based on symmetric and isotropic systems.
The first example can be the reconnection physics of the sawtooth instability, and driven reconnection experiments in laboratory [4,5] may have more relevance to the theoretical models than the naturally occurring reconnection events in tokamaks and solar plasmas. The dynamics of the magnetic flux rope of a solar flare and associated corona mass ejection (CME) [134] through reconnection events have a strong resemblance to the reconnection event of the 1/1 kink instability in a sawtooth oscillation. A major difference could be that the 1/1 kink instability in a sawtooth oscillation is bounded by the q = 1 surface whereas the flux rope in solar flare (i.e. half loop as shown in Figure 32) has no boundary and kinking of the current-carrying flux rope can be freely twisted. Therefore, the reconnection process can have more freedom than that of the 1/1 kink mode in tokamak plasmas [135,136]. The common factor is that kink modes in a solar flare and a sawtooth oscillation in a tokamak carry current and its shear will be important factor for growth of the instability as shown in simulation of the flux rope in solar flare [137]. An example of modeling using an ideal MHD code with zero beta plasmas is given in Figure 32.
Here, it is assumed that the initial magnetic rope (i.e. half of toroidal plasma) is created by strong current under the chromosphere and the rope carries a current which is the source of the kink instability similar to the case of sawtooth kink instability.
The second example can be the reconnection timescale which appears to be different whether the reconnection event is in a symmetric system or asymmetric system. The dominant fast reconnection time in sawtooth crash event has been observed mainly in a localized reconnection zone (non-axisymmetric system), and the initiation was proceeded at the pointy region in a ballooning type bulge (i.e. highly distorted 1/1 kink mode) as demonstrated in 2-D images prior to the crash. It is not surprising to find that the characteristic reconnection time (τ c ) estimated based on the Sweet-Parker model in a helically symmetric system does not agree with the measured timescale. On the other hand, this characteristic time (τ c ) is consistent with that of the axisymmetric reconnection zone ('post-cursor oscillation' case). Therefore, the Petschek model may be more suitable for the fast reconnection with a narrow reconnection zone and thin current sheet. The reconnection events in the magnetic rope of solar plasma can have various twists of the kink due to various possible current flow patterns in open boundary, and the timescale can be vastly different. The lessons from reconnection physics in the study of the sawtooth instability with a 2-D visualization tool suggests that application of the theoretical model should be done with caution, since little is known about the currents in the flux rope and boundary conditions such as the axisymmetry.

Summary
Among various concepts of fusion plasma devices, the tokamak concept has been at the front of the line and the progress on the stellarator concept is rapidly catching up. The future success of the ITER project warrants challenging a design of advanced fusion reactor and modeling with predictive capability is essential. New insights and multi-dimensional boundary constraints are critical in validating the theoretical models for improved control system for the critical MHD and Alfvén instabilities. New findings of these instabilities by 2-D ECEi system have not only provided an opportunity to improve the relevant theoretical modelings but also a broadened understanding of the underlying physical mechanisms. Validation of the long-standing q 0 problem (i.e. MHD aspect of the 1/1 kink stability condition) and new insights of the reconnection timescale (i.e. the importance of axisymmetry of the 1/1 kink mode in determining the reconnection time) in the sawtooth instability vindicated the discarded 'full reconnection model'. Furthermore, the newly uncovered physics may find applications to similar instabilities and reconnection events in solar and geomagnetic plasmas. The modeling capability of the Alfvén eigenmode instability is improved owing to detailed 2-D images of the spatially distributed RSAE. In study of the NTM instability, it is demonstrated that 2-D imaging data can be effective in determining the critical stability parameters of the 2/1 mode (i.e. Δ' and ω c ). Simultaneous measurement of macroscopic MHD (i.e. 2/1 island) and the non-intuitive anisotropic micro-turbulence and its flow shear in the presence of a 2/1 island is justified by theoretical models. Here, the anisotropy of microturbulence and strong flow shear are from a competition between growth rates of the TEM and ITG and E x B shearing rates near the X and Opoints of the 2/1 island. The 2-D imaging data of the ELMs and ELM-crash has provided unprecedented new physics information. Here, validation of the mode structure with the reduced MHD code (M3D-C 1 ) provided a confidence of the measured ELM structures. The 2-D dynamics of the ELMs such as temporal evolution to the crash and asymmetry of mode behaviors are new information of the ELMs. In the ELM-crash suppression experiment, the simultaneously measured ELMs and turbulence induced by the RMP provided an opportunity to understand fundamentals of the ELM-crash suppression physics. The slow growth of the turbulence on the magnetic diffusion timescale with the slow RMP current ramp-up suggests the nature of the turbulence may be magnetic rather than electrostatic. Also, bifurcation of the perpendicular electron flow was observed through direct measurement of the turbulence. No correlation between the growth of the turbulence and flow speed raises an important issue on the E r x B contribution to the perpendicular flow speed and enhanced transport at the pedestal region. The inter-coupling of the growth and decay of the turbulence, the timescales of the ELM-crash amplitude (D α spikes) into and from the suppression period and perpendicular rotation are extremely important for the suppression mechanism of the ELM-crashes.