Effects of hydrothermal factors and human activities on the vegetation coverage of the Qinghai

Scientific Reports volume 13, Article number: 12488 (2023) Cite this article

22 Accesses

Metrics details

A systematic understanding of the spatio-temporal changes and driving factors in the Qinghai-Tibet Plateau holds significant scientific reference value for the future of ecological sustainable development. This paper utilizes MODIS normalized difference vegetation index (NDVI) and meteorological data to investigate the spatio-temporal changes and driving factors of vegetation coverage in the Qinghai-Tibet Plateau from 2001 to 2020. Methods employed include the dimidiate pixel model, trend analysis, partial correlation analysis, and residual analysis. The results demonstrate a generally fluctuating upward trend in vegetation coverage across the Tibetan Plateau over the past two decades, with spatial expansion occurring from northwest to southeast. Vegetation coverage exhibits a positive correlation with climate factors. Approximately 60.7% of the area showed a positive correlation between vegetation fractional cover (FVC) and precipitation, with 8.66% of the area demonstrating extremely significant (p < 0.05) and significant (p < 0.01) positive correlation. Human activities, on the whole, have contributed to the enhancement of vegetation cover in the Qinghai-Tibet Plateau. The areas where human activities have positively impacted vegetation cover are primarily situated in north-central Qinghai and north of Ngari, while areas experiencing degradation include certain grassland regions in central-eastern Yushu, Nagqu, and Lhasa.

As an important part of terrestrial ecosystems, vegetation coverage and growth condition significantly influences the ecological environment and the capacity for sustainable development1,2. The fractional vegetation cover (FVC), which is defined as the ratio of the vertical projection of ground vegetation to the total area of a specific region, serves as a valuable indicator for assessing vegetation dynamics, ecological changes, and surface vegetation growth3,4,5. Hydrothermal conditions are the primary nonbiological factors that determine vegetation characteristics6. The underlying connection between FVC and climate factors indicates that FVC may indirectly reflect climate change. Hence, the spatiotemporal variation in FVC proves to be an effective means of elucidating vegetation growth dynamics, monitoring drought stress, and evaluating ecosystem quality in a warming world7.

The Qinghai-Tibet Plateau (QTP), situated in Asia, is the largest and highest plateau in the world, often referred to as the "roof of the world" and the "third pole". With a total area exceeding 2.5 million km28,9, it serves as an inland plateau in China8,9. QTP is a sensitive region in terms of climate change and an ecologically fragile region, and its ecological change has an important impact on Asia and even the global climate10. The distinctive climate and geographical characteristics of the QTP have facilitated the development of diverse ecosystem types, such as forests, shrubs, alpine grasslands, alpine meadows, and alpine deserts. The acquisition of vegetation coverage data and the exploration of its spatial and temporal variations are crucial for assessing ecosystem quality on the QTP.

Previous studies have consistently reported an upward trend in vegetation coverage on the QTP. Li et al. observed an increasing trend in QTP's vegetation coverage from 2001 to 201011. Duan et al. analyzed NDVI data from 2000 to 2018 and found that alpine meadows, alpine grasslands, and overall vegetation on the QTP exhibited an increasing trend in NDVI during the growing season12. Similarly, Zhu et al. conducted a comprehensive analysis of time series NDVI data from 2000 to 2018 and confirmed a positive trend in vegetation greenness across the QTP13. Additionally, these studies have demonstrated the influence of climate change and anthropogenic factors on vegetation growth. Zhong et al. reported an overall increase in vegetation density on the QTP from 1999 to 2014, attributed to short-term warming and increased precipitation14. Zhu et al. indicated that the positive effect of climate change on QTP’s vegetation is weakening while the negative impact is increasing, with human activities playing an increasingly negative role15. Han et al. identified a significant correlation between average temperature, total precipitation, and vegetation coverage during the growing season on the QTP, with precipitation being the main controlling factor for vegetation growth16. Huang et al. found that the dominant influence of climate and human activity varied spatiotemporally, with human-dominated regions being smaller than climate-dominated areas17. Sun et al. used the geodetector model to analyze the effects of data on the intensity of human activity and weather on the spatial distribution of vegetation NDVI in the QTP region18.

In conclusion, vegetation changes are driven by climate change and human activities. Climate change manifests primarily through significant alterations in temperature and precipitation, impacting plant photosynthesis, respiration, and overall growth19. Human activities exert both negative and positive influences on vegetation, encompassing urbanization, deforestation, overgrazing, agricultural practices, afforestation, and ecological initiatives20.

Thus, by utilizing climate data and long-term series of MODIS NDVI datasets, among other sources, along with the residual trend analysis method, this study has the following objectives: (1) to analyze the spatiotemporal dynamics of vegetation variation on the QTP over the past two decades; (2) to elucidate the primary factors contributing to vegetation changes since 2001; (3) to quantify the influence of climate change and anthropogenic activities on vegetation variation in the QTP. This study enhances our comprehension of the underlying mechanisms driving vegetation changes in plateau regions, while offering valuable insights for vegetation restoration efforts on the QTP.

