The global build-up to intrinsic edge localized mode bursts seen in divertor full flux loops in JET

A global signature of the build-up to an intrinsic edge localized mode (ELM) is found in the temporal analytic phase of signals measured in full flux azimuthal loops in the divertor region of JET. Toroidally integrating, full flux loop signals provide a global measurement proportional to the voltage induced by changes in poloidal magnetic flux; they are electromagnetically induced by the dynamics of spatially integrated current density. We perform direct time-domain analysis of the high time-resolution full flux loop signals VLD2 and VLD3. We analyze plasmas where a steady H-mode is sustained over several seconds during which all the observed ELMs are intrinsic; there is no deliberate intent to pace the ELMing process by external means. ELM occurrence times are determined from the Be II emission at the divertor. We previously [Chapman et al., Phys. Plasmas 21, 062302 (2014); Chapman et al., in 41st EPS Conference on Plasma Physics, Europhysics Conference Abstracts (European Physical Society, 2014), Vol. 38F, ISBN 2-914771-90-8] found that the occurrence times of intrinsic ELMs correlate with specific temporal analytic phases of the VLD2 and VLD3 signals. Here, we investigate how the VLD2 and VLD3 temporal analytic phases vary with time in advance of the ELM occurrence time. We identify a build-up to the ELM in which the VLD2 and VLD3 signals progressively align to the temporal analytic phase at which ELMs preferentially occur, on a ∼2−5ms timescale. At the same time, the VLD2 and VLD3 signals become temporally phase synchronized with each other, consistent with the emergence of coherent global dynamics in the integrated current density. In a plasma that remains close to a global magnetic equilibrium, this can reflect bulk displacement or motion of the plasma. This build-up signature to an intrinsic ELM can be extracted from a time interval of data that does not extend beyond the ELM occurrence time, so that these full flux loop signals could assist in ELM prediction or mitigation.

The global build-up to intrinsic edge localized mode bursts seen in divertor full flux loops in JET A global signature of the build-up to an intrinsic edge localized mode (ELM) is found in the temporal analytic phase of signals measured in full flux azimuthal loops in the divertor region of JET. Toroidally integrating, full flux loop signals provide a global measurement proportional to the voltage induced by changes in poloidal magnetic flux; they are electromagnetically induced by the dynamics of spatially integrated current density. We perform direct time-domain analysis of the high time-resolution full flux loop signals VLD2 and VLD3. We analyze plasmas where a steady Hmode is sustained over several seconds during which all the observed ELMs are intrinsic; there is no deliberate intent to pace the ELMing process by external means. ELM occurrence times are determined from the Be II emission at the divertor. We previously [Chapman et ISBN 2-914771-90-8] found that the occurrence times of intrinsic ELMs correlate with specific temporal analytic phases of the VLD2 and VLD3 signals. Here, we investigate how the VLD2 and VLD3 temporal analytic phases vary with time in advance of the ELM occurrence time. We identify a build-up to the ELM in which the VLD2 and VLD3 signals progressively align to the temporal analytic phase at which ELMs preferentially occur, on a $2 À 5ms timescale. At the same time, the VLD2 and VLD3 signals become temporally phase synchronized with each other, consistent with the emergence of coherent global dynamics in the integrated current density. In a plasma that remains close to a global magnetic equilibrium, this can reflect bulk displacement or motion of the plasma. This build-up signature to an intrinsic ELM can be extracted from a time interval of data that does not extend beyond the ELM occurrence time, so that these full flux loop signals could assist in ELM prediction or mitigation. Enhanced confinement (H-mode) regimes in tokamak plasmas are characterized by intense, short duration relaxation events known as edge localized modes (ELMs). [1][2][3][4][5] Prevention of large amplitude ELMs is essential for ITER as each ELM releases particles and energy which load the plasma facing components; scaled up to ITER, 6 the largest such loads would be unacceptable. Theoretical 7,8 and observational 9 works suggest that the peeling-ballooning MHD instability of the plasma edge underlies ELM initiation, but as yet there is no comprehensive understanding of the sequence of physical processes involved in ELMing in terms of self-consistent nonlinear plasma physics.
Quantitative characterization of the dynamics of ELMing processes via their time domain properties, such as inter-ELM time intervals and ELM event waiting times, is relatively novel [10][11][12][13][14][15] and has provided evidence of unexpected structure in the sequence of ELM occurrence times. Recently, 16,17 we found that the signals from a system scale diagnostic, the toroidally integrating full flux loops in the divertor region of JET, contain statistically significant information on the occurrence times of intrinsic ELMs: the ELMs tend to preferentially occur when the full flux loop signals are at a specific analytic phase of their timeseries. Since the full flux loop signals capture aspects of the global plasma dynamics including large scale plasma motion, this may suggest, as first proposed in Ref. 17, a nonlinear feedback on a global scale where the control system and plasma behave as a single nonlinearly coupled system, rather than as driver and response. This feedback may act to pace the intrinsic ELMs. a) Electronic mail: S.C.Chapman@warwick.ac.uk b) In this paper, we investigate the time dynamics of the full flux loop signal temporal analytic phases, and directly test whether these signal phases contain information on the build-up to an intrinsic ELM. We perform direct time domain analysis of high time resolution signals from the full flux loops in the divertor region in JET. These full flux loop VLD2 and VLD3 currents are proportional to the voltage induced by changes in poloidal magnetic flux. We use a simultaneous high time resolution Be II signal to determine the intrinsic ELM occurrence times. We focus on a sequence of JET plasmas that have a steady flat top for $5s and which all exhibit intrinsic ELMing in that there is no deliberate intent to control the ELMing process by external means. Importantly, the full flux loop signals have sufficiently large signal dynamic range, compared to the noise, to allow the time evolving instantaneous temporal analytic phase to be determined on timescales between one ELM and the next. ELMs tend to occur preferentially at a specific temporal analytic phase in the VLD2 and VLD3 signals. Here, we find that the temporal analytic phases become progressively more strongly ordered from about 2-5 ms before the ELM up to the ELM time. Furthermore, the VLD2 and VLD3 signals become temporally phase synchronized with each other during this build-up time. Global synchronized plasma dynamics is thus part of the build-up to an intrinsic ELM. The organization of the paper is as follows: in Section II we introduce the data used in this study, in Section III we describe how the full flux loop instantaneous temporal analytic phases are determined, our main results are given in Section IV, in Section V we quantify the statistical significance of these results, and in Section VI we present a possible interpretation of these results following the scenario of Ref. 17. We provide significance tests against null hypotheses, that is, phase alignment by chance coincidence, in the Appendix.

