趙文婷　姜曉晗　李萌萌　焦菊英　嚴晰芹　祁泓錕

摘要：植被對暴雨的攔蓄是調(diào)控全球變暖進程中頻發(fā)的雨洪災害的重要手段。為探究黃丘區(qū)大規(guī)模退耕后自然植被坡面對暴雨的攔蓄作用，針對陜西省延安市安塞區(qū)坊塌小流域2019—2022年全部發(fā)生的6場暴雨事件，使用土壤水分自動監(jiān)測系統(tǒng)對1個灌叢和3個草本群落坡面10 m土層進行逐小時監(jiān)測。分析結(jié)果表明：暴雨對草本群落最大補給深度可達200 cm，但對灌叢群落僅為100 cm;灌叢的攔蓄效果優(yōu)于草本群落，草本群落的攔蓄效果隨退耕年限增長而提高;自然灌草植被坡面對低雨量短歷時暴雨的攔蓄率可達60.8%，而對高雨量長歷時暴雨僅可攔蓄40.8%。因此，建議在黃丘區(qū)合理保護和布局灌-草生態(tài)系統(tǒng)，以提升植被-土壤對暴雨的容蓄減災作用。

關(guān)鍵詞：暴雨事件;灌草植被;攔蓄作用;土壤入滲;黃丘區(qū)

中圖分類號：S152

文獻標志碼：A

文章編號：1001-6791（2023）05-0731-13

暴雨是水土流失及洪水災害的主要誘因，嚴重威脅生態(tài)環(huán)境安全［1］。植物-土壤生態(tài)系統(tǒng)能夠通過冠層和枯枝落葉層的攔擋及“土壤水庫”等方式來降低雨洪災害風險［2］。自退耕還林/草工程實施以來，黃丘區(qū)植被覆蓋度顯著提高［3］，土壤水分的入滲性能和蓄持能力顯著提升［4］。然而，在全球變暖導致暴雨發(fā)生頻率大幅增加的背景下［5］，退耕恢復植被對暴雨的攔蓄作用及抵御能力如何還不是很明確，但這對提升植被-土壤生態(tài)系統(tǒng)對暴雨的減災作用、優(yōu)化植被恢復布局具有重要的指導意義。因此，暴雨條件下不同退耕植被對水分攔蓄和土壤入滲能力的影響仍需進一步研究［6］，明確退耕恢復植被對暴雨產(chǎn)流的影響及其對降雨的調(diào)控，以盡最大可能提升植被對降雨的調(diào)控作用，避免降雨-洪水災害鏈的形成。

近年來，植被對降雨的攔蓄作用受到國際廣泛關(guān)注［7-8］。然而，由于退耕還林/草工程為中國特色生態(tài)恢復工程，故國內(nèi)研究聚焦于退耕造成的影響，多集中在不同植被類型與恢復方式對降雨的響應(yīng)以及干化土層形成特征與恢復等［9］，而針對不同暴雨條件下退耕植被蓄水減洪能力的研究較少。在黃丘區(qū)，一般降雨的土壤補給深度為100～200 cm，僅有較大的降雨事件才能滲透至深層土壤［10］。如在2013年7月延安長歷時極端降雨事件中（總雨量高達792.9 mm），自然退耕草本群落的補給深度可達500 cm，但人工退耕林地僅可補給至300 cm［11］。雖然極端降雨事件發(fā)生頻率相對較低，但降雨強度在50～100 mm/d的暴雨即可導致黃丘區(qū)在短時間內(nèi)產(chǎn)生地表徑流，在造成雨水資源浪費的同時，也產(chǎn)生了嚴重的土壤侵蝕［11］。

人工植被掠奪性地利用有限的土壤水資源，形成了明顯的土壤干層，影響生態(tài)系統(tǒng)的可持續(xù)發(fā)展［4，12］。而自然恢復的植被由于與生態(tài)環(huán)境相適宜而具有較高的生態(tài)功能［13］。因此，本文采用多探頭土壤水分自動監(jiān)測儀進行長期連續(xù)定位觀測，對2019—2022年陜西省延安市安塞區(qū)坊塌小流域不同退耕年限的自然恢復植被坡面土壤水分對暴雨的動態(tài)響應(yīng)特征進行分析，以揭示黃丘區(qū)坡面植被自然恢復后對暴雨的攔蓄作用，為區(qū)域合理利用土壤水資源和優(yōu)化植被布局提供科學支撐。

1 材料與方法

1.1 研究區(qū)概況

第5期趙文婷，等：黃丘區(qū)自然植被對暴雨的攔蓄作用

水科學進展第34卷

