Bare land exposed by glacier retreat provides new opportunities for ecosystem development. Investigating primary vegetation succession in deglaciated regions can provide significant insights for ecological restoration, particularly for future climate change scenarios. Nonetheless, research on this topic in the Qinghai-Tibet Plateau has been exceedingly limited. This study aimed to investigate vegetation succession in the deglaciated area of the Zepu glacier during the Little Ice Age in southeastern Tibet. Quadrat surveys were performed on arboreal communities, and trends in vegetation change were assessed utilizing multi-year (1986-2024) remote sensing data. The findings indicate that vegetation succession in the Zepu glacier deglaciated area typically adheres to a sequence of bare land-shrub-tree, divided into four stages: (1) shrub (species include Larix griffithii Mast., Hippophae rhamnoides subsp. yunnanensis Rousi, Betula utilis D. Don, and Populus pseudoglauca C. Wang & P. Y. Fu); (2) broadleaf forest primarily dominated by Hippophae rhamnoides subsp. yunnanensis Rousi; (3) mixed coniferous-broadleaf forest with Hippophae rhamnoides subsp. yunnanensis Rousi and Populus pseudoglauca C. Wang & P. Y. Fu as the dominant species; and (4) mixed coniferous-broadleaf forest dominated by Picea likiangensis (Franch.) E. Pritz. Soil depth and NDVI both increase with succession. Species diversity is significantly higher in the third stage compared to other successional stages. In addition, soil moisture content is significantly greater in the broadleaf-dominated communities than in the conifer-dominated communities. An analysis of NDVI from 1986 to 2024 reveals an overall positive trend in vegetation recovery in the area, with 93% of the area showing significant vegetation increase. Temperature is the primary controlling factor for this recovery, showing a positive correlation with vegetation cover. The results indicate that Key ecological indicators-including species composition, diversity, NDVI, soil depth, and soil moisture content-exhibit stage-specific patterns, reflecting distinct phases of primary succession. These findings enhance our comprehension of vegetation succession in deglaciated areas and their influencing factors in deglaciated areas, providing theoretical support for vegetation restoration in climate change.

本文使用1988—2021年较高空间分辨率的Landsat卫星遥感数据细致分析了位于青藏高原腹地的长江源区植被绿度变化情况，并使用地表长期形变速率、热融湖塘及热融滑塌分布定量描述多年冻土退化状态并由此明晰植被对多年冻土退化的响应。研究表明：(1)长江源区NDVI自1988年以来平均增加速率为0.003 1/a，绿化区域占总源区的91.9%，褐化区域占总源区的7.1%；(2)随着地表沉降速率增大，植被绿化速率加快，但沉降速率高于25 mm/a时，NDVI增长速率逐渐减缓，当沉降速率进一步增加时，甚至一些地区表现出了明显的NDVI减小情况；(3)热融滑塌的出现显著破坏了发生区域内植被，但周围120 m范围内植被发育情况好于全源区平均水平，热融湖塘区及其周围60 m范围内的植被绿化程度落后于全源区平均水平。长江源植被区域褐化比例为7.1%[只统计褐化植被/植被区域（NDVI>0.1部分）]，其中多年冻土退化直接影响了29.6%的褐化植被（只统计褐化植被中有明显冻土影响信号的部分占褐化区域比重）。

本文使用1988—2021年较高空间分辨率的Landsat卫星遥感数据细致分析了位于青藏高原腹地的长江源区植被绿度变化情况，并使用地表长期形变速率、热融湖塘及热融滑塌分布定量描述多年冻土退化状态并由此明晰植被对多年冻土退化的响应。研究表明：(1)长江源区NDVI自1988年以来平均增加速率为0.003 1/a，绿化区域占总源区的91.9%，褐化区域占总源区的7.1%；(2)随着地表沉降速率增大，植被绿化速率加快，但沉降速率高于25 mm/a时，NDVI增长速率逐渐减缓，当沉降速率进一步增加时，甚至一些地区表现出了明显的NDVI减小情况；(3)热融滑塌的出现显著破坏了发生区域内植被，但周围120 m范围内植被发育情况好于全源区平均水平，热融湖塘区及其周围60 m范围内的植被绿化程度落后于全源区平均水平。长江源植被区域褐化比例为7.1%[只统计褐化植被/植被区域（NDVI>0.1部分）]，其中多年冻土退化直接影响了29.6%的褐化植被（只统计褐化植被中有明显冻土影响信号的部分占褐化区域比重）。

