Известия РАН. Серия географическая, 2022, T. 86, № 5, стр. 731-745

Механизмы формирования ландшафтно-зональных систем Волжского бассейна (эмпирико-статистическое моделирование)

Э. Г. Коломыц *

Пущинский научный центр РАН, Институт фундаментальных проблем биологии
Пущино, Россия

* E-mail: egk2000@mail.ru

Поступила в редакцию 31.03.2022
После доработки 15.05.2022
Принята к публикации 12.07.2022

Полный текст (PDF)

Аннотация

Изложено численное обоснование концептуального положения отечественной географической экологии о генетическом единстве природной зоны и географического ландшафта. На основе известных методических подходов к полисистемному анализу географической среды проведены количественные оценки ландшафтной организации бореального экотона – ансамбля зональных переходов от бореального пояса к поясу суббореальному на Восточно-Европейском субконтиненте. Задача решалась с помощью эмпирико-статистических моделей ландшафтно-зональных связей. В качестве рабочего инструмента анализа использована разработанная автором кибернетическая модель природных комплексов, описывающая иерархическую систему их организации. Осуществлены три методических подхода к полисистемному анализу: 1) функциональный, с выделением структур экологических ниш ландшафтных единиц в их факторном пространстве; 2) типологический, с дифференциацией территории по группам ландшафтов и с анализом их взаимодействия; 3) иерархический, с использованием кластерных (таксономических) моделей ландшафтной классификации. Получены меры тесноты связей между различными геокомпонентами во множестве ландшафтных единиц, а по этим мерам выделены геосинэргические “ядра” параметрической гомогенности зональных геопространств. Установлены конкретные парциальные ландшафтные структуры, образованные наиболее тесно связанными геокомпонентами и определяющие “лицо” ландшафтно-зональных систем. С помощью информационного моделирования проведена дифференцированная оценка ландшафтной организации природных зон (подзон) путем операций с матрицами эквивалентности геокомпонентных блоков. По теоретико-множественным и кластерным моделям геокомплексных связей описано территориальное пространство бореального экотона и его зональных частей. На основе количественных критериев оценен вклад различных факторов в ландшафтную организацию зонального геопространства. Это позволяет определить меру чувствительности ландшафтной структуры данной территории к изменениям зонально-климатических условий, что имеет непосредственное отношение к оценкам устойчивости ландшафтно-зональных систем в меняющемся климате.

Ключевые слова: природная зона, бореальный экотон, региональный ландшафт, иерархическая природно-территориальная организация, геокомпоненты, эмпирико-статистическое моделирование

ВВЕДЕНИЕ

В учении о географическом способе организации систем земной природы центральное место занимает выдвинутая в свое время В.В. Докучаевым (1881, 1948), а затем всесторонне развитая и теоретически обоснованная А.А. Григорьевым (1966, 1970) концепция о цельности и неразрывности географической среды, получившая наиболее полное развитие в теории географической зональности. Это концептуальное положение комплексной физической географии звучало уже в трудах А. Гумбольдта и П.Н. Савицкого. Дальнейшее развитие оно получило также в фундаментальных работах Л.С. Берга, Г.Д. Рихтера, В.Р. Волобуева, Ф.Н. Милькова, М.И. Будыко о зональных ландшафтных комплексах и вещественно-энергетических факторах и процессах, определяющих географическую зональность. Предложенная В.В. Докучаевым концепция природной зональности превратилась в один из теоретических инвариантов современной географии (Преображенский и др., 1997). Л.С. Берг (1958, с. 210) был, по-видимому, первым, кто охарактеризовал географические зоны как ландшафтные, подчеркнув глубокое единство понятий “ландшафт” и “природная зона”.

Следует особо подчеркнуть тесную связь последующего совершенствования теории зональности как общепланетарного биоклиматического явления с развитием ландшафтоведения, а затем и геоэкологии, что отобразило общую тенденцию сближения географии и экологии (Сочава, 1970, 2005). Экологический подход позволил существенно расширить рамки такого, уже традиционного объекта физической географии, как природная зональность (Д.Л. Арманд, 1975; Базилевич и др., 1986). В структуре и функционировании зональных типов ландшафтов и региональных экосистем, прежде всего, в сложности строения фитобиоты и ее продуктивности, видится проявление высшей формы организации биосферы. Тем самым, создается наиболее оптимальная для биоты форма ее устойчивого функционирования (Герасимов, 1985).

Дальнейшее развитие данного направления географической (в том числе ландшафтной) экологии связано с разработкой проблем иерархии природно-территориальных образований, которая обусловлена отношениями и связями между геокомпонентами и геопотоками, с одной стороны, и пространством создающихся ими структур – с другой (Forman, 2006; Phillips, 2002; Pickett and Cadenasso, 2002). Утверждается, что ландшафтное исследование должно быть полимасштабным, поскольку большинство ландшафтообразующих процессов и явлений может быть понято только при анализе смежных иерархических уровней (Marceau, 1999; Turner et al., 2001; Wu and Hobbs, 2002). При этом главное внимание необходимо обратить на свойства ландшафта как пространственно-гетерогенной системы: сложность, нелинейную динамику, эмерджентность, критические состояния, самоорганизацию.

Данное направление получило существенное развитие в отечественном ландшафтоведении (Арманд, 1988; Зейдис и др., 2001; Пузаченко, 2006; и др.). При освещении проблемы иерархии и масштаба установлены “резонансные” сочетания масштаба межкомпонентных отношений, создающих целостные ландшафтные структуры (Хорошев и др., 2010). Разномасштабная система связей “вмещаемых” и “вмещающих” геокомпонентов прослежена на примерах бинарных отношений “рельеф (геологическое строение)–растительность” (Хорошев, 2010). В обобщающей работе (Хорошев, 2016) представлена концепция, развивающая моносистемную ландшафтную парадигму Докучаева, с построением количественных (статистических) моделей межкомпонентных связей в полимасштабном географическом пространстве и с достаточно полной ландшафтно-экологической интерпретацией этих моделей. Раскрыты механизмы ландшафтной организации не только на различных структурных уровнях, но и связи межуровневые, характеризующие системообразующее участие геокомпонентов с различными характерными временами и соответствующей размерностью.