II. ELM AND FULL FLUX LOOP TIME SIGNATURES
We analysed the sequence of JET plasmas 83769-83775 discussed in Ref. 16. These are a subset of plasmas 83630-83794 analysed in Ref. 14. Each has a flat-top H-mode duration of $5s. These all exhibit intrinsic ELMing in that there is no attempt to precipitate ELMs; the only externally applied time varying fields are those produced by the control system. ELM occurrence times are inferred from the Be II signal, which we will compare with measurements of the inductive voltage in the full flux loops VLD2 and VLD3. These circle the JET tokamak toroidally at a location just below and outside the divertor coils. These signals are sampled at 100 ls time resolution, which can be compared to the shear Alfven transit time at the top of the edge pedestal which is of order a few microseconds.
The configuration of these diagnostics on JET is shown in Figure 1. The signal voltage is induced by changes in poloidal magnetic flux through the surface encompassed by the loops. We determine the ELM occurrence times t k by identifying the peak of the Be II signal within each ELM using the method described in Ref. 16.
From the occurrence times t k of these peaks, the time intervals between successive ELMs, the ELM waiting times, Dt k ¼ t k À t kÀ1 are found. In these plasmas, there is time structure in the probability density of ELM waiting times. There is a lower cutoff in the ELM waiting time at Dt $ 10ms, and there are time intervals where ELMs occur less often (Ref. 16, and in other plasmas 12 ). Large ensemble statistical studies across many JET plasmas have also revealed 16,17 that the ELM waiting time probability distribution shows time structure, that is, some ELM waiting times are more likely than others.
Signal traces for a representative pair of successive ELMs with a waiting time of $30ms are shown in Figure 2 for plasma 83771. The absolute polarity of these signals is set by convention. Following each ELM, the figures show a characteristic large amplitude oscillatory response in both of the full flux loop signals, the first cycle of which is on a timescale of $10ms. We previously identified 16 a class of prompt ELMs which are clustered approximately within 10 < Dt < 15ms and appear to be directly paced by this response to the previous ELM. These prompt ELMs will be excluded from the current analysis, here we consider ELMs that occur on longer timescales such that this initial flux loop signal response to an ELM is seen to decay. Intervals of quasiperiodic oscillations can be seen in the VLD2 and VLD3 signals throughout the time between one ELM and the next. We will now directly obtain the instantaneous temporal analytic phase of these signals in order to test for information in these oscillations.

