王振昌，郭相平，吳夢(mèng)洋，欒雅珺，陳 盛，邱讓建，劉春偉，操信春，虞穎蕾，曹 宇(.河海大學(xué)水利水電學(xué)院，南京 0098；.南京信息工程大學(xué)江蘇省農(nóng)業(yè)氣象重點(diǎn)實(shí)驗(yàn)室，南京 0044)

株高是水稻重要生長(zhǎng)指標(biāo)之一，也是構(gòu)成水稻理想株型的重要因素，適宜的株高是協(xié)調(diào)產(chǎn)量結(jié)構(gòu)、群體光合效率和抗倒伏能力的重要性狀[1,2]。株高除受遺傳因素影響外，也受外界環(huán)境的影響。前人研究表明，施肥量、CO2濃度、O3濃度以及地溫變化等因素均會(huì)影響水稻株高等生長(zhǎng)指標(biāo)[3-6]。除上述因素外，水稻不同生育階段進(jìn)行的灌溉處理，也能顯著影響水稻株高等生長(zhǎng)性狀[4,7,8]。在水稻灌排實(shí)踐中，人們主要關(guān)注各生育階段的適宜耐淹深度和耐淹(旱)歷時(shí)對(duì)水稻產(chǎn)量、品質(zhì)以及環(huán)境效益的影響[9-12]，然而對(duì)于不同灌溉方式，尤其是不同生育階段旱澇交替脅迫處理下水稻的抗倒伏性能研究較少[13]。水稻隨著株高以及地上鮮重的增大，基部莖節(jié)的彎曲力矩變大，水稻倒伏幾率增大[13-16]，嚴(yán)重影響水稻的產(chǎn)量、品質(zhì)及收割效率[14]。研究不同生育階段旱澇交替脅迫處理對(duì)水稻株高的動(dòng)態(tài)變化，以及應(yīng)用生長(zhǎng)曲線模型對(duì)株高生長(zhǎng)進(jìn)行模擬，可為估算水稻生長(zhǎng)后期的倒伏風(fēng)險(xiǎn)提供數(shù)據(jù)支持。近年來(lái)研究人員已構(gòu)建了多種用于對(duì)植物生長(zhǎng)曲線進(jìn)行描述的數(shù)學(xué)模型, 并發(fā)現(xiàn)多數(shù)植物的生長(zhǎng)均符合 Logistic, Gompertz 和 von Bertalanffy 等非線性模型，且上述模型多用于模擬不同品種、基質(zhì)、氮素水平以及種植密度條件下植物的生長(zhǎng)動(dòng)態(tài)[16-21], 而對(duì)于單生育階段的旱澇交替脅迫對(duì)水稻株高生長(zhǎng)的模擬鮮見(jiàn)報(bào)道。前人研究表明，一定生育階段適度的水分脅迫可使復(fù)水后短期內(nèi)作物的生理生長(zhǎng)得以改善[22-25]。然而，植物生長(zhǎng)曲線模型是否適合模擬不同生育階段旱澇交替脅迫及其補(bǔ)償效應(yīng)產(chǎn)生的株高生長(zhǎng)動(dòng)態(tài)變化，尚未見(jiàn)報(bào)道。因此本文通過(guò)采用盆栽試驗(yàn)，研究水稻株高在分蘗期、拔節(jié)期旱澇交替脅迫下的動(dòng)態(tài)變化，并通過(guò) Logistic, Gompertz與von Bertalanffy 3種典型的非線性生長(zhǎng)曲線模型對(duì)常規(guī)的淺水勤灌方式及不同生育階段旱澇交替脅迫下水稻生長(zhǎng)規(guī)律進(jìn)行生長(zhǎng)曲線擬合，并經(jīng)過(guò)分析不同灌溉處理的株高生長(zhǎng)速度以及不同旱澇脅迫的生長(zhǎng)差異，揭示水稻在不同灌溉方式下株高生長(zhǎng)規(guī)律，優(yōu)選出適宜旱澇交替脅迫處理下的株高生長(zhǎng)模型，以期為進(jìn)一步模擬不同水分調(diào)控方式對(duì)水稻株高影響，以及估算水稻生長(zhǎng)后期的倒伏風(fēng)險(xiǎn)提供理論及數(shù)據(jù)參考。

