book_5_20.html

<!DOCTYPE html> 
<html>
<head>
      <title>Моделирование трансформации землепользования</title>
      <meta charset="utf-8">
<meta name="viewport" content="width=device-width, initial-scale=1">
<meta name="description" content="Трансформация землепользования в сельской местности, Моделирование степени зарастания сельскохозяйственных угодий" />
    <link rel="canonical" 
href="https://www.example.com/keywords.html" >
    <link rel="stylesheet" type="text/css" href="style.css">
   <link rel="preload" href="c:/fonts/montserrat/montserrat.woff2" as="font">
    <link href= "https://fonts.fontstorage.com/import/montserrat.css">

<style>
    
* {
  box-sizing: border-box;

}

figure {
    width: 47%; /* Ширина если надо расположить 2 картинки в ряд*/
    float: left; /* Выстраиваем элементы по горизонтали */
    margin: 0 0 0 0%; /* Отступ слева */
   text-indent: 0px; /* убираем отступ для картинки как ни странно */
   /* background: #f0f0f0; /* Цвет фона */
    border-radius: 1px; /* Радиус скругления */
    padding: 1%; /* Поля */
   }
   figure:first-child {
    margin-left: 0; /* Убираем отступ для первого элемента */
   }

picture {
 width: 30%; /* Ширина если надо расположить 3 картинки в ряд*/
    float: left; /* Выстраиваем элементы по горизонтали */
    margin: 0 0 0 0%; /* Отступ слева */
   text-indent: 0px; /* убираем отступ для картинки как ни странно */
   /* background: #f0f0f0; /* Цвет фона */
    border-radius: 1px; /* Радиус скругления */
    padding: 1%; /* Поля */
}
picture:first-child {
    margin-left: 0; /* Убираем отступ для первого элемента */
   }

   image {
    width: 100%; /* Ширина если надо расположить 1 картинки в ряд*/
    float: left; /* Выстраиваем элементы по горизонтали */
    margin: 0 0 0 0%; /* Отступ слева */
   text-indent: 0px; /* убираем отступ для картинки как ни странно */
   /* background: #f0f0f0; /* Цвет фона */
    border-radius: 1px; /* Радиус скругления */
    padding: 1%; /* Поля */
   }
</style>
</head>
<body>

<div class="sidenav">

    <a href="#h1">20. ДИАГНОСТИКА ТРАНСФОРМАЦИИ ЗЕМЛЕПОЛЬЗОВАНИЯ В СЕЛЬСКОЙ МЕСТНОСТИ С ИСПОЛЬЗОВАНИЕМ ГИС-МОДЕЛИРОВАНИЯ</a>
  

    <a href="#h3_1">20.1. Трансформация землепользования в сельской местности и ГИС-моделирование</a>
    <a href="#h3_2">20.2. Гипотеза и состав модели для исследования процессов экореабилитации</a>
    <a href="#h3_3">20.3. Моделирование степени зарастания сельскохозяйственных угодий</a>
    <a href="#h3_4">20.4. Оценка вклада основных природных и социально-географических факторов в трансформацию землепользования</a>
    <a href="#h3_5">20.5. Оценка изменения функций и типологии сельских поселений в связи с трансформацией землепользования</a>
    
</div>

<div class="content">
  
  <h1 id="h1">V. ДИАГНОСТИКА СЕЛЬСКОЙ МЕСТНОСТИ С ИСПОЛЬЗОВАНИЕМ ГИС-МОДЕЛИРОВАНИЯ</h1>     

  <h2 id="h2">20. ДИАГНОСТИКА ТРАНСФОРМАЦИИ ЗЕМЛЕПОЛЬЗОВАНИЯ В СЕЛЬСКОЙ МЕСТНОСТИ С ИСПОЛЬЗОВАНИЕМ ГИС-МОДЕЛИРОВАНИЯ</h2> 

<h3 id="h3_1">20.1. Трансформация землепользования в сельской местности и ГИС-моделирование</h3>

<p>Трансформация землепользования и земельного покрова (<span class="monospace">land use/land cover change</span>) - одно из важнейших следствий изменения характера взаимодействия природы и общества. В данном разделе демонстрируются возможности исследования причинно-следственной связи трансформации землепользования средствами ГИС-моделирования  территории российского Нечерноземья в границах Смоленской области, испытывающей наряду с обычными факторами, определяющими современные тренды изменения системы сельского расселения, дополнительное воздействие приграничного (с республикой Беларусь) положения.</p> 

<p>Для территории Нечерноземья в целом характерны концентрация сельскохозяйственного производства вокруг более-менее крупных населенных пунктов, с одновременной деградацией периферийной части системы расселения, сопровождающейся запустением и забросом освоенных ранее сельскохозяйственных земель. Одним из наиболее достоверных (и независимых от официальной статистики способов) мониторинга такого рода процессов является оценка зарастания сельскохозяйственных угодий на основе временных серий летних спутниковых изображений (<span class="monospace">Landsat</span>  или <span class="monospace">Sentinel</span>) с определением вероятности перехода одного типа землепользования и/или сукцессии в другой.</p>

<p>В целом ряде работ (см., например,  <a href="BIBLIO.html#Прищепов и др., 2013">[Прищепов и др., 2013]</a>)  доказано, что практически на всей территории Российского Нечерноземья площадь открытых поверхностей без растительности монотонно сокращается, при этом темпы и характер экореабилитации определяются набором факторов, в первую очередь, мезо- и микроклиматом, влажностью почв, локальным рельефом местности и породным составом ближайших лесов. Не последнюю роль играет и характер предшествующего землепользования.</p>

<p>Само по себе использование подобных подходов сделалось возможным и популярным в значительной степени благодаря широкому распространению и постоянному уточнению специально созданной и поддерживаемой модели глобальных данных <span class="monospace">land use/land cover (LULC)</span> и прежде всего так называемой Глобальной карты лесного покрова университета Мериленда (Global Forest Change) <a href="BIBLIO.html#Potapov_et_al., 2015">[Potapov_et_al., 2015]</a>. Весьма популярным в такого рода моделировании является использование различных ландшафтных индексов, в частности <span class="monospace">NDVI</span>, а также сравнение результатов полученных по ДДЗ со сведениями публикуемых в официальных статистических сводках - Докладах о состоянии и использовании земель.</p>

<p>Трансформация землепользования - "пограничная тема", имеющая проекцию в области физической географии, экологии и экономической географии: это обстоятельство должно приниматься во внимание при выстраивании предварительной гипотезы любого моделирования в данной сфере. Так исследованиями эконом-географов доказано <a href="BIBLIO.html#Люри_et_al._2010">[Люри и др., 2010</a>; <a href="BIBLIO.html#Казьмин_2016">Казьмин, 2016]</a>, что в пределах Европейской России происходит масштабный процесс концентрации сельскохозяйственных функций, приводящий, с одной стороны, к сокращению ареалов экономической активности с последующей интенсификацией землепользования в границах таких ареалов, с другой стороны -  к утрате функций даже экстенсивного сельского хозяйства на  периферии с забросом и запустением масштабных площадей бывших сельхозугодий. Таким образом, географическая локация и позиция на шкале "Центр-Периферия" - как в макрорегиональном, так и региональном и локальных масштабах приобретают особое значение.</p>

<p>Значительное сокращение общей площади сельхозугодий отмечено и для принятой здесь в качестве сцены моделирования Смоленской области: с начала 2000-х гг. область ежегодно теряет от 0,1 до 1,5 тыс. га., при этом активно происходят процессы вторичного заболачивания почв: площади заболоченных сенокосов и пастбищ составляют около 200 тыс. га, пашни — более 160 тыс. Проблемность региона обусловлена его макроэкономическим положением  между российским Центром и развитой в сельскохозяйственном отношении республикой Беларусь, что является причиной агро-рекреационной трансформации всей системы сельского расселения ("дачное" давление Московского региона) и подавления производственной функции в результате "неподъемной" (для региона) рыночной конкуренции с соседней страной. Дополнительное усложнение всей картины землепользования вносит наличие на территории крупной ООПТ - национального парка Смоленское Поозерье, основная функция которого - сохранение исторических культурных и природных ландшафтов - также не способствует развитию современного сельского хозяйства.</p>