III. DETERMINATION OF FULL FLUX LOOP INSTANTANEOUS TEMPORAL ANALYTIC PHASE
A time series S(t) has a corresponding analytic signal defined by SðtÞ þ iHðtÞ ¼ A exp½i/ðtÞ, where H(t) is the Hilbert transform of S(t), defined in Refs. 18, 19, and 21, see also Refs. 20 and 22. This defines an instantaneous temporal analytic amplitude A(t) and phase /ðtÞ ¼ xðtÞt where the instantaneous frequency is xðtÞ for the real signal S(t). We compute the analytic signal by Hilbert transform over each waiting time Dt k between each pair of ELMs. The procedure is summarized in the schematic shown in Figure 3, which shows the domain over which the Hilbert transform is calculated relative to a pair of ELMs occurring at t kÀ1 ¼ t ELM1 and t k ¼ t ELM2 . We will obtain the temporal analytic phase for the full flux loop signals for a sequence of times dt preceding the second ELM of each pair, that is, at times t ELM2 À dt. We will need to choose a zero time t 0 to define a temporal analytic phase difference in the full flux loop signals D/ ¼ /ðtÞ À /ðt 0 Þ; here t 0 will be an estimate of the occurrence time of the first ELM.
The full flux loop signals are sufficiently above the noise that we can use this method to determine their instantaneous temporal analytic phase. The instantaneous temporal analytic phase cannot be directly extracted for the Be II signal because its noise level is usually too high. We first perform a 3 point spline smoothing on the VLD2 and VLD3 time series to remove noise fluctuations on the sampling timescale. The signal analyzed must oscillate about zero in order for the instantaneous temporal analytic phase to be well determined from the analytic signal, we can ensure this locally by subtracting a locally determined mean specified as shown in Figure 2. The signal local mean is determined within a window that is shifted back in time by dt, in the results shown here we used a window T A ¼ ½t k À 0:01; t k À 0:0025s À dt, relative to each ELM occurrence time t k (so that on the schematic dt m ¼ 2:5ms). The Hilbert transform has a singlesided Fourier transform which is approximated via fast Fourier transform over the finite time window of the data. We therefore need to use a time window that is larger than that of the time domain of interest to avoid edge effects, so that we only calculate the instantaneous temporal analytic phase at times within a window edge time interval dt E of the ends of the time window of data. We have found that in these time-series a dt E > 1ms is sufficient to give stable values of the instantaneous temporal analytic phase and all results presented here use this value of dt E . We have varied T A , dt m , and dt E to check the robustness of our results.
The above methods are only effective if the full flux loop signals have good signal/noise, do not have too large a dynamic range in response to all the ELMs, and if the mean of the signal does not vary too rapidly. The high rate of change of instantaneous temporal analytic phase with time of the full flux loop signals requires well defined ELM occurrence times in order to cleanly determine any phase relationship.

IV. FULL FLUX LOOP INSTANTANEOUS TEMPORAL ANALYTIC PHASE AND BUILD-UP TO AN ELM
We will first discuss results using Hilbert transform window (a) which extends beyond the time t ELM2 of the second ELM, so that we can obtain the instantaneous temporal analytic phase at times up to the ELM occurrence time. In Figure 4, we plot the instantaneous temporal analytic phase of the full flux loop signal versus time for all the ELMs in JET plasma 83771. The main figure panel plots time from t 0 , that is, Dt ¼ t À t 0 versus the instantaneous temporal analytic phase difference D/ ¼ /ðtÞ À /ðt 0 Þ of the VLD2 signal. In Figure 4, we set t 0 ¼ t ELM1 the time of the first ELM as determined from the Be II signal. The first (red circle) and second To facilitate comparison, we have standardized the signal amplitudes by dividing by a multiple (10 for Be II, 2 for the VLD2 and VLD3) of their respective means over the flat-top H-mode duration, and then subtracted a local mean calculated over the interval denoted by the pair of vertical dotted-dashed blue lines. The sign convention of the VLD2 and VLD3 signals is chosen, such that they have opposite polarity. The ELM occurrence times are indicated by vertical red and green lines. For reference the time interval between 0 and 5ms before the second ELM is shaded in grey.
FIG. 3. Schematic (not to scale) of the procedure to determine the time dependent full flux loop difference in temporal analytic phase D/ k ðdtÞ as a function of the time interval dt measured back from the time t ELM2 of the second ELM in an ELM pair. D/ k ðdtÞ is with respect to the temporal analytic phase at the time t ELM1 of the first ELM. The temporal analytic phase can only be determined at times that are within the Hilbert transform window edge dt E . One can then choose a Hilbert transform window (a) that goes beyond t ELM2 , so that dt > 0 or (b) that stops at t ELM2 , so that dt > dt E .
(green circle) ELM times, as determined from the Be II signal, are overplotted on each corresponding VLD2 trace. On the plot, the first ELM has coordinates Dt ¼ 0 and D/ ¼ 0 by definition. The coordinates of the second ELM are the waiting time Dt k ¼ t ELM2 À t ELM1 and corresponding temporal analytic phase difference D/ k ¼ /ðt ELM2 Þ À /ðt ELM1 Þ. Histograms are shown of the waiting times Dt k (top panel) and differences in temporal analytic phase D/ k (right panel) for all the k ¼ 1; :::; N ELM pairs. There is a group of prompt ELMs 16 with Dt < 15ms, indicated by pink bars, which are distinct in both arrival time and temporal analytic phase. We have previously identified these prompt ELM events as being directly correlated with the response to the previous ELM and will exclude them from the analysis to follow by only considering ELM pairs with waiting times Dt > 15ms. These ELMs, with Dt > 15ms, are bunched in temporal analytic phase with a peak around zero. We obtain the same results for the VLD3 and for the other plasmas in the sequence.
ELMs are thus more likely to occur when the full flux loop signals are at a specific temporal analytic phase with respect to that of the preceding ELM. Prompt ELMs occur within the coherent (in amplitude and temporal phase) response to the previous ELM which can clearly be seen in the full flux loop signals. 16 For all other, non-prompt ELMs, the full flux loop signals do not remain coherent in both amplitude and temporal phase throughout the inter-ELM time interval. The question is then whether there is detectable temporal phase coherence at all times (implying that the system retains a memory of the preceding ELM) or whether temporal phase coherence is lost, and then re-emerges as part of the build-up to the next ELM. Figure 5 shows polar plots of the histogram of the temporal analytic phase differences D/ k ðdtÞ for all of the ELM pairs in plasma 83771. The temporal analytic phase difference is determined at time t ¼ t ELM2 À dt, that is, at time dt before the second ELM. As the preceding ELM generates a large amplitude response in the full flux loop signals on a timescale $10ms, we will exclude ELM pairs with waiting times Dt < 15ms þ dt; hence the number N of samples in the histogram decreases with increasing dt. As in Figure 4, we use Hilbert transform window (a) which extends beyond the time t ELM2 of the second ELM, so that we can obtain the instantaneous temporal analytic phase at times up to the ELM occurrence time. The bottom panels in Figure 5 are at the time of the ELM, dt ¼ 0, so that the bottom left hand plot is a polar histogram of the same data as in the right hand panel of Figure 4. The time before the ELM dt increases moving up the plot. We then see that the temporal analytic phase difference qualitatively becomes progressively more ordered from dt $ 5ms and that there is a clear temporal phase bunching after dt $ 2ms. This suggests that there is a signature of the build-up to an ELM in the full flux loop signals and we will quantify the degree of temporal phase bunching in Sec. V. Figures 4 and 5 required a Hilbert transform window (a) which extended beyond the time t ELM2 of the second ELM, so that the instantaneous temporal analytic phase at times up to the ELM occurrence time could be calculated. As a check on the robustness of the ELM build-up signature, we repeat the analysis with Hilbert transform window (b) which stops at the time t ELM2 of the second ELM, so that only information before the ELM occurrence time is used. The resulting polar histograms are shown in Figure 6, where apart from the different Hilbert transform window, the data and analysis are the same as that used to produce Figure 5. Now, we can only consider times before the second ELM dt > dt E ¼ 1ms. Nevertheless, we still see in these histograms a clear temporal phase bunching on the same timescales as in Figure 5, where information from times beyond the ELM time t ELM2 was used.
ELM occurrence times are marked on each VLD2 trace with yellow-filled red circles (first ELM) and yellow-filled green circles (second ELM). Right hand panel: histogram of VLD2 D/ at the time of all the second ELMs (blue), overplotted (pink) for the prompt ELMs with waiting times Dt < 15ms. Top panel: histogram of ELM occurrence times Dt ¼ t À t 0 for the first ELM (red) and the second ELMs (green), overplotted (pink) for the prompt ELMs. The frequency N of first ELM times has been rescaled by 1/10.
The above results test for temporal coherence in the build-up to an ELM, that is, over what time interval do we always see the same temporal analytic phase in the VLD2 or VLD3 just before an ELM. The VLD2 and VLD3 full flux loops are both in the divertor region of JET and in Figure 2 we can see that they are very similar in their time variation, however they are not identical. The time evolving difference in temporal analytic phase between these signals provides a measure of spatial coherence, that is, coherent large scale plasma motion in the region of these full flux loops will tend to make their phases align. We test this idea in Figure 7 where we plot polar histograms of the instantaneous temporal analytic phase difference between the VLD2 and VLD3 signals at times dt before each of the ELMs in plasma 83771, in the same format as Figure 5. From the top panel, we see that their temporal analytic phase difference at all times shows some alignment, it is within $660 of its mean at dt ¼ 10ms before the ELM. However, again for times dt < 5ms, that is, just before the ELM, we see that the phase difference in these two signals tends to zero, that is, they become temporally phase synchronized.

