EN
The wintering Great Grey Shrike Lanius excubitor was censused in Poland more than two decades ago (during the winters of 1988/89 and 1989/90). Single censuses were made during each winter on sample plots varying from 4.1 to 35.1 km2 (mean: 16.4 km2). Altogether, data from 404 plots, covering in total more than 6.6 thousands km2 were used for the analysis. Based on this data and environmental information gathered in GIS databases (Corine land cover "CLC1990" database, digital elevation model "GTOP030" dataset) we modelled habitat- and spatial-related variation in shrike distribution. Birds were recorded on 45% of study plots. The mean density was 4.8 individuals/100 km2. In both seasons density in meadows was roughly two times higher than in arable fields. We modelled spatial distribution of wintering birds using combination of GLM and three-dimensional local regression. Models were cross-validated to check their temporal consistency. Repeatable spatial pattern of population distribution allowed to make predictive distribution maps. As a general rule, wintering shrikes avoid regions with severe, continental climate and prefer lowlands dominated by agriculture. On the basis of our models, population of the Great Grey Shrike wintering in Poland during late 1980s can be roughly estimated at about 7.7 thousands of individuals.
PL
Dla wielu gatunków ptaków zimowanie jest krytycznym okresem ich cyklu życiowego, decydującym o późniejszych zmianach liczebności populacji. Srokosz jest gatunkiem osiadłym lub podejmującym niewielkie wędrówki. Z dotychczasowych badań wynika, że wielkość populacji zimującej przede wszystkim zależy od pogody i dostępności pokarmu. Zatem, te właśnie czynniki są uważane za decydujące o wyborze miejsca zimowania. Do tej pory informacje dotyczące zimowania srokosza pochodziły jednak z niewielkich powierzchni badawczych. Brak informacji z dużych obszarów wynikał z faktu, iż ich pozyskanie jest niezwykle pracochłonne i wymaga bardzo dobrze zaplanowanych badań. Wymagania dotyczące odpowiedniej wielkości próby i metodycznie poprawnego sposobu zbierania danych nie zawsze jednak są spełnione, co w szczególności dotyczy prac wykonywanych (w niedalekiej nawet) przeszłości. Błędy takie powodują, że często dane te nie nadają się bezpośrednio do analizy. Z drugiej strony takie historyczne źródła informacji są niezwykle cenne, gdyż umożliwiają porównania i śledzenie procesów zachodzących na przestrzeni wielu lat. W przypadku srokosza ma to szczególne znaczenie, gdyż gatunek ten wykazuje spadki liczebności na terenie całej zachodniej Europy Stabilne lub nawet wzrastające liczebnie populacje utrzymują się w krajach Europu Środkowej, co pozwala przypuszczać, iż przyczyn zmian liczebności należy dopatrywać się w zmianach zachodzących w krajobrazie rolniczym. Praca niniejsza jest próbą zastosowania stosunkowo nowych technologii (teledetekcja satelitarna, modelowanie predyktywne) do analizy danych zebranych przed ponad 20 laty. Starano się w niej pokazać, że nawet niezbyt dobrej jakości dane historyczne mogą posłużyć do modelowania wybiórczości środowiskowej i zmienności przestrzennej rozmieszenia populacji. Zaprezentowano predyktywne mapy rozmieszczenia, a także oszacowania dotyczące zagęszczeń i wielkości zimującej populacji srokosza w Polsce na przełomie lat 80. i 90. XX wieku. Dane terenowe zostały zebrane podczas dwóch zim (1988/1989 i 1989/1990). Łącznie policzono ptaki na, odpowiednio, 192 i 212 powierzchniach pokrywających obszar 3073.5 i 3551.0 km2. Średnia wielkość powierzchni wynosiła 16.4 km2 (zakres: 4.1-35.1). Powierzchnie były zlokalizowane w całym kraju (Fig. 1). Obserwatorzy wybierali je samodzielnie w krajobrazie rolniczym i wykonywali pojedyncze liczenia między połową grudnia a połową stycznia. Zastosowano dane teledetekcyjne o pokryciu terenu pochodzące z bazy Corine "CLC1990" oraz numeryczny model terenu z bazy "GTOP030". Wszystkie analizy były wykonywane w rozdzielczości 4 km2. Dla każdej powierzchni dane były: liczba ptaków oraz powierzchniowe udziały klas pokrycia terenu. Zagęszczenie na każdej powierzchni może być więc opisane jako średnia ważona zagęszczeń w poszczególnych środowiskach (równanie 1). Estymacja parametrów modelu (równanie 2) wykonana została metodą największej wiarygodności przy wykorzystaniu logarytmicznej funkcji wiążącej i założeniu, że rozkład błędów jest rozkładem Poissona. Przedziały ufności zostały obliczone metodą bootstrap. Następnie, dla każdej powierzchni obliczono zagęszczenie resztowe (RD), czyli różnicę między zagęszczeniem obserwowanym a oczekiwanym. Wielkość ta mówi, o ile zagęszczenie na danej powierzchni badawczej różni się od zagęszczenia oczekiwanego jedynie na podstawie struktury środowisk. Przestrzenna zmienność RD była modelowana poprzez dopasowanie do trójwymiarowych współrzędnych (współrzędne geodezyjne i wysokość n.p.m.) modelu regresji lokalnie ważonej (loess). Zdolność modeli do generalizacji była oceniana na podstawie testów krzyżowych: dla każdego sezonu był budowany model i jego skuteczność była sprawdzana na znanych wartościach z innego sezonu. Według tego samego algorytmu wyprodukowane zostały mapy rozmieszczenia. Dla każdego sezonu oddzielnie, w kwadratach 4x4 km oszacowano RD na postawie współrzędnych geodezyjnych i wysokości. Potem, predykcje z modelu przestrzennego zostały dodane do predykcji z modelu środowiskowego. Tak powstałe predykcje z obu sezonów zostały uśrednione, czego efektem jest mapa zaprezentowana na Fig. 3. Wszystkie obliczenia wykonano w środowisku R 2.8.1. Srokosz był obserwowany na 45% powierzchni ze średnim zagęszczeniem 4.8 (95% CI: 4.1-5.6) os./100 km2. Średnie zagęszczenia nie różniły się między latami badań, natomiast stwierdzono istotne różnice między środowiskami (Tab. 1). Wielkość zimującej populacji została oszacowana na 7670 osobników (95% CI: 5700-10600, Tab. 2). Stwierdzono istotne interakcje między położeniem geograficznym a wysokością (Fig. 2). Modele dopasowane niezależnie do danych z dwóch różnych sezonów wykazują podobne wzorce, co wskazuje na spójność predykcji. Także mapy predyktywne dla obu sezonów są podobne. Generalnie, srokosz unika zimą miejsc o surowym klimacie kontynentalnym (NE Polska oraz wyżyny, pogórza i góry). Preferuje łagodny klimat morski i nisko położone równiny, zwłaszcza doliny rzeczne (Fig. 3). Informacje na temat wybiórczości miejsc zimowania srokosza, a zwłaszcza wiedza o rozmieszczeniu odpowiednich środowisk w krajobrazie rolniczym, mogą pozwolić na opracowanie skutecznych metod ochrony tego gatunku w większej skali przestrzennej.