Термин «мониторинг» со второй половины ХХ в. все активнее используется в научной литературе в связи с постоянным ростом хозяйственного освоения территорий и акваторий земного шара. Согласно [1, c. 6] под мониторингом понимается «система повторных целенаправленных наблюдений за исследуемыми объектами в пространстве и времени». Географические информационные системы (ГИС) и космические снимки высокого разрешения в последние годы стали все активнее использоваться в целях оперативного мониторинга различных природных и антропогенных объектов [2–5]. Использование космических снимков для целей мониторинга представляет особенно актуальную задачу сейчас, при постоянно увеличивающихся потоках информации и необходимости принятия решений. Особенно остро это проявляется для территорий, которые активно используются и вовлекаются в хозяйственное развитие.
Крымский полуостров после 2014 г. испытывает особенно мощный процесс преобразования, развития и модернизации инфраструктуры, масштабы которого уже сложно оценивать традиционными методами. В связи с этим постоянно нарастает нагрузка на экосистемы и ландшафты Крымского полуострова.
Цель исследования – показать возможность применения мультиспектральных космических снимков для целей мониторинга ландшафтов в пределах Бахчисарайского района Республики Крым.
Материалы и методы исследования
Для мониторинга за вегетационным состоянием растений и определения температуры поверхности использованы космические снимки Landsat-8, которые имеют разрешение 30х30 м/пиксель, и содержат набор мультиспектральных каналов, которые являются хорошим средством для ведения мониторинга. Существенным недостатком, при высоком разрешении, является периодичность съемки, которая в среднем составляет два снимка в месяц. Однако с учетом того, что над территорией Крымского полуострова наблюдается значительная облачность, особенно в зимние месяцы, подобрать серию снимков за год не представляется возможным. В связи с этим в работе использованы выборочные космические снимки Landsat-8 за весну – осень 2020 г.
Для целей мониторинга по космическим снимкам возможно рассчитывать множество показателей, характеризующих ту или иную территорию [6, 7]. В работе анализируется вегетационный индекс NDVI и температура поверхности, рассчитанные на основе космических снимков Landsat-8. Методика исследования изложена в работах [8, 9].
Расчет вегетационного индекса NDVI происходит по формуле (1) [10]:
NDVI = 
(1)
где NIR – отражение в ближней инфракрасной области спектра, RED – отражение в красной области спектра.
Расчет температуры поверхности производится с использованием методики, описанной в работах [11–13], по формуле (2), с использованием 10-го канала мультиспектрального архива снимка Landsat-8, ввиду лучших результатов по сравнению с 11-м каналом, как показано в работе [14].
T = (1321.08 / (Ln ((774.89 / ((«B10.TIF» *
*0.0003342) + 0.1)) + 1))) – 273.15, (2)
где «B10.TIF» – 10-й канал мультиспектрального архива снимка Landsat-8 за определенную дату.
Обработка и расчеты произведены с использованием программного комплекса QGIS. Ландшафтная структура приводится по ландшафтно-типологической карте Г.Е. Гришанкова [15] (рис. 1).
Рис. 1. Фрагмент ландшафтно-типологической карты Г.Е. Гришанкова на территорию Бахчисарайского района Республики Крым (номера ландшафтов приводятся согласно [15])
Результаты исследования и их обсуждение
В результате исследования для территории Бахчисарайского района Республики Крым получены значения вегетационного индекса NDVI за 2020 г. (рис. 2) на выбранные даты, имеющие минимальное облачное покрытие.
Рис. 2. Значение вегетационного индекса NDVI за 2020 г. в пределах Бахчисарайского района Республики Крым
Как видно из рис. 2, максимальные значения NDVI характерны для юго-восточной части Бахчисарайского района, где сосредоточены природные лесные экосистемы, в то время как минимальные – для северо-западной части, где активно ведется хозяйственная деятельность и на сельскохозяйственных полях еще не наступила активная фаза вегетации и центральной части, где проходят основные магистрали и расположены крупные населенные пункты. В то же время, если рассматривать сезонную динамику, то в летние месяцы наблюдаются максимальные значения NDVI, которые достигают в юго-восточной части Бахчисарайского района значений более 0,8. Анализ значений NDVI показывает, что максимальные значения наблюдаются в пределах групп местной с естественной непреобразованной растительностью.
Особый интерес представляет распределение значений NDVI в пределах групп местностей ландшафтов в пределах Бахчисарайского района Республики Крым. На рис. 3 представлена карта средних значений вегетационного индекса NDVI в пределах Бахчисарайского района Республики Крым для групп местностей.
Рис. 3. Средние значения NDVI в пределах групп местностей на ландшафтно-типологической карте в пределах Бахчисарайского района Республики Крым
Рис. 4. Динамика изменения вегетационного индекса NDVI в пределах группы местностей № 39 (денудационно-останцовые и аккумулятивные равнины с ковыльно-типчаковыми степями) и № 62 (эрозионное овражно-балочное низкогорье с дубовыми лесами и лесокустарниковыми зарослями)
Как можно видеть из рис. 2, в пределах групп местностей наибольшие средние значения вегетационного индекса NDVI приходятся на группы местностей среднегорно-склонового пояса дубовых, можжевелово-дубовых и смешанных широколиственных лесов. Если рассматривать значение NDVI в пределах восстановленных ландшафтов, представленных на территории Бахчисарайского района, то этот вегетационный индекс косвенно показывает преобразованность ландшафтного контура, ввиду того, что значительно увеличивается неоднородность распределения NDVI. Например, группа местностей денудационно-останцовые равнины с разнотравными степями, зарослями типа «дубки» и колючекустарниковыми зарослями типа «шибляк» характеризуется средними значениями NDVI, которые меньше, чем в большинстве других групп местностей в пределах пояса лесостепи на останцово-денудационных, структурных денудационных и аккумулятивных равнинах, куэстовых возвышенностях, ввиду того что по ней проходит большая часть дорожной инфраструктуры, связывающая Симферополь и Севастополь, а также расположены населенные пункты. На рис. 4 представлен график изменения вегетационного индекса NDVI за рассматриваемый период, из которого видно, что максимальные значения присутствуют практически в течение всего года для контуров с лесной растительностью, а минимальные – для контуров, где активно идет хозяйственное использование территории и выращивают сельскохозяйственные культуры.
Мониторинг температуры поверхности, выполненный на основании космических снимков Landsat-8, представлен на рис. 5.
Рис. 5. Значение температуры поверхности (LST) за 2020 г. в пределах Бахчисарайского района Республики Крым
Для анализа температуры поверхности, как и вегетационного индекса NDVI, снимки за 30 апреля 2020 г., 17 июня 2020 г., 03 июля 2020 г. и 07 октября 2020 г. подходят частично, так как имеют минимальное перекрытие облаками. Как видно из рис. 4, в пределах Бахчисарайского района Республики Крым максимальные температуры поверхности характерны для северо-запада, где имеется большое количество открытых участков, а минимальные – для юго-восточной части района, где преобладают леса. К тому же существенный вклад в уменьшение температуры поверхности с северо-востока на юго-запад вносит увеличение абсолютных высот. Наибольшую наглядность показывает рис. 6, где представлено распределение средних температур поверхности в пределах групп местностей, где вышеуказанная закономерность отчетливо видна. Расчет коэффициента корреляции между средними значениями абсолютной высоты и температуры поверхности в пределах групп местностей показал значения -0,85 для 04 августа 2020 г. и -0,82 для 21 сентября 2020 г., что свидетельствует о наличии высокой отрицательной связи между абсолютной высотой и температурой поверхности.
Рис. 6. Среднее значение температуры поверхности в пределах групп местностей на ландшафтно-типологической карте в пределах Бахчисарайского района Республики Крым
В то же время корреляционный анализ показывает наличие очень высокой обратной отрицательной связи между средними значениями NDVI и температурой поверхности. Для 04 августа 2020 г. коэффициент корреляции составляет -0,92, а для 21 сентября 2020 г. составляет -0,91.
Таким образом можно говорить, что распределение вегетационного индекса NDVI и температуры поверхности (LST) носит квазизональный характер и нарушается только в пределах участков, где территория активно преобразована и используется в хозяйственной деятельности (населенные пункты, транспортная инфраструктура, земли сельскохозяйственного использования) и отчетливо выделяется при мониторинге.
Заключение
Использование космических снимков высокого разрешения позволяет вести мониторинг различных территорий, в том числе и ландшафтов и экосистем, при условии отсутствия данных полевых натурных исследований и слабом покрытии территории Крымского полуострова сетью метеостанций. Полученные данные репрезентативны и позволяют оценивать современное состояние и динамику различных показателей для исследуемых территорий. В результате исследования установлено, что на территории Бахчисарайского района Республики Крым значения вегетационного индекса NDVI увеличиваются, а температуры поверхности уменьшаются с северо-запада на юго-восток, что связано с увеличением абсолютных высот и уменьшением антропогенной преобразованности территории.
Исследование выполнено в рамках темы НИР «Изучение пространственно-временной организации водных и сухопутных экосистем с целью развития системы оперативного мониторинга на основе данных дистанционного зондирования и ГИС-технологий. Регистрационный номер: 121040100327-3».