西北風積沙區采煤擾動下土壤侵蝕與養分演變特征

聶小軍，高 爽，陳永亮，張合兵

?

西北風積沙區采煤擾動下土壤侵蝕與養分演變特征

聶小軍1，高 爽1，陳永亮2，張合兵1

（1. 河南理工大學測繪與國土信息工程學院，焦作 454000； 2. 中國神華神東哈拉溝煤礦，榆林 719315）

為更好地理解西北風積沙礦區生態環境演變規律，以神東哈拉溝與上灣礦區為例，利用137Cs示蹤法分析了未采區、自恢復沉陷區（1、2、4、8 a沉陷區）與植被修復區（13 a沉陷區）的土壤侵蝕與養分特征，研究了礦區土壤侵蝕與養分的演變規律。結果表明：采煤擾動可以導致未采區土壤侵蝕強度增大與有機碳、微生物量碳、全氮、全磷、堿解氮養分明顯損失。開采沉陷后，沉陷區土壤侵蝕強度隨著時間的推移呈現出先增加后降低的趨勢；其中，地表沉陷后的最初2 a是土壤侵蝕急劇增大的時段。采煤擾動下，土壤有機碳、微生物量碳、全氮、全磷、堿解氮養分的演變規律與土壤侵蝕演變密切相關。植被修復可以有效降低沉陷區土壤侵蝕強度與提高土壤養分含量。西北風積沙區采煤沉陷地表的生態恢復應該及時開展地表沉陷后的早起侵蝕防治，兼顧考慮植被修復與土壤微生物聯合修復，以此促進土壤、植被正向演替。

土壤；侵蝕；養分；采煤塌陷；生態環境演變；137Cs示蹤；風積沙區

0 引 言

西北地區是中國重要的煤炭基地，煤炭資源儲量豐富、賦存穩定、基巖薄、煤層厚、埋深淺、地質條件簡單，極適于高強度開采。隨著國家煤炭開采戰略的西移，集群化、高強度開采方式使得地表生態環境的破壞范圍廣、速度快、形式多、時空差異大。加之，西北礦區多處于風積沙區，生態脆弱，突出表現為干旱缺水、植被覆蓋低、抗擾動能力差，在采礦擾動和自然營力的共同作用下，煤炭開采導致的地表生態環境惡化問題極其突出，嚴重制約了礦區的可持續發展。

采煤擾動導致的土壤質量變化是西北風積沙礦區生態環境研究關注的一個重要科學議題。受采煤擾動影響，土壤質量呈動態變化，因此開展該方面的研究應該關注土壤質量的演變。先前大量的研究從靜態的角度評價了地表沉陷后土壤水分、養分等理化性質的變化，探討了地表破壞（裂縫、沉陷）對土壤質量的影響[1-6]。最近，臺曉麗等[7-8]從動態的角度研究了地表沉陷后的土壤水分、養分演變。目前的研究極大地增加了對礦區土壤質量變化的認識，但是這些研究通常以沉陷區附近的未采區為對照，忽略了煤炭開采對未采區生態環境的影響（例如：已采區地下含水層破壞導致整個礦區土壤干旱加劇），這不利于全面理解西北風積沙礦區土壤質量的演變規律。

土壤侵蝕是導致全球土地退化的一個主要原因[9]。西北風積沙區風蝕荒漠化嚴重，采煤擾動導致該環境問題更加突出，但是有關該礦區土壤侵蝕演變的研究鮮有報道。目前用于評價土壤侵蝕的方法主要有修正的通用土壤流失方程（RUSLE，revised universal soil loss equation）與同位素示蹤法。在RUSLE中，SL（坡長、坡度）是評價土壤侵蝕強度的地形因子。由于西北風積沙礦區地貌類型主要為波狀沙丘，地下煤層近水平分布，開采沉陷導致的地形因子SL變化不明顯，這將限制利用RUSLE對該區土壤侵蝕演變的準確評價。與通用土壤流失方程不同，137Cs同位素示蹤法不受地形因子的限制，只需要比較研究區與背景區的137Cs含量差異，即可獲取研究區的侵蝕或沉積狀況；而且，137Cs示蹤法具備簡單、評價精度高、研究費用低等特點，因此該方法目前被證明是研究土壤侵蝕（水蝕、風蝕、耕作侵蝕）的可靠方法[10]。此外，137Cs示蹤法能夠很好地建立土壤侵蝕-養分間的聯系[11-13]，為研究土壤養分侵蝕動態提供了方法論上的支持。

基于以上背景，本研究以神東哈拉溝與上灣礦區為研究區，利用137Cs示蹤法調查未采區、自恢復沉陷區（地表沉陷后穩定期分別為1、2、4、8 a的沉陷區）的土壤侵蝕與養分特征，揭示其演變規律并探討土壤侵蝕對土壤養分演變的影響，同時也檢驗植被修復區（地表沉陷后穩定期為(13+12) a經濟林恢復）的土壤保持效果，以期為西北風積沙礦區的生態恢復與重建提供理論支持。

