寧夏引黃灌區(qū)不同灌水處理下春小麥光響應(yīng)曲線模型研究*

金建新，李株丹，黃建成，何進(jìn)勤，何振嘉

(1. 寧夏農(nóng)林科學(xué)院農(nóng)業(yè)資源與環(huán)境研究所，銀川市，750002； 2. 寧夏大學(xué)土木與水利工程學(xué)院，銀川市，750021； 3. 陜西省土地工程建設(shè)集團(tuán)有限責(zé)任公司，西安市，710075)

0 引言

作物光合強(qiáng)度的大小決定其生產(chǎn)有機(jī)物能力的大小，凈光合速率是反映作物利用CO2、水等原料生產(chǎn)有機(jī)物質(zhì)的重要指標(biāo)，為研究不同有效輻射下作物光合作用的定量變化，建立了許多光響應(yīng)模型，包括直角雙曲線模型、非直角雙曲線模型、指數(shù)模型、葉子飄模型等[1]，并用在各種作物上，通過不斷改進(jìn)，取得了較為可靠的模擬效果，在不同施肥水平和灌水量條件下擬合水稻、小麥、玉米等作物的光響應(yīng)過程[2-4]，各模型精度雖有差異，但總體來說，能較好反應(yīng)不同有效凈輻射下各作物光合作用變化趨勢(shì)，其中葉子飄模型對(duì)栓條櫟和刺槐光響應(yīng)過程的模擬也具有較好的效果[5]。不同模型對(duì)不同條件、不同作物的模擬效果不同，因此，部分學(xué)者對(duì)特定條件下大豆、棉花、碧桃、豇豆、向日葵、毛竹、咖啡等進(jìn)行了光響應(yīng)曲線模型的比較研究[6-10]，提出各作物適宜的光響應(yīng)曲線模型，用于指導(dǎo)實(shí)際生產(chǎn)，以取得最佳經(jīng)濟(jì)效益。方保婷等[11]探討了施氮量對(duì)小麥灌漿后期光響應(yīng)的影響，通過改善不同光合有效輻射強(qiáng)度下旗葉光合速率，提出最佳施氮量。張立志[12]利用冬小麥光響應(yīng)曲線模型計(jì)算和分析了光合特征參數(shù)，結(jié)合產(chǎn)量和土壤鹽分變化特征提出了較為合理的微咸水灌水方案及灌溉制度。黃明等[13]分析了水分脅迫對(duì)小麥不同抗旱品種光響應(yīng)、產(chǎn)量影響，為不同生態(tài)條件下小麥品種的選擇提供參考。小麥光響應(yīng)結(jié)合其他參數(shù)，可進(jìn)一步優(yōu)化小麥田間管理、栽培措施等，但是針對(duì)寧夏引黃灌區(qū)滴灌條件，不同灌水量下春小麥適宜的光響應(yīng)模型報(bào)道較少。

因此，本文設(shè)置滴灌不同灌水量，分析評(píng)價(jià)直角雙曲線模型、非直角雙曲線模型、指數(shù)模型、葉子飄模型、動(dòng)力學(xué)模型5種模型的模擬效果和差異性，以確定不同灌溉定額下最優(yōu)的春小麥光響應(yīng)模型，并對(duì)相應(yīng)的光響應(yīng)參數(shù)進(jìn)行計(jì)算，以期為寧夏引黃灌區(qū)春小麥水分的高效利用提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

試驗(yàn)于2019年3—7月在寧夏農(nóng)林科學(xué)院永寧試驗(yàn)基地進(jìn)行，該試驗(yàn)站位于永寧縣望洪鎮(zhèn)農(nóng)豐村，地處北緯38°22′，東經(jīng)106°21′，海拔1 120 m，試驗(yàn)區(qū)屬典型的中溫帶大陸性氣候，光熱資源豐富，多年平均氣溫為8.7 ℃，多年平均積溫為3 245.6 ℃，年日照時(shí)數(shù)為2 866.7 h，無霜期平均167 d，多年平均降水量為201.4 mm，多年平均蒸發(fā)量為1 470.1 mm，試驗(yàn)?zāi)攴荽盒←溕趦?nèi)降雨量見表1。試驗(yàn)區(qū)土壤為粘壤土，耕層平均田間持水量為31.24%(體積百分比)，飽和含水量為38.75%，凋萎系數(shù)為12.78%，土壤平均容重為1.42 g/cm3。堿解氮75 mg/kg，速效磷24.7 mg/kg，速效鉀158 mg/kg，有機(jī)質(zhì)含量12.8 g/kg。

