李 昊，譚方穎，王建林**，譚凱炎，徐 英，王志偉

(1.哈爾濱師范大學(xué)，哈爾濱 150025；2.國家氣象中心，北京 100081；3.中國氣象科學(xué)研究院，北京 100081；4.山東省夏津縣氣象局，夏津 253200；5.山西省氣候中心，太原 030006)

?

華北地區(qū)夏玉米干物質(zhì)分配系數(shù)的模擬***

李 昊1，譚方穎2，王建林2**，譚凱炎3，徐 英4，王志偉5

(1.哈爾濱師范大學(xué)，哈爾濱 150025；2.國家氣象中心，北京 100081；3.中國氣象科學(xué)研究院，北京 100081；4.山東省夏津縣氣象局，夏津 253200；5.山西省氣候中心，太原 030006)

摘要：干物質(zhì)分配系數(shù)是驅(qū)動玉米生長發(fā)育模型的關(guān)鍵參數(shù)。利用2013年、2014年連續(xù)兩年在山東夏津、河北固城、山西運城進行的田間試驗觀測資料，采用比值法、線性回歸等訂正方法，獲取完整的玉米全生育期內(nèi)生物量序列；在此基礎(chǔ)上，以有效積溫模擬的發(fā)育進程為自變量，構(gòu)建了華北夏玉米干物質(zhì)分配的動態(tài)變化模型。結(jié)果表明：（1）三站玉米干物質(zhì)分配系數(shù)有相同的時間動態(tài)變化特征。葉片干物質(zhì)分配系數(shù)從出苗開始持續(xù)減少；莖稈干物質(zhì)分配系數(shù)先增后減，最大值出現(xiàn)在抽雄前后；穗棒干物質(zhì)分配系數(shù)在玉米抽雄后持續(xù)增加，抽雄后20d左右達到1，即干物質(zhì)不再向葉莖分配。（2）華北夏玉米生育期內(nèi)葉、莖的干物質(zhì)分配系數(shù)均可用分段式非線性模型模擬。葉的干物質(zhì)分配系數(shù)以抽雄后10～15d為界，之前干物質(zhì)分配系數(shù)隨發(fā)育進程可用三次多項式動態(tài)模型模擬，之后變?yōu)?；莖的干物質(zhì)分配系數(shù)以抽雄后20～25d為界，之前干物質(zhì)分配系數(shù)隨發(fā)育進程可用四次多項式動態(tài)模型模擬，之后變?yōu)?；穗棒干物質(zhì)分配系數(shù)通過依據(jù)任意發(fā)育階段葉、莖、穗棒的干物質(zhì)分配系數(shù)之和為1的原則計算求得。檢驗結(jié)果表明華北地區(qū)夏玉米干物質(zhì)分配系數(shù)動態(tài)模型模擬效果良好。

關(guān)鍵詞：干物質(zhì)分配系數(shù)；模型；田間試驗；積溫；發(fā)育進程

李昊,譚方穎,王建林,等.華北地區(qū)夏玉米干物質(zhì)分配系數(shù)的模擬[J].中國農(nóng)業(yè)氣象,2016,37(3):335-342

玉米是自然界貯存能量最多的植物之一，是重要的糧食、飼料和經(jīng)濟兼用作物，其產(chǎn)量的高低直接影響國家糧食安全和玉米產(chǎn)業(yè)的發(fā)展[1]。華北是中國玉米的主要產(chǎn)區(qū)之一，播種面積和總產(chǎn)量均占全國的30%[2]，在華北地區(qū)開展玉米田間試驗研究對保障玉米產(chǎn)業(yè)發(fā)展具有重要意義。

