范廷玉 鐘建 王順

摘 要：為探討高潛水位采煤拉張裂隙區(qū)土壤中元素遷移特征，選擇淮北礦區(qū)朱莊礦3522工作面地表拉張裂隙區(qū)為研究對象，按地勢坡度從高到低劃分為Ⅰ、Ⅱ、Ⅲ區(qū)，在裂隙兩側(cè)間隔50cm布設(shè)21個觀測點(diǎn)，每個觀測點(diǎn)自上而下每20cm取樣，共計(jì)105個樣品，分析含水率、pH、有機(jī)質(zhì)（SOM）、有效磷（AP）和速效鉀（AK）等指標(biāo)。結(jié)果表明：拉張裂隙區(qū)土壤剖面含水率隨深度增大;SOM在Ⅱ區(qū)臺階0～40cm地表聚集。AP在垂向剖面上分布無顯著差異，在水平方向上在Ⅱ區(qū)、Ⅲ區(qū)顯著積累;裂隙對0～20cm AK影響顯著，土壤整體肥力屬于瘦瘠到一般水平，斜坡水土流失與裂隙的隔斷攔截作用，導(dǎo)致SOM在臺階區(qū)積累，AP在低區(qū)積累，AK在高區(qū)積累。

關(guān)鍵詞：高潛水位;地表拉張裂隙;土壤特征

Abstract：In order to explore the characteristics of element migration in the soil of tension fracture area in high phreatic water level coal mining， the surface tension fracture area of 3522 working face in Zhuzhuang Mine of Huaibei mining area was selected as the research object. According to the terrain slope， it was divided into zone I， zone II and zone III. 21 observation points were set at 50 cm intervals on both sides of the fracture， and 105 samples were taken from each observation point every 20 cm from top to bottom. Organic matter （SOM）， available phosphorus （AP）， available potassium （AK）and other indicators were monitored. The results were showed as following：the water content of soil profile increased with the depth in the tension fracture area， and SOM accumulated on the 0～40cm surface of the step in area II. There was no significant difference in the vertical distribution of AP， but AP accumulated significantly in area II and Area III in the horizontal direction. The surface tension cracks had a significant effect to the AK concentration of 0～20cm soil depth， while the soil fertility belonged to the general level. The partition and interception of slope soil erosion and fissures lead to the accumulation of SOM in the bench area， AP in the low area and AK in the high area.

Key words：High phreatic water level; surface tension fractures; soil characteristics

煤炭是我國的主要能源，2018年我國能源生產(chǎn)結(jié)構(gòu)中，原煤占比68.3%，能源消耗結(jié)構(gòu)中煤炭占57.7%。井工開采是主要的采煤方式，導(dǎo)致地表塌陷并產(chǎn)生地表拉張裂隙[1-2]。根據(jù)采煤工作面生產(chǎn)與對應(yīng)地表拉張裂隙發(fā)育情況，可將地表拉張裂隙分為橫向裂隙（與切眼平行）與縱向裂隙（與風(fēng)巷、機(jī)巷平行），橫向裂隙位于沉陷盆地底部，沿采煤推進(jìn)方向不斷拉張與閉合，規(guī)模小且淺;縱向裂隙則分布在沉陷盆地區(qū)，形成后不斷拉張，規(guī)模相對大且深。地表沉陷裂隙導(dǎo)致了沉陷區(qū)土壤呈現(xiàn)高度空間異質(zhì)性[3-4]，水土流失[5]、優(yōu)先流[6-8]、微地形改造[9]中粘粒、粉粒等易隨著地表徑流遷移，眾多學(xué)者圍繞土壤質(zhì)量改變與植被群落結(jié)構(gòu)退化[10]，圍繞微生物群落多樣性[11-12]、酶活性[13]等恢復(fù)土壤活性進(jìn)行了研究，但目前針對沉陷區(qū)土壤性質(zhì)研究并未區(qū)分拉張裂隙類型，進(jìn)行土壤水分、營養(yǎng)元素的精細(xì)化空間分布研究。論文以淮北礦業(yè)集團(tuán)朱莊煤礦縱向拉張裂隙為研究區(qū)域，討論該區(qū)域土壤營養(yǎng)元素分布特征，為沉陷區(qū)土壤環(huán)境保護(hù)與土地質(zhì)量提升提供理論依據(jù)，對于提高農(nóng)業(yè)產(chǎn)量具有重要的現(xiàn)實(shí)意義。

