Effects of Climate Variability and Human Activities on Vegetation Dynamics across the Qinghai–Tibet Plateau from 1982 to 2020

.


Introduction
Due to increasing greenhouse gas concentrations, global warming, an incontrovertible reality, has resulted in alterations of precipitation patterns and increased temperatures, affecting ecosystems both directly and indirectly [1].Vegetation plays a pivotal role in terrestrial ecosystems by linking the biosphere, pedosphere, atmosphere, lithosphere, and hydrosphere [2][3][4].The QTP, recognized as the "Roof of the World" and the "Third Pole", holds the highest elevation among global plateaus and is often referred to as the "Asian Water Tower".It covers approximately 70% of the total plateau area and supports the most extensive alpine grassland ecosystem globally.Owing to its unique geography and complex environment, the QTP is one of the primary ecologically vulnerable regions in the world, serving as a responsive and sensitive zone to climate change, influencing not only Asia but also the broader Northern Hemisphere [5,6].
Over the past 50 years, warming in the QTP has doubled the global warming rate, particularly in winter.Hence, the QTP is becoming one of the foremost indicators of significant global warming [5][6][7].Within the global warming trend, temperature escalations contribute to glacier shrinkage, snowpack melting, and a concomitant reduction in the total area of the glaciers and permafrost zones.This intricate interplay can trigger a series of unpredictable socioenvironmental and ecological crises [4].These factors directly influence the hydrological patterns of the QTP, notably causing divergences in precipitation between the northern and southern regions.The enhanced variability in precipitation is an evident aspect of climate change on this plateau [8].Overall, climate change in the QTP is characterized by the intertwined effects of warming and increased moisture despite regional variations [9].This climate profile indicates a complex interplay of factors, necessitating a detailed examination and understanding of the interconnected factors that determine the ecological future of the QTP.
By 2015, the QTP, a significant national grazing zone, saw its actual livestock capacity rise to 158 million sheep units, surpassing the theoretical capacity by over 1.6 times.Overgrazing, especially in the northwest QTP, has increased since the 1990s and peaked in the 2000s [10].Moreover, the QTP also sustains approximately 5.3 million pastoralists who had largely relied on grass-fed livestock farming by 2015 [11].With population growth and infrastructure development, socioeconomic advancements have intensified pressures on ecosystems [6,12,13], mainly in the western and eastern areas of the QTP.The latest data indicate that alpine arid desert and sparse vegetation account for 34.9% of the QTP, exhibiting severe desertification, soil erosion, and freeze-thaw erosion.In addition, forest and shrub degradation reached 59%, primarily in the Hengduan Valley, while grassland degradation reached 80%, mainly in the northwest of the QTP [5,14].There are also adverse effects, such as soil texture changes, nutrient depletion, and shifts in soil communities [15,16].These changes can be attributed to climate change and human activities.
Grassland degradation on the Qinghai-Tibet Plateau has become an urgent ecological issue.The driving forces behind this can be subject to various perspectives in the academic community.Although several studies suggest that climatic factors play a dominant role [17], others believe that human activities exert a more pronounced influence than climate change [18,19].Therefore, examining the roles of climate change and human activities in vegetation dynamics is essential for understanding the implications of climate warming and for offering insights into the ecological sustainability of the QTP.
Currently, remote sensing data with a high spatial-temporal resolution are extensively employed for monitoring surface vegetation dynamics and biomass assessment [20,21].The NDVI and Enhanced Vegetation Index (EVI) are commonly applied indices for assessing the vegetation status in the QTP.Shen et al. [9] utilized three satellite-derived NDVI datasets (AVHRR NDVI, SPOT NDVI, and MODIS NDVI) to investigate temporal changes in QTP vegetation.The AVHRR NDVI data indicated a greening trend from 1982 to 2010, while all datasets highlighted a decline in the growing season NDVI in the southwest of the plateau during 2000-2010.In contrast, greening persisted in the northeastern region.Based on the GIMMS NDVI, Wei et al. [19] identified temperature rises and decreases in the high-altitude austere climate as factors expanding vegetation and improving high-altitude vegetation, particularly in the southwest, albeit with degradation in the northern and eastern QTP.Wang et al. [22] observed a general NDVI increase since the 1980s, albeit with notable spatial variability and localized decreases.
Many scholars have explored the relationship between vegetation and climate on the QTP and found that climate variability has a direct impact on the growth of vegetation [23].Li et al. [13] utilized MODIS NDVI data to observe that vegetation in the northeastern and southwestern regions of the QTP responded prominently to precipitation, while the southern and central-eastern regions were more responsive to temperature.Utilizing the random forest algorithm, Wang et al. [22] assessed vegetation changes in the QTP and identified precipitation variability and the DEM (Digital Elevation Model) as primary drivers for both greening and browning across various vegetation types.Certain studies have indicated that, in addition to largely driving NDVI growth, rising temperatures can also intensify water shortages in low-precipitation areas of the QTP [24].In addition, Lou et al. [25] employed a novel deep learning method and identified precipitation and temperature as the principal factors behind the greening of the QTP.
Numerous researchers have explored the primary drivers affecting vegetation growth in different subregions of the QTP.Bai et al. [26] discovered a modest upward trend in the NDVI in Sanjiangyuan from 1982 to 2015.Although the annual average temperature significantly increased, the precipitation remained largely unchanged.They concluded that the temperature rise facilitated vegetation growth and recovery in this area.Generally, on the QTP, heightened precipitation during warmer periods favors alpine vegetation development.Conversely, in arid periods temperature acts as the dominant constraint on vegetation growth [27].These results indicate that, at the regional level, the influence of human activities on vegetation changes may not surpass that of climate change [13,17].
Numerous scholars have explored vegetation changes and their driving mechanisms across various time periods and study areas within the QTP, using diverse satellite datasets [18,19,28,29].Based on the research above, we have used the NOAA CDR AVHRR NDVI long-term time series data to assess vegetation changes on the QTP from 1982 onwards, covering over four decades.In addition, scholars have mostly used methods such as machine learning and remote sensing model simulation to explore the factors affecting vegetation change.In this paper, utilizing a residual-based method, we evaluated the effects of climate change and human activity on vegetation changes for nearly 40 years.Our approach identified the impacts of climate change and human activities on NDVI trends over this extended period, capturing the relative contributions of climatic and anthropogenic factors to vegetation changes.Such differentiation is crucial for providing a scientific foundation for formulating effective vegetation restoration strategies for the QTP.The main objectives of this research were as follows: (1) to examine both temporal and spatial variations in vegetation across the plateau and its specific subregions; (2) to quantify the effects of climate change (CC) and human activities (HAs) on vegetation and to determine their individual contributions; and (3) to investigate the influencing mechanisms of meteorological elements, such as temperature increases and precipitation shifts, on vegetation, as well as the implications of contemporary human activities, including livestock counts and population fluctuations, for these dynamics.

Study Area
Located between 26 • -39 • N and 73 • -104 • E (Figure 1b), the QTP has an average elevation of 4400 m.It extends from the Pamir Plateau in the west to the Hengduan Mountains in the east, bordered by the Kunlun, Alps, and Qilian Mountains to the north and the Himalayas to the south.Measuring 2800 km in length from east to west, with a width varying between 300 and 1500 km, the QTP covers an area of roughly 2.57 million km 2 .This region includes, either entirely or in part, six provinces: Tibet, Qinghai, Sichuan, Yunnan, Xinjiang, and Gansu [30].

