余恩旭, 张明芳,2, 蒋志魏, 徐亚莉, 邓诗宇

(1.电子科技大学 资源与环境学院, 成都 611731; 2.电子科技大学 长三角研究院(湖州), 浙江 湖州 313001)

土壤侵蚀会引起土地退化、土壤肥力下降、泥沙淤积和生态环境恶化等问题[1]。土壤侵蚀已成为全球的重要生态环境问题，尤其在地形起伏较大的山区尤为严重[2]。岷江上游流域地处青藏高原到四川盆地的过渡带，地形起伏大，地质破碎，夏季多暴雨，易发生滑坡、崩塌、泥石流等地质灾害；加之早期森林资源破坏严重，水土保持能力下降，局部土地沙化和石质化等，导致该地区水土流失严重，被列为国家级水土保持重点预防和防治区域之一[3]。与此同时，岷江上游流域作为长江上游重要水源涵养和水土保持区，是成都平原的生态围栏。该区域土壤侵蚀变化直接影响长江上游和成都平原的生态环境质量。因此，开展岷江上游土壤侵蚀的时空特征以及其成因研究十分必要，能够为相关部门治理水土流失、构筑长江上游生态安全屏障提供科学依据。

部分学者围绕岷江上游土壤侵蚀分布及驱动因素开展了相关研究[2,4-8]。根据龚雪梅的研究，1995—2014年岷江上游流域土壤侵蚀敏感性以轻度和中度为主，高度和极敏感区域主要分布在干旱河谷和高山裸土区[7]。姜琳等进一步指出2000—2010年岷江上游流域侵蚀较严重的区域主要分布在25°～45°坡度带、海拔<2 000 m和海拔4 000～5 000 m区域[2]。岷江上游流域土壤侵蚀主要受到地形、降水、土地利用、植被覆盖和人类活动的影响。据孟兆鑫等对岷江流域1995—2005年土壤侵蚀研究表明林地、草地面积之和与土壤侵蚀面积呈显著线性相关，且降水是土壤侵蚀的驱动因素之一[4]。刘金山等对岷江上游2015年土壤侵蚀风险评估结果表明中度侵蚀区域分布在水网附近，并指出土壤侵蚀与土地利用类型和人类活动关系密切[6]。但是上述研究多针对岷江上游流域某一年或短期土壤侵蚀空间分布特点和部分影响因素，尚不足以全面和系统地揭示区域土壤侵蚀的长期动态(趋势性和平稳性)、空间分布特征及其成因的区域分异规律。本研究应用USLE模型，结合ArcGIS，SPSS等专业软件，运用相关性分析、平稳性分析、趋势分析、偏相关等方法研究岷江上游流域2001—2017年土壤侵蚀时空变化的趋势性和平稳性特征、时空分异规律及其对驱动因子尤其是植被变化的响应。本研究可为岷江上游流域土壤侵蚀分类治理、提升植被水源涵养功能、构筑长江上游生态安全屏障提供科学依据。

1 研究区概况

岷江上游流域(102°35′—103°56′E，30°45′—33°09′N)位于青藏高原东南缘，横断山脉北端与川西高山峡谷东侧的结合部，包含黑水县、理县、汶川县、茂县、松潘县和都江堰市部分地区。流域面积约为2.287万km2，干流全长约337 km(图1)。该流域地处四川盆地丘陵山地过渡地带，西北高东南低，流域海拔在724～5 893 m，地形起伏较大，地貌类型较复杂。岷江上游水资源充沛，年均流量约158亿m3，丰水期(6—10月)流量占全年的 75%。受地形影响，气候垂直分异明显，沿海拔梯度分布亚热带、温带、亚高山和高山寒带等多种气候[9]。夏季温凉，冬春寒冷，年均气温在5.7～13.5℃，干湿季分明，空间上降水分配不均，年均降水量约为500～850 mm。土壤种类较丰富，有高山草甸土、暗棕壤等[10]。受气候和地形影响，植被类型多样且呈现垂直分布规律[11]，森林、灌木和草甸的面积比重分别为35.38%，40.18%和24.13%[12]。岷江上游流域是典型贫困区，经济发展落后，生产生活多以传统林业和牧业为主[5]。冷杉、云杉为该地区主要的针叶林树种。建国初期，原始针叶林遭受到严重砍伐，天然林资源几近枯竭，生态系统退化严重；直到20世纪末，随着“退耕还林”和“天然林保护工程”等生态保护项目的实施，该流域生态环境逐渐恢复[13]。

