по истории и философии науки История изучения динамики популяций (стр. 2 из 5)

пробрели всемирную известность и были многократно воспроизведены в разных модификациях на компьютерах. Красота фрактальных изображений завораживает. (См. "Красота фракталов", М., 1995 - перевод с англ. книги: H.-O.Peitgen, P.H.Richer " The Beauty of Fractals",Springer, 1986).

К дискретному логистическому уравнению мы обратимся позднее, а сейчас вспомним тот биологический факт, что в природе популяции имеют не только максимальную численность, определяемую величиной экологической ниши K, но и минимальную критическую численность L. При падении численности популяции ниже этой критической величины из-за неблагоприятных условий, или в результате хищнического промысла, восстановление популяции становится невозможным.

Величина нижней критической плотности различна для разных видов. Исследования биологов показали, что она может составлять всего лишь пару особей на тысячу квадратных километров в случае ондатры, и сотни тысяч особей для американского странствующего голубя. Заранее трудно было предположить, что столь многочисленный вид уже перешел через критическую границу своей численности и обречен на вымирание. Например, для голубых китов критическая граница численности оказалась равной десяткам - сотням. Хищническое истребление этих гигантских животных привело к тому, что их осталось слишком мало в Мировом океане. И хотя охота на них давно запрещена, надежд на восстановление популяции голубых китов практически нет. Кривые показателей численности для трех видов китов приведены на рис. 2.

Рис.2. Динамика численности трех видов китов в мировом океане. По оси ординат отложен индекс численности - число убитых китов на 1 тыс. судо-тонно-суток. (Gulland, 1971)

Модели, описывающие как внутривидовую конкуренцию, определяющую верхнюю границу численности популяции, так и нижнюю критическую численность популяции, имеют два устойчивых стационарных решения. Одно из них - нулевое для начальных численностей, которые ниже наименьшей критической численности популяции. Другое равно K - емкости экологической ниши в случае, когда начальная численность выше наименьшей критической величины. Такими "триггерными" свойствами обладает нелинейное уравнение, предложенное А.Д.Базыкиным [1]

(6)

В формуле (6) первый член в правой части описывает размножение двуполой популяции, скорость которого пропорциональна квадрату численности (вероятности встреч особей разного пола) для малых плотностей, и пропорциональна числу самок в популяции - для больших плотностей популяции. Второй член описывает смертность, пропорциональную численности, а третий - внутривидовую конкуренцию, подобно тому, как это было в логистическом уравнении (4).

Зависимости численности от времени и скорости прироста от численности представлены на рис. 3 (а,б). Кривые 1,-5 соответствуют различным начальным численностям. x=0 и x=K -устойчивые стационарные состояния, x=L -неустойчивое, разделяющее области влияния устойчивых состояний равновесия. Величины L и K различны для разных популяций и могут быть определены из наблюдений и экспериментов.

Из рисунка 3а видно, что скорость восстановления популяции после ее падения в силу промысла или неблагоприятных условий зависит от того, насколько близка новая начальная численность к опасной границе L. Если ущерб, нанесенный популяции невелик (меньше половины емкости экологической ниши) популяция быстро восстанавливается по кривой 1, не имеющей точки перегиба. В случае, когда численность оставшейся популяции близка к критической, восстановление происходит сначала очень медленно, популяция надолго "застревает" вблизи опасной границы", а затем уже, "набрав силы", более быстро выходит на устойчивый стационарный уровень K (кривая 3). Кривая 2 представляет промежуточный случай. Кривые 4, 5 иллюстрирует вырождение популяции в случае, когда начальная численность опустилась ниже критической границы. обращает на себя внимание сходства начальных участков кривых 3 и 5. Близость к опасной границе со стороны больших значений (3) и меньших (5) выражается в долгом пребывании системы в неопределенном состоянии, когда малые флуктуации могут легко "перебросить" систему через опасную границу в "благополучную" область возврата к стационарному значению K, или, наоборот, в область вымирания. В это время сторонний наблюдатель не сможет определить по форме кривой динамики численности, какая судьба ожидает систему. Для самих участников жизненной драмы - нахождения системы вблизи опасной границы - исход не очевиден. Важно понимать, что в этой ситуации чрезвычайно важны любые, даже очень малые усилия, направленные на преодоление критического барьера.

