Сравнение пространственно-временной динамики промысловых видов животных, обитающих на территориях Среднего Приамурья и Свердловской области

Большинство популяций животных, включая охотничьи виды, очень часто демонстрирует многолетние колебания, вызываемые эндо- и экзогенными факторами. Как правило, на локальном уровне циклические изменения численности наблюдаются вследствие воздействия плотностно-зависимой регуляции и особенностей жизненного цикла особей [2, 5, 6, 9, 11]. В случае взаимодействия нескольких группировок одного вида или разных видов, обитающих в одном ареале (например, миг-рационно связанные субпопуляции, система хищник-жертва, общий хищник и т.д.), режимы динамики численности существенно усложняются, могут возникать такие явления как синхронизация колебаний и колебания с запаздыванием [8, 9, 11, 13, 14]. Так, резкий рост численности популяции, обитающей в некоторой части ареала, может приводить к синхронному росту численности вида по всему ареалу, вызванному миграцией на сопредельные территории [7, 10, 12]. Кроме того, известны случаи, когда наблюдается сдвиг фазы колебаний (колебания с запаздыванием) численностей разных видов в одном местообитании. Так, в работе [13] показано, что в рамках одного региона (Юкон, Канада) плотности молоди тонкорогого барана (Ovis dalli) и зайца-беляка (Lepus americanus) имеют колебания с запаздыванием в 1–2 года, данный эффект авторы объясняют косвенным воздействием общих хищников. Действительно, в сообществе, состоящем из трех популяций разных видов, взаимодействующих по принципу хищник – две жертвы, возможны колебания численностей с запаздыванием между жертвами, в силу того, что хищник переключается с одной жертвы на другую.

Также в природе отмечаются связанные колебания численности в географически разобщенных популяциях, когда воздействие некоторых факторов синхронизирует изменения численности популяций, имеющих схожую систему регуляции [2, 3, 12]. Наблюдаемая синхронность на обширных территориях обычно интерпретируется как следствие воздействия внешних синхронизирующих эффектов (погода, корма, хищники и т.д.) [11, 13].

Данная статья посвящена анализу и сопоставлению динамики численности охотничье-промыс-ловых видов млекопитающих, обитающих на территориях Еврейской автономной и Свердловской областей (ЕАО и СО). Целью работы является анализ колебаний численности популяций и поиск связей между ними как на региональном, так и межрегиональном уровнях. В рамках заявленной цели решаются следующие задачи: 1) анализ наличия цикличности в динамике численности популяций отдельных видов для каждой из территорий (автокорреляционный анализ); 2) выявление синхронных колебаний численности разных видов, обитающих в одном регионе (территории Еврейской автономной и Свердловской областей) (корреляционный анализ); 3) анализ синхронизирован- ности колебаний численности популяций одного вида, обитающих в географически разобщенных ареалах (кросс-корреляционный анализ).

Материал и методы

В работе использованы временные ряды оценок общей численности популяций за период с 1988 по 2013 гг., полученные методом зимнего маршрутного учета (ЗМУ), проведенного на территориях Еврейской автономной и Свердловской областей. Данные для Свердловской области предоставлены Департаментом по охране, контролю и регулированию использования животного мира Свердловской области, для Еврейской автономной области – управлением по охране и использованию объектов животного мира правительства Еврейской автономной области. Сравнивали и анализировали динамику численности следующих видов животных: рысь ( Felis lynx Linnaeus, 1758), обыкновенная лисица ( Vulpes vulpes) , колонок ( Mustela sibirica ), лось ( Alces alces ), дикий кабан ( Sus scrofa ), сибирская косуля ( Capreolus pygargus Pallas, 1871), обыкновенная белка ( Sciurus vulgaris ), заяц-беляк ( Lepus timidus ), соболь ( Martes zibellina ).

Для решения первой задачи был использован автокорреляционный анализ. Для анализа колебаний численности популяций и поиска связей между ними как на территориальном, так и межтерриториальном уровнях использованы автокорреляционный, корреляционный и кросс-корреляционный анализы.

Выявление циклических колебаний численности в локальных популяциях, обитающих на территориях Еврейской автономной и Свердловской областей

