基于ERT的黃土高原小流域典型坡地和壩地土壤儲水量反演

段國秀，賈小旭,3*，白 曉，劉成功,3，魏孝榮,3

基于ERT的黃土高原小流域典型坡地和壩地土壤儲水量反演

段國秀1，賈小旭1,3*，白 曉2，劉成功1,3，魏孝榮2,3

（1.中國科學院 地理科學與資源研究所/生態系統網絡觀測與模擬重點實驗室，北京 100101；2.西北農林科技大學，陜西 楊凌 712100；3.中國科學院大學 資源與環境學院，北京 100190）

【目的】反演黃土高原小流域典型坡地和壩地的土壤儲水量（），為小流域水資源管理與生態建設提供科學依據。【方法】在陜北六道溝小流域選擇典型坡地和壩地，利用電阻率層析成像法（ERT）測定其土壤電阻率（），結合實測土壤體積含水率（v），構建與v的冪函數模型，量化0～12.5 m土層剖面的，并對比分析坡地與壩地的空間分布差異。【結果】對于黃土高原坡地和壩地2種典型地貌類型，與v之間均具有顯著的冪函數關系（坡地：2=0.66，<0.01；壩地：2=0.64，<0.01）；坡地分布受坡位影響顯著，沿坡頂至坡底逐漸降低，壩地在垂直方向上呈“高-低-高”的分布格局，具有明顯的飽和層和非飽和層，而坡地剖面均為非飽和狀態；面積為1.8×104m2的壩地0～6.5 m非飽和層的可儲水量為1.49×104m3，6.5～12.5 m飽和層可儲水量為5.10×104m3；而相近面積的坡地0～12.5 m非飽和層可儲水量為2.59×104m3。【結論】ERT可用于黃土區深層土壤水資源時空分布研究，六道溝小流域壩地土壤儲水量豐富，是坡地土壤儲水量的2.5倍，對小流域的水循環和生態系統服務具有重要作用。

黃土高原；電阻率層析成像法；土壤電阻率；坡地；壩地；土壤含水率

0 引 言

【研究意義】黃土高原氣候干旱，降水分布不均，水土流失嚴重，是我國典型的生態脆弱區和氣候敏感區[1]，其生態建設是鞏固北方生態屏障、促進西部經濟社會可持續發展的重要保障。植被恢復和淤地壩建設分別是該地區水土流失治理的重要生物措施和工程措施[2]。隨著退耕還林（草）工程的大規模實施，小流域坡地土地利用發生了顯著改變，植被覆蓋度顯著增加，加劇了深層土壤水分的消耗[3]。目前，土壤水分供給不足已成為制約該區生態建設和植被恢復的主要因素[4]。與坡地不同，淤地壩內淤積而成的壩地水肥條件優越，是該區建設基本農田的重要地理單元。當前，黃土高原地區共有淤地壩約12萬座，淤成壩地約31萬hm2[5-6]，在提高農業生產力、防治水土流失和改善生態環境等方面發揮著重要作用。量化小流域典型坡地和壩地的土壤儲水量（）及其空間分布格局，對小流域水資源優化管理、生態恢復和可持續的農業生產具有重要科學意義[1,7]。

【研究進展】當前，關于土壤水分時空分布的研究大多采用烘干法、中子儀法、原位傳感器等方法測定土壤含水率（）[8-9]，但這些方法監測范圍較小，存在空間上的不連續性且費時、耗力等缺點，限制了對深層土壤水分時空特征的探索。宇宙射線快中子法和遙感反演方法雖然能獲取更大尺度的時空分布，但獲取深度和空間分辨率有限[10]。發展大尺度和高空間分辨率的監測技術是解決這一問題的關鍵。探地雷達（Ground Penetrating Radar，GPR）、核磁共振（Nuclear Magnetic Resonance Imaging，NMR）、電阻率層析成像（Electrical Resistivity Tomography，ERT）等地球物理技術是將土壤剖面特征（如土壤結構、土壤水分）與地球物理信號相關聯，能實現多尺度高分辨率的土壤特征監測，極大地促進了水文地理的發展[7,11]。研究表明，ERT可實現多尺度土壤水文狀況的原位連續監測，在土壤水和地下水的空間格局研究中具有較高的應用潛力和價值[12-14]。然而，在黃土高原地區利用ERT測量土壤剖面電阻率（），進而反演的研究較少[15-16]，缺乏小流域坡地和淤地壩不同水分條件的對比研究。岳寧等[15]利用ERT和水分傳感器測定2 m深度土壤剖面的和，建立了隴中半干旱區農田與之間的定量關系。Sun等[17]在黃土高原地區分別建立了多因素影響下的和的關系模型，但模型受眾多因素限制，無法大范圍推廣應用。Celano等[18]利用實測數據建立了與的指數模型，較好地評估了果園土壤水分的季節變化。張世斌等[19]利用室內模擬試驗建立了黃土與飽和度之間的關系式，重塑了黃土水分遷移特性。袁鴻猷等[16]利用ERT測定了黃土高原壩地，基于和的關系模型估算了壩地淺層地下水儲量。以往研究對利用ERT技術反演土壤水分或地下水儲量進行了探索，但對于小流域坡地和壩地的對比研究尚屬缺乏，對不同地類下的土壤儲水特征認識不足，且和的反演模型也存在較大差異，已有研究結果難以被廣泛應用。因此，本研究利用ERT開展黃土區小流域典型坡地和淤地壩的及時空分布的反演研究，對于小流域坡地植被建設和壩地農業生產管理具有重要的科學意義。

