СЦЕНАРНАЯ МОДЕЛЬ ЭФФЕКТА ВРЕМЕННОГО РЕЗКОГО СОКРАЩЕНИЯ ЧИСЛЕННОСТИ ПОПУЛЯЦИИ С БОЛЬШИМ РЕПРОДУКТИВНЫМ ПАРАМЕТРОМ

Полный текст

Введение

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

Вестник Самарского университета. Естественнонаучная серия. 2021. Том 27, № 2. С. 80–90

Vestnik of Samara University. Natural Science Series. 2021, vol. 27, no. 2, pp. 80–90 81

Например, вселение бабочки Cydalima perspectalis в 2012 г. разрушило экосистему самшитовых рощ Черноморского побережья. Популяции видов-вселенцев очень различно ведут себя при попадании в разные новые ареалы. Например, вредоносный гребневик Mnemiopsis leidyi в Азовском и Каспийском морях демонстрировал разную и несхожую динамику. Бывает так, что вид, который в исходном ареале был малочисленный, в новой среде демонстрирует уверенный рост от малой исходной группы и даже с переходом к вспышке численности. Неудержимый рост случался с обычными европейскими речными моллюсками рода Dreissena при попадании в водоемы Северной Америки [1]. Интересно, что если в водоем попадал один из видов моллюсков, то наблюдался взрывообразный рост. Однако, если вселенцев было несколько, то итоговая их совместная биомасса оказывалась намного меньше той, которую мог бы достичь при захвате ареала один вселенец.

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

В серии предыдущих работ [2; 3] мы предложили феноменологические сценарии развития процесса расселения агрессивного чужеродного вида в модификациях уравнений с запаздыванием. Модели рассматривают случаи: № 1 — гибель агрессивного нового вида c min N∗(t; τr) = 0 при независимой регуляции, № 2 — спорадические и колебательные режимы изменений, № 3 — закрепление вселенца

на некотором допустимом малом пороге численности «бутылочного горлышка», который не угнетает среду, № 4 — разрушение циклической траектории N∗(t, rτ ) в (1) при N (tˆ) > K и появление псевдопериодической траектории в момент t > tˆ+ τ2. Данные сценарии далеко не полностью описывают

наблюдаемые ситуации при инвазиях [4]. Например, при распространении вредителей кукурузного мотылька и непарного шелкопряда в Северной Америке [5] наблюдаются резко возникающие и медленно затухающие пилообразные осцилляции к минимальному значению численности limt→∞ max N∗(t, rτ ) = L. Одной из нерешенных проблем является получение в модели инвазивного сценария № 5 — уверенного прохождения на первом этапе инвазии глубокого кризиса с повторной фазой быстрого роста численности. Цель работы — разработка новой специфической модели переходной кризисной динамики с быстрым восстановлением. Актуальность работы обусловлена задачей поиска способов биологического контроля опасных вселенцев, разработки методики целенаправленной интродукции естественных врагов. Задача сложна тем, что быстро или через некоторое время атакуемая популяция приспосабливается к врагу, вырабатывает ответную адаптацию к паразиту или иммунитет к патогенному микроорганизму. Время реакции среды и адаптации вида — важные характеристики процессов, которые можно отразить феноменологическим уравнением. Новизной работы будет разделение величин запаздывания по типам в регуляции в предложенной модели. Данный вариант кризиса существенно отличается от рассмотренного нами ранее сценария промысловой деградации популяций осетровых рыб Каспийского моря в гибридной

модели [6].

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

dN (

dt = r1N 1 −

N (t − τ1))

K

(H − N (t − τ2)) . (1)

Траектория при увеличении амплитуды неустоявшихся колебаний не будет в итоге притягиваться к замкнутому подмножеству фазового пространства, как показано на рис. 1, где не происходит изменений параметров по ходу вычислительного эксперимента. В декабре подобная динамика наблюдается в динамике пандемии COVID. После циклической динамики из серии трех эпидемических волн со сравнимыми максимумами в середине 2021 года число случаев заражения COVID в день стремительно увеличилось (рис. 2.) Так наша гипотетическая модель получила реальную интепретацию.