Как правило, динамика численности большинства видов животных стремится к равновесию, со- вершая при этом периодические или апериодические флуктуации с различной амплитудой. В математической биологии широко известно понятие популяционных циклов [4, 9, 11], для анализа которых удобно использовать автокорреляционную функцию, позволяющую выявлять циклическую составляющую временных рядов [1]. Автокорреляционная функция временного ряда представляет собой последовательность коэффициентов автокорреляции 1, 2 и т.д. порядков. Найденные автокорреляционные функции для большинства видов, обитающих на рассматриваемых территориях, говорят о наличии тенденции, т.е. существует значимая связь между численностью популяции данного вида в этом году от ее численности в предыдущем. Значимые коэффициенты автокорреляции первого порядка на территории ЕАО имеют следующие виды: кабан, косуля, соболь, кабарга, изюбр, лось, бурый медведь, выдра. На территории Свердловской области – косуля, куница, заяц-беляк, колонок, лисица. На рис. 1 для некоторых видов, демонстрирующих наличие тенденции, представлены коррелограммы – графики зависимости значений коэффициентов автокорреляции от величины лага (порядка коэффициента автокорреляции).

Особый интерес представляют случаи, когда наиболее высоким оказывается коэффициент автокорреляции порядка m >  1 . Это означает, что ряд содержит циклические колебания с периодом m . Из многообразия рассмотренных временных рядов периодические флуктуации наблюдаются в популяциях белки (Свердловская область), косули (ЕАО), соболя (ЕАО) и куницы (Свердловская область) (рис. 1). На рис. 1 видно, что в популяции белки, обитающей на территории Свердловской области, существует противофазный восьмилетний

Рис. 1. Коррелограммы временных рядов: а–б с тенденцией и в–е с циклическими колебаниями

цикл, скорее всего связанный с миграционной активностью данного вида. Численность куницы испытывает трехлетние колебания, которые могут быть объяснены трехлетним циклом численности многих видов мышевидных грызунов. Автокорреляционный анализ выявил статистически значимую цикличность колебаний численности соболя, которые могут быть объяснены периодическими изменениями кормовой базы, а именно колебаниями численности мышевидных грызунов и урожая орехов. Наиболее интересная ситуация связана с популяцией косули в ЕАО, как видно, значимые коэффициенты автокорреляции наблюдаются для лага 2, 3 и 4 года. Однако их значения меньше коэффициента 1 порядка, поэтому можно принять гипотезу о том, что численность косули, обитающей на территории ЕАО, зависит от ее численности в предыдущем году. Для всех остальных подвидов коэффициенты автокорреляции не являются статистически значимыми, то есть временные ряды не содержат линейной тенденции и, следовательно, для исследования циклической составляющей требуются дополнительные исследования.

Синхронные изменения численности популяций, обитающих в рамках одного ареала (территории Еврейской автономной и Свердловской областей)

Для более подробного анализа колебаний численности разных видов, обитающих в одном регионе, и выявления связей между ними использовался параметрический корреляционный анализ. Такой подход позволяет исследовать синхронность колебаний численности популяций в рамках одного ареала. Найденные линейные коэффициенты корреляции для всевозможных сочетаний пар временных рядов общей численности промысловых животных, обитающих внутри территорий Свердловской (10 видов) и Еврейской автономной областей (14 видов) (общих видов 8), представлены на рис. 2.

На основе имеющихся оценок (рис. 2) можно сделать вывод о том, что для большинства видов, обитающих на территории ЕАО, характерны синфазные (синхронные) колебания. Действительно, для экологических систем достаточно часто наблюдается локальная когерентность популяционной динамики различных видов, существующих на одной территории, объясняемая воздействием одних и тех же внешних факторов (в частности, погодных). Популяции лисицы, колонка, куницы и горностая в Свердловской области испытывают противофазные колебания. Наиболее ярко выделяются синфазные колебания численности косули и кабана, что наблюдается на обеих территориях.

Дополнительно проводился расчет коэффициентов корреляции между временными рядами численностей одинаковых видов животных, обитающих на территориях Еврейской автономной и Свердловской областей. Как оказалось, синхронные изменения численности для этих регионов не характерны. Скорее всего, это связано с их отдаленностью друг от друга (между центрами Свердловской и Еврейской автономной областей – более 4500 км).

Рис. 2. Значимые корреляционные связи между численностями промысловых животных, обитающих внутри локальной территории (уровень значимости 0,01)

Анализ пространственных волн численности популяций одного вида, обитающих в географически разобщенных ареалах

(территории Еврейской автономной и Свердловской областей)