<p>Исследования социально-географических причин трансформации земель на территории Смоленской области  позволили выявить сокращение площади сельскохозяйственных угодий, сжатие и стягивание к районным центрам системы расселения, а также несоответствие характера землепользования в приграничных районах Смоленской области и Беларуси. Показано, что доля земель, которые были выведены из сельскохозяйственного оборота в 1990—2000 годы и сейчас заняты кустарниковой растительностью и мелколесьем, существенно варьирует по районам Смоленской области, вследствие особенностей территориальной структуры хозяйства, плодородия почв, удаленностью от центров первичной переработки и потребления продукции сельского хозяйства и пищевой промышленности <a href="BIBLIO.html#Барановский_2017">[Барановский, 2017]</a>.</p>


<p>Еще одна особенность территории Смоленской области как сцены моделирования - свойственное ей высокое ландшафтное разнообразие <a href="BIBLIO.html#Шкаликов_2004">[Шкаликов, 2004]</a>, определяемое свежестью всего комплекса гляциального, флювиогляциального, лимно-гляциального и лимно-аллювиального рельефа, а также положением на экотоне двух природно-ландшафтных зон - южной тайги и хвойно-широколиственных лесов. Это обстоятельство, в свою очередь, сказывается на разнообразии локальных условий поддержания эффективности сельского хозяйственного производства и развития процессов экореабилитации.</p>

<p>С учетом охарактеризованной выше специфики региона геоинформационное моделирование может быть полезным для ответа на следующие вопросы:</p>
<ol>
<li>В какой мере общедоступные глобальные данные могут использоваться для изучения динамики землепользования на региональном уровне и насколько актуальной является  параллельная оценка их валидности и возможности корректировки?</li>
<li>Каков реальный вклад основных природных и социально-географических факторов в трансформацию землепользования, и в какой степени этот вклад может быть оценен средствами ГИС-моделирования? </li>
<li>Можно ли в рамках ГИС-модели получить корректную картину пространственной дифференциации процессов трансформация землепользования по территории региона используя одну из общедоступных ОТЕ - например, матрицу сельских поселений Смоленской области?</li>
</ol>
<br>


<h3 id="h3_2">20.2. Гипотеза и состав модели для исследования процессов экореабилитации</h3>

<p>По природным условиям Смоленская область типична для регионов Центральной России. Большая ее часть занята Смоленско-Московской возвышенностью, на северо-западе и юго-востоке – низменные участки. Около <span class="greencursiv">40%</span> территории области занимают леса, преимущественно из мелколиственных и хвойных пород, почвы со средним плодородием. В структуре сельского хозяйства по стоимости продукции животноводство незначительно превышает растениеводство. По производству большинства зерновых, зернобобовых и других культур область находится во второй половине рейтинга российских регионов. По целому ряду признаков можно назвать Смоленскую область типичным регионом северо-запада России в части основных тенденций социально-экономического развития сельских территорий.</p>

<p>С учетом вышеприведенных обстоятельств структура геоинформационной модели (не претендующая на комплексную оценку трансформации сельского землепользования) может быть выстроена следующим образом. Зарастание полей - одно из проявлений процессов экореабилитации и изменения сельской системы расселения, которые контролируются сложным набором факторов. Тогда в качестве гипотезы можно принять что существуют разные причины, связанные с размещением конкретного сельскохозяйственного угодья относительно населенных пунктов, транспортных коммуникаций, лесных опушек и т.д., и эти особенности размещения можно считать <span class="bolditalic">детерминантами (объясняющими факторами) локального уровня</span>. Кроме того, на локальном уровне зарастание может контролироваться чисто <span class="bolditalic">физико-географическими факторами</span> - положением конкретного поля на форме мезорельефа, близостью к эрозионной сети, характером почвенного покрова. К <span class="bolditalic">детерминантам мезоуровня</span> можно отнести размещение территории сельского поселения относительно муниципального центра и крупных транспортных магистралей. Наконец, <span class="bolditalic">детерминантами макроуровня</span> (в региональном масштабе) могут считаться положение сельского района относительно Смоленска, границы с республикой Беларусь и западной периферии Московской агломерации. Согласно нашей гипотезе <span class="bolditalic">существует некое сочетание указанных факторов, приводящее к разному уровню развития процессов экореабилитации: (завершенному, значительному, слабому, невыраженному), индикаторами которых в данном модели выступает степень зарастания конкретных сельскохозяйственных угодий</span> <a href="BIBLIO.html#Колб_и_др__2018">[Колбовский, Климанова, Бавшин, 2018]</a>.</p>
<br>

<a id="00_Схема_модели"><img src="Pict_5_20/00_Схема_модели.png" width="100%"></a><br>
<br><span class="imgtitle">Рис. 20.1 Общая схема гипотезы моделирования 

<p>В рамках данной модели в качестве исследуемого феномена - характер экореабилитации бывших сельскохозяйственных угодий (СУ), выраженный через степень зарастания; в роли предикторов, в той или иной степени влияющих на процесс экореабилитации, рассматриваются три группы факторов:</p>
<ul>
<li>природные факторы, определяющие качество СУ (геоморфометрические параметры, типы ландшафтных местоположений, типы почв);</li>
<li>факторы положения СУ относительно природных объектов (близость к опушкам лесных массивов и древесно-кустарниковых групп,  болотам, озерам и водохранилищам, рекам);</li>
<li>факторы размещения, определяющие доступность СУ (близость к населенным пунктам разного типа и людности, близость к дорогам).</li>
</ul>

<p>Очевидно, что процессы экореабилитации могут  оцениваться в рамках различных операционно-территориальных единиц как природного (например, речные бассейны), так и социального (единицы АТД) порядка. По гипотезе, принятой в данном исследовании, <span class="bolditalic">трансформация землепользования прежде всего ощущается на локальном и среднемасштабном уровнях, приводя к изменению функциональной типологии сельских поселений,  а затем, – в региональном масштабе, вызывая общее перераспределение нагрузки между ареалами эффективной экономической деятельности и периферией, несущей средостабилизирующие (природоохранные) и рекреационные функции</span>.</p>
</p>

<p>Потенциальный характер взаимосвязи между исследуемым феноменом (степень экореабилитации выбывших из оборота сельскохозяйственных угодий) и предикторами определил необходимость развертывания трех самостоятельных блоков геоинформационной модели:</p>

<ol>
<li>Исследование особенностей применения общедоступных глобальных данных [14] для изучения динамики землепользования на региональном уровне с параллельной оценкой их валидности и возможности корректировки;</li>
<li>Оценка вклада основных природных и социально-географических факторов в трансформацию землепользования посредством построения каузальной ГИС-модели;</li>
<li>Обобщение закономерностей пространственной дифференциации процессов трансформация землепользования по территории региона в рамках матрицы сельских поселений Смоленской области.</li>
</ol>

<p>В соответствии с содержанием каждого блока моделирования строится специальная комплексная ГИС-модель сельскохозяйственных угодий, включающая около 20 тематических растров, принадлежащих к трем группам:  растры геоморфометрических параметров, ландшафтных местоположений, растры характеризующие отношения близости и пространственные закономерности освоения. Для всех контуров сельскохозяйственных угодий, привлеченных из официальных источников (сайта Министерства Сельского хозяйства РФ) необходимо предварительно получить показатели фактического зарастания. Далее показатели экореабилитации и объясняющих факторов извлекаются в сетку операционно-территориальных единиц модели  - сельских поселений Смоленской области.</p>
<br>


<h3 id="h3_3">20.3. Моделирование степени зарастания сельскохозяйственных угодий</h3>

<p>Выявление степени зарастания СУ по данным официальной статистики на сегодняшний день затруднено отсутствием надежной базы (векторных слоев) сельскохозяйственных угодий  даже для центральных областей Российской Федерации.</p>

<p>Для данной модели использовались контура сельскохозяйственных угодий <a href="http://atlas.mcx.ru/">Атласа земель сельскохозяйственного назначения</a> по Смоленской области и растр лесного покрова <a href="http://glad.geog.umd.edu/dataset/eastern-europe-forset-cover-dynamics-1985-2012/">Eastern Europe forest cover dynamics</a>. В унифицированной легенде <span class="monospace">Eastern Europe forest cover dynamics</span>, выделяется семь типов ареалов: 1 – стабильные нелесные земли; 2 – стабильные лесные земли; 3 –  лесные земли, сформированные на ранее (по состоянию на 1985 г.) нелесных землях; 4 – утраты лесов; 5 – утраты лесов с последующим ростом леса; 6 – участки с чередованием утрат и прироста леса;  7  - утраты лесов на участках, которые приобрели лесной покров после нелесного состояния в 1985 г.. Построенные на основе методики выявления изменений (<span class="monospace">Change Detection</span>), пять из семи данных категорий отражают последствия трансформации в бинарной классификации – лес- не лес, не различая причин подобных изменений.</p>