V. CIRCULAR STATISTICS AND THE RAYLEIGH TEST
We use the Rayleigh test and associated circular statistics (see, e.g., Refs. 23 and 24, and references therein) to quantify the extent to which the temporal analytic phase FIG. 5. Polar plots of histograms of VLD2 (left panels) and VLD3 (right panels) temporal analytic phase difference D/ just before an ELM in the flat-top of JET plasma 83771. The phase difference is calculated from the time of the first ELM to a time dt before the second ELM, so that D/ ¼ /ðt ELM2 À dtÞ À /ðt ELM1 Þ. From bottom to top dt ¼ ½0, 1, 2, 5, 10ms. The histograms include all ELM pairs with waiting times Dt > 15ms À dt, the number N in the histogram decreases with increasing dt. Hilbert transform time window (a) is used to determine the VLD temporal analytic phases and it extends dt E ¼ 1ms beyond t ELM2 , the time of the second ELM. The interval used to determine the VLD signal means just before the ELM moves back with dt.
FIG. 6. Polar plots of histograms of VLD2 (left panels) and VLD3 (right panels) temporal analytic phase difference D/ just before an ELM in the flat-top of JET plasma 83771. The same data and analysis are used as in Figure 5 except that now Hilbert transform time window (b) is used to determine the VLD temporal analytic phases and it stops at t ELM2 , the time of the second ELM. Phase differences can then only be determined for times dt > dt E ¼ 1ms before t ELM2 . From bottom to top dt ¼ ½1, 2, 5, 10ms. differences are aligned, and the statistical significance of any such alignment. Using the procedure described above, we determine the temporal analytic phase differences D/ k for the k ¼ 1; :::; N ELM pairs in given plasma. If each temporal phase is represented by a unit vector r k ¼ ðx k ; y k Þ ¼ ðcos D/ k ; sin D/ k Þ, then a measure of their alignment is given by the magnitude of the vector sum, normalized to N. This is most easily realized if we use unit magnitude complex variables to represent the r k ¼ e iD/ k . Then if the Rayleigh number is the magnitude of the sum and the mean temporal analytic phase angle is /. Clearly, if R ¼ 1 the temporal phases are completely aligned, however R ¼ 0 does not distinguish random alignment from ordered anti-alignment. We will consider two other statistics here. The first is an estimate of how closely aligned the temporal phases are with the mean phase angle. We can calculate centred trigonometric moments relative to the mean phase angle / We will consider q ¼ 2, then the temporal analytic phase angle of m 1 2 2 , that is, d/ 2 =2 is a measure of the angular variance around the mean /; this can take values ½0 À 6p. We will plot this quantity as a standardized, positive definite, angular variance r / ¼ jd/ 2 j=2p, so that r / is in the range of ½0; 1.
The second is an estimate of the p-value under the null hypothesis that the vectors are uniformly distributed around the circle which is given by so that a small value of p indicates significant departure from uniformity, i.e., the null hypothesis can be rejected with 95% confidence for p < 0.05. We now calculate the Rayleigh R, the standardized angular variance, and p values as a function of the time before the ELM dt, corresponding to the polar histograms above. Figure 8 corresponds to the analysis of Figure 5, where we have used Hilbert time window (a) to obtain the FIG. 7. Polar plots of histograms of the difference instantaneous temporal analytic phase between the VLD2 and VLD3 signals just before an ELM in the flat-top of JET plasma 83771. The temporal phase difference is calculated directly between the two signals, so that D/ ¼ /ðVLD2ðtÞÞ À /ðVLD3ðtÞÞ. The format is the same as in Figure 5. Hilbert transform window (a) is used to determine the phases.
FIG. 8. Rayleigh statistics for VLD2 difference in temporal analytic phase D/ just before an ELM in the flat-top of JET plasma 83771. The temporal analytic phase difference is calculated from time t 0 to a time dt before the second ELM, so that D/ ¼ /ðt ELM2 À dtÞ À /ðt 0 Þ. The figure plots variation with dt (x axis) of Rayleigh R for t 0 ¼ t ELM1 (green) and for t 0 ¼ t VLDmin (blue). The corresponding standardized angular variance is plotted as light blue and green shading. The corresponding p-values are indicated by the dark and light grey shading, respectively. The p ¼ 0.05 level is indicated by the horizontal dashed black line. The sample includes all ELM pairs with inter-ELM time intervals Dt > 15ms þ dt, the number N in the sample decreases with increasing dt; the fraction NðdtÞ=Nðdt ¼ 0Þ of ELM pairs in the sample is plotted (black). Hilbert transform window (a) is used. VLD2 difference in temporal analytic phase D/ k at times dt just before each ELM. The temporal analytic phase difference is again calculated from the time of the first ELM to a time dt before the second ELM, so that D/ k ¼ /ðt ELM2 À dtÞ À/ðt 0 Þ. In Figure 8, the green line is the Rayleigh R for the analysis of Figure 5 where we calculate the temporal analytic phase differences from the zero time at the first ELM t 0 ¼ t ELM1 . The times of the extrema of the characteristic initial large amplitude oscillatory response to an ELM, which is seen in both the full flux loop signals, have been found 1 to provide a better determined zero time t 0 . The blue line in Figure 8 is the R obtained for t 0 ¼ t VLDmin , the first minimum of the VLD2 signal following the preceding ELM. We can then see that R > 0.3 for dt < 5ms before the ELM occurs, and systematically increases as we approach the ELM occurrence time. Within this time interval, the standardized angular variance r / is small, it gradually increases with dt as the temporal analytic phases become more disordered. The p-statistic remains small for times dt < 15ms indicating that the distribution of temporal analytic phases remains far from circular. However, this is not a smooth trend, there are short intervals for example, around dt $ 6ms where p $ 0:05. We have found that for dt > 5ms, the details of where shortlived fluctuations in R, r / , and p-value occur are not robust, they vary with the dataset and with the detailed parameters of how the Hilbert transform is computed. However, the overall trends are robust, in particular, alignment of the temporal analytic phases around a single value for dt < 5ms, that is, large R and small r / and p-value.
In Figure 6, we only used signals up to, and not beyond, the time of the ELM in order to test for the ELM build-up signature in the full flux loop temporal analytic phases. The corresponding circular statistics are plotted in Figure 9 where the temporal analytic phase differences are obtained for t 0 ¼ t VLDmin . The blue line in Figure 9 replots that in Figure 8, it is calculated using Hilbert transform time window (a) which extends to times beyond the ELM occurrence time. The red line in Figure 9 is obtained using the same analysis and data, but with temporal analytic phases calculated using Hilbert transform time window (b) which stops at the ELM occurrence time. We can see that the build-up to an ELM in the full flux loops can still be resolved only using information from before the ELM occurrence time.
Finally, in Figure 10 we plot the Rayleigh statistics for the difference in instantaneous temporal analytic phase between the VLD2 and VLD3 signals that was shown in Figure 7. These signals are very similar in their time variation as can be seen in Figure 2, however they are not identical. From top panel of Figure 10, we see that their temporal analytic phase difference at all times shows some alignment, it is within $660 of its mean so that R $ 0:5 in Figure 10. However, again for times dt < 5ms, that is, just before the ELM, we see that the difference in temporal analytic phase between these two signals tends to zero, that is, they become temporally phase synchronized.
We have quantified the values that these circular statistics can take for these time-series due to chance coincidence. Chance coincidence can occur between time-series that have non-trivial time-structure, for example, roughly periodic ELM occurrence times may preferentially occur at specific temporal phases of a roughly sinusoidal signal. We have constructed a set of surrogate time-series and repeated the above analysis to explore this possibility. This is described in detail in the Appendix, and establishes that the temporal phase alignments seen for times dt < 5ms can be distinguished as statistically distinct from chance occurrence and thus are evidence for correlation.
Our main result is that there is a signature of the buildup a non-prompt ELM in the temporal analytic phases of the FIG. 9. Rayleigh test for VLD2 difference in temporal analytic phase D/ just before an ELM in the flat-top of JET plasma 83771. The format is as in the previous figure. The temporal analytic phase difference is calculated from time t 0 to a time dt before the second ELM, so that D/ ¼ /ðt ELM2 À dtÞ À/ðt 0 Þ where t 0 ¼ t VLDmin . The blue line is the same as the previous figure, Hilbert transform window (a) is used to determine the temporal analytic phases. The red line is the same analysis using Hilbert transform window (b). These time windows end at t ELM2 þ 1ms (blue) and t ELM2 (red), respectively. Normalized angular variance shown in shaded light blue and red, respectively. The corresponding p-values are indicated by the light and dark grey shading, respectively.
FIG. 10. Rayleigh test for the difference in instantaneous temporal analytic phase between the VLD2 and VLD3 signals just before an ELM in the flattop of JET plasma 83771. The format is as in the previous figure. The temporal analytic phase difference is calculated directly between the two signals, so that D/ ¼ /ðVLD2ðtÞÞ À /ðVLD3ðtÞÞ. The green line is the R value, Hilbert transform window (a) is used to determine the temporal analytic phases. Normalized angular variance shown in shaded light green and the p-values by the light grey shading, here it is too small to be visible on the plot. full flux loop signals. ELMs tend to occur preferentially at a specific temporal analytic phase in the VLD2 and VLD3 signals. The R, r / , and p-values are all consistent with alignment of the temporal analytic phases from about 2-5 ms just before the ELM occurs. Furthermore, the VLD2 and VLD3 signals become temporally phase synchronized with each other during this build-up time. Global, spatio-temporally synchronized plasma dynamics is thus part of the build-up to an intrinsic ELM. We cannot detect statistically significant temporal phase coherence at all times, whereas we do detect temporal phase coherence re-emerging as part of the buildup to the next ELM.

