Linear mixed effects models have been widely used in the spatial analysis of environmental processes. However, parameter estimation and spatial predictions involve the inversion and determinant of the n times n dimensional spatial covariance matrix of the data process, with n being the number of...