水氮耦合對黑土稻作產量與氮素吸收利用的影響

秦子元 張忠學 孫 迪 宋 健 張作合 李鐵成

(1.東北農業大學水利與土木工程學院， 哈爾濱 150030；2.東北農業大學農業農村部農業水資源高效利用重點實驗室， 哈爾濱 150030；3.綏化學院農業與水利工程學院， 綏化 152061)

0 引言

水稻是我國主要糧食作物，2021年我國水稻總產量預計達到2.14億t，約占我國糧食總產量的1/3[1]。在耕地數量有限、人口持續增長的背景下，到2030年，我國水稻產量需提高20%，才能滿足國內糧食需求[2]。農民通常施用高量氮肥，以求最大限度地提升水稻產量，但持續的高量投入水氮資源也帶來了諸多負面影響[3]。一方面，由于過量施用氮肥、施肥方式不當以及田間水分管理不協調等因素，使得我國氮肥利用率僅為30%～35%，大量氮肥通過氨揮發、硝化與反硝化、地表徑流和淋溶等方式損失，導致土壤酸化、河流湖泊富營養化、地下水污染、溫室氣體排放加劇等環境問題[4]；另一方面，我國水資源短缺嚴重，而有約70%的農業用水用于稻田灌溉，傳統淹水灌溉模式不斷增加耗水，隨著水資源供應緊缺，稻田生態系統生產力將受到嚴重威脅[5]。目前，我國正在引導農業生產從單純追求高產向穩定生產且注重資源節約和環境友好的新型農業發展，因此，推進節水灌溉和優化施肥技術的研究，對我國農業可持續發展具有重要意義。

近30年，我國東北平原稻田平均施氮量為159 kg/hm2，遠高于99.5 kg/hm2的世界平均水平，減氮潛力巨大[6]。同時，隨著我國糧食生產重心北移，黑龍江省逐漸成為我國主要糧食生產基地和糧食戰略儲備基地，農業水資源消耗日益增加不僅加重了潛在旱災風險而且削弱了水體自凈能力[7]。近年來，隨著農業水資源供需矛盾加劇、環境污染事件頻發，以節水灌溉和優化施肥來實現水稻穩產和高產的理論和技術研究成為熱點[8]。目前，國內外專家學者關于水氮耦合對水稻產量和氮素吸收利用已有諸多研究[9-14]。以往研究大多關注于水稻產量、氮素吸收和水氮利用效率等方面，而對水稻干物質和氮素累積速率系統的研究較少，且適宜水氮耦合方式能否進一步提升水稻產量與植株氮素吸收轉運以及干物質和氮素累積速率的協同性還需探討。

本文基于大田試驗，研究黑土地區不同水氮耦合條件下的水稻產量、水氮利用效率、植株干物質累積量及累積速率、氮素累積量及累積速率、氮素在不同器官中的轉運規律以及產量相關性分析，旨在為建立區域性節水節肥、增產、減排的高效水氮耦合方式提供理論依據與技術參考。

1 材料與方法

1.1 試驗區概況

試驗于2020年5—10月在黑龍江省慶安國家灌溉試驗重點站(46°58′8″N，127°40′2″E)進行。試驗田位于慶安縣平安鎮，地處松嫩平原呼蘭河流域中下游，是典型的寒地黑土分布區，屬寒溫帶大陸性季風氣候，夏季溫熱多雨，冬季寒冷干燥。全年無霜期約128 d，多年平均氣溫2.9℃，多年平均蒸發量為1 213.4 mm，多年平均降水量559.8 mm，降水多集中在7、8月。供試土壤類型為黑土型水稻土，試驗前耕層土壤(0～20 cm)基本理化性質如下：pH值6.47，有機質質量比43.1 g/kg，全氮質量比1.69 g/kg，全磷質量比0.67 g/kg，全鉀質量比19.99 g/kg，堿解氮質量比159.21 mg/kg，速效磷質量比27.56 mg/kg，速效鉀質量比158.3 mg/kg。

1.2 試驗設計

