Giovanna, space-time cross covariance has never seemed to me something that can be done "simple", but package gstat at least offers you: 1. to include time as a third dimension, and model 3D anisotropy, 2. to calculate cross variograms between different moments in time, and proceed with cokriging. -- Edzer