毛烏素沙漠與黃土高原過渡帶氣候變化對玉米產(chǎn)量的影響*

高 娟，馮靈芝，張建康，王 云，曹 梅，常善剛

(1.陜西省榆林市氣象局，榆林 719000； 2.陜西省西安市氣象局,西安 710016)

·糧食安全·

毛烏素沙漠與黃土高原過渡帶氣候變化對玉米產(chǎn)量的影響*

高 娟1，馮靈芝1，張建康1，王 云1，曹 梅2，常善剛1

(1.陜西省榆林市氣象局，榆林 719000； 2.陜西省西安市氣象局,西安 710016)

基于1971～2012年榆林市氣象數(shù)據(jù)和玉米產(chǎn)量數(shù)據(jù)，分析過去42年毛烏素沙漠與黃土高原過渡帶玉米生長期內(nèi)氣候因子的變化趨勢，并利用一階差分法、相關(guān)分析、多元逐步回歸法分析各生育階段氣候因子對玉米產(chǎn)量的影響。研究結(jié)果表明：過去42年，玉米生長期內(nèi)熱量資源呈顯著增加的趨勢，生長期平均溫度、平均最高溫度和平均最低溫度分別增加0.24℃/10 a、0.26℃/10 a 和0.32℃/10 a，降水量、光照時數(shù)無顯著變化。一階差分、回歸分析結(jié)果表明，玉米產(chǎn)量與全生長期平均溫度、平均最高溫度、積溫呈極顯著的負(fù)相關(guān)，其中與平均溫度相關(guān)系數(shù)最大，全生長期溫度每升高1℃，玉米減產(chǎn)835.5kg/hm2； 與降水量變化呈顯著正相關(guān)，降水量每增加100mm，玉米增產(chǎn)364.7kg/hm2； 與日照時數(shù)變化呈負(fù)相關(guān)，日照時數(shù)每增加100h，玉米減產(chǎn)325.8kg/hm2。各生育階段中，出苗-三葉期至灌漿乳熟期的平均最高溫度、平均溫度及苗期日照對玉米產(chǎn)量產(chǎn)生負(fù)影響，七葉-拔節(jié)期的降水和成熟期的平均最高溫度對玉米產(chǎn)量產(chǎn)生正影響； 多元逐步回歸模型分析結(jié)果表明，顯著影響玉米產(chǎn)量的是抽雄-開花期和灌漿-乳熟期平均最高溫度、拔節(jié)期降水量、出苗-三葉期日照，模型能夠解釋毛烏素沙漠與黃土高原過渡帶玉米56%的產(chǎn)量變化。

毛烏素沙漠 黃土高原 一階差分 氣候變化 玉米產(chǎn)量

0 引言

氣候變暖已成為國際公認(rèn)事實，政府間氣候變化專門委員會第5次評估報告再次指出氣候變暖，特別是從20世紀(jì)50年代以來，變暖趨勢尤為明顯，并且較IPCC第4次評估報告結(jié)果有進(jìn)一步升高之勢[1]。氣候變暖會影響到農(nóng)業(yè)、生態(tài)、經(jīng)濟(jì)等一系列問題[2-3]。國內(nèi)外大量研究表明，氣候變化給農(nóng)業(yè)帶來巨大挑戰(zhàn)，對糧食產(chǎn)量、種植布局、生產(chǎn)管理等方面會造成不利影響，增大未來農(nóng)業(yè)生產(chǎn)的不穩(wěn)定性[4-5]。各氣候因子中，溫度、降雨、光照等因子變化對玉米種植的影響較大，玉米更易受到氣候變化的負(fù)影響，影響幅度具有區(qū)域差異性。王柳等的研究指出，黃土高原及周邊地區(qū)平均溫度上升1℃，玉米平均減產(chǎn)21.6%[6]，呂碩等的研究表明，對于同一品種，溫度升高可使玉米產(chǎn)量潛力下降22%～26%[7]，侯艷林等通過試驗研究表明，氣候變化對玉米單產(chǎn)呈負(fù)面影響[8]，李言照的研究表明光、溫因子對玉米產(chǎn)量影響極大[9]，吉相臣等人研究中提到，在不澇不旱范圍內(nèi)，玉米產(chǎn)量與降水量呈正相關(guān)[10]，杜軍等研究表明，氣候變暖加劇了氣候暖干化發(fā)展，使玉米播種期提前[11]，屈洋等人研究指出，平均最高溫度變化對關(guān)中西部玉米產(chǎn)量影響最大[12]，扈艷萍等的研究指出，在發(fā)育期內(nèi)的降水對遼寧地區(qū)玉米產(chǎn)量影響最大[13]，馬雅麗等研究則指出，生長前期降水對山西地區(qū)玉米的影響最顯著[14]。