采用灌溉模式和施氮量2因素全面試驗。設置常規淹灌(F)、淺濕灌溉(W)和控制灌溉(C)3種，常規淹灌在水稻返青期后保持30～50 mm水層；淺濕灌溉采用“前水不見后水”的灌溉方式，待田面呈濕潤狀態，再灌下次水；控制灌溉田面不再長時間建立水層，以根層土壤含水率及土壤表相確定灌水時間、灌水次數及灌水定額。不同灌溉模式各生育期水分管理見表1。全生育期施氮量設置4個水平：0、85、110、135 kg/hm2(N0、N1、N2、N3)。試驗共計12個處理，每個處理3次重復，共36個小區，隨機區組排列，各試驗小區面積為100 m2(10 m×10 m)，田埂高20 cm，寬25 cm，小區四側布置塑料板和水泥埂，以減少串流和側滲。各小區單獨灌排，進水管接裝小型水表。

表1 水稻各生育期水分管理Tab.1 Water management during each growth period of rice field

供試水稻品種為當地廣泛種植的“綏粳18”，水稻插秧每穴定3株，株行距為10 cm×30 cm。供試氮肥為尿素(含N 46%)，分基肥(45%)、蘗肥(20%)和穗肥(35%)施用；磷肥為過磷酸鈣(含P2O512%)，P2O5施入量為45 kg/hm2，全部作為基肥施入；鉀肥為硫酸鉀(含K2O 60%)，K2O施入量為80 kg/hm2，按基肥(50%)和8.5葉齡(50%)分施。2020年5月19日大田施基肥，5月20日水稻幼苗移栽大田，6月11日施蘗肥，7月18日施穗肥，9月25日收獲。稻田田間管理如病蟲草害防治等與當地大田生產相同。

1.3 測定項目與方法

1.3.1氣象數據測定

降水量和氣溫等氣象數據由自動氣象站(DZZ 2型，天津氣象儀器廠)記錄(圖1)。

1.3.2灌溉水量測定

水稻移栽大田后，每日于固定時刻用毫米刻度尺讀取田面水層深度(當田面有水層時)，或用土壤水分檢測儀測定大田土壤含水率(當田面無水層時)，當田面水層深度或土壤含水率達到各生育期水分管理下限時，則灌水至水分管理上限，記錄灌溉日期，并用小型水表計算各次灌溉水量。

1.3.3干物質量與氮素含量測定

各小區分別于分蘗期、拔節孕穗期、抽穗開花期、乳熟期和成熟期隨機選取水稻5穴，將整株水稻挖起，剪去根部后裝入自封袋運回實驗室清洗，泥土清洗干凈后，將植株分為葉、莖鞘和穗(抽穗后)分別裝入袋中，置于干燥箱中105℃殺青0.5 h，然后80℃干燥至質量恒定，冷卻至室溫用精度0.01 g電子天平測定各部分器官干物質量。稱量后用粉碎機將植株各部分器官粉碎后過80目(0.18 mm)網篩，經H2SO4-H2O2法消煮后，取待測液用連續流動分析儀(AutoAnalyzer-3型，德國 Bran+Luebbe公司)測定氮素含量。

1.3.4土壤無機氮含量測定

各小區分別于水稻施基肥前和收獲后，用土鉆按照五點取樣法取0～20 cm、20～40 cm、40～60 cm共3層土壤樣品，裝入自封袋混合均勻后，-20℃保存，各層土壤用100 mL環刀取土，密封后立即帶回實驗室，105℃干燥至質量恒定，計算土壤容重和含水率。土壤樣品解凍后，稱取5 g過2 mm篩孔新鮮土樣于250 mL廣口瓶中，并加1 mol/L KCl 100 mL充分浸提后，經定性濾紙過濾，待測液用連續流動分析儀(AutoAnalyzer-3型，德國 Bran+Luebbe公司)測定銨態氮和硝態氮含量。

1.3.5產量測定

水稻成熟期，各試驗小區隨機選取1 m2水稻人工收割，用脫粒機脫粒，稻谷攤開晾曬至含水率約14%時稱量測產。

1.4 計算方法

采用Logistic生長方程進行水稻干物質、氮素累積量的非線性回歸擬合[15]，其函數表達式為

y=a/(1+be-kt)

(1)

式中y——水稻地上部干物質或氮素累積量，kg/hm2

a——干物質或氮素累積量的潛在最大值，kg/hm2

b——與干物質或氮素累積量有關的回歸參數，b>0

k——干物質或氮素累積量的增長率，k<0

t——水稻移栽后時間，d