<p>С помощью инструмента <span class="blackbold">Зональная Статистика</span> <span class="red">ArcMAP10.x</span> или утилиты <span class="blackbold">Add Grid Values to Shapes</span> <span class="blue">SAGA GIS</span> определяется доля, занимаемая древесным покровом  в пределах конкретного поля. Следует иметь ввиду, что в соответствии с правилами создания слоев <span class="monospace">forest cover</span> в данной процедуре невозможно учесть самые начальные фазы зарастания угодья кустарниками или подростом мелколиственных и хвойных пород.</p>
 
<p>Далее в зависимости от доли (%), занимаемой древесным покровом вся совокупность контуров угодий была ранжирована  по степени зарастания (способом <span class="cursive">Естественные границы|Natural break</span>) на 10 классов. В результате этого сопоставления была получена карта ареалов трансформации землепользования, легенда которой включает имеющиеся в настоящий момент лесные массивы разной степени трансформации (контура взяты с карты <span class="monospace">Global Forest Change</span>), селитебные земли (контура взяты с карты <span class="monospace">Open Street Map</span>) и сельскохозяйственные угодья разной степени трансформации – от обрабатываемых "чистых" полей до полей в начальной стадии зарастания  и практически полностью заросших.</p>


<a id="СХ_зарастание_Ранг_Цвет"><img src="Pict_5_20/СХ_зарастание_Ранг_Цвет.png" width="55%" ></a><br>
<br><span class="imgtitle">Рис. 20.2  Ранжирование сельскохозяйственных угодий по степени зарастания от 1 до 10 (1 – поля, возделываемые в настоящее время; 10 – полностью заросшие поля)


<p>Сопоставление этих данных с результатами полевых исследований, проведенных в  2012-15 гг. в Смоленском и Демидовском районах области, показало, что факторы, вызвавшие подобную трансформацию, могут быть разбиты на четыре группы:</p>

<ul>
<li>демутационные сукцессии, вызванные прежде всего забросом и зарастанием сельскохозяйственных угодий, экореабилитацией мелиорированных болотных массивов, а также постпирогенные и пострубочные сукцессии;</li>
<li>антропогенные воздействия (рубки, прокладывание инженерных линеаментов, строительство дорог);</li>
<li>экстремальные природные явления (для Смоленской области это, прежде всего, обширные ветровалы в зрелых и перестойных лесных насаждениях),</li>
<li>результаты естественной многолетней динамики ландшафтов (изменение типа и характера питания болотных массивов, эвтрофикация дистрофных озер, зарастание вновь образованны пойменных грив и иные)</li>
</ul>

<p>Это обстоятельство не позволяет напрямую использовать глобальные данные для верифицированных оценок трансформации землепользования на региональном уровне и построения причинно-следственных моделей с возможностью последующего прогноза. С другой стороны, чрезвычайно высокая трудоемкость получения результатов дешифрирования на обширные территории делает невозможным повторение исходного алгоритма в условиях ограниченных средств и вычислительных возможностей. В подобных случаях можно попытаться  диверсифицировать стандартные категории <span class="monospace">change detection</span> на ключевых ареалах экспертно (вручную участок за участком) привлекая крупнодетальные космические снимков (например "подгрузить" вектора ключевых ареалов в программа SASPLANETA). Результат можно верифицировать дополнительно расчетом индексом <span class="monospace">NDVI</span> по снимкам <span class="monospace">Landsat</span> с учетом региональных особенностей вегетации и проведения основных стадий сельскохозяйственных работ.</p>

<p>Разумеется, оптимальное решение - проведение полевых исследований на ключевых участках. Для данной модели такие исследования были проведены   период с 2012 по 2016 гг. на территории двух муниципальных районов с резко различными природными и социально-экономическими условиями Смоленского (прицентрального) и Демидовского (периферийного с высокой долей  природных экосистем и заброшенных сельскохозяйственных угодий)</p>

<p>Результаты, отраженные в <a href="#Таблица 20.1">Таблице 20.1</a>, свидетельствуют о значительной диверсификации  действительного содержания одинаковых стандартных глобальных индексов по четырем группам факторов. Так, например, категория <span class="greencursiv">4</span> – утраты лесов, соответствует как непосредственным техногенным воздействиям (застройка, устройство карьеров, прокладка линеаментов), так и последствиям экстремальных явлений (гари, ветровалы).  Категория <span class="greencursiv">3</span> (лесные земли, сформированные на ранее по состоянию на 1985 г. нелесных землях) попадает в ареалы разнотипных демутационных сукцессий  - зарастание полей разных  ландшафтных местоположений, а также пойменных лугов, ранее использовавшихся под сенокосы. Таким образом, выявлено, что для внутрирайонного уровня использование данных дистанционного зондирования позволяет достаточно достоверно определить тренды, но не позволяет установить причины и ландшафтно-экологические последствия процессов трансформации.</p>

<a id="Таблица 20.1"><span class="imgtitle">Таблица 20.1   Диверсификация типов трансформации земель на территории Смоленского и Демидовского районов Смоленской области по данным мануального дешифрирования и полевых исследований

<div class="table">
<table id="customers">
        <table border="1">

<tr>
  <th>№</th>
  <th>Тип  трансформации</th>
  <th colspan="2">Доля от общей площади земель района, испытавших трансформацию в %</th>
  <th>Тип (индексы) трансформации forest change</th>
</tr>
<tr>
<td></td>
<td></td>
<td>Демидовский</td>
<td>Смоленский</td>
<td></td>
</tr>
<tr>
<td>1</td>
<td colspan="4">Антропогенные воздействия</td>
</tr>
<tr>
<td>1.1.</td>
<td>застройка</td>
<td>0,5</td>
<td>2,43</td>
<td>4, 1, 2, 5</td>
</tr>
<tr>
<td>1.2.</td>
<td>карьеры</td>
<td>0,5</td>
<td>0,7</td>
<td>4</td>
</tr>
<tr>
<td>1.3.</td>
<td>прокладка линеаментов</td>
<td>0,05</td>
<td>1,67</td>
<td>4</td>
</tr>
<tr>
<td>1.4.</td>
<td>рубки</td>
<td>2,11</td>
<td>0,58</td>
<td>4</td>
</tr>
<tr>
<td>1.5.</td>
<td>гари</td>
<td>0,31</td>
<td>-</td>
<td>4</td>
</tr>
<tr>
<td></td>
<td>Итого</td>
<td>3,47</td>
<td>5,38</td>
<td>1, 2, 4, 5</td>
</tr>
<tr>
<td>2</td>
<td colspan="4">Демутационные сукцессии</td>
</tr>
<tr>
<td>2.1.</td>
<td>зарастание селитебных земель</td>
<td></td>
<td></td>
<td>3, 5</td>
</tr>
<tr>
<td>2.2.</td>
<td>зарастание нарушенных территорий</td>
<td>0,06</td>
<td>0,36 </td>
<td>3, 5</td>
</tr>
<tr>
<td>2.3.</td>
<td>зарастание полевых контуров на границе с лесом</td>
<td>23,48</td>
<td>41,22</td>
<td>3</td>
</tr>
<tr>
<td>2.4.</td>
<td>зарастание полей среди леса</td>
<td>14,51</td>
<td>2,55</td>
<td>3</td>
</tr>
<tr>
<td>2.5.</td>
<td>зарастание  полей в массивах сельхозземель</td>
<td>36,96</td>
<td>36,64</td>
<td>3</td>
</tr>
<tr>
<td>2.6.</td>
<td>зарастание лугов речных пойм</td>
<td>3,85</td>
<td>5,63</td>
<td>3, 6</td>
</tr>
<tr>
<td>2.7.</td>
<td>восстановление естественной растительности в местах рубок и ветровалов</td>
<td>10,96</td>
<td>6,36</td>
<td>3</td>
</tr>
<tr>
<td></td>
<td>Итого</td>
<td>89,82</td>
<td>92,76</td>
<td>3, 5, 6</td>
</tr>
<tr>
<td>3</td>
<td colspan="4">Экстремальные природные явления</td>
</tr>
<tr>
<td>3.1.</td>
<td>образование ареалов ветровалов</td>
<td>0,1</td>
<td>-</td>
<td>4, 5</td>
</tr>
<tr>
<td>4</td>
<td colspan="4">Среднемноголетняя динамика ландшафтов</td>
</tr>
<tr>
<td>4.1</td>
<td>развитие сообществ водных растений в дистрофных озерах</td>
<td>0,18</td>
<td>-</td>
<td>3</td>
</tr>
<td>4.2.</td>
<td>зарастание верховых болот</td>
<td>4,55</td>
<td>-</td>
<td>3</td>
</tr>
<tr>
<td></td>
<td>Итого</td>
<td>4,73</td>
<td>-</td>
<td>3</td>
</tr>

