The Influences of Climate Change and Human Activities on Vegetation Dynamics in the Qinghai-Tibet Plateau

Grasslands occupy nearly three quarters of the land surface of the Qinghai-Tibet plateau (QTP) and play a critical role in regulating the ecological functions of the QTP. Ongoing climate change and human interference have greatly affected grasslands on the QTP. Differentiating human-induced and climate-driven vegetation changes is vital for both ecological understanding and the management of husbandry. In this study, we employed statistical analysis of annual records, various sources of remote sensing data, and an ecosystem process model to calculate the relative contribution of climate and human activities to vegetation vigor on the QTP. The temperature, precipitation and the intensity and spatial pattern of livestock grazing differed between the periods prior to and after the year 2000, which led to different vegetation dynamics. Overall, increased temperature and enhanced precipitation favored vegetation growth. However, their combined effects exhibited strong spatial heterogeneity. Specifically, increased temperature restrained vegetation growth in dry steppe regions during a period of slightly increasing precipitation from 1986 to 2000 and in meadow regions during a period of precipitation decline during 2000–2011, thereby making precipitation a dominant factor. An increase in precipitation tended to enhance vegetation growth in wet meadow regions during warm periods, and temperature was the limiting factor in Tibet during dry periods. The dominant role played by climate and human activities differed with location and targeted time period. Areas dominated by human activities are much smaller than those dominated by climate. The effects of grazing on grassland pasture were more obvious under unfavorable climate conditions than under suitable ones.


Introduction
As a dominant terrestrial ecosystem on the Qinghai-Tibet Plateau (QTP), grasslands play a pivotal role in linking the pedospheric, atmospheric, and hydrospheric systems in the region, and even the whole of Asia [1][2][3][4][5].Meanwhile, effects of global climate change have been reported to be more obvious in regions of high altitude [6][7][8][9].With a mean elevation of approximately 4000 m above sea level, temperature on the plateau has been increasing at a much faster rate than in most other regions in the world [10,11].In addition, as a traditional mode of land use on the plateau [11], livestock husbandry has developed rapidly over the past several decades, contributing to significant environmental and grassland degradation [12].In order to recover and protect the environments on the plateau, during the late 20th century, the Chinese government implemented a series of projects, including grazing exclusion and establishment of national conservation areas.To maximize the effectiveness of these projects, regular monitoring of grassland dynamics and identifying its underlying mechanisms are necessary.Knowledge gained therefrom can also be used as a basis for sound policy making.
Due to its capability in measuring chlorophyll and energy absorption, the Normalized Difference Vegetation Index (NDVI) has been widely used as an indicator of numerous vegetative biophysical parameters, such as biomass and productivity [13][14][15].For regional scale studies, the NDVI derived from satellite data can be readily accessed by a wide array of instruments.Recently, the Advanced Very High Resolution Radiometer (AVHRR) developed by the National Oceanic and Atmospheric Administration (NOAA) generated the third-generation Global Inventory Modeling and Mapping Studies (GIMMS3g) NDVI dataset [16].Due to its superior performance relative to prior versions and its capability for providing the longest time-series record from the early 80s, the GIMMS3g dataset has been widely used to identify vegetation activity trends and relationships with climatic factors [17,18].However, studies have suggested that the GIMMS NDVI time-series should be dealt with cautiously when being applied to semiarid and sparsely vegetated areas such as the Sahel and the Tibetan Plateau, since the data have been identified as being biased in these regions [19,20].To overcome this limitation, it is necessary to better integrate optional data sources [21].Since the late 90s, improved sensors are being designed specifically for vegetation monitoring, exemplified by the Terra MODIS (Moderate Resolution Imaging Spectro-radiometer), whose NDVI time-series have been successfully utilized as surrogates for the GIMMS3g dataset in regions where GIMMS3g data is biased [22,23].
To date, numerous studies have analyzed inter-annual variations of the NDVI and their relationship to temperature and precipitation throughout the QTP [24,25].It was found that grassland productivity has been significantly influenced by climate change on the plateau [26,27].Elevated temperature has favored grassland growth [21].Precipitation has been identified as a potential factor driving annual maximum NDVI variability for meadows and grasslands with medium vegetation cover [28].The growing season NDVI of the steppe is driven by both temperature and precipitation [29].In addition, effects of precipitation and temperature on vegetation growth have been found to exhibit strong spatial heterogeneity throughout the arid and semi-arid regions [30,31].In the headstream region of Yangtze River, overgrazing has resulted in severe grassland degradation [32], whereas at the same time some researchers have pointed out that light-intensity grazing improves alpine meadow productivity and promotes adaptation to environment change on the plateau [33].The contradictory findings indicate that livestock husbandry exerts significant effects on grassland ecosystems contributing to high spatial heterogeneity over the entire plateau.
To the authors' best knowledge, previous studies have not distinguished anthropogenic drivers and climate drivers at the county level on a geographic scale.The respective effects of human activity and climate on vegetation growth still remain unclear.In particular, knowledge related to the contributions arising from anthropogenic factors remains extremely scarce due to the absence of data on the intensity of anthropogenic activities.A thorough understanding of the relative contribution of each driving factor is essential for developing sustainable management capable of coping with consequences of climate change.In this study, we used the satellite-based vegetation index NDVI and the climatic dataset along with data on livestock grazing intensity to explore co-variations between vegetation dynamics and environmental variables across the QTP grasslands.Since statistical data on stocking intensity is available only for the years 1986-2011, this study is focused mainly on this 26-year period.The principal objectives of this research have been to: (1) Quantify the spatial and temporal variations of vegetation dynamics in response to climate change and anthropogenic activities on the plateau; (2) Identify the dominant factors regulating vegetation dynamics of the QTP grasslands prior to and post 2000.

Study Area
The Qinghai-Tibet Plateau lies in the central part of the "Third Pole".The plateau is located in the semi-arid/arid alpine climate zone, with an annual precipitation of about 250 mm.For the coldest, the mean temperature is lower than −5 • C and for the warmest month it is lower than 10 • C, respectively.The plateau is separated by the Tanggula Mountains into two zones, along a rainfall gradient from 50 to 200 mm in the western Tibet Autonomous Region (Tibet) to 200-500 mm in the east, including the province of Qinghai [34].There is also a precipitation gradient from 50 to 150 mm in the southeast to 300-450 mm in the northwest.Grassland, occupying nearly 75% of the plateau, is the dominant vegetation type on the QTP (Figure 1).For alpine grasslands, pasture for sheep and yaks is the most common anthropogenic activity [35].Since the early 1990s, governments have launched a series of projects to cope with environmental degradation.Nearly one third of the grasslands on the plateau is protected under the delineation of National Conservation Areas (NCAs).There are three NCAs on the plateau, including the Qiangtang Wildlife Sanctuary, the Kekexili National Reserve and the Sanjiangyuan National Nature Reserve, which were established in 1993, 1995, and 2000, respectively.In addition, each of the counties' boundary in the two provinces were regrouped into seven major subzones with reference to the climatic zone, topography, and the national reserve of the Qinghai-Tibet Plateau (Figure 1).