干物質(zhì)生產(chǎn)是形成作物經(jīng)濟產(chǎn)量的基礎(chǔ)。玉米與其它綠色植物一樣，都是通過光合作用產(chǎn)生干物質(zhì)，并按照一定的生物學(xué)規(guī)律將干物質(zhì)分配至各器官[3]。而干物質(zhì)分配系數(shù)是隨發(fā)育進程變化的一個動態(tài)參數(shù)，能精確表達干物質(zhì)向莖、葉和儲存器官的分配情況，被廣泛應(yīng)用于作物生長發(fā)育模擬模型中。如目前被國內(nèi)廣泛運用的荷蘭WOFOST、ORYZA等作物模型[4-7]。由于干物質(zhì)分配系數(shù)表征玉米干物質(zhì)積累與分配規(guī)律[8-12]，因此確定干物質(zhì)分配系數(shù)要基于大量的干物質(zhì)積累與分配試驗。國外對此進行了較多的相關(guān)研究。但受品種、管理措施、環(huán)境等因素影響，國外模型中的玉米干物質(zhì)分配系數(shù)并不適于國內(nèi)玉米生長特性[13-14]，直接影響了其在國內(nèi)的模擬效果，也成為在國內(nèi)推廣應(yīng)用的瓶頸。因此，基于國內(nèi)玉米種植試驗，構(gòu)建適宜中國區(qū)域玉米生長發(fā)育的干物質(zhì)分配系數(shù)，有利于提高國外相關(guān)作物生長發(fā)育模型在國內(nèi)的應(yīng)用效果。本文以華北地區(qū)夏玉米田間試驗為基礎(chǔ)，通過對試驗數(shù)據(jù)的分析訂正，構(gòu)建華北夏玉米干物質(zhì)分配系數(shù)模型，以供相關(guān)研究參考使用。

1 資料與方法

1.1 試驗介紹

（1）試驗點概況

在華北夏玉米種植區(qū)，根據(jù)空間代表性并參考夏玉米種植面積，選擇山東夏津、河北固城、山西運城試驗站，進行大田玉米觀測試驗。其中夏津試驗地為山東省德州市夏津縣氣象局下屬氣象觀測站，固城試驗地為中國氣象科學(xué)研究院農(nóng)業(yè)氣象實驗基地，運城試驗地為山西省運城市氣象局下屬氣象觀測站。

（2）試驗時間

分別在2013年、2014年進行夏玉米田間試驗，與當?shù)卮筇镒魑锿诓シN，同期收獲。詳情見表1。

表1 各試驗點玉米播種和收獲日期（月-日）Table 1 Sawing and harvest date of maize in the three experimental stations（mm-dd）

（3）試驗設(shè)計

每年在山東夏津、河北固城、山西運城3個試驗點分別設(shè)立2個試驗小區(qū)，單個小區(qū)面積約30m2，共計6個小區(qū)。每個試驗小區(qū)設(shè)3個重復(fù)處理，均選用鄭單958為試驗品種，在播種前施底肥，拔節(jié)-孕穗前追肥1次，保證土壤養(yǎng)分供應(yīng)。采取遮雨和灌溉的措施，全生育期土壤相對濕度保持在60%～100%，無旱澇災(zāi)害發(fā)生。采取人工播種方式，按當?shù)仄毡椴シN時間進行播種，行距50cm，株距35cm，深度約為7cm，密度約為6株·m-2。

（4）觀測項目

生物量觀測：在七葉、拔節(jié)、拔節(jié)后10d、抽雄后10d、抽雄后20d以及收獲（成熟）期，按玉米長勢好、中、差在每個小區(qū)分別隨機選擇1個整株玉米，連根到地上部分整株取下，放入袋中，帶回試驗室。將葉、莖、穗等器官分離，放入烘干機烘干，按器官分別稱重，記錄干物質(zhì)質(zhì)量。各器官質(zhì)量取3株的平均值。2013年觀測中發(fā)現(xiàn)部分觀測時間間隔較大，故2014年各站在玉米拔節(jié)前10d、拔節(jié)后20d加測兩次生物量數(shù)據(jù)，

生育期觀測：根據(jù)農(nóng)業(yè)氣象觀測規(guī)范，判斷并記錄玉米播種期、出苗期、三葉期、七葉期、拔節(jié)期、抽雄期、開花期、吐絲期、乳熟期、收獲（成熟）期等發(fā)育期資料。

土壤濕度觀測：采用烘干稱重法。根據(jù)玉米發(fā)育特征，設(shè)定在拔節(jié)前和拔節(jié)后取樣深度不同，拔節(jié)前為0～20cm，拔節(jié)后為0～50cm。土壤濕度觀測時間從播種前一天開始直至收獲，每7日測定1次，取樣方式為人工鉆土法。

