1. Introduction
Vegetation, which plays a significant role in controlling desertification and conservation of soil and water [
1], is an important component of terrestrial ecosystems and acts a pivotal part in climate change by affecting carbon storage, hydrological cycle and energy balance [
2,
3]. Therefore, vegetation coverage is an important indicator of ecological environment and global climate change [
4,
5,
6]. For nearly half a century, the continuous impact of climate change and anthropogenic disturbances on global scale has become the main driving force for ecosystem change and has brought a huge impact on the vegetation ecosystem [
7,
8,
9]. Undoubtedly, there is a close internal relationship between vegetation evolution and climate factors, but human activities can also interfere with vegetation evolution [
4]. Understanding vegetation evolution and its response to climate change and human activities is a key task in predicting future ecosystem evolution, providing a basis for ecological management, ecosystem protection and decision-making [
10,
11,
12].
The relationship between climate and vegetation has been widely discussed and studied by many scholars. Xu et al. [
13] investigated the dynamic evolution of vegetation and its relations with climatic factors during 1982–2011 in China. They discovered that the spatiotemporal variations of vegetation dynamic evolution are controlled primarily by temperature and secondly by precipitation. However, the combined effects of temperature and precipitation exhibit strong spatial heterogeneity [
14], and the complexity of the climate–vegetation relationship is also spatially and temporally variable [
15]. For example, in northwestern China, the precipitation might be the key driving factor of vegetation growth [
16], while, in Chinese Loess Plateau, the temperature is a main control factor of the seasonal change of vegetation and precipitation is an important factor for vegetation variation [
17]. In addition, Xie et al. [
18] showed temperature has a positive effect on vegetation in most periods in the semi-humid region but has no significant effect in the arid region, and precipitation has a positive effect on vegetation in summer in the arid region and in autumn in the semi-arid region. Because of the close relationship between climate and vegetation, the climate change inevitably has a certain impact on the vegetation ecosystem. Xu et al. [
19] found that, from 1982 to 2000, global climate change has contributed to an increase in vegetation cover in the Qinghai-Tibet Plateau, and precipitation is the major climatic factor influencing interannual variation of average vegetation cover. Kong et al. [
20] investigated vegetation response to climate change at Northern Hemisphere (NH) scale. Their research results show that factors potentially influencing vegetation growth in different parts of NH were complex and varied; for instance, temperature was recognized as the critical factor behind vegetation greenness in high latitudes especially for spring and autumn temperature in North America and Siberia. Additionally, some studies showed that vegetation changes and the response of vegetation to climate varied with seasons [
18,
21,
22,
23]. The relationship between vegetation and climate is not only the response of vegetation to climate, but also the impact of vegetation on regional climate. Some studies also indicated that vegetation impacts the regional climate by modulating the land–atmosphere exchanges of heat, water and momentum [
24,
25,
26].
Human activities, which is another important factor affecting vegetation, should not be neglected in analyzing the vegetation dynamic evolution under the changing environment [
18]. Zhang et al. [
27] used the Carnegie–Ames–Stanford approach model to assess the status of vegetation in the Three-River Source Region across different periods from 1982 to 2012 and found that human activities had a weak negative impact from 1982 to 2000 and a favorable impact from 2001 to 2012 on vegetation growth or recovery. Hua et al. [
28] found that land use change was the dominant factor driving long-term changes in vegetation greenness in China. Brandt et al. [
29] reported that, in Sub-Saharan Africa, the increases in woody cover were associated with low population growth and were driven by increases in CO
2 in the humid zones and by increase in precipitation in drylands, whereas the decreases in woody cover were associated with high population growth. In addition, Li et al. [
30] analyzed the main characteristics, spatial-temporal distribution and driving forces of vegetation restoration in the Shaanxi–Gansu–Ningxia Region, reporting that human activities are the main driving forces in vegetation restoration. For example, the “Grain for Green Project”, which turns cultivated land into forest land and aims to protect the ecological environment in China, was identified as the main cause leading to gradual vegetation stabilization in the Shaanxi–Gansu–Ningxia Region [
31,
32]. In addition to environmental protection projects, water conservancy project construction may also has an impact on vegetation. For example, Zhang et al. [
33] analyzed the impact of the Three Gorges Water Conservancy Project on environment and reported that cropland, woodland and grassland areas reduced continuously, while river and built-up area increased from 2000 to 2005, and significant changes in land use and vegetation cover have occurred in the Three Gorges Reservoir (TGR) Area.
Climate change and human activities are the two main driving forces of vegetation cover change [
4]. It is a key task to distinguish the effects of human activities and climate change on vegetation evolution for vegetation ecosystem protection and human activities impact assessment. The residual analysis method was firstly applied to discriminate between climate or human-induced dryland degradation by Evans et al. [
34]. Then, this method began to be widely used to distinguish the impact of climate and human activities on vegetation [
35,
36]. For example, Jiang et al. [
37] used a residual analysis trend method to distinguish the effects of climatic change and human activities on vegetation evolution dynamics in Central Asia. Their research works highlighted that sparse vegetation and the degradation of some shrubs in the southern part of the Karakum Desert, the southern Ustyurt Plateau and the wetland delta of the Large Aral Sea were mainly triggered by human activities. Similarly, by using residual analysis, Sun et al. [
38] found that human activities played a major role in vegetation variation of North China. In addition, based on residual analysis of NDVI variation, Wang et al. [
39] pointed out that human activities had either improved or degraded vegetation cover in some parts in southern China. In theory, the residual analysis method is feasible, and establishing an ideal model of vegetation–climate is the key to identify the impact of human activities on vegetation. However, in practice, it is difficult to find the data without human intervention. Therefore, to overcome the shortcomings of residual analysis method, this study only investigated the impact of the TGR impoundment and did not consider other human activities. Then, a vegetation–climate model was established, which is not affected by impoundment, and the impact of water storage was studied by residual analysis.
The Three Gorges Water Conservancy Project at the upper end of the Yangtze River is currently the world’s largest water conservancy project. Many studies in the Three Gorges Reservoir Region focused on the land use/cover changes [
33] and the response of vegetation to natural and anthropogenic driving factors [
40]. These studies did not pay attention to the impact of the Three Gorges Reservoir’s impoundment. At present, there are few studies on the impact of TGR water impoundment on vegetation–climate relationship. Therefore, this study focused on vegetation–climate relationship and the impact of impoundment in different zones. The NDVI was chosen as a valid indictor for vegetation coverage variation, while precipitation and temperature were chosen as climate factors. The aims of this study were as follows: (1) to establish a vegetation–climate model and explore the relationship between them; (2) to propose a new method for analyzing the impact of water impoundment; (3) to quantify the impact of impoundment on climate–vegetation response relationship; and (4) to compare and verify the proposed new method with the existing method. Hopefully, knowledge of the vegetation variation and vegetation–climate relationship will promote the protection of the ecological environment in this region and manage ecosystem under the impacts of climate change and human activities.
4. Discussion
Vegetation in most areas showed an increasing trend (
Figure 2), while, in Sichuan Basin or the lower reaches of the Yangtze River, vegetation decreased insignificantly or significantly. Urban construction land increased in Sichuan Basin, Chongqing, the middle and lower reaches of the Yangtze River and other places (
Figure 3). Combining
Figure 2 and
Figure 3 to speculate, urbanization and intensive population activities may be the reasons for the heterogeneous distribution of vegetation evolution. Furthermore, the relationship between NDVI and precipitation was weaker than that between NDVI and temperature, and the relationships varied by region.
Combining the vegetation–climate regression model and the sensitivity of vegetation to temperature (
Figure 10,
Figure 11 and
Figure 12), it can be seen that NDVI will increase when the temperature increases, but it will decrease when temperature exceeds a certain limit. The response mechanism of vegetation activity to temperature is mainly reflected in the degree of warming influence on the processes of photosynthesis and respiration. Moderate warming can have a positive effect on the enhancement of vegetation activity process. However, excessively high temperature will adversely affect the vegetation activity process. As excessive increase in temperature may accelerate the evaporation of soil and form dryness trends, the vegetation will prevent itself from losing water by reducing leaf area and light saturation point, resulting in a corresponding reduction in vegetation coverage and a limited photosynthesis rate [
71]. On the other hand, the rising temperature increases the rate of autotrophic respiration and transpiration of vegetation, accelerates the consumption of organic matter and leads to a reduction in the net productivity, resulting in the inhibition of vegetation activity.
According to the vegetation–climate regression model and the sensitivity of vegetation to precipitation (
Figure 13), the observation that precipitation has promotion effect and inhibition effect on vegetation growth can be made. Water participates in physiological and biochemical processes such as photosynthesis and transpiration of vegetation. Many nutrients and minerals in the soil can only be absorbed by plants when dissolved in water. Therefore, increasing precipitation will lead to an increase in photosynthetic rate and organic matter production in vegetation and promote vegetation activities such as growth and cover [
72]. However, when the precipitation exceeds the requirement for vegetation, it may also adversely affect vegetation activities such as growth and development indirectly by reducing radiation and increasing relative humidity.
The sensitivity analysis results in
Section 3.5 show that the sensitivity of vegetation to temperature was higher than that of vegetation to precipitation. From a geographical perspective, the study area is located in the subtropical monsoon climate and humid zone, with abundant precipitation, suitable temperature, high soil moisture and small evaporation. Thus, vegetation activity is not restricted by water, and the increase in temperature is conducive to the extension of vegetation growth season and accumulation of dry matter quality [
73].
Residual analysis method and SI and DI indicator method were used to distinguish and quantify the impact of TGR impoundment on the vegetation–climate response relationship. The results in
Figure 14 show that the vegetation–climate relationship may be affected by impoundment, and impact degree decreased as the distance from subarea to TGR increased. Judging from the results, the impact of TGR impoundment may be the reason for vegetation–climate response relationship variation. However, how TGR impoundment affected the relationship is still unknown. In addition, residual analysis method and SI and DI indicator method have certain limitations when they are used to study the impact of human activities. It is more difficult to build a perfectly ideal model without the influence of human activities, but focusing on a specific impact will make it more feasible to build a model that is not affected by this specific impact. Improving methods to make it more widely used is also one of the directions that can be studied in the future.
5. Conclusions
This study used partial correlation analysis, grid point analysis, residual analysis and Mann–Kendall test methods to analyze and quantify the impact of TGR impoundment on the vegetation–climate relationship in the TGRR and its 100-km buffer zone, based on SPOT/VEGETATION NDVI and ERA5 datasets during 1998–2018. Two types of index were proposed and compared with residual analysis method to quantify the impact of TGR impoundment on the vegetation–climate relationship. Finally, the conclusion can be drawn as follows:
In the TGRR and its 100-km buffer zones, NDVI in most areas showed a significant increasing trend. However, in Sichuan Basin or the lower reaches of the Yangtze River, intensive human activities might be the reason for NDVI decreasing significantly. The partial correlation coefficients of NDVI–temperature were higher than those of NDVI–precipitation, and the dynamic response of vegetation cover to temperature and precipitation changes had strong spatial heterogeneity. More importantly, temperature was the main driving factor of vegetation cover change.
The multiple polynomial regression, which simplified relationship and explained the physical mechanism, could effectively describe the response of vegetation–climate before and after water impoundment. The residuals between predicted NDVI and the observed value after impoundment were mostly negative and decreased significantly. In other words, observed NDVI was higher than predicted. The significant decreasing residuals that cannot be explained by independent variables (climatic factors: temperature and precipitation) may be caused by impoundment of TGR. The trends of the residuals were not completely similar in different zones, but this difference was not intuitive.
In this study, SI and DI, which are more intuitive and clearer when displaying the results, were proposed, respectively, to describe and quantify the vegetation’s sensitivity to climate and the difference between before and after water impoundment. In terms of the sensitivity of vegetation to climate, SI and SVI indices were defined, which can effectively quantify the difference of vegetation–climate response between before and after water impoundment. The sensitivity of vegetation to temperature was affected by precipitation and temperature, while the sensitivity of vegetation to precipitation was affected by temperature. The variation of sensitivity after impoundment was slightly more obvious than that before impoundment, and the direction of this sensitivity variation was related to the temperature threshold. In addition, SI results indicate that vegetation was more sensitive to changes in temperature. Hot and humid conditions will be more conducive to vegetation growth. In terms of the difference in the response of NDVI to climate between before and after water impoundment, the DI results show that water impoundment might have an impact on the response relationship of vegetation to climate, and the impact degree decreased with increasing the distance between subarea and TGR.
Comparing residual analysis method and SI and DI indicator method, it can be found that these two methods are essential to quantify the difference of the vegetation–climate relationship between before and after impoundment to reflect the impact on vegetation. They can both quantify and show the difference before and after the impact, but the residual analysis method showed the impact of impoundment through the significant change trend of the residual, while SI and DI indicators directly showed this impact degree. The results of the residual analysis method are less intuitive compared with those of the SI and DI indicators. In addition, it seems that the SI and DI index methods are more effective when comparing the impact of impoundment on vegetation in different areas. Therefore, it may be desirable to combine the residual analysis method with SI and DI indicator method to better and more comprehensively analyze the impact of human activities. A new idea is provided for the study on the impact of human activities in the future.
This study distinguished and quantified the impact of TGR impoundment on the vegetation–climate response relationship. However, its impact is a complex and far-reaching process, which needs to be further studied to find the mechanism and provide a basis for management and decision-making.