Океан поёт на разные голоса ✎

Оригинал: "Rotated spectral principal component analysis (rsPCA) for identifying dynamical modes of variability in climate systems"
arXiv:2004.11411 · 2020-04-23 · arXiv non-exclusive · собственный разбор · 2 мин
Как учёные научились разделять похожие колебания климата, поворачивая оси анализа.

Климат Земли похож на хор, в котором каждый участник поёт свою партию. Одни колебания медленные, как многолетние изменения температуры океана, другие быстрые, как короткие атмосферные волны. Учёным нужно расслышать каждый голос, но голоса перекрываются, и ноты сливаются в кашу.

Спектральный анализ главных компонент — это способ рассортировать голоса по частотам, но и он путается, когда два певца берут одну ноту с одинаковой силой.

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

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

Затем метод применили к реальным данным температуры поверхности Тихого океана за 1940–2018 годы. На низких частотах он выделил моду Эль-Ниньо с периодом 2–7 лет, и её совпадение с известным индексом составило 0.72. На более высоких частотах вращение помогло разделить несколько мод, включая волну в южной части Тихого океана, которая движется на запад со скоростью около 0.4 метра в секунду. Для сравнения: это скорость неторопливой прогулки человека. Раньше такие моды сливались в неразличимый гул, а теперь их можно слушать по отдельности.

Вейвлет Морле назван в честь французского геофизика Жана Морле, который разработал его в 1980-х для анализа сейсмических сигналов, а не климатических данных.

Понятия
анализ главных компонент NOAA CPC нелинейный спектральный анализ Лапласиана ПЦР тихоокеанская декадная осцилляция модель осциллятора с подзарядкой энсо резонансная частота Ниньо 4 сезонный цикл естественная изменчивость спектроскопия с подавлением шума реанализ угловая частота южная осцилляция спектроскопия многомерный индекс ЭНСО спектральная картина спектральный анализ спектроскопия атмосфер экзопланет