книги2 / konf_15-24
.pdf
Электронный архив УГЛТУ
Научная статья УДК 630*521
МОДЕЛЬ ОБРАЗУЮЩЕЙ СТВОЛОВ ДЕРЕВЬЕВ СОСНЫ КОСТРОМСКОЙ ОБЛАСТИ
Владимир Викторович Гостев1, Александр Вячеславович Лебедев2
1, 2 Российский государственный аграрный университет – МСХА имени К. А. Тимирязева, Москва, Россия
1v.gostev@rgau-msha.ru
2alebedev@rgau-msha.ru
Аннотация. Изменение диаметра дерева на любой высоте можно аппроксимировать с помощью уравнений образующей древесного ствола. Рассматривается четырехпараметрическая функция образующей древесного ствола, апробированная на региональных данных, наиболее корректно аппроксимирующая изменение диаметра с высотой в сосновых древостоях Костромской области.
Ключевые слова: форма ствола, образующая ствола, регрессионные модели, сосновые древостои, Костромская область
Благодарности: исследование выполнено за счет гранта Российского научного фонда № 23-76-01016, https://rscf.ru/project/23-76-01016/.
Original article
MODEL OF THE FORMAT OF PINE TREE TRUNKS
IN THE KOSTROMA REGION
Vladimir V. Gostev1, Aleksandr V. Lebedev2
1, 2 Russian State Agrarian University – Moscow Timiryazev Agricultural Academy, Moscow, Russia
1v.gostev@rgau-msha.ru
2alebedev@rgau-msha.ru
Abstract. The change in tree diameter at any height can be approximated using the equations of the tree stem taper. A four-parameter function of the stem taper is considered, tested on regional data, which most correctly approximates the change in diameter with height in pine stands of the Kostroma region.
Keywords: stem shape, stem taper, regression models, pine stands, Kostroma region
© Гостев В. В., Лебедев А. В., 2024
71
Электронный архив УГЛТУ
Acknowledgments: this study has been supported by the grants the Russian Science Foundation, RSF 23-76-01016, https://rscf.ru/en/project/23-76-01016/.
Модели образующей древесного ствола получили широкое распространение в лесохозяйственных исследованиях [1–3]. Уравнения образующей древесных стволов являются надежной основой для выявления сортиментной структуры стволов и древостоев. Для Костромской области ранее не производились работы, направленные на моделирование образующей стволов [4–6]. Поэтому целью исследования являлась разработка модели образующей древесного ствола с одним уравнением, наиболее корректно описывающей изменение диаметра с высотой в древостоях сосны обыкновенной, произрастающих в Костромской области.
В качестве экспериментальных данных использовались данные 10064 измерений диаметров для 692 деревьев сосны. Диапазон диаметров на высоте 1,3 м составляет от 3,9 до 54,1 см, высот – от 5,6 до 34,5 м. Возраст деревьев составил от 10 до 148 лет. Полученные данные характеризуют господствующие и согосподствующие деревья из полога насаждения. Расчеты производились в среде Python, статистические выводы сделаны на уровне значимости 5 %.
Наиболее адекватным уравнением образующей стволов деревьев сосны в условиях Костромской области является четырехпараметрическая модель [8], имеющая следующий вид:
|
|
|
|
|
|
|
|
hi |
|
0,25 |
|
|
hi |
|
0,5 |
|
D |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
b |
b |
|
b |
|
b |
H |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
H h |
|
0 |
1 |
H |
|
2 |
H |
|
3 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
di |
D |
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
H 1,3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где di – диаметр дерева на высоте hi |
; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
D – диаметр дерева на высоте 1,3 м;– высота дерева, м;– параметры модели.
Для полученной модели рассчитаны критерии достоверности [7]: квадратный корень из среднеквадратической ошибки (RMSE = 1,497), средний абсолютный процент ошибки (PE = 8,655), средняя абсолютная ошибка (MAE = 0,926), среднее смещение остатков (MBE = 0,051), коэффициент детерминации (R2 = 0,982). Анализ распределения остатков и метрик качества показывает, что модель характеризуется высокой обобщающей способностью.
На рис. 1 представлена прогнозируемая по уравнению образующая стволов сосны обыкновенной, построенная по предсказанным значениям диаметров, в сравнении с фактическими для дерева средней величины из
72
Электронный архив УГЛТУ
выборки с таксационным диаметром 21,9 см и высотой 23,0 м. Модель адекватно описывает образующую комлевой, центральной частей ствола и вершины.
Предсказанные уравнением значения диаметров были сопоставлены с данными из таблиц сбега древесных стволов по разрядам высот [9]. На рис. 2 представлен график зависимости диаметров из таблиц сбега и рассчитанных с помощью модели по разрядам высот. На рисунке можно отметить занижение таблицами сбега значений диаметров по сравнению с данными модели для 2–4 разрядов высот, что в конечном итоге приводит к занижению объемов стволов и их частей. Предложенная модель образующей получена на основе региональных данных и, учитывая местные особенности древостоев, точнее предсказывает значения диаметров на любой высоте.
|
30 |
|
|
|
|
|
, см |
25 |
|
|
|
|
|
|
|
|
|
|
|
|
в коре |
20 |
|
|
|
|
|
15 |
|
|
|
|
|
|
Диаметр |
|
|
|
|
|
|
10 |
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
0 |
5 |
10 |
15 |
20 |
25 |
|
|
|
Высота, м |
|
|
|
|
|
|
Исходные данные |
Модель 14 |
|
|
Рис. 1. Кривая образующей ствола, аппроксимированная при помощи модели |
||||||
|
и фактические значения диаметров для дерева средней величины |
|
||||
|
|
|
(D = 21,9 см, H = 23,0 м) |
|
|
|
Диаметр по модели, см
80 |
|
|
70 |
|
|
60 |
1б |
|
|
||
50 |
1а |
|
|
||
40 |
1 |
|
2 |
||
|
||
30 |
3 |
|
|
||
20 |
4 |
|
|
||
10 |
5 |
0 |
|
|
|
|
|
|
|
5a |
0 |
10 |
20 |
30 |
40 |
50 |
60 |
70 |
80 |
|
|
Диаметр по таблицам сбега, см |
|
|
||||
Рис. 2. Сопоставление диаметров из таблиц сбега c диаметрами, рассчитанными по модели, по разрядам высот
73
Электронный архив УГЛТУ
Полученная модель образующей может использоваться для оценки сортиментной структуры отдельных стволов и древостоев в сочетании с моделями распределения деревьев по толщине. Разработанная модель позволяет вычислять различные виды сбега ствола [10] и его отдельных частей. Требуются дальнейшие исследования, направленные на разработку моделей сортиментной структуры стволов, основывающихся на предложенном уравнении образующей, что позволит получить качественные региональные таксационные нормативы.
Список источников
1.Забавская Л. Н., Вайс А. А. Параметры образующей функции “Harris” и форма нижней части деревьев сосны // Хвойные бореальной зоны. 2021. № 39 (2). С. 95–101.
2.Батурин К. В., Стариков А. В. Анализ математических моделей образующих древесных стволов // Актуальные направления научных исследований XXI века: теория и практика. 2015. № 3 (5). С. 23–26. DOI: 10.12737/16198
3.Kozak A. My last words on taper equations // Forestry Chronicle. 2004.
№80. P. 507–515. DOI: 10.5558/tfc80507-4
4.Кузьмичев В. В., Лебедев А. В. Закономерности изменения размеров и качества древесины деревьев в лесах Европейской России (по материалам А. А. Крюденера). Кологрив : Государственный природный заповедник «Кологривский лес», 2022. 96 с.
5.Кузьмичев В. В., Лебедев А. В., Гостев В. В. Образующая стволов деревьев сосны в таблицах А. А. Крюденера // Материалы V Национальной конференции по итогам научной и производственной работы преподавателей и студентов в области лесного дела, ландшафтной архитектуры, мелиорации и экологии, посвященной 100-летию со дня рождения профессора М. А. Дудорева : Сборник материалов конференции (Саратов, 15–19 мая 2023 г.). Саратов : ООО «Амирит», 2023. С. 153–156.
6.Петровский В. С., Малышев В. В., Мурзинов Ю. В. Моделирование параметров древесных стволов в насаждении // Лесотехнический журнал.
2012. № 4 (8). С. 18–22.
7.Дубенок Н. Н., Лебедев А. В., Гостев В. В. Регрессионные модели смешанных эффектов зависимости высоты от диаметра ствола в сосновых древостоях европейской части России // Лесной вестник – Forestry Bulletin. 2023. Т. 27, № 5. С. 37–47. DOI: 10.18698/2542-1468-2023-5-37-47
8.Zeng W. S., Liao Z. Y. A study on taper equation // Scientia Silvae Sinica. 1997. № 33. P. 127–132.
9.Общесоюзные нормативы для таксации лесов : справочник / [Загреев В. В. и др.]. Москва : Колос, 1992. 494 с.
10.Сальникова И. С., Анчугова Г. В., Нагимов З. Я. Таксация леса : учебное пособие. Екатеринбург : УГЛТУ, 2017. 72 с.
74
Электронный архив УГЛТУ
Научная статья УДК 630.568
ФИТОМАССА СОСНЫ СИБИРСКОЙ (PINUS SIBIRICA DU TOUR) НА ВЕРХНЕЙ ГРАНИЦЕ ЛЕСА В УСЛОВИЯХ АЛТАЕ-САЯНСКОГО ГОРНО-ТАЕЖНОГО РАЙОНА
Антон Максимович Громов1, Павел Александрович Моисеев2, Андрей Андреевич Григорьев3, Дмитрий Сергеевич Балакин4
1, 4 Уральский государственный лесотехнический университет, Екатеринбург, Россия 2, 3 Институт экологии растений и животных УрО РАН, Екатеринбург, Россия
1heytonny@yandex.ru
2moiseev@ipae.uran.ru
3grigoriev.a.a@ipae.uran.ru
4dmitrijbalakin047@gmail.com
Аннотация. В XX веке, в связи с наблюдениями д-ра биол. наук С. Г. Шиятова об оценке и анализе климатических функций отклика древесных растений, стало понятно, что леса играют одну из важнейших ролей в стабилизации климата по средствам поглощения углерода из атмосферы и с последующим депонированием его в фитомассе. А. А. Григорьев в своей работе по изучению Южного Урала, основываясь на количественной оценке трансформации высокогорных экосистем, показал, что потепление климата обусловливает экспансию древесно-кустарниковой растительности в горные тундры и луга, а также увеличение густоты и продуктивности ранее произраставших древостоев. Основываясь на фактах о глобальном изменении климата и о данных по продвижению растительности в горные тундры, изучение взаимосвязей фитомассы древостоев с климатическими условиями приобретает особую актуальность.
Впервые произведены расчеты запаса фитомассы на верхнем пределе произрастания древостоя Горного Алтая и Западного Саяна, а также произведена оценка запасов по структурным частям фитомассы деревьев.
Ключевые слова: горная тундра, редина, редколесье, фитомасса, Сосна сибирская (Pinus sibirica Du Tour)
© Громов А. М., Моисеев П. А., Григорьев А. А., Балакин Д. С., 2024
75
Электронный архив УГЛТУ
Original article
PHYTOMASS OF SIBERIAN PINE (PINUS SIBIRICA DU TOUR) ON THE UPPER BORDER OF THE FOREST IN THE ALTAI-SAYAN MOUNTAIN TAIGA REGION
Anton M. Gromov1, Pavel A. Moiseev2, Andrey A. Grigoriev3, Dmitry S. Balakin4
1, 4 Ural State Forest Engineering University, Yekaterinburg, Russia
2, 3 Institute of Plant and Animal Ecology, Ural Branch of the Russian Academy of Sciences, Yekaterinburg, Russia
1heytonny@yandex.ru
2moiseev@ipae.uran.ru
3grigoriev.a.a@ipae.uran.ru
4dmitrijbalakin047@gmail.com
Abstract. In the XX century, in connection with the observations of Doctor of Biological Sciences S. G. Shiyatov on the assessment and analysis of the climatic response functions of woody plants, it became clear that forests play one of the most important roles in climate stabilization by means of carbon absorption from the atmosphere and subsequent deposition in phytomass. Grigoriev, in his work on the study of the Southern Urals, based on a quantitative assessment of the transformation of high-altitude ecosystems, showed that climate warming causes the expansion of tree and shrub vegetation into mountain tundra and meadows, as well as an increase in the density and productivity of previously growing stands. Based on the facts about global climate change and data on the promotion of vegetation in the mountain tundra, the study of the relationship of phytomass of stands with climatic conditions is of particular relevance.
For the first time, calculations of the phytomass reserve at the upper limit of the growth of the stand of the Altai Mountains and Western Sayan were made, as well as an assessment of reserves for structural parts of the phytomass of trees.
Keywords: mountain tundra, redina, woodlands, phytomass, Siberian pine (Pinus sibirica Du Tour)
Район исследований – в пределах Алтае-Саянского горно-таежного района, а именно хребет Холодный белок (50°11-44' с. ш. и 83°15-40" в. д. 2250 м над ур. м) и гора Кохош (51°45-15' с. ш. и 89°43-14" в. д. 2524 м над ур. м), расположенные в пределах центральных возвышенностей горной цепи (рис. 1).
76
Электронный архив УГЛТУ
Рис. 1. Район исследования – хр. Холодный белок (Горный Алтай), г. Кохош (Западный Саян)
На склонах восточной экспозиции хребта Холодный белок и горы Кохош в переходной зоне лес-горная тундра летом 2021–2022 гг. заложены высотные профили. В пределах каждого профиля зафиксированы четыре высотных уровня (ВУ) (рис. 2).
Рис. 2. Район заложения высотных профилей:
а – (гонный Алтай) хр. Холодный белок, б – (Западный Саян) г. Кохош
Верхний уровень (1 ур.) представлен еденичными деревьями; (2 ур). – у верхней границы редин, средний уровень; (3 ур.) – у границы редколесий, нижний; (4 ур.) – у верхней границы сомкнутых лесов. На протяжении всех высотных уровней были заложены от 4 до 6 пробных площадок. Средняя площадь пробной площади составляет 200 м2 [1].
После заложения площадок были произведены работы по замеру таксационных показателей по каждому дереву. Определялись: высота дерева, диаметр ствола у основания и на высоте груди (1,3 м), диаметр проекции
77
Электронный архив УГЛТУ
кроны по двум направлениям, жизненное состояние дерева или оценка жизнеспособности. Для определения возраста дерева, с помощью возрастного бурава, из деревьев, диаметром превышающих 3 см у основания, были извлечены буровые образцы древесины (керны). В лабораторных условиях каждый полученный образец приклеивался на деревянную основу, и для точного определения возраста производилась зачистка лезвием. Для определения ширины годичных колец все образцы (керны) были измерены установкой LINTABVI (Frank Rinn-Rinntech, Германия). В измерении датировки самого ближнего к центру кольца использовалось программное обес-
печение TSAP−3.0.
Для изучения структуры и запасов наземной фитомассы в летний период был собран материал на переходной зоне лес-горная тундра. Отбор деревьев происходил по ступеням толщины (шаг ступени 4 см) с предварительным определением размеров. Такой принцип отбора наиболее подходящий для исследования зависимостей массы различных частей деревьев от их размерных показателей (диаметра, высоты и др.). Так как программа исследования подразумевает определение общей надземной фитомассы и отдельных ее компонентов, то в учет принимались не только такие показатели, как диаметр у основания, на высоте груди и т. д, но и происходило разделение на составляющие фракции, а именно ствол в коре, ветви, хвоя, корни и генеративные органы [2].
При измерении фитомассы ствола использовались электронные весы (безмен) с точностью ±50 г. Также происходил отбор дисков на высоте каждого метра от высоты ствола, включая основание и высоту груди. После выпиливания производилось взвешивание дисков в коре и без коры в полевой лаборатории. При дальнейшем исследовании диски проходили высушивание в сушильном шкафу (ШСП−0,25−100, Россия) при температуре 105 ºC до постоянного веса. Это позволило переделить содержание абсолютно сухой массы фракции ствола путем взвешивания в сырой навеске и после сушки. Были отобраны и обработаны в двух разных районах суммарно 18 модельных деревьев сосны сибирской (Pinus sibirica Du Tour) по изложенной ранее методике.
На основе результатов, полученных при изучении заложенных пробных площадей, создана характеристика древостоев по высотным уровням. Данные представлены в табл. ниже.
Преобладающей породой в составе на нижних уровнях является сосна сибирская (Pinus sibirica Du Tour). На четвертых высотных уровнях различия в возрасте минимальны, но уже на третьем уровне породный процент различия составляет 35 %. На первом же уровне он составил 63 %.
Различия наблюдаются и в диаметре на высоте груди (1,3 м). По мере продвижения вниз по склону визуально можно увидеть существенное улучшение микроклиматических условий. Происходит резкое увеличение количества деревьев, в среднем в 8 раз на вторых высотных уровнях и на третьих – в 25 раз.
78
Электронный архив УГЛТУ
Можно заключить, что в зависимости от высоты над уровнем моря наблюдаются морфометрические различия лесообразующих пород.
Как ранее и утверждалось, объем деревьев и их структурных частей зависит от огромного количества различных факторов. В таких исследованиях важное значение имеет обоснованный выбор факторов (таксационных показателей). Именно они объясняют максимальную долю изменчивости массы стволов и крон. При этом существует множество различных методов определения корреляции фитомассы структурных частей дерева. Одной из таких методик является нахождение фитомассы путем расчета многофакторных уравнений. Отмечается, что в уравнениях парной связи целесообразно использовать показатель D2H. Именно этим методом сделан акцент в данной работе [3].
Большинство ученых, проводивших изучение этой тематики, признают, что для выражения зависимостей наиболее обоснованным является применение степенной (алгоритмической) функции, статистические параметры которой корректно объясняются с биологической позиции (рис. 3) [4].
Рис. 3. Зависимость общей фитомассы стволов модельных деревьев от диаметра
Благодаря разработанным уравнениям зависимости и материалам распределения деревьев, на пробных площадках был произведен расчет запаса фитомассы на единицу площади (рис. 4). Для полной картины состава древостоя на склонах была добавлена Лиственница сибирская (Larix sibirica ladeb).
Рис. 4. Запас фитомассы структурных частей по высотным уровням
79
Электронный архив УГЛТУ
Характеристика древостоев на исследуемых профилях, разбитых на высотные уровни
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
|
Профиль уровеньи |
Экспозиция |
Порода |
Диаметр 1,3 м, |
Высота, м |
Возраст, лет |
Диаметр кроны, м |
Проективное покрытие м,кроны2/га |
шт,Отпад./га |
Плотность |
шт./га |
|||||
|
|
|
|
см |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Средн. |
|
Макс |
Средн. |
Макс |
Средн. |
Макс |
Средн. |
Макс |
|
|
<1,5 |
>1,5 |
|
|
|
|
|
|
|
м |
м |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
Горный Алтай, хр. Холодный белок |
|
|
|
|
|
||||
|
1 |
E |
Ps |
5,8±2,9 |
|
13,0 |
2,4±0,7 |
3,7 |
28±18 |
96 |
1,7±0,5 |
2,9 |
500 |
33 |
470 |
191 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
E |
Ps |
7,4±5,2 |
|
22,0 |
3,6±1,5 |
7,3 |
38±22 |
116 |
2,0±0,9 |
5,7 |
1046 |
33 |
154 |
529 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
E |
Ps |
12,7±8,2 |
|
34,1 |
5,5±2,7 |
13,0 |
48±18 |
104 |
2,8±1,3 |
6,1 |
3001 |
11 |
55 |
760 |
|
4 |
E |
Ps |
15,3±8,8 |
|
42,0 |
9,8±4,5 |
20,9 |
91±19 |
137 |
3,4±1,5 |
6,8 |
7610 |
143 |
0 |
881 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
80 |
|
|
|
|
|
|
|
Западные Саяны, хр. Кохош |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
E |
Ps |
3,1±1,6 |
|
8,3 |
2,3±0,8 |
6,3 |
– |
– |
0,9±0,4 |
1,7 |
467 |
97 |
678 |
687 |
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
E |
Ps |
10,1±7,5 |
|
37,0 |
5,3±2,2 |
9,6 |
– |
– |
1,9±1,4 |
5,5 |
3348 |
44 |
70 |
987 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
E |
Ps |
19,3±9,1 |
|
43,0 |
8,3±2,9 |
14,0 |
– |
– |
3,6±1,5 |
6,2 |
8987 |
44 |
132 |
907 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
E |
Ps |
17,5±19 |
|
90,0 |
7,4±4,7 |
20,0 |
– |
– |
2,6±1,8 |
8,3 |
6740 |
62 |
343 |
890 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Примечание. Ps – Сосна сибирская (Pinus sibirica Du Tour), E – восточная экспозиция.
80