張澤宇,馬波,2,李占斌,2*,曾建輝,孫寶洋,趙玉泉

凍融交替作用對黃土高原三種典型土壤可蝕性影響

張澤宇1,馬波1,2,李占斌1,2*,曾建輝1,孫寶洋3,趙玉泉4

1. 西北農(nóng)林科技大學水土保持研究所,黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室, 陜西 楊凌 712100 2. 中國科學院水利部水土保持研究所, 陜西 楊凌 712100 3. 長江科學院水土保持研究所, 湖北 武漢 430010 4. 西北農(nóng)林科技大學資源環(huán)境學院, 陜西 楊凌 712100

為研究凍融交替作用對黃土高原北部土壤可蝕性的影響，選取黃土高原北部季節(jié)性凍融區(qū)的三種典型土壤（灌淤土、黃綿土和風沙土）作為研究對象，設計室內(nèi)凍融模擬試驗，利用基于機械組成和有機質(zhì)含量的EPIC模型與Dg模型估算土壤可蝕性值，系統(tǒng)分析凍融交替作用對土壤機械組成、有機質(zhì)含量及土壤可蝕性值的影響。結(jié)果表明：（1）凍融交替作用使三種供試土壤粘粒含量顯著減少（<0.05），土壤質(zhì)地粗化；（2）土壤有機質(zhì)含量隨凍融交替次數(shù)的增加總體呈先增加后減小再增加的趨勢；（3）相同凍融條件下，利用EPIC模型與Dg模型估算出的土壤可蝕性值具有顯著差異（<0.05），利用EPIC模型估算時，土壤可蝕性值大小表現(xiàn)為風沙土>灌淤土>黃綿土，利用Dg模型的估算結(jié)果與之相反。

凍融交替; 高原土壤; 可蝕性

黃土高原北部季節(jié)性凍融區(qū)的土壤在秋末及初春季節(jié)會發(fā)生“夜凍晝?nèi)凇钡默F(xiàn)象，該現(xiàn)象稱為凍融交替作用，其通過影響土壤的理化性質(zhì)，使土壤更易成為水力侵蝕、風力侵蝕等的侵蝕物質(zhì)來源，加劇土壤侵蝕，嚴重影響該區(qū)域土壤侵蝕狀況及農(nóng)業(yè)生產(chǎn)活動[1-3]。對此，國內(nèi)外學者做了大量研究，研究發(fā)現(xiàn)凍融交替作用能夠減小土壤容重[4,5]，增大土體滲透性、使松散土的孔隙減小、密實土的孔隙增大[6]，改變土壤粘聚力和抗剪強度[7]，影響土壤團聚體含量及其穩(wěn)定性[8]，改變土壤顆粒和機械組成[9]，有機質(zhì)含量的變化受凍融程度、土壤質(zhì)地等因素的綜合影響[10]。

土壤可蝕性（Soil Erodibility）是指土壤對侵蝕的敏感性，是評價土壤對侵蝕敏感程度的重要定量化參數(shù)，國際上一般采用K值即土壤可蝕性因子來表示[11]。經(jīng)過幾十年國內(nèi)外學者的大量研究，目前關(guān)于土壤可蝕性的測定方法已經(jīng)比較成熟，土壤理化性質(zhì)的測定、放水沖刷試驗、模擬降雨和小區(qū)試驗被用于土壤水蝕可蝕性的研究中[12]，并提出諸多模型，其中使用較廣泛的模型有USLE、RUSLE、EPIC模型、Dg模型、Torri.D模型等。張科利等[13]將USLE、EPIC模型和Dg模型應用于我國各地區(qū)土壤侵蝕，并提出適用于我國的三種模型修正公式。

本研究旨在通過室內(nèi)模擬凍融試驗分析凍融交替作用對土壤機械組成和有機質(zhì)含量的影響，利用模型及修正式對土壤可蝕性值進行估算，進而研究凍融對土壤可蝕性的影響，對比利用兩種模型估算得到的結(jié)果，分析其優(yōu)缺點，為該地區(qū)土壤可蝕性研究提供一定的科學依據(jù)。

1 材料與方法

1.1 供試土壤

本研究供試土壤為黃土高原北部季節(jié)性凍融區(qū)的三種典型土壤：陜西省安塞黃綿土、內(nèi)蒙古包頭市灌淤土和內(nèi)蒙古鄂爾多斯達拉特旗風沙土，于撂荒地采取0～30 cm土壤。將供試土壤自然風干后去除植物根系秸稈及石塊等雜質(zhì)，過5 mm篩備用。如表1所示為供試土壤基本性質(zhì)。