表1 春小麥生育期內(nèi)降雨量Tab. 1 Rainfall during the growth period of spring wheat

1.2 試驗(yàn)設(shè)計(jì)

春小麥品種為“寧春55號(hào)”，于2019年3月2日播種，播種方式為條播，播種量為375 kg/hm2，行距為10 cm，小麥于7月10日收獲，為保障出苗，試驗(yàn)地冬灌量為 1 350 m3/hm2。本試驗(yàn)為灌水量單因子試驗(yàn)，灌水方式為滴灌，滴頭采用內(nèi)鑲貼片式滴頭，滴頭間距為30 cm，滴灌帶鋪設(shè)間距為50 cm。根據(jù)滴灌小麥研究結(jié)論[14]，試驗(yàn)設(shè)置5個(gè)不同灌溉定額梯度，分別為 1 350 m3/hm2、1 950 m3/hm2、2 400 m3/hm2、2 850 m3/hm2、3 150 m3/hm2，在兩葉一心、分蘗期、拔節(jié)期、孕穗期、灌漿期和乳熟期6個(gè)時(shí)期進(jìn)行灌溉，如表2所示。試驗(yàn)共計(jì)5個(gè)處理，每個(gè)處理3個(gè)重復(fù)，共計(jì)15個(gè)小區(qū)，小區(qū)長(zhǎng)10 m，寬6 m，用水表嚴(yán)格控制灌水量，總施肥量N、P2O5、K2O分別為240 kg/hm2、135 kg/hm2、60 kg/hm2，全部利用滴灌隨水追入，其他田間操作均一致，定期鋤草和病蟲害防治。

表2 試驗(yàn)處理表Tab. 2 Test treatmentTable m3/hm2

1.3 測(cè)定指標(biāo)與方法

在孕穗期和灌漿期灌水后，選擇晴朗的天氣，每個(gè)小區(qū)隨機(jī)選擇3株，觀測(cè)時(shí)間為早上9：00至11： 00，選擇長(zhǎng)勢(shì)正常、完全展開的葉片進(jìn)行觀測(cè)。采用CIRAS-3型便攜式光合儀自帶紅藍(lán)光源測(cè)定不同光合有效輻射梯度下的小麥葉片凈光合速率，開機(jī)后預(yù)熱20 min，設(shè)置氣源CO2濃度固定為400 μmol/mol，光照范圍設(shè)置11個(gè)光強(qiáng)梯度，依次為0、50、100、300、500、800、1 000、1 200、1 500、1 800、2 000 μmol/(m2·s)。儀器自動(dòng)記錄胞間CO2濃度、氣孔導(dǎo)度、葉片水蒸氣壓虧缺、凈光合速率、蒸騰速率和光合水分利用效率等指標(biāo)。

各小區(qū)隨機(jī)選擇5株小麥，利用LI-3000C葉面積儀測(cè)定所有植株綠色葉片的面積，取平均值后計(jì)算小麥生育期葉面積指數(shù)。利用TRIME-PICO管式TDR系統(tǒng)測(cè)定土壤體積含水率，每個(gè)小區(qū)在滴灌帶處埋設(shè)1根TDR管，深度為1 m，在灌水前后按照20 cm一層分層測(cè)定土壤含水率，降雨前后加測(cè)，首次測(cè)量時(shí)，利用烘干稱重法校核TDR測(cè)定的土壤含水率。

1.4 光響應(yīng)曲線模型