氣象資料觀測：逐日最高溫度和最低溫度數(shù)據(jù)均來自各試驗站的自動氣象觀測站。

1.2 生物量數(shù)據(jù)訂正

田間試驗受天氣因素、人為操作及取樣的隨機性等因素的影響，導(dǎo)致生物量觀測資料出現(xiàn)缺測、個別數(shù)據(jù)異常等現(xiàn)象，需要對觀測數(shù)據(jù)進行訂正處理。常用的訂正方法有比值法、差值法、均值法、線性回歸法等[15-21]。本文根據(jù)不同情況，采用如下方法進行資料訂正。

（1）比值法。比值法是通過尋找數(shù)據(jù)間的特定比值關(guān)系，對問題數(shù)據(jù)進行訂正的方法。如2013年山西運城兩個小區(qū)七葉期觀測得到單株葉干重分別為1.70g、1.90g，莖干重數(shù)據(jù)為0，分析判斷認為，葉、莖可能未分開稱重?？紤]到相同條件下，同一發(fā)育階段各站玉米生長情況大致相同，參考計算出的同期山東、河北葉、莖干重之比（分別為1.55和1.67）的平均值1.61，用比值法訂正得到兩個小區(qū)七葉期單株葉干重分別為1.04g、1.17g，莖干重分別為0.66g、0.73g。研究期運用比值法共處理相似問題數(shù)據(jù)3個。

（2）線性回歸法。在觀測數(shù)據(jù)中，個別數(shù)據(jù)出現(xiàn)異常，明顯不符合作物生長規(guī)律，此類情況下，采用線性回歸模型進行訂正。如2014年河北固城一個小區(qū)拔節(jié)后20d觀測得到單株葉干重數(shù)據(jù)為9.65g，低于前一次觀測（拔節(jié)后10d）的9.75g，判斷為異常觀測值。本文根據(jù)2014年河北固城另一重復(fù)小區(qū)的觀測數(shù)據(jù)建立單株葉干重與發(fā)育進程的回歸方程，即

式中，y為單株葉干重（g），x為發(fā)育進程。方程通過了0.05水平的顯著性檢驗（n=7，R0.05=0.754）。根據(jù)發(fā)育進程訂正拔節(jié)后20d觀測值。研究期用線性回歸法處理相似問題數(shù)據(jù)共7次。

1.3 夏玉米發(fā)育進程計算

干物質(zhì)分配系數(shù)是隨發(fā)育進程變化的函數(shù)，而國內(nèi)外作物模型中普遍使用有效積溫模擬作物發(fā)育進程[22]。本文采用由荷蘭瓦赫寧根大學(xué)開發(fā)的WOFOST作物模型[23]中的發(fā)育進程計算公式，利用氣象觀測資料將生物量觀測時對應(yīng)的生育期量化為發(fā)育進程，以便建立干物質(zhì)分配系數(shù)模擬模型。即

式中，Tei為玉米生育期內(nèi)第i日的有效溫度（℃）；iT為該日的平均氣溫（℃），取日最高、最低氣溫的平均值；Tb為玉米發(fā)育的下限溫度（℃），本文取值為6；Tmax.e為玉米發(fā)育的上限溫度（℃），本文取值為30；j是發(fā)育進程，1表示出苗-抽雄期，2表示抽雄-收獲期；TSUM為玉米發(fā)育到某一時刻所需的有效積溫（℃·d）；fred為光周期影響因子，本文取值為1.00；DVS為作物發(fā)育進程，其中出苗期DVS為0，抽雄期DVS為1.00，成熟期DVS 為2.00。

根據(jù)以上公式，按照作物處于不同生長發(fā)育時期，依據(jù)所對應(yīng)的有效積溫，可以將作物生長發(fā)育的時期進行量化。如山東夏津2013年拔節(jié)后10d對應(yīng)的發(fā)育進程為0.80。3站各發(fā)育階段有效積溫詳情見表2。

