Magnetic reconnection (MR) and the associated concurrently occurring waves have been extensively studied at large-scale plasma boundaries, in quasi-symmetric and asymmetric configurations in the terrestrial magnetotail and at the magnetopause. Recent high-resolution observations by MMS (Magnetospheric Multiscale) spacecraft indicate that MR can occur also in the magnetosheath where the conditions are highly turbulent when the upstream shock geometry is quasi-parallel. The strong turbulent motions make the boundary conditions for evolving MR complicated. In this paper it is demonstrated that the wave observations in localized regions of MR can serve as an additional diagnostic tool reinforcing our capacity for identifying MR events in turbulent plasmas. It is shown that in a close resemblance with MR at large-scale boundaries, turbulent reconnection associated whistler waves occur at separatrix/outflow regions and at the outer boundary of the electron diffusion region, while lower hybrid drift waves are associated with density gradients during the crossing of the current sheet. The lower hybrid drift instability can make the density inhomogeneities rippled. The identification of MR associated waves in the magnetosheath represents also an important milestone for developing a better understanding of energy redistribution and dissipation in turbulent plasmas.

In the present paper, we investigate the power-law behaviour of the magnetic field spectra in the Earths magnetosheath region using Cluster spacecraft data under solar minimum condition. The power spectral density of the magnetic field data and spectral slopes at various frequencies are analysed. Propagation angle and compressibility are used to test the nature of turbulent fluctuations. The magnetic field spectra have the spectral slopes between -1.5 to 0 down to spatial scales of 20 ion gyroradius and show clear evidence of a transition to steeper spectra for small scales with a second power-law, having slopes between -2.6 to -1.8. At low frequencies, f_sc<0.3f_ci(where f_ci is ion gyro-frequency), propagation angle approximately 90 degrees to the mean magnetic field, B_0, and compressibility shows a broad distribution, 0.1 < R > 0.9. On the other hand at f_sc>10f_ci, the propagation angle exhibits a broad range between 30-90 degree while 'R' has a small variation: 0.2 < R > 0.5. We conjecture that at high frequencies, the perpendicularly propagating Alfven waves could partly explain the statistical analysis of spectra. To support our prediction of kinetic Alfven wave-dominated spectral slope behaviour at high frequency, we also present a theoretical model and simulate the magnetic field turbulence spectra due to the nonlinear evolution of kinetic Alfven waves. The present study also shows the analogy between the observational and simulated spectra.

In weakly collisional space plasmas, the turbulent cascade provides most of the energy that is dissipated at small scales by various kinetic processes. Understanding the characteristics of such dissipative mechanisms requires the accurate knowledge of the fluctuations that make energy available for conversion at small scales, as different dissipation processes are triggered by fluctuations of a different nature. The scaling properties of different energy channels are estimated here using a proxy of the local energy transfer, based on the third-order moment scaling law for magnetohydrodynamic turbulence. In particular, the sign-singularity analysis was used to explore the scaling properties of the alternating positive-negative energy fluxes, thus providing information on the structure and topology of such fluxes for each of the different type of fluctuations. The results show the highly complex geometrical nature of the flux, and that the local contributions associated with energy and cross-helicity nonlinear transfer have similar scaling properties. Consequently, the fractal properties of current and vorticity structures are similar to those of the Alfv\'enic fluctuations.

The degree of deviation from the thermodynamic equilibrium in the ion velocity distribution functions (VDFs), measured by the Magnetospheric Multiscale (MMS) mission in the Earth's turbulent magnetosheath, is quantitatively investigated. Taking advantage of MMS ion data, having a resolution never reached before in space missions, and of the comparison with Vlasov-Maxwell simulations, this analysis aims at relating any deviation from Maxwellian equilibrium to typical plasma parameters. Correlations of the non-Maxwellian features with plasma quantities such as electric fields, ion temperature, current density and ion vorticity are very similar in both magnetosheath data and numerical experiments, and suggest that distortions in the ion VDFs occur close to (but not exactly at) peaks in current density and ion temperature. Similar results have also been found during a magnetopause crossing by MMS. This work could help clarifying the origin of distortion of the ion VDFs in space plasmas.

Using high resolution Cluster satellite observations, we show that the turbulent solar wind is populated by magnetic discontinuities at different scales, going from proton down to electron scales. The structure of these layers resembles the Harris equilibrium profile in plasmas. Using a multi-dimensional intermittency technique, we show that these structures are connected through the scales. Supported by numerical simulations of magnetic reconnection, we show that observations are consistent with a scenario where many current layers develop in turbulence, and where the outflow of these reconnection events are characterized by complex sub-proton networks of secondary islands, in a self-similar way. The present work establishes that the picture of "reconnection in turbulence" and "turbulent reconnection", separately invoked as ubiquitous, coexist in space plasmas.