曹亦兵，黃收兵，王媛媛，夏雨晴，孟慶鋒，陶洪斌，王璞

?

玉米群體生長與光截獲的動態(tài)模擬及應(yīng)用

曹亦兵，黃收兵，王媛媛，夏雨晴，孟慶鋒，陶洪斌，王璞

（中國農(nóng)業(yè)大學(xué)農(nóng)學(xué)院，北京 100193）

【目的】群體光截獲對玉米生長發(fā)育及產(chǎn)量形成均有重要的影響，合理的群體結(jié)構(gòu)直接決定著玉米產(chǎn)量的增加與穩(wěn)定。群體結(jié)構(gòu)的測量較為繁瑣，費時費力，且不能與當(dāng)?shù)靥栞椛淝闆r充分結(jié)合，從而無法對群體結(jié)構(gòu)和功能做出整體的評價。因此，有必要構(gòu)建合理、簡便的玉米群體光截獲動態(tài)模型，結(jié)合群體不同生育期的形態(tài)指標(biāo)與光輻射情況，對群體結(jié)構(gòu)做出快速而精準(zhǔn)的評價，為不同區(qū)域玉米品種選擇、播種密度確定、以及后期管理等提供參考意見。【方法】將玉米群體結(jié)構(gòu)模型與光分布模型相結(jié)合，建立簡易可靠的群體光截獲動態(tài)模型。該模型可模擬玉米整個生育期群體光截獲的動態(tài)變化，并對不同生育期群體結(jié)構(gòu)做出相應(yīng)的評價。同時，在中國農(nóng)業(yè)大學(xué)吳橋?qū)嶒炚鹃_展了兩個田間試驗，對模型的準(zhǔn)確性和應(yīng)用功能進行了驗證。試驗均為三因素，分別為品種、密度和氮肥處理，其中試驗一：設(shè)置2個品種（鄭單958，丹玉405），3個密度（4.5×104，6.0×104，7.5×104株/hm2）和3個氮肥處理（不施氮；施氮180 kg·hm-2，于播前、大喇叭口期和吐絲期以1﹕4﹕1的比例施入；施氮270 kg·hm-2，于播前、大喇叭口期和吐絲期以4﹕4﹕1的比例施入）；試驗二：設(shè)置1個品種（鄭單958）、1個密度（8.25×104株/hm2）和6個播期（4月20日，5月5日與20日，6月4日與18日，7月3日）。在兩個試驗中，對玉米不同生育期的群體形態(tài)指標(biāo)、干物重、及產(chǎn)量均進行了測定。【結(jié)果】吐絲期和灌漿中期的群體光截獲率模擬值與測量值的相關(guān)系數(shù)分別為0.91和0.85，均達到顯著相關(guān)。運用該模型對兩個田間試驗數(shù)據(jù)的模擬結(jié)果顯示，光截獲率（量）隨灌漿期的推進呈先增后減的趨勢，吐絲后50 d光截獲衰減加速。在試驗地區(qū)，密度超過6×104株/hm2時，光截獲隨密度變化不大；播期為6月中上旬時形成的玉米群體具有較為理想的光截獲能力，結(jié)構(gòu)功能潛力較大。產(chǎn)量與灌漿期有效光輻射呈正相關(guān)，但相關(guān)性隨生育期推進逐漸減弱，表明灌漿后期產(chǎn)量的增加對群體光截獲的需求減弱。在玉米生長發(fā)育后期物質(zhì)轉(zhuǎn)移對產(chǎn)量的貢獻可能發(fā)揮著更加重要的作用。產(chǎn)量與群體光截獲相關(guān)性的動態(tài)變化與品種有關(guān)，緊湊株型優(yōu)于平展型?！窘Y(jié)論】選擇株型緊湊且生育后期群體結(jié)構(gòu)維持較好的玉米品種，通過增加密度可進一步增加夏玉米產(chǎn)量。調(diào)整播期增強灌漿初期光合效率也是有效的增產(chǎn)方法之一。

玉米；群體結(jié)構(gòu)；模型；光截獲; 密度；播期；產(chǎn)量

0 引言