本文采用5種光響應(yīng)曲線模型進(jìn)行模擬計(jì)算[15]，利用DPSv18.10中非線性回歸分析對(duì)光響應(yīng)數(shù)據(jù)進(jìn)行擬合。

1) 直角雙曲線模型

(1)

式中：I——有效輻射，μmol/(m2·s)；

Pn——凈光合速率，μmol/(m2·s)；

α——表觀量子效率或者初始量子效率，表示光響應(yīng)曲線的初始斜率；

Rd——暗呼吸速率，μmol/(m2·s)，表示有效輻射為0時(shí)的呼吸速率；

Pmax——最大凈光合速率，μmol/(m2·s)，表示光響應(yīng)曲線的極值。

2) 非直角雙曲線模型

(2)

式中：θ——區(qū)角。

3) 指數(shù)模型

指數(shù)模型自1991年提出以來，發(fā)展出的形式較多，但均為沒有極值的函數(shù)，本文采用BASSMAN和ZWIER提出的指數(shù)模型[12]，其表達(dá)式為

Pn(I)=Pmax(1-e-αI/Pmax)-Rd

(3)

4) 葉子飄模型

(4)

式中：β——光抑制項(xiàng)修正系數(shù)；

γ——光飽和項(xiàng)。

5) 動(dòng)力學(xué)模型

(5)

式中：Km——當(dāng)Pn是最大凈光合速率一半時(shí)的光合有效輻射值；

Ic——作物在光補(bǔ)償點(diǎn)時(shí)的光合有效輻射值。

1.5 評(píng)價(jià)指標(biāo)

本研究采用決定系數(shù)R2、一致性指數(shù)IA、偏差BIAS、均方根誤差RMSE、平均絕對(duì)誤差MAE等指標(biāo)對(duì)不同模型擬合得到的春小麥光響應(yīng)曲線進(jìn)行評(píng)價(jià)，對(duì)各模型模擬的精確性進(jìn)行評(píng)估，表達(dá)式如式(6)～式(10)所示。

(6)

(7)

(8)

(9)

(10)

式中：Oi——實(shí)測(cè)值；

Si——模擬值；

R2值、IA越接近1，RMSE、MAE、BIAS越小說明模擬值誤差越小，模擬效果越好。

1.6 數(shù)據(jù)處理方法

試驗(yàn)數(shù)據(jù)利用Excel2009進(jìn)行數(shù)據(jù)處理和繪圖，利用origin2019進(jìn)行模型回歸和參數(shù)率定，利用DPSv18.10進(jìn)行差分析。

2 結(jié)果與分析

2.1 土壤儲(chǔ)水量和小麥葉面積指數(shù)變化

圖1為各處理小麥全生育期0～30 cm內(nèi)土壤儲(chǔ)水量變化情況?？梢钥闯?，在5月19日孕穗之前，同一時(shí)期各處理土壤儲(chǔ)水量變化較大，特別是灌水前處理5、處理4和處理3土壤儲(chǔ)水量大于處理2和處理1，如在4月27日土壤儲(chǔ)水量最大值的處理5較最小的處理2大33.47%。在小麥孕穗后，除了處理1外，其他處理土壤儲(chǔ)水量差異逐漸縮小，灌水前各處理土壤儲(chǔ)水量最大值較最小值大20.71%～25.44%。

圖1 土壤儲(chǔ)水量變化Fig. 1 Changes in soil water storage

圖2為孕穗期小麥葉面積指數(shù)LAI?？梢钥闯?，隨著灌水量的增加，小麥葉面積指數(shù)也增加，各處理間表現(xiàn)為極顯著性差異(P<0.01)。尤其當(dāng)灌溉定額超過3 000 m3/hm2時(shí)，小麥葉面積指數(shù)迅速增加，處理4較處理3大13.64%，而處理5較處理4則大32.65%，說明當(dāng)增加灌水量時(shí)，能促進(jìn)小麥葉片細(xì)胞分裂和伸長(zhǎng)，增大葉面積，為小麥光合作用提高和后期光合同化物積累奠定基礎(chǔ)。

