長期低投入的稻魚共作對土壤養(yǎng)分循環(huán)及有機(jī)碳庫穩(wěn)定性的影響

摘要： 【目的】研究長期低投入稻魚共作模式下土壤有機(jī)碳庫和養(yǎng)分化學(xué)計(jì)量比等的變化，深入解析稻魚模式的碳庫變化，提升稻魚模式生態(tài)功能。【方法】以湖南傳統(tǒng)稻魚共作代表縣辰溪縣的長期低投入稻魚共作系統(tǒng)為研究對象，分析共作15 年（RF15） 和50 年（RF50） 稻魚模式下土壤養(yǎng)分含量、理化性質(zhì)、養(yǎng)分儲量、有機(jī)碳組分、化學(xué)計(jì)量比和土壤碳庫指數(shù)的變化，探討關(guān)鍵影響因子?！窘Y(jié)果】較之水稻單作（RM），RF15 的有機(jī)碳、全氮、全磷、硝態(tài)氮、銨態(tài)氮、速效磷含量，有機(jī)碳儲量、氮儲量、磷儲量，易氧化有機(jī)碳含量及其在有機(jī)碳中的占比均顯著降低；C∶N、pH、顆粒有機(jī)碳含量、顆粒有機(jī)碳占比、礦質(zhì)結(jié)合態(tài)有機(jī)碳含量占比、碳庫活度和碳庫活度指數(shù)分別顯著增加47.98%、13.15%、35.47%、72.24%、31.68%、58.07% 和58.07%。RF50 的全氮含量、C∶P、N∶P、硝態(tài)氮含量分別顯著增加13.13%、33.08%、14.31% 和51.52%，全磷含量、容重、有機(jī)碳儲量、磷儲量、易氧化有機(jī)碳含量、礦質(zhì)結(jié)合態(tài)有機(jī)碳含量及其占比均顯著降低。RF15 和RF50 的碳庫管理指數(shù)整體保持穩(wěn)定，且均大于100。在長期稻魚模式下，與土壤有機(jī)碳庫穩(wěn)定指標(biāo)顯著相關(guān)的環(huán)境因子為全氮含量、C∶N、有機(jī)碳儲量?！窘Y(jié)論】低投入的稻魚模式系統(tǒng)雖然會導(dǎo)致養(yǎng)分流失，但自我維持時(shí)間的延長會使系統(tǒng)具備養(yǎng)分儲量自我恢復(fù)的功能。長期稻魚種養(yǎng)可以提升土壤質(zhì)量，進(jìn)而改善稻田生態(tài)環(huán)境。研究結(jié)果為優(yōu)化稻魚共作模式提供了理論依據(jù)，有助于提升該模式的生態(tài)功能和可持續(xù)性。

關(guān)鍵詞： 稻魚共作；養(yǎng)分含量；養(yǎng)分儲量；化學(xué)計(jì)量比；碳庫管理指數(shù)

中圖分類號： S181；S154.1 文獻(xiàn)標(biāo)志碼： A 文章編號： 1001-411X（2024）06-0856-09

稻魚共作模式是利用稻、魚互利共生的關(guān)系[1]，實(shí)現(xiàn)農(nóng)田提質(zhì)增效的一種綠色、高效的生態(tài)種養(yǎng)模式[2-3]，機(jī)理包括魚排泄物的增肥效應(yīng)和魚活動對土壤養(yǎng)分礦化的促進(jìn)作用。稻田水產(chǎn)動物和水稻互利共生，增加生物多樣性，提高土地、光、水等資源的利用效率，在鄉(xiāng)村振興的大背景下得到大面積的示范推廣。截至2019 年，我國使用稻魚共生模式的稻田面積已達(dá)到200 萬hm2[4]，總收益穩(wěn)定提升[5]。

大量研究表明，稻田復(fù)合種養(yǎng)模式可以不同程度地提高稻田土壤中有機(jī)質(zhì)以及氮、磷、鉀等養(yǎng)分含量，有效提高土壤肥力，使稻田在少施或不施化肥的情況下也能為水稻生長提供足量營養(yǎng)物質(zhì)[6-10]。水生動物在稻田里的活動有利于稻田土壤通氣，排泄物可以為稻田提供有機(jī)肥料[11]，水稻也能吸收利用魚飼料中的一部分氮元素，這些都有助于減少農(nóng)藥和化肥的使用，促進(jìn)稻田中物質(zhì)的相互利用和能量的流動[12]。如紫云英還田+稻魚共生模式不僅有效提高土壤有機(jī)質(zhì)、全氮、堿解氮、速效磷、速效鉀含量，還明顯提高土壤pH，加速土壤養(yǎng)分轉(zhuǎn)化和利用[13]。有研究發(fā)現(xiàn)養(yǎng)魚稻田土壤過氧化氫酶、磷酸酶、脲酶、蔗糖酶及脫氫酶活性分別高于對照樣地33.3%、19.3%、20.2%、39.8% 和72.7%[14]。稻田養(yǎng)魚提高土壤微生物的豐富度，增加微生物的功能種類，但對微生物數(shù)量、均勻度和代謝強(qiáng)度的影響不大[15]。上述研究[13-14] 多為大量飼料投入的稻魚綜合種養(yǎng)模式。大量投入飼料可能造成飼料殘留，水質(zhì)惡化，破壞生態(tài)平衡[16]，而傳統(tǒng)稻魚模式的外源投飼很少，魚類主要攝食雜草、底棲生物、水稻害蟲等[17]，對外源物質(zhì)的依賴減少，有助于實(shí)現(xiàn)農(nóng)業(yè)的可持續(xù)發(fā)展。如湖南辰溪、新化等地的稻魚模式均以系統(tǒng)自我運(yùn)行為主。