1 材料與方法

1.1 區(qū)域概況

淮北礦區(qū)位于安徽省北部，總面積4 650km2[14]，采煤沉陷地面積約272km2[15]，約占總面積5.9%。研究區(qū)位于淮北市杜集區(qū)朱莊礦3522工作面的地表沉陷區(qū)，位于人民中路北側(cè)，東臨岱河，地理坐標(biāo)為116°51′02″～116°51′09″E、33°57′20″～33°57′25″N，3522工作面長150m，寬130m。研究區(qū)屬于暖溫帶半濕潤季風(fēng)氣候，年均降水量823.4mm，年極端最大降水量1 441.4mm。沉陷土壤類型主要為半水成土綱的潮土與砂姜黑土，土地利用類型為林地，潛水位埋深在1～3m。

1.2 點(diǎn)位布設(shè)與樣品采集

選取平行于3522工作面風(fēng)巷的兩條縱向拉張裂隙（C1、C2）作為研究對象，C1、C2中間為沉陷臺階，東西走向長約70m，裂隙寬度2cm～15cm，北高南低，坡度約2°～3°（見圖1）。從高（北）到低（南）布設(shè)一條長25m采樣線，以裂隙為界從平面上將沉陷斜坡分為Ⅰ（高）、Ⅱ（中）、Ⅲ（低）區(qū)。在裂隙兩側(cè)每隔0.5m布設(shè)一個采樣點(diǎn)，共21個點(diǎn)位。每個點(diǎn)位以20cm為間隔使用土壤分層采樣器采集天然土壤樣品，采至1m深，共105個樣品。采樣時間為2019年4月6日，采集的土樣去除石礫和動植物殘?bào)w后混勻裝入寫好標(biāo)簽的密封袋帶回實(shí)驗(yàn)室。

1.3 樣品前處理與測試

將帶回新鮮土樣，部分測定土壤含水率。剩余用無菌手套捏碎混勻，去除雜草和石礫等。均勻的鋪在干凈的紙上，于室內(nèi)自然風(fēng)干，采用四分法用研缽研磨，分別過10目、20目和60目篩，裝入聚乙烯塑料袋保存。按照文獻(xiàn)[16]分別測定土壤pH、有機(jī)質(zhì)、有效磷、速效鉀等指標(biāo)。其中pH測定取10目土樣8g于離心管中，采用液土比2.5混合振蕩靜置30min后測定上清液pH值，采用重鉻酸鉀加熱法測定有機(jī)質(zhì)，采用鉬銻抗比色法測定有效磷，采用四苯硼鈉比色法測定速效鉀。測試過程中，每批樣品設(shè)置兩個空白試驗(yàn)，兩個測定結(jié)果的相對偏差小于等于50%。每20個樣品測定10%的平行雙樣，兩個測定結(jié)果的相對偏差小于等于10%，校準(zhǔn)曲線的相關(guān)系數(shù)均大于0.999。

1.4 評價方法

土壤肥力系數(shù)采用修正的內(nèi)梅羅綜合指數(shù)法[17]評價，土壤各屬性分級指標(biāo)（見表1）參照第二次全國土壤普查分級標(biāo)準(zhǔn)和《土壤分析技術(shù)規(guī)范》[18]，綜合考慮研究區(qū)土壤各指標(biāo)值范圍確定。

采用標(biāo)準(zhǔn)化計(jì)算公式消除各參數(shù)之間的量綱差異，將計(jì)算得出的肥力系數(shù)與肥力分級標(biāo)準(zhǔn)對比可知肥力情況，肥力分級標(biāo)準(zhǔn)：P≥2.7，很肥沃;1.80≤P<2.70，肥沃;0.90≤P<1.80，一般;P<0.90，瘦瘠。