表1 三種供試土壤基本性質(zhì)

1.2 室內(nèi)模擬凍融試驗

將風干過篩后的土樣放入絕熱泡沫盒（內(nèi)徑長72.5×寬40×深34.9 cm3）中，根據(jù)調(diào)查和實測容重與預設土壤前期含水量（4%、6%、8%、15%，風沙土增設2%）分層裝土并壓實。填土完成后，將土壤表面刮平。用保鮮膜覆蓋并包裹保溫盒，密封放置24 h，以使水分在土體中均勻分布滲透。將配好土樣放入調(diào)好溫度的冰柜中，打開盒蓋，用保鮮膜覆蓋表層，防止冰柜上蓋板冷凝水進入土體，在冰柜中凍結(jié)12 h。凍結(jié)后取出將保鮮膜揭開，在室溫下融化12 h。本研究中，土體凍結(jié)12 h然后融化12 h的過程稱為一次凍融交替。凍融交替次數(shù)設置為：0、1、3、5、7、10、15、20次。土樣在室溫條件下融化12 h后，采取土體表層（0 ～ 5 cm）土壤于自封袋中，將采集好的土壤自然風干，利用馬爾文激光粒度儀測定其機械組成，使用重鉻酸鉀外加熱法測定有機質(zhì)含量，每組設置3次重復實驗。

1.3 土壤可蝕性K值估算與修正

本研究中，利用基于土壤機械組成和有機質(zhì)含量的EPIC模型和Dg模型估算土壤可蝕性值。1990年Williams等[14]提出的侵蝕-生產(chǎn)力模型EPIC中的值作為衡量土壤可蝕性的指標，見式（1）：

式中：SAN為砂粒含量(%)；SIL為粉粒含量(%)；CLA為粘粒含量(%)；C為有機碳含量(%)；SN1=1-SAN/100。

1984年Shirazi等[15]提出在土壤理化性質(zhì)資料有限的情況下，建議只考慮土壤幾何平均粒徑（Dg）來計算土壤可蝕性值，計算公式（2）為：

Dg=e(0.01?lnm)(3)