1 模型選擇與數(shù)據(jù)分析

1.1 模型選擇

植物生長(zhǎng)曲線法作為趨勢(shì)外推法的一種重要方法，在描述及預(yù)測(cè)生物個(gè)體的生長(zhǎng)發(fā)育等領(lǐng)域中已得到廣泛應(yīng)用[19,20]。常用的數(shù)學(xué)模型有Logistics、Gompertz 和von Bertaffany 等模型[21,22]，其中Logistics和Gompertz 模型拐點(diǎn)固定，且均具有飽和增長(zhǎng)特征，而von Bertalanffy模型的拐點(diǎn)可變[11,12,23]。Logistic, Gompertz與von Bertalanffy等3種非線性模型的表達(dá)式及特征參數(shù)如表1所示。

表1 非線性生長(zhǎng)模型及特征參數(shù)Tab.1 Three types of nonlinear models and characteristic parameters for fitting the growth curve

注：t為水稻移栽后天數(shù)；Y是在時(shí)間t的株高；A為水稻最大株高；K為相對(duì)生長(zhǎng)速率參數(shù)；B為常規(guī)參數(shù)。

1.2 試驗(yàn)設(shè)計(jì)及株高數(shù)據(jù)測(cè)定

本試驗(yàn)于2013年5-10月在河海大學(xué)南方地區(qū)高效灌排與農(nóng)業(yè)水土環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室玻璃溫室內(nèi)進(jìn)行(南京, 31°57′N, 118° 50′E, 海拔14.4 m)，試驗(yàn)水稻品種為當(dāng)?shù)馗弋a(chǎn)品種“南粳44”(Oryza sativa L. Nanjing44)。試驗(yàn)桶高為25.5 cm，上部和底部直徑分別為 24和18 cm。試驗(yàn)用土采自附近稻田耕作層，土壤類型為黏壤土，土壤中速效氮與速效磷含量分別為50.2 和 9.58 mg/kg。所用土壤經(jīng)曬干、打碎、過(guò)篩和風(fēng)干后，均勻施肥，施肥量為每千克風(fēng)干土尿素[CO(NH2)2]0.22 g、硫酸鉀(K2SO4)0.09 g、磷酸二氫鉀(KH2PO4)0.25 g、有機(jī)肥1.67 g，其中有機(jī)肥氮素含量質(zhì)量分?jǐn)?shù)為3.75%，裝桶后土壤容重為1.29 g/cm3，田間持水率(Field capacity，F(xiàn)C)為 33.0%(質(zhì)量含水率)。試驗(yàn)于2013年5月8日育秧，6月21日選擇大小基本一致的3葉1心的秧苗移栽，每桶3穴，每穴移栽秧苗2株。本試驗(yàn)共設(shè)置5個(gè)灌溉處理，分別為對(duì)照處理(CK)，分蘗期輕旱-澇-輕旱處理(T-LD)，分蘗期重旱-澇-重旱處理(T-HD)，拔節(jié)期輕旱-澇-輕旱處理(J-LD)和拔節(jié)期重旱-澇-重旱處理(J-HD)。試驗(yàn)中對(duì)照處理(CK)為淺水勤灌，除黃熟期外，各生育階段均保持5 cm水層不變。試驗(yàn)中干旱脅迫通過(guò)人工稱重法實(shí)現(xiàn)：在試驗(yàn)期間，通過(guò)一天2次的頻率(07∶00和18∶00 各稱量1 次)稱取試驗(yàn)桶重，當(dāng)土壤含水量低于或接近于灌水下限時(shí)，人工灌水至灌水上限，維持桶內(nèi)土壤含水量處于相應(yīng)生育階段的灌水上限和灌水下限之間；試驗(yàn)中的澇脅迫通過(guò)將水稻試驗(yàn)桶統(tǒng)一放入對(duì)應(yīng)水位的水箱來(lái)實(shí)現(xiàn)。水分脅迫開(kāi)始前，各處理均為淺水勤灌。水稻分蘗期旱脅迫于2013年7月21日開(kāi)始，歷時(shí) 5 d，于7月26日結(jié)束，隨即由旱脅迫轉(zhuǎn)入澇脅迫，并持續(xù)5 d，于8月7日再次轉(zhuǎn)入旱脅迫，8月12日分蘗期脅迫結(jié)束，后續(xù)轉(zhuǎn)入淺水勤灌；水稻拔節(jié)期旱脅迫于2013年8月20日開(kāi)始，脅迫歷時(shí)5 d；8月25日由旱脅迫轉(zhuǎn)入澇脅迫，脅迫歷時(shí)5 d；9月2日再次轉(zhuǎn)入旱脅迫，9月7日拔節(jié)期脅迫結(jié)束后隨即轉(zhuǎn)入淺水勤灌，并于10月29日(DAT 130)收獲水稻。試驗(yàn)水分處理旱脅迫分為重旱(土壤含水率上下限分別為 60%FC和70%FC)和輕旱(土壤含水率上下限分別為70%FC和80%FC)，澇脅迫為輕澇(在分蘗期和拔節(jié)期分別保持桶中土面距水面10和15 cm)。水稻各水分控制方案如圖1所示。上述各灌溉處理隨機(jī)排列，每種灌溉處理設(shè)40個(gè)重復(fù)。水稻植株移栽后，分別于DAT 24、35、40、51、59、65、70、77、94、108、118和130，隨機(jī)抽取6株水稻用直尺測(cè)定株高。水稻抽穗前株高為土面至每叢最高葉尖的高度，抽穗后為土面至最高穗頂(不連芒)的高度。

