黃河中游煤礦區沉陷黃土坡面土壤養分空間變化規律研究

劉 具，宋世杰，杜 麟，楊 磊，梁躍強，程 坤

(1.中國煤炭工業協會咨詢中心，北京市朝陽區，100013；2. 西安科技大學地質與環境學院，陜西省西安市，710054)

煤炭作為我國的主體能源，長時期內對保障國家能源安全發揮著重大作用[1]。黃河流域不僅是我國重要生態屏障，更是煤炭資源主產區，素有“能源流域”之稱，我國煤炭資源的開采主要集中在黃河中上游，這些區域大多處于干旱、半干旱地區，生態環境脆弱[2-4]。煤炭開采不可避免地引發生態環境損傷，加劇水土沉陷流失、耕地損失、植被退化。有學者曾提出自然與人工干預相結合的修復方式[5]，并指出采煤擾動下植被修復可有效降低土壤有機碳及氮磷養分含量[6-9]。雖然開采沉陷及其衍生的生態環境損害復雜多樣，但土壤質量損傷是生態環境負效應的關鍵要素，而有機質、銨態氮、磷和鉀等養分含量是表征土壤質量的核心指標，且極大地制約著地表植被的生長狀況和修復方式。在治理開采沉陷過程中，黃土坡面會產生各種變形，導致坡面土壤養分產生顯著空間異質性[10-12]。

本文以黃河中游煤礦區陜北張家峁礦區井田為研究對象，主要研究陜北煤礦區沉陷黃土坡面土壤養分的空間變化特征與規律，用于指導沉陷區生態修復工作，具有重要的科學價值和實踐意義。

1 研究區概況

研究區位于陜西省榆林市神木縣北部張家峁井田，地處陜西省榆林市神木縣北部，屬典型中溫帶半干旱大陸性氣候，冬季寒冷、夏季炎熱，晝夜溫差懸殊，多年平均氣溫為8.5℃，最高氣溫為38.9℃，最低氣溫為-28.5℃，平均降水量為436.7 mm，多集中于7～9月，平均蒸發量為1 907～2 122 mm，是降水量的4～5倍，多年平均風速為2.5 m/s、最大風速為25 m/s，凍土最大深度為146 cm。研究區除西南角外均為黃土溝壑地貌，地表為第四系中上更新統黃土，厚度一般為50～100 m，坎陡溝深、梁峁相間。井田全區主采煤層為3-1煤層、埋深為0～108 m、平均厚度為2.84 m，在長壁綜采開采方式下開采沉陷發育明顯，地表下沉系數約為0.7。

2 材料與方法

2.1 布點與采樣

在研究區選擇沉陷穩定(2 a左右)的3處形態、大小相近的典型黃土坡面，分別在坡頂、坡中和坡腳的3個部位隨機布設3個1 m×1 m的采樣方，采用“五點法”分別采集0～10、10～20、20～40、40～60 cm等4個深度的土壤；將每個黃土坡面不同部位、相同土層深度土壤充分混合后，分別裝入純棉布袋；將同區域未沉陷且形態相似的黃土坡面作為對照，按相同方法采集土壤樣品；共采集48個樣品。研究區沉陷黃土坡面土壤樣品采集示意圖如圖1所示。

圖1 研究區沉陷黃土坡面土壤樣品采集示意圖

2.2 樣品處理與檢測方法

樣品經晾曬、除雜、研磨、過篩等處理，采用總有機碳分析儀(Vario TOC)測定有機質含量，全自動流動分析儀(AA3)測定銨態氮含量，電感耦合等離子體(PerkinElme2000，ICP)測定速效鉀含量，可見分光光度計測定速效磷含量，每項平行測定3次，取算數平均值作為該項測試的最終值。

2.3 數據處理

實驗數據采用SPSS 21.0軟件進行顯著性分析和相關性分析，所得圖件采用Origin 9.1軟件繪制。

3 實驗與數據分析

3.1 室內實驗

通過室內實驗，測定了土壤有機質、銨態氮、速效磷、速效鉀等養分指標，陜北煤礦區沉陷黃土坡面土壤養分含量測定結果見表1。

表1 陜北煤礦區沉陷黃土坡面土壤養分含量測定結果

3.2 數據分析

根據檢測結果繪制了不同部位及土層深度下各土壤養分含量情況，陜北煤礦區沉陷黃土坡面不同部位及土層深度的土壤養分對比如圖2所示。

3.2.1 沉陷黃土坡面土壤有機質空間變化規律

由圖2(a)可以看出：

(1)沉陷黃土坡面土壤有機質含量在不同土層垂直深度上存在明顯差異，坡頂0～10 cm土壤有機質含量為550 mg/kg，是40～60 cm土壤的4.6倍；坡中0～10 cm土壤有機質含量為730 mg/kg，是40～60 cm土壤的5.6倍；坡腳10～20 cm土壤有機質含量為542 mg/kg，分別是0～10 cm、40～60 cm土壤的1.6倍和1.4倍。對于坡頂和坡中，沉陷黃土坡面土壤有機質含量均呈現隨土層深度增加而逐漸降低的變化趨勢，且在坡中最明顯，有機質含量主要集中于0～20 cm土層土壤中；對于坡腳，有機質含量隨土層垂直深度增加而表現為先增加后下降的變化趨勢，有機質含量主要集中于10～40 cm土層土壤。