СТРАТЕГИЯ НАУЧНОГО ПОИСКА

Представленные в настоящем сообщении результаты исследования механизмов формирования ландшафтно-зональных геопространств опираются на изложенные научно-методические положения современной географической экологии. Цель научного поиска состояла в установлении характера и степени влияния гидро-климатических факторов, определяющих “лицо” “вмещающих” зональных геопространств, на территориальную упорядоченность “вмещаемых” низших, согласно (Исаченко, 1965), региональных классификационных единиц – ландшафтов, выделяемых в качестве элементов этих пространств. Необходимо было выявить ландшафтообразующую роль факторов зональности, а также территориальную дифференциацию ландшафтных единиц по сочетанию определенных значений этих факторов. Следовало рассмотреть также обратную связь – иерархическую систему ландшафтообразования, которая обеспечивает определенную устойчивости зональных геопространств.

Таким образом, объект исследования представлял собой достаточно сложное географическое образование, поэтому в качестве методической основы исследования использованы предложенные А.К. Черкашиным (2005, 2020) методические подходы к полисистемному анализу. Изучаемая полисистема (множество ландшафтно-зональных единств основного водосбора Волжского бассейна) расчленялась на три разнокачественных подсистемы (подмножества) – функциональную, объектно-типологическую и иерархическую, по каждой из которых создавались эмпирико-статистические модели ландшафтно-зональных связей. Этот подход созвучен трем известным в физической географии методическим положениям: межкомпонентной корреляции, сопряженности между геокомплексами и их соподчиненности в иерархическом ряду (Сочава, 2005).

Организация зональных геопространств тесно связана с пространством экологическим – системой координат геокомпонетных признаков, в которой формируется множество экологических ниш фитоценологических элементов природных комплексов. В данном случае использованы два из трех известных (Одум, 1975) аспектов экологической ниши: 1) пространственный, близкий к понятиям местообитания и местоположения, и 2) ресурсный (он же факторный), определяемый положением ниши объекта относительно градиентов того или иного вещественно-энергетического фактора (температуры, влажности и др.). Соответственно выделяeтся экологическое пространство двух категорий: факторное, или ресурсное (геокомпонентное) и территориальное (геокомплексное). Размещение экологических ниш биогеосистем в зональном геопространстве соответствует климатическим нишам доминирующих видов-эдификаторов и образуемых ими растительных сообществ. Различные сочетания экологических ниш ландшафтов образуют множество экологических пространств вышестоящих по рангу геосистем зонально-регионального уровня.

Исследование проведено на примере ансамбля природных зон Средней полосы Русской равнины, основную часть которой составляет Волжский бассейн. Эта территория в силу своих природных особенностей традиционно была объектом изучения зональных форм организации географической среды. На ее примере были заложены основы концепции о зональных ландшафтных комплексах и о вещественно-энергетических факторах и процессах, определяющих географическую зональность.

ИСХОДНЫЕ МАТЕРИАЛЫ И МЕТОДЫ АНАЛИЗА

Средняя полоса Русской равнины, определяемая как бореальный экотон (Коломыц, 2005), является обширной полосой перехода от бореального пояса к суббореальному. В пределах основного водосбора Волжского бассейна она охватывает среднюю и южную тайгу, подтайгу, лесостепь и северную степь. “Сердцевину” бореального экотона образует Главный ландшафтный рубеж Русской равнины (Мильков, 1981), проходящий по речным долинам Оки–Волги–Камы и отделяющий подтаежную зону от зоны лесостепи, с резко выраженным переходом соотношения тепла и влаги через 1 (Базилевич и др., 1986).

Исходная информация для анализа снималась с построенной в масштабе 1 : 2 500 000 ландшафтной карты (автор – В.П. Юнина) и 25 созданных нами карт ландшафтно-геофизических условий основного водосбора Волжского бассейна и его окружения (Коломыц, 2005). Ландшафтная карта создавалась автором на основе известных мелкомасштабных геоморфологических, почвенных и геоботанических карт, а также собственных выборочных полевых исследований. Следует отметить, что геоботанические карты 35–40-летней давности (Исаченко, Лавренко, 1974; Котова, 1987) по своему содержанию не потеряли актуальности до настоящего времени. Установлено (Исаев, 2012), что за период с 1968 по 2008 г. лесистость Европейской части России изменились всего на 6–8%. Несущественно изменилась и породная структура лесов: доля хвойных пород снизилась на 4%, а доля мягколиственных – увеличилась на 5%. Осталось практически таким же соотношение площадей с различной возрастной структурой лесов.

Новейшая “Карта лесных экосистем Северной Евразии”, составленная по спутниковым данным SPOT-Vegetation (Барталев и др., 2011; Исаев, 2012) не могла быть использована. Выделенные на ней “типы растительности” (например, хвойные вечнозеленые леса, хвойные листопадные, лиственные леса вообще, смешанные с преобладанием хвойных, и т.п.) имеют скорее лесопромышленное, нежели лесоведческое значение, что затрудняет интерпретацию этой карты в свете классических законов лесной биогеоценологии, согласно (Сукачев, 1972). Такие категории лесов совершенно не соответствуют принятому в традиционной геоботанике содержательному смыслу понятий “классы растительных формаций” и “типы растительности” (Сочава, 1979), на основе которых и выделяется система природных зон и подзон Волжского бассейна. Более того, в этой легенде не находят отражения важнейшие принципы фитоценологии и биогеоценологии (Базилевич и др., 1986; Вальтер, 1982; Сочава, 1979; Сукачев, 1972), в том числе понятия коренных и производных ассоциаций, эндоэкогенетические и экзогенные динамические смены растительного покрова. Соответственно отсутствует разделение лиственных лесов на широколиственные (коренные), по которым выделяются, как известно, зональные типы и подтипы суббореальной растительности, и леса мелколиственные (производные). Наконец, нет разделения пространств хвойных вечнозеленых пород на ареалы темно- и светлохвойных, что лишает возможности выделять зональные и подзональные биогеосистемы в самом бореальном поясе.

