寧夏馬鈴薯產(chǎn)量形成的關(guān)鍵影響氣象因素分析

鄭 方,徐 青,馬力文,張學(xué)藝

(中國氣象局旱區(qū)特色農(nóng)業(yè)氣象災(zāi)害監(jiān)測預(yù)警與風(fēng)險(xiǎn)管理重點(diǎn)實(shí)驗(yàn)室/寧夏氣象防災(zāi)減災(zāi)重點(diǎn)實(shí)驗(yàn)室,寧夏 銀川 750002)

馬鈴薯是寧夏中南部山區(qū)面積最大的農(nóng)作物和戰(zhàn)略性主導(dǎo)產(chǎn)業(yè)，種植面積近230 Khm2，總產(chǎn)量42.5萬t。該區(qū)域處于我國西北農(nóng)牧交錯(cuò)帶的黃土高原，丘陵溝壑縱橫，土層深厚，通透性好，十分適合種植馬鈴薯，旱作馬鈴薯是貧困山區(qū)農(nóng)民抵御干旱的重要農(nóng)產(chǎn)品。由于該區(qū)域十年九旱，災(zāi)害頻繁，農(nóng)業(yè)生產(chǎn)受氣候影響很大[1]，旱作馬鈴薯產(chǎn)量低而不穩(wěn)，年際波動(dòng)大。2016年之后寧夏把馬鈴薯主食化作為脫貧濟(jì)困、促進(jìn)增收的重大舉措，大力發(fā)展馬鈴薯特色產(chǎn)業(yè)。寧夏馬鈴薯播種面積和總產(chǎn)在逐年提高，但單產(chǎn)年際間波動(dòng)劇烈。中部干旱帶1981-2000年馬鈴薯單產(chǎn)由4 000 kg·hm-2的水平上升至9 000 kg·hm-2的水平，翻了1.25倍，但自21世紀(jì)初開始下降至目前的6 300 kg·hm-2的水平；南部山區(qū)1981-2005年馬鈴薯單產(chǎn)由5 000 kg·hm-2的水平逐漸上升至13 000 kg·hm-2的水平，單產(chǎn)翻了1.6倍，近5 a來略下降至12 000 kg·hm-2的水平。中部干旱帶的馬鈴薯的生長受氣候條件的影響更大，單產(chǎn)增產(chǎn)幅度比南部山區(qū)低，特別是近10 a來呈減產(chǎn)趨勢，干旱造成單產(chǎn)的區(qū)域差距大，產(chǎn)量低而不穩(wěn)。為此，筆者研究在寧夏西海固貧困區(qū)開展大田分期播種試驗(yàn)和適地播種，分析氣溫、降水等氣象因子對(duì)寧夏貧困區(qū)旱地馬鈴薯產(chǎn)量的影響，確定該地區(qū)影響產(chǎn)量的關(guān)鍵氣象因子，提出適應(yīng)氣候變化的對(duì)策。

氣象條件對(duì)作物的影響是一個(gè)復(fù)雜的過程，某氣象要素在特定的區(qū)域總是能滿足作物生長時(shí)，該地該要素對(duì)產(chǎn)量的影響就變得不太重要。當(dāng)某一氣象要素作為資源，在該區(qū)域變異性較大且總是不太能滿足作物生長時(shí)，這個(gè)氣象因子就變得非常重要，我們就稱之為關(guān)鍵因子，或在作物的某一發(fā)育階段總是不能滿足作物生長，這個(gè)發(fā)育期便被稱為關(guān)鍵期。因此每一個(gè)地區(qū)，每一種作物的關(guān)鍵氣象因子和關(guān)鍵發(fā)育期并不完全相同。如河南冬小麥返青-抽穗期≥10℃有效積溫相對(duì)于其他形式的熱量資源波動(dòng)變化大，與產(chǎn)量的相關(guān)系數(shù)高，是冬小麥生長的關(guān)鍵氣象因子[2]。在陜西省寶雞地區(qū)，冬前10-12月和次年5月的氣溫對(duì)冬小麥產(chǎn)量形成起著至關(guān)重要的作用[3]。種植在不同地區(qū)的馬鈴薯關(guān)鍵期和關(guān)鍵氣象因子各有不同。李華鵬等在四川臥龍地區(qū)的試驗(yàn)研究表明，在塊莖增長期，較高的氣溫對(duì)馬鈴薯塊莖增大有利[4]。吳炫柯等通過分期播種試驗(yàn)發(fā)現(xiàn)，在廣西柳州馬鈴薯結(jié)薯-成熟期是馬鈴薯關(guān)鍵生育期，期間的氣象因子對(duì)產(chǎn)量影響較大[5]。在青海馬鈴薯全生育期內(nèi)氣溫對(duì)其產(chǎn)量的影響并不大，而生育期內(nèi)降水對(duì)馬鈴薯產(chǎn)量的影響十分明顯[6]。在遼寧沈陽地區(qū)，年降水量增加，馬鈴薯產(chǎn)量呈增加趨勢，而年平均氣溫減少則造成馬鈴薯減產(chǎn)[7]。朱赟赟等認(rèn)為雖然在寧夏不同地區(qū)影響馬鈴薯產(chǎn)量的氣象因子不同，但降水量始終是影響馬鈴薯產(chǎn)量的第一主要因子[8]。