The average FVC, precipitation, and temperature of the Tibetan Plateau exhibited a general fluctuating upward trend from 2001 to 2020, as shown in Fig. 1a–c.

Annual average FVC (a), annual precipitation (b) and annual average temperature (c) change of the QTP from 2001 to 2020.

The annual average FVC reached its maximum value of 0.410 in 2020, while the minimum value was 0.381 in 2015. Precipitation and temperature experienced lower values around 2015 and higher values around 2019. A clear relationship can be observed between vegetation growth and hydrothermal factors, with vegetation coverage varying under different hydrothermal conditions. Overall, the average annual vegetation coverage of the QTP has exhibited a gradual increasing trend over the past 20 years; however, the growth rate has been relatively slow at 0.0006 per year.

To provide a more intuitive representation of vegetation coverage, the FVC was categorized into five levels14 (Table 1). The vegetation coverage in the QTP exhibited a distinct spatial pattern of increasing from northwest to southeast in a step-like manner (Fig. 2).

Annual FVC grade distribution of QTP from 2001 to 2020.

Regions with high vegetation coverage were predominantly located in the southeastern areas, characterized by favorable water and heat conditions, encompassing Sichuan Province, Yunnan Province, the southeastern part of Tibet Autonomous Region, the eastern and southern parts of Qinghai Province, and the southern part of Gansu Province. These regions were predominantly covered by coniferous forests and mixed coniferous and broad-leaved forests. Conversely, the regions with low vegetation coverage were situated in the northwestern part of the plateau, characterized by unfavorable water and heat conditions, encompassing the western part of Tibet Autonomous Region, Xinjiang Uygur Autonomous Region, the northwestern part of Qinghai Province, and the western part of Gansu Province. These regions were predominantly covered by grasslands and deserts. This observation indicated a significant correlation between the spatial distribution pattern of vegetation coverage and the hydrothermal conditions.

The statistical analysis revealed that the high and the low vegetation covered area constituted 25.02% and 52.85%,of the total research area, respectively, while the medium and high, the medium, and the medium and low covered area accounted for 7.30%, 6.30%, and 8.53%,respectively.

In order to investigate the spatial evolution characteristics of vegetation coverage on the QTP in more detail, we employed the Theil-Sen Median trend analysis method and the Mann–Kendall significance test method to analyze the trend of variation in vegetation coverage.

The results, as depicted in Fig. 3, revealed that the area with a highly significant increase and significant increase in vegetation coverage constituted 11.94%. This area was predominantly located in the northern region, encompassing the Xinjiang Uygur Autonomous Region (including the Kunlun Mountain system and Aljinshan Mountain system), Gansu Province (Qilian Mountain system), and the northern part of Qinghai Province (Bayankra Mountain system and the vicinity of Qinghai Lake). Conversely, the area with a highly significant decrease and significant decrease accounted for 3.94% of the total research area and was primarily concentrated in the central Tibet Autonomous Region and southwestern Qinghai Province.

Change trend of FVC on QTP from 2001 to 2020.

The spatial distribution of precipitation in the QTP exhibited a gradual increase from northwest to southeast, aligning closely with the distribution pattern of vegetation coverage (Fig. 4).

Average annual precipitation change on the QTP from 2001 to 2020.

Following the partial correlation analysis between vegetation coverage and precipitation, we derived a spatial distribution map illustrating the partial correlation coefficient and the significance of the correlation between FVC and precipitation on the QTP from 2001 to 2020 (Fig. 5).

Partial correlation coefficient and significance between FVC and precipitation on the QTP from 2001 to 2020.

According to the statistical analysis, 60.7% of the research area exhibited a positive correlation between FVC and precipitation. Among these, 8.66% of the area displayed a highly significant and significant positive correlation, encompassing the grassland areas in Gansu Province, the grassland and meadow areas in southern Qinghai Province, as well as the alpine vegetation, grassland, and meadow areas in central Tibet Autonomous Region. The remaining 39.3% of the area exhibited a negative correlation. Of this portion, a small proportion (1.7%) exhibited a highly significant and significant negative correlation, while 33.5% exhibited an insignificant negative correlation. This included the southeast of Qinghai Province, the Chai Damu Basin, Yunnan Province, and the southeast of Tibet Autonomous Region.

The temperature in the QTP was strongly influenced by altitude (Fig. 6), with temperatures decreasing as altitude increases.

Distribution map of annual average temperature and DEM on Qinghai-Tibet Plateau.

By analyzing the spatial distribution of the partial correlation coefficient and significance between vegetation coverage and temperature from 2001 to 2020 (Fig. 7), it revealed that 52.8% of the area exhibited a positive correlation between vegetation coverage and temperature. This positive correlation was primarily observed in most parts of Qinghai Province and the northwest of Tibet Autonomous Region (adjacent to Xinjiang). Among these regions, 2.22% exhibited an extremely significant and significant positive correlation, while 50.58% exhibited an insignificant positive correlation. Additionally, 47.2% of the area displayed a negative correlation, predominantly in the southwest of Tibet Autonomous Region. Of this portion, 2.86% exhibited an extremely significant and significant negative correlation, while 44.34% exhibited an insignificant negative correlation.

