杜兆國　肖洋　張瑞豪　徐金忠　邵社剛　倪棟　張振燁

［關(guān)鍵詞］土壤酶；侵蝕溝；坡耕地；土壤理化性質(zhì)；東北黑土區(qū)

［摘要］東北黑土區(qū)作為我國重要的糧食產(chǎn)區(qū)，侵蝕溝發(fā)育導(dǎo)致土壤地力下降，嚴(yán)重制約了黑土區(qū)農(nóng)業(yè)發(fā)展。為了分析黑土區(qū)侵蝕溝發(fā)育對(duì)表層土壤酶活性的影響，探究酶活性對(duì)侵蝕-沉積作用的響應(yīng)機(jī)制，選取黑龍江省哈爾濱市延壽縣雙奎河雙安村河段不同發(fā)育程度的侵蝕溝作為研究對(duì)象，測定不同發(fā)育程度侵蝕溝溝頭、溝中、溝尾段0～40 cm土層的堿性磷酸酶、脲酶、蔗糖酶、過氧化氫酶活性和土壤理化性質(zhì)，并進(jìn)行相關(guān)性分析和冗余分析，結(jié)果表明：①隨著侵蝕發(fā)育，土壤酶活性總體呈下降趨勢(shì)，隨著土層加深，土壤酶活性降低；②4類酶活性和土壤養(yǎng)分指標(biāo)呈顯著正相關(guān)、和土壤密度呈顯著負(fù)相關(guān)，有機(jī)質(zhì)含量是影響土壤酶活性最主要的因子；③在侵蝕-沉積作用下，土壤酶活性基本呈現(xiàn)溝頭、溝中段隨侵蝕發(fā)育而下降明顯，溝尾段因沉積作用而下降減緩的規(guī)律。

［中圖分類號(hào)］ S154；S157.1［文獻(xiàn)標(biāo)識(shí)碼］ A［文章編號(hào)］ 1000－0941（2023）08－0038－06

微生物活動(dòng)釋放的酶參與了土壤有機(jī)質(zhì)和養(yǎng)分的形成[1]、循環(huán)[2]、分解[3]、礦化[4]等諸多環(huán)節(jié)，其催化活動(dòng)是土壤生物化學(xué)反應(yīng)的重要一環(huán)。酶的活性與土壤環(huán)境息息相關(guān)，因土壤酶對(duì)外界氣候、地形、植被等環(huán)境因素具備較強(qiáng)敏感性，故成為表征土壤理化性質(zhì)狀況的重要指標(biāo)[5]。測定酶的活性，可以了解特定養(yǎng)分在土壤中的轉(zhuǎn)化速率，掌握土壤生物化學(xué)反應(yīng)的相對(duì)強(qiáng)度[6]。侵蝕作為土壤退化的重要表現(xiàn)形式之一，不僅在宏觀上生成溝道改變地形，而且直接影響了土壤生態(tài)的微觀環(huán)境[7-8]。大量研究表明侵蝕溝的演變對(duì)土壤理化性質(zhì)有顯著影響，造成土壤粒徑分布和土壤內(nèi)部溫度、含水量等水熱條件的變化[9]，進(jìn)而改變土壤酶活性[10-12]。

侵蝕作用會(huì)打破土壤養(yǎng)分和微生物之間的良性循環(huán)，導(dǎo)致土壤地力下降，影響糧食產(chǎn)量[13-14]。東北黑土區(qū)作為我國重要的糧食產(chǎn)區(qū)，每年因侵蝕造成的糧食損失多達(dá)280萬t，侵蝕溝的發(fā)育導(dǎo)致土壤地力下降，嚴(yán)重制約了黑土區(qū)農(nóng)業(yè)發(fā)展。本研究以東北黑土區(qū)不同發(fā)育程度的侵蝕溝為研究對(duì)象，通過量化分析不同侵蝕發(fā)育程度侵蝕溝溝頭、溝中、溝尾段0～40 cm表層土壤酶活性分布狀況，探討侵蝕作用對(duì)表層黑土土壤酶活性的影響，以期為黑土區(qū)侵蝕溝治理工作提供支撐。

