深松結合秸稈還田對黑土孔隙結構的影響

楊建君，蓋浩，張夢璇，蔡育蓉，王力艷，王立剛

深松結合秸稈還田對黑土孔隙結構的影響

1中國農業科學院農業資源與農業區劃研究所/農業農村部面源污染控制重點實驗室，北京 100081；2青岡縣氣象局，黑龍江青岡 151600

【目的】東北黑土實施深松結合秸稈還田對土壤孔隙結構影響的研究缺乏明確性結論，為此開展本研究，旨在研究深松結合秸稈還田對黑土結構的影響機制，為合理耕層創建提供科學依據。【方法】利用在東北典型黑土區——黑龍江省綏化市青岡縣開展的5年田間定位試驗為平臺，設置農民常規耕作（FP）、單獨深松25 cm（T2）、深松25 cm結合秸稈還田（T3）和深松35 cm結合秸稈還田（T4）等處理，采用CT掃描技術開展土壤孔隙結構可視化和定量化研究，并結合田間持水量和容重等指標，探究深松結合秸稈還田對黑土孔隙結構的影響。【結果】通過土壤孔隙二維和三維圖像可以清晰看出，各處理0—20 cm土層孔隙分布均明顯少于20—40 cm土層，深松結合秸稈還田（T3、T4）的孔隙分布明顯多于FP，增加了結構更為復雜的大孔隙。定量化分析表明，相較于FP，單獨深松25 cm（T2）顯著提高20—30 cm土層總孔隙度103.0%（＜0.05），主要通過顯著提高小孔隙（孔隙直徑d≤0.50 mm）孔隙度91.3%和中孔隙（0.50 mm＜d≤1.00 mm）孔隙度143.5%來實現（＜0.05）；而深松結合秸稈還田（T3、T4）則可顯著提高0—30 cm土層總孔隙度109.8%—382.7%（＜0.05），主要通過顯著提高大孔隙（d＞1.00 mm）孔隙度221.5%—661.7%和中孔隙孔隙度105.4%—544.9%來實現（＜0.05）。T3、T4還顯著提高了0—30 cm土層孔隙的分形維數9.9%—17.7%（＜0.05），降低了歐拉數32.4%—66.4%（＜0.05），顯著提高田間持水量24.2%—40.6%（＜0.05）。進一步分析得出不同孔徑孔隙度、總孔隙度與田間持水量、分形維數呈極顯著正相關，與歐拉數呈極顯著負相關（＜0.01）。【結論】深松結合秸稈還田能夠提高黑土大中孔隙的孔隙度、改善孔隙結構和連通性、增加田間持水量，尤以深松35 cm結合秸稈還田的效果最為顯著，可作為東北黑土合理耕層構建的推薦措施。

東北黑土；深松結合秸稈還田；CT掃描技術；土壤孔隙結構

0 引言

【研究意義】東北黑土區作為我國重要商品糧生產基地，自開墾以來，隨著連年高強度的開發利用，機械作業的連年壓實導致其結構變差，嚴重影響了黑土的肥力供給和產能提升，甚至威脅到我國的糧食安全[1]。因此，如何保護黑土地，提高黑土耕地質量，成為黑土地保護亟待解決的首要問題[2]【前人研究進展】眾多學者針對于黑土結構變差的問題提出了不同解決措施，包括改進耕作措施[3-4]、秸稈還田[5-6]、施用生物炭[7]等，結果證實均可以不同程度改善黑土結構，有利于黑土合理耕層的構建，其中深松和秸稈還田措施是公認為具有降低土壤容重[8]、增加耕層厚度[9]、改善孔隙分布[5]的效應，若將兩者結合，則可以充分發揮兩者的優勢，改善土壤理化性狀[10-11]，有效改良黑土結構。當前研究多通過測定土壤有機碳含量、土壤容重、貫入阻力等基礎理化性狀[6,10,12]來探究深松和秸稈還田措施對黑土的影響，鮮有從土壤孔隙結構的角度進行研究，特別是缺乏可視化和定量化相結合的方法，對其改善土壤性狀的效果難以有科學直觀的判斷與評價。土壤孔隙結構指土壤孔隙的形態大小、數量搭配及其分布狀況，包括孔隙度、孔隙半徑、孔隙大小分布等形態和數量特征，以及孔隙的相互連通狀況和孔隙之間相關性等空間分布特征，是土壤結構的重要方面[13]，孔隙結構直接影響土壤持水能力[5,14]和土壤肥力[15]。因此，對土壤孔隙結構的研究具有重要的實踐意義和指導作用。CT技術，即X射線計算機斷層掃描（X-ray computed tomography）技術[16]，可以快速獲取原狀土壤內部結構，是非破壞性檢測土壤孔隙3D結構的一種新興手段[17]，可以實現土壤孔隙結構的可視化與定量化[18-19]，與傳統基本理化性狀指標監測技術相比具有不可比擬的優勢。隨著CT掃描技術及后續分析手段的發展進步，近年來已經有不少學者利用這一方式進行土壤孔隙可視化與定量化研究。邱琛等[5]2019年利用CT掃描技術研究有機物料還田深度對黑土孔隙結構影響，結果顯示秸稈淺混還田、秸稈深混還田和秸稈配合有機肥深混還田處理對0—15 cm土層＞500 μm孔隙數量和孔隙度分別顯著增加了18.1%—179.9%和69.2%—256.0%（＜0.05），表明施用有機物料能夠增加相應土層的孔隙數量，改善孔隙分布，增加孔隙結構的復雜性，并且在土壤中出現了交叉孔隙和細長孔隙，提高了孔隙的連通性。王憲玲等[20]基于CT技術研究發現相比單施化肥，有機無機肥配施可以提高20—40 cm土層大孔隙度91.7%（＜0.05），增加大孔隙數量54.8%，土壤持水和導水性也能明顯提高。房煥等[21]基于CT技術研究發現稻麥輪作區進行全量秸稈還田能夠降低土壤容重19.4%，增大土壤大孔隙度110.6%，改善水稻土物理結構。【本研究切入點】目前基于CT掃描技術的研究多是探究免耕、有機肥替代化肥、秸稈還田與常規耕作措施之間的土壤孔隙結構差異[20-25]，而關于深松結合秸稈還田措施對黑土孔隙結構的研究鮮有報道，且缺乏系統性。【擬解決的關鍵問題】針對上述現狀和問題，本研究利用在東北典型黑土區開展的田間定位試驗，通過CT掃描技術開展土壤孔隙結構可視化與定量化的研究，探究深松結合秸稈還田措施對東北黑土土壤孔隙結構的影響，以期為深松結合秸稈還田措施對黑土物理結構的影響機制研究及合理耕層的構建提供科學依據。

1 材料與方法

1.1 試驗地概況

試驗地點設在黑龍江省綏化市青岡縣民政鎮進化村（126°11′3.14″N，46°42′42.11″E）。該地處于松嫩平原腹地，是典型東北黑土區域，溫帶大陸性季風氣候，平均海拔457 m，年平均氣溫2.4—2.6℃，全年無霜期130 d左右。試驗開展期間（2017—2021年）年降水量為614—753 mm，且降雨主要集中在6—8月，占全年降雨量60%—80%；2020年試驗地點遭遇極端臺風天氣，9月又發生多次強降雨（月總降水量達220.7 mm），全年降雨量達到753 mm。供試土壤為黏壤質黑鈣土，0—100 cm土層基礎理化性狀如表1所示。