Data and Materials
We utilized the AVHRR NDVI 5th edition from the NOAA CDR as the remote sensing NDVI data source, featuring gridded daily NDVI measurements from 1981 to 10 days before the present using NOAA satellite data (7, 9, 11, 14, 16, 17, 18, and 19) and presented In the northwestern regions of the QTP, an alpine desert steppe dominated by an arid climate is characterized by expanses of the Gobi Desert and visible rock formations (Figure 1a).Isolated pockets of alpine swamp meadows punctuate the landscape and cluster around water bodies, such as lakes and rivers.The vegetation of the QTP has shifted from evergreen broad-leaf forests in the southeast to desert grasslands in the northwest (Figure 1a), representing one of the most diverse vertical mountainous landscapes in China.This unique ecological structure indicates that the QTP plays a significant role as a critical carbon sink [31].Agricultural activities in the QTP are tailored to the challenging environment.Due to the harsh climate and elevated altitude, pastoralism and animal husbandry dominate.Agricultural practices largely focus on the subsistence and cultivation of cold-resistant crops, such as barley and potatoes.Populations on the QTP cluster in valleys and river basins, which offer relatively mild conditions.Key administrative and economic centers include Lhasa and Xining, with smaller settlements dispersed throughout the plateau.

Data and Materials
We utilized the AVHRR NDVI 5th edition from the NOAA CDR as the remote sensing NDVI data source, featuring gridded daily NDVI measurements from 1981 to 10 days before the present using NOAA satellite data (7, 9, 11, 14, 16, 17, 18, and 19) and presented on a 0.05 • × 0.05 • global grid [32].These data used AVHRR bands 1 and 2 for NDVI calculation (From Table 1).This high-quality, long-term dataset facilitates global vegetation growth studies and is applicable to areas such as arid region monitoring, agricultural forecasting, crop management, and biology.This dataset is one of the Land Surface CDR Version 5 products produced by the NASA Goddard Space Flight Center (GSFC) and the University of Maryland (UMD).Improvements for Version 5 include using the improved surface reflectance data; correcting the data for known errors in time, latitude, and longitude variables; and improvements in the global and variable attribute definitions.The dataset is in the netCDF-4 file format following ACDD and CF conventions.The dataset is accompanied by algorithm documentation, a data flow diagram, and the source code for the NOAA CDR Program.To complement the AVHRR CDR NDVI for QTP research, we selected the Chinese meteorological forcing dataset (CMFD) spanning from 1979 to 2018, sourced from the National Tibetan Plateau Data Center [33].This dataset integrates various observational data, including remote sensing products, meteorological station records, and reanalysis data, and is primarily designed to support surface change research in China.The dataset captured seven near-surface meteorological elements from January 1979 to December 2018, including the 2 m air temperature, surface pressure, specific humidity, 10 m wind speed, downward short-wave and long-wave radiation, and precipitation rate.It has a 3 h temporal resolution and a spatial resolution of 0.1 To harmonize our survey with the focus of the study, DEM data featuring a resolution of 90 × 90 m were purchased from NASA (Figure 1b).Vegetation type data were derived from ESA's Climate Change Initiative Land Cover dataset V2015 at a resolution of 300 × 300 m (Figure 1a).Socioeconomic data, including year-end livestock inventories, the Gross Domestic Product (GDP), and population figures, were primarily extracted from recent statistical yearbooks of Qinghai Province, the Tibet Autonomous Region, and other relevant provinces.These data were subsequently restructured at the county level for analysis.The eco-climatic zones followed the framework (as shown in Table 2) downloaded from Resource and Environmental Sciences Data Centre (https://www.resdc.cn).

Methods (1) Maximum Value Composition Method
Firstly, there is data preprocessing.Because the NDVI data are daily gridded NOAA CDR AVHRR NDVI data, we applied the maximum value composition method to transform daily data into a monthly NDVI to mitigate certain effects, such as the solar angle, cloud interference, and atmospheric residuals.This approach has been widely employed in NDVI product preprocessing [19].The calculation procedure is as follows: where NDV I i is the daily NDVI and NDV I j is the monthly NDVI.Through the above equation, we aggregated the daily NDVI data into monthly NDVI datasets.Subsequently, based on these monthly NDVI datasets, we synthesized four-season NDVI and annual NDVI datasets, from which the annual NDVImax is obtained for all the NDVIs involved in the following.
(2) Trend Analysis In this study, we utilized the Theil-Sen Median method to assess temporal patterns of regional factors, including the annual NDVI, mean annual temperature, and annual precipitation at the pixel level.The Theil-Sen Median method, commonly referred to as Sen's slope estimation, is adept at identifying trends in time-series data.Notably, it can handle datasets with missing information, offers computational efficiency, and remains robust to measurement inaccuracies and outliers.Its widespread application in prolonged time-series trend analysis is well recognized.The calculation formula is as follows: where x i and x j denote the time-series data.If β > 0, the NDVI exhibits an upward trend; otherwise, it shows a downward trend.To validate this trend, the Mann-Kendall method was employed [34], which has the advantages of independence from assumptions about a normal distribution or linear trends and resilience against missing values and outliers.This nonparametric approach has been extensively utilized to analyze long-term time-series data trends.With this treatment, we obtained the spatial trends in the NDVI, temperature, precipitation, and other elements on the QTP.
(3) Residual Analysis Residual analysis has been widely employed to isolate the effects of both human activities and climate variables on vegetation change [35].Typically, multiple linear regressions can establish relationships between the NDVI and climate variables, such as the mean temperature, precipitation, and solar radiation over a consistent period.This facilitates the estimation of the impacts of climate variables on vegetation, indirectly highlighting the crude impacts of human activities on the NDVI [36].In this study, we delineated the linear relationship between the NDVI and monthly climate variables as follows: where climate elements include the Temp, representing the average annual temperature; Pre, representing the annual precipitation; and Solar, representing downward short-wave radiation.The NDVI is the real NDVI value, NDV I cc denotes the effect of climate factors on vegetation, and NDV I ha denotes the residual, which is a proxy for the impacts of human activities on vegetation (including pastoralism, population growth, and road infrastructure development).In this paper, we use the residual method to separate the effects of climatic elements and human activities on NDVI dynamics.In general, the annual NDVI of vegetation is mainly determined by climatic conditions and the intensity of human activities.Consequently, we disentangle the contribution of climate change to the NDVI by establishing a regression model based on the relationship between climate conditions and the NDVI, predicting the annual NDVI determined by climatic conditions.Under the assumption of there being no other nondeterministic factors, the residual variation between the observed NDVI values and the predicted values based on climate change represents the portion attributed to human activities.In the absence of a human activity influence, the interannual variation of residuals should exhibit random characteristics centered around zero.If there is a significant decreasing trend in the interannual variation of residuals, it indicates vegetation degradation induced by human activities.Conversely, it suggests that human activities have improved the ecological environment.Subsequently, a trend analysis was conducted on the NDV I ha and NDV I cc over an extended time series.F NDV I denotes the linear slope of the NDVI, F CC denotes the linear slope of the NDVI under the impacts of climate variables, and F H A denotes the linear slope of the NDVI under the influence of human activities.Furthermore, the relative contributions of climate change and human activities were calculated using the equations shown in Table 3 [37,38].
The Pearson correlation coefficient r is a standard metric used to measure the correlation between two variables.Therefore, to assess the response of vegetation to climatic and anthropogenic factors, correlation analysis was conducted.The correlation coefficients between the NDVI and meteorological data were determined as follows: where x i and y i represent the NDVI and meteorological factor at the pixel scale during 1982-2020, respectively, and r xy denotes the correlation coefficient between indicators x i and y i .In this paper, we calculate the spatially distributed correlations of the NDVI with temperature and NDVI with precipitation, respectively.In addition, in the data validation section we also calculate the spatial correlation between NOAA CDR AVHRR NDVI and MODIS NDVI data and between NOAA CDR AVHRR NDVI and GIMMS-3gNDVI data to validate their usability.

