張帥杰，段 萌，張寶忠，胡振華

（1.山西農(nóng)業(yè)大學(xué) 林學(xué)院，山西 太谷 030801；2.中國水利水電科學(xué)研究院 流域水循環(huán)模擬與調(diào)控國家重點實驗室，北京 100038）

研究玉米冠層結(jié)構(gòu)對探索其生長發(fā)育規(guī)律、空間分布特性等具有重要意義。玉米冠層結(jié)構(gòu)影響玉米冠層內(nèi)光分布特性及葉片的光合作用速率[1]，進而改變冠層內(nèi)部物質(zhì)和能量循環(huán)，調(diào)節(jié)物質(zhì)分配和作物產(chǎn)量[2-3]。為此，本研究探究了玉米冠層各器官的結(jié)構(gòu)特征，構(gòu)建玉米各器官的形態(tài)結(jié)構(gòu)模型（主要包括節(jié)間、葉片、雄穗、雌穗、葉鞘等器官的模擬模型），為進一步研究玉米冠層內(nèi)生理生態(tài)過程與形態(tài)結(jié)構(gòu)建成提供基礎(chǔ)。

在作物冠層結(jié)構(gòu)模擬方面，國內(nèi)外學(xué)者已經(jīng)開展了大量研究工作。國外學(xué)者對作物器官形態(tài)建成從不同尺度構(gòu)建了作物生長動態(tài)模擬模型。WATANABE等[4]結(jié)合軟件分析水稻的形態(tài)建成，對水稻植株結(jié)構(gòu)進行了三維可視化研究；DORNBUSCH等[5]以春大麥為研究對象，基于器官3D功能-結(jié)構(gòu)模型描述了禾本科作物植株的結(jié)構(gòu)；MCCOWN等[6]分析了日照對小麥葉面積的影響，構(gòu)建了小麥葉面積變化模型。在玉米冠層結(jié)構(gòu)方面，F(xiàn)OURNIER等[7]研究了玉米單個節(jié)間的伸長過程及規(guī)律；GREYSON等[8]用葉位的變化對葉片長寬進行了規(guī)律性的總結(jié)；PARVEZ等[9]探索了葉片長與寬的比例關(guān)系；SALAH等[10]研究了葉鞘的伸長過程。國內(nèi)學(xué)者在作物形態(tài)模擬方面取得了較多的研究進展，在節(jié)間生長模擬方面，孟亞利等[11]、常麗英等[12]討論了分蘗發(fā)生與衰亡的生理生態(tài)過程及其相互關(guān)系、生長度日與小麥節(jié)間的相互關(guān)系，建立了水稻莖蘗生長模型；劉鐵梅等[13]、陳國慶等[14]分別探尋了節(jié)間長度及節(jié)間寬度等形態(tài)指標(biāo)隨生育進程的變化規(guī)律，構(gòu)建了小麥個體和群體莖蘗動態(tài)綜合模型和基于生長度日的節(jié)間生長動態(tài)模擬模型；楊恒山等[15]、宋有洪等[16]分別通過植株生理發(fā)育時間和生物量與玉米節(jié)間的變化規(guī)律，建立了玉米的動態(tài)生長模型；葉片和葉鞘擴展方面，常麗英等[17]、馮陽春等[18]分別分析了水稻葉與生長度日的相互關(guān)系和水稻主莖葉片位置與有效積溫的關(guān)系，建立了水稻葉片生長動態(tài)模型和葉片葉齡模型；嚴(yán)美春等[19]、陳國慶等[20]分別擬合了葉熱間距和生長度日與小麥葉片、葉鞘的相互關(guān)系，構(gòu)建了小麥葉片和葉鞘動態(tài)生長模型；ZHAO等[21]、馬 韞 韜 等[22]、郭 新 宇 等[23]、宋 有 洪 等[24]分 別 通過三維空間、三維數(shù)學(xué)描述以及生物量分析了玉米葉片和葉鞘的形態(tài)變化，從不同的維度建立了玉米的葉片生長模型；穗的發(fā)育方面，劉桃菊等[25]、譚子輝等[26]分別考慮了小麥葉齡、積溫與麥穗的生長發(fā)育的關(guān)系，建立了小麥麥穗生長模型；楊恒山等[15]研究了氮素對玉米生長的影響，建立了基于生理發(fā)育時間的玉米穗生長模型。由于玉米冠層結(jié)構(gòu)較復(fù)雜以及形態(tài)測量數(shù)據(jù)的難獲取性，使得前人研究中多是以單一器官為研究對象，同時以生長度日為基礎(chǔ)對玉米形態(tài)模擬的研究較少，使其器官分配及形態(tài)建成模擬研究成為作物生長與形態(tài)建成模擬研究中的薄弱環(huán)節(jié)。