2 結(jié)果與討論

2.1 土壤理化性質(zhì)分布特征

如圖2（a）～（b）所示，土壤含水率隨著土層深度增加，沿著地勢降低，距離裂隙的50～100cm范圍內(nèi)的0～20cm土層含水率最低。Ⅰ區(qū)、 Ⅱ區(qū)、 Ⅲ區(qū)的含水率逐級下降（見圖2（c）～（d）），均值分別為13.73%、11.4%、8.53%，最小值均在表層土壤，分別為5%、5%、4%，Ⅰ-Ⅱ區(qū)、Ⅱ-Ⅲ區(qū)含水率分別下降18%、25%。pH值8.21～9.04，均值8.66，不同深度的pH值變化規(guī)律相似，且裂隙兩側(cè)50cm范圍的pH值最大，伴隨到裂隙距離增加緩慢降低。土壤粘粒含量（見表2）隨著深度增加而增大，Ⅰ區(qū)、Ⅱ區(qū)和Ⅲ區(qū)的粘粒含量均值分別為5.42%、3.41%和3.03%。

淺層土（0～40cm）中有機(jī)質(zhì)在Ⅱ區(qū)累積， 有效磷含量在水平分布上（見圖3）呈現(xiàn)顯著差異， Ⅰ區(qū)含量顯著低于Ⅱ區(qū)和Ⅲ區(qū)。 在0～40cm， Ⅱ區(qū)大于Ⅲ區(qū);在40～100cm，Ⅱ區(qū)小于Ⅲ區(qū)。有效磷含量在距裂隙50cm處（見圖4）Ⅲ區(qū)表層土中的最高，達(dá)5.35mg·kg-1，隨距裂隙長度增加，Ⅰ區(qū)、Ⅲ區(qū)有效磷均呈現(xiàn)先下降再增加的趨勢。速效鉀含量11.82～222.51mg·kg-1，均值86.41mg·kg-1。在水平分布上呈現(xiàn)與有效磷相反的趨勢，在0～40cm土層，Ⅱ區(qū)含量最高，Ⅲ區(qū)含量最低;40～100cm土層，最高的為Ⅰ區(qū)，Ⅱ區(qū)和Ⅲ區(qū)差異并不顯著。

降雨期間的優(yōu)先流與非降雨期間的裂隙增加土壤蒸發(fā)，均導(dǎo)致表層土壤含水率不同程度下降[14-20]，其中以0～20cm表層土蒸散發(fā)量最高，距離裂隙較近位置的土壤含水率會明顯小于較遠(yuǎn)處的土壤，影響范圍在50～100cm左右，有研究表明在室內(nèi)模擬實(shí)驗(yàn)中得出裂隙影響范圍為距裂隙距離45cm[21]，Ⅲ區(qū)更加接近沉陷盆底，裂隙更加密集，土壤持水能力下降幅度最大。

土壤鹽分和可溶性養(yǎng)分與水分運(yùn)移規(guī)律一致，拉張裂隙區(qū)土壤含水率的下降導(dǎo)致其他水溶性指標(biāo)分布的差異性，如pH值上升、有效磷隨著地勢降低升高等。沉陷裂隙的阻截作用[22]、縱向拉張裂隙區(qū)中水土流失與裂隙優(yōu)先流[23-24]并存，造成土壤有機(jī)質(zhì)、有效磷和速效鉀空間分布差異性。土壤中有機(jī)物的主要載體為粘粒[25]，在降雨等外力作用下會隨地表徑流遷移[26]，向裂隙處聚集、隨裂隙優(yōu)先流向下遷移[27-28]。土壤表層的磷和鉀等元素具有向底土層轉(zhuǎn)移和聚集的趨勢[29]。

