В 1960 году метеоролог на станции в Перу записывает температуру воды, а фермер в Индии смотрит на небо в ожидании муссона. Между этими двумя событиями нет прямой связи, но через несколько месяцев одно отзовётся в другом. Климат Земли похож на разговор, в котором собеседники отвечают не сразу, а с паузой: океан нагревается, и лишь спустя месяцы на другом континенте меняется режим дождей. Старые методы слышали только одновременные реплики, поэтому половина разговора оставалась неразборчивой.

Авторы новой работы решили научить статистику слышать паузы. Каждый временной ряд температуры и осадков они превратили в аналитический сигнал с помощью преобразования Гильберта, которое сдвигает колебания на четверть периода. Это позволило учесть, что максимум тепла в одном регионе может наступить раньше, чем пик дождей в другом. Чтобы избежать искажений на краях шестидесятилетнего ряда, данные продлили в прошлое и будущее с помощью модели Тета, а затем обрезали лишнее. После этого вычислили ковариационную матрицу между двумя комплексными полями и применили сингулярное разложение. Первые 150 мод, объясняющие 99.82% ковариации, подвергли вращению Варимакса, чтобы сделать паттерны более локальными и понятными.

Результат оказался похож на расшифровку многоголосого хора. Метод выделил шесть основных мод связанной изменчивости океана и осадков. Первая мода (62.7% ковариации) описывает сезонный цикл: фазовая функция позволила вычислить, что сезонный максимум температуры в северной части Тихого океана следует за южной с задержкой 196 дней, а в Индийском океане осадки отстают от температуры на 81 день. Вторая мода (7.0%) соответствует Эль-Ниньо: корреляция с индексом ONI составила 0.941. Третья мода (6.9%) отражает глобальное потепление с корреляцией 0.781 с индексом потепления океана. Четвёртая (1.6%) воспроизводит Тихоокеанскую декадную осцилляцию (корреляция 0.565), пятая (1.6%) — Эль-Ниньо Модоки (0.506), шестая (1.1%) — Атлантическую меридиональную моду (0.741).

Раньше эти телесвязи находили по отдельности в разрозненных региональных исследованиях, и каждый раз приходилось заранее предполагать, какая задержка важна. Новый метод впервые объединяет комплексное представление сигнала с вращением Варимакса в рамках анализа максимальной ковариации и одной процедурой извлекает все основные моды вместе с их временными сдвигами. Он не требует предварительного удаления трендов, что критично в условиях меняющегося климата.

Теперь разговор океана и суши можно слушать целиком, со всеми паузами. Фермер в Индии и метеоролог в Перу связаны невидимой нитью задержанного сигнала, и новый метод учит нас различать её среди шума. Впереди — применение к осцилляции Маддена-Джулиана и другим парам переменных, от солёности океана до влажности почвы. Возможно, мы услышим новые реплики в этом древнем диалоге.