Vegetation Growth Status and Topographic Effects in Frozen Soil Regions on the Qinghai–Tibet Plateau

: The Qinghai–Tibet Plateau (QTP), which is known as Earth’s “Third Pole”, is a driver of global climate change, and long-term monitoring of QTP vegetation can reveal changes attributable to climatic and human inﬂuences. Previous research monitoring vegetation on the QTP focused primarily on spatiotemporal variations of vegetation indices, while few studies have considered fractional vegetation cover (FVC) in relation to topographic and frozen soil factors. We used MODIS-EVI, digital elevation models, and frozen soil data to investigate topographic effects on vegetation growth status in different soil types on the QTP during 2000–2020. (1) FVC showed a trend of increase during 2000–2020, and the FVC on the QTP decreased from the southeast to the northwest in spatial distribution. FVC in permafrost regions was the lowest, followed by seasonal frozen soil areas; FVC in unfrozen areas was the highest. (2) With increasing elevation, FVC of permafrost, seasonal frozen, and unfrozen soil areas showed downward trends for each aspect. In seasonal frozen soil areas, at elevation ≤ 4000 m (>4000 m), FVC of sunny (shady) slopes was greater than that of shady (sunny) slopes. In permafrost regions, except at elevations of 3000–4000 m, FVC of shady slopes was greater than that of sunny slopes. In unfrozen soil areas, at elevation >4000 m, FVC of sunny slopes was obviously greater than that of shady slopes. (3) With increasing slope, FVC in seasonal frozen and permafrost soil (unfrozen soil) regions showed a trend of increase (decrease). In seasonal frozen soil areas, FVC of sunny (shady) slopes was greater than that of shady (sunny) slopes on slopes ≤ 6 ◦ (>6 ◦ ). In permafrost regions, FVC of sunny slopes was less than that of shady slopes. With increasing slope, the inﬂuence of aspect became more obvious. In unfrozen soil areas, FVC of sunny slopes was slightly greater than that of shady slopes. Topographic effects especially the elevation and slope effects might signiﬁcantly affect the spatiotemporal variations of vegetation growth status in frozen soil regions on the QTP.


Introduction
The Qinghai-Tibet Plateau (QTP) is the highest plateau in the world and an ecological security barrier in China. It plays an important role in biodiversity protection [1], soil and water conservation [2], and climate change [3]. Frozen soil, which is distributed widely on the QTP and is very sensitive to climate change, can affect not only the local ecology and climatic system but also the global carbon cycle and climate change [4,5]. Thawing of the active layer of permafrost increases soil moisture that can trigger dramatic change from sensible heat flux to latent heat flux in areas of both alpine meadow and bare land [6]. The permafrost layer is considered an impermeable layer, and freezing and thawing processes can change the water phase and water cycle. Adaption or succession will occur in vegetation relying on shallow groundwater with change of the hydrological cycle [7].

