Два шага к точному прогнозу ливней ✎

Оригинал: "A two-step approach to model precipitation extremes in California based on max-stable and marginal point processes"
arXiv:1204.0286 · 2012-04-02 · arXiv non-exclusive · собственный разбор · 2 мин
Как ежедневные записи метеостанций превращаются в карту риска редких наводнений.

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

Первый шаг похож на настройку каждого микрофона в студии перед записью оркестра. Для каждой из 35 метеостанций Калифорнии по отдельности подбирают кривую, описывающую хвост распределения осадков. В хвосте живут редкие, но разрушительные ливни. Здесь важны не только средние дожди, но и то, насколько быстро убывает вероятность экстремальных значений. Настройка идёт по ежедневным превышениям порога, а не по годовым пикам, поэтому кривая опирается на сотни точек вместо 55.

Ключевая идея: сначала очистить каждый сигнал по отдельности, а уж потом настраивать связи между ними. Так шум одного датчика не портит настройку соседнего.

Второй шаг наводит резкость на связи между станциями. Дождь в Напе и дождь в Уинтерсе часто приходят вместе, потому что их рождает один атмосферный фронт. Чтобы измерить эту связь, метод берёт уже не ежедневные записи, а годовые максимумы на парах станций. Попарное сравнение похоже на прослушивание дуэтов: если два инструмента всегда звучат громко одновременно, значит, между ними есть общая причина. Так оценивают, насколько синхронны экстремумы в разных точках карты.

Результат виден на карте рисков. Традиционный подход, опиравшийся только на годовые максимумы, находил значимое влияние Эль-Ниньо на 14 станциях. Новый метод увидел этот эффект на 22 станциях. Доверительные интервалы для совместных 50-летних уровней возврата сузились на 27–51%. Это значит, что страховщик или городской планировщик получает не размытое пятно, а чёткую линию, где заканчивается зона приемлемого риска.

Оценка эффективности традиционного метода составила от 23% до 95% в зависимости от параметра. Для формы хвоста, самого капризного параметра, эффективность падала до 45%.

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

Совместный 50-летний уровень возврата для пары станций Напа и Уинтерс при значении индекса Южного колебания -1 составил около 10–15 см осадков за сутки, в зависимости от метода оценки.

G(z) = \exp\left\{-\left[1 + \xi\left(\frac{z-\mu}{\sigma}\right)\right]^{-1/\xi}\right\}
Распределение максимумов блоков, где μ — параметр положения, σ — масштаба, ξ — формы.
\Pr\{Z(x_1) \le z, \ldots, Z(x_p) \le z\} = \exp(-\zeta/z)
Мера зависимости экстремумов: ζ=1 при полной зависимости, ζ=p при независимости.
Понятия
многомерный индекс ЭНСО NOAA CPC максимальное правдоподобие океанический индекс нино ERSST Ниньо 1.2 вероятностный прогноз энсо аномалия осадков относительный они ансамблевый прогноз байесовское свидетельство модель Зебиака–Кейна критерий отношения правдоподобия весенний барьер предсказуемости вероятностная модель вложенная выборка методы Монте-Карло аналоговый прогноз оценка параметров