Journal topic
Ann. Geophys., 36, 497–507, 2018
https://doi.org/10.5194/angeo-36-497-2018

Special issue: Space weather connections to near-Earth space and the...

Ann. Geophys., 36, 497–507, 2018
https://doi.org/10.5194/angeo-36-497-2018

Regular paper 23 Mar 2018

Regular paper | 23 Mar 2018

# Non-Gaussianity and cross-scale coupling in interplanetary magnetic field turbulence during a rope–rope magnetic reconnection event

Non-Gaussianity and cross-scale coupling in interplanetary magnetic field turbulence during a rope–rope magnetic reconnection event
Rodrigo A. Miranda1,2, Adriane B. Schelin2, Abraham C.-L. Chian3,4,5, and José L. Ferreira2 Rodrigo A. Miranda et al.
• 1UnB-Gama Campus, University of Brasília (UnB), Brasília DF 70910-900, Brazil
• 2Plasma Physics Laboratory, Institute of Physics, University of Brasília (UnB), Brasília DF 70910-900, Brazil
• 3School of Mathematical Sciences, University of Adelaide, Adelaide SA 5005, Australia
• 4Institute of Aeronautical Technology (ITA), São José dos Campos, SP 12228-900, Brazil
• 5National Institute for Space Research (INPE), P.O. Box 515, São José dos Campos, SP 12227-010, Brazil

Correspondence: Rodrigo A. Miranda (rmiracer@unb.br)

Abstract

In a recent paper it was shown that magnetic reconnection at the interface region between two magnetic flux ropes is responsible for the genesis of interplanetary intermittent turbulence. The normalized third-order moment (skewness) and the normalized fourth-order moment (kurtosis) display a quadratic relation with a parabolic shape that is commonly observed in observational data from turbulence in fluids and plasmas, and is linked to non-Gaussian fluctuations due to coherent structures. In this paper we perform a detailed study of the relation between the skewness and the kurtosis of the modulus of the magnetic field $|\mathbit{B}|$ during a triple interplanetary magnetic flux rope event. In addition, we investigate the skewness–kurtosis relation of two-point differences of $|\mathbit{B}|$ for the same event. The parabolic relation displays scale dependence and is found to be enhanced during magnetic reconnection, rendering support for the generation of non-Gaussian coherent structures via rope–rope magnetic reconnection. Our results also indicate that a direct coupling between the scales of magnetic flux ropes and the scales within the inertial subrange occurs in the solar wind.

Keywords. Space plasma physics (turbulence)

1 Introduction

The solar wind can be regarded as a network of entangled magnetic flux tubes and Alfvénic fluctuations propagating within each flux tube . Flux tubes can emerge locally in the solar wind as a consequence of the magnetohydrodynamic turbulent cascade . An alternative view describes coherent structures as “fossile” structures that emanate from the solar surface and are advected by the solar wind .

The probability distribution functions (PDFs) of turbulent space plasmas display sharp peaks and fat tails on small scales within the inertial subrange , as well as departures from self-similarity and monofractality . These features are due to the presence of rare, large-amplitude coherent structures which dominate the statistics of fluctuations on small scales and can be quantified by the computation of statistical moments.

A robust parabolic dependence between the normalized third-order moment (skewness) and the normalized fourth-order moment (kurtosis) has been found in local concentrations of contaminants in atmospheric turbulence as found by . also found a similar skewness–kurtosis parabolic relation using global data of sea-surface temperature fluctuations. reported a similar skewness–kurtosis dependence in electron density fluctuations in plasma confinement experiments. obtained a skewness–kurtosis parabolic relation for datasets of human reaction times for visual stimuli. Since then, the presence of a skewness–kurtosis relation in different physical scenarios has attracted much attention and has been associated with the presence of non-Gaussian fluctuations due to coherent structures . The skewness–kurtosis parabolic relation was also found in time series of two-point differences of the modulus of the magnetic field by . They demonstrated that the parabolic relation is due to nonlocal interaction between large-scale structures and small-scale intermittency.

In this paper we investigate the skewness–kurtosis relation during a triple interplanetary magnetic flux rope (IMFR) event detected by Cluster-1 in the solar wind. This event was recently characterized by . They demonstrated the occurrence of magnetic reconnection at the interface region of two IMFRs and that this reconnection can be the origin of interplanetary intermittent turbulence. Our results show that the skewness–kurtosis parabolic relation is enhanced during the reconnection between flux ropes, and that is a natural consequence of the interaction between flux ropes.

This paper is organized as follows. Section 2 presents the statistical tools employed for the data analysis, including the equations to compute the skewness and the kurtosis. Section 3 describes the triple-IMFR event. The skewness–kurtosis relation is analyzed in detail in Sect. 4. The interpretations of these results are presented in Sect. 5. Finally, we conclude in Sect. 6.

2 Data analysis tools

Let θi, $i=\mathrm{1},\mathrm{\dots },N$ be the time series of a quantity of interest (e.g., the modulus of the magnetic field $|\mathbit{B}|$). The skewness of θi can be computed as follows:

$\begin{array}{}\text{(1)}& S=\frac{\mathrm{1}}{N}\sum _{i=\mathrm{1}}^{N}{\left(\frac{{\mathit{\theta }}_{i}-〈{\mathit{\theta }}_{i}〉}{\mathit{\sigma }}\right)}^{\mathrm{3}},\end{array}$

where θi represents the average of θi, N represents the number of data points, and σ is the SD of θi. The flatness of θi is given by

$\begin{array}{}\text{(2)}& F=\frac{\mathrm{1}}{N}\sum _{i=\mathrm{1}}^{N}{\left(\frac{{\mathit{\theta }}_{i}-〈{\mathit{\theta }}_{i}〉}{\mathit{\sigma }}\right)}^{\mathrm{4}},\end{array}$

from which the kurtosis can be obtained by

$\begin{array}{}\text{(3)}& K=F-\mathrm{3}.\end{array}$

For a Gaussian function $S=K=\mathrm{0}$. The skewness quantifies the degree of asymmetry of the PDF of θi, whereas the kurtosis quantifies the departure of the flatness of the PDF of θi from the flatness of a Gaussian distribution which is equal to 3. The definition of kurtosis in Eq. (3) is sometimes called “excess kurtosis” .