表2 三站各發(fā)育階段有效積溫統(tǒng)計（℃·d）Table 2 The effective accumulated temperature of different developmental stages in three stations（℃·d）

1.4 干物質(zhì)分配系數(shù)的計算

基于發(fā)育進程（DVS），葉、莖和儲藏器官的分配系數(shù)[21]分別由以下公式計算

式中，F(xiàn)LTB為葉片干物質(zhì)分配系數(shù)，F(xiàn)STB為莖的干物質(zhì)分配系數(shù)，F(xiàn)OTB為儲藏器官（穗棒）的干物質(zhì)分配系數(shù)。WL、WS、WO分別為葉、莖和儲藏器官（穗棒）的干重（g），TAGP為同期植株地上部分總干重（g）。i表示生物量觀測次序，i=1為第一次觀測，葉、莖生物量第一次觀測為七葉期，穗棒生物量第一次觀測為抽雄后10d；i＞1為其后的觀測，葉、莖生物量的觀測時間依次為拔節(jié)、拔節(jié)后10d、抽雄后10d，穗棒生物量的觀測時間依次為抽雄后20d、成熟期。2014年在拔節(jié)前10d和拔節(jié)后20d有兩次加測。同一發(fā)育進程時FLTB、FSTB、FOTB之和為1。

2 結(jié)果與分析

2.1 各試驗點玉米干物質(zhì)分配系數(shù)計算結(jié)果

由表3-表5可見，兩年的試驗均表明，3個試驗站玉米生育進程中葉片的干物質(zhì)分配系數(shù)最大值出現(xiàn)在七葉期-拔節(jié)期，拔節(jié)以后隨發(fā)育進程逐步減?。磺o稈的干物質(zhì)分配系數(shù)最大值出現(xiàn)在拔節(jié)后10～20d，隨后逐漸減小；抽雄后穗棒的分配系數(shù)增加迅速，抽雄后20d左右變?yōu)?，即干物質(zhì)不再向葉莖分配。

表3 2013和2014年山東夏津站實測玉米干物質(zhì)分配系數(shù)Table 3 Partition coefficient of dry matter during maize growing stage in Xiajin station,Shandong province, in 2013 and 2014

表4 河北固城站實測玉米干物質(zhì)分配系數(shù)Table 4 Partition coefficient of dry matter in Gucheng station,Hebei province

表5 山西運城站實測玉米干物質(zhì)分配系數(shù)Table 5 Partition coefficient of dry matter in Yucheng station,Shanxi province

2.2 玉米干物質(zhì)分配系數(shù)模擬模型的建立

將3個試驗點兩年玉米干物質(zhì)分配系數(shù)系列間進行相關(guān)分析，結(jié)果見表6。由表中可以看出，它們相互間的相關(guān)系數(shù)均在0.9以上，均通過了0.05水平的顯著性檢驗（n=8，R0.05=0.706）。2013年與2014年各試驗站的葉、莖及穗棒干物質(zhì)分配系數(shù)相關(guān)系數(shù)也在0.9以上（P＜0.05），說明在水分充分保證的前提下，華北地區(qū)夏玉米干物質(zhì)平均分配系數(shù)具有較穩(wěn)定的年際變化和一致的時間動態(tài)變化特征，可以基于三地的試驗結(jié)果，構(gòu)建適于華北地區(qū)夏玉米主要器官干物質(zhì)分配系數(shù)動態(tài)模型。

