This paper is aimed at the problem of predicting the land subsidence or upheave in an area, using GNSS position time series. Since machine learning algorithms have presented themselves as strong prediction tools in different fields of science, we employ them to predict the next values of the GNSS position time series. For this reason, we present an algorithm that takes advantage of the machine learning algorithms for the prediction of positions in a GNSS time series. The proposed algorithm has two steps-preprocessing and prediction. In the preprocessing phase, the periodic tidal and atmospheric signals in the time series are removed and coordinates are transferred to the local coordinate system. In the prediction phase, eight different machine learning algorithms are used, namely, multilayer perceptron, Bayesian neural network, radial basis functions, Gaussian processes, k-nearest neighbor, generalized regression neural network, classification and regression trees, and support vector regression. We show the superiority of the Gaussian processes algorithm, compared to other methods, in 14 different real GNSS time series studies. The proposed algorithm can achieve up to 4 millimeters in accuracy, with the average accuracy as 2 centimeters across all time series.