【研究意義】群體結(jié)構(gòu)是影響玉米光合特性和產(chǎn)量的重要指標(biāo)之一，理想的群體結(jié)構(gòu)有助于光線在冠層合理分布，顯著提高作物群體光合特性和生物學(xué)產(chǎn)量[1-3]。群體結(jié)構(gòu)是一個復(fù)雜的綜合指標(biāo)，包括葉面積指數(shù)、葉面積分布、莖葉夾角、葉向值和株高等，且諸多指標(biāo)均隨生育進程的推進不斷變化[4]。群體結(jié)構(gòu)是一個動態(tài)變化的過程，任何一個指標(biāo)均可影響群體光分布的情況。另外，太陽高度和方位也是影響作物群體光空間分布的重要因素之一[5]。因此，運用傳統(tǒng)方法測定和描述群體結(jié)構(gòu)及其光分布特征是一個復(fù)雜、耗時且精確度難以控制的過程?！厩叭搜芯窟M展】基于試驗數(shù)據(jù)和模型模擬的研究手段能夠在一定程度上簡化群體結(jié)構(gòu)特征和光分布，同時提高效率和精確度?，F(xiàn)有一系列模型，例如群體結(jié)構(gòu)-功能模型[6-7]、3D結(jié)構(gòu)模型[8-9]和生長模型[10]等，均將光分布和光截獲作為重要的模擬部分，并試圖將不同的模型相結(jié)合，以此準(zhǔn)確地模擬群體光合特性。其中，群體葉面積指數(shù)或葉片空間分布是以上模型考慮的重要參數(shù)之一。例如，在群體結(jié)構(gòu)-功能模型中，群體葉片的空間分布，尤其是穗位層葉片可決定冠層總體的光合性能和物質(zhì)生產(chǎn)[11]。在3D結(jié)構(gòu)模型中，冠層不同層次的群體結(jié)構(gòu)，單個葉片的空間取向均是重點模擬的指標(biāo)，以期準(zhǔn)確地模擬群體光分布情況，并對群體結(jié)構(gòu)做出評價[9]。但是，上述模型均是模擬某個生育期或某個時間段的群體結(jié)構(gòu)，而對其進行整個生育期動態(tài)變化的模擬卻鮮有報道。在模擬群體光分布的過程中，太陽高度和方位隨玉米生育進程的動態(tài)變化和對玉米群體光分布的影響沒有重視。這些因素均是模擬群體整個生育期光合性能的重要因子，并最終影響產(chǎn)量模擬的準(zhǔn)確性。其次，優(yōu)化群體光合性能的目的是增加干物質(zhì)累積，促進產(chǎn)量的提高，但是上述模型沒有建立群體光分布與產(chǎn)量的具體關(guān)系，從而難以評價不同類型群體結(jié)構(gòu)及其光分布對產(chǎn)量的貢獻。同時，現(xiàn)有的群體光分布模型逐漸趨于復(fù)雜，所需參數(shù)較多，任何一個參數(shù)的微小變動均可影響模擬結(jié)果的準(zhǔn)確性和穩(wěn)定性，使模擬過程更加繁瑣?！颈狙芯壳腥朦c】如何在保證模擬結(jié)果準(zhǔn)確性的同時減少輸入?yún)?shù)，提高模型的簡便性和可操作性，對于指導(dǎo)農(nóng)業(yè)生產(chǎn)非常重要。在簡單的田間試驗或以往的試驗數(shù)據(jù)基礎(chǔ)上，借助簡易模型能夠快速地對玉米群體結(jié)構(gòu)做出評價，對于不同區(qū)域玉米品種選擇、播種密度、以及后期管理同樣有非常重要的意義?！緮M解決的關(guān)鍵問題】本研究利用玉米群體結(jié)構(gòu)模型與光分布模型，擬建立一個簡便可靠的玉米群體光截獲動態(tài)模型，以期對不同生育階段群體結(jié)構(gòu)做出評價。

1 材料與方法

1.1 模型設(shè)計

該模型將Chelle等[12]和Maddonni等[13]提出的光輻射模型與Wilson[14]和Welles等[15]提出的冠層結(jié)構(gòu)模型進行改進和整合。Maddonni等[13]提出群體有效光截獲公式（the fraction of incident photosynthetically active radiation intercepted by canopy，fIPAR）可表達為：

f=1--K×LAI（1）

式中，為消光系數(shù)，為葉面積指數(shù)。

Welles等[15]提出光線在群體中的穿透率為：

(,)=exp[-(,)(,)] （2）

式中，(,)為光輻射的群體結(jié)構(gòu)區(qū)域，為太陽天頂角，為太陽方位角，為葉片密度（每立方米群體中的葉片面積），(,)為光線在冠層中的路徑。在測量過程中，太陽方位角將被平均，可在公式中忽略不記，所以公式（2）可簡化為：

()=exp[-()()] （3）

由公式（1）、（3）可得，

()=1-f=1--K×LAI（4）

進而，由公式（3）、（4）可得，

exp[-()()]=-K×LAI（5）

最終得到公式（6）：

其中，在單一的群體結(jié)構(gòu)中，葉片密度與冠層高度和葉面積指數(shù)有關(guān)，可表示為：

（7）

光線在冠層群體中的路徑()與冠層高度和太陽天頂角相關(guān)，可表示為：

由公式（6）、（7）、（8）可知，消光系數(shù)可轉(zhuǎn)換為：

（9）

根據(jù)Wilson[13]和Welles等[14]提出的冠層結(jié)構(gòu)模型，光線在冠層群體中的輻射區(qū)域可表示為：

式中，葉莖夾角，為太陽天頂角。

由公式（9）、（10），將消光系數(shù)進一步轉(zhuǎn)換為：

因此，冠層群體光截獲模型最終表達為：

（12）

1.2 模型驗證與應(yīng)用

