CC BY
13
/?\Л научный электронный журнал
\+У ПРИНЦИПЫ экологии
петрозаводский государственный
VI 1ИВСРСИТЕТ
О
http://ecopri.ru http://petrsu.ru
УДК УДК 574.34:595.789:57.087.1
Опыт оценки абсолютной численности популяции по распределению числа отловов в эксперименте с мечением на примере имаго перламутровок рода БОЬОПА ЦШЕСТА,
ЬЕРЮОРТЕВА, МУМРНАЬЮАЕ)
горбач
вячеслав Васильевич
Петрозаводский государственный университет, gorbach@psu.karelia.ru
Ключевые слова:
метод мечения с повторным отловом особей
распределение Пуассона геометрическое распределение стохастическая модель динамики численности Джолли - Себера Boloria freija Boloria aquilonaris
Получена: 10 апреля 2018
Аннотация: Абсолютную численность популяций животных оценивают по результатам экспериментов с мечением особей. Популяции дневных чешуекрылых обычно изучают методом Джолли - Себера, использующем в качестве математической основы стохастическую модель. Однако не всегда удается получить исходные данные, удовлетворяющие минимальным требованиям к объему выборок, периодичности и пропорциональности отловов. В таких случаях остается возможность оценить абсолютную численность по одному из известных законов распределения частот. В данной работе рассмотрены способы оценки абсолютной численности популяций. Исходные выборки получены методом мечения с повторными отловами по схеме Джолли - Себера в пространственно обособленных группировках бабочек-перламутровок рода Во1опа. Исследования проводили на болотах Близкое и Осоковое в заповеднике «Кивач» летом 1995, 1996 и 2016 гг. Всего помечено 261 имаго Во1опа ^е'а и 3628 имаго В. aquilonaris, число повторных отловов составило 274 и 968 соответственно. В качестве теоретических платформ взяты распределение Пуассона и геометрическое распределение. Общая задача состояла в том, чтобы по доле меченых особей в отловах определить число особей, которые не попались ни разу. Предложен скрипт для расчета абсолютной численности и ее доверительных интервалов по Пуассонову и геометрическому распределениям и оценки адекватности эмпирических частот теоретическим моделям на языке R. Установлено, что Пуассонова модель, предполагающая равную вероятность отлова любой особи, более адекватна экспериментальным данным. Соответствие эмпирических частот этому распределению можно рассматривать в качестве меры изолированности популяции. Сходство Пуассоновых значений абсолютной численности с оценками, сделанными по модели Джолли - Себера, которая предполагает наличие входящих и выходящих потоков особей, ставит под сомнение рекомендации по изначальному выбору определенных методов для «закрытых» и «открытых» популяций. Геометрическую модель используют для описания процессов с непостоянной вероятностью совершения событий, т. е. она должна лучше соответствовать «открытой» системе Джолли - Себера, но на практике эта модель дала сильно завышенные оценки, особенно в приложении к большим по объему выборкам.
© Петрозаводский государственный университет
года Подписана к печати: 25 сентября 2018 года
Введение
В основе одного из наиболее распространенных подходов к оценке абсолютной численности популяций животных лежит метод мечения с повторным отловом особей (MRR-техника: mark - release - recapture). К настоящему времени описано множество модификаций постановки MRR-эксперимента для различных таксономических и экологических групп объектов и выработано немало рекомендаций по использованию математических моделей для аппроксимации результатов исследований. Из руководств на русском языке наиболее популярным среди зоологов остается сводка Г. Коли (1979), из англоязычных изданий - монографии Ч. Кребса (Krebs, 1999), Т. Саусвуда и П. Хен-дерсона (Southwood, Henderson, 2000). Общий принцип всех предлагаемых алгоритмов состоит в том, чтобы по доле меченых особей в отловах определить число особей, которые не были отловлены ни разу, непременным требованием к выборке является случайность отлова любой особи.
Популяции бабочек обычно изучают методом Джолли - Себера, использующем в качестве математической основы стохастическую модель. Пожалуй, это наиболее информативный подход, который по серии многократно повторяемых отловов позволяет оценить целый ряд скрытых демографических переменных, в частности абсолютную численность, численность для каждой даты исследования, скорость роста и т. д. Вместе с тем этот метод обладает рядом особенностей, которые снижают его эффективность. Практика применения указывает на чувствительность используемой модели к малым выборкам, периодичности и пропорциональности отловов (Горбач, 2013). В первом случае бывает невозможно провести расчеты из-за пропуска повторных отловов меченых особей в малочисленных популяциях, другие затруднения обычно возникают из-за неблагоприятных для лета бабочек погодных условий, которые нередко увеличивают временные интервалы между запланированными этапами (датами) мечения или вынуждают досрочно завершить очередной этап, не дав отловить достаточное число особей. В случаях подобных сбоев в ходе эксперимента остается возможность оценить абсолютную численность по одному из известных законов распределения частот. Правильная стратегия в этом случае состоит не в том, чтобы использовать для расчетов
какую-либо модель, бездумно принимая на веру лежащую в ее основе гипотезу, а в том, чтобы из имеющихся вариантов выбрать наилучший по адекватности экспериментальным данным (Коли, 1979).
Цель настоящей работы состоит в сравнении оценок абсолютной численности и в обобщении опыта использования дискретных распределений для описания частот отловов особей в экспериментах с перламутровками рода Boloria. В качестве теоретической платформы мы выбрали две модели -распределение Пуассона и геометрическое распределение, которое в случае с бабочками лучше других описывает результаты экспериментов.
Материалы
Объектами изучения стали два вида бабочек-перламутровок рода Boloria -B. freija (Thunberg, 1791) и B. aquilonaris (Stichel, 1908). Это типичные ледниковые реликты лесной зоны, населяющие в Карелии сфагновые болота. Поскольку болотные массивы представляют собой местообитания островного типа, то и популяции изучаемых видов распадаются на пространственно разобщенные группировки, приуроченные к отдельным выделам. Исследования проводили в заповеднике «Кивач» на сопредельных болотах Близкое и Осоковое, ранее рассматриваемых нами как единый массив (Горбач, 1998, 2011). Данные для анализа выбирали, учитывая дифференциацию объектов по численности: B. aquilonaris является обычным болотным видом, тогда как B. freija - одним из малочисленных. То же касается и рассматриваемых сезонов: 1995 г. стал самым удачным для B. aquilonaris за последние десятилетия, численность же B. freija была настолько низкой, что не имело смысла начинать эксперимент, в сезоне 1996 г. она достигла локального максимума, в 2016 г. соответствовала обычному для этого вида уровню. Общий объем выборок составляет 3889 меченых бабочек (табл. 1 и 2).
Методы
Методология получения выборок подробно описана ранее (Горбач, 2013), она стандартна для подхода Джолли - Себера и полностью соответствует требованиям MRR-эксперимента. Общая задача при оценке абсолютной численности состояла в том, чтобы, зная частоту повторных отловов меченых особей (к = 1, 2, 3, 4, ...), найти частоту нулевого класса (к = 0), т. е. определить
Таблица 1. Распределение частот отлова перламутровок B. freija и оценки абсолютной численности вида на исследуемой территории Table 1. Distribution of the Freija fritillary captures and the estimates of the absolute number of the butterflies in the study area
1996 2016
k -
самцы самки имаго самцы самки имаго
1 61 42 103 19 18 37
2 37 29 66 17 11 28
3 10 3 13 7 3 10
4 1 1 1 1 2
5 1 1
Распределение Пуассона
Np, экз. 178 131 307 60 51 110
lim, экз. 159-206 114-5-157 276-338 54-70 43-69 100-125
X2 2.638 7.123 8.445 1.424 0.244 1.081
p 0.450 0.028 0.037 0.840 0.970 0.897
Геометрическое распределение
Ng, экз. 305 228 524 97 86 182
lim, экз. 266-365 194-282 469-599 84-118 69-121 159-214
X2 10.197 14.307 23.625 6.125 1.975 7.477
p 0.017 <0.001 <0.001 0.190 0.578 0.113
Стохастическая модель динамики численности Джолли - Себера
Njs, экз. 216 114 314 119 102 198
SE, экз. 121 61 141 64 81 93
Примечание. k - частотные классы эмпирического распределения (число отловов особи), Ng, Np, Njs - абсолютная численность, lim - доверительный интервал, SE - стандартная ошибка, х2 и p - оценка и значимость отличий эмпирического и теоретического распределений (критерий Пирсона), жирным шрифтом указаны значимые различия p<0 = 0.05.
Note. k are classes of the empirical frequencies (number of captures of the individual); Ng, Np, Njs - absolute numbers by Poison, geometric and Jolly - Seber models, accordingly; lim - confidence interval; SE - standard error; х2 and p - estimatе and significance of differences between empirical and theoretical distributions (Pearson's criterion), significant differences (р<Ш = 0.05) are indicated in bold.
число особей, которые в ходе эксперимента не попались ни разу. Размер этого класса, а следовательно, и величина искомой переменной (суммы всех частот) зависит от характера распределения.
Распределение Пуассона предполагает постоянную вероятность отлова особей, которая возможна только в «закрытой» популяции при очень низкой, в идеале - нулевой смертности за время эксперимента. Распределение частот представляет собой аппроксимацию биномиального распределения для случаев, когда р<<д (Коросов, 2007), где
р - вероятность найти объект нужного качества, в нашем случае - отловить меченую особь, ц - вероятность не найти такой объект, 0 < р < 1 и ц = 1 - р. Случайная величина принимает при этом целочисленные значения с вероятностью Рк=е-Л-Лк/к!, где Л - параметр распределения, к! - факториал числа к, е - основание натурального логарифма. Абсолютная численность рассчитывается по формуле Np = э/л, где 5 - общее число отловов, ¡л - средняя арифметическая распределения с нулевым классом, которая задается равенством М = ¡/(1 - е-^). Нахождение зна-
Таблица 2. Распределение частот отлова перламутровок B. aquilonaris и оценки абсолютной численности вида на исследуемой территории Table 2. Distribution of the Cranberry fritillary captures and the estimates of the absolute number of the
butterflies in the study area
1, 1995 1996
К самцы самки имаго самцы самки имаго
1 941 985 1926 690 514 1237
2 109 132 241 106 82 188
3 12 9 21 8 5 13
4 2 2
Распределение Пуассона
Np, экз. 4963 4847 9789 3198 2374 5813
lim, экз. 4340^5847 4285^5622 8919^10892 2795^3747 2054^2849 5239^6538
0.723 2.342 0.328 0.898 1.653 2.146
p 0.697 0.504 0.954 0.638 0.199 0.342
Геометрическое распределение
Ng, экз. 9534 9277 18778 6096 4521 11094
lim, экз. 8285^11305 8152^10832 17031^20990 5278^7196 3879^5474 9942^12548
X2 0.256 3.315 2.784 4.765 5.409 9.418
p 0.880 0.346 0.426 0.092 0.067 0.009
Стохастическая модель динамики численности Джолли - Себера
Njs, экз. 3762 3420 7203 2041 1590 3761
SE, экз. 1643 1344 2106 1028 783 1442
Примечание. Обозначения см. в табл. 1. For notations see Table 1.
чения переменной ¡л составляет определенную сложность, поскольку данное уравнение имеет прямое решение лишь в случаях, когда средняя арифметическая усеченного (эмпирического) распределения M > 2. Тогда используют интерполяционные многочлены Лагранжа, в иных случаях искомое значение подбирают методом проб и ошибок (Коли, 1979). Проще эта задача решается средствами оптимизации, встроенными в компьютерные программы, например с помощью пакета «Поиск решения» в среде MS Excel или функции optimize в среде R.
Геометрическое распределение используют для описания процессов, в которых вероятность совершения тех или иных событий не постоянна (Коли, 1979). Считают, что эта модель будет адекватной экспериментальным данным при неравной предрасполо-
женности особей к попаданию в ловушки. Геометрическое распределение является предельным случаем отрицательного биномиального распределения с функцией вероятности Рк = р^к. Абсолютная численность рассчитывается как Ng = п- п)/^ - 1), где п - число меченых особей.
Погрешности оценок определяли простым непараметрическим бутстрепом с числом итераций В = 5000 (Шитиков, 2012; Шитиков, Розенберг, 2013). Построив вариационный ряд бутстрепированных значений средней, методом процентилей находили границы ее доверительных интервалов (рис. 1). Пределы изменчивости Np и Ng оценивали по этим экстремумам, приняв 5 = п-М,. , где М - минимальное или макси-
¡¡т ^ит „
мальное значение бутстрепированной средней. Адекватность моделей исходным дан-
ным определяли по критерию Пирсона (х2), сравнивая эмпирические частоты с теоретическими распределениями. Частоты соответствующих теоретических распределений (Ак) вычисляли по общей формуле Ак = N-9^ где N - абсолютная численность. Вклад отличий нулевых классов в значение критерия принимали равным нулю, допустив, что число непойманных особей а0 действительно соответствует рассчитанному значению А0. Параметры распределений определены как Л = М, р = 1-(s - n)/(s - 1). Использован алгоритм расчетов на языке R (см. скрипт). Число
частотных классов теоретического распределения приравнено к числу классов эмпирической выборки, остатки их «хвостов» слиты с частотой Ак последнего класса. При расчете числа степеней свободы учитывали наличие нулевых классов. Оценка абсолютной численности методом Джолли-Себера выполнена в соответствии с ранее опубликованной методикой (Горбач, 2013), в том числе была реализована процедура оптимизации путем подгонки модельных параметров под биологически разумные значения.
q
¡П -1
{Ч
о
1П ч
1.5 1.5 1.7 1.8 1.9 2.0 2.1
А
о О
гч
о in
о
<=>
о m
м
г
-лпг Jill
Т"
т
baJiann-
I
1.10 1.11 1.12 1.13 1.14 1.15 1.16
Рис. 1. Распределение значений средней арифметической числа отловов M, полученных бутстреп-методом для имаго B. freija 2016 г. (а) и B. aqulonaris 1995 г. (б). А - число выборок (IA = B = 5000), красной линией указано значение M исходной выборки, синей - значение М, скорректированное бут-стрепом, зелеными - границы доверительного интервала M -M
Fig. 1. Distribution of the arithmetic mean value (M) of capture numbers obtained by the bootstrap method for all Freija fritillaries in 2016 (а) and all Cranberry fritillaries in 1995 (б). A is the number of samples (IA = B = 5000), the red line indicates the M value of the original sample, the blue line is the M corrected by bootstrap, the green lines are the confidence interval M . +M
r' ° min max
Результаты
Результаты оценки абсолютной численности и последующего анализа распределений частот отловов бабочек в MRR-эксперименте позволили прийти к следующим заключениям: 1) наиболее высокие значения численности дает геометрическая модель, они в полтора-два раза выше Пуассоновых оценок; 2) значения, полученные методом Пуассона, лучше сходятся с модельными оценками Джолли - Себера; 3) распределение Пуассона более адекватно экспериментальным данным. Вместе с тем имеется и несколько отклонений от общих трендов. Так, оценки абсолютной численности, сделанные методом Джолли - Себера по выборкам 2016 г., более сходны с геометрическими оценками (см. табл. 1); отловы самцов В. aquilonaris в 1995 г. геометрическая модель аппроксими-
рует лучше, чем распределение Пуассона (см. табл. 2).
Общая закономерность в распределении эмпирических частот состоит в недостатке особей, отловленных один раз, и избытке особей, отловленных два раза, по сравнению с теоретическими предсказаниями (рис. 2). Предельный вариант наблюдали у самок В. /ге^а в 1996 г.: отмеченная непропорциональность привела к существенным расхождениям с теоретическими частотами и значимо сказалась на распределении обобщенных частот для всех особей вида в этот сезон. Лишь у упомянутых выше самцов В. aquilonaris эмпирические частоты а1 = 941, а2 = 109 и а3 = 12 смогли попасть в коридор из предсказанных частот Ар-Ад: 931-944, 113-105 и 10-13.
ЭКЗ.
300
250 -
ZOO
150
100 -
50
□ J
2
— 3
0 12 3 4 5 6
-•—•—« i » i
7 8 5 10
Рис. 2. Частота отловов особей для выборок имаго B. freija 1996 (а) и 2016 г (б). k - число отловов особи (частотные классы), 1 - эмпирические частоты а, 2 - частоты геометрического распределения Ag, 3 -
частоты распределения Пуассона Ap Fig. 2. The frequencies of captures for all Freija fritillaries in 1996 (а) and 2016 (б). k - the capture number of the specimens (classes of frequencies), 1 - empirical frequencies a, 2 - frequencies of the geometric
distribution Ag, 3 - the Poisson's frequencies Ap
Обсуждение
Стохастическая модель динамики численности Джолли - Себера предсказывает весьма высокий уровень элиминации изученных группировок бабочек за счет гибели и эмиграции особей (Горбач, 2013). Следовательно, требование равной вероятности отловов нарушается, означая, что оценки абсолютной численности, рассчитанные методом Пуассона, не могут быть вполне адекватными и геометрическая модель, допускающая непостоянство вероятности, должна описывать полученные выборки лучше. Между тем результаты сравнения распределений расходятся с этим предположением. Сходство оценок по модели Джолли - Себера, предусматривающей неустойчивость хода демографических процессов, и модели Пуассона свидетельствует о вполне хорошем описательном потенциале последней, указывая на ее вполне приемлемую толерантность к нарушениям равновероятности отловов в ходе эксперимента. Отдавая предпочтение Пуассоновым оценкам, следует подчеркнуть, что здесь речь идет о пространственно обособленных
группировках с естественным образом ограниченными входящими и выходящими потоками особей. Соответствие частот распределению Пуассона помимо прочего может служить своего рода индикатором «закрытости» подобных популяционных систем.
К оценкам на основе геометрической модели следует относиться осторожно, как к неоправданно высоким, особенно в случаях больших выборок. С ростом населения снижается вероятность повторных отловов, чем меньше средняя, тем больше прирост нулевого класса и выше оценка численности. Например, если судить по рассчитанной численности В. aqui¡onaris (см. табл. 2), то получается, что в 1995 г. средняя плотность населения в скоплениях должна составлять около 500 имаго за сезон на учетную площадку 50 х 50 м (Горбач, 2011), но общее число отловленных на этих площадках особей было существенно меньше: оно не превышало 300 при среднем показателе в 50 экз. Дрейф эмпирических частот в сторону геометрического распределения можно трактовать как проявление разного рода демографических возмущений, в частности роста миграцион-
ных потоков. Однако в упомянутом выше случае с самцами В. aqui¡onaris причина была иная, связанная с условиями получения выборки. Дело в том, что неординарное поведение этого вида позволяло продолжать эксперимент при ухудшении погодных условий, когда бабочки переставали летать, - они оставались на соцветиях сабельника, поэтому, обходя скопления этого растения, можно было собирать, метить и возвращать их обратно. При этом существенно уменьшалась доля самцов, многие из которых, в отличие от самок, скрывались в нижнем ярусе растительности, становясь недоступными. Дней с неустойчивой или пасмурной погодой в сезоне 1995 г. было больше половины, поэтому указанная особенность повлияла на общую выборку, и соответствие геометрическому распределению следует рассматривать здесь как следствие изменчивости вероятности отлова самцов в различных условиях.
Ответ на вопрос, почему значения абсолютной численности В. /к^а в 2016 г., рассчитанные методом Джолли - Себера, сильно отклоняются от Пуассоновых оценок, заключен в нестабильности стохастической модели, прилагаемой к малым выборкам. Участие случайных событий в формировании вариант существенно возрастает, усиливая диспропорции в распределении числа отловов. В результате модельные параметры, используемые для расчета переменных, принимают невозможные с биологической точки зрения значения. Несоответствие их реальности давало заведомо заниженные итоговые оценки, вплоть до того, что рассчитанная абсолютная численность была меньше числа меченых особей. Введение в модель биологически разумных ограничений сделало эти параметры более адекватными, но при этом значения численности стали неоправданно высокими. Понять справедливость этого заключения можно, сравнив равноценные по объему выборки В. /к^а, полученные для самок в 1996 г. и всех имаго в 2016 г. (см. табл. 1). Во втором случае метод Джолли - Себера дал значение абсолютной численности в 1.7 раза больше. Если же сопоставить обобщенные выборки этого вида, а в том и другом сезоне удавалось отлавливать почти всех замеченных бабочек, то видно, что в 2016 г. обилие было по меньшей мере в два раза ниже, чем в 1996 г. Таким образом, остается констатировать, что и здесь Пуассонова модель более реалистично отражает наблюдаемую действительность.
Отмеченная диспропорция между частотами отловов в выборках и их предсказанными вероятностями (см. рис. 2) говорит в пользу «закрытости» исследованных группировок. Обычно уменьшение числа встреч немеченых особей на фоне роста повторных отловов наблюдается при довольно стабильном составе населения, когда особи подолгу остаются на тех участках, где были помечены, а прирост за счет выхода свежих бабочек из куколок и иммиграции несущественен. В качестве иллюстрации можно предложить выборку В. а. В сезоне 1996 г. плотность населения этой перламутровки была сравнительно высокой, но тем не менее удавалось отлавливать всех встреченных бабочек. В результате число впервые отловленных особей с каждым посещением местообитания уменьшалось, а встречаемость меченых росла. У самок оседлость проявилась в большей мере, поскольку они, по сравнению с самцами, не так активны поначалу и начинают расселяться, лишь отложив значительную часть яиц в месте своего отрождения. При росте населения, как в случае с В. aqui¡onaris) возможность отлова каждой замеченной бабочки падает, и, вследствие того, что многие из них остаются немечеными, влияние оседлости на частоту отловов минимизируется. При уменьшении объема выборок, подобно В. /к^а в 2016 г., перестает работать закон больших чисел, и распределение частот оказывается сильно зависимым от случайности каждой встречи.
Заключение
Частоту отловов бабочек в MRR-эксперименте при изучении пространственно обособленных группировок лучше всего описывает распределение Пуассона. Оценки абсолютной численности, рассчитанные на основе этой модели, представляются наиболее адекватными: в целом они довольно хорошо согласуются со значениями, полученными методом Джолли - Себера. В случае малых выборок аппроксимация экспериментальных данных распределением Пуассона, по-видимому, остается единственно возможным способом получения адекватных оценок численности населения. Смещение эмпирических частот в сторону геометрического распределения может быть индикатором динамичности демографических процессов, - рост элиминации неизбежно уменьшает долю повторных отловов, а активное пополнение популяции ведет к увеличению доли впервые отлов-
ленных особей. Рассматривая соответствие распределению Пуассона в качестве меры обособленности группировки, нужно понимать, что локальная оседлость особей может менять характер эмпирического распределения, увеличивая частоту повторных
Библиография
отловов, так, что отличия от теоретического распределения становятся значимыми. Расхождение с Пуассоновой моделью в таком случае будет аргументом в пользу большей «закрытости» изучаемой популяционной системы.
Горбач В. В. Сезонная динамика численности и половой состав популяции перламутровки Boloria aquilonaris (Lepidoptera, Nymphalidae) // Зоологический журнал. 1998. Т. 77. № 5. С. 576581.
Горбач В. В. Пространственная структура популяции и подвижность имаго перламутровки Boloria aquilonaris (Lepidoptera, Nymphalidae) // Экология. 2011. № 4. С. 289-296. doi: 10.1134/ S1067413611040060.
Горбач В. В. Изучение динамики численности методом Джолли - Себера на примере имаго булавоу-сых чешуекрылых (Insecta, Lepidoptera: Hesprioidea et Papilionoidea) // Принципы экологии. 2013. Т. 2. № 2. С. 14-28. doi: 10.15393/j1.art.2013.2601.
Коли Г. Анализ популяций позвоночных . М.: Мир, 1979. 362 с.
Коросов А. В. Специальные методы биометрии . Петрозаводск: Изд-во ПетрГУ, 2007. 364 с.
Шитиков В. К. Использование рандомизации и бутстрепа при обработке результатов экологических наблюдений // Принципы экологии. 2012. № 1. С. 4-24. doi: 10.15393/j1.art.2012.481.
Шитиков В. К., Розенберг Г. С. Рандомизация и бутстреп: статистический анализ в биологии и экологии с использованием R . Тольятти: Кассандра, 2013. 314 с. URL: http://www.ievbras.ru/ ecostat/Kiril/Article/A32/Stare.htm (дата обращения: 21.03.2018).
Krebs C. J. Ecological Methodology. Menlo Park: Addison-Wesley, 1999. 620 p. URL: http://www.zoology. ubc.ca/~krebs/books.html (дата обращения 31.03.2018).
Southwood T. R. E., Henderson P. A. Ecological methods. Oxford: Wiley-Blackwell, 2000. 592 p.
ON THE ESTIMATION OF THE ABSOLUTE POPULATION SIZE FROM THE DISTRIBUTION OF CAPTURES IN THE MRR-EXPERIMENT ON THE EXAMPLE OF THE FRITILLARY BUTTERFLIES OF GENUS BOLORIA {INSECTA, LEPIDOPTERA,
NYMPHALIDAE)
GORBACH Vyacheslav Vasil'evich
Petrozavodsk State University, gorbach@psu.karelia.ru
Key words:
mark-recapture method Poisson's distribution geometric distribution stochastic model of population dynamics by Jolly - Seber Boloria freija Boloria aquilonaris
Summary: The absolute size of an animal population is estimated from the experiments with marked specimens. Populations of butterflies are usually studied by the Jolly - Seber method, which uses the stochastic model as a mathematical basis. However, the initial data that correspond to the minimum requirements for sample size, periodicity and proportionality of captures are not always possible to obtain. In such cases, it remains to estimate the absolute number by one of the known laws of frequency distribution. In this article we consider the techniques of estimation of the absolute population size. The original samples were obtained by the mark-recapture method according to the Jolly-Seber scheme in insular local populations of fritillary butterflies of the genus Boloria. Studies were carried out in the mires Blizkoe and Osokovoe in Kivach Reserve in the summer of 1995, 1996 and 2016. In total, 261 imagos of B. freija and 3628 of B. aquilonaris were marked; the number of recaptures was 274 and 968, respectively. Poisson and geometric distributions were taken as theoretic platforms. The general task was to determine the number of individuals that were not captured based on the proportion of marked specimens in the sample. A script is proposed for calculating the absolute size and its confidence intervals for Poisson and geometric distributions as well as estimating the adequacy of empirical frequencies to theoretic models in the language R. It was established that the Poisson model assuming an equal probability of capture of any specimen is more adequate to the experimental data. The conformity of the empirical frequencies to the Poisson distribution can be considered as a measure of isolation of the population. The Poisson values of the absolute size were similar to the estimates by the Jolly - Seber model, which presumes the existence of inbound and outbound streams of individuals. It calls in question the recommendations for initial choosing certain methods for "closed" and "open" populations. The geometric model is used to describe the processes with non-constant probability of events, i. e. it should better correspond to the "open" Jolly - Seber system. But in practice this model gave highly inflated values, especially when applied to larger samples.
Received on: 10 April 2018 Published on: 25 September 2018
References
Gorbach V. V. Seasonal dynamics and sex ration in a population of the butterfly Boloria aquilonaris (Lepidortera, Nymphalidae), Zoologicheskiy zhurnal. 1998. T. 77. No. 5. P. 576-581.
Gorbach V. V. Spatial distribution and mobility of butterflies in a population of the Cranberry Fritillary Boloria aquilonaris (Lepidoptera, Nymphalidae), Ekologiya. 2011. No. 4. P. 289-296. doi: 10.1134/ S1067413611040060.
Gorbach V. V. The study of population dynamics by the Jolly - Seber method in the butterflies (Insecta, Lepidoptera: Hesprioidea et Papilionoidea), Principy ekologii. 2013. T. 2. No. 2. P. 14-28. doi: 10.15393/j1.art.2013.2601.
Koli G. Analysis of vertebrate populations. M.: Mir, 1979. 362 p.
Korosov A. V. Special methods of biometry. Petrozavodsk: Izd-vo PetrGU, 2007. 364 p.
Krebs C. J. Ecological Methodology. Menlo Park: Addison-Wesley, 1999. 620 p. URL: http://www.zoology.
ubc.ca/~krebs/books.html (data obrascheniya 31.03.2018). Shitikov V. K. Rozenberg G. S. Randomization and bootstrap: statistical analysis in biology and ecology using R. Tol'yatti: Kassandra, 2013. 314 p. URL: http://www.ievbras.ru/ecostat/Kiril/Article/A32/ Stare.htm (data obrascheniya: 21.03.2018). Shitikov V. K. Use of randomization and bootstrap in processing the result of ecological observations,
Principy ekologii. 2012. No. 1. P. 4-24. doi: 10.15393/j1.art.2012.481. Southwood T. R. E., Henderso n P. A. Ecological methods. Oxford: Wiley-Blackwell, 2000. 592 p.