VI. DISCUSSION
Our results rely upon a new approach to the analysis of an existing JET diagnostic, the full flux loops VLD2 and VLD3 signals, alongside ELM timings from the Be II signal. We have performed direct time-domain analysis of high time resolution toroidally integrating full flux loop signals arising from the dynamics of spatially integrated current density, with high time resolution determination of the ELM timings. In addition to the main results of this paper, in this section we will develop our recent conjecture 17 in the light of these results. Our aim is to frame a testable hypothesis for future work.
It is well established experimentally that ELMs can be triggered by applied magnetic "kicks" delivered by the vertical stabilization control coils which drive vertical plasma movement, including global changes in the divertor region. In such triggering experiments in JET, ELMs preferentially occur when the plasma is in a specific temporal phase in its vertical motion (downwards), and delays of $2-3 ms typically are observed between the start of the kick and the ELM. 25,26 Similar behaviour, i.e., ELM occurrence when the plasma is in a specific temporal phase in its vertical motion, is seen in other devices, e.g., Refs. 27 and 28, and references therein. Furthermore, the velocity perturbation associated with intrinsic ELMs is found to set a minimum threshold value that must be exceeded in order to trigger ELMs with the vertical coils. 27 The build-up to an ELM that has been magnetically kicked thus involves global plasma motion at a specific temporal phase. This global plasma displacement can then modify conditions at the plasma edge, such that peeling-ballooning and perhaps other instabilities become active, leading to the ELM burst. The details are complex and may be device dependent; 30 but the essential point here is that the kicked ELM burst follows a global perturbation in the plasma dynamics and occurs at a specific temporal phase thereof.
We have presented evidence for the emergence of coherent global dynamics in the integrated current density in the $2-5 ms build-up to an intrinsic ELM. This signature of the build up to an ELM is global, albeit perhaps poloidally localised, in that (i) it can be found in timeseries that are toroidally integrating and (ii) it is seen at both full flux loops which are at different radial positions in the divertor region. In a plasma that remains close to a global magnetic equilibrium, this can reflect bulk displacement or motion of the plasma. We see this build-up in the full flux loop signals which track the dynamics of the integrated current density in the divertor region. The VLD2 and VLD3 signals become temporally phase synchronized during this build-up, suggesting a spatially coherent large-scale plasma perturbation. The intrinsic ELMs are found to preferentially occur at a specific temporal analytic phase in the full flux loop signals, that is, at a specific temporal phase in this global perturbation in the plasma. If this global perturbation is sufficient to modify conditions at the plasma edge to favour instability, then an intrinsic ELM can occur.
Our results suggest one possible scenario for intrinsic ELMing where the plasma and its interacting environment together self-generate a global plasma perturbation, such that the plasma is magnetically "self-kicked," which then leads to an ELM. Self-generation of global motion could occur via nonlinear feedback between the multiscale dynamics of the plasma and its interacting environment, including the control system, as we first suggested in Ref. 17. The steady state of the JET flat top plasmas is actively maintained by perturbations from the control system reacting to plasma motion. Integrated over the largest spatial scales, the reaction of the plasma to these perturbations is seen in the full flux loop signals. These signals reflect the control system and plasma behaving as a single nonlinearly coupled system, rather than as driver and response. If there were coupling between the global plasma environment, including the control system, and each of several growing modes in the plasma, these modes could become synchronized, 20-22 through their individual interactions with the global plasma/control system environment, without the need of coupling between the modes themselves. Large scale plasma motion would then develop on timescales characteristic of the dynamics of the global plasma environment. We have found an ELM buildup timescale of $2-5 ms, which is similar both to the $2ms time constant of the known unstable mode in the vertical control system on JET, 29 and the $2-3 ms response time to generate global plasma motion from active kicks in the vertical stabilization control coils. 25,26 The VLD2 and VLD3 full flux loops also capture the initial integrated plasma and control system response to an ELM. 1 If this integrated plasma and control system response again corresponds to global plasma motion, it may be expected to act as a "kick" to directly trigger an ELM, if this global perturbation is sufficient to modify conditions at the plasma edge for instability. We found 16 that prompt ELMs sometimes occur at a specific temporal phase within this initial response to the previous ELM. This suggests an additional testable hypothesis: that compound ELMs are a pattern of successive prompt ELMs and again arise from global plasma motion emerging as above. This is consistent with the observation 12 of a narrow spread in the time intervals between successive component ELMs in a compound ELM sequence. We would then expect to see a well-defined temporal phase relationship between high time resolution full flux loop signals and the burst occurrence times within compound ELMs.
Although the above is a conjecture, it frames hypotheses that are testable by direct time-domain analysis of the relevant signals if they can be obtained at sufficiently high time resolution, pointing to future work that may further the understanding of the ELMing process.