A common way to characterize asymmetry and non-Gaussianity of θi as a function of scale τ is through the time series of two-point differences:

$\mathit{\delta }{\mathit{\theta }}_{i}\left(\mathit{\tau }\right)={\mathit{\theta }}_{i+\mathit{\tau }}-{\mathit{\theta }}_{i}.$

The skewness of θi on scale τ is then

$\begin{array}{}\text{(4)}& S\left(\mathit{\tau }\right)=\frac{\mathrm{1}}{N}\sum _{i=\mathrm{1}}^{N}{\left(\frac{\mathit{\delta }{\mathit{\theta }}_{i}-〈\mathit{\delta }{\mathit{\theta }}_{i}〉}{{\mathit{\sigma }}_{\mathit{\tau }}}\right)}^{\mathrm{3}},\end{array}$

and the flatness is

$\begin{array}{}\text{(5)}& F\left(\mathit{\tau }\right)=\frac{\mathrm{1}}{N}\sum _{i=\mathrm{1}}^{N}{\left(\frac{\mathit{\delta }{\mathit{\theta }}_{i}-〈\mathit{\delta }{\mathit{\theta }}_{i}〉}{{\mathit{\sigma }}_{\mathit{\tau }}}\right)}^{\mathrm{4}},\end{array}$

where στ is the SD of δθi(τ). From Eq. (5) the kurtosis as a function of scale is obtained by

$\begin{array}{}\text{(6)}& K\left(\mathit{\tau }\right)=F\left(\mathit{\tau }\right)-\mathrm{3}.\end{array}$

A functional relation between the skewness and the kurtosis of θi as defined by Eqs. (1)–(3) has been observed in a variety of scenarios . This relation is given by

$\begin{array}{}\text{(7)}& K=\mathit{\alpha }{S}^{\mathrm{2}}+\mathit{\beta },\end{array}$

where α and β are the coefficients that characterize a parabolic curve.

We compute the α and β coefficients by applying a least-square fit between (S,K) values obtained from the observational data and Eq. (7) following the Levenberg–Marquardt algorithm , which is a popular method to fit a dataset into nonlinear equations. In order to quantify how well the computed (S,K) values are fitted into Eq. (7) we employ the correlation index r which measures the correlation between two datasets Xi and Yi, $i=\mathrm{1},\mathrm{\dots },N$:

$\begin{array}{}\text{(8)}& r=\frac{\mathrm{1}}{{\mathit{\sigma }}_{X}{\mathit{\sigma }}_{Y}}\sum _{i=\mathrm{1}}^{N}\frac{\left({X}_{i}-〈{X}_{i}〉\right)\left({Y}_{i}-〈{Y}_{i}〉\right)}{N},\end{array}$

where σX and σY represent the SD of Xi and Yi, respectively. The correlation index $r\in \left[-\mathrm{1},\mathrm{1}\right]$. If r=1 there is complete correlation between Xi and Yi, whereas $r=-\mathrm{1}$ indicates anticorrelation. The value r=0 represents absence of correlation.

In summary, the analysis is described by the following steps:

• Compute S and K from the modulus of magnetic field $|\mathbit{B}|$ using Eqs. (1)–(3).

• Apply the Levenberg–Marquardt algorithm to find α and β in Eq. (7) that best fit the (S,K) values.

• Use α and β obtained from the previous step in Eq. (7) to obtain empirical values of K as a function of S.

• Compute the correlation index r between the values of K from the previous step and the values of K from the observational data. The r index will measure how close the K values computed by Eq. (3) are to the K values obtained empirically from Eq. (7).

We repeat these steps for S(τ) and K(τ) of two-point differences using Eqs. (4) and (6) in the first step. There are several computational programs for data analysis that implement the Levenberg–Marquardt algorithm. Here we use the implementation available in the GNU Octave program .

We note that several papers regarding the relation between skewness and kurtosis have employed the definition of what we refer to as flatness (Eq. 2). Throughout this paper we will focus on the kurtosis defined by Eq. (3).

3 A triple-IMFR event

Figure 1a shows the time series of the modulus of magnetic field $|\mathbit{B}|$ obtained by the FGM instrument onboard Cluster-1 from 00:00 to 12:00 UT on 2 February 2002. During this interval Cluster-1 was in the solar wind upstream of the Earth's bow shock . The magnetic field data are collected by Cluster-1 at a resolution of 22 Hz . Figure 1 also presents an overview of other in situ plasma parameters for the selected interval, namely, the three components of B in the GSE coordinates, the angles ΦB and ΘB of the solar wind magnetic field B relative to the Sun–Earth x axis in the ecliptic plane, and out of the ecliptic, respectively, in the polar GSE coordinates; the modulus of the ion bulk flow velocity $|{\mathbf{V}}_{\mathrm{i}}|$, the ion number density ni, the ion temperature perpendicular to the magnetic field Ti and the ion plasma βi, which is the ratio between plasma kinetic pressure and magnetic pressure. The Cluster-1 plasma measurements are given by the ion spectrometry experiment CIS .

This event is characterized by the presence of three interplanetary magnetic flux ropes. Magnetic flux ropes are magnetic structures described as bundles of twisted, current-carrying magnetic field lines bent into a tube-like shape, spiralling around a common axis . During this event three IMFRs were identified by using a combination of criteria for large-scale magnetic cloud boundary layers and small-scale IMFRs . The interval of each IMFR is indicated by horizontal arrows in Fig. 1a, and their timings are shown in Table 1.

Figure 1Cluster-1 magnetic field and plasma parameters from 00:00 to 12:00 UT on 2 February 2002. From top to bottom: modulus of magnetic field $|\mathbit{B}|$ (nT), three components of B (nT) in the GSE coordinates, azimuth angle ΦB (), latitude angle Θ (), modulus of ion bulk velocity $|{\mathbf{V}}_{\mathrm{i}}|$ (km s−1), ion number density ni (cm−3), ion temperature Ti (eV) and ion plasma beta βi. Horizontal arrows indicate the interval of IMFR-1 (black), IMFR-2 (red) and IMFR-3 (blue). The front and rear boundary layers of each IMFR are indicated by the vertical dotted lines.

