CC BY
62
Математические структуры и моделирование 2015. № 3(35). С. 56-60
УДК 004.9:631.4
математическая модель взаимосвязи «растительность-почва» в лесных экосистемах
А.К. Гуц
д.ф.-м.н., профессор, e-mail: guts@omsu.ru Л.А. Володченкова
к.б.н., старший преподаватель, e-mail: Volodchenkova2007@yandex.ru Омский государственный университет им. Ф.М. Достоевского
Аннотация. В статье предложена математическая модель плодородия почвы, основанная на дифференциальном уравнении, допускающем катастрофы сборки при изменении влажности почвы. Дана также модель взаимосвязи продукции фитомассы четырёхъярусного леса и плодородия почвы.
Ключевые слова: модель, почва, растительность, лесная экосистема, теория катастрофа, катастрофа сборки.
Введение
Модели продукции фитомассы многоярусного леса, построенные в [1], существенным образом учитывают влажность почвы, но никак не связаны со структурой почвы, с таким её важнейшим фактором как плодородие.
Цель данной статьи состоит, с одной стороны, в построении математической модели плодородия почвы, которое меняется в зависимости от влажности почвы и типа почвообразующей породы, а с другой стороны, построение модели взаимосвязи «растительность-почва».
В стационарном случае (при неизменяющемся плодородии почвы), получаем алгебраическое уравнение, из которого можно выразить влажность почвы W = W(у,р) как функцию меры плодородия почвы у и типа почвообразующей породы р. Подставляя W = W(у,р) в дифференциальное уравнение для продукции фитомассы лесной экосистемы, мы получаем уравнение, которое определяет динамику продукции фитомассы в зависимости от меры плодородия и типа почвообразующей породы.
Результаты данной статьи были анонсированы в материалах Международной научной конференции «Математическое и компьютерное моделирование» (Омск, 2014) [2,3].
Математические структуры и моделирование. 2015. №3(35)
57
1. Модель плодородия почвы с точки зрения катастрофы «сборка»
Плодородие почвы определяется содержанием в ней необходимых для растений питательных веществ, уровнем устойчивой влажности, наличием в почве воздуха, строением и структурой почвы, определяющими возможность развития корней. Уровень плодородия зависит от процента содержания гумуса — комплекса органических соединений, образовавшихся в результате разложения микроорганизмами остатков растительного и животного происхождения.
Высоким плодородием обладают почвы (суглинок), возникшие в результате разрушения мягких горных пород, низким — почвы (подзолистые почвы), образованные при разрушении твёрдых пород.
Большинство характеристик почвы связано со свойствами почвообразующей породы [4]. Условием плодородия почвы является оптимальная и постоянная влажность почвы. Вода — обязательное условие почвообразования и формирования почвенного плодородия. Если слишком сухо, почва становится плотнее, а жизнь в ней почти замирает — прекращается разложение органики и фиксация азота. При избытке воды всё в почве задыхается и начинается вредное бескислородное сбраживание органических остатков (силосование).
Учитывая сказанное, мы ограничимся при построении модели почвы только двумя её основополагающими характеристиками: типом почвообразующей породы и влажностью. Их изменение может привести к скачкообразному изменению плодородия почвы, и это мы смоделировали, вводя в правую часть уравнения, описывающего динамику плодородия почвы, функции, соответствующие катастрофам типа «сборка» (см. рис. 1).
Модель представляет собой дифференциальное уравнение вида:
ddy = Y ■ [(Р - Ро) - y2]y - S ■ (W - W_)(W - W+), (1)
0 <W- < W+,
где y — мера плодородия почвы, p — мера типа почвообразующей породы, W — влажность почвы, W- — значение влажности почвы, которое характеризует нехватку воды, и, соответственно, W+ — ее избыток, у, S — положительные константы.
В результате для определённых почв в точке (po,W-) происходит катастрофа падения (от точки С к точке D) плодородия, связанная с нехваткой воды в почве, а в точке (p0,W+) — катастрофа падения (от точки А к точке В) плодородия при избытке влаги.
2. Описание взаимосвязи «растительность-почва» в лесных экосистемах
Почва относится к числу основных факторов, определяющих условия произрастания деревьев. Растительность и почвы связаны потоками энергии и вещества и совместно с животными и микроорганизмами формируют целостный
58 A.K. Гуц, Л.А. Володченкова. Модель системы «растительность-почва»
Рис. 1. Слева изображена катастрофа плодородия почвы, вызванная нехваткой влаги в почве,
а справа — избытком влаги
лесной биогеоценотический покров. В [1] была предложена ярусно-мозаичная модель лесной экосистемы в виде дифференциального уравнения для продукции фитомассы с правой частью, описываемой в рамках математической теории катастроф. Взаимосвязь с почвой учитывалась посредством введения в правую часть в качестве внешнего параметра влажности почвы.
Предлагается следующая модель взаимосвязи «растительность-почва», расширяющая модель четырёхъярусного леса, предложенную в [1],
dx д
~dt = - 0xV (x,k,m,a,w), dt = Y ■ [(p - po) - y2]y - S ■ (W - W_)(W - W+),
V(x, k, m, a,w) = — (x - x0)6 + k(x - x0)4 + m(x - x0)3 + a(x - x0)2 + w(x 6
k = -ck(CI - CIo), m = cm(--1),
У
a = -ca(УАН - УАН0), w = Aw(W - W0),
2
(2)
xo),
0 < W- < Wo < W+,
где x — продукция фитомассы (т/га за год), у — мера плодородия почвы, CI — индекс конкуренции Вайса [1]; s2/у — коэффициент дисперсии, являющийся показателем равномерности распределения деревьев в пространстве; если s2 /у близко к нулю, то распределение регулярное, к единице — случайное, а чем
Математические структуры и моделирование, 2015. №3(35)
59
более единицы, — тем мозаичнее; УАН — уровень антропогенной нагрузки на район [1], p — мера типа почвообразующей породы, W - влажность почвы, W-— значение влажности почвы, которое характеризует нехватку воды, и, соответственно, W+ — ее избыток, 7,5 - положительные константы, коэффициент a = а1а2аза4, где aj — доля фитомассы j-го яруса в фитомассе всего леса, ck,cm,ca,cw — постоянные коэффициенты.
Величины CI0, yAH0,Wo — это критические значения факторов, обозначающие границы экологической устойчивости фитоценоза [1].
Рассматривая стационарные равновесия для у, находим W = W(y,p). Подставляя это в первое уравнение, концентрируем внимание либо на равновесных состояниях для продукции, либо на её динамике.
Первое уравнение описывает четырёхъярусный лес. Добавление второго уравнения таким способом, как это мы предлагаем, сохраняет все результаты, полученные нами для теоретико-катастрофического описания четырёхъярусного леса и подробно изложенные в [1].
Второе уравнение — это уравнение (1), введённое в § 1, реализующее упрощённое представление о плодородии почвы и учитывающее только два фактора: тип почвообразующей породы и влажность почвы. Их изменение может привести к скачкообразному изменению плодородия почвы, и это мы смоделировали, вводя в правую часть уравнения катастрофу типа «сборка». В точке (p0,W-) происходит катастрофа падения плодородия, связанная с нехваткой воды в почве, а в точке (p0, W+) — катастрофа падения плодородия при избытке влаги. 3
3. Пример двойной экологической катастрофы
Переувлажнение почвы, как видно из рис. 1, влечёт катастрофу (типа «сборка») плодородия, т.е. величина у скачком переходит к меньшему значению у. В при этом мы будем наблюдать, как было показано в [5], скачкообразный переход (катастрофа типа «бабочка») лесной экосистемы к новому равновесию, характеризующемуся крайне низкой продуктивностью (рис. 2).
Рис. 2. График потенциала V(x) при сильной влажности почвы W = W0 + 60% (лес стоит в
воде).
60 A.K. Гуц, Л.А. Володченкова. Модель системы «растительность-почва»
Литература
1. Гуц А.К., Володченкова Л.А. Кибернетика катастроф лесных экосистем, Омск : Изд-во КАН, 2012. 220 с.
2. Володченкова Л.А. Модель плодородия почвы с точки зрения катастрофы «сборка» // Математическое и компьютерное моделирование: сборник материалов Международной научной конференции (Омск, 21 ноября 2014 г.). Омск : изд-во Ом. гос. ун-та, 2014. С. 25-26.
3. Гуц А.К., Володченкова Л.А. Теоретико-катастрофическая модель взаимосвязи «Растительность-почва» в лесных экосистемах // Математическое и компьютерное моделирование: сборник материалов Международной научной конференции (Омск, 21 ноября 2014 г.). Омск : изд-во Ом. гос. ун-та, 2014. С. 23-24.
4. Солодова Н.И., Христофорова Л.Б., Малахова Н.И. Основы устойчивого управления лесным хозяйством. СПб. 2006. 111 с.
5. Гуц А.К., Володченкова Л.А. Динамика вымокающего лесного фитоценоза // Вестник Омского университета. 2013. Вып. 4. С. 19-22.
mathematical model of interrelation "vegetation-soil"
in forest ecosystems
A.K. Guts
Dr.Se.(Phys.-Math.}, Professor, e-mail: guts@omsu.ru L.A. Volodehenkova
Senior Lecturer, e-mail: Volodehenkova2007@yandex.ru Omsk State University n.a. F.M. Dostoevskiy
Abstract. In the article the mathematical model of fertility of the soil based on the differential equation allowing the cusp catastrophes at change of soil humidity is offered. Also the model of an interrelation of phytomass production for forest and soil fertility is given.
Keywords: model, soil, vegetation, forest ecosystem, theory of catastrophe, cusp catastrophe.