Study Area
The Qinghai-Tibet Plateau lies in the central part of the "Third Pole".The plateau is located in the semi-arid/arid alpine climate zone, with an annual precipitation of about 250 mm.For the coldest, the mean temperature is lower than −5 °C and for the warmest month it is lower than 10 °C, respectively.The plateau is separated by the Tanggula Mountains into two zones, along a rainfall gradient from 50 to 200 mm in the western Tibet Autonomous Region (Tibet) to 200-500 mm in the east, including the province of Qinghai [34].There is also a precipitation gradient from 50 to 150 mm in the southeast to 300-450 mm in the northwest.Grassland, occupying nearly 75% of the plateau, is the dominant vegetation type on the QTP (Figure 1).For alpine grasslands, pasture for sheep and yaks is the most common anthropogenic activity [35].Since the early 1990s, governments have launched a series of projects to cope with environmental degradation.Nearly one third of the grasslands on the plateau is protected under the delineation of National Conservation Areas (NCAs).There are three NCAs on the plateau, including the Qiangtang Wildlife Sanctuary, the Kekexili National Reserve and the Sanjiangyuan National Nature Reserve, which were established in 1993, 1995, and 2000, respectively.In addition, each of the counties' boundary in the two provinces were regrouped into seven major subzones with reference to the climatic zone, topography, and the national reserve of the Qinghai-Tibet Plateau (Figure 1).

Meteorological Data
Monthly precipitation and air temperature datasets during the period of 1986-2011 were collected from the Chinese Meteorological Data Service [36].In total, data from 75 meteorological stations were utilized (Figure 1).The original meteorological data was interpolated to continuous surface data with an 8 km spatial resolution using the AUSPLINE interpolation program [37].We used the growing season accumulative precipitation (GSP) and growing season mean temperature (GST) as the two main factors influencing climate for each county.

Meteorological Data
Monthly precipitation and air temperature datasets during the period of 1986-2011 were collected from the Chinese Meteorological Data Service [36].In total, data from 75 meteorological stations were utilized (Figure 1).The original meteorological data was interpolated to continuous surface data with an 8 km spatial resolution using the AUSPLINE interpolation program [37].We used the growing season accumulative precipitation (GSP) and growing season mean temperature (GST) as the two main factors influencing climate for each county.

NDVI Time Series
Recognizing that the data may reflect biases in the GIMMS3g NDVI dataset on the Qinghai-Tibet Plateau after the year 2000, we used 8 km and 15-day composite GIMMS3g data [38] for 1986-2000, and both GIMMS3g data and the 1 km and 16-day composite MODIS dataset (MOD13A2) [39] for the period of 2000-2011.Seven MODIS tiles-h24v05, h25v05, h25v06, h26v05, h26v06, h27v05, h27v06-are required for our study region.The MODIS tiles were mosaicked and then reprojected from the Sinusoidal projection to the Albers Conical Equal Area projection.Subset images for the study region were clipped from the GIMMS3g data and projected MODIS mosaics by the QTP boundary.Neither quality control data nor quality flag data were used in this study.The maximum value composite (MVC) method was used to select the higher value of bimonthly NDVI images for obtaining the monthly NDVI.Then, the monthly NDVI from May to September was averaged to yield the growing season NDVI (GSNDVI).GSNDVI images from 1986 to 2011 were developed in stacks for time-series analysis.We then extracted mean GSNDVI values from all pixels within each county's grassland for further analysis.

Livestock Grazing Intensity Data
The socio-economic data was acquired from the yearly statistical book, which provides statistical data at a county resolution.We obtained the livestock type and herd number data for 57 counties in Tibet and 21 counties around Sanjiangyuan in Qinghai Province from 1986 to 2011.The number of different animals from the Yearbooks was converted into sheep equivalent unit number (SUN) to quantify the stocking intensity within each county.Conventionally, 1 sheep and 1 goat equals 1 sheep unit, and other animals such as horses and yaks are equivalent to 5 sheep units [40,41].We used the stocking intensity in 1986 and 2000 as the baseline stocking instensity for each county for 1986-2000 and 2000-2011.The time series for grazing intensity for these two periods was constructed using each year's grazing intensity relative to the year 1986 or 2000.

Trend and Shift Trend Analysis
We utilized linear regression to detect the change rates of vegetation greenness (growing-season NDVI, GSNDVI), GST, GSP and livestock grazing intensity (LI) on the QTP during 1986-2000 and 2000-2011 at both regional and county scale.For example, the change rate for vegetation greenness was calculated as: where when θ slope > 0, the vegetation greenness increased and vice versa.The sequential Mann-Kendall (SQ-MK) test [42] was used to statistically assess whether there was a shift in trends of the climate-driving factors at a regional scale.The SQ-MK test, as a non-parametric "trend shift analysis" technology, has been increasingly applied to detect the change point in time series of vegetation greenness, temperature, vegetation phenology, etc. [43].This test sets up two series, a progressive one uf (t) and a backward one ub(t).If they cross each other and diverge beyond a specific threshold value, then there is a statistically significant trend shift.For the GSP time series during 1986-2011, the application of the SQ-MK test has the following four steps in sequence: 1 For each comparison, the number of cases GSP i > GSP j is counted and indicated by n i at the i-th year, where GSP i (i = 1, 2, . . ., n) and GSP j (j = 1, . . ., i − 1) are the sequential values in the inter-annual GSP series. 2 We defined the SQ-MK test statistic as t i , which can be calculated as follows: 3 The mean and variance of the test statistic are calculated using the following two equations, respectively.
A sequential forward series uf (t) value can be calculated via Equation ( 5): Similarly, a sequential backward ub(t) analysis of the SQ-MK test is calculated starting from the end of the time series data.The interception point of uf (t) and ub(t) can be compared with the threshold value to distinguish whether it is the significant trend shift point.

Explanation of Climatic Factors on Vegetation Dynamics
First, partial correlation analysis was conducted between both climatic and anthropogenic driving factors and GSNDVI at a county resolution to assess the interaction effects of individual environment drivers on grassland dynamics while controlling for other factors.In addition, to separate the relative contribution of each driving factor, the fraction in variation partitioning of the multiple standardized regression was applied.The outputs are the variations explained exclusively by each factor, the variations explained by the overlapping effects of each pair of factors, and the amount explained by the overlapping effects of all the explanatory factors.The amount of variations not explained by the complete model was also provided.
where ED i (i = 1, 2, . . ., n) are environmental drivers, including GSP, GST and LI; GSNDV I and ED i are inter-annual mean values of GSNDVI and environmental drivers, respectively; and σ GSNDV I and σ i are standard deviations of GSNDVI and climatic drivers, respectively.Each time series was normalized by subtracting its mean value and then dividing by its standard deviation.Thus any unit change in each variable has the same statistical meaning.Previous studies have shown that the regression coefficients β i can reflect the relative effect or contribution of each climate driver to GSNDVI and can be compared directly for the independent variables.With the above multi-regression, driving factors with larger standard regression coefficients can be assumed to contribute more to the variation of GSNDVI, and hence identified as the dominant driving factor with respect to the GSNDVI variation.

