Based on the classical Gaussian process (GP) model, we propose a multi-scale Gaussian process (MGP) model to predict the existence of chaotic time series. The MGP employs a covariance function that is constructed by a scaling function with its different dilations and translations, ensuring that the optimal hyperparameter is easy to determine.