圖1 灌溉處理示意圖Fig.1 Diagram of experimental setup

1.3 數(shù)據(jù)統(tǒng)計(jì)與參數(shù)擬合

利用SPSS 13.0 (SPSS, Chicago, IL, USA) 進(jìn)行方差分析與鄧肯多重比較對(duì)不同灌溉處理的株高數(shù)據(jù)進(jìn)行分析。選用表1中3種非線性增長(zhǎng)模型，取各灌溉方式下65%的實(shí)際測(cè)定株高，利用SPSS 13.0 中的Nonlinear Regression功能，通過(guò)比較擬合優(yōu)度的擬合優(yōu)度最大估算法分別對(duì)不同灌溉處理水稻株高增長(zhǎng)過(guò)程進(jìn)行擬合，求出不同灌溉方式相應(yīng)模型下的表達(dá)式及特征參數(shù)，建立不同生育階段旱澇脅迫下水稻株高生長(zhǎng)曲線模型。

1.4 模型評(píng)價(jià)標(biāo)準(zhǔn)

取65%的株高數(shù)據(jù)用于參數(shù)估計(jì)后，取剩下的35%的株高數(shù)據(jù)用于上述3種非直線模型實(shí)測(cè)值和模擬值的比較。修正效率系數(shù)E1、修正一致性指標(biāo)d1和平均絕對(duì)誤差MAE為3種模型評(píng)價(jià)指標(biāo)[26]。E1反應(yīng)模型預(yù)測(cè)值和實(shí)測(cè)值的符合程度，其值在-∞與1之間，數(shù)值越大表明符合程度越好，d1的值在0與1之間，數(shù)值越大，表明其一致性優(yōu)越。上述3種模型評(píng)價(jià)指標(biāo)的計(jì)算方法分如下：

(3)

2 結(jié)果分析與模型擬合

2.1 不同生育階段旱澇交替脅迫下水稻株高動(dòng)態(tài)變化分析

