Волны денге: что связывает лихорадку с океаном и листвой

Оригинал: "Wavelet Analysis of Dengue Incidence and its Correlation with Weather and Vegetation Variables in Costa Rica"
arXiv:2107.05740 · 2021-07-02 · CC BY 4.0 · 3 мин
Вейвлет-анализ выявил, что вспышки денге в Коста-Рике пульсируют с периодами 1, 2 и 3 года, синхронно с растительностью, температурой поверхности и океаническими индексами.

Контекст

Денге в Коста-Рике эндемичен с 1993 года, и с тех пор вспышки повторяются ежегодно, но их ритм неодинаков в разных частях страны. Рядом с этой работой лежат исследования, которые смотрят на город как на организм и различают эпидемии не по месту, а по времени. Настоящее исследование встраивается в это русло, но берёт не один город и не одну переменную, а 32 кантона и шесть спутниковых индексов, чтобы увидеть, как денге дышит в масштабе всей страны.

Методы

С января 2000 по декабрь 2019 года собирали еженедельные данные о заболеваемости денге по 32 кантонам, агрегировали их помесячно, извлекли квадратный корень и стандартизировали. Затем для каждого ряда вычислили вейвлет-спектр мощности с вейвлетом Морле, чтобы выделить доминирующие периоды. Для связи с климатом использовали вейвлет-когерентность и фазовые сдвиги между заболеваемостью и индексами: EVI, NDVI, NDWI, дневной и ночной температурой поверхности, TNA и четырьмя индексами Эль-Ниньо. Значимость оценивали по 1000 бутстреп-серий.

Результаты

Доминирующие периоды денге во всех 32 кантонах лежат в полосах 1, 2 и 3 года, причём до 2008 года преобладали циклы 1 и 2 года, а после 2008 года усилился 3-летний цикл. Растительные индексы EVI, NDVI и NDWI значимо коррелируют с заболеваемостью в годовом цикле в кантонах центральной и северной тихоокеанской части, с фазовым сдвигом около трёх месяцев. Температура поверхности коррелирует с денге во всех кантонах, но сигналы противофазны: пик заболеваемости отстаёт от температурного минимума на 3–6 месяцев. TNA связан с денге в периодах 1 и 3 года, особенно после 2011 года. Индексы Эль-Ниньо 3, 3.4 и 4 показывают сильную корреляцию в 3-летнем периоде в большинстве кантонов.

Значение

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

Развитие

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

Влияние

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

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

Ближайший шаг — включить в анализ осадки и влажность, а также проверить устойчивость связей на более длинных рядах после 2019 года.

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

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

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

W_x(s,\tau) = \frac{1}{\sqrt{s}} \int_{-\infty}^{\infty} x(t) \Psi^*\left(\frac{t-\tau}{s}\right) dt
Разложение сигнала по масштабам и временным сдвигам
R_{x,y}(s,\tau) = \frac{|\langle W_{x,y}(s,\tau) \rangle|^2}{|\langle W_x(s,\tau) \rangle|^2 |\langle W_y(s,\tau) \rangle|^2}
Мера синхронности двух временных рядов на каждом масштабе и времени

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

  • кантонов: 32
  • лет наблюдений: 20
  • доминирующие периоды: 1, 2, 3 года
  • фазовый сдвиг растительности: 3 месяца
  • фазовый сдвиг температуры: 3–6 месяцев
Понятия
океанический индекс нино энсо многомерный индекс ЭНСО NOAA CPC Ниньо 1.2 декадная изменчивость относительный они межгодовая изменчивость модель Зебиака–Кейна северное экваториальное противотечение сезонный цикл атлантическая мультидекадная осцилляция