We analyse velocity fluctuations in the solar wind at magneto-fluid scales in two datasets, extracted from Wind data in the period 2005-2015, that are characterised by strong or weak expansion. Expansion affects measurements of anisotropy because it breaks axisymmetry around the mean magnetic field. Indeed, the small-scale three-dimensional local anisotropy of magnetic fluctuations ({\delta}B) as measured by structure functions (SF_B) is consistent with tube-like structures for strong expansion. When passing to weak expansion, structures become ribbon-like because of the flattening of SFB along one of the two perpendicular directions. The power-law index that is consistent with a spectral slope -5/3 for strong expansion now becomes closer to -3/2. This index is also characteristic of velocity fluctuations in the solar wind. We study velocity fluctuations ({\delta}V) to understand if the anisotropy of their structure functions (SF_V ) also changes with the strength of expansion and if the difference with the magnetic spectral index is washed out once anisotropy is accounted for. We find that SF_V is generally flatter than SF_B. When expansion passes from strong to weak, a further flattening of the perpendicular SF_V occurs and the small-scale anisotropy switches from tube-like to ribbon-like structures. These two types of anisotropy, common to SF_V and SF_B, are associated to distinct large-scale variance anisotropies of {\delta}B in the strong- and weak-expansion datasets. We conclude that SF_V shows anisotropic three-dimensional scaling similar to SF_B, with however systematic flatter scalings, reflecting the difference between global spectral slopes.

We study the anisotropy with respect to the local magnetic field of turbulent magnetic fluctuations at magnetofluid scales in the solar wind. Previous measurements in the fast solar wind obtained axisymmetric anisotropy, despite that the analysis method allows nonaxisymmetric structures. These results are probably contaminated by the wind expansion that introduces another symmetry axis, namely, the radial direction, as indicated by recent numerical simulations. These simulations also show that while the expansion is strong, the principal fluctuations are in the plane perpendicular to the radial direction. Using this property, we separate 11 yr of Wind spacecraft data into two subsets characterized by strong and weak expansion and determine the corresponding turbulence anisotropy. Under strong expansion, the small-scale anisotropy is consistent with the Goldreich & Sridhar critical balance. As in previous works, when the radial symmetry axis is not eliminated, the turbulent structures are field-aligned tubes. Under weak expansion, we find 3D anisotropy predicted by the Boldyrev model, that is, turbulent structures are ribbons and not tubes. However, the very basis of the Boldyrev phenomenology, namely, a cross-helicity increasing at small scales, is not observed in the solar wind: the origin of the ribbon formation is unknown.

The spectral anisotropy of turbulent structures has been measured in the solar wind since 1990, relying on the assumption of axisymmetry about the mean magnetic field, B0. However, several works indicate that this hypothesis might be partially wrong, thus raising two questions: (i) is it correct to interpret measurements at 1 AU (the so-called Maltese cross) in term of a sum of slab and 2D turbulence? (ii) what information is really contained in the Maltese cross? We solve direct numerical simulations of the MHD equations including the transverse stretching exerted by the solar wind flow and study the genuine 3D anisotropy of turbulence as well as that one resulting from the assumption of axisymmetry about B0. We show that the evolution of the turbulent spectrum from 0.2 to 1 AU depends strongly on its initial anisotropy. An axisymmetric spectrum with respect to B0 keeps its axisymmetry, i.e., resists stretching perpendicular to radial, while an isotropic spectrum becomes essentially axisymmetric with respect to the radial direction. We conclude that close to the Sun, slow-wind turbulence has a spectrum that is axisymmetric around B0 and the measured 2D component at 1 AU describes the real shape of turbulent structures. On the contrary, fast-wind turbulence has a more isotropic spectrum at the source and becomes radially symmetric at 1 AU. Such structure is hidden by the symmetrization applied to the data that instead returns a slab geometry.

