Коммиттор как компас: где заканчивается предсказуемость климата ✎

Оригинал: "Committor Functions for Climate Phenomena at the Predictability Margin: The example of El Niño Southern Oscillation in the Jin and Timmerman model"
arXiv:2106.14990 · 2021-06-28 · arXiv non-exclusive · собственный разбор · 3 мин
Функция коммиттора раскрывает, когда вероятность события Эль-Ниньо предсказуема, а когда — нет.

Контекст

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

Методы

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

Результаты

Главный результат — выделение четырёх областей фазового пространства по предсказуемости события сильного Эль-Ниньо: две области детерминированной предсказуемости (вероятность 0 или 1), область вероятностной предсказуемости (вероятность между 0 и 1, слабо зависящая от начальных условий) и область вероятностной непредсказуемости (сильная чувствительность вероятности к начальным условиям). При промежуточной амплитуде шума (σ = 5×10⁻⁵) эти области сосуществуют, причём вероятностно предсказуемая область возникает из-за локальной неустойчивости, связанной с положительными конечновременными показателями Ляпунова. Кроме того, обнаружено, что среднее время перехода между аттракторами не подчиняется закону Аррениуса, а растёт как степенная функция σ⁻¹ в широком диапазоне малых шумов, что противоречит ожиданиям из теории больших уклонений.

Значение

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

Развитие

Тема может развиваться в направлении применения коммиттора к более реалистичным климатическим моделям и реальным данным. Авторы уже анонсируют использование методов машинного обучения для вычисления коммиттора в модели CESM2.1 для прогноза экстремальных тепловых волн. Также перспективно исследование математической природы нарушения закона Аррениуса, возможно, связанного с фрактальной структурой границ бассейнов притяжения.

Влияние

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

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

Следующие шаги включают разработку эффективных алгоритмов на основе машинного обучения для вычисления коммиттора в высокоразмерных системах и теоретическое объяснение наблюдаемого нарушения закона Аррениуса.

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

Работа напрямую связана с фундаментальной проблемой предсказуемости хаотических систем: как долго и с какой точностью можно прогнозировать редкие события при наличии шума. Она также затрагивает вопрос об универсальности закона Аррениуса для переходов между сложными аттракторами, который остаётся открытым в теории динамических систем.

Среднее время перехода между режимами в модели может достигать 500 000 безразмерных единиц времени, что соответствует тысячелетиям, при этом распределение времён перехода оказывается экспоненциальным, как у радиоактивного распада.

q(x) = P( max_{0 ≤ t ≤ T} x(t) > -1 | X(0) = x )
Вероятность того, что траектория, стартующая из точки x, превысит порог -1 (сильное Эль-Ниньо) в течение времени T.
E[τ_c] ~ A exp(ΔV/σ²)
Ожидаемая экспоненциальная зависимость среднего времени перехода от обратного квадрата амплитуды шума, которая нарушается в модели.

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

  • период предельного цикла: 186 безразмерных единиц (около 50 лет)
  • среднее время перехода при σ=5×10⁻⁵: 1039 безразмерных единиц
  • амплитуда шума для промежуточного режима: 5×10⁻⁵
  • диапазон времён перехода: от 25 до 3000 лет
Понятия
эль-ниньо модель Джина–Тиммермана вероятностный прогноз весенний барьер предсказуемости ла-нинья межгодовая изменчивость прогноз сезонный цикл южная осцилляция эль-ниньо энсо NOAA CPC предсказуемость эффект конденсатора Индийского океана осадки многомерный индекс ЭНСО предсказание сезонный форсинг начальные условия аномалия осадков IRI