1. Стабильность, колебания и запаздывание в моделях популяций

   Цикличность в биоситемах — интересный и все еще загадочный феномен, который наблюдается экспериментально и для изолированных популяций. Обычно используют насекомых или аквариумных

   Переварюха А.Ю. Сценарная модель эффекта временного резкого сокращения численности популяции...

   82Perevaryukha A.Yu. Scenarios model of the effect of a temporary sharp reduction of population with a large...

    

   

   Рис. 1. Разрушение релаксационного цикла N∗(t) → ∞ в (1) r1 = 0.000033

   Fig. 1. Destruction of the relaxation cycleN∗(t) → ∞ в (1) r1 = 0.000033

    

   

    

   Рис. 2. График динамики заболеваемости COVID по регистрируемым случам во всех странах Fig. 2. Graph of the dynamics of COVID incidence by reported cases in all countries

    

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

   Гидробиолог и «отец американской лимнологии» Джордж Эвелин Хатчинсон один из первых начал применять математические методы в экологии и в [7] указал на эффект влияния от предшествующих состояний, a Pоберт Мэй позднее выписал феноменологическое расширение известной модели Ферхюлста

   N˙ = rN (1 − N/K) с t − τ :

   dN (

   −

    

   = rN 1

   dt

   N (t − τ )) K

    

   . (2)

   В уравнении (2) возникает бифуркация рождения цикла N∗(t; rτ ) при rτ > π/2, исследованная

   во многих работах [8]. Емкость насыщения экологической ниши K тут играет роль точки центра: ∀N (0), | max N∗(t) − min N∗(t)| ≈ K [9]. Изначальная интерпретация параметра K как

   уравновешивающего предела ниши в других модификациях может меняться. Часто уравнение исследуют

   в однопараметрической форме: x˙

   = λx(t)(1 − x(t − τ )). Включать запаздывание rN (t − τ )f (N (t); K)

   не оправданно из экологических соображений. Для модели Xaтчинcонa есть несколько обобщений [10]:

   dN (

   −

    

   = rN 1

   dt

   c1N (t − τ ) + c2N (t − τ1))

   K

    

   . (3)

   Модель Гoпoлсэми [11] с вариативной регулирующей емкостью от уровня K интересна для моделирования серии пиков вредителей:

   dN

   = rN (t)

   dt

   ( K − N (t − τ ) )

   (K + ςN (t − τ ))

    

   . (4)

   Модели бывают схожими по поведению их решений, но отличными по принципам интерпретации. Идея регуляции требует включения в (4) фиксированной характеристики — емкости насыщения

   Вестник Самарского университета. Естественнонаучная серия. 2021. Том 27, № 2. С. 80–90

   Vestnik of Samara University. Natural Science Series. 2021, vol. 27, no. 2, pp. 80–90 83

   экологической ниши, которая теоретически может достигаться популяцией ∀t → ∞. Подобный подход применим далеко не ко всем ситуациям. Нельзя определить, что представляет собой балансовая емкость N = K, f (N ) = 0 в контексте взрывообразно развивающейся инвазии вселенца с неконтролируемым репродуктивным потенциалом в новой среде. Актуальный пример — острая вирусная инфекция. Вирус в организме иногда способен захватывать «нишу» и проникать во все доступные клетки, но тогда болезнь быстро заканчивается летально и до того как будут поражены абсолютно все клетки. Подобное верно для многих ситуаций, когда предел K-емкости может быть достигнут, но только на ограниченное

   время. Дo достижения состояния N (t) → K − ϵ раньше начинаются необратимые изменения. Другой

   важный аспект: чтобы достичь предельно возможного K, N (0) должно быть выше порога J . Важнее

   не предельная теоретическая емкость K, а положение неразрушительного порога J .

   Модель с регуляцией rf (N (t−τ )) без участия K-ёмкости ниши в форме известной в математической биологии функции Mорана–Pикеpa f (x) = rxe−bx предложена при анализе лабораторных экспериментов

   A. Hиколсoнa:

   dN

   

   dt = rN (t − τ ) exp(−bN (t − τ )) − δN (t). (5)

   Модель при увеличении rτ демонстрирует релаксационные колебания [12].

   У метода получения осциллирующих решений в моделях c t − τ есть существенная проблема с точки зрения экологической обоснованности поведения решения. Потому простое включение t−τ в уже

   известные модели выглядит не совсем биологически убедительно при увеличении значений ключевых параметров rτ .

   Проблемы возникают для корректного описания в уравнениях c t−τ колебаний с большой амплитудой

   численности. Если увеличить rτ в (4), то цикл быстро приобретет форму негармонических колебаний.

   Релаксационный цикл с увеличивающейся амплитудой по мере возрастания rτ получает вид все дальше отстающих друг от друга во времени Λ-образных пиков max N∗(t) ≫ K. Tраектория в минимумах

   цикла начинает сколь угодно близко приближаться к нулю [13]. Известно [14], что для модели (2) (в однопараметрической форме) выполняется:

   ( λ 1 + (1 + λ) ln λ

   [ ln2 λ ] )

   min N∗(t, λ) = exp

   −e + 2λ +

   

   λ + O λ2

   − 1 .

   Окрестности минимумов колебаний, где траектория пробегает «околонулевые» значения, становятся чрезвычайно длинными. Для популяционной интерпретации эта картина нереалистична, потому актуальна специфическая модификация модели для краткого популяционного кризиса.

   В [15] рассматривалась модель колебаний:

   dN

   

   dt = λN (t)f (N (t − 1)) , (6)

   где предполагается λ ≫ 1, a f (x) — дифференцируемая функция, разложимая в асимптотический ряд,

   для которой выполняются:

   ∞

   f (0) = 1, f (x) = −a0 + ∑

   k=1

   ak

   

   x

    

   , x → ∞, a0 > 0. (7)

   k

   Свойствам (7) соответствует функция регуляции из уравнения Гополсэми (4):

   (1 − x)

   f (x) =

   (1 + ςx)

   . (8)

   ς в (4) тоже влияет на характеристики релаксационного цикла.

   Из дополнительных ограничений при преобразовании исходного уравнения (4) к сингулярно

   возмущенному

   x˙ = F (x(t − 1), ϵ) и к предельному релейному уравнению

   x˙ = R(x(t − 1)) для (4) c (7)

   было подтверждено существование единственного орбитально устойчивого цикла N∗(t; λ). Однако из

   утверждения о характеристиках колебаний в [15] получается:

   min N∗(t, λ) ∼ C1 exp (−λa0), C1, a0 = const > 0.

   Соответствие свойств цикла популяционному процессу проблематично, если λ изначально указывается достаточно большим значением. Проблема популяционной интерпретируемости циклического режима остается актуальной.

   Серийные вспышки численности вселенцев нельзя описать при свойстве минимумов:

   ∗

    

   min0<t<T

   N∗(t, r) → ϵ, и ϵ ≪ 1, слишком низких для быстрого повторного достижения экстремально

   высоких значений численности. Известны экстремальные популяционные колебания y арктических

   зайцев и грызунов — «лемминговые циклы» с большой амплитудой и малыми минимумами, но это неустойчивые режимы существования и цикличности нарушаются при неблагоприятном внешнем возмущении.

   Переварюха А.Ю. Сценарная модель эффекта временного резкого сокращения численности популяции...

   84Perevaryukha A.Yu. Scenarios model of the effect of a temporary sharp reduction of population with a large...

    