Data Quality Assessment
Owing to missing NOAA CDR AVHRR NDVI raw data for October, November, and December in 1994 and 2020 in China, a value of 0 was substituted during data processing.To ensure reliability, the acquired NOAA CDR AVHRR NDVI data were assessed for validation.During 1982-2015, the correlation coefficient between the NOAA CDR AVHRR NDVI and GIMMS3g NDVI was significant (r = 0.95, p < 0.05, n = 408).Given that both datasets rely on the same satellite sensor, a high correlation is anticipated.During 2001-2020, the correlation between the NOAA CDR AVHRR NDVI and MODIS NDVI was marginally lower (r = 0.88, p < 0.05, n = 240).During 2001-2015, the correlation coefficient between the MODIS NDVI and GIMMS3g NDVI was notably high (r = 0.93, p < 0.05, n = 180).Over the period 2001-2015, a high correlation coefficient was observed between the NOAA CDR AVHRR NDVI and MODIS NDVI (r = 0.93, p < 0.05, n = 180).For the period 2001-2015, the consistency between the NOAA CDR AVHRR NDVI and MODIS NDVI surpassed that of the GIMMS3g NDVI.Similarly, between 2001 and 2020 the consistency between the NOAA CDR AVHRR NDVI and MODIS NDVI was higher than that of the GIMMS3g NDVI.The reduced consistency of the MODIS NDVI can be attributed to the absence of NOAA CDR AVHRR NDVI data during 2015-2020, particularly the missing values in 2020.
To provide additional validation of the above process, we also calculated the spatial correlation between the individual NDVI data (Figure 2a,b).The NOAA CDR AVHRR NDVI exhibits a robust spatial correlation with the other two NDVI datasets across the entire QTP.Although slight variations are observed in different subregions, they do not significantly impact our ability to dynamically monitor vegetation changes across the entire plateau.Hence, the NOAA CDR AVHRR NDVI data, with their compatibility with other NDVI datasets in the QTP and extensive time series, can be well suited for subsequent analytical endeavors.To provide additional validation of the above process, we also calculated the spatial correlation between the individual NDVI data (Figure 2a,b).The NOAA CDR AVHRR NDVI exhibits a robust spatial correlation with the other two NDVI datasets across the entire QTP.Although slight variations are observed in different subregions, they do not significantly impact our ability to dynamically monitor vegetation changes across the entire plateau.Hence, the NOAA CDR AVHRR NDVI data, with their compatibility with other NDVI datasets in the QTP and extensive time series, can be well suited for subsequent analytical endeavors.

Time Series of Major Change
As shown in Figure 3a, the NDVI over the QTP has exhibited a significant increase since 1982 at a rate of 0.0013/a.By employing different monitoring data sources and diverse vegetation indices (e.g., the Leaf Area Index, Net Primary Production), numerous researchers have similarly reported increasing vegetation trends on the QTP, with NDVI rates ranging between 0.002 and 0.0020/a [39,40].We divided the period into three distinct phases: an increase of 0.0012/a from to 1982 to 2000 and a decline of 0.0010/a between 2000 and 2015, followed by a pronounced increase at 0.0130/a thereafter.In the QTP, alpine grassland and alpine meadows primarily dominate the vegetation distribution, typically exhibiting low NDVI values.Therefore, NDVI values below 0.10 were categorized as bare ground [35].Figure 3b illustrates that winter and spring NDVI values have remained relatively stable in recent years at approximately 0.2.This stability

Time Series of Major Change
As shown in Figure 3a, the NDVI over the QTP has exhibited a significant increase since 1982 at a rate of 0.0013/a.By employing different monitoring data sources and diverse vegetation indices (e.g., the Leaf Area Index, Net Primary Production), numerous researchers have similarly reported increasing vegetation trends on the QTP, with NDVI rates ranging between 0.002 and 0.0020/a [39,40].We divided the period into three distinct phases: an increase of 0.0012/a from to 1982 to 2000 and a decline of 0.0010/a between 2000 and 2015, followed by a pronounced increase at 0.0130/a thereafter.To provide additional validation of the above process, we also calculated the spatial correlation between the individual NDVI data (Figure 2a,b).The NOAA CDR AVHRR NDVI exhibits a robust spatial correlation with the other two NDVI datasets across the entire QTP.Although slight variations are observed in different subregions, they do not significantly impact our ability to dynamically monitor vegetation changes across the entire plateau.Hence, the NOAA CDR AVHRR NDVI data, with their compatibility with other NDVI datasets in the QTP and extensive time series, can be well suited for subsequent analytical endeavors.

Time Series of Major Change
As shown in Figure 3a, the NDVI over the QTP has exhibited a significant increase since 1982 at a rate of 0.0013/a.By employing different monitoring data sources and diverse vegetation indices (e.g., the Leaf Area Index, Net Primary Production), numerous researchers have similarly reported increasing vegetation trends on the QTP, with NDVI rates ranging between 0.002 and 0.0020/a [39,40].We divided the period into three distinct phases: an increase of 0.0012/a from to 1982 to 2000 and a decline of 0.0010/a between 2000 and 2015, followed by a pronounced increase at 0.0130/a thereafter.In the QTP, alpine grassland and alpine meadows primarily dominate the vegetation distribution, typically exhibiting low NDVI values.Therefore, NDVI values below 0.10 were categorized as bare ground [35].Figure 3b illustrates that winter and spring NDVI values have remained relatively stable in recent years at approximately 0.2.This stability In the QTP, alpine grassland and alpine meadows primarily dominate the vegetation distribution, typically exhibiting low NDVI values.Therefore, NDVI values below 0.10 were categorized as bare ground [35].Figure 3b  The QTP, with its high altitude; intricate topography; extensive latitudinal and longitudinal ranges; and diverse eco-climatic and vegetation zones, exhibits pronounced spatial heterogeneity in NDVI change (Figure 4).Notably, the NDVI progressively decreased from east to west owing to the particular vegetation distribution on the QTP.The alpine grasslands and meadows on the QTP were further influenced by global changes [42,43].Figure 4 illustrates that only 17.6% of the vegetation NDVI demonstrated a decline along the west-north and northern edge of the QTP from 1982 to 2020.Conversely, the regions east of 95 • E on the QTP demonstrate a more pronounced upward trend.The inset of the NDVI spatial distribution map indicates an NDVI trend increase of 0.00004-0.00005/aper 1 • longitude, while the adjacent graph on the right represents an NDVI trend increase of 0.00003-0.00005/aper 1 • latitude.These graphs suggest that an average NDVI trend exceeding 0.002/a is primarily situated south of 29 The QTP, with its high altitude; intricate topography; extensive latitudinal and longitudinal ranges; and diverse eco-climatic and vegetation zones, exhibits pronounced spatial heterogeneity in NDVI change (Figure 4).Notably, the NDVI progressively decreased from east to west owing to the particular vegetation distribution on the QTP.The alpine grasslands and meadows on the QTP were further influenced by global changes [42,43].Figure 4 illustrates that only 17.6% of the vegetation NDVI demonstrated a decline along the west-north and northern edge of the QTP from 1982 to 2020.Conversely, the regions east of 95°E on the QTP demonstrate a more pronounced upward trend.The inset of the NDVI spatial distribution map indicates an NDVI trend increase of 0.00004-0.00005/aper 1° longitude, while the adjacent graph on the right represents an NDVI trend increase of 0.00003-0.00005/aper 1° latitude.These graphs suggest that an average NDVI trend exceeding 0.002/a is primarily situated south of 29°N and east of 98°E.

