不同生育期干旱對(duì)玉米生長(zhǎng)及產(chǎn)量的影響模擬

黃巖 李晶 王瑩 張曉月 李雨鴻 焦敏 周斌 王賀然

摘要 利用作物生長(zhǎng)模型模擬苗期、拔節(jié)期和灌漿期發(fā)生不同時(shí)長(zhǎng)干旱（用連續(xù)無(wú)雨日數(shù)表示）時(shí)玉米的生長(zhǎng)過(guò)程，分析單一發(fā)育期干旱發(fā)生干旱對(duì)玉米產(chǎn)量的影響。結(jié)果表明，苗期發(fā)生5～40 d的持續(xù)干旱，對(duì)玉米籽粒灌漿無(wú)明顯影響，最終產(chǎn)量無(wú)明顯變化。拔節(jié)期或灌漿期分別發(fā)生5～40 d的持續(xù)干旱，均會(huì)對(duì)玉米籽粒灌漿產(chǎn)生負(fù)面影響，并最終導(dǎo)致產(chǎn)量下降，且干旱持續(xù)時(shí)間越長(zhǎng)對(duì)灌漿的影響越早，減產(chǎn)越嚴(yán)重。當(dāng)玉米拔節(jié)期持續(xù)干旱5～10 d時(shí)，2013年和2016年的最終產(chǎn)量均無(wú)明顯變化;持續(xù)干旱15～20 d時(shí)，2013年減產(chǎn)10%以內(nèi)，2016年減產(chǎn)20%以內(nèi)。持續(xù)干旱25～40 d時(shí)，2013年減產(chǎn)15%～62%，2016年減產(chǎn)30%～67%。當(dāng)玉米灌漿期持續(xù)干旱5～30 d時(shí)，2013年玉米產(chǎn)量基本無(wú)明顯影響;當(dāng)玉米灌漿期持續(xù)干旱5～10 d時(shí)，2016年玉米產(chǎn)量基本無(wú)明顯影響;當(dāng)玉米灌漿期持續(xù)干旱15～30 d時(shí)，2016年玉米減產(chǎn)1.8%～8.20%?？梢?jiàn)，拔節(jié)期同等程度干旱造成的玉米減產(chǎn)幅度均大于苗期和灌漿期干旱。

關(guān)鍵詞 干旱;作物模型;產(chǎn)量;影響模擬

中圖分類號(hào)：S513 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：2095-3305（2019）06-047-03

DOI： 10.19383/j.cnki.nyzhyj.2019.06.017

Impact Simulation of Drought on Maize Growth and Yield in Different Growth Stages

HUANG Yanet al（Ecological Meteorology and Satellite Remote Sensing Center of Liaoning Province，Shenyang，Liaoning 110166）

Abstract In order to explore the impacts of drought occurred in single development stage respectively on the growth formation and yield of maize，the growth formation and yield were simulated by WOFOST under different drought grades at seedling stage，jointing stage and filling stage. The drought grades were indicated by the number of no rain days. The results showed that 5 days to 40 days continuous drought occurred at single development stage affected on the growth and yield formation. The longer the drought period was，the larger yield reduction rate got. The reduction rate was not obvious changed when continuous drought for 5-40 days，occurred at the seedling stage. And the reduction rate was also not obvious changed when continuous drought for 5-10 days，occurred at the jointing stage. The reduction rate was 1.37% and 6.54%，7.55% and 18.78%，16.81% and 32.96%，33.28% and 52.96%，53.22% and 60.17%，61.09% and 66.84% when continuous drought for 15 days，20 days，25 days，30 days，35 days，40 days，occurred at the jointing stage of 2013 and 2016，respectively. The reduction rate was not obvious changed when continuous drought for 5-40 days，occurred at the filling stage of 2013. The yield reduction rate was 1.80%-8.20% when continuous drought for 15-30 days occurred at the filling stage. The results indicated that yield reduction was more significant when drought occurred at the jointing stage than that at the other two stages. The results could provide scientific basis for maize planting and drought disaster mitigation in the northeast of China.

Key words Drought;Crop model;Yield;Impact simulation