淮北礦區(qū)屬于典型的高潛水礦區(qū)，土壤含水量與潛水位之間關(guān)系緊密。有研究表明淮北平原土壤含水率變化除受氣象、人類活動因素影響外，潛水埋深影響也十分顯著[30]。豐水期潛水位大概在1m左右，對深層土壤（60～100cm）中的土壤水分產(chǎn)生較大影響，土壤粘粒、養(yǎng)分等可以伴隨高潛水位向下遷移，潛水位隨著季節(jié)波動，也會導(dǎo)致深層土壤各指標(biāo)變化的復(fù)雜性。

2.2 土壤肥力綜合評價

由表3可知，Ⅰ區(qū)土壤 pH 系數(shù)為2.74，屬很肥沃水平，Ⅱ和Ⅲ土壤pH< 2.70，屬于肥沃水平。Ⅰ、Ⅱ、Ⅲ區(qū)速效鉀肥力系數(shù)分別為1.47～1.75，屬于一般到肥沃水平。有機(jī)質(zhì)和有效磷的肥力系數(shù)分別為0.19～0.25和0.28～0.71，屬于瘦瘠水平，有機(jī)質(zhì)、磷養(yǎng)分是肥力較低的限制因子。整體而言，Ⅰ、Ⅱ、Ⅲ裂隙土壤肥力綜合系數(shù)介于0.89～0.94，屬于瘦瘠到一般水平。Ⅰ、Ⅱ、Ⅲ區(qū)土壤綜合肥力系數(shù)隨土壤深度的變化如圖5所示，Ⅰ、Ⅱ、Ⅲ區(qū)土壤綜合肥力系數(shù)最大值均在0～20cm的表層土，在0～40cm的土層土壤肥力系數(shù)Ⅱ區(qū)>Ⅰ區(qū)>Ⅲ區(qū)，在40～100cm土層土壤肥力系數(shù)為Ⅲ區(qū)>I區(qū)> Ⅱ區(qū)，0～40cm和40～100cm土壤肥力呈相反趨勢。

將肥力系數(shù)和各指標(biāo)做相關(guān)性分析，如表4所示，土壤肥力系數(shù)與土壤的pH無顯著相關(guān)，但是與有機(jī)質(zhì)、有效磷及速效鉀均極顯著相關(guān)，其相關(guān)系數(shù)分別為0.661、0.501、0.84。

土壤養(yǎng)分流失主要分為伴隨徑流流失和伴隨泥沙流失兩種形式，其中有機(jī)質(zhì)和N主要通過徑流流失，而P和K則主要伴隨泥沙流失[31]，該結(jié)論與本研究中P隨深度增加流失量不顯著的現(xiàn)象正好相反，這可能是因?yàn)樵撗芯繀^(qū)域土壤P嚴(yán)重缺乏而K含量較高所致[32]。其中K主要富集在0.01～0.05mm粒徑中，隨徑流遷移以粘粒居多[33]。I和Ⅱ區(qū)土壤綜合肥力系數(shù)大致隨裂隙距離的增加而增加，這可能是由于裂隙的產(chǎn)生，增加了土壤與空氣的接觸面積，使得水土流失現(xiàn)象加劇，致使肥力損失;而Ⅲ區(qū)隨著距采動地裂隙距離的增加，土壤肥力綜合系數(shù)大致呈下降趨勢，這可能與當(dāng)?shù)氐牡匦蔚貏萦嘘P(guān)，I區(qū)地勢較高，Ⅱ區(qū)位于中間較平坦，Ⅲ區(qū)地勢較低，土壤肥力綜合系數(shù)大致隨地勢的降低而減小。

2.3 裂隙區(qū)土壤營養(yǎng)元素遷移特征

