基于同位素的神東礦區植物吸水深度及生態水位研究

赫云蘭，楊英明，王璽凱，李可新，房魯鷺，劉卓鳴

(1.中國礦業大學(北京) 煤炭精細勘探與智能開發全國重點實驗室，北京 100083；2.國家能源集團 煤炭開采水資源保護與利用國家重點實驗室，北京 100011；3.中國礦業大學(北京) 地球科學與測繪工程學院，北京 100083)

神東礦區開采煤層埋深小、采高大，高強度的煤炭開采使覆巖的破壞和變形直接傳導到地表，對近地表的水資源、土地資源和地表生態環境產生了不利影響[1-3]。其中地表地形的變化、采煤造成的地裂縫發育和隔水層變化會破壞潛水含水層的結構，使泉水消失、河流徑流量減小、地下水埋深增大，甚至地下水環境發生變化，使依賴地下水生存和受地下水補給影響的地表植被發生退化甚至死亡[4-8]。礦區地下水大規模的下降將引發系統性的生態風險，在干旱半干旱的神東礦區，水資源成為制約植被生長的限制性因素，由于該區域降水季節變化大，持續穩定的地下水源成為許多植物的重要水分來源，因此研究該區域生態安全地下水位具有重要的意義[9，10]。

地下水生態水位是指既能滿足地表陸生植物生長發育對地下水的需求，同時又避免造成土壤鹽漬化和土地荒漠化的地下水動態區間[11]。就地處干旱半干旱區的神東礦區而言，地表多為稀疏荒漠植被，地表植被、土壤水和地下水位相互制約達到區域水資源變化的均衡，增強生態系統的穩定性[12，13]，但采煤造成的地下水位下降可能會對該區域生態系統的穩定性造成威脅。傳統的對于地下水生態水位的研究方法多集中于數值模擬方法[14]、生態統計方法[15]和遙感方法[16]。大多數研究傾向于揭示植物生長狀態與地下水水位之間的關系，但從植物的角度出發，植物根系對于多孔介質中水的吸收能力和水分利用策略也是地下水生態水文研究的一個重要組成部分。近年來興起的水穩定同位素方法成為劃分植物水分來源的有力工具[17]。土壤水穩定同位素廣泛應用于各個學科，如水文學、生態學、土壤科學和生物地球化學等，但水穩定同位素作為一種示蹤工具來描述水在不同的生態、水文和地質區劃中的起源、流動路徑、運輸過程和停留時間的全部潛力尚未被開發。土壤水分是連接不同水體的紐帶，利用水穩定同位素示蹤技術研究干旱半干旱礦區土壤水分蒸發、運移對了解區域水文過程、促進干旱區環境修復、合理生態用水具有重要意義[18，19]。

對于干旱半干旱煤礦區而言，煤炭開采會誘發地裂縫、地下水位下降和地面沉降等環境問題，對地表生態產生一定的擾動。厘清植物水分利用的深度和策略是煤礦區生態修復與治理的一個關鍵因素。Zhao等通過分析黑河流域降水、土壤水和地下水中穩定氧同位素組成，識別和量化了植物的水分吸收，分析了不同氣候條件下植物根吸水所需主要的水源[20]。本文以該區域分布最為廣泛的沙蒿為研究對象[21]，通過人工投放同位素示蹤劑(2H)的方法，根據標記沙蒿植株的莖葉同位素特征，分析該區域井工礦未采動區植物吸水深度[22]和用水策略[23]，從植物學角度探討神東礦區的生態水位。為該研究區生態修復與治理提供了一定的科學依據。

1 研究區概況

神東礦區地理位置位于榆林市神木縣以北，府谷市以西，伊克昭盟的伊金霍洛旗以及鄂爾多斯城的東南，其經緯度坐標為109°51′—110°46′E、38°52′—39°41′N之間，北有毛烏素沙地，南有黃土高原，不僅是兩大高原的過渡地帶，還是水力侵蝕區和風力侵蝕區的過渡帶。

該礦區海拔1200 m左右，地表形態以風沙丘陵和黃土丘陵為主，地表土質類型主要為風沙土。地表植被類型喬木以小葉楊為主，灌木以沙蒿、沙柳、檸條為主，草本植物多為一年生，如硬質早熟禾、沙鞭、狗尾草等。氣候屬于北溫帶半干旱大陸氣候，常年干旱少雨，多年平均降雨量436.7 mm，年蒸發量是降雨量的5～6倍。一年當中降水主要集中在6—8月，最低氣溫出現在12月和1月期間，如圖1所示。研究區既是大型現代化煤炭基地，同時也是典型的生態脆弱區。

圖1 研究區降水與溫度情況Fig.1 Precipitation and temperature in the research area

2 研究方法

本次研究選擇神東礦區井工礦未采動區為研究對象，地表土質類型主要為風沙土。土壤厚度不均，總體厚度在5～8 m。地表植被覆蓋率不高，主要優勢植物為沙蒿。2020年10月，采用人工投放一定豐度同位素水的方法，結合現場原位示蹤試驗探討神東礦區井工礦未采動區植物吸水深度與生態水位。