Partial correlation coefficient and significance between FVC and temperature on the QTP from 2001 to 2020.

Residual analysis was employed to isolate the influence of human activities on FVC (Fig. 8).

FVC change trend affected by human activities.

The statistical analysis revealed that 39.19% of the area experienced improvement due to human activities, while 34.12% remained relatively unchanged, and 26.69% experienced degradation. Consequently, over the past 20 years, human activities have enhanced the overall vegetation coverage of the QTP, and fostered vegetation growth. With respect to land use type (Fig. 9), the areas that experienced improvement were primarily located in the grasslands of Xinjiang, the central and northern areas of Qinghai, cultivated areas in Xining, the northern and southwestern areas of Ali, the southwestern areas of Shigatse, and Nujiang. Areas with minimal change predominantly consisted of bare lands and forest lands, such as the Qaidam Basin in the Hercynian region, and forest lands in the southern parts of Shannan and Nyingchi. Regions exhibiting a detrimental impact were mainly concentrated in grasslands of Nagqu and Yushu, the eastern part of Guoluo, and parts of Qamdo, Diqing, Liangshan, etc. It indicated a relatively severe degradation trend in human-altered surfaces.

Map of land use types on the QTP.

The vegetation coverage in the QTP exhibited a consistent overall improvement, with a spatial distribution that increased from the northwest to the southeast. This distribution pattern aligns with the prevailing climatic conditions in the QTP. The northwestern region experienced less favorable climatic conditions for vegetation growth compared to the southeastern region. The trend analysis revealed that the vegetation coverage in the QTP has generally improved, but local areas also experienced degradation. These findings aligned with a previous study21,22,23. Additionally, using the Google Earth Engine (GEE) platform, we acquired the total surface water areas of Tibet and Qinghai from 2001 to 2020 (Fig. 10) through processing the v1.3 dataset of JRC Annual Water Classification History.

Interannual variation of total surface water area in Tibet and Qinghai from 2001 to 2020.

It was observed that the areas experiencing a significant or substantial decline corresponded closely to the lake-concentrated regions. In addition, by comparing the vegetation type maps, it showed that the dominant vegetation type in these areas was mainly grasslands. Over time, the expansion of the lake's size led to the inundation of vegetation in that particular area, resulting in a decrease in vegetation coverage. These findings align with previous studies24,25. A stable area existed in the southernmost part, characterized by subtropical and tropical mountain coniferous forest vegetation that maintains its dense coverage throughout the year. Notably, the regions exhibiting an improving trend in QTP's vegetation were primarily concentrated in the western and northern parts, while the eastern and southern regions displayed a degrading trend. This discrepancy may be attributed to higher population density and relatively greater human activity intensity compared to the western and northern regions.

Through partial correlation analysis, a significant positive correlation was observed between vegetation coverage of QTP and both precipitation and air temperature across the majority of the study area. Precipitation and temperature exhibited complementary driving effects, with precipitation having a greater impact on vegetation coverage compared to temperature. The northern regions of QTP predominantly consisted of deserts, where vegetation coverage exhibited a positive correlation with temperature, but a negative correlation with precipitation. The data demonstrated the sensitivity of vegetation coverage to the changes in water and thermal conditions. Our findings revealed that the changes in precipitation, particularly in the northern regions, were the primary driver of vegetation degradation, aligning with previous research26. The occurrence of low temperature and precipitation levels around 2015 provided a plausible explanation for the observed minimum vegetation coverage during that year. Additionally, the relevant research has noted the occurrence of a super El Niño event in 2015 on QTP, which led to a severe drought across almost the entire plateau27. This extreme event further contributed to the lowest vegetation coverage recorded on QTP in 2015. The results demonstrated that favorable hydrothermal conditions played a pivotal role in promoting vegetation growth.

Apart from natural climate factors, human activities significantly influence vegetation growth. Human activities have a more localized impact on vegetation compared to the broad-scale effects of climatic factors, particularly in densely populated regions28. The distinctive geographical conditions of QTP render the majority of the plateau unsuitable for human habitation, particularly in the western regions, resulting in human activities being predominantly concentrated in the eastern part of the plateau. Based on the statistical findings, the areas exhibiting improved vegetation coverage due to human activities accounted for 39.19%, primarily concentrated in Xinjiang, central and northern Qinghai, cultivated areas of Xining, northern and southwestern Ali, as well as the eastern and southwestern regions of Shigatse, among others. The degraded area constituted 26.69%, primarily distributed in the grassland regions of central and eastern Yushu, areas with man-made surface coverage in central Tibet (Nagqu, Lhasa), as well as the eastern regions of Guoluo and Xining, among others. The above conclusions were consistent with those of previous studies28,29. The Northern Steppe region encompasses numerous snow-covered mountains and glaciers, often situated at high altitudes. The vegetation in this area primarily consists of tundra and other alpine plants. Furthermore, the region has a sparse population, limited urban development, and negligible grazing activities, all of which have minimal impact on vegetation coverage. The degraded areas primarily result from intensive animal husbandry, predominantly affecting grassland and meadow vegetation. The expansion of local animal husbandry has accompanied improved living standards. However, the disorderly management of grazing practices and the absence of timely remedial measures have expedited the degradation of grassland vegetation30. Moreover, Xining, the leading contributor to Qinghai Province's GDP, exhibits a significant proportion of man-made surface area. The rapid economic growth, dense population, expanding road construction, and infrastructural development in Xining have resulted in notable human activity footprints that impact vegetation coverage. The favorable impact of human activities on vegetation in the northeastern region of the Tibetan Plateau since 2000 can primarily be attributed to the implementation of four significant ecological projects in China: the Natural Forest Resources Protection Project, the Three Northern Protection Forests Project, the Return of Cultivated Land to Forests Project, and the Grazing Ban Project15,31.