</table>

*Индексы глобальных данных forest change: 1 - stable non-forest; 2 stable forest; 3 - forest gain over non-forest in 1985; 4 - forest loss; 5 - forest loss followed by forest gain; 6 -  repeated forest loss separated by forest gain; 7  - forest loss on areas which gain forest cover after non-forest state in 1985

<p>При этом значительными являются различия в трансформации земель между максимально урбанизированным (относительно условий Смоленской области) и максимально сохранившим природные экосистемы и экстенсивное сельское хозяйство периферийным районом. Для Смоленского района свойственна высокая доля трансформации землепользования, связанная с застройкой, прокладкой дорог и инженерных линеаментов, при этом отсутствуют типы трансформации связанные с естественной динамикой ландшафтов, характерные для заболоченного и заозеренного Демидовского района. Однако (и это можно считать особенностью региона) прицентральное положение Смоленского района практически не сказывается на снижении интенсивности зарастания сельхозугодий.</p>

<p>Для решения второй задачи – корректировки метода выявления трендов трансформации землепользования  -   спутниковые изображения <span class="monospace">Landsat 4,5,7,8</span> за период с 1984 по 2016 годы обрабатывались для ареалов квадратной формы размерностью <span class="greencursiv">100*100 км</span> с предварительным выбором алгоритма расчета индекса <span class="monospace">NDVI</span>, который учитывал бы особенности регионального цикла сельскохозяйственных работ (для полевых угодий) и разработкой классификации демутационных сукцессии (для лесопокрытых земель) <a href="BIBLIO.html#Колб_и_др__2018">[Колбовский и др., 2018]</a>.</p>

<p>Полученные тематические карты, отображающие стадии восстановительных сукцессий  на бывших сельскохозяйственных угодьях, показали, что в течении всего периода наблюдения площади открытых участков без растительности монотонно сокращались.  Обнаружено, что данный процесс несет выраженный унаследованный  (от позднесоветского и раннего перестроечного периода) характер: так к 1990 г. около <span class="greencursiv">38%</span> нелесных территорий уже были покрыты древесной растительностью. За следующие 25 лет прирост площади лесов (в основном за счет лиственных пород) составил около <span class="greencursiv">18%</span>, из которых половина приходится на 1990-е годы. В 2015 году <span class="greencursiv">46,0%</span> бывших нелесных территорий было занято березовыми и <span class="greencursiv">9,6%</span> - сосновыми молодняками.</p>

<p>Отдельный интерес в моделях трансформации землепользования и ландшафтов (<span class="monospace">LULC</span>) представляет сравнительный анализ выявляемых изменений с данными лесоустройства по территориям лесного фонда региона. В данном случае оценка точности карты растительности и динамики изменения лесного покрова  проводились  проводилась по двум параметрам: по покрытой лесом площади и по группам пород (по хвойным и лиственным лесам).</p>

<p>Предварительно, верификации границ лесного фонда, использованных в оценке точности, проводится анализ официальных статистических показателей площади лесного фонда с площадями ареалов векторных карт лесного фонда ГУЛФ. В результате сравнения площади лесного фонда по данным ГУЛФ 2008 года (<span class="greencursiv">896700 га</span>) и по цифровой карте лесного фонда 2010 года (<span class="greencursiv">1954160 га</span>) было выявлено значительное отклонение в оценках площадей Смоленской области (<span class="greencursiv">54,1%</span>), что по сути является еще одним подтверждением масштабного процесса экореабилитации увеличения  общей лесопокрытой площади Смоленской области.</p>

<p>Сравнение площадей лесов, полученных по данным дистанционного зондирования (Landsat), растра <a href="https://glad.earthengine.app/view/global-forest-change#dl=1;old=off;bl=off;lon=100.5078125;lat=-13.209155847813578;zoom=3;">Global Forest Change</a>  Университета Мериленда, с данными статистической отчетности показало хорошую сходимость площадей лесов на уровне субъектов (среднеквадратическое отклонение <span class="greencursiv">3,4%</span>).</p>

<p>Третья задача данного блока модели - определение особенностей экореабилитации на исследуемой территории Смоленской области. Для моделирования изменений лесного покрова был использован ключевой участок – территория национального парка Смоленское Поозерье, для которого был создан пространственно-временной набор данных, отражающий динамику лесного покрова. Чтобы преодолеть несоответствие понятий и методов оценок лесохозяйственных характеристик, использовался архив данных космической съемки. Лесной покров в анализе определялся набором биофизических признаков - <span class="cursive">сомкнутость крон деревьев</span>  и <span class="cursive">высота насаждений</span>, параметры дешифрирования, которые (в силу понятных причин) не всегда соответствует данным лесоустройства.</p>

<p>Дополнительно с использованием редактора JOSM проекта OpenStreetMap выполнялась верификация нелесных территорий национального парка. В качестве подложки использовались крупномасштабные снимки от <span class="monospace">DigitalGlobe, Mapbox satellite </span> и др., с разрешением <span class="greencursiv">1 м</span> и лучше, дата съемки 2010 -2016 гг. Сотрудником НП Смоленское Поозерье <b>И.М.Бавшиным </b>было верифицировано более <span class="greencursiv">300</span> полигонов болот, <span class="greencursiv">2337</span>  полигонов необрабатываемых полей, <span class="greencursiv">860</span> полигонов обрабатываемых полей и <span class="greencursiv">704</span> полигона обрабатываемых земельных участков в границах населенных пунктов (огороды, приусадебные хозяйства).</p>

<p>Сравнительный анализ позволил выявить значительные расхождения между традиционными материалами лесоустройства ФГБУ Национальный парк Смоленское Поозерье по состоянию на 2015 год и атрибутированными тематическими растрами, полученными по данным Landsat. Так, при сравнении площадей, занятых различными преобладающими породами процент совпадений составил: для сосновых лесов <span class="greencursiv">82,7%</span>, еловых <span class="greencursiv">76,4%</span>, березовых <span class="greencursiv">77,6%</span>, осиновых <span class="greencursiv"><23,2%</span> сероольшанников <span class="greencursiv">71,8%</span>, черноольшанников <span class="greencursiv">50,4%</span>. Аналогичные показатели получены и при сравнении соотношения лиственных и хвойных пород в лесонасаждениях; совпадение составило: хвойные  <span class="greencursiv">53,7%</span>,  смешанные, с преобладанием хвойных <span class="greencursiv">43,8%</span>, смешанные  <span class="greencursiv">41,1%</span>, смешанные, с преобладанием лиственных <span class="greencursiv">19,4%</span>, лиственные  <span class="greencursiv">88,2%</span>. Столь существенное расхождение серьезно затрудняет совместную интерпретацию результатов состояния растительности по данным ДДЗ и лесотаксации, вынуждая выбрать первые в качестве более надежного источника. </p>
<br>



<h3 id="h3_4">20.4. Оценка вклада основных природных и социально-географических факторов  в трансформацию землепользования</h3>


<p>Оценка вклада основных природных и социально-географических факторов в  трансформацию землепользования осуществлялась в рамках блока ГИС-модели, включающего общедоступные слои <span class="monospace">OpenStreetMap</span> на Смоленскую область, контура сельскохозяйственных угодий (СУ), а также специальные тематические (пофакторные) слои, принадлежащие к трем группам факторов:</p>

<ol>
<li>природные факторы, определяющие качество СУ - геоморфометрические параметры, типы ландшафтных местоположений, типы почв;</li>
<li>факторы положения СУ относительно природных объектов - близость к опушкам лесных массивов и древесно-кустарниковых групп, болотам, рекам, озерам и водохранилищам; </li>
<li>факторы локализации, определяющие доступность СУ - близость к населенным пунктам разного типа и людности, близость к дорогам.</li>
</ol>

<p>На начальном этапе необходимо получить нормализованные атрибутированные растры основных геоморфометрических параметров:</p> 

<div class="script"> 

Normalised Height|Нормализованная высота, <br>
Height Above River|Высота над уровнем рек, <br>
Slope|Уклон,<br>
Plain Curvature|Плановая кривизна,<br>
Profile Curvature|Профильная кривизна,<br>
Aspect|Экспозиция, <br>
TPI|Индекс топографической позиции
</div>

