Systematic analysis of numerical simulations of two-dimensional magnetohydrodynamic turbulence reveals the presence of a large number of $X$-type neutral points where magnetic reconnection occurs. We examine the statistical properties of this ensemble of reconnection events that are spontaneously generated by turbulence. The associated reconnection rates are distributed over a wide range of values and scales with the geometry of the diffusion region. Locally, these events can be described through a variant of the Sweet-Parker model, in which the parameters are externally controlled by turbulence. This new perspective on reconnection is relevant in space and astrophysical contexts, where plasma is generally in a fully turbulent regime.

An overview is given of important properties of spatial and temporal intermittency, including evidence of its appearance in fluids, magnetofluids and plasmas, and its implications for understanding of heliospheric plasmas. Spatial intermittency is generally associated with formation of sharp gradients and coherent structures. The basic physics of structure generation is ideal, but when dissipation is present it is usually concentrated in regions of strong gradients. This essential feature of spatial intermittency in fluids has been shown recently to carry over to the realm of kinetic plasma, where the dissipation function is not known from first principles. Spatial structures produced in intermittent plasma influence dissipation, heating, and transport and acceleration of charged particles. Temporal intermittency can give rise to very long time correlations or a delayed approach to steady-state conditions, and has been associated with inverse cascade or quasi-inverse cascade systems, with possible implications for heliospheric prediction.

We examine statistics of rapid spatial variations of the magnetic field in simulations of magnetohydrodynamic (MHD) turbulence, by analyzing intermittency properties, and by using classical methods for identifying discontinuities. The methods identify similar structures, and give very similar event distribution functions. When the results are scaled to the correlation length, the average waiting times agree with typically reported waiting times between solar wind discontinuities. Thus discontinuities may be related to flux tube boundaries and intermittent structures that appear spontaneously in MHD turbulence.

Using direct numerical simulations of a hybrid Vlasov-Maxwell model, kinetic processes are investigated in a two-dimensional turbulent plasma. In the turbulent regime, kinetic effects manifest through a deformation of the ion distribution function. These patterns of non-Maxwellian features are concentrated in space nearby regions of strong magnetic activity: the distribution function is modulated by the magnetic topology, and can elongate along or across the local magnetic field. These results open a new path on the study of kinetic processes such as heating, particle acceleration, and temperature~anisotropy, commonly observed in astrophysical and laboratory plasmas.