1 材料與方法

1.1 研究區概況

研究區選擇神東哈拉溝（39°23′54″N，110°12′30″E）與上灣（39°16′52″N，110°09′32″E）礦區（圖1a），地處毛烏素沙漠東南邊緣，海拔1 084～1 350 m；地貌類型主要為波狀固定、半固定沙丘，地表多為第四系風積沙等松散層覆蓋。研究區屬溫帶半干旱半沙漠的高原大陸性季風氣候，降雨稀少且年內變化大，年均降雨量為375 mm（1957－2016），一年的降水主要集中在7～9月中旬，約占全年降水量的60%以上，且多以暴雨的形式出現。研究區蒸發強烈，年均蒸發量為2 382 mm，是年降雨量的6.4倍。春、冬兩季為研究區風季，西北風盛行，風速一般在3.30 m/s左右，最大可達24 m/s（1979年11月）。生態系統主要為荒漠草原，植被以沙柳（）、沙蒿（）、紫穗槐（）和檸條（）等典型荒漠植被為主。土壤類型為風沙土或黃綿土，表層土壤質地為砂土或砂壤，結構較松散，土壤肥力低，風蝕嚴重。研究區煤層埋深較淺，高強度井工開采導致的地表破壞（沉陷、裂縫）、水土流失、植被退化等生態環境問題突出。

圖1 研究區與試驗樣地位置示意圖

1.2 樣地選擇

根據哈拉溝與上灣礦采煤工作面背景資料，2013年在研究區內選擇未采區、自恢復沉陷區（地表沉陷后穩定期分別為1、2、4、8 a的沉陷區，下文分別簡稱為1、2、4、8 a沉陷區）、植被修復沉陷區（地表沉陷后穩定期為(13+12) a經濟林恢復）來開展研究（圖1b）。未采區、4個自恢復沉陷區以及植被修復沉陷區的供試樣地面積均為0.15 km2。未采區、自恢復沉陷區與植被修復沉陷區地貌、土壤類型一致，地貌類型均為波狀沙丘，沙丘均為淺丘，土壤類型均為風沙土，土壤侵蝕主要為風蝕。自恢復沉陷區的植被以自然恢復為主，植被類型為落葉闊葉灌叢和沙生類，覆蓋度30%～60%，主要物種包括沙柳（）、沙蒿（）、紫穗槐（）、檸條（）、沙米（）、針茅（）、狗尾草（()）等。植被修復沉陷區的修復措施為水平溝與魚鱗坑整地+杏樹（）種植，修復林地的郁閉度為0.6。本研究調查發現，未采區雖然還沒有發生地表沉陷，但已經受到了礦業活動的明顯影響，如：周邊地表沉陷對未采區地表-地下水文環境的影響、運煤公路對未采區的切割以及對其植被的破壞（植被覆蓋度50%～60%；未采區植被類型、主要物種與自恢復沉陷區一樣）等，因此，像沉陷區一樣，未采區也屬于采煤擾動區，反映地表沉陷之前的采煤擾動信息。另外，根據神東礦區井田分布范圍資料，本研究在神東礦區邊界外（距離上灣礦西部25 km左右）選擇植被良好（覆蓋度80%）且未受采煤擾動的區域作為研究區的對照區（CP，圖1a），以此反映礦區未受采煤擾動前的背景信息。而且，對照區的地貌類型、土壤類型、植被類型、主要植物物種與未采區、自恢復沉陷區完全一樣，從而確保對照區的代表性。對照區的供試樣地面積同樣為0.15 km2。

未采區、沉陷區、對照區土壤概況見表1。

表1 未采區、沉陷區、對照區土壤（0～15 cm）概況

表1中，對照區土壤容重小，主要是由于對照區植被情況良好，覆蓋度高達80%（遠高于未采區50%～60%與沉陷區30%～60%的植被覆蓋度），植物根系發達，因而導致對照區土壤疏松，結果表現出土壤容重小的現象。此外，對照區土壤容重小與土壤侵蝕強度有關。相比較未采區、自恢復沉陷區，對照區土壤侵蝕弱（見下文中的土壤侵蝕研究結果），因此，對照區土壤中保留的黏粒含量相對高（表1），所以這也是導致對照區土壤容重小的一個原因。而未采區、自恢復沉陷區土壤侵蝕強，由于黏粒相較于砂粒易受侵蝕（風蝕、水蝕）搬運，所以土壤中的黏粒含量降低、砂粒含量增加（表1），最終導致土壤容重增加。

137Cs示蹤法已被成功應用于風蝕區土壤侵蝕研究[14-16]。應用該示蹤技術開展土壤侵蝕評價的前提是獲取研究區域的137Cs背景值。本研究中，為利用137Cs示蹤法評價采煤擾動區內的未采區、沉陷區土壤侵蝕變化特征，在對照區附近選擇一處具有60 a歷史的林地（楊樹林，樣地面積為0.05 km2，圖1a）作為137Cs沉降本底區來獲取137Cs背景值。該本底區地表凋落物層次清晰，為不同分解階段的植物殘體，而且未發現泥沙混合現象，表明該土壤無人為擾動、無明顯侵蝕和堆積，因此是較理想的137Cs本底值采樣區。

