Magnetospheric Multiscale (MMS) observations of plasma turbulence generated by a Kelvin-Helmholtz (KH) event at the Earth's magnetopause are compared with a high-resolution two-dimensional (2D) hybrid direct numerical simulation (DNS) of decaying plasma turbulence driven by large-scale balanced Alfv\'enic fluctuations. The simulation, set up with four observation-driven physical parameters (ion and electron betas, turbulence strength, and injection scale) exhibits a quantitative agreement on the spectral, intermittency, and cascade-rate properties with in situ observations, despite the different driving mechanisms. Such agreement demonstrates a certain universality of the turbulent cascade from magnetohydrodynamic (MHD) to sub-ion scales, whose properties are mainly determined by the selected parameters, also indicating that the KH instability-driven turbulence has a quasi-2D nature. The validity of the Taylor hypothesis in the sub-ion spatial range suggests that the fluctuations at sub-ion scales have predominantly low frequencies, consistent with a kinetic Alfv\'en wave-like nature or with quasi-static structures. Finally, the third-order structure function analysis indicates that the cascade rate of the turbulence generated by a KH event in the magnetopause is an order of magnitude larger than in the ambient magnetosheath.

The relationship between a decaying plasma turbulence and proton fire hose instabilities in a slowly expanding plasma is investigated using three-dimensional (3-D) hybrid expanding box simulations. We impose an initial ambient magnetic field along the radial direction, and we start with an isotropic spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations with zero cross-helicity. A turbulent cascade rapidly develops and leads to a weak proton heating that is not sufficient to overcome the expansion-driven perpendicular cooling. The plasma system eventually drives the parallel and oblique fire hose instabilities that generate quasi-monochromatic wave packets that reduce the proton temperature anisotropy. The fire hose wave activity has a low amplitude with wave vectors quasi-parallel/oblique with respect to the ambient magnetic field outside of the region dominated by the turbulent cascade and is discernible in one-dimensional power spectra taken only in the direction quasi-parallel/oblique with respect to the ambient magnetic field; at quasi-perpendicular angles the wave activity is hidden by the turbulent background. These waves are partly reabsorbed by protons and partly couple to and participate in the turbulent cascade. Their presence reduces kurtosis, a measure of intermittency, and the Shannon entropy but increases the Jensen-Shannon complexity of magnetic fluctuations; these changes are weak and anisotropic with respect to the ambient magnetic field and it's not clear if they can be used to indirectly discern the presence of instability-driven waves.

By means of three dimensional high-resolution hybrid simulations we study the properties of the magnetic field spectral anisotropies near and beyond ion kinetic scales. By using both a Fourier analysis and a local analysis based on multi-point 2nd-order structure function techniques, we show that the anisotropy observed is less than what expected by standard wave normal modes turbulence theories although the non linear energy transfer is still in the perpendicular direction, only advected in the parallel direction as expected balancing the non-linear energy transfer time and the decorrelation time. Such result can be explained by a phenomenological model based on the formation of strong intermittent two-dimensional structures in the plane perpendicular to the local mean field that have some prescribed aspect ratio eventually depending on the scale. This model support the idea that small scales structures, such as reconnecting current sheets, contribute significantly to the formation of the turbulent cascade at kinetic scales.

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.

This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence.

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.

A dynamical vectorial equation for homogeneous incompressible Hall-MHD turbulence together with the exact scaling law for third-order correlation tensors, analogous to that for the incompressible MHD, is rederived and applied to the results of two-dimensional hybrid simulations of plasma turbulence. At large (MHD) scales the simulations exhibits a clear inertial range where the MHD dynamic law is valid. In the sub-ion range the cascade continues via the Hall term but the dynamic law derived in the framework of incompressible Hall MHD equations is obtained only in a low plasma beta simulation. For a higher beta plasma the cascade rate decreases in the sub-ion range and the change becomes more pronounced as the plasma beta increases. This break in the cascade flux can be ascribed to non thermal (kinetic) features or to others terms in the dynamical equation that are not included in the Hall-MHD incompressible approximation.