1.2 試驗設計

試驗開始于2017年，采用大區設計，每個試驗處理面積為1 200 m2（50 m×24 m）。共設置4個處理，分別是農民常規耕作（FP，旋耕）、單獨深松25 cm（T2）、深松25 cm結合秸稈還田（T3）、深松35 cm結合秸稈還田（T4），T3、T4處理均為秸稈全量還田，試驗田及取樣點布置如圖1所示。

每年秋季玉米收獲后，深松結合秸稈還田處理（T3、T4）地塊首先用秸稈碎混機進行秸稈原位碎混（碎混至秸稈長度小于5 cm），每兩年進行一次相應深度的深松（每次深松位置基本一致），之后使用翻地機器將粉碎后的秸稈翻到20—25 cm土層，用耙地機器打碎土塊后旋耕起壟（作業機械類型及流程如圖2所示）。單獨深松處理（T2）：玉米收獲后將秸稈移出試驗田，每兩年進行一次深松，之后旋耕起壟。農民常規耕作（FP）：玉米收獲后將秸稈移出試驗田，進行旋耕起壟。所有深松作業均使用德國LEMKEN公司的Karat 9深松機進行，深松深度通過儀器數字電控；所有旋耕作業深度一致，均約15 cm。T2、T3、T4處理于2017、2019、2021年秋季分別進行相應的深松，2020年秋季由于受到臺風和土壤濕度太大的影響，未進行翻地、起壟等措施，于2021年春季進行了相應的耕作措施。

圖1 試驗田及取樣點布置示意圖

表1 2017年試驗區土壤基礎理化性狀

圖2 深松結合秸稈還田耕作流程

各處理均于5月初播種，10月初收獲。供試品種為天澤703，播種行間距66 cm，株間距25 cm，種植密度約為6萬株/hm2。

1.3 樣品采集與測定

本研究于2021年8月進行原狀土柱取樣（圖2）。目前關于CT技術的前期取樣尚未有統一的技術規范，但都是以獲取原狀土體，盡量減少取樣過程對所取原狀土體的擾動為原則。一些學者使用PVC管進行原狀土壤取樣[23,26]，本題組在2019年也嘗試利用該方法進行取樣，但效果不理想，表明在東北質地較為黏重的黑土使用PVC管取樣難度較大且取樣過程易發生較大的擾動，進而影響測定的準確性，因此本課題組采用定制的白鋼取樣盒進行剖面取樣（取樣盒長、寬、高參數為50 cm×10 cm×10 cm，50 cm×10 cm的一面可開合，如圖3所示）。

每個試驗大區取3個原狀土柱作為重復，取樣點布置如圖1所示。由于本研究采用大區試驗，耕作措施均由中型及大型機械進行，機械壓實及深松尺作業均主要作用于壟溝，故取樣點均布設在壟溝。原狀土柱取樣時，先在取樣地塊挖出長約60 cm，寬約40 cm，深約60 cm的長方體剖面，將取樣盒扣在取樣剖面（寬40 cm，深60 cm的一面，要保證此面較平滑，豎直方向垂直地面），然后在取樣盒底部墊上厚木板（減少擾動以及防止取樣盒受敲擊變形），用錘子擊打厚木板直至將取樣盒完全扣入土體，在側面用鐵鍬將裝有原狀土柱的取樣盒挖出，挖出后用刮刀將多余土壤去掉，扣上蓋子，用膠帶密封，盡量減少水分蒸發散失，使用氣柱包裝袋將裝有原狀土體的取樣盒包好，立即送到實驗室進行CT掃描。取樣及運輸過程盡量減少對原狀土柱的擠壓、碰撞等擾動。掃描樣品完成后，用體積100 cm3的環刀分別在原狀土柱0—10、10—20、20—30、30—40 cm土層取原狀土，用于田間持水量、容重指標的測定[27]，并取0—10、10—20、20—30、30—40 cm土層土樣，風干測土壤有機碳含量。

圖3 取樣盒示意圖

原狀土柱CT掃描委托英華檢測（上海）有限公司分析測試中心測定，采用GE公司的phoenix v | tome x m多功能X射線微聚焦CT系統，掃描參數：掃描電壓為220 kV，掃描電流為230 μA，分辨率為100 μm，為避免取樣及運輸對兩端的擾動，原狀土柱最上端3 cm及最下端7 cm不進行掃描，僅掃描其余40 cm原狀土柱，掃描結束每個原狀土柱獲取約4 000張16位tiff格式圖像。

1.4 圖像處理及數據獲取

圖像處理采用image J軟件進行，主要處理過程如下[17]：（1）導入圖像。為更直觀整體展現原狀土柱孔隙結構，可視化對象為整個原狀土柱。基于計算機計算能力的限制，可視化過程采取隔一選一的方式（Increment=2），每個土樣選取約2 000張灰度圖像進行處理；定量化過程分層進行處理（10 cm/層），每層選取對應約1 000張灰度圖像進行處理（Increment=1）；（2）圖像增強。調節圖像亮度、對比度；（3）圖像去噪。通過中值濾波平滑圖像；（4）選擇分析區域。為了避免邊界部分的影響，利用方形工具選取未受擾動區域進行分析，本研究中可視化與定量化均選取200像素×200像素的分析區域，換算實際土樣切面面積大小為2 cm×2 cm；（5）二值分割。目視法確定圖像的分割閾值，對圖像進行二值化處理；（6）孔隙可視化及定量化。對于分割后的二值圖像，進行顆粒分析獲取每個土層的孔隙面積和孔隙度，利用Bone J插件獲取孔隙大小分布、分形維數及歐拉數，孔隙的可視化通過3D viewer插件實現。

1.5 數據分析

采用Microsoft Excel 2019對試驗數據進行整理，采用SPSS 21.0進行單因素方差分析，分析比較深松結合秸稈還田措施對土壤孔隙結構參數及物理性質的影響，并用Duncan法進行多重比較（α=0.05）。采用Origin 2021作圖。采用SPSS 21.0 軟件進行土壤物理性質與土壤結構參數之間的Pearson相關性分析。

2 結果

2.1 不同深松結合秸稈還田對土壤孔隙結構的影響

2.1.1 土壤孔隙二維和三維結構可視化特征 不同處理土壤孔隙二維結構示意如圖4所示，農民常規耕作處理（FP）整體大孔隙較少，更多是較小的孔隙，且較為分散，相互之間缺乏連通性，而經過深松的T2、T3、T4處理二維圖像都直觀地反映出土壤大、中孔隙數量有不同程度的增加，且深松結合秸稈還田處理（T3、T4）孔隙改良效果更明顯，在二維圖像中體現為土壤大孔隙變多，孔隙性狀變復雜，孔隙和孔隙之間的連通也清晰可見。