1.2.1 模型驗證 試驗一于2010—2011年在中國農(nóng)業(yè)大學(xué)吳橋?qū)嶒炚荆?7°41′N，116°37′E）進行，其數(shù)據(jù)用于模型建立及驗證。試驗采取裂區(qū)設(shè)計，品種和密度為主因素，氮肥為副因素，每個處理3次重復(fù)，小區(qū)面積為60 m2。供試品種為鄭單958（ZD 958）和丹玉405（DY 405）；種植密度為4.5×104，6.0×104，7.5×104株/hm2，行距均為60 cm；氮肥處理為不施氮、施氮180 kg·hm-2（分別在播前、大喇叭口期和吐絲期以1﹕4﹕1的比例施入）、施氮270 kg·hm-2（分別在播前、大喇叭口期和吐絲期以4﹕4﹕1的比例施入）。同時，施用磷肥（重過磷酸鈣）105 kg P2O5·hm-2、鉀肥（K2SO4）120 kg K2O·hm-2和ZnSO415 kg·hm-2作為基肥。本試驗于6月20播種，10月7日收獲。

1.2.2模型應(yīng)用 試驗二于2009年在中國農(nóng)業(yè)大學(xué)吳橋?qū)嶒炚具M行，其數(shù)據(jù)以及部分試驗一數(shù)據(jù)用于驗證模型的應(yīng)用效果。試驗設(shè)置6個播期，具體播種時間及相應(yīng)吐絲期和收獲期見表1，采取隨機區(qū)組設(shè)計，每個處理3次重復(fù)，小區(qū)面積為50 m2。供試品種為鄭單958，設(shè)置密度為8.25×104株/hm2，行距為60 cm；氮肥用量180 kg N·hm-2（于播前和大喇叭口期分別以1﹕2的比例施入），其他肥料用量和施入時間同試驗一。

表1 試驗二的播種、吐絲和收獲時間（月/日，2009）

1.2.3 測量指標(biāo)和數(shù)據(jù)來源 葉面積指數(shù)（LAI）：分別于吐絲期、灌漿期由傳統(tǒng)長寬法測量全株葉面積。每個小區(qū)取樣3株，單葉葉面積＝長×寬×0.75，單株葉面積為單葉葉面積的總和。LAI=單株葉面積×單位土地面積內(nèi)株數(shù)/單位土地面積。群體綠葉面積衰老速度用曲線方程=ae-cx/(1+e-cx)模擬，其中為灌漿期某一時刻相對吐絲期的綠葉面積（RGLA，%），為吐絲后天數(shù)，為RGLA的理論初始值（100%），為葉片衰老的啟動參數(shù)，為衰老參數(shù)[16]。該方程可用于模擬玉米灌漿期葉面積指數(shù)的動態(tài)變化。

葉莖夾角：于吐絲期利用量角器側(cè)量葉脈與莖邊緣的夾角，每個小區(qū)測量10株。

透光率：于吐絲期、灌漿期采用LAI-2000（Li-Cor，Lincoln，NE，USA）在群體底層（地面）測量所得，每個小區(qū)不同位置測量4次。

光截獲率（%）=1-透光率。

產(chǎn)量：每小區(qū)未取樣區(qū)域收獲玉米5行，每行4 m，稱量所收果穗總鮮重。按平均鮮穗重從所收果穗中隨機選取10穗，測定籽粒含水率，然后換算成實際產(chǎn)量（按14%折算含水率）。

太陽天頂角：數(shù)據(jù)來自NOAA ESRL（http://www. esrl.noaa.gov/gmd/grad/solcalc/azel.html）。同時測定吐絲期和灌漿期取樣當(dāng)天中午12點的太陽天頂角。

光合有效輻射：數(shù)據(jù)來自“中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)”泊頭站（http://cdc.cma.gov.cn/），站點編號為54 618，38°05′N, 116°33′E，海拔132 m，距吳橋試驗站約50 km。

光合有效輻射截獲量：群體光合有效輻射截獲量（MJ·m-2，下文用有效光截獲表述）=f×，灌漿期光合有效輻射累計量（MJ·m-2）=?(day)，灌漿期群體光合有效輻射累計截獲量（MJ·m-2）=? f(day)×(day)。式中，f為群體截光率（%），f(day)為灌漿期每日群體光截獲率，為光合有效輻射量（MJ·m-2），(day)為灌漿期每日光合有效輻射量，f和f(day)均為模擬值。

1.3 數(shù)據(jù)處理與分析

本研究運用SAS 9.4程序“Proc. Corr.”分析模型模擬值與田間試驗測量值之間的相關(guān)性以及產(chǎn)量與光截獲的相關(guān)性，整個模型在Excel 2010中運行。

2 結(jié)果

2.1 光截獲模型的驗證

為驗證模型的準(zhǔn)確性，使用試驗一的研究結(jié)果分 別對玉米吐絲期和灌漿中期各處理光截獲率的測量值與模擬值進行相關(guān)分析（圖1）。結(jié)果表明，兩者呈顯著相關(guān)關(guān)系，且吐絲期相關(guān)系數(shù)高于灌漿中期，分別為0.91和0.85。