玉米是我國(guó)三大糧食作物和重要飼料作物之一，玉米播種面積約占我國(guó)耕地面積的20%，產(chǎn)量約占全國(guó)糧食總產(chǎn)量的25%[1]。1951—2000年?yáng)|北地區(qū)南部向干旱發(fā)展，這顯然與全球氣候變暖的大背景有關(guān)[2]。干旱對(duì)玉米的生長(zhǎng)發(fā)育產(chǎn)生影響從而最終導(dǎo)致產(chǎn)量下降，一般可使玉米減產(chǎn)20%～30%，是影響玉米生產(chǎn)的重要因素[3]。玉米產(chǎn)量下降的程度不但取決于干旱的嚴(yán)重程度，還取決于干旱發(fā)生時(shí)玉米的生長(zhǎng)階段。干旱發(fā)生的程度、持續(xù)時(shí)間及生育進(jìn)程的不同，對(duì)玉米產(chǎn)量的影響有所不同[4]。不同生育時(shí)期干旱的研究結(jié)果表明，任何生育時(shí)期的土壤干旱均會(huì)導(dǎo)致玉米減產(chǎn)，其中拔節(jié)期干旱減產(chǎn)較重，苗期相對(duì)較輕[5]。初夏旱對(duì)春玉米的產(chǎn)量影響較大，春旱和伏旱對(duì)春玉米產(chǎn)量、初夏旱和伏旱對(duì)夏玉米產(chǎn)量影響均較小[6]。目前，玉米在不同發(fā)育階段受到不同程度干旱影響下產(chǎn)量的定量響應(yīng)特征研究較少，運(yùn)用作物模型進(jìn)行玉米干旱模擬研究能夠較好地彌補(bǔ)田間水分控制試驗(yàn)的不足，同時(shí)能夠更加細(xì)致地模擬不同發(fā)育階段受到干旱影響的玉米產(chǎn)量的定量變化。因此，筆者擬基于WOFOST作物生長(zhǎng)模型，模擬分析不同持續(xù)日數(shù)干旱對(duì)東北地區(qū)玉米產(chǎn)量的影響，為區(qū)域防災(zāi)減災(zāi)提供科學(xué)依據(jù)。

1 資料與方法

1.1 研究區(qū)概況

試驗(yàn)地為遼寧省錦州農(nóng)業(yè)氣象試驗(yàn)站。錦州所在區(qū)域?qū)俚湫蜏貛Ъ撅L(fēng)型氣候，年平均氣溫9.5℃，1月平均氣溫-8.0℃，7月平均氣溫24.4℃，年均降水量565.9 mm，主要作物玉米的生育期為5—9月，土壤為典型棕壤。土壤凋萎濕度為0.105 cm3/cm3，田間持水量為0.361 cm3/cm3，飽和含水量為0.519 cm3/cm3。

1.2 干旱模擬試驗(yàn)研究

供試玉米品種為丹玉39。干旱模擬時(shí)段選擇在玉米的苗期、拔節(jié)期和灌漿期，模擬時(shí)間從苗期、拔節(jié)期和灌漿期的第5 d開(kāi)始，設(shè)置玉米生長(zhǎng)季內(nèi)單一發(fā)育期發(fā)生干旱3種處理情景。單一發(fā)育期發(fā)生干旱指在正常年景下苗期、拔節(jié)期和灌漿期分別發(fā)生持續(xù)時(shí)間為5、10、15、20、25、30、35和40 d的干旱，苗期分別記為M05、M10、M15、M20、M25、M30、M35、M40，拔節(jié)期分別記為B05、B10、B15、B20、B25、B30、B35、B40，灌漿期分別記為G05、G10、G15、G20、G25、G30、G35、G40。該研究將干旱定義為在正常降水條件下，從控水開(kāi)始，將氣象要素中的降水量均設(shè)定為0，其余氣象要素不變。

1.3 減產(chǎn)率的計(jì)算

以2013、2016年為模擬年份，同時(shí)為了使其更接近正常年份，模擬時(shí)將2013、2016年氣象要素中的逐日降水量以多年平均值（1981—2010年）替代，其他要素均采用2013、2016年的逐日實(shí)測(cè)值，并將其記為CK，相應(yīng)地其模擬產(chǎn)量值記為Yck（kg/hm2）。干旱條件下模擬的產(chǎn)量與正常年份的模擬產(chǎn)量作對(duì)比得到減產(chǎn)率，即