由圖2知，在分蘗期第一次旱脅迫結(jié)束時(shí)(DAT 35)，T-LD和T-HD的株高從40 cm(DAT24)分別增長(zhǎng)到48.75和47.5 cm，其中干旱脅迫重的T-HD比脅迫輕的T-LD處理株高低了1.05 cm，經(jīng)統(tǒng)計(jì)分析，兩者無(wú)顯著差異(p>0.05)，但兩處理均值均顯著低于無(wú)干旱脅迫的CK處理(p<0.01)；其后的澇脅迫處理 (水層深度持續(xù)保持10 cm)，使得該時(shí)期T-HD的增長(zhǎng)速度顯著大于T-LD處理，但澇脅迫5 d后(DAT40)T-LD和T-HD的平均株高仍顯著低于CK處理(p<0.01)；隨后的澇轉(zhuǎn)旱過(guò)程以及其后5d的旱脅迫，使得T-HD和T-LD處理的株高直至分蘗期結(jié)束時(shí)，仍顯著低于CK處理(p<0.01)。水稻拔節(jié)期旱脅迫于2013年8月20日開(kāi)始(DAT 60，圖1)，脅迫歷時(shí)5 d后(DAT 65)，J-HD和J-LD處理的株高雖低于CK處理，但三者無(wú)顯著差異(p<0.05)，且均顯著高于T-HD和T-LD處理(p<0.01)；8月25日(DAT 65)J-HD和J-LD處理由旱脅迫轉(zhuǎn)入澇脅迫，脅迫歷時(shí)5 d，并于DAT73再次轉(zhuǎn)入旱脅迫; J-HD和J-LD處理經(jīng)過(guò)5 d持續(xù)的干旱脅迫后(DAT 78)，J-HD株高顯著低于其余4個(gè)處理(p<0.01)，且T-LD處理株高已經(jīng)趕超J-LD處理(74.8 cm vs 73.9 cm);拔節(jié)期脅迫結(jié)束后，隨即轉(zhuǎn)入淺水勤灌，DAT94測(cè)定的數(shù)據(jù)表明，T-HD(93.35 cm)和T-LD處理株高(92.35 cm)與CK(91.85 cm)無(wú)顯著差異，且顯著大于J-HD(83.05 cm)和J-LD處理(81.47 cm，p<0.01)。

圖2 不同時(shí)期旱澇交替脅迫作用下水稻株高的變化動(dòng)態(tài)Fig.2 Changes of rice plant heights under different irrigation regimes

2.2 植株生長(zhǎng)模型的擬合與分析

利用Logistic, Gompertz和von Bertalanffy 模型對(duì)不同生育階段旱澇交替脅迫的水稻株高進(jìn)行擬合。從表2和圖3可知，3種非直線生長(zhǎng)模型的擬合度均較好(R2均大于0.83)，說(shuō)明3種模型均能較好地?cái)M合各種灌溉條件下水稻生長(zhǎng)曲線，其中全生育期淺水勤灌和分蘗期旱澇交替脅迫條件下的R2達(dá)到0.90以上，而對(duì)于拔節(jié)期的交替重旱-輕澇處理(J-HD)，3種模型的R2均值卻不足0.83。另外，從表2可知，3種模型預(yù)測(cè)各灌溉條件水稻株高生長(zhǎng)過(guò)程線的最大株高、拐點(diǎn)高度、拐點(diǎn)時(shí)間均存在較大差異：Logistic 模型預(yù)測(cè)的拐點(diǎn)株高和拐點(diǎn)時(shí)間均顯著大于Gompertz 和von Bertalanffy模型，而Gompertz和von Bertalanffy模型預(yù)測(cè)的水稻最大株高均顯著大于Logistic 模型。結(jié)合該粳稻品種(Nanjing 44)的生長(zhǎng)特性，其最大株高不應(yīng)超過(guò)115 cm, Logistic模型預(yù)測(cè)值在此范圍內(nèi)，而Gompertz 和 von Bertalanffy模擬T-LD和T-HD處理下的最大株高已經(jīng)超過(guò)128 cm，不符合水稻實(shí)際生長(zhǎng)狀況。

表2 2種擬合非線性模型參數(shù)估計(jì)值和擬合度Tab.2 Fitting degree and parameter evaluation ofthree fitting nonlinear models

2.3 實(shí)際觀測(cè)值與擬合曲線估計(jì)值比較