**, P＜0.01顯著。下同

2.2 模型應(yīng)用

2.2.1 密度與群體光截獲的關(guān)系 施氮量為180 kg·hm-2（播前﹕大喇叭口期﹕吐絲期=1﹕4﹕1）是吳橋地區(qū)較為常用的施氮方式，因此以本處理來反映吐絲后群體光截獲隨生育進程的變化情況。不同密度下，兩品種群體光截獲率（圖2）和光截獲量（圖3）在整個灌漿期均呈現(xiàn)先升高后降低的趨勢。光截獲率在吐絲后30 d達到峰值。由于受陰雨寡照影響，有效光輻射在吐絲后20—30 d處于較低水平，導(dǎo)致光截獲量波動較大。群體光截獲率和截獲量均在吐絲后50 d開始加速衰減。在4.5×104株/hm2密度下，光截獲率（量）明顯低于6.0×104和7.5×104株/hm2種植密度；而在6.0×104和7.5×104株/hm2密度下光截獲量非常接近（圖2—3）。

2.2.2 播期與群體光截獲的關(guān)系 利用已經(jīng)建立的模型，對試驗二播期試驗的數(shù)據(jù)進行分析整理。模擬結(jié)果顯示，不同播期群體有效光截獲率（量）在整個灌漿期呈現(xiàn)先升高后降低的趨勢（圖4），其中，光截獲率隨播期的推遲逐漸升高，即播期1—3群體處于較低的光截獲水平，而播期4—6群體的光截獲水平則較高，且播期4—6的數(shù)值非常接近（圖4-A）。受不同階段有效光輻射的影響，雖然整體規(guī)律表現(xiàn)一致（先升后降），但單個播期群體在灌漿期不同階段的光截獲量（MJ·m-2）波動較大（圖4-B）。整個灌漿期光合有效輻射累積量（MJ·m-2）隨播期推遲呈逐漸降低的趨勢，而光合有效輻射累積截獲量（MJ·m-2）隨播期推遲呈先升高后降低的趨勢，其中播期3的數(shù)值最高（圖5）。

圖2 灌漿期鄭單958（ZD 958）和丹玉405（DY 405）在3個密度下的群體光截獲率模擬值

圖3 灌漿期鄭單958（ZD 958）和丹玉405（DY 405）在3個密度下的群體光合有效輻射截獲量模擬

圖4 不同播期下鄭單958群體灌漿期光合有效輻射截獲率模擬值（A）與十天平均光截獲量模擬值（B）隨生育進程的變化

圖5 灌漿期不同播期群體光合有效輻射累計量與累計截獲量

2.2.2 產(chǎn)量與群體光截獲的關(guān)系 運用本模型對兩個不同試驗的產(chǎn)量和光合有效輻射數(shù)據(jù)進行綜合分析后發(fā)現(xiàn)，籽粒產(chǎn)量與光合有效輻射均呈正相關(guān)關(guān)系（圖6—7）。其中，播期試驗（試驗二）的相關(guān)性高于密度和氮肥試驗（試驗一）。由圖6可知，產(chǎn)量與灌漿期光合有效輻射累計量的相關(guān)性（0.86）高于與光輻射截獲量的相關(guān)性（0.80），且均達到顯著相關(guān)。由圖7可知，產(chǎn)量與灌漿期不同階段（10 d平均值，試驗一）群體有效光輻射截獲量同樣呈現(xiàn)正相關(guān)性，但相關(guān)系數(shù)較低且不顯著，并隨灌漿期的推進逐漸降低。鄭單958相關(guān)系數(shù)的變動范圍為0.31—0.21，丹玉405為0.55—0.33，后者在不同灌漿階段均高于前者。

3 討論

3.1 模型驗證與分析

由不同時期群體光截獲率模擬值與測量值的顯著相關(guān)性可知（圖1），該模型的準(zhǔn)確性已達到較為理想的水平，可以用于評價玉米群體結(jié)構(gòu)。對吐絲期和灌漿中期進行比較發(fā)現(xiàn)，測量值均高于模擬值，其主要原因是測量值未能排除底部衰老葉片對光截獲的影響[4, 17]。底部衰老或干枯葉片可以攔截部分光照，從而影響光截獲的測量值。兩者的相關(guān)性在灌漿中期低于吐絲期，有以下四點原因：隨著生育期的推進，葉片出現(xiàn)下披的現(xiàn)象，導(dǎo)致葉莖夾角略有增大[18]，由于模型中未涉及此參數(shù)的校正導(dǎo)致模擬值略有偏差；隨著底部衰老葉片增多，其對光截獲率測量值的影響增大，而在模型中葉面積參數(shù)是有效綠葉面積更接近真實狀態(tài)；未考慮雄穗對光輻射的截獲作用[19-20]，在一定程度上可能影響模型的準(zhǔn)確性；另外，模型中葉面積指數(shù)為模擬值，并非每個生育期的真實測量值，因此造成了一定的模擬偏差?？傮w而言，與群體光截獲的實測方法相比，該模擬方法不僅節(jié)省大量的時間和人力，而且可通過對群體不同生育期光截獲的模擬分析，對群體結(jié)構(gòu)的動態(tài)變化進行評價。