圖中白色部分為孔隙，黑色部分為土壤基質，實際尺寸為2 cm×2 cm

從不同處理土壤孔隙三維結構來看，0—20 cm土層孔隙數量及分布明顯少于20—40 cm土層，這在未進行秸稈還田的FP處理和單獨深松25 cm的T2處理表現尤為突出（圖5），證實了東北黑土存在耕層結構變差的問題。T2處理0—20 cm土層孔隙大小分布及復雜程度相較于FP處理差異較小，但在20—40 cm土層，T2處理土壤孔隙復雜程度有所提高。而T3和T4處理對比FP處理，在0—30 cm土層（T3、T4），甚至30—40 cm土層（T4），孔徑較大的孔隙都明顯增多，且結構看起來更為復雜，分布更為密集，表明深松結合秸稈還田措施對其作業深度及以上土層孔隙作用效果明顯，有效增加了大孔隙分布，提高了孔隙結構的復雜程度。

圖中長方體內紅色部分為孔隙，其余部分為土壤基質, 實際尺寸為2 cm×2 cm×40 cm

2.1.2 土壤不同直徑孔隙分布特征 從總孔隙度分析表明，相較于FP處理，T3和T4處理對0—30 cm土層總孔隙度產生較為顯著的影響，分別顯著提高128.4%—382.7%和109.8%—339.3%（＜0.05），T4處理對30—40 cm土層孔隙也產生較大影響，較FP處理顯著提高總孔隙度152.3%（＜0.05，圖6）。單獨深松25 cm（T2）僅對深松作業深度土層（20—30 cm）產生顯著影響，較FP處理顯著提高總孔隙度103.0%（＜0.05）。從中可以看出，深松結合秸稈還田措施可以有效提高其作業深度土層孔隙度，為土體適宜的“三相比”的形成奠定了土壤結構基礎[28-29]。

圖6 不同深松結合秸稈還田對不同深度土壤總孔隙度的影響

不同處理在不同土層對大、中、小孔隙的孔隙度具有不同的影響效果（圖7）。在0—30 cm土層，與農民常規耕作（FP）相比，深松結合秸稈還田處理（T3、T4）都能顯著提高中孔隙（0.50 mm＜d≤1.00 mm）和大孔隙（d＞1.00 mm）的孔隙度，但T3和T4處理之間差異不顯著；在30—40 cm土層，只有T4較FP處理顯著提高小孔隙（d≤0.50 mm）和中孔隙的孔隙度，T3處理作用效果則不顯著。T2處理主要作用效果體現在20—30 cm土層的小孔隙和中孔隙，且效果弱于T3和T4處理。在0—10 cm土層，T3、T4較FP處理分別顯著提高小孔隙、中孔隙和大孔隙孔隙度231.3%、544.9%、661.7%和246.5%、406.2%、586.7%（＜0.05）；在10—20 cm土層，T3、T4處理分別較FP處理顯著提高中孔隙和大孔隙孔隙度161.4%、571.6%和138.8%、523.3%（＜0.05），小孔隙孔隙度有提高趨勢，但差異不顯著（＞0.05）；在20—30 cm土層，相較于FP處理，T3、T4處理分別顯著提高小孔隙、中孔隙和大孔隙孔隙度143.1%、155.2%、287.1%和158.2%、105.4%、221.5%（＜0.05）。T4處理在30—40 cm土層主要作用效果體現在小孔隙和中孔隙孔隙度，較FP處理分別顯著提高63.0%和91.3%（＜0.05）。T2處理主要作用于10—20 cm土層的中孔隙孔隙度，以及20—30 cm土層小孔隙和中孔隙孔隙度，較FP處理分別顯著提高151.6%、91.3%和143.5%（＜0.05）。

圖7 不同深松結合秸稈還田對不同深度土層各孔徑孔隙度的分布

2.1.3 土壤孔隙分形維數和歐拉數 孔隙分形維數綜合反映了土壤孔隙的大小與形狀，其值越大，孔隙的狀況越為優越[5,30]。與農民常規耕作（FP）相比，深松結合秸稈還田處理（T3、T4）土壤孔隙分形維數有顯著提高，且主要體現在機械作業深度及以上的土層（0—30 cm），而單獨深松25 cm處理T2與FP處理相比雖然有提高趨勢，但并未表現出顯著性差異（表2）。在0—30 cm土層，T3、T4處理的分形維數較FP顯著提高9.9%—17.7%（＜0.05），但兩處理間未表現出顯著差異。在30—40 cm土層，僅有T4處理的分形維數顯著高于FP處理9.5%（＜0.05）。

表2 不同深松結合秸稈還田對土壤孔隙分形維數和歐拉數的影響

歐拉數可以用來表征土壤孔隙的連通性，其值越低，表示孔隙連通性越好[31]。本研究中0—40 cm土層范圍內，相較于FP處理，T2、T3、T4處理的歐拉數均顯著降低（表2），其中以T3處理降低幅度最大，降低了40.1%—66.4%（＜0.05）；T4處理次之，降低了32.4%—56.3%（＜0.05），T2處理降低了32.2%—38.3%（＜0.05）。這也表明深松結合秸稈還田不僅提高了相應土層孔隙度的大中孔隙孔隙度，還能改善孔隙的結構特征，提高孔隙連通性，有利于土壤水分的運移和保存[17]。

2.2 不同深松結合秸稈還田對土壤物理性狀的影響

對比FP處理，T2、T3、T4處理均有降低土壤容重的趨勢，但均尚未達到顯著水平（圖8）。T2、T3、T4處理均不同程度提高了土壤田間持水量，其中T3和T4處理效果最為顯著。在0—30 cm土層，相較于FP處理，T3、T4田間持水量分別顯著提高了20.7%—40.6%和24.2%—35.0%（＜0.05），T2處理較FP處理顯著提高了0—10和10—20 cm土層田間持水量18.9%和26.3%（＜0.05）。在30—40 cm土層深度，T2、T3、T4處理較FP處理均有提高趨勢，但未達到顯著水平（＞0.05）。

2.3 土壤孔隙結構與物理性狀的相關性

對各土壤孔隙結構參數和土壤物理性狀參數進行相關分析表明，不同孔徑孔隙度、總孔隙度均與歐拉數呈現極顯著負相關（＜0.01），與分形維數呈現極顯著正相關（＜0.01）（圖9）。小孔隙、中孔隙、大孔隙孔隙度、總孔隙度、分形維數均與土壤容重呈現極顯著負相關（＜0.01），與田間持水量呈現極顯著正相關（＜0.01），即土壤孔隙分布及結構會影響到土壤容重及田間持水量。

3 討論

3.1 深松結合秸稈還田措施下黑土孔隙結構的可視化