對Logistic函數分別求一階、二階導數，并求其特征值方程，可得Logistic曲線3個關鍵點，其橫坐標分別對應始盛期(T1，d)、高峰期(T2，d)和盛末期(T3，d)，以及最快增長速率(Vmax，kg/(hm2·d))，利用始盛期和盛末期可將Logistic曲線擬合的生長過程分為漸增期(0～T1)、快增期(T1～T3)和緩增期(T3～∞)。特征值方程分別為

T1=(lnb-1.317)/k

(2)

T2=lnb/k

(3)

T3=(lnb+1.317)/k

(4)

ΔT=T3-T1

(5)

Vmax=ak/4

(6)

式中 ΔT——快增期持續時間，d

作物耗水量(ET，mm)計算公式為

ET=P+I+G+ΔW-R-D

(7)

式中P——降水量，mm

I——單位面積灌水量，mm

G——地下水補給量，mm，由于試驗區地下水埋深較大，故G取零

ΔW——水稻移栽和收獲時0～60 cm土壤儲水變化量，mm

R——水稻生育期內排水量，mm

D——深層滲漏量，mm

灌溉水分利用效率(IWUE，kg/m3)計算公式為

IWUE=Y/I

(8)

式中Y——水稻產量，kg/hm2

水分生產效率(WUE，kg/m3)計算公式為

WUE=Y/ET

(9)

氮肥農學利用效率(AEN，kg/kg)計算公式為

AEN=(Ya-Yb)/Na

(10)

式中Ya——施氮處理水稻產量，kg/hm2

Yb——不施氮處理水稻產量，kg/hm2

Na——施氮量，kg/hm2

氮素籽粒生產效率(NGPE，kg/kg)計算公式為

NGPE=Y/Tn

(11)

式中Tn——水稻氮素總累積量，kg/hm2

氮肥偏生產力(PFP，kg/kg)計算公式為

PFP=Ya/Na

(12)

氮素收獲指數(NHI，%)計算公式為

NHI=Pn/Tn×100%

(13)

式中Pn——水稻穗部氮素累積量，kg/hm2

百千克籽粒吸氮量(NUG，kg)計算公式為

NUG=Tn/Y×100

(14)

氮素吸收量(Nu，kg/hm2)計算公式為

Nu=NDm

(15)

式中N——水稻植株氮素含量，%

Dm——水稻干物質量，kg/hm2

氮素轉運量(Nt，kg/hm2)計算公式為

Nt=Nh-Nm

(16)

式中Nh——抽穗期水稻莖鞘或葉氮素累積量，kg/hm2

Nm——成熟期水稻莖鞘或葉氮素滯留量，kg/hm2

氮素轉運率(Nte，%)計算公式為

Nte=Nt/Nh×100%

(17)

氮素轉運貢獻率(Ntce，%)計算公式為

Ntce=(Nlt+Nsst)/Npt×100%

(18)

式中Nlt——水稻葉氮素轉運量，kg/hm2

Nsst——水稻莖鞘氮素轉運量，kg/hm2

Npt——抽穗期至成熟期水稻穗氮素增加量，kg/hm2

土壤氮素表觀盈虧量(Sapl，kg/hm2)計算公式為

Sapl=Na+Nsi-Tn-Nsr

(19)

其中

Nsi(sr)=0.1dρω(N)

(20)

式中Nsi——土壤初始無機氮量，kg/hm2

Nsr——土壤殘留無機氮量，kg/hm2

d——土層厚度，cm

ρ——土壤容重，g/cm3

ω(N)——土壤中無機氮質量比，mg/kg

1.5 數據分析

采用Excel 2019和Origin 2017處理數據及作圖，利用SPSS 22.0進行方差分析，運用Origin 2017軟件進行Logistic生長函數擬合，模型有效性通過效率系數(EF)、一致性指數(AI)[16-17]和決定系數(R2)評價。

2 結果與分析

2.1 水氮耦合對水稻收獲后地上部干物質、氮素累積量的影響

