邊遠(yuǎn) 滕飛 平冰宇 王紅宇

摘要 為研究氣候因素對(duì)吉林市春玉米產(chǎn)量的影響，基于2000—2020年的農(nóng)業(yè)和產(chǎn)量數(shù)據(jù)，對(duì)吉林市的氣候因子變化和玉米單產(chǎn)變化進(jìn)行相關(guān)分析，并運(yùn)用逐步回歸方法建立氣象單產(chǎn)和氣候因素的數(shù)學(xué)模型。結(jié)果表明：近21年來吉林市春玉米生育期內(nèi)除日照時(shí)數(shù)外，年平均氣溫、總降水量、年平均風(fēng)速和濕度的線性均呈逐年上升趨勢(shì)；玉米實(shí)際單產(chǎn)和氣象單產(chǎn)線性在2010年以前呈下降趨勢(shì)，在2010年之后為上升趨勢(shì)，氣象單產(chǎn)數(shù)據(jù)波動(dòng)較大；日照時(shí)數(shù)與玉米單產(chǎn)存在顯著的正相關(guān)關(guān)系。研究旨在探尋對(duì)吉林市春玉米產(chǎn)量影響顯著的因子，從而合理規(guī)劃種植方案，以實(shí)現(xiàn)糧食增產(chǎn)。

關(guān)鍵詞 氣候因素；玉米產(chǎn)量；相關(guān)分析；逐步回歸

中圖分類號(hào)：S513 文獻(xiàn)標(biāo)識(shí)碼：B 文章編號(hào)：2095–3305（2023）09–0-04

糧食安全是“國之大者”。黨的二十大報(bào)告中提出，“全方位夯實(shí)糧食安全根基”“確保中國人的飯碗牢牢端在自己手中”?？梢?，糧食對(duì)人民和國家有著舉足輕重的影響[1]。玉米是重要的糧食產(chǎn)物，然而影響玉米產(chǎn)量高低的因素很多，如土壤狀況、生產(chǎn)力水平和氣候條件等，面對(duì)當(dāng)今嚴(yán)峻的全球變暖問題，研究玉米產(chǎn)量如何受氣候變化影響成為重要的課題。

近年來，對(duì)玉米產(chǎn)量與氣候因子之間的關(guān)系研究，許多學(xué)者進(jìn)行探討。羅瑞林[2]對(duì)內(nèi)蒙古玉米種植區(qū)的溫度和降水量變化趨勢(shì)進(jìn)行分析，并進(jìn)行肯德爾檢驗(yàn)，計(jì)算干燥度指數(shù)，分析溫度、降水量和玉米生長(zhǎng)發(fā)育的相關(guān)性及其對(duì)玉米產(chǎn)量的影響。侯英雨等[3]提出作物生長(zhǎng)日可利用降水量的推算方法和氣象指標(biāo)時(shí)間插值算法，構(gòu)建了東北地區(qū)日尺度氣候適宜度模型。趙錦等[4]學(xué)者利用APSIM-Maize模型模擬東三省春玉米可能種植區(qū)域內(nèi)氣象站點(diǎn)逐年的雨養(yǎng)產(chǎn)量，并結(jié)合統(tǒng)計(jì)學(xué)方法，分析春玉米雨養(yǎng)產(chǎn)量高產(chǎn)性和穩(wěn)產(chǎn)性的變化特征。王占彪等[5]通過分析華北平原的氣候數(shù)據(jù)，分析夏玉米營(yíng)養(yǎng)生長(zhǎng)期、并進(jìn)期、生殖生長(zhǎng)期和全生育期的生長(zhǎng)度日、高溫度日、降水量及氣候傾向率的時(shí)空分布。