毛烏素沙漠和黃土高原過渡帶，生態(tài)環(huán)境脆弱，榆林境內(nèi)西北部為毛烏素風(fēng)沙草灘區(qū)，東南為黃土丘陵溝壑區(qū)。玉米是榆林現(xiàn)代特色農(nóng)業(yè)主導(dǎo)作物，常年播種面積10.666 7萬hm2，玉米年產(chǎn)量占全市糧食總產(chǎn)量的近50%，成為全市總產(chǎn)量第一的大宗農(nóng)作物，玉米種植收入是當(dāng)?shù)胤N糧農(nóng)民的重要經(jīng)濟(jì)來源。氣候變暖對生態(tài)環(huán)境脆弱區(qū)極為不利[15]。大量研究表明，氣候變暖，積溫增多，有使榆林農(nóng)耕期、作物生長期延長的趨勢[16]，而關(guān)于榆林氣候變化對玉米生長及產(chǎn)量的影響研究仍較少見，因此，文章基于1971～2012年氣候因子和玉米產(chǎn)量數(shù)據(jù)，研究過去42年玉米生長期內(nèi)溫度、降水、光照的變化趨勢，并分析氣候因子對玉米產(chǎn)量的影響，為合理利用氣候資源、保障玉米生產(chǎn)安全提供科學(xué)依據(jù)。

1 資料與方法

1.1 研究區(qū)域與數(shù)據(jù)來源

圖1 榆林市12縣區(qū)氣象觀測站點和農(nóng)氣觀測站點分布

毛烏素沙漠與黃土高原過渡帶主要分布于陜西省北部榆林市境內(nèi)，因此，該研究將榆林地區(qū)作為研究區(qū)域，圖1為榆林市12縣區(qū)氣象觀測站點和農(nóng)氣觀測站點分布圖。該區(qū)域面積4.357 8萬km2，屬于暖溫帶和溫帶半干旱大陸性季風(fēng)氣候，年平均溫度9.6℃，年降水量404.5mm。

氣象數(shù)據(jù)來源于榆林12個縣級地面氣象觀測站1971～2012年逐日觀測資料，其中平均溫度采用逐日4次定時觀測氣溫的平均值。對逐日平均溫度、日最高溫度、日最低溫度，日降水量、日照時數(shù)統(tǒng)計得到各要素逐月資料，即月平均溫度、月平均最高溫度、月平均最低溫度、月總降水量、月總?cè)照諗?shù)據(jù)等。

玉米生育期數(shù)據(jù)來源于榆陽區(qū)和綏德縣2個農(nóng)氣觀測站1999～2012年的農(nóng)氣觀測記錄，包括播種，出苗，三葉，七葉，拔節(jié)，抽雄，開花，吐絲，灌漿乳熟，成熟10個生育期。由于區(qū)域和年際差異性，該研究以5～9月為玉米生長期，榆陽區(qū)和綏德縣的生育期數(shù)據(jù)作為玉米生育期數(shù)據(jù)，并根據(jù)玉米生長特性，將全生育期劃分為5個生育階段，分別為出苗-三葉期，七葉-拔節(jié)期，抽雄-開花期，灌漿-乳熟期，成熟期。

玉米產(chǎn)量數(shù)據(jù)來源于榆林市統(tǒng)計局統(tǒng)計年鑒1971～2012年玉米產(chǎn)量數(shù)據(jù)。

1.2 研究方法與指標(biāo)

1.2.1 氣候傾向率

采用氣候傾向率作為氣候因子年際變化趨勢指標(biāo)，采用一元線性回歸方程描述：

Xi=a+bti(i=1,2,3,4,5，…，n)

(1)

其中，Xi為因變量，ti為自變量，i為年序號，ɑ為回歸常數(shù)，b代表氣候傾向率，ɑ、b用最小二乘法估計。

1.2.2 0-1標(biāo)準(zhǔn)化方法

(2)

其中，Xi表示因子數(shù)值，i為年序號，max(X)為樣本數(shù)據(jù)的最大值，min(X)為樣本數(shù)據(jù)的最小值。min-max標(biāo)準(zhǔn)化對基礎(chǔ)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，對原始數(shù)據(jù)作線性變換，使結(jié)果落到[0,1]區(qū)間。

1.2.3 一階差分法