<p>Далее значения тематических растров извлекаются в векторные контура сельскохозяйственных угодий инструментом <span class="blackbold">Add Grid Values to Shapes</span> <span class="blue">SAGA GIS</span>.</p>

<p>На основе полученных тематических растров с помощью усовершенствованного алгоритма <span class="blackbold">Fuzzy Overlay /span>  с добавлением переменной <span class="monospace"></span>TPI третьей размерности и растра <span class="monospace">Height Above River|Высота над уровнем рек</span> строится карта ландшафтных местоположений.</p>

<a id="00_Landscape_Siites_Black"><img src="Pict_5_20/00_Landscape_Siites_Black.png" width="55%"></a><br>
<br><span class="imgtitle">Рис. 20.3 Векторная карта ландшафтных местоположений – модель, построенная на основе набора геоморфометрических переменных по ЦРМ масштаба 30 м в пикселе

<p>Моделирование природных факторов позволило провести кластеризацию угодий по всему набору условий локализации полей в пределах конкретной мезоформы рельефа, в пределах определенного ландшафтного местоположения и типа почв. Кластерный анализ выявил семь отчетливо выделяющихся классов, каждый из которых характеризуется определенным высотным уровнем, локализацией на катене гребне-килевой дифференциации и приуроченностью к рассеивающим\концентрирующим, склонам различной крутизны, а также позволил констатировать чрезвычайное разнообразие и высокую дисперсность распределения полей разных типов по территории Смоленской области. </p>

<div class="script_01">
<figure> 
<a id="Кластер_Геоморфа"><img src="Pict_5_20/Кластер_Геоморфа.png" width="100%" ></a>
</figure>
<figure> 
<a id="Диаграмма_Геоморфа"><img src="Pict_5_20/Диаграмма_Геоморфа.png" width="61%" ></a>
</figure>
</div>
<br>
<span class="imgtitle">Рис. 20.4 a)  Классы сельскохозяйственных угодий Смоленской области по геоморфометрическим параметрам b) Диаграмма доли, занимаемой каждым классом

<p>В структуре сельхозугодий Смоленской области преобладают поля на невысоких и низких расположенных мезоформах рельефа, расположенных в относительной близости к эрозионной сети, пологие и слабопокатые (<span class="bolditalic">класс 6</span>), а также на средних частях некрутых склонов (<span class="bolditalic">класс 1</span>). Склонами значительной крутизны от этих двух типов отличается <span class="bolditalic">класс 5</span>. <span class="bolditalic">Класс 3</span> – поля на привершинных поверхностях, близок к нем <span class="bolditalic">класс 4</span>, приуроченный к  верхним крутым частям склонов сводов возвышенностей. Малочисленный <span class="bolditalic">класс 2</span> – невысокие водоразделы, слабодренированные и плоские. <span class="bolditalic">Класс 7</span> –  поля на невысоких, но хорошо расчлененных  поверхностях (эродированных вторичных моренных равнинах).</p>

<p>Моделирование положения сельскохозяйственных угодий по отношению к природным объектам отражает один из аспектов потенциала экореабилитации, поскольку очевидно, что поля граничащие с лесными массивами имеют больше шансов подвергнуться зарастанию, а поля, расположенные недалеко от болот и когда-то подвергнутые мелиорации посредством дренажных систем – вторичному заболачиванию и закочкариванию. Стандартный набор слоев <span class="monospace">OSM</span> позволяет провести расчет параметров близости (инструмент <span class="blackbold">Near</span> <span class="red">ArcMAP10.x</span>) по отношению к 10 видам природных объектов, которые предварительно должны быть извлечены выборкой в отдельные (новые) векторные слои:</p>

<div class="script">

леса (vegetation-polygon, "NATURAL" = 'wood')<br>
кустарники и поросль (кустарники и поросль, "NATURAL" = 'scrub')<br>
болота (water-polygon, "NATURAL" = 'wetland')<br>
крупные реки (water-polygon, "NATURAL" = 'water')<br>
средние реки (water-line, "WATERWAY" = 'river')<br>
малые реки (water-line, "WATERWAY" = 'stream', Lenght>=10 км)<br>
ручьи (water-line, "WATERWAY" = 'stream', Lenght<=10 км <br>
озера (water-polygon, "WATERWAY" = 'lake')<br>
водохранилища (water-polygon, "WATERWAY" 'reservoir')<br>
каналы и дрены (water-line, "WATERWAY" = 'drain' OR "WATERWAY" = 'canal')<br>
</div>

<p>Значения близости для каждого контура поля по отношению к набору факторов позволяют провести классификацию сельскохозяйственных угодий: <span class="blackbold">Анализ группирования</span> <span class="red">ArcMAP10.x</span> определил оптимально число "групп" - <span class="greencursiv">7</span>; кластеризация СУ действительно выявляет хорошо интерпретируемые классы,   различающиеся по условиям локализации.</p>
 
<a id="Кластеры_Близость_Природы"><img src="Pict_5_20/Кластеры_Близость_Природы.png" width="60%"></a><br>
<br><span class="imgtitle">Рис. 20.5 Классы сельскохозяйственных угодий Смоленской области выделившиеся по близости к природным объектам

<p><span class="bolditalic">Класс 7</span> преобладает количественно среди прочих угодий и объединяет поля, расположенные в относительной близости от лесных массивов, болот и озер – условия характерные для дренированных небольшими водотоками (средними и малыми реками) заболоченных и заозеренных озерно-водноледниковых равнин, занимающих значительную площадь на территории Смоленской области. <span class="bolditalic">Класс 6</span> - это СУ с усредненными значениями всех параметров. <span class="bolditalic">Класс 2</span> фиксирует СУ, находящиеся в пределах фрагментированных ландшафтов с обилием древесно-кустарниковых групп, лесополос и ремизов, кроме того это местность с обилием озер. <span class="bolditalic">Класс 3</span> – также в пределах слабодренированных равнин в свое время мелиорированных системой осушительных каналов и дрен, с близким расположением лесов и болот. Наконец, четко выделяется <span class="bolditalic">класс 5</span> - СУ, находящиеся в хорошо дренированных местоположениях относительно далеко от лесных опушек но, в непосредственной близости от всех элементов речной сети (крупных и малых рек).</p>

<p>Следующий шаг - моделирование положения контуров полей по отношению к социально-экономическим факторам (населенным пунктам и дорогам) выявило отчетливую анизотропию современной картины распределения сельскохозяйственных угодий. Отдельный интерес представляет определение "силы влияния" параметра близости  моделируемых факторов. Для этой цели по гистограмме параметра близости фиксировались пороговые значения пяти классов (способ - Natural Breaks), далее определяется число объектов (т.е., полей) попадающих в первые два класса близости: логично предположить, что если таких объектов много, значит фактор имеет высокую значимость, и наоборот, если множество объектов расположены далеко от предполагаемого "центра притяжения" это может означать, что фактор не работает или работает как "отталкиватель": в частности, например, поля не могут располагаться в непосредственной близости от границ крупного города. Кроме того, графики близости позволяют  оценить характер распределения (нормальное, экспоненциальное) и знак связи.</p>

<div class="script_01">
<figure> 
<a id="Дист_Город"><img src="Pict_5_20/Дист_Город.png" width="99%"></a>
</figure>
<figure> 
<a id="Дист_Деревни"><img src="Pict_5_20/Дист_Деревни.png" width="100%"></a>
</figure>
</div>
<br>
<span class="imgtitle">Рис. 20.6 a)  Гистограмма  распределения полей по расстоянию от городов, b) Гистограмма  распределения полей по расстоянию от  деревень

<p>Оценка, проведенная отдельно для разных типов населенных пунктов -  городов, поселков городского типа, деревень, мелких хуторских пунктов и дачных поселков - выявила два резко различающихся типа распределения: условно городской и деревенский. Города Смоленской области формирует кольца освоенного пространства <a href="#Дист_Город">(Рис. 20.6 a)</a>, в пределах которых на ближних расстояниях предпочтение отдается иным видам землепользования (здесь расположены лесопарковые зоны лесничеств, рекреационные зоны, дачные кооперативы, промышленные предприятия, элементы транспортной инфраструктуры), поэтому наибольший процент СУ находится в полосе  приблизительно 15-30 км от городского центра (в целом распределение имеет нормальный характер) и далее кривая идет резко на убыль – это значит что  на расстоянии около 50 км всякие признаки сельскохозяйственного использования  земель исчезают. Совершенно  иное (логнормальное) распределение обнаруживается  для деревень,  где большинство СУ укладывается в полосу 500 м – 2,5 км, а за пределами 5 км зона освоенности уступает либо лесным массивам, либо зоне тяготения соседнего населенного пункта <a href="#Дист_Деревни">(Рис. 20.6 b)</a>.</p>

