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

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

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

中圖分類號(hào) S513文獻(xiàn)標(biāo)識(shí)碼 A文章編號(hào) 1007-7731（2021）09-0147-03

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

1 資料來源

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

2 分析方法

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

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

2.2 EOF分析方法的計(jì)算步驟 計(jì)算每個(gè)特征向量的方差貢獻(xiàn)（貢獻(xiàn)率）：

及前m個(gè)特征向量的累積方差貢獻(xiàn)：

2 結(jié)果與分析

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

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

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

從時(shí)間系數(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個(gè)階段發(fā)生冷害的年份并去掉重復(fù)的年份，共22年發(fā)生低溫冷害，本文的研究范圍為48年，因此，低溫冷害基本上為2年一遇。由此可知：東北地區(qū)玉米出苗到七葉階段發(fā)生冷害的次數(shù)較少，且大多發(fā)生在60年代和70年代;出苗到抽雄、出苗到乳熟2個(gè)階段，70年代和80年代出現(xiàn)冷害的頻率較大;時(shí)間系數(shù)的最小值出現(xiàn)在出苗到成熟階段的1972年和1976年，而這2年前3個(gè)階段的時(shí)間系數(shù)均為正值，表明這2年在玉米即將成熟階段出現(xiàn)了較為嚴(yán)重的低溫冷害。

3 結(jié)論

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

參考文獻(xiàn)

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

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

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

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

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