井工采煤作業(yè)過程導(dǎo)致地表沉陷，水土流失是營養(yǎng)元素遷移的首要驅(qū)動力，降雨過程中大雨強(qiáng)下的徑流泥沙攜帶、小雨強(qiáng)下的徑流水?dāng)y帶以及超滲產(chǎn)流[34]，可礦化泥沙結(jié)合態(tài)的養(yǎng)分、可溶性溶解態(tài)養(yǎng)分[35-36]均會伴隨徑流遷移，沉陷邊界頂部與底部土壤有機(jī)質(zhì)含量比中部明顯下降，有機(jī)質(zhì)在沉陷臺階處聚集;磷具有很強(qiáng)的吸附性[37]，淋失量小，短歷暴雨事件產(chǎn)生的侵蝕泥沙是磷流失的主要載體[38-39] ，約60%磷流失是通過泥沙攜帶遷移，在沉陷斜坡上，有效磷在沉陷盆底聚集;速效鉀易溶于水，當(dāng)斜坡降雨強(qiáng)度小于土壤的入滲能力時，部分速效鉀會溶于水并全部入滲深層土壤;當(dāng)降雨強(qiáng)度大于入滲能力時，會產(chǎn)生徑流，速效鉀隨著徑流流失，研究區(qū)的高潛水位對沉陷盆底土壤影響更大，速效鉀隨著潛水向下遷移，導(dǎo)致速效鉀的含量在沉陷斜坡頂部最高。結(jié)合土壤性質(zhì)分析和肥力評價，伴隨降雨的水土流失、裂隙攔截、高潛水位三者共同作用導(dǎo)致了縱向裂隙區(qū)的養(yǎng)分空間分布差異性，營養(yǎng)元素的遷移特征如圖6所示。

3 結(jié)論

（1）在高潛水位采煤拉張裂隙區(qū)土壤含水率、pH、有機(jī)質(zhì)、有效磷和速效鉀進(jìn)行了研究，發(fā)現(xiàn)在沉陷坡度和裂隙的共同作用下，垂直方向上含水率會隨著地勢的降低而降低，隨著土層厚度增加而增大?？傮w上Ⅲ區(qū)含水率較Ⅱ區(qū)下降25%，Ⅱ區(qū)較Ⅰ區(qū)下降18%，水平方向上在裂隙附近的土壤含水率明顯低于距離裂隙較遠(yuǎn)處的含水率，裂隙影響范圍在100cm左右。

（2）在縱向拉張裂隙區(qū)，土壤的pH無顯著差異，有機(jī)質(zhì)在地勢較平坦的Ⅱ區(qū)0～40cm地表聚集，在地勢較高的地方，有機(jī)質(zhì)含量隨著距裂隙距離成正比，在地勢較低處則相反。土壤有效磷在地勢較低處聚集，在地表土距離裂隙不同位置的變化和有機(jī)質(zhì)有相同的規(guī)律。土壤速效鉀在高區(qū)聚集，且在地表0～20cm處含量較高。

（3）通過修正的內(nèi)梅羅指數(shù)法分析，該研究區(qū)土壤整體肥力偏低，屬于瘦瘠到一般水平，有機(jī)質(zhì)是肥力較低的主要限制因子。

（4）高潛水位、水土流失和裂隙攔截作用共同導(dǎo)致了局部營養(yǎng)元素的空間分布差異性，在對沉陷區(qū)土壤提質(zhì)增效時需在不同區(qū)域采取差異性措施，也為沉陷區(qū)土壤環(huán)境保護(hù)提供了理論參考。

參考文獻(xiàn)：

[1] 張兵.開采沉陷動態(tài)預(yù)計(jì)模型構(gòu)建與算法實(shí)現(xiàn)[D].北京：中國礦業(yè)大學(xué)，2017.

[2] 徐平.采動沉陷影響下埋地管道與土相互作用及力學(xué)響應(yīng)研究[D].北京：中國礦業(yè)大學(xué)，2015.

[3] 張淑媛，秦勝金，宗國，等.遼東山區(qū)次生林土壤主要化學(xué)性質(zhì)空間異質(zhì)性[J].生態(tài)學(xué)雜志，2018，37（7）：2 076-2 082.

