# On the Statistical Properties of Turbulent Energy Transfer Rate in the Inner Heliosphere

Luca Sorriso-Valvo, Francesco Carbone, Silvia Perri, Antonella Greco, Raffaele Marino, Roberto Bruno

► To cite this version:

Luca Sorriso-Valvo, Francesco Carbone, Silvia Perri, Antonella Greco, Raffaele Marino, et al. On the Statistical Properties of Turbulent Energy Transfer Rate in the Inner Heliosphere. Solar Physics, 2018, 293 (1), pp.10.1007/s11207-017-1229-6. hal-03866844

HAL Id: hal-03866844

https://hal.science/hal-03866844v1

Submitted on 20 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# On the Statistical Properties of Turbulent Energy Transfer Rate in the Inner Heliosphere

Luca Sorriso-Valvo<fnref n="1" /> · Francesco Carbone<fnref n="2" />

Silvia Perri<fnref n="3" /> · Antonella Greco<fnref n="3" /> · Raffaele Marino<fnref n="4" />

Roberto Bruno<fnref n="5" />

**Abstract** The transfer of energy from large to small scales in solar wind turbulence is an important ingredient of the long-standing question of the mechanism of the interplanetary plasma heating. Previous studies have shown that magnetohydrodynamic (MHD) turbulence is statistically compatible with the observed solar wind heating as it expands in the heliosphere. However, in order to understand which processes contribute to the plasma heating, it is necessary to have a local description of the energy flux across scales. To this aim, it is customary to use indicators such as the magnetic field partial variance of increments (PVI), which is associated with the local, relative, scale-dependent magnetic energy. A more complete evaluation of the energy transfer should also include other terms, related to velocity and cross-helicity. This is achieved here by introducing a proxy for the local, scale-dependent turbulent energy transfer rate $\epsilon_{\Delta t}(t)$, based on the third-order moment scaling law for MHD turbulence. Data from *Helios 2* are used to determine the statistical properties of such a proxy in comparison with the magnetic and velocity fields PVI, and the correlation with local solar wind heating is computed. PVI and $\epsilon_{\Delta t}(t)$ are generally well correlated; however, $\epsilon_{\Delta t}(t)$ is a very sensitive proxy that can exhibit large amplitude values, both positive and negative, even for low amplitude peaks in the PVI. Furthermore, $\epsilon_{\Delta t}(t)$ is very well correlated with local increases of the temperature when large amplitude bursts of energy transfer are localized, thus suggesting an important role played by this proxy in the study of plasma energy dissipation.

**Keywords** Solar wind · Turbulence · Intermittency

<figure><img src="image_2.png" /></figure>

<footnote>L. Sorriso-Valvo
lucasorriso@gmail.com</footnote>

<footnote n="1">1 Nanotec – UOS di Cosenza, CNR, Ponte P. Bucci, Cubo 31C, 87036, Rende, Italy</footnote>

<footnote n="2">2 IIA – U.O.S. di Cosenza, CNR, c/o: Unical-Polifunzionale, 87036 Rende, Italy</footnote>

<footnote n="3">3 Dipartimento di Fisica, Università della Calabria, 87036 Rende (CS), Italy</footnote>

<footnote n="4">4 Laboratoire de Mécanique des Fluides et d'Acoustique, École Centrale de Lyon, CNRS – Université de Lyon, 69134 Écully, France</footnote>

<footnote n="5">5 IAPS-INAF, via Fosso del Cavaliere 100, 00133 Roma, Italy</footnote>
# 1. Introduction

The solar wind is the most important example of natural plasma turbulence that can be probed using satellite-born instrumentation (Tu and Marsch, 1995; Bruno and Carbone, 2013). This means space missions dedicated to *in situ* measurements of plasma parameters and electromagnetic fields provide a unique chance to obtain direct experimental observation of the turbulent dynamics of space plasmas. Solar wind plasma has a very low density (of the order of a few particles per cubic centimeter) and high temperature (of the order of $10^5$ K), and is embedded in a radially decreasing background magnetic field of the order of $B_{SW} \sim 5$ nT near the Earth, resulting in a weakly collisional, magnetized flow. Strong acceleration mechanisms push the wind away from the Sun to a typical speed $v_{SW} \sim 350-750$ km s$^{-1}$, making the wind supersonic and superalfvénic. Observations also show that the temperature decreases with the distance from the Sun more slowly than expected for an adiabatically expanding plasma (Schwenn, 1983; Freeman, Totten, and Arya, 1992; Goldstein et al., 1996). Understanding the heating mechanism providing the non-adiabatic cooling of the expanding solar wind is a long-standing open question in astrophysics (Matthaeus et al., 1994; Richardson et al., 1995). One of the possible sources of heat is the dissipation of the kinetic and magnetic energy available in the form of large-scale fluctuations, which can be traced back to the Sun and the solar corona. This requires that the energy is transported from large scale to smaller scales by a turbulent cascade, where kinetic plasma processes can convert it into particle heating. Observations of magnetic and velocity power spectral density have shown that solar wind fluctuations follow the typical Kolmogorov-like power-law energy decay (Kolmogorov, 1941; Frisch, 1995; Marsch and Tu, 1997; Bruno and Carbone, 2013), in a range between the typical correlation length (corresponding to a few hours) (Matthaeus et al., 2005), and the typical scale where kinetic processes arise (a few seconds) (Leamon et al., 2000). This gives about three decades of inertial range where non-linear energy transfer occurs, resulting in high Reynolds number turbulence (Sorriso-Valvo et al., 2006; Weygand et al., 2007; Sorriso-Valvo, Yordanova, and Carbone, 2010). Turbulence is therefore the major ingredient to make a connection between the large-scale fluctuations and the small-scale microphysics processes (Alexandrova et al., 2013).