We investigate the spectral properties of plasma turbulence from fluid to sub-ion scales by means of high-resolution three-dimensional (3D) numerical simulations performed with the hybrid particle-in-cell (HPIC) code CAMELIA. We produce extended turbulent spectra with well-defined power laws for the magnetic, ion bulk velocity, density, and electric fluctuations. The present results are in good agreement with previous two-dimensional (2D) HPIC simulations, especially in the kinetic range of scales, and reproduce several features observed in solar wind spectra. By providing scaling tests on many different architectures and convergence studies, we prove CAMELIA to represent a very efficient, accurate and reliable tool for investigating the develpoment of the turbulent cascade in the solar wind, being able to cover simultaneously several decades in wavenumber, also in 3D.

Spectral properties of the turbulent cascade from fluid to kinetic scales in collisionless plasmas are investigated by means of large-size three-dimensional (3D) hybrid (fluid electrons, kinetic protons) particle-in-cell simulations. Initially isotropic Alfv\`enic fluctuations rapidly develop a strongly anisotropic turbulent cascade, mainly in the direction perpendicular to the ambient magnetic field. The omnidirectional magnetic field spectrum shows a double power-law behavior over almost two decades in wavenumber, with a Kolmogorov-like index at large scales, a spectral break around ion scales, and a steepening at sub-ion scales. Power laws are also observed in the spectra of the ion bulk velocity, density, and electric field, both at magnetohydrodynamic (MHD) and at kinetic scales. Despite the complex structure, the omnidirectional spectra of all fields at ion and sub-ion scales are in remarkable quantitative agreement with those of a two-dimensional (2D) simulation with similar physical parameters. This provides a partial, a-posteriori validation of the 2D approximation at kinetic scales. Conversely, at MHD scales, the spectra of the density and of the velocity (and, consequently, of the electric field) exhibit differences between the 2D and 3D cases. Although they can be partly ascribed to the lower spatial resolution, the main reason is likely the larger importance of compressible effects in a full geometry. Our findings are also in remarkable quantitative agreement with solar wind observations.

Motivated by recent observational searches of sub-structure in starless molecular cloud cores, we investigate the evolution of density perturbations on scales smaller than the Jeans length embedded in contracting isothermal clouds, adopting the same formalism developed for the expanding Universe and the solar wind. We find that initially small amplitude, Jeans-stable perturbations (propagating as sound waves in the absence of a magnetic field), are amplified adiabatically during the contraction, approximately conserving the wave action density, until they either become nonlinear and steepen into shocks at a time $t_{\rm nl}$, or become gravitationally unstable when the Jeans length decreases below the scale of the perturbations at a time $t_{\rm gr}$. We evaluate analytically the time $t_{\rm nl}$ at which the perturbations enter the non-linear stage using a Burgers' equation approach, and we verify numerically that this time marks the beginning of the phase of rapid dissipation of the kinetic energy of the perturbations. We then show that for typical values of the rms Mach number in molecular cloud cores, $t_{\rm nl}$ is smaller than $t_{\rm gr}$, and therefore density perturbations likely dissipate before becoming gravitational unstable. Solenoidal modes grow at a faster rate than compressible modes, and may eventually promote fragmentation through the formation of vortical structures.

A new path for the generation of a sub-ion scale cascade in collisionless plasma turbulence, triggered by magnetic reconnection, is uncovered by means of high-resolution two-dimensional hybrid-kinetic simulations employing two complementary approaches, Lagrangian and Eulerian, and different driving mechanisms. The simulation results provide clear numerical evidences that the development of power-law energy spectra below the so-called ion break occurs as soon as the first magnetic reconnection events take place, regardless of the actual state of the turbulent cascade at MHD scales. In both simulations, the reconnection-mediated small-scale energy spectrum of parallel magnetic fluctuations exhibits a very stable spectral slope of ~-2.8, whether or not a large-scale turbulent cascade has already fully developed. Once a quasi-stationary turbulent state is achieved, the spectrum of the total magnetic fluctuations settles towards a spectral index of -5/3 in the MHD range and of ~-3 at sub-ion scales.

The relationship between a decaying strong turbulence and the mirror instability in a slowly expanding plasma is investigated using two-dimensional hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations which have energy equipartition between kinetic and magnetic fluctuations and vanishing correlation between the two fields. A turbulent cascade rapidly develops, magnetic field fluctuations exhibit a Kolmogorov-like power-law spectrum at large scales and a steeper spectrum at sub-ion scales. The imposed expansion (taking a strictly transverse ambient magnetic field) leads to generation of an important perpendicular proton temperature anisotropy that eventually drives the mirror instability. This instability generates large-amplitude, nonpropagating, compressible, pressure-balanced magnetic structures in a form of magnetic enhancements/humps that reduce the perpendicular temperature anisotropy.

We investigate properties of the ion-scale spectral break of solar wind turbulence by means of two-dimensional high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field perpendicular to the simulation box and add a spectrum of in-plane, large-scale, magnetic and kinetic fluctuations. We perform a set of simulations with different values of the plasma beta, distributed over three orders of magnitude, from 0.01 to 10. In all the cases, once turbulence is fully developed, we observe a power-law spectrum of the fluctuating magnetic field on large scales (in the inertial range) with a spectral index close to -5/3, while in the sub-ion range we observe another power-law spectrum with a spectral index systematically varying with $\beta$ (from around -3.6 for small values to around -2.9 for large ones). The two ranges are separated by a spectral break around ion scales. The length scale at which this transition occurs is found to be proportional to the ion inertial length, $d_i$, for $\beta \ll 1$ and to the ion gyroradius, $\rho_i = d_i\sqrt{\beta}$, for $\beta \gg 1$, i.e., to the larger between the two scales in both the extreme regimes. For intermediate cases, i.e., $\beta \sim 1$, a combination of the two scales is involved. We infer an empiric relation for the dependency of the spectral break on $\beta$ that provides a good fit over the whole range of values. We compare our results with in situ observations in the solar wind and suggest possible explanations for such a behavior.

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.

Proton temperature anisotropies between the directions parallel and perpendicular to the mean magnetic field are usually observed in the solar wind plasma. Here, we employ a high-resolution hybrid particle-in-cell simulation in order to investigate the relation between spatial properties of the proton temperature and the peaks in the current density and in the flow vorticity. Our results indicate that, although regions where the proton temperature is enhanced and temperature anisotropies are larger correspond approximately to regions where many thin current sheets form, no firm quantitative evidence supports the idea of a direct causality between the two phenomena. On the other hand, quite a clear correlation between the behavior of the proton temperature and the out-of-plane vorticity is obtained.

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$.

The relationship between a decaying strong turbulence and kinetic instabilities in a slowly expanding plasma is investigated using two-dimensional (2-D) hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfv\'enic fluctuations which have energy equipartition between kinetic and magnetic fluctuations and vanishing correlation between the two fields. A turbulent cascade rapidly develops, magnetic field fluctuations exhibit a Kolmogorov-like power-law spectrum at large scales and a steeper spectrum at ion scales. The turbulent cascade leads to an overall anisotropic proton heating, protons are heated in the perpendicular direction, and, initially, also in the parallel direction. The imposed expansion leads to generation of a large parallel proton temperature anisotropy which is at later stages partly reduced by turbulence. The turbulent heating is not sufficient to overcome the expansion-driven perpendicular cooling and the system eventually drives the oblique firehose instability in a form of localized nonlinear wave packets which efficiently reduce the parallel temperature anisotropy. This work demonstrates that kinetic instabilities may coexist with strong plasma turbulence even in a constrained 2-D regime.

We investigate properties of plasma turbulence from magneto-hydrodynamic (MHD) to sub-ion scales by means of two-dimensional, high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field, perpendicular to the simulation box, and we add a spectrum of large-scale magnetic and kinetic fluctuations, with energy equipartition and vanishing correlation. Once the turbulence is fully developed, we observe a MHD inertial range, where the spectra of the perpendicular magnetic field and the perpendicular proton bulk velocity fluctuations exhibit power-law scaling with spectral indices of -5/3 and -3/2, respectively. This behavior is extended over a full decade in wavevectors and is very stable in time. A transition is observed around proton scales. At sub-ion scales, both spectra steepen, with the former still following a power law with a spectral index of ~-3. A -2.8 slope is observed in the density and parallel magnetic fluctuations, highlighting the presence of compressive effects at kinetic scales. The spectrum of the perpendicular electric fluctuations follows that of the proton bulk velocity at MHD scales, and flattens at small scales. All these features, which we carefully tested against variations of many parameters, are in good agreement with solar wind observations. The turbulent cascade leads to on overall proton energization with similar heating rates in the parallel and perpendicular directions. While the parallel proton heating is found to be independent on the resistivity, the number of particles per cell and the resolution employed, the perpendicular proton temperature strongly depends on these parameters.

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.

We present results from a high-resolution and large-scale hybrid (fluid electrons and particle-in-cell protons) two-dimensional numerical simulation of decaying turbulence. Two distinct spectral regions (separated by a smooth break at proton scales) develop with clear power-law scaling, each one occupying about a decade in wave numbers. The simulation results exhibit simultaneously several properties of the observed solar wind fluctuations: spectral indices of the magnetic, kinetic, and residual energy spectra in the magneto-hydrodynamic (MHD) inertial range along with a flattening of the electric field spectrum, an increase in magnetic compressibility, and a strong coupling of the cascade with the density and the parallel component of the magnetic fluctuations at sub-proton scales. Our findings support the interpretation that in the solar wind large-scale MHD fluctuations naturally evolve beyond proton scales into a turbulent regime that is governed by the generalized Ohm's law.

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 long-term evolution of large-amplitude Alfven waves propagating in the solar wind is investigated by performing two-dimensional MHD simulations within the expanding box model. The linear and nonlinear phases of the parametric decay instability are studied for both circularly polarized waves in parallel propagation and for arc-polarized waves in oblique propagation. The non-monochromatic case is also considered. In the oblique case, the direct excitation of daughter modes transverse to the local background field is found for the first time in an expanding environment, and this transverse cascade seems to be favored for monochromatic mother waves. The expansion effect reduces the instability growth rate, and it can even suppress its onset for the lowest frequency modes considered here, possibly explaining the persistence of these outgoing waves in the solar wind.

Although the mechanisms responsible for heating the Sun's corona and accelerating the solar wind are still being actively investigated, it is largely accepted that photospheric motions provide the energy source and that the magnetic field must play a key role in the process. \citet{2010ApJ...708L.116V} presented a model for heating and accelerating the solar wind based on the turbulent dissipation of Alfv\'en waves. We first use a time-dependent model of the solar wind to reproduce one of \citeauthor{2010ApJ...708L.116V}'s solutions; then we extend its application to the case when the energy equation includes thermal conduction and radiation losses, and the upper chromosphere is part of the computational domain. Using this model, we explore parameter space and describe the characteristics of a fast-solar-wind solution. We discuss how this formulation may be applied to a 3D MHD model of the coron a and solar wind \citep{2009ApJ...690..902L}.

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.

Context: Supersonic nonthermal motions in molecular clouds are often interpreted as long-lived magnetohydrodynamic (MHD) waves. The propagation and amplitude of these waves is affected by local physical characteristics, most importantly the gas density and the ionization fraction. Aims: We study the propagation, reflection and dissipation of Alfv\'en waves in molecular clouds deriving the behavior of observable quantities such as the amplitudes of velocity fluctuations and the rate of energy dissipation. Methods: We formulated the problem in terms of Els\"asser variables for transverse MHD waves propagating in a one-dimensional inhomogeneous medium, including the dissipation due to collisions between ions and neutrals and to a nonlinear turbulent cascade treated in a phenomenological way. We considered both steady-state and time-dependent situations and solved the equations of the problem numerically with an iterative method and a Lax-Wendroff scheme, respectively. Results: Alfv\'en waves incident on overdense regions with density profiles typical of cloud cores embedded in a diffuse gas suffer enhanced reflection in the regions of the steepest density gradient, and strong dissipation in the core's interior. These effects are especially significant when the wavelength is intermediate between the critical wavelength for propagation and the typical scale of the density gradient. For larger wave amplitudes and/or steeper input spectra, the effects of the perpendicular turbulent cascade result in a stronger energy dissipation in the regions immediately surrounding the dense core. Conclusions: The results may help to interpret the sharp decrease of line width observed in the environments of low-mass cloud cores in several molecular transitions.

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).

We report observations of quasi-periodic pulsations (QPPs) during the X2.2 flare of 2011 February 15, observed simultaneously in several wavebands. We focus on fluctuations on time scale 1-30 s and find different time lags between different wavebands. During the impulsive phase, the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) channels in the range 25-100 keV lead all the other channels. They are followed by the Nobeyama RadioPolarimeters at 9 and 17 GHz and the Extreme Ultra-Violet (EUV) channels of the Euv SpectroPhotometer (ESP) onboard the Solar Dynamic Observatory (SDO). The Zirconium and Aluminum filter channels of the Large Yield Radiometer (LYRA) onboard the Project for On-Board Autonomy (PROBA2) satellite and the SXR channel of ESP follow. The largest lags occur in observations from the Geostationary Operational Environmental Satellite (GOES), where the channel at 1-8 {\AA} leads the 0.5-4 {\AA} channel by several seconds. The time lags between the first and last channels is up to 9 s. We identified at least two distinct time intervals during the flare impulsive phase, during which the QPPs were associated with two different sources in the Nobeyama RadioHeliograph at 17 GHz. The radio as well as the hard X-ray channels showed different lags during these two intervals. To our knowledge, this is the first time that time lags are reported between EUV and SXR fluctuations on these time scales. We discuss possible emission mechanisms and interpretations, including flare electron trapping.

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

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.

A model is presented for generation of fast solar wind in coronal holes, relying on heating that is dominated by turbulent dissipation of MHD fluctuations transported upwards in the solar atmosphere. Scale-separated transport equations include large-scale fields, transverse Alfvenic fluctuations, and a small compressive dissipation due to parallel shears near the transition region. The model accounts for proton temperature, density, wind speed, and fluctuation amplitude as observed in remote sensing and in situ satellite data.

We study the formation and evolution of a turbulent spectrum of Alfv\'en waves driven by reflection off the solar wind density gradients, starting from the coronal base up to 17 solar radii, well beyond the Alfv\'enic critical point. The background solar wind is assigned and 2D shell models are used to describe nonlinear interactions. We find that the turbulent spectra are influenced by the nature of reflected waves. Close to the base, these give rise to a flatter and steeper spectrum for the outgoing and reflected waves respectively. At higher heliocentric distance both spectra evolve toward an asymptotic Kolmogorov spectrum. The turbulent dissipation is found to account for at least half of the heating required to sustain the background imposed solar wind and its shape is found to be determined by the reflection-determined turbulent heating below 1.5 solar radii. Therefore reflection and reflection-driven turbulence are shown to play a key role in the accelerationof the fast solar wind and origin of the turbulent spectrum found at 0.3 AU in the heliosphere.

We solve the problem of propagation and dissipation of Alfvenic turbulence in a model solar atmosphere consisting of a static photosphere and chromosphere, transition region, and open corona and solar wind, using a phenomenological model for the turbulent dissipation based on wave reflection. We show that most of the dissipation for a given wave-frequency spectrum occurs in the lower corona, and the overall rms amplitude of the fluctuations evolves in a way consistent with observations. The frequency spectrum, for a Kolmogorov-like slope, is not found to change dramatically from the photosphere to the solar wind, however it does preserve signatures of transmission throughout the lower atmospheric layers, namely oscillations in the spectrum at high frequencies reminiscent of the resonances found in the linear case. These may disappear once more realistic couplings for the non-linear terms are introduced, or if time-dependent variability of the lower atmospheric layer is introduced.