CT掃描技術可以實現對土壤孔隙結構的可視化[18]，以直觀了解土壤耕層孔隙狀況。本研究中，從二維和三維圖像都可以清晰看出，在0—20 cm土層，各處理均不同程度出現孔隙分布相對較少的情況，而20—40 cm土層孔隙分布更多，這與邱琛等[5]在黑土研究中0—15 cm土層孔隙度明顯多于15—30 cm土層孔隙度的結果有所不同，其原因在于本研究采用大區試驗，各耕作管理措施均通過機械進行，符合目前黑土農機耕種的生產實際情況，同時取樣點布設在壟溝，更充分體現了機械對土壤一定程度壓實的作用[32]，即機械作業對上層土壤直接的碾壓作用導致了上層土壤孔隙，特別是大孔隙明顯減少[33]，呈現出上層土壤孔隙少于下層土壤孔隙的現象。而邱琛等開展的研究主要是通過傳統的小區試驗來體現有機物料還田深度對孔隙結構的影響，未充分體現農機作業的壓實影響。本研究中，經過5年的試驗，單獨深松處理（T2）在0—20 cm土層孔隙分布與FP處理相比未呈現出明顯差異，而深松結合秸稈還田處理（T3、T4）在0—20 cm出現一些較大且結構復雜的孔隙，與FP處理相比呈現出較為明顯的差異，表明深松結合秸稈還田可以在一定程度上緩解機械碾壓造成上層土壤孔隙減少、土壤結構變差的問題，這與房煥等[21]得出的秸稈還田顯著提高表層土壤孔隙度及孔隙復雜程度的研究結果基本一致。

圖8 不同深松結合秸稈還田對土壤容重和田間持水量的影響

圖中SP、MP、LP、TP、FD、EN、BD、FC分別代表小孔隙孔隙度、中孔隙孔隙度、大孔隙孔隙度、總孔隙度、分形維數、歐拉數、容重、田間持水量。**表示在0.01水平上極顯著

關于原狀土體取樣，目前尚未有針對原狀土體取樣范圍及方式建立統一的規范，學者們多根據試驗情況采取適合自己研究的方法，例如PVC管[5,21,26]等。本課題組采用定制的白鋼取樣盒進行剖面取樣，整體效果有所改善，但仍存在原狀土柱邊緣受到擾動破損的狀況，故本研究在最終處理時僅保留中心2 cm×2 cm×40 cm的土柱，最大限度減少了取樣擾動的影響，保證了原狀土柱樣品分析的準確性與科學性，當然，若能在取樣方式上進一步完善改進，實現更好的原狀土體取樣效果，將有力促進CT掃描技術在土壤孔隙結構研究上的應用，有利于更好呈現土壤結構特征。

3.2 深松結合秸稈還田措施下黑土孔隙結構的定量化

本研究中，單獨深松25 cm處理（T2）通過顯著提高了作業深度（20—30 cm）土層小孔隙和中孔隙孔隙度來提高了總孔隙度，這與李嵩等[9]提出的深松提高犁底層通氣孔隙度和總孔隙度的結論基本一致。深松在不翻轉土壤的情況下對相應土層進行橫向切割，活化耕層土壤結構，增加土壤孔隙度[34]，但深松后土壤結構不穩定，土壤強度較低，易被再次壓實[35]，因此本研究中單獨深松25 cm處理（T2）作用效果主要體現在作業深度的小孔隙和中孔隙，并未對大孔隙產生顯著影響，這一結論在深松35 cm結合秸稈還田處理（T4）的30—40 cm土層也有體現（圖7）。但兩者在孔隙結構參數方面有所不同，T4處理深松作業深度30—40 cm土層孔隙分形維數較FP處理顯著提高，而T2處理深松作業深度20—30 cm土層分形維數較FP處理并未呈現出顯著差異，這可能是機械碾壓對上層土壤壓實作用效果強于下層土壤[36]，從而導致上層土壤深松效果的弱化。

深松結合秸稈還田作業處理（T3、T4）不僅有效提高深松作業深度的小孔隙和中孔隙孔隙度，還顯著提高了秸稈還田作業深度及以上耕層（0—30 cm）的大孔隙孔隙度，其孔隙結構參數表現也優于單獨深松25 cm處理（T2），這與叢聰等[10]2021年在黑土區研究得出的深松配施有機物料顯著增加土壤孔隙度的結論基本一致。一方面還田的秸稈在微生物和酶的共同作用下向土壤提供了大量小于土壤密度的有機物質，有助于土壤有機碳的固存，并與土壤顆粒結合形成穩定疏松的結構，增大了土壤孔隙度[37-38]，促進原有孔隙向大孔隙轉化[39]，這在本研究中也有所體現，T3、T4處理較FP處理顯著提高了0—40 cm土層土壤有機碳含量6.6%—27.5%（＜0.05）（表3）。另一方面，秸稈還田還有利于農田土壤生物的生存環境的維護[40]，顯著增加土壤動物優勢類群[41]，而土壤動物在土壤中的活動促進了土壤孔隙的形成，改善土壤結構[42]，這可能是T3、T4處理大孔隙孔隙度顯著提高的另一個原因。此外，T3、T4處理是通過翻耕將碎混的秸稈還到20—25 cm深度土層，每年秸稈還田都會伴隨著土層翻轉與混合，因此經過5年的深松結合秸稈還田作業后，秸稈還田對大孔隙的作用效果可體現在整個0—30 cm土層，而不僅僅是秸稈還田深度土層，加之深松及翻耕過程對土壤的擾動，進一步促進了孔隙的形成，改善了孔隙狀況[28,43-44]，即深松結合秸稈還田有效發揮了深松和秸稈還田兩措施優勢，具有改良土壤孔隙結構的效果，促進黑土結構更加趨向于合理化，這就給農作物的生長發育創造了良好的土壤微環境。

表3 各處理0—40 cm土層土壤有機碳含量

3.3 深松結合秸稈還田對黑土物理性狀的影響

本研究中，T2、T3、T4處理在不同深度耕層容重較FP處理均有降低趨勢，但未達到顯著水平，這可能是機械的壓實作用削弱了對容重的改良效果[45]，或者是5年的作用年限相對較短，仍未充分發揮出對容重的作用效果，還需要更長年限試驗來體現[46]。相較于FP處理，T3、T4處理在0—30 cm土層田間持水量均顯著提高，結合前文定量化與可視化結果中T3、T4處理大中孔隙顯著增加的結果，即土壤中混入打碎的秸稈起到“楔子”的作用，配合深松及翻耕對土壤的擾動，改善了土壤的物理結構，優化了土壤“三相比”[47]，增加了可以儲存水分的孔隙[48]，從而提高了土壤的持水保水性能，這與鞠忻倪[49]等提出的土壤孔隙的數量、大小和空間結構對土壤的持水保水性能有著重要影響的結論相吻合；另一方面，還田的秸稈轉換為土壤中的有機質，具有巨大的比表面積、親水基團和較強的吸水性，提高了土壤毛管孔隙度，增強了土壤的持水能力[1,50]。對土壤孔隙結構參數與土壤物理性狀進行相關分析顯示，土壤孔隙結構參數和土壤容重、田間持水量具有極顯著的相關關系，即深松結合秸稈還田措施通過提高不同孔徑的孔隙度，改善孔隙連接性及孔隙狀況，從而降低了土壤容重，增加了田間持水量，以實現土壤物理性狀的改良，這與阮仁杰[51]研究結論基本一致。