在不同灌溉處理下，均用65%實(shí)測(cè)株高數(shù)據(jù)模擬出的模型預(yù)測(cè)株高，并用剩余35%的實(shí)測(cè)數(shù)據(jù)進(jìn)行驗(yàn)證。用1∶1 作圖法，比較模擬值與實(shí)測(cè)值在標(biāo)準(zhǔn)對(duì)稱線兩側(cè)的分布，即分別以實(shí)測(cè)值為橫坐標(biāo)，以模擬值為縱坐標(biāo)，作等間隔相關(guān)圖，同時(shí)計(jì)算二者的修正效率系數(shù)E1、修正一致性指標(biāo)d1和平均絕對(duì)誤差MAE。從表3和圖4可知，實(shí)測(cè)值與模擬值所建立的回歸方程的決定系數(shù)達(dá)到0.85以上，且回歸方程統(tǒng)計(jì)檢驗(yàn)均達(dá)到極顯著水平(p<0.01)；經(jīng)分析CK，T-LD和T-HD處理下的實(shí)測(cè)株高值和模擬株高值，可知 Logistic 模型條件下的E1、d1均大于Gompertz 和von Bertalanffy，而MAE小于Gompertz和von Bertalanffy模型，表明在CK，T-LD與T-HD灌溉方式下，用logistic模型模擬的株高生長(zhǎng)狀況與實(shí)際符合程度好，一致性優(yōu)越；經(jīng)分析J-LD和J-LD的實(shí)測(cè)株高和模擬株高在3種模型下E1、d1和MAE無(wú)差異，表明在J-LD和J-LD灌溉方式下3種模型模擬的生長(zhǎng)狀況無(wú)差異，但整體來(lái)看，各模型在模擬CK，T-LD與T-HD條件下株高生長(zhǎng)狀況比J-LD和J-HD處理的效果好(CK，T-LD與T-HD處理在各模型下的E1、d1均顯著大于J-LD和J-HD，而MAE則小于J-LD和J-HD處理)。

2.4 不同灌溉方式下株高增長(zhǎng)速率隨時(shí)間變化狀況

經(jīng)實(shí)際觀測(cè)值與擬合曲線估計(jì)值比較可知，上述3種非曲線模型中，Logistic 模型模擬5種灌溉方式下株高生長(zhǎng)狀態(tài)最好。通過(guò)對(duì)Logistic 模型求一階導(dǎo)數(shù)，得出株高增長(zhǎng)速率隨移栽后天數(shù)變化的方程為：

(4)

式中：v(t)為株高增長(zhǎng)速率隨時(shí)間的變化方程；A為水稻最大株高；K為相對(duì)生長(zhǎng)速率參數(shù)；B為常規(guī)參數(shù)。

圖3 不同灌溉方式下水稻株高生長(zhǎng)模型擬合Fig.3 Fitting nonliear models of rice height under different irrigation regimes

從圖5可知，水稻秧苗移入大田后，由于苗期的旱澇交替脅迫，在DAT 40之前，CK，J-LD和J-HD的增長(zhǎng)速率大于T-LD和T-HD處理；水稻進(jìn)入拔節(jié)期以后，J-LD和J-HD由于受到旱澇交替脅迫，其株高增長(zhǎng)速率明顯低于CK，T-LD和T-HD處理，甚至在DAT 60以后，T-LD和T-HD處理的株高增長(zhǎng)速率反超CK處理，表現(xiàn)出較強(qiáng)的補(bǔ)償效應(yīng)。

3 討 論