Model Simulation Analysis
Similar to NDVI, Net Primary Productivity (NPP) is also an effective indicator of vegetation growth condition.The Terrestrial Ecosystem Model (TEM) was used to simulate the climatic Potential Net Primary Productivity (NPP P ) [44], which represents the maximum ecosystem NPP under certain climatic conditions and can be assumed to be driven solely by climatic factors.The Carnegie-Ames-Stanford Approach (CASA) model was used to simulate actual NPP (NPP A ), which is the actually occurring NPP.The NPP A can be assumed to be influenced by both climatic factors and human activities.Thus, the human-induced NPP (NPP H ) affected exclusively by human activities can be simulated as follows: We used the Equation ( 6) to obtain the spatial pattern of NPP and its two components [44].Then, the total amount for each component was calculated for each county, and the contributions of climatic and anthropogenic factors were calculated via Equation (8).The adjusted−R 2 is the contribution of climate to vegetation vigor and 1−R 2 is the unexplained variation possibly caused by human activities, including grazing.If the R 2 is greater than 0.5, then climate drivers are regarded as the dominant driving factor and vice versa.The result of the ecosystem productivity model was comparable to the contribution value from the statistical model's result derived from the Equation ( 6) at the county level.

Variation in Enviromental Variables and Vegetation Vigor Dynamics
Figure 2 shows that the GSP over the whole QTP did not exhibit significant trends or any trend change.The uf(t) and ub(t) curves for GST intersected three times between 1993 and 2000.All of the intersection points were statistically significant (exceeding the confidence level of 95%).Then, the entire plateau was separated into Tibet and Qinghai, for which some trend change points for temperature and precipitation can be identified.Over the period 1986-2011 in Tibet, the uf(t) and ub(t) curves intersected in 1998 for GSP and three times between 1994 and 1998 for GST during 1986-2011.All the intersected points were statistically significant.In contrast, for GSP in Qinghai, no significant trends or trend change points were identified.Points of intersection were just slightly below the confidence level of 95%.For GST in Qinghai, intersection points were also identified between 1995 and 1998.Thus, the year of 2000 was treated as the trend change point (year) for climate.To keep in line with climate driving factors, we separated the livestock grazing intensity time series into segments prior to and post 2000.For remote sensing NDVI data, we also separated GIMMS data into segments prior to and post 2000.MODIS data is available only after 2000.Thus, the separation point for remote sensing data is also in accord with that for the climate data.
We used the Equation ( 6) to obtain the spatial pattern of NPP and its two components [44].Then, the total amount for each component was calculated for each county, and the contributions of climatic and anthropogenic factors were calculated via Equation (8).The adjusted−R 2 is the contribution of climate to vegetation vigor and 1−R 2 is the unexplained variation possibly caused by human activities, including grazing.If the R 2 is greater than 0.5, then climate drivers are regarded as the dominant driving factor and vice versa.The result of the ecosystem productivity model was comparable to the contribution value from the statistical model's result derived from the Equation ( 6) at the county level.

Variation in Enviromental Variables and Vegetation Vigor Dynamics
Figure 2 shows that the GSP over the whole QTP did not exhibit significant trends or any trend change.The uf(t) and ub(t) curves for GST intersected three times between 1993 and 2000.All of the intersection points were statistically significant (exceeding the confidence level of 95%).Then, the entire plateau was separated into Tibet and Qinghai, for which some trend change points for temperature and precipitation can be identified.Over the period 1986-2011in Tibet, the uf(t) and ub(t) curves intersected in 1998 for GSP and three times between 1994 and 1998 for GST during 1986-2011.All the intersected points were statistically significant.In contrast, for GSP in Qinghai, no significant trends or trend change points were identified.Points of intersection were just slightly below the confidence level of 95%.For GST in Qinghai, intersection points were also identified between 1995 and 1998.Thus, the year of 2000 was treated as the trend change point (year) for climate.To keep in line with climate driving factors, we separated the livestock grazing intensity time series into segments prior to and post 2000.For remote sensing NDVI data, we also separated GIMMS data into segments prior to and post 2000.MODIS data is available only after 2000.Thus, the separation point for remote sensing data is also in accord with that for the climate data.For Qinghai, the GST increased significantly years both prior to 2000 (0.05 • C•year −1 ) and afterwards (0.03 • C•year −1 ).In contrast, the GSP exhibited a non-significant decreasing trend (−0.98 mm•year −1 ) over the 26-year period for grasslands across the plateau.For the period of 1986-2000, the GSP exhibited insignificant trends in Qinghai and a significant increasing trend in Tibet (3.9 mm•year −1 ).Then, from 2000 to 2011, GSP exhibited an increasing trend in Qinghai while displaying no significant change in Tibet.With respect to the stocking intensity, livestock numbers in Tibet have increased steadily since 1986 and shown a significant increasing trend from 2000 to 2011.In SJY, the livestock number showed a clear decreasing trend, especially after the Restoration of Grazing Project post 2000.
The vegetation vigor, as indicated by GSNDVI, exhibited increasing trends over the whole plateau for both the periods prior to and post 2000, as shown in Figure 3. Recalling that the entire 26-year period was divided into two segments at the year 2000, the GIMMS3g data showed that the increasing rate was 0.006/decade and 0.012/decade for the periods prior to and post 2000, respectively.After 2000, the MODIS data-based vegetation vigor exhibited an increasing rate of 0.0068/decade.Since GIMMS3g-derived GSNDVI showed an unreliable increase from 2008 to 2011 and MODIS GSNDVI did not, we only detect the county scale pattern of vegetation dynamics from 2000-2011 based on MODIS GSNDVI.Averaged across grasslands for all counties, GST increased significantly (0.05 °C•year −1 , p < 0.05) for the entire plateau during 1986-2011.Across Tibet, the averaged inter-annual GST increased significantly (p < 0.05) for the periods of 1986-2000 (0.06 °C•year −1 ) and 2000-2011 (0.04 °C•year −1 ).For Qinghai, the GST increased significantly years both prior to 2000 (0.05 °C•year −1 ) and afterwards (0.03 °C•year −1 ).In contrast, the GSP exhibited a non-significant decreasing trend (−0.98 mm•year −1 ) over the 26-year period for grasslands across the plateau.For the period of 1986-2000, the GSP exhibited insignificant trends in Qinghai and a significant increasing trend in Tibet (3.9 mm•year −1 ).Then, from 2000 to 2011, GSP exhibited an increasing trend in Qinghai while displaying no significant change in Tibet.With respect to the stocking intensity, livestock numbers in Tibet have increased steadily since 1986 and shown a significant increasing trend from 2000 to 2011.In SJY, the livestock number showed a clear decreasing trend, especially after the Restoration of Grazing Project post 2000.
The vegetation vigor, as indicated by GSNDVI, exhibited increasing trends over the whole plateau for both the periods prior to and post 2000, as shown in Figure 3. Recalling that the entire 26-year period was divided into two segments at the year 2000, the GIMMS3g data showed that the increasing rate was 0.006/decade and 0.012/decade for the periods prior to and post 2000, respectively.After 2000, the MODIS data-based vegetation vigor exhibited an increasing rate of 0.0068/decade.Since GIMMS3g-derived GSNDVI showed an unreliable increase from 2008 to 2011 and MODIS GSNDVI did not, we only detect the county scale pattern of vegetation dynamics from 2000-2011 based on MODIS GSNDVI.For all counties across the plateau, GST showed steady increases during both periods.Thus, GST spatial patterns were not included in Figure 4.For vegetation vigor from 1986 to 2000, only six counties in the southern plateau showed a decreasing trend, while half of the counties in Tibet displayed a decreasing trend from 2000 to 2011.GSP over most counties in Qinghai decreased from 1986 to 2000, while the rest of the plateau experienced an increase in water supply from precipitation in the growing season.During 2000-2011, most counties in Tibet showed a decrease in GSP, the reverse was true for all counties in Qinghai.Prior to 2000, nearly all counties covered by alpine meadow in Tibet exhibited an increasing trend of grazing pressure, while counties in Sanjiangyuan showed a decreasing trend.Intensity of grazing also increased around the border of the Qiangtang Wildlife Sanctuary.The counties across the eastern part of the Qiangtang plateau even exhibited an increasing trend in grazing pressure after 2000.For all counties across the plateau, GST showed steady increases during both periods.Thus, GST spatial patterns were not included in Figure 4.For vegetation vigor from 1986 to 2000, only six counties in the southern plateau showed a decreasing trend, while half of the counties in Tibet displayed a decreasing trend from 2000 to 2011.GSP over most counties in Qinghai decreased from 1986 to 2000, while the rest of the plateau experienced an increase in water supply from precipitation in the growing season.During 2000-2011, most counties in Tibet showed a decrease in GSP, the reverse was true for all counties in Qinghai.Prior to 2000, nearly all counties covered by alpine meadow in Tibet exhibited an increasing trend of grazing pressure, while counties in Sanjiangyuan showed a decreasing trend.Intensity of grazing also increased around the border of the Qiangtang Wildlife Sanctuary.The counties across the eastern part of the Qiangtang plateau even exhibited an increasing trend in grazing pressure after 2000.