В природе неоднократно отмечалась и отмечается крупномасштабная пространственно-временная когерентность популяционной динамики видов и экологических процессов на территориях, отстоящих друг от друга на сотни километров, когда прямое взаимодействие наблюдаемых видов друг на друга исключено. Следует заметить, что впервые на существование такого феномена глобальной когерентности экологических процессов указывал А.Л. Чижевский еще в начале прошлого века. В рамках данного раздела анализ колебаний численности между популяциями одного вида, обитающими на территориях Свердловской и Еврейской автономной областей, проведен на основе кросс-корреляционного анализа. Наиболее яркие результаты, демонстрирующие пространственновременную когерентность динамики видов, заселяющих географически разобщенные ареалы, представлены на рис. 3.

Как видно, у шести видов из девяти удалось выявить взаимосвязанные колебания с запаздыванием (оценки линейного коэффициента кросс-корреляции, уровень значимости 0,05). Обращают на себя внимание синхронные колебания численности популяций лисицы и колонка в направлении ЕАО-СО, при этом сдвиг фазы составляет один год.

Примечательность данного результата заключается в том, что эти виды имеют схожую кормовую базу; они преимущественно питаются мышевидными грызунами. Численности популяций рыси также демонстрируют синхронные колебания: изменения, наблюдаемые в ЕАО, запаздывают на три года относительно флуктуации численности рыси, обитающей в Свердловской области, т.е. если в данном году в Свердловской области наблюдается рост численности, то через 3 года с определенной вероятностью следует ожидать некоторый рост численности популяции рыси в ЕАО. В соответствии с полученными оценками коэффициентов кросс-корреляции в популяции белки в направлении ЕАО – Свердловская область существует запаздывание колебаний с фазой задержки 5 лет. Это может быть связано с миграционной активностью вида в поисках обильной кормовой базы, когда особи мигрируют на сопредельные территории. Помимо синхронных колебаний с запаздыванием наблюдаются противофазные колебания между численностями популяций зайца и кабана, обитающими на территориях Еврейской автономной и Свердловской областей. В целом полученные оценки позволяют говорить о существовании пространственных волн численности популяций одного вида.

Заключение

Показано, что большинство популяций охотни-чье-промысловых видов животных, обитающих на территориях Свердловской области и Среднего год, соответствующий численности в ЕАО 1995 1997 1999 2001 2003 2005

год, соответствующий численности в ЕАО

1994     1996     1998    2000

в в

в год, соответствующий численности в СО

4.5

О

год, соответствующий численности в СО

год, соответствующий численности в СО год, соответствующий численности в СО

1989   1991   1993   1995   1997

год, соответствующий численности в ЕАО 1994 1996 1998 2000 2002 2004

К в к

2002   2004   2006 2008   2010

год, соответствующий численности в ЕАО

3 1995 1997 1999 2001 2003 2005

год, соответствующий численности в СО

год, соответствующий численности в ЕАО 1981 1983 1985 1987   1989 1991

1996 1998 2000 2002 2004 2006 год, соответствующий численности в СО

Рис. 3. Пространственные волны численности популяций одного вида, обитающих в географически разобщенных ареалах; численность представлена в виде натуральных логарифмов

Приамурья, имеют схожие тенденции изменения численности, совершая при этом периодические или апериодические флуктуации с различной амплитудой. Для большинства видов животных характерно наличие инерциальной тенденции, т.е. численность популяции в данном году существенно зависит от ее состояния в предыдущем году (Свердловская область: косуля, куница, заяц-беляк, колонок, лисица; Среднее Приамурье: кабан, косуля, соболь, кабарга, изюбрь, лось, бурый медведь, выдра). Периодические флуктуации наблюдаются в популяциях белки (Свердловская область), соболя (ЕАО) и куницы (Свердловская область).

В ходе исследования показано, что колебания численности большинства видов оказались в значительной мере синхронизированными. Согласованность прослеживается как в колебаниях численности одного вида в разных регионах, так и в динамике численности разных видов охотничьих животных в одном регионе. Это позволяет с высокой вероятность полагать, что на территориях Среднего Приамурья и Свердловской области происходят закономерные циклические колебания численности многих животных, синхронно охватывающие значительные части ареалов.

Исследование выполнено при частичной финансовой поддержке Комплексной программы фундаментальных исследований «Дальний Восток» и проекта 12-С-4-1012 УрО РАН «Пространственно-временная когерентность экологических процессов на территории Евразии».