Лотерейный билет климата: почему Эль-Ниньо греет втрое сильнее, чем кажется ✎

Оригинал: "Direct and Indirect Effects -- An Information Theoretic Perspective"
arXiv:1912.10508 · 2019-12-22 · arXiv non-exclusive · собственный разбор · 2 мин
Новые меры причинного влияния, привязанные к конкретным значениям, показывают: эффект Эль-Ниньо на температуру Тихоокеанского Северо-Запада втрое больше среднего.

Выигрыш в лотерею меняет жизнь победителя, но для статистики всех игроков почти ничего не значит: вероятность выигрыша ничтожна, и средний выигрыш на билет остаётся крошечным. В климате есть похожие события: фаза Эль-Ниньо случается не каждый год, зато когда случается, температура Тихоокеанского Северо-Запада меняется заметно. Классические меры причинного эффекта усредняют влияние по всем возможным фазам, поэтому редкий, но мощный эффект Эль-Ниньо тонет в среднем. Авторы работы предложили способ измерять влияние именно конкретного значения причины, как если бы статистик смотрел не на средний выигрыш всех билетов, а на выигрыш именно этого, вытянутого из лотерейного барабана.

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

Чтобы отделить эффект конкретного значения от среднего, авторы сравнивают два распределения вероятностей исхода. Первое — распределение температуры после того, как мы принудительно задали причину, например фазу Эль-Ниньо. Второе — «естественное» распределение, усреднённое по всем возможным фазам. Различие между этими распределениями измеряется дивергенцией Кульбака-Лейблера, то есть мерой того, насколько одно распределение отличается от другого. Если после вмешательства температура сдвигается сильно, дивергенция велика, и эффект конкретного значения оказывается большим.

На данных климатической модели CESM2 за 1850–2015 годы эффект Эль-Ниньо на температуру Тихоокеанского Северо-Запада оказался примерно втрое больше среднего эффекта по всем фазам. Прямой эффект, который передаётся через атмосферную циркуляцию, превышает косвенный, идущий через персистентность температуры. Это согласуется с физическим механизмом: Эль-Ниньо влияет на регион напрямую, а не через память океана.

Самое сильное влияние обнаружено для условного эффекта прошлой температуры на текущую при фазе Эль-Ниньо: вмешательство, задающее холодную аномалию в год Эль-Ниньо, даёт эффект примерно в пять раз сильнее среднего.

Такой подход сближает два мира: статистический причинный вывод, который умеет оценивать эффекты вмешательств, и теорию информации, которая улавливает нелинейные взаимодействия. Раньше информационные меры не различали вклад конкретных значений причины; теперь они отвечают на вопрос «каков эффект именно этого значения?». Метод можно применить не только к климату, но и к нейронауке, где важно оценить влияние конкретного стимула на активность нейронов, или к экономике, где интересует эффект определённого политического решения. Выигрышный билет наконец-то можно измерить отдельно от всех остальных.

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

STE(x \to Y) = D(p(Y | \hat{x}) \| \sum_{x'} p(x') p(Y | \hat{x'}))
Дивергенция Кульбака-Лейблера между распределением Y после принудительного задания X=x и усреднённым по всем возможным x распределением Y.
SCDE(x \to Y; z) = D(p(Y | \hat{x}, \hat{z}) \| \sum_{x'} p(x') p(Y | \hat{x'}, \hat{z}))
Эффект X=x на Y при фиксированном медиаторе Z=z.
SNIE(x \to Y) = D(p(Y | \hat{x}) \| \sum_{x',z'} p(x') p(z' | \hat{x'}) p(Y | \hat{x}, z'))
Эффект X=x на Y, передающийся только через медиатор Z.
Понятия
IRI причинность по Грейнджеру индекс Эль-Ниньо информационный критерий эль-ниньо относительный они причинность эпидемия индекс южной осцилляции эффект конденсатора Индийского океана межгодовая изменчивость сезонный форсинг NOAA CPC телеконнекция теория информации сезонный цикл прогноз североатлантическая осцилляция линейная инверсная модель ERSST