(2)與對照組相比較，沉陷黃土坡面任何部位和任何土層垂直深度中的土壤有機質含量均有所損失。從不同坡面部位角度看，垂直深度為0～60 cm以內的不同坡面部位土壤有機質平均損失率由大到小的排序為坡頂(55.2%)>坡中(26.6%)>坡腳(0.4%)；從不同土層深度角度看，3個坡面部位的各層土壤有機質平均損失率由大到小的排序為20～40 cm(41.0%)>40～60 cm(36.9%)>0～10 cm(29.8%)>10～20 cm(28.8%)。由此可見，開采沉陷造成了全坡面土壤有機質含量的損失，而坡頂20～40 cm土壤損失最為嚴重。

圖2 陜北煤礦區沉陷黃土坡面不同部位及土層深度的土壤養分對比

3.2.2 沉陷黃土坡面土壤銨態氮的空間變化規律

由圖2(b)可以看出：

(1)沉陷黃土坡面土壤銨態氮含量在不同土層垂直深度上存在明顯差異，坡頂10～20 cm、40～60 cm土壤銨態氮含量為3.66 mg/kg、3.43 mg/kg，分別是其它2層土壤平均含量的1.8倍和1.7倍；坡中10～20 cm土壤銨態氮含量為3.49 mg/kg，分別是0～10 cm、40～60 cm土壤的1.5倍和1.4倍；坡腳0～10 cm土壤銨態氮含量為4.11 mg/kg，是40～60 cm土壤的1.9倍。由此可見，坡頂銨態氮含量主要集中于10～20 cm及40～60 cm土層土壤中，坡中銨態氮含量主要集中于10～20 cm土層土壤中，坡腳銨態氮含量主要集中于0～20 cm土層土壤中。

(2)與對照組相比較，沉陷黃土坡面任何部位和任何土層垂直深度中的土壤銨態氮含量均有所損失。具體而言，從不同坡面部位角度看，垂直深度0～60 cm以內的不同坡面部位土壤銨態氮平均損失率由大到小的排序為坡腳(47.1%)>坡頂(37.7%)>坡中(15.1%)；從不同土層深度角度看，3個坡面部位的各層土壤銨態氮平均損失率由大到小的排序為0～10 cm土層(52.0%)>20～40 cm土層(41.5%)>10～20 cm土層(28.1%)>40～60 cm土層(13.2%)。由此可見，開采沉陷造成了全坡面土壤銨態氮含量的損失，而坡腳0～10 cm土層土壤損失最為嚴重。

3.2.3 沉陷黃土坡面土壤速效磷的空間變化規律

由圖2(c)可以看出：

(1)沉陷黃土坡面土壤速效磷含量在不同部位上存在明顯差異，坡腳0～10 cm土壤速效磷平均含量為270 mg/kg，均為坡頂、坡中的1.4倍；坡腳10～60 cm土壤速效磷平均含量為253 mg/kg，分別是坡頂、坡中的1.2倍和1.6倍。由此可見，在坡頂向坡中，再向坡腳轉換過程中，土層深度10～60 cm土壤速效磷含量呈現出先降低后增加的變化趨勢，且速效磷含量主要集中于坡腳部位。

(2)沉陷黃土坡面土壤速效磷含量在各個部位上不同土層垂直深度的差異不明顯，坡頂0～10 cm土壤速效磷含量為190 mg/kg，與其它各層土壤含量相當；坡中0～10 cm土壤速效磷含量為200 mg/kg，僅是40～60 cm土壤的1.2倍；坡腳0～10 cm土壤速效磷含量為270 mg/kg，僅是40～60 cm土壤的1.3倍。由此可見，無論在坡頂、坡中和坡腳，沉陷黃土坡面土壤速效磷含量在各土層上的分布較為均勻。

(3)與對照組相比較，沉陷黃土坡面的任何部位和任何土層垂直深度上，土壤速效磷含量均有所損失。具體而言，從不同坡面部位角度看，垂直深度0～60 cm以內的不同坡面部位土壤速效磷平均損失率由大到小的排序為坡中(32.7%)>坡頂(25.7%)>坡腳(11.2%)；從不同土層深度角度看，3個坡面部位的各層土壤有機質平均損失率由大到小的排序為40～60 cm土層(33.9%)>20～40 cm土層(24.2%)>10～20 cm土層(17.1%)>0～10 cm土層(14.4%)。由此可見，開采沉陷造成了全坡面土壤速效磷含量的損失，而坡中40～60 cm土壤損失最為嚴重。

3.2.4 沉陷黃土坡面土壤速效鉀的空間變化規律

由圖2(d)可以看出：

