Введение
Одним из последствий возрастающей антропогенной нагрузки и увеличения воздействия человека на окружающую среду является глобальное изменение климата на планете. Ведущими международными группами ученых разработано несколько моделей и сценариев, которые пытаются предсказать, как в будущем изменятся различные биоклиматические параметры на планете. Несомненно, изменение климата повлечет за собой изменение видового состава и распространения различных групп животных и растений в каждом регионе. В связи с этим прогнозирование и моделирование будущего пространственного распределения видов (особенно индикаторных) является одной из приоритетных задач. Луговые дневные чешуекрылые (Lepidoptera) являются биологическими индикаторами состояния луговых сообществ, в том числе практики землепользования. Например, увеличение интенсификации в использовании ресурсов или оставление земель оказывает прямое влияние на местные виды луговых бабочек [5].
Цель данной работы – создать модель пространственного распределения луговых видов бабочек Брянской области и предсказать изменения пространственного распределения их местообитаний на основе модели будущего климата - CGCM3.1_CCCMA.
Материалы и методы
В качестве исследуемой группы нами были выбраны дневные булавоусые чешуекрылые, характерные для луговых биомов. В данную группу были включены 11 видов, для которых луга являются основным местообитанием на территории Европы [5]. По составу она не однородна. Можно выделить широко распространенные виды (толстоголовка лесовик Ochlodes sylvanus, зорька обыкновенная Anthocharis cardamines, червонец пятнистый Lycaena phlaeas, голубянка Икар Polyommatus icarus, сенница Памфил Coenonympha pamphilus, крупноглазка воловий глаз Maniola jurtina) и специализированные виды (бурокрылка Тагет Erynnis tages, голубянка малая Cupido minimus, голубянка Арион Maculinea arion, голубянка лесная Cyaniris semiargus, голубянка Коридон Polyommatus coridon).
Сбор данных проводился в течение вегетационных периодов 2010- 2012 гг. на трансектах по общепринятой методике учета дневных бабочек [1]. Для построения модели пространственного распределения дневных чешуекрылых на территории Брянской области и выявления факторов, обуславливающих распространение данной группы насекомых, использовался статистический метод – MaximumEntropy (MaxEnt) и пакет математических программ с одноименным названием. Данный подход основан на принципе нахождения распределения вероятности максимальной энтропии с учетом ограничений, налагаемых доступной информацией о наблюдаемом распространении видов и распределении экологических условий в районе исследования [4]. Метод строит вероятностную модель c заданием условной вероятности:
,
где Z(b) – нормирующий множитель.
Для моделирования распространения видов методом максимальной энтропии использовались ГИС-данные: стандартный набор данных Биоклим (Bioclim) [6], который включает 19 параметров, а также климатическая модель CGCM3.1_CCCMA, содержащая предполагаемые значения параметров в 2030, 2060 и 2080 гг. [2]. Расчеты проводились с учетом прогноза SRES A2, который предсказывает существенное увеличение численности населения, интенсивное использование энергии и эмиссию парниковых газов. Это наиболее пессимистичный сценарий из предложенных Межправительственной группой экспертов по изменению климата (IPCC) [3].
Результаты исследования и их обсуждение
В 2010 – 2012 гг. была сформирована база данных по 27 трансектам, на которых проводился учет чешуекрылых. Было отмечено 223 точки, где были зарегистрированы луговые виды дневных бабочек. С использованием данных Биоклим (BioClim) была построена достоверная модель современного пространственного распределения чешуекрылых по территории Брянской области (рис. 1).
Рис. 1. Модель пространственного распределения дневных булавоусых чешуекрылых на территории Брянской области
(красный цвет означает наиболее пригодные для луговых бабочек местообитания, зеленый цвет - наименее пригодные местообитания для луговых бабочек)
Большая часть территории, наиболее подходящей для группы луговых бабочек, располагается в западной части Брянской области. Согласно расчетам на неё приходится 43,34% от общей площади региона. Сюда входят территория левобережья Десны, Брянский лесной массив, а также самые южные районы Брянской области.
Среди использованных в модели биоклиматических факторов были выявлены те, которые оказывают наибольшее влияние на пространственное распределение луговых видов чешуекрылых. К ним относятся: коэффициент вариации осадков (25,9%), годовая сумма осадков (14,8%), сумма осадков самой холодной четверти года (14,7%), сумма осадков самого влажного месяца года (13,8), максимальная температура самого тёплого месяца года (7,7%). В скобках указана доля в процентах вклада фактора в общее влияние на распространение бабочек.
С использованием данных BioClim для 2030, 2060 и 2080 гг. из климатической модели CGCM3.1_CCCMA были построены соответствующие достоверные модели распространения бабочек, приуроченных к луговым местообитаниям. Было установлено, что площадь местообитаний, наиболее подходящих для данной группы, будет изменяться с течением времени, что представлено на рис. 2.
Рис. 2. Изменение площади местообитаний, пригодных для группы луговых дневных булавоусых чешуекрылых с 2012 по 2080 гг.
К 2030 г. площадь подходящих местообитаний сократится на 3,15% по сравнению с 2012 г., к 2060 г. она, наоборот, увеличится и составит 49,40% от общей площади Брянской области. В 2080 г. площадь наиболее благоприятных территорий сократится на 11,3% по сравнению с уровнем 2060 г. и составит 38,10%. Таким образом, к 2080 г. площадь подходящих местообитаний сократится на 4,1% по сравнению с 2012 г. и составит 13295,88 км2.
Также с течением времени изменяется влияние различных биоклиматических переменных на пространственное распределение дневных чешуекрылых. Так, по сравнению с 2012 г. в 2030 г. влияние коэффициента вариации осадков и годовой суммы осадков ослабевает, при этом резко возрастает влияние годовой амплитуды температуры, сумма осадков самой сухой четверти года, средней температуры самой влажной четверти года. К 2060 г. сильно возрастает влияние годовой суммы осадков по сравнению с 2012 и 2030 гг., а также возрастает влияние таких параметров, как стандартное отклонение температур, максимальная температура самого теплого месяца года, средняя температура самой теплой четверти года. В 2060 г. резко возрастает влияние суммы осадков самой тёплой четверти года, но при этом основное влияние оказывает годовая сумма осадков. В 2080 г. еще больше увеличивается влияние суммы осадков самой тёплой четверти года по сравнению с другими годами, а также увеличивается вклад суммы осадков самой сухой четверти года, а влияние годовой суммы осадков сокращается, но при этом остается одним из ведущих факторов в 2080 г.
Выводы
Таким образом, в результате проведенных исследований можно сделать следующие выводы.
1. Изменение климата будет оказывать влияние на распространение дневных булавоусых чешуекрылых на территории Брянской области.
2. Площадь территории наиболее пригодных для обитания бабочек, приуроченных к луговым местообитаниям, к 2080 г. сократится на 4,1% по сравнению с 2012 г.
3. В настоящее время в качестве фактора, оказывающего наибольшее влияние на распространение бабочек, выступает количество осадков.
4. Влияние вклада различных биоклиматических компонент в пространственное распределение луговых бабочек в будущем будет изменяться, но ведущим фактором будет оставаться количество осадков.
5. Поскольку глобальные модели изменения климата имеют ряд допущений, то составленный прогноз изменения пространственного распределения дневных чешуекрылых на территории Брянской области не идеален и будет уточняться в последующей доработке и совершенствовании климатической модели изменения климата.
Резензенты:
Зайцева Е.В., д.б.н, профессор, заведующая кафедрой зоологии и анатомии ФГОБУ ВПО «Брянский государственный университет им. ак. И.Г. Петровского», г. Брянск.
Анищенко Л.Н., д.с-х.н., профессор кафедры экологии и рационального природопользования ФГОБУ ВПО «Брянский государственный университет им. ак. И.Г. Петровского», г. Брянск.