不同水氮耦合處理下，水稻收獲后地上部干物質、氮素累積量見圖2(圖中不同小寫字母表示處理間差異顯著(P<0.05)，下同)。試驗結果表明，常規淹灌和淺濕灌溉模式下，水稻地上部各器官干物質累積量隨施氮量增加而增大；控制灌溉模式下，地上部各器官干物質累積量隨施氮量增加先增大后減小。除CN0處理莖鞘干物質累積量低于WN0，CN0處理穗干物質累積量低于WN0和FN0外，相同施氮量下，各處理不同器官干物質累積量，按灌溉模式由大到小依次為：控制灌溉、淺濕灌溉、常規淹灌。各處理穗干物質累積量均高于葉和莖鞘干物質累積量，葉、莖鞘和穗的干物質累積量分別占地上部干物質累積量的8.72%～11.24%、34.50%～38.17%、50.59%～56.15%，CN2處理穗部干物質累積量最大，但與CN3處理不存在顯著性差異(P>0.05)。3種灌溉模式下，水稻地上部不同器官氮素累積量均隨施氮量的增加而增大。相同施氮量下，常規淹灌和淺濕灌溉模式水稻地上部不同器官氮素累積量均低于控制灌溉；相同灌溉模式下，除FN0和FN1莖鞘氮素累積量不存在顯著性差異外(P>0.05)，其余各施氮處理莖鞘和穗的氮素累積量顯著高于不施氮處理(P<0.05)。從不同器官來看，葉、莖鞘和穗的干物質累積量分別占地上部氮素累積量的7.24%～9.26%、21.73%～26.67%、65.99%～70.55%，葉、莖鞘和穗的平均氮素累積量分別為12.03、33.43、98.11 kg/hm2，分別占地上部平均氮素累積量的8.38%、23.28%、68.34%，其中CN3處理葉、莖鞘和穗的氮素累積量最大，分別為19.71、49.62、145.81 kg/hm2。控制灌溉處理相同施氮量下的葉、莖鞘和穗的氮素累積量較常規淹灌提高了27.80%～43.42%、18.32%～24.97%、13.85%～24.25%，較淺濕灌溉提高了0.96%～13.18%、10.73%～12.86%、10.53%～12.61%。

2.2 水氮耦合對水稻地上部干物質、氮素積累速率的影響

不同水氮耦合處理下，水稻地上部干物質、氮素積累量的Logistic生長方程擬合及模型有效性評估如表2所示。不同處理水稻地上部干物質、氮素累積量Logistic生長模型一致性指數分別為0.998 4～0.999 8、0.996 1～0.999 2，效率系數分別為0.955 0～0.995 3、0.890 1～0.977 9，決定系數R2均不低于0.975 6，模型有效性評估表明，利用Logistic生長方程擬合干物質、氮素累積過程效果較好，模型的預測值與實測值吻合度較高。擬合結果顯示，不同水氮耦合條件下各處理水稻地上部干物質、氮素累積的理論最大值分別為11 073.51～17 738.81 kg/hm2、73.49～213.89 kg/hm2，其中CN2處理干物質累積理論最大值較其他處理高0.05%～63.87%，CN3處理氮素累積理論最大值較其他處理高10.15%～191.04%。

對Logistic生長方程求一階導數，得到水稻地上部干物質、氮素累積速率隨移栽后時間的變化曲線(圖3)。水稻地上部干物質、氮素累積速率的變化曲線為單峰曲線，隨移栽后時間的推進呈現先增加后減小的趨勢，不同水氮耦合方式下，氮素累積速率達到峰值時間早于干物質累積速率達到峰值的時間。相同灌溉模式下，水稻地上部干物質、氮素累積速率隨施氮量的增加而增大，但當施氮量為135 kg/hm2時，干物質、氮素累積速率明顯減緩。相同施氮量下，3種灌溉模式干物質、氮素累積速率由大到小依次為控制灌溉、淺濕灌溉、常規淹灌。

對Logistic生長方程求二階導數，得到水稻地上部干物質、氮素累積過程主要特征參數(表3)。不同水氮耦合處理地上部干物質、氮素累積速率平均快增期分別為34.37 d和30.72 d，地上部干物質、氮素累積達到最大增長速率分別為第71.44天和第60.05天，平均氮素累積速率達到峰值時間比平均干物質累積速率達到峰值時間提前11.39 d。3種灌溉模式下，干物質、氮素累積始盛期隨施氮量提升而提前，施氮處理平均干物質、氮素累積速率始盛期比不施氮處理提前3.08、6.20 d。在快增期內，施氮量為135 kg/hm2時，各灌水模式下的地上部干物質、氮素最大增長速率達到最大值，其中CN3處理地上部干物質、氮素最大增長速率為340.95、4.97 kg/(hm2·d)，與其他處理相比分別提升了2.55%～63.84%、12.56%～229.56%。