图1 岷江上游流域区位

2 数据与方法

2.1 数据来源与预处理

研究数据包括土地覆盖、土壤、降水、植被指数、地形等数据，数据信息及来源见表1。MCD12Q1土地覆被数据用于研究不同土地利用类型/植被类型土壤侵蚀变化状况；中国科学院南京土壤研究所发布的1∶100万中国土壤数据库所包含的沙粒、黏粒、粉粒、有机质含量等土壤理化信息用于土壤侵蚀模数中土壤可蚀性因子计算。基于SPOT/VEGETATION NDVI卫星遥感数据，采用最大值合成法生成2001—2017年 5—9月的空间分辨率为 1 km NDVI月值数据用于降雨侵蚀力因子计算；DEM数据空间分辨率为30 m，用于地形因子计算以及不同地形条件下土壤侵蚀状况分析。降水栅格数据采用ANUSPLIN专业插值软件对获取的降水数据(txt文件)进行插值得到。投影坐标统一为WGS_1984_UTM_Zone_47N，空间分辨率统一为500 m×500 m。

表1 研究数据介绍

2.2 研究方法

2.2.1 土壤侵蚀模数计算 使用通用土壤流失方程计算岷江上游流域逐年土壤侵蚀模数(A)，涉及降雨侵蚀力因子(R)、土壤可蚀性因子(K)、植被覆盖因子(C)、地形因子(LS)和水土保持措施因子(P)，求取土壤侵蚀模数后按照土壤侵蚀模数分级标准(表2)[14]进行分级。其中LS因子空间分布图由SAGA GIS软件中Hydrology工具生成，P因子空间分布图根据表3在ArcMap中生成，其余各因子的计算公式和相关说明如下所示[9,15]：

A=R×K×LS×C×P

fc=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoil)

式中:A为土壤侵蚀模数〔t/(km2·a)〕;R为降雨侵蚀力因子〔MJ/(km2·a)〕;K为土壤可蚀性因子(t/MJ);LS为地形因子；C为植被覆盖因子；P为水土保持措施因子,均无量纲;i为月份;Ji为月降水量；J为年降水量(mm);SD为沙粒含量；CI为黏粒含量；SI为粉粒含量；CA为有机碳含量,单位均为(%);fc为植被覆盖度；NDVIveg为全植被覆盖像元的NDVI值；NDVIsoil表示无植被覆盖像元的NDVI值；将NDVI统计直方图5%所对应的值确定为NDVIsoil；NDVI统计直方图95%所对应的值定为NDVIveg。

2.2.2 趋势、相关性以及平稳性检验 Pearson相关系数方法常用来计算线性相关、相互独立、成对观测以及符合/接近正态分布的两组变量之间直线相关性[16]。Kendall相关系数和Spearman相关系数是一种衡量两个变量依赖性的非参数指标。这两种方法对于数据的要求较少，不考虑数据的分布形态，适用性较广[17-18]。本研究采用上述方法检验土壤侵蚀模数与海拔、高程、年降水量等环境因子的时空相关性。使用偏相关方法分析不同类型植被(林地、灌木、草地)的土壤侵蚀模数和植被覆盖、降水之间的相关性，分别排除降水和植被覆盖度的影响。

单根检验(DF检验)是用来检验序列平稳性的方法，本文所使用的ADF检验是对DF检验的扩充。通过建立时间序列的自回归方程，引入滞后算子构建特征方程，根据特征根的绝对值来判断序列是否平稳，其准确性受到自回归方程精确性影响[19]。本文使用该方法检验土壤侵蚀轻度及以上面积在时间序列上的平稳性。

表2 土壤侵蚀模数分级标准

表3 岷江上游地区不同土地利用类型的P值[15]

3 结果与分析

3.1 土壤侵蚀时空变化特征

据岷江上游流域2001年、2006年、2011年、2016年土壤侵蚀分级图显示，岷江上游流域土壤侵蚀多发生在西部、西南部、东北部和东南部，汶川县、理县、黑水县和松潘县土壤侵蚀较明显(图2)。其中强度、极强度和剧烈土壤侵蚀主要发生在松潘县东部，黑水县东部，理县中部、北部及西部，汶川西部等地区。

