Все выпуски

Рассматривается классическая для нелинейной динамики модель «брюсселятор», дополненная третьей переменной, играющей роль быстро диффундирующего ингибитора. Модель исследуется в одномерном случае в области параметров, где проявляются два типа диффузионной неустойчивости однородного стационарного состояния системы: волновая неустойчивость, приводящая к самопроизвольному формированию автоволн, и неустойчивость Тьюринга, приводящая к самопроизвольному формированию стационарных диссипативных структур, или структур Тьюринга. Показано, что благодаря субкритическому характеру бифуркации Тьюринга взаимодействие двух неустойчивостей в данной системе приводит к самопроизвольному формированию стационарных диссипативных структур еще до прохождения бифуркации Тьюринга. В ответ на различные случайные шумовые возмущения пространственно-однородного стационарного состояния в исследуемой параметрической области в окрестности точки двойной бифуркации в системе могут устанавливаться различные режимы: как чистые, состоящие только из стационарных или только автоволновых диссипативных структур, так и смешанные, при которых разные режимы проявляются в разных участках расчетного пространства. В рассматриваемой параметрической области система является мультистабильной и проявляет высокую чувствительность к начальным шумовым условиям, что приводит к размытию границ между качественно разными режимами. При этом даже в зоне доминирования смешанных режимов с преобладанием структур Тьюринга значительную вероятность имеет установление чистого автоволнового режима. В случае установившихся смешанных режимов достаточно сильное локальное возмущение в участке расчетного пространства, где проявляется автоволновой режим, может инициировать локальное формирование новых стационарных диссипативных структур. Локальное возмущение стационарного однородного состояния в исследуемой области параметрического пространства приводит к качественно схожей карте устоявшихся режимов, при этом зона доминирования чистых автоволновых режимов расширяется с увеличением амплитуды локального возмущения. В двумерном случае в системе не устанавливаются смешанные режимы. При эволюции системы в случае появления локальных структур Тьюринга под воздействием автоволнового режима со временем они заполняют все расчетное пространство.

Предложена трехкомпонентная модельпланк тонного сообщества с дискретным временем. Сообщество представлено зоопланктоном и двумя конкурирующими за ресурсы видами фитопланктона: токсичным и нетоксичным. Модельдв ух связанных уравнений Рикера, ориентированная на описание динамики конкурентного сообщества, используется для описания взаимодействия двух видов фитопланктона и позволяет неявно учитывать ограничение роста биомассы каждого из видов-конкурентов доступностью внешних ресурсов. Изъятие фитопланктона за счет питания зоопланктоном описывается трофической функцией Холлинга II типа с учетом насыщения хищника. Способность фитопланктона защищаться от хищничества и избирательность питания хищника учтены в виде ограничения потребления: зоопланктон питается только нетоксичным фитопланктоном.

Анализ сценариев перехода от стационарной динамики к колебаниям численности сообщества показал, что потеря устойчивости нетривиального равновесия, соответствующего сосуществованию двух видов фитопланктона и зоопланктона, может происходитьч ерез каскад бифуркаций удвоения периода, также возникает бифуркация Неймарка – Сакера, ведущая к возникновению квазипериодических колебаний. Вариация внутрипопуляционных параметров фито- или зоопланктона может приводитьк выраженным изменениям динамического режима в сообществе: резким переходам от регулярной к квазипериодической динамике и далее к точным циклам с небольшим периодом или даже стационарной динамике. В областях мультистабильности возможна кардинальная смена как динамического режима, так и состава сообщества за счет изменения начальных условий или же текущего состава сообщества. Предложенная в данной работе трехкомпонентная модель динамики сообщества с дискретным временем, являясь достаточно простой, позволяет получитьадекв атную динамику взаимодействующих видов: возникают динамические режимы, отражающие основные свойства экспериментальной динамики. Так, наблюдается динамика характерная для модели «хищник–жертва» без учета эволюции — с отставанием динамики хищника от жертвы примерно на четвертьперио да. Рассмотрение генетической неоднородности фитопланктона, даже в случае выделения всего двух генетически различных форм: токсичного и нетоксичного, позволяет наблюдатьв модели как длиннопериодические противофазные циклы хищника и жертвы, так и скрытые циклы, при которых плотностьч исленности жертв остается практически постоянной, а плотность численности хищников колеблется, демонстрируя влияние быстрой эволюции, маскирующей трофическое взаимодействие видов.