為了更加系統(tǒng)、詳實地對玉米冠層各器官結(jié)構(gòu)進行表征，本研究結(jié)合玉米各器官發(fā)育規(guī)律，運用Logistic方程構(gòu)建各器官形態(tài)結(jié)構(gòu)模型，系統(tǒng)分析了夏玉米各器官的動態(tài)發(fā)育變化特征，構(gòu)建了基于生理生態(tài)過程的玉米單位器官模擬模型，旨在為在玉米生長模擬過程中實現(xiàn)植株冠層不同器官形態(tài)建成的可視化奠定基礎(chǔ)。

1 材料和方法

1.1 試驗地概況

試驗區(qū)位于北緯39°37′25″、東經(jīng)116°25′51″的中國水利水電科學(xué)研究院國家節(jié)水灌溉北京工程技術(shù)研究中心（北京）大興試驗基地。試驗區(qū)屬于溫帶大陸性季風(fēng)氣候，夏季炎熱多雨，春冬干旱少雨，年平均氣溫12℃左右，年日照時間約2 600 h，多年平均降水量約550 mm，全年水面蒸發(fā)量在1 800 mm以上。試驗基地面積約4 hm2，土質(zhì)為沙壤土，田間持水量為33.4%左右。光照、氣溫及土壤水分等氣象條件適宜夏玉米正常生長。

1.2 試驗材料

供試夏玉米品種為紀(jì)元168，購買于河北廊坊。

1.3 試驗設(shè)計

試驗于2020年6月播種夏玉米，種植密度為6萬株/hm2（株行距30 cm×55 cm），分12個試驗小區(qū)，每個小區(qū)面積56 m2。夏玉米整個生育期內(nèi)施肥2次，基肥施復(fù)合肥（總養(yǎng)分為45%，N∶P2O5∶K2O=15∶15∶15），8月8日追施尿素。每個小區(qū)取代表性玉米3株進行重復(fù)測定，生育期內(nèi)每3 d進行一次測定。生育期內(nèi)降水基本能滿足夏玉米需水，故各處理除播前灌水60 mm以保證出苗外，整個生育期內(nèi)均不再進行灌溉，其他田間管理措施一致。

1.4 測定參數(shù)及方法

將玉米冠層器官分為五大類，分別為節(jié)、葉、葉鞘、雄穗和雌穗。節(jié)間器官主要表征莖稈的形態(tài)結(jié)構(gòu)特點，包括節(jié)間的長度和寬度；葉片主要包括玉米葉長和葉寬動態(tài)模擬；葉鞘主要包括葉鞘長度和寬度的模擬，實現(xiàn)整株玉米的葉鞘器官模擬；雄穗和雌穗分別對其穗長和穗寬進行分析，實現(xiàn)穗器官的模擬。

1.4.1 節(jié)間測定 自拔節(jié)期開始，每3 d定株測量各節(jié)間長度和粗度，測量部位為節(jié)間中間最寬及最窄處直徑，計算二者的平均值。

1.4.2 葉片測定 從第1葉開始，每3 d定株測量葉片長度和寬度，直至葉片定長、定寬。

1.4.3 雄穗測定 自雄穗出現(xiàn)開始，每3 d定株測量雄穗長度和寬度，長度為雄穗和穗稈的總長度，寬度測量部位為最寬處直徑。

1.4.4 雌穗測定 自雌穗出現(xiàn)開始，每3 d定株測量雌穗長度和寬度，寬度使用游標(biāo)卡尺測量。