1研究區(qū)概況

研究區(qū)雙奎河雙安村河段位于黑龍江省哈爾濱市延壽縣，屬松嫩平原和小興安嶺余脈的過渡區(qū)，以低丘地形為主，海拔145～300 m。屬寒溫帶大陸性氣候區(qū)，年均氣溫2.8 ℃，年均降水量571.7 mm。土壤類型為黑土，有機(jī)質(zhì)含量1%～8%，土壤質(zhì)地為壤土和黏壤土，底層成土母質(zhì)為第四紀(jì)黃土狀亞土。經(jīng)現(xiàn)場實(shí)測，地表下垂直高度20～25 cm處存在明顯的犁底層，平均厚度約4 cm。

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

2.1土樣采集

為保證室外環(huán)境下試驗(yàn)變量的一致性，在雙奎河雙安村河段選取一條支流（流域面積0.13～1.3 km2），采集時(shí)間為2021年5月春耕前，此時(shí)無積雪、降水少，氣候、人為因素對(duì)土壤的擾動(dòng)較小。此外要確保采集土樣的海拔、坡向、坡度等地形因素相近。采樣時(shí)侵蝕溝地貌中主要分布的植物種類包括蔞蒿、薹草、杞柳、紅足蒿、青蒿等。

依據(jù)黑龍江省地方標(biāo)準(zhǔn)《侵蝕溝分級(jí)與分類》（DB23/T 2412—2019），選取流域內(nèi)3條不同發(fā)育程度的侵蝕溝，分別記為淺溝A1、小型半穩(wěn)定溝A2、中型穩(wěn)定切溝A3（見圖1），侵蝕發(fā)育程度為A1＜A2＜A3。各溝道的坡面沿縱斷面線分為溝頭、溝中、溝尾3個(gè)部位，分別在每個(gè)部位的坡面位置設(shè)置2個(gè)溝坡采樣區(qū)，在溝底集水區(qū)設(shè)置1個(gè)溝底采樣區(qū)（見圖2），采樣區(qū)大小為0.5 m×0.5 m。溝坡采樣時(shí)需垂直于坡面，溝底采樣前需去除厚3 cm的上游沖刷浮土。每個(gè)采樣區(qū)在采集時(shí)以地表為基準(zhǔn)，按照機(jī)械分層法垂直向下分別采集0～20 cm、＞20～40 cm兩個(gè)土層，采集樣方尺寸為長15 cm×寬15 cm×深20 cm，重復(fù)采樣3次。在溝道溝頭間設(shè)置未發(fā)生侵蝕的空白對(duì)照組CK，同樣在0～20 cm、＞20～40 cm兩個(gè)土層采樣3次后進(jìn)行混樣處理。采樣點(diǎn)基本情況見表1。

2.2土樣測定

采集的土樣使用塑料盒妥善保存，3份重復(fù)樣品中：1份存儲(chǔ)于-20 ℃的冰柜中用于測定土壤微生物量，參與酶活性計(jì)算；另外2份在避光處攤開陰干、過篩，分別用于測定土壤理化性質(zhì)和酶活性。土壤理化性質(zhì)和酶活性的測定方法見表2，測定得到的土壤理化性質(zhì)結(jié)果見表3。

2.3數(shù)據(jù)分析

完成土樣采集和測定后，采用SPSS 25、Excel 2020軟件對(duì)測定數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析和相關(guān)性分析，采用Canoco 5.0軟件進(jìn)行冗余分析。

3結(jié)果

3.1不同發(fā)育程度侵蝕溝的土壤酶活性特征

3.1.1堿性磷酸酶

