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