3D Magnetohydrodynamic (MHD) turbulent flows with initially magnetic and kinetic energies at equipartition spontaneously develop a magnetic excess (or residual energy), as well in numerical simulations and in the solar wind. Closure equations obtained in 1983 describe the residual spectrum as being produced by a dynamo source proportional to the total energy spectrum, balanced by a linear Alfv\'en damping term. A good agreement was found in 2005 with incompressible simulations; however, recent solar wind measurements disagree with these results. The previous dynamo-Alfv\'en theory is generalized to a family of models, leading to simple relations between residual and total energy spectra. We want to assess these models in detail against MHD simulations and solar wind data. The family of models is tested against compressible decaying MHD simulations with low Mach number, low cross-helicity, zero mean magnetic field, without or with expansion terms (EBM or expanding box model). A single dynamo-Alfv\'en model is found to describe correctly both solar wind scalings and compressible simulations without or with expansion. It is equivalent to the 1983-2005 closure equation but with critical balance of nonlinear turnover and linear Alfv\'en times, while the dynamo source term remains unchanged. The discrepancy with previous incompressible simulations is elucidated. The model predicts a linear relation between the spectral slopes of total and residual energies $m_R = -1/2 + 3/2 m_T$. Examining the solar wind data as in \cite{2013ApJ...770..125C}, our relation is found to be valid whatever the cross-helicity, even better so at high cross-helicity, with the total energy slope varying from $1.7$ to $1.55$.

We study the anisotropy of II-order structure functions defined in a frame attached to the local mean field in three-dimensional (3D) direct numerical simulations of magnetohydrodynamic turbulence, including or not the solar wind expansion. We simulate spacecraft flybys through the numerical domain by taking increments along the radial (wind) direction that forms an angle of $45^o$ with the ambient magnetic field. We find that only when expansion is taken into account, do the synthetic observations match the 3D anisotropy observed in the solar wind, including the change of anisotropy with scales. Our simulations also show that the anisotropy changes dramatically when considering increments oblique to the radial directions. Both results can be understood by noting that expansion reduces the radial component of the magnetic field at all scales, thus confining fluctuations in the plane perpendicular to the radial. Expansion is thus shown to affect not only the (global) spectral anisotropy, but also the local anisotropy of second-order structure functions by influencing the distribution of the local mean field, which enters this higher-order statistics.