不同發(fā)育程度溝道間堿性磷酸酶活性變化顯著（p＜0.05）（見圖3）。表4是不同發(fā)育程度侵蝕溝溝頭、溝中、溝尾段0～40 cm土層的堿性磷酸酶活性，在0～20 cm土層，A1、A2、A3的堿性磷酸酶活性分別為0.51、0.54、0.20 μg/（g·h），各溝道堿性磷酸酶活性相較于CK分別下降17.74%、12.90%和67.74%，總體呈隨侵蝕發(fā)育波動(dòng)下降趨勢(shì)；在＞20～40 cm土層，A1、A2、A3的堿性磷酸酶活性分別為0.30、0.24、0.12 μg/（g·h），各溝道堿性磷酸酶活性呈CK＞A1＞A2＞A3，隨著侵蝕發(fā)育，堿性磷酸酶活性不斷下降。各溝道0～20 cm土層堿性磷酸酶活性總體高于＞20～40 cm土層，且二者之間差異極顯著（p＜0.01），表明隨土層加深，堿性磷酸酶活性下降。從空間分布上來看，A1、A2、A3中0～20 cm土層和＞20～40 cm土層的堿性磷酸酶活性最大值均出現(xiàn)在溝尾段，且在沉積作用下，堿性磷酸酶活性在溝尾段隨侵蝕發(fā)育而降幅減緩，而在溝頭、溝中段則隨著侵蝕溝發(fā)育降幅明顯。

3.1.2脲酶

不同發(fā)育程度溝道間土壤脲酶活性變化顯著（p＜0.05）（見圖4）。表5是不同發(fā)育程度侵蝕溝溝頭、溝中、溝尾段0～40 cm土層的脲酶活性，在0～20 cm土層，A1、A2、A3的脲酶活性分別為1.59、1.08、1.02 μg/（g·h），分別較CK下降了4.22%、34.94%、38.55%，在＞20～40 cm土層，A1、A2、A3的脲酶活性分別為1.07、0.65、0.61 μg/（g·h），分別較CK下降了15.08%、48.41%、51.59%。0～20 cm和＞20～40 cm土層脲酶活性均呈現(xiàn)CK＞A1＞A2＞A3，脲酶活性隨著侵蝕發(fā)育呈遞減趨勢(shì)。各溝道0～20 cm土層脲酶活性總體高于＞20～40 cm土層，且二者之間差異極顯著（p＜0.01），表明隨土層加深，脲酶活性下降；且兩個(gè)土層間的脲酶活性變化率A3＞A2＞A1，表明侵蝕溝的發(fā)育加速了脲酶活性的下降。各溝道＞20～40 cm土層脲酶活性的變異系數(shù)高于0～20 cm土層，表明在＞20～40 cm土層脲酶活性受侵蝕發(fā)育的影響更大。從空間分布上看，A1、A2、A3溝尾段脲酶活性較CK的降幅最小，表明隨著侵蝕發(fā)育，溝尾段脲酶活性下降最少，溝頭段下降最多。

3.1.3蔗糖酶

不同發(fā)育程度溝道間土壤蔗糖酶活性變化不顯著（見圖5）。表6是不同發(fā)育程度侵蝕溝溝頭、溝中、溝尾段0～40 cm土層的蔗糖酶活性，在0～20 cm土層，A1、A2、A3的蔗糖酶活性分別為47.37、48.55、44.28 μg/（g·h），呈CK＞A2＞A1＞A3的趨勢(shì)；在＞20～40 cm土層，A1、A2、A3的蔗糖酶活性分別為26.68、24.31、31.71 μg/（g·h），分別較CK下降了18.78%、26.00%、3.44%。各溝道0～20 cm土層蔗糖酶活性總體高于＞20～40 cm土層，且二者之間差異極顯著（p＜0.01），表明隨土層加深，蔗糖酶活性下降。從空間分布上看，在0～20 cm土層，溝頭、溝中段蔗糖酶活性隨侵蝕發(fā)育呈下降趨勢(shì)，而溝尾段呈上升趨勢(shì)，且A2、A3溝尾段的蔗糖酶活性均高于CK，表明蔗糖酶活性在溝頭、溝中段受到養(yǎng)分遷移的影響，在侵蝕區(qū)活性下降，而在溝尾段養(yǎng)分沉積，蔗糖酶活性隨侵蝕發(fā)育而增加。