生長(zhǎng)曲線法是趨勢(shì)外推法的一種重要方法，廣泛應(yīng)用于描述及預(yù)測(cè)生物個(gè)體的生長(zhǎng)發(fā)育以及經(jīng)濟(jì)特性等領(lǐng)域中，其常用的數(shù)學(xué)模型有Logistics、Gompertz 和 von Bertaffany 等模型，其中Logistics和Gompertz 模型均具有飽和增長(zhǎng)特征，在預(yù)測(cè)植物生長(zhǎng)發(fā)育進(jìn)程等方面應(yīng)用最多[ 17-21]。前人研究表明Logistic生長(zhǎng)曲線模型能較好地描述植物生長(zhǎng)的模型[27,28]，而馬秋月等的研究發(fā)現(xiàn)von Bertalanffy模型擬合植物生長(zhǎng)的效果比較好，并認(rèn)為von Bertalanffy模型的拐點(diǎn)可變，而Logistic與Gompertz的拐點(diǎn)固定，是造成Logistic與Gompertz擬合效果差于von Bertalanffy的主要原因[19]。本研究運(yùn)用Logistics、Gompertz和von Bertaffany 模型模擬了不同灌溉方式下的株高變化，結(jié)果表明，不同灌溉處理株高增長(zhǎng)均符合漸增期-快增期-緩增期3個(gè)階段的“S”形曲線生長(zhǎng)的一般規(guī)律，說(shuō)明株高作為生物積累量的重要指標(biāo)可以反映植物的生長(zhǎng)狀況。比較水稻株高的拐點(diǎn)時(shí)間、拐點(diǎn)株高和最大株高的模型預(yù)測(cè)值與實(shí)測(cè)值的接近程度、綜合預(yù)測(cè)值和實(shí)測(cè)值的線性相關(guān)性、修正效率系數(shù)E1、修正一致性指標(biāo)d1以及平均絕對(duì)誤差MAE，可知Logistics 模型能更好地模擬不同灌溉方式下水稻株高變化動(dòng)態(tài)，這與Wu等[29]的研究結(jié)果相一致，而不同于馬秋月等的研究成果[19]。

圖4 不同灌溉方式下水稻株高模擬值與實(shí)測(cè)值的比較Fig.4 Comparison of measured rice height and estimated rice height

表3 3種非線性模型實(shí)測(cè)值和模擬值相關(guān)性比較Tab.3 Correlation between simulated values and observed values under three different nonlinear growth models

圖5 不同灌溉方式下株高增長(zhǎng)速率圖Fig.5 Growth rates of rice height under different irrigation regimes

根據(jù)Logistics模型曲線擬合不同灌溉方式下的株高生長(zhǎng)的數(shù)據(jù)(表2、表3、圖3和圖4 )，可知Logistic模型曲線模擬常規(guī)灌溉(CK)以及分蘗期旱澇交替脅迫的灌溉處理(T-LD和T-HD處理)的效果優(yōu)于拔節(jié)期旱澇交替脅迫的處理(J-LD和J-HD處理)。這可能解釋為：分蘗期旱澇交替脅迫下的株高可以通過(guò)后期補(bǔ)償效應(yīng)來(lái)保持 Logistic 模型推導(dǎo)過(guò)程[17]相對(duì)生長(zhǎng)率K和常規(guī)參數(shù)B在生育期前后保持協(xié)調(diào)的狀況有關(guān)；而對(duì)于水稻拔節(jié)期旱澇脅迫，由于株高生長(zhǎng)受限且后期的補(bǔ)償效應(yīng)較弱，不能保證上述參數(shù)的前后協(xié)調(diào)。

從表2和圖3可知，水稻最大株高以及拐點(diǎn)時(shí)間受旱澇交替脅迫實(shí)施的時(shí)期影響較大，并且不同時(shí)期的旱澇交替脅迫下拐點(diǎn)株高存在顯著差異，在分蘗期旱澇交替脅迫處理的拐點(diǎn)株高高于全生育期淺水勤灌的CK處理，而拔節(jié)期旱澇交替脅迫處理的值為最低，這說(shuō)明不同時(shí)期的旱澇交替脅迫顯著影響植物的生長(zhǎng)狀況。通過(guò)對(duì)Logistic 模型求一階導(dǎo)數(shù)，得出株高增長(zhǎng)速率隨水稻移栽天數(shù)的方程：分蘗期水稻經(jīng)歷旱澇交替脅迫后，能在生育后期產(chǎn)生一定的補(bǔ)償效應(yīng)，促進(jìn)株高的增長(zhǎng)；拔節(jié)期的旱澇交替脅迫后進(jìn)行復(fù)水，其株高反而小于整個(gè)生育期淺水勤灌的處理，這與前人的研究表明水稻分蘗期旱后復(fù)水的后效性大于拔節(jié)期的結(jié)論相一致[7,8,22,23]。近期我們的研究成果也表明，生育早期旱澇交替脅迫后水稻比常規(guī)灌溉較早形成通氣組織、較大的導(dǎo)管面積和通氣面積以及較大的光合速率是水稻出現(xiàn)明顯補(bǔ)償效應(yīng)的原因[24,25]。

