Let ZN, N >= 1, denote the integer lattice points in the N-dimensional Euclidean space. Asymptotic normality of kernel estimators of the multivariate density of stationary random fields indexed by ZN is established. Appropriate choices of the bandwiths are found. The random fields are assumed to satisfy some mixing conditions. The results apply to many spatial models.