Сеть предупреждает об Эль-Ниньо за год

Оригинал: "Using Network Theory and Machine Learning to predict El Niño"
arXiv:1803.10076 · 2018-03-27 · CC BY 4.0 · 4 мин
Гибридная модель из статистики и нейросети предсказывает Эль-Ниньо за 12 месяцев, преодолевая весенний барьер предсказуемости.

Контекст

Эль-Ниньо — это тёплый язык воды в восточной экваториальной части Тихого океана, который возникает примерно раз в четыре года. Он влияет на погоду по всему миру: от засух в Австралии до наводнений в Перу. Рядом с этой работой в архиве лежат исследования, которые тоже пытаются заглянуть в будущее климатической системы: одни используют нейросети для раннего предупреждения о засухе, другие строят гравитационные карты галактик с помощью машинного обучения. Общее направление — поиск скрытых паттернов в больших массивах данных. Данная статья стоит в этом же русле, но фокусируется на конкретной задаче: продлить горизонт предсказания Эль-Ниньо за пределы шести месяцев, где классические модели и ансамбли теряют навык.

Методы

Гибридная модель раскладывает временной ряд индекса NINO3.4 (средняя аномалия температуры поверхности в центральной части Тихого океана) на линейную и нелинейную части. Линейную часть предсказывает классическая модель ARIMA(12,1,0) или ARIMA(12,1,1) — она улавливает инерцию и сезонность. Остаток, то есть нелинейные отклонения, предсказывает нейронная сеть прямого распространения с одним или двумя скрытыми слоями. В сеть подают не только сам индекс, но и три дополнительных атрибута: объём тёплой воды (WWV) в экваториальной полосе, вторую главную компоненту остаточного ветрового напряжения (PC2), связанную с порывами западного ветра, и синусоиду с годовым периодом для учёта сезонного цикла. Для прогнозов на 12 месяцев вместо WWV используют сетевой показатель c2, вычисленный по корреляциям аномалий уровня моря. Сеть строят на пороге корреляции 0.9 со скользящим окном в один год. Обучение идёт на 80% временного ряда, тестирование — на оставшихся 20%. Дополнительно проводят кросс-валидацию с 200 случайными точками старта тестового периода.

Результаты

Гибридная модель показала лучшее качество, чем среднее ансамблевого прогноза CFSv2, на сроках 3–6 месяцев. Например, в январе 2010 года все члены ансамбля CFSv2 переоценили индекс NINO3.4, тогда как гибридная модель дала более точный прогноз. Нормированная среднеквадратичная ошибка (NRMSE) у гибридной модели была ниже, чем у CFSv2, при прогнозе на 4 месяца вперёд (0.016 против 0.021 в среднем по тестовому периоду). На сроке 6 месяцев гибридная модель всё ещё обходила ансамбль, хотя и недооценила Эль-Ниньо 2010 года. Самое важное: прогноз на 12 месяцев с использованием сетевого показателя c2 оказался не хуже, а иногда и лучше прогноза на 6 месяцев. Средняя NRMSE для 12-месячного прогноза составила около 0.018, что сопоставимо с 0.017 для 6-месячного. Это первый случай, когда статистическая модель достигает такого горизонта без потери качества. Кросс-валидация подтвердила устойчивость: при разных разбиениях на обучающую и тестовую выборки пики распределения NRMSE совпадали, хотя разброс рос с увеличением горизонта прогноза.

Значение

Работа показывает, что сетевые метрики, отражающие глобальные корреляции в климатической системе, могут служить ранними индикаторами Эль-Ниньо. В отличие от соседних исследований, которые предсказывают засуху или плазменные неустойчивости, здесь впервые объединены физически мотивированные атрибуты с топологией климатической сети в рамках гибридной модели. Это не просто ещё одна нейросеть: выбор атрибутов опирается на модель Zebiak–Cane, что даёт интерпретируемость. Результат бросает вызов представлению о непреодолимом весеннем барьере предсказуемости — возможно, барьер связан не с хаосом, а с недостатком информации в используемых переменных.

Развитие

Метод можно улучшить, заменив линейную часть ARIMA на динамическую модель, учитывающую волновую динамику океана. Также перспективно использовать нелинейные меры связи, такие как взаимная информация с запаздыванием или перенос энтропии, для отбора атрибутов. Расширение сети за пределы Тихого океана, например включение Индийского океана, может дать ещё более ранние сигналы. Наконец, можно изменить функцию потерь, чтобы модель лучше предсказывала экстремальные события, а не только среднюю ошибку.

Влияние

Результаты повлияют на оперативное прогнозирование Эль-Ниньо в метеорологических центрах, позволяя предупреждать о засухах и наводнениях за год. Метод также применим к другим климатическим модам, таким как Индийский океанский диполь.

Следующие шаги

Проверить модель на независимых данных после 2014 года и на других индексах, например NINO1+2. Исследовать, какие именно корреляции в сети дают предсказательную силу, чтобы связать c2 с конкретными физическими процессами.

Ключевые проблемы

В науках о климате остаётся открытым вопрос о природе весеннего барьера предсказуемости: является ли он следствием хаотической динамики или недостатка наблюдательных данных. Эта работа показывает, что добавление сетевых переменных, улавливающих глобальные корреляции, может ослабить барьер, но не объясняет его механизм. Также неясно, насколько универсален сетевой показатель c2 для разных типов Эль-Ниньо (восточного и центрального).

Сетевой показатель c2 — это просто доля узлов, которые образуют пары. В климатической сети он растёт, когда локальные корреляции усиливаются перед Эль-Ниньо, и падает после. По сути, океан начинает «шептаться» парами, прежде чем закричать.

NRMSE(y_A, y_B) = \frac{1}{\max(y_A, y_B) - \min(y_A, y_B)} \sqrt{\frac{\sum_{k=1}^{n} (y_A^k - y_B^k)^2}{n}}
Ошибка прогноза, делённая на размах наблюдаемых значений; чем меньше, тем лучше.
c_s = \frac{s n_s}{N}
Доля узлов, входящих в кластеры размера s; для s=2 это пары связанных узлов.
Z_t = Y_t + N_t
Наблюдаемый ряд раскладывается на линейную часть Y_t (ARIMA) и нелинейный остаток N_t (нейросеть).

Ключевые числа

  • горизонт прогноза: 12 месяцев
  • NRMSE для 12-месячного прогноза: 0.018
  • NRMSE для 6-месячного прогноза: 0.017
  • порог корреляции для сети: 0.9
  • длина скользящего окна: 1 год
Понятия
весенний барьер предсказуемости южная осцилляция эль-ниньо регрессия нейронной сети энсо модель Джина–Тиммермана индекс Эль-Ниньо ла-нинья многомерный индекс ЭНСО нейронные сети прогностического обучения распространение экваториальных волн модель Зебиака–Кейна ансамблевый прогноз машинное обучение навык прогноза NOAA CPC ERSST NOAA PMEL TAO относительный они буйковые станции TAO/TRITON предсказание