The question of the relative importance of coherent structures and waves has for a long time attracted a great deal of interest in astrophysical plasma turbulence research, with a more recent focus on kinetic scale dynamics. Here we utilize high-resolution observational and simulation data to investigate the nature of waves and structures emerging in a weakly collisional, turbulent kinetic plasma. Observational results are based on in situ solar wind measurements from the Cluster and MMS spacecraft, and the simulation results are obtained from an externally driven, three-dimensional fully kinetic simulation. Using a set of novel diagnostic measures we show that both the large-amplitude structures and the lower-amplitude background fluctuations preserve linear features of kinetic Alfven waves to order unity. This quantitative evidence suggests that the kinetic turbulence cannot be described as a mixture of mutually exclusive waves and structures but may instead be pictured as an ensemble of localized, anisotropic wave packets or "eddies" of varying amplitudes, which preserve certain linear wave properties during their nonlinear evolution.

We present results from a three-dimensional particle-in-cell simulation of plasma turbulence, resembling the plasma conditions found at kinetic scales of the solar wind. The spectral properties of the turbulence in the subion range are consistent with theoretical expectations for kinetic Alfv\' en waves. Furthermore, we calculate the local anisotropy, defined by the relation $k_{\parallel}(k_{\perp})$, where $k_{\parallel}$ is a characteristic wave number along the local mean magnetic field at perpendicular scale $l_{\perp}\sim 1/k_{\perp}$. The subion range anisotropy is scale dependent with $k_{\parallel}<k_{\perp}$ and the ratio of linear to nonlinear time scales is of order unity, suggesting that the kinetic cascade is close to a state of critical balance. Our results compare favorably against a number of \emph{in situ} solar wind observations and demonstrate---from first principles---the feasibility of plasma turbulence models based on a critically balanced cascade of kinetic Alfv\' en waves.

We report the results of a direct comparison between different kinetic models of collisionless plasma turbulence in two spatial dimensions. The models considered include a first principles fully-kinetic (FK) description, two widely used reduced models [gyrokinetic (GK) and hybrid-kinetic (HK) with fluid electrons], and a novel reduced gyrokinetic approach (KREHM). Two different ion beta ($\beta_i$) regimes are considered: 0.1 and 0.5. For $\beta_i=0.5$, good agreement between the GK and FK models is found at scales ranging from the ion to the electron gyroradius, thus providing firm evidence for a kinetic Alfv\'en cascade scenario. In the same range, the HK model produces shallower spectral slopes, presumably due to the lack of electron Landau damping. For $\beta_i=0.1$, a detailed analysis of spectral ratios reveals a slight disagreement between the GK and FK descriptions at kinetic scales, even though kinetic Alfv\'en fluctuations likely still play a significant role. The discrepancy can be traced back to scales above the ion gyroradius, where the FK and HK results seem to suggest the presence of fast magnetosonic and ion Bernstein modes in both plasma beta regimes, but with a more notable deviation from GK in the low-beta case. The identified practical limits and strengths of reduced-kinetic approximations, compared here against the fully-kinetic model on a case-by-case basis, may provide valuable insight into the main kinetic effects at play in turbulent collisionless plasmas, such as the solar wind.