[4] 許傳陽，馬守臣，張合兵，等.煤礦沉陷區(qū)沉陷裂縫對土壤特性和作物生長的影響[J].中國生態(tài)農(nóng)業(yè)學(xué)報(bào)，2015，23（5）：597-604.

[5] 郭文兵，白二虎，趙高博.高強(qiáng)度開采覆巖地表破壞及防控技術(shù)現(xiàn)狀與進(jìn)展[J].煤炭學(xué)報(bào)，2020，45（2）：509-523.

[6] 張信寶.關(guān)于中國水土流失研究中若干理論問題的新見解[J].水土保持通報(bào)，2019，39（6）：302-306.

[7] 趙令，蘇濤，李杰衛(wèi)，等.采煤沉陷區(qū)地表水土流失時空演化研究[J].煤礦安全，2019，50（4）：37-40.

[8] 常文科，李成杰，郭曉明，等.采煤沉陷區(qū)林地優(yōu)先流現(xiàn)象下的土壤特性分布特征研究[J].河南理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2020，39（6）：51-57.

[9] 李斯佳，王金滿，萬德鵬，等.采煤沉陷地微地形改造及其應(yīng)用研究進(jìn)展[J].生態(tài)學(xué)雜志，2018，37（6）：1 612-1 619.

[10] 畢銀麗，鄒慧，彭超，等.采煤沉陷對沙地土壤水分運(yùn)移的影響[J].煤炭學(xué)報(bào)，2014，39（S2）：490-496.

[11] 劉英，雷少剛，程林森，等.采煤塌陷影響下土壤含水率變化對檸條氣孔導(dǎo)度、蒸騰與光合作用速率的影響[J].生態(tài)學(xué)報(bào)，2018，38（9）：3 069-3 077.

[12] SUN SHIYONG，SUN HUI，ZHANG DESHUN，et al.Response of Soil Microbes to Vegetation Restoration in Coal Mining Subsidence Areas at Huaibei Coal Mine，China[J].International Journal of Environmental Research and Public Health，2019，16（10）：1 757-1 771.

[13] 李媛媛.采煤塌陷地泥漿泵復(fù)墾土壤微生物多樣性及土壤酶活性研究[D].北京：中國礦業(yè)大學(xué)，2015.

[14] 孫浩，周春財(cái)，徐仲雨，等.淮北礦區(qū)土壤重金屬空間分布與環(huán)境評價[J].中國科學(xué)技術(shù)大學(xué)學(xué)報(bào)，2018，48（7）：560-566.

[15] 鄧士強(qiáng).關(guān)于淮北市采煤沉陷區(qū)調(diào)查的幾點(diǎn)認(rèn)識[J].西部資源，2020（2）：65-67.

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

[17] 耿慶龍，李娜，賴寧，等.基于修正內(nèi)梅羅指數(shù)法的果園土壤肥力評價[J].新疆農(nóng)業(yè)科學(xué)，2019，56（11）：2 104-2 111.

[18] 全國土壤普查辦公室.中國土壤普查技術(shù)[M].北京：農(nóng)業(yè)出版社，1992：2-200.

[19] 全占軍，程宏，于云江，等.煤礦井田區(qū)地表沉陷對植被景觀的影響--以山西省晉城市東大煤礦為例[J].植物生態(tài)學(xué)報(bào)，2006，30（3）：414-420.

[20] 黎煒，陳龍乾，周天建，等.我國采煤沉陷地土壤質(zhì)量研究進(jìn)展[J].煤炭科學(xué)術(shù)，2011，39（5）：125-128.

[21] 畢銀麗，伍越，張健，等.采用HYDRUS模擬采煤沉陷地裂縫區(qū)土壤水鹽運(yùn)移規(guī)律[J].煤炭學(xué)報(bào)，2020，45（1）：360-367.

[22] 喬岡，徐友寧，陳華清，等.寧東煤礦區(qū)地裂縫對植被生態(tài)環(huán)境的影響[J].地質(zhì)通報(bào)，2018，37（12）：2 176-2 183.