圖2 孕穗期小麥LAIFig. 2 Wheat LAI in booting stage

2.2 不同灌溉定額下春小麥的光響應(yīng)曲線變化特征

由光響應(yīng)曲線(圖3)可知，不同滴灌灌水量處理下春小麥凈光合速率均隨著光合有效輻射的增加而增加，當(dāng)光合有效輻射超過一定值后，其凈光合速率趨于平緩，甚至表現(xiàn)為下降的趨勢(shì)，如處理1在孕穗期有效輻射為2 000 μmol/(m2·s)時(shí)春小麥凈光合速率較1 800 μmol/(m2·s)時(shí)低5.97%，其他處理均有類似的規(guī)律。

(a) 孕穗期

(b) 灌漿期 圖3 不同灌溉定額下灌漿前后期春小麥光響應(yīng)曲線Fig. 3 Light response curves of spring wheat before and after filling under different irrigation quotas

在同一有效輻射下，各處理凈光合速率表現(xiàn)為隨著灌水量的增加而增加，這種現(xiàn)象在灌漿期當(dāng)有效輻射超過800 μmol/(m2·s)時(shí)更為明顯，但在孕穗期當(dāng)有效輻射超過300 μmol/(m2·s)時(shí)就已經(jīng)凸顯，說明孕穗期春小麥耗水比較大，對(duì)水分更為敏感。在孕穗期當(dāng)?shù)喂喙嗨砍^2 400 m3/hm2時(shí)，隨著光合有效輻射的增加，春小麥凈光合速率迅速增加，說明存在一個(gè)灌水閾值，當(dāng)超過該值后，春小麥能迅速提高光合速率，增加同化物積累量。

2.3 5種模型對(duì)春小麥灌漿期光響應(yīng)的模擬

5種模型的春小麥灌漿期光響應(yīng)曲線如圖4所示。由圖4(a)可以看出，處理1當(dāng)光合有效輻射PAR大于100 μmol/(m2·s)時(shí)，各模型模擬結(jié)果表現(xiàn)出較大的差異，在500 μmol/(m2·s)范圍內(nèi)，動(dòng)力學(xué)模型和直角雙曲線模型模擬效果較好。當(dāng)PAR大于1 200 μmol/(m2·s)時(shí)，除了葉子飄模型外，其余模型模擬Pn均較實(shí)測(cè)值小。圖4(b)是處理2的春小麥光響應(yīng)模擬曲線，可以看出非直角雙曲線模型模擬值偏低，當(dāng)PAR較小時(shí)，葉子飄模型模擬效果最好，但是當(dāng)PAR大于1 000 μmol/(m2·s)時(shí)，葉子飄模型模擬值偏大，其余指數(shù)模型、直角雙曲線模型、動(dòng)力學(xué)模型整體模擬效果較好。圖4(c)中5種模型模擬結(jié)果較一致，和實(shí)測(cè)值也較為接近，特別是當(dāng)PAR小于1 000 μmol/(m2·s)時(shí)模擬效果較好。圖4(d)中對(duì)處理4模擬結(jié)果可以看出，當(dāng)PAR小于800 μmol/(m2·s)時(shí)，5種模型模擬結(jié)果和實(shí)測(cè)值均較接近，PAR在800～1 500 μmol/(m2·s)之間時(shí)，模擬值均較實(shí)測(cè)值低。圖4(e)中指數(shù)模型對(duì)處理5的模擬效果較好，葉子飄模型和非直角雙曲線模型模擬值較大，而動(dòng)力學(xué)模型和直角雙曲線模型模擬值較低。整體來看，各模型對(duì)處理3的模擬結(jié)果較為一致，說明灌水量過大或者過小，不同模型模擬結(jié)果會(huì)產(chǎn)生不同程度的誤差。當(dāng)超過光飽和點(diǎn)時(shí)，葉子飄模型能較好地模擬春小麥實(shí)際光響應(yīng)過程，即隨著PAR的增加Pn減小，其余模型沒有明顯的下降趨勢(shì)，不同滴灌定額灌水處理下直角雙曲線模型和動(dòng)力學(xué)模型模擬結(jié)果均較為接近。