Spatial Changes in Different Zones
We conducted an in-depth analysis of vegetation trends across distinct zones to elucidate the spatial heterogeneity of vegetation changes in the QTP.Our study evaluated various eco-climatic zones, as illustrated in Figure 5a.The zones with minimal changes mainly included the eco-climatic regions H1D1, H2D3, H2D1, and H2D2.Notably, these zones were positioned within the arid expanse and western periphery of the QTP and were characterized by sparse vegetation cover.As shown in Figure 5d, the distribution across the eco-climatic zones was relatively uniform.However, the humid zones, particularly VA5 and VA6, had the lowest proportion.Although the H1C2 zone, centrally located in this domain, exhibited a mean NDVI trend just below 0.001/a, regions indicating significant enhancement (with NDVI trends surpassing 0.002/a) were predominantly identified on the eastern edge and central area of the QTP.Zones such as H2C1, VA5, VA6, H2AB1, and H1B1 were situated in lower alpine temperate areas or moist to subhumid climate regions.

Spatial Changes in Different Zones
We conducted an in-depth analysis of vegetation trends across distinct zones to elucidate the spatial heterogeneity of vegetation changes in the QTP.Our study evaluated various eco-climatic zones, as illustrated in Figure 5a.The zones with minimal changes mainly included the eco-climatic regions H1D1, H2D3, H2D1, and H2D2.Notably, these zones were positioned within the arid expanse and western periphery of the QTP and were characterized by sparse vegetation cover.As shown in Figure 5d, the distribution across the eco-climatic zones was relatively uniform.However, the humid zones, particularly VA5 and VA6, had the lowest proportion.Although the H1C2 zone, centrally located in this domain, exhibited a mean NDVI trend just below 0.001/a, regions indicating significant enhancement (with NDVI trends surpassing 0.002/a) were predominantly identified on the eastern edge and central area of the QTP.Zones such as H2C1, VA5, VA6, H2AB1, and H1B1 were situated in lower alpine temperate areas or moist to subhumid climate regions.We examined NDVI dynamics across diverse land cover types, including 10 primary categories.The analysis utilized land use data primarily sourced from the land use data shown in Figure 1a, where the original 23 classifications were combined into 10 primary categories."Tree cover, flooded" and "other forest" collectively accounted for only 0.002%, and other vegetation types, such as wetland, settlement, and shrubland, each represent less than 1%.Consequently, NDVI trend analysis was not conducted for categories below 1%. Figure 5b distinctly highlights the increased growth in "agriculture," "forest," and "other forest" zones.Both cropland and tree cover demonstrated increases exceeding 0.002/a.In the second tier, grassland, wetland, and settlement exhibited increments of approximately 0.001/a, with shrubland and water reflecting modest increases of approximately 0.0008/a.The "water" category included surface features, such as glaciers and snow.The increases in snow and tree lines due to global warming have induced vegetated area expansion, thereby increasing the NDVI for aquatic regions.Conversely, bare land (including the desert and Gobi categories) lacking significant NDVI changes was consistent with the existing knowledge.
Figure 5c illustrates the NDVI trends per 500 m elevation increment.As shown in Figure 5f, the elevation profile of the QTP primarily laid between 4500 and 5500 m, representing a significant 70.3% plateau.With an increasing altitude, there was a noticeable decline in the NDVI.The QTP was divided into 15 subregions, with only four smaller zones demonstrating notable NDVI increases of approximately 0.002/a.These zones cover less than 1% of the entire QTP.Below 2500 m elevation, an NDVI increase of 0.002/a was evident.In contrast, above 6000 m the NDVI change was minimal at 0.0006/a.Notably, between 3500 and 4500 m there was a substantial NDVI increase of 0.0018/a.Compared with other elevational bands, the vegetation within the 3500-4500 m range demonstrates increased responsiveness to environmental changes, especially climatic variations.

Contribution of HA and CC to NDVI Change
As shown in Figure 6, we estimated the slope of human activities (HAs) and climate change (CC, encompassing Pre, Temp, and Srad) to the NDVI from 1982 to 2020 for each grid of the QTP, assessed using residual analysis.Figure 6a illustrates that HA positively influenced NDVI changes in 61.8% of the QTP.Specifically, the contribution of HA exceeded 0.0005/a in 30.5% of the regions and 0.001/a in 14.2%, predominantly in the southern and southeastern QTP.Conversely, the northwestern and northern parts of the QTP exhibited a negative contribution to HA.On average, HA contributed 0.0003/a to the We examined NDVI dynamics across diverse land cover types, including 10 primary categories.The analysis utilized land use data primarily sourced from the land use data shown in Figure 1a, where the original 23 classifications were combined into 10 primary categories."Tree cover, flooded" and "other forest" collectively accounted for only 0.002%, and other vegetation types, such as wetland, settlement, and shrubland, each represent less than 1%.Consequently, NDVI trend analysis was not conducted for categories below 1%. Figure 5b distinctly highlights the increased growth in "agriculture," "forest," and "other forest" zones.Both cropland and tree cover demonstrated increases exceeding 0.002/a.In the second tier, grassland, wetland, and settlement exhibited increments of approximately 0.001/a, with shrubland and water reflecting modest increases of approximately 0.0008/a.The "water" category included surface features, such as glaciers and snow.The increases in snow and tree lines due to global warming have induced vegetated area expansion, thereby increasing the NDVI for aquatic regions.Conversely, bare land (including the desert and Gobi categories) lacking significant NDVI changes was consistent with the existing knowledge.
Figure 5c illustrates the NDVI trends per 500 m elevation increment.As shown in Figure 5f, the elevation profile of the QTP primarily laid between 4500 and 5500 m, representing a significant 70.3% plateau.With an increasing altitude, there was a noticeable decline in the NDVI.The QTP was divided into 15 subregions, with only four smaller zones demonstrating notable NDVI increases of approximately 0.002/a.These zones cover less than 1% of the entire QTP.Below 2500 m elevation, an NDVI increase of 0.002/a was evident.In contrast, above 6000 m the NDVI change was minimal at 0.0006/a.Notably, between 3500 and 4500 m there was a substantial NDVI increase of 0.0018/a.Compared with other elevational bands, the vegetation within the 3500-4500 m range demonstrates increased responsiveness to environmental changes, especially climatic variations.

