СТАТИСТИЧЕСКАЯ ОБРАБОТКА РЕЗУЛЬТАТОВ МОНИТОРИНГА ОРНИТОКОМПЛЕКСОВ НА ТЕРРИТОРИИ БОТИЕВСКОГО ВЕТРОПАРКА (ЗАПОРОЖСКАЯ ОБЛАСТЬ)

Научная статья
DOI:
https://doi.org/10.60797/IRJ.2025.159.78
Выпуск: № 9 (159), 2025
Предложена:
10.04.2025
Принята:
23.07.2025
Опубликована:
17.09.2025
370
3
XML
PDF

Аннотация

В статье рассмотрено применение методов статистической обработки данных для анализа динамики изменения численности птиц на территории Ботиевской ветровой электростанции и в буферных зонах в период её эксплуатации с 2013 г. по 2021 г. включительно. На первом этапе анализа проверялась возможность появления ошибочных измерений или аномальных величин, которые выпадали из общих статистических закономерностей, для чего использовался метода Ирвина. В результате проверки часть данных была изъята из дальнейшего рассмотрения. На втором этапе проверялась возможность корреляционной зависимости скорректированной численности птиц на территории ВЭС от времен её эксплуатации в различные сезоны. Значимость коэффициентов корреляции проверялась с помощью Т-критерия Стьюдента. Вычисленное значение критерия для зимнего сезона равно -0,062, что свидетельствует об отсутствии зависимости количества птиц от времени эксплуатации. Критерии Стьюдента для весеннего, гнездового, осеннего сезонов и общей численности птиц за весь год находились в пределах от 2,77 до 5,12, что больше критической величины 2,45 для 5%-го уровня значимости. Следовательно, в этом случае можно говорить о зависимости, которая была представлена в виде уравнений регрессии: (весенний сезон), (гнездовой сезон), (осенний сезон), (общая численность).

1. Введение

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

,
. Ботиевская ветровая электростанция (ВЭС) расположена на побережье Азовского моря, поэтому изучение её влияния на миграционные процессы представляет большой интерес. Результаты мониторинга Ботиевского ветрового парка и прилегающей территории получены группой сотрудников Научно-исследовательского центра «Биоразнообразие» Мелитопольского государственного университета имени Б. Хмельницкого под руководством кандидата биологических наук В. Д. Сиохина в 2013–2021 годах и опубликованы в отчёте
. В статье Сидоренко В.С., Горлова П.И. и Сиохина В.Д. 
приведены данные об изменении численности птиц во время весенней миграции на территории ВЭС. Настоящая работа посвящена статистическому анализу динамики миграционных процессов в зимний, весенний, гнездовой и осенний периоды миграции на территории Ботиевской ВЭС и буферных зонах.

2. Структурные характеристики территории Ботиевской ВЭС и динамика птиц в различные сезоны в 2013-2021 годах

Раздел составлен по данным работы

Ботиевская ВЭС расположена на побережье Азовского моря в пределах Приазовского района Запорожской области (рис. 1). В антропогенном аспекте площадка ВЭС с прилегающей территорией и буферными зонами занята только сельскохозяйственными угодьями и лесополосами. Лесонасаждения представлены лесополосами, которые отличается составом древесных насаждений, шириной и кустарниковым ярусом.

Схема расположения Ботиевской ВЭС

Рисунок 1 - Схема расположения Ботиевской ВЭС

Учёт сезонных комплексов птиц проводился маршрутным методом. Для получения информации о динамике птичьих перемещений проводились наблюдения в отдельных мониторинговых точках и на выделенных участках, что позволяло определять угрозы столкновения птиц с ветровыми установками. Численность учтённых птиц на территории ВЭС и буферных зонах в зимний, весенний, гнездовой и осенний периоды в 2013–2021 годах представлялась в виде таблиц. Пример информации о видовом разнообразии и численности птиц в период весенней миграции в 2013–2021 годах приведён в таблице 1.

Таблица 1 - Сравнительная характеристика видового разнообразия и численности птиц в период весенней миграции в 2013-2021 годах

Вид

2013

2014

2015-2019

2020

2021

2013-2021

1

Поганка черношейная

 

 

25

25

 

50

2

Поганка большая

-

101

698

4

11

814

3

Баклан большой

87

65

11735

169

617

12673

4

Цапля большая белая

11

-

-

4

-

15

...

122

Коноплянка

-

-

80

-

11

91

123

Просянка

5

-

31

19

8

63

124

Овсянка обыкновенная

-

-

8

-

-

8

125

Овсянка садовая

-

-

1

-

-

1

Всего

Количество видов

28

38

89

52

41

92

Численность птиц

7921

6773

32667

2402

2612

52575

Обработка полученных данных показала, что орнитокомплекс на территории Ботиевской ВЭС в период 2013-2021 гг. состоял из 92 видов птиц. Наименьшее количество видов зафиксировано в 2018 г. (23 вида). Наибольшее количество видов зафиксировано в 2020 г. (52 вида).

