Работа стоит на отшибе карты исследований: ближайшие соседи касаются сверхпроводимости, квантовых парадоксов и нанорезонаторов, и лишь одна работа затрагивает куперовские пары и сверхтекучесть. Вокруг уже изучены квантовые измерения и численные методы для уравнения Дирака, но данное исследование переносит спектральный анализ марковских операторов в климатическую модель. Это первый случай применения редукции Рюэля-Полликотта к высокоразмерной стохастической системе океана и атмосферы.
Модель Кейна-Зебиака размерностью 930 переменных возмущали винеровским процессом в поле ветрового стресса с уровнем шума 0.01. Наблюдения проецировали на плоскость из восточной температуры поверхности океана и западной глубины термоклина, разбитую на 100×100 ячеек. По 6000-летним симуляциям с шагом 5 дней оценивали матрицы переходов с лагом 4 года, а собственные значения пересчитывали в резонансы через логарифм отношения матриц при лагах 4 года и 4 года плюс полмесяца.
До бифуркации (параметр связи ξ=2.80 и 2.85) резонансы образуют треугольную решётку, порождаемую собственными значениями якобиана стационарной точки. После бифуркации (ξ=2.90 и 2.95) они выстраиваются в параболы, где мнимая часть кратна основной частоте, а вещественная растёт квадратично с номером гармоники. Спектральная щель между ведущими резонансами и мнимой осью сужается при подходе к бифуркации, что обостряет пики в спектре мощности и замедляет спад корреляций. Собственные векторы имеют вид произведения гармоники на полиномы Лагерра или Эрмита, а их фазовые изолинии наклонены из-за ненулевого коэффициента закрутки β.
В отличие от соседних работ, которые анализируют низкоразмерные или квантовые системы, здесь редукция применена к климатической модели промежуточной сложности с тысячами степеней свободы. Это первая проверка того, что параболические структуры резонансов, предсказанные для нормальной формы Хопфа, возникают в реалистичной модели Эль-Ниньо. Работа показывает, что фазовая диффузия, а не только устойчивость предельного цикла, управляет потерей предсказуемости.
Метод можно распространить на неавтономные системы, чтобы учесть сезонный цикл, который фазирует события Эль-Ниньо к декабрю-январю. Также перспективно вычислить коэффициенты фазовой диффузии напрямую из анализа Флоке предельного цикла модели Кейна-Зебиака, что потребует численного продолжения цикла и фундаментальной матрицы.
Результаты повлияют на диагностику предсказуемости Эль-Ниньо в климатических моделях и на теорию стохастических бифуркаций в высокоразмерных системах.
Следующий шаг — распространить редукцию на модели с сезонным форсингом и сравнить коэффициенты фазовой диффузии с наблюдаемыми индексами Эль-Ниньо.
Для климатической динамики остаётся открытым вопрос, насколько фазовая диффузия, вызванная быстрыми атмосферными процессами, ограничивает практическую предсказуемость Эль-Ниньо за пределами нескольких месяцев. Работа связывает эту проблему с измеримым спектральным зазором, но не решает её для полной неавтономной системы с сезонным циклом.
Шумовое воздействие с уровнем всего 0.01 (безразмерная амплитуда ветрового стресса) способно возбудить колебания с периодом 3–4 года в модели, которая без шума просто затухает к стационарному состоянию.