Для съема информации с ландшафтной карты Волжского бассейна применялся известный метод ортогональных сеток. Шаг между узлами сетки был, как правило, меньше среднего поперечника видового ландшафтного контура. Вся территория была покрыта 1467 точками. Аналогичным образом использовались опубликованные картографические материалы по Русской равнине.

Исходная ландшафтная карта включала 67 видов ландшафтов, которые дифференцируются в значительной мере по азональным (геоморфологическим) признакам, согласно (Исаченко, 1965). В целях более четкой зональной “экологизации” ландшафтной карты было проведено повышение ранга ландшафтных единиц до уровня типологических групп ландшафтов, в трактовке А.В. Николаева (1978), путем введения полыновской системы мезомасштабных ландшафтных сопряжений. Тем самым удалось сблизить ранги рассматриваемых нами ландшафтных единиц с единицами растительного и почвенного покровов. В пределах лесных зон и лесостепи выделено 20 типологических групп ландшафтов (рис. 1).

Рис. 1.

Карта типологических групп ландшафтов территории Волжского бассейна и его окружения: а – территории и акватории речных долин, озер и водохранилищ; б – заболоченные территории; в – границы групп ландшафтов; г – границы природных зон и подзон. 1–20 – порядковые номера типологических групп ландшафтов.

Решение поставленных задач осуществлялось с помощью построения серии взаимосвязанных эмпирико-статистических моделей:

1) теоретико-информационных (А.Д. Арманд, 1975; Пузаченко, Скулкин, 1981), описывающих пространственную геокомпонентную сопряженность, т.е. функциональное расслоение, по терминологии (Черкашин, 2005), с помощью которой был создан региональный банк экологических ниш видов ландшафтов, характеризующий их факторное (параметрическое) пространство;

2) теоретико-множественных (Андреев, 1980), вскрывающих территориально дифференцирующий аспект ландшафтной организации (объектно-типологическое расслоение) путем операций с экологическими нишами (дескриптивными векторами) ландшафтных систем, с получением различных мер их сходства и включения;

3) кластерных (таксономических) моделей (Дюран, Оделл, 1977; Жуков и др., 1980), характеризующих соподчиненную структуру ландшафтно-зональных связей (иерархическое расслоение).

Многомерный ландшафтно-экологический анализ зональных геопространств Русской равнины проведен в рамках понятийной концепции о структурных уровнях природно-территориальной организации, опирающейся на кибернетическую модель природного комплексафон–каркас–узоркак иерархической системы управления (Коломыц, 1988). Дадим некоторые пояснения. Территориальное многообразие природной среды представлено функциями зонального фона и пространственного каркаса, а также расслоенным пространством ландшафтного узора (эпигенетической структуры) по элементам каркаса в конкретных условиях. Физико-географический фон характеризует общий материально-энергетический уровень природно-территориальной организации. Фоновая функция – это непрерывное распределение признака без заметных скачков. Каркас ландшафта представляет собой совокупность наиболее пространственно протяженных и наименее изменчивых во времени структурно-функциональных элементов, которые соответствуют территориальному масштабу данной геосистемы и определяют соответствующую относительно замкнутую сеть переноса вещества и энергии, а также узлы и линии точек поворота геопотоков. Ландшафтный узор – это ландшафтно-типологическая дифференциация (расслоение) зонального фонового пространства на базе территориального и факторального каркаса. Геополя и геопотоки при граничных условиях физико-географического фона и данного каркаса формируют процессор; его характеристики являются факторами прямого ландшафтообразующего действия.

РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ

Факторное пространство и общая таксономическая модель природной зональности

В ландшафтно-зональном геопространстве каждый геокомпонент рассматривался в системе структурных уровней “фон–каркас–узор”, где фон и каркас с процессором (геопотоками) образуют входные переменные, а группы ландшафтов выступают как выход системы (узор). Все геокомпонентные признаки распределены по блокам кибернетической модели (табл. 1).

Таблица 1.  

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

Геокомпонентные признаки Параметры связей
Но-мер Наименование и обозначение
Физико-географический фон и ландшафтный каркас К(А;В) К(В/А)
1
2
3
4
5
6
7
8
9
10
11
12
13
Годовая суммарная радиация, QС
Средняя температура января, t (янв)
Коэффициент зимней континентальности Полозовой, Кз(Пол)
Годовой поверхностный сток, S(год)
Максимальная высота снежного покрова, h (сп)
Снежно-температурный коэффициент Рихтера, СТК
Коэффициент снежности Осокина, К(сн)
Типы морфоструктур, МСТ
Морфоструктура и морфоскульптура, МСТ + МСК
Ступени абсолютных высот, Н(абс)
Современные тектонические движения, F
Механический состав почвообразующх пород, Спп
Роды ландшафтов, РЛ
0.120
0.085
0.168
0.145
0.099
0.104
0.202
0.199
0.140
0.174
0.067
0.125
0.249
0.238
0.224
0.357
0.350
0.251
0.267
0.251
0.273
0.321
0.200
0.176
0.183
0.321
Нерасчлененная система блоков “каркас–процессор”
14
15
16
17
18
19
Годовой радиационный баланс, R(год)
Годовая сумма осадков, r(год)
Годовой подземный сток, U(год)
Коэффициент стока, f
Глубина залегания грунтовых вод, Z(гв)
Литология и влажность грунтов, SW(гр)
0.166
0.122
0.183
0.129
0.274
0.175
0.343
0.298
0.395
0.295
0.275
0.253
Процессор (внутренние геопотоки) 
20
21
22
23
24
25
26
27
28
Средняя температур июля, t (июля)
Сумма биологически активных температур,∑ t(cр-cут) ≥ 100
Годовая испаряемость, Е0
Сумма осадков теплого периода, r(тп)
Летние запасы влаги в почве, W(почв)
Радиационный индекс сухости Будыко, I(Буд)
Коэффициент увлажнения Высоцкого−Иванова, К(увл)
Гидротермический коэффициент Селянинова, ГТК
Первичная биопродуктивность, Р(перв)
0.187
0.206
0.210
0.140
0.296
0.162
0.169
0.184
0.081
0.379
0.404
0.362
0.317
0.303
0.361
0.376
0.379
0.166
Ландшафтный узор
29
30
31
Группа видов почв, ГВП
Химический состав грунтовых вод, J(гв)
Почвенно-геохимические комплексы, SG
0.125
0.183
0.142
0.309
0.212
0.379