圖6 籽粒產(chǎn)量與灌漿期光合有效輻射累積量和累積截獲量的相關(guān)性(試驗二)

圖7 不同品種群體產(chǎn)量與灌漿期不同階段光合有效輻射截獲量的相關(guān)性

與現(xiàn)有的玉米群體光截獲模型相比[6-10, 21]，本模型將群體葉面積作為一個整體，模擬了整個群體的光截獲狀態(tài)，簡化了模擬步驟，省略了大量參數(shù)，使模擬過程更加簡潔，并且實現(xiàn)了群體光截獲隨生育期變化的動態(tài)模擬。在接下來的研究中，可以進一步細化光輻射在群體不同層次的分布，為深入研究光照、不同部位葉片對干物質(zhì)生產(chǎn)及產(chǎn)量的影響奠定基礎(chǔ)。

3.2 模型應(yīng)用及分析

光照在群體內(nèi)的分布和截獲是評價作物群體結(jié)構(gòu)的重要參數(shù)，但至今仍沒有研究提出針對不同群體結(jié)構(gòu)的光截獲評價范圍[6, 22-23]，也鮮有光截獲模型涉及籽粒產(chǎn)量[4]?？梢娋_模擬群體光截獲具有一定難度，主要因為涉及群體光截獲的影響因素較多，例如品種、密度、播期等。品種是調(diào)控群體結(jié)構(gòu)的重要措施之一，本文中丹玉405比鄭單958更為緊湊（18.0°和22.4°，吐絲期整株平均值），在高密度下有進一步提高光截獲的潛力（圖2—3）。但是，當(dāng)密度大于6×104株/hm2時，群體光截獲率（量）變化很小，尤其是華北平原地區(qū)最常用的品種鄭單958，表明在當(dāng)?shù)赝ㄟ^增加密度進一步提升產(chǎn)量的空間有限，這點與徐麗娜等[24]的研究結(jié)果一致。在播期試驗中，產(chǎn)量與群體光輻射截獲量的相關(guān)性處于較低的水平=0.8（圖6），表明品種鄭單958在8.25×104株/hm2（播期試驗的種植密度）密度下的群體結(jié)構(gòu)并不理想；可能由于群體結(jié)構(gòu)過于郁閉，截獲的有效光輻射不能被充分利用。另外，產(chǎn)量與群體的光合能力（群體功能）、后期物質(zhì)轉(zhuǎn)運同樣密切相關(guān)[25]，這些因素也可影響產(chǎn)量與群體光輻射的相關(guān)性。

通過早播或者晚收能夠提高灌漿期光合有效輻射累計量，從而達到高產(chǎn)的目的[26-27]。本研究中播期早于5月20日的群體光截獲率處于較低的水平，說明群體結(jié)構(gòu)的優(yōu)勢并沒有得到充分發(fā)揮；播期晚于6月4日的群體光截獲率則較高，且不同播期間數(shù)值非常接近。上述結(jié)果說明，晚于6月4日播種更利于群體結(jié)構(gòu)功能的充分發(fā)揮。鑒于當(dāng)?shù)匾荒陜墒斓母髦贫龋ǘ←?夏玉米），夏玉米整個生育期被控制在6月中下旬至9月底，即100 d左右[27-29]。夏玉米播種太早將影響冬小麥的產(chǎn)量，而播種太晚會顯著降低夏玉米的產(chǎn)量[30]。因此，大幅度調(diào)整夏玉米播期的可能性不大。另外，吐絲開花期高溫、陰雨等逆境條件同樣制約著夏玉米播期的選擇[26]。適當(dāng)早播（6月中上旬）在提高群體光輻射利用效率的同時，也提高了避免高溫、寡照等逆境的可能性。

由不同階段光截獲與產(chǎn)量的相關(guān)性可知（圖7），灌漿初期光照對產(chǎn)量的影響大于灌漿的其他時期，可能與當(dāng)?shù)赝陆z期前后的陰雨寡照環(huán)境有關(guān)。徐麗娜等[24]、陳建忠等[31]指出在該試驗區(qū)域（黑龍港流域）7—8月的降雨量可達全年的70%—80%，日照時數(shù)為5.8 h左右。這段時間正值夏玉米吐絲、灌漿初期，“源”的物質(zhì)生產(chǎn)能力不足，將直接影響夏玉米產(chǎn)量的進一步增加。

