黄河源区高寒退化草地空间分布特征

李积兰 王苑 李希来 刘慧 王晓波

摘 要 为了解黄河源地区“黑土滩”退化草地的空间分布规律和分布的关键地形特征，利用2019年黄河源区玛沁县、河南县、泽库县遥感影像解译数据，结合实地调查，分析高寒“黑土滩”型退化草地空间分布特征。结果表明：黄河源区“黑土滩”型退化草地海拔为3 500～4 000 m，且在海拔3 700～3 800 m时退化最严重；从坡向分析来看，“黑土滩”型退化草地主要发生于阳坡、半阳坡，从坡度分析来看，“黑土滩”型退化草地主要发生于缓坡地和滩地。聚类分析发现研究区草地类型为阳坡极度退化草地、半阳坡重度退化草地、阳坡-滩地轻度退化草地、半阳坡轻度退化草地、湿地轻度退化草地、阴坡未退化草地等6类。随着退化程度的加剧，植被盖度、地上生物量减少；同一土层的全氮、全磷、碱解氮、速效磷、有机质含量在湿地轻度退化草地和阴坡轻度退化草地中最高，阳坡极度退化草地最低，且差异显著（P<0.05）；丰富度指数、多样性指数、均匀度指数以及土壤速效钾含量在重度、轻度退化草地较高，湿地轻度退化草地显著最低。阳坡极度退化草地全钾含量、鼠害破坏率和鼠洞数最高，湿地轻度退化草地全钾含量最低，与阳坡极度退化草地差异显著（P<0.05）；湿地轻度退化草地和阴坡轻度退化草地均无鼠类活动。

关键词 黄河源区；退化草地；空间分布；植物群落

黄河源区位于青藏高原腹地，三江源自然保护区的核心区域，是黄河流域的主要产流区、水源涵养区，也是中国重要的生态屏障［1］。近几十年来，黄河源的局地气候发生了变化，加之超载过牧、乱砍滥垦、随意采挖、车辆行驶等人为因素的干扰，使得源区生态环境恶化，出现草地持续退化［2-3］。最常见的问题为草毡层破坏、草地破碎化、鼠害、冻融、侵蚀，毒杂草蔓延，秃斑块产生和扩大，最终导致“黑土滩”的形成［4-6］。秃斑化是“黑土滩”最典型的表现，且在高寒草地退化中呈中、重度以上，也是鼠洞密度最高的区域［7-9］。

近年来，诸多学者进行了“黑土滩”退化草地的成因分析和治理、恢复措施的研究，发现  2000-2008年三江源区草地覆盖度降低［10］，  2008-2009年的遥感信息判断玛多县的“黑土滩”面积增加至  1 592 km2［11］，也有研究认为1994-2009年草地退化面积逐渐增加，2005-2011年间黄河源区土壤水蚀增加，退化草地尚未得到恢复［12-13］，2007年以前“黑土滩”扩张态势，2011年以后由于退牧还草和生态保护工程减轻了草地退化趋势，使得退化速度降低［14-15］。黄河源区1996-2015年植被覆盖度呈增加趋势的区域占57.25%，基本不变的区域占16.02%，植被覆盖度呈下降趋势的区域占26.73%。植被覆盖度下降的主要原因是黄河源头及一些河谷地带、环湖地区受人类影响较大，且东南部海拔较低地区受到过度放牧的影响，而且毒杂草的面积也自1996年至2015年，  20 a间增加了6 882 km2［16］，表明有些地区“黑土滩”退化草地面积仍处于扩张态势。

许多学者在退化草地时空分布、成因分析、恢复治理等方面做了大量研究工作，但遥感影像解译与实地验证相结合进行退化草地空间分布的分析研究报道较少。为进一步了解黄河源地区“黑土滩”退化草地的空间分布规律和分布的关键地形特征，本研究利用LandSat卫星影像，采用目视解译结合实地调查的方式，对黄河源地区—玛沁县、泽库县、河南县2019年遥感影像数据和高寒草地退化草地群落结构、土壤和啮齿动物活动特征进行分析，以期阐释研究区草地退化的空间分布规律与退化原因，为黄河源地区草地退化的综合治理与区域可持续发展提供数据支持。

1 材料与方法

1.1 研究地自然概况

研究地区位于青藏高原黄河源区内，主要包括玛沁县、河南县、泽库县。玛沁县位于东经  98°～100°56′，北纬33°43'～35°16′。河南县位于东经100°53′～102°12′，北纬34°13′～34°52′。泽库县地理坐标在东经100°34′～102°08′，北纬  34°45′～35°32′。气候特点为高原大陆性季风气候，属高原亚寒带湿润气候区。年均降雨量  535 mm，雨热同季，降水量较丰富且集中；年均气温  1 ℃，气温较低，热量不足，无霜期短，日照时间长，太阳辐射强;冷季漫长干冷，多大风；暖季短促凉爽。土壤类型为高寒草地土，草地类型为高寒草地，主要建群种为高山嵩草（Kobresia pygmaea）、矮嵩草（Kobresia humilis）、线叶嵩草（Kobresia capillifolia）、藏嵩草（Kobresia tibetica）、垂穗披碱草（Elymus nutans）、紫花针茅（Stipa purpurea）等［17-18］。

1.2 样地选择