綜上所述，學(xué)者可進(jìn)一步完善細(xì)化氣候產(chǎn)量的分離研究，為適應(yīng)氣候條件使農(nóng)業(yè)產(chǎn)量達(dá)到最大化，選取2000—2020年吉林市的氣象數(shù)據(jù)和農(nóng)業(yè)數(shù)據(jù)。首先，將統(tǒng)計(jì)年鑒中吉林市歷年玉米總產(chǎn)量和種植面積數(shù)據(jù)進(jìn)行預(yù)先處理計(jì)算玉米單產(chǎn)，為后續(xù)氣象產(chǎn)量求解奠定基礎(chǔ)；其次，將各氣候因子先單獨(dú)考慮，分析單一因素與氣象產(chǎn)量的相關(guān)性；最后，分析各氣候因子對(duì)產(chǎn)量的綜合作用并求解逐步回歸方程，為農(nóng)業(yè)生產(chǎn)部門和政府相關(guān)部門制定相關(guān)風(fēng)險(xiǎn)防范措施、提高農(nóng)作物產(chǎn)量、研究氣候變化特征提供有力的論證依據(jù)。

1 研究區(qū)域概況

吉林市地處43°31′～44°40′N，125°40′～127°56′E之間，位于吉林省中部偏東方向，是吉林省的第二大城市，位于長(zhǎng)白山區(qū)向松遼平原過渡地帶，地貌類型復(fù)雜，有“遠(yuǎn)迎長(zhǎng)白，近繞松花”之勢(shì)，土地肥沃，適宜水稻、旱糧農(nóng)事耕作。吉林市屬于溫帶大陸性季風(fēng)氣候，四季分明，春季少雨干燥，夏季溫暖多雨，秋季涼爽晴朗，冬季漫長(zhǎng)而寒冷，氣溫受地形影響，由西、西北向東、東南氣溫逐漸下降[6]。

2 數(shù)據(jù)來源與處理

2.1 數(shù)據(jù)來源

氣象數(shù)據(jù)來自中國氣象網(wǎng)的V3.0地面氣象觀測(cè)數(shù)據(jù)集，包含吉林市基本氣象站2000—2020年的氣象數(shù)據(jù)，而且對(duì)缺測(cè)數(shù)據(jù)進(jìn)行線性內(nèi)插，同時(shí)刪除異常值，最后選取生育期內(nèi)年平均氣溫T、總降水量P、日照時(shí)數(shù)SH、年平均風(fēng)速F、年平均濕度H作為氣候因子變量加以研究。玉米產(chǎn)量等農(nóng)業(yè)數(shù)據(jù)來自吉林省統(tǒng)計(jì)年鑒中2000—2020年各個(gè)市、縣的主要農(nóng)作物總播種面積和產(chǎn)量統(tǒng)計(jì)等。

2.2 數(shù)據(jù)處理

重點(diǎn)探討氣候因子對(duì)氣象產(chǎn)量的影響，從吉林省統(tǒng)計(jì)年鑒中整理2000—2020年吉林市的玉米總產(chǎn)量X和播種面積S，使用Excel計(jì)算每年的玉米實(shí)際單產(chǎn)Y，然后使用一階差分法計(jì)算相鄰兩年春玉米產(chǎn)量的差值，即可得到春玉米氣象單產(chǎn)Yw（t/ha）：

Y=X / S（1）

Yw=Yn+1-Yn（2）

式（1）、（2）中，n=1，2，…，20。

3 結(jié)果與分析

3.1 氣象數(shù)據(jù)分析

分析并繪制吉林市2000—2020逐年的氣象數(shù)據(jù)折線圖（圖1），以時(shí)間序列作為自變量，將各氣候因子的觀測(cè)值作為因變量，觀察氣象數(shù)據(jù)的線性變化趨勢(shì)[7]。

