Small area estimation is studied under a spatio-temporal Fay–Herriot model. Model fitting based on restricted maximum likelihood is described and empirical best linear unbiased predictors are derived under the model. A parametric bootstrap procedure is proposed for the estimation of the mean...