据图3A，整个流域土壤侵蚀以微度侵蚀为主，其次是轻度侵蚀、中度侵蚀和重度侵蚀，轻度及以上侵蚀面积总体变化趋势不显著(p>0.10)且变化平稳(p<0.10)，但年际波动明显，其中2001年、2004年、2008年、2012年、2014年、2015年和2017年土壤侵蚀相对严重。2001—2017年流域土壤侵蚀变化趋势的空间分异明显。17年间土壤侵蚀模数呈显著减小的区域主要分布在流域的西部和西南部，显著增加的区域主要分布在北部和东南部(图3B)。

图2 2001-2017年岷江上游流域土壤侵蚀分级

2001—2017年土壤侵蚀模数变化显著的区域主要集中在海拔2 500～5 000 m和坡度大于15°的地区。其中49.74%的土壤侵蚀模数显著减小(p<0.10)区域位于2 500～4 000 m海拔带，51.18%的土壤侵蚀模数显著增大(p<0.10)区域位于4 000～5 000 m海拔带。89.00%的土壤侵蚀模数显著减小(p<0.10)区域位于坡度大于15°的地区，与此同时，81.80%土壤侵蚀模数显著增大(p<0.10)也分布在坡度大于15°的地区(图4)。就年均降水量看，土壤侵蚀显著减少的区域主要分布在800～1 200 mm降水带，土壤侵蚀显著增加的区域主要分布在400～800 mm降水带，面积占比分别为58.09%，53.81%。就土地利用类型看，土壤侵蚀显著减少的区域主要分布在林地，而显著增加的区域主要分布在草地。

3.2 土壤侵蚀时空分异规律

土壤侵蚀模数与海拔、坡度、年降水量、植被类型空间上显著相关(p<0.05)(表4)，土壤侵蚀的时空分异规律明显。

3.2.1 年降水量与土壤侵蚀 2001—2017年岷江上游年降水量呈上升趋势(图5A)，2001—2007年整体趋势下降，波动性较强；2008—2017年降水量增加，增速为170.33 mm/10 a。就各降水带内土壤侵蚀强度分布来看(图5B)， 400～800 mm，800～1 200 mm、大于1 200 mm降水带下的轻度及以上土壤侵蚀面积占比依次增加，说明降水量越高的地带其高强度土壤侵蚀的比例越高。

17年间，整个流域内土壤侵蚀主要发生在降水量低的区域。400～800 mm降水带土壤侵蚀面积占流域土壤侵蚀总面积的67.19%，并且以微度侵蚀为主，轻度及以上侵蚀比例仅为10.73%。800～1 200 mm降水带土壤侵蚀面积占流域土壤侵蚀总面积的32.14%，其中43.77%为轻度及以上侵蚀。>1 200 mm降水带土壤侵蚀面积最少(0.67%)，其中轻度及以上侵蚀面积占79.41%(表5)。

图3 2001-2017年不同侵蚀级别趋势及土壤侵蚀模数变化趋势空间分布

图4 土壤侵蚀模数变化趋势分布特征

3.2.2 海拔与土壤侵蚀 据图6A，不同海拔带小于1 000 m(高平原)，1 000～2 500 m(低山)，2 500～4 000 m(中山)、和大于4 000 m(高山亚高山)，中度以上侵蚀面积分别占各海拔带面积的0.44%，0.64%，1.93%，29.05%。即高海拔带高强度土壤侵蚀发生比例更高。由图6B，2001—2017年高平原和低山地带土壤侵蚀面积无明显增减，中山地带轻度及以上侵蚀面积17年内减少0.67%，相应地，高山亚高山地带增加1.05%。2001—2007年，中山、亚高山地区轻度及以上土壤侵蚀面积呈减少趋势，分别减少4.05%，0.77%。2008—2017年，中山、亚高山地带轻度及以上侵蚀面积分别增加1.66%，1.86%。说明2008年地震后，中山和高山亚高山区域土壤侵蚀面积呈持续增加趋势。

表4 土壤侵蚀模数与环境因子的空间相关性

图5 2001-2017年岷江上游流域年际降水量变化及不同降水带土壤侵蚀强度

表5 2001-2017年岷江上游流域降水带与不同级别土壤侵蚀状况

图6 2001-2017年岷江上游流域不同海拔带土壤侵蚀强度及比例变化