觀察圖1a可知，受全球氣候變暖的影響，吉林市在2000—2020年玉米生育期內(nèi)的年平均溫度的線性趨勢(shì)呈現(xiàn)逐年遞增狀態(tài)，氣候傾向率為0.127 ℃/10年。21年間，受高溫少雨旱情影響，平均氣溫最高出現(xiàn)在2000年，而年平均氣溫最低出現(xiàn)在2009年；觀察圖1b可知，總降水量變化幅度最大，線性趨勢(shì)呈現(xiàn)上升狀態(tài)，增長(zhǎng)速率分別為11.888 mm/年，春季降水偏少，夏季降水較多。2009年降水量因旱災(zāi)出現(xiàn)突破性的低值，而2019年由于強(qiáng)降水導(dǎo)致降水量達(dá)到最大值，降水量相對(duì)適中但總體趨勢(shì)變化不穩(wěn)定；觀察圖1c可知，總?cè)照諘r(shí)數(shù)呈遞減趨勢(shì)，遞減速率分別為3.100 5 h/年，但高日照時(shí)數(shù)的年份較多，有17年的總?cè)照諘r(shí)數(shù)超過1 000 h。2005年日照時(shí)數(shù)是歷史同期的最低值，為879.08 h，而2002年日照時(shí)數(shù)最高可達(dá)1 208.16 h；觀察圖1d可知，年平均風(fēng)速的線性趨勢(shì)變化不大，風(fēng)速基本為1.7～2.2 m/s，吉林市常年主導(dǎo)風(fēng)向?yàn)槲髂巷L(fēng)，平均風(fēng)速最大的年份是2019年，風(fēng)速值為2.20 m/s，平均風(fēng)速最小的年份為2012年，風(fēng)速值為1.80 m/s；觀察圖1e可知，年平均濕度的線性趨勢(shì)呈現(xiàn)逐年遞增的狀態(tài)，速率為0.163 1%/年。平均濕度最大的在2020年，濕度值為77.85%，平均濕度最小的為2007年，濕度值為68.32%。

3.2 玉米產(chǎn)量數(shù)據(jù)分析

觀察2000—2020年吉林市春玉米的實(shí)際單產(chǎn)和氣象單產(chǎn)分成2個(gè)階段（圖2），第一階段是2000—2009年，在此期間春玉米實(shí)際單產(chǎn)和氣象單產(chǎn)均呈現(xiàn)下降趨勢(shì)，實(shí)際單產(chǎn)以每年0.037 t/ha的速率下降，同時(shí)氣象產(chǎn)量以每年0.058 7 t/ha的速率下降；第二階段是2010—2020年，春玉米實(shí)際單產(chǎn)和氣象單產(chǎn)均呈現(xiàn)上升趨勢(shì)，實(shí)際單產(chǎn)以每年0.047 3 t/hm2的速率上升，而氣象產(chǎn)量以每年0.015 1 t/hm2的速率上升，兩者同步變化。綜上，氣象產(chǎn)量的波動(dòng)變化間接證明了氣候因子極大地影響玉米的產(chǎn)量[8]。

3.3 單一氣候因子對(duì)吉林市春玉米產(chǎn)量的相關(guān)性分析

首先，對(duì)產(chǎn)量和氣象數(shù)據(jù)進(jìn)行正態(tài)性檢驗(yàn)，經(jīng)驗(yàn)證，所有要素全部具備正態(tài)性特質(zhì)。其次，采用一階差分法進(jìn)一步處理氣象數(shù)據(jù)，一階差分既可以清除序列所包含的隨機(jī)性趨勢(shì)，又可以降低氣象數(shù)據(jù)對(duì)時(shí)間的依賴性，即減弱時(shí)間對(duì)氣象數(shù)據(jù)的影響[9]。最后，將玉米單產(chǎn)的變化值設(shè)為?yw，將氣候因子的變化量設(shè)為?xi（i=1，2，3，4，5，分別代表年平均氣溫、總降水量、日照時(shí)數(shù)、年平均風(fēng)速和年平均濕度），對(duì)單一氣候因子和玉米產(chǎn)量分別進(jìn)行Pearson相關(guān)分析，結(jié)果見表1，并對(duì)顯著性高的氣候因子求出如下線性回歸方程：

?yw=k?xi+b（3）

式（3）中，k代表一次回歸方程的系數(shù)，b代表一次回歸方程的常數(shù)項(xiàng)。

觀察表1可知，2000—2020年吉林市春玉米單產(chǎn)變化?yw只和日照時(shí)數(shù)變化?x3之間存在顯著的正相關(guān)關(guān)系，相關(guān)系數(shù)為0.555（P＜0.05），表示中等程度的相關(guān)，兩者線性回歸結(jié)果見表2。

