In this study, we explore the spectral properties of the distribution matrices of the mapping method and its relation to the distributive mixing of passive scalars. The spectral (or eigenvector-eigenvalue) decomposition of these matrices constitutes discrete approximations to the eigenmodes of the continuous advection operator in periodic flows. The eigenvalue spectrum always lies within the unit circle and due to mass conservation, always accommodates an eigenvalue equal to one with trivial (uniform) eigenvector. The asymptotic state of a fully chaotic mixing flow is dominated by the eigenmode corresponding with the eigenvalue closest to the unit circle ( dominant eigenmode ). This eigenvalue determines the decay rate; its eigenvector dete...