Problems with likelihood estimation of covariance functions of spatial Gaussian processes

Maximum likelihood has frequently been suggested as a way to estimate covariance parameters in spatial Gaussian processes. The results of some simple examples based on simulated data are discussed. It is shown by an example with real data that the profile likelihood can be multimodal, and that the global maximum may not correspond to a sensible value of the parameters.