碳、氮、磷的變化特征及其化學(xué)計(jì)量比（C∶N、C∶P、N∶P） 可以有效衡量微生物營養(yǎng)需求與土壤有效養(yǎng)分的關(guān)系[18]。當(dāng)C∶N 較高時(shí)，植物所需的氮素將被微生物爭奪，滿足其生長；而C∶N 較低時(shí)，植物的氮素養(yǎng)分將不會被微生物消耗[19]。C∶P 較高時(shí)，微生物將吸收磷元素補(bǔ)充養(yǎng)分生長，植物對磷元素的有效性吸收將變差，這將在一定程度上約束其生長，而比值較低則反之[20]。主要依靠種群互利而維持的稻魚系統(tǒng)，其有機(jī)碳組分、有機(jī)碳庫穩(wěn)定性和土壤養(yǎng)分化學(xué)計(jì)量比如何變化，相關(guān)報(bào)道還比較少。

因此，本研究以湖南傳統(tǒng)稻魚養(yǎng)殖縣懷化市辰溪縣的長期低投入稻魚模式為研究對象，比較不同年限稻魚共作下稻田土壤養(yǎng)分儲量、有機(jī)碳組分、化學(xué)計(jì)量比的差異，探究長期稻魚模式下土壤特性變化的響應(yīng)因子，為稻田種養(yǎng)模式的系統(tǒng)運(yùn)行和自我維持提供理論依據(jù)，為優(yōu)化提升稻魚系統(tǒng)的生態(tài)功能提供參考。

1 材料與方法

1.1 研究地點(diǎn)

選取辰溪縣仙人灣清水塘村（27°44′58′′N，110°21′9′′E） 為研究對象，該地區(qū)屬于中亞熱帶季風(fēng)濕潤氣候，四季分明，雨量充足，熱量和水資源豐富，無霜期長。

1.2 樣品收集與處理

當(dāng)?shù)厮境Ｒ?guī)種植和稻魚共作種植在每年6 月進(jìn)行，均為一季稻種植。2 種模式的水稻秸稈均全量還田，且冬閑不泡水，水稻品種為‘玉針香’。

常規(guī)水稻種植施用有機(jī)肥和農(nóng)藥，田間水分進(jìn)行干濕交替管理；稻魚共作模式在水稻移栽10～20 d 后放魚苗，每公頃0.75～1.20 萬條，放魚后淹水20 cm 直至水稻成熟，且不施用化肥、農(nóng)藥，不投喂魚飼料。稻魚共作模式和水稻單作模式均施用2 250 kg·hm?2 有機(jī)肥，w（N+P2O5+K2O）≥8%，按基肥一次性全部施用。

根據(jù)對當(dāng)?shù)胤N養(yǎng)大戶的調(diào)研結(jié)果，選取種植密度為20 cm×30 cm、位置相鄰的稻魚共作15 和50 年（RF15、RF50） 的2 種傳統(tǒng)稻田養(yǎng)魚模式，以鄰近50 m 左右的水稻單作模式（RM） 為對照，于2023 年10 月在稻田排水、收獲后取樣。在每個試驗(yàn)點(diǎn)采集10 個約 3 m × 25 m 的試驗(yàn)塊作為重復(fù)，均按“S”形布設(shè)9 個采樣點(diǎn)，用土鉆自上而下分別采集土層0～20 cm 的土壤樣品。

1.3 指標(biāo)測定

土壤銨態(tài)氮、硝態(tài)氮、速效磷、全氮、全磷、有機(jī)碳含量，pH 按照《土壤農(nóng)業(yè)化學(xué)分析方法》[21] 測定；容重根據(jù)環(huán)刀法測定；土壤中礦質(zhì)結(jié)合態(tài)有機(jī)碳和顆粒有機(jī)碳含量參照章曉芳等[22] 的方法測定；易氧化有機(jī)碳含量參考顧曉娟等[23] 的方法測定。有機(jī)碳、全氮和全磷含量用于計(jì)算碳、氮、磷化學(xué)計(jì)量比（C∶N、C∶P、N∶P）。土壤養(yǎng)分儲量計(jì)算公式[24]如下：