VII. CONCLUSIONS
We have performed direct time domain analysis of ELMing in JET plasmas where a steady H-mode is sustained over several seconds, during which there is no deliberate intent to control the ELMing process by external means. We identified the ELM occurrence times from the Be II signal and have determined their relationship with the temporal analytic phase of the VLD2 and VLD3 toroidally integrating full flux loop current signals, which are a high time resolution global measurement proportional to the voltage induced by changes in poloidal magnetic flux in the divertor region.
We have established that there is a signature of the build-up to an ELM in the temporal analytic phases of the full flux loop signals. Just before an ELM, the full flux loop temporal analytic phases progressively align such that at the ELM, they have the same value as at the previous ELM. This alignment is seen to develop over the $2-5 ms before the ELM. It is sufficiently strong that it can be distinguished from temporal phase relationships that could occur by coincidence in these quasi-oscillatory signals. We are able to recover this build-up signature using only data from before the ELM occurrence time. It thus possesses predictive power. While the full flux loops track each other at all times, that is, they have a temporal phase relationship with each other that is distinct from random, they become strongly temporally phase synchronized within this build-up time before an ELM, consistent with globally spatially coherent plasma dynamics in the divertor region.
These results may assist ELM prediction and mitigation, in that real time knowledge of the full flux loop signal temporal analytic phases indicates future times when ELM occurrence is statistically more likely. The full flux loop signals capture aspects of the global dynamics of the plasma, including large scale plasma motion, plasma dynamics in the divertor region, and mutual interaction with the control system. Our result may thus provide new insight into the ELMing process. We suggest a possible scenario that unifies our understanding of intrinsic ELMing, and magnetic pacing of ELMs that uses the vertical stabilization coils to drive bulk plasma motion. The full flux loop signals can be seen in the time-series to have intervals where there is a clear sinusoidal component, with a characteristic period of $10ms and the ELM waiting times have time structure; they are not random. The analysis is performed on a restricted sized sample. We now test a series of null hypotheses that capture scenarios where the temporal phase alignment that we report above could occur by coincidence. We will use the same circular statistics as above to distinguish the likelihood of coincidental occurrence in a quantitative manner. We use Hilbert transform window (a) and the same dataset as in the main paper, plasma 83771, to construct the surrogates. We have repeated this analysis for all the other JET plasmas in this sequence, and we obtain the same results.

