The space-wise approach to GOCE data reduction exploits the spatial correlation of the observations by “projecting” them on a spherical grid at mean satellite altitude; to this aim a local collocation prediction based on a global covariance function of the potential T can be used. Since the implicit hypothesis of rotational invariance is not fulfilled by the gravitational potential, the global covariance function cannot describe the local characteristics of the field in the different interpolation areas. As a consequence the estimation error of the gridded values is not homogeneously distributed all over the reference sphere, but it is much higher over the Himalaya, the Alps, the Andes, etc., i.e. over areas where the random field presents ...