Relationship between NDVI and Individual Environmental Factors
Effects of climate and grazing intensity on vegetation vigor were evaluated for the period of 1986-2000 using GIMMS3g and for 2000-2011 using MODIS (See Figure 5).Over the two periods, GSNDVI was significantly correlated with precipitation in more than 70% of the plateau.In less than 20% of the grassland areas, GSNDVI was significantly and negatively correlated with temperature.The total area where each driving factor exerted significant effects on vegetation vigor differed between 1986-2000 and 2000-2011.Prior to 2000, the vegetation dynamic was positively correlated with GST and GSP for nearly 95% and 60% of the plateau, respectively.GST displayed a negative correlation with GSNDVI for less than 50% of the plateau during 1986-2000 and 2000-2011.LI exerted a greater negative effect on GSNDVI from 2000 to 2011 than during 1986-2000.Variations in growing season NDVI over most counties of the plateau were positively correlated with GSP both prior to and after 2000.There were only a few counties situated on the wetter southern plateau whose grassland vegetation vigor showed an insignificant negative correlation with GSP (Figure 6a).During 2000-2011, the extremely dry counties around the Qiangtang Wildlife Sanctuary and the Qaidam Desert displayed insignificant negative correlation (Figure 6d).Grassland productivity is generally low in those counties where the correlation between GSP and vegetation was insignificant.In the arid alpine steppe area of the Qiangtang Plateau, vegetation vigor was

Relationship between NDVI and Individual Environmental Factors
Effects of climate and grazing intensity on vegetation vigor were evaluated for the period of 1986-2000 using GIMMS3g and for 2000-2011 using MODIS (See Figure 5).Over the two periods, GSNDVI was significantly correlated with precipitation in more than 70% of the plateau.In less than 20% of the grassland areas, GSNDVI was significantly and negatively correlated with temperature.The total area where each driving factor exerted significant effects on vegetation vigor differed between 1986-2000 and 2000-2011.Prior to 2000, the vegetation dynamic was positively correlated with GST and GSP for nearly 95% and 60% of the plateau, respectively.GST displayed a negative correlation with GSNDVI for less than 50% of the plateau during 1986-2000 and 2000-2011.LI exerted a greater negative effect on GSNDVI from 2000 to 2011 than during 1986-2000.

