Память океана гасит хаос Эль-Ниньо ✎

Оригинал: "Chaotic dynamics of fractional Vallis system for El-Nino"
arXiv:1808.08075 · 2018-08-24 · arXiv non-exclusive · собственный разбор · 3 мин
В дробной модели Валлиса хаос Эль-Ниньо исчезает при уменьшении порядка дробной производной, а критический порядок вычислен точно.

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

В новой работе уравнения Валлиса переписали с дробными производными. Обычная производная показывает, как быстро меняется величина прямо сейчас. Дробная производная добавляет к этому память: она учитывает, что происходило с системой в прошлом, причём недавние события важнее давних. Это похоже на то, как человек принимает решение, опираясь не только на текущий момент, но и на весь свой опыт. Океан тоже помнит: тёплая вода, пришедшая в прошлом году, продолжает влиять на течения и ветры.

Дробный порядок α в уравнениях работает как регулятор громкости хаоса. При α=1 память отсутствует, и система ведёт себя как классическая модель Валлиса. Чем меньше α, тем сильнее океан помнит прошлое и тем спокойнее становится его поведение.

Исследователи численно решили дробную систему для разных значений параметра p, который описывает среднее влияние экваториальных ветров на океан. При α=1 хаос наблюдается для |p|≤2.45. Когда α уменьшается до 0.985, диапазон хаоса сужается до |p|≤1.5. При α=0.975 хаос остаётся лишь в узкой полосе вблизи p=0. А при α=0.965 хаос исчезает полностью. Это значит, что достаточно сильная память океана способна погасить хаотические колебания Эль-Ниньо.

Критический порядок, ниже которого хаос пропадает, вычислен точно. Для одного набора параметров он равен 0.975, то есть суммарный порядок трёх уравнений составляет 2.925. Для несоразмерной системы, где разные уравнения имеют разные порядки производных, критический суммарный порядок оказался ещё ниже: 2.71 при α1=0.71, α2=α3=1. Это показывает, что память в одном из уравнений может быть решающей.

Наибольший показатель Ляпунова при α=0.98 и p=0.35 равен 0.4181. Положительное значение подтверждает хаос, но при дальнейшем уменьшении α этот показатель становится отрицательным, и траектории перестают разбегаться.

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

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

Модель Валлиса получается из системы Лоренца добавлением всего одного параметра p, который описывает среднее влияние экваториальных ветров на океан. При p=0 и B=1 она в точности совпадает с системой Лоренца.

D^{\alpha_1} x = By - C(x+p), \quad D^{\alpha_2} y = xz - y, \quad D^{\alpha_3} z = -xy - z + 1
Три уравнения описывают эволюцию тока x, разности температур y и средней температуры z с дробными производными порядка α_i.
|\arg(\lambda)| > \frac{\alpha \pi}{2}
Точка равновесия асимптотически устойчива, если аргументы всех собственных значений матрицы Якоби превышают απ/2.
Понятия
теория хаоса сезонный цикл межгодовая изменчивость смешанные моды колебаний индекс Эль-Ниньо случайная флуктуация показатель ляпунова