Как из хаоса рождается порядок: строгий вывод моделей Эль-Ниньо ✎

Оригинал: "Rigorous Derivation of Stochastic Conceptual Models for the El Niño-Southern Oscillation from a Spatially-Extended Dynamical System"
· Nan Chen, Yinling Zhang
arXiv:2212.14179 · 2022-12-29 · arXiv non-exclusive · собственный разбор · 2 мин
Учёные математически строго вывели простые модели Эль-Ниньо из сложной пространственной системы.

Каждые несколько лет тёплые воды Тихого океана поднимаются к поверхности у берегов Перу, и погода на планете сходит с привычного расписания. Это Эль-Ниньо — нерегулярное явление, которое повторяется раз в 3–7 лет, но до сих пор ускользает от точного прогноза. Учёные описывают его с помощью моделей разной сложности: от простых «коробочных» уравнений до гигантских симуляций всей планеты. Простые модели удобны, но обычно их строят по наитию, а не выводят из более полных систем. Новая работа заполняет этот пробел: она показывает, как из сложного пространственного хаоса строго получить ту самую простую мелодию.

Климатическая система напоминает оркестр, где каждый инструмент играет свою партию. Ветры, течения и температура океана связаны между собой, и чтобы услышать главную тему, нужно выделить её из общего шума.

Исследователи взяли пространственную модель, где атмосфера и океан взаимодействуют на межгодовом масштабе. Её основа линейна и устойчива, а случайные порывы ветра, которые приходят из внутрисезонных процессов, добавляют шум, сила которого зависит от температуры океана на западе Тихого океана. Чтобы упростить систему, они применили разложение по собственным векторам — это как найти особые направления, вдоль которых колебания происходят независимо. Две ведущие моды с периодом около 4.6 года и временем затухания 1.5 года объясняют почти всю изменчивость. Проецируя шум на эти моды, получают два уравнения, которые и есть концептуальная модель.

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

Результат удивляет точностью: реконструкция по двум модам почти неотличима от исходного решения. Спектр мощности показывает пик в нужном диапазоне 3–7 лет, а распределение температуры поверхности океана на востоке имеет положительную асимметрию и тяжёлый хвост, что означает редкие, но сильные Эль-Ниньо. Корреляция с наблюдениями превышает 0.85. В рамках общей схемы строго выводятся известные модели: разрядно-зарядный осциллятор, где тепло накапливается и сбрасывается, и запаздывающий осциллятор, где текущее изменение зависит от прошлого состояния. Оказалось, что фазовый сдвиг между температурой востока и глубиной термоклина на западе составляет 122 градуса, а не 90, как у простого маятника. Эта деталь важна для правильной задержки.

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

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

\frac{d\tilde{a}}{dt} = -d_o \tilde{a} + \tilde{c}_{11} \tilde{a} + \tilde{c}_{12} \tilde{b} + 2g_{a,Re}, \quad \frac{d\tilde{b}}{dt} = -d_o \tilde{b} + \tilde{c}_{21} \tilde{a} + \tilde{c}_{22} \tilde{b} + 2g_{b,Re}
Общий вид низкоразмерной концептуальной модели после проецирования на две ведущие моды; коэффициенты выражаются через собственные значения и собственные векторы.
\frac{dT_E(t)}{dt} = -d_o T_E(t) + \tilde{c}_{11} T_E(t) + \tilde{c}_{12} c_y e^{-d_o \phi / \omega_o} T_E(t - \phi/\omega_o) + \alpha_{T_E} a_p(t) + \text{интегральный член}
Стохастическая версия запаздывающего осциллятора с экспоненциальным затуханием запаздывающего члена и вкладом шума с памятью.
Понятия
энсо южная осцилляция эль-ниньо модель Джина–Тиммермана индекс Эль-Ниньо многомерный индекс ЭНСО модель осциллятора с подзарядкой весенний барьер предсказуемости ла-нинья южная осцилляция буйковые станции TAO/TRITON ERSST естественная изменчивость стохастическая динамическая система объём тёплой воды вспышка западных ветров относительный они NOAA CPC эффект конденсатора Индийского океана северное экваториальное противотечение распространение экваториальных волн