東北地區(qū)玉米冷害指數(shù)時空分析

崔一鳴 王祿祿 李曉婷 姜新 許馨月 李銘鑒 辛得

摘 要：選取東北地區(qū)氣象觀測站1966—2013年的氣象資料以及玉米生育期資料，計算玉米4個生育階段的冷害指數(shù)，并對冷害指數(shù)進行EOF時空分解，分析了東北地區(qū)玉米低溫冷害的時空分布規(guī)律。

關(guān)鍵詞：低溫冷害;冷害指數(shù);EOF分解

中圖分類號 S513文獻標識碼 A文章編號 1007-7731（2021）09-0147-03

EOF分析方法最早是在1902年由英國統(tǒng)計學(xué)家Karl Pearson提出，在20世紀中期才由Lorenz初次引用至大氣學(xué)科當中。20世紀70年代中后期以來，EOF分析方法隨著計算機行業(yè)的迅猛發(fā)展，在大氣科學(xué)學(xué)科中的應(yīng)用變得更為簡單順暢。本文利用EOF分析方法，對東北地區(qū)玉米低溫冷害進行研究，以期為當?shù)赜衩變龊Ψ乐翁峁﹨⒖肌?/p>

1 資料來源

選取東北地區(qū)1966—2013年間36個氣象站逐日平均溫度及月平均溫度以及東北地區(qū)36個氣象臺站1966—2013年玉米播種、出苗、三葉、七葉、拔節(jié)、抽雄、開花、吐絲、乳熟及成熟日期。

2 分析方法

2.1 冷害指數(shù)的計算[1] 計算公式如下：

式中：T為某日的日平均溫度，T0為該日所在的生育期的適宜溫度，T1為該日所在的生育期的下限溫度，T2為該日所在的生育期的上限溫度。

2.2 EOF分析方法的計算步驟 計算每個特征向量的方差貢獻（貢獻率）：

及前m個特征向量的累積方差貢獻：

2 結(jié)果與分析

EOF分析特征向量的方差貢獻率越大，說明該特征向量的空間分布越典型[2]。對東北地區(qū)所選36個站點1966—2013年4個發(fā)育階段冷害指數(shù)矩陣分別進行EOF分析，得出特征向量及相應(yīng)的特征值、時間系數(shù)。4個階段的前3個模態(tài)的累積方差貢獻率分別為81.51%、80.55%、80.24%、87.66%，都超過了80%，說明前3個模態(tài)足可以代表東北地區(qū)低溫冷害的時空變化特征。東北玉米4個生育階段第1模態(tài)的方差貢獻率分別達到總方差的56.98%、61.09%、62.67%、77.01%，是方差貢獻最大的1個模態(tài)，說明第1模態(tài)可以代表東北玉米發(fā)育階段冷害最主要空間變化的模態(tài)，所以本文主要分析第1模態(tài)。

2.1 東北玉米4個階段冷害指數(shù)第1模態(tài)分析 第1模態(tài)中，玉米4個發(fā)育階段的特征向量符號一致，均為正值，說明東北玉米4個發(fā)育階段冷害的空間變化是基本一致[3]，且4個生育階段冷害的一致特征均占總方差的50%以上，表明東北地區(qū)冷害受大尺度天氣系統(tǒng)的影響較為明顯[4]。圖1為1966—2013年東北玉米4個階段冷害指數(shù)第1模態(tài)分布圖，從圖中可以看出：低于0.1的低值區(qū)從出苗七葉到出苗成熟，范圍不斷增加，出苗七葉階段主要分布在遼寧南部、中西部地區(qū)，到出苗成熟階段幾乎整個遼寧都位于低值區(qū)，說明遼寧地區(qū)冷害程度發(fā)生較輕，且冷害主要發(fā)生在生育階段的前期。0.1到0.2的中值區(qū)在出苗到七葉階段幾乎占據(jù)整個研究范圍的2/3，然后開始逐漸減弱，出苗到成熟階段只剩吉林西部、黑龍江小部分地區(qū)位于中值區(qū)。高于0.2的高值區(qū)隨時間不斷增加，出苗到七葉階段僅在吉林地區(qū)有小范圍的高值區(qū)，出苗到抽雄階段擴展到整個吉林省東部地區(qū)，出苗到乳熟階段開始明顯向北延伸，從出苗到成熟階段，幾乎整個黑龍江省以及吉林的中部、東部地區(qū)都位于高值區(qū)。吉林東北部地區(qū)在玉米生長發(fā)育的4個階段始終位于中高值區(qū)，說明吉林東北部對冷害最為敏感。