Примечание. В общем списке геокомпонентных признаков отсутствуют группы растительных формаций, поскольку по данному признаку выделялись сами виды ландшафтов в пределах каждого рода, а также зональных типов и подзональных подтипов.

Согласно (А.Д. Арманд, 1975; Пузаченко, Скулкин, 1981), рассчитывались две меры пространственных связей явления А (типологической группы ландшафтов) и фактора В (геокомпонентного признака), в долях единицы:

1) нормированный коэффициент сопряженности K(А;В):

(1)
$K\left( {А;В} \right) = \frac{{{{2}^{{T(AB)}}} - 1}}{{{{2}^{{H({\text{min }}A,B)}}} - 1}},$

2) коэффициент приема информации K(B/A):

(2)
$K({B \mathord{\left/ {\vphantom {B A}} \right. \kern-0em} A}) = \frac{{T(AB)}}{{H(B)}}.$
Здесь 2Т(АВ) – число состояний явления А и фактора В, H(A) и H(B) – общие меры их разнообразия (априорная негэнтропия), Т(АВ) – общее количество информации, переданное от фактора к явлению. Значение К(А;В) = 0.19 примерно соответствует коэффициенту корреляции 0.7 (Пузаченко, Скулкин, 1981).

Результаты информационно-картографического анализа пространственной сопряженности групп ландшафтов (ГЛ) с характеристиками различных функциональных блоков природного комплекса (см. табл. 1) выявили ведущую роль не только первичных входных (фоново-каркасных), но и процессорных вещественно-энергетических параметров, определяющих общую зональную структуру Русской равнины. В качестве доминирующих начальных входных переменных, обладающих максимальной независимостью, оказались: годовая суммарная радиация (Qc), годовая сумма осадков (r(год)), типы морфоструктур и морфоскульптур (МСТ + МСК) и механический состав почвообразующих пород (Сnn). Следуя методу, изложенному в (Пузаченко, Скулкин, 1981), зависимость распределения площадей групп ландшафтов по всему зональному спектру Волжского бассейна от указанных входных переменных была выражена в виде следующего линейного полинома, имеющего чисто иллюстративный характер:

(3)
$\begin{gathered} {\text{ГЛ}} = 0.24Qc + 0.30r(год) + {\text{ }}0.32(МСТ + MСК) + \\ + \,\,0.12Cnn + {\text{ }}0.02Х. \\ \end{gathered} $
Здесь коэффициентами при аргументах служат коэффициенты приема информации (К(В/А)) (см. табл. 1). Сумма всех коэффициентов, включая коэффициент при неизвестном аргументе Х, равна 1.

Как видно из уравнения, дифференциация групп ландшафтов бореального экотона почти целиком определяется влиянием четырех указанных факторов. При этом, роль климатических и литогенных входных переменных соразмерна, с некоторым “перевесом” (до 54%) факторов климатической группы. В составе последней влияние поступающей солнечной энергии и атмосферной влаги так же примерно одинаково, между тем как в литогенной группе факторов ведущую роль играют генетические типы рельефа (МСТ + МСК). Значение механического состава почвообразующих пород на региональном уровне оказалось гораздо менее значимым.

Ландшафтная структура зональных геопространств выражена таксономическими моделями групп ландшафтов по отношениям сходства их экологических ниш в многомерном пространстве ландшафтно-геофизических факторов. Такая модель построена по всему набору 21 обменно-транзитных признаков (1–7, 14–18 и 20–28, см. табл. 1), с применением метода “вроцлавской таксономии” (Жуков и др., 1980). Для этого сначала строилась матрица распределения всех ландшафтно-геофизических параметров (суммарной радиации, радиационного баланса, средней температуры июля и т.д.) по 21 группе ландшафтов. В каждую ячейку этой матрицы помещалась “норма” параметра, свойственная данной ландшафтной группе. При вычислении нормы использовались центральные значения всех градаций данного ландшафтно-геофизического параметра, а также элементы матрицы К(А;В), которые были введены в расчеты в качестве “весовых” коэффициентов. Затем вычислялись среднеквадратичные отклонения и нормированные значения признаков, по которым составлялась матрица нормированных значений по всем ГЛ. Наконец, с использованием евклидовой метрики и “весовых” коэффициентов рассчитывались парные таксономические расстояния между ГЛ по всем признакам. В качестве “весовых” коэффициентов использованы приведенные к 1 значения параметра К(А;В).

По матрице расстояний построено граф-дерево минимальной длины связи между ГЛ (рис. 2). В основании графа лежит связка среднетаежных геосистем (1–3), охватывающая почти все ландшафтное разнообразие этой подзоны: эрозионно-денудационные, моренные, денудационно- и моренно-зандровые пихто-ельники и сосняки на элювиальных, транзитных и трансэлювиальных местоположениях. От среднетаежной плеяды отходят две ветви последовательного сходства, каждая из которых соответствует генетическим типам рельефа – родам ландшафтов. Левая ветвь состоит из трансзональной цепи ландшафтных групп: 1–5–9–13, 14–18–19, – которые имеют большие таксономические расстояния. Эта цепь характеризует водораздельные природные комплексы и пересекает весь бореальный экотон от средней тайги (ГЛ 1), через южную тайгу (5) и подтайгу (9), далее к широколиственным лесам (13, 14) и, наконец, к лесостепи (18) и северной степи (ГЛ 21).

Рис. 2.