Table 1Beginning and end of the intervals depicted in Fig. 1, corresponding to the boundary layers of three interplanetary magnetic flux ropes (IMFRs) on 2 February 2002.

4 Skewness–kurtosis relation

## 4.1 Time series of $|\mathbit{B}|$

Figure 2a shows the time series of $|\mathbit{B}|$ detected by Cluster-1 on 2 February 2002 (Julian day 32) from 00:32 to 03:18 UT. Five regions were defined during this interval and are indicated using arrows. These regions represent the interior region of IMFR-1 (R1), the interface of IMFR-1 and IMFR-2 (I12), the interior of IMFR-2 (R2), the interface of IMFR-2 and IMFR-3 (I23), and the interior of IMFR-3 (R3). Their timings are indicated in Table 2. Each region has a duration of 30 min, which gives 40 358 data points. During this event current sheets were detected at the front boundary layer of IMFR-1 and at the interface region between IMFR-2 and IMFR-3. This interface region was identified as a source of intermittent turbulence by . A current sheet was detected at the leading edge of IMFR-1 using data from ACE and Cluster-1, and a current sheet was detected at the interface region between IMFR-2 and IMFR-3 using data from Cluster-1, ACE and Wind .

Table 2Timing of the five selected regions during the triple-IMFR event on 2 February 2002.

Figure 2(a) Time series of $|\mathbit{B}|$ from 00:32 to 03:18 UT on 2 February 2002. Five regions of 30 min each are highlighted using different colors: interior of IMFR-1 (R1, black), interface region of IMFR-1 and IMFR-2 (I12, green), the interior of IMFR-2 (R2, red), the interface of IMFR-2 and IMFR-3 (I23, violet), and the interior of IMFR-3 (R3, blue). The interval of each IMFR is indicated by horizontal arrows as in Fig. 1. (b, c) Time series of the skewness S and the kurtosis K computed using a sliding overlapping window of size 10 000 data points and a window shift of 400 data points. The SD computed in each window is represented by a gray area.

The SK parabolic relation described by Eq. (7) can be verified by computing S and K from a number of datasets corresponding to different realizations of an experiment. In the case of a time series, the parabolic relation can be tested by computing S and K using datasets extracted from the time series with sliding windows. The size of the sliding window is a critical parameter for this type of analysis. Since S and K are higher statistical moments, the number of data points inside the window should be large enough to guarantee the robust estimation of S and K. However, if the time series is divided into sliding windows with a large number of data points, then the number of (S, K) values may be insufficient to verify the parabolic relation of Eq. (7). This can be solved by defining overlapping windows; nevertheless, the overlapping cannot be too large in order to obtain a set of independent (S, K) values. To determine the optimal window size, we applied a procedure to estimate the maximum order of the statistical moment in a time series . We computed the maximum statistical order in each sliding window of size 5000 data points across the time series of Fig. 2, and a window shift of 400 data points. Then, we increased the size of the window by 1000 data points (keeping the same window shift), computed the maximum order in each window and then repeated the procedure. We found that a sliding window of size 10 000 data points is large enough for a robust estimation of moments up to the sixth order in all windows and at the same time allows a sufficient number of estimations of S an K to be obtained to test the parabolic relation of Eq. (7). Figure 2b and c show the resulting time series of S and K, respectively. The SD gives an estimation of the uncertainty of the computed S and K inside each window, and is represented using a gray area. From this figure we observe that from 02:26 to 02:35 UT the uncertainty of S and K increases due to the large variation in $|\mathbit{B}|$ at the interface between IMFR-2 and IMFR-3. A similar behavior was observed in magnetic field data during an interplanetary shock event by . The uncertainty within sliding windows that contain the large variations in $|\mathbit{B}|$ increases due to nonstationarity. Following , we exclude these windows from further analysis.

Figure 3Kurtosis K as a function of skewness S computed using overlapping windows of size 10 000 data points and a window shift of 400 data points, for (a) the interior region of IMFR-1, (b) the interface of IMFR-1 and IMFR-2, (c) the interior of IMFR-2, (d) the interface of IMFR-2 and IMFR-3, and (e) the interior of IMFR-3. In each panel, the least-square fit with the parabolic function $K=\mathit{\alpha }{S}^{\mathrm{2}}+\mathit{\beta }$ is displayed as a dashed line (see Table 3).

Figure 3 shows K as a function of S for the five regions previously defined. A least-square fit with Eq. (7) is displayed as a dashed line. Table 3 shows the resulting fit for each region, as well as the correlation index r between the points in the scatter plot and the fitted parabolic function computed using Eq. (8). Since the interpretation of α and β is under debate (see the discussion in Sect. 5) we will focus on the computed value of r.

Table 3The least-square fits of Eq. (7) computed from the scatter plots of Fig. 3, and the correlation index r for the five regions defined.

The correlation index r shown in the last column of Table 3 measures how well the data points can be adjusted by the parabolic function given by Eq. (7). All regions display r>0.5. The lowest correlation is obtained for the interval corresponding to R2, in agreement with a visual inspection of Fig. 3c. For this interval, most of the points in Fig. 3c tend to accumulate around $\left(S,K\right)=\left(\mathrm{0},\mathrm{0}\right)$, which is the value obtained for a Gaussian distribution (i.e., in the absence of coherent structures). Therefore, the interior of IMFR-2 is characterized by a low degree of non-Gaussianity and intermittency in comparison with the other intervals.

The highest value of the correlation is obtained during I23 (see Table 3). Figure 3d shows that points spread near the fitted parabola and far from the (0, 0) Gaussian point. This indicates that this interval is characterized by a higher degree of non-Gaussianity. These results are in agreement with the results of , which found that the interior of IMFR-2 has lower degrees of non-Gaussianity and phase coherence, and a nearly monofractal scaling when compared with other intervals. For the interface of IMFR-2 and IMFR-3 they observed higher degrees of non-Gaussianity and phase synchronization, and a strong departure from monofractality.

Figure 4The power spectral density (PSD, left panels) and the compensated PSD (right panels) for (a) the time series of $|\mathbit{B}|$ from 00:32 UT until 08:40 UT, (b) the IMFR-1 interior region, (c) the interface between IMFR-1 and IMFR-2, (d) the IMFR-2 interior region, (e) the interface between IMFR-2 and IMFR-3, and (f) the IMFR-3 interior region. Vertical dashed lines indicate the beginning and the end of the inertial subrange.