【切入點】以往研究大多關注黃土高原坡面土壤水分及區域尺度土壤水分的時空變異性，由于監測技術的限制，相關土壤水分時空分布及儲水量估算的研究大多集中在5 m土層內[8-9, 20-24]，僅有少數研究關注深層[24-27]。深層土壤水分是黃土區人工植被生長和應對極端干旱的重要水源，一旦耗盡將難以恢復[1,3]，其信息的缺失限制了黃土高原深層土壤水的可持續管理。【擬解決的關鍵問題】鑒于此，本研究以陜北神木六道溝小流域土壤質地較為均勻的典型坡地和壩地為研究對象，利用ERT測定的二維空間分布，基于實測剖面土壤體積含水率（v），建立和v的定量關系模型，探究ERT在反演坡地和壩地方面的適用性，并估算深層，分析坡地和壩地深層土壤儲水特征及空間分布格局，以期為黃土高原小流域不同地貌單元土壤水定量評估與管理提供科學依據和技術支撐。

1 材料與方法

1.1 研究區概況

研究區位于陜西省神木市以西14 km處的六道溝小流域（圖1（a）），海拔1 080～1 275 m。該流域地處黃土高原與毛烏素沙地的過渡地帶，是黃土高原水蝕、風蝕交錯帶強烈的侵蝕區。年平均氣溫為8.4 ℃，年平均降水量為427 mm，屬于典型的半干旱大陸性季風氣候區。由于水蝕和風蝕的交替危害及人為濫砍濫伐和過度放牧，天然喬、灌叢植被幾乎破壞殆盡。經人工植被恢復后，目前小流域主要有草地（44%）、灌木林地（26%）和農田（16%）3種土地利用類型[27]。草地包括以長芒草為主的天然草地和以苜蓿為主的人工草地，灌木地以檸條為主，農用地以1 a生農作物為主[28]。

注 B為壩地，P為坡地；數字代表每個樣地中監測電阻率時布設的測線編號。

本研究所選坡地為小流域內遠離侵蝕溝的一個典型坡面（圖1（b）），坡長350 m，海拔落差42 m。坡上植被以人工植被為主，其中坡地上部110 m內是種植32 a的苜蓿地，坡中部是種植16 a的杏樹林，坡底部為種植16 a的杏樹林和零星分布的刺槐和鐵桿蒿[29]。壩地位于所選坡地底部的西北面，總長度為550 m[10]，面積約為1.8×104m2。自壩頭至壩尾分布有檸條灌木地、撂荒灌草地、苜蓿草地和玉米為主的農耕地，兩側零星分布有少許榆樹和旱柳。沉積泥沙主要為粉粒泥沙，質地較為均一[30]。

1.2 數據獲取

于2019年7月7日，利用美國AGI（Advanced Geosciences, Inc.）公司研制的SuperSting-R8/IP高密度電導儀，采用對電阻率垂向變異較為敏感的Wenner電極陣列，對坡地和壩地剖面土壤電阻率（m）進行了測定[31-32]。對于坡地，沿坡面布設3條間隔16.5 m的平行測線（P1—P3）。對于壩地，從壩頭至壩尾布設1條測線B4，同時在垂直于B4測線方向另外布設3條測線（B1—B3），測線具體信息和位置分別見表1和圖1。

表1 坡地和壩地各樣線基本情況

測線P1—P3布設在距坡面原有的3列中子管0.5 m距離處，每列共有28根。測線B4沿原有的22根中子管布設，所有中子管測深均為5.2 m。在測定土壤剖面m的同時，利用中子儀（型號CNC503DR）測定相應位置0～5 m不同土層深度的v，其中0～1 m土層間隔10 cm，1～5 m土層間隔20 cm。此外，在壩地中部安裝自記式水位計（型號ZYHC-300B）以監測地下水位狀況。利用環刀采集坡地和壩地剖面原狀土壤樣品，在室內測定飽和土壤含水率（s）。

1.3 數據處理與分析

土壤真實由田間測定的m反演得到。數據反演前，需對m進行質量控制，即刪除滿足下列條件的m值：①m高于5 000 Ω m的數值；②重復測量誤差大于3%的數值[17,33-34]；③負值。使用EarthImager 2D軟件，采用平滑模型對m進行反演，利用m值和模型反演的值計算得到均方根誤差（），對數據擬合程度進行評價，并通過EarthImager 2D軟件中的擬合誤差直方圖，在考慮誤差分布的基礎上適當調整擬合閾值，去除誤差高于100%的數據點再次反演[17]。

一般來說，ERT樣線兩端的靈敏度和分辨率較低[17]。因此，選擇位于樣線中部的中子管實測v及其對應位置的建立數據集，隨機選取其中80%的數據進行建模，其余20%用于模型驗證。利用決定系數（2）和對模型精度進行評價。

