CC BY
26
УДК 629.075
МОДЕЛИРОВАНИЕ ПОРОГОВОГО ЭФФЕКТА В ЭФФЕКТИВНОСТИ ВОСПРОИЗВОДСТВА ВОЛЖСКОЙ ПОПУЛЯЦИИ ОСЕТРА
©2013 А.Ю. Переварюха
Санкт-Петербургский институт информатики и автоматизации РАН, г. Санкт-Петербург
Поступила 20.11.2013
Допустимый баланс между воспроизводством и смертностью от разных факторов может резко нарушаться из-за непостоянства репродуктивного потенциала, что показал ретроспективный анализ сведений о нересте осетра в Волге. Статья рассматривает моделирование популяцион-ного процесса при существовании различных уровней эффективности воспроизводства. Предложена система ОДУ, формализующая влияние скорости роста на смертность в ювенальный период развития и учитывающая возможность потерь неоплод отворенной икры. В вычислительной среде получена дискретно-непрерывная динамическая система, позволяющая оценивать последствия резкого уменьшения эффективности воспроизводства, что сказывается в наличии минимально допустимой для выживания вида Ааретег gueldenstaedtii численности нерестовой части популяции.
Ключевые слова: пороговые эффекты в динамике популяций, деградация волжских осетровых, дискретно-непрерывные динамические системы.
ВВЕДЕНИЕ
Методы математического моделирования по-пуляционных процессов развиваются с 20-х годов прошлого века и при существенном интересе исследователей к данной области накоплен набор общедоступных подходов. Можно выделить в классификации большие группы «вольтерров-ских» моделей, формализующих взаимодействие конкурирующих популяций на основе систем ОДУ, методы трофо-динамического подхода предназначены для описания переноса энергии по пищевым цепям. В ихтиологии и энтомологии развивались модели, ориентированные на определение механизмов саморегуляции численности при ограничивающем влиянии выделенных факторов среды. При наличии арсенала моделей практическая организация рационального использования биоресурсов на основе прогностической математического аппарата оценивается как сложная. Поставивший под угрозу исчезновения пере-лов осетровых рыб Каспийского моря можно объяснить ошибками и развалом системы регулирования эксплуатации, но коллапс рыбных запасов отмечался в странах экологически ответвленным обществом и действующей научно-обоснованной системой регулирования промысла, как в Норвегии [1]. Вынужденными решениями об остановке уже и так оскудевшего промысла удается лишь более или менее своевременно предотвращать дальнейшие необратимые последствия деградации биоресурсов.
Переварюха Андрей Юрьевич, кандидат технических наук, старший научный сотрудник, [email protected]
К причинам отсутствия универсальной отработанной методики можно отнести то, что построение математического аппарата моделей изначально было слабо связано с данными наблюдений. Теоретическая экология формулируется на уровне качественных принципов, а не строгих законов. Заимствованные из механики методы описания процессов в отсутствии строго формализованных законов оказались не столь эффективны в новой области применения. Начиная с уравнения П. Ферхюльста предполагалось, что популяция характеризуется некоторым репродуктивным потенциалом, г-параметром, что выразилось в том числе в представлении о К- и г-стратегиях выживания видов. Однако, наблюдения свидетельствуют, что г-параметр постоянен только в некоторых усредненных условиях [2].
В статье на основе анализа наблюдаемых эффектов предлагается модель формирования поколений, описывающая скорость убыли численности в раннем онтогенезе рыб. Рассматриваются свойства новой дискретно-непрерывно мульти-стабильной динамической системы. Кратко обсуждается сложноинтерпретируемая динамика применяемых в промысловых расчетах функциональных итераций с появлением хаотического режима.
ДИНАМИЧЕСКАЯ СЛОЖНОСТЬ МОДЕЛЕЙ ИХТИОЛОГИИ
Службами учета Канады неоднократно отмечался нестабильный характер воспроизводства облавливаемых популяций лососевых рыб, мигрирующих для размножения в реки. Фиксировалось отсутствие пропорциональной зависимости К ~ кН между нерестовым запасом и сезонным
2287
пополнением. Возникла теория о роли различных факторов смертности в зависимости от плотности рыб на нерестилищах [3]. Личиночные стадии крупных осетровых и лососевых протекают при высокой плотности, даже если взрослых рыб не так много. Канадский биолог У.Е. Рикер, рассматривая причины уменьшения пополнения при возрастании нерестового стада, предложил теоретическую концепцию о существовании зависимости запаса и пополнения, и роли её свойств при долговременной эксплуатации [4]. Рикер предложил простую, как казалось, формулу расчета величины пополнения R = f(S)
f(S) = aS exp(-bS), (1)
где S - от английского «stock» величина нерестового запаса, Ъ - коэффициент, отражающий величину, обратную количеству производителей, при котором число выжившей молоди максимально, имеет смысл Ъ= 1; а - параметр характеристики репродуктивного потенциала. График зависимости пополнения от численности производителей называется кривой пополнения. (1) задает куполообразную кривую с асимптотой lims f{S) = 0 и единственным нетривиальным пересечением с биссектрисой координатного угла R=S.
Очевидно, что при применении модели (1) возникает ряд сложностей. При увеличении количества отложенной икры выживаемость молоди должна уменьшаться до неразличимых значений, что противоречит наблюдениям над аквариумными популяциями. Пополнение уменьшается до некоторого Q>0 = const. (1) при моделировании на ЭВМ сообщества гидробионтов часто использовалась как кусочно-гладкая, до некоторого критического количества Sc. при S > Sc количество молоди предикативно устанавливалось Q = const.
При низких численностях производителей (1) предсказывает увеличение скорости воспроизводства. При стремящемся к нулю количестве отложенной икры, эффективность пополнения стремится к предельному значению: dR
— = аехр(-й5')(1 - bS), lim acxp(-bS)(\ -bS) = а.
dS s^o
(la)
Подобное свойство сделало сомнительным использование (1) для моделирования деградировавших из-за нерационального промысла популяций рыб, в частности интересующей нас динамики численности волжских осетровых. Более того, (1а) не соответствует концепции экологии, известному как принцип Олли (W. Allee) [1].
Модель (1) обладает интересными свойствами, о которых не имел представления ее разработчик
и потому не предложил для них адекватной биологической интерпретации.
Представим процесс изменения состояния популяции, определяемый зависимостью запас-пополнение в виде динамической системы, математического объекта, для которого можно указать набор величин: динамических переменных, характеризующих состояние системы. Значения таких переменных в последующий момент времени рассчитываются из текущих значений по закону, называемому оператором эволюции. Динамическая система представляется тройкой
состоящей из фазового пространства М, времени Т, оператора эволюции v|/. Для Vx е М и Vi, s е Т выполняется условие:
V(y(x,t),s) = \\f(x,s + t).
Множество {\|/(i)(x)}ieJ, фазовая траектория
точки х. Графически эволюция диссипативной динамической системы во времени представляется движением точек в фазовом пространстве к аттрактору, подмножеству фазового пространства А CIМ инвариантное относительно эволюции:
(А) = А для всех / е / и такое, что существует окрестность U множества А. в которой для всех у gU выполняется lim i|/'j(у) = А. Регу-
t—>СО
лярными аттракторами для отображения интервала \|/ : I —> I считаются состояние равновесия с
неподвижной точкой х*: lim '(}') = х* и устой-
t—>СО
чивый цикл, соответствующий режиму периодических автоколебаний.
Стремящиеся к одному аттрактору множество точек есть его область притяжения, которая может быть непрерывным подинтвалом Q d I, так и образовывать фрактальное множество [5]. Рассмотрим динамическую систему как полугруппу
итераций {\|/;)};а0, и пусть Rn., R], ... последовательность точек эволюции системы, определённых условием R. | = \|/(Д ) при всех / > 0.
Качественное поведение зависит от параметра а, он является управляющим. До определённого значения а, не превышающего бифуркационное, система стремится к точечному аттрактору R* = In а/Ъ . Первый метаморфоз поведения системы происходит, когда производная, вычисленная в неподвижной точке, перестаёт удовлетворять критерию устойчивости | v|у '(R )|<1. Для (1) это происходит при выполнении условия
а > е2:
2288
\|/ \К) = ае м - ЬЯае ш,
* - 1п (2 -ь^^ а(\ — 1п а) \\s\R ) = ае ъ -Ъ-ае ъ = —-- = 1 - 1па.
Ъ еЬа
1 -\па = —1, а = е
Параметр Ь не влияет на топологические характеристики фазового портрета. Динамическая система стремится в устойчивое циклическое состояние с периодом 2: глобальный аттрактор, состоящий из двух периодических точек. Факт, что
хп+1 = ахпъ ы'п имеет в большом диапазоне параметра цикл, говорит о том, что не стоит от построенного по эмпирическим данным графика ожидать характерной купол о образной кривой.
Далее увеличивая а, будем наблюдать рост амплитуды колебаний, и по достижению порогового значения а> 12.51, произойдёт бифуркация удвоения периода и установится цикл р = 4. При дальнейшем увеличении а все быстрее реализуется каскад бифзфкаций удвоения периода. Около а> 14,8 невозможно выделить устойчивых точек или замкнутого цикла. Происходит детерминированный хаос, напоминающий стохастический процесс очень чувствительный к начальным условиям. Траектория притягивается к подмножеству фазового пространства, получившему из-за своих фрактальных свойств название «странный аттрактор». Важно знать, что сценарием хаотиза-ции через бесконечный каскад удвоений, исследов аным на примере хи+1 = г\хп (1 - хп ) М. Фэйгенбаумом в [6], спектр нелинейных эффектов не ограничивается. В ограниченных диапазонах при а> 14,8 наблюдаются «окна периодичности» устойчивых циклов р Ф 2г, далее вновь испытывающих удвоения.
В хп+1 = ах„е
-Ьх„
времени почти регулярных колебаний (ламинарная фаза) с интервалами хаотического поведения (турбулентная фаза). Перед появлением циклических точек ветви графика третьей итерации функции /3(х) образуют так называемые каналы, в строгой последовательности по которым траектория движется довольно долго, что соответствует ламинарной фазе перемежаемости. Выход точки траектории из канала определяет короткую турбулентную фазу - отчетливо видимый на графике «всплеск». Такое поведение может служить индикатором скорых резких изменений.
Рис. 1. Перемежаемость с циклом р = 3 в
наблюдаем последовательность касательных бифуркаций с появлением устойчивых циклов периодов р Ф 21. Ли и Йорк в 1975 г. опубликовали статью «Period three implies chaos», где показали, что если одномерное отображение R +l = R. ) при некотором значении
параметра а — Ь имеет цикл периода р = 3, то
оно при i имеет и бесконечное множество циклов всех других периодов. Ранее А. Шарковский в [7] доказал, если одномерное отображение f : I —> I имеет цикл р = п, то оно также имеет
циклы со всевозможными периодами предшествующими числу п, среди целых чисел, выписанных в порядке Шарковского, который завершает число три.
Появлению окон предшествует явление «перемежаемости» (отмечен стрелкой на рис. 1). Перемежаемостью считается режим чередования во
-Ьх„
ПРОБЛЕМА ИНТЕРПРЕТАЦИИ НЕЛИНЕЙНОСТИ
Появление хаоса интерпретировалось в биологическом контексте как последствие повышения репродуктивного потенциала популяции. В (1) это а, но для динамики подобных унимодальных
отображений при мы увидим набор резких
изменений поведения, для которых нам трудно найти рациональное толкование. При использовании более сложных нелинейных дискретных моделей практические выводы будут определяться свойствами критерия устойчивости аттрактора как функции нескольких переменных. Рассматривая различные воздействия на популяцию, подчиняющуюся (1), мы получим, например, что искусственное воспроизводство хаотизирует динамику, а промысел наоборот уменьшает амплитуду колебаний. Однако, выводы будут универсальны для целого семейства итераций функций, удовлетворяющих критериям теоремы Д. Синджера, называемых ^^/-отображения.
Для построения куполообразных кривых запас-пополнение были предложены довольно сложные преобразования исходных данных наблюдений. Возникает вопрос: имеют ли смысл методы построения кривой, если эмпирические данные о воспроизводстве не прошедшей все стадии деградации популяций и для которых справедлива (1), будут представлены в виде сгущений точек на графике, а не ожидаемой куполообразной кривой с явным максимумом.
Не для всех рыб актуальна теория пополнения, в гораздо большей степени ее методы обоснованы
2289
для анадромных видов, заходящих из морей для нереста на речном грунте.
ПРОБЛЕМА РЕЗКОГО СОКРАЩЕНИЯ ВОСПРОИЗВОДСТВА
Изучение ретроспективных данных о популяциях анадромных рыб Нижней Волги позволило выявить нелинейные особенности в формировании пополнения. Общей задачей наших исследований является анализ на основе вычислительных модельных сценариев антропогенно-
обусловленных процессов деградации биоресурсов экосистемы Волго-Каспия. К прогнозированию уловов осетровых методы динамической теории пополнения У. Рикера ранее не применялись, основным фактором считался гидрологический режим реки.
Среди разнообразной опубликованной информации интересны сведения о количестве молоди осетра, полученной от некоторого запаса, обработанные нами методом скользящей средней (рис. 2). Количество пропущенных на нерест в Волгу производителей за период наблюдений изменялась существенно: от 2500 тыс. экз. в 1979-81 гг. до 25 тыс. в 2000-х гг. [8]. В настоящее время мы не имеем достоверной информации о нынешнем состоянии популяции Ааретег gueldenstaedtii и численности нерестовой части. С 2010 г. волжские популяции осетровых внесены в «Красную книгу» и промысел запрещен. Имеющиеся сведения с 1968 г. об общей численности осетровых Каспия по данным тралового учета [9] многими специалистами оцениваются как сомнительные [10] и менее достоверные, чем информация о проходе производителей на нерестилища в створе реки.
^оллкено праизэцдителЕй- "осярв" тыс аи
Рис. 2. Данные о воспроизводстве осетра
Ааретег guldenstaedtii.
Кривая воспроизводства осетра (на рис. 2 среднее по трем ближайшим) характеризуется как изломанная куполообразная. Мы можем восстановить отчетливо гладкую кривую, а не сгущение точек. Следовательно, популяционная динамика определялась промысловой смертностью более,
чем зависимостью /(5). Стрелкой отмечен момент: при пропуске менее 350-10Д шт. непропорционально резко уменьшается количество скатывающейся в море молоди. Выявлен статистически достоверный пороговый эффект скачкообразного перехода эффективности воспроизводства на некий отличный от нуля минимальный уровень.
Подобное явление оказалось характерно не только для единственной волжской популяции. Сам У. Рикер не рассматривал (1) как единственно возможную нелинейную форму зависимости между запасом и пополнением, описав в работе [11] неожиданный для промысла эффект резкого падения воспроизводства, нарисовав для его объяснения кривую воспроизводства как подобие ступенчатой функции. Рикер назвал эффект «лезвием ножа», т.к. даже небольшой перелов 5% может подорвать запасы такой популяции и потому оптимальная стратегия промысла, изымающая весь излишний прирост нерестового стада, крайне опасна для практики в плане стабильной долговременной эксплуатации.
Зоологами неоднократно подтверждалось для не заботящихся о потомстве видов, что существует оптимальный для воспроизводства диапазон численности репродуктивной части популяции. При уменьшении численности эффективность воспроизводства рыб должна снижаться, так как уменьшается вероятность встречи особей разного пола, пригодных для размножения участках рек, но ожидаемое снижение может оказаться непредвиденно резким.
Модельное описание амплитуды изменения репродуктивного потенциала, с сохранением гладкого характера зависимости, целесообразно провести с применением современного подхода к вычислительному исследованию гибридных динамических систем.
РАЗРАБОТКА НОВОЙ РЕПРОДУКТИВНОЙ МОДЕЛИ
Наличие ограниченных пищевых ресурсов в неявном виде учитывалось известными моделями, но общим недостатком функциональных итераций нужно считать игнорирование проявления декомпенсационного фактора смертности. Не учитывался фактор изменения пищевых потребностей по мере развития молоди, которое происходит синхронно с началом сезонного уменьшения кормовой базы. Декомпенсационные факторы - увеличивающие смертность при уменьшении плотности особей снижают эффективность нереста, уменьшая количество икры реально вступившей в репродуктивный процесс при $ = К . Прямой учет такой зависимости становится очевидной необходимостью при моделировании динамики популяций подвергающихся перелову. Задачей разработки новой модели запас-
2290
пополнение стало создание гибкого математического аппарата для согласования качественного характера поведения имитационной модели промысловой популяции с ретроспективными данными о численности молоди волжских анадром-ных рыб в различные годы и объеме вылова.
Очевидно, что поведение модели зависит от математических особенностей выбранной функции воспроизводства, особенно вне диапазона оптимальной численности нерестового стада. Важны не только наличие минимума и максимума зависимости, но и точки перегиба.
Формализуем увеличение пищевых потребностей молоди исходя из факта, что скорость роста находится в обратной зависимости от численности поколения, но не в обратно пропорциональной. Это согласуется с данными наблюдений биологов в частности, с результатами экспериментов над ростом рыб при различной их плотности. Согласно наблюдениям, при увеличении плотности возникает асимметричное распределение размерной структуры популяции в сторону преобладания особей с меньшими размерами. Описывать убыль численности поколения N (!) на интервале модельного времени [0,7] будут следующие объединённые в систему дифференциальные уравнения (а, ß, с - константы):
— = -(aw(t)N(t) + ®(S)$)N(t) dt
^¡r = , / , w(0) = wn, Щ0) = IS.
dt ^(0+C
где: S — величина нерестового запаса; w(t) -отражает уровень размерного развития поколения, влияющий на увеличение пищевых потребностей; g - параметр, учитывающий ограниченность количества доступных для молоди кормовых объектов. С, — параметр ограничивающий темп развития вне зависимости от N(1); X -средняя плодовитость особей; а - коэффициент компенсационной смертности; ß - коэффициент декомпенсационной смертности; k е (1,2); интервал уязвимости t G [0, / ] специфичен для вида.
Убывающая функция
ОД = 1 + ехр(-с£2), 0(Я) -И отра-
жает действие известного в экологии эффекта Олли (W. Allee) или более точнее «эффекта агрегированной группы». Отмечалось неоднократно, что для животных ведущих стадный образ жизни выживаемость молоди в группе повышается, хотя такое существование увеличивает пищевую конкуренцию у взрослых. Когда численность нерестового стада низка, то большая протяженность русловых нерестилищ становится уже негативным фактором.
Функция эффекта Олли 0(Л) введена в уравнение для N(1), а не в переопределение начальных условий N(0) = т.к. даже не вступившие в репродуктивный процесс кладки икры после их гибели ухудшает токсикологическую обстановку на нерестилище, с -характеризует степень действия эффекта, который не проявляется если численность взрослой части популяции осетровых достаточно велика. При незначительном количестве повторно нерестующих рыб численно получаем Я = N(7).
Графиком исследованной с применением численного решения в инструментальной среде последовательности задачи Коши (2) для вычисления численности выживших особей поколения к моменту времени Т новой модели является унимодальная кривая с уменьшающимся наклоном ниспадающей правой ветви. Кривая имеет ненулевую горизонтальную асимптоту, две нетривиальные точки пересечения с биссектрисой координатного угла Я =5" и характеризуется при двух точках перегиба незнакопостоянным дифференциальным инвариантом Шварца (рис. 3). Фазовое пространство в случае динамической системы
\|/" (Ни) на основе системы (2) разделяется на две
области притяжения: О] и ГК Границей между ними является неустойчивая стационарная точка
7?^ервого пересечения кривой с биссектрисой
тл*
Второе равновесие л, устойчиво.
Анализ устойчивости неподвижных точек дискретно-непрерывной системы в инструментальной среде можно проводить с использованием свойства второй итерации \|/(\|/(х)). Необходимым и достаточным условием устойчивости неподвижной точки одномерного отображения
х*является неравенство: \|/2(х) > X при х<х*, \|Г(х) < Хприх>х*.
О Z ООО 4 ООО 6 ООО 0 ООО 10 ООО 12 ООО 14 ООО 16 ООО
Нерестовый запас .V. шт.
Рис. 3. Зависимость = /(Б) полученная из системы (2).
(2) исследовалась с применением отечественной инструментальной среды моделирования AnyLogic-6, предоставляющей выбор численных методов с изменяющимся шагом интегрирования
2291
доя работы с системами ОДУ первого порядка записанными в форме Коши.
ЗАКЛЮЧЕНИЕ
Поведение траектории динамической системы, с использованием в качестве оператора эволюции решение задачи Коши (2), качественно отличается от итераций (1) возможностью притяжения к двум аттракторам. Соответственно присутствуют две области притяжения, границей которых служит репеллер, неустойчивая особая точка первого пересечения кривой с биссектрисой координатного угла. Траектории с начальными условиями, разделёнными репеллером, покидают его окрестность и приближаются к разным аттракторам. Один из аттракторов точка (0,0) на плоскости 1{/Н. Если численность популяции соответствует области притяжения этого аттрактора, произойдёт вымирание. Данное свойство отражает наличие минимально допустимой численности для существования вида.
Для промысла кривая подобного вида на рис. 3 означает, что популяция долго способна поддерживать высокий процент изъятия, но потом произойдет резкое снижение уловов, которое очень сложно предвидеть экспертам. Деградация может растянуться на длительный период, когда численность популяции будет колебаться около неустойчивого равновесия. В подобной ситуации главное вовремя оценить, что действительно происходит с запасами и своевременно ввести мораторий на промысел, не уповая на хороший уровень половодья.
Работа выполнена при поддержке РФФИ (Проект 13-07-00925).
СПИСОК ЛИТЕРАТУРЫ
1. Гиске Я., Шулдаль Х.Р., Шлагстаг Д. Экологиче-
ское моделирование для рыболовства // Модели многовидового управления [под ред. Т. Редсета]. М.: Изд-во ВНИРО, 2002. С. 17-79.
2. ОдумЮ. Экология. М.: Мир, 1986. 376 с.
3. Neave F. Principles affecting the size of pink and chum salmon population in British Columbia // Journal of the fisheries research board of Canada, 1953. Vol. 9. №9. P. 450-491.
4. Richer W. Stock and recruitment // J. Fisheries re-
search board of Canada, 1954. №.11. C. 559-623.
5. Переварюха А.Ю. Об определении фрактальных объектов в динамике моделей управления биоресурсами // Труды СПИИРАН, 2013. Т.24. С. 211221
6. Feigenbaum М. Universal behavior in nonlinear sys-
tems // Physica D., 1983. Vol. 7. №1-3. P. 16-39.
7. Шарковский A.H. Сосуществование циклов у непрерывного преобразования прямой в себя // Украинский математический журнал. 1964. Т.26. №1. С. 61-71.
8. Вещее П.В., Гутенева Г.И. Эффективность есте-
ственного воспроизводства осетровых в Низовьях Волги в современных условиях // Экология. 2012. №2. С. 123-127.
9. Панин Г.Н., Мамедое P.M., Митрофанов И.В. Со-
временное состояние Каспийского моря. М.: Наука, 2005. С. 281.
10. Легеза М.И., Маилян Р.А. К методике исследования численности осетровых в Каспийском море // Тезисы отчетной сессии ЦНИОРХ. Астрахань, 1974. С. 82-83.
11. Richer W.E. Big effect from small causes: two examples from fish population dynamics // Journal Fisheries research board of Canada. 1963. Vol. 20. №2. P. 257-264.
SIMULATION OF THE THRESHOLD EFFECT IN EFFICIENCY OF REPRODUCTION OF RUSSIAN STURGEON IN THE VOLGA RIVER
© 2013 A.U. Perevaryukha
St. Petersburg Institute for Informatics and Automation, St. Petersburg
Acceptable balance between reproduction and mortality from various factors could dramatically disturbed because of the impermanence in reproductive potential, that showed a retrospective analysis of data on spawning sturgeon in the Volga. Article considers development of the model of population process with the existence of different levels of reproductive efficiency. We proposed a system of ODE, which formalizes the influence of the growth rate on mortality in juvenile period and take into account the loss of unfertilized eggs. In a computing environment was received discrete-continuous dynamical system which allowed assessing the effects of a sharp decrease in reproductive efficiency. This effect leads to the presence minimum acceptable size for the survival of Acipenser gueldenstaed-tii spawning population.
Key words, threshold effects in the dynamics of populations, degradation of the Volga sturgeon, discrete-continuous dynamical systems
Perevaryukha Andrei Yrievich, candidate of technical sciences, senior Researcher, [email protected]
2292