This paper investigated the temporal and spatial variation characteristics of FVC influenced by hydrothermal factors and human activities, providing valuable scientific references for the sustainable ecological development of the Qinghai-Tibetan Plateau (QTP). The findings of this study revealed that vegetation growth on the QTP was influenced by various factors, including temperature, precipitation, and human activities. The main conclusions are as follows:

Over the past 20 years, the overall vegetation coverage of the QTP exhibited an upward trend, with greening predominantly observed in the northern regions and degradation primarily observed in the southwestern regions.

Vegetation response to climatic factors displayed notable spatial variations: vegetation in the northeastern and southwestern regions exhibited negative correlations with temperature and positive correlations with precipitation, while vegetation in the southeast demonstrated negative correlations with precipitation and positive correlations with temperature. Additionally, vegetation in the northwest displayed positive correlations with both climate factors.

Over the past 20 years, human activities had an overall positive effect on the vegetation coverage of the QTP, although significant regional variations existed. In the northeastern part of the plateau, human activities exerted a beneficial impact on vegetation. Nevertheless, the adverse effects of human activities on vegetation in the southeastern part of the plateau can not be ignored, and deserve serious attention to prevent further disturbances.

This study focused on two primary climate factors, namely, annual average temperature and annual precipitation. However, given the complexity of geographical systems and human activities, it is necessary to consider other climatic factors. The factors such as relative humidity, solar radiation, and evapotranspiration play significant roles in climate change and are vital for plant development. Climate conditions represent a highly complex system, and the unique geographical location of the Tibetan Plateau renders altitude a crucial factor influencing plant distribution. Subsequent investigations will incorporate altitude, slope, and diverse meteorological conditions to extensively explore the geographic distribution of FVC under the influence of climate change. Moreover, due to the challenges associated with data collection in the QTP, this study solely examined the impact of human activities based on population density, land use, and livestock areas. However, the inclusion of additional data, such as the distribution of tourist attractions and social points of interest, would significantly enhance the accuracy of analyzing the influence of human activities. Consequently, future research can offer more comprehensive multidimensional analysis outcomes, leading to a clearer understanding of the influencing factors that impact vegetation coverage in the region.