(a) 處理1

(b) 處理2

(c) 處理3

(d) 處理4

(e) 處理5 圖4 5種模型的春小麥灌漿期凈光合速率光響應(yīng)曲線Fig. 4 Light response curves of net photosynthetic rate of spring wheat in the early stage of grain filling

利用不同光響應(yīng)曲線模型對(duì)各處理小麥孕穗期的葉片光響應(yīng)曲線進(jìn)行擬合，擬合所得光響應(yīng)數(shù)據(jù)見表3。由表3分析可知，各模型對(duì)不同滴灌灌水量下春小麥孕穗期光響應(yīng)過程模擬的決定系數(shù)較大，均達(dá)到了0.98以上，但是僅根據(jù)決定系數(shù)并不能判斷模型模擬值的可靠程度，還需要借助其他指標(biāo)綜合判斷。最大凈光合速率是表示春小麥光合同化能力的指標(biāo)，是最大光合速率和呼吸速率之差，處理1、處理2、處理3、處理4和處理5，利用光合儀實(shí)測(cè)得到各處理春小麥最大凈光合速率分別為27.97 μmol/(m2·s)、27.87 μmol/(m2·s)、30.57 μmol/(m2·s)、30.97 μmol/(m2·s)和31.07 μmol/(m2·s)，葉子飄模型模擬結(jié)果與實(shí)測(cè)值較為接近，非直角雙曲線模型對(duì)處理1和處理2的模擬結(jié)果較好，其他模型對(duì)不同滴灌量處理下春小麥最大凈光合速率的模擬值均偏大，不能反映實(shí)際同化能力。

作物暗呼吸速率是在無光條件下，作物的呼吸速率，因此，其值和作物生理活性密切相關(guān)，暗呼吸速率Rd非直角雙曲線模型和葉子飄模型模擬效果較好，其他3種模型規(guī)律性較差。

表觀量子效率是反應(yīng)作物在光照強(qiáng)度比較弱時(shí)，作物吸收和利用太陽輻射能力的指標(biāo)，實(shí)測(cè)結(jié)果為0.022 3、0.024 8、0.024 1、0.029 8和0.030 9，隨著灌水量的增加而增加，表明春小麥隨著滴灌灌溉定額的增加，其吸收和利用弱光強(qiáng)的能力也提高，這可能是因?yàn)橥寥浪殖渥銜r(shí)小麥生理活性較高。模擬結(jié)果來看，直角雙曲線模型、指數(shù)模型、葉子飄模型和動(dòng)力學(xué)模型規(guī)律較好，但是除了動(dòng)力學(xué)模型外，其余模型模擬得到的表觀量子效率值均偏大。葉子飄等指出，表觀量子效率不能充分體現(xiàn)作物光能利用效率，而用光補(bǔ)償點(diǎn)更加科學(xué)，5種模型中只有非直角雙曲線模型對(duì)光補(bǔ)償點(diǎn)的模擬規(guī)律較好。

表3 5種光響應(yīng)模型擬合的主要光合參數(shù)Tab. 3 Photosynthetic parameters obtained by fitting five light response models to the light response data of spring wheat

2.4 5種模型對(duì)春小麥光響應(yīng)模擬的精度分析

利用RMSE、MAE、BIAS、IA和R2指標(biāo)計(jì)算各模型對(duì)不同滴灌定額處理下孕穗期春小麥光響應(yīng)過程進(jìn)行計(jì)算，評(píng)估各模型的模擬精度，結(jié)果見表4。

分析可知，不同模型對(duì)各處理春小麥光響應(yīng)曲線模擬的IA和R2均為0.99以上，因此，主要根據(jù)RMSE、MAE和BIAS這3個(gè)指標(biāo)進(jìn)行評(píng)價(jià)。