土壤養(yǎng)分儲量=Ci×VW×H／10， （1）

式中，Ci 代表土壤營養(yǎng)成分含量，i 分別代表全氮、全磷或有機(jī)碳；VW 表示土壤容重，g/cm3；H 代表土層深度，cm，本研究為20 cm。

易氧化有機(jī)碳含量占比=易氧化有機(jī)碳含量=有機(jī)碳含量×100%， （2）

顆粒有機(jī)碳和礦質(zhì)結(jié)合態(tài)有機(jī)碳含量占比的計(jì)算方式與公式（2） 一致，將式中易氧化有機(jī)碳含量改為顆粒有機(jī)碳含量和礦質(zhì)結(jié)合態(tài)有機(jī)碳含量。

土壤碳庫管理指數(shù)（Carbon pool managementindex，CPMI） 的計(jì)算公式[25] 如下：

CPMI = CPI×CAI×100%， （3）

CPI = 樣品有機(jī)碳含量／對照有機(jī)碳含量， （4）

CAI = 樣品CPA／對照CPA， （5）

CPA = 活性有機(jī)碳含量／非活性有機(jī)碳含量， （6）

活性有機(jī)碳含量=易氧化有機(jī)碳含量+顆粒有機(jī)碳含量， （7）

非活性有機(jī)碳含量= 有機(jī)碳含量-活性有機(jī)碳含量，（8）

式中，CPI 代表碳庫指數(shù)（Carbon pool index），CAI代表碳庫活度指數(shù)（Carbon pool activity index），CPA 代表碳庫活度（Carbon pool activity）。

1.4 統(tǒng)計(jì)分析

使用Microsoft Excel 2019 進(jìn)行數(shù)據(jù)錄入和整理；使用IBM SPSS Statistics 24 進(jìn)行單因素方差分析，利用LSD（Least-significant difference） 進(jìn)行差異顯著性比較，Plt;0.05 時(shí)認(rèn)為存在顯著差異；所有圖用Origin 2024 繪制；使用Canoco 5 進(jìn)行冗余分析。

2 結(jié)果與分析

2.1 土壤碳、氮、磷養(yǎng)分含量以及化學(xué)計(jì)量比

如圖1 所示，土壤有機(jī)碳含量表現(xiàn)為RMgt;RF50gt;RF15。與RM 相比，RF15 的土壤有機(jī)碳含量降低19.14%，差異顯著（Plt;0.05）；RF50 略有降低，未達(dá)到顯著性差異（圖1A）。全氮含量為RF50gt;RMgt;RF15。較之RM，RF15 顯著降低47.84%，RF50顯著提高13.13%（Plt;0.05）（圖1B）。全磷含量表現(xiàn)為RMgt;RF15gt;RF50。較之RM，RF15、RF50 分別顯著降低11.35%、24.76%（Plt;0.05）（圖1C）。從土壤化學(xué)計(jì)量比來看，與RM 相比，RF15 的C∶N 顯著提高47.98%（圖1D），RF50 的C∶P 和N∶P 分別顯著提高33.08%、14.31%（Plt;0.05）（圖1E、1F）。

2.2 土壤速效養(yǎng)分含量和理化性質(zhì)

如表1 所示，不同稻魚共作年限和單作模式下，土壤速效養(yǎng)分含量和理化性質(zhì)差異顯著（Plt;0.05）。較之RM，RF15 的硝態(tài)氮、銨態(tài)氮和速效磷含量分別顯著降低52.56%、46.96% 和66.51%，pH 顯著增加13.15%（Plt;0.05），容重?zé)o顯著變化；RF50 的銨態(tài)氮、速效磷含量和pH 無明顯變化，硝態(tài)氮含量顯著提高51.52%，容重顯著降低11.11%（Plt;0.05）。稻魚共作15 年顯著降低土壤速效養(yǎng)分含量，但顯著提高pH；共作更久的RF50 的速效養(yǎng)分含量和pH 逐漸恢復(fù)到與RM 無顯著差異，但顯著降低容重。

2.3 土壤養(yǎng)分儲量

如圖2 所示，有機(jī)碳儲量表現(xiàn)為RMgt;RF50gt;RF15，與RM 相比，RF50 顯著降低19.37%（Plt;0.05），但比RF15 略有提升；氮儲量表現(xiàn)為RF50gt;RMgt;RF15，RF15 較之RM 顯著降低47.98%（Plt;0.05）， RF50 與RM 差異不顯著；磷儲量表現(xiàn)為RMgt;RF15gt;RF50，其中RF15 顯著降低11.69%，RF50 顯著降低33.08%（Plt;0.05）。結(jié)果表明，稻魚共作會顯著降低磷儲量，共作更久的RF能更好地維持碳儲量和氮儲量。

2.4 土壤有機(jī)碳組分含量及其占有機(jī)碳比例

