CC BY
53
Владикавказский математический журнал 2017, Том 19, Выпуск 4, С. 58-69
УДК 681.3.06
СЦЕНАРИИ НЕВЫНУЖДЕННОИ ДЕСТРУКЦИИ ПОПУЛЯЦИИ В МОДИФИКАЦИИ УРАВНЕНИЯ ХАТЧИНСОНА1
А. Ю. Переварюха
Рассматривается проблема моделирования резких изменений в режиме автоколебаний, присущих видам, которые способны воздействовать на среду своего обитания. Актуальность работы обусловлена необходимостью совершенствования методов математической биологии для все чаще проявляющихся нестационарных и экстремальных типов популяционной динамики. Стремительные переходы к резким флуктуациям численности возникают при инвазиях активно размножающихся видов вредителей. Предложена модификация уравнения Хатчинсона с учетом существенной роли достижения предпороговой численности, меньшей предельной емкости экологической ниши К из уравнения Ферхюльста, и существенно большей нижней пороговой численности Ь из уравнения Базыкина: Ь -С Н < К. В нашем уравнении при изменении действующего запаздывания регуляции т описывается атипичный сценарий развития опасной вспышки насекомых. Как следует из экологических примеров, популяционные циклы с большой амплитудой часто оказываются неустойчивы. Часто цикл — переходный режим. Не всегда происходит плавное затухание осцилляций Ж* (г, ^ К. В новой модели после бифуркации Андронова — Хопфа при Т = т* + £ и появления автоколебаний негармонической формы с увеличением их амплитуды резко происходит потеря диссипативного свойства траектории. Вычислительный сценарий с внезапным выходом неустановившегося цикла Ж* (Тг, £) из области допустимых значений численности интерпретируется как специфическое нарушение функционирования среды, ведущее к деструкции биосистемы в очаге вспышки насекомых или безвозвратной гибели в случае островной популяции млекопитающих.
Ключевые слова: уравнение Хатчинсона, динамика насекомых вредителей, модель особой вспышки численности, бифуркация Андронова — Хопфа, неустойчивый цикл.
Варианты моделей асимптотической динамики популяций для стационарного или периодически возмущенного репродуктивного потенциала развиты достаточно подробно, но существует проблема анализа ряда экстремальных плохо прогнозируемых явлений в экологических процессах, которые входят в противоречие с базовыми принципами регулирования численности.
Ограниченность пределов роста численности популяции — один из таких ключевых принципов для системной экологии [1]. Первой популяционной моделью ограниченного роста стала модель Ферхюльста 1838 г. Асимптотика решения уравнения подразумевала плавный выход численности к стационарному состоянию К, соответствующему балансу со средой:
© 2017 Переварюха А. Ю.
1 Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований, проект № 17-07-00125.
1. Введение
(1)
щг) =
1 + N(0)(ен — 1 )/К
Ферхюльст проводил оценку репродуктивного параметра г уравнения (1) по демографическим данным европейских стран и уже тогда предлагал модификации. Считающаяся хрестоматийной квадратичная саморегуляция численности в (1) не рассматривалась Ферхюльстом как единственно возможная для модели ограниченного роста. В 1845 г. он выписал такое обобщенное уравнение:
и г Ф //
М — = тр — пр оЬ
г К К г
появились значительно позднее. Позднее он предположил [2], что динамика народонаселения Европы развивается согласно (1) только до некоторого порога, так как не наблюдалось замедления темпов прироста населения, которое для (1) должно происходить после N > Щ.
В правую часть (1) можно добавить сомножитель — Ь), увеличив число ее нулей или раздельно учесть смертность, зависящую от плотности скопления §N2 и естественную яN, много увеличив число параметров:
сЙ 7 + сгЖ ' v ;
Ь
для видов с минимально допустимой численностью группы [3]. Отметим, что для формы правой части (1) /(^— Ь) в вычислительных моделях еще необходима компенсирующая сверхбыстрый прирост перенормировка репродуктивного параметра: Г1 ^ г, но вымирание N(Ь) ^ 0 выглядит слишком стремительным и безвозвратным явлением.
Многие популяции могут брать начало от нескольких сохранивших репродуктивную активность особей и через некоторое время демонстрировать режимы сложных флукту-аций. С вселения единичных особей начались регулярные нашествия европейского кукурузного мотылька ОзЬпта пиЫШгз в Северной Америке и затянувшаяся с 2012 г. вспышка дальневосточной самшитовой огневки СуйаИта регзресЬаЫз в Краснодарском крае. Исчезновение ряда популяций (характерно для островных сообществ, где нет разреженной границы ареала) может происходить как раз не из состояния близкого к неустой-
ЬК получения в модели режима, соответствующего частному, но значимому сценарию попу-ляционной динамики. Например, из-за истощения пригодной растительности следующее крайне многочисленное поколение вредителей вдруг оказывается последним, и так часто заканчивается многолетняя вспышка численности.
2. Проблема популяционных автоколебаний
Решения (1) и (2) плохо согласовывались с наблюдавшимися разнообразного вида флуктуациями популяций самых разных животных. Флуктуации регистрируются даже для слонов в Национальном парке Крюгера, потому приходится проводить отстрелы с целью удержать популяцию в оптимальном по расчетам для сохранения растительной среды состоянии [4]. Слоны способны стремительно достигать пиковых значений N > 9000, тогда так для них рассчитано К = 7000. Хотя согласно экологической тео-
К
стремящейся к асимптотическому равновесию N ^ К±е. Признаков уравновешивания у расчетного предела К не наблюдается и после изъятия особей. Через некоторое время после отстрела слоны начинают вновь оказывать чрезмерное давление на растительность. Лемминги считаются примером г-стратегии с нерегулярностью своих больших пиков, но вспышки грызунов для арктических экосистем с редкой растительностью из полярных кустарничков тундры эквиваленты воздействию слонов на леса. Колебательная динамика популяции может оказываться неустойчивой. Для гренландского лемминга давно были известны популяционные циклы значительной амплитуды, которые внезапно завершились после серии пиков в конце XX века, и в настоящее время популяция леммингов мало изменяется на низком уровне численности [5].
J - --г- —- Р- -V---1
1985 1990 1995 2000 2005 2010 Year
Рис. 1. Неустойчивость популяционного цикла лемминга.
Модель взаимодействия двух видов Вольтерра, как хорошо известно, позволяет получить циклическое решение системы двух дифференциальных уравнений. Казалось бы, уравнения позволяют формализовать колебания с отставанием по фазе, наблюдаемые в сообществе «зайцы рыси». Как было установлено позднее, цикличность численности канадских зайцев обнаружена и в регионах, где их враги не обитают. Более того, аналогичные квазигармонические колебания свойственны многим другим северным популяциям грызунов, в том числе обособленным островным группам, и которые по экологическим причинам не могут быть описаны исходя из модели взаимодействия двух видов типа «хищник жертва». Не подтвердились выводы об устойчивых колебаниях из модели Вольтерра в лабораторных экспериментах Г. Гаузе с простейшими организмами. Колебания в лабораторных экспериментах с насекомыми получил С. Утида, но для взаимоотношений вида «паразит хозяин». Паразитизм имеет отличия от трофического взаимодействия «хищник жертва», так как текущая динамика численности паразитической осы зависит от количества доступных яиц (но не взрослых жуков) в предшествующем поколении.
Из лабораторных экспериментов энтомолога Никольсона [6] стало понятно, что колебания численности могут появляться у изолированных популяций, обитающих при постоянных условиях и получающих фиксированное количество корма.
3. Развитие уравнения (1) с запаздыванием
Хатчинсоном [7] была предложена модификация (1), исходя из запаздывающего действия регуляции в размножении дафний, что привело к уравнению с отклоняющимся но времени аргументом «logistic DDE»:
SWi-ВД.
Введенное в уравнение запаздывание т служит характеристикой регуляции — использования и восстановления ресурсов или накопления/распада отравляющих продуктов метаболизма.
Возникающие в уравнении режимы поведения траектории были подробно исследованы во многих работах, в том числе отечественных [8]. При малых значениях запаздывания т динамика модели опишет затухающие колебания N(Ь) ^ К. В (3) установлена возможность возникновение бифуркации Андронова Хопфа с появлением устойчивого предельного цикла N (Ь,г). Нарушение устойчивости состояния равновесия зависит от величины гт. Дальнейшее увеличение гт > вызывает переход в режим релаксационных колебаний. Быстрое возрастание амплитуды колебаний выраженной негармонической формы при увеличении временного промежутка между максимумами и минимумами, стремящимся к неотличимым от нуля значениям, приводят такой релаксационный цикл (рис. 2) к затруднениям при экологическом обосновании.
Рис. 2. Релаксационный цикл в уравнении (3), K = 15000.
Многим насекомым свойственны продолжительные вспышки численности, но они начинаются после преодоления некоторого порогового значения, когда паразиты не могут эффективно подавлять размножение. Формы колебаний для автохтонных и инвазионных популяций могут быть различными. Помимо затухающих и гармонических колебаний может возникать и противоположное явление возникновение флуктуации с затяжными пиками численности. Такая ситуация возможна, например, при переходе некоторых популяций к каннибализму и или может быть обусловлена изменением направления естественного отбора в особых условиях существования. Дополнение «репродуктивного» сомножителя запаздыванием rN(t — т)f [N(t),N(t — т)] не несут экологического смысла.
В [9] рассматривалась модификация модели Хатчинсона (с приведенными коэффициентами) :
^L = XN(t)f(N(t-l)), (4)
где f (■) — дифференцируемая функция, разложимая в асимптотический ряд, для которой выполняются условия
те
/(0) = 1, f{x) = -a0 + J2— > Ж 00, а0 > 0. (4')
k=1 Xk
Условиям (4;) соответствует, например, предложенная ранее в [10] для тех же целей функция
/И = (4//)
(1 + cx)
В подобной модификации (известной как «food limited DDE») расширяется число параметров, определяющих характеристики релаксационного цикла неклассической формы, так как при увеличении с в (4") сжимается амплитуда.
Из дополнительных ограничений при преобразовании исходного уравнения (4) к сингулярно возмущенному x = F(x(t — 1),е) и далее к предельному релейному уравнению с запаздыванием x = R(x(t — 1)) предполагавтся А ^ 1. Для (4) с (4;) утверждается существование единственного орбитально экспоненциально устойчивого цикла. По-прежнему непонятно какой именно нелинейный популяционный процесс может помочь предсказать свойство (4), так как из утверждения о свойствах цикла в [9] получено
min N*(t, А) ~ C\ exp (—Аао), C,a0 = const > 0,
А
использовании функции (4") наши вычислительные эксперименты показывают влияние
с
обсуждается в [9], и требуется условие с ^ 1, так так при малых с свойства цикла (4) будут значительно менее интерпретируемы, чем в (3).
Для (3) и (4) установлено свойство диссппативностп траектории:
N(t) < exp(rr). (4*)
4. Модель динамики вспышек листовертки
Некоторым популяциям вредителей свойственно специфическое нерегулярное явление пилообразной вспышки численности «sawtooth oscillations» [11], которую мы можем описать при дополнении модели (4), (4"). После равновесного состояния при численности бабочек, не вызывающей дефолиации, вдруг появляется серия непоследовательных очень многочисленных поколений, на порядок превышающих среднемноголетнюю численность за предшествующую декаду в очаге вспышки. Так внезапно проявляет себя еловая листовертка Choristoneura fumiferana в лесах Северной Америки (на рис. 3 изменение площади погибшего леса с тремя выраженными пиками но данным [12]), где более 30 лет постоянно наблюдалась полная потеря хвои (дефолиация) на значительных площадях вплоть до прекращения вспышки в конце 1980-х гг. Завершение атипичной вспышки может происходить так же внезапно, как и ее развитие.
Рис. 3. Пилообразная вспышка листовертки в масштабах гибели леса.
Мы можем предложить математическое описание такого частного, но весьма неблагоприятного для лесного хозяйства сценария популяционной динамики, увеличив степень нелинейности правой части в (4"):
^ = гдгт (
(5)
Так удается сгладить недостаток у минимумов шшо<г<т* Ж* (Ь,г) ^ е, е ^ 1, слишком низких для возможностей повторного быстрого роста численности при свойственных вредителям больших г, то отойдя от уеловий (4'): /(х) ^ 0 х ^ В вычислительных исследованиях для (5) мы после бифуркации можем получить реализацию сценария пилообразной вспышки (рис. 4). Основной значимый аспект (5) очередной стремительный неконтролируемый рост численности начинается от значения близкого к существовавшему перед бифуркацией равновесию, оптимального баланса со средой. Не хватает у всех аналогичных модификаций уравнений с запаздыванием некоторых триггерных свойств, чтобы рассматривать аспекты быстрого сценария самозавершения вспышек.
Рис. 4. Переход к релаксационному циклу в новой модели (5).
Замечание. Известно, что в итерациях нелинейных функций Я?+1 = ^(Я?), имеющих максимум ф'(Кт) = 0 ,ф"(Ят) = 0, поведение траектории даже более разнообразно. В ихтиологии применяется:
Я?+1 = аЯ? е-ш1, а> 1, 0 < Ь < 1,
где при а > е2 возникают циклы разнообразных периодов, и дискретизация задачи поможет найти оптимальное решение. Однако, помимо периода циклы итераций отличаются взаимным расположением точек. При бифуркациях удвоения периода вокруг теряющей устойчивость циклической точки возникнут две новые, сохраняя симметрию ветвей на бифуркационной диаграмме сценария Фейгенбаума. Короткие циклы грызунов в таких терминах это монотонные перестановки с пиком в конце, которые можно найти в «окнах периодичности» с циклами р = 2 очень узких диапазонах значений управляющего параметра [13]. Существует ряд причин, но которым возникновение хаотического аттрактора в точке накопления каскада бифуркаций удвоения не является хорошим свойством для описания популяционных процессов.
5. Модификация с предкритическим порогом
Из развития уравнения (3) можно получить описание еще одного специфического экологического явления, которое не было рассмотрено в известных работах. В случае вспышек насекомых-фитофагов понятие емкости экологической ниши из (1) утрачивает смысл, а равновесие при малой их численности имеет другую природу. Из анализа дилеммы парка Крюгера предположим существование некоторого нредкритического порогового уровня Н < К. Для генерирующих вспышки популяций нижний критический порог полагается весьма малым Ь ^ Н. Положим, что достижение значения численно-К
обитания. Переход через мягкий порог имеет значение для механизмов контроля внутри-популяционной структуры. Тогда на динамику системы оказывает влияние отклонение [Н — N(Ь — т)], притом величина отклонения может быть как положительной, так и отрицательной. Модифицируем (3) следующим образом:
(6)
Можно считать, что при смене знака отклонения члены правой части (воспроизводства и регуляции) меняются своими функциональными ролями в модели. Для правой части (6) не выполняются условия (4;), указанные для f (ж) при обосновании утверждения о единственном орбитально устойчивом цикле, потому не обязательно теорема (2.1) из [9] будет полностью справедлива в нашем случае.
6. Вычислительный анализ новой модели
При малом значении запаздывания в (6) получаем затухающие осцилляции с N ^ Н. При увеличении т ми т\ в таком уравнении (с Т\ = 10-3г) возникнет притягивающий цикл (рис. 5, К = 15000, Н = 5000).
Рис. 5. 1 — цикл после бифуркации в (6), 2 — динамика (3) при аналогичных параметрах.
На рис. 5 динамика после плавного прохождения бифуркационного изменения но сценарию Андронова — Хопфа при переходе Т\ через критическое значение (в вычислительном эксперименте рис. 5 все параметры модели сохраняются), показывающая установление цикла при N(0) = H + е, H = 5000, K = 15000.
Цикл быстро становится релаксационным с возрастанием Т\, из имеющейся, например, в [14] теории о модификациях logistic DDE мы не можем подтвердить его глобальную устойчивость. Амплитуда колебаний в отличие аналогичного случая модели (3) не
возрастает до нереалистичных величин, но остается в разумных для биологии границах. Скошенная форма пиков даже более похожа на реальные данные нонуляционной динамики. В (6) при N(0) < Н переход к установившимся флуктуациям происходит не плавно, но через скачкообразный переходный режим.
При дальнейшем увеличении значения т\т произойдет другое резкое изменение поведения траектории, которая перестанет притягиваться к замкнутому подмножеству фазового пространства. Изменение репродуктивного параметра можно трактовать как следствие ослабления давления паразитов на выживаемость личинок. Если давление ослабнет совсем, то при такой бифуркации траектория системы до установления цикла со своей увеличивающейся огромной амплитудой может быть резко выброшена за пределы допустимых для ее существования значений при N(Ь — т) > К. Такая потеря установившегося режима считается жесткой (катастрофической). На рис. 6 ситуация после бифуркации, когда в результате изменения параметра Т\ (в сравнении с аналогичным изменением т в (3) в момент Ь = 150) траектория в релаксационных колебаниях преодолевает значение К, и далее N(Ь) ^ то реализует сценарий кризиса перенасыщения, появления последнего сверхмногочисленного поколения и деструкцию всей биосистемы.
Рис. 6. 1 — выброс за границы емкости среды в (6), 2 — динамика (3) при аналогичном изменении параметра.
Цикл сжимается и выбрасывается, вычислительный эксперимент заканчивается сообщением программы о переполнении при вычислениях с плавающей точкой. Свойство диссипативности (4*) более не выполняется.
7. Проблема обобщения всех модификаций (3)
Выглядит перспективным объединение идей двух предложенных модификаций, но не все варианты, возникающие в уравнениях, с запаздыванием имеют прикладной экологический смысл. Можно представить обобщенное уравнение в следующем виде:
увеличивая нелинейность в правой части до кубической. Тогда образующая для N(0) < Н релаксационный цикл траектории при уменьшении с начнет достаточно параболически плавно и неограниченно возрастать (рис. 7, со = 0.11, бифуркационное с = 0.1 т = 8).
Рис. 7. Медленный сценарий ухода траектории из цикла в (7).
Подобный путь потерн решением свойства диссинативности видится экологически не обоснованным и сценарий нонуляционной динамики вряд ли реализуется. Мы рассматривали ситуацию изолированной саморегулируемой популяции, потому должно присутствовать ограниченное сверху предельное состояние (цикл или равновесие), либо происходит деструкция биосистемы скачкообразно.
Авторы [10], пытаясь преодолеть отмеченные недостатки, в дальнейшем расширили модель (4), (4") до системы, где параметры А, с становятся изменяющимися с периодом ш функциями (обусловленными сезонностью характеристик среды обитания) и /(Ж(Ь — шш)). Однако в лабораторных экспериментах с насекомыми Никольсона и У гады не вносились никакие вынуждающие воздействия. Для кукурузного мотылька ОяЫта пиЬИаИв на Юге России свойственна динамика с серией резких пиков численности, перемежающаяся с годами длительной депрессии популяции, которую иногда связывают с цикличностью активности магнитного поля Солнца, но физического объяснения эффекту пока не получено. Потому актуально модельное описание именно внезапного возникновения специфических популяционных автоколебаний вредителей.
8. Заключение
Предложена модель возникновения популяционных колебаний релаксационной формы, где вместо дальнейшего увеличения амплитуды происходит выброс траектории в режим за рамками ограниченного прироста численности. Будем считать введенное в уравнение значение Н мягким пороговым состоянием отложенного действия, в отличие от явной нижней критической точки Ь в (2) и жесткого уровня остановки прироста К в (1). Емкость среды не выступает более асимптотическим предельным значением для затухающих осцилляций. В новой модели критическим, но спонтанно преодолеваемым состоянием служит К. N > К на время 6Ь > т подразумевает стремительное катастрофическое разрушение популяцией своей среды обитания. Пример катастрофы успешно размножавшейся искусственно вселенной популяции судьба исчезнувших оленей на острове Беринга, куда в 1882 г. было завезено 15 голов, через 20 лет их численность колебалась 1,5 3 тыс. с учетом промысла. Без регуляции к 1917 г. на острове не осталось оленей, обильные лишайниковые пастбища были полностью истощены. В демографии самый известный случай кризиса перенаселения загадка падения цивилизации острова Пасхи.
Таким образом, мы получили вычислительную модель, описывающую переход численности из окрестности неявного порога в устойчивый режим колебаний. В современном обзоре по применению уравнений с запаздыванием [15] в математической биологии
сценарий с разрушением автоколебаний не рассматривался. Бифуркация, для которой требовалось изменение условий регуляции, приводит систему не к изменению амплитуды/периода цикла или переходу к хаотическому движению, но к скачку за границы экологического баланса в виде резкой потери замкнутого колебательного контура. Полученные при исследовании уравнения характеристики достаточно особых случаев попу-ляционного процесса не могут, по нашему мнению, быть распространены на обобщенный случай динамики распространения вида-вселенца. Модификация модели (6) с медленной потерей переходного циклического режима не получает содержательной популяционной интерпретации.
Литература
1. Jeffers J. An Introduction to Systems Analysis: with Ecological Applications.—London: Edward Arnold, 1978.
2. Bacaer N. A Short History of Mathematical Population Dynamics.—London: Springer-Verlag, 2011. DOI: 10.1007/978-0-85729-115-8.
3. Вазыкин А. Д. Теоретическая и математическая экология: проблема опасных границ и критериев приближения к ним // Математика и моделирование.—Пущино: ОНТИ НЦБИ АН СССР, 1990.— С. 232-237.
4. van Aarde R. Culling and the dynamics of the Kruger National Park African elephant population // Animal Conservation.-1999.-Vol. 2, № 2.-P. 287-294. DOI: 10.1111/j.1469-1795.1999.tb00075.x.
5. Gilg O., Sittler Т., Hanski I. Climate change and cyclic predator-prey population dynamics in the high Arctic // Global Change Biology.-2009.-Vol. "l5, № ll.-P. 2634-2652. DOI: 10.1111/j.l365-2486.2009.01927.x.
6. Nicholson A. An outline of the dynamics of animal populations // Austral. J. Zool.—1954.—Vol. 2, № l.-P. 9-65. DOI: 10.1071/Z09540009.
7. Hutchinson G. An Introduction to Population Ecology.—New Haven: Yale University Press, 1978.
8. Ворздыко В. И. Об исследовании популяционной модели Хатчинсона // Диф. уравнения.—1985.— Т. 21.-С. 316-318.
9. Колесов А. Ю., Мищенко Е. Ф., Розов Н. X. Об одной модификации уравнения Хатчинсона // Журн. вычисл. матем. и мат. физики.—2010.—Т. 50, № 12.—С. 2099-2112.
DOI: 10.1134/S0965542510120031.
10. Gopalsamy К., Kulenovic М., Ladas G. Time lags in a «food-limited» population model // Applicable Analysis.-1988.-Vol. 31, № 3.-P. 225-237. DOI: 10.1080/00036818808839826.
11. Cooke В., NealiS V., Regniere J. Insect Defoliators as Periodic Disturbances in Northern Forest Ecosystems.—Plant Disturbance Ecology: the Process and the Response.—Elsevier: Burlington, 2007.— P. 487-525.
12. Gray D. R. Historical spruce budworm defoliation records adjusted for insecticide protection in New Brunswick // J. Acad. Entomol. Soc.—2007.—Vol. 115, № 1.—P. 1-6.
13. Шарковский A. 11.. Коляда С. Ф., ('пник А. Г., Федоренко В. В. Динамика одномерных отображений.—Киев: Наукова думка, 1989.
14. Ladas G., Qian С. Oscillation and global stability in a delay logistic equation // Dynamics and Stability of Systems.-1994.-Vol. 9, № 2.-P. 153-162. DOI: 10.1080/02681119408806174."
15. Ruan S. Delay differential equations in single species dynamics // Delay Differential Equations and Appl.-Berlin: Springer, 2006.-P. 477-517. DOI: 10.1007/1-4020-3647-7^11.
Статья поступила 26 июня 2016 г.
Переварюха Андрей Юрьевич Санкт-Петербургский институт информатики и автоматизации РАН,
старший научный сотрудник лаборатории прикладной информатики РОССИЯ, 199178, Санкт-Петербург, 14-линия 39 E-mail: [email protected]
SCENARIO OF INVOLUNTARY DESTRUCTION OF A POPULATION IN A MODIFIED HUTCHINSON EQUATION
Perevarukha A. Yu.
The problem of simulating abrupt changes in the mode of self-oscillations inherent in species that are capable of affecting their habitat is considered. The relevance of this work is the need to improve methods of mathematical biology to study non-stationary and extreme types of population dynamics, which often occur in practice. Rapid transitions to sharp fluctuations in the number of infestations occur during invasions of actively breeding pest species as Ostrinia nubilalis. Modification of the Hutchinson equation suggested in the article, taking into account a significant role of achieving the subthreshold point number that is less than the limiting capacity of the ecological niche K in Verhulst equation, but a significantly higher number of lower threshold L in Bazykin equation: L C H < K. In our equation, the atypical scenario of the development of a dangerous outbreak of insects is described with the change in the acting delay of regulation r value. As follows from ecological examples, population cycles with large amplitude are often unstable. Often the cycle is a transitional mode. The smooth damping of the oscillations N*(r,t) ^ K does not always occur. In the new model, after the Hopf bifurcation, with the value f = t* + £ and the appearance of auto-oscillations of the nonharmonic form with increasing amplitude, the loss of the dissipative property of the trajectory sharply occurs. The computational scenario with the sudden output of the transient cycle N* (fr, t) from the range of admissible values of abundance is interpreted as a specific disturbance of the functioning of the habitat. The loss of a compact attracting set leads to the destruction of the biosystem in the locus of an outbreak of insects or irretrievable death in the case of an island population of mammals.
Key words: Hutchinson equation, dynamics of pest insects, model of exclusive outbreak, Hopf bifurcation, unstable cycle.
References
1. Jeffers J. An Introduction to Systems Analysis: with Ecological Applications. London: Edward Arnold, 1978.
2. Bacaer N. A Short History of Mathematical Population Dynamics. London: Springer-Verlag, 2011. DOI: 10.1007/978-0-85729-115-8.
3. Bazykin A. D. Theoretical and mathematical ecology: dangerous boundaries and criteria of approach them, Mathematics and Modelling. N. Y.: Nova Sci. Publ. Inc., 1993. P. 321-328.
4. van Aarde R. Culling and the dynamics of the Kruger National Park African elephant population, Animal Conservation, 1999, vol. 2, no. 2, pp. 287-294. DOI: 10.1111/j.1469-1795.1999.tb00075.x.
5. Gilg O., Sittler T., and Hanski I. Climate change and cyclic predator-prey population dynamics in the high Arctic, Global Change Biology, 2009, vol. 15, no. 11, pp. 2634-2652. DOI: 10.1111/j. 1365-2486.2009.01927.x.
6. Nicholson A. An outline of the dynamics of animal populations, Austral. J. Zool., 1954, vol. 2, no. 1, pp. 9-65. DOI: 10.1071 /Z09540009.
7. Hutchinson G. An Introduction to Population Ecology. New Haven: Yale University Press, 1978.
8. Borzdyko V. I. Investigation of the Hutchinson Population Model. Differentsial'nye Uravneniya [Differential EquationsJ, 1985, vol. 21, no. 2, pp. 316-318 (in Russian).
9. Kolesov A. Y., Mishchenko E. F., and Rozov N. K. A modification of Hutchinson's equation, Comput. Math, and Math. PPhys., 2010, vol. 50, no. 12, pp. 1990-2002.
10. Gopalsamy K., Kulenovic M., and Ladas G. Time lags in a «food-limited» population madel, Applicable Analysts, 1988, vol. 31, no. 3, pp. 225-237. DOI: 10.1080/00036818808839826.
11. Cooke B., Nealis V., and Regniere J. Insect Defoliators as Periodic Disturbances in Northern Forest Ecosystems, Plant Disturbance Ecology: the Process and the Response. Elsevier: Burlington, 2007. P. 487-525.
12. Gray D. R. Historical spruce budworm defoliation records adjusted for insecticide protection in New Brunswick, J. Acad. Entomol. Soc., 2007, vol. 115, no. 1, pp. 1-6.
13. Sharkovsky A. N., Kolyada S. F., Sivak A. G., and Fedorenko V. V. Dynamics of One-Dimensional Maps. Dorbrecht: Kluwer Acad. Publ., 1997. 272 pp.
14. Ladas G., Qian С. Oscillation and global stability in a delay logistic equation, Dynamics and Stability of Systems, 1994, vol. 9, no. 2, pp. 153-162. DOI: 10.1080/02681119408806174.
15. Ruan S. Delay differential equations in single species dynamics, Delay Differential Equations and Appl. Berlin: Springer, 2006, pp. 477-517. DOI: 10.1007/1-4020-3647-7^11.
Received June 26, 2016
Perevaryukha Andrey Yurievich
St. Petersburg Institute for Informatics and Automation
of the Russian Academy of Sciences,
Senior researcher of the Laboratory of Applied Informatics 14-th Linia, V.l., 39, St. Petersburg, 199178, Russia E-mail: [email protected]