表2 水稻地上部干物質、氮素累積量的Logistic生長方程擬合及模型有效性Tab.2 Logistic growth equation fitting and model validity of dry matter and nitrogen accumulation in shoots of rice

表3 水稻地上部干物質、氮素累積過程主要特征值Tab.3 Main characteristic values of dry matter and nitrogen accumulation process in shoots of rice

2.3 水氮耦合對水稻抽穗后氮素轉運的影響

不同水氮耦合處理下，水稻抽穗期至成熟期的氮素轉運見表4。試驗結果表明，灌溉模式和施氮量對葉片和莖鞘的氮素轉運量、氮素轉運率，穗部氮素增加量和氮素轉運貢獻率均有極顯著影響(P<0.01)，灌溉模式和施氮量的交互效應對葉片和莖鞘的氮素轉運量、氮素轉運率有顯著影響(P<0.05)，對穗部氮素增加量和氮素轉運貢獻率有極顯著影響(P<0.01)。水稻抽穗期至成熟期葉片和莖鞘氮素轉運量變化趨勢一致，在0～110 kg/hm2施氮量范圍內，相同灌溉模式下，氮素轉運量隨施氮量的增加而增大，施氮量為135 kg/hm2時，氮素轉運量略有降低，但與施氮量為110 kg/hm2時的氮素轉運量差異不顯著(P>0.05)。水稻抽穗期至成熟期葉片氮素轉運率高于莖鞘氮素轉運率，相同灌溉模式下，葉片和莖鞘氮素轉運率以及氮素轉運貢獻率隨施氮量的增加而減小，穗部氮素增加量隨施氮量增加而增大。相同施氮量，3種灌溉模式對葉片和莖鞘氮素轉運量以及穗部氮素增加量的影響，均以控制灌溉處理最高，表明適宜水分脅迫可促進葉片和莖鞘氮素轉運以及籽粒氮素累積。

表4 不同水氮處理水稻抽穗期至成熟期的氮素轉運Tab.4 Nitrogen transfer from heading stage to maturity stage of rice treated with different water and nitrogen treatments

2.4 水氮耦合對水稻產量和水氮利用效率的影響

不同水氮耦合處理下，水稻產量和水分利用效率如圖4所示。試驗結果表明，相同灌溉模式，不施氮處理水稻產量顯著低于施氮處理(P<0.05)，常規淹灌、淺濕灌溉和控制灌溉增產幅度分別為24.45%～76.31%、31.03%～62.37%和66.30%～99.78%。常規淹灌和淺濕灌溉處理水稻產量隨施氮量的增加而提升，施氮量超過110 kg/hm2時，水稻增產幅度相對較小；控制灌溉處理水稻產量隨施氮量的增加先增加后小幅度下降，CN2處理水稻產量最高，為10 272.57 kg/hm2，但CN2和CN3處理水稻產量差異不顯著(P>0.05)。相同施氮量，除CN0外，控制灌溉模式各處理水稻產量均高于常規淹灌和淺濕灌溉模式水稻產量。施氮量為110、135 kg/hm2時，3種灌溉模式的IWUE和WUE差異不顯著(P>0.05)，常規淹灌和淺濕灌溉處理IWUE和WUE隨施氮量增加而增大，控制灌溉處理IWUE和WUE隨施氮量增加先增后降，CN2處理IWUE和WUE均最高，分別為3.14、2.04 kg/m3，與CN0相比提高了96.25%和70.57%。

不同水氮耦合處理對水稻氮素利用的影響如表5所示。試驗結果表明，灌溉模式、施氮量以及灌溉模式和施氮量的交互效應對氮肥農學利用率、氮肥偏生產力、百千克籽粒吸氮量、氮素籽粒生產效率和氮素收獲指數均有極顯著影響(P<0.01)。CN2處理氮肥農學利用率顯著高于其他處理(P<0.05)，當施氮量相同時，常規淹灌和淺濕灌溉模式各處理氮肥農學利用率均小于控制灌溉模式。相同灌溉模式下，氮肥偏生產力隨施氮量的增加而降低，施氮量相同時，3種灌溉模式的氮肥偏生產力由大到小依次為控制灌溉、淺濕灌溉、常規淹灌，控制灌溉下氮肥偏生產力比常規淹灌和淺濕灌溉分別提升了4.77%～28.03%和3.79%～14.41%。控制灌溉模式下水稻平均百千克籽粒吸氮量為1.84 kg，較淺濕灌溉和常規淹灌分別提高了6.84%和7.97%。施氮量為110、135 kg/hm2時，常規淹灌和淺濕灌溉模式下氮素籽粒生產效率無顯著差異(P>0.05)，控制灌溉模式下氮素籽粒生產效率差異顯著(P<0.05)。不同處理間氮素收獲指數介于65.99%～70.55%之間，表明水稻吸收的氮素大部分用于形成籽粒。

