Uncertainty estimation for time series forecasting via Gaussian process regression surrogates