## 4.2 Time series of $\mathit{\delta }|\mathbit{B}|$

Next, we investigate the SK parabolic relation as a function of scale within the inertial subrange. The left side of Fig. 4a shows the power spectral density (PSD) as a function of frequency f of the time series of $|\mathbit{B}|$ from the beginning of IMFR-1 at 00:32 UT until the end of IMFR-3 at 08:40 UT. The right side of Fig. 4a shows the compensated PSD which is the original PSD multiplied by f5∕3 . The inertial subrange should appear as a frequency range in which the compensated PSD is almost horizontal. The following panels in Fig. 4 show the PSD and the compensated PSD for R1, I12, R2, I23 and R3. A common frequency range in which the compensated PSD is almost horizontal for all regions is indicated by two vertical dashed lines. From Fig. 4, the inertial subrange starts at f=0.01Hz and ends at f=0.1Hz, which correspond to scales τ=100s and τ=10s, respectively.

Figure 5Probability distribution functions (PDFs) of ΔB(τ) for τ=10s (continuous line) and τ=100s (dashed line). (a) The interior region of IMFR-1, (b) the interface of IMFR-1 and IMFR-2, (c) the interior of IMFR-2, (d) the interface of IMFR-2 and IMFR-3, and (e) the interior of IMFR-3. A Gaussian distribution function is represented by the gray area.

The intermittent aspect of interplanetary magnetic field turbulence can be demonstrated by constructing the PDF of the normalized magnetic-field differences

$\mathrm{\Delta }B\left(\mathit{\tau }\right)=\frac{\mathit{\delta }B\left(\mathit{\tau }\right)-〈\mathit{\delta }B〉}{{\mathit{\sigma }}_{B}},$

where $\mathit{\delta }B\left(\mathit{\tau }\right)=|\mathbit{B}\left(t+\mathit{\tau }\right)|-|\mathbit{B}\left(t\right)|$, and the brackets denote the average value. Figure 5 shows the PDFs of ΔB constructed from the magnetic field fluctuations of the five regions, for τ=10s and τ=100s. From this figure it is clear that the PDFs are closer to a Gaussian distribution (represented by the gray area in Fig. 5) at τ=100s (large scale), and become non-Gaussian at τ=10s (small scale), exhibiting sharp peaks and fat tails. This figure demonstrates that magnetic field fluctuations become more intermittent as the scale τ becomes smaller.

Next, we analyze the SK relation of ΔB(τ) at τ=10s and τ=100s. Figure 6a shows the time series of ΔB(τ=100 s). Figure 6b and c show the time series of S and K computed using a sliding overlapping window as in Sect. 4.1. The gray area indicates the uncertainty of the S and K values. As in Fig. 2, we observe a large uncertainty from 02:26 to 02:35 UT due to the interface between IMFR-2 and IMFR-3; therefore these S and K values are excluded from further analysis.

Figure 6(a) Time series of δB(τ=100 s) from 00:00 to 04:00 UT on 2 February 2002. Five regions of 30 min each are highlighted using different colors: interior of IMFR-1 (R1, black), interface region of IMFR-1 and IMFR-2 (I12, green), the interior of IMFR-2 (R2, red), the interface of IMFR-2 and IMFR-3 (I23, violet), and the interior of IMFR-3 (R3, blue). The interval of each IMFR is indicated by horizontal arrows as in Fig. 1. (b, c) Time series of the skewness S and the kurtosis K computed using a sliding overlapping window with the same parameters as in Fig. 2. The SD computed in each window is represented by a gray area.

Figures 7 and 8 show the SK scatter plots for τ=100s and τ=10s, respectively. From these figures, we note that R2 does not display a parabolic shape on the two selected scales. The low value of the correlation index of R2 shown in Tables 4 and 5 confirms that the data points fit the parabolic shape poorly. This indicates that magnetic field fluctuations during R2 are nearly Gaussian even on the smallest scale.

Figure 7Kurtosis K as a function of skewness S computed from the time series of $\mathit{\delta }|\mathbit{B}|\left(\mathit{\tau }\right)$, where τ=100s. (a) The interior region of IMFR-1, (b) the interface of IMFR-1 and IMFR-2, (c) the interior of IMFR-2, (d) the interface of IMFR-2 and IMFR-3, and (e) the interior of IMFR-3. In each panel, the least-square fit with the parabolic function $K=\mathit{\alpha }{S}^{\mathrm{2}}+\mathit{\beta }$ is displayed as a dashed line (see Table 4).

Figure 8Kurtosis K as a function of skewness S computed from the time series of $\mathit{\delta }|\mathbit{B}|\left(\mathit{\tau }\right)$, where τ=10s. (a) The interior region of IMFR-1, (b) the interface of IMFR-1 and IMFR-2, (c) the interior of IMFR-2, (d) the interface of IMFR-2 and IMFR-3, and (e) the interior of IMFR-3. In each panel, the least-square fit with the parabolic function $K=\mathit{\alpha }{S}^{\mathrm{2}}+\mathit{\beta }$ is displayed as a dashed line (see Table 5).

Table 4The least-square fits of Eq. (7) computed from the scatter plots of Fig. 7 (τ=100s). The fitting function of IMFR-2 was not applicable (n/a) due to the small correlation value.

Table 5Same as in Table 3 for τ=10s. The fitting function of IMFR-2 was not applicable (n/a) due to the small correlation value.

Except for R2, all other regions show a parabolic shape at τ=100s that is enhanced at τ=10s, in agreement with the intermittent nature of magnetic field turbulence. Magnetic field fluctuations in the solar wind turbulence display a scale dependence in which they become intermittent as the scale becomes smaller, within the inertial subrange, due to rare, large-amplitude coherent structures. As a consequence, statistics of magnetic field fluctuations such as the PDFs of the ΔB (Fig. 5) departure from Gaussian statistics as τ decreases. By comparing the values of the correlation index shown in Table 4 for τ=100s with those of Table 5 for τ=10s we note that, for each region, the correlation r increases on the smallest scale, confirming that the SK parabolic relation displays scale dependence within the inertial subrange.

