In time series analysis, it is well-known that the differencing operator ∇d may transform a non-stationary series, {Z(t)} say, to a stationary one, {W(t)} = ∇dZ(t)}; and there are many procedures for analysing and modelling {Z(t)} which exploit this transformation. Rather differently, Matheron (1973) introduced a set of measures on Rn that transform an appropriate non-stationary spatial process to stationarity, and Cressie (1988) then suggested that specialized low-order analogues of these measures, called increment-vectors, be used in time series analysis. This paper develops a general theory of increment-vectors which provides a more powerful transformation tool than mere simple differencing. The methodology gives a handle on the second-moment structure and divergence behaviour of homogeneously non-stationary series which leads to many important applications such as determining the correct degree of differencing, forecasting and interpolation.

The intrinsic random functions (IRF) are a particular case of the Guelfand generalized processes with stationary increments. They constitute a much wider class than the stationary RF, and are used in practical applications for representing non-stationary phenomena. The most important topics are: existence of a generalized covariance (GC) for which statistical inference is possible from a unique realization; theory of the best linear intrinsic estimator (BLIE) used for contouring and estimating problems; the turning bands method for simulating IRF; and the models with polynomial GC, for which statistical inference may be performed by automatic procedures.