Logistics模型相對(duì)于其他2個(gè)生長(zhǎng)模型可以較好地模擬出不同灌溉方式下植物的生長(zhǎng)狀況(表2、表3和圖3、圖4)；通過(guò)對(duì)Logistics模型一階求導(dǎo)得出的株高生長(zhǎng)速率隨時(shí)間變化的方程，也可以較好地模擬出分蘗期旱澇交替脅迫在生育后期對(duì)株高的補(bǔ)償效應(yīng)。然而通過(guò)圖5可知，生育后期模型模擬出T-LD和T-HD生長(zhǎng)速率與實(shí)際情況不符，因?yàn)樵谠摃r(shí)期，尤其是DAT 100以后，實(shí)際生產(chǎn)中株高的生長(zhǎng)已經(jīng)停止，而模型模擬出的株高生長(zhǎng)速率還有較大的值，這表明Logistics模型在模擬不同階段的旱澇交替脅迫對(duì)株高的影響時(shí)還存在一些缺陷，還需要進(jìn)一步改進(jìn)模型，使得該模型能夠更準(zhǔn)確地模擬不同旱澇交替脅迫下植株的生長(zhǎng)動(dòng)態(tài)，為估算水稻生長(zhǎng)后期的倒伏風(fēng)險(xiǎn)等情況提供理論及數(shù)據(jù)參考。

□

[1] 袁隆平.雜交水稻超高產(chǎn)育種[J].雜交水稻, 1997,12(6):1-3.

[2] 楊守仁,張步龍,王進(jìn)民, 等. 水稻理想株型育種的理論和方法初論[J]. 中國(guó)農(nóng)業(yè)科學(xué), 1984,(3):6-13.

[3] 趙軼鵬, 邵在勝, 王云霞, 等．大氣 CO2和O3濃度升高對(duì)汕優(yōu)63生長(zhǎng)動(dòng)態(tài)、物質(zhì)生產(chǎn)和氮素吸收的影響[J]. 生態(tài)學(xué)報(bào)，2015，35(24):8 128-8 138.

[4] 楊長(zhǎng)明, 楊林章，顏廷梅，等. 不同養(yǎng)分和水分管理模式對(duì)水稻抗倒伏能力的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào)， 2004，15(4):646-650.

[5] Shi G Y，Yang L X，Wang Y X，et al． Impact of elevated ozone concentration on yield formation of four Chinese rice cultivars under fully open-air field conditions[J]. Agriculture, Ecosystems and Environment, 2009,131(3 /4):178-184.

[6] Zhu C, Ziska L H, Sakai H, et al. Vulnerability of lodging risk to elevated CO2and increased soil temperature differs between rice cultivars[J]. Eur. J. Agron., 2013,46:20-24.

[7] 郝樹榮,郭相平,王文娟.旱后復(fù)水對(duì)水稻生長(zhǎng)的后效影響[J].農(nóng)業(yè)機(jī)械學(xué)報(bào) 2010，41(7):76- 79．

[8] 肖夢(mèng)華, 胡秀君, 褚琳琳. 水稻株高生長(zhǎng)對(duì)旱澇交替脅迫的動(dòng)態(tài)響應(yīng)研究[J]. 節(jié)水灌溉，2015,(9):15-22.

[9] 呂 露, 馮常萍, 崔遠(yuǎn)來(lái). 水稻不同節(jié)水灌溉模式的水分利用效率比較-以荊門、桂林為例[J]. 節(jié)水灌溉，2011,(11):15-18.

[10] 郭相平, 袁 靜, 郭 楓, 等. 水稻蓄水-控灌技術(shù)初探[J].農(nóng)業(yè)工程學(xué)報(bào)，2009,25(4):70-73.