The highest correlation value for τ=100s (Table 4) corresponds to I23. This indicates that the ongoing magnetic reconnection occurring in this region can act as a source of non-Gaussianity and intermittent turbulence even on the largest scale. At τ=10s, Table 5 shows that r=0.99 at I23 and r=0.98 at R1. Small-scale current sheets were detected in these two intervals by and are responsible for intermittency and non-Gaussian fluctuations. Our result demonstrates that they are also responsible for the enhancement of the SK parabolic relation. Note that there are points in Fig. 8d that are further away from the (0, 0) Gaussian point, compared to Fig. 8a. This means that while the scatter plots of R1 and I23 are highly correlated with Eq. (7), the numerical values of S and K, which measure the degree of asymmetry and non-Gaussianity respectively, can be higher at I23.

5 Discussion

A theoretical explanation of the parabolic relation between the skewness and kurtosis of turbulent fluids and plasmas is still an open question. proposed a nonlinear Langevin equation with external forcing that can account for the parabolic relation between S and K. extended this model to include self-generated internal instabilities in plasmas. argued that a parabolic relation can be obtained as a natural consequence of a number of constraints expected to be met for most physical systems. proposed a simplified model of a synthetic intermittent time series, constructed from a random number of coherent structures with random amplitudes embedded in a background Gaussian noise, and demonstrated that their model can predict a SK parabolic relation. A similar study was performed by using a model of coherent plasma flux events.

Although a theoretical explanation of the SK relation is still unclear, there is a consensus that the parabolic shape is due to non-Gaussianity related to coherent structures, whereas points near $\left(S,K\right)=\left(\mathrm{0},\mathrm{0}\right)$ correspond to Gaussian fluctuations. This is confirmed by models of synthetic time series. For example, proposed a model of intermittent time series which consists of a superposition of Gaussian and non-Gaussian random fluctuations. Their model includes a parameter that measures the deviation from Gaussianity. The resulting PDF derived from their model displays asymmetric long tails that reproduce measured distributions of plasma density fluctuations in plasma magnetic confinement devices as well as distributions of X-ray emissions detected from accretion disks . Their model also leads to a parabolic relation between S and K. observed a transition from a parabolic shape to the $\left(S,K\right)=\left(\mathrm{0},\mathrm{0}\right)$ point by increasing the intensity of the Gaussian noise in their model of synthetic time series, constructed by adding deterministic fluctuations and Gaussian noise. However, a quantification of the parabolic shape is needed for an objective comparison between different datasets. We have found that the computation of the correlation index r allows time series dominated by either Gaussian and non-Gaussian fluctuations to be clearly distinguished. Despite the simplicity of this approach, it represents an alternative way to compare the degree of non-Gaussianity due to asymmetry and fat tails in the PDFs of different datasets, and can be applied to observational data and results from numerical simulations.

The stochastic model of a time series proposed by assumes that the non-Gaussian fluctuations arise from a quadratic nonlinear term. By increasing the degree of non-Gaussianity the skewness and the kurtosis converge to extreme values: $S=±\mathrm{2}\sqrt{\mathrm{2}}$ and K=12. This means that experimental data governed by nonlinear processes of quadratic order should lead to SK scatter plots with $S\in \left[-\mathrm{2}\sqrt{\mathrm{2}},\mathrm{2}\sqrt{\mathrm{2}}\right]$ and K<12. The scatter plots shown in Fig. 3a, c and e seem to agree with these limits; however, in Fig. 3d there are some points in which $S<-\mathrm{2}\sqrt{\mathrm{2}}$. also propose that processes described by higher-order nonlinearities can result in SK parabolic shapes with S outside the interval $\left[-\mathrm{2}\sqrt{\mathrm{2}},\mathrm{2}\sqrt{\mathrm{2}}\right]$, which can explain the behavior of $|\mathbit{B}|$ during the magnetic reconnection occurring in the I23 interval.

In the previous sections we showed and discussed the value of the correlation index measuring how well the SK scatter plots fit with a parabola. As mentioned before, there is no agreement on the interpretation of the coefficients α and β in Eq. (7). argues that the coefficients are not likely to offer relevant information about the underlying process. However, discussed an interpretation of the α and β coefficients based on their model of a synthetic time series. The value of the α coefficient depends on the statistics of the fluctuations due to coherent structures and is not necessarily constant in time. For the β coefficient, if the number of coherent structures in a time series can be represented as random independent variables that follow a Poisson distribution function (which models the occurrence of rare events), then β=3. Deviations from this value can be interpreted as a departure from the independence assumption, which means that there is interaction among coherent structures . Since we define kurtosis to be the flatness minus three, the previous statement is equivalent to say that deviations from β=0 are due to interacting coherent structures. From Table 3, we note that all intervals have nonzero values of β. Recall that this event is characterized by a rope–rope magnetic reconnection involving IMFR-2 and IMFR-3, with formation of a bifurcated current sheet acting as a source of intermittent turbulence . The interaction between the small-scale IMFR-2 and the medium-scale IMFR-3 occurring during this event gives support for the interpretation of the β parameter by .

demonstrated that the SK parabolic relation is also observed for time series of two-point differences of $|\mathbit{B}|$ in the solar wind. They showed that this relation is enhanced in the presence of large-scale events such as interplanetary shocks, whereas for nonshock intervals, the parabolic relation is not observed. In this case the SK parabolic relation represents a signature of direct coupling between large-scale structures (interplanetary shocks) and small-scale intermittency. Our results indicate that the SK parabolic relation is present during reconnection between a small-scale IMFR with a duration of ∼60min and a medium-scale IMFR with a duration of ∼7 h (see Table 1). The only region in which the parabolic relation is not observed is in the interior of IMFR-2. This region was found to have a low degree of intermittency and nearly monofractal scaling. Therefore, our results are in accordance with cross-scale coupling between IMFR scales and scales within the inertial subrange.

6 Conclusions