2. Модификация с квадратичной функцией противодействия биотического окружения

   Существует концепция минимальной критической численности группы, которая теоретически необходима для выживаемости локальной популяции. Локальные группы леммингов изредка действительно погибают. В известной модели Бaзыкинa [16] введен L-порог непосредственно как

   априори существующее ∀N (0), t неустойчивое равновесие ∀N (t) < L, limt→∞ N (t) = 0:

   dN (

   

   −

    

   = rN 1

   dt

   N )

   

   K (N − L). (9)

   В эту версию бистабильной логистической модели можно было бы включить и запаздывание (N (t −

   — τL) − L), но тут включение не получит разумной сущностной интерпретации. В реальности оценить

   заранее величину L-порога практически невозможно. Несколько малочисленных реликтовых видов

   существуют миллионы лет. Многие агрессивные инвазии начинались с единичных проникших особей.

   Для начала инфекционного процесса достаточно одного проникшего в организм вириона ВИЧ. Если верна гипотеза, для популяции критически опасно N (t) → L, то жесткая квадратичная саморегуляция rf (N 2) в уравнении (9) c (N (t − τL) − L) явно выглядит избыточной. Многомиллионные колонии

   социальных насекомых не страдают от высокой плотности, но минимально необходимое количество

   рабочих насекомых действительно известно для пчелиных семей. Существует менее известная модель с квадратичным фактором внешнего сопротивления −δN 2, которая лучше подходит для модификаций с N (t − τ ):

   dN

   

   dt = r1

   γN 2

    

   

   γ + σN

   − ςN − δN 2. (10)

   Изначально уравнения N˙

   = rN (t)f (N (t − τ )) предложили с целью описать колебания y популяций

   в стационарных условиях лаборатории, где τ предполагали характеристикой жизненного цикла вида.

   Интерпретация временных параметров — проблема. Часто длина жизненного цикла вида и период между пиками естественных популяционных колебаний просто несопоставимые между собой величины, что верно и для актуальной проблемы динамики волн эпидемии. Мы предлагаем рассматривать запаздывание τ как системную величину, относящуюся к многокомпонентному взаимодействию противоборствующих факторов.

   Предложенные ранее модели c f (N (t−τ )) рассматривали специфические колебательные режимы, но

   среди проблем математической экологии есть сценарии активного противоборства без периодических

   флуктуаций [17].

    