137Cs示蹤土壤侵蝕的基本原理為：如果采樣區的137Cs濃度小于137Cs背景值，則該樣區被斷定為侵蝕區；反之，為沉積區。在侵蝕區，土壤侵蝕強度與137Cs濃度呈反比關系，即137Cs濃度越低，土壤侵蝕強度越大[17-18]。

1.3 樣品采集與分析

樣品采集方案為：首先，采集137Cs本底樣品。在137Cs本底區，設置7個137Cs背景點，用取土鉆（直徑=8 cm，鉆頭長度為15 cm；荷蘭Eijkelkamp公司生產）在每個樣點按照5 cm的土壤深度間隔采集至地表以下40 cm。其次，在上述未采區、5個沉陷區以及對照區分別隨機布設5個重復樣點，總計35個樣點。根據137Cs本底值的土壤深度分布特征（主要分布在0～25 cm），每個樣點按照15 cm的土壤深度間隔，分層采集0～15、>15～30 cm的土壤樣品。樣品采集結束，帶回實驗室測定。

待測樣品在實驗室自然風干、過篩，0～15 cm深度的土壤樣品用于測定土壤顆粒組成、容重、pH值、有機碳、全氮、全磷、堿解氮、有效磷、微生物量碳指標；0～30 cm深度的分層土壤樣品用于測定137Cs濃度與土壤容重。137Cs樣品的測定利用美國堪培拉（CANBERRA）公司生產的高純鍺（HPGe）探頭多道能譜儀。每個樣品測試質量不小于350 g，測試時間為40 000～60 000 s，測試結果相對誤差<5%。經儀器測定的137Cs含量為質量濃度（Bq/kg），結合土壤容重將其轉化為面積濃度（Bq/m2）。土壤顆粒組成、容重、pH、有機碳、全氮、全磷、堿解氮、有效磷采用常規分析方法來測定[19]，即：土壤顆粒組成通過吸管法來測定，然后把顆粒組分歸類成2～0.05 （砂粒）、0.05～0.002（粉粒）、<0.002 mm（黏粒）3部分，根據美國農部制土壤質地判定方法來劃定本研究中的土壤質地；土壤容重的測定是通過對本研究中按照5與15 cm土壤深度間隔采集的全部樣品利用烘干法（105 ℃）來進行（與環刀法原理一樣）；pH值采用電位法來測定（水土比為2.5:1）；有機碳用重鉻酸鉀油浴外加熱法測定；全氮與堿解氮分別用凱氏定氮法與堿解擴散法來測定；全磷用HClO4-H2SO4酸溶-鉬銻抗比色法測定，有效磷用碳酸氫鈉浸提-鉬銻抗比色法測定。土壤微生物量碳通過氯仿熏蒸－K2SO4提取法測定[20]。

1.4 數據分析

所有測定數據采用SPSS19.0統計軟件來處理。未采區、沉陷區以及對照區間的各指標均值差異采用ONE－ANOVA單因素方差分析來進行，在<0.05的顯著性水平下認為均值具有顯著性差異。

2 結果與討論

2.1 土壤侵蝕

本研究調查的137Cs本底值為726（SD:±34）Bq/m2。137Cs本底值深度分布格局（圖2）表明：隨著土壤深度的增加，137Cs面積濃度降低，且137Cs含量集中分布在0～25 cm的土壤深度范圍內，總體上呈現出單峰+負指數分布形態，這種格局符合137Cs在區域背景樣地土壤中的分布規律[21]。

圖2 137Cs本底值的土壤剖面分布格局

未采區、沉陷區、對照區（CP）的137Cs面積濃度（326～581 Bq/m2）均低于137Cs本底值（726 Bq/m2）（圖3），降幅為32%～55%，表明未采區、沉陷區與對照區的土壤遭受了嚴重的侵蝕，也反映了西北風積沙區風蝕嚴重的背景特征。

Note: CP, control plot; UEP, unexploited plot; 1ySP, 1-year subsided plot; 2ySP, 2-year subsided plot; 4ySP, 4-year subsided plot; 8ySP, 8-year subsided plot; 13ySP, 13-year subsided plot, the same below.

對照區137Cs含量（494 Bq/m2）反映研究區未受采煤擾動下的土壤侵蝕水平；與對照區相比，未采區137Cs含量（454 Bq/m2）降低8%，表明采煤擾動加劇了未采區土壤侵蝕，也反映出采煤擾動影響下礦區地表沉陷之前的土壤侵蝕變化。未采區土壤侵蝕增強可能與井下開采導致的地下水水文格局破壞有關。由于研究區煤系上覆巖層之上的第四系河谷沖積層潛水為土壤、植被的涵養水層，井下工作面開采產生的導水裂縫帶導通了第四系潛水含水層，使得采空區形成地下漏斗，從而誘發未采區的地下水向采空區匯集，進而加劇未采區地表土壤干旱、沙化，最終導致未采區土壤侵蝕強度增加。