目前相關(guān)文獻(xiàn)中，消除玉米產(chǎn)量趨勢效應(yīng)常采取的方法主要有直線滑動平均法分離趨勢產(chǎn)量、試驗田定位觀測、一階差分法等[17]。Liu Y等人的研究表明，一階差分法與直線滑動平均法在消除技術(shù)改進(jìn)、品種更替等緩慢變化對作物產(chǎn)量的影響時，其最終結(jié)果大致相同[18]。由于一階差分法較直線滑動平均法計算方法簡便，因此該文對玉米產(chǎn)量采用一階差分法進(jìn)行處理，其計算公式為：

ΔYi+1=Yi+1-Yi(i=1,2,3…n)

(3)

其中，Yi+1為第二年產(chǎn)量值，Yi為當(dāng)年產(chǎn)量值，i為年序號。采用第二年產(chǎn)量值減去當(dāng)年產(chǎn)量值，趨勢產(chǎn)量在求差值計算中已包含在實際產(chǎn)量中被消除，得到的差值即為氣象因子變化引起的產(chǎn)量變化。

1.2.4 多元逐步回歸模型

Y=a1X1+a2X2+a3X3+a4X4+a5X5+…anXn+b0

(4)

式中，Y代表玉米產(chǎn)量，X1，X2，X3…Xn分別代表各氣候因子，n表示氣候因子個數(shù)。a代表回歸系數(shù)，b0為常數(shù)項。采用SPSS19.0軟件進(jìn)行統(tǒng)計分析，選用多元逐步回歸法，設(shè)置F<0.05 進(jìn)入，F(xiàn)>0.1剔除。

研究首先基于1971～2012年榆林市氣象數(shù)據(jù)和玉米產(chǎn)量數(shù)據(jù)，分析過去42年毛烏素沙漠與黃土高原過渡帶玉米發(fā)育期內(nèi)氣候因子的變化趨勢； 其次，利用一階差分法消除趨勢產(chǎn)量，同時對玉米全生長期、各關(guān)鍵發(fā)育階段平均溫度、平均最高溫度、平均最低溫度、降水量、日照時數(shù)進(jìn)行一階差分處理，用相關(guān)分析法研究玉米產(chǎn)量變化與生長期各氣候因子變化的相關(guān)性，定量分析各氣候因子變化對玉米產(chǎn)量的影響； 研究各關(guān)鍵生育階段對應(yīng)月的氣候因子變化與產(chǎn)量變化的相關(guān)關(guān)系； 再次，對玉米產(chǎn)量數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，構(gòu)建各發(fā)育期的逐步回歸模型，進(jìn)一步分析各發(fā)育期氣候因子對玉米產(chǎn)量的綜合影響。

2 結(jié)果與討論

2.1 玉米生長期內(nèi)氣候變化

圖2為1971～2012年榆林玉米生長期氣候因子年際變化，由圖2a可見,1971～2012年榆林玉米生長期平均溫度、平均最高溫度和平均最低溫度均呈現(xiàn)明顯上升趨勢，上升幅度分別為0.24℃/10a、0.26℃/10a和0.32℃/10a，均達(dá)到了極顯著水平(P<0.01)，最低溫度上升速度高于最高溫度和平均溫度，生長期(≥10℃)積溫增多趨勢達(dá)到極顯著水平，上升幅度達(dá)68.63 ℃·d/10a[16]。

圖2 1971～2012年榆林玉米生長期氣候因子年際變化

由圖2b可見，過去42年，玉米生長期內(nèi)降水量有緩慢上升趨勢(5.3mm/10a)，總?cè)照諘r數(shù)有下降趨勢(17h/10a)，但變化趨勢均未達(dá)到顯著。表1為1971～2012年各縣區(qū)玉米生長期平均溫度和平均最高溫度變化趨勢，由表1可見，玉米生長期內(nèi)的平均溫度和平均最高溫度一致呈上升趨勢，各縣區(qū)變化趨勢差異較大，空間上表現(xiàn)為溫度增加速率由東南向西北遞增，其中東南部的佳縣、綏德、吳堡、清澗平均溫度上升趨勢未通過顯著性檢驗，西北部的神木、榆陽、府谷、橫山、靖邊、定邊增加0.19～0.42℃/10a； 榆陽、佳縣平均最高溫度未通過顯著性檢驗，其余10個縣增加0.21～0.35℃/10a，說明平均最高溫度升高更為普遍，西北部的神木、府谷、橫山、靖邊、定邊增加0.24～0.35℃/10a。