3.2.3 坡度与土壤侵蚀 据图7A，坡度小于5°，5°～8°，8°～15°，15°～25°，25°～35°，大于35°地区轻度及以上侵蚀面积占各坡度带比例分别为4.22%，7.76%，14.10%，19.38%，22.38%，26.45%，其中轻度及以上侵蚀集中发生在坡度大于15°地区，并随着坡度增大高强度土壤侵蚀占比逐渐加大。

据图7B，2001—2017年，坡度小于15°地带土壤侵蚀变化不明显，坡度大于15°地带土壤侵蚀变化趋势明显，其中15°～25°，25°～35°坡度带轻度及以上侵蚀面积分别增加0.40%，0.10%，大于35°坡度带轻度及以上侵蚀面积减少0.75%。2001—2007年15°～25°，25°～35°，大于35°坡度带土壤侵蚀状况呈转好趋势，轻度及以上侵蚀面积分别减少0.75%，1.63%，3.04%。2008—2017年15°～25°，25°～35°，大于35°坡度带土壤侵蚀加剧，各坡度带轻度及以上侵蚀面积分别增加0.86%，0.90%，0.05%。

图7 2001-2017年岷江上游流域不同坡度土壤侵蚀强度及比例变化

3.2.4 土地利用与土壤侵蚀 据图8A，岷江上游流域土壤侵蚀主要发生在未利用地、草地、林地，而耕地、建设用地、水域湿地面积在该流域占比较少，且以微度侵蚀为主(99%以上)，其中耕地土壤侵蚀区域仅有0.68%发生轻度及以上退化。就每种土地利用类型内发生较强土壤侵蚀的情况看，轻度及以上侵蚀程度未利用地(87.11%)>草地(39.75%)>林地(9.49%)。由图8B看出，17年间，草地轻度及以上土壤侵蚀呈上升趋势，2017年较2001年增加1.87%；2001—2007年草地土壤侵蚀有明显改善，轻度及以上侵蚀面积下降0.55%，但2008—2017年侵蚀面积明显上升，增加2.25%。2001—2017年林地、未利用地土壤侵蚀面积无明显增减。

由图9看出，不同林地类型发生轻度及以上土壤侵蚀面积占比由大到小分别为疏林地(15.04%)>针叶林(13.50%)>混交林(4.41%)>阔叶林(0.97%)，说明疏林地易发生土壤侵蚀，针叶林、混交林次之，而阔叶林土壤侵蚀程度相对较低。

图8 2001-2017年岷江上游流域不同土地利用类型土壤侵蚀强度及面积比例变化

图9 不同林地类型土壤侵蚀强度

3.3 土壤侵蚀对长期植被变化响应

由图10看出, 2001—2017年岷江上游流域林地、灌木和草地的NDVI总体均呈现显著增加的趋势(p<0.05)，且变化过程大致相同。2001—2007年3种植被类型NDVI呈上升趋势，受2008年汶川地震影响，2008—2010年NDVI有明显下降，2010年后植被逐渐恢复，3种植被类型NDVI呈不同程度的上升趋势。

分别对各植被类型的土壤侵蚀模数A与年降水量PRE，NDVI进行偏相关分析，结果见表6。林地土壤侵蚀模数A与NDVI呈显著负相关(r=-0.66,p<0.05)；林地、灌木和草地土壤侵蚀模数A与年降水量PRE呈显著正相关(r=0.73,p<0.05;r=0.64,p<0.05;r=0.95,p<0.05;r=0.92,p<0.05)。说明林地土壤侵蚀主要受到植被覆盖和降水的影响，而灌木和草地土壤侵蚀则主要受到降水的影响。

图10 不同植被类型NDVI变化趋势

表6 土壤侵蚀模数与年降水量和NDVI偏相关分析

4 讨 论

4.1 土壤侵蚀空间分布规律及成因

研究结果表明岷江上游流域土壤侵蚀模数多年平均值为801.60 t/(km2·a)，与四川省水土重点防治区土壤侵蚀模数多年均值接近[20]；并且本研究2010年、2017年流域平均土壤侵蚀模数计算结果与刘佳[21]研究结果也基本一致。岷江上游流域土壤侵蚀以微度和轻度侵蚀为主，其次是轻度侵蚀、中度侵蚀和重度侵蚀。总体而言土壤侵蚀呈恶化趋势，总体结果与相似区域研究基本相近[9,22-25]。尽管使用数据的来源、空间分辨率以及计算方法的不同会造成相同区域同一时间段内识别结果的差异[2,9,21]，本研究识别出的土壤侵蚀分级空间分布与姜琳、邓兵、祝聪等[2,15,26]得出的总体规律基本一致。该流域断裂发育较多，地质构造活动较活跃，地表破碎程度较大，并且分布有大面积典型干旱河谷生态脆弱区[27]。受地形地貌、降水、植被覆盖、地质活动以及人类活动的影响，土壤侵蚀分布的时空异质性较高。