陜西省延安市安塞區(qū)坊塌小流域位于黃丘區(qū)腹地，屬第二副區(qū)，具有溝壑縱橫、地形破碎等典型黃丘區(qū)特征，也曾遭受嚴重的水土流失問題。該流域在1999年被選作退耕還林/草工程的“先行者”進行示范及推廣［14］，擁有耕地、果園、林地、草地等多種土地利用類型，有壩系、水平梯田等水土保持工程措施，布設(shè)有大量徑流小區(qū)和土壤水分監(jiān)測儀器，且廣泛開展了植被恢復、土壤恢復、土壤碳循環(huán)、侵蝕產(chǎn)沙、泥沙連通性等研究［15-17］，已成為典型的農(nóng)業(yè)生態(tài)小流域和研究小流域。因此，坊塌小流域具有成為后續(xù)黃丘區(qū)植被恢復建設(shè)“由綠變穩(wěn)”、生態(tài)保護與暴雨洪水減災工作的一把“金釘子”的潛力。坊塌小流域地理坐標介于36°47′19″N—36°49′35″N、109°14′40″E—109°17′09″E，面積為10.5 km2，海拔為997～1 731 m;屬于暖溫帶半濕潤氣候向半干旱氣候過渡區(qū)，以水力侵蝕為主，多年平均氣溫為8.8 ℃，年平均降水量為542.6 mm，降雨集中在7—10月且多暴雨，多年平均無霜期為160 d［18］;屬于暖溫帶落葉闊葉林向暖溫帶草原區(qū)過渡的森林草原區(qū)，土壤類型以黃綿土為主，黃綿土顆粒組成以細砂粒和粉粒為主，疏松多孔，滲透性較強［19］。該流域優(yōu)勢植物種與黃丘區(qū)也基本相似［20］，主要有刺槐（Robinia pesudoacacia）、檸條（Caragana korshinskii）、狼牙刺（Sophora davidii）、鐵桿蒿（Artemisia gmelinii）、草木樨狀黃芪（Astragalus melilotoides）、白羊草（Bothriochloa ischaemum）、長芒草（Stipa bungeana）、茭蒿（Artemisia giraldii）和達烏里胡枝子（Lespedeza daurica）等［21］，典型自然植物群落主要為鐵桿蒿群系和狼牙刺群落。

1.2 樣地選擇、土壤水分與降雨監(jiān)測

自然植被恢復演替時間長，而完整的固定樣地演替觀測很少，因此“時空互代”法被廣泛使用在自然植被恢復研究中。黃丘區(qū)的植被恢復演替需經(jīng)歷一年生草本群落、多年生蒿禾類草本群落、灌木群落和喬木群落階段［20，22］。退耕還林/草工程已大規(guī)模實行24 a，已由一年生草本群落演替到多年生草本群落及灌木群落階段，但尚未演替到喬木群落階段。其中，半灌木鐵桿蒿是植被演替到多年生蒿禾類草本階段的主要群落的建群種或亞優(yōu)勢種;而灌木狼牙刺在多年生篙類群落階段侵入，因處干燥貧瘠生境，基質(zhì)恢復慢，可存留20～40 a，是黃丘區(qū)干旱生境中的主要植被類型［23-24］。依據(jù)課題組于2022年對黃丘區(qū)森林帶、森林草原帶和草原帶的42個坡面物種組成的調(diào)查結(jié)果，鐵桿蒿群系出現(xiàn)頻率高達95.2%，狼牙刺群落約占灌木群落的60.0%?？梢?，鐵桿蒿群系和狼牙刺群落是該區(qū)域的典型植物群落。同時，2019—2022年坊塌小流域的79個草本樣方內(nèi)鐵桿蒿群系出現(xiàn)頻率為81.0%。其中，鐵桿蒿群落占鐵桿蒿群系的36.9%，鐵桿蒿+草木樨狀黃芪群落占32.7%，鐵桿蒿+長芒草群落占15.6%。在26個灌木樣方中，狼牙刺群落占42.3%?？梢姡@4種植被群落類型在坊塌小流域也具有代表性。因此，本研究選取了立地條件相似且具有這4種典型植物群落的連貫坡面，通過多次野外考證與農(nóng)戶走訪確定退耕年限，最終選擇截至2022年退耕47 a的鐵桿蒿群落、32 a的鐵桿蒿+草木樨狀黃芪群落、12 a的鐵桿蒿+長芒草群落和47 a的灌叢狼牙刺群落作為研究對象，從不同演替階段和不同群落結(jié)構(gòu)差異角度進行對比分析，以深入了解退耕還林/草工程實施后不同恢復年限間同一群系的自然草本群落、灌木群落和草本群落對土壤水分調(diào)控效果的差異性，以評估自然退耕坡面對暴雨的抵御效果，為后續(xù)制定合理的植被建設(shè)方向提供科學參考，以實現(xiàn)退耕還林/草工程成果的可持續(xù)性。

退耕坡面上各布設(shè)有1個徑流小區(qū)（圖1，表1），規(guī)格為2.5 m×10.0 m。同時，坡面均設(shè)有CS650-CR1000土壤水分自動監(jiān)測系統(tǒng)（美國Campbell公司），可對0～1 000 cm土層深度的土壤體積含水量進行長期原位、多層次的土壤水分動態(tài)監(jiān)測，水分探頭安裝深度分別為10、30、50、70、100、200、300、500、700和1 000 cm，監(jiān)測頻率為1次/h，體積含水率變化范圍為4.0%～34.1%。土壤水分觀測年限為2019—2022年，4 a的小時尺度0～1 000 cm土層的土壤水分監(jiān)測數(shù)據(jù)精細地刻畫了降雨過程中土壤水分的動態(tài)變化情況，為退耕還林/草工程生態(tài)效益評價以及雨洪災害預防提供了寶貴的數(shù)據(jù)支撐。

