CC BY
73
УДК 581.526
ИМИТАЦИОННАЯ МОДЕЛЬ ДИНАМИКИ ФИТОМАССЫ ДЕРЕВЬЕВ СОСНЫ ОБЫКНОВЕННОЙ
А. И. Старцев
Нижегородская государственная сельскохозяйственная академия, 603107, Российская Федерация, Нижний Новгород, пр. Гагарина, 97 E-mail: astartsev2005@rambler.ru
Разработана имитационная модель роста деревьев сосны обыкновенной. Модель имеет девять фазовых переменных: масса (объём) ствола, фитомасса ветвей, хвои, тонких, толстых корней, коры, валовый годичный фотосинтез, фонд запасных углеводов и высота. Скорость роста дерева зависит от общего количества углеводов, синтезируемых растением в течение вегетационного периода. Разработанная модель позволяет вычислять биометрические показатели и фитомассу отдельных компонентов деревьев в возрасте от двух до 120-150 лет.
Ключевые слова: модель роста дерева; фитомасса ствола, хвои, ветвей, корней; валовый годичный фотосинтез; фонд запасных углеводов.
Введение. Разработка моделей роста и динамики фитомассы древостоев является актуальной задачей для решения проблемы оценки биологической продуктивности лесных экосистем. Для её решения широко используются различные математические модели, наиболее перспективными из которых являются имитационные модели роста древостоев и отдельных деревьев.
Модель О. Г. Чертова [1] предназначена для описания экосистемы одного растения (ЭОР). Под ней понимается пространство, занимаемое надземной частью высшего растения, и объём почвы, занятый его корневой системой. Лимитирующими факторами модели приняты фото-синтетически активная радиация (ФАР), вода и элементы питания.
В модели В. В. Бугровского и др. [2, 3] прирост биомассы определяется разностью между фотосинтезом и расходами на дыхание с учётом опада. В отличие от других моделей полог леса разбит на пять равных по массе слоёв листьев. Каждый слой характеризуется своей кривой фотосинтеза, наклон, высота плато и точка
компенсации которой зависят от условий освещённости. Модель предназначена для расчёта годичной продуктивности насаждений дуба черешчатого.
Наиболее полной экофизиологиче-ской моделью является модель РЦЦ-1 К. Кулля и др. [4] для свободно растущего дерева ели. Она содержит девять фазовых переменных (масса хвои, ствола, тонких и толстых корней, фонды углеводов в листьях, стволе и корнях, фонд азота и высота дерева). Модель позволяет рассчитывать сезонную (по декадам) и многолетнюю динамику прироста, однако не учитывает процессы конкуренции в древостое за свет и элементы питания.
Для деревьев дуба черешчатого, растущих в чистых одновозрастных древо-стоях, разработана модель роста с учётом взаимозатенения деревьев, влияния половодий, осадков, испарения и транспира-ции на содержание влаги в почве, влияние температур, влажности воздуха, продолжительности солнечного сияния и давления влаги в почве на интенсивность фотосинтеза [5]. Модель содержит 15 фазовых переменных и описывает динамику при-
© Старцев А. И., 2016.
Для цитирования: Старцев А. И. Имитационная модель динамики фитомассы деревьев сосны обыкновенной // Вестник Поволжского государственного технологического университета. Сер.: Лес. Экология. Природопользование. 2016. № 1 (29). С. 30-41.
роста биомассы ствола, листьев, ветвей, толстых и тонких корней, фонда азота и углеводов, фонда запасаемых углеводов.
Модели для распространения осиновых древостоев составлялись с использованием коэффициентов расселения и плотности популяции [6].
В работе Ю.П. Демакова и др. [7] использовались аллометрические зависимости для оценки фитомассы различных фракций деревьев и моделирования их динамики.
Статистические модели также разрабатывались для изучения закономерностей накопления углерода всеми компонентами древостоев в процессе роста и для составления таблиц биологической продуктивности насаждений различных классов бонитета [8-10]. Подробные обзоры моделей роста приведены в ряде монографий [4, 11, 12].
Все перечисленные модели имеют преимущественно теоретический характер. В лучшем случае они применяются для описания динамики роста в высоту или по объёму стволов, поэтому разработка моделей, описывающих динамику компонентов фитомассы деревьев и древосто-ев, продолжает оставаться актуальной. Одному из подходов к решению этой задачи и посвящена данная работа.
Целью работы являлась разработка модели, отражающей фактический ход роста деревьев сосны по массе и объёму, а также динамику фитомассы хвои (листьев), ветвей, тонких и толстых корней.
Для достижения поставленной цели необходимо было создать математическую модель роста дерева, описывающую основные процессы роста отдельных структурных элементов фитомассы деревьев.
Математическое моделирование.
Разработанная модель предназначена для имитации роста деревьев сосны обыкновенной. Основу модели составляет система уравнений, часто использующаяся в теории катастроф для описания систем
«хищник-жертва» [13]. Однако эти уравнения можно использовать и для систем «ресурс-потребитель», что вполне приемлемо для моделирования роста растительных организмов. В этом случае ресурсом выступает продукция фотосинтеза, а потребителем - прирост биомассы растений с учётом дыхания. Общий вид использованных уравнений следующий:
=а( х) х-V (х) у
Cdy=kxV ( х) у - ту,
(1)
где а(х) - можно интерпретировать как скорость продуцирования ассимилятов при фотосинтезе единицы массы хвои (листьев); V(x) - скорость потребления ассимилятов на рост единицей поверхности флоэмы ствола, ветвей и толстых корней; kl - кпд переработки биомассы асси-милятов в биомассу ствола, корней и ветвей; т - коэффициент естественной гибели клеток тканей растения, определяемый продолжительностью жизни клеток заболони древесины; dx/dt - скорость изменения величины фонда углеводов в растении; ёуШ - скорость роста биомассы.
Обычно а(х) задают наиболее простым логистическим уравнением (рост с самоограничением)
а
(х)=а0(1-х / xk X
(2)
где щ - максимальная скорость продуцирования; х] - максимально возможная масса ассимилирующего аппарата. Однако проведённые численные эксперименты показали, что в молодом возрасте при быстром нарастании массы хвои соотношение (2) не обеспечивает стабильности решения, поэтому в качестве функции, саморегулирующей продукцию фотосинтеза (и как следствие рост хвои), было использовано уравнение, применяемое для описания эффектов взаимозатенения в GAP-моделях для светолюбивых видов [12]:
Кр =ехр[-1,136(Мь /Мтах)-0,08], (3)
где Кр - относительная (в долях единицы от максимальной) интенсивность фотосинтеза единицы массы хвои с учётом эффектов взаимозатенения; Мь - фактическая масса хвои, кг; Мтах - максимально возможная масса хвои для данного вида в заданных условиях произрастания, кг. Такая замена вполне правомерна, так как в данной системе важно только то, чтобы эти функции были качественно похожими [13].
Скорость роста, пропорциональную скорости потребления ассимилятов V(x), в моделях растений часто задают S-образной функцией. С учётом того, что скорость роста зависит от концентрации углеводов в тканях, её можно записать следующим образом:
Дтах^Г
V (C) =
(4)
(кп + [С]п)
где [С] - концентрация углеводов в растущих тканях; к и п - параметры; Атах -максимальная скорость прироста биомассы, приходящаяся на 1 дм2 поверхности флоэмы ствола, ветвей и толстых корней
[14].
Концентрацию углеводов в растущих тканях можно принять равной их концентрации во флоэме, т. е. обратно пропорциональной площади поверхности
[С]. = С, L h S
(5)
где С - общее количество углеводов в растении, кг; £ - площадь поверхности, дм ; /=1...3 - ветви, ствол и толстые корни соответственно. Тогда, подставляя соотношение (5) в (4) и умножая его на соответствующие площади, после преобразования получим уравнения для расчёта годичного прироста ветвей (dMcу/di), ствола (dMs/di) и толстых корней (dMA/di) соответственно:
dMc dt
А CnS
Amax C1 Scy ,
knSny + C1n'
(6)
dMst = Amax C^st ,
dt knSnst + C? '
dMA_ АтаХС^А ,
(7)
(8)
dt кп^пА + с
где С1 - количество углеводов в растении, образованное за вегетационный период, за исключением расходов на дыхание и формирование тонких корней в весенне-летний период, кг. Общая продукция фотосинтеза в вегетационном периоде и количество углеводов, которое может быть израсходовано на рост ствола, ветвей и толстых корней (С1), рассчитывается по формулам
с = крккрмь + С^, (9)
с !=(с -яр (му +мк) -
- (1+яг ^мг / dt),
где С - общее количество углеводов, образованное в течение вегетационного периода, кг; С^ - остаток углеводов из прошлого вегетационного периода, расходуемый на рост в текущем вегетационном периоде, кг; dMr/dt - весенний прирост тонких корней, кг; Яр - коэффициент дыхания поддержания; Я, - коэффициент дыхания роста; Му - общая живая биомасса ветвей, ствола и толстых корней, кг; МЯ - масса тонких корней, кг; крн - максимально возможная величина годичной продукции фотосинтеза единицы массы хвои, зависящая от условий произрастания (класса бонитета) и густоты древостоя, определяющей интенсивность конкуренции за свет, кг кг-1 г.-1
Интегрирование системы (1) проводилось в следующей последовательности. Сначала по формуле (9) определялась общая продукция фотосинтеза в предстоящем вегетационном периоде. Затем, учитывая приоритеты расхода растения на дыхание живых тканей и рост тонких корней, по уравнению (10) определялось начальное количество углеводов (С1), которое расходуется на рост ствола, ветвей и толстых корней, рассчитываемых по уравнениям (6)-(8). Расчёты ведутся по
годам. При этом год разбивается на пять шагов интегрирования. В процессе расчёта на каждом шаге величина С1 корректируется с учётом снижения количества углеводов вследствие роста биомассы дерева ^С). Такая корректировка осуществляется по формулам:
dC=~R„ (Mv + Mr ) - (1+Rr )
dM„
dt
-(1+Rr )
/ ddMsL++ dMA dt dt dt
(11)
С =С + dC. (12)
Общая масса тонких корней зависит от влажности почвы, обеспеченности её минеральными элементами (от условий произрастания), климата, интенсивности корневой конкуренции и др. Для сходных местообитаний можно предположить, что масса корней пропорциональна массе хвои:
Mr = krML,
dM
Ro _
dt
= kjkRM f 2 ,
(14)
где dMRc> - осенний прирост корней за счёт фонда запасных питательных веществ, кг; kl - коэффициент эффективности переработки углеводов в биомассу; Мр - часть фонда запасных питательных веществ, расходуемая на рост, кг; ^ - доля фонда запасных углеводов, идущая на рост тонких корней. Принято, что за счёт фонда запасных веществ осенью формируется 50 % всей массы тонких корней. Остальная часть образуется весной-летом следующего года в течение вегетационного периода. В этом случае массу корней, которую растению необходимо сформировать весной - летом текущего вегетационного периода (dMr/dt), можно определить по формуле:
dMr = M - dMRo dt r
dt
(15)
(13)
где &г - коэффициент, зависящий от региона, лесорастительных условий и др.
Это предположение можно обосновать и тем, что при недостатке влаги в почве и недостаточном количестве корней рост всех частей растения замедляется или прекращается, а корни начинают интенсивно расти, пока не будет достигнут баланс между системой «корни - листья», обеспечивающий устойчивое функционирование растения в создавшихся условиях [15]. В данной модели коэффициент &г принят равным единице.
Продолжительность жизни тонких корней довольно изменчива и составляет от нескольких недель до одного года, а иногда и более. В модели она принята равной одному году. Рост тонких корней проходит в два этапа. Первый начинается осенью в конце вегетационного периода и осуществляется за счёт фонда запасных питательных веществ и оттока углеводов из отмирающей хвои
Изменение величины фонда запасных углеводов определяется как остаток неизрасходованных на рост и дыхание углеводов к концу вегетационного периода и оттока части углеводов из отмирающих органов растения (хвоя и тонкие корни)
dMf = С, + к М , (16)
f 1 от оп' \ /
где &от - коэффициент оттока подвижных углеводов из листьев, принятый равным 0,3 [15]; Моп - масса опада хвои, кг; С1 - остаток неизрасходованных углеводов к концу вегетационного периода, кг; dMf - увеличение фонда запасных углеводов к концу текущего вегетационного периода, кг.
Общая фитомасса дерева к концу вегетационного периода составит
МА = М А_1 +£ М, (17) 1=1 dt
где /=1...5 обозначает прирост ветвей, ствола, толстых, тонких корней и генеративных органов; А - возраст, лет.
В осенний период после прекращения роста идёт накопление запасных питательных веществ, из которых весной формируются хвоя и часть раннего (весеннего) прироста. При этом углеводы откладываются в клетках в виде зёрен крахмала и других сложных неподвижных форм угле-
водов. При избытке питательных веществ предположено, что они идут на осенний рост корней. Используя это предположение, потенциальную «ёмкость» тканей для запасных питательных веществ можно определить из соотношения:
Fp = кг(Му -М,)-Мл, (18)
где Мр - величина фонда запасных веществ в живых тканях на середину вегетационного периода, кг; кр - коэффициент «ёмкости» живой биомассы для запасных питательных веществ; М¥ - масса живых тканей к концу вегетационного периода с учётом прироста, кг.
Общее количество запасных веществ, которое может отложиться в клетках тканей их запасающих, составляет:
F1 = kpdMf, (19)
где кр - коэффициент трансформации из подвижных углеводов в запасаемые, а dMf определяется по (16).
Если потенциальная «ёмкость» тканей больше общего количества запасаемых углеводов (^), то масса фонда к концу вегетационного периода составит:
М = 1 + F2. (20)
В случае, когда количество запасаемых углеводов больше потенциальной «ёмкости», фонд запасных питательных веществ полностью заполняется до его максимальной величины, а их избыток добавляется к остатку углеводов (С^Д переходящих на следующий вегетационный период:
Мр= Мл + Fр, (21)
С1аМ = кр ( dMf - Fp ) , (22)
где кр - коэффициент трансформации подвижных углеводов в тканях растений.
Динамика фонда запасных углеводов с конца вегетационного периода текущего года к началу следующего описывается следующим образом. Часть запасных веществ остаётся в виде резерва, а часть переходит в подвижные углеводы и расходуется осенью текущего года на рост тон-
ких корней и весной следующего года на рост хвои
dF2 = кгМ/, (23)
= Мг - dF2, (24)
Мр 2 = kpdF2, (25)
где dF2 - расходуемая часть запасных веществ, кг; Мр - масса углеводов с учётом затрат на их трансформацию из запасных в подвижные, которая расходуется на рост тонких корней осенью, рассчитываемый по (14), и новой хвои весной следующего года, кг; кр - доля от общего фонда запасных веществ, трансформирующаяся в подвижные углеводы. Коэффициент кр введён ввиду того, что часть запасного фонда остаётся в «законсервированном» виде как резерв для компенсации потерь и повреждений живых тканей дерева в результате воздействия различных факторов, а также расходуемый при плодоношении.
Из работ по физиологии древесных растений известно, что в весенний период формирование листьев до достижения ими примерно 50-60 % конечной массы идёт за счёт запасных питательных веществ [15], поэтому в модели масса новой хвои следующего года полностью определяется величиной фонда запасных углеводов предыдущего года
Ми = кркьМр 2/ к 2, (26)
где кь - доля подвижных углеводов, расходуемая на рост новой хвои; Мь1 - масса хвои будущего года, формируемой за счёт запасного фонда, кг; к2 - доля роста хвои за счёт запасных веществ.
Для сосны в зоне хвойно-широколиственных лесов средняя продолжительность жизни хвои составляет три года, поэтому масса хвои в следующем году составит:
Мь(А +1) = ^МЩ1) -Моп + Мы, (27)
г=1
Моп = Мп(А +1 - п), (28)
где А - возраст дерева, лет; п - максимальный возраст хвои, лет; Моп - масса хвои в возрасте трёх лет, кг. Таким обра-
зом, к началу следующего года хвоя в возрасте трёх лет опадает, а общая масса увеличивается за счёт образования новой хвои. Для сосны на территории России п может изменяться от одного-двух лет в загрязнённых районах до шести лет в зоне лесостепи.
Остаток подвижных запасных углеводов, переходящий на рост дерева в следующий вегетационный период, определяется из соотношения
СЬа* = (1 - К - К )Mf 2 + СЬаМ. (29) Полученный остаток прибавляется к общей продукции фотосинтеза следующего года, определяемого по (6), и весь цикл расчётов повторяется для следующего вегетационного периода.
В качестве факторов, регулирующих количество живой биомассы в растении на начало следующего вегетационного периода, выступают продолжительность жизни клеток в древесине ствола, ветвей, толстых корней и обеспеченность их запасными углеводами в осенне-зимний период. Баланс живых тканей рассчитывается по следующей схеме. Сначала определяется общий вес живой биомассы, где может разместиться накопленный к концу вегетационного периода фонд запасных углеводов (МГ11у)
М,
MFllv =
kt
(30)
где kf - коэффициент «ёмкости» живых тканей для запасаемых углеводов. После этого рассчитывается ежегодно отмирающая биомасса (Mmd) на основании предположения, что с заданной вероятностью погибнут живые ткани, которым не хватит запасных углеводов на зимний период
Мтй =(М11У - М¥1п )ку, (31)
где Мцу - масса живых тканей в древесине ствола, ветвей и толстых корней, кг; ку -вероятность гибели тканей, оставшихся без запасных углеводов.
Если за период, равный продолжительности жизни клеток (4), общая погибшая масса, рассчитанная по (31), меньше прироста древесины ^ лет назад,
то масса погибших клеток в текущем году принимается равной приросту живых тканей ^ лет назад
, 1=1
^^ death
dt
(32)
'(A-tk)
где dM1 / dt (1=1,3 ) - приросты ствола, ветвей и толстых корней ^ лет назад, кг.
Если же за этот период общая погибшая масса больше прироста древесины ^ лет назад (т. е. гибель из-за недостатка запасных углеводов превышала гибель в результате старения), масса погибших тканей в осенне-зимний период принималась равной массе, рассчитанной по (31).
Таким образом, количество живой биомассы в растении на конец вегетационного периода текущего года или на начало следующего (Му(А+1)) рассчитывается по балансовому уравнению
Mliv = Mliv - M death
Л dM,
+ > -
La
i =1
dt
(33)
где Y^dM1/dt - сумма приростов сучьев, ствола и толстых корней в текущем году, кг.
Масса ствола (М^) и его объём (К^) к концу вегетационного периода составят:
М (34)
Mst = Mst +■
dt
Vst= Mst / p,
(35)
где р - плотность древесины, принятая для сосны в расчётах равной 0,52 кг м .
Масса сучьев определяется из соотношения:
dM
My = My + cy cy dt
cy
- M
opad ?
(36)
где Mopad - масса отмерших сучьев за вегетационный период, принятая пропорциональной общей их массе
МоРа* = КраМсу , (37)
где корас1 - скорость отмирания ветвей, зависящая от густоты древостоя, его потенциальной производительности (класса бонитета) и генетических особенностей популяции, кг-г"1. Для толстых корней опад не рассчитывался ввиду большой неопределённости его величины.
Таким образом, общая величина ежегодного опада определялась как сумма биомассы хвои в возрасте трёх лет, общей массы годового прироста тонких корней и опада ветвей.
Диаметр ствола и площадь его поверхности рассчитываются по эмпирическим формулам, коэффициенты которых установлены на основании обмеров модельных деревьев
d6K = exp(0,5255ln Vst --0,5354ln h+1,5431), dK = exp(-0,3573 ln A-0,07749 ln h+ +0,4045 ln Vst -0,2622)+dKK, Sa = exp(0,5596 ln Vst +
(38)
(39)
(40)
+0,36121п h+2,4026), где dвк и dбк - диаметры в коре и без коры на высоте 1,3 м соответственно, см; А -возраст дерева, лет; И - его высота, м; Vst -объём ствола без коры, дм ; Sst - площадь поверхности ствола, дм2.
Площадь поверхности ветвей ^су) рассчитывалась по эмпирической формуле
77М
=---. (41)
^ 0,027А + 0,33
Контрольные расчёты по модели показали, что использование уравнения (6) для расчётов приростов массы ветвей с площадью поверхности, вычисляемой по (41), даёт завышенные значения, поэтому для концентраций был введён поправочный коэффициент ксу, изменяющийся в пределах 0,25...0,35, и уравнение (6) было заменено уравнением
dMcУ_Amax(kcУC1) nSy
dt knSn + (kC)n
dH dt
= к
cb
где кИ - коэффициент скорости роста, соответствующий максимально возможному годичному приросту в высоту, м г-1; к2 -концентрация, при которой скорость роста достигает 1/2 её максимальной величины, кг дм-3; Ь - эмпирический коэффициент.
Концентрация в начале вегетационного периода зависит от величины остатка фонда углеводов, сформированного в предыдущем периоде (С^), и объёма ствола без коры
С
(44)
C2 =
' last
V„
(42)
'су ' \ксус\,
Для описания прироста в высоту в имитационной модели используется уравнение, в котором скорость роста зависит от концентрации углеводов в стволе дерева в начале вегетационного периода (С2)
кь + Cb
(43)
где кг - количество углеводов, расходуемых на весенний рост корней, кг. Расчёт прироста проводится одновременно с вычислением динамики концентрации и биомассы дерева в основном блоке модели.
Техника эксперимента и методы обработки результатов. Проверка модели роста показала, что для растений в возрасте до 4-10 лет формула (43) даёт завышенные результаты. Связано это с тем, что в этом возрасте объём ствола в зависимости от условий роста составляет всего 0,1.0,7 дм-3, поэтому удельный объём ветвей в общем объёме биомассы растения достаточно высок и игнорирование их массы при расчётах концентрации углеводов на единицу объёма может существенно завышать оценки концентраций, а значит и приростов в высоту. В связи с этим введён ограничитель предельного прироста в высоту в этом возрасте на основе предположения, что в молодых растениях с объёмом ствола до 0,7 дм-3 на прирост в высоту расходуется около 40 % общего прироста ствола (<М^Ш). Тогда объём прироста верхушечного побега можно определить делением этой части биомассы на плотность его древесины. Разделив объём прироста на площадь сечения побега, можно установить его длину. В этом случае в молодом возрасте дополнительно рассчитывается максимально возможный прирост в высоту по формуле
dH 0AdM„t / dt
dt
(45)
р1кЯ2 '
где р1 - плотность древесины верхушечного побега, принятая в модели равной 0,4 кгдм" ; R - его радиус, дм. Если прирост в высоту, вычисленный по (43), превышает максимальный по (45), то его величина принимается равной dHm/dt. А если прирост dH/dt по (43) меньше максимально возможного, то он принимается за истинный.
Таким образом, рост дерева зависит от общего количества углеводов, синтезируемых растением в течение вегетационного периода с учётом остатка углеводов, расходуемых из запасного фонда, созданного в прошедшие годы. В свою очередь эта величина определяется потенциальным плодородием почв, густотой древостоя и положением дерева в пологе, от которого зависит степень затенения кроны, а значит и интенсивность фотосинтеза. Кроме того, использование такой схемы позволяет учитывать при моделировании роста влияние на величину валовой продукции фотосинтеза количество осадков и температуры, т. е. учитывать и влияние климата.
Для расчётов по модели, помимо общих для всех деревьев одной породы постоянных коэффициентов, необходимо для каждого растения установить всего шесть параметров. Это коэффициент скорости роста в высоту кц и коэффициент к2, представляющий собой концентрацию, при которой скорость роста достигает 1/2 его максимальной величины из уравнения (43), концентрация углеводов во флоэме (к), при которой скорость прироста древесины ствола составляет 1/2 её максимальной величины (4), интенсивность опада сучьев и ветвей (корас) в формуле (37), радиус верхушечного побега молодых деревьев (К) объёмом до 0,7 дм-3, поправочный коэффициент ксу для корректировки концентрации углеводов во флоэме ветвей при расчётах прироста их массы (42).
Основным параметром, определяющим скорость роста дерева, является ко-
эффициент крЦ. По своей сути он представляет годичный нетто-фотосинтез единицы массы хвои данного дерева, зависящий от интенсивности солнечной радиации, доходящей до дерева с учётом её поглощения в пологе древостоя, доступности влаги и обеспеченности элементами питания. Таким образом, при имитации роста дерева в древостое достаточно задать шесть индивидуальных параметров, перечисленных выше, а основным управляющим параметром будет коэффициент фотосинтеза, изменяющийся в онтогенезе дерева и зависящий от динамики густоты древостоя вследствие самоизреживания и проведения различных видов рубок. Средняя величина его для древостоя зависит от типа леса, почвы и условий увлажнения. Другими словами, рассматриваемая модель позволяет учитывать различия в типах леса и потенциальной производительности древостоев (классах бонитета).
При расчётах прироста ветвей по (42) используется поправочный коэффициент ксу к концентрации углеводов во флоэме ветвей. Его величина довольно тесно связана с параметрами роста массы ствола
ку =-1,2503к+16,1978к2 -* (46)
-0,6005кср + 0,2997,
ср
где кср - среднее значение коэффициента к для данного древостоя, которое рассчитывается по формуле
кср = 0,2117ехр[0,3077(Ьои + 2)]. (47)
При расчёте динамики массы ветвей по (36) - (37) используется коэффициент, характеризующий скорость их отмирания (кораД который устанавливается по уравнению
Км =1,3604ехр(-19,716к - 2,106) -
(48)
-0,0041.
Для расчёта прироста в высоту в молодом возрасте до 5-10 лет по формулам (43) - (45) используется радиус верхушечного побега (К), который зависит от скорости роста в высоту, интенсивности солнечной радиации и густоты молодняка,
т. е. напряжённости конкурентных отношений. Учитывая, что такие данные часто отсутствуют, его величина определяется по формуле
R=0,0341кИ +0,0462р-0,0068 (49)
где р - полнота, устанавливаемая по таксационному описанию прежнего лесоустройства или других архивных источников. При отсутствии таких данных можно использовать современную полноту древостоя, хотя такое решение будет приблизительным.
Таким образом, после определения шести индивидуальных коэффициентов для каждого отдельного дерева выполняется цикл расчётов биомассы и размеров (диаметра и высоты) ствола, биомассы хвои, ветвей, толстых и тонких корней в возрасте от двух до 120-150 лет, что вполне достаточно для практических целей.
Использование рассмотренной модели роста деревьев сосны возможно в нескольких вариантах. Первый вариант предполагает использование данных перечёта на пробной площади и рубку моделей с полным анализом ствола, определением компонентов надземной фитомассы и годичной продукции деревьев. В этом случае для каждого такого дерева рассчитываются шесть индивидуальных коэффициентов (кь, к, к2, ксу, кора^ R) по приведённым выше формулам.
После этого выполняются расчёты по модели. Полученные результаты сравниваются с фактическими данными. Эта процедура используется для проверки адекватности модели, её настройки, а также при уточнении значений используемых параметров и формул для расчётов индивидуальных коэффициентов.
При втором варианте модель используется в качестве имитационной системы для расчётов таксационных и биометрических характеристик отдельных деревьев и древостоев. В этом случае по формулам (46)-(49) рассчитываются параметры для 30 или более деревьев в зависимости от
моделируемой площади древостоя и его полноты. К полученным значениям коэффициентов добавляют нормально распределённые случайные величины, умноженные на среднеквадратические отклонения (±8о). Аналогично рассчитываются коэффициенты интенсивности фотосинтеза. Нормально распределённые случайные величины получают с использованием метода Монте-Карло, встроенных датчиков случайных чисел или иными способами. После обработки результатов расчётов статистическими методами получают средние таксационные характеристики моделируемого древостоя, его надземную и общую фитомассу, массу опада и другие характеристики.
Результаты. Для проверки адекватности модели были использованы данные измерений стволов 48 деревьев сосны, срубленных на семи пробных площадях. Пробные площади были заложены в Богородском, Сосновском лесничествах Нижегородской области и в Костромском лесничестве Костромской области. Возраст древостоев составлял 45-120 лет, класс бонитета насаждений изменялся от Ia до III. У каждого срубленного модельного дерева измерялись объём ствола, текущий прирост объёма по десятилетиям, масса хвои и ветвей. Масса стволов рассчитывалась через их объём и удельную плотность абсолютно сухой древесины.
С использованием разработанной модели для каждого дерева вычислялись объём ствола, масса хвои, ветвей и другие показатели по приведённым выше формулам. Расчёты проводились с возраста один год до возраста, в котором эти модели были срублены. Полученные результаты для объёмов стволов достаточно близки к фактическим, а коэффициент детерминации модели составляет 0,849.
Для массы хвои разброс фактических и расчётных значений несколько выше. Объясняется это значительным варьированием фитомассы ассимиляционного аппарата даже у деревьев одного возраста и
размера. Связано это с тем, что масса хвои и ветвей существенно зависит от истории роста и развития древостоя, его густоты, интенсивности изреживаний в прошлые годы. Тем не менее, результаты расчётов имеют хорошую точность (за исключением двух модельных деревьев), а коэффициент детерминации составляет 0,714.
Для ветвей точность модели аналогична точности расчётов для фитомассы хвои. Коэффициент детерминации здесь равен 0,72.
В процессе расчётов для древостоев каждой пробной площади вычислялась средняя высота как среднее арифметическое по моделируемым деревьям. Отклонения от фактических средних высот колебались от +3,6 до -9,0 %, что находится в пределах ошибок измерений при инвентаризации насаждений.
Таким образом, модель вполне адекватно описывает динамику роста деревьев сосны в разных условиях произрастания, начиная с возраста два года и заканчивая возрастом 145 лет.
Выводы. Предлагаемая модель при наличии данных измерений модельных деревьев позволяет реконструировать их ход роста с получением расчётных данных о фитомассе хвои и ветвей в различных возрастных периодах. Кроме того, имеется возможность прогнозирования роста отдельных деревьев и древостоя на ближайшие десять лет. Если дополнительно привлечь модель, имитирующую динамику густоты насаждения, то продолжительность такого прогноза можно увеличить до возраста 100-150 лет, т. е. до возраста спелости сосновых древостоев.
Достоинством рассмотренной модели является то, что она имеет не только тео-
ретический характер, но и описывает фактический рост деревьев по объёму ствола, динамику фитомассы хвои, ветвей с момента появления всходов и до предельного возраста их жизни. При этом учитывается разнообразие условий местопроизрастания, индивидуальные (наследственные) особенности деревьев, а также внутриценотические отношения. Точность расчётов для отдельных деревьев будет значительно выше, если объединить модель отдельных деревьев в модель древостоя. В такой системе взаимодействие и взаимовлияние деревьев друг на друга в процессе роста и развития определяются также густотой насаждений в разных возрастах. Кроме того, ход роста зависит от интенсивности и повторяемости различных видов рубок, выполняемых в лесу.
При отсутствии данных измерений модельных деревьев можно, задавая тип леса и соответствующий ему класс бонитета, начальную густоту, возраст смыкания крон (Ать), распределение количества деревьев по относительным высотам и диаметрам, рассчитывать биометрические показатели деревьев в различных возрастах. Однако для получения информации о параметрах всего древостоя необходимо дополнительно использовать модели са-моизреживания, аналогичные применяемым в GAP-моделях или других моделях динамики роста и густоты древостоев.
Рассмотренная модель достаточно универсальна и при наличии данных полевых измерений может быть легко адаптирована для описания роста и динамики фитомассы отдельных компонентов деревьев как хвойных, так и лиственных пород.
Список литературы
1. Чертов О. Г. Математическая модель экосистемы одного растения // Журнал общей биологии. 1983. Т. 44, № 3. С. 406-414.
2. Бугровский В. В., Дудин Е. Б., Мелли-на Е. Г., Цельникер Ю. Л. Моделирование продукционных процессов в чистых древостоях // Журнал общей биологии. 1982. Т. 43, № 4. С. 480-488.
3. Бугровский В. В., Лютов Л. И., Мелли-на Е. Г. Моделирование развития древостоев с учетом климатических условий // Журнал общей биологии. 1987. Т. 48, № 1. С. 84-91.
4. Кулль К., Кулль О. Динамическое моделирование роста деревьев. Таллинн: Валгус, 1989. 232 с.
5. Старцев А. И., Корепанов А.А. Анализ причин усыхания пойменных дубрав р. Урал с использованием динамической модели роста // Дуб -порода третьего тысячелетия: сб. науч. тр. -Гомель: ИЛ НАН Беларуси, 1998. Вып. 48. С. 384-391.
6. Демаков Ю. П., Смыков А.Е., Гаврицко-ва Н.Н. Структура, продуктивность и динамика осинников Республики Марий Эл // Вестник Марийского государственного технического университета. Сер.: Лес. Экология. Природопользование. 2011. № 2 (12). С. 24-38.
7. Демаков Ю. П., Пуряев А.С., Черных В.Л., Черных Л.В. Использование аллометрических зависимостей для оценки фитомассы различных фракций деревьев и моделирования их динамики // Вестник Поволжского государственного технологического университета. Сер.: Лес. Экология. Природопользование. 2015. № 2 (26). С. 19-36.
8. Усольцев В. А. Фитомасса лесов Северной Евразии: нормативы и элементы географии. Екатеринбург: УрО РАН, 2002. 762 с.
9. Курбанов Э.А. Бюджет углерода сосновых экосистем Волго-Вятского района. Йошкар-Ола: МарГТУ, 2002. 298 c
10. Курбанов Э.А. Углерододепонирующие насаждения Киотского протокола. Йошкар-Ола: МарГТУ, 2007. 184 с.
11. Chertov O. G., Komarov A. S., Karev G. P. Modern Approaches in Forest Ecosystem Modelling: EFI Research Report 8. Brill// Leiden, Boston, Köln., 1999. 130 p.
12. Shugart H. H. A theory of forest dynamics. New York: Springer-Verlag, 1984. 278 p.
13. Свирежев Ю. М. Нелинейные волны, дис-сипативные структуры и катастрофы в экологии. М.: Наука. Гл. ред. физ.-мат. лит., 1987. 369 с.
14. Сиротенко О. Д., Варчева С. Е. Прикладная авторегуляционная модель роста растений // Математическое моделирование в агрометеорологии. Тр. ВНИИСМ. Вып. 26 Л.: Гидрометеоиздат, 1990. С. 45-52.
15. Крамер П. Д., Козловски Т. Т. Физиология лесных растений. М.: Лесная промышленность, 1983. 464 с.
Статья поступила в редакцию 24.10.15.
Информация об авторе
СТАРЦЕВ Алексей Иванович - кандидат сельскохозяйственных наук, доцент кафедры лесной таксации и лесоустройства, Нижегородская государственная сельскохозяйственная академия. Область научных интересов - фитомасса деревьев и древостоев, математическое моделирование роста деревьев и древостоев, информационные технологии в лесном хозяйстве. Автор 50 публикаций.
UDC 581.526
SIMULATION MODEL FOR DYNAMICS ESTIMATION OF SCOTS PINE BIOMASS
А. I. Startsev Nizhny Novgorod State Agricultural Academy, 97, Gagarin St., Nizhny Novgorod, 603107, Russian Federation E-mail: astartsev2005@rambler.ru
Key words: tree growth model; stem biomass; needles; branches; roots; gross annual photosynthesis; carbohydrate reserve.
ABSTRACT
The paper is aimed at the development of a simulation model for estimation of Scots Pine growth. The model features nine phase variables: stem mass (volume), biomass of branches, needles, tender and thick roots, bark, gross annual photosynthesis, carbohydrate reserve and increment. The growth rate of a woody plant depends on the total amount of carbohydrates synthesized by the tree throughout the vegetation period. This value is based on the level of soil fertility and humidity, density of the wood stand, tree position in the canopy which defines the level of crown shading. Individual characteristics of the growth rate arising as a result of hereditary properties are taken into consideration. These include height growth rate, biomass increment rate, etc. Conclusion. The model developed allows calculating biometric indicators and biomass of certain tree components aged from 2 to 120-150 years. Model verification was carried out using the data from 48 cut down pine trees aged from 37 to 145 years. The determination coefficient (R2) amounted 0.849 for stems, 0.714 for needles and 0.72 for branches. The model is well adaptable for other tree species.
REFERENCES
1. Chertov O. G. Matematicheskaya model ekosistemy odnogo rasteniya [Mathematical model of a plant ecosystem]. Zhurnal obshchej biologii [General biology journal]. 1983. Vol. 44, No 3. Pp. 406-414.
2. Bugrovskiy V. V., Dudin E. B., Mellina E. G., Celniker Yu. L. Modelirovanie produkcionnyh pro-cessov v chistyh drevostoyah [Production processes modeling in pure sands]. Zhurnal obshchej biologii [General biology joural]. 1982. Vol. 43, No 4. Pp. 480-488.
3. Bugrovskij V. V., Lyutov L. I., Mellina E. G. Modelirovanie razvitiya drevostoev s uchetom klimat-icheskih usloviy [Modelling of stand development based on the climatic region]. Zhurnal obshchej biologii [General biology journal]. 1987. Vol. 48. No 1. Pp. 84-91.
4. Kull K., Kull O. Dinamicheskoe modelirovanie rosta derev'ev [Dynamic modeling of tree growth]. Tallinn: Valgus, 1989. 232 p.
5. Startsev A. I., Korepanov A.A. Analiz prichin usyhaniya pojmennyh dubrav r. Ural s ispolzovaniem dinamicheskoj modeli rosta. Dub-poroda tret'ego ty-syacheletiya: sb. nauch. tr. - Gomel [Analysis of reasons of floodplain forest drying out by River Ural using the dynamic growth model]: IL NAN Belarus, 1998. Iss. 48. Pp. 384-391.
6. Demakov Yu. P., Smykov A.E., Gavritsko-va N.N. Struktura, produktivnost i dinamika osinnikov Respubliki Mari El [Structure, productivity and dynamics of aspen forest in Mari El]. Vestnik Marijskogo gosudarstvennogo tekhnicheskogo univer-siteta. Seriya: Les. Ekologia. Prirodopolzovanie [Vestnik of Mari State Technical University. Series: Forest, Ecology, Nature Management]. 2011. No 2 (12). Pp. 24-38.
7. Demakov Yu. P., Puryaev A.S., Chernyh V.L., Chernyh L.V. Ispolzovanie allometricheskih zavisi-mostej dlya otsenki fitomassy razlichnyh fraktsiy dereviev i modelirovaniya ih dinamiki [Application of
allometric dependencies for the purpose of biomass estimation of different tree fractions and modeling their dynamics]. Vestnik Povolzhskogo gosudarstven-nogo tekhnologicheskogo universiteta. Seriya: Les. Ekologiya. Prirodopolzovanie [Vestnik of Volga State University of Technology. Series: Forest Ecology, Nature Management]. 2015. No 2 (26). Pp. 19-36.
8. Usoltsev V. A. Fitomassa lesov Severnoj Ev-razii: normativy i ehlementy geografii [Biomass of forests in North Eurasia: regulatory standard and geography elements]. Yekaterinburg: UrO RAN, 2002. 762 p.
9. Kurbanov E.A. Byudzhet ugleroda sosnovyh ehkosistem Volgo-Vyatskogo raiona [Carbon budget of pine ecosystem in Volga-Vyatka regions]. Yoshkar-Ola: MarGTU, 2002. 298 p.
10. Kurbanov E.A. Uglerododeponiruyushchie nasazhdeniya Kiotskogo protokola [Carbon sequestration forest plantation of the Kyoto protocol]. Yoshkar-Ola: MarGTU. 2007. 184 p.
11. Chertov O. G., Komarov A. S., Karev G. P. Modern Approaches in Forest Ecosystem Modelling: EFI Research Report 8. Brill// Leiden, Boston, Köln., 1999. 130 p.
12. Shugart H. H. A theory of forest dynamics. New York: Springer-Verlag, 1984. 278 p.
13. Svirezhev Yu. M. Nelinejnye volny, dissipa-tivnye struktury i katastrofy v ehkologii [Nonlinear waves, dissipative constructions and catastrophes in ecology]. Moscow: Nauka. Gl. red. fiz.-mat. lit., 1987. 369 p.
14. Sirotenko O. D., Varcheva S. E. Prikladnaya avtoregulyatsionnaya model rosta rastenij [Applied autoregulating model of plant growth]. Matematich-eskoe modelirovanie v agrometeorologii [Mathematic modeling in agricultural metrology] Tr. VNIISM. Iss. 26. Leningrad: Gidrometeoizdat, 1990. p. 45-52.
15. Kramer P. D., Kozlovski T. T. Fiziologiya le-snyh rastenij [Physiology of forest plants]. Moscow: Lesnaya promyshlennost, 1983. 464 p.
The article was received 24.10.15.
Citation for an article: Startsev A. I. Simulation Model for Dynamics Estimation of Scots Pine Biomass.
Vestnik of Volga State University of Technology. Ser.: Forest. Ecology. Nature Management. 2016. No 1(29).
Pp. 30-41.
Information about the author
STARTSEV Aleksey Ivanovich - Candidate of Agricultural Sciences, Associate Professor of the Department of Forest Inventory and Forest Management at Nizhny Novgorod State Academy of Agriculture. Research interests - tree and stand biomass, mathematical modeling and tree and stand growth, information technology in forestry. Author of 50 publications.