神木縣的平均溫度和平均最高溫度上升趨勢最大，平均溫度、平均最高溫度分別增加0.42℃/10a,0.35℃/10a，均達(dá)到極顯著水平(P<0.01)； 佳縣的平均溫度和榆陽區(qū)平均最高溫度升高速率最小，未通過顯著性檢驗。

表1 1971～2012年各縣區(qū)玉米生長期平均溫度和平均最高溫度變化趨勢 ℃/10 a

2.2 氣候因子變化對玉米產(chǎn)量的影響

圖3 1971～2012年榆林市玉米單產(chǎn)變化

圖3為1971～2012年榆林玉米單產(chǎn)變化圖，由圖3可見，過去42年榆林市玉米單產(chǎn)在波動中有急劇上升的趨勢，平均逐年增加92.41kg/hm2。研究采用一階差分法消除技術(shù)進(jìn)步、種植制度、品種更替等因素的影響，分析玉米單產(chǎn)變化量與生長期溫度、降水、日照變化量的關(guān)系。

2.2.1 生長期各氣候因子變化對玉米產(chǎn)量的影響

基于1971～2012年的玉米單產(chǎn)數(shù)據(jù)與生長期各氣候因子，通過一階差分、相關(guān)分析法計算逐年玉米單產(chǎn)變化和各氣候因子變化的相關(guān)關(guān)系，采用回歸分析法研究各氣候因子變化對產(chǎn)量的影響，結(jié)果如圖4所示。

由圖4可見，榆林地區(qū)玉米產(chǎn)量對生長期內(nèi)的平均溫度、平均最高溫度、積溫等氣候因子變化較敏感，總體來說表現(xiàn)為氣候變暖對玉米產(chǎn)量有負(fù)影響，其次降水量、日照時數(shù)的變化對玉米也有顯著影響，日照時數(shù)增多對玉米產(chǎn)量不利，降水量增多則利于玉米增產(chǎn)。各因子對玉米產(chǎn)量的影響具有差異性。其中，對玉米產(chǎn)量變化產(chǎn)生極顯著負(fù)影響的是平均溫度變化(P<0.01)，平均溫度每升高1℃，產(chǎn)量減少835.5kg/hm2(圖4a)； 其次為平均最高溫度(P<0.01)，生長期內(nèi)最高溫度每升高1℃，產(chǎn)量減少617.1kg/hm2(圖4c)； 再次為積溫變化(P<0.01)，生長期積溫每增加100℃·d，玉米減產(chǎn)255.3kg/hm2(圖4b)； 最低溫度變化對玉米產(chǎn)量變化影響未達(dá)到顯著水平(圖4d)。

其次為日照時數(shù)變化對玉米產(chǎn)量變化有顯著負(fù)影響(P<0.05)，日照時數(shù)每增加100h，玉米減產(chǎn)325.8kg/hm2(圖4f)。

對玉米產(chǎn)量變化產(chǎn)生顯著正影響的是生長期平均降水量(P<0.05)，降水量每增加100mm，使玉米產(chǎn)量增加364.7kg/hm2(圖4e)。

此結(jié)果與王柳、呂碩、吉相臣等人的研究結(jié)果趨向一致，KucharikCj等人對于美國威斯康星州30年氣候變化與玉米產(chǎn)量的分析也表明，夏季溫度每升1℃，會使得產(chǎn)量下降近13%，而氣候變涼、降水增多可提高玉米生產(chǎn)[19]，王春春等人研究也指出，溫度升高導(dǎo)致玉米減產(chǎn)，一方面可能因溫度升高，使生育期變短，致灌漿不足有關(guān)，另一方面也可能因溫度升高使水分散失，干旱風(fēng)險增大，影響玉米生長[20]。榆林玉米生長期內(nèi)，如何適應(yīng)氣候變暖成為玉米穩(wěn)產(chǎn)增產(chǎn)急需關(guān)注的問題。

圖4 1971～2012年榆林市玉米產(chǎn)量變化與生長期氣候因子變化的關(guān)系

2.2.2 不同發(fā)育階段各氣候因子變化對玉米產(chǎn)量的影響

基于榆林地區(qū)1971～2012年的玉米單產(chǎn)數(shù)據(jù)與玉米出苗-三葉、七葉-拔節(jié)、抽雄-開花、灌漿-乳熟、成熟5個關(guān)鍵生育期氣候因子，利用一階差分法建立逐年玉米單產(chǎn)變化值與各階段氣候因子變化值的對應(yīng)序列，采用相關(guān)分析、回歸分析研究不同發(fā)育階段各氣候因子對玉米產(chǎn)量的影響。