采样时间为8月初，样区选择在玛沁县、河南县、泽库县境内，与影像解译数据样区一致，每个样区面积400  km2（图1），每个样区内设置5个样地，样地尽可能选在具有代表性的区域且均匀分布，面积20 m×20 m，在每个样地的对角线上均匀设置3个1 m×1 m的样方，共有15个样地，45个样方，样地位置图见图1。

1.3 观测指标与测定方法

1.3.1 遥感影像解译 人工目视解译广泛应用于植被类型的解译中［19］。采用LandSat8（OLI）卫星影像，对2019年8月玛沁县、河南县、泽库县全域影像数据进行人工目视解译。影像预处理过程：对影像进行几何纠正，然后影像融合，形成  15 m分辨率的工作影像。在遥感解译中，“黑土滩”退化草地为植被覆盖度小于60％以下的草地［20-21］。

1.3.2 野外验证调查 利用青海省地理国情监测外业调绘核查系统进行野外验证调查，地面样方调查采用高精度手持GPS进行定位，可精确定位每一个样点，对应到15 m遥感影像上的空间位置。为了更加客观反映“黑土滩”的基本情况，在选取研究地时，根据原生植被和秃斑地占研究地的比例，划分为原生、轻度退化、重度退化和极度退化四种程度［22-23］。根据玛沁县、河南县、泽库县影像解译情况，利用青海省地理国情监测外业调绘核查系统分别对3个县的样区进行实地调查（图2）。玛沁县样区海拔起伏较大，到达的样地基本为海拔较低、坡度较小的区域，样区内阳坡、半阳坡、缓坡地和滩地的草地退化较严重，草地植被类型主要以黄帚橐吾（Ligularia virgaurea）、白苞筋骨草（Ajuga lupulina）、圆萼刺参（Morina chinensis）、铁棒锤（Aconitum pendulum）等为主；河南县样区海拔起伏变化不大，样区内阳坡、半阳坡、缓坡地和滩地的草地退化较严重，草地植被类型主要以密花香薷（Elsholtzia densa Benth.）、黄帚橐吾（Ligularia virgaurea）、鵝绒委陵菜（Potentilla anserina）、白苞筋骨草（Ajuga lupulina）等为主；泽库县样区实际地理位置包含了同德县，由于样区内大部分区域为同德牧场，退化严重的区域为湿地、缓坡地和滩地，植物以狼毒大戟（Euphorbia fischeriana）、兰石草（Lancea tibetica）、冰草（Agropyron cristatum）、白苞筋骨草（Ajuga lupulina）、细叶亚菊（Ajania tenuifolia）、鹅绒委陵菜（Potentilla anserina）等为主。

1.3.3 土壤理化性质的测定 在同一样方内，按照对角线用土钻分别采集0～10 cm、10～20 cm、20～30 cm的土样，将同一样方内同一土层的土样混匀后编号装入自封袋内，带回实验室自然风干，进行土壤养分的测定。土壤全氮、全磷、速效磷含量测定用SMART CHEM 450全自动间断化学分析仪测定；土壤碱解氮含量测定用水解-扩散法；土壤全钾、速效钾含量用火焰光度计法测定［24］；pH用土壤原位pH计（IQ150）测定。

1.3.4 植物群落的测定 调查并记录样方内的海拔、经纬度及每个样方的植物种类组成、植物高度、盖度、生物量（地上）等［25］。

盖度：目测法测定样方总盖度和各物种分盖度。

相对盖度=样方内某种植物的分盖度/样方所有物种分盖度之和×100%

高度：样方中每种植物选取10株用钢卷尺测定其自然高度（不足10株的按实际株数测定）。

相对高度=样方内某种植物的平均高度/样方内所有物种的平均高度之和×100%

鲜质量：0.5 m×0.5 m的样方中按物种分别用剪刀齐地面剪下，用电子天平（千分之一）测其鲜质量。

相对生物量=样方内某种植物的质量/样方内所有物种的质量之和×100%

重要值=（相对高度+相对盖度+相对生物量）/3

群落多样性测定：物种的多样性指数选用以下4个指标表示：Patrick丰富度指数、Pielou均匀度指数、Shannon-Wiener 多样性指数和SimPson多样性指数［26］。

1.3.5 鼠害调查 在每个20 m×20 m的样地内，采用对角线取样，分别测定洞口、秃斑地和秃丘的长度。破坏率=各项所截长度的总和/样线长×100%［27］。

1.3.6 数据处理 采用ArcGIS软件进行矢量数据采集；Excel 2019对数据进行整理；用  Spss.20软件LSD法进行多重比较，结果以“平均数±标准差”表示，以P<0.05为差异显著的判断标准；用Hierarchical Cluster法进行聚类分析；绘图用SigmaPlot 10。

2 结果与分析

2.1 黄河源区“黑土滩”型退化草地空间分布

2.1.1 “黑土滩”型退化草地面积 根据2019年玛沁、河南、泽库3 个县解译的遥感影像图及数据（图3-a，3-b），玛沁县“黑土滩”型退化草地面积最大，19.87 km2，占所有草地面积的6.66%；泽库县“黑土滩”型退化草地面积次之，3.62 km2，占所有草地面积的1.11%，河南县“黑土滩”型退化草地面积最小，2.99 km2，占所有草地面积的  0.80%，表明黄河源区玛多县草地退化最严重，其次为泽库县，河南县较轻。