如表2 所示，與RM 相比，稻魚共作（RF15、RF50） 顯著改變有機(jī)碳組分含量及其占土壤有機(jī)碳的比例（Plt;0.05），其中RF15 的易氧化有機(jī)碳含量及其占比分別顯著降低32.61%、16.34%，顆粒有機(jī)碳含量及其占比以及礦質(zhì)結(jié)合態(tài)有機(jī)碳占比分別顯著增加35.47%、72.24% 和31.68%；RF50 的礦質(zhì)結(jié)合態(tài)有機(jī)碳含量及其占比分別顯著降低36.00%、32.85%。結(jié)果表明RF15 能夠顯著提高穩(wěn)定性碳的占比。

2.5 土壤碳庫指數(shù)

如表3 所示，2 個年限的稻魚共作模式未顯著改變碳庫管理指數(shù)，但顯著改變碳庫活度、碳庫活度指數(shù)和碳庫指數(shù)（Plt;0.05）。與RM 相比，RF15 的碳庫活度和碳庫活度指數(shù)均顯著增加58.07%，碳庫指數(shù)顯著降低19.14%（Plt;0.05）；RF50 相關(guān)的碳庫指數(shù)均無顯著變化。稻魚共作15、50 年的碳庫管理指數(shù)整體保持穩(wěn)定，均大于100。

2.6 土壤有機(jī)碳庫穩(wěn)定性指標(biāo)與環(huán)境因子的關(guān)系50 較之RF15

對有機(jī)碳庫穩(wěn)定性相關(guān)指標(biāo)和環(huán)境因子指標(biāo)進(jìn)行冗余分析，結(jié)果如圖3 所示。其前兩軸的特征值分別為0.545 4 和0.234 1，累計(jì)解釋變異貢獻(xiàn)率為77.95%。其中易氧化有機(jī)碳含量占比與全氮、速效磷、硝態(tài)氮、銨態(tài)氮、易氧化有機(jī)碳、土壤有機(jī)碳含量，有機(jī)碳儲量，氮儲量，C∶P，N∶P 呈正相關(guān)關(guān)系；與顆粒有機(jī)碳、礦質(zhì)結(jié)合態(tài)有機(jī)碳、全磷含量，磷儲量，容重，pH，C∶N 呈負(fù)相關(guān)關(guān)系。顆粒有機(jī)碳含量占比、礦質(zhì)結(jié)合態(tài)有機(jī)碳含量占比和碳庫管理指數(shù)與全氮、速效磷、硝態(tài)氮、銨態(tài)氮、易氧化有機(jī)碳、土壤有機(jī)碳含量，有機(jī)碳儲量， C∶ P，N∶P 呈負(fù)相關(guān)，與全磷含量、磷儲量、容重、pH、C∶N 呈正相關(guān)。碳庫管理指數(shù)與顆粒有機(jī)碳含量占比、礦質(zhì)結(jié)合碳含量占比呈正相關(guān)，與易氧化有機(jī)碳含量占比呈負(fù)相關(guān)。其中，全氮含量、C∶N、有機(jī)碳儲量與土壤有機(jī)碳庫穩(wěn)定指標(biāo)顯著相關(guān)（Plt;0.05）。

3 討論與結(jié)論

3.1 長期低投入稻魚共作對土壤養(yǎng)分含量的影響

稻田土壤為水稻和魚生長提供養(yǎng)分和空間[26]。除了成土母質(zhì)、氣候、溫度等天然環(huán)境因素影響土壤理化性質(zhì)外，種植模式、肥料、農(nóng)藥、田間管理等人為因素也至關(guān)重要[27]。研究表明，較之常規(guī)純化肥投入，有機(jī)肥投入[28]、秸稈還田[29]、飼料投喂[30]也均可改變土壤養(yǎng)分含量。本研究為低飼料投入的稻魚共作模式，有機(jī)碳、硝態(tài)氮、銨態(tài)氮含量均在共作15 年時(shí)顯著降低，到共作50 年時(shí)又逐漸恢復(fù)到與對照差異不大的水平，全氮含量表現(xiàn)為先顯著降低又顯著增加。張琴等[31] 發(fā)現(xiàn)稻魚模式下有機(jī)碳含量在水稻進(jìn)入齊穗期后的第11 天顯著降低；李成芳[32] 同樣發(fā)現(xiàn)在稻鴨和稻魚種養(yǎng)模式下，土壤全氮含量在整個水稻生育期內(nèi)隨水稻生長而逐漸降低。而陳佳等[33] 發(fā)現(xiàn)有飼料投喂的稻魚共生模式使得土壤中全氮含量增加；周榆淇等[34] 經(jīng)連續(xù)兩年定點(diǎn)試驗(yàn)發(fā)現(xiàn)，稻魚模式下土壤速效養(yǎng)分含量變化無穩(wěn)定規(guī)律。本研究RF15 模式碳、氮含量降低可能是傳統(tǒng)農(nóng)戶稻田養(yǎng)魚模式投喂飼料極少，外源投入少而每年大量有機(jī)物（稻米和魚） 被帶離系統(tǒng)，使得土壤碳、氮含量逐漸降低。另外，稻魚共作長期淹水導(dǎo)致CO2 排放增加，也可能導(dǎo)致有機(jī)碳流失[35]，養(yǎng)分流失也更多[36]。較之RF15，RF50 的有機(jī)碳和全氮含量又逐漸提升，可能是更長期的稻魚模式將水中浮游動植物中的養(yǎng)分截留在系統(tǒng)中的量增大[37]，其次，共作年限增加逐步提高土壤中的酶活性[14， 34]、微生物多樣性和豐度，促進(jìn)土壤對養(yǎng)分的吸收能力[15， 38]。在本研究的稻魚共作模式（RF15、RF50） 中，全磷含量顯著降低，原因可能在于在南方酸性紅壤背景下，磷素的轉(zhuǎn)化效率低，更容易流失，與其他學(xué)者[39-41] 的研究結(jié)論相似。