While at spatial scales smaller than the typical ion scales, kinetic processes must be included in the dynamics, at larger scales the magnetohydrodynamics (MHD) framework (Biskamp, 1997) is a good approximation to describe the turbulent motion of the solar wind plasma. The properties of solar wind turbulence have been studied for more than 40 years using both experimental data and, more recently, numerical simulations. It is now understood that the spectral properties of the magnetic and velocity fluctuations depend on the wind speed, distance from the Sun, solar activity, correlation between velocity and magnetic field, and other local plasma parameters, making the solar wind a complex environment with high variability (Bruno and Carbone, 2013). Solar wind turbulence is also characterized by anisotropy and intermittency, which have been deeply studied in the past (Bavassano et al., 1982; Marsch and Tu, 1997). Intermittency, in particular, is related to the appearance of small-scale structures typical of turbulence. In the solar wind, these are mostly current sheets, magnetic discontinuities, vorticity structures or similar features, and they are usually identified using field increments or wavelet-based detection techniques (Veltri and Mangeney, 1999; Bruno et al., 2001; Greco et al., 2009; Greco and Perri, 2014; Zhdankin, Boldyrev, and Chen, 2016). One example of identification technique recently introduced is the partial variance of increments (PVI), based on the evaluation of the intensity of the field gradients at a given position and scale (Greco et al., 2009; Greco and Perri, 2014). Recent investigations have confirmed that the intermittent structures are associated with enhanced
plasma heating. Both ions (Osman et al., 2012; Tessein et al., 2013) and electrons (Chasapis et al., 2015) display energization in the proximity of the most intense current sheets. This has also been confirmed through the analysis of numerical simulations of the Vlasov-Maxwell equations (Servidio et al., 2012). The processes responsible for such heating may involve magnetic reconnection, plasma instabilities and enhancement of collisions, and they are still poorly understood (Chen, 2016). In this work, a different data analysis technique is proposed to identify the regions of space that are carrying energy towards the small scales, in order to understand the link between the presence of strong turbulent fluctuations and a local enhancement of the temperature. In analogy with the von Karman-Howart law for Navier-Stokes turbulence, the MHD turbulent energy flux across the scales is regulated by a relation, often referred to as the Politano-Pouquet law (PP) (Politano and Pouquet, 1998). This is a statistical prescription for the scaling law of the mixed third-order moment of the Elsasser fields increments, and is obtained directly from the MHD equations under the assumptions of stationarity, isotropy, incompressibility and vanishing dissipation coefficients (i.e. within the inertial range). Although some of the above assumptions are only marginally satisfied in the solar wind, the validity of the PP law has been successfully verified in numerical simulations (Sorriso-Valvo et al., 2002) and in the solar wind (MacBride, Forman, and Smith, 2005; Sorriso-Valvo et al., 2007; Marino et al., 2008; MacBride, Smith, and Forman, 2008; Marino et al., 2012). Subsequently, a variety of extensions of the PP law to more complex, realistic systems have been introduced, where the approximations of incompressibility (Carbone et al., 2009; Banerjee and Galtier, 2013) and isotropy (Stawarz et al., 2011; Osman et al., 2011) have been relaxed, and in some cases verified in experimental data (Banerjee et al., 2016; Hadid, Sahraoui, and Galtier, 2017). The effects of the solar wind expansion have also been considered (Gogoberidze, Perri, and Carbone, 2013; Hellinger et al., 2013), and an attempt to include the small-scale effects described by Hall MHD has been performed (Galtier, 2008). In this work, we chose not to take into account any of the above modifications, in order to provide a first-order estimation of the energy transfer rate. The evaluation of the contributions emerging when approximations are relaxed represents an interesting possible improvement that we leave to future work.

The basic version of the PP law reads

$$Y^{\pm}(\Delta t) = \langle |\Delta z_{\Delta t}^{\pm}(t)|^2 \Delta z_{\parallel, \Delta t}^{\mp}(t) \rangle = -\frac{4}{3} \langle \epsilon^{\pm} \rangle \Delta t \langle v \rangle. \quad (1)$$

Here $\Delta z_{\Delta t}^{\pm}(t) = z^{\pm}(t + \Delta t) - z^{\pm}(t)$ indicates the increment of the field $z^{\pm}$ across a temporal scale $\Delta t$. The same notation applies to the other fields used in this work. The subscript $\parallel$ indicates the longitudinal component, i.e. parallel to the bulk speed in solar wind time series; $z^{\pm} = v \pm B/\sqrt{4\pi\rho}$ are the Elsasser variables that couple the solar wind velocity $v$ and the magnetic field $B$, transformed in velocity units using the solar wind density $\rho$; $Y^{\pm}(\Delta t)$ are the mixed third-order moments, and $\langle \epsilon^{\pm} \rangle$ is the mean energy transfer rate, estimated over the whole domain. In order to study spacecraft time series, all spatial scales $\ell$ were customarily transformed in the time lags $\Delta t = \ell/|\langle v \rangle|$ through the bulk flow speed $\langle v \rangle$ averaged over the entire dataset. This is allowed by the Taylor hypothesis (Taylor, 1938), which is robustly valid for solar wind fluctuations in the inertial range (Perri et al., 2017). The PP law in Equation 1 thus indicates that the nonlinear transport of energy across the time scales is proportional to the time scale via the mean energy transfer rate.

In MHD numerical simulations, the statistical properties of the local energy dissipation can be studied directly (Zhdankin, Boldyrev, and Chen, 2016). However, when the plasma is weakly collisional, MHD viscous and resistive dissipative terms are not defined. In such cases, if a sufficient scale separation exists between the inertial range and the dissipative
scales, even though the dissipation mechanisms are unknown the PP law can provide an estimate of the mean energy transfer rate. This has recently been measured from solar wind data, providing values compatible with the energy necessary for the observed non-adiabatic cooling (Marino et al., 2008; Smith et al., 2009; Carbone et al., 2009; Marino et al., 2011; Coburn et al., 2012; Banerjee et al., 2016; Hadid, Sahraoui, and Galtier, 2017).

Although the PP law is only valid in a statistical sense, dimensional considerations suggest that it may be possible to use its local values as a proxy of the local energy transfer rate (Marsch and Tu, 1997). This proxy was recently used to validate a multifractal model of the statistical properties of the turbulent fluctuations (Sorriso-Valvo et al., 2015). In this paper, it will be used as a tool to identify places where energy is being transferred towards the small scales.

By analogy with the definition of mean energy transfer rate, we define a “local” pseudo-energy transfer rate proxy (LET) by