岷江上游流域土壤侵蚀主要分布在年降水量低的低海拔干旱河谷区。该流域约20.8%的区域为干旱河谷(1 300～2 200 m)，植被类型以干旱河谷灌丛为主，植被覆盖度低，裸露地表比例较大，雨季极易发生大面积土壤侵蚀[15,21]。由于总体降水量低，该区域土壤侵蚀主要以微度和轻度侵蚀为主。而高海拔区域，降水量高，植被类型以阔叶林和暗针叶林为主，植被覆盖度较高，林冠层和枯落物层对降水具有较强的截留作用，能够有效地减少土壤侵蚀，土壤侵蚀分布面积较小。但是在高海拔、高降水区域，土壤侵蚀以高强度土壤侵蚀为主，与王猛等[28]研究结果一致。水土流失通常与侵蚀雨量呈极显著的线性关系[29]，降水通过影响径流大小从而影响土壤侵蚀程度[6]。由于该流域高海拔地区山高坡陡，集中高强度降水(集中发生在6—10月，多为大到暴雨)快速汇集产生径流，极易发生泥石流、滑坡。地表若缺乏有效植被覆盖，极易引发高强度土壤侵蚀[30-31]。

就地形而言，随着海拔和坡度上升，土壤侵蚀程度越大，该结果与何兴元[9]、邱国玉[32]等的研究结果一致。岷江上游流域植被分布随海拔变化呈现垂直差异，海拔由低到高主要分布的植被为干旱森林和灌丛(1 300～2 600 m)、针阔混交林和针叶林(1 800～3 900 m)、高山灌丛(3 600～4 300 m)和高山草甸(>4 300 m)[33]。据研究表明，林地、草地、灌木土壤保持能力/强度由大到小分别为林地>灌木>草地[34-36]。高海拔4 000～5 000 m地区由于多高山裸土以及荒草地，土壤保持能力差，加之海拔3 800～4 200 m地区大部分地层属于第四系，其结构松散，颗粒较细，在地质运动和集中强降雨冲刷下高海拔地区更易发生侵蚀[2,31]。与此同时，岷江上游流域约50%区域坡度大于45°，地势陡峭，河谷深切。由于土壤入渗量随坡度增大而减小，地表产流量随之增大，加之重力作用随坡度增加不断增强，水力侵蚀效果越明显，进而导致坡度大的区域土壤侵蚀量越大[6,15,37]。

土地利用类型的差异也是造成流域土壤侵蚀空间分异明显的重要因素之一。岷江上游流域发生较强土壤侵蚀的3种土地利用类型主要为未利用地、草地和林地。未利用地多为裸地，其植被稀疏、极易受到降水的侵蚀。虽然林地、草地均具有固结土壤表层，滞留水分，减少水土流失的能力[8]，但是森林植被覆盖度相对较高，且根系较草地更发达粗壮，使得森林生态系统具有更高的稳定性[38]，林地土壤保持能力更强。相关研究也表明，林地、草地、未利用地土壤保持能力由大到小依次为林地>草地>未利用地(裸地)[35-36,39]。因此，岷江上游流域各土地利用类型发生土壤侵蚀程度由大到小依次为未利用地>草地>林地，基本与多数研究结果一致[2,21,40]。森林植被枯落物层能延缓坡面径流产生，降低坡面径流流速，在一定程度上减少坡面产流产沙[41]。岷江上游流域不同林分枯落物储量大致表现为阔叶林>针阔混交林>针叶林[42-45]，枯落物储量越大，在一定程度上减少地表径流能力越强。岷江上游流域针叶林主要分布在海拔3 500～4 500 m，阔叶林、混交林分布在海拔2 500～3 500 m[10]，呈现明显地带差异，高海拔地区更易发生土壤侵蚀。此外，由于不同林地类型下表层土壤抗蚀性表现为阔叶林>针阔混交林>针叶林[46-49]，疏林地水土保持能力很低，同时由于林冠层稀疏，强降雨下容易发生溅蚀，因此研究区不同林型土壤侵蚀程度差异表现为疏林地>针叶林>混交林>阔叶林。