結(jié)合表1中結(jié)論和對(duì)玉米的現(xiàn)有研究可知，玉米作為一種喜光的作物，在種子發(fā)芽期，因內(nèi)部此時(shí)沒有生成葉綠體，其生長(zhǎng)只需空氣、溫度和水分。從種子萌發(fā)后，幼苗細(xì)胞中出現(xiàn)葉綠體，開始進(jìn)行光合作用。玉米出苗后，若長(zhǎng)期處于短日照條件下，則發(fā)育加快、植株矮小、提前抽雄開花而降低產(chǎn)量；若長(zhǎng)期處于長(zhǎng)日照條件下，則植株增高、莖葉繁茂、抽雄開花期延遲，甚至不能開花結(jié)果[10]。綜上所述，從相關(guān)分析結(jié)果和作物生長(zhǎng)特性角度均驗(yàn)證日照時(shí)數(shù)變化對(duì)氣象單產(chǎn)變化有顯著的影響。

其他氣候因素對(duì)氣象單產(chǎn)的影響不顯著，不代表對(duì)產(chǎn)量沒有影響，只是在年平均尺度溫度、降水、風(fēng)速和濕度的變化不顯著，所以相關(guān)系數(shù)相對(duì)較低。實(shí)際上，從玉米播種到收獲，溫度影響力也不容忽視，春玉米的出苗時(shí)間和生長(zhǎng)發(fā)育的速度都離不開溫度的作用，溫度太高或者太低都不適合春玉米的種植。玉米葉片在進(jìn)行光合作用時(shí)，降水是合成生長(zhǎng)所需有機(jī)物質(zhì)和分配礦物元素的條件之一，包括玉米體溫的調(diào)節(jié)和在拔節(jié)、抽穗、開花階段都離不開水分。風(fēng)速對(duì)玉米的產(chǎn)量也有一定的影響，有助于玉米的光合作用，還可幫助其授粉，過大的風(fēng)速會(huì)使玉米葉片掉落，縮小光合作用的面積，正常生長(zhǎng)受到阻礙。玉米對(duì)環(huán)境中濕度的變化較為敏感，合適的濕度可以促進(jìn)其生長(zhǎng)發(fā)育。發(fā)芽期是玉米種植的重要階段，此時(shí)胚根開始生長(zhǎng)，合適的濕度值應(yīng)保持為70%～80%，有利于提供充足的水分；生長(zhǎng)期是玉米種植的主要階段，合適的濕度有利于玉米植株生長(zhǎng)；抽穗期是玉米開始開花的階段，合適的濕度對(duì)其傳粉和受精有益；成熟期是玉米種植的最后階段，合適的濕度可保證玉米的質(zhì)量和儲(chǔ)存，以防潮濕和發(fā)霉。

由表2可知，將日照變化作為自變量，而將氣象單產(chǎn)變化作為因變量進(jìn)行線性回歸分析，公式如下：

?yw=0.085+0.010*?x3（4）

式（4）中，模型R2=0.308代表日照變化可以解釋氣象單產(chǎn)變化30.8%的變化原因，模型通過F檢驗(yàn)（F=7.550，P=0.014＜0.05），即說明日照變化會(huì)對(duì)氣象單產(chǎn)變化產(chǎn)生影響關(guān)系，且日照變化的回歸系數(shù)值為0.010（t=2.748，P=0.014＜0.05），意味著日照變化會(huì)對(duì)氣象單產(chǎn)變化產(chǎn)生顯著的正向影響關(guān)系。

3.4 氣候因子對(duì)吉林市春玉米產(chǎn)量的回歸分析

研究單一氣候因子對(duì)春玉米氣象單產(chǎn)的影響，初步發(fā)現(xiàn)日照時(shí)數(shù)和氣象單產(chǎn)呈顯著正相關(guān)。然而，實(shí)際生活中各氣候因素間互相作用，表1中降水與溫度也呈顯著正相關(guān)，即存在多重共線性[11]。因此，建立氣象單產(chǎn)與氣候因子的逐步回歸方程，t檢驗(yàn)判斷顯著性時(shí)，須滿足P＜0.05，變量才能在方程中存在，方程如下：