Среди представленных видов по численности преобладает Турухтан (Philomachuspugnax), Бережливая ворона (Corvusfrugilegus) и Пелагаз (Tadornatadorna). Наименьшее количество наблюдается у Коровайки (Plegadisfalcinellus), Аиста белого (Ciconiaciconia), Утки серой (Anasstrepera), Ястреба малого, Веретенника малого (Limosalapponica), Сыча домового (Athenenoctua), Дрозда певчего (Turdusphilomelos)  и Овсянки садовой — по одной особи.

В течение всего весеннего периода наблюдений зафиксировано 52575 особей 92 видов. Численность птиц в различные сезоны представлена в таблице 2.

Таблица 2 - Численность птиц, учтённых в 2013-2021 годах на территории Ботиевской станции и буферных зонах в различные периоды

Численность птиц

Количество видов

Зимний период

8 541

49

Весенняя миграция

52 575

92

Гнездовой период

18 316

108

Осенняя миграция

50 552

87

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

В таблице 3 приведена численность зарегистрированных птиц по данным отчёта

.

Таблица 3 - Количество птиц на территории Ботиевской станции и буферных зонах

Год

2013

2014

2015

2016

2017

2018

2019

2020

2021

Зимний сезон

741

511

624

2644

619

365

1622

868

547

Весенний сезон

7921

6773

16297

9526

2593

2480

1971

2402

2612

Гнездовой сезон

3021

5364

2714

2301

1145

792

238

1206

1544

Осенний сезон

6330

7845

1255

16097

7869

2777

4162

2472

1645

Всего

18013

20493

20890

30568

12226

6414

7993

6948

6348

3. Статистическая обработка изменения численности птиц в период эксплуатации ВЭС в различные сезоны от 2013 г. до 2021 г.

Статистический анализ данных, представленных в таблице 3, проводился в соответствии с рекомендациями, изложенными в работе

. На первом этапе с использованием метода Ирвина проверялась возможность появления ошибочных измерений или аномальных величин, которые выпадают из общих статистических закономерностей полученных результатов. Метод основан на определении коэффициентов λi=|хii-1|/sх для каждой строки из таблицы 3 в отдельности, где хi  — численность птиц в i-м году (i=0 относится к 2013 г., i=1 — к 2014 г. и т.д., i=9 — к 2021 г.), sх — среднеквадратическое отклонение, вычисляемое по формуле

(1)

Величина хср определяет среднее значение численности птиц в анализируемом сезоне. Полученные значения λi сравнивались с критическим параметром λα. Если λi>λα , то больший из двух элементов с номерами i-1 и i в статистическом ряде отбрасывался из анализа. Параметр λα для объёма выборки, равного 9, и 5%-го уровня значимости приблизительно равен 1,9. Следует отметить, что отбракованный элемент может не являться результатом ошибки измерения, а отражать влияние известного или неизвестного фактора. Но здесь математика бессильна. Решение по отбраковке принимает исследователь. Вычисленные коэффициенты λi приведены в таблице 4.

Таблица 4 - Параметр Ирвина λi в случае учёта птиц на территории Ботиевской ВЭС

i

2

3

4

5

6

7

8

9

Год

2014

2015

2016

2017

2018

2019

2020

2021

Зимний сезон

0,31

0,15

2,76

2,77

0,33

1,66

1,0

0,42

Миграция весной 

0,24

1,97

1,4

1,43

0,05

0,23

0,22

0,08

Гнездовой сезон

1,52

1,72

0,27

0,75

0,98

2,61

2,33

0,66

Миграция осенью

0,32

1,41

3,17

1,76

1,57

0,64

0,72

0,34

Общая численность

0,29

0,05

1,13

2,15

1,41

0,41

0,13

0,08

Из таблицы 4 видно, что результаты ряда измерений выпадают из общих статистических закономерностей и подлежат отбраковке. К ним относятся следующие данные:

Численность птиц в зимний период миграции в 2016 г. (2644 птиц), параметр λ4=2,76.

Численность птиц в весенний период миграции в 2015 г. (2387 птиц), параметр λ3=1,97.

Численность птиц в гнездовой период миграции в 2019 г. (238 птиц), параметр λ7=2,61.

Численность птиц в осенний период миграции в 2016 г. (16097 птиц), параметр λ4=3,17.

Общая численность в 2016 г. (30568 птиц), параметр λ5=2,15.

Скорректированные данные по численности птиц после отбраковки возможных ошибочных измерений или аномальных величин представлены в таблице 5.

Таблица 5 - Скорректированная численность птиц на территории Ботиевской станции и буферных зонах

Год

2013

2014

2015

2016

2017

2018

2019

2020

2021

Зимний сезон

741

511

624

 

619

365

1622

868

547

Весенний сезон

7921

6773

 

9526

2593

2480

1971

2402

2612

Гнездовой сезон

3021

5364

2714

2301

1145

792

 

1206

1544

Осенний сезон

6330

7845

1255

 

7869

2777

4162

2472

1645

Всего

18013

20493

20890

 

12226

6414

7993

6948

6348

На втором этапе анализировалась возможность корреляционной зависимости численности птиц xот года учёта tв период 2013–2021 г. Наличие линейной связи между ними определялось по значению коэффициента корреляции