In this paper we investigated the relation between the skewness and the kurtosis during a triple-IMFR event on 2 February 2002. This event was divided into five regions, namely, the interior of IMFR-1, the interface of IMFR-1 and IMFR-2, the interior of IMFR-2, the interface of IMFR-2 and IMFR-3, and the interior of IMFR-3. We then computed the skewness S and the kurtosis K of $|\mathbit{B}|$ using a sliding window, and showed that the scatter plots of K as a function of S display a parabolic shape for all regions. The highest value of the correlation index computed by a least-square fit between the (S,K) values and Eq. (7) occurs at the interface of IMFR-2 and IMFR-3. This region was found to be the source of intermittent turbulence due to a magnetic reconnection between the small-size IMFR-2 and the medium-size IMFR-3 . Therefore, the enhanced SK parabolic relation is related to non-Gaussian fluctuations due to coherent structures emerging from intermittent turbulence generated via magnetic reconnection. The lowest value of the correlation index was obtained at the interior of IMFR-2, in agreement with the results of , who found that this region is characterized by a low degree of non-Gaussianity and phase synchronization, and nearly monofractal scaling.

We also analyzed the SK relation using two-point differences of $|\mathbit{B}|$ on two different scales within the inertial subrange. By computing the compensated PSD we selected an interval of frequencies in which all regions exhibit $-\mathrm{5}/\mathrm{3}$ scaling corresponding to the inertial subrange and selected two timescales representing the largest scale (τ=100s) and the smallest scale (τ=10s) within the inertial subrange. We found that the scatter plot of IMFR-2 on the largest scale (τ=100s) and on the smallest scale (τ=10s) accumulate around the $\left(S,K\right)=\left(\mathrm{0},\mathrm{0}\right)$ point. The least-square fit with Eq. (7) results in a low correlation index, which confirms that magnetic field fluctuations in this region are nearly Gaussian. All other regions displayed parabolic shapes. At τ=100s, the correlation index is high for the interface of IMFR-2 and IMFR-3, indicating that the magnetic reconnection that occurs in this region can generate non-Gaussian fluctuations on the largest scale. On the smallest scale, the correlation index is higher for two regions, namely, the interior of IMFR-1 and the interface of IMFR-2 and IMFR-3. This result can be due to non-Gaussian fluctuations resulting from small-scale current sheets detected within these regions . Our analysis indicates that the SK parabolic relation observed in interplanetary magnetic field turbulence is enhanced on small scales within the inertial subrange.

Our findings give support to the conclusion by that rope–rope magnetic reconnection acts as a source of interplanetary intermittent turbulence and suggest that magnetic reconnection is responsible for non-Gaussian PDFs with asymmetric shapes and fat tails. The results are also in agreement with the results of in that the SK parabolic relation is a signature of direct coupling between IMFR scales and small-scale intermittency.

Code and data availability
Code and data availability.

All data analyzed in this paper are publicly available via the Cluster Science Archive at http://www.cosmos.esa.int/web/csa (ESA, 2018). Numerical codes are also freely available at https://github.com/rmiracer (Miranda, 2018).

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Special issue statement
Special issue statement.

This article is part of the special issue “Space weather connections to near-Earth space and the atmosphere”. It is a result of the 6 Simpósio Brasileiro de Geofísica Espacial e Aeronomia (SBGEA), Jataí, Brazil, 26–30 September 2016.

Acknowledgements
Acknowledgements.

The authors are grateful to the reviewer for valuable comments. The authors would like to thank Heng Qiang Feng for providing the estimated times of the boundary layers for the three IMFRs observed by Cluster-1. Rodrigo A. Miranda acknowledges support from FAPDF (Brazil) under grant 0193.000984/2015. Adriane B. Schelin acknowledges support from FAPDF under grant 0193.000.884/2015. Abraham C.-L. Chian acknowledges the award of a PVE Distinguished Visiting Professor Fellowship by CAPES (grant no. 88881.068051/2014-01) and the hospitality of Erico Rempel of ITA. José L. Ferreira acknowledges support from the UNIESPAÇO program of the Brazilian Space Agency (AEB), the National Council of Technological and Scientific Development (CNPq), and FAPDF.
The topical editor, Alisson Dal Lago, thanks the two anonymous referees for help in evaluating this paper.

References

Antar, G. Y., Krasheninnikov, S. I., Devynck, P., Doerner, R. P., Hollmann, E. M., Boedo, J. A., Luckhardt, S. C., and Conn, R. W.: Experimental evidence of intermittent convection in the edge of magnetic confinement devices, Phys. Rev. Lett., 87, 065001, https://doi.org/10.1103/PhysRevLett.87.065001, 2001. a

Antar, G. Y., Counsell, G., Yu, Y., Labombard, B., and Devynck, P.: Universality of intermittent convective transport in the scrape-off layer of magnetically confined devices, Phys. Plasmas, 10, 419, https://doi.org/10.1063/1.1536166, 2003. a

Bale, S. D., Kellogg, P. J., Mozer, F. S., Horbury, T. S., and Rème, H.: Measurement of the electric fluctuation spectrum of magnetohydrodynamic turbulence, Phys. Rev. Lett., 94, 215002, https://doi.org/10.1103/PhysRevLett.94.215002, 2005.

Balogh, A., Carr, C. M., Acuña, M. H., Dunlop, M. W., Beek, T. J., Brown, P., Fornacon, K.-H., Georgescu, E., Glassmeier, K.-H., Harris, J., Musmann, G., Oddy, T., and Schwingenschuh, K.: The Cluster Magnetic Field Investigation: overview of in-flight performance and initial results, Ann. Geophys., 19, 1207–1217, https://doi.org/10.5194/angeo-19-1207-2001, 2001. a, b

Bard, Y.: Nonlinear Parameter Estimation, Academic Press, New York, 1974. a

Bergsaker, A. S., Fredriksen, Å., Pécseli, H. L., and Trulsen, J. K.: Models for the probability densities of the turbulent plasma flux in magnetized plasmas, Phys. Scripta, 90, 108005, https://doi.org/10.1088/0031-8949/90/10/108005, 2015. a, b, c, d

Bershadskii, A. and Sreenivasan, K. R.: Intermittency and the passive nature of the magnitude of the magnetic field, Phys. Rev. Lett., 93, 064501, https://doi.org/10.1103/PhysRevLett.93.064501, 2004. a