Study Area
The QTP is located in the southwest of China (27.5°-37.5°N, 75°-105°E). It ran from the Kunlun Mountains in the north to the Himalayas in the south, and from the K rakoram Mountains in the west to the Hengduan Mountains in the east. The QTP, wh includes Tibet, Qinghai Provinces, and parts of Xinjiang, Gansu, Sichuan, and Yunn Provinces, covers an area of 2.57 × 10 6 km 2 . The QTP is the highest plateau in the wor with average elevation of 4300 m. Influenced by the East Asian, Indian Southwest, a Westerly monsoon circulation systems, the distribution of both water and heat on the Q has distinct variation [7]. Annual average temperature is −2.0 °C. Annual average prec itation is 340 mm, and its distribution is very uneven, with summer rainfall account for more than 70% of the total. Frozen soil is distributed widely on the QTP, cover approximately 96% of the total area, which makes the QTP the largest low-latitude p mafrost region in the world [31]. The main types of soil on the QTP are seasonal froz soil, permafrost soil, and unfrozen soil ( Figure 1). The areas of permafrost soil and s sonal frozen soil account for 40% and 56% of the total area, respectively. There are vario types of vegetation on the QTP, including coniferous forest, coniferous and broadleav mixed forest, broadleaved forest, shrub, desert, grassland, grass, meadow, swamp, a cultivated vegetation. Among them, grassland covers the largest area (27.

Data
To continuously and completely realize dynamic monitoring of vegetation grow on the QTP, MOD13Q1-EVI images acquired during 2000-2020 (March-October) w used in this study. The data were obtained from NASA's Land Processes Distributed A tive Archive Center (LP DAAC) located at the USGS Earth Resources Observation a Science (EROS) Center (https://e4ftl01.cr.usgs.gov/, accessed on 25 March 2021). The s tial resolution of MODI13Q1 data is 250 m and the temporal resolution is 16 d. The MOD Reprojection Tool was used for mosaic and projection transformation of the images, a the projection type adopted was Albers Equal Area Conic projection (central longitu 105°E, latitude: 25°N and 47°N). The international maximum value composites meth was used to calculate the largest EVI image sequences during 2000-2020.
Digital elevation model (DEM) data were used to extract the topography of the Q The DEM data with the 30 m spatial resolution were obtained from the Geospatial D Cloud site, Computer Network Information Center of the Chinese Academy of Scien (http: www.gscloud.cn, accessed on 6 December 2021). Slope and aspect were calcula from the DEM data. The elevation of the QTP was divided into six levels: ≤2000, 20

Data
To continuously and completely realize dynamic monitoring of vegetation growth on the QTP, MOD13Q1-EVI images acquired during 2000-2020 (March-October) were used in this study. The data were obtained from NASA's Land Processes Distributed Active Archive Center (LP DAAC) located at the USGS Earth Resources Observation and Science (EROS) Center (https://e4ftl01.cr.usgs.gov/, accessed on 25 March 2021). The spatial resolution of MODI13Q1 data is 250 m and the temporal resolution is 16 d. The MODIS Reprojection Tool was used for mosaic and projection transformation of the images, and the projection type adopted was Albers Equal Area Conic projection (central longitude: 105 • E, latitude: 25 • N and 47 • N). The international maximum value composites method was used to calculate the largest EVI image sequences during 2000-2020.

Fractional Vegetation Coverage
Fractional vegetation coverage (FVC) refers to the percent of the vertical projection area covered by vegetation as a fraction of the total area [33]. Generally, FVC is calculated using remote sensing techniques that include mixed-pixel decomposition models, regression models, machine learning methods, and vegetation canopy models [34]. Among those, mixed-pixel decomposition models are used most widely, and the pixel binary model is the most representative of that class. The pixel binary model is relatively simple and has obvious physical advantages in that it is better able to weaken the influence of the atmosphere, soil background, and vegetation type [25,35]. The pixel binary model assumes that each pixel can be decomposed into full vegetation cover and bare soil [36]. The Gutman model and the Carlson model are the pixel binary models used most widely, owing to their applicability to areas with high and low FVC, respectively [37,38]. Owing to the characteristics of sparse vegetation on the QTP, the Carlson model was used to estimate FVC in this study. As the basic principle of the Carlson model is to assume that each pixel can be decomposed into pure vegetation and pure soil, the information of pixel (S) can be expressed as the sum of the vegetation cover (Sv) and the soil cover (Ss): For a pixel containing only soil and vegetation, the area occupied by the vegetation is the FVC, and the proportion of the area occupied by soil is 1 − FVC; therefore, FVC can be expressed as follows:

Fractional Vegetation Coverage
Fractional vegetation coverage (FVC) refers to the percent of the vertical projection area covered by vegetation as a fraction of the total area [33]. Generally, FVC is calculated using remote sensing techniques that include mixed-pixel decomposition models, regression models, machine learning methods, and vegetation canopy models [34]. Among those, mixed-pixel decomposition models are used most widely, and the pixel binary model is the most representative of that class. The pixel binary model is relatively simple and has obvious physical advantages in that it is better able to weaken the influence of the atmosphere, soil background, and vegetation type [25,35]. The pixel binary model assumes that each pixel can be decomposed into full vegetation cover and bare soil [36]. The Gutman model and the Carlson model are the pixel binary models used most widely, owing to their applicability to areas with high and low FVC, respectively [37,38]. Owing to the characteristics of sparse vegetation on the QTP, the Carlson model was used to estimate FVC in this study. As the basic principle of the Carlson model is to assume that each pixel can be decomposed into pure vegetation and pure soil, the information of pixel (S) can be expressed as the sum of the vegetation cover (S v ) and the soil cover (S s ): For a pixel containing only soil and vegetation, the area occupied by the vegetation is the FVC, and the proportion of the area occupied by soil is 1 − FVC; therefore, FVC can be expressed as follows: where EVI soil and EVI veg represent the EVI of pure soil pixels and pure vegetation pixels, respectively. For bare soil pixels, EVI soil should be close to 0, but because of the influence of atmospheric, surface wetland, and sunlight conditions, values range from −0.1 to 0.2 [37,39]. For pure vegetation pixels, EVI veg varies depending on vegetation type, composition, and seasonal change in growth. Therefore, we used a 0.5% confidence level, and we selected EVI accumulation values of 0.5% and 99.5% as EVI soil and EVI veg , respectively. Both EVI soil and EVI veg were used in Equation (2) to calculate FVC during 2000-2020. In accordance with the Standard for Classification and Gradation of Soil Erosion of China (SL190-2007), and with consideration of the actual situation on the QTP, we divided FVC into five levels: bare land (<10%), low FVC (10%-30%), medium FVC (30%-50%), high FVC (50%-70%), and extremely high FVC (≥70%).

Sen Trend Analysis and Mann-Kendall Test
Sen trend analysis and the Mann-Kendall test are important methods used widely for determining the trend of long time series data in vegetation analysis [40,41]. In comparison with the least squares method, the Sen trend can avoid interference from missing data, elimination of outliers, and form of data distribution on the analysis results. The Sen trend is expressed as follows: where Slope FVC is the trend of the time series variables; a positive (negative) value of SlopeFVC indicates a trend of increase (decrease) in FVC. The Mann-Kendall test is a nonparametric statistical test recommended by the World Meteorological Organization, which is used widely to determine the statistical significance of a trend. The Mann-Kendall test does not require samples to obey a certain distribution and the result is not disturbed by a few outliers. The method has high quantification, wide detection range, and convenient calculation, making it suitable for objective assessment of the trend of time series of sequential variables and type variables [42,43]. The Mann-Kendall test is formulated as described in the following.
Let a stationary independent sequence be X t (t = 1..., n, n is the sequence length), where S is defined as follows: where sgn is a symbolic function, when n ≥ 10, S approximately obeys the normal distribution, and the mean E(S) = 0 and variance V(S) = n (n − 1)(2n + 5)/18, respectively, without considering the equivalent data points in the sequence. The value of test statistic Z of a positive distribution can be calculated using Equation (6): when Z > 0 (Z < 0), the sequence has an upward (downward) trend. Considering the symmetry of a normal distribution, the significance level is given as α; when |Z| > Z(1 − α/2),

Coefficient of Variation
The coefficient of variation (CV) is a normalized measure of the dispersion degree of a probability distribution and is the ratio of the variation to the mean, which can better evaluate the variability of vegetation change associated with topographic effects [29,44].
where MD represents the difference in mean FVC between sunny and shady slopes and Mean represents the mean FVC for each aspect. When CV is close to 0, the difference in FVC between shady and sunny slopes is smaller. A negative (positive) CV value indicates that the sunny (shady) slope FVC is smaller than the shady (sunny) slope FVC.

Correlation Analysis
Correlation analysis is a method for analyzing the degree of correlation between two or more variables. In this study, correlation analysis was used to analyze the relationship between FVC and temperature or precipitation: where r xy is the correlation coefficient that has values in the range of −1 and 1 (an absolute value closer to 1 represents stronger correlation), and x and y represent the mean of variables x and y, respectively.

Spatial Variation of FVC
The annual mean FVC on the QTP during 2000-2020 was 0.494 ( Figure 3), and it showed a trend of decrease from the southeast to the northwest. The area of low coverage (FVC < 0.3) was largest, accounting for 34.84% of the total area, and it was distributed mainly in southwestern, central, and northern parts of the QTP. The area with extremely high coverage (FVC > 0.7) accounted for 33.62% of the total area, and it was distributed mainly in southern and eastern parts of the QTP. The area of high coverage (0.5 ≤ FVC < 0.7) accounted for 15.94% of the total area, and it was distributed mainly in central and southern parts of the QTP. The area of medium coverage (0.3 ≤ FVC < 0.5) accounted for 15.60% of the total area, and it was distributed mainly in southwestern and central parts of the QTP.
Remote Sens. 2022, 14, x FOR PEER REVIEW 7 of 1 mainly in southern and eastern parts of the QTP. The area of high coverage (0.5 ≤ FVC 0.7) accounted for 15.94% of the total area, and it was distributed mainly in central and southern parts of the QTP. The area of medium coverage (0.3 ≤ FVC < 0.5) accounted fo 15.60% of the total area, and it was distributed mainly in southwestern and central part of the QTP. The spatial change trend of mean annual FVC on the QTP during 2000-2020 (Figur 4a) showed that the trend value of FVC was −0.0481 to 0.047. The area with an increasing trend of FVC accounted for 54.32% of the total area, distributed mainly in central and eastern parts of the QTP, including southern Xinjiang Province, and northeastern, north western, and central parts of Gansu and Qinghai provinces. The area with a decreasing trend of FVC accounted for 45.68% of the total area, distributed mainly in southern part of the QTP, including western Sichuan Province, northwestern Yunnan Province, centra and southeastern Tibet, and southwestern and southeastern regions of Qinghai Province The significance test (Figure 4b) of FVC on the QTP showed that the change in mean FVC was not significant in most areas (81.65% of the total area), exhibited significant and ex tremely significant increase in 6.62% and 5.54% of the area, respectively, and exhibited significant and extremely significant decrease in 4.43% and 2.21% of the area, respectively  For different soil types, the mean trend values of spatial change rate ( Figure 4a) fo seasonal frozen soil, permafrost soil, and unfrozen soil regions were 0.0002, 0.0004, and 0.001, respectively, indicating that the FVC in different soil type regions showed gradua increasing trends and the mean trend of unfrozen soil was relatively higher than those o permafrost and seasonal frozen soil. In terms of vegetation growth trends, seasonal frozen soil exhibited a trend of gradual decline in 52.48% of the total area, distributed mainly in southeastern and southern Xinjiang Province, northern Tibet, and central and western The spatial change trend of mean annual FVC on the QTP during 2000-2020 ( Figure 4a) showed that the trend value of FVC was −0.0481 to 0.047. The area with an increasing trend of FVC accounted for 54.32% of the total area, distributed mainly in central and eastern parts of the QTP, including southern Xinjiang Province, and northeastern, northwestern, and central parts of Gansu and Qinghai provinces. The area with a decreasing trend of FVC accounted for 45.68% of the total area, distributed mainly in southern parts of the QTP, including western Sichuan Province, northwestern Yunnan Province, central and southeastern Tibet, and southwestern and southeastern regions of Qinghai Province. The significance test (Figure 4b) of FVC on the QTP showed that the change in mean FVC was not significant in most areas (81.65% of the total area), exhibited significant and extremely significant increase in 6.62% and 5.54% of the area, respectively, and exhibited significant and extremely significant decrease in 4.43% and 2.21% of the area, respectively.  The spatial change trend of mean annual FVC on the QTP during 2000-2020 ( Figure  4a) showed that the trend value of FVC was −0.0481 to 0.047. The area with an increasing trend of FVC accounted for 54.32% of the total area, distributed mainly in central and eastern parts of the QTP, including southern Xinjiang Province, and northeastern, northwestern, and central parts of Gansu and Qinghai provinces. The area with a decreasing trend of FVC accounted for 45.68% of the total area, distributed mainly in southern parts of the QTP, including western Sichuan Province, northwestern Yunnan Province, central and southeastern Tibet, and southwestern and southeastern regions of Qinghai Province. The significance test (Figure 4b) of FVC on the QTP showed that the change in mean FVC was not significant in most areas (81.65% of the total area), exhibited significant and extremely significant increase in 6.62% and 5.54% of the area, respectively, and exhibited significant and extremely significant decrease in 4.43% and 2.21% of the area, respectively.  For different soil types, the mean trend values of spatial change rate (Figure 4a) for seasonal frozen soil, permafrost soil, and unfrozen soil regions were 0.0002, 0.0004, and 0.001, respectively, indicating that the FVC in different soil type regions showed gradual increasing trends and the mean trend of unfrozen soil was relatively higher than those of permafrost and seasonal frozen soil. In terms of vegetation growth trends, seasonal frozen soil exhibited a trend of gradual decline in 52.48% of the total area, distributed mainly in southeastern and southern Xinjiang Province, northern Tibet, and central and western For different soil types, the mean trend values of spatial change rate (Figure 4a) for seasonal frozen soil, permafrost soil, and unfrozen soil regions were 0.0002, 0.0004, and 0.001, respectively, indicating that the FVC in different soil type regions showed gradual increasing trends and the mean trend of unfrozen soil was relatively higher than those of permafrost and seasonal frozen soil. In terms of vegetation growth trends, seasonal frozen soil exhibited a trend of gradual decline in 52.48% of the total area, distributed mainly in southeastern and southern Xinjiang Province, northern Tibet, and central and western Qinghai Province. The area with a trend of gradual rise accounted for 47.52% of the total seasonal frozen soil area, distributed mainly in northwestern Yunnan Province, western Sichuan Province, southwestern Gansu Province, eastern Qinghai Province, and southern Tibet. In permafrost regions, the area with a trend of gradual decrease accounted for 44.77% of the total area, distributed mainly in southwestern parts of Tibet. The area with a trend of gradual increase accounted for 55.23% of the total permafrost area, distributed mainly in eastern-central parts of the border between Qinghai Province and Gansu Province, southwestern parts of Qinghai Province, and northern parts of Tibet. The unfrozen soil area showed a trend of decrease (increase) that accounted for 45.28% (54.72%) of the total area.

Effect of Topographic Variations
The topography of the QTP is complex, and the distribution and growth status of vegetation differ depending on terrain. Among various topographic factors, elevation, slope, and aspect are considered to best reflect the variation of vegetation conditions, and they are used widely in evaluations of vegetation change [27,45]. Thus, this study investigated the effects of elevation, slope, and aspect on vegetation growth on the QTP.

Effect of Elevation and Aspect
The areal proportions of FVC at different elevation levels on the QTP varied greatly ( Figure 5). When elevation was <4000 m, the areal proportion of extremely high coverage was highest, i.e., the values at elevations of ≤2000 and 3000-4000 m were 91.7% and 62.1%, respectively. No area had extremely high coverage at elevation >6000 m. The area of extremely high coverage showed a trend of gradual decrease with increasing elevation. The areal proportion of high coverage at different elevation levels was relatively small; the highest proportion (19.42%) was distributed at elevations of 4000-5000 m, and the lowest proportion (3.99%) was distributed at elevations ≤2000 m. No area had high coverage at elevation >6000 m. The areal proportion of high coverage exhibited a trend of gradual increase and then decrease with increasing elevation. The areal proportion of middle coverage was relatively low at each elevation level; the maximum was 20.72% at elevations of 5000-6000 m and the minimum was 1.92% when elevation was 2000 m. No area had middle coverage at elevation >6000 m. The areal proportion of low coverage was maximum when elevation was >6000 m. With increasing elevation, the areal proportion of low coverage gradually increased. With increase of elevation, the areal proportion of extremely high and low coverage on the QTP decreased gradually and increased, respectively, which might be related to the decrease in temperature with increasing elevation. Thus, temperature might be an important limiting factor for plant growth on the QTP [21].

Effect of Topographic Variations
The topography of the QTP is complex, and the distribution and growth status of vegetation differ depending on terrain. Among various topographic factors, elevation, slope, and aspect are considered to best reflect the variation of vegetation conditions, and they are used widely in evaluations of vegetation change [45,27]. Thus, this study investigated the effects of elevation, slope, and aspect on vegetation growth on the QTP.

Effect of Elevation and Aspect
The areal proportions of FVC at different elevation levels on the QTP varied greatly ( Figure 5). When elevation was <4000 m, the areal proportion of extremely high coverage was highest, i.e., the values at elevations of ≤2000 and 3000-4000 m were 91.7% and 62.1%, respectively. No area had extremely high coverage at elevation >6000 m. The area of extremely high coverage showed a trend of gradual decrease with increasing elevation. The areal proportion of high coverage at different elevation levels was relatively small; the highest proportion (19.42%) was distributed at elevations of 4000-5000m, and the lowest proportion (3.99%) was distributed at elevations ≤2000 m. No area had high coverage at elevation >6000 m. The areal proportion of high coverage exhibited a trend of gradual increase and then decrease with increasing elevation. The areal proportion of middle coverage was relatively low at each elevation level; the maximum was 20.72% at elevations of 5000-6000 m and the minimum was 1.92% when elevation was 2000 m. No area had middle coverage at elevation >6000 m. The areal proportion of low coverage was maximum when elevation was >6000 m. With increasing elevation, the areal proportion of low coverage gradually increased. With increase of elevation, the areal proportion of extremely high and low coverage on the QTP decreased gradually and increased, respectively, which might be related to the decrease in temperature with increasing elevation. Thus, temperature might be an important limiting factor for plant growth on the QTP [21].   Figure 6). In unfrozen soil areas, FVC was largest at elevation ≤5000 m; the maximum value (0.929) was at elevation ≤2000 m. As elevation increased, FVC decreased gradually. In seasonal frozen soil areas, FVC was highest (0.812) at elevation ≤2000 m. With increasing elevation, the variation of FVC showed a trend of decline, followed by increase, and then decline. In permafrost regions, FVC mainly decreased with increasing elevation. At all elevations except 2000-3000 m, FVC values of permafrost regions were smaller than those of seasonal frozen soil regions. At elevation >6000 m, there was only permafrost and FVC was relatively low (0.038). The variation of FVC with different elevation and aspect is presented in Table 1. The FVC of both permafrost soil and unfrozen soil exhibited a downward trend with increasing elevation. The FVC of seasonal frozen soil generally decreased gradually with increasing elevation but with a slight increase at elevations of 3000-4000 m. In seasonal frozen soil, when elevation was ≤4000 m, FVC of sunny slopes was larger than that of shady slopes, and FVC of semi-sunny slopes was larger than that of semi-shady slopes. At elevation >4000 m, FVC of shady slopes was larger than that of sunny slopes, and FVC of semi-shady slopes was larger than that of semi-sunny slopes. In permafrost regions, FVC of semi-shady slopes was larger than that of semi-sunny slopes. Except at elevations of 3000-4000 m, FVC of shady slopes was larger than that of sunny slopes. In unfrozen soil areas, at elevations of 2000-4000 m, FVC of shady slopes was larger than that of sunny slopes, and at elevations ≤2000 or ≥4000 m, FVC of sunny slopes was larger than that of shady slopes. The unfrozen and seasonal frozen soil regions on the QTP are distributed mainly at elevations of ≤5000 and ≤6000 m, respectively, whereas the permafrost soil is distributed at all elevation levels. FVC varied obviously with elevation for the different types of soil ( Figure 6). In unfrozen soil areas, FVC was largest at elevation ≤5000 m; the maximum value (0.929) was at elevation ≤2000 m. As elevation increased, FVC decreased gradually. In seasonal frozen soil areas, FVC was highest (0.812) at elevation ≤2000 m. With increasing elevation, the variation of FVC showed a trend of decline, followed by increase, and then decline. In permafrost regions, FVC mainly decreased with increasing elevation. At all elevations except 2000-3000 m, FVC values of permafrost regions were smaller than those of seasonal frozen soil regions. At elevation >6000 m, there was only permafrost and FVC was relatively low (0.038). The variation of FVC with different elevation and aspect is presented in Table 1. The FVC of both permafrost soil and unfrozen soil exhibited a downward trend with increasing elevation. The FVC of seasonal frozen soil generally decreased gradually with increasing elevation but with a slight increase at elevations of 3000-4000 m. In seasonal frozen soil, when elevation was ≤4000 m, FVC of sunny slopes was larger than that of shady slopes, and FVC of semi-sunny slopes was larger than that of semi-shady slopes. At elevation >4000 m, FVC of shady slopes was larger than that of sunny slopes, and FVC of semi-shady slopes was larger than that of semi-sunny slopes. In permafrost regions, FVC of semi-shady slopes was larger than that of semi-sunny slopes. Except at elevations of 3000-4000 m, FVC of shady slopes was larger than that of sunny slopes. In unfrozen soil areas, at elevations of 2000-4000 m, FVC of shady slopes was larger than that of sunny slopes, and at elevations ≤2000 or ≥4000 m, FVC of sunny slopes was larger than that of shady slopes.    Changes of CV with elevation on the QTP were analyzed (Table 2). In seasonal frozen soil areas, FVC of sunny slopes was larger than that of shady slopes at elevation ≤4000 m, and CV was largest (32.55%) at elevations of 2000-3000 m; thus, aspect had obvious impact on FVC. At elevation >4000 m, FVC of shady slopes was larger than that of sunny slopes, and CV was negative. At elevations of 5000-6000 m, CV was smallest (−29.32%), and it was obviously affected by aspect. In permafrost regions, except at elevations of 3000-4000 m, FVC of shady slopes was larger than that of sunny slopes, indicating that aspect had marked influence on FVC at elevation >4000 m. The influence of aspect on FVC in unfrozen soil areas was small when elevation was ≤4000 m; CV was maximum (27.17%) at elevations of 4000-5000 m. FVC of sunny slopes was greater than that of shady slopes at elevation ≤2000 m, and FVC of shady slopes was greater than that of sunny slopes at elevations of 2000-4000 m.