3.2 長期低投入稻魚共作對土壤養(yǎng)分儲量和化學(xué)計(jì)量比的影響

土壤碳、氮、磷儲量是反映其生態(tài)系統(tǒng)服務(wù)功能的重要指標(biāo)[40-41]。養(yǎng)分儲量庫是所有養(yǎng)分的來源，其高低能夠反映出當(dāng)?shù)氐咎锿寥鲤B(yǎng)分流失、積累情況。本研究中RF15 的有機(jī)碳儲量、氮儲量和磷儲量均顯著降低；RF50 的有機(jī)碳和磷儲量顯著降低，氮儲量變化和全氮含量變化保持一致，其中RF50 有機(jī)碳儲量顯著降低的原因可能在于土壤容重顯著降低。本研究RF50 的氮儲量沒有顯著變化，原因可能是長期稻魚共作有機(jī)模式下，水體和土壤中浮游生物、微生物的多樣性增加，微生物殘?bào)w返回系統(tǒng)內(nèi)，逐步提升RF50 土壤氮儲量。

C ∶N 升高表明有機(jī)質(zhì)積累增加[ 4 2 ] ，本研究RF15 的C∶N 增加，但其有機(jī)碳和全氮含量均顯著降低，說明全氮含量的降低比例更大。而RF50 的C∶P 和N∶P 比均顯著增加，這會降低土壤中的磷素有效性[43]。

3.3 長期低投入稻魚共作對土壤有機(jī)碳庫組分及其穩(wěn)定性的影響

活性有機(jī)碳是指土壤中可在一定時(shí)間內(nèi)發(fā)生轉(zhuǎn)化或周轉(zhuǎn)、被微生物和植物直接利用、影響碳平衡的有機(jī)碳[44]。目前廣泛認(rèn)為有機(jī)碳是影響土壤肥力的關(guān)鍵因素，有機(jī)碳含量高的農(nóng)田結(jié)構(gòu)更好，被侵蝕風(fēng)險(xiǎn)更低[45]。土壤有機(jī)碳因其對全球可持續(xù)發(fā)展的多種服務(wù)和益處而日益受到關(guān)注[46]。碳庫管理指數(shù)被認(rèn)為是定量評價(jià)土壤質(zhì)量指數(shù)的最有效工具，可以作為衡量農(nóng)業(yè)管理措施影響土壤質(zhì)量的早期指標(biāo)，并有效評估農(nóng)業(yè)管理措施對可以在土壤中通過微生物作用相互轉(zhuǎn)化的惰性有機(jī)碳和礦質(zhì)結(jié)合態(tài)有機(jī)碳的影響[47-48]。礦質(zhì)結(jié)合態(tài)有機(jī)碳和活性有機(jī)碳的分配可以反映有機(jī)碳被微生物分解的難易程度[49]。本研究RF15 的易氧化有機(jī)碳含量及其占比顯著降低，顆粒有機(jī)碳含量及其占比以及礦質(zhì)結(jié)合態(tài)有機(jī)碳占比顯著增加；RF50 的礦質(zhì)結(jié)合態(tài)有機(jī)碳含量及其占比顯著降低?？傮w來看，稻魚共作不會造成土壤碳庫管理指數(shù)的顯著改變，其可能的原因是RF15 的礦質(zhì)結(jié)合態(tài)有機(jī)碳含量增加的同時(shí)，顆粒有機(jī)碳含量也在增加，碳庫管理指數(shù)也呈增加趨勢；RF50 的土壤碳庫管理指數(shù)沒有顯著改變，可能是還有其他活性有機(jī)碳含量也在同時(shí)降低。土壤中顆粒有機(jī)碳和礦質(zhì)結(jié)合態(tài)有機(jī)碳占比可在一定程度上反映土壤有機(jī)碳的穩(wěn)定程度，占比越大，土壤有機(jī)碳越易礦化、周轉(zhuǎn)期越短或活性越高，反之，土壤有機(jī)質(zhì)就越穩(wěn)定[ 5 0 ]。與RM 相比，RF15 和RF50 的有機(jī)碳含量降低，但是顆粒有機(jī)碳、礦質(zhì)結(jié)合態(tài)有機(jī)碳含量及其占比以及大部分土壤碳庫指數(shù)指標(biāo)升高，說明稻魚模式可以加速土壤有機(jī)質(zhì)周轉(zhuǎn)，并提高其活性。原因可能是稻魚模式引入了水生動物糞便作為碳源，以支持土壤微生物的生長和繁殖，從而提高土壤碳周轉(zhuǎn)率[51]。長期低投入外源性營養(yǎng)物質(zhì)的稻魚模式能夠保持良好的碳庫穩(wěn)定性，在保證土壤質(zhì)量的前提下，提高稻田系統(tǒng)的物質(zhì)產(chǎn)出。