Biskamp, D., Schwarz, E., Zeiler, A., Celani, A., and Drake, J. F.: Electron magnetohydrodynamic turbulence, Phys. Plasmas, 6, 751, https://doi.org/10.1063/1.873312, 1999. a

Borovsky, J. E.: Flux tube texture of the solar wind: Strands of the magnetic carpet at 1 AU?, J. Geophys. Res., 113, A08110, https://doi.org/10.1029/2007JA012684, 2008. a, b

Bruno, R. and Carbone, V.: The solar wind as a turbulence laboratory, Living Rev. Sol. Phys., 2, 4, https://doi.org/10.12942/lrsp-2005-4, 2005.

Bruno, R., Carbone, V., Veltri, P., Pietropaolo, E., and Bavassano, B.: Identifying intermittency events in the solar wind, Planet. Space Sci., 49, 1201–1210, 2001. a, b, c

Bruno, R., Carbone, V., Bavassano, B., and Sorriso-Valvo, L.: Observations of magnetohydrodynamic turbulence in the 3-D heliosphere, Adv. Space Res., 35, 939–950, 2005.

Bruno, R., Carbone, V., Chapman, S., Hnat, B., Noullez, A., and Sorriso-Valvo, L.: Intermittent character of interplanetary magnetic field fluctuations, Phys. Plasmas, 14, 032901, https://doi.org/10.1063/1.2711429, 2007. a

Burlaga, L. F. and Viñas, A. F.: Multi-scale probability distributions of solar wind speed fluctuations at 1 AU described by a generalized Tsallis distribution, Geophys. Res. Lett., 31, L16807, https://doi.org/10.1029/2004GL020715, 2004.

Chian, A. C.-L. and Miranda, R. A.: Cluster and ACE observations of phase synchronization in intermittent magnetic field turbulence: a comparative study of shocked and unshocked solar wind, Ann. Geophys., 27, 1789–1801, https://doi.org/10.5194/angeo-27-1789-2009, 2009. a, b

Chian, A. C.-L., Feng, H. Q., Hu, Q., Loew, M. H., Miranda, R. A., Muñoz, P. R., Sibeck, D. G., and Wu, D. J.: Genesis of interplanetary intermittent turbulence: A case study of rope-rope magnetic reconnection, Astrophys. J., 832, 179, https://doi.org/10.3847/0004-637X/832/2/179, 2016. a, b, c, d, e, f, g, h, i, j, k, l, m

de Wit, T. D.: Can high-order moments be meaningfully estimated from experimental turbulence measurements?, Phys. Rev. E, 70, 055302, https://doi.org/10.1103/PhysRevE.70.055302, 2004. a

Eaton, J. W.: GNU Octave and reproducible research, J. Process. Contr., 22, 1433, https://doi.org/10.1016/j.jprocont.2012.04.006, 2012. a

Eaton, J. W., Bateman, D., Hauberg, S., and Wehbring, R.: GNU Octave version 3.8.1 manual: a high-level interactive language for numerical computations, CreateSpace Independent Publishing Platform, ISBN 441413006, 2014. a

ESA: Cluster Science Archive, available at: http://www.cosmos.esa.int/web/csa, last access: 19 March 2018.

Feng, H. Q., Wu, D. J., and Chao, J. K. J.: Size and energy distributions of interplanetary magnetic flux ropes, Geophys. Res., 112, A02102, https://doi.org/10.1029/2006JA011962, 2007. a

Greco, A., Chuychai, P., Matthaeus, W. H., Servidio, S., and Dmitruk, P.: Intermittent MHD structures and classical discontinuities, Geophys. Res. Lett., 35, L19111, https://doi.org/10.1029/2008GL035454, 2008. a

Greco, A., Matthaeus, W. H., Servidio, S., Chuychai, P., and Dmitruk, P.: Statistical analysis of discontinuities in solar wind ACE data and comparison with intermittent MHD turbulence, Astrophys. J., 691, L111–L114, 2009. a

Guszejnov, D., Lazányi, N., Bencze, A., and Zoletnik, S.: On the effect of intermittency of turbulence on the parabolic relation between skewness and kurtosis in magnetized plasmas, Phys. Plasmas, 20, 112305, https://doi.org/10.1063/1.4835535, 2013. a, b, c, d, e, f

Kamide, Y. and Chian, A. C.-L. (Eds.): Handbook of the Solar-Terrestrial Environment, Springer, Berlin, 2007.

Koga, D., Chian, A. C.-L., Miranda, R. A., and Rempel, E. L.: Intermittent nature of solar wind turbulence near the Earth's bow shock: phase coherence and non-Gaussianity, Phys. Rev. E, 75, 046401, https://doi.org/10.1103/PhysRevE.75.046401, 2007. a

Krommes, J. A.: The remarkable similarity between the scaling of kurtosis with squared skewness for TORPEX density fluctuations and sea-surface temperature fluctuations, Phys. Plasmas, 15, 030703, https://doi.org/10.1063/1.2894560, 2008. a, b

Labit, B., Furno, I., Fasoli, A., Diallo, A., Müller, S. H., Plyushchev, G., Podestà, M., and Poli, F. M.: Universal Statistical Properties of Drift-Interchange Turbulence in TORPEX Plasmas, Phys. Rev. Lett., 98, 255002, https://doi.org/10.1103/PhysRevLett.98.255002, 2007. a, b, c

Leamon, R. J., Smith, C. W., Ness, N. F., and Matthaeus, W. H.: Observational constraints on the dynamics of the interplanetary magnetic field dissipation range, J. Geophys. Res., 103, 4475–4787, 1998.

Lepping, R. P., Burlaga, L. F., Szabo, A., Ogilvie, K. W., Mish, W. H., Vassiliadis, D., Lazarus, A. J., Steinberg, J., Farrugia, C. J., Janoo, L., and Mariani, F.: The Wind magnetic cloud and events of October 18–20, 1995: Interplanetary properties and as triggers for geomagnetic activity, J. Geophys. Res., 102, 14049, https://doi.org/10.1029/97JA00272, 1997. a

Levenberg, K.: A method for the solution of certain non-linear problems in least squares, Q. Appl. Math., 2, 164–168, 1944. a

Marquardt, D. W.: An algorithm for least-squares estimation of nonlinear parameters, J. Soc. Ind. Appl. Math., 11, 2, https://doi.org/10.1137/0111030, 1963. a