2.1.2 高寒退化草地海拔梯度分布 玛沁县退化草地主要分布于海拔3 600～4 100 m和  4 200～4 800 m，退化面积为1 986.97 hm2，其中在海拔3 700～4 100 m退化较严重，退化面积占所有海拔区域的97.90%（图4-a，4-d）。河南县退化草地分布于海拔3 500～3 800 m和4 401～4 500 m，退化面积为3.59 km2，其中3 500～  3 700 m退化比较严重，退化面积占所有海拔区域的72.24%（图4-b，4-d）。泽库县退化草地分布于海拔3 500～3 800 m，退化面积为3.63 km2，其中在海拔3 600～3 800 m草地退化比较严重，退化面积占所有海拔区域的88.76%（图4-c，4-d）。因此，黄河源区“黑土滩”型退化草地的海拔主要在3 501～4 000 m。

2.1.3 高寒退化草地坡度分布 由图5可以看出，3个县的“黑土滩”型退化草地面积均在缓坡地最大，其中玛沁县“黑土滩”型退化草地在缓坡地占50%，滩地占45%，陡坡地占4%（图5-a，5-b）；河南县“黑土滩”型退化草地在缓坡地占72%，陡坡地占15%，滩地占12%（图5-c，5-d）；泽库县“黑土滩”型退化草地在缓坡地占71%，陡坡地占23%，滩地占5%（图5-e，5-f）。因此，黄河源区“黑土滩”型退化草地主要分布于缓坡地、滩地。

2.1.4 高寒退化草地坡向分布 由图6可知，3 个县的“黑土滩”型退化草地面积均由阳坡、半阳坡、半阴坡、阴坡逐渐减小，其中玛沁县“黑土滩”型退化草地面积较大，且阳坡和半阳坡面积最大，占所有坡向面积的52.89%，阴坡和半阴坡稍低，占所有坡向面积的47.12%（图6-a，6-b）。河南县阳坡地“黑土滩”型退化草地面积最大，占45%，半阳坡占27%，半阴坡占16%；阴坡占12%（图6-c，  6-d）。泽库县“黑土滩”型退化草地坡向分布趋势与河南县一致，阳坡地面积占58%，半阳坡地占29%，半阴坡地占28%，阴坡地占5%（图6-e，6-f）。因此，“黑土滩”型退化草地主要分布于阳坡地、半阳地。

2.2 高寒草地实地调查群落结构基本特征

根據实地调查及植物、土壤相关指标的测定结果，对黄河源区玛沁县、河南县、泽库县3个县15个样地45个样方的31个植物群落、土壤相关指标进行聚类分析（图6）。经综合考虑，标尺值选择5，草地退化类型分为阳坡极度退化草地、半阳坡重度退化草地、阳坡-滩地轻度退化草地、半阳坡轻度退化草地、湿地轻度退化草地、阴坡未退化草地6类（表1）。根据聚类分析结果，极度退化草地分布于阳坡，且主要在玛沁县，这与遥感影像解译数据显示的阳坡地退化较严重和玛沁县退化草地面积较大的结果一致，说明解译数据是可靠的。

2.2.1 高寒草地植被特征 根据聚类分析结果对不同退化类型草地的植物群落进行分析（表2）。植被盖度随着草地退化程度的加剧而逐渐降低，阳坡极度退化草地与湿地轻度退化草地、阴坡未退化草地间差异显著（P<0.05）。丰富度随着草地退化程度的增加呈现先增大后减小的趋势，湿地轻度退化草地与半阳坡轻度退化草地差异显著（P<0.05），其他地形差异不显著（P>0.05）。多样性指数Shannon-wiener、Simpson随着草地退化程度的加剧，变化不明显，湿地轻度退化草地与半阳坡轻度退化草地差异显著（P<0.05）。均匀度指数Pielou随着草地退化程度的加剧，各地形差异不显著（P>0.05）。地上生物量随着退化程度的加剧而减小，各地形差异不显著（P>  0.05）。因此，阳坡极度退化草地的植被盖度、地上生物量较其他草地类型低，且主要发生在海拔3 800 m左右；湿地轻度退化草地丰富度指数、多样性指数和均匀度指数均较其他草地类型低，但地上生物量较其他地形高。

2.2.2 土壤理化性质特征 由表3可知，在不同退化类型的同一土层中，全氮、全磷、碱解氮、有机质含量变化趋势基本一致，均在湿地轻度退化草地、阴坡未退化草地含量最高，阳坡极度退化草地最低，阴坡未退化草地与湿地轻度退化草地差异不显著（P>0.05），与其他退化草地类型差异显著（P<0.05）。而全钾、速效钾含量变化趋势与其他指标相反，均在阳坡极度退化草地、半阳坡重度退化草地、半阳坡轻度退化草地、阳坡-滩地轻度退化草地含量较高，湿地轻度退化草地、阴坡未退化草地含量较低，且湿地退化草地的全钾含量与其他草地类型差异显著（P<0.05），其余草地类型间差异不显著（P>0.05）；阴坡退化草地速效钾含量与阳坡极度退化草地差异不显著（P>  0.05），与其他类型间差异显著（P<0.05）。速效磷含量在半阳坡轻度退化草地含量最高，阳坡-滩地轻度退化草地次之，二者差异不显著（P>  0.05），与其他草地类型差异显著（P<0.05）；半阳坡重度退化草地含量最低，且与其他草地类型差异显著（P<0.05），其余草地类型间差异不显著（P>0.05）。

