ECOLOGICAL AND BIOLOGICAL EXPLANATION OF MATHEMATICAL MODELING OF CTENOPHORE MNEMIOPSIS LEIDYI POPULATION IN THE CASPIAN SEA
Authors:
1.
Astrakhan State Technical University
2.
Astrakhan State Technical University
Russian Federation
Russian Federation
3.
Caspian Research Institute of Fishery
Type:
Article
Pages:
from 47
to 61
Status:
Published
Received:
05.03.2015
Accepted:
25.03.2015
Published:
25.03.2015
Subject area:
UDK 593.8 (262.81)
GRNTI 34.39 Физиология человека и животных
GRNTI 62.13 Биотехнологические процессы и аппараты
GRNTI 69.01 Общие вопросы рыбного хозяйства
GRNTI 69.25 Аквакультура. Рыбоводство
GRNTI 69.31 Промышленное рыболовство
GRNTI 69.51 Технология переработки сырья водного происхождения
GRNTI 87.19 Загрязнение и охрана вод суши, морей и океанов
GRNTI 34.39 Физиология человека и животных
GRNTI 62.13 Биотехнологические процессы и аппараты
GRNTI 69.01 Общие вопросы рыбного хозяйства
GRNTI 69.25 Аквакультура. Рыбоводство
GRNTI 69.31 Промышленное рыболовство
GRNTI 69.51 Технология переработки сырья водного происхождения
GRNTI 87.19 Загрязнение и охрана вод суши, морей и океанов
Language:
Russian
Keywords:
invader, hydrobionts, ctenophores, zooplankton, invasion, forage base, model, trophism, environmental factors, ecosystem
Abstract (English):
Ctenophore Mnemiopsis leidyi was found in the Caspian Sea in 1999 and it immediately proved itself as an invader with aggressive specific r -strategy, which in conditions of a new reservoir (unlimited food supply, lack of predator, parasites and diseases and so on) is increasing its number and biomass very quickly. At low concentrations of food organisms, the physiological adaptation "mechanism of switching of somatic growth on generative" is becoming active. Such mechanism of redistribution of products within the population makes possible to use effectively the limited food resource. Setting up a simulation model of the spatial-temporal distribution of Mnemiopsis leidyi population, based on the cybernetic approach, is impossible without taking into account seasonal and interannual trophic dynamics for the entire water area of the Caspian Sea. The invasion of Mnemiopsis leidyi has resulted in serious irreversible changes at all levels, from populations to ecosystems of the whole sea. Therefore, modeling invasions is a promising and important way of studying of the ecosystems of the Caspian Sea, being of great practical and theoretical interest. The paper describes one of the first experiences of the researches in modeling of spatial-temporal dynamics of settling of the Caspian population of Mnemiopsis leidyi from wintering areas of the sea. The scheme of the imitation model of the population of Mnemiopsis leidyi and the process of its migration after wintering, using the refined submodels of trophic dynamics and water cycle, is presented. It allows to determine the rates and terms of the spring-summer distribution of Mnemiopsis leidyi ; the succession of inhabiting the parts of the sea; the degree of influence on the development of the nutritive base; the degree of the influence on the number of populations, commercial stocks of fish species of the Caspian Sea. Besides, the simulation results can be used to predict and develop new invasion models of invading species into other water bodies.
Ctenophore Mnemiopsis leidyi was found in the Caspian Sea in 1999 and it immediately proved itself as an invader with aggressive specific r -strategy, which in conditions of a new reservoir (unlimited food supply, lack of predator, parasites and diseases and so on) is increasing its number and biomass very quickly. At low concentrations of food organisms, the physiological adaptation "mechanism of switching of somatic growth on generative" is becoming active. Such mechanism of redistribution of products within the population makes possible to use effectively the limited food resource. Setting up a simulation model of the spatial-temporal distribution of Mnemiopsis leidyi population, based on the cybernetic approach, is impossible without taking into account seasonal and interannual trophic dynamics for the entire water area of the Caspian Sea. The invasion of Mnemiopsis leidyi has resulted in serious irreversible changes at all levels, from populations to ecosystems of the whole sea. Therefore, modeling invasions is a promising and important way of studying of the ecosystems of the Caspian Sea, being of great practical and theoretical interest. The paper describes one of the first experiences of the researches in modeling of spatial-temporal dynamics of settling of the Caspian population of Mnemiopsis leidyi from wintering areas of the sea. The scheme of the imitation model of the population of Mnemiopsis leidyi and the process of its migration after wintering, using the refined submodels of trophic dynamics and water cycle, is presented. It allows to determine the rates and terms of the spring-summer distribution of Mnemiopsis leidyi ; the succession of inhabiting the parts of the sea; the degree of influence on the development of the nutritive base; the degree of the influence on the number of populations, commercial stocks of fish species of the Caspian Sea. Besides, the simulation results can be used to predict and develop new invasion models of invading species into other water bodies.
Keywords:
invader, hydrobionts, ctenophores, zooplankton, invasion, forage base, model, trophism, environmental factors, ecosystem
invader, hydrobionts, ctenophores, zooplankton, invasion, forage base, model, trophism, environmental factors, ecosystem
Введение Биологическое загрязнение является одной из серьезных мировых экологических проблем, т. к. в настоящее время балластные воды морских судов - наиболее часто встречающийся способ инвазии водных вселенцев [1]. С точки зрения биологического загрязнения моделирование инвазий - перспективное и актуальное направление морской экологии, имеющее для экологов большой практический и теоретический интерес. Существует острая необходимость в разработке имитационной модели всей экосистемы Каспийского моря. Для такой сложной задачи предусматривается поэтапное решение. Моделирование такой крупной структурной единицы, как морская экосистема должно базироваться на более мелких моделях, таких как подмодели окружающей среды и подмодели популяций основных гидробионтов Каспийского моря, в том числе - популяции гребневика Mnemiopsis leidyi (рис. 1, 2). Гребневик был обнаружен в Каспийском море в 1999 г. [2] и сразу же зарекомендовал себя как вселенец с агрессивной видовой r-стратегией [3]. Ранее представители данного типа гидробионтов (тип Ctenophora) в Каспийском море не встречались [4]. Появление гребневика привело к необратимым серьёзным изменениям на всех уровнях - от отдельных популяций до экосистемы всего моря. Рис. 1. Внешний вид вселенца в Каспийское море - гребневика Mnemiopsis leidyi (фото В. Б. Ушивцева, А. М. Камакина) Рис. 2. Внутреннее строение Mnemiopsis leidyi: 1 - аборальный орган; 2 и 3 - гребные пластинки; 4 - аурикулюс; 5-7 - каналы; 8 - боковая лопасть (рис. Т. А. Шигановой) Видовые особенности выживания Mnemiopsis leidyi Чем же определена агрессивность вида Mnemiopsis leidyi? Во-первых, гребневики являются гермафродитами, т. е. для создания новой популяции достаточно одной особи, а не двух. Во-вторых, это вид с простым жизненным циклом - без метаморфоз [5], т. е. простой короткоцикличный вид, с необычайно высокой плодовитостью - до 10 000 яиц/сут, что позволяет ему лавинообразно, в геометрической прогрессии, наращивать свою численность. На первом этапе моделирования необходимо знать энергетические затраты на продуцирование яиц и среднюю плодовитость каждой размерной группы популяции или хотя бы среднепопуляционную в каждом временном шаге модели. Это учитывается в приведенных ниже формулах (8)-(11). Необходимо отметить, что особенностью данного вида является педогенез - размножение организма на стадии личинки (когда личинка размножается, доля половозрелой части популяции в модели соответственно увеличивается). У личинок «родительского» поколения с дальнейшим ростом тела гонады рассасываются, но через несколько недель, уже у взрослого гребневика, они вновь формируются, и он опять готов к репродукции [6]. Мнемиопсис является активным хищником, захватывает добычу с помощью лопастей (рис. 2), покрытых клейкой слизью. Естественно, площадь лопастей у групп старшего возраста больше и в течение жизни при дальнейшем росте организма происходит ее увеличение, за счет чего процесс питания протекает эффективнее и интенсивнее. Соответственно, это биологическое свойство в моделировании также необходимо учитывать в функции «Зависимость линейных размеров гребневика и интенсивность потребления зоопланктона». В отличие от медуз Mnemiopsis leidyi является макрофагом, т. е. способен поедать относительно крупную добычу - длиной до 1 см и более [7]. Например, у особей каспийской популяции в желудках попадались личинки макрозообентоса, икра и личинки рыб (рис. 3, 4). Таким образом, с увеличением размеров организма хищника размерный состав пищи также меняется, т. е. вероятность поедания крупной добычи увеличивается. Рис. 3. Головная часть личинки анчоусовидной кильки (Clupeonella engrauliformes, Borodin) (фото А. М. Камакина) Рис. 4. Внешний вид пелагической икры каспийских килек (р. Clupeonella) (фото А. М. Камакина) Благодаря вышеперечисленным эколого-биологическим свойствам вида популяция мнемиопсиса в нативном ареале (залив Наррагансетт, штат Род-Айленд, США) в течение 1-2 месяцев может увеличивать свою численность примерно в 100 000 раз (на 5 порядков) [8-11]. С началом массового развития в 2003 г. Mnemiopsis leidyi в Каспийском море наблюдалось резкое снижение уловов его пищевых конкурентов. В наибольше степени это отразилось на запасах анчоусовидной кильки (Clupeonella engrauliformes, Borodin), т. к. данный вид каспийских килек (р. Clupeonella) имеет наиболее сходный спектр питания [12, 13]. Из этого следует, что ретроспективное поэтапное моделирование, т. е., в нашем случае, с момента инвазии Mnemiopsis leidyi и до настоящего времени, является одним из направлений современных научных исследований в экологии по выяснению последствий антропогенного влияния биологического загрязнения на промысловые запасы гидробионтов. Для наиболее полного изучения процесса инвазии гребневика необходимо создать пространственно-временную имитационную модель экосистемы всего Каспийского моря. Однако столь сложную модель можно воплотить только через детализацию и включение в нее подмодели отдельных более мелких геоморфологических элементов (Северный, Средний и Южный Каспий) [14] и экосистемных структурных элементов (популяция, субпопуляции, виды, сообщества и т. п.). В исследованиях на Каспии, как и в аналогичных исследованиях Азово-Черноморского бассейна [15], на первом этапе создания модели возможно применение агрегированной 2-компартментальной модели экосистемы, в частности Северного Каспия. При моделировании кормовой базы популяции Mnemiopsis leidyi в динамике в качестве исходных берутся результаты агрегированной модели [16], а также натурные и экспертные оценки распределения трофности всей акватории ареала гребневика с районированием моря на сегменты (рис. 5) [14]. Рис. 5. Районирование акватории Северного Каспия на отдельные районы-сегменты Основные принципы моделирования экосистемы Каспия Трофическая модель динамики популяции Mnemiopsis leidyi является частным случаем трофодинамики всей экосистемы пелагиали Каспийского моря, она базируется на знании структуры и связей пищевых цепей, представленных в виде совокупности основных элементов, между которыми осуществляется обмен веществом и энергией в соответствии с их направлением и объемом [15]: dx/dt = (P - Pt) e + u - v, (1) где x = (x1, x2, …, xn) - вектор состояния, исходными данными являются значения биомассы основных элементов экосистемы на единицу площади акватории в основных районах-сегментах моря; P = (Pij) - матрица прямых энергетических потоков (от j-го элемента экосистемы к i-му); u, v - векторы входящих и исходящих положительных экзогенных потоков (нерест, естественная и промысловая смертность, миграция и т. п.). Матричное представление формулы (1) основано на допущении аддитивности энергетических потоков. При этом все компоненты удобнее выражать в единицах эквивалента энергии органического вещества (1 г Сорг = 10 ккал) [15]. Так как рассматриваемая система характеризуется диссипацией энергии на всех уровнях экосистемы, то используется неявная схема [15], в которой частные потоки энергии выражаются формулой Pij = gij · xj, (2) где gij - скорость трансформации j-го элемента в i-й. Далее, применяя неявную схему Эйлера для интегрирования формулы (1) с шагом t, получаем устойчивый алгоритм, который выступает в качестве самостоятельной дискретной динамической подмодели [16, 17]: xt + t = B (xt + P + U) ; (3a) P = A (xt + P + U); (3b) где xt = (xt1, xt2, ..., xtn)T; P = (P1, P2, ..., Pn)T; ; где T - символ транспонирования. Этот шаг построения модели заключается в параметризации скоростей gij, характеризующих основные эколого-биологические процессы в системе (выедание, отмирание, метаболизм и т. п.). При рассмотрении пространственной неоднородности необходимо также учитывать при водообмене между районами моря физический перенос гидробионтов (активную и пассивную миграцию). Наряду с гладкими функциональными зависимостями типа Холлинга, Контуа и т. д., трофические функции (скорость питания и роста) отражают один из основополагающих принципов экологии - принцип минимума Либиха. В соответствии с ним скорость питания аппроксимируется выражением [15, 18]: m (x, d) = min(1, x, d), (4) где x = y/y*, d = g/g*, - приведенная концентрация пищи (с учётом весовых коэффициентов предпочтения wj); g = y/x - обеспеченность пищей; x*, y* - «пороговые» концентрации хищника и жертвы соответственно; g* = y*/x*. Использование негладкой функции (4) и безразмерного показателя m как индикатора условий питания позволяет дифференцировать трофическое пространство, что облегчает идентификацию и анализ модели [15]. При лимитирующей концентрации пищи x минимум достигается на 2-м члене, при низкой обеспеченности d - на 3-м. Соответственно, при ухудшении трофических условий рацион гидробионтов в среднем уменьшается, а при ограничении линейным приближением в окрестности порога насыщения m = 1 получается формула [15, 18]: С = M [1 - a (1 - m)] x, (5) где M - максимальная скорость питания; 0 < a< 1 - коэффициент чувствительности рациона к изменению условий питания. Так как темп весового роста отдельных особей не может сильно отклоняться от некоторого условно среднестандартного, то коэффициент a лежит в пределах 0 < a ≤ 0,2. При дефиците пищи снижается темп роста и, соответственно, растёт популяционная смертность. Такую добавочную смертность D можно описать линеаризованной зависимостью [15] в окрестности m = 1: D = M b (1 - m) сx, (6) где b - коэффициент чувствительности популяционной смертности к дефициту пищи. С учетом уравнений (4)-(6) можно ввести модифицированную трофическую функцию (зависимость наиболее эффективного удельного рациона от оптимальных трофических условий нагула): f(m) = (C - D)/x = Mmin {0,1 - (a + b) (1 - m)}. (7) Если a + b > 1, то f(m) приобретает S-образный вид вследствие образования «трофического» порога [15, 18]. Имитационная модель каспийской популяции Mnemiopsis leidyi Mnemiopsis leidyi является типичным зоопланктофагом, и, соответственно, его можно отнести по ряду параметров к группе видов с агрессивной стратегией развития (r-стратегией) [3]. Такой вид, попадая в благоприятные условия нового водоема (неограниченные пищевые ресурсы, отсутствие хищников, паразитов и болезней и т. п.) не только выживает, но и очень быстро наращивает численность и биомассу популяции [19]. Данный вид-вселенец обладает специализированным адаптационным механизмом, позволяющим ему выживать в неблагоприятных условиях нового водоема. Например, при низкой концентрации кормовых организмов запускается адаптационный физиологический механизм переключения соматического роста на генеративный. Кроме этого, у рыб разных возрастов существуют различия в рационе питания. Самые младшие по возрасту группы этого вида гребневика питаются микропланктоном: простейшими, фитопланктоном, науплиальными стадиями веслоногих рачков (отр. Copepoda), встречающимися в достаточном количестве. Кроме вышеперечисленного, повышению выживаемости способствует увеличение численности популяции за счет педогенеза [1, 5, 20-22]. Данный эколого-биологический популяционный процесс можно описать с помощью следующей математической подмодели [15]: dx1/dt = (1 + s) Pg - e1x1 - dx1; (8) dx2/dt = P - Pg - e2x2 + dx1; (9) dy/dt = r K(t) y / {K(t) + [(r/e k) - 1] y} - e y - C; (10) Pg = r (y) P, (11) где x1 - биомасса неполовозрелой части популяции; x2 - биомасса половозрелых особей; y - приведенная биомасса кормовых (планктонных) организмов; С - выедание Mnemiopsis leidyi зоопланктона; e1, e2 - скорость отмирания Mnemiopsis leidyi; ek - скорость отмирания зоопланктона; d - скорость перехода организмов Mnemiopsis leidyi в половозрелую группу; P и Pg - чистая и генеративная продукция Mnemiopsis leidyi; s - весовой прирост неполовозрелой части популяции; r (y) - генеративная часть общей продукции, зависящая от концентрации зоопланктона; r - скорость-брутто воспроизводства биомассы зоопланктона; K(t) - динамика биомассы зоопланктона при отсутствии Mnemiopsis leidyi. При расчете продукции популяции Mnemiopsis leidyi в качестве трофической функции используется формула (7) [15], в итоге получаем следующую формулу: P = k1 k2 f(m) x2, (12) где k1, k2 - коэффициенты эффективности процессов ассимиляции и роста, а выедание зоопланктона С, в свою очередь, задается формулой (5). Учитывая тесную функциональную связь численности Mnemiopsis leidyi с факторами среды, можно предположить, что в формуле (7) a = 0, b = 1. С учетом педогенеза неполовозрелая стадия непродолжительна (d >> e), то x1 << x2. Если для упрощения модели принять, что e1 = e2 = e, то, суммируя уравнения (8)-(9) и с учётом формул (11)-(12), в первом приближении получаем dx/dt = M1m [1 + sr (y)] x - ex, (13) где x = x1 + x2 - биомасса всей популяции Mnemiopsis leidyi: M1 = Mk1 k2. Согласно вышеописанной гипотезе, получаем dr(y) / dy < 0. (14) При снижении концентрации пищи вышеприведенные формулы описывают механизм активации процесса генеративного роста, что позволяет сохранить высокую скорость образования продукции [15]. В гипотезе медленного изменения концентрации K(t) (квазистационарности), выступающей в данной модели в качестве экзогенной переменной, система уравнений (10) и (13) может быть описана методами качественной теории дифференциальных уравнений. Указанная система имеет 3 состояния равновесия, одно из которых тривиальное (x = 0; y = K(t)). Таким образом, условием вселения Mnemiopsis leidyi является неустойчивость тривиального равновесия [15]. Данное условие достигается, если в момент вселения имеется достаточное количество кормовых организмов K(t). Критическую концентрацию корма можно найти из уравнения (15): min (1, K/y) [1 + sr (K)] = e /M1. (15) Необходимо указать, что в уравнении (14) единственное нетривиальное устойчивое равновесие сдвигается в сторону уменьшения концентрации y и роста биомассы x. Это свидетельствует о том, что внутрипопуляционный механизм перераспределения продукции позволяет наиболее эффективно использовать пищевой ресурс, имеющийся на данный момент времени [15]. Пространственная подмодель каспийской популяции Mnemiopsis leidyi Представители типа Ctenophora в поисках пищи практически не совершают продолжительных направленных перемещений, за исключением вертикальных [5]. На Каспии Mnemiopsis leidyi совершает сезонные перемещения по акватории море за счет пассивного переноса течениями [23], поэтому, при разработке пространственно-временной модели популяции гребневика, миграцию можно рассматривать как пассивный элемент системы [15]. В основе пространственной модели пелагиали Каспийского моря необходимо положить кибернетический подход, как и для других бассейнов, например, Азово-Черноморского [15]. Акватория Северного Каспия поделена на 8 двухслойных районов-сегментов (рис. 5), водообмен между которыми осуществляется согласно определенной схеме [24, 25]. С точки зрения гидрологии толщина верхнего, динамически активного слоя моря равна 20-25 м. С точки зрения биологии вида минимальная нижняя граница вертикального распространения Mnemiopsis leidyi в летний период ограничена термоклином, т. е. глубиной 20-30 м [23]. Зимой, при отсутствии термоклина, граница проникновения гребневика опускается до 100-метровой глубины [26]. Летом и осенью нижняя граница биотопа вида Mnemiopsis leidyi достигает дна в мелководных районах юго-восточного шельфа Южного Каспия, северо-западного шельфа Среднего Каспия, а также почти всей акватории Северного Каспия. В схеме описания циркуляции вод Каспийского моря имеется два антициклонических круговорота в Южнокаспийской и Дербентской котловинах, осуществляющих перенос водных масс в направлении против часовой стрелки вдоль береговой линии. Общепринято, что схема круговых течений неизменна в течение всего теплого времени года, а зимой, в присутствии вертикальной гомотермии, на восточном мелководье Южного Каспия образуется относительно теплое течение, идущее вдоль всего восточного побережья и достигающее в Северном Каспии границ льдообразования [24, 25]. На остальной акватории моря изменения касаются только интенсивности водообмена. Динамика вертикального водообмена тесно связана со сроками и пространственными особенностями сезонного формирования слоя термоклина. Согласно гидрологическим данным, подъем вод происходит в центрах круговоротов [24, 25]. Ветровые течения считаются стохастической компонентой, которая учитывается при определении водообмена между Средним и Северным Каспием, а также внутри мелководного Северного Каспия - между соседними квадратами-сегментами. На макроуровне водообмен отображает диффузионные процессы [14]. Масштаб процесса диффузии определяется размерами каждого из 8 районов-сегментов (рис. 5). Географический размер квадрата составляет 10"×10", на широте Каспийского моря в метрических единицах его стороны соответственно составили 10 и 7 морских миль или 18,52 и 13,00 км. Настройку и верификацию подмодели водообмена в Северном Каспии можно проводить по среднемноголетним показателям солёности [15]. Идентифицированная таким образом подмодель используется для моделирования процесса пассивного расселения Mnemiopsis leidyi по акватории моря из районов зимовки. Имитационное моделирование процесса распространения популяции Mnemiopsis leidyi в Каспийском море Для оценки экологических параметров популяции вида Mnemiopsis leidyi (табл.) в качестве калибровочного используется коэффициент смертности [15, 27]. Выявление картины пространственно-временного распределения популяции Mnemiopsis leidyi невозможно без учета сезонной и межгодовой трофодинамики. Проводится сравнительный анализ развития кормовой базы, соответствующей двум основным периодам - до инвазии гребневика Mnemiopsis leidyi (1980-1999 гг.) [28, 29], и после (2000-2014 гг.) [30]. Например, в эксперименте без вселения Mnemiopsis leidyi в Азово-Черноморском бассейне расчетные данные были получены с помощью компьютерной модели MT Base 1.1 [15]. За точку отсчета, т. е. год первого появления гребневика в Каспийском море, должен приниматься 1999 г. - год первого официального обнаружения этого вселенца [2]. По данным проф. С. П. Воловика [1], начальное (стартовое) значение биомассы должно задаваться величиной равной 50 кг сырого веса на 1 м3, что в основных районах ареала в расчёте на весь верхний слой соответствует концентрации 10-6 кал/м3 [27]. Основные эколого-биологические параметры вида популяции Mnemiopsis leidyi в Черном и Каспийском морях [15, 28, 29] Показатель Черное море Каспийское море Единица измерения Интервал оптимальности по температуре 14-25 16-28 оС Интервал толерантности по температуре 5-30 6-32 оС Интервал оптимальности по солёности 10-20 9,0-14,0 ‰ Интервал толерантности по солёности 6-30 3,5-14,5 ‰ Максимальная скорость потребления пищи 50 50 1/год Коэффициент эффективности усвоения k1 0,75 0,75 1/год Коэффициент эффективности роста k2 0,60 0,60 1/год Минимальная смертность 4,2 4,7 1/год Максимальная смертность 9,0 9,5 1/год Гребневик Mnemiopsis leidyi в Каспии попал, естественно, в новые для себя трофические условия [2]. Соответственно, параметр «пороговая» концентрация пищи» y, при которой начинается лимитирование при дальнейшем росте популяции, определяется эмпирически и рассматривается в качестве калибровочного [15]. Так, при больших значениях данного параметра (> 300 кал/м3) Mnemiopsis leidyi поступал в Северный Каспий раньше, а при низких значениях (< 50 кал/м3) его вселение происходит позже, например в 2002 г. [26]. В итоге теоретически определено адекватное для модели значение этого параметра - 200 кал/м3 [15]. Многолетняя динамика. В зависимости от района моря максимальное значение численности Mnemiopsis leidyi наблюдается в августе - сентябре: до 820 экз./м3 (рис. 6). В северной части Каспийского моря в октябре - ноябре начинает прослеживаться снижение (в 10-20 раз) биомассы. В последние годы ежегодный максимум постоянно наблюдается в районе мелководий (северо-западная часть Среднего Каспия и юго-восточная - Южного Каспия [23]. Рис. 6. Многолетняя динамика (2001-2012 гг.) Mnemiopsis leidyi, август - сентябрь Многолетняя динамика распределения Mnemiopsis leidyi имеет следующие закономерности: - в первые годы после вселения (предположительно 1996-1998 гг.) развитие популяции Mnemiopsis leidyi носит скрытый характер (латентная фаза), и в этот период она сосредоточена в наиболее благоприятной по условиям среды обитания южной части Каспийского моря; - в Южном Каспии - в исходной области образования новой популяции Mnemiopsis leidyi; максимум численности наблюдался на 5-6 год после вселения, т. е. в 2002-2003 гг. (рис. 6); - на 7-8 год отмечается снижение ее численности на порядок (в 15-16 раз) с последующей стабилизацией на этом уровне [23]. Сезонная динамика. Анализ карт-схем сезонного распределения Mnemiopsis leidyi по акватории моря (рис. 7) выявил область постоянного местоположения зимнего ареала - это открытая часть Южного Каспия [26]. а б Рис. 7. Карта распределения численности Mnemiopsis leidyi каспийской популяции (слой 0-100 м) в аномальные зимы, экз./м3: а - холодная, 2002 г.; б - теплая, 2004 г. Особенности летнего распределения (рис. 8) несколько иные - в южной части Каспийского моря субпопуляция Mnemiopsis leidyi образует 2 ядра. Первое ядро (биомассы) формируется на северо-западе Среднего Каспия. Оно формируется при высоких концентрациях зоопланктона и предопределено внутрипопуляционным процессом соматическим роста. Второе (ядро численности) образуется на западе и юго-востоке Южного Каспия. Образование второго ядра является следствием лимитирования кормовой базы, сдерживающей увеличение биомассы хищника-зоопланктофага, но достаточной для поддержания высокого уровня прироста за счет численности. Соответственно, из этой области осуществляется пассивная миграция (за счет течений) в малозаселенные районы, например в восточные и северо-восточные Среднего Каспия [23]. В Южном Каспии зимнее распределение Mnemiopsis leidyi как хищника соответствует распределению его жертвы - Larvae Lamellibranchiata. Данное явление наблюдается только в холодный период года (с ноября по июнь), что является сезонной особенностью распределения популяции Mnemiopsis leidyi (рис. 9). Судя по картам, максимальные концентрации биомассы гребневика Mnemiopsis leidyi наблюдались в центральном глубоководном районе моря, а минимальные - в районе восточного мелководья у о-ва Огурчинский [12, 13]. Есть и другие особенности распределения мнемиопсиса, например пути весенних миграций идут не вдоль восточного побережья, где можно использовать циклоническое течение, а вдоль западного, т. е. за счет системы небольших противотечений, направленных против основной - волжской струи. а б Рис. 8. Карта распределения численности каспийской популяции Mnemiopsis leidyi (слой 0-50 м) в 2006 г., экз./м3: а - в июне; б - в августе а б Рис. 9. Карта зимнего распределения биомассы в Каспийском море (слой 0-100 м) в январе 2002 г., г/м3: а - хищника (Mnemiopsis leidyi); б - жертвы (Larvae Lamellibranchiata) Многолетние данные показали, что в течение всего вегетационного сезона у западного побережья Среднего и Северного Каспия популяция Mnemiopsis leidyi имеет более высокие значения концентрации, чем у восточного [23]. Заключение Имитационная модель популяции Mnemiopsis leidyi основана на закономерностях сезонной динамики процесса его весенних и осенних перемещений. В совокупности с использованием детализированной подмодели трофической динамики и подмодели водообмена эта модель позволит определять: - темпы и сроки весенних и осенних миграций Mnemiopsis leidyi; - последовательность заселения районов моря; - степень влияния на развитие кормовой базы; - степень влияния на численность и промысловые запасы популяций рыб Каспийского моря. Выделяются четыре основных этапа годового цикла развития популяции Mnemiopsis leidyi в Каспийском море: - 1-й этап - пассивная миграция из районов зимовки; - 2-й этап - наращивание численности и биомассы; - 3-й этап - достижение летнего пика развития; - 4-й этап - осенний спад численности и биомассы и зимнее отмирание части популяции в Северном Каспии и в прибрежной мелководной части Среднего и Южного Каспия. В заключение можно сказать, что, кроме исследований в других морях [15, 16, 31], данная работа является одним из первых опытов исследований по моделированию пространственно-временной динамики годового цикла развития каспийской популяции Mnemiopsis leidyi. Необходимо отметить, что в дальнейшем через определенный период постоянно будет требоваться корректировка модели, т. е. исходная информация постоянно будет нуждаться в уточнении в ходе натурных экспедиционных исследований. Опыт исследований в других бассейнах показывает, что результаты работы модели достаточно объективны и адекватны и, соответственно, достаточно близки к натурным исследованиям морских экосистем [15, 16]. Сами результаты и методика моделирования экосистемы Каспийского моря и отдельных ее элементов могут быть использованы для прогнозирования численности и биомассы как гидробиологических сообществ в целом, так и отдельно - составляющих элементов (виды, популяции и т. п.), в частности - для прогнозирования инвазий вселенцев в новые водоемы.
1. Harbison G. R. The ctenophore, Mnemiopsis leidyi, in the Black Sea: a holoplanktonic organism transported in the ballast water of ships / G. R. Harbison, S. P. Volovik // Nonindigenous Estuarine and Marine Organisms (NEMO). Proceedings of the Conference and Workshop, Seattle, Washington, April 1993, U. S. Department of Commerce, National Oceanic and Atmospheric Administration, Office of the Chief Scientist (September 1994). U. S. Government Document N C55.2. N 73. P. 25-36.
2. Ivanov V. P. Simultaneous invasion of the Caspian Sea by two jellies Mnemiopsis end Aurelia / V. P. Ivanov, A. M. Kamakin, V. B. Ushivtzev, T. A. Shiganova, O. Zhukova, N. Aladin, I. Susan, S. I. Wilson, G. R. Harbison, H. J. Dumont // J. of Invasion. In: Biological Invasions. Kluwer Academic Publishers, Printed in the Netherlands. 2000. N 2. P. 255-258.
3. Grin N. Biologiya / N. Grin, U. Staut, D. Teylor: v 3 t. / pod red. R. Sopera. M.: Mir, 1990. T. 2. 325 s.
4. Atlas bespozvonochnyh Kaspiyskogo morya / pod red. prof. Ya. A. Birshteyna, L. G. Vinogradova i dr. M.: Pisch. prom-st', 1968. 416 s.
5. Dogel' V. A. Zoologiya bespozvonochnyh: ucheb. dlya un-v / V. A. Dogel': pod red. prof. Yu. I. Polyanskogo. M.: Vyssh. shk., 1975. S. 130-134.
6. Backer L. D. Laboratory culture of lobate ctenophore Mnemiopsis mccradyi with notes on feeding and fecundity / L. D. Backer, M. R. Reeve // Marine Biol. 1974. Vol. 26, no. 1. P. 57-62.
7. Main R. J. Observations of the feeding mechanism of a ctenophore, Mnemiopsis leidyi / R. J. Main // Biol. Bull. 1928. Vol. 55, no. 2. P. 69-78.
8. Kremer P. M. Predation by the ctenophore Mnemiopsis leidyi in Narragansett Bay, Rhode Island / P. M. Kremer // Estuaries. 1979. Vol. 2, no. 2. P. 97-105.
9. Oviatt C. A. Perdition on the ctenophore Mnemiopsis leidyi by butterfish Peprilus triacanthus in Narragansett Bay, Rhode Island / C. A. Oviatt, P. M. Kremer // Chesapeake Sci. 1977. Vol. 18, no. 2. P. 236-240.
10. Deason E. E. Mnemiopsis leidyi (Ctenophora) in Narragansett Bay, 1975-1979: Abundance, size composition and estimation of grazing / E. E. Deason // Estuar. Coast. Shelf Sci. 1982. Vol. 15, no. 2. P. 121-134.
11. Deason E. E. Experimental evaluation of herbivory in the ctenophore Mnemiopsis leidyi relevant to ctenophore - zooplankton - phytoplankton interactions in Narragansett Bay, Rhode Islend, USA / E. E. Deason, T. J. Smauda // J. Plankton Res. 1982. Vol. 4, no. 2. P. 219-236.
12. Kamakin A. M. Vliyanie Mnemiopsis leidyi na nekotorye troficheskie zven'ya Kaspiyskogo morya / A. M. Kamakin, O. A. Chizhenkova, V. F. Zaycev // Yug Rossii: Ekologiya, razvitie. 2010. № 2. S. 33-42.
13. Kamakin A. M. The impact of invader Mnemiopsis leidyi on aboriginal fauna of the Caspian Sea / A. M. Kamakin, U. A. Parirskiy, L. V. Nikulina // IV International Symposium "Invasion of alien species in Holarctic" (Borok-4), 2013, Russia: Programme & Book of abstracts. Yaroslavl: Publisher's bureau «Filigran'», 2013. P. 69.
14. Kamakin A. M. Elementy imitacionnoy gidrologicheskoy modeli Severnogo Kaspiya kak sredy obitaniya subpopulyacii grebnevika Mnemiopsis leidyi / A. M. Kamakin, D. N. Katunin // Voprosy sovremennoy nauki i praktiki. Problemy, tehnologii, upravlenie: sb. nauch. st. Bugul'ma: NO «FƏN-NAUKA», 2012. S. 20-25.
15. Berdnikov S. V. Imitacionnaya model' osnovnyh komponentov ekosistemy Ohotskogo morya / S. V. Berdnikov, Yu. A. Dombrovskiy, A. G. Ostrovskaya, M. V. Prihod'ko, L. I. Titova, Yu. V. Tyutyunov // Morskoy gidrofizicheskiy zhurnal. 1989. № 3. S. 52-57.
16. Lebedeva L. P. Izmenchivost' chislennosti i biomassy grebnevika-mnemiopsisa v Chernom more (model'noe issledovanie) / L. P. Lebedeva //Okeanologiya. 1998. T. 38, № 5. C. 727-733.
17. Selyutin V. V. Krugovorot veschestva i potok energii v ekologicheskih sistemah: ot modeli sistemy k sisteme modeley / V. V. Selyutin // Obozrenie prikladnoy i promyshlennoy matematiki. 1994. T. 1, № 6. C. 957-973.
18. Gil'derman Yu. I. Modeli L-sistem (sistemy s limitiruyuschimi faktorami) / Yu. I. Gil'derman, K. N. Kudrina, I. A. Poletaev // Issledovaniya po kibernetike. M.: Sov. radio, 1970. S. 165-209.
19. Karpevich A. F. Izbr. tr. / A. F. Karpevich: v 2 t. T. 2. Akklimatizaciya gidrobiontov i nauchnye osnovy akvakul'tury. M.: Izd-vo VNIRO, 1998. 870 s.
20. Mianzan H. W. Estudio preliminar sobre distribucion y abendancia de Mnemiopsis mccradyi en el estuario de Bahna Blanca, Argentina (Ctenophora) / H. W. Mianzan, M. L. Sabatini // Spheniscus. 1985. N 1. P. 53-58.
21. Harbison G. R. On the natural history end distribution of oceanic ctenophores / G. R. Harbison, L. P. Madin, S. N. Rwanberg // Deep-Sea Res. 1978. Vol. 25. P. 233-256.
22. Harbison G. Not all ctenophores are hermaphrodites. Studies on the systematics, distribution, sexuality and development of two species of Ocyropsis / G. Harbison, R. Miller. Mar. Biol. 1986. Vol. 90, iss. 3. P. 413-424.
23. Kamakin A. M. Zakonomernosti mnogoletney i mezhsezonnoy dinamiki populyacii grebnevika Mnemiopsis leidyi v Kaspiyskom more / A. M. Kamakin, V. F. Zaycev // Yug Rossii: ekologiya, razvitie. 2012. № 1. S. 96-103.
24. Gidrometeorologiya i gidrohimiya morey. T. VI. Kaspiyskoe more. Vyp. 1. Gidrometeorologicheskie usloviya / otv. red. F. S. Terziev. SPb.: Gidrometeoizdat, 1992. 360 s.
25. Kaspiyskoe more: gidrologiya i gidrohimiya / pod red. S. S. Baydina, A. N. Kosareva. M.: Nauka, 1986. 261 s.
26. Kamakin A. M. Sostoyanie populyacii Mnemiopsis leidyi v Kaspiyskom more v zimniy period / A. M. Kamakin, O. A. Chizhenkova, V. F. Zaycev // Vestn. Astrahan. gos. tehn. un-ta. Ser.: Rybnoe hozyaystvo. 2009. № 1. S. 58-64.
27. Reeve M. R. Growth dynamics of a ctenophore (Mnemiopsis) in relation to a variable food supply. I. Carbon biomass, feeding, egg production, growth and assimilation efficiency / M. R. Reeve, M. A. Syms, P. Kremer // J. of Plankton Research. 1989. Vol. 11, no. 3, pp. 535-552.
28. Kaspiyskoe more: fauna i biologicheskaya produktivnost' / otv. red. E. A. Yablonskaya i dr. M.: Nauka, 1985. 278 s.
29. Gidrometeorologiya i gidrohimiya morey. T. VI. Kaspiyskoe more. Vyp. 2. Gidrohimicheskie usloviya i okeanologicheskie osnovy formirovaniya biologicheskoy produktivnosti / otv. red. F. S. Terziev. SPb.: Gidrometeoizdat, 1996. 324 s.
30. Polyaninova A. A. Vidy-vselency v Kaspii i ih rol' v ekosisteme morya / A. A. Polyaninova. Astrahan': Izd-vo KaspNIRH, 2007. 104 s.
31. Lebedeva L. P. Model'noe issledovanie vliyaniya grebnevika Mnemiopsis na planktonnoe soobschestvo Chernogo morya / L. P. Lebedeva, E. A. Shushkina // Okeanologiya. 1994. T. 34, № 1. C. 79-87.