表5 水氮耦合對水稻氮素利用的影響Tab.5 Effect of water and nitrogen interaction on rice nitrogen utilization

2.5 水稻產量相關性分析

水稻產量與植株成熟期地上部干物質、氮素累積量關系如圖5所示。不同灌溉模式下，水稻產量和成熟期地上部干物質、氮素累積量之間存在顯著或極顯著的正相關(R2為0.937 1～0.998 4)，水稻成熟期干物質、氮素累積量均可作為衡量產量的評價依據或指標，稻作控制灌溉模式下，水稻產量隨地上部氮素累積量的增加先增大后減小，表明水稻植株對氮素吸收存在閾值，過多氮素累積不利于產量形成。水稻產量和水氮利用效率間相關分析見表6。水稻產量與IWUE、WUE、AEN、NUG之間呈極顯著正相關(P<0.01)，與NGPE之間呈極顯著負相關(P<0.01)，與其他指標之間相關性不顯著(P>0.05)；IWUE與WUE、AEN之間呈極顯著正相關(P<0.01)，與NUG之間呈顯著正相關(P<0.05)；WUE與AEN之間呈極顯著正相關(P<0.01)，與PFP之間呈顯著正相關(P<0.05)；AEN與PFP之間呈顯著正相關(P<0.05)；PFP與NUG之間呈極顯著負相關(P<0.01)，與NGPE之間呈極顯著正相關(P<0.01)；NUG與NGPE之間呈極顯著負相關(P<0.01)，與NHI之間呈顯著正相關(P<0.05)。

3 討論

水稻地上部干物質、氮素積累量能反映作物的生產能力[18]，是影響水稻產量的重要指標，水氮條件作為水稻生長發育的主要限制因子，適宜的水氮耦合模式不僅可以提高水稻產量和水氮利用效率，而且能夠影響水稻不同器官中干物質、氮素累積和分配，對提高水稻干物質、氮素累積量產生正效應[19]，但過量的水氮供應也會使作物“徒長”，不利于產量形成[20]。本研究表明，稻作控制灌溉模式下，施氮處理地上部不同器官干物質累積量高于常規淹灌和淺濕灌溉，控制灌溉地上部各器官干物質累積量隨施氮量增加先增大后減小，這是因為控制灌溉模式下葉片光合能力提升，有利于有機質累積，同時控制灌溉可以改變土壤水分狀況，改善根區通氣條件，為水稻根系生長提供適宜環境，不僅能夠促進根系發育，還能延緩后期根系衰老和退化，使得水稻在生長中后期能吸收更多養分[21]，這與孟翔燕等[22]的研究結果一致，但當施氮量為135 kg/hm2，控制灌溉模式下水稻出現“貪青”現象，過量氮肥施用使得水稻晚熟，生育期推遲，地上部各器官干物質累積量小幅度下降，影響水稻產量。控制灌溉和淺濕灌溉模式下水稻地上部氮素累積量高于常規淹灌，水稻不同器官氮素累積量與施氮量呈正相關，控制灌溉和淺濕灌溉處理相同施氮量下的穗氮素累積量較常規淹灌提高了13.85%～24.25%、2.27%～10.09%，說明節水灌溉模式有利于氮素吸收轉化并促進植株氮素向籽粒轉運，這與崔遠來等[23]的研究結果相似。

表6 水稻產量和水氮利用效率的相關系數Tab.6 Correlation coefficients between rice yield and water and nitrogen use efficiency