Effect of Slope and Aspect
On the QTP, FVC exhibited notable differences with slope ( Figure 7). The areal proportion of low coverage decreased gradually with increasing slope. The areal proportion of medium coverage decreased gradually with increasing slope except for slopes of for 2 • -6 • . The areal proportions of high and extremely high coverage increased gradually with increasing slope. Changes of CV with elevation on the QTP were analyzed (Table 2). In seasonal frozen soil areas, FVC of sunny slopes was larger than that of shady slopes at elevation ≤4000 m, and CV was largest (32.55%) at elevations of 2000-3000 m; thus, aspect had obvious impact on FVC. At elevation >4000 m, FVC of shady slopes was larger than that of sunny slopes, and CV was negative. At elevations of 5000-6000 m, CV was smallest (−29.32%), and it was obviously affected by aspect. In permafrost regions, except at elevations of 3000-4000 m, FVC of shady slopes was larger than that of sunny slopes, indicating that aspect had marked influence on FVC at elevation >4000 m. The influence of aspect on FVC in unfrozen soil areas was small when elevation was ≤4000 m; CV was maximum (27.17%) at elevations of 4000-5000 m. FVC of sunny slopes was greater than that of shady slopes at elevation ≤2000 m, and FVC of shady slopes was greater than that of sunny slopes at elevations of 2000-4000 m.