Contribution of HA and CC to NDVI Change
As shown in Figure 6, we estimated the slope of human activities (HAs) and climate change (CC, encompassing Pre, Temp, and Srad) to the NDVI from 1982 to 2020 for each grid of the QTP, assessed using residual analysis.Figure 6a illustrates that HA positively influenced NDVI changes in 61.8% of the QTP.Specifically, the contribution of HA exceeded 0.0005/a in 30.5% of the regions and 0.001/a in 14.2%, predominantly in the southern and southeastern QTP.Conversely, the northwestern and northern parts of the QTP exhibited a negative contribution to HA.On average, HA contributed 0.0003/a to the NDVI changes across the QTP, representing 41.2% of the entire region.The cumulative impact of the three climatic factors under CC was positive for 75.7% of the QTP, with most areas exceeding 0.0005/a.The mean CC contribution to the NDVI change of the QTP was 0.0008/a, accounting for 44.1% of the total NDVI change.As depicted in Figure 6c,d, we computed the relative contribution rates of CC and HA to NDVI changes.The total contributions of CC to NDVI were dominant (with relative contribution rates surpassing 50%) in 70.6% of the areas, mainly distributed in the central and northeastern regions of the QTP (Figure 6c).The dominant contribution rates of HA to NDVI change were observed in 29.4% of the QTP, primarily in the northwestern and southeastern regions (Figure 6d). Figure 7 illustrates the magnitude and relative contribution rate of each factor to the NDVI changes at the city scale.Most cities exhibited a CC contribution surpassing 0.001/a, with only three cities, Haixi, Rikaze, and Naqu, registering contributions below this value.Generally, HA contributions did not exceed 0.001/a, with the Ngari Prefecture notably presenting a negative HA contribution.Figure 7b shows that the contribution of CC to vegetation exceeded that of HA, with CC contributions averaging around 60% in cities, reaching a peak of 85.8%.As depicted in Figure 6c,d, we computed the relative contribution rates of CC and HA to NDVI changes.The total contributions of CC to NDVI were dominant (with relative contribution rates surpassing 50%) in 70.6% of the areas, mainly distributed in the central and northeastern regions of the QTP (Figure 6c).The dominant contribution rates of HA to NDVI change were observed in 29.4% of the QTP, primarily in the northwestern and southeastern regions (Figure 6d). Figure 7 illustrates the magnitude and relative contribution rate of each factor to the NDVI changes at the city scale.Most cities exhibited a CC contribution surpassing 0.001/a, with only three cities, Haixi, Rikaze, and Naqu, registering contributions below this value.Generally, HA contributions did not exceed 0.001/a, with the Ngari Prefecture notably presenting a negative HA contribution.Figure 7b shows that the contribution of CC to vegetation exceeded that of HA, with CC contributions averaging around 60% in cities, reaching a peak of 85.8%.As depicted in Figure 6c,d, we computed the relative contribution rates of CC and HA to NDVI changes.The total contributions of CC to NDVI were dominant (with relative contribution rates surpassing 50%) in 70.6% of the areas, mainly distributed in the central and northeastern regions of the QTP (Figure 6c).The dominant contribution rates of HA to NDVI change were observed in 29.4% of the QTP, primarily in the northwestern and southeastern regions (Figure 6d). Figure 7 illustrates the magnitude and relative contribution rate of each factor to the NDVI changes at the city scale.Most cities exhibited a CC contribution surpassing 0.001/a, with only three cities, Haixi, Rikaze, and Naqu, registering contributions below this value.Generally, HA contributions did not exceed 0.001/a, with the Ngari Prefecture notably presenting a negative HA contribution.Figure 7b shows that the contribution of CC to vegetation exceeded that of HA, with CC contributions averaging around 60% in cities, reaching a peak of 85.8%.

Vegetation and Climate Change
Over the past four decades, the QTP has experienced an increase in both temperature and precipitation [24].The overall temperature of the QTP has exhibited a relatively uniform and significant warming trend (from Figure 8a), particularly in zones H2D3, H1C2, H2D1, H2AB1, H1C1, and H1B1, where average temperatures have climbed by approximately 1.5-2.6 • C from the 1980s to the 2020s.However, areas such as H1D1 and H2D2, situated on the northwestern edge of the QTP, recorded a slight temperature decrease of approximately 0.9-1.1 • C. Precipitation patterns exhibited considerable temporal and spatial variation, with pronounced spatial heterogeneity [5,44].Figure 8b highlights regions within a blue polygon where the coefficient of precipitation variation surpassed 35%.Overall, precipitation in the plateau exhibited an upward trend, especially in the central and western areas, including H1C1, H1C2, H1D1, H2C1, H2D2, and H2D3.In these regions, the average precipitation increased by 24-130 mm compared with 1982.Moreover, the southwestern and southeastern parts of the plateau experienced decreases in precipitation levels by an average of 20-35 mm compared with 1982.