3. Моделирования фазы резкого кризиса в сценарии агрессивной инвазии

Далее мы рассмотрим новую модель для важной ситуации, которая не связана с возникновением режима осцилляций — актуального сценария преодоления внезапного кризиса в момент, когда еще не был исчерпан потенциал быстрого роста. Предложим модель на основе логарифмической формы регуляции для обоснования лабораторных экспериментов Г. Гаузе с гибелью обоих конкурирующих видов инфузорий в искусственной экосистеме после интродукции хищника [18]. В уравнение включен

отложенный фактор. В следующей форме уравнения N˙

= rN (t)f (N (t−τ ))−F (N ) мы получим сценарий

гибели вида после колебаний при убыли с запаздыванием δN (t − τ ):

dN

= N (t)r ln

dt

( K )

N (t − τ )

− δN (t − τ1). (11)

Этот вариант не совсем тот сценарий, что мы хотим описать в модели. Для f (N ) = rN ln(K/N )

точка перегиба Np кривой решения N˙

= f (N ) лежит ниже K/2, так как f′(Np) = 0, Np = K/e.

Параметр ограниченной среды в данной модели не тождествен понятию «ёмкость ниши» y Ферхюлста

или Хатчинсона. Потому мы используем в записи K. В вычислительном сценарии с (11) наблюдается гибель популяции агрессивного вселенца после двух максимумов осцилляций, как было в опытах Георга Гаyзе с инфузориями — хищником и жертвой, это неустойчивая система. При уменьшении

r-параметра траектория демонстрирует обычные затухающие N (t) → K гармонические колебания, но в

(11) пожертвовали удобным свойством ∀t, N (0) � 0N (t) � 0.

Включение нелинейности подавления δNk (t − τ1), τ � τ1 обосновано тем, что текущее воздействие

паразитов может определяться предшествующим состоянием популяции хозяина, но такое воздействие в

Вестник Самарского университета. Естественнонаучная серия. 2021. Том 27, № 2. С. 80–90

Vestnik of Samara University. Natural Science Series. 2021, vol. 27, no. 2, pp. 80–90 85

свою очередь связано зависимостью с концентрацией доступных для заражения жертв. Формализованное в модели воздействие может в реальности отражать целенаправленное биологическое подавление вредных видов, которое осуществляют при выпуске их паразитических естественных врагов, например наездников семейств Braconidae, Ichneumonidae, которых специально разводят. Против непарного шелкопряда Lymantria dispar применяют распыление энтомопатогенных грибков Beauveria bassiana. Давление паразитов-наездников на популяцию фитофагов зависит от количества доступных им жертв в прошлом поколении, что отличается от хрестоматийного варианта «хищник–жертва». Однако методами биологического контроля полного подавления достичь удается крайне редко, так как их воздействие регулируется более сложным образом. Большинство жертв в ходе коэволюции вырабатывают средства адаптации к паразиту, потому паразит должен эволюционировать быстрее. Известно из наблюдений, что эффективность атак паразитов возрастает при скученности жертв [19], а при малой численности жертвы уходят от обнаружения. Патогенные микроорганизмы аналогично интенсивно распространяются только при повышенной плотности.