Effect of Slope and Aspect
On the QTP, FVC exhibited notable differences with slope ( Figure 7). The areal proportion of low coverage decreased gradually with increasing slope. The areal proportion of medium coverage decreased gradually with increasing slope except for slopes of for 2°-6°. The areal proportions of high and extremely high coverage increased gradually with increasing slope.  had extremely high coverage for all slope levels, although it decreased gradually with increasing slope. Variations of FVC with different slopes and aspect were analyzed (Table 3). Results showed that FVC exhibited little difference in relation to slope and aspect for each soil type. For seasonal frozen soil, when slope was ≤15°, FVC exhibited little difference in relation to slope and aspect; however, when slope was >15°, the slope and aspect effect became more obvious with increasing slope. For slope >25°, FVC values of shady and semishady slopes were larger than those of sunny and semi-sunny slopes. For permafrost regions, when slope was ≤2°, the slope and aspect effect was not obvious. When slope was 2°-6°, the FVC value of semi-shady slopes was larger than middle coverage, whereas it was low coverage for other slopes and aspects. When slope was >6°, the FVC values of shady and semi-shady slopes were larger than those of sunny and semi-sunny slopes. For unfrozen soil, when slope was ≤2°, the FVC values of semi-shady and semi-sunny slopes were relatively large, while those of shady and sunny slopes were relatively small. When slope was 2°-6°, the FVC values of sunny and shady slopes were greater than those of semi-sunny and semi-shady slopes. When slope was 6°-15°, the effect of aspect was substantial and for sunny slopes, FVC was larger than that for other aspects; when slope was >15°, the effect of aspect was relatively weak.  Variations of FVC with different slopes and aspect were analyzed (Table 3). Results showed that FVC exhibited little difference in relation to slope and aspect for each soil type. For seasonal frozen soil, when slope was ≤15 • , FVC exhibited little difference in relation to slope and aspect; however, when slope was >15 • , the slope and aspect effect became more obvious with increasing slope. For slope >25 • , FVC values of shady and semi-shady slopes were larger than those of sunny and semi-sunny slopes. For permafrost regions, when slope was ≤2 • , the slope and aspect effect was not obvious. When slope was 2 • -6 • , the FVC value of semi-shady slopes was larger than middle coverage, whereas it was low coverage for other slopes and aspects. When slope was >6 • , the FVC values of shady and semi-shady slopes were larger than those of sunny and semi-sunny slopes. For unfrozen soil, when slope was ≤2 • , the FVC values of semi-shady and semi-sunny slopes were relatively large, while those of shady and sunny slopes were relatively small. When slope was 2 • -6 • , the FVC values of sunny and shady slopes were greater than those of semi-sunny and semi-shady slopes. When slope was 6 • -15 • , the effect of aspect was substantial and for sunny slopes, FVC was larger than that for other aspects; when slope was >15 • , the effect of aspect was relatively weak. Obvious differences were found in the variations of CV with different slopes for the different soil types (Table 4). For seasonal frozen soil, when slope was ≤6 • , FVC of sunny slopes was larger than that of shady slopes, and when slope was >6 • , FVC of shady slopes was larger than that of sunny slopes. For permafrost, CV was negative for all slopes; for sunny slopes, FVC was less than that for shady slopes, and CV increased gradually with increasing slope. For unfrozen soil, CV was positive for all slopes; for sunny slopes, FVC was larger than that for shady slopes.