表2為5～9月逐月對應(yīng)的關(guān)鍵生育期平均溫度、平均最高溫度、平均最低溫度、降水量、日照時數(shù)變化與產(chǎn)量變化的相關(guān)系數(shù)表，由表2可見，榆林地區(qū)玉米產(chǎn)量與5個關(guān)鍵發(fā)育期中出苗-三葉期的平均溫度、日照時數(shù)、七葉-拔節(jié)期的平均溫度、平均最高溫度、抽雄-開花期的平均溫度、平均最高溫度呈極顯著負(fù)相關(guān)(P<0.01)，與出苗-三葉期平均最高溫度、灌漿-乳熟期平均溫度、平均最高溫度、成熟期降水量呈顯著負(fù)相關(guān)(P<0.05)； 與七葉-拔節(jié)期降水量呈極顯著正相關(guān)(P<0.01)，與成熟期平均最高溫度呈顯著正相關(guān)(P<0.05)。各發(fā)育階段的平均最低溫度與產(chǎn)量的相關(guān)性不明顯。

利用多元逐步回歸法分析標(biāo)準(zhǔn)化處理的玉米產(chǎn)量逐年變化數(shù)據(jù)和5個關(guān)鍵生育階段5個氣象因子共25個因子逐年變化值，分別用X1、X2、X3、X4、X5、…X25表示25個因子，構(gòu)建模型如式(5)所示，以進(jìn)一步分析多因子對玉米產(chǎn)量的綜合影響。

ΔY=0.385ΔX7-0.267ΔX19-0.331ΔX3-0.273ΔX14+0.890

(5)

其中,5個生育階段分別為出苗-三葉、七葉-拔節(jié)、抽雄-開花、灌漿-乳熟、成熟期,5個氣象因子分別為平均溫度、平均最高溫度、平均最低溫度、降水量、日照時數(shù)。多元逐步回歸模型中，ΔY代表玉米逐年產(chǎn)量差，ΔX7表示七葉-拔節(jié)期降水變化值，ΔX19表示灌漿-乳熟期平均最高溫度變化值，ΔX3代表出苗-三葉期日照變化值，ΔX14代表抽雄開花期的平均最高溫度變化值。F=11.235，P<0.01，模型達(dá)到極顯著。

由式(5)可知，玉米各發(fā)育階段中，拔節(jié)期的降水量，灌漿-乳熟期、抽雄-開花期平均最高溫度、苗期的日照是影響玉米產(chǎn)量的主導(dǎo)因子。R2=0.56，說明這些因子變化對玉米產(chǎn)量的影響貢獻(xiàn)可達(dá)56%。

綜合各發(fā)育階段氣候因子變化與產(chǎn)量的關(guān)系，榆林地區(qū)玉米從出苗到乳熟期，平均溫度、平均最高溫度變化對產(chǎn)量的影響不容忽視，特別是抽雄-開花期，正值榆林高溫頻發(fā)期，持續(xù)高溫極易導(dǎo)致玉米“卡脖旱”，影響玉米授粉； 灌漿-乳熟期的高溫會影響?zhàn)B分運輸，使顆粒瘦秕，粒重減輕； 七葉-拔節(jié)期是玉米需水關(guān)鍵期，而此時榆林進(jìn)入初夏，如苗期蓄水不足，易產(chǎn)生春夏連旱、初夏干旱，應(yīng)適時灌溉，以滿足玉米生長的水分需求； 玉米成熟期，籽粒干物質(zhì)累積對熱量需求高，而此時期榆林多連陰雨、秋淋天氣，會間接導(dǎo)致溫度、日照下降，影響籽粒質(zhì)量。關(guān)注各發(fā)育階段的氣象因子變化，應(yīng)對并減緩各關(guān)鍵期的不利因子對玉米種植的負(fù)面影響是玉米種植管理需重點關(guān)注的問題。

采用2013、2014年氣象數(shù)據(jù)及產(chǎn)量數(shù)據(jù)驗證模型準(zhǔn)確性，結(jié)果及誤差見表3?？梢姌?gòu)建回歸模型檢驗2013、2014年玉米產(chǎn)量的誤差分別為5.1%、4.5%，但由于玉米產(chǎn)量變化的影響因素眾多，該文僅針對氣候變化宏觀影響進(jìn)行探討，對氣象災(zāi)害、耕作變化等的影響未能考慮在內(nèi)，有待于今后更進(jìn)一步研究。