與對照區相比，沉陷歷史為1、2、4、8 a的自恢復沉陷區137Cs面積濃度在未采區8%的降幅基礎上繼續降低（圖3），分別降低16%、29%、34%、10%。這表明：開采沉陷進一步加劇了研究區的土壤侵蝕強度；而且，自地表沉陷后，隨著時間的推移，土壤侵蝕強度呈現出先增加后降低的趨勢。綜合未采區、1、2、4、8 a沉陷區137Cs含量變化，也可以看出：地表沉陷后的最初2 a，土壤侵蝕急劇增大；4 a以后，土壤侵蝕呈現降低趨勢。該結果可能歸因于植被的破壞與恢復。實地調查發現，高強度開采沉陷后的最初2 a，由于地表裂縫多（尤其是垂直于井下工作面推進方向的動態裂縫），灌喬木受裂縫拉斷根系影響而枯萎，禾本、草本科植物因裂縫增加了土壤水分蒸發通道而大面積干旱死亡，該時段內植被對土壤的保持作用極差，因此，土壤侵蝕急劇增大；及至第4年，裂縫多被風沙掩埋，植被開始好轉，植被固土作用也相應轉強，因此該時段后，土壤侵蝕開始減弱。從1、2、4、8a沉陷區都是自恢復沉陷區這個角度來看，8年沉陷區相較于1、2、4a沉陷區，其137Cs含量呈現增加趨勢。這個結果說明：西北風積沙區沉陷地表隨著自恢復期限的延長，土壤侵蝕水平可以降低。

植被修復沉陷區由于植被修復期長（12 a）且取得了良好的修復效果，137Cs面積濃度相較于對照區增加（圖3），增幅為18%，表現出最低的土壤侵蝕強度，由此可見，植被修復是防治西北風積沙區采煤沉陷地表土壤侵蝕的有效手段。

2.2 土壤養分

圖4顯示了未采區、沉陷區、對照區之間土壤養分變化。與對照區相比，未采區土壤有機碳、全氮、全磷、堿解氮含量顯著降低（0.05），降低幅度分別為55%、64%、29%、36%，表明地表沉陷之前這些養分的明顯損失。但是，未采區土壤有效磷含量相較于對照區顯著增加（0.05），這可能是因為井下開采對土壤的涵養水層（第四系潛水）的破壞加劇了土壤干旱，從而增強了土壤有機磷向無機磷的礦化。

圖4 未采區、沉陷區、對照區的土壤養分含量

對于自恢復沉陷區而言，土壤有機碳、全氮、堿解氮含量隨著時間的推移而增加（圖4），該結果與劉哲榮等[22-23]的研究結果一致，表明西北風積沙區開采沉陷后土壤養分的自恢復演變趨勢。具體來說，1、2 a沉陷區土壤有機碳、全氮、堿解氮含量最低，且低于未采區（0.05）；4、8 a沉陷區這3個養分的含量分別接近（0.05）、顯著高于未采區，但仍低于對照區。這表明：地表沉陷后的最初2 a，土壤有機碳、全氮、堿解氮損失最嚴重；4~8 a后，它們的含量將恢復至甚或超出地表沉陷前的水平，但很難恢復到采煤擾動前的水平。另外，土壤全磷含量僅在4 a沉陷區顯著低于未采區（<0.05），但相較于其他3個自恢復沉陷區（1、2、8 a）無顯著差異（>0.05），這表明：地表沉陷后的第4年，全磷損失明顯；其后隨著時間的推移，全磷含量可以恢復到地表沉陷前的水平。土壤有效磷在1、2、4、8 a沉陷區的含量均低于未采區（<0.05），大體上表現出隨著時間的推移而逐漸降低至采煤擾動前水平的趨勢。該結果可能歸因于地表沉陷后，隨著植被自恢復效果的逐漸好轉，植被對有效磷養分的利用不斷增強。此外，這4個自恢復沉陷區土壤有效磷含量均不同程度地高于對照區。土壤有效磷虧缺被證實是限制許多礦區土壤、植被正向演替的影響因子[24]。但是，本文結果表明本研究區并沒有發生因采煤擾動導致的土壤有效磷虧缺，因此本礦區的生態修復應該合理利用土壤有效磷本底優勢，促進土壤、植被正向演替。

在植被修復沉陷區（13 a沉陷區），土壤有機碳、堿解氮含量顯著高于對照區（<0.05），全氮、全磷、有效磷含量接近對照區（>0.05），這表明植被修復可以有效提升西北風積沙區開采沉陷地表退化土壤的養分。