進一步基于2013年、2014年山東夏津、河北固城、山西運城試驗站試驗資料，以發(fā)育進程（x）為自變量，以干物質(zhì)分配系數(shù)為因變量，分別建立葉、莖分配系數(shù)隨發(fā)育進程變化的關(guān)系模型，結(jié)果見圖1。由圖1a可見，當生育進程在0～1.46范圍時，葉片干物質(zhì)分配系數(shù)隨生育進程以較快速度下降，二者關(guān)系可用一元三次方程表示，方程復(fù)相關(guān)系數(shù)為0.981，通過了0.05水平的顯著性檢驗（n=37，R0.05=0.320）；當生育進程達到1.46（抽雄后10～15d）后，葉片干物質(zhì)分配系數(shù)減小為0。由圖1b可見，當生育進程在0～1.58范圍時，莖稈干物質(zhì)分配系數(shù)隨生育進程先上升后下降，二者關(guān)系可用一元四次方程表示，方程相關(guān)系數(shù)0.900，通過了0.05水平的顯著性檢驗（n=37，R0.05=0.320）；當生育進程達到1.46（抽雄后20～25d）后，莖稈干物質(zhì)分配系數(shù)減小為0。穗棒的分配系數(shù)由于抽雄后才開始出現(xiàn)，觀測次數(shù)較少，難以準確模擬，所以依據(jù)任意發(fā)育階段葉、莖、穗棒的干物質(zhì)分配系數(shù)之和為1的原則（式4-式6），直接計算求得穗棒干物質(zhì)分配系數(shù)。最終確定玉米不同階段各器官干物質(zhì)分配系數(shù)模型為

（1）葉片干物質(zhì)分配系數(shù)模型

（2）莖稈干物質(zhì)分配系數(shù)模型

（3）穗棒干物質(zhì)分配系數(shù)模型

式中，y1為葉片干物質(zhì)分配系數(shù)；y2為莖稈干物質(zhì)分配系數(shù)；y3為穗棒干物質(zhì)分配系數(shù)；x為發(fā)育進程。

表6 各站間不同年份玉米干物質(zhì)分配系數(shù)的相關(guān)性檢驗Table 6 The correlation coefficients of the partition coefficient between the three experimental stations in 2013 and 2014

圖1 玉米葉、莖干物質(zhì)分配系數(shù)的動態(tài)模型Fig. 1 Dynamic model of maize dry matter partition with development stage(DVS)

2.3 玉米干物質(zhì)分配系數(shù)動態(tài)模擬模型的檢驗

以2012年在山東夏津開展的夏玉米（鄭單958）田間試驗資料為基礎(chǔ)，對玉米干物質(zhì)分配系數(shù)動態(tài)模擬模型進行驗證，結(jié)果見圖2。由圖中可見，模擬值與實測值散點均勻分布在1:1對角線的兩側(cè)，葉、莖、穗棒干物質(zhì)分配系數(shù)模擬值與實測值的RMSE值分別為0.028、0.028、0.033，小于實測數(shù)據(jù)平均值的10%～15%，表明葉、莖干物質(zhì)分配系數(shù)動態(tài)模型模擬效果良好。

圖2 玉米干物質(zhì)分配系數(shù)動態(tài)模型檢驗（n=6）Fig. 2 Dynamic model test results of maize dry matter partition coefficient（n=6）

3 討論與結(jié)論

（1）干物質(zhì)分配系數(shù)是驅(qū)動作物生長模型的重要參數(shù)。現(xiàn)有作物生長模型中用以模擬同化物分配的干物質(zhì)分配系數(shù)，多是非連續(xù)的階段性取值，缺乏統(tǒng)一性與準確性。本文通過在華北地區(qū)夏津、固城、運城等3個地點進行的連續(xù)兩年田間試驗，對試驗數(shù)據(jù)分析訂正基礎(chǔ)上，構(gòu)建一套完整的華北夏玉米葉、莖干物質(zhì)分配系數(shù)動態(tài)模型，實現(xiàn)了計算玉米發(fā)育期內(nèi)任意階段干物質(zhì)分配系數(shù)的可能。經(jīng)檢驗干物質(zhì)分配系數(shù)動態(tài)模型模擬效果較好，基本能夠反映華北玉米生長發(fā)育過程中干物質(zhì)分配的變化情況。

（2）在水分條件滿足情況下，葉、莖隨生育進程呈非線性變化，模型為分階段性函數(shù)。葉的干物質(zhì)分配系數(shù)以發(fā)育進程1.42（抽雄后10～15d）為界，之前為三次多項式模型，之后停止分配；莖的干物質(zhì)分配系數(shù)以發(fā)育進程1.58（抽雄后20～25d）為界，之前為四次多項式模型，之后停止分配。國內(nèi)對玉米干物質(zhì)分配系數(shù)的深入研究報道較少，在影響作物干物質(zhì)分配的諸多因素中，水分條件是影響華北夏玉米最關(guān)鍵的因素。張淑杰等從不同方面研究了水分脅迫對玉米生長和干物質(zhì)轉(zhuǎn)移的影響[24-26]。因此，進一步研究水分脅迫條件下，玉米干物質(zhì)分配系數(shù)的動態(tài)變化規(guī)律，是未來深入研究的方向之一。