表2 玉米產(chǎn)量變化與各發(fā)育期氣候因子變化的相關(guān)分析

發(fā)育期平均溫度平均最高溫度平均最低溫度降水量日照時數(shù)出苗-三葉期-0.446**-0.377*-0.3020.135-0.433**七葉-拔節(jié)期-0.438**-0.446**-0.2700.539**-0.294抽雄-開花期-0.493**-0.546**-0.1950.094-0.304灌漿-乳熟期-0.359*-0.342*-0.2160.029-0.30成熟期0.2580.396*-0.046-0.368*0.144 注:*表示通過0.05顯著性檢驗,**表示通過0.01顯著性檢驗

表3 回歸模型精度檢驗

年份ΔX7(mm)ΔX19(℃)ΔX3(h)ΔX14(℃)ΔY(kg/hm2)預(yù)測產(chǎn)量(kg/hm2)實際產(chǎn)量(kg/hm2)誤差(%)2013-20.5-1.7-34.4-3.81130737870205.12014-28.5-1.836.41.6598761872904.5

3 結(jié)論與討論

(1)1971～2012年，毛烏素沙漠與黃土高原過渡帶玉米生長期內(nèi)平均溫度、平均最高溫度、平均最低溫度等氣候因子年際變化均呈現(xiàn)極顯著的上升趨勢，上升幅度分別為0.24℃/10a、0.26℃/10a和0.32℃/10a，在空間上表現(xiàn)為升溫速率由東南向西北遞增。生長期內(nèi)降水量有緩慢上升趨勢，總?cè)照諘r數(shù)有下降趨勢，變化不明顯。

(2)過去42年，玉米產(chǎn)量與生長期各氣候因子呈極顯著負(fù)相關(guān)的是平均溫度、平均最高溫度、積溫，與平均溫度的相關(guān)系數(shù)最高，平均溫度每升高1℃，造成產(chǎn)量減少835.5kg/hm2，其次為平均最高溫度，平均最高溫度每升高1℃，產(chǎn)量減少617.1kg/hm2。玉米產(chǎn)量與生長期降水量變化呈正相關(guān)，降水量每增加100mm，使玉米產(chǎn)量增加364.7kg/hm2； 與日照時數(shù)變化呈負(fù)相關(guān)，日照時數(shù)每增加100h，玉米減產(chǎn)325.8kg/hm2。

(3)在各個發(fā)育階段，玉米產(chǎn)量與出苗-三葉期平均溫度、日照、與七葉-拔節(jié)期平均溫度、最高溫度、與抽雄-開花期平均溫度、最高溫度均呈極顯著負(fù)相關(guān)，與七葉-拔節(jié)期降水量呈極顯著的正相關(guān)，其次與成熟期平均最高溫度顯著正相關(guān)。

(4)各發(fā)育階段的各氣候因子變化綜合作用時，苗期日照、拔節(jié)期降水量，抽雄開花期的平均最高溫度和灌漿-乳熟期的平均最高溫度是影響玉米產(chǎn)量的主要因子，這些因子變化對玉米產(chǎn)量的影響貢獻(xiàn)可達(dá)56%。用2013、2014年的實際產(chǎn)量和氣象因子數(shù)據(jù)驗證模型效果良好。

全球氣候變化背景下，榆林玉米生長期氣候變化與全球氣候變暖趨勢一致，榆林作為西北干旱區(qū)，降水本來就少，溫度升高可致干旱加劇，自然降水更難以補(bǔ)足作物生長所需水分，對玉米產(chǎn)量安全提出挑戰(zhàn)。特別是抽雄-開花、灌漿-乳熟期等生育關(guān)鍵期缺少水分供給可造成玉米雄穗敗育，粒重減輕，而成熟期降水過多不利于干物質(zhì)積累，造成產(chǎn)量下降[21]。盡管施肥能促進(jìn)增產(chǎn)，但也只有在降水量充足的情況下才能有明顯的增產(chǎn)成效[22]，可從一定程度上降低種植成本[23]。靳英華等人研究指出，早播可以促進(jìn)玉米根系生長[24]，使根系扎入墑情較好層次。不同播種期能給玉米生長配置特定的氣候條件組合，段金省研究指出，灌漿-乳熟期若與一年中降水最多時段吻合，有助于提高灌漿效果[25]。輻射增多雖然能促進(jìn)有機(jī)物質(zhì)轉(zhuǎn)移并增產(chǎn)，但榆林光照資源豐富，輻射過強(qiáng)時引起蒸騰加劇，容易形成高溫干旱使玉米減產(chǎn)。因此應(yīng)對氣候變暖，避免干旱風(fēng)險，是榆林玉米生產(chǎn)需重點關(guān)注的問題，對穩(wěn)定玉米生產(chǎn)尤為重要。