Di=×100（1）

式中，Di為第i年不同程度干旱下（干旱持續(xù)日數(shù)為5、10、15、20、25、30、35、40 d等）的減產(chǎn)率（%），若Di為正數(shù)，則為增產(chǎn)率;Ydi為第i年不同程度干旱下模擬的產(chǎn)量（kg/hm2）。

1.4 作物模型及其適宜性檢驗(yàn)

1.4.1 模型簡(jiǎn)介 WOFOST模型是荷蘭開(kāi)發(fā)的一個(gè)著名的機(jī)理性模型，描述作物的基本生理過(guò)程，包括光合作用、呼吸作用等，并描述這些過(guò)程如何受環(huán)境影響，模擬步長(zhǎng)為日。模型采用“積溫法”，以溫度為驅(qū)動(dòng)，模擬作物發(fā)育過(guò)程。

1.4.2 模型參數(shù)的調(diào)試與確定 準(zhǔn)確的參數(shù)是模擬結(jié)果準(zhǔn)確和穩(wěn)定的前提，需根據(jù)歷史資料，進(jìn)行模型參數(shù)的本地化。作物模型參數(shù)主要包括作物參數(shù)和土壤參數(shù)。作物參數(shù)又包括發(fā)育參數(shù)和生長(zhǎng)參數(shù)。發(fā)育參數(shù)主要包括不同發(fā)育階段所需的有效積溫和光周期影響因子等，這些參數(shù)是作物品種的固有屬性，該研究主要根據(jù)前人研究成果和實(shí)際發(fā)育期計(jì)算得到。生長(zhǎng)參數(shù)主要包括光合速率、呼吸速率、比葉面積、光合產(chǎn)物轉(zhuǎn)化系數(shù)、干物質(zhì)分配系數(shù)及葉片衰老系數(shù)等，這些參數(shù)的確定主要來(lái)自于田間試驗(yàn)數(shù)據(jù)，并采用“試錯(cuò)法”進(jìn)行訂正和確定。土壤參數(shù)主要包括與土壤本身特性相關(guān)的物理參數(shù)和初始條件，如凋萎濕度、田間持水量、飽和含水量、飽和導(dǎo)水率、下滲速率以及初始土壤水分含量等，這些參數(shù)主要根據(jù)實(shí)際田間土壤成分測(cè)定結(jié)果而獲得，部分參數(shù)依據(jù)前人研究成果而定。

為確保參數(shù)本地化后的模型不受典型氣象災(zāi)害年、氣候變化等的影響，以及以后區(qū)域化應(yīng)用，選擇2010、2011年的氣象數(shù)據(jù)和農(nóng)業(yè)氣象觀測(cè)數(shù)據(jù)進(jìn)行參數(shù)本地化。結(jié)果顯示，玉米開(kāi)花天數(shù)模擬值與實(shí)測(cè)值相差2 d，玉米成熟天數(shù)模擬值與實(shí)測(cè)值相差2 d，產(chǎn)量模擬值與實(shí)測(cè)值均方根誤差為166.88 kg/hm2。由此可見(jiàn)，模型對(duì)玉米生長(zhǎng)發(fā)育及產(chǎn)量模擬效果較好，參數(shù)本地化具有較好的初步效果。

1.4.3 模型適宜性檢驗(yàn) 通過(guò)模型模擬結(jié)果與實(shí)測(cè)結(jié)果的圖形比較以及各項(xiàng)評(píng)價(jià)指標(biāo)來(lái)評(píng)價(jià)模型在錦州的適宜性。采用以下統(tǒng)計(jì)量作為檢驗(yàn)?zāi)Ｐ偷闹笜?biāo)：模擬值與實(shí)測(cè)值間的決定系數(shù)（R2）、歸一化均方根誤差（NRMSE）。其中，R2可反映模擬值與實(shí)測(cè)值間的一致性，其值愈接近1說(shuō)明模擬效果愈好;NRMSE值反映了模擬值與實(shí)測(cè)值間的絕對(duì)誤差和相對(duì)誤差，其值愈小，表明誤差值愈小[7]。以2012—2017年連續(xù)6年的發(fā)育期和同期產(chǎn)量對(duì)確定的模型參數(shù)進(jìn)行適宜性檢驗(yàn)，結(jié)果見(jiàn)表1和圖1。

