Волго-Ахтубинская пойма является значимым природным комплексом, в котором леса играют важную роль [1]. Пойменные леса являются жизненным пространством для множества растений и животных. Они предоставляют укрытие, пищу и места гнездования для многих видов, включая редкие и уязвимые. Биологическое разнообразие пойменных лесов способствует поддержанию экосистемного баланса и предоставляет ценные экологические услуги [2]. Однако пойменные леса также подвержены угрозам, особенно в контексте гидрологических и климатических изменений. Снижение частоты и уровней половодья приводит к стабилизации грунтовых вод, накоплению солей в зоне аэрации, что приводит к деградации аборигенных растительных сообществ. Особенно губительны эти изменения для дубрав и ивняков, которые замещаются инвазивными видами [3]. Изменения климата, такие как повышение температуры и засухи, могут снизить устойчивость лесов и способствовать распространению пожаров и вредителей. Это может привести к деградации и потере пойменных лесов, что имеет негативные последствия для биоразнообразия.
Для сохранения и защиты пойменных лесов в Волго-Ахтубинской пойме необходимо проводить научные исследования по оценке их актуального состояния, разрабатывать эффективные стратегии управления и принимать меры по адаптации к гидрологическим и климатическим изменениям. Важно осознавать ценность и уникальность пойменных лесов и принимать меры для их сохранения для экологической устойчивости региона. Технологии дистанционного зондирования Земли из космоса широко применяются для картографирования лесных насаждений [4, 5], тем не менее даже в пределах Волгоградской области методы спутникового мониторинга пойменных лесов недостаточно проработаны. Имеющиеся исследования чаще основаны на очень трудоемком экспертном дешифрировании данных сверхвысокого пространственного разрешения [6–8], которые обновляются достаточно редко в рамках общедоступных сервисов либо имеют высокую стоимость.
Наличие доступных электронных карт типов земного покрова, полученных по спутниковым данным высокого пространственного разрешения, открывает возможность картографирования различных категорий земель на региональном и муниципальном уровнях. Первым информационным продуктом типов земного покрова разрешения 10 м/пикс. стал FROM-GLC10, изначально основанный на спутниковых снимках Landsat разрешения 30 м/пикс. Появление миссии Sentinel-2 позволило получить данные пространственного разрешения 10 м/пикс., что сделало возможным улучшение пространственного разрешения данных FROM-GLC10 до 10 м/пикс. [9]. Данные ESRI Land Cover, основанные на спутниковых изображениях Sentinel-2, имеют временной охват 2017–2022 гг. [10], что открывает возможности для анализа изменений в землепользовании на этом отрезке времени. Данные ESA WorldCover доступны только на 2020 г., но основаны как на оптических спутниковых изображениях Sentinel-2, так и на радиолокационных данных Sentinel-1 [11], что повышает их точность при картографировании защитных лесных насаждений. Также существует еще один глобальный информационный продукт разрешения 30 м/пикс. – Global Land Cover 30 (GLC30), но он показал значительно меньшую точность при картографировании лесных насаждений, а леса на пойменных землях по этим данным относятся к водно-болотным угодьям без выделения древостоев [12]. По этим причинам в данной работе GLC30 не рассматривается.
Основные типы земного покрова, которые выделяются по данным высокого пространственного разрешения: пашни, леса, травянистые экосистемы (grassland), кустарники, водоемы, искусственные поверхности (застройка), открытые почвы, вечные льды и снега. Также некоторые информационные продукты выделяют водно-болотные угодья, в том числе мангры. Достоинствами указанных глобальных данных о типах земного покрова являются их общеземной охват, свободный доступ, высокое пространственное разрешение. К недостаткам этих данных можно отнести относительно небольшое количество классов (8–10), ограниченность глубины архива (1–5 лет), нерегулярное обновление, существенные различия в точности у разных источников данных при картографировании разных типов земного покрова.
Кроме информационных продуктов типов земного покрова для картографирования лесов используются специализированные данные. Одним из наиболее распространенных источников является Global Forest Change (GFC) пространственного разрешения 30 м/пикс. [13]. Достоинством этого информационного продукта является большая глубина архива: данные о лесопокрытой площади на 2000 г. с учетом ее изменений за 2001–2022 гг.
Цель данного исследования – провести сравнительный анализ результатов картографирования лесов в пойменных ландшафтах Нижней Волги в пределах Волгоградской области по данным глобальных информационных продуктов Global Forest Change, FROM-GLC10, ESRI, ESA. Для достижения цели решались следующие задачи: разработка карт лесов по данным разных информационных продуктов, создание результирующего растра по данным всех источников, сопоставление лесистости в разрезе муниципальных образований с определением точности на основе коэффициента детерминации, относительной и среднеквадратической ошибок.
Материалы и методы исследования
Исследование проводится в северной части Волго-Ахтубинской поймы в границах Волгоградской области (рис. 1). Площадь исследуемой территории составляет 185,23 тыс. га. В работе кроме данных о типах земного покрова ESA, ESRI, FROM-GLC10 были использованы данные Global Forest Change – 2022. Продукт основан на автоматическом распознавании лесов по мультиспектральным снимкам Landsat в видимом и инфракрасном диапазонах спектра и имеет пространственное разрешение 30 м. GFC-2022 – это набор данных, который разделен на тайлы размером 10х10 градусов в географической проекции в трех растровых слоях: «treecover2000» – слой леса, «loss» – потеря и «gain» – прирост. Древесный покров на 2000 г. определялся как закрытие полога для всей растительности высотой более 5 м. В слое «treecover2000» от 1 до 100 указывается процент древесного покрова в пикселе. За покрытую лесом площадь принимались пиксели GFC с проективным покрытием 30 и более процентов, вычитались пиксели исчезнувшего в 2000–2022 гг. леса, добавлялись пиксели восстановившихся насаждений.
Все данные предоставляются в формате GeoTiFF и общеземной географической системе координат, только данные ESRI спроецированы в универсальную поперечную проекцию Меркатора (UTM). Это позволяет свободно использовать данные с использованием популярных геоинформационных программ, например QGIS.
Рис. 1. Лесной покров по данным разных информационных продуктов (а – GFC, б – ESA, в – ESRI, г – FROM-GLC10)
Типы земного покрова по данным различных информационных продуктов в растровых файлах закодированы различными цифровыми значениями: у ESRI Land Cover пиксели леса имеют значение 2, у ESA World Cover лес обозначен значением 10, а FROM-GLC10 – 20. Поэтому были созданы бинарные маски лесного покрова согласно каждому источнику данных, где лесным пикселям присваивалось значение 1, а всем остальным – 0.
Поскольку каждый информационный продукт характеризуется различными ошибками, как пропусками, так и ложными определениями лесов [12], был разработан результирующий растровый слой, пиксели которого содержали количество информационных продуктов, относивших данный пиксель к лесопокрытой площади. За достоверные пиксели леса в этом случае принимались пиксели по данным двух и более информационных продуктов. После чего было проведено сравнение площадей леса по данным отдельных продуктов и результирующего растра в разрезе муниципальных образований, что позволило оценить точность картографирования леса на основе разных источников сведений о типах земного покрова. Для сравнительного анализа использовался показатель лесистости – отношение площади леса к общей площади территории, т.е. доля лесных насаждений в общей площади исследования. Статистическая обработка выполнена в Microsoft Excel.
Результаты исследования и их обсуждение
В результате сравнения на территории Волго-Ахтубинской поймы общая площадь лесного покрова составляет от 19,9 тыс. га по данным ESRI до 41,1 тыс. га по данным GFC, что соответствует данным экспертного дешифрирования, согласно которым площадь леса в Волго-Ахтубинской пойме без волжских островов и правобережной прибортовой поймы составляет 34,4 тыс. га [7]. Согласно результирующему растру по данным всех информационных продуктов площадь лесов составила 34,1 тыс. га, что меньше оценки по экспертному дешифрированию, тем более учитывая площади островов.
Средняя лесистость полуострова Сарептский по данным информационных продуктов типов земного покрова составила от 60 до 74 %, что соответствует данным О.Ю. Кошелевой по результатам экспертного дешифрирования, согласно которым лесистость составляет 64,2 % [8]. Лесистость на территории острова Сарпинский по данным информационных продуктов составила от 22 до 33 %, по данным В.Н. Бодровой на 2018 г. средняя лесистость острова составляла около 40–45 % [6]. Такая разница, возможно, связана с недостаточной точностью информационных продуктов или вспышкой численности ясеневой узкотелой златки после 2018 г., из-за которой множество деревьев на острове погибло [14, 15].
Сравнительный анализ значений лесистости по данным различных информационных продуктов и составленного на их основе результирующего растра (рис. 2) показал, что наиболее тесную связь по величине коэффициента детерминации имеют данные ESRI (R2 = 0,87), также у этих данных угловой коэффициент наиболее близок к единице. Тем не менее данные ESRI существенно занижают лесистость в сравнении с результирующим растром (относительная ошибка – -42 %), а также имеют достаточно большую среднеквадратическую ошибку (RMSE = 8,27 %). Следующим информационным продуктом по величине коэффициента детерминации является FROM-GLC10 (R2 = 0,82), в среднем по этим данным лесистость завышена на 6,3 % при среднеквадратической ошибке RMSE = 3,94 %. Информационные продукты ESA и GFC имеют сопоставимые коэффициенты детерминации (R2 = 0,62, R2 = 0,59 соответственно) и среднеквадратическую ошибку (RMSE = 6,48 и 6,3 % соответственно) при сравнении с результирующими данными. Тем не менее относительная ошибка при картографировании по данным GFC составляет 19,6 %, а по данным ESA – -13,4 %.
Таким образом, данные FROM-GLC10 обладают наилучшей точностью среди исследованных продуктов. Данные ESRI хорошо коррелированы с результирующим растром, но существенно занижают площадь лесов, что может свидетельствовать о важности правильного определения лесных насаждений и минимального количества ложных определений. По данным GFC площади лесов существенно завышены, что может быть связано с большим количеством зарослей из тростника и других высокостебельных трав, которые зачастую ошибочно определяются как лес при использовании автоматизированных методов обработки данных дистанционного зондирования Земли [12]. Данные ESA также имеют недостаточно высокую точность, занижая площадь леса по сравнению с результирующим растром.
Рис. 2. Сравнение долей лесопокрытой площади по данным информационных продуктов (вертикальная ось) с результирующей величиной по данным двух и более источников (горизонтальная ось)
Заключение
В результате исследований была определена лесистость пойменных ландшафтов Нижней Волги в пределах Волгоградской области, полученная величина площади лесов 34,3 тыс. га согласуется с результатами экспертного дешифрирования. Наибольшей точностью характеризуются данные FROM-GLC10 (R2 = 0,82, RMSE = 3,94 %), в то время как остальные исследованные информационные продукты по отдельности не обеспечивают достаточную точность при картографировании пойменных лесов. Тем не менее построение результирующего растра на основе всех источников данных позволило получить электронную карту лесов в пойменных ландшафтах необходимой точности.
Разработанные электронные карты лесных насаждений Волго-Ахтубинской поймы позволят провести мониторинг состояния древостоев как на основе данных дистанционного зондирования Земли из космоса, так и с помощью традиционных наземных методов лесной таксации. При регулярном обновлении исходных информационных продуктов типов земного покрова станет возможным постоянный мониторинг площади лесных насаждений, что даст возможность определить ключевые факторы динамики их состояния.