為有效應(yīng)對氣候變暖，榆林地區(qū)一方面應(yīng)盡可能順應(yīng)自然，積極研究、嘗試引進(jìn)抗逆抗旱及晚熟品種，另一方面應(yīng)加強(qiáng)水利建設(shè)，擇時灌溉，開展人工增雨彌補(bǔ)玉米生長所需水分。同時更要加強(qiáng)田間管理技術(shù)研究，根據(jù)溫度、降水的監(jiān)測和預(yù)報信息，選擇合理播期，為玉米種植贏得有利氣候條件，配以播前催芽、合理的種植密度、覆膜等栽培技術(shù)，可規(guī)避低溫冷害、高溫、干旱風(fēng)險。

[1]IntergovernmentalPanelOnClimateChange.FifthAssessmentReport(AR5),2014-07-08

[2] 徐雨晴,苗秋菊，沈永平.2008年：氣候持續(xù)變暖，極端事件頻發(fā).氣候變化研究進(jìn)展,2009, 5(1): 56～60

[3] 劉穎杰,林而達(dá).氣候變暖對中國不同地區(qū)農(nóng)業(yè)的影響.氣候變化研究進(jìn)展,2007, 3(4): 229～233

[4] 趙鴻,肖國舉，王潤元，等.氣候變化對半干旱雨養(yǎng)農(nóng)業(yè)區(qū)春小麥生長的影響.地球科學(xué)進(jìn)展,2007, 22(3): 322～327

[5] 秦大河. 氣候變化的事實與影響及對策.中國科學(xué)基金,2003,17(1)： 1～3.59～63

[6] 王柳,熊偉，溫小樂，等.溫度降水等氣候因子變化對中國玉米產(chǎn)量的影響.農(nóng)業(yè)工程學(xué)報,2014, 30(21): 138～146

[7] 呂碩,楊曉光，趙錦，等.氣候變化和品種更替對東北地區(qū)春玉米產(chǎn)量潛力的影響.農(nóng)業(yè)工程學(xué)報,2013,29(18)： 17～190

[8] 侯艷林,姜文來.氣候變化背景下主要氣象因素對作物單產(chǎn)影響的模擬研究——以吉林公主嶺玉米為例.中國農(nóng)業(yè)資源與區(qū)劃,2014, 35(6): 58～64

[9] 李言照,東先旺，劉光亮，等.光溫因子對玉米產(chǎn)量及產(chǎn)量構(gòu)成因素值的影響.中國生態(tài)農(nóng)業(yè)學(xué)報,2002, 10(2): 86～89

[10]吉相臣,劉昂，康玉俠，等.鐵嶺地區(qū)降水量與玉米產(chǎn)量相關(guān)分析及預(yù)測模型研究.現(xiàn)代農(nóng)業(yè)科技,2012，(5)： 306～308,313

[11]杜軍,劉英.氣候變化對隴東塬區(qū)界限溫度和作物種植結(jié)構(gòu)的影響分析.中國農(nóng)學(xué)通報,2015, 31(27): 253～257

[12]屈洋,孟慶立，趙寧娟，等.岐山關(guān)中西部半濕潤區(qū)氣候變化對玉米產(chǎn)量的影響.種子世界,2014,(10)： 31～33

[13]扈艷萍,曹敏建，劉敏.遼寧省玉米主產(chǎn)區(qū)氣候因子與玉米產(chǎn)量相關(guān)性的研究.玉米科學(xué),2008, 16(3): 140～146

[14]馬雅麗,王志偉，欒青，等.玉米產(chǎn)量與生態(tài)氣候因子的關(guān)系.中國農(nóng)業(yè)氣象,2009, 30(4): 565～568

[15]劉曉瓊,劉彥隨，延軍平，等.生態(tài)脆弱區(qū)多年氣候變化特征分析.干旱區(qū)資源與環(huán)境,2008, 22(1): 54～59

[16]高娟,馬鋒，曹梅，等.榆林市近42年農(nóng)業(yè)界限溫度變化特征分析.陜西氣象, 2015，(3)： 32～37

[17]史文嬌,陶福祿，張朝，等.基于統(tǒng)計模型識別氣候變化對農(nóng)業(yè)產(chǎn)量貢獻(xiàn)的研究進(jìn)展.地理學(xué)報,2012,67(9)： 1213～1222

