the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
NonGaussianity and crossscale coupling in interplanetary magnetic field turbulence during a rope–rope magnetic reconnection event
Adriane B. Schelin
Abraham C.L. Chian
José L. Ferreira
In a recent paper (Chian et al., 2016) 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 thirdorder moment (skewness) and the normalized fourthorder 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 nonGaussian 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 $\left\mathit{B}\right$ during a triple interplanetary magnetic flux rope event. In addition, we investigate the skewness–kurtosis relation of twopoint differences of $\left\mathit{B}\right$ 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 nonGaussian 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)
The solar wind can be regarded as a network of entangled magnetic flux tubes and Alfvénic fluctuations propagating within each flux tube (Bruno et al., 2001; Borovsky, 2008). Flux tubes can emerge locally in the solar wind as a consequence of the magnetohydrodynamic turbulent cascade (Matthaeus and Montgomery, 1980; Veltri, 1999; Greco et al., 2008, 2009; Telloni et al., 2016). An alternative view describes coherent structures as “fossile” structures that emanate from the solar surface and are advected by the solar wind (Borovsky, 2008; Bruno et al., 2001).
The probability distribution functions (PDFs) of turbulent space plasmas display sharp peaks and fat tails on small scales within the inertial subrange (SorrisoValvo et al., 2001; Bruno et al., 2001; Koga et al., 2007; Chian and Miranda, 2009), as well as departures from selfsimilarity and monofractality (Bershadskii and Sreenivasan, 2004; Bruno et al., 2007; Miranda et al., 2013). These features are due to the presence of rare, largeamplitude 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 thirdorder moment (skewness) and the normalized fourthorder moment (kurtosis) has been found in local concentrations of contaminants in atmospheric turbulence as found by Mole and Clarke (1995). Sura and Sardeshmukh (2007) also found a similar skewness–kurtosis parabolic relation using global data of seasurface temperature fluctuations. Labit et al. (2007) reported a similar skewness–kurtosis dependence in electron density fluctuations in plasma confinement experiments. Medina and Díaz (2016) 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 (Krommes, 2008; Sattin et al., 2009; Sandberg et al., 2009; Guszejnov et al., 2013; Bergsaker et al., 2015) and has been associated with the presence of nonGaussian fluctuations due to coherent structures (Labit et al., 2007; Sandberg et al., 2009; Guszejnov et al., 2013; Bergsaker et al., 2015). The skewness–kurtosis parabolic relation was also found in time series of twopoint differences of the modulus of the magnetic field by Vörös et al. (2006). They demonstrated that the parabolic relation is due to nonlocal interaction between largescale structures and smallscale intermittency.
In this paper we investigate the skewness–kurtosis relation during a triple interplanetary magnetic flux rope (IMFR) event detected by Cluster1 in the solar wind. This event was recently characterized by Chian et al. (2016). 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 tripleIMFR 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.
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 $\left\mathit{B}\right$). The skewness of θ_{i} can be computed as follows:
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
from which the kurtosis can be obtained by
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” (Sattin et al., 2009).
A common way to characterize asymmetry and nonGaussianity of θ_{i} as a function of scale τ is through the time series of twopoint differences:
The skewness of θ_{i} on scale τ is then
and the flatness is
where σ_{τ} is the SD of δθ_{i}(τ). From Eq. (5) the kurtosis as a function of scale is obtained by
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 (e.g., Mole and Clarke, 1995; Labit et al., 2007; Medina and Díaz, 2016). This relation is given by
where α and β are the coefficients that characterize a parabolic curve.
We compute the α and β coefficients by applying a leastsquare fit between (S,K) values obtained from the observational data and Eq. (7) following the Levenberg–Marquardt algorithm (Levenberg, 1944; Marquardt, 1963; Bard, 1974), 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 X_{i} and Y_{i}, $i=\mathrm{1},\mathrm{\dots},N$:
where σ_{X} and σ_{Y} represent the SD of X_{i} and Y_{i}, respectively. The correlation index $r\in [\mathrm{1},\mathrm{1}]$. If r=1 there is complete correlation between X_{i} and Y_{i}, 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 $\left\mathit{B}\right$ 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 twopoint 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 (Eaton, 2012; Eaton et al., 2014).
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).
Figure 1a shows the time series of the modulus of magnetic field $\left\mathit{B}\right$ obtained by the FGM instrument onboard Cluster1 (Balogh et al., 2001) from 00:00 to 12:00 UT on 2 February 2002. During this interval Cluster1 was in the solar wind upstream of the Earth's bow shock (Chian and Miranda, 2009). The magnetic field data are collected by Cluster1 at a resolution of 22 Hz (Balogh et al., 2001). 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 $\left{\mathbf{V}}_{\mathrm{i}}\right$, the ion number density n_{i}, the ion temperature perpendicular to the magnetic field T_{i} and the ion plasma β_{i}, which is the ratio between plasma kinetic pressure and magnetic pressure. The Cluster1 plasma measurements are given by the ion spectrometry experiment CIS (Rème et al., 2001).
This event is characterized by the presence of three interplanetary magnetic flux ropes. Magnetic flux ropes are magnetic structures described as bundles of twisted, currentcarrying magnetic field lines bent into a tubelike shape, spiralling around a common axis (Russell and Elphic, 1979; Telloni et al., 2016; Chian et al., 2016). During this event three IMFRs were identified by Chian et al. (2016) using a combination of criteria for largescale magnetic cloud boundary layers (Lepping et al., 1997; Wei et al., 2003) and smallscale IMFRs (Moldwin et al., 2000; Feng et al., 2007). The interval of each IMFR is indicated by horizontal arrows in Fig. 1a, and their timings are shown in Table 1.
4.1 Time series of $\left\mathit{B}\right$
Figure 2a shows the time series of $\left\mathit{B}\right$ detected by Cluster1 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 IMFR1 (R_{1}), the interface of IMFR1 and IMFR2 (I_{12}), the interior of IMFR2 (R_{2}), the interface of IMFR2 and IMFR3 (I_{23}), and the interior of IMFR3 (R_{3}). 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 IMFR1 and at the interface region between IMFR2 and IMFR3. This interface region was identified as a source of intermittent turbulence by Chian et al. (2016). A current sheet was detected at the leading edge of IMFR1 using data from ACE and Cluster1, and a current sheet was detected at the interface region between IMFR2 and IMFR3 using data from Cluster1, ACE and Wind (Chian et al., 2016).
The S–K 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 (de Wit, 2004; Miranda et al., 2013). 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 $\left\mathit{B}\right$ at the interface between IMFR2 and IMFR3. A similar behavior was observed in magnetic field data during an interplanetary shock event by Vörös et al. (2006). The uncertainty within sliding windows that contain the large variations in $\left\mathit{B}\right$ increases due to nonstationarity. Following Vörös et al. (2006), we exclude these windows from further analysis.
Figure 3 shows K as a function of S for the five regions previously defined. A leastsquare 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.
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 R_{2}, in agreement with a visual inspection of Fig. 3c. For this interval, most of the points in Fig. 3c tend to accumulate around $(S,K)=(\mathrm{0},\mathrm{0})$, which is the value obtained for a Gaussian distribution (i.e., in the absence of coherent structures). Therefore, the interior of IMFR2 is characterized by a low degree of nonGaussianity and intermittency in comparison with the other intervals.
The highest value of the correlation is obtained during I_{23} (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 nonGaussianity. These results are in agreement with the results of Chian et al. (2016), which found that the interior of IMFR2 has lower degrees of nonGaussianity and phase coherence, and a nearly monofractal scaling when compared with other intervals. For the interface of IMFR2 and IMFR3 they observed higher degrees of nonGaussianity and phase synchronization, and a strong departure from monofractality.
4.2 Time series of $\mathit{\delta}\left\mathit{B}\right$
Next, we investigate the S–K 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 $\left\mathit{B}\right$ from the beginning of IMFR1 at 00:32 UT until the end of IMFR3 at 08:40 UT. The right side of Fig. 4a shows the compensated PSD which is the original PSD multiplied by f^{5∕3} (Biskamp et al., 1999). 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 R_{1}, I_{12}, R_{2}, I_{23} and R_{3}. 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.01 Hz and ends at f=0.1 Hz, which correspond to scales τ=100 s and τ=10 s, respectively.
The intermittent aspect of interplanetary magnetic field turbulence can be demonstrated by constructing the PDF of the normalized magneticfield differences
where $\mathit{\delta}B\left(\mathit{\tau}\right)=\left\mathit{B}\right(t+\mathit{\tau}\left)\right\left\mathit{B}\right(t\left)\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 τ=10 s and τ=100 s. From this figure it is clear that the PDFs are closer to a Gaussian distribution (represented by the gray area in Fig. 5) at τ=100 s (large scale), and become nonGaussian at τ=10 s (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 S–K relation of ΔB(τ) at τ=10 s and τ=100 s. 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 IMFR2 and IMFR3; therefore these S and K values are excluded from further analysis.
Figures 7 and 8 show the S–K scatter plots for τ=100 s and τ=10 s, respectively. From these figures, we note that R_{2} does not display a parabolic shape on the two selected scales. The low value of the correlation index of R_{2} shown in Tables 4 and 5 confirms that the data points fit the parabolic shape poorly. This indicates that magnetic field fluctuations during R_{2} are nearly Gaussian even on the smallest scale.
Except for R_{2}, all other regions show a parabolic shape at τ=100 s that is enhanced at τ=10 s, 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, largeamplitude 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 τ=100 s with those of Table 5 for τ=10 s we note that, for each region, the correlation r increases on the smallest scale, confirming that the S–K parabolic relation displays scale dependence within the inertial subrange.
The highest correlation value for τ=100 s (Table 4) corresponds to I_{23}. This indicates that the ongoing magnetic reconnection occurring in this region can act as a source of nonGaussianity and intermittent turbulence even on the largest scale. At τ=10 s, Table 5 shows that r=0.99 at I_{23} and r=0.98 at R_{1}. Smallscale current sheets were detected in these two intervals by Chian et al. (2016) and are responsible for intermittency and nonGaussian fluctuations. Our result demonstrates that they are also responsible for the enhancement of the S–K 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 R_{1} and I_{23} are highly correlated with Eq. (7), the numerical values of S and K, which measure the degree of asymmetry and nonGaussianity respectively, can be higher at I_{23}.
A theoretical explanation of the parabolic relation between the skewness and kurtosis of turbulent fluids and plasmas is still an open question. Sura and Sardeshmukh (2007) proposed a nonlinear Langevin equation with external forcing that can account for the parabolic relation between S and K. Krommes (2008) extended this model to include selfgenerated internal instabilities in plasmas. Sattin et al. (2009) 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. Guszejnov et al. (2013) 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 S–K parabolic relation. A similar study was performed by Bergsaker et al. (2015) using a model of coherent plasma flux events.
Although a theoretical explanation of the S–K relation is still unclear, there is a consensus that the parabolic shape is due to nonGaussianity related to coherent structures, whereas points near $(S,K)=(\mathrm{0},\mathrm{0})$ correspond to Gaussian fluctuations. This is confirmed by models of synthetic time series. For example, Sandberg et al. (2009) proposed a model of intermittent time series which consists of a superposition of Gaussian and nonGaussian 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 (Antar et al., 2001, 2003) as well as distributions of Xray emissions detected from accretion disks (Sandberg et al., 2009). Their model also leads to a parabolic relation between S and K. Bergsaker et al. (2015) observed a transition from a parabolic shape to the $(S,K)=(\mathrm{0},\mathrm{0})$ 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 nonGaussian fluctuations to be clearly distinguished. Despite the simplicity of this approach, it represents an alternative way to compare the degree of nonGaussianity 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 Sandberg et al. (2009) assumes that the nonGaussian fluctuations arise from a quadratic nonlinear term. By increasing the degree of nonGaussianity the skewness and the kurtosis converge to extreme values: $S=\pm \mathrm{2}\sqrt{\mathrm{2}}$ and K=12. This means that experimental data governed by nonlinear processes of quadratic order should lead to S–K scatter plots with $S\in [\mathrm{2}\sqrt{\mathrm{2}},\mathrm{2}\sqrt{\mathrm{2}}]$ 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}}$. Sandberg et al. (2009) also propose that processes described by higherorder nonlinearities can result in S–K parabolic shapes with S outside the interval $[\mathrm{2}\sqrt{\mathrm{2}},\mathrm{2}\sqrt{\mathrm{2}}]$, which can explain the behavior of $\left\mathit{B}\right$ during the magnetic reconnection occurring in the I_{23} interval.
In the previous sections we showed and discussed the value of the correlation index measuring how well the S–K scatter plots fit with a parabola. As mentioned before, there is no agreement on the interpretation of the coefficients α and β in Eq. (7). Sattin et al. (2009) argues that the coefficients are not likely to offer relevant information about the underlying process. However, Guszejnov et al. (2013) 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 (Guszejnov et al., 2013). 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 IMFR2 and IMFR3, with formation of a bifurcated current sheet acting as a source of intermittent turbulence (Chian et al., 2016). The interaction between the smallscale IMFR2 and the mediumscale IMFR3 occurring during this event gives support for the interpretation of the β parameter by Guszejnov et al. (2013).
Vörös et al. (2006) demonstrated that the S–K parabolic relation is also observed for time series of twopoint differences of $\left\mathit{B}\right$ in the solar wind. They showed that this relation is enhanced in the presence of largescale events such as interplanetary shocks, whereas for nonshock intervals, the parabolic relation is not observed. In this case the S–K parabolic relation represents a signature of direct coupling between largescale structures (interplanetary shocks) and smallscale intermittency. Our results indicate that the S–K parabolic relation is present during reconnection between a smallscale IMFR with a duration of ∼60 min and a mediumscale 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 IMFR2. This region was found to have a low degree of intermittency and nearly monofractal scaling. Therefore, our results are in accordance with crossscale coupling between IMFR scales and scales within the inertial subrange.
In this paper we investigated the relation between the skewness and the kurtosis during a tripleIMFR event on 2 February 2002. This event was divided into five regions, namely, the interior of IMFR1, the interface of IMFR1 and IMFR2, the interior of IMFR2, the interface of IMFR2 and IMFR3, and the interior of IMFR3. We then computed the skewness S and the kurtosis K of $\left\mathit{B}\right$ 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 leastsquare fit between the (S,K) values and Eq. (7) occurs at the interface of IMFR2 and IMFR3. This region was found to be the source of intermittent turbulence due to a magnetic reconnection between the smallsize IMFR2 and the mediumsize IMFR3 (Chian et al., 2016). Therefore, the enhanced S–K parabolic relation is related to nonGaussian 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 IMFR2, in agreement with the results of Chian et al. (2016), who found that this region is characterized by a low degree of nonGaussianity and phase synchronization, and nearly monofractal scaling.
We also analyzed the S–K relation using twopoint differences of $\left\mathit{B}\right$ 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 (τ=100 s) and the smallest scale (τ=10 s) within the inertial subrange. We found that the scatter plot of IMFR2 on the largest scale (τ=100 s) and on the smallest scale (τ=10 s) accumulate around the $(S,K)=(\mathrm{0},\mathrm{0})$ point. The leastsquare 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 τ=100 s, the correlation index is high for the interface of IMFR2 and IMFR3, indicating that the magnetic reconnection that occurs in this region can generate nonGaussian fluctuations on the largest scale. On the smallest scale, the correlation index is higher for two regions, namely, the interior of IMFR1 and the interface of IMFR2 and IMFR3. This result can be due to nonGaussian fluctuations resulting from smallscale current sheets detected within these regions (Chian et al., 2016). Our analysis indicates that the S–K 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 Chian et al. (2016) that rope–rope magnetic reconnection acts as a source of interplanetary intermittent turbulence and suggest that magnetic reconnection is responsible for nonGaussian PDFs with asymmetric shapes and fat tails. The results are also in agreement with the results of Vörös et al. (2006) in that the S–K parabolic relation is a signature of direct coupling between IMFR scales and smallscale intermittency.
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).
The authors declare that they have no conflict of interest.
This article is part of the special issue “Space weather connections to nearEarth 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.
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
Cluster1. 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/201401) 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.
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 scrapeoff 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 inflight performance and initial results, Ann. Geophys., 19, 1207–1217, https://doi.org/10.5194/angeo1912072001, 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/00318949/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/lrsp20054, 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 SorrisoValvo, L.: Observations of magnetohydrodynamic turbulence in the 3D heliosphere, Adv. Space Res., 35, 939–950, 2005.
Bruno, R., Carbone, V., Chapman, S., Hnat, B., Noullez, A., and SorrisoValvo, 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.: Multiscale 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/angeo2717892009, 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 roperope magnetic reconnection, Astrophys. J., 832, 179, https://doi.org/10.3847/0004637X/832/2/179, 2016. a, b, c, d, e, f, g, h, i, j, k, l, m
de Wit, T. D.: Can highorder 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 highlevel 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 SolarTerrestrial 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 nonGaussianity, 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 seasurface 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 DriftInterchange 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 nonlinear problems in least squares, Q. Appl. Math., 2, 164–168, 1944. a
Marquardt, D. W.: An algorithm for leastsquares 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 fullydeveloped 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.: Smallscale 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.: Wavenumber 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., CoeurJoly, 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., BavassanoCattaneo, 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/angeo1913032001, 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 delCastilloNegrete, 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/00318949/79/04/045006, 2009. a, b, c, d
SorrisoValvo, 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 NonGaussian 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/0004637X/826/2/205, 2016. a, b
Veltri, P.: MHD turbulence in the solar wind: selfsimilarity, intermittency and coherent structures, Plasma Phys. Contr. F., 41, A787–A795, 1999. a
Vörös, Z., Leubner, M. P., and Baumjohann, W. J.: Crossscale couplinginduced 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/angeo25132007, 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