Четырехкомпонентная (N–P–Z–D) динамическая модель планктонной пищевой цепи
УДК 577.3:574.583
К. А. Подгорный
ИБВВ им. И. Д. Папанина РАН, Борок, Россия kap@ibiw.yaroslavl.ru
A FOUR-COMPONENT
(NUTRIENT–PHYTOPLANKTON–ZOOPLANKTON–DETRITUS)
DYNAMIC MODEL OF PLANKTON FOOD CHAIN
K. A. Podgornyj
I. D. Papanin Institute for the Biology of Inland Waters of the Russian Academy of Sciences,
Borok, Russia, kap@ibiw.yaroslavl.ru
Рассматривается пространственно однородная модель планктонной пищевой цепи, описывающая временную динамику концентраций биогенных элементов (N), фитопланктона (P), зоопланктона (Z) и планктонного детрита (D). Модель представляет собой систему четырех обыкновенных дифференциальных уравнений:
(1)
где – максимальная скорость роста фитопланктона; b – коэффициент поглощения света водой; c – коэффициент самозатенения фитопланктона; e, – константы полунасыщения; r – скорость дыхания фитопланктона; – максимальная скорость роста зоопланктона; q – коэффициент выедания зоопланктона организмами высших трофических уровней; k – скорость поступления биогенных элементов; – концентрация биогенных элементов, поступающих в систему (если рассматривается верхний квазиоднородный слой океана, то – это концентрация биогенов ниже слоя перемешивания); s, – скорости оседания фитопланктона и детрита соответственно; – скорость минерализации детрита; – порядок ферментативной реакции ; – коэффициент эффективности потребления фитопланктона зоопланктоном; – коэффициент экскреции пищи зоопланктоном; – коэффициент регенерации пищи, не усвоенной организмами высших трофических уровней; – скорость выделения метаболитов (токсических компонентов) клетками фитопланктона. Все переменные модели неотрицательны и выражены в граммах углерода на м3. Параметры модели неотрицательны, причем , и . Последнее условие означает, что фитопланктон не вымирает полностью при очень низкой концентрации биогенных элементов.
В отличие от модели (Edwards, 2001) в данной модели убыль зоопланктона описывается не квадратичной, а более простой линейной зависимостью вида . Для описания выедания фитопланктона зоопланктоном используется функция Холлинга II типа вместо функции Холлинга III-го типа . Важным отличием данной модели является использование нелинейной зависимости при описании процесса минерализации детрита. Детрит представляет собой не индивидуальное соединение, а смесь различных по своей природе групп веществ, которые имеют разную реакционную способность по отношению к распаду. Поэтому суммарная скорость процесса, описывающая распад всех компонентов, зависит от их соотношения в смеси, а порядок ферментативной реакции может изменяться в достаточно широких пределах (Долгоносов, Губернаторова, 2005). Еще одной отличительной особенностью данной модели является включение в нее слагаемого , которое позволяет учитывать отрицательное влияние метаболитов фитопланктона на развитие зоопланктона.
Цель работы состояла в аналитическом и предварительном численном исследовании динамического поведения модели в зависимости от изменения порядка ферментативной реакции, описывающей разложение детрита, а также скорости выделения метаболитов клетками фитопланктона.
Для системы уравнений (1) возможны следующие стационарные состояния: , , . Если , где , то стационарное состояние – устойчивый узел, в противном случае – неустойчивая особая точка типа седло. Для обеспечения устойчивости состояния необходимо одновременное выполнение двух неравенств: и . Устойчивость особой точки определяется следующими условиями: , , , где
а – элементы матрицы Якоби для системы (1).
Численные эксперименты с моделью (1) проводились методом Рунге-Кутта-4 при следующих значениях параметров: гС/м3; (м сут)–1; м–1; м2/гС; сут–1; ; сут–1; гС/м3; ; сут–1; сут–1; гС/м3; сут–1; ; м3/(гС сут); сут–1; ; сут–1. Ниже для примера приведены результаты четырех вариантов вычислений.
Вариант 1. . – узел; – устойчивый фокус; – предельный цикл; – бифуркация удвоения периода, появление двухобходного предельного цикла; – серия бифуркаций удвоения периода; – хаотические колебания; – возврат к периодическому режиму, но с более сложной, чем прежде, структурой колебаний. При этом возрастают как амплитуда, так и период колебаний; – предельный цикл.
Вариант 2. . – узел; – устойчивый фокус; – предельный цикл; – бифуркация удвоения периода, появление двухобходного предельного цикла; – серия бифуркаций удвоения периода; – хаотические колебания; – возврат к периодическому режиму, но с более сложной структурой колебаний, возрастанием амплитуды и периода.
Вариант 3. . – узел; – устойчивый фокус; предельный цикл; – бифуркация удвоения периода, появление двухобходного предельного цикла; – серия бифуркаций удвоения периода; – хаотические колебания.
Вариант 4. . – узел; – устойчивый фокус; предельный цикл. Во всех четырех вариантах при значениях происходит вымирание зоопланктона и переход системы из состояния в .
Таким образом, видно, что при возрастании порядка ферментативной реакции, описывающей разложение детрита в планктонной системе, происходит увеличение характерного времени распада детрита и уменьшение возможных типов динамического поведения системы. В этом смысле можно говорить о том, что детрит представляет собой важнейший стабилизирующий элемент экосистемы. Влияние метаболических выделений фитопланктона на зоопланктон также представляет собой важный элемент регулирования поведения планктонной системы. За счет этого влияния возможно ослабление пресса зоопланктона на фитопланктон, перераспределение потока вещества в системе и направление его по цепочке . В результате чего фитопланктон получает определенные преимущества для своего развития. Это влияние особенно велико при уменьшении порядка ферментативной реакции .
Рассмотренная в данной работе модель не учитывает того, что зоопланктон на самом деле может питаться не только фитопланктоном, но и детритом. При этом у зоопланктона существуют определенные пищевые предпочтения в питании, которые могут изменяться с течением времени в зависимости от тех или иных условий. Рассмотрение этого круга вопросов является предметом дальнейших исследований.
Zoocenosis — 2009
Біорізноманіття та роль тварин в екосистемах: Матеріали V Міжнародної наукової конференції. – Дніпропетровськ: Ліра, 2009. – С. 80-82.