Statistical properties of scale-dependent velocity and magnetic field fluctuations decomposed into components parallel and perpendicular to the local magnetic field in different ISM phases. Fluctuations are computed as the square root of the second-order structure function. They are shown separately for the WNM, UNM, CNM, and total ISM (including all phases). The shaded areas represent the standard deviation over the realizations. Dashed and dashed–dotted gray lines represent reference slopes of 1/3 and 1/2, corresponding to Kolmogorov and Burgers turbulence scalings, respectively. The numerical dissipation scale r_dis is around 0.5 pc.