本文使用1988—2021年较高空间分辨率的Landsat卫星遥感数据细致分析了位于青藏高原腹地的长江源区植被绿度变化情况，并使用地表长期形变速率、热融湖塘及热融滑塌分布定量描述多年冻土退化状态并由此明晰植被对多年冻土退化的响应。研究表明：(1)长江源区NDVI自1988年以来平均增加速率为0.003 1/a，绿化区域占总源区的91.9%，褐化区域占总源区的7.1%；(2)随着地表沉降速率增大，植被绿化速率加快，但沉降速率高于25 mm/a时，NDVI增长速率逐渐减缓，当沉降速率进一步增加时，甚至一些地区表现出了明显的NDVI减小情况；(3)热融滑塌的出现显著破坏了发生区域内植被，但周围120 m范围内植被发育情况好于全源区平均水平，热融湖塘区及其周围60 m范围内的植被绿化程度落后于全源区平均水平。长江源植被区域褐化比例为7.1%[只统计褐化植被/植被区域（NDVI>0.1部分）]，其中多年冻土退化直接影响了29.6%的褐化植被（只统计褐化植被中有明显冻土影响信号的部分占褐化区域比重）。

本文使用1988—2021年较高空间分辨率的Landsat卫星遥感数据细致分析了位于青藏高原腹地的长江源区植被绿度变化情况，并使用地表长期形变速率、热融湖塘及热融滑塌分布定量描述多年冻土退化状态并由此明晰植被对多年冻土退化的响应。研究表明：(1)长江源区NDVI自1988年以来平均增加速率为0.003 1/a，绿化区域占总源区的91.9%，褐化区域占总源区的7.1%；(2)随着地表沉降速率增大，植被绿化速率加快，但沉降速率高于25 mm/a时，NDVI增长速率逐渐减缓，当沉降速率进一步增加时，甚至一些地区表现出了明显的NDVI减小情况；(3)热融滑塌的出现显著破坏了发生区域内植被，但周围120 m范围内植被发育情况好于全源区平均水平，热融湖塘区及其周围60 m范围内的植被绿化程度落后于全源区平均水平。长江源植被区域褐化比例为7.1%[只统计褐化植被/植被区域（NDVI>0.1部分）]，其中多年冻土退化直接影响了29.6%的褐化植被（只统计褐化植被中有明显冻土影响信号的部分占褐化区域比重）。

本文使用1988—2021年较高空间分辨率的Landsat卫星遥感数据细致分析了位于青藏高原腹地的长江源区植被绿度变化情况，并使用地表长期形变速率、热融湖塘及热融滑塌分布定量描述多年冻土退化状态并由此明晰植被对多年冻土退化的响应。研究表明：(1)长江源区NDVI自1988年以来平均增加速率为0.003 1/a，绿化区域占总源区的91.9%，褐化区域占总源区的7.1%；(2)随着地表沉降速率增大，植被绿化速率加快，但沉降速率高于25 mm/a时，NDVI增长速率逐渐减缓，当沉降速率进一步增加时，甚至一些地区表现出了明显的NDVI减小情况；(3)热融滑塌的出现显著破坏了发生区域内植被，但周围120 m范围内植被发育情况好于全源区平均水平，热融湖塘区及其周围60 m范围内的植被绿化程度落后于全源区平均水平。长江源植被区域褐化比例为7.1%[只统计褐化植被/植被区域（NDVI>0.1部分）]，其中多年冻土退化直接影响了29.6%的褐化植被（只统计褐化植被中有明显冻土影响信号的部分占褐化区域比重）。