Дерево минимальной длины между группами ландшафтов Волжского бассейна по 25 ландшафтно-геофизическим признакам (см. табл. 1). Группы ландшафтов: Средняя тайга: 1 – эрозионно-денудационные и денудационно-зандровые, Э, пихтово-еловые и елово-сосновые; 2 – моренные и моренно-зандровые, Э и Т, пихтово-елово-сосновые и сосновые; 3 – зандровые, ТА, елово-сосновые и сосновые; 4 – озерно-ледниковые и долинно-зандровые, еловые и сосновые. Южная тайга: 5 – эрозионно-денудационные, Э, пихтово-еловые; 6 – моренные и озерно-ледниковые, Э и ТЭ, еловые и елово-сосновые; 7 – денудационно- и моренно-зандровые; Т, пихтово-еловые; 8 – зандровые и долинно-зандровые, ТА и А, елово-сосновые. Смешанные леса: 9 – эрозионно-денудационные, Э и ТЭ, широколиственно-пихтово-еловые и дубово-липовые; 10 – моренные и озерно-ледниковые, Э и ТА, еловые и елово-сосновые; 11 ‒ денудационно- и моренно-зандровые, Т, пихтово-еловые и елово-сосновые; 12 – зандровые и долинно-зандровые, ТА и А, сосновые и елово-сосновые. Широколиственные леса: 13 –эрозионно-денудационные, Э, ясенево-дубово-липовые и сосново-дубовые; 14 – эрозионно-денудационные, Э и ТЭ, дубово-липовые и ясенево-липово-дубовые; 15а – эрозионные и озерно-ледниковые, Т, дубово-липовые и сосново-широколиственные; 15б – денудационно- и моренно-зандровые, Т, хвойно-широколиственные и сосновые; 16 – зандровые и долинно-зандровые, ТА и А, широколиственно-хвойные и сосновые. Средняя и южная лесостепь: 17 – эрозионно-денудационные, Э и ТЭ, ясенево-липово-дубовые, с луговыми степями; 18 – эрозионные и “сыртовые”, ТЭ и Э, луговостепные с липово-дубовыми лесами; 19–20 – озерно-ледниковые и древнеаллювиальные, Т и А, лугово-степные, широколиственно-сосновые. Северная степь: 21 – эрозионно-денудационные, “сыртовые” и древнеаллювиальные, ТЭ и ТА, разнотравно-типчаково-ковыльные. Типы мезоместоположений: Э – элювиальный, ТЭ – трансэлювиальный, Т – транзитный, ТА – трансаккумулятивный, А – аккумулятивный.

Правая ветвь графа соответствует средним и нижним звеньям мезомасштабных ландшафтных сопряжений. Это речные долины и понижения с прилегающими к ним нижними и средними склонами междуречий. Данная ветвь прослеживается только до зоны смешанных лесов включительно и отсутствует на южном крыле бореального экотона. В отличие от левой (водораздельной) ветви, она имеет меньшие таксономические расстояния между своими звеньями. Это указывает на более высокую однородность биоклиматических условий низменных равнин и речных долин по сравнению с водоразделами.

Таким образом, вырисовываются два различных по своей ландшафтно-геофизической структуре трансзональных экологических пространства. К истокам такой макродифференциации восходят основные причинные механизмы территориальной организации региональных геосистем. Одно из ярких последействий этих форм причинности – характерное для Русской равнины явление прямой и обратной вертикальной дифференциации ландшафтов (Мильков, 1981). Первое, водораздельное, экопространство полностью трансзонально, однако более “размытое”; оно отличается меньшей однородностью ландшафтно-геофизических условий между звеньями цепи зональных переходов. Второе, низменно-долинное, пространство ограничено только северным крылом экотона и не простирается южнее зоны смешанных лесов, поэтому его можно именовать собственно бореальным экопространством низменных равнин.

Функциональная организация ландшафтно-зональных геопространств

Качественная определенность каждой группы ландшафтов и особенности размещения ее ареалов должны определяться характером контактной геокомпонентной структуры и степенью ее геосинэргической однородности (Хаазе, 1980). Степень “контактности и геосинэргичности” ландшафтной структуры природных зон (подзон) была выражена через меру тесноты связей (МТС) между признаками различных геокомпонентов. Такой подход основан на представлениях о преобладании в географической оболочке не функциональных, а стохастических связей, с принципиально вероятностным характером географического детерминизма (Арманд и др., 1969). С этой целью все рассматриваемые геокомпонентные признаки были сгруппированы в 7 блоков (табл. 2), которые характеризуют три входных системообразующих фактора (приземную атмосферу, континентальную гидросферу и поверхностную литосферу), а также две выходных характеристики ландшафтов – первичную биопродуктивность и почвенно-геохимические комплексы.

Таблица 2.  

Распределение структурно-функциональных признаков геосистем по основным геокомпонентным блокам

Блок Обобщенное наименование Номерa признаков (см. табл. 1)
A
B
C
D
E
F
G
Радиационно-термический режим атмосферы
Составляющие водного баланса территории
Комплексные гидротермические показатели
Геолого-геоморфологические условия
Гидроэдафические условия местоположений
Первичная биологическая продуктивность
Почвенно-геохимические комплексы
1, 2, 14, 20–22
4, 5, 15–17, 23
3, 6, 7, 25–27
8–13
18, 19, 24, 30
28
31

На основе матриц частот совместной встречаемости каждой пары А и В признаков (по их соответственно i-м и j-м градациям) рассчитывались парциальные коэффициенты связи, согласно (А.Д. Арманд, 1975):

(4)
${{K}_{{ij}}} = 2{{P}_{{ij}}}/\left( {{{P}_{i}} + {{P}_{j}}} \right),$
где Pij – условные, а Pi и Pj – априорные вероятности. Очевидно, что 0 < Kij < 1. Для каждой точки измерительной сети составлялась матрица Kij между всеми 30 признаками. В каждой клетке матрицы помещались те значения Kij, которые отвечают парному сочетанию градаций признаков, принадлежащему данной точке. Объединением точек по ареалам групп ландшафтов получали значения МТС для соответствующих объектов.

