# Scaling laws of turbulence and heating of fast solar wind: the role of density fluctuations

## Abstract

Incompressible and isotropic magnetohydrodynamic turbulence in plasmas can be described by an exact relation for the energy flux through the scales. This Yaglom-like scaling law has been recently observed in the solar wind above the solar poles observed by the Ulysses spacecraft, where the turbulence is in an Alfvénic state. An analogous phenomenological scaling law, suitably modified to take into account compressible fluctuations, is observed more frequently in the same dataset. Large scale density fluctuations, despite their low amplitude, play thus a crucial role in the basic scaling properties of turbulence. The turbulent cascade rate in the compressive case can moreover supply the energy dissipation needed to account for the local heating of the non-adiabatic solar wind.

###### pacs:

96.50.Ci; 47.27.Gs; 96.50.Tf; 52.35.RaThe interplanetary space is permeated by the solar wind (1), a magnetized, supersonic flow of charged particles originating in the high solar atmosphere and blowing away from the sun. Low frequency fluctuations of solar wind variables are often described in the framework of fully developed hydromagnetic (MHD) turbulence (2); (3). The large range of scales involved, spanning from 1 AU ( km) down to a few kilometers, make the solar wind the largest “laboratory” where MHD turbulence can be investigated using measurements collected in situ by instruments onboard spacecraft (3). MHD turbulence is often investigated through the Elsässer variables , computed from the local plasma velocity and magnetic field , being the plasma mass density. In terms of such variables, MHD equations can be rewritten as , where is the total hydromagnetic pressure, and diss indicates dissipative terms involving the viscosity and the magnetic diffusivity. As in the Navier-Stokes equations for neutral fluids, the nonlinear terms cause the turbulent energy transfer between different scales, at high Reynolds numbers where dissipative terms can be neglected. However, in the MHD case, they couple the two Elsässer variables, so that the Alfvénic MHD fluctuations , propagating along the background magnetic field, are advected by fluctuations propagating in the opposite direction. The presence of strong correlations (or anti-correlations) between velocity and magnetic fluctuations, along with a nearly constant magnetic intensity and low amplitude density fluctuations, is usually referred to as Alfvénic state of turbulence, and implies that one of the two modes should be negligible, making the nonlinear term of MHD equations vanish for pure Alfvénic fluctuations. In that case, the turbulent energy transfer should also disappear (4). Alfvénic turbulence is observed almost ubiquitously in fast wind. This holds both in the ecliptic fast streams, and in the high latitude wind blowing directly from the sun coronal holes (5); (6); (3). As pointed out in (4), the observation of Alfvénic state turbulence in the solar wind represents therefore a paradoxical “contradiction in terms”.

MHD turbulence however satisfies an important analytical relation, which is the equivalent for magnetized fluids of the Kolmogorov or the Yaglom relations. Under suitable hypotheses, it has been shown (7); (8) that the pseudo-energy fluxes through the scale of the increments of the Elsässer fields follow a linear scaling relation

(1) |

Here, represents the component of the increment along the direction , and are the dissipation rates per unit mass of the pseudo-energies ( indicates space averages). The scaling law (1) has been recently observed experimentally in polar wind (8) and in the ecliptic plane (9); (10). The confirmation of the scaling law (1) is an important step towards the solution of the apparent paradox of the Alfvénic turbulent state, because it unambiguously shows that an MHD cascade is present, maybe with a weak transfer rate, despite the strong velocity-magnetic fields correlations.