Даже беглый взгляд на впечатляющее множество современных работ по математическому моделированию популяционной динамики позволяет заключить, что основной интерес авторов сосредоточен вокруг двух-трех ключевых направлений исследований, связанных с описанием и анализом динамики, либо отдельных структурированных популяций, либо систем однородных популяций, взаимодействующих между собой в экологическом сообществе или (и) в физическом пространстве. В рамках данной работы приводится обзор и систематизируются научные исследования и результаты, полученные на сегодняшний день в ходе развития идей и подходов математического моделирования динамики структурированных и взаимодействующих популяций. В вопросах моделирования динамики численности изолированных популяций описана эволюция научных идей по пути усложнения моделей — от классической модели Мальтуса до современных моделей, учитывающих множество факторов, влияющих на популяционную динамику. В частности, рассматриваются динамические эффекты, к которым приводит учет экологической емкости среды, плотностно-зависимая регуляция, эффект Олли, усложнение возрастной и стадийной структуры. Особое внимание уделяется вопросам мультистабильности популяционной динамики. Кроме того, представлены исследования, в которых анализируется влияние промыслового изъятия на динамику структурированных популяций и возникновение эффекта гидры. Отдельно рассмотрены вопросы возникновения и развития пространственных диссипативных структур в пространственно разобщенных популяциях и сообществах, связанных миграциями. Здесь особое внимание уделяется вопросам частотной и фазовой мультистабильности популяционной динамики, а также возникновению пространственных кластеров. В ходе систематизации и обзора задач, посвященных моделированию динамики взаимодействующих популяций, основное внимание уделяется сообществу «хищник–жертва». Представлены ключевые идеологические подходы, применяемые в современной математической биологии при моделировании систем типа «хищник–жертва», в том числе с учетом структуры сообщества и промыслового изъятия. Кратко освещены вопросы возникновения и сохранения мозаичной структуры в пространственно распределенных и миграционно связанных сообществах.

Проводится исследование вольтерровской модели, описывающей конкуренцию трех видов. Соответствующая система дифференциальных уравнений первого порядка с квадратичной правой частью после замены переменных сводится к системе с восемью параметрами. Два из них характеризуют скорости роста популяций, для первого вида этот параметр принят равным единице. Остальные шесть коэффициентов задают матрицу взаимодействий видов. Ранее при аналитическом исследовании так называемых симметричной модели [May, Leonard, 1975] и асимметричной модели [Chi, Wu, Hsu, 1998] с коэффициентами роста, равными единице, были установлены соотношения на коэффициенты взаимодействия, при которых система имеет однопараметрическое семейство предельных циклов. В данной работе проведено численно-аналитическое исследование полной системы на основе косимметричного подхода, позволившего определить соотношения на параметры, которым отвечают семейства равновесий. Получены различные варианты однопараметрических семейств и показано, что они могут состоять как из устойчивых, так и из неустойчивых равновесий. В случае матрицы взаимодействий с единичными коэффициентами найдены мультикосимметрия системы и двухпараметрическое семейство равновесий, существующее при любых коэффициентах роста. Для различных коэффициентов взаимодействия найдены значения параметров роста, при которых реализуются периодические режимы. Их принадлежность семейству предельных циклов подтверждена расчетом мультипликаторов. В широком диапазоне значений, нарушающих соотношения, при которых обеспечивается существование циклов, получается типичное при разрушении косимметрии медленное колебательное установление. Приведены примеры, когда фиксированному значению одного параметра роста отвечают два значения другого параметра, так что существуют разные семейства периодических режимов. Таким образом, установлена вариативность сценариев развития трехвидовой системы.

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

Рассматривается пространственно-временная модель тритрофической системы, описывающая взаимодействие жертвы, хищника и суперхищника в среде с неоднородным распределением ресурса. Учитываются всеядность суперхищника (Intraguild Predation, IGP), диффузия и направленная миграция (таксис), который моделируется с помощью логарифмической функции от ресурса и плотности жертвы. Основное внимание уделено анализу мультистабильности системы и роли косимметрии в формировании континуальных семейств стационарных решений. С использованием численно-аналитического подхода изучаются пространственно-однородные и неоднородные стационарные решения. Установлено, что при выполнении дополнительных соотношений между параметрами, характеризующими локальное взаимодействие хищников, и коэффициентами диффузии система обладает косимметрией, что приводит к возникновению семейства устойчивых стационарных решений, пропорциональных функции ресурса. Показано, что косимметрия не зависит от функции ресурса в случае неоднородной среды. Проведено исследование устойчивости стационарных распределений с помощью спектрального метода. Нарушение условий косимметрии приводит к разрушению семейства и появлению изолированных стационарных состояний, а также к длительным переходным процессам, отражающим память системы об исчезнувшем семействе. В зависимости от начальных условий и параметров в системе реализуются переходы к режимам с одним хищником (выживание хищника или суперхищника) или к сосуществованию хищников. Численные эксперименты на основе метода прямых (разностная схема по пространственной переменной и метод Рунге – Кутты для интегрирования по времени) подтверждают мультистабильность системы и иллюстрируют исчезновение семейства решений при разрушении косимметрии.

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