3.4 結(jié)論

1） 本研究通過傳統(tǒng)化學(xué)方法比較水稻單作、稻魚共作15 年、稻魚共作50 年硝態(tài)氮、銨態(tài)氮和速效磷含量和全量養(yǎng)分含量的變化，證明了長期低投入循環(huán)農(nóng)業(yè)模式穩(wěn)定全量養(yǎng)分含量的能力。

2） 稻魚共作15、50 年的碳庫管理指數(shù)整體保持穩(wěn)定，且均大于100，表明長期稻魚種養(yǎng)提升了土壤質(zhì)量，進(jìn)而改善了稻田生態(tài)環(huán)境。

3） 對有機(jī)碳穩(wěn)定性相關(guān)指標(biāo)和環(huán)進(jìn)行冗余分析，揭示在長期稻魚模式下，與土壤有機(jī)碳庫穩(wěn)定指標(biāo)顯著相關(guān)的環(huán)境因子為全氮含量、C∶N、有機(jī)碳儲量。

參考文獻(xiàn)：境因子指標(biāo)

[1]鄭華斌， 賀慧， 姚林， 等. 稻田飼養(yǎng)動物的生態(tài)經(jīng)濟(jì)效應(yīng)及其應(yīng)用前景[J]. 濕地科學(xué)， 2015， 13（4）： 510-517.

[2]ISBELL F， CRAVEN D， CONNOLLY J， et al. Biodiversityincreases the resistance of ecosystem productivityto climate extremes[J]. Nature， 2015， 526（7574）： 574-577.

[3]RENARD D， TILMAN D. National food production stabilizedby crop diversity[J]. Nature， 2019， 571（7764）：257-260.

[4]于秀娟， 郝向舉， 黨子喬， 等. 中國稻漁綜合種養(yǎng)產(chǎn)業(yè)發(fā)展報(bào)告（2023）[J]. 中國水產(chǎn)， 2023（8）： 19-26.

[5]WANG B， SUN Y， JIAO W. Ecological benefit evaluationof agricultural heritage system conservation： A casestudy of the Qingtian Rice-Fish Culture System[J]. Journalof Resources and Ecology， 2021， 12（4）： 489-497.

[6]丁姣龍. 稻魚共生年限對耕層土壤養(yǎng)分與微生物影響研究[D]. 長沙： 湖南農(nóng)業(yè)大學(xué)， 2021.

[7]佀國涵， 袁家富， 彭成林， 等. 長期稻蝦共作模式提高稻田土壤生物肥力的機(jī)理[J]. 植物營養(yǎng)與肥料學(xué)報(bào)，2020， 26（12）： 2168-2176.

[8]張苗苗， 宗良綱， 謝桐洲. 有機(jī)稻鴨共作對土壤養(yǎng)分動態(tài)變化和經(jīng)濟(jì)效益的影響[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào)，2010， 18（2）： 256-260.

[9]朱練峰， 房偉平， 莊雪浩， 等. 稻蛙共作對土壤理化特性和水稻產(chǎn)量的影響[J]. 中國稻米， 2023， 29（5）： 23-27.

[10]陳曉云， 孫文濤， 于鳳泉， 等. 稻蟹生態(tài)種養(yǎng)模式對稻田土壤肥力及生產(chǎn)效益的影響[J]. 土壤通報(bào)， 2021， 52（5）：1165-1172.

[11]劉元生， 孟慶紅， 何騰兵， 等. 稻田生態(tài)養(yǎng)魚水質(zhì)動態(tài)與水稻生長及經(jīng)濟(jì)效益研究[J]. 耕作與栽培， 2003（5）： 5-6.

[12]XIE J， HU L， TANG J， et al. Ecological mechanisms underlyingthe sustainability of the agricultural heritagerice-fish coculture system[J]. Proceedings of the NationalAcademy of Sciences of the United States of America，2011， 108（50）： E1381-E1387.

[13]呂廣動， 黃璜， 梁玉剛， 等. 紫云英還田+稻魚共生對水稻土壤養(yǎng)分及產(chǎn)量的影響[J]. 西南農(nóng)業(yè)學(xué)報(bào)， 2020，33（8）： 1729-1735.