Discussion
Within the context of global climate change, temperature and precipitation are considered the main climatic factors leading to vegetation change. Additionally, human population, agriculture, and animal husbandry production activities are also important factors driving vegetation change [46]. The interaction between these natural and human factors has had substantial impact on the surface processes of the QTP, which has been reflected in changes to regional vegetation growth.

Natural Factors
Temperature and precipitation are important indicators of climate change and represent the most important climatic factors affecting plant distribution and vegetation growth on the QTP. In this study, temperature and precipitation data were spatially interpolated onto grids, and correlation analysis was performed between FVC and the two climatic factors (p < 0.05) (Figure 9). Generally, FVC had positive correlation with temperature in the range of 0-0.868 (r = 0.009) over most of the QTP. Areas with positive correlation were distributed mainly in Qinghai Province, northwestern Gansu Province, and western Sichuan Province. The area with negative correlation in the range of −0.860 to 0 was relatively small and distributed mainly in central Qinghai Province, southern Xinjiang Province, and southwestern parts of Tibet with high elevation (>4000 m). In the alpine environment, the rise in temperature was conducive to vegetation growth. However, with increasing elevation, the marked decrease in temperature became an important limiting factor of plant growth.  When the temperature in the frozen soil area is higher than 0 • C, the ice in the frozen soil starts to melt, which can provide certain water for the soil and vegetation. Vegetation on the QTP might be sensitive to temperature, especially the accumulated temperature over 0 • C. Relationship between accumulated temperature over 0 • C and FVC was analyzed, and the result showed FVC had negative correlation with accumulated temperature in the range of −1-0 (r = −0.018) on the most of the QTP. Areas with negative correlation were distributed mainly in southwestern parts of Qinghai Province, southern Xinjiang Province, southern and southwestern parts of Tibet, Yunnan and Sichuan Provinces. The area with positive correlation in the range of 0 to 1 was relatively small and distributed mainly in the central and eastern parts of Qinghai Province, northern part of Tibet and Gansu Provinces.
Over most of the QTP, positive correlation in the range of 0-0.937 was found to exist between FVC and precipitation (r = 0.089). Areas of negative correlation (−0.894 to 0) were distributed in small parts of the QTP. The regions with high correlation were distributed mainly in southern and northwestern Gansu Province, northern and central Qinghai Province, southern Xinjiang Province, and southwestern Tibet. Ample amounts of precipitation encourage vigorous growth of local vegetation; thus, precipitation is an important factor affecting the growth status of local vegetation on the QTP. The areas with negative correlation were distributed mainly in southeastern parts of the QTP, where precipitation showed a trend of decrease, which inhibited plant growth. Owing to the influence of the monsoons and westerly winds, projections indicate that annual precipitation in northeastern, central, and southwestern parts of the QTP will increase in 2020, 2050 and 2080, whereas annual precipitation in southeastern parts will decrease [47]. The precipitation distribution will further affect vegetation growth on the QTP.
During 2000-2020, the climate of the QTP was mainly warm and humid ( Figure 10). Temperature rise could meet the heat conditions required for plant growth in alpine areas. Additionally, it could increase the melting degree of frozen soil and thereby increase water availability for vegetation. The effects of temperature rise in areas with different topography and vegetation types would be notably different. The increase in temperature would stimulate the growth of alpine meadow, but temperature-induced drought would inhibit such growth [48]. Moreover, thawing lakes and thawing collapses attributable to thawing of frozen soil would also degrade or eradicate some vegetation. Overall, increased temperature and enhanced precipitation favor vegetation growth. However, on the QTP, their combined effects exhibited strong spatial heterogeneity. An increase in precipitation tended to enhance vegetation when the temperature tended to increase, but when precipitation decreased, temperature was the limiting factor [49]. The warm and wet areas of the QTP accounted for 32.78% of the total area; cold and wet areas accounted for 32.59%; warm and dry areas accounted for 19.07%, and cold and dry areas accounted for 15.56%. The tendency for warm and wet conditions was the main reason for the improvement of vegetation on the QTP. For different frozen soil types, the effects of temperature and precipitation change were very different. In seasonal frozen soil areas, the area with warm and wet conditions accounted for the largest proportion of the total area (36.40%), followed by areas with cold and wet conditions (29.98%), warm and dry conditions (17.82%), and cold and dry conditions (15.8%). In the permafrost region, the area with warm and wet conditions accounted for the largest proportion of the total area (36.93%), followed by areas with cold and wet conditions (27.73%), warm and dry conditions (20.6%), and cold and dry conditions (14.74%). The area with warm and wet conditions in both seasonal frozen soil areas and permafrost soil regions was largest, which played an important role in promoting the improvement of local vegetation growth. In unfrozen soil regions, the area with cold and dry conditions was largest (31.59%), followed by areas with warm and dry conditions (29.36%), areas with cold and wet conditions (20.62%), and areas with warm and wet conditions (18.43%). The area with warm and wet conditions in unfrozen soil regions accounted for the smallest proportion, whereas the areas with cold and dry conditions and warm and dry conditions accounted for 60.95% of the total. The drying trend of the environment showed a certain inhibitory effect on the local vegetation.