$$ \epsilon_{\Delta t}^{\pm}(t) = \frac{|\Delta z_{\Delta t}^{\pm}(t)|^2 \Delta z_{\parallel, \Delta t}^{\mp}(t)}{\Delta t \langle \mathbf{v} \rangle}, \quad (2) $$

so that the local energy transfer rate at the scale $\Delta t$ is computed as $\epsilon_{\Delta t}(t) = (\epsilon_{\Delta t}^{+}(t) + \epsilon_{\Delta t}^{-}(t))/2$. At a given scale, each field increment can thus be associated with the local value of $\epsilon_{\Delta t}(t)$ (Marsch and Tu, 1997; Sorriso-Valvo et al., 2015). In terms of velocity and magnetic field, the LET is $\epsilon_{\Delta t}(t) \propto 2\Delta v_{\parallel}(\Delta v^2 + \Delta b^2) - 4\Delta b_{\parallel}(\Delta \mathbf{v} \cdot \Delta \mathbf{b})$, where the first term is associated with the energy advected by the velocity, and the second to the velocity–magnetic field correlations coupled to the longitudinal magnetic field. The statistical properties of this proxy will be explored here using *Helios* 2 measurements in the inner heliosphere. This article is organized as follows: in Section 2 we describe the data and the diagnostic variables used for the analysis; in Section 3 we study the scale-dependent statistical properties of the local energy transfer rate, and compare them with the PVI, a standard identification tool; finally, in Section 4 the correlation with the local proton temperature is studied.

## 2. Description of Data

This work presents a study on *Helios* 2 data, which have been thoroughly analyzed for about 40 years and still represent a milestone in the study of the inner heliosphere. More in particular, the data selected here are for 11 intervals taken during the first 4 months of 1976, at low solar activity, while *Helios* 2 orbit spanned between 1 AU, on day 17, to 0.29 AU on day 108. Each interval includes 2178 data points at 81-s cadence, covering about 2 days of measurements. All intervals were extracted during relatively stationary wind conditions, i.e. far from the inter-stream interaction regions, and include time series of magnetic field $\mathbf{B}(t)$, velocity $\mathbf{v}(t)$, proton number density $n_p(t)$ and proton temperature $T_p(t)$. Five samples refer to slow solar wind, with average bulk speed $V_{SW} \lesssim 450$ km s$^{-1}$, and six to fast solar wind, with $V_{SW} \gtrsim 550$ km s$^{-1}$. General parameters of the 11 intervals are given in Table 1.

In the next sections, we will use the *Helios* 2 data to characterize the statistical properties of LET.

## 3. Statistical Properties of the Local Energy Transfer Rate