Vegetation and Climate Change
Over the past four decades, the QTP has experienced an increase in both temperature and precipitation [24].The overall temperature of the QTP has exhibited a relatively uniform and significant warming trend (from Figure 8a), particularly in zones H2D3, H1C2, H2D1, H2AB1, H1C1, and H1B1, where average temperatures have climbed by approximately 1.5-2.6 °C from the 1980s to the 2020s.However, areas such as H1D1 and H2D2, situated on the northwestern edge of the QTP, recorded a slight temperature decrease of approximately 0.9-1.1 °C.Precipitation patterns exhibited considerable temporal and spatial variation, with pronounced spatial heterogeneity [5,44].Figure 8b highlights regions within a blue polygon where the coefficient of precipitation variation surpassed 35%.Overall, precipitation in the plateau exhibited an upward trend, especially in the central and western areas, including H1C1, H1C2, H1D1, H2C1, H2D2, and H2D3.In these regions, the average precipitation increased by 24-130 mm compared with 1982.Moreover, the southwestern and southeastern parts of the plateau experienced decreases in precipitation levels by an average of 20-35 mm compared with 1982.Figure 8c illustrates that the NDVI predominantly exhibited a positive correlation with the overall temperature.Specifically, in the southwestern and eastern parts of Qinghai Province, notably northeast of Qinghai Lake, and along the Gansu-Sichuan border the correlation coefficient between the temperature and NDVI exceeded 0.65.These areas, characterized by lower elevations and higher temperatures, facilitated vegetation growth.Conversely, east of 85°E in Tibet, the temperature-NDVI correlation was mildly positive with a coefficient of approximately 0.45.However, certain regions within the Ngari Prefecture showed a negative correlation, with a coefficient of −0.35.In arid regions with sparse vegetation, during this stage, a decrease in temperature and a significant increase in precipitation promoted vegetation growth in the area, resulting in an increase in the NDVI [45].
In most parts of Tibet, the eastern region of Qinghai Province, and the western part adjoining Tibet, the correlation coefficient between the annual precipitation and NDVI was approximately 0.6 (Figure 8d).However, a slight negative correlation existed between the precipitation and NDVI in regions such as the Xinjiang Uygur Autonomous Region adjacent to Tibet, eastern Tibet, and certain parts of Yunnan.In precipitation-rich zones such as the southern Tibetan Plateau, vegetation exhibited a minimal response to Figure 8c illustrates that the NDVI predominantly exhibited a positive correlation with the overall temperature.Specifically, in the southwestern and eastern parts of Qinghai Province, notably northeast of Qinghai Lake, and along the Gansu-Sichuan border the correlation coefficient between the temperature and NDVI exceeded 0.65.These areas, characterized by lower elevations and higher temperatures, facilitated vegetation growth.Conversely, east of 85 • E in Tibet, the temperature-NDVI correlation was mildly positive with a coefficient of approximately 0.45.However, certain regions within the Ngari Prefecture showed a negative correlation, with a coefficient of −0.35.In arid regions with sparse vegetation, during this stage, a decrease in temperature and a significant increase in precipitation promoted vegetation growth in the area, resulting in an increase in the NDVI [45].
In most parts of Tibet, the eastern region of Qinghai Province, and the western part adjoining Tibet, the correlation coefficient between the annual precipitation and NDVI was approximately 0.6 (Figure 8d).However, a slight negative correlation existed between the precipitation and NDVI in regions such as the Xinjiang Uygur Autonomous Region adjacent to Tibet, eastern Tibet, and certain parts of Yunnan.In precipitation-rich zones such as the southern Tibetan Plateau, vegetation exhibited a minimal response to precipitation changes.In the heavily degraded northern margin of the QTP, the NDVI indicated a negative correlation with both the temperature and precipitation, whereas in regions with robust vegetation growth the NDVI typically exhibited a positive correlation with both the temperature and precipitation.
Based on the above analysis, Qinghai and Central Tibet exhibited pronounced warming and humidification trends.Notably, in areas near the Qaidam Basin, there was a significant increase in the NDVI.This suggests that regional climate warming and humidification are conducive to surface greening and vegetation growth.The semi-arid and arid regions of the QTP house most alpine grasslands, which are susceptible to climate warming and are highly sensitive to precipitation changes.Considering the correlation between temperature and precipitation, an increase in temperature enhanced vegetation growth across the entire plateau.Additionally, north of 33 • N, increasing precipitation further promoted vegetation growth.Conversely, south of 33 • N, a decline in precipitation constrained vegetation growth, indicating an intricate relationship between the plateau ecology and precipitation patterns.
Subsequently, six counties in the western Naqu region were selected (as shown in Figure 1b).As a typical study area, the selected areas were situated in northern Tibet, constituting a key pastureland of the QTP.Research has indicated a prevalent greening trend in northern Tibet [3,46], possibly due to its historically sparse vegetation coverage and recent climate "warming and humidification."This can facilitate the growth of alpine vegetation.Figure 9 highlights the noticeable vegetation enhancement in this region, albeit with varying degradations.MODIS data pinpoint degradation is mainly south of 32 • N, particularly in the southeast.In contrast, the AVHRR data indicated degradation predominantly in the central and southern parts, with less accurate southeast monitoring.After 2000, the climate displayed a trend towards "warmer and drier" conditions, further intensifying vegetation degradation.
precipitation changes.In the heavily degraded northern margin of the QTP, the NDVI indicated a negative correlation with both the temperature and precipitation, whereas in regions with robust vegetation growth the NDVI typically exhibited a positive correlation with both the temperature and precipitation.
Based on the above analysis, Qinghai and Central Tibet exhibited pronounced warming and humidification trends.Notably, in areas near the Qaidam Basin, there was a significant increase in the NDVI.This suggests that regional climate warming and humidification are conducive to surface greening and vegetation growth.The semi-arid and arid regions of the QTP house most alpine grasslands, which are susceptible to climate warming and are highly sensitive to precipitation changes.Considering the correlation between temperature and precipitation, an increase in temperature enhanced vegetation growth across the entire plateau.Additionally, north of 33°N, increasing precipitation further promoted vegetation growth.Conversely, south of 33°N, a decline in precipitation constrained vegetation growth, indicating an intricate relationship between the plateau ecology and precipitation patterns.
Subsequently, six counties in the western Naqu region were selected (as shown in Figure 1b).As a typical study area, the selected areas were situated in northern Tibet, constituting a key pastureland of the QTP.Research has indicated a prevalent greening trend in northern Tibet [3,46], possibly due to its historically sparse vegetation coverage and recent climate "warming and humidification."This can facilitate the growth of alpine vegetation.Figure 9 highlights the noticeable vegetation enhancement in this region, albeit with varying degradations.MODIS data pinpoint degradation is mainly south of 32°N, particularly in the southeast.In contrast, the AVHRR data indicated degradation predominantly in the central and southern parts, with less accurate southeast monitoring.After 2000, the climate displayed a trend towards "warmer and drier" conditions, further intensifying vegetation degradation.Temperature manipulation experiments conducted by Wang et al. [24] indicated that the influence of temperature on the NDVI depends on precipitation levels.Specifically, areas with annual precipitation of approximately 500 mm demonstrated pronounced NDVI growth.Increasing evidence has emphasized the crucial influence of precipitation on species richness, alpine grassland ecosystem diversity, primary production, and carbon and water cycle spatial distribution.

Human Drivers
Since the launch of the Qinghai-Tibet Railway and the Grand Western Development Program in 2000, the central and southern QTP has experienced accelerated socioeconomic development (see Figure 10).A surge in sectors such as mining, manufacturing, services, and tourism has led to population migration to the region, resulting in consistent population growth [19].As depicted in Figure 11a, the QTP population has more than doubled, increasing from 4.  Temperature manipulation experiments conducted by Wang et al. [24] indicated that the influence of temperature on the NDVI depends on precipitation levels.Specifically, areas with annual precipitation of approximately 500 mm demonstrated pronounced NDVI growth.Increasing evidence has emphasized the crucial influence of precipitation on species richness, alpine grassland ecosystem diversity, primary production, and carbon and water cycle spatial distribution.