同一退化类型不同土层中，随着土层的加深，全氮、全磷、全钾、碱解氮、速效磷、速效钾、有机质含量逐渐降低，全氮含量除了在阳坡极度退化草地的各土层间差异不显著外（P>0.05），在其余草地类型的各土层间差异显著（P<0.05）；全磷含量在所有草地类型的各土层间差异均不显著（P>0.05）。碱解氮含量0～10 cm土层与20～30 cm土层差异显著（P<0.05）。在所有退化类型中，0～10 cm土层的速效磷、速效钾、有机质含量与其他土层差异显著（P<0.05）。

2.2.3 鼠害情况特征 根据表4，鼠害破坏率、洞口数随着草地退化程度的减弱而减弱，阳坡极度退化草地的鼠洞数量和鼠害破坏率最高，半阳坡重度退化草地次之，半阳坡轻度退化草地与阳坡-滩地轻度退化草地中较少，湿地轻度退化草地和阴坡未退化草地中均未有鼠害破坏情况发生。鼠洞数量、鼠害破坏率在极度退化草地与重度退化草地和轻度退化草地差异显著（P<0.05），在半阳坡重度退化草地与阳坡-滩地轻度退化草地差异不显著（P>0.05），与半阳坡轻度退化草地差异显著（P<0.05），半阳坡轻度退化草地与阳坡-滩地轻度退化草地差异不显著（P>0.05）。

3 讨  论

3.1 黄河源区“黑土滩”型退化草地的空间分布特征

“黑土滩”指退化的高寒草地，是青藏高原独特的自然条件下高寒草地严重退化经风蚀和水蚀后形成大面积裸地产物，只有在特定生态区域的高寒草地才能退化成“黑土滩”［28］。从海拔分布来看，本研究认为黄河源区“黑土滩”型退化草地主要分布于海拔3 500～4 000 m，在海拔3 700～  3 800 m时退化最严重。这与已有研究认为“黑土滩”的海拔分布在3 000～5 000 m［29］和   3 600～4 500 m［30］的范围是一致的，这可能与随着海拔的升高，家畜和人类活动减少有关［31］。阳坡地的退化面积最大，其次为半阳坡、半阴坡，阴坡地退化面积最小，缓坡地的“黑土滩”型草地面积大于滩地和陡坡地；因此从坡向分析来看，“黑土滩”型退化草地主要发生于阳坡、半阳坡，从坡度分析来看，“黑土滩”型退化草地主要发生于缓坡地和滩地，这与已有研究发现的“黑土滩”多出现于青藏高原山地阳坡和半阳坡山麓及山前滩地一致［14］，这可能是由于家畜和人类更喜欢在阳坡、半阳坡、缓坡地和滩地活动，因此“黑土滩”的发生可能与该地区特殊的生态系统结构和生态过程有关。

3.2 高寒退化草地植物群落结构特征

经实地调查及聚类分析，草地类型划分为阳坡极度退化草地、半阳坡重度退化草地、阳坡-滩地轻度退化草地、半阳坡轻度退化草地、湿地轻度退化草地、阴坡未退化草地6类。草地植被变化是草地退化最明显的表现，研究区植被盖度和地上生物量随着退化程度的加剧而减少，且阳坡极度退化草地最低，湿地轻度退化草地最高。物种丰富度指数是反映群落中物种多少的参数，本研究中丰富度指数在轻度退化草地较高，而湿地较低，表明轻度退化草地由于阔叶类植物的增加［32］，使得植物较丰富，而湿地主要以藏嵩草和苔草为主，植物种类较单一，丰富度指数较低。多样性是反映一個群落内物种的数量和均匀程度的测度指标［33］，生态系统的退化，不仅使群落在物种组成上发生重要变化，而且也导致物种多样性发生了改变［34］。在本研究不同退化草地类型中，多样性指数和均匀度指数变化不同，其中湿地的多样性指数较低，轻度退化草地较高，与许世洋等［35］研究得出中轻度退化草地多样性指数较高的结果一致。

3.3 高寒退化草地土壤理化性质特征

随着退化程度的加剧，土壤全氮、全磷、碱解氮、速效磷、有机质含量呈下降趋势，与有关研究［36-38］结果相吻合。这是由于土壤有机质是土壤肥力的物质基础，土壤全氮的95%来源于有机质，因而有机质在土壤中的积累会直接影响氮、磷含量的变化。另外，当植被发生退化时，土壤中的固氮作用减弱甚至停止，并且退化加速了生物残体的分解，有机氮转化为无机氮，进而在水土流失中损失［39］。全钾、速效钾含量在重极度退化草地含量较高，且0～10 cm土层的速效磷、速效钾、碱解氮、有机质含量相较其他土层差异显著（P<  0.05）。原因可能是磷和钾含量受土壤风化的影响，草地退化加速了土壤的风化，使得土壤中的磷和钾释放更多，因此全磷和全钾的含量有所升高，而退化草地中土壤微生物减少，使得植物残体的分解速率降低，导致速效磷降低［40］，湿地的速效钾含量低的原因可能是湿地含水量较高，不利于速效钾的分解。

3.4 高寒退化草地鼠害活动特征