本文使用1988—2021年较高空间分辨率的Landsat卫星遥感数据细致分析了位于青藏高原腹地的长江源区植被绿度变化情况，并使用地表长期形变速率、热融湖塘及热融滑塌分布定量描述多年冻土退化状态并由此明晰植被对多年冻土退化的响应。研究表明：(1)长江源区NDVI自1988年以来平均增加速率为0.003 1/a，绿化区域占总源区的91.9%，褐化区域占总源区的7.1%；(2)随着地表沉降速率增大，植被绿化速率加快，但沉降速率高于25 mm/a时，NDVI增长速率逐渐减缓，当沉降速率进一步增加时，甚至一些地区表现出了明显的NDVI减小情况；(3)热融滑塌的出现显著破坏了发生区域内植被，但周围120 m范围内植被发育情况好于全源区平均水平，热融湖塘区及其周围60 m范围内的植被绿化程度落后于全源区平均水平。长江源植被区域褐化比例为7.1%[只统计褐化植被/植被区域（NDVI>0.1部分）]，其中多年冻土退化直接影响了29.6%的褐化植被（只统计褐化植被中有明显冻土影响信号的部分占褐化区域比重）。

本文使用1988—2021年较高空间分辨率的Landsat卫星遥感数据细致分析了位于青藏高原腹地的长江源区植被绿度变化情况，并使用地表长期形变速率、热融湖塘及热融滑塌分布定量描述多年冻土退化状态并由此明晰植被对多年冻土退化的响应。研究表明：(1)长江源区NDVI自1988年以来平均增加速率为0.003 1/a，绿化区域占总源区的91.9%，褐化区域占总源区的7.1%；(2)随着地表沉降速率增大，植被绿化速率加快，但沉降速率高于25 mm/a时，NDVI增长速率逐渐减缓，当沉降速率进一步增加时，甚至一些地区表现出了明显的NDVI减小情况；(3)热融滑塌的出现显著破坏了发生区域内植被，但周围120 m范围内植被发育情况好于全源区平均水平，热融湖塘区及其周围60 m范围内的植被绿化程度落后于全源区平均水平。长江源植被区域褐化比例为7.1%[只统计褐化植被/植被区域（NDVI>0.1部分）]，其中多年冻土退化直接影响了29.6%的褐化植被（只统计褐化植被中有明显冻土影响信号的部分占褐化区域比重）。

本文使用1988—2021年较高空间分辨率的Landsat卫星遥感数据细致分析了位于青藏高原腹地的长江源区植被绿度变化情况，并使用地表长期形变速率、热融湖塘及热融滑塌分布定量描述多年冻土退化状态并由此明晰植被对多年冻土退化的响应。研究表明：(1)长江源区NDVI自1988年以来平均增加速率为0.003 1/a，绿化区域占总源区的91.9%，褐化区域占总源区的7.1%；(2)随着地表沉降速率增大，植被绿化速率加快，但沉降速率高于25 mm/a时，NDVI增长速率逐渐减缓，当沉降速率进一步增加时，甚至一些地区表现出了明显的NDVI减小情况；(3)热融滑塌的出现显著破坏了发生区域内植被，但周围120 m范围内植被发育情况好于全源区平均水平，热融湖塘区及其周围60 m范围内的植被绿化程度落后于全源区平均水平。长江源植被区域褐化比例为7.1%[只统计褐化植被/植被区域（NDVI>0.1部分）]，其中多年冻土退化直接影响了29.6%的褐化植被（只统计褐化植被中有明显冻土影响信号的部分占褐化区域比重）。

本文使用1988—2021年较高空间分辨率的Landsat卫星遥感数据细致分析了位于青藏高原腹地的长江源区植被绿度变化情况，并使用地表长期形变速率、热融湖塘及热融滑塌分布定量描述多年冻土退化状态并由此明晰植被对多年冻土退化的响应。研究表明：(1)长江源区NDVI自1988年以来平均增加速率为0.003 1/a，绿化区域占总源区的91.9%，褐化区域占总源区的7.1%；(2)随着地表沉降速率增大，植被绿化速率加快，但沉降速率高于25 mm/a时，NDVI增长速率逐渐减缓，当沉降速率进一步增加时，甚至一些地区表现出了明显的NDVI减小情况；(3)热融滑塌的出现显著破坏了发生区域内植被，但周围120 m范围内植被发育情况好于全源区平均水平，热融湖塘区及其周围60 m范围内的植被绿化程度落后于全源区平均水平。长江源植被区域褐化比例为7.1%[只统计褐化植被/植被区域（NDVI>0.1部分）]，其中多年冻土退化直接影响了29.6%的褐化植被（只统计褐化植被中有明显冻土影响信号的部分占褐化区域比重）。