群體光截獲與產(chǎn)量的相關(guān)性在整個灌漿期一直處于較低的水平，表明光照不足只是制約當(dāng)?shù)叵挠衩自霎a(chǎn)的因素之一。由于灌漿后期群體結(jié)構(gòu)維持功能衰減加快（吐絲后50 d左右，圖2）以及光輻射量的減小，灌漿初期的相關(guān)性大于灌漿后期，表明在灌漿中后期物質(zhì)生產(chǎn)對產(chǎn)量的貢獻逐漸減弱，而物質(zhì)轉(zhuǎn)移隨生育期進行變得愈加重要。

4 結(jié)論

通過利用玉米群體結(jié)構(gòu)模型與光分布模型建立了簡便有效的玉米群體光截獲動態(tài)模型，該模型準(zhǔn)確的描述玉米冠層的光輻射截獲能力，并可結(jié)合當(dāng)?shù)氐墓廨椛洵h(huán)境，對吐絲后整個生育期的群體光截獲做出動態(tài)模擬，細化光輻射在群體內(nèi)的分布。該模型結(jié)合吳橋地區(qū)大田試驗分析發(fā)現(xiàn)，根據(jù)不同玉米品種形態(tài)特征，該試驗地區(qū)夏玉米最適種植密度為6.0×104—7.5×104株/hm2，最適播期為6月中上旬。在現(xiàn)有的耕作措施和品種條件下，通過增加密度進一步增加夏玉米產(chǎn)量的空間有限，但通過調(diào)整播期增強灌漿初期光合效率，或選擇株型緊湊且能在生育后期維持較好群體結(jié)構(gòu)的品種是有效的增產(chǎn)方法。

致謝：感謝徐麗娜、朱金城博士在數(shù)據(jù)分析與處理方面給予的幫助與支持！

[1] 董樹亭, 胡昌浩, 岳壽松, 王群瑛, 高榮岐, 潘子龍.夏玉米群體光合速率特性及其與冠層結(jié)構(gòu)、生態(tài)條件的關(guān)系. 植物生態(tài)學(xué)與地植物學(xué)學(xué)報, 1992, 16(4): 372-378.

DONG S T, HU C H, YUE S S, WANG Q Y, GAO R Q, PAN Z L. The characteristics of canopy photosynthesis of summer corn and its relationship with canopy structure and ecological conditions., 1992, 16(4): 372-378. (in Chinese)

[2] 呂麗華, 陶洪斌, 夏來坤, 張雅杰, 趙明, 趙久然, 王璞. 不同種植密度下的夏玉米群體結(jié)構(gòu)及光合特性. 作物學(xué)報, 2008, 34(3): 447-455.

Lü L H, TAO H B, XIA L K, ZHANG Y J, ZHAO M, ZHAO J R, WANG P. Canopy structure and photosynthesis traits of summer maize under different planting densities., 2008, 34(3): 447-455. (in Chinese)

[3] 薛吉全, 梁宗鎖, 馬國勝, 路海東, 任建宏. 玉米不同株型耐密性的群體生理指標(biāo)研究. 應(yīng)用生態(tài)學(xué)報, 2002, 13(1): 55-59.

XUE J Q, LIANG Z S, MA G S, LU H D, REN J H. Population physiological indices on density-tolerance of maize in different plant type., 2002, 13(1): 55-59. (in Chinese)

[4] BIRCH C J, ANDRIEU B, FOURNIER C, VOS J, ROOM P. Modelling kinetics of plant canopy architecture concepts and applications., 2003, 19(4): 519-533.

[5] CHELLE M, EVERS J B, COMBES D, VARLET-GRANCHER C, VOS J, ANDRIEU B. Simulation of the three-dimensional distribution of the red: Far-red ratio within crop canopies., 2007, 176(1): 223-234.

[6] SARLIKIOTI V, DE VISSER P H, MARCELIS L F. Exploring the spatial distribution of light interception and photosynthesis of canopies by means of a functional-structural plant model., 2011, 107(5): 875-883.

[7] 王謙, 陳景玲, 孫治強. 用LI-2000冠層分析儀確定作物群體外紅多面高度. 農(nóng)業(yè)工程學(xué)報, 2005, 21(8): 70-74.

WANG Q, CHEN J L, SUN Z Q. Determination of the height of active surface of the crop colonies using Li-2000 Canopy Analyzer., 2005, 21(8): 70-74. (in Chinese)

[8] DROUET J L. MODICA and MODANCA: Modelling the three- dimensional shoot structure of graminaceous crops from two methods of plant description., 2003, 83(2): 215-222.

[9] 孔婭, 勞彩蓮, 曹素云. 利用3D 模型模擬天空與葉面散射對玉米冠層截光率的影響. 農(nóng)業(yè)工程學(xué)報, 2011, 27(5): 248-252.

KONG Y, LAO C L, CAO S Y. Effect of sky radiation and leaf scattering on maize canopy light interception by 3D modeling., 2011, 27(5): 248-252. (in Chinese)

[10] DROUET J L, PAGES L. GRAAL: A model of GRowth, Architecture and carbon ALlocation during the vegetative phase of the whole maize plant : Model description and parameterisation., 2003, 165(2/3): 147-173.