(2)

где sx=√sx2, st=√st2 — среднеквадратические отклонения.

Значимость коэффициента (2) проверялась с использованием Т-критерия Стьюдента 

для степеней свободы:

(3)

Если значение Т больше критической величины Tkr при выбранном уровне значимости q, то между параметрами количества птиц х и временем эксплуатации t существует линейная связь, которую можно представить в виде линейного уравнения регрессии

(4)

где a=xср-rtxtasx/st,   b=rtxsx/st.

В противном случае, при Т<Tkr подобная связь отсутствует. Вычисленные коэффициенты корреляции для выборок, представленных в таблице 5, соответствующие значения Т-критерия Стьюдента для 5%-го уровня значимости для n-2=8-2=6 степеней свободы и критические значения Tkr приведены в таблице 6.

Таблица 6 - Параметры корреляционной связи

Параметры

rtx

T

Tkr

a

b

Зимний сезон

-0,062

0,163

2,45

570,09

32,59

Весенний сезон

-0,847

4,193

2,45

8935,53

-838,24

Гнездовой сезон

-0,815

3,736

2,45

4040,82

-374,73

Осенний сезон

-0,742

2,766

2,45

6954,79

-519,11

Общая численность

-0,889

5,123

2,45

22782,32

-2022,77

Примечание: rtx – коэффициенты корреляции; T-критерии Стьюдента; Tkr – критические значения критериев Стьюдента; b – коэффициенты уравнения регрессии

Из таблицы 6 видно, что во всех случаях, за исключением зимнего сезона, коэффициент Стьюдента (3) больше критического значения. Следовательно, корреляционная связь существует для большинства сезонов, которую в соответствии с формулой (4) можно представить в виде уравнений регрессии.

Весенний сезон:

(5)

Гнездовой сезон:

(6)

Осенний сезон:

(7)

Общая численность:

(8)

Корреляционные зависимости численности птиц от времени эксплуатации ВЭС, рассчитанные по формулам (5)–(8), приведены на рисунках 2–5. Наиболее сильно выраженные зависимости относятся к весеннему сезону, осеннему сезону и к общей численности птиц.

Корреляционная зависимость численности птиц от времени эксплуатации ВЭС для весеннего сезона, вычисленная по формуле (5)

Рисунок 2 - Корреляционная зависимость численности птиц от времени эксплуатации ВЭС для весеннего сезона, вычисленная по формуле (5)

Корреляционная зависимость численности птиц от времени эксплуатации ВЭС для гнездового сезона, вычисленная по формуле (6)

Рисунок 3 - Корреляционная зависимость численности птиц от времени эксплуатации ВЭС для гнездового сезона, вычисленная по формуле (6)

Корреляционная зависимость численности птиц от времени эксплуатации ВЭС для осеннего сезона, вычисленная по формуле (7)

Рисунок 4 - Корреляционная зависимость численности птиц от времени эксплуатации ВЭС для осеннего сезона, вычисленная по формуле (7)

Корреляционная зависимость общей численности птиц от времени эксплуатации ВЭС, вычисленная по формуле (8)

Рисунок 5 - Корреляционная зависимость общей численности птиц от времени эксплуатации ВЭС, вычисленная по формуле (8)

Достоверность полученных зависимостей подтверждается значениями коэффициентов детерминации rtx2 и результатом анализа адекватности уравнений регрессии (5-8) с использованием критерия Фишера для 5%-го уровня значимости, которые приведены в таблице 7. Критерий Фишера вычислялся по формуле
(9)

где — число параметров в уравнении регрессии, в данном случае равное 2, n=8 — количество измерений птиц, которое независимо от сезона было одинаковым и согласно таблице 5 равным 8.

Таблица 7 - Значения коэффициента корреляции rtx, коэффициента детерминации rtx2, критерия Фишера F и критического значения критерия Фишера Fкр

Параметры

rtx

rtx2

F

Fкр

Весенний сезон

-0,847

0,723

15,66

5,59

Гнездовой сезон

-0,815

0,664

11,86

5,59

Осенний сезон

-0,742

0,551

7,36

5,59

Общая численность

-0,889

0,790

22,57

5,59

Примечание: по источнику [6]

4. Заключение

Поскольку временная зависимость количества птиц, учтённых в зимний сезон отсутствует, то в этом случае можно говорить только о численности птиц в рамках статистического разброса относительно среднего значения хср=737. Величина доверительного интервала определялась в предположении, что выборочные данные подчиняются нормальному закону распределения. Разброс численности птиц Δ относительно среднего значения оценивался по формуле Δ=0,67sx для 50%-го доверительного интервала и по формуле Δ=1,64sx для 90%-го интервала

. Минимальное хmin и максимальное хmax значения интервалов вычислялись по формулам хminср-Δ, хmax=хср+Δ. Нижний и верхний пределы для доверительных интервалов оказались равными: 495 и 979 для 50%-го интервала, 272÷1200 для 90%-го интервала. 

Метрика статьи

Просмотров:370
Скачиваний:3
Просмотры
Всего:
Просмотров:370