降雨數(shù)據(jù)通過布置在本研究退耕植被坡面的翻斗式自記雨量計獲取，監(jiān)測頻率為4次/h。降雨過程依據(jù)最小降雨間歇6 h進行劃分，即降雨時間間隔小于6 h則為1次降雨過程，反之則為2次［25］。并對降雨事件的最大30 min雨強（I30）和平均雨強等參數(shù)進行計算。

1.3 土壤對降雨攔蓄作用的表征指標

土壤對降雨的攔蓄即為降雨的就地入滲［26］，其中降雨進入土壤后在土層內(nèi)的分布及其隨時間的變化是土壤攔蓄降雨的重要過程，故本研究選用土壤入滲量、土壤水分入滲率/消退率對暴雨條件下土壤的攔蓄過程進行描述，并采用土壤水分補給率指示土壤對暴雨的攔蓄效益。

1.4 數(shù)據(jù)分析

采用CS650-CR1000土壤水分自動監(jiān)測系統(tǒng)儀器的配套軟件LoggerNet和Microsoft Excel 2016進行土壤水分數(shù)據(jù)提取及分析;基于SPSS 26.0開展土壤水分指標與降雨指標、雨前土壤含水量及土壤深度的K-S正態(tài)性檢驗和Spearman相關(guān)性分析，以探究不同指標之間的相關(guān)關(guān)系;本文中的圖件均通過Origin 2022和Adobe Photoshop 2020繪制。

2 結(jié)果分析

2.1 暴雨事件特征

根據(jù)國家氣象局降雨等級劃分標準，研究區(qū)在2019—2022年共發(fā)生6場暴雨事件。其中，暴雨5場，大暴雨1場，累計降水量為666.5 mm，占4 a內(nèi)共發(fā)生的244次降雨事件總降水量的63.7%。同時，這些暴雨事件的徑流深累計35.8 mm，占比高達總徑流深的67.3%。這6場暴雨事件的降雨與產(chǎn)流特征如圖2和表2所示。根據(jù)各場暴雨的歷時和雨量特征將以上6場暴雨劃分為長歷時、高雨量暴雨（暴雨3和5）和短歷時、低雨量暴雨（暴雨1、2、4和6）。同時結(jié)合強度劃分標準［30］將其進一步劃分為長歷時、中雨強、高雨量暴雨（暴雨3和5），短歷時、中雨強、低雨量暴雨（暴雨2和6）以及短歷時、高雨強、低雨量暴雨（暴雨1和4）。

暴雨事件對植物群落坡面的補給深度最深可達200 cm。其中，暴雨對草本群落坡面的補給深度為70～200 cm，而對灌叢坡面的補給深度僅為70～100 cm。且不同暴雨類型對植物群落坡面的補給深度也存在差異。在短歷時、低雨量的暴雨事件中，各植物群落充分入滲后的平均補給深度為70～100 cm，其中，土壤水分入滲深度（d）對暴雨1的響應(yīng)最小，經(jīng)充分入滲，除鐵桿蒿群落的補給深度達到100 cm，其余3個群落的最大入滲深度均為70 cm;而對于長歷時、高雨量暴雨事件，各群落土壤水分補給深度為100～200 cm（圖3，其中，W₀為雨前土壤含水量；S₀為雨后0 h土壤入滲量；S₉₆為雨后96 h土壤入滲量），其中，狼牙刺群落的補給深度為100 cm，其余群落均為200 cm，這表明長歷時、高雨量的暴雨事件對植物群落坡面的入滲深度更深。

植被群落坡面土壤水分入滲總量對暴雨的響應(yīng)較為復雜，主要受到降水量和I₃₀的主導（p<0.05），而平均雨強對土壤水分入滲總量的影響卻并不顯著（p>0.05）（表3）。各植被群落坡面在短歷時、低雨量暴雨中的平均土壤水分入滲總量為48.7 mm（圖3），其中，暴雨1和暴雨2的降水量相當，但暴雨2的I₃₀是暴雨1的1.4倍，理應(yīng)暴雨1的土壤入滲總量高于暴雨2，但由于暴雨2的雨前土壤含水量高于暴雨1，其各群落坡面的降雨入滲總量較暴雨1低6.06%～16.00%，平均入滲深度向下延伸30～50 cm。在長歷時、高雨量暴雨事件中，各植被群落坡面的平均土壤水分入滲總量為72.9 mm，遠高于短歷時、低雨量暴雨事件（圖3）。其中，暴雨5的降水量低于暴雨3，但平均雨前土壤含水量較暴雨3多43.1 mm，其在植物群落坡面的平均土壤水分入滲量較暴雨3少50.3 mm。可見，除了降水量和I₃₀外，雨前土壤含水量也會對土壤入滲產(chǎn)生影響。

土壤水分入滲率（V）隨著平均雨強和I₃₀的增加而增加（p<0.01）（表3）。其中，短歷時、高雨強、低雨量暴雨的土壤水分入滲率為2.11 mm/h;其次為短歷時、中雨強、低雨量暴雨，V=1.11 mm/h;而長歷時、中雨強、大雨量暴雨的V僅為0.96 mm/h（圖4）。