<p>Отдельный интерес в рамках данного блока общей модели представляет выявление гипотетической закономерности размещения СУ относительно сельских НП различной размерности. Для выявления возможной связи все сельские населенные пункты  были разделены по людности на 5 классов (пороговые значения - естественные границы):</p>

<span class="monospace">
1-й класс 0 - 9 человек,<br>  
2-й класс 10 - 30,<br> 
3-й класс 32 -  89,<br> 
4-й класс 90 -  241,<br> 
5-й класс 242 - 458 человек. </span>

<p>Затем определялся параметр близости полей к населенным пунктам различной людности.</p>

<a id="Таблица 20.2"><span class="imgtitle">Таблица 20.2 Расположение сельскохозяйственных угодий Смоленской области  относительно сельских населенных пунктов различной людности 


<div class="table">
<table id="customers">
        <table border="1" width=100%>


<tr>
<th rowspan="2">Классы людности</th>
<th colspan="5">Кластеры полей по близости к населенным пунктам (расстояние в м)</th>
</tr>
<tr>
<th>Ближайшие</th>
<th>Близкие</th>
<th>Средне удаленные</th>
<th>Удаленные</th>
<th>Самые дальние</th>
</tr>
<tr>
<td>1</td>
<td>659,8</td>
<td>1787,1</td>
<td>3147,6</td>
<td>5013,4</td>
<td>9911,2</td> 
</tr>
<tr>
<td>2</td>
<td>861,1</td>
<td>2232,8</td>
<td>3878,8</td>
<td>6190,9</td>
<td>9992,4</td> 
</tr>
<tr>
<td>3</td>
<td>979,3</td>
<td>2547,9</td>
<td>4273,3</td>
<td>6469,3</td>
<td>9998,5</td> 
</tr>
<tr>
<td>4</td>
<td>1097,0</td>
<td>2900,9</td>
<td>4861,7</td>
<td>7096,9</td>
<td>9998,9</td> 
</tr>
<tr>
<td>5</td>
<td>1018,4</td>
<td>2821,8</td>
<td>4742,9</td>
<td>7056,1</td>
<td>9957,3</td> 
</tr>

</table>


<p>В общем случае для территории Смоленской области можно констатировать что чем меньше населенный пункт, тем меньше расстояние до СУ, принадлежащих к одному и тому же классу близости, но эта тенденция работает лишь до населенных пунктов 4 класса людности.</p>

<a id="Таблица 20.3"><span class="imgtitle">Таблица 20.3 Доля сельскохозяйственных угодий (%), попадающих в разные классы близости относительно населенных пунктов различной людности 


<div class="table">
<table id="customers">
        <table border="1" width=100%>

<tr>
<th rowspan="2">Классы людности</th>
<th colspan="5">Доля сельхозугодий, попадающих в разные классы близости к НП</th>
</tr>
<tr>
<th>Ближайшие</th>
<th>Близкие</th>
<th>Средне удаленные</th>
<th>Удаленные</th>
<th>Самые дальние</th>
</tr>
<tr>
<td>1</td>
<td>43,3</td>
<td>25,6</td>
<td>18,8</td>
<td>10,1</td>
<td>3,3</td> 
</tr>
<tr>
<td>2</td>
<td>38,7</td>
<td>26,6</td>
<td>20,0</td>
<td>11,2</td>
<td>3,2</td> 
</tr>
<tr>
<td>3</td>
<td>34,7</td>
<td>24,2</td>
<td>21,4</td>
<td>13,7</td>
<td>6,0</td> 
</tr>
<tr>
<td>4</td>
<td>36,0</td>
<td>22,9</td>
<td>18,9</td>
<td>13,9</td>
<td>8,3</td> 
</tr>
<tr>
<td>5</td>
<td>48,0</td>
<td>17,1</td>
<td>15,7</td>
<td>11,7</td>
<td>7,5</td> 
</tr>
</table>

<p>Тенденция может быть лучше понята, если выявить долю полей, попадающих в классы близости к населенным пунктам различной людности: здесь зависимость несет более сложный характер. Процент полей в ближайшей группе выше для ненаселенных и самых населенных деревень, в деревнях средней людности этот процент уменьшается.</p>

<a id="Доля_СУ_Классы_близости"><img src="Pict_5_20/Доля_СУ_Классы_близости.png" width="70%"></a><br>
<br><span class="imgtitle">Рис. 20.7 Доля сельхозугодий попадающих в разные классы близости относительно населенных пунктов различной людности: ряд 1 - ближайшие поля, ряд 2 - близкие, ряд 3 - среднеудаленные, ряд 4 - удаленные, ряд 5 - самые дальние


<p>Таким образом, система расселения Смоленской области формирует свое силовое поле распределения сельскохозяйственных угодий с центрами-аттракторами различной силы и дистанции воздействия. В сочетании с геоморфологическими факторами и фактором расстояния до природных объектов общая мозаика агроэкосистем приобретает характер сложного феномена, исторически формировавшегося в многомерном пространстве факторов.</p>


<p>Выявление характера  зависимости процессов экореабилитации от трех групп факторов (ландшафтно-геоморфологических, переменных близости к природным объектам и переменных близости к населенным пунктам) проводилось в два этапа. На первом этапе внутри каждой группы с помощью <span class="blackbold">Explanatory regression|Исследовательской регрессии</span>  <span class="red">ArcMAP10.x</span>  верифицировалась достаточность взаимной независимости привлеченных для анализа факторов. В группе геоморфологических факторов все переменные (несмотря на то, что все они были получены на основе высотных данных SRTM) продемонстрировали достаточный уровень автономности (т.е., прошли критерий пространственной автокорелляции). Гипотеза автономности была подтверждена для  группы факторов социально-экономического характера, связанных с системой расселения. В группе факторов рисунка освоения ожидаемо связанными оказались факторы близости к различным звеньям речной сети, поэтому в итоговой модели все они были объединены в качестве признака расстояния от СУ до ближайшего (любого) водотока.</p>