本研究可視化結果明確了不同耕層的孔隙分布及目視特征，定量化結果實現了對孔徑分布及孔隙結構的定量化，后續的研究中，若能將可視化、定量化結果與更多物理化學及生物學性狀相結合，細化分析各土壤孔隙特征對土壤理化性狀等影響機制及貢獻程度，將極大有助于我們深入理解耕作對土壤微生態的影響。

4 結論

4.1 從二維和三維圖像中可以明顯顯示出東北黑土農民常規耕作（FP）受到機械作業壓實的影響，0—20 cm土層孔隙分布明顯少于20—40 cm土層，深松結合秸稈還田處理（T3、T4）可以在一定程度改善這一現象，且增加了結構更為復雜的大孔隙。

4.2 相較于FP處理，單獨深松25 cm處理（T2）顯著提高20—30 cm土層總孔隙度103.0%（＜0.05），主要通過顯著提高小孔隙（d≤0.50 mm）孔隙度91.3%和中孔隙（0.50 mm＜d≤1.00 mm）孔隙度143.5%實現（＜0.05）；而深松結合秸稈還田處理（T3、T4）則可顯著提高0—30 cm土層的總孔隙度109.8%—382.7%（＜0.05），主要通過顯著提高大孔隙（d＞1.00 mm）孔隙度221.5%—661.7%和中孔隙孔隙度105.4%—544.9%實現（＜0.05）。

4.3 深松結合秸稈還田處理（T3、T4）較FP處理顯著提高了0—30 cm土層孔隙的分形維數9.9%—17.7%（＜0.05），降低了歐拉數32.4%—66.4%（＜0.05），顯著提高了田間持水量24.2%—40.6%（＜0.05），在改善孔隙狀況，提高了孔隙連通性的同時，也有增強土壤持水能力的體現。

4.4 相關分析顯示，田間持水量和不同孔徑孔隙度、總孔隙度、分形維數呈極顯著正相關，與歐拉數呈極顯著負相關，容重和不同孔徑孔隙度、總孔隙度、分形維數呈極顯著負相關，與歐拉數呈極顯著正相關（＜0.01），即田間持水量和容重受到土壤孔隙直徑分布及孔隙結構的影響。

綜上所述，深松結合秸稈還田能夠提高土壤大中孔徑的孔隙度、改善孔隙結構和連通性、增加田間持水量，可作為東北黑土合理耕層構建的推薦技術。

[1] 郭孟潔, 李建業, 李健宇, 齊佳睿, 張興義. 實施16年保護性耕作下黑土土壤結構功能變化特征. 農業工程學報, 2021, 37(22): 108-118.

GUO M J, LI J Y, LI J Y, QI J R, ZHANG X Y. Changes of soil structure and function after 16-year conservation tillage in black soil. Transactions of the Chinese Society of Agricultural Engineering, 2021, 37(22): 108-118. (in Chinese)

[2] 韓曉增, 鄒文秀. 我國東北黑土地保護與肥力提升的成效與建議. 中國科學院院刊, 2018, 33(2): 206-212. doi:10.16418/j.issn.1000- 3045.2018.02.011.

HAN X Z, ZOU W X. Effects and suggestions of black soil protection and soil fertility increase in northeast China. Bulletin of Chinese Academy of Sciences, 2018, 33(2): 206-212. doi:10.16418/j.issn. 1000-3045.2018.02.011. (in Chinese)

[3] 張博文, 楊彥明, 張興隆, 李金龍, 陳新宇, 李志新. 連續深松對黑土結構特性和有機碳及碳庫指數影響. 中國土壤與肥料, 2019(2): 6-13.

ZHANG B W, YANG Y M, ZHANG X L, LI J L, CHEN X Y, LI Z X. Effects of continuous deep loosening on soil physical characteristics, organic carbon content and carbon pool index in black soil. Soil and Fertilizer Sciences in China, 2019(2): 6-13. (in Chinese)

[4] 李景, 吳會軍, 武雪萍, 王碧勝, 姚宇卿, 呂軍杰. 長期免耕和深松提高了土壤團聚體顆粒態有機碳及全氮含量. 中國農業科學, 2021, 54(2): 334-344.

LI J, WU H J, WU X P, WANG B S, YAO Y Q, Lü J J. Long-term conservation tillage enhanced organic carbon and nitrogen contents of particulate organic matter in soil aggregates. Scientia Agricultura Sinica, 2021, 54(2): 334-344. (in Chinese)

[5] 邱琛, 韓曉增, 陳旭, 陸欣春, 嚴君, 馮玉鈿, 甘佳偉, 鄒文秀, 劉國輝. CT掃描技術研究有機物料還田深度對黑土孔隙結構影響.農業工程學報, 2021, 37(14): 98-107.

QIU C, HAN X Z, CHEN X, LU X C, YAN J, FENG Y T, GAN J W, ZOU W X, LIU G H. Effects of organic amendment depths on black soil pore structure using CT scanning technology. Transactions of the Chinese Society of Agricultural Engineering, 2021, 37(14): 98-107. (in Chinese)

[6] 高洪軍, 彭暢, 張秀芝, 李強, 朱平, 王立春. 秸稈還田量對黑土區土壤及團聚體有機碳變化特征和固碳效率的影響. 中國農業科學, 2020, 53(22): 4613-4622.

GAO H J, PENG C, ZHANG X Z, LI Q, ZHU P, WANG L C. Effects of corn straw returning amounts on carbon sequestration efficiency and organic carbon change of soil and aggregate in the black soil area. Scientia Agricultura Sinica, 2020, 53(22): 4613-4622. (in Chinese)

[7] 張偉明, 陳溫福, 孟軍, 金梁, 郭偉, 趙洪亮. 東北地區秸稈生物炭利用潛力、產業模式及發展戰略研究. 中國農業科學, 2019, 52(14): 2406-2424.

ZHANG W M, CHEN W F, MENG J, JIN L, GUO W, ZHAO H L. Study of straw-biochar on utilization potential, industry model and developing strategy in northeast China. Scientia Agricultura Sinica, 2019, 52(14): 2406-2424. (in Chinese)

[8] 翟振, 李玉義, 郭建軍, 王婧, 董國豪, 郭智慧, 逄煥成. 耕深對土壤物理性質及小麥-玉米產量的影響. 農業工程學報, 2017, 33(11): 115-123.

ZHAI Z, LI Y Y, GUO J J, WANG J, DONG G H, GUO Z H, PANG H C. Effect of tillage depth on soil physical properties and yield of winter wheat-summer maize. Transactions of the Chinese Society of Agricultural Engineering, 2017, 33(11): 115-123. (in Chinese)

[9] 李嵩, 韓巍, 張凱, 依艷麗. 不同耕作方式對遼西褐土物理性狀及玉米根系分布的影響. 玉米科學, 2020, 28(6): 101-106. doi:10. 13597/j.cnki.maize.science.20200615.

LI S, HAN W, ZHANG K, YI Y L. Effect of different tillage methods on soil structure and maize root distribution in cinnamon soil area in western Liaoning. Journal of Maize Sciences, 2020, 28(6): 101-106. doi:10.13597/j.cnki.maize.science.20200615. (in Chinese)

