CC BY
519
УДК 519:577.4+502.7
МОДЕЛЬ ЯРУСНО-МОЗАИЧНОГО ЛЕСА И МОДЕЛИРОВАНИЕ СУКЦЕССИИ
А.К. Гуц, Е.О. Хлызов
Предлагается математическая модель мозаичного многоярусного леса с целью описания серий сукцессии.
Введение
В этой статье ставится задача построения математической динамической модели мозаичного многоярусного леса, с помощью которой можно было бы прогнозировать состояния лесной экосистемы, подверженной влиянию внешних управляющих факторов.
Модель основывается на четырёх управляющих внешних факторах, задающих среду экосистемы. Это влажность почвы w, мозаичность т, наличие конкуренции к и антропогенное вмешательство а в лесную экосистему (вырубка леса, пожары и т.д.). Модель может быть усложнена за счёт введения дополнительных управляющих внешних факторов, но при этом она становится менее наглядной и требует при её использовании уже гораздо более серьёзных математических знаний.
1. Понятие ярусности леса
Ярус - это неоднородность в вертикальном распределении фитомассы леса.
Дюрье выделяет крупные ярусы деревьев, кустарников, трав, напочвенного покрова, а в пределах этих ярусов - подъярусы:
- верхний, первый ярус А, образуют деревья. Его подъярус А. 1 состоит из высоких деревьев; подъярус А.2 - деревья второй величины - рябина обыкновенная, черёмуха обыкновенная, ива козья, дикая яблоня.
- второй ярус В состоит из кустарников, образующих подлесок, - лещина обыкновенная, жимолость лесная, крушина ломкая, бересклет европейский.
- третий ярус С леса состоит из трав. Подъярус С.1 высоких трав - чистец лесной, бор развесистый, борцы; подъярус С.2 низких трав - сныть обыкновенная, осока волосистая, пролестник многолетний и др.
- четвёртый ярус D - мхи, грибы, лишайники.
Copyright © 2011 А.К. Гуц, Е.О. Хлызов
Омский государственный университет им. Ф.М. Достоевского
E-mail: [email protected], [email protected]
Ярусное расположение растений связано с неодинаковой освещённостью. Количество света уменьшается от яруса к ярусу. Много света получают деревья первого яруса и очень мало - мхи и лишайники. В еловом лесу кустарники не растут: ветви елей задерживают очень много света, в таком лесу всегда сумрачно.
«Многоярусные биоценозы, представленные большим количеством видов растений, животных и микроорганизмов, связанных между собой разнообразными пищевыми и пространственными отношениями, называются сложными. Они наиболее устойчивы к неблагоприятным воздействиям. Исчезновение какого-либо вида существенно не отражается на судьбе таких биоценозов. В них происходит лишь незначительная перестройка организации, при которой популяции одного и даже нескольких видов могут заменяться экологически близкими видами, а стабильность сообщества определяется количественной регуляцией численности одних видов другими» [1].
2. Понятие мозаичности леса
Мозаичность - это неоднородность в горизонтальном (по площади) распределении фитомассы леса.
Мозаичность связана с неравномерным распределением деревьев в лесу. На рис. 1 представлены регулярное, случайное и контагиозное (пятнистое, мозаичное) распределения.
Рис. 1. Типы мозаичности леса: а) регулярное распределение, б) случайное распределение, в) контагиозное (пятнистое, мозаичное) распределение
Мозаичность в лесу может быть установлена с помощью, например, следующих двух способов измерения: коэффициента дисперсии и индекса размещения.
3. Шестиярусная модель леса
«Наиболее широко принятой в современной в экологии лесных сообществ является ярусно-мозаичная концепция леса как сложной системы» [2].
Лесной фитоценоз, его состояние, характеризуем с помощью такого понятия, как первичная биологическая продуктивность леса. Степень продуктивности фитоценоза в момент времени I определяем как функцию и х = х(£).
Первичная биологическая продуктивность характеризуется образованием биомассы (первичной продукции) в процессе фотосинтеза зелёными растениями (автотрофами), которые образуют первый трофический уровень экосистемы и служат началом всех цепей питания.
Продуктивность лесного биоценоза в момент времени £ будет определять продуктивность леса в следующий момент времени £ + (Й.
Иначе говоря,
х(Ь + ей) = х{£) + А(уЬ, х)(й, (1)
где А(£,ж)сЙ - величина, описывающая отклонения в продуктивности, произошедшие на отрезке времени ей.
Из (1) имеем то, что называется дифференциальным уравнением
<1х .. ч /г>ч
— = А(1,х). (2)
В основе продуктивного процесса растений лежит фотосинтез. Растения под воздействием солнечной энергии, поглощая листьями из атмосферы углекислый газ и корневой системой из почвы воду, создают органическое вещество. Недостаток влаги (—ги) > 0 является фактором, не способствующим благополучию леса.
В таком случае следует в правую часть дифференциального уравнения (2) добавить член (—(—«;)):
^ = А(Ь, х) - (-т). (3)
В самом простом случае можно принять, что А(Ь,х) = к\х, т.е.
с1х , ...
— = к\х + ь). (4)
ей
Коэффициент к\ можно посчитать постоянным. Но тогда доброкачественность леса будет нарастать как геометрическая прогрессия, и это делает бессмысленной нашу модель. Поэтому начнём её усложнять.
Коэффициент к\ отвечает за «прирост доброкачественности леса». Учтём, что для здорового леса обязательной чертой является наличие ярусности.
Выделим шесть ярусов леса: два яруса деревьев, один ярус кустарников, два яруса трав, один ярус напочвенного покрова (см. § 1).
«Каждый ярус, входящий в состав фитоценоза, оказывает влияние на другие ярусы и в свою очередь подвергается их влиянию. Поэтому фитоценоз необходимо рассматривать как нечто целое, а ярусы фитоценоза - как его структурные части, которые в некоторых случаях могут быть относительно самостоятельными» [1].
Взаимовлияние ярусов должно дать вклад в доброкачественность леса вида
01\Х ■ си^х ■ а^х ■ а^ха^ха^ = ах6,
01 —
(5)
Каждый коэффициент сц > 0 характеризует степень участия г-го яруса во взаимодействии ярусов. Доля фитомассы г-го яруса характеризуется величиной сцх. Взаимовлияние г-го яруса и j-ro яруса - это произведение щх ■ cxjx.
Если a,i —> 0, то мы констатируем отсутствие г-го яруса, означающее меньшую ярусность фитоценоза и, следовательно, его меньшую степень биоразнообразия, меньшую устойчивость к неблагоприятным воздействиям.
Примем, что
к\ = ах6 — р(х), (6)
где а = а\а2аза4 = const - «сила» взаимодействия четырёх ярусов леса, а величина р(х) - это то, что снижает продуктивность леса.
Тогда следует принять, что
р(х) = — к4х, (7)
вырубка леса, пожары конкуренция мозаичность
где квадратичный член к^х2 характеризует конкуренцию растений в лесу (часть фитоцеоза - Р\х влияет на другую его часть - [32х. Взаимовлияние конкурирующих растений есть произведение j3\x ■ [32х = ках2). Мозаичность леса пропорциональна, на наш взгляд, всей фитомассе х леса - отсюда член к4х. Коэффициент к4> 0 характеризует степень влияния горизонтальной структуры (мозаичности) на динамику прироста фитомассы.
Объединяя уравнения (4) - (7), получаем следующую модель лесной экосистемы:
7 1 7 3 7 2
— = ах — к2х — кзх + к4х + w = at
О ( OL g ^3 4 . ^4 з . к2 2
= Й18;'' +—Х +Т;Г +—1 +ГО;Г/ (8)
ИЛИ
Нт п
~di = dxV^X,k,m,a,W^ ^
где
V(х, к, rri, a, w) = —х8 + кх4 + тх3 + ах2 + wx, (10)
к3 к4 к2
к = —т=—, а =--------------------.
4 ’ 3 ’ 2
Сопоставляем лесному фитоценозу потенциальную функцию
V(х, к, т, а, и>) =
= +?тгж^+ чаж^/ + ^гих^
наличие вырубка леса, влажность
конкуренции пожары, почвы
ураган,
разлив нефти
Это частный случай катастрофы типа «звезда» [4]:
V(х,р, д, к,т, а, ио) = —х8 + рх6 + дх5 + кх4 + тх3 + ах2 + тх. (11)
Катастрофа «звезда» описывает, вообще-то говоря, семь равновесных состояний сразу, из которых не более четырёх являются устойчивыми (наличие сразу четырёх локальных минимумов у функции V).
Таким образом, имеем следующую модель шестиярусного лесного фитоценоза:
Нт г)
= -т^У(х,к,т,а,т), (12)
ГД6 а
V(х, к, т, а, т) = —х8 + кх4 + тх3 + ах2 + -хих. (13)
8
Отметим, что продуктивность характеризуется неравенством х > 0, наличие конкуренции неравенством к < 0, действенность оконной динамики неравенством т > 0, вырубка лесов, пожары неравенством а < 0, недостаток влаги неравенством т < 0.
При моделировании лесных экологических катастроф с помощью построенной модели важно дать ответ на следующий вопрос: если, например, при пожаре сгорели, скажем, два верхних яруса леса, то можно ли считать, что допустимо дальнейшее использование 6-ярусной модели (12) - (13)? Можно, особенно когда мы ограничиваемся качественным описанием сукцессии. Действительно, при выгоревших ярусах 1,2 мы имеем ситуацию, когда а\ —> 0 и
а,2 —> 0. Следовательно, уменьшается сила взаимодействия ярусов ко, но мы имеем ту же потенциальную функцию (13), но с другим значением а. Для качественного анализа изменение значения величины а не играет существенной роли.
4. Равновесия биоценоза;
сукцессия как последовательная смена равновесий
Сообщество - это группа организмов различных видов, проживающих на общей территории и взаимодействующих между собой.
Конкретное значение продуктивности лесной экосистемы, находящейся в состоянии равновесия, будем интерпретировать как конкретное растительное сообщество. Смена равновесия, сопровождающаяся изменением значения продуктивности лесной экосистемы, - это смена одних сообществ другими (динамика растительного покрова по Сукачёву).
Равновесия лесного биоценоза, описываемого уравнением (12) - (13), находятся как решения х = х{к,т,а,ги) уравнения
дУ (х, к, т, а, т) = р ^
дх
Лесную экосистему можно вывести из состояния равновесия многими способами: пожаром, наводнением или засухой. После такого нарушения равновесия
новая экосистема сама себя восстанавливает, и этот процесс носит регулярный характер, называется сукцессией и повторяется в самых разных ситуациях.
Сукцессия - это последовательный ряд смены серийных (временно существующих) растительных сообществ на конкретном местообитании после выведения конкретной экосистемы из состояния равновесия [1].
Можно также сказать, что сукцессия, как процесс, представляет собой последовательную (необратимую) смену биоценозов, преемственно возникающих на одной и той же территории в результате влияния природных или антропогенных факторов.
Различают множество форм сукцессии: фитогенная, зоогенная, ландшафтная, антропогенная, пирогенная, катастрофическая и др.
Каждое сообщество, промежуточное в серии, называемое стадией, существует достаточно долго и слабо меняется во времени. Поэтому стадию можно также считать равновесным состоянием.
Равновесие в природе на самом деле зависит от окружающей среды, а среда эта постоянно подвержена изменениям. Пожары, наводнения, колебания количества атмосферных осадков оказывают влияние на среду, в которой произрастает лес. И растения, конечно же, не могут не реагировать на эти изменения. Получается, что экосистема все время пытается либо сохранить равновесие, либо попасть в новое равновесие. Вмешательство человека - всего лишь ещё один способ изменить окружающую среду и таким образом повлиять на направление развития экосистемы.
Смена сообщества происходит под влиянием факторов, которые связаны с экосистемой и являются её характеристиками.
На роль таких факторов могут претендовать такие характеристики фитоценоза, как внутривидовая и межвидовая конкуренция, мозаичность, оконная динамика. Изменения этих факторов отражают в себе микроэволюцию фитосреды, т.е. накопление изменений фитосреды, сначала благоприятных для фитоценоза, а затем неблагоприятных, и, значит, характеризуют наступающую смену сообщества, смену равновесия.
Следовательно, сукцессия может моделироваться в рамках математической теории катастроф, созданная для описания резких изменений, смен, т.е. катастроф равновесных состояний. Катастрофы при изменении одних факторов могут интерпретироваться как изменение климаксного сообщества, а других, при фиксированных первых, как смена сообщества в серии.
При смене равновесия в нашей модели меняется значение продуктивности леса х(Ь). Она нами определена так, что характеризует разные её значения, соответствующие разным лесным сообществам, т.е. разным стадиям сукцессии. Смена стадии - это смена равновесия, а в рамках теории катастроф Тома
- это бифуркация, означающая то, что в физике называется фазовым переходом. Следовательно, смена равновесия, смена стадии может описываться с привлечением аппарата теории фазовых переходов. Такой подход применён в работе А.С. Исаева и др. [3]. Поскольку теория катастроф Тома включает в себя теорию фазовых переходов [5,6], то наш подход является более общим, допускающим самые различные обобщения.
Крупное нарушение в окружающей среде1
Со<?6щеСТ5$ серии ___________I_________
Лишайники
Голая земля
І
+
їоро
*
+
ВОДОРОСЛИ
Травы Кустарники
ІШМ -■гіМВі.
Луговая Кустарниковая растительность заросль
Пионерное
сообщество
Лес
Кллыйкснов
сообщество
Рис. 2. Схема сукцессии [7]
Темпы сукцессий неодинаковы. На участке, лишённом растительного покрова, первые стадии сменяют друг друга через год несколько лет, а далее процесс смен замедляется, и более позднее стадии восстановления или формирования растительного покрова занимают десятилетия, а затем и столетия [1].
«Время смены одного сообщества другим сильно различается. Типичная последовательность сукцессий, приводящих к появлению дубрав или сосновых лесов в средней полосе, занимает около 200 лет; при этом скорости ранних сукцессий гораздо выше, чем скорости поздних.
В сходных условиях развиваются сходные сукцессии. Факторами, определяющими состав климаксного сообщества, могут быть климат, рельеф, дренаж почвы и т.д.
Нередко результатом сукцессии на таких территориях является восстановление исходного биогеоценоза. Применительно к лесным сообществам такие процессы получили название демутация» [7].
«При этом конкретная экосистема возвращается в своё исходное состояние и пребывает в нем до тех пор, пока не изменятся климат, рельеф, гидрологический режим, пока вновь не пройдёт пожар или не случится какая-то другая катастрофа. И вновь начнётся новая сукцессия, которая либо приведёт к восстановлению исходного сообщества, либо нет. В результате сукцессии на конкретном местообитании восстанавливается исходное растительное сообщество, называемое геоботаниками климаксовым, или коренным. Коренное сообщество растений устойчиво и в данных климатических условиях не изменяется» [1].
Если внешние факторы к,т,а,т, изменяясь в 4-мерном пространстве с осями к,т,а,и>, пересекают так называемое бифуркационное множество Ву, то происходит резкая смена равновесия экосистемы, происходит то, что математики называют бифуркацией (катастрофой). Это соответствует экологической катастрофе, постигшей лесной биоценоз. Продолжающиеся изменения факторов - это новые бифуркции, новые переходы к новым равновесиям, которые
рассматриваем уже как серийные (временно существующие) растительные сообщества. Возврат факторов к исходным значениям - это создание условий к полному восстановлению биоценза.
Бифуркационное множество Ву находят, исключая х и решая систему уравнений:
дУ(х,щ, ...,ип) _ д2У(х,щ,...,ип) _
дх ’ дх2
5. Смены климаксных сообществ
Последнее сообщество серии (цепи), стабильное и находящееся в равновесии с окружающей средой, называется климаксным сообществом. Дальнейшее изменение климаксного сообщества возможно только при изменении окружающих условий.
Мы видим, что климаксное сообщество - это равновесное состояние экосистемы, которое соответствует набору конкретных значений внешних факторов (влажность почвы, климат, антропогенные воздействия).
Смена климаксных сообществ происходит при действии экзогенных факторов, таких как пожары, разливы нефти, засуха и т.д.
В рамках нашей модели для описания смены климаксных сообществ мы должны рассмотреть проекции бифуркационного множества на плоскость «пожары (вырубка) а - влажность т» или на плоскость «атмосфера т - пожары (вырубка) а».
5.1. Бифуркации в плоскости
«пожары (вырубка) а - влажность т»
Бифуркации в плоскости пожары (вырубка) а - влажность т при фиксированных р = д = 0, т = О, к = —20 (сильная конкуренция) даны на рис. 3.
Взятые точки: (а, из) = (-7, 0); (20, 0); (60, 0); (30, -11); (30,11); (24,-18); (24,18).
5.2. Бифуркации в плоскости «атмосфера т - пожары (вырубка) а»
Фактор т можно трактовать как состояние атмосферы, т.е. как экзогенный фактор, существенно влияющий на состояние лесной экосистемы.
Бифуркации в плоскости атмосфера т - пожары (вырубка) а дана при фиксированных р = д = 0; = 0, к = —20 на рис. 4.
Взятые точки: (т, а) = (-10, 0); (10, 0); (70, 0); (35, -21); (35, 21).
6. Смены стадий
Смены стадий происходят под воздействием эндогенных факторов, таких как конкуренция и мозаичность.
-34.0
-4.0
73.0
.20.0
-34.0
V
173.0
20.0
-34.0
V
'73.0
-34.0
73.0
-34.0
В рамках нашей модели таковыми факторами являются бифуркации в плоскости «конкуренция к - мозаичность га» дана при фиксированных р = q = 0; w = 0, а = 20 на рис. 5. Взятые точки: (к, га) = (—20, 0); (1, 0); (0,10); (0, —10).
Литература
1. Москалюк Т.А. Курс лекций по биогеоценологии. URL: http: //www.botsad. ru/p_papers .htm (дата обращения: 12.05.2009).
2. Ризниченко Г.Ю. Экология математическая. URL: http: //www. library,
biophys .msu. ru/MathMod/EM. HTML#e9 (дата обращения: 16.05.09)
3. Исаев А.С., Суховольский В.Г., Бузыкин А.И., Овчинникова Т.М. Сукцессионные процессы в лесных сообществах: модели фазовых переходов // Хвойные бореальной зоны. 2008. Т.25, № 1-2. С.9-15.
4. Woodcock A., Poston Т. A geometrical study of the elementary catastrophes // Lectures Notes in Math. № 373. Springer, 1974. 257 p.
5. Постон Т., Стюарт И. Теория катастроф и её приложения. М. : Мир, 1980.
6. Гилмор Р. Прикладная теория катастроф. М. : Мир, 1984.
7. Динамика природных сообществ // Биология. Электронный учебник. URL: http: //www.ebio.ru/eko07.html (дата обращения: 16.05.09).