[11] PALMER A F E, HEICHEL G H, MUSGRAVE R B. Patterns of translocation, respiratory loss, and redistribution of14C in maize labeled after flowering., 1973, 13(3): 371-376.

[12] CHELLE M, ANDRIEU B. The nested radiosity model for the distribution of light within plant canopies., 1998, 111(1): 75-91.

[13] MADDONNI G A, CHELLE M, DROUET J L, ANDRIEU B. Light interception of contrasting azimuth canopies under square and rectangular plant spatial distributions: Simulations and crop measurements., 2001, 70(1): 1-13.

[14] WILSON J W. Inclined point quadrats., 1960, 59(1): 1-7.

[15] WELLES J, NORMAN J. Instrument for indirect measurement of canopy architecture., 1991, 83(5): 818-825.

[16] 劉開昌, 董樹亭, 趙海軍, 王慶成, 李宗新, 劉霞, 張慧. 我國玉米自交系葉片保綠性及其與產(chǎn)量的關(guān)系. 作物學(xué)報, 2009, 35(9): 1662-1671.

LIU K C, DONG S T, ZHAO H J, WANG Q C, LI Z X, LIU X, ZHANG H. Leaf stay-green traits in Chinese maize inbred lines and their relationship with grain yield., 2009, 35(9): 1662-1671. (in Chinese)

[17] 張子山, 楊程, 高輝遠, 李耕, 劉鵬. 保綠玉米與早衰玉米葉片衰老過程中葉綠素降解與光合作用光化學(xué)活性的關(guān)系. 中國農(nóng)業(yè)科學(xué), 2012, 45(23): 4794-4800.

ZHANG Z S, YANG C, GAO H Y, LI G, LIU P. Relationship between photochemistry activity and decrease in chlorophyll content during senescence in leaves of stay green and quick-leaf-senescence inbred line of maize., 2012, 45(23): 4794-4800. (in Chinese)

[18] TIAN F, BRADBURY P J, BROWN P J, HUNG H, SUN Q, FLINT-GARCIA S, ROCHEFORD T R, MCMULLEN M D, HOLLAND J B, BUCKLER E S. Genome-wide association study of leaf architecture in the maize nested association mapping population., 2011, 43(2): 159-162.

[19] 黃華國, 竇寶成, 胡妮. 雄穗對玉米冠層熱輻射方向性的影響分析. 紅外與毫米波學(xué)報, 2011, 30(2): 120-123.

HUANG H G, DOU B C, HU N. Tassel effect on the thermal emission directionality of corn canopies., 2011, 30(2): 120-123. (in Chinese)

[20] ANDRIEU B, IVANOV N, BOISSARD P. Simulation of light interception from a maize canopy model constructed by stereo plotting., 1995, 75(1): 103-119.

[21] CABRERA-BOSQUET L, FOURNIER C, BRICHET N, WELCKER C, SUARD B, TARDIEU F. High-throughput estimation of incident light, light interception and radiation-use efficiency of thousands of plants in a phenotyping platform., 2016, 212(1): 269-281.

[22] BORRAS L, MADDONNI G A, OTEGUI M E. Leaf senescence in maize hybrids: Plant population, row spacing and kernel set effects., 2003, 82(1): 13-26.

[23] Koetz B, Baret F, Poilvé H, Hill J. Use of coupled canopy structure dynamic and radiative transfer models to estimate biophysical canopy characteristics., 2005, 95(1): 115-124.

[24] 徐麗娜, 陶洪斌, 黃收兵, 明博, 王璞. 黑龍港流域夏玉米產(chǎn)量提升限制因素. 生態(tài)學(xué)報, 2013, 33(1): 222-228.

XU L N, TAO H B, HUANG S B, MING B, WANG P. Analysis on the limiting factors to further improve yield of summer maize in Heilonggang River Valley., 2013, 33(1): 222-228. (in Chinese)

[25] ZHOU B, YUE Y, SUN X, WANG X, WANG Z, MA W, ZHAO M. Maize grain yield and dry matter production responses to variations in weather conditions., 2016, 108(1): 196-204.

[26] 李向嶺, 李從鋒, 侯玉虹, 侯海鵬, 葛均筑, 趙明. 不同播期夏玉米產(chǎn)量性能動態(tài)指標(biāo)及其生態(tài)效應(yīng). 中國農(nóng)業(yè)科學(xué), 2012, 45(6): 1074-1083.

Li X L, Li C F, Hou Y H, Hou H P, Ge J Z, Zhao M. Dynamic characteristics of summer maize yield performance in different planting dates and its effect of ecological factors., 2012, 45(6): 1074-1083. (in Chinese)

[27] 明博, 朱金城, 陶洪斌, 徐麗娜, 郭步慶, 王璞. 黑龍港流域玉米不同生育階段氣象因子對產(chǎn)量性狀的影響. 作物學(xué)報, 2013, 39(5): 919-927.

