Спектральный анализ временного ряда среднегодовых значений аномалии глобально осредненной приповерхностной температуры показывает наличие гармоник лунного нодального цикла с периодом 18.6 года, что может использоваться при прогнозе значений ряда. Рассматривались три модели ряда: авторегрессия AR(p), комбинированная модель авторегрессии – проинтегрированного скользящего среднего ARIMA(p,d,q) и искусственная нейронная сеть. Показано, что наилучшие результаты прогноза аномалии глобальной температуры на три года дает модель ARIMA(4,1,4).