Relation (1), which is of general validity within MHD turbulence, is not always realized in the solar wind observations (8). Indeed, it is possible that local characteristics of the solar wind plasma do not always satisfy the assumptions required for (1) to be valid, namely large-scale homogeneity, isotropy, and incompressibility. The role of anisotropy of solar wind turbulence, which is expected to be important, is not considered in this letter. Density fluctuations in solar wind have a low amplitude, so that nearly incompressible MHD framework is usually considered (11); (12). However, compressible fluctuations are observed, typically convected structures characterized by anti-correlation between kinetic pressure and magnetic pressure (13). Properties and interaction of the basic MHD modes in the compressive case have also been considered in the past (14); (15). In the present paper, we show that density fluctuations, despite their very low amplitude, play a central role in the turbulent energy transfer, and that a phenomenological scaling law obtained by taking into account density fluctuations is observed in a much larger proportion of fast solar wind. We also show that the turbulent dissipation can account for a large fraction of the local heating causing a slower than expected decrease of temperature with distance. A first attempt to include density fluctuations in the framework of fluid turbulence was due to Lighthill (16). He pointed out that in a compressible energy cascade, the mean energy transfer rate per unit volume should be constant in a statistical sense ( being the characteristic velocity fluctuations at the scale ), obtaining . Fluctuations of a density-weighted velocity field should thus follow the usual Kolmogorov scaling . The same phenomenological conjecture can be introduced in MHD turbulence by considering the pseudo-energy dissipation rates per unit volume , and introducing density-weighted Elsässer fields, defined as . The equivalent of the Yaglom-type relation

(2) |

should then hold for the density-weighted increments . Note that we have defined the flux so that it reduces to in the case of constant density, allowing for comparisons between the compressible scaling (2) and the purely incompressible one (1). Despite its simple phenomenological derivation, the introduction of the density fluctuations in the Yaglom-type scaling (2) seems to describe correctly the turbulent cascade for compressible fluid (or magnetofluid) turbulence. The law for the velocity field has been observed in recent numerical simulations (17); (18).

We will now study the cascade properties of compressive MHD turbulence from solar wind data collected by spacecraft Ulysses. In order to avoid as far as possible variations due to solar activity, or other ecliptic disturbances such as slow wind sources, coronal mass ejection, current sheets, we concentrate our analysis on pure Alfvénic state turbulence observed in high speed polar wind. We use here measurements from Ulysses spacecraft in the first six months of 1996. Such period was characterized by low solar activity, so that solar origin disturbances werebv almost absent. Moreover, the spacecraft orbit was at high and slowly decreasing heliolatitude, from about to , and presented small variations of the heliocentric distance , from 3 AU to 4 AU. Since the mean wind speed in the spacecraft frame is much larger than the typical velocity fluctuations, and is nearly aligned with the radial direction , space scales can be viewed as time scales , related through the Taylor hypothesis by . We then used 8 minutes averaged time series of both Elsässer variables and density (obtained as the sum of proton density and 4 times He density), to compute the density-weighted time series . From this time series we calculate the increments for different time lags , and the third-order mixed structure functions by time averaging over windows of fixed duration . The same procedure has been also been used to calculate the quantities using the time series of the Elsässer fields . In order to eliminate instationarities, heliolatitude and heliocentric distance changes, and to explore the wind properties locally, averages are computed over a moving window of about 11 days, consisting of 2048 data points. Accuracy of third order moments estimate (19) was tested with such sample size (20). We found that the third-order structure functions computed from the Ulysses data show a linear scaling

(3) |

during a considerable fraction of the period under study. In particular, we observed linear scaling of in about half of the signal, while displays scaling on about a quarter of the sample. As comparison, the corresponding incompressive scaling law for was only observed in a third of the whole period, considerably smaller than the compressible case (8). The portions of wind where the scaling is present are distributed in the whole period, and their extensions span from 6 hours up to 10 days. The linear scaling law generally extends on about 2 decades, from a few minutes to one day or more. For the compressible scaling, the two fluxes coexist in a large number of cases. This does not hold for the incompressive scaling, where in general the scaling periods for the two fluxes are disjoint.

