Большие системы. Связность, сложность и катастрофы - Касти Дж.
Скачать (прямая ссылка):
короткие промежутки времени значительные колебания биомассы фитопланктона
в небольших эвтрофных озерах. Поскольку некоторые особенности наблюдаемых
изменений фитоплактона свидетельствуют в пользу возможности применения
методов теории катастроф, ниже предлагается динамическая модель этого
процесса с использованием катастрофы типа сборки.
Будем считать общую биомассу ') всех видов планктонных водорослей в озере
выходной переменной х системы. В большинстве озер, особенно в случаях
доминирования сине-зеленых водорослей, лимитирующим биогенным элементом
является фосфор. Исходя из этого, концентрацию растворимого фосфата2)
следует включить в модель как один из параметров управления аь В качестве
эквивалента уровня освещенности и температуры в озере будем рассматривать
изменения одного из представителей фитопланктона - род Anabaena.
Наблюдения показывают, что отмиранию Апа-Ьаепа предшествует ее скопление
на поверхности, вызванное способностью этого вида к образованию газовых
вакуолей и поднятию на поверхность в условиях небольшой освещенности.
Более того, колебания концентрации растворенного фосфата оказываются
противоположными колебаниям концентрации фитопланктона. Таким образом, в
качестве второго параметра управления выбираем биомассу1) Anabaena.
') В единицах концентрации хлорофилла, мкг/мл. - Прим. ред s) В единицах
концентрации фосфора, мкг/мл. - Прим. ред.
Устойчивость, катастрофы и адаптируемость больших систем 171
Для вывода уравнения, описывающего изменения концентрации фосфата,
используется тот факт, что скорость накопления фосфора планктонными
водорослями пропорциональна концентрации последних во время "цветения"
воды. Когда наступает отмирание, фосфат осаждается на дно вместе с
клетками и может быть "освобожден" только после окончания периода
"цветения" и следующего за ним процесса отмирания. Таким образом,
концентрация растворенного фосфата будет по экспоненте приближаться к
равновесному уровню для данного озера.
Чтобы получить уравнение, определяющее концентрацию Anabaena, заметим,
что скорость роста Anabaena будет зависеть только от численности ее живых
клеток, в то время как снижение характерной для этого вида скорости роста
будет зависеть также и от концентрации клеток конкурирующих видов,
потребляющих тот же биогенный элемент. В конечном итоге уравнение для
общей концентрации фитопланктона вытекает из логистической модели,
которая подробно описывается в следующем разделе.
Учитывая все сказанное выше и используя каноническое уравнение сборки для
выходной переменной х, получим следующую систему уравнений, моделирующую
рассматриваемый процесс:
х = - (с^3 - с2щх + с3а2),
", = - с4х (а, - а0), а2 = с5а,а2 - с6а 2х,
где с .....с6-константы скорости и а0 - равновесная кон-
центрация фосфата.
С учетом экспериментальных данных, полученных в шт. Алабама (США) при
исследовании пруда для разведения зубатки, приведенные выше уравнения
были численно проинтегрированы для описания цикла эвтрофикации. При этом
использовались следующие значения констант:
с, = 3,10 (мкг хлорофилла/мл)-2(время)-1, с2 = 0,60 (мкг фосфора/мл)-1
(время)-1, сг = 0,05 (время)-1,
с4 = 1,00 (мкг хлорофилла/мл)-1 (время)-1, с5 = 1,35 (мкг фосфора/мл)-1
(время)-1, с6 = 1,95 (мкг хлорофилла/мл)-1 (время)-1, а0 == 0.
Начальными значениями были *о = 0,11 мкг хлорофилла/мл, а°=1,28 мкг
фосфора/мл, а^ = 0,02 мкг хлорофилла/мл.
172 Устойчивость, катастрофы и адаптируемость больших систем
Результаты эксперимента показаны на рис. 5.8-5.10. Для концентрации
Anabaena и общей концентрации фитопланктона данные эксперимента
достаточно хорошо согласуются с приведенной выше моделью, за исключением
большого спада в пятый интервал времени. Этот спад соответствует
возрастанию средней силы ветра в соответствии с измерениями и, возможно,
является следствием усиленного вертикального
0 1 2 3 4 5 Я 7 8 9 10 11 12
Время (за единицу дыбран период В 1,75 дня)
Рис. 5.8. Траектории переменных управления и переменных состояния для
суммарной концентрации (х) водорослей [Parks et al., 1975].
•-• результаты наблюдений; X-X теоретические значения.
или горизонтального перемешивания, в результате которого в месте взятия
проб клетки могли переместиться из поверхностного слоя пруда.
Таким образом, моделирование дает возможность получить детерминированный
диапазон, в пределах которого на цикл цветение - отмирание водорослей
могут влиять случайные явления, например погода.
Совпадение теоретических и экспериментальных данных о концентрации
растворенного фосфата оказывается менее точным. Подобное расхождение
может быть вызвано несколькими причинами. Например, при измерении
концентрации фосфата могли не учитываться растворенные удобрения из-за
временных и пространственных разрывов в измерениях, которые проводились
только раз в 1-3 дня и только в одной точке пруда. Кроме того, модель
может быть
Время (за единицу Выбран период 61,75 дня)
Рис. 5.9. Траектории переменных управления и переменных состояния для
