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

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


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.

96.50.Ci; 47.27.Gs; 96.50.Tf; 52.35.Ra

The 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


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


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


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: One example of the mixed third order compressible pseudo-energy flux  as computed from the Ulysses data during days 23 to 32 of 1996. The incompressible flux  in the same time window and a linear fit are also indicated. In this case, both compressible and incompressible fluxes obey a Yaglom-like law.
Figure 2: Top panel: an example of the third order compressible pseudo-energy flux  during days 1 to 10 of 1996. Bottom panel: for days 66 to 75 of the same year. In both panels, the corresponding incompressible fluxes  (no scaling present) and a linear fit are displayed.

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: Radial profile of the pseudo-energy transfer rates obtained from the turbulent cascade rate through the Yaglom relation, for both the compressive and incompressive case. The solid lines represent the radial profiles of the heating rate required to obtain the observed temperature profile.

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.


  1. A.J. Hundhausen, Coronal Expansion and Solar Wind, Springer, New York (1972)
  2. C.-Y. Tu and E. Marsch, Space Sci. Rev. 73, 1 (1995)
  3. R. Bruno and V. Carbone, Living Rev. Solar Phys. 2, 4 (2005)
  4. M. Dobrowolny, A. Mangeney and P. Veltri, Phys. Rev. Lett. 45, 144 (1980)
  5. J.M. Belcher and L. Davis Jr, J. Geophys. Res. 76, 3534 (1971)
  6. R. Bruno, B. Bavassano and U. Villante, J. Geophys. Res. 90, 4373 (1985)
  7. H. Politano, and A. Pouquet, J. Geophys. Lett. 25, 273 (1998)
  8. L. Sorriso-Valvo, et al., Phys. Rev. Lett. 99, 115001 (2007)
  9. B.J. Vasquez, et al. J. Geophys. Res. 112, A07101 (2007)
  10. B.T. MacBride, C.W. Smith, and M.A. Forman, Astrophys J. 679, 1644 (2008)
  11. 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)
  12. W.H. Matthaeus, et al., J. Geophys. Res. 96, 5421 (1991); B. Bavassano, and R.Bruno, J. Geophys. Res. 100, 9475 (1995)
  13. C.-Y. Tu, and E. Marsch, J. Geophys. Res. 99, 21481 (1994)
  14. P. Goldreich, and S. Sridhar, Astrophys. J. 438, 763 (1995)
  15. J. Cho, and A. Lazarian, Phys. Rev. Lett. 88, 245001 (2002)
  16. M.J. Lighthill, in IAU Symp. 2, Gas Dynamics of Cosmic Clouds (Amsterdam: North Holland), 121 (1955)
  17. A.G. Kritsuk, et al., Astrophys. J. 665, 416 (2007)
  18. G. Kowal, and A. Lazarian, Astrophys. J. 666, L69 (2007)
  19. J. J. Podesta, et al., Nonlin. Process. in Geophys. 16, 99 (2009)
  20. T. Dudok de Wit, Phys. Rev. E 70, 055302R (2004)
  21. R. Marino, et al., Astrophys J. 677, L71 (2008)
  22. M.-M. Mac Low, Astrophys. J. 524, 169 (1999)
  23. C. W. Smith, et al., J. Geophys. Res. 106, 8253 (2001)
  24. M.K. Verma, D.A. Roberts, and M.L. Goldstein, J. Geophys. Res. 100, 19839 (1995)
  25. J.W. Bieber, W. Wanner, and W.H. Matthaeus, J. Geophys. Res. 101, 2511 (1996)
  26. T.S. Horbury, M.A. Forman, and S. Oughton, Plasma phys. Control. Fusion 47, B703 (2005)
  27. B.E. Goldstein, et al., Astron. Astrophys. 316, 296 (1996)
  28. J. J. Podesta, M. A. Forman and C. W. Smith, Phys. Plasmas 14, 092305 (2007)
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description