Оптимальное управление однородной популяцией, находящейся под действием естественного отбора

Полный текст

Оптимизация промысла, являясь ключевой задачей математической биологии, не теряет актуальности. Под оптимальной понимается такая доля изъятия из популяции, которая обеспечивает стабильный максимальный равновесный уровень промыслового изъятия при условии невырождения популяции [1, 3], что в теории позволяет балансировать между экономической выгодой и сохранением биоресурсов. Актуальность проблемы возрастает  из-за того, что антропогенное воздействие затрагивает жизнедеятельность большинства природных популяций и ставит под угрозу их существование; здесь можно отметить как видоизменение и фрагментацию среды обитания, так и переэксплуатацию отдельных популяций. Это требует уточнения прогнозов динамики популяций [5, 7]. Современные методы включают машинное обучение и анализ больших данных для управления промыслом [например, 8-9], однако такие методы более ориентированы на внешние проявления изменения динамики численности популяций и сообществ, но не стремятся выявить внутрипопуляционные процессы, являющиеся причиной этих изменений.

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

Рассмотрено влияние промысла с постоянной долей изъятия (u) на динамику численности и генетического состава популяции с неперекрывающимися поколениями, подверженной действию плотностно-зависимого естественного отбора:

 

где n – это номер поколения,   - средняя приспособленность популяции в n-м поколении, x_n – численность n-ого поколения, измеряемая в относительных единицах; q_n – частота аллеля A в n-м поколении, w_AA(n), w_Aa(n), w_aa(n) – приспособленности генотипов AA, Aa и aa соответственно в n-м поколении; приспособленность зависит от численности популяции линейно:.  Тогда каждый генотип характеризует его ресурсный (K_ij>0) и мальтузианский (R_ij>0) параметры.

Отметим, что в случае равенства приспособленностей всех генотипов:  модель (1) представляет собой дискретный аналог модели Ферхюльста:.  Корректность данной модели (x ≥ 0) сохраняется при выполнении следующих условий a  [3]. Таким образом, это накладывает на параметры и начальные условия модели (1) следующие ограничения: R_ij < 3 и .

Величина изъятия («урожай») в n-м поколении составляет . В соответствии с концепцией оптимального промысла необходимо определить уровень изъятия из популяции, который обеспечивает стабильный максимальный равновесный промысловый «урожай» при условии невырождения популяции. Система (1) имеет три равновесных решения: два мономорфных ({q=0, x=x₁} и {q = 1, x = x₂}) и одно полиморфное {q = q₃, x = x₃}. Из первого уравнения системы (1) можно выразить долю изъятия u:  и получить выражение, определяющее объем изъятия в виде: .

Для того чтобы промысел с постоянной долей изъятия был оптимальным, необходимо выполнение условия: , следовательно, . Далее были определены значения численности и оптимальные доли изъятия в каждом нетривиальном равновесии системы (1), а также проанализированы условия их устойчивости.

Оказалось, что стационарные значения частот аллелей в состояниях, обеспечивающих максимальный уровень изъятия, остаются такими же, как и в случае, когда промысел отсутствует, а равновесные значения численности уменьшаются в два раза. Кроме того, аналитически показано, что при оптимальном изъятии, соответствующем отдельному равновесию, данная неподвижная точка всегда имеет одно собственное значение, не превосходящее единицы по модулю. Второе определяется взаимным расположением приспособленности гетерозиготы и гомозигот при соответствующем данному равновесию уровне численности x_i. Мономорфная стационарная точка оказывается устойчивой, если в данном промысловом равновесии приспособленность соответствующей гомозиготы больше приспособленности гетерозиготы; а полиморфное равновесие оказывается устойчивым, если приспособленность гетерозиготы больше приспособленностей обеих гомозигот при соответствующем промысловом равновесном значении численности популяции. С этим связано проявление эффекта смены типа генетического равновесия в результате неизбирательного промысла (пример на рис. 1). Так полиморфная популяция в результате промысла может утратить своё генетическое разнообразие, а мономорфная, наоборот, сохранить. 

Отметим, что каждая из найденных долей изъятия u_ij является оптимальной только для соответствующего ей равновесия. Для других неподвижных точек данный уровень изъятия является произвольным. В результате динамика численности популяции стабилизируется на равновесном уровне, если система оказывается в области притяжения этого равновесия. Поскольку для других двух неподвижных точек системы (1) текущий уровень изъятия в общем случае не удовлетворяет условиям оптимальности при попадании в их область притяжения, динамика численности популяции может продолжить флуктуировать, поскольку одно из собственных чисел текущего равновесия системы (1) может превосходить единицу по модулю. Поэтому вариация начальных условий может менять как результат эволюции промысловой популяции, так и её динамический режим.

Об авторах

О. Л. Жданова

Институт автоматики и процессов управления ДВО РАН

Автор, ответственный за переписку.
Email: axanka@iacp.dvo.ru
ORCID iD: 0000-0002-3090-986X
Россия, ул. Радио 5, г. Владивосток, 690041

Е. А. Колбина

Институт автоматики и процессов управления ДВО РАН

Email: pavlova@iacp.dvo.ru
ORCID iD: 0000-0001-9565-5605
Россия

E. Я. Фрисман

Институт комплексного анализа региональных проблем ДВО РАН

Email: frisman@mail.ru
ORCID iD: 0000-0003-1629-2610
Россия, ул. Шолом-Алейхема 4, г. Биробиджан, 679016

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

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