Обширное развитие природных пожаров, охвативших различные районы Северной Евразии в последнее десятилетие, привело к осознанию необходимости более детального изучения причин их развития, экологических и социально-экономических последствий [1–4]. Значительный современный интерес к изучению травяных пожаров в степных регионов во многом обусловлен:
а) открывшимися возможностями доступа к материалам дистанционного зондирования Земли (ДЗЗ);
б) проблемами осуществления противопожарных мероприятий на заповедных территориях;
в) необходимостью оценки разнообразных экологических последствий;
г) увеличением ущерба населенным пунктам от распространения огня с прилегающих территорий.
Открытый огонь до сих пор используется сельхозпроизводителями для уничтожения растительных остатков и в качестве доступного способа улучшения кормовых качеств пастбищных и сенокосных угодий [4, 5]. Следует констатировать, что по степным регионам полностью отсутствуют какие-либо официальные сведения о количестве и площадях природных пожаров, что свидетельствует о недостаточном внимании к ним со стороны природоохранных служб и органов власти, несмотря на достаточно разработанную нормативно-правовую базу. Добавим, что травяные пожары, помимо непосредственного воздействия на степные экосистемы и их длительного восстановления, служат поставщиком углерода в атмосферу и влияют на изменение процессов его депонирования почвами [6].
Исходя из того, что природные пожары в настоящее время являются одним из наиболее значимых факторов формирования современных степных ландшафтов и экосистем, наиболее актуальным остается получение достоверных фактических данных, на основе которых возможно проведение анализа пространственно-временных закономерностей их развития и оценки экологических последствий. В рамках исследования получены результаты о возможности применения космоснимков, полученных значительно позже времени развития пожаров, и рассмотрены доступные варианты использования различных данных ДЗЗ.
Материалы и методы исследования
В качестве основного метода получения сведений о площадях гарей было принято визуальное дешифрирование комбинированных снимков Landsat, основанное на экспертной оценке временной серии предпожарных и постпожарных снимков с учетом характерной разницы световых тонов. Методология апробирована на серии ключевых участков Заволжско-Уральского региона; наибольший временной охват был получен по пяти ключевым участкам общей площадью 5098 км². Для этих территорий ранее [4] были дешифрированы гари за 1984–2014 гг., в рамках настоящего исследования временной ряд дополнен гарями за последние годы (2015–2017 гг.). Анализ полученной базы данных позволил выявить максимально возможный многолетний ход развития пожаров и констатировать тенденцию их активизации с начала 2000-х гг. Помимо приведенных ключевых участков визуальное дешифрирование гарей на основе снимков MODIS было проведено для обширного региона (221,4 тыс. км2), охватывающего сектор степной и прилегающих к ней лесостепной и полупустынной зон в пределах Заволжья, Южного Урала, Предуралья и Зауралья, Северного Прикаспия и Приаралья за характерные (2005, 2010, 2017) годы.
Данные визуального дешифрирования были сопоставлены с другими возможными источниками информации:
1) ареалы, полученные в ходе автоматизированной классификации снимков Landsat ГИС-средствами (ENVI, ArcGIS, QGIS) с использованием спектральных преобразований NBR, MIRBI, NDVI и др.;
2) данные «теплового» канала Landsat;
3) архив данных по термоточкам FIRMS (The Fire Information for Resource Management System) за 2001–2017 гг. – результат идентификации тепловых аномалий по снимкам MODIS;
4) изображения, представляющие собой интерпретацию данных по тепловым аномалиям, выполненную SFMS (ScanEx Fire Monitoring Service) по собственному алгоритму для отображения примерных площадей природных пожаров.
Перечисленные источники информации анализировались на различных пространственных уровнях, от крупных мезорегионов до отдельных ареалов гарей. Оценка временных серий (протяженностью до 2–3 лет) классифицированных изображений и данных тепловых каналов Landsat производилась с целью выявления продолжительности и особенностей восстановления степных экосистем. Также эти исходные и преобразованные изображения Landsat рассмотрены в качестве альтернативы визуальному дешифрированию гарей, выявлена степень их достоверности при идентификации площадей пожаров. Анализ продуктов MODIS по тепловым аномалиям проведен, прежде всего, для получения обзорных сведений о степени подверженности пожарам отдельных географических регионов и в качестве потенциальных данных о многолетней динамике пожаров.
Результаты исследования и их обсуждение
Отметим, что все без исключения материалы ДЗЗ в аспекте поставленной цели имеют одни и те же потенциальные недостатки, связанные:
а) с моментальностью и определенной периодичностью съемки (отсутствие возможности непрерывного наблюдения за развитием пожаров);
б) особенностями распространения степных пожаров (кратковременность воздействия огня и, как правило, быстротечность его распространения);
в) неблагоприятными метеоусловиями (облачность) на время получения снимков (необходимость анализа серий снимков, неопределенная достоверность данных по тепловым аномалиям);
г) недостаточно высоким пространственным разрешением космоснимков (в первую очередь MODIS), не позволяющим идентифицировать небольшие по площади гари, преобладающие в количественном выражении.
Несмотря на означенные проблемы, каждый из перечисленных источников информации может быть использован как самостоятельно (для решения каких-либо частных вопросов), так и комплексно при задаче получить максимально полную и достоверную информацию о пространственно-временных аспектах распространения природных пожаров, особенно для обширных территорий.
Опыт визуального дешифрирования
Визуальное дешифрирование космических изображений позволяет получить наиболее детальные для идентификации степных пожаров результаты, но остается весьма трудоемким процессом. С учетом редкой периодичности снимков Landsat одним из этапов оценки достоверности полученных результатов стало выявление максимальной удаленности времени получения космических изображений от времени прохождения пожара. Для анализа была выбрана одна из наиболее продолжительных временных серий, состоящая из 15 снимков Landsat за 2009 г. и одного весеннего снимка за 2010 г. для оценки возможности использования для позднеосенних гарей.
Анализировалось распределение яркости изображения по смежным горелым и негорелым участкам по пяти пожарам, наблюдавшимся в различных местах сцены в разные месяцы пожароопасного периода. Рассматриваемые пары фрагментов подбирались с условием, чтобы каждый из них располагался на однотипной поверхности. Определялось значение яркости (средние, медиана и отклонения) по гистограмме в графическом редакторе по стандартной градиентной шкале, от черного (0) до белого (255).
Исходя из того, что анализу подвергались готовые композитные изображения с различными настройками сведения, для расчетов была принята разность значений по негорелой и горелой частям (рис. 1); дополнительно приведена шкала с шагом в 5 единиц, позволяющая примерно оценить степень контрастности между рассматриваемыми негорелыми и горелыми участками и оценить полученные графики. Отмеченные различия на основе значений яркости лишь частично описывают эти различия, так как на композитных изображениях контрастно заметны и цветовые отличия – естественные участки имеют зеленовато-бежевую окраску, а гари от недавних пожаров – от черной до темно-бурой. Идентификации гарей также способствует, как правило, языковатый характер расположения границ.
Рис. 1. Динамика разностей средних значений яркости (график) и разностей отклонений от средних значений (диаграмма) смежных негорелых и горелых территорий. Даты пожаров: А – 110 день (20 апреля, степь), Б – 123 (3 мая, степь), В – 170 (19 июня, степь), Г – 170 (19 июня, залесенная пойма), Д – 268 (25 сентября). Е – фрагмент шкалы различий яркости, по шкале от 0 (черный) до 255 (белый)
Полученные результаты свидетельствуют, что композитные изображения достаточно информативны для выявления гарей визуальным способом. Естественно, чем ближе снимок по дате к дню прохождения пожара, тем контрастнее различия тонов гарей и фоновых поверхностей. Наиболее существенное снижение различий между горелыми и негорелыми территориями наблюдается в период активного нарастания растительной фитомассы, вплоть до конца мая – начала июня. Соответственно, ранне-весенние (апрель – начало мая) гари могут быть достаточно объективно идентифицированы по снимкам, сделанным не позднее середины мая. Помимо возможного отсутствия снимков за этот короткий период весьма сложна идентификация характерных для региона сельскохозяйственных палов (сжигание стерневых остатков на пахотных угодьях), так как практически сразу же вслед за выжиганием проводится механизированная обработка. Гари летних пожаров достаточно хорошо визуально дешифрируются практически вплоть до начала залегания снежного покрова. Позднеосенние пожары часто не могут быть выявлены из-за отсутствия достаточного количества пригодных для анализа снимков, в том числе из-за неустойчивой облачной погоды и выпадения снежного покрова. Пожары в этот период, как и наблюдаемые в апреле, большей частью связаны с действиями сельхозпроизводителей. В этих случаях вполне возможно использование ранневесенних снимков последующего года.
Показатели отклонений значений яркости от среднего значения до прохождения пожаров имеют значения, близкие к нулевым, так как для анализа подбираются компактные и относительно однородные участки. Вместе с тем выгорание травяной растительности «обнажает» подстилающую поверхность, что в совокупности с горелыми остатками создает более пеструю структуру изображения по яркости (рис. 1).
Исходя из вышеизложенного, минимальное количество снимков для достоверного визуального дешифрирования гарей должно быть не менее 3–4 за год. Дешифрирование одной и той же гари по материалам ДЗЗ с разным пространственным разрешением – высокого (снимки сервиса Google Earth), среднего (Landsat) и низкого (MODIS) – показывает на отсутствие значимых расхождений в ее площади, обычно не ниже 90 %.
Краткие выводы о возможностях использования интерпретированных данных Landsat и MODIS
При сопоставлении многолетних данных визуального дешифрирования снимков Landsat (совокупные площади гарей по пяти ключевым участкам) и данных архива термоточек (суммарное количество, для территории Заволжско-Прикаспийского региона) выявлена высокая степень сходства этих рядов (коэффициент корреляции 0,64) за общий период – 2001–2017 гг. (рис. 2). Это означает, что интерпретация данных архива термоточек остается одним из относительно нетрудоемких способов получения информации о пространственно-временных особенностях развития пожаров. В сочетании с наличием атрибутивной информации по дате и времени каждой термоточки этот источник данных можно использовать при сопоставлении с другими данными (в том числе и в суточном выражении), например для изучения роли погодно-климатических условий.
Рис. 2. Данные многолетней (1984–2017) динамики совокупных площадей степных пожаров и численности термоточек по Заволжско-Прикаспийскому региону
Следует заметить, что быстрое прохождение степных пожаров и немного менее скоротечное остывание поверхности приводит к тому, что тепловые аномалии из снимков MODIS отображают лишь положение огневого фронта на момент съемки. Соответственно, длительно развивающиеся пожары обычно отображены в виде двух и более (по количеству снимков) групп термоточек, соответствующих положению фронта пожара. «Моментальность» съемки и получаемых данных приводит к тому, что в основном фиксируются обширные пожары (площадью более 5 км2), хотя и по таким возгораниям встречаются «пропуски» данных. По этой же причине менее значимые по площади и одновременно наиболее часто отмечающиеся пожары зачастую вообще не фиксируются. Эти недостатки наследуются и интерпретированным данным – площадям активных пожаров, полученным алгоритмами SFMS, значения которых в среднем составляют 25–35 % от совокупной фактической площади гарей и до 40–45 % для отдельных крупных гарей.
Данные теплового канала Landsat и изображения, полученные в ходе проведения спектральных преобразований (NBR, MIRBI, NDVI и др.), позволяют решить две крупных задачи:
1) автоматизированно получить ареалы гарей, близкие к фактическим;
2) проследить продолжительность восстановительных процессов в различных ландшафтно-географических условиях.
В этих целях нами рассмотрена 3-летняя серия снимков, позволившая посредством анализа температурного фона выявить продолжительность и пространственные закономерности постпожарных сукцессий [7]. Было выявлено, что температура поверхности гарей в теплый период в среднем выше на 2,5–4,0 °С; на следующий год эти различия повсеместно остаются, снижаясь до 1,5–2,0 °С; и лишь на третий год после пожара (к окончанию полного второго вегетационного периода) различия практически нивелируются (в среднем 0,2 °С).
Серией исследований нами также выявлено, что выгоревшие участки могут быть эффективно дешифрированы применением специализированных индексов NBR и MIRBI; менее достоверно использование нормализованного относительного индекса растительности NDVI и метода классификации ISODATA [8]. При этом в период активного нарастания фитомассы (со второй половины мая до конца июня) достоверность резко снижается, но затем вновь достигает 80–95 % от фактических площадей. Главным недостатком использования автоматизированных методов анализа материалов ДЗЗ в степных регионах является необходимость проведения дополнительных работ по выявлению возделываемых полей, поверхность которых дает «отклик», неотличимый от сгоревших территорий.
Заключение
Визуальное дешифрирование гарей для степных регионов остается одним из наиболее достоверных способов получения фактической информации, особенно для регионов с разреженным травянистым покровом [9]. Методология автоматизированной идентификации гарей геоинформационными средствами постоянно развивается [10], но с учетом особенностей распространения степных пожаров применима лишь для отдельных аспектов интересующей проблемы. Идентификация гарей для обширных степных регионов возможна при условии разработки алгоритма многоступенчатой обработки материалов ДЗЗ, в том числе учитывающего ландшафтную неоднородность и характер землепользования.
Исследование выполнено в рамках гранта РФФИ № 18-05-00088 «Природные пожары в степных регионах Евразии: анализ пространственно-временных закономерностей и геоэкологических последствий».