1.4.5 葉鞘測定 從第1葉開始，每3 d定株測量葉鞘長度和寬度，直至葉鞘定長、定寬。

1.5 數(shù)據(jù)分析

試驗運用Matlab結(jié)合Richard方程進行參數(shù)的率定驗證，運用Excel軟件對數(shù)據(jù)進行分析。

2 結(jié)果與分析

2.1 模型構(gòu)建

2.1.1 玉米生長發(fā)育階段模擬 根據(jù)紀(jì)元168玉米各器官形成日期與生長度日的關(guān)系推算出玉米各器官開始生長的生長度日以及同器官相鄰部位之間節(jié)和葉熱間距。

式中，Temp(i)為一天各時段中每個時段的平均溫度，Tmin為日最低溫度，Tmax為日最高溫度，Tfac(i)為溫度變化因子，Tb為基點溫度，DTT為每天的積溫，GDD為每天積溫的累積形成的生長度日。

2.1.2 玉米葉片動態(tài)模擬 玉米葉片伸長和增寬符合S型曲線生長，因此，采用Logistic方程對葉片進行模擬得到其方程公式（5）～（9）。

式中，YLenn（GDD）為第n片葉片在GDD時刻下 的 葉 片 長 度，YWenn（GDD）為 第n片 葉 片 在GDD時刻下的葉片寬度；YLenn為第n片葉片的葉片長度潛在最大值，YWenn為第n片葉片的葉片寬度 潛 在 最 大 值；t1、s1、m1、t2、s2、m2為 模 型 參 數(shù)；YGDDn為第n片葉片發(fā)育時的初始GDD，GDD為某時刻葉片的生長度日；121.71為本研究玉米品種從播種到出苗這段時間所需的GDD；PHYLL表示相鄰葉片出現(xiàn)的熱時間積溫間隔，a1、b1、c1為葉長潛在最大值方程參數(shù)，a2、b2、c2為葉寬潛在最大值方程參數(shù)。

2.1.3 玉米葉鞘生長動態(tài)模擬 玉米各葉位葉鞘從玉米各節(jié)末端開始向上伸長和增粗直至葉片末端，與葉片的生長具有同時性。采用Logistic方程對葉鞘進行模擬。

式中，YQLenn（GDD）為第n片葉片在GDD時刻葉鞘長度，YQWenn（GDD）為第n片葉片葉鞘寬度；YQLenn為第n片葉片的葉鞘長度潛在最大值，YQWenn為第n片葉片的葉鞘寬度潛在最大值；t3、s3、m3、t4、s4、m4為模 型 參 數(shù)；YQGDDn為 第n片 葉鞘發(fā)育時的初始GDD，GDD為某時刻葉鞘的生長度日；PHYLL表示相鄰葉鞘出現(xiàn)的熱時間積溫間隔。a3、b3、c3為葉鞘長度潛在最大值方程參數(shù)，a4、b4、c4為葉鞘寬度潛在最大值方程參數(shù)。

2.1.4 玉米節(jié)間動態(tài)模擬 玉米各節(jié)間從玉米基部開始向上伸長和增粗，生長具有同時性且符合S型曲線生長。因此，本研究用Logistic方程對玉米節(jié)間進行模擬。

式中，INLenn（GDD）為第n節(jié)節(jié)間在生長度日時刻的長度，INWenn（GDD）為第n節(jié)節(jié)間在生長度日時刻的寬度；INLenn為第n節(jié)節(jié)間長度的潛在最大值，INWenn為第n節(jié)節(jié)間寬度的潛在最大值，t5、s5、m5、t6、s6、m6為模型參數(shù)；INGDDn為第n節(jié)開始伸長和增粗的GDD，GDD為某時刻節(jié)間的生長度日，678.69為本研究玉米品種從播種到出苗這段時間所需的GDD；PHJLL則表示相鄰節(jié)間出現(xiàn)的熱時間積溫間隔。a5、b5、c5為節(jié)間長度潛在最大值方程參數(shù)，a6、b6、c6為節(jié)間寬度潛在最大值方程參數(shù)。