本研究調查了采煤擾動區土壤微生物量碳的變化。土壤微生物量碳作為土壤有機碳庫中最為活躍的組分，盡管只占總有機碳的1%～5%，但是它積極參與有機質的分解、養分的循環與轉化，是反映土壤微生物活性強弱的最重要的指標[25]。本調查（圖4）發現，采煤擾動區土壤微生物量碳的演變趨勢類似于土壤有機碳、全氮、堿解氮。與土壤有機碳、全氮、堿解氮、全磷含量在植被修復沉陷區高于或接近對照區不同的是，土壤微生物量碳含量在植被修復沉陷區低于（<0.05）對照區，表明植被修復對沉陷區土壤微生物活性的改善效果弱于有機碳與氮磷養分。因此，研究區生態修復應兼顧考慮植被修復與土壤微生物修復措施。

2.3 微生物熵、碳氮比、碳磷比

為進一步探討采煤擾動區土壤養分的動態，本研究分析了土壤養分比值指標：微生物熵（土壤微生物量碳與有機碳的比值；下文用MBC/SOC來表示）、碳氮比（C/N）與碳磷比（C/P）。圖5顯示：與對照區（3.43%）相比，未采區、沉陷區的土壤MBC/SOC均出現了不同程度的降低（1.87%～3.18%），其中，8 a自恢復沉陷區與植被修復沉陷區的MBC/SOC最低。與MBC/SOC的變化相反，未采區、沉陷區的C/N（10.4～15.2）相較于對照區呈現不同程度的增加，其中，未采區、8 a自恢復沉陷區與植被修復沉陷區的C/N最高。相較于對照區，未采區的C/P降低，植被修復沉陷區的C/P增加；但在自恢復沉陷區，C/P有增有減，其中，C/P在1、2、4 a自恢復沉陷區降低，而在8 a自恢復沉陷區增加。

圖5 未采區、沉陷區、對照區土壤微生物熵、碳氮比及碳磷比

土壤微生物熵是衡量土壤有機質積累水平的重要指標，它比單獨用有機碳或微生物量碳能更好地反映有機質質量的變化[26]。土壤C/N、C/P通常被用來表征土壤氮、磷素的礦化能力[27-29]。土壤C/N、C/P低，有利于促進微生物分解有機質釋放養分，增加土壤有效氮、磷供給；反之，會導致微生物在分解有機質的過程中氮、磷受限，從而產生微生物與植物對土壤無機氮、磷的競爭，不利于植物的生長。本研究土壤MBC/SOC變化表明，土壤有機質積累水平在地表沉陷前與沉陷后可以出現不同程度的降低，尤其是在沉陷后。8 a自恢復沉陷區與13 a植被修復沉陷區表現出最低的MBC/SOC，該結果歸因于這兩個沉陷歷史長的沉陷區土壤微生物活性恢復程度遠低于有機碳。因此，開展礦區土壤有機質質量評價應該綜合考慮土壤有機碳與微生物量碳。土壤C/N結果表明：地表沉陷前與沉陷后土壤氮素礦化能力降低，土壤有效氮供應不足。土壤C/P變化表明：受采煤擾動影響，礦區地表在其沉陷前及沉陷后的短期內（1～4 a），土壤磷素礦化能力增強。前述結果“未采區、1、2、4 a沉陷區土壤有效磷含量最高（見表2）”也證實了礦區地表沉陷前及沉陷后短期內的土壤磷素礦化能力增強這一結果。8 a自恢復沉陷區與植被修復沉陷區C/N、C/P最高，表明隨著時間推移，土壤氮磷素礦化能力降低的趨勢，這可能是因為該階段（地表沉陷后的8～13 a）沉陷區植被恢復效果轉好，進入土壤中的有機物（枯枝落葉）迅速增多，但土壤微生物活性未能同步增強，所以限制了微生物對有機N、P底物的礦化分解。

2.4 土壤養分與侵蝕的聯系

在山地、丘陵等地貌區，土壤侵蝕是導致土壤養分損失的主要原因[9,30]。通過2.1與2.2節結果可以看出：與對照區相比，未采區、1、2、4、8 a自恢復沉陷區、植被修復區的土壤有機碳、全氮、全磷、微生物量碳、堿解氮變化相似于137Cs變化，表明西北風積沙區采煤擾動地表土壤有機碳及氮、磷養分的演變規律與土壤侵蝕演變之間存在密切聯系。相關分析也顯示：土壤有機碳、全氮、全磷、微生物量碳、堿解氮含量均與137Cs含量呈正相關關系（0.76～0.87，<0.05；見圖6），這個結果與程靜霞等[31-32]的研究結果一致，表明礦區土壤有機碳及氮、磷養分含量隨采煤擾動下土壤侵蝕強度的增加而降低。因此，土壤侵蝕防治是西北風積沙區采煤擾動地表土壤養分保持的首要任務。

圖6 土壤養分與137Cs之間的相關關系

3 結 論

本研究對西北風積沙礦區采煤擾動下土壤侵蝕與養分演變特征調查的結果表明：

1）受采煤擾動影響，未采區土壤侵蝕強度增大，土壤有機碳、微生物量碳、全氮、全磷、堿解氮養分損失明顯，土壤有機質積累水平、氮素礦化能力降低，土壤磷素礦化能力增加。