灌叢群落的入滲能力強于草本群落，而草本群落的入滲能力隨退耕年限的增加而增加。暴雨條件下灌叢和草本植物群落坡面0～200 cm的土壤水分入滲總量介于22.1～102.9 mm（圖3）。其中，灌叢狼牙刺群落坡面對暴雨響應(yīng)最為強烈，入滲總量均值高達61.4 mm;而草本群落鐵桿蒿+長芒草、鐵桿蒿+草木樨狀黃芪和鐵桿蒿群落坡面的入滲總量分別為56.6、51.3和65.7 mm。不同植物群落坡面的土壤水分入滲率也展現(xiàn)出同樣趨勢，灌叢群落較草本群落的土壤水分入滲率高10.9%～41.2%，鐵桿蒿群落較鐵桿蒿+長芒草群落和鐵桿蒿+草木樨狀黃芪群落的土壤水分入滲率分別增加0.50和0.15 mm/h（圖4）。

2.3 暴雨條件下不同植物群落坡面土壤水分的補給特征

降雨結(jié)束初期退耕時間長的群落坡面對暴雨的攔蓄作用較退耕時間短的群落更強，但坡面的攔蓄作用會隨時間的推移呈降低趨勢，而且退耕年限長的植物群落坡面對暴雨的攔截作用隨時間下降的幅度較退耕年限短的群落更大。因為植物群落坡面的土壤水分補給率與暴雨事件的降雨歷時呈現(xiàn)極顯著負相關(guān)關(guān)系（p<0.01）（表3）。在短歷時、低雨量暴雨中，退耕47 a的鐵桿蒿群落和狼牙刺群落坡面雨后0 h的平均土壤水分補給率大于退耕時間較短的鐵桿蒿+草木樨狀黃芪和鐵桿蒿+長芒草群落;而在降雨結(jié)束96 h后，退耕年限較長的狼牙刺群落和鐵桿蒿群落的平均土壤水分補給率分別降低了26.3%、21.1%，但退耕年限較短的群落僅降低了9.3%和4.0%（表4）。

高雨強暴雨較中雨強暴雨的平均土壤水分補給能力更強。其中，短歷時、高雨強、低雨量暴雨的雨后0 h平均土壤水分補給率為64.8%，雨后96 h降低至53.8%；而短歷時、中雨強、低雨量暴雨平均土壤水分補給率為56.8%，雨后96 h縮減為33.3%；在長歷時、中雨強、高雨量暴雨中，雨后0 h土壤水分補給率為40.8%，遠低于短歷時、低雨量暴雨條件下的土壤水分補給率（表4）。

不同植物群落各土層的土壤水分消退率（R）隨著土壤深度的增加快速下降（圖5）。其中，狼牙刺、鐵桿蒿群落10～50 cm土層土壤水分消退率相對較高，平均值分別為0.13 mm/h和0.11 mm/h;經(jīng)96 h充分入滲后，土壤水分向深層運移，70～200 cm土層的土壤水分增加，所以70～100 cm土層的土壤水分消退率較低，平均值分別為0.04 mm/h和0.03 mm/h。鐵桿蒿+草木樨狀黃芪、鐵桿蒿+長芒草群落10～30 cm土層的土壤水分消退率相對較高，平均值均為0.11 mm/h，而50～100 cm土層的土壤水分消退率平均值均降至0.03 mm/h。這4種群落坡面各土層的平均土壤水分消退率分別為0.11、0.09、0.08、0.08 mm/h。對不同暴雨事件而言，長歷時、高雨量暴雨事件的土壤水分消退率平均值為0.06 mm/h;而短歷時、低雨量暴雨事件的消退率平均值為0.09 mm/h，其中的高雨強暴雨事件的土壤水分消退率高達0.11 mm/h，特別是暴雨1，土壤水分消退率高達0.15mm/h，各植物群落在10～30 cm土層的土壤水分消退率為其他場次暴雨相應(yīng)土層土壤水分消退率的0.5～2.3倍。

3 討論

黃丘區(qū)坊塌小流域自實行退耕還林/草措施以來，植被覆蓋度及地上生物量大幅增加，土壤水分物理性質(zhì)和結(jié)構(gòu)得以有效改善［31］，但退耕后土壤水分對暴雨的響應(yīng)研究仍有不足，所以本研究采用由自計雨量計獲取的2019—2022年流域內(nèi)發(fā)生的6場暴雨進行了研究，雖然暴雨次數(shù)較少，但通過對2007—2020年位于黃丘區(qū)的綏德縣、志丹縣和吳起縣的暴雨事件統(tǒng)計發(fā)現(xiàn)，這些站點在14 a內(nèi)均僅發(fā)生了7次暴雨和1次大暴雨，而在2017—2020年4 a內(nèi)發(fā)生了3～4次暴雨和1次大暴雨?？梢?，黃丘區(qū)暴雨事件的出現(xiàn)愈加頻繁。因此，坊塌小流域在4 a監(jiān)測到的5場暴雨和1場大暴雨事件，在數(shù)量上具有典型性。同時，6場暴雨事件可被劃分為“長歷時、中雨強、高雨量暴雨”、“短歷時、中雨強、低雨量暴雨”以及“短歷時、高雨強、低雨量暴雨”3種黃丘區(qū)常出現(xiàn)的類型［32］，且均有2場重復。因此，雖然本研究時段中僅發(fā)生了6場暴雨事件，但其對黃丘區(qū)暴雨事件的發(fā)生是具有代表性的。