壩地土壤剖面包含非飽和層和飽和層。通過和v之間的關系模型計算各測點的土壤含水率vi，然后獲得非飽和層的平均土壤含水率`vi，利用式（1）計算非飽和層的。飽和層的通過s計算獲得，計算方法見式（2）[20]。飽和層的上下邊界利用的空間分布并結合地下水位觀測數據進行判斷。

式中：1和2分別為非飽和層和飽和層的（m3）；s為土體積（m3）；`θ和s分別為非飽和層土壤含水率（%）和飽和層土壤含水率（%）。

為了與壩地進行對比分析，坡地則是在整個坡面選擇與壩地面積相近、深度一致的土體進行反演和計算。根據數據點相對位置在坡面建立矩形網格，計算每個矩形網格的v，共有1 710個矩形單元網格。中心點坐標定義為（，，），為3 m，代表相鄰兩電極間距的50%（電極間距6 m），代表與當前電極所在樣線相鄰平行樣線之間水平距離的50%，代表矩形網格中心點垂向深度[15]。現場作業的命令文件在編寫時根據測距、電極間距和最大測深將土壤剖面自動劃分為5層，各層深度分別為0.38、3.10、6.01、9.11 m和12.43 m。

2 結果與分析

2.1 坡地和壩地土壤電阻率空間分布特征

坡地剖面的空間分布如圖2所示。在沿坡面方向上，沿坡面由坡頂至坡底呈逐漸下降趨勢，具體為坡頂（0～100 m）>坡中（100～250 m）>坡底（250～330 m）。這種分布特征與坡面降水再分配有關，部分降水以徑流的形式沿坡面向坡底匯聚，導致坡底剖面v較高，較低。在垂直方向，地表至20 m深度內的土層相比更深土層的值高，特別是在坡頂位置，這可能是由于坡頂降水補給量較小且深根系人工植被過度消耗土壤水分，導致土壤干燥化[3,20]。

注 Iteration表示迭代次數，RMSE表示ERT實測視電阻率數據與軟件反演時重建模型計算數據之間的均方根誤差；L2為視電阻率值加權誤差的平方和；EP（Electrode Spacing）為電極間距。

壩地剖面在垂直方向上呈“高-低-高”的層狀分布格局（圖3）。淺層高阻區厚度自壩尾向壩頭方向逐漸增加，這主要是由于植被類型的不同。淺層高阻區以下至15 m深度范圍為連續帶狀低阻區，值在16～50 Ω m范圍內分布，這與地下水的觀測結果基本一致，是壩地飽和含水層的分布區域，與袁鴻猷等[16]研究結果相近。根據壩地剖面分布規律和值大小，選取該低阻區同一測點下值變化差異最大的2個點，分別提取這2個點各自對應的深度，即為該測點處飽和含水層上、下邊界的深度，取各測點的平均值作為壩地飽和含水層的上、下邊界深度，分別為6.5、12.5 m。飽和含水層上邊界深度與壩地地下水觀測的水位結果基本一致，觀測當月地下水平均水位深度也為6.5 m。另外，飽和含水層以下的高阻區可能是由基巖所致[10]。

注 Iteration=4，RMSE=2.83%，L2=0.89，EP=6 m。

與壩地縱向（即壩尾—壩頭方向）不同，壩地橫向（即與壩尾—壩頭的垂直方向）剖面空間分布較為復雜，橫向剖面B1、B2和B3分別位于壩頭、壩中和壩尾，各測線空間變異較大，特別是B1和B2這2條測線（圖4）。土壤淺層高電阻區由壩尾至壩頭由1 m（B3）增至7 m（B1），這可能是由于水流挾沙能力和重力等多因素影響，使粗顆粒在壩頭先沉積，壩頭土壤砂粒量較高，使得剖面較高[16,30]。同時，壩頭至壩尾地勢逐漸降低，部分降水在壩尾匯集，淺層地下水埋深逐漸變淺，高阻區厚度隨之減小，如圖4（c）中位于壩尾的B3剖面地勢最低，具有較為連續的飽和帶分布特征。此外，3條樣線地表植被類型有顯著差異，根系耗水量不同也會導致剖面的差異。如圖4（a）位于壩頭的B1中部是以高耗水檸條為主的灌木林地，而兩側零星分布有榆樹、旱柳等喬木和少許撂荒草本。以檸條為主的人工灌木林地根系耗水深度可達21 m[27]。B2測線為以苜蓿、鐵桿蒿為主的撂荒草地，樣線剖面分布的高阻區也是由于苜蓿根系過度耗水導致v降低所致。此外，B2測線四周分布有小的洼地，可能存在測量時相應的電極與地面接觸不良，且為人工苜蓿地和撂荒草地復合區域，植被類型和微地形均較為復雜，其和迭代次數比其余剖面均較大[17]（圖4（b））。而B3測線則是以玉米為主的農耕地，作物根系主要分布在1 m土層內，對深層v影響很小，故B3測線剖面分布較為均勻[27]。

圖4 壩地橫向（與壩尾-壩頭垂直方向）剖面土壤電阻率空間分布

2.2 土壤含水率與電阻率的關系模型