yw=αx1+βx2+γx3+φx4+Φx5+c（5）

式（5）中，α、β、γ、φ、Φ代表方程中對(duì)應(yīng)變量的系數(shù)，c代表常數(shù)，結(jié)果見表3。

將全部氣象因子作為自變量，而將氣象產(chǎn)量作為因變量進(jìn)行逐步回歸分析，經(jīng)過模型自動(dòng)識(shí)別，最終余下日照一項(xiàng)在模型中，R2=0.028，意味著日照可以解釋氣象產(chǎn)量的22.8%變化原因。而且模型通過F檢驗(yàn)（F=5.306，P=

0.033＜0.05），說明模型有效[12]。公式如下：

yw= -7.378+0.007*x3（6）

式（6）中，日照的回歸系數(shù)值為0.007

（t=2.303，P=0.033＜0.05），意味著日照時(shí)數(shù)與吉林市春玉米氣象單產(chǎn)存在顯著的正向影響關(guān)系[13]。

4 結(jié)論

分析了吉林市2000—2020年春玉米生育期內(nèi)氣候因素和產(chǎn)量的線性趨勢(shì)，并對(duì)兩者的相關(guān)性和回歸方程進(jìn)行重點(diǎn)求解，研究結(jié)論和多位學(xué)者的成果也基本吻合。但是，玉米產(chǎn)量的影響因素不僅僅只有氣候，還包括品種等，雖然采取一階差分法得到氣象單產(chǎn)，仍有意外情況如臨時(shí)的農(nóng)業(yè)政策調(diào)整和玉米價(jià)格波動(dòng)等，需要進(jìn)一步提高方法準(zhǔn)確度；相關(guān)分析和逐步回歸方程結(jié)果顯示氣象單產(chǎn)只與日照時(shí)數(shù)呈顯著正相關(guān)關(guān)系，而與其他氣候因素?zé)o關(guān)，并不代表真實(shí)情況，可能其他氣象數(shù)據(jù)的年際差異較小影響不顯著。針對(duì)玉米產(chǎn)量的研究下一步可以從自然因素和非自然因素入手，應(yīng)用更全面精準(zhǔn)的數(shù)學(xué)模型控制變量助力農(nóng)業(yè)生產(chǎn)。本研究主要結(jié)論如下：

（1）從氣候因素變化趨勢(shì)角度分析可知，吉林市近21年春玉米生育期的氣候變化除日照時(shí)數(shù)外，年平均氣溫、總降水量、年平均風(fēng)速和濕度的線性變化趨勢(shì)均呈逐年上升，上升速率分別為0.012 7 ℃/年、11.888 mm/年、0.009 9 m/s/年和0.163 1%/年，而日照時(shí)數(shù)的下降速率為3.100 5 h/年。

（2）從春玉米產(chǎn)量角度分析可知，吉林市的玉米實(shí)際單產(chǎn)和氣象單產(chǎn)在2010前呈下降趨勢(shì)，從2010—2020年呈上升趨勢(shì)，氣象單產(chǎn)數(shù)據(jù)波動(dòng)較大，說明氣候因素確實(shí)對(duì)玉米產(chǎn)量產(chǎn)生不可忽視的影響。

（3）從全部氣候因素對(duì)春玉米影響角度分析可知，應(yīng)用Pearson相關(guān)分析和逐步回歸方程求解，發(fā)現(xiàn)日照時(shí)數(shù)是對(duì)吉林市春玉米單產(chǎn)有顯著正影響的氣候變量，為合理規(guī)劃農(nóng)業(yè)生產(chǎn)提供有力的科學(xué)依據(jù)。

參考文獻(xiàn)

[1] 盧景波.深入學(xué)習(xí)貫徹習(xí)近平總書記關(guān)于糧食安全的重要論述 不斷提升國家糧食安全保障能力和水平[J].中國糧食經(jīng)濟(jì)，2022，35（12）：4-8.