Human Drivers
Since the launch of the Qinghai-Tibet Railway and the Grand Western Development Program in 2000, the central and southern QTP has experienced accelerated socioeconomic development (see Figure 10).A surge in sectors such as mining, manufacturing, services, and tourism has led to population migration to the region, resulting in consistent population growth [19].As depicted in Figure 11a, the QTP population has more than doubled, increasing from 4.3 million in 1970 to 9.4 million in 2021.With societal and economic advancements in China, the primary industry has witnessed substantial growth since 2000.In addition, the secondary and tertiary industries exhibit sharp increases in a linear trend.economic advancements in China, the primary industry has witnessed substantial growth since 2000.In addition, the secondary and tertiary industries exhibit sharp increases in a linear trend.The QTP primarily belongs to two different provinces, Qinghai and Tibet (Figure 11d).Due to the sociocultural differences between both regions, it should be evaluated whether their trajectories have diverged over recent decades.Therefore, this study sourced data on cattle, sheep, horses, and other livestock from the Tibet and Qinghai Province Statistical Yearbooks.These data were then uniformly converted into sheep units to analyze the livestock population trends in both provinces.
Based on the observed trends from 1982 to 2020, we categorized the period into three distinct phases.As depicted in Figure 11c, from 1982 to 2000 Tibet's NDVI exhibited a positive trajectory, growing at 0.001/a (p > 0.05).Concurrently, livestock counts experienced a modest increase of 0.13 million sheep/a (Figure 11b).Both the temperature and precipitation significantly increased during this period (p < 0.05), making the conditions conducive to vegetation growth.From 2000 to 2005, the NDVI declined at a rate of 0.0076/a (p < 0.05).In contrast, livestock numbers increased at a rapid pace of 1.23 million sheep/a.The QTP primarily belongs to two different provinces, Qinghai and Tibet (Figure 11d).Due to the sociocultural differences between both regions, it should be evaluated whether their trajectories have diverged over recent decades.Therefore, this study sourced data on cattle, sheep, horses, and other livestock from the Tibet and Qinghai Province Statistical Yearbooks.These data were then uniformly converted into sheep units to analyze the livestock population trends in both provinces.
Based on the observed trends from 1982 to 2020, we categorized the period into three distinct phases.As depicted in Figure 11c, from 1982 to 2000 Tibet's NDVI exhibited a positive trajectory, growing at 0.001/a (p > 0.05).Concurrently, livestock counts experienced a modest increase of 0.13 million sheep/a (Figure 11b).Both the temperature and precipitation significantly increased during this period (p < 0.05), making the conditions conducive to vegetation growth.From 2000 to 2005, the NDVI declined at a rate of 0.0076/a (p < 0.05).In contrast, livestock numbers increased at a rapid pace of 1.23 million sheep/a.The QTP primarily belongs to two different provinces, Qinghai and Tibet (Figure 11d).Due to the sociocultural differences between both regions, it should be evaluated whether their trajectories have diverged over recent decades.Therefore, this study sourced data on cattle, sheep, horses, and other livestock from the Tibet and Qinghai Province Statistical Yearbooks.These data were then uniformly converted into sheep units to analyze the livestock population trends in both provinces.
Based on the observed trends from 1982 to 2020, we categorized the period into three distinct phases.As depicted in Figure 11c, from 1982 to 2000 Tibet's NDVI exhibited a positive trajectory, growing at 0.001/a (p > 0.05).Concurrently, livestock counts experienced a modest increase of 0.13 million sheep/a (Figure 11b).Both the temperature and precipitation significantly increased during this period (p < 0.05), making the conditions conducive to vegetation growth.From 2000 to 2005, the NDVI declined at a rate of 0.0076/a (p < 0.05).In contrast, livestock numbers increased at a rapid pace of 1.23 million sheep/a.Although temperature increased during this phase, precipitation remained relatively stable, potentially intensifying drought conditions and hampering vegetation greening.During 2005-2020, the livestock population experienced a dramatic decline at the rate of 0.63 million sheep/a (p < 0.05).Simultaneously, both the temperature and precipitation demonstrated no significant alterations, while the NDVI increased at 0.0034/a (p < 0.05).Overall, the livestock count in Tibet in the 2020s was close to that in the 1980s, suggesting a potential continued decline.
The NDVI and livestock population in Qinghai Province were divided into three distinct periods.During 1982-1991, as illustrated in Figure 11f, the NDVI in Qinghai exhibited minimal variation.Livestock figures displayed a negligible increase, peaking in 1991 (Figure 11e).As the temperature increased and precipitation remained relatively constant, both climatic and livestock changes were detrimental to vegetation growth.From 1992 to 2002, the NDVI in Qinghai demonstrated a positive trend, increasing at 0.005/a.Between 1992 and 1999, livestock counts declined substantially at a rate of 1.08 million sheep/a (p < 0.05).During this phase, both the temperature and precipitation increased.From 2003 to 2020, the NDVI in Qinghai showed a marginal growth of 0.001/a (p > 0.05).Between 2000 and 2020, the livestock population decreased (p > 0.05).During this period, both the temperature and precipitation increased, which enhanced vegetation growth.In summary, as inferred from Figure 11e, the livestock numbers in the 2020s were lower than those in the 1980s in Qinghai.
In general, "warming and humidification" contributed to vegetation greening [17].On the one hand, a persistent rise in temperature might amplify extreme heat events, potentially mitigating the positive feedback between temperature and vegetation growth [40], as evidenced in Tibet from 2000 to 2005.On the other hand, the spatially dispersed grazing activities on the plateau suggest that changes in livestock numbers may not immediately influence vegetation feedback, as illustrated by vegetation trends in Qinghai from 1982 to 1990.
Restoring degraded grassland ecosystems can be an intricate and long-term endeavor.In response to highland grassland degradation, the Chinese government has implemented multiple measures since the 1950s, including initiatives such as the "Three-North Shelterbelt," reforestation of neglected farmland, conversion of grazing areas to forests, and targeted grassland restoration projects.By 2012, the data confirmed that the fenced region in northern Qinghai Province reached 3.32 million hectares.In certain degraded alpine grasslands, fencing has been proven to be beneficial for vegetation recovery [47].The "Grassland Ecosystem Compensation Policy" (GECP), introduced in 2010, has not only amplified livestock off-take rates, decreasing livestock counts, but also increased incomes for herders.This policy, in turn, has alleviated human-induced pressure on grasslands, enhancing their quality, although the extent of the improvement varies [46,48].

Limitations
By 2012, fencing enclosed approximately one-third of the northern Tibetan Plateau.Such initiatives have demonstrated the ability to promote vegetation consistency and may even improve plant diversity.Nonetheless, it is crucial to recognize that fencing can impede wildlife mobility, intensify grazing in unfenced areas, and pose significant financial challenges to regional and national authorities [41].The effects of fencing differ across contexts, indicating the need for detailed evaluation in future research.Several studies have highlighted a reduction in vegetation near infrastructure, such as roads, rivers, and lakes [48].Therefore, future research should consider the effects of transportation systems, herder settlement, and urban expansion on grassland dynamics.

Conclusions
Across the QTP, the NDVI exhibited an increase of 0.0013/a without significant variations during winter and spring.The eastern QTP experienced a more pronounced increase in NDVI, whereas only 17.6% of grasslands, primarily located in the northwest and northern periphery, demonstrated a declining trend.Residual analysis indicated that climate change was the primary driver of vegetation shifts for most of the plateau, accounting for approximately 70.6%.In contrast, human activities exhibited a more focal impact in the southeastern section of the plateau, accounting for 29.4% of the changes.Prospectively, the vegetation on the QTP can maintain the greening due to warming and increased humidity.
After 2000, the positive impacts of human activities increased, reflecting the results of various ecological initiatives.Future studies should incorporate the consequences of such interventions on vegetation changes to accurately quantify the influence of human activity on vegetation at the pixel scale and to provide new solutions for the harmonious and friendly coexistence of human beings and nature in the context of climate variability.
Remote Sens. 2023, 15, x FOR PEER REVIEW 4 of 18 river basins, which offer relatively mild conditions.Key administrative and economic centers include Lhasa and Xining, with smaller settlements dispersed throughout the plateau.

Figure 1 .
Figure 1.The vegetation type (a) and DEM (b) on the QTP in China.

Figure 1 .
Figure 1.The vegetation type (a) and DEM (b) on the QTP in China.
Remote Sens. 2023, 15, x FOR PEER REVIEW 8 of 18 absence of NOAA CDR AVHRR NDVI data during 2015-2020, particularly the missing values in 2020.