Relationship between NDVI and Individual Environmental Factors
Effects of climate and grazing intensity on vegetation vigor were evaluated for the period of 1986-2000 using GIMMS3g and for 2000-2011 using MODIS (See Figure 5).Over the two periods, GSNDVI was significantly correlated with precipitation in more than 70% of the plateau.In less than 20% of the grassland areas, GSNDVI was significantly and negatively correlated with temperature.The total area where each driving factor exerted significant effects on vegetation vigor differed between 1986-2000 and 2000-2011.Prior to 2000, the vegetation dynamic was positively correlated with GST and GSP for nearly 95% and 60% of the plateau, respectively.GST displayed a negative correlation with GSNDVI for less than 50% of the plateau during 1986-2000 and 2000-2011.LI exerted a greater negative effect on GSNDVI from 2000 to 2011 than during 1986-2000.Variations in growing season NDVI over most counties of the plateau were positively correlated with GSP both prior to and after 2000.There were only a few counties situated on the wetter southern plateau whose grassland vegetation vigor showed an insignificant negative correlation with GSP (Figure 6a).During 2000-2011, the extremely dry counties around the Qiangtang Wildlife Sanctuary and the Qaidam Desert displayed insignificant negative correlation (Figure 6d).Grassland productivity is generally low in those counties where the correlation between GSP and vegetation was insignificant.In the arid alpine steppe area of the Qiangtang Plateau, vegetation vigor was Variations in growing season NDVI over most counties of the plateau were positively correlated with GSP both prior to and after 2000.There were only a few counties situated on the wetter southern plateau whose grassland vegetation vigor showed an insignificant negative correlation with GSP (Figure 6a).During 2000-2011, the extremely dry counties around the Qiangtang Wildlife Sanctuary and the Qaidam Desert displayed insignificant negative correlation (Figure 6d).Grassland productivity is generally low in those counties where the correlation between GSP and vegetation was insignificant.In the arid alpine steppe area of the Qiangtang Plateau, vegetation vigor was positively and significantly correlated with GSP from 1986 to 2000.From 2000 to 2011, the semi-arid alpine meadow area in northeast plateau also displayed a positive correlation.
Vegetation dynamics of more than half of the counties over the entire plateau showed a positive correlation with GST (Figure 6b,e).For the arid Qiangtang Wildlife Sanctuary and the wet area of Qinghai, vegetation vigor was negatively correlated with GST before 2000.Then, from 2000 to 2011, the negative correlation was mainly distributed in the alpine meadow counties over the Tibet.Compared with climate factors, from 1986 to 2000, those counties whose vegetation dynamics showed a significant negative correlation with livestock intensity were distributed mainly in the alpine steppe region of the Qiangtang Plateau (Figure 6c).The counties across the arid alpine desert and steppe region show significant negative correlation of GST with vegetation vigor.GST in counties across the Qiangtang plateau showed a significant positive correlation with vegetation vigor just as GSP did.In most regions of the alpine meadow ecosystem, grazing intensity failed to show significant correlations with vegetation dynamics (Figure 6c,f).
Remote Sens. 2016, 8, 876 9 of 18 positively and significantly correlated with GSP from 1986 to 2000.From 2000 to 2011, the semi-arid alpine meadow area in northeast plateau also displayed a positive correlation.Vegetation dynamics of more than half of the counties over the entire plateau showed a positive correlation with GST (Figure 6b,e).For the arid Qiangtang Wildlife Sanctuary and the wet area of Qinghai, vegetation vigor was negatively correlated with GST before 2000.Then, from 2000 to 2011, the negative correlation was mainly distributed in the alpine meadow counties over the Tibet.Compared with climate factors, from 1986 to 2000, those counties whose vegetation dynamics showed a significant negative correlation with livestock intensity were distributed mainly in the alpine steppe region of the Qiangtang Plateau (Figure 6c).The counties across the arid alpine desert and steppe region show significant negative correlation of GST with vegetation vigor.GST in counties across the Qiangtang plateau showed a significant positive correlation with vegetation vigor just as GSP did.In most regions of the alpine meadow ecosystem, grazing intensity failed to show significant correlations with vegetation dynamics (Figure 6 c,f).

The Climatic and Anthropogenic Contribution to Vegetation Vigor
Similar to results of the partial correlation coefficients, the spatial patterns of the dominant driving factors clearly differed between the periods prior to and post 2000.From 1986 to 2000 (Figure 7a), dynamics of grassland vegetation over all counties in QTP were driven primarily by precipitation.Increased precipitation favored vegetation vigor in the alpine steppe more than that in the arid region and the alpine meadow, except for the southern plateau and Sanjiangyuan.Increased precipitation and temperature provided a more favorable environment for alpine meadows across the South plateau and Qinghai.For the semi-arid ecosystem in the southern plateau and Sanjiangyuan, increased temperature was dominant.There are mainly three limited regions where grazing exerted great effects on vegetation dynamics: the alpine steppe around the border of the Qiangtang Wildlife Sanctuary, the southern plateau, and the Sanjiangyuan.With respect to the alpine steppe, where grazing intensity is mainly at a low or medium level, decreased grazing intensity around the border of the Qiangtang Wildlife Sanctuary tended to favor growth of grasslands in these counties (Figures A1a and 7a).For counties across the Sanjiangyuan and southern plateau in which grazing intensity is low or medium, warming and wetter climate improved vegetation growth

The Climatic and Anthropogenic Contribution to Vegetation Vigor
Similar to results of the partial correlation coefficients, the spatial patterns of the dominant driving factors clearly differed between the periods prior to and post 2000.From 1986 to 2000 (Figure 7a), dynamics of grassland vegetation over all counties in QTP were driven primarily by precipitation.Increased precipitation favored vegetation vigor in the alpine steppe more than that in the arid region and the alpine meadow, except for the southern plateau and Sanjiangyuan.Increased precipitation and temperature provided a more favorable environment for alpine meadows across the South plateau and Qinghai.For the semi-arid ecosystem in the southern plateau and Sanjiangyuan, increased temperature was dominant.There are mainly three limited regions where grazing exerted great effects on vegetation dynamics: the alpine steppe around the border of the Qiangtang Wildlife Sanctuary, the southern plateau, and the Sanjiangyuan.With respect to the alpine steppe, where grazing intensity is mainly at a low or medium level, decreased grazing intensity around the border of the Qiangtang Wildlife Sanctuary tended to favor growth of grasslands in these counties (Figures A1a and 7a).For counties across the Sanjiangyuan and southern plateau in which grazing intensity is low or medium, warming and wetter climate improved vegetation growth maintaining the higher livestock carrying capacity of the alpine meadow (Figures A1a and 7a).For counties with high and extremely high grazing intensity, decreasing grazing intensity still exerted a negative impact on the vegetation after severe grazing disturbance (Figures A1a and 7a).
The dominant factors also displayed different patterns during 2000-2011 as compared with 1986-2000 (Figure 7b).As such, 81.6% of the areas in counties with declining vegetation vigor was mainly influenced by temperature or grazing.Meanwhile, precipitation exerted a greater impact on alpine meadows in Qinghai and the Qiangtang Wildlife Sanctuary.In the latter and on the southern edge of southwest plateau where the Himalaya Mountains intrude, precipitation still played a major role in vegetation dynamics.In the south plateau and Qiangtang plateau, temperatures acted as the dominant factor over a much larger area than during 1986-2000.Across the whole Qinghai-Tibet Plateau, the number of counties where grazing had been the dominant factor shrunk from that for the period prior to 2000.All of the counties that experienced enhanced vegetation, driven by the dominant factor of grazing, were found to also experience on average a high or medium decrease in grazing pressure from 2000 to 2011 (Figures A1a and 7b).
Remote Sens. 2016, 8, 876 10 of 18 maintaining the higher livestock carrying capacity of the alpine meadow (Figures A1a and 7a).For counties with high and extremely high grazing intensity, decreasing grazing intensity still exerted a negative impact on the vegetation after severe grazing disturbance (Figures A1a and 7a).
The dominant factors also displayed different patterns during 2000-2011 as compared with 1986-2000 (Figure 7b).As such, 81.6% of the areas in counties with declining vegetation vigor was mainly influenced by temperature or grazing.Meanwhile, precipitation exerted a greater impact on alpine meadows in Qinghai and the Qiangtang Wildlife Sanctuary.In the latter and on the southern edge of southwest plateau where the Himalaya Mountains intrude, precipitation still played a major role in vegetation dynamics.In the south plateau and Qiangtang plateau, temperatures acted as the dominant factor over a much larger area than during 1986-2000.Across the whole Qinghai-Tibet Plateau, the number of counties where grazing had been the dominant factor shrunk from that for the period prior to 2000.All of the counties that experienced enhanced vegetation, driven by the dominant factor of grazing, were found to also experience on average a high or medium decrease in grazing pressure from 2000 to 2011 (Figures A1 and 7b).