[18]LiuY，WangE，YangXG，etal.ContributionsofclimaticandcropvarietalchangestocropproductionintheNorthChinaPlain，since1980s.GlobalChangeBiology,2010, 16(1): 2289～2299

[19]KucharikCJ，SerbinSP.ImpactsofrecentclimatechangeonWisconsincornandsoybeanyieldtrend.EnvironmentalResearchLetters,2008,3： 034003

[20]王春春,黃山，鄧艾興，等.東北雨養(yǎng)農(nóng)區(qū)氣候變暖趨勢與春玉米產(chǎn)量變化的關(guān)系分析.玉米科學(xué),2010,18(6)： 64～68

[21]葛均筑. 氣象資源特性對玉米產(chǎn)量形成的影響及長江中游玉米高產(chǎn)關(guān)鍵技術(shù)研究.武漢：華中農(nóng)業(yè)大學(xué),2015

[22]冶明珠,郭建平，袁彬，等.氣候變化背景下東北地區(qū)熱量資源及玉米溫度適宜度.應(yīng)用生態(tài)學(xué)報,2012, 23(10): 2786～2794

[23]范成方,史建民.山東省糧食種植成本影響因素的實證分析——以玉米、小麥為例.中國農(nóng)業(yè)資源與區(qū)劃,2014, 35(2): 67～74

[24]靳英華,周道瑋，秦麗杰，等.適應(yīng)氣候變化的吉林省半干旱區(qū)玉米播種期.應(yīng)用生態(tài)學(xué)報,2012, 23(10): 2795～2802

[25]段金省,牛國強(qiáng).氣候變化對隴東塬區(qū)玉米播種期的影響.干旱地區(qū)農(nóng)業(yè)研究,2007, 25(2): 235～238

EFFECT OF CLIMATE CHANGE ON CORN YIELD IN YULIN, NORTH SHAANXI PROVINCE OF CHINA*

Gao Juan1，F(xiàn)eng Lingzhi1，Zhang Jiankang1，Wang Yun1, Cao Mei2，Chang Shangang1

(1.Yulin Meteorological Bureau of Shaanxi province，Yulin 719000, China； 2.Xi′an Meteorological Bureau of Shaanxi province，Xi′an 710016, China)

The study analyzed the effects of climate change on corn yield in Yulin, the typical ecotone between Mu Us desert and Loess Plateau in north Shaanxi province of China, based on the meteorological data and statistics data from 1971 to 2012. The methods of first order difference, correlation analysis and multiple liner regression were used. The results showed that the heat resource for whole corn growth period increased significantly during the past 42a. The average temperature, average maximum temperature and mean minimum temperature increased by 0.24℃/10a, 0.26℃/10a and 0.32℃/10a, respectively. The changes of precipitation and sunlight hours were not as obvious as the temperature during the corn growth period. The average temperature, average maximum temperature, accumulated temperature had significant negative correlation with maize yield.The correlation coefficient between maize yield and the average temperature was the biggest. Whenthe average temperature increased by 1℃, maize yield would decrease by 835.5 kg/hm2; However, the precipitation increased by 100mm, maize yield would increase by 364.7kg/hm2;and when the sunlight hours increased by 100h, maize yield would decrease by 325.8kg/hm2. For the each growth stage, the average maximum temperature, average temperature from three leaves stage to milking stage, and the sunlight hours at seedling stage had negative influence on the yield of maize. The precipitation at the seven leaves-jointing stage and the average maximum temperature at the maturity showed positive effect on maize yield. The average maximum temperature during pumping male flowering and maturity stage, the precipitation of jointing stage, and the sunlight hours during the seedling-three leaves stage had significant influence on the maize yield.The model can explain the 56% change of maize yield in the typical ecotone between Mu Us desert and Loess Plateau.

Mu Us Desert; Loess Plateau; first order difference; climate change; corn yield

10.7621/cjarrp.1005-9121.20161218

2016-03-10

高娟(1980—)，女，陜西榆林人，碩士、高級工程師、副主任。研究方向：氣象為農(nóng)服務(wù)管理及氣象減災(zāi)服務(wù)工作。Email： 158749792@qq.com

*資助項目：榆林市科技局2015年技術(shù)研究與發(fā)展計劃項目“榆林氣候變化對玉米產(chǎn)量的影響及增產(chǎn)技術(shù)研究”(2015JH-Z0)

S16; S513

A

1005-9121[2016]12-0118-08