[14]孫剛， 房巖， 王欣， 等. 稻魚復(fù)合種養(yǎng)對水田土壤酶活性的影響[J]. 農(nóng)業(yè)與技術(shù)， 2009， 29（1）： 23-26.

[15]展茗， 曹湊貴， 汪金平， 等. 稻鴨、稻魚復(fù)合生態(tài)系統(tǒng)土壤微生物多樣性特征分析[J]. 土壤學(xué)報(bào)， 2008， 45（6）：1179-1183.

[16]JIA B， TANG Y， TIAN L， et al. Impact of fish farmingon phosphorus in reservoir sediments[J]. Scientific Reports，2015， 5： 16617.

[17]閔寬洪， 郁桐炳. 浙江青田“稻魚共生”系統(tǒng)發(fā)展的新模式： 從傳統(tǒng)田魚生產(chǎn)到現(xiàn)代漁業(yè)文化產(chǎn)業(yè)[J]. 中國漁業(yè)經(jīng)濟(jì)， 2009， 27（1）： 25-28.

[18]ROSINGER C， ROUSK J， SANDéN H. Can enzymaticstoichiometry be used to determine growth-limiting nutrientsfor microorganisms？： A critical assessment in twosubtropical soils[J]. Soil Biology and Biochemistry，2019， 128： 115-126.

[19]CHAPIN F S， MATSON P A， MOONEY H A， et al.Principles of terrestrial ecosystem ecology[M]. New York：Springer， 2002.

[20]DENG Q， CHENG X， HUI D， et al. Soil microbial communityand its interaction with soil carbon and nitrogendynamics following afforestation in central China[J]. Scienceof the Total Environment， 2016， 541： 230-237.

[21]魯如坤. 土壤農(nóng)業(yè)化學(xué)分析方法[M]. 北京： 中國農(nóng)業(yè)科技出版社， 2000.

[22]章曉芳， 鄭生猛， 夏銀行， 等. 紅壤丘陵區(qū)土壤有機(jī)碳組分對土地利用方式的響應(yīng)特征[J]. 環(huán)境科學(xué)， 2020，41（3）： 1466-1473.

[23]顧曉娟， 于耀泓， 劉悅， 等. 華南典型人工林對土壤易氧化有機(jī)碳的影響[J]. 森林與環(huán)境學(xué)報(bào)， 2023， 43（2）： 145-151.

[24]LUO X， HOU E， CHEN J， et al. Dynamics of carbon， nitrogen，and phosphorus stocks and stoichiometry resultingfrom conversion of primary broadleaf forest to plantationand secondary forest in subtropical China[J].CATENA， 2020， 193： 104606.

[25]BLAIR G J， LEFROY R D， LISLE L. Soil carbon fractionsbased on their degree of oxidation， and the developmentof a carbon management index for agricultural systems[J]. Australian Journal of Agricultural Research，1995， 46（7）： 1459-1466.

[26]林秋燕， 鐘全林， 李寶銀， 等. 氮磷添加對刨花楠幼林生長與葉性狀及土壤養(yǎng)分關(guān)系的影響[J]. 應(yīng)用與環(huán)境生物學(xué)報(bào)， 2023， 29（5）： 1068-1076.

[27]李科江， 馬俊永， 曹彩云， 等. 長期定位施用不同種類有機(jī)肥對作物產(chǎn)量及潮土理化性狀的影響[J]. 河北農(nóng)業(yè)科學(xué)， 2007（1）： 60-63.

[28]何偉， 王會， 韓飛， 等. 長期施用有機(jī)肥顯著提升潮土有機(jī)碳組分[J]. 土壤學(xué)報(bào)， 2020， 57（2）： 425-434.