<p>На следующем этапе в программном комплексе span <class="blue">SAGA GIS</span> проводился собственно регрессионный анализ  с использованием инструмента <span class="blackbold">Multiple Linear Regression Analysis (Shapes)</span>;   степень зарастания определялась как зависимая величина (<span class="cursive">`Dependent Variable</span>), а признаки трех групп факторов - как предикторы (<span class="cursive">predictors</span>).</p>  

<p>Полученные результаты позволяют констатировать выраженный стохастический  характер  процесса зарастания сельхозугодий. Проверка статистической значимости и силы зависимостей по группам выглядит следующим образом. Наименее выражены связи между экореабилитацией СУ и геоморфологическими факторами – здесь достоверная (отрицательная по знаку) связь получена только для  <span class="monospace">индекса Топографической позиции (TPI)</span>,  характеризующего расположение участка на склоне (верхняя, средняя, нижняя части) и на шкале гребне-килевой дифференциации (тальвеги - гребни). В общем случае можно утверждать, что <span class="bolditalic">участки расположенные ближе к эрозионной сети и в нижней части склона имеют большую вероятность подвергнуться зарастанию</span>. В двух других группах были выявлены 8 значимых факторов с положительной связью, расположившихся в следующем (по значимости) порядке:</p> 
<span class="monospace">
близость к лесному массиву, <br>
близость к городу,<br> 
близость к поселку городского типа,<br> 
близость к дороге,<br> 
близость к реке,<br> 
близость к деревне,<br> 
размер поля.<br>
</span>
<p> Из выявленных восьми факторов интерпретируемы практически все, седьмой (по значимости) – близость к деревне объяснить сложнее: определенную роль могло сыграть обстоятельство передачи ближайших полей в паевую собственность бывших членов совхозов и использование их под частные  огороды, застройку и т.д.., В последующие годы при обезлюдении многих мелких населенных пунктов окраинные поля  подверглись забросу и зарастанию.</p>
<br>


<h3 id="h3_5">20.5. Оценка изменения функций и типологии сельских поселений в связи с трансформацией землепользования</h3>

<p>Полученные данные о характере трансформации свидетельствуют о широком развитии процессов экореабилитации, даже на землях Смоленского района, прилегающего к областному центру, что в свою очередь вызывает необходимость постановки вопросов  о смене набора и приоритета функций, выполняемых сельскими поселениями (производственные, агрорекреационные, рекреационные, экосервисные, иные) и, следовательно – о  перспективах дальнейшего развития сельских районов области, а также о сохранении средостабилизирующей способностей региональной сети охраняемых природных территорий. </p>

<p>Сельское поселение - подходящая операционно-территориальная единица оценки трансформации землепользования по ряду причин. Во-первых, сельские поселения (СП) соразмерны по средней площади единицам ландшафтной структуры территории, приблизительно соответствующим уровню "местности" или "ландшафта" (в системе индивидуального ландшафтного картографирования, и, следовательно - характеризуются совокупностью более-менее однородных природных параметров (рельеф, почвенный покров, характер растительности и т.д.). Во-вторых, СП в значительной степени соответствуют ареалам, на которые проецировалась в советскую эпоху производственно-экономическая активность советских и коллективных хозяйств (колхозов и совхозов); в этом смысле СП являются преемниками некой общей истории освоения, продолжавшейся несколько десятилетий. На момент разработке моделей в пределах Смоленской области насчитывалось <span class="greencursiv">234</span>  сельских поселения.</p>

<p>Оценка преобладающих функций сельских поселений может быть основана на моделировании трех групп факторов, характеризующих ландшафтную структуру, рисунок освоения и актуальный уровень трансформации землепользования.</p> 


<p>В рамках данного блока модели ландшафтная структура определена для <span class="greencursiv">234</span> сельских поселений области: сетка сельских поселений Смоленской области была наложена на  векторную карту ландшафтных местоположений с использованием инструмента <span class="blackbold">Tabulate Intersection</span>  <span class="red">ArcMAP10.x</span> (<span class="monospace">Analysis Tools >> Statistics</span>), который позволяет вычислить площадь каждого типа в пределах конкретного сельского поселения.  Поскольку общее число местоположений зависит не только от ландшафтной структуры местности, но и от размеров сельского поселения (потенциально, чем больше площадь поселений, тем большее число местоположений оно способно вместить), для нормализации результата использован показатель удельного разнообразия - общее число типов местоположений, отнесенное к площади СП.</p>

<p>Для выявления классов сельских поселений, близких  по ландшафтной структуре  использовался кластерный анализ (вариант без пространственных ограничений). Обработка проводилась для различного числа кластеров с экспертной оценкой результатов. Оптимальным был признан результат для <span class="greencursiv">7</span> кластеров. Для получения содержательной интерпретации результатов классификации представлялось важным получить некие усредненные спектры ландшафтных местоположений каждой группы. Для этой цели шейп-файл кластеров обрабатывался инструментом  <span class="blackbold">Summury statistic|Суммарная статистика</span> с номерами кластеров в качестве значений поля  <span class="cursive">case</span> и подсчетом суммарных значений по колонкам каждого вида местоположений. </p>

<a id="01_Кластер_Ланд_Белый"><img src="Pict_5_20/01_Кластер_Ланд_Белый.png" width="65%"></a><br>
<br><span class="imgtitle">Рис. 20.8 Классы СП Смоленской области, выделившиеся по особенностям ландшафтной структуры


<p>Обработка и построение графиков-диаграмм позволили выявить два вида различий. Прежде всего классы сельских поселений различаются по числу относящихся к ним поселений (от <span class="greencursiv">1-2</span> до <span class="greencursiv">16-54</span>). Появление немногочисленных классов объясняется попаданием небольших (как правило, небольших по площади) поселений в уникальные ландшафтные условия (например, поймы и низкие террасы крупных долин) <a href="#Таблица 20.4">Таблица 20.4</a>.</p>


<a id="Таблица 20.4"><span class="imgtitle">Таблица 20.4 Преобладающие типы ландшафтных комплексов в семи классах сельских поселений
<div class="table">
<table id="customers">
        <table border="1" width=100%>

<tr>
<th>ID кластеров</th>
<th>Преобладающие комплекс местоположений</th>
<th>Число СП</th>
</tr>

<tr>
<td>0</td>
<td>Расчлененные верховьями долин склоны высоких морен</td>
<td>50</td>
</tr>
<tr>
<td>1</td>
<td>Привершинные слабодренированные поверхности конечных морен</td>
<td>41</td>
</tr>
<tr>
<td>2</td>
<td>Привершинные поверхности и склоны вторичных моренных равнин, эрозионно-расчлененные</td>
<td>36</td>
</tr>
<tr>
<td>3</td>
<td>Крупные долины в окружении долинных зандров и  склонов моренных равнин</td>
<td>18</td>
</tr>
<tr>
<td>4</td>
<td>Тектонические низменности с озерно-ледниковыми равнинами и низкими террасами</td>
<td>9</td>
</tr>
<tr>
<td>5</td>
<td>Всхолмленные озерно-ледниковые и озерно-водноледниковые равнины</td>
<td>32</td>
</tr>
<tr>
<td>6</td>
<td>Широкие речные долины в окружении комплекса террас и долинных зандров</td>
<td>48</td>
</tr>

</table>

<p>Довольно значительными оказались различия между ландшафтной структурой групп  <a href="#Гистгр_Ланд_Структура">(Рис. 20.9)</a>.</p>

<a id="Гистгр_Ланд_Структура"><img src="Pict_5_20/Гистгр_Ланд_Структура.png" width="70%"></a>
<br><span class="imgtitle">Рис. 20.9 Ландшафтная структура семи классов сельских поселений 

<p>Классификация сельских поселений Смоленской области по характеру освоенности проведена с учетом группы факторов, которые можно считать, с одной стороны – условиями, с другой – результатами освоения. В последнем случае можно констатировать что эти факторы в совокупности формируют так называемый рисунок освоения: соотношение освоенных и природных ареалов. К факторам-условиям можно отнести, заболоченность, заозеренность (для Смоленской области как территории свежего оледенения характерна высокая степень присутствия в ландшафте послеледниковых озер), плотность речной сети. К факторам-результатам освоения – прежде всего площадь СУ, площадь селитьбы, плотность дорожной сети.  Залесенность – уникальный фактор, который является одновременно и фактором-условием, и фактором результатом. Вся совокупность параметров оценивались картометрически с помощью инструмента <span class="blackbold">Table Intersection,</span> который позволяет выявить абсолютные (площадь) и относительные (процент покрытия) показатели; таким образом был выявлен индивидуальный спектра  освоения для каждого СП  - соотношение освоенных и природных ареалов, плотность речной и дорожной сети, плотность населенных пунктов.</p>

<a id="Таблица 20.4"><span class="imgtitle">Таблица 20.4 Преобладающие типы ландшафтных комплексов в семи классах сельских поселений
<div class="table">
<table id="customers">
        <table border="1" width=100%>

<tr>
<th rowspan="2">Номер кластера</th>
<th colspan="7">Доля территории кластера, занятая различными типами использования земель, (%)</th>
</tr>
<tr>
<th>Поля </th>
<th>Леса</th>
<th>Акватории</th>
<th>Болота</th>
<th>Города</th>
<th>Деревни</th>
<th>Дороги</th>
</tr>
<tr>
<td>0</td>
<td>36,33</td>
<td>45,37</td>
<td>2,24</td>
<td>6,99</td>
<td>0,00</td>
<td>1,89</td>
<td>0,56</td>
</tr>
<tr>
<td>1</td>
<td>47,31</td>
<td>28,29</td>
<td>0,27</td>
<td>2,97</td>
<td>0,15</td>
<td>2,40</td>
<td>0,47</td>
</tr>
<tr>
<td>2</td>
<td>39,25</td>
<td>46,30</td>
<td>0,11</td>
<td>3,71</td>
<td>0,24</td>
<td>1,88</td>
<td>0,49</td>
</tr>
<tr>
<td>3</td>
<td>36,48</td>
<td>45,01</td>
<td>2,20</td>
<td>8,48</td>
<td>0,10</td>
<td>1,81</td>
<td>0,58</td>
</tr>
<tr>
<td>4</td>
<td>40,61</td>
<td>41,63</td>
<td>2,03</td>
<td>10,63</td>
<td>0,10</td>
<td>1,67</td>
<td>0,83</td>
</tr>
<tr>
<td>6</td>
<td>33,41</td>
<td>48,79</td>
<td>0,53</td>
<td>6,24</td>
<td>0,26</td>
<td>1,55</td>
<td>0,57</td>
</tr>
<tr>
<td>6</td>
<td>38,15</td>
<td>48,06</td>
<td>0,40</td>
<td>5,17</td>
<td>0,38</td>
<td>1,58</td>
<td>0,59</td>
</tr>
<tr>
<td>7</td>
<td>31,11</td>
<td>54,01</td>
<td>0,55</td>
<td>8,87</td>
<td>0,08</td>
<td>1,31</td>
<td>0,80</td>
</tr>
<tr>
<td>8</td>
<td>34,92</td>
<td>45,76</td>
<td>1,74</td>
<td>9,44</td>
<td>0,32</td>
<td>1,72</td>
<td>0,74</td>
</tr>


</table>

<p>Возможность классификации спектров и оптимальное число групп  определялось с использованием инструмента <span class="blackbold">Анализ группирования</span>    <span class="red">ArcMAP10.x</span>. Затем посредством кластерного анализа были выявлены 9 кластеров сельских поселений, различающихся по пространственному характеру освоения, причем выявленные спектры обнаруживают достаточное четкое расхождение типов с ID от <span class="greencursiv">0</span> до <span class="greencursiv">8</span> <a href="#01_Спектры_Освоения">(Рис. 20.10)</a>.</p>

<a id="01_Спектры_Освоения"><img src="Pict_5_20/01_Спектры_Освоения.png" width="70%"></a>
<br><span class="imgtitle">Рис. 20.10 Спектры пространственного характера освоения сельских поселений Смоленской области

<p><span class="bolditalic">Тип 1</span> и <span class="bolditalic">тип 2</span> – это СП с наибольшей площадью распашки и значительной долей деревень, высокой плотностью дорог, сравнительно низкой заболоченностью и небольшой площадью акваторий значительных рек и водохранилищ. К ним сравнительно близок <span class="bolditalic">тип 0</span>, которому также свойственна высокая доля сельхозугодий, но при этом он отличается более высокой плотностью дорог, присутствием большей поверхности  болот. <span class="bolditalic">Тип 7</span> при средней распаханности является наиболее залесенным и заболоченным, и обладает весьма низкой плотностью дорог и низкой плотностью деревень. К нему приближается <span class="bolditalic">тип 5</span>  также с высокой залесенностью, малой плотностью дорог, но с чуть большей площадью селитьбы. <span class="bolditalic">Тип 6</span> и <span class="bolditalic">тип 8</span> формируют близкую группу со средней распаханностью, высокой лесистостью, средней обводненностью, но отличаются друг от друга плотностью дорог и площадью селитьбы. <span class="bolditalic">Тип 3</span> и <span class="bolditalic">тип 4</span>  обладают средними значениями всех факторов, но при этом – низкой площадью селитьбы. Логичным было бы ожидать проявление связи между долей сельхозугодий и относительною площадью селитьбы, но это условие выполняется не для всех типов СП, таким образом можно констатировать что <span class="bolditalic">традиционная деревня, как основной элемент системы сельского расселения перестает быть организующим звеном системы сельскохозяйственного производства</span>.</p>


<p>Наблюдающиеся тренды в изменении землепользования наряду с унаследованными свойствами и особенностями экономико-географического положения существенно дифференцируют ткань административно-территориального деления уже на нижней ступени ее иерархии – сети сельской поселений. Несмотря на то, что кластеризация специально проведена методом игнорирующим соседство (<span class="blackbold">Hill Climbing</span> Ruben 1967), полученные результаты позволяют констатировать дифференциацию территории на меридионально вытянутые полосы - от западной (внешней) границы области к восточной, граничащей с Московской областью.</p> 

<a id="02_Типолог_СП_Белый_Фон"><img src="Pict_5_20/02_Типолог_СП_Белый_Фон.png" width="60%"></a>
<br><span class="imgtitle">Рис. 20.11 Типы сельских поселений Смоленской области по характеру освоенности территории 

<p>Характерно, что к границе с республикой Беларусь примыкают полосы двух типов СП с низкой площадью селитьбы, что фиксирует депопуляцию сельской местности в этих районах. К ним примыкают типы 3,  4 и 6 тип также характеризующиеся невысокой долей сельской селитьбы, в совокупности формирующие своего рода осевую зону области. Далее на восток расположены полосы типа 5 и 7, с максимальной лесистостью и заболоченностью, но при этом сохраняющие средние показатели распашки и селитьбы. Краевым к Москве является полоса СП 8-го типа с относительно высокой долей водно-болотных угодий и низкой плотностью дорог на фоне сохранившейся селитьбы. Данные результаты можно считать своего рода <span class="bolditalic">status quo</span> - достигнутым положением для середины второго десятилетия XXI века.</p>


<p>Мерой трансформации сельскохозяйственных функций в указанный период стала интенсивность процессов зарастания полей в границах сельских поселений </p>


<a id="Степень_Экореабилитации_СП"><img src="Pict_5_20/Степень_Экореабилитации_СП.png" width="60%"></a>
<br><span class="imgtitle">Рис. 20.12 Степень экореабилитации и зарастания полей в сельских поселениях Смоленской области: 1 – отсутствует; 2 - в начальной стадии; 3 – незначительная; 4 –средняя; 5 – значительная; 6 – высокая; 7 – критическая

<p>Сопоставление полученных данных с выделенными ранее типами поселений по характеру освоенности демонстрирует четко выраженные пространственные закономерности. Поселения, где сельскохозяйственные угодья практически не испытывают процессов зарастания, расположены вдоль двух субширотных полос. Одна следует через центральную часть области вдоль трассы М1, другая – в ее южной части - вдоль трассы А130. Ареалы же сельских поселений с ярко-выраженными процессами зарастания (категории 5-7) образуют такие же полосы в 30 км к северу и к югу от трассы М1. Подобные закономерности отчасти соотносятся с типами поселений по характеру использования земель. Так, процессы зарастания в большей степени характерны для поселений, отнесенных к типам 5,8,9. В то же время утрата сельскохозяйственных функций отмечается и в сельских поселениях, отнесенных к типам с высокой освоенностью – 1 и 2.</p>

<p>Таким образом, несмотря на ограниченность набора привлеченных данных ГИС-моделирование процесса трансформации землепользования может обеспечивать вполне содержательные результаты и способствовать выявлению реальных и трендов.</p>

<p>На современном уровне развития гис-технологий глобальные данные по изменению ландшафтного покрова (Landcover) и трансформации лесного покрова  корректно отражают основные тренды в изменении землепользования и как правило точно локализованы, но требуют региональной верификации относительно региональных проявлений  четырех групп факторов, «ответственных» за их реализацию демутационные сукцессии, антропогенные воздействия, экстремальные природные явления, результаты естественной многолетней динамики ландшафтов. Конкретные региональные изменения также как и вызвавшие их причины могут быть установлены в результате осуществления самостоятельных экспертных сравнительных оценок по снимкам Landsat и расчетов индексов NDVI и параллельного проведения полевых исследований на ключевых участках</p>

<p>Оценка вклада основных природных и социально-географических факторов в трансформацию землепользования может быть осуществлена посредством построения специальной геоинформационной модели, позволяющей учесть три группы факторов – ландшафтно-географические, позиционные, определяющие положение сельхозугодий относительно природных объектов, и ситуационные, описывающие близость к населенным пунктам различного типа, границам различного ранга и  дорожной сети.</p>
<p>Следует особо отметить обнаруженную в ходе гис-моделирования зависимость между сохранностью сельхозугодий и принадлежностью к территории одного хозяйствующего субъекта, что говорит о значимости управленческого (т.е. «субъективного») фактора. Увеличение доли земель занятых селитьбой в Смоленской области (также как и в других областях Нечерноземья), вызвано не ростом численности населения  и не реальным расширением зоны застройки, а изменением правового статуса – переводом земель сельскохозяйственного назначения (паевых участков) в земли населенных пунктов с целью агрегирования, увеличение кадастровой стоимости и дальнейшей перепродажи. В целом проявление связи между долей сельхозугодий и относительной площадью селитьбы выполняется не для всех типов сельских поселений, следовательно можно констатировать что традиционная деревня, как основной элемент системы сельского расселения перестает быть организующим звеном системы сельскохозяйственного производства. </p>

<p>Наблюдающиеся тренды в изменении землепользования наряду с унаследованными свойствами и особенностями экономико-географического положения существенно дифференцируют «ткань» административно-территориального деления уже на нижней ступени ее иерархии – сети сельской поселений. Итоговые типы, отражающие актуальную ситуацию освоения сельских поселений, укладываются на территории Смоленской области в меридионально вытянутые полосы, отражающие «напряженность» социально-экономического «силового» поля, полюсами которого являются с одной стороны столичный московский регион, с другой – республика Беларусь.</p>


<br>
<footer id="main-footer">Пространственный анализ в геоэкологии &copy; Е.Ю.Колбовский, 2022 </footer>
</body>
</html>