筆者研究以歷年各市縣馬鈴薯種植面積為權(quán)重進(jìn)行氣象因子合成，利用氣象適宜度和作物水分虧缺和產(chǎn)量進(jìn)行擬合，找出產(chǎn)量擬合的關(guān)鍵影響氣象因素，為農(nóng)業(yè)生產(chǎn)提供支持。

1 資料和方法

1.1 資料來源

產(chǎn)量、面積資料來源于1961-2015年寧夏中南部山區(qū)各市縣馬鈴薯的產(chǎn)量調(diào)查數(shù)據(jù)。氣象資料源于歷年各氣象觀測站歷史觀測記錄，并根據(jù)本研究的方法將其進(jìn)行面積合成。

1.2 氣象資料面積權(quán)重合成

將區(qū)域內(nèi)多個(gè)氣象觀測站的歷年氣象資料，以歷年各市縣馬鈴薯種植面積為權(quán)重進(jìn)行合成來反映馬鈴薯生長氣象條件：

(1)

式中，Mj為j日綜合權(quán)重的氣象要素值(主要為日平均氣溫、日降水量、日照時(shí)數(shù)、日最高氣溫、日最低氣溫等)，Mij為i站j日的氣象要素，Si為i縣當(dāng)年馬鈴薯種植面積。按照該年各縣馬鈴薯種植面積分別合成當(dāng)年的氣象資料，作為當(dāng)年的區(qū)域綜合氣象資料。按照上述方法，將1954-2016年的逐站氣象資料合成為馬鈴薯區(qū)域綜合氣象資料。

1.3 相對(duì)氣象產(chǎn)量

為明確氣象因子的作用，需濾除品種、灌溉條件改善等生產(chǎn)力水平提高的趨勢產(chǎn)量，提取受氣候影響的產(chǎn)量波動(dòng)。筆者以單產(chǎn)與趨勢產(chǎn)量的比值表示相對(duì)氣象產(chǎn)量，該比值大于1表示氣象產(chǎn)量高于該年平均產(chǎn)量，即增產(chǎn)年，比值小于1則當(dāng)年產(chǎn)量在趨勢產(chǎn)量以下，為減產(chǎn)年。相對(duì)氣象產(chǎn)量表示為：

y/yt=1+yW/yt

(2)

參考山區(qū)各地固海揚(yáng)水工程、鹽環(huán)定揚(yáng)黃灌區(qū)建設(shè)發(fā)揮效用的時(shí)間，以1976年為界限，分段模擬了馬鈴薯的趨勢產(chǎn)量(圖1)。

(3)

1.4 作物水分滿足率和水分盈虧量

采用彭曼公式計(jì)算參考作物蒸散量，公式如下：

(4)

式中，T為日平均氣溫，Rn為冠層凈輻射(MJm-2d-1)，G為土壤熱通量(MJm-2d-1)，U2為2 m處平均風(fēng)速(ms-1)，VPD為水氣壓差(KPa),Δ為飽和水氣壓曲線斜率(Kpa℃-1)，γ為干濕表常數(shù)。

參考FAO提供的不同生育階段作物Kc系數(shù)，估算歷年馬鈴薯不同生育期需水量。作物生長第i日的需水量ETci為：

ETci=Kci×EToi