4.2 土壤侵蚀长期动态及成因

2001—2017年岷江上游流域土壤侵蚀呈缓慢增加趋势，与年降水量的变化大致相同。因此降水是影响不同植被覆盖区域土壤侵蚀的主要因子[50]。但不同时段影响因子有所差异。2001—2007年，各个海拔带、坡度带和土地利用类型的土壤侵蚀变化基本与降水年际波动一致。尽管该时段林地、灌丛、草地植被指数均呈增加趋势，但并未对土壤侵蚀起到明显的减缓作用[2]。2008—2017年该地区大暴雨极端天气频发，降水量呈显著增加的趋势，同时岷江上游流域地处龙门山断裂带，地质运动较活跃，尤其是受2008年汶川地震影响，大面积植被遭受破坏，同时在流域高海拔区域孕育了相当丰富的固体松散物和堆积物，叠加陡峻的地势条件，为岷江上游流域次生地质灾害(泥石流、滑坡等)频繁发生提供了有利条件。受地质灾害和降水等因素的共同作用，流域土壤侵蚀呈逐年上升趋势[51]，但趋势性弱于降水量的变化。其中以高海拔带(>2 500 m)、坡度大于25°以及草地为主的地区土壤侵蚀增加趋势最为明显。值得注意的是，2008—2017年林地和灌丛NDVI呈明显增加趋势，而草地NDVI增加速度弱于上述两种土地利用类型。由此可见，天然林保护工程、退耕还林以及震后植被修复工程，显著提高了林地和灌丛的植被覆盖度[52]，对林地和灌丛区土壤侵蚀的控制起到了较好的效果，有效减缓了降水和地质灾害对土壤侵蚀的增加作用。由于草地大量向林地转移[21]，且草地结构简单易受干扰，因而草地植被恢复力度较林地和灌丛弱，其土壤侵蚀防止作用提升效果则不明显。

综合而言，岷江上游流域一系列植被恢复工程对土壤侵蚀控制起到了一定积极的作用，但并未从总体上遏制土壤侵蚀增加的趋势。因此，未来该区域生态恢复的重点为进一步研发林地、灌丛和草地水源涵养功能提升的恢复措施。尤其是需要加强对高海拔带(>2 500 m)、坡度大于25°以及草地区域的植被恢复研究，进一步提高各类型植被水源涵养功能，减少区域土壤侵蚀。

不过本研究针对性分析单因子与土壤侵蚀之间关系，对于不同因子之间的交互作用的考虑尚有欠缺。根据相关人员[53-55]在渭河流域、子午岭地区、天山北坡中段山区和杨凌示范区的研究表明植被覆盖度、造林面积、降雨、海拔、坡度、土壤类型、人类活动等影响因子变化对土壤保持变化均起交互增强作用。王猛等对西南地区土壤侵蚀影响因子进行交互作用分析也得出两两因子之间交互协同作用表现为非线性增强或相互增强，其中坡度和盖度交互作用最大[28]。因此，在今后的研究中使用更高分辨率的数据并且可加强因子之间交互作用的研究，提升研究结果的完整性和精确度。

5 结 论

岷江上游流域轻度及以上土壤侵蚀区域主要分布在理县、汶川县、黑水县和松潘县。土壤侵蚀模数与海拔、坡度、年降水量和植被类型在空间上显著相关，表现为随海拔(坡度或年降水量)增加，土壤侵蚀程度增加。林地、草地、未利用地为岷江上游发生土壤侵蚀的主要土地利用类型，轻度及以上侵蚀程度表现为未利用地>草地>林地，不同林地类型土壤侵蚀程度表现为疏林地>针叶林>混交林>阔叶林。2001—2007年，各个海拔带、坡度带和土地利用类型的土壤侵蚀变化基本与降水年际波动一致。2008—2017年受地质灾害和降水等因素的共同作用，流域土壤侵蚀呈逐年上升趋势，其中以高海拔带(>2 500 m)、坡度大于25°以及草地为主的地区土壤侵蚀增加趋势最为明显。天然林保护工程、退耕还林以及震后植被修复工程，显著提高了林地和灌丛的植被覆盖度，有效减缓了降水和地质灾害对土壤侵蚀的增加作用。而草地植被恢复力度较林地和灌丛弱，其土壤侵蚀防止作用提升效果则不明显。