The measure of the third-order structure function, Y, is employed in the solar wind to compute the cascade rate of turbulence. In the absence of a mean field B0=0, Y is expected to be isotropic (radial) and independent of the direction of increments, so its measure yields directly the cascade rate. For turbulence with mean field, as in the solar wind, Y is expected to become more two dimensional (2D), that is, to have larger perpendicular components, loosing the above simple symmetry. To get the cascade rate one should compute the flux of Y, which is not feasible with single-spacecraft data, thus measurements rely upon assumptions about the unknown symmetry. We use direct numerical simulations (DNS) of magneto-hydrodynamic (MHD) turbulence to characterize the anisotropy of Y. We find that for strong guide field B0=5 the degree of two-dimensionalization depends on the relative importance of shear and pseudo polarizations (the two components of an Alfv\'en mode in incompressible MHD). The anisotropy also shows up in the inertial range. The more Y is 2D, the more the inertial range extent differs along parallel and perpendicular directions. We finally test the two methods employed in observations and find that the so-obtained cascade rate may depend on the angle between B0 and the direction of increments. Both methods yield a vanishing cascade rate along the parallel direction, contrary to observations, suggesting a weaker anisotropy of solar wind turbulence compared to our DNS. This could be due to a weaker mean field and/or to solar wind expansion.

Forced, weak MHD turbulence with guide field is shown to adopt different regimes, depending on the magnetic excess of the large forced scales. When the magnetic excess is large enough, the classical perpendicular cascade with $5/3$ scaling is obtained, while when equipartition is imposed, an isotropic $3/2$ scaling appears in all directions with respect to the mean field (\cite{2010PhRvE..82b6406G} or GM10). We show here that the $3/2$ scaling of the GM10 regime is not ruled by a small-scale cross-helicity cascade, and propose that it is a 3D extension of a perpendicular weak Iroshnikov-Kraichnan (IK) cascade. We analyze in detail the structure functions in real space and show that they closely follow the critical balance relation both in the local frame and the global frame: we show that there is no contradiction between this and the isotropic $3/2$ scaling of the spectra. We propose a scenario explaining the spectral structure of the GM10 regime, that starts with a perpendicular weak IK cascade and extends to 3D by using quasi-resonant couplings. The quasi-resonance condition happens to reduce the energy flux in the same way as is done in the weak perpendicular cascade, so leading to a $3/2$ scaling in all directions. We discuss the possible applications of these findings to solar wind turbulence.

We study the evolution of turbulence in the solar wind by solving numerically the full 3D magneto-hydrodynamic (MHD) equations embedded in a radial mean wind. The corresponding equations (expanding box model or EBM) have been considered earlier but never integrated in 3D simulations. Here, we follow the development of turbulence from 0.2 AU up to about 1.5 AU. Starting with isotropic spectra scaling as $k^{-1}$, we observe a steepening toward a $k^{-5/3}$ scaling in the middle of the wavenumber range and formation of spectral anisotropies. The advection of a plasma volume by the expanding solar wind causes a non-trivial stretching of the volume in directions transverse to radial and the selective decay of the components of velocity and magnetic fluctuations. These two effects combine to yield the following results. (i) Spectral anisotropy: gyrotropy is broken, and the radial wavevectors have most of the power. (ii) Coherent structures: radial streams emerge that resemble the observed microjets. (iii) Energy spectra per component: they show an ordering in good agreement with the one observed in the solar wind at 1 AU. The latter point includes a global dominance of the magnetic energy over kinetic energy in the inertial and $f^{-1}$ range and a dominance of the perpendicular-to-the-radial components over the radial components in the inertial range. We conclude that many of the above properties are the result of evolution during transport in the heliosphere, and not just the remnant of the initial turbulence close to the Sun.

The magnetopause boundary implies two kinds of variations: a density/ temperature gradient and a magnetic field rotation. These two kinds are always observed in a close vicinity of each other, if not inseparably mixed. We present a case study from the Cluster data where the two are clearly separated and investigate the natures of both layers. We evidence that the first one is a slow shock while the second is a rotational discontinuity. The interaction between these two kinds of discontinuities is then studied with the help of 1,5-D magnetohydrodynamics simulations. The comparison with the data is quite positive and leads to think that most of the generic properties of the magnetopause may be interpreted in this sense.

The dynamical response of the solar corona to surface and sub-surface perturbations depends on the chromospheric stratification, and specifically on how efficiently these layers reflect or transmit incoming Alfv\'en waves. While it would be desirable to include the chromospheric layers in the numerical simulations used to study such phenomena, that is most often not feasible. We defined and tested a simple approximation allowing the study of coronal phenomena while taking into account a parametrised chromospheric reflectivity. We addressed the problems of the transmission of the surface rotation to the corona and that of the generation of polar plumes by Alfv\'en waves (Pinto et al., 2010, 2011). We found that a high (yet partial) effective chromospheric reflectivity is required to properly describe the angular momentum balance in the corona and the way the surface differential rotation is transmitted upwards. Alfv\'en wave-driven polar plumes maintain their properties for a wide range of values for the reflectivity, but they become bursty (and eventually disrupt) when the limit of total reflection is attained.

The description of the turbulent spectrum of magnetic fluctuations in the solar wind in the kinetic range of scales is not yet completely established. Here, we perform a statistical study of 100 spectra measured by the STAFF instrument on the Cluster mission, which allows to resolve turbulent fluctuations from ion scales down to a fraction of electron scales, i.e. from $\sim 10^2$ km to $\sim 300$ m. We show that for $k_{\perp}\rho_e \in[0.03,3]$ (that corresponds approximately to the frequency in the spacecraft frame $f\in [3,300]$ Hz), all the observed spectra can be described by a general law $E(k_\perp)\propto k_\perp^{-8/3}\exp{(-k_\perp \rho_e)}$, where $k_{\perp}$ is the wave-vector component normal to the background magnetic field and $\rho_e$ the electron Larmor radius. This exponential tail found in the solar wind seems compatible with the Landau damping of magnetic fluctuations onto electrons.

We propose a phenomenological theory of incompressible magnetohydrodynamic turbulence in the presence of a strong large-scale magnetic field, which establishes a link between the known anisotropic models of strong and weak MHD turbulence We argue that the Iroshnikov-Kraichnan isotropic cascade develops naturally within the plane perpendicular to the mean field, while oblique-parallel cascades with weaker amplitudes can develop, triggered by the perpendicular cascade, with a reduced flux resulting from a quasi-resonance condition. The resulting energy spectrum $E(k_\parallel,k_\bot)$ has the same slope in all directions. The ratio between the extents of the inertial range in the parallel and perpendicular directions is equal to $b_{rms}/B_0$. These properties match those found in recent 3D MHD simulations with isotropic forcing reported in [R. Grappin and W.-C. M\"uller, Phys. Rev. E \textbf{82}, 26406 (2010)].

The transition from weak to strong turbulence when passing from large to small scales in magnetohydrodynamic (MHD) turbulence with guide field is a cornerstone of anisotropic turbulence theory. We present the first check of this transition, using the Shell-RMHD which combines a shell model of perpendicular nonlinear coupling and linear propagation along the guide field. This model allows us to reach Reynolds numbers around $10^6$. We obtain surprisingly good agreement with the theoretical predictions, with a reduced perpendicular energy spectrum scaling as $k^{-2}$ at large scales and as $k^{-5/3}$ at small scales, where critical balance between nonlinear and propagation time is reached. However, even in the strong regime, a high level of excitation is found in the weak coupling region of Fourier space, which is due to the rich frequency spectrum of large eddies. A corollary is that the reduced parallel spectral slope is not a definite test of the spectral anisotropy, contrary to standard belief.

We present a mechanism for the formation of the low frequency 1/f magnetic spectrum based on numerical solutions of a shell reduced-MHD model of the turbulent dynamics inside the sub-Alfv\'enic solar wind. We assign reasonably realistic profiles to the wind speed and the density along the radial direction, and a radial magnetic field. Alfv\'en waves of short periodicity (600 s) are injected at the base of the chromosphere, penetrate into the corona and are partially reflected, thus triggering a turbulent cascade. The cascade is strong for the reflected wave while it is weak for the outward propagating waves. Reflection at the transition region recycles the strong turbulent spectrum into the outward weak spectrum, which is advected beyond the Alfv\'enic critical point without substantial evolution. There, the magnetic field has a perpendicular power-law spectrum with slope close to the Kolmogorov -5/3. The parallel spectrum is inherited from the frequency spectrum of large (perpendicular) eddies. The shape is a double power-law with slopes of -1 and -2 at low and high frequencies respectively, the position of the break depending on the injected spectrum. We suggest that the double power-law spectrum measured by Helios at 0.3 AU, where the average magnetic field is not aligned with the radial (contrary to our assumptions) results from the combination of such different spectral slopes. At low frequency the parallel spectrum dominates with its characteristic 1/f shape, while at higher frequencies its steep spectral slope (-2) is masked by the more energetic perpendicular spectrum (slope -5/3).

The turbulent spectrum of magnetic fluctuations in the solar wind displays a spectral break at ion characteristic scales. At electron scales the spectral shape is not yet completely established. Here, we perform a statistical study of 102 spectra at plasma kinetic scales, measured by the Cluster/STAFF instrument in the free solar wind. We show that the magnetic spectrum in the high frequency range, [1,400] Hz, has a form similar to what is found in hydrodynamics in the dissipation range ~Ak^(-\alpha)exp(-kl_d). The dissipation scale l_d is found to be correlated with the electron Larmor radius \rho_e. The spectral index \alpha varies in the range [2.2,2.9] and is anti-correlated with l_d, as expected in the case of the balance between the energy injection and the energy dissipation. The coefficient A is found to be proportional to the ion temperature anisotropy, suggesting that local ion instabilities may play some role for the solar wind turbulence at plasma kinetic scales. The exponential spectral shape found here indicates that the effective dissipation of magnetic fluctuations in the solar wind has a wave number dependence similar to that of the resistive term in collisional fluids ~k^2\delta B.

Coronal loops act as resonant cavities for low frequency fluctuations that are transmitted from the deeper layers of the solar atmosphere and are amplified in the corona, triggering nonlinear interactions. However trapping is not perfect, some energy leaks down to the chromosphere, thus limiting the turbulence development and the associated heating. We consider the combined effects of turbulence and leakage in determining the energy level and associated heating rate in models of coronal loops which include the chromosphere and transition region. We use a piece-wise constant model for the Alfven speed and a Reduced MHD - Shell model to describe the interplay between turbulent dynamics in the direction perpendicular to the mean field and propagation along the field. Turbulence is sustained by incoming fluctuations which are equivalent, in the line-tied case, to forcing by the photospheric shear flows. While varying the turbulence strength, we compare systematically the average coronal energy level (E) and dissipation rate (D) in three models with increasing complexity: the classical closed model, the semi-open corona model, and the corona-chromosphere (or 3-layer) model, the latter two models allowing energy leakage. We find that: (i) Leakage always plays a role (even for strong turbulence), E and D are systematically lower than in the line-tied model. (ii) E is close to the resonant prediction, i.e., assuming effective turbulent correlation time longer than the Alfven coronal crossing time (Ta). (iii) D is close to the value given by the ratio of photospheric energy divided by Ta (iv) The coronal spectra exibits an inertial range with 5/3 spectral slope, and a large scale peak of trapped resonant modes that inhibit nonlinear couplings. (v) In the realistic 3-layer model, the two-component spectrum leads to a damping time equal to the Kolmogorov time reduced by a factor u_rms/Va_corona

We study the connections between the sun's convection zone and the evolution of the solar wind and corona. We let the magnetic fields generated by a 2.5D axisymmetric kinematic dynamo code (STELEM) evolve in a 2.5D axisymmetric coronal isothermal MHD code (DIP). The computations cover an 11 year activity cycle. The solar wind's asymptotic velocity varies in latitude and in time in good agreement with the available observations. The magnetic polarity reversal happens at different paces at different coronal heights. Overall sun's mass loss rate, momentum flux and magnetic braking torque vary considerably throughout the cycle. This cyclic modulation is determined by the latitudinal distribution of the sources of open flux and solar wind and the geometry of the Alfv\'en surface. Wind sources and braking torque application zones also vary accordingly.

We present the first results of simulations of giant polar plumes and coronal flows. We use a 2.5D axisymmetric MHD numerical model of an isothermal corona and slow solar wind. A plume is generated just above asmall magnetic bipole embedded in an unipolar flux region which is perturbed by Alfv\'en waves injected from the coronal base. The boundary conditions are transparent. The results are compared to those obtained previously with a 1D wind model in which plumes are generated as a consequence of variations of the heating and flux-tube expansion parameters.

In a recent work Grappin et al. [1] have shown that low- frequency movements can be transmitted from one footpoint to the other along a magnetic loop, thus mimicking a friction effect of the corona on the photosphere, and invalidating the line-tying approximation. We consider here successively the effect of high frequencies and turbulent damping on the process. We use a very simple atmospheric model which allows to study analytically the laminar case, and to study the turbulent case both using simple phenomenological arguments and a more sophisticated turbulence model [2]. We find that, except when turbulent damping is such that all turbulence is damped during loop traversal, coupling still occurs between distant footpoints, and moreover the coronal field induced by photospheric movements saturates at finite values.

We investigate how wave propagation is modified by the presence of heat sources and sinks, in the simple 1D, hydrodynamical case, including chromosphere and solar wind. We integrate the time-dependent hydrodynamic equations of the solar wind with spherical symmetry, including conduction, radiative cooling and a prescribed mechanical heat flux. Once a quasi-stationary wind is established, we study the response of the system to pressure oscillations at the photospheric boundary. We use transparent boundary conditions. We find that wavepackets with high enough amplitude propagating upward from the photosphere implode just below the transition region. This implosion is due to the radiative cooling term generating pressure holes close to the wave crests of the wave, which make the wave collapse. In the case where heat sources and sinks are not present in the equations, the wave remains stable whatever the initial wave amplitude, which is compatible with published work. Instability should be observable when and where the TR is high enough above the optically thick regions.