Human Factors
Human factors are also an important driving force affecting vegetation gro the QTP. The surge in population and the large amount of greenhouse gas emissio had obvious consequences regarding global warming [50]. According to the sixth of the Intergovernmental Panel on Climate Change issued in 2021, global temperat has reached 1.1 °C (above the preindustrial benchmark temperature) and is in the of reaching the proposed upper limit (1.5 °C) for sustainable human habitability [ der this background of global warming, degradation of glaciers and frozen soi QTP has changed the surface water cycle process, and affected the distributi growth of vegetation. Following recent rapid growth of the population on the Q impact of human activities on the ecological environment has increased, which h affected local vegetation growth. Human activities could lead to increase in local tion coverage through the Grain-to-Green Program, afforestation, and promoting tural technological progress. However, damage caused to the ecological environm unsustainable farming practices, urban expansion, overgrazing, and mining cou lead to substantial deterioration in regional vegetation coverage. Taking Qinghai P as an example, following implementation of China's western development prog

Human Factors
Human factors are also an important driving force affecting vegetation growth on the QTP. The surge in population and the large amount of greenhouse gas emissions have had obvious consequences regarding global warming [50]. According to the sixth report of the Intergovernmental Panel on Climate Change issued in 2021, global temperature rise has reached 1.1 • C (above the preindustrial benchmark temperature) and is in the process of reaching the proposed upper limit (1.5 • C) for sustainable human habitability [51]. Under this background of global warming, degradation of glaciers and frozen soil on the QTP has changed the surface water cycle process, and affected the distribution and growth of vegetation. Following recent rapid growth of the population on the QTP, the impact of human activities on the ecological environment has increased, which has also affected local vegetation growth. Human activities could lead to increase in local vegetation coverage through the Grain-to-Green Program, afforestation, and promoting agricultural technological progress. However, damage caused to the ecological environment by unsustainable farming practices, urban expansion, overgrazing, and mining could also lead to substantial deterioration in regional vegetation coverage. Taking Qinghai Province as  .04 million ha, accounting for 26.50% of the entire area of the province. The land areas with desertification and obvious desertification trend accounted for 12.46 million and 4.13 million ha, respectively. Desertification on the QTP remains a serious problem that has major impact on local vegetation coverage.
Although increases in the human population and livestock numbers might cause substantial pressure on regional vegetation, implementation of active policies could also provide favorable conditions for vegetation restoration. Programs for relocation of herdsmen and adjustment of the production structure in ecologically fragile areas have led to continuous reduction in the proportion of the population involved in agriculture and animal husbandry, and further reduced the grassland pressure in agricultural and pastoral areas of the QTP. Additionally, implementation of a series of ecological engineering policies such as the Construction of Ecological Security Barriers has curbed the trend of vegetation degradation or even led to local improvement [52]. In terms of reserve construction, Qinghai Province has established 11 nature reserves (http://www.mee.gov.cn/, accessed on 9 October 2021). All these nature reserves have played a very positive role in vegetation restoration in Qinghai Province. Moreover, Qinghai Province began to implement the Grainto-Green Program in 2000, and it had completed afforestation of an area of 8.91 × 10 5 ha by 2011. By the end of 2016, the area of farmland returned to forests was 7.40 × 10 5 ha, and the total forested area had reached 4.52 × 10 7 ha. In terms of desertification control projects, Qinghai Province started the "Three-North Shelterbelt" Project in 1978. Since 2016, Qinghai Province has implemented the Qinghai Desertification Prevention and Control Plan (2011-2020), and it had completed desertification land control tasks over an area of 4.24 × 10 5 ha by the end of 2019. According to the results of the Fifth Desertification Land Monitoring in Qinghai Province, the areas of desertification land and sand land were reduced by 5.1 × 10 4 and 5.7 × 10 4 ha, respectively, between 2009 and 2014. The ecological environment of the sand land area has improved substantially, and an initial green ecological barrier has been established. The active environmental protection policies outlined above have played a very important role in vegetation restoration and improved vegetation growth on the QTP.

Conclusions
Using MODIS-EVI and DEM data, this study investigated the spatiotemporal variation characteristics of vegetation growth status in areas of the QTP with different soil types during 2000-2020 and analyzed FVC variations attributable to different topographic conditions. The principal conclusion was that the variation trend of FVC showed a gradual trend of growth. The multiyear mean value of FVC was 0.421, and the vegetation index gradually decreased from the southeast to the northwest. The FVC value in the permafrost soil area was the smallest, followed by the seasonal frozen soil area, and the unfrozen soil area. Topographic effects varied for all soil types depending on elevation and slope. Correlation analysis between vegetation coverage and temperature or precipitation revealed positive correlation between FVC and precipitation or temperature. Analysis of temperature and humidity revealed that the QTP climate has tended to become warmer and wetter, which is the principal natural reason for the improvement of FVC on the QTP. Although increases in the human population and numbers of livestock might cause considerable pressure on local vegetation, implementation of active policies such as the Grain-to-Green Program, afforestation, closing mountains for afforestation, and promoting agricultural technological progress has increased vegetation coverage.