Далее, из начальных (полных) матриц мер тесноты межкомпонентных связей удалялись все значения Kij между признаками внутри каждого геокомпонентного блока и оставлены только межблочные отношения, которые образовали вторую группу матриц. Здесь параметр Kij характеризовал уже меру связей между пятью входными абиотическими геокомпонентными блоками (А, В, С, D и Е), а также между двумя выходными характеристиками (P и G, см. табл. 2). Тем самым был получен показатель степени взаимной эквивалентности (конформности) коррелятных геокомпонентных признаков в пределах ареалов данной ГЛ, а также в любом выбранном для анализа вышестоящем по рангу геопространстве. По данному показателю были установлены семь классов эквивалентности структурно-функциональных блоков в ландшафтах Волжского бассейна (табл. 3).

Таблица 3.  

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

Показатель  Класс эквивалентности
I II III IV V VI VII
Коэффициент связи 0.504–
–0.307
0.306–
–0.255
0.254–
–0.216
0.215–
–0.185
0.184–
–0.146
0.145–
–0.110
0.109–
–0.037

На мерах тесноты связей между ландшафтообразующими признаками построены кластерные модели геокомпонентных структур всего бореальнoго экотона и его зональных подразделений. На общей кластер-дендрограмме (рис. 3) выделяются пять основных кластеров, объединяющие три первых уровня эквивалентности признаков. Наиболее высокий уровень образуют следующие связки признаков, которые назовем геосинэргическимиядрамигеокомпонентной структуры природных зон (обозначения см. в табл. 1): в кластере a – [t(янв)–S(год)]–Е0–t(июл); в кластере b – [h(сп)–Кл(Пол)][∑t(cр-cут) ≥100–r(хп)–С(nn)– r(год)–J(гв); в кластере c–[К(сн) – ГТК]; в кластере d – [( f – СТК)][( I(Буд) – F)]–SW(гр)[QС(МСТ + + МСК)]–К(увл)–R(год)–Р(перв); в кластере e [SG–r(тп)]–Z(гв).

Рис. 3.

Кластер-дендрограмма тесноты связей между ландшафтообразующими признаками (см. табл. 1) в целом по всему бореальному экотону. 0–0.1–0.2–0.3 – евклидово расстояние. I–V – уровни эквивалентности признаков. ае – основные кластеры.

Характерно множество связей между гидро-климатическими факторами, которое существенно повышаeт геосинэргическую однородность зональных геопространств, а также вносит приоритетный вклад в разнообразие ландшафтов. Признаки же литогенной группы не имеют аналогичного количества перекрестных взаимосвязей и поэтому не могут создавать подобного разнообразия. По кластер-дендрограммам выявлены следующие свойства параметрического пространства бореального экотона Русской равнины.

1. Признаки ландшафтного каркаса и фона обладают большей геосинэргической гомогенностью, нежели процессорные параметры. Локально-региональные геопотоки преломляют начальные связи между входными переменными в сторону их определенной разбалансированности, создавая мозаику местных отклонений от фоново-каркасной “нормы” и повышая уровень разнообразия ландшафтов.

2. По геосинэргетике ландшафтного каркаса средне- и южнотаежные геопространства наиболее близки к обобщенной ситуации на экотоне. По мере движения в сторону широколиственных лесов и южной лесостепи наблюдается усиление тесноты взаимосвязей каркасных признаков. Это указывает на общее повышение единства первичных входных переменных на южном крыле экотона, где значительное их число образуют кластеры уже на I уровне эквивалентности. На северном же крыле экотона, особенно в зоне смешанных лесов, объединение каркасных признаков в кластеры начинается только со II-го уровня.

3. Среди процессорных признаков аномально высокой геосинэргичной однородностью отличается средняя тайга. В этом “ядре типичности” бореального пояса взаимообусловленность геокомпонентов явно усиливается локально-региональными геопотоками уже на II–III уровнях эквивалентности, благодаря чему обеспечивается относительно высокое генетическое единство собственно бореальных ландшафтов. В более южных зональных геосистемах число кластеров в блоке процессора остается весьма большим при всех уровнях общности. Это означает, что местные геопотоки здесь дробят параметрическое пространство на группы слабо связных ландшафтообразующих факторов. Таким образом, в пределах бореального экотона формируется комплекс процессорных признаков (прямого ландшафтообразующего воздействия) с гораздо более низкой взаимной эквивалентностью, нежели в “ядре типичности” бореального пояса. В этом видится одна из причин повышенной чувствительности зональных подразделений экотона и границ между ними к фоновым возмущениям гидро-климатической системы.

Ландшафтная типология зональных геопространств

Типологическоe расслоение ландшафтно-зональной полисистемы Волжского бассейна проведено по шести блокам структурно-функциональных признаков ландшафтов (см. табл. 2). Выполнялось теоретико-множественное моделирование сначала групп ландшафтов как однородных единиц, а затем и зональных геопространств, с их гетерогенной ландшафтной структурой.

Оценивались бинарные отношения (меры сходства и включения) ГЛ по множеству факторов их формирования, объединенных в указанные группы (см. табл. 2). Каждая пара групп ландшафтов А и В характеризовалась одним и тем же n числом элементов ai и bi из данной группы факторов, с определенными численными значениями. Меры сходства N(ai, bi), а также взаимного включения J(ai; bi) и J(bi; ai) элемента ai в bi и элемента bi в ai находились из соотношений (Андреев, 1980):

(5)
$\begin{gathered} N\left( {{{a}_{i}}{{b}_{i}}} \right) = \frac{{2({{a}_{i}} \cap {{b}_{i}})}}{{{{a}_{i}} + {{b}_{i}}}};\,\,\,\,J\left( {{{a}_{i}};{{b}_{i}}} \right) = \frac{{{{a}_{i}} \cap {{b}_{i}}}}{{{{a}_{i}}}}; \\ J\left( {{{b}_{i}};{{a}_{i}}} \right) = \frac{{{{a}_{i}} \cap {{b}_{i}}}}{{{{b}_{i}}}}. \\ \end{gathered} $

По мерам сходства проведена оценка парциального влияния различных гидро-климатических факторов на ландшафтную структуру зональных геопространств. На графах сходства (рис. 4) были выделены корреляционные плеяды ландшафтных объединений путем задания двух порогов однородности: при N(A, B) > 0.40 и при N(A, B) > 0.60.

Рис. 4.

Корреляционные плеяды групп видов ландшафтов Волго-Вятского региона в их зависимости от годового радиационного баланса (а), суммы биологически активных температур (б) и годового количества осадков (в). 1, 2, 3, … – группы ландшафтов (см. рис. 1). Мера сходства: 1 – 0.01–0.020; 2 – 0.21–0.40; 3 – 0.41–0.60; 4 – 0.61–0.80; 5 – 0.080–0.100. 6 – границы между плеядами первого уровня сходства. I–V – корреляционные плеяды.

По величине N(A, B) отчетливо выражена климатическая граница между бореальным и суббореальным поясами (см. рис. 4а). При этом лесостепь (вместе с широколиственными лесами) менее однородна, чем средняя и южная тайга и подтайга. По сумме биологически активных температур (см. рис. 4б), годовой испаряемости и радиационному индексу сухости множество ГЛ разбивается на пять зональных типов (подтипов) уже на первом уровне сходства, с четким обособлением средней и южной лесостепи. Наиболее резок переход от южной тайги к смешанным лесам. Эту границу, по-видимому, и следует считать основным климатическим рубежом между бореальным и суббореальным поясами. Годовое количества осадков предопределяет весьма “размытую” дифференциацию зональных геопространств (см. рис. 4в).

Меры включения позволяют упорядочить рассматриваемые геосистемы в ряды экологических групп, или экотипов, по (Одум, 1975), которые характеризуют уровень структурно-функционального развития каждой из них, а также эффективность использования биотой ресурсов абиотической среды. В общем случае J(ai;bi) ≠ J(bi;ai), т.е. бинарные отношения включения нетранзитивны. При J(A;Bi) > J(B;A) можно говорить о том, что геосистема А проще по структуре и функционированию и менее устойчива, чем система В, в отношении данного фактора, а при J(A;Bi) < J(B;A) – напротив, она оказывается более высокоразвитой и более устойчивой.

В качестве примера приведем теоретико-множественную модель ландшафтно-типологической организации природных зон (подзон) Волжского бассейна по фактору среднеиюльской температуры (рис. 5). На орграфе отчетливо виден, прежде всего, резкий контраст в отношениях включения ГЛ между средней и южной тайгой, с одной стороны, и смешанными лесами – с другой. Таежные геосистемы отличаются большими пересечениями экологических ниш, т.е., более высокой повторяемостью в пространстве одних и тех же термических условий летнего сезона. Судя по количеству и “мощности” входных стрелок (ребер орграфа), наибольшим разнообразием по фактору среднеиюльской температуры отличаются: в подзоне средней тайги трансэлювиальные (2) и транзитные (3) ландшафтные группы (моренные, и зандровые пихто-ельники и сосняки), а в южнотаежной подзоне – ландшафты 5-й элювиальной группы. Наименее разнообразны среднетаежные эрозионно-денудационные и денудационно-зандровые пихто-ельники и сосняки (ГЛ 1), с максимальными числом и “мощностью” выходных стрелок).

