Дирижёр оркестра постепенно теряет контроль над темпом. Сначала вторая скрипка чуть отстаёт, потом виолончели ускоряются, и вот уже каждый музыкант играет в своём ритме, хотя ещё минуту назад все следовали одной мелодии. Примерно так в климатической системе Эль-Ниньо две независимые частоты — собственный ритм обратной связи и навязанный сезонный цикл — могут сосуществовать годами, а затем резко рассогласоваться.
Внутри этой кухни работает модель ГЗТ: простое уравнение с запаздыванием, где скорость изменения глубины термоклина зависит от её значения в прошлом. Когда параметры попадают в резонансный язык 2:7 (два цикла обратной связи за семь лет), в фазовом пространстве возникают торы — поверхности, по которым траектория движется, как нить по бублику, наматывая витки с двумя разными периодами.
Теоретики давно предсказали, что при изменении параметров два таких тора могут сблизиться и исчезнуть, оставив после себя область сложной динамики — пузырь Шенсинера. Но одно дело нарисовать схему на бумаге, и совсем другое — увидеть, как она разворачивается в уравнении с бесконечномерным фазовым пространством. Авторы работы использовали программный пакет DDE-Biftool, чтобы продолжать периодические решения и вычислять их мультипликаторы Флоке — числа, показывающие, насколько устойчив каждый цикл.
Внутри пузыря для резонанса 2:7 обнаружилась структура, хорошо согласующаяся с теоретической картиной Баесенса и Маккея: кривые нейтральных сёдел, гомоклинических и гетероклинических переходов, а также кривые складывания торов. Но нашлось и отличие: бифуркация тора оказалась субкритической, то есть рождался седловой тор, а не устойчивый. Это значит, что при пересечении критической линии система не плавно переходит на новый режим, а срывается в резкое переключение, как если бы оркестр не постепенно расстроился, а разом перешёл на какофонию.
Раньше было неясно, как именно разрушаются торы в дифференциальных уравнениях с запаздыванием, и совпадает ли бифуркационная структура с теоретической схемой для отображений. Теперь впервые для запаздывающего уравнения численно воспроизведена полная картина, и показано, что бифуркация тора субкритическая. Это стало возможным благодаря продолжению периодических решений и вычислению неустойчивых многообразий седловых орбит в двумерных проекциях, связанных с сечениями Пуанкаре. Такой подход позволил отследить все ключевые бифуркации внутри резонансного языка 2:7.
Результаты могут помочь в интерпретации климатических временных рядов и в поиске предвестников резких переходов в системах с несколькими частотами, хотя прямого практического применения пока нет. Дирижёр, потерявший контроль, — это не только красивая метафора, но и предупреждение: если мы научимся распознавать момент, когда торы готовы сложиться, то, возможно, успеем расслышать первые фальшивые ноты в климатическом оркестре.
Резонансный язык 2:7 в модели означает, что за 7 лет совершается ровно 2 цикла обратной связи, и именно этот язык оказался самым крупным в рассматриваемой области параметров, что облегчило численный анализ.
Идея резких климатических переходов из-за взаимодействия нескольких частот напоминает сюжет фильма «Послезавтра», где нарушение океанической циркуляции приводит к катастрофическим изменениям, хотя в реальности такие переходы могут быть более постепенными и связаны с квазипериодической динамикой.