УДК 630*181 (182.5)
DOI: 10.33764/2618-981X-2020-3-1-241-247
ОЦЕНКА ПОТЕНЦИАЛЬНОЙ ПРОДУКТИВНОСТИ ДРЕВОСТОЕВ
(НА ПРИМЕРЕ ГОСУДАРСТВЕННОГО ПРИРОДНОГО ЗАПОВЕДНИКА «АЗАС»)
Сергей Кимович Фарбер
Институт леса им. В. Н. Сукачева СО РАН - обособленное подразделение ФИЦ КНЦ СО РАН, 660036, Россия, г. Красноярск, Академгородок, 50/28, доктор сельскохозяйственных наук, ведущий научный сотрудник, тел. (391)249-46-35, e-mail: sfarber@ksc.krasn.ru
Наталья Сергеевна Кузьмик
Институт леса им. В. Н. Сукачева СО РАН - обособленное подразделение ФИЦ КНЦ СО РАН, 660036, Россия, г. Красноярск, Академгородок, 50/28, кандидат сельскохозяйственных наук, ученый секретарь, тел. (391)249-44-37, e-mail: natalia_5791@mail.ru
Нина Иосифовна Молокова
Федеральное государственное бюджетное учреждение «Государственный природный заповедник «Азас», 668530, Россия, Республика Тыва, Тоджинский р-н, с. Тоора-Хем, ул. Агбаан, 20, кандидат биологических наук, ведущий научный сотрудник, тел. (394)502-12-99, e-mail: azasmolokova@mail.ru
Потенциальная продуктивность древостоев рассматривается в зависимости от лесорас-тительных условий. Формируются множественные уравнения регрессии. В качестве зависимой переменной используется класс бонитета. Почва, влага и тепло - основные факторы, влияющие на развитие насаждений. Теплообеспеченность местоположений рассматривается, как производное от показателей рельефа. Показатели влажности и плодородия местоположений конструируются на основе характеристик таксационных выделов. Достоверное снижение класса бонитета по мере набора высоты наблюдается только у лиственничных древосто-ев. Влияние экспозиции и величины уклона на продуктивность древостоев не прослеживается. Ранжированный ряд влажности местоположений демонстрируется относительно обобщенного качественного показателя лесорастительных условий - типа леса. Наблюдается наличие статистически значимой связи продуктивности древостоев с типами почв.
Ключевые слова: рельеф, продуктивность насаждений, типы лесов и нелесных сообществ, ландшафт, картографирование.
EVALUATION OF POTENTIAL PRODUCTIVITY OF TREES
(ON THE EXAMPLE OF THE STATE NATURAL RESERVE «AZAS»)
Sergey K. Farber
Federal Research Center Krasnoyarsk Scientific Center, Russian Academy of Sciences, Siberian Branch V. N. Sukachev Institute of Forest, Russian Academy of Sciences, Siberian Branch, 50/28, Akademgorodok, Krasnoyarsk, 660036, Russia, D. Sc., Leading Researcher, phone: (391)249-46-35, e-mail: sfarber@ksc.krasn.ru
Natalia S. Kuzmik
Federal Research Center Krasnoyarsk Scientific Center, Russian Academy of Sciences, Siberian Branch V. N. Sukachev Institute of Forest, Russian Academy of Sciences, Siberian Branch, 50/28, Akademgorodok, Krasnoyarsk, 660036, Russia, Ph. D., Scientific Secretary of the Academic Council, phone: (391)249-44-37, e-mail: natalia_5791@mail.ru
Nina I. Molokova
State Wildlife Reserve «Azas», 20, Agbaan, Toora-Khem, 668530, Russia, Tyva Republic, Ph. D., Leading Researcher, phone: (394-50)2-12-99, e-mail: azasmolokova@mail.ru
The potential productivity of forest stands is considered depending on forest conditions. Multiple regression equations are formed. As a dependent variable, the bonitet class is used. Soil, moisture and heat are the main factors affecting the development of plantings. The heat supply of locations is considered as a derivative of the topography. The moisture and fertility indicators of locations are constructed based on the characteristics of taxation stands. A significant decrease in the bonitet class as climbing is observed only in larch stands. The influence of exposure and slope on the productivity of stands is not traced. A ranked range of location moisture is demonstrated relative to a generalized qualitative indicator of forest conditions such as forest. There is a statistically significant relationship between the productivity of the stands and soil types.
Key words: relief, stand productivity, types of forests and non-forest communities, landscape, mapping.
Введение
Государственный природный заповедник «Азас» расположен в республике Тыва. По лесорастительному районированию территория заповедника входит в Тоджинский округ подтаежных сосново-лиственничных и горнотаежных лиственничных, сосновых и кедровых лесов [1]. Лесоустройство в заповеднике проводилось в 2015 г. по 3-му разряду. Количество таксационных выделов в заповеднике 11739. Из них сосновых - 56, лиственничных -2424, кедровых - 5931, еловых - 237. Анализ выполнен с привлечением всех таксационных выделов заповедника. Кроме данных массовой таксации и ЦМР в работе использовались нормативно-справочные материалы лесоустройства -схемы типов леса и нелесных сообществ, планово-картографические и фондовые материалы заповедника.
Продуктивность древостоев определяется посредством многофакторного регрессионного анализа. В качестве зависимой переменной принимается класс бонитета. В качестве аргументов - показатели лесорастительных условий. Почва, влага и тепло - основные факторы, влияющие на рост и состояние насаждений. Уровень значимости уравнений принят равным 0,95. Далее для каждого таксационного выдела производился расчет продуктивности древосто-ев. Результаты были внесены в соответствующее поле атрибутивной таблицы ГИС, что определяет возможность последующего картографирования.
Цель работы - определение и картографирование потенциальной продуктивности (классов бонитета) древостоев, в т. ч. на не покрытых лесом землях.
Тепло. Температурный режим местоположений находится в зависимости от показателей рельефа. Отсюда следует, что при выявлении сопряженности характеристик лесного покрова с лесорастительными условиями количество тепла можно заменить показателями рельефа местности [2].
При наличии связи формируется уравнение регрессии
,
где Вг - класс бонитета, зависящий от рельефа(количества тепла). Чтобы избежать необходимости формирования отдельных зависимостей для северных и южных экспозиций румбам симметричным относительно оси север - юг присваивается одноименный индекс экспозиции.
Методы и материалы
Анализ проводился средствами ArcGISSpatialAnalyst. Для получения показателей рельефа использована ЦМР SRTM. Принята следующая градация показателей рельефа:
- H (абсолютная высота) - интервал 200 м;
- Asp (экспозиция) - интервал 45°;
- Sl (уклон) - интервал 2°.
Рельеф заповедника горный. Однако вопреки ожиданиям корреляция классов бонитета древостоев с показателями рельефа в основном отсутствует. Положительная корреляция с абсолютной высотой наблюдается только у лиственничных древостоев R=0,402. Влияние экспозиции и уклонов на продуктивность древостоев не прослеживается. Влияние рельефа (тепла) на продуктивность лиственничных древостоев описывается уравнением:
3, = 3.77-:- 0.027 ■ И-- 0.290 ■ = 0.170>.
Влага. Оценка степени влажности местоположений (почв) производится по описаниям лесотаксационных выделов. Оценка комплексная, учитывающая в т. ч. информацию из макета дополнительных сведений. По отношению к типам леса принимается последовательное увеличение от сухих к более увлажненным. Порядковый номер (индекс) в ранжированном ряду отождествляется со степенью влажности местоположений (почв) [2]. При наличии связи формируется уравнение регрессии. Общий вид зависимости Ви = f{JAО, где W порядковый номер (индекс) в ранжированном ряду влажности почв, Вw - класс бонитета, зависящий от влажности местоположений (почв).
Сформированы уравнения:
- сосняки -
з,., = -..б:: - о.з-и - - o.osi);
- лиственничники -
£Fff = 4,064 + 0079 * W2 - 0,568* Wf (Л2 = 0,418);
- кедровники -
Bw = 5,706+ 0,050 * W - 0,460 * VK2, (Jt2 = 0,270);
- ельники -
£?. = з.з:*- 0.030- vr-.-jr- = С'-оОл
Теснота связи влажности почв с продуктивностью лиственничных, кедровых и еловых древостоев умеренная, с продуктивностью сосновых - слабая. Насаждения заповедника произрастают при определенной влажности. Сосняки занимают местоположения с баллами влажности 2 и 3; лиственничники от 1 до 9; кедровники от 1 до 9; ельники от 3 до 9. Еще более влажные местоположения занимают заболоченные кустарники. Получается, что влажность почв лимитирует местоположение пород деревьев, и этот факт далее используется при картографировании потенциальной продуктивности древесных пород. Например, местоположения с индексом влажности 1 в качестве потенциально возможных для произрастания ели выпадают.
Плодородие почв. Принимается, что класс бонитета древостоя производное лесорастительных условий, в перечень которых входит в т. ч. плодородие лесных почв. Можно записать Bs = /Т5), где: Bs - класс бонитета древостоя, зависящий от плодородия почв; S - показатель лесорастительных условий. Класс бонитета древостоя и наименование (тип) почв фигурируют в описании таксационного выдела. Для получения зависимости^ = /(^остается определится с конструкцией показателя лесорастительных условий S.
Возможность сравнения появляется при использовании заранее оговоренной линии отсчёта. В качестве показателя лесорастительных условий ранее было предложено отношение i£/dö, где ä9 = /(A}; d, h - диаметр и высота главной породы древостоя [3]. Вообще же преимущество за показателем, варьирующим менее других. Поэтому в качестве измерителя лесорастительных условий лучше использовать не отношение диаметров iÉ/f£a, а отношение высот h/hQ. Тогда в качестве линии отсчета будет выступать усредненная линия роста по высоте (уравнение регрессии = /(А), где А - возраст древостоя, лет). Посредством введения дополнительного аргумента d точность оценок ho повышается. Тогда = /{Д,А), при этом началом отсчета будет уже не линия, а плоскость.
Сопоставление плодородия почв возможно относительно отдельных древесных пород, в нашем случае относительно уравнений h0 = fid, АX которые индивидуальны для сосны лиственницы кедра и ели. Для определенной древесной породы и типа почвы среднее значение показателя лесорастительных условий = 2 7я гДе п ~ число таксационных выделов. Общее для всех пород деревьев и типа почвы средневзвешенное значение показателя лесорастительных условий
где iV — 2 Щ ■
Сформированы уравнения:
- сосняки -
ft0 = 0Д32+ 1,084 * d - 0,010* d2 — 0,033*4, (R2 = 0,856);
- лиственничники -
A0 = 2,053+ 1,373 * d- 0,014 *d2- 0,104*4 № = 0,?S9>;
- кедровники -
= 0,04?+ 0,809* d - 0,008* d2 + 0,014 * Д (ft2 = 0,705);
- ельники -
ha = —2,791 + 1Д30 * d - 0,011 * d2 + 0,012 * А (tf3 = 0,858).
На основе уравнений ti$ = / (f£, дня каждого таксационного выдела рассчитаны значения S = Далее получены показатели лесорастительных условий Si и S0. Оценка плодородия почв (и картирование) на основе Si возможна только для таксационных выделов i - той древесной породы. Показатель Se относится ко всем древесным породам. Соответственно оценка плодородия почв (и картирование) на основе Se уже возможна только на основании наименования почвы для всех таксационных выделов, включая выдела лиственных древостоев (березняки и осинники), выдела погибших насаждений, вырубки.
Аналитический вид уравнений Bs = f СУ)получен по данным классов бонитета древостоев и показателям условий произрастания Se таксационных выделов:
- сосняки - связи нет;
- лиственничники -
Bs = -7&580- 105141 + 18?,55&*5г,(й2 = 0,443);
- кедровники -
Bs = 4,822 - 8,220*+ 8,375*£0, (Д2 = 0,226);
- ельники -
Bs = 14 JS9- 9,624- SS,[Д3 = 0,335).
Результаты
Потенциальная продуктивность древостоев. Формируются многофакторные уравнения регрессии. В качестве предикторов можно использовать значения Вг, Вм„ Bs таксационных выделов. Тогда общий вид зависимости В$ = ff_EriBwlBs}, где Во - класс бонитета, зависящий от тепла, влажности и плодородияпочв. Количество вычислительной работы существенно сокращается если в качестве предикторов используются непосредственно показатели рельефа и почв (для тестового участка Д, =
Сформированы уравнения:
- сосняки -
Bq = 4,622 - 0,311 * W, (К2 = 0,081);
- лиственничники -
В9 = -63,303 - 78,929 * Si + 0,052 * W2 - 0,002 * Н - 0,419 * W Н-
+ 0,006* - 0,520);
- кедровники -
Bq = 20,339 +12,385 * 5} + 0,033 * W2 - 0,286 * W - 2743 2 * SQt [JR2 =
- ельники -
Rö =3,528-1 0,030 * W2, (fi2 = 0,480).
Картографирование потенциальной продуктивности произведено по значениям полей атрибутивной таблицы ГИС абсолютная высота Н, влажность почв Ж и плодородие почв Бо (таблица).
Данные атрибутивной таблицы ГИС (фрагмент)
Показатели таксационного выдела Лесорастительные условия, индекс Потенциальная продуктивность древостоев Во, класс бонитета
Тип леса A, лет h, м d, см Класс бонитета H W So Сосна Лиственница Кедр Ель
Лбрзм 170 20 28 4 9 3 1.03 4.06 3.03 4.71 3.79
Кбрзм 180 19 28 4 10 3 1.03 4.06 3.14 4.71 3.79
Кбрзм 180 20 30 4 10 3 1.03 4.06 3.14 4.71 3.79
Лбрзм 170 23 30 3 10 3 1.03 4.06 3.14 4.71 3.79
Лртбр 180 22 30 4 10 3 1.01 4.06 3.44 4.76 3.79
Лртбр 170 23 30 3 10 3 1.01 4.06 3.44 4.76 3.79
Кртбр 10 1 4 10 3 1.01 4.06 3.44 4.76 3.79
Лбрзм 170 22 28 4 10 3 1.03 4.06 3.14 4.71 3.79
Кглмх 180 17 30 5 10 7 0.97 4.25 5.05 4.98
Кермх 190 20 36 4 10 7 0.97 4.25 5.05 4.98
Лермх 180 21 28 4 10 7 0.97 4.25 5.05 4.98
Кбрзм 200 20 32 4 10 3 1.03 4.06 3.14 4.71 3.79
Кглмх 190 16 26 5 11 7 0.97 4.37 5.05 4.98
БКУл- 10 1 5А 10
Лглмх 180 22 28 4 11 7 0.97 4.37 5.05 4.98
Лермх 160 16 24 5 11 7 0.97 4.37 5.05 4.98
Кбрзм 200 19 32 4 11 3 1.03 4.06 3.26 4.71 3.79
БКУм 10 1 5А 9 0.96
Еермх 110 16 20 5 9 7 0.97 4.15 5.05 4.98
Условные обозначения: Лбрзм - лиственничники бруснично-зеленомошные; Кбрзм -кедровники бруснично-зеленомошные; Лртбр - лиственничники разнотравно-брусничные; Кртбр - кедровники разнотравно-брусничные; Кглмх - кедровники голубнично-моховые; Кермх - кедровники ерниково-моховые; Лермх - лиственничники ерниково-моховые; БКУл-гтр-береза кустарниковая лугово-разнотравная; БКУмхбл - береза кустарниковая мохово-болотная; Еермх - ельники ерниково-моховые.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Типы лесов гор Южной Сибири. - Новосибирск: Наука, 1980. - 334 с.
2. Фарбер С. К., Кошкарова В. Л., Кузьмик Н. С. Картографирование лесных формаций голоцена с использованием основных показателей климата - тепла и влаги // Сиб. лесн. журн. 2017. № 6. - С. 26-40.
3. Фарбер С. К. Лесные измерения по среднемасштабным аэроснимкам. Красноярск: Изд-во СОРАН, 1997. -106 с.
© С. К. Фарбер, Н. С. Кузьмик, Н. И. Молокова, 2020