處理1和處理2中葉子飄模型對(duì)光響應(yīng)曲線模擬值的RMSE、MAE、BIAS這3個(gè)指標(biāo)最小，2個(gè)處理中該3指標(biāo)較其他模型分別小15.9%～33.4%、23.01%～44.5%和21.3%～36.8%。同時(shí)IA、R2指標(biāo)也最大，2個(gè)處理中葉子飄模型較其他模型大0.13%～0.3%、0.4%～1.2%，差異較小。說明葉子飄模型對(duì)處理1和處理2春小麥孕穗期的光響應(yīng)模擬效果最好。處理3和處理5中葉子飄模型、指數(shù)模型的模擬結(jié)果較好，如RMSE指標(biāo)，兩模型較其他3個(gè)模型均大2.0%～3.0%，兩模型之間相差僅為0.13%。處理4中葉子飄模型、指數(shù)模型和非直角雙曲線模型的模擬效果較好。葉子飄模型、指數(shù)模型和非直角雙曲線模型對(duì)處理3、處理4、處理5的模擬精度較高，各模型對(duì)灌水量較小處理的光響應(yīng)過程模擬效果較差，寧夏引黃灌區(qū)春小麥灌溉定額小于2 400 m3/hm2時(shí)，可以利用葉子飄模型來模擬春小麥孕穗期的光響應(yīng)過程。

表4 5種模型對(duì)光響應(yīng)曲線的模擬精度比較分析Tab. 4 Comparison and analysis of simulation accuracy of 5 models to light response curve

2.5 不同灌水處理下小麥產(chǎn)量

滴灌條件下隨著灌水量的增加穗長(zhǎng)、千粒重、穗粒數(shù)、穗數(shù)等產(chǎn)量構(gòu)成要素也增加，各處理表現(xiàn)為顯著性差異。如表5所示，穗長(zhǎng)和穗粒數(shù)均是處理4最大，處理1最小，最大值較最小值分別大25.4%和11%。千粒重和穗數(shù)均為處理5最大，處理1最小，最大值分別較最小值大13.2%和11.3%。產(chǎn)量為處理4最大，其次是處理3和處理5。整體來看，當(dāng)灌水量大于2 400 m3/hm2(處理3)時(shí)，各指標(biāo)均顯著增加，說明水量不足嚴(yán)重抑制了小麥產(chǎn)量及相關(guān)要素生長(zhǎng)發(fā)育。

表5 不同處理小麥產(chǎn)量及產(chǎn)量構(gòu)成要素Tab. 5 Wheat yield and yield components under different treatments

3 討論

作物生長(zhǎng)期土壤水分主要來源為灌溉、降雨和地下水補(bǔ)給，水分是作物生長(zhǎng)重要的原料，參與光合、同化、蒸騰等多個(gè)生理過程，和小麥生長(zhǎng)密切相關(guān)。小麥孕穗期之前棵間蒸發(fā)較大，消耗過多的土壤水分，造成灌水量較小的處理土壤儲(chǔ)水量迅速減小，隨著小麥葉面積的增加，棵間蒸發(fā)消耗減少，前期水分不足促進(jìn)小麥根系生長(zhǎng)，吸收更多的深層土壤水分，同時(shí)也使小麥生育期提前，造成各處理間土壤儲(chǔ)水量差異減小。在畦灌條件下過大的灌水量會(huì)導(dǎo)致群體內(nèi)通風(fēng)透光條件變差，降低小麥葉面積指數(shù)[16]，在灌水量適當(dāng)時(shí)，小麥葉面積指數(shù)和灌水具有顯著的相關(guān)性[17]。滴灌條件下，小麥葉面積指數(shù)隨著灌水量增大而增大[18]，這和本文研究結(jié)果一致，當(dāng)灌水量一致時(shí)，增加灌水次數(shù)可以提高小麥葉面積指數(shù)，在灌溉定額為360～405 mm時(shí)增幅尤為顯著[19]?？傊?，隨著灌水量的增加小麥葉面積指數(shù)也增大，因小麥為淺根系作物，較難利用深層土壤水分，所以在灌溉定額不變的條件下，增加灌水次數(shù)也可以顯著提高小麥葉面積指數(shù)。