2）地表沉陷后，隨著時間的推移，土壤侵蝕強度呈現出先增加后降低的趨勢，土壤有機碳、微生物量碳、全氮、全磷、堿解氮含量呈現出先降低后增加的趨勢。地表沉陷后的最初2 a是土壤侵蝕急劇增大、碳氮養分損失最嚴重的時段。

3）采煤擾動下，礦區土壤有機碳、微生物量碳、全氮、全磷、堿解氮養分演變與土壤侵蝕演變密切相關。

4）植被修復可以有效降低采煤沉陷區土壤侵蝕強度與提高土壤有機碳與氮磷養分含量。

基于本研究結果，西北風積沙礦區的生態恢復應該以侵蝕防治為首要任務，合理利用土壤有效磷本底優勢，兼顧考慮植被修復與土壤微生物修復，以此促進土壤、植被正向演替。此外，本研究也啟示研究者在今后開展礦區生態環境演變研究時，如果把未采區作為對照區，可能存在不能準確評價沉陷區土壤、植被、水等生態環境要素變化的風險。

[1] 王銳，馬守臣，張合兵，等.干旱區高強度開采地表裂縫對土壤微生物學特性和植物群落的影響[J]. 環境科學研究，2016，29(9)：1249－1255.

Wang Rui, Ma Shouchen, Zhang Hebing, et al. Effects of surface cracks caused by high intensity coal mining on soil microbial characteristics and plant communities in arid regions[J]. Research of Environmental Sciences, 2016, 29(9): 1249－1255. (in Chinese with English abstract)

[2] 趙紅梅，張發旺，宋亞新，等. 大柳塔采煤塌陷區土壤含水量的空間變異特征分析[J]. 地球信息科學學報，2010，12(6)：753－760.

Zhao Hongmei, Zhang Fawang, Song Yaxin, et al. Spatial variation of soil moisture content in mining subsidence areas of Da Liuta, Shenmu County, Shanxi Province[J]. Journal of Geo-information Science, 2010,12(6): 753－760. (in Chinese with English abstract)

[3] 臧蔭桐，汪季，丁國棟，等.采煤沉陷后風沙土理化性質變化及其評價研究[J].土壤學報，2010，47(2)：262－269.

Zang Yintong, Wang Ji, Ding Guodong, et al. Sand soil physical and chemical properties change after mining subsidence and its evaluation research[J]. Acta Pedologica Sinica, 2010, 47(2): 262－269. (in Chinese with English abstract)

[4] 卞正富，雷少剛，常魯群，等. 基于遙感影像的荒漠化礦區土壤含水率的影響因素分析[J]. 煤炭學報，2009，34(4)：520－525.

Bian Zhengfu, Lei Shaogang, Chang Luqun, et al. Based on the influence factors of mining area desertification soil moisture content of remote sensing imag analysis[J]. Journal of China Coal Society, 2009, 34(4): 520－525. (in Chinese with English abstract)

[5] 魏江生，賀曉，胡春元，等.干旱半干旱地區采煤塌陷對沙質土壤水分特性的影響[J].干旱區資源與環境，2006，20(5)：84－88.

Wei Jiangsheng, He Xiao, Hu Chunyuan, et al. Influence of ground collapse caused by coal mining activities on the water characteristics of sandy soil in arid and semi-arid area[J]. Journal of arid Land Resources and Environment, 2006, 20(5): 84－88. (in Chinese with English abstract)

[6] 王健，高永，魏江生，等.采煤塌陷對風沙區土壤理化性質影響的研究[J].水土保持學報，2006，20(5)：52－55.

Wang Jian, Gao Yong, Wei Jiangsheng, et al. Research on the effects of coal mining subsidence on sandy soil physical and chemical properties[J]. Journal of Soil and Water Conservation, 2006, 20(5): 52－55. (in Chinese with English abstract)

[7] 臺曉麗，胡振琪，陳超. 西部風沙區不同采煤沉陷區位土壤水分中子儀監測[J]. 農業工程學報，2016，32(15)：225－231.