Comparison of Statistical Results and Productivity Model Results
The statistical regression model using data from annual records indicated that in most of grasslands in the QTP-82.3%from 1986 to 2000 and 90.6% from 2000 to 2011-vegetation dynamics were driven by climate (Figure 8a,b).The regression model and the process model, as well as the remote sensing model showed that the region dominated by anthropogenic activity comprised less than 15.7% and 39.5% of the vegetated land prior to and post 2000, respectively.Areas in this region as detected by the regression model were smaller and more dispersed than those identified in the ecosystem productivity model.During 2000-2011, these latter areas were found mainly in counties inside Tibet along the Qinghai-Tibet Railway (Figure 8d).This was not revealed by a statistical model since this kind of anthropogenic activity was not recorded in the livestock yearbook.Furthermore, in the ecosystem productivity model, the relationships were first calculated by each pixel, while in the statistical model the analysis was conducted by each county.These scale differences are likely the main cause of these spatial inconsistencies.The two methods both revealed that the total area of climate-dominated grasslands was far larger than that dominated by anthropogenic activities.

Comparison of Statistical Results and Productivity Model Results
The statistical regression model using data from annual records indicated that in most of grasslands in the QTP-82.3%from 1986 to 2000 and 90.6% from 2000 to 2011-vegetation dynamics were driven by climate (Figure 8a,b).The regression model and the process model, as well as the remote sensing model showed that the region dominated by anthropogenic activity comprised less than 15.7% and 39.5% of the vegetated land prior to and post 2000, respectively.Areas in this region as detected by the regression model were smaller and more dispersed than those identified in the ecosystem productivity model.During 2000-2011, these latter areas were found mainly in counties inside Tibet along the Qinghai-Tibet Railway (Figure 8d).This was not revealed by a statistical model since this kind of anthropogenic activity was not recorded in the livestock yearbook.Furthermore, in the ecosystem productivity model, the relationships were first calculated by each pixel, while in the statistical model the analysis was conducted by each county.These scale differences are likely the main cause of these spatial inconsistencies.The two methods both revealed that the total area of climate-dominated grasslands was far larger than that dominated by anthropogenic activities.

Discussion
The significant warming and non-significant precipitation trends across the Qinghai-Tibet Plateau as identified in this study were in accord with findings of previous studies [46,47].However, pronounced shifts in the temporal climate trend occurred around the year 2000.Geographically, Tibet and Qinghai are characterized by different thermal heating systems and separate water vapor sources [48,49].The higher elevation of Tibet than Qinghai can lead to more rapid temperature increases in Tibet.In addition, more dense vegetation in Qinghai relative to that in Tibet can cause stronger evapotranspiration, thereby attenuating surface warming by altering the surface latent heat flux and sensible heat flux [50,51].Furthermore, there is a higher percentage of glacial and lake surface area in Tibet than in Qinghai [52].Under a warming climate, changes in glacier and snow cover altered the proportion of received shortwave radiation (surface albedo) and outgoing longwave radiation (surface emissivity), leading to a greater rate of warming in Tibet than in Qinghai [53][54][55][56].The spatial pattern of precipitation change is much more complicated than that of temperature.The precipitable water in the atmosphere has been shown to have increased significantly since the 1990s due to increased atmospheric moisture from melting glaciers in Tibet [57].Across the Qiangtang plateau and southeast plateau, summer rainfall is controlled mainly by 'up-and-over' moisture transport, which is deep convection over the Indian subcontinent [58].Increased aerosol accumulation over northern India will warm and dry the upper troposphere over the QTP, and strongly modulate ENSO influence [59,60], ultimately leading to a decline in growing season precipitation in the Qiangtang and southeast plateaus.
Under duel pressure from the ongoing changing climate and social-economic development, the Chinese government has taken various measures to tune the local economy and industry in the direction of sustainable development [61].Relatively, the economic structure in Tibet is simpler than