草地秃斑、裸露和破碎化过程是从阳坡-滩地开始的［41］，且在鼠类活动区域首先形成“黑土滩”［42］，中、重度“黑土滩”退化草地是鼠洞密度非常高的地区［7，43］，“黑土滩”分布区的关键地形特征与高原啮齿动物栖息地条件、放牧家畜的适宜放牧区相吻合。本研究发现阳坡极度退化草地的破坏率最大、鼠洞数最多，说明鼠类喜活动于阳坡、滩地；湿地轻度退化草地、阴坡未退化草地无鼠类活动。也有研究认为鼠兔相对密度与土壤湿度之间呈显著的负相关，河漫滩湿地地势平坦，土壤湿度大，鼠洞密度最小［44］，与本研究结果一致。

4 结论

4.1 黄河源区“黑土滩”型退化草地的空间分布

黃河源区“黑土滩”型退化草地主要发生在海拔3 500～4 000 m的阳坡、半阳坡的缓坡地和滩地，以海拔3 800 m左右退化最为严重。

4.2 高寒退化草地群落结构特征

退化草地类型划分为阳坡极度退化草地、半阳坡重度退化草地、阳坡-滩地轻度退化草地、半阳坡轻度退化草地、湿地轻度退化草地、阴坡未退化草地6类。植被盖度、地上生物量随着退化程度的加剧而减少，丰富度指数、多样性指数、均匀度指数均在轻度退化草地较高，在湿地轻度退化草地较低。

4.3 高寒退化草地土壤理化性质基本特征

同一土层不同退化类型草地的全氮、全磷、碱解氮、速效磷、有机质含量总体随着退化程度的加剧而降低，其中湿地轻度退化草地、阴坡轻度退化草地中的全氮、全磷、碱解氮、有机质含量高于其他退化草地；阳坡重度退化草地速效磷含量最低；全钾含量在阳坡极度退化草地最高，在湿地轻度退化草地最低；速效钾含量在半阳坡重度退化草地、半阳坡轻度退化草地、阳坡-滩地轻度退化草地和阴坡轻度退化草地较高，在阳坡极度退化草地和湿地轻度退化草地较低。同一退化类型，随着土层深度的加深，各指标总体呈下降趋势，0～10 cm土层的速效磷、速效钾、碱解氮、有机质含量均显著高于其他土层。

4.4 高寒退化草地鼠害活动特征

研究区鼠害破坏率、鼠洞数随着草地退化程度的加剧而增加。阳坡极度退化草地的破坏率最大、鼠洞数最多，重度退化草地次之，半阳坡轻度较低；湿地轻度退化草地、阴坡未退化草地无鼠类活动。

参考文献 Reference：

［1］ 赵志平，吴晓莆，李 果，等.黄河源区高寒草地 NDVI 格局与梯度变化［J］.草业科学，2013，30（12）：1917-1925.

ZHAO ZH P，WU X P，LI G，et al. Spatial pattern and change gradient of alpine grassland in the Source Region of Yellow River［J］. Patacultural Science，2013，30（12）：1917-1925.

［2］ 王 莺，李耀辉，孙旭映.气候变化对黄河源区生态环境的影响［J］.草业科学，2015，32（5）：539-551.

WANG Y，LI Y H，SUN X Y. Impact of climate change on the eco-environment in the Yellow River Source［J］. Patacultural Science，2015，32（5）：539-551.

［3］ 翟颖佳，李耀辉，陈玉华.全球及中国区域气候变化预估研究主要进展简述［J］.干旱气象，2013，31（4）：803-813.

ZHAI Y J，LI Y H，CHEN Y H. Major progress of global and China regional climate change projection［J］.Journal of Arid Meteorology，2013，31（4）：803-813.

［4］ DONG Q M ，ZHAO X Q ，WU G L，et al. A review of formation mechanism and restoration measures of “black-soil-type” degraded grassland in the Qinghai-Tibetan Plateau［J］.Environmental Geology，2013，70（5）：2359-2370.

［5］ 李子好.青藏高原三江源区“黑土滩”退化高寒草地秃斑块分布格局及土壤养分空间变异［D］.兰州：兰州大学，2020.

LI Z H. Spatial pattern and variation of bald patches and soil nutrition in ‘bare land degraded apine grassland in Three-river source region on the Tibetan Plateau［D］. Lanzhou：Lanzhou University，2020.

［6］ LIN L，LI Y K，XU X L，et al. Predicting parameters of degradation succession processes of Tibetan Kobresia  grasslands［J］. Solid Earth，2015，6（4）：1237-1246.

［7］ 韩立辉，尚占环，任国华，等.青藏高原“黑土滩”退化草地植物和土壤对秃斑块面积变化的响应［J］.草业学报，2011，20（1）：1-6.

HAN L H，SHANG ZH H，REN G H，et al. The response of plants and soil on black soil patch of the Qinghai-Tibetan Plateau to variation of bare-patch areas［J］. Acta Prtacuitur Sinica，2011，  20（1）：1-6.

［8］ 林慧龙，王钊齐，尚占环.江河源区“黑土滩”退化草地秃斑与鼠洞的分形特征［J］.草地学报，2010，18（4）：477-484.

LIN H L，WANG ZH Q，SHANG ZH H. Features on fractal dimension of barren patch and mouse hole among different degenerated succession stages on alpine meadow in the source region of the Yangtze and Yellow river，Qinghai-Tibetan Plateau，China［J］. Acta Agrestla Sinica，2010，  18（4）：477-484.