[11] 郭相平，張展羽，殷國(guó)璽．稻田控制排水對(duì)減少氮磷損失的影響[J]. 上海交通大學(xué)學(xué)報(bào)：農(nóng)業(yè)科學(xué)版，2006,24(3):307-310.

[12] 郭相平, 袁 靜, 郭 楓, 等. 水稻蓄水-控灌技術(shù)初探[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2009,(4):70-73.

[13] 郭相平, 甄 博, 王振昌. 旱澇交替脅迫增強(qiáng)水稻抗倒伏性能[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2013,29(12): 130-135.

[14] 李小坤, 李云春, 魯劍巍, 等. 強(qiáng)降雨致洪澇災(zāi)害下不同因素對(duì)水 稻倒伏的影響[J].自然災(zāi)害學(xué)報(bào), 2012,21(6):99-103.

[15] 彭世彰, 張正良, 龐桂斌. 控制灌溉條件下寒區(qū)水稻莖稈抗倒伏力學(xué)評(píng)價(jià)及成因分析[J]．農(nóng)業(yè)工程學(xué)報(bào)，2009,25(1):6-10.

[16] Setter T L, Laureles E V, Mazaredo A M. Lodging reduces yield of rice by self-shading and reductions in canopy photosynthesis[J]. Field Crop Res, 1997,49:95-106.

[17] 莫惠棟. Logistic方程及其應(yīng)用[J].江蘇農(nóng)學(xué)院學(xué)報(bào), 1983,4(2):53-57.

[18] 錢永強(qiáng)，孫振元，韓 蕾，等. 基于Logistic 與Gompertz 非線性模型的野牛草克隆生長(zhǎng)模擬與分析[J].中國(guó)農(nóng)業(yè)科學(xué), 2011,44(11):2 252-2 259.

[19] 馬秋月, 陳贏男, 渠紀(jì)騰, 等. 簸箕柳種內(nèi)雜交F1群體株高生長(zhǎng)曲線的擬合[J]. 南京林業(yè)大學(xué)學(xué)報(bào), 2013,37(4):13-16.

[20] 姚克敏,鄒江石,買 苗,等.兩系法雜交水稻株高變化規(guī)律及其與氣象條件的關(guān)系[J]. 雜交水稻(增刊), 1999,14:49-51.

[21] 趙雨明，盧桂賓，賀義才．帥棗系品種枝條果實(shí)生長(zhǎng)發(fā)育動(dòng)態(tài)[J]．東北林業(yè)大學(xué)學(xué)報(bào)，2011，39(9):45-48．

[22] 郝樹榮,郭相平,王為木,等. 水稻分蘗期水分脅迫及復(fù)水對(duì)根系生長(zhǎng)的影響[J].干旱地區(qū)農(nóng)業(yè)研究, 2007,25(1):149-152.

[23] 郝樹榮, 郭相平, 張展羽, 等. 水稻根冠功能對(duì)水分脅迫及復(fù)水的補(bǔ)償響應(yīng)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2010, 41(5):52-55.

[24] 甄 博, 郭相平, 陸紅飛. 旱澇交替脅迫對(duì)水稻分蘗期根解剖結(jié)構(gòu)的影響[J].農(nóng)業(yè)工程學(xué)報(bào), 2015,31(9):107-113.

[25] 陸紅飛, 郭相平, 甄 博, 等. 旱澇交替脅迫條件下粳稻葉片光合特性[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2016,32(8):105-112.

[26] Legates DR, McCabe G J. Evaluating the use of 'goodness-of-fit'measures in hydrologic and hydroclimatic model validation[J]. Water Resour Res, 1999,35:233-241.

[27] Niklas K L．Plant allometry：the scaling of form and process[M]. Chicago：University of Chicago Press，1994．

[28] West G B，Brown J H，Enquist B J．A general model for ontogenetic growth[J]. Nature, 2001,413:628-631.

[29] Wu R L，Ma C X，Yang M C K，et al．Quantitative trait loci for growth trajectories in Populus[J]．Genetical Research, 2003,81(1):51-64．