Discussion
The significant warming and non-significant precipitation trends across the Qinghai-Tibet Plateau as identified in this study were in accord with findings of previous studies [46,47].However, pronounced shifts in the temporal climate trend occurred around the year 2000.Geographically, Tibet and Qinghai are characterized by different thermal heating systems and separate water vapor sources [48,49].The higher elevation of Tibet than Qinghai can lead to more rapid temperature increases in Tibet.In addition, more dense vegetation in Qinghai relative to that in Tibet can cause stronger evapotranspiration, thereby attenuating surface warming by altering the surface latent heat flux and sensible heat flux [50,51].Furthermore, there is a higher percentage of glacial and lake surface area in Tibet than in Qinghai [52].Under a warming climate, changes in glacier and snow cover altered the proportion of received shortwave radiation (surface albedo) and outgoing longwave radiation (surface emissivity), leading to a greater rate of warming in Tibet than in Qinghai [53][54][55][56].The spatial pattern of precipitation change is much more complicated than that of temperature.The precipitable water in the atmosphere has been shown to have increased significantly since the 1990s due to increased atmospheric moisture from melting glaciers in Tibet [57].Across the Qiangtang plateau and southeast plateau, summer rainfall is controlled mainly by "up-and-over" moisture transport, which is deep convection over the Indian subcontinent [58].Increased aerosol accumulation over northern India will warm and dry the upper troposphere over the QTP, and strongly modulate ENSO influence [59,60], ultimately leading to a decline in growing season precipitation in the Qiangtang and southeast plateaus.
Under duel pressure from the ongoing changing climate and social-economic development, the Chinese government has taken various measures to tune the local economy and industry in the direction of sustainable development [61].Relatively, the economic structure in Tibet is simpler than in Qinghai, resulting in a greater dependency on husbandry in Tibet.Due to the rapid husbandry development of Tibet, over the last three decades, total livestock numbers increased faster in Tibet than in Qinghai.Sanjiangyuan, the home to 300,000 Tibetan pastoralists, encompasses the headwaters of the Yangtze, Yellow and Mekong Rivers.The fundamental environmental and ecological significance of Sanjiangyuan to China underscores the mounting efforts being paid in every aspect, and Sanjiangyuan is one of the best-protected natural reserves in China [62].All these efforts have obviously slowed the growth rate of livestock numbers in Qinghai.Reduced grazing intensity can lead to a series of positive environmental consequences, such as reducing compaction of surface soil and improving soil water retention and hydraulic conductivity [63].
Both climate change and anthropogenic activity could have a non-negligible influence on vegetation growth [64,65].Consistent with trends of climatic and human activity-related factors, the vigor of grassland vegetation across the plateau displayed different trends for the periods prior to and post 2000.During the last decade of the 20th century, an increase in vegetation vigor occurred across almost all of the plateau.After 2000, the decline in vegetation vigor was detected mainly across the previously highly productive high-cold steppe and meadows of the Qiangtang and south plateaus [66-68], which are the primary fodder-providing regions.This shift created huge challenges for husbandry development in these regions.
For ecosystems other than grasslands, including forests and agricultural lands, effects of human activities are various and complex and thereby hard to quantify.Normally, contributions of human activities to variations in vegetation were calculated from the residuals of the NDVI-soil moisture or NDVI-rainfall regression models by applying residual trend analysis (RESTREND) [69,70].However, the RESTREND can only quantify effects stemming from climatic variables.The Qinghai-Tibet Plateau represents a unique system, where intensity of human interference is quite low and pasture has been a traditional landscape for over a thousand years long [71].Grazing is the primary type of human activity on the plateau, which provides a rare opportunity to simplify the process of quantifying effects of anthropogenic activity.
Similar to previous studies, we found that vegetation was greatly affected by precipitation [72,73].Temperature clearly led to a decline in vegetation vigor in Tibet during the hotter years and in drier regions, which is consistent with the findings of other similar studies [66,74].The spatially-varied dominant effects caused by precipitation and temperature embody the heterogeneous climate sensitivity on the plateau under different climatic conditions [75].Temperature is the main limiting factor under conditions of adequate moisture in areas such as alpine meadows in Qinghai and across the warming meadows in Tibet.Increasing precipitation benefits vegetation growth by providing a sufficient water supply under a continually warming climate and throughout water-limited ecosystems, which has also been confirmed along grassland transects in Inner Mongolia [76].In addition to the direct precipitation from the atmosphere, permafrost soils and melting of snow at high elevations could also serve as an indirectly supply of water.These effects might counterbalance the decline in precipitation at the southwestern border of the Qiangtang Wildlife Sanctuary.
Providing food and habitat for livestock is a vital ecosystem service of grasslands.Traditionally, for the rangeland contracted to individual pastoralists, pastoral residents have tended to maximize their livestock numbers during years of favorable climate and abundance of grassland forage [77].However, the local nomads are not willing to decrease their livestock herds even under situations when fodder is scarce [78].As a result, enhanced livestock intensity, along with increasing temperatures and slightly decreasing precipitation, has clearly resulted in an alarming reduction of vegetation vigor in some regions of Tibet.Impacted vegetation can change the energy balance [79] and moisture conditions [80] on the plateau, which will in turn impact soil and its capacity for water conservation; hence, the "Asian water tower" and "ecological hurdle" functions of the Tibetan Plateau.
In contrast to the results of some studies conducted at a local scale [81], this study identified greater effects on vegetation vigor from climate than from grazing.Findings consistent with ours have been suggested in studies of other similar ecosystems [41,82] and research on mechanisms underlying the scope of climatic driving forces [83].Also, in our research, the decreased grazing pressure can benefit the restoration of grassland in the alpine desert region and regions where grazing pressure is not severe, which is also reflected in Wang's research at a finer spatial resolution [84].In the long run, climate inevitably causes impacts on the vigor of grassland vegetation in the Qinghai-Tibet Plateau, just as it impacted the greenness of the Sahel during the years following the drought in the 1980s [85][86][87].Moderate utilization of pasture by humans affects vegetation vigor in the short run.Corresponding measures may be suggested to be taken by both local governments and herdsman grounded in a thorough understanding of their combined effects and respective contributions.Such knowledge is critical to achieving sustainable management of the pastoral ecosystems of the Qinghai-Tibet Plateau.
The accuracy of this study is constrained by the fact that spatially explicit data on grazing intensity at a regional scale in terms of grazing locations and grazing period is not available.Currently, data from annual records is provided county by county.In the future, more comprehensive and spatially explicit survey data and remote sensing data at finer spatial resolutions are needed to enhance research capabilities.With more refined datasets on enclosures (i.e., fencing) and artificial grassland dataset, we can conduct further research on the lagging effects of grazing on the recovery of grasslands, in particular geographic regions.
Another limitation of the study is that the relationships between vegetation dynamics and potential driving forces have been quantified, but the quantified relationship cannot be directly used to represent underlying causal relationships.Also, there might be lagged effects caused by overgrazing.The underlying mechanisms need to be further explored using a comprehensive model that can incorporate most of the relevant factors [83].

Conclusions
(1) On an overall basis, the grassland ecosystems of the Qinghai Plateau have experienced a climate that is clearly changing along with varying intensities of grazing disturbance.Regions where climate has played a dominant role in affecting vegetation growth are much larger than those similarly affected by grazing, which demonstrates that grassland on the Qinghai-Tibet is a climate-driven system.(2) Prior to and following the year 2000, the climate conditions on the Qinghai-Tibet Plateau displayed different trends.From 1986 to 2000, a wetter and warmer climate improved grassland growth through most of the Plateau.From 2000 to 2011, the drier and hotter climate disfavored grassland growth, especially in the arid regions across Tibet.(3) Under favorable climate conditions, grassland ecosystems such as the alpine meadow in the Qinghai-Tibet Plateau can sustain the traditional grazing intensity.Under the harsher climate conditions, government and herders themselves did not manage and lighten the grazing intensity, causing a decline and even degradation of grassland vegetation.

Appendix A
Figure A1 illustrated the spatial pattern of the grazing pressure over the Qinghai-Tibet Plateau prior to 2000 and post 2000.The grazing pressure is the livestock number (sheep equivalent unit number, SUN) per grassland area for a particular year.The high grazing pressure was mainly distributed on alpine meadows of the Qiangtang Plateau and southern Tibet and Sanjiangyuan.During 1986-2000, the grazing pressure showed a continuous increasing trend mainly over the Southwest Plateau and the South Plateau, where the grazing pressure was at a high level (Figure A1a).Meanwhile, the grazing pressure in the extremely high grazing pressure counties of Qiangtang Pleateau showed a continuous decreasing trend.In Sanjiangyuan, the grazing pressure mostly was decreasing.During 2000-2011, most of the high and extremely high grazing pressure counties showed an upward trend over southwest plateau, south plateau and Qiangtang Plateau (Figure A1b).In Sanjiangyuan, half of the extremely high grazing pressure counties showed an upward trend.Some other studies found that climate is the principal driving force for grassland degradation, whereas human activities are the dominant factor in grassland restoration [84].
Remote Sens. 2016, 8, 876 14 of 18 grazing pressure showed a continuous increasing trend mainly over the Southwest Plateau and the South Plateau, where the grazing pressure was at a high level (Figure A1a).Meanwhile, the grazing pressure in the extremely high grazing pressure counties of Qiangtang Pleateau showed a continuous decreasing trend.In Sanjiangyuan, the grazing pressure mostly was decreasing.During 2000-2011, most of the high and extremely high grazing pressure counties showed an upward trend over southwest plateau, south plateau and Qiangtang Plateau (Figure A1b).In Sanjiangyuan, half of the extremely high grazing pressure counties showed an upward trend.Some other studies found that climate is the principal driving force for grassland degradation, whereas human activities are the dominant factor in grassland restoration [84].