當前，黃丘區(qū)的自然植被坡面土壤水分對降雨的響應(yīng)過程較為復雜，主要受降水量、I₃₀、降雨歷時的主導調(diào)控。特別是I₃₀較平均雨強影響的解釋力更為顯著，因為平均雨強容易掩蓋自然降雨過程中的“強度爆發(fā)”與“間歇性”等降雨波動特征［27］。隨著降雨強度增加，土壤攔蓄率隨之增加，且在4.8和7.9 mm/h的高雨強暴雨事件中坡面土壤水分入滲率達到2.11 mm/h，遠超當?shù)氐耐寥婪€(wěn)滲率1.15～1.30 mm/h［33］，但并未發(fā)生超滲產(chǎn)流事件，可能黃丘區(qū)的植被建設(shè)已使降雨-產(chǎn)流模式發(fā)生自“超滲產(chǎn)流”向“蓄滿產(chǎn)流”的轉(zhuǎn)變［34］。在降雨過程中，表層土壤水分先受到補充;由于不同土層間土水勢的差異，水分逐漸向深層運移;但隨著降雨歷時的增加，不同土層間的土水勢差值逐漸減小，導致各自然植物群落坡面的土壤水分入滲率和補給率隨降雨歷時的延伸出現(xiàn)不同程度的下降［35］。

水分在土壤剖面內(nèi)的運移及分配受到植被類型、土壤質(zhì)地和雨前土壤含水量等下墊面因素的調(diào)控。在黃丘區(qū)，對土壤水分的補給深度基本集中在地下70～200 cm，其中草地坡面的補給深度最深可達200 cm。但由于灌叢蒸散量大，表層土壤水分消耗劇烈，因此灌叢群落坡面100 cm以下土層的水分難以得到補充［36］。此外，植被恢復與重建過程中，土壤有機質(zhì)、孔隙度、團粒結(jié)構(gòu)及水穩(wěn)性團聚體增加，土壤重度降低［37-38］，土壤質(zhì)地發(fā)生改善，所以呈現(xiàn)退耕年限久的植被對暴雨的攔蓄能力強和灌叢群落坡面較草本群落坡面對暴雨的攔蓄能力更強的規(guī)律。具體在本研究中，草本群落的土壤有機質(zhì)含量自退耕12 a的7.62 g/kg提升至12.03 g/kg，水穩(wěn)性團聚體含量由50.64%提升至66.12%，毛管孔隙度由50.81%提升至52.47%，土壤重度因此降低了0.04 g/cm（表1），促使水分的下滲與貯蓄［29，31］。故草本群落中，退耕47 a的入滲量較退耕32 a和12 a分別提升14.0%和22.0%。灌木群落階段的毛管孔隙度、有機質(zhì)含量比同退耕年限的草本群落更高，土壤重度更低，故而擁有更強的入滲能力［27，39］。同時，退耕年限長的狼牙刺和鐵桿蒿群落坡面由于土壤孔隙度的增加，對其土層內(nèi)入滲的土壤水分再分配迅速，也可迅速排空大部分進入淺土層內(nèi)的降雨，如若短期內(nèi)再次降雨，仍可容納較多的雨水。此外，來自降雨的入滲水分首先要補足表層因植被導致的水分虧缺才能繼續(xù)下滲，所以雨前土壤含水量高的坡面會較含水量低的坡面先補足表層土壤水分進行深層入滲，導致入滲深度更深。然而，坡面在補足表層水分后，濕潤前沿的水力梯度也隨之減小，下滲的過程相對更為困難，因此入滲量更少［27］。

綜上，退耕植被恢復的確有利于增強土壤持蓄水分的能力，減少徑流匯集，錯峰消洪，促成“偽蓄滿產(chǎn)流”。然而，土壤只能攔蓄部分降雨，其產(chǎn)生的徑流仍可誘發(fā)土壤侵蝕，增大洪水風險。因此，在黃丘區(qū)仍應(yīng)在加強自然恢復植被坡面攔蓄能力的同時，合理配置引排措施，以期更好地發(fā)揮植被-土壤對暴雨的攔蓄作用。

4 結(jié)論

本研究通過對陜西省延安市安塞區(qū)坊塌小流域在2019—2021年所發(fā)生的6次暴雨事件進行逐小時土壤水分監(jiān)測，定量分析了不同自然灌草植被坡面在不同暴雨類型下的攔蓄效果，探究了黃丘區(qū)退耕自然灌草植被坡面對暴雨的攔蓄作用。主要結(jié)論如下：

（1） 6場暴雨對草本群落的補給深度最深可達200 cm，但對灌叢狼牙刺群落的補給深度最深僅為100 cm。灌叢狼牙刺群落的土壤水分入滲量、入滲率均較草本群落高;草本群落中，退耕年限長的鐵桿蒿群落的入滲量、入滲率均高于退耕年限較短的鐵桿蒿+草木樨狀黃芪群落和鐵桿蒿+長芒草群落。

