Собственные микросостояния озона: как климат расслаивается по высоте

Оригинал: "Eigen Microstates and Their Evolution of Global Ozone at Different Geopotential Heights"
arXiv:2107.00843 · 2021-07-02 · CC BY 4.0 · 2 мин
Метод собственных микросостояний раскладывает глобальное поле озона на независимые моды, каждая из которых отвечает за свой климатический процесс.

Контекст

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

Методы

Использованы ежедневные данные отношения массы озона к массе воздуха из реанализа ERA5 с 1 января 1979 по 30 июня 2020 года, на 37 уровнях давления от 1000 до 1 гПа, с горизонтальным разрешением 2,5°×2,5°. Каждая точка сетки рассматривалась как агент системы, а её временной ряд после нормировки образовывал микросостояние. Сингулярное разложение ансамбля микросостояний дало собственные микросостояния и их эволюцию. Для связи мод с физическими процессами применяли регрессионную модель LOTUS, которая учитывает индексы КДО, Эль-Ниньо, аэрозольную оптическую толщину, солнечный поток и высоту тропопаузы.

Результаты

Без удаления сезонности первое собственное микросостояние объясняет от 10% до 47% изменчивости в зависимости от высоты: максимум у поверхности и в стратосфере, минимум около 300 гПа. Пространственная структура показывает два кластера: северное и южное полушария с противоположными знаками, а граница между ними смещается к северу с высотой. Фаза годового цикла сдвигается: минимум наступает в феврале у поверхности и в июне на 50 гПа. После удаления сезонности первая мода отражает долгосрочный тренд, показывая снижение озона до середины 1990-х и восстановление после 2000 года на 50 гПа, что связывают с Монреальским протоколом. Третья мода на 100 гПа коррелирует с давлением тропопаузы (r = 0,57). Моды, связанные с Эль-Ниньо, появляются между 400 и 150 гПа, их корреляция с индексом Эль-Ниньо растёт с высотой от 0,24 до 0,47. Сигнал квазидвухлетней осцилляции выделяется на 50, 30 и 10 гПа с корреляциями до 0,67.

Значение

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

Развитие

Метод может быть распространён на другие поля: распределение углекислого газа, океанические течения. Авторы также видят потенциал для обнаружения точек невозврата в климатической системе, таких как коллапс Атлантической меридиональной циркуляции или таяние арктического льда.

Влияние

Результаты полезны для климатологии, атмосферной химии и мониторинга озонового слоя.

Следующие шаги

Применить метод к более высокому разрешению и другим переменным, например, температуре или концентрации парниковых газов.

Ключевые проблемы

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

Вес первой моды озона повторяет профиль температуры атмосферы: обе кривые имеют излом около 300 гПа.

s^k_i(t) = \frac{\delta O^k_i(t)}{\sqrt{\langle [\delta O^k_i(t)]^2 \rangle}}
Отклонение от среднего, делённое на стандартное отклонение
A^k = U^k \cdot \Lambda^k \cdot (V^k)^T = \sum_j \lambda^k_j u^k_j \otimes (v^k_j)^T
Сингулярное разложение матрицы данных
\frac{C^k}{C^k_0} \cdot u^k_j = (\lambda^k_j)^2 u^k_j
Собственные значения дают веса мод

Ключевые числа

  • период данных: 1979–2020, 15188 дней
  • уровней давления: 37 (от 1000 до 1 гПа)
  • разрешение сетки: 2,5°×2,5°, 10368 точек
  • максимальная корреляция с КДО: 0,67 на 30 гПа
  • вес первой моды: от 10% до 47% в зависимости от высоты
Понятия
озоновый слой тропосфера тропосферная двухлетняя осцилляция химия атмосферы спутниковый рекорд UAH южная осцилляция эль-ниньо энсо многомерный индекс ЭНСО муссон спутниковые слои геопотенциальная высота углеродный цикл модель Джина–Тиммермана атмосферное давление весенний барьер предсказуемости водородная атмосфера температура поверхности моря парниковый эффект межгодовая изменчивость метан