[2] 羅瑞林.氣候變化對(duì)內(nèi)蒙古春玉米產(chǎn)量影響的研究[D].呼和浩特：內(nèi)蒙古農(nóng)業(yè)大學(xué)，2013.

[3] 侯英雨，張艷紅，王良宇，等.東北地區(qū)春玉米氣候適宜度模型[J].應(yīng)用生態(tài)學(xué)報(bào)，2013，24（11）：3207-3212.

[4] 趙錦，楊曉光，劉志娟，等.全球氣候變暖對(duì)中國種植制度的可能影響Ⅹ：氣候變化對(duì)東北三省春玉米氣候適宜性的影響[J].中國農(nóng)業(yè)科學(xué)，2014，47（16）： 3143-3156.

[5] 王占彪，王猛，尹小剛，等.氣候變化背景下華北平原夏玉米各生育期水熱時(shí)空變化特征[J].中國生態(tài)農(nóng)業(yè)學(xué)報(bào)， 2015，23（4）：473-481.

[6] 鮑新華.吉林—長(zhǎng)春—四平城市經(jīng)濟(jì)帶表土環(huán)境質(zhì)量評(píng)價(jià)與生態(tài)地球化學(xué)分區(qū)[D].長(zhǎng)春：吉林大學(xué)，2011.

[7] 王柳，熊偉，溫小樂，等.溫度降水等氣候因子變化對(duì)中國玉米產(chǎn)量的影響[J].農(nóng)業(yè)工程學(xué)報(bào)，2014，30（21）：138-146.

[8] 車曉翠，李洪麗，張春燕，郭聃，沈海鷗.1980年代以來氣候變化對(duì)吉林省玉米產(chǎn)量的影響[J].水土保持研究，2021， 28（2）：230-234，241.

[9] 李輝，姚鳳梅，張佳華，等.東北地區(qū)玉米氣候產(chǎn)量變化及其對(duì)氣候變化的敏感性分析[J].中國農(nóng)業(yè)氣象，2014，35（4）： 423-428.

[10] 李秀穎.探析農(nóng)業(yè)氣象對(duì)玉米種植的影響及作用[J].中國農(nóng)業(yè)信息，2015，27 （22）：113-114.

[11] Liu Y， Wang E， Yang X G， et al. Contributions of climatic and crop varietal changes to crop production in the North China Plain， since 1980s[J]. Global Change Biology， 2010， 16（1）： 2289—2299.

[12] 曹永強(qiáng)，李玲慧，路潔，等.氣候變化對(duì)遼寧省玉米產(chǎn)量的影響[J].中國農(nóng)村水利水電，2020，62（11）：132-137.

[13] 梁烜赫，徐晨，王冰，等.吉林省不同生態(tài)區(qū)氣象因子對(duì)玉米產(chǎn)量影響的評(píng)價(jià)[J].吉林農(nóng)業(yè)科學(xué)，2015，40（4）：17-20.

Statistical Analysis of Climatic Factors on Spring Corn Yield in Jilin City

Bian Yuan et al（Jilin Normal University， Siping， Jilin 136000）

Abstract In order to study the impact of climate factors on spring corn yield in Jilin City， based on agricultural and yield data from 2000 to 2020， a correlation analysis was conducted between climate factor changes and corn yield changes in Jilin City， and a mathematical model for meteorological yield and climate factors was established using stepwise regression method. The results Showed that in the past 21 years， except for sunshine hours， the linear changes in annual average temperature， total precipitation， annual average wind speed， and humidity have all increased year by year during the growth period of spring corn in Jilin City; The actual yield of corn and the linear trend of meteorological yield showed a downward trend before 2010 and an upward trend after 2010， with significant fluctuations in meteorological yield data; There was a significant positive correlation between sunshine hours and maize yield per unit area. The research aimed to explore the factors that have a significant impact on the yield of spring corn in Jilin City， in order to plan planting plans reasonably to achieve grain yield increase.

Key words Climate factors; Corn yield; Correlation analysis; Stepwise regression