The data was obtained from the National Aeronautics and Space Administration (https://ladsweb.modaps.eosdis.nasa.gov/). The dataset comprise MOD13Q1 grid NDVI data products spanning from 2001 to 2020, featuring a spatial resolution of 250 m and a temporal resolution of 16 days. We performed data extraction, splicing, reprojection, cropping, and other preprocessing using the Modis Reprojection Tool and ArcGIS software. The Maximum Value Composite (MVC) method was applied to mitigate the negative influences of cloud, cloud shadow, aerosol, water vapor, visual angle, and solar altitude angle on vegetation coverage calculation32,33.

The data was obtained from the National Earth System Science Data Center (http://www.geodata.cn). This dataset34 was derived from the global 0.5 climate data provided by CRU (Climatic Research Unit) and the global high-resolution climate data from WorldClim. The dataset was produced in China using the Delta spatial downscaling scheme and validated with 496 independent meteorological observation points. The validation results demonstrated the reliability and accuracy of the climate data. The downscaling process and accuracy evaluation were conducted as follows:

The annual mean, maximum, and minimum variables of temperature (TMPs) and precipitation (PRE) for each month were derived from the CRU time series data. The derived climatological dataset has a spatial resolution of 30′, consistent with the CRU dataset.

The anomalous time series data for each climate variable were calculated using the CRU time series data and the derived climatological datasets.

where \({TMP}_{an\left(yr,m\right)}\) and \({PRE}_{an\left(yr,m\right)}\) are the anomalies for temperatures and precipitation, respectively; \({TMP}_{\left(yr,m\right)}\) and \({PRE}_{\left(yr,m\right)}\) are the absolute temperatures and precipitation values, respectively; \({CRUTMP}_{\left(m\right)}\) and \({CRUPRE}_{\left(m\right)}\) are the 30′ climatology for temperatures and precipitation, respectively; and \(m\) and \(yr\) correspond to month (January–December) and year, respectively.

The 30′ anomaly time series data were spatially interpolated to achieve a higher spatial resolution consistent with the reference dataset from WorldClim.

The high-resolution anomaly time series data were converted to absolute climate time series data using WorldClim's reference dataset.

where \(m\) and \(yr\) are defined as above; \(res\) represents spatial resolution, i.e., 10′, 5′, 2.5′, and 0.5′; \({TMP}_{\left(yr,m,res\right)}\) and \({PRE}_{\left(yr,m,res\right)}\) are the absolute temperatures and precipitation values with a spatial resolution of res, respectively; \({TMP}_{an\left(yr,m,res\right)}\) and \({PRE}_{an\left(yr,m,res\right)}\) represent anomalies with a spatial resolution of res for temperatures and precipitation, respectively; and \({WorldTMP}_{\left(m,res\right)}\) and \({WorldPRE}_{\left(m,res\right)}\) represent climatology datasets from WorldClim at a spatial resolution of res for temperatures and precipitation, respectively.

The original CRU dataset and the downscaled datasets were evaluated using four statistical metrics: Pearson's correlation coefficient (Cor), mean absolute error (MAE), root mean square error (RMSE), and Nash–Sutcliffe efficiency coefficient (NSE). Additionally, the WorldClim data were evaluated at various spatial resolutions by comparing the climatological values of the WorldClim data with observations from the corresponding geographical locations, using the MAE and Cor indices. The sample size corresponds to the number of independent stations.

Downloaded from the National Tibetan Plateau Scientific Data Center35,36,37,38,39 (http://data.tpdc.ac.cn).

Land use type data is the 2020 version of 30 m global land cover data released by the Ministry of Natural Resources of China.

NDVI exhibits a positive correlation with vegetation coverage. Hence, a pixel binary model wass employed in conjunction with NDVI to estimate vegetation coverage40,41. In this paper, we adopted a widely used method to eliminate noise errors. According to the cumulative frequency histogram of NDVI on the whole image, we intercepted the upper and lower threshold values with 5% and 95% confidence levels, and approximated the NDVI values representing bare land and pure vegetation pixels, respectively42. The specific formula is as follows:

where \(FVC\) is the Fractional Vegetation Cover; \(NDVI\) is the mixed image element value; \(NDV{I}_{\text{soil}}\) is the bare soil image element value; \(NDV{I}_{\text{veg}}\) is the pure vegetation image element value.

The Theil Sen Medium trend analysis method is a robust nonparametric estimation algorithm that effectively mitigates the influence of outliers in long time series analysis. It represents an improvement over the least squares linear regression method43,44. The Mann Kendall test is a nonparametric statistical method introduced by Mann in 1945 and subsequently refined by Kendal and Sneyers45,46. One of its advantages is its independence from the assumption of normal distribution for the measured values and linearity of the trend. Moreover, it is unaffected by missing values and outliers47.

Sen's slope is estimated by the formula:

where \(\beta\) is the slope; \(i\) and \(j\) are the time series; \({\text{FV}}{\text{C}}_{i}\) and \(FV{C}_{j}\) are the FVC values in year \(i\) and year \(j\), respectively; \({\text{median}}\) is the median function. When \(\beta\) >0, FVC has an upward trend; when \(\beta\) < 0, FVC has a downward trend.

The results of the Sen trend analysis were tested for significance using the Mann–Kendall method, and the test statistic S was calculated as follows:

We used the test statistic Z for the trend test.

where \(S{\text{gn}}\) is the symbolic function, \(k\) and \(j\) are the time series, and \(n\) is the number of monitoring years.

In this paper, using significance levels of P = 0.05 and 0.01, the trend was considered statistically significant when the absolute value of Z exceeded 1.96 and 2.58, corresponding to the 95% and 99% significance levels, respectively. Table 2 presented the methodology used to assess the significance of the observed trend.

Partial correlation analysis is a valuable method for examining the relationship between two specific variables while controlling for the influence of a third variable. In partial correlation analysis, the influence of a third variable on the relationship between two variables is removed, allowing for an examination of the correlation between the remaining two variables48,49. Using available data on vegetation coverage, precipitation, and temperature, the partial correlation coefficients between vegetation coverage and each of the two climate variables (precipitation and temperature) are computed, and their statistical significance is assessed using t-tests. The specific formula used for calculating the partial correlation coefficients is as follows:

where \({r}_{xy,z}\) represents the partial correlation coefficient between x and y when the z variable is fixed; \({r}_{xy,zw}\) is the second-order partial correlation coefficient between x and y when the z and w variables are fixed; \({r}_{xy}\), \({r}_{xz}\), and \({r}_{yz}\) respectively represent the correlation coefficients between x and y, x and z, y and z; n is the number of years of monitoring; \({x}_{i}\) and \({y}_{i}\) are respectively the values of the x and y variables in year i, and \(\overline{x }\) and \(\overline{y }\) are respectively the annual means of x and y.

Following the computation of each correlation coefficient, a t-test was performed to assess the significance:

where n is the number of monitoring years and q is the number of controllable variables.

Using predetermined significance levels of P = 0.05 and 0.01, the t-values exceeding 2.120 and 2.921, respectively, were considered significant at the 95% and 99% confidence levels.

This method is an effective approach for distinguishing the effects of climate change and human activities on vegetation dynamics50. The impact of human factors is isolated by removing the effects of precipitation and temperature from the long-term time series of vegetation coverage51. Considering that the variations in vegetation coverage are primarily driven by climate and human activities, a multiple regression model is developed in this study to estimate vegetation coverage using precipitation and temperature as predictors. The residual represents the discrepancy between the observed and predicted vegetation coverage, capturing the portion attributed to human activities. Lastly, the linear regression method is employed to analyze the temporal trend of the annual residuals. The formula is expressed as:

where: \(pre\) is annual precipitation data; \(tem\) is annual temperature data; a, b, c are regression coefficients; \(FV{C}_{p}\) is the predicted value; \(FV{C}_{a}\) is the actual Fractional Vegetation Cover calculated from MOD13Q1-NDVI values; \(\delta\) is residual. \(\delta\) > 0 means human activities play a positive role; \(\delta\) < 0 means human activities play a negative role.

The trends of vegetation change under the influence of human activities were classified into five categories: significant degradation52 (< − 0.006), slight degradation (− 0.006 to − 0.0006), basically unchanged (− 0.0006 to 0.0006), slight improvement (0.0006 to 0.006), and significant improvement (> 0.006).

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Wenfu, P. et al. Quantifying influences of natural factors on vegetation NDVI changes based on geographical detector in Sichuan, western China. J. Clean. Prod. 233, 353–367 (2019).

Article Google Scholar

Yang, C. et al. Spatiotemporal changes in vegetation cover and its influencing factors in the Loess plateau of China based on the geographically weighted regression model. Forests 12, 673. https://doi.org/10.3390/f12060673 (2021).

Article Google Scholar

Mu, B. et al. Vegetation cover change and its attribution in China from 2001 to 2018. Remote Sens. 13, 496. https://doi.org/10.3390/rs13030496 (2021).

Article ADS Google Scholar

Miaomiao, L et al. Estimaion of vegetation fraction in the Upper Basin of Miyun Reservoir by remote sensing. Resour. Sci. 26, 153–159 (2004).

Changlong, L. et al. Spatial and temporal characteristics of forest and grass cover in the potential range of desertification in China from 2000 to 2020. Acta Geogr. Sin. 77, 2803–2816 (2022).

Google Scholar

Jiao, K. et al. Precipitation drives the NDVI distribution on the Tibetan Plateau while high warming rates may intensify its ecological droughts. Remote Sens. 13, 1305. https://doi.org/10.3390/rs13071305 (2001).

Article ADS Google Scholar

Fang, H., Li, S., Zhang, Y., Wei, S. & Wang, Y. New insights of global vegetation structural properties through an analysis of canopy clumping index, fractional vegetation cover, and leaf area index. Sci. Remote Sens. 4, 100027 (2021).

Article Google Scholar

Wang, Z. et al. Vegetation expansion on the Tibetan Plateau and its relationship with climate change. Remote Sens. 12, 4150. https://doi.org/10.3390/rs12244150 (2020).

Article ADS Google Scholar

Zhong, L. et al. Assessment of vegetation dynamics and their response to variations in precipitation and temperature in the Tibetan Plateau. Clim. Change 103, 519–535 (2010).

Article ADS Google Scholar

Hua, T. et al. Effectiveness of protected areas edges on vegetation greenness, cover and productivity on the Tibetan Plateau, China. Landsc. Urban Plan. 224, 104421 (2022).

Article Google Scholar

Juncheng, L. et al. Spatiotemporal changes of vegetation cover in response to climate change on the Tibetan Plateau. Acta Geol. Sin. 88, 974–983 (2014).

Article Google Scholar

Duan, H. et al. Spatial and temporal differences in alpine meadow, alpine steppe and all vegetation of the Qinghai-Tibetan Plateau and their responses to climate change. Remote Sens. 13, 669 (2021).

Article ADS Google Scholar

Yuying Z et al. Changes of vegetation greenness and its response to drought-wet variation on the Qingzang Plateau. Chin. J. Plant Ecol. 47(01), 51–64 (2023).

Zhong, L., Ma, Y., Xue, Y. & Piao, S. Climate change trends and impacts on vegetation greening over the Tibetan Plateau. J. Geophys. Res. Atmos. 124, 7540–7552. https://doi.org/10.1029/2019JD030481 (2019).

Article ADS Google Scholar

Zhu, B., Zhang, Z., Tian, J., Kong, R. & Chen, X. Increasing negative impacts of climatic change and anthropogenic activities on vegetation variation on the Qinghai-Tibet plateau during 1982–2019. Remote Sens. 14, 4735. https://doi.org/10.3390/rs14194735 (2022).

Article ADS Google Scholar

Han, H., Yin, Y., Zhao, Y. & Qin, F. Spatiotemporal variations in fractional vegetation cover and their responses to climatic changes on the Qinghai-Tibet plateau. Remote Sens. 15, 2662. https://doi.org/10.3390/rs15102662 (2023).

Article ADS Google Scholar

Huang, K. et al. The influences of climate change and human activities on vegetation dynamics in the Qinghai-Tibet plateau. Remote Sens. 8, 876. https://doi.org/10.3390/rs8100876 (2016).

Article ADS Google Scholar

Sun, L. et al. Impacts of climate change and human activities on NDVI in the Qinghai-Tibet plateau. Remote Sens. 15, 587. https://doi.org/10.3390/rs15030587 (2023).

Article ADS Google Scholar

Wang, Y., Zhang, Z. & Chen, X. Quantifying influences of natural and anthropogenic factors on vegetation changes based on geodetector: A case study in the Poyang Lake Basin, China. Remote Sens. 13, 5081. https://doi.org/10.3390/rs13245081 (2021).

Article ADS Google Scholar

Zhang, F. et al. Changes in forest net primary productivity in the Yangtze river basin and its relationship with climate change and human activities. Remote Sens. 11, 1451. https://doi.org/10.3390/rs11121451 (2019).

Article ADS Google Scholar

Zhang, Y. et al. Spatio-temporal changes of vegetation net primary productivity and its driving factors on the Qinghai-Tibetan Plateau from 2001 to 2017. Remote Sens. 13, 1566 (2021).

Article ADS Google Scholar

Zhang, L., Guo, H., Lei, L. & Yan, D. Monitoring vegetation greenness variations in Qinghai-Tibet Plateau with MODIS vegetation index. In 2011 IEEE International Geoscience and Remote Sensing Symposium, 760–762 (2011)https://doi.org/10.1109/IGARSS.2011.6049241.

Guo, B. et al. Determining the contributions of climate change and human activities to the vegetation NPP dynamics in the Qinghai-Tibet Plateau, China, from 2000 to 2015. Environ. Monit. Assess. 192, 663. https://doi.org/10.1007/s10661-020-08606-6 (2020).

Article PubMed Google Scholar

Guoqing, Z. et al. Progress in remote sensing monitoring of lake area, water level, and volume changes on the Tibetan Plateau. Natl. Remote Sens. Bull. 26, 115–125. https://doi.org/10.11834/jrs.20221171 (2022).

Article Google Scholar

Zhang, G. et al. 100 years of lake evolution over the Qinghai-Tibet Plateau. Earth Syst. Sci. Data. 13, 3951–3966. https://doi.org/10.5194/essd-13-3951-2021 (2021).

Article ADS Google Scholar

Deng, X., Wu, L., He, C. & Shao, H. Study on spatiotemporal variation pattern of vegetation coverage on Qinghai-Tibet Plateau and the analysis of its climate driving factors. Int. J. Environ. Res. Public Health. 19, 8836. https://doi.org/10.3390/ijerph19148836 (2022).

Article PubMed PubMed Central Google Scholar

Lei, Y. et al. Extreme lake level changes on the Tibetan Plateau associated with the 2015/2016 El Niño. Geophys. Res. Lett. 46, 5889–5898. https://doi.org/10.1029/2019GL081946 (2019).

Article ADS Google Scholar

Chen, J. et al. Spatiotemporal variation of vegetation on the Qinghai-Tibet Plateau and the influence of climatic factors and human activities on vegetation trend (2000–2019). Remote Sens. 12, 3150. https://doi.org/10.3390/rs12193150 (2020).

Article ADS Google Scholar

Sun, Y. et al. Grazing intensity and human activity intensity data sets on the Qinghai-Tibetan Plateau during 1990–2015. Geosci. Data J. 9, 140–153. https://doi.org/10.1002/gdj3.127 (2022).

Article ADS MathSciNet Google Scholar

Wei, Y. et al. Dual influence of climate change and anthropogenic activities on the spatiotemporal vegetation dynamics over the Qinghai-Tibetan Plateau from 1981 to 2015. Earth’s Future. 10, e2021EF002566. https://doi.org/10.1029/2021EF002566 (2022).

Article ADS Google Scholar

Guo, Z., Nie, H., Zang, R. & Zhang, Y. Eco-benefit evaluation on natural forest protection in southwest China. J. Inn. Mong. Agric. Univ. 32, 65–72 (2011).

Google Scholar

Holben, B. N. Characteristics of maximum-value composite images from temporal AVHRR data. Int. J. Remote Sens. 7, 1417–1434. https://doi.org/10.1080/01431168608948945 (1986).

Article Google Scholar

Liu, Y. et al. Temporal and spatial variation in vegetation coverage and its response to climatic change in marshes of Sanjiang Plain, China. Atmosphere 13, 2077. https://doi.org/10.3390/atmos13122077 (2022).

Article ADS Google Scholar

Shouzhang, P. et al. 1 km monthly temperature and precipitation dataset for China from 1901 to 2017. Earth Syst. Sci. Data. 11, 1931–1946. https://doi.org/10.5194/essd-11-1931-2019 (2019).

Article Google Scholar

Yili, Z. et al. Integration dataset of Tibet Plateau boundary. (National Tibetan Plateau/Third Pole Environment Data Center, 2019)https://doi.org/10.11888/Geogra.tpdc.270099.CSTR:18406.11.Geogra.tpdc.270099.

Yili, Z. et al. Datasets of the boundary and area of the Tibetan Plateau. Digit. J. Glob. Change Data Repos. https://doi.org/10.3974/geodb.2014.01.12.V1 (2014).

Article Google Scholar

Yili, Z. et al. Boundary data of the Tibetan Plateau (2021 Version). Digit. J. Glob. Change Data Repos. https://doi.org/10.3974/geodb.2021.07.10.V1 (2021).

Article Google Scholar

Yili, Z. et al. Redetermine the region and boundaries of Tibetan Plateau. Geogr. Res. 40, 1543–1553 (2021).

Google Scholar

Yili, Z. et al. Comparison of boundary datasets covering Tibetan Plateau between 2021 and 2014 versions. J. Glob. Change Data Discov. 5, 32–42. https://doi.org/10.3974/geodp.2021.04.04.CSTR:20146.14.2021.04.04 (2021).

Article Google Scholar

Li, X. et al. Changes in vegetation coverage and migration characteristics of center of gravity in the arid desert region of northwest China in 30 recent years. Land. 11, 1688. https://doi.org/10.3390/land11101688 (2022).

Article Google Scholar

Anees, S. A. et al. Estimation of fractional vegetation cover dynamics and its drivers based on multi-sensor data in Dera Ismail Khan, Pakistan. J. King Saud Univ.-Sci. 34, 102217 (2022).

Article Google Scholar

Guangyi, D. et al. Response of vegetation variation to climate change and human activities in semi-arid swamps. Front. Plant Sci. https://doi.org/10.3389/fpls.2022.990592 (2022).

Article Google Scholar

Achite, M. et al. Precipitation and runoff trend analysis in the Wadi Mina Basin (Northern Algeria) using non-parametric tests and the ITA method. Sustainability. 14, 9892. https://doi.org/10.3390/su14169892 (2022).

Article Google Scholar

Xiong, Q. et al. Temporal and spatial variation of land surface temperature in recent 20 years and analysis of the effect of land use in Jiangxi Province, China. Atmosphere 13, 1278. https://doi.org/10.3390/atmos13081278 (2022).

Article ADS Google Scholar

Mann, H. B. Non-parametric tests against trend. Econometrica 13, 245–259 (1945).

Article MathSciNet MATH Google Scholar

Kendall, M. G. Rank Correlation Measures (Charles Griffin, 1975).

Google Scholar

Le, H. T. et al. Characterizing spatiotemporal patterns of mangrove forests in can gio biosphere reserve using sentinel-2 imagery. Appl. Sci. 10, 4058. https://doi.org/10.3390/app10124058 (2020).

Article CAS Google Scholar

Yang, J. et al. Vegetation distribution and vegetation activity changes in the Hengduan Mountains from 1992 to 2020. Acta Geogr. Sin. 77, 2787–2802 (2022).

Google Scholar

Liu, H. et al. Effects of climate and human activities on vegetation cover changes in Danjiangkou water source areas. Trans. Chin. Soc. Agric. Eng. 36, 97–105 (2020).

CAS Google Scholar

Berdimbetov, T. et al. Climatic change and human activities link to vegetation dynamics in the Aral sea basin using NDVI. Earth Syst Environ. 5, 303–318. https://doi.org/10.1007/s41748-021-00224-7 (2021).

Article ADS Google Scholar

Yansong, J. et al. Impacts of climate change and human activities on NDVI change in Eastern Coastal Areas of China. Environ. Sci. https://doi.org/10.13227/j.hjkx.202207039 (2023).

Article Google Scholar

Zhihui, T. et al. Driving mechanism of the spatiotemporal evolution of vegetation in the Yellow River Basin from 2000 to 2020. Environ. Sci. 43, 743–751. https://doi.org/10.13227/j.hjkx.202105213 (2022).

Article Google Scholar

Download references

This work was supported by “The Fundamental Research Funds for Central Public-Interest Scientific Institution, grant number: AR2203”. Acknowledgement for the data support from “National Tibetan Plateau Data Center (http://data.tpdc.ac.cn)”, “Resource and Environmental Science and Data Center of Chinese Academy of Sciences (http://www.resdc.cn)” and “National Earth System Science Data Center (http://www.geodata.cn)”.

Key Laboratory of Surveying and Mapping Science and Geospatial Information Technology of MNR, Chinese Academy of Surveying & Mapping, Beijing, 100036, China

Jianxiao Guo, Liang Zhai & Huiyong Sang

College of Resource Environment and Tourism, Capital Normal University, Beijing, 100048, China

Siyuan Cheng

CPC Central Party School (Chinese Academy of Governance), Beijing, 100089, China

Hongwei Li

You can also search for this author in PubMed Google Scholar

J.G. and H.S. designed the research; J.G. and L.Z. performed the analysis; J.G., L.Z. and H.S. wrote the draft. L.Z., J.G., H.S., S.C., and H.L. contributed to the interpretation of the results and the writing of the paper.

Correspondence to Liang Zhai or Huiyong Sang.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Guo, J., Zhai, L., Sang, H. et al. Effects of hydrothermal factors and human activities on the vegetation coverage of the Qinghai-Tibet Plateau. Sci Rep 13, 12488 (2023). https://doi.org/10.1038/s41598-023-39761-8

Download citation

Received: 25 January 2023

Accepted: 30 July 2023

Published: 01 August 2023

DOI: https://doi.org/10.1038/s41598-023-39761-8

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.