（2） 灌木狼牙刺群落的土壤水分入滲總量、土壤水分入滲率以及土壤水分補給率較草本群落更優(yōu)，即灌木狼牙刺群落對暴雨的攔蓄效果更強，草本群落對暴雨的攔蓄效果隨退耕年限的增長而提升;但雨前土壤含水量較高會降低植物群落對暴雨的攔蓄作用。

（3） 自然灌草植被坡面在短歷時、高雨強、低雨量暴雨前期可平均攔蓄約64.8%的降雨，在短歷時、中雨強、低雨量暴雨條件下可攔蓄56.8%，但長歷時、中雨強、高雨量暴雨僅可攔蓄40.8%。

參考文獻：

［1］汪星，張敬曉，呂望，等.黃土丘陵區(qū)深層干化土壤對降雨入滲的響應(yīng)［J］.干旱地區(qū)農(nóng)業(yè)研究，2021，39（4）：29-38，63.（WANG X，ZHANG J X，LYU W，et al.Infiltration response of deep dry soil to rainfall in the loess hilly region［J］.Agricultural Research in the Arid Areas，2021，39（4）：29-38，63.（in Chinese））

［2］趙炯昌，衛(wèi)偉，段興武.模擬降雨下黃土坡面水沙過程對3種灌草植被垂直結(jié)構(gòu)變化的響應(yīng)［J］.生態(tài)學報，2021，41（21）：8602-8611.（ZHAO J C，WEI W，DUAN X W.Response of the runoff and sediment process on loess slope to the vertical structure changes of three shrub and grass vegetations under simulated rainfall［J］.Acta Ecologica Sinica，2021，41（21）：8602-8611.（in Chinese））

［3］MAO R C，SONG J X，TANG B，et al.Vegetation variation regulates soil moisture sensitivity to climate change on the Loess Plateau［J］.Journal of Hydrology，2023，617：128763.

［4］HAN L，CHANG Y Q，CHEN R，et al.Response of soil moisture to vegetation and trade-off analysis in the hilly area of the Loess Plateau，China［J］.Ecological Indicators，2022，142：109273.

［5］GRNDEMANN G J，van de GIESEN N，BRUNNER L，et al.Rarest rainfall events will see the greatest relative increase in magnitude under future climate change［J］.Communications Earth & Environment，2022，3：235.

［6］金釗.黃土高原小流域退耕還林還草的生態(tài)水文效應(yīng)與可持續(xù)性［J］.地球環(huán)境學報，2022，13（2）：121-131.（JIN Z.Ecohydrological effects and sustainability of the grain for green project on the Chinese Loess Plateau［J］.Journal of Earth Environment，2022，13（2）：121-131.（in Chinese））

［7］PARONUZZI P，BOLLA A.Rainfall infiltration and slope stability of alpine colluvial terraces subject to storms （NE Italy）［J］.Engineering Geology，2023，323：107199.

［8］BAHDDOU S，OTTEN W，WHALLEY W R，et al.Changes in soil surface properties under simulated rainfall and the effect of surface roughness on runoff，infiltration and soil loss［J］.Geoderma，2023，431：116341.

［9］白曉，賈小旭，邵明安，等.黃土高原北部土地利用變化對長期土壤水分平衡影響模擬［J］.水科學進展，2021，32（1）：109-119.（BAI X，JIA X X，SHAO M A，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.（in Chinese））

［10］YU X N，HUANG Y M，LI E G，et al.Effects of rainfall and vegetation to soil water input and output processes in the Mu Us Sandy Land，Northwest China［J］.Catena，2018，161：96-103.

［11］尹秋龍，焦菊英，寇萌.極端強降雨條件下黃丘區(qū)不同植被類型土壤水分特征［J］.自然資源學報，2015，30（3）：459-469.（YIN Q L，JIAO J Y，KOU M.The soil moisture characteristics under different vegetation types after extremely heavy rainfall on the hilly-gullied Loess Plateau［J］.Journal of Natural Resources，2015，30（3）：459-469.（in Chinese））

［12］GE F C，XU M X，LI B B，et al.Afforestation reduced the deep profile soil water sustainability on the semiarid Loess Plateau［J］.Forest Ecology and Management，2023，544：121240.

［13］HUANG L M，SHAO M A.Advances and perspectives on soil water research in China′s Loess Plateau［J］.Earth-Science Reviews，2019，199：102962.

［14］李宗善，楊磊，王國梁，等.黃土高原水土流失治理現(xiàn)狀、問題及對策［J］.生態(tài)學報，2019，39（20）：7398-7409.（LI Z S，YANG L，WANG G L，et al.The management of soil and water conservation in the Loess Plateau of China：present situations，problems，and counter-solutions［J］.Acta Ecologica Sinica，2019，39（20）：7398-7409.（in Chinese））

［15］XU Q，LI M M，JIANG X H，et al.Response of rill erosion to rainfall types and maintenance on the Loess Plateau：implications for road erosion control［J］.Catena，2022，219：106642.

［16］YAN X Q，JIAO J Y，LI M M，et al.Lateral sediment connectivity of landslides occurred under a heavy rainstorm and its influence on sediment yield of slope-channel cascade on the Loess Plateau［J］.Catena，2022，216：106378.

［17］SUN P C，PAN Y W，WU Y P，et al.The effects of land use and cover changes on lateral carbon losses from an ungagged headwater basin on the Chinese Loess Plateau［J］.Journal of Hydrology，2023，623：129751.