2.2 東北玉米4個階段冷害指數(shù)時間系數(shù)分析 圖2為東北玉米4個發(fā)育階段第1特征向量所對應(yīng)的時間系數(shù)變化曲線，代表東北玉米4個發(fā)育階段冷害的年際變化。出苗到七葉、出苗到抽雄和出苗到乳熟3個階段時間系數(shù)波動較小，表明這3個階段東北全區(qū)玉米低溫冷害時間分布較為一致;而出苗到成熟階段時間系數(shù)的波動較大，波動較大的年份為1971年和1975年。由于4個生育階段第1特征向量均為正值，根據(jù)EOF原理，若某年的時間系數(shù)為負，則該年的時間函數(shù)和空間函數(shù)的乘積也為負值，即冷害指數(shù)也為負值，則表明該年發(fā)生冷害，且值越小說明冷害越嚴重。從圖中可以看出，第1模態(tài)4個階段的時間系數(shù)均呈上升趨勢，這表明了東北地區(qū)玉米各發(fā)育階段冷害的強度呈減弱的趨勢。

從時間系數(shù)圖可以看出，出苗到七葉階段，發(fā)生冷害的年份有1969年、1973年、1974年和1997年;出苗到抽雄階段，發(fā)生冷害的年份有1969年、1972年、1973年、1974年、1983年、1987年、1992年和1995年;出苗到乳熟階段，發(fā)生冷害的年份有1969年、1971年、1972年、1974年、1976年、1977年、1981、1983、1986、1987、1989、1992、1993、1995年;出苗到成熟階段，發(fā)生冷害的年份有1966、1967、1968、1969、1971、1972、1974、1975、1976、1977、1981、1983、1987、1989、1990、1992、1995、1997、2009年。綜合4個階段發(fā)生冷害的年份并去掉重復(fù)的年份，共22年發(fā)生低溫冷害，本文的研究范圍為48年，因此，低溫冷害基本上為2年一遇。由此可知：東北地區(qū)玉米出苗到七葉階段發(fā)生冷害的次數(shù)較少，且大多發(fā)生在60年代和70年代;出苗到抽雄、出苗到乳熟2個階段，70年代和80年代出現(xiàn)冷害的頻率較大;時間系數(shù)的最小值出現(xiàn)在出苗到成熟階段的1972年和1976年，而這2年前3個階段的時間系數(shù)均為正值，表明這2年在玉米即將成熟階段出現(xiàn)了較為嚴重的低溫冷害。

3 結(jié)論

EOF分析結(jié)果表明：東北玉米4個發(fā)育階段冷害的空間變化是基本一致。遼寧地區(qū)冷害程度發(fā)生較輕，且冷害主要發(fā)生在生育階段的前期;吉林東北部地區(qū)對冷害最為敏感。當遼寧地區(qū)發(fā)生冷害時，黑龍江地區(qū)普遍偏暖。東北玉米第1模態(tài)時間系數(shù)表明：出苗到七葉、出苗到抽雄和出苗到乳熟3個階段東北全區(qū)玉米低溫冷害時間分布較為一致;而出苗到成熟階段時間系數(shù)的波動較大，波動較大的年份為1972年和1976年，說明這2年在乳熟到成熟階段發(fā)生了較為嚴重的冷害。低溫冷害基本上為2～3年一遇。

參考文獻

[1]崔一鳴，張明哲，王祿祿.1966～2013年東北地區(qū)玉米低溫冷害指數(shù)分析[J].現(xiàn)代農(nóng)業(yè)科技，2018，15：65～70

[2]高曉容，王春乙，張繼權(quán).東北玉米低溫冷害時空分布與多時間尺度變化規(guī)律分析[J].災(zāi)害學(xué)，2012，27（4）：65～70

[3]賈燕.基于EOF分析遼寧省極端氣溫時空分布的研究.安徽農(nóng)業(yè)科學(xué)[J]，2008，36（11）：4589-4590，4620

[4]張劍明，黃晴.近46年來湖南省氣溫變化的若干特點.干旱氣象[J]，2012，30（3）：387-403

（責(zé)編：王慧晴）