Рис. 5.

Орграф отношений включения групп видов ландшафтов таежной и подтаежной зон по средней температуре июля. Мера включения: а – ≤ 0.30; б – 0.31–0.40; в – 0.41–0.50; г – 0.51–0.60; д – 0.61–0.70; е – 0.71–0.80; ж – 0.81–0.92. 1, 2, 3, … – группы ландшафтов (см. рис. 1).

ЗАКЛЮЧЕНИЕ

1. Проведено численное обоснование концептуального положения отечественного естествознания о генетическом единстве природной зоны и географического ландшафта. С помощью полисистемного моделирования описаны причинные механизмы формирования ландшафтно-зональных геопространств при взаимодействии разнокачественных природных компонентов – гидро-климатических, литогенных и почвенно-биотических. Впервые дана интерпретация послойной функции факторного воздействия в разных математических формулах. С их помощью выражены парные количественные и качественные взаимосвязи: евклидовое расстояние, частота совместной встречаемости, мера сходства множеств градаций признаков. В результате раскрыты те причинно-следственные связи между структурными элементами ландшафта и факторами его формирования, которые объединяют эти элементы в единое зональное целое.

2. Выявлено региональное разнообразие зональных способов организации географической среды и соответствующее многообразие формирования дискретных природно-территориальных образований (региональных ландшафтов) в системе непрерывных гидро-климатических полей. Эти силовые поля создают вещественно-энергетический каркас ландшафтных комплексов, который, взаимодействуя с литогенным каркасом, создает целую серию парциальных ландшафтных структур.

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