2 結(jié)果與分析

2.1 苗期干旱對(duì)玉米產(chǎn)量的影響模擬

由圖2可知，當(dāng)玉米苗期持續(xù)干旱5～40 d時(shí)，最終產(chǎn)量比CK略增?？梢?jiàn)，如果苗期發(fā)生干旱，對(duì)玉米最終產(chǎn)量沒(méi)有產(chǎn)生負(fù)面影響，這與紀(jì)瑞鵬等[1]的研究一致。

2.2 拔節(jié)期干旱對(duì)玉米產(chǎn)量的影響模擬

由圖3可知，當(dāng)玉米拔節(jié)期持續(xù)干旱5～10 d時(shí)，2013年和2016年的最終產(chǎn)量均無(wú)明顯變化;持續(xù)干旱15、20、25、30、35、40 d時(shí)，2013年分別減產(chǎn)1.37%、7.55%、16.81%、33.28%、53.22%和61.09%，2016年分別減產(chǎn)6.54%、18.78%、32.96%、52.96%、60.17%和66.84%?？梢?jiàn)，拔節(jié)期發(fā)生15～40 d的持續(xù)干旱，會(huì)對(duì)玉米籽粒灌漿產(chǎn)生負(fù)面影響，并最終導(dǎo)致產(chǎn)量降低，且干旱持續(xù)時(shí)間越長(zhǎng)，對(duì)灌漿的影響越早，減產(chǎn)越嚴(yán)重。

2.3 灌漿期干旱對(duì)玉米產(chǎn)量的影響模擬

由圖4可知，當(dāng)玉米灌漿期持續(xù)干旱5～20 d時(shí)，2013年的最終產(chǎn)量無(wú)明顯變化;持續(xù)干旱25、30 d時(shí)，2013年分別減產(chǎn)0.09%和0.12%，基本無(wú)影響;當(dāng)玉米灌漿期持續(xù)干旱5～10 d時(shí)，2016年的最終產(chǎn)量無(wú)明顯變化;持續(xù)干旱15、20、25和30 d時(shí)，2016年分別減產(chǎn)1.75%、3.87%、6.06%和8.10%?？梢?jiàn)，灌漿期發(fā)生干旱對(duì)玉米的影響與前期生長(zhǎng)狀況等多方面因素有關(guān)?？傮w看，灌漿期發(fā)生15～30 d的持續(xù)干旱，會(huì)對(duì)玉米籽粒灌漿產(chǎn)生負(fù)面影響，并最終導(dǎo)致產(chǎn)量降低，且干旱持續(xù)時(shí)間越長(zhǎng)，對(duì)灌漿的影響越早，減產(chǎn)越嚴(yán)重。與拔節(jié)期相比，灌漿期發(fā)生相同持續(xù)時(shí)間的干旱，其受影響程度和減產(chǎn)幅度均小于拔節(jié)期。

3 結(jié)論與討論

玉米苗期、拔節(jié)期和灌漿期分別發(fā)生持續(xù)日數(shù)相同的干旱時(shí)，拔節(jié)期干旱對(duì)玉米產(chǎn)量的影響程度最大，苗期影響最小。表明玉米苗期較為耐旱，拔節(jié)以后對(duì)水分虧缺的敏感性增強(qiáng)。具體表現(xiàn)為當(dāng)玉米苗期持續(xù)干旱5～40 d時(shí)，2013、2016年模擬最終產(chǎn)量略增。當(dāng)玉米拔節(jié)期持續(xù)干旱5～10 d時(shí)，2013年和2016年玉米的最終產(chǎn)量均無(wú)明顯變化;持續(xù)干旱15～20 d時(shí)，2013年減產(chǎn)10%以內(nèi)，2016年減產(chǎn)20%以內(nèi)。持續(xù)干旱25～40 d時(shí)，2013年減產(chǎn)15%～62%，2016年減產(chǎn)30%～67%。當(dāng)玉米灌漿期持續(xù)干旱5～10 d時(shí)，2013年和2016年玉米的最終產(chǎn)量基本無(wú)明顯影響;當(dāng)玉米灌漿期持續(xù)干旱15～30 d時(shí)，2013年玉米產(chǎn)量無(wú)明顯影響，2016年玉米減產(chǎn)1.8%～8.2%。