坡地和壩地剖面變異系數（）分別為40%和26%，v的分別為30%和56%，均具有較高的變異性。壩地的中值和平均值（46.26、46.20 Ω m）較坡地（79.1、74.61 Ω m）低，而v的中值和平均值（19.40%、18.53%）比坡地（14.45%、11.43%）高。袁水龍等[23]研究發現，壩地土壤水分呈層狀分布特征，0～0.4 m和1.4 m以下壩地的v明顯高于坡地，而0.4～1.4 m坡地v高于壩地，這與壩地植被根系耗水有關。

利用在壩地和坡地獲得的與v數據集，分別建立適于坡地和壩地的與v定量關系模型，并驗證其準確性，結果如圖5和圖6所示。坡地和壩地v實測值與模擬值之間均具有顯著的相關關系（<0.01），2分別為0.61和0.75，分別為3.66和3.15，表明建立的-v關系模型可以較準確地反演坡地和壩地土壤水分狀況。Celano等[18]在2種管理措施（耕作和覆蓋）的橄欖園地也建立了v與顯著相關的指數函數模型，用于評價2種果園秋春期v的變化，通過高分辨率的土壤結構二維分布與水分運移過程的圖像，研究剖面v分布的動態差異，評估深層土壤（>1.0 m）的v，發現覆蓋具有顯著的土壤儲水性能。岳寧等[15]在隴中半干旱區玉米田對降雨前后土壤二維剖面進行監測，建立了v與的線性關系（2=0.65，=96），探究了不同條件下土壤含水率變化的可能原因。

圖5 坡地土壤電阻率與土壤含水率的回歸模型與模型評價

圖6 壩地土壤電阻率與土壤含水率的回歸模型與模型評價

2.3 坡地和壩地土壤儲水量估算

小流域壩地的估算分為土壤水飽和帶和非飽和帶2部分。其中，土壤水飽和帶上、下邊界深度平均值分別為6.5、12.5 m。依據式（2）估算獲得面積為1.8×104m2、厚度為6.0 m的壩地土壤飽和含水層的為5.10×104m3。利用壩地-v定量關系模型獲得非飽和層的v，其中0～1.44、1.44～3.03 m和3.03～6.50 m土層的v分別為10.58%、13.76%和13.89%。在0～6.5 m非飽和層，隨土層深度的增加，基于反演獲取的θ逐漸增加。依據式（1）估算獲得面積為1.8×104m2、厚度為6.5 m的壩地非飽和層為1.49×104m3。與壩地不同，所選坡地0～12.5 m均為非飽和狀態。因此，利用坡地-v定量關系模型計算獲得坡地剖面v，依據式（1）估算獲得面積1.75×104m2、厚度為12.5 m的坡地剖面為2.59×104m3。壩地單位面積0～12.5 m剖面的達到3.66 m3，而相同深度的坡地單位面積剖面的為1.48 m3，壩地約為坡地的2.5倍，表明陜北六道溝小流域壩地儲水量豐富，是小流域水資源重要的存儲場所。

3 討 論

受土壤溫度、質地、溶質濃度等多種因素影響，而在特定環境中，往往受個別因素主導[35]。例如坡地P1測線因其v受坡位和植被影響，導致P1測線受坡位和植被影響顯著：坡頂植被以人工苜蓿植被為主，植被耗水量大，加之受風速、太陽輻射等強烈影響，蒸散量大，降水補給到深層土壤中的水分更少，導致剖面土壤較高[20]。與P1測線對應的0～4 m剖面v也具有相似的空間分布特征，即坡頂高阻區對應的0～80 m距離范圍內低v區域比坡中和坡底面積更大，且深度更深（圖7（a））。可見，雖然人工高耗水植被分布于整個坡面，但高阻區主要集中在坡頂位置，坡地分布特征受植被和坡位的共同影響[20]。

然而，由于壩地大部分中子管進水導致v測定深度較淺，僅為0～80 cm，圖3壩地B4測線對應的剖面v難以與剖面進行對比分析（圖7（b））。但圖3水平方向自壩尾至壩頭，淺層值逐漸增大，植被類型依次為農耕地、苜蓿草地、撂荒灌草地和檸條灌木地，表明其分布特征可能與植被耗水量增加有關，這與劉成功等[10]的研究結果一致。垂直方向壩地呈明顯的“高-低-高”分層分布格局，袁水龍等[23]長期監測、分析了黃土高原一壩地土壤水分時空分布特征，結果也表明，壩地土壤水分具有明顯的分層現象，表層變化劇烈，隨著深度的增加變化程度減弱，且壩尾各層（0.6～1.4 m）均明顯高于壩中和壩頭。此外，壩地中層呈現連續帶狀低阻區，表明壩地中部存在土壤水飽和層。張翔等[36]利用同位素示蹤技術分析陜北綏德韭園溝流域一壩地土壤水的來源及分布，發現其土壤水垂直方向賦存規律從“較為平均分布的地表徑流為主”逐漸向“地下水為主、地表徑流為輔”轉變，這也表明部分壩地中層存在土壤水飽和層。此外，鄒俊亮等[37]在本研究的壩地采集壩淤土時，在5.3 m深度已見水，可見多年來該壩地中部保持有較多的土壤水，這些研究結果也表明壩地v存在分層分布規律，具有一定的蓄水潛能。目前現有觀測方法難以監測確定壩地地下水位下邊界[38-40]，而ERT可用于估計土壤水飽和帶上下邊界深度進而估算地下水資源儲量。