[23] 郭巧玲，蘇寧，丁斌，等.采煤塌陷區(qū)裂隙優(yōu)先流特征[J].中國水土保持科學(xué)，2018，16（3）：112-120.

[24] LV GANG，LI JUN，LI YEXIN，et al.Preferential flow characteristics of reclaimed mine soils in a surface coal mine dump[J].Environmental Monitoring and Assessment，2017，189（6）：266-287.

[25] 蔡岸冬，張文菊，申小冉，等.長期施肥土壤不同粒徑顆粒的固碳效率[J].植物營養(yǎng)與肥料學(xué)報(bào)，2015，21（6）：1 431-1 438.

[26] 嚴(yán)家平，陳孝楊，程方奎，等.礦區(qū)土壤裂隙優(yōu)先流對土壤銨態(tài)氮遷移及土壤結(jié)構(gòu)的影響[J].農(nóng)業(yè)工程學(xué)報(bào)，2018，34（2）：120-126.

[27] ZHANG HUIJUAN，MA SHOUCHEN，LIU WENKAI，et al.Three-Dimensional Spatial Simulation and Distribution Characteristics of Soil Organic Matter in Coal Mining Subsidence Area[J].Materials Science Forum，2020，6038：437-448.

[28] JING ZHAORUI，WANG JINMAN，ZHU YUCHENG，et al.Effects of land subsidence resulted from coal mining on soil nutrient distributions in a loess area of China[J].Journal of Cleaner Production，2018，177：350-361.

[29] 馬康，楊帆，張玉秀.西北干旱半干旱區(qū)煤炭井工開采對土壤肥力質(zhì)量的影響研究進(jìn)展[J].中國科學(xué)院大學(xué)學(xué)報(bào)，2020，37（4）：442-449.

[30] 賈青，余鐘波，楊傳國，等.淮北平原土壤濕度變化特征及其與潛水埋深響應(yīng)關(guān)系[J].水電能源科學(xué)，2017，35（11）：123-126.

[31] 陳欣，姜曙千，張克中，等.紅壤坡地磷素流失規(guī)律及其影響因素[J].土壤侵蝕與水土保持學(xué)報(bào)，1999（3）：38-41.

[32] 李仁山.廣西南部桉樹人工林土壤性質(zhì)與水土流失特征[D].南寧：廣西大學(xué)，2014.

[33] SHIGAKI F，SHARPLEY A，PROCHNOW L I .Rainfall intensity and phosphorus source effects on phosphorus transport in surface runoff from soil trays[J].The ence of the Total Environment，2007，373（1）：334-343.

[34] 沈紫燕，王輝，平李娜，等.前期土壤含水率對粘性紅壤產(chǎn)流產(chǎn)沙及溶質(zhì)運(yùn)移的影響[J].水土保持學(xué)報(bào)，2014，28（1）：58-62.

[35] 張冠華，劉國彬，易亮.植被格局對坡面流阻力影響的試驗(yàn)研究[J].水土保持學(xué)報(bào)，2014，28（4）：55-59.

[36] DREWRY J J，NEWHAM L T H，CROKE B F W.Suspended sediment，nitrogen and phosphorus concentrations and exports during storm-events to the Tuross estuary，Australia[J].Journal of environmental management，2009，90（2）：879-887.

[37] 王帥兵，宋婭麗，王克勤，等.不同雨型下反坡臺階減少紅壤坡耕地氮磷流失的效果[J].農(nóng)業(yè)工程學(xué)報(bào)，2018，34（13）：160-169.

[38] M L RODRIGUEZ-BLANCO， M M TABOADA-CASTRO， M T Taboada-Castro. Phosphorus transport into a stream draining from a mixed land use catchment in Galicia （NW Spain）：Significance of runoff events[J].Journal of Hydrology，2013，481：12-21.

[39] 喬閃閃，吳磊，彭夢玲.人工模擬降雨條件下黃土坡面水-沙-氮磷流失特征[J].環(huán)境科學(xué)研究，2018，31（10）：1 728-1 735.

（責(zé)任編輯：李 麗）