ELM time and instantaneous temporal analytic phase of one full flux loop signal
We need to quantify the temporal phase alignment that could occur due to coincidence in comparing the ELM arrival times with a single signal, one of the full flux loops. If, for example, the full flux loop signals were simply monochromatic sinusoids and the ELMs occurred in a sufficiently periodic fashion, one would see ELMs preferentially occurring at particular temporal analytic phases in the full flux loop signals whether or not the sequence of ELM occurrence times and the full flux loop signals were related to each other.
We therefore test the statistical significance of the above results against some alternative hypotheses. We can represent these alternative hypotheses by constructing surrogate time-series that retain some, but not all, of the properties of the original data. We aim to test that the above results are significant compared to a random process. We also aim to quantify trivial correlation, that is, coincidences between ELM arrival time and full flux loop temporal analytic phase. Coincidences could arise in a finite dataset where both the sequence of ELM arrival times, and the full flux loop signals contain time structure that includes periodicity. Here, the ELM waiting times have a mean period and a "comb like" multi-periodic structure, and the full flux loop signals exhibit intervals of oscillatory behaviour. We will calculate the same circular statistics in exactly the same manner as above for the following surrogate datasets, the results are shown in the three panels of Figure 11. (C) Full flux loops are single constant frequency sinusoids ( Figure 11, bottom panel): We replace the full flux loop signal with a sinusoid with period T ¼ 10ms, approximately the characteristic period of the oscillatory response seen following an ELM. We trial two surrogates: the first is a single sinusoid running through the entire time-series, and second, we reset the temporal analytic phase of the sinusoid to zero at the time if the first ELM in each pair, that is, at the start of each ELM waiting time.
From these surrogates, we can conclude the following. First, surrogate A establishes the value of R $ 0:1 that occurs from the temporal phases in a random signal. Here, p > 0.05 so that the distribution of temporal phases is indistinguishable from circular, they randomly occur at all angles. Surrogate B preserves both the full structure of the VLD2 signal and the probability distribution of ELM waiting times. Now, the ELM waiting time distribution has time structure, some waiting times occur more frequently than others. In a finite sized sample, randomly permuting them cannot generate coincidences with all possible temporal analytic phases of the VLD2 signal and this will lead to some alignment. As we move through the VLD2 time series by varying dt the degree of alignment will fluctuate. This can indeed be seen to give R $ 0:3 which is larger than the random signal surrogate A. In the dt $ 5ms before the surrogate ELM time, and both the angular variance r / and p-value are small so that there is some alignment. This sets an upper bound for R and a lower bound for r / which can occur by such coincidences. Finally, surrogate C produces p > 0.05 everywhere except dt < 3ms. For dt > 3ms, the distribution of temporal phases is indistinguishable from circular, they randomly occur at all angles. At smaller dt, there is again some alignment, which reaches a similar alignment, that is R and angular variance r / , as in surrogate B, and for the same reason, the ELM waiting times have preferred values and these preferentially coincide with some temporal phases of the single sinusoid surrogate.
Comparing these surrogates with our result of Figure 8, we conclude that the alignment seen for dt < 5ms is statistically significant and cannot be accounted for by chance coincidence between the sequence of ELM occurrence times and the temporal analytic phase of the full flux loop signals. The alignment in 5 < dt < 15ms is stronger than that of a random process (R $ 0:1: surrogate A) but is comparable with that arising from temporal phase coincidence (R $ 0:3, surrogates B and C) and thus cannot be distinguished from it.