參考文獻References

[1]董樹亭,張吉旺.建立玉米現(xiàn)代產(chǎn)業(yè)技術(shù)體系,加快玉米生產(chǎn)發(fā)展[J].玉米科學(xué),2008,16(4):18-20.

Dong S T,Zhang J W.The establishment of maize modern industrial technology system,accelerate the development of maize production[J].Journal of Maize Sciences,2008,16(4): 18-20.(in Chinese)

[2]趙先麗,張玉書,紀瑞鵬,等.遼寧春玉米出苗期水分脅迫試驗初探[J].氣象與環(huán)境學(xué)報,2010,26(4):35-39.

Zhao X L,Zhang Y S,Ji R P,et al.Preliminary discussion on water stress of spring maize during seedling in Liaoning province[J].Journal of Meteorology and Environment,2010, 26(4):35-39.(in Chinese)

[3]王建林.現(xiàn)代農(nóng)業(yè)氣象業(yè)務(wù)[M].北京:氣象出版社,2010: 50-51.

Wang J L.Modern agricultural meteorological operations[M]. Beijing:China Meteorological Press,2010:50-51. (in Chinese)

[4]王銳,李亞飛,張麗娟,等.土壤濕度驅(qū)動WOFOST模型及其適應(yīng)性[J].中國農(nóng)業(yè)氣象,2015,36(3):263-271.

Wang R,Li Y F,Zhang L J,et al.WOFOST model based on soil moisture driven and its adaptability[J].Chinese Journal of Agrometeorology,2015,36(3):263-271.(in Chinese)

[5]馬玉平,王石立,張黎.針對華北小麥越冬的WOFOST模型改進[J].中國農(nóng)業(yè)氣象,2005,26(3):145-149.

Ma Y P,Wang S L,Zhang L.Study on improvement of WOFOST against overwinter of wheat in north China[J]. Chinese Journal of Agrometeorology,2005,26(3):145-149.(in Chinese)

[6]韓湘云,景元書,浩宇,等.基于田間試驗的水稻模型ORYZA2000區(qū)域參數(shù)比較[J].干旱氣象,1991,2(3):294.

Han X Y,Jing Y S, Hao Y,et al.The comparison of regional parameters Of ORYZA2000 based on the field test[J].Arid Meteorology,1991,2(3):294.(in Chinese)

[7]薛昌穎,楊曉光,Bouman B A M,等.ORYZA2000模型模擬北京地區(qū)旱稻的適應(yīng)性初探[J].作物學(xué)報,2005,12(31): 1567-1571.

Xue C Y,Yang X G, Bouman B A M,et al.Preliminary approach on adaptability of ORYZA2000 model for aerobic rice in Beijing region[J].Acta Agronomica Sinica,2005,12(31): 1567-1571.(in Chinese)

[8]Grant R F.Simulation of maize phenology[J].Ageron J,1989, 81:451-457.

[9]Penning de Vries F W T,Jansen D M,Ten Berge H F M,et al.Simulation of ecophysiological processes of growth in several annual crops[M].Wageningen:Centre for Agricultural Publishing and Documentation,1989:111-155.

[10]Karlen D L,Lflannery R,Sadler E J.Dry matter nitrogen phosphorus and potassium accumulation rate by corn on Norfolk Loamy Sand[J].Agron,1987,79:649-656.

[11]尹紅征,呂冰清,鄭國清,等.玉米產(chǎn)量形成模擬模型研究[J].華北農(nóng)學(xué)報,2004,19(3):73-76.

Yi H Z,Lv B Q,Zheng G Q,et al.Simulation models of the maize yield formation[J].Acta Agriculturae Boreali-Sinica,2004,19(3):73-76.(in Chinese)