[29]齊智娟， 徐敬文， 張忠學(xué)， 等. 秸稈還田配施氮肥對黑土玉米田土壤CO2 排放與碳平衡的影響[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2024， 55（6）： 284-293.

[30]沈建凱， 黃璜， 傅志強(qiáng)， 等. 稻鴨生態(tài)種養(yǎng)系統(tǒng)直播水稻根表和根際土壤營養(yǎng)特性研究[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào)，2010， 18（6）： 1151-1156.

[31]張琴， 孟祥杰， 陳燦， 等. 不同稻、鴨、魚共生模式對土壤質(zhì)量及水稻生長發(fā)育的影響[J]. 江蘇農(nóng)業(yè)科學(xué)，2023， 51（20）： 58-67.

[32]李成芳. 稻田生態(tài)種養(yǎng)模式氮素轉(zhuǎn)化規(guī)律的研究[D].武漢： 華中農(nóng)業(yè)大學(xué)， 2008.

[33]陳佳， 趙璐峰， 戴然欣， 等. 稻魚共生系統(tǒng)的土壤產(chǎn)甲烷和甲烷氧化微生物群落[J]. 生態(tài)學(xué)雜志， 2023， 42（12）：2961-2971.

[34]周榆淇， 鄒冬生， 王安崠， 等. 長期淹水條件下稻魚復(fù)合種養(yǎng)對土壤養(yǎng)分和酶活性的影響[J]. 農(nóng)業(yè)現(xiàn)代化研究，2022， 43（5）： 911-920.

[35]GREENWAY H， ARMSTRONG W， COLMER T D.Conditions leading to high CO2 （gt; 5 kPa） in waterloggedfloodedsoils and possible effects on root growth andmetabolism[J]. Annals of Botany， 2006， 98（1）： 9-32.

[36]馮國祿， 李書迪， 許尤厚， 等. 撒施液體復(fù)合肥后不同蓄水深度的水分管理對稻田養(yǎng)分流失潛力的影響[J]. 中國土壤與肥料， 2018（1）： 83-86.

[37]張丹， 閔慶文， 成升魁， 等. 應(yīng)用碳、氮穩(wěn)定同位素研究稻田多個物種共存的食物網(wǎng)結(jié)構(gòu)和營養(yǎng)級關(guān)系[J]. 生態(tài)學(xué)報(bào)， 2010， 30（24）： 6734-6740.

[38]田興， 鄒利， 鄧時(shí)銘， 等. 養(yǎng)殖禾花鯉對稻田土壤養(yǎng)分、微生物多樣性的影響[J]. 江蘇農(nóng)業(yè)科學(xué)， 2023， 51（17）：191-197.

[39]馬微微， 陳燦， 黃璜， 等. 壟作稻魚雞共生對稻田土壤養(yǎng)分含量及水稻產(chǎn)量的影響[J]. 河南農(nóng)業(yè)科學(xué)， 2021，50（8）： 9-17.

[40]宋希娟， 楊成德， 陳秀蓉， 等. 東祁連山高寒草地生態(tài)系統(tǒng)N、P 養(yǎng)分含量研究[J]. 草原與草坪， 2008（6）： 46-49.

[41]鐘建軍， 李林， 魏識廣， 等. 漓江流域喀斯特森林土壤碳氮磷儲量分布特征及其影響因素[J]. 水土保持學(xué)報(bào)，2023， 37（6）： 180-186.

[42]WANG M， CHEN H， ZHANG W， et al. Soil nutrientsand stoichiometric ratios as affected by land use and lithologyat county scale in a karst area， southwestChina[J]. Science of the Total Environment， 2018， 619：1299-1307.

[43]HUANG J， LIU L， QI K， et al. Differential mechanismsdrive changes in soil C pools under N and P enrichment in a subalpine spruce plantation[J]. Geoderma， 2019， 340：213-223.

[44]沈宏， 曹志洪， 胡正義. 土壤活性有機(jī)碳的表征及其生態(tài)效應(yīng)[J]. 生態(tài)學(xué)雜志， 1999， 18（3）： 32-38.

[45]LIU Y， GE T， VAN GROENIGEN K J， et al. Rice paddysoils are a quantitatively important carbon store accordingto a global synthesis[J]. Communications Earth amp; Environment，2021， 2（1）： 154.

[46]MANLAY R J， FELLER C， SWIFT M. Historical evolutionof soil organic matter concepts and their relationshipswith the fertility and sustainability of cropping systems[J]. Agriculture Ecosystems amp; Environment， 2007，119（3/4）： 217-233.

[47]舒業(yè)勤， 彭復(fù)細(xì)， 雷文碩， 等. 稻油復(fù)種不同措施下土壤有機(jī)碳組分積累及其穩(wěn)定性特征[J/OL]. 土壤學(xué)報(bào)（2024-04-08） [2024-05-10]. https：//link.cnki.net/urlid/32.1119.p.20240402.1444.002.

[48]薛志婧， 李霄云， 焦磊， 等. 土壤礦質(zhì)結(jié)合態(tài)有機(jī)碳形成及穩(wěn)定機(jī)制的研究進(jìn)展[J]. 水土保持學(xué)報(bào)， 2023， 37（5）：12-23.

[49]COTRUFO M F， RANALLI M G， HADDIX M L， et al.Soil carbon storage informed by particulate and mineralassociatedorganic matter[J]. Nature Geoscience， 2019，12（12）： 989-994.

[50]藍(lán)賢瑾， 呂真真， 劉秀梅， 等. 長期施肥對紅壤性水稻土顆粒有機(jī)質(zhì)和礦物結(jié)合態(tài)有機(jī)質(zhì)含量與化學(xué)組成的影響[J]. 土壤， 2021， 53（1）： 140-147.

[51]WANG C， YANG Q， CHEN J， et al. Variations in soilorganic carbon fractions and microbial community in ricefields under an integrated cropping system[J]. Agronomy，2024， 14（1）： 81.

【責(zé)任編輯 李慶玲】

基金項(xiàng)目：國家自然科學(xué)基金區(qū)域創(chuàng)新發(fā)展聯(lián)合基金重點(diǎn)項(xiàng)目（U21A20184）