Модели климата средней сложности часто содержат петли обратной связи с большими временами запаздывания, связанными с переносом энергии через океан. В случае Эль-Ниньо задержки возникают из-за времени, которое требуется волнам Кельвина и Россби, чтобы пересечь Тихий океан. Традиционные исследования опирались на длительное численное интегрирование, что не позволяет увидеть всю структуру возможных режимов. Данная работа помещает модель ЭНЮК в контекст современной теории динамических систем с запаздыванием, где бифуркационный анализ становится стандартным инструментом.
Авторы используют модель запаздывающего осциллятора, предложенную Циперманом и коллегами, которая сводит всю физику к одному скалярному уравнению с двумя запаздываниями и параметрическим годовым возбуждением. С помощью пакета DDE-Biftool они численно продолжают положения равновесия и периодические орбиты как решения краевых задач, отслеживая их устойчивость через спектр характеристической функции и мультипликаторы Флоке. Резонансные языки вычисляются как поверхности в расширенном пространстве, что позволяет обойти медленную сходимость вблизи порога возбуждения.
Построена полная бифуркационная диаграмма в плоскости параметров (dk, k0), где dk — амплитуда годовой вариации связи океана и атмосферы, а k0 — её среднее значение. Обнаружены все резонансные языки k:ℓ с ℓ ≤ 11, включая основные 1:3 и 1:4, а также кривые бифуркаций тора и удвоения периода. Показано, что тривиальное решение устойчиво в обширной области ниже кривой тора, а выше неё возникают устойчивые резонансные периодические орбиты. Точка параметров, использованная в исходной работе Ципермана, лежит чуть левее языка 1:3, что объясняет наблюдавшуюся там непериодическую динамику, близкую к трёхлетнему циклу.
Работа демонстрирует, что численный бифуркационный анализ дифференциальных уравнений с запаздыванием сегодня вполне осуществим и даёт информацию, недоступную при простом моделировании. В отличие от соседних исследований, которые либо ограничиваются симуляциями, либо рассматривают более простые модели без запаздывания, здесь впервые систематически прослежены резонансные поверхности в полной плоскости параметров. Это позволяет строго определить границы предсказуемости ЭНЮК и чувствительность модели к плохо известным параметрам.
Методику можно распространить на более реалистичные модели ЭНЮК с несколькими запаздываниями и пространственной структурой. Также перспективно исследование перекрытия резонансных языков при больших амплитудах возбуждения, что может приводить к хаотической динамике. Другой путь — включение стохастических возмущений для сравнения с линейными моделями, управляемыми шумом.
Результаты имеют прямое значение для климатологии, поскольку дают строгий инструмент анализа чувствительности моделей ЭНЮК. Они также полезны в других областях, где встречаются дифференциальные уравнения с запаздыванием, например в лазерной физике и популяционной биологии.
Следующим шагом может стать бифуркационный анализ более детальных моделей, таких как модель Зебиака-Кейна, с учётом пространственной изменчивости. Также важно изучить влияние стохастического шума на резонансные языки и сравнить с наблюдаемыми данными.
В климатологии остаётся открытым вопрос о том, является ли Эль-Ниньо детерминированным хаосом или линейной системой, возбуждаемой шумом. Данная работа предоставляет бифуркационную основу для различения этих гипотез, показывая, какие области параметров соответствуют устойчивым периодическим режимам, а какие — квазипериодическим или хаотическим.
Волны Кельвина, переносящие тёплые аномалии через Тихий океан, движутся со скоростью около 2–3 метров в секунду, преодолевая 10 000 километров примерно за два месяца.
В романе Кима Стэнли Робинсона «Сорок знаков дождя» климатические модели с обратными связями и задержками играют ключевую роль в предсказании резких изменений погоды, что перекликается с идеей чувствительности ЭНЮК к параметрам.