3.1.4過氧化氫酶

不同發(fā)育程度溝道間土壤過氧化氫酶活性變化顯著（p＜0.05）（見圖6）。表7是不同發(fā)育程度侵蝕溝溝頭、溝中、溝尾段0～40 cm土層的過氧化氫酶活性，在0～20 cm土層，A1、A2、A3的過氧化氫酶活性分別為2.85、2.28、2.32 μg/（g·h），分別較CK下降了8.65%、26.92%、25.64%，表明過氧化氫酶活性隨侵蝕發(fā)育呈波動(dòng)下降趨勢(shì)；在＞20～40 cm土層，A1、A2、A3的過氧化氫酶活性分別為2.05、1.28、1.68 μg/（g·h），過氧化氫酶活性整體呈現(xiàn)CK＞A1＞A3＞A2，過氧化氫酶活性最小值出現(xiàn)在A2。各溝道0～20 cm土層過氧化氫酶活性總體高于＞20～40 cm土層，且二者之間差異顯著（p＜0.05），隨土層加深，過氧化氫酶活性下降。

3.2酶活性的影響因素分析

3.2.1酶活性與土壤理化性質(zhì)的關(guān)系

將各項(xiàng)酶活性和土壤理化性質(zhì)指標(biāo)采用SPSS軟件進(jìn)行探索性因子分析，KMO（Kaiser-Meyer-Olkin）檢驗(yàn)統(tǒng)計(jì)量為0.759（p＜0.001），表明變量之間相關(guān)性較強(qiáng)。表8是酶活性和土壤理化性質(zhì)指標(biāo)的相關(guān)性分析結(jié)果，堿性磷酸酶、脲酶、蔗糖酶、過氧化氫酶4類酶活性均與有機(jī)質(zhì)、全氮、全磷、全鉀、速效鉀、有效磷、堿解氮7項(xiàng)土壤養(yǎng)分指標(biāo)呈顯著正相關(guān)（p＜0.05），表明隨著侵蝕發(fā)育，土壤養(yǎng)分逐漸減少，酶活性隨之下降。4類酶活性與土壤密度呈顯著負(fù)相關(guān)（p＜0.05），隨著土壤密度的增加，酶活性下降。堿性磷酸酶活性與土壤含水量呈顯著正相關(guān)（p＜0.01），其余3類酶活性與土壤含水量之間的關(guān)系不顯著。

酶活性和土壤理化性質(zhì)的冗余分析結(jié)果見圖7。在0～40 cm土層中，第一軸與第二軸分別解釋了土壤酶活性和土壤理化性質(zhì)變化的51.8%和4.34%，總解釋度為56.14%。在各項(xiàng)理化性質(zhì)指標(biāo)中，土壤有機(jī)質(zhì)含量是影響土壤酶活性的最主要因子，其解釋度為46.0%（p=0.002），表明有機(jī)質(zhì)的分布極大影響著侵蝕溝道內(nèi)生物酶的活性；其次是全氮含量與速效鉀含量，解釋度分別為3.2%（p=0.044）和2.8%（p=0.048）。同時(shí)堿性磷酸酶、脲酶、蔗糖酶、過氧化氫酶4類酶活性指標(biāo)均與有機(jī)質(zhì)含量、全氮含量與速效鉀含量呈正相關(guān)關(guān)系，表明隨著侵蝕溝發(fā)育，有機(jī)質(zhì)、全氮、速效鉀等土壤肥力的遷移進(jìn)一步影響了土壤中生物酶的活性。

3.2.2酶活性與土壤侵蝕的關(guān)系