Существуют интересные сценарии экстремальной популяционной динамики с пороговой регуляцией. В экосистемах с высоким индексом биоразнообразия наблюдаются сложные многоступенчатые механизмы регуляции численности видов. Длительно доминирующим даже при высоком репродуктивном потенциале ни один из видов-конкурентов стать не может. Тем не менее иногда наблюдаются нерегулярные короткие вспышки массовых размножений у уязвимых для естественных врагов вредителей в вечнозеленом лесу с очаговой дефолиацией леса. Псиллиды вдруг переходят к вспышке численности, хотя у них много врагов и нет химической защиты. В отряде перепончатокрылых существуют семейства ос-паразитов, которые атакуют других перепончатокрылых. Жизненный цикл паразитической осы должен быть точно приспособлен для развития личинки внутри конкретных хозяев. Действия всех паразитов специфичны, им подходят только определенные стадия развития их видов-хозяев. Обычно численность сверхпаразитов невысока при малом для насекомых r-параметре. Насекомых фитофагов эффективно регулируют паразитические наездники, но первичные паразиты сами становятся жертвами атак паразитов второго порядка Coccidoctonus psyllae, если их концентрация становится велика. Сверхпаразиты бывают даже третьего порядка, но у насекомых далеко не самая сложная система контроля. В результате межвидового соперничества паразитов у листоблошек появляется возможность для ограниченной во времени фазы роста численности, если их популяция преодолеет порог и ускользнет от контроля первичных паразитов. Фаза роста скоро сменяется резким падением. Псиллиды не уничтожают растительность непосредственно, но в липких выделениях листоблошек развиваются патогенные для поврежденных листьев деревьев микроорганизмы.

Рассмотрим ситуацию активного противодействия, которое формируется спустя некоторое время против быстро размножающейся популяции. Биотическое сопротивление в свою очередь имеет сложную регуляцию и с запаздыванием F (Nk, τ1). Через некоторое время популяция, внутренне регулируемая с собственным запаздыванием, начинает встречать сопротивление со стороны биотического окружения, и данный интервал времени тоже становится важной характеристикой процесса. Актуальный для всех нас сегодня пример — некоторые люди после контакта с вирусом эффективно мобилизуют адаптивный иммунитет и их B-лимфоциты скоро начинают вырабатывать иммуноглобулины IgG с высокой аффинностью к коротким белковым фрагментам вируса. У других иммунная активация происходит замедленно, вирус получает время на активное развитие в организме. В худшем для нас сценарии инфекция переходит в хроническое состояние с накоплением дальнейших мутаций. Механизм иммунной активации — процесс многоступенчатый и триггерный, c событийной компонентой, где не все события при выборе эпитопов и активации наивных лимфоцитов априори предопределены. В полностью непрерывной модели для сущностного описания механизм запуска ответа недоступен, потому моделирование инфекционной динамики тут останется иллюстративным. Опишем противоборство

для саморегулируемого размножения воздействующего на среду вселенца в форме уравнения N˙ =

= rN (t)f (N (t − τ )) − F (Nk (t − τ1); J ).

Представим, что развивается ситуация внезапного обострения конкурентного противоборства

для популяций при приближении к некоторому порогу численности J , что значительно меньше теоретически допустимой емкости. При приближении к порогу конкуренция при внешнем воздействии обостряется, однако уровень давления зависит от состояния популяции нелинейно. Отразим пороговый эффект противодействия активному росту численности с возникновением кризиса в уравнении при логарифмической саморегуляции:

dN

= N (t)r ln

dt

( K )

N (t − τ )

Nm(t − τ1)

2

δ

− (J − N (t))

− qN (t), δ > q, m � 2, N (0) < J < K. (12)

В вычислительном сценарии для модели (12) на первом этапе происходит увеличение численности небольшой группы N (0) < J по варианту логистической кривой, но рост останавливается до

Переварюха А.Ю. Сценарная модель эффекта временного резкого сокращения численности популяции...

86Perevaryukha A.Yu. Scenarios model of the effect of a temporary sharp reduction of population with a large...

прохождения Np. Вместо обычной стабилизации после краткого превышения равновесия «овершута» далее начинается стадия резкой убыли численности с возрастанием F (N 2; J−1) при N → J до того,

как потенциал роста будет исчерпан внутривидовой lnK -регуляцией из-за предела исчерпания ресурсов. Воздействие специфических биотических внешних, но зависимых от численности N (t − τ1) факторов

создает порог, при приближении к которому популяция переходит в демографический кризис (рис. 3).

Рис. 3. Сценарий прохождения этапа кризиса в новой модели (12) при

m = 2, r = 7.17 · 10−3, J = 103, K = 15 · 103, τ = τ1 = 48, δ = 10, q = 4.2 · 10−3, N (0) = 10.