2.1.5 玉米雄穗生長動態(tài)模擬 玉米雄穗在玉米拔節(jié)初期開始生長，根據(jù)雄穗的生長結(jié)合模型得到玉米雄穗伸長、增粗方程。

式中，LT,t為玉米在任一生長發(fā)育時刻的雄穗長度；WT,t為玉米在任一生長發(fā)育時刻的雄穗寬度；LT,max為玉米雄穗長度潛在最大值，WT,max為玉米雄穗寬度潛在最大值，是品種遺傳參數(shù)；t7、s7、m7、t8、s8、m8為模型參數(shù)；XGDDn為雄穗開始生長發(fā)育的生長度日；GDD為雄穗某一時刻的生長發(fā)育的生長度日。

2.1.6 玉米雌穗生長動態(tài)模擬 玉米雌穗在玉米拔節(jié)初期開始出現(xiàn)，在抽穗期開始突然迅速生長，對雌穗進行模擬得到玉米雌穗伸長、增粗方程。

式中，LE,t為玉米在任一生長發(fā)育時刻的雌穗長度，WE,t為玉米在任一生長發(fā)育時刻的雌穗寬度；LE,max為玉米雌穗長度潛在最大值，WE,max為玉米雌穗寬度潛在最大值數(shù)值，是品種遺傳參數(shù)；t9、s9、m9、t10、s10、m10為模型參數(shù)；CGDDn為雌穗 開始生長發(fā)育的生長度日；GDD為雌穗某一時刻的生長發(fā)育的生長度日。

2.2 模型率定與驗證

2.2.1 模型參數(shù)率定 本研究采用2020年實測數(shù)據(jù)進行模型的率定與驗證，并通過縮張算法進行非線性擬合，進一步對玉米各器官模型參數(shù)進行參數(shù)率定。表1為各器官生長模型率定所得到的參數(shù)取值。其中，節(jié)間長度、葉片長度、葉片寬度、葉鞘寬度與葉位存在二次函數(shù)關(guān)系，節(jié)間寬度、葉鞘長度與葉位存在一次函數(shù)關(guān)系，因此，本研究通過二次方程擬合得到各葉位對應(yīng)的潛在最大值，如表2所示。

表1 模型率定參數(shù)取值Tab.1 The model calibration parameter values

表2 潛在最大值方程參數(shù)Tab.2 The equation parameter of potential maximum

2.2.2 模型驗證 利用試驗資料對所建成的玉米節(jié)間、葉片、雄穗、雌穗和葉鞘的長度和寬度模型進行了驗證，結(jié)果表明（圖1～5），本模型在模擬玉米節(jié)間生長在不同GDD時的節(jié)間、葉片、雄穗、雌穗、葉鞘長度和寬度時具有較好的一致性。統(tǒng)計分析表明，模型對葉片長度和寬度預(yù)測的RMSE分別為0.532、0.508；對葉鞘長度和寬度預(yù)測的RMSE分別為0.845、0.198；對節(jié)間長度和寬度預(yù)測的RMSE分別為0.734、0.585；對雄穗長度和寬度預(yù)測的RMSE分別為1.850、0.019；對雌穗長度和寬度預(yù)測的RMSE分別為1.689、0.236。

圖1 玉米葉片在不同GDD時的觀測值與模擬值的比較Fig.1 Comparison of observed and simulated values of maize leaves at different GDD

圖2 玉米葉鞘在不同GDD時的觀測值與模擬值的比較Fig.2 Comparison of observed and simulated values of maize leaf sheaths at different GDD

圖3 玉米節(jié)間在不同GDD時的觀測值與模擬值的比較Fig.3 Comparison of observed and simulated values of maize internodes at different GDD

圖4 玉米雄穗在不同GDD時的觀測值與模擬值的比較Fig.4 Comparison of observed and simulated values of maize tassels at different GDD

圖5 玉米雌穗在不同GDD時的觀測值與模擬值的比較Fig.5 Comparison of observed and simulated values of maize ears at different GDD

3 結(jié)論與討論

