Введение
Согласно федеральному законодательству [1] основной целью системы государственного экологического мониторинга является обеспечение охраны окружающей среды. В зависимости от объектов и задач мониторинга в Единой системе государственного экологического мониторинга выделено 16 подсистем, одна из которых связана с государственным мониторингом земель. Согласно положениям земельного законодательства [2], данный вид мониторинга служит для осуществления рационального использования и охраны земельных ресурсов.
В настоящем исследовании затронуты вопросы применения методов обработки данных дистанционного зондирования Земли (многозональных космических снимков) для мониторинга земель лесного фонда, к которому в силу положений лесного законодательства [3] относятся как лесные, так и нелесные земли. Некоторые профильные задачи по мониторингу лесных земель решаются в рамках государственного мониторинга воспроизводства лесов и государственного лесопатологического мониторинга, являющихся подсистемами государственного экологического мониторинга, а также при проведении мониторинга пожарной опасности в лесах и лесных пожаров.
Вне зависимости от подсистемы государственного экологического мониторинга одним из ключевых источников информации выступают данные дистанционного зондирования Земли, которые позволяют оперативно проектировать и обновлять картографические материалы при помощи использования геоинформационных систем и программных продуктов по дешифрированию снимков.
Цель исследования – апробация методик дешифрирования многозональных космических снимков для целей мониторинга земель лесного фонда и разработки соответствующих картографических материалов на основе применения геоинформационных технологий.
Материалы и методы исследования
В качестве объекта исследования выступила территория Темниковского территориального лесничества Республики Мордовия, находящегося на северо-западе региона в Темниковском, Теньгушевском и Зубово-Полянском районах. Лесничество граничит на севере с Нижегородской областью, на северо-востоке с Мордовским государственным заповедником им. П.Г. Смидовича, на западе с Рязанской областью, на востоке – с Краснослободским и на юге – с Виндрейским территориальными лесничествами. В составе территориального лесничества в соответствии с лесохозяйственным регламентом [4] выделены Барашевское (13163 га), Кечемировское (19360 га), Темниковское (15761 га), Теньгушевское (11777 га), Харинское (15483 га) участковые лесничества.
Базовые картографические слои (границы субъекта, муниципальных районов, кадастровых районов, населенных пунктов, территориального и участковых лесничеств, железные и автомобильные дороги, гидрография) были спроектированы в геоинформационной системе Maplnfo Professional с помощью оцифровки картографических материалов лесохозяйственного регламента лесничества [4] и схемы территориального планирования Республики Мордовия [5].
В работе использованы многозональные космические снимки Landsat-8, 9, полученные с геопортала Геологической службы США (USGS): 2014 г. – от 6 июля; 2018 г. – от 6 июля; 2023 г. – от 27 апреля, 8 июля и 4 октября. Дешифрирование многозональных космических снимков проводилось в программном комплексе QGIS. На первом этапе изображения были радиометрически откалиброваны в соответствии с методикой поставщика данных ДЗЗ [6]. Для этого попиксельно для каждого снимка «сырое» безразмерное нормализованное значение DN (Digital Number) пересчитано в параметр отражательной способности (Reflectance)
В результате проведенного исследования были рассчитаны и проанализированы пространственно-временные изменения числовых расчетных показателей вегетационных индексов – параметров, определяемых в результате операций с спектральными каналами. Для целей анализа динамики растительного покрова были вычислены нормализованный разностный индекс растительности (Normalized Difference Vegetation Index, NDVI), нормализованный разностный водный индекс (Normalized Difference Water Index, NDWI) и хлорофилльный индекс (Green Chlorophyll Index, Clgreen) (табл. 1). Указанные индексы успешно используются исследователями в работах последних лет для мониторинга состояния лесной растительности.
Таблица 1
Формулы для расчета вегетационных индексов
Вегетационный индекс |
Формула |
Переменные |
NDVI |
(NIR – Red) / (NIR + Red) |
NIR – значение спектрального отражения в ближнем инфракрасном канале SWIR – значение спектрального отражения в среднем инфракрасном канале Red – значение спектрального отражения в красном канале Green – значение спектрального отражения в зеленом канале |
NDWI |
(NIR – SWIR) / (NIR + SWIR) |
|
Clgreen |
(NIR / Green) – 1 |
NDVI – один из самых часто используемых вегетационных индексов в современных работах по мониторингу лесных насаждений и оценке фитомассы [7–9]. Впервые его применение было описано в 1973 г. в работе J.W. Rouse, R.H. Haas, J.A. Schell и D.W. Deering. Для растительного покрова значения данного индекса имеют положительные значения. Показатели данного индекса варьируются от –1 до +1 в зависимости от состава представленной растительности, состояния растительного покрова и других факторов. Чем оно выше, тем больше фитомассы регистрируется в одном пикселе космического снимка. Для лесной растительности в середине вегетационного периода показатель, как правило, выше 0,7.
NDWI, предложенный B. Goa в 1996 г., используется с целью мониторинга изменений содержания воды в растениях и поверхностных водах [10, 11]. Амплитуда показаний данного индекса от –1 до +1. С помощью него оценивается влажность в растительности. Залесные территории имеют значение около 0,4.
Индекс Clgreen, активно используемый в сельскохозяйственных исследованиях с начала 2000-х гг., в современных работах также применяется и для мониторинга состояния лесного покрова [12, 13]. Высокие значения индекса связаны с наибольшим содержанием хлорофилла в листьях растений. Значения выше 1 свидетельствуют о здоровой растительности с большим содержанием хлорофилла в листовых пластинах.
Результаты исследования и их обсуждение
На основе полученных картографических материалов (рис. 1–3) и статистических сведений (табл. 2) сделаны выводы об изменении показателей растительного покрова Темниковского территориального лесничества Республики Мордовия.
Оценка общего количества фитомассы по NDVI показывает, что территория Темниковского территориального лесничества соответствует классу «густой растительности» общепринятой шкалы данного индекса. При этом значения с 2014 по 2023 г. имеют тренд на увеличение (с 0,72332 до 0,74273). Данная тенденция характерна для всех участковых лесничеств, при этом для Харинского участкового лесничества значение NDVI изменилось с класса «разреженная растительность» (0,69637) на класс «густая растительность» (0,72985). В 2023 г. во всех пяти участковых лесничествах наблюдаются значения выше 0,7.
Рис. 1. Динамика значений вегетационного индекса NDVI на территории Темниковского территориального лесничества Республики Мордовия в 2023 г.: № 1 – 27 апреля, № 2 – 8 июля и № 3 – 4 октября
Рис. 2. Динамика значений вегетационного индекса NDWI на территории Темниковского территориального лесничества Республики Мордовия в 2023 г.: № 1 – 27 апреля, № 2 – 8 июля и № 3 – 4 октября
Рис. 3. Динамика значений вегетационного индекса Clgreen на территории Темниковского территориального лесничества Республики Мордовия в 2023 г.: № 1 – 27 апреля, № 2 – 8 июля и № 3 – 4 октября
Таблица 2
Результаты расчета вегетационных индексов на территории Темниковского территориального лесничества Республики Мордовия в 2014–2023 гг.
NDVI |
NDWI |
Clgreen |
|||||||
06.07. 2014 г. |
06.07. 2018 г. |
08.07. 2023 г. |
06.07. 2014 г. |
06.07. 2018 г. |
08.07. 2023 г. |
06.07. 2014 г. |
06.07. 2018 г. |
08.07. 2023 г. |
|
Территориальное лесничество |
0,72332 |
0,74004 |
0,74273 |
0,45563 |
0,45317 |
0,4503 |
3,40933 |
3,64539 |
3,40933 |
Участковые лесничества |
|||||||||
Барашевское |
0,74066 |
0,76489 |
0,76153 |
0,45642 |
0,45432 |
0,45269 |
3,63577 |
3,98541 |
3,63577 |
Кочемировское |
0,72535 |
0,72549 |
0,74503 |
0,45215 |
0,4515 |
0,45215 |
3,35746 |
3,47216 |
3,35746 |
Темниковское |
0,72727 |
0,74995 |
0,75089 |
0,45098 |
0,45329 |
0,45098 |
3,526 |
3,79055 |
3,52600 |
Теньгушевское |
0,73129 |
0,74293 |
0,74626 |
0,4597 |
0,45508 |
0,45197 |
3,47542 |
3,60161 |
3,47542 |
Харинское |
0,69637 |
0,72489 |
0,72985 |
0,45141 |
0,45483 |
0,45141 |
3,11561 |
3,45980 |
3,11561 |
Целевое назначение и категории защитных лесов |
|||||||||
Леса, расположенные в защитных полосах лесов |
0,70384 |
0,72186 |
0,7341 |
0,4522 |
0,45168 |
0,4522 |
3,19715 |
3,41789 |
3,39725 |
Запретные полосы лесов, расположенные вдоль водных объектов |
0,73363 |
0,75329 |
0,75706 |
0,45946 |
0,45497 |
0,45946 |
3,66835 |
3,84497 |
3,75406 |
Леса, расположенные в лесопарковых зонах |
0,68204 |
0,7157 |
0,71713 |
0,4502 |
0,45441 |
0,45022 |
2,93606 |
3,24830 |
3,12461 |
Противоэрозионные леса |
0,70246 |
0,71156 |
0,74393 |
0,45129 |
0,4557 |
0,45129 |
3,35788 |
3,43475 |
3,61548 |
Эксплуатационные леса |
0,72461 |
0,74086 |
0,74638 |
0,45211 |
0,45354 |
0,45211 |
3,40354 |
3,64874 |
3,52477 |
Оценка общего количества фитомассы по NDVI показывает, что территория Темниковского территориального лесничества соответствует классу «густой растительности» общепринятой шкалы данного индекса. При этом значения с 2014 по 2023 г. имеют тренд на увеличение (с 0,72332 до 0,74273). Данная тенденция характерна для всех участковых лесничеств, при этом для Харинского участкового лесничества значение NDVI изменилось с класса «разреженная растительность» (0,69637) на класс «густая растительность» (0,72985). В 2023 г. во всех пяти участковых лесничествах наблюдаются значения выше 0,7.
Аналогичная ситуация регистрируется при анализе территории в разрезе целевого назначения. Как для категорий защитных лесов, так и для эксплуатационных лесов показатель NDVI демонстрирует рост. Наибольшие значения в 2023 г. характерны для эксплуатационных лесов (0,74638) и запретных полос лесов, расположенных вдоль водных объектов (0,75706).
NDWI ежегодно имеет высокие показатели на территории всего территориального лесничества (более 0,45). Анализ значений по участковым лесничествам показывает, что за все периоды наблюдений сохраняются стабильные значения, соответствующие уровню «зеленая растительность, насыщенная водой» общепринятой шкалы данного индекса.
В формуле NDWI используется канал SWIR, который реагирует не только на изменения содержания воды в листовых пластинах, но и на условия увлажнения в верхнем слое почвы. Закономерно, что для запретных полос, расположенных вдоль водных объектов, данный показатель несколько выше, чем для остальных категорий. Повышенные значения NDWI в апреле по отношению к октябрю обусловлены распусканием листьев, накоплением в них влаги, а также сохраняющейся влажностью почвы после весеннего снеготаяния.
Вегетационный индекс Clgreen показывает оценку фотосинтетической активности растительного покрова. Чем она выше, тем больше в листьях растений содержится хлорофилла.
Рис. 4. Средние значения вегетационных индексов по лесным кварталам
В целом для территории Темниковского территориального лесничества индекс Clgreen выше 3,41, что свидетельствует о высокой фотосинтетической активности лесных земель. При этом в сравнении 2018 и 2023 гг. показатель несколько снизился (с 3,64539 до 3,40933). Значения по участковым лесничествам в 2018 г.: Барашевское – 3,99, Кочемировское – 3,47, Темниковское – 3,79, Теньгушевское – 3,60, Харинское – 3,46. Результаты анализа лесной растительности по целевому назначению и категориям защитных лесов показали, что наименьшее значение вегетационный индекс принимает в категории лесов, расположенных в лесопарковых зонах, а наибольшее регистрируется для запретных полос лесов, расположенных вдоль водных объектов.
Для более эффективной интерпретации пространственной мониторинговой информации о параметрах залесенных территорий и принятия управленческих решений создан полигональный слой с лесными кварталами, в котором каждой единице присвоено усредненное значение вегетационных индексов (рис. 4).
Результаты расчета вегетационных индексов и автоматизированной контролируемой кластеризации космических снимков методом максимального подобия позволили провести анализ динамики лесопокрытой территории в границах Темниковского территориального лесничества за последние 10 лет в период 2014–2023 гг. Были определены территории, расчищенные от леса в 2014–2023 гг.; территории, залесенные в 2014–2023 гг., и территории без изменений лесного покрова (рис. 5).
По данным геоинформационного мониторинга площадь территорий, расчищенных от леса в 2014–2023 гг., составила 2561,80 га. Площадь территорий, залесенных в 2014–2023 гг., – 2027,10 га. Площадь лесопокрытых территорий в 2014–2023 гг. уменьшилась на 534,70 га. Наибольшие изменения зарегистрированы для Барашевского и Харинского лесничеств.
Рис. 5. Изменение лесистости на территории Темниковского территориального лесничества Республики Мордовия в 2014–2023 гг.
Заключение
В результате проведенного исследования выявлены важнейшие качественные и количественные особенности лесных земель в составе Темниковского территориального лесничества. Спроектированы электронные слои и базы данных, отражающие пространственные характеристики лесопокрытых территорий. Результаты могут быть использованы как для мониторинга земель лесного фонда в рамках государственного мониторинга земель, так и при реализации программ других видов государственного экологического мониторинга, например государственного мониторинга воспроизводства лесов.