PLONOWANIE PSZENŻYTA OZIMEGO
A WARUNKI METEOROLOGICZNE OKRESU JESIENNO-ZIMOWEGO W REJONIE WARMIŃSKO-MAZURSKIM
Jan Grabowski, Katarzyna Pożarska, Krystyna Grabowska
Katedra Gospodarki Wodnej, Klimatologii i Kształtowania Środowiska Uniwersytet Warmińsko-Mazurski w Olsztynie
Pl. Łódzki 1, 10-727 Olsztyn e-mail: jangrab@uwm.edu.pl
S t r e s z c z e n i e: Wielkość i jakość plonu zbóż w znacznym stopniu jest uzależniona od warun-ków pogodowych okresu wegetacji, jak i kształtowania się wskaźniwarun-ków meteorologicznych w okre-sie jeokre-sienno-zimowym. Warunki te, zwłaszcza w Polsce, wpływają na dużą zmienność plonowania wynoszącą średnio 30-40%. Celem niniejszej pracy była ocena wpływu wskaźników meteorolo-gicznych w okresie jesienno-zimowym na plonowanie pszenżyta ozimego w latach 1993-2013 w rejonie Warmińsko-Mazurskim. Pszenżyto jest jedną z głównych roślin zbożowych uprawianych w Polsce. Jego udział w strukturze zasiewów wynosi ok. 10%, a średnie plony wahają się w
grani-cach 3-4 t·ha–1. W badaniach uwzględniono następujące wskaźniki meteorologiczne: liczbę dni od
siewu do okresu zimowego, średnie dobowe temperatury od siewu do początku zimy, sumę opadów od siewu do okresu zimowego, długość okresu zimowego, sumę średnich dobowych temperatur zimy, sumę opadów okresu zimowego, długość zalegania pokrywy śnieżnej, liczbę dni z odwilżami atmosferycznymi. Wpływ badanych wskaźników meteorologicznych okresu jesienno-zimowego na plonowanie pszenżyta ozimego oceniono statystycznie, wykorzystując analizę korelacji i regresji dla jednej lub wielu zmiennych niezależnych. W analizowanym 21-letnim okresie badań plony
pszen-żyta ozimego wynosiły średnio 5,97 t·ha–1. Najmniejsze były w 1998 roku – 3,78 t·ha–1, największe
8,83 t·ha–1 w 2004 r. Z przeprowadzonych analiz wynika, iż wielkość plonu pszenżyta ozimego
istotnie zależała tylko od liczby dni z odwilżami atmosferycznymi. O łącznym wpływie na plon decydowały głównie: liczba dni z odwilżami atmosferycznymi w okresie zimy, wartości średniej dobowej temperatury od siewu do okresu zimowego i sumy opadów okresu zimowego.
WSTĘP
Plonowanie roślin w warunkach klimatycznych Polski charakteryzuje się dużą zmiennością czasową i przestrzenną. Związane to jest między innymi z warunka-mi glebowywarunka-mi, agrotechnicznywarunka-mi oraz dużym zróżnicowaniem przebiegu pogody w poszczególnych latach.
Pszenżyto ozime jest uprawiane w Polsce od początku lat osiemdziesiątych jako zboże paszowe, a jego udział w strukturze zasiewów wynosi nieco ponad 10% (GUS2013). Na wzrost, rozwój i plonowanie roślin ozimych wpływają za-równo wskaźniki meteorologiczne okresu wegetacyjnego, jak również okresu jesienno-zimowego (Rymuza i in. 2012). W literaturze zdecydowanie więcej opracowań dotyczy wpływu warunków meteorologicznych okresu wegetacyjnego na plonowanie pszenżyta, niż zależności plonowania od warunków okresu jesien-no-zimowego (Jaczewska-Kalicka 2008, Pisulewska i in. 1998).
Warunki meteorologiczne okresu jesieni i zimy zasadniczo wpływają na prze-zimowanie roślin, odporność na choroby i szkodniki, a także na niesprzyjające ekstremalne warunki pogodowe w okresie wiosennej wegetacji – co istotnie wpływa na wielkość plonu.
Od warunków meteorologicznych w tym okresie zależy ponadto wartość bio-logiczna nasion oraz zawartość składników pokarmowych (Stankiewicz i in. 2003). W organizmach zwierzęcych składniki mineralne wpływają na prawidłowy przebieg różnych procesów metabolicznych. Zasadniczą drogą ich uzupełniania są odpowiedniej jakości pasze (Stankiewicz i in. 2003).
O zawartości składników organicznych i mineralnych w ziarnie decydują ta-kie czynniki jak: gatunek, odmiana, warunki glebowe, zabiegi agrotechniczne – szczególnie nawożenie i ochrona roślin – a także warunki pogodowe.
Celem niniejszej pracy było zbadanie wpływu wybranych wskaźników mete-orologicznych w okresie jesienno-zimowym na plonowanie pszenżyta ozimego w latach 1993-2013.
MATERIAŁY I METODY BADAŃ
Badania z pszenżytem ozimym prowadzono w Tomaszkowie, oddalonym o ok. 6 km na zachód od Olsztyna. Doświadczenia ścisłe realizowano na glebie brunatnej właściwej, zaliczanej do kompleksu pszennego dobrego, metodą blo-ków losowych w 4 powtórzeniach. Co kilka lat aktualizowano odmiany pszenżyta zalecane do uprawy w rejonie Warmińsko-Mazurskim (Malno, Bogo, Woltario, Grenado).
Plony ziarna pszenżyta (Y) pozyskiwano z obiektów standardowych o śred-nim poziomie nawożenia (N 120 kg·ha–1, P 31 kg·ha–1, K 83 kg·ha–1).
Okres zimowy ustalono, przyjmując za jego początek po wystąpieniu średniej dobowej wartości temperatury powietrza ≤0°C w ciągu pięciu kolejnych dni (Atlas... 2001). Jako odwilż atmosferyczną przyjęto co najmniej dwudniowe ocie-plenie w okresie jesienno-zimowym, ze średnią dobową temperaturą powietrza > 0°C (Atlas... 1990).
Warunki klimatyczne w okresie jesienno-zimowym scharakteryzowano na podstawie kształtowania się następujących wskaźników meteorologicznych:
x1 – liczba dni od siewu do początku okresu zimowego,
x2 – średnia dobowa temperatura powietrza od siewu do okresu zimowego, x3 – suma opadów od siewu do okresu zimowego,
x4 – długość okresu zimowego,
x5 – suma temperatur średnich dobowych zimy, x6 – suma opadów okresu zimowego,
x7 – długość zalegania pokrywy śnieżnej, x8 – liczba dni z odwilżami atmosferycznymi.
Do oceny zależności plonu pszenżyta ozimego od poszczególnych wskaźni-ków meteorologicznych wykorzystano analizę korelacji (Pearsona i Spearmana).
Łączny wpływ badanych wskaźników meteorologicznych na wielkość plonu pszenżyta ozimego oceniono, stosując analizę korelacji i regresji wielowymiaro-wej z wyborem optymalnego podzbioru zmiennych niezależnych.
Do obliczeń statystycznych i prezentacji wyników wykorzystano programy komputerowe Excel oraz Statistica ver. 12 (StatSoft… 2014).
WYNIKI I DYSKUSJA
Plony pszenżyta ozimego w 21-letnim okresie badań wahały się od 3,78 t·ha–1 w 1998 r. do 8,83 t·ha–1 w 2004 r. i wyniosły średnio 5,97 t·ha–1. Zmienność uzy-skanych wyników (SD = 1,44) pozwoliła na przyjęcie hipotezy o normalności rozkładu plonów. Ponadto, wielkość plonów w zależności od lat badań charakte-ryzowała się trendem istotnym statystycznie (p < 0,05) (rys. 1).
Warunki jesienno-zimowe dwudziestojednoletniego okresu badań charaktery-zowały się także dużą zmiennością, zbadano zatem normalność rozkładów po-szczególnych wskaźników w tym okresie (tab. 1). Z testu Shapiro-Wilka wynika, że spośród 8 wskaźników meteorologicznych uwzględnionych w badaniach tylko połowa (x2, x6, x7, x8) charakteryzuje się rozkładem normalnym.
Rys. 1. Plonowanie pszenżyta ozimego w Tomaszkowie w latach 1993-2013 Fig. 1. Yielding of winter triticale in Tomaszkowo in 1993-2013
Tabela 1. Warunki jesienno-zimowe wegetacji pszenżyta ozimego w Tomaszkowie w latach 1993-2013 Table 1. Autumn-winter conditions of winter triticale vegetation in Tomaszkowo in 1993-2013
Warunki wegetacji Conditions of vegetation
Średnia
Average Min Max SD
Normalność rozkładu; test S-W Normality of distribution; test S-W x1 – Liczba dni od siewu do okresu zimowego
Number of days from sowing to winter period 75,57 39 124 25,05 0,069
x2 – Średnia dobowa temperatura od siewu do okresu zimowego / Average daily temperature from sowing to winter period (°C)
6,54 4,5 10,1 1,35 0,328 (N)
x3 – Suma opadów od siewu do okresu zimowego
Total precipitation from sowing to winter period (mm) 101,95 37 246 46,81 0,034
x4 – Długość okresu zimowego (dni) Duration of
winter period (days) 86,10 14 128 32,46 0,037
x5 – Suma temp. średnich dobowych zimy
Total average daily temperature in winter period (°C) –228,10 –680 –10 169,06 0,010
x6 – Suma opadów okresu zimowego
Total precipitation in winter period (mm) 75,52 7 151 38,57 0,902 (N)
x7 – Długość zalegania pokrywy śnieżnej (dni)
Duration of snow cover retention (days) 66,51 13 116 29,68 0,398 (N)
x8 – Liczba dni z odwilżami atmosferycznymi
Number of days with atmospheric thaws 29,00 4 70 18,56 0,112 (N)
Na to, że pozostałe zmienne x1, x3, x4 i x5 odbiegały od rozkładu normalnego w znacznej mierze wpływały warunki atmosferyczne zaistniałe w 1993 i 2004 roku. Okres jesienno-zimowy w roku 1993 charakteryzował się specyficznym prze-biegiem warunków meteorologicznych. Dotyczy to długiego okresu od siewu do początku zimy (91 dni), gwałtownego spadku średnich dobowych temperatur po-wietrza poniżej 0°C, wskutek czego rośliny weszły w stan anabiozy nieodpowied-nio zahartowane. Ponadto, okres jesienny tego roku charakteryzował się niskimi temperaturami powietrza (4,5°C), dużą sumą opadów (140 mm) oraz częstymi odwilżami atmosferycznymi (35 dni). Natomiast warunki meteorologiczne oma-wianego okresu w 2004 roku charakteryzowały się stosunkowo niewielką liczbą dni od siewu do okresu zimowego – 52 dni, dużą sumą opadów – 128 mm, a w okresie zimowym krótko zalegającą pokrywą śnieżną – 65 dni oraz małą licz-bą dni z odwilżami atmosferycznymi – 12 dni. Warunki te skrajnie niekorzystnie wpływały na wegetację i zimowy spoczynek roślin, co z kolei negatywnie wpły-wało na wielkość plonu. Potwierdza to m.in. Czarnecka i Kalbarczyk (2004).
Ocenę zależności plonu pszenżyta od analizowanych wskaźników meteorolo-gicznych przedstawiono w tabeli 2.
Tabela 2. Zależność plonu pszenżyta ozimego od zmiennych określających wskaźniki jesienno-zimowe w latach 1993-2013. Współczynniki korelacji liniowej r lub korelacji rang R
Table 2. Winter triticale yield dependence on variables determining autumn-winter indicators in 1993-2013. Linear correlation coefficients r or rank correlation R
Zmienne – warunki wegetacji Variables – conditions of vegetation
Korelacja / Correlation r
(Pearsona)
R (Spearmana) x1 – Liczba dni od siewu do okresu zimowego
Number of days from sowing to winter period 0,033
x2 – Średnia dobowa temperatura powietrza od siewu do okresu zimowego
Average daily temperature from sowing to winter period (oC) 0,318
x3 – Suma opadów od siewu do okresu zimowego - Total
precipita-tion from sowing to winter period (mm) 0,085
x4 – Długość okresu zimowego (dni)
Duration of winter period (days) –0,270
x5 – Suma temp. średnich dobowych zimy
Total average daily temperature in winter period (oC) –0,121
x6 – Suma opadów okresu zimowego
Total precipitation in winter period (mm) –0,265
x7 – Długość zalegania pokrywy śnieżnej (dni)
Duration of snow cover retention (days) –0,271
x8 – Liczba dni z odwilżami atmosferycznymi
Number of days with atmospheric thaws –0,558*
Z danych zawartych w tabeli 2 wynika, iż plon pszenżyta zależy istotnie tylko od liczby dni z odwilżami atmosferycznymi i zależność ta jest ujemna. Tłumaczyć to można tym, że występowanie odwilży atmosferycznych, przechodzących często w odwilże glebowe powoduje między innymi demineralizację próchnicy, zuboże-nie kompleksu sorpcyjnego, wypłukiwazuboże-nie przyswajalnych składników pokarmo-wych przy często występujących opadach. Wskutek tego gleba traci właściwą strukturę i teksturę, co prowadzi do obniżenia jej żyzności i urodzajności (Gorlach i Mazur 2001). Zależność plonu od pozostałych zmiennych nie jest istotna.
Tabela 3. Łączny wpływ cech określających warunki jesienno-zimowe w Tomaszkowie w latach 1993-2013 na plon pszenżyta ozimego. Analiza wielozmiennej korelacji i regresji liniowej
Table 3. Total impact of characteristics defining autumn-winter conditions in 1993-2013 in Tomaszkowo on the yield of winter triticale. Multivariate analysis of correlation and linear regression
A n al iz o w an e z m ie n n e V ar ia b le s an al y se d L ic zb a d n i o d s ie w u d o o k re su z i-m o w eg o – N u m b e r o f d ay s fr o m so w in g t o w in te r p er io d Ś re d n ia d o b o w a te m p . p o w ie tr z a o d si ew u d o o k re su z im o w eg o – A v er a -g e d ai ly t em p er at u rę f ro m s o w in g t o w in te r p er io d ( °C ) S u m a o p ad ó w o d s ie w u d o o k re su zi m o w eg o – T o ta l p re ci p it at io n f ro m so w in g t o w in te r p er io d ( m m ) D łu g o ść o k re su z im o w eg o – D u ra ti o n o f w in te r p er io d ( d n i-d ay s) S u m a te m p . śr ed n ic h d o b o w y ch z im y – T o ta l av er ag e d ai ly t em p er at u re i n w in te r p er io d ( °C ) S u m a o p ad ó w o k re su z im o w eg o – T o ta l p re ci p it at io n i n w in te r p er io d ( m m ) D łu g o ść z al eg an ia p o k ry w y ś n ie żn ej – D u ra ti o n o f sn o w c o v er r et en ti o n ( d n i-d ay s) L ic zb a d n i z o d w il ża m i at m o sf e -ry cz n y m i – N u m b er o f d ay s w it h a tm o s-p h er ic t h aw s Y x1 x2 x3 x4 x5 x6 x7 x8
Plon / Yield (t·ha–1) Model pełny / Full model
Równanie regresji Regression equation: Y = f(x)
Y= 5,2480 – 0,0004·1 + 0,3863·2 + 0,0024·3 + 0,0196·4 + 0,0017·5 – 0,0152·6 – 0,0041·7 – 0,0669·8
Ocena regresji /
Regres-sion evaluation: R2; p R 2 = 0,624; p < 0,075* Znormalizowane wsp. regresji / Standarised regression coefficients: b’ –0,007 0,361 0,078 0,441 0,198 –0,407 –0,085 –0,859 p 0,982 0,113 0,786 0,319 0,378 0,073 0,731 0,010**
Plon / Yield (t·ha–1) Model optymalny – Optimal model
Równanie regresji / Regression equation: Y = f(x)
Y = 5,8977 + 0,3509·2 – 0,0112·6 – 0,0477·8 Ocena regresji /
Regres-sion evaluation: R2; p R 2 = 0,441; p < 0,005** Znormalizowane wsp. regresji / Standarised regression coefficients: b’ 0,328 –0,298 –0,913 p 0,068* 0,095* 0,002**
Łączny wpływ badanych wskaźników meteorologicznych okresu jesienno-zimowego na plonowanie pszenżyta oszacowano, wykorzystując analizę regresji wielozmiennej (tab. 3). Istotne równanie regresji uwzględnia wszystkie analizo-wane warunki wegetacji (tzw. model pełny) – przy czym najsilniej ujemnie na plon wpływają: suma opadów okresu zimowego i liczba dni z odwilżami atmosfe-rycznymi.
Zastosowanie procedur eliminacji nieistotnych zmiennych doprowadziło do oszacowania nowego istotnego równania regresji, uwzględniającego minimalny zestaw zmiennych (model optymalny). Na plon pszenżyta wpływają łącznie trzy zmienne: dodatnio średnia temperatura dobowa od siewu do okresu zimowego, a ujemnie suma opadów okresu zimowego oraz liczba dni z odwilżami atmosfe-rycznymi – przy czym najsilniej obniżająco wpływa liczba dni z odwilżami at-mosferycznymi.
Trudno jest zilustrować zależność plonu pszenżyta łącznie od trzech zmien-nych. Pokazano zatem, jak zmienia się plonowanie od każdej z tych trzech zmiennych na osobnych rysunkach.
Zależność plonu pszenżyta od średniej dobowej temperatury powietrza w okresie jesiennym (rys. 2) wykazuje trend dodatni, ale nieistotny statystycznie (r = 0,318), tzn. wraz ze wzrostem temperatury w tym okresie wzrasta wielkość plonu. Odpowiednie temperatury sprzyjają właściwemu wzrostowi i rozwojowi roślin, wskutek tego rośliny mają dostatecznie wykształconą część nadziemną i podziemną, co warunkuje właściwe przezimowanie roślin. Niskie temperatury w tym okresie obniżają wielkość plonu, potwierdzają to również badania Czar-neckiej i Kalbarczyka (2004).
Rys. 2. Zależność plonu pszenżyta ozimego (t·ha–1) od średniej temperatury dobowej od siewu do okresu zimowego (°C)
Fig. 2. Winter triticale yield dependence (t ha–1) on average daily temperature from sowing to winter period (°C)
Na rysunku 3 przedstawiono zależność wielkości plonu pszenżyta ozimego od sumy opadów okresu zimowego. Trend tej zależności ma wartość ujemną, tzn. wzrost opadów w tym okresie generuje niewielki spadek plonu.
Rys. 3. Zależność plonu pszenżyta ozimego (t·ha–1) od sumy opadów okresu zimowego (mm) Fig. 3. Winter triticale yield dependence (t ha–1) on total precipitation in winter period (mm)
Spośród badanych wskaźników meteorologicznych wpływających na plono-wanie pszenżyta najwyższy spadek wielkości plonu powodują odwilże atmosfe-ryczne w okresie zimowym (rys. 4). Trend spadku plonu w zależności od odwilży atmosferycznych – wysoki (r = –0,558).
Rys. 4. Zależność plonu pszenżyta ozimego (t·ha–1) od liczby dni z odwilżami atmosferycznymi Fig. 4. Winter triticale yield dependence (t ha–1) on number of days with atmospheric thaws
WNIOSKI
1. Plonowanie pszenżyta ozimego oraz przebieg warunków meteorologicznych w latach 1993-2013 były znacznie zróżnicowane w poszczególnych latach badań.
2. Spośród analizowanych łącznie wskaźników meteorologicznych istotnie dodatnio na plonowanie pszenżyta ozimego wpływały średnie dobowe temperatu-ry powietrza od siewu do początku okresu zimowego, a ujemnie sumy opadów i odwilże atmosferyczne w okresie zimy.
3. Liczba dni z odwilżami atmosferycznymi okazała się jedyną istotną cechą spośród analizowanych wskaźników meteorologicznych wpływających bezpo-średnio na plon pszenżyta ozimego.
PIŚMIENNICTWO
Atlas klimatyczny elementów i zjawisk szkodliwych dla rolnictwa w Polsce. 1990. Red. Koźmiński C., Górski T., Michalska B. Wyd. IUNG – Puławy, AR, Szczecin.
Atlas klimatycznego ryzyka uprawy roślin w Polsce. 2001. Red. Koźmiński C., Michalska B. Wyd. AR, Szczecin.
Czarnecka M., Kalbarczyk E., 2004. Wpływ warunków meteorologicznych okresu jesiennego na zimowanie pszenżyta w Polsce. Acta Agrophysica, 3(2), 221-234.
Gorlach E., Mazur T., 2001. Chemia rolna. Podstawy żywienia i zasady nawożenia roślin. Wyd. PWN Warszawa.
Jaczewska-Kalicka A., 2008. Wpływ zmian klimatycznych na plonowanie i ochronę zbóż w Polsce. Postępy w Ochronie Roślin, 48(2), 415-425.
Pisulewska E., Zając T., Oleksy A., 1998. Skład mineralny ziarna wybranych odmian pszenżyta ozimego w warunkach zróżnicowanego nawożenia azotem. Biul. IHAR, 205/206, 179-188. GUS 2013. Roczniki Statystyczne. Zakład Wydawnictw Statystycznych, Warszawa.
Rymuza K., Marciniuk-Kluska A., Bombik A., 2012. Plonowanie zbóż ozimych w zależności od warunków termiczno-opadowych na polach produkcyjnych rolniczej stacji doświadczalnej w Zawadach. Woda-Środowisko-Obszary Wiejskie, 12, 2(38), 207-220.
Stankiewicz C., Gruszecka D., Steć E., Mitrus J., Walo P., Kamecki M., 2003. Effect of interaction of parental components on the content of protein and amino acids in triticale grain. EJPAU, 6(2), 6. StatSoft, Inc. (2014). STATISTICA, version 12. www.statsoft.com.
YIELDING OF WINTER TRITICALE IN RELATION TO METEOROLOGICAL CONDITIONS IN AUTUMN-WINTER PERIOD
IN THE WARMIA AND MAZURY REGION
Jan Grabowski, Katarzyna Pożarska, Krystyna Grabowska
Department of Water, Climate and Environmental Management University of Warmia and Mazury in Olsztyn
Pl. Łódzki 1, 10-727 Olsztyn e-mail: jangrab@uwm.edu.pl
A b s t r a c t: The quantity and quality of yields of cereals are dependent to a large extent on the weather conditions during the vegetation season, as well as on the development of the meteorological indicators in the autumn-winter period. These conditions, especially in Poland, cause a large variability of yielding, on average at the level of 30-40%. The aim of this study was to evaluate the impact of meteorological indicators which occurred in the autumn-winter period on the yielding of winter tritica-le in 1993-2013 in the region of Warmia and Mazury. Triticatritica-le is one of the main cereal plants culti-vated in Poland. Its participation in the total crop structure is approx. 10%, and average yields are in
the range of 3-4 t ha–1. The following meteorological indicators were included in the study: number of
days from sowing to winter period, average daily temperature from sowing to winter period, total precipitation from sowing to winter period, duration of winter period, total average daily temperature in winter period, total precipitation in winter period, duration of snow cover retention, number of days with atmospheric thaws. The impact of the examined meteorological indicators in autumn-winter period on the yielding of winter triticale was analysed statistically using the analysis of correlation and regression for one or more independent variables. In the analysed 21-year period the yields of winter
triticale averaged 5.97 t ha–1. The lowest were in 1998 – 3.78 t ha–1, the highest – 8.83 t ha–1 – in 2004.
The conducted analyses indicate that the yield of winter triticale significantly depended only on the number of days with atmospheric thaws. The total impact on yield was mainly determined by the number of days with atmospheric thaws in winter period, average daily temperature from sowing to winter period and total precipitation in winter period.
