In this talk we will briefly introduce models accounting for the correlation over the spacetime domain.

We consider recent results for a class of non-separable spacetime models and to outline a computational implementation including a simple way to introduce non-stationarity.

We discuss some practical details through a working example considering a new R package dedicated to this work.