Figure 1 shows one example of both mixed third-order structure functions and computed in the same 11 days windows where the scaling was observed. Figure 2 shows two more examples of scaling, observed both for and , in two different time windows. The scaling extends over 2 decades, while behaves linearly on the whole range of scales considered here (3 decades). In the last example, the scaling is not present for the incompressible fluxes . This example shows that the inclusion of compressible effect through the density-weigthed fluctuations improves the scaling (2) and modify the energy cascade. The scaling relation (2) also allows a direct estimate of the pseudo-energy transfer rates in the compressible case. A fit of the linear law (3) provides the local values of the amount of pseudo-energy transferred from large to small scales by the turbulent MHD cascade. This was already measured in the incompressive case (8); (21), so that it is possible to compare the transfer rates in the two cascades. The mean values, computed over the 46 observed scaling cases at different radial distances from the sun, ( their dispersion, in [J kg sec]) for the compressible cascade are (29 cases) and (17 cases). Both values are considerably larger than the corresponding values for the incompressive case (, 24 cases, and , 11 cases (21)). This result shows again that the cascade in the solar wind is strongly enhanced by density fluctuations, despite their small amplitude. Note that the new variables are built by coupling the Elsässer fields with the density, before computing the scale dependent increments. Moreover, the third order moments are very sensitive to intense field fluctuations (intermittency), that could arise when density fluctuations are correlated with velocity and magnetic field. Similar results, but with considerably smaller effect, were found in numerical simulations of compressive MHD (22). We should point out that experimental values of energy transfer rate in the incompressive case had been also estimated with different tecniques from different datasets (9); (10). Those values are not in agreement with the present (incompressive case) results. However, the different nature of wind (ecliptic vs polar, fast vs slow, at different radial distances from the sun) makes such comparison only indicative.

An interesting open question is the problem of the solar wind heating. The first models of solar wind assumed an adiabatic cooling due to spherical expansion of plasma blowing out of the sun. This would result in a radial decrease of the proton temperature with . On the contrary, spacecraft measurements (27) have shown that the temperature decay is slower than the adiabatic prescription, with . This implies that some local heating mechanism is present. One standing hypothesis is that the heating could be provided by energy dissipation occurring at the small scales of a turbulent cascade (23); (9); (24). Using equation (1), the rates at which the incompressible turbulent pseudo-energy is transported down the scales, and eventually dissipated at small scale, can be measured directly from data. This has recently be used to investigate whether or not a turbulent cascade can heat the solar wind. Results were however not conclusive. In fact, the measured incompressible dissipation rate of pseudo-energies can only account for up to 50 % of the solar wind heating (10); (21).

Figure 3 shows the radial profiles of the pseudo-energy transfer rates for both the compressive and incompressive cascades. In the same figure, we show the profiles of the heating rates needed to obtain the observed temperatures, as estimated from heating models (9); (24); (21) and from the measured temperatures (the two different values refer to the different estimates of the temperature obtained from Ulysses instruments). It is evident that, while the incompressive cascade cannot provide all the energy needed to heat the wind, the density fluctuations coupled with magnetohydrodynamic turbulence can supply the amount of energy required. This evidence shows the importance of the density fluctuations in polar, fast solar wind turbulence, confirming that it should be considered as an example of compressive fully developed MHD turbulence. Note that, since in a few samples we measured both and in the same period, the values of the energy and cross-helicity transfer rates can be disentangled. From the values obtained, it is clear that the cross-helicity contribution, indicating the importance of the Alfvénic state of turbulence, can vary from a negligible fraction (less than 1 %) to a considerable 25 % of the energy contribution. Since its amplitude does not appear to be correlated with the observation of the cascade, Alfvénicity seems not to play a crucial role in the cascade at the observed scales. This would be in agreement with previous analysis of solar wind turbulence anisotropy, where the Alfvénic contribution to the field fluctuations is small (25); (26).

