ИСПОЛЬЗОВАНИЕ РОССИЙСКИХ ДАННЫХ ДЗЗ ДЛЯ ИЗУЧЕНИЯ ДИНАМИКИ ВОССТАНОВЛЕНИЯ ЭКОЛОГО-РЕСУРСНОГО ПОТЕНЦИАЛА ЛЕСНЫХ РЕГИОНОВ ПОСЛЕ ВОЗДЕЙСТВИЯ ПОЖАРОВ

Полный текст

Вопрос прогнозирования и мониторинга пожароопасной обстановки широко изучается на территории каждой страны. Разрабатываются наиболее функциональные и систематизированные средства мониторинга в целях более раннего обнаружения точек возгорания и предупреждения чрезвычайных экологических ситуаций, связанных с лесными пожарами [3]. В связи с непрерывным развитием и модернизацией космических технологий, в частности в сфере дистанционного зондирования Земли, и невозможностью полного предупреждения риска возникновения таких масштабных экологических бедствий и катастроф, как лесные пожары, актуально совершенствование метода космического мониторинга пожароопасных ситуаций для поиска новых более эффективных подходов к изучению причин возникновения пожаров и оценки их последствий [2]. Цели работы - определить скорость восстановления биомассы территории пострадавшей в результате воздействия лесных пожаров и оценить изменения эколого-ресурсного потенциала с помощью космических средств дистанционного зондирования. При том, что сам принцип мониторинга процесса зарастания гарей после пожаров и оценки биомассы с помощью вегетационных индексов не является новинкой (П.А. Украинский [3], Рожков Ю.Ф., Кондакова М.Ю. [4]), впервые делается попытка объединить их для определения последствий пожаров и скорости восстановления территории (в том числе и с пересчетом в тонны образовавшейся фитомассы). МАТЕРИАЛЫ И МЕТОДЫ ИССЛЕДОВАНИЯ Исследуемый объект располагается на территории Борского и Лысковского района Нижегородской области (рис. 1). Данная область в значительной степени повреждена в результате воздействия катастрофически крупных пожаров в 2010 г. Лесные ресурсы данного района занимают 40-45% территории, отличаются преобладанием хвойных пород деревьев и составляют 2% от всех лесных запасов России, что служит высоким показателем, так как их большая часть сосредоточена на востоке страны (Сибирь и Дальний Восток). Рис. 1. Географическое положение исследуемой территории [Fig. 1. Geographical location of the study area] В соответствии с принятой классификацией леса исследуемых районов области относятся к высокому классу пожарной опасности (средний класс природной пожарной опасности по области - 2,5: Борский район - 2 класс, Лысковский - 3 класс пожарной опасности) [5], о чем свидетельствуют последствия пожаров (рис. 2). Рис. 2. Последствия пожаров 2010 г. на исследуемой территории [Fig. 2. Consequences of fires in 2010 in the study area] Для расчета количества восстанавливающейся растительной биомассы на территории в рамках исследования применен нормализованный относительный индекс растительности NDVI (Normalized Difference Vegetation Index), основанный на спектральных свойствах растительности [6]. Индекс может принимать значения от -1 до 1. Для растительности индекс NDVI принимает положительные значения, обычно от 0,2 до 0,8 (табл. 1) и рассчитывается по формуле [1]: NDVI = rNIR -rred . rNIR + rred Следовательно, вегетационный индекс способен дать относительную оценку растительного покрова территории, которая в дальнейшем, при использовании наземной информации и полевых данных, может быть пересчитана в абсолютные значения [8]. Дискретная шкала NDVI [7] Таблица 1 Тип объекта Отражение в красной области спектра Отражение в инфракрасной области спектра Значение NDVI Густая растительность 0,1 0,5 0,7 Разряженная растительность 0,1 0,3 0,5 Открытая почва 0,25 0,3 0,025 Discrete NDVI scale Table 1 Object type Reflection in the red spectral Reflection in the infrared spectral range NDVI value Thick vegetation 0,1 0,5 0,7 Discharged vegetation 0,1 0,3 0,5 Open soil 0,25 0,3 0,025 В качестве исходных использованы данные, полученные в июле-сентябре 2009, 2010, 2011, 2012 и 2013 гг. (рис. 3) аппаратурой КМСС со спутника «Метеор-М» № 1 (предоставлены Научным центром оперативного мониторинга Земли) АО «Российские космические системы») и данные Landsat-5 (http://glovis.usgs.gov/ index.stml). Данные обработаны с использованием программ ERDAS Imagine, ENVI, ArcGIS. Карта и гистограмма на основе индекса NDVI № 1. Landsat-5 (10 августа 2009 года) Map and histogram based on NDVI № 1. Landsat-5 (August 10, 2009) Карта и гистограмма на основе индекса NDVI № 2. «Метеор-М» № 1 (15 августа 2011 года) Map and histogram based on NDVI № 2. «Meteor-M» № 1 (August 15, 2011) Карта и гистограмма на основе индекса NDVI № 3. «Метеор-М» № 1 (17 сентября 2012 года) Map and histogram based on NDVI № 3. «Meteor-M» № 1 (September 17, 2012) Карта и гистограмма на основе индекса NDVI № 4. «Метеор-М» № 1 (14 июля 2013 года) Map and histogram based on NDVI № 4. «Meteor-M» № 1 (July 14, 2013) Рис. 3. Построение карт и гистограмм на основе NDVI исследуемой территории [Fig. 3. Mapping and histograms building based on the NDVI of the study area] РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЯ И ИХ ОБСУЖДЕНИЕ В результате анализа гистограмм по картам на основе индекса NDVI с 2009 по 2013 гг. видна определенная тенденция: в 2009 г. фитомасса исследуемой территории и общее количество растительности находилось на своем пике. После пожаров 2010 г. значительная часть фитомассы была потеряна, и с 2011 по 2013 гг. отмечается ее не равномерная, но стойкая тенденция к восстановлению. Используя результаты анализа гистограмм на основе NDVI (см. рис. 2), и справочные данные официальных документов Рослесхоза, IISS (Международный институт прикладного системного анализа), лесного плана Нижегородской области, паспортов Борского и Лысковского районов Нижегородской области, в несколько этапов произведен усредненный количественный расчет прироста биомассы на выгоревшей территории с 2011 по 2013 гг. 1. Расчет общей площади лесов и фитомассы исследуемой территории. Согласно материалам МИПСА [6] на территорию РФ приходится около 69 671 150 000 т фитомассы. Нижегородская область обладает 418 220 000 т фитомассы страны. Исходя из площади Нижегородской области и площадей обоих исследуемых районов - Борского и Лысковского, был произведет расчет количества фитомассы, приходящийся на каждый из районов (табл. 2). Таблица 2 Table 2 Площадь и объем фитомассы Нижегородской области и Борского, Лысковского районов The area and volume of phytomass of the Nizhny Novgorod region and Borsky, Lyskovsky districts Территория Площадь территории S, км2 Фитомасса территории Ф, т Нижегородская область 76 900 418 220 000 Борский район 3584 19 491 553,71 Лысковский район 2134 11 605 740,96 Territory Area of territory S, km2 Phytomass of the territory Ф, t Nizhny Novgorod region 76 900 418 220 000 Borsky district 3584 19 491 553,71 Lyskovsky district 2134 11 605 740,96 90000 80000 70000 Площадь, км2 60000 50000 40000 76900 37580 Условные обозначения: o общая площадь территории; o площадь лесной территории 30000 20000 10000 0 3584 2509 2134 960 Нижегородская область Борский район Лысковский район Рис. 4. Общая площадь и площадь лесов Нижегородской области, Борского и Лысковского районов (Источники: расчеты авторов и работа [10]) [Fig. 4. Total area and area of forests of the Nizhny Novgorod region, Borsky and Lyskovsky districts (Source: calculations of the authors and [10])] Леса занимают приблизительно 53% территории области, причем на севере территории лесистость достигает 80%, а в юго-восточный районах снижается до 1% (рис. 4) [9]. По количеству фитомассы, приходящейся на каждый из районов, и площади лесов каждого района было рассчитано количество фитомассы на единицу площади каждого района (табл. 3). Количество фитомассы на единицу площади каждого района Таблица 3 Район Площадь лесов, км2 Фитомасса на единицу площади, т/км2 Общая фитомасса, т Общая фитомасса на единицу площади, т/км2 Борский 2509 (70%) 7768,6 31 097 294,67 9 928,9 Лысковский 960 (45%) 12089,3 The amount of phytomass per area of each district Table 3 District Forest area, km2 Phytomass per area, t/km2 Total phytomass, t Total phytomass per area, t/km2 Borsky 2509 (70%) 7768,6 31 097 294,67 9 928,9 Lyskovsky 960 (45%) 12089,3 2. Расчет количества восстанавливающийся фитомассы с использованием данных гистограмм, построенных по картам на основе вегетационного индекса NDVI. При подсчете прироста фитомассы с использованием гистограмм, было необходимо учесть, что значения, отложенные по оси ординат на гистограммах, соответствуют количеству пикселей. Аппаратура КМСС спутника «Метеор-М» № 1 имеет разрешение 60 м, т.е. один пиксель соответствует 3600 м2. Тогда фитомасса территории: Ф3600 9928,9 т×3600 м 2 = 10000 м 2 = 3574,4 т. (1) Таким образом, один пиксель снимка (3600 м2) содержит приблизительно 3 574,4 т фитомассы. Важно учитывать, что данное значение характерно для идеального состояния растительности на территории, следовательно, необходимо провести градацию количества фитомассы, в зависимости от значения NDVI: - 0,1-0,3 - <100/3600 т/м2; - 0,3-0,4 - 1500/3600 т/м2; - 0,5-1 - 2500-3574,4/3600 т/м2. В качестве эталонного значения NDVI, на уровне которого в дальнейшем производилось сравнение по годам, взято значение 0,4 (в соответствии с анализом неповрежденной территории на основе данных за 2009 год), что соответствует удовлетворительному состоянию растительного покрова и количеству фитомассы, равному 1500 т на пиксель. 2011 год. Согласно гистограмме (см. рис. 3), количество пикселей, соответствующих значению 0,4 индекса NDVI, составляет 6 252, соответственно: Ф2011 = 6252 · 0,0036 км2 · 1500 т = 33760 т, (2) т.е. в 2011 г. общая масса растительности на исследуемой территории составляла 33 760 т. 2012 год. Согласно гистограмме (см. рис. 3), количество пикселей, соответствующих значению 0,4 индекса NDVI, составляет 77. Для данного года был введен поправочный коэффициент в связи с осенним сезоном [6], который составил 0,09 единиц индекса NDVI. Таким образом, в августе 2012 г. на значение 0,4 индекса NDVI приходилось порядка 7106 пикселей, следовательно: Ф2012 = 7106 · 0,0036 км2 · 1500 т = 38372 т., (3) т.е. в 2012 г. общая масса растительности на исследуемой территории составляла 38 372 т. 2013 год. Согласно гистограмме (см. рис. 3), количество пикселей, соответствующих значению 0,4 индекса NDVI, составляет 8 335, соответственно: Ф2013 = 8335 · 0,0036 км2 · 1500 т = 45900 т., (4) т.е. в 2013 г. общая масса растительности на исследуемой территории составляла 45 900 т. Обобщая расчеты (табл. 4) можно сделать вывод о том, что увеличение фитомассы с 2011 по 2012 гг. составило 4612 т, а с 2012 по 2013 гг. - 7528 т, т.е. общее количество фитомассы за 2 года увеличилось приблизительно на 36%. Объем фитомассы исследуемого объекта по годам [Table 4. The volume of phytomass of the object by years] Таблица 4 Год [Year] Объем фитомассы, т [Phytomass volume, t] Прирост фитомассы по сранению с прошлым годом [Growth of phytomass compared to last year] 2011 33 760 - 2012 38 372 14% 2013 45 900 19% 50000 45000 Фитомасса, т 40000 35000 30000 25000 2011 2012 2013 Год Рис. 5. График изменения количества фитомассы с 2011 по 2013 гг. [Fig. 5. Schedule of changes in the amount of phytomass from 2011 to 2013] Таким образом, на основе анализа полученных результатов (рис. 5) можно выявить увеличение биомассы лесной растительности на исследуемой территории за 3 года, прошедших после катастрофически крупных пожаров 2010 г., что свидетельствует о положительной динамике восстановления эколого-ресурсного потенциала лесных хозяйств на территории Нижегородской области. Предлагаемая методика позволяет оценить состояние и эколого-ресурсный потенциал, биомассу (в частности фитомассу) территории и ее продуктивность с использованием российских данных дистанционного зондирования.

Об авторах

Елизавета Андреевна Григорец

Российский университет дружбы народов

Автор, ответственный за переписку.
Email: 5977749@mail.ru

аспирант кафедры экологического мониторинга и прогнозирования Российского университета дружбы народов

Российская Федерация, 113093, Москва, Подольское шоссе, 8/5

Дарья Олеговна Капралова

Российский университет дружбы народов

Email: dorris@yandex.ru

старший преподаватель кафедры экологического мониторинга и прогнозирования Российского университета дружбы народов

Российская Федерация, 113093, Москва, Подольское шоссе, 8/5

Лариса Ивановна Пермитина

НЦ ОМЗ АО «Российские космические системы»

Email: permitina@ntsomz.ru

ведущий научный сотрудник НЦ ОМЗ АО «Российские космические системы»

Российская Федерация, 127490, Москва, ул. Декабристов, 51/25

Список литературы