MING B, ZHU J C, TAO H B, XU L N, GUO B Q, WANG P. Effects of meteorological factors at different growth stages on yield traits of maize (L) in Heilonggang Basin., 2013, 39(5): 919-927. (in Chinese)

[28] 付雪麗, 張惠, 賈繼增, 杜立豐, 付金東, 趙明. 冬小麥-夏玉米“雙晚”種植模式的產(chǎn)量形成及資源效率研究. 作物學(xué)報, 2009, 35(9): 1708-1714.

FU X L, ZHANG H, JIA J Z, DU L F, FU J D, ZHAO M. Yield performance and resources use efficiency of winter wheat and summer maize in double late-cropping system., 2009, 35(9): 1708-1714. (in Chinese)

[29] WANG J, WANG E, YANG X, ZHANG F, YIN H. Increased yield potential of wheat-maize cropping system in the North China Plain by climate change adaptation., 2012, 113(3/4): 825-840.

[30] 薛慶禹, 王靖, 曹秀萍, 馬薇, 馮利平. 不同播期對華北平原夏玉米生長發(fā)育的影響. 中國農(nóng)業(yè)大學(xué)學(xué)報, 2012, 17(5): 30-38.

XUE Q Y, WANG J, CAO X P, MA W, FENG L P. Effect of sowing date and variety on growth and population characteristics of summer maize in North China Plain., 2012, 17(5): 30-38. (in Chinese)

[31] 陳建忠, 肖荷霞. 黑龍港流域氣象生態(tài)因子對夏玉米穗粒數(shù)的影響. 河北農(nóng)業(yè)科學(xué), 2000, 4(3): 23-27.

Chen J Z, Xiao H X. Effects of ecoclimatical factors on grain number of summer maize in Heilonggang Area., 2000, 4(3): 23-27. (in Chinese)

（責(zé)任編輯 楊鑫浩）

Dynamic Simulation of Relationship Between Light Interception and Growth of Maize Population and Its Application

CAO YiBing, Huang ShouBing, Wang YuanYuan, Xia YuQing, Meng QingFeng, Tao HongBin, Wang Pu

(College of Agronomy and Biotechnology, China Agricultural University, Beijing 100193)

【Objective】The population light interception (PLI) has a great impact on growth, development, and grain production of maize (L.) in interaction with canopy structure of maize population. It costs amounts of time and labor to measure the morphology index in the evaluation of maize canopy structure. And it is impossible to integrate solar radiation into the evaluation of canopy structure. Therefore, it is quite necessary to develop a simple but feasible model that combine canopy structure and solar radiation to improve the evaluation of maize canopy structure, and thus to guide farmers in crop management.【Method】Based on the canopy structure model and light distribution model, a simple and reliable light interception model of maize canopy was developed. In order to evaluate the prediction accuracy of this model, two field trials were conducted at Wuqiao Experimental Station of China Agricultural University in 2009-2011. Trial 1 was set up as a split-plot design with varieties (ZD 958 and DY 405) and plant densities (45 000, 60 000, and 75 000 plant/hm2) as main factors, and nitrogen management (No nitrogen; 180 kg·hm-2N applied at sowing﹕13-leaf stage﹕silking stage=1﹕4﹕1; 270 kg·hm-2N applied at sowing﹕13-leaf stage﹕silking stage =4﹕4﹕1) as second factor. Trial 2 was in a randomized complete block design, including one variety (ZD 958), one plant density (82 000 plant/hm2), and six sowing dates (April 20, May 5 and 20, June 4 and 18, and July 3, 2009). 【Result】The simulated PLI rates were significantly correlated with measured PLI rates at silking stage and mid kernel-filling stage, withvalues of 0.91 and 0.85, respectively. During the entire kernel-filling stage, the PLI values first increased and then decreased, and the reduction became more obvious on 50 days after silking. Moreover, the PLI values remained at the highest level when plant density was greater than 60 000 plant/hm2. Additionally, population structure under sowing date of June 4 achieved the highest PLI, demonstrating its great potential of structure-function for the yield improvement. PLI and yield were positively correlated, while their correlation coefficient declined over time during grain-filling period, indicating that the dry matter allocation might play more important roles on yield during the late grain-filling period. The decline of correlation coefficient also depended on variety; variety with erecter leaves could remain a higher coefficient value over time. 【Conclusion】The potential to further improve maize yield through increasing plant density (＞75 000 plant/hm2) is limited at experimental site. Some strategies should be concerned, such as (1) applying varieties with erect leaves and long-lasting canopy structure persistence during late grain-filling period, and (2) adjusting sowing date according to climate situation to increase population photosynthesis at early grain-filling period.

maize (L.); canopy structure; model; light interception; plant density; sowing date; yield

2016-08-25；

2017-03-14

國家玉米產(chǎn)業(yè)技術(shù)體系（201203031）、“糧食豐產(chǎn)增效科技創(chuàng)新”重點專項（SQ2017YFNC050062）

黃收兵，E-mail：huangshoubing@cau.edu.cn

聯(lián)系方式：曹亦兵，E-mail：784298412@qq.com。