We analyzed the higher harmonics of temperature perturbation in the modulated electron cyclotron heating experiment on TJ-II plasma. The higher harmonics (e.g., 5th harmonic) exhibited significantly weaker decay in amplitude as they propagated in radius as compared with the prediction by diffusive model. The change in the time derivative of temperature at the onset (and turning-off) of the heating power propagates in radius with very little temporal smoothening. Recently, a new method for analyzing anomalous transport in magnetically confined plasmas was developed. This new method directly measures the gradient-flux relation by studying the temperature response with respect to modulations of the electron cyclotron heating (ECH) power. Using this method, hysteresis in the gradient-flux relation was observed in LHD plasmas [1,2]. Because heat flux is a multivalued function of gradient, the temperature perturbation dynamics are far from a simple diffusive response with constant diffusion coefficient. Therefore, higher harmonics can exhibit smaller or larger spatial damping rates in comparison with the diffusive response. To observe the violation of the simple diffusive relation q = −nχ∇T with homogeneous χ, measuring the higher harmonics of temperature perturbation in the heating modulation experiment was proposed [3]. In this rapid communication, we report the application of the abovementioned method to the modulated electron cyclotron heating (MECH) experiment on the TJ-II [4].The power modulation experiment with line-averaged values around 0.5 × 10 19 m −3 was performed on a lowdensity plasma using continuous on-axis heating with 250 kW and modulated 50 kW (square waveform, 50% duty cycle) microwave beams. The modulation frequency of ECH ( f mod ) is 180 Hz, mostly absorbed in the central core r/a ≤ 0.3 (plasma radius a = 0.2 m), as determined in previous heat wave propagation experiments [5]. The electron temperature profile was measured at ten points in the high-field side region of 0.0 < r/a < 0.9 using an ECE author's e-mail: inagaki@riam.kyushu-u.ac.jp heterodyne radiometer with a time resolution of 10 µs [6]. The comparison with the T e measured using a He beam indicates that ECE can measure T e even at the edge [7]. We applied the conditional averaging technique [1]. In this procedure, the ECH turn-on time (t on ) is detected in each period and temporal evolutions of the electron temperature for time intervals of −3 ms < t − t on < 3 ms were extracted; the extracted signals were averaged over 25 modulations. Post averaging, the noise level significantly decreased. Figure 1 illustrates the power spectrum of the conditionally averaged temperature perturbation and Fig. 1 Power spectrum of the conditionally averaged temperature perturbation (black solid line) and that calculated from the raw signal (red dashed line). Higher harmonics up to the 5th are presented. The black dashed line is for visual aid.