干旱對(duì)作物的危害程度與其發(fā)生的生育階段、持續(xù)時(shí)間長(zhǎng)短以及作物的自身品種特性如品種類型等有關(guān)。玉米生育期需水量與品種、氣候等多方面因素有關(guān)，在不同生長(zhǎng)發(fā)育期，其需水量也不盡相同。苗期水分脅迫結(jié)束后，植株可快速地部分或全部彌補(bǔ)前期干旱所減少的生長(zhǎng)量，而且苗期適當(dāng)干旱可促進(jìn)玉米根系的發(fā)育，從而增強(qiáng)抗旱能力[5]。拔節(jié)期后，需水逐漸增多，抽雄吐絲期是玉米的水分臨界期，如果這時(shí)供水不足或不及時(shí)，對(duì)產(chǎn)量影響很大[8]。該研究對(duì)苗期、拔節(jié)期、灌漿期不同階段干旱進(jìn)行模擬試驗(yàn)，從模擬結(jié)果上看，當(dāng)玉米在苗期、拔節(jié)期和灌漿期分別發(fā)生干旱，且發(fā)生持續(xù)日數(shù)相同的干旱時(shí)，拔節(jié)期干旱對(duì)產(chǎn)量形成的影響程度最大。即苗期較為耐旱，拔節(jié)以后玉米對(duì)水分虧缺越來(lái)越敏感，灌漿期影響較小。這一研究結(jié)論與白向歷等[5，9]的研究結(jié)果相一致，符合玉米全生育期需水規(guī)律。因此，玉米拔節(jié)至抽孕穗期要及時(shí)補(bǔ)充水分，以滿足其正常生長(zhǎng)發(fā)育需要。

該研究采用數(shù)值模擬的方法，研究結(jié)論是針對(duì)單站的結(jié)果，且干旱的影響效應(yīng)還與選擇的玉米品種熟型有關(guān)，相關(guān)研究有待進(jìn)一步深入。

參考文獻(xiàn)

[1] 紀(jì)瑞鵬，車宇勝，朱永寧，等. 干旱對(duì)東北春玉米生長(zhǎng)發(fā)育和產(chǎn)量的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2012（11）：3021-3026.

[2] 謝安，孫永罡，白人海. 中國(guó)東北近50年干旱發(fā)展及對(duì)全球氣候變暖的響應(yīng)[J]. 地理學(xué)報(bào)，2003（S1）：75-82.

[3] 齊偉，張吉旺，王空軍，等. 干旱脅迫對(duì)不同耐旱性玉米雜交種產(chǎn)量和根系生理特性的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2010，21（1）：48-52.

[4] 白莉萍，隋方功，孫朝暉，等. 土壤水分脅迫對(duì)玉米形態(tài)發(fā)育及產(chǎn)量的影響[J]. 生態(tài)學(xué)報(bào)，2004（7）：1556-1560.

[5] 白向歷，孫世賢，楊國(guó)航，等. 不同生育時(shí)期水分脅迫對(duì)玉米產(chǎn)量及生長(zhǎng)發(fā)育的影響[J]. 玉米科學(xué)，2009（2）：60-63.

[6] 欒慶祖，葉彩華，莫志鴻，等. 基于WOFOST模型的玉米干旱損失評(píng)估：以北京為例[J]. 中國(guó)農(nóng)業(yè)氣象，2014，35（3）：311-316.

[7] 劉志娟，楊曉光，王靜，等. APSIM玉米模型在東北地區(qū)的適應(yīng)性[J]. 作物學(xué)報(bào)，2012，38（4）：740-746.

[8] 張建平，何永坤，王靖，等. 不同發(fā)育期干旱對(duì)玉米籽粒形成與產(chǎn)量的影響模擬[J]. 中國(guó)農(nóng)業(yè)氣象，2015，36（1）：43-49.

[9] 張淑杰，張玉書，紀(jì)瑞鵬，等. 東北地區(qū)玉米干旱時(shí)空特征分析[J]. 干旱地區(qū)農(nóng)業(yè)研究，2011（1）：231-236.

責(zé)任編輯：鄭丹丹