［18］李玉進，胡澍，焦菊英，等.黃土丘陵區(qū)不同侵蝕環(huán)境下土壤有機碳對植被恢復的響應(yīng)［J］.生態(tài)學報，2017，37（12）：4100-4107.（LI Y J，HU S，JIAO J Y，et al.Response of soil organic carbon to vegetation restoration in different erosion environments in the hilly-gullied region of the Loess Plateau［J］.Acta Ecologica Sinica，2017，37（12）：4100-4107.（in Chinese））

［19］高海東，李占斌，李鵬，等.基于土壤侵蝕控制度的黃土高原水土流失治理潛力研究［J］.地理學報，2015，70（9）：1503-1515.（GAO H D，LI Z B，LI P，et al.The capacity of soil loss control in the Loess Plateau based on soil erosion control degree［J］.Acta Geographica Sinica，2015，70（9）：1503-1515.（in Chinese））

［20］焦菊英，張振國，賈燕鋒，等.陜北丘陵溝壑區(qū)撂荒地自然恢復植被的組成結(jié)構(gòu)與數(shù)量分類［J］.生態(tài)學報，2008，28（7）：2981-2997.（JIAO J Y，ZHANG Z G，JIA Y F，et al.Species composition and classification of natural vegetation in the abandoned lands of the hilly-gullied region of North Shaanxi Province［J］.Acta Ecologica Sinica，2008，28（7）：2981-2997.（in Chinese））

［21］嚴方晨，焦菊英，曹斌挺，等.黃丘區(qū)撂荒地不同演替階段植物群落的土壤抗蝕性：以坊塌流域為例［J］.應(yīng)用生態(tài)學報，2016，27（1）：64-72.（YAN F C，JIAO J Y，CAO B T，et al.Soil anti-erodibility of abandoned lands during different succession stages of plant community in hilly-gullied region of the Loess Plateau：take Fangta small watershed as an example［J］.Chinese Journal of Applied Ecology，2016，27（1）：64-72.（in Chinese））

［22］曾光，楊勤科，張信寶.黃丘區(qū)退耕地植被自然恢復過程：以吳起縣雙樹溝流域為例［J］.中國水土保持科學，2008，6（3）：48-52.（ZENG G，YANG Q K，ZHANG X B.Natural succession process of vegetation on the abanded land in the hilly-gully region of Loess Plateau：a case study of Shuangshugou watershed in Wuqi County［J］.Science of Soil and Water Conservation，2008，6（3）：48-52.（in Chinese））

［23］朱志誠.陜北黃土高原灌木林的類型及其動態(tài)特性［J］.陜西林業(yè)科技，1992（1）：36-42.（ZHU Z C.The types and dynamic characteristics of shrubs on Loess Plateau in North Shaanxi［J］.Shaanxi Forest Science and Technology，1992（1）：36-42.（in Chinese））

［24］徐學華，張文輝，張慧，等.黃土丘陵區(qū)狼牙刺群落恢復過程中物種多樣性動態(tài)研究［J］.河北農(nóng)業(yè)大學學報，2006，29（3）：1-5.（XU X H，ZHANG W H，ZHANG H，et al.Study on the species diversity dynamics of Sophora davidii community in loess hilly region［J］.Journal of Agricultural University of Hebei，2006，29（3）：1-5.（in Chinese））

［25］DUNKERLEY D.The importance of incorporating rain intensity profiles in rainfall simulation studies of infiltration，runoff production，soil erosion，and related landsurface processes［J］.Journal of Hydrology，2021，603：126834.

［26］陳寶根，王仕琴，宋獻方.一維土壤水分運動模擬在土壤水分特征研究中的應(yīng)用：以華北平原衡水實驗站為例［J］.水文，2011，31（3）：64-70.（CHEN B G，WANG S Q，SONG X F.Application of 1-D soil water movement simulation in research on soil water characteristics：a case of Hengshui experimental station［J］.Journal of China Hydrology，2011，31（3）：64-70.（in Chinese））

［27］靳曉輝，陳敏建，閆龍，等.不同灌溉方式對降雨入滲深度的影響［J］.排灌機械工程學報，2018，36（11）：1168-1174.（JIN X H，CHEN M J，YAN L，et al.lnfluences of different irrigation methods on rainfall infiltration depth［J］.Journal of Drainage and Irrigation Machinery Engineering，2018，36（11）：1168-1174.（in Chinese））

［28］JIA Y H，SHAO M G.Temporal stability of soil water storage under four types of revegetation on the northern Loess Plateau of China［J］.Agricultural Water Management，2013，117：33-42.

［29］CUI Z，HUANG Z，WANG Y B，et al.Soil water deficit was effectively alleviated by higher water infiltration after the short-term forestland-to-farmland conversion in semi-arid area［J］.Journal of Hydrology，2022，610：127893.

［30］羅琳，王忠靜，劉曉燕，等.黃河流域中游典型支流汛期降雨特性變化分析［J］.水利學報，2013，44（7）：848-855.（LUO L，WANG Z J，LIU X Y，et al.Changes in characteristics of precipitation in flood season over five typical basins of middle reaches of the Yellow River in China［J］.Journal of Hydraulic Engineering，2013，44（7）：848-855.（in Chinese））

