Эль-Ниньо — Южное колебание (ЭНСО) — это нерегулярное чередование тёплых (Эль-Ниньо) и холодных (Ла-Нинья) фаз температуры поверхности океана в тропической части Тихого океана. Оно повторяется каждые 3–7 лет и сбивает привычный ход ветров и осадков по всей планете. Существует целая иерархия моделей ЭНСО: от простейших концептуальных «коробочных» моделей с несколькими уравнениями до сложных моделей общей циркуляции. Концептуальные модели ценны своей математической прозрачностью и малыми вычислительными затратами, но обычно строятся на физической интуиции, а не выводятся строго из более полных систем. Данная работа заполняет этот пробел: она показывает, как из пространственно-распределённой стохастической динамической системы с линейным устойчивым ядром и мультипликативным шумом можно строго получить низкоразмерные концептуальные модели.
Исходная пространственно-распределённая система описывает взаимодействие атмосферы, океана и температуры поверхности океана (SST) на межгодовом масштабе. Её динамическое ядро линейно, детерминированно и устойчиво, а внутрисезонные порывы ветра (включая влияние осцилляции Маддена–Джулиана) параметризованы простым стохастическим процессом с мультипликативным шумом. Амплитуда шума зависит от SST в западной части Тихого океана через двухуровневый марковский скачкообразный процесс. Для понижения размерности применяется анализ главных компонент на основе разложения по собственным векторам матрицы линеаризованной системы. Выделяются две ведущие комплексно-сопряжённые моды с периодом около 4.6 года и временем затухания 1.5 года, которые объясняют подавляющую часть изменчивости. Проецируя стохастическое воздействие на эти моды, получают двумерную систему обыкновенных дифференциальных уравнений.
Проекция полного решения на две ведущие моды почти неотличима от исходного решения по пространственно-временным паттернам и ключевым статистикам: спектр мощности имеет пик в диапазоне 3–7 лет, а функция плотности вероятности SST в восточной части Тихого океана сильно не-гауссова с положительной асимметрией и тяжёлым хвостом, соответствующим экстремальным Эль-Ниньо. Корреляция между реконструированной и наблюдаемой SST в восточной части Тихого океана превышает 0.85 на всех долготах. В рамках общей схемы строго выводятся стохастические версии разрядно-зарядного осциллятора и запаздывающего осциллятора. Для разрядно-зарядного осциллятора коэффициенты связи положительны для положительной обратной связи Бьеркнеса и механизма подзарядки, а отрицательный коэффициент описывает разрядку. Фазовый сдвиг между SST востока и термоклином запада составляет 122 градуса, а не 90, что согласуется с наблюдениями. В запаздывающем осцилляторе появляется экспоненциальное затухание запаздывающего члена и два типа вклада шума: мгновенный и с памятью.
Работа стоит особняком среди соседних исследований: она не просто предлагает очередную концептуальную модель, а даёт строгий мост между пространственно-распределёнными и низкоразмерными описаниями ЭНСО. В отличие от многих существующих моделей, которые фокусируются на малом наборе переменных или осреднённых величинах, здесь низкоразмерная модель сохраняет полные пространственно-временные паттерны для атмосферы, океана и SST через собственные векторы. Это позволяет напрямую связывать различные компоненты системы и проверять гипотезы о причинно-следственных связях. Строгий вывод также проясняет роль мгновенных и запаздывающих эффектов стохастического шума в крупномасштабной динамике ЭНСО.
Текущая схема ориентирована на восточно-тихоокеанские события Эль-Ниньо. Для учёта разнообразия ЭНСО, включая центрально-тихоокеанские Эль-Ниньо, необходимо добавить механизм зональной адвекции океана и, возможно, нелинейность в исходную пространственно-распределённую модель. Полученные концептуальные модели могут служить полигоном для тестирования стратегий включения стохастических параметризаций в более сложные модели промежуточной сложности и модели общей циркуляции. Также возможно построение иерархии моделей с разной сложностью с помощью методов, управляемых данными.
Результаты могут улучшить понимание механизмов ЭНСО и предоставить простые, но физически обоснованные модели для прогнозирования индексов Эль-Ниньо и изучения предсказуемости с точки зрения ансамблевых методов.
Включить механизм центрально-тихоокеанских событий и нелинейность в исходную модель, а затем проверить, как низкоразмерные модели воспроизводят разнообразие ЭНСО. Использовать концептуальные модели как тестовые стенды для разработки стохастических параметризаций в более сложных климатических моделях.
Работа напрямую связана с проблемой разнообразия Эль-Ниньо: почему события возникают в разных районах Тихого океана и имеют разную амплитуду и длительность. Строгий вывод низкоразмерных моделей даёт основу для изучения того, как нелинейность и адвекция порождают центрально-тихоокеанские события, которые пока не воспроизводятся линейной схемой.
Фазовый сдвиг между температурой поверхности океана на востоке и глубиной термоклина на западе составляет 122 градуса, а не 90, как можно было бы ожидать для простого гармонического осциллятора. Эта деталь важна для правильного воспроизведения временной задержки в запаздывающем осцилляторе.