(5)

其中Kci為作物系數(shù)。FAO推薦的馬鈴薯作物的Kc系數(shù)見表1。

表1 主要作物的Kc值(FAO)

農(nóng)田蒸散量ET由作物蒸騰和土壤蒸發(fā)構(gòu)成，在干旱條件下，可忽略水分下滲量和毛管水上升量，其簡易平衡方程為：

ET=V1-V2+(1-r)∑R+P

(5)

式中，V1、V2分別為某時(shí)段開始和結(jié)束時(shí)一定深度土層內(nèi)的土壤貯水量[mm]，可由土壤重量含水量與土壤容重的乘積來確定。r為平均徑流系數(shù)，∑R為某時(shí)段累計(jì)降水量，P為灌溉量。

水分滿足率是某時(shí)段內(nèi)作物耗水量與需水量的比值，其計(jì)算公式為：

(6)

水分盈虧量是某時(shí)段內(nèi)作物耗水量與需水量的差值，其計(jì)算公式為：

βK=ET-ETC

(7)

(6)、(7)式中，βs、βk分別為水分滿足率和水分盈虧，ET為作物耗水量，ETc為作物需水量。

1.5 水分利用效率(WUE)

水分利用效率是評(píng)價(jià)水分虧缺下植物生長適宜度的綜合指標(biāo)之一。其計(jì)算方法為[9]：

式中：Yd為經(jīng)濟(jì)產(chǎn)量，ETc為作物需水量。

1.6 氣象適宜度

作物氣候適宜度是農(nóng)作物在各生長階段內(nèi)所需要水分、溫度、光照等氣象條件的滿足程度，在評(píng)價(jià)作物生長氣象條件優(yōu)劣等方面有獨(dú)特的作用，評(píng)價(jià)的客觀性和科學(xué)性較強(qiáng)。利用氣候適宜度以及寧夏中南部山區(qū)馬鈴薯生長氣象條件適宜性，建立產(chǎn)量氣象模型，可評(píng)價(jià)氣象條件對(duì)中南部山區(qū)馬鈴薯生長發(fā)育、產(chǎn)量形成的適宜性。筆者研究溫度適宜度函數(shù)的計(jì)算采用馬樹慶等基于溫度三基點(diǎn)溫度不對(duì)稱提出的溫度適宜度函數(shù)[10]，其計(jì)算公式如下：

(9)

式中t1為馬鈴薯某發(fā)育時(shí)段的下限溫度；t2為上限溫度；t0為馬鈴薯發(fā)育的最適溫度；t為某旬的實(shí)際溫度。馬鈴薯的三基點(diǎn)溫度指標(biāo)如表2。

表2 馬鈴薯不同生育期的三基點(diǎn)溫度

改進(jìn)魏瑞江等在旱作農(nóng)田的降水適宜度公式[11]，以上旬降水量在本旬的存留比例與本旬降水量的累加值占本旬需水量的70%和2倍為界，如果本旬水分收入項(xiàng)少于該作物旬需水量的70%，適宜度為降水量除以需水量，如果高于本旬作物需水量的2倍，則水分過多對(duì)作物生長產(chǎn)生負(fù)作用，用作物需水量除以水分收入項(xiàng)，在需水量的70%至2倍的范圍內(nèi)，適宜度為1，即：

(10)

式中Ri-1j、Ri分別為上旬和本旬降水量，r為存留系數(shù)，即上旬降水用于本旬的量。根據(jù)利用固原農(nóng)業(yè)氣象試驗(yàn)站土壤水分監(jiān)測與估算結(jié)果的比較，統(tǒng)一用0.5，效果較好。ETci為利用彭曼公式與FAO推薦的馬鈴薯作物系數(shù)推算的需水量[12]。

大多數(shù)文獻(xiàn)研究認(rèn)為，日照時(shí)數(shù)達(dá)到可照時(shí)數(shù)的70%(日照百分率)為臨界點(diǎn)[13],其計(jì)算公式如下：

(11)

S為某階段日照時(shí)數(shù)，S0為太陽可照時(shí)數(shù)，可通過地理緯度、太陽赤緯計(jì)算，b為系數(shù)，可參照相關(guān)研究結(jié)果。

在確定最適宜的光溫水單因子氣候適宜度模型后，可利用立方根法可建立光溫水綜合氣候適宜度：