(1)沉陷黃土坡面土壤速效鉀含量在不同部位上存在明顯差異，坡腳土層深度0～10 cm土壤速效鉀含量為242.8 mg/kg，分別是坡頂、坡中的1.3倍和1.5倍；坡中土層深度10～20 cm土壤速效鉀含量為206.4 mg/kg，分別是坡頂、坡腳的1.2倍和1.1倍；坡頂土層深度20～60 cm土壤速效鉀平均含量為263.4 mg/kg，分別是坡中、坡腳的1.1倍和1.4倍。由此可見，在坡頂向坡中，再向坡腳轉換過程中，土層深度10～20 cm土壤速效鉀含量呈現出先增加后降低變化趨勢，且速效鉀含量主要集中于坡中部位；土層深度20～60 cm土壤速效鉀主要集中于坡頂部位。

(2)沉陷黃土坡面土壤速效鉀含量在不同土層垂直深度上存在明顯差異，坡頂20～40 cm土壤速效鉀含量最高，為271.2 mg/kg，10～20 cm土壤速效鉀含量最低，為172.3 mg/kg，前者是后者的1.6倍；坡中40～50 cm土壤速效鉀含量為249.9 mg/kg，是0～10 cm土壤的1.5倍；坡腳0～10 cm土壤速效鉀含量最高，為242.8 mg/kg，10～20 cm土壤速效鉀含量最低，為197.5 mg/kg，前者是后者的1.2倍。由此可見，對于坡頂和坡腳，沉陷黃土坡面速效鉀含量主要集中于坡頂40～60 cm、坡腳0～10cm土層土壤中。

(3)與對照組相比較，除20～40 cm土層及坡中40～60 cm土層外，沉陷黃土坡面土壤速效鉀含量均有所損失。具體而言，從不同坡面部位角度看，垂直深度0～60 cm以內的不同坡面部位土壤速效鉀平均損失率由大到小的排序為坡頂(7.3%)>坡腳(6.4%)>坡中(-1.6%)；從不同土層深度角度看，3個坡面部位的各層土壤速效鉀平均損失率由大到小的排序0～10 cm(14.6%)>10～20 cm土層(13.3%)>40～60 cm土層(-1.2%)>20～40 cm土層(-11.0%)。由此可見，開采沉陷造成了全坡面0～20 cm土層、坡頂及坡腳40～60 cm土層土壤速效鉀含量的損失，其中以全坡面0～10 cm土層損失最嚴重。

3.3 沉陷黃土坡面土壤養分修復策略

基于實驗及數據分析，陜北煤礦區開采沉陷發生后，會造成沉陷黃土坡面土壤有機質、銨態氮、速效磷等養分全坡面性損失以及速效鉀局部性損失和局部性富集[10，12]。因此，修復陜北采煤沉陷區生態環境，應根據區域生態環境質量定位以及修復后土地利用類型與植被類型，有選擇、有目標開展沉陷黃土坡面土壤養分修復。具體策略包括以下3個方面。

(1)采取人工修復與自然恢復相結合的方式實現全坡面土壤有機質、銨態氮、速效磷、速效鉀(0～20 cm土層)的修復。

(2)黃土沉陷坡面的坡頂、坡中土壤養分最為貧瘠，自然恢復潛力小，應以人工修復為主，而坡腳土壤養分較為豐富，自然恢復潛力大，應以自然恢復為主。

(3)坡頂20～40 cm土層、坡腳0～10 cm土層、坡中40～60 cm土層、全坡面0～10 cm土層應依次作為沉陷黃土坡面土壤有機質、銨態氮、速效磷、速效鉀修復的重點部位。

4 結論

(1)陜北煤礦區沉陷黃土坡面土壤有機質、銨態氮、速效磷、速效鉀含量在不同坡面部位上存在明顯差異，且因土層垂直深度不同而表現出不同的空間變化規律。

(2)陜北煤礦區沉陷黃土坡面土壤有機質、銨態氮、速效鉀含量在不同土層垂直深度上存在明顯差異，且因部位不同而表現出不同的空間變化規律；土壤速效磷含量在各個部位上不同土層垂直深度的差異不明顯。

(3)陜北煤礦區開采沉陷會造成黃土坡面土壤有機質、銨態氮、速效磷等養分全坡面性的損失以及速效鉀的局部性損失和局部性富集。其中坡頂20～40 cm土壤有機質含量損失最為嚴重；坡腳0～10 cm土壤銨態氮含量損失最為嚴重；坡中40～60 cm土壤速效磷含量損失最為嚴重；全坡面0～10 cm土壤速效鉀損失最為嚴重，但開采沉陷也使得全坡面20～40 cm土壤的速效鉀含量出現明顯富集現象。

(4)陜北煤礦區沉陷黃土坡面土壤養分修復應采取人工與自然相結合方式，坡頂、坡中應以人工修復為主，坡腳應以自然恢復為主，并分析了沉陷黃土坡面土壤有機質修復重點部位在坡頂20～40 cm土層，銨態氮修復重點部位在坡腳0～10 cm土層、速效磷修復重點部位在坡中40～60 cm土層、速效鉀修復重點部位在全坡面0～10 cm土層。