由于輸入參數的數量和種類在理論和實踐上都影響模型的精度，而輸入參數具有高度的空間異質性和尺度依賴性[10]，且受多種因素的共同影響，很難建立高精度、多因子模型。考慮到特定環境中，往往受個別因素主導[31]，以及本研究結果在田間尺度上的應用性，本文選取土壤質地較為均一的地區開展監測[10, 41]，降低了其余因素對的影響，故認為v是影響該區的主要因子。本研究建立的特定環境下-v定量關系模型也具有較高的模擬精度，可用于黃土區坡地和壩地深層的估算。然而，本研究尚未探討與土壤不同理化性質之間的作用機制，在未來研究中應考慮氣象、地形、環境等因子對的影響，室內外試驗相結合，并加強土壤影響因子與之間作用機理的研究，建立多因子、高精度、普適性的物理耦合模型。

壩地在防止水土流失、糧食增產、改善區域生態環境等方面發揮著顯著的經濟、社會和生態效益[4]。研究表明，淤地壩建設在小流域系統中表現為土壤的碳和氮匯[42]，也有研究發現壩地是增加小流域地下水資源量的重要來源[38]，這部分水資源是黃土高原小流域不可忽視的寶貴水源，對生態系統功能的持續發揮具有重要作用。然而，賦存于不同壩地的水資源形式和儲量較為復雜，受諸如植被、氣候、土壤、排水等多種因素的影響，對于壩地水資源賦存特征和具體數量的研究還處在探索階段。未來應加強壩地水資源形成過程、賦存特征、轉化機制、儲量及經濟、生態效益等方面的研究，更深入了解小流域坡地和壩地系統水循環過程，以評價小流域水量平衡。本研究基于-v定量關系模型估算的壩地也存在著一定的不確定性。例如，由于水流挾沙能力和重力等多因素的影響，沉積顆粒的分選過程使壩地在淤積形成過程中粗顆粒在壩頭先沉積，垂直方向其土壤質地也存在分層現象，而本研究并未考慮土壤質地在壩地水平和垂直方向的異質性。盡管如此，土壤電學特性較其他土壤性質更易測得，ERT技術仍是高效獲取深層的有效方法，在多尺度土壤水文過程和土壤電學性質的無損監測中具有較好的潛力。

4 結 論

1）小流域典型坡地剖面土壤沿坡頂至坡底逐漸降低，主要受坡位和植被的影響；而壩地土壤空間分布較為復雜，具有明顯的水平和垂直分異特征，垂直方向呈“高-低-高”的分布格局，水平方向淺層土壤受植被類型影響顯著。壩地具有明顯的土壤水非飽和層和飽和層，而坡地剖面均為非飽和層。

2）小流域典型坡地和壩地土壤與v之間均具有顯著的冪函數關系，利用ERT可較為準確地反演深層v并獲得飽和含水層厚度，可用于黃土區不同地貌單元土壤水非飽和層和飽和層儲水量及空間分布研究。

3）神木六道溝小流域面積1.8×104m2的壩地和坡地0～12.5 m的分別為6.59×104和2.59×104m3，表明壩地比坡地儲存了更多的水資源，水儲量是坡地的2.5倍，對小流域水循環與生態系統服務功能具有重要作用。

[1] JIA Xiaoxu, SHAO Mingan, WEI Xiaorong, et al. Policy development for sustainable soil water use on China’s Loess Plateau[J]. Science Bulletin (Beijing), 2020, 65(24): 2 053-2 056.

[2] 劉曉燕, 高云飛, 馬三保, 等. 黃土高原淤地壩的減沙作用及其時效性[J]. 水利學報, 2018, 49(2): 145-155.

LIU Xiaoyan, GAO Yunfei, MA Sanbao, et al. Sediment reduction of warping dams and its timeliness in the Loess Plateau[J]. Journal of Hydraulic Engineering, 2018, 49(2): 145-155.

[3] SHAO Mingan, WANG Yunqiang, XIA Yongqiu, et al. Soil drought and water carrying capacity for cegetation in the critical zone of the Loess Plateau: A review[J]. Vadose Zone Journal, 2018, 17(1): 1-8.

[4] 賈小旭, 邵明安, 張晨成, 等. 黃土高原南北樣帶不同土層土壤水分變異與模擬[J]. 水科學進展, 2016, 27(4): 520-528.

JIA Xiaoxu, SHAO Mingan, ZHANG Chencheng, et al. Variation and simulation of soil water content within different soil depths along the south-north transect of the Loess Plateau[J]. Advances in Water Science, 2016, 27(4): 520-528.

[5] 艾開開. 黃土高原淤地壩發展變遷研究[D]. 楊凌: 西北農林科技大學, 2019.

AI Kaikai, Study on development and change of warping dams on Loess Plateau[D]. Yangling: Northwest A & F University, 2019.

[6] 王云強, 張興昌, 韓鳳朋. 黃土高原淤地壩土壤性質剖面變化規律及其功能探討[J]. 環境科學, 2008, 29(4): 1 020-1 026.

WANG Yunqiang, ZHANG Xingchang, HAN Fengpeng. Profile variability of soil properties in Check Dam on the Loess Plateau and Its Functions[J]. Environmental Science, 2008, 29(4): 1 020-1 026.