Temporal analytic phase difference between VLD2 and VLD3 signals
The full flux loop signals both contain time structure that includes periodicity, on roughly the same period T ¼ 10ms. We now test against the coincidence that could occur in the temporal analytic phase difference between sinusoidal signals sampled at a sequence of times (the ELM arrival times) that have time structure.
FIG. 11. Rayleigh statistics for VLD2 temporal analytic phase difference D/ just before an ELM for three different surrogate time-series for the flattop of JET plasma 83771. Format is as in previous figures. Top: the individual data points of the VLD2 signal have been randomly permuted within each ELM waiting time. Middle: the sequence of ELM waiting times is randomly permuted. Bottom: the VLD2 signal is replaced by a single sinusoid of period T ¼ 10ms throughout the entire sequence of ELMs (blue); the sin temporal analytic phase is reset to zero at the time of the previous ELM (red). Red and blue shading are the normalized angular variance. Light and dark grey are the p-value of single and reset sin waves, respectively.
(A) ELM arrives at a random time (Figure 12, top panel): For each ELM pair, we randomly select a time within the time interval to the next ELM, that is, the second ELM arrives at a random time. (B) No pattern in the sequence of ELM waiting times ( Figure 12, bottom panel): We randomly permute the time sequence of ELM waiting times fDt 1 ; Dt 2 ; :::;Dt j ::Dt N g as in surrogate B above. (C) One of the full flux loops is a single constant frequency sinusoid ( Figure 13): We replace the VLD2 full flux loop signal with a sinusoid with period T ¼ 10ms, the characteristic period of the oscillatory response seen following an ELM. We trial two surrogates, the first is a single sinusoid running through the entire time-series, and second, we reset the temporal analytic phase of the sinusoid to zero at the time if the first ELM in each pair, that is, at the start of each ELM waiting time. The results are similar, one case is shown.
These surrogates establish that the full flux loops are similar in temporal analytic phase at all times, surrogates A and B have an R $ 0:5. Comparing Figure 10, we see that this is the value at dt > 5ms. Here, we use the unchanged VLD3 and ELM arrival times. The VLD2 signal is replaced by a single sinusoid of period T ¼ 10ms and the sin temporal analytic phase is reset to zero at the time of the previous ELM.