(12)

式中，CS為氣候適宜度指數(shù)，CS(Ta)、CS(Pr)、CS(S)分別為溫度、水分、光照適宜度。

2 結(jié)果

2.1 影響相對(duì)氣象產(chǎn)量的因子相關(guān)分析

根據(jù)固原馬鈴薯發(fā)育期觀測，一般在4月下旬至5月上旬播種，6月中旬現(xiàn)蕾，6月下旬至7月中旬為花序形成、開花階段，7月下旬至9月底為塊莖膨大階段，9月下旬末或10月初采挖。由于各地馬鈴薯播種時(shí)間前后可達(dá)一個(gè)月，收獲期有些年份可延遲到10月下旬，為了比較不同年份馬鈴薯生育氣象條件，筆者研究統(tǒng)一將1/4～20/6期間定為苗期、21/6～20/7定為花序形成至開花期、21/7～30/9定為塊莖膨大期。以1954-2014年的馬鈴薯相對(duì)氣象產(chǎn)量與不同生育時(shí)段內(nèi)的溫度、降水、日照適宜度、水分盈虧量和水分滿足率等進(jìn)行相關(guān)分析。氣象因子與相對(duì)氣象產(chǎn)量的相關(guān)篩選結(jié)果見表3。

表3 馬鈴薯相對(duì)氣象產(chǎn)量與氣象因子相關(guān)分析(1954-2014年)

2.2 水分滿足率及盈虧量與產(chǎn)量的關(guān)系

盡管水分盈虧量和水分滿足度與氣象產(chǎn)量的相關(guān)因子并沒有通過0.05的檢驗(yàn)，但其與相對(duì)氣象產(chǎn)量的擬合方程的復(fù)相關(guān)系數(shù)依然通過了R0.05的顯著性檢驗(yàn)(圖2a、b)，苗期水分滿足率一般在30%～150%之間，61 a中有47 a在100%以下，水分盈虧以負(fù)值居多，表示苗期水分多數(shù)年份得不到滿足。

花序形成至開花期水分滿足率和盈虧量與產(chǎn)量的擬合方程的復(fù)相關(guān)系數(shù)也通過R0.001的顯著性檢驗(yàn)(圖3a、b)。在61 a中，完全能夠滿足馬鈴薯生長的水分條件的僅有2 a，其余59 a水分滿足率在-5.8%～-92.4%之間，多年平均狀況在-60%左右，水分盈虧量在-164～-10 mm，多年平均缺水60 mm以上，表明絕大多數(shù)年份6月下旬至7月中旬降水量嚴(yán)重不足，基本不能滿足馬鈴薯水分需求，嚴(yán)重限制了產(chǎn)量的提高，需要在此階段補(bǔ)充灌溉才能大幅度提高產(chǎn)量。

塊莖膨大期水分滿足率和盈虧量與產(chǎn)量的關(guān)系不顯著(圖4a、b)，表明后期馬鈴薯不缺水，水分能夠滿足膨大需求，有些年份雨水過多，引起晚疫病而減產(chǎn)。61 a中，全生育期水分能夠安全滿足需水的僅2 a，水分滿足率一般在36%～92%之間，缺水量一般在20～370 mm之間。

2.3 水分利用效率(MUE)與產(chǎn)量的關(guān)系

馬鈴薯全生育期的WUE與產(chǎn)量關(guān)系極其密切，相對(duì)氣象產(chǎn)量隨著WUE的增大而升高，通過R0.001的顯著性檢驗(yàn)(圖5)。通過不同時(shí)段降水量與WUE的相關(guān)普查，發(fā)現(xiàn)馬鈴薯薯塊膨大期降水量與WUE的負(fù)相關(guān)達(dá)到極顯著，通過了R0.001的顯著性檢驗(yàn)(圖6)，說明膨大期降水量對(duì)產(chǎn)量的影響大。究其原因是馬鈴薯生長后期出現(xiàn)持續(xù)陰雨寡照的年份很多，雨水過多往往造成馬鈴薯晚疫病嚴(yán)重發(fā)生，特別是近5 a晚疫病大面積發(fā)生，田間病株率幾乎100%，導(dǎo)致秋雨偏好的年份馬鈴薯反而減產(chǎn)，甚至嚴(yán)重減產(chǎn)。