Tai Xiaoli, Hu Zhenqi, Chen Chao. Monitoring soil moisture at different subsidence areas of mining in western windy and sandy region with neutron instrument[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(15): 225－231. (in Chinese with English abstract)

[8] 鄒慧，畢銀麗，朱郴韋，等. 采煤沉陷對沙地土壤水分分布的影響[J]. 中國礦業大學學報，2014，43(3)：496－501.

Zou Hui, Bi Yinli, Zhu Chenwei, et al. Effect of mining subsidence on soil moisture dynamic changes of sandy land[J]. Journal of China University of Mining & Technology, 2014, 43(3): 496－501. (in Chinese with English abstract)

[9] Oldeman L R.The Global Extent of Soil Degradation. In: Soil Resilience and Sustainable Land Use (eds Greenland D J, Szabolcs I)[M]. CAB International, Wallingford, UK, 1994, 99－118.

[10] Li Sheng, Lobb D A, Kachanoski R G, et al. Comparing the use of the traditional and repeated sampling-approach of the137Cs technique in soil erosion estimation[J]. Geoderma, 2011, 160: 324－335.

[11] Nie Xiaojun, Zhang Jianhui, Cheng Jingxia, et al. Effect of soil redistribution on various organic carbons in a water- and tillage-eroded soil[J]. Soil & Tillage Research, 2016, 155(1): 1－8.

[12] Martinez C, Hancock R, Kalma J D. Relationships between137Cs and soil organic carbon (SOC) in cultivated and never-cultivated soils: An Australian example[J]. Geoderma, 2010, 158(3): 137－147.

[13] Ni Shijun, Zhang Jianhui. Variation of chemical properties as affected by soil erosion on hillslopes and terraces[J]. European Journal of Soil Science, 2007, 58(6): 1285－1292.

[14] 嚴平，董光榮，張信寶，等.137Cs 法測定青藏高原土壤風蝕的初步結果[J]. 科學通報，2000，45(2)：199－204.

[15] 張春來，鄒學勇.干草原地區土壤137Cs沉積特征[J].科學通報，2002，47(3)：221－225.

[16] 胡云鋒，劉紀遠，莊大方，等. 風蝕土壤剖面137Cs的分布及侵蝕速率的估算[J]. 科學通報，2005，50(9)：933－937.

[17] Van Pelt R S, Zobeck T M, Ritchie J C, et al. Validating the use of Cs-137 measurements to estimate rates of soil redistribution by wind[J]. Catena, 2007, 70(3): 455－464.

[18] Walling D E, Quine T A. Use of caesium-137 as a tracer of erosion and sedimentation: Handbook for the application of the caesium-137 technique[J]. Estudio Fao Montes,1993.

[19] 張甘霖，龔子同.土壤調查實驗室分析方法[M]. 北京：科學出版社，2012：8－49.

[20] Vance E D, Brookes P C, Jenkinson D S. An extraction method for measuring soil microbial biomass C[J]. Soil Biology and Biochemistry, 1987, 19(6): 703－707.

[21] Walling D E, Quine T A. Calibration of caesium-137 measurements to provide quantitative erosion rate data[J]. Land Degradation and Rehabilitation, 1990, 2(3): 161－175.

[22] 劉哲榮，燕玲，賀曉，等. 采煤沉陷干擾下土壤理化性質的演變-以大柳塔礦采區為例[J]. 干旱區資源與環境，2014，28(4)：133－138.

Liu Zherong, Yan Ling, He Xiao, et al. Effects of mining subsidence on physical and chemical properties of soil in the subsided land of the Daliuta Mining Area[J]. Journal of Arid Land Resources and Environment, 2014, 28(4): 133－138. (in Chinese with English abstract)

[23] 王琦，全占軍，韓煜，等. 采煤塌陷對風沙區土壤性質的影響[J]. 中國水土保持科學，2013，11(6)：110－118.

Wang Qi, Quan Zhanjun, Han Yu, et al. Effects of mining subsidence on soil properties in windy desert area[J]. Science of Soil and Water Conservation, 2013, 11(6): 110－118. (in Chinese with English abstract)

[24] Chen H, Zhang C, Zhu Y. Phosphorus: A limiting factor for restoration of soil fertility in a newly reclaimed coal mined site in Xuzhou, China[J]. Land Degradation & Development, 1998, 9(2): 115－121.

[25] Marinari S, Mancinelli R, Campiglia E, et al. Chemical and biological indicators of soil quality in organic and conventional farming systems in Central Italy[J]. Ecological Indictors, 2006, 6(4): 701－711.

[26] Anderson TH, Domsch KH. Ratios of microbial biomass carbon to total organic carbon in arable soils[J]. Soil Biology and Biochemistry, 1989, 21(4): 471－479.

[27] Wang J L, Zhong Z M, Wang Z H, et al. Soil C/P distribution characteristics of alpine steppe ecosystems in the Qinhai-Tibetan Plateau[J]. Acta Prataculturae Sinica, 2014, 23(2): 9－19.

[28] Wiesmeier M, Hübner R, Barthold F, et al. Amount, distribution and driving factors of soil organic carbon and nitrogen in cropland and grassland soils of southeast Germany (Bavaria)[J]. Agriculture, Ecosystems and Environment, 2013, 176: 39－52.

[29] Luo Y Q, Currie W S, Dukes J S, et al. Progressive nitrogen limitation of ecosystem responses to rising atmospheric carbon dioxide[J]. Bioscience, 2004, 54(8): 731－739.

[30] Yan Yuchun, Xin Xiaoping, Xu Xingliang, et al. Quantitative effects of wind erosion on the soil texture and soil nutrients under different vegetation coverage in a semiarid steppe of northern China[J]. Plant & Soil, 2013, 369(1): 585－598.

[31] 程靜霞，聶小軍，劉昌華. 煤炭開采沉陷區土壤有機碳空間變化[J]. 煤炭學報，2014，39(12)：2495－2500.

Cheng Jingxia, Nie Xiaojun, Liu Changhua. Soil organic carbon spatial change of coal mining subsidence area[J]. Journal of China Coal Society,2014, 39(12): 2495－2500. (in Chinese with English abstract)

[32] 張合兵，聶小軍，程靜霞.137Cs 示蹤采煤沉陷坡土壤侵蝕及其對土壤養分的影響[J].農業工程學報，2015，31(4)：137－143.

Zhang Hebing, Nie Xiaojun, Cheng Jingxia.137Cs tracing of soil erosion and its impact on soil nutrients across subsidence slope induced by coal mining[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(4): 137－143. (in Chinese with English abstract)

Characteristics of soil erosion and nutrients evolution under coal mining disturbance in aeolian sand area of Northwest China

Nie Xiaojun1, Gao Shuang1, Chen Yongliang2, Zhang Hebing1

(,,454000,; 2.,,719315,)

Deterioration of the local ecological environment due to coal mining is serious in the aeolian sand area, Northwest China. A better understanding of the evolution of aeolian sandy soil under coal mining disturbance is essential prior to the initiation of restoration of local ecological environment. In this study, the137Cs tracing technique was used to investigate the characteristics of soil erosion and nutrients in coal mining disturbed lands of Shendong mining area, Northwest China. A plot without being disturbed by coal mining activities near Shendong mining area was selected as the control plot. As for the mining disturbed lands, an unexploited plot and 5 subsidence plots were selected as the study sites. The 5 subsidence plots include 4 self-recovery subsidence plots and one vegetation restoration plot. The 4 self-recovery subsidence plots (1ySP, 2ySP, 4ySP, 8ySP) have undergone stabilization for 1, 2, 4, and 8 years, respectively after surface subsidence. The vegetation restoration plot has experienced a 13-year duration of stabilization after surface subsidence, and has been planted with almond trees for 12 years. It was found that137Cs inventories from the unexploited plot, subsidence plots and control plot were 32%-55% lower than local137Cs reference inventory (726 Bq/m2), which demonstrates the background of severe wind erosion desertification in the study area. In comparison with the control plot, the137Cs inventories in the unexploited plot, 1ySP, 2ySP, 3ySP and 4ySP decreased by 8%, 16%, 29%, 34%, and 10%, respectively. Changes in the nutrients, including soil organic carbon (SOC), microbial biomass carbon (MBC), total nitrogen (TN), available nitrogen (AN) and total phosphorus (TP), showed a similar trend to the137Cs variation. Unlike SOC, MBC, TN, AN, and TP, available phosphorus showed higher concentrations in the unexploited plot and the 4 self-recovery subsidence plots than that in the control plot. In addition, the137Cs, SOC and TN concentrations increased in the vegetation restoration plot in comparison with that in the control plot. Overall, there were significant correlations between SOC, MBC, TN, AN, TP and137Cs within the mining disturbed lands. It was also found that the coal-mining disturbed lands exhibited smaller ratios of MBC to SOC but greater ratios of SOC to TN (i.e. C/N) than the control plot. These results suggested that the intensity of soil erosion firstly increased and then decreased after surface subsidence in the aeolian sand area. Soil erosion increased sharply in the first 2 years after surface subsidence. The evolution of soil nutrients is closely associated with soil erosion evolution under coal mining disturbance. The present study provides evidences for the acceleration of soil erosion and nutrients depletion in unexploited lands of mining area by coal mining disturbance. As a consequence, there would be a considerable underestimation of the impact of mining subsidence on soil evolution when using unexploited plots as the control. Our study also suggested that vegetation restoration may effectively weaken soil erosion and improve soil nutrients in subsidence lands. Therefore, ecological restoration in subsidence lands of the aeolian sand area, Northwest China, should combat against earlier soil erosion after surface subsidence, and a plant-microorganism integrated remediation approach should be adopted to promote the positive succession of soil and vegetation.

soils; erosion; nutrients; mining subsidence; ecological environment evolution;137Cs tracing method; aeolian sand area

10.11975/j.issn.1002-6819.2018.02.017

S157.1; P208

A

1002-6819(2018)-02-0127-08

2017-08-31

2017-12-06

國家自然科學基金委員會與神華集團有限責任公司聯合資助項目（U1261206）；河南省高校創新團隊支持計劃（18IRTSTHN008）

聶小軍，男，山西曲沃人，博士，副教授，主要從事土壤侵蝕與環境生態研究。Email：niexj2005@126.com

聶小軍，高 爽，陳永亮，張合兵. 西北風積沙區采煤擾動下土壤侵蝕與養分演變特征[J]. 農業工程學報，2018，34(2)：127－134. doi：10.11975/j.issn.1002-6819.2018.02.017 http://www.tcsae.org

Nie Xiaojun, Gao Shuang, Chen Yongliang, Zhang Hebing. Characteristics of soil erosion and nutrients evolution under coal mining disturbance in aeolian sand area of Northwest China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(2): 127－134. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2018.02.017 http://www.tcsae.org