[10] 叢聰, 王天舒, 岳龍凱, 周璇, 李玉明, 堯水紅. 深松配施有機物料還田對黑土區坡耕地土壤物理性質的改良效應. 中國土壤與肥料, 2021(3): 227-236.

CONG C, WANG T S, YUE L K, ZHOU X, LI Y M, YAO S H. Amendment effect of subsoiling with organic materials application on soil physical properties of slope cropland in mollisol region. Soil and Fertilizer Sciences in China, 2021(3): 227-236. (in Chinese)

[11] 劉平奇. 不同耕作措施下東北黑土平/坡耕地有機碳平衡及影響因素研究[D]. 北京: 中國農業科學院, 2020.

LIU P Q. Study on organic carbon balance and influencing factors of flat/sloping farmlands with different tillage measures in black soil region of northeast China[D]. Beijing: Chinese Academy of Agricultural Sciences, 2020. (in Chinese)

[12] 賀美, 王迎春, 王立剛, 李成全, 王利民, 李玉紅, 劉平奇. 深松施肥對黑土活性有機碳氮組分及酶活性的影響. 土壤學報, 2020, 57(2): 446-456.

HE M, WANG Y C, WANG L G, LI C Q, WANG L M, LI Y H, LIU P Q. Effects of subsoiling combined with fertilization on the fractions of soil active organic carbon and soil active nitrogen, and enzyme activities in black soil in northeast China. Acta Pedologica Sinica, 2020, 57(2): 446-456. (in Chinese)

[13] 程亞南, 劉建立, 張佳寶. 土壤孔隙結構定量化研究進展. 土壤通報, 2012, 43(4): 988-994. doi:10.19336/j.cnki.trtb.2012.04.038.

CHENG Y N, LIU J L, ZHANG J B. Advance in the study on quantification of soil pore structure. Chinese Journal of Soil Science, 2012, 43(4): 988-994. doi:10.19336/j.cnki.trtb.2012.04.038. (in Chinese)

[14] POESEN J, INGELMO-SANCHEZ F. Runoff and sediment yield from topsoils with different porosity as affected by rock fragment cover and position. CATENA, 1992, 19(5): 451-474. doi:10.1016/ 0341-8162(92)90044-C.

[15] EDWARDS W M, NORTON L D, REDMOND C E. Characterizing macropores that affect infiltration into nontilled soil. Soil Science Society of America Journal, 1988, 52(2): 483-487. doi:10.2136/ sssaj1988.03615995005200020033x.

[16] 周虎, 李文昭, 張中彬, 彭新華. 利用X射線CT研究多尺度土壤結構. 土壤學報, 2013, 50(6): 1226-1230.

ZHOU H, LI W Z, ZHANG Z B, PENG X H. Characterization of multi-scale soil structure with X-ray computed tomography. Acta Pedologica Sinica, 2013, 50(6): 1226-1230. (in Chinese)

[17] 張宏媛, 逄煥成, 盧闖, 劉娜, 張曉麗, 李玉義. CT掃描分析秸稈隔層孔隙特征及其對土壤水入滲的影響. 農業工程學報, 2019, 35(6): 114-122.

ZHANG H Y, PANG H C, LU C, LIU N, ZHANG X L, LI Y Y. Pore characteristics of straw interlayer based on computed tomography images and its influence on soil water infiltration. Transactions of the Chinese Society of Agricultural Engineering, 2019, 35(6): 114-122. (in Chinese)

[18] ELLIOT T R, HECK R J. A comparison of optical and X-ray CT technique for void analysis in soil thin section. Geoderma, 2007, 141(1/2): 60-70. doi:10.1016/j.geoderma.2007.05.001.

[19] MUNKHOLM L J, HECK R J, DEEN B. Soil pore characteristics assessed from X-ray micro-CT derived images and correlations to soil friability. Geoderma, 2012, 181/182: 22-29. doi:10.1016/j.geoderma. 2012.02.024.

[20] 王憲玲, 趙志遠, 馬艷婷, 鄭朝霞, 鄭偉, 翟丙年, 趙政陽. 基于CT掃描技術研究有機無機肥長期配施對土壤物理特征的影響. 植物營養與肥料學報, 2020, 26(9): 1647-1655.

WANG X L, ZHAO Z Y, MA Y T, ZHENG Z X, ZHENG W, ZHAI B N, ZHAO Z Y. Study on the effects of long-term application of chemical fertilizer combined with manure on soil physical properties of apple orchard based on CT scanning technology. Journal of Plant Nutrition and Fertilizers, 2020, 26(9): 1647-1655. (in Chinese)

[21] 房煥, 李奕, 周虎, 顏曉元, 彭新華. 稻麥輪作區秸稈還田對水稻土結構的影響. 農業機械學報, 2018, 49(4): 297-302.

FANG H, LI Y, ZHOU H, YAN X Y, PENG X H. Effects of straw incorporation on paddy soil structure in rice-wheat rotation system. Transactions of the Chinese Society for Agricultural Machinery, 2018, 49(4): 297-302. (in Chinese)

[22] GARBOUT A, MUNKHOLM L J, HANSEN S B. Tillage effects on topsoil structural quality assessed using X-ray CT, soil cores and visual soil evaluation. Soil and Tillage Research, 2013, 128: 104-109. doi:10.1016/j.still.2012.11.003.

[23] 蔡太義, 張佳寶, 張叢志, 黃會娟, 白玉紅, 趙占輝, 李太魁. 基于顯微CT研究施肥方式對砂姜黑土大孔隙結構的影響. 干旱區資源與環境, 2017, 31(12): 143-149. doi:10.13448/j.cnki.jalre.2017.393.

CAI T Y, ZHANG J B, ZHANG C Z, HUANG H J, BAI Y H, ZHAO Z H, LI T K. Effects of fertilization mode on macropore characteristics of the Shajiang black soil (vertisol) based on computed tomography (CT) images. Journal of Arid Land Resources and Environment, 2017, 31(12): 143-149. doi:10.13448/j.cnki.jalre. 2017.393. (in Chinese)

[24] 張琪. 長期免耕措施下土壤結構定量化研究[D]. 北京: 中國地質大學(北京), 2020.

ZHANG Q. Quantitative description of soil structure under long-term no-tillage farming[D]. Beijing: China University of Geosciences, 2020. (in Chinese)

[25] GUO Y F, FAN R Q, ZHANG X P, ZHANG Y, WU D H, MCLAUGHLIN N, ZHANG S X, CHEN X W, JIA S X, LIANG A Z. Tillage-induced effects on SOC through changes in aggregate stability and soil pore structure. Science of the Total Environment, 2020, 703: 134617. doi:10.1016/j.scitotenv.2019.134617.

[26] 甘磊, 張靜舉, 黃太慶, 陳曉冰, 馬蕊, 張金蓮, 陳廷速. 基于CT技術的甘蔗地不同耕作措施下土壤孔隙結構研究. 西南農業學報, 2017, 30(8): 1843-1848. doi:10.16213/j.cnki.scjas.2017.8.025.