2.4 氣象適宜度與產(chǎn)量的關(guān)系

馬鈴薯相對(duì)氣象產(chǎn)量與5月、9月份的溫度適宜度為明顯正相關(guān)，播種至出苗階段、塊莖膨大后期溫度往往偏低，兩頭熱量不足限制了產(chǎn)量的提高(圖7a、7b)。全生育期溫度累積適宜度指數(shù)與產(chǎn)量呈2次曲線關(guān)系，最適值為158，熱量不足或過多產(chǎn)量均會(huì)降低(圖7c)。說明雖然適宜程度高，對(duì)馬鈴薯生長越有有利，但不一定對(duì)產(chǎn)量形成有利。

6月下旬至8月上旬是馬鈴薯花序形成、開花和塊莖膨大初期，對(duì)水分敏感，其與產(chǎn)量的決定系數(shù)為0.1846，達(dá)到R0.001的極顯著水平，此時(shí)氣溫處于全年最高時(shí)段，正值馬鈴薯開花至塊莖膨大初期的水分敏感期，雨水偏少影響小薯膨大數(shù)量，干旱抑制塊莖膨大，對(duì)結(jié)薯量、產(chǎn)量影響很大(圖8a)。薯塊膨大期水分適宜度與產(chǎn)量呈2次曲線關(guān)系，分別通過了R0.001的顯著性檢驗(yàn)，薯塊膨大期是馬鈴薯干物質(zhì)積累的重要階段，此時(shí)水分缺少影響產(chǎn)量形成，但水分適宜度高往往意味著陰雨天多影響馬鈴薯光合積累造成產(chǎn)量降低(圖8b)。

8月中旬至9月下旬的塊莖膨大時(shí)段日照適宜度與產(chǎn)量正相關(guān)(圖9a)，說明膨大期日照不足往往與陰雨寡照聯(lián)系緊密，多雨造成晚疫病危害嚴(yán)重。從7-9月份和光溫水綜合氣象適宜度與產(chǎn)量呈正相關(guān)關(guān)系，通過了R0.001的顯著性檢驗(yàn)，隨著綜合累積適宜度指數(shù)的增大，產(chǎn)量提高。

3 結(jié)論與討論

寧夏中南部山區(qū)馬鈴薯產(chǎn)量與氣象條件關(guān)系緊密，受自然條件限制，這些地區(qū)以旱作栽培為主，農(nóng)業(yè)生產(chǎn)與降水多寡和降水出現(xiàn)的季節(jié)息息相關(guān)，光熱資源是進(jìn)一步提高產(chǎn)量的基礎(chǔ)。限制馬鈴薯產(chǎn)量的關(guān)鍵氣象因素有苗期的溫度，花序形成-開花期為馬鈴薯生長發(fā)育的關(guān)鍵時(shí)期，水分適宜度、日照適宜度和綜合適宜度均對(duì)產(chǎn)量有影響，說明此階段馬鈴薯對(duì)氣象條件敏感，特別水對(duì)水分敏感，是馬鈴薯的水分關(guān)鍵期。薯塊膨大期是馬鈴薯另一個(gè)關(guān)鍵期，溫度、降水量和光照條件對(duì)產(chǎn)量均有影響，其中以日照尤為顯著。

從馬鈴薯實(shí)際生產(chǎn)調(diào)查中也可發(fā)現(xiàn)，苗期寧夏中南部山區(qū)馬鈴薯需水量較小，水分資源不是限制因子，但熱量資源略顯不夠，所以苗期氣溫為產(chǎn)量形成的關(guān)鍵因子。花序形成-開花期，隨著氣溫升高，熱量資源已不是限制因子了，此時(shí)馬鈴薯進(jìn)入了需水旺盛期，自然降水嚴(yán)重缺乏，此階段水分是關(guān)鍵因子，以開花前后的降水量影響最大。薯塊膨大期寧夏中南部山區(qū)往往出現(xiàn)低溫、陰雨寡照天氣，誘發(fā)嚴(yán)重的晚疫病，導(dǎo)致減產(chǎn)。苗期在水分適宜條件下，溫度偏高有利于早發(fā)，使開花期提早，有利于避開花期高溫少雨的不利氣象條件，從而有利于增產(chǎn)。