Гляциологическая карта Девисона (Davisson map) — это один из самых точных методов оценки ледникового покрова и его изменений во времени. Она позволяет визуализировать среднее распределение толщины льда и выявлять аномальные участки, где таяние или накопление ледяной массы происходит быстрее или медленнее, чем в среднем по региону.
Что такое «средняя карта Девисона»?
В отличие от традиционных гляциологических съёмок, карта Девисона строится не на точечных измерениях, а на усреднённых данных спутниковой лазерной альтиметрии и радиолокационной интерферометрии за длительный период (обычно 5–10 лет). Это позволяет «сгладить» сезонные колебания и получить репрезентативную картину.
С чего начать: подготовка исходных данных
Источники данных
- NASA ICESat-2 ATLAS — лазерная альтиметрия с точностью до 2 см по высоте поверхности льда.
- Sentinel-1 / RADARSAT-2 — синтетическая апертура радара (SAR) для измерения скорости движения ледника.
- MODIS / Landsat 8-9 — мультиспектральные снимки для контроля поверхности (облака, трещины, озёра талой воды).
Временной охват
Для корректного усреднения минимально требуется 5 полных гидрологических лет (1 октября по 30 сентября следующего года). Такой интервал учитывает два полных цикла «зима—лето» и позволяет устранить годовые выбросы.
Пошаговый алгоритм построения
Шаг 1. Создание цифровой модели рельефа (DEM)
- Загрузите исходные тайлы высот из проекта ArcticDEM или REMA.
- Выполните коррекцию смещения относительно WGS-84 и приведите к единой системе координат EPSG:3413.
- Удалите выбросы, вызванные облачностью или тенью, методом MAD-медианного фильтра.
Шаг 2. Нормализация высотных срезов
Для каждого спутникового прохода строим:
- «Снежную поверхность» — усреднённое значение высоты за ±3 дня.
- «Базовую поверхность» — высоту твёрдого грунта под льдом, полученную из радиолокационного зондирования.
Разность этих двух слоёв даёт мгновенную толщину льда.
Шаг 3. Усреднение и интерполяция
- Сложите все мгновенные толщины и разделите на количество проходов.
- Используйте кривые Безье для плавной интерполяции между траекториями спутников.
- Заполните пробелы методом kriging 3-D с учётом анизотропии (вдоль потока ледника).
Шаг 4. Оптимизация цветовой шкалы
| Диапазон толщины, м | Цвет (HEX) | Интерпретация |
|---|---|---|
| 0 – 50 | #ffcccc | Морской прибрежный лёд, высокий риск таяния |
| 50 – 250 | #b3d9ff | Типовой континентальный ледник |
| 250 – 800 | #0073e6 | Колоссальные запасы, устойчивые к климатическим колебаниям |
| 800+ | #003366 | Основные ледяные щиты Гренландии или Антарктиды |
Расшифровка полученной карты
Основные метрики
- Mean Ice Thickness (MIT) — усреднённая толщина льда по всей площади.
- Ice Volume Anomaly (IVA) — отклонение текущего объёма от долгосрочного среднего.
- Balance Velocity Field (BVF) — скорость перемещения льда, рассчитанная по направлениям течения и влиянию рельефа.
Как читать «тёплые» и «холодные» пятна
Тёплые пятна (отрицательная аномалия) — участки, где лёд стал тоньше на ≥25 м за 5-летний период. Чаще всего совпадают с:
- Зонами поверхностного таяния (озёра талой воды).
- Участками повышенной тектонической активности (гейзеры, теплые источники).
Холодные пятна (положительная аномалия) — участки, где лёд утолщился на ≥25 м. Причины:
- Увеличение годового снежного выпадения (сильные снежные бураны).
- Закупорка ледяной дамбой оттока талой воды, что приводит к «запечатыванию» ледяной массы.
Как использовать карту Девисона на практике
Мониторинг изменения климата
Регулярное обновление средней карты (раз в 2–3 года) позволяет отслеживать ускорение таяния и оценивать вклад ледников в повышение уровня Мирового океана.
Планирование научных экспедиций
Области с высокой положительной аномалией выбирают для установки ледовых ядер, там сохраняется самый древний и чистый лёд для анализа изотопов.
Моделирование ледоходов и навигация судов
Карта интегрируется в навигационные системы ледоколов и позволяет прокладывать маршруты, минимизирующие риск застрять в многолетнем льду.
Краткий FAQ
Можно ли строить карту без спутниковых данных?
В принципе да, но для этого потребуется сеть наземных радарных станций и ежегодное аэрогеофизическое зондирование. Это в 8–10 раз дороже и доступно лишь крупным научным консорциумам.
Какой минимальный размер карты?
Рекомендуемый размер тайла, не менее 50×50 км. Меньшие области дают слишком большую статистическую погрешность (>±15 м).
Построение средней карты Девисона — это сочетание точной науки и современных технологий. При правильном подходе вы получаете беспрецедентно точный снимок ледяной толщи, который можно использовать как для академических исследований, так и для практической логистики в Арктике и Антарктике.
Для автоматизации процесса рекомендуется использовать Python-библиотеки icepyx и xDEM, а для визуализации — QGIS 3.34 LTR с плагином Copernicus Glacier Tools.