4. Полисистемное моделирования природных комплексов ландшафтного и зонального уровней позволяет на основе количественных показателей определить вклад различных факторов в ландшафтную организацию зонального геопространства, оценить его чувствительность к их возможным изменениям. Широкая трактовка оценочных показателей дает возможность усовершенствовать методологию полисистемного расслоения с применением универсальных уравнений обработки пространственных данных.

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

  1. Андреев В.Л. Классификационные построения в экологии и систематике. М.: Наука, 1980. 142 с.

  2. Арманд А.Д. Информационные модели природных комплексов. М.: Наука, 1975. 126 с.

  3. Арманд А.Д. Самоорганизация и саморегулирование географических систем. М.: Наука, 1988. 260 с.

  4. Арманд Д.Л. Наука о ландшафте. М.: Мысль, 1975. 287 с.

  5. Арманд Д.Л., Преображенский В.С., Арманд А.Д. Природные комплексы и современные методы их изучения // Изв. АН СССР. Сер. геогр. 1969. № 5. С. 5–16.

  6. Базилевич Н.И., Гребенщиков О.С., Тишков А.А. Географические закономерности структуры и функционирования экосистем. М.: Наука, 1986. 297 с.

  7. Барталев С.А., Егоров В.А., Ершов Д.В. и др. Спутниковое картографирование растительного покрова России по данным спектрорадиометра MODIS // Совр. пробл. дистанционного зондирования Земли из космоса. 2011. Т. 8. № 4. С. 285–302.

  8. Берг Л.С. Избранные труды. М.: Изд-во АН СССР, 1958. Т. II. Физическая география. 426 с.

  9. Вальтер Г. Растительность Земного Шара. Эколого-физиологическая характеристика. М.: Прогресс, 1974. Т. II. 423 с.

  10. Герасимов И.П. Экологические проблемы в прошлой, настоящей и будущей географии Мира. М.: Наука, 1985. 247 с.

  11. Григорьев А.А. Закономерности строения и развития географической среды. М.: Мысль, 1966. 382 с.

  12. Григорьев А.А. Типы географической среды. Избранные теоретичепские работы. М.: Мысл, 1970. 550 с.

  13. Докучаев В.В. О законности известного географического распределения наземно-растительных почв на территории Европейской России // Тр. С.-Петербург. общ-ва естествоиспытателей. СПб., 1881. Т. XII. С. 65–66.

  14. Докучаев В.В. Учение о природных зонах. М.: Географгиз, 1948. 63 с.

  15. Дюран Б., Оделл П. Кластерный анализ. М.: Статистика, 1977. 128 с.

  16. Жуков В.Т., Сербенюк С.Н., Тикунов В.С. Математико-картографическое моделирование в географии. М.: Мысль, 1980. 224 с.

  17. Зейдис И.М., Кружалин В.И., Симонов Ю.Г. и др. Общие свойства динамики геосистем // Вестн. Моск. ун-та. Сер. 5. География. 2001. № 4. С. 3–8.

  18. Исаченко А.Г. Основы ландшафтоведения и физико-географическое районирование. М.: Высшая школа, 1965. 327 с.

  19. Карта растительности Европейской части СССР. М-б 1 : 2 500 000 / ред. Т.И. Исаченко, Е.М. Лавренко. М.: ГУГК, 1974. 6 л.

  20. Коломыц Э.Г. Бореальный экотон и географическая зональность. Прогнозный Атлас-монография. М.: Наука, 2005. 390 с.

  21. Коломыц Э.Г. Экотон как объект физико-географического исследования // Изв. РАН. Сер. геогр. 1988. № 5. С. 24–36.

  22. Мильков Ф.Н. Физическая география: современное состояние, закономерности, проблемы. Воронеж: Изд-во Воронеж. ун-та, 1981. 398 с.

  23. Николаев В.А. Классификация и мелкомасштабное картографирование ландшафтов. М.: Изд-во Моск. ун-та, 1978. 62 с.

  24. Одум Ю. Основы экологии. М.: Мир, 1975. 740 с.

  25. Преображенский В.С., Александрова Т.Д., Максимова Л.В. География в меняющемся мире. Век ХХ. М.: Ин-т географии РАН, 1997. 273 с.

  26. Пузаченко Ю.Г. Глобальное биологическое разнообразие и его пространственно-временная изменчивость // Современные глобальные изменения природной среды. М.: Научный мир, 2006. Т. 2. С. 306–377.

  27. Пузаченко Ю.Г., Скулкин В.С. Структура растительности лесной зоны СССР: Системный анализ. М.: Наука, 1981. 275 с.

  28. Разнообразие и динамика лесных экосистем России. Кн. 2 / ред. А.С. Исаев. М.: КМК, 2012б. 478 с.

  29. Растительность Европейской части СССР и Кавказа. М-б 1 : 2 000 000 / ред. Т.В. Котова. М.: ГУГК, 1987. 4 л.

  30. Сочава В.Б. География и экология: Материалы 5-го Съезда Географического общества СССР. Л.: ВГО, 1970. 42 с.

  31. Сочава В.Б. Растительный покров на тематических картах. Новосибирск: Наука, 1979. 189 с.

  32. Сочава В.Б. Избранные труды. Теоретическая и прикладная география. Новосибирск: Наука, 2005. 288 с.

  33. Сукачев В.Н. Избранные труды. Л.: Наука, 1972. Т. 1. Основы лесной типологии и биогеоценологии. 418 с.

  34. Хаазе Г. Изучение топических и хорических структур, их динами и развития в ландшафтных системах // Структура, динамика и развитие ландшафтов / отв. ред. В.С. Преображенский. М.: Ин-т географии АН СССР, 1980. С. 57–81.

  35. Хорошев А.В. Рельеф как фактор полимасштабной организации межкомпонентных связей в лесных ландшафтах Восточно-Европейской равнины // Вестн. Моск. ун-та. Сер. 5. География. 2010. № 3. С. 35–42.

  36. Хорошев А.В. Полимасштабная организация географического ландшафта. М.: КМК, 2016. 416 с.

  37. Хорошев А.В., Мерекалова К.А., Алещенко Г.М. Полимасштабная организация межкомпонентных отношений в ландшафте // Изв. РАН. Сер. геогр. 2010. № 1 С. 26–36.

  38. Черкашин А.К. Полисистемное моделирование. Новосибирск: Наука, 2005. 280 с.

  39. Черкашин А.К. Геокартографическое мышление в современной науке // Геодезия и картография. 2020. № 7. С. 27–36.

  40. Forman R.T.T. Landscape Mosaics. Cambridge Univ. Press, 2006. 632 p.

  41. Marceau D.J. The Scale issue in social and natural sciences // Canad. J. Remote Sens. 1999. Vol. 25. № 4. P. 347–356.

  42. Phillips J.D. Global and local factors in earth surface systems // Ecol. Modeling. 2002. Vol. 149. № 3. P. 257–272.

  43. Pickett S.T.F., Cadenasso M.L. The ecosystem as a multidimemsional concept6 meaning, model, and metaphor // Ecosystems. 2002. Vol. 5. P. 1–10.

  44. Turner M., Gardner R.Y., O’Neill R.V. Landscape Ecology in Theory and Practice: Pattern and Process. Springer Verlag, 2001. 352 p.

  45. Wu J., Hobbs R. Key issues and research priorities in landscape ecology. An idiosyncratic synthesis // Landscape Ecol. 2002. Vol. 17. № 4. P. 355–365.

Дополнительные материалы отсутствуют.