Figure 3 .
Figure 3. Annual NDVI (a) and NDVI in each season (b) from 1982 to 2020.The spring is March to May, the summer is June to August, the autumn is September to November, and the winter is December to February in the next year.
Remote Sens. 2023, 15, x FOR PEER REVIEW 8 of 18 absence of NOAA CDR AVHRR NDVI data during 2015-2020, particularly the missing values in 2020.

Figure 3 .
Figure 3. Annual NDVI (a) and NDVI in each season (b) from 1982 to 2020.The spring is March to May, the summer is June to August, the autumn is September to November, and the winter is December to February in the next year.

Figure 3 .
Figure 3. Annual NDVI (a) and NDVI in each season (b) from 1982 to 2020.The spring is March to May, the summer is June to August, the autumn is September to November, and the winter is December to February in the next year.
illustrates that winter and spring NDVI values have remained relatively stable in recent years at approximately 0.2.This stability could be attributed to the relatively dry winter climate induced by global warming, which counteracts the vegetation greening trends.However, NDVI trends during summer and autumn closely align with the interannual fluctuations of the QTP because most vegetation Remote Sens. 2023, 15, 4988 9 of 17 in this region consists of seasonal deciduous plants, exhibiting peak growth in July and August within a single growth cycle [41].3.3.Spatial Change 3.3.1.Spatial Variability of Major Change

Figure 4 .
Figure 4. Spatial distribution of interannual NDVI trends during 1982-2020 on the QTP.The upper panel shows NDVI trend for each 1° longitude bin, right panel shows the longitudinal gradient.

Figure 4 .
Figure 4. Spatial distribution of interannual NDVI trends during 1982-2020 on the QTP.The upper panel shows NDVI trend for each 1 • longitude bin, right panel shows the longitudinal gradient.

Figure 5 .
Figure 5.The statistical trends of vegetation in different zones: (a) eco-climatic zones; (b) vegetation land covers (agriculture: all kinds of cropland, forest: broad-leaved and needle-leaved tree cover, grassland: mosaic herbaceous covered grassland and grassland, wetland: shrub or herbaceous covered water, settlement: urban, sparse vegetation: liches and sparse vegetation, bare: bare areas, water: water, other forest: tree cover with flooded); and (c) 500 m interval altitude zones, as well as statistical proportions of areas for different regions (d-f).

Figure 5 .
Figure 5.The statistical trends of vegetation in different zones: (a) eco-climatic zones; (b) vegetation land covers (agriculture: all kinds of cropland, forest: broad-leaved and needle-leaved tree cover, grassland: mosaic herbaceous covered grassland and grassland, wetland: shrub or herbaceous covered water, settlement: urban, sparse vegetation: liches and sparse vegetation, bare: bare areas, water: water, other forest: tree cover with flooded); and (c) 500 m interval altitude zones, as well as statistical proportions of areas for different regions (d-f).
Remote Sens. 2023, 15, x FOR PEER REVIEW 11 of 18 NDVI changes across the QTP, representing 41.2% of the entire region.The cumulative impact of the three climatic factors under CC was positive for 75.7% of the QTP, with most areas exceeding 0.0005/a.The mean CC contribution to the NDVI change of the QTP was 0.0008/a, accounting for 44.1% of the total NDVI change.

Figure 6 .
Figure 6.(a,b) Spatial distribution of interannual  and  trends during 1982-2020.(c,d) Spatial distribution of relative contribution rates in HA and CC to the annual NDVI change on the QTP.

Figure 7 .
Figure 7. (a,b) Contributions and contribution rates of CC and HA to NDVI change at the city scale.

Figure 6 .
Figure 6.(a,b) Spatial distribution of interannual NDV I ha and NDV I cc trends during 1982-2020.(c,d) Spatial distribution of relative contribution rates in HA and CC to the annual NDVI change on the QTP.
Remote Sens. 2023, 15, x FOR PEER REVIEW 11 of 18 NDVI changes across the QTP, representing 41.2% of the entire region.The cumulative impact of the three climatic factors under CC was positive for 75.7% of the QTP, with most areas exceeding 0.0005/a.The mean CC contribution to the NDVI change of the QTP was 0.0008/a, accounting for 44.1% of the total NDVI change.

Figure 6 .
Figure 6.(a,b) Spatial distribution of interannual  and  trends during 1982-2020.(c,d) Spatial distribution of relative contribution rates in HA and CC to the annual NDVI change on the QTP.

Figure 7 .
Figure 7. (a,b) Contributions and contribution rates of CC and HA to NDVI change at the city scale.Figure 7. (a,b) Contributions and contribution rates of CC and HA to NDVI change at the city scale.

Figure 7 .
Figure 7. (a,b) Contributions and contribution rates of CC and HA to NDVI change at the city scale.Figure 7. (a,b) Contributions and contribution rates of CC and HA to NDVI change at the city scale.

Figure 8 .
Figure 8. (a,b) Temporal trends during 1979-2018 of mean annual temperature and annual precipitation.(c,d) Spatial distributions of correlation coefficient between NDVI and mean annual temperature and annual precipitation on the QTP.

Figure 8 .
Figure 8. (a,b) Temporal trends during 1979-2018 of mean annual temperature and annual precipitation.(c,d) Spatial distributions of correlation coefficient between NDVI and mean annual temperature and annual precipitation on the QTP.

Figure 9 .
Figure 9. Spatial distribution of interannual NDVI trends during 2000-2020 on the Naqu, from NOAA CDR AVHRR NDVI (a) and MODIS NDVI (b).The smaller plots at the top right of the figure show the slopes of the annual changes in NDVI for the same region.
3 million in 1970 to 9.4 million in 2021.With societal and

Figure 9 .
Figure 9. Spatial distribution of interannual NDVI trends during 2000-2020 on the Naqu, from NOAA CDR AVHRR NDVI (a) and MODIS NDVI (b).The smaller plots at the top right of the figure show the slopes of the annual changes in NDVI for the same region.

Figure 11 .
Figure 11.(a) The population and the production of three industries over the QTP during 1970-2020.NDVI, livestock population, temperature, and precipitation of the Tibet (b,c) and Qinghai provinces (e,f) over time series.(d) The overview of the administrative boundary between Tibet and Qinghai provinces.

Figure 11 .
Figure 11.(a) The population and the production of three industries over the QTP during 1970-2020.NDVI, livestock population, temperature, and precipitation of the Tibet (b,c) and Qinghai provinces (e,f) over time series.(d) The overview of the administrative boundary between Tibet and Qinghai provinces.

Figure 11 .
Figure 11.(a) The population and the production of three industries over the QTP during 1970-2020.NDVI, livestock population, temperature, and precipitation of the Tibet (b,c) and Qinghai provinces (e,f) over time series.(d) The overview of the administrative boundary between Tibet and Qinghai provinces.

Table 1 .
Overview of the AVHRR bands.

Table 3 .
Calculation of the relative contribution rates of CC and HA to NDVI change.
[41]and east of 98 • E. Sens. 2023,15,x FOR PEER REVIEW 9 of 18 could be attributed to the relatively dry winter climate induced by global warming, which counteracts the vegetation greening trends.However, NDVI trends during summer and autumn closely align with the interannual fluctuations of the QTP because most vegetation in this region consists of seasonal deciduous plants, exhibiting peak growth in July and August within a single growth cycle[41].