GAN L, ZHANG J J, HUANG T Q, CHEN X B, MA R, ZHANG J L, CHEN T S. Pore structure in sugarcane soil under different tillage managements based on CT scanning. Southwest China Journal of Agricultural Sciences, 2017, 30(8): 1843-1848. doi:10.16213/j.cnki. scjas.2017.8.025. (in Chinese)

[27] 許靜. 不同試驗方法測定田間持水量的對比研究[D]. 長春: 吉林大學, 2018.

XU J. Comparative study on the field capacity based on different test methods[D]. Changchun: Jilin University, 2018. (in Chinese)

[28] 鄒文秀, 韓曉增, 嚴君, 陳旭, 陸欣春, 邱琛, 郝翔翔. 耕翻和秸稈還田深度對東北黑土物理性質的影響. 農業工程學報, 2020, 36(15): 9-18.

ZOU W X, HAN X Z, YAN J, CHEN X, LU X C, QIU C, HAO X X. Effects of incorporation depth of tillage and straw returning on soil physical properties of black soil in Northeast China. Transactions of the Chinese Society of Agricultural Engineering, 2020, 36(15): 9-18. (in Chinese)

[29] 鄒洪濤, 王勝楠, 閆洪亮, 馬迎波, 范慶鋒, 黃毅, 張玉龍. 秸稈深還田對東北半干旱區土壤結構及水分特征影響. 干旱地區農業研究, 2014, 32(2): 52-60.

ZOU H T, WANG S N, YAN H L, MA Y B, FAN Q F, HUANG Y, ZHANG Y L. Effects of straw deep returning on soil structure moisturein semiarid region of Northeast China. Agricultural Research in the Arid Areas, 2014, 32(2): 52-60. (in Chinese)

[30] 馮杰, 于紀玉. 利用CT掃描技術確定土壤大孔隙分形維數. 灌溉排水學報, 2005, 24(4): 26-28, 40. doi:10.13522/j.cnki.ggps.2005.04. 007.

FENG J, YU J Y. Determination fractal dimension of soil macropore using comput ed tomography. Journal of Irrigation and Drainage, 2005, 24(4): 26-28, 40. doi:10.13522/j.cnki.ggps.2005.04.007. (in Chinese)

[31] ZHANG Z B, LIU K L, ZHOU H, LIN H, LI D M, PENG X H. Linking saturated hydraulic conductivity and air permeability to the characteristics of biopores derived from X-ray computed tomography. Journal of Hydrology, 2019, 571: 1-10. doi:10.1016/j.jhydrol.2019. 01.041.

[32] GRONLE A, LUX G, B?HM H, SCHMIDTKE K, WILD M, DEMMEL M, BRANDHUBER R, WILBOIS K P, HE? J. Effect of ploughing depth and mechanical soil loading on soil physical properties, weed infestation, yield performance and grain quality in sole and intercrops of pea and oat in organic farming. Soil and Tillage Research, 2015, 148: 59-73. doi:10.1016/j.still.2014.12. 004.

[33] 周艷麗, 盧秉福. 農田機械壓實對土壤物理特性的影響. 中國農機化學報, 2018, 39(9): 66-70. doi:10.13733/j.jcam.issn.2095-5553.2018. 09.015.

ZHOU Y L, LU B F. Influence on soil physical properties for farmland mechanical compaction. Journal of Chinese Agricultural Mechanization, 2018, 39(9): 66-70. doi:10.13733/j.jcam.issn.2095- 5553.2018.09.015. (in Chinese)

[34] ZHANG Y J, WANG R, WANG S L, WANG H, XU Z G, JIA G C, WANG X L, LI J. Effects of different sub-soiling frequencies incorporated into no-tillage systems on soil properties and crop yield in dryland wheat-maize rotation system. Field Crops Research, 2017, 209: 151-158. doi:10.1016/j.fcr.2017.05.002.

[35] LIU Z D, QIN A Z, ZHAO B, ATA-UL-KARIM S T, XIAO J F, SUN J S, NING D F, LIU Z G, NAN J Q, DUAN A W. Yield response of spring maize to inter-row subsoiling and soil water deficit in Northern China. PLoS ONE, 2016, 11(4): e0153809. doi:10.1371/journal.pone. 0153809.

[36] 齊紅志, 余天雨, 劉天學. 農田機械碾壓對土壤物理特性及玉米生長和產量的影響. 南方農業學報, 2021, 52(4): 951-958.

QI H Z, YU T Y, LIU T X. Effects of mechanical compaction on soil physical properties, maize growth and yield in farmland. Journal of Southern Agriculture, 2021, 52(4): 951-958. (in Chinese)

[37] 張聰, 慕平, 尚建明. 長期持續秸稈還田對土壤理化特性、酶活性和產量性狀的影響. 水土保持研究, 2018, 25(1): 92-98. doi:10. 13869/j.cnki.rswc.2018.01.016.

ZHANG C, MU P, SHANG J M. Effects of continuous returning corn straw on soil chemical properties, enzyme activities and yield trait. Research of Soil and Water Conservation, 2018, 25(1): 92-98. doi:10. 13869/j.cnki.rswc.2018.01.016. (in Chinese)

[38] FERRO N D, CHARRIER P, MORARI F. Dual-scale micro-CT assessment of soil structure in a long-term fertilization experiment. Geoderma, 2013, 204/205: 84-93. doi:10.1016/j.geoderma.2013.04. 012.

[39] 丁奠元, 馮浩, 趙英, 杜璇. 氨化秸稈還田對土壤孔隙結構的影響. 植物營養與肥料學報, 2016, 22(3): 650-658.

DING D Y, FENG H, ZHAO Y, DU X. Effect of ammoniated straw returning on soil pore structure. Journal of Plant Nutrition and Fertilizers, 2016, 22(3): 650-658. (in Chinese)

[40] 楊旭, 高梅香, 張雪萍, 林琳, 沙迪, 張利敏. 秸稈還田對耕作黑土中小型土壤動物群落的影響. 生態學報, 2017, 37(7): 2206-2216.

YANG X, GAO M X, ZHANG X P, LIN L, SHA D, ZHANG L M. Effect of straw-returning management on meso-micro soil fauna in a cultivated black soil area. Acta Ecologica Sinica, 2017, 37(7): 2206-2216. (in Chinese)

[41] 朱強根, 朱安寧, 張佳寶, 張煥朝, 楊淑莉, 王意錕. 保護性耕作下土壤動物群落及其與土壤肥力的關系. 農業工程學報, 2010, 26(2): 70-76.

ZHU Q G, ZHU A N, ZHANG J B, ZHANG H C, YANG S L, WANG Y K. Relation of agricultural soil fauna and soil fertility under conservation tillage systems. Transactions of the Chinese Society of Agricultural Engineering, 2010, 26(2): 70-76. (in Chinese)

[42] 劉鵬飛. 玉米秸稈還田量對黑土區農田土壤動物群落的影響[D]. 呼和浩特: 內蒙古農業大學, 2019.

LIU P F. Effect of quantity of returned corn straw on the cropland soil fauna community in a black soil area[D]. Hohhot: Inner Mongolia Agricultural University, 2019. (in Chinese)