[7] LIU Wenzhao, ZHANG Xunchang, DANG Tinghui, et al. Soil water dynamics and deep soil recharge in a record wet year in the southern Loess Plateau of China[J]. Agricultural Water Management, 2010, 97(8): 1 133-1 138.

[8] 孟婷婷, 王歡元, 劉金寶, 等. 黃土丘陵區退耕還林草對深層土壤水分動態的影響[J]. 灌溉排水學報, 2021, 40(12): 78-84.

MENG Tingting, WANG Huanyuan, LIU Jinbao, et al. Changes in subsoil water content after withdrawing crop production in the hilly loess in Northwestern China[J]. Journal of Irrigation and Drainage, 2021, 40(12): 78-84.

[9] 孟秦倩, 王健, 吳發啟. 延安丘陵溝壑區坡面果園土壤水庫特征[J]. 灌溉排水學報, 2008, 27(1): 93-95, 99.

MENG Qinqian, WANG Jian, WU Faqi. Characteristics of slope-orchard soil reservoir in Yan’an loess hilly region[J]. Journal of Irrigation and Drainage, 2008, 27(1): 93-95, 99.

[10] 劉成功, 賈小旭, 趙春雷, 等. 黃土區坡地和壩地土壤電阻率分布特征及應用[J]. 土壤學報, 2022, 59(1): 139-147.

LIU Chenggong, JIA Xiaoxu, ZHAO Chunlei, et al. Distribution and application of soil electrical resistivity in slope land and dam land in the Loess Region[J]. Acta Pedologica Sinica, 2022, 59(1): 139-147.

[11] 段國秀, 賈小旭, 白曉, 等. 電阻率層析成像法在土壤水文學中的應用：基于CiteSpace的計量分析[J]. 土壤通報, 2021, 52(6): 1 447-1 459.

DUAN Guoxiu, JIA Xiaoxu, BAI Xiao, et al. Application and progress of electrical resistivity tomography in pedohydrology：Citespace-based quantitative analysis[J]. Chinese Journal of Soil Science, 2021, 52(6): 1 447-1 459.

[12] BEFF Laure, GüNTHER Thomas, VANDOORNE Bertrand, et al. Three-dimensional monitoring of soil water content in a maize field using Electrical Resistivity Tomography[J]. Hydrology and Earth System Sciences, 2013, 17(7): 595-609.

[13] BRUNET Pascal, CLéMENT Rémi, BOUVIER Christophe. Monitoring soil water content and deficit using Electrical Resistivity Tomography (ERT): A case study in the Cevennes area, France[J]. Journal of Hydrology (Amsterdam), 2010, 380(1-2): 146-153.

[14] MICHOT Didier, BENDERITTER Yves, DORIGNY Abel, et al. Spatial and temporal monitoring of soil water content with an irrigated corn crop cover using surface electrical resistivity tomography[J]. Water Resources Research, 2003, 39(5): 1 138.

[15] 岳寧, 董軍, 李玲, 等. 基于高密度電阻率成像法的隴中半干旱區土壤含水量監測研究[J]. 中國生態農業學報, 2016, 24(10): 1 417-1 427.

YUE Ning, DONG Jun, LI Ling, et al. A soil water content monitoring in Longzhong semi-arid region by high-density electrical resistivity tomography[J]. Chinese Journal of Eco-Agriculture, 2016, 24(10): 1 417-1 427.

[16] 袁鴻猷, 樊軍, 金沐, 等. 黃土高原淤地壩土壤水分和淺層地下水時空分布特征解析[J]. 灌溉排水學報, 2020, 39(10): 50-56.

YUAN Hongyou, FAN Jun, JIN Mu, et al. Spatiotemporal distribution of soil water and shallow groundwater in check dams in the loess plateau of China[J]. Journal of Irrigation and Drainage, 2020, 39(10): 50-56.

[17] SUN Hui, WANG Yunqiang, ZHAO Yali, et al. Assessing the value of electrical resistivity derived soil water content: Insights from a case study in the Critical Zone of the Chinese Loess Plateau[J]. Journal of Hydrology, 2020, 589: 125 132.

[18] CELANO Giuseppe, PALESE Assunta Maria, CIUCCI A, et al. Evaluation of soil water content in tilled and cover-cropped olive orchards by the geoelectrical technique[J]. Geoderma, 2011, 163(3-4): 163-170.

[19] 張世斌, 朱才輝, 袁繼國. 降雨條件下重塑黃土中水分遷移模型試驗研究[J]. 水利學報, 2019, 50(5): 621-630.

ZHANG Shibin, ZHU Caihui, YUAN Jiguo. Laboratory model tests on moisture migration in remolded loess under rainfall conditions[J]. Journal of Hydraulic Engineering, 2019, 50(5): 621-630.

[20] MEI Xuemei, MA Lan, ZHU Qingke, et al. The variability in soil water storage on the loess hillslopes in China and its estimation[J]. Catena, 2019, 172: 807-818.

[21] 白曉, 賈小旭, 邵明安, 等. 黃土高原北部土地利用變化對長期土壤水分平衡影響模擬[J]. 水科學進展, 2021, 32(1): 109-119.