［31］晏清洪，原翠萍，雷廷武，等.降雨類型和水土保持對黃土區(qū)小流域水土流失的影響［J］.農(nóng)業(yè)機械學報，2014，45（2）：169-175.（YAN Q H，YUAN C P，LEI T W，et al.Effact of rainstorm patterns and soil erosion control practies on soil and water loss in small watershed on Loess Plateau［J］.Transactions of the Chinese Society for Agricultural Machinery，2014，45（2）：169-175.（in Chinese））

［32］蔣定生，黃國俊.黃土高原土壤入滲速率的研究［J］.土壤學報，1986，23（4）：299-305.（JIANG D S，HUANG G J.Study on the filtration rate of soils on the Loess Plateau of China［J］.Acta Pedologica Sinica，1986，23（4）：299-305.（in Chinese））

［33］JIN Z，GUO L，YU Y L，et al.Storm runoff generation in headwater catchments on the Chinese Loess Plateau after long-term vegetation rehabilitation［J］.Science of the Total Environment，2020，748：141375.

［34］梁越，焦菊英.黃土高原小流域產(chǎn)沙性降雨標準分析［J］.中國水土保持科學，2019，17（3）：8-14.（LIANG Y，JIAO J Y.Analysis on rainfall standards while sediment occurring in small watersheds on the Loess Plateau［J］.Science of Soil and Water Conservation，2019，17（3）：8-14.（in Chinese））

［35］PARIZEK B，ROSTAGNO C M，SOTTINI R.Soil erosion as affected by shrub encroachment in Northeastern Patagonia［J］.Journal of Range Management，2002，55（1）：43-48.

［36］楊新民，楊文治.灌木林地的水分平衡研究［J］.水土保持研究，1998，5（1）：109-118.（YANG X M，YANG W Z.Study on the water balance in shrub land［J］.Research of Soil and Water Conservation，1998，5（1）：109-118.（in Chinese））

［37］陳文媛，張少妮，華瑞，等.黃土丘陵區(qū)林草恢復進程中土壤入滲特征研究［J］.北京林業(yè)大學學報，2017，39（1）：62-69.（CHEN W Y，ZHANG S N，HUA R，et al.Effects of forestland and grassland restoration process on soil infiltration characteristics in loess hilly region［J］.Journal of Beijing Forestry University，2017，39（1）：62-69.（in Chinese））

［38］顧朝軍，穆興民，孫文義，等.極端暴雨洪水及侵蝕產(chǎn)沙對延河流域植被恢復響應(yīng)的比較研究［J］.自然資源學報，2017，32（10）：1755-1767.（GU C J，MU X M，SUN W Y，et al.Comparative analysis of the responses of rainstorm flood and sediment yield to vegetation rehabilitation in the Yanhe River basin［J］.Journal of Natural Resources，2017，32（10）：1755-1767.（in Chinese））

［39］肖婧，王兵.黃丘區(qū)坡面尺度撂荒草地入滲特征影響因素試驗研究［J］.水科學進展，2020，31（6）：820-831.（XIAO J，WANG B.Experimental study on the influencing factors of infiltration characteristics of abandoned grassland on slope scale in the loess hilly and gully region［J］.Advances in Water Science，2020，31（6）：820-831.（in Chinese））

Interception and storage of heavy rainfall by natural vegetations

in the loess hilly and gully area

The study is financially supported by the National Natural Science Foundation of China （No.U2243213;No.42077078）.

ZHAO Wenting JIANG Xiaohan LI Mengmeng JIAO Juying YAN Xiqin QI Hongkun

（1. State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau，Northwest A & F University，

Yangling 712100，China;

2. Institute of Soil and Water Conservation，Chinese Academy of Sciences and

Ministry of Water Resource，Yangling 712100，China）

Abstract：Vegetative trapping of heavy rainfall is an important mechanism for mitigating frequent rain-induced flood disasters in the context of global warming.To investigate the role of natural vegetation in trapping heavy rainfall following the widespread implementation of the “the Grain for Green” project in the Loess Hilly and Gully area，an automatic soil moisture monitoring system was employed to continuously monitor the soil moisture content in the top 10 m on one shrub-covered slope and three herbaceous slopes during six heavy rainfall events occurring in the Fangta watershed of the Ansai District，Shaanxi Province，from 2019 to 2022.The results showed that the maximum infiltration depth of the herbaceous slopes was approximately 200 cm during the heavy rainfall events，while for the shrub-covered slope，it was only 100 cm.Furthermore，the shrub-covered slope exhibited superior rainfall interception compared to the herbaceous slopes.The trapping efficianncy of the herbaceous slopes improved with each successive year of “Grain for Green” project implementation.Additionally，natural shrub-covered slopes displayed a trapping efficiency of up to 60.8% during low-intensity，short-duration rainfall events，and 40.8% during high-intensity，long-duration rainfall events.Therefore，it is advisable to strategically conserve and allocate shrub-herbaceous ecosystems within the Loess Hilly and Gully area to enhance the storage and mitigation capacity of vegetation and soil in response to heavy rainfall.

Key words：heavy rainfall event;scrub and herbaceous;heavy rainfall interception and storage;soil infiltration;the loess hilly and gully area