# Апанас: предварительная проверка тепловой аномалии по Landsat

**Дата анализа:** 12 сентября 2026 года  
**Статус:** воспроизводимый точечный дистанционный анализ; требуется площадное картирование и независимая сверка.

## 1. Метод

Через публичный STAC-каталог Microsoft Planetary Computer отобраны малоблачные сцены Landsat 8/9 Collection 2 Level 2 за январь–март 2015–2026 годов для района Апанаса.

Для каждой доступной сцены сравнивались пиксели:

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

Использован продукт **Surface Temperature, lwir11/ST_B10**. Пересчёт выполнен по формуле USGS:

`ST (K) = DN × 0.00341802 + 149.0`, затем из K вычиталось 273,15.

Пиксели проверялись по QA_PIXEL; облачные и затенённые облаками значения исключались.

Официальное описание и коэффициент USGS: https://www.usgs.gov/landsat-missions/landsat-collection-2-surface-temperature

Описание QA_PIXEL: https://www.usgs.gov/landsat-missions/landsat-collection-2-quality-assessment-bands

Документация Planetary Computer Data API: https://planetarycomputer.microsoft.com/docs/quickstarts/using-the-data-api/

## 2. Предварительные результаты

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

| Дата Landsat | Центр, °C | Север, °C | Юг, °C | Соседняя нарушенная территория, °C | Центр минус нарушенный контроль |
|---|---:|---:|---:|---:|---:|
| 15.02.2016 | −9,56 | −4,29 | −20,42 | −19,91 | **+10,35** |
| 17.03.2018 | +9,24 | +11,52 | −1,17 | +1,45 | **+7,80** |
| 13.03.2020 | +4,85 | +9,51 | −2,43 | −2,14 | **+6,99** |
| 28.02.2021 | −16,35 | −11,35 | −22,56 | −23,11 | **+6,76** |
| 31.01.2022 | −12,58 | −10,31 | −20,54 | −21,02 | **+8,45** |
| 08.02.2022 | −10,55 | −7,95 | −19,73 | −18,81 | **+8,25** |
| 22.03.2023 | +5,45 | +6,72 | −6,56 | −1,21 | **+6,66** |
| 23.03.2023 | +9,72 | +11,57 | −3,62 | +1,75 | **+7,97** |
| 28.01.2024 | −11,10 | −10,11 | −18,10 | −18,00 | **+6,90** |
| 09.03.2024 | +0,32 | +3,63 | −6,79 | −7,09 | **+7,41** |
| 28.03.2025 | +10,76 | +10,76 | −0,81 | +2,22 | **+8,54** |

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

## 3. Что это означает

Повторение аномалии в разные годы, в морозные даты и относительно другого открытого нарушенного участка существенно ослабляет простое объяснение только тёмным цветом породы. Особенно показательно, что:

- северная и центральная части почти всегда наиболее тёплые;
- южная часть в те же минуты значительно холоднее;
- разность между центром и соседней нарушенной территорией удерживается примерно в диапазоне **+6,7…+10,4 °C**.

Это является **сильным дистанционным свидетельством устойчивой локальной поверхностной тепловой аномалии** в центральной/северной части массива по меньшей мере в 2016–2025 годах.

## 4. Чего анализ пока не доказывает

Точечные значения сами по себе ещё не доказывают:

- физическую причину аномалии;
- площадь подземного горения;
- температуру внутренних очагов;
- химический состав газов;
- что каждый пиксель относится именно к федеральному объекту площадью 81,8096 га;
- тождество всего видимого массива одному юридическому объекту.

Альтернативные факторы — различия экспозиции склонов, структуры поверхности, влажности и эмиссивности — требуют проверки. Однако они должны объяснить многолетнюю пространственно устойчивую разницу относительно близких контрольных участков.

## 5. Связь со снегом

QA_PIXEL в ряде сцен классифицирует пиксели массива как снежные, несмотря на существенно более высокую рассчитанную температуру поверхности. Landsat имеет разрешение порядка 30 м для этого продукта; один пиксель может смешивать снег, открытую породу и локально нагретую поверхность. Поэтому сочетание «снег + повышенная температура» нельзя интерпретировать буквально как полностью покрытый снегом горячий участок без просмотра исходных каналов.

Главный результат зимнего анализа — не визуальное отсутствие снега, а **устойчивая разность температур между одновременными соседними пикселями**.

## 6. Файлы воспроизводимости

- `Источники/Landsat_Aпанас/thermal_points.csv`;
- `Источники/Landsat_Aпанас/thermal_points.json`;
- `tools/landsat_thermal_points.mjs`.

## 7. Следующий шаг

1. Построить площадные температурные карты наиболее ясных сцен 15.02.2016, 31.01.2022 и 28.01.2024. **Выполнено.**
2. Проверить, совпадают ли максимумы трёх дат пространственно. **Предварительно выполнено: основные максимумы устойчиво повторяются в северной и центрально-южной частях одного и того же вытянутого массива.**
3. Сравнить аномалию с формой северной и южной частей массива.
4. Наложить приблизительный контур федерального объекта и будущую карту КЕМ 00936 ТЭ.
5. Сверить зоны максимумов с наземными фото/видео дыма, провалов и измерений температуры.

Площадные карты сохранены:

- `Источники/Landsat_Aпанас/Апанас_Landsat_ST_2016-02-15.png`;
- `Источники/Landsat_Aпанас/Апанас_Landsat_ST_2022-01-31.png`;
- `Источники/Landsat_Aпанас/Апанас_Landsat_ST_2024-01-28.png`.

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

## Промежуточный вывод

После спутниковой проверки формулировку можно усилить с «объект визуально похож на горящий» до следующей:

> В центральной и северной части крупного отвального массива у Апанаса Landsat систематически фиксирует локальную температурную аномалию зимой и ранней весной; в проверенных ясных сценах 2016–2025 годов центральный пиксель был примерно на 6,7–10,4 °C теплее соседнего нарушенного горными работами контрольного участка. Причина аномалии должна быть подтверждена площадным анализом и наземными данными, однако результат согласуется с сообщениями о длительном внутреннем горении.