BAI Xiao, JIA Xiaoxu, SHAO Mingan, et al. Simulating long-term soil water balance in response to land use change in the Northern China’s Loess Plateau[J]. Advances in Water Science, 2021, 32(1): 109-119.

[22] LIU Yu, MIAO Haitao, HUANG Ze, et al. Soil water depletion patterns of artificial forest species and ages on the Loess Plateau (China)[J]. Forest Ecology and Management, 2018(417): 137-143.

[23] 袁水龍, 李占斌, 張揚, 等. 黃土高原小流域壩地水分時空分布特征[J]. 水土保持研究, 2018, 25(3): 29-34.

YUAN Shuilong, LI Zhanbin, ZHANG Yang, et al. Impact of layered deposition on temporal and spatial distribution characterisitic of soil moisture of check dam land[J]. Research of Soil and Water Conservation, 2018, 25(3): 29-34.

[24] 劉中奇, 朱清科, 秦偉, 等. 半干旱黃土區典型林地土壤水分消耗與補給動態研究[J]. 灌溉排水學報, 2010, 29(3):36-40.

LIU Zhongqi, ZHU Qingke, QIN Wei, et al. Soil water consumption and supply of typical forest lands in the Semiarid Loess Region[J]. Journal of Irrigation and Drainage, 2010, 29(3):36-40.

[25] JIA Xiaoxu, SHAO Mingan, WEI Xiaorong, et al. Hillslope scale temporal stability of soil water storage in diverse soil layers[J]. Journal of Hydrology (Amsterdam), 2013, 498: 254-264.

[26] HU Wei, SHAO Mingan, REICHARDT Klaus. Using a new criterion to identify sites for mean soil water storage evaluation[J]. Soil ence Society of America Journal, 2010, 74(3): 762-773.

[27] WANG Yunqiang, HU Wei, ZHU Yuanjun, et al. Vertical distribution and temporal stability of soil water in 21-m profiles under different land uses on the Loess Plateau in China[J]. Journal of Hydrology, 2015, 527: 543-554.

[28] 劉丙霞. 黃土區典型灌草植被土壤水分時空分布及其植被承載力研究[D]. 北京: 中國科學院大學, 2015.

LIU Bingxia. Experimental study of soil water spatial-temporal distribution and soil water carrying capacity for vegetation of typical shrub and grass on the Northern Loess Plateau[D]. Beijing: University of Chinese Academy of Sciences, 2015.

[29] 白一茹, 邵明安. 黃土高原水蝕風蝕交錯帶不同土地利用方式坡面土壤水分特性研究[J]. 干旱地區農業研究, 2009, 27(1): 122-129.

BAI Yiru, SHAO Mingan. Soil water properties in a slope of different land use in the wind-water erosion crisscross region on the Loess Plateau[J]. Agricultural Research in the Arid Areas, 2009, 27(1): 122-129.

[30] 張瑋, 楊明義, 張風寶, 等. 黃丘區小流域壩地沉積泥沙粒徑剖面分布特征[J]. 水土保持研究, 2015, 22(2): 17-21.

ZHANG Wei, YANG Mingyi, ZHANG Fengbao, et al. Profile Distribution of Particle Size of Sediment at a Check Dam in a Small Watershed of the Loess Plateau[J]. Research of Soil and Water Conservation, 2015, 22(2): 17-21.

[31] 馬東豪, 張佳寶, 吳忠東, 等. 電阻率成像法在土壤水文學研究中的應用及進展[J]. 土壤學報, 2014, 51(3): 439-447.

MA Donghao, ZHANG Jiabao, WU Zhongdong, et al. Application of electrical resistivity tomography to study on soil hydrology and its advance[J]. Acta Pedologica Sinica, 2014, 51(3): 439-447.

[32] LOKE M H. Electrical Imaging Surveys for Environmental and Engineering Studies, a practical guide to 2-D and 3-D surveys[J/OL]. Heritage Geophysics, 2000. http://www.terrajp.co.jp/lokenote.pdf.

[33] JOHONSON Timothy C, SLATER Lee, NTARLAGIANNIS Dimitrios, et al. Monitoring groundwater-surface water interaction using time-series and time-frequency analysis of transient three-dimensional electrical resistivity changes[J]. Water Resources Research, 2012, 48(7): 506.

[34] BRILLANTE Luca, BOIS Benjamin, MATHIEU Olivier, et al. Monitoring soil volume wetness in heterogeneous soils by electrical resistivity. A field-based pedotransfer function[J]. Journal of Hydrology (Amsterdam), 2014, 516: 56-66.

[35] WANG Yunqiang, SHAO Mingan, LIU Zhipeng, et al. Regional spatial pattern of deep soil water content and its influencing factors[J]. Hydrological Sciences Journal, 2012, 57(2): 265-281.

[36] 張翔, 李占斌, 李鵬. 淤地壩土壤水來源與分布研究[J]. 西安理工大學學報, 2013(1): 1-7.

ZHANG Xiang, LI Zhanbin, LI Peng. The study of soil water source and its distribution in soil-saving dam[J]. Journal of Xi’an University of Technology, 2013(1): 1-7.

[37] 鄒俊亮, 邵明安, 龔時慧. 不同植被和土壤類型下土壤水分剖面的分異[J]. 水土保持研究, 2011, 18(6): 12-17.