時間步長是作物的生理年齡，它與作物的形態(tài)存在聯(lián)系，是能準(zhǔn)確描述發(fā)育狀態(tài)的先決條件。因此，模型要對玉米形態(tài)結(jié)構(gòu)模擬，時間坐標(biāo)的確定就不可或缺。從農(nóng)業(yè)角度來看，時間坐標(biāo)一般通過葉熱間距、生物量、生理生長時間（PDT）以及生長度日（GDD）4種形式進行計算。葉熱間距的計算目前在國內(nèi)運用較少且還未有統(tǒng)一的方法，因此，模型運用時具有不確定性和波動性。宋有洪等[24]建立了節(jié)間、葉片、穗的長寬與生物量的異速生長關(guān)系，能夠相對穩(wěn)定地模擬玉米形態(tài)的建成，該方法主要應(yīng)用于溫室內(nèi)的盆栽作物。楊恒山等[15]通過考慮光周期和溫度等因素以生長發(fā)育日為時間坐標(biāo)建立了模型，并且考慮了氮肥效應(yīng)因子對形態(tài)建成的影響，穩(wěn)定地反映了生理發(fā)育效應(yīng)與作物生長過程的關(guān)系。本研究將試驗觀測與模擬研究相結(jié)合，以玉米形態(tài)特征的變化規(guī)律為基礎(chǔ)，以GDD為時間坐標(biāo)結(jié)合熱間距建立了可系統(tǒng)描述玉米各器官動態(tài)生長的模型，并利用實測數(shù)據(jù)對所建模型進行了初步檢驗，結(jié)果表明，模型能夠較穩(wěn)定地反映玉米的形態(tài)生長過程。

本模型中確定了GDD與各器官形態(tài)參數(shù)的關(guān)系，即葉片開始生長發(fā)育的GDD為121.71，相鄰葉片出現(xiàn)的熱間距間隔為29.42；節(jié)間開始生長發(fā)育的GDD為678.69，相鄰節(jié)間部位出現(xiàn)的熱間距間隔為11.28；雄穗開始生長發(fā)育的GDD為686.69；雌穗開始生長發(fā)育的GDD為798.82。GDD與生長發(fā)育階段對應(yīng)關(guān)系為模型構(gòu)建奠定了基礎(chǔ)。同時通過確定葉片、葉鞘、節(jié)間潛在最大值與其位置的關(guān)系，得出葉片長度、葉片寬度、葉鞘寬度、節(jié)間長度與其位置二次函數(shù)關(guān)系式，葉鞘長度、節(jié)間寬度與其位置一次函數(shù)關(guān)系式，計算不同部位的潛在最大值，繼而計算整個生長發(fā)育生育期內(nèi)的各器官發(fā)育程度。本研究通過模型對于各器官的R2、RMSE以及離散程度進行綜合性分析得出，對于葉片、葉鞘和雄穗發(fā)育擬合程度較好，模擬值與實測值離散程度較小；但節(jié)間和雌穗擬合程度相對較差，模擬值與實測值離散程度相對較大，節(jié)間長度模擬中一部分模擬值略低于實測值，雌穗因發(fā)育的突發(fā)性和迅速性而前期較為不穩(wěn)定，導(dǎo)致整體擬合程度略差。本研究建立了以GDD為時間步長的玉米器官生長模型，對玉米節(jié)間、葉片、雄穗、雌穗以及葉鞘的形態(tài)變化過程進行了定量化模擬，較好地體現(xiàn)了玉米器官建成的動態(tài)特征。用Logistic模型描述了器官的發(fā)育過程，符合玉米器官發(fā)育的生物學(xué)規(guī)律，能更好地量化玉米生長動態(tài)。根據(jù)器官潛在最大值與其位置關(guān)系建立聯(lián)系，既減少了模型的參數(shù)又提升了模型的實用性。但由于玉米形態(tài)數(shù)據(jù)的難獲取性，使模型在不同的生長環(huán)境下適用性仍需做進一步研究。本研究為玉米生長進程預(yù)測和實現(xiàn)生理生態(tài)過程模擬提供一定的技術(shù)支撐，為作物形態(tài)結(jié)構(gòu)研究作物器官冠層可視化的構(gòu)建提供了思路和方法。