前人研究表明，水稻基因型、水分狀況、施肥水平和氣象因素等均會影響植株干物質和氮素累積過程，體現在Logistic模型上就是擬合參數的變化[24]。本研究結果表明，從水稻干物質、氮素累積速率的變化過程看，不同水氮耦合模式基本不影響水稻干物質、氮素累積速率的變化趨勢，各處理水稻干物質、氮素累積速率隨移栽后時間的變化仍為單峰曲線，且隨移栽后時間的推進先增后減，利用Logistic生長模型擬合效果較好，但不同水氮耦合模式下，水稻干物質和氮素累積對應的漸增期、快增期、緩增期、累積量和累積速率均有差異。葉廷紅等[25]研究表明，水稻氮素累積快增期明顯早于干物質累積快增期，干物質累積最大速率出現在拔節孕穗期，這與本文研究結果相似，這是由于水稻植株在分蘗期以生根、生葉和分蘗為主，需要大量的氮素來形成氮化物，當養分累積足夠多時進行大量物質生產，因此，水稻氮素累積速率到達高峰期后，干物質累積速率逐漸進入始盛期。有關研究認為，植株干物質和氮素累積速率高峰期和最大增長速率出現時間越早，快增期持續時間越短，越有利于干物質和氮素累積以及優質群體結構的形成[26]，這與本文結論相似又有所不同，本研究表明，3種灌溉模式下，干物質和氮素累積速率始盛期、高峰期和盛末期隨施氮量的增加而提前，但快增期持續時間并沒有隨施氮量的增加明顯縮短，甚至增加，相比于常規淹灌模式，淺濕灌溉和控制灌溉模式下水稻平均干物質累積速率快增期縮短1.59 d和1.42 d，平均氮素累積速率快增期縮短3.35 d和6.00 d，從干物質和氮素累積特征值看，節水灌溉處理具有快增期持續時間短和增長速率大的特征，更有利于水稻植株生長中期干物質和氮素累積。稻作控制灌溉模式，從抽穗開花期開始N3水平下植株干物質和氮素累積速率逐漸開始低于N2或N1水平，而成熟期水稻干物質累積量隨施氮量的增加先增大后減小，氮素累積量隨施氮量的增大而增大，說明水稻植株對氮素養分的空間競爭使其對氮素吸收存在閾值，植株氮素吸收過多并不會使水稻干物質量持續增加，反而會導致水稻減產和品質下降[27]。相同灌溉模式下，當施氮量為135 kg/hm2時，干物質和氮素累積速率明顯減緩，高氮條件下，水稻對氮素的需求與氮肥養分釋放速率不匹配，氮肥的高投入并沒有帶來高回報。此外，本文Logistic模型是以移栽后時間為自變量對水稻生長過程進行模擬分析，而年際間氣候狀況有所差異，單純通過移栽后時間表示植株不同生長過程可能會產生偏差，由于水稻生育期內不同階段植株干物質和氮素累積所需熱量固定，一些專家學者用有效積溫作為自變量對不同植株生長過程進行Logistic模擬[28]，并取得良好效果，但目前尚無統一方法進行植株Logistic生長過程模擬分析，這將是今后進一步的研究重點。

水稻產量是植株干物質積累、分配、運輸和轉化的結果，水氮利用效率則可以反映作物生產過程中的能量轉化效率，也能衡量作物水氮限制條件下的生長適宜程度[29]。已有研究表明，一般情況下隨著施氮量的增加，水稻產量、水分利用效率增加，而氮肥利用率降低[30]。本研究表明，施氮量在0～135 kg/hm2范圍內，常規淹灌和淺濕灌溉模式下水稻產量和水分利用效率隨著施氮量的增加而提升，控制灌溉模式下，當施氮量為135 kg/hm2時，水稻產量和水分利用效率呈下降趨勢，同時，控制灌溉模式下施氮處理水稻產量和水分利用效率顯著高于常規淹灌和淺濕灌溉模式，其原因在于控制灌溉模式可增強水稻根區通氣性，改善稻田土壤條件，加速土壤有機質礦化，促使水稻產量提高，同時，控制灌溉模式降低了灌水頻次和灌水量，田面不再長期建立水層，使得稻田滲漏量、蒸發量降低，田間蓄雨能力較常規淹灌和淺濕灌溉更高，徑流量減小，雨水利用率提升，進而提升水分利用效率[31]。相同灌溉模式下，不施氮或低氮處理的氮肥偏生產力和氮素籽粒生產效率較高，但水稻產量和水分利用效率卻顯著低于適當的高施肥處理。因此，協調水稻產量和水氮利用效率之間的矛盾，應在維持較高產量的基礎上，提升水氮利用效率，同時降低水稻營養體對氮素的過多吸收。