堿性磷酸酶、脲酶、蔗糖酶、過氧化氫酶與土壤養(yǎng)分的變化均具有相對(duì)專一性，可以比較準(zhǔn)確地反映土壤中有機(jī)質(zhì)和氮、磷、鉀的含量。酶活性主要受土壤養(yǎng)分中的有機(jī)質(zhì)、全氮和速效鉀含量的影響，而侵蝕溝內(nèi)各部位的酶活性也受土壤養(yǎng)分遷移的影響。從空間分布上來看，在侵蝕-沉積作用下，溝頭、溝中段被侵蝕的土壤裹挾土壤養(yǎng)分遷移，在溝尾段沉積；受土壤養(yǎng)分遷移的影響，土壤酶活性基本呈現(xiàn)出溝頭、溝中段隨侵蝕發(fā)育而下降明顯，溝尾段因沉積作用而下降減緩的規(guī)律。

黑土區(qū)侵蝕溝的發(fā)育主要分為3個(gè)方面：一是溯源侵蝕。主要是由地表徑流驅(qū)動(dòng)，當(dāng)降水強(qiáng)度超過土壤入滲速度時(shí)形成徑流，徑流在集水區(qū)匯集，沖蝕沿途的土壤，造成水土流失，同時(shí)土壤養(yǎng)分也隨徑流遷移，溝頭段土壤隨侵蝕發(fā)育逐漸“粗骨化”“貧瘠化”，使得微生物失去營養(yǎng)來源，造成溝頭段酶活性的下降，堿性磷酸酶、脲酶、蔗糖酶活性在溝頭段隨侵蝕發(fā)育呈下降趨勢(shì)均能證明這個(gè)規(guī)律。二是溝壁崩塌。溝壁崩塌在淺溝尚不明顯，當(dāng)侵蝕溝發(fā)育為切溝時(shí)，坡面不同部位的土壤理化性質(zhì)差異逐漸明顯，坡面頂部受雨滴濺蝕，土壤機(jī)械位移較少，而坡面上、中部同時(shí)受到徑流沖刷，土壤養(yǎng)分被徑流剝離并向下沉積，造成坡面土壤養(yǎng)分的流失和坡腳處土壤養(yǎng)分的累積，使得土壤養(yǎng)分在坡面呈現(xiàn)多—少—多的分布規(guī)律，酶活性則呈現(xiàn)相同分布規(guī)律。三是溝底下切。侵蝕溝發(fā)育時(shí)，表層土壤質(zhì)地疏松，且微生物活動(dòng)劇烈的肥沃土壤被沖蝕，裸露出相對(duì)貧瘠緊實(shí)的下層土層，加上黑土分層明顯，使得酶活性在不同土層間差異顯著，本研究中堿性磷酸酶、脲酶、蔗糖酶、過氧化氫酶4類酶活性在0～20 cm和＞20～40 cm土層間均差異顯著。

4結(jié)論

1）除蔗糖酶外，不同發(fā)育程度侵蝕溝間土壤堿性磷酸酶活性、脲酶活性、過氧化氫酶活性差異顯著（p<0.05），4類酶活性在0～20 cm和＞20～40 cm土層間均差異顯著（p<0.01）。隨著侵蝕發(fā)育，土壤酶活性總體呈下降趨勢(shì)；隨著土層加深，土壤酶活性降低。

2）4類酶活性均與有機(jī)質(zhì)、全氮、全磷、全鉀、速效鉀、有效磷、堿解氮7項(xiàng)土壤養(yǎng)分指標(biāo)呈顯著正相關(guān)（p＜0.05），與土壤密度呈顯著負(fù)相關(guān)（p＜0.05）。冗余分析結(jié)果表明，有機(jī)質(zhì)含量是影響土壤酶活性最主要的因子，其解釋度為46.0%（p=0.002），其次是全氮含量和速效鉀含量，解釋度分別為3.2%（p=0.044）和2.8%（p=0.048）。

3）在侵蝕-沉積作用下，溝頭、溝中段被侵蝕的土壤裹挾土壤養(yǎng)分遷移，在溝尾段沉積；受到土壤養(yǎng)分遷移的影響，土壤酶活性基本呈現(xiàn)出溝頭、溝中段隨侵蝕發(fā)育而下降明顯，溝尾段因沉積作用而下降減緩的規(guī)律。