The statistical properties of LET can provide quantitative information as regards the characteristics of the turbulent cascade occurring in the “fluid” range in solar wind plasma. For
<table><caption>**Table 1** List of the 11 49-hour intervals used for this work. For each interval: the day of the year 1976 when the sample begins, DOY; the heliocentric distance, R; the mean speed, $V_{SW}$ are indicated.</caption><thead><tr><th>DOY</th><th>R (AU)</th><th>V<sub>SW</sub> (km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>22</td><td>0.98</td><td>676</td></tr><tr><td>28</td><td>0.97</td><td>348</td></tr><tr><td>32</td><td>0.96</td><td>587</td></tr><tr><td>46</td><td>0.90</td><td>433</td></tr><tr><td>49</td><td>0.88</td><td>643</td></tr><tr><td>72</td><td>0.70</td><td>411</td></tr><tr><td>75</td><td>0.67</td><td>632</td></tr><tr><td>81</td><td>0.59</td><td>343</td></tr><tr><td>85</td><td>0.54</td><td>543</td></tr><tr><td>99</td><td>0.35</td><td>431</td></tr><tr><td>105</td><td>0.30</td><td>727</td></tr></tbody></table>

each of the samples indicated in Table 1, we have calculated the proxy of LET, $\epsilon_{\Delta t}(t)$, at different time scales $\Delta t$, using Equation 2 as described in Section 1. Figure 1 shows one example of $\epsilon_{\Delta t}(t)$, at two different scales (two upper panels). The scale-dependent nature of the signal is evident from the comparison of the two panels, and is the typical signature of intermittency, resulting in the strongly bursty field observed at small $\Delta t$. In order to compare LET with the standard PVI, we estimated the latter, including both velocity and magnetic structure contributions, as

$$PVI_{\Delta t}^2(t) = \frac{\Delta \mathbf{B}(t)^2}{\Delta \mathbf{B}_{\text{rms}}^2} + \frac{\Delta \mathbf{v}(t)^2}{\Delta \mathbf{v}_{\text{rms}}^2}, \quad (3)$$

where $\Delta \mathbf{B}_{\text{rms}}$ and $\Delta \mathbf{v}_{\text{rms}}$ indicate the standard deviation of the magnetic field and velocity increments at the scale $\Delta t$, computed over each interval. The main difference between the two proxies is that, while the PVI account for the amplitude of magnetic and velocity gradients (so being an estimate of electric current and vorticity structures), the LET carries information on the magnetic and kinetic energy coupled to the advecting velocity field, and on the cross-helicity coupled to the longitudinal magnetic field. Furthermore, contrary to PVI, the LET is signed, and might contain additional information as regards the local direction of the energy flux. Examples of $PVI_{\Delta t}^2(t)$ for the same interval and for the two same scales are given in Figure 1, in the third and fourth panels from the top. A similar scale-dependent burstiness is observed for the PVI, although qualitative difference is present. In particular, the LET appears burstier than PVI. For a direct visual comparison, in the bottom panel of Figure 1 we show the two proxies at the smallest scale $\Delta t = 81$ s, in a short time interval, as to compare the correspondence of energy bursts. As is evident, there are times of good correspondence between LET and PVI, but also times when only one of the two proxies has one or more peaks. Figure 2 shows scatter plots of the two variables for one fast (top panel) and one slow wind (bottom panel) interval, which demonstrates a good qualitative agreement between them. This is confirmed by the large associated Spearman correlation coefficient, $\rho_S \sim 0.9$. However, it should be noticed that the correlation is less evident when large values of $\epsilon_{\Delta t}$ and $PVI_{\Delta t}(t)$ are considered, i.e. where the energy flow is larger and at the most intense current structures. In particular, for the fast wind of DOY 85 (top panel) there is an evident presence of points with larger $\epsilon_{\Delta t}$ and small $PVI_{\Delta t}(t)$ (the isolated population lying above the correlated points), indicating times when an enhanced energy flux does not
<figure><img src="image_3.png" /><figcaption>**Figure 1** Examples of $\epsilon_{\Delta t}(t)$ and $\text{PVI}_{\Delta t}(t)$ at two different scales $\Delta t$, for the slow wind interval of DOY 28. A magnification of both proxies at the resolution scale $\Delta t = 81$ s is shown in the *bottom panel*, for a shorter time interval.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 2** Scatter plots of $|\epsilon_{\Delta t}(t)|$ versus $PVI_{\Delta t}(t)$ at the smallest scale $\Delta t = 81$ s, for the fast wind interval of DOY 85 (top) and for the slow wind interval of DOY 99 (bottom). The corresponding Spearman correlation coefficients $\rho_S$ are indicated.</figcaption></figure>

necessarily correspond to comparably strong current sheets. This effect is still present, but it is less relevant in slow wind intervals, where the smaller correlations between velocity and magnetic field reduce the difference between the two proxies.

In order to fully describe the statistical properties of the two proxies, we compute for each interval the probability distribution function (PDF) of both variables, at different scales. For the LET, we show the distribution $P(|\epsilon_{\Delta t}(t)|)$, having verified that the positive and negative parts of the variable have very similar statistics. Figure 3 illustrates the difference between the two proxies. It can be noted that at large values the functional form of the PDFs of the two proxies change. Although most of the data are in the core of the distribution (small values), the interesting large bursts represented in the right tails are particularly relevant for this work. Both proxies clearly show scale-dependent PDF, with the typical increase of the tails as the scale decreases, indicating the increasing presence of bursts of energy transfer (Sorriso-Valvo et al., 2015), typical of intermittency. The scaling behavior of the PDFs of $|\epsilon_{\Delta t}|$ and $PVI_{\Delta t}(t)$ can be described through standard models of the turbulent cascade. The distribution of $|\epsilon_{\Delta t}|$ appears to be a stretched exponential function (Sorriso-Valvo et al., 2015)

$$P(|\epsilon_{\Delta t}|) \sim \exp(-b|\epsilon_{\Delta t}|^c), \quad (4)$$

where b is a parameter related with the most probable value of the distribution, and $c \equiv c(\Delta t)$ describes the shape of the curve. In particular, for $c = 2$ a Gaussian PDF is
<figure><img src="image_5.png" /><figcaption>**Figure 3** Examples of PDFs of $|\epsilon_{\Delta t}(t)|$ (top panel) and $\text{PVI}_{\Delta t}(t)$ (bottom panel) at three different scales $\Delta t$, for the fast wind interval of DOY 22. The error bars are obtained by propagating the standard Poisson statistical uncertainty of each bin count.</figcaption></figure>

obtained, while $c=1$ corresponds to an exponential distribution. Values $0 < c < 1$ are associated with distribution whose tails can be more and more approximated by power laws as $c$ decreases. The stretched exponential distribution can be interpreted in the framework of the extreme deviations theory (EDT) applied to the fragmentation process occurring in the turbulent cascade, and more specifically when the statistics is controlled by a small number of extremely intense events (Frisch and Sornette, 1997). As turbulence is intermittent, this condition applies when the scale decreases, and EDT can be invoked to describe the statistics. Moreover, the PDF tail’s “flatness” can be quantitatively represented through the parameter $c$, so that smaller $c$ corresponds to higher-tailed distributions and then to a higher probability of occurrence of extreme intense events. Thus, this is the most relevant scale-dependent parameter to describe the intermittency of the system.

The fit of $P(|\epsilon_{\Delta t}|)$ at one time scale for the interval of DOY 22 is displayed in the top panel of Figure 4. Similar fitting quality is achieved for all scales and all datasets. For each data interval, it is thus possible to describe the scale dependence of $P(|\epsilon_{\Delta t}|)$ through $c(\Delta t)$. The top panel of Figure 5 shows one example of scaling of $c(\Delta t)$ for the same interval. In the same plot, the vertical bar indicates the turbulence correlation time estimated for this interval, $\tau_c \sim 1000$ s (Bruno et al., 2009). Note that time scales smaller than the correlation time lie in the inertial range, while larger time scales are often associated with the $1/f$ spectral region (Bruno and Carbone, 2013). A double power-law scaling $c(\Delta t) \sim \Delta t^\gamma$ is evident. Indeed, the break occurs near the correlation time, so that the scaling exponent $\gamma$ is different for the inertial range and for the $1/f$ range. Power laws are often associated with the pres-
<figure><img src="image_6.png" /><figcaption>**Figure 4** Examples of fitted PDF of $|\epsilon_{\Delta t}(t)|$ for the fast wind interval of DAY 22 (top panel) and $PVI_{\Delta t}(t)$ for the fast wind interval of DOY 32 (bottom panel) at the smallest time scale $\Delta t = 81$ s. The error bars are obtained by propagating the standard Poisson statistical uncertainty of each bin count.</figcaption></figure>

ence of correlations between bursts, indicating the presence of a non-stochastic process. The values found from power-law fits are $\gamma_{kol} = 0.076 \pm 0.01$ and $\gamma_{1/f} = 0.30 \pm 0.02$. According to EDT, $\gamma$ may be inversely proportional to the typical number of fragmentation steps $N_c$ occurring during the cascade. Therefore, larger $\gamma$ values are associated with a smaller number of steps in the cascade. For the example given here $N_c \simeq 12$ in the inertial range, while $N_c \simeq 3$ in the $1/f$ range. Although the EDT model might not be fully adapted to describe solar wind turbulence, it allows us to estimate the properties of intermittency based on the dissipation, at variance with standard approaches based on the field increments (PVI being an example of the latter). Moreover, the appearance of a power law in the scaling of the parameter $c$ in the inertial range is indicative that $|\epsilon_{\Delta t}(t)|$ is a suitable variable for the description of the turbulent cascade.

On the other hand, as shown in the bottom panel of Figure 4, the PVI distribution can be well described by a log-normal function:

$$P(PVI_{\Delta t}) = \frac{1}{(\sqrt{2\pi}\lambda PVI_{\Delta t})} e^{-\frac{\log(PVI_{\Delta t} - \langle PVI_{\Delta t} \rangle)^2}{2\lambda^2}}, \quad (5)$$

where $\langle PVI_{\Delta t} \rangle$ is the mean of PVI at the scale $\Delta t$ and $\lambda \equiv \lambda(\Delta t)$ is its scale-dependent variance, which determines the width of the distribution. The log-normal statistics of PVI can be naively understood in terms of the multiplicative process underlying the intermittent turbulent cascade (Frisch, 1995): at each position, for a given scale, the field increment is the result
**Figure 5** Examples of scaling of the fitting parameters *c*(*Δt*) for the fast wind interval of DOY 22 (*top panel*) and *λ*(*Δt*) for the fast wind interval of DOY 32 (*bottom panel*), obtained from the fit of the PDFs of |*ε*<sub>Δt</sub>(*t*)| and *PVI*<sub>Δt</sub>(*t*), respectively. The *vertical dashed line* indicates the turbulence correlation scale for those intervals (Bruno *et al.*, 2009). For the parameter *c*, two power-law fits are also superposed. The scaling exponent is *γ*<sub>kol</sub> = 0.076 ± 0.01 in the inertial range (corresponding to *N*<sub>c</sub> ≃ 12 steps in the intermittent cascade) and *γ*<sub>1/f</sub> = 0.30 ± 0.02 in the 1/f range (corresponding to only *N*<sub>c</sub> ≃ 3 steps in the cascade, if any exists). In both panels, *error bars* represent the uncertainty of the parameters obtained from the χ<sup>2</sup> minimization of the fitting procedure.

<figure><img src="image_7.png" /></figure>

of all the previous fragmentation steps, which can be expressed in terms of multiplicative random factors. After a large number of fragmentations, the logarithm of the field increment will be the sum of the logarithms of randomly distributed multiplicative factors, so that the central limit theorem applies and the final log(*PVI*) value will obey Gaussian statistics. The fitting procedure provides, again, a quantitative estimation of the non-randomness of the fields, in terms of the presence of high tails. The parameter *λ*(*Δt*) controls the width of the distribution, so that, for larger *λ*, the distribution is broader, and the tails include the increasingly larger bursts of *PVI*<sub>Δt</sub>(*t*) arising at small scales because of intermittency. As for the stretched exponential parameter *c*(*Δt*), the decrease of the parameter *λ*(*Δt*) with the scale, shown in the bottom panel of Figure 5 for DOY 32, suggests the non-self-similar nature of the fluctuations, although a clear power law is not identified in this case. This confirms that a common mechanism, the nonlinear energy cascade, underlies the generation of the energy
bursts in the turbulent field. Note that a similar scale dependence of the distributions was obtained analyzing the coarse-grained energy dissipation in MHD numerical simulations, and a proxy of dissipation, similar to PVI, in solar wind data (Zhdankin, Boldyrev, and Chen, 2016).

Finally, the estimation of a proxy for the local energy transfer rate permits the study of the multifractal properties of the turbulent cascade, which are related to intermittency (Frisch, 1995). In the framework of the multifractal cascade models, the energy dissipation rate should be distributed in space as a multifractal object (Paladin and Vulpiani, 1987). Thus, the inhomogeneous character of the time series $\epsilon_{\Delta t}(t)$ highlighted by our analysis can be interpreted in terms of the multifractal properties of the field. It is therefore interesting to use this approach to characterize the degree of intermittency of the solar wind data, as opposed to the standard approach based on the scaling properties of velocity and magnetic field increments (Burlaga, 1992; Macek, 2006; Macek, Wawrzaszek, and Carbone, 2011). In order to do so, we have evaluated the multifractal spectrum of the proxy ($\epsilon_{\Delta t}(t)$), estimated at the resolution scale $\Delta t = 81$ s, using a standard box-counting procedure (see for example the details given in Macek, Wawrzaszek, and Carbone, 2011; Sorriso-Valvo et al., 2017). The presence of singular structures (e.g. the bursts of local energy transfer rate) is revealed by the power-law scaling of the $q$th order partition functions ($\chi_q(\delta t) \propto \delta t^{\tau_q}$) of a suitably defined coarse-grained probability measure associated to the LET, for each scale $\delta t$. Note that $\delta t$ is the time scale over which the coarse graining is computed and is not related to the scale $\Delta t$ used for field increments computation. The set of scaling exponents ($\tau_q$) describes the inhomogeneity of the singularity strength, and thus the multifractal properties of the field (Frisch and Parisi, 1983). In particular, the exponents are expected to depend linearly on the order $q$ for mono-fractal objects, where only one singularity exponent is present. On the contrary, deviation from linearity indicates multifractality, i.e. a broader set of singularity exponents (Grassberger and Procaccia, 1983). Such deviation can be estimated using theoretical models, e.g. the p-model (Meneveau and Sreenivasan, 1987), which is the one adopted in this work. The p-model was originally developed for the description of the energy cascade in Navier-Stokes turbulence. It is a simple representation of the cascade in which the energy at one given position and scale is redistributed unevenly to two smaller scale (or “daughter”) structures. The fraction of energy transferred at each step to each daughter structure is given by a cascade of multipliers randomly extracted from a binomial distribution, i.e. $p$ or $1-p$, where $0 \le p \le 0.5$ is the parameter that determines the characteristics of the cascade. In this simplistic view, the scaling exponents $\tau_q$ are directly related to the value of $p$ through $\tau_q = -\log_2[p^q + (1-p)^q]$ (Meneveau and Sreenivasan, 1987). The parameter $p$ is thus a good quantitative measure of the deviation from self-similar (or fractal) scaling, i.e. of the degree of multifractality of the system. In particular, $p \approx 0.5$ is an indication of mono-fractal fields, while smaller values are associated with greater multifractality.

For each solar wind $|\epsilon_{\Delta t}(t)|$ sample, the probability measures and their partition functions $\chi_q(\delta t)$ have been computed, for $q \in [-3, 3]$ with step $dq = 0.05$. The values of $\tau_q$ have then been evaluated by a fit of the partition functions to power laws, in the range $\delta t \in [162, 2000]$ s. Then the exponents $\tau_q$ have been fitted to the p-model and the parameter $p$ has been estimated. In Figure 6 we show one example of the scaling exponents $\tau_q$ for DOY 22. A fit with the p-model is also indicated. The values of $p$ obtained for all intervals lie in the range $0.8 \le p \le 0.9$, and they are compatible with the usual strong intermittency parameters obtained using the fields increments (Horbury et al., 1997; Sorriso-Valvo et al., 2017). Thus, the strong intermittent character of the field is well captured by such a multifractal analysis, based on a simple multiplicative model.
<figure><img src="image_8.png" /><figcaption>**Figure 6** Scaling exponents $\tau_q$ (symbols) for the *Helios 2* interval of DOY 22 with the p-model fit (red line). For the interval in this example, $p = 0.87$.</figcaption></figure>

**4. Conditioned Analysis of Temperature**

After describing the general statistical properties of the LET, it is interesting to point out its relationship with solar wind temperature. To this aim, the same analysis carried out by Osman et al. (2012) has been performed on *Helios 2* data using both the PVI and the LET at the resolution scale $\Delta t = 81$ s, for the localization of the small-scale structures. In this case, all 11 intervals were analyzed together, in order to increase the statistical significance of the conditioning procedure. For this reason, for each sample the temperature time series was previously normalized to its mean value, allowing the comparison between different heliocentric distances and different wind type. Upon identification of five incremental thresholds of the PVI values, and ten (five positive, five negative) of the LET, the conditioned average proton temperature has been estimated as a function of the distance *D* from each PVI or LET structure, i.e. where their values overcome a given threshold, $\theta_{PVI}$ and $\theta_{\epsilon}$, respectively. This procedure provides the secular temperature profile around structures, $\langle T_p | \theta_{PVI} \rangle$ and $\langle T_p | \theta_{\epsilon} \rangle$. Figure 7 shows these temperature profiles as a function of the distance *D* from the structures, conditioned to the indicated thresholds $\theta_{\epsilon}$ (coded in different colors and symbols). Note that all the curves have been arbitrarily shifted vertically for clarity. While the low-threshold curve ($\theta_{\epsilon} = 1$) is approximately constant (no local temperature increase), there is a striking evidence of temperature increase localized near the energetic LET structures, which is more evident as the conditioning value of $\epsilon_{\Delta t}$ is increased. In particular, for ($\theta_{\epsilon} = 5$) the amplitude of the central peak indicates approximately 8% higher temperature where the energy transfer rate is higher. Such evidence results from the robust presence of hotter plasma near the structures, whereas the random fluctuations of the temperature are statistically canceled out far from these. There is also the appearance of a typical size of the higher-temperature site, which is approximately 160 s around the structure. This might depend on the scale under study, so that a deeper analysis is left to an investigation in progress, based on higher resolution data. For the data associated with negative energy flux (lower part of the panel), the temperature profile has an evident threshold-dependent, incremental decrease approaching
<figure><img src="image_9.png" /><figcaption>**Figure 7** Conditionally averaged normalized temperature $\langle T_p | \theta_{\epsilon} \rangle$ as a function of the distance from the structure center, for different positive and negative values of the threshold. All curves are arbitrarily vertically shifted for clarity.</figcaption></figure>

the structure, suggesting that whenever the energy flows from smaller to larger scale, the plasma is heated less than on average. This effect is less localized than the possible heating observed at positive $\epsilon_{\Delta t}(t)$, and may be similar to the effect observed by Osman *et al.* (2012) for low PVI data. Therefore, it is evident that there is a strong localization of higher plasma temperature near the sites of larger energy flux towards the small scales, possibly associated with local plasma heating, while the times with larger, negative energy flux are associated to colder plasma. For comparison, the same analysis was carried out using PVI, i.e. repeating the Osman *et al.* (2012) procedure, and it is depicted in Figure 8. In this case, indication of higher temperature at the PVI structures is strongly reduced with respect to LET.

The discrepancy between LET and PVI conditioning is related to the difference, already observed in Figure 2, between the two proxies for large values, which are the most relevant for this analysis. The better performance of the LET shows that it is a more sensitive proxy, able to highlight the possible turbulent heating properties even when using a limited size dataset.

# 5. Conclusions

Unveiling the connection between the processes occurring at fluid scales and at kinetic scales ranges is important to understand how weakly collisional space plasmas dissipate the energy
**Figure 8** Conditionally averaged normalized temperature $\langle T_p | \theta_{PVI} \rangle$ as a function of the distance from the structure center, for different positive values of the threshold. All curves are arbitrarily vertically shifted for clarity.

<figure><img src="image_10.png" /></figure>

cascading from large to small scales. In particular, the solar wind represents an example of a collisionless plasma with a clear indication of heating due to turbulent energy dissipation. At the same time, it has the important advantage of the availability of in situ measurements, at variance with other astrophysical plasmas. In this article we have presented one possible tool for the identification of the local transfer of turbulent energy across scales, the LET $\epsilon_{\Delta t}(t)$, based on the third-order moment scaling law for MHD turbulence (Politano and Pouquet, 1998). This proxy differs from the usual tools such as LIM and PVI, as it includes, besides the direct kinetic and magnetic energy contributions, also cross-terms representing for example the cross-helicity contribution. The statistical analysis of the proxy provides insight into the scaling properties of the turbulent cascade, consistent with the standard turbulence analysis of solar wind plasmas. We have used *Helios 2* data to describe the properties of LET as compared to the PVI. Good correlations are found between PVI and LET, confirming that both indicators are suitable for the description of the turbulent cascade of energy. However, for large energy flux important differences arise. Moreover, the signed variable $\epsilon_{\Delta t}(t)$ carries information as regards the possible direction of cross-scale energy flow, which is hidden in the positive-defined variables LIM and PVI. This could be useful for better understanding the coupling mechanisms occurring near the MHD break scale. The study of this particular aspect is being addressed in a different work. The LET has also been studied through multifractal analysis based on the dissipation, rather than on the field increments. Such an alternative analysis has confirmed the highly intermittent character of solar wind MHD turbulence. Finally, a convincing correspondence between times of enhanced energy transfer rate and local temperature increase has been clearly demonstrated, indicating that LET is a useful tool for the identification of regions of interest for the study of turbulent energy dissipation.

Because of its ability to track the link between the two ranges of scales, the LET could be useful for interpreting data from numerical simulations of the Vlasov-Maxwell system for the description of kinetic processes in collisionless plasmas. Similarly, it could have important implications for the analysis and interpretation of data from space missions providing high-resolution plasma measurements, such as MMS (Burch *et al.*, 2016) and the ESA candidate mission THOR (Vaivads *et al.*, 2016), but also for the forthcoming *Parker Solar Probe* (Fox *et al.*, 2016) and *Solar Orbiter* (Müller *et al.*, 2013).

**Acknowledgements** SP acknowledges support by the Agenzia Spaziale Italiana under the contract ASI-INAF 2015-039-R.O "Missione M4 di ESA: Partecipazione Italiana alla fase di assessment della missione
THOR". RM acknowledges financial support from the program PALSE (Programme Avenir Lyon Saint-Etienne) of the University of Lyon, in the frame of the program Investissements d'Avenir (No. ANR-11-IDEX-0007). The authors declare they have no conflict of interest arising from the above funding.

## References

Alexandrova, O., Chen, C.H.K., Sorriso-Valvo, L., Horbury, T.S., Bale, S.D.: 2013, *Space Sci. Rev.* **178**, 101. DOI.

Banerjee, S.: 2013, *Phys. Rev. E* **87**, 013019. DOI.

Banerjee, S., Hadid, L.Z., Sahraoui, F., Galtier, S.: 2016, *Astrophys. J. Lett.* **829**, L27. DOI.

Bavassano, B., Dobrowolny, M., Fanfoni, G., Mariani, F., Ness, N.F.: 1982, *Solar Phys.* **78**, 373. DOI.

Biskamp, D.: 1997, *Nonlinear Magnetohydrodynamics*, Cambridge University Press, Cambridge.

Bruno, R., Carbone, V.: 2013, *Living Rev. Solar Phys.* **10**, 2. DOI.

Bruno, R., Carbone, V., Veltri, P., Pietropaolo, E., Bavassano, B.: 2001, *Planet. Space Sci.* **49**, 1201. DOI.

Bruno, R., Carbone, V., Vörös, Z., D'Amicis, R., Bavassano, B., Cattaneo, M.B., Mura, A., Milillo, A., Orsini, S., Veltri, P., Sorriso-Valvo, L., Zhang, T., Biernat, H., Rucker, H., Baumjohann, W., Jankovičová, D., Kovács, P.: 2009, *Earth Moon Planets* **104**, 101. DOI.

Burch, J.L., Moore, T.E., Torbert, R.B., Giles, B.L.: 2016, *Space Sci. Rev.* **199**, 5. DOI.

Burlaga, L.F.: 1992, *J. Geophys. Res.* **97**, 4283. DOI.

Carbone, V., Marino, R., Sorriso-Valvo, L., Noullez, A., Bruno, R.: 2009, *Phys. Rev. Lett.* **103**, 061102. DOI.

Chasapis, A., Retinò, A., Sahraoui, F., Vaivads, A., Khotyaintsev, Yu.V., Sundkvist, D., Greco, A., Sorriso-Valvo, L., Canu, P.: 2015, *Astrophys. J. Lett.* **804**, L1. DOI.

Chen, C.H.K.: 2016, *J. Plasma Phys.* **82**, 53582060. DOI.

Coburn, J.T., Smith, C.W., Vasquez, B.J., Stawarz, J.E., Forman, M.A.: 2012, *Astrophys. J.* **754**, 93. DOI.

Fox, N.J., Velli, M.C., Bale, S.D., Decker, R., Driesman, A., Howard, R.A., Kasper, J.C., Kinnison, J., Kusterer, M., Lario, D., Lockwood, M.K., McComas, D.J., Raouafi, N.E., Szabo, A.: 2016, *Space Sci. Rev.* **204**, 7. DOI.

Freeman, J.W., Totten, T., Arya, S.: 1992, A determination of the polytropic index of the free streaming solar wind using improved temperature and density radial power-law indices. *Eos Trans. AGU* **73**, 238.

Frisch, U.: 1995, *Turbulence: The Legacy of A.N. Kolmogorov*, Cambridge University Press, Cambridge.

Frisch, U., Parisi, G.: 1983, In: Ghil, M., Benzi, R., Parisi, G. (eds.) *Turbulence and Predictability of Geophysical Flows and Climate Dynamics* **84**, North-Holland, Amsterdam.

Frisch, U., Sornette, D.: 1997, *J. Phys. 17*, 1155. DOI.

Galtier, S.: 2008, *Phys. Rev. E* **77**, 015302(R). DOI.

Gogoberidze, G., Perri, S., Carbone, V.: 2013, *Astrophys. J.* **769**, 111. DOI.

Goldstein, B.E., Neugebauer, M., Phillips, J.L., Bame, S., Gosling, J.T., McComas, D., Wang, Y.-M., Sheeley, N.R., Guess, S.T.: 1996, *Astron. Astrophys.* **316**, 296.

Grassberger, P., Procaccia, I.: 1983, *Phys. Rev. Lett.* **50**, 346. DOI.

Greco, A., Perri, S.: 2014, *Astrophys. J.* **784**, 163. DOI.

Greco, A., Matthaeus, W.H., Servidio, S., Dmitruk, P.: 2009, *Phys. Rev. E* **80**, 046401. DOI.

Hadid, L.Z., Sahraoui, F., Galtier, S.: 2017, *Astrophys. J.* **838**, 9. DOI.

Hellinger, P., Trávníček, P.M., Štverák, Š., Matteini, L., Velli, M.: 2013, *J. Geophys. Res.* **118**, 1351. DOI.

Horbury, T.S., Balogh, A., Forsyth, R.J., Smith, E.J.: 1997, *Adv. Space Res.* **19**, 847. DOI.

Kolmogorov, A.N.: 1941, *Dokl. Akad. Nauk SSSR* **30**, 301.

Leamon, R.J., Matthaeus, W.H., Smith, C.W., Zank, G.P., Mullan, D.J., Oughton, S.: 2000, *Astrophys. J.* **537**, 1054. DOI.

MacBride, B.T., Forman, M.A., Smith, C.W.: 2005, In: Fleck, B., Zurbuchen, T.H., Lacoste, H. (eds.) *Solar Wind 11/SOHO 16, Connecting Sun and Heliosphere*, ESA Special Publication **592**, 613.

MacBride, B.T., Smith, C.W., Forman, M.A.: 2008, *Astrophys. J.* **679**, 1644. DOI.

Macek, W.M.: 2006, *Space Sci. Rev.* **122**, 329. DOI.

Macek, W.M., Wawrzaszek, A., Carbone, V.: 2011, *Geophys. Res. Lett.* **38**, L19103. DOI.

Marino, R., Sorriso-Valvo, L., Carbone, V., Noullez, A., Bruno, B., Bavassano, B.: 2008, *Astrophys. J.* **677**, L71. DOI.

Marino, R., Sorriso-Valvo, L., Carbone, V., Veltri, P., Noullez, A., Bruno, R.: 2011, *Planet. Space Sci.* **59**, 592. DOI.

Marino, R., Sorriso-Valvo, L., D'Amicis, R., Carbone, V., Bruno, R., Veltri, P.: 2012, *Astrophys. J.* **750**, 41. DOI.

Marsch, E., Tu, C.: 1997, *Nonlinear Process. Geophys.* **4**, 101. DOI.
Matthaeus, W.H., Dasso, S., Weygand, J.M., Milano, L.J., Smith, C.W., Kivelson, M.J.: 2005, *Phys. Rev. Lett.* **95**, 231101. DOI.

Matthaeus, W.H., Oughton, S., Pontius, D.H. Jr., Zhou, Y.: 1994, *J. Geophys. Res.* **99**, 19267. DOI.

Meneveau, C., Sreenivasan, K.R.: 1987, *Phys. Rev. Lett.* **59**, 1424. DOI.

Müller, D., Marsden, R.G., St. Cyr, D.C., Gilbert, H.R.: 2013, *Solar Phys.* **285**, 27. DOI.

Osman, K.T., Wan, M., Matthaeus, W.H., Weygand, J.M., Dasso, S.: 2011, *Phys. Rev. Lett.* **107**, 165001. DOI.

Osman, K.T., Matthaeus, W.H., Wan, M., Rappazzo, A.F.: 2012, *Phys. Rev. Lett.* **108**, 261102. DOI.

Paladin, G., Vulpiani, A.: 1987, *Phys. Rep.* **156**, 147. DOI.

Perri, S., Servidio, S., Vaivads, A., Valentini, F.: 2017, *Astrophys. J., Suppl. Ser.* **231**, 4. DOI.

Politano, H., Pouquet, A.: 1998, *Geophys. Res. Lett.* **25**, 273. DOI.

Richardson, J.D., Paularena, K.I., Lazarus, A.J., Belcher, J.W.: 1995, *Geophys. Res. Lett.* **22**, 325. DOI.

Servidio, S., Valentini, F., Califano, F., Veltri, P.: 2012, *Phys. Rev. Lett.* **108**, 045001. DOI.

Schwenn, R.: 1983 In: *Solar Wind Five, NASA Conf. Publ., United States, CP-2280*, 489.

Smith, C.W., Stawarz, J.E., Vasquez, B.J., Forman, M.A., MacBride, B.T.: 2009, *Phys. Rev. Lett.* **103**, 201101. DOI.

Sorriso-Valvo, L., Yordanova, E., Carbone, V.: 2010, *Europhys. Lett.* **90**, 59001. DOI.

Sorriso-Valvo, L., Carbone, V., Bruno, R., Veltri, P.: 2006, *Europhys. Lett.* **75**, 832. DOI.

Sorriso-Valvo, L., Carbone, F., Leonardis, E., Chen, C.H.K., Šafránková, J., Němeček, Z.: 2017, *Adv. Space Res.* **59**, 1642. DOI.

Sorriso-Valvo, L., Carbone, V., Noullez, A., Politano, H., Pouquet, A., Veltri, P.: 2002, *Phys. Plasmas* **9**, 89. DOI.

Sorriso-Valvo, L., Marino, R., Carbone, V., Noullez, A., Lepreti, F., Veltri, P., Bruno, R., Bavassano, B., Pietropaolo, P.: 2007, *Phys. Rev. Lett.* **99**, 115001. DOI.

Sorriso-Valvo, L., Marino, R., Lijoi, L., Perri, S., Carbone, V.: 2015, *Astrophys. J.* **807**, 86. DOI.

Stawarz, J.E., Vasquez, B.J., Smith, C.W., Forman, M.A., Klewicki, J.: 2011, *Astrophys. J.* **736**, 44. DOI.

Taylor, G.I.: 1938, *Proc. Roy. Soc. London Ser. A, Math. Phys. Sci.* **164**(919), 476.

Tessein, J.A., Matthaeus, W.H., Wan, M., Osman, K.T., Ruffolo, D., Giacalone, J.: 2013, *Astrophys. J. Lett.* **776**, L8. DOI.

Tu, C., Marsch, E.: 1995, *Space Sci. Rev.* **73**, 1. DOI.

Vaivads, A., et al.: 2016, *J. Plasma Phys.* **82**, 905820501. DOI.

Veltri, P., Mangeney, A.: 1999, In: Habbai, S.R., Hesser, R., Hollweg, J.V., Isenberg, P.A. (eds.) *Solar Wind Nine, The American Institute of Physics, College Park*, 543.

Weygand, J.M., Matthaeus, W.H., Dasso, S., Kivelson, M.G., Walker, R.J.: 2007, *J. Geophys. Res.* **112**, A10201. DOI.

Zhdankin, V., Boldyrev, S., Chen, C.H.K.: 2016, *Mon. Not. Roy. Astron. Soc.* **457**, L69. DOI.