ZOU Junliang, SHAO Mingan, GONG Shihui. Effects of different vegetation and soil types on profile variability of soil moisture[J]. Research of Soil and Water Conservation, 2011, 18(6): 12-17.

[38] LUO Zhanbin, YONG Chenxu, FAN Jun, et al. Precipitation recharges the shallow groundwater of check dams in the loessial hilly and gully region of China[J]. Science of The Total Environment, 2020, 742: 140 625.

[39] 白樂, 殷會娟, 李昭悅, 等. 黃土丘陵區淤地壩作用下地下水位動態演化特征[J]. 人民黃河, 2021, 43(3): 74-77.

BAI Le, YIN Huijuan, LI Zhaoyue, et al. Dynamic evolution characteristics of groundwater level under the action of warping dam in loess hilly area[J]. Yellow River, 2021, 43(3): 74-77.

[40] 周立花, 延軍平, 徐小玲, 等. 黃土高原淤地壩對土壤水分及地表徑流的影響：以綏德縣辛店溝為例[J]. 干旱區資源與環境, 2006, 20(3): 112-115.

ZHOU Lihua, YAN Junping, XU Xiaoling, et al. Effect of soil saving dams on soil water content and runoff in loess plateau: Case study on Xindian watershed in Suide Country[J]. Journal of Arid Land Resources and Environment, 2006, 20(3): 112-115.

[41] 趙培培. 黃土高原小流域典型壩地土壤水分和泥沙空間分布特征[D]. 北京: 中國科學院大學, 2010.

ZHAO Peipei. Spatial distribution of soil water content and sediment in the dam farmlands in a small catchment of the Loess Plateau[D]. Beijing: University of Chinese Academy of Sciences, 2010.

[42] LIU Chun, LI Zhongwu, DONG Yuting, et al. Do land use change and check-dam construction affect a real estimate of soil carbon and nitrogen stocks on the Loess Plateau of China[J]. Ecological Engineering, 2017, 101: 220-226.

Variation of Soil Water over Slopes and Retained Lands in Loess Region: Investigated Using Electrical Resistivity Tomography

DUAN Guoxiu1, JIA Xiaoxu1,3*, BAI Xiao2, LIU Chenggong1,3, WEI Xiaorong2,3

(1. Key laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research, Beijing 100101, China; 2. Northwest A&F University, Yangling 712100, China;3. College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100190, China)

【Objective】Crop growth and ecological functions in arid and semi-arid loess regions in northwestern China are limited not only by topsoil water directly available to crop but by deep soil water which functions as a reservoir. Slopes and lands formed by artificial retaining are two typical geographical units in small watersheds in the loess plateau. This paper presents a method to estimate soil water distribution in them up to 12.5 m deep.【Method】The method was based on electrical resistivity tomography (ERT). We measured electrical resistivity of the soil in typical slope and retained land in the small Liudaogou watershed in northern Shaanxi province. Using the measured soil volumetric water content, a power function relating the resistivity to soil water content was established, from which we calculated water distribution and water storage in the 0～12.5 m profile in both the slope and the retained land.【Result】The distribution of electrical resistivity over the slope was significantly affected by slope position, with the resistivity decreasing gradually from the slope top to the slope toe. Such changes were associated with vegetation consumption of the topsoil water and redistribution of the infiltrated precipitation over the slope. The vertical distribution of the resistivity in the retained land showed a high-low-high variation; this was also related to root water uptake from the topsoil and precipitation recharge to the subsoil. ERT inversion showed co-existence of a saturated zone and an unsaturated zone in the retained land, while the slope was partly saturated. In a 1.8×104m2of retained land, there was 1.49×104m3of water in the 0～6.5 m unsaturated layer, and 5.10×104m3of water in the saturated layer. In a same area but on the slope, there was only 2.59×104m3of water in the 0～12.5 m soil layer.【Conclusion】Retained land contains more water than slope, and the deep soil in it functions as a reservoir banking infiltration water in wet seasons. ERT is suitable for measuring spatiotemporal variation in soil moisture in both slopes and flatten plains in the loess plateau.

loess plateau;electrical resistivity tomography (ERT); soil resistivity; slope land; dam land; soil water content

段國秀, 賈小旭, 白曉, 等. 基于ERT的黃土高原小流域典型坡地和壩地土壤儲水量反演[J]. 灌溉排水學報, 2023, 42(3): 104-111.

DUAN Guoxiu, JIA Xiaoxu, BAI Xiao, et al. Variation of Soil Water over Slopes and Retained Lands in Loess Region: Investigated Using Electrical Resistivity Tomography[J]. Journal of Irrigation and Drainage, 2023, 42(3): 104-111.

1672 - 3317（2023）03 - 0104 - 08

P343.9

A

10.13522/j.cnki.ggps.2022274

2022-05-11

國家自然科學基金項目（42022048）

段國秀（1996-），女。碩士研究生，主要從事土壤物理與水文生態研究。E-mail: duangx@nwafu.edu.cn

賈小旭（1985-），男。研究員，主要從事土壤物理與農業生態研究。E-mail: jiaxx@igsnrr.ac.cn

責任編輯：韓 洋