[43] 徐永剛, 馬強, 周樺, 姜春明, 宇萬太. 秸稈還田與深松對土壤理化性狀和玉米產量的影響. 土壤通報, 2015, 46(2): 428-432. doi:10.19336/j.cnki.trtb.2015.02.026.

XU Y G, MA Q, ZHOU H, JIANG C M, YU W T. Effects of straw returning and deep loosening on soil physical and chemical properties and maize yields. Chinese Journal of Soil Science, 2015, 46(2): 428-432. doi:10.19336/j.cnki.trtb.2015.02.026. (in Chinese)

[44] XU D, MERMOUD A. Topsoil properties as affected by tillage practices in North China. Soil and Tillage Research, 2001, 60(1/2): 11-19. doi:10.1016/S0167-1987(01)00167-2.

[45] GAN L, PENG X, PETH S, HORN R. Effects of grazing intensity on soil thermal properties and heat flux underandvegetation in Inner Mongolia, China. Soil and Tillage Research, 2012, 118: 147-158. doi:10.1016/j.still.2011.11.005.

[46] 崔正果. 不同年限玉米秸稈還田對黑土土壤理化性狀以及土壤微生物的影響[D]. 長春: 吉林大學, 2019.

CUI Z G. Effect of different years straw returning on physicochemical and microorganism characteristics of black soil[D]. Changchun: Jilin University, 2019. (in Chinese)

[47] 周彥莉, 吳海梅, 周彥棟, 尚旭民, 逄蕾. 短期秸稈不同還田方式對土壤結構和水分影響. 干旱區研究, 2022, 39(2): 502-509. doi:10.13866/j.azr.2022.02.18.

ZHOU Y L, WU H M, ZHOU Y D, SHANG X M, PANG L. Effects of different short-term straw returning methods on soil structure and water content. Arid Zone Research, 2022, 39(2): 502-509. doi:10. 13866/j.azr.2022.02.18. (in Chinese)

[48] MATEO-MARíN N, BOSCH-SERRA à D, MOLINA M G, POCH R M. Impacts of tillage and nutrient management on soil porosity trends in dryland agriculture. European Journal of Soil Science, 2022, 73(1): e13139. doi:10.1111/ejss.13139.

[49] 鞠忻倪, 賈玉華, 甘淼, 金珊, 肖波. 黃土溝壑區不同地形部位土壤大孔隙特征研究. 土壤學報, 2018, 55(5): 1098-1107.

JU X N, JIA Y H, GAN M, JIN S, XIAO B. Characteristics of soil macropores in the gully area of loess plateau as affected by terrain. Acta Pedologica Sinica, 2018, 55(5): 1098-1107. (in Chinese)

[50] 王瑋璐, 賀康寧, 張潭, 王先棒, 張震中. 青海高寒區水源涵養林土壤機械組成和理化性質對其飽和導水率和持水能力的影響. 植物資源與環境學報, 2020, 29(2): 69-77.

WANG W L, HE K N, ZHANG T, WANG X B, ZHANG Z Z. Effects of mechanical components and physical and chemical properties of soil in water conservation forests in cold highland area of Qinghai on its saturated hydraulic conductivity and water holding capacity. Journal of Plant Resources and Environment, 2020, 29(2): 69-77. (in Chinese)

[51] 阮仁杰. 土壤孔隙特征對氧化亞氮排放和硝態氮釋放的影響[D]. 合肥: 安徽農業大學, 2021.

RUAN R J. Nitrate released and nitrous oxide emission under different soil pore characteristics[D]. Hefei: Anhui Agricultural University, 2021. (in Chinese)

Effect of Subsoiling Combined with Straw Returning Measure on Pore Structure of Black Soil

YANG JianJun,1GAI Hao1, ZHANG MengXuan1, CAI YuRong1, WANG LiYan2, WANG LiGang1

1Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences/Key Laboratory of Non-point Source Pollution Control, Ministry of Agriculture and Rural Affairs, Beijing 100081;2Qinggang Meteorological Bureau, Qinggang 151600, Heilongjiang

【Objective】The effect of subsoiling combined with straw returning on soil pore structure of black soil in Northeast China is lack of definite judgment. Aimed at such problem, this research was conducted in order to provide scientific basis for the influence mechanism of this measure on soil structure of black soil and the establishment of reasonable tillage.【Method】In this study, a 5-year field positioning experiment conducted in Qinggang County, Suihua City, Heilongjiang Province, a typical black soil area in Northeast China was used as a platform. Farmers’ conventional treatment (FP), 25 cm subsoiling alone treatment (T2), 25 cm subsoiling combined with straw returning treatment (T3) and 35 cm subsoiling combined with straw returning treatment (T4) were set. The visualization and quantification of soil pore structure were studied using CT scanning technology, and combined with field capacity and bulk density to explore the effect of subsoiling combined with straw returning on the pore structure of black soil. 【Result】The results showed that the two-dimensional and three-dimensional images of soil pores clearly showed that the pore distribution at 0-20 cm soil layer was significantly less than that at 20-40 cm soil layer in all treatments, while the pore distribution of suboiling combined with straw returning treatment (T3 and T4) was significantly higher than that under FP treatment, and their macropores with more complex structure were increased. Quantitative analysis showed that compared with FP treatment, the total porosity of 20-30 cm soil layer under T2 was significantly increased by 103.0% (＜0.05) , which was achieved by significantly increasing the micropores porosity(pore diameter d≤0.50 mm) by 91.3% and mesopores porosity (0.50 mm＜d≤1.00 mm) by 143.5% (＜0.05). While subsowing combined with straw returning treatments (T3 and T4) significantly increased the total porosity of 0-30 cm soil layer by 109.8%-382.7% (＜0.05), which was achieved by significantly increasing the macropores porosity (d＞1.00 mm) by 221.5%-661.7% and the mesopores porosity by 105.4%-544.9% (＜0.05). In addition, compared with FP treatment, subsoiling combined with straw returning (T3 and T4) significantly increased the fractal dimension of soil pores at 0-30cm soil layer by 9.9%-17.7% (＜0.05), decreased the Euler number by 32.4%-66.4% (＜0.05), and significantly increased the field water capacity by 24.2%-40.6% (＜0.05). Further analysis showed that the different pore size porosity and total porosity was significantly positively correlated with field capacity and fractal dimension, but significantly negatively correlated with Euler number (＜0.01). 【Conclusion】Subsoiling combined with straw returning could improve porosity of macropores and mesopores pore of black soil, improve pore structure and connectivity, and increase field capacity, especially the effect of 35 cm subsoiling combined with straw returning treatment is the most significant, which could be recommended measure for rational tillage construction in black soil of Northeast China.

Northeast black soil; subsoiling combined with straw returning; CT scanning technology; pore structure in soil

10.3864/j.issn.0578-1752.2023.05.007

2022-07-11；

2022-09-14

國家自然科學基金（32171564）、中國農業科學院科技創新工程項目（CAASZDRW202202）、現代農業產業技術體系北京市創新團隊（BAIC08-2022）

楊建君，E-mail：officeyjj@tom.com。通信作者王立剛，E-mail：wangligang@caas.cn

（責任編輯 李云霞）