示蹤試驗開始之前，采集試驗點兩個500 cm深度的土壤剖面分析土壤含水率和氫氧穩定同位素特征。具體實驗方法如圖2所示，首先在現場鉆取不同深度的鉆孔，鉆孔深度選擇30，50，70，100，150，200，300，500 cm，將2H同位素示蹤劑投放到鉆孔底部。研究區植物2H天然同位素本底值在0.015%左右，為了使示蹤效果明顯，人工配比的同位素水的豐度設定15%，是天然值的1000倍左右。投放同位素水后，連續7天采集不同深度鉆孔周邊沙蒿莖葉，冷藏保存并送至實驗室提取、測試同位素值，判斷植物吸水最大深度，采用MixSIAR同位素混合模型計算分析土壤水穩定同位素，劃分植物吸收不同深度段土壤水分的比例。

圖2 現場原位示蹤實驗Fig.2 On site in-situ tracing experimental diagram

3 研究結果與討論

3.1 土壤剖面含水率

不同地區、不同氣候、不同土壤結構其在垂向上土壤含水率都不同，土壤含水率對于地表生態環境與植物的生長至關重要。研究區土壤剖面上含水率如圖3所示，雨季過后，地表到60 cm深度含水率上升并保持較高水平，60～160 cm深度段土壤含水率快速下降到2%左右，160～340 cm深度段土壤含水率緩慢上升達到3%，340～500 cm深度段土壤含水率快速波動上升到5%左右。

圖3 土壤剖面含水率Fig.3 Soil moisture profile

3.2 植物吸水極限深度

利用植物莖葉2H同位素特征，判斷植物根系吸水的最大深度，植物吸收到同位素水，莖葉的同位素值會較本底值大。本次現場原位示蹤實驗選取了不同深度投放同位素水，并在地表連續7天采集鉆孔周邊沙蒿的莖葉進行測試，其中為確保沙蒿根系相當，選擇高度相近的沙蒿采集樣品，如圖4所示。

圖4 研究區植物莖葉δ2H同位素特征Fig.4 δ2H isotope characteristics in plant stems and leaves in the research area

實驗結果顯示，同位素水投放不同深度情況下沙蒿連續7天采集的樣品δ2H同位素豐度不同，這與同位素水投放深度和植物根系吸水利用率有關。

從圖4中可以看出，近地表30 cm到200 cm植物莖葉的同位素值基本均高于本底值。而投放在300 cm深度時，沙蒿δ2H同位素只有第2天高于本底值。投放500 cm深度下，植物已吸收不到高豐度的同位素水。在300 cm深度條件下，第1天沙蒿δ2H同位素低于本底值，第2天高于于本底值，后續采樣均低于本底值。由于人工投放同位素水時，鉆孔底部會形成局部過飽和的狀態，同位素水會同時向下和向上產生一定的運移。在300 cm投放同位素水時，從土壤剖面上同位素富集特征看，分別向上和向下運移了50 cm左右，所以判斷植物沙蒿吸水的極限深度為2.5～3 m，如圖5所示。

圖5 300 cm深度土壤水δ2H和δ18O同位素特征Fig.5 δ2H and δ18O isotope characteristics in the soil water at a depth of 300 cm

3.3 土壤剖面水分利用

植物根系對土壤水分利用率在垂向上具有差異，這與植物類型、根系對不同深度土壤水分吸收效率、季節有關。利用天然δ2H和δ18O同位素示蹤可以對植物吸水來源和土壤剖面利用深度進行判斷，太湖流域河岸植物銀杏主要利用0～60 cm深度的土壤水和降水，而桑樹則主要利用120～150 cm的土壤水和地下水[24]。MixSIAR同位素混合模型可以量化不同深度土壤水分對植物吸收水分的貢獻，這對于分析植物主要利用那個剖面深度的水分提供了有效的方法和手段。

本研究通過人工投放高豐度的δ2H同位素示蹤劑，提高了示蹤效果。植物樣品采集于2020年10月，利用MixSIAR同位素混合模型分析了研究區沙蒿對不同深度土壤水的利用率。將240 cm深度的土壤剖面劃分了12個層段，經模型計算得出12個層段土壤水對沙蒿吸水的貢獻率相差不大，如圖6所示，其中，60～80 cm，80～100 cm和100～120 cm深度段相對水分貢獻率較低，其他層段貢獻率都接近9%。由圖6可以看出：60～120 cm深度段土壤含水率快速下降，此時正值雨季過后，推測雨季的降水濕潤鋒剛入滲到該深度段，之前該深度段為水分比較匱乏的區域，沙蒿在此段的水分吸收比例較低；沙蒿水分吸收雖然受到突然含水率高低的影響，但各深度段的水分吸收比例差別不大，顯示出其廣泛而較均勻的水分吸收策略。

圖6 不同深度土壤水分利用情況Fig.6 Soil moisture utilized by plants at different depths

4 結 語

氫氧同位素本身作為水的一部分，在研究水循環方面是理想的示蹤劑。采用該方法在神東礦區未采動區域開展了現場原位示蹤實驗研究，獲得了較好的示蹤效果。根據優勢植物沙蒿和土壤水δ2H同位素特征，研究區沙蒿根系最大吸水深度在2.5～3 m左右。雨季過后研究區域土壤含水率淺處和深處高，中間低。10月份沙蒿對土壤剖面水分利用率受到土壤含水率高低的影響但總體表現相差不多，60～120 cm層段利用率相對較低，其他層段利用率接近9%，顯示出其較均勻的水分吸收策略。本研究揭示了神東礦區沙蒿的根系吸水深度和水分吸收策略，從植物的角度測量出該區域地下水補充植物水分的生態水位下限為2.5～3 m。