To summarize, we used the density-weighted Elsässer fields to show for the first time that a phenomenological compressive Yaglom-like relation is verified to a large extent within the solar wind turbulence. This implies that low amplitude density fluctuations play a crucial role for scaling laws of solar wind turbulence (22). This observation also confirm the recent results for the Kolmogorov -law from numerical simulations of compressible turbulence (17), while no experimental evidences from real fluids had been found so far. This could be attributed to the incompressible nature of flows in ordinary fluids accessible to laboratory experiments. Here in fact, we present the first experimental observation of relation (2) in real systems. Using solar wind data, we have access to a sample of weakly compressible MHD turbulence in nature. Scaling law is found to be quite common and extends on a large range of scales, indicating not only that a nonlinear MHD cascade for pseudo-energies is active in the solar wind turbulence, but also that compressible effects are an important ingredient of the cascade. We point out that the observed departures from the scaling law could be due to presence of inhomogeneity and anisotropy in the solar wind (28). The compressive corrections to the cascade also cause the transfer of a considerably larger amount of energy toward the small scales, where it can be dissipated to heat the plasma locally. The role of anisotropy in the solar wind turbulent cascade still remains an open question.

### References

- A.J. Hundhausen, Coronal Expansion and Solar Wind, Springer, New York (1972)
- C.-Y. Tu and E. Marsch, Space Sci. Rev. 73, 1 (1995)
- R. Bruno and V. Carbone, Living Rev. Solar Phys. 2, 4 (2005)
- M. Dobrowolny, A. Mangeney and P. Veltri, Phys. Rev. Lett. 45, 144 (1980)
- J.M. Belcher and L. Davis Jr, J. Geophys. Res. 76, 3534 (1971)
- R. Bruno, B. Bavassano and U. Villante, J. Geophys. Res. 90, 4373 (1985)
- H. Politano, and A. Pouquet, J. Geophys. Lett. 25, 273 (1998)
- L. Sorriso-Valvo, et al., Phys. Rev. Lett. 99, 115001 (2007)
- B.J. Vasquez, et al. J. Geophys. Res. 112, A07101 (2007)
- B.T. MacBride, C.W. Smith, and M.A. Forman, Astrophys J. 679, 1644 (2008)
- D. Montgomery, M.R. Brown, and W.H. Matthaeus, J. Geophys. Res. 92, 282 (1987); W.H. Matthaeus, and M.R. Brown, Phys. Fluids 31, 3634 (1988); G.P. Zank, and W.H. Matthaeus, Phys. Fluids 3, 69 (1991); G.P. Zank, and W.H. Matthaeus, Phys. Fluids 5, 257 (1993)
- W.H. Matthaeus, et al., J. Geophys. Res. 96, 5421 (1991); B. Bavassano, and R.Bruno, J. Geophys. Res. 100, 9475 (1995)
- C.-Y. Tu, and E. Marsch, J. Geophys. Res. 99, 21481 (1994)
- P. Goldreich, and S. Sridhar, Astrophys. J. 438, 763 (1995)
- J. Cho, and A. Lazarian, Phys. Rev. Lett. 88, 245001 (2002)
- M.J. Lighthill, in IAU Symp. 2, Gas Dynamics of Cosmic Clouds (Amsterdam: North Holland), 121 (1955)
- A.G. Kritsuk, et al., Astrophys. J. 665, 416 (2007)
- G. Kowal, and A. Lazarian, Astrophys. J. 666, L69 (2007)
- J. J. Podesta, et al., Nonlin. Process. in Geophys. 16, 99 (2009)
- T. Dudok de Wit, Phys. Rev. E 70, 055302R (2004)
- R. Marino, et al., Astrophys J. 677, L71 (2008)
- M.-M. Mac Low, Astrophys. J. 524, 169 (1999)
- C. W. Smith, et al., J. Geophys. Res. 106, 8253 (2001)
- M.K. Verma, D.A. Roberts, and M.L. Goldstein, J. Geophys. Res. 100, 19839 (1995)
- J.W. Bieber, W. Wanner, and W.H. Matthaeus, J. Geophys. Res. 101, 2511 (1996)
- T.S. Horbury, M.A. Forman, and S. Oughton, Plasma phys. Control. Fusion 47, B703 (2005)
- B.E. Goldstein, et al., Astron. Astrophys. 316, 296 (1996)
- J. J. Podesta, M. A. Forman and C. W. Smith, Phys. Plasmas 14, 092305 (2007)