Fig. 3. Scenario of passing through the crisis stage in the new model (12) with

m = 2, r = 7.17 · 10−3, J = 103, K = 15 · 103, τ = τ1 = 48, δ = 10, q = 4.2 · 10−3, N (0) = 10.

Популяция проходит глубокий кризис, так как реакция ее биотических противников аналогично ослабляется, и в результате популяционный процесс после переходных колебаний стабилизируется на уровне limt→∞ = P < J . Параметр K тут теоретическая предельная емкость пищевых ресурсов,

которая для изначально малой локальной группы не достигается. В иммунологической интерпретации

только очень высокая доза заражения N (0) способна стать летальной. При уменьшении коэффициента воздействия δ, от которого зависит глубина кризиса, циклический режим сохраняется. Модель

применима для исходно малой группы N (0) < J/4, при N (0) = 2/3J, m > 2 развитие покажет v-образный эффект одного из вариантов прохождения стадии «бутылочного горлышка» N (tm) ≈ 0+ϵ, limt→∞ N (t) =

= P . Очевидно, многие популяции не смогли преодолеть состояние «bottleneck» при увеличении q либо из-за внутренних причин, связанных с эффектом Олли [20]. В рассмотренном сценарии адаптируется именно внешнее воздействие как давление ос-паразитов. Для учета адаптации самой популяции (например, формирование иммунитета в группе [21]) можно использовать корректировку в

регуляции ln(√(K − N (t − ν))2/N (t − τ )). Ограничивающий показатель воздействия J, N (t) < J, ∀t, 0 <

N (0) < J в уравнении (12) зеркально симметричен по смыслу критическому порогу L из уравнения

с критическим порогом как нижней грани для существования локальной популяции 0 < N (0) < L,

limt→∞ N (t) = 0. Пороговый уровень J¯ можно включать и в составляющую логарифмической регуляции.

Выводы

Экстремальная динамика инвазивной популяции в основном зависит от двух факторов: 1) истощение необходимых ресурсов и разрушение среды; 2) активация биотического сопротивления со стороны автохтонного окружения либо из-за специально интродуцированных конкурентов или патогенов. Для популяций крупных долгоживущих животных эти два фактора могут накладываться. Полученный в модели сценарий хорошо описывает динамику локальных популяций оленей на островах Канады (рис. 4), куда они были интродуцированы. Популяция при отсутствии хищников размножается, но олени при большой скученности уничтожают растительность и разрушают свою среду. Питание становится недостаточным, смертность молодых особей увеличивается. Снижается прирост биомассы [22]. Популяция входит в кризис. Восстановление после кризиса зависит от способности к регенерации их пищевых ресурсов, что происходит с запаздыванием. Полной регенерации не происходит и так создается порог, который значительно меньше предельно допустимой емкости экологической ниши. В системе присутствуют и конкурируют два запаздывания: репродуктивный цикл оленей и темпы восстановления среды. Это создаст колебания со снижением среднего значения. Северная растительность восстанавливается медленно после вытаптывания. Интересно, что при хоте и искусственном изъятии

Вестник Самарского университета. Естественнонаучная серия. 2021. Том 27, № 2. С. 80–90

Vestnik of Samara University. Natural Science Series. 2021, vol. 27, no. 2, pp. 80–90 87

популяция оленей стабилизируется, но без изъятия флуктуации продолжаются, что и подтверждает разработанная нами модель.

Рис. 4. Динамика популяций оленей при их интродукции на острова в Канаде по данным сайта https://www.msudeer.msstate.edu/deer-habitat-carrying-capacity.php

Fig. 4. Dynamics of deer populations during their introduction to the islands in Canada according to the website https://www.msudeer.msstate.edu/deer-habitat-carrying-capacity.php

Разработанная модель может стать основой для серии сценарных модификаций с использованием разных вариантов функций регуляции:

dN

= rN

dt

3 (N − L) − F (N (t − τ ); J ). (13)

(

√

( N )Θ)

−

1

K

Подобные варианты интересны для моделирования замедленного иммунного ответа.

Об авторах

А. Ю. Переварюха

Автор, ответственный за переписку.
Email: madelf@pisem.net
ORCID iD: 0000-0002-1049-0096

кандидат технических наук, руководитель группы Биоинформатики и математической геофизики

Список литературы

Дополнительные файлы