［9］ 王雪璐.青藏高原三江源高寒草地生态系统土壤侵蚀研究［D］.兰州：兰州大学，2016.

WANG X L. Soil esrosion of alpine grassland ecosystem on Tibetan Plateau Sanjiangyuan region［D］.Lanzhou：Lanzhou University，2016.

［10］ [ZK（#]刘晓东，刘荣堂，刘爱军，等.三江源区草地覆盖遥感信息提取方法及动态研究［J］.草地学报，2010，18（2）：154-159.

LIU X D，LIU R T，LIU A J，et al. Study on information extraction and the dynamic  monitoring of grassland coverage in three river source area［J］. Acta Agrestia Sinica，2010，18（2）：154-159.

［11］ 安 如，姜丹萍，李晓雪，等.基于地面实测高光谱数据的三江源中东部草地植被光谱植被特征研究［J］.遥感技术与应用，2014，29（2）：202-211.

AN R，JIANG D P，LI X X，et al. Using Hyperspectral data to determine spectral characteristics of grassland vegetation in general and eastern parts of Three-river source［J］. Remote Sensing Technology and Application，2014，29（2）：202-211.

［12］ 魏卫东，李希来.三江源区高寒草地退化草地土壤侵蚀模型与模拟研究［J］.环境科学与管理，2013，38（7）：26-30.

WEI W D，LI X L. Soil erosion model and simulation of degraded grassland of alpine meadow in Sanjiangyuan region［J］. Environmental Science and Management，2013，38（7）：26-30.

［13］ 蔣 冲，高艳妮，李 芬，等.1956-2010年三江源区水土流失状况演变［J］.环境科学研究，2017，30（1）：20-29.

JIANG CH，GAO Y N，LI F，et al. Soil erosion variation and its attribution in the Three-Rivers Headwater region，China，during the period of 1956-2010［J］. Research of Environment Science，2017，30（1）：20-29.

［14］ 朱 霞，钞振华，杨永顺，等.三江源区“黑土滩”型退化草地时空变化［J］.草业科学2014，31（9）：1628-1636.

ZHU X，CHAO ZH H. YANG Y SH，et al. Spatial and temporal in the Three River source type degraded grassland in the Three River source region［J］. Pratacultural Science，2014，31（9）：1628-1636.

［15］ 徐剑波，宋立生，赵之重，等.近15a来黄河源地区玛多县草地植被退化的遥感动态监测［J］.干旱区地理，2014，  35（4）：615-622.

XU J B，SONG L SH，ZHAO ZH ZH，et al.  Monitoring grassland degradation dynamically at Maduo county in source region of Yellow River in past 15 years based on remote sensing［J］. Arid Land Geography，2014，35（4）：615-622.

［16］ 王俊奇，王广军，梁四海，等.1996-2015年黄河源区植被覆盖度提取和时空变化分析［J］.冰川冻土，2021，43（2）：662-674.

WANG J Q，WANG G J，LIANG S H，et al. Extraction of vegetation coverage and analysis of temporal and spatial changes in the source area of the Yellow River from 1996 to 2015［J］. Journal of Glaciology and Geocryology，2021，43（2）：662-674.

［17］ 柴 瑜，李希来，于金峰，等.有机肥施用量对黄河源不同坡向退化高寒草地土壤团聚体及有机碳的影响［J］.草地学报，2022，30（7）：1613-1620.

CHAI  Y，LI X L，YU J F，et al.Effects of different applications of organic fertilizer in degraded alpine meadowon soil aggregates and organic carbon in the Source Zone of Yellow River［J］. Acta Agrestia Sinica，2022，30（7）：1613-1620.

［18］ 张群英，李 捷，郝力壮，等.体外法评价青海省泽库县高寒草地草场不同物候期牧草的营养价值［J］.动物营养学报，2020，32（3）：1415-1432.

ZHANG Q Y，LI J，HAO L ZH，et al.Evaluation of forage nutritional value of alpine meadow grassland in different phenological periods in Zeku County of Qinghai Province by in vitro method［J］.Chinese Journal of Animal Nutrition，2020，32（3）：1415-1432.

［19］ 韩 俊. 不同的遥感解译方法在地表覆盖分类信息提取中的应用研究［D］.昆明：昆明理工大学，2016.

HAN J.Application of different remote sensing interpretation methods in surface cover classification information extraction［D］. Kunming：Kunming University of Science and Technology，2016.

［20］ 刘海军，魏晓琴.3S技术在三江源区“黑土滩”退化草地调查中的应用［J］.青海国土经略，2007（3）：29-31.

LIU H J，WEI X Q. Application of 3S technology in the investigation of “black soil beach” degraded grassland in the source area of Three-Rivers［J］.Qinghai Guotu Jinglue，2007（3）：29-31.

［21］ 李永花，王 苑，赵成福.三江源区“黑土滩”的遥感图像解译［J］.测绘技术装备，2008，10（4）：38-41.

LI Y H，WANG Y，ZHAO CH F. Remote sensing image interpretation of black soil beach in the source area of Three Rivers［J］.Geomatics Technology and Equipment，2008，10（4）：38-41.

［22］ 董全民，馬玉寿，许长军，等.三江源区“黑土滩”退化草地分类分级体系及分类恢复研究［J］.草地学报，2015，23（3）：441-447.

DONG Q M，MA Y SH，XU CH J，et al. Study of classification and gradation，restoration of Black-Soil Beach degraded grassland in the Headwaters of Three rivers［J］. Acta Agrestia Sinica，2015，23（3）：441-447.

［23］ 李积兰，李希来，唐 燕，等. 三江源头“黑土滩”草地分布现状及退化特征研究［J］.草业与畜牧，2009（2）：6-13.

LI J L，LI X L，TANG Y，et al. Study on grassland distribution and degradation characteristics of “black soil beach” in the source of Three-rivers［J］. Caoye Yu Xumu，2009（2）：6-13.

［24］ 鲍士旦.土壤农化分析［M］.北京：中国农业出版社，2010.

BAO SH D. Soil Agrochemical Analysis ［M］.Beijing：China Agricultural Press，2010.

［25］ 李希来，李积兰. 三江源区优势种植物矮嵩草繁殖策略与环境适应［M］.武汉：长江出版社，2019.

LI X L，LI J L.Reproductive Strategy and Environmental Adaptation of Dominant Species Kobresia humiliss in the Source Area of Three-rivers［M］.Wuhan：Changjiang Press，2019.

［26］ 马玉寿，尚占环，施建军，等.黄河源区“黑土滩”退化草地群落类型多样性及其群落结构研究［J］.草业科学，2006，23（12）：6-11.

MA Y SH，SHANG ZH H，SHI J J，et al. Studies on communities diversity and their structure of “bleak-soil-land” degraded grsassland in the headwater of Yellow River［J］. Patacultural Science，2006，23（12）：6-11.

［27］ 刘 若.草原保护学［M］.北京：中国农业出版社，2005.

LIU R. Protection of Steppe［M］. Beijing：China Agriculture Press，2005.

［28］ 尚占环，龙瑞军.青藏高原“黑土型”退化草地成因与恢复［J］.生态学杂志，2005，24（6）：652-656.

SHANG ZH H，LONG R J. Formation reason and recovering problem of the ‘black soil type degradation alpine in Qinghai-Tibetan plateau［J］. Chinese of Journal of Ecology，2005，24（6）：652-656.

［29］ 尚占环，董全民，施建军，等. 青藏高原“黑土滩”退化草地及其生态恢复近10年研究进展--兼论三江源生态恢复问题［J］.草地学报，2018，26（1）：1-21.

SHANG ZH H，DONG Q M，SHI J J，et al. Research progress in recent ten years of ecological restoration for ‘Black soil land degraded grassland on Tibetan plateau concurrently discuss of ecological restoration in Sangjiangyuan region［J］. Acta Agrestia Sinica，2018，26（1）：1-21.

［30］ 李希来.青藏高原“黑土滩”形成的自然因素与生物学机制［J］.草业科学，2002，19（1）：20-22.

LI X L. Natural factors and biological mechanism of the formation of “black soil beach” in the Qinghai Tibet Plateau［J］. Pratacultural Science，2002，19（1）：20-22.

［31］ 魏瑞琪，邵懷勇，李林峰，等.基于MODIS的川西北江河源区草地退化状况时空分析［J］.物探化探计算技术，2018，40（2）：262-268.

WEI R Q，SHAO H Y，LI L  F，et al.Temporal and spatial analysis of grassland degradation in yangtze and Yellow River Headwater region in Northwest Sichuan Based on MODIS［J］. Computing Techniques for Geophysical and Geoche，2018，40（2）：262-268.

［32］ 刘育红，魏卫东，杨元武，等. 高寒草地退化草地植被与土壤因子关系冗余分析［J］.西北农业学报，2018，27（4）：480-490.

LIU Y H，WEI W D，YANG Y W，et al. Redundancy   analysis on relationships between grassland vegetation and soil factors on degraded alpine meadow［J］. Agriculturae Boreali-occidentalis Sinica，2018，27（4）：480-490.

［33］ 杨振宇. 玛曲草原“黑土滩”分布特征及其植被演替规律的研究［D］. 兰州：兰州大学，2007.

YANG ZH Y. Study on distribution characteristics and vegetation succession rules of black soil beach in Maqu grassland［D］. Lanzhou：Lanzhou University，2007.

［34］ 李秋年. 高寒草地退化草地植物群落结构特征及物种多样性的初步分析［J］.青海环境，2004，14（1）：30-33.

LI Q N. Preliminary analysis of plant community structure and species diversity of degraded grassland in alpine meadow［J］.Qinghai Environment，2004，14（1）：30-33.

［35］ 许世洋，李雪萍，李敏权，等. 青藏高原高寒草地退化对草地群落多样性及土壤特性的影响［J］.中国草地学报，2022，44（8）：20-27.

XU SH Y，LI X P，LI M Q，et al. Effects of alpine mesdow degration on the grassland communities diversity and soil  properties on the Tibetan plateau［J］.Chinese Journal of Grassland，2022，44（8）：20-27.

［36］ 崔宁洁，张丹桔，刘 洋，等.不同林龄马尾松人工林林下植物多样性与土壤理化性质［J］.生态学杂志，2014，  33（10）：2610-2617.

CUI N J，ZHANG D J，LIU Y，et al. Plant diversity and soil physicochemical properties under different aged Pinus massoniana plantations［J］. Chinese of Journal of Ecology，2014，33（10）：2610-2617.

［37］ 楊红善，那·巴特尔，周学辉，等.不同放牧强度对肃北高寒草原土壤肥力的影响［J］.水土保持学报，2009，23（1）：150-153.

YANG H SH，NA·BATEER，ZHOU X H，et al. Influence of different grazing intensity to soil fertility in the Subei alpine steppes［J］. Journal of Soil and Water Conservation，2009，23（1）：150-153.

［38］ 刘育红，李希来，魏卫东，等.高寒草地“黑土滩”对土壤养分的影响［J］.西北农业学报，2009，18（3）：304-308.

LIU Y H，LI X L，WEI W D，et al. The influence of “Black Beach” to soils nutrient on Alpine Meadow［J］.Agriculturae Boreali-occidentalis Sinica，2009，18（3）：304-308.

［39］ 袁知洋，邓邦良，郭晓敏，等.武功山山地草地土壤全量氮磷钾分布格局及对不同退化程度的响应［J］.西北林学院学报，2015，30（3）：14-20.

YUAN ZH Y，DENG B L，GUO X M，et al. Distribution pattern of total nitrogen，phosphorus and potassium in mountain meadow of Wugong Mountain and its response to different degradation degrees［J］.Journal of Northwest Forestry University，2015，30（3）：14-20.

［40］ 王玉琴，尹亚丽，李世雄.不同退化程度高寒草甸土壤理化性质及酶活性分析［J］.生态环境学报，2019，28（6）：1108-1116.

WANG Y Q，YIN Y L，LI SH X.Physicochemical properties and enzymatic activities of alpine meadow  at different degradation degrees［J］.Ecology and Environmental Sciences，2019，28（6）：1108-1116.

［41］ LI X L，PERRY G L W，BRIERLEY G J，et al.Quantitative assessment of degradation classifications for degraded alpine meadows（Heitutan） ，Sanjiangyuan，western China［J］.Land Degradation ＆ Development，2014，25：417-427.

［42］ ZHANG Y，DONG S，GAO Q，et al. Responses of alpine vegetation and soils to the disturbance of plateau pike （Ochotona Curzon-iae） at burrow level on the Qinghai-Tibetan plateau of China［J］.Ecological Engineering，2016，88：232-236.

［43］ 林慧龙，王钊齐，尚占环.江河源区“黑土滩”退化草地秃斑与鼠洞的分形特征［J.草地学报，2010，18（4）：477-484.

LIN H L，WANG ZH Q，SHANG ZH  H.Features on fractal dimension of barren patch and mouse hole among different degenerated succession stages on alpine meadow in the source region of the Yangtze and Yellow river，Qinghai-Tibetan plateau，China［J］. Acta Agrestia Sinica，2010，18（4）：477-484.

［44］ 张海娟，谢久祥，李希来. 青海河南县高原鼠兔的分布及其与土壤特性的相关性分析［J］.青海大学学报（自然科学版），2016，34（1）：14-50.

ZHANG H J，XIE J X，LI X L.The distribution of the plateau pika and correlation withsoil properties of Henan County in Qinghai［J］. Journal of Qinghai University（Natural Science Edition），2016，34（1）：14-50.

Spatial Distribution Characteristics of Alpine Degraded Grassland in Source Region of Yellow River

Abstract In order to understand the  regularity of spatial distribution and the key topographic features of the degraded grassland in the black soil beach，the data of remote sensing image interpretation in the year of 2019 in counties of Maqin，Henan and Zeku in the source region of the Yellow River were used and field investigation was conducted to analyze the spatial distribution characteristics of alpine degraded grassland in black soil beach.The results showed that the degraded grassland in the black soil  beach of the source region of the Yellow River was distributed at the altitude of 3 500-4 000 m，and themost serious degradation was at the altitude of 3 700-3 800 m. From the slope analysis，the  degraded grassland in black soil beach mainly occurred on the sunny  and semi-sunny slopes，from the   analysis of slope degree，the degraded grassland in black soil  beach mainly occurred  on the gentle slope and beach.After the cluster analysis，the degraded grassland could be divided into six types，namely  extremely degraded grassland on the sunny slope，severely degraded grassland on the semi-sunny slope，slightly degraded grassland on the beach and sunny slope，slightly degraded grassland on the semi-sunny slope，slightly degraded grassland on the wetland，and non-degraded grassland on the shady slope. With the deterioration of degradation，the vegetation coverage and aboveground biomass decreased; the contents of total nitrogen，total phosphorus，alkali-hydrolysable nitrogen，available phosphorus and organic matter in the same soil layer of the different degradations were highest in  the slightly degraded grassland on the wetland and onthe shady slope，while it was the lowest in extremely degraded grassland on the sunny slope，and the differences were significant （P<0.05）. The richness index，diversity index，evenness index and soil available potassium content were higher in severely and slightly degraded grassland，howerver，it was lower in slightly degraded grassland on thewetland，and there were significant differences between slightly degraded grassland on the wetland  and other types of degraded grassland （P<0.05）. In the extremely degraded grassland on the sunny slope，the total potassium content and the rat damage rate were highest，the number of rat holes were the most ，the total potassium content was the lowest in slightly   degraded grassland on the the wetland，which was extremely different from the extremely degraded grassland on the sunny slope（P<0.05），and  there was no any rodent activity on the slightly degraded grassland on the shady  slope，so the characteristics of terrain characteristics was consistent with the conditions of rodents  habitat.

Key words Source region of the Yellow River; Degraded grassland; Spatial distribution; Plant community