[12]王忠孝,高學(xué)曾,滕世玉.玉米生理[M].北京:農(nóng)業(yè)出版社, 1997:270-300.

Wang Z X,Gao X Z,Teng S Y.Corn physiology[M]. Beijing:China Agricultural Press,1997:270-300.(in Chinese)

[13]王石立,馬玉平.作物生長模擬模型在我國農(nóng)業(yè)氣象業(yè)務(wù)中的應(yīng)用研究進展及思考[J].氣象,2008,34(6): 3-10.

Wang S L,Ma Y P,Teng S Y.The progress in application of crop growth simulation models to agro- meteorological services in China[J]. Meteorological Monthly,2008,34(6): 3-10.(in Chinese)

[14]劉克禮,劉景輝.春玉米干物質(zhì)積累分配與轉(zhuǎn)移規(guī)律的研究[J].內(nèi)蒙古農(nóng)牧學(xué)院學(xué)報,1994,15(1):1-9.

Liu K L,Liu J H.Accumulation of dry matter distribution and transfer laws of spring corn research[J].Journal of Inner Mongolia College of Agriculture,1994,15(1):1-9.(in Chinese)

[15]姚萍.普通物理實驗中的數(shù)據(jù)處理[J].黔西南民族師范高等?？茖W(xué)校學(xué)報,2007,(2):105-108.

Yao P.Data processing in the general physics experiment[J]. Journal of Southwest Guizhou Teachers College for Nationalities,2007(2):105-108.(in Chinese)

[16]幺枕生,丁裕國.氣候統(tǒng)計[M].北京:氣象出版社,1990: 52-58.

Yao Z S,Ding Y G.Climate statistics[M].Beijing:China Meteorological Press,1990:52-58.(in Chinese)

[17]馬力文,李鳳霞,梁旭.寧夏干旱及其對農(nóng)業(yè)生產(chǎn)的影響[J].干旱地區(qū)農(nóng)業(yè)研究,2001,19(4);102-109.

Ma L W,Li F X,Liang X.Characteristics of drought and its nfluence to agriculture in Ningxia[J]. Agricultural Research in the Arid Areas,2001,19(4);102-109.(in Chinese)

[18]李軍,黃敬峰,王秀珍.山區(qū)月平均氣溫的短序列訂正方法研究[J].浙江大學(xué)學(xué)報（農(nóng)業(yè)與生命科學(xué)版）,2005, 31 (2):165-170.

Li J,Hang J F,Wang X Z.Research on adjusting methods of monthly average temperature series in mountainous area[J]. Journal of Zhejiang University(Agriculture & Life Sciences), 2005,31(2):165-170.(in Chinese)

[19]李慶楊,王能超,易大義.數(shù)值分析(第4版) [M].北京:施普林格出版社和清華大學(xué)出版社,2001:70-72.

Li Q Y,Wang N C,Yi D Y.Numerical analysis[M].Beijing: Tsinghua University Press and Springer Verlag, 2001:70-72.(in Chinese)

[20]李莉,朱建躍.T213降水預(yù)報訂正系統(tǒng)的建立與研究[J].應(yīng)用氣象學(xué)報,2006,8:131-133.

Li L,Zhu J Y.The establishment and research of T213 precipitation calibration system[J].Quarterly Journal of Applied Meteorology, 2006,(8):131-133.(in Chinese)

[21]楊松,杞明輝,姚德寬.誤差訂正在預(yù)報集成中的應(yīng)用研究[J].氣象,2003,29(12):22-24.

Yang S,Qi M H,Yao D K.Application of error correction in forecast ensemble[J].Meteorological Monthly, 2003,29(12): 22-24.(in Chinese)

[22]謝云,James R Kiniry.國外作物生長模型發(fā)展綜述[J].作物學(xué)報,2002,(2):190-195.

Xie Y,James R K.Review on crop growth model development abroad[J].Acta Agronomica Sinica,2002, (2): 190-195.(in Chinese)

[23]Hijmans R J,Guiking-Lens I M,van Diepen C A.WOFOST, user guide for the WOFOST 6.0 crop growth simulation model[M].Wageningen:Netherlands Technical document 12, DLOW in and Staring Centre, 1994:145.