式中，f為原土壤中第個粒徑級的等級質(zhì)量分數(shù)，%；m為第個粒徑級兩端數(shù)值的算術(shù)平均值，mm；公式中值的單位是美國制，計算后將值乘0.1317，轉(zhuǎn)化為國際制單位[t·hm2·h /(MJ·m·hm2）]。

本研究利用張科利等在研究我國土壤可蝕性值估算提出的修正式對土壤可蝕性值進行修正轉(zhuǎn)換[13]，修正式如下：EPIC=-0.01383+0.51575williams(=0.613，=0.106) (4)

Dg=-0.00911+0.55066shirazi(=0.705，=0.051) (5)

1.4 相對變化量

本研究中三種土壤質(zhì)地不同，為方便統(tǒng)一比較凍融交替作用對土壤理化性質(zhì)的影響，對數(shù)據(jù)進行歸一化，用相對變化量來表征變化程度，計算公式：

其中，未凍融與凍融后指土壤凍融交替前后測得的指標數(shù)據(jù)。

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

本研究試驗數(shù)據(jù)利用軟件Excel 2019和SPSS 26進行數(shù)理統(tǒng)計分析及圖表制作，進行方差分析和配對樣本檢驗，用最小顯著性差異法（LSD）對試驗結(jié)果進行多重比較。

2 結(jié)果分析

2.1 凍融交替對土壤機械組成的影響

本研究的三種供試土壤中，灌淤土粘粒（<0.002 mm）含量最高，為28.55%，其次是黃綿土粘粒含量為17.59%，風沙土粘粒含量最少，為6.73%；黃綿土粉粒（0.002～0.05 mm）含量最高，為68.63%，其次是灌淤土粉粒含量為44.78%，風沙土粉粒含量最少，為31.96%；風沙土砂粒（0.05～2 mm）含量最高，為61.32%，其次是灌淤土砂粒含量為26.67%，黃綿土砂粒含量最少，為13.78%。在凍融交替過程中三種供試土壤粘粒、粉粒和砂粒含量的變化呈現(xiàn)明顯的規(guī)律：灌淤土，黃綿土和風沙土粘粒（< 0.002 mm）含量隨凍融交替次數(shù)的增加均顯著減小（<0.05），如圖1所示，風沙土粘粒含量降低幅度最大，最大降幅達48.37%，灌淤土和黃綿土粘粒含量的最大降幅分別為19.08%和28.71%。初始含水量越大，凍融交替后粘粒含量減小得越多。黃綿土粉粒（0.002～0.05 mm）含量總體呈減小趨勢，灌淤土和風沙土粉粒含量總體呈增加趨勢；土壤砂粒（0.05～2 mm）含量總體呈增加趨勢，均無顯著性差異（> 0.05）。在凍融交替15次后，三種土壤粘粒、粉粒和砂粒含量的總體趨于穩(wěn)定。

圖1 三種土壤粘粒含量相對變化量隨凍融交替次數(shù)的變化

Fig.1 The relative change of three kinds of soil clay content with the number of freeze-thaw cycles

2.2 凍融交替對土壤有機質(zhì)含量的影響

本研究中，三種土壤在未進行凍融處理時的有機質(zhì)含量表現(xiàn)為：灌淤土（12.33 g/kg）>黃綿土（5.38 g/kg）>風沙土（3.31 g/kg）。在凍融過程中，灌淤土有機質(zhì)含量的變化范圍為12.19～12.76 g/kg，黃綿土有機質(zhì)含量的變化范圍為5.25～5.83 g/kg，風沙土有機質(zhì)含量的變化范圍為3.25～3.52 g/kg。隨著凍融交替次數(shù)的增加，灌淤土、黃綿土和風沙土的有機質(zhì)含量變化差別較大，但總體呈現(xiàn)先增大后減小再增大的變化規(guī)律。黃綿土和風沙土的有機質(zhì)含量變化幅度逐漸增大，相反，灌淤土的有機質(zhì)含量變化幅度逐漸減小。如圖2所示，灌淤土、黃綿土和風沙土有機質(zhì)含量分別在第3次和第20次凍融交替后達到最高水平，較未凍融土壤分別增加2.64%、8.28%和5.26%。在凍融交替10次和15次時，不同初始含水量條件下三種土壤的有機質(zhì)含量均達到最低水平。除初始含水量為4%的風沙土外，不同初始含水量條件下三種土壤經(jīng)過20次凍融交替后，土壤有機質(zhì)含量均增加。

圖2 三種土壤有機質(zhì)含量相對變化量隨凍融交替次數(shù)的變化

Fig.2 The relative change of three kinds of soil organic matter content with the number of freeze-thaw cycles

2.3 凍融交替對土壤可蝕性K值的影響

利用基于土壤機械組成和有機質(zhì)含量的EPIC模型和Dg模型計算三種土壤在凍融交替條件下的土壤可蝕性值，并用公式（4）和（5）進行修正。配對樣本檢驗結(jié)果表明，相同條件下，兩種方法估算出的值具有顯著差異（<0.05）。利用EPIC模型計算的凍融前后3種土壤值（EPIC）變異系數(shù)均大于Dg模型（Dg）。利用 EPIC 模型估算的土壤可蝕性值表現(xiàn)為風沙土>灌淤土>黃綿土，與利用Dg模型的估算結(jié)果相反。在凍融過程中，灌淤土、黃綿土和風沙土EPIC在凍融過程中變化范圍為0.0068 ～ 0.0091、0.0054 ～ 0.0065和0.0240 ～ 0.0314 t·hm2·h /( MJ·mm·hm2），KDg變化范圍為0.0167 ～ 0.0185、0.0179 ～ 0.0187和0.0110 ～ 0.0135 t·hm2·h /( MJ·mm·hm2）。

圖3 三種土壤KEPIC和KDg相對變化量隨凍融交替次數(shù)的變化

Fig.3 The relative changes of KEPIC and KDg of three soils with the number of freeze-thaw cycles

如圖3所示，不同含水量條件下，隨凍融交替次數(shù)的增加，灌淤土和黃綿土EPIC和Dg均呈增加趨勢,無顯著差異(>0.05)，在15次凍融交替后，EPIC和Dg趨于穩(wěn)定。風沙土EPIC隨凍融交替次數(shù)的增加總體呈增加趨勢，相反Dg隨凍融交替次數(shù)的增加呈減小趨勢，均無顯著差異(>0.05)。灌淤土在凍融交替15次后變化量達到最大，EPIC相對變化量為24.95%，Dg相對變化量為7.69%；黃綿土在凍融交替20次后變化量達到最大，EPIC相對變化量為16.06%，Dg相對變化量為5.29%；風沙土在凍融交替10次后變化量達到最大，EPIC相對變化量為12.51%，Dg相對變化量為18.33%。分析可知，EPIC模型和Dg模型應用于灌淤土和黃綿土時得到的變化趨勢總體相近，應用于風沙土時，兩種模型得到的結(jié)果差別巨大。

3 討 論

凍融交替過程中，發(fā)生“凍后聚熵”的現(xiàn)象，溫度變化間接導致土體中水勢變化，從而驅(qū)動水分運移。在凍結(jié)過程中，土壤溫度由上至下依次降低，上層土壤溫度下降時，水汽壓降低，低于下層土壤水汽壓，形成汽壓梯度，下層土壤中的水汽向上層移動；在融化過程中，土壤溫度由上至下依次升高，上層土壤中的水分向下移動。同時，凍融作用會使土壤的孔隙度增大，這為細小顆粒的移動提供了通道。粘粒融化后在水分運動和重力作用下，逐漸向下層移動，使得表層土樣的粘粒含量均顯著減少，且減小幅度受凍融交替次數(shù)和初始含水量影響，凍融交替次數(shù)越多，初始含水量越大，凍融程度越大，粘粒含量減少越多。已有研究發(fā)現(xiàn)，凍融作用會使土壤的干容重減小，且略小于下層土壤的干容重[16]。究其原因可能為水分在向下運移過程中，會挾帶部分細顆粒。表層土壤粉粒和砂粒含量增大，是由于土壤顆粒含量表示為質(zhì)量分數(shù)，當粘粒含量減少，粉粒和砂粒所占比重增加，故含量增大。

使用模型估算土壤可蝕性值時，不同公式的理化性質(zhì)指標不同，其敏感度也不同，Dg模型以土壤機械組成為基礎(chǔ)估算K值，而EPIC模型還需土壤有機質(zhì)含量，凍融交替對有機質(zhì)含量有影響，因此，兩種模型估算的土壤可蝕性K值結(jié)果存在差異。Dg模型公式是土壤可蝕性值關(guān)于幾何平均粒徑Dg的函數(shù)關(guān)系式，研究發(fā)現(xiàn)，當Dg=0.0211 mm時，值最大；當Dg<0.0211 mm時，值與Dg呈正相關(guān)關(guān)系；當Dg>0.0211 mm時，值與Dg呈負相關(guān)關(guān)系。在本研究中，黃綿土和灌淤土的幾何平均粒徑Dg均略小于0.0211 mm，風沙土的幾何平均粒徑Dg遠大于0.0211 mm，故出現(xiàn)風沙土的土壤可蝕性值小于黃綿土和灌淤土的結(jié)果。這與王彬等[17]和Romkens等[18]研究相同，土壤過于粘重（Dg減?。┗蛸|(zhì)地過于粗化（Dg增加），均會使土壤可蝕性減小，最大的土壤可蝕性出現(xiàn)在壤質(zhì)土情況。因此，Dg模型的適用條件受到一定限制，EPIC模型更適合黃土高原季節(jié)性凍融區(qū)的侵蝕預報與研究。

4 結(jié) 論

（1）凍融交替過程中，三種土壤的粘粒含量顯著減小（<0.05），土壤質(zhì)地粗化，灌淤土、黃綿土和風沙土凍融前后相對變化量達到-19.08%、-28.71%和-48.73%；

（2）凍融交替過程中，三種土壤的有機質(zhì)含量總體呈先增加后減小再增加的趨勢；

（3）相同凍融交替條件下，EPIC模型與Dg模型估算出的土壤可蝕性值具有顯著差異（<0.05），利用EPIC模型估算時，土壤可蝕性值表現(xiàn)為風沙土>灌淤土>黃綿土，相反，利用Dg模型的估算結(jié)果表現(xiàn)為風沙土<灌淤土<黃綿土。

[1] 范昊明,蔡強國.凍融侵蝕研究進展[J].中國水土保持科學,2003,1(4):50-55

[2] 李占斌,李社新,任宗萍,等.凍融作用對坡面侵蝕過程的影響[J].水土保持學報,2015(5):56-60

[3] 張科利,劉宏遠.東北黑土區(qū)凍融侵蝕研究進展與展望[J].中國水土保持科學,2018,16(1):17-24

[4] 肖俊波.季節(jié)性凍融對土壤可蝕性影響的試驗研究[D].楊凌:西北農(nóng)林科技大學,2017

[5] 劉佳,范昊明,周麗麗,等.凍融循環(huán)對黑土容重和孔隙度影響的試驗研究[J].水土保持學報,2009,23(6):186-189

[6] 倪萬魁,師華強.凍融循環(huán)作用對黃土微結(jié)構(gòu)和強度的影響[J].冰川凍土,2014,36(4):922-927

[7] Wang L, Zuo X, Zheng F,. The effects of freeze-thaw cycles at different initial soil water contents on soil erodibility in Chinese Mollisol region [J]. Catena, 2020,193:104615

[8] Edwards LM. The effect of alternate freezing and thawing on aggregate stability and aggregate size distribution of some Prince Edward Island soils [J]. Journal of Soil Science, 1991,42:193-204

[9] Zhang Z, Wei MA, Feng W,. Reconstruction of soil Particle composition during freeze-thaw cycling: a review [J]. Pedosphere, 2016,26(2):167-179

[10] Koponen HT, Martikainen PJ. Soil water content and freezing temperature affect freeze–thaw related N2O production in organic soil [J]. Nutrient Cycling in Agroecosystems, 2004,69(3):213-219

[11] Bouyoucos GJ. The clay ratio as a criterion of susceptibility of soils to erosion [J]. Journal of the American Society of Agronomy, 1935,27:738-741

[12] 孫寶洋,李占斌,肖俊波,等.凍融作用對土壤理化性質(zhì)及風水蝕影響研究進展[J].應用生態(tài)學報,2019,30(1):337-347

[13] 張科利,彭文英,楊紅麗.中國土壤可蝕性值及其估算[J].土壤學報,2007,44(1):7-13

[14] Williams JR. EPIC-erosion/productivity impact calculator: 1. Model documentation [J]. Technical Bulletin - United States Department of Agriculture, 1990,4(4):206-207

[15] Shirazi MA, Hart JW, Boersma L. A unifying quantitative analysis of soil texture: improvemnt of precision and extension of scale [J]. Soil Science Society of America Journal, 1988,52(1):181-190

[16] 溫美麗,劉寶元,魏欣,等.凍融作用對東北黑土容重的影響[J].土壤通報,2009,40(3):492-495

[17] 王彬,鄭粉莉, R?mkans MJM.水蝕過程的土壤可蝕性研究述評[J].水土保持研究,2013,20(1):277-286

[18] R?mkens MJM, Poesen J, Wang JY. Relationship between the USLE Soil Erodibility Factor and soil properties [C]//International soil conservation conference on Land Conservation For Future Generations, II. System Analysis for erosion evaluation. Bangkok, Thailand, 1988:371-385

Influence of Freeze-Thaw Alternation on the Erosibility of Three Typical Soils in the Loess Plateau

ZHANG Ze-yu1, MA Bo1,2, LI Zhan-bin1,2*, ZENG Jian-hui1, SUN Bao-yang3, ZHAO Yu-quan4

1.712100,2.712100,3.430010,4.712100,

In order to study the influence of the alternation of freezing and thawing on the soil erodibility in the northern Loess Plateau, three typical soils (irrigated silt soil, loess soil and aeolian sand soil) in the seasonal freezing and thawing area of the northern Loess Plateau were selected as the research objects, and an indoor freeze-thaw simulation was designed in the experiment, the EPIC model and the Dg model based on mechanical composition and organic matter content are used to estimate the soil erodibilityvalue, and the influence of freeze-thaw alternation on soil mechanical composition, organic matter content and soil erodibilityvalue is systematically analyzed. The results showed that: (1) The freeze-thaw alternation effect significantly reduced the clay content of the three types of soils (<0.05), and the soil texture became coarser; (2) The soil organic matter content increased first and then decreased with the increase in the number of freeze-thaw alternations. (3) Under the same freezing and thawing conditions, the soil erodibilityvalue estimated by the EPIC model and the Dg model are significantly different (<0.05). When the EPIC model is used to estimate the soil erodibility, thevalue of soil erodibility is in the order of aeolian sandy soil> irrigated silt soil> loessal soil, and the estimation result using the Dg model is contrary to this.

Freeze-thaw alternation; plateau soil; erodibility

S155.2+94

A

1000-2324(2021)05-0813-06

2020-12-25

2021-05-12

國家自然科學基金面上項目:黃土高原凍融對水蝕過程的作用機理研究(41771311);中央級公益性科研院所基本業(yè)務費項目:凍融作用對長江源區(qū)土壤可蝕性影響機制研究(CKSF2019179/TB);中央級公益性科研院所基本業(yè)務費項目:江源地區(qū)水循環(huán)及生態(tài)境演變與適應性保護研究(CKSF2019292/SH);西北農(nóng)林科技大學大學生創(chuàng)新訓練項目:凍融交替對土壤抗蝕性的影響研究(S202010712569)

張澤宇(1996-),男,碩士研究生,主要從事土壤侵蝕研究. E-mail:m18392426454@163.com

通訊作者:Author for correspondence.E-mail:zhanbinli@126.com