Figure 1 .
Figure 1.The study region and meteorological stations for different grassland types.

Figure 1 .
Figure 1.The study region and meteorological stations for different grassland types.

Figure 2 .
Figure 2. The sequential Mann-Kendall (SQ-MK) test to detect the trend shift points during 1986-2011 of (a) Growing season mean temperature (GST) for Qinghai-Tibet Plateau; (b) Growing season mean precipitation (GSP) for Tibet; (c) GSP for Qinghai; (d) GST for Qinghai-Tibet Plateau'; (e) GST for Tibet; and (f) GST for Qinghai.

Figure 2 .
Figure 2. The sequential Mann-Kendall (SQ-MK) test to detect the trend shift points during 1986-2011 of (a) Growing season mean temperature (GST) for Qinghai-Tibet Plateau; (b) Growing season mean precipitation (GSP) for Tibet; (c) GSP for Qinghai; (d) GST for Qinghai-Tibet Plateau; (e) GST for Tibet; and (f) GST for Qinghai.

Figure 3 .
Figure 3.The trend of vegetation vigor in the Qinghai-Tibet Plateau during the period of 1986-2011, prior to 2000, and post 2000, GIMMS NDVI3g (Global Inventory Modeling and Mapping Studies Normalized Difference Vegetation Index) and MODIS NDVI (Moderate Resolution Imaging Spectroradiometer Normalized Difference Vegetation Index).

Figure 3 .
Figure 3.The trend of vegetation vigor in the Qinghai-Tibet Plateau during the period of 1986-2011, prior to 2000, and post 2000, GIMMS NDVI3g (Global Inventory Modeling and Mapping Studies Normalized Difference Vegetation Index) and MODIS NDVI (Moderate Resolution Imaging Spectro-radiometer Normalized Difference Vegetation Index).

Figure 4 .
Figure 4.The change of Growing season mean precipitation (GSP), Livestock grazing intensity (LI) and vegetation vigor (based on growing season NDVI from Global Inventory Modeling and Mapping Studies Normalized Difference Vegetation Index (GIMMS3g) and Moderate Resolution Imaging Spectro-radiometer (MODIS)) over counties on the Qinghai-Tibet Plateau (Tibet and Sanjiangyuan for LI) during (a) the periods prior to 2000 and (b) the periods post 2000.

Figure 5 .
Figure 5.The percentage of pixels with positive and negative partial correlation coefficients indicating the determination of each environmental factor on the growing season NDVI (GSNDVI).

Figure 4 .
Figure 4.The change of Growing season mean precipitation (GSP), Livestock grazing intensity (LI) and vegetation vigor (based on growing season NDVI from Global Inventory Modeling and Mapping Studies Normalized Difference Vegetation Index (GIMMS3g) and Moderate Resolution Imaging Spectro-radiometer (MODIS)) over counties on the Qinghai-Tibet Plateau (Tibet and Sanjiangyuan for LI) during (a) the periods prior to 2000 and (b) the periods post 2000.

Figure 4 .
Figure 4.The change of Growing season mean precipitation (GSP), Livestock grazing intensity (LI) and vegetation vigor (based on growing season NDVI from Global Inventory Modeling and Mapping Studies Normalized Difference Vegetation Index (GIMMS3g) and Moderate Resolution Imaging Spectro-radiometer (MODIS)) over counties on the Qinghai-Tibet Plateau (Tibet and Sanjiangyuan for LI) during (a) the periods prior to 2000 and (b) the periods post 2000.

Figure 5 .
Figure 5.The percentage of pixels with positive and negative partial correlation coefficients indicating the determination of each environmental factor on the growing season NDVI (GSNDVI).

Figure 5 .
Figure 5.The percentage of pixels with positive and negative partial correlation coefficients indicating the determination of each environmental factor on the growing season NDVI (GSNDVI).

Figure 7 .
Figure 7.The dominant factors-the growing season temperature GST (T), the growing season precipitation GSP (P), and livestock grazing intensity LI (G)-driving variation in vegetation vigor across the Qinghai-Tibet Plateau during (a) the periods prior to 2000 and (b) the periods post 2000.The "+" indicates areas of decreasing trend in vegetation vigor; while the "−" represents areas of declining vegetation vigor.

Figure 7 .
Figure 7.The dominant factors-the growing season temperature GST (T), the growing season precipitation GSP (P), and livestock grazing intensity LI (G)-driving variation in vegetation vigor across the Qinghai-Tibet Plateau during (a) the periods prior to 2000 and (b) the periods post 2000.The "+" indicates areas of decreasing trend in vegetation vigor; while the "−" represents areas of declining vegetation vigor.

Figure 8 .
Figure 8.The dominant factors driving the dynamics of vegetation vigor across Qinghai-Tibet Plateau (a,b): Dominant driving factors during 1986-2000 and 2000-2011 as identified from statistical model; (c,d): Dominant driving factors during 1986-2000 and 2000-2011 with higher variation explained via ecosystem production model.Bar charts in upper-left corner show the percentage of total area impacted by driving factors.Note: Contribution of climate change in (a,b) to vegetation vigor = ∑ GST + GSP contributions.

Figure 8 .
Figure 8.The dominant factors driving the dynamics of vegetation vigor across Qinghai-Tibet Plateau (a,b): Dominant driving factors during 1986-2000 and 2000-2011 as identified from statistical model; (c,d): Dominant driving factors during 1986-2000 and 2000-2011 with higher variation explained via ecosystem production model.Bar charts in upper-left corner show the percentage of total area impacted by driving factors.Note: Contribution of climate change in (a,b) to vegetation vigor = ∑ GST + GSP contributions.

Figure A1 .
Figure A1.The multi-year average grazing pressure and the temporal trend (a) During 1986-2000; (b) During 2000-2011 at a county scale.

Figure A1 .
Figure A1.The multi-year average grazing pressure and the temporal trend (a) During 1986-2000; (b) During 2000-2011 at a county scale.