[24]張淑杰,張玉書.水分脅迫對玉米生長發(fā)育及產(chǎn)量形成的影響研究[J].中國農(nóng)學(xué)通報,2011,27(12):68-72.

Zhang S J,Zhang Y S.Influences of water stress on growth and development of maize and yield[J]. Chinese Agricultural Science Bulletin,2011,27(12):68-72. (in Chinese)

[25]于志青,于衛(wèi)衛(wèi),譚秀山,等.水分脅迫對夏玉米干物質(zhì)積累與分配的影響[J].華北農(nóng)學(xué)報,2009,24(S2): 149-154.

Yu Z Q,Yu W W,Tan X S.Effect of water stress on dry-matter partition of summer maize[J].Acta Agriculturae Boreali-Sinica,2009,24(S2):149-154. (in Chinese)

[26]張玉書,米娜.土壤水分脅迫對玉米生長發(fā)育的影響研究進展[J].中國農(nóng)學(xué)通報,2012,28(3):1-7.

Zhang Y S,Mi N.Influences of soil water stress on growth and development of maize[J].Chinese Agricultural Science Bulletin,2012,28(3):1-7.(in Chinese)

Simulation on Dry Matter Distribution Coefficient for Summer Maize in North China

LI Hao1, TAN Fang-ying2, WANG Jian-lin2, TAN Kai-yan3, XU Ying4, WANG Zhi-wei5
(1.Harbin Normal University, Harbin 150025, China; 2.National Meteorological Center of China Meteorological Administration, Beijing 100081; 3.Chinese Academy of Meteorological Sciences, Beijing 100081; 4.Xiajin County Meteorological Bureau of Shandong Province, Xiajin, 253200; 5.Climate Center of Shanxi Province, Taiyuan, 030006)

Abstract:Dry matter distribution coefficient was a key parameter for driving maize growth model. Based on the field experimental data of two consecutive years (2013 and 2014) in Xiajin, Shandong province, Gucheng, Hebei province and Yuncheng, Shanxi province, the biomass sequence in the whole growth duration of maize was obtained by using the ratio method and linear regression correction method. On such a basis, the dynamic change model of dry matter distribution for summer maize in North China was established by taking the development stage of effective simulation of accumulated temperature as the independent variable. The results showed that: (1) There were same dynamic change characteristics for dry matter distribution coefficient of maize from three stations. The dry matter distribution coefficient for leaves kept reducing from seedling process, and the dry matter distribution coefficient for stems increased at the beginning and then decreased with the maximum value occurring at tasseling stage. The dry matter partitioning coefficient for ears increased after the tasseling and reached 1 after 20 days, which meant the dry matter was no longer allocated to leaves and stems. (2) Both dry matter distribution coefficients for leaves and stems of summer maize in North China could be simulated by piecewise nonlinear model. The dry matterdistribution coefficient for leaves was distinguished at 10-15 days after tasseling, which could be simulated by three times polynomial dynamic model before tasseling along with the development process, and became 0 after tasseling. The dry matter distribution coefficient for stems was distinguished at 20-25 days after tasseling, which could be simulated by four times polynomial dynamic model before tasseling along with the development process, and became 0 after tasseling. The dry matter distribution coefficient for ears could calculated by the sum of the dry matter distribution coefficient for leaves, stems and ears, which was 1. The results indicated that the simulation effect of dynamic model for dry matter distribution coefficient of summer maize was proficient in North China.

Key words:Dry matter distribution coefficient; Model; Field experiment; Accumulated temperature; Development stage

doi:10.3969/j.issn.1000-6362.2016.03.009

* 收稿日期：2015-12-07**通訊作者。E-mail：wangjl@cma.gov.cn

基金項目：公益性行業(yè)（氣象）科研專項 “干旱對華北玉米生長發(fā)育和產(chǎn)量影響的定量評估技術(shù)研究（GYHY201306038）”

作者簡介：李昊（1990-），碩士生，主要從事農(nóng)業(yè)氣象研究。E-mail：killerlight9023@163.com