Именно популяции, численность которых близка к нижней критической численности, занесены в Красную книгу. Удастся ли перенести каждый конкретный вид на "Зеленые страницы", куда переносят виды, исчезновение которых удалось предотвратить, зависит от многих обстоятельств, в частности как от репродуктивных усилий вида, так и от усилий людей, спасающих эти виды.

Влияние запаздывания.

Уравнения, которые мы рассматривали до сих пор, предполагают, что процессы размножения и гибели происходят одновременно и популяция мгновенно реагирует на любое изменение внешних условий. Однако в реальности это не так. Всегда имеется некоторое запаздывание, которое вызвано несколькими причинами.

Развитие любой взрослой особи из оплодотворенного яйца требует определенного времени T. Поэтому если какое-нибудь изменение в окружающей среде, например, увеличение ресурса, вызовет внезапное повышение продуктивности взрослых особей, то соответствующе изменение численности произойдет лишь по прошествии времени T. Это означает, что уравнение

dx/dt=f(x), (7)

где x - численность взрослых особей, следует заменить уравнением:

dx/dt=f(xt--T) (8)

где xt--T - численность половозрелых особей в момент t-T.

В реальных популяциях интенсивность размножения и гибели различны в разных возрастных группах. Например, у насекомых откладывают яйца взрослые особи, а конкуренция наиболее выражена на личиночной стадии. Такие процессы, как отравление среды продуктами метаболизма, каннибализм и т.п. в наиболее сильной степени воздействуют на ранние возрастные стадии, а их интенсивность зависит от численности взрослых особей, т.е. отрицательное влияние на коэффициент естественного прироста оказывают особи предыдущего поколения. С учетом этих обстоятельств, логистическое уравнение (4) перепишется в виде:

dx/dt=x(r-xt-T) (9)

Наиболее распространенное и изученное в динамике популяций уравнение Хатчинсона учитывает тот факт, что особи размножаются лишь с определенного возраста, и имеет вид:

(10)

Смысл модели (10) заключается в том, что уровень лимитирования системы зависит не только от общей численности популяции в данный момент времени t, определяемой емкостью среды, но и от количества половозрелых особей в момент времени t-T. Еще более точное уравнение, учитывающее распределение времени запаздывания:

Вид функции распределения времен запаздывания w(t-s) представлен на рис. 4. Такого типа уравнения могут иметь колебательные решения. Это легко проверить для простого линейного уравнения:

которое имеет периодическое решение

в широком диапазоне значений скоростей роста r и времени запаздывания T.

В технике хорошо известно, что запаздывание в регуляции системы может привести к возникновению колебаний переменных. Если система регулируется петлей обратной связи, в которой происходит существенная задержка, то весьма вероятно возникновение колебаний. Если продолжительность задержки в петле обратной связи больше собственного времени системы, могут возникнуть колебания с нарастающей амплитудой, нарушаются их период и фаза.

Рис.4. Модель динамики популяции с учетом распределения времен запаздывания. Типичный вид весовой функции w(t)

Дискретные модели популяций с неперекрывающимися поколениями

Даже в таких популяциях, где особи размножаются несколько лет подряд (млекопитающие и птицы, многолетние растения), наличие сезонов размножения вносит некоторое запаздывание в процессы регуляции численности. Если же взрослые особи, размножающиеся в данном году, редко или никогда не доживают до того, чтобы размножиться в будущем году, как, например, у однолетних растений, мелких грызунов, многих насекомых, это оказывает существенное влияние на динамику их численности. В этом случае уравнение (7) следует заменить уравнением

Nn+1=N(xn), (11)

где Nn - численность популяции в году n.

Наблюдения над динамикой численности показывают, что в таких системах при малых численностях N растет от одной генерации к другой, а при высоких - падает. Это свойство - резко расти при малых N и падать при больших, проявляется в экономике как закон "бумов и спадов". В таких случаях функция F - одноэкстремальная, вид ее изображен на рис. 5а.

Функция такого типа может быть описана с помощью различных формул. Ниболее широко распространена версия дискретного логистического уравнения, предложенная Мораном для численности насекомых (1950) и Рикером для рыбных популяций (1954):