水稻植株中來自土壤的氮素與土壤肥力關系密切，衡量最佳水氮耦合處理，不僅要考慮水稻產量、氮素吸收和水氮利用效率，還需考慮土壤氮素盈虧狀況。研究區黑土層厚度一般在53～58 cm之間，黑土區水稻吸收肥料氮素的比例遠低于土壤氮素，本文將0～60 cm土體無機氮作為整體，研究土壤氮素表觀平衡，不同水氮處理下土壤氮素盈虧量見圖6。研究結果表明，相同施氮量下，不同灌溉模式對土壤氮素盈虧量的影響不顯著(P>0.05)，在施氮量0、85 kg/hm2條件下，土壤氮素虧損，在施氮量110、135 kg/hm2條件下，土壤氮素盈余，施氮量為85～110 kg/hm2時土壤氮素盈虧出現臨界點。在水稻生產中，土壤氮素高盈余意味著高的土壤氮素損失風險，這些氮素一部分殘留在土壤中供后季作物吸收利用，另一部分則通過淋溶、氨揮發和反硝化等途徑損失，增加了環境污染風險；土壤氮素長期虧損則會造成土壤退化，土壤氮素耗竭阻礙稻田系統可持續發展，因此，根據黑土區環境特征選擇適宜水氮耦合模式，將土壤氮素盈虧控制在合理范圍內極其重要。

4 結論

(1)水稻控制灌溉模式下，施氮處理地上部不同器官干物質累積量高于常規淹灌和淺濕灌溉，控制灌溉地上部各器官干物質累積量隨施氮量的增加先增大后減小，施氮量為135 kg/hm2時抑制水稻干物質累積。相同灌溉模式下，水稻地上部不同器官氮素累積量隨施氮量的增加而增大，相同施氮量，控制灌溉模式下葉、莖鞘和穗的氮素累積量較常規淹灌提高了27.80%～43.42%、18.32%～24.97%、13.85%～24.25%，較淺濕灌溉提高了0.96%～13.18%、10.73%～12.86%、10.53%～12.61%。

(2)水稻地上部干物質、氮素累積速率隨移栽后時間的推進先增加后減小，平均氮素累積速率達到峰值時間比平均干物質累積速率達到峰值時間提前11.39 d。相同灌溉模式下，干物質、氮素累積始盛期隨施氮量增加而提前，干物質、氮素累積速率隨施氮量的增加而增大，施氮量為135 kg/hm2時，干物質、氮素累積速率明顯減緩。相同施氮量下，控制灌溉模式干物質、氮素累積速率高于淺濕灌溉和常規淹灌。

(3)常規淹灌和淺濕灌溉模式水稻產量和水分利用效率隨施氮量的增加而增大，控制灌溉模式水稻產量和水分利用效率隨施氮量的增加先增大后減小，相同施氮量，控制灌溉模式水稻產量(CN0處理除外)和水分利用效率高于常規淹灌和淺濕灌溉模式，其中CN2處理水稻產量、IWUE和WUE均最大，分別為10 272.57 kg/hm2、3.14 kg/m3和2.04 kg/m3。控制灌溉模式可顯著提升水稻氮肥農學利用效率和氮肥偏生產力。

(4)相同灌溉模式下，水稻抽穗期至成熟期葉片和莖鞘氮素轉運量隨施氮量的增加先增大后減小，葉片氮素轉運率高于莖鞘氮素轉運率，葉片和莖鞘氮素轉運率以及氮素轉運貢獻率隨施氮量的增加而減小。相同施氮量，控制灌溉模式葉片和莖鞘氮素轉運量以及穗部氮素增加量高于常規淹灌和控制灌溉。

(5)植株干物質、氮素累積量能夠作為不同灌溉模式下衡量水稻產量的指示指標。水稻產量與灌溉水分利用效率、水分生產效率、氮肥農學利用效率、百千克籽粒吸氮量之間呈極顯著正相關(P<0.01)，與氮素籽粒生產效率之間呈極顯著負相關(P<0.01)。綜合考慮水稻產量、氮素吸收、水氮利用效率和土壤氮素盈虧等方面，CN2處理為最佳水氮耦合模式。