Matthaeus, W. H. and Montgomery, D.: Selective decay hypothesis at high mechanical and magnetic Reynolds numbers, Ann. NY Acad. Sci., 357, 203–222, 1980. a

Matthaeus, W. H., Goldstein, M. L., and Smith, C.: Evaluation of magnetic helicity in homogeneous turbulence, Phys. Rev. Lett., 48, 1256–1259, 1982.

Medina, J. M. and Díaz, J. A.: Extreme reaction times determine fluctuation scaling in human color vision, Phys. A, 461, 125–132, 2016. a, b

Miranda, R. A.: Numerical tools for statistical analysis, available at: https://github.com/rmiracer, last access: 21 March 2018.

Miranda, R. A., Chian, A. C.-L., and Rempel, E. L.: Universal scaling laws for fully-developed magnetic field turbulence near and far upstream of the Earth's bow shock, Adv. Space Res., 51, 1893–1901, 2013. a, b

Moldwin, M. B., Ford, S., Lepping, R., Slavin, J., and Szabo, A.: Small-scale magnetic flux ropes in the solar wind, Geophys. Res. Lett., 27, 57, https://doi.org/10.1029/1999GL010724, 2000. a

Mole, N. and Clarke, E. D.: Relationships between higher moments of concentration and of dose in turbulent dispersion, Bound.-Lay. Meteorol., 73, 35–52, 1995. a, b

Narita, Y., Glassmeier, K.-H., and Treumann, R. A.: Wave-number spectra and intermittency in the terrestrial foreshock region, Phys. Rev. Lett., 97, 191101, https://doi.org/10.1103/PhysRevLett.97.191101, 2006.

Rème, H., Aoustin, C., Bosqued, J. M., Dandouras, I., Lavraud, B., Sauvaud, J. A., Barthe, A., Bouyssou, J., Camus, Th., Coeur-Joly, O., Cros, A., Cuvilo, J., Ducay, F., Garbarowitz, Y., Medale, J. L., Penou, E., Perrier, H., Romefort, D., Rouzaud, J., Vallat, C., Alcaydé, D., Jacquey, C., Mazelle, C., d'Uston, C., Möbius, E., Kistler, L. M., Crocker, K., Granoff, M., Mouikis, C., Popecki, M., Vosbury, M., Klecker, B., Hovestadt, D., Kucharek, H., Kuenneth, E., Paschmann, G., Scholer, M., Sckopke, N., Seidenschwang, E., Carlson, C. W., Curtis, D. W., Ingraham, C., Lin, R. P., McFadden, J. P., Parks, G. K., Phan, T., Formisano, V., Amata, E., Bavassano-Cattaneo, M. B., Baldetti, P., Bruno, R., Chionchio, G., Di Lellis, A., Marcucci, M. F., Pallocchia, G., Korth, A., Daly, P. W., Graeve, B., Rosenbauer, H., Vasyliunas, V., McCarthy, M., Wilber, M., Eliasson, L., Lundin, R., Olsen, S., Shelley, E. G., Fuselier, S., Ghielmetti, A. G., Lennartsson, W., Escoubet, C. P., Balsiger, H., Friedel, R., Cao, J.-B., Kovrazhkin, R. A., Papamastorakis, I., Pellat, R., Scudder, J., and Sonnerup, B.: First multispacecraft ion measurements in and near the Earth's magnetosphere with the identical Cluster ion spectrometry (CIS) experiment, Ann. Geophys., 19, 1303–1354, https://doi.org/10.5194/angeo-19-1303-2001, 2001.  a

Russell, C. T. and Elphic, R. C.: Observation of magnetic flux ropes in the Venus ionosphere, Nature, 279, 616, https://doi.org/10.1038/279616a0, 1979. a

Sandberg, I., Benkadda, S., Garbet, X., Ropokis, G., Hizanidis, K., and del-Castillo-Negrete, D.: Universal probability distribution function for bursty transport in plasma turbulence, Phys. Rev. Lett., 103, 165001, https://doi.org/10.1103/PhysRevLett.103.165001, 2009. a, b, c, d, e, f

Sattin, F., Agostini, M., Cavazzana, R., Serianni, G., Scarin, P., and Vianello, N.: About the parabolic relation existing between the skewness and the kurtosis in time series of experimental data, Phys. Scripta, 79, 045006, https://doi.org/10.1088/0031-8949/79/04/045006, 2009. a, b, c, d

Sorriso-Valvo, L., Carbone, V., Giuliani, P., Veltri, P., Bruno, R., Antoni, V., and Martines, E.: Intermittency in plasma turbulence, Planet. Space Sci., 49, 1193–1200, 2001. a

Sura, P. and Sardeshmukh, P. D.: A Global View of Non-Gaussian SST Variability, J. Phys. Oceanogr., 38, 638, https://doi.org/10.1175/2007JPO3761.1, 2007. a, b

Telloni, D., Carbone, V., Perri, S., Bruno, R., Lepreti, F., and Veltri, P.: Relaxation processes within flux ropes in solar wind, Astrophys. J., 826, 205, https://doi.org/10.3847/0004-637X/826/2/205, 2016. a, b

Veltri, P.: MHD turbulence in the solar wind: self-similarity, intermittency and coherent structures, Plasma Phys. Contr. F., 41, A787–A795, 1999. a

Vörös, Z., Leubner, M. P., and Baumjohann, W. J.: Cross-scale coupling-induced intermittency near interplanetary shocks, J. Geophys. Res., 111, A02102, https://doi.org/10.1002/2015JA021257, 2006. a, b, c, d, e

Vörös, Z., Baumjohann, W., Nakamura, R., Runov, A., Volwerk, M., Takada, T., Lucek, E. A., and Rème, H.: Spatial structure of plasma flow associated turbulence in the Earth's plasma sheet, Ann. Geophys., 25, 13–17, https://doi.org/10.5194/angeo-25-13-2007, 2007.

Wei, F. S., Liu, R., Fan, Q., and Feng, X. S.: Identification of the magnetic cloud boundary layers, J. Geophys. Res., 108, A1263, https://doi.org/10.1029/2002JA009511, 2003. a