光合作用是作物進(jìn)行能量轉(zhuǎn)換的重要過程，水不僅是光合作用的介質(zhì)，也是重要原料之一。隨著灌水量的增加，各光合有效輻射下小麥凈光合速率和灌水量成正比[20]，這和本研究結(jié)果一致，但是當(dāng)灌水量過大時(shí)，小麥凈光合速率反而有所降低，說明小麥在開花后受到水分脅迫或者土壤水分不足，尤其是小麥前期也受到水分虧缺時(shí)，小麥PSII反應(yīng)中心失去活性，損傷了光合機(jī)構(gòu)，光合作用受到抑制，造成光飽和點(diǎn)降低，光抑制現(xiàn)象提前[21]。5種模型對(duì)小麥光響應(yīng)曲線的擬合均具有一定的精度，但是當(dāng)小麥?zhǔn)艿剿置{迫時(shí)，各模型的模擬精度差異較大，可能是因?yàn)楣嗨蛔?，影響了春小麥孕穗期生理過程，抑制了其生長(zhǎng)發(fā)育和光合同化作用，造成模擬值偏差較大。直角雙曲線模型和非直角雙曲線模型、指數(shù)模型對(duì)高灌水處理下的春小麥光響應(yīng)曲線的模擬精度顯著高于低灌水處理[15]，在充分灌溉條件下，指數(shù)模型模擬結(jié)果較實(shí)測(cè)值偏大，非直角雙曲線模型較模擬值偏小，除了葉子飄模型外，其余模型不收斂，不能很好地模擬光飽和后光響應(yīng)曲線的走勢(shì)情況[19]，這和本文結(jié)果一致。

4 結(jié)論

1) 各處理土壤儲(chǔ)水量和孕穗期LAI均與灌水量呈正相關(guān)，處理5較處理1的土壤儲(chǔ)水量大4.11%～37.23%，在小麥孕穗后差異逐漸減小。處理5孕穗期平均LAI較其他處理大32.34%～50.91%，特別當(dāng)灌溉定額超過3 000 m3/hm2時(shí)，小麥葉面積指數(shù)迅速增加。

2) 隨著光合有效輻射的增加，各處理春小麥凈光合速率先增加后平緩或減小，該變化閾值介于1 200～1 800 μmol/(m2·s)之間。孕穗期當(dāng)灌水量超過2 400 m3/hm2時(shí)，春小麥凈光合速率迅速增加，說明在一定范圍內(nèi)適當(dāng)提高灌水量，可以促進(jìn)春小麥同化速率。

3) 各模型對(duì)各處理春小麥孕穗期光響應(yīng)過程模擬的R2均達(dá)到0.98以上，但是各光合參數(shù)不同模型差異較大，最大凈光合速率是葉子飄模型模擬效果較好，表觀量子效率是動(dòng)力學(xué)模型規(guī)律暗呼吸速率Rd非直角雙曲線模型和葉子飄模型模擬效果較好。

4) 葉子飄模型、指數(shù)模型和非直角雙曲線模型對(duì)處理3、處理4、處理5的模擬精度較高，各處理RMSE和MAE最大值為1.25 μmol/(m2·s)、0.98 μmol/(m2·s)，其他3個(gè)指標(biāo)差異不大。當(dāng)寧夏引黃灌區(qū)春小麥灌溉定額小于2 400 m3/hm2時(shí)，可以利用葉子飄模型來模擬春小麥孕穗期的光響應(yīng)過程。

5) 穗長(zhǎng)、穗粒數(shù)和產(chǎn)量為處理4最大，較其他處理分別大6.1%～25.4%、2.5%～11.1%、1.4%～14.8%，千粒重和穗數(shù)為處理5最大，較其他處理分別大0.2%～12.2%和2.3%～9.2%。整體來看，當(dāng)灌水量大于2 400 m3/hm2(處理3)時(shí)，春小麥產(chǎn)量及各指標(biāo)均顯著增加。