[參考文獻(xiàn)]

[1] 史吉平，張夫道，林葆.長期施肥對(duì)土壤有機(jī)質(zhì)及生物學(xué)特性的影響[J].土壤肥料，1998（3）：7-11.

[2] 沈仁芳，趙學(xué)強(qiáng).土壤微生物在植物獲得養(yǎng)分中的作用[J].生態(tài)學(xué)報(bào)，2015，35（20）：6584-6591.

[3] 王戰(zhàn)臺(tái)，曹旻霞，蔡文昌，等.土地利用類型對(duì)侵蝕黃壤腐殖質(zhì)及土壤酶的影響[J].水土保持通報(bào)，2017，37（5）：27-33.

[4] 辛芝紅，李君劍，趙小娜，等.煤礦區(qū)不同復(fù)墾年限的土壤有機(jī)碳礦化和酶活性特征[J].環(huán)境科學(xué)研究，2017，30（10）：1580-1586.

[5] 孟立君，吳鳳芝.土壤酶研究進(jìn)展[J].東北農(nóng)業(yè)大學(xué)學(xué)報(bào)，2004（5）：622-626.

[6] 郭曼，鄭粉莉，和文祥，等.黃土丘陵區(qū)不同退耕年限植被多樣性變化及其與土壤養(yǎng)分和酶活性的關(guān)系[J].土壤學(xué)報(bào)，2010，47（5）：979-986.

[7] 左自泉，張斌，覃發(fā)超，等.東北黑土區(qū)侵蝕溝長度規(guī)模-頻率關(guān)系與空間分布規(guī)律[J].西華師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2022，43（3）：322-327.

[8] 姜蕓，王軍，張莉.東北典型黑土區(qū)侵蝕溝形態(tài)及分布特征[J].農(nóng)業(yè)工程學(xué)報(bào)，2020，36（7）：157-165.

[9] 張瑞豪.延壽縣典型坡耕地侵蝕溝土壤物理性質(zhì)特征研究[D].哈爾濱：黑龍江大學(xué)，2022：8-9.

[10] 劉華征.典型黑土緩坡耕地土壤侵蝕特征解析[D].沈陽：沈陽農(nóng)業(yè)大學(xué)，2022：7-14.

[11] 安娟.東北黑土區(qū)土壤侵蝕過程機(jī)理和土壤養(yǎng)分遷移研究[D].北京：中國科學(xué)院研究生院（教育部水土保持與生態(tài)環(huán)境研究中心），2012：9-12.

[12] 李海強(qiáng).東北黑土區(qū)侵蝕小流域土壤質(zhì)量空間分異特征及影響因素研究[D].楊凌：西北農(nóng)林科技大學(xué)，2021：82-86.

[13] 張孝存.東北典型黑土區(qū)流域侵蝕—沉積對(duì)土壤質(zhì)量的影響[D].西安：陜西師范大學(xué)，2013：81-86.

[14] 趙旭.侵蝕退化黑土地力提升研究[D].楊凌：西北農(nóng)林科技大學(xué)，2021：1-2.

收稿日期： 2023-06-30

基金項(xiàng)目： 國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目子課題（2021YFD150080503）；黑龍江省自然科學(xué)基金項(xiàng)目（LH2020D018）；公路交通環(huán)境保護(hù)技術(shù)交通行業(yè)重點(diǎn)實(shí)驗(yàn)室開放課題（2014—8801）；中央級(jí)公益性科研院所基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目（2017—9058）

第一作者： 杜兆國（1971—），男，黑龍江哈爾濱人，高級(jí)工程師，學(xué)士，主要從事水土保持生態(tài)建設(shè)工作。

通信作者： 肖洋（1978—），男，黑龍江哈爾濱人，副教授，博士，主要從事農(nóng)業(yè)水土保持研究。

E-mail： xiaoyang@hlju.edu.cn

（責(zé)任編輯李佳星）