Given a general homogeneous non-stationary autoregressive integrated moving average process ARIMA(p,d,q), the corresponding model for the subseries obtained by a systematic sampling is derived. The article then shows that the sampled subseries approaches approximately to an integrated moving average process IMA(d,l), l≤(d-l